Extended Brown Hamiltonian Dynamics

• The extended Brown Hamiltonian is a secular, integrable model for weakly hierarchical three-body systems that incorporates nonlinear short-period effects into a closed-form framework.
• It combines the classical quadrupole disturbing function with both Brown’s correction and a novel extended term to accurately describe modified von Zeipel–Lidov–Kozai oscillations.
• The model provides precise analytical predictions and diagnostics, validated by numerical simulations, for applications such as the dynamics of irregular satellites.

Searching arXiv for the specified papers and closely related work on the extended Brown Hamiltonian and Brown corrections. The extended Brown Hamiltonian is a secular, integrable Hamiltonian model for weakly hierarchical three-body systems in the test-particle limit, developed to describe modified von Zeipel–Lidov–Kozai (ZLK) oscillations when classical double-averaged quadrupole theory and the classical Brown correction are no longer sufficiently accurate. In the formulation introduced by Lei and Grishin, the model retains the quadrupole-order disturbing function but incorporates nonlinear short-period effects accumulated from both the outer and inner orbits, yielding a closed-form Hamiltonian and a closed-form transformation between mean and osculating elements (Lei et al., 19 May 2025). Subsequent work analyzed the model’s phase-space structure, the associated modified Lidov integral $C_{\rm ZLK}$, and applications to irregular satellites of the giant planets, where the model provides an efficient criterion for identifying ZLK libration in the weak-hierarchy regime (Lei et al., 15 Jul 2025, Lei et al., 13 May 2026).

1. Conceptual setting and scope

The extended Brown Hamiltonian is formulated for a hierarchical triple in the test-particle limit, with an inner body of negligible mass orbiting a central mass $m_0$ and an outer perturber of mass $m_p$ moving on a fixed Kepler orbit (Lei et al., 19 May 2025). The relevant dynamical context is a low-hierarchy or weak-hierarchy triple, where nonlinear perturbations associated with periodic oscillations of both the inner and outer binaries accumulate secularly and materially affect long-term dynamics (Lei et al., 19 May 2025, Lei et al., 15 Jul 2025).

The motivating systems explicitly include high-altitude lunar satellites influenced by the Earth, planetary satellites perturbed by the Sun, and stellar binaries affected by a supermassive black hole (Lei et al., 19 May 2025). In such systems, the classical double-averaged quadrupole Hamiltonian captures the leading ZLK mechanism, while Brown’s 1936 correction accounts for first-order nonlinear short-period effects of the outer orbit; however, these contributions alone can become insufficient when the hierarchy is mild and the inner-orbital period is also comparable to the ZLK timescale (Lei et al., 19 May 2025).

In this framework, the extended Brown Hamiltonian is not a higher-multipole expansion of the disturbing function. Rather, it is obtained by truncating the disturbing function at quadrupole order while retaining all nonlinearities up to $\mathcal O\!\left((n_p/n)^2\right)$ produced by successive von Zeipel averaging (Lei et al., 19 May 2025). A plausible implication is that its improvement over classical secular theory arises from a more complete treatment of short-period accumulation rather than from adding octupole or higher spatial harmonics.

2. Hamiltonian structure and closed-form terms

The basic canonical variables are Delaunay variables. In one normalization, these are written as

$$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$

with $\mu=Gm_0$ (Lei et al., 19 May 2025). In the normalized formulation used in the analytical study, one works with

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$

where $\{g,G\}$ and $\{h,H\}$ are canonically conjugate pairs and $H$ is an exact integral because $m_0$0 is cyclic (Lei et al., 15 Jul 2025).

After double averaging and retention of nonlinear corrections through second order in the mean-motion ratio, the secular Hamiltonian can be written, up to an additive or overall constant factor, as

$m_0$1

where

$m_0$2

and the three dimensionless components are given in closed form (Lei et al., 19 May 2025).

The classical quadrupole term is

$m_0$3

and the classical Brown correction is

