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.

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

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

A thermostat for molecular dynamics of complex fluids

Uncovering novel phase transitions in dense dry polar active fluids using a lattice Boltzmann method

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

Canonical Temperature Control by Molecular Dynamics

Advanced Thermostats for Molecular Dynamics

A Novel Monte Carlo Approach to the Dynamics of Fluids --- Single Particle Diffusion, Correlation Functions and Phase Ordering of Binary Fluids

Universal two-component dynamics in supercritical fluids

Microcanonical Statistical Mechanics of a Nonideal Fluid: A Pedagogical Approach to Obtain the Combined First and Second Law of Thermodynamics

Achievements and Prospects of Molecular Dynamics Simulations in Thermofluid Sciences

Monte Carlo Simulation of Thermodynamic Properties for Two-Dimensional Lennard-Jones Fluids

Particle fluid interactivity deteriorates buoyancy driven thermal transport in nanosuspensions : A multi component lattice Boltzmann approach

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

Molecular dynamics study of nanoconfined TIP4P/2005 water: how confinement and temperature affect diffusion and viscosity

Comparison of Dissipative Particle Dynamics and Langevin thermostats for out-of-equilibrium simulations of polymeric systems

Mesoscopic Characterization of Bubble Dynamics in Flow Boilling following A Pseudopotential-based Approach

Asymptotic methods for confined fluids

Thermal convection of viscoelastic fluids in concentric rotating cylinders: Elastic turbulence and kinetic energy budget analysis

Generalized Langevin dynamics of a nanoparticle using a finite element approach: thermostating with correlated noise

Discrete Laplacian Thermostat for Spin Systems with Conserved Dynamics