Langevin dynamics with constraints and computation of free energy differences

Lelièvre, Tony; Rousset, Mathias; Stoltz, Gabriel

doi:10.1090/S0025-5718-2012-02594-4

Langevin dynamics with constraints and computation of free energy differences
HTML articles powered by AMS MathViewer

by Tony Lelièvre, Mathias Rousset and Gabriel Stoltz PDF

Math. Comp. 81 (2012), 2071-2125 Request permission

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 splitting strategy. We show how to correct the scheme so that it samples exactly the canonical measure restricted on a submanifold, using a Metropolis-Hastings correction 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. (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 discretization error is proposed, and its overdamped limit is studied. Numerical illustrations are provided for (ii) and (iii).

References

Elena Akhmatskaya and Sebastian Reich, GSHMC: an efficient method for molecular simulation, J. Comput. Phys. 227 (2008), no. 10, 4934–4954. MR 2414842, DOI 10.1016/j.jcp.2008.01.023
Luigi Ambrosio, Nicola Fusco, and Diego Pallara, Functions of bounded variation and free discontinuity problems, Oxford Mathematical Monographs, The Clarendon Press, Oxford University Press, New York, 2000. MR 1857292
V. I. Arnol′d, Mathematical methods of classical mechanics, 2nd ed., Graduate Texts in Mathematics, vol. 60, Springer-Verlag, New York, 1989. Translated from the Russian by K. Vogtmann and A. Weinstein. MR 997295, DOI 10.1007/978-1-4757-2063-1
Nawaf Bou-Rabee and Houman Owhadi, Long-run accuracy of variational integrators in the stochastic context, SIAM J. Numer. Anal. 48 (2010), no. 1, 278–297. MR 2608370, DOI 10.1137/090758842
Raphaël Chetrite and Krzysztof Gawȩdzki, Fluctuation relations for diffusion processes, Comm. Math. Phys. 282 (2008), no. 2, 469–518. MR 2421485, DOI 10.1007/s00220-008-0502-9
C. Chipot and A. Pohorille (eds.), Free energy calculations, Springer Series in Chemical Physics, vol. 86, Springer, 2007.
N. Chopin, T. Lelièvre, and G. Stoltz, Free energy methods for bayesian inference: Efficient exploration of univariate gaussian mixture posteriors, Stat. Comput., arXiv:1003.0428v4, (2011).
G. Ciccotti, R. Kapral, and E. Vanden-Eijnden, Blue Moon sampling, vectorial reaction coordinates, and unbiased constrained dynamics, Chem. Phys. Chem 6 (2005), no. 9, 1809–1814.
Giovanni Ciccotti, Tony Lelievre, and Eric Vanden-Eijnden, Projection of diffusions on submanifolds: application to mean force computation, Comm. Pure Appl. Math. 61 (2008), no. 3, 371–408. MR 2376846, DOI 10.1002/cpa.20210
Gavin E. Crooks, Nonequilibrium measurements of free energy differences for microscopically reversible Markovian systems, J. Statist. Phys. 90 (1998), no. 5-6, 1481–1487. MR 1628273, DOI 10.1023/A:1023208217925
—, Entropy production fluctuation theorem and the nonequilibrium work relation for free-energy differences, Phys. Rev. E 60 (1999), no. 3, 2721–2726.
E. Darve, Thermodynamic integration using constrained and unconstrained dynamics, Free Energy Calculations (C. Chipot and A. Pohorille, eds.), Springer, 2007, pp. 119–170.
C. Dellago, P. G. Bolhuis, and D. Chandler, On the calculation of reaction rate constants in the transition path ensemble, J. Chem. Phys. 110 (1999), no. 14, 6617–6625.
W. K. den Otter, Thermodynamic integration of the free energy along a reaction coordinate in Cartesian coordinates, J. Chem. Phys. 112 (2000), no. 17, 7283–7292.
P. A. M. Dirac, Generalized Hamiltonian dynamics, Canad. J. Math. 2 (1950), 129–148. MR 43724, DOI 10.4153/cjm-1950-012-1
S. Duane, A. D. Kennedy, B. J. Pendleton, and D. Roweth, Hybrid Monte-Carlo, Phys. Lett. B 195 (1987), no. 2, 216–222.
Marie Duflo, Random iterative models, Applications of Mathematics (New York), vol. 34, Springer-Verlag, Berlin, 1997. Translated from the 1990 French original by Stephen S. Wilson and revised by the author. MR 1485774, DOI 10.1007/978-3-662-12880-0
Weinan E and Eric Vanden-Eijnden, Metastability, conformation dynamics, and transition pathways in complex systems, Multiscale modelling and simulation, Lect. Notes Comput. Sci. Eng., vol. 39, Springer, Berlin, 2004, pp. 35–68. MR 2089952, DOI 10.1007/978-3-642-18756-8_{3}
Stewart N. Ethier and Thomas G. Kurtz, Markov processes, Wiley Series in Probability and Mathematical Statistics: Probability and Mathematical Statistics, John Wiley & Sons, Inc., New York, 1986. Characterization and convergence. MR 838085, DOI 10.1002/9780470316658
Lawrence C. Evans and Ronald F. Gariepy, Measure theory and fine properties of functions, Studies in Advanced Mathematics, CRC Press, Boca Raton, FL, 1992. MR 1158660
M. Fixman, Simulation of polymer dynamics. I. General theory, J. Chem. Phys. 69 (1978), 1527–1537.
Ernst Hairer, Christian Lubich, and Gerhard Wanner, Geometric numerical integration, 2nd ed., Springer Series in Computational Mathematics, vol. 31, Springer-Verlag, Berlin, 2006. Structure-preserving algorithms for ordinary differential equations. MR 2221614
Carsten Hartmann, An ergodic sampling scheme for constrained Hamiltonian systems with applications to molecular dynamics, J. Stat. Phys. 130 (2008), no. 4, 687–711. MR 2387561, DOI 10.1007/s10955-007-9470-2
Carsten Hartmann and Christof Schütte, A constrained hybrid Monte-Carlo algorithm and the problem of calculating the free energy in several variables, ZAMM Z. Angew. Math. Mech. 85 (2005), no. 10, 700–710. MR 2172065, DOI 10.1002/zamm.200410218
Carsten Hartmann and Christof Schütte, A geometric approach to constrained molecular dynamics and free energy, Commun. Math. Sci. 3 (2005), no. 1, 1–20. MR 2132822, DOI 10.4310/CMS.2005.v3.n1.a1
Carsten Hartmann and Christof Schütte, Comment on two distinct notions of free energy, Phys. D 228 (2007), no. 1, 59–63. MR 2334508, DOI 10.1016/j.physd.2007.02.006
Lars Hörmander, Hypoelliptic second order differential equations, Acta Math. 119 (1967), 147–171. MR 222474, DOI 10.1007/BF02392081
A. M. Horowitz, A generalized guided Monte Carlo algorithm, Phys. Lett. B 268 (1991), 247–252.
C. Jarzynski, Nonequilibrium equality for free energy differences, Phys. Rev. Lett. 78 (1997), no. 14, 2690–2693.
Wolfgang Kliemann, Recurrence and invariant measures for degenerate diffusions, Ann. Probab. 15 (1987), no. 2, 690–707. MR 885138
J. Latorre, C. Hartmann, and Ch. Schütte, Free energy computation by controlled Langevin processes, Procedia Computer Science 1 (2010), 1591–1600.
B. J. Leimkuhler and S. Reich, Simulating Hamiltonian dynamics, Cambridge Monographs on Applied and Computational Mathematics, vol. 14, Cambridge University Press, 2005.
Benedict J. Leimkuhler and Robert D. Skeel, Symplectic numerical integrators in constrained Hamiltonian systems, J. Comput. Phys. 112 (1994), no. 1, 117–125. MR 1277499, DOI 10.1006/jcph.1994.1085
Tony Lelièvre, Mathias Rousset, and Gabriel Stoltz, Computation of free energy differences through nonequilibrium stochastic dynamics: the reaction coordinate case, J. Comput. Phys. 222 (2007), no. 2, 624–643. MR 2313418, DOI 10.1016/j.jcp.2006.08.003
Tony Lelièvre, Mathias Rousset, and Gabriel Stoltz, Free energy computations, Imperial College Press, London, 2010. A mathematical perspective. MR 2681239, DOI 10.1142/9781848162488
T. Lelièvre, M. Rousset, and G. Stoltz, Langevin dynamics with constraints and computation of free energy differences, arXiv preprint 1006.4914 (2010).
P. B. Mackenzie, An improved hybrid Monte Carlo method, Phys. Lett. B 226 (1989), no. 3-4, 369–371.
Jerrold E. Marsden and Tudor S. Ratiu, Introduction to mechanics and symmetry, 2nd ed., Texts in Applied Mathematics, vol. 17, Springer-Verlag, New York, 1999. A basic exposition of classical mechanical systems. MR 1723696, DOI 10.1007/978-0-387-21792-5
G. N. Milstein and M. V. Tretyakov, Quasi-symplectic methods for Langevin-type equations, IMA J. Numer. Anal. 23 (2003), no. 4, 593–626. MR 2011342, DOI 10.1093/imanum/23.4.593
G. N. Milstein and M. V. Tretyakov, Stochastic numerics for mathematical physics, Scientific Computation, Springer-Verlag, Berlin, 2004. MR 2069903, DOI 10.1007/978-3-662-10063-9
D. D. L. Minh and A. B. Adib, Optimized free energies from bidirectional single-molecule force spectroscopy, Phys. Rev. Lett. 100 (2008), 180602.
S. Park, F. Khalili-Araghi, E. Tajkhorshid, and K. Schulten, Free energy calculation from steered molecular dynamics simulations using Jarzynski’s equality, J. Chem. Phys. 119 (2003), no. 6, 3559–3566.
D. C. Rapaport, The art of molecular dynamics simulations, Cambridge University Press, 1995.
Sebastian Reich, Smoothed Langevin dynamics of highly oscillatory systems, Phys. D 138 (2000), no. 3-4, 210–224. MR 1744627, DOI 10.1016/S0167-2789(99)00200-6
J. Schlitter and M. Klähn, A new concise expression for the free energy of a reaction coordinate, J. Chem. Phys. 118 (2003), no. 5, 2057–2060.
J. E. Straub, M. Borkovec, and B. J. Berne, Molecular-dynamics study of an isomerizing diatomic in a Lennard-Jones fluid, J. Chem. Phys. 89 (1988), no. 8, 4833–4847.
E. Vanden-Eijnden and G. Ciccotti, Second-order integrators for Langevin equations with holonomic constraints, Chem. Phys. Lett. 429 (2006), no. 1-3, 310–316.