$m_0$4

The new extended contribution is

$m_0$5

These formulae are stated to be closed form with respect to the eccentricities of the inner and outer orbits, and the generating function used for the mean–osculating transformation is likewise closed form and converges for any $m_0$6, unlike power-series expansions (Lei et al., 19 May 2025).

A second, algebraically equivalent representation emphasizes the constant

$m_0$7

where $m_0$8 is the mutual inclination at zero eccentricity, and introduces

$m_0$9

In that representation,

$m_p$0

with

$m_p$1

$m_p$2

$m_p$3

which is the form used in the applications to irregular satellites (Lei et al., 13 May 2026).

3. Ordering, derivation, and approximation regime

The derivation assumes a negligible-mass inner body, a fixed outer Kepler perturber, and a hierarchical semimajor-axis ratio $m_p$4, while truncating the disturbing function at quadrupole order $m_p$5 (Lei et al., 19 May 2025). A bookkeeping parameter $m_p$6 is introduced, and the Hamiltonian is expanded to second order in $m_p$7 (Lei et al., 19 May 2025).

Two small parameters control the corrections: $m_p$8 with $m_p$9 associated with the classical Brown correction and $\mathcal O\!\left((n_p/n)^2\right)$0 with the new extended term (Lei et al., 19 May 2025). In the irregular-satellite scaling used later, these parameters are expressed as

$\mathcal O\!\left((n_p/n)^2\right)$1

where $\mathcal O\!\left((n_p/n)^2\right)$2 (Lei et al., 13 May 2026).

The derivation proceeds by two successive von Zeipel canonical transformations (Lei et al., 19 May 2025). In Step I, one removes the inner fast angle $\mathcal O\!\left((n_p/n)^2\right)$3 via a generating function $\mathcal O\!\left((n_p/n)^2\right)$4, producing the single-averaged Hamiltonian and then removing its short-period remainder. In Step II, one removes the outer fast angle $\mathcal O\!\left((n_p/n)^2\right)$5 via $\mathcal O\!\left((n_p/n)^2\right)$6, producing the double-averaged Hamiltonian $\mathcal O\!\left((n_p/n)^2\right)$7 and then removing its short-period remainder, which generates $\mathcal O\!\left((n_p/n)^2\right)$8. The second-order terms in the final averaging produce two secular corrections: the classical Brown term $\mathcal O\!\left((n_p/n)^2\right)$9, from accumulated outer short-period effects, and the new extended term $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$0, from accumulated inner short-period effects (Lei et al., 19 May 2025).

The final secular problem is integrable, with a single degree of freedom in $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$1 and with $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$2 and $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$3 conserved (Lei et al., 19 May 2025). This is central to the later analytical development of fixed points, separatrices, maximum eccentricities, and critical inclinations (Lei et al., 15 Jul 2025).

The reported validity regime is “mildly to weakly hierarchical triples,” with truncation at $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$4 stated to be valid up to $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$5 (Lei et al., 13 May 2026). This suggests that the model is intended precisely for systems that are not asymptotically deep in the hierarchy, but are still hierarchical enough for the quadrupole truncation to remain meaningful.

4. Mean–osculating transformation and secular equations

A distinctive feature of the framework is that it provides not only the secular Hamiltonian but also a closed-form canonical generating function for transforming between osculating elements $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$6 and secular or mean elements $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$7 (Lei et al., 19 May 2025). To first order in $$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$8, the generating function is written as

$$L=\sqrt{\mu a},\quad G=L\sqrt{1-e^2},\quad H=G\cos i,\quad (l,g,h)=(M,\omega,\Omega),$$9

where $\mu=Gm_0$0 removes inner-period terms and $\mu=Gm_0$1 removes outer-period terms (Lei et al., 19 May 2025).

The resulting change of variables is given in Lagrange-planetary-equations form, evaluated at the mean elements: $\mu=Gm_0$2

$\mu=Gm_0$3

