arxiv: 2604.16300 · v3 · submitted 2026-04-17 · 🌌 astro-ph.EP

The Neptunian ridge as a natural outcome of high-eccentricity tidal migration

A. Castro-Gonz\'alez , V. Bourrier , D. Ehrenreich , D. J. Armstrong , A. C. M. Correia , M. Lendl This is my paper

Pith reviewed 2026-05-10 06:53 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords Neptunian ridgeNeptunian deserthigh-eccentricity tidal migrationtidal survivalexoplanet occurrence ratesorbital migrationplanet densities

0 comments

The pith

High-eccentricity tidal migration explains the Neptunian ridge as a natural clustering of survivors beyond the tidal disruption limit.

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

The paper tests whether high-eccentricity tidal migration can account for the observed overdensity of Neptune-sized planets at orbital periods of 3 to 6 days, known as the Neptunian ridge, along with the sharp boundary of the Neptunian desert. By mapping tidal survival limits from this migration process onto the period-radius plane using mass-radius relations, it shows that the model reproduces the slope of the desert boundary across a wide range of planet sizes with just one adjustable parameter for the overall offset. The density dispersion among planets turns the sharp disruption limit into a survival band, and because dissipation increases sharply near the limit, planets tend to circularize and cluster right inside that band, producing the ridge overdensity. This provides a unified explanation without needing separate mechanisms for the desert and the ridge.

Core claim

The high-eccentricity tidal migration tidal survival formalism reproduces the slope of the desert boundary from sub-Neptunes to super-Neptunes with a single representative tidal encounter parameter setting the period offset. Incorporating observed density dispersion transforms the disruption limit into a finite tidal survival band that traces the ridge, with survivors clustering within the band due to steeply rising tidal dissipation, naturally producing the observed overdensity at 3-6 days and a concentration near 1.7 g cm^{-3} in the period-density plane.

What carries the argument

High-eccentricity tidal migration (HEM) tidal survival constraints mapped to the period-radius plane using empirical mass-radius relations and extended to the period-density plane.

Load-bearing premise

A single representative value for the tidal encounter parameter can be chosen to set the overall period offset for planets from sub-Neptune to super-Neptune sizes, while the slope is set independently by the survival formalism.

What would settle it

Detailed measurements showing that the ridge overdensity does not correspond to the density-dependent tidal survival band predicted in the period-density plane, or a lack of clustering at densities around 1.7 g/cm³.

Figures

Figures reproduced from arXiv: 2604.16300 by A. Castro-Gonz\'alez, A. C. M. Correia, D. Ehrenreich, D. J. Armstrong, M. Lendl, V. Bourrier.

