Showing 1–50 of 133 results for author: Walsh, C

How Land-Mass Distribution Influences the Atmospheric Dynamics of Tidally Locked Terrestrial Exoplanets

Authors: F. Sainsbury-Martinez, C. Walsh, G. J. Cooke, D. R. Marsh

Abstract: Interpretation of the ongoing efforts to simulate the atmospheres of potentially-habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally-locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terres… ▽ More Interpretation of the ongoing efforts to simulate the atmospheres of potentially-habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally-locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terrestrial planets: the inclusion of an Earth-like land-ocean distribution with orography. To do this we use the Earth-system model WACCM6/CESM2 to run a pair of TRAPPIST-1e models with N$_2$-O$_2$ atmospheres and with the sub-stellar point fixed over either land or ocean. The presence of orography shapes atmospheric transport, and in the case of Earth-like orography, breaks the symmetry between the northern and southern hemispheres which was previously found in slab ocean models. For example, peak zonal jet speeds in the southern hemisphere are $50\rightarrow100\%$ faster than similar jets in the northern hemisphere. This also affects the meridional circulation, transporting equatorial material towards the south-pole. As a result we also find significant changes in the atmospheric chemistry, including the accumulation of potentially lethal quantities of ozone at both the south pole and the surface. Future studies which investigate the effects of land-mass distribution on the dynamics of exoplanetary atmospheres should pay close attention to both the day-side land-fraction as well as the orography of the land. Simply modelling a flat land-mass will not give a complete picture of its dynamical impact. △ Less

Submitted 2 July, 2024; v1 submitted 1 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ

Hybrid approach predicts a lower binding energy for benzene on water ice

Authors: Victoria H. J. Clark, David M. Benoit, Marie Van de Sande, Catherine Walsh

Abstract: In this paper we provide a highly accurate value for the binding energy of benzene to proton-ordered crystalline water ice (XIh), as a model for interstellar ices. We compare our computed value to the latest experimental data available from temperature programmed desorption (TPD) experiments and find that our binding energy value agrees well with data obtained from binding to either crystalline or… ▽ More In this paper we provide a highly accurate value for the binding energy of benzene to proton-ordered crystalline water ice (XIh), as a model for interstellar ices. We compare our computed value to the latest experimental data available from temperature programmed desorption (TPD) experiments and find that our binding energy value agrees well with data obtained from binding to either crystalline or amorphous ice. Importantly, our new value is lower than that used in most astrochemical networks by about nearly half its value. We explore the impact of this revised binding energy value for both an AGB outflow and a protoplanetary disk. We find that the lower value of the binding energy predicted here compared with values used in the literature (4050 K versus 7587 K) leads to less depletion of gas-phase benzene in an AGB outflow, and leads to a shift outwards in the benzene snowline in the midplane of a protoplanetary disk. Using this new value, the AGB model predicts lower abundances of benzene in the solid phase throughout the outflow. The disk model also predicts a larger reservoir of gas-phase benzene in the inner disk, which is consistent with the recent detections of benzene for the first time in protoplanetary disks with JWST. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Modelling predicts a molecule-rich disk around the AGB star L2 Puppis

Authors: M. Van de Sande, C. Walsh, T. Danilovich, F. De Ceuster, T. Ceulemans

Abstract: The nearby oxygen-rich AGB star L2 Pup hosts a well-studied nearly edge-on disk. To date, disks around AGB stars have not been chemically studied in detail. By combining a parameterisation commonly used for protoplanetary disks and archival ALMA observations, we retrieved an updated density and temperature structure of this disk. This physical model was then used as input to the first chemical mod… ▽ More The nearby oxygen-rich AGB star L2 Pup hosts a well-studied nearly edge-on disk. To date, disks around AGB stars have not been chemically studied in detail. By combining a parameterisation commonly used for protoplanetary disks and archival ALMA observations, we retrieved an updated density and temperature structure of this disk. This physical model was then used as input to the first chemical model of an AGB disk. The model shows that the physical structure of the disk has a large impact on its chemistry, with certain species showing large changes in column density relative to a radial outflow, indicating that chemistry could be used as a tracer of disks that cannot be directly imaged. Despite its oxygen-rich nature, the daughter species formed within the disk are surprisingly carbon-rich. Two chemical regimes can be distinguished: cosmic-ray induced chemistry in the midplane and photochemistry induced by the interstellar radiation field in the outer regions. Certain complex organic molecules are formed in the midplane. This occurs via gas-phase chemistry only, as the disk is too warm for dust-gas chemistry. The photochemistry in the outer regions leads to the efficient formation of (long) carbon-chains. The predictions of the model allow us to tentatively put the disk's age $\lesssim 10^5$ yr. Additional observations are necessary to better constrain the physical structure of L2 Pup's disk and are essential to test the predictions made by the chemical model. Our exploratory work paves the way for a more general study of the chemistry of AGB disks. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: Accepted for publication in MNRAS

Lethal surface ozone concentrations are possible on habitable zone exoplanets

Authors: G. J. Cooke, D. R. Marsh, C. Walsh, F. Sainsbury-Martinez

