Projects per year
Abstract / Description of output
In this article, we focus on the sampling of the configurational Gibbs-Boltzmann distribution, that is, the calculation of averages of functions of the position coordinates of a molecular $N$-body system modelled at constant temperature. We show how a formal series expansion of the invariant measure of a Langevin dynamics numerical method can be obtained in a straightforward way using the Baker-Campbell-Hausdorff lemma. We then compare Langevin dynamics integrators in terms of their invariant distributions and demonstrate a superconvergence property (4th order accuracy where only 2nd order would be expected) of one method in the high friction limit; this method, moreover, can be reduced to a simple modification of the Euler-Maruyama method for Brownian dynamics involving a non-Markovian (coloured noise) random process. In the Brownian dynamics case, 2nd order accuracy of the invariant density is achieved. All methods considered are efficient for molecular applications (requiring one force evaluation per timestep) and of a simple form. In fully resolved (long run) molecular dynamics simulations, for our favoured method, we observe up to two orders of magnitude improvement in configurational sampling accuracy for given stepsize with no evident reduction in the size of the largest usable timestep compared to common alternative methods.
Original language | English |
---|---|
Pages (from-to) | 34-56 |
Number of pages | 21 |
Journal | Applied Mathematics Research Express |
Volume | 2013 |
Issue number | 1 |
Early online date | 29 Jun 2012 |
DOIs | |
Publication status | Published - 2013 |
Keywords / Materials (for Non-textual outputs)
- molecular dynamics
- statistical mechanics
- sampling
- Langevin dynamics
- Brownian dynamics
- stochastic dynamics
- diffusion equation
Fingerprint
Dive into the research topics of 'Rational Construction of Stochastic Numerical Methods for Molecular Sampling'. Together they form a unique fingerprint.Projects
- 1 Finished
-
Science and Innovation: Numerical Algorithms and Intelligent Software for the Evolving HPC Platform
1/08/09 → 31/07/14
Project: Research