DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Global sensitivity analysis in stochastic simulators of uncertain reaction networks

Abstract

Stochastic models of chemical systems are often subjected to uncertainties in kinetic parameters in addition to the inherent random nature of their dynamics. Uncertainty quantification in such systems is generally achieved by means of sensitivity analyses in which one characterizes the variability with the uncertain kinetic parameters of the first statistical moments of model predictions. In this work, we propose an original global sensitivity analysis method where the parametric and inherent variability sources are both treated through Sobol’s decomposition of the variance into contributions from arbitrary subset of uncertain parameters and stochastic reaction channels. The conceptual development only assumes that the inherent and parametric sources are independent, and considers the Poisson processes in the random-time-change representation of the state dynamics as the fundamental objects governing the inherent stochasticity. Here, a sampling algorithm is proposed to perform the global sensitivity analysis, and to estimate the partial variances and sensitivity indices characterizing the importance of the various sources of variability and their interactions. The birth-death and Schlögl models are used to illustrate both the implementation of the algorithm and the richness of the proposed analysis method. The output of the proposed sensitivity analysis is also contrasted with a local derivative-based sensitivity analysismore » method classically used for this type of systems.« less

Authors:
 [1]; ORCiD logo [2];  [3]
  1. King Abdullah Univ. of Science and Technology (KAUST), Thuwal (Saudi Arabia). CEMSE Division
  2. King Abdullah Univ. of Science and Technology (KAUST), Thuwal (Saudi Arabia). CEMSE Division; Univ. Paris-Saclay, Gif-sur-Yvette (France). LIMSI and CNRS
  3. King Abdullah Univ. of Science and Technology (KAUST), Thuwal (Saudi Arabia). CEMSE Division; Duke Univ., Durham, NC (United States). Dept. of Mechanical Engineering and Materials Science
Publication Date:
Research Org.:
Duke Univ., Durham, NC (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); King Abdullah Univ. of Science and Technology (KAUST), Thuwal (Saudi Arabia)
OSTI Identifier:
1423885
Alternate Identifier(s):
OSTI ID: 1337433
Grant/Contract Number:  
SC0008789
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 145; Journal Issue: 24; Journal ID: ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Stochastic Process; Poisson's equation; Decomposition reactions; nonlinear dynamics; monte carlo methods; modeling; probability theory; real functions; mathematical physics; partial differential equations; reaction mechanisms

Citation Formats

