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

Recognition: 2 theorem links

· Lean Theorem

The ALMA Survey of Gas Evolution of PROtoplanetary Disks (AGE-PRO): Constraints on disk turbulence, fragmentation velocity, and inner pebble fluxes

Lilian Luo , Paola Pinilla , Camila Pulgar\'es , Laura M. P\'erez , Miguel Vioque , Nicol\'as T. Kurtovic , Anibal Sierra , Carolina Agurto-Gangas

show 9 more authors

Rossella Anania John Carpenter Lucas A. Cieza Dingshan Deng James Miley Ilaria Pascucci Giovanni P. Rosotti Beno\^it Tabone Ke Zhang

Authors on Pith no claims yet

Pith reviewed 2026-05-15 16:41 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords protoplanetary disksdust evolutionturbulencefragmentation velocitypebble fluxALMA observationspressure trapssubstructures

0 comments

The pith

Dust evolution models with low turbulence and low fragmentation velocity best match observations for nearly half of 30 protoplanetary disks, with inner pebble fluxes correlating more strongly with disk age than with substructures.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper models dust evolution across a sample of 30 protoplanetary disks observed by ALMA to test how turbulence and fragmentation affect the transport of solids into planet-forming regions. Four model setups are run using the DustPy code, each combining a different turbulent viscosity value with a different fragmentation velocity threshold, while adding pressure traps by adjusting the gas surface density profile from the data. Nearly half the disks are reproduced best by the lowest-turbulence, lowest-fragmentation case. Models of smooth disks consistently underpredict the observed dust mass, pointing to unresolved traps. Pebble delivery rates into the inner disk track disk age more closely than the presence of visible substructures.

Core claim

Using the DustPy code on the AGE-PRO sample, no single combination of turbulent viscosity and fragmentation velocity reproduces all 30 disks, yet the configuration with α = 10^{-4} and v_frag = 1 m s^{-1} provides the best match for nearly half. Pressure traps are included by perturbing the gas surface density according to the 1.3 mm continuum profiles. Synthetic images generated with RADMC-3D are compared directly to the data. Smooth disks without detected substructures are underpredicted in dust mass, while the inward pebble flux shows a clearer dependence on evolutionary stage than on radial substructure presence.

What carries the argument

DustPy dust evolution simulations that track grain growth, fragmentation, and radial drift for four pairs of α and v_frag values, with pressure traps added through direct perturbation of the gas surface density profile derived from continuum observations.

If this is right

The low-turbulence, low-fragmentation model supplies a workable description for a large fraction of disks across different ages.
Apparently smooth disks likely contain unresolved pressure traps that retain more dust than current models allow.
Time-dependent radial transport of pebbles dominates over static substructure effects in setting inner-disk solid delivery.
The resulting pebble fluxes provide a reference for interpreting upcoming multi-wavelength and JWST water-vapor data.

Where Pith is reading between the lines

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

Higher-resolution imaging could test whether the underpredicted dust masses in smooth disks arise from hidden traps.
Planet-formation timelines may shift systematically with disk age because of the age-linked change in pebble supply.
Disks with different morphologies may need coupled gas-dust hydrodynamic calculations for more accurate trap strengths.

Load-bearing premise

Perturbing the gas surface density profile from the continuum intensity is enough to capture the dust-trapping effect of substructures without running full hydrodynamic models of the gas.

What would settle it

High-resolution observations that confirm the absence of substructures in currently smooth disks while still showing higher dust masses than the smooth models predict, or direct measurements of inner-disk pebble fluxes that fail to increase with disk age.

read the original abstract

