pith. sign in

arxiv: 2606.08485 · v1 · pith:VDWATIVJnew · submitted 2026-06-07 · 🧮 math.DS

Linking Averaged and Unaveraged Three-Body Dynamics Near Smaller Primaries: Symmetric Periodic Orbits

Beom Park , Kathleen C. Howell This is my paper

Pith reviewed 2026-06-27 18:10 UTC · model grok-4.3

classification 🧮 math.DS
keywords three-body problemperiodic orbitsaveraged dynamicsbifurcation analysisrestricted three-body problemfrequency frameworksymmetric orbitsinvariant tori
0
0 comments X

The pith

Averaged equilibria in three-body dynamics correspond to symmetric periodic orbits in the unaveraged models.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper aims to connect simplified averaged models of three-body motion near a small primary with the more accurate unaveraged versions such as the Hill Restricted Three-Body Problem and Circular Restricted Three-Body Problem. It establishes this link by showing how equilibria from the averaged system map to families of symmetric periodic orbits in the full models. A unified frequency framework is used to characterize how invariant tori correspond across the different dynamical descriptions. Resonance ratio parity supplies an initialization scheme that predicts the number and symmetry types of admissible solutions in advance. Bifurcation and frequency analysis then trace the global structure of the resulting orbit families.

Core claim

Averaged equilibria are explicitly linked to symmetric periodic orbits in the unaveraged three-body systems through a unified frequency framework that characterizes the mapping of invariant tori, with an initialization scheme based on resonance ratio parity allowing a priori prediction of solution multiplicity and symmetry types.

What carries the argument

The unified frequency framework for mapping invariant tori across averaged and unaveraged dynamical models, together with the resonance-ratio-parity initialization scheme for apse configurations.

If this is right

  • Symmetric periodic orbit families can be traced through their global evolution using bifurcation and frequency analysis.
  • Archetypical bifurcation diagrams provide an atlas of the symmetric periodic orbit web within the HR3BP and CR3BP.
  • The topological origins of complex periodic orbit families become visible from the averaged starting points.
  • The resulting atlas supplies a practical reference for selecting initial conditions in multi-body trajectory design.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same resonance-parity rule could be tested for orbit families that are not required to be symmetric.
  • Extending the frequency mapping to include small additional perturbations would check robustness outside the exact restricted problems.
  • The initialization scheme might reduce the computational cost of generating starting guesses for long-period families in higher-fidelity ephemeris models.

Load-bearing premise

The parity of the resonance ratio enables an a priori initialization scheme to identify admissible apse configurations and predict solution multiplicity and symmetry types.

What would settle it

A numerical continuation search in the CR3BP for a chosen resonance ratio that yields a different count or symmetry set of symmetric periodic orbits than the multiplicity predicted by the parity-based initialization scheme.

Figures

Figures reproduced from arXiv: 2606.08485 by Beom Park, Kathleen C. Howell.