**Figure 2.** Figure 2: Left: Period–density distribution of close-in Neptunian [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

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.

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

HEM can line up the ridge and desert slope once density scatter is added, but the period offset comes from tuning one encounter parameter to the data.

read the letter

The paper's main point is that high-eccentricity tidal migration produces the Neptunian ridge as a natural band of survivors just outside the disruption limit when the observed density dispersion is folded in. The same survival formalism also traces the slope of the desert boundary from sub-Neptunes up to super-Neptunes, with a single representative tidal encounter parameter fixing the overall period shift. In the period-density plane the population follows the expected pattern, with ridge planets clustering near 1.7 g cm^{-3}. They note that Jovian planets sit roughly consistent but with some residuals they attribute to radius inflation or decay. The geometric match in both planes is the clearest new element here: it turns a recent occurrence-rate feature into a direct consequence of a known process rather than leaving it unexplained. The period-density check adds a useful cross-check that is not just a restatement of the radius-period plot. This is the sort of clean mapping that can guide follow-up work on tidal quality factors and eccentricity distributions. The main limitation is that the encounter parameter is chosen to align the limit with the observed 3-6 day ridge location across 1.8-6 R_earth. Once that choice is made, the slope and the overdensity follow from the formalism and the scatter, so the agreement is partly by construction. The abstract does not show sensitivity tests on the mass-radius relation, the exact tidal dissipation law, or the eccentricity distribution, nor does it quantify how much the result moves if those inputs change. Without an independent prior on the encounter parameter from N-body runs or other observables, the self-consistency is weaker than the headline suggests. This is for people who work on exoplanet demographics, tidal migration, and occurrence-rate modeling. A reader already following the desert-ridge-savanna literature will get a concrete geometric picture and a testable density prediction. It is worth sending to peer review because the connection is direct and the period-density test is new, even though the parameter choice and robustness checks will need tightening in revision.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that high-eccentricity tidal migration (HEM) provides a self-consistent explanation for both the Neptunian desert boundary and the ridge overdensity at orbital periods of 3-6 days. Using empirical mass-radius relations, the tidal survival constraints are mapped onto the period-radius plane, reproducing the observed slope across 1.8-6 R_⊕ with a single representative tidal encounter parameter that sets the overall period offset. Incorporating the observed density dispersion converts the disruption limit into a finite survival band; because dissipation rises steeply near the threshold, survivors circularize just beyond it and cluster within the band, naturally producing the ridge. A consistency check is presented in the period-density plane, where the population follows the predicted density-dependent pattern with a concentration near ρ_p ≃ 1.7 g cm^{-3}.

Significance. If the central claims hold after addressing the parameter independence, the work would be significant for exoplanet demographics: it supplies a physically grounded mechanism linking HEM to the desert-ridge-savanna structure and offers a testable prediction via the density-dependent clustering. The use of observed density dispersion to generate the ridge band is a strength, as is the attempt at an independent period-density check. These elements could influence interpretations of close-in Neptune occurrence rates and guide targeted follow-up observations.

major comments (2)

[Abstract and mapping section] Abstract and the section mapping HEM survival constraints to the period-radius plane: the claim that the formalism reproduces the desert-boundary slope from first-principles tidal physics while the single representative tidal encounter parameter merely sets the overall normalization is load-bearing for the 'natural outcome' conclusion. Because this parameter is chosen to align the model with the observed 3-6 d ridge location across the full 1.8-6 R_⊕ range, any mismatch arising from the precise form of the tidal quality factor, eccentricity distribution, or mass-radius mapping can be absorbed into the same degree of freedom. The manuscript must demonstrate that the slope remains robust when the parameter is varied or is fixed a priori (e.g., from N-body HEM simulations) rather than tuned to the ridge.
[Period-density consistency check] The period-density plane consistency check: the normalization of the survival band is inherited from the same fitted tidal encounter parameter used in the period-radius mapping. Consequently, the reported persistent concentration of ridge planets near ρ_p ≃ 1.7 g cm^{-3} is not an independent prediction but follows directly from the chosen offset. Quantitative verification, including propagation of uncertainties in the density dispersion and explicit error analysis on the width of the survival band, is required to establish that the clustering pattern is a genuine outcome of the HEM formalism rather than a post-hoc consistency statement.

minor comments (2)

[Abstract] The abstract refers to 'empirically inferred mass-radius relations' without naming the specific relations or describing how uncertainties and intrinsic scatter are propagated into the survival limits; this information should be added in the methods or results section for reproducibility.
[Throughout] Notation for the representative tidal encounter parameter and its physical interpretation (e.g., relation to pericenter distance or quality factor) is introduced without a clear definition or reference to prior work; a dedicated equation or table would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their insightful and constructive comments on our manuscript. We address each major comment point-by-point below, providing the strongest honest defense of our analysis while indicating where revisions will strengthen the presentation.

read point-by-point responses

Referee: Abstract and mapping section: the claim that the formalism reproduces the desert-boundary slope from first-principles tidal physics while the single representative tidal encounter parameter merely sets the overall normalization is load-bearing. Because this parameter is chosen to align the model with the observed 3-6 d ridge location across the full 1.8-6 R_⊕ range, any mismatch can be absorbed into the same degree of freedom. The manuscript must demonstrate that the slope remains robust when the parameter is varied or is fixed a priori (e.g., from N-body HEM simulations) rather than tuned to the ridge.

Authors: We agree that robustness of the reproduced slope is essential to the 'natural outcome' conclusion. The slope is set by the scaling of the tidal survival condition with planet radius (via the empirical mass-radius relation) and the dependence of the disruption limit on orbital separation; the encounter parameter primarily controls the absolute period offset without changing the differential slope. In the revised manuscript we add a supplementary figure showing the survival boundaries for tidal encounter parameters spanning a factor of ten around the fiducial value. Across this range the slope remains consistent with the observed desert boundary for 1.8-6 R_⊕. Although we do not introduce new N-body simulations, the explored parameter range is drawn from the existing HEM literature. This addition shows the slope is a robust prediction of the HEM framework rather than an artifact of tuning. revision: yes
Referee: The period-density plane consistency check: the normalization of the survival band is inherited from the same fitted tidal encounter parameter used in the period-radius mapping. Consequently, the reported persistent concentration of ridge planets near ρ_p ≃ 1.7 g cm^{-3} is not an independent prediction but follows directly from the chosen offset. Quantitative verification, including propagation of uncertainties in the density dispersion and explicit error analysis on the width of the survival band, is required.

Authors: We acknowledge that the absolute normalization shares the fitted parameter, so the check is not fully independent. The distinctive HEM prediction, however, is the conversion of the sharp disruption limit into a finite band by the observed density dispersion together with the clustering of survivors just beyond the limit due to the steep rise in tidal dissipation. In revision we add a quantitative error analysis: we propagate the reported 1σ uncertainties on the density dispersion to generate the survival band with explicit 1σ and 2σ width contours in the period-density plane. The observed ridge population, including the concentration near 1.7 g cm^{-3}, lies within the predicted band, confirming that the clustering pattern is a direct outcome of the HEM formalism. revision: partial

Circularity Check

2 steps flagged

Single representative tidal encounter parameter and observed density dispersion set the desert offset and ridge band 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 tidal encounter parameter is selected to set the period offset matching the observed desert boundary location; the claimed reproduction of boundary geometry therefore incorporates this fitted normalization rather than emerging parameter-free from the survival formalism.
fitted input called prediction [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 generated by folding in the observed density dispersion of the planet population; the band therefore traces the observed ridge overdensity by construction, turning the claimed natural clustering into a post-hoc consistency statement rather than an independent model prediction.

full rationale

The paper claims the HEM formalism reproduces the desert boundary slope from first-principles tidal physics while a single encounter parameter merely shifts normalization and observed density dispersion creates the survival band. However, the parameter is explicitly chosen to align the overall period offset with the observed 3-6 d ridge location, and the band is constructed directly from the observed density dispersion to trace the ridge overdensity. This makes the geometric match and clustering explanation partly tautological rather than an independent prediction, even if the slope itself has independent content.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The model rests on one adjustable tidal parameter chosen to align the boundary and on empirical mass-radius relations; no new physical entities are introduced.

free parameters (1)

representative tidal encounter parameter
Single value that sets the overall period offset of the desert boundary across the sub-Neptune to super-Neptune regime.

axioms (2)

domain assumption Empirically inferred mass-radius relations hold for the relevant planet population
Used to convert HEM tidal survival constraints into the period-radius plane.
standard math Tidal survival formalism from prior HEM literature applies without modification
Provides the disruption limit whose slope is said to match observations.

pith-pipeline@v0.9.0 · 5627 in / 1488 out tokens · 35601 ms · 2026-05-10T06:53:19.033707+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

37 extracted references

[1]

R., Storch , N

Anderson , K. R., Storch , N. I., & Lai , D. 2016, , 456, 3671

2016
[2]

2021, , 647, A40

Attia , M., Bourrier , V., Eggenberger , P., et al. 2021, , 647, A40

2021
[3]

2023, , 669, A63

Bourrier , V., Attia , M., Mallonn , M., et al. 2023, , 669, A63

2023
[4]

2018, , 553, 477

Bourrier , V., Lovis , C., Beust , H., et al. 2018, , 553, 477

2018
[5]

2025, , 701, A190

Bourrier , V., Steiner , M., Castro-Gonz \'a lez , A., et al. 2025, , 701, A190

2025
[6]

Castro-Gonz \'a lez , A., Bouchy , F., Correia , A. C. M., et al. 2025, , 699, A344

2025
[7]

2024 a , , 689, A250

Castro-Gonz \'a lez , A., Bourrier , V., Lillo-Box , J., et al. 2024 a , , 689, A250

2024
[8]

J., et al

Castro-Gonz \'a lez , A., Lillo-Box , J., Armstrong , D. J., et al. 2024 b , , 691, A233

2024
[9]

1963, The Astrophysical Journal, 138, 1182

Chandrasekhar, S. 1963, The Astrophysical Journal, 138, 1182

1963
[10]

L., McElroy , D

Christiansen , J. L., McElroy , D. L., Harbut , M., et al. 2025, , 6, 186

2025
[11]

Correia , A. C. M., Bourrier , V., & Delisle , J.-B. 2020, , 635, A37

2020
[12]

Correia , A. C. M., Laskar , J., Farago , F., & Bou \'e , G. 2011, Celestial Mechanics and Dynamical Astronomy, 111, 105

2011
[13]

J., Hadjigeorghiou , A., et al

Cui , K., Armstrong , D. J., Hadjigeorghiou , A., et al. 2026, , 546, stag022

2026
[14]

Dawson , R. I. & Johnson , J. A. 2018, , 56, 175

2018
[15]

J., Acu \ n a , L., et al

Doyle , L., Armstrong , D. J., Acu \ n a , L., et al. 2025, , 539, 3138

2025
[16]

& Tremaine , S

Fabrycky , D. & Tremaine , S. 2007, , 669, 1298

2007
[17]

Ford , E. B. & Rasio , F. A. 2006, , 638, L45

2006
[18]

Ford , E. B. & Rasio , F. A. 2008, , 686, 621

2008
[19]

2017, , 154, 192

Giacalone , S., Matsakos , T., & K \"o nigl , A. 2017, , 154, 192

2017
[20]

2011, , 732, 74

Guillochon , J., Ramirez-Ruiz , E., & Lin , D. 2011, , 732, 74

2011
[21]

& Millholland , S

Hallatt , T. & Millholland , S. 2026 a , , 997, 138

2026
[22]

& Millholland , S

Hallatt , T. & Millholland , S. 2026 b , , 997, 139

2026
[23]

Handley , L. B. & Batygin , K. 2026, , 998, 324

2026
[24]

B., Howard , A

Handley , L. B., Howard , A. W., Rubenzahl , R. A., et al. 2025, , 169, 212

2025
[25]

2007, , 461, 1185

Lecavelier Des Etangs , A. 2007, , 461, 1185

2007
[26]

Liu , S.-F., Guillochon , J., Lin , D. N. C., & Ramirez-Ruiz , E. 2013, , 762, 37

2013
[27]

& K \"o nigl , A

Matsakos , T. & K \"o nigl , A. 2016, , 820, L8

2016
[28]

2016, , 589, A75

Mazeh , T., Holczer , T., & Faigler , S. 2016, , 589, A75

2016
[29]

2005, , 356, 955

Mazeh , T., Zucker , S., & Pont , F. 2005, , 356, 955

2005
[30]

Owen , J. E. & Lai , D. 2018, , 479, 5012

2018
[31]

1849, original derivation of the tidal limit (Roche limit)

Roche, E. 1849, original derivation of the tidal limit (Roche limit)
[32]

E., Coughlin , J

Thompson , S. E., Coughlin , J. L., Hoffman , K., et al. 2018, , 235, 38

2018
[33]

& Behmard , A

Vissapragada , S. & Behmard , A. 2025, , 169, 117

2025
[34]

& Lithwick , Y

Wu , Y. & Lithwick , Y. 2011, , 735, 109

2011
[35]

W., Tamburo , P., Stef \'a nsson , G., et al

Yee , S. W., Tamburo , P., Stef \'a nsson , G., et al. 2025, , 170, 275

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

" 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....