How substructures and disk properties affect dust evolution and the delivery of solids and volatiles into planet-forming regions remains an open question. We present results from tailored dust evolution modeling of the AGE-PRO ALMA large program, a sample of 30 protoplanetary disks spanning different evolutionary stages. Visibility fitting of the AGE-PRO ALMA data (at 1.3\,mm) reveals that approximately half of the disks exhibit radial substructures. Combined with stellar properties, disk inclinations, and gas mass estimates from CO isotopologues and N$_2$H$^+$, this well-characterized set of disks provides an ideal testbed to constrain dust evolution models across different ages and disk morphologies. Using the dust evolution code \texttt{DustPy}, we simulate dust evolution in each disk under four model configurations, varying two key free parameters: the turbulent viscosity ($\alpha = 10^{-4}, 10^{-3}$) and fragmentation velocity ($v_{\rm{frag}} = 1 \mathrm{m\,s^{-1}}, 10 \mathrm{m\,s^{-1}}$). Pressure traps are incorporated by perturbing the gas surface density based on the continuum intensity profiles, and synthetic observations generated with \texttt{RADMC-3D} are compared to these profiles. While no single model fits all disks, nearly half are best reproduced by the configuration with low turbulence and low fragmentation velocity ($\alpha = 10^{-4}, v_{\rm{frag}} = 1\,\mathrm{m\,s^{-1}}$). Models of smooth disks underpredict dust mass, possibly indicating unresolved substructures. Pebble fluxes into inner disk regions correlate more strongly with disk age than with the presence of substructures, highlighting time-dependent dust transport as a key factor in shaping inner disk composition. Our results also provide a comparative baseline for interpreting multiwavelength and JWST water vapor observations.

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.

Desk Editor's Note private letter to a colleague

The paper gives concrete model preferences for low turbulence and low fragmentation in half the AGE-PRO disks plus an age correlation for inner pebble flux, but the pressure trap modeling is a simplified approximation.

read the letter

The main thing to know is that this paper finds low turbulence and low fragmentation velocity best match nearly half the disks in their sample, while inner pebble fluxes correlate more with disk age than with substructures. They apply the DustPy code to the full AGE-PRO set of 30 disks. Each gets its own run using measured gas masses, inclinations, and substructures identified from ALMA 1.3 mm data. The four configurations test alpha of 10^{-4} versus 10^{-3} and v_frag of 1 versus 10 m/s. Substructures enter by perturbing the gas surface density profile to match the continuum intensity, then RADMC-3D produces synthetic observations for comparison. This gives a useful set of constraints because the sample spans evolutionary stages with solid ancillary data. The preference for the low-alpha, low-v_frag case in many disks and the stronger age dependence for pebble delivery are new quantifications for this dataset. Smooth disk models falling short on dust mass also flags that unresolved structure may be common. The soft spot is the pressure trap treatment. Directly adjusting gas density from observed intensity avoids full hydrodynamic simulations, so the traps do not evolve with dust feedback or planet-induced flows. If this overestimates trapping, the low-turbulence models will look stronger than they should, and the age-versus-substructure correlation could shift. The abstract mentions no single model fits all disks, which is honest but leaves the overall robustness dependent on how the 'best' fits are scored. This work is for researchers modeling dust evolution and solid delivery in protoplanetary disks. It supplies reference points that planet formation calculations can use. The data handling looks careful enough that it should go to peer review rather than get desk rejected.

Referee Report

2 major / 1 minor

Summary. The paper reports tailored DustPy dust evolution simulations for the 30-disk AGE-PRO ALMA sample. Four discrete model grids are run by varying turbulent viscosity α (10^{-4} or 10^{-3}) and fragmentation velocity v_frag (1 or 10 m s^{-1}); substructures are added by perturbing the gas surface density profile directly from the observed 1.3 mm continuum intensity. Synthetic RADMC-3D images are compared to the data. No single configuration reproduces the full sample, but nearly half the disks are reported as best matched by the low-α/low-v_frag run. Smooth-disk models under-predict dust mass, and inner pebble fluxes are found to correlate more strongly with disk age than with the presence of substructures.

Significance. If the modeling framework is shown to be robust, the work supplies the first systematic, observationally anchored constraints on turbulence and fragmentation across a statistically useful range of disk ages and morphologies. The finding that time-dependent pebble delivery dominates over substructure trapping for inner-disk fluxes would be a useful baseline for interpreting JWST water-vapor and multi-wavelength data. The use of a well-characterized sample with CO-derived gas masses is a clear strength.

major comments (2)

[modeling description] Modeling section: pressure traps are implemented by directly perturbing Σ_gas from the observed continuum intensity profile. This bypasses self-consistent hydrodynamics, dust back-reaction, and planet-induced flows; if the imposed perturbation overestimates trap depth or lifetime, the low-turbulence (α=10^{-4}) models will be artificially favored and the reported age-versus-substructure correlation for pebble flux will be biased.
[results] Results section: the claim that 'nearly half' the disks are best reproduced by α=10^{-4}, v_frag=1 m s^{-1} is presented without quantitative fit metrics (e.g., reduced χ², residual maps, or formal model-selection statistics). It is therefore unclear how robust the ranking is to noise, calibration uncertainties, or the choice of radial weighting.

minor comments (1)

[abstract] The abstract states that smooth-disk runs 'underpredict dust mass' but does not quantify the typical factor or its dependence on stellar mass or age; a short table or histogram would clarify this statement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

Thank you for the constructive feedback on our manuscript. We have carefully considered each major comment and provide point-by-point responses below. Where appropriate, we have made revisions to strengthen the paper.

read point-by-point responses

Referee: Modeling section: pressure traps are implemented by directly perturbing Σ_gas from the observed continuum intensity profile. This bypasses self-consistent hydrodynamics, dust back-reaction, and planet-induced flows; if the imposed perturbation overestimates trap depth or lifetime, the low-turbulence (α=10^{-4}) models will be artificially favored and the reported age-versus-substructure correlation for pebble flux will be biased.

Authors: We agree that our modeling approach involves simplifications by directly perturbing the gas surface density profile to match observed substructures, which does not include full hydrodynamical simulations or dust back-reaction effects. This choice was made to anchor the models closely to the observations without introducing additional assumptions about planet masses or migration. We acknowledge the potential for overestimating trap depths, which could bias towards lower α values. In the revised manuscript, we have expanded the discussion in Section 3 to explicitly address these limitations and their possible impact on the results. Additionally, we have performed a sensitivity test on trap depths and included it in the appendix to assess robustness of the age-substructure correlation. revision: partial
Referee: Results section: the claim that 'nearly half' the disks are best reproduced by α=10^{-4}, v_frag=1 m s^{-1} is presented without quantitative fit metrics (e.g., reduced χ², residual maps, or formal model-selection statistics). It is therefore unclear how robust the ranking is to noise, calibration uncertainties, or the choice of radial weighting.

Authors: The classification of 'best reproduced' was primarily based on qualitative assessment of how well the synthetic radial intensity profiles matched the observed ones across the sample. We recognize that this lacks the rigor of quantitative metrics. To address this concern, we have added quantitative comparisons using reduced chi-squared statistics for the radial profiles in the results section of the revised manuscript. We have also included example residual maps and discussed the impact of noise and calibration uncertainties on the model selection. This strengthens the robustness of our claim that nearly half the disks favor the low-α, low-v_frag configuration. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper runs forward DustPy simulations for four discrete (α, v_frag) combinations, perturbs Σ_gas from observed 1.3 mm continuum to impose traps, generates RADMC-3D synthetics, and selects the best-matching configuration per disk by direct comparison. The reported 'constraints' (nearly half best-fit by α=10^{-4}, v_frag=1 m s^{-1}) and the age-vs-substructure correlation for pebble fluxes are outputs of this model-selection process applied to independent dust-evolution physics. No step reduces by construction to its inputs: the match is not guaranteed, the perturbation is an explicit modeling choice rather than a self-definition, and no load-bearing self-citation or uniqueness theorem is invoked. This is standard parameter-constraint modeling against external data, not a circular derivation.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The results rest on two varied parameters chosen to match observations and on standard assumptions within the DustPy dust-evolution framework plus an ad-hoc treatment of pressure traps.

free parameters (2)

turbulent viscosity alpha = 10^{-4} or 10^{-3}
Varied between 10^{-4} and 10^{-3} as one of the two key free parameters in the four model configurations
fragmentation velocity v_frag = 1 or 10 m/s
Varied between 1 and 10 m/s as one of the two key free parameters in the four model configurations

axioms (2)

domain assumption Dust coagulation, fragmentation, and radial drift physics as implemented in the DustPy code
Standard assumptions of the publicly available dust-evolution code used for all simulations
ad hoc to paper Pressure traps can be approximated by perturbing gas surface density profiles derived from continuum intensity
Method used to incorporate observed substructures without running full hydrodynamic simulations

pith-pipeline@v0.9.0 · 5735 in / 1531 out tokens · 46337 ms · 2026-05-15T16:41:22.009496+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Using the dust evolution code DustPy, we simulate dust evolution in each disk under four model configurations, varying two key free parameters: the turbulent viscosity (α=10^{-4},10^{-3}) and fragmentation velocity (v_frag=1 m s^{-1},10 m s^{-1}). Pressure traps are incorporated by perturbing the gas surface density based on the continuum intensity profiles
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

nearly half are best reproduced by the configuration with low turbulence and low fragmentation velocity (α=10^{-4}, v_frag=1 m s^{-1})

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.