Relative Periodic Orbits in the Gutzwiller-type Anisotropic Kepler Problem and n-body Problem
Published 13 Aug 2026 in math.DS, math.CA, and math.SG | (2608.12901v1)
Abstract: We study the relative periodic orbits in the Gutzwiller-type anisotropic Kepler problem, which is a generalized model derived from the classical Gutzwiller anisotropic Kepler problem. By reducing the rotational symmetry of the z-axis, we obtain a reduced system with two degrees of freedom and its energy surface forms a compact and regular three-sphere in a certain parameter range. Combining the estimation of the Seifert rotation number of the planar Kepler orbit and the CHHL formula introduced in \cite{CHHL23}, we prove that this system admits infinitely many periodic orbits on every compact regular energy surface. This model can be applied to the (1+2n)-body problem to find infinitely many relative periodic orbits with hip-hop symmetry as well as relative periodic orbits in the n-pyramidal problem.
The paper proves infinitely many periodic orbits on every compact regular energy surface when the effective anisotropy satisfies β=D/C−1≥0, using contact-volume estimates, rotation numbers, and the CHHL two-or-infinity theorem.
The analysis shows that the planar Kepler orbit’s rotation number strictly exceeds the volume threshold except in the integrable spatial Kepler case, where equality occurs.
The result guarantees infinitely many hip-hop relative periodic orbits for all central mass ratios and extends the n-pyramidal conclusion to all mass ratios for 2≤n≤472, while leaving the D<C regime open.
Overview and main result
This paper by Hu, Ou, Qiao, and Yang (2608.12901) studies the reduced Gutzwiller-type anisotropic Kepler problem, a two-degree-of-freedom Hamiltonian system of the form
with A0≥0 and Ai,Bi>0. The classical Gutzwiller anisotropic Kepler problem corresponds to A0=0, n=1. The main theorem states that, under the condition i=1∑nAiBi≥i=0∑nAi (equivalently β=D/C−1≥0 with D=∑AiBi, C=∑Ai), every compact regular energy surface admits infinitely many periodic orbits. The proof combines a contact-volume estimate, an estimate of the Seifert rotation number of the planar Kepler orbit via Maslov-type index theory, and the CHHL formula from the resolution of the two-or-infinity conjecture [CHHL23, CHHL26].
The strategy is a comparison argument: if the CHHL formula held with only two simple Reeb orbits, one would have vol(M,λ)=Tp2/ρ^p for the planar Kepler orbit A0≥00. The paper shows instead that A0≥01 strictly (with equality only in the integrable spatial Kepler case), forcing infinitely many Reeb orbits.
Topology of energy surfaces and the reduction framework
A key structural parameter is A0≥02, where A0≥03 is the energy level and A0≥04 the conserved angular momentum about the A0≥05-axis. The paper classifies the topology completely: when A0≥06, the energy surface A0≥07 is compact and homeomorphic to A0≥08; at A0≥09 it degenerates to a single point (the non-degenerate global minimum); below this value it is empty; and for Ai,Bi>00 it is unbounded in the Ai,Bi>01-direction. In the compact range, Ai,Bi>02 is a contact-type hypersurface in Ai,Bi>03, so the Hamiltonian flow coincides with a Reeb flow up to reparametrization.
The planar Kepler orbit Ai,Bi>04 is explicit, with eccentricity Ai,Bi>05, and serves as the distinguished Reeb orbit whose rotation number is compared against the contact volume.
Contact volume estimates
The volume analysis reduces the general system to a comparison anisotropic Kepler problem with a single effective anisotropy Ai,Bi>06. Using strict convexity of Ai,Bi>07 and Jensen's inequality, the authors show
Ai,Bi>08
with equality if and only if Ai,Bi>09 and all A0=00 coincide — i.e., precisely when the system is the reduced spatial Kepler problem. This volume bound is sharp in exactly the integrable case, which is what makes the subsequent strict rotation-number inequality decisive.
Rotation number via the Hill stability equation
The linearized flow along A0=01 decouples into two subsystems. The transverse one, after a time-dependent symplectic change of variables, becomes the Hill stability equation
A0=02
a special case of Ince's equation that also governs the stability of homographic solutions in the A0=03-body problem. Its Seifert rotation number A0=04 is analyzed through the A0=05-degenerate curves A0=06 of the associated self-adjoint Fredholm operator.
For A0=07, the paper extends a previously known inequality (valid for A0=08) to all A0=09 by proving, via a blow-up technique near collision singularities and Sturm comparison applied to the coefficient function n=10, that
n=11
where the explicit threshold n=12 satisfies n=13 on n=14. Combining this with the volume estimate yields n=15 whenever n=16, and also in the degenerate case n=17 with genuine anisotropy (n=18 or unequal n=19), where i=1∑nAiBi≥i=0∑nAi0 but the volume strictly exceeds i=1∑nAiBi≥i=0∑nAi1. By the corollary of the CHHL formula, each compact regular energy surface carries infinitely many periodic orbits. Equality holds only for the integrable Kepler case, consistent with the earlier results of Hu–Ou–Qiao covering i=1∑nAiBi≥i=0∑nAi2; together with a conformal symplectic rescaling handling i=1∑nAiBi≥i=0∑nAi3, infinitely many periodic orbits are now established for the classical anisotropic Kepler problem for alli=1∑nAiBi≥i=0∑nAi4.
Applications to the hip-hop i=1∑nAiBi≥i=0∑nAi5-body problem
The hip-hop symmetric invariant subsystem of the i=1∑nAiBi≥i=0∑nAi6-body problem — i=1∑nAiBi≥i=0∑nAi7 equal masses forming two regular i=1∑nAiBi≥i=0∑nAi8-gons reflected across planes perpendicular to the i=1∑nAiBi≥i=0∑nAi9-axis, plus a central mass of ratio β=D/C−1≥00 fixed at the center of mass — reduces after fixing angular momentum to a Hamiltonian of the Gutzwiller type with β=D/C−1≥01, β=D/C−1≥02, β=D/C−1≥03, β=D/C−1≥04, and β=D/C−1≥05, β=D/C−1≥06. The planar Kepler orbit here is the elliptic relative equilibrium (ERE) of the regular β=D/C−1≥07-gon central configuration.
Verifying the hypothesis β=D/C−1≥08 amounts to showing β=D/C−1≥09, where D=∑AiBi0 and D=∑AiBi1. The authors prove this for all D=∑AiBi2 and D=∑AiBi3 using the lower bound D=∑AiBi4, which exceeds D=∑AiBi5 already at D=∑AiBi6 and diverges as D=∑AiBi7. Consequently, every compact regular energy surface of the hip-hop system carries infinitely many relative periodic orbits with hip-hop symmetry under Newtonian gravity, for all central mass ratios. This strengthens Terracini–Venturelli's variational existence result for weak-force potentials: their construction did not control which energy surface the orbits lie on, whereas the present result gives orbits on every compact regular energy surface. The paper notes, however, that the total rotation angle about the D=∑AiBi8-axis over one period cannot be determined for any given orbit — a limitation inherited from the method.
Applications to the D=∑AiBi9-pyramidal problem
For the C=∑Ai0-pyramidal problem (C=∑Ai1 masses, C=∑Ai2 identical masses in a regular C=∑Ai3-gon symmetric about a fixed axis along which the first body moves), the reduced Hamiltonian fits the framework with C=∑Ai4, C=∑Ai5, C=∑Ai6, C=∑Ai7, where C=∑Ai8 is the mass ratio. Here the hypothesis C=∑Ai9 becomes vol(M,λ)=Tp2/ρ^p0, i.e. vol(M,λ)=Tp2/ρ^p1. Since vol(M,λ)=Tp2/ρ^p2 grows like vol(M,λ)=Tp2/ρ^p3, this holds for vol(M,λ)=Tp2/ρ^p4 — a concrete numerical threshold. Within this range, every compact regular energy surface admits infinitely many periodic orbits for all mass ratios vol(M,λ)=Tp2/ρ^p5, extending prior results that required either sufficiently large mass ratio or sufficiently small vol(M,λ)=Tp2/ρ^p6. For vol(M,λ)=Tp2/ρ^p7 the hypothesis fails and the question remains open.
Limitations and open questions
Several restrictions are intrinsic to the approach. The main theorem requires vol(M,λ)=Tp2/ρ^p8; systems with vol(M,λ)=Tp2/ρ^p9 are not covered, and the corresponding rotation-number comparison is not established. The A0≥000-pyramidal application is limited to A0≥001 by the same condition. For the hip-hop problem, although infinitely many relative periodic orbits exist on each compact energy surface, the winding number about the symmetry axis remains undetermined. Finally, the result concerns only compact regular energy surfaces in the window A0≥002; unbounded energy surfaces (A0≥003) fall outside the scope, though analogous infinite-orbit statements have been obtained for non-compact surfaces in related problems such as the spatial isosceles three-body problem.
Conclusion
The paper establishes that the Gutzwiller-type anisotropic Kepler problem admits infinitely many periodic orbits on every compact regular energy surface whenever the effective anisotropy parameter satisfies A0≥004, with equality in the governing volume–rotation-number inequality characterizing exactly the integrable spatial Kepler problem. The proof synthesizes Jensen-based contact volume comparison, Maslov-index computation of the Seifert rotation number through the Hill stability equation, and the CHHL two-or-infinity dichotomy. As applications, it yields infinitely many hip-hop-symmetric relative periodic orbits in the A0≥005-body problem for all A0≥006 and all central mass ratios, and infinitely many periodic orbits in the A0≥007-pyramidal problem for all mass ratios when A0≥008.