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.

Pressure control using stochastic cell rescaling

Resolving the Dynamic Properties of Entangled Linear Polymers in Non-Equilibrium Coarse Grain Simulation with a Priori Scaling Factors.

Application of Berendsen Barostat in Dissipative Particle Dynamics for Nonequilibrium Dynamic Simulation.

Constant pressure and temperature discrete-time Langevin molecular dynamics

Pressure Control in Dissipative Particle Dynamics and Its Application in Simulating Micro- and Nano-Bubbles

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

Benchmark of Schemes for Multiscale Molecular Dynamics Simulations

Stochastic boundary conditions for molecular dynamics simulations

Stochastic Dynamics with Correct Sampling for Constrained Systems

Constant Pressure Molecular Dynamics Simulation for Ionic System

Ensembles Are Required to Handle Aleatoric and Parametric Uncertainty in Molecular Dynamics Simulation

A Reversible Unwrapping Algorithm for Constant Pressure Molecular Dynamics Simulations

Density Fluctuations, Solvation Thermodynamics and Coexistence Curves in Grand Canonical Molecular Dynamics Simulations

Efficient particle control in systems with large density gradients

A molecular dynamics method for calculating molecular volume changes appropriate for biomolecular simulation

Direct Simulation and Mass Flow Stochastic Algorithms to Solve a Sintering-Coagulation Equation

Exact lattice-based stochastic cell culture simulation algorithms incorporating spontaneous and contact-dependent reactions

Feedback control for stochastic gas flow

Stochastic Approximation to MBAR and TRAM: Batchwise Free Energy Estimation

The Roles of Space and Stochasticity in Computational Simulations of Cellular Biochemistry: Quantitative Analysis and Qualitative Insights

A local resampling trick for focused molecular dynamics