$\mu=Gm_0$4

with $\mu=Gm_0$5 and all partial derivatives taken at the mean elements (Lei et al., 19 May 2025).

The secular equations of motion under $\mu=Gm_0$6 are

$\mu=Gm_0$7

and may be decomposed as

$\mu=Gm_0$8

with each contribution available in closed form (Lei et al., 19 May 2025).

The paper states that these expressions may be inserted directly into a 1D ODE integrator or used to find analytical Jacobian-elliptic solutions (Lei et al., 19 May 2025). This indicates that the model is not only conceptually integrable but also practically deployable as a reduced dynamical system.

5. Modified Lidov integral and phase-space organization

A central analytical development in later papers is the modified Lidov integral $\mu=Gm_0$9, introduced to characterize the modified ZLK dynamics under the extended Hamiltonian (Lei et al., 15 Jul 2025, Lei et al., 13 May 2026). In the classical quadrupole problem, the invariant is

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$0

with

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$1

and equivalently

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$2

in the normalized notation (Lei et al., 15 Jul 2025).

In the extended model,

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$3

with

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$4

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$5

$$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$6

and one checks directly that $$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$7 and $$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$8, so $$g=\omega,\quad G=\sqrt{1-e^2},\quad h=\Omega,\quad H=G\cos i,$$9 is a second integral of the 1-degree-of-freedom system (Lei et al., 15 Jul 2025).

In the alternative $\{g,G\}$0-based representation, the same quantity satisfies

$\{g,G\}$1

which implies $\{g,G\}$2 exactly (Lei et al., 13 May 2026). The later applications paper explicitly calls $\{g,G\}$3 a “practical diagnostic index” and a “decisive parameter” for identifying the ZLK resonance (Lei et al., 13 May 2026).

The separatrix between libration and circulation is characterized by

$\{g,G\}$4

and, in the $\{g,G\}$5 phase portrait, the level curve $\{g,G\}$6 passes through the unstable fixed point at $\{g,G\}$7 (Lei et al., 15 Jul 2025, Lei et al., 13 May 2026). Trajectories with $\{g,G\}$8 correspond to full circulation of $\{g,G\}$9 from $\{h,H\}$0 to $\{h,H\}$1, while trajectories with $\{h,H\}$2 are confined to libration around $\{h,H\}$3 or $\{h,H\}$4 (Lei et al., 13 May 2026). Accordingly, the analytic criterion

$\{h,H\}$5

is used to identify ZLK libration in the weak-hierarchy regime (Lei et al., 13 May 2026).

A recurrent theme is that the Brown correction breaks the prograde–retrograde symmetry of classical quadrupole ZLK dynamics. Under the classical model, $\{h,H\}$6 is even in $\{h,H\}$7, implying invariance under $\{h,H\}$8. The second-order correction $\{h,H\}$9 is also even in $H$0, but the first-order Brown term $H$1 is odd in $H$2, so the full Hamiltonian is asymmetric under $H$3 (Lei et al., 15 Jul 2025). The same asymmetry appears in the Lidov integrals, where only $H$4 is odd in $H$5 (Lei et al., 15 Jul 2025). This is the formal reason that ZLK properties in prograde and retrograde regimes are no longer symmetric.

6. Analytical results for modified ZLK dynamics

The analytical study in Paper II develops perturbative expressions for the principal geometric and dynamical features of modified ZLK oscillations under the extended Brown Hamiltonian (Lei et al., 15 Jul 2025). The four highlighted characteristics are the location of the ZLK center, the maximum eccentricity on the separatrix, the critical inclination for the onset of libration, and the separatrix boundary itself.

For the stable fixed point, stationarity with $H$6 implies $H$7, and $H$8 leads to a quartic in

$H$9

The perturbative solution is written as

$m_0$00

with

$m_0$01

and the corresponding fixed-point eccentricity and inclination given by

$m_0$02

(Lei et al., 15 Jul 2025).

