Mobility estimation for Langevin dynamics using control variates

06/20/2022
by   G. A. Pavliotis, et al.
0

The scaling of the mobility of two-dimensional Langevin dynamics in a periodic potential as the friction vanishes is not well understood for non-separable 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 previous numerical evidence on the scaling of the mobility in the low friction regime for a simple non-separable potential.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/29/2018

Stochastic Gradient Hamiltonian Monte Carlo with Variance Reduction for Bayesian Inference

Gradient-based Monte Carlo sampling algorithms, like Langevin dynamics a...
research
08/24/2017

Randomized Dimension Reduction for Monte Carlo Simulations

We present a new unbiased algorithm that estimates the expected value of...
research
09/21/2022

Particle-based, rapid incremental smoother meets particle Gibbs

The particle-based, rapid incremental smoother (PARIS) is a sequential M...
research
12/14/2019

Control variates and Rao-Blackwellization for deterministic sweep Markov chains

We study control variate methods for Markov chain Monte Carlo (MCMC) in ...
research
05/04/2022

Asymptotic analysis of diabatic surface hopping algorithm in the adiabatic and non-adiabatic limits

Surface hopping algorithms, as an important class of quantum dynamics si...
research
06/14/2022

Unbiased Estimation using the Underdamped Langevin Dynamics

In this work we consider the unbiased estimation of expectations w.r.t. ...
research
01/02/2023

State and parameter learning with PaRIS particle Gibbs

Non-linear state-space models, also known as general hidden Markov model...

Please sign up or login with your details

Forgot password? Click here to reset