We present a method for carrying out long time scale dynamics simulations within the harmonic transition state theory approximation. For each state of the system, characterized by a local minimum on the potential energy surface, multiple searches for saddle points are carried out using random initial directions. The dimer method is used for the saddle point searches and the rate for each transition mechanism is estimated using harmonic transition state theory. Transitions are selected and the clock advanced according to the kinetic Monte Carlo algorithm. Unlike traditional applications of kinetic Monte Carlo, the atoms are not assumed to sit on lattice sites and a list of all possible transitions need not be specified beforehand. Rather, the relevant transitions are found on the fly during the simulation. A multiple time scale simulation of Al(100) crystal growth is presented where the deposition event, occurring on the time scale of picoseconds, is simulated by ordinary classical dynamics, but the time interval in between deposition events, on the order of milliseconds, is simulated by the long time scale algorithm. The Al(100) surface is found to grow remarkably smooth, even at 30 K because of concerted displacements of multiple atoms with significantly lower activation energy than adatom diffusion on the flat terrace.
Skip Nav Destination
Article navigation
1 December 2001
Research Article|
December 01 2001
Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table Available to Purchase
Graeme Henkelman;
Graeme Henkelman
Department of Chemistry, University of Washington, Seattle, Washington 98195-1700
Search for other works by this author on:
Hannes Jónsson
Hannes Jónsson
Department of Chemistry, University of Washington, Seattle, Washington 98195-1700
Faculty of Science, VR-II, University of Iceland, 107 Reykjavı́k, Iceland
Search for other works by this author on:
Graeme Henkelman
Department of Chemistry, University of Washington, Seattle, Washington 98195-1700
Hannes Jónsson
Department of Chemistry, University of Washington, Seattle, Washington 98195-1700
Faculty of Science, VR-II, University of Iceland, 107 Reykjavı́k, Iceland
J. Chem. Phys. 115, 9657–9666 (2001)
Article history
Received:
June 11 2001
Accepted:
September 13 2001
Connected Content
This is a companion to:
Optimization of hyperplanar transition states
Citation
Graeme Henkelman, Hannes Jónsson; Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table. J. Chem. Phys. 1 December 2001; 115 (21): 9657–9666. https://doi.org/10.1063/1.1415500
Download citation file:
Pay-Per-View Access
$40.00
Sign In
You could not be signed in. Please check your credentials and make sure you have an active account and try again.
Citing articles via
DeePMD-kit v2: A software package for deep potential models
Jinzhe Zeng, Duo Zhang, et al.
CREST—A program for the exploration of low-energy molecular chemical space
Philipp Pracht, Stefan Grimme, et al.
GPAW: An open Python package for electronic structure calculations
Jens Jørgen Mortensen, Ask Hjorth Larsen, et al.
Related Content
Novel predefined-time stability theory and its application in sliding mode control of synchronizing chaotic systems
Rev. Sci. Instrum. (December 2024)
Transition path sampling for non-equilibrium dynamics without predefined reaction coordinates
J. Chem. Phys. (January 2020)
Transition state theory: Variational formulation, dynamical corrections, and error estimates
J. Chem. Phys. (November 2005)
Numerical tests of coherence-corrected surface hopping methods using a donor-bridge-acceptor model system
J. Chem. Phys. (May 2019)
Probability distributions of molecular observables computed from Markov models
J. Chem. Phys. (June 2008)