Multilevel Hybrid Split Step Implicit Tau-Leap

02/23/2020
by   Chiheb Ben Hammouda, et al.
0

In biochemically reactive systems with small copy numbers of one or more reactant molecules, the dynamics are dominated by stochastic effects. To approximate those systems, discrete state-space and stochastic simulation approaches have been shown to be more relevant than continuous state-space and deterministic ones. These stochastic models constitute the theory of Stochastic Reaction Networks (SRNs). In systems characterized by having simultaneously fast and slow timescales, existing discrete space-state stochastic path simulation methods, such as the stochastic simulation algorithm (SSA) and the explicit tau-leap (explicit-TL) method, can be very slow. In this talk, we propose a novel implicit scheme, split-step implicit tau-leap (SSI-TL), to improve numerical stability and provide efficient simulation algorithms for those systems. Furthermore, to estimate statistical quantities related to SRNs, we propose a novel hybrid Multilevel Monte Carlo (MLMC) estimator in the spirit of the work by Anderson and Higham (SIAM Multiscal Model. Simul. 10(1), 2012). This estimator uses the SSI-TL scheme at levels where the explicit-TL method is not applicable due to numerical stability issues, and then, starting from a certain interface level, it switches to the explicit scheme. We present numerical examples that illustrate the achieved gains of our proposed approach in this context.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/14/2019

Importance sampling for a robust and efficient multilevel Monte Carlo estimator for stochastic reaction networks

The multilevel Monte Carlo (MLMC) method for continuous time Markov chai...
research
10/01/2018

Multilevel Adaptive Sparse Grid Quadrature for Monte Carlo models

Many problems require to approximate an expected value by some kind of M...
research
03/07/2023

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

We develop new multilevel Monte Carlo (MLMC) methods to estimate the exp...
research
05/20/2021

A flexible split-step scheme for MV-SDEs

We present an implicit Split-Step explicit Euler type Method (dubbed SSM...
research
07/20/2019

Drift-preserving numerical integrators for stochastic Hamiltonian systems

The paper deals with numerical discretizations of separable nonlinear Ha...
research
09/07/2022

Multilevel Path Branching for Digital Options

We propose a new Monte Carlo-based estimator for digital options with as...
research
02/28/2021

Accelerated Jarzynski Estimator with Deterministic Virtual Trajectories

The Jarzynski estimator is a powerful tool that uses nonequilibrium stat...

Please sign up or login with your details

Forgot password? Click here to reset