Resonance Capture of a Test Particle by an Eccentric Planet in the Presence of Externally-Driven Apsidal Precession

Dong Lai; JT Laune

REVIEW 3 major objections 4 minor 20 references

Resonance Capture of a Test Particle by an Eccentric Planet in the Presence of Externally-Driven Apsidal Precession

T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read External differential apsidal precession can defeat mean-motion resonance capture even at rates far below the resonance-overlap threshold, through two distinct disruption channels.

desk verdict New two-channel picture of MMR capture failure under differential precession is worth refereeing; the low-frequency secular-resonance lobe is plausible but rests on a small-e equilibrium formula used outside its validity range. read the letter →

arxiv 2505.09937 v1 pith:KUOF3DXA submitted 2025-05-15 astro-ph.EP

JT Laune , Dong Lai This is my paper · ORCID

classification astro-ph.EP MSC 70F1537N05

keywords mean-motionresonancecaptureapsidalprecessionrestrictedthree-bodyproblemsecularoverlapKuiperbeltNeptunemigration

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that capture into a first-order $j:j+1$ mean-motion resonance is not determined by the resonance alone: externally driven differential apsidal precession, with the planet and test particle precessing at different rates, can suppress capture even when that precession is far weaker than the usual resonance-overlap threshold. Working in the eccentric restricted three-body problem with a slowly outward-migrating planet, the authors identify two critical precession frequencies at which capture fails: the resonance-overlap frequency $\Delta\omega_{\rm crit}$ of equation (25) and the secular apsidal-resonance frequency $\Delta\omega_{\rm crit,sec}$ of equation (34). The second channel is notable because it destroys capture at planet eccentricities as small as $e_p \sim 0.001$. If correct, the result means resonant trans-Neptunian object capture during Neptune's outward migration depends sensitively on Neptune's eccentricity and on the precession imposed by the inner giant planets.

What carries the argument

The central objects are the two resonance angles $\theta$ (associated with the test particle's eccentricity) and $\theta_p$ (associated with the planet's eccentricity), whose standard reduction to a single mixed angle, valid at $\Delta\omega_{\rm ext}=0$, fails once differential apsidal precession splits them. The load-bearing relations are the resonance-overlap threshold $\Delta\omega_{\rm crit}=\Delta n_c$, the width of the $\theta$-resonance separatrix at first appearance (equations 22--25), and the secular equilibrium eccentricity $e_{\rm eq,sec}$ of equation (33), whose denominator vanishes at $\Delta\omega_{\rm crit,sec}\simeq \frac{1}{4}(j/(j+1))\alpha b^{(1)}_{3/2}(\alpha_0)\mu_p n_p$ (equation 34). The vanishing denominator marks the apsidal precession resonance where the planet's precession frequency matches the test particle's net secular precession, and the divergence is the mechanical reason capture fails at the low-frequency lobe.

What would settle it

Run a full three-body numerical simulation (no secular disturbing-function truncation, no small-e expansion) for the same restricted setup with $\mu_p=5\times10^{-5}$, $e_p=0.001$, and $\Delta\omega_{\rm ext}\simeq\Delta\omega_{\rm crit,sec}$; if the particle is nonetheless captured into the 2:3 resonance, the predicted low-frequency disruption lobe fails, whereas a full simulation that reproduces the eccentricity excitation and escape would confirm it.

Watch

Extended reading notes

Core claim

In the restricted three-body problem of a test particle outside an eccentric planet that migrates outward, the two first-order resonance angles $\theta = (j+1)\lambda - j\lambda_p - \varpi$ and $\theta_p = (j+1)\lambda - j\lambda_p - \varpi_p$ are split by a differential apsidal precession $\Delta\omega_{\rm ext} = \omega_p - \omega$. The paper shows numerically, integrating equations (13)--(16), that capture into the 2:3 resonance fails in two distinct lobes of the $(\Delta\omega_{\rm ext}, e_p)$ plane. At $\Delta\omega_{\rm ext}$ near $\Delta\omega_{\rm crit}$, the widened separatrices of the two sub-resonances overlap and the particle exits through chaotic libration--circulation with low eccentricity. At the much smaller frequency $\Delta\omega_{\rm crit,sec}$, a secular apsidal resonance locks the apsidal angle $\varpi_p - \varpi$ near $\pi$, the particle's eccentricity is excited to roughly 0.23, and capture is skipped even for $e_p$ as small as about 0.001. The comparison with the giant-planet-induced precession rate $\Delta\omega_{\rm GP} \simeq 6.85\times 10^{-5}\, n_p$ suggests these effects are relevant to resonant Kuiper Belt objects captured by migrating Neptune.

