Laura Martínez Pérez

Atacama Large Millimeter/submillimeter (ALMA) high-resolution surveys have so far focused mainly on large, bright protoplanetary disks, where substructures are routinely detected, leaving compact disks poorly characterized. Yet, the survival of millimeter emission in Class II sources implies that radial drift must have been halted by pressure traps, previously undetected due to limited angular resolution. We present ALMA Band 7 (0.9 mm) continuum observations of 11 of the smallest ( R ₆₈ < 10 au) Class II disks in Lupus, reaching 0 . ″ 02 (2 au) resolution, the highest currently available with ALMA at this wavelength. We detect substructures in four disks, including a clear inner cavity of radius 1.6 au, and hints of substructure in three more. Two disks also host faint outer rings beyond their bright inner core. Combining these results with demographic information on Class II disks observed with ALMA, we examine how the substructure detection rate varies with disk size, restricting to disks whose diameter is resolved by at least three resolution elements. We find that nearly all disks (≳95%) larger than 30 au show substructures, compared to a still substantial ∼50% between 3 and 30 au, with the lower rate likely affected by optical depth. The planet masses required to open gaps within 5 au fall in the super-Earth to mini-Neptune regime. These results demonstrate that substructures are ubiquitous across the disk population, and that planet formation conditions can be probed down to a few astronomical units from the host star.

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₂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 DustPy, we simulate dust evolution in each disk under four model configurations, varying two key free parameters: the turbulent viscosity (α = 10⁻⁴, 10⁻³) and fragmentation velocity (v_frag = 1m s⁻¹, 10m s⁻¹). Pressure traps are incorporated by perturbing the gas surface density based on the continuum intensity profiles, and synthetic observations generated with 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 (α = 10⁻⁴, v_frag = 1 m s⁻¹). 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.

Despite observational progress in planet formation, the stage in which planetesimals grow into planets remains poorly understood. During this phase, protoplanets may develop gaseous envelopes that are warmer than the surrounding disk gas, potentially providing observable signatures through molecules otherwise depleted in cold regions. In this Letter, we report the detection of the silicon sulfide isotopologues ²⁸ SiS J = 16−15 and ³⁰ SiS J = 18−17 in the protoplanetary disk around PDS 66 (MP Mus) at a significance of ∼5 σ −6 σ , using the Atacama Large Millimeter/submillimeter Array. These constitute the second and first detections of ²⁸ SiS and ³⁰ SiS in a protoplanetary disk, respectively. The emission appears as a compact source at r = 60 au in the southwestern region of the disk, unresolved with a ∼ 0 . ″ 5 beam, and shows a velocity consistent with Keplerian rotation, suggesting a protoplanetary origin. By modeling the line fluxes, we constrain the emitting radius to ∼0.5−4 au and estimate a SiS mass of 10 ²² –10 ²³ g, corresponding to at least ∼10% of the silicon contained in local dust grains. Because complete sublimation of a substantial fraction of dust grains by local processes is difficult to achieve, this result instead implies an accumulation of silicon from a larger region. We propose that a circumplanetary envelope surrounding a low-mass protoplanet, where pebble accretion and subsequent sublimation of grains may enhance gaseous silicon abundance with respect to observable dust grains around it, can account for the observed characteristics.

Context . High-contrast imaging relies on advanced coronagraphs and adaptive optics (AO) to attenuate the starlight. However, residual aberrations, especially non-common path aberrations between the AO channel and the coronagraph channel, limit the instrument performance. While post-processing techniques such as spectral or angular differential imaging (ADI) can partially address those issues, they suffer from self-subtraction and inefficiencies at small angular separations or when observations are conducted far from transit. Aims . We previously demonstrated the on-sky performance of coherent differential imaging (CDI), which offers a promising alternative. It allows for isolating coherent starlight residuals through speckle modulation, which can then be subtracted from the raw images during post-processing. This work aims to validate a CDI method on real science targets using VLT/SPHERE, demonstrating its effectiveness in imaging almost face-on circumstellar disks, which are typically challenging to retrieve with ADI. Methods . We temporally modulated the speckle field in VLT/SPHERE images, applying small phase offsets on the AO deformable mirror while observing stars surrounded by circumstellar material: HR 4796A, CPD-36 6759, HD 169142, and HD 163296. We hence separated the astrophysical scene from the stellar speckle field, whose lights are mutually incoherent. Results . Combining a dozen of data frames and reference coronagraph point spread functions through a Karhunen–Loève image projection framework, we recover the circumstellar disks without the artifacts that are usually introduced by common post-processing algorithms (e.g., self-subtraction). Conclusions . The CDI method therefore represents a promising strategy for calibrating the effect of static and quasi-static aberrations in future direct imaging surveys. Indeed, it is efficient, does not require frequent telescope slewing, and does not introduce image artifacts to first order.

