Journal
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 14, Issue 10, Pages 5262-5272Publisher
AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.8b00446
Keywords
-
Funding
- Center for Gas Separations Relevant to Clean Energy Technologies, an Energy Frontier Research Center - U.S. Department of Energy, Office of Science, Basic Energy Sciences [DE-SC0001015]
- Deutsche Forschungsgemeinschaft (DFG) [SPP 1570]
- Office of Science of the U.S. Department of Energy [DE-AC02-05CH11231]
Ask authors/readers for more resources
The flying ice cube effect is a molecular dynamics simulation artifact in which the use of velocity rescaling thermostats sometimes causes violation of the equipartition theorem, affecting both structural and dynamic properties. The reason for this artifact and the conditions under which it occurs have not been fully understood. Since the flying ice cube effect was first demonstrated, a new velocity rescaling algorithm (the CSVR thermostat) has been developed and become popular without its effects on the equipartition theorem being truly known. Meanwhile, the use of simple velocity rescaling and Berendsen (weak coupling) thermostat algorithms has not abated but has actually continued to grow. Here, we have calculated the partitioning of the kinetic energy between translational, rotational, and vibrational modes in simulations of diatomic molecules to explicitly determine whether the equipartition theorem is violated under different thermostats and while rescaling velocities to different kinetic energy distributions. We have found that the underlying cause of the flying ice cube effect is a violation of balance leading to systematic redistributions of kinetic energy under simple velocity rescaling and the Berendsen thermostat. When velocities are instead rescaled to the canonical ensemble's kinetic energy distribution, as is done with the CSVR thermostat, the equipartition theorem is not violated, and we show that the CSVR thermostat satisfies detailed balance. The critical necessity for molecular dynamics practitioners to abandon the use of popular yet incorrect velocity rescaling algorithms is underscored with an example demonstrating that the main result of a highly cited study is entirely due to artifacts resulting from the study's use of the Berendsen thermostat.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available