Search

Measuring the gas mass of protoplanetary disks, the reservoir available for giant planet formation, has proven to be difficult. We currently lack a far-infrared observatory capable of observing HD, and the most common gas mass tracer, CO, suffers from a poorly constrained CO-to-H$_2$ ratio. Expanding on previous work, we investigate if N2H+, a chemical tracer of CO poor gas, can be used to observationally measure the CO-to-H$_2$ ratio and correct CO-based gas masses. Using disk structures obtained from the literature, we set up thermochemical models for three disks, TW Hya, DM Tau and GM Aur, to examine how well the CO-to-H$_2$ ratio and gas mass can be measured from N2H+ and C18O line fluxes. Furthermore, we compare these gas masses to independently gas masses measured from archival HD observations. The N2H+ (3-2)/C18O (2-1) line ratio scales with the disk CO-to-H$_2$ ratio. Using these two lines, we measure 4.6e-3 Msun < Mdisk < 1.1e-1 Msun for TW Hya, 1.5e-2 Msun < Mdisk < 9.6e-2 Msun for GM Aur and 3.1e-2 Msun < Mdisk < 9.6e-2 Msun for DM Tau. These gas masses agree with values obtained from HD within their respective uncertainties. The uncertainty on the N2H+ + C18O gas mass can be reduced by observationally constraining the cosmic ray ionization rate in disks. These results demonstrate the potential of using the combination of N2H+ and C18O to measure gas masses of protoplanetary disks., Comment: 10 pages, 5 figures, accepted for publication in ApJL

[Abridged] Planet formation is expected to be enhanced around snowlines in protoplanetary disks, in particular around the water snowline. However, the close proximity of the water snowline to the host star and water in the Earth's atmosphere makes a direct detection of the water snowline in protoplanetary disks challenging. Following earlier work on protostellar envelopes, the aim of this research is to investigate the validity of HCO$^+$ and H$^{13}$CO$^+$, as tracers of the water snowline in protoplanetary disks, as HCO$^+$ is destroyed by gas-phase water. Two small chemical networks are used to predict the HCO$^+$ abundance in a typical Herbig Ae disk. Subsequently, the corresponding emission profiles are modelled for H$^{13}$CO$^+$ and HCO$^+$ $J=2-1$, which provides the best balance between brightness and optical depth effects of the continuum emission. The HCO$^+$ abundance jumps by two orders of magnitude just outside the water snowline at 4.5 AU. We find that the emission of H$^{13}$CO$^+$ and HCO$^+$ is ring-shaped due to three effects: destruction of HCO$^+$ by gas-phase water, continuum optical depth, and molecular excitation effects. The presence of gas-phase water causes an additional drop of only $\sim$13% and 24% in the center of the disk, for H$^{13}$CO$^+$ and HCO$^+$, respectively. For the much more luminous outbursting source V883Ori, our models predict that the effect of dust and excitation are not limiting if the snowline is located outside $\sim$40 AU. Our analysis of ALMA observations of HCO$^+$ $J=3-2$ is consistent with the water snowline located around 100 AU. The HCO$^+$ abundance drops steeply around the water snowline, but dust and excitation can conceal the drop in HCO$^+$ emission due to the water snowline. Therefore, locating the water snowline with HCO$^+$ in Herbig disks is very difficult, but it is possible for outbursting sources like V883Ori., Comment: Accepted for publication in Astronomy and Astrophysics. 13 pages, 11 figures + appendices

The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the disk gaseous outer radius are based on the extent of the CO rotational emission, which, during the evolution, is also affected by the changing physical and chemical conditions in the disk. We use physical-chemical DALI models to study how the extent of the CO emission changes with time in a viscously expanding disk and investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. We find that the gas outer radius (R_co) measured from our models matches the expectations of a viscously spreading disk: R_co increases with time and for a given time R_co is larger for a disk with a higher viscosity alpha_visc. However, in the extreme case where the disk mass is low (less than 10^-4 Msun) and alpha_visc is high (larger than 10^-2), R_co will instead decrease with time as a result of CO photodissociation in the outer disk. For most disk ages R_co is up to 12x larger than the characteristic size R_c of the disk, and R_co/R_c is largest for the most massive disk. As a result of this difference, a simple conversion of R_co to alpha_visc will overestimate the true alpha_visc of the disk by up to an order of magnitude. We find that most observed gas outer radii in Lupus can be explained using a viscously evolving disk that starts out small (R_c = 10 AU) and has a low viscosity (alpha_visc = 10^-4 - 10^-3)., Comment: 19 pages, 17 figures, accepted in A&A

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., Comment: 17 pages, 11 figures, accepted in A&A

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., Comment: Accepted for publication in A&A

The extent of the gas in protoplanetary disks is observed to be universally larger than the extent of the dust. This is often attributed to radial drift and grain growth of the mm grains, but line optical depth produces a similar observational signature. We investigate in what parts of the disk structure parameter space dust evolution and line optical depth are the dominant drivers of the observed gas and dust size difference. Using the thermochemical model DALI with dust evolution included we ran a grid of models aimed at reproducing the observed gas and dust size dichotomy. The relation between R_dust and dust evolution is non-monotonic and depends on the disk structure. R_gas is directly related to the radius where the CO column density drops below 10^15 cm^-2 and CO becomes photodissociated. R_gas is not affected by dust evolution but scales with the total CO content of the disk. R_gas/R_dust > 4 is a clear sign for dust evolution and radial drift in disks, but these cases are rare in current observations. For disks with a smaller R_gas/R_dust, identifying dust evolution from R_gas/R_dust requires modelling the disk structure including the total CO content. To minimize the uncertainties due to observational factors requires FWHM_beam < 1x the characteristic radius and a peak SNR > 10 on the 12CO emission moment zero map. For the dust outer radius to enclose most of the disk mass, it should be defined using a high fraction (90-95%) of the total flux. For the gas, any radius enclosing > 60% of the 12CO flux will contain most of the disk mass. To distinguish radial drift and grain growth from line optical depth effects based on size ratios requires disks to be observed at high enough angular resolution and the disk structure should to be modelled to account for the total CO content of the disk., Comment: 18 pages, 27 figures, accepted in A&A

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., Comment: 22 pages, 14 figures, A&A accepted

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., Comment: 12 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

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., Comment: 22 pages, 10 figures; accepted to ApJ