Load-bearing premise

The low-frequency disruption channel is predicted from a far-from-resonance equilibrium calculation (equation 33) that assumes the resonant terms average out and that the test particle's eccentricity stays small, while in the simulated disruption trajectory the eccentricity grows to about 0.23 before resonance crossing; if the divergence in equation (33) is an artifact of the truncated secular expansion rather than a real apsidal resonance, that channel needs a different explanation, though the high-frequency overlap channel would remain.

Editorial extensions

If this is right

A finite external apsidal precession must be included in any resonance-capture calculation, because the standard single-angle integrable model becomes invalid at $\Delta\omega_{\rm ext} \neq 0$, with chaos appearing where the sub-resonance separatrices cross.
In the early Solar System, the inner-giant-planet precession $\Delta\omega_{\rm GP}\simeq 6.85\times10^{-5} n_p$ is close enough to $\Delta\omega_{\rm crit,sec}\simeq 1.54\,\mu_p n_p$ that Neptune's eccentricity at the time of migration determines whether the 2:3 resonance captures TNOs or fails.
The secular-resonance failure channel operates at $e_p\sim0.001$, so even a nearly circular migrating planet can fail to capture if the differential precession matches the secular resonance frequency.
Because successful capture requires either very small eccentricity or precession rates far away from the two critical values, the resonant-to-nonresonant ratio of the Kuiper Belt is not simply a migration-smoothness diagnostic; it also encodes differential precession and Neptune's eccentricity.
Faster migration makes capture more robust near the critical frequencies, since less time is spent in the chaotic or secular-excitation zones, so the migration timescale and the disruption channels trade off.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The secular-resonance disruption channel gives a natural population pathway from the Kuiper Belt to the scattered disk that does not require close encounters with Neptune itself: failed capture near $\Delta\omega_{\rm crit,sec}$ pumps eccentricity to orbit-crossing values, and the particle is then ejected or scattered by the planet on a subsequent pass; the paper raises this possibility but does n
Because $\Delta\omega_{\rm crit,sec}$ scales linearly with $\mu_p n_p$, equivalent low-frequency disruption lobes should appear for more massive planets and at other $j:j+1$ resonances, shifted to correspondingly larger precession rates; this is a testable prediction of the equilibrium-divergence argument.
The boundary between capture and disruption near $\Delta\omega_{\rm crit,sec}$ is chaotic, so the observable prediction is not a sharp line but a capture probability; an ensemble of integrations with random initial $\varpi$ at fixed $\Delta\omega_{\rm ext}$ and $e_p$ would map this probability as a function of migration speed.
If the small-e expansion in equation (33) is improved to finite eccentricity, the apparent divergence at $\Delta\omega_{\rm crit,sec}$ may broaden into a finite-width resonance; one test is to compare the location of the failure lobe in full N-body integrations with the analytic frequency.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Briefly, this paper gives a restricted three-body result worth knowing: when a migrating eccentric planet's resonance sweeps past a test particle, externally driven differential apsidal precession can suppress capture at two distinct frequencies. The high-frequency channel is ordinary resonance overlap between the θ and θp sub-resonances; the threshold (Eq 25) matches the captured/disrupted boundary in the numerics and the chaotic signature in the FLI maps. The lower-frequency channel is the novel piece: a secular apsidal resonance (Eq 34) that disrupts capture even for planet eccentricities down toward 1e-3. The diagnostic runs in Fig 7 do a lot of work—turning off the secular forcing kills the low-frequency spike, turning off the θp resonant term does not. That is good experimental hygiene.

