A new fast multiple-shooting method for computing periodic orbits in symplectic maps leveraging simultaneous Floquet vector computation to avoid large linear systems
Abstract
Given a 4D symplectic map that has a normally hyperbolic invariant cylinder foliated by invariant tori, those with rational rotation numbers are themselves foliated by subharmonic periodic orbits (SPOs). If is part of a perturbative family , one is often interested in computing those SPOs which persist for . Assuming that a persisting SPO of has been identified, in this paper, we develop a quasi-Newton method which solves for the SPO simultaneously with its Floquet vectors and multipliers. This in turn enables continuation by the perturbation parameter . The resulting SPO and Floquet vectors are then used to compute Taylor parameterizations of the SPO's weak stable and unstable manifolds, if they exist. Our quasi-Newton method is based on an adaptation of the parameterization method for invariant tori, with this paper being the first-ever to apply such a framework to directly compute periodic orbit points themselves. The new algorithm improves on efficiency compared to prior multi-shooting methods for SPOs, and notably applies to the case of stroboscopic maps of 2.5 DOF Hamiltonian flows resulting from periodic perturbations of 2 DOF systems. The tools have been successfully used for studies of resonant orbits in perturbed real-life celestial systems, the results of which are summarized as a demonstration of the methods' utility.
Keywords
Cite
@article{arxiv.2601.00149,
title = {A new fast multiple-shooting method for computing periodic orbits in symplectic maps leveraging simultaneous Floquet vector computation to avoid large linear systems},
author = {Bhanu Kumar},
journal= {arXiv preprint arXiv:2601.00149},
year = {2026}
}
Comments
32 pages, 7 figures