Navarro Jimenez, M., Le Maître, O. P., and Knio, O. M. Global sensitivity analysis in stochastic simulators of uncertain reaction networks. United States: N. p., 2016. Web. doi:10.1063/1.4971797.
Navarro Jimenez, M., Le Maître, O. P., & Knio, O. M. Global sensitivity analysis in stochastic simulators of uncertain reaction networks. United States. https://doi.org/10.1063/1.4971797
Navarro Jimenez, M., Le Maître, O. P., and Knio, O. M. Fri . "Global sensitivity analysis in stochastic simulators of uncertain reaction networks". United States. https://doi.org/10.1063/1.4971797. https://www.osti.gov/servlets/purl/1423885.
@article{osti_1423885,
title = {Global sensitivity analysis in stochastic simulators of uncertain reaction networks},
author = {Navarro Jimenez, M. and Le Maître, O. P. and Knio, O. M.},
abstractNote = {Stochastic models of chemical systems are often subjected to uncertainties in kinetic parameters in addition to the inherent random nature of their dynamics. Uncertainty quantification in such systems is generally achieved by means of sensitivity analyses in which one characterizes the variability with the uncertain kinetic parameters of the first statistical moments of model predictions. In this work, we propose an original global sensitivity analysis method where the parametric and inherent variability sources are both treated through Sobol’s decomposition of the variance into contributions from arbitrary subset of uncertain parameters and stochastic reaction channels. The conceptual development only assumes that the inherent and parametric sources are independent, and considers the Poisson processes in the random-time-change representation of the state dynamics as the fundamental objects governing the inherent stochasticity. Here, a sampling algorithm is proposed to perform the global sensitivity analysis, and to estimate the partial variances and sensitivity indices characterizing the importance of the various sources of variability and their interactions. The birth-death and Schlögl models are used to illustrate both the implementation of the algorithm and the richness of the proposed analysis method. The output of the proposed sensitivity analysis is also contrasted with a local derivative-based sensitivity analysis method classically used for this type of systems.},
doi = {10.1063/1.4971797},
journal = {Journal of Chemical Physics},
number = 24,
volume = 145,
place = {United States},
year = {Fri Dec 23 00:00:00 EST 2016},
month = {Fri Dec 23 00:00:00 EST 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 8 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Polynomial chaos expansion for sensitivity analysis
journal, July 2009

  • Crestaux, Thierry; Le Maıˆtre, Olivier; Martinez, Jean-Marc
  • Reliability Engineering & System Safety, Vol. 94, Issue 7
  • DOI: 10.1016/j.ress.2008.10.008

Approximation of Population Processes
book, January 1981


Analysis of variance designs for model output
journal, March 1999


Uncertainty Quantification in MD Simulations. Part I: Forward Propagation
journal, January 2012

  • Rizzi, F.; Najm, H. N.; Debusschere, B. J.
  • Multiscale Modeling & Simulation, Vol. 10, Issue 4
  • DOI: 10.1137/110853169

Automatic identification of model reductions for discrete stochastic simulation
journal, July 2012

  • Wu, Sheng; Fu, Jin; Li, Hong
  • The Journal of Chemical Physics, Vol. 137, Issue 3
  • DOI: 10.1063/1.4733563

Preconditioned Bayesian Regression for Stochastic Chemical Kinetics
journal, July 2013

  • Alexanderian, Alen; Rizzi, Francesco; Rathinam, Muruhan
  • Journal of Scientific Computing, Vol. 58, Issue 3
  • DOI: 10.1007/s10915-013-9745-5

Variance decomposition in stochastic simulators
journal, June 2015

  • Le Maître, O. P.; Knio, O. M.; Moraes, A.
  • The Journal of Chemical Physics, Vol. 142, Issue 24
  • DOI: 10.1063/1.4922922

Binomial distribution based τ-leap accelerated stochastic simulation
journal, January 2005

  • Chatterjee, Abhijit; Vlachos, Dionisios G.; Katsoulakis, Markos A.
  • The Journal of Chemical Physics, Vol. 122, Issue 2
  • DOI: 10.1063/1.1833357

Hybrid Chernoff Tau-Leap
journal, January 2014

  • Moraes, Alvaro; Tempone, Raul; Vilanova, Pedro
  • Multiscale Modeling & Simulation, Vol. 12, Issue 2
  • DOI: 10.1137/130925657

Efficient computation of parameter sensitivities of discrete stochastic chemical reaction networks
journal, January 2010

  • Rathinam, Muruhan; Sheppard, Patrick W.; Khammash, Mustafa
  • The Journal of Chemical Physics, Vol. 132, Issue 3
  • DOI: 10.1063/1.3280166

Uncertainty quantification in MD simulations of concentration driven ionic flow through a silica nanopore. I. Sensitivity to physical parameters of the pore
journal, May 2013

  • Rizzi, F.; Jones, R. E.; Debusschere, B. J.
  • The Journal of Chemical Physics, Vol. 138, Issue 19
  • DOI: 10.1063/1.4804666

Sensitivity Analysis in Chemical Kinetics
journal, October 1983


Reduction of chemical reaction networks through delay distributions
journal, March 2013

  • Barrio, Manuel; Leier, André; Marquez-Lago, Tatiana T.
  • The Journal of Chemical Physics, Vol. 138, Issue 10
  • DOI: 10.1063/1.4793982

Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels
journal, March 2000

  • Gibson, Michael A.; Bruck, Jehoshua
  • The Journal of Physical Chemistry A, Vol. 104, Issue 9
  • DOI: 10.1021/jp993732q

Uncertainty quantification in MD simulations of concentration driven ionic flow through a silica nanopore. II. Uncertain potential parameters
journal, May 2013

  • Rizzi, F.; Jones, R. E.; Debusschere, B. J.
  • The Journal of Chemical Physics, Vol. 138, Issue 19
  • DOI: 10.1063/1.4804669

Multilevel Monte Carlo for Continuous Time Markov Chains, with Applications in Biochemical Kinetics
journal, January 2012

  • Anderson, David F.; Higham, Desmond J.
  • Multiscale Modeling & Simulation, Vol. 10, Issue 1
  • DOI: 10.1137/110840546

A pathwise derivative approach to the computation of parameter sensitivities in discrete stochastic chemical systems
journal, January 2012

  • Sheppard, Patrick W.; Rathinam, Muruhan; Khammash, Mustafa
  • The Journal of Chemical Physics, Vol. 136, Issue 3
  • DOI: 10.1063/1.3677230

Spectral stochastic uncertainty quantification in chemical systems
journal, September 2004

  • Reagan, M. T.; Najm 4, H. N.; Debusschere, B. J.
  • Combustion Theory and Modelling, Vol. 8, Issue 3
  • DOI: 10.1088/1364-7830/8/3/010

Error analysis of tau-leap simulation methods
journal, December 2011

  • Anderson, David F.; Ganguly, Arnab; Kurtz, Thomas G.
  • The Annals of Applied Probability, Vol. 21, Issue 6
  • DOI: 10.1214/10-AAP756

Constructing stochastic models from deterministic process equations by propensity adjustment
journal, January 2011

  • Wu, Jialiang; Vidakovic, Brani; Voit, Eberhard O.
  • BMC Systems Biology, Vol. 5, Issue 1
  • DOI: 10.1186/1752-0509-5-187

Stochastic Finite Elements: A Spectral Approach
book, January 1991


PC analysis of stochastic differential equations driven by Wiener noise
journal, March 2015


Hybrid pathwise sensitivity methods for discrete stochastic models of chemical reaction systems
journal, January 2015

  • Wolf, Elizabeth Skubak; Anderson, David F.
  • The Journal of Chemical Physics, Vol. 142, Issue 3
  • DOI: 10.1063/1.4905332

Efficient step size selection for the tau-leaping simulation method
journal, January 2006

  • Cao, Yang; Gillespie, Daniel T.; Petzold, Linda R.
  • The Journal of Chemical Physics, Vol. 124, Issue 4
  • DOI: 10.1063/1.2159468

An adaptive multi-level simulation algorithm for stochastic biological systems
journal, January 2015

  • Lester, C.; Yates, C. A.; Giles, M. B.
  • The Journal of Chemical Physics, Vol. 142, Issue 2
  • DOI: 10.1063/1.4904980

Exact stochastic simulation of coupled chemical reactions
journal, December 1977

  • Gillespie, Daniel T.
  • The Journal of Physical Chemistry, Vol. 81, Issue 25
  • DOI: 10.1021/j100540a008

Multilevel Monte Carlo Path Simulation
journal, June 2008


Spectral Methods for Uncertainty Quantification: With Applications to Computational Fluid Dynamics
book, January 2010


Making best use of model evaluations to compute sensitivity indices
journal, May 2002


The slow-scale linear noise approximation: an accurate, reduced stochastic description of biochemical networks under timescale separation conditions
journal, January 2012

  • Thomas, Philipp; Straube, Arthur V.; Grima, Ramon
  • BMC Systems Biology, Vol. 6, Issue 1
  • DOI: 10.1186/1752-0509-6-39

Derivative based global sensitivity measures and their link with global sensitivity indices
journal, June 2009


Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method
journal, December 2003

  • Rathinam, Muruhan; Petzold, Linda R.; Cao, Yang
  • The Journal of Chemical Physics, Vol. 119, Issue 24
  • DOI: 10.1063/1.1627296

Master Equations, Stochastic Differential Equations, and Bifurcations
journal, February 1979


Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index
journal, February 2010

  • Saltelli, Andrea; Annoni, Paola; Azzini, Ivano
  • Computer Physics Communications, Vol. 181, Issue 2
  • DOI: 10.1016/j.cpc.2009.09.018

The subtle business of model reduction for stochastic chemical kinetics
journal, February 2009

  • Gillespie, Dan T.; Cao, Yang; Sanft, Kevin R.
  • The Journal of Chemical Physics, Vol. 130, Issue 6
  • DOI: 10.1063/1.3072704

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions
journal, December 1976


A modified next reaction method for simulating chemical systems with time dependent propensities and delays
journal, December 2007

  • Anderson, David F.
  • The Journal of Chemical Physics, Vol. 127, Issue 21
  • DOI: 10.1063/1.2799998

Sensitivity Analysis for Chemical Models
journal, July 2005

  • Saltelli, Andrea; Ratto, Marco; Tarantola, Stefano
  • Chemical Reviews, Vol. 105, Issue 7
  • DOI: 10.1021/cr040659d

A finite difference method for estimating second order parameter sensitivities of discrete stochastic chemical reaction networks
journal, December 2012

  • Wolf, Elizabeth Skubak; Anderson, David F.
  • The Journal of Chemical Physics, Vol. 137, Issue 22
  • DOI: 10.1063/1.4770052

Works referencing / citing this record:

Sensitivity-Driven Adaptive Construction of Reduced-space Surrogates
journal, December 2018

  • Vohra, Manav; Alexanderian, Alen; Safta, Cosmin
  • Journal of Scientific Computing, Vol. 79, Issue 2
  • DOI: 10.1007/s10915-018-0894-4