Upper limits on protolunar disc masses using ALMA observations of directly-imaged exoplanets
Pith reviewed 2026-05-25 13:59 UTC · model grok-4.3
The pith
ALMA non-detections of four young directly-imaged exoplanets place upper limits on circumplanetary dust below the mass needed to form regular satellite systems.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The non-detections of 1.3 mm continuum emission around PZ Tel, AB Pic, 51 Eri, and κ And imply dust masses in their circumplanetary discs below 10^{-4} times the planet mass for both extended (0.4 Hill radius) and compact (0.04 Hill radius) disc models, which the authors take as evidence that dust has grown beyond metre sizes to form moonetesimals within less than 10 Myr.
What carries the argument
Radiative transfer models of protolunar discs with two representative sizes (0.4 and 0.04 times the exoplanet's Hill radius) that convert the ALMA non-detections into upper limits on dust mass.
If this is right
- The upper limits are lower than the dust mass required to explain regular satellite systems.
- Dust growth to moonetesimals occurs in timescales under 10 Myr.
- The sample of ALMA non-detections of young companions increases by 50%.
- Comparison with viscous evolution and debris disc models supports the interpretation of rapid solid growth.
Where Pith is reading between the lines
- Similar non-detections in additional systems could indicate that moon formation around giant planets is generally complete by 10 Myr.
- If future ALMA observations detect dust above these limits in comparable planets, it would point to longer-lived primordial dust or different disc conditions in those cases.
- These upper limits may suggest that the assembly of regular satellites around extrasolar gas giants follows a timeline comparable to the Solar System's Galilean moons.
Load-bearing premise
The two chosen disc sizes and standard dust properties adequately represent all possible circumplanetary material distributions, and that 10^{-4} planet mass is the correct dust threshold for forming regular satellite systems.
What would settle it
Detection of 1.3 mm emission from any of these four systems at a level corresponding to more than 10^{-4} planet mass of dust, or finding a young directly-imaged planet with sufficient dust to match the threshold without evidence of moon formation.
Figures
read the original abstract
The Solar System gas giants are each surrounded by many moons, with at least 50 prograde satellites thought to have formed from circumplanetary material. Just like the Sun is not the only star surrounded by planets, extrasolar gas giants are likely surrounded by satellite systems. Here, we report on ALMA observations of four <40 Myr old stars with directly-imaged companions: PZ Tel, AB Pic, 51 Eri, and $\kappa$ And. Continuum emission at 1.3 mm is undetected for any of the systems. Since these are directly-imaged companions, there is knowledge of their temperatures, masses and locations. These allow for upper limits on the amount of circumplanetary dust to be derived from detailed radiative transfer models. These protolunar disc models consider two disc sizes: 0.4 and 0.04 times the exoplanet's Hill radius. The former is representative of hydrodynamic simulations of circumplanetary discs while the latter a case with significant radial drift of solids. The more compact case is also motivated by the semi-major axis of Callisto, enclosing Jupiter's Galilean satellites. All upper limits fall below the expected amount of dust required to explain regular satellite systems ($10^{-4}$ times the mass of their central planet). Upper limits are compared with viscous evolution and debris disc models. Our analysis suggests that the non detections can be interpreted as evidence of dust growth beyond metre sizes to form moonetesimals in timescales <10 Myr. This sample increases by 50% the number of ALMA non-detections of young companions available in the literature.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports ALMA 1.3 mm non-detections toward four young (<40 Myr) directly-imaged exoplanets (PZ Tel b, AB Pic b, 51 Eri b, κ And b). Using radiative-transfer models of protolunar discs with two fixed radial extents (0.4 R_Hill and 0.04 R_Hill) and standard dust properties, the authors derive 3σ upper limits on circumplanetary dust mass that lie below 10^{-4} M_p for all targets. They compare these limits to viscous-evolution and debris-disc models and conclude that the non-detections constitute evidence for grain growth beyond metre sizes (moonetesimal formation) on timescales <10 Myr. The work increases the published sample of such ALMA non-detections by 50%.
Significance. If the upper limits remain robust under reasonable variations in disc geometry and opacity, the result supplies the first direct observational indication that the rapid satellite-formation channel inferred for the Solar System also operates around extrasolar gas giants. The addition of four new systems to the small existing sample is a concrete, incremental contribution to the empirical timeline of circumplanetary dust evolution.
major comments (3)
- [Abstract and §5] Abstract and §5: the central interpretive claim—that the non-detections imply dust growth beyond metre sizes within <10 Myr—requires that the radiative-transfer mass limits lie below the 10^{-4} M_p satellite-formation threshold for the adopted disc sizes. No quantitative exploration is presented of how the limits change for disc radii outside the 0.04–0.4 R_Hill range or for grain-size distributions that reduce mm-wave opacity relative to the assumed ISM-like properties; either change would map the same ALMA flux upper bound to a higher total dust mass and could remove the tension with the formation models.
- [§3] §3 (radiative-transfer setup): the two disc sizes are motivated by hydrodynamical simulations and by Callisto’s orbit, respectively, but the paper does not demonstrate that these choices bracket the plausible range of circumplanetary material distributions. A brief sensitivity test (e.g., a factor-of-two change in outer radius or a different surface-density power-law index) would directly address whether the reported limits are conservative or optimistic.
- [§4] §4 (results): the upper limits are stated to fall below 10^{-4} M_p, yet the precise numerical values, their dependence on the assumed planet parameters, and the conversion from flux to mass for each of the two disc sizes are not tabulated or plotted; without these data it is difficult to assess how close the limits sit to the threshold or how they would shift under modest model variations.
minor comments (2)
- [Abstract] The abstract states that the 0.04 R_Hill case is “also motivated by the semi-major axis of Callisto”; a one-sentence clarification of the numerical correspondence would help readers unfamiliar with the Galilean system.
- [Throughout] Figure captions and text use “Hill radius” without an explicit definition or reference to the standard expression R_H = a (M_p / 3 M_*)^{1/3}; adding the definition once would improve accessibility.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive report. We address each major comment below. Where the comments identify areas for improvement, we have revised the manuscript accordingly to strengthen the presentation and robustness of the results.
read point-by-point responses
-
Referee: [Abstract and §5] Abstract and §5: the central interpretive claim—that the non-detections imply dust growth beyond metre sizes within <10 Myr—requires that the radiative-transfer mass limits lie below the 10^{-4} M_p satellite-formation threshold for the adopted disc sizes. No quantitative exploration is presented of how the limits change for disc radii outside the 0.04–0.4 R_Hill range or for grain-size distributions that reduce mm-wave opacity relative to the assumed ISM-like properties; either change would map the same ALMA flux upper bound to a higher total dust mass and could remove the tension with the formation models.
Authors: The two adopted disc radii are motivated by specific hydrodynamical results and by the Callisto orbit as representative cases; we agree that a quantitative sensitivity study would better support the robustness of the <10^{-4} M_p conclusion. For grain-size distributions with lower mm-wave opacity, the same flux upper limit would indeed correspond to a higher dust-mass upper limit. We will add a dedicated paragraph in §5 that (i) states the adopted opacities explicitly, (ii) notes that reduced opacity would relax the mass limits, and (iii) discusses why the ISM-like assumption remains the appropriate baseline for testing the rapid-growth hypothesis. A full parameter sweep is beyond the scope of the present work but will be flagged as future work. revision: partial
-
Referee: [§3] §3 (radiative-transfer setup): the two disc sizes are motivated by hydrodynamical simulations and by Callisto’s orbit, respectively, but the paper does not demonstrate that these choices bracket the plausible range of circumplanetary material distributions. A brief sensitivity test (e.g., a factor-of-two change in outer radius or a different surface-density power-law index) would directly address whether the reported limits are conservative or optimistic.
Authors: We will perform and report a brief sensitivity test in a new subsection of §3. This will include (a) varying the outer radius by a factor of two around each fiducial value while keeping total mass fixed, and (b) changing the surface-density power-law index from the adopted value. The resulting changes in the derived 3σ dust-mass upper limits will be quantified and compared with the 10^{-4} M_p threshold. revision: yes
-
Referee: [§4] §4 (results): the upper limits are stated to fall below 10^{-4} M_p, yet the precise numerical values, their dependence on the assumed planet parameters, and the conversion from flux to mass for each of the two disc sizes are not tabulated or plotted; without these data it is difficult to assess how close the limits sit to the threshold or how they would shift under modest model variations.
Authors: We agree that the numerical values and their dependence on planet mass, radius, and distance should be presented explicitly. We will add a new table in §4 that lists, for each target and each disc size: the 3σ flux upper limit, the corresponding dust-mass upper limit, the assumed planet parameters, and the conversion factor from flux to mass. A supplementary figure showing the mass limits versus planet mass will also be included. revision: yes
Circularity Check
No circularity: upper limits computed from external planet parameters and standard RT models
full rationale
The derivation computes ALMA flux upper limits into dust-mass upper bounds via radiative-transfer models whose inputs (planet mass, temperature, location, adopted disc radii 0.4/0.04 R_Hill, and standard dust opacities) are taken from direct-imaging literature and hydrodynamic simulations, not fitted to the present non-detections. The comparison to the 10^{-4} M_p satellite-formation threshold and the subsequent interpretive suggestion of rapid dust growth are downstream inferences, not self-definitional or fitted-input predictions. No load-bearing self-citation, uniqueness theorem, or ansatz smuggling is present in the reported chain; the result remains falsifiable against independent mass thresholds and disc-size assumptions.
Axiom & Free-Parameter Ledger
free parameters (1)
- disc size scaling factors
axioms (1)
- domain assumption Expected dust mass required to form regular satellite systems is 10^{-4} times the central planet mass
Reference graph
Works this paper leans on
- [1]
-
[2]
Bonnefoy, M., Chauvin, G., Rojo, P., et al.\ 2010, , 512, A52
work page 2010
-
[3]
Bonnefoy, M., Currie, T., Marleau, G.-D., et al.\ 2014, , 562, A111
work page 2014
-
[4]
Bonnefoy, M., Zurlo, A., Baudino, J. L., et al.\ 2016, , 587, A58
work page 2016
-
[5]
Booth, M., Jord \'a n, A., Casassus, S., et al.\ 2016, , 460, L10
work page 2016
- [6]
- [7]
- [8]
- [9]
-
[10]
Chauvin, G., Lagrange, A.-M., Dumas, C., et al.\ 2004, , 425, L29
work page 2004
-
[11]
Chauvin, G., Lagrange, A.-M., Zuckerman, B., et al.\ 2005, , 438, L29
work page 2005
-
[12]
Chauvin, G., Lagrange, A.-M., Dumas, C., et al.\ 2005, , 438, L25
work page 2005
-
[13]
Christiaens V., et al., 2019, ApJ, 877, L33
work page 2019
-
[14]
Cilibrasi, M., Szul \'a gyi, J., Mayer, L., et al.\ 2018, , 480, 4355
work page 2018
-
[15]
Crida, A., & Charnoz, S.\ 2012, Science, 338, 1196
work page 2012
-
[16]
Dr a \.z kowska J., Szul \'a gyi J., 2018, ApJ, 866, 142
work page 2018
-
[17]
P., et al.\ 2012, Astrophysics Source Code Library, ascl:1202.015
Dullemond, C. P., et al.\ 2012, Astrophysics Source Code Library, ascl:1202.015
work page 2012
-
[18]
Gillon, M., Triaud, A. H. M. J., et al.\ 2017, , 542, 456
work page 2017
-
[19]
Hayashi, C.\ 1981, Progress of Theoretical Physics Supplement, 70, 35
work page 1981
-
[20]
Hinkley, S., Pueyo, L., Faherty, J. K., et al.\ 2013, , 779, 153
work page 2013
-
[21]
Submillimeter emission associated with candidate protoplanets
Isella A., Benisty M., Teague R., Bae J., Keppler M., Facchini S., P \'e rez L. M., 2019, arXiv e-prints, arXiv:1906.06308
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[22]
J., Quinn, S., et al.\ 2016, , 822, L3
Jones, J., White, R. J., Quinn, S., et al.\ 2016, , 822, L3
work page 2016
-
[23]
Kennedy G. M., Wyatt M. C., Su K. Y. L., Stansberry J. A., 2011, MNRAS, 417, 2281
work page 2011
- [24]
-
[25]
Lagrange, A.-M., Gratadour, D., Chauvin, G., et al.\ 2009, , 493, L21
work page 2009
-
[26]
Lagrange, A.-M., Boccaletti, A., Langlois, M., et al.\ 2018, arXiv:1809.08354
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [27]
- [28]
-
[29]
R., Barman, T., et al.\ 2015, Science, 350, 64
Macintosh, B., Graham, J. R., Barman, T., et al.\ 2015, Science, 350, 64
work page 2015
-
[30]
Maire, A.-L., Bonnefoy, M., Ginski, C., et al.\ 2016, , 587, A56
work page 2016
-
[31]
Marino S., et al., 2016, MNRAS, 460, 2933
work page 2016
- [32]
-
[33]
Marois, C., Macintosh, B., Barman, T., et al.\ 2008, Science, 322, 1348
work page 2008
-
[34]
Matr \`a , L., Wyatt, M. C., et al.\ 2019, arXiv e-prints , arXiv:1902.04081
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[35]
McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K.\ 2007, Astronomical Data Analysis Software and Systems XVI, 376, 127
work page 2007
-
[36]
Morbidelli, A., Szul \'a gyi, J., Crida, A., et al.\ 2014, , 232, 266
work page 2014
- [37]
- [38]
-
[39]
Ricci L., Rome H., Pinilla P., et al.\ 2017, ApJ, 846, 19
work page 2017
-
[40]
Ricci L., et al., 2017, AJ, 154, 24
work page 2017
- [41]
- [42]
-
[43]
Song, I., Schneider, G., Zuckerman, B., et al.\ 2006, , 652, 724
work page 2006
-
[44]
Su, K. Y. L., MacGregor, M. A., Booth, M., et al.\ 2017, , 154, 225
work page 2017
- [45]
- [46]
-
[47]
G., M \'e nard, F., Caceres, C., et al.\ 2017, , 154, 26
Wolff, S. G., M \'e nard, F., Caceres, C., et al.\ 2017, , 154, 26
work page 2017
- [48]
- [49]
- [50]
- [51]
- [52]
- [53]
- [54]
-
[55]
M., & Isella, A.\ 2018, , 479, 1850
Zhu, Z., Andrews, S. M., & Isella, A.\ 2018, , 479, 1850
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.