MOBILITY ESTIMATION FOR LANGEVIN DYNAMICS USING CONTROL VARIATES

The scaling of the mobility of two-dimensional Langevin dynamics in a periodic potential as the friction vanishes is not well understood for nonseparable potentials. Theoretical results are lacking, and numerical calculation of the mobility in the underdamped regime is challenging because the computational cost of standard Monte Carlo methods is inversely proportional to the friction coefficient, while deterministic methods are ill-conditioned. In this work, we propose a new variance-reduction method based on control variates for efficiently estimating the mobility of Langevin-type dynamics. We provide bounds on the bias and variance of the proposed estimator and illustrate its efficacy through numerical experiments, first in simple one-dimensional settings and then for two-dimensional Langevin dynamics. Our results corroborate prior numerical evidence that the mobility scales as γ^-σ, with 0 < σ ≤ 1, in the low friction regime for a simple nonseparable potential.