期刊
JOURNAL OF MOLECULAR LIQUIDS
卷 365, 期 -, 页码 -出版社
ELSEVIER
DOI: 10.1016/j.molliq.2022.120116
关键词
Thermostat; Barostat; Molecular dynamics simulation; Physical property; Isothermal-isochoric(NVT) ensemble; Isothermal-isobaric(NPT) ensemble
资金
- National Natural Science Foundation of China [22078104]
- Science and Technology Key Project of Guangdong Province [2020B010188002]
- Guangdong Provincial Key Lab of Green Chemical Product Technology [GC202111]
The choice of thermostats/barostats has a significant impact on the accuracy of simulated properties in molecular dynamics simulations. Nose-Hoover/V-rescale thermostat and Parrinello-Rahman barostat with moderate coupling strength are recommended for common NVT/NPT production simulations.
Thermostats/barostats that maintain temperature/pressure as constant (on average) in molecular dynamics (MD) simulations are essential for simulating isothermal-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 Nose-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 sup-presses 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, Nose-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 simu-lated temperature close to the target. (c) 2022 Elsevier B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据