Showing 1–50 of 96 results for author: Loomis, R

Maser Activity of Organic Molecules toward Sgr B2(N)

Authors: Ci Xue, Anthony Remijan, Alexandre Faure, Emmanuel Momjian, Todd R. Hunter, Ryan A. Loomis, Eric Herbst, Brett McGuire

Abstract: At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH$_2$NH, HNCNH, and HCOOCH$_3$. However, the lack of spatial distribution information of these new maser species has prevented us from assessing the excitation conditions of the maser emission as… ▽ More At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH$_2$NH, HNCNH, and HCOOCH$_3$. However, the lack of spatial distribution information of these new maser species has prevented us from assessing the excitation conditions of the maser emission as well as their pumping mechanisms. Here, we present a mapping study toward Sgr B2 North (N) to locate the region where the complex maser emission originates. We report the first detection of the Class I methanol (CH$_3$OH) maser at 84 GHz and the first interferometric map of the methanimine (CH$_2$NH) maser at 5.29 GHz toward this region. In addition, we present a tool for modeling and fitting the unsaturated molecular maser signals with non-LTE radiative transfer models and Bayesian analysis using the Markov-Chain Monte Carlo approach. These enable us to quantitatively assess the observed spectral profiles. The results suggest a two-chain-clump model for explaining the intense CH$_3$OH Class I maser emission toward a region with low continuum background radiation. By comparing the spatial origin and extent of maser emission from several molecular species, we find that the 5.29 GHz CH$_2$NH maser has a close spatial relationship with the 84 GHz CH$_3$OH Class I masers. This relationship serves as observational evidence to suggest a similar collisional pumping mechanism for these maser transitions. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 17 figures

Rotational Spectrum and First Interstellar Detection of 2-Methoxyethanol Using ALMA Observations of NGC 6334I

Authors: Zachary T. P. Fried, Samer J. El-Abd, Brian M. Hays, Gabi Wenzel, Alex N. Byrne, Laurent Margulès, Roman A. Motiyenko, Steven T. Shipman, Maria P. Horne, Jes K. Jørgensen, Crystal L. Brogan, Todd R. Hunter, Anthony J. Remijan, Andrew Lipnicky, Ryan A. Loomis, Brett A. McGuire

Abstract: We use both chirped-pulse Fourier transform and frequency modulated absorption spectroscopy to study the rotational spectrum of 2-methoxyethanol in several frequency regions ranging from 8.7-500 GHz. The resulting rotational parameters permitted a search for this molecule in Atacama Large Millimeter/submillimeter Array (ALMA) observations toward the massive protocluster NGC 6334I as well as source… ▽ More We use both chirped-pulse Fourier transform and frequency modulated absorption spectroscopy to study the rotational spectrum of 2-methoxyethanol in several frequency regions ranging from 8.7-500 GHz. The resulting rotational parameters permitted a search for this molecule in Atacama Large Millimeter/submillimeter Array (ALMA) observations toward the massive protocluster NGC 6334I as well as source B of the low-mass protostellar system IRAS 16293-2422. 25 rotational transitions are observed in the ALMA Band 4 data toward NGC 6334I, resulting in the first interstellar detection of 2-methoxyethanol. A column density of $1.3_{-0.9}^{+1.4} \times 10^{17}$ cm$^{-2}$ is derived at an excitation temperature of $143_{-39}^{+31}$ K. However, molecular signal is not observed in the Band 7 data toward IRAS 16293-2422B and an upper limit column density of $2.5 \times 10^{15}$ cm$^{-2}$ is determined. Various possible formation pathways--including radical recombination and insertion reactions--are discussed. We also investigate physical differences between the two interstellar sources that could result in the observed abundance variations. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: Accepted for publication in The Astrophysical Journal Letters

The First Spatially-resolved Detection of $^{13}$CN in a Protoplanetary Disk and Evidence for Complex Carbon Isotope Fractionation

Authors: Tomohiro C. Yoshida, Hideko Nomura, Kenji Furuya, Richard Teague, Charles J. Law, Takashi Tsukagoshi, Seokho Lee, Christian Rab, Karin I. Öberg, Ryan A. Loomis

Abstract: Recent measurements of carbon isotope ratios in both protoplanetary disks and exoplanet atmospheres have suggested a possible transfer of significant carbon isotope fractionation from disks to planets. For a clearer understanding of the isotopic link between disks and planets, it is important to measure the carbon isotope ratios in various species. In this paper, we present a detection of the… ▽ More Recent measurements of carbon isotope ratios in both protoplanetary disks and exoplanet atmospheres have suggested a possible transfer of significant carbon isotope fractionation from disks to planets. For a clearer understanding of the isotopic link between disks and planets, it is important to measure the carbon isotope ratios in various species. In this paper, we present a detection of the $^{13}$CN $N=2-1$ hyperfine lines in the TW Hya disk with the Atacama Large Millimeter/submillimeter Array. This is the first spatially-resolved detection of $^{13}$CN in disks, which enables us to measure the spatially resolved $^{12}$CN/$^{13}$CN ratio for the first time. We conducted non-local thermal equilibrium modeling of the $^{13}$CN lines in conjunction with previously observed $^{12}$CN lines to derive the kinetic temperature, ${\rm H_2}$ volume density, and column densities of $^{12}$CN and $^{13}$CN. The ${\rm H_2}$ volume density is found to range between $ (4 - 10)\times10^7 \ {\rm cm^{-3}}$, suggesting that CN molecules mainly reside in the disk upper layer. The $^{12}$CN/$^{13}$CN ratio is measured to be $ 70^{+9}_{-6}$ at $30 < r < 80$ au from the central star, which is similar to the $\rm ^{12}C/^{13}C$ ratio in the interstellar medium. However, this value differs from the previously reported values found for other carbon-bearing molecules (CO and HCN) in the TW Hya disk. This could be self-consistently explained by different emission layer heights for different molecules combined with preferential sequestration of $\rm ^{12}C$ into the solid phase towards the disk midplane. This study reveals the complexity of the carbon isotope fractionation operating in disks. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 16 pages, 7 figures, accepted for publication in ApJ

An Automated Chemical Exploration of NGC 6334I at 340 au Resolution

Authors: Samer J. El-Abd, Crystal L. Brogan, Todd R. Hunter, Kin Long Kelvin Lee, Ryan A. Loomis, Brett A. McGuire

Abstract: Much of the information gleaned from observations of star-forming regions comes from the analysis of their molecular emission spectra, particularly in the radio regime. The time-consuming nature of fitting synthetic spectra to observations interactively for such line-rich sources, however, often results in such analysis being limited to data extracted from a single-dish observation or a handful of… ▽ More Much of the information gleaned from observations of star-forming regions comes from the analysis of their molecular emission spectra, particularly in the radio regime. The time-consuming nature of fitting synthetic spectra to observations interactively for such line-rich sources, however, often results in such analysis being limited to data extracted from a single-dish observation or a handful of pixels from an interferometric observation. Yet, star-forming regions display a wide variety of physical conditions that are difficult, if not impossible, to accurately characterize with such a limited number of spectra. We have developed an automated fitting routine that visits every pixel in the field of view of an ALMA data cube and determines the best-fit physical parameters, including excitation temperature and column densities, for a given list of molecules. In this proof-of-concept work, we provide an overview of the fitting routine and apply it to 0".26, 1.1 km s$^{-1}$ resolution ALMA observations of two sites of massive star-formation in NGC 6334I. Parameters were found for 21 distinct molecules by generating synthetic spectra across 7.48 GHz of spectral bandwidth between 280 and 351 GHz. Spatial images of the derived parameters for each of the > 8000 pixels are presented with special attention paid to the C$_2$H$_4$O$_2$ isomers and their relative variations. We highlight the greater scientific utility of the column density and velocity images of individual molecules compared to traditional moment maps of single transitions. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 40 pages, 71 figures, accepted for publication in The Astrophysical Journal

JWST-MIRI Spectroscopy of Warm Molecular Emission and Variability in the AS 209 Disk

Authors: Carlos E. Muñoz-Romero, Karin I. Öberg, Andrea Banzatti, Klaus M. Pontoppidan, Sean M. Andrews, David J. Wilner, Edwin A. Bergin, Ian Czekala, Charles J. Law, Colette Salyk, Richard Teague, Chunhua Qi, Jennifer B. Bergner, Jane Huang, Catherine Walsh, Viviana V. Guzmán, L. Ilsedore Cleeves, Yuri Aikawa, Jaehan Bae, Alice S. Booth, Gianni Cataldi, John D. Ilee, Romane Le Gal, Feng Long, Ryan A. Loomis , et al. (2 additional authors not shown)

Abstract: We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $μ$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal det… ▽ More We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $μ$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal detections of CO$_2$, HCN, and a possible colder water vapor component. Using slab models with a detailed treatment of opacities and line overlap, we retrieve the column density, emitting area, and excitation temperature of water vapor and OH, and provide upper limits for the observable mass of other molecules. Compared to MIRI spectra of other T-Tauri disks, the inner disk of AS 209 does not appear to be atypically depleted in CO$_2$ nor HCN. Based on \textit{Spitzer IRS} observations, we further find evidence for molecular emission variability over a 10-year baseline. Water, OH, and CO$_2$ line luminosities have decreased by factors 2-4 in the new MIRI epoch, yet there are minimal continuum emission variations. The origin of this variability is yet to be understood. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ

High turbulence in the IM Lup protoplanetary disk: Direct observational constraints from CN and C$_2$H emission

Authors: T. Paneque-Carreño, A. F. Izquierdo, R. Teague, A. Miotello, E. A. Bergin, R. Loomis, E. F. van Dishoeck

