Convergence of empirical distributions in an interpretation of quantum mechanics.
Ontology highlight
ABSTRACT:
SUBMITTER: McKeague IW
PROVIDER: S-EPMC5542025 | biostudies-other | 2016 Aug
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.
| S-EPMC8477969 | biostudies-literature
Project description: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.
| S-EPMC6909101 | biostudies-literature
Project description:We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.
| S-EPMC10107847 | biostudies-literature
Project description:A quantum mechanics/molecular mechanics (QM/MM) implementation that uses the Gaussian electrostatic model (GEM) as the MM force field is presented. GEM relies on the reproduction of electronic density by using auxiliary basis sets to calculate each component of the intermolecular interaction. This hybrid method has been used, along with a conventional QM/MM (point charges) method, to determine the polarization on the QM subsystem by the MM environment in QM/MM calculations on 10 individual H(2)O dimers and a Mg(2+)-H(2)O dimer. We observe that GEM gives the correct polarization response in cases when the MM fragment has a small charge, while the point charges produce significant over-polarization of the QM subsystem and in several cases present an opposite sign for the polarization contribution. In the case when a large charge is located in the MM subsystem, for example, the Mg(2+) ion, the opposite is observed at small distances. However, this is overcome by the use of a damped Hermite charge, which provides the correct polarization response.
| S-EPMC2656107 | biostudies-literature
Project description:We raise fundamental questions about the very meaning of conservation laws in quantum mechanics, and we argue that the standard way of defining conservation laws, while perfectly valid as far as it goes, misses essential features of nature and has to be revisited and extended.
| S-EPMC7817196 | biostudies-literature
Project description:Non-relativistic quantum mechanics is reformulated here based on the idea that relational properties among quantum systems, instead of the independent properties of a quantum system, are the most fundamental elements to construct quantum mechanics. This idea, combining with the emphasis that measurement of a quantum system is a bidirectional interaction process, leads to a new framework to calculate the probability of an outcome when measuring a quantum system. In this framework, the most basic variable is the relational probability amplitude. Probability is calculated as summation of weights from the alternative measurement configurations. The properties of quantum systems, such as superposition and entanglement, are manifested through the rules of counting the alternatives. Wave function and reduced density matrix are derived from the relational probability amplitude matrix. They are found to be secondary mathematical tools that equivalently describe a quantum system without explicitly calling out the reference system. Schrödinger Equation is obtained when there is no entanglement in the relational probability amplitude matrix. Feynman Path Integral is used to calculate the relational probability amplitude, and is further generalized to formulate the reduced density matrix. In essence, quantum mechanics is reformulated as a theory that describes physical systems in terms of relational properties.
| S-EPMC6127378 | biostudies-literature
Project description:Network analysis has become an increasingly prevalent research tool across a vast range of scientific fields. Here, we focus on the particular issue of comparing network statistics, i.e. graph-level measures of network structural features, across multiple networks that differ in size. Although "normalized" versions of some network statistics exist, we demonstrate via simulation why direct comparison is often inappropriate. We consider normalizing network statistics relative to a simple fully parameterized reference distribution and demonstrate via simulation how this is an improvement over direct comparison, but still sometimes problematic. We propose a new adjustment method based on a reference distribution constructed as a mixture model of random graphs which reflect the dependence structure exhibited in the observed networks. We show that using simple Bernoulli models as mixture components in this reference distribution can provide adjusted network statistics that are relatively comparable across different network sizes but still describe interesting features of networks, and that this can be accomplished at relatively low computational expense. Finally, we apply this methodology to a collection of ecological networks derived from the Los Angeles Family and Neighborhood Survey activity location data.
| S-EPMC5052017 | biostudies-literature
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.
| S-EPMC6233882 | biostudies-literature
Project description:For exploration of chemical and biological systems, the combined quantum mechanics and molecular mechanics (QM/MM) and machine learning (ML) models have been developed recently to achieve high accuracy and efficiency for molecular dynamics (MD) simulations. Despite its success on reaction free energy calculations, how to identify new configurations on insufficiently sampled regions during MD and how to update the current ML models with the growing database on the fly are both very important but still challenging. In this article, we apply the QM/MM ML method to solvation free energy calculations and address these two challenges. We employ three approaches to detect new data points and introduce the gradient boosting algorithm to reoptimize efficiently the ML model during ML-based MD sampling. The solvation free energy calculations on several typical organic molecules demonstrate that our developed method provides a systematic, robust, and efficient way to explore new chemistry using ML-based QM/MM MD simulations.
| S-EPMC6448400 | biostudies-literature
Project description:Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long-range effects, to polarization and advanced embedding potentials.
| S-EPMC8341087 | biostudies-literature