Showing 1–50 of 122 results for author: Booth, A

Supernova Pointing Capabilities of DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electr… ▽ More The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on $^{40}$Ar and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called ``brems flipping'', as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE's burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage. △ Less

Submitted 14 July, 2024; originally announced July 2024.

Comments: 25 pages, 16 figures

Report number: FERMILAB-PUB-24-0319-LBNF

Viscous circumbinary protoplanetary discs -- I. Structure of the inner cavity

Authors: Anna B. T. Penzlin, Richard A. Booth, Richard P. Nelson, Christoph M. Schäfer, Wilhelm Kley

Abstract: Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficul… ▽ More Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficult problem, for which there is no complete analytical solution, that predicts the shape of the observed disc accurately. We use numeric simulation to evolve circumbinary discs with varying disc viscosities and investigate the size and shape of the inner cavities in such protoplanetary discs. We have simulated over 140 locally isothermal 2D grid-based disc models for > 3e4 binary orbits each and mapped out the parameter space relevant for protoplanetary discs. With this, it becomes possible to create parametrised profiles for individual discs to compare to observation and find limits to their binary eccentricity or internal viscosity from the simulation data. In the long-term simulations larger cavity sizes than previously considered are possible within the parameter space (< 6 binary separations). As an example, we find that the eccentricity of the disc around HD 142527 suggests the impact of the binary dynamics on the disc. However, even considering the larger cavity sizes, the large size of the cavity in HD 142527 remains unexplained by the simulations considering the most recent orbital constraints. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 13 pages, 14 figures (appendix: 2 pages, 3 figures); accepted by mnras

BOWIE-ALIGN: How formation and migration histories of giant planets impact atmospheric compositions

Authors: Anna B. T. Penzlin, Richard A. Booth, James Kirk, James E. Owen, Eva-Maria Ahrer, Duncan A. Christie, Alastair B. Claringbold, Emma Esparza-Borges, M. López-Morales, N. J. Mayne, Mason McCormack, Annabella Meech, Vatsal Panwar, Diana Powell, Denis E. Sergeev, Jake Taylor, Peter J. Wheatley, Maria Zamyatina

Abstract: Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the composition… ▽ More Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the compositions of 8 hot Jupiters around F stars, 4 with orbits aligned with the stellar rotation axis and 4 misaligned. Using the alignment as an indicator for planets that underwent disc migration or high-eccentricity migration, one can determine whether migration history produces notable differences in composition between the two samples of planets. This paper describes the planet formation model that motivates our observing programme. Our model traces the accretion of chemical components from the gas and dust in the disc over a broad parameter space to create a full, unbiased model sample from which we can estimate the range of final atmospheric compositions. For high metallicity atmospheres (O/H > 10 times solar), the C/O ratios of aligned and misaligned planets diverge, with aligned planets having lower C/O (< 0.25) due to the accretion of oxygen-rich silicates from the inner disc. However, silicates may rain out instead of releasing their oxygen into the atmosphere. This would significantly increase the C/O of aligned planets (C/O > 0.6), inverting the trend between the aligned and misaligned planets. Nevertheless, by comparing statistically significant samples of aligned and misaligned planets, we expect atmospheric composition to constrain how planets form. △ Less

Submitted 4 July, 2024; v1 submitted 3 July, 2024; originally announced July 2024.

Comments: 11pages 10 figures, (appendix: 6 page, 4 figures), submitted to mnras

BOWIE-ALIGN: A JWST comparative survey of aligned vs misaligned hot Jupiters to test the dependence of atmospheric composition on migration history

Authors: James Kirk, Eva-Maria Ahrer, Anna B. T. Penzlin, James E. Owen, Richard A. Booth, Lili Alderson, Duncan A. Christie, Alastair B. Claringbold, Emma Esparza-Borges, Chloe E. Fisher, Mercedes López-Morales, N. J. Mayne, Mason McCormack, Annabella Meech, Vatsal Panwar, Diana Powell, Jake Taylor, Denis E. Sergeev, Daniel Valentine, Hannah R. Wakeford, Peter J. Wheatley, Maria Zamyatina

Abstract: A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft… ▽ More A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft break with different orbital alignments. Under the assumption that aligned planets migrate through the inner disc, while misaligned planets migrate after disc dispersal, the act of migrating through the inner disc should lead to a measurable difference in the C/O between aligned and misaligned planets. We expect the amplitude and sign of this difference to depend on the amount of planetesimal accretion and whether silicates accreted from the inner disc release their oxygen. Here, we identify all known exoplanets that are suitable for testing this hypothesis, describe our JWST survey, and use noise simulations and atmospheric retrievals to estimate our survey's sensitivity. With the selected sample of four aligned and four misaligned hot Jupiters, we will be sensitive to the predicted differences in C/O between aligned and misaligned hot Jupiters for a wide range of model scenarios. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 13 pages, 8 figures, submitted to RASTI

C/O Ratios and the formation of wide separation exoplanets

Authors: Edwin A. Bergin, Richard A. Booth, Maria Jose Colmenares, John D. Ilee