Abstract: Constraining turbulence in disks is key towards understanding their evolution through the transport of angular momentum. Until now measurements of high turbulence have remained elusive and methods for estimating turbulence relay mostly on complex radiative transfer models of the data. Using the disk emission from IM Lup, a source proposed to be undergoing Magneto-Rotational Instabilities (MRI) and… ▽ More Constraining turbulence in disks is key towards understanding their evolution through the transport of angular momentum. Until now measurements of high turbulence have remained elusive and methods for estimating turbulence relay mostly on complex radiative transfer models of the data. Using the disk emission from IM Lup, a source proposed to be undergoing Magneto-Rotational Instabilities (MRI) and possibly have high turbulence values in the upper disk layers, we present a new way of directly measuring turbulence without need of radiative transfer or thermochemical models. Through the characterization of the CN and C$_2$H emission in IM Lup, we aim to connect the information on the vertical and thermal structure of a particular disk region to derive turbulence at that location. By using an optically thin tracer it is possible to directly measure turbulence from the non-thermal broadening of the line. The vertical layers of the CN and C$_2$H emission are traced directly from the channel maps using ALFAHOR. By comparing their position to that of optically thick CO observations we are able to characterize the kinetic temperature of the emitting region. Using a simple parametric model of the line intensity with DISCMINER we accurately measure the emission linewidth and separate the thermal and non-thermal components. Assuming that the non-thermal component is fully turbulent, we are able to directly estimate the turbulent motions at the studied radial and vertical location of CN emission. IM Lup shows high turbulence of Mach 0.4-0.6 at $z/r \sim$ 0.25. Considering previous estimates of low turbulence near the midplane, this may indicate a vertical gradient in the disk turbulence, which is a key prediction of MRI studies. CN and C$_2$H are both emitting from a localized upper disk region at $z/r =$0.2-0.3, in agreement with thermochemical models. △ Less

Submitted 6 February, 2024; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A

MAPS: Constraining Serendipitous Time Variability in Protoplanetary Disk Molecular Ion Emission

Authors: Abygail R. Waggoner, L. Ilsedore Cleeves, Ryan A. Loomis, Yuri Aikawa, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Jenny K. Calahan, Gianni Cataldi, Charles J. Law, Romane Le Gal, Feng Long, Karin I. Öberg, Richard Teague, David J. Wilner

Abstract: Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged t… ▽ More Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged the line and continuum data for each individual day the five sources were observed (HD 163296, AS 209, GM Aur, MWC 480, and IM Lup, with between 3 to 6 unique visits per source). Significant enhancement (>3σ) was not observed, but we find variations in the spectral profiles in all five disks. Variations in AS 209, GM Aur, and HD 163296 are tentatively attributed to variations in HCO+ flux, while variations in IM Lup and MWC 480 are most likely introduced by differences in the \textit{uv} coverage, which impact the amount of recovered flux during imaging. The tentative detections and low degree of variability are consistent with expectations of X-ray flare driven HCO+ variability, which requires relatively large flares to enhance the HCO+ rotational emission at significant (>20%) levels. These findings also demonstrate the need for dedicated monitoring campaigns with high signal to noise ratios to fully characterize X-ray flare driven chemistry. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ, 18 pages, 9 figures

CoCCoA: Complex Chemistry in hot Cores with ALMA. Selected oxygen-bearing species

Authors: Y. Chen, M. L. van Gelder, P. Nazari, C. L. Brogan, E. F. van Dishoeck, H. Linnartz, J. K. Jørgensen, T. R. Hunter, O. H. Wilkins, G. A. Blake, P. Caselli, K. -J. Chuang, C. Codella, I. Cooke, M. N. Drozdovskaya, R. T. Garrod, S. Ioppolo, M. Jin, B. M. Kulterer, N. F. W. Ligterink, A. Lipnicky, R. Loomis, M. G. Rachid, S. Spezzano, B. A. McGuire

Abstract: Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protos… ▽ More Complex organic molecules (COMs) have been observed to be abundant in the gas phase toward protostars. Deep line surveys have been carried out only for a limited number of well-known high-mass star forming regions using the Atacama Large Millimeter/submillimeter Array (ALMA), which has unprecedented resolution and sensitivity. Statistical studies on oxygen-bearing COMs (O-COMs) in high-mass protostars using ALMA are still lacking. With the recent CoCCoA survey, we are able to determine the column density ratios of six O-COMs with respect to methanol (CH$_3$OH) in a sample of 14 high-mass protostellar sources to investigate their origin through ice and/or gas-phase chemistry. The selected species are: acetaldehyde (CH$_3$CHO), ethanol (C$_2$H$_5$OH), dimethyl ether (DME, CH$_3$OCH$_3$), methyl formate (MF, CH$_3$OCHO), glycolaldehyde (GA, CH$_2$OHCHO), and ethylene glycol (EG, (CH$_2$OH)$_2$). DME and MF have the highest and most constant ratios within one order of magnitude, while the other four species have lower ratios and exhibit larger scatter by 1-2 orders of magnitude. We compare the O-COM ratios of high-mass CoCCoA sources with those of 5 low-mass protostars available from the literature, along with the results from experiments and simulations. We find that the O-COM ratios with respect to methanol are on the same level in both the high- and low-mass samples, which suggests that these species are mainly formed in similar environments during star formation, probably in ice mantles on dust grains during early pre-stellar stages. Current simulations and experiments can reproduce most observational trends with a few exceptions, and hypotheses exist to explain the differences between observations and simulations/experiments, such as the involvement of gas-phase chemistry and different emitting areas of molecules. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 17 pages and 7 figures in the main body, 15 pages and 8 figures in the appendix

Journal ref: A&A 678, A137 (2023)

The ALMA Interferometric Pipeline Heuristics

Authors: Todd R. Hunter, Remy Indebetouw, Crystal L. Brogan, Kristin Berry, Chin-Shin Chang, Harold Francke, Vincent C. Geers, Laura Gómez, John E. Hibbard, Elizabeth M. Humphreys, Brian R. Kent, Amanda A. Kepley, Devaky Kunneriath, Andrew Lipnicky, Ryan A. Loomis, Brian S. Mason, Joseph S. Masters, Luke T. Maud, Dirk Muders, Jose Sabater, Kanako Sugimoto, László Szűcs, Eugene Vasiliev, Liza Videla, Eric Villard , et al. (3 additional authors not shown)

Abstract: We describe the calibration and imaging heuristics developed and deployed in the ALMA interferometric data processing pipeline, as of ALMA Cycle 9. The pipeline software framework is written in Python, with each data reduction stage layered on top of tasks and toolkit functions provided by the Common Astronomy Software Applications package. This framework supports a variety of tasks for observator… ▽ More We describe the calibration and imaging heuristics developed and deployed in the ALMA interferometric data processing pipeline, as of ALMA Cycle 9. The pipeline software framework is written in Python, with each data reduction stage layered on top of tasks and toolkit functions provided by the Common Astronomy Software Applications package. This framework supports a variety of tasks for observatory operations, including science data quality assurance, observing mode commissioning, and user reprocessing. It supports ALMA and VLA interferometric data along with ALMA and NRO45m single dish data, via different stages and heuristics. In addition to producing calibration tables, calibrated measurement sets, and cleaned images, the pipeline creates a WebLog which serves as the primary interface for verifying the data quality assurance by the observatory and for examining the contents of the data by the user. Following the adoption of the pipeline by ALMA Operations in 2014, the heuristics have been refined through annual development cycles, culminating in a new pipeline release aligned with the start of each ALMA Cycle of observations. Initial development focused on basic calibration and flagging heuristics (Cycles 2-3), followed by imaging heuristics (Cycles 4-5), refinement of the flagging and imaging heuristics with parallel processing (Cycles 6-7), addition of the moment difference analysis to improve continuum channel identification (2020 release), addition of a spectral renormalization stage (Cycle 8), and improvement in low SNR calibration heuristics (Cycle 9). In the two most recent Cycles, 97% of ALMA datasets were calibrated and imaged with the pipeline, ensuring long-term automated reproducibility. We conclude with a brief description of plans for future additions, including self-calibration, multi-configuration imaging, and calibration and imaging of full polarization data. △ Less

Submitted 25 July, 2023; v1 submitted 12 June, 2023; originally announced June 2023.

Comments: Accepted for publication by Publications of the Astronomical Society of the Pacific, 65 pages, 20 figures, 10 tables, 2 appendices. Small corrections and additions from the proof stage have been applied

Journal ref: Publications of the Astronomical Society of the Pacific, volume 135, number 1049, 074501, 2023 July 24

Detection of Interstellar $E$-1-cyano-1,3-butadiene in GOTHAM Observations of TMC-1

Authors: Ilsa R. Cooke, Ci Xue, P. Bryan Changala, Hannah Toru Shay, Alex N. Byrne, Qi Yu Tang, Zachary T. P. Fried, Kin Long Kelvin Lee, Ryan A. Loomis, Thanja Lamberts, Anthony Remijan, Andrew M. Burkhardt, Eric Herbst, Michael C. McCarthy, Brett A. McGuire

Abstract: We report the detection of the lowest energy conformer of $E$-1-cyano-1,3-butadiene ($E$-1-C$_4$H$_5$CN), a linear isomer of pyridine, using the fourth data reduction of the GOTHAM deep spectral survey toward TMC-1 with the 100 m Green Bank Telescope. We performed velocity stacking and matched filter analyses using Markov chain Monte Carlo simulations and find evidence for the presence of this mol… ▽ More We report the detection of the lowest energy conformer of $E$-1-cyano-1,3-butadiene ($E$-1-C$_4$H$_5$CN), a linear isomer of pyridine, using the fourth data reduction of the GOTHAM deep spectral survey toward TMC-1 with the 100 m Green Bank Telescope. We performed velocity stacking and matched filter analyses using Markov chain Monte Carlo simulations and find evidence for the presence of this molecule at the 5.1$σ$ level. We derive a total column density of $3.8^{+1.0}_{-0.9}\times 10^{10}$ cm$^{-2}$, which is predominantly found toward two of the four velocity components we observe toward TMC-1. We use this molecule as a proxy for constraining the gas-phase abundance of the apolar hydrocarbon 1,3-butadiene. Based on the three-phase astrochemical modeling code NAUTILUS and an expanded chemical network, our model underestimates the abundance of cyano-1,3-butadiene by a factor of 19, with a peak column density of $2.34 \times 10^{10}\ \mathrm{cm}^{-2}$ for 1,3-butadiene. Compared to the modeling results obtained in previous GOTHAM analyses, the abundance of 1,3-butadiene is increased by about two orders of magnitude. Despite this increase, the modeled abundances of aromatic species do not appear to change and remain underestimated by 1--4 orders of magnitude. Meanwhile, the abundances of the five-membered ring molecules increase proportionally with 1,3-butadiene by two orders of magnitudes. We discuss implications for bottom-up formation routes to aromatic and polycyclic aromatic molecules. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Astronomical Detection of the Interstellar Anion C10H- towards TMC-1 from the GOTHAM Large Program on the GBT

