Theoretical and numerical comparison of some sampling methods for molecular dynamics
Open Access
- 1 March 2007
- journal article
- research article
- Published by EDP Sciences in ESAIM: Mathematical Modelling and Numerical Analysis
- Vol. 41 (2), 351-389
- https://doi.org/10.1051/m2an:2007014
Abstract
The purpose of the present article is to compare different phase-space sampling methods, such as purely stochastic methods (Rejection method, Metropolized independence sampler, Importance Sampling), stochastically perturbed Molecular Dynamics methods (Hybrid Monte Carlo, Langevin Dynamics, Biased Random Walk), and purely deterministic methods (Nosé-Hoover chains, Nosé-Poincaré and Recursive Multiple Thermostats (RMT) methods). After recalling some theoretical convergence properties for the various methods, we provide some new convergence results for the Hybrid Monte Carlo scheme, requiring weaker (and easier to check) conditions than previously known conditions. We then turn to the numerical efficiency of the sampling schemes for a benchmark model of linear alkane molecules. In particular, the numerical distributions that are generated are compared in a systematic way, on the basis of some quantitative convergence indicators.Keywords
This publication has 59 references indexed in Scilit:
- Quasi-symplectic methods for Langevin-type equationsIMA Journal of Numerical Analysis, 2003
- Sur quelques algorithmes récursifs pour les probabilités numériquesESAIM: Probability and Statistics, 2001
- Understanding Modern Molecular Dynamics: Techniques and ApplicationsThe Journal of Physical Chemistry B, 1999
- Explicit reversible integrators for extended systems dynamicsMolecular Physics, 1996
- Nosé–Hoover chains: The canonical ensemble via continuous dynamicsThe Journal of Chemical Physics, 1992
- A generalized guided Monte Carlo algorithmPhysics Letters B, 1991
- Second-order discretization schemes of stochastic differential systems for the computation of the invariant lawStochastics and Stochastic Reports, 1990
- An improved hybrid Monte Carlo methodPhysics Letters B, 1989
- Hybrid Monte CarloPhysics Letters B, 1987
- Smooth Transition Densities for One-Dimensional DiffusionsBulletin of the London Mathematical Society, 1985