Voir la notice de l'article provenant de la source Numdam
For a stationary Markov process the detailed balance condition is equivalent to the time-reversibility of the process. For stochastic differential equations (SDE's), the time discretization of numerical schemes usually destroys the time-reversibility property. Despite an extensive literature on the numerical analysis for SDE's, their stability properties, strong and/or weak error estimates, large deviations and infinite-time estimates, no quantitative results are known on the lack of reversibility of discrete-time approximation processes. In this paper we provide such quantitative estimates by using the concept of entropy production rate, inspired by ideas from non-equilibrium statistical mechanics. The entropy production rate for a stochastic process is defined as the relative entropy (per unit time) of the path measure of the process with respect to the path measure of the time-reversed process. By construction the entropy production rate is nonnegative and it vanishes if and only if the process is reversible. Crucially, from a numerical point of view, the entropy production rate is an a posteriori quantity, hence it can be computed in the course of a simulation as the ergodic average of a certain functional of the process (the so-called Gallavotti-Cohen (GC) action functional). We compute the entropy production for various numerical schemes such as explicit Euler-Maruyama and explicit Milstein's for reversible SDEs with additive or multiplicative noise. In addition we analyze the entropy production for the BBK integrator for the Langevin equation. The order (in the time-discretization step Δt) of the entropy production rate provides a tool to classify numerical schemes in terms of their (discretization-induced) irreversibility. Our results show that the type of the noise critically affects the behavior of the entropy production rate. As a striking example of our results we show that the Euler scheme for multiplicative noise is not an adequate scheme from a reversibility point of view since its entropy production rate does not decrease with Δt.
@article{M2AN_2014__48_5_1351_0,
author = {Katsoulakis, Markos and Pantazis, Yannis and Rey-Bellet, Luc},
title = {Measuring the {Irreversibility} of {Numerical} {Schemes} for {Reversible} {Stochastic} {Differential} {Equations}},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis },
pages = {1351--1379},
publisher = {EDP-Sciences},
volume = {48},
number = {5},
year = {2014},
doi = {10.1051/m2an/2013142},
mrnumber = {3264357},
language = {en},
url = {http://geodesic.mathdoc.fr/articles/10.1051/m2an/2013142/}
}
TY - JOUR AU - Katsoulakis, Markos AU - Pantazis, Yannis AU - Rey-Bellet, Luc TI - Measuring the Irreversibility of Numerical Schemes for Reversible Stochastic Differential Equations JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2014 SP - 1351 EP - 1379 VL - 48 IS - 5 PB - EDP-Sciences UR - http://geodesic.mathdoc.fr/articles/10.1051/m2an/2013142/ DO - 10.1051/m2an/2013142 LA - en ID - M2AN_2014__48_5_1351_0 ER -
%0 Journal Article %A Katsoulakis, Markos %A Pantazis, Yannis %A Rey-Bellet, Luc %T Measuring the Irreversibility of Numerical Schemes for Reversible Stochastic Differential Equations %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2014 %P 1351-1379 %V 48 %N 5 %I EDP-Sciences %U http://geodesic.mathdoc.fr/articles/10.1051/m2an/2013142/ %R 10.1051/m2an/2013142 %G en %F M2AN_2014__48_5_1351_0
Katsoulakis, Markos; Pantazis, Yannis; Rey-Bellet, Luc. Measuring the Irreversibility of Numerical Schemes for Reversible Stochastic Differential Equations. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 48 (2014) no. 5, pp. 1351-1379. doi: 10.1051/m2an/2013142
Cité par Sources :