In this paper, we present preliminary results illustrating the effect of cosmic rays on solid-phase chemistry in models of both TMC-1 and several sources with physical conditions identical to TMC-1 ...except for hypothetically enhanced ionization rates. Using a recent theory for the addition of cosmic-ray-induced reactions to astrochemical models, we calculated the radiochemical yields, called G values, for the primary dust grain ice-mantle constituents. We show that the inclusion of this nonthermal chemistry can lead to the formation of complex organic molecules from simpler ice-mantle constituents, even under cold core conditions. In addition to enriching ice mantles, we find that these new radiation-chemical processes can lead to increased gas-phase abundances as well, particularly for HOCO, NO2, HC2O, methyl formate (HCOOCH3), and ethanol (CH3CH2OH). These model results imply that HOCO-and perhaps NO2-might be observable in TMC-1. Future detections of either of these two species in cold interstellar environments could provide strong support for the importance of cosmic-ray-driven radiation chemistry. The increased gas-phase abundance of methyl formate can be compared with abundances achieved through other formation mechanisms such as pure gas-phase chemistry and surface reactions.
Abstract
The Molecules with ALMA at Planet-forming Scales Large Program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high ...angular resolution (0.″15 and 0.″30 beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3–1.3 mJy beam
−1
and 0.4–1.9 mJy beam
−1
for Bands 6 and 3, respectively). In this article, we describe the multistage workflow—built around the CASA
tclean
image deconvolution procedure—that we used to generate the core data product of the MAPS LP: the position–position–velocity image cubes for each spectral line. Owing to the expansive nature of the survey, we encountered a range of imaging challenges: some are familiar to the submillimeter protoplanetary disk community, like the need to use an accurate CLEAN mask, and others are less well known, like the incorrect default flux scaling of the CLEAN residual map first described by Jorsater & van Moorsel (the “JvM effect”). We distill lessons learned into recommended workflows for synthesizing image cubes of molecular emission. In particular, we describe how to produce image cubes with accurate fluxes via “JvM correction,” a procedure that is generally applicable to any image synthesized via CLEAN deconvolution but is especially critical for low signal-to-noise ratio (S/N) emission. We further explain how we used visibility tapering to promote a common, fiducial beam size and contextualize the interpretation of S/N when detecting molecular emission from protoplanetary disks. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.
Abstract
Planets form and obtain their compositions in dust- and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The ...distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks—around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480—in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in four spectral setups, which together cover ∼50 lines from over 20 different species. This paper introduces the Astrophysical Journal Supplement’s MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical substructures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation.
The nature and abundance of sulfur chemistry in protoplanetary disks (PPDs) may impact the sulfur inventory on young planets and therefore their habitability. PPDs also offer an interesting test bed ...for sulfur chemistry models, since each disk shows a diverse set of environments. In this context, we present new sulfur molecule observations in PPDs and new S-disk chemistry models. With the Atacama Large Millimeter/submillimeter Array we observed the CS 5-4 rotational transition toward five PPDs (DM Tau, DO Tau, CI Tau, LkCa 15, MWC 480) and the CS 6-5 transition toward three PPDs (LkCa 15, MWC 480, and V4046 Sgr). Across this sample, CS displays a range of radial distributions, from centrally peaked to gaps and rings. We also present the first detection in PPDs of 13CS 6-5 (LkCa 15 and MWC 480), C34S 6-5 (LkCa 15), and H2CS 817-716, 919-818, and 918-817 (MWC 480) transitions. Using LTE models to constrain column densities and excitation temperatures, we find that either 13C and 34S are enhanced in CS or CS is optically thick despite its relatively low brightness temperature. Additional lines and higher spatial resolution observations are needed to distinguish between these scenarios. Assuming that CS is optically thin, CS column density model predictions reproduce the observations within a factor of a few for both MWC 480 and LkCa 15. However, the model underpredicts H2CS by 1-2 orders of magnitude. Finally, comparing the H2CS/CS ratio observed toward the MWC 480 disk and toward different interstellar medium sources, we find the closest match with prestellar cores.
Abstract
We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of ...1.″4 from the star at a position angle of 161°), isolated via
13
CO
J
= 2−1 emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a 117 × 82 mas beam and manifests as a point source in
13
CO, indicating that its diameter is ≲14 au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using
12
CO and near-infrared scattered-light observations and is associated with localized velocity perturbations in
12
CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the
13
CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is ≳35 K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is ≳0.095
M
Jup
≃ 30
M
⊕
adopting a standard
13
CO abundance. From the nondetection of millimeter continuum emission at the location of the CPD (3
σ
flux density ≲26.4
μ
Jy), we infer that the CPD dust mass is ≲0.027
M
⊕
≃ 2.2 lunar masses, indicating a low dust-to-gas mass ratio of ≲9 × 10
−4
. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion.
Shocks are a crucial probe for understanding the ongoing chemistry within ices on interstellar dust grains, where many complex organic molecules (COMs) are believed to be formed. However, previous ...work has been limited to the initial liberation into the gas phase through nonthermal desorption processes such as sputtering. Here we present results from the adapted three-phase gas-grain chemical network code nautilus, with the inclusion of additional high-temperature reactions, nonthermal desorption, collisional dust heating, and shock physics parameters. This enhanced model is capable of reproducing many of the molecular distributions and abundance ratios seen in our prior observations of the prototypical shocked outflow L1157. In addition, we find that, among others, NH2CHO, HCOOCH3, and CH3CHO have significant post-shock chemistry formation routes that differ from those of many other COMs observed in shocks. Finally, a number of selected species and phenomena are studied here with respect to their usefulness as shock tracers in various astrophysical sources.
Abstract
Here we present high-resolution (15–24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs ...observations of the (
J
= 2–1) and (1–0) lines of
13
CO and C
18
O and the (
J
= 1–0) line of C
17
O for five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermochemical modeling framework based on the CO data, the broadband spectral energy distribution, and the millimeter continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C
17
O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H
2
abundance 10–100 times lower than the interstellar medium ratio. The CO gas distributions between 150 and 400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet–disk interactions, but some CO gaps are 5–10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.
Abstract
Constraining dust properties of planet-forming disks via high-angular-resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may ...be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distribution of solid mass and maximum grain size in each disk, including or excluding dust scattering. In the nonscattering model, the dust surface density and maximum grain size profiles decrease from the inner disks to the outer disks, with local maxima at the bright ring locations, as expected from dust trapping models. The inferred maximum grain sizes from the inner to outer disks decrease from 1 cm to 1 mm. For IM Lup, HD 163296, and MWC 480 in the scattering model, two solutions are compatible with their observed inner disk emission: one solution corresponding to a maximum grain size of a few millimeters (similar to the nonscattering model), and the other corresponding to a size of a few hundred micrometers. Based on the estimated Toomre parameter, only IM Lup—which shows a prominent spiral morphology in millimeter dust—is found to be gravitationally unstable. The estimated maximum Stokes number in all the disks lies between 0.01 and 0.3, and the estimated turbulence parameters in the rings of AS 209 and HD 163296 are close to the threshold where dust growth is limited by turbulent fragmentation. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.
H2CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an Atacama Large Millimeter/submillimeter Array ...survey of H2CO toward 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H2CO is detected toward 13 disks and tentatively detected toward a fourteenth. We find both centrally peaked and centrally depressed emission morphologies, and half of the disks show ring-like structures at or beyond expected CO snowline locations. Together these morphologies suggest that H2CO in disks is commonly produced through both gas-phase and CO-ice-regulated grain-surface chemistry. We extract disk-averaged and azimuthally-averaged H2CO excitation temperatures and column densities for four disks with multiple H2CO line detections. The temperatures are between 20-50 K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H2CO emission in disks generally emerges from the warm molecular layer, with some contributions from the colder midplane. Applying the same H2CO excitation temperatures to all disks in the survey, we find that H2CO column densities span almost three orders of magnitude (∼5 × 1011-5 × 1014 cm−2). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H2CO compared to T Tauri disks, possibly because of less CO freeze-out. More H2CO observations toward Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H2CO and other oxygen-bearing organics efficiently form during planet formation.
We report the first detection of the molecule cyanomethyl, CH2CN, in a protoplanetary disk. Until now, CH2CN had only been observed at earlier evolutionary stages, in the molecular clouds TMC-1, ...Sgr2, and L483, in the prestellar core L1544, and toward the protostar L1527. We detect six transitions of ortho-CH2CN toward the disk around nearby T Tauri star TW Hya. An excitation analysis reveals that the disk-averaged column density, N , for ortho-CH2CN is (6.3 ± 0.5) × 1012 cm−2, which is rescaled to reflect a 3:1 ortho-para ratio, resulting in a total column density, N tot, of (8.4 ± 0.7) × 1012 cm−2. We calculate a disk-average rotational temperature, T rot = 40 ± 5 K, while a radially resolved analysis shows that T rot remains relatively constant across the radius of the disk. This high rotation temperature suggests that in a static disk and if vertical mixing can be neglected, CH2CN is largely formed through gas-phase reactions in the upper layers of the disk, rather than solid-state reactions on the surface of grains in the disk midplane. The integrated intensity radial profiles show a ring structure consistent with molecules such as CN and DCN. We note that this is also consistent with previous lower-resolution observations of centrally peaked CH3CN emission toward the TW Hya disks, since the observed emission gap disappears when convolving our observations with a larger beam size. We obtain a CH2CN/CH3CN ratio ranging between 4 and 10. This high CH2CN/CH3CN is reproduced in a representative chemical model of the TW Hya disk that employs standard static disk chemistry model assumptions, i.e., without any additional tuning.