pith. machine review for the scientific record. sign in

arxiv: 2603.01534 · v2 · submitted 2026-03-02 · 🌌 astro-ph.EP

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

Authors on Pith no claims yet

Pith reviewed 2026-05-15 17:11 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords giant planetsexoplanet multiplicityradial velocityplanet formationHARPSCalifornia Legacy SurveySun-like starsoccurrence rates
0
0 comments X

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.

This paper derives the intrinsic multiplicity distribution of giant planets from two cleaned radial velocity surveys. It reports that about 8 percent of Sun-like stars host exactly one giant planet within 10 au, with 2 to 7 percent hosting two and even smaller fractions hosting three or more. The low rate of multiples differs from theoretical models based on planet-planet scatterings, which produce too many systems with several giants. The difference between the two surveys is attributed to their distinct stellar metallicity distributions, and the results constrain formation pathways by showing that most giant planets occur singly.

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

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

  • 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

Figures reproduced from arXiv: 2603.01534 by Jiayin Li, Wei Zhu.

Figure 1
Figure 1. Figure 1: The distribution of 822 HARPS sample stars in the minimum mass vs. semi–major axis plane. Planets with different observed multiplicities are differentiated with different labels and colors. The contours represent the sur￾vey completeness, with lighter shades indicating higher sen￾sitivity and darker shades indicating lower sensitivity. Both based on the data from M. Mayor et al. (2011). The vertical dashed… view at source ↗
Figure 3
Figure 3. Figure 3: The distribution of the remaining 351 CLS sam￾ple stars in the minimum mass vs. semimajor axis plane after selection. The meanings of the markers, background contours, and various lines are the same as described in the caption of [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: This figure shows in blue dots the planets de￾tected in the less biased CLS sample whose properties fall within our chosen parameter space. The values shown in the grid and their meanings are the same as described in the caption of [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Illustrations of the derived intrinsic planet multiplicity distributions for the three chosen planet classes from the cleaned CLS (blue) and HARPS (orange) datasets. The maximum-likelihood solutions and 1σ–3σ confidence intervals are indicated by the colored bars with increasing transparency. np mp Fp 0.0 0.5 1.0 1.5 2.0 CLS / HARPS [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Ratios of CLS to HARPS values for average num￾ber of planets per star (¯np), average multiplicity ( ¯mp), and the fraction of stars with planets (Fp). The values are taken for the “All giant” population. Error bars indicate 1σ un￾certainties. The dashed line at unity marks equal values between the two samples. even more metal-rich environment than do cold giants, thus putting constraints on the in situ for… view at source ↗
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.

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

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

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

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

1 steps flagged

Minor self-citation to Zhu (2022) method paper; central inference is data-driven fitting with no reduction by construction

specific steps
  1. 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

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on statistical inference from RV detections that assumes standard completeness models and representative samples after cleaning; no new entities are introduced.

free parameters (1)
  • multiplicity fractions per survey
    Fitted values (e.g., 7.8%, 2.3%) chosen to match observed planet counts after completeness correction.
axioms (2)
  • domain assumption RV detection completeness is known and can be inverted to recover intrinsic occurrence rates
    Invoked to convert observed detections into intrinsic multiplicity percentages.
  • domain assumption Cleaned survey samples are unbiased representatives of Sun-like stars
    Required for generalizing the fractions to the broader population.

pith-pipeline@v0.9.0 · 5611 in / 1260 out tokens · 47564 ms · 2026-05-15T17:11:11.665903+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

59 extracted references · 59 canonical work pages · 2 internal anchors

  1. [1]

    L., McArthur, B

    Bean, J. L., McArthur, B. E., Benedict, G. F., et al. 2007, AJ, 134, 749, doi: 10.1086/519956

  2. [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. [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. [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. [5]

    L., Knutson, H

    Bryan, M. L., Knutson, H. A., Howard, A. W., et al. 2016, ApJ, 821, 89, doi: 10.3847/0004-637X/821/2/89

  6. [6]

    P., Vogt, S

    Butler, R. P., Vogt, S. S., Marcy, G. W., et al. 2000, ApJ, 545, 504, doi: 10.1086/317796

  7. [7]

    N., & Davies, M

    Carrera, D., Raymond, S. N., & Davies, M. B. 2019, A&A, 629, L7, doi: 10.1051/0004-6361/201935744

  8. [8]

    B., Matsumura, S., & Rasio, F

    Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580, doi: 10.1086/590227

  9. [9]

    L., McElroy, D

    Christiansen, J. L., McElroy, D. L., Harbut, M., et al. 2025, PSJ, 6, 186, doi: 10.3847/PSJ/ade3c2

  10. [10]

    P., Marcy, G

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

  12. [12]

    C., & Chiang, E

    Duffell, P. C., & Chiang, E. 2015, ApJ, 812, 94, doi: 10.1088/0004-637X/812/2/94

  13. [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. [14]

    B., Mulders, G

    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. [15]

    A., Marcy, G

    Fischer, D. A., Marcy, G. W., Butler, R. P., et al. 2001, ApJ, 551, 1107, doi: 10.1086/320224

  16. [16]

    A., & Valenti, J

    Fischer, D. A., & Valenti, J. 2005, The Astrophysical Journal, 622, 1102, doi: 10.1086/428383

  17. [17]

    B., & Rasio, F

    Ford, E. B., & Rasio, F. A. 2008, ApJ, 686, 621, doi: 10.1086/590926

  18. [18]

    W., & Morton, T

    Foreman-Mackey, D., Hogg, D. W., & Morton, T. D. 2014, ApJ, 795, 64, doi: 10.1088/0004-637X/795/1/64

  19. [19]

    A., & Petrovich, C

    Frelikh, R., Jang, H., Murray-Clay, R. A., & Petrovich, C. 2019, ApJL, 884, L47, doi: 10.3847/2041-8213/ab4a7b

  20. [20]

    J., Rosenthal, L

    Fulton, B. J., Rosenthal, L. J., Hirsch, L. A., et al. 2021, ApJS, 255, 14, doi: 10.3847/1538-4365/abfcc1

  21. [21]

    X., Wang, S., et al

    Gan, T., Wang, S. X., Wang, S., et al. 2023, AJ, 165, 17, doi: 10.3847/1538-3881/ac9b12

  22. [22]

    Y., Ford, E

    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. [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. [24]

    A., Fulton, B

    Knutson, H. A., Fulton, B. J., Montet, B. T., et al. 2014, ApJ, 785, 126, doi: 10.1088/0004-637X/785/2/126

  25. [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. [26]

    H., & Peale, S

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

  27. [27]

    Lin, D. N. C., & Ida, S. 1997, ApJ, 477, 781, doi: 10.1086/303738

  28. [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. [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. [30]

    The HARPS search for southern extra-solar planets XXXIV. Occurrence, mass distribution and orbital properties of super-Earths and Neptune-mass planets

    Mayor, M., Marmier, M., Lovis, C., et al. 2011, arXiv e-prints, arXiv:1109.2497. https://arxiv.org/abs/1109.2497

  31. [31]

    Morbidelli, A., & Raymond, S. N. 2016, Journal of Geophysical Research (Planets), 121, 1962, doi: 10.1002/2016JE005088

  32. [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. [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. [34]

    W., Jha, S., Korzennik, S

    Noyes, R. W., Jha, S., Korzennik, S. G., et al. 1997, ApJL, 483, L111, doi: 10.1086/310754

  35. [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. [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. [37]

    A., Marcy, G

    Petigura, E. A., Marcy, G. W., Winn, J. N., et al. 2018, AJ, 155, 89, doi: 10.3847/1538-3881/aaa54c

  38. [38]

    A., & Ford, E

    Rasio, F. A., & Ford, E. B. 1996, Science, 274, 954, doi: 10.1126/science.274.5289.954

  39. [39]

    J., Howard, A

    Rosenthal, L. J., Howard, A. W., Knutson, H. A., & Fulton, B. J. 2024, ApJS, 270, 1, doi: 10.3847/1538-4365/acffc0

  40. [40]

    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

  41. [41]

    C., Israelian, G., & Mayor, M

    Santos, N. C., Israelian, G., & Mayor, M. 2001a, A&A, 373, 1019, doi: 10.1051/0004-6361:20010648

  42. [42]

    C., Israelian, G., & Mayor, M

    Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153, doi: 10.1051/0004-6361:20034469

  43. [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. [44]

    G., Butler, R

    Tinney, C. G., Butler, R. P., Marcy, G. W., et al. 2001, ApJ, 551, 507, doi: 10.1086/320097

  45. [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. [46]

    2000, A&A, 356, 590

    Udry, S., Mayor, M., Naef, D., et al. 2000, A&A, 356, 590

  47. [47]

    J., & Marzari, F

    Weidenschilling, S. J., & Marzari, F. 1996, Nature, 384, 619, doi: 10.1038/384619a0

  48. [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. [49]

    , keywords =

    Wittenmyer, R. A., Wang, S., Horner, J., et al. 2020, MNRAS, 492, 377, doi: 10.1093/mnras/stz3436

  50. [50]

    Wright, J. T. 2005, PASP, 117, 657, doi: 10.1086/430369

  51. [51]

    T., Marcy, G

    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. [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. [53]

    2019, ApJ, 873, 8, doi: 10.3847/1538-4357/ab0205

    Zhu, W. 2019, ApJ, 873, 8, doi: 10.3847/1538-4357/ab0205

  54. [54]

    2022, AJ, 164, 5, doi: 10.3847/1538-3881/ac6f59

    Zhu, W. 2022, AJ, 164, 5, doi: 10.3847/1538-3881/ac6f59

  55. [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. [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. [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. [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. [59]

    K., & Howard, A

    Zink, J. K., & Howard, A. W. 2023, ApJL, 956, L29, doi: 10.3847/2041-8213/acfdab