Phase space volume preserving dynamics for non-Hamiltonian systems

TL;DR

This paper introduces a volume-preserving linearized dynamics framework for non-Hamiltonian systems, enabling invariant phase space measures and improved analysis of chaotic and dissipative systems.

Contribution

It proposes a novel volume-preserving operator derived from the stability matrix, extending classical density matrix theory to non-Hamiltonian dynamics.

Findings

Invariant phase space volume achieved for non-Hamiltonian systems.

Efficient computation of Lyapunov exponents and entropy flow.

Applicable to chaotic, dissipative, and driven systems without re-orthogonalization.

Abstract

Infinitesimal volumes stretch and contract as they coevolve with classical phase space trajectories according to linearized dynamics. Unless these tangent-space dynamics are modified, chaotic evolution causes the volume spanned by evolving tangent vectors to collapse. However, this collapse is unphysical and due to their exponential alignment along the most expanding direction, independent of the compressibility of the phase-space volume. Here, we propose an alternative linearized dynamics and rectify the generalized Liouville equation to preserve phase space volume, even for non-Hamiltonian systems. Within a classical density matrix theory, we define the time evolution operator from the anti-symmetric part of the stability matrix so that phase space volume is time-invariant. The operator generates orthogonal transformations without distorting volume elements, providing an invariant measure for dissipative dynamics and a evolution equation for the density matrix akin to the quantum mechanical Liouville-von Neumann equation. The compressibility of volume elements is determined by a non-orthogonal operator made from the symmetric part of the stability matrix. We analyze complete sets of basis vectors for the tangent space dynamics of chaotic systems, which may be dissipative, transient or driven, without re-orthogonalization of tangent vectors. The linear harmonic oscillator, the Lorenz-Fetter model, and the H\'enon-Heiles system demonstrate the computation of the instantaneous Lyapunov exponent spectrum and the local Gibbs entropy flow rate using these bases and show that it is numerically convenient.