Authors: Anthony Remijan, Haley N. Scolati, Andrew M. Burkhardt, P. Bryan Changala, Steven B. Charnley, Ilsa R. Cooke, Martin A. Cordiner, Harshal Gupta, Eric Herbst, Kin Long Kelvin Lee, Ryan Loomis, Christopher N. Shingledecker, Mark A. Siebert, Ci Xue, Michael C. McCarthy, Brett A. McGuire

Abstract: Using data from the GOTHAM (GBT Observations of TMC-1: Hunting for Aromatic Molecules) survey, we report the first astronomical detection of the C10H- anion. The astronomical observations also provided the necessary data to refine the spectroscopic parameters of C10H-. From the velocity stacked data and the matched filter response, C10H- is detected at >9σ confidence level at a column density of 4… ▽ More Using data from the GOTHAM (GBT Observations of TMC-1: Hunting for Aromatic Molecules) survey, we report the first astronomical detection of the C10H- anion. The astronomical observations also provided the necessary data to refine the spectroscopic parameters of C10H-. From the velocity stacked data and the matched filter response, C10H- is detected at >9σ confidence level at a column density of 4.04e11 cm-2. A dedicated search for the C10H radical was also conducted towards TMC-1. In this case, the stacked molecular emission of C10H was detected at a ~3.2σ confidence interval at a column density of 2.02e11 cm-2. However, since the determined confidence level is currently <5σ, we consider the identification of C10H as tentative. The full GOTHAM dataset was also used to better characterize the physical parameters including column density, excitation temperature, linewidth, and source size for the C4H, C6H and C8H radicals and their respective anions, and the measured column densities were compared to the predictions from a gas/grain chemical formation model and from a machine learning analysis. Given the measured values, the C10H-/C10H column density ratio is ~2.0 - the highest value measured between an anion and neutral species to date. Such a high ratio is at odds with current theories for interstellar anion chemistry. For the radical species, both models can reproduce the measured abundances found from the survey; however, the machine learning analysis matches the detected anion abundances much better than the gas/grain chemical model, suggesting that the current understanding of the formation chemistry of molecular anions is still highly uncertain. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 38 Pages, 24 Figures, 12 Tables, 8 Appendices

MSC Class: 85-11 ACM Class: A.1

A kinematically detected planet candidate in a transition disk

Authors: Jochen Stadler, Myriam Benisty, Andrés F. Izquierdo, Stefano Facchini, Richard Teague, Nicolas Kurtovic, Paola Pinilla, Jaehan Bae, Megan Ansdell, Ryan Loomis, Satoshi Mayama, Laura M. Pérez, Leonardo Testi

Abstract: Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps an… ▽ More Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps and moment maps of the line of sight velocity and peak intensity. We fitted a Keplerian model of the channel-by-channel emission to study line profile differences and produced deprojected radial profiles for all velocity components. The $^{12}$CO emission shows a cavity inwards of $\sim$56 au and within the dust continuum ring at 81 au. The azimuthal brightness variations in the $^{12}$CO line and dust continuum are broadly aligned with the shadows detected in scattered-light observations. We find a strong localized non-Keplerian feature toward the west within the continuum ring (at $R=41\pm 10$ au and $PA=280\pm 2 ^\circ$). A tightly wound spiral is also detected which extends over 300$^\circ$ in azimuth, possibly connected to the localized non-Keplerian feature. Finally, a bending of the iso-velocity contours within the gas cavity indicates a highly perturbed inner region, possibly related to the presence of a misaligned inner disk. While broadly aligned with the scattered-light shadows, the localized non-Keplerian feature cannot be solely due to changes in temperature. Instead, we interpret the kinematical feature as tracing a massive companion located at the edge of the dust continuum ring. We speculate that the spiral is caused by buoyancy resonances driven by planet-disk interactions. However, this potential planet at $\sim$41 au cannot explain the gas-depleted cavity, the low accretion rate, and the misaligned inner disk, which suggests the presence of another companion closer in. △ Less

Submitted 11 January, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: 13 pages, 14 figures. To be published in Astronomy & Astrophysics Letters

Journal ref: A&A 670, L1 (2023)

Mapping Protoplanetary Disk Vertical Structure with CO Isotopologue Line Emission

Authors: Charles J. Law, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Myriam Benisty, Stefano Facchini, Kevin Flaherty, Andrea Isella, Sheng Jin, Jun Hashimoto, Jane Huang, Ryan A. Loomis, Feng Long, Carlos E. Muñoz-Romero, Teresa Paneque-Carreño, Laura M. Pérez, Chunhua Qi, Kamber R. Schwarz, Jochen Stadler, Takashi Tsukagoshi, David J. Wilner, Gerrit van der Plas

Abstract: High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$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… ▽ More High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$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 2D (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 (${\lesssim}$0."1 to ${\approx}$0."2; ${\approx}$15-30 au) observations of CO, $^{13}$CO, and C$^{18}$O in either or both J=2-1 and J=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 in $^{13}$CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C$^{18}$O in any disk, which is instead consistent with C$^{18}$O emission originating from the midplane. Both the J=2-1 and J=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 that $^{13}$CO 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. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 31 pages, 15 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/7430257

Regularized Maximum Likelihood Image Synthesis and Validation for ALMA Continuum Observations of Protoplanetary Disks

Authors: Brianna Zawadzki, Ian Czekala, Ryan A. Loomis, Tyler Quinn, Hannah Grzybowski, Robert C. Frazier, Jeff Jennings, Kadri M. Nizam, Yina Jian

Abstract: Regularized Maximum Likelihood (RML) techniques are a class of image synthesis methods that achieve better angular resolution and image fidelity than traditional methods like CLEAN for sub-mm interferometric observations. To identify best practices for RML imaging, we used the GPU-accelerated open source Python package MPoL, a machine learning-based RML approach, to explore the influence of common… ▽ More Regularized Maximum Likelihood (RML) techniques are a class of image synthesis methods that achieve better angular resolution and image fidelity than traditional methods like CLEAN for sub-mm interferometric observations. To identify best practices for RML imaging, we used the GPU-accelerated open source Python package MPoL, a machine learning-based RML approach, to explore the influence of common RML regularizers (maximum entropy, sparsity, total variation, and total squared variation) on images reconstructed from real and synthetic ALMA continuum observations of protoplanetary disks. We tested two different cross-validation (CV) procedures to characterize their performance and determine optimal prior strengths, and found that CV over a coarse grid of regularization strengths easily identifies a range of models with comparably strong predictive power. To evaluate the performance of RML techniques against a ground truth image, we used MPoL on a synthetic protoplanetary disk dataset and found that RML methods successfully resolve structures at fine spatial scales present in the original simulation. We used ALMA DSHARP observations of the protoplanetary disk around HD 143006 to compare the performance of MPoL and CLEAN, finding that RML imaging improved the spatial resolution of the image by up to a factor of 3 without sacrificing sensitivity. We provide general recommendations for building an RML workflow for image synthesis of ALMA protoplanetary disk observations, including effective use of CV. Using these techniques to improve the imaging resolution of protoplanetary disk observations will enable new science, including the detection of protoplanets embedded in disks. △ Less

Submitted 16 June, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 27 pages, 12 figures, accepted for publication in PASP

Discovery of Interstellar 2-Cyanoindene (2-C$_9$H$_7$CN) in GOTHAM Observations of TMC-1

Authors: Madelyn L. Sita, P. Bryan Changala, Ci Xue, Andrew M. Burkhardt, Christopher N. Shingledecker, Kin Long Kelvin Lee, Ryan A. Loomis, Emmanuel Momjian, Mark A. Siebert, Divita Gupta, Eric Herbst, Anthony J. Remijan, Michael C. McCarthy, Ilsa R. Cooke, Brett A. McGuire

Abstract: We present laboratory rotational spectroscopy of five isomers of cyanoindene (2-, 4-, 5-, 6-, and 7-cyanoindene) using a cavity Fourier-transform microwave spectrometer operating between 6-40 GHz. Based on these measurements, we report the detection of 2-cyanoindene (1H-indene-2-carbonitrile; 2-C$_9$H$_7$CN) in GOTHAM line survey observations of the dark molecular cloud TMC-1 using the Green Bank… ▽ More We present laboratory rotational spectroscopy of five isomers of cyanoindene (2-, 4-, 5-, 6-, and 7-cyanoindene) using a cavity Fourier-transform microwave spectrometer operating between 6-40 GHz. Based on these measurements, we report the detection of 2-cyanoindene (1H-indene-2-carbonitrile; 2-C$_9$H$_7$CN) in GOTHAM line survey observations of the dark molecular cloud TMC-1 using the Green Bank Telescope at centimeter wavelengths. Using a combination of Markov Chain Monte Carlo (MCMC), spectral stacking, and matched filtering techniques, we find evidence for the presence of this molecule at the 6.3$σ$ level. This provides the first direct observation of the ratio of a cyano-substituted polycyclic aromatic hydrocarbon (PAH) to its pure hydrocarbon counterpart, in this case indene, in the same source. We discuss the possible formation chemistry of this species, including why we have only detected one of the isomers in TMC-1. We then examine the overall hydrocarbon:CN-substituted ratio across this and other simpler species, as well as compare to those ratios predicted by astrochemical models. We conclude that while astrochemical models are not yet sufficiently accurate to reproduce absolute abundances of these species, they do a good job at predicting the ratios of hydrocarbon:CN-substituted species, further solidifying -CN tagged species as excellent proxies for their fully-symmetric counterparts. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted in ApJL

ALMA Detection of Dust Trapping around Lagrangian Points in the LkCa 15 Disk

Authors: Feng Long, Sean M. Andrews, Shangjia Zhang, Chunhua Qi, Myriam Benisty, Stefano Facchini, Andrea Isella, David J. Wilner, Jaehan Bae, Jane Huang, Ryan A. Loomis, Karin I. Öberg, Zhaohuan Zhu

Abstract: We present deep high-resolution ($\sim$50 mas, 8 au) ALMA 0.88 and 1.3 mm continuum observations of the LkCa 15 disk. The emission morphology shows an inner cavity and three dust rings at both wavelengths, but with slightly narrower rings at the longer wavelength. Along a faint ring at 42 au, we identify two excess emission features at $\sim$10$σ$ significance at both wavelengths: one as an unreso… ▽ More We present deep high-resolution ($\sim$50 mas, 8 au) ALMA 0.88 and 1.3 mm continuum observations of the LkCa 15 disk. The emission morphology shows an inner cavity and three dust rings at both wavelengths, but with slightly narrower rings at the longer wavelength. Along a faint ring at 42 au, we identify two excess emission features at $\sim$10$σ$ significance at both wavelengths: one as an unresolved clump and the other as an extended arc, separated by roughly 120 degrees in azimuth. The clump is unlikely to be a circumplanetary disk (CPD) as the emission peak shifts between the two wavelengths even after accounting for orbital motion. Instead, the morphology of the 42 au ring strongly resembles the characteristic horseshoe orbit produced in planet--disk interaction models, where the clump and the arc trace dust accumulation around Lagrangian points $L_{4}$ and $L_{5}$, respectively. The shape of the 42 au ring, dust trapping in the outer adjacent ring, and the coincidence of the horseshoe ring location with a gap in near-IR scattered light, are all consistent with the scenario of planet sculpting, with the planet likely having a mass between those of Neptune and Saturn. We do not detect point-like emission associated with a CPD around the putative planet location ($0.''27$ in projected separation from the central star at a position angle of $\sim$60\degr), with upper limits of 70 and 33 $μ$Jy at 0.88 and 1.3 mm, respectively, corresponding to dust mass upper limits of 0.02--0.03 $M_{\oplus}$. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: ApJL in press

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

Authors: Jaehan Bae, Richard Teague, Sean M. Andrews, Myriam Benisty, Stefano Facchini, Maria Galloway-Sprietsma, Ryan A. Loomis, Yuri Aikawa, Felipe Alarcon, Edwin Bergin, Jennifer B. Bergner, Alice S. Booth, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzman, Jane Huang, John D. Ilee, Nicolas T. Kurtovic, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Francois Menard, Karin I. Oberg , et al. (7 additional authors not shown)

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^\circ$), 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… ▽ More 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^\circ$), 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\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ 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 $\gtrsim35$ 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 $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-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. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted for publication in the ApJ Letters (July 7, 2022), 19 pages, 13 figures, interactive figures (Figure 7, 8, 9) are available at http://jaehanbae.com/as209/

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

Authors: Charles J. Law, Sage Crystian, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Kevin Flaherty, Viviana V. Guzmán, Jane Huang, John D. Ilee, Joel H. Kastner, Ryan A. Loomis, Feng Long, Laura M. Pérez, Sebastián Pérez, Chunhua Qi, Giovanni P. Rosotti, Dary Ruíz-Rodríguez, Takashi Tsukagoshi, David J. Wilner

Abstract: High spatial resolution CO observations of mid-inclination (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 o… ▽ More High spatial resolution CO observations of mid-inclination (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 ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: 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$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$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\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute 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 ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/6410045

A Search for Heterocycles in GOTHAM Observations of TMC-1

Authors: Timothy J. Barnum, Mark A. Siebert, Kin Long Kelvin Lee, Ryan A. Loomis, P. Bryan Changala, Steven B. Charnley, Madelyn L. Sita, Ci Xue, Anthony J. Remijan, Andrew M. Burkhardt, Brett A. McGuire, Ilsa R. Cooke

Abstract: We have conducted an extensive search for nitrogen-, oxygen- and sulfur-bearing heterocycles toward Taurus Molecular Cloud 1 (TMC-1) using the deep, broadband centimeter-wavelength spectral line survey of the region from the GOTHAM large project on the Green Bank Telescope. Despite their ubiquity in terrestrial chemistry, and the confirmed presence of a number of cyclic and polycyclic hydrocarbon… ▽ More We have conducted an extensive search for nitrogen-, oxygen- and sulfur-bearing heterocycles toward Taurus Molecular Cloud 1 (TMC-1) using the deep, broadband centimeter-wavelength spectral line survey of the region from the GOTHAM large project on the Green Bank Telescope. Despite their ubiquity in terrestrial chemistry, and the confirmed presence of a number of cyclic and polycyclic hydrocarbon species in the source, we find no evidence for the presence of any heterocyclic species. Here, we report the derived upper limits on the column densities of these molecules obtained by Markov Chain Monte Carlo (MCMC) analysis and compare this approach to traditional single-line upper limit measurements. We further hypothesize why these molecules are absent in our data, how they might form in interstellar space, and the nature of observations that would be needed to secure their detection. △ Less

Submitted 10 April, 2022; originally announced April 2022.

Comments: Accepted in the Journal of Physical Chemistry A

CH3-Terminated Carbon Chains in the GOTHAM Survey of TMC-1: Evidence of Interstellar CH3C7N

Authors: Mark A. Siebert, Kin Long Kelvin Lee, Anthony J. Remijan, Andrew M. Burkhardt, Ryan A. Loomis, Michael C. McCarthy, Brett A. McGuire

Abstract: We report a systematic study of all known methyl carbon chains toward TMC-1 using the second data release of the GOTHAM survey, as well as a search for larger species. Using Markov-Chain Monte Carlo simulations and spectral line stacking of over 30 rotational transitions, we report statistically significant emission from methylcyanotriacetylene (CH$_3$C$_7$N) at a confidence level of 4.6$σ$, and u… ▽ More We report a systematic study of all known methyl carbon chains toward TMC-1 using the second data release of the GOTHAM survey, as well as a search for larger species. Using Markov-Chain Monte Carlo simulations and spectral line stacking of over 30 rotational transitions, we report statistically significant emission from methylcyanotriacetylene (CH$_3$C$_7$N) at a confidence level of 4.6$σ$, and use it to derive a column density of ${\sim}$10$^{11}$ cm$^{-2}$. We also searched for the related species, methyltetraacetylene (CH$_3$C$_8$H), and place upper limits on the column density of this molecule. By carrying out the above statistical analyses for all other previously detected methyl-terminated carbon chains that have emission lines in our survey, we assess the abundances, excitation conditions, and formation chemistry of methylpolyynes (CH3C$_{2n}$H) and methylcyanopolyynes (CH3C$_{2n-1}$N) in TMC-1, and compare those with predictions from a chemical model. Based on our observed trends in column density and relative populations of the A and E nuclear spin isomers, we find that the methylpolyynes and methylcyanopolyynes families exhibit stark differences from one another, pointing to separate interstellar formation pathways, which is confirmed through gas-grain chemical modeling with nautilus. △ Less

Submitted 14 January, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 15 pages, 7 figures. Accepted for publication in ApJ

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

Authors: Jennifer B. Bergner, Karin I. Oberg, Viviana V. Guzman, Charles J. Law, Ryan A. Loomis, Gianni Cataldi, Arthur D. Bosman, Yuri Aikawa, Sean M. Andrews, Edwin A. Bergin, Alice S. Booth, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, John D. Ilee, Romane Le Gal, Feng Long, Hideko Nomura, Francois Menard, Chunhua Qi, Kamber R. Schwarz, Richard Teague, Takashi Tsukagoshi, Catherine Walsh, David J. Wilner , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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 towards five disk systems. All disks show bright CN emission within $\sim$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. △ Less

Submitted 16 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Molecules with ALMA at Planet-forming Scales (MAPS) XVI: Characterizing the impact of the molecular wind on the evolution of the HD 163296 system

Authors: Alice S. Booth, Benoit Tabone, John D. Ilee, Catherine Walsh, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Jennifer B. Bergner, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzman, Jane Huang, Charles J. Law, Romane Le Gal, Feng Long, Ryan A. Loomis, Francois Menard, Karin I. Oberg, Chunhua Qi, Kamber R. Schwarz, Richard Teague , et al. (4 additional authors not shown)

Abstract: 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) di… ▽ More 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 through 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. △ Less

Submitted 15 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: Accepted ApJ July 30th 2021 This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS). X. Studying deuteration at high angular resolution toward protoplanetary disks

Authors: Gianni Cataldi, Yoshihide Yamato, Yuri Aikawa, Jennifer B. Bergner, Kenji Furuya, Viviana V. Guzmán, Jane Huang, Ryan A. Loomis, Chunhua Qi, Sean M. Andrews, Edwin A. Bergin, Alice S. Booth, Arthur D. Bosman, L. Ilsedore Cleeves, Ian Czekala, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Hideko Nomura, Karin I. Öberg, Kamber R. Schwarz, Richard Teague , et al. (4 additional authors not shown)

