Transient subtraction: A control variate method for computing transport coefficients

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Monmarché, Pierre, Spacek, Renato, Stoltz, Gabriel
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913694789664768
author Monmarché, Pierre
Spacek, Renato
Stoltz, Gabriel
author_facet Monmarché, Pierre
Spacek, Renato
Stoltz, Gabriel
contents In molecular dynamics, transport coefficients measure the sensitivity of the invariant probability measure of the stochastic dynamics at hand with respect to some perturbation. They are typically computed using either the linear response of nonequilibrium dynamics, or the Green--Kubo formula. The estimators for both approaches have large variances, which motivates the study of variance reduction techniques for computing transport coefficients. We present an alternative approach, called the \emph{transient subtraction technique} (inspired by early work by Ciccotti and Jaccucci in 1975), which amounts to simulating a transient dynamics started off equilibrium and relaxing towards the equilibrium state, from which we subtract a sensibly coupled equilibrium trajectory, resulting in an estimator with smaller variance. We present the mathematical formulation of the transient subtraction technique, give error estimates on the bias and variance of the associated estimator, and demonstrate the relevance of the method through numerical illustrations for various systems.
format Preprint
id arxiv_https___arxiv_org_abs_2410_00212
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Transient subtraction: A control variate method for computing transport coefficients
Monmarché, Pierre
Spacek, Renato
Stoltz, Gabriel
Numerical Analysis
Statistical Mechanics
82C31, 65C40, 46N30
In molecular dynamics, transport coefficients measure the sensitivity of the invariant probability measure of the stochastic dynamics at hand with respect to some perturbation. They are typically computed using either the linear response of nonequilibrium dynamics, or the Green--Kubo formula. The estimators for both approaches have large variances, which motivates the study of variance reduction techniques for computing transport coefficients. We present an alternative approach, called the \emph{transient subtraction technique} (inspired by early work by Ciccotti and Jaccucci in 1975), which amounts to simulating a transient dynamics started off equilibrium and relaxing towards the equilibrium state, from which we subtract a sensibly coupled equilibrium trajectory, resulting in an estimator with smaller variance. We present the mathematical formulation of the transient subtraction technique, give error estimates on the bias and variance of the associated estimator, and demonstrate the relevance of the method through numerical illustrations for various systems.
title Transient subtraction: A control variate method for computing transport coefficients
topic Numerical Analysis
Statistical Mechanics
82C31, 65C40, 46N30
url https://arxiv.org/abs/2410.00212