1
|
Steffen J. Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package. J Chem Theory Comput 2023; 19:5334-5355. [PMID: 37555628 DOI: 10.1021/acs.jctc.3c00568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.
Collapse
Affiliation(s)
- Julien Steffen
- Chair of Theoretical Chemistry, Friedrich-Alexander University Erlangen-Nürnberg, Egerlandstraße 3, 91058 Erlangen, Bavaria, Germany
| |
Collapse
|
2
|
Alkorta I, Elguero J, Del Bene JE. Perturbing the O-H …O Hydrogen Bond in 1-oxo-3-hydroxy-2-propene. Molecules 2021; 26:3086. [PMID: 34064185 PMCID: PMC8196739 DOI: 10.3390/molecules26113086] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Revised: 05/06/2021] [Accepted: 05/12/2021] [Indexed: 11/16/2022] Open
Abstract
Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to identify and characterize equilibrium structures and transition structures on the 1-oxo-3-hydroxy-2-propene: Lewis acid potential energy surfaces, with the acids LiH, LiF, BeH2, and BeF2. Two equilibrium structures, one with the acid interacting with the C=O group and the other with the interaction occurring at the O-H group, exist on all surfaces. These structures are separated by transition structures that present the barriers to the interconversion of the two equilibrium structures. The structures with the acid interacting at the C=O group have the greater binding energies. Since the barriers to convert the structures with interaction occurring at the O-H group are small, only the isomers with interaction occurring at the C=O group could be experimentally observed, even at low temperatures. Charge-transfer energies were computed for equilibrium structures, and EOM-CCSD spin-spin coupling constants 2hJ(O-O), 1hJ(H-O), and 1J(O-H) were computed for equilibrium and transition structures. These coupling constants exhibit a second-order dependence on the corresponding distances, with very high correlation coefficients.
Collapse
Affiliation(s)
- Ibon Alkorta
- Instituto de Química Médica, CSIC, Juan de la Cierva, 3, E-28006 Madrid, Spain;
| | - José Elguero
- Instituto de Química Médica, CSIC, Juan de la Cierva, 3, E-28006 Madrid, Spain;
| | - Janet E. Del Bene
- Department of Chemistry, Youngstown State University, Youngstown, OH 44555, USA
| |
Collapse
|
3
|
Steffen J. A new class of reaction path based potential energy surfaces enabling accurate black box chemical rate constant calculations. J Chem Phys 2019; 150:154105. [DOI: 10.1063/1.5092589] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Julien Steffen
- Institut für Physikalische Chemie, Christian-Albrechts-Universität, Olshausenstraße 40, D–24098 Kiel, Germany
| |
Collapse
|
4
|
Shannon RJ, Hornung B, Tew DP, Glowacki DR. Anharmonic Molecular Mechanics: Ab Initio Based Morse Parametrizations for the Popular MM3 Force Field. J Phys Chem A 2019; 123:2991-2999. [PMID: 30793911 DOI: 10.1021/acs.jpca.8b12006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.
Collapse
Affiliation(s)
- R J Shannon
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K.,Department of Mechanical Engineering , Stanford University , 452 Escondido Mall , Stanford , California 94305 , United States
| | - B Hornung
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K
| | - D P Tew
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K
| | - D R Glowacki
- School of Chemistry, Cantock's Close , University of Bristol , Bristol BS8 1TS , U.K.,Department of Computer Science , University of Bristol , Bristol BS8 1UB , U.K
| |
Collapse
|
5
|
Naskar P, Talukder S, Chaudhury P, Ghosh S. The effect of stochastic barrier fluctuation on semiclassical transmission probability and Shannon entropy of a symmetric double well potential. INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY 2018; 118:e25667. [DOI: 10.1002/qua.25667] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/19/2023]
Affiliation(s)
- Pulak Naskar
- Department of Chemistry; University of Calcutta, 92 A P C Road; Kolkata 700 009 India
| | - Srijeeta Talukder
- Department of Physical Chemistry; Indian Association for the Cultivation of Science, Jadavpur; Kolkata 700 032 India
| | - Pinaki Chaudhury
- Department of Chemistry; University of Calcutta, 92 A P C Road; Kolkata 700 009 India
| | - Subhasree Ghosh
- Department of Chemistry; Serampore College; Serampore, Hooghly, West Bengal 712201 India
| |
Collapse
|
6
|
Robertson C, Habershon S. Harmonic-phase path-integral approximation of thermal quantum correlation functions. J Chem Phys 2018; 148:102316. [PMID: 29544325 DOI: 10.1063/1.5002189] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an approximation to the thermal symmetric form of the quantum time-correlation function in the standard position path-integral representation. By transforming to a sum-and-difference position representation and then Taylor-expanding the potential energy surface of the system to second order, the resulting expression provides a harmonic weighting function that approximately recovers the contribution of the phase to the time-correlation function. This method is readily implemented in a Monte Carlo sampling scheme and provides exact results for harmonic potentials (for both linear and non-linear operators) and near-quantitative results for anharmonic systems for low temperatures and times that are likely to be relevant to condensed phase experiments. This article focuses on one-dimensional examples to provide insights into convergence and sampling properties, and we also discuss how this approximation method may be extended to many-dimensional systems.
Collapse
Affiliation(s)
- Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing, University Of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University Of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
7
|
Karandashev K, Xu ZH, Meuwly M, Vaníček J, Richardson JO. Kinetic isotope effects and how to describe them. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:061501. [PMID: 29282447 PMCID: PMC5729036 DOI: 10.1063/1.4996339] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/16/2017] [Accepted: 08/23/2017] [Indexed: 06/01/2023]
Abstract
We review several methods for computing kinetic isotope effects in chemical reactions including semiclassical and quantum instanton theory. These methods describe both the quantization of vibrational modes as well as tunneling and are applied to the ⋅H + H2 and ⋅H + CH4 reactions. The absolute rate constants computed with the semiclassical instanton method both using on-the-fly electronic structure calculations and fitted potential-energy surfaces are also compared directly with exact quantum dynamics results. The error inherent in the instanton approximation is found to be relatively small and similar in magnitude to that introduced by using fitted surfaces. The kinetic isotope effect computed by the quantum instanton is even more accurate, and although it is computationally more expensive, the efficiency can be improved by path-integral acceleration techniques. We also test a simple approach for designing potential-energy surfaces for the example of proton transfer in malonaldehyde. The tunneling splittings are computed, and although they are found to deviate from experimental results, the ratio of the splitting to that of an isotopically substituted form is in much better agreement. We discuss the strengths and limitations of the potential-energy surface and based on our findings suggest ways in which it can be improved.
Collapse
Affiliation(s)
- Konstantin Karandashev
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Zhen-Hao Xu
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jeremy O Richardson
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, Eidgenössische Technische Hochschule Zürich (ETHZ), CH-8093 Zürich, Switzerland
| |
Collapse
|
8
|
Steffen J, Hartke B. Cheap but accurate calculation of chemical reaction rate constants from ab initio data, via system-specific, black-box force fields. J Chem Phys 2017; 147:161701. [DOI: 10.1063/1.4979712] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Changing pictures of molecular faces and depths of potential acting on an electron in molecule for intramolecular proton transfer reactions of formic acid and malonaldehyde. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.05.035] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
10
|
Wang W, Zhao Y. The dissociation and recombination rates of CH 4 through the Ni(111) surface: The effect of lattice motion. J Chem Phys 2017; 147:044703. [PMID: 28764359 DOI: 10.1063/1.4995299] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Methane dissociation is a prototypical system for the study of surface reaction dynamics. The dissociation and recombination rates of CH4 through the Ni(111) surface are calculated by using the quantum instanton method with an analytical potential energy surface. The Ni(111) lattice is treated rigidly, classically, and quantum mechanically so as to reveal the effect of lattice motion. The results demonstrate that it is the lateral displacements rather than the upward and downward movements of the surface nickel atoms that affect the rates a lot. Compared with the rigid lattice, the classical relaxation of the lattice can increase the rates by lowering the free energy barriers. For instance, at 300 K, the dissociation and recombination rates with the classical lattice exceed the ones with the rigid lattice by 6 and 10 orders of magnitude, respectively. Compared with the classical lattice, the quantum delocalization rather than the zero-point energy of the Ni atoms further enhances the rates by widening the reaction path. For instance, the dissociation rate with the quantum lattice is about 10 times larger than that with the classical lattice at 300 K. On the rigid lattice, due to the zero-point energy difference between CH4 and CD4, the kinetic isotope effects are larger than 1 for the dissociation process, while they are smaller than 1 for the recombination process. The increasing kinetic isotope effect with decreasing temperature demonstrates that the quantum tunneling effect is remarkable for the dissociation process.
Collapse
Affiliation(s)
- Wenji Wang
- College of Chemistry & Pharmacy, Northwest A&F University, Yangling, 712100 Shaanxi Province, People's Republic of China
| | - Yi Zhao
- State Key Laboratory for Physical Chemistry of Solid Surfaces and Fujian Provincial Key Lab of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
11
|
Wong KY, Xu Y, Xu L. Pitfall in Free-Energy Simulations on Simplest Systems. ChemistrySelect 2017. [DOI: 10.1002/slct.201601160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Kin-Yiu Wong
- Department of Physics; High Performance Cluster Computing Centre; Institute of Computational and Theoretical Studies; Hong Kong Baptist University; 224 Waterloo Road Kowloon Tong Hong Kong
- Institute of Research and Continuing Education; Hong Kong Baptist University (Shenzhen); Shenzhen China
| | - Yuqing Xu
- Department of Physics; High Performance Cluster Computing Centre; Institute of Computational and Theoretical Studies; Hong Kong Baptist University; 224 Waterloo Road Kowloon Tong Hong Kong
- Institute of Research and Continuing Education; Hong Kong Baptist University (Shenzhen); Shenzhen China
| | - Liang Xu
- Department of Physics; High Performance Cluster Computing Centre; Institute of Computational and Theoretical Studies; Hong Kong Baptist University; 224 Waterloo Road Kowloon Tong Hong Kong
| |
Collapse
|
12
|
Mazzuca JW, Schultz CP. Quantum Mechanical Enhancement of Rate Constants and Kinetic Isotope Effects for Water-Mediated Proton Transfer in a Model Biological System. J Phys Chem A 2017; 121:819-826. [DOI: 10.1021/acs.jpca.6b10337] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- James W. Mazzuca
- Chemistry Department, Alma College, Alma, Michigan 48801, United States
| | - Chase P. Schultz
- Chemistry Department, Alma College, Alma, Michigan 48801, United States
| |
Collapse
|
13
|
Ceriotti M, Fang W, Kusalik PG, McKenzie RH, Michaelides A, Morales MA, Markland TE. Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges. Chem Rev 2016; 116:7529-50. [DOI: 10.1021/acs.chemrev.5b00674] [Citation(s) in RCA: 339] [Impact Index Per Article: 42.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Michele Ceriotti
- Laboratory
of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Wei Fang
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Peter G. Kusalik
- Department
of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Ross H. McKenzie
- School
of Mathematics and Physics, University of Queensland, Brisbane, 4072 Queensland Australia
| | - Angelos Michaelides
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Thomas E. Markland
- Department
of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305, United States
| |
Collapse
|
14
|
Mori Y, Masuda Y. Effect of solvent on proton location and dynamic behavior in short intramolecular hydrogen bonds studied by molecular dynamics simulations and NMR experiments. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2015.07.002] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
15
|
Cecchini M. Quantum Corrections to the Free Energy Difference between Peptides and Proteins Conformers. J Chem Theory Comput 2015; 11:4011-22. [DOI: 10.1021/acs.jctc.5b00260] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marco Cecchini
- Laboratoire d’Ingénierie
des Fonctions Moléculaires Institut de Science et d’Ingénierie
Supramoléculaires, Université de Strasbourg, 8 allée
Gaspard Monge, F-67083 Strasbourg Cedex, France
| |
Collapse
|
16
|
Kubař T, Welke K, Groenhof G. New QM/MM implementation of the DFTB3 method in the gromacs package. J Comput Chem 2015; 36:1978-89. [PMID: 26238364 DOI: 10.1002/jcc.24029] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Accepted: 07/09/2015] [Indexed: 01/07/2023]
Abstract
The approximate density-functional tight-binding theory method DFTB3 has been implemented in the quantum mechanics/molecular mechanics (QM/MM) framework of the Gromacs molecular simulation package. We show that the efficient smooth particle-mesh Ewald implementation of Gromacs extends to the calculation of QM/MM electrostatic interactions. Further, we make use of the various free-energy functionalities provided by Gromacs and the PLUMED plugin. We exploit the versatility and performance of the current framework in three typical applications of QM/MM methods to solve biophysical problems: (i) ultrafast proton transfer in malonaldehyde, (ii) conformation of the alanine dipeptide, and (iii) electron-induced repair of a DNA lesion. Also discussed is the further development of the framework, regarding mostly the options for parallelization.
Collapse
Affiliation(s)
- Tomáš Kubař
- Institute of Physical Chemistry and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, 76131, Germany
| | - Kai Welke
- Department of Chemistry, Nagoya University, Nagoya, 464-8602, Japan
| | - Gerrit Groenhof
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, Jyväskylä, 40014, Finland
| |
Collapse
|
17
|
Glowacki DR, Orr-Ewing AJ, Harvey JN. Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD3CN treated with a parallel multi-state EVB model. J Chem Phys 2015; 143:044120. [DOI: 10.1063/1.4926996] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David R. Glowacki
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
- PULSE Institute and Department of Chemistry, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | | | - Jeremy N. Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Heverlee, Belgium
| |
Collapse
|
18
|
Ghosh S, Talukder S, Sen S, Chaudhury P. Optimised polychromatic field-mediated suppression of H-atom tunnelling in a coupled symmetric double well: two-dimensional malonaldehyde model. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1068393] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
19
|
Review of computer simulations of isotope effects on biochemical reactions: From the Bigeleisen equation to Feynman's path integral. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2015; 1854:1782-94. [PMID: 25936775 DOI: 10.1016/j.bbapap.2015.04.021] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2015] [Revised: 04/13/2015] [Accepted: 04/22/2015] [Indexed: 12/18/2022]
Abstract
Enzymatic reactions are integral components in many biological functions and malfunctions. The iconic structure of each reaction path for elucidating the reaction mechanism in details is the molecular structure of the rate-limiting transition state (RLTS). But RLTS is very hard to get caught or to get visualized by experimentalists. In spite of the lack of explicit molecular structure of the RLTS in experiment, we still can trace out the RLTS unique "fingerprints" by measuring the isotope effects on the reaction rate. This set of "fingerprints" is considered as a most direct probe of RLTS. By contrast, for computer simulations, oftentimes molecular structures of a number of TS can be precisely visualized on computer screen, however, theoreticians are not sure which TS is the actual rate-limiting one. As a result, this is an excellent stage setting for a perfect "marriage" between experiment and theory for determining the structure of RLTS, along with the reaction mechanism, i.e., experimentalists are responsible for "fingerprinting", whereas theoreticians are responsible for providing candidates that match the "fingerprints". In this Review, the origin of isotope effects on a chemical reaction is discussed from the perspectives of classical and quantum worlds, respectively (e.g., the origins of the inverse kinetic isotope effects and all the equilibrium isotope effects are purely from quantum). The conventional Bigeleisen equation for isotope effect calculations, as well as its refined version in the framework of Feynman's path integral and Kleinert's variational perturbation (KP) theory for systematically incorporating anharmonicity and (non-parabolic) quantum tunneling, are also presented. In addition, the outstanding interplay between theory and experiment for successfully deducing the RLTS structures and the reaction mechanisms is demonstrated by applications on biochemical reactions, namely models of bacterial squalene-to-hopene polycyclization and RNA 2'-O-transphosphorylation. For all these applications, we used our recently-developed path-integral method based on the KP theory, called automated integration-free path-integral (AIF-PI) method, to perform ab initio path-integral calculations of isotope effects. As opposed to the conventional path-integral molecular dynamics (PIMD) and Monte Carlo (PIMC) simulations, values calculated from our AIF-PI path-integral method can be as precise as (not as accurate as) the numerical precision of the computing machine. Lastly, comments are made on the general challenges in theoretical modeling of candidates matching the experimental "fingerprints" of RLTS. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.
Collapse
|
20
|
Wang W, Zhao Y. The direct and precursor mediated dissociation rates of H2 on a Ni(111) surface. Phys Chem Chem Phys 2015; 17:5901-12. [PMID: 25630487 DOI: 10.1039/c4cp05624h] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The dissociation and recombination rates of physisorbed H2, and the direct and steady state dissociation (i.e., the precursor mediated dissociation) rates of gas phase H2 on Ni(111), as well as the corresponding kinetic isotope effects, are calculated using the quantum instanton method, together with path integral Monte Carlo and adaptive umbrella sampling techniques. All these rates except the recombination one first decrease and then increase with the increasing temperature, and their minimum values appear at about 250, 300 and 250 K, respectively. These non-monotonic behaviors reveal that the quantum effect of H2 should be very remarkable at low temperatures. The steady state rates are smaller than the direct rates at low temperatures, however, they become larger than the direct ones at high temperatures, these two kinds of rates become equal at about 400 and 300 K on the rigid and quantum lattices, respectively. The quantum motion of the lattice can enhance the direct and steady state rates, and it increases the steady state rate much more than the direct one, for instance, the direct and steady state rates on the quantum lattice are 1.30 and 2.08 times larger than that on the rigid one at 300 K. The calculated kinetic isotope effects are much larger than 1, which reveals that H2 always has a larger rate than that of D2, and the direct process predicts much larger kinetic isotope effects than the steady state process at low temperatures. In addition, the kinetic isotope effects are not affected by the lattice motion.
Collapse
Affiliation(s)
- Wenji Wang
- College of Science, Northwest A&F University, Yangling, 712100, Shaanxi Province, P. R. China.
| | | |
Collapse
|
21
|
Glover WJ, Casey JR, Schwartz BJ. Free Energies of Quantum Particles: The Coupled-Perturbed Quantum Umbrella Sampling Method. J Chem Theory Comput 2014; 10:4661-71. [DOI: 10.1021/ct500661t] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- William J. Glover
- Department of Chemistry and
Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Jennifer R. Casey
- Department of Chemistry and
Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Benjamin J. Schwartz
- Department of Chemistry and
Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| |
Collapse
|
22
|
Affiliation(s)
- Joshua P. Layfield
- Department of Chemistry, 600 South Mathews Avenue, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| | - Sharon Hammes-Schiffer
- Department of Chemistry, 600 South Mathews Avenue, University of Illinois at Urbana-Champaign, Urbana, IL 61801
| |
Collapse
|
23
|
Nagy T, Yosa Reyes J, Meuwly M. Multisurface Adiabatic Reactive Molecular Dynamics. J Chem Theory Comput 2014; 10:1366-75. [PMID: 26580356 DOI: 10.1021/ct400953f] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Adiabatic reactive molecular dynamics (ARMD) simulation method is a surface-crossing algorithm for modeling chemical reactions in classical molecular dynamics simulations using empirical force fields. As the ARMD Hamiltonian is time dependent during crossing, it allows only approximate energy conservation. In the current work, the range of applicability of conventional ARMD is explored, and a new multisurface ARMD (MS-ARMD) method is presented, implemented in CHARMM and applied to the vibrationally induced photodissociation of sulfuric acid (H2SO4) in the gas phase. For this, an accurate global potential energy surface (PES) involving 12 H2SO4 and 4 H2O + SO3 force fields fitted to MP2/6-311G++(2d,2p) reference energies is employed. The MS-ARMD simulations conserve total energy and feature both intramolecular H-transfer reactions and water elimination. An analytical treatment of the dynamics in the crossing region finds that conventional ARMD can approximately conserve total energy for limiting cases. In one of them, the reduced mass of the system is large, which often occurs for simulations of solvated biomolecular systems. On the other hand, MS-ARMD is a general approach for modeling chemical reactions including gas-phase, homogeneous, heterogeneous, and enzymatic catalytic reactions while conserving total energy in atomistic simulations.
Collapse
Affiliation(s)
- Tibor Nagy
- Department of Chemistry, University of Basel , 4056 Basel, Switzerland
| | | | - Markus Meuwly
- Department of Chemistry, University of Basel , 4056 Basel, Switzerland
| |
Collapse
|
24
|
Huang J, Buchowiecki M, Nagy T, Vaníček J, Meuwly M. Kinetic isotope effect in malonaldehyde determined from path integral Monte Carlo simulations. Phys Chem Chem Phys 2014; 16:204-11. [DOI: 10.1039/c3cp53698j] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
25
|
Wang W, Zhao Y. Dissociation rates of H2on a Ni(100) surface: the role of the physisorbed state. Phys Chem Chem Phys 2014; 16:13318-28. [DOI: 10.1039/c4cp01705f] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
26
|
Wang W, Zhao Y. Quantum instanton calculation of rate constant for CH4 + OH → CH3 + H2O reaction: torsional anharmonicity and kinetic isotope effect. J Chem Phys 2013; 137:214306. [PMID: 23231230 DOI: 10.1063/1.4768874] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Thermal rate constants for the title reaction are calculated by using the quantum instanton approximation within the full dimensional Cartesian coordinates. The results reveal that the quantum effect is remarkable for the reaction at both low and high temperatures, and the obtained rates are in good agreement with experimental measurements at high temperatures. Compared to the harmonic approximation, the torsional anharmonic effect of the internal rotation has a little influence on the rates at low temperatures, however, it enhances the rate by about 20% at 1000 K. In addition, the free energy barriers for the isotopic reactions and the temperature dependence of kinetic isotope effects are also investigated. Generally speaking, for the title reaction, the replacement of OH with OD will reduce the free energy barrier, while substituting D for H (connected to C) will increase the free energy barrier.
Collapse
Affiliation(s)
- Wenji Wang
- College of Science, Northwest A&F University, Shaanxi Province, Yangling, 712100, People's Republic of China.
| | | |
Collapse
|
27
|
Abstract
Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data directly from condensed phase ab initio simulations to develop reactive molecular dynamics models that do not require predefined empirical functions. Instead, the interactions used in the reactive model are expressed as linear combinations of interpolating functions that are optimized by using a linear least-squares algorithm. One notable benefit of the procedure outlined here is the capability to minimize the number of parameters requiring nonlinear optimization. The method presented can be generally applied to multiscale problems and is demonstrated by generating reactive models for the hydrated excess proton and hydroxide ion based directly on condensed phase ab initio molecular dynamics simulations. The resulting models faithfully reproduce the water-ion structural properties and diffusion constants from the ab initio simulations. Additionally, the free energy profiles for proton transfer, which is sensitive to the structural diffusion of both ions in water, are reproduced. The high fidelity of these models to ab initio simulations will permit accurate modeling of general chemical reactions in condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus facilitating a proper statistical sampling of the coupling to slow, large-scale motions of the system.
Collapse
Affiliation(s)
- Chris Knight
- Computing, Environment, and Life Sciences, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | | | | |
Collapse
|
28
|
Paasche A, Schirmeister T, Engels B. Benchmark Study for the Cysteine-Histidine Proton Transfer Reaction in a Protein Environment: Gas Phase, COSMO, QM/MM Approaches. J Chem Theory Comput 2013; 9:1765-77. [PMID: 26587634 DOI: 10.1021/ct301082y] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Proton transfer reactions are of crucial interest for the investigation of proteins. We have investigated the accuracy of commonly used quantum chemical methods for the description of proton transfer reactions in different environments (gas phase, COSMO, QM/MM) using the proton transfer between the catalytic dyad residues cysteine 145 and histidine 41 of SARS coronavirus main protease as a case study. The test includes thermodynamic, kinetic, and structural properties. The study comprises computationally demanding ab initio approaches (HF, CC2, MP2, SCS-CC2, SCS-MP2, CCSD(T)), popular density functional theories (BLYP, B3LYP, M06-2X), and semiempirical methods (MNDO/d, AM1, RM1, PM3, PM6). The approximated coupled cluster approach LCCSD(T) is taken as a reference method. We find that the robustness of the tested methods with respect to the environment correlates well with the level of theory. As an example HF, CC2, MP2, and their SCS variants show similar errors for gas phase, COSMO, or QM/MM computations. In contrast for semiempirical methods, the errors strongly diversify if one goes from gas phase to COSMO or QM/MM. Particular problems are observed for the recent semiempirical methods PM6 and RM1, which show the best performance for gas phase calculations but possess larger errors in conjunction with COSMO. Finally, a combination of SCS-MP2 and B3LYP or M06-2X allows reliable estimates about remaining errors.
Collapse
Affiliation(s)
- Alexander Paasche
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| | - Tanja Schirmeister
- Institut für Pharmazie und Biochemie, Johannes Gutenberg-Universität Mainz, Staudinger Weg 5, 55128 Mainz, Germany
| | - Bernd Engels
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
29
|
|
30
|
Quantum Instanton Evaluations of the Thermal Rate Constants for
Complex Systems. ADVANCES IN CHEMICAL PHYSICS 2012. [DOI: 10.1155/2012/483504] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Quantum instanton (QI) approximation is recently proposed for the evaluations of the chemical reaction rate constants with use of full dimensional potential energy surfaces. Its strategy is to use the instanton mechanism and to approximate time-dependent quantum dynamics to the imaginary time propagation of the quantities of partition function. It thus incorporates the properties of the instanton idea and the quantum effect of partition function and can be applied to chemical reactions of complex systems. In this paper, we present the QI approach and its applications to several complex systems mainly done by us. The concrete systems include, (1) the reaction of H+CH4→H2+CH3, (2) the reaction of H+SiH4→H2+SiH3, (3) H diffusion on Ni(100) surface; and (4) surface-subsurface transport and interior migration for H/Ni. Available experimental and other theoretical data are also presented for the purpose of comparison.
Collapse
|
31
|
Terranova ZL, Corcelli SA. Monitoring Intramolecular Proton Transfer with Two-Dimensional Infrared Spectroscopy: A Computational Prediction. J Phys Chem Lett 2012; 3:1842-1846. [PMID: 26291870 DOI: 10.1021/jz300714t] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Proton transfer processes are ubiquitous and play a vital role in a broad range of chemical and biochemical phenomena. The ability of two-dimensional infrared (2D IR) spectroscopy with a carbon-deuterium (C-D) reporter to monitor the kinetics of proton transfer in the model compound malonaldehyde was demonstrated computationally. One of the two carbonyl/enol carbon atoms in malonaldehyde was labeled with a C-D bond. The C-D stretch vibrational frequency provides ∼150 cm(-1) of sensitivity to the two tautomers of malonaldehyde. Mixed quantum mechanics/molecular mechanics simulations employing the self-consistent-charge density functional tight binding (SCC-DFTB) method were used to compute 2D IR line shapes for the C-D stretch of labeled malonaldehyde in aqueous solution. The 2D IR spectra reveal cross peaks from the chemical exchange of the proton. The kinetics for the growth of the cross-peaks (and the decay of the diagonal peaks) precisely match the proton transfer rate observed in the SCC-DFTB simulations.
Collapse
Affiliation(s)
- Z L Terranova
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - S A Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
32
|
Head-Gordon M. Publications. Mol Phys 2012. [DOI: 10.1080/00268976.2012.693252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
33
|
BUCHOWIECKI MARCIN. THE RATE CONSTANT OF THE O(3P) + HCl → OH + Cl REACTION WITHIN QUANTUM INSTANTON APPROXIMATION. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633612500095] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Temperature dependence of the rate constant of the O (3 P ) + HCl → OH + Cl reaction was calculated for 200–700 K range within the quantum instanton approximation using the 3A″ potential energy surface. Curvature of the Arrhenius plot arising from quantum effects in the hydrogen exchange reaction was observed. Results of the present method were found to significantly improve upon those obtained within the ICVT/μOMT method. The errors pertaining to the present method were estimated.
Collapse
Affiliation(s)
- MARCIN BUCHOWIECKI
- Institute of Physics, University of Szczecin, Wielkopolska 15, Szczecin, 70-451, Poland
| |
Collapse
|
34
|
Buchowiecki M. Quantum calculations of the temperature dependence of the rate constant and the equilibrium constant for the NH3+H⇌NH2+H2 reaction. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.01.075] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
35
|
Poma A, Monteferrante M, Bonella S, Ciccotti G. The quantum free energy barrier for hydrogen vacancy diffusion in Na3AlH6. Phys Chem Chem Phys 2012; 14:15458-63. [DOI: 10.1039/c2cp42536j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
36
|
Glowacki DR, Orr-Ewing AJ, Harvey JN. Product energy deposition of CN + alkane H abstraction reactions in gas and solution phases. J Chem Phys 2011; 134:214508. [DOI: 10.1063/1.3595259] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Stare J, Mavri J, Grdadolnik J, Zidar J, Maksić ZB, Vianello R. Hydrogen Bond Dynamics of Histamine Monocation in Aqueous Solution: Car–Parrinello Molecular Dynamics and Vibrational Spectroscopy Study. J Phys Chem B 2011; 115:5999-6010. [DOI: 10.1021/jp111175e] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Jernej Stare
- National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
| | - Janez Mavri
- National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
- EN→FIST Centre of Excellence, Dunajska 156, SI-1000 Ljubljana, Slovenia
| | - Jože Grdadolnik
- National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
- EN→FIST Centre of Excellence, Dunajska 156, SI-1000 Ljubljana, Slovenia
| | - Jernej Zidar
- National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
- EN→FIST Centre of Excellence, Dunajska 156, SI-1000 Ljubljana, Slovenia
| | | | - Robert Vianello
- National Institute of Chemistry, Hajdrihova 19, SI-1000 Ljubljana, Slovenia
- Ruđer Bošković Institute, Bijenička cesta 54, HR-10000 Zagreb, Croatia
| |
Collapse
|
38
|
Azuri A, Engel H, Doron D, Major DT. Path-Integral Calculations of Nuclear Quantum Effects in Model Systems, Small Molecules, and Enzymes via Gradient-Based Forward Corrector Algorithms. J Chem Theory Comput 2011; 7:1273-86. [DOI: 10.1021/ct100716c] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Asaf Azuri
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Hamutal Engel
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dvir Doron
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Dan Thomas Major
- Department of Chemistry and the Lise Meitner-Minerva Center of Computational Quantum Chemistry, Bar-Ilan University, Ramat-Gan 52900, Israel
| |
Collapse
|
39
|
Rommel JB, Goumans TPM, Kästner J. Locating Instantons in Many Degrees of Freedom. J Chem Theory Comput 2011; 7:690-8. [DOI: 10.1021/ct100658y] [Citation(s) in RCA: 99] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Judith B. Rommel
- Computational Biochemistry Group, Institute of Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - T. P. M. Goumans
- Gorlaeus Laboratories, LIC, Leiden University, Leiden, The Netherlands
| | - Johannes Kästner
- Computational Biochemistry Group, Institute of Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|