Robust and efficient configurational molecular sampling via Langevin dynamics

B. Leimkuhler, C. Matthews

Research output: Contribution to journalArticlepeer-review

Abstract

A wide variety of numerical methods are evaluated and compared for solving the stochastic differential equations encountered in molecular dynamics. The methods are based on the application of deterministic impulses, drifts, and Brownian motions in some combination. The Baker-Campbell-Hausdorff expansion is used to study sampling accuracy following recent work by the authors, which allows determination of the stepsize-dependent bias in configurational averaging. For harmonic oscillators, configurational averaging is exact for certain schemes, which may result in improved performance in the modelling of biomolecules where bond stretches play a prominent role. For general systems, an optimal method can be identified that has very low bias compared to alternatives. In simulations of the alanine dipeptide reported here (both solvated and unsolvated), higher accuracy is obtained without loss of computational efficiency, while allowing large timestep, and with no impairment of the conformational exploration rate (the effective diffusion rate observed in simulation). The optimal scheme is a uniformly better performing algorithm for molecular sampling, with overall efficiency improvements of 25 or more in practical timestep size achievable in vacuum, and with reductions in the error of configurational averages of a factor of ten or more attainable in solvated simulations at large timestep.
Original languageEnglish
JournalJournal of Chemical Physics
Volume138
Issue number17
DOIs
Publication statusPublished - 7 May 2013

Fingerprint

Dive into the research topics of 'Robust and efficient configurational molecular sampling via Langevin dynamics'. Together they form a unique fingerprint.

Cite this