Metadynamics of Paths.

We present a method to sample reactive pathways via biased molecular dynamics simulations in trajectory space. We show that the use of enhanced sampling techniques enables unconstrained exploration of multiple reaction routes. Time correlation functions are conveniently computed via reweighted averages along a single trajectory and kinetic rates are accessed at no additional cost. These abilities are illustrated analyzing a model potential and the umbrella inversion of NH_{3} in water. The algorithm allows a parallel implementation and promises to be a powerful tool for the study of rare events.

[1]  Massimiliano Bonomi,et al.  PLUMED 2: New feathers for an old bird , 2013, Comput. Phys. Commun..

[2]  Massimiliano Bonomi,et al.  Reconstructing the equilibrium Boltzmann distribution from well‐tempered metadynamics , 2009, J. Comput. Chem..

[3]  Finding multiple reaction pathways via global optimization of action , 2017, Nature communications.

[4]  A. V. van Duin,et al.  ReaxFF reactive force field development and applications for molecular dynamics simulations of ammonia borane dehydrogenation and combustion. , 2010, The journal of physical chemistry. A.

[5]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[6]  Michele Ceriotti,et al.  Iterative unbiasing of quasi-equilibrium sampling. , 2019, Journal of chemical theory and computation.

[7]  Peter G Bolhuis,et al.  Rate constants for diffusive processes by partial path sampling. , 2004, The Journal of chemical physics.

[8]  M. Parrinello,et al.  Metadynamics with Adaptive Gaussians. , 2012, Journal of chemical theory and computation.

[9]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[10]  Thijs J. H. Vlugt,et al.  On the efficient sampling of pathways in the transition path ensemble , 2001 .

[11]  Daniel Sebastiani,et al.  Generalized variational density functional perturbation theory , 2000 .

[12]  M. Parrinello,et al.  A time-independent free energy estimator for metadynamics. , 2015, The journal of physical chemistry. B.

[13]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[14]  Noam Bernstein,et al.  Exploration, Sampling, And Reconstruction of Free Energy Surfaces with Gaussian Process Regression. , 2016, Journal of chemical theory and computation.

[15]  Ron Elber,et al.  Calculation of classical trajectories with a very large time step: Formalism and numerical examples , 1996 .

[16]  M Parrinello,et al.  Action-derived molecular dynamics in the study of rare events. , 2001, Physical review letters.

[17]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[18]  Hiroshi Fujisaki,et al.  Onsager-Machlup action-based path sampling and its combination with replica exchange for diffusive and multiple pathways. , 2010, The Journal of chemical physics.

[19]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[20]  L Donati,et al.  Girsanov reweighting for path ensembles and Markov state models. , 2017, The Journal of chemical physics.

[21]  Lars Onsager,et al.  Fluctuations and Irreversible Processes , 1953 .

[22]  Ernesto E. Borrero,et al.  Avoiding traps in trajectory space: Metadynamics enhanced transition path sampling , 2016 .

[23]  M. Parrinello,et al.  Study of an F center in molten KCl , 1984 .

[24]  Gregory A. Voth,et al.  Hyper-parallel algorithms for centroid molecular dynamics: application to liquid para-hydrogen , 1996 .

[25]  Peter G. Bolhuis,et al.  A novel path sampling method for the calculation of rate constants , 2003 .

[26]  William H. Miller,et al.  Quantum mechanical transition state theory and a new semiclassical model for reaction rate constants , 1974 .

[27]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[28]  M. Parrinello,et al.  From metadynamics to dynamics. , 2013, Physical review letters.

[29]  Jörg Behler,et al.  High order path integrals made easy. , 2016, The Journal of chemical physics.

[30]  William H. Miller,et al.  Quantum mechanical rate constants for bimolecular reactions , 1983 .

[31]  Michele Parrinello,et al.  Well-tempered metadynamics converges asymptotically. , 2014, Physical review letters.

[32]  Ron Elber,et al.  A method for determining reaction paths in large molecules: application to myoglobin , 1987 .

[33]  Matteo Salvalaglio,et al.  Time-independent free energies from metadynamics via mean force integration. , 2019, The Journal of chemical physics.

[34]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[35]  Bin Zhang,et al.  Path-accelerated stochastic molecular dynamics: Parallel-in-time integration using path integrals. , 2019, The Journal of chemical physics.

[36]  Michele Parrinello,et al.  Rethinking Metadynamics: from Bias Potentials to Probability Distributions. , 2020, The journal of physical chemistry letters.

[37]  Michele Parrinello,et al.  Enhancing Important Fluctuations: Rare Events and Metadynamics from a Conceptual Viewpoint. , 2016, Annual review of physical chemistry.

[38]  D. Chandler,et al.  Introduction To Modern Statistical Mechanics , 1987 .

[39]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

[40]  Peter G. Bolhuis,et al.  Nested Transition Path Sampling. , 2018, Physical review letters.

[41]  Michele Parrinello,et al.  Variational Flooding Study of a SN2 Reaction. , 2017, The journal of physical chemistry letters.

[42]  B. Keller,et al.  Girsanov reweighting for metadynamics simulations. , 2018, The Journal of chemical physics.