DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations

Journal Article · · Computer Physics Communications
 [1];  [2]; ORCiD logo [3];  [4]
  1. Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Mechanical Engineering
  2. Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Materials Science and Engineering
  3. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Chemical Sciences Division
  4. Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Mechanical Engineering; Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Materials Science and Engineering

A window-based steady-state detection algorithm has been developed for application to kinetic Monte Carlo simulation data. The algorithm, termed F-t-Pj-RG sequentially applies an F-test, a t-test, and a projection test on adjacent windows of the data while rolling (or shifting) and growing the windows when any of the tests fail. In aggregate, the algorithm is able to (a) automatically reject the warm-up period as not being at steady-state, as well as (b) determine an appropriate window size for converged statistics when sampling the data, which is necessary for detection of steady-state, and (c) detect steady-state within a particular tolerance. The last step, the projection test, is actually an oscillating-slope projection test, and is performed on j sequential data windows (i.e., more than two adjacent windows). It requires more than simply being within the user defined tolerance: the oscillating-slope projection test includes a condition that the slope must oscillate around zero when 2, which is an additional indication of steady-state. When all three tests are passed, the F-t-Pj test is passed, indicating that the prerequisites of steady-state detection have been met and also that conditions consistent with the definition of steady-state have been realized. This algorithm is applied to a variety of data sets that correspond to the diverse type of data trends that can be produced by kinetic Monte Carlo simulations. The algorithm is shown to be robust in its ability to handle differing functional forms, and is able to detect steady-state with low computational cost. Finally, the low computational cost of this method and its robustness towards varied data trends make it suitable for on-the-fly use in kinetic Monte Carlo simulations.

Research Organization:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Sponsoring Organization:
USDOE
Grant/Contract Number:
AC05-00OR22725; LOIS 8457
OSTI ID:
1495967
Alternate ID(s):
OSTI ID: 1548150
Journal Information:
Computer Physics Communications, Vol. 232, Issue C; ISSN 0010-4655
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 4 works
Citation information provided by
Web of Science

References (26)

Experiment-Based Kinetic Monte Carlo Simulations: CO Oxidation over RuO 2 (110) journal December 2011
First-Principles-Based Monte Carlo Simulation of Ethylene Hydrogenation Kinetics on Pd journal December 2000
Adlayer inhomogeneity without lateral interactions: Rationalizing correlation effects in CO oxidation at RuO 2 (110) with first-principles kinetic Monte Carlo journal February 2011
First-Principles-Based Kinetic Monte Carlo Simulation of Nitric Oxide Reduction over Platinum Nanoparticles under Lean-Burn Conditions journal November 2010
Unraveling the Complexity of Catalytic Reactions via Kinetic Monte Carlo Simulation: Current Status and Frontiers journal November 2012
Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO 2 journal April 2018
An efficient method for on-line identification of steady state journal December 1995
A steady-state detection (SSD) algorithm to detect non-stationary drifts in processes journal March 2013
Methodology of on-line optimisation applied to a chemical reactor journal February 2008
Acceleration and sensitivity analysis of lattice kinetic Monte Carlo simulations using parallel processing and rate constant rescaling journal October 2017
kmos: A lattice kinetic Monte Carlo framework journal July 2014
Theoretical foundations of dynamical Monte Carlo simulations journal July 1991
A general method for numerically simulating the stochastic time evolution of coupled chemical reactions journal December 1976
Exact stochastic simulation of coupled chemical reactions journal December 1977
The steady state and equilibrium approximations: A geometrical picture journal April 1988
SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations journal October 2017
Below-Room-Temperature C–H Bond Breaking on an Inexpensive Metal Oxide: Methanol to Formaldehyde on CeO 2 (111) journal November 2017
Equivalence of on-lattice stochastic chemical kinetics with the well-mixed chemical master equation in the limit of fast diffusion journal December 2011
Fluctuation-dominated kinetics in diffusion-controlled reactions journal July 1985
Role of density fluctuations in bimolecular reaction kinetics journal February 1978
A + B → 0 Reaction with Different Initial Patterns journal January 1996
Fractal Reaction Kinetics journal September 1988
Ion clustering in aqueous solutions probed with vibrational energy transfer journal March 2011
Vibrational Energy Transfer: An Angstrom Molecular Ruler in Studies of Ion Pairing and Clustering in Aqueous Solutions journal February 2015
Ion Segregation in Aqueous Solutions journal November 2012
A local superbasin kinetic Monte Carlo method journal April 2013

Cited By (2)

A Practical Guide to Surface Kinetic Monte Carlo Simulations journal April 2019
A Practical Guide to Surface Kinetic Monte Carlo Simulations text January 2019

Figures / Tables (20)