Probabilistic interpretation and random walk on spheres algorithms for the Poisson-Boltzmann equation in molecular dynamics
Motivated by the development of efficient Monte Carlo methods for PDE models in molecular dynamics, we establish a new probabilistic interpretation of a family of divergence form operators with discontinuous coefficients at the interface of two open subsets of . This family of operators includes the case of the linearized Poisson-Boltzmann equation used to compute the electrostatic free energy of a molecule. More precisely, we explicitly construct a Markov process whose infinitesimal generator...