Figure 1
Figure 1. Figure 1: Geometry of the basis vectors for the inertial and rotating frames. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: DADM solution space within the reduced domain, [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: η = νs/νm for the DADM constructed for µ = 0 and µ ≈ 0.01215 (Earth-Moon system). the unaveraged dynamics. These gaps include: (1) the kinematic transformation from the orbital elements in the inertial frame to the states within the rotating frames, and (2) the dynamical disparity between the integrable averaged model and the non-integrable R3BPs. To bridge these gaps, symmetry serves as a critical criteri… view at source ↗
Figure 4
Figure 4. Figure 4: Symmetric apse configurations within the HRF (prograde, [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Potential symmetric apse configurations within the PIF. Abbreviations A, D, N, S denote ascending [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Representation of the modular evolution over a quarter period ( [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Process for transitioning DADM equilibria to analog symmetric POs within the unaveraged dynamics. [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Doubly symmetric (XOZ/YOZ) circular PO in the DADM (p [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Doubly symmetric (XOZ/YOZ) circular PO in the HR3BP (initialized from Fig. 8). [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Singly symmetric (XOZ) circular PO in the CR3BP (initialized from Fig. 8). [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Singly symmetric (XOZ) frozen PO in the DADM (p [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Singly symmetric (XOZ) frozen PO in the HR3BP (initialized from Fig. 11). [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Singly symmetric (XOZ) frozen PO in the CR3BP (initialized from Fig. 11). [PITH_FULL_IMAGE:figures/full_fig_p024_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Analysis of circular periodic orbits at i = 0◦ . 27 [PITH_FULL_IMAGE:figures/full_fig_p027_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Analysis of circular periodic orbits at i = 180◦ . 5.1.3 Period-Multiplying Bifurcations and Number of Branches Building on the planar-limit connections in Sections 5.1.1–5.1.2, this subsection characterizes how the spatial circular PO families branch from the planar g/LPO and f/DRO families. Owing to the circular geometry (e = 0), these spatial families emerge through period-multiplying bifurcations; the… view at source ↗
Figure 16
Figure 16. Figure 16: Analysis of frozen periodic orbits near i = 90◦ . 32 [PITH_FULL_IMAGE:figures/full_fig_p032_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Broken bifurcation between circular and frozen POs for [PITH_FULL_IMAGE:figures/full_fig_p034_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Archetypical bifurcation diagrams for unaveraged dynamics with different parities. [PITH_FULL_IMAGE:figures/full_fig_p035_18.png] view at source ↗
read the original abstract

Within a three-body system comprised of two celestial bodies and a spacecraft, the dynamical environment near a smaller primary is significantly perturbed, motivating a balance between global insight and model fidelity. While averaged dynamics offer an integrable model to classify solution landscapes, they inherently lack the accuracy of the unaveraged dynamics, such as the Hill Restricted Three-Body Problem and Circular Restricted Three-Body Problem. This work establishes a systematic bridge between the averaged and unaveraged regimes by explicitly linking averaged equilibria to symmetric periodic orbits in the unaveraged three-body systems. A unified frequency framework is introduced to characterize the mapping of invariant tori across the dynamical models. Leveraging the parity of the resonance ratio, an initialization scheme is developed to identify admissible apse configurations, enabling the a priori prediction of solution multiplicity and symmetry types. Furthermore, the global evolution of families derived from averaged equilibria is traced via bifurcation and frequency analysis. These findings are synthesized into archetypical bifurcation diagrams, providing a comprehensive atlas of the symmetric periodic orbit web within the HR3BP and CR3BP. The resulting framework not only clarifies the topological origins of complex periodic orbit families but also offers a versatile tool for trajectory design in cislunar and multi-body environments.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper claims to establish a systematic bridge between averaged and unaveraged three-body dynamics near smaller primaries by explicitly linking averaged equilibria to symmetric periodic orbits in the HR3BP and CR3BP. It introduces a unified frequency framework to characterize the mapping of invariant tori, develops an initialization scheme leveraging the parity of the resonance ratio to predict solution multiplicity and symmetry types, traces the global evolution of families via bifurcation and frequency analysis, and synthesizes the results into archetypical bifurcation diagrams providing an atlas of the symmetric periodic orbit web.

Significance. If the central mappings and diagrams hold, the work offers a practical atlas and versatile tool for trajectory design in cislunar and multi-body environments by clarifying topological origins of periodic orbit families. The systematic linking of averaged and unaveraged regimes via frequency analysis is a strength when supported by explicit constructions.

minor comments (2)
  1. The abstract asserts the existence of the bridge, framework, and bifurcation diagrams but supplies no derivations, error analysis, or verification steps; consider adding a dedicated section or appendix with explicit numerical verification of at least one mapped orbit family to strengthen the claims.
  2. The initialization scheme based on resonance parity is described at a high level; clarify the precise definition of 'parity' and the admissible apse configurations with an example in §3 or §4 to improve reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript on linking averaged and unaveraged three-body dynamics and for recommending minor revision. The summary accurately captures the central contributions of the unified frequency framework, resonance-parity initialization, and archetypical bifurcation diagrams.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's core claim is a systematic continuation from averaged equilibria to symmetric periodic orbits in the HR3BP/CR3BP using a frequency map and resonance-parity initialization. This is a standard dynamical-systems construction resting on external models (averaged three-body dynamics, continuation methods) rather than any self-definition, fitted-input prediction, or self-citation chain that reduces the result to its inputs. No load-bearing step is shown to be equivalent to its own assumptions by construction; the parity device is presented as an a-priori selection rule, not a derived output. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities can be extracted. The work appears to rest on standard assumptions of the restricted three-body problem.

pith-pipeline@v0.9.1-grok · 5744 in / 1175 out tokens · 29745 ms · 2026-06-27T18:10:33.731993+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

45 extracted references · 40 canonical work pages · 1 internal anchor

  1. [1]

    Fuller, E

    S. Fuller, E. Lehnhardt, C. Zaid, and K. Halloran. Gateway program status and overview.Journal of Space Safety Engineering, 9(4):625–628, 2022. doi:10.1016/j.jsse.2022.07.008

    work page doi:10.1016/j.jsse.2022.07.008 2022

  2. [2]

    Grasset, M

    O. Grasset, M. K. Dougherty, A. Coustenis, E. J. Bunce, C. Erd, D. Titov, M. Blanc, A. Coates, P. Drossart, L. N. Fletcher, et al. JUpiter ICy moons Explorer (JUICE): An ESA mission to orbit Ganymede and to characterise the Jupiter system.Planetary and Space Science, 78:1–21, 2013. doi:10.1016/j.pss.2012.12.002

    work page doi:10.1016/j.pss.2012.12.002 2013

  3. [3]

    J. R. Spencer and F. Nimmo. Enceladus: An active ice world in the Saturn system.Annual Review of Earth and Planetary Sciences, 41(1):693–717, 2013. doi:10.1146/annurev-earth-050212-124025

    work page doi:10.1146/annurev-earth-050212-124025 2013

  4. [4]

    Szebehely.Theory of orbits: The restricted problem of three bodies

    V . Szebehely.Theory of orbits: The restricted problem of three bodies. Elsevier, 2012

    2012

  5. [5]

    E., et al

    H. von Zeipel. Sur l’application des séries de m. lindstedt à l’étude du mouvement des comètes périodiques.Astronomische Nachrichten, 183:345–418, 1910. doi:10.1002/asna.19091832202

    work page doi:10.1002/asna.19091832202 1910

  6. [6]

    M. L. Lidov. Evolution of the orbits of artificial satellites of planets as affected by gravitational perturbation from external bodies.Artificial Earth Satellites, 8:5–45, 1961. Originally in Russian. English translation:Planetary and Space Science,9, 719–759, 1962, https://doi.org/10.1016/00 32-0633(62)90129-0

    work page doi:10.1016/00 1961

  7. [7]

    Y . Kozai. Secular perturbations of asteroids with high inclination and eccentricity.The Astronomical Journal, 67:591–598, 1962. doi:10.1086/108790

    work page doi:10.1086/108790 1962

  8. [8]

    Ito and K

    T. Ito and K. Ohtsuka. The Lidov-Kozai oscillation and Hugo von Zeipel.arXiv preprint arXiv:1911.03984, 2019. doi:10.48550/arXiv.1911.03984

    work page doi:10.48550/arxiv.1911.03984 1911

  9. [9]

    D. J. Scheeres, M. D. Guman, and B. F. Villac. Stability analysis of planetary satellite orbiters: application to the Europa orbiter.Journal of Guidance, Control, and Dynamics, 24(4):778–787, 2001. doi:10.2514/2.4778

    work page doi:10.2514/2.4778 2001

  10. [10]

    A. F. B. de Almeida Prado. Third-body perturbation in orbits around natural satellites.Journal of Guidance, Control, and Dynamics, 26(1):33–40, 2003. doi:10.2514/2.5042

    work page doi:10.2514/2.5042 2003

  11. [11]

    R. A. Broucke. Long-term third-body effects via double averaging.Journal of Guidance, Control, and Dynamics, 26(1):27–32, 2003. doi:10.2514/2.5041

    work page doi:10.2514/2.5041 2003

  12. [12]

    T. A. Ely. Stable constellations of frozen elliptical inclined lunar orbits.The Journal of the Astronautical Sciences, 53:301–316, 2005. doi:10.1007/BF03546355

    work page doi:10.1007/bf03546355 2005

  13. [13]

    Lunar Frozen Orbits

    D. Folta and D. Quinn. Lunar frozen orbits. InAIAA/AAS Astrodynamics Specialist Conference and Exhibit, page 6749, 2006. doi:10.2514/6.2006-6749

    work page doi:10.2514/6.2006-6749 2006

  14. [14]

    T. A. Ely and E. Lieb. Constellations of elliptical inclined lunar orbits providing polar and global coverage.The Journal of the Astronautical Sciences, 54(1):53–67, 2006. doi:10.1007/BF03256476

    work page doi:10.1007/bf03256476 2006

  15. [15]

    Lara and R

    M. Lara and R. Russell. Computation of a science orbit about Europa.Journal of Guidance, Control, and Dynamics, 30(1):259–263, 2007. doi:10.2514/1.22493

    work page doi:10.2514/1.22493 2007

  16. [16]

    R. P. Russell and A. T. Brinckerhoff. Circulating eccentric orbits around planetary moons.Journal of Guidance, Control, and Dynamics, 32(2):424–436, 2009. doi:10.2514/1.38593

    work page doi:10.2514/1.38593 2009

  17. [17]

    McArdle and R

    S. McArdle and R. P. Russell. Circulating, eccentric periodic orbits at the Moon.Celestial Mechanics and Dynamical Astronomy, 133(4):18, 2021. doi:10.1007/s10569-021-10013-z

    work page doi:10.1007/s10569-021-10013-z 2021

  18. [18]

    C. J. Franz and R. P. Russell. Database of planar and three-dimensional periodic orbits and families near the Moon.The Journal of the Astronautical Sciences, 69(6):1573–1612, 2022. doi:10.1007/s40295-022- 00361-9. 36

    work page doi:10.1007/s40295-022- 2022

  19. [19]

    K. C. Howell and J. V . Breakwell. Almost rectilinear halo orbits.Celestial mechanics, 32(1):29–52,

  20. [20]

    doi:10.1007/BF01358402

    work page doi:10.1007/bf01358402

  21. [21]

    E. M. Zimovan-Spreen, K. C. Howell, and D. C. Davis. Near rectilinear halo orbits and nearby higher-period dynamical structures: orbital stability and resonance properties.Celestial Mechanics and Dynamical Astronomy, 132(5):28, 2020. doi:10.1007/s10569-020-09968-2

    work page doi:10.1007/s10569-020-09968-2 2020

  22. [22]

    R. P. Russell. Global search for planar and three-dimensional periodic orbits near Europa.The Journal of the Astronautical Sciences, 54(2):199–226, 2006. doi:10.1007/BF03256483

    work page doi:10.1007/bf03256483 2006

  23. [23]

    R. L. Restrepo and R. P. Russell. A database of planar axisymmetric periodic orbits for the solar system. Celestial Mechanics and Dynamical Astronomy, 130(7):49, 2018. doi:10.1007/s10569-018-9844-6

    work page doi:10.1007/s10569-018-9844-6 2018

  24. [24]

    Moreno, C

    A. Moreno, C. Aydin, O. van Koert, U. Frauenfelder, and D. Koh. Bifurcation graphs for the CR3BP via symplectic methods: On the Jupiter–Europa and Saturn–Enceladus systems.The Journal of the Astronautical Sciences, 71(6):51, 2024. doi:10.1007/s40295-024-00462-7

    work page doi:10.1007/s40295-024-00462-7 2024

  25. [25]

    C. Aydin. Exploration of vertical self-resonant bifurcations from DRO in the Earth-Moon CR3BP. arXiv preprint arXiv:2508.17286, 2025. doi:10.48550/arXiv.2508.17286

    work page doi:10.48550/arxiv.2508.17286 2025

  26. [26]

    Aydin and A

    C. Aydin and A. Batkhin. Studying network of symmetric periodic orbit families of the Hill prob- lem via symplectic invariants.Celestial Mechanics and Dynamical Astronomy, 137(2):1–77, 2025. doi:10.1007/s10569-025-10241-7

    work page doi:10.1007/s10569-025-10241-7 2025

  27. [27]

    K. C. Howell, D. J. Grebow, and Z. P. Olikara. Design using Gauss’ perturbing equations with applications to lunar south pole coverage. InAAS/AIAA Space Flight Mechanics Meeting, Sedona, Arizona, January 2007

    2007

  28. [28]

    R. P. Russell and M. Lara. Long-lifetime lunar repeat ground track orbits.Journal of Guidance, Control, and Dynamics, 30(4):982–993, 2007. doi:10.2514/1.27104

    work page doi:10.2514/1.27104 2007

  29. [29]

    Lara and J

    M. Lara and J. F. San Juan. Dynamic behavior of an orbiter around Europa.Journal of Guidance, Control, and Dynamics, 28(2):291–297, 2005. doi:10.2514/1.5686

    work page doi:10.2514/1.5686 2005

  30. [30]

    M. Lara, R. Russell, and B. Villac. Classification of the distant stability regions at Europa.Journal of Guidance, Control, and Dynamics, 30(2):409–418, 2007. doi:10.2514/1.22372

    work page doi:10.2514/1.22372 2007

  31. [31]

    D. C. Koblick and P. Kelly. Novel three-body tulip-shaped orbit families for lunar missions.The Journal of the Astronautical Sciences, 72(4):32, 2025. doi:10.1007/s40295-025-00510-w

    work page doi:10.1007/s40295-025-00510-w 2025

  32. [32]

    L. Peng, P. Shi, Y . Liang, and N. Pushparaj. Analog of lunar Sun-synchronous orbits based on spatial distant retrograde orbits.Journal of Guidance, Control, and Dynamics, 48(6):1439–1448, 2025. doi:10.2514/1.G008375

    work page doi:10.2514/1.g008375 2025

  33. [33]

    Zhang, X

    J. Zhang, X. Jiang, Y . Yuan, and H. Dai. Time-regularized bifurcation framework for constructing tulip orbits in the CRTBP.Nonlinear Dynamics, 114(9):620, 2026. doi:10.1007/s11071-026-12465-0

    work page doi:10.1007/s11071-026-12465-0 2026

  34. [34]

    D. A. Vallado.Fundamentals of astrodynamics and applications, volume 12. Springer Science & Business Media, 2001

    2001

  35. [35]

    J. M. Longuski, F. R. Hoots, and G. E. Pollock IV .Introduction to Orbital Perturbations, volume 40. Springer Nature, 2022. doi:10.1007/978-3-030-89758-1

    work page doi:10.1007/978-3-030-89758-1 2022

  36. [36]

    Lunar frozen orbits revisited.Celestial Mechanics and Dynamical Astronomy, 130(10): 61, September 2018

    T. Nie and P. Gurfil. Lunar frozen orbits revisited.Celestial Mechanics and Dynamical Astronomy, 130 (10):61, 2018. doi:10.1007/s10569-018-9858-0

    work page doi:10.1007/s10569-018-9858-0 2018

  37. [37]

    V . I. Arnold.Mathematical methods of classical mechanics, volume 60. Springer Science & Business Media, 2013. doi:10.1007/978-1-4757-1693-1

    work page doi:10.1007/978-1-4757-1693-1 2013

  38. [38]

    Celletti

    A. Celletti. From infinite to finite time stability in celestial mechanics and astrodynamics.Astrophysics and Space Science, 368(12):106, 2023. doi:10.1007/s10509-023-04264-5

    work page doi:10.1007/s10509-023-04264-5 2023

  39. [39]

    C. Aydin. The linear symmetries of Hill’s lunar problem.Archiv der Mathematik, 120(3):321–330,

  40. [40]

    doi:10.1007/s00013-022-01822-1

    work page doi:10.1007/s00013-022-01822-1

  41. [41]

    I. A. Robin and V . V . Markellos. Numerical determination of three-dimensional periodic orbits generated from vertical self-resonant satellite orbits.Celestial Mechanics, 21(4):395–434, 1980. doi:10.1007/BF01231276. 37

    work page doi:10.1007/bf01231276 1980

  42. [42]

    Springer Berlin Heidelberg, 2002.isbn: 9783540459491

    M. Hénon.Generating families in the restricted three-body problem. Springer, 2002. doi:10.1007/3- 540-69650-4

    work page doi:10.1007/3- 2002

  43. [43]

    Strömgren

    E. Strömgren. Connaissance actuelle des orbites dans le problème des trois corps.Bulletin astronomique, Observatoire de Paris, 9(1):87–130, 1933. doi:10.3406/bastr.1933.14090. English translation by Maxime Murray and J. D. Mireles James (2018)

    work page doi:10.3406/bastr.1933.14090 1933

  44. [44]

    Guzzetti, N

    D. Guzzetti, N. Bosanac, A. Haapala, K. C. Howell, and D. C. Folta. Rapid trajectory design in the Earth–Moon ephemeris system via an interactive catalog of periodic and quasi-periodic orbits.Acta Astronautica, 126:439–455, 2016. doi:10.1016/j.actaastro.2016.06.029

    work page doi:10.1016/j.actaastro.2016.06.029 2016

  45. [45]

    Bifurcations of highly inclined near halo orbits using Moser regularization

    C. Joung, D. Koh, and O. van Koert. Bifurcations of highly inclined near halo orbits using Moser regularization.arXiv preprint arXiv:2512.03849, 2025. doi:10.48550/arXiv.2512.03849. 38

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2512.03849 2025