Abstract:We aim to construct higher order tau-leaping methods for numerically simulating stochastic chemical kinetic systems in this paper. By adding a random correction to the primitive tau-leaping scheme in each time step, we greatly improve the accuracy of the tau-leaping approximations. This gain in accuracy actually comes from the reduction in the local truncation error of the scheme in the order of tau, the marching time step size. While the local truncation error of the primitive tau-leaping method is O(tau(2)) for all moments, our Poisson random correction tau-leaping method, in which the correction term is a Poisson random variable, can reduce the local truncation error for the mean to O(tau(3)), and both Gaussian random correction tau-leaping methods, in which the correction term is a Gaussian random variable, can reduce the local truncation error for both the mean and covariance to O(tau(3)). Numerical results demonstrate that these novel methods more accurately capture crucial properties such as the mean and variance than existing methods for simulating chemical reaction systems. This work constitutes a first step to construct high order numerical methods for simulating jump processes. With further refinement and appropriately modified step-size selection procedures, the random correction methods should provide a viable way of simulating chemical reaction systems accurately and efficiently.

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

An Approach to Numerically Solving the Poisson Equation

Numerical Solutions of Jump Diffusions with Markovian Switching

The truncated EM method for stochastic differential equations with Poisson jumps

From Monte Carlo to neural networks approximations of boundary value problems

Explicit and implicit error inhibiting schemes with post-processing

Numerical approximation of a hybrid Poisson-jump Ait-Sahalia-type interest rate model with delay

First order strong approximation of Ait-Sahalia-type interest rate model with Poisson jumps

Approximation schemes for McKean-Vlasov and Boltzmann type equations (error analyses in total variation distance)

Numerical method to solve impulse control problems for partially observed piecewise deterministic Markov processes

On the asymptotic stability and numerical analysis of solutions to nonlinear stochastic differential equations with jumps

Estimation of domain truncation error for a system of PDEs arising in option pricing

On the Poisson equation for nonreversible Markov jump processes

Highly Accurate Tau-Leaping Methods with Random Corrections.

Probabilistic limit behaviors of numerical discretizations for time-homogeneous Markov processes

Reliable numerical solution of a class of nonlinear elliptic problems generated by the Poisson-Boltzmann equation

Convergence of the Semi-implicit Euler Numerical Solutions of Stochastic Jump-diffusion Equations with Distributed Memory

Recursive Optimal Stopping with Poisson Stopping Constraints

Strongly convergent homogeneous approximations to inhomogeneous Markov jump processes and applications

Numerical Method In Reproducing Kernel Space For An Inverse Source Problem For The Fractional Diffusion Equation

On Brownian motion approximation of compound Poisson processes with applications to threshold models