pith. machine review for the scientific record. sign in

arxiv: 2604.25704 · v1 · submitted 2026-04-28 · 🌌 astro-ph.EP

Recognition: unknown

Suppression of Resonant Overstability at Sharp Migration Gradients

Konstantin Batygin , Ian R. Brunton , Alessandro Morbidelli
Authors on Pith no claims yet

Pith reviewed 2026-05-07 14:34 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords mean-motion resonancesresonant overstabilityprotoplanetary disksorbital migrationtorque reversalsplanet-disk interactionsN-body simulations
0
0 comments X

The pith

Resonant overstability is quenched when migration gradients are sufficiently sharp.

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

Mean-motion resonances are expected to form at the inner edges of protoplanetary disks through large-scale orbital convergence. Dissipative forces that enable capture can also drive librations overstable and break the resonances. This paper develops an analytical theory reducing the dynamics to a damped driven oscillator and shows that overstability is suppressed once the local gradient steepness parameter exceeds the ratio of semi-major axis to eccentricity evolution timescales. The result is verified with N-body integrations. A reader would care because this determines whether resonant chains persist or dissolve, shaping the final architectures of planetary systems.

Core claim

Adopting the dissipative circular restricted three-body problem as a paradigm, a WKB-style analysis reduces the resonant dynamics to a damped, driven harmonic oscillator. Within this framework an effective frictional term arises that is proportional to the local migration-rate gradient, encoded in the dimensionless coefficient beta. The analytical theory predicts that overstability is quenched once beta is at least as large as tau_a over tau_e. Direct N-body integrations verify and refine the analytic results, while type-I scalings indicate that the competition is a borderline outcome sensitive to the detailed structure of planet-disk interactions.

What carries the argument

WKB-style reduction of the resonant dynamics to a damped driven harmonic oscillator whose effective friction is set by the dimensionless steepness parameter beta of the local migration gradient.

If this is right

  • Overstability is subdued once the torque reversal becomes sufficiently sharp.
  • The stability threshold is set by the ratio of the semi-major axis and eccentricity evolution timescales.
  • N-body integrations confirm the analytic quenching condition.
  • The outcome remains sensitive to the precise radial structure of disk torques.

Where Pith is reading between the lines

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

  • Resonance survival probabilities could shift under disk models with different torque profiles.
  • The slow-variation assumption could be relaxed to treat abrupt transitions directly.
  • Analogous gradient effects may appear in other dissipative resonance problems in celestial mechanics.

Load-bearing premise

The migration gradient varies slowly compared to the resonant libration period.

What would settle it

N-body runs that scan the gradient steepness across the predicted threshold and measure whether libration amplitudes grow or decay exactly when beta falls below tau_a over tau_e.

Figures

Figures reproduced from arXiv: 2604.25704 by Alessandro Morbidelli, Ian R. Brunton, Konstantin Batygin.

Figure 1
Figure 1. Figure 1: Schematic of the dissipative circular restricted three-body problem. A massless test particle resides in an interior k : k − 1 resonance (here, taken to be k = 3) with a circular exterior perturber, and evolves subject to disk￾driven evolution of a and e. The resonance width (cor￾responding to a perturber of mass m′ = 10−4 M) is illus￾trated with bold black curves while the resonant equilibrium is shown wi… view at source ↗
Figure 2
Figure 2. Figure 2: Numerical verification of overstability quench￾ing. Time-series of eccentricity (top) and semi-major axis (bottom) are shown for integrations of the dissipative circu￾lar restricted three-body problem with m′ = 10−4M⋆ and τa/τe = 10, for a range of torque-gradient parameters β. In all runs, outward migration first produces capture into the interior 3:2 resonance and growth of e toward the equi￾librium valu… view at source ↗
read the original abstract

Mean-motion resonances are expected to frequently arise at the inner edges of protoplanetary disks, where planet-disk interactions facilitate large-scale orbital convergence. Under certain conditions, however, the same dissipative forces that promote resonant capture can drive resonant librations overstable, ultimately breaking commensurabilities. Here we examine the onset of overstability near disk torque reversals and show that it can be subdued when the transition is sufficiently sharp. Adopting the dissipative circular restricted three-body problem as a paradigm, we present a WKB-style analysis that reduces the resonant dynamics to a damped, driven harmonic oscillator. Within this framework, we obtain an effective frictional term that is proportional to the local migration-rate gradient, parameterized by a dimensionless coefficient $\beta$ that encodes the steepness of the local torque reversal. Our analytical theory predicts that overstability is quenched once $\beta \gtrsim \tau_a/\tau_e$, where $\tau_a$ and $\tau_e$ denote the characteristic disk-driven evolution timescales of semi-major axis and eccentricity. We verify and refine our analytic results with direct $N$-body integrations. Simple estimates based on conventional type-I scalings suggest that the competition between overstability and its mitigation at disk inner edges is a borderline outcome that is sensitive to the detailed structure of planet-disk interactions.

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 manuscript examines the suppression of resonant overstability in mean-motion resonances near sharp migration gradients at protoplanetary disk edges. Using the dissipative circular restricted three-body problem, it applies a WKB-style reduction of the resonant equations to a damped driven harmonic oscillator, deriving an effective frictional term proportional to the dimensionless gradient parameter β. The central analytic result is that overstability is quenched for β ≳ τ_a/τ_e, where τ_a and τ_e are the disk-driven semi-major axis and eccentricity evolution timescales. This prediction is verified and refined via direct N-body integrations, with estimates suggesting the outcome is borderline under conventional type-I torque scalings.

Significance. If the result holds, the work supplies a first-principles mechanism linking local torque-reversal steepness to resonance stability, with direct relevance to capture and survival of resonances at disk inner edges. The clean analytic reduction to an effective oscillator and the parameter-free relation between β and the timescale ratio are strengths, as is the direct N-body verification that refines the threshold. These elements provide a falsifiable prediction that can be tested against more realistic disk models.

major comments (1)
  1. [analytic derivation] WKB-style reduction (analytic derivation): the mapping to a damped driven oscillator assumes the migration gradient varies slowly compared to the resonant libration period, allowing β to be treated as locally constant over one cycle. However, the quenching regime β ≳ τ_a/τ_e corresponds precisely to sharp transitions (high β) where this slow-variation premise is violated by construction. The manuscript states that N-body runs verify the prediction, but the analytic threshold itself rests on an approximation whose validity range shrinks in the regime where the effect is claimed to operate; a quantitative estimate of the breakdown scale (e.g., ratio of gradient length to libration period) is needed to bound the applicability of the derived friction term.
minor comments (2)
  1. [abstract] The abstract and introduction use τ_a and τ_e without an early explicit definition or reference to their conventional expressions in type-I migration theory; a brief parenthetical or equation reference would improve readability.
  2. [figures] Figure captions for the N-body results should explicitly state the range of β values explored and how the measured libration amplitudes compare to the analytic threshold, to make the verification quantitative.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful and constructive review. The single major comment identifies a legitimate question about the range of validity of the WKB reduction. We address it directly below and have revised the manuscript to incorporate a quantitative bound on the approximation.

