Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences.

The chemical composition of planets is inherited by the distribution of the various molecular species in the protoplanetary disk at thetime of their formation. As of today, only a handful of disks has been imaged in multiple spectral lines with high spatial resolution. As part of a small campaign devoted to the chemical characterization of disk-outflow sources in Taurus, we report on new ALMA Band 6 observations with 20 au resolution toward the embedded young star DG Tau B. Images of the continuum emission reveals a dust disk with rings and, putatively, a leading spiral arm. The disk, as well as the prominent outflow cavities, are detected in CO, H2CO, CS, and CN while they remain undetected in SO2, HDO, and CH3OH. From the absorption of the back-side outflow, we inferred that the disk emission is optically thick in the inner 50 au. This morphology explains why no line emission is detected from this inner region and poses some limitations toward the calculation of the dust mass and the characterization of the inner gaseous disk. The H2CO and CS emission from the inner 200 au is mostly from the disk and their morphology is very similar. The CN emission significantly differs from the other two molecules as it is observed only beyond 150 au. This ring-like morphology is consistent with previous observations and the predictions of thermochemical disk models. Finally, we constrained the disk-integrated column density of all molecules. In particular, we found that the CH3OH/H2CO ratio must be smaller than 2, making the methanol non-detection still consistent with the only such a ratio available from the literature (1.27 in TW Hya).

The UV radiation field is a critical regulator of gas-phase chemistry in surface layers of disks around young stars. In an effort to understand the relationship between photocatalyzing UV radiation fields and gas emission observed at infrared and sub-mm wavelengths, we present an analysis of new and archival HST, Spitzer, ALMA, IRAM, and SMA data for five targets in the Lupus cloud complex and 14 systems in Taurus-Auriga. The HST spectra were used to measure LyA and FUV continuum fluxes reaching the disk surface, which are responsible for dissociating relevant molecular species (e.g. HCN, N2). Semi-forbidden C II] 2325 and UV-fluorescent H2 emission were also measured to constrain inner disk populations of C+ and vibrationally excited H2. We find a significant positive correlation between 14 micron HCN emission and fluxes from the FUV continuum and C II] 2325, consistent with model predictions requiring N2 photodissociation and carbon ionization to trigger the main CN/HCN formation pathways. We also report significant negative correlations between sub-mm CN emission and both C II] and FUV continuum fluxes, implying that CN is also more readily dissociated in disks with stronger FUV irradiation. No clear relationships are detected between either CN or HCN and LyA or UV-H2 emission. This is attributed to the spatial stratification of the various molecular species, which span several vertical layers and radii across the inner and outer disk. We expect that future observations with JWST will build on this work by enabling more sensitive IR surveys than were possible with Spitzer.

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time.

We present new 890 $\mu m$ continuum ALMA observations of 5 brown dwarfs (BDs) with infrared excess in Lupus I and III -- which, in combination with 4 BDs previously observed, allowed us to study the mm properties of the full known BD disk population of one star-forming region. Emission is detected in 5 out of the 9 BD disks. Dust disk mass, brightness profiles and characteristic sizes of the BD population are inferred from continuum flux and modeling of the observations. Only one source is marginally resolved, allowing for the determination of its disk characteristic size. We conduct a demographic comparison between the properties of disks around BDs and stars in Lupus. Due to the small sample size, we cannot confirm or disprove if the disk mass over stellar mass ratio drops for BDs, as suggested for Ophiuchus. Nevertheless, we find that all detected BD disks have an estimated dust mass between 0.2 and 3.2 $M_{\bigoplus}$; these results suggest that the measured solid masses in BD disks can not explain the observed exoplanet population, analogous to earlier findings on disks around more massive stars. Combined with the low estimated accretion rates, and assuming that the mm-continuum emission is a reliable proxy for the total disk mass, we derive ratios of $\dot{M}_{\mathrm{acc}} / M_{\mathrm{disk}}$ significantly lower than in disks around more massive stars. If confirmed with more accurate measurements of disk gas masses, this result could imply a qualitatively different relationship between disk masses and inward gas transport in BD disks.

