skip to: onlinetools | mainnavigation | content | footer

Stephen D. Bond

Citation:

Stephen D. Bond and Benedict J. Leimkuhler, Molecular Dynamics and the Accuracy of Numerically Computed Averages, Acta Numerica 16 (2007) 1-65.

Abstract:

Molecular dynamics is discussed from a mathematical perspective. The recent history of method development is briefly surveyed with an emphasis on the use of geometric integration as a guiding principle. The recovery of statistical mechanical averages from molecular dynamics is then introduced, and the use of backward error analysis as a technique for analysing the accuracy of numerical averages is described. This article gives the first rigorous estimates for the error in statistical averages computed from molecular dynamics simulation based on backward error analysis. It is shown that molecular dynamics introduces an appreciable bias at stepsizes which are below the stability threshold. Simulations performed in such a regime can be corrected by use of a stepsize-dependent reweighting factor. Numerical experiments illustrate the efficacy of this approach. In the final section, several open problems in dynamics-based molecular sampling are considered.

Full Text:

DOI: 10.1017/S0962492906280012

Errata:

Errata (68K, PDF)

Preprint:

Molecular Dynamics and the Accuracy of Numerically Computed Averages (1.1M, PDF)

Bibtex:
@article{BoLe2007,
  author  = {Stephen D. Bond and Benedict J. Leimkuhler},
  title   = {Molecular Dynamics and the Accuracy of Numerically Computed
             Averages},
  journal = {Acta Numerica},
  volume  = 16,
  year    = 2007,
  pages   = {1--65},
  doi     = {10.1017/S0962492906280012}
}

Stephen Bond