Recognition: no theorem link
The Intrinsic Multiplicity Distribution of Exoplanets Revealed from the Radial Velocity Method. II. Constraints on Giant Planet Multiplicity from Different Surveys
Pith reviewed 2026-05-15 17:11 UTC · model grok-4.3
The pith
Giant planet multiplicity around Sun-like stars falls short of planet-planet scattering predictions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
From the HARPS sample, (7.8±1.4%, 2.3±1.2%, 0.5+0.8-0.3%) of Sun-like stars have one, two, or three giant planets within 10 au; the California Legacy Survey yields (7.3±2.8%, 7.2±2.3%, <1.3%, 1.0+1.0-0.6%) for one through four. After correcting for detection completeness, the total fraction of stars with at least one giant planet is 10.6±1.2% and 15.8±2.1% respectively, and this distribution is inconsistent with most planet-planet scattering models, which predict too many multi-giant systems.
What carries the argument
The intrinsic multiplicity distribution recovered by applying detection completeness corrections to cleaned radial velocity survey samples of giant planets above Saturn mass.
If this is right
- Scattering is disfavored as the dominant process shaping giant planet systems after formation.
- Most giant planets form and remain in single-planet configurations rather than packed multiples.
- Metallicity primarily affects overall occurrence rates rather than the shape of the multiplicity distribution.
- Theoretical models must be revised to produce mostly isolated giant planets within 10 au.
Where Pith is reading between the lines
- The pattern supports core-accretion formation with limited subsequent dynamical stirring for most systems.
- Extending the same completeness-corrected approach to sub-Saturn planets could test whether the single-planet dominance holds at lower masses.
- Independent constraints from transit or direct-imaging surveys on outer giants would provide a cross-check on the 10 au cutoff.
Load-bearing premise
The cleaned RV survey samples represent Sun-like stars without bias once detection completeness and metallicity differences are accounted for.
What would settle it
A new volume-complete RV survey that directly counts all giant planets above Saturn mass within 10 au around Sun-like stars and finds a substantially higher fraction of multi-planet systems than the reported 2-7 percent would falsify the inferred distribution.
Figures
read the original abstract
Compared to the commonly used planet occurrence rates, the multiplicity distribution of planets can be more useful in constraining the formation and evolution pathways of planetary systems. This work follows an earlier work of Zhu (2022) and derive the intrinsic multiplicity distribution of giant planets (with masses above Saturn mass) from two independent radial velocity (RV) surveys. In particular, we find that $(7.8\pm1.4\%, 2.3\pm1.2\%, 0.5^{+0.8}_{-0.3}\%)$ of Sun-like stars in the HARPS sample have $(1, 2, 3)$ giant planets within 10\,au, whereas $(7.3\pm2.8\%, 7.2\pm2.3\%, <1.3\%, 1.0^{+1.0}_{-0.6}\%)$ of Sun-like stars in the California Legacy Survey (CLS) have $(1, 2, 3, 4)$ giant planets within 10\,au. Here we have further cleaned the CLS sample and removed planet detections that were not discovered in the survey mode. The total fraction of Sun-like stars with at least one giant planet within 10\,au from the two samples are $10.6\pm1.2\%$ and $15.8\pm2.1\%$, respectively, and the difference may be accounted for by their different metallicity distributions. We briefly discuss the theoretical implications of our results. In particular, the inferred giant planet multiplicity distribution is inconsistent with most of the proposed theoretical models involving planet--planet scatterings, which predict too many multi-giant systems.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives the intrinsic multiplicity distribution of giant planets (masses above Saturn mass) within 10 au from two independent radial velocity surveys (HARPS and the cleaned California Legacy Survey). It reports specific fractions of Sun-like stars hosting 1, 2, 3 (and 4 for CLS) such planets, with totals of 10.6±1.2% and 15.8±2.1% respectively, attributing sample differences to metallicity distributions, and concludes that the low multiplicity is inconsistent with most planet-planet scattering models, which overpredict multi-giant systems.
Significance. If the completeness corrections and sample cleaning prove robust, the quantitative multiplicity fractions provide stronger constraints on giant planet formation and evolution than occurrence rates alone, as they directly test predictions from scattering scenarios. The use of two independent surveys and explicit uncertainty reporting adds value for theoretical comparisons, though the result hinges on unverified aspects of the inversion method.
major comments (2)
- [Abstract and methods] Abstract and methods section: the reported fractions (e.g., HARPS 7.8±1.4% single, 2.3±1.2% double) are presented as the outcome of fitting observed detections to an intrinsic distribution, but no explicit equations, likelihood function, or step-by-step completeness modeling are shown to demonstrate how detection biases for 2+ planet systems within 10 au (including metallicity effects) are inverted without circularity.
- [CLS sample description] CLS sample description: the cleaning step that removes non-survey-mode detections is load-bearing for the claim of low multi-giant fractions (7.3% single, 7.2% double), yet no quantitative test is provided on whether the cleaned sample remains representative of Sun-like stars or whether completeness for higher-multiplicity systems is systematically underestimated due to unmodeled interactions or noise.
minor comments (1)
- The abstract would benefit from a one-sentence outline of the statistical method used to recover the intrinsic distribution from the surveys.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight areas where the methods and sample validation can be clarified, which will strengthen the manuscript. We address each major comment below and will revise the paper accordingly.
read point-by-point responses
-
Referee: [Abstract and methods] Abstract and methods section: the reported fractions (e.g., HARPS 7.8±1.4% single, 2.3±1.2% double) are presented as the outcome of fitting observed detections to an intrinsic distribution, but no explicit equations, likelihood function, or step-by-step completeness modeling are shown to demonstrate how detection biases for 2+ planet systems within 10 au (including metallicity effects) are inverted without circularity.
Authors: We agree that additional explicit details would improve clarity. The inversion follows the Bayesian framework of Zhu (2022), in which the observed detection counts are modeled as the intrinsic multiplicity distribution convolved with survey-specific completeness functions derived from injection-recovery tests. The likelihood is a Poisson likelihood on the binned detections, with the completeness for multi-planet systems computed by simulating the joint detectability of multiple signals (accounting for period and mass correlations). Metallicity dependence enters only through the occurrence prior, not the completeness itself, avoiding circularity. We will add the explicit likelihood expression, a flowchart of the completeness calculation, and a dedicated subsection in Methods describing these steps. revision: yes
-
Referee: [CLS sample description] CLS sample description: the cleaning step that removes non-survey-mode detections is load-bearing for the claim of low multi-giant fractions (7.3% single, 7.2% double), yet no quantitative test is provided on whether the cleaned sample remains representative of Sun-like stars or whether completeness for higher-multiplicity systems is systematically underestimated due to unmodeled interactions or noise.
Authors: The cleaning removes detections obtained outside the uniform survey cadence to maintain consistent selection and observing strategy. We will add a table comparing the distributions of stellar mass, metallicity, and spectral type before and after cleaning, confirming that the cleaned sample remains representative of the original Sun-like star population. For completeness, our injection tests already incorporate realistic RV noise and sampling; however, planet-planet dynamical interactions are not explicitly simulated in the recovery step. We will add a quantitative estimate of the potential bias (based on N-body simulations of a subset of systems) and include this as a systematic uncertainty in the final error budget. revision: partial
Circularity Check
Minor self-citation to Zhu (2022) method paper; central inference is data-driven fitting with no reduction by construction
specific steps
-
self citation load bearing
[Abstract]
"This work follows an earlier work of Zhu (2022) and derive the intrinsic multiplicity distribution of giant planets (with masses above Saturn mass) from two independent radial velocity (RV) surveys."
The core inference procedure is imported from the overlapping-author prior paper; while the current application to new samples yields independent numerical results, the self-citation is the sole justification for the fitting methodology used to obtain the reported fractions.
full rationale
The paper applies the multiplicity inference framework from the author's prior work (Zhu 2022) to cleaned HARPS and CLS RV data, fitting assumed intrinsic fractions (1-planet, 2-planet, etc.) against observed detections after applying survey-specific completeness corrections. This is standard forward modeling; the output fractions are not equivalent to any prior fitted constants or self-cited inputs by construction. The self-citation supports the method but is not load-bearing for the new numerical results or the comparison to scattering models. No self-definitional loops, fitted inputs renamed as predictions, or ansatz smuggling occur. The derivation remains self-contained against the external RV catalogs and completeness functions.
Axiom & Free-Parameter Ledger
free parameters (1)
- multiplicity fractions per survey
axioms (2)
- domain assumption RV detection completeness is known and can be inverted to recover intrinsic occurrence rates
- domain assumption Cleaned survey samples are unbiased representatives of Sun-like stars
Reference graph
Works this paper leans on
-
[1]
Bean, J. L., McArthur, B. E., Benedict, G. F., et al. 2007, AJ, 134, 749, doi: 10.1086/519956
-
[2]
2023, A&A, 674, A178, doi: 10.1051/0004-6361/202245040
Bitsch, B., & Izidoro, A. 2023, A&A, 674, A178, doi: 10.1051/0004-6361/202245040
-
[3]
2020, A&A, 643, A66, doi: 10.1051/0004-6361/202038856
Bitsch, B., Trifonov, T., & Izidoro, A. 2020, A&A, 643, A66, doi: 10.1051/0004-6361/202038856
-
[4]
2025, ApJL, 995, L43, doi: 10.3847/2041-8213/ae2603
Brahm, R., Trifonov, T., Jord´ an, A., et al. 2025, ApJL, 995, L43, doi: 10.3847/2041-8213/ae2603
-
[5]
Bryan, M. L., Knutson, H. A., Howard, A. W., et al. 2016, ApJ, 821, 89, doi: 10.3847/0004-637X/821/2/89
-
[6]
Butler, R. P., Vogt, S. S., Marcy, G. W., et al. 2000, ApJ, 545, 504, doi: 10.1086/317796
-
[7]
Carrera, D., Raymond, S. N., & Davies, M. B. 2019, A&A, 629, L7, doi: 10.1051/0004-6361/201935744
-
[8]
Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580, doi: 10.1086/590227
-
[9]
Christiansen, J. L., McElroy, D. L., Harbut, M., et al. 2025, PSJ, 6, 186, doi: 10.3847/PSJ/ade3c2
-
[10]
Cumming, A., Butler, R. P., Marcy, G. W., et al. 2008, Publications of the Astronomical Society of the Pacific, 120, 531554, doi: 10.1086/588487
-
[11]
doi:10.1146/annurev-astro-081817-051853 , eprint =
Dawson, R. I., & Johnson, J. A. 2018, ARA&A, 56, 175, doi: 10.1146/annurev-astro-081817-051853
work page internal anchor Pith review doi:10.1146/annurev-astro-081817-051853 2018
-
[12]
Duffell, P. C., & Chiang, E. 2015, ApJ, 812, 94, doi: 10.1088/0004-637X/812/2/94
-
[13]
2021, A&A, 656, A69, doi: 10.1051/0004-6361/202038553
Emsenhuber, A., Mordasini, C., Burn, R., et al. 2021, A&A, 656, A69, doi: 10.1051/0004-6361/202038553
-
[14]
Fernandes, R. B., Mulders, G. D., Pascucci, I., Mordasini, C., & Emsenhuber, A. 2019, The Astrophysical Journal, 874, 81, doi: 10.3847/1538-4357/ab0300
-
[15]
Fischer, D. A., Marcy, G. W., Butler, R. P., et al. 2001, ApJ, 551, 1107, doi: 10.1086/320224
-
[16]
Fischer, D. A., & Valenti, J. 2005, The Astrophysical Journal, 622, 1102, doi: 10.1086/428383
-
[17]
Ford, E. B., & Rasio, F. A. 2008, ApJ, 686, 621, doi: 10.1086/590926
-
[18]
Foreman-Mackey, D., Hogg, D. W., & Morton, T. D. 2014, ApJ, 795, 64, doi: 10.1088/0004-637X/795/1/64
-
[19]
Frelikh, R., Jang, H., Murray-Clay, R. A., & Petrovich, C. 2019, ApJL, 884, L47, doi: 10.3847/2041-8213/ab4a7b
-
[20]
Fulton, B. J., Rosenthal, L. J., Hirsch, L. A., et al. 2021, ApJS, 255, 14, doi: 10.3847/1538-4365/abfcc1
-
[21]
Gan, T., Wang, S. X., Wang, S., et al. 2023, AJ, 165, 17, doi: 10.3847/1538-3881/ac9b12
-
[22]
He, M. Y., Ford, E. B., Ragozzine, D., & Carrera, D. 2020, AJ, 160, 276, doi: 10.3847/1538-3881/abba18 Giant Planet Intrinsic Multiplicities from HARPS & CLS11
-
[23]
Ida, S., Lin, D. N. C., & Nagasawa, M. 2013, ApJ, 775, 42, doi: 10.1088/0004-637X/775/1/42 Juri´ c, M., & Tremaine, S. 2008, ApJ, 686, 603, doi: 10.1086/590047
-
[24]
Knutson, H. A., Fulton, B. J., Montet, B. T., et al. 2014, ApJ, 785, 126, doi: 10.1088/0004-637X/785/2/126
-
[25]
2023, A&A, 677, A71, doi: 10.1051/0004-6361/202346165
Lagrange, A.-M., Philipot, F., Rubini, P., et al. 2023, A&A, 677, A71, doi: 10.1051/0004-6361/202346165
-
[26]
Lee, M. H., & Peale, S. J. 2002, ApJ, 567, 596, doi: 10.1086/338504
-
[27]
Lin, D. N. C., & Ida, S. 1997, ApJ, 477, 781, doi: 10.1086/303738
-
[28]
J., Ragozzine, D., Fabrycky, D
Lissauer, J. J., Ragozzine, D., Fabrycky, D. C., et al. 2011, ApJS, 197, 8, doi: 10.1088/0067-0049/197/1/8
-
[29]
2004, A&A, 415, 391, doi: 10.1051/0004-6361:20034250
Mayor, M., Udry, S., Naef, D., et al. 2004, A&A, 415, 391, doi: 10.1051/0004-6361:20034250
-
[30]
Mayor, M., Marmier, M., Lovis, C., et al. 2011, arXiv e-prints, arXiv:1109.2497. https://arxiv.org/abs/1109.2497
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[31]
Morbidelli, A., & Raymond, S. N. 2016, Journal of Geophysical Research (Planets), 121, 1962, doi: 10.1002/2016JE005088
-
[32]
D., Pascucci, I., Apai, D., & Ciesla, F
Mulders, G. D., Pascucci, I., Apai, D., & Ciesla, F. J. 2018, AJ, 156, 24, doi: 10.3847/1538-3881/aac5ea
-
[33]
2024, ApJ, 969, 133, doi: 10.3847/1538-4357/ad3046
Nagpal, V., Goldberg, M., & Batygin, K. 2024, ApJ, 969, 133, doi: 10.3847/1538-4357/ad3046
-
[34]
Noyes, R. W., Jha, S., Korzennik, S. G., et al. 1997, ApJL, 483, L111, doi: 10.1086/310754
-
[35]
2023, in Astronomical Society of the Pacific Conference Series, Vol
Paardekooper, S., Dong, R., Duffell, P., et al. 2023, in Astronomical Society of the Pacific Conference Series, Vol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 685, doi: 10.48550/arXiv.2203.09595
-
[36]
2025, ApJ, 985, 7, doi: 10.3847/1538-4357/adc7a9
Pan, M., Liu, B., Jiang, L., et al. 2025, ApJ, 985, 7, doi: 10.3847/1538-4357/adc7a9
-
[37]
Petigura, E. A., Marcy, G. W., Winn, J. N., et al. 2018, AJ, 155, 89, doi: 10.3847/1538-3881/aaa54c
-
[38]
Rasio, F. A., & Ford, E. B. 1996, Science, 274, 954, doi: 10.1126/science.274.5289.954
-
[39]
Rosenthal, L. J., Howard, A. W., Knutson, H. A., & Fulton, B. J. 2024, ApJS, 270, 1, doi: 10.3847/1538-4365/acffc0
-
[40]
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
-
[41]
Santos, N. C., Israelian, G., & Mayor, M. 2001a, A&A, 373, 1019, doi: 10.1051/0004-6361:20010648
-
[42]
Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153, doi: 10.1051/0004-6361:20034469
-
[43]
C., Mayor, M., Naef, D., et al
Santos, N. C., Mayor, M., Naef, D., et al. 2001b, A&A, 379, 999, doi: 10.1051/0004-6361:20011366
-
[44]
Tinney, C. G., Butler, R. P., Marcy, G. W., et al. 2001, ApJ, 551, 507, doi: 10.1086/320097
-
[45]
2012, AJ, 143, 94, doi: 10.1088/0004-6256/143/4/94
Tremaine, S., & Dong, S. 2012, AJ, 143, 94, doi: 10.1088/0004-6256/143/4/94
- [46]
-
[47]
Weidenschilling, S. J., & Marzari, F. 1996, Nature, 384, 619, doi: 10.1038/384619a0
-
[48]
doi:10.1146/annurev-astro-082214-122246 , eprint =
Winn, J. N., & Fabrycky, D. C. 2015, ARA&A, 53, 409, doi: 10.1146/annurev-astro-082214-122246
-
[49]
Wittenmyer, R. A., Wang, S., Horner, J., et al. 2020, MNRAS, 492, 377, doi: 10.1093/mnras/stz3436
-
[50]
Wright, J. T. 2005, PASP, 117, 657, doi: 10.1086/430369
-
[51]
Wright, J. T., Marcy, G. W., Howard, A. W., et al. 2012, The Astrophysical Journal, 753, 160, doi: 10.1088/0004-637X/753/2/160
-
[52]
2011, ApJ, 735, 109, doi: 10.1088/0004-637X/735/2/109
Wu, Y., & Lithwick, Y. 2011, ApJ, 735, 109, doi: 10.1088/0004-637X/735/2/109
-
[53]
2019, ApJ, 873, 8, doi: 10.3847/1538-4357/ab0205
Zhu, W. 2019, ApJ, 873, 8, doi: 10.3847/1538-4357/ab0205
-
[54]
2022, AJ, 164, 5, doi: 10.3847/1538-3881/ac6f59
Zhu, W. 2022, AJ, 164, 5, doi: 10.3847/1538-3881/ac6f59
-
[55]
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
-
[56]
2021, ARA&A, 59, 291, doi: 10.1146/annurev-astro-112420-020055
Zhu, W., & Dong, S. 2021, ARA&A, 59, 291, doi: 10.1146/annurev-astro-112420-020055
-
[57]
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
-
[58]
2018, AJ, 156, 92, doi: 10.3847/1538-3881/aad22a
Zhu, W., & Wu, Y. 2018, AJ, 156, 92, doi: 10.3847/1538-3881/aad22a
-
[59]
Zink, J. K., & Howard, A. W. 2023, ApJL, 956, L29, doi: 10.3847/2041-8213/acfdab
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.