Abstract: Deuterium fractionation is dependent on various physical and chemical parameters. Thus, the formation location and thermal history of material in the solar system is often studied by measuring its D/H ratio. This requires knowledge about the deuteration processes operating during the planet formation era. We aim to study these processes by radially resolving the DCN/HCN (at 0.3" resolution) and N… ▽ More Deuterium fractionation is dependent on various physical and chemical parameters. Thus, the formation location and thermal history of material in the solar system is often studied by measuring its D/H ratio. This requires knowledge about the deuteration processes operating during the planet formation era. We aim to study these processes by radially resolving the DCN/HCN (at 0.3" resolution) and N$_2$D$^+$/N$_2$H$^+$ (0.3 to 0.9") column density ratios toward the five protoplanetary disks observed by the Molecules with ALMA at Planet-forming scales (MAPS) Large Program. DCN is detected in all five sources, with one newly reported detection. N$_2$D$^+$ is detected in four sources, two of which are newly reported detections. We derive column density profiles that allow us to study the spatial variation of the DCN/HCN and N$_2$D$^+$/N$_2$H$^+$ ratios at high resolution. DCN/HCN varies considerably for different parts of the disks, ranging from $10^{-3}$ to $10^{-1}$. In particular, the inner disk regions generally show significantly lower HCN deuteration compared with the outer disk. In addition, our analysis confirms that two deuterium fractionation channels are active, which can alter the D/H ratio within the pool of organic molecules. N$_2$D$^+$ is found in the cold outer regions beyond $\sim$50 au, with N$_2$D$^+$/N$_2$H$^+$ ranging between $10^{-2}$ and 1 across the disk sample. This is consistent with the theoretical expectation that N$_2$H$^+$ deuteration proceeds via the low-temperature channel only. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 28 November, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: Accepted for publication in the Astrophysical Journal Supplement. 55 pages, 30 figures. Replacement of earlier version with updated references

Journal ref: The Astrophysical Journal Supplement Series, 2021, Volume 257, Issue 1

Molecules with ALMA at Planet-forming Scales (MAPS) XIV: Revealing disk substructures in multi-wavelength continuum emission

Authors: Anibal Sierra, Laura M. Pérez, Ke Zhang, Charles J. Law, Viviana V. Guzmán, Chunhua Qi, Arthur D. Bosman, Karin I. Öberg, Sean M. Andrews, Feng Long, Richard Teague, Alice S. Booth, Catherine Walsh, David J. Wilner, François Ménard, Gianni Cataldi, Ian Czekala, Jaehan Bae, Jane Huang, Jennifer B. Bergner, John D. Ilee, Myriam Benisty, Romane Le Gal, Ryan A. Loomis, Takashi Tsukagoshi , et al. (3 additional authors not shown)

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 distributi… ▽ More 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 few hundred micrometer sizes. 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. △ Less

Submitted 20 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

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

Authors: Yuri Aikawa, Gianni Cataldi, Yoshihide Yamato, Ke Zhang, Alice S. Booth, Kenji Furuya, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Jennifer B. Bergner, Arthur D. Bosman, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzmán, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Ryan A. Loomis, Francois Ménard, Hideko Nomura, Karin I. Öberg, Chunhua Qi, Kamber R. Schwarz, Richard Teague , et al. (3 additional authors not shown)

Abstract: We observed HCO$^+$ $J=1-0$ and H$^{13}$CO$^+$ $J=1-0$ emission towards 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\arcsec\,resolution in all five disks, while H$^{13}$CO$^+$ is detected (SNR$>6 σ$) towards GM Aur and HD 163296 and tentatively detected (SNR$>3 σ$) towards the other disks by a… ▽ More We observed HCO$^+$ $J=1-0$ and H$^{13}$CO$^+$ $J=1-0$ emission towards 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\arcsec\,resolution in all five disks, while H$^{13}$CO$^+$ is detected (SNR$>6 σ$) towards GM Aur and HD 163296 and tentatively detected (SNR$>3 σ$) towards the other disks by a matched filter analysis. Inside a radius of $R\sim 100$ au, the HCO$^+$ column density is flat or shows a central dip. At outer radii ($\gtrsim 100$ au), the HCO$^+$ column density decreases outwards, while the column density ratio of HCO$^+$/CO is mostly in the range of $\sim 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\gtrsim 100$ au, the HCO$^+$ abundance is $\sim 3 \times 10^{-11} - 3\times 10^{-10}$, which is consistent with a template disk model with X-ray ionization. At the smaller radii, the abundance decreases inwards, 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$^+$, N$_2$H$^+$, and N$_2$D$^+$ indicates that the midplane ionization rate is $\gtrsim 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. △ Less

Submitted 20 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: accepted to ApJS, 33 pages, 20 figures

Molecules with ALMA at Planet-forming Scales (MAPS) VI: Distribution of the small organics HCN, C2H, and H2CO

Authors: Viviana V. Guzmán, Jennifer B. Bergner, Charles J. Law, Karin I. Oberg, Catherine Walsh, Gianni Cataldi, Yuri Aikawa, Edwin A. Bergin, Ian Czekala, Jane Huang, Sean M. Andrews, Ryan A. Loomis, Ke Zhang, Romane Le Gal, Felipe Alarcón, John D. Ilee, Richard Teague, L. Ilsedore Cleeves, David J. Wilner, Feng Long, Kamber R. Schwarz, Arthur D. Bosman, Laura M. Pérez, François Ménard, Yao Liu

Abstract: 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… ▽ More 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. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 23 pages, 10 figures, 4 tables. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Molecules with ALMA at Planet-forming Scales (MAPS). IX. Distribution and Properties of the Large Organic Molecules HC$_3$N, CH$_3$CN, and $c$-C$_3$H$_2$

Authors: John D. Ilee, Catherine Walsh, Alice S. Booth, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Jennifer B. Bergner, Arthur D. Bosman, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzmán, Jane Huang, Charles J. Law, Romane Le Gal, Ryan A. Loomis, François Ménard, Hideko Nomura, Karin I Öberg, Chunhua Qi, Kamber R. Schwarz, Richard Teague, Takashi Tsukagoshi, David J. Wilner , et al. (2 additional authors not shown)

Abstract: 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 HC$_3$N, CH$_3$CN, and $c$-C$_3$H$_2$ in five protoplanetary disks observed as part… ▽ More 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 HC$_3$N, CH$_3$CN, and $c$-C$_3$H$_2$ 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$-C$_3$H$_2$ and CH$_3$CN in IM Lup. We observe a range of morphologies -- central peaks, single or double rings -- with no clear correlation in morphology between molecule nor disk. Emission is generally compact and on scales comparable with the millimetre 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 CH$_3$CN originates from lower relative heights in the disks when compared with HC$_3$N, 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 Series. △ Less

Submitted 15 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 24 pages, 11 figures, 5 tables. Accepted for publication in ApJSS. Updated to cross-reference other MAPS publications

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

Authors: Romane Le Gal, Karin I. Öberg, Richard Teague, Ryan A. Loomis, Charles J. Law, Catherine Walsh, Edwin A. Bergin, Francois Menard, David J. Wilner, Sean M. Andrews, Yuri Aikawa, Alice S. Booth, Gianni Cataldi, Jennifer B. Bergner, Arthur D. Bosman, L. Ilsedore Cleeves, Ian Czekala, Kenji Furuya, Viviana V. Guzmán, Jane Huang, John D. Ilee, Hideko Nomura, Chunhua Qi, Kamber R. Schwarz, Takashi Tsukagoshi , et al. (2 additional authors not shown)

Abstract: 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$) observatio… ▽ More 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 $\sim0."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 C$_2$S and SO were also serendipitously covered but only upper limits are found. For MWC 480, we present complementary ALMA observations at $\sim0."5$, of CS, $^{13}$CS, C$^{34}$S, H$_2$CS, OCS, and SO$_2$. We find a column density ratio N(H$_{2}$CS)/N(CS)$\sim2/3$, suggesting that a substantial part of the sulfur reservoir in disks is in organic form (i.e., C$_x$H$_y$S$_z$). 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 super-solar 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. △ Less

Submitted 17 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement (27 pages, 13 figures, 5 tables)

Molecules with ALMA at Planet-forming Scales (MAPS) I: Program Overview and Highlights

Authors: Karin I. Oberg, Viviana V. Guzman, Catherine Walsh, Yuri Aikawa, Edwin A. Bergin, Charles J. Law, Ryan A. Loomis, Felipe Alarcon, Sean M. Andrews, Jaehan Bae, Jennifer B. Bergner, Yann Boehler, Alice S. Booth, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Kenji Furuya, Jane Huang, John D. Ilee, Nicolas T. Kurtovic, Romane Le Gal, Yao Liu, Feng Long , et al. (13 additional authors not shown)

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 fro… ▽ More 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 4 spectral setups, which together cover ~50 lines from over 20 different species. This paper introduces the ApJS 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 sub-structures, 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. △ Less

Submitted 16 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: Accepted for publication in the ApJS MAPS Special Issue. v2 has updated MAPS references and a correction to Fig. 3

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

Authors: Felipe Alarcón, Arthur Bosman, Edwin Bergin, Ke Zhang, Richard Teague, Jaehan Bae, Yuri Aikawa, Sean M. Andrews, Alice Booth, Jenny Calahan, Gianni Cataldi, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Ryan A. Loomis, François Ménard, Karin Öberg, Kamber R. Schwarz, Merel L. R. Van't Hoff, Catherine Walsh, David J. Wilner

Abstract: 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 C$^{18}$O J=2-1 emission exhibit a… ▽ More 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 C$^{18}$O 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 C$_2$H column density traces changes in the C/O ratio rather than the H$_2$ gas surface density. However, the presence of a massive planet (> 0.2 M$_{Jup}$) 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 H$_2$ surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series. △ Less

Submitted 17 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 21 pages, 18 figures

Molecules with ALMA at Planet-forming Scales (MAPS) V: CO gas distributions

