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.

The Andersen Thermostat In Molecular Dynamics

The Andersen thermostat in molecular dynamics

Understanding Molecular Dynamics with Stochastic Processes via Real or Virtual Dynamics

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

Advanced Thermostats for Molecular Dynamics

Molecular Dynamics Simulation of Self-Diffusion Coefficient and Its Relation with Temperature Using Simple Lennard-Jones Potential

Constant pressure and temperature discrete-time Langevin molecular dynamics

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

Calculation of transport properties from molecular dynamics simulation

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

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

Thermostats and Thermostat Strategies for Molecular Dynamics Simulations of Nanofluidics.

Stochastic Dynamics with Correct Sampling for Constrained Systems

The Effects Of Lowe-Andersen Temperature Controlling Method On The Polymer Properties In Mesoscopic Simulations

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

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

A thermostat for molecular dynamics of complex fluids

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

A Simple and Accurate Algorithm for Path Integral Molecular Dynamics with the Langevin Thermostat.

Studies in molecular dynamics. XV. High temperature description of the transport coefficients

Molecular dynamics in a grand ensemble: Bergmann–Lebowitz model and adaptive resolution simulation