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.

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

Using Stein’s Method to Analyze Euler–Maruyama Approximations of Regime-Switching Jump Diffusion Processes

Highly Accurate Tau-Leaping Methods with Random Corrections.

Variable-Step Euler–Maruyama Approximations of Regime-Switching Jump Diffusion Processes

The Weak Convergence Analysis of Tau-Leaping Methods: Revisited

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

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

Robust Stochastic Chemical Reaction Networks and Bounded Tau-Leaping

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

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

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

Compensated projected Euler-Maruyama method for stochastic differential equations with superlinear jumps

Stochastic analysis of chemical reactions in multi-component interacting systems at criticality

Effects of different discretisations of the Laplacian upon stochastic simulations of reaction-diffusion systems on both static and growing domains

Local convergence analysis of L1/finite element scheme for a constant delay reaction-subdiffusion equation with uniform time mesh

Convergence of a Metropolized Integrator for Stochastic Differential Equations with Variable Diffusion Coefficient

Two-scale Large Deviations for Chemical Reaction Kinetics Through Second Quantization Path Integral

How Discrete and Continuous Diffusion Meet: Comprehensive Analysis of Discrete Diffusion Models via a Stochastic Integral Framework

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

Convergence analysis of an explicit method and its random batch approximation for the McKean-Vlasov equations with non-globally Lipschitz conditions