Optimal explicit stabilized postprocessed τ-leap method for the simulation of chemical kinetics

06/17/2021
by   Assyr Abdulle, et al.
0

The simulation of chemical kinetics involving multiple scales constitutes a modeling challenge (from ordinary differential equations to Markov chain) and a computational challenge (multiple scales, large dynamical systems, time step restrictions). In this paper we propose a new discrete stochastic simulation algorithm: the postprocessed second kind stabilized orthogonal τ-leap Runge-Kutta method (PSK-τ-ROCK). In the context of chemical kinetics this method can be seen as a stabilization of Gillespie's explicit τ-leap combined with a postprocessor. The stabilized procedure allows to simulate problems with multiple scales (stiff), while the postprocessing procedure allows to approximate the invariant measure (e.g. mean and variance) of ergodic stochastic dynamical systems. We prove stability and accuracy of the PSK-τ-ROCK. Numerical experiments illustrate the high reliability and efficiency of the scheme when compared to other τ-leap methods.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

07/24/2019

Approximate Numerical Integration of the Chemical Master Equation for Stochastic Reaction Networks

Numerical solution of the chemical master equation for stochastic reacti...
08/31/2021

Sensitivity Approximation by the Peano-Baker Series

In this paper we develop a new method for numerically approximating sens...
09/10/2019

A class of robust numerical methods for solving dynamical systems with multiple time scales

In this paper, we develop a class of robust numerical methods for solvin...
05/18/2021

Learning stochastic dynamical systems with neural networks mimicking the Euler-Maruyama scheme

Stochastic differential equations (SDEs) are one of the most important r...
09/01/2020

Variance Reduction with Array-RQMC for Tau-Leaping Simulation of Stochastic Biological and Chemical Reaction Networks

We explore the use of Array-RQMC, a randomized quasi-Monte Carlo method ...
09/28/2020

Parameter Experimental Analysis of the Reservoirs Observers using Echo State Network Approach

Dynamical systems has a variety of applications for the new information ...
10/20/2020

Algorithmic Reduction of Biological Networks With Multiple Time Scales

We present a symbolic algorithmic approach that allows to compute invari...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.