Dataset Information

Quantum-classical mechanics as an alternative to quantum mechanics in molecular and chemical physics.

ABSTRACT: In quantum mechanics, the theory of quantum transitions is grounded on the convergence of a series of time-dependent perturbation theory. In nuclear and atomic physics, this series converges because the dynamics of quantum transitions (quantum jumps) are absent by definition. In molecular and chemical physics, on the contrary, the dynamics of "quantum" transitions, being determined by the joint motion of a light electron (or electrons) and very heavy nuclei, are present by definition, and the series of time-dependent perturbation theory becomes singular. An exception is the dynamic problem for stationary states in the Born-Oppenheimer adiabatic approximation, when the electronic subsystem turns out to be "off" from the general dynamic process and therefore is not dynamically full-fledged: it only forms an electric potential in which the nuclei oscillate. Removing the aforementioned singularity can be accomplished in two ways. The first method was consisted of introducing an additional postulate in the form of the Franck-Condon principle into molecular quantum mechanics, in which the adiabatic approximation is used. The second method was proposed by the author and consisted of damping the singular dynamics of the joint motion of an electron and nuclei in the intermediate (transient) state of molecular "quantum" transitions by introducing chaos. This chaos arises only during molecular quantum transitions and is called dozy chaos. Formally, the damping is carried out by replacing an infinitely small imaginary addition in the spectral representation of the complete Green's function of the system with its finite quantity. The damping chaos (dozy chaos) leads to the continuity of the energy spectrum in the molecular transient state, which is a sign of classical mechanics. Meanwhile, the initial and final states of the molecule obey quantum mechanics in the adiabatic approximation. Molecular quantum mechanics, which takes into account the chaotic dynamics of the transient state of molecular "quantum" transitions, can be called quantum-classical (dozy-chaos) mechanics. The efficacy of the damping for the aforementioned singularity was previously shown by dozy-chaos mechanics of elementary electron transfers in condensed matter, which is the simplest case of dozy-chaos mechanics, and its applications to a whole number of problems, especially to the optical spectra in polymethine dyes and their aggregates. This paper provides a regular exposition of this dozy-chaos (quantum-classical) mechanics of the elementary electron transfers. The main results of its applications presented in the introduction are also described.

SUBMITTER: Egorov VV

PROVIDER: S-EPMC6909101 | biostudies-literature | 2019 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Quantum-classical mechanics as an alternative to quantum mechanics in molecular and chemical physics.

Egorov Vladimir V VV

Heliyon 20191204 12

In quantum mechanics, the theory of quantum transitions is grounded on the convergence of a series of time-dependent perturbation theory. In nuclear and atomic physics, this series converges because the dynamics of quantum transitions (quantum jumps) are absent by definition. In molecular and chemical physics, on the contrary, the dynamics of "quantum" transitions, being determined by the joint motion of a light electron (or electrons) and very heavy nuclei, are present by definition, and the se ...[more]

PMID: 31872096

Similar Datasets

Project description:Direct molecular dynamics (MD) simulation with ab initio quantum mechanical and molecular mechanical (QM/MM) methods is very powerful for studying the mechanism of chemical reactions in a complex environment but also very time-consuming. The computational cost of QM/MM calculations during MD simulations can be reduced significantly using semiempirical QM/MM methods with lower accuracy. To achieve higher accuracy at the ab initio QM/MM level, a correction on the existing semiempirical QM/MM model is an attractive idea. Recently, we reported a neural network (NN) method as QM/MM-NN to predict the potential energy difference between semiempirical and ab initio QM/MM approaches. The high-level results can be obtained using neural network based on semiempirical QM/MM MD simulations, but the lack of direct MD samplings at the ab initio QM/MM level is still a deficiency that limits the applications of QM/MM-NN. In the present paper, we developed a dynamic scheme of QM/MM-NN for direct MD simulations on the NN-predicted potential energy surface to approximate ab initio QM/MM MD. Since some configurations excluded from the database for NN training were encountered during simulations, which may cause some difficulties on MD samplings, an adaptive procedure inspired by the selection scheme reported by Behler [ Behler Int. J. Quantum Chem. 2015 , 115 , 1032 ; Behler Angew. Chem., Int. Ed. 2017 , 56 , 12828 ] was employed with some adaptions to update NN and carry out MD iteratively. We further applied the adaptive QM/MM-NN MD method to the free energy calculation and transition path optimization on chemical reactions in water. The results at the ab initio QM/MM level can be well reproduced using this method after 2-4 iteration cycles. The saving in computational cost is about 2 orders of magnitude. It demonstrates that the QM/MM-NN with direct MD simulations has great potentials not only for the calculation of thermodynamic properties but also for the characterization of reaction dynamics, which provides a useful tool to study chemical or biochemical systems in solution or enzymes.

Project description:The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.

Project description:Structure and dynamics are both critical to RNA's vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies. We have expanded on atomic-resolution crystal structures of RNAs in functionally relevant states by applying explicit solvent MD simulations to explore their dynamics and conformational changes on the submicrosecond time scale. MD relies on simplified atomistic, pairwise additive interaction potentials (force fields). Because of limited sampling, due to the finite accessible simulation time scale and the approximated force field, high-quality starting structures are required. Despite their imperfection, we find that currently available force fields empower MD to provide meaningful and predictive information on RNA dynamics around a crystallographically defined energy minimum. The performance of force fields can be estimated by precise QM calculations on small model systems. Such calculations agree reasonably well with the Cornell et al. AMBER force field, particularly for stacking and hydrogen-bonding interactions. A final verification of any force field is accomplished by simulations of complex nucleic acid structures. The performance of the Cornell et al. AMBER force field generally corresponds well with and augments experimental data, but one notable exception could be the capping loops of double-helical stems. In addition, the performance of pairwise additive force fields is obviously unsatisfactory for inclusion of divalent cations, because their interactions lead to major polarization and charge-transfer effects neglected by the force field. Neglect of polarization also limits, albeit to a lesser extent, the description accuracy of other contributions, such as interactions with monovalent ions, conformational flexibility of the anionic sugar-phosphate backbone, hydrogen bonding, and solute polarization by solvent. Still, despite limitations, MD simulations are a valid tool for analyzing the structural dynamics of existing experimental structures. Careful analysis of MD simulations can identify problematic aspects of an experimental RNA structure, unveil structural characteristics masked by experimental constraints, reveal functionally significant stochastic fluctuations, evaluate the structural role of base ionization, and predict structurally and potentially functionally important details of the solvent behavior, including the presence of tightly bound water molecules. Moreover, combining classical MD simulations with QM calculations in hybrid QM/MM approaches helps in the assessment of the plausibility of chemical mechanisms of catalytic RNAs (ribozymes). In contrast, the reliable prediction of structure from sequence information is beyond the applicability of MD tools. The ultimate utility of computational studies in understanding RNA function thus requires that the results are neither blindly accepted nor flatly rejected, but rather considered in the context of all available experimental data, with great care given to assessing limitations through the available starting structures, force field approximations, and sampling limitations. The examples given in this Account showcase how the judicious use of basic MD simulations has already served as a powerful tool to help evaluate the role of structural dynamics in biological function of RNA.

Dataset Information

Quantum-classical mechanics as an alternative to quantum mechanics in molecular and chemical physics.

Publications

Quantum-classical mechanics as an alternative to quantum mechanics in molecular and chemical physics.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets