摘要

We investigate some unusual behaviour observed while performing molecular dynamics simulations of small molecular clusters using a constrained Langevin thermostat. Atoms appear to be thermalised to different temperatures that depend on their mass and on the total number of particles in the system. The deviation from the zeroth law of thermodynamics can be considerable for small systems of heavy and light particles. We trace this behaviour to the absence of thermal noise acting on the centre of mass of the system. This is demonstrated by solving the stochastic dynamics for the constrained thermostat and comparing the results with simulation data. By removing the constraint, the Langevin thermostat may be restored to its intended behaviour. We also investigate a Langevin thermostat constrained to have zero total force acting on its centre of mass, and find similar deficiencies.

  • 出版日期2014-11-17