Disc lifetime distribution as a function of the mass of host star
Pith reviewed 2026-05-07 14:30 UTC · model grok-4.3
The pith
Protoplanetary disc lifetimes depend on host star mass, peaking later around low-mass stars.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By binning disc-fraction data from star clusters according to host-star mass and fitting Weibull distributions, the analysis finds that disc-lifetime distributions peak at 3.72 Myr for stars between 1 and 3 solar masses and at 7.20 Myr for stars between 0.01 and 0.2 solar masses, with both distributions broad and the low-mass distribution somewhat wider; the initial disc fraction is lower for higher-mass stars.
What carries the argument
Mass-binned Weibull distribution fitted to observed disc fractions versus cluster age, with shape and scale parameters allowed to vary between low-mass and high-mass stellar bins.
Load-bearing premise
Observations of disc fractions in clusters give an unbiased record of how long discs survive around stars of different masses.
What would settle it
Mass-resolved disc-fraction measurements in a single young cluster that lie outside the predicted Weibull curves for either the low-mass or high-mass bin.
Figures
read the original abstract
The lifetime of protoplanetary discs is a critical factor for planet formation. Although the mean disc lifetime provides an estimate of the typical period available for planet formation, it does not capture the substantial variability in individual disc lifetimes or their dependence on host star mass. This study addresses these limitations by deriving the disc lifetime distribution as a function of stellar mass. Our results reveal a pronounced mass-dependence. Performing a phenomenological fit using a Weibull distribution, we find the maxima of the distributions at $t_{max}^H =$3.72 Myr for high-mass stars ($\approx$ 1.00--3.00 $M_{\odot}$) and $t_{max}^L =$ 7.20 Myr for low-mass stars ($\approx$ 0.01--0.20 $M_{\odot}$) assuming an initial disc fraction of $f_{init} = 0.8$. All distributions are broad (typically 3.2 Myr $< \sigma <$ 4.7 Myr), with the distribution for low-mass stars being somewhat broader. Our analysis indicates that not all stars are initially surrounded by a disc (60% $< f_{init} <$ 90% at cluster zero age), and that the initial disc fraction is even lower ($f_{init} \approx$ 40%) for higher-mass stars. The potential mechanisms responsible for the observed spread and mass-dependence of disc lifetime distributions and initial disc fractions are discussed. Our primary aim is to demonstrate the methodology; more robust constraints will require improved data on mass-dependent disc fractions. Nevertheless, the derived mass-dependent disc lifetime distributions can already serve as a valuable input or a benchmark for planet-formation synthesis models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives mass-dependent distributions of protoplanetary disc lifetimes by performing phenomenological Weibull fits to observed disc fractions as a function of cluster age, separately for high-mass (≈1–3 M⊙) and low-mass (≈0.01–0.2 M⊙) stars. It reports distribution maxima of t_max^H = 3.72 Myr and t_max^L = 7.20 Myr (with f_init = 0.8), broad widths (3.2–4.7 Myr), lower initial disc fractions for higher-mass stars, and discusses possible mechanisms; the work is framed as a demonstrative methodology requiring better data.
Significance. If the reported mass dependence and broad distributions are robust, the results supply a useful, observationally motivated input for planet-formation synthesis models that currently rely on mean lifetimes. The explicit separation into mass bins and the allowance for f_init < 1 address a clear limitation of single-parameter lifetime estimates.
major comments (3)
- [Results section (phenomenological fit description)] The numerical values t_max^H = 3.72 Myr, t_max^L = 7.20 Myr and the quoted σ range are direct outputs of the Weibull fits, yet the manuscript provides no description of the underlying disc-fraction data points, their binning in age or mass, the fitting procedure (likelihood, priors, treatment of non-detections), or error estimation. Without these details the claimed mass dependence cannot be independently verified.
- [Discussion (mechanisms and assumptions)] The central claim of pronounced mass dependence assumes that the binned disc-fraction-versus-age data accurately reflect intrinsic lifetime distributions. Potential systematics—luminosity-dependent IR-excess completeness for high-mass stars and possible differences in cluster age calibrations—are not quantified or corrected for, so any offset in the high-mass bin would be absorbed into the reported t_max difference.
- [Abstract and Results (Weibull fit)] Because the Weibull parameters (scale, shape, f_init) are fitted directly to the same disc-fraction data that define the lifetime distribution, the reported maxima and widths are largely a reparameterization of the input observations rather than an independent derivation; this circularity limits the strength of the mass-dependence conclusion.
minor comments (3)
- [Abstract] The abstract states both f_init = 0.8 and the broader range 60 % < f_init < 90 %; the text should clarify which value applies to which mass bin and how the range was obtained.
- [Methods/Results] The symbol σ is used for the distribution width; its precise relation to the Weibull shape and scale parameters should be stated explicitly, preferably with the functional form of the distribution.
- [Throughout] A table or figure listing the input disc-fraction data points (age, fraction, mass bin, reference) would greatly improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments, which have helped improve the clarity and rigor of the manuscript. We have revised the paper to address the lack of methodological details and to expand the discussion of observational systematics. We maintain that the phenomenological Weibull modeling provides a valid way to quantify mass-dependent lifetime distributions from the available data, while acknowledging the demonstrative nature of the study.
read point-by-point responses
-
Referee: The numerical values t_max^H = 3.72 Myr, t_max^L = 7.20 Myr and the quoted σ range are direct outputs of the Weibull fits, yet the manuscript provides no description of the underlying disc-fraction data points, their binning in age or mass, the fitting procedure (likelihood, priors, treatment of non-detections), or error estimation. Without these details the claimed mass dependence cannot be independently verified.
Authors: We agree that these details were insufficient. The revised manuscript includes a new subsection in Results that describes the compilation of disc-fraction data from the literature for the two mass bins, the adopted age binning, the maximum-likelihood fitting of the Weibull model (with the explicit likelihood function accounting for binomial errors and non-detections treated as censored data), and bootstrap resampling for parameter uncertainties. The analysis code has also been made available as supplementary material to enable independent verification. revision: yes
-
Referee: The central claim of pronounced mass dependence assumes that the binned disc-fraction-versus-age data accurately reflect intrinsic lifetime distributions. Potential systematics—luminosity-dependent IR-excess completeness for high-mass stars and possible differences in cluster age calibrations—are not quantified or corrected for, so any offset in the high-mass bin would be absorbed into the reported t_max difference.
Authors: This is a valid concern. We have added a dedicated paragraph in the Discussion that addresses these systematics, noting that luminosity-dependent completeness for high-mass stars could bias disc fractions low at older ages and that cluster age uncertainties could shift the data. We explicitly state that no quantitative corrections were applied due to the absence of detailed completeness functions in the source studies, and that these effects may contribute to the reported t_max difference. We retain the mass-dependence conclusion as qualitative and consistent with other evidence (e.g., photoevaporation), while highlighting it as a limitation requiring future work. revision: partial
-
Referee: Because the Weibull parameters (scale, shape, f_init) are fitted directly to the same disc-fraction data that define the lifetime distribution, the reported maxima and widths are largely a reparameterization of the input observations rather than an independent derivation; this circularity limits the strength of the mass-dependence conclusion.
Authors: We respectfully disagree that this is problematic circularity. The binned disc fractions versus age constitute the empirical survival function; fitting a flexible Weibull form is a standard parametric approach in survival analysis to recover the underlying lifetime density, from which the mode t_max and width σ are obtained. The mass dependence arises from performing separate fits to the high-mass and low-mass subsets, yielding statistically different parameters. We have clarified this motivation in the revised Results and Discussion sections, emphasizing the phenomenological character of the model and its utility for comparison with theory. revision: no
Circularity Check
No significant circularity in the derivation of mass-dependent disc lifetime distributions
full rationale
The paper collects observational disc-fraction data from clusters of different ages, bins the data by stellar mass, and fits a Weibull distribution (with free parameters including f_init) to model the observed decline. The reported t_max^H = 3.72 Myr, t_max^L = 7.20 Myr and sigma ranges are the directly optimized parameters of this fit under the explicit phenomenological assumption of a Weibull form. Because the input data are external observations and the functional form is stated as an ansatz rather than derived, the output is an inference from independent inputs rather than a re-expression by construction. No self-citations, uniqueness theorems, or renamings of known results are invoked as load-bearing steps. The paper explicitly frames the work as a methodology demonstration requiring better data, confirming the chain remains non-circular.
Axiom & Free-Parameter Ledger
free parameters (4)
- t_max^H =
3.72 Myr
- t_max^L =
7.20 Myr
- sigma =
3.2-4.7 Myr
- f_init =
0.8 (range 0.6-0.9; ~0.4 for high-mass)
axioms (2)
- domain assumption The fraction of stars retaining discs declines with cluster age according to a Weibull survival function whose parameters depend only on stellar mass.
- domain assumption Observed disc fractions in young clusters directly trace the underlying lifetime distribution without major contamination from dynamical ejection or observational bias.
Reference graph
Works this paper leans on
-
[1]
Akeson, R. L., Jensen, E. L. N., Carpenter, J., et al. 2019, ApJ, 872, 158
work page 2019
-
[2]
Andrews, S. M. 2020, ARA&A, 58, 483
work page 2020
-
[3]
M., Kamp, I., Henning, T., et al
Arabhavi, A. M., Kamp, I., Henning, T., et al. 2024, Science, 384, 1086
work page 2024
- [4]
-
[5]
Baraffe, I., Homeier, D., Allard, F ., & Chabrier, G. 2015, A&A, 577, A42
work page 2015
-
[6]
Barenfeld, S. A., Carpenter, J. M., Sargent, A. I., et al. 2019, ApJ, 878, 45
work page 2019
- [7]
-
[8]
Blandford, R. D. & P ayne, D. G. 1982, MNRAS, 199, 883 Briceño, C., Calvet, N., Hernández, J., et al. 2019, AJ, 157, 85
work page 1982
-
[9]
D., Kamp, I., Meeus, G., Oudmaijer, R
Brittain, S. D., Kamp, I., Meeus, G., Oudmaijer, R. D., & Waters, L. B. F . M. 2023, SSRv, 219, 7
work page 2023
-
[10]
Coleman, G. A. L. & Haworth, T. J. 2022, MNRAS, 514, 2315–2332
work page 2022
-
[11]
Colmenares, M. J., Bergin, E. A., Salyk, C., et al. 2024, ApJ, 977, 173 Dickson-V andervelde, D. A., Wilson, E. C., & Kastner, J. H. 2021, AJ, 161, 87 Duchêne, G. & Kraus, A. 2013, ARA&A, 51, 269
work page 2024
- [12]
-
[13]
Ercolano, B. & P ascucci, I. 2017, Royal Society Open Science, 4, 170114
work page 2017
-
[14]
2014, in Protostars and Planets VI, ed
Espaillat, C., Muzerolle, J., Najita, J., et al. 2014, in Protostars and Planets VI, ed. H. Beuther, R. S. Klessen, C. P . Dullemond, & T. Henning, 497–520
work page 2014
-
[15]
E., Henning, T., Jayawardhana, R., & Oliveira, J
Fedele, D., van den Ancker, M. E., Henning, T., Jayawardhana, R., & Oliveira, J. M. 2010, A&A, 510, A72
work page 2010
- [16]
-
[17]
Galli, P . A. B., Bouy, H., Olivares, J., et al. 2021b, A&A, 646, A46 Gárate, M., Birnstiel, T., Pinilla, P ., et al. 2023, A&A, 679, A15
work page 2023
-
[18]
Gomez, J. L., Guilera, O. M., Miller Bertolami, M. M., Castro-Martínez, E., & Ronco, M. P . 2025, arXiv e-prints, arXiv:2510.06978
- [19]
-
[20]
Harris, R. J., Andrews, S. M., Wilner, D. J., & Kraus, A. L. 2012, ApJ, 751, 115 12 Hernández, J., Hartmann, L., Calvet, N., et al. 2008, ApJ, 686, 1195
work page 2012
-
[21]
Kaeufer, T., Waters, R., Gasman, D., et al. 2026, MNRAS, 545, staf2056
work page 2026
-
[22]
M., Matrà, L., F acchini, S., et al
Kennedy, G. M., Matrà, L., F acchini, S., et al. 2019, Nat. Astron., 3, 230
work page 2019
-
[23]
Kenyon, S. J. & Hartmann, L. 1995, ApJS, 101, 117
work page 1995
-
[24]
Komaki, A., Fukuhara, S., Suzuki, T. K., & Y oshida, N. 2023, arXiv e-prints, arXiv:2304.13316
-
[25]
Kral, Q., Matrà, L., Kennedy, G. M., Marino, S., & Wyatt, M. C. 2020, MNRAS, 497, 2811
work page 2020
-
[26]
Kraus, A. L., Ireland, M. J., Hillenbrand, L. A., & Martinache, F . 2012, ApJ, 745, 19
work page 2012
-
[27]
Kraus, A. L., Ireland, M. J., Huber, D., Mann, A. W., & Dupuy, T. J. 2016, AJ, 152, 8
work page 2016
- [28]
- [29]
- [30]
-
[31]
Lada, C. J. & Lada, E. A. 2003, ARA&A, 41, 57
work page 2003
- [32]
-
[33]
Luhman, K. L. 2021, AJ, 163, 24
work page 2021
-
[34]
Luhman, K. L. 2022, AJ, 163, 25
work page 2022
-
[35]
Luhman, K. L. & Esplin, T. L. 2020, AJ, 160, 44
work page 2020
- [36]
-
[37]
Manara, C. F ., Ansdell, M., Rosotti, G. P ., et al. 2023, in Astronomical Society of the P acific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y. Aikawa, T. Muto, K. T omida, & M. T amura, 539
work page 2023
-
[38]
F ., Prusti, T., Comeron, F ., et al
Manara, C. F ., Prusti, T., Comeron, F ., et al. 2018, A&A, 615, L1 Manzo-Martínez, E., Calvet, N., Hernández, J., et al. 2020, ApJ, 893, 56
work page 2018
-
[39]
Michel, A., van der Marel, N., & Matthews, B. C. 2021, ApJ, 921, 72
work page 2021
- [40]
-
[41]
Murphy, S. J., Mamajek, E. E., & Bell, C. P . M. 2018, MNRAS, 476, 3290
work page 2018
-
[42]
Offner, S. S. R., Moe, M., Kratter, K. M., et al. 2023, in Astronomical Society of the P acific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y. Aikawa, T. Muto, K. T omida, & M. T amura, 275 P ascucci, I., T esti, L., Herczeg, G. J., et al. 2016, ApJ, 831, 125 P ecaut, M. J. & Mamajek, E. E. 2016, MNRAS, 461, 794
work page 2023
- [43]
- [44]
- [45]
- [46]
-
[47]
Pfalzner, S., Steinhausen, M., & Menten, K. 2014, ApJL, 793, L34
work page 2014
-
[48]
2022, European Physical Journal Plus, 137, 1206
Pinilla, P . 2022, European Physical Journal Plus, 137, 1206
work page 2022
- [49]
-
[50]
Pinilla, P ., Sierra, A., Kurtovic, N. T., et al. 2025, MNRAS, 543, 2723 P olnitzky, F . A., Ratzenböck, S., Großschedl, J., & Alves, J. 2026, A&A, 707, A216 Ribas, Á., Bouy, H., & Merín, B. 2015, A&A, 576, A52 Ribas, Á., Merín, B., Bouy, H., & Maud, L. T. 2014, A&A, 561, A54
work page 2025
-
[51]
Richert, A. J. W., Getman, K. V ., Feigelson, E. D., et al. 2018, MNRAS, 477, 5191
work page 2018
-
[52]
Ronco, M. P ., Guilera, O. M., Cuadra, J., et al. 2021, ApJ, 916, 113
work page 2021
-
[53]
Bertolami, M. M. 2024, A&A, 682, A155
work page 2024
-
[54]
Sana, H., Le Bouquin, J.-B., Lacour, S., et al. 2014, ApJS, 215, 15
work page 2014
-
[55]
Guerra-Alvarado, O. M. 2024, A&A, 689, A236
work page 2024
- [56]
-
[57]
Sellek, A. D. & van Dishoeck, E. F . 2025, A&A, 701, A239
work page 2025
-
[58]
Silverberg, S. M., Kuchner, M. J., Wisniewski, J. P ., et al. 2016, ApJL, 830, L28 van der Marel, N. 2023, European Physical Journal Plus, 138, 225 van der Marel, N. & Mulders, G. D. 2021, AJ, 162, 28 van T erwisga, S. E., Hacar, A., van Dishoeck, E. F ., Oonk, R., & P ortegies Zwart, S. 2022, A&A, 661, A53
work page 2016
- [59]
- [60]
-
[61]
Wang, J., Fischer, D. A., Xie, J.-W., & Ciardi, D. R. 2014, ApJ, 791, 111
work page 2014
-
[62]
Wilhelm, M. J. C. & P ortegies Zwart, S. 2022, MNRAS, 509, 44
work page 2022
-
[63]
Xie, C., P ascucci, I., Long, F ., et al. 2023, ApJL, 959, L25 Stellar mass-dependent disc lifetime distributions13 APPENDIX T able A.1.Disc fractions of star clusters closer than 200 au. IdentificationdAgeN ⋆ fd Limit Median mass logρ c Source pc Myr [M ⊙] [M ⊙/pc3] Alessi 13 108 30 162 0.049 a 0.04M ⊙ a) UCL/LCC 150 15–20 3665 0.09±0.01 0.15 -0.85–(-1.0...
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.