The TW Hya Rosetta Stone Project IV: A Hydrocarbon-rich Disk Atmosphere

Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like C2H and C3H2 have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved c–C3H2 observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at 15–20 au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, Eu = 29–97 K. We find the ortho-to-para ratio of c–C3H2 is consistent with 3 throughout extent of the emission, and the total abundance of both c–C3H2 isomers is (7.5–10) × 10−11 per H atom, or 1%–10% of the previously published C2H abundance in the same source. We find c–C3H2 comes from a layer near the surface that extends no deeper than z/r = 0.25. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside ~30 au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridional flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry

The TW Hya Rosetta Stone Project. II. Spatially Resolved Emission of Formaldehyde Hints at Low-temperature Gas-phase Formation

Formaldehyde (H2CO) is an important precursor to organics like methanol (CH3OH). It is important to understand the conditions that produce H2CO and prebiotic molecules during star and planet formation. H2CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold (lesssim20 K) dust grains. To understand which pathway dominates, gaseous H2CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and <3 linked to cold formation in the solid state. We present spatially resolved Atacama Large Millimeter/submillimeter Array observations of multiple ortho- and para-H2CO transitions in the TW Hya protoplanetary disk to test H2CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of 33 ± 2 K, (1.1 ± 0.1) × 1012 cm−2 and 25 ± 2 K, (4.4 ± 0.3) × 1011 cm−2 for ortho- and para-H2CO, respectively, and an OPR of 2.49 ± 0.23. A radially resolved analysis shows that the observed H2CO emits mostly at rotational temperatures of 30–40 K, corresponding to a layer with z/R ≥ 0.25. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to 2.0 ± 0.5. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, led us to speculate that gas-phase formation is responsible for the observed H2CO across the TW Hya disk

UV-driven chemistry as a signpost of late-stage planet formation

The chemical reservoir within protoplanetary disks has a direct impact on planetary compositions and the potential for life. A long-lived carbon- and nitrogen-rich chemistry at cold temperatures (≤ 50 K) is observed within cold and evolved planet-forming disks. This is evidenced by bright emission from small organic radicals in 1–10 Myr aged systems that would otherwise have frozen out onto grains within 1 Myr. We explain how the chemistry of a planet-forming disk evolves from a cosmic-ray/X-ray-dominated regime to a ultraviolet-dominated chemical equilibrium. This, in turn, will bring about a temporal transition in the chemical reservoir from which planets will accrete. This photochemical dominated gas phase chemistry develops as dust evolves via growth, settling and drift, and the small grain population is depleted from the disk atmosphere. A higher gas-to-dust mass ratio allows for deeper penetration of ultraviolet photons is coupled with a carbon-rich gas (C/O > 1) to form carbon-bearing radicals and ions. This further results in gas phase formation of organic molecules, which then would be accreted by any actively forming planets present in the evolved disk

The TW Hya Rosetta Stone Project IV: A Hydrocarbon-rich Disk Atmosphere

Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like C2H and C3H2 have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved c–C3H2 observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at 15–20 au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, Eu = 29–97 K. We find the ortho-to-para ratio of c–C3H2 is consistent with 3 throughout extent of the emission, and the total abundance of both c–C3H2 isomers is (7.5–10) × 10^−11 per H atom, or 1%–10% of the previously published C2H abundance in the same source. We find c–C3H2 comes from a layer near the surface that extends no deeper than z/r = 0.25. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside ~30 au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridional flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry

The TW Hya Rosetta Stone Project. II. Spatially Resolved Emission of Formaldehyde Hints at Low-temperature Gas-phase Formation

Formaldehyde (H2CO) is an important precursor to organics like methanol (CH3OH). It is important to understand the conditions that produce H2CO and prebiotic molecules during star and planet formation. H2CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold (≲20 K) dust grains. To understand which pathway dominates, gaseous H2CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and <3 linked to cold formation in the solid state. We present spatially resolved Atacama Large Millimeter/submillimeter Array observations of multiple ortho- and para-H2CO transitions in the TW Hya protoplanetary disk to test H2CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of 33 ± 2 K, (1.1 ± 0.1) × 1012 cm−2 and 25 ± 2 K, (4.4 ± 0.3) × 1011 cm−2 for ortho- and para-H2CO, respectively, and an OPR of 2.49 ± 0.23. A radially resolved analysis shows that the observed H2CO emits mostly at rotational temperatures of 30–40 K, corresponding to a layer with z/R ≥ 0.25. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to 2.0 ± 0.5. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, led us to speculate that gas-phase formation is responsible for the observed H2CO across the TW Hya disk

The TW Hya Rosetta Stone Project. III. Resolving the Gaseous Thermal Profile of the Disk

The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide-reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J = 2–1 and C18O J = 2–1 as well as archival ALMA observations of 12CO J = 3–2, 13CO J = 3–2 and 6–5, and C18O J = 3–2 and 6–5. Additionally, we reproduce a Herschel observation of the HD J = 1–0 line flux and the spectral energy distribution and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best-fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 M⊙ and a thin upper layer of gas depleted of small dust with a thickness of ~1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer, as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge, and, in this work, HD provides the additional constraint to derive the gas mass and support the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future

