Abstract:Thermostats/barostats that maintain temperature/pressure as constant (on average) in molecular dynamics (MD) simulations are essential for simulatingisothermal-isochoric (NVT) or isothermal-isobaric (NPT) systems. In order to simulate accurate physical properties, thermostats/barostats need not only sample a correct NVT/NPT ensemble but also minimally disturb the particles’ (Newtonian) dynamics. In other words, the NVT/NPT simulations need to yield accurate physical properties (dynamics properties included, e.g. diffusivity and viscosity) and also their fluctuations. However, few studies have studied in-depth the effects of thermostats/barostats on a comprehensively wide range of simulated properties. In this work, commonly used thermostats, e.g. Andersen, stochastic dynamics (SD), Berendsen, V-rescale and Nosé-Hoover thermostats, and barostats, e.g. Berendsen and Parrinello-Rahman barostats, with different coupling strengths were studied to see whether they can yield accurate simulated properties and fluctuations. The theoretical values of physical properties’ fluctuations were calculated via statistical mechanics to provide a comparison benchmark for MD simulations. Particularly, the accuracy of thermostats was studied in the context of NPT or non-equilibrium molecular dynamics (NEMD) simulations, which has long been overlooked. Berendsen thermostat/barostat suppresses the fluctuations of energy/volume (due to the exponential decay of deviation of the simulated temperature/pressure from the target value) and yields inaccurate simulated properties. In addition, Andersen and SD thermostats perturb the particles’ dynamics so violently (due to the random component in the algorithm) that they fail to accurately simulate dynamics properties. Overall, Nosé-Hoover/V-rescale thermostat, and Parrinello-Rahman barostat with moderate coupling strength are recommended for common NVT/NPT production simulations. NPT or NEMD simulations would require more efficient (or stronger coupling) thermostats than NVT equilibrium molecular dynamics (EMD) to maintain the simulated temperature close to the target.

On the Actual Difference between the Nosé and the Nosé–Hoover Thermostats: A Critical Review of Canonical Temperature Control by Molecular Dynamics

Canonical Temperature Control by Molecular Dynamics

Discrete Newtonian dynamics with Nosé-Hoover thermostats

Influence of Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids.

A thermostat for molecular dynamics of complex fluids

Non-ergodicity of the Nose-Hoover Thermostatted Harmonic Oscillator

Unification of Nonequilibrium Molecular Dynamics and the Mode-Resolved Phonon Boltzmann Equation for Thermal Transport Simulations

Study of Critical Dynamics in Fluid via Molecular Dynamics in Canonical Ensemble

Constant-temperature molecular-dynamics algorithms for mixed hard-core/continuous potentials

Effects of thermostats/barostats on physical properties of liquids by molecular dynamics simulations

Non-Ergodicity of Nose-Hoover chain thermostat in computationally achievable time

Analyses of the perturbative motion in the Coupled Nosé-Hoover equation for assessing the system parameter values

Nonequilibrium Temperature and Thermometry in Heat-Conducting Phi-4 Models

Fast-Forward Langevin Dynamics with Momentum Flips

A Unified Thermostat Scheme for Efficient Configurational Sampling for Classical/quantum Canonical Ensembles Via Molecular Dynamics.

From Hard Spheres and Cubes to Nonequilibrium Maps with Thirty-Some Years of Thermostatted Molecular Dynamics

Advanced Thermostats for Molecular Dynamics

A Unified Phonon Interpretation for the Non-Fourier Heat Conduction by Non-equilibrium Molecular Dynamics Simulations

Numerical and Monte Carlo Bethe Ansatz Method: 1d Heisenberg Model

Managing Temperature in Open Quantum Systems Strongly Coupled with Structured Environments