Adaptive Kinetic Monte Carlo combines the simplicity of Kinetic Monte Carlo (KMC) with a saddle point search algorithm based on Molecular Dynamics (MD) in order to simulate metastable systems. Key to making Adaptive KMC effective is a stopping criterion for the saddle point search. In this work, we examine a criterion of S. T. Chill and G. Henkelman (J. Chem. Phys. 140 (2014), no. 21, 214110), which is based on the fraction of total reaction rate found instead of the fraction of observed saddles. The criterion uses the Eyring-Kramers law to estimate the reaction rate at the MD search temperature. We also consider a related criterion that remains valid when the Eyring-Kramers law is not. We examine the mathematical properties of both estimators and prove their mean square errors are well behaved, vanishing as the simulation continues to run.
ANALYSIS OF ESTIMATORS FOR ADAPTIVE KINETIC MONTE CARLO
Creators
David Aristoff - Colorado State University
Samuel T. Chill - QuantumWise A/S, Austin, TX 78749, United States
Gideon Simpson - Drexel University
Publication Details
Communications in applied mathematics and computational science, v 11(2), pp 171-186
Publisher
Mathematical Science Publ
Number of pages
16
Grant note
NSF-DMS-1522398 / National Science Foundation; National Science Foundation (NSF)
DE-SC0012733 / US Department of Energy; United States Department of Energy (DOE)
Resource Type
Journal article
Language
English
Academic Unit
Mathematics
Web of Science ID
WOS:000399134300002
Scopus ID
2-s2.0-85013229416
Other Identifier
991019168414004721
InCites Highlights
Data related to this publication, from InCites Benchmarking & Analytics tool: