Explicit Asymptotics on First Passage Times of Diffusion Processes
Angelos Dassios, Luting Li
TL;DR
This work develops a unified perturbation framework for explicit first passage time densities of time-homogeneous diffusions by formulating a killed Dirichlet problem in the Laplace domain and solving a recursive sequence in the frequency domain. The method yields closed-form asymptotic FPT densities for single-sided level crossings, with an exact probabilistic representation of truncation error and rigorous convergence properties that hold beyond small perturbations. Application to Ornstein–Uhlenbeck and Bessel processes provides concrete Nth-order densities expressed via parabolic-cylinder functions, along with detailed left/right tail behavior and error terms. Numerical experiments with downward OU FPTD demonstrate strong accuracy and substantial computational advantages over conventional inverse-Laplace techniques, while highlighting limitations in long-time tail regions for larger orders. The framework further extends to time-changed diffusions (e.g., exponential-Shiryaev/CIR-type models) and offers practical tools for survival analysis and finance where closed-form FPTs are valuable.
Abstract
We introduce a unified framework for solving first passage times of time-homogeneous diffusion processes. According to the killed version potential theory and the perturbation theory, we are able to deduce closed-form solutions for probability densities of single-sided level crossing problem. The framework is applicable to diffusion processes with continuous drift functions, and a recursive system in the frequency domain has been provided. Besides, we derive a probabilistic representation for error estimation. The representation can be used to evaluate deviations in perturbed density functions. In the present paper, we apply the framework to Ornstein-Uhlenbeck and Bessel processes to find closed-form approximations for their first passage times; another successful application is given by the exponential-Shiryaev process. Numerical results are provided at the end of this paper.
