A method for accelerating the molecular dynamics simulation of infrequent events

For infrequent-event systems, transition state theory (TST) is a powerful approach for overcoming the time scale limitations of the molecular dynamics (MD) simulation method, provided one knows the locations of the potential-energy basins (states) and the TST dividing surfaces (or the saddle points) between them. Often, however, the states to which the system will evolve are not known in advance. We present a new, TST-based method for extending the MD time scale that does not require advanced knowledge of the states of the system or the transition states that separate them. The potential is augmented by a bias potential, designed to raise the energy in regions other than at the dividing surfaces. State to state evolution on the biased potential occurs in the proper sequence, but at an accelerated rate with a nonlinear time scale. Time is no longer an independent variable, but becomes a statistically estimated property that converges to the exact result at long times. The long-time dynamical behavior is ex...

[1]  M. R. Marcelin,et al.  Contribution à l'étude de la cinétique physico-chimique , 1915 .

[2]  H. Eyring.  The Activated Complex in Chemical Reactions , 1935 .

[3]  Juro Horiuti,et al.  On the Statistical Mechanical Treatment of the Absolute Rate of Chemical Reaction , 1938 .

[4]  H. Kramers.  Brownian motion in a field of force and the diffusion model of chemical reactions , 1940 .

[5]  S. Chandrasekhar.  Stochastic problems in Physics and Astronomy , 1943 .

[6]  G. Vineyard.  Frequency factors and isotope effects in solid state rate processes , 1957 .

[7]  James C. Keck,et al.  Statistical investigation of dissociation cross-sections for diatoms , 1962 .

[8]  A. Nowick,et al.  Diffusion in solids: recent developments , 1975 .

[9]  Ralph E. Christoffersen,et al.  Algorithms for Chemical Computations , 1977 .

[10]  John A. Montgomery,et al.  Trajectory analysis of a kinetic theory for isomerization dynamics in condensed phases , 1979 .

[11]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[12]  B. C. Garrett,et al.  Variational Transition State Theory , 1980 .

[13]  J. Doll.  A unified theory of dissociation , 1980 .

[14]  Surface self‐diffusion constants at low temperature: Monte Carlo transition state theory with importance sampling , 1984 .

[15]  M. Baskes,et al.  Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals , 1984 .

[16]  Arthur F. Voter,et al.  Dynamical corrections to transition state theory for multistate systems: Surface self-diffusion in the rare-event regime , 1985 .

[17]  Graeme E. Murch,et al.  Diffusion in solids : recent developments , 1985 .

[18]  Flynn,et al.  Theory of classical diffusion jumps in solids. , 1985, Physical review. B, Condensed matter.

[19]  A. Voter.  A Monte Carlo method for determining free‐energy differences and transition state theory rate constants , 1985 .

[20]  D. C. Clary,et al.  The Theory of Chemical Reaction Dynamics , 1986 .

[21]  A. Voter,et al.  Classically exact overlayer dynamics: Diffusion of rhodium clusters on Rh(100). , 1986, Physical review. B, Condensed matter.

[22]  Arthur F. Voter,et al.  Accurate Interatomic Potentials for Ni, Al and Ni3Al , 1986 .

[23]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[24]  John E. Straub,et al.  Classical and modern methods in reaction rate theory , 1988 .

[25]  Voter.  Computing classically exact diffusion constants using short-time trajectories. , 1989, Physical review letters.

[26]  J. Doll,et al.  Using multistate dynamical corrections to compute classically exact diffusion constants at arbitrary temperature , 1989 .

[27]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[28]  Kellogg,et al.  Surface self-diffusion on Pt(001) by an atomic exchange mechanism. , 1990, Physical review letters.

[29]  Mark E. Tuckerman,et al.  Molecular dynamics algorithm for condensed systems with multiple time scales , 1990 .

[30]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[31]  Chen,et al.  Displacement distribution and atomic jump direction in diffusion of Ir atoms on the Ir(001) surface. , 1990, Physical review letters.

[32]  P. Feibelman,et al.  Diffusion path for an Al adatom on Al(001). , 1990, Physical review letters.

[33]  J. B. Adams,et al.  EAM study of surface self-diffusion of single adatoms of fcc metals Ni, Cu, Al, Ag, Au, Pd, and Pt , 1991 .

[34]  B. C. Garrett,et al.  Centroid‐density quantum rate theory: Dynamical treatment of classical recrossing , 1993 .

[35]  E. Sevick,et al.  A chain of states method for investigating infrequent event processes occurring in multistate, multidimensional systems , 1993 .

[36]  Murray S. Daw,et al.  The embedded-atom method: a review of theory and applications , 1993 .

[37]  Abhinandan Jain,et al.  Protein simulations using techniques suitable for very large systems: The cell multipole method for nonbond interactions and the Newton‐Euler inverse mass operator method for internal coordinate dynamics , 1994, Proteins.

[38]  Sebastian Doniach,et al.  Long‐time overdamped Langevin dynamics of molecular chains , 1994, J. Comput. Chem..

[39]  Jónsson,et al.  Low-temperature homoepitaxial growth of Pt(111) in simulated vapor deposition. , 1994, Physical review. B, Condensed matter.

[40]  R. Stephen Berry,et al.  Statistical interpretation of topographies and dynamics of multidimensional potentials , 1995 .

[41]  G. Huber,et al.  Weighted-ensemble Brownian dynamics simulations for protein association reactions. , 1996, Biophysical journal.

[42]  Pawan Kumar,et al.  Smart Monte Carlo for accurate simulation of rare‐event dynamics: Diffusion of adsorbed species on solid surfaces , 1996 .

[43]  M. Nadeau.  Proteins : Structure , Function , and Genetics , .