Recognition: no theorem link
Survival of Protoplanetary Disks in Upper Scorpius from Population Synthesis Models with External Photoevaporation
Pith reviewed 2026-05-15 19:21 UTC · model grok-4.3
The pith
Population synthesis models show viscous evolution plus mild external FUV photoevaporation produces 3-7 Myr disk lifetimes matching Upper Scorpius observations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Population synthesis models of viscous protoplanetary disks subject to mild external FUV radiation fields of 1-100 G0 produce gas disk lifetimes of 3-7 Myr that are consistent with observed dispersal timescales in Upper Scorpius. Low-mass stars are more susceptible to dispersal, with their fraction among survivors dropping from 76 percent to 51 percent by 10 Myr. The majority of long-lived disks have low viscosity alpha below 10 to the minus 3.5 and initial characteristic radius below 125 AU, while initial disk-to-star mass ratio has little effect. Surviving disks show a sharp initial decline in median mass and radius followed by a slight increase due to survivorship bias.
What carries the argument
Population synthesis models that couple viscous disk evolution, parameterized by the viscosity coefficient alpha, with external photoevaporation driven by uniform FUV fields of 1-100 G0.
If this is right
- Disk lifetimes fall in the 3-7 Myr range for viscosity parameters between 10^{-4} and 10^{-2}.
- Low-mass stars of 0.1 solar masses experience faster disk dispersal due to weaker gravitational binding.
- Long-lived disks are preferentially those with low viscosity below 10^{-3.5} and initial characteristic radii below 125 AU.
- Median gas disk mass and size drop sharply in the first 0.2 Myr then rise slightly due to survivorship bias.
- Surviving disk properties correlate with stellar mass and FUV field strength.
Where Pith is reading between the lines
- External photoevaporation may shape disk populations even in regions with only moderate radiation, suggesting it is not limited to intense environments.
- Disk fraction surveys in other clusters with independently measured FUV levels could directly test whether the same viscosity and radius preferences hold.
- The sharp early drop in median disk mass followed by an upturn implies that current observed samples are biased toward initially compact, low-viscosity systems.
- Varying the FUV field across a simulated cluster rather than holding it uniform would reveal how spatial differences affect the final survivor population.
Load-bearing premise
The external FUV radiation fields are assumed to be mild, uniform, and fixed in the 1-100 G0 range, while the explored initial disk conditions and viscosity range are taken as representative enough to achieve consistency with observations.
What would settle it
Direct measurements showing that the typical FUV field strength in Upper Scorpius lies well outside 1-100 G0, or an observed distribution of disk lifetimes that falls significantly outside the 3-7 Myr window, would falsify the models.
read the original abstract
We present population synthesis models of viscous protoplanetary disks subject to mild external far-ultraviolet (FUV) radiation fields ($F_{\rm UV}=1\text{-}100\,$G$_0$). Our simulations focus on gas disk evolution, exploring stellar masses drawn from an Initial Mass Function and a range of initial disk conditions. We quantify the fraction of surviving disks across $10\,\mathrm{Myr}$ of evolution, track the evolution of gas disk mass and size, and compare our results with observations of protoplanetary disks in the Upper Scorpius region, including the ten targets studied by the AGE-PRO ALMA Large Program. We find that models combining viscous evolution with external photoevaporation yield disk lifetimes of $3\text{-}7\,\mathrm{Myr}$, consistent with observed dispersal timescales, particularly for $10^{-4} \leq \alpha \leq 10^{-2}$. Low-mass stars ($0.1\,$M$_\odot$) are more susceptible to disk dispersal due to their weaker gravitational binding, with their fraction among all surviving disks dropping from $76\%$ at birth to $51\%$ by $10\,\mathrm{Myr}$. The majority of the long-lived disks are those with low viscosity $\alpha < 10^{-3.5}$ and initial characteristic radius $R_c < 125\,\mathrm{AU}$, while the initial disk-to-star mass ratio does not play an important role. The median gas disk mass and radius of the surviving disks exhibit a sharp decline in the first $0.2\,\mathrm{Myr}$ of evolution, followed by a slight increase that reflects survivorship bias. We also explore correlations between gas disk mass and size vs. stellar mass and FUV strength. Our findings highlight the critical role of external photoevaporation in shaping disk populations even at moderate levels of FUV radiation fields.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents population synthesis simulations of protoplanetary disks evolving under viscous accretion and mild external FUV photoevaporation (1-100 G0). Stellar masses are drawn from an IMF, and initial disk parameters (including characteristic radius Rc and disk-to-star mass ratio) are varied across a grid. The models track the fraction of surviving disks, gas mass, and radius evolution over 10 Myr, with direct comparison to Upper Scorpius observations including the ten AGE-PRO ALMA targets. Key results are that viscous + external-photoevaporation models produce disk lifetimes of 3-7 Myr consistent with observed dispersal timescales specifically for 10^{-4} ≤ α ≤ 10^{-2}, that low-mass stars are more susceptible to dispersal, and that long-lived disks preferentially have low α (<10^{-3.5}) and small initial Rc (<125 AU). Median properties of survivors show an initial sharp decline followed by a slight rise due to survivorship bias.
Significance. If the central results hold after addressing the parameter-anchoring concern, the work provides a clear demonstration that moderate external FUV fields can shape disk lifetimes and population statistics in regions like Upper Scorpius, without needing extreme radiation or ad-hoc dispersal mechanisms. The forward-modeling approach, explicit tracking of survivorship bias, and comparison to specific ALMA data constitute strengths that would make the findings useful for interpreting disk surveys and for input into planet-formation population models.
major comments (1)
- [Abstract] Abstract and results sections: The headline claim that models yield 3-7 Myr lifetimes consistent with Upper Scorpius data specifically for 10^{-4} ≤ α ≤ 10^{-2} is obtained by exploring this viscosity range together with initial disk conditions. No independent observational or theoretical anchor for the viscosity distribution in Upper Scorpius (e.g., from turbulence measurements, disk-size evolution independent of lifetime, or accretion-rate statistics) is provided. If the true α distribution lies outside this interval, the predicted survival fractions and median lifetimes would shift, undermining the claimed consistency.
minor comments (3)
- Methods: The exact functional forms and parameter distributions for initial disk mass, Rc, and stellar-mass sampling should be stated explicitly (including any truncation or IMF weighting) so that the runs are reproducible.
- Results: Survival fractions, median masses, and radii are reported without error bars or uncertainty ranges arising from the finite number of simulated systems or from variations in the assumed FUV field; these should be added to all quantitative statements and figures.
- Figure captions and text: The assumed uniformity of the FUV field (1-100 G0) and its spatial variation across the cluster should be justified or tested with a sensitivity run, as the abstract presents it as a fixed mild background.
Simulated Author's Rebuttal
We thank the referee for their constructive review and for highlighting the need to better contextualize the viscosity parameter. We address the major comment below and have revised the manuscript to clarify the basis for the explored α range while preserving the core results on disk lifetimes under mild external photoevaporation.
read point-by-point responses
-
Referee: [Abstract] Abstract and results sections: The headline claim that models yield 3-7 Myr lifetimes consistent with Upper Scorpius data specifically for 10^{-4} ≤ α ≤ 10^{-2} is obtained by exploring this viscosity range together with initial disk conditions. No independent observational or theoretical anchor for the viscosity distribution in Upper Scorpius (e.g., from turbulence measurements, disk-size evolution independent of lifetime, or accretion-rate statistics) is provided. If the true α distribution lies outside this interval, the predicted survival fractions and median lifetimes would shift, undermining the claimed consistency.
Authors: We agree that the manuscript explores a range of α values (10^{-4} to 10^{-2}) rather than deriving a specific distribution from independent Upper Scorpius observations. This range is motivated by standard values adopted in the viscous disk literature, informed by accretion-rate measurements and turbulence constraints from other regions (e.g., Taurus and Lupus). Our central result is that, within this commonly used interval, the combination of viscous evolution and mild external FUV photoevaporation reproduces the observed 3-7 Myr disk lifetimes in Upper Scorpius. We have revised the abstract and results sections to emphasize that consistency is shown for the explored α range, added a short discussion paragraph citing prior observational anchors for α, and noted the sensitivity of survival fractions to the α distribution. If future data constrain α outside this interval, the predicted lifetimes would indeed shift, but the models demonstrate how moderate external photoevaporation shapes populations for plausible viscosities. revision: partial
Circularity Check
No significant circularity; forward population synthesis compared to observations
full rationale
The paper runs forward population synthesis simulations of viscous disk evolution plus external photoevaporation across explored ranges of α (10^{-4} to 10^{-2}), initial disk masses, radii, and stellar masses drawn from an IMF. Resulting survival fractions, lifetimes (3-7 Myr), and median properties are then compared to Upper Scorpius observations (including AGE-PRO targets). No equation or claim reduces a derived quantity to a fitted input by construction, no self-citation chain bears the central result, and no ansatz or uniqueness theorem is smuggled in. The match for particular α values is an output of the physics-based runs, not a tautology. This is standard external validation of a model grid.
Axiom & Free-Parameter Ledger
free parameters (2)
- viscosity parameter alpha =
10^{-4} to 10^{-2}
- initial characteristic radius Rc =
<125 AU
axioms (2)
- domain assumption Disks evolve according to the standard viscous alpha prescription
- domain assumption External photoevaporation mass loss is set by FUV flux in the 1-100 G0 range
Forward citations
Cited by 1 Pith paper
-
The dispersal of compact protoplanetary discs
Compact protoplanetary discs disperse inside-out when photoevaporation is limited to their cut-off radius, unlike the outside-in dispersal seen in extended discs.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.