Some Remarks on Preconditioning Molecular Dynamics
The SMAI journal of computational mathematics, Volume 4 (2018) , pp. 57-80.

We consider a Preconditioned Overdamped Langevin algorithm that does not alter the invariant distribution (up to controllable discretisation errors) and ask whether preconditioning improves the standard model in terms of reducing the asymptotic variance and of accelerating convergence to equilibrium. We present a detailed study of the dependence of the asymptotic variance on preconditioning in some elementary toy models related to molecular simulation. Our theoretical results are supported by numerical simulations.

Published online: 2018-03-28
DOI: https://doi.org/10.5802/smai-jcm.29
Keywords: Preconditioned Overdamped Langevin algorithm, Asymptotic Variance, Central Limit Theorem, Model Hamiltonians, Lattice Model
@article{SMAI-JCM_2018__4__57_0,
author = {Houssam AlRachid and Letif Mones and Christoph Ortner},
title = {Some Remarks on Preconditioning Molecular Dynamics},
journal = {The SMAI journal of computational mathematics},
publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
volume = {4},
year = {2018},
pages = {57-80},
doi = {10.5802/smai-jcm.29},
language = {en},
url={smai-jcm.centre-mersenne.org/item/SMAI-JCM_2018__4__57_0/}
}
AlRachid, Houssam; Mones, Letif; Ortner, Christoph. Some Remarks on Preconditioning Molecular Dynamics. The SMAI journal of computational mathematics, Volume 4 (2018) , pp. 57-80. doi : 10.5802/smai-jcm.29. https://smai-jcm.centre-mersenne.org/item/SMAI-JCM_2018__4__57_0/

[1] A. Barducci; G. Bussi; M. Parrinello Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method, Phys. Rev. Lett., Volume 100(2) (2008), 020603 pages

[2] A. Beskos; G. Roberts; A. Stuart Optimal scalings for local Metropolis–Hastings chains on nonproduct targets in high dimensions, The Annals of Applied Probability, Volume 19(3) (2009), p. 863-98

[3] R. Bhattacharya A central limit theorem for diffusions with periodic coefficients, The Annals of Probability, Volume 13(2) (1985), pp. 385-396

[4] H.J. Brascamp; E.H. Lieb On extensions of the Brunn-Minkowski and Prékopa-Leindler theorems, including inequalities for log concave functions, and with an application to the diffusion equation, Journal of functional analysis, Volume 22(4) (1976), pp. 366-389

[5] A. Brünger; C. L. Brooks; M. Karplus Stochastic boundary conditions for molecular dynamics simulations of ST2 water, Chemical Physics Letters, Volume 105(5) (1984), pp. 495-500

[6] P. Cattiaux; D. Chafai; A. Guillin Central limit theorems for additive functionals of ergodic Markov diffusions processes, ALEA, Volume 9(2) (2012), pp. 337-382

[7] A. B. Duncan; T. Lelièvre; G. A. Pavliotis Variance Reduction using Nonreversible Langevin Samplers, Journal of Statistical Physics, Volume 163(3) (2016), pp. 457-491

[8] M. Girolami; B. Calderhead Riemann Manifold Langevin and Hamiltonian Monte Carlo methods, Journal of the Royal Statistical Society: Series B (Statistical Methodology), Volume 73(2) (2011), pp. 123-214

[9] P. W. Glynn; S.P. Meyn A Liapounov bound for solutions of the Poisson equation, The Annals of Probability (1996), pp. 916-931

[10] E. Hairer; C. Lubich; G. Wanner Geometric numerical integration illustrated by the Stormer-Verlet method, Acta numerica, Volume 12(12) (2003), pp. 399-450

[11] I. S. Helland Central limit theorems for martingales with discrete or continuous time, Scandinavian Journal of Statistics, Volume 9(2) (1982), pp. 79-94

[12] M. Hitrik; K. Pravda-Starov Spectra and semigroup smoothing for non-elliptic quadratic operators, Math. Ann., Volume 344(4) (2009), pp. 801-846

[13] W. Hornak; R. Abel; A. Okur; B. Strockbine; A. Roitberg; C. Simmerling Comparison of multiple Amber force fields and development of improved protein backbone parameters, Prot. Struct. Funct. Bioinf., Volume 65(3) (2006), pp. 712-725

