1
|
Bracher CE, Allen CJ, Singleton DA. Nuclear Quantum Effects on the Nature of Hydroboration Selectivity: Experimental Effects of First-Collision Tunneling. J Am Chem Soc 2024; 146:25907-25911. [PMID: 39284010 PMCID: PMC11440546 DOI: 10.1021/jacs.4c09306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/26/2024]
Abstract
The understanding of selectivity in reactions exhibiting nonstatistical dynamics is impeded by the limitations of trajectory studies with regard to nuclear quantum effects, especially tunneling. We described here the use of ring-polymer molecular dynamics (RPMD) to account for an unusual regiochemical isotope effect on the regioselectivity of hydroborations of alkenes with BH3/BD3. RPMD is able to account for the experimental observation, while statistical approaches and classical trajectories fail. The combination of experiment and RPMD trajectories suggests that tunneling in the initial collision of reactants is the major source of the nonstatistical selectivity.
Collapse
Affiliation(s)
- Christoph E Bracher
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| | - Connor J Allen
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| | - Daniel A Singleton
- Department of Chemistry, Texas A&M University, P.O. Box 30012, College Station, Texas 77842, United States
| |
Collapse
|
2
|
Ojha AA, Thakur S, Ahn SH, Amaro RE. DeepWEST: Deep Learning of Kinetic Models with the Weighted Ensemble Simulation Toolkit for Enhanced Sampling. J Chem Theory Comput 2023; 19:1342-1359. [PMID: 36719802 DOI: 10.1021/acs.jctc.2c00282] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Recent advances in computational power and algorithms have enabled molecular dynamics (MD) simulations to reach greater time scales. However, for observing conformational transitions associated with biomolecular processes, MD simulations still have limitations. Several enhanced sampling techniques seek to address this challenge, including the weighted ensemble (WE) method, which samples transitions between metastable states using many weighted trajectories to estimate kinetic rate constants. However, initial sampling of the potential energy surface has a significant impact on the performance of WE, i.e., convergence and efficiency. We therefore introduce deep-learned kinetic modeling approaches that extract statistically relevant information from short MD trajectories to provide a well-sampled initial state distribution for WE simulations. This hybrid approach overcomes any statistical bias to the system, as it runs short unbiased MD trajectories and identifies meaningful metastable states of the system. It is shown to provide a more refined free energy landscape closer to the steady state that could efficiently sample kinetic properties such as rate constants.
Collapse
Affiliation(s)
- Anupam Anand Ojha
- Department of Chemistry, University of California San Diego, La Jolla, California92093, United States
| | - Saumya Thakur
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, Maharashtra400076, India
| | - Surl-Hee Ahn
- Department of Chemical Engineering, University of California Davis, Davis, California95616, United States
| | - Rommie E Amaro
- Department of Chemistry, University of California San Diego, La Jolla, California92093, United States
| |
Collapse
|
3
|
Abstract
The kinetics of a dynamical system dominated by two metastable states is examined from the perspective of the activated-dynamics reactive flux formalism, Markov state eigenvalue spectral decomposition, and committor-based transition path theory. Analysis shows that the different theoretical formulations are consistent, clarifying the significance of the inherent microscopic lag-times that are implicated, and that the most meaningful one-dimensional reaction coordinate in the region of the transition state is along the gradient of the committor in the multidimensional subspace of collective variables. It is shown that the familiar reactive flux activated dynamics formalism provides an effective route to calculate the transition rate in the case of a narrow sharp barrier but much less so in the case of a broad flat barrier. In this case, the standard reactive flux correlation function decays very slowly to the plateau value that corresponds to the transmission coefficient. Treating the committor function as a reaction coordinate does not alleviate all issues caused by the slow relaxation of the reactive flux correlation function. A more efficient activated dynamics simulation algorithm may be achieved from a modified reactive flux weighted by the committor. Simulation results on simple systems are used to illustrate the various conceptual points.
Collapse
Affiliation(s)
- Benoît Roux
- Department of Biochemistry and Molecular Biology, Department of Chemistry, The University of Chicago, 5735 S Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
4
|
Garcia-Meseguer R, Carpenter BK. Re-Evaluating the Transition State for Reactions in Solution. European J Org Chem 2018. [DOI: 10.1002/ejoc.201800841] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Barry K. Carpenter
- School of Chemistry; Cardiff University; CF10 3AT Cardiff United Kingdom
| |
Collapse
|
5
|
Shannon RJ, Amabilino S, O’Connor M, Shalishilin DV, Glowacki DR. Adaptively Accelerating Reactive Molecular Dynamics Using Boxed Molecular Dynamics in Energy Space. J Chem Theory Comput 2018; 14:4541-4552. [DOI: 10.1021/acs.jctc.8b00515] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Robin J. Shannon
- Mechanical Engineering, Stanford University, Stanford, California 94305, United States
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Silvia Amabilino
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Mike O’Connor
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, U.K
| | | | - David R. Glowacki
- Mechanical Engineering, Stanford University, Stanford, California 94305, United States
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
- Department of Computer Science, University of Bristol, Bristol BS8 1UB, U.K
| |
Collapse
|
6
|
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
|
7
|
Martins‐Costa MTC, Anglada JM, Ruiz‐López MF. Computational Insights into the CH
3
Cl+OH Chemical Reaction Dynamics at the Air–Water Interface. Chemphyschem 2017; 18:2747-2755. [DOI: 10.1002/cphc.201700437] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2017] [Revised: 06/02/2017] [Indexed: 11/11/2022]
Affiliation(s)
| | - Josep M. Anglada
- Departament de Química Biològica i Modelització MolecularIQAC, CSIC c/ Jordi Girona 18 E-08034 Barcelona Spain
| | - Manuel F. Ruiz‐López
- SRSMCUniversity of Lorraine and CNRS Boulevard des Aiguillettes, BP 70239 54506 Vandoeuvre-lès-Nancy France
| |
Collapse
|
8
|
Warshel A, Bora RP. Perspective: Defining and quantifying the role of dynamics in enzyme catalysis. J Chem Phys 2017; 144:180901. [PMID: 27179464 DOI: 10.1063/1.4947037] [Citation(s) in RCA: 142] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Enzymes control chemical reactions that are key to life processes, and allow them to take place on the time scale needed for synchronization between the relevant reaction cycles. In addition to general interest in their biological roles, these proteins present a fundamental scientific puzzle, since the origin of their tremendous catalytic power is still unclear. While many different hypotheses have been put forward to rationalize this, one of the proposals that has become particularly popular in recent years is the idea that dynamical effects contribute to catalysis. Here, we present a critical review of the dynamical idea, considering all reasonable definitions of what does and does not qualify as a dynamical effect. We demonstrate that no dynamical effect (according to these definitions) has ever been experimentally shown to contribute to catalysis. Furthermore, the existence of non-negligible dynamical contributions to catalysis is not supported by consistent theoretical studies. Our review is aimed, in part, at readers with a background in chemical physics and biophysics, and illustrates that despite a substantial body of experimental effort, there has not yet been any study that consistently established a connection between an enzyme's conformational dynamics and a significant increase in the catalytic contribution of the chemical step. We also make the point that the dynamical proposal is not a semantic issue but a well-defined scientific hypothesis with well-defined conclusions.
Collapse
Affiliation(s)
- Arieh Warshel
- Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, USA
| | - Ram Prasad Bora
- Department of Chemistry, University of Southern California, SGM 418, 3620 McClintock Avenue, Los Angeles, California 90089, USA
| |
Collapse
|
9
|
Dahlen O, van Erp TS. Mesoscopic modeling of DNA denaturation rates: Sequence dependence and experimental comparison. J Chem Phys 2015; 142:235101. [DOI: 10.1063/1.4922519] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Oda Dahlen
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, Realfagbygget D3-117 7491 Trondheim, Norway
| | - Titus S. van Erp
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), Høgskoleringen 5, Realfagbygget D3-117 7491 Trondheim, Norway
| |
Collapse
|
10
|
Klenin KV. Efficient calculation of rate constants: Downhill versus uphill sampling. J Chem Phys 2014; 141:074103. [DOI: 10.1063/1.4892565] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Divi S, Chatterjee A. Accelerating rare events while overcoming the low-barrier problem using a temperature program. J Chem Phys 2014; 140:184115. [DOI: 10.1063/1.4875476] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Laner M, Horta BA, Hünenberger PH. Phase-transition properties of glycerol-monopalmitate lipid bilayers investigated by molecular dynamics simulation: influence of the system size and force-field parameters. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2012.755526] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
13
|
Mauguière FAL, Collins P, Ezra GS, Wiggins S. Bond breaking in a Morse chain under tension: Fragmentation patterns, higher index saddles, and bond healing. J Chem Phys 2013; 138:134118. [DOI: 10.1063/1.4798641] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Li Y, Suleimanov YV, Li J, Green WH, Guo H. Rate coefficients and kinetic isotope effects of the X + CH4 → CH3 + HX (X = H, D, Mu) reactions from ring polymer molecular dynamics. J Chem Phys 2013; 138:094307. [DOI: 10.1063/1.4793394] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Meadley SL, Escobedo FA. Thermodynamics and kinetics of bubble nucleation: Simulation methodology. J Chem Phys 2012; 137:074109. [DOI: 10.1063/1.4745082] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
WHITE GEORGEWN, GOLDMAN SAUL, GRAY CG. Test of rate theory transmission coefficient algorithms. An application to ion channels. Mol Phys 2012. [DOI: 10.1080/00268970009483391] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
17
|
Van Erp TS. Dynamical Rare Event Simulation Techniques for Equilibrium and Nonequilibrium Systems. ADVANCES IN CHEMICAL PHYSICS 2012. [DOI: 10.1002/9781118309513.ch2] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
18
|
Lau D, Lam R. Atomistic Prediction of Nanomaterials: Introduction to Molecular Dynamics Simulation and a Case Study of Graphene Wettability. IEEE NANOTECHNOLOGY MAGAZINE 2012. [DOI: 10.1109/mnano.2011.2181739] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
19
|
Zheng J, Truhlar DG. Multi-path variational transition state theory for chemical reaction rates of complex polyatomic species: ethanol + OH reactions. Faraday Discuss 2012; 157:59-88; discussion 113-40. [DOI: 10.1039/c2fd20012k] [Citation(s) in RCA: 105] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Anderson JB, Anderson LE, Kussmann J. Monte Carlo simulations of single- and multistep enzyme-catalyzed reaction sequences: effects of diffusion, cell size, enzyme fluctuations, colocalization, and segregation. J Chem Phys 2010; 133:034104. [PMID: 20649305 DOI: 10.1063/1.3459111] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Following the discovery of slow fluctuations in the catalytic activity of an enzyme in single-molecule experiments, it has been shown that the classical Michaelis-Menten (MM) equation relating the average enzymatic velocity and the substrate concentration may hold even for slowly fluctuating enzymes. In many cases, the average velocity is that given by the MM equation with time-averaged values of the fluctuating rate constants and the effect of enzyme fluctuations is simply averaged out. The situation is quite different for a sequence of reactions. For colocalization of a pair of enzymes in a sequence to be effective in promoting reaction, the second must be active when the first is active or soon after. If the enzymes are slowly varying and only rarely active, the product of the first reaction may diffuse away before the second enzyme is active, and colocalization may have little value. Even for single-step reactions the interplay of reaction and diffusion with enzyme fluctuations leads to added complexities, but for multistep reactions the interplay of reaction and diffusion, cell size, compartmentalization, enzyme fluctuations, colocalization, and segregation is far more complex than for single-step reactions. In this paper, we report the use of stochastic simulations at the level of whole cells to explore, understand, and predict the behavior of single- and multistep enzyme-catalyzed reaction systems exhibiting some of these complexities. Results for single-step reactions confirm several earlier observations by others. The MM relationship, with altered constants, is found to hold for single-step reactions slowed by diffusion. For single-step reactions, the distribution of enzymes in a regular grid is slightly more effective than a random distribution. Fluctuations of enzyme activity, with average activity fixed, have no observed effects for simple single-step reactions slowed by diffusion. Two-step sequential reactions are seen to be slowed by segregation of the enzymes for each step, and results of the calculations suggest limits for cell size. Colocalization of enzymes for a two-step sequence is seen to promote reaction, and rates fall rapidly with increasing distance between enzymes. Low frequency fluctuations of the activities of colocalized enzymes, with average activities fixed, can greatly reduce reaction rates for sequential reactions.
Collapse
Affiliation(s)
- James B Anderson
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, USA.
| | | | | |
Collapse
|
21
|
Collins P, Ezra GS, Wiggins S. Phase space structure and dynamics for the Hamiltonian isokinetic thermostat. J Chem Phys 2010; 133:014105. [DOI: 10.1063/1.3455712] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Allen RJ, Valeriani C, Rein Ten Wolde P. Forward flux sampling for rare event simulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2009; 21:463102. [PMID: 21715864 DOI: 10.1088/0953-8984/21/46/463102] [Citation(s) in RCA: 229] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Rare events are ubiquitous in many different fields, yet they are notoriously difficult to simulate because few, if any, events are observed in a conventional simulation run. Over the past several decades, specialized simulation methods have been developed to overcome this problem. We review one recently developed class of such methods, known as forward flux sampling. Forward flux sampling uses a series of interfaces between the initial and final states to calculate rate constants and generate transition paths for rare events in equilibrium or nonequilibrium systems with stochastic dynamics. This review draws together a number of recent advances, summarizes several applications of the method and highlights challenges that remain to be overcome.
Collapse
Affiliation(s)
- Rosalind J Allen
- SUPA, School of Physics and Astronomy, University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, UK
| | | | | |
Collapse
|
23
|
Hernández-Ortiz JP, Chopra M, Geier S, de Pablo JJ. Hydrodynamic effects on the translocation rate of a polymer through a pore. J Chem Phys 2009; 131:044904. [DOI: 10.1063/1.3184798] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Ezra GS, Waalkens H, Wiggins S. Microcanonical rates, gap times, and phase space dividing surfaces. J Chem Phys 2009; 130:164118. [DOI: 10.1063/1.3119365] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Olsson MHM, Parson WW, Warshel A. Dynamical contributions to enzyme catalysis: critical tests of a popular hypothesis. Chem Rev 2007; 106:1737-56. [PMID: 16683752 DOI: 10.1021/cr040427e] [Citation(s) in RCA: 254] [Impact Index Per Article: 14.9] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mats H M Olsson
- Department of Chemistry, University of Southern California, 3620 McClintock Avenue, Los Angeles, California 90089-1062, USA.
| | | | | |
Collapse
|
26
|
Borrero EE, Escobedo FA. Reaction coordinates and transition pathways of rare events via forward flux sampling. J Chem Phys 2007; 127:164101. [PMID: 17979313 DOI: 10.1063/1.2776270] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ernesto E Borrero
- School of Chemical and Biomolecular Engineering, Cornell University, Ithaca, New York 14853, USA
| | | |
Collapse
|
27
|
van Erp TS. Efficiency analysis of reaction rate calculation methods using analytical models I: The two-dimensional sharp barrier. J Chem Phys 2006; 125:174106. [PMID: 17100428 DOI: 10.1063/1.2363996] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We analyze the efficiency of different methods for the calculation of reaction rates in the case of a simple two-dimensional analytical benchmark system. Two classes of methods are considered: the first is based on the free energy calculation along a reaction coordinate and the calculation of the transmission coefficient, the second on the sampling of dynamical pathways. We give scaling rules for how this efficiency depends on barrier height and width, and we hand out simple optimization rules for the method-specific parameters. We show that the path sampling methods, using the transition interface sampling technique, become exceedingly more efficient than the others when the reaction coordinate is not the optimal one.
Collapse
Affiliation(s)
- Titus S van Erp
- Centrum voor Oppervlaktechemie en Katalyse, K.U. Leuven, Kasteelpark Arenberg 23, B-3001 Leuven, Belgium.
| |
Collapse
|
28
|
Nangia S, Jasper AW, Miller TF, Truhlar DG. Army ants algorithm for rare event sampling of delocalized nonadiabatic transitions by trajectory surface hopping and the estimation of sampling errors by the bootstrap method. J Chem Phys 2006; 120:3586-97. [PMID: 15268520 DOI: 10.1063/1.1641019] [Citation(s) in RCA: 64] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The most widely used algorithm for Monte Carlo sampling of electronic transitions in trajectory surface hopping (TSH) calculations is the so-called anteater algorithm, which is inefficient for sampling low-probability nonadiabatic events. We present a new sampling scheme (called the army ants algorithm) for carrying out TSH calculations that is applicable to systems with any strength of coupling. The army ants algorithm is a form of rare event sampling whose efficiency is controlled by an input parameter. By choosing a suitable value of the input parameter the army ants algorithm can be reduced to the anteater algorithm (which is efficient for strongly coupled cases), and by optimizing the parameter the army ants algorithm may be efficiently applied to systems with low-probability events. To demonstrate the efficiency of the army ants algorithm, we performed atom-diatom scattering calculations on a model system involving weakly coupled electronic states. Fully converged quantum mechanical calculations were performed, and the probabilities for nonadiabatic reaction and nonreactive deexcitation (quenching) were found to be on the order of 10(-8). For such low-probability events the anteater sampling scheme requires a large number of trajectories ( approximately 10(10)) to obtain good statistics and converged semiclassical results. In contrast by using the new army ants algorithm converged results were obtained by running 10(5) trajectories. Furthermore, the results were found to be in excellent agreement with the quantum mechanical results. Sampling errors were estimated using the bootstrap method, which is validated for use with the army ants algorithm.
Collapse
Affiliation(s)
- Shikha Nangia
- Department of Chemistry, University of Minnesota, Minneapolis 55455-0431, USA.
| | | | | | | |
Collapse
|
29
|
Chen LY, Horing NJM. Transition rate prefactors for systems of many degrees of freedom. J Chem Phys 2006; 124:164102. [PMID: 16674124 DOI: 10.1063/1.2188943] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
When a minimum on the potential energy surface is surrounded by multiple saddle points with similar energy barriers, the transition pathways with greater prefactors are more important than those that have similar energy barriers but smaller prefactors. In this paper, we present a theoretical formulation for the prefactors, computing the probabilities for transition paths from a minimum to its surrounding saddle points. We apply this formulation to a system of 2 degrees of freedom and a system of 14 degrees of freedom. The first is Brownian motion in a two-dimensional potential whose global anharmonicities play a dominant role in determining the transition rates. The second is a Lennard-Jones (LJ) cluster of seven particles in two dimensions. Low lying transition states of the LJ cluster, which can be reached directly from a minimum without passing through another minimum, are identified without any presumption of their characteristics nor of the product states they lead to. The probabilities are computed for paths going from an equilibrium ensemble of states near a given minimum to the surrounding transition states. These probabilities are directly related to the prefactors in the rate formula. This determination of the rate prefactors includes all anharmonicities, near or far from transition states, which are pertinent in the very sophisticated energy landscape of LJ clusters and in many other complex systems.
Collapse
Affiliation(s)
- L Y Chen
- Department of Physics, University of Texas at San Antonio, San Antonio, Texas 78249-0697, USA.
| | | |
Collapse
|
30
|
Vanden-Eijnden E, Tal FA. Transition state theory: Variational formulation, dynamical corrections, and error estimates. J Chem Phys 2005; 123:184103. [PMID: 16292895 DOI: 10.1063/1.2102898] [Citation(s) in RCA: 104] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Transition state theory (TST) is revisited, as well as evolutions upon TST such as variational TST in which the TST dividing surface is optimized so as to minimize the rate of recrossing through this surface and methods which aim at computing dynamical corrections to the TST transition rate constant. The theory is discussed from an original viewpoint. It is shown how to compute exactly the mean frequency of transition between two predefined sets which either partition phase space (as in TST) or are taken to be well-separated metastable sets corresponding to long-lived conformation states (as necessary to obtain the actual transition rate constants between these states). Exact and approximate criterions for the optimal TST dividing surface with minimum recrossing rate are derived. Some issues about the definition and meaning of the free energy in the context of TST are also discussed. Finally precise error estimates for the numerical procedure to evaluate the transmission coefficient kappaS of the TST dividing surface are given, and it is shown that the relative error on kappaS scales as 1/square root(kappaS) when kappaS is small. This implies that dynamical corrections to the TST rate constant can be computed efficiently if and only if the TST dividing surface has a transmission coefficient kappaS which is not too small. In particular, the TST dividing surface must be optimized upon (for otherwise kappaS is generally very small), but this may not be sufficient to make the procedure numerically efficient (because the optimal dividing surface has maximum kappaS, but this coefficient may still be very small).
Collapse
Affiliation(s)
- Eric Vanden-Eijnden
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA.
| | | |
Collapse
|
31
|
Zhao Y, Mil'nikov G, Nakamura H. Evaluation of canonical and microcanonical nonadiabatic reaction rate constants by using the Zhu–Nakamura formulas. J Chem Phys 2004; 121:8854-60. [PMID: 15527348 DOI: 10.1063/1.1801971] [Citation(s) in RCA: 34] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We consider a problem of calculating both thermal and microcanonical rate constants for nonadiabatic chemical reactions. Instead of using the conventional transition state theory, we use a generalized seam surface and introduce a concept of a coordinate dependent effective nonadiabatic transition probability based on the Zhu-Nakamura theory which can treat the nonadiabatic tunneling properly. The present approach can be combined with Monte Carlo method so as to be applicable to chemical reactions in complicated systems. The method is demonstrated to work well in wide energy and temperature range. Numerical tests also show that it is very essential for accurate evaluation of the thermal rate constant to use the generalized seam surface and take into account the nonadiabatic tunneling effect.
Collapse
Affiliation(s)
- Yi Zhao
- Department of Theoretical Studies, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| | | | | |
Collapse
|
32
|
|
33
|
Schenter GK, Garrett BC, Truhlar DG. Generalized transition state theory in terms of the potential of mean force. J Chem Phys 2003. [DOI: 10.1063/1.1597477] [Citation(s) in RCA: 98] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|
34
|
TOLOKH IGORS, WHITE GEORGEWN, GOLDMAN SAUL, GRAY CG. Prediction of ion channel transport from Grote—Hynes and Kramers theories. Mol Phys 2002. [DOI: 10.1080/00268970210124828] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
35
|
Bolhuis PG, Chandler D, Dellago C, Geissler PL. Transition path sampling: throwing ropes over rough mountain passes, in the dark. Annu Rev Phys Chem 2002; 53:291-318. [PMID: 11972010 DOI: 10.1146/annurev.physchem.53.082301.113146] [Citation(s) in RCA: 1305] [Impact Index Per Article: 59.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This article reviews the concepts and methods of transition path sampling. These methods allow computational studies of rare events without requiring prior knowledge of mechanisms, reaction coordinates, and transition states. Based upon a statistical mechanics of trajectory space, they provide a perspective with which time dependent phenomena, even for systems driven far from equilibrium, can be examined with the same types of importance sampling tools that in the past have been applied so successfully to static equilibrium properties.
Collapse
Affiliation(s)
- Peter G Bolhuis
- Department of Chemical Engineering, Nieuwe Achtergracht 166, 1018 WV Amsterdam, The Netherlands.
| | | | | | | |
Collapse
|
36
|
HU JIANGUO, GOLDMAN SAUL, GRAY CG, GUY HROBERT. Calculation of the conductance and selectivity of an ion-selective potassium channel (IRK1) from simulation of atomic scale models. Mol Phys 2000. [DOI: 10.1080/00268970009483320] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
37
|
Li GS, Martins Costa MT, Millot C, Ruiz-López MF. Effect of solvent fluctuations on proton transfer dynamics: a hybrid AM1/MM molecular dynamics simulation on the [H3N–H–NH3]+ system. Chem Phys 1999. [DOI: 10.1016/s0301-0104(98)00370-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
38
|
den Otter WK, Briels WJ. The reactive flux method applied to complex isomerization reactions: Using the unstable normal mode as a reaction coordinate. J Chem Phys 1997. [DOI: 10.1063/1.473573] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Voter AF. A method for accelerating the molecular dynamics simulation of infrequent events. J Chem Phys 1997. [DOI: 10.1063/1.473503] [Citation(s) in RCA: 505] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
40
|
Strnad M, Martins-Costa MTC, Millot C, Tuñón I, Ruiz-López MF, Rivail JL. Molecular dynamics simulations of elementary chemical processes in liquid water using combined density functional and molecular mechanics potentials. II. Charge separation processes. J Chem Phys 1997. [DOI: 10.1063/1.473458] [Citation(s) in RCA: 59] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Neria E, Karplus M. A position dependent friction model for solution reactions in the high friction regime: Proton transfer in triosephosphate isomerase (TIM). J Chem Phys 1996. [DOI: 10.1063/1.472889] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
42
|
Neria E, Fischer S, Karplus M. Simulation of activation free energies in molecular systems. J Chem Phys 1996. [DOI: 10.1063/1.472061] [Citation(s) in RCA: 1020] [Impact Index Per Article: 36.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|