Analytic and Monte Carlo approximations to the distribution of the first-passage time of drifted diffusion with stochastic resetting and mixed boundary conditions

Abstract

This paper introduces two techniques for computing the distribution of the absorption or first passage time of the drifted Wiener diffusion subject to Poisson resetting times to an upper hard wall barrier and to a lower absorbing barrier. The first method, which we call the Padé-partial fraction approximation, starts with the Padé approximation to the Laplace transform of the first passage time distribution, which is then exactly inverted by means of the partial fraction decomposition. The second method, which we call the multiresolution algorithm, is a Monte Carlo technique that exploits the properties of the Wiener process to generate Brownian bridges at increasing levels of resolution. Our numerical study reveals that the multiresolution algorithm has higher efficiency than standard Monte Carlo, whereas the faster Padé-partial fraction method is accurate in various circumstances and provides an analytical formula. Also, a closed-form exact expression for the expected first passage time is derived.