skip to main content

Title: Stochastic linear multistep methods for the simulation of chemical kinetics

In this paper, we introduce the Stochastic Adams-Bashforth (SAB) and Stochastic Adams-Moulton (SAM) methods as an extension of the τ-leaping framework to past information. Using the Θ-trapezoidal τ-leap method of weak order two as a starting procedure, we show that the k-step SAB method with k ≥ 3 is order three in the mean and correlation, while a predictor-corrector implementation of the SAM method is weak order three in the mean but only order one in the correlation. These convergence results have been derived analytically for linear problems and successfully tested numerically for both linear and non-linear systems. A series of additional examples have been implemented in order to demonstrate the efficacy of this approach.
Authors:
 [1] ;  [2] ;  [3] ;  [4]
  1. Departamento de Informática, University of Valladolid, Valladolid (Spain)
  2. Department of Computer Science, University of Oxford, Oxford (United Kingdom)
  3. (Australia)
  4. School of Mathematical Sciences, Queensland University of Technology, Brisbane (Australia)
Publication Date:
OSTI Identifier:
22416091
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 142; Journal Issue: 6; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; CHEMICAL REACTION KINETICS; CHEMICAL REACTIONS; CORRELATIONS; IMPLEMENTATION; NONLINEAR PROBLEMS; STOCHASTIC PROCESSES