Sticky coupling as a control variate for sensitivity analysis

Shiva Darshan,Andreas Eberle,Gabriel Stoltz

2024-09-24

Abstract:We present and analyze a control variate strategy based on couplings to reduce the variance of finite difference estimators of sensitivity coefficients, called transport coefficients in the physics literature. We study the bias and variance of a sticky-coupling and a synchronous-coupling based estimator as the finite difference parameter $\eta$ goes to zero. For diffusions with elliptic additive noise, we show that when the drift is contractive outside a compact the bias of a sticky-coupling based estimator is bounded as $\eta \to 0$ and its variance behaves like $\eta^{-1}$, compared to the standard estimator whose bias and variance behave like $\eta^{-1}$ and $\eta^{-2}$, respectively. Under the stronger assumption that the drift is contractive everywhere, we additionally show that the bias and variance of the synchronous-coupling based estimator are both bounded as $\eta \to 0$. Our hypotheses include overdamped Langevin dynamics with many physically relevant non-convex potentials. We illustrate our theoretical results with numerical examples, including overdamped Langevin dynamics with a highly non-convex Lennard-Jones potential to demonstrate both failure of synchronous coupling and the effectiveness of sticky coupling in the not globally contractive setting.

Probability,Numerical Analysis,Statistics Theory,Computation

What problem does this paper attempt to address?

The paper aims to address the problem of calculating transport coefficients in molecular dynamics. Specifically, the authors propose a control variable method based on a coupling strategy to reduce the variance of finite difference estimators when calculating sensitivity coefficients (i.e., transport coefficients in physical literature). The core contributions of the paper are:

1. **Proposing the sticky coupling strategy**: This coupling method can effectively reduce estimation bias and variance, especially in non-globally contracting scenarios.
2. **Comparing with synchronous coupling**: Synchronous coupling performs well when dynamics are strongly contracting in all regions; however, in non-globally contracting scenarios, sticky coupling is more effective.

Through theoretical analysis and numerical experiments, the authors demonstrate the performance of these two coupling methods under different conditions and prove that under certain assumptions, sticky coupling can effectively reduce bias while maintaining low variance. These results are significant for improving the efficiency of transport coefficient calculations in molecular dynamics simulations.