Molecules with ALMA at Planet-forming Scales (MAPS). VIII. CO Gap in AS 209—Gas Depletion or Chemical Processing?

Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial-resolution data from the Molecules with ALMA at Planet-forming Scales (MAPS) program. The observations of C18O J = 2–1 emission exhibit a strong depression at 88 au overlapping with the positions of multiple gaps in millimeter dust continuum emission. We find that the observed CO column density is consistent with either gas surface-density perturbations or chemical processing, while C2H column density traces changes in the C/O ratio rather than the H2 gas surface density. However, the presence of a massive planet (>0.2 MJup) would be required to account for this level of gas depression, which conflicts with constraints set by the dust emission and the pressure profile measured by gas kinematics. Based on our models, we infer that a local decrease of CO abundance is required to explain the observed structure in CO, dominating over a possible gap-carving planet present and its effect on the H2 surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales. XX. The Massive Disk around GM Aurigae

Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopologue HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to H2 and the disk temperature structure. Using new Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program as well as archival ALMA observations, we construct a disk physical/chemical model of the protoplanetary disk GM Aur. Our model is in good agreement with the spatially resolved CO isotopologue emission from 11 rotational transitions with spatial resolution ranging from 0farcs15 to 0farcs46 (24–73 au at 159 pc) and the spatially unresolved HD J = 1–0 detection from Herschel. Our best-fit model favors a cold protoplanetary disk with a total gas mass of approximately 0.2 M⊙, a factor of 10 reduction in CO gas inside roughly 100 au and a factor of 100 reduction outside of 100 au. Despite its large mass, the disk appears to be on the whole gravitationally stable based on the derived Toomre Q parameter. However, the region between 70 and 100 au, corresponding to one of the millimeter dust rings, is close to being unstable based on the calculated Toomre Q of <1.7. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS). XVI. Characterizing the Impact of the Molecular Wind on the Evolution of the HD 163296 System

During the main phase of evolution of a protoplanetary disk, accretion regulates the inner-disk properties, such as the temperature and mass distribution, and in turn, the physical conditions associated with planet formation. The driving mechanism behind accretion remains uncertain; however, one promising mechanism is the removal of a fraction of angular momentum via a magnetohydrodynamic (MHD) disk wind launched from the inner tens of astronomical units of the disk. This paper utilizes CO isotopologue emission to study the unique molecular outflow originating from the HD 163296 protoplanetary disk obtained with the Atacama Large Millimeter/submillimeter Array. HD 163296 is one of the most well-studied Class II disks and is proposed to host multiple gas-giant planets. We robustly detect the large-scale rotating outflow in the 12CO J = 2 − 1 and the 13CO J = 2 − 1 and J = 1 − 0 transitions. We constrain the kinematics, the excitation temperature of the molecular gas, and the mass-loss rate. The high ratio of the rates of ejection to accretion (5–50), together with the rotation signatures of the flow, provides solid evidence for an MHD disk wind. We find that the angular momentum removal by the wind is sufficient to drive accretion though the inner region of the disk; therefore, accretion driven by turbulent viscosity is not required to explain HD 163296's accretion. The low temperature of the molecular wind and its overall kinematics suggest that the MHD disk wind could be perturbed and shocked by the previously observed high-velocity atomic jet. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS). XVII. Determining the 2D Thermal Structure of the HD 163296 Disk

Understanding the temperature structure of protoplanetary disks is key to interpreting observations, predicting the physical and chemical evolution of the disk, and modeling planet formation processes. In this study, we constrain the two-dimensional thermal structure of the disk around the Herbig Ae star HD 163296. Using the thermochemical code RAC2D, we derive a thermal structure that reproduces spatially resolved Atacama Large Millimeter/submillimeter Array observations (∼0farcs12 (13 au)–0farcs25 (26 au)) of 12CO J = 2 − 1, 13CO J = 1 − 0, 2 − 1, C18O J = 1 − 0, 2 − 1, and C17O J = 1 − 0, the HD J = 1 − 0 flux upper limit, the spectral energy distribution (SED), and continuum morphology. The final model incorporates both a radial depletion of CO motivated by a timescale shorter than typical CO gas-phase chemistry (0.01 Myr) and an enhanced temperature near the surface layer of the the inner disk (z/r ≥ 0.21). This model agrees with the majority of the empirically derived temperatures and observed emitting surfaces derived from the J = 2 − 1 CO observations. We find an upper limit for the disk mass of 0.35 M⊙, using the upper limit of the HD J = 1 − 0 and J = 2 − 1 flux. With our final thermal structure, we explore the impact that gaps have on the temperature structure constrained by observations of the resolved gaps. Adding a large gap in the gas and small dust additionally increases gas temperature in the gap by only 5%–10%. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement