Acceleration scheme for particle transport in kinetic Monte Carlo methods.

Kinetic Monte Carlo (kMC) simulations are frequently used to study (electro-)chemical processes within science and engineering. kMC methods provide insight into the interplay of stochastic processes and can link atomistic material properties with macroscopic characteristics. Significant problems concerning the computational demand arise if processes with large time disparities are competing. Acceleration algorithms are required to make slow processes accessible. Especially, the accelerated superbasin kMC (AS-kMC) scheme has been frequently applied within chemical reaction networks. For larger systems, the computational overhead of the AS-kMC is significant as the computation of the superbasins is done during runtime and comes with the need for large databases. Here, we propose a novel acceleration scheme for diffusion and transport processes within kMC simulations. Critical superbasins are detected during the system initialization. Scaling factors for the critical rates within the superbasins, as well as a lower bound for the number of sightings, are derived. Our algorithm exceeds the AS-kMC in the required simulation time, which we demonstrate with a 1D-chain example. In addition, we apply the acceleration scheme to study the time-of-flight (TOF) of charge carriers within organic semiconductors. In this material class, time disparities arise due to a significant spread of transition rates. The acceleration scheme allows a significant acceleration up to a factor of 65 while keeping the error of the TOF values negligible. The computational overhead is negligible, as all superbasins only need to be computed once.

[1]  D. Gillespie A rigorous derivation of the chemical master equation , 1992 .

[2]  H. Resat,et al.  Probability-Weighted Dynamic Monte Carlo Method for Reaction Kinetics Simulations , 2001 .

[3]  Johannes Popp,et al.  Generalized Kinetic Monte Carlo Framework for Organic Electronics , 2018, Algorithms.

[4]  D. Vlachos A Review of Multiscale Analysis: Examples from Systems Biology, Materials Engineering, and Other Fluid–Surface Interacting Systems , 2005 .

[5]  A. Jansen,et al.  Dynamic Monte-Carlo simulations of reactions in heterogeneous catalysis , 1999 .

[6]  Markos A. Katsoulakis,et al.  Coarse-grained stochastic processes and kinetic Monte Carlo simulators for the diffusion of interacting particles , 2003 .

[7]  Jianguo Dai,et al.  Coarse-grained lattice kinetic Monte Carlo simulation of systems of strongly interacting particles. , 2008, The Journal of chemical physics.

[8]  L. Zepeda-Ruiz,et al.  Monte Carlo Simulations of Crystal Growth , 2015 .

[9]  K. Kaski,et al.  Kinetic Monte Carlo simulation of nucleation on patterned substrates , 2000 .

[10]  R. Drautz,et al.  Kinetic Monte Carlo simulations of vacancy diffusion in nondilute Ni- X ( X = Re, W, Ta) alloys , 2018, Physical Review Materials.

[11]  B. Engels,et al.  First-principles calculations of anisotropic charge-carrier mobilities in organic semiconductor crystals , 2011, 1102.4289.

[12]  W M Young,et al.  Monte Carlo studies of vacancy migration in binary ordered alloys: I , 1966 .

[13]  K. Reuter,et al.  Charge Transport in Molecular Materials: An Assessment of Computational Methods. , 2017, Chemical reviews.

[14]  Tiejun Li,et al.  Analysis of Explicit Tau-Leaping Schemes for Simulating Chemically Reacting Systems , 2007, Multiscale Model. Simul..

[15]  A. Gagliardi,et al.  Charge carrier mobility of disordered organic semiconductors with correlated energetic and spatial disorder. , 2018, Physical chemistry chemical physics : PCCP.

[16]  Abhijit Chatterjee,et al.  Spatially adaptive lattice coarse-grained Monte Carlo simulations for diffusion of interacting molecules. , 2004, The Journal of chemical physics.

[17]  Michail Stamatakis,et al.  Unraveling the Complexity of Catalytic Reactions via Kinetic Monte Carlo Simulation: Current Status and Frontiers , 2012 .

[18]  Wonryull Koh,et al.  An accelerated algorithm for discrete stochastic simulation of reaction-diffusion systems using gradient-based diffusion and tau-leaping. , 2011, The Journal of chemical physics.

[19]  Abhijit Chatterjee,et al.  An overview of spatial microscopic and accelerated kinetic Monte Carlo methods , 2007 .

[20]  M. Neurock,et al.  Generalized Temporal Acceleration Scheme for Kinetic Monte Carlo Simulations of Surface Catalytic Processes by Scaling the Rates of Fast Reactions. , 2017, Journal of chemical theory and computation.

[21]  R. Smith,et al.  Dynamic simulation of crystal growth by Monte Carlo method—I. Model description and kinetics , 1992 .

[22]  A. Chatterjee,et al.  Net-event kinetic Monte Carlo for overcoming stiffness in spatially homogeneous and distributed systems , 2005, Comput. Chem. Eng..

[23]  J. B. Adams,et al.  Kinetic lattice Monte Carlo model for oxygen vacancy diffusion in praseodymium doped ceria: Applications to materials design , 2011 .

[24]  H. Bässler Charge Transport in Disordered Organic Photoconductors a Monte Carlo Simulation Study , 1993 .

[25]  Muruhan Rathinam,et al.  Consistency and Stability of Tau-Leaping Schemes for Chemical Reaction Systems , 2005, Multiscale Model. Simul..

[26]  E. Abrahams,et al.  Impurity Conduction at Low Concentrations , 1960 .

[27]  S. Piana,et al.  Simulating micrometre-scale crystal growth from solution , 2005, Nature.

[28]  Peifeng Zhang,et al.  Kinetic Monte Carlo simulation of Cu thin film growth , 2004 .

[29]  Abhijit Chatterjee,et al.  Accurate acceleration of kinetic Monte Carlo simulations through the modification of rate constants. , 2010, The Journal of chemical physics.

[30]  D G Vlachos Temporal coarse-graining of microscopic-lattice kinetic Monte Carlo simulations via tau leaping. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[31]  Desmond J. Higham,et al.  Modeling and Simulating Chemical Reactions , 2008, SIAM Rev..

[32]  Abhijit Chatterjee,et al.  Multiscale spatial Monte Carlo simulations: multigriding, computational singular perturbation, and hierarchical stochastic closures. , 2006, The Journal of chemical physics.

[33]  A. Gagliardi,et al.  Machine Learning–Based Charge Transport Computation for Pentacene , 2018, Advanced Theory and Simulations.

[34]  D. Gillespie Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .