Skip to content

Canonical dynamics: Equilibrium phase-space distributions

Why this mattered

Hoover’s paper helped turn Nosé’s extended-system thermostat from an elegant statistical-mechanical construction into a practical molecular-dynamics method. Nosé had shown that deterministic equations of motion could generate canonical ensembles, but his formulation used a scaled time variable that made direct simulation and interpretation awkward. Hoover reformulated the dynamics without explicit time scaling, replacing it with equations in which a thermostat variable, later known through the Nosé-Hoover thermostat, acts as a dynamical friction coefficient. This made canonical sampling accessible within ordinary time evolution rather than as a post-processed or reparameterized trajectory.

The paradigm shift was that temperature control no longer had to be imposed by stochastic collisions, velocity rescaling, or weak coupling schemes external to Hamiltonian-like dynamics. Instead, equilibrium thermodynamic constraints could be represented by additional deterministic degrees of freedom whose stationary distribution was analytically tied to the canonical phase-space density. That opened a route to simulations that were simultaneously dynamical, reproducible, and ensemble-aware, especially for nonequilibrium molecular dynamics where one wants to study transport under shear, heat flow, or external driving while maintaining controlled thermodynamic conditions.

The paper also clarified a limitation that became central to later work: deterministic thermostats do not automatically guarantee ergodic sampling, as illustrated by Hoover’s analysis of the one-dimensional harmonic oscillator. That observation shaped subsequent developments such as Nosé-Hoover chains, improved barostats, and more sophisticated deterministic and stochastic thermostatting schemes. In this sense, the paper mattered both as an enabling simplification and as a diagnostic landmark: it made canonical molecular dynamics routine, while exposing the deep connection between ensemble generation, chaos, and ergodicity that later breakthroughs had to address.

Abstract

Nos\'e has modified Newtonian dynamics so as to reproduce both the canonical and the isothermal-isobaric probability densities in the phase space of an N-body system. He did this by scaling time (with s) and distance (with ${V}^{1/D}$ in D dimensions) through Lagrangian equations of motion. The dynamical equations describe the evolution of these two scaling variables and their two conjugate momenta ${p}{s}$ and ${p} act as thermodynamic friction coefficients. We find that these friction coefficients have Gaussian distributions. From the distributions the extent of small-system non-Newtonian behavior can be estimated. We illustrate the dynamical equations by considering their application to the simplest possible case, a one-dimensional classical harmonic oscillator.}$. Here we develop a slightly different set of equations, free of time scaling. We find the dynamical steady-state probability density in an extended phase space with variables x, ${p}_{x}$, V, \ensuremath{\epsilon}\ifmmode \dot{}\else .{}\fi{}, and \ensuremath{\zeta}, where the x are reduced distances and the two variables \ensuremath{\epsilon}\ifmmode \dot{}\else .{}\fi{} and \ensuremath{\zeta

Sources