Authors: Ke Zhang, Alice S. Booth, Charles J. Law, Arthur D. Bosman, Kamber R. Schwarz, Edwin A. Bergin, Karin I. Öberg, Sean M. Andrews, Viviana V. Guzmán, Catherine Walsh, Chunhua Qi, Merel L. R. van 't Hoff, Feng Long, David J. Wilner, Jane Huang, Ian Czekala, John D. Ilee, Gianni Cataldi, Jennifer B. Bergner, Yuri Aikawa, Richard Teague, Jaehan Bae, Ryan A. Loomis, Jenny K. Calahan, Felipe Alarcón , et al. (10 additional authors not shown)

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 $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling f… ▽ More 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 $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-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 of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-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. △ Less

Submitted 23 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. 36 pages, 21 figures, accepted for publication in ApJS

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

Authors: Kamber R. Schwarz, Jenny K. Calahan, Ke Zhang, Felipe Alarcón, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Ryan A. Loomis, Enrique Macías, Melissa McClure, François Ménard, Karin I. Öberg, Richard Teague , et al. (3 additional authors not shown)

Abstract: 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 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk tempe… ▽ More 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 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ 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 isotopolog emission from eleven rotational transitions with spatial resolution ranging from 0.15'' to 0.46'' (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 solar masses, 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. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 18 pages, 12 figures

Molecules with ALMA at Planet-forming Scales (MAPS) XIX. Spiral Arms, a Tail, and Diffuse Structures Traced by CO around the GM Aur Disk

Authors: Jane Huang, Edwin A. Bergin, Karin I. Öberg, Sean M. Andrews, Richard Teague, Charles J. Law, Paul Kalas, Yuri Aikawa, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, John D. Ilee, Romane Le Gal, Viviana V. Guzmán, Feng Long, Ryan A. Loomis, François Ménard, Hideko Nomura, Chunhua Qi, Kamber R. Schwarz , et al. (6 additional authors not shown)

Abstract: The concentric gaps and rings commonly observed in protoplanetary disks in millimeter continuum emission have lent the impression that planet formation generally proceeds within orderly, isolated systems. While deep observations of spatially resolved molecular emission have been comparatively limited, they are increasingly suggesting that some disks interact with their surroundings while planet fo… ▽ More The concentric gaps and rings commonly observed in protoplanetary disks in millimeter continuum emission have lent the impression that planet formation generally proceeds within orderly, isolated systems. While deep observations of spatially resolved molecular emission have been comparatively limited, they are increasingly suggesting that some disks interact with their surroundings while planet formation is underway. We present an analysis of complex features identified around GM Aur in $^{12}$CO $J=2-1$ images at a spatial resolution of $\sim40$ au. In addition to a Keplerian disk extending to a radius of $\sim550$ au, the CO emission traces flocculent spiral arms out to radii of $\sim$1200 au, a tail extending $\sim1800$ au southwest of GM Aur, and diffuse structures extending from the north side of the disk up to radii of $\sim1900$ au. The diffuse structures coincide with a "dust ribbon" previously identified in scattered light. The large-scale asymmetric gas features present a striking contrast with the mostly axisymmetric, multi-ringed millimeter continuum tracing the pebble disk. We hypothesize that GM Aur's complex gas structures result from late infall of remnant envelope or cloud material onto the disk. The morphological similarities to the SU Aur and AB Aur systems, which are also located in the L1517 cloud, provide additional support to a scenario in which interactions with the environment are playing a role in regulating the distribution and transport of material in all three of these Class II disk systems. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 16 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 34 pages, 21 figures, in press at ApJS, cross-references updated

Molecules with ALMA at Planet-forming Scales (MAPS). XV. Tracing protoplanetary disk structure within 20 au

Authors: Arthur D. Bosman, Edwin A. Bergin, Ryan A. Loomis, Sean M. Andrews, Merel L. R. van 't Hoff, Richard Teague, Karin I. Öberg, Viviana V. Guzmán, Catherine Walsh, Yuri Aikawa, Felipe Alarcón, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Hideko Nomura , et al. (3 additional authors not shown)

Abstract: Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large p… ▽ More Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large program to construct radial surface brightness profiles with a ~3 au effective resolution for the CO isotopologue J=2-1 lines using the line velocity profile. IM Lup reveals a central depression in 13CO and C18O that is ascribed to a pileup of ~500 $M_\oplus$ of dust in the inner 20 au, leading to a gas-to-dust ratio of around <10. This pileup is consistent with efficient drift of grains ($\gtrsim$ 100 $M_\oplus$ Myr$^{-1}$) and a local gas-to-dust ratio that suggests that the streaming instability could be active. The CO isotopologue emission in the GM Aur disk is consistent with a small (~15 au), strongly depleted gas cavity within the ~40 au dust cavity. The radial surface brightness profiles for both the AS 209 and HD 163296 disks show a local minimum and maximum in the C18O emission at the location of a known dust ring (~14 au) and gap (~10 au), respectively. This indicates that the dust ring has a low gas-to-dust ratio ($>$ 10) and that the dust gap is gas-rich enough to have optically thick C18O. △ Less

Submitted 16 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 22 pages, 14 figures, accepted by ApJS. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Update with correct references to other MAPS papers

Molecules with ALMA at Planet-forming Scales (MAPS). VII. Sub-stellar O/H and C/H and super-stellar C/O in planet feeding gas

Authors: Arthur D. Bosman, Felipe Alarcón, Edwin A. Bergin, Ke Zhang, Merel L. R. van 't Hoff, Karin I. Öberg, Viviana V. Guzmán, Catherine Walsh, Yuri Aikawa, Sean M. Andrews, Jennifer B. Bergner, Alice S. Booth, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Kenji Furuya, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, Ryan A. Loomis, François Ménard, Hideko Nomura , et al. (6 additional authors not shown)

Abstract: The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models, wi… ▽ More The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models, with variable C/O ratios and small-grain abundances, from the DALI code with CO and C2H column densities derived from the high-resolution observations of the disks of AS 209, HD 163296, and MWC 480. Elevated C/O ratios (~2.0), even within the CO ice line, are necessary to match the inferred C2H column densities, over most of the pebble disk. Combined with constraints on the CO abundances in these systems, this implies that both the O/H and C/H ratios in the gas are substellar by a factor of 4-10, with the O/H depleted by a factor of 20-50, resulting in the high C/O ratios. This necessitates that even within the CO ice line, most of the volatile carbon and oxygen is still trapped on grains in the midplane. Planets accreting gas in the gaps of the AS 209, HD 163296, and MWC 480 disks will thus acquire very little carbon and oxygen after reaching the pebble isolation mass. In the absence of atmosphere-enriching events, these planets would thus have a strongly substellar O/H and C/H and superstellar C/O atmospheric composition. △ Less

Submitted 16 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 19 pages, 8 figures This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Updates references for other MAPS papers

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

Authors: Richard Teague, Jaehan Bae, Yuri Aikawa, Sean M. Andrews, Edwin A. Bergin, Jennifer B. Bergner, Yann Boehler, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, Ian Czekala, Viviana V. Guzmán, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Feng Long, Ryan A. Loomis, François Ménard, Karin I. Öberg, Laura M. Pérez, Kamber R. Schwarz, Anibal Sierra, Catherine Walsh, David J. Wilner , et al. (2 additional authors not shown)

Abstract: 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 $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Kepleria… ▽ More 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 $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Keplerian rotation in the projected velocity fields of both disks ($\lesssim\!5\%$ of the local Keplerian velocity), suggestive of large-scale (10s of au 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_φ$, $v_r$ and $v_z$. Using the optically thick $^{12}$CO emission as a probe of the gas temperature, local variations of $\approx\! 3$ K ($\approx\! 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, $\sim\! 1$ $M_{\rm Jup}$, at a radius of $\approx 245$ au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. (2018a) and Teague et al. (2018a) advocating for multiple $\sim\! 1$ $M_{\rm Jup}$ 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. △ Less

Submitted 20 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 32 pages, 21 figures, accepted for publication in ApJS, MAPS cross-references updated

Molecules with ALMA at Planet-forming Scales (MAPS) IV: Emission Surfaces and Vertical Distribution of Molecules

Authors: Charles J. Law, Richard Teague, Ryan A. Loomis, Jaehan Bae, Karin I. Öberg, Ian Czekala, Sean M. Andrews, Yuri Aikawa, Felipe Alarcón, Edwin A. Bergin, Jennifer B. Bergner, Alice S. Booth, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Kenji Furuya, Viviana V. Guzmán, Jane Huang, John D. Ilee, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Hideko Nomura , et al. (10 additional authors not shown)

Abstract: The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a funct… ▽ More The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a function of radius ($r$). Using this method, we measure emitting surfaces for a suite of CO isotopologues, HCN, and C$_2$H. We find that $^{12}$CO emission traces the most elevated regions with $z/r > 0.3$, while emission from the less abundant $^{13}$CO and C$^{18}$O probes deeper into the disk at altitudes of $z/r \lesssim 0.2$. C$_2$H and HCN have lower opacities and SNRs, making surface fitting more difficult, and could only be reliably constrained in AS 209, HD 163296, and MWC 480, with $z/r \lesssim 0.1$, i.e., relatively close to the planet-forming midplanes. We determine peak brightness temperatures of the optically thick CO isotopologues and use these to trace 2D disk temperature structures. Several CO temperature profiles and emission surfaces show dips in temperature or vertical height, some of which are associated with gaps and rings in line and/or continuum emission. These substructures may be due to local changes in CO column density, gas surface density, or gas temperatures, and detailed thermo-chemical models are necessary to better constrain their origins and relate the chemical compositions of elevated disk layers with those of planet-forming material in disk midplanes. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 20 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 31 pages, 20 figures, accepted for publication in ApJS, MAPS cross-references updated

Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

Authors: Charles J. Law, Ryan A. Loomis, Richard Teague, Karin I. Öberg, Ian Czekala, Sean M. Andrews, Jane Huang, Yuri Aikawa, Felipe Alarcón, Jaehan Bae, Edwin A. Bergin, Jennifer B. Bergner, Yann Boehler, Alice S. Booth, Arthur D. Bosman, Jenny K. Calahan, Gianni Cataldi, L. Ilsedore Cleeves, Kenji Furuya, Viviana V. Guzmán, John D. Ilee, Romane Le Gal, Yao Liu, Feng Long, François Ménard , et al. (10 additional authors not shown)

Abstract: The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We i… ▽ More The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 May, 2022; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 62 pages, 31 figures, accepted for publication in ApJS, MAPS cross-references updated, corrected Figure 21, updated gas disk sizes (Table 2, Figures 15-16) from associated Erratum

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

Authors: Jenny K. Calahan, Edwin A. Bergin, Ke Zhang, Kamber R. Schwarz, Karin I. Oberg, Viviana V. Guzman, Catherine Walsh, Yuri Aikawa, Felipe Alarcon, Sean M. Andrews, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Arthur D. Bosman, Gianni Cataldi, Ian Czekala, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Feng Long, Ryan A. Loomis, Francois Menard, Hideko Nomura, Chunhua Qi , et al. (4 additional authors not shown)

Abstract: 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 Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spa… ▽ More 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 Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spatially resolved ALMA observations (~0.12 arcsec (13 au) - 0.25 arcsec (26 au)) of CO 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 time scale 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 Msun, 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. △ Less

Submitted 24 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: 15 pages + 11 pages of appendix, accepted to ApJS, part of MAPS collaboration

Molecules with ALMA at Planet-forming Scales (MAPS) II: CLEAN Strategies for Synthesizing Images of Molecular Line Emission in Protoplanetary Disks

Authors: Ian Czekala, Ryan A. Loomis, Richard Teague, Alice S. Booth, Jane Huang, Gianni Cataldi, John D. Ilee, Charles J. Law, Catherine Walsh, Arthur D. Bosman, Viviana V. Guzmán, Romane Le Gal, Karin I. Öberg, Yoshihide Yamato, Yuri Aikawa, Sean M. Andrews, Jaehan Bae, Edwin A. Bergin, Jennifer B. Bergner, L. Ilsedore Cleeves, Nicolas T. Kurtovic, François Ménard, Hideko Nomura, Laura M. Pérez, Chunhua Qi , et al. (5 additional authors not shown)

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 and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our mu… ▽ More 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 and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our multi-stage 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 sub-mm protoplanetary disk community, like the benefits of using an accurate CLEAN mask, and others less well-known, like the incorrect default flux scaling of the CLEAN residual map first described in Jorsater & van Moorsel 1995 (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 the "JvM correction," a procedure that is generally applicable to any image synthesized via CLEAN deconvolution but is especially critical for low S/N emission. We further explain how we used visibility tapering to promote a common, fiducial beam size and contextualize the interpretation of signal to noise ratio when detecting molecular emission from protoplanetary disks. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 24 September, 2021; v1 submitted 13 September, 2021; originally announced September 2021.

Comments: ApJS accepted. Part of the MAPS ALMA large program series: http://www.alma-maps.info/. Updated bibliography with MAPS LP arXiv references

Detection of Interstellar H$_2$CCCHC$_3$N

Authors: C. N. Shingledecker, K. L. K. Lee, J. T. Wandishin, N. Balucani, A. M. Burkhardt, S. B. Charnley, R. Loomis, M. Schreffler, M. Siebert, M. C. McCarthy, B. A. McGuire

Abstract: The chemical pathways linking the small organic molecules commonly observed in molecular clouds to the large, complex, polycyclic species long-suspected to be carriers of the ubiquitous unidentified infrared emission bands remain unclear. To investigate whether the formation of mono- and poly-cyclic molecules observed in cold cores could form via the bottom-up reaction of ubiquitous carbon-chain s… ▽ More The chemical pathways linking the small organic molecules commonly observed in molecular clouds to the large, complex, polycyclic species long-suspected to be carriers of the ubiquitous unidentified infrared emission bands remain unclear. To investigate whether the formation of mono- and poly-cyclic molecules observed in cold cores could form via the bottom-up reaction of ubiquitous carbon-chain species with, e.g. atomic hydrogen, a search is made for possible intermediates in data taken as part of the GOTHAM (GBT Observations of TMC-1 Hunting for Aromatic Molecules) project. Markov-Chain Monte Carlo (MCMC) Source Models were run to obtain column densities and excitation temperatures. Astrochemical models were run to examine possible formation routes, including a novel grain-surface pathway involving the hydrogenation of C$_6$N and HC$_6$N, as well as purely gas-phase reactions between C$_3$N and both propyne (CH$_3$CCH) and allene (CH$_2$CCH$_2$), as well as via the reaction CN + H$_2$CCCHCCH. We report the first detection of cyanoacetyleneallene (H$_2$CCCHC$_3$N) in space toward the TMC-1 cold cloud using the Robert C. Byrd 100 m Green Bank Telescope (GBT). Cyanoacetyleneallene may represent an intermediate between less-saturated carbon-chains, such as the cyanopolyynes, that are characteristic of cold cores and the more recently-discovered cyclic species like cyanocyclopentadiene. Results from our models show that the gas-phase allene-based formation route in particular produces abundances of H$_2$CCCHC$_3$N that match the column density of $2\times10^{11}$ cm$^{-2}$ obtained from the MCMC Source Model, and that the grain-surface route yields large abundances on ices that could potentially be important as precursors for cyclic molecules. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: Accepted in A&A Letters

Discovery of the Pure Polycyclic Aromatic Hydrocarbon Indene ($c$-C$_9$H$_8$) with GOTHAM Observations of TMC-1

Authors: Andrew M. Burkhardt, Kin Long Kelvin Lee, P. Bryan Changala, Christopher N. Shingledecker, Ilsa R. Cooke, Ryan A. Loomis, Hongji Wei, Steven B. Charnley, Eric Herbst, Michael C. McCarthy, Brett A. McGuire

Abstract: Polycyclic Aromatic Hydrocarbons (PAHs) have long been invoked in the study of interstellar and protostellar sources, but the unambiguous identification of any individual PAH has proven elusive until very recently. As a result, the formation mechanisms for this important class of molecules remain poorly constrained. Here we report the first interstellar detection of a pure hydrocarbon PAH, indene… ▽ More Polycyclic Aromatic Hydrocarbons (PAHs) have long been invoked in the study of interstellar and protostellar sources, but the unambiguous identification of any individual PAH has proven elusive until very recently. As a result, the formation mechanisms for this important class of molecules remain poorly constrained. Here we report the first interstellar detection of a pure hydrocarbon PAH, indene (C$_9$H$_8$), as part of the GBT Observations of TMC-1: Hunting for Aromatic Molecules (GOTHAM) survey. This detection provides a new avenue for chemical inquiry, complementing the existing detections of CN-functionalized aromatic molecules. From fitting the GOTHAM observations, indene is found to be the most abundant organic ring detected in TMC-1 to date. And from astrochemical modeling with NAUTILUS, the observed abundance is greater than the model's prediction by several orders of magnitude suggesting that current formation pathways in astrochemical models are incomplete. The detection of indene in relatively high abundance implies related species such as cyanoindene, cyclopentadiene, toluene, and styrene may be detectable in dark clouds. △ Less

Submitted 27 May, 2021; v1 submitted 30 April, 2021; originally announced April 2021.

Comments: 18 pages, 2 figures and 3 tables in the main text. 1 table and 1 figure in the Appendix. Published in The Astrophysical Journal Letters. Supplementary data available in the DataVerse entry provided in text

Journal ref: ApJL 2021, 913 L18

Detection of Two Interstellar Polycyclic Aromatic Hydrocarbons via Spectral Matched Filtering

Authors: Brett A. McGuire, Ryan A. Loomis, Andrew M. Burkhardt, Kin Long Kelvin Lee, Christopher N. Shingledecker, Steven B. Charnely, Ilsa R. Cooke, Martin A. Cordiner, Eric Herbst, Sergei Kalenskii, Mark A. Siebert, Eric R. Willis, Ci Xue, Anthony J. Remijan, Michael C. McCarthy

Abstract: Ubiquitous unidentified infrared emission bands are seen in many astronomical sources. Although these bands are widely, if not unanimously, attributed to the collective emission from polycyclic aromatic hydrocarbons, no single species from this class has been detected in space. We present the discovery of two -CN functionalized polycyclic aromatic hydrocarbons, 1- and 2-cyanonaphthalene, in the in… ▽ More Ubiquitous unidentified infrared emission bands are seen in many astronomical sources. Although these bands are widely, if not unanimously, attributed to the collective emission from polycyclic aromatic hydrocarbons, no single species from this class has been detected in space. We present the discovery of two -CN functionalized polycyclic aromatic hydrocarbons, 1- and 2-cyanonaphthalene, in the interstellar medium aided by spectral matched filtering. Using radio observations with the Green Bank Telescope, we observe both bi-cyclic ring molecules in the molecular cloud TMC-1. We discuss potential in situ gas-phase formation pathways from smaller organic precursor molecules. △ Less

Submitted 17 March, 2021; originally announced March 2021.

Journal ref: Science (2021) vol 371, pg 1265-1269

Interstellar Detection of 2-Cyanocyclopentadiene, C$_5$H$_5$CN, a Second Five-Membered Ring Toward TMC-1

Authors: Kin Long Kelvin Lee, P. Bryan Changala, Ryan A. Loomis, Andrew M. Burkhardt, Ci Xue, Martin A. Cordiner, Steven B. Charnley, Michael C. McCarthy, Brett A. McGuire

Abstract: Using radio observations with the Green Bank Telescope, evidence has now been found for a second five-membered ring in the dense cloud Taurus Molecular Cloud-1 (TMC-1). Based on additional observations of an ongoing, large-scale, high-sensitivity spectral line survey (GOTHAM) at centimeter wavelengths toward this source, we have used a combination of spectral stacking, Markov chain Monte Carlo (MC… ▽ More Using radio observations with the Green Bank Telescope, evidence has now been found for a second five-membered ring in the dense cloud Taurus Molecular Cloud-1 (TMC-1). Based on additional observations of an ongoing, large-scale, high-sensitivity spectral line survey (GOTHAM) at centimeter wavelengths toward this source, we have used a combination of spectral stacking, Markov chain Monte Carlo (MCMC), and matched filtering techniques to detect 2-cyanocyclopentadiene, a low-lying isomer of 1-cyanocyclopentadiene, which was recently discovered there by the same methods. The new observational data also yields a considerably improved detection significance for the more stable isomer and evidence for several individual transitions between 23 - 32 GHz. Through our MCMC analysis, we derive total column densities of $8.3\times10^{11}$ and $1.9\times10^{11}$ cm$^{-2}$ for 1- and 2-cyanocyclopentadiene respectively, corresponding to a ratio of 4.4(6) favoring the former. The derived abundance ratios point towards a common formation pathway - most likely being cyanation of cyclopentadiene by analogy to benzonitrile. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 6 pages, 3 figures and 1 table in the main text. 2 tables and 2 figures in the Appendix. Accepted for publication in The Astrophysical Journal Letters. Supplementary data available in the DataVerse entry provided in text

The TW Hya Rosetta Stone Project IV: A hydrocarbon rich disk atmosphere

Authors: L. Ilsedore Cleeves, Ryan A. Loomis, Richard Teague, Edwin A. Bergin, David J. Wilner, Jennifer B. Bergner, Geoffrey A. Blake, Jenny K. Calahan, Paolo Cazzoletti, Ewine F. van Dishoeck, Viviana V. Guzman, Michiel R. Hogerheijde, Jane Huang, Mihkel Kama, Karin I. Oberg, Chunhua Qi, Jeroen Terwisscha van Scheltinga, Catherine Walsh

Abstract: 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 $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project… ▽ More 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 $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ 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, $E_u=29-97$ K. We find the ortho-to-para ratio of $\rm C_3H_2$ is consistent with 3 throughout extent of the emission, and the total abundance of both $\rm C_3H_2$ isomers is $(7.5-10)\times10^{-11}$ per H atom, or $1-10$% of the previously published $\rm C_2H$ abundance in the same source. We find $\rm C_3H_2$ 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 $\sim30$ au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridonial 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. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures, Accepted in ApJ

Discovery of Interstellar trans-cyanovinylacetylene (HCCCH=CHCN) and vinylcyanoacetylene (H$_2$C=CHC$_3$N) in GOTHAM Observations of TMC-1

Authors: Kin Long Kelvin Lee, Ryan A. Loomis, Andrew M. Burkhardt, Ilsa R. Cooke, Ci Xue, Mark A. Siebert, Christopher N. Shingledecker, Anthony Remijan, Steven B. Charnley, Michael C. McCarthy, Brett A. McGuire

Abstract: We report the discovery of two unsaturated organic species, trans-(E)-cyanovinylacetylene and vinylcyanoacetylene, using the second data release of the GOTHAM deep survey towards TMC-1 with the 100 m Green Bank Telescope. For both detections, we performed velocity stacking and matched filter analyses using Markov chain Monte Carlo simulations, and for trans-(E)-cyanovinylacetylene, three rotationa… ▽ More We report the discovery of two unsaturated organic species, trans-(E)-cyanovinylacetylene and vinylcyanoacetylene, using the second data release of the GOTHAM deep survey towards TMC-1 with the 100 m Green Bank Telescope. For both detections, we performed velocity stacking and matched filter analyses using Markov chain Monte Carlo simulations, and for trans-(E)-cyanovinylacetylene, three rotational lines were observed at low signal-to-noise (${\sim}$3$σ$). From this analysis, we derive column densities of $2\times10^{11}$ and $3\times10^{11}$ cm$^{-2}$ for vinylcyanoacetylene and trans-(E)-cyanovinylacetylene, respectively, and an upper limit of $<2\times10^{11}$ cm$^{-2}$ for trans-(Z)-cyanovinylacetylene. Comparisons with G3//B3LYP semi-empirical thermochemical calculations indicate abundances of the [H$_3$C$_5$N}] isomers are not consistent with their thermodynamic stability, and instead their abundances are mainly driven by dynamics. We provide discussion into how these species may be formed in TMC-1, with reference to related species like vinyl cyanide (CH$_2$=CHCN). As part of this discussion, we performed the same analysis for ethyl cyanide (CH$_3$CH$_2$CN), the hydrogenation product of CH$_2$=CHCN. This analysis provides evidence -- at 4.17$σ$ significance -- an upper limit to the column density of $<4\times10^{11}$ cm$^{-2}$; an order of magnitude lower than previous upper limits towards this source. △ Less

Submitted 12 February, 2021; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: Four figures, two tables. Five supplementary figures and four supplementary tables. Published in The Astrophysical Journal Letters

Journal ref: ApJL 2021, 908 L11

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

Authors: Jenny Calahan, Edwin Bergin, Ke Zhang, Richard Teague, Ilsedore Cleeves, Jennifer Bergner, Geoffrey A. Blake, Paolo Cazzoletti, Viviana Guzman, Michiel R. Hogerheijde, Jane Huang, Mihkel Kama, Ryan Loomis, Karin Oberg, Ewine F. van Dishoeck, Jeroen Terwisscha van Scheltinga, Catherine Walsh, David Wilner, Charlie Qi

Abstract: 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 o… ▽ More 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, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, 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 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 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 supports 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. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 28 pages, 12 main figures, 14 appendix figures. Accepted to ApJ

The TW Hya Rosetta Stone Project II: Spatially resolved emission of formaldehyde hints at low-temperature gas-phase formation

Authors: Jeroen Terwisscha van Scheltinga, Michiel R. Hogerheijde, L. Ilsedore Cleeves, Ryan A. Loomis, Catherine Walsh, Karin I. Öberg, Edwin A. Bergin, Jennifer B. Bergner, Geoffrey A. Blake, Jenny K. Calahan, Paolo Cazzoletti, Ewine F. van Dishoeck, Viviana V. Guzmán, Jane Huang, Mihkel Kama, Chunhua Qi, Richard Teague, David J. Wilner

Abstract: Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pa… ▽ More Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pathway dominates, gaseous H$_2$CO'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 ALMA observations of multiple ortho- and para-H$_2$CO transitions in the TW Hya protoplanetary disk to test H$_2$CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of $33\pm2$ K, ($1.1\pm0.1)\times10^{12}$ cm$^{-2}$ and $25\pm2$ K, $(4.4\pm0.3)\times10^{11}$ cm$^{-2}$ for ortho- and para-H$_2$CO, respectively, and an OPR of $2.49\pm0.23$. A radially resolved analysis shows that the observed H$_2$CO emits mostly at rotational temperatures of 30-40 K, corresponding to a layer with $z/R\ge0.25$. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to $2.0\pm0.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, lead us to speculate that gas-phase formation is responsible for the observed H$_2$CO across the TW Hya disk. △ Less

Submitted 13 November, 2020; originally announced November 2020.

The TW Hya Rosetta Stone Project I: Radial and vertical distributions of DCN and DCO+

Authors: Karin I. Oberg, L. Ilsedore Cleeves, Jennifer B. Bergner, Joseph Cavanaro, Richard Teague, Jane Huang, Ryan A. Loomis, Edwin A. Bergin, Geoffrey A. Blake, Jenny Calahan, Paolo Cazzoletti, Viviana Veloso Guzman, Michiel R. Hogerheijde, Mihkel Kama, Jeroen Terwisscha van Scheltinga, Chunhua Qi, Ewine van Dishoeck, Catherine Walsh, David J. Wilner

Abstract: 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 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project,… ▽ More 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 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards 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 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. There appears to be 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. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: Accepted for publication in The Astronomical Journal

Ubiquitous Aromatic Carbon Chemistry at the Earliest Stages of Star Formation

Authors: Andrew M. Burkhardt, Ryan A. Loomis, Christopher N. Shingledecker, Kin Long Kelvin Lee, Anthony J. Remijan, Michael C. McCarthy, Brett A. McGuire

Abstract: Benzonitrile ($c$-C$_6$H$_5$CN), a polar proxy for benzene ($c$-C$_6$H$_6$}), has the potential to serve as a highly convenient radio probe for aromatic chemistry, provided this ring can be found in other astronomical sources beyond the molecule-rich prestellar cloud TMC-1 where it was first reported by McGuire et al. in 2018. Here we present radio astronomical evidence of benzonitrile in four add… ▽ More Benzonitrile ($c$-C$_6$H$_5$CN), a polar proxy for benzene ($c$-C$_6$H$_6$}), has the potential to serve as a highly convenient radio probe for aromatic chemistry, provided this ring can be found in other astronomical sources beyond the molecule-rich prestellar cloud TMC-1 where it was first reported by McGuire et al. in 2018. Here we present radio astronomical evidence of benzonitrile in four additional pre-stellar, and possibly protostellar, sources: Serpens 1A, Serpens 1B, Serpens 2, and MC27/L1521F. These detections establish benzonitrile is not unique to TMC-1; rather aromatic chemistry appears to be widespread throughout the earliest stages of star formation, likely persisting at least to the initial formation of a protostar. The abundance of benzonitrile far exceeds predictions from models which well reproduce the abundances of carbon chains, such as HC$_7$N, a cyanpolyyne with the same heavy atoms, indicating the chemistry responsible for planar carbon structures (as opposed to linear ones) in primordial sources is favorable but not well understood. The abundance of benzonitrile relative to carbon-chain molecules displays sizable variations between sources within the Taurus and Serpens clouds, implying the importance of physical conditions and initial elemental reservoirs of the clouds themselves. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: To Appear in Nature Astronomy, 53 Pages, 15 Figures, 12 Tables