skip to main content


Title: New Constraints on Protoplanetary Disk Gas Masses in Lupus
Abstract

Gas mass is a fundamental quantity of protoplanetary disks that directly relates to their ability to form planets. Because we are unable to observe the bulk H2content of disks directly, we rely on indirect tracers to provide quantitative mass estimates. Current estimates for the gas masses of the observed disk population in the Lupus star-forming region are based on measurements of isotopologues of CO. However, without additional constraints, the degeneracy between H2mass and the elemental composition of the gas leads to large uncertainties in such estimates. Here, we explore the gas compositions of seven disks from the Lupus sample representing a range of CO-to-dust ratios. With Band 6 and 7 ALMA observations, we measure line emission for HCO+, HCN, and N2H+. We find a tentative correlation among the line fluxes for these three molecular species across the sample, but no correlation with13CO or submillimeter continuum fluxes. For the three disks where N2H+is detected, we find that a combination of high disk gas masses and subinterstellar C/H and O/H are needed to reproduce the observed values. We find increases of ∼10–100× previous mass estimates are required to match the observed line fluxes. This work highlights how multimolecular studies are essential for constraining the physical and chemical properties of the gas in populations of protoplanetary disks, and that CO isotopologues alone are not sufficient for determining the mass of many observed disks.

 
more » « less
Award ID(s):
1910106
NSF-PAR ID:
10364360
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
927
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 229
Size(s):
["Article No. 229"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instabilities or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successfully, we need to understand what the observed gas disk size actually traces. We studied the relation betweenRCO,90%, the radius that encloses 90% of the12CO flux, andRc, the radius that encodes the physical disk size, in order to provide simple prescriptions for conversions between these two sizes. For an extensive grid of thermochemical models, we calculateRCO,90%from synthetic observations and relate properties measured at this radius, such as the gas column density, to bulk disk properties, such asRcand the disk massMdisk. We found an empirical correlation between the gas column density atRCO,90%and disk mass:Ngas(RCO,90%)3.73×1021(Mdisk/M)0.34cm2. Using this correlation we derive an analytical prescription ofRCO,90%that only depends onRcandMdisk. We deriveRcfor disks in Lupus, Upper Sco, Taurus, and the DSHARP sample, finding that disks in the older Upper Sco region are significantly smaller (〈Rc〉 = 4.8 au) than disks in the younger Lupus and Taurus regions (〈Rc〉 = 19.8 and 20.9 au, respectively). This temporal decrease inRcgoes against predictions of both viscous and wind-driven evolution, but could be a sign of significant external photoevaporation truncating disks in Upper Sco.

     
    more » « less
  2. Abstract

    High-spatial-resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (≈30°–75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we can map the two-dimensional (r,z) disk structure from the disk upper atmosphere, as traced by CO, to near the midplane, as probed by less abundant isotopologues. Here, we present high-angular-resolution (≲0.″1 to ≈0.″2; ≈15–30 au) observations of CO,13CO, and C18O in either or bothJ= 2–1 andJ= 3–2 lines in the transition disks around DM Tau, Sz 91, LkCa 15, and HD 34282. We derived line emission surfaces in CO for all disks and in13CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C18O in any disk, which is instead consistent with C18O emission originating from the midplane. Both theJ= 2–1 andJ= 3–2 lines show similar heights. Using the derived emission surfaces, we computed radial and vertical gas temperature distributions for each disk, including empirical temperature models for the DM Tau and LkCa 15 disks. After combining our sample with literature sources, we find that13CO line emitting heights are also tentatively linked with source characteristics, e.g., stellar host mass, gas temperature, disk size, and show steeper trends than seen in CO emission surfaces.

     
    more » « less
  3. Abstract High spatial resolution CO observations of midinclination (≈30°–75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically thick, spatially resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of Atacama Large Millimeter/submillimeter Array archival data yields CO emission surfaces, dynamically constrained stellar host masses, and disk atmosphere gas temperatures for the disks around the following: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50–2.10 M ⊙ ), ages (∼0.3–23 Myr), and CO gas radial emission extents (≈200–1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ( z / r ≈ 0.1 to ≳0.5) hinted at in previous studies. We compute the radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ≈2–5× gas pressure scale heights (H g ) and could potentially be calibrated as empirical tracers of H g . 
    more » « less
  4. Abstract

    We present Keck Planet Imager and Characterizer (KPIC) high-resolution (R∼35,000)K-band thermal emission spectroscopy of the ultrahot Jupiter WASP-33b. The use of KPIC’s single-mode fibers greatly improves both blaze and line-spread stabilities relative to slit spectrographs, enhancing the cross-correlation detection strength. We retrieve the dayside emission spectrum with a nested-sampling pipeline, which fits for orbital parameters, the atmospheric pressure–temperature profile, and the molecular abundances. We strongly detect the thermally inverted dayside and measure mass-mixing ratios for CO (logCOMMR=1.10.6+0.4), H2O (logH2OMMR=4.10.9+0.7), and OH (logOHMMR=2.11.1+0.5), suggesting near-complete dayside photodissociation of H2O. The retrieved abundances suggest a carbon- and possibly metal-enriched atmosphere, with a gas-phase C/O ratio of0.80.2+0.1, consistent with the accretion of high-metallicity gas near the CO2snow line and post-disk migration or with accretion between the soot and H2O snow lines. We also find tentative evidence for12CO/13CO ∼ 50, consistent with values expected in protoplanetary disks, as well as tentative evidence for a metal-enriched atmosphere (2–15 × solar). These observations demonstrate KPIC’s ability to characterize close-in planets and the utility of KPIC’s improved instrumental stability for cross-correlation techniques.

     
    more » « less
  5. null (Ed.)
    Context. Recent years have seen building evidence that planet formation starts early, in the first ~0.5 Myr. Studying the dust masses available in young disks enables us to understand the origin of planetary systems given that mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. Aims. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. Methods. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (1.1–1.3 mm) continuum observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka -band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities measured in the image plane. Results. We find a strong linear correlation between the ALMA and VLA fluxes, demonstrating that emission at both wavelengths is dominated by dust emission. For a subsample of optically thin sources, we find a median spectral index of 2.5 from which we derive the dust opacity index β = 0.5, suggesting significant dust growth. Comparison with ALMA surveys of Orion shows that the Class I dust disk mass distribution between the two regions is similar, but that the Class 0 disks are more massive in Perseus than those in Orion. Using the DIANA opacity model including large grains, with a dust opacity value of κ 9 mm = 0.28 cm 2 g −1 , the median dust masses of the embedded disks in Perseus are 158 M ⊕ for Class 0 and 52 M ⊕ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M ⊕ and 12 M ⊕ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value κ 1.3 mm = 2.3 cm 2 g −1 . The dust masses of young Class 0 and I disks are larger by at least a factor of ten and three, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. Conclusions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of ~15%. A higher efficiency of ~30% is necessary if the planet formation is set to start in Class I disks. 
    more » « less