The Neptunian ridge as a natural outcome of high-eccentricity tidal migration
Pith reviewed 2026-05-21 09:06 UTC · model grok-4.3
The pith
High-eccentricity tidal migration produces the Neptunian ridge as a density-dependent survival band after circularization.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The HEM tidal survival formalism reproduces the slope of the desert boundary across the 1.8 to 6 Earth-radius regime with a single tidal encounter parameter. In the Jovian regime the boundary stays broadly consistent, with small deviations attributable to radius inflation or decay. Incorporating observed density dispersion converts the disruption limit into a tidal survival band that traces the ridge, and because dissipation increases sharply toward the threshold, HEM survivors circularize just beyond it and cluster within the band, producing the observed overdensity. In the period-density plane the population follows the predicted density-dependent survival pattern with a persistent ridge-1
What carries the argument
The high-eccentricity tidal migration survival constraints mapped onto the period-radius plane via empirical mass-radius relations and broadened by observed density dispersion into a finite survival band.
If this is right
- The desert boundary slope is reproduced from sub-Neptunes through super-Neptunes with one fixed tidal encounter parameter.
- In the Jovian regime the boundary remains consistent with the survival limit aside from possible radius inflation effects.
- Ridge planets are expected to cluster near a density of 1.7 g cm^{-3} in the period-density plane.
- Steeply rising tidal dissipation near the limit forces survivors to circularize inside the survival band.
Where Pith is reading between the lines
- Similar ridges could appear in other radius or period ranges if comparable density dispersions exist in those populations.
- Tighter constraints on individual planet densities in the 3-6 day range would directly test whether the ridge sits inside the predicted survival band.
- The model implies that revisions to the adopted mass-radius relation would shift the expected ridge location in a predictable way.
Load-bearing premise
The observed dispersion in planet densities is large enough to turn the sharp tidal disruption limit into a survival band whose position and width after circularization match the ridge overdensity.
What would settle it
A high-precision census showing that ridge planets do not concentrate near 1.7 g cm^{-3} or that the ridge overdensity lies outside the density-dependent tidal survival band predicted by HEM.
Figures
read the original abstract
Recent occurrence-rate analyses have shown that the transition between the Neptunian desert and the savanna is not smooth but instead exhibits an overdensity of planets at $P_{\rm orb}\simeq3$-$6$ d, known as the Neptunian ridge. We confronted the high-eccentricity tidal migration (HEM) scenario with this updated desert-ridge-savanna landscape. We mapped the HEM tidal survival constraints onto the period-radius plane using empirically inferred mass-radius relations and provided an independent consistency check in the period-density plane. The HEM tidal survival formalism reproduces the slope of the desert boundary across the sub-Neptune to super-Neptune/sub-Saturn regime ($1.8\,\rm R_\oplus \lesssim R_{\rm p} \lesssim 6\,\rm R_\oplus$), with a single representative tidal encounter parameter setting the overall period offset. In the Jovian regime, the boundary remains broadly consistent with the survival limit, with residual deviations likely due to radius inflation or orbital decay. Incorporating the observed density dispersion transforms the disruption limit into a finite tidal survival band that traces the ridge. Because tidal dissipation rises steeply towards the disruption threshold, HEM survivors are expected to circularise just beyond this limit, clustering within the band and naturally producing the ridge overdensity. In the period-density plane, the population follows the predicted density-dependent survival and clustering pattern, with a persistent concentration of ridge planets near $\rho_{\rm p}\simeq1.7\,\mathrm{g\,cm^{-3}}$. High-eccentricity tidal migration thus provides a self-consistent explanation for the ridge and desert boundary geometry.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that high-eccentricity tidal migration (HEM) provides a self-consistent explanation for both the slope of the Neptunian desert boundary and the ridge overdensity at P_orb ≈ 3-6 d. Using empirically inferred mass-radius relations, the authors map HEM tidal survival constraints onto the period-radius plane, reproducing the observed desert slope across 1.8 R_⊕ ≲ R_p ≲ 6 R_⊕ with a single representative tidal encounter parameter that sets the overall period offset. They further argue that the observed dispersion in planet densities converts the sharp tidal disruption limit into a finite survival band whose location and width trace the ridge; because tidal dissipation increases steeply near the threshold, survivors circularize and cluster inside this band. An independent consistency check is presented in the period-density plane, where the population follows the predicted density-dependent pattern with ridge planets concentrated near ρ_p ≈ 1.7 g cm^{-3}.
Significance. If the quantitative details of the density-to-band mapping can be shown to align with the observed ridge amplitude and width without further tuning, the work would supply a dynamical mechanism that simultaneously accounts for the desert boundary geometry and the ridge overdensity, strengthening the case for HEM relative to in-situ or disk-migration scenarios. The independent period-density consistency check is a methodological strength that allows falsifiability. However, the reliance on an empirically chosen tidal parameter to fix the period offset limits the predictive power of the model.
major comments (3)
- [period-radius mapping and abstract] The reproduction of the desert boundary slope (abstract and the period-radius mapping section) is achieved by selecting a single representative tidal encounter parameter whose value is set to match the observed boundary location. This choice makes the overall period offset a tuned quantity rather than an a-priori prediction, directly affecting the central claim that HEM is self-consistent with the observed landscape.
- [density dispersion to survival band] The transformation of the sharp tidal disruption limit into a finite survival band that traces the ridge (abstract and the density-dispersion discussion) invokes the observed density dispersion but provides no explicit mapping: which quantiles of the density distribution define the band edges, how circularization shifts the final periods, or whether the resulting overdensity amplitude quantitatively matches the observed ridge strength. This step is load-bearing for the ridge explanation yet remains unverified in the presented results.
- [Jovian regime discussion] In the Jovian regime the boundary is described as broadly consistent but with residual deviations attributed to radius inflation or orbital decay. These exceptions indicate that the same HEM formalism is not uniformly predictive across the full radius range, weakening the assertion of a self-consistent explanation for the entire desert-ridge-savanna structure.
minor comments (2)
- [abstract] The abstract states that the population follows the predicted pattern 'with a persistent concentration of ridge planets near ρ_p ≃ 1.7 g cm^{-3}' but does not reference the specific figure or table that demonstrates this concentration; adding an explicit cross-reference would improve clarity.
- [methods / formalism] Notation for the tidal encounter parameter is introduced without an equation label or explicit functional form in the summary sections; defining it once with a numbered equation would aid readers tracing the single-parameter choice.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed report. We address each major comment below and indicate where revisions will be made to the manuscript.
read point-by-point responses
-
Referee: [period-radius mapping and abstract] The reproduction of the desert boundary slope (abstract and the period-radius mapping section) is achieved by selecting a single representative tidal encounter parameter whose value is set to match the observed boundary location. This choice makes the overall period offset a tuned quantity rather than an a-priori prediction, directly affecting the central claim that HEM is self-consistent with the observed landscape.
Authors: The slope of the boundary in the period-radius plane is the central prediction of the HEM survival criterion once the empirical mass-radius relation is adopted; this slope emerges without any adjustment of the tidal parameter. The representative tidal encounter parameter is introduced only to fix the absolute period offset for direct comparison with observations and corresponds to a characteristic value expected from high-eccentricity dynamical histories. We have revised the abstract and mapping section to clarify this distinction between the untuned slope prediction and the normalized offset. revision: partial
-
Referee: [density dispersion to survival band] The transformation of the sharp tidal disruption limit into a finite survival band that traces the ridge (abstract and the density-dispersion discussion) invokes the observed density dispersion but provides no explicit mapping: which quantiles of the density distribution define the band edges, how circularization shifts the final periods, or whether the resulting overdensity amplitude quantitatively matches the observed ridge strength. This step is load-bearing for the ridge explanation yet remains unverified in the presented results.
Authors: We agree that an explicit quantitative mapping strengthens the ridge explanation. In the revised manuscript we now specify that the band edges are set by the 16th and 84th percentiles of the observed density distribution at each radius. Circularization is modeled by shifting planets that reach the disruption threshold to the circular period corresponding to their post-damping semi-major axis. We also include a direct comparison showing that the predicted overdensity amplitude lies within the uncertainties of the measured ridge strength. These details have been added to the density-dispersion section. revision: yes
-
Referee: [Jovian regime discussion] In the Jovian regime the boundary is described as broadly consistent but with residual deviations attributed to radius inflation or orbital decay. These exceptions indicate that the same HEM formalism is not uniformly predictive across the full radius range, weakening the assertion of a self-consistent explanation for the entire desert-ridge-savanna structure.
Authors: The manuscript's primary claim of self-consistency applies to the Neptunian regime (1.8–6 R_⊕). In the Jovian regime we already state that the boundary is only broadly consistent and explicitly attribute the residuals to radius inflation and orbital decay—processes that become relevant at larger radii and are outside the basic HEM survival calculation. We have revised the text to make the limited scope of the uniform prediction clearer while retaining the transparent discussion of the exceptions. revision: partial
Circularity Check
Single tidal encounter parameter tuned to observed boundary offset; density dispersion then defines survival band that traces ridge by construction
specific steps
-
fitted input called prediction
[Abstract]
"The HEM tidal survival formalism reproduces the slope of the desert boundary across the sub-Neptune to super-Neptune/sub-Saturn regime (1.8 R⊕ ≲ Rp ≲ 6 R⊕), with a single representative tidal encounter parameter setting the overall period offset."
The representative parameter is chosen to align the overall period offset with the observed desert boundary; the subsequent claim that the formalism reproduces the boundary therefore incorporates a fitted quantity as if it were an independent prediction.
-
self definitional
[Abstract]
"Incorporating the observed density dispersion transforms the disruption limit into a finite tidal survival band that traces the ridge. Because tidal dissipation rises steeply towards the disruption threshold, HEM survivors are expected to circularise just beyond this limit, clustering within the band and naturally producing the ridge overdensity."
The finite survival band is constructed by mapping the observed density dispersion onto the disruption limit; the resulting band is then asserted to trace the ridge and to produce the observed overdensity, rendering the explanation equivalent to the input data by construction.
full rationale
The derivation uses empirically inferred mass-radius relations to map HEM survival limits, but selects one representative tidal encounter parameter to set the overall period offset so that the formalism reproduces the desert boundary location. Incorporating the observed density dispersion then broadens the sharp disruption limit into a finite band whose location and width are stated to trace the ridge overdensity. Because the parameter choice and the dispersion mapping are both anchored directly to the same observed features the model seeks to explain, the central claim that HEM naturally produces the ridge reduces to a fitted input relabeled as a prediction. The slope match and period-density consistency check retain some independent content, preventing a higher score.
Axiom & Free-Parameter Ledger
free parameters (1)
- tidal encounter parameter =
representative value chosen to match observed boundary
axioms (1)
- domain assumption Empirically inferred mass-radius relations accurately map tidal survival constraints from mass to radius for planets between 1.8 and 6 Earth radii.
Reference graph
Works this paper leans on
- [1]
-
[2]
Attia , M., Bourrier , V., Eggenberger , P., et al. 2021, , 647, A40
work page 2021
- [3]
- [4]
-
[5]
Bourrier , V., Steiner , M., Castro-Gonz \'a lez , A., et al. 2025, , 701, A190
work page 2025
-
[6]
Castro-Gonz \'a lez , A., Bouchy , F., Correia , A. C. M., et al. 2025, , 699, A344
work page 2025
-
[7]
Castro-Gonz \'a lez , A., Bourrier , V., Lillo-Box , J., et al. 2024 a , , 689, A250
work page 2024
- [8]
-
[9]
1963, The Astrophysical Journal, 138, 1182
Chandrasekhar, S. 1963, The Astrophysical Journal, 138, 1182
work page 1963
-
[10]
Christiansen , J. L., McElroy , D. L., Harbut , M., et al. 2025, , 6, 186
work page 2025
-
[11]
Correia , A. C. M., Bourrier , V., & Delisle , J.-B. 2020, , 635, A37
work page 2020
-
[12]
Correia , A. C. M., Laskar , J., Farago , F., & Bou \'e , G. 2011, Celestial Mechanics and Dynamical Astronomy, 111, 105
work page 2011
-
[13]
J., Hadjigeorghiou , A., et al
Cui , K., Armstrong , D. J., Hadjigeorghiou , A., et al. 2026, , 546, stag022
work page 2026
-
[14]
Dawson , R. I. & Johnson , J. A. 2018, , 56, 175
work page 2018
-
[15]
Doyle , L., Armstrong , D. J., Acu \ n a , L., et al. 2025, , 539, 3138
work page 2025
- [16]
-
[17]
Ford , E. B. & Rasio , F. A. 2006, , 638, L45
work page 2006
-
[18]
Ford , E. B. & Rasio , F. A. 2008, , 686, 621
work page 2008
- [19]
- [20]
- [21]
- [22]
-
[23]
Handley , L. B. & Batygin , K. 2026, , 998, 324
work page 2026
-
[24]
Handley , L. B., Howard , A. W., Rubenzahl , R. A., et al. 2025, , 169, 212
work page 2025
- [25]
-
[26]
Liu , S.-F., Guillochon , J., Lin , D. N. C., & Ramirez-Ruiz , E. 2013, , 762, 37
work page 2013
- [27]
- [28]
- [29]
-
[30]
Owen , J. E. & Lai , D. 2018, , 479, 5012
work page 2018
-
[31]
1849, original derivation of the tidal limit (Roche limit)
Roche, E. 1849, original derivation of the tidal limit (Roche limit)
-
[32]
Thompson , S. E., Coughlin , J. L., Hoffman , K., et al. 2018, , 235, 38
work page 2018
- [33]
- [34]
-
[35]
W., Tamburo , P., Stef \'a nsson , G., et al
Yee , S. W., Tamburo , P., Stef \'a nsson , G., et al. 2025, , 170, 275
work page 2025
-
[36]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[37]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.