Abstract: Ozone ($\textrm{O}_3$) is important for the survival of life on Earth because it shields the surface from ionising ultraviolet (UV) radiation. However, the existence of $\textrm{O}_3$ in Earth's atmosphere is not always beneficial. Resulting from anthropogenic activity, $\textrm{O}_3$ exists as a biologically harmful pollutant at the surface when it forms in the presence of sunlight and other poll… ▽ More Ozone ($\textrm{O}_3$) is important for the survival of life on Earth because it shields the surface from ionising ultraviolet (UV) radiation. However, the existence of $\textrm{O}_3$ in Earth's atmosphere is not always beneficial. Resulting from anthropogenic activity, $\textrm{O}_3$ exists as a biologically harmful pollutant at the surface when it forms in the presence of sunlight and other pollutants. As a strong oxidiser, $\textrm{O}_3$ can be lethal to several different organisms; thus, when assessing the potential habitability of an exoplanet, a key part is determining whether toxic gases could be present at its surface. Using the Whole Atmosphere Community Climate Model version 6 (WACCM6; a three-dimensional chemistry-climate model), twelve atmospheric simulations of the terrestrial exoplanet TRAPPIST-1 e are performed with a variety of $\textrm{O}_2$ concentrations and assuming two different stellar spectra proposed in the literature. Four atmospheric simulations of the exoplanet Proxima Centauri b are also included. Some scenarios for both exoplanets exhibit time-averaged surface $\textrm{O}_3$ mixing ratios exceeding harmful levels of 40 ppbv, with 2200 ppbv the maximum concentration found in the cases simulated. These concentrations are toxic and can be fatal to most life on Earth. In other scenarios $\textrm{O}_3$ remains under harmful limits over a significant fraction of the surface, despite there being present regions which may prove inhospitable. In the case that $\textrm{O}_3$ is detected in a terrestrial exoplanet's atmosphere, determining the surface concentration is an important step when evaluating a planet's habitability. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 24 pages, 8 figures. Accepted in The Planetary Science Journal

Chemistry in externally FUV irradiated disks in the outskirts of the Orion Nebula

Authors: Javiera K. Díaz-Berríos, Viviana V. Guzmán, Catherine Walsh, Karin I. Öberg, L. Ilsedore Cleeves, Elizabeth Artur de la Villarmois, John Carpenter

Abstract: Most stars are born in stellar clusters and their protoplanetary disks, which are the birthplaces of planets, can therefore be affected by the radiation of nearby massive stars. However, little is known about the chemistry of externally irradiated disks, including whether or not their properties are similar to the so-far better-studied isolated disks. Motivated by this question, we present ALMA Ba… ▽ More Most stars are born in stellar clusters and their protoplanetary disks, which are the birthplaces of planets, can therefore be affected by the radiation of nearby massive stars. However, little is known about the chemistry of externally irradiated disks, including whether or not their properties are similar to the so-far better-studied isolated disks. Motivated by this question, we present ALMA Band 6 observations of two irradiated Class II protoplanetary disks in the outskirts of the Orion Nebula Cluster (ONC) to explore the chemical composition of disks exposed to (external) FUV radiation fields: the 216-0939 disk and the binary system 253-1536A/B, which are exposed to radiation fields of $10^2-10^3$ times the average interstellar radiation field. We detect lines from CO isotopologues, HCN, H$_2$CO, and C$_2$H toward both protoplanetary disks. Based on the observed disk-integrated line fluxes and flux ratios, we do not find significant differences between isolated and irradiated disks. The observed differences seem to be more closely related to the different stellar masses than to the external radiation field. This suggests that these disks are far enough away from the massive Trapezium stars, that their chemistry is no longer affected by external FUV radiation. Additional observations towards lower-mass disks and disks closer to the massive Trapezium stars are required to elucidate the level of external radiation required to make an impact on the chemistry of planet formation in different kinds of disks. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 32 pages, 22 figures, Accepted for publication in ApJ

The Impact of Cometary 'impacts' on the Chemistry, Climate, and Spectra of Hot Jupiter Atmospheres

Authors: Felix Sainsbury-Martinez, Catherine Walsh

