No Arabic abstract
We present the first high-resolution sub-mm survey of both dust and gas for a large population of protoplanetary disks. Characterizing fundamental properties of protoplanetary disks on a statistical level is critical to understanding how disks evolve into the diverse exoplanet population. We use ALMA to survey 89 protoplanetary disks around stars with $M_{ast}>0.1~M_{odot}$ in the young (1--3~Myr), nearby (150--200~pc) Lupus complex. Our observations cover the 890~$mu$m continuum and the $^{13}$CO and C$^{18}$O 3--2 lines. We use the sub-mm continuum to constrain $M_{rm dust}$ to a few Martian masses (0.2--0.4~$M_{oplus}$) and the CO isotopologue lines to constrain $M_{rm gas}$ to roughly a Jupiter mass (assuming ISM-like $rm {[CO]/[H_2]}$ abundance). Of 89 sources, we detect 62 in continuum, 36 in $^{13}$CO, and 11 in C$^{18}$O at $>3sigma$ significance. Stacking individually undetected sources limits their average dust mass to $lesssim6$ Lunar masses (0.03~$M_{oplus}$), indicating rapid evolution once disk clearing begins. We find a positive correlation between $M_{rm dust}$ and $M_{ast}$, and present the first evidence for a positive correlation between $M_{rm gas}$ and $M_{ast}$, which may explain the dependence of giant planet frequency on host star mass. The mean dust mass in Lupus is 3$times$ higher than in Upper Sco, while the dust mass distributions in Lupus and Taurus are statistically indistinguishable. Most detected disks have $M_{rm gas}lesssim1~M_{rm Jup}$ and gas-to-dust ratios $<100$, assuming ISM-like $rm {[CO]/[H_2]}$ abundance; unless CO is very depleted, the inferred gas depletion indicates that planet formation is well underway by a few Myr and may explain the unexpected prevalence of super-Earths in the exoplanet population.
We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply Keplerian masking to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys.
We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.
We perform a comprehensive demographic study of the CO extent relative to dust of the disk population in the Lupus clouds, in order to find indications of dust evolution and possible correlations with other properties. We increase up to 42 the number of disks of the region with measured CO and dust sizes ($R_{mathrm{CO}}$, $R_{mathrm{dust}}$) from observations with the Atacama Large Millimeter/submillimeter Array (ALMA). The sizes are obtained from modeling the ${^{12}}$CO $J = 2-1$ line emission and continuum emission at $sim 0.89$ mm with an empirical function (Nuker profile or Gaussian function). The CO emission is more extended than the dust continuum, with a $R_{68%}^{mathrm{CO}}$/$R_{68%}^{mathrm{dust}}$ median value of 2.5, for the entire population and for a sub-sample with high completeness. 6 disks, around $15%$ of the Lupus disk population have a size ratio above 4. Based on thermo-chemical modeling, this value can only be explained if the disk has undergone grain growth and radial drift. These disks do not have unusual properties in terms of stellar mass ($M_{star}$), disk mass ($M_{mathrm{disk}}$), CO and dust sizes ($R_{mathrm{CO}}$, $R_{mathrm{dust}}$), and mass accretion. We search for correlations between the size ratio and $M_{star}$, $M_{mathrm{disk}}$, $R_{mathrm{CO}}$ and $R_{mathrm{dust}}$: only a weak monotonic anti-correlation with the $R_{mathrm{dust}}$ is found. The lack of strong correlations is remarkable and suggests that the bulk of the population may be in a similar evolutionary stage, independent of the stellar and disk properties. These results should be further investigated, since the optical depth difference between CO and dust continuum may play a role in the inferred size ratios. Lastly, the CO emission for the majority of the disks is consistent with optically thick emission and an average CO temperature of around 30 K.
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 ($tausim$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.
Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $alpha simeq 10^{-4}$. The value of $alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.