Efficient modified ``walk on spheres'' algorithm for the linearized Poisson-Bolzmann equation
Abstract
A discrete random walk method on grids was proposed and used to solve the linearized Poisson-Boltzmann equation (LPBE) [R. Ettelaie, J. Chem. Phys. 103, 3657 (1995)]. Here, we present an efficient grid-free random walk method. Based on a modified "walk on spheres" algorithm [B. S. Elepov and G. A. Mihailov, Sov. Math. Dokl. 14, 1276 (1973)] for the LPBE, this Monte Carlo algorithm uses a survival probability distribution function for the random walker in a continuous and free diffusion region. This simulation method is illustrated by computing four analytically solvable problems. In all cases, excellent agreement is observed.
- Publication:
-
Applied Physics Letters
- Pub Date:
- February 2001
- DOI:
- arXiv:
- arXiv:math-ph/0010032
- Bibcode:
- 2001ApPhL..78..787H
- Keywords:
-
- 05.60.Cd;
- 02.60.Lj;
- 02.50.Ng;
- Classical transport;
- Ordinary and partial differential equations;
- boundary value problems;
- Distribution theory and Monte Carlo studies;
- Mathematical Physics;
- Mathematics - Numerical Analysis
- E-Print:
- 12 pages, 5 figures