Abstract: The gas and solid-state C/O ratios provide context to potentially link the atmospheric composition of planets to that of the natal disk. We provide a synthesis of extant estimates of the gaseous C/O and C/H ratios in planet-forming disks obtained primarily through analysis of Atacama Large Millimeter Array observations. These estimates are compared to atmospheric abundances of wide separation (> 1… ▽ More The gas and solid-state C/O ratios provide context to potentially link the atmospheric composition of planets to that of the natal disk. We provide a synthesis of extant estimates of the gaseous C/O and C/H ratios in planet-forming disks obtained primarily through analysis of Atacama Large Millimeter Array observations. These estimates are compared to atmospheric abundances of wide separation (> 10 au) gas giants. The resolved disk gas C/O ratios, from seven systems, generally exhibit C/O > 1 with sub-solar, or depleted, carbon content. In contrast, wide separation gas giants have atmospheric C/O ratios that cluster near or slightly above the presumed stellar value with a range of elemental C/H. From the existing disk composition, we infer that the solid-state mm/cm-sized pebbles have a total C/O ratio (solid cores and ices) that is solar (stellar) in content. We explore simple models that reconstruct the exoplanet atmospheric composition from the disk, while accounting for silicate cloud formation in the planet atmosphere. If wide separate planets formed via the core-accretion mechanism, they must acquire their metals from pebble or planetesimal accretion. Further, the dispersion in giant planet C/H content is best matched by a disk composition with modest and variable factors of carbon depletion. An origin of the wide separation gas giants via gravitational instability cannot be ruled out as stellar C/O ratios should natively form in this scenario. However, the variation in planet metallicity with a stellar C/O ratio potentially presents challenges to these models. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 17 pages, 6 figures, accepted by Astrophysical Journal Letters

Seeing the unseen: a method to detect unresolved rings in protoplanetary disks

Authors: Chiara E. Scardoni, Richard A. Booth, Cathie J. Clarke, Giovanni P. Rosotti, Alvaro Ribas

Abstract: While high resolution ALMA observations reveal a wealth of substructure in protoplanetary discs, they remain incapable of resolving the types of small scale dust structures predicted, for example, by numerical simulations of the streaming instability. In this Letter, we propose a method to find evidence for unresolved, optically thick dusty rings in protoplanetary disks. We demonstrate that, in pr… ▽ More While high resolution ALMA observations reveal a wealth of substructure in protoplanetary discs, they remain incapable of resolving the types of small scale dust structures predicted, for example, by numerical simulations of the streaming instability. In this Letter, we propose a method to find evidence for unresolved, optically thick dusty rings in protoplanetary disks. We demonstrate that, in presence of unresolved rings, the brightness of an inclined disc exhibits a distinctive emission peak at the minor axis. Furthermore, the azimuthal brightness depends on both the geometry of the rings and the dust optical properties; we can therefore use the azimuthal brightness variations to both detect unresolved rings and probe their properties. By analyzing the azimuthal brightness in the test-case of ring-like substructures formed by streaming instability, we show that the resulting peak is likely detectable by ALMA for typical disc parameters. Moreover, we present an analytic model that not only qualitatively but also quantitatively reproduces the peak found in the simulations, validating its applicability to infer the presence of unresolved rings in observations and characterize their optical properties and shape. This will contribute to the identification of disk regions where streaming instability (and thus planet formation) is occurring. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Accepted by ApJ, 10 pages, 6 figures

Faraday tomography with CHIME: the `tadpole' feature G137+7

Authors: Nasser Mohammed, Anna Ordog, Rebecca A. Booth, Andrea Bracco, Jo-Anne C. Brown, Ettore Carretti, John M. Dickey, Simon Foreman, Mark Halpern, Marijke Haverkorn, Alex S. Hill, Joseph W Kania, Roland Kothes, T. L. Landecker, Joshua MacEachern, Kiyoshi W. Masui, Aimee Menard, Ryan R. Ransom, Wolfgang Reich, Patricia Reich, J. Richard Shaw, Seth R. Siegel, Mehrnoosh Tahani, Alec J. M. Thomson, Tristan Pinsonneault-Marotte , et al. (4 additional authors not shown)

Abstract: A direct consequence of Faraday rotation is that the polarized radio sky does not resemble the total intensity sky at long wavelengths. We analyze G137+7, which is undetectable in total intensity but appears as a depolarization feature. We use the first polarization maps from the Canadian Hydrogen Intensity Mapping Experiment. Our $400-729$ MHz bandwidth and angular resolution, $17'$ to $30'$, all… ▽ More A direct consequence of Faraday rotation is that the polarized radio sky does not resemble the total intensity sky at long wavelengths. We analyze G137+7, which is undetectable in total intensity but appears as a depolarization feature. We use the first polarization maps from the Canadian Hydrogen Intensity Mapping Experiment. Our $400-729$ MHz bandwidth and angular resolution, $17'$ to $30'$, allow us to use Faraday synthesis to analyze the polarization structure. In polarized intensity and polarization angle maps, we find a "tail" extending $10^\circ$ from the "head" and designate the combined object the "tadpole". Similar polarization angles, distinct from the background, indicate that the head and tail are physically associated. The head appears as a depolarized ring in single channels, but wideband observations show that it is a Faraday rotation feature. Our investigations of H I and H$α$ find no connections to the tadpole. The tail suggests motion of either the gas or an ionizing star through the ISM; the B2(e) star HD 20336 is a candidate. While the head features a coherent, $\sim -8$ rad m$^2$ Faraday depth, Faraday synthesis also identifies multiple components in both the head and tail. We verify the locations of the components in the spectra using QU fitting. Our results show that $\sim$octave-bandwidth Faraday rotation observations at $\sim 600$ MHz are sensitive to low-density ionized or partially-ionized gas which is undetectable in other tracers. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: Resubmitted to ApJ after referee report

Chemistry across dust and gas gaps in protoplanetary disks: modelling the co-spatial molecular rings in the HD 100546 disk

Authors: M. Leemker, A. S. Booth, E. F. van Dishoeck, L. Wölfer, B. Dent

Abstract: High-resolution observations show that typically both the dust and the gas in nearby extended protoplanetary disks are structured, possibly related to radial and azimuthal variations in the disk density and/or chemistry. The aim of this work is to identify the expected location and intensity of rings seen in molecular line emission of HCN, CN, C$_2$H, NO, [CI], and HCO$^+$ in gapped disks while ex… ▽ More High-resolution observations show that typically both the dust and the gas in nearby extended protoplanetary disks are structured, possibly related to radial and azimuthal variations in the disk density and/or chemistry. The aim of this work is to identify the expected location and intensity of rings seen in molecular line emission of HCN, CN, C$_2$H, NO, [CI], and HCO$^+$ in gapped disks while exploring a range of physical conditions across the gap. In particular, we model HD 100546 disk where molecular rings are co-spatial with the dust rings at $\sim$20 and $\sim$200 au, in contrast to most other gapped disks. The fiducial model of a gapped disk with a 15 au gas cavity, a 20 au dust cavity, and a shallow (a factor of $\lesssim10$) gas and deep dust gap at 40-175 au provides a good fit to the continuum, CO isotopologues, HCN, and HCO$^+$ in the HD 100546 disk. However, the predictions for [CI], CN, C$_2$H and NO do not match the intensity nor the morphology of the observations. An exploration of the parameter space shows that in general the molecular emission rings are only co-spatial with the dust rings if the gas gap between the dust rings is depleted by at least four orders of magnitude in gas or if the C/O ratio of the gas is varying as a function of radius. For shallower gaps the decrease in the UV field roughly balances the effect of a higher gas density for UV tracers such as CN, C$_2$H, and NO. Therefore, these radicals are not good tracers of the gas gap depth. The C/O ratio primarily effects the intensity of the lines without changing the morphology much. The co-spatial rings observed in the HD 100546 disk could be indicative of a radially varying C/O ratio in the HD 100546 disk with a C/O above 1 in a narrow region across the dust rings, together with a shallow gas gap that is depleted by a factor of $\sim$10 in gas, and a reduced background UV field. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A

Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs

Authors: Alison K. Young, Maggie Celeste, Richard A. Booth, Ken Rice, Adam Koval, Ethan Carter, Dimitris Stamatellos

Abstract: The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate… ▽ More The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths respectively. The second, the "Modified Lombardi" method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the Modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted to MNRAS. 10 pages

Searching for planet-driven dust spirals in ALMA visibilities

Authors: Edward T. Stevenson, Álvaro Ribas, Jessica Speedie, Richard A. Booth, Cathie J. Clarke

Abstract: ALMA (Atacama Large Millimetre/submillimetre Array) observations of the thermal emission from protoplanetary disc dust have revealed a wealth of substructures that could evidence embedded planets, but planet-driven spirals, one of the more compelling lines of evidence, remain relatively rare. Existing works have focused on detecting these spirals using methods that operate in image space. Here, we… ▽ More ALMA (Atacama Large Millimetre/submillimetre Array) observations of the thermal emission from protoplanetary disc dust have revealed a wealth of substructures that could evidence embedded planets, but planet-driven spirals, one of the more compelling lines of evidence, remain relatively rare. Existing works have focused on detecting these spirals using methods that operate in image space. Here, we explore the planet detection capabilities of fitting planet-driven spirals to disc observations directly in visibility space. We test our method on synthetic ALMA observations of planet-containing model discs for a range of disc/observational parameters, finding it significantly outperforms image residuals in identifying spirals in these observations and is able to identify spirals in regions of the parameter space in which no gaps are detected. These tests suggest that a visibility-space fitting approach warrants further investigation and may be able to find planet-driven spirals in observations that have not yet been found with existing approaches. We also test our method on six discs in the Taurus molecular cloud observed with ALMA at 1.33 mm, but find no evidence for planet-driven spirals. We find that the minimum planet masses necessary to drive detectable spirals range from $\approx$ 0.03 to 0.5 $M_{\text{Jup}}$ over orbital radii of 10 to 100 au, with planet masses below these thresholds potentially hiding in such disc observations. Conversely, we suggest that planets $\gtrsim$ 0.5 to 1 $M_{\text{Jup}}$ can likely be ruled out over orbital radii of $\approx$ 20 to 60 au on the grounds that we would have detected them if they were present. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 16 pages, 14 figures, 3 tables. Accepted 2024 April 18 for publication in MNRAS

Introducing cuDisc: a 2D code for protoplanetary disc structure and evolution calculations

Authors: Alfie Robinson, Richard A. Booth, James E. Owen

Abstract: We present a new 2D axisymmetric code, cuDisc, for studying protoplanetary discs, focusing on the self-consistent calculation of dust dynamics, grain size distribution and disc temperature. Self-consistently studying these physical processes is essential for many disc problems, such as structure formation and dust removal, given that the processes heavily depend on one another. To follow the evolu… ▽ More We present a new 2D axisymmetric code, cuDisc, for studying protoplanetary discs, focusing on the self-consistent calculation of dust dynamics, grain size distribution and disc temperature. Self-consistently studying these physical processes is essential for many disc problems, such as structure formation and dust removal, given that the processes heavily depend on one another. To follow the evolution over substantial fractions of the disc lifetime, cuDisc uses the CUDA language and libraries to speed up the code through GPU acceleration. cuDisc employs a second-order finite-volume Godonuv solver for dust dynamics, solves the Smoluchowski equation for dust growth and calculates radiative transfer using a multi-frequency hybrid ray-tracing/flux-limited-diffusion method. We benchmark our code against current state-of-the-art codes. Through studying steady-state problems, we find that including 2D structure reveals that when collisions are important, the dust vertical structure appears to reach a diffusion-settling-coagulation equilibrium that can differ substantially from standard models that ignore coagulation. For low fragmentation velocities, we find an enhancement of intermediate-sized dust grains at heights of ~ 1 gas scale height due to the variation in collision rates with height, and for large fragmentation velocities, we find an enhancement of small grains around the disc mid-plane due to collisional ''sweeping'' of small grains by large grains. These results could be important for the analysis of disc SEDs or scattered light images, given these observables are sensitive to the vertical grain distribution. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 18 pages, accepted for publication in MNRAS. GitHub repository link: https://github.com/cuDisc/cuDisc

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 Role of Drag and Gravity on Dust Concentration in a Gravitationally Unstable Disc

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru, James Wurster, Hossam Aly, Richard Alexander, Ken Rice, Richard A. Booth

Abstract: We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If w… ▽ More We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If we include both, grains that would otherwise be partially decoupled will become well coupled and trace the spirals. For the dust grains most influenced by drag (with Stokes numbers near unity), the dust disc quickly becomes gravitationally unstable and rapidly forms clumps with masses between 0.15 - 6 Earth masses. A large fraction of clumps are below the threshold where runaway gas accretion can occur. However, if dust self-gravity is neglected, the dust is unable to form clumps, despite still becoming trapped in the gas spirals. When large dust grains are unable to feel either gas gravity or drag, the dust is unable to trace the gas spirals. Hence, full physics is needed to properly simulate dust in gravitationally unstable discs. Dust trapping of large grains in spiral arms of discs stable to gas fragmentation could explain planet formation in very young discs by a population of planetesimals formed due to the combined roles of drag and gravity in the earliest stages of a disc's evolution. Furthermore, it highlights that gravitationally unstable discs are not just important for forming gas giants quickly, it can also rapidly form Earth mass bodies. △ Less

Submitted 24 January, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. Updated references

Constraining the gas mass of Herbig disks using CO isotopologues

Authors: L. M. Stapper, M. R. Hogerheijde, E. F. van Dishoeck, L. Lin, A. Ahmadi, A. S. Booth, S. L. Grant, K. Immer, M. Leemker, A. F. Pérez-Sánchez

Abstract: The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing fo… ▽ More The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing for the possibility of using CO as a more effective total gas mass tracer. Using ALMA archival data and new NOEMA data of 12CO, 13CO, and C18O transitions of 35 Herbig disks within 450 pc, the masses are determined using the thermo-chemical code Dust And LInes (DALI). The majority of Herbig disks for which 13CO and C18O are detected are optically thick in both. Computing the gas mass using a simple optically thin relation between line flux and column density results in an underestimate of the gas mass of at least an order of magnitude compared to the masses obtained with DALI. The inferred gas masses with DALI are consistent with a gas-to-dust ratio of at least 100. These gas-to-dust ratios are two orders of magnitude higher compared to those found for T Tauri disks using similar techniques, even over multiple orders of magnitude in dust mass, illustrating the importance of chemical conversion of CO in colder T Tauri disks. Similar high gas-to-dust ratios are found for Herbig group I and II disks. Since group II disks have dust masses comparable to T Tauri disks, their higher CO gas masses illustrate the determining role of temperature. Compared to debris disks, Herbig disks have four orders of magnitude higher gas masses. At least one Herbig disk, HD 163296, has a detected molecular disk wind, but our investigation has not turned up other detections of the CO disk wind in spite of similar sensitivities. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy and Astrophysics. 25 pages, 19 figures, plus appendices

Observing planetesimal formation under streaming instability in the rings of HD 163296

Authors: Francesco Zagaria, Cathie J. Clarke, Richard A. Booth, Stefano Facchini, Giovanni P. Rosotti

Abstract: We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a… ▽ More We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a factor of two. We combine our results with literature measurements of the dust surface density and grain size to determine the dust-to-gas ratio and Stokes number in the 67 au and 100 au rings. Our estimates suggest that particle clumping is taking place under the effect of streaming instability (SI) in the 100 au ring. Even though in the presence of external isotropic turbulence this process might be hindered, we provide evidence that turbulence is non-isotropic in both rings and likely originating from mechanisms (such as ambipolar diffusion) that could ease particle clumping under SI. Finally, we determine the mass accretion rate under the assumption that the disc is in steady state and turbulence regulates angular momentum transport. Our results are in tension with spectroscopic measurements and suggest that other mechanisms might be responsible for accretion, in qualitative agreement with the detection of a magneto-centrifugal wind in this system. Applying our method to larger samples can be used to statistically assess if SI is a viable mechanism to form planetesimals in bright rings. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 13 pages, 4 figures; accepted for publication on ApJL

Formation of misaligned second-generation discs through flyby encounters

Authors: Jeremy L. Smallwood, Rebecca Nealon, Nicolás Cuello, Ruobing Dong, Richard A. Booth

Abstract: Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary di… ▽ More Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary disc, respectively. In both our $N$--body and hydrodynamic simulations, we find that the captured, second-generation disc forms at a tilt twice the initial flyby tilt. This relationship is robust to variations in the flyby's tilt, position angle, periastron, and mass. We extend this concept by also simulating the case where the flyby has a disc of material prior to the encounter but we do not find the same trend. An inclined disc with respect to the primary disc around a misaligned flyby is tilted by a few degrees, remaining close to its initial disc tilt. Therefore, if a disc is present around the flyby before the encounter, the disc may not tilt up to twice the perturber tilt depending on the balance between the angular momentum of the circumsecondary disc and captured particles. In the case where the perturber has no initial disc, analyzing the orientation of these second-generation discs can give information about the orbital properties of the flyby encounter. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 17 pages, 15 figures, accepted to MNRAS

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

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

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

Submitted 22 August, 2023; originally announced August 2023.

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

Tracing snowlines and C/O ratio in a planet-hosting disk: ALMA molecular line observations towards the HD169142 disk

Authors: Alice S. Booth, Charles J. Law, Milou Temmink, Margot Leemker, Enrique Macias

Abstract: The composition of a forming planet is set by the material it accretes from its parent protoplanetary disk. Therefore, it is crucial to map the chemical make-up of the gas in disks to understand the chemical environment of planet formation. This paper presents molecular line observations taken with the Atacama Large Millimeter/submillimeter Array of the planet-hosting disk around the young star HD… ▽ More The composition of a forming planet is set by the material it accretes from its parent protoplanetary disk. Therefore, it is crucial to map the chemical make-up of the gas in disks to understand the chemical environment of planet formation. This paper presents molecular line observations taken with the Atacama Large Millimeter/submillimeter Array of the planet-hosting disk around the young star HD 169142. We detect N2H+, CH3OH, [CI], DCN, CS, C34S, 13CS, H2CS, H2CO, HC3N and c-C3H2 in this system for the first time. Combining these data with the recent detection of SO and previously published DCO+ data, we estimate the location of H2O and CO snowlines and investigate radial variations in the gas phase C/O ratio. We find that the HD 169142 disk has a relatively low N2H+ flux compared to the disks around Herbig stars HD 163296 and MWC 480 indicating less CO freeze-out and place the CO snowline beyond the millimetre disk at ~150 au. The detection of CH3OH from the inner disk is consistent with the H2O snowline being located at the edge of the central dust cavity at ~20 au. The radially varying CS/SO ratio across the proposed H2O snowline location is consistent with this interpretation. Additionally, the detection of CH3OH in such a warm disk adds to the growing evidence supporting the inheritance of complex ices in disks from the earlier, colder stages of star formation. Finally, we propose that the giant HD 169142 b located at 37 au is forming between the CO2 and H2O snowlines where the local elemental make of the gas is expected to have C/O=1.0. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: Accepted A&A 13th August 2023

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

SO and SiS Emission Tracing an Embedded Planet and Compact $^{12}$CO and $^{13}$CO Counterparts in the HD 169142 Disk

Authors: Charles J. Law, Alice S. Booth, Karin I. Öberg

Abstract: Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in mol… ▽ More Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in molecular gas observations. Using high-angular-resolution ALMA archival data of the HD 169142 disk, we identify compact SO J=8$_8$-7$_7$ and SiS J=19-18 emission coincident with the position of a ${\sim}$2 M$_{\rm{Jup}}$ planet seen as a localized, Keplerian NIR feature within a gas-depleted, annular dust gap at ${\approx}$38 au. The SiS emission is located along an azimuthal arc and has a similar morphology as a known $^{12}$CO kinematic excess. This is the first tentative detection of SiS emission in a protoplanetary disk and suggests that the planet is driving sufficiently strong shocks to produce gas-phase SiS. We also report the discovery of compact $^{12}$CO and $^{13}$CO J=3-2 emission coincident with the planet location. Taken together, a planet-driven outflow provides the best explanation for the properties of the observed chemical asymmetries. We also resolve a bright, azimuthally-asymmetric SO ring at ${\approx}$24 au. While most of this SO emission originates from ice sublimation, its asymmetric distribution implies azimuthal temperature variations driven by a misaligned inner disk or planet-disk interactions. Overall, the HD 169142 disk shows several distinct chemical signatures related to giant planet formation and presents a powerful template for future searches of planet-related chemical asymmetries in protoplanetary disks. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 22 pages, 12 figures, accepted for publication in ApJL

Deprojecting and constraining the vertical thickness of exoKuiper belts

Authors: James Terrill, Sebastian Marino, Richard A. Booth, Yinuo Han, Jeff Jennings, Mark C. Wyatt

Abstract: Constraining the vertical and radial structure of debris discs is crucial to understanding their formation, evolution and dynamics. To measure both the radial and vertical structure, a disc must be sufficiently inclined. However, if a disc is too close to edge-on, deprojecting its emission becomes non-trivial. In this paper we show how Frankenstein, a non-parametric tool to extract the radial brig… ▽ More Constraining the vertical and radial structure of debris discs is crucial to understanding their formation, evolution and dynamics. To measure both the radial and vertical structure, a disc must be sufficiently inclined. However, if a disc is too close to edge-on, deprojecting its emission becomes non-trivial. In this paper we show how Frankenstein, a non-parametric tool to extract the radial brightness profile of circumstellar discs, can be used to deproject their emission at any inclination as long as they are optically thin and axisymmetric. Furthermore, we extend Frankenstein to account for the vertical thickness of an optically thin disc ($H(r)$) and show how it can be constrained by sampling its posterior probability distribution and assuming a functional form (e.g. constant $h=H/r$), while fitting the radial profile non-parametrically. We use this new method to determine the radial and vertical structure of 16 highly inclined debris discs observed by ALMA. We find a wide range of vertical aspect ratios, $h$, ranging from $0.020\pm0.002$ (AU Mic) to $0.20\pm0.03$ (HD 110058), which are consistent with parametric models. We find a tentative correlation between $h$ and the disc fractional width, as expected if wide discs were more stirred. Assuming discs are self-stirred, the thinnest discs would require the presence of at least 500 km-sized planetesimals. The thickest discs would likely require the presence of planets. We also recover previously inferred and new radial structures, including a potential gap in the radial distribution of HD 61005. Finally, our new extension of Frankenstein also allows constraining how $h$ varies as a function of radius, which we test on 49 Ceti, finding that $h$ is consistent with being constant. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS. 17 pages. 16 figures

Investigating the asymmetric chemistry in the disk around the young star HD 142527

Authors: Milou Temmink, Alice S. Booth, Nienke van der Marel, Ewine F. van Dishoeck

Abstract: The atmospheric composition of planets is determined by the chemistry of the disks in which they form. Studying the gas-phase molecular composition of disks thus allows us to infer what the atmospheric composition of forming planets might be. Recent observations of the IRS 48 disk have shown that (asymmetric) dust traps can directly impact the observable chemistry, through radial and vertical tran… ▽ More The atmospheric composition of planets is determined by the chemistry of the disks in which they form. Studying the gas-phase molecular composition of disks thus allows us to infer what the atmospheric composition of forming planets might be. Recent observations of the IRS 48 disk have shown that (asymmetric) dust traps can directly impact the observable chemistry, through radial and vertical transport, and the sublimation of ices. The asymmetric HD 142527 disk provides another good opportunity to investigate the role of dust traps in setting the disk's chemical composition. In this work, we use archival ALMA observations of the HD 142527 disk to obtain an as large as possible molecular inventory, which allows us to investigate the possible influence of the asymmetric dust trap on the disk's chemistry. We present the first ALMA detections of [C I], 13C18O, DCO+, H2CO and additional transition of HCO+ and CS in this disk. In addition, we have acquired upper limits for non-detected species such as SO and CH3OH. For the majority of the observed molecules, a decrement in the emission at the location of the dust trap is found. For the main CO isotopologues continuum over-subtraction likely causes the observed asymmetry, while for CS and HCN we propose that the observed asymmetries are likely due to shadows cast by the misaligned inner disk. As the emission of the observed molecules is not co-spatial with the dust trap and no SO or CH3OH are found, thermal sublimation of icy mantles does not appear to play a major role in changing the gas-phase composition of the outer disk in HD 142527 disk. Using our observations of 13C18O and DCO+ and a RADMC-3D model, we determine the CO snowline to be located beyond the dust traps, favouring cold gas-phase formation of H2CO, rather than the hydrogenation of CO-ice and subsequent sublimation. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A on 12/04/2023

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

A major asymmetric ice trap in a planet-forming disk IV. Nitric oxide gas and a lack of CN tracing sublimating ices and a C/O ratio $< 1$

Authors: M. Leemker, A. S. Booth, E. F. van Dishoeck, N. van der Marel, B. Tabone, N. F. W. Ligterink, N. G. C. Brunken, M. R. Hogerheijde

Abstract: [Abridged] Most well-resolved disks observed with ALMA show signs of dust traps. These dust traps set the chemical composition of the planet forming material in these disks, as the dust grains with their icy mantles are trapped at specific radii and could deplete the gas and dust of volatiles at smaller radii. In this work we analyse the first detection of nitric oxide (NO) in a protoplanetary dis… ▽ More [Abridged] Most well-resolved disks observed with ALMA show signs of dust traps. These dust traps set the chemical composition of the planet forming material in these disks, as the dust grains with their icy mantles are trapped at specific radii and could deplete the gas and dust of volatiles at smaller radii. In this work we analyse the first detection of nitric oxide (NO) in a protoplanetary disk. We aim to constrain the nitrogen chemistry and the gas-phase C/O ratio in the highly asymmetric dust trap in the Oph-IRS 48 disk. We use ALMA observations of NO, CN, C$_2$H, and related molecules and model the effect of the dust trap on the physical and chemical structure using the thermochemical code DALI. Furthermore, we explore how ice sublimation contributes to the observed emission lines. NO is only observed at the location of the dust trap but CN and C$_2$H are not detected in the Oph-IRS 48 disk. This results in an CN/NO column density ratio of $< 0.05$ and thus a low C/O ratio at the location of the dust trap. The main gas-phase formation pathways to NO through OH and NH in the fiducial model predict NO emission that is an order of magnitude lower than is observed. The gaseous NO column density can be increased by factors ranging from 2.8 to 10 when the H$_2$O and NH$_3$ gas abundances are significantly boosted by ice sublimation. However, these models are inconsistent with the upper limits on the H$_2$O and OH column densities derived from observations. We propose that the NO emission in the Oph-IRS 48 disk is closely related to the nitrogen containing ices sublimating in the dust trap. The non-detection of CN constrains the C/O ratio both inside and outside the dust trap to be $< 1$ if all nitrogen initially starts as N$_2$ and $\leq 0.6$, consistent with the Solar value, if (part of) the nitrogen initially starts as N or NH$_3$. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 673, A7 (2023)

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)

A New Distance to the Supernova Remnant DA 530 Based on HI Absorption of Polarized Emission

Authors: Rebecca A. Booth, Roland Kothes, Tom Landecker, Jo-Anne Brown, Andrew Gray, Tyler Foster, Eric Greisen

Abstract: Supernova remnants (SNRs) are significant contributors of matter and energy to the interstellar medium. Understanding the impact and the mechanism of this contribution requires knowledge of the physical size, energy, and expansion rate of individual SNRs, which can only come if reliable distances can be obtained. We aim to determine the distance to the SNR DA 530 (G93.3+6.9), an object of low surf… ▽ More Supernova remnants (SNRs) are significant contributors of matter and energy to the interstellar medium. Understanding the impact and the mechanism of this contribution requires knowledge of the physical size, energy, and expansion rate of individual SNRs, which can only come if reliable distances can be obtained. We aim to determine the distance to the SNR DA 530 (G93.3+6.9), an object of low surface brightness. To achieve this, we used the Dominion Radio Astrophysical Observatory Synthesis Telescope and the National Radio Astronomy Observatory Very Large Array to observe the absorption by intervening HI of the polarized emission from DA 530. Significant absorption was detected at velocities $-28$ and -67 km/s (relative to the local standard of rest), corresponding to distances of 4.4 and 8.3 kpc, respectively. Based on the radio and X-ray characteristics of DA 530, we conclude that the minimum distance is 4.4$^{+0.4}_{-0.2}$ kpc. At this minimum distance, the diameter of the SNR is 34$^{+4}_{-1}$ pc, and the elevation above the Galactic plane is 537$^{+40}_{-32}$ pc. The $-67$ km/s absorption likely occurs in gas whose velocity is not determined by Galactic rotation. We present a new data processing method for combining Stokes $Q$ and $U$ observations of the emission from an SNR into a single HI absorption spectrum, which avoids the difficulties of the noise-bias subtraction required for the calculation of polarized intensity. The polarized absorption technique can be applied to determine distances to many more SNRs. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Disentangling protoplanetary disk gas mass and carbon depletion through combined atomic and molecular tracers

Authors: J. A. Sturm, A. S. Booth, M. K. McClure, M. Leemker, E. F. van Dishoeck

Abstract: The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$C… ▽ More The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$CO, C$^{18}$O and optically thin C$^{17}$O $J$=2-1 lines, and use additional high angular resolution Atacama Large Millimeter Array millimeter continuum and CO data to construct a representative model of LkCa 15. The transitions that constrain the gas mass and carbon abundance most are C$^{17}$O 2-1, N${_2}$H$^+$ 3-2 and HD 1-0. Using these three molecules we find that the gas mass in the LkCa 15 disk is $M_\mathrm{g}=0.01 ^{+0.01}_{-0.004} M_{\odot}$, a factor of six lower than estimated before. The carbon abundance is C/H = ($3 \pm 1.5) \times10^{-5}$, implying a moderate depletion of elemental carbon by a factor of 3-9. All other analyzed transitions also agree with these numbers, within a modeling uncertainty of a factor of two. Using the resolved \ce{C2H} image we find a C/O ratio of $\sim$1, which is consistent with literature values of H$_2$O depletion in this disk. The lack of severe carbon depletion in the LkCa 15 disk is consistent with the young age of the disk, but contrasts with the higher depletions seen in older cold transition disks. Combining optically thin CO isotopologue lines with N$_2$H$^+$ is promising to break the degeneracy between gas mass and CO abundance. The moderate level of depletion for this source with a cold, but young disk, suggests that long carbon transformation timescales contribute to the evolutionary trend seen in the level of carbon depletion among disk populations, rather than evolving temperature effects and presence of dust traps alone. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 16 pages, 8 figures. Accepted for publication in A&A

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/

Extreme Pebble Accretion in Ringed Protoplanetary Discs

Authors: Daniel P. Cummins, James E. Owen, Richard A. Booth

Abstract: Axisymmetric dust rings containing tens to hundreds of Earth masses of solids have been observed in protoplanetary discs with (sub-)millimetre imaging. Here, we investigate the growth of a planetary embryo in a massive (150M$_\oplus$) axisymmetric dust trap through dust and gas hydrodynamics simulations. When accounting for the accretion luminosity of the planetary embryo from pebble accretion, th… ▽ More Axisymmetric dust rings containing tens to hundreds of Earth masses of solids have been observed in protoplanetary discs with (sub-)millimetre imaging. Here, we investigate the growth of a planetary embryo in a massive (150M$_\oplus$) axisymmetric dust trap through dust and gas hydrodynamics simulations. When accounting for the accretion luminosity of the planetary embryo from pebble accretion, the thermal feedback on the surrounding gas leads to the formation of an anticyclonic vortex. Since the vortex forms at the location of the planet, this has significant consequences for the planet's growth: as dust drifts towards the pressure maximum at the centre of the vortex, which is initially co-located with the planet, a rapid accretion rate is achieved, in a distinct phase of ``vortex-assisted'' pebble accretion. Once the vortex separates from the planet due to interactions with the disc, it accumulates dust, shutting off accretion onto the planet. We find that this rapid accretion, mediated by the vortex, results in a planet containing $\approx$ 100M$_\oplus$ of solids. We follow the evolution of the vortex, as well as the efficiency with which dust grains accumulate at its pressure maximum as a function of their size, and investigate the consequences this has for the growth of the planet as well as the morphology of the protoplanetary disc. We speculate that this extreme formation scenario may be the origin of giant planets which are identified to be significantly enhanced in heavy elements. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 20 pages, 14 figures. Accepted for publication in MNRAS. Supplementary material can be found at https://doi.org/10.5281/zenodo.6761702

Super-resolution trends in the ALMA Taurus survey: Structured inner discs and compact discs

Authors: Jeff Jennings, Marco Tazzari, Cathie J. Clarke, Richard A. Booth, Giovanni P. Rosotti

Abstract: The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank')… ▽ More The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank') to identify new substructure in three compact and seven extended sources. Among the new features we identify three trends: a higher occurrence rate of substructure in the survey's compact discs than previously seen, underresolved (potentially azimuthally asymmetric) substructure in the innermost disc of extended sources, and a 'shoulder' on the trailing edge of a ring in discs with strong depletion at small radii. Noting the shoulder morphology is present in multiple discs observed at higher resolution, we postulate it is tracing a common physical mechanism. We further demonstrate how a super-resolution frank brightness profile is useful in motivating an accurate parametric model, using the highly structured source DL Tau in which frank finds two new rings. Finally we show that sparse (u, v) plane sampling may be masking the presence of substructure in several additional compact survey sources. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted to MNRAS

Inward and outward migration of massive planets: moving towards a stalling radius

Authors: C. E. Scardoni, C. J. Clarke, G. P. Rosotti, R. A. Booth, R. D. Alexander, E. Ragusa

Abstract: Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simu… ▽ More Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simulations exploring a range of $K'$ and disc mass values to test and extend these results. The excitation of planet eccentricity in live-planet simulations breaks the direct dependence of migration rate (rate of change of semi-major axis) on the torques imposed, an effect that `fixed-planet' simulations cannot treat. By disentangling the contribution to the torque due to the semi-major axis evolution from that due to the eccentricity evolution, we recover the relation between the magnitude and sign of migration and $K'$ and argue that this relation may be better expressed in terms of the related gap depth parameter $K$. We present a toy model in which the sign of planetary migration changes at a limiting value of $K$, through which we explore planets' migration in viscously evolving discs. The existence of the torque reversal shapes the planetary system's architecture by accumulating planets either at the stalling radius or in a band around it (defined by the interplay between the planet migration and the disc evolution). In either case, planets pile up in the area $1-10$ au, disfavouring hot Jupiter formation through Type II migration in the planet-dominated regime. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

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

Dust formation in the outflows of catastrophically evaporating planets

Authors: Richard A. Booth, James E. Owen, Matthäus Schulik

Abstract: Ultra-short period planets offer a window into the poorly understood interior composition of exoplanets through material evaporated from their rocky interiors. Among these objects are a class of disintegrating planets, observed when their dusty tails transit in front of their host stars. These dusty tails are thought to originate from dust condensation in thermally-driven winds emanating from the… ▽ More Ultra-short period planets offer a window into the poorly understood interior composition of exoplanets through material evaporated from their rocky interiors. Among these objects are a class of disintegrating planets, observed when their dusty tails transit in front of their host stars. These dusty tails are thought to originate from dust condensation in thermally-driven winds emanating from the sublimating surfaces of these planets. Existing models of these winds have been unable to explain their highly variable nature and have not explicitly modelled how dust forms in the wind. Here we present new radiation-hydrodynamic simulations of the winds from these planets, including a minimal model for the formation and destruction of dust, assuming that nucleation can readily take place. We find that dust forms readily in the winds, a consequence of large dust grains obtaining lower temperatures than the planet's surface. As hypothesised previously, we find that the coupling of the planet's surface temperature to the outflow properties via the dust's opacity can drive time-variable flows when dust condensation is sufficiently fast. In agreement with previous work, our models suggest that these dusty tails are a signature of catastrophically evaporating planets that are close to the end of their lives. Finally, we discuss the implications of our results for the dust's composition. More detailed hydrodynamic models that self-consistently compute the nucleation and composition of the dust and gas are warranted in order to use these models to study the planet's interior composition. △ Less

Submitted 26 October, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: 15 pages. Authors accepted version, to be published in MNRAS

Iceline Variations Driven by Protoplanetary Disc Gaps

Authors: Madelyn Broome, Oliver Shorttle, Mihkel Kama, Richard A. Booth

Abstract: The composition of forming planets is strongly affected by the protoplanetary disc's thermal structure. This thermal structure is predominantly set by dust radiative transfer and viscous (accretional) heating and can be impacted by gaps - regions of low dust and gas density that can occur when planets form. The effect of variations in dust surface density on disc temperature has been poorly unders… ▽ More The composition of forming planets is strongly affected by the protoplanetary disc's thermal structure. This thermal structure is predominantly set by dust radiative transfer and viscous (accretional) heating and can be impacted by gaps - regions of low dust and gas density that can occur when planets form. The effect of variations in dust surface density on disc temperature has been poorly understood until now. In this work, we use the radiative transfer code MCMax to model the 2D dust thermal structure with individual gaps corresponding to planets with masses of 0.1 M$_J$ - 5 M$_J$ and orbital radii of 3, 5, and 10 AU. Low dust opacity in the gap allows radiation to penetrate deeper into the disc and warm the midplane by up to 16 K, but only for gaps located in the region of the disc where stellar irradiation is the dominant source of heat (here, a$\gtrsim$4 AU). In viscously-heated regions (a$\lesssim$4 AU), the midplane of the gap is relatively cooler by up to 100 K. Outside of the gap, broad radial oscillations in heating and cooling are present due to changes in disc flaring. These thermal features affect local segregation of volatile elements (H$_2$O, CH$_4$, CO2, CO) between the dust and gas. We find that icelines experience dramatic shifts relative to gapless models: up to 6.5 AU towards the star and 4.3 AU towards the midplane. While quantitative predictions of iceline deviations will require more sophisticated models which include transport and sublimation/condensation kinetics, our results provide evidence that planet-induced iceline variations represent a potential feedback from the planet onto the composition of material it is accreting. △ Less

Submitted 14 September, 2023; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 11 pages, 5 figures. Accepted to MNRAS, published 20 April 2023

Journal ref: MNRAS, 522, 3378 (2023)

Gas temperature structure across transition disk cavities

Authors: M. Leemker, A. S. Booth, E. F. van Dishoeck, A. F. Pérez-Sánchez, J. Szulágyi, A. D. Bosman, S. Bruderer, S. Facchini, M. R. Hogerheijde, T. Paneque-Carreño, J. A. Sturm

Abstract: [Abridged] Most disks observed at high angular resolution show substructures. Knowledge about the gas surface density and temperature is essential to understand these. The aim of this work is to constrain the gas temperature and surface density in two transition disks: LkCa15 and HD 169142. We use new ALMA observations of the $^{13}$CO $J=6-5$ transition together with archival $J=2-1$ data of… ▽ More [Abridged] Most disks observed at high angular resolution show substructures. Knowledge about the gas surface density and temperature is essential to understand these. The aim of this work is to constrain the gas temperature and surface density in two transition disks: LkCa15 and HD 169142. We use new ALMA observations of the $^{13}$CO $J=6-5$ transition together with archival $J=2-1$ data of $^{12}$CO, $^{13}$CO and C$^{18}$O to observationally constrain the gas temperature and surface density. Furthermore, we use the thermochemical code DALI to model the temperature and density structure of a typical transition disk. The $6-5/2-1$ line ratio in LkCa15 constrains the gas temperature in the emitting layers inside the dust cavity to be up to 65 K, warmer than in the outer disk at 20-30 K. For the HD 169142, the peak brightness temperature constrains the gas in the dust cavity of HD 169142 to be 170 K, whereas that in the outer disk is only 100 K. Models also show that a more luminous central star, a lower abundance of PAHs and the absence of a dusty inner disk increase the temperature of the emitting layers and hence the line ratio in the gas cavity. The gas column density in the LkCa15 dust cavity drops by a factor >2 compared to the outer disk, with an additional drop of an order of magnitude inside the gas cavity at 10 AU. In the case of HD 169142, the gas column density drops by a factor of 200$-$500 inside the gas cavity, which could be due to a massive companion of several M$_{\mathrm{J}}$. The broad dust-depleted gas region from 10-68 AU for LkCa15 may imply several lower mass planets. This work demonstrates that knowledge of the gas temperature is important to determine the gas surface density and thus whether planets, and if so what kind of planets, are the most likely carving the dust cavities. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Astronomy and astrophysics

Journal ref: A&A 663, A23 (2022)

A major asymmetric ice trap in a planet-forming disk: III. First detection of dimethyl ether

Authors: Nashanty G. C. Brunken, Alice S. Booth, Margot Leemker, Pooneh Nazari, Nienke van der Marel, Ewine F. van Dishoeck

Abstract: The complex organic molecules (COMs) detected in star-forming regions are the precursors of the prebiotic molecules that can lead to the emergence of life. By studying COMs in more evolved protoplanetary disks we can gain a better understanding of how they are incorporated into planets. This paper presents ALMA band 7 observations of the dust and ice trap in the protoplanetary disk around Oph IRS… ▽ More The complex organic molecules (COMs) detected in star-forming regions are the precursors of the prebiotic molecules that can lead to the emergence of life. By studying COMs in more evolved protoplanetary disks we can gain a better understanding of how they are incorporated into planets. This paper presents ALMA band 7 observations of the dust and ice trap in the protoplanetary disk around Oph IRS 48. We report the first detection of dimethyl ether (CH3OCH3) in a planet-forming disk and a tentative detection of methyl formate (CH3OCHO). We determined column densities for the detected molecules and upper limits on non-detected species using the CASSIS spectral analysis tool. The inferred column densities of CH3OCH3 and CH3OCHO with respect to methanol (CH3OH) are of order unity, indicating unusually high abundances of these species compared to other environments. Alternatively, the 12CH3OH emission is optically thick and beam diluted, implying a higher CH3OH column density and a smaller emitting area than originally thought. The presence of these complex molecules can be explained by thermal ice sublimation, where the dust cavity edge is heated by irradiation and the full volatile ice content is observable in the gas phase. This work confirms the presence of oxygen-bearing molecules more complex than CH3OH in protoplanetary disks for the first time. It also shows that it is indeed possible to trace the full interstellar journey of COMs across the different evolutionary stages of star, disk, and planet formation. △ Less

Submitted 6 March, 2022; originally announced March 2022.

Comments: Accepted to A&A - 20th Jan 2022

Observing planet-driven dust spirals with ALMA

Authors: Jessica Speedie, Richard A. Booth, Ruobing Dong

Abstract: ALMA continuum observations of thermal emission from the dust component of protoplanetary disks have revealed an abundance of substructures that may be interpreted as evidence for embedded planets, but planet-driven spiral arms -- perhaps one of the most compelling lines of evidence -- have proven comparatively elusive. In this work, we test the capabilities of ALMA to detect the planet-driven spi… ▽ More ALMA continuum observations of thermal emission from the dust component of protoplanetary disks have revealed an abundance of substructures that may be interpreted as evidence for embedded planets, but planet-driven spiral arms -- perhaps one of the most compelling lines of evidence -- have proven comparatively elusive. In this work, we test the capabilities of ALMA to detect the planet-driven spiral signal in continuum emission. Carrying out hydrodynamic simulations and radiative transfer calculations, we present synthetic Band 7 continuum images for a wide range of disk and observing conditions. We show that thermal mass planets at tens of au typically drive spirals detectable within a few hours of integration time, and the detectable planet mass may be as low as $\sim$Neptune mass ($0.3 \, M_{\rm th}$). The grains probed by ALMA form spirals morphologically identical to the underlying gas spiral. The temperature of the dust spiral is crucial in determining its contrast, and spirals are easier to detect in disks with an adiabatic equation of state and longer cooling times. Resolving the spiral is not necessary for its detection; with the help of residual maps, the optimal beam size is a few times the spiral width at a constant noise level. Finally, we show how the presence of gaps and rings can impair our ability to recognize co-located spirals. Our work demonstrates the planet-finding potential of the current design specification of ALMA, and suggests that observing capability is not the bottleneck in searching for spirals induced by thermal mass planets. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 11 figures, 20 pages. Accepted for publication in The Astrophysical Journal. Image repository available at https://doi.org/10.6084/m9.figshare.19148912

Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

Authors: Jessica Speedie, Ralph E. Pudritz, Alex J. Cridland, Farzana Meru, Richard A. Booth

Abstract: ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed thro… ▽ More ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed throughout the disk, assuming either a high or low turbulent $α$ viscosity. We have discovered that there is a bifurcation in the migration history of forming planets as a consequence of varying the disk viscosity. In our high viscosity disks, inward migration prevails and yields compact planetary systems, tempered only by planet trapping at the water iceline around 5 au. In our lower viscosity models however, low mass planets can migrate outward to twice their initial orbital radii, driven by a radially extended region of strong outward-directed corotation torques located near the heat transition (where radiative heating of the disk by the star is comparable to viscous heating) - before eventually migrating inwards. We derive analytic expressions for the planet mass at which the corotation torque dominates, and find that this "corotation mass" scales as $M_{\rm p, corot} \sim α^{2/3}$. If disk winds dominate the corotation torque, the corotation mass scales linearly with wind strength. We propose that the observed bifurcation in disk demographics into a majority of compact dust disks and a minority of extended ring/gap systems is a consequence of a distribution of viscosity across the disk population. △ Less

Submitted 4 February, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: 25 pages, 15 figures, 10 movies (links embedded). Published in MNRAS

Journal ref: MNRAS, Volume 510, Issue 4, March 2022, Pages 6059-6084

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). 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