Abstract: Impacts from icy and rocky bodies have helped shape the composition of solar system objects, for example the Earth-Moon system, or the recent impact of comet Shoemaker-Levy 9 with Jupiter. It is likely that such impacts also shape the composition of exoplanetary systems. Here we investigate how cometary impacts might affect the atmospheric composition/chemistry of hot Jupiters, which are prime tar… ▽ More Impacts from icy and rocky bodies have helped shape the composition of solar system objects, for example the Earth-Moon system, or the recent impact of comet Shoemaker-Levy 9 with Jupiter. It is likely that such impacts also shape the composition of exoplanetary systems. Here we investigate how cometary impacts might affect the atmospheric composition/chemistry of hot Jupiters, which are prime targets for characterisation. We introduce a parametrised cometary impact model that includes thermal ablation and pressure driven breakup, which we couple with the 1D `radiative-convective' atmospheric model ATMO, including disequilibrium chemistry. We use this model to investigate a wide range of impactor masses and compositions, including those based on observations of Solar System comets, and interstellar ices (with JWST). We find that even a small impactor (R = 2.5 km) can lead to significant short-term changes in the atmospheric chemistry, including a factor $>10$ enhancement in H$_2$O, CO, CO$_2$ abundances, and atmospheric opacity more generally, and the near complete removal of observable hydrocarbons, such as CH$_4$, from the upper atmosphere. These effects scale with the change in atmospheric C/O ratio and metallicity. Potentially observable changes are possible for a body that has undergone significant/continuous bombardment, such that the global atmospheric chemistry has been impacted. Our works reveals that cometary impacts can significantly alter or pollute the atmospheric composition/chemistry of hot Jupiters. These changes have the potential to mute/break the proposed link between atmospheric C/O ratio and planet formation location relative to key snowlines in the natal protoplanetary disc. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ

An ALMA molecular inventory of warm Herbig Ae disks: II. Abundant complex organics and volatile sulphur in the IRS 48 disk

Authors: Alice S. Booth, Milou Temmink, Ewine F. van Dishoeck, Lucy Evans, John D. Ilee, Mihkel Kama, Luke Keyte, Charles J. Law, Margot Leemker, Nienke van der Marel, Hideko Nomura, Shota Notsu, Karin Öberg, Catherine Walsh

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in the… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in these disks. The IRS 48 disk shows a particularly rich chemistry that can be directly linked to its asymmetric dust trap. Here, we present ALMA observations of the IRS 48 disk where we detect 16 different molecules and make the first robust detections of H213CO, 34SO, 33SO and c-H2COCH2 (ethylene oxide) in a protoplanetary disk. All of the molecular emissions, aside from CO, are colocated with the dust trap and this includes newly detected simple molecules such as HCO+, HCN and CS. Interestingly, there are spatial offsets between different molecular families, including between the COMs and sulphur-bearing species, with the latter being more azimuthally extended and located radially further from the star. The abundances of the newly detected COMs relative to CH3OH are higher than the expected protostellar ratios, which implies some degree of chemical processing of the inherited ices during the disk lifetime. These data highlight IRS 48 as a unique astrochemical laboratory to unravel the full volatile reservoir at the epoch of planet and comet formation and the role of the disk in (re)setting chemical complexity. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted to AJ, 21 pages, 7 figures

An ALMA molecular inventory of warm Herbig Ae disks: I. Molecular rings, asymmetries and complexity in the HD 100546 disk

Authors: Alice S. Booth, Margot Leemker, Ewine F. van Dishoeck, Lucy Evans, John D. Ilee, Mihkel Kama, Luke Keyte, Charles J. Law, Nienke van der Marel, Hideko Nomura, Shota Notsu, Karin Öberg, Milou Temmink, Catherine Walsh

Abstract: Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three t… ▽ More Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three transition disks around Herbig Ae stars all show oxygen-rich gas traced via the unique detections of the molecules SO and CH3OH. We present the first results of an ALMA line survey at 337 to 357 GHz of such disks and focus this paper on the first Herbig Ae disk to exhibit this chemical signature - HD 100546. In these data, we detect 19 different molecules including NO, SO and CH3OCHO (methyl formate). We also make the first tentative detections of H213CO and 34SO in protoplanetary disks. Multiple molecular species are detected in rings, which are, surprisingly, all peaking just beyond the underlying millimeter continuum ring at 200 au. This result demonstrates a clear connection between the large dust distribution and the chemistry in this flat outer disk. We discuss the physical and/or chemical origin of these sub-structures in relation to ongoing planet formation in the HD 100546 disk. We also investigate how similar and/or different the molecular make up of this disk is to other chemically well-characterised Herbig Ae disks. The line-rich data we present motivates the need for more ALMA line surveys to probe the observable chemistry in Herbig Ae systems which offer unique insight into the composition of disk ices, including complex organic molecules. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted to AJ, 25 pages, 11 figures

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

The UMIST Database for Astrochemistry 2022

Authors: T J Millar, C Walsh, M Van de Sande, A J Markwick

Abstract: Detailed astrochemical models are a key component to interpret the observations of interstellar and circumstellar molecules since they allow important physical properties of the gas and its evolutionary history to be deduced. We update one of the most widely used astrochemical databases to reflect advances in experimental and theoretical estimates of rate coefficients and to respond to the large i… ▽ More Detailed astrochemical models are a key component to interpret the observations of interstellar and circumstellar molecules since they allow important physical properties of the gas and its evolutionary history to be deduced. We update one of the most widely used astrochemical databases to reflect advances in experimental and theoretical estimates of rate coefficients and to respond to the large increase in the number of molecules detected in space since our last release in 2013. We present the sixth release of the UMIST Database for Astrochemistry (UDfA), a major expansion of the gas-phase chemistry that describes the synthesis of interstellar and circumstellar molecules. Since our last release, we have undertaken a major review of the literature which has increased the number of reactions by over 40% to a total of 8767 and increased the number of species by over 55% to 737. We have made a particular attempt to include many of the new species detected in space over the past decade, including those from the QUIJOTE and GOTHAM surveys, as well as providing references to the original data sources. We use the database to investigate the gas-phase chemistries appropriate to O-rich and C-rich conditions in TMC-1 and to the circumstellar envelope of the C-rich AGB star IRC+10216 and identify successes and failures of gas-phase only models. This update is a significant improvement to the UDfA database. For the dark cloud and C-rich circumstellar envelope models, calculations match around 60% of the abundances of observed species to within an order of magnitude. There are a number of detected species, however, that are not included in the model either because their gas-phase chemistry is unknown or because they are likely formed via surface reactions on icy grains. Future laboratory and theoretical work is needed to include such species in reaction networks. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 26 pages, 8 figures, accepted for publication in A&A

Degenerate interpretations of O$_3$ spectral features in exoplanet atmosphere observations due to stellar UV uncertainties: a 3D case study with TRAPPIST-1e

Authors: Gregory Cooke, Dan Marsh, Catherine Walsh, Allison Youngblood

Abstract: TRAPPIST-1e is a potentially habitable terrestrial exoplanet orbiting an ultra-cool M Dwarf star and is a key target for observations with the James Webb Space Telescope (JWST). One-dimensional photochemical modelling of terrestrial planetary atmospheres has shown the importance of the incoming stellar UV flux in modulating the concentration of chemical species, such as O$_3$ and H$_2$O. In additi… ▽ More TRAPPIST-1e is a potentially habitable terrestrial exoplanet orbiting an ultra-cool M Dwarf star and is a key target for observations with the James Webb Space Telescope (JWST). One-dimensional photochemical modelling of terrestrial planetary atmospheres has shown the importance of the incoming stellar UV flux in modulating the concentration of chemical species, such as O$_3$ and H$_2$O. In addition, three-dimensional (3D) modelling has demonstrated anisotropy in chemical abundances due to transport in tidally locked exoplanet simulations. We use the Whole Atmosphere Community Climate Model Version 6 (WACCM6), a 3D Earth System Model, to investigate how uncertainties in the incident UV flux, combined with transport, affect observational predictions for TRAPPIST-1e (assuming an initial Earth-like atmospheric composition). We use two semi-empirical stellar spectra for TRAPPIST-1 from the literature. The UV flux ratio between them can be as large as a factor of 5000 in some wavelength bins. Consequently, the photochemically-produced total O$_3$ columns differ by a factor of 26. Spectral features of O$_3$ in both transmission and emission spectra vary between these simulations (e.g. differences of 19 km in transmission spectra effective altitude for O$_3$ at 0.6 $μ$m). This leads to potential ambiguities when interpreting observations, including overlap with scenarios that assume alternative O$_2$ concentrations. Hence, to achieve robust interpretations of terrestrial exoplanetary spectra, characterisation of the UV spectra of their host stars is critical. In the absence of such stellar measurements, atmospheric context can still be gained from other spectral features (e.g. H$_2$O), or by comparing direct imaging and transmission spectra in conjunction. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 20 pages, 11 figures, accepted in ApJ

Higher Water Loss on Earth-like Exoplanets in Eccentric Orbits

Authors: Binghan Liu, Dan Marsh, Catherine Walsh, Gregory Cooke

Abstract: The climate of a terrestrial exoplanet is controlled by the type of host star, the orbital configuration and the characteristics of the atmosphere and the surface. Many rocky exoplanets have higher eccentricities than those in the Solar System, and about 18% of planets with masses $< 10 \mathrm{M}_{\oplus}$ have $e>0.1$. Underexplored are the implications of such high eccentricities on the atmosph… ▽ More The climate of a terrestrial exoplanet is controlled by the type of host star, the orbital configuration and the characteristics of the atmosphere and the surface. Many rocky exoplanets have higher eccentricities than those in the Solar System, and about 18% of planets with masses $< 10 \mathrm{M}_{\oplus}$ have $e>0.1$. Underexplored are the implications of such high eccentricities on the atmosphere, climate, and potential habitability on such planets. We use WACCM6, a state-of-the-art fully-coupled Earth-system model, to simulate the climates of two Earth-like planets; one in a circular orbit ($e=0$), and one in an eccentric orbit ($e=0.4$) with the same mean insolation. We quantify the effects of eccentricity on the atmospheric water abundance and loss given the importance of liquid water for habitability. The asymmetric temperature response in the eccentric orbit results in a water vapour mixing ratio in the stratosphere ($> 20$ ppmv) that is approximately five times greater than that for circular orbit ($\sim 4$ ppmv). This leads to at most $\sim 3$ times increases in both the atmospheric hydrogen loss rate and the ocean loss rate compared with the circular case. Using the Planetary Spectrum Generator, we simulate the idealised transmission spectra for both cases. We find that the water absorption features are stronger at all wavelengths for the $e=0.4$ spectrum than for the circular case. Hence, highly-eccentric Earth-like exoplanets may be prime targets for future transmission spectroscopy observations to confirm, or otherwise, the presence of atmospheric water vapour. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 23 pages, 6 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, stad1828, Published: 16 June 2023

Disentangling physics and chemistry in AGB outflows: revealing degeneracies when adding complexity

Authors: Marie Van de Sande, Catherine Walsh, Tom J. Millar

Abstract: Observations of the outflows of asymptotic giant branch (AGB) stars continue to reveal their chemical and dynamical complexity. Spherical asymmetries, such as spirals and disks, are prevalent and thought to be caused by binary interaction with a (sub)stellar companion. Furthermore, high density outflows show evidence of dust-gas interactions. The classical chemical model of these outflows - a gas-… ▽ More Observations of the outflows of asymptotic giant branch (AGB) stars continue to reveal their chemical and dynamical complexity. Spherical asymmetries, such as spirals and disks, are prevalent and thought to be caused by binary interaction with a (sub)stellar companion. Furthermore, high density outflows show evidence of dust-gas interactions. The classical chemical model of these outflows - a gas-phase only, spherically symmetric chemical kinetics model - is hence not appropriate for a majority of observed outflows. We have included several physical and chemical advancements step-by-step: a porous density distribution, dust-gas chemistry, and internal UV photons originating from a close-by stellar companion. Now, we combine these layers of complexity into the most chemically and physically advanced chemical kinetics model of AGB outflows to date. By varying over all model parameters, we obtain a holistic view of the outflow's composition and how it (inter)depends on the different complexities. A stellar companion has the largest influence, especially when combined with a porous outflow. We compile sets of gas-phase molecules that trace the importance of dust-gas chemistry and allow us to infer the presence of a companion and porosity of the outflow. This shows that our new chemical model can be used to infer physical and chemical properties of specific outflows, as long as a suitable range of molecules is observed. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Faraday Discussions 2023, accepted manuscript. 15 pages, 7 figures, 4 tables. Supplementary information can be found at https://pubs.rsc.org/en/content/articlelanding/2023/fd/d3fd00039g

Molecules with ALMA at Planet-forming Scales (MAPS). Complex Kinematics in the AS 209 Disk Induced by a Forming Planet and Disk Winds

Authors: Maria Galloway-Sprietsma, Jaehan Bae, Richard Teague, Myriam Benisty, Stefano Facchini, Yuri Aikawa, Felipe Alarcón, Sean M. Andrews, Edwin Bergin, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzmán, Jane Huang, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Karin I. Öberg, Catherine Walsh, David J. Wilner

Abstract: We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a cohere… ▽ More We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a coherent upward flow arising from the gap. The upward gas flow is as fast as $150~{\rm m~s}^{-1}$ in the regions traced by $^{12}$CO emission, which corresponds to about 50% of the local sound speed or $6\%$ of the local Keplerian speed. Such an upward gas flow is difficult to reconcile with an embedded planet alone. Instead, we propose that magnetically driven winds via ambipolar diffusion are triggered by the low gas density within the planet-carved gap, dominating the kinematics of the gap region. We estimate the ambipolar Elsasser number, Am, using the HCO$^+$ column density as a proxy for ion density and find that Am is ~0.1 at the radial location of the upward flow. This value is broadly consistent with the value at which numerical simulations find ambipolar diffusion drives strong winds. We hypothesize the activation of magnetically-driven winds in a planet-carved gap can control the growth of the embedded planet. We provide a scaling relationship which describes the wind-regulated terminal mass: adopting parameters relevant to 100 au from a solar-mass star, we find the wind-regulated terminal mass is about one Jupiter mass, which may help explain the dearth of directly imaged super-Jovian-mass planets. △ Less

Submitted 12 May, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

Comments: This paper has been accepted for publication in the Astrophysical Journal (ApJ)

Azimuthal C/O Variations in a Planet-Forming Disk

Authors: Luke Keyte, Mihkel Kama, Alice S. Booth, Edwin A. Bergin, L. Ilsedore Cleeves, Ewine F. van Dishoeck, Maria N. Drozdovskaya, Kenji Furuya, Jonathan Rawlings, Oliver Shorttle, Catherine Walsh

Abstract: The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and… ▽ More The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and between disks. This is mostly tied to the icelines of major volatile carriers such as CO and H2O. Using ALMA observations of CS and SO, we have unearthed evidence for an entirely novel type of C/O variation in the protoplanetary disk around HD 100546: an azimuthal variation from a typical, oxygen-dominated ratio (C/O=0.5) to a carbon-dominated ratio (C/O>1.0). We show that the spatial distribution and peculiar line kinematics of both CS and SO molecules can be well-explained by azimuthal variations in the C/O ratio. We propose a shadowing mechanism that could lead to such a chemical dichotomy. Our results imply that tracing the formation history of giant exoplanets using their atmospheric C/O ratios will need to take into account time-dependent azimuthal C/O variations in a planet's accretion zone. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: Accepted in Nature Astronomy

UV-driven Chemistry as a Signpost for Late-stage Planet Formation

Authors: Jenny K. Calahan, Edwin A. Bergin, Arthur D. Bosman, Evan Rich, Sean M. Andrews, Jennifer B. Bergner, L. Ilsedore Cleeves, Viviana V. Guzman, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Karin I. Oberg, Richard Teague, Catherine Walsh, David J. Wilner, Ke Zhang

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

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted to Nature Astronomy, Published Dec 8th 2022

Sulphur monoxide emission tracing an embedded planet in the HD 100546 protoplanetary disk

Authors: Alice S. Booth, John D. Ilee, Catherine Walsh, Mihkel Kama, Luke Keyte, Ewine F. van Dishoeck, Hideko Nomura

Abstract: Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming… ▽ More Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming planets. We report high angular resolution ALMA Band 7 observations of sulphur monoxide (SO) in the nearby planet-hosting disk around Herbig star HD 100546. SO is rarely detected in evolved protoplanetary disks but in other environments, it is most often utilised as a tracer of shocks. The SO emission from the HD 100546 disk is primarily originating from gas within the approx. 20 au mm-dust cavity and shows a clear azimuthal brightness asymmetry of a factor of 2. In addition, we see a significant difference in the line profile shape when comparing these new Cycle 7 data to Cycle 0 data of the same SO transitions. We discuss the different physical/chemical mechanisms that could be responsible for this asymmetry and time variability including disk winds, disk warps, and a shock triggered by a (forming) planet. We propose that the SO is enhanced in the cavity due to the presence of a giant planet. The SO asymmetry complements evidence for hot circumplanetary material around the giant planet HD 100546 c traced via CO ro-vibrational emission. This work sets the stage for further observational and modelling efforts to detect and understand the chemical imprint of a forming planet on its parent disk. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: Accepted A&A 24th October

Journal ref: A&A 669, A53 (2023)

Variability due to climate and chemistry in observations of oxygenated Earth-analogue exoplanets

Authors: Gregory Cooke, Dan Marsh, Catherine Walsh, Sarah Rugheimer, Geronimo Villanueva

Abstract: The Great Oxidation Event was a period during which Earth's atmospheric oxygen (O$_2$) concentrations increased from $\sim 10^{-5}$ times its present atmospheric level (PAL) to near modern levels, marking the start of the Proterozoic geological eon 2.4 billion years ago. Using WACCM6, an Earth System Model, we simulate the atmosphere of Earth-analogue exoplanets with O$_2$ mixing ratios between 0.… ▽ More The Great Oxidation Event was a period during which Earth's atmospheric oxygen (O$_2$) concentrations increased from $\sim 10^{-5}$ times its present atmospheric level (PAL) to near modern levels, marking the start of the Proterozoic geological eon 2.4 billion years ago. Using WACCM6, an Earth System Model, we simulate the atmosphere of Earth-analogue exoplanets with O$_2$ mixing ratios between 0.1% and 150% PAL. Using these simulations, we calculate the reflection/emission spectra over multiple orbits using the Planetary Spectrum Generator. We highlight how observer angle, albedo, chemistry, and clouds affect the simulated observations. We show that inter-annual climate variations, as well short-term variations due to clouds, can be observed in our simulated atmospheres with a telescope concept such as LUVOIR or HabEx. Annual variability and seasonal variability can change the planet's reflected flux (including the reflected flux of key spectral features such as O$_2$ and H$_2$O) by up to factors of 5 and 20, respectively, for the same orbital phase. This variability is best observed with a high-throughput coronagraph. For example, HabEx (4 m) with a starshade performs up to a factor of two times better than a LUVOIR B (6 m) style telescope. The variability and signal-to-noise ratio of some spectral features depends non-linearly on atmospheric O$_2$ concentration. This is caused by temperature and chemical column depth variations, as well as generally increased liquid and ice cloud content for atmospheres with O$_2$ concentrations of $<$1% PAL. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS

The gas-phase reaction of NH2 with formaldehyde (CH2O) is not a source of formamide (NH2CHO) in interstellar environments

Authors: Kevin M. Douglas, Daniel Lucas, Catherine Walsh, Niclas A. West, Mark A. Blitz, Dwayne E. Heard

Abstract: The first experimental study of the low-temperature kinetics of the gas-phase reaction of NH2 with formaldehyde (CH2O) has been performed. This reaction has previously been suggested as a source of formamide (NH2CHO) in interstellar environments. A pulsed Laval nozzle equipped with laser-flash photolysis and laser-induced fluorescence spectroscopy was used to create and monitor the temporal decay… ▽ More The first experimental study of the low-temperature kinetics of the gas-phase reaction of NH2 with formaldehyde (CH2O) has been performed. This reaction has previously been suggested as a source of formamide (NH2CHO) in interstellar environments. A pulsed Laval nozzle equipped with laser-flash photolysis and laser-induced fluorescence spectroscopy was used to create and monitor the temporal decay of NH2 in the presence of CH2O. No loss of NH2 could be observed via reaction with CH2O and we place an upper-limit on the rate coefficient of <6x10-12 cm3 molecule-1 s-1 at 34K. Ab initio calculations of the potential energy surface were combined with RRKM calculations to predict a rate coefficient of 6.2x10-14 cm3 molecule-1 s-1 at 35K, consistent with the experimental results. The presence of a significant barrier, 18 kJ mol-1, for the formation of formamide as a product, means that only the H-abstraction channel producing NH3 + CHO, in which the transfer of an H-atom can occur by quantum mechanical tunnelling through a 23 kJ mol-1 barrier, is open at low temperatures. These results are in contrast with a recent theoretical study which suggested that the reaction could proceed without a barrier and was therefore a viable route to gas-phase formamide formation. The calculated rate coefficients were used in an astrochemical model which demonstrated that this reaction produces only negligible amounts of gas-phase formamide under interstellar and circumstellar conditions. The reaction of NH2 with CH2O is therefore not an important source of formamide at low temperatures in interstellar environments. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: Manuscript, 14 pages, 4 figures. Supporting Information, 8 pages, 2 figures. Accepted for publication in The Astrophysical Journal Letters

The Molecular Composition of Shadowed Protosolar Disk Midplanes beyond the Water Snowline

Authors: Shota Notsu, Kazumasa Ohno, Takahiro Ueda, Catherine Walsh, Christian Eistrup, Hideko Nomura

Abstract: The disk midplane temperature is potentially affected by the dust traps/rings. The dust depletion beyond the water snowline will cast a shadow. In this study, we adopt a detailed gas-grain chemical reaction network, and investigate the radial gas and ice abundance distributions of dominant carbon-, oxygen-, and nitrogen-bearing molecules in disks with shadow structures beyond the water snowline ar… ▽ More The disk midplane temperature is potentially affected by the dust traps/rings. The dust depletion beyond the water snowline will cast a shadow. In this study, we adopt a detailed gas-grain chemical reaction network, and investigate the radial gas and ice abundance distributions of dominant carbon-, oxygen-, and nitrogen-bearing molecules in disks with shadow structures beyond the water snowline around a protosolar-like star. In shadowed disks, the dust grains at around $3-8$ au are predicted to have more than around $5-10$ times amounts of ices of organic molecules such as H$_{2}$CO, CH$_{3}$OH, and NH$_{2}$CHO, saturated hydrocarbon ices such as CH$_{4}$ and C$_{2}$H$_{6}$, in addition to H$_{2}$O, CO, CO$_{2}$, NH$_{3}$, N$_{2}$, and HCN ices, compared with those in non-shadowed disks. In the shadowed regions, we find that hydrogenation (especially of CO ice) is the dominant formation mechanism of complex organic molecules. The gas-phase N/O ratios show much larger spatial variations than the gas-phase C/O ratios, thus the N/O ratio is predicted to be a useful tracer of the shadowed region. N$_{2}$H$^{+}$ line emission is a potential tracer of the shadowed region. We conclude that a shadowed region allows the recondensation of key volatiles onto dust grains, provides a region of chemical enrichment of ices that is much closer to the star than within a non-shadowed disk, and may explain to some degree the trapping of O$_{2}$ ice in dust grains that formed comet 67P/Churyumov-Gerasimenko. We discuss that, if formed in a shadowed disk, Jupiter does not need to have migrated vast distances. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 52 pages, 22 Figures, 2 Tables, Accepted for publication in The Astrophysical Journal (ApJ) on August 7th, 2022

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/

Modeling snowline locations in protostars: The impact of the structure of protostellar cloud cores

Authors: Nadia M. Murillo, Tien-Hao Hsieh, Catherine Walsh

Abstract: Abridged Context: Snowlines during star and disk formation are responsible for a range of effects during the evolution of protostars, such as setting the chemical composition of the envelope and disk. This in turn influences the formation of planets by changing the elemental compositions of solids and affecting the collisional properties and outcomes of dust grains. Snowlines can also reveal acc… ▽ More Abridged Context: Snowlines during star and disk formation are responsible for a range of effects during the evolution of protostars, such as setting the chemical composition of the envelope and disk. This in turn influences the formation of planets by changing the elemental compositions of solids and affecting the collisional properties and outcomes of dust grains. Snowlines can also reveal accretion bursts, providing insight into the formation process of stars. Methods: A numerical chemical network coupled with a grid of cylindrical-symmetric physical models was used to identify what physical parameters alter the CO and H$_2$O snowline locations. The investigated parameters are the initial molecular abundances, binding energies of CO and H$_2$O, heating source, cloud core density, outflow cavity opening angle, and disk geometry. Simulated molecular line emission maps were used to quantify the change in the snowline location with each parameter. Conclusions: The models presented in this work show that the CO and H$_2$O snowline locations do not occur at a single, well-defined temperature as is commonly assumed. Instead, the snowline position depends on luminosity, cloud core density, and whether a disk is present or not. Inclination and spatial resolution affect the observability and successful measurement of snowline locations. We note that N$_2$H$^+$ and HCO$^+$ emission serve as good observational tracers of CO and H$_2$O snowline locations. However, constraints on whether or not a disk is present, the observation of additional molecular tracers, and estimating envelope density will help in accurately determining the cause of the observed snowline position. Plots of the N$_2$H$^+$ and HCO$^+$ peak emission radius versus luminosity are provided to compare the models with observations of deeply embedded protostars aiming to measure the CO and H$_2$O snowline locations. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 29 pages, 21 figures, accepted for publication in A&A

Journal ref: A&A 665, A68 (2022)

Testing White Dwarf Age Estimates using Wide Double White Dwarf Binaries from Gaia EDR3

Authors: Tyler M. Heintz, J. J. Hermes, Kareem El-Badry, Charlie Walsh, Jennifer L. van Saders, C. E. Fields, Detlev Koester

Abstract: White dwarf (WD) stars evolve simply and predictably, making them reliable age indicators. However, self-consistent validation of the methods for determining WD total ages has yet to be widely performed. This work uses 1565 wide ( > 100 au) WD+WD binaries and 24 new triples containing at least two WDs to test the accuracy and validity of WD total age determinations. For these 1589 wide double-WD b… ▽ More White dwarf (WD) stars evolve simply and predictably, making them reliable age indicators. However, self-consistent validation of the methods for determining WD total ages has yet to be widely performed. This work uses 1565 wide ( > 100 au) WD+WD binaries and 24 new triples containing at least two WDs to test the accuracy and validity of WD total age determinations. For these 1589 wide double-WD binaries and triples, we derive total ages of each WD using photometric data from all-sky surveys, in conjunction with Gaia parallaxes and current hydrogen-atmosphere WD models. Ignoring initial-to-final-mass relations and considering only WD cooling ages, we find that roughly 21-36% of the more massive WDs in a system have a shorter cooling age. Since more massive WDs should be born as more massive main-sequence stars, we attribute this unphysical disagreement as evidence of prior mergers or the presence of an unresolved companion, suggesting that roughly 21-36% of wide WD+WD binaries were once triples. Among the 423 wide WD+WD pairs that pass high-fidelity cuts, we find that 25% total age uncertainties are generally appropriate for WDs with masses > 0.63 Msun and temperatures < 12,000 K, and provide suggested inflation factors for age uncertainties for higher-mass WDs. Overall, WDs return reliable stellar ages, but we detail cases where total ages are least reliable, especially for WDs < 0.63 Msun. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Comments: 24 pages, 10 figures. Submitted to ApJ. Catalog available at https://sites.bu.edu/buwd/files/2022/05/Table_A1.csv

Unveiling the outer dust disc of TW Hya with deep ALMA observations

Authors: John D. Ilee, Catherine Walsh, Jeff Jennings, Richard A. Booth, Giovanni P. Rosotti, Richard Teague, Takashi Tsukagoshi, Hideko Nomura

Abstract: The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALM… ▽ More The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALMA observations at 1 mm of the nearest protoplanetary disc, TW Hya. Using the visibility analysis tool `frank', we reveal a structured millimetre intensity distribution out to 100 au, well beyond previous estimates of 60-70 au. Our analysis suggests the presence of a new millimetre continuum gap at 82 au, coincident with similar features seen in optical/near-infrared scattered light and millimetre molecular line observations. Examination of the fit residuals confirms the presence of the previously reported au-scale continuum excess at 52 au (P.A. = 242.5 degrees). Our results demonstrate the utility of combining deep, moderate resolution observations with super-resolution analysis techniques to probe the faintest regions of protoplanetary discs. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

