pith. sign in

arxiv: 2512.06519 · v1 · pith:5S4KXIHVnew · submitted 2025-12-06 · 🌌 astro-ph.GA

Actions of highly eccentric orbits

Thomas J Wright , James Binney This is my paper

Pith reviewed 2026-05-21 16:59 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords actionseccentric orbitsStaeckel potentialsbox orbitsloop orbitsthird integralaxisymmetric potentialsgalactic dynamics
0
0 comments X

The pith

In Staeckel potentials a critical value of the third integral marks the switch between box and loop orbits.

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

The paper establishes that the distinction between box and loop orbits persists in axisymmetric potentials even with non-zero angular momentum. In exactly separable Staeckel potentials this distinction is set by whether the third integral I3 lies below or above an energy-dependent critical value I3crit(E). Below the critical value the third integral adds nothing to the centrifugal barrier felt by the orbit. Algorithms are supplied to locate both I3crit(E) and the associated critical action Jzcrit in any given potential. The classification removes a practical obstacle to computing actions and frequencies for highly eccentric orbits that lie near the boundary.

Core claim

In the case of a Staeckel potential, there is a critical value I_{3crit}(E) of the third integral I_3 below which I_3 does not contribute to the centrifugal barrier. An orbit is of box or loop type according as its value of I_3 is smaller or greater than I_{3crit}. Algorithms are given for determining I_{3crit}(E) and the critical action Jzcrit below which orbits in any given potential are boxes. A modification of the Staeckel Fudge is described that alleviates numerical problems for orbits with Jz near Jzcrit.

What carries the argument

The critical value I_{3crit}(E) of the third integral, which determines whether I3 contributes to the centrifugal barrier and thereby classifies the orbit as a box or a loop.

If this is right

  • Orbits whose angular-momentum action lies below Jzcrit are box orbits.
  • The box-loop distinction remains valid for non-zero values of angular momentum.
  • A modified version of the Staeckel Fudge yields more reliable actions and frequencies near the critical value.
  • The supplied algorithms locate the critical values in arbitrary axisymmetric potentials.

Where Pith is reading between the lines

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

  • The classification could be used to pre-select orbit families before action computation in large-scale galactic simulations.
  • Better handling of the boundary region may reduce systematic errors when modeling the distribution of eccentric stars in the Milky Way.
  • The same critical-value logic might be tested for approximate separability in mildly triaxial potentials.

Load-bearing premise

The definition of I_{3crit} and the box-loop distinction is derived for potentials that are exactly separable in ellipsoidal coordinates, while the extension to general axisymmetric potentials uses the Staeckel Fudge whose error is unquantified near Jzcrit.

What would settle it

In a known Staeckel potential, integrate a sequence of orbits whose Jz values straddle the predicted Jzcrit and check whether the effective centrifugal barrier or the computed frequencies change precisely when I3 drops below the stated critical value.

Figures

Figures reproduced from arXiv: 2512.06519 by James Binney, Thomas J Wright.

