Search

This work aims at providing fundamental general tools for the analysis of water spectra as observed in protoplanetary disks with JWST-MIRI. We analyze 25 high-quality spectra from the JDISC Survey reduced with asteroid calibrators as presented in Pontoppidan et al. 2024. First, we present a spectral atlas to illustrate the clustering of water transitions from different upper level energies ($E_u$) and identify single (un-blended) lines that provide the most reliable measurements. With the atlas, we demonstrate two important excitation effects: one related to the opacity saturation of ortho-para line pairs that overlap, and the other to the sub-thermal excitation of $v=1-1$ lines scattered across the $v=0-0$ rotational band. Second, from this larger line selection we define a list of fundamental lines spanning $E_u$ from 1500 to 6000 K to develop simple line-ratio diagrams as diagnostics of temperature components and column density. Third, we report the detection of disk-rotation Doppler broadening of molecular lines, which demonstrates the radial distribution of water emission at different $E_u$ and confirms from gas kinematics a radially-extended $\approx 170$--220 K reservoir, close to the ice sublimation front. We also report the detection of narrow blue-shifted absorption from an inner disk wind in ro-vibrational H2O and CO lines, which may be observed in disks at inclinations $> 50$ deg. We summarize these findings and tools into a general recipe that should be beneficial to community efforts to study water in planet-forming regions., Comment: Community input and feedback is very welcome. Referee report received and asking for minor revisions. Some analysis updates are in progress. The rest of the sample in Appendix H will be included later

The mid-infrared water vapor emission spectrum provides a novel way to characterize the delivery of icy pebbles towards the innermost ($<5$ au) regions of planet-forming disks. Recently, JWST MIRI-MRS showed that compact disks exhibit an excess of low-energy water vapor emission relative to extended multi-gapped disks, suggesting that icy pebble drift is more efficient in the former. We carry out detailed emission line modeling to retrieve the excitation conditions of rotational water vapor emission in a sample of four compact and three extended disks within the JDISC Survey. We present two-temperature H$_2$O slab model retrievals and, for the first time, constrain the spatial distribution of water vapor by fitting parametric radial temperature and column density profiles. Such models statistically outperform the two-temperature slab fits. We find a correlation between the observable hot water vapor mass and stellar mass accretion rate, as well as an anti-correlation between cold water vapor mass and sub-mm dust disk radius, confirming previously reported water line flux trends. We find that the mid-IR spectrum traces H$_2$O with temperatures down to 180-300 K, but the coldest 150-170 K gas remains undetected. Furthermore the H$_2$O temperature profiles are generally steeper and cooler than the expected `super-heated' dust temperature in passive irradiated disks. The column density profiles are used to estimate icy pebble mass fluxes, which suggest that compact and extended disks may produce markedly distinct inner-disk exoplanet populations if local feeding mechanisms dominate their assembly., Comment: Accepted for publication in ApJ

Detecting planet signatures in protoplanetary disks is fundamental to understanding how and where planets form. In this work, we report dust and gas observational hints of planet formation in the disk around 2MASS-J16120668-301027, as part of the ALMA Large Program "AGE-PRO: ALMA survey of Gas Evolution in Protoplanetary disks". The disk was imaged with the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 6 (1.3 mm) in dust continuum emission and four molecular lines: $^{12}$CO(J=2-1), $^{13}$CO(J=2-1), C$^{18}$O(J=2-1), and H$_2$CO(J=3$_{(3,0)}$-2$_{(2,0)}$). Resolved observations of the dust continuum emission (angular resolution of $\sim 150$ mas, 20 au) show a ring-like structure with a peak at $0.57 ^{\prime \prime}$ (75 au), a deep gap with a minimum at 0.24$^{\prime \prime}$ (31 au), an inner disk, a bridge connecting the inner disk and the outer ring, along with a spiral arm structure, and a tentative detection (to $3\sigma$) of a compact emission at the center of the disk gap, with an estimated dust mass of $\sim 2.7-12.9$ Lunar masses. We also detected a kinematic kink (not coincident with any dust substructure) through several $^{12}$CO channel maps (angular resolution $\sim$ 200 mas, 30 au), located at a radius of $\sim 0.875^{\prime \prime}$ (115.6 au). After modeling the $^{12}$CO velocity rotation around the protostar, we identified a tentative rotating-like structure at the kink location with a geometry similar to that of the disk. We discuss potential explanations for the dust and gas substructures observed in the disk, and their potential connection to signatures of planet formation., Comment: 24 pages, 15 figures

The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$\alpha$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs., Comment: 23 pages, 23 figures, accepted by A&A (abstract shortened)

Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport varying the initial conditions of the disk and substructure to find the best match to the observed distributions. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure, namely a perturbation strong enough to be able to trap particles, and that is formed early in the evolution of the disk, that is within 0.4Myr. Agreement is reached by relatively high initial disk masses ($10^{-2.3}M_{\star}\leqslant M_{disk}\leqslant10^{-0.5}M_{\star}$) and moderate levels of turbulence ($10^{-3.5}\leqslant\alpha\leqslant 10^{-2.5}$). Other disk parameters play a weaker role. Only opacities with high absorption efficiency can reproduce the observed spectral indices. Our results extend to the whole population that substructure is likely ubiquitous, so far assessed only in individual disks and implies that most "smooth" disks hide unresolved substructure.

The composition and distribution of the gas in a protoplanetary disk plays a key role in shaping the outcome of the planet formation process. Observationally, the recovery of information such as the emission height and brightness temperature from interferometric data is often limited by the imaging processes. To overcome the limitations of image-reconstruction when analyzing gas emission from interferometric observations, we have introduced a parametric model to fit the main observable properties of the gaseous disk component in the visibility plane. This approach is also known as parametric visibility modeling. We applied our parametric visibility modeling to the gas brightness distribution of the molecular line emission from 12CO J=3-2 and 13CO J=3-2 in the disk around MHO 6, a very-low-mass star in the Taurus star-forming Region. To improve the flux fidelity of our parametric models, we combined models with different pixel resolution before the computation of their visibilities, referred to as ``nesting images.'' When we apply our parametric visibility modeling to MHO 6, with independent fits to the emission from its CO isopotologues, the models return the same consistent results for the stellar mass, disk geometry, and central velocity. The surface height and brightness temperature distribution are also recovered. When compared to other disks, MHO 6 surface height is among the most elevated surfaces, consistent with the predictions for disks around very-low-mass stars. This work demonstrates the feasibility of running rapidly iterable parametric visibility models in moderate resolution and sensitivity interferometric observations. More importantly, this methodology opens the analysis of disk's gas morphology to observations where image-based techniques are unable to robustly operate, as in the case of the compact disk around MHO 6., Comment: Accepted for publication in A&A

Aims. To comprehend the efficiency of dust evolution within protoplanetary disks, it is crucial to conduct studies of these disks using high-resolution observations at multiple wavelengths with the Atacama Large Millimeter/submillimeter Array (ALMA). Methods. In this work, we present high-frequency ALMA observations of the HL Tau disk using its Band 9 centered at a wavelength of 0.45 mm. These observations achieve the highest angular resolution in a protoplanetary disk to date, 12 milliarcseconds (mas), allowing the study of the dust emission at scales of 2 au. We use these data to extend the previously published multi-wavelength analysis of the HL Tau disk. Results. Our new 0.45 mm data traces mainly optically thick emission, providing a tight constraint to the dust temperature profile. We derive maximum particle sizes of $\sim$1 cm from the inner disk to $\sim$60 au. Beyond this radius, we find particles between 300 $\mu$m and 1 mm. Moreover, an intriguing asymmetry is observed at 32 au in the northeast inner part of the HL Tau disk at 0.45 mm. We propose that this asymmetry is the outcome of a combination of factors including the optically thick nature of the emission, the orientation of the disk, and a relatively large dust scale height of the grains. To validate this, we conducted a series of radiative transfer models using the RADMC-3D software. If this scenario is correct, our measured dust mass within 32 au would suggest a dust scale height H/R> 0.08 for the inner disk. Finally, the unprecedented resolution allowed us to probe for the first time the dust emission down to a few au scales. We observed an increase in brightness temperature inside the estimated water snowline and speculate whether it could indicate the presence of a traffic jam effect in the inner disk. Abridge, Comment: Accepted for publication in A&A, 13 pages, 8 figures

The K7 T Tauri star PDS 70 remains the best laboratory for investigating the influence of giant planet formation on the structure of the parental disk. One of the most intriguing discoveries is the detection of a resolved inner disk from ALMA observations that extends up to the orbit of PDS 70b. It is challenging to explain this inner disk because most of the dust particles are expected to be trapped at the outer edge of the gap opened by PDS 70b and PDS 70c. By performing dust evolution models in combination with radiative transfer simulations that match the gas disk masses obtained from recent thermo-chemical models of PDS 70, we find that when the minimum grain size in the models is larger than 0.1$\mu$m, there is an efficient filtration of dust particles, and the inner disk is depleted during the first million year of dust evolution. To maintain an inner disk, the minimum grain size in the models therefore needs to be smaller than 0.1$\mu$m. Only when grains are that small are they diffused and dragged along with the gas throughout the gap opened by the planets. The small grains transported in the inner disk grow and drift into it, but the constant reservoir of dust particles that are trapped at the outer edge of the gap and that continuously fragment allows the inner disk to refill on million-year timescales. Our flux predictions at millimeter wavelength of these models agree with ALMA observations. These models predict a spectral index of 3.2 in the outer and 3.6 in the inner disk. Our simple analytical calculations show that the water emission in the inner disk that was recently observed with the James Webb Space Telescope may originate from these ice-coated small grains that flow through the gap, grow, and drift toward the innermost disk regions to reach the water snowline., Comment: Accepted for publication in A&A. Final version after language edition

High-angular resolution imaging by ALMA has revealed the near-universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M-dwarfs in the galaxy. Here we present high-resolution (~50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ~10 au for 2M0434, 2M0436 and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared SEDs to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet-disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ~10 au cavities under current observations and theories. Currently, disks around mid-to-late M stars in Taurus show a higher detection frequency of cavities than earlier type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass., Comment: 25 pages, 6 figures. Accepted for publication in ApJ

The dust in planet-forming disks evolve rapidly through growth and radial drift, and external photoevaporation also contributes to this evolution in massive star-forming regions. We test whether the presence of substructures can explain the survival of the dust component and observed millimeter continuum emission in protoplanetary disks located within massive star-forming regions. We also characterize the dust content removed by the photoevaporative winds. For this, we performed hydrodynamical simulations of protoplanetary disks subject to irradiation fields of $F_{UV} = 10^2$, $10^3$, and $10^4\, G_0$, with different dust trap locations. We used the FRIED grid to derive the mass loss rate for each irradiation field and disk properties, and then measure the evolution of the dust mass over time. For each simulation we estimate continuum emission at $\lambda = 1.3\, \textrm{mm}$ along with the radii encompassing $90\%$ of the continuum flux, and characterize the dust size distribution entrained in the photoevaporative winds, along with the resulting far-ultraviolet (FUV) cross section. Our simulations show that the presence of dust traps can extend the lifetime of the dust component of the disk to a few millionyears if the FUV irradiation is $F_{UV} \lesssim 10^3 G_0$, but only if the dust traps are located inside the photoevaporative truncation radius. The dust component of a disk quickly disperse if the FUV irradiation is strong ($10^4\, G_0$) or if the substructures are located outside the photoevaporation radius. We do find however, that the dust grains entrained with the photoevaporative winds may result in an absorption FUV cross section of $\sigma \approx 10^{-22}\, \textrm{cm}^2$ at early times of evolution ($<$0.1 Myr), which is enough to trigger a self-shielding effect that reduces the total mass loss rate, and slow down the disk dispersal in a negative feedback loop process., Comment: Accepted for publication in A&A

(Exo-)planets inherit their budget of chemical elements from a protoplanetary disk. The disk temperature determines the phase of each chemical species, which sets the composition of solids and gas available for planet formation. We investigate how gap structures, which are widely seen by recent disk observations, alter the thermal and chemical structure of a disk. Planet-disk interaction is a leading hypothesis of gap formation and so such changes could present a feedback that planets have on planet-forming material. Both the planet gap-opening process and the disk thermal structure are well studied individually, but how the gap-opening process affects disk thermal structure evolution remains an open question. We develop a new modelling method by iterating hydrodynamical and radiative transfer simulations to explore the gap-opening feedback on disk thermal structure. We carry out parameter studies by considering different planet locations rp and planet masses Mp. We find that for the same rp and Mp, our iteration method predicts a wider and deeper gap than the non-iteration method. We also find that the inner disk and gap temperature from the iteration method can vary strongly from the non-iteration or disk without planets, which can further influence dust-trap conditions, iceline locations, and distribution of various ices, such as H2O, CO2, and CO on large dust grains ("pebbles"). Through that, a gap-opening planet can complicate the canonical picture of the non-planet disk C/O ratio and influence the composition of the next generation of planetesimals and planets., Comment: 17 pages, 15 figures, accepted for publication in MNRAS

We present an analysis of the JWST/MIRI-MRS spectrum of Sz114, an accreting M5 star surrounded by a large dust disk with a gap at ~39au. The spectrum is molecular-rich: we report the detection of water, CO, CO2, HCN, C2H2, and H2. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water lines with the 17.22 micron emission surpassing that of most late M-star disks by an order of magnitude in flux and aligning instead with disks of earlier-type stars. Moreover, flux rations of C2H2/H2O and HCN/H2O in Sz114 also resemble those of earlier-type disks, with a slightly elevated CO2/H2O ratio. While accretion heating can boost all infrared lines, the unusual properties of Sz114 could be explained by the young age of the source, its formation under unusual initial conditions (a large massive disk), and/or the presence of dust substructures. The latter delay the inward drift of icy pebbles and thus preserve a lower C/O ratio over an extended period. In contrast, late M-star disks-which are typically faint, small in size, and likely lack significant substructures-may have more quickly depleted the outer icy reservoir and already evolved out of a water-rich inner disk phase. Our findings underscore the unexpected diversity within mid-infrared spectra of late M-star disks, highlighting the need to expand the observational sample for a comprehensive understanding of their variations and thoroughly test pebble drift and planet formation models., Comment: 15 pages, 7 figures, submitted to ApJL

While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A., Comment: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

Planet formation models rely on knowledge of the physical conditions and evolutionary processes in protoplanetary disks, in particular the grain size distribution and dust growth timescales. In theoretical models, several barriers exist that prevent grain growth to pebble sizes and beyond, such as the radial drift and fragmentation. Pressure bumps have been proposed to overcome such barriers. In the past decade ALMA has revealed observational evidence for the existence of such pressure bumps in the form of dust traps, such as dust rings, gaps, cavities and crescents through high-resolution millimeter continuum data originating from thermal dust emission of pebble-sized dust grains. These substructures may be related to young protoplanets, either as the starting point or the consequence of early planet formation. Furthermore, disk dust masses have been measured for complete samples of young stars in clusters, which provide initial conditions for the solid mass budget available for planet formation. However, observational biases exist in the selection of high-resolution ALMA observations and uncertainties exist in the derivation of the disk dust mass, which both may affect the observed trends. This chapter describes the latest insights in dust evolution and disk continuum observations. Specifically, disk populations and evolutionary trends are described, as well as the uncertainties therein, and compared with exoplanet demographics., Comment: Submitted version, to be published in: Handbook of Exoplanets, 2nd Edition, Hans Deeg and Juan Antonio Belmonte (Eds. in Chief), Springer International Publishing AG, part of Springer Nature. Comments welcome. Comments welcome

Diverse protoplanetary disk morphology can result from planet-disk interaction, suggesting planetary presence. To date, most scattered light imaging campaigns have probed polarized light, which is only a fraction of the total light and not very sensitive to planets. To observe and characterize protoplanetary disk systems in the near-infrared in both polarized and total intensity light, we carried out an unprecedented study of scattering properties of disks, as well as of any planetary companions. Using SPHERE with star-hopping at the Very Large Telescope, we observed 29 disk hosts and their reference stars in $K_s$-band polarized light. We extracted disks in total intensity by adopting the data imputation concept with sequential non-negative matrix factorization (DI-sNMF). We obtained high-quality disk images in total intensity for 15 systems and in polarized light for 23. For well-recovered disks in polarized light and total intensity, we parameterized the polarization fraction phase functions using scaled beta distribution: the peak of polarization fraction tentatively correlates with the peak scattering angle, which could be reproduced using certain compact dust, yet more detailed modeling studies are needed. We investigated the empirical DI-sNMF detectability of disks using logistic regression: total intensity detectability of disks primarily depends on host star brightness. For disks with SPHERE data in $Y$-/$J$-/$H$-band, we summarized their polarized color at ~90 deg scattering angle: most of disks are blue in polarized $J-K_s$ color, and they are relatively redder as stellar luminosity increases, indicating larger scatterers. High-quality disk imagery in both total intensity and polarized light thus allows for disk characterization in polarization fraction, and reduces the confusion between disk and planetary signals., Comment: 25 pages, 16 figures, 3 tables, A&A accepted. Data files in FITS format will be publicly available

[abridged]Photoevaporation and dust-trapping are individually considered to be important mechanisms in the evolution and morphology of protoplanetary disks. We studied how the presence of early substructures affects the evolution of the dust distribution and flux in the millimeter continuum of disks that are undergoing photoevaporative dispersal. We also tested if the predicted properties resemble those observed in the population of transition disks. We used the numerical code Dustpy to simulate disk evolution considering gas accretion, dust growth, dust-trapping at substructures, and mass loss due to X-ray and EUV (XEUV) photoevaporation and dust entrainment. Then, we compared how the dust mass and millimeter flux evolve for different disk models. We find that, during photoevaporative dispersal, disks with primordial substructures retain more dust and are brighter in the millimeter continuum than disks without early substructures, regardless of the photoevaporative cavity size. Once the photoevaporative cavity opens, the estimated fluxes for the disk models that are initially structured are comparable to those found in the bright transition disk population ($F_\textrm{mm} > 30\, \textrm{mJy}$), while the disk models that are initially smooth have fluxes comparable to the transition disks from the faint population ($F_\textrm{mm} < 30\, \textrm{mJy}$), suggesting a link between each model and population. Our models indicate that the efficiency of the dust trapping determines the millimeter flux of the disk, while the gas loss due to photoevaporation controls the formation and expansion of a cavity, decoupling the mechanisms responsible for each feature. In consequence, even a planet with a mass comparable to Saturn could trap enough dust to reproduce the millimeter emission of a bright transition disk, while its cavity size is independently driven by photoevaporative dispersal., Comment: Accepted for publication in A&A

PDS 70 is so far the only young disc where multiple planets have been detected by direct imaging. The disc has a large cavity when seen at sub-mm and NIR wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disc shaped by two giant planets, and in particular to test whether disc winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O I] 6300AA line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O I] line luminosity is high with respect to the accretion luminosity when compared to a large sample of discs with cavities in nearby star-forming regions. The FWHM and blue-shifted peak of the [O I] line suggest an emission in a region very close to the star, favouring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He I 10830AA line, which is unusual for low-accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner disc wind could also have had a significant contribution to clearing the inner-disc., Comment: 13 pages, 7 figures, accepted for publication in ApJ

Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles drift inward from the outer disk and sublimate after crossing the snowline. Previous analyses of IRS water spectra, however, were uncertain due to the low spectral resolution that blended lines together. We present new JWST-MIRI spectra of four disks, two compact and two large with multiple radial gaps, selected to test the scenario that water vapor inside the snowline is regulated by pebble drift. The higher spectral resolving power of MIRI-MRS now yields water spectra that separate individual lines, tracing upper level energies from 900 K to 10,000 K. These spectra clearly reveal excess emission in the low-energy lines in compact disks, compared to the large disks, demonstrating an enhanced cool component with $T \approx$ 170-400 K and equivalent emitting radius $R_{\rm{eq}}\approx$ 1-10 au. We interpret the cool water emission as ice sublimation and vapor diffusion near the snowline, suggesting that there is indeed a higher inwards mass flux of icy pebbles in compact disks. Observation of this process opens up multiple exciting prospects to study planet formation chemistry in inner disks with JWST., Comment: Accepted for publication in ApJ Letters

Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H$_2$O) and their abundance within the water snow line region (within a few au). In this study, we aim to understand what effect the presence of gaps in the outer gas disk may have on water vapor enrichment in the inner disk. Building on previous work, we employ a volatile-inclusive multi-Myr disk evolution model that considers an evolving ice-bearing drifting dust population, sensitive to dust-traps, which loses its icy content to sublimation upon reaching the snow line. We find that vapor abundance in the inner disk is strongly affected by fragmentation velocity (v$_{\rm f}$) and turbulence, which control how intense vapor enrichment from pebble delivery is, if present, and how long it may last. Generally, for disks with low to moderate turbulence ($\alpha$ $\le$ 1 $\times$ 10$^{-3}$) and for a range of v$_{\rm f}$, radial location, and gap depth (especially that of the innermost gaps), can significantly alter enrichment. Shallow inner gaps may continuously leak material from beyond it, despite the presence of additional deep outer gaps. We finally find that the for realistic v$_{\rm f}$ ($\le$ 10 m s$^{-1}$), presence of gaps is more important than planetesimal formation beyond the snow line in regulating pebble and volatile delivery into the inner disk., Comment: 24 Pages, 15 Figures, Accepted to ApJ

[Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the steady-state MRI-driven accretion, by employing and improving our previous 1+1D MRI-driven turbulence model. Second, we relax the steady-state accretion assumption in this disk accretion model, and partially couple it to a dust evolution model in order to investigate how the evolution of dust (dynamics and grain growth processes combined) and MRI-driven accretion are intertwined on million-year timescales. Results. Dust coagulation and settling lead to a higher gas ionization degree in the protoplanetary disk, resulting in stronger MRI-driven turbulence as well as a more compact dead zone. On the other hand, fragmentation has an opposite effect because it replenishes the disk in small dust particles. Since the dust content of the disk decreases over million years of evolution due to radial drift, the MRI-driven turbulence overall becomes stronger and the dead zone more compact until the disk dust-gas mixture eventually behaves as a grain-free plasma. Furthermore, our results show that dust evolution alone does not lead to a complete reactivation of the dead zone. Conclusions. The MRI activity evolution (hence the temporal evolution of the MRI-induced $\alpha$-parameter) is controlled by dust evolution and occurs on a timescale of local dust growth, as long as there is enough dust particles in the disk to dominate the recombination process for the ionization chemistry. Once it is no longer the case, it is expected to be controlled by gas evolution and occurs on a viscous evolution timescale., Comment: 23 pages, 13 figures, Accepted for publication in A&A

Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $\sigma$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3\sigma$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5\sigma$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST., Comment: 15 pages, 5 figures. Accepted for publication in ApJ

Circumstellar disks do not evolve in isolation, as about half of solar-type stars were born in binary or multiple systems. Resolving disks in binary systems provides the opportunity to examine the influence of stellar companions on the outcomes of planet formation. We aim to investigate and compare disks in stellar multiple systems with near-infrared scattered-light imaging as part of the Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS) program. We used polarimetric differential imaging with SPHERE/IRDIS at the VLT to search for scattered light from the circumstellar disks in three multiple systems, CHX 22, S CrA, and HP Cha. We performed astrometric and orbit analyses for the stellar companions using archival HST, VLT/NACO, and SPHERE data. Combined with the age and orbital constraints, the observed disk structures provide insights into the evolutionary history and the impact of the stellar companions. The small grains in CHX 22 form a tail-like structure surrounding the close binary, which likely results from a close encounter and capture of a cloudlet. S CrA shows intricate structures (tentative ringed and spiral features) in the circumprimary disk as a possible consequence of perturbations by companions. The circumsecondary disk is truncated and connected to the primary disk via a streamer, suggesting tidal interactions. In HP Cha, the primary disk is less disturbed and features a tenuous streamer, through which the material flows towards the companions. The comparison of the three systems spans a wide range of binary separation (50 - 500 au) and illustrates the decreasing influence on disk structures with the distance of companions. This agrees with the statistical analysis of exoplanet population in binaries, that planet formation is likely obstructed around close binary systems, while it is not suppressed in wide binaries., Comment: 19 pages, 6 figures, accpeted for publication in A&A

Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps and moment maps of the line of sight velocity and peak intensity. We fitted a Keplerian model of the channel-by-channel emission to study line profile differences and produced deprojected radial profiles for all velocity components. The $^{12}$CO emission shows a cavity inwards of $\sim$56 au and within the dust continuum ring at 81 au. The azimuthal brightness variations in the $^{12}$CO line and dust continuum are broadly aligned with the shadows detected in scattered-light observations. We find a strong localized non-Keplerian feature toward the west within the continuum ring (at $R=41\pm 10$ au and $PA=280\pm 2 ^\circ$). A tightly wound spiral is also detected which extends over 300$^\circ$ in azimuth, possibly connected to the localized non-Keplerian feature. Finally, a bending of the iso-velocity contours within the gas cavity indicates a highly perturbed inner region, possibly related to the presence of a misaligned inner disk. While broadly aligned with the scattered-light shadows, the localized non-Keplerian feature cannot be solely due to changes in temperature. Instead, we interpret the kinematical feature as tracing a massive companion located at the edge of the dust continuum ring. We speculate that the spiral is caused by buoyancy resonances driven by planet-disk interactions. However, this potential planet at $\sim$41 au cannot explain the gas-depleted cavity, the low accretion rate, and the misaligned inner disk, which suggests the presence of another companion closer in., Comment: 13 pages, 14 figures. To be published in Astronomy & Astrophysics Letters

Dense stellar environments as hosts of ongoing star formation increase the probability of gravitational encounters among stellar systems during the early stages of evolution. Stellar interaction may occur through non-recurring, hyperbolic or parabolic passages (a so-called 'fly-by'), through secular binary evolution, or through binary capture. In all three scenarios, the strong gravitational perturbation is expected to manifest itself in the disc structures around the individual stars. Here, we present near-infrared polarised light observations that were taken with the SPHERE/IRDIS instrument of three known interacting twin-disc systems: AS 205, EM* SR 24, and FU Orionis. The scattered light exposes spirals likely caused by the gravitational interaction. On a larger scale, we observe connecting filaments between the stars. We analyse their very complex polarised intensity and put particular attention to the presence of multiple light sources in these systems. The local angle of linear polarisation indicates the source whose light dominates the scattering process from the bridging region between the two stars. Further, we show that the polarised intensity from scattering with multiple relevant light sources results from an incoherent summation of the individuals' contribution. This can produce nulls of polarised intensity in an image, as potentially observed in AS 205. We discuss the geometry and content of the systems by comparing the polarised light observations with other data at similar resolution, namely with ALMA continuum and gas emission. Collective observational data can constrain the systems' geometry and stellar trajectories, with the important potential to differentiate between dynamical scenarios of stellar interaction., Comment: 23 pages, 17 figures, accepted for publication in MNRAS

Brown dwarfs and very low mass stars are a significant fraction of stars in our galaxy, and they are interesting laboratories to investigate planet formation in extreme conditions of low temperature and densities. In addition, the dust radial drift of particles is expected to be a more difficult barrier to overcome during the first steps of planet formation in these disks. ALMA high-angular resolution observations of few protoplanetary disks around BDs and VLMS have shown substructures as in the disks around Sun-like stars. Such observations suggest that giant planets embedded in the disks are the most likely origin of the observed substructures. However, this type of planets represent less than 2% of the confirmed exoplanets so far around all stars, and they are difficult to form by different core accretion models (either pebble or planetesimal accretion). Dedicated deep observations of disks around BDs and VLMS with ALMA and JWST will provide significant progress on understanding the main properties of these objects (e.g., disk size and mass), which is crucial for determining the physical mechanisms that rule the evolution of these disks and the effect on the potential planets that may form in these environments., Comment: Invited contribution as an article in the European Physics Journal + (accepted for publication)

Over the last years, large (sub-)millimetre surveys of protoplanetary disks have well constrained the demographics of disks, such as their millimetre luminosities, spectral indices, and disk radii. Additionally, several high-resolution observations have revealed an abundance of substructures in the disks dust continuum. The most prominent are ring like structures, likely due to pressure bumps trapping dust particles. The origins and characteristics of these bumps, nevertheless, need to be further investigated. The purpose of this work is to study how dynamic pressure bumps affect observational properties of protoplanetary disks. We further aim to differentiate between the planetary- versus zonal flow-origin of pressure bumps. We perform one-dimensional gas and dust evolution simulations, setting up models with varying pressure bump features. We subsequently run radiative transfer calculations to obtain synthetic images and the different quantities of observations. We find that the outermost pressure bump determines the disks dust size across different millimetre wavelengths. Our modelled dust traps need to form early (< 0.1 Myr), fast (on viscous timescales), and must be long lived (> Myr) to obtain the observed high millimetre luminosities and low spectral indices of disks. While the planetary bump models can reproduce these observables irrespectively of the opacity prescription, the highest opacities are needed for the zonal flow bump model to be in line with observations. Our findings favour the planetary- over the zonal flow-origin of pressure bumps and support the idea that planet formation already occurs in early class 0-1 stages of circumstellar disks. The determination of the disks effective size through its outermost pressure bump also delivers a possible answer to why disks in recent low-resolution surveys appear to have the same sizes across different millimetre wavelengths., Comment: 22 pages, 15 figures. To be published in Astronomy & Astrophysics

The radial velocities and direct imaging observations of exoplanets have suggested that the frequency of giant planets may decrease for intermediate-mass stars ($2.5-8\,M_\odot$). The key mechanism that could hinder their formation remains unclear. From a theoretical point of view, planet formation around intermediate-mass stars may take place on longer timescales, which -- coupled with fast migration and efficient photoevaporation -- may prevent planetary formation in these environments. In this letter, we investigate the temporal evolution of the radial drift for dust particles in disks when stellar evolution is taken into account. We demonstrate that the particle drift velocity around intermediate-mass stars sharply increases after 1$-$2 Myr, potentially forming a difficult barrier to overcome in the first steps of planet formation. This high radial drift could explain the lack of disk detections around intermediate-mass stars older than 3$-$4 Myr, as opposed to low-mass stars ($<2.5\,M_\odot$), where the drift may not be the most impactful factor for the disk evolution. Future high-resolution images of these disks can help us to explain why planets around intermediate-mass stars may be rare. In addition, we can explore whether the role of efficient dust radial drift does in fact hinder planet formation around intermediate-mass stars -- or otherwise., Comment: Accepted for publication in A&A (Letters). Minor changes after language edition

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales., Comment: 42 pages, 20 figures, accepted by ApJ, reduced data available at https://zenodo.org/record/6408903 (typos from v1 fixed)

The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emission span from 50 to 1000\,au. This range could be explained by viscous evolution models with different $\alpha$ values (mostly of $10^{-4}-10^{-3}$) and/or a spread of initial conditions. The CO sizes for most disks are also consistent with MHD wind models that directly remove disk angular momentum, but very large initial disk sizes would be required to account for the very extended CO disks in the sample. As no CO size evolution is observed across stellar ages of 0.5--20\,Myr in this sample, determining the dominant mechanism of disk evolution will require a more complete sample for both younger and more evolved systems. We find that the CO emission is universally more extended than the continuum emission by an average factor of $2.9\pm1.2$. The ratio of the CO to continuum sizes does not show any trend with stellar mass, mm continuum luminosity, or the properties of substructures. The GO Tau disk has the most extended CO emission in this sample, with an extreme CO to continuum size ratio of 7.6. Seven additional disks in the sample show high size ratios ($\gtrsim4$) that we interpret as clear signs of substantial radial drift., Comment: Accepted for publication in ApJ

Observations of young stars hosting transition disks show that several of them have high accretion rates, despite their disks presenting extended cavities in their dust component. This represents a challenge for theoretical models, which struggle to reproduce both features. We explore if a disk evolution model, including a dead zone and disk dispersal by X-ray photoevaporation, can explain the high accretion rates and large gaps (or cavities) measured in transition disks. We implement a dead zone turbulence profile and a photoevaporative mass loss profile into numerical simulations of gas and dust. We perform a population synthesis study of the gas component, and obtain synthetic images and SED of the dust component through radiative transfer calculations. This model results in long lived inner disks and fast dispersing outer disks, that can reproduce both the accretion rates and gap sizes observed in transition disks. For a dead zone of turbulence $\alpha_{dz} = 10^{-4}$ and extent $r_{dz}$ = 10 AU, our population synthesis study shows that $63\%$ of our transition disks are accreting with $\dot{M}_g > 10^{-11} M_\odot/yr$ after opening a gap. Among those accreting transition disks, half display accretion rates higher than $5\times10^{-10} M_\odot/yr$ . The dust component in these disks is distributed in two regions: in a compact inner disk inside the dead zone, and in a ring at the outer edge of the photoevaporative gap, which can be located between 20 AU and 100 AU. Our radiative transfer calculations show that the disk displays an inner disk and an outer ring in the millimeter continuum, a feature observed in some transition disks. A disk model considering X-ray photoevaporative dispersal in combination with dead zones can explain several of the observed properties in transition disks including: the high accretion rates, the large gaps, and long-lived inner disks at mm-emission., Comment: 16 pages, 15 figures. Accepted in A&A

The mass of protoplanetary disks is arguably one of their most important quantities shaping their evolution toward planetary systems, but it remains a challenge to determine this quantity. Using the high spatial resolution now available on telescopes such as the Atacama Large Millimeter/submillimeter Array (ALMA), recent studies derived a relation between the disk surface density and the location of the "dust lines". This is a new concept in the field, linking the disk size at different continuum wavelengths with the radial distribution of grain populations of different sizes. We aim to use a dust evolution model to test the dependence of the dust line location on disk gas mass. In particular, we are interested in the reliability of the method for disks showing radial substructures, as recent high-resolution observations revealed. Our models show that the determination of the dust line location is a promising approach to the mass estimate of protoplanetay disks, but the exact relation between the dust line location and disk mass depends on the structure of the particular disk. We calibrated the relation for disks without evidence of radial structures, while for more complex structures we ran a simple dust evolution model. However, this method fails when there is evidence of strong dust traps. It is possible to reveal when dust evolution is dominated by traps, providing the necessary information for when the method should be applied with caution., Comment: 14 pages, 18 figures, accepted by Astronomy & Astrophysics

[abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $\alpha$ is determined self-consistently under the framework of viscously-driven accretion from detailed considerations of the MRI with non-ideal MHD effects (Ohmic resistivity and ambipolar diffusion), accounting for disk heating by stellar irradiation, non-thermal sources of ionization, and dust effects on the ionization chemistry. Additionally, the magnetic field strength is constrained and adopted to maximize the MRI activity. We demonstrate the use of our framework by investigating steady-state MRI-driven accretion in a fiducial protoplanetary disk model around a solar-type star. We find that the equilibrium solution displays no pressure maximum at the dead zone outer edge, except if a sufficient amount of dust particles have accumulated there before the disk reaches a steady-state accretion regime. Furthermore, the steady-state accretion solution describes a disk that displays a spatially extended long-lived inner disk gas reservoir (the dead zone) accreting a few $10^{-9}\, M_{\odot}.\rm{yr}^{-1}$. By conducting a detailed parameter study, we find that the extend to which the MRI can drive efficient accretion is primarily determined by the total disk gas mass, the representative grain size, the vertically-integrated dust-to-gas mass ratio, and the stellar X-ray luminosity. A self-consistent time-dependent coupling between gas, dust, stellar evolution models and our general framework on million-year timescales is required to fully understand the formation of dead zones and their potential to trap dust particles., Comment: Accepted for publication in A&A

Context. The protoplanetary disk around the star HD 100546 displays prominent substructures in the form of two concentric rings. Recent observations with the Atacama Large Millimeter/sub-millimeter Array (ALMA) have revealed these features with high angular resolution and have resolved the faint outer ring well. This allows us to study the nature of the system further. Aims. Our aim is to constrain some of the properties of potential planets embedded in the disk, assuming that they induce the observed rings and gaps. Methods. We present the self-calibrated $0.9\,$mm ALMA observations of the dust continuum emission from the circumstellar disk around HD 100546. These observations reveal substructures in the disk that are consistent with two rings, the outer ring being much fainter than the inner one. We reproduced this appearance closely with a numerical model that assumes two embedded planets. We varied planet and disk parameters in the framework of the planet-disk interaction code FARGO3D and used the outputs for the gas and dust distribution to generate synthetic observations with the code RADMC-3D. Results. From this comparison, we find that an inner planet located at $r_1 = 13\,$au with a mass $M_1 = 8 M_{\rm{Jup}}$ and an outer planet located at $r_2 = 143\,$au with a mass $M_2 = 3 M_{\rm{Jup}}$ leads to the best agreement between synthetic and ALMA observations (deviation less than $3\sigma$ for the normalized radial profiles). To match the very low brightness of the outer structure relative to the inner ring, the initial disk gas surface density profile needs to follow an exponentially tapered power law (self-similar solution), rather than a simple power-law profile., Comment: Accepted for publication in A&A

PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$\alpha$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 855$\mu$m at high angular resolution ($\sim$20mas, 2.3au) that aim to resolve the circumplanetary disks and constrain their dust masses. Our observations confirm the presence of a compact source of emission co-located with PDS70c, spatially separated from the circumstellar disk and less extended than $\sim$1.2au in radius, a value close to the expected truncation radius of the cicumplanetary disk at a third of the Hill radius. The emission around PDS70c has a peak intensity of $\sim$86$\pm$16 $\mu \mathrm{Jy}~\mathrm{beam}^{-1}$ which corresponds to a dust mass of $\sim$0.031M$_{\oplus}$ or $\sim$0.007M$_{\oplus}$, assuming that it is only constituted of 1 $\mu$m or 1 mm sized grains, respectively. We also detect extended, low surface brightness continuum emission within the cavity near PDS70b. We observe an optically thin inner disk within 18au of the star with an emission that could result from small micron-sized grains transported from the outer disk through the orbits of b and c. In addition, we find that the outer disk resolves into a narrow and bright ring with a faint inner shoulder., Comment: ApJ Letters, in press; 19 pages, 9 figures

Recent ALMA observations have found many protoplanetary discs with rings that can be explained by gap-opening planets less massive than Jupiter. Meanwhile, recent studies have suggested that protoplanetary discs should have low levels of turbulence. Past computational work on low-viscosity discs has hinted that these two developments might not be self-consistent because even low-mass planets can be accompanied by vortices instead of conventional double rings. We investigate this potential discrepancy by conducting hydrodynamic simulations of growing planetary cores in discs with various aspect ratios ($H/r=0.04$, 0.06, 0.08) and viscosities ($1.5 \times 10^{-5} \lesssim \alpha \lesssim 3 \times 10^{-4}$), having these cores accrete their gas mass directly from the disc. With $\alpha < 10^{-4}$, we find that sub-Saturn-mass planets in discs with $H/r \le 0.06$ are more likely to be accompanied by dust asymmetries compared to Jupiter-mass planets because they can trigger several generations of vortices in succession. We also find that vortices with $H/r = 0.08$ survive $>6000$ planet orbits regardless of the planet mass or disc mass because they are less affected by the planet's spiral waves. We connect our results to observations and find that the outward migration of vortices with $H/r \ge 0.08$ may be able to explain the cavity in Oph IRS 48 or the two clumps in MWC 758. Lastly, we show that the lack of observed asymmetries in the disc population in Taurus is unexpected given the long asymmetry lifetimes in our low viscosity simulations ($\alpha \sim 2 \times 10^{-5}$), a discrepancy we suggest is due to these discs having higher viscosities., Comment: 24 pages, 20 figures, 2 tables; Accepted to MNRAS

We aim to estimate if structures, such as cavities, rings, and gaps, are common in disks around VLMS and to test models of structure formation in these disks. We also aim to compare the radial extent of the gas and dust emission in disks around VLMS, which can give us insight about radial drift. We studied six disks around VLMS in the Taurus star-forming region using ALMA Band 7 ($\sim 340\,$GHz) at a resolution of $\sim0.1''$. The targets were selected because of their high disk dust content in their stellar mass regime. Our observations resolve the disk dust continuum in all disks. In addition, we detect the $^{12}$CO ($J=3-2$) emission line in all targets and $^{13}$CO ($J=3-2$) in five of the six sources. The angular resolution allows the detection of dust substructures in three out of the six disks, which we studied by using UV-modeling. Central cavities are observed in the disks around stars MHO\,6 (M5.0) and CIDA\,1 (M4.5), while we have a tentative detection of a multi-ringed disk around J0433. Single planets of masses $0.1\sim0.4\,M_{\rm{Jup}}$ would be required. The other three disks with no observed structures are the most compact and faintest in our sample. The emission of $^{12}$CO and $^{13}$CO is more extended than the dust continuum emission in all disks of our sample. When using the $^{12}$CO emission to determine the gas disk extension $R_{\rm{gas}}$, the ratio of $R_{\rm{gas}}/R_{\rm{dust}}$ in our sample varies from 2.3 to 6.0, which is consistent with models of radial drift being very efficient around VLMS in the absence of substructures. Our observations do not exclude giant planet formation on the substructures observed. A comparison of the size and luminosity of VLMS disks with their counterparts around higher mass stars shows that they follow a similar relation.

We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $<$10 au). We find that two discs around embedded protostars lack the clear gaps and rings that are ubiquitous in more evolved sources with Class II SEDs. Our sample includes 5 objects with previously known large dust cavities (r $>$20 au). We find that the 1.3 mm radial profiles of these objects are in good agreement with those produced by numerical simulations of dust evolution and planet-disc interactions, which predict the accumulation of mm-sized grains at the edges of planet-induced cavities. Our long-baseline observations resulted in the largest sample of discs observed at $\sim$3-5 au resolution in any given star-forming region (15 objects when combined with Ophiuchus objects in the DSHARP Large Program) and allow for a demographic study of the brightest $\sim5\%$ of the discs in Ophiuchus (i.e. the most likely formation sites of giant planets in the cloud). We use this unique sample to propose an evolutionary sequence and discuss a scenario in which the substructures observed in massive protoplanetary discs are mainly the result of planet formation and dust evolution. If this scenario is correct, the detailed study of disc substructures might provide a window to investigate a population of planets that remains mostly undetectable by other techniques., Comment: 21 pages, 10 figures. Appendix with 3 additional figures. Version 2, identical to previous one, but now accepted for publication (in MNRAS)

[abridged] Recent laboratory experiments indicate that destructive collisions of icy dust particles occur with much lower velocities than previously thought. When these new velocities are considered from laboratory experiments in dust evolution models, a growth to pebble sizes in protoplanetary disks (PPDs) is difficult. This may contradict (sub-)mm observations and challenge the formation of planetesimals and planets. We investigate the conditions that are required in dust evolution models for growing and trapping pebbles in PPDs when the fragmentation speed is 1ms$^{-1}$ in the entire disk. We distinguish the parameters controlling the effects of turbulent velocities, vertical stirring, radial diffusion, and gas viscous evolution, always assuming that particles cannot diffuse faster (radially or vertically) than the gas. To form pebbles and produce effective particle trapping, the parameter that controls the particle turbulent velocities must be small ($\delta_t\lesssim10^{-4}$). In these cases, the vertical settling can limit the formation of pebbles, which also prevents particle trapping. Therefore the parameter that sets the vertical settling of the grains must be $\delta_z<10^{-3}$. Our results suggest that different combinations of the particle and gas diffusion parameters can lead to a large diversity of millimeter fluxes and dust-disk radii. When pebble formation occurs and trapping is efficient, gaps and rings have higher contrast at mm-emission than in the NIR. In the case of inefficient trapping, structures are also formed at the two wavelengths, producing deeper and wider gaps in the NIR. Our results highlight the importance of obtaining observational constraints of gas and particle diffusion parameters and the properties of gaps at short and long wavelengths to better understand basic features of PPDs and the origin of the structures that are observed in these objects., Comment: Accepted for publication in A&A

ISO-Oph 2 is a wide-separation (240 au) binary system where the primary star harbors a massive (M$_{dust}$ $\sim$40 M$_{\oplus}$) ring-like disk with a dust cavity $\sim$50 au in radius and the secondary hosts a much lighter (M$_{dust}$ $\sim$0.8 M$_{\oplus}$) disk. As part of the high-resolution follow-up of the "Ophiuchus Disk Survey Employing ALMA" (ODISEA) project, we present 1.3 mm continuum and $^{12}$CO molecular line observations of the system at 0''02 (3 au) resolution. We resolve the disk around the primary into two non-axisymmetric rings and find that the disk around the secondary is only $\sim$7 au across and also has a dust cavity (r $\sim$2.2 au). Based on the infrared flux ratio of the system and the M0 spectral type of the primary, we estimate the mass of the companion to be close to the brown dwarf limit. Hence, we conclude that the ISO-Oph 2 system contains the largest and smallest cavities, the smallest measured disk size, and the resolved cavity around the lowest mass object (M$_{\star}$ $\sim$0.08 M$_\odot$) in Ophiuchus. From the $^{12}$CO data, we find a bridge of gas connecting both disks. While the morphology of the rings around the primary might be due to an unseen disturber within the cavity, we speculate that the bridge might indicate an alternative scenario in which the secondary has recently flown by the primary star causing the azimuthal asymmetries in its disk. The ISO-Oph 2 system is therefore a remarkable laboratory to study disk evolution, planet formation, and companion-disk interactions., Comment: 10 pages, 4 figures, 1 table

We present a synergic study of protoplanetary disks to investigate links between inner disk gas molecules and the large-scale migration of solid pebbles. The sample includes 63 disks where two types of measurements are available: i) spatially-resolved disk images revealing the radial distribution of disk pebbles (mm-cm dust grains), from millimeter observations with ALMA or the SMA, and ii) infrared molecular emission spectra as observed with Spitzer. The line flux ratios of H2O with HCN, C2H2, and CO2 all anti-correlate with the dust disk radius R$_{dust}$, expanding previous results found by Najita et al. (2013) for HCN/H2O and the dust disk mass. By normalization with the dependence on accretion luminosity common to all molecules, only the H2O luminosity maintains a detectable anti-correlation with disk radius, suggesting that the strongest underlying relation is between H2O and R$_{dust}$. If R$_{dust}$ is set by large-scale pebble drift, and if molecular luminosities trace the elemental budgets of inner disk warm gas, these results can be naturally explained with scenarios where the inner disk chemistry is fed by sublimation of oxygen-rich icy pebbles migrating inward from the outer disk. Anti-correlations are also detected between all molecular luminosities and the infrared index n$_{13-30}$, which is sensitive to the presence and size of an inner disk dust cavity. Overall, these relations suggest a physical interconnection between dust and gas evolution both locally and across disk scales. We discuss fundamental predictions to test this interpretation and study the interplay between pebble drift, inner disk depletion, and the chemistry of planet-forming material., Comment: Accepted for publication on ApJ

The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $\mu$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a part of the disk. The polarization vectors are parallel to the disk minor axis, and the polarization fraction was derived to be $1-2$\%. These characteristics are consistent with models of self-scattering of submillimeter-wave emission, which indicate a maximum grain size of $\sim100$ $\mu$m. However, this is inconsistent with the previous interpretation of millimeter/centimeter dust particles being trapped by a vortex. To explain both, ALMA polarization and previous ALMA and VLA observations, we suggest that the thermal emission at 860 $\mu$m wavelength is optically thick ($\tau_{\rm abs}\sim7.3$) at the dust trap with the maximum observable grain size of $\sim100$ $\mu$m rather than an optically thin case with $\sim$ cm dust grains. We note that we cannot rule out that larger dust grains are accumulated near the midplane if the 860 $\mu$m thermal emission is optically thick., Comment: 22 pages, 17 figures, accepted for publication in ApJ

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($\alpha_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($\tau\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation., Comment: accepted for publication in ApJ

Demographic surveys of protoplanetary disks, carried out mainly with ALMA, have provided access to a large range of disk dust masses ($M_{\rm{dust}}$) around stars with different stellar types and in different star-forming regions. These surveys found a power-law relation between $M_{\rm{dust}}$ and $M_{\star}$ that steepens in time, but which is also flatter for transition disks (TDs). We performed dust evolution models, which included perturbations to the gas surface density with different amplitudes to investigate the effect of particle trapping on the $M_{\rm{dust}}-M_{\star}$ relation. These perturbations were aimed at mimicking pressure bumps that originated from planets. We focused on the effect caused by different stellar and disk masses based on exoplanet statistics that demonstrate a dependence of planet mass on stellar mass and metallicity. Models of dust evolution can reproduce the observed $M_{\rm{dust}}-M_{\star}$ relation in different star-forming regions when strong pressure bumps are included and when the disk mass scales with stellar mass (case of $M_{\rm{disk}}=0.05\,M_\star$ in our models). This result arises from dust trapping and dust growth beyond centimeter-sized grains inside pressure bumps. However, the flatter relation of $M_{\rm{dust}}-M_{\star}$ for TDs and disks with substructures cannot be reproduced by the models unless the formation of boulders is inhibited inside pressure bumps. In the context of pressure bumps originating from planets, our results agree with current exoplanet statistics on giant planet occurrence increasing with stellar mass, but we cannot draw a conclusion about the type of planets needed in the case of low-mass stars. This is attributed to the fact that for $M_\star<1\,M_\odot$, the observed $M_{\rm{dust}}$ obtained from models is very low due to the efficient growth of dust particles beyond centimeter-sizes inside pressure bumps., Comment: Accepted for publication in A&A, minor changes after language edition

We utilize ALMA archival data to estimate the dust disk size of 152 protoplanetary disks in Lupus (1-3 Myr), Chamaeleon I (2-3 Myr), and Upper-Sco (5-11 Myr). We combine our sample with 47 disks from Tau/Aur and Oph whose dust disk radii were estimated, as here, through fitting radial profile models to visibility data. We use these 199 homogeneously derived disk sizes to identify empirical disk-disk and disk-host property relations as well as to search for evolutionary trends. In agreement with previous studies, we find that dust disk sizes and millimeter luminosities are correlated, but show for the first time that the relationship is not universal between regions. We find that disks in the 2-3 Myr-old Cha I are not smaller than disks in other regions of similar age, and confirm the Barenfeld et al. (2017) finding that the 5-10 Myr USco disks are smaller than disks belonging to younger regions. Finally, we find that the outer edge of the Solar System, as defined by the Kuiper Belt, is consistent with a population of dust disk sizes which have not experienced significant truncation., Comment: ApJ accepted, 38 pages, 16 figures, 68k compatible

We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~\mu$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths., Comment: ApJ accepted, 25 pages, 11 figures

We report new Atacama Large Millimeter/sub-millimeter Array (ALMA) Band 3 observations at 2.75 mm of the TD around SR 24S with an angular resolution of $\sim$0.11''$\times$ 0.09'' and a peak signal-to-noise ratio of $\sim24$. We detect an inner disk and a mostly symmetric ring-like structure that peaks at $\sim$0.32'', that is $\sim$37 au at a distance of $\sim$114.4 pc. The full width at half maximum of this ring is $\sim$28 au. We analyze the observed structures by fitting the dust continuum visibilities using different models for the intensity profile, and compare with previous ALMA observations of the same disk at 0.45 mm and 1.30 mm. We qualitatively compare the results of these fits with theoretical predictions of different scenarios for the formation of a cavity or large gap. The comparison of the dust continuum structure between different ALMA bands indicates that photoevaporation and dead zone can be excluded as leading mechanisms for the cavity formation in SR 24S disk, leaving the planet scenario (single or multiple planets) as the most plausible mechanism. We compared the 2.75 mm emission with published (sub-)centimeter data and find that the inner disk is likely tracing dust thermal emission. This implies that any companion in the system should allow dust to move inwards throughout the gap and replenish the inner disk. In the case of one single planet, this puts strong constraints on the mass of the potential planet inside the cavity and the disk viscosity of about $\lesssim$5 $M_{\rm{Jup}}$ and $\alpha\sim10^{-4}-10^{-3}$, respectively., Comment: Accepted to ApJ

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet., Comment: Accepted for publication in A&A, 11 pages, 8 figures