Fevering Interstellar Ices Have More CH3OD

Authors: Beatrice M. Kulterer, Maria N. Drozdovskaya, Stefano Antonellini, Catherine Walsh, Tom J. Millar

Abstract: Mono-deuterated methanol is thought to form during the prestellar core stage of star formation. Observed variations in the CH2DOH/CH3OD ratio suggest that its formation is strongly dependent on the surrounding cloud conditions. Thus, it is a potential tracer of the physical conditions before the onset of star formation. A single-point physical model representative of a typical prestellar core is c… ▽ More Mono-deuterated methanol is thought to form during the prestellar core stage of star formation. Observed variations in the CH2DOH/CH3OD ratio suggest that its formation is strongly dependent on the surrounding cloud conditions. Thus, it is a potential tracer of the physical conditions before the onset of star formation. A single-point physical model representative of a typical prestellar core is coupled to chemical models to investigate potential formation pathways towards deuterated methanol at the prestellar stage. Simple addition reactions of H and D are not able to reproduce observed abundances. The implementation of an experimentally verified abstraction scheme leads to the efficient formation of methyl-deuterated methanol, but lacks sufficient formation of hydroxy-deuterated methanol. CH3OD is most likely formed at a later evolutionarymstage, potentially from H-D exchange reactions in warm ices between HDO (and D2O) and CH3OH. The CH2DOH/CH3OD ratio is not an appropriate tracer of the physical conditions during the prestellar stage, but might be better suited as a tracer of ice heating. △ Less

Submitted 3 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ACS Earth and Space Chemistry; 56 pages, 7 figures

Surface astrochemistry: a computational chemistry perspective

Authors: H. M. Cuppen, A. Fredon, T. Lamberts, E. M. Penteado, M. Simons, C. Walsh

Abstract: Molecules in space are synthesized via a large variety of gas-phase reactions, and reactions on dust-grain surfaces, where the surface acts as a catalyst. Especially, saturated, hydrogen-rich molecules are formed through surface chemistry. Astrochemical models have developed over the decades to understand the molecular processes in the interstellar medium, taking into account grain surface chemist… ▽ More Molecules in space are synthesized via a large variety of gas-phase reactions, and reactions on dust-grain surfaces, where the surface acts as a catalyst. Especially, saturated, hydrogen-rich molecules are formed through surface chemistry. Astrochemical models have developed over the decades to understand the molecular processes in the interstellar medium, taking into account grain surface chemistry. However, essential input information for gas-grain models, such as binding energies of molecules to the surface, have been derived experimentally only for a handful of species, leaving hundreds of species with highly uncertain estimates. Moreover, some fundamental processes are not well enough constrained to implement these into the models. The proceedings gives three examples how computational chemistry techniques can help answer fundamental questions regarding grain surface chemistry. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: (12 pages, 8 figures)

