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

Title: Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification

Abstract

Here, a massively parallel method to build large transition rate matrices from temperature-accelerated molecular dynamics trajectories is presented. Bayesian Markov model analysis is used to estimate the expected residence time in the known state space, providing crucial uncertainty quantification for higher-scale simulation schemes such as kinetic Monte Carlo or cluster dynamics. The estimators are additionally used to optimize where exploration is performed and the degree of temperature acceleration on the fly, giving an autonomous, optimal procedure to explore the state space of complex systems. The method is tested against exactly solvable models and used to explore the dynamics of C15 interstitial defects in iron. Our uncertainty quantification scheme allows for accurate modeling of the evolution of these defects over timescales of several seconds.

Authors:
ORCiD logo [1]; ORCiD logo [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1458961
Report Number(s):
LA-UR-18-21562
Journal ID: ISSN 2475-9953; PRMHAR; TRN: US1901538
Grant/Contract Number:  
AC52-06NA25396
Resource Type:
Accepted Manuscript
Journal Name:
Physical Review Materials
Additional Journal Information:
Journal Volume: 2; Journal Issue: 5; Journal ID: ISSN 2475-9953
Publisher:
American Physical Society (APS)
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE; Material Science

Citation Formats

Swinburne, Thomas David, and Perez, Danny. Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification. United States: N. p., 2018. Web. doi:10.1103/PhysRevMaterials.2.053802.
Swinburne, Thomas David, & Perez, Danny. Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification. United States. https://doi.org/10.1103/PhysRevMaterials.2.053802
Swinburne, Thomas David, and Perez, Danny. Mon . "Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification". United States. https://doi.org/10.1103/PhysRevMaterials.2.053802. https://www.osti.gov/servlets/purl/1458961.
@article{osti_1458961,
title = {Self-optimized construction of transition rate matrices from accelerated atomistic simulations with Bayesian uncertainty quantification},
author = {Swinburne, Thomas David and Perez, Danny},
abstractNote = {Here, a massively parallel method to build large transition rate matrices from temperature-accelerated molecular dynamics trajectories is presented. Bayesian Markov model analysis is used to estimate the expected residence time in the known state space, providing crucial uncertainty quantification for higher-scale simulation schemes such as kinetic Monte Carlo or cluster dynamics. The estimators are additionally used to optimize where exploration is performed and the degree of temperature acceleration on the fly, giving an autonomous, optimal procedure to explore the state space of complex systems. The method is tested against exactly solvable models and used to explore the dynamics of C15 interstitial defects in iron. Our uncertainty quantification scheme allows for accurate modeling of the evolution of these defects over timescales of several seconds.},
doi = {10.1103/PhysRevMaterials.2.053802},
journal = {Physical Review Materials},
number = 5,
volume = 2,
place = {United States},
year = {Mon May 21 00:00:00 EDT 2018},
month = {Mon May 21 00:00:00 EDT 2018}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 30 works
Citation information provided by
Web of Science

Figures / Tables:

FIG. 1 FIG. 1: TAMMBER workflow. TAD MD produces interstate transition trajectories which are analyzed by Bayesian rate estimators and static calculation. An absorbing Markov chain gives then gives the expected residence time and optimally allocates resources and the degree of temperature acceleration. The cycle is then repeated until the target residencemore » time is achieved.« less

Save / Share:

Works referenced in this record:

Atomistic Simulations of Activated Processes in Materials
journal, July 2017


Discrete path sampling
journal, October 2002


Temperature-accelerated dynamics for simulation of infrequent events
journal, June 2000

  • So/rensen, Mads R.; Voter, Arthur F.
  • The Journal of Chemical Physics, Vol. 112, Issue 21
  • DOI: 10.1063/1.481576

Unsupervised Calculation of Free Energy Barriers in Large Crystalline Systems
journal, March 2018


The Modern Temperature-Accelerated Dynamics Approach
journal, June 2016


A new algorithm for Monte Carlo simulation of Ising spin systems
journal, January 1975


Kinetic activation-relaxation technique
journal, October 2011


Dynamical integration of a Markovian web: A first passage time approach
journal, August 2007

  • Boulougouris, Georgios C.; Theodorou, Doros N.
  • The Journal of Chemical Physics, Vol. 127, Issue 8
  • DOI: 10.1063/1.2753153

Monte Carlo Sampling of a Markov Web
journal, March 2005

  • Boulougouris, Georgios C.; Frenkel, Daan
  • Journal of Chemical Theory and Computation, Vol. 1, Issue 3
  • DOI: 10.1021/ct049900m

The nature of high-energy radiation damage in iron
journal, February 2013


Adaptive temperature-accelerated dynamics
journal, February 2011

  • Shim, Yunsic; Amar, Jacques G.
  • The Journal of Chemical Physics, Vol. 134, Issue 5
  • DOI: 10.1063/1.3541823

Structural Relaxation Made Simple
journal, October 2006


Stability of self-interstitial clusters with C15 Laves phase structure in iron
journal, June 2014


Practical graph isomorphism, II
journal, January 2014


Hyperdynamics: Accelerated Molecular Dynamics of Infrequent Events
journal, May 1997


Formation of prismatic loops from C15 Laves phase interstitial clusters in body-centered cubic iron
journal, March 2015


Fast Parallel Algorithms for Short-Range Molecular Dynamics
journal, March 1995


Commentary: The Materials Project: A materials genome approach to accelerating materials innovation
journal, July 2013

  • Jain, Anubhav; Ong, Shyue Ping; Hautier, Geoffroy
  • APL Materials, Vol. 1, Issue 1
  • DOI: 10.1063/1.4812323

Long-Time Dynamics through Parallel Trajectory Splicing
journal, December 2015

  • Perez, Danny; Cubuk, Ekin D.; Waterland, Amos
  • Journal of Chemical Theory and Computation, Vol. 12, Issue 1
  • DOI: 10.1021/acs.jctc.5b00916

Tutorial on maximum likelihood estimation
journal, February 2003


Reaction-rate theory: fifty years after Kramers
journal, April 1990

  • Hänggi, Peter; Talkner, Peter; Borkovec, Michal
  • Reviews of Modern Physics, Vol. 62, Issue 2
  • DOI: 10.1103/RevModPhys.62.251

Multiscale modelling of defect kinetics in irradiated iron
journal, December 2004

  • Fu, Chu-Chun; Torre, Jacques Dalla; Willaime, François
  • Nature Materials, Vol. 4, Issue 1
  • DOI: 10.1038/nmat1286

Big Data of Materials Science: Critical Role of the Descriptor
journal, March 2015


Uncertainty in a Markov state model with missing states and rates: Application to a room temperature kinetic model obtained using high temperature molecular dynamics
journal, September 2015

  • Chatterjee, Abhijit; Bhattacharya, Swati
  • The Journal of Chemical Physics, Vol. 143, Issue 11
  • DOI: 10.1063/1.4930976

Benchmarks for Characterization of Minima, Transition States, and Pathways in Atomic, Molecular, and Condensed Matter Systems
journal, November 2014

  • Chill, Samuel T.; Stevenson, Jacob; Ruehle, Victor
  • Journal of Chemical Theory and Computation, Vol. 10, Issue 12
  • DOI: 10.1021/ct5008718

Scalable kernel polynomial method for calculating transition rates
journal, June 2013


Comparison of empirical interatomic potentials for iron applied to radiation damage studies
journal, November 2010


A climbing image nudged elastic band method for finding saddle points and minimum energy paths
journal, December 2000

  • Henkelman, Graeme; Uberuaga, Blas P.; Jónsson, Hannes
  • The Journal of Chemical Physics, Vol. 113, Issue 22, p. 9901-9904
  • DOI: 10.1063/1.1329672

Monte Carlo Algorithms with Absorbing Markov Chains: Fast Local Algorithms for Slow Dynamics
journal, January 1995


Irradiation-Induced Formation of Nanocrystallites with C 15 Laves Phase Structure in bcc Iron
journal, January 2012


Molecular dynamics saddle search adaptive kinetic Monte Carlo
journal, June 2014

  • Chill, Samuel T.; Henkelman, Graeme
  • The Journal of Chemical Physics, Vol. 140, Issue 21
  • DOI: 10.1063/1.4880721

An invariant form for the prior probability in estimation problems
journal, September 1946

  • Jeffreys, Harold
  • Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, Vol. 186, Issue 1007, p. 453-461
  • DOI: 10.1098/rspa.1946.0056

Everything you wanted to know about Markov State Models but were afraid to ask
journal, September 2010


Parallel replica method for dynamics of infrequent events
journal, June 1998


A new class of enhanced kinetic sampling methods for building Markov state models
journal, October 2017

  • Bhoutekar, Arti; Ghosh, Susmita; Bhattacharya, Swati
  • The Journal of Chemical Physics, Vol. 147, Issue 15
  • DOI: 10.1063/1.4984932

Big Data of Materials Science - Critical Role of the Descriptor
text, January 2014


Works referencing / citing this record:

Atomistic-object kinetic Monte Carlo simulations of irradiation damage in tungsten
journal, May 2019

  • Mason, D. R.; Sand, A. E.; Dudarev, S. L.
  • Modelling and Simulation in Materials Science and Engineering, Vol. 27, Issue 5
  • DOI: 10.1088/1361-651x/ab1a1e

Atomistic-Object Kinetic Monte Carlo simulations of irradiation damage in tungsten
text, January 2019


Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.