English

A new fast multiple-shooting method for computing periodic orbits in symplectic maps leveraging simultaneous Floquet vector computation to avoid large linear systems

Dynamical Systems 2026-01-05 v1 Earth and Planetary Astrophysics Chaotic Dynamics

Abstract

Given a 4D symplectic map F0F_0 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 F0F_0 is part of a perturbative family FεF_\varepsilon, one is often interested in computing those SPOs which persist for ε>0\varepsilon >0. Assuming that a persisting SPO of F0F_0 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 ε\varepsilon. 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

R2 v1 2026-07-01T08:47:33.153Z