Monte Carlo simulation with fixed steplength for diffusion processes in nonhomogeneous media
arXiv:1212.0362 · doi:10.1016/j.jcp.2012.12.029
Abstract
Monte Carlo simulation is one of the most important tools in the study of diffusion processes. For constant diffusion coefficients, an appropriate Gaussian distribution of particle's steplengths can generate exact results, when compared with integration of the diffusion equation. It is important to notice that the same method is completely erroneous when applied to non-homogeneous diffusion coefficients. A simple alternative, jumping at fixed steplegths with appropriate transition probabilities, produces correct results. Here, a model for diffusion of calcium ions in the neuromuscular junction of the crayfish is used as a test to compare Monte Carlo simulation with fixed and Gaussian steplegth.
References in corpus (7)
- Kinetic Equations for Diffusion in the Presence of Entropic Barriers
- Biased diffusion in confined media: Test of the Fick-Jacobs approximation and validity criteria
- Entropic particle transport: higher order corrections to the Fick-Jacobs diffusion equation
- Entropic particle transport in periodic channels
- Anderson localization as position-dependent diffusion in disordered waveguides
- Antiprotons from Dark Matter: Effects of a Position-Dependent Diffusion Coefficient
- Position-dependent stochastic diffusion model of ion channel gating