Molecules with ALMA at Planet-forming Scales (MAPS). VI. Distribution of the Small Organics HCN, C2H, and H2CO

Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high-angular-resolution (10–50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find that the HCN emission arises in a temperate (20–30 K) layer in the disk, while C2H is present in relatively warmer (20–60 K) layers. In the case of HD 163296, we find a decrease in column density for HCN and C2H inside one of the dust gaps near ∼83 au, where a planet has been proposed to be located. We derive H2CO column density profiles assuming temperatures between 20 and 50 K, and find slightly higher column densities in the colder disks around T Tauri stars than around Herbig Ae stars. The H2CO column densities rise near the location of the CO snowline and/or millimeter dust edge, suggesting an efficient release of H2CO ices in the outer disk. Finally, we find that the inner 50 au of these disks are rich in organic species, with abundances relative to water that are similar to cometary values. Comets could therefore deliver water and key organics to future planets in these disks, similar to what might have happened here on Earth. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Exploring HNC and HCN line emission as probes of the protoplanetary disk temperature

The distributions and abundances of molecules in protoplanetary disks are powerful tracers of the physical and chemical disk structures. The abundance ratios of HCN and its isomer HNC are known to be sensitive to gas temperature. Their line ratios might therefore offer a unique opportunity to probe the properties of the emitting gas. Using the 2D thermochemical code DALI, we ran a set of disk models accounting for different stellar and disk properties, with an updated chemical network for the nitrogen chemistry. These modeling results were then compared with observations, including new ALMA observations of HNC $J=3-2$ for the TW Hya disk and HNC $J=1-0$ for 29 disks in Lupus. Similar to CN, HCN and HNC have brighter line emission in models with larger disk flaring angles and higher UV fluxes. HNC and HCN are predicted to be abundant in the warm surface layer and outer midplane region, which results in ring-shaped emission patterns. However, the precise emitting regions and emission morphology depend on the probed transition, as well as on other parameters such as C and O abundances. The modeled HNC-to-HCN line intensity ratio increases from $<0.1$ in the inner disk to up to 0.8 in the outer disk regions, which can be explained by efficient HNC destruction at high temperatures. Disk-integrated HNC line fluxes from current scarce observations and its radial distribution in the TW Hya disk are broadly consistent with our model predictions. The HNC-to-HCN flux ratio robustly increases with radius (decreasing temperature), but its use as a chemical thermometer in disks is affected by other factors, including UV flux and C and O abundances. High-spatial resolution ALMA disk observations of HNC and HCN that can locate the emitting layers would have the great potential to constrain both the disk thermal and UV radiation structures, and also to verify our understanding of the nitrogen chemistry

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. I. Radial and Vertical Distributions of DCN and DCO+

Molecular D/H ratios are frequently used to probe the chemical past of solar system volatiles. Yet it is unclear which parts of the solar nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0farcs2–0farcs4 Atacama Large Millimeter/submillimeter Array (ALMA) observations of DCO+ and DCN 2–1, 3–2, and 4–3 toward the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project, augmented with archival data. DCO+ is characterized by an excitation temperature of ~40 K across the 70 au radius pebble disk, indicative of emission from a warm, elevated molecular layer. Tentatively, DCN is present at even higher temperatures. Both DCO+ and DCN present substantial emission cavities in the inner disk, while in the outer disk the DCO+ and DCN morphologies diverge: most DCN emission originates from a narrow ring peaking around 30 au, with some additional diffuse DCN emission present at larger radii, while DCO+ is present in a broad structured ring that extends past the pebble disk. Based on a set of simple parametric disk abundance models, these emission patterns can be explained by a near-constant DCN abundance exterior to the cavity, and an increasing DCO+ abundance with radius. In conclusion, the ALMA observations reveal an active deuterium fractionation chemistry in multiple disk regions around TW Hya, but not in the cold planetesimal-forming midplane and in the inner disk. More observations are needed to explore whether deuterium fractionation is actually absent in these latter regions, and if its absence is a common feature or something peculiar to the old TW Hya disk

Molecules with ALMA at Planet-forming Scales (MAPS). XI. CN and HCN as Tracers of Photochemistry in Disks

UV photochemistry in the surface layers of protoplanetary disks dramatically alters their composition relative to previous stages of star formation. The abundance ratio CN/HCN has long been proposed to trace the UV field in various astrophysical objects; however, to date the relationship between CN, HCN, and the UV field in disks remains ambiguous. As part of the ALMA Large Program MAPS (Molecules with ALMA at Planet-forming Scales), we present observations of CN N = 1–0 transitions at 0farcs3 resolution toward five disk systems. All disks show bright CN emission within ∼50–150 au, along with a diffuse emission shelf extending up to 600 au. In all sources we find that the CN/HCN column density ratio increases with disk radius from about unity to 100, likely tracing increased UV penetration that enhances selective HCN photodissociation in the outer disk. Additionally, multiple millimeter dust gaps and rings coincide with peaks and troughs, respectively, in the CN/HCN ratio, implying that some millimeter substructures are accompanied by changes to the UV penetration in more elevated disk layers. That the CN/HCN ratio is generally high (>1) points to a robust photochemistry shaping disk chemical compositions and also means that CN is the dominant carrier of the prebiotically interesting nitrile group at most disk radii. We also find that the local column densities of CN and HCN are positively correlated despite emitting from vertically stratified disk regions, indicating that different disk layers are chemically linked. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS). XII. Inferring the C/O and S/H Ratios in Protoplanetary Disks with Sulfur Molecules

