Systematic errors in diffusion coefficients from long-time molecular dynamics simulations at constant pressure
arXiv:2003.09205 · doi:10.1063/5.0008316
Abstract
In molecular dynamics simulations under periodic boundary conditions, particle positions are typically wrapped into a reference box. For diffusion coefficient calculations using the Einstein relation, the particle positions need to be unwrapped. Here, we show that a widely used heuristic unwrapping scheme is not suitable for long simulations at constant pressure. Improper accounting for box-volume fluctuations creates, at long times, unphysical trajectories and, in turn, grossly exaggerated diffusion coefficients. We propose an alternative unwrapping scheme that resolves this issue. At each time step, we add the minimal displacement vector according to periodic boundary conditions for the instantaneous box geometry. Here and in a companion paper [J. Chem. Phys. XXX, YYYYY (2020)], we apply the new unwrapping scheme to extensive molecular dynamics and Brownian dynamics simulation data. We provide practitioners with a formula to assess if and by how much earlier results might have been affected by the widely used heuristic unwrapping scheme.
arXiv comment6 pages, 5 figures. The following article has been accepted for publication at The Journal of Chemical Physics
References in corpus (2)
Cited by in corpus (8)
- Optimal estimates of diffusion coefficients from molecular dynamics simulations
- Small ionic radii limit time step in Martini 3 molecular dynamics simulations
- Bayesian Estimation of the Hydroxyl Radical Diffusion Coefficient at Low Temperature and High Pressure from Atomistic Molecular Dynamics
- Hydrodynamic Interactions in Ion Transport -- Theory and Simulation
- An OrthoBoXY-Method for Various Alternative Box Geometries
- Flow Matching for Accelerated Simulation of Atomic Transport in Crystalline Materials
- H-bonds in Crambin: Coherence in an alpha helix
- A Reversible Unwrapping Algorithm for Constant Pressure Molecular Dynamics Simulations