Journal ref: Astrochemistry VII -- Through the Cosmos from Galaxies to Planets", proceedings of the IAU Symposium No. 332, 2017, Puerto Varas, Chile; eds. M. Cunningham, T. Millar and Y. Aikawa,p 293 - 304

Inside-Out Planet Formation. VII. Astrochemical Models of Protoplanetary Disks and Implications for Planetary Compositions

Authors: Arturo Cevallos Soto, Jonathan C. Tan, Xiao Hu, Chia-Jung Hsu, Catherine Walsh

Abstract: Inside-Out Planet Formation (IOPF) proposes that the abundant systems of close-in Super-Earths and Mini-Neptunes form in situ at the pressure maximum associated with the Dead Zone Inner Boundary (DZIB). We present a model of physical and chemical evolution of protoplanetary disk midplanes that follows gas advection, radial drift of pebbles and gas-grain chemistry to predict abundances from 300~au… ▽ More Inside-Out Planet Formation (IOPF) proposes that the abundant systems of close-in Super-Earths and Mini-Neptunes form in situ at the pressure maximum associated with the Dead Zone Inner Boundary (DZIB). We present a model of physical and chemical evolution of protoplanetary disk midplanes that follows gas advection, radial drift of pebbles and gas-grain chemistry to predict abundances from 300~au down to the DZIB near 0.2 au. We consider typical disk properties relevant for IOPF, i.e., accretion rates 1E-9 < dM/dt / (Msun/yr) < 1E-8 and viscosity parameter alpha = 1E-4, and evolve for fiducial duration of t = 1E5 years. For outer, cool disk regions, we find that C and up to 90% of O nuclei start locked in CO and O2 ice, which keeps abundances of CO2 and H2O one order of magnitude lower. Radial drift of icy pebbles is influential, with gas-phase abundances of volatiles enhanced up to two orders of magnitude at ice-lines, while the outer disk becomes depleted of dust. Disks with decreasing accretion rates gradually cool, which draws in icelines closer to the star. At <~1 au, advective models yield water-rich gas with C/O ratios <~ 0.1, which may be inherited by atmospheres of planets forming here via IOPF. For planetary interiors built by pebble accretion, IOPF predicts volatile-poor compositions. However, advectively-enhanced volatile mass fractions of ~10% can occur at the water ice line. △ Less