read point-by-point responses

  1. Referee: WKB-style reduction (analytic derivation): the mapping to a damped driven oscillator assumes the migration gradient varies slowly compared to the resonant libration period, allowing β to be treated as locally constant over one cycle. However, the quenching regime β ≳ τ_a/τ_e corresponds precisely to sharp transitions (high β) where this slow-variation premise is violated by construction. The manuscript states that N-body runs verify the prediction, but the analytic threshold itself rests on an approximation whose validity range shrinks in the regime where the effect is claimed to operate; a quantitative estimate of the breakdown scale (e.g., ratio of gradient length to libration period) is needed to bound the applicability of the derived friction term.

    Authors: We agree that the slow-variation assumption underlying the WKB reduction must be checked explicitly in the high-β regime. In the revised manuscript we have added a dedicated paragraph (new Section 3.3) that derives the relevant length-scale ratio. Defining the local gradient length as L_grad ≡ a/β and the characteristic libration length as L_lib ≡ a (e / |d e / d t|_lib)^{1/2} (where the libration frequency is taken from the resonant Hamiltonian), we obtain L_grad / L_lib ≈ (τ_a / τ_e) / (β τ_lib / τ_a). At the quenching threshold β ∼ τ_a / τ_e this ratio is O(1) for the fiducial parameters explored in the N-body survey. We therefore state that the leading-order frictional term remains a useful approximation provided β is not orders of magnitude larger than τ_a / τ_e; beyond that point the analytic threshold should be regarded as an order-of-magnitude estimate only. The N-body integrations, which impose no slow-variation assumption, continue to show quenching near the predicted value, indicating that the effective damping captures the dominant physics even when the WKB premise is only marginally satisfied. We have also added a brief caveat in the abstract and conclusions clarifying the regime of applicability. revision: yes

Circularity Check

0 steps flagged

No circularity: first-principles WKB reduction with independent verification

full rationale

The central result (overstability quenched for β ≳ τ_a/τ_e) follows directly from reducing the dissipative circular restricted three-body equations to an effective damped driven oscillator via WKB, where β parameterizes the local torque gradient as an input from disk structure. This is not a fit to resonant outcomes, not a self-definition, and not dependent on self-citation chains. N-body runs are presented only as verification and refinement after the analytic threshold is obtained. No load-bearing step reduces to its own inputs by construction; the derivation remains self-contained under the stated slow-variation assumption.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The analysis rests on the standard dissipative circular restricted three-body problem plus a WKB approximation for slowly varying coefficients; β is introduced as a dimensionless measure of gradient steepness rather than derived from first principles.

free parameters (1)
  • β
    Dimensionless coefficient that encodes the steepness of the local torque reversal; its value is set by the assumed disk structure.
axioms (2)
  • domain assumption Dissipative forces in the circular restricted three-body problem can be modeled with constant τ_a and τ_e timescales.
    Invoked to reduce the resonant dynamics to a linear damped oscillator.
  • domain assumption WKB approximation applies when the migration gradient varies slowly relative to the libration frequency.
    Used to obtain the effective frictional term proportional to the gradient.

pith-pipeline@v0.9.0 · 5533 in / 1447 out tokens · 92516 ms · 2026-05-07T14:34:30.934716+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

40 extracted references · 39 canonical work pages

  1. [1]

    arXiv , author =:2102.08612 , journal =

    Ataiee, S. & Kley, W.\ 2021, , 648, A69. doi:10.1051/0004-6361/202038772

    work page doi:10.1051/0004-6361/202038772 2021

  2. [2]

    & Morbidelli, A.\ 2013, , 556, A28

    Batygin, K. & Morbidelli, A.\ 2013, , 556, A28. doi:10.1051/0004-6361/201220907

    work page doi:10.1051/0004-6361/201220907 2013

  3. [3]

    & Morbidelli, A.\ 2020, , 894, 2, 143

    Batygin, K. & Morbidelli, A.\ 2020, , 894, 2, 143. doi:10.3847/1538-4357/ab8937

    work page doi:10.3847/1538-4357/ab8937 2020

  4. [4]

    2026, , in review

    Batygin, K., & Morbidelli, A. 2026, , in review

    2026

  5. [5]

    D., Drew, J

    Beaug \'e , C., Michtchenko, T. A., & Ferraz-Mello, S.\ 2006, , 365, 4, 1160. doi:10.1111/j.1365-2966.2005.09779.x

    work page doi:10.1111/j.1365-2966.2005.09779.x 2006

  6. [6]

    Brunton, I. R. & Batygin, K.\ 2025, , 991, 1, 15. doi:10.3847/1538-4357/adf432

    work page doi:10.3847/1538-4357/adf432 2025

  7. [7]

    Brunton, I. R. & Batygin, K.\ 2025, , 990, 1, L11. doi:10.3847/2041-8213/adfa94

    work page doi:10.3847/2041-8213/adfa94 2025

  8. [8]

    A.\ 1976, American Journal of Physics, 44, 10, 944

    Burns, J. A.\ 1976, American Journal of Physics, 44, 10, 944. doi:10.1119/1.10237

    work page doi:10.1119/1.10237 1976

  9. [9]

    2023, AJ, 165, 33, doi: 10.3847/1538-3881/aca327

    Dai, F., Masuda, K., Beard, C., et al.\ 2023, , 165, 2, 33. doi:10.3847/1538-3881/aca327

    work page doi:10.3847/1538-3881/aca327 2023

  10. [10]

    Deck, K. M. & Batygin, K.\ 2015, , 810, 2, 119. doi:10.1088/0004-637X/810/2/119

    work page doi:10.1088/0004-637x/810/2/119 2015

  11. [11]

    Delisle, J.-B., Correia, A. C. M., & Laskar, J.\ 2015, , 579, A128. doi:10.1051/0004-6361/201526285

    work page doi:10.1051/0004-6361/201526285 2015

  12. [12]

    F.\ 1968, , 141, 363, doi:10.1093/mnras/141.3.363

    Dermott, S. F.\ 1968, , 141, 363, doi:10.1093/mnras/141.3.363

    work page doi:10.1093/mnras/141.3.363 1968

  13. [13]

    doi:10.48550/arXiv.2510.12895

    Epstein-Martin, M., Stone, N., & Becker, J.\ 2025, arXiv:2510.12895. doi:10.48550/arXiv.2510.12895

    work page doi:10.48550/arxiv.2510.12895 2025

  14. [14]

    Goldreich, P.\ 1965, , 130, 159, doi:10.1093/mnras/130.3.159

    work page doi:10.1093/mnras/130.3.159 1965

  15. [15]

    , volume =

    Goldreich, P. & Tremaine, S.\ 1980, , 241, 425. doi:10.1086/158356

    work page doi:10.1086/158356 1980

  16. [16]

    arXiv , author =:1308.4688 , journal =

    Goldreich, P. & Schlichting, H. E.\ 2014, , 147, 2, 32. doi:10.1088/0004-6256/147/2/32

    work page doi:10.1088/0004-6256/147/2/32 2014

  17. [17]

    doi:10.1093/mnrasl/slv156

    Go \'z dziewski, K., Migaszewski, C., Panichi, F., et al.\ 2016, , 455, 1, L104. doi:10.1093/mnrasl/slv156

    work page doi:10.1093/mnrasl/slv156 2016

  18. [18]

    Greenberg, R.\ 1973, , 165, 305, doi:10.1093/mnras/165.3.305

    work page doi:10.1093/mnras/165.3.305 1973

  19. [19]

    , keywords =

    Grishin, E., Gilbaum, S., & Stone, N. C.\ 2024, , 530, 2, 2114. doi:10.1093/mnras/stae828

    work page doi:10.1093/mnras/stae828 2024

  20. [20]

    doi:10.1007/BF01228946

    Henrard, J.\ 1982, Celestial Mechanics, 27, 1, 3. doi:10.1007/BF01228946

    work page doi:10.1007/bf01228946 1982

  21. [21]

    Lee, M. H. & Peale, S. J.\ 2002, , 567, 1, 596. doi:10.1086/338504

    work page doi:10.1086/338504 2002

  22. [22]

    Lin, D. N. C. & Papaloizou, J.\ 1986, , 309, 846. doi:10.1086/164653

    work page doi:10.1086/164653 1986

  23. [23]

    https://doi.org/10.1038/s41586-023-06692-3

    Luque, R., Osborn, H. P., Leleu, A., et al.\ 2023, , 623, 7989, 932. doi:10.1038/s41586-023-06692-3

    work page doi:10.1038/s41586-023-06692-3 2023

  24. [24]

    , keywords =

    Masset, F. S., Morbidelli, A., Crida, A., et al.\ 2006, , 642, 1, 478. doi:10.1086/500967

    work page doi:10.1086/500967 2006

  25. [25]

    M., Fabrycky, D

    Mills, S. M., Fabrycky, D. C., Migaszewski, C., et al.\ 2016, , 533, 7604, 509. doi:10.1038/nature17445

    work page doi:10.1038/nature17445 2016

  26. [26]

    doi:10.1017/CBO9781139174817

    Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cambridge Univ. Press), doi:10.1017/CBO9781139174817

    work page doi:10.1017/cbo9781139174817 1999

  27. [27]

    doi:10.1016/0021-8928(75)90060-X

    Neishtadt, A.\ 1975, Journal of Applied Mathematics and Mechanics, 39, 4, 594. doi:10.1016/0021-8928(75)90060-X

    work page doi:10.1016/0021-8928(75)90060-x 1975

  28. [28]

    L., & Haiman, Z

    Paardekooper, S.-J. & Papaloizou, J. C. B.\ 2009, , 394, 4, 2283. doi:10.1111/j.1365-2966.2009.14511.x

    work page doi:10.1111/j.1365-2966.2009.14511.x 2009

  29. [29]

    J.\ 1976, , 14, 215

    Peale, S. J.\ 1976, , 14, 215. doi:10.1146/annurev.aa.14.090176.001243

    work page doi:10.1146/annurev.aa.14.090176.001243 1976

  30. [30]

    2023, A&A, 670, A148, doi: 10.1051/0004-6361/202245196

    Pichierri, G., Bitsch, B., & Lega, E.\ 2023, , 670, A148. doi:10.1051/0004-6361/202245196

    work page doi:10.1051/0004-6361/202245196 2023

  31. [31]

    , keywords =

    Rein, H. & Liu, S.-F.\ 2012, , 537, A128. doi:10.1051/0004-6361/201118085

    work page doi:10.1051/0004-6361/201118085 2012

  32. [32]

    Roy, A. E. & Ovenden, M. W.\ 1954, , 114, 232. doi:10.1093/mnras/114.2.232

    work page doi:10.1093/mnras/114.2.232 1954

  33. [33]

    T.\ 1972, , 160, 169, doi:10.1093/mnras/160.2.169

    Sinclair, A. T.\ 1972, , 160, 169, doi:10.1093/mnras/160.2.169

    work page doi:10.1093/mnras/160.2.169 1972

  34. [34]

    , keywords =

    Tamayo, D., Rein, H., Shi, P., et al.\ 2020, , 491, 2, 2885. doi:10.1093/mnras/stz2870

    work page doi:10.1093/mnras/stz2870 2020

  35. [35]

    & Hadden, S.\ 2025, , 986, 1, 11

    Tamayo, D. & Hadden, S.\ 2025, , 986, 1, 11. doi:10.3847/1538-4357/adc1c4

    work page doi:10.3847/1538-4357/adc1c4 2025

  36. [36]

    The Astrophysical Journal , author =

    Tanaka, H., Takeuchi, T., & Ward, W. R.\ 2002, , 565, 2, 1257. doi:10.1086/324713

    work page doi:10.1086/324713 2002

  37. [37]

    R.\ 1997, , 126, 2, 261

    Ward, W. R.\ 1997, , 126, 2, 261. doi:10.1006/icar.1996.5647

    work page doi:10.1006/icar.1996.5647 1997

  38. [38]

    & Lai, D.\ 2017, , 468, 3, 3223

    Xu, W. & Lai, D.\ 2017, , 468, 3, 3223. doi:10.1093/mnras/stx668

    work page doi:10.1093/mnras/stx668 2017

  39. [39]

    Yap, T. E. & Batygin, K.\ 2026, arXiv:2601.00786

    work page arXiv 2026

  40. [40]

    M., & Calvet, N

    Zhu, Z., Stone, J. M., & Calvet, N.\ 2024, , 528, 2, 2883. doi:10.1093/mnras/stad3712

    work page doi:10.1093/mnras/stad3712 2024