The main soft spot is the analytic support for the secular channel. Eq 33 is a far-from-resonance equilibrium in which resonant terms are averaged away, and it diverges at Eq 34. But the representative disrupting trajectory in Fig 5 shows e climbing to about 0.23 while n/np is still approaching the 2:3 value. That is exactly where O(e^2) secular theory and the averaging assumption are unreliable. The paper itself notes the equilibrium e at the resonance location is about 0.09 at those parameters—a factor of about 2.5 below the simulated peak. So the divergence in Eq 33 could be an artifact of the truncation rather than a real dynamical singularity. That does not kill the claim—the numerical lobe is there and the diagnostic runs support a secular effect—but it means the analytic threshold is not yet established. A direct test, such as a full N-body integration or a non-perturbative treatment, would settle it.

Two smaller points. Each point in the capture maps is a single integration; in the chaotic boundary regions the lobe edges could be fuzzy. The authors show two trajectories at one boundary point, which is suggestive but not statistics. No code or data is shipped, which makes it harder to probe the suspicious region. The Neptune/TNO discussion is appropriately hedged and does not oversell.

Overall: this is a solid, clearly written extension of El Moutamid and Murray et al. The high-frequency lobe is on firm ground; the low-frequency lobe is plausible and important but needs the analytic rough edge cleaned up. I would send it to a serious referee, and ask for a targeted test of the secular-resonance threshold plus some release of artifacts.

Referee Report

3 major / 4 minor

Summary. The paper studies capture of a massless test particle into first-order j:j+1 mean-motion resonances with an eccentric, outwardly migrating planet, in the presence of externally imposed differential apsidal precession Δωext. Using a resonant Hamiltonian, FLI maps, and direct integrations of a four-variable ODE system, the authors identify two critical precession rates that disrupt capture: resonance overlap at Δωcrit (Eq. 25) and a secular apsidal resonance at Δωcrit,sec (Eq. 34). Capture maps in the ep–Δωext plane show disruption lobes near both frequencies, and the authors argue that this is relevant to resonance capture of trans-Neptunian objects by migrating Neptune. The central claim is that even a differential precession rate far below the resonance-overlap criterion can strongly modify capture probability, and that the low-frequency channel survives down to planet eccentricities ep ≈ 0.001.

Significance. If the low-frequency disruption channel survives scrutiny, the result is significant: it predicts a non-monotonic dependence of capture probability on differential apsidal precession and planet eccentricity, with parameter-free analytic thresholds that can be compared with N-body simulations of Neptune migration. The paper has clear strengths: the critical frequencies are analytic predictions rather than fitted constants; the capture maps are direct numerical integrations that could in principle disagree with those predictions; the FLI maps and diagnostic runs (turning off the secular forcing, or setting βp=0) isolate the two proposed channels; and the appendices supply the underlying derivations. The paper is well organized and the conclusions are stated precisely enough to be tested. The main weakness is that the low-frequency secular-resonance channel rests on a truncated, far-from-resonance equilibrium calculation applied in a regime where the simulated trajectory reaches eccentricities where that calculation is quantitatively unreliable.

major comments (3)

