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

doi:10.1016/j.cpc.2018.05.013

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

Full Record
Other Related Research

Abstract

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 varietymore »« less

Authors:

Nellis, Chris ^[1]; Danielson, Thomas ^[2];

^[3]; Hin, Celine ^[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
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 Mechanical Engineering; Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States). Dept. of Materials Science and Engineering

Publication Date:: 2018-06-23

Research Org.:: Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)

Sponsoring Org.:: USDOE

OSTI Identifier:: 1495967

Alternate Identifier(s):: OSTI ID: 1548150

Grant/Contract Number:: AC05-00OR22725; LOIS 8457

Resource Type:: Accepted Manuscript

Journal Name:: Computer Physics Communications

Additional Journal Information:: Journal Volume: 232; Journal Issue: C; Journal ID: ISSN 0010-4655

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; Steady-state detection; Kinetic Monte Carlo simulations; Complex chemical reactions

Citation Formats


                    Nellis, Chris, Danielson, Thomas, Savara, Aditya, and Hin, Celine. The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations.  United States: N. p., 2018. 
        Web.  doi:10.1016/j.cpc.2018.05.013.

Copy to clipboard


                    Nellis, Chris, Danielson, Thomas, Savara, Aditya, & Hin, Celine. The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations.  United States.  doi:10.1016/j.cpc.2018.05.013.

Copy to clipboard


                    Nellis, Chris, Danielson, Thomas, Savara, Aditya, and Hin, Celine. Sat .  
        "The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations".  United States.  doi:10.1016/j.cpc.2018.05.013.  https://www.osti.gov/servlets/purl/1495967.

Copy to clipboard


                    
@article{osti_1495967,

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

  author       = {Nellis, Chris and Danielson, Thomas and Savara, Aditya and Hin, Celine},

  abstractNote = {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.},

  doi          = {10.1016/j.cpc.2018.05.013},

  journal      = {Computer Physics Communications},

  number       = C,

  volume       = 232,

  place        = {United States},

  year         = {2018},

  month        = {6}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

DOI: 10.1016/j.cpc.2018.05.013

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 1 work

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

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

Journal Article Nellis, Chris ; Danielson, Thomas ; Savara, Aditya ; ... - Computer Physics Communications
Cited by 1
DOI: 10.1016/j.cpc.2018.05.013
SQERT-T: alleviating kinetic Monte Carlo (KMC)-stiffness in transient KMC simulations [SQERT-T: Alleviating KMC-stiffness in Transient Kinetic Monte Carlo Simulations]

Journal Article Savara, Aditya Ashi ; Sutton, Jonathan E. - Journal of Physics. Condensed Matter

Lattice based kinetic Monte Carlo (KMC) is often used for simulating the dynamics of systems at a supramolecular scale, based on molecular scale transitions. A common challenge in KMC simulations is rapid 'back-and-forth' reactions, which dominate the events executed during simulations and inhibit the ability for simulations to reach longer time scales. Such processes are fast frivolous processes (FFPs) and are one manifestation of a phenomenon referred to as KMC-stiffness. Here, an algorithm for staggered quasi-equilibrium rank-based throttling geared towards transient kinetics (SQERT-T) is presented. Within the SQERT-T methodology, a pace-restrictor reaction and an FFP floor are utilized along withmore »« less
DOI: 10.1088/1361-648X/aacb6d

Full Text Available
SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations

Journal Article Danielson, Thomas ; Sutton, Jonathan E. ; Hin, Céline ; ... - Computer Physics Communications

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 time-steps during KMC simulations, with the simulation spending an inordinate amount of KMC steps / cpu-time simulating fast frivolous processes (FFPs) without progressing the system (reaction network). In order tomore »« less
Cited by 7
DOI: 10.1016/j.cpc.2017.05.016

Full Text Available
Kinetic Activation-Relaxation Technique and Self-Evolving Atomistic Kinetic Monte Carlo: Comparison of on-the-fly kinetic Monte Carlo algorithms

Journal Article Beland, Laurent Karim ; Osetskiy, Yury N. ; Stoller, Roger E. ; ... - Computational Materials Science

Here, we present a comparison of the Kinetic Activation–Relaxation Technique (k-ART) and the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC), two off-lattice, on-the-fly Kinetic Monte Carlo (KMC) techniques that were recently used to solve several materials science problems. We show that if the initial displacements are localized the dimer method and the Activation–Relaxation Technique nouveau provide similar performance. We also show that k-ART and SEAKMC, although based on different approximations, are in agreement with each other, as demonstrated by the examples of 50 vacancies in a 1950-atom Fe box and of interstitial loops in 16,000-atom boxes. Generally speaking, k-ART’s treatment ofmore »« less
Cited by 5
DOI: 10.1016/j.commatsci.2014.12.001

Full Text Available
Steady state likelihood ratio sensitivity analysis for stiff kinetic Monte Carlo simulations

Journal Article Núñez, M. ; Vlachos, D. G. - Journal of Chemical Physics
Cited by 6
DOI: 10.1063/1.4905957

Similar Records