Synchronous parallel kinetic Monte Carlo for continuum diffusion-reaction systems

  • Authors:

  • E. Martínez;J. Marian;M. H. Kalos;J. M. Perlado

  • Affiliations:

  • Lawrence Livermore National Laboratory, Livermore, CA 94551, USA and Instituto de Fusión Nuclear, Universidad Politécnica de Madrid, 28006 Madrid, Spain;Lawrence Livermore National Laboratory, Livermore, CA 94551, USA;Lawrence Livermore National Laboratory, Livermore, CA 94551, USA;Instituto de Fusión Nuclear, Universidad Politécnica de Madrid, 28006 Madrid, Spain

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2008

Quantified Score

Hi-index 31.46

Visualization

Abstract

A novel parallel kinetic Monte Carlo (kMC) algorithm formulated on the basis of perfect time synchronicity is presented. The algorithm is intended as a generalization of the standard n-fold kMC method, and is trivially implemented in parallel architectures. In its present form, the algorithm is not rigorous in the sense that boundary conflicts are ignored. We demonstrate, however, that, in their absence, or if they were correctly accounted for, our algorithm solves the same master equation as the serial method. We test the validity and parallel performance of the method by solving several pure diffusion problems (i.e. with no particle interactions) with known analytical solution. We also study diffusion-reaction systems with known asymptotic behavior and find that, for large systems with interaction radii smaller than the typical diffusion length, boundary conflicts are negligible and do not affect the global kinetic evolution, which is seen to agree with the expected analytical behavior. Our method is a controlled approximation in the sense that the error incurred by ignoring boundary conflicts can be quantified intrinsically, during the course of a simulation, and decreased arbitrarily (controlled) by modifying a few problem-dependent simulation parameters.