[14] F. Hou; J. Goodman; D.W. Hogg; J. Weare; C. Schwab An affine-invariant sampler for exoplanet fitting and discovery in radial velocity data, The Astrophysical Journal, Volume 745(2) (2012)

[15] A. Iserles A First Course in the Numerical Analysis of Differential Equations, Cambridge Texts in Applied Mathematics, Cambridge University Press (2008)

[16] A. M. Stuart J. C. Mattingly; D. J. Higham Ergodicity for SDEs and approximations: locally Lipschitz vector fields and degenerate noise, Stochastic processes and their applications, Volume 101(2) (2002), pp. 185-232

[17] J. Kästner; W. Thiel Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: Umbrella integration, J. Chem. Phys., Volume 123 (2005)

[18] C. Kipnis; S.R.S. Varadhan OCentral limit theorem for additive functionals of reversible Markov processes and applications to simple exclusions, Communications in Mathematical Physics, Volume 104(1) (1986), pp. 1-19

[19] A. Laio; M. Parrinello Escaping free-energy minima, Proceedings of the National Academy of Sciences, Volume 99 (2002), pp. 12562-12566

[20] T. Lelièvre Personal communication (2017)

[21] T. Lelièvre; F. Nier; G. Pavliotis Optimal non-reversible linear drift for the convergence to equilibrium of a diffusion, Journal of Statistical Physics, Volume 152(2) (2013), pp. 237-274

[22] J. Martin; L. C. Wilcox; C. Burstedde; O. Ghattas A stochastic Newton MCMC method for large-scale statistical inverse problems with application to seismic inversion, SIAM Journal on Scientific Computing, Volume 34(3) (2012), p. A1460-A1487

[23] S. P. Meyn; R. L. Tweedie Stability of Markovian processes III: Foster-Lyapunov criteria for continuous-time processes, Advances in Applied Probability, Volume 25(3) (1993), pp. 518-548

[24] S. P. Meyn; R. L. Tweedie A survey of Foster-Lyapunov techniques for general state space Markov processes, In Proceedings of the Workshop on Stochastic Stability and Stochastic Stabilization, Metz, France. Citeseer (1993)

[25] L. Mones; N. Bernstein; G. Csanyi Exploration, Sampling, And Reconstruction of Free Energy Surfaces with Gaussian Process Regression, Journal of Chemical Theory and Computation, Volume 12(10) (2016), pp. 5100-5110

[26] L. Mones; G. Csanyi; C. Ortner in preparation (2017)

[27] D. Olson; C. Ortner Regularity and locality of point defects in multilattices, to appear in AMRX, arXiv:1608.08930 (2016)

[28] M. Ottobre; G. A. Pavliotis; K. Pravda-Starov Exponential return to equilibrium for hypoelliptic quadratic systems, Journal of Functional Analysis, Volume 262(9) (2012), pp. 4000-4039

[29] D. Packwood; J. Kermode; L. Mones; N. Bernstein; N. I. M. Gould J. Woolley; C. Ortner; G. Csanyi universal preconditioner for simulating condensed phase materials, J. Chem. Phys., Volume 144(16) (2016)

[30] N. S. Pillai; A. M. Stuart; A. H. Thiéry Optimal scaling and diffusion limits for the Langevin algorithm in high dimensions, The Annals of Applied Probability, Volume 22(6) (2017), pp. 2320-2356

[31] C. P. Robert; G. Casella Monte Carlo Statistical Methods, 2nd ed. Springer, New York (2004)

[32] G. O. Roberts; R. L. Tweedie Exponential convergence of Langevin distributions and their discrete approximations, Bernoulli, Volume 2(4) (1996), pp. 341-363

[33] V. Villani Hypocoercivity, Bull. Am. Math. Soc., Volume 46(2) (2009)

[34] T. Xifara; C. Sherlock; S. Livingstone; S. Byrne; M. Girolami Langevin diffusions and the Metropolis-adjusted Langevin algorithm, Statistics & Probability Letters, Volume 91 (2014), pp. 14-19