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

Title: A derivation and scalable implementation of the synchronous parallel kinetic Monte Carlo method for simulating long-time dynamics

Abstract

Kinetic Monte Carlo (KMC) simulations are used to study long-time dynamics of a wide variety of systems. Unfortunately, the conventional KMC algorithm is not scalable to larger systems, since its time scale is inversely proportional to the simulated system size. A promising approach to resolving this issue is the synchronous parallel KMC (SPKMC) algorithm, which makes the time scale size-independent. This paper introduces a formal derivation of the SPKMC algorithm based on local transition-state and time-dependent Hartree approximations, as well as its scalable parallel implementation based on a dual linked-list cell method. The resulting algorithm has achieved a weak-scaling parallel efficiency of 0.935 on 1024 Intel Xeon processors for simulating biological electron transfer dynamics in a 4.2 billion-heme system, as well as decent strong-scaling parallel efficiency. The parallel code has been used to simulate a lattice of cytochrome complexes on a bacterial-membrane nanowire, and it is broadly applicable to other problems such as computational synthesis of new materials.

Authors:
 [1];  [1];  [1];  [1];  [1]
  1. University of Southern California, Los Angeles, CA (United States)
Publication Date:
Research Org.:
Univ. of Southern California, Los Angeles, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES). Chemical Sciences, Geosciences & Biosciences Division (CSGB); Defense Threat Reduction Agency (DTRA); USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1538203
Alternate Identifier(s):
OSTI ID: 1550232
Grant/Contract Number:  
SC0010609; SC0014607; HDTRA1-14-1-0074; FG02-13ER16415
Resource Type:
Accepted Manuscript
Journal Name:
Computer Physics Communications
Additional Journal Information:
Journal Volume: 219; Journal Issue: C; Journal ID: ISSN 0010-4655
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; kinetic Monte Carlo simulation; divide-and-conquer algorithm; parallel computing

Citation Formats

Byun, Hye Suk, El-Naggar, Mohamed Y., Kalia, Rajiv K., Nakano, Aiichiro, and Vashishta, Priya. A derivation and scalable implementation of the synchronous parallel kinetic Monte Carlo method for simulating long-time dynamics. United States: N. p., 2017. Web. doi:10.1016/j.cpc.2017.05.028.
Byun, Hye Suk, El-Naggar, Mohamed Y., Kalia, Rajiv K., Nakano, Aiichiro, & Vashishta, Priya. A derivation and scalable implementation of the synchronous parallel kinetic Monte Carlo method for simulating long-time dynamics. United States. https://doi.org/10.1016/j.cpc.2017.05.028
Byun, Hye Suk, El-Naggar, Mohamed Y., Kalia, Rajiv K., Nakano, Aiichiro, and Vashishta, Priya. Tue . "A derivation and scalable implementation of the synchronous parallel kinetic Monte Carlo method for simulating long-time dynamics". United States. https://doi.org/10.1016/j.cpc.2017.05.028. https://www.osti.gov/servlets/purl/1538203.
@article{osti_1538203,
title = {A derivation and scalable implementation of the synchronous parallel kinetic Monte Carlo method for simulating long-time dynamics},
author = {Byun, Hye Suk and El-Naggar, Mohamed Y. and Kalia, Rajiv K. and Nakano, Aiichiro and Vashishta, Priya},
abstractNote = {Kinetic Monte Carlo (KMC) simulations are used to study long-time dynamics of a wide variety of systems. Unfortunately, the conventional KMC algorithm is not scalable to larger systems, since its time scale is inversely proportional to the simulated system size. A promising approach to resolving this issue is the synchronous parallel KMC (SPKMC) algorithm, which makes the time scale size-independent. This paper introduces a formal derivation of the SPKMC algorithm based on local transition-state and time-dependent Hartree approximations, as well as its scalable parallel implementation based on a dual linked-list cell method. The resulting algorithm has achieved a weak-scaling parallel efficiency of 0.935 on 1024 Intel Xeon processors for simulating biological electron transfer dynamics in a 4.2 billion-heme system, as well as decent strong-scaling parallel efficiency. The parallel code has been used to simulate a lattice of cytochrome complexes on a bacterial-membrane nanowire, and it is broadly applicable to other problems such as computational synthesis of new materials.},
doi = {10.1016/j.cpc.2017.05.028},
journal = {Computer Physics Communications},
number = C,
volume = 219,
place = {United States},
year = {Tue Jun 27 00:00:00 EDT 2017},
month = {Tue Jun 27 00:00:00 EDT 2017}
}

Journal Article:

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

Save / Share:

Works referenced in this record:

Rigorous synchronous relaxation algorithm for parallel kinetic Monte Carlo simulations of thin film growth
journal, March 2005


Long-range electron transport in Geobacter sulfurreducens biofilms is redox gradient-driven
journal, September 2012

  • Snider, R. M.; Strycharz-Glaven, S. M.; Tsoi, S. D.
  • Proceedings of the National Academy of Sciences, Vol. 109, Issue 38
  • DOI: 10.1073/pnas.1209829109

Quantum Molecular Dynamics in the Post-Petaflops Era
journal, November 2015

  • Romero, Nichols A.; Nakano, Aiichiro; Riley, Katherine M.
  • Computer, Vol. 48, Issue 11
  • DOI: 10.1109/MC.2015.337

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


Kinetic Monte Carlo simulations with minimal searching
journal, February 2002


