Harmonic oscillators in the Nosé-Hoover environment

We study the dynamics of an ensemble of noninteracting harmonic oscillators in a nonlinear dissipative environment described by the Nosé-Hoover model, and find the histogram for energy regions of phase space against visiting time by employing numerical simulation. The results agree with the analysis of the Nosé-Hoover equations effected with the method of averaging for small values of the dissipative parameter α of the thermostat. We find oscillations at frequencies proportional to α/m,m being the characteristic mass of the particle, about the stationary state corresponding to equilibrium, for sufficiently small α . In this region of α the histogram does not correspond to Gibbs’ canonical distribution. For larger values of α the motion becomes irregular. The phenomena could have an important bearing upon simulating molecular dynamics in the Nosé-Hoover thermostat.