While planets are commonly discovered around main-sequence stars, the processes leading to their formation are still far from being understood. Current planet population synthesis models, which aim to describe the planet formation process from the protoplanetary disk phase to the time exoplanets are observed, rely on prescriptions for the underlying properties of protoplanetary disks where planets form and evolve. The recent development in measuring disk masses and disk-star interaction properties, i.e., mass accretion rates, in large samples of young stellar objects demand a more careful comparison between the models and the data. We performed an initial critical assessment of the assumptions made by planet synthesis population models by looking at the relation between mass accretion rates and disk masses in the models and in the currently available data. We find that the currently used disk models predict mass accretion rate in line with what is measured, but with a much lower spread of values than observed. This difference is mainly because the models have a smaller spread of viscous timescales than what is needed to reproduce the observations. We also find an overabundance of weakly accreting disks in the models where giant planets have formed with respect to observations of typical disks. We suggest that either fewer giant planets have formed in reality or that the prescription for planet accretion predicts accretion on the planets that is too high. Finally, the comparison of the properties of transition disks with large cavities confirms that in many of these objects the observed accretion rates are higher than those predicted by the models. On the other hand, PDS70, a transition disk with two detected giant planets in the cavity, shows mass accretion rates well in line with model predictions.

ALMA surveys have shown that CO emission in protoplanetary disks is much fainter than expected. Accordingly, CO-based gas masses and gas/dust ratios are orders of magnitude lower than previously thought. This may be explained either as fast gas dispersal, or as chemical evolution and locking up of volatiles in larger bodies leading to the low observed CO fluxes. The latter processes lead to enhanced C/O ratios in the gas, which may be reflected in enhanced abundances of molecules like C2H. The goal of this work is to employ C2H observations to understand whether low CO fluxes are caused by volatile depletion, or by fast gas dissipation. We present ALMA Cycle 4 C2H observations of a subsample of nine sources in Lupus. The integrated C2H emission is determined and compared to previous CO observations and model predictions. Seven out of nine disks are detected in C2H, whose line emission is almost as bright as 13CO. All detections are significantly brighter than the typical sensitivity of the observations, hinting at a bimodal distribution of the C2H line intensities. When compared with DALI models, the observed C2H fluxes can be reproduced only if some level of volatile carbon and oxygen depletion is allowed and C/O>1 in the gas. Models with reduced gas/dust ratios fail instead to reproduce the observations. A steeper than linear correlation between C2H and CN emission line is found. This is linked to the fact that C2H emission lines are affected more strongly by C/O variations than CN lines. Ring-like structures are detected both in C2H and in continuum emission but, as for CN, they do not seem to be connected. Sz 71 shows ring shaped emission in both C2H and CN with the location of the peak intensity coinciding. Our new ALMA C2H observations favour volatile carbon and oxygen depletion rather than fast gas dispersal to explain the faint CO observations for most of the disks.

The mass accretion rate is a crucial parameter for the study of the evolution of accretion discs around young low-mass stellar and substellar objects (YSOs). We revisit the case of MYLup, an object where VLT/X-Shooter data suggested a negligible mass accretion rate, and show it to be accreting on a level similar to other ClassII YSOs in Lupus based on HST observations. In our HST-COS and -STIS spectra, we find many emission lines, as well as substantial FUV continuum excess emission, which can be ascribed to active accretion. The total luminosity of the CIV 1549A doublet is 4.1$\times10^{-4}$ Lsun. Using scalings between accretion luminosity, Lacc, and CIV luminosity from the literature, we derive Lacc $\sim$2$\times10^{-1}$ Lsun, which is more than an order of magnitude higher than the upper limit estimated from the X-Shooter observations. We discuss possible reasons for the X-Shooter-HST discrepancy, the most plausible being that the low contrast between the continuum excess emission and the photospheric+chromospheric emission at optical wavelengths in MYLup hampered detection of excess emission. The luminosity of the FUV continuum and CIV lines, strong H$_{2}$ fluorescence, and a "1600 A Bump" place MYLup in the class of accreting objects with gas-rich discs. So far, MYLup is the only peculiar case in which a significant difference between the HST and X-Shooter Macc estimates exists that is not ascribable to variability. The mass accretion rate inferred from the revisited Lacc estimate is Macc $\sim$ 1( $^{+1.5}_{-0.5}$ )$\times10^{-8}$ Msun/yr. This value is consistent with the typical value derived for accreting YSOs of similar mass in Lupus and points to less clearing of the inner disc than indicated by near- and mid-infrared observations. This is confirmed by ALMA data, which show that the gaps and rings seen in the sub-millimetre are relatively shallow.

The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.

A good constraint of when the growth of dust grains from sub-micrometer to millimeter sizes occurs, is crucial for planet formation models. This provides the first step towards the production of pebbles and planetesimals in protoplanetary disks. Currently, it is well established that Class II objects have large dust grains. However, it is not clear when in the star formation process this grain growth occurs. We use multi-wavelength millimeter observations of a Class I protostar to obtain the spectral index of the observed flux densities $\alpha_\mathrm{mm}$ of the unresolved disk and the surrounding envelope. Our goal is to compare our observational results with visibility modeling at both wavelengths simultaneously. We present data from NOEMA at 2.7 mm and SMA at 1.3 mm of the Class I protostar, Per-emb-50. We model the dust emission with a variety of parametric and radiative transfer models to deduce the grain size from the observed emission spectral index. We find a spectral index in the envelope of Per-emb-50 of $\alpha_{\rm env}$=$3.3\pm0.3$, similar to the typical ISM values. The radiative transfer modeling of the source confirms this value of $\alpha_{\rm env}$ with the presence of dust with a $a_\mathrm{max}$$\leq$100 $\mu$m. Additionally, we explore the backwarming effect, where we find that the envelope structure affects the millimeter emission of the disk. Our results reveal grains with a maximum size no larger than $100$ $\mu$m in the inner envelope of the Class I protostar Per-emb-50, providing an interesting case to test the universality of millimeter grain growth expected in these sources.

The formation of planets occurs within protoplanetary disks surrounding young stars, resulting in perturbation of the gas and dust surface densities. Here, we report the first evidence of spatially resolved gas surface density ($\Sigma_{g}$) perturbation towards the AS~209 protoplanetary disk from the optically thin C$^{18}$O ($J=2-1$) emission. The observations were carried out at 1.3~mm with ALMA at a spatial resolution of about 0.3$\arcsec$ $\times$ 0.2$\arcsec$ (corresponding to $\sim$ 38 $\times$ 25 au). The C$^{18}$O emission shows a compact ($\le$60~au), centrally peaked emission and an outer ring peaking at 140~au, consistent with that observed in the continuum emission and, its azimuthally averaged radial intensity profile presents a deficit that is spatially coincident with the previously reported dust map. This deficit can only be reproduced with our physico-thermochemical disk model by lowering $\Sigma_{gas}$ by nearly an order of magnitude in the dust gaps. Another salient result is that contrary to C$^{18}$O, the DCO$^{+}$ ($J=3-2$) emission peaks between the two dust gaps. We infer that the best scenario to explain our observations (C$^{18}$O deficit and DCO$^{+}$ enhancement) is a gas perturbation due to forming-planet(s), that is commensurate with previous continuum observations of the source along with hydrodynamical simulations. Our findings confirm that the previously observed dust gaps are very likely due to perturbation of the gas surface density that is induced by a planet of at least 0.2~M$\rm_{Jupiter}$ in formation. Finally, our observations also show the potential of using CO isotopologues to probe the presence of saturn mass planet(s).

The cyanide radical CN is abundant in protoplanetary disks, with line fluxes often comparable to those of $^{13}$CO. It is known to be sensitive to UV irradiation of the upper disk atmosphere, with models predicting ring-shaped emission. We seek to characterize the CN emission from 94 Class II disks in the Lupus star forming region, compare it to observations in other regions, and interpret our observations with a grid of models. The CN emission morphology is discussed for two primordial disks, Sz 71 and Sz 98, and modeled in more detail. ALMA observed CN N=3-2 in Lupus disks down to sensitivities better than previous surveys. Models constructed with the physico-chemical code DALI are used to study the integrated fluxes of the disks and resolved emission of CN in disks without (dust) substructures. We find that CN N=3-2 is bright, and detected in 38% of sources, but its disk-integrated flux is not strongly correlated to either $^{13}$CO or continuum flux. Compared to pre-ALMA single-dish surveys, no significant difference in the CN flux distributions in Lupus and Taurus-Auriga is found, although $\rho$ Ophiuchus disks may be fainter on average. We find ring-shaped CN emission with peak radii of ~50 AU in two resolved disks. We conclude that a large fraction of sources is faint in CN; only exponential gas surface density cutoffs at $R_{c} \leq 15$ AU can reconcile observations with models. This is the first observational evidence of such a compact gas disk population in Lupus. Absolute intensities and the emission morphology of CN are reproduced by DALI models without the need for any continuum substructure; they are unrelated to the CO snowline location. These observations and the successful modeling of these rings provide a new probe of the structure and conditions in disks, and particularly their incident UV radiation field, if disk size is determined by the data.

