Adaptive Runge–Kutta integration for stiff systems: Comparing Nose and Nose–Hoover dynamics for the harmonic oscillator

We describe the application of adaptive (variable time step) integrators to stiff differential equations encountered in many applications. Linear harmonic oscillators subject to nonlinear thermal can exhibit either stiff or smooth dynamics. Two closely related examples, Nose’s dynamics and Nose–Hoover dynamics, are both based on Hamiltonian mechanics and generate microstates consistent with Gibbs’ canonical ensemble. Nose’s dynamics is stiff and can present severe numerical difficulties. Nose–Hoover dynamics, although it follows exactly the same trajectory, is smooth and relatively trouble-free. We emphasize the power of adaptive integrators to resolve stiff problems such as the Nose dynamics for the harmonic oscillator. The solutions also illustrate the power of computer graphics to enrich numerical solutions.