First-Passage Kinetic Monte Carlo method

Abstract

We present a new efficient method for Monte Carlo simulations of diffusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditional small diffusion hops and propagates the diffusing particles over long distances through a sequence of super-hops, one particle at a time. By partitioning the simulation space into non-overlapping protecting domains each containing only one or two particles, the algorithm factorizes the N-body problem of collisions among multiple Brownian particles into a set of much simpler single-body and two-body problems. Efficient propagation of particles inside their protective domains is enabled through the use of time-dependent Green's functions (propagators) obtained as solutions for the first-passage statistics of random walks. The resulting Monte Carlo algorithm is event-driven and asynchronous; each Brownian particle propagates inside its own protective domain and on its own time clock. The algorithm reproduces the statistics of the underlying Monte-Carlo model exactly. Extensive numerical examples demonstrate that for an important class of diffusion-reaction models the new algorithm is efficient at low particle densities, where other existing algorithms slow down severely.

Elaborates on Phys. Rev. Lett., 97:230602, 2006. Small revisions from prior version only

Cited by in corpus (26)

• A First-Passage Kinetic Monte Carlo Algorithm for Complex Diffusion-Reaction Systems
• From Molecular Dynamics to Brownian Dynamics
• Efficient Reactive Brownian Dynamics
• Event Driven Langevin simulations of Hard Spheres
• Particle-based Multiscale Modeling of Calcium Puff Dynamics
• First-Passage Exponents of Multiple Random Walks
• A First-Passage Kinetic Monte Carlo Method for Reaction-Drift-Diffusion Processes
• Coupling all-atom molecular dynamics simulations of ions in water with Brownian dynamics
• An efficient multi-scale Green's Functions Reaction Dynamics scheme
• First-passage time to clear the way for receptor-ligand binding in a crowded environment
• Efficient kinetic Monte Carlo method for reaction-diffusion processes with spatially varying annihilation rates
• Diffusive Exit Rates Through Pores in Membrane-Enclosed Structures
• Slow Kinetics of Brownian Maxima
• Numerical analysis of homogeneous and inhomogeneous intermittent search strategies
• Narrow escape problem in two-shell spherical domains
• Path integral approach to theories of diffusion-influenced reactions
• Absorption kinetics of vacancies by cavities in Aluminum: numerical characterization of sink strengths and first-passage statistics through Krylov subspace projection and eigenvalue deflation
• Kinetic Monte Carlo methods for three-dimensional diffusive capture problems in exterior domains
• Ordering statistics of 4 random walkers on the line
• Spatial-Stochastic Simulation of Reaction-Diffusion Systems
• Extinction and Survival in Two-Species Annihilation
• On the formalization of Asynchronous First Passage Algorithms
• Diffusive search and trajectories on spatial networks: a propagator approach
• Stochastic single-particle based simulations of cellular signaling embedded into computational models of cellular morphology
• Varying the resolution of the Rouse model on temporal and spatial scales: application to multiscale modelling of DNA dynamics
• Smoluchowski reaction kinetics for reactions of any order