For the separatrix maximum eccentricity, one defines

$m_0$03

and imposes

$m_0$04

or equivalently $m_0$05, obtaining a quartic in $m_0$06 (Lei et al., 15 Jul 2025). The perturbation solution is

$m_0$07

from which

$m_0$08

follow (Lei et al., 15 Jul 2025).

For the critical inclination, the onset of libration is determined by the small-eccentricity saddle conditions

$m_0$09

leading to a quartic in

$m_0$10

The perturbative solution is

$m_0$11

with zeroth-order term

$m_0$12

where $m_0$13 corresponds to prograde and $m_0$14 to retrograde motion, and

$m_0$15

(Lei et al., 15 Jul 2025).

The perturbation scheme itself is developed in two variants. “Approach I” treats both $m_0$16 and $m_0$17 as small and expands in a double series. “Approach II” absorbs $m_0$18 exactly by starting from $m_0$19 and then treats only $m_0$20 as small (Lei et al., 15 Jul 2025). The reported quoted results are those of Approach II, expanded to third order in $m_0$21, and the paper states that analytical predictions show excellent agreement with numerical results (Lei et al., 15 Jul 2025).

7. Applications to irregular satellites and empirical performance

The principal numerical tests reported for Paper I concern four high-altitude moons of Jupiter: Pasiphae, Kore, Callirrhoe, and Philophrosyne (Lei et al., 19 May 2025). For these systems, the reported parameters are approximately $m_0$22 and $m_0$23 (Lei et al., 19 May 2025). The paper compares four models for Pasiphae and Kore: the classical $m_0$24 only, $m_0$25, the full extended model $m_0$26, and direct $m_0$27-body integration (Lei et al., 19 May 2025). It reports that only the full extended model captures both the period and amplitude well, with errors $m_0$28, whereas the truncated models have $m_0$29 errors and can even predict the wrong libration-versus-circulation behavior; similar agreement is reported for all four satellites over thousands of years (Lei et al., 19 May 2025).

Paper II extends this numerical validation to extensive $m_0$30-body integrations over grids in $m_0$31, up to $m_0$32 secular ZLK periods (Lei et al., 15 Jul 2025). The analytic center, boundary, and critical-inclination curves are stated to lie on top of numerically identified libration islands with typical relative errors

$m_0$33

for Approach II, and

$m_0$34

even into the weak-hierarchy regime $m_0$35 (Lei et al., 15 Jul 2025).

Paper III applies the modified Lidov integral $m_0$36 to the known population of irregular satellites of the giant planets (Lei et al., 13 May 2026). Using mean orbital elements computed from direct $m_0$37-body averaging and then forming $m_0$38, $m_0$39, $m_0$40, $m_0$41, and $m_0$42, the study reports that exactly 27 objects satisfy $m_0$43 (Lei et al., 13 May 2026). Direct $m_0$44-body simulations confirm 26 of these as sustained librators, with the sole exception of S/2019 S1, whose discrepancy is attributed to its proximity to the separatrix, $m_0$45 (Lei et al., 13 May 2026). The same work states that the red contour $m_0$46 marks the separatrix between circulating and librating phases in individual $m_0$47 portraits, and that the extended Brown trajectories “hug” the osculating curves almost perfectly (Lei et al., 13 May 2026).

These applications clarify both the utility and the limitations of the framework. The utility lies in providing a high-accuracy secular model and a simple analytical resonance indicator in regimes where classical averaging is quantitatively inadequate. The limitation is embedded in the assumptions: test-particle dynamics, a fixed outer Kepler perturber, quadrupole truncation of the disturbing function, and neglect of terms beyond $m_0$48 or $m_0$49 (Lei et al., 19 May 2025, Lei et al., 13 May 2026). A common misconception would be to treat the model as a generic three-body theory valid without these restrictions; the published formulation is explicitly tailored to weakly hierarchical, quadrupole-order, test-particle secular dynamics.