Submitted 13 September, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: Submitted to MNRAS; comments welcome; for summary of Inside-Out Planet Formation papers see http://www.cosmicorigins.space/iopf

ALMA High-resolution Multiband Analysis for the Protoplanetary Disk around TW Hya

Authors: Takashi Tsukagoshi, Hideko Nomura, Takayuki Muto, Ryohei Kawabe, Kazuhiro D. Kanagawa, Satoshi Okuzumi, Shigeru Ida, Catherine Walsh, Tom J. Millar, Sanemichi Z. Takahashi, Jun Hashimoto, Taichi Uyama, Motohide Tamura

Abstract: We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstr… ▽ More We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstruct the images; multi-frequency synthesis (MFS) and the fiducial image-oriented method, where each band is imaged separately and the frequency dependence is fitted pixel by pixel. We find that the MFS imaging with the second order of Taylor expansion can reproduce the frequency dependence of the continuum emission between Bands 3 and 7 in a manner consistent with previous studies and is a reasonable method to reconstruct the spectral index map. The image-oriented method provides a spectral index map consistent with the MFS imaging, but with a two times lower resolution. Mock observations of an intensity model were conducted to validate the images from the two methods. We find that the MFS imaging provides a high-resolution spectral index distribution with an uncertainty of $<10$~\%. Using the submillimeter spectrum reproduced from our MFS images, we directly calculated the optical depth, power-law index of the dust opacity coefficient ($β$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $β$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Comments: 17pages, 12 figures, Accepted for publication in The Astrophysical Journal

Astrochemical modelling of infrared dark clouds

Authors: Negar Entekhabi, Jonathan C. Tan, Giuliana Cosentino, Chia-Jung Hsu, Paola Caselli, Catherine Walsh, Wanggi Lim, Jonathan D. Henshaw, Ashley T. Barnes, Francesco Fontani, Izaskun Jiménez-Serra

Abstract: Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. We aim to determine the astrochemical conditions, especially cosmic ray ionisatio… ▽ More Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. We aim to determine the astrochemical conditions, especially cosmic ray ionisation rate (CRIR) and chemical age, in different regions of the massive IRDC G28.37+00.07 by comparing observed abundances of multiple molecules and molecular ions with the predictions of astrochemical models. We have computed a series of single-zone astrochemical models with a gas-grain network that systematically explores the parameter space of density, temperature, CRIR, and visual extinction. We have also investigated the effects of choices of CO ice binding energy and temperatures achieved in transient heating of grains when struck by cosmic rays. We selected 10 positions across the IRDC that are known to have a variety of star formation activity. We utilised mid-infrared (MIR) extinction maps and sub-mm emission maps to measure the mass surface densities of these regions, needed for abundance and volume density estimates. The sub-mm emission maps were also used to measure temperatures. We then used IRAM-30m observations of various tracers to estimate column densities and thus abundances. Using estimates of the abundances of CO, HCO$^+$ and N$_2$H$^+$ we find consistency with astrochemical models that have relatively low CRIRs of $ζ\sim10^{-18}$ to $\sim10^{-17}\:{\rm s}^{-1}$, with no evidence for systematic variation with the level of star formation activity. Astrochemical ages are found to be < 1 Myr. We discuss potential sources of systematic uncertainties in these results and the overall implications for IRDC evolutionary history and astrochemical models.(abridged for arXiv) △ Less

Submitted 23 March, 2022; v1 submitted 7 November, 2021; originally announced November 2021.

Comments: Accepted by A&A

Journal ref: A&A 662, A39 (2022)

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

Different molecular filament widths as tracers of accretion onto filaments

Authors: Gilberto C. Gómez, Catherine Walsh, Aina Palau

Abstract: We explore how dense filament widths, when measured using different molecular species, may change as a consequence of gas accretion toward the filament. As a gas parcel falls into the filament, it will experience different density, temperature, and extinction values. The rate at which this environment changes will affect differently the abundance of different molecules. So, a molecule that forms q… ▽ More We explore how dense filament widths, when measured using different molecular species, may change as a consequence of gas accretion toward the filament. As a gas parcel falls into the filament, it will experience different density, temperature, and extinction values. The rate at which this environment changes will affect differently the abundance of different molecules. So, a molecule that forms quickly will better reflect the local physical conditions a gas parcel experiences than a slower-forming molecule. Since these differences depend on how the respective timescales compare, the different molecular distributions should reflect how rapidly the environment changes, i.e., the accretion rate toward the filament. We find that the filament widths measured from time-dependent abundances for C2H, CO, CN, CS, and C3H2, are sensitive the most to this effect, being those molecules the ones presenting also the wider filament widths. On the contrary, molecules such as N2H+, NH3, H2CO, HNC or CH3OH are not so sensitive to accretion and present the narrowest filament widths. We propose that ratios of filament widths for different tracers could be a useful tool to estimate the accretion rate onto the filament. △ Less

Submitted 30 March, 2022; v1 submitted 12 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS