ENERGY FLUCTUATIONS INDUCED BY THE NOSE THERMOSTAT

The Nose thermostat is a method for simulating a canonical ensemble in molecular-dynamics calculations. It is generally agreed that for large values of Q, the parameter determining the response of the thermostat to the system, the approach to equilibrium is slow and for intermediate values it is rapid. However, for small Q Nose found the thermostat oscillates rapidly and does not equilibrate with the system; we found it equilibrates but yields incorrect values for energy fluctuations, whereas Cho and Joan-nopoulos found it equilibrates very slowly, eventually yielding correct energy fluctuations. Here we demonstrate how these differences arise.