Brownian Motion of Radioactive Particles: Derivation and Monte Carlo Test of Spatial and Temporal Distributions

Stochastic processes such as diffusion can be analyzed by means of a partial differential equation of the Fokker-Planck type (FPE), which yields a transition probability density, or by a stochastic differential equation of the Langevin type (LE), which yields the time evolution of a statistical process variable. Provided the stochastic process is continuous and certain boundary conditions are met, the two approaches yield equivalent information. However, Brownian motion of radioactively decaying particles is not a continuous process because the Brownian trajectories abruptly terminate when the particle decays. Recent analysis of the Brownian motion of decaying particles by both approaches has led to different mean-square displacements. In this paper, we demonstrate the complete equivalence of the two approaches by 1) showing quantitatively and operationally how the probability densities and statistical moments predicted by the FPE and LE relate to one another, 2) verifying that both approaches lead to identical statistical moments at all orders, and 3) confirming that the analytical solution to the FPE accurately describes the Brownian trajectories obtained by Monte Carlo simulations based on the LE. The analysis in this paper addresses both the spatial distribution of the particles (i.e. the question of displacement as a function of diffusion time) and the temporal distribution (i.e. the question of first-passage time to fixed absorbing boundaries).