[§4.3, Eq. (33)–(34), Figure 5] The low-frequency disruption threshold is derived by assuming the system is far from resonance, so the resonant terms average out and ṅ=0 and ė=0, and by truncating the secular disturbing function at O(e²). The representative disrupting trajectory in Figure 5 reaches e≈0.23 while n/np approaches the 2:3 resonant value, which is exactly the regime where neither assumption is safe. At the parameters of Figure 5 (μp=5×10⁻⁵, ep=0.03, Δωext=10⁻⁴np), Eq. (33) gives e_eq≈0.09 at the resonance location, roughly a factor of 2.5 below the simulated peak eccentricity. Because Eq. (34) and the blue line in Figure 6 are the only quantitative prediction tying the low-frequency lobe to the secular resonance, I ask for a direct test that this divergence is not an artifact of the truncated expansion. For example, the authors could recompute the capture map for a range of μp and verify that the lobe center follows the predicted 1.54 μp np scaling, or replace the O(e²) equilibrium estimate with the full (untruncated) secular disturbing function and show that the divergence at Δωcrit,sec persists. Without such a test, the existence of the low-frequency disruption channel remains insufficiently supported, even though the high-frequency resonance-overlap lobe at Δωcrit is well supported.
[Figures 6 and 7] Each cell in the capture maps appears to represent a single integration, with no repeated runs or uncertainty quantification. The text itself describes the capture/disruption boundary as fuzzy and initial-condition sensitive (§5.2), and Figure 2 shows a chaotic resonant zone with a mixture of regular islands and chaotic regions. With one trajectory per parameter cell, the plotted lobe boundaries, and especially the claim that the low-frequency lobe reaches ep≈0.001, are not yet robust. I recommend repeating each integration over a modest ensemble of initial ϖ and θp values near the two critical frequencies and reporting the capture probability per cell, or at least stating explicitly how many initial conditions per cell were used and how the boundary depends on integration time and error tolerance.
[§4.1 and §6] There is an internal inconsistency about the eccentricity range over which capture is disrupted. Section 4.1 states that 'capture is secure for very small ep's, regardless of the Δωext values,' while the abstract, Section 5, and Section 6 claim that MMR capture can fail 'even for very small values of ep≈0.001.' These statements cannot both be correct as written. The authors should define 'very small' quantitatively and reconcile the two statements with the actual content of Figure 6, where the low-frequency lobe apparently does extend to small ep. This is a central claim of the paper and needs to be stated unambiguously.

minor comments (4)

[§2.2, Eq. (17)] The derivation of ω/ωp≃α₀^{7/2} would be easier to follow if the steps n/np≈α^{3/2} and (ap/a)²=α² were written out before the final expression.
[§3.2, Figure 1] In the third row, the description of the pink trajectory as 'stretched' is not immediately visible in the printed figure; marking the relevant trajectory with an arrow or label would help the reader connect the text to the panel.
[§4.1, Figure 6 caption] The caption does not state how many initial conditions were integrated per parameter cell or whether the random choices of initial ϖ and θp differ from cell to cell; please add this information to the caption or to §2.1.
[§5.4] The comparison with Murray et al. (2022) is qualitative; a brief statement about how much eccentricity damping is needed to suppress the two disruption channels would help the reader judge the astrophysical relevance for comparable-mass planets.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the two critical precession frequencies are parameter-free analytic predictions, and the capture maps are independent numerical integrations that could falsify them.

full rationale

The paper's central claims are the two disruption channels at Delta_omega_crit (Eq. 25) and Delta_omega_crit,sec (Eq. 34). Both are obtained analytically: Eq. (25) follows from the pendulum/Andoyer resonance width (Eqs. 22-24) and Eq. (34) from setting the secular apsidal precession equilibrium denominator to zero in Eq. (33). Neither formula contains parameters fitted to the capture-outcome maps; the maps in Figures 6-7 are produced by direct numerical integration of equations (13)-(16) over a grid of initial conditions, and the green/blue lines are overlaid after the fact as comparisons. The paper itself states the comparison rather than the fit: 'We see in Figure 6 and the middle panel of Figure 7 that this value fits the spike of disruption well.' A numerical experiment can disagree with the analytic threshold, and the authors' own Figure 5 shows e growing to about 0.23, outside the small-e regime, so the secular formula is not tautologically enforced by the integration. The self-citations (Laune et al. 2022) are used for qualitative comparison of apsidal alignment behavior and are not load-bearing; the integrable Delta_omega_ext=0 reduction is re-derived in Appendix C, not imported purely by citation. The only mild concern is that the analytic secular threshold and the equations of motion share the same truncated disturbing function (12), so a systematic error in that truncation would affect both; this is a model-validity/robustness caveat, not a circular reduction of the kind where the prediction is defined as the input or fitted to the output. Accordingly, no circular step meets the evidentiary bar of the review instructions.

Assumptions & free parameters 2 free parameters · 7 assumptions · 0 invented entities

