Abstract:Abstract We construct a higher-order adaptive method for strong approximations of exit times of Itô stochastic differential equations (SDEs). The method employs a strong Itô–Taylor scheme for simulating SDE paths, and adaptively decreases the step size in the numerical integration as the solution approaches the boundary of the domain. These techniques complement each other nicely: adaptive timestepping improves the accuracy of the exit time by reducing the magnitude of the overshoot of the numerical solution when it exits the domain, and higher-order schemes improve the approximation of the state of the diffusion process. We present two versions of the higher-order adaptive method. The first one uses the Milstein scheme as the numerical integrator and two step sizes for adaptive timestepping: $h$ when far away from the boundary and $h^2$ when close to the boundary. The second method is an extension of the first one using the strong Itô–Taylor scheme of order 1.5 as the numerical integrator and three step sizes for adaptive timestepping. Under some regularity assumptions, we show that for any $\xi>0$, the strong error is ${\mathcal{O}}(h^{1-\xi })$ and ${\mathcal{O}}(h^{3/2-\xi })$ for the first and second method, respectively. Provided quite restrictive commutativity conditions hold for the diffusion coefficient, we further show that the expected computational cost for both methods is ${\mathcal{O}}(h^{-1} \log (h^{-1}))$. This results in a near doubling/trebling of the strong error rate compared to the standard Euler–Maruyama-based approach, while the computational cost rate is kept close to order one. Numerical examples that support the theoretical results are provided, and we discuss the potential for extensions that would further improve the strong convergence rate of the method.

Highly Accurate Tau-Leaping Methods with Random Corrections.

Analysis of Explicit Tau-Leaping Schemes for Simulating Chemically Reacting Systems.

A Weak Second Order Tau-Leaping Method for Chemical Kinetic Systems.

Quasi-Monte Carlo methods applied to tau-leaping in stochastic biological systems

New Methods of Simulating Lévy Processes

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

The Weak Convergence Analysis of Tau-Leaping Methods: Revisited

Estimation of parameter sensitivities for stochastic reaction networks using tau-leap simulations

Robust Stochastic Chemical Reaction Networks and Bounded Tau-Leaping

A hybrid tau-leap for simulating chemical kinetics with applications to parameter estimation

$S$-Leaping: An adaptive, accelerated stochastic simulation algorithm, bridging $τ$-leaping and $R$-leaping

Beyond the adiabatic limit in systems with fast environments: a $τ$-leaping algorithm

Multi-level methods and approximating distribution functions

Central Limit Theorems of a Recursive Stochastic Algorithm with Applications to Adaptive Designs

A Numerical Truncation Approximation with A Posteriori Error Bounds for the Solution of Poisson's Equation

Chebyshev Methods with Discrete Noise: the Tau-Rock Methods

Adaptive Weak Approximation of Diffusions with Jumps

Stochastically accelerated perturbative triples correction in coupled cluster calculations

Numerical Solutions of Jump Diffusions with Markovian Switching

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions

Higher-order adaptive methods for exit times of Itô diffusions