The First Interferometric Measurements of NH2D/NH3 Ratio in Hot Corinos

The chemical evolution of nitrogen during star and planet formation is still not fully understood. Ammonia (NH3) is a key specie in the understanding of the molecular evolution in star-forming clouds and nitrogen isotope fractionation. In this paper, we present high-spatial-resolution observations of multiple emission lines of NH3 toward the protobinary system NGC1333 IRAS4A with the Karl G. Jansky Very Large Array. We spatially resolved the binary (hereafter, 4A1 and 4A2) and detected compact emission of NH3 transitions with high excitation energies (≳100 K) from the vicinity of the protostars, indicating the NH3 ice has sublimated at the inner hot region. The NH3 column density is estimated to be ∼1017-1018 cm−2. We also detected two NH2D transitions, allowing us to constrain the deuterium fractionation of ammonia. The NH2D/NH3 ratios are as high as ∼0.3-1 in both 4A1 and 4A2. From comparisons with the astrochemical models in the literature, the high NH2D/NH3 ratios suggest that the formation of NH3 ices mainly started in the prestellar phase after the formation of bulk water ice finished, and that the primary nitrogen reservoir in the star-forming cloud could be atomic nitrogen (or N atoms) rather than nitrogen-bearing species such as N2 and NH3. The implications on the physical properties of IRAS4A’s cores are discussed as well

The First Interferometric Measurements of NH₂D/NH₃ Ratio in Hot Corinos

The chemical evolution of nitrogen during star and planet formation is still not fully understood. Ammonia (NH_{3}) is a key specie in the understanding of the molecular evolution in star-forming clouds and nitrogen isotope fractionation. In this paper, we present high-spatial-resolution observations of multiple emission lines of NH_{3} toward the protobinary system NGC1333 IRAS4A with the Karl G. Jansky Very Large Array. We spatially resolved the binary (hereafter, 4A1 and 4A2) and detected compact emission of NH3 transitions with high excitation energies (≳100 K) from the vicinity of the protostars, indicating the NH_{3} ice has sublimated at the inner hot region. The NH3 column density is estimated to be ∼10^{17}–10^{18} cm^{−2}. We also detected two NH_{2}D transitions, allowing us to constrain the deuterium fractionation of ammonia. The NH_{2}D/NH_{3} ratios are as high as ∼0.3–1 in both 4A1 and 4A2. From comparisons with the astrochemical models in the literature, the high NH_{2}D/NH_{3} ratios suggest that the formation of NH3 ices mainly started in the prestellar phase after the formation of bulk water ice finished, and that the primary nitrogen reservoir in the star-forming cloud could be atomic nitrogen (or N atoms) rather than nitrogen-bearing species such as N_{2} and NH_{3}. The implications on the physical properties of IRAS4A's cores are discussed as well

Early Planet Formation in Embedded Disks (eDisk). VIII. A Small Protostellar Disk around the Extremely Low-Mass and Young Class 0 Protostar, IRAS 15398-3359

Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting to find the deconvolved size and $2\sigma$ radius of the dust disk to be $4.5\times2.8\,\mathrm{au}$ and $3.8\,\mathrm{au}$, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be $0.6-1.8\,M_\mathrm{jup}$, indicating a very low-mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the PV diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are $0.022\,M_\odot$ and $31.2\,\mathrm{au}$ from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be $0.1\,M_\odot$. The protostellar mass-accretion rate and the specific angular momentum at the protostellar disk edge are found to be between $1.3-6.1\times10^{-6}\,M_\odot\,\mathrm{yr^{-1}}$ and $1.2-3.8\times10^{-4}\,\mathrm{km\,s^{-1}\,pc}$, respectively, with an age estimated between $0.4-7.5\times10^{4}\,$yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars.Comment: 28 pages, 16 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Progra

Molecules with ALMA at Planet-forming Scales (MAPS). XI. CN and HCN as tracers of photochemistry in disks

Funding: I.C. was supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. A.S.B. acknowledges the studentship funded by the Science and Technology Facilities Council of the United Kingdom (STFC). J.D.I. acknowledges support from the Science and Technology Facilities Council of the United Kingdom (STFC) under ST/T000287/1. C.W. acknowledges financial support from the University of Leeds, STFC, and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, MR/T040726/1).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 0.″3 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.Publisher PDFPeer reviewe

Molecules with ALMA at Planet-forming Scales (MAPS). XIV. Revealing disk substructures in multiwavelength continuum emission

Funding: I.C. was supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. A.S.B acknowledges the studentship funded by the Science and Technology Facilities Council of the United Kingdom (STFC). C.W. acknowledges financial support from the University of Leeds, STFC and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, MR/T040726/1).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.Publisher PDFPeer reviewe

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

Funding: K.Z., K.R.S., J.H., J.B., J.B.B., and I.C. acknowledge the support of NASA through Hubble Fellowship grants HST-HF2-51401.001, HST-HF2-51419.001, HST-HF2-51460.001-A, HSTHF2-51427.001-A, HST-HF2-51429.001-A, and HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. A.S.B. acknowledges the studentship funded by the Science and Technology Facilities Council of the United Kingdom (STFC). C.W. acknowledges financial support from the University of Leeds, STFC, and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, MR/T040726/1).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 ~0"15 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.Publisher PDFPeer reviewe

