SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations
Abstract
Lattice based Kinetic Monte Carlo (KMC) simulations offer a powerful simulation technique for investigating large reaction networks while retaining spatial configuration information, unlike ordinary differential equations. However, large chemical reaction networks can contain reaction processes with rates spanning multiple orders of magnitude. This can lead to the problem of “KMC stiffness” (similar to stiffness in differential equations), where the computational expense has the potential to be overwhelmed by very short timesteps during KMC simulations, with the simulation spending an inordinate amount of KMC steps / cputime simulating fast frivolous processes (FFPs) without progressing the system (reaction network). In order to achieve simulation times that are experimentally relevant or desired for predictions, a dynamic throttling algorithm involving separation of the processes into speedranks based on event frequencies has been designed and implemented with the intent of decreasing the probability of FFP events, and increasing the probability of slow process events  allowing rate limiting events to become more likely to be observed in KMC simulations. This Staggered QuasiEquilibrium Rankbased Throttling for Steadystate (SQERTSS) algorithm designed for use in achieving and simulating steadystate conditions in KMC simulations. Lastly, as shown in this work, the SQERTSS algorithm also works for transient conditions: themore »
 Authors:
 Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Materials Science and Engineering
 Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Chemical Sciences Division
 Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Materials Science and Engineering; Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Mechanical Engineering
 Publication Date:
 Research Org.:
 Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
 Sponsoring Org.:
 USDOE Laboratory Directed Research and Development (LDRD) Program
 OSTI Identifier:
 1376303
 Grant/Contract Number:
 AC0500OR22725
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Computer Physics Communications
 Additional Journal Information:
 Journal Volume: 219; Journal Issue: C; Journal ID: ISSN 00104655
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 97 MATHEMATICS AND COMPUTING; Kinetic Monte Carlo; Stiffness; Steady state; Nonequilibrium
Citation Formats
Danielson, Thomas, Sutton, Jonathan E., Hin, Céline, and Savara, Aditya. SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations. United States: N. p., 2017.
Web. doi:10.1016/j.cpc.2017.05.016.
Danielson, Thomas, Sutton, Jonathan E., Hin, Céline, & Savara, Aditya. SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations. United States. doi:10.1016/j.cpc.2017.05.016.
Danielson, Thomas, Sutton, Jonathan E., Hin, Céline, and Savara, Aditya. 2017.
"SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations". United States.
doi:10.1016/j.cpc.2017.05.016.
@article{osti_1376303,
title = {SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations},
author = {Danielson, Thomas and Sutton, Jonathan E. and Hin, Céline and Savara, Aditya},
abstractNote = {Lattice based Kinetic Monte Carlo (KMC) simulations offer a powerful simulation technique for investigating large reaction networks while retaining spatial configuration information, unlike ordinary differential equations. However, large chemical reaction networks can contain reaction processes with rates spanning multiple orders of magnitude. This can lead to the problem of “KMC stiffness” (similar to stiffness in differential equations), where the computational expense has the potential to be overwhelmed by very short timesteps during KMC simulations, with the simulation spending an inordinate amount of KMC steps / cputime simulating fast frivolous processes (FFPs) without progressing the system (reaction network). In order to achieve simulation times that are experimentally relevant or desired for predictions, a dynamic throttling algorithm involving separation of the processes into speedranks based on event frequencies has been designed and implemented with the intent of decreasing the probability of FFP events, and increasing the probability of slow process events  allowing rate limiting events to become more likely to be observed in KMC simulations. This Staggered QuasiEquilibrium Rankbased Throttling for Steadystate (SQERTSS) algorithm designed for use in achieving and simulating steadystate conditions in KMC simulations. Lastly, as shown in this work, the SQERTSS algorithm also works for transient conditions: the correct configuration space and final state will still be achieved if the required assumptions are not violated, with the caveat that the sizes of the timesteps may be distorted during the transient period.},
doi = {10.1016/j.cpc.2017.05.016},
journal = {Computer Physics Communications},
number = C,
volume = 219,
place = {United States},
year = 2017,
month = 6
}
Web of Science

Progress has been rapid in increasing the efficiency of energy conversion in nanoparticles. However, extraction of the photogenerated charge carriers remains challenging. Encouragingly, the charge mobility has been improved recently by driving nanoparticle (NP) films across the metalinsulator transition (MIT). To simulate MIT in NP films, we developed a hierarchical Kinetic Monte Carlo transport model. Electrons transfer between neighboring NPs via activated hopping when the NP energies differ by more than an overlap energy, but transfer by a nonactivated quantum delocalization, if the NP energies are closer than the overlap energy. As the overlap energy increases, emerging percolating clusters supportmore »

Hydrogenation of AcetyleneEthylene Mixtures over Pd and PdAg Alloys: FirstPrinciples Based Kinetic Monte Carlo Simulations
The kinetics for the selective hydrogenation of acetyleneethylene mixtures over model Pd(111) and bimetallic PdAg alloy surfaces were examined using first principles based kinetic Monte Carlo (KMC) simulations to elucidate the effects of alloying as well as process conditions (temperature and hydrogen partial pressure). The mechanisms that control the selective and unselective routes which included hydrogenation, dehydrogenation and CC bond breaking pathways were analyzed using firstprinciple density functional theory (DFT) calculations. The results were used to construct an intrinsic kinetic database that was used in a variable time step kinetic Monte Carlo simulation to follow the kinetics and the molecularmore » 
Evaluation of the interindividual human variation in bioactivation of methyleugenol using physiologically based kinetic modeling and Monte Carlo simulations
The present study aims at predicting the level of formation of the ultimate carcinogenic metabolite of methyleugenol, 1′sulfooxymethyleugenol, in the human population by taking variability in key bioactivation and detoxification reactions into account using Monte Carlo simulations. Depending on the metabolic route, variation was simulated based on kinetic constants obtained from incubations with a range of individual human liver fractions or by combining kinetic constants obtained for specific isoenzymes with literature reported human variation in the activity of these enzymes. The results of the study indicate that formation of 1′sulfooxymethyleugenol is predominantly affected by variation in i) P450 1A2catalyzed bioactivationmore » 
Monte Carlo Simulations of HighSpeed, TimeGated MCPbased Xray Detectors: Saturation Effects in DC Pulsed Modes and Detector Dynamic Range
We present here results of continued efforts to understand the performance of microchannel plate (MCP)–based, highspeed, gated, xray detectors. This work involves the continued improvement of a Monte Carlo simulation code to describe MCP performance coupled with experimental efforts to better characterize such detectors. Our goal is a quantitative description of MCP saturation behavior in both static and pulsed modes. A new model of charge buildup on the walls of the MCP channels is briefly described. The simulation results agree favorably with experimental data obtained with a shortpulse, highintensity ultraviolet (UV) laser. These results indicate that a weak saturation canmore » 
Monte Carlo simulations of highspeed, timegated microchannelplatebased xray detectors: Saturation effects in dc and pulsed modes and detector dynamic range
We present here results of continued efforts to understand the performance of microchannel plate (MCP)based, highspeed, gated, xray detectors. This work involves the continued improvement of a Monte Carlo simulation code to describe MCP performance coupled with experimental efforts to better characterize such detectors. Our goal is a quantitative description of MCP saturation behavior in both static and pulsed modes. A new model of charge buildup on the walls of the MCP channels is briefly described. The simulation results are compared to experimental data obtained with a shortpulse, highintensity ultraviolet laser, and good agreement is found. These results indicate thatmore »