2MASS J16120668–3010270 (hereafter 2MJ1612) is a young M0 star that hosts a protoplanetary disk in the Upper Scorpius star-forming region. Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations of 2MJ1612 show a mildly inclined disk (i = 37°) with a large dust-depleted gap (R _cav ≈ 0 . ″ 4 or 53 au). We present high-contrast Hα observations from MagAO-X on the 6.5 m Magellan telescope and new high-resolution submillimeter dust continuum observations with ALMA of 2MJ1612. On both 2025 April 13 and 16, we recovered a point source with Hα excess with a signal-to-noise ratio ≳5 within the disk gap in our MagAO-X angular and spectral differential images at a separation of 141.96 ± 2.10 mas (23.45 ± 0.29 au deprojected) from the star and a position angle ​​​​​of 159 . ° 00 ± 0 . ° 55. Furthermore, this Hα source is within close proximity to a K-band point source in the SPHERE/IRDIS observation taken on 2023 July 21. The astrometric offset between the K band and Hα source can be explained by orbital motion of a bound companion. Thus, our observations can be best explained by the discovery of an accreting protoplanet, 2MJ1612 b, with an estimated mass of 4 M _Jup and a Hα line flux ranging from (29.7 ± 7.5) × 10⁻¹⁶ erg s cm² to (8.2 ± 3.4) × 10⁻¹⁶ erg s cm². 2MJ1612 b is likely the third example of an accreting Hα protoplanet responsible for carving the gap in its host disk, joining PDS 70 b and c. Further study is necessary to confirm and characterize this protoplanet candidate and to identify any additional protoplanets that may also play a role in shaping the gap.

We perform visibility fitting to the dust continuum Band 6 1.3 mm data of the 30 protoplanetary disks in the Atacama Large Millimeter/submillimeter Array Survey of Gas Evolution of PROtoplanetary Disks (AGE-PRO) Large Program. We obtain disk geometries, dust-disk radii, and azimuthally symmetric radial profiles of the intensity of the dust continuum emission. We examine the presence of continuum substructures in the AGE-PRO sample by using these radial profiles and their residuals. We detect substructures in 15 out of 30 disks. We report five disks with large (>15 au) inner dust cavities. The Ophiuchus Class I disks show dust-disk substructures in ∼80% of the resolved sources. This evidences the early formation of substructures in protoplanetary disks. A spiral is identified in IRS 63, hinting to gravitational instability in this massive disk. We compare our dust-disk brightness radial profiles with gas-disk brightness radial profiles and discuss colocal substructures in both tracers. In addition, we discuss the evolution of dust-disk radii and substructures across Ophiuchus, Lupus, and Upper Scorpius. We find that disks in Lupus and Upper Scorpius with large inner dust cavities have typical gas-disk masses, suggesting an abundance of dust cavities in these regions. The prevalence of pressure dust traps at later ages is supported by a potential trend with time with more disks with large inner dust cavities (or transition disks) in Upper Scorpius and the absence of evolution of dust-disk sizes with time in the AGE-PRO sample. We propose this is caused by an evolutionary sequence with a high fraction of protoplanetary disks with inner protoplanets carving dust cavities.

Protoplanetary disk evolution can be deeply influenced by the UV radiation emitted by neighboring massive stars (mainly of spectral types O and B). We show that the process of external photoevaporation, which causes an outside-in depletion of disk material due to environmental UV radiation, can lead to a significant decrease in disk size, disk mass, and lifetime even at moderate irradiation levels (1–10 G₀). In this work, we investigate the role of external photoevaporation in shaping the masses and sizes of the 10 AGE-PRO disks in the Upper Scorpius (Upper Sco) region, which we estimate to be subject to far-ultraviolet (FUV) fluxes ranging between ∼2 and ∼12 G₀, on average. We compare the disk masses and sizes resulting from 1D numerical viscous evolution simulations, in which the effect of external photoevaporation is included, to the values retrieved from the AGE-PRO observations. While the pure viscous framework fails in adequately explaining the observed disk properties in Upper Sco, with the inclusion of external photoevaporation, we can successfully reproduce gas disk sizes for seven out of 10 sources within a factor <2, when the initial disk mass is 1%–10% of the stellar mass. We emphasize the importance of accounting for the environmental irradiation when comparing star-forming regions of different ages, even when moderate FUV irradiation fields are experienced, as in the case of Upper Sco.

The architecture of planetary systems depends on the evolution of the disks in which they form. In this work, we develop a population synthesis approach to interpret the Atacama Large Millimeter/submillimeter Array survey of Gas Evolution of PROtoplanetary Disks (AGE-PRO) measurements of disk gas mass and size considering two scenarios: turbulence-driven evolution with photoevaporative winds and MHD wind-driven evolution. A systematic method is proposed to constrain the distribution of disk parameters from the disk fractions, accretion rates, disk gas masses, and CO gas sizes. We find that turbulence-driven accretion with initially compact disks (R ₀ ≃ 5–20 au), low mass-loss rates, and relatively long viscous timescales (t _ν,0 ≃ 0.4–3 Myr or α _SS ≃ 2–4 × 10⁻⁴) can reproduce the disk fractions and gas sizes. However, the distribution of apparent disk lifetimes defined as the M D / M ̇ * ratio is severely overestimated by turbulence-driven models. On the other hand, MHD wind-driven accretion can reproduce the bulk properties of disk populations from Ophiuchus to Upper Scorpius assuming compact disks with an initial magnetization of about β ≃ 10⁵ (α _DW ≃ 0.5–1 × 10⁻³) and a magnetic field that declines with time. More studies are needed to confirm the low masses found by AGE-PRO, notably for compact disks that question turbulence-driven accretion. The constrained synthetic disk populations can now be used for realistic planet population models to interpret the properties of planetary systems on a statistical basis.

The potential for planet formation of a circumstellar disk depends on the dust and gas reservoirs, which evolve as a function of the disk age. The Atacama Large Millimeter/submillimeter Array AGE-PRO Large Program has measured several disk properties across three star-forming regions of different ages, and in this study, we compare the observational results to dust evolution simulations. Using DustPy for the dust evolution, and RADMC-3D for the radiative transfer, we ran a large grid of models spanning stellar masses of 0.25, 0.50, 0.75, and 1.0 M _⊙, with different initial conditions, including: disk sizes, disk gas masses, and dust-to-gas ratio, and viscosity. Our models are performed assuming smooth, weakly, or strongly substructured disks, aiming to investigate if any observational trend can favor or exclude the presence of dust traps. The observed gas masses in the disks of the AGE-PRO sample are not reproducible with our models, which only consider viscous evolution with constant α, suggesting that additional physical mechanisms play a role in the evolution of the gas mass of disks. When comparing the dust continuum emission fluxes and sizes at 1.3 mm, we find that most of the disks in the AGE-PRO sample are consistent with simulations that have either weak or strong dust traps. The evolution of spectral index in the AGE-PRO sample is also suggestive of an unresolved population of dust traps. Future observations at high angular resolution are still needed to test several hypotheses that result from comparing the observations to our simulations, including that more massive disks in gas mass have the potential to form dust traps at larger disk radii.

The evolution of the gas mass of planet-forming disks around young stars is crucial for our understanding of planet formation, yet it has proven hard to constrain observationally, due both to the difficulties of measuring gas masses and the lack of a homogeneous sample. Here we present a large grid of thermochemical models that we use to measure protoplanetary gas disk masses of AGE-PRO, the Atacama Large Millimeter/submillimeter Array survey of Gas Evolution in PROtoplanetary disks. AGE-PRO covers a sample of 30 disks around similar spectral type (M3-K6) stars with ages between 0.1 and 10 Myr. Our approach is to simultaneously fit observations of CO isotopologues and N₂H⁺, a complementary molecule produced when CO freezes out. We find that the median gas mass of the three regions decreases over time, from 7 . 0 − 2.6 + 4.4 × 1 0 − 3 M ⊙ in Ophiuchus (≲1 Myr) to 9 . 4 − 3.4 + 5.4 × 1 0 − 4 M ⊙ for Lupus (∼1–3 Myr) and 6 . 8 − 2.8 + 5.1 × 1 0 − 4 M ⊙ for Upper Sco (∼2–6 Myr), with ∼1 dex scatter in gas mass in each region. We note that the gas mass distributions for Lupus and Upper Sco look very similar, which could be due to survivorship bias for the latter. The median bulk CO abundance in the CO emitting layer is found to be a factor ∼10 lower than the interstellar medium value but does not significantly change between Lupus and Upper Sco. From Lupus to Upper Sco, the median gas-to-dust mass ratio increases by a factor ∼3 from ∼40 to ∼120, suggesting efficient inward pebble drift and/or the formation of planetesimals.