Sulfur-bearing molecules play an important role in prebiotic chemistry and planet habitability. They are also proposed probes of chemical ages, elemental C/O ratio, and grain chemistry processing. Commonly detected in diverse astrophysical objects, including the solar system, their distribution and chemistry remain, however, largely unknown in planet-forming disks. We present CS (2 − 1) observations at ∼0farcs3 resolution performed within the ALMA MAPS Large Program toward the five disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. CS is detected in all five disks, displaying a variety of radial intensity profiles and spatial distributions across the sample, including intriguing apparent azimuthal asymmetries. Transitions of C2S and SO were also serendipitously covered, but only upper limits are found. For MWC 480, we present complementary ALMA observations at ∼ 0farcs5 of CS, 13CS, C34S, H2CS, OCS, and SO2. We find a column density ratio N(H2CS)/N(CS) ∼ 2/3, suggesting that a substantial part of the sulfur reservoir in disks is in organic form (i.e., CxHySz). Using astrochemical disk modeling tuned to MWC 480, we demonstrate that N(CS)/N(SO) is a promising probe for the elemental C/O ratio. The comparison with the observations provides a supersolar C/O. We also find a depleted gas-phase S/H ratio, suggesting either that part of the sulfur reservoir is locked in solid phase or that it remains in an unidentified gas-phase reservoir. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS). XVIII. Kinematic Substructures in the Disks of HD 163296 and MWC 480

We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet-forming Scales (MAPS) large program. Using the J = 2–1 transitions of 12CO, 13CO, and C18O imaged at spatial resolutions of ∼0farcs15 and with a channel spacing of 200 m s−1, we find perturbations from Keplerian rotation in the projected velocity fields of both disks (≲5% of the local Keplerian velocity), suggestive of large-scale (tens of astronomical units in size), coherent flows. By accounting for the azimuthal dependence on the projection of the velocity field, the velocity fields were decomposed into azimuthally averaged orthogonal components, vϕ, vr, and vz. Using the optically thick 12CO emission as a probe of the gas temperature, local variations of ≈3 K (≈5% relative changes) were observed and found to be associated with the kinematic substructures. The MWC 480 disk hosts a suite of tightly wound spiral arms. The spirals arms, in conjunction with the highly localized perturbations in the gas velocity structure (kinematic planetary signatures), indicate a giant planet, ∼1 MJup, at a radius of ≈245 au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. and Teague et al. advocating for multiple ∼1 MJup planets embedded in the disk. These results demonstrate that molecular line observations that characterize the dynamical structure of disks can be used to search for the signatures of embedded planets. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement