HAL will be down for maintenance from Friday, June 10 at 4pm through Monday, June 13 at 9am. More information
Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

Langevin dynamics with constraints and computation of free energy differences

Abstract : In this paper, we consider Langevin processes with mechanical constraints. The latter are a fundamental tool in molecular dynamics simulation for sampling purposes and for the computation of free energy differences. The results of this paper can be divided into three parts. (i) We propose a simple discretization of the constrained Langevin process based on a standard splitting strategy. We show how to correct the scheme so that it samples {\em exactly} the canonical measure restricted on a submanifold, using a Metropolis rule in the spirit of the Generalized Hybrid Monte Carlo (GHMC) algorithm. Moreover, we obtain, in some limiting regime, a consistent discretization of the overdamped Langevin (Brownian) dynamics on a submanifold, also sampling exactly the correct canonical measure with constraints. The corresponding numerical methods can be used to sample (without any bias) a probability measure supported by a submanifold. (ii) For free energy computation using thermodynamic integration, we rigorously prove that the longtime average of the Lagrange multipliers of the constrained Langevin dynamics yields the gradient of a rigid version of the free energy associated with the constraints. A second order time discretization using the Lagrange multipliers is proposed. (iii) The Jarzynski-Crooks fluctuation relation is proved for Langevin processes with mechanical constraints evolving in time. An original numerical discretization without time-step error is proposed. Numerical illustrations are provided for (ii) and (iii).
Document type :
Preprints, Working Papers, ...
Complete list of metadata

https://hal.archives-ouvertes.fr/hal-00495517
Contributor : Gabriel Stoltz Connect in order to contact the contributor
Submitted on : Monday, June 28, 2010 - 9:36:06 AM
Last modification on : Wednesday, March 23, 2022 - 3:50:58 PM

Links full text

Identifiers

  • HAL Id : hal-00495517, version 1
  • ARXIV : 1006.4914

Collections

Citation

Tony Lelièvre, Mathias Rousset, Gabriel Stoltz. Langevin dynamics with constraints and computation of free energy differences. 2010. ⟨hal-00495517⟩

Share

Metrics

Record views

173