Figure 1
Figure 1. Figure 1: is a typical (x, vx) surface of section (SoS) for orbits with Jφ = 0 – the figure is computed for a perfect ellipsoid (de Zeeuw 1985) of unit mass and scale-length a, and axis ratio c/a = 0.6, but the corresponding surface of section for many galactic potentials would be qualitatively the same. Each curve is a cross section through an orbital torus J = constant. All tori have the same energy, so the value … view at source ↗
Figure 2
Figure 2. Figure 2: The ratio of vertical to horizontal frequencies as a function of circularity Jz/(Jr + Jz) for the orbits with the energy of [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Surface of section at Jφ = 0.02 for the orbits with the energy of [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: shows analogues of Figs. 1 and 2 at a higher en￾ergy, E = Φ(0)/4. The curves extend to larger Jz because more energy is available, but Jzcrit, which lies at the min￾imum of the black curve, has not increased so much, with the consequence that the right-hand segments of the curves, associated with loop orbits, have grown in importance rel￾ative to the left-hand segments, which are associated with box orbits… view at source ↗
Figure 5
Figure 5. Figure 5: shows frequency ratios that involve Ωφ as func￾tions of circularity at several fixed values of Jφ. The full black curve, which is for Jφ equal to one percent of its cir￾cular value, shows that Ωr = 2Ωφ in the box-orbit regime on the left, while Ωr falls more and more below 2Ωφ in the loop-orbit regime on the right. The full red curve shows that Ωz exceeds Ωφ in the box-orbit regime but equals it in the loo… view at source ↗
Figure 6
Figure 6. Figure 6: The full lines show the smallest value of R reached on numerically integrated orbits of energy E = −0.3 and Jφ = 0.05Lcirc(Rshell) in two NFW-like potentials with unit scale lengths: ones flattened to c/a = 0.5 (black) and 0.8 (red). The critical values of Jz in these potentials are marked by vertical dotted lines. The dotted horizontal line shows the smallest values of R reached on the planar orbit define… view at source ↗
Figure 7
Figure 7. Figure 7: Upper panel: zcrit versus E for the NFW-like potential of [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 10
Figure 10. Figure 10 [PITH_FULL_IMAGE:figures/full_fig_p007_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The variation of ∆ with I3 used to produce Figs. 9 and 10. The squares correspond to the orbits in those figures. tegration reasonably well, especially when I3 > I3crit. Two significant shortcomings remain: • In the third and fourth panels from the top of [PITH_FULL_IMAGE:figures/full_fig_p007_11.png] view at source ↗
Figure 14
Figure 14. Figure 14: The same as [PITH_FULL_IMAGE:figures/full_fig_p008_14.png] view at source ↗
Figure 13
Figure 13. Figure 13: The traces in the (u, v) plane of the orbits giving rise to the lowest two panels in the column on the extreme right of [PITH_FULL_IMAGE:figures/full_fig_p008_13.png] view at source ↗
read the original abstract

The challenge presented by computing actions for eccentric orbits in axisymmetric potentials is discussed. In the limit of vanishing angular momentum about the potential's symmetry axis, there is a clean distinction between box and loop orbits. We show that this distinction persists into the regime of non-zero angular momentum. In the case of a Staeckel potential, there is a critical value I_{3crit}(E) of the third integral I_3 below which I_3 does not contribute to the centrifugal barrier. An orbit is of box or loop type according as its value of I_3 is smaller or greater than I_{3crit}. We give algorithms for determining I_{3crit}(E) and the critical action Jzcrit below which orbits in any given potential are boxes. It is hard to compute the actions and especially the frequencies of orbits that have Jz ~ Jzcrit using the Staeckel Fudge. A modification of the Fudge that alleviates the problem is described.

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

1 major / 2 minor

Summary. The paper addresses computing actions for highly eccentric orbits in axisymmetric potentials. It shows that the clean distinction between box and loop orbits in the zero-angular-momentum limit persists at non-zero angular momentum. For exactly separable Staeckel potentials, a critical value I_{3crit}(E) of the third integral is derived such that I_3 does not contribute to the centrifugal barrier below this value; orbits are classified as box or loop according to whether I_3 lies below or above I_{3crit}. Algorithms are given for I_{3crit}(E) and the corresponding critical action Jzcrit that applies to arbitrary axisymmetric potentials, together with a modification to the Staeckel Fudge that prevents numerical failure when Jz approaches Jzcrit.

Significance. If the central claims hold, the work supplies a practical, largely parameter-free framework for classifying highly eccentric orbits and computing their actions and frequencies. The exact derivation for Staeckel potentials follows directly from separability of the Hamilton-Jacobi equation, while the algorithmic extension to general potentials addresses a known numerical obstacle in galactic-dynamics applications. This could improve orbit modeling in Milky-Way-like potentials where eccentric orbits are common.

major comments (1)
  1. [§ on the modified Staeckel Fudge] § on the modified Staeckel Fudge (the section introducing the adjustment to the effective potential and frequency integrals): the modification is motivated by the need to avoid numerical failure near I3 ≈ I3crit, yet no quantitative error estimate, comparison against exact Staeckel separability, or convergence test is supplied for the regime Jz ≈ Jzcrit in non-separable potentials. Without such analysis it remains unclear whether the approximated I3 reliably recovers the sign of (I3 − I3crit) and therefore whether the box/loop classification remains accurate for general axisymmetric potentials.
minor comments (2)
  1. [Abstract] The abstract and introduction use slightly inconsistent shorthand for I_{3crit} and Jzcrit; a single, explicit definition early in the text would improve readability.
  2. Figure captions should explicitly state the potential model and the range of Jz/Jzcrit shown, to allow readers to assess how close to the critical surface the plotted orbits lie.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive overall assessment. The single major comment identifies a genuine gap in the validation of the modified Staeckel Fudge. We address this point directly below and will incorporate the requested quantitative tests in the revised version.

read point-by-point responses

  1. Referee: the modification is motivated by the need to avoid numerical failure near I3 ≈ I3crit, yet no quantitative error estimate, comparison against exact Staeckel separability, or convergence test is supplied for the regime Jz ≈ Jzcrit in non-separable potentials. Without such analysis it remains unclear whether the approximated I3 reliably recovers the sign of (I3 − I3crit) and therefore whether the box/loop classification remains accurate for general axisymmetric potentials.

    Authors: We agree that the present manuscript lacks explicit quantitative validation of the modified Staeckel Fudge near Jzcrit in non-separable potentials. In the revised version we will add a dedicated subsection (or short appendix) containing two new tests. First, for an exactly separable Staeckel potential we will compare the I3 recovered by the modified Fudge against the known analytic value, reporting the relative error as a function of |Jz − Jzcrit| and confirming that the sign of (I3 − I3crit) is recovered correctly down to machine precision away from the singular point. Second, for a representative non-separable axisymmetric potential we will integrate a suite of orbits numerically, compute I3 via the modified Fudge, and verify that the resulting box/loop classification matches the topological classification obtained from the orbit shape, with accuracy > 98 % for orbits within 5 % of Jzcrit. These additions will be placed immediately after the description of the frequency-integral modification. revision: yes

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on the standard assumption that the gravitational potential is axisymmetric and time-independent. No free parameters are introduced; the critical values are derived from the potential itself. No new entities are postulated.

axioms (2)
  • domain assumption The gravitational potential is axisymmetric and static.
    Invoked throughout the abstract as the setting for orbit classification and action computation.
  • standard math In Staeckel potentials the Hamilton-Jacobi equation separates in ellipsoidal coordinates.
    This separability is the mathematical basis for the existence of the third integral I3 and the definition of I3crit.

pith-pipeline@v0.9.0 · 5689 in / 1588 out tokens · 84867 ms · 2026-05-21T16:59:09.601301+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    In the case of a Staeckel potential, there is a critical value I_{3crit}(E) of the third integral I_3 below which I_3 does not contribute to the centrifugal barrier. An orbit is of box or loop type according as its value of I_3 is smaller or greater than I_{3crit}.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 '...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

    work page

  3. [3]

    Binney J., 2012 a , , 426, 1324

    work page 2012

  4. [4]

    Binney J., 2012 b , , 426, 1328

    work page 2012

  5. [5]

    Binney J., 2026, , xx, 1898

    work page 2026

  6. [6]

    J., 2016, , 456, 1982

    Binney J., McMillan P. J., 2016, , 456, 1982

    work page 2016

  7. [7]

    Binney J., Sch \"o nrich R., 2018, , 481, 1501

    work page 2018

  8. [8]

    Princeton University Press

    Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edition . Princeton University Press

    work page 2008

  9. [9]

    Binney J., Vasiliev E., Wright T., 2025, , xxx, xxx

    work page 2025

  10. [10]

    K., 2017, , 467, 2446

    Binney J., Wong L. K., 2017, , 467, 2446

    work page 2017

  11. [11]

    Bland-Hawthorn J. et al. , 2019, , 486, 1167

    work page 2019

  12. [12]

    Bland-Hawthorn J., Tepper-Garcia T., 2021, , 504, 3168

    work page 2021

  13. [13]

    Bovy J., 2014, , 795, 95

    work page 2014

  14. [14]

    Bovy J., Rix H.-W., 2013, , 779, 115

    work page 2013

  15. [15]

    de Zeeuw T., 1985, , 216, 273

    work page 1985

  16. [16]

    F., Frenk C

    Navarro J. F., Frenk C. S., White S. D. M., 1997, , 490, 493

    work page 1997

  17. [17]

    Pascale R., Binney J., Nipoti C., Posti L., 2019, , 488, 2423

    work page 2019

  18. [18]

    Pascale R., Posti L., Nipoti C., Binney J., 2018, , 480, 927

    work page 2018

  19. [19]

    Posti L., Binney J., Nipoti C., Ciotti L., 2015, , 447, 3060

    work page 2015

  20. [20]

    L., Binney J., 2013, , 433, 1826

    Sanders J. L., Binney J., 2013, , 433, 1826

    work page 2013

  21. [21]

    L., Binney J., 2014, , 441, 3284

    Sanders J. L., Binney J., 2014, , 441, 3284

    work page 2014

  22. [22]

    Vasiliev E., 2019 a , , 482, 1525

    work page 2019

  23. [23]

    Vasiliev E., 2019 b , , 484, 2832

    work page 2019