Orthosymplectic integration of linear Hamiltonian systems

Abstract

The authors describe a continuous, orthogonal and symplectic factorization procedure for integrating unstable linear Hamiltonian systems. The method relies on the development of an orthogonal, symplectic change of variables to block triangular Hamiltonian form. Integration is thus carried out within the class of linear Hamiltonian systems. Use of an appropriate timestepping strategy ensures that the symplectic pairing of eigenvalues is automatically preserved. For long-term integrations, as are needed in the calculation of Lyapunov exponents, the favorable qualitative properties of such a symplectic framework can be expected to yield improved estimates. The method is illustrated and compared with other techniques in numerical experiments on the Henon-Heiles and spatially discretized Sine-Gordon equations.