Search papers, labs, and topics across Lattice.
The paper introduces a thermodynamic extrapolation method to remove discretization errors in molecular dynamics simulations arising from finite time steps. By modeling the time-step dependence of thermodynamic observables (energy, volume, temperature) as linear functions of Δt², the method allows for extrapolation to the zero time step limit, effectively recovering Boltzmann-consistent statistics. This approach enables the use of larger time steps without compromising thermodynamic accuracy, particularly beneficial for enhanced sampling methods.
MD simulations can now achieve Boltzmann-consistent statistics independent of the chosen time step, unlocking faster simulations without sacrificing thermodynamic accuracy.
The integration time step is a critical determinant of performance in molecular dynamics simulations, governing the trade-off between speed and fidelity. Although 2 fs remains the standard in atomistic biomolecular simulations, the push for performance has popularized a 4 fs time step with hydrogen mass repartitioning, often combined with multiple time stepping or mass rescaling. However, it is often unclear whether a chosen protocol is overly aggressive, as the apparent numerical stability of a trajectory can mask underlying thermodynamic inaccuracies. Increasing the time step will exacerbate systematic discretization errors, inherent to all numerical integration algorithms. In the widely used Verlet family of integrators, these errors manifest as $\mathcal{O}(Δt^2)$ deviations in thermodynamic observables such as potential energy and volume, and for common Langevin splitting schemes, even temperature. We demonstrate that these deviations follow a simple, linear thermodynamic model, allowing for their rigorous removal by extrapolation to the zero time step limit. In turn, the time-step dependence provides us with estimates of the system heat capacity, compressibility, and thermal expansion coefficient. This framework allows us to construct consistent probability distributions of energy and volume across thermodynamic states, effectively recovering Boltzmann-consistent statistics at a target condition independent of time step. These considerations are particularly important for enhanced sampling methods such as replica exchange and umbrella sampling, which rely on rigorous Boltzmann sampling and require accurate energies and temperatures for valid replica exchange probabilities and statistical reweighting.
