Adaptive Runge–Kutta integration for stiff systems: Comparing Nosé and Nosé–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 constraints can exhibit either stiff or smooth dynamics. Two closely related examples, Nosé's dynam...

Full description

Saved in:
Bibliographic Details
Published inAmerican journal of physics Vol. 84; no. 10; pp. 786 - 794
Main Authors Graham Hoover, William, Clinton Sprott, Julien, Griswold Hoover, Carol
Format Journal Article
LanguageEnglish
Published Woodbury American Institute of Physics 01.10.2016
Subjects
Comparative analysis
Differential equations
Harmonic analysis
Oscillators
Online AccessGet full text

Cover

More Information
Summary: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 constraints can exhibit either stiff or smooth dynamics. Two closely related examples, Nosé's dynamics and Nosé–Hoover dynamics, are both based on Hamiltonian mechanics and generate microstates consistent with Gibbs' canonical ensemble. Nosé's dynamics is stiff and can present severe numerical difficulties. Nosé–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 Nosé dynamics for the harmonic oscillator. The solutions also illustrate the power of computer graphics to enrich numerical solutions.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0002-9505
1943-2909
DOI:10.1119/1.4959795