The model uses no invented physical entities. No parameter is fitted to the capture maps; the critical frequencies are analytic predictions. The main hand-chosen inputs are the migration timescale and initial test-particle eccentricity, plus standard constants evaluated at resonance. The heaviest unverified input is the small-e secular equilibrium used to predict the low-frequency disruption lobe.

free parameters (2)

migration timescale tau_m = 10^6 P_p,0 (main runs), 10^5 P_p,0 (Figure 7 bottom)
Chosen by hand to set slow, adiabatic migration. Section 5.1 shows capture outcomes depend on it near the critical lines, so it is a free parameter of the model.
test particle initial eccentricity e0 = 0.001
Chosen initial condition for all capture runs; the maps are not averaged over e0. The claim about very small planet eccentricity is made for this low test-particle e.

assumptions (7)

domain assumption The restricted three-body Hamiltonian (7) truncated to first-order resonant terms and the secular disturbing function (12) captures the capture dynamics.
This truncation defines the model; higher-order terms and finite test-particle mass are neglected. Invoked in Section 2.
domain assumption External apsidal precession rates omega and omega_p are constant, with companion perihelion angle equal to omega_p t.
Section 2 states this for simplicity; real precession varies during migration and disk dispersal, which the parameter survey does not cover.
domain assumption Laplace coefficients beta, beta_p, b3_2^(1), and b3_2^(2) are evaluated at alpha_0 and held constant.
Used in equations (8)-(9) and (12) and in the secular resonance threshold (34); the paper notes this simplification but does not test sensitivity.
domain assumption Migration is slow and adiabatic: tau_m n_p is much larger than 1.
Section 2 sets tau_m = 10^6 P_p,0; this assumption underlies the capture-process interpretation.
domain assumption Far from resonance, resonant terms average out, so the secular equilibrium (33) and divergence (34) describe the pre-resonance eccentricity excitation.
Section 4.3 uses this to derive the secular-resonance channel; the assumption of small e is violated in Figure 5 (e about 0.23).
domain assumption The external precession follows a quadrupole interior-potential scaling omega/omega_p approximately alpha_0^(7/2) (equation 17).
Used to set omega and Delta_omega_ext in the capture maps in Section 4.1; valid for interior giant planets, not for general external sources.
standard math Hamiltonian perturbation theory, canonical transformations, Lagrange planetary equations, and the FLI variational method are valid as used.
Background tools used throughout; no machine-checked verification.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Resonance Capture of a Test Particle by an Eccentric Planet in the Presence of Externally-Driven Apsidal Precession." pith.science (2026). https://pith.science/paper/KUOF3DXA

@misc{pith2026250509937,
  author       = {Pith},
  title        = {Pith review of: Resonance Capture of a Test Particle by an Eccentric Planet in the Presence of Externally-Driven Apsidal Precession},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KUOF3DXA}},
  note         = {Machine review of arXiv:2505.09937}
}

abstract

Planets undergoing convergent migration can be captured into mean-motion resonance (MMR), in which the planets' periods are related by integer ratios. The dynamics of MMR are typically considered in isolation, including only the forces between the planets and with the central star. However, the planets are often subjected to external forces that induce apsidal precessions, which may split the MMR into two sub-resonances and give rise to chaotic motion due to resonance overlap. In this study, we investigate how such externally induced differential apsidal precession affects capture into first-order $j:j+1$ MMRs. We study the restricted three-body problem for a test particle outside of an eccentric planet, with the planet undergoing outward migration. We find that capture can be sensitive to the differential apsidal precession, even when the precession rate is much smaller than the resonance overlap criterion. We identify two critical precession frequencies -- related to resonance overlap and secular apsidal resonance -- around which capture is disrupted even for very small planet eccentricity. Our results may help clarify the capture process of resonant trans-Neptunian objects by outward-migrating Neptune in early Solar System.

Figures

Figures reproduced from arXiv: 2505.09937 by the authors.

**Figure 1.** Phase-space trajectories of the 2:3 resonant system at varying degrees of differential apsidal precession, ∆ωext, for µp = 5×10−5 and ep = 0.03. The secular forcing from mp on the test mass is neglected. Each color-coded trajectory corresponds to a different initial SMA. The left panels show θ vs. a; the right panels show θp vs. a. The system is initialized with k = −0.1, which sets the initial e via equation (20). … view at source ↗

**Figure 2.** FLI (Section 2.3) as a function of the initial a and θ for the resonant system (without the secular forcing from mp and planet migration) for µp = 5 × 10−5 and ep = 0.03. The value of k = −0.1 is held fixed and e is chosen via equation (20). We initialize θp = θ + π. Left: ∆ωext = 0.1np; the separatrices are clearly visible and well-separated. Right: ∆ωext = 0.03np; the separatrices overlap and widespread chaos perm… view at source ↗

**Figure 3.** Capture of the test particle into a 2:3 MMR with a planet of mass µp = 5 × 10−5 and ep = 0.03, under a weak differential apsidal precession (∆ωext = 10−5np,0). The initial conditions of the test particle are θp = 0, e = 0.001, ϖ = π and a/ap ≃ 1.48 (n/np,0 = 1/1.8). As the planet migrates outward, n/np increases and approaches 2/3. When k (equation 20) crosses ksep (equation 22, red line), θ begins to librate and e … view at source ↗

Figures from the paper (4 more)

**Figure 4.** Figure 4: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]

**Figure 5.** Figure 5: Same as [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]

**Figure 6.** Figure 6: Capture outcomes for the 2 : 3 MMR as a function of differential apsidal precession rate ∆ωext and planet eccentricity ep. Yellow dots indicate disruption; purple dots indicate successful capture. The resonance overlap threshold ∆ωcrit (equation 25) and secular resonan…

**Figure 7.** Figure 7: Same as [PITH_FULL_IMAGE:figures/full_fig_p012_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 19 canonical work pages

[1]

Batygin, K., Ngo, H., and Meshkat, T. (2019). An Excess of Jupiter Analogs in Super-Earth Systems. The Astronomical Journal, 157(2):52

work page 2019
[2]

M., Payne, M., and Holman, M

Deck, K. M., Payne, M., and Holman, M. J. (2013). First-order resonance overlap and the stability of close two-planet systems. The Astrophysical Journal , 774(2):129. El Moutamid, M., Sicardy, B., and Renner, S. (2014). Coupling between corotation and Lindblad resonances in the presence of secular precession rates. Celestial Mechanics and Dynamical Astron...

work page 2013
[3]

Fernandez, J. A. and Ip, W.-H. (1984). Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals. Icar, 58(1):109–120

work page 1984
[4]

Gomes, R. S. (2000). Planetary Migration and Plutino Orbital Inclinations. The Astronomical Journal , 120:2695–2707

work page 2000
[5]

Gomes, R. S. (2003). The origin of the Kuiper Belt high-inclination population. Icarus, 161(2):404–418

work page 2003
[6]

Hahn, J. M. and Malhotra, R. (2005). Neptune’s Migration into a Stirred-Up Kuiper Belt: A Detailed Comparison of Simulations to Observations. AJ, 130(5):2392–2414

work page 2005
[7]

and Lemaitre, A

Henrard, J. and Lemaitre, A. (1983). A second fundamental model for resonance. Celestial Mechanics, 30(2):197–218

work page 1983
[8]

Henrard, J., Lemaitre, A., Milani, A., and Murray, C. D. (1986). The reducing transformation and Apocentric Librators. Celestial mechanics, 38(4):335–344

work page 1986

Show all 20 references

[9]

T., Rodet, L., and Lai, D

Laune, J. T., Rodet, L., and Lai, D. (2022). Apsidal alignment and anti-alignment of planets in mean-motion resonance: Disc-driven migration and eccentricity driving. Monthly Notices of the Royal Astronomical Society, 517(3):4472–4488

2022
[10]

Malhotra, R. (1993). The origin of Pluto’s peculiar orbit. Nature, 365(6449):819–821

1993
[11]

Malhotra, R. (1995). The Origin of Pluto’s Orbit: Implications for the Solar System Beyond Neptune. AJ, 110:420

1995
[12]

Murray, C. D. and Dermott, S. F. (2000). Solar System Dynamics. Cambridge University Press, Cambridge

2000
[13]

Murray, Z., Hadden, S., and Holman, M. J. (2022). The Effects of Disk Induced Apsidal Precession on Planets Captured into Mean Motion Resonance. The Astrophysical Journal, 931(1):66. 14

2022
[14]

Murray-Clay, R. A. and Schlichting, H. E. (2011). USING KUIPER BELT BINARIES TO CONSTRAIN NEPTUNE’S MIGRATION HISTORY. The Astrophysical Journal, 730(2):132. Nesvorn´ y, D. (2015). Evidence for Slow Migration of Neptune from the Inclination Distribution of Kuiper Belt Objects....

2011
[15]

Peale, S. J. (1976). Orbital resonance in the solar system. ARA&A, 14:215–246

1976
[16]

C., Petigura, E

Petit, A. C., Petigura, E. A., Davies, M. B., and Johansen, A. (2020). Resonance in the K2-19 system is at odds with its high reported eccentricities. arXiv:2003.04931 [astro-ph]

2020 arXiv
[17]

A., and Laskar, J., editors (2016)

Skokos, C., Gottwald, G. A., and Laskar, J., editors (2016). Chaos Detection and Predictability , volume 915 of Lecture Notes in Physics . Springer, Berlin, Heidelberg. Van Zandt, J., Petigura, E. A., Lubin, J., Weiss, L. M.,

2016
[18]

W., and Zink, J

Tyler, D., Yee, S. W., and Zink, J. (2025). The TESS–Keck Survey. XXIV. Outer Giants May Be More Prevalent in the Presence of Inner Small Planets. AJ, 169(5):235

2025
[19]

Wisdom, J. (1986). Canonical solution of the two critical argument problem. Celestial Mechanics, 38:175–180

1986
[20]

old” canonical variables are related to the “new

Zhu, W. and Wu, Y. (2018). The Super Earth-Cold Jupiter Relations. The Astronomical Journal , 156:92. 15 APPENDIX A. VARIATIONAL EQUATIONS FOR THE FAST LYAPUNOV INDICATOR METHOD We scale the problem using np and ap, and we denote npt = τ. We follow the method in Skokos et al. ...

2018

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.

[1] [1]

Batygin, K., Ngo, H., and Meshkat, T. (2019). An Excess of Jupiter Analogs in Super-Earth Systems. The Astronomical Journal, 157(2):52

work page 2019

[2] [2]

M., Payne, M., and Holman, M

Deck, K. M., Payne, M., and Holman, M. J. (2013). First-order resonance overlap and the stability of close two-planet systems. The Astrophysical Journal , 774(2):129. El Moutamid, M., Sicardy, B., and Renner, S. (2014). Coupling between corotation and Lindblad resonances in the presence of secular precession rates. Celestial Mechanics and Dynamical Astron...

work page 2013

[3] [3]

Fernandez, J. A. and Ip, W.-H. (1984). Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals. Icar, 58(1):109–120

work page 1984

[4] [4]

Gomes, R. S. (2000). Planetary Migration and Plutino Orbital Inclinations. The Astronomical Journal , 120:2695–2707

work page 2000

[5] [5]

Gomes, R. S. (2003). The origin of the Kuiper Belt high-inclination population. Icarus, 161(2):404–418

work page 2003

[6] [6]

Hahn, J. M. and Malhotra, R. (2005). Neptune’s Migration into a Stirred-Up Kuiper Belt: A Detailed Comparison of Simulations to Observations. AJ, 130(5):2392–2414

work page 2005

[7] [7]

and Lemaitre, A

Henrard, J. and Lemaitre, A. (1983). A second fundamental model for resonance. Celestial Mechanics, 30(2):197–218

work page 1983

[8] [8]

Henrard, J., Lemaitre, A., Milani, A., and Murray, C. D. (1986). The reducing transformation and Apocentric Librators. Celestial mechanics, 38(4):335–344

work page 1986