Assessing numerical methods for molecular and particle simulation

Xiaocheng Shang, Martin Kröger, Benedict Leimkuhler

Research output: Contribution to journalArticlepeer-review

10 Citations (Scopus)
139 Downloads (Pure)

Abstract

We discuss the design of state-of-the-art numerical methods for molecular dynamics, focusing on the demands of soft matter simulation, where the purposes include sampling and dynamics calculations both in and out of equilibrium. We discuss the characteristics of different algorithms, including their essential conservation properties, the convergence of averages, and the accuracy of numerical discretizations. Formulations of the equations of motion which are suited to both equilibrium and nonequilibrium simulation include Langevin dynamics, dissipative particle dynamics (DPD), and the more recently proposed “pairwise adaptive Langevin” (PAdL) method, which, like DPD but unlike Langevin dynamics, conserves momentum and better matches the relaxation rate of orientational degrees of freedom. PAdL is easy to code and suitable for a variety of problems in nonequilibrium soft matter modeling; our simulations of polymer melts indicate that this method can also provide dramatic improvements in computational efficiency. Moreover we show that PAdL gives excellent control of the relaxation rate to equilibrium. In the nonequilibrium setting, we further demonstrate that while PAdL allows the recovery of accurate shear viscosities at higher shear rates than are possible using the DPD method at identical timestep, it also outperforms Langevin dynamics in terms of stability and accuracy at higher shear rates.
Original languageEnglish
Pages (from-to)8565-8578
JournalSoft Matter
Volume13
Issue number45
DOIs
Publication statusPublished - 25 Oct 2017

Fingerprint

Dive into the research topics of 'Assessing numerical methods for molecular and particle simulation'. Together they form a unique fingerprint.

Cite this