Partially Reflected Diffusion

Abstract

The radiation (reactive or Robin) boundary condition for the diffusion equation is widely used in chemical and biological applications to express reactive boundaries. The underlying trajectories of the diffusing particles are believed to be partially absorbed and partially reflected at the reactive boundary; however, the relation between the reaction constant in the Robin boundary condition and the reflection probability is not well defined. In this paper we define the partially reflected process as a limit of the Markovian jump process generated by the Euler scheme for the underlying Ito dynamics with partial boundary reflection. Trajectories that cross the boundary are terminated with probability $P\sqrt{\Delta t}$ and otherwise are reflected in a normal or oblique direction. We use boundary layer analysis of the corresponding master equation to resolve the nonuniform convergence of the probability density function of the numerical scheme to the solution of the Fokker–Planck equation in a half-space, wi...