Exascale computing and big data
journal, June 2015

  • Reed, Daniel A.; Dongarra, Jack
  • Communications of the ACM, Vol. 58, Issue 7
  • DOI: 10.1145/2699414

Thermally activated long range electron transport in living biofilms
journal, January 2015

  • Yates, Matthew D.; Golden, Joel P.; Roy, Jared
  • Physical Chemistry Chemical Physics, Vol. 17, Issue 48
  • DOI: 10.1039/C5CP05152E

Dynamic correlations in electron liquids. I. General formalism
journal, March 1989


Billion-atom synchronous parallel kinetic Monte Carlo simulations of critical 3D Ising systems
journal, February 2011

  • Martínez, E.; Monasterio, P. R.; Marian, J.
  • Journal of Computational Physics, Vol. 230, Issue 4
  • DOI: 10.1016/j.jcp.2010.11.006

Rapid electron exchange between surface-exposed bacterial cytochromes and Fe(III) minerals
journal, March 2013

  • White, G. F.; Shi, Z.; Shi, L.
  • Proceedings of the National Academy of Sciences, Vol. 110, Issue 16
  • DOI: 10.1073/pnas.1220074110

Enhanced sampling in molecular dynamics: use of the time-dependent Hartree approximation for a simulation of carbon monoxide diffusion through myoglobin
journal, December 1990

  • Elber, R.; Karplus, M.
  • Journal of the American Chemical Society, Vol. 112, Issue 25
  • DOI: 10.1021/ja00181a020

Dielectric Genome of van der Waals Heterostructures
journal, June 2015


Synchronous parallel kinetic Monte Carlo for continuum diffusion-reaction systems
journal, April 2008


Kinetic Monte Carlo Simulations and Molecular Conductance Measurements of the Bacterial Decaheme Cytochrome MtrF
journal, September 2014

  • Byun, Hye Suk; Pirbadian, Sahand; Nakano, Aiichiro
  • ChemElectroChem, Vol. 1, Issue 11
  • DOI: 10.1002/celc.201402211

A framework for stochastic simulations and visualization of biological electron-transfer dynamics
journal, August 2015


Suppressing Roughness of Virtual Times in Parallel Discrete-Event Simulations
journal, January 2003


Shewanella oneidensis MR-1 nanowires are outer membrane and periplasmic extensions of the extracellular electron transport components
journal, August 2014

  • Pirbadian, S.; Barchinger, S. E.; Leung, K. M.
  • Proceedings of the National Academy of Sciences, Vol. 111, Issue 35
  • DOI: 10.1073/pnas.1410551111

Identification of Relaxation and Diffusion Mechanisms in Amorphous Silicon
journal, August 1998


Multiresolution molecular dynamics algorithm for realistic materials modeling on parallel computers
journal, December 1994


\mathcal{O}(N) methods in electronic structure calculations
journal, February 2012


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

Multiresolution load balancing in curved space: the wavelet representation
journal, June 1999


The ‘porin-cytochrome’ model for microbe-to-mineral electron transfer: Microbe-to-mineral electron transfer
journal, May 2012


A combined electrochemical and optical trapping platform for measuring single cell respiration rates at electrode interfaces
journal, June 2015

  • Gross, Benjamin J.; El-Naggar, Mohamed Y.
  • Review of Scientific Instruments, Vol. 86, Issue 6
  • DOI: 10.1063/1.4922853

Theoretical foundations of dynamical Monte Carlo simulations
journal, July 1991

  • Fichthorn, Kristen A.; Weinberg, W. H.
  • The Journal of Chemical Physics, Vol. 95, Issue 2
  • DOI: 10.1063/1.461138

Electrical transport along bacterial nanowires from Shewanella oneidensis MR-1
journal, October 2010

  • El-Naggar, M. Y.; Wanger, G.; Leung, K. M.
  • Proceedings of the National Academy of Sciences, Vol. 107, Issue 42
  • DOI: 10.1073/pnas.1004880107

Current Status of Transition-State Theory
journal, January 1996

  • Truhlar, Donald G.; Garrett, Bruce C.; Klippenstein, Stephen J.
  • The Journal of Physical Chemistry, Vol. 100, Issue 31
  • DOI: 10.1021/jp953748q

Pathfinder: A parallel search algorithm for concerted atomistic events
journal, February 2007


Nanoscopic mechanisms of singlet fission in amorphous molecular solid
journal, April 2013

  • Mou, Weiwei; Hattori, Shinnosuke; Rajak, Pankaj
  • Applied Physics Letters, Vol. 102, Issue 17
  • DOI: 10.1063/1.4795138

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions
journal, December 1976


Exponential multiplicity of inherent structures
journal, January 1999


A fast algorithm for particle simulations
journal, December 1987


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


A space–time-ensemble parallel nudged elastic band algorithm for molecular kinetics simulation
journal, February 2008


A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations
journal, May 2014

  • Shimojo, Fuyuki; Hattori, Shinnosuke; Kalia, Rajiv K.
  • The Journal of Chemical Physics, Vol. 140, Issue 18
  • DOI: 10.1063/1.4869342

Faster Monte Carlo simulations
journal, February 1995


Cloud-based simulations on Google Exacycle reveal ligand modulation of GPCR activation pathways
journal, December 2013

  • Kohlhoff, Kai J.; Shukla, Diwakar; Lawrenz, Morgan
  • Nature Chemistry, Vol. 6, Issue 1
  • DOI: 10.1038/nchem.1821