Molecules with ALMA at Planet-forming Scales (MAPS). XIII. HCO+ and disk ionization structure

Funding: K.Z., K.R.S., J.H., J.B., J.B.B., and I.C. acknowledge the support of NASA through Hubble Fellowship grants HST-HF2-51401.001, HST-HF2-51419.001, HST-HF2-51460.001-A, HST-HF2-51427.001-A, HST-HF2-51429.001-A, and HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. J.D.I. acknowledges support from the Science and Technology Facilities Council of the United Kingdom (STFC) under ST/T000287/1. C.W. acknowledges financial support from the University of Leeds, STFC and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, and MR/T040726/1).We observed HCO+ J = 1 - 0 and H13CO+ J = 1 - 0 emission toward the five protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 as part of the MAPS project. HCO+ is detected and mapped at 0"3 resolution in all five disks, while H13CO+ is detected (S/N > 6σ) toward GM Aur and HD 163296 and tentatively detected (S/N > 3σ) toward the other disks by a matched filter analysis. Inside a radius of R ~ 100 au, the HCO+ column density is flat or shows a central dip. At outer radii (≳100 au), the HCO+ column density decreases outward, while the column density ratio of HCO+/CO is mostly in the range of ~10-5-10-4. We derived the HCO+ abundance in the warm CO-rich layer, where HCO+ is expected to be the dominant molecular ion. At R ≳ 100 au, the HCO+ abundance is ~3 × 10-11 - 3 × 10-10, which is consistent with a template disk model with X-ray ionization. At the smaller radii, the abundance decreases inward, which indicates that the ionization degree is lower in denser gas, especially inside the CO snow line, where the CO-rich layer is in the midplane. Comparison of template disk models with the column densities of HCO+, N2H+, and N2D+ indicates that the midplane ionization rate is ≳10-18 s-1 for the disks around IM Lup, AS 209, and HD 163296. We also find hints of an increased HCO+ abundance around the location of dust continuum gaps in AS 209, HD 163296, and MWC 480. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.Publisher PDFPeer reviewe

Molecules with ALMA at Planet-Forming Scales (MAPS). IX. Distribution and properties of the large organic molecules HC3N, CH3CN, and c-C3H2

Funding: J.D.I. acknowledges support from the Science and Technology Facilities Council of the United Kingdom (STFC) under ST/T000287/1. C.W. acknowledges financial support from the University of Leeds, STFC and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, MR/T040726/1). A.S.B acknowledges the studentship funded by the Science and Technology Facilities Council of the United Kingdom (STFC). I.C. was supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555.The precursors to larger, biologically relevant molecules are detected throughout interstellar space, but determining the presence and properties of these molecules during planet formation requires observations of protoplanetary disks at high angular resolution and sensitivity. Here, we present 0"3 observations of HC3N, CH3CN, and c-C3H2 in five protoplanetary disks observed as part of the Molecules with ALMA at Planet-forming Scales (MAPS) Large Program. We robustly detect all molecules in four of the disks (GM Aur, AS 209, HD 163296, and MWC 480) with tentative detections of c-C3H2 and CH3CN in IM Lup. We observe a range of morphologies-central peaks, single or double rings-with no clear correlation in morphology between molecule or disk. Emission is generally compact and on scales comparable with the millimeter dust continuum. We perform both disk-integrated and radially resolved rotational diagram analysis to derive column densities and rotational temperatures. The latter reveals 5-10 times more column density in the inner 50-100 au of the disks when compared with the disk-integrated analysis. We demonstrate that CH3CN originates from lower relative heights in the disks when compared with HC3N, in some cases directly tracing the disk midplane. Finally, we find good agreement between the ratio of small to large nitriles in the outer disks and comets. Our results indicate that the protoplanetary disks studied here are host to significant reservoirs of large organic molecules, and that this planet- and comet-building material can be chemically similar to that in our own solar system. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.Publisher PDFPeer reviewe

Molecules with ALMA at Planet-forming Scales (MAPS). XII. Inferring the C/O and S/H ratios in protoplanetary disks with sulfur molecules

Funding: I.C. was supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. C.W. acknowledges financial support from the University of Leeds, STFC, and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, MR/T040726/1). J.D.I. acknowledges support from the Science and Technology Facilities Council of the United Kingdom (STFC) under ST/T000287/1.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 ~0"3 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 ~ 0"5 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.Publisher PDFPeer reviewe

Molecules with ALMA at Planet-forming Scales (MAPS). XVII. Determining the 2D thermal structure of the HD 163296 disk

Funding: J.D.I. acknowledges support from the Science and Technology Facilities Council of the United Kingdom (STFC) under ST/T000287/1. C.W. acknowledges financial support from the University of Leeds, STFC, and UKRI (grant Nos. ST/R000549/1, ST/T000287/1, MR/T040726/1). I.C. was supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51405.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555.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 (~0"12 (13 au)-0"25 (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.Publisher PDFPeer reviewe