Accelerated molecular dynamics with the bond-boost method

We present a new method for accelerating molecular-dynamics simulations of infrequent events. The method targets simulation of systems that spend most of the time in local energy minima, with slow transitions in between, as is the case with low-temperature surface diffusion. The potential-energy surface is modified by adding a boost potential in regions close to the local minima, such that all transition rates are increased while relative rates are preserved. The boost potential is an empirical function determined by the deviation of the bond lengths of a specified set of atoms from equilibrium. The method requires no previous knowledge of the processes involved and it can be applied to a wide variety of interaction potentials. Application to the diffusion of Cu atoms on the Cu(100) surface using an embedded-atom potential yields correct rates for adatom hopping, exchange, as well as vacancy and dimer diffusion with speed-ups up to several orders of magnitude.

[1]  John C. Tully,et al.  Molecular dynamics of infrequent events: Thermal desorption of xenon from a platinum surface , 1981 .

[2]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[3]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[4]  Foiles,et al.  Embedded-atom-method functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys. , 1986, Physical review. B, Condensed matter.

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

[6]  A. Voter.  A method for accelerating the molecular dynamics simulation of infrequent events , 1997 .

[7]  Laurent J. Lewis,et al.  Self-diffusion of adatoms, dimers, and vacancies on Cu(100) , 1997 .

[8]  A. Voter.  Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events , 1997 .

[9]  John W. Wilkins,et al.  Simple bias potential for boosting molecular dynamics with the hyperdynamics scheme , 1998 .

[10]  K. Fichthorn,et al.  Accelerated molecular dynamics of infrequent events , 1999 .

[11]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[12]  A. Voter,et al.  Temperature-accelerated dynamics for simulation of infrequent events , 2000 .

[13]  G. Henkelman,et al.  Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table , 2001 .

[14]  A. Voter,et al.  Closing the gap between experiment and theory: crystal growth by temperature accelerated dynamics. , 2001, Physical review letters.

[15]  Simulations of long time scale dynamics using the dimer method , 2001 .

[16]  A. Voter,et al.  Closing the Gap between Experiment and Theory , 2001 .

[17]  K. Fichthorn,et al.  Accelerated molecular dynamics of rare events using the local boost method , 2001 .

[18]  R. Miron,et al.  The Step and Slide method for finding saddle points on multidimensional potential surfaces , 2001 .

[19]  A. Voter,et al.  Extending the Time Scale in Atomistic Simulation of Materials Annual Re-views in Materials Research , 2002 .