pith. sign in

arxiv: 2606.12524 · v1 · pith:4OGFHAE3new · submitted 2026-06-10 · 🌌 astro-ph.EP

An early look at how gas giants shape small planet bulk compositions

Joseph Y. Tang , Marta L. Bryan , Eve J. Lee This is my paper

Pith reviewed 2026-06-27 08:00 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords exoplanetsgas giantsinner small planetsmetal-rich systemsplanet densitycore massoccurrence ratesplanet formation
0
0 comments X

The pith

In metal-rich systems gas giants are preferentially found with lower-density inner small planets of similar core mass.

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

The authors assemble 43 systems containing 68 inner small planets with measured masses and radii, then compute the conditional gas-giant occurrence rate after correcting each system for its own detection sensitivity to distant giants. Across the full sample they find no statistically significant trend with inner-planet density, envelope mass fraction, or core mass. In the metal-rich subset, however, a hint emerges that gas giants accompany lower-density planets while core masses remain comparable. The pattern is consistent with metal-rich disks enabling rapid core growth that triggers early gas accretion, and the lack of a strong core-mass trend may reflect later photoevaporation.

Core claim

After correcting for heterogeneous detection sensitivities, the gas giant occurrence rate P(GG|ISP) shows no significant dependence on inner small planet density, envelope mass fraction, or core mass across the full sample. In metal-rich systems, however, gas giants are preferentially found with lower density planets that have similar core masses. The result is consistent with metal-enriched disks catalyzing rapid core assembly and kickstarting gas accretion early, while muted differences with respect to core mass may indicate contamination by post-formation photoevaporation.

What carries the argument

The conditional gas-giant occurrence rate P(GG|ISP) computed after system-by-system sensitivity corrections, evaluated separately inside the metal-rich subset.

If this is right

  • Metal-enriched disks enable rapid core assembly that triggers early gas accretion.
  • Inner small planets in such systems retain lower densities while core masses stay comparable.
  • Photoevaporation after formation can mask true core-mass trends in occurrence statistics.
  • The same preference appears when larger radius samples are split across the radius valley.
  • The same preference appears when larger mass samples are split at 10 Earth masses.

Where Pith is reading between the lines

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

  • If confirmed, the trend would imply that disk metallicity sets a timing offset between inner and outer planet growth phases.
  • Future occurrence studies could test whether the density offset disappears in low-metallicity hosts.
  • Models that tie core growth speed directly to solid surface density could be calibrated against the observed density split.

Load-bearing premise

The heterogeneous detection-sensitivity corrections applied to the 43 systems accurately recover the true underlying gas-giant occurrence rate without residual bias.

What would settle it

A larger sample of metal-rich systems, after uniform high-precision sensitivity corrections, shows no preference for lower-density inner planets around those that host gas giants.

Figures

Figures reproduced from arXiv: 2606.12524 by Eve J. Lee, Joseph Y. Tang, Marta L. Bryan.

Figure 2
Figure 2. Figure 2: The median completeness map for all 43 systems. The 13 detected cold gas giants (0.5–20 MJup, 1–10 AU) are plotted as blue dots. Note that all gas giants lie in a region of high sensitivity parameter space for our sample. 4. A SANITY CHECK TO CONFIRM OCCURRENCE RATE P(GG|ISP) FOR OUR SAMPLE With our sample of 43 small planet systems that have both measured masses and radii, we can explore the question of w… view at source ↗
Figure 1
Figure 1. Figure 1: Top: Mass-radius diagram for small planets in￾cluded in our sample. Small planets in systems with and without gas giants are plotted in blue and orange, respec￾tively. We include a density curve for an Earth-like planet in red from L. Zeng et al. (2019). Middle: Stellar mass distributions for our sample, split into systems with gas gi￾ants (blue) and without gas giants (red). Bottom: Stellar metallicity di… view at source ↗
Figure 3
Figure 3. Figure 3: A histogram of planet densities (top), envelope mass fractions (EMFs) (middle), and core masses (bot￾tom) for our sample of small planets. Blue and red colors represent systems with and without gas giants, respectively. The gray zone in the middle panel represents the empty re￾gion below the lowest curve in E. D. Lopez & J. J. Fortney (2014) where no EMF is provided. Instead, all H/He frac￾tions are set to… view at source ↗
Figure 4
Figure 4. Figure 4: Left: Small planet radius vs. mass for our sample, where small planets in systems with outer gas giants are outlined in red. Interpolated H/He envelope fractions from E. D. Lopez & J. J. Fortney (2014) are overplotted, with solid curves corresponding to Fplanet = 0.1F⊕, and dashed lines to planets with 10F⊕ assuming an age of 1 Gyr. The solid red line shows an Earth-like planet composition (L. Zeng et al. … view at source ↗
Figure 5
Figure 5. Figure 5: Conditional occurrence rates P(GG|ISP) split by density (left column), EMF (middle column), and core mass (right column) for the entire sample (top row), metal-rich systems (middle row), and metal-poor systems (bottom row). The shaded region indicates the 68% confidence interval for each distribution. Densities were divided at 1ρ⊕, EMF at a value of 1%, and core mass at 7 M⊕ [PITH_FULL_IMAGE:figures/full_… view at source ↗
Figure 6
Figure 6. Figure 6: Left column: Conditional occurrence rates of outer gas giants in systems with inner super-Earths (SE) or sub-Nep￾tunes (SN). SE/SN are categorized using the C. S. K. Ho & V. Van Eylen (2023) equation 11 whereby SE fall below the valley and SN are above. Top, middle, and bottom panels show this frequency for all, high, and low metallicity systems respectively. Middle column: Separating out planets that fall… view at source ↗
Figure 7
Figure 7. Figure 7: The same histograms as [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
read the original abstract

Gas giants may shape the reservoir of solids and gas in the inner disk in which the small planets assemble. To test this possibility, we collect a sample of 43 exoplanetary systems containing 68 inner small planets (ISP) with both measured masses (1-20 M$_{\oplus}$) and radii (1-4 R$_{\oplus}$). After correcting for heterogeneous individual system sensitivities to distant gas giants, we calculate the gas giant occurrence rate in ISP systems P(GG$|$ISP) as a function of inner small planet density, envelope mass fraction (EMF), and core mass. While we find no significant difference between P(GG$|$ISP) given high/low small planet density, EMF, or core mass, we see hints of a trend when only looking at the metal-rich systems. Despite the substantial limitations due to small sample sizes, we find that gas giants in metal-rich systems are preferentially found with lower density planets with similar core masses. We find consistent hints of trends using larger samples of inner planets with measured radii divided across the radius valley or with measured masses divided across 10 $M_\oplus$. Our result is consistent with more metal-enriched disks catalyzing rapid core assembly and kickstarting the gas accretion early, while the muted difference in the outer giant occurrence rate with respect to core mass may indicate contamination by post-formation photoevaporation.

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

3 major / 1 minor

Summary. The paper assembles a sample of 43 systems hosting 68 inner small planets (1-20 M⊕, 1-4 R⊕) with measured masses and radii. After applying heterogeneous sensitivity corrections for distant gas giants, it computes the conditional occurrence rate P(GG|ISP) versus inner-planet density, envelope mass fraction, and core mass. No overall difference is found, but a post-hoc split into metal-rich hosts yields a reported hint that gas giants prefer lower-density inner planets of comparable core mass; the authors interpret this as evidence that metal-rich disks enable rapid core growth and early gas accretion, with photoevaporation possibly contaminating the core-mass trend. Larger radius- or mass-selected samples are cited as yielding consistent hints.

Significance. If the reported trend survives larger samples and validated sensitivity corrections, it would link outer giant-planet occurrence to inner-planet bulk composition in a metallicity-dependent way, offering a testable signature of disk solid/gas reservoir modification by gas giants. The work supplies no machine-checked proofs or parameter-free derivations, and the small-sample, post-hoc nature limits immediate impact.

major comments (3)
  1. [Abstract] Abstract and main text: the central claim is a 'hint' of a trend only after an overall null result and only within the metal-rich subset; no p-values, bootstrap uncertainties, or quantitative significance are supplied for this subdivision, so the strength of evidence for the stated preference cannot be evaluated.
  2. [Abstract] Abstract and methods: the calculation of P(GG|ISP) after heterogeneous sensitivity corrections is presented without any test or quantification showing that the correction factors are uncorrelated with the binning variables (density, metallicity); residual bias correlated with these quantities would produce the observed pattern without physical origin.
  3. [Abstract] Abstract: the metal-rich trend rests on post-hoc subdivision of an already modest sample of 43 systems; the authors themselves note 'substantial limitations due to small sample sizes,' yet the manuscript offers no power analysis or minimum-sample requirement to support reporting the trend.
minor comments (1)
  1. [Abstract] The abstract states 'we find consistent hints of trends using larger samples' but does not specify how those samples avoid the same sensitivity-correction issues or how they are selected relative to the primary 43-system sample.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed comments. We address each major comment point by point below, providing the strongest honest defense of the manuscript while agreeing to revisions where the concerns are valid and addressable.

read point-by-point responses

  1. Referee: [Abstract] Abstract and main text: the central claim is a 'hint' of a trend only after an overall null result and only within the metal-rich subset; no p-values, bootstrap uncertainties, or quantitative significance are supplied for this subdivision, so the strength of evidence for the stated preference cannot be evaluated.

    Authors: We agree that quantitative measures would allow readers to better assess the evidence. Although the result is presented as a 'hint' due to the exploratory post-hoc nature and small sample, we will add bootstrap resampling to report uncertainties on P(GG|ISP) in the metal-rich subset and compute associated p-values or significance levels in the revised manuscript. revision: yes

  2. Referee: [Abstract] Abstract and methods: the calculation of P(GG|ISP) after heterogeneous sensitivity corrections is presented without any test or quantification showing that the correction factors are uncorrelated with the binning variables (density, metallicity); residual bias correlated with these quantities would produce the observed pattern without physical origin.

    Authors: This is a fair point on potential residual bias. In the revised methods section, we will include explicit quantification (e.g., Spearman rank correlations or scatter plots) between the individual sensitivity correction factors and the binning variables (density, metallicity, core mass). We will discuss any detected correlations and their possible impact on the trends. revision: yes

  3. Referee: [Abstract] Abstract: the metal-rich trend rests on post-hoc subdivision of an already modest sample of 43 systems; the authors themselves note 'substantial limitations due to small sample sizes,' yet the manuscript offers no power analysis or minimum-sample requirement to support reporting the trend.

    Authors: We already emphasize the small-sample limitations in the text. To address the request for justification of reporting the trend, we will add a simple power analysis (e.g., using binomial or bootstrap-based simulations) estimating the sample size required to detect the observed effect size at higher significance, while still framing the current result as a hint. revision: yes

Circularity Check

0 steps flagged

No circularity: direct statistical comparison of occurrence rates

full rationale

The paper assembles an observational sample of 43 systems, applies heterogeneous sensitivity corrections to compute P(GG|ISP), and bins the resulting rates by measured inner-planet density, EMF, and core mass. No equation reduces any reported trend to a quantity defined by a fitted parameter, no self-citation supplies a load-bearing uniqueness theorem or ansatz, and the central result is an empirical pattern extracted from the corrected data rather than a renaming or self-referential construction. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis depends on the accuracy of detection-sensitivity corrections and on the assumption that the 43-system sample remains representative after those corrections; these are standard domain assumptions rather than new entities or free parameters introduced by the paper.

axioms (2)
  • domain assumption Heterogeneous individual system sensitivities to distant gas giants can be corrected to yield an unbiased occurrence rate P(GG|ISP).
    Explicitly invoked when the authors calculate occurrence rates after correction.
  • domain assumption The distinction between metal-rich and metal-poor systems is sufficiently clean that post-hoc subdivision does not introduce selection bias.
    Used when reporting the trend only in the metal-rich subset.

pith-pipeline@v0.9.1-grok · 5781 in / 1547 out tokens · 34076 ms · 2026-06-27T08:00:55.322034+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

42 extracted references · 40 canonical work pages · 3 internal anchors

  1. [1]

    2025, arXiv

    Ahrer, E.-M., Radica, M., Piaulet-Ghorayeb, C., et al. 2025, arXiv. https://arxiv.org/abs/2504.20428

    arXiv 2025

  2. [2]

    2015, Proceedings of the National Academy of Science, 112, 4214, doi: 10.1073/pnas.1423252112

    Batygin, K., & Laughlin, G. 2015, Proceedings of the National Academy of Science, 112, 4214, doi: 10.1073/pnas.1423252112

    work page doi:10.1073/pnas.1423252112 2015

  3. [3]

    S., Naponiello, L., Pezzetta, E., et al

    Bonomo, A. S., Naponiello, L., Pezzetta, E., et al. 2025, A&A, 700, A126, doi: 10.1051/0004-6361/202452523

    work page doi:10.1051/0004-6361/202452523 2025

  4. [4]

    No excess of cold Jupiters in small planet systems

    Bonomo, A. S., Dumusque, X., Massa, A., et al. 2023, Astronomy and Astrophysics, 677, A33, doi: 10.1051/0004-6361/202346211

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

  5. [5]

    L., Knutson, H

    Bryan, M. L., Knutson, H. A., Lee, E. J., et al. 2019, The Astronomical Journal, 157, 52, doi: 10.3847/1538-3881/aaf57f

    work page doi:10.3847/1538-3881/aaf57f 2019

  6. [7]

    L., & Lee, E

    Bryan, M. L., & Lee, E. J. 2024b, The Astrophysical Journal Letters, 968, L25, doi: 10.3847/2041-8213/ad5013

    work page doi:10.3847/2041-8213/ad5013 2041

  7. [8]

    L., & Lee, E

    Bryan, M. L., & Lee, E. J. 2025, Resolving the Super-Earth/Gas Giant Connection in Stellar Mass and Metallicity, arXiv, doi: 10.48550/arXiv.2502.01748

    work page doi:10.48550/arxiv.2502.01748 2025

  8. [9]

    Chachan, Y., & Lee, E. J. 2023, The Astrophysical Journal Letters, 952, L20, doi: 10.3847/2041-8213/ace257

    work page doi:10.3847/2041-8213/ace257 2023

  9. [10]

    A., Knutson, H

    Chachan, Y., Dalba, P. A., Knutson, H. A., et al. 2022, ApJ, 926, 62, doi: 10.3847/1538-4357/ac3ed6

    work page doi:10.3847/1538-4357/ac3ed6 2022

  10. [11]

    Sub-Neptunes Show a Stronger Correlation with Cold Jupiters than Super-Earths Especially in Metal-rich Systems

    Chen, D.-C., Dai, F., Ma, B., Liu, S.-F., & Yu, C. 2026, arXiv e-prints, arXiv:2604.20203, doi: 10.48550/arXiv.2604.20203

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2604.20203 2026

  11. [12]

    I., Chiang, E., & Lee, E

    Dawson, R. I., Chiang, E., & Lee, E. J. 2015, Monthly Notices of the Royal Astronomical Society, 453, 1471, doi: 10.1093/mnras/stv1639

    work page doi:10.1093/mnras/stv1639 2015

  12. [13]

    J., & Petigura, E

    Fulton, B. J., & Petigura, E. A. 2018, The Astronomical Journal, 156, 264, doi: 10.3847/1538-3881/aae828

    work page doi:10.3847/1538-3881/aae828 2018

  13. [14]

    J., Petigura, E

    Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, The Astronomical Journal, 154, 109, doi: 10.3847/1538-3881/aa80eb

    work page doi:10.3847/1538-3881/aa80eb 2017

  14. [15]

    2022, Astronomy and Astrophysics, 658, A176, doi: 10.1051/0004-6361/202141429

    Heidari, N., Boisse, I., Orell-Miquel, J., et al. 2022, Astronomy and Astrophysics, 658, A176, doi: 10.1051/0004-6361/202141429

    work page doi:10.1051/0004-6361/202141429 2022

  15. [16]

    Ho, C. S. K., & Van Eylen, V. 2023, Monthly Notices of the Royal Astronomical Society, 519, 4056, doi: 10.1093/mnras/stac3802

    work page doi:10.1093/mnras/stac3802 2023

  16. [17]

    R., Yalcinkaya, S., Osborn, H

    Kane, S. R., Yalcinkaya, S., Osborn, H. P., et al. 2020, VizieR Online Data Catalog, 516, J/AJ/160/129, doi: 10.26093/cds/vizier.51600129

    work page doi:10.26093/cds/vizier.51600129 2020

  17. [18]

    A., et al

    Lambrechts, M., Morbidelli, A., Jacobson, S. A., et al. 2019, A&A, 627, A83, doi: 10.1051/0004-6361/201834229

    work page doi:10.1051/0004-6361/201834229 2019

  18. [19]

    Lee, E. J. 2019, ApJ, 878, 36, doi: 10.3847/1538-4357/ab1b40

    work page doi:10.3847/1538-4357/ab1b40 2019

  19. [20]

    , keywords =

    Lee, E. J., & Connors, N. J. 2021, ApJ, 908, 32, doi: 10.3847/1538-4357/abd6c7

    work page doi:10.3847/1538-4357/abd6c7 2021

  20. [21]

    , keywords =

    Lee, E. J., Karalis, A., & Thorngren, D. P. 2022, ApJ, 941, 186, doi: 10.3847/1538-4357/ac9c66

    work page doi:10.3847/1538-4357/ac9c66 2022

  21. [22]

    W., Lee, E

    Lin, J. W., Lee, E. J., & Chiang, E. 2018, MNRAS, 480, 4338, doi: 10.1093/mnras/sty2159

    work page doi:10.1093/mnras/sty2159 2018

  22. [23]

    Liu, B., & Ormel, C. W. 2018, A&A, 615, A138, doi: 10.1051/0004-6361/201732307

    work page doi:10.1051/0004-6361/201732307 2018

  23. [24]

    D., & Fortney, J

    Lopez, E. D., & Fortney, J. J. 2014, The Astrophysical Journal, 792, 1, doi: 10.1088/0004-637X/792/1/1

    work page doi:10.1088/0004-637x/792/1/1 2014

  24. [25]

    2023, Astronomy and Astrophysics, 680, A84, doi: 10.1051/0004-6361/202347744

    Martioli, E., H´ ebrard, G., de Almeida, L., et al. 2023, Astronomy and Astrophysics, 680, A84, doi: 10.1051/0004-6361/202347744

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

  25. [26]

    2023, AJ, 165, 265, doi: 10.3847/1538-3881/acd175

    Ment, K., & Charbonneau, D. 2023, AJ, 165, 265, doi: 10.3847/1538-3881/acd175

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

  26. [27]

    2015, Icarus, 258, 418, doi: 10.1016/j.icarus.2015.06.003

    Morbidelli, A., Lambrechts, M., Jacobson, S., & Bitsch, B. 2015, Icarus, 258, 418, doi: 10.1016/j.icarus.2015.06.003

    work page doi:10.1016/j.icarus.2015.06.003 2015

  27. [28]

    J., et al

    Ohno, K., Schlawin, E., Bell, T. J., et al. 2025, ApJL, 979, L7, doi: 10.3847/2041-8213/ada02c

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

  28. [29]

    Ormel, C. W. 2017, in Astrophysics and Space Science

    2017

  29. [30]

    445, Formation, Evolution, and Dynamics of Young Solar Systems, ed

    Library, Vol. 445, Formation, Evolution, and Dynamics of Young Solar Systems, ed. M. Pessah & O. Gressel, 197, doi: 10.1007/978-3-319-60609-5 7

    work page doi:10.1007/978-3-319-60609-5

  30. [31]

    , keywords =

    Owen, J. E., & Jackson, A. P. 2012, MNRAS, 425, 2931, doi: 10.1111/j.1365-2966.2012.21481.x

    work page doi:10.1111/j.1365-2966.2012.21481.x 2012

  31. [32]

    Kepler planets: a tale of evaporation

    Owen, J. E., & Wu, Y. 2013, The Astrophysical Journal, 775, 105, doi: 10.1088/0004-637X/775/2/105

    work page internal anchor Pith review doi:10.1088/0004-637x/775/2/105 2013

  32. [33]

    E., & Wu, Y

    Owen, J. E., & Wu, Y. 2017, ApJ, 847, 29, doi: 10.3847/1538-4357/aa890a

    work page doi:10.3847/1538-4357/aa890a 2017

  33. [34]

    , keywords =

    Parc, L., Bouchy, F., Venturini, J., Dorn, C., & Helled, R. 2024, A&A, 688, A59, doi: 10.1051/0004-6361/202449911

    work page doi:10.1051/0004-6361/202449911 2024

  34. [35]

    2022, A&A, 664, A65, doi: 10.1051/0004-6361/202142828

    Pinamonti, M., Sozzetti, A., Maldonado, J., et al. 2022, A&A, 664, A65, doi: 10.1051/0004-6361/202142828

    work page doi:10.1051/0004-6361/202142828 2022

  35. [36]

    2025, arXiv e-prints, arXiv:2510.11397, doi: 10.48550/arXiv.2510.11397

    Rochon, A., Artigau, ´E., Weisserman, D., et al. 2025, arXiv e-prints, arXiv:2510.11397, doi: 10.48550/arXiv.2510.11397

    work page doi:10.48550/arxiv.2510.11397 2025

  36. [37]

    J., Fulton, B

    Rosenthal, L. J., Fulton, B. J., Hirsch, L. A., et al. 2021, The Astrophysical Journal Supplement Series, 255, 8, doi: 10.3847/1538-4365/abe23c Van Eylen, V., Agentoft, C., Lundkvist, M. S., et al. 2018, MNRAS, 479, 4786, doi: 10.1093/mnras/sty1783

    work page doi:10.3847/1538-4365/abe23c 2021

  37. [38]

    J., Morbidelli, A., Raymond, S

    Walsh, K. J., Morbidelli, A., Raymond, S. N., O’Brien, D. P., & Mandell, A. M. 2011, Nature, 475, 206, doi: 10.1038/nature10201

    work page internal anchor Pith review doi:10.1038/nature10201 2011

  38. [39]

    M., Isaacson, H., Howard, A

    Weiss, L. M., Isaacson, H., Howard, A. W., et al. 2024, The Astrophysical Journal Supplement Series, 270, 8, doi: 10.3847/1538-4365/ad0cab 18

    work page doi:10.3847/1538-4365/ad0cab 2024

  39. [40]

    B., Sasselov, D

    Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, Proceedings of the National Academy of Sciences, 116, 9723, doi: 10.1073/pnas.1812905116

    work page doi:10.1073/pnas.1812905116 2019

  40. [41]

    2024, Research in Astronomy and Astrophysics, 24, 045013, doi: 10.1088/1674-4527/ad3132

    Zhu, W. 2024, Research in Astronomy and Astrophysics, 24, 045013, doi: 10.1088/1674-4527/ad3132

    work page doi:10.1088/1674-4527/ad3132 2024

  41. [42]

    2018, ApJ, 860, 101, doi: 10.3847/1538-4357/aac6d5

    Zhu, W., Petrovich, C., Wu, Y., Dong, S., & Xie, J. 2018, ApJ, 860, 101, doi: 10.3847/1538-4357/aac6d5

    work page doi:10.3847/1538-4357/aac6d5 2018

  42. [43]

    2018, AJ, 156, 92, doi: 10.3847/1538-3881/aad22a

    Zhu, W., & Wu, Y. 2018, The Astronomical Journal, 156, 92, doi: 10.3847/1538-3881/aad22a

    work page doi:10.3847/1538-3881/aad22a 2018