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

Title: Uncertainty quantification for generalized Langevin dynamics

Abstract

Here, we present efficient finite difference estimators for goal-oriented sensitivity indices with applications to the generalized Langevin equation (GLE). In particular, we apply these estimators to analyze an extended variable formulation of the GLE where other well known sensitivity analysis techniques such as the likelihood ratio method are not applicable to key parameters of interest. These easily implemented estimators are formed by coupling the nominal and perturbed dynamics appearing in the finite difference through a common driving noise or common random path. After developing a general framework for variance reduction via coupling, we demonstrate the optimality of the common random path coupling in the sense that it produces a minimal variance surrogate for the difference estimator relative to sampling dynamics driven by independent paths. In order to build intuition for the common random path coupling, we evaluate the efficiency of the proposed estimators for a comprehensive set of examples of interest in particle dynamics. These reduced variance difference estimators are also a useful tool for performing global sensitivity analysis and for investigating non-local perturbations of parameters, such as increasing the number of Prony modes active in an extended variable GLE.

Authors:
ORCiD logo [1];  [1]; ORCiD logo [1]
  1. Univ. of Massachusetts, Amherst, MA (United States). Dept. of Mathematics and Statistics
Publication Date:
Research Org.:
Univ. of Massachusetts, Amherst, MA (United States). Dept. of Mathematics and Statistics
Sponsoring Org.:
USDOE Office of Science (SC); National Science Foundation (NSF)
OSTI Identifier:
1465368
Alternate Identifier(s):
OSTI ID: 1336482
Grant/Contract Number:  
SC0010723; DMS-1515712
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Chemical Physics
Additional Journal Information:
Journal Volume: 145; Journal Issue: 22; Journal ID: ISSN 0021-9606
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING

Citation Formats

Hall, Eric J., Katsoulakis, Markos A., and Rey-Bellet, Luc. Uncertainty quantification for generalized Langevin dynamics. United States: N. p., 2016. Web. doi:10.1063/1.4971433.
Hall, Eric J., Katsoulakis, Markos A., & Rey-Bellet, Luc. Uncertainty quantification for generalized Langevin dynamics. United States. https://doi.org/10.1063/1.4971433
Hall, Eric J., Katsoulakis, Markos A., and Rey-Bellet, Luc. Wed . "Uncertainty quantification for generalized Langevin dynamics". United States. https://doi.org/10.1063/1.4971433. https://www.osti.gov/servlets/purl/1465368.
@article{osti_1465368,
title = {Uncertainty quantification for generalized Langevin dynamics},
author = {Hall, Eric J. and Katsoulakis, Markos A. and Rey-Bellet, Luc},
abstractNote = {Here, we present efficient finite difference estimators for goal-oriented sensitivity indices with applications to the generalized Langevin equation (GLE). In particular, we apply these estimators to analyze an extended variable formulation of the GLE where other well known sensitivity analysis techniques such as the likelihood ratio method are not applicable to key parameters of interest. These easily implemented estimators are formed by coupling the nominal and perturbed dynamics appearing in the finite difference through a common driving noise or common random path. After developing a general framework for variance reduction via coupling, we demonstrate the optimality of the common random path coupling in the sense that it produces a minimal variance surrogate for the difference estimator relative to sampling dynamics driven by independent paths. In order to build intuition for the common random path coupling, we evaluate the efficiency of the proposed estimators for a comprehensive set of examples of interest in particle dynamics. These reduced variance difference estimators are also a useful tool for performing global sensitivity analysis and for investigating non-local perturbations of parameters, such as increasing the number of Prony modes active in an extended variable GLE.},
doi = {10.1063/1.4971433},
journal = {Journal of Chemical Physics},
number = 22,
volume = 145,
place = {United States},
year = {Wed Dec 14 00:00:00 EST 2016},
month = {Wed Dec 14 00:00:00 EST 2016}
}

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

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

Save / Share:

Works referenced in this record:

Sensitivity Anaysis as an Ingredient of Modeling
journal, November 2000

  • Campolongo, F.; Tarantola, S.; Saltelli, A.
  • Statistical Science, Vol. 15, Issue 4
  • DOI: 10.1214/ss/1009213004

Mori–Zwanzig formalism as a practical computational tool
journal, January 2010

  • Hijón, Carmen; Español, Pep; Vanden-Eijnden, Eric
  • Faraday Discuss., Vol. 144
  • DOI: 10.1039/B902479B

Applications of the generalized Langevin equation: Towards a realistic description of the baths
journal, January 2015


Parametric sensitivity analysis for stochastic molecular systems using information theoretic metrics
journal, July 2015

  • Tsourtis, Anastasios; Pantazis, Yannis; Katsoulakis, Markos A.
  • The Journal of Chemical Physics, Vol. 143, Issue 1
  • DOI: 10.1063/1.4922924

Convergence Rates of Finite-Difference Sensitivity Estimates for Stochastic Systems
journal, August 1993


Anomalous diffusion: Exact solution of the generalized Langevin equation for harmonically bounded particle
journal, January 2006


Factorial Sampling Plans for Preliminary Computational Experiments
journal, May 1991


Monte Carlo methods for security pricing
journal, June 1997


Applications of Malliavin calculus to Monte Carlo methods in finance
journal, August 1999

  • Fournié, Eric; Lasry, Jean-Michel; Lebuchoux, Jérôme
  • Finance and Stochastics, Vol. 3, Issue 4
  • DOI: 10.1007/s007800050068

A multivariate generalization of Prony's method
journal, February 2016


SPSens: a software package for stochastic parameter sensitivity analysis of biochemical reaction networks
journal, October 2012


A relative entropy rate method for path space sensitivity analysis of stationary complex stochastic dynamics
journal, February 2013

  • Pantazis, Yannis; Katsoulakis, Markos A.
  • The Journal of Chemical Physics, Vol. 138, Issue 5
  • DOI: 10.1063/1.4789612

Path-Space Information Bounds for Uncertainty Quantification and Sensitivity Analysis of Stochastic Dynamics
journal, January 2016

  • Dupuis, Paul; Katsoulakis, Markos A.; Pantazis, Yannis
  • SIAM/ASA Journal on Uncertainty Quantification, Vol. 4, Issue 1
  • DOI: 10.1137/15M1025645

Some Guidelines and Guarantees for Common Random Numbers
journal, June 1992


Approximation by exponential sums revisited
journal, March 2010


Accelerating the convergence of path integral dynamics with a generalized Langevin equation
journal, February 2011

  • Ceriotti, Michele; Manolopoulos, David E.; Parrinello, Michele
  • The Journal of Chemical Physics, Vol. 134, Issue 8
  • DOI: 10.1063/1.3556661

Colored-Noise Thermostats à la Carte
journal, March 2010

  • Ceriotti, Michele; Bussi, Giovanni; Parrinello, Michele
  • Journal of Chemical Theory and Computation, Vol. 6, Issue 4
  • DOI: 10.1021/ct900563s

Numerical Experiments on the Stochastic Behavior of a Lennard-Jones Gas System
journal, September 1973


Likelihood ratio gradient estimation for stochastic systems
journal, October 1990


Nonergodicity Mimics Inhomogeneity in Single Particle Tracking
journal, June 2008


Nuclear Quantum Effects in Solids Using a Colored-Noise Thermostat
journal, July 2009


Subdiffusive behavior in a trapping potential: Mean square displacement and velocity autocorrelation function
journal, August 2009


On the Convergence Rates of IPA and FDC Derivative Estimators
journal, August 1994


A simple framework to justify linear response theory
journal, March 2010


Direct Computation of Dynamical Response by Molecular Dynamics: The Mobility of a Charged Lennard-Jones Particle
journal, September 1975


Efficient multiple time scale molecular dynamics: Using colored noise thermostats to stabilize resonances
journal, January 2011

  • Morrone, Joseph A.; Markland, Thomas E.; Ceriotti, Michele
  • The Journal of Chemical Physics, Vol. 134, Issue 1
  • DOI: 10.1063/1.3518369

Statistical challenges in microrheology: STATISTICAL CHALLENGES IN MICRORHEOLOGY
journal, May 2012


Malliavin Greeks without Malliavin calculus
journal, November 2007


Accelerated Sensitivity Analysis in High-Dimensional Stochastic Reaction Networks
journal, July 2015


Langevin Equation with Colored Noise for Constant-Temperature Molecular Dynamics Simulations
journal, January 2009


On approximation of functions by exponential sums
journal, July 2005


Asymptotic analysis for the generalized Langevin equation
journal, April 2011


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

Comparison of finite difference based methods to obtain sensitivities of stochastic chemical kinetic models
journal, February 2013

  • Srivastava, Rishi; Anderson, David F.; Rawlings, James B.
  • The Journal of Chemical Physics, Vol. 138, Issue 7
  • DOI: 10.1063/1.4790650

Numerical integration of the extended variable generalized Langevin equation with a positive Prony representable memory kernel
journal, July 2013

  • Baczewski, Andrew D.; Bond, Stephen D.
  • The Journal of Chemical Physics, Vol. 139, Issue 4
  • DOI: 10.1063/1.4815917

Optical Measurements of Frequency-Dependent Linear Viscoelastic Moduli of Complex Fluids
journal, February 1995


Problem reduction, renormalization, and memory
journal, January 2006

  • Chorin, Alexandre; Stinis, Panagiotis
  • Communications in Applied Mathematics and Computational Science, Vol. 1, Issue 1
  • DOI: 10.2140/camcos.2006.1.1

Multilevel hybrid split-step implicit tau-leap
journal, June 2016


Generalized Langevin Equation with Fractional Gaussian Noise: Subdiffusion within a Single Protein Molecule
journal, October 2004


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

On the Malliavin approach to Monte Carlo approximation of conditional expectations
journal, January 2004


Goal-oriented sensitivity analysis for lattice kinetic Monte Carlo simulations
journal, March 2014

  • Arampatzis, Georgios; Katsoulakis, Markos A.
  • The Journal of Chemical Physics, Vol. 140, Issue 12
  • DOI: 10.1063/1.4868649

Efficient stochastic thermostatting of path integral molecular dynamics
journal, September 2010

  • Ceriotti, Michele; Parrinello, Michele; Markland, Thomas E.
  • The Journal of Chemical Physics, Vol. 133, Issue 12
  • DOI: 10.1063/1.3489925

Incorporation of memory effects in coarse-grained modeling via the Mori-Zwanzig formalism
journal, December 2015

  • Li, Zhen; Bian, Xin; Li, Xiantao
  • The Journal of Chemical Physics, Vol. 143, Issue 24
  • DOI: 10.1063/1.4935490

An effective screening design for sensitivity analysis of large models
journal, October 2007

  • Campolongo, Francesca; Cariboni, Jessica; Saltelli, Andrea
  • Environmental Modelling & Software, Vol. 22, Issue 10
  • DOI: 10.1016/j.envsoft.2006.10.004

An Efficient Finite Difference Method for Parameter Sensitivities of Continuous Time Markov Chains
journal, January 2012

  • Anderson, David F.
  • SIAM Journal on Numerical Analysis, Vol. 50, Issue 5
  • DOI: 10.1137/110849079

Parameter estimation for exponential sums by approximate Prony method
journal, May 2010


Applications of Malliavin calculus to Monte-Carlo methods in finance. II
journal, April 2001

  • Fournié, Eric; Lasry, Jean-Michel; Lebuchoux, Jérôme
  • Finance and Stochastics, Vol. 5, Issue 2
  • DOI: 10.1007/PL00013529

Model Comparison and Assessment for Single Particle Tracking in Biological Fluids
journal, October 2016

  • Lysy, Martin; Pillai, Natesh S.; Hill, David B.
  • Journal of the American Statistical Association, Vol. 111, Issue 516
  • DOI: 10.1080/01621459.2016.1158716

Convergence rates for steady-state derivative estimators
journal, December 1992

  • L'Ecuyer, Pierre
  • Annals of Operations Research, Vol. 39, Issue 1
  • DOI: 10.1007/bf02060938

Parametric sensitivity analysis for biochemical reaction networks based on pathwise information theory
journal, January 2013

  • Pantazis, Yannis; Katsoulakis, Markos A.; Vlachos, Dionisios G.
  • BMC Bioinformatics, Vol. 14, Issue 1
  • DOI: 10.1186/1471-2105-14-311

Computing generalized Langevin equations and generalized Fokker-Planck equations
journal, June 2009

  • Darve, E.; Solomon, J.; Kia, A.
  • Proceedings of the National Academy of Sciences, Vol. 106, Issue 27
  • DOI: 10.1073/pnas.0902633106

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems
journal, April 2014


Efficient estimators for likelihood ratio sensitivity indices of complex stochastic dynamics
journal, March 2016

  • Arampatzis, Georgios; Katsoulakis, Markos A.; Rey-Bellet, Luc
  • The Journal of Chemical Physics, Vol. 144, Issue 10
  • DOI: 10.1063/1.4943388

Generalized Langevin equation for solids. I. Rigorous derivation and main properties
journal, September 2008


Pathwise Sensitivity Analysis in Transient Regimes
book, January 2015


Open Classical Systems
book, January 2006


Works referencing / citing this record:

Homogenization for Generalized Langevin Equations with Applications to Anomalous Diffusion
journal, February 2020


Fine tuning classical and quantum molecular dynamics using a generalized Langevin equation
journal, March 2018

  • Rossi, Mariana; Kapil, Venkat; Ceriotti, Michele
  • The Journal of Chemical Physics, Vol. 148, Issue 10
  • DOI: 10.1063/1.4990536