In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $\alpha_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $\nu$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $\beta$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $\beta =1.7$ gives better results than $\beta=1$. A wide distribution of disk masses is also needed to produce some cases with $\alpha <2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account.

It is key to constrain the gas surface density distribution, Sigma_gas, as function of disk radius in protoplanetary disks. In this work we investigate if spatially resolved observations of rarer CO isotopologues may be good tracers of Sigma_gas. Physical-chemical models with different input Sigma_gas(R) are run. The input disk surface density profiles are compared with the simulated 13CO intensity radial profiles to check if and where the two follow each other. There is always an intermediate region in the disk where the slope of the 13CO radial emission profile and Sigma_gas(R) coincide. At small radii the line radial profile underestimates Sigma_gas, as 13CO emission becomes optically thick. The same happens at large radii where the column densities become too low and 13CO is not able to efficiently self-shield. If the gas surface density profile is a simple power-law of the radius, the input power-law index can be retrieved within 20% uncertainty if one choses the proper radial range. If instead Sigma_gas(R) follows the self-similar solution for a viscously evolving disk, retrieving the input power-law index becomes challenging, in particular for small disks. Nevertheless, it is found that the power-law index can be in any case reliably fitted at a given line intensity contour around 6 K km/s, and this produces a practical method to constrain the slope of Sigma_gas(R). Application of such a method is shown in the case study of the TW Hya disk. Spatially resolved 13CO line radial profiles are promising to probe the disk surface density distribution, as they directly trace Sigma_gas(R)profile at radii well resolvable by ALMA. There, chemical processes like freeze-out and isotope selective photodissociation do not affect the emission, and, assuming that the volatile carbon does not change with radius, no chemical model is needed when interpreting the observations.

A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $\alpha_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed.

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.

In recent ALMA surveys, the gas distributions and velocity structures of most of the protoplanetary disks can still not be imaged at high S/N due to the short integration time. In this work, we re-analyzed the ALMA 13CO (3-2) and C18O (3-2) data of 88 young stellar objects in Lupus with the velocity-aligned stacking method to enhance S/N and to study the kinematics and disk properties traced by molecular lines. This method aligns spectra at different positions in a disk based on the projected Keplerian velocities at their positions and then stacks them. This method enhances the S/N ratios of molecular-line data and allows us to obtain better detections and to constrain dynamical stellar masses and disk orientations. We obtain 13CO detections in 41 disks and C18O detections in 18 disks with 11 new detections in 13CO and 9 new detections in C18O after applying the method. We estimate the disk orientations and the dynamical stellar masses from the 13CO data. Our estimated dynamical stellar masses correlate with the spectroscopic stellar masses, and in a subsample of 16 sources, where the inclination angles are better constrained, the two masses are in a good agreement within the uncertainties and with a mean difference of 0.15 Msun. With more detections of fainter disks, our results show that high gas masses derived from the 13CO and C18O lines tend to be associated with high dust masses estimated from the continuum emission. Nevertheless, the scatter is large (0.9 dex), implying large uncertainties in deriving the disk gas mass from the line fluxes. We find that with such large uncertainties it is expected that there is no correlation between the disk gas mass and the mass accretion rate with the current data. Deeper observations to detect disks with gas masses <1E-5 Msun in molecular lines are needed to investigate the correlation between the disk gas mass and the mass accretion rate.

This work aims to understand which midplane conditions are probed by the DCO$^+$ emission in the disk around the Herbig Ae star HD 169142. We explore the sensitivity of the DCO$^+$ formation pathways to the gas temperature and the CO abundance. The DCO$^+$ $J$=3-2 transition was observed with ALMA at a spatial resolution of 0.3". The HD 169142 DCO$^+$ radial intensity profile reveals a warm, inner component at radii <30 AU and a broad, ring-like structure from ~50-230 AU with a peak at 100 AU just beyond the millimeter grain edge. We modeled DCO$^+$ emission in HD 169142 with a physical disk structure adapted from the literature, and employed a simple deuterium chemical network to investigate the formation of DCO$^+$ through the cold deuterium fractionation pathway via H$_2$D$^+$. Contributions from the warm deuterium fractionation pathway via CH$_2$D$^+$ are approximated using a constant abundance in the intermediate disk layers. Parameterized models show that alterations to the midplane gas temperature and CO abundance of the literature model are both needed to recover the observed DCO$^+$ radial intensity profile. The best-fit model contains a shadowed, cold midplane in the region z/r < 0.1 with an 8 K decrease in gas temperature and a factor of five CO depletion just beyond the millimeter grain edge, and a 2 K decrease in gas temperature for r > 120 AU. The warm deuterium fractionation pathway is implemented as a constant DCO$^+$ abundance of 2.0$\times$10$^{-12}$ between 30-70 K. The DCO$^+$ emission probes a reservoir of cold material in the HD 169142 outer disk that is not revealed by the millimeter continuum, the SED, nor the emission from the 12CO, 13CO, or C18O $J$=2-1 lines.

We present a study of molecular gas in the inner disk $\left(r < 20 \, \text{AU} \right)$ around RY Lupi, with spectra from HST-COS, HST-STIS, and VLT-CRIRES. We model the radial distribution of flux from hot gas in a surface layer between $r = 0.1-10$ AU, as traced by Ly$\alpha$-pumped H$_2$. The result shows H$_2$ emission originating in a ring centered at $\sim$3 AU that declines within $r < 0.1$ AU, which is consistent with the behavior of disks with dust cavities. An analysis of the H$_2$ line shapes shows that a two-component Gaussian profile $\left(\text{FWHM}_{broad, H_2} = 105 \pm 15 \, \text{km s}^{-1}; \, \text{FWHM}_{narrow, H_2} = 43 \pm 13 \, \text{km s}^{-1} \right)$ is statistically preferred to a single-component Gaussian. We interpret this as tentative evidence for gas emitting from radially separated disk regions $\left( \left \langle r_{broad, H_2} \right \rangle \sim 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, H_2} \right \rangle \sim 3 \pm 2 \, \text{AU} \right)$. The 4.7 $\mu$m $^{12}$CO emission lines are also well fit by two-component profiles $\left( \left \langle r_{broad, CO} \right \rangle = 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, CO} \right \rangle = 15 \pm 2 \, \text{AU} \right)$. We combine these results with 10 $\mu$m observations to form a picture of gapped structure within the mm-imaged dust cavity, providing the first such overview of the inner regions of a young disk. The HST SED of RY Lupi is available online for use in modeling efforts.

Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities hinting at recently formed, giant planets. However, many of these studies are biased towards the brightest disks in the nearby star forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (>20 AU radius) from a complete disk survey of the Lupus star forming region, using ALMA Band 7 observations at 0.3" (22-30 AU radius) resolution of the 345 GHz continuum, 13CO and C18O 3-2 observations and the Spectral Energy Distribution of each source. Gas and dust surface density profiles are derived using the physical-chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star forming region. The dust cavity sizes range from 20-90 AU radius and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is ~11%, which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities.

The formation of planets strongly depends on the total amount as well as on the spatial distribution of solids in protoplanetary disks. Thanks to the improvements in resolution and sensitivity provided by ALMA, measurements of the surface density of mm-sized grains are now possible on large samples of disks. Such measurements provide statistical constraints that can be used to inform our understanding of the initial conditions of planet formation. We analyze spatially resolved observations of 36 protoplanetary disks in the Lupus star forming complex from our ALMA survey at 890 micron, aiming to determine physical properties such as the dust surface density, the disk mass and size and to provide a constraint on the temperature profile. We fit the observations directly in the uv-plane using a two-layer disk model that computes the 890 micron emission by solving the energy balance at each disk radius. For 22 out of 36 protoplanetary disks we derive robust estimates of their physical properties. The sample covers stellar masses between ~0.1 and ~2 Solar masses, and we find no trend between the average disk temperatures and the stellar parameters. We find, instead, a correlation between the integrated sub-mm flux (a proxy for the disk mass) and the exponential cut-off radii (a proxy of the disk size) of the Lupus disks. Comparing these results with observations at similar angular resolution of Taurus-Auriga/Ophiuchus disks found in literature and scaling them to the same distance, we observe that the Lupus disks are generally fainter and larger at a high level of statistical significance. Considering the 1-2 Myr age difference between these regions, it is possible to tentatively explain the offset in the disk mass/disk size relation with viscous spreading, however with the current measurements other mechanisms cannot be ruled out.

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 $>3\sigma$ 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.

Protoplanetary disks around young stars are the sites of planet formation. While the dust mass can be estimated using standard methods, determining the gas mass - and thus the amount of material available to form giant planets - has proven to be very difficult. Hydrogen deuteride (HD) is a promising alternative to the commonly-used gas mass tracer, CO. We aim to examine the robustness of HD as tracer of the disk gas mass, specifically the effect of gas mass on the HD FIR emission and its sensitivity to the vertical structure. Deuterium chemistry reactions relevant for HD were implemented in the thermochemical code DALI and models were run for a range of disk masses and vertical structures. The HD J=1-0 line intensity depends directly on the gas mass through a sublinear power law relation with a slope of ~0.8. Assuming no prior knowledge about the vertical structure of a disk and using only the HD 1-0 flux, gas masses can be estimated to within a factor of 2 for low mass disks (M$_{\rm disk} < 10^{-3}$ M$_\odot$). For more massive disks, this uncertainty increases to more than an order of magnitude. Adding the HD 2-1 line or independent information about the vertical structure can reduce this uncertainty to a factor of ~3 for all disk masses. For TW Hya, using the radial and vertical structure from Kama et al. 2016b the observations constrain the gas mass to $6\cdot10^{-3}$ M$_\odot$ < M$_{\rm disk} < 9\cdot10^{-3}$ M$_\odot$. Future observations require a 5$\sigma$ sensitivity of $1.8\cdot10^{-20}$ W m$^{-2}$ ($2.5\cdot10^{-20}$ W m$^{-2}$) and a spectral resolving power R > 300 (1000) to detect HD 1-0 (HD 2-1) for all disk masses above $10^{-5}$ M$_\odot$ with a line-to-continuum ratio > 0.01. These results show that HD can be used as an independent gas mass tracer with a relatively low uncertainty and should be considered as an important science goal for future FIR missions.

The $\sigma$ Orionis cluster is important for studying protoplanetary disk evolution, as its intermediate age ($\sim$3-5 Myr) is comparable to the median disk lifetime. We use ALMA to conduct a high-sensitivity survey of dust and gas in 92 protoplanetary disks around $\sigma$ Orionis members with $M_{\ast}\gtrsim0.1 M_{\odot}$. Our observations cover the 1.33 mm continuum and several CO $J=2-1$ lines: out of 92 sources, we detect 37 in the mm continuum and six in $^{12}$CO, three in $^{13}$CO, and none in C$^{18}$O. Using the continuum emission to estimate dust mass, we find only 11 disks with $M_{\rm dust}\gtrsim10 M_{\oplus}$, indicating that after only a few Myr of evolution most disks lack sufficient dust to form giant planet cores. Stacking the individually undetected continuum sources limits their average dust mass to 5$\times$ lower than that of the faintest detected disk, supporting theoretical models that indicate rapid dissipation once disk clearing begins. Comparing the protoplanetary disk population in $\sigma$ Orionis to those of other star-forming regions supports the steady decline in average dust mass and the steepening of the $M_{\rm dust}$-$M_{\ast}$ relation with age; studying these evolutionary trends can inform the relative importance of different disk processes during key eras of planet formation. External photoevaporation from the central O9 star is influencing disk evolution throughout the region: dust masses clearly decline with decreasing separation from the photoionizing source, and the handful of CO detections exist at projected separations $>1.5$ pc. Collectively, our findings indicate that giant planet formation is inherently rare and/or well underway by a few Myr of age.

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. To study the interaction between the protoplanet and the disk itself observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3\,mm, $^{12}$CO, $^{13}$CO and C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 56 - 83\,$au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity ($R \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 180\,$au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution).

An era has started in which gas and dust can be observed independently in protoplanetary disks, thanks to the recent surveys with ALMA. The first near-complete high-resolution disk survey in both dust and gas in a single star-forming region has been carried out in Lupus, finding surprisingly low gas/dust ratios. The goal of this work is to fully exploit CO isotopologues observations in Lupus, comparing them with physical-chemical model results, in order to obtain gas masses for a large number of disks. We have employed physical-chemical models to analyze continuum and CO isotopologues observations of Lupus disks, including isotope-selective processes and freeze-out. Employing also the ALMA 13CO-only detections, disk gas masses have been calculated for a total of 34 sources, expanding the sample of 10 disks studied by Ansdell et al. (2016), where also C18O was detected. We confirm that overall gas-masses are very low, often smaller than 1 $M_{\rm J}$, if volatile carbon is not depleted. Accordingly, global gas/dust ratios predominantly between 1 and 10. Low CO-based gas masses and gas/dust ratios may indicate rapid loss of gas, or alternatively chemical evolution, e.g. via sequestering of carbon from CO to more complex molecules, or carbon locked up in larger bodies. Current ALMA observations cannot distinguish between these two hypotheses. We have simulated both scenarios, but chemical model results do not allow us to rule out one of the two. Assuming that all Lupus disks have evolved mainly due to viscous processes over the past few Myr, the observed correlation between the current mass accretion rate and dust mass found by Manara et al. (2016) implies a constant gas-to-dust ratio, which is close to 100 based on the observed $M_{\rm disk}/\dot{M}_{\rm acc}$ ratio. This in turn points to a scenario in which carbon depletion is responsible for the low CO isotopologue line luminosities.

Despite intensive studies of protoplanetary disks, there is still no reliable way to determine their total mass and their surface density distribution, quantities that are crucial for describing both the structure and the evolution of disks up to the formation of planets. The goal of this work is to use less abundant CO isotopologues, whose detection is routine for ALMA, to infer the gas mass of disks. Isotope-selective effects need to be taken into account in the analysis, because they can significantly modify CO isotopologues line intensities. CO isotope-selective photodissociation has been implemented in the physical-chemical code DALI and 800 disk models have been run for a range of disk and stellar parameters. Dust and gas temperature structures have been computed self-consistently, together with a chemical calculation of the main species. Both disk structure and stellar parameters have been investigated. Total fluxes have been ray-traced for different CO isotopologues and for various transitions for different inclinations. A combination of 13CO and C18O total intensities allows inference of the total disk mass, although with non-negligible uncertainties. These can be overcome by employing spatially resolved observations, i.e. the disk's radial extent and inclination. Comparison with parametric models shows differences at the factor of a few level, especially for extremely low and high disk masses. Finally, total line intensities for different CO isotopologue and for various low-J transitions are provided and are fitted to simple formulae. The effects of a lower gas-phase carbon abundance and different gas/dust ratios are investigated as well, and comparison with other tracers is made. Disk masses can be determined within a factor of a few by comparing CO isotopologue lines observations with the simulated line fluxes, modulo the uncertainties in the volatile elemental abundances.

A relation between the mass accretion rate onto the central young star and the mass of the surrounding protoplanetary disk has long been theoretically predicted and observationally sought. For the first time, we have accurately and homogeneously determined the photospheric parameters, mass accretion rate, and disk mass for an essentially complete sample of young stars with disks in the Lupus clouds. Our work combines the results of surveys conducted with VLT/X-Shooter and ALMA. With this dataset we are able to test a basic prediction of viscous accretion theory, the existence of a linear relation between the mass accretion rate onto the central star and the total disk mass. We find a correlation between the mass accretion rate and the disk dust mass, with a ratio that is roughly consistent with the expected viscous timescale when assuming an interstellar medium (ISM) gas-to-dust ratio. This confirms that mass accretion rates are related to the properties of the outer disk. We find no correlation between mass accretion rates and the disk mass measured by CO isotopologues emission lines, possibly owing to the small number of measured disk gas masses. This suggests that the mm-sized dust mass better traces the total disk mass and that masses derived from CO may be underestimated, at least in some cases.

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $\leq0.1\,M_{\odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{\rm gas}=(0.1-1.5)\times 10^{-4}$. HD 100546 itself is a $\lambda\,$Bo\"{o}tis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{\rm gas}=(0.2-5.0)\times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{\rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems.

Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

One of the methods for deriving disk masses relies on direct observations of the gas, whose bulk mass is in the outer cold ($T\lesssim30$K) regions. This zone can be well traced by rotational lines of less abundant CO isotopologues, that probe the gas down to the midplane. The total CO gas mass is then obtained with the isotopologue ratios taken to be constant at the elemental isotope values found in the local ISM. This approach is however imprecise, because isotope selective processes are ignored. The aim of this work is an isotopologue selective treatment of CO isotopologues, in order to obtain a more accurate determination of disk masses. The isotope-selective photodissociation, the main process controlling the abundances of CO isotopologues in the CO-emissive layer, is properly treated for the first time in a full disk model (DALI, Bruderer et al. 2012; Bruderer 2013). The chemistry, thermal balance, line and continuum radiative transfer are all considered together with a chemical network that treats $^{13}$CO, C$^{18}$O, C$^{17}$O, isotopes of all included atoms, and molecules, as independent species. Isotope selective processes lead to regions in the disk where the isotopologues abundance ratios are considerably different from the elemental ratio. The results of this work show that considering CO isotopologue ratios as constants can lead to an underestimate of disk masses by up to almost two orders of magnitude if grains have grown to larger sizes. This may explain observed discrepancies in mass determinations from different tracers. The dependence of the various isotopologues emission on stellar and disk parameters is investigated. Including CO isotope selective processes is crucial to determine the gas mass of the disk accurately (through ALMA observations) and thus to provide the amount of gas which may eventually form planets or change the dynamics of forming planetary systems.

We present new 3 mm ATCA data of two Class I Young Stellar Objects in the Ophiucus star forming region: Elias29 and WL12. For our analysis we compare them with archival 1.1 mm SMA data. In the (u,v) plane the two sources present a similar behavior: a nearly constant non-zero emission at long baselines, which suggests the presence of an unresolved component and an increase of the fluxes at short baselines, related to the presence of an extended envelope. Our data analysis leads to unusually low values of the spectral index $\alpha_{\rm 1.1-3mm}$, which may indicate that mm-sized dust grains have already formed both in the envelopes and in the disk-like structures at such early stages. To explore the possible scenarios for the interpretation of the sources we perform a radiative transfer modeling using a Monte Carlo code, in order to take into account possible deviations from the Rayleigh-Jeans and optically thin regimes. Comparison between the model outputs and the observations indicates that dust grains may form aggregates up to millimeter size already in the inner regions of the envelopes of Class I YSOs. Moreover, we conclude that the embedded disk-like structures in our two Class Is are probably very compact, in particular in the case of WL12, with outer radii down to tens of AU.

Deep HST broad-band images taken with ACS and WFPC2 of the giant ($\sim 1000$ AU diameter) dark silhouette proplyd 114-426 in the Orion Nebula show that this system is tilted, asymmetric, warped and photoevaporated. The exquisite angular resolution of ACS allows us to map the distribution of dust grains at the northern translucent edge of the disk, dominated by the photoevaporative flow. Using the Mie theory for standard circumstellar disk grains, we find evidence for a spatial gradient in grain size. The typical dust radius, $\simeq 0.2-0.7 \mu$m (less than what reported by previous studies) becomes smaller as the distance from the disk center increases, consistent with the expectations for the dynamic of dust entrained in a gaseous photoevaporative wind. Our analysis of the disk morphology and location within the nebula indicates that this system is photoevaporated by the diffuse radiation field of the Orion Nebula, while being shielded from the radiation coming directly from the central Trapezium stars. We estimate the mass-loss rate from the disk surface and the time-scale for total disk dissipation, which turns out to be of the order of $10^4$yr. Such a short time, of the order of 1/100 of the cluster age, indicates that this system is seen on the verge of destruction. This is compatible with the exceptional nature of the disk, namely its combination of huge size and low mass. Finally, we briefly discuss the viability of possible mechanisms that may lead to the peculiar morphology of this system: external UV flux, binary star and past close encounter.