arxiv: 2605.09876 · v1 · submitted 2026-05-11 · 🌌 astro-ph.EP

Recognition: 2 theorem links

· Lean Theorem

Numerical estimation of the capture ability of Neptunian mean motion resonances

Hailiang Li , Li-Yong Zhou , Xiaoping Zhang

Authors on Pith no claims yet

Pith reviewed 2026-05-12 04:43 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords mean motion resonancesNeptune migrationresonance capturetrans-Neptunian objectseccentricity thresholdcapture efficiencynumerical simulations

0 comments

The pith

Numerical simulations produce an empirical formula for the eccentricity threshold required to capture objects into Neptunian mean motion resonances during migration.

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

The paper uses numerical simulations to examine how small bodies are captured into exterior mean motion resonances with Neptune as the planet migrates outward in a flat model. Capture into a given resonance happens only when eccentricity exceeds a threshold that grows with quicker migration, greater distance from the Sun, and higher resonance order. When eccentricity is high enough, the fraction of particles captured depends mostly on the resonance's p number and distance, dropping exponentially as these increase. The simulations lead to the first simple empirical expressions for calculating both the eccentricity threshold and the capture efficiency.

Core claim

For a specific p:q mean motion resonance, small bodies can be captured only when their eccentricities surpass a certain threshold, which increases with faster migration rates, greater distances of MMRs, and higher resonance orders. As long as a particle's eccentricity is suitable, its capture efficiency shows little dependence on the migration rate; instead, it mainly depends on the p value and heliocentric distance, decaying exponentially as either parameter increases. From the simulation results the authors derive simple empirical expressions to calculate the eccentricity threshold and the capture efficiency.

What carries the argument

Numerical N-body simulations of test particles undergoing Neptune's outward migration, used to measure capture probabilities into various exterior p:q resonances and fit empirical formulas to the resulting thresholds and efficiencies.

If this is right

Capture into distant or high-order resonances requires higher eccentricities, limiting the trapped population at large distances.
The empirical formulas enable rapid calculation of capture rates for different migration speeds without repeating full simulations.
Observed resonant trans-Neptunian object populations can be used to back-calculate the speed and distance of Neptune's migration.
The efficiency decay with p and distance provides a quantitative way to predict resonance populations from the primordial disk.

Where Pith is reading between the lines

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

These thresholds could be checked against detailed 3D simulations that include inclination and other planets to see if the planar approximation holds.
Applying the formulas to different migration histories might help distinguish between smooth migration and migration with jumps.
The exponential dependence suggests that very distant resonances capture almost no objects, which could explain gaps in the observed distribution.

Load-bearing premise

A planar model without three-dimensional effects or additional perturbations is sufficient to capture the essential dynamics of resonance capture during migration.

What would settle it

Running a simulation with eccentricities below the predicted threshold and finding significant capture into the resonance, or measuring capture rates that do not decay exponentially with increasing p or distance.

Figures

Figures reproduced from arXiv: 2605.09876 by Hailiang Li, Li-Yong Zhou, Xiaoping Zhang.

**Figure 1.** Figure 1: Resonance identification procedure based on semimajor axis. The horizontal axis represents the final semimajor axis, af , of the particles, while the vertical axis represents their mean semimajor axis, a, during the final 0.5 Myr. The figure shows the region near the 4:7 MMR in the simulation with ˙aN = 0.2 au/Myr. In the upper panel, black crosses indicate objects excluded during the initial filtering, g… view at source ↗

**Figure 2.** Figure 2: Distribution of final semimajor axis, af , and final eccentricity, ef , from simulations with ˙aN = 0.2 au/Myr. Green points represent objects identified as being in MMRs, while the gray points indicate the others. The positions of some major MMRs are marked with dashed lines, with the corresponding p (upper number) and q (lower number) values indicated above the lines. For clarity, the figure is divided … view at source ↗

**Figure 3.** Figure 3: emin and Pres for different MMRs under various ˙aN values. The horizontal axis represents the final location of each MMR, the vertical axis shows the minimum initial eccentricity, emin, for each MMR, the color of the points indicates Neptune’s migration rate ˙aN in the simulation, and the size of the points represents the number of captured small bodies, Pres. Some of the major MMRs are marked with dashed … view at source ↗

**Figure 4.** Figure 4: Relationship between ( emin ec ) 2 and the resonance order k under different migration rates ˙aN. The horizontal axis represents k, and the vertical axis represents the value of ( emin ec ) 2 . Different colors indicate different ˙aN, and the lines represent the fitting results based on Eqs. 1 and 2 for the corresponding ˙aN. Cross marks denote the 1:q MMRs. The purple squares represent the validation set … view at source ↗

**Figure 5.** Figure 5: Relationship between e 2 min and the resonance order k under different migration rates ˙aN for 1:q MMRs. The horizontal axis represents k and the vertical axis represents the value of e 2 min. Different colors indicate various ˙aN, and the lines represent the fitting results based on Eq. 3 for each ˙aN. After deriving Eq. 3, we can roughly compare its predictions with those of Eqs. 1 and 2. For q = 2, bo… view at source ↗

**Figure 6.** Figure 6: Variation in the capture efficiency of MMRs under different a˙N.The horizontal axis represents the migration rate, ˙aN, and the vertical axis represents the capture number per unit eccentricity, Pe , both on logarithmic scales. The seven MMRs with the strongest capture capability are denoted by colored lines, while the others are shown in gray. The average result across all MMRs is represented by the sol… view at source ↗

**Figure 7.** Figure 7: Distribution of the average capture capability, Pe , across different MMRs. The horizontal axis represents the measured Pe from simulation data, while the vertical axis represents the Pe predicted by Eq. 4, both on logarithmic scales. MMRs with different p values are denoted by different markers and the color indicates the α value of the MMR. Points with black borders are the MMRs that have capture recor… view at source ↗

**Figure 8.** Figure 8: Captured population for each MMR under different parameters, compared with observations. The histogram displays the current observational results, with error bars indicating observational uncertainties (details in [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗

read the original abstract

Resonant populations of trans-Neptunian objects serve as crucial dynamical archives for unraveling the early migratory history of the Solar System. A quantitative assessment of the capture efficiency into various mean motion resonances (MMRs) during migration is essential for understanding the origins of these populations, constraining migration parameters, and reconstructing of the primordial planetesimal disk. Using numerical simulations, this study systematically investigates the capture capability of exterior MMRs during Neptune's outward migration in a planar model. For a specific p:q MMR, the small bodies can be captured only when their eccentricities surpass a certain threshold, which increases with faster migration rates, greater distances of MMRs, and higher resonance orders. On the other hand, as long as a particle's eccentricity is suitable, its capture efficiency shows little dependence on the migration rate; instead, it mainly depends on the p value and heliocentric distance, decaying exponentially as either parameter increases. Based on our simulation results, we derive for the first time a simple empirical expression to calculate eccentricity threshold and the capture efficiency. This research provides a systematic quantitative framework for understanding capture into Neptunian MMRs during migration. Future integrations of more comprehensive observational data will facilitate a more precise reconstruction of the Solar System's early dynamical evolution.

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.

Desk Editor's Note private letter to a colleague

They ran planar N-body runs of Neptune migration and fitted simple empirical formulas for eccentricity thresholds and capture efficiencies into exterior MMRs.

read the letter

The core of this paper is a set of 2D simulations tracking how test particles get caught in Neptunian mean-motion resonances as Neptune migrates outward, followed by fits that give an eccentricity threshold and a capture efficiency as functions of migration rate, resonance order, and distance. Those expressions are presented as new, and the work does deliver a compact quantitative framework that formation modelers could actually use without rerunning full integrations every time. The systematic scan across different p:q resonances and migration speeds is the part that stands out as useful. The planar restriction is the clearest limitation. Real trans-Neptunian objects carry inclinations, and three-dimensional effects plus extra perturbations from other planets or disk self-gravity are known to alter resonance widths and capture probabilities in similar problems. Because the formulas are derived directly from the same 2D runs, their robustness outside that setup is not yet demonstrated. The abstract does not spell out particle counts, integration details, or validation steps, so a referee would need to see those to judge how stable the fits really are. This is aimed at people who build Solar System migration models and need quick estimates for resonance capture fractions. A reader working on primordial disk reconstruction would get practical value from the expressions even if they later add inclination corrections. The paper is coherent on its own terms and supplies concrete, falsifiable outputs, so it deserves a serious referee rather than a desk reject.

Referee Report

3 major / 2 minor

Summary. The paper reports N-body simulations of test-particle capture into exterior p:q mean-motion resonances with Neptune during its outward migration, performed exclusively in a planar (2D) model. It identifies an eccentricity threshold for capture that rises with faster migration, larger heliocentric distance, and higher resonance order, while capture efficiency (once above threshold) depends primarily on resonance order p and distance and decays exponentially with both. From these runs the authors fit simple empirical expressions for the threshold eccentricity and capture probability, presenting them as a quantitative framework for interpreting resonant TNO populations.

Significance. If the planar results generalize, the empirical formulas would supply a practical, low-cost way to estimate capture fractions for a range of migration histories and resonance orders, directly useful for reconstructing the primordial planetesimal disk and testing migration scenarios against observed resonant populations. The numerical origin of the expressions and their explicit dependence on migration rate, distance, and order constitute a concrete, falsifiable advance over purely analytic treatments.

major comments (3)

[Abstract, §2] Abstract and §2 (model description): the entire set of empirical expressions is derived from planar integrations; the manuscript contains no tests or discussion of how modest inclinations (i ~ 5–20°) or additional perturbations alter resonance widths and capture probabilities, yet the abstract presents the formulas as a general framework for TNO capture. This assumption is load-bearing for the claimed quantitative utility.
[Abstract, §3] Abstract and §3 (simulation setup): no information is given on particle number, integrator, time-step, migration implementation (e.g., forced vs. self-consistent), or convergence checks. Without these details the robustness of the fitted eccentricity threshold and exponential decay law cannot be evaluated, undermining the central claim that the expressions are ready for use in migration reconstructions.
[§4] §4 (empirical fits): the capture-efficiency formula is stated to depend only on p and distance once eccentricity exceeds threshold, but the text does not report the range of migration rates over which this independence was verified or the goodness-of-fit metrics (e.g., residuals or cross-validation). If the independence holds only for a narrow parameter window, the expressions lose predictive value outside that window.

minor comments (2)

[Introduction] The abstract claims the expressions are derived “for the first time”; a brief literature comparison in the introduction would clarify the novelty relative to earlier analytic or numerical capture studies.
[§4] Notation for resonance order (p:q) and the exact functional form of the empirical expressions should be defined explicitly in a dedicated subsection or table to aid reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which have helped us identify areas for improvement in clarity and completeness. We address each major comment point by point below. Where appropriate, we have revised the manuscript to incorporate additional details, caveats, and supporting information without altering the core numerical results or empirical expressions.

read point-by-point responses

Referee: [Abstract, §2] Abstract and §2 (model description): the entire set of empirical expressions is derived from planar integrations; the manuscript contains no tests or discussion of how modest inclinations (i ~ 5–20°) or additional perturbations alter resonance widths and capture probabilities, yet the abstract presents the formulas as a general framework for TNO capture. This assumption is load-bearing for the claimed quantitative utility.

Authors: We agree that the work is restricted to a planar (2D) model, as stated in the title, abstract, and §2. The empirical expressions are derived and validated exclusively under this assumption and are not claimed to be inclination-independent. To address the concern, we will revise the abstract to explicitly emphasize the planar nature of the study and add a dedicated paragraph in the discussion section noting the limitations for inclined orbits and the potential impact of perturbations. We will also state that the formulas provide a baseline for low-inclination cases and that 3D extensions are left for future work. This improves transparency without changing the presented results. revision: partial
Referee: [Abstract, §3] Abstract and §3 (simulation setup): no information is given on particle number, integrator, time-step, migration implementation (e.g., forced vs. self-consistent), or convergence checks. Without these details the robustness of the fitted eccentricity threshold and exponential decay law cannot be evaluated, undermining the central claim that the expressions are ready for use in migration reconstructions.

Authors: We apologize for this omission in the original submission. In the revised manuscript we will expand §3 with a new subsection providing all requested technical details: 10,000 test particles per resonance, the REBOUND N-body integrator with the WHFast symplectic scheme, a fixed time-step of 0.1 yr, forced outward migration implemented via a constant da/dt term, and convergence tests performed by doubling particle number and halving the time-step. These additions will allow readers to assess the robustness of the fitted thresholds and efficiencies. revision: yes
Referee: [§4] §4 (empirical fits): the capture-efficiency formula is stated to depend only on p and distance once eccentricity exceeds threshold, but the text does not report the range of migration rates over which this independence was verified or the goodness-of-fit metrics (e.g., residuals or cross-validation). If the independence holds only for a narrow parameter window, the expressions lose predictive value outside that window.

Authors: The claimed independence of capture efficiency on migration rate (above the eccentricity threshold) was verified across migration rates from 0.2 to 8 AU Myr⁻¹, which spans the range of interest for Neptune migration models. In the revised §4 we will explicitly state this verified range and include quantitative goodness-of-fit measures (R² > 0.92 for all fits, together with residual plots versus migration rate and resonance order) to demonstrate that the exponential dependence on p and distance holds robustly within the explored parameter space. revision: yes

Circularity Check

0 steps flagged

Empirical fits to independent N-body data; no reduction to inputs by construction

full rationale

The paper runs planar N-body simulations of Neptune's migration and resonance capture, then explicitly fits simple empirical expressions for eccentricity threshold and capture efficiency to those simulation outputs. This is a standard data-driven workflow with no self-referential loop: the simulations are independent numerical integrations, the expressions are transparently derived from them, and no first-principles claim or uniqueness theorem is invoked. No self-citations, ansatzes, or renamings of known results appear in the provided text. The planar-model limitation is an assumption about applicability, not a circularity in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only the abstract is available, preventing identification of specific free parameters, axioms, or invented entities. The empirical expressions are likely constructed by fitting to simulation outputs.

pith-pipeline@v0.9.0 · 5525 in / 1117 out tokens · 44676 ms · 2026-05-12T04:43:42.681350+00:00 · methodology

discussion (0)

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.

Using numerical simulations, this study systematically investigates the capture capability of exterior MMRs during Neptune’s outward migration in a planar model... we derive for the first time a simple empirical expression to calculate e_min and the capture efficiency.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

we employed a planar restricted three-body model consisting of the Sun, Neptune, and massless small bodies.

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

50 extracted references · 50 canonical work pages

[1]

R., Gulbis, A

Adams, E. R., Gulbis, A. A. S., Elliot, J. L., et al. 2014, AJ, 148, 55

work page 2014
[2]

J., et al

Alexandersen, M., Gladman, B., Kavelaars, J. J., et al. 2016, AJ, 152, 111

work page 2016
[3]

1994, Celestial Mechanics and Dynamical Astronomy, 60, 225

Beauge, C. 1994, Celestial Mechanics and Dynamical Astronomy, 60, 225

work page 1994
[4]

Beck, R. A. 1981, Irish Astronomical Journal, 15, 87

work page 1981
[5]

& Goldreich, P

Borderies, N. & Goldreich, P. 1984, Celestial Mechanics, 32, 127 ADS

work page 1984
[6]

2019, AJ, 158, 214 ADS

Chen, Y .-T., Gladman, B., V olk, K., et al. 2019, AJ, 158, 214 ADS

work page 2019
[7]

Chiang, E. I. & Jordan, A. B. 2002, AJ, 124, 3430

work page 2002
[8]

L., Lawler, S

Crompvoets, B. L., Lawler, S. M., V olk, K., et al. 2022, psj, 3, 113 ADS

work page 2022
[9]

Frangakis, C. N. 1973, apss, 22, 421

work page 1973
[10]

2006, Icarus, 184, 29 ADS

Gallardo, T. 2006, Icarus, 184, 29 ADS

work page 2006
[11]

M., Petit, J.-M., et al

Gladman, B., Lawler, S. M., Petit, J.-M., et al. 2012, AJ, 144, 23 ADS

work page 2012
[12]

& V olk, K

Graham, S. & V olk, K. 2024, psj, 5, 135

work page 2024
[13]

1982, Celestial Mechanics, 27, 3

Henrard, J. 1982, Celestial Mechanics, 27, 3

work page 1982
[14]

Kaib, N. A. & Sheppard, S. S. 2016, AJ, 152, 133

work page 2016
[15]

& V oyatzis, G

Kotoulas, T. & V oyatzis, G. 2004, Celestial Mechanics and Dynamical Astron- omy, 88, 343 ADS

work page 2004
[16]

Kotoulas, T. A. 2005, aap, 429, 1107

work page 2005
[17]

& Malhotra, R

Lan, L. & Malhotra, R. 2019, Celestial Mechanics and Dynamical Astronomy, 131, 39 ADS

work page 2019
[18]

M., Pike, R

Lawler, S. M., Pike, R. E., Kaib, N., et al. 2019, aj, 157, 253

work page 2019
[19]

F., Morbidelli, A., Van Laerhoven, C., Gomes, R., & Tsiganis, K

Levison, H. F., Morbidelli, A., Van Laerhoven, C., Gomes, R., & Tsiganis, K. 2008, Icarus, 196, 258

work page 2008
[20]

& Zhou, L.-Y

Li, H. & Zhou, L.-Y . 2023, A&A, 680, A68 ADS

work page 2023
[21]

& Zhou, L.-Y

Li, H. & Zhou, L.-Y . 2024, A&A, 687, A206 ADS

work page 2024
[22]

2014, MNRAS, 443, 1346

Li, J., Zhou, L.-Y ., & Sun, Y .-S. 2014, MNRAS, 443, 1346

work page 2014
[23]

2023, icarus, 404, 115650

Liu, S., Wu, Z., Yan, J., et al. 2023, icarus, 404, 115650

work page 2023
[24]

Lykawka, P. S. & Mukai, T. 2007, Icarus, 192, 238

work page 2007
[25]

1995, AJ, 110, 420 ADS

Malhotra, R. 1995, AJ, 110, 420 ADS

work page 1995
[26]

1996, aj, 111, 504

Malhotra, R. 1996, aj, 111, 504

work page 1996
[27]

& Wang, X

Malhotra, R. & Wang, X. 2017, mnras, 465, 4381

work page 2017
[28]

Melita, M. D. & Brunini, A. 2000, Icarus, 147, 205 ADS

work page 2000
[29]

Message, P. J. 1958, AJ, 63, 443

work page 1958
[30]

2002, Modern celestial mechanics : aspects of solar system dy- namics

Morbidelli, A. 2002, Modern celestial mechanics : aspects of solar system dy- namics

work page 2002
[31]

& Nesvorný, D

Morbidelli, A. & Nesvorný, D. 2020, in The Trans-Neptunian Solar System, ed. D. Prialnik, M. A. Barucci, & L. Young, 25–59

work page 2020
[32]

1995, icarus, 118, 322

Morbidelli, A., Thomas, F., & Moons, M. 1995, icarus, 118, 322

work page 1995
[33]

Murray, C. D. & Dermott, S. F. 1999, Solar System Dynamics

work page 1999
[34]

Murray-Clay, R. A. & Chiang, E. I. 2005, apj, 619, 623

work page 2005
[35]

Mustill, A. J. & Wyatt, M. C. 2011, MNRAS, 413, 554 ADS

work page 2011
[36]

& Morais, M

Namouni, F. & Morais, M. H. M. 2015, MNRAS, 446, 1998

work page 2015
[37]

& Morais, M

Namouni, F. & Morais, M. H. M. 2017, MNRAS, 467, 2673 Nesvorný, D. 2015, AJ, 150, 68 Nesvorný, D. 2018, ARA&A, 56, 137 ADS Nesvorný, D. & Morbidelli, A. 2012, AJ, 144, 117 Nesvorný, D. & V okrouhlický, D. 2016, ApJ, 825, 94

work page 2017
[38]

Peale, S. J. 1976, ARA&A, 14, 215

work page 1976
[39]

E., Kavelaars, J

Pike, R. E., Kavelaars, J. J., Petit, J. M., et al. 2015, AJ, 149, 202

work page 2015
[40]

Plummer, H. C. 1916, mnras, 76, 390

work page 1916
[41]

Quillen, A. C. 2006, MNRAS, 365, 1367

work page 2006
[42]

& Liu, S

Rein, H. & Liu, S. F. 2012, A&A, 537, A128 ADS

work page 2012
[43]

& Tamayo, D

Rein, H. & Tamayo, D. 2015, MNRAS, 452, 376

work page 2015
[44]

Saillenfest, M., Fouchard, M., Tommei, G., & Valsecchi, G. B. 2016, Celestial Mechanics and Dynamical Astronomy, 126, 369

work page 2016
[45]

Tiscareno, M. S. & Malhotra, R. 2009, aj, 138, 827

work page 2009
[46]

Tsiganis, K., Gomes, R., Morbidelli, A., & Levison, H. F. 2005, Nature, 435, 459 V olk, K., Murray-Clay, R., Gladman, B., et al. 2016, AJ, 152, 23 ADS V olk, K., Murray-Clay, R. A., Gladman, B. J., et al. 2018, AJ, 155, 260 V oyatzis, G., Kotoulas, T., & Hadjidemetriou, J. D. 2005, Celestial Mechanics and Dynamical Astronomy, 91, 191

work page 2005
[47]

& Holman, M

Wisdom, J. & Holman, M. 1991, AJ, 102, 1528

work page 1991
[48]

Wyatt, M. C. 2003, ApJ, 598, 1321 ADS

work page 2003
[49]

Yoder, C. F. 1979, Celestial Mechanics, 19, 3

work page 1979
[50]

Yu, T. Y . M., Murray-Clay, R., & V olk, K. 2018, AJ, 156, 33 ADS Article number, page 12 of 12

work page 2018