Showing 1–50 of 266 results for author: Andrews, S

Constraints on the gas-phase C/O ratio of DR Tau's outer disk from CS, SO, and C$_2$H observations

Authors: Jane Huang, Edwin A. Bergin, Romane Le Gal, Sean M. Andrews, Jaehan Bae, Luke Keyte, J. A. Sturm

Abstract: Millimeter wavelength observations of Class II protoplanetary disks often display strong emission from hydrocarbons and high CS/SO values, providing evidence that the gas-phase C/O ratio commonly exceeds 1 in their outer regions. We present new NOEMA observations of CS $5-4$, SO $7_6-6_5$ and $5_6-4_5$, C$_2$H $N=3-2$, HCN $3-2$, HCO$^+$ $3-2$, and H$^{13}$CO$^+$ $3-2$ in the DR Tau protoplanetary… ▽ More Millimeter wavelength observations of Class II protoplanetary disks often display strong emission from hydrocarbons and high CS/SO values, providing evidence that the gas-phase C/O ratio commonly exceeds 1 in their outer regions. We present new NOEMA observations of CS $5-4$, SO $7_6-6_5$ and $5_6-4_5$, C$_2$H $N=3-2$, HCN $3-2$, HCO$^+$ $3-2$, and H$^{13}$CO$^+$ $3-2$ in the DR Tau protoplanetary disk at a resolution of $\sim0.4''$ (80 au). Estimates for the disk-averaged CS/SO ratio range from $\sim0.4-0.5$, the lowest value reported thus far for a T Tauri disk. At a projected separation of $\sim180$ au northeast of the star, the SO moment maps exhibit a clump that has no counterpart in the other lines, and the CS/SO value decreases to $<0.2$ at its location. Thermochemical models calculated with DALI indicate that DR Tau's low CS/SO ratio and faint C$_2$H emission can be explained by a gas-phase C/O ratio that is $<1$ at the disk radii traced by NOEMA. Comparisons of DR Tau's SO emission to maps of extended structures traced by $^{13}$CO suggest that late infall may contribute to driving down the gas-phase C/O ratio of its disk. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: Accepted by ApJ

Exploring the Complex Ionization Environment of the Turbulent DM Tau Disk

Authors: Deryl E. Long, L. Ilsedore Cleeves, Fred C. Adams, Sean Andrews, Edwin A. Bergin, Viviana V. Guzmán, Jane Huang, A. Meredith Hughes, Chunhua Qi, Kamber Schwarz, Jacob B. Simon, David Wilner

Abstract: Ionization drives important chemical and dynamical processes within protoplanetary disks, including the formation of organics and water in the cold midplane and the transportation of material via accretion and magneto-hydrodynamic (MHD) flows. Understanding these ionization-driven processes is crucial for understanding disk evolution and planet formation. We use new and archival ALMA observations… ▽ More Ionization drives important chemical and dynamical processes within protoplanetary disks, including the formation of organics and water in the cold midplane and the transportation of material via accretion and magneto-hydrodynamic (MHD) flows. Understanding these ionization-driven processes is crucial for understanding disk evolution and planet formation. We use new and archival ALMA observations of HCO+, H13CO+, and N2H+ to produce the first forward-modeled 2D ionization constraints for the DM Tau protoplanetary disk. We include ionization from multiple sources and explore the disk chemistry under a range of ionizing conditions. Abundances from our 2D chemical models are post-processed using non-LTE radiative transfer, visibility sampling, and imaging, and are compared directly to the observed radial emission profiles. The observations are best fit by a modestly reduced CR ionization rate ($ζ_{CR}$ ~ 10$^{-18}$ s$^{-1}$) and a hard X-ray spectrum (hardness ratio [HR] = 0.3), which we associate with stellar flaring conditions. Our best-fit model under-produces emission in the inner disk, suggesting that there may be an additional mechanism enhancing ionization in DM Tau's inner disk. Overall, our findings highlight the complexity of ionization in protoplanetary disks and the need for high resolution multi-line studies. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, accepted to be published in The Astrophysical Journal (June 25, 2024)

A Dust-Trapping Ring in the Planet-Hosting Disk of Elias 2-24

Authors: Adolfo S. Carvalho, Laura M. Perez, Anibal Sierra, Maria Jesus Mellado, Lynne A. Hillenbrand, Sean Andrews, Myriam Benisty, Tilman Birnstiel, John M. Carpenter, Viviana V. Guzman, Jane Huang, Andrea Isella, Nicolas Kurtovic, Luca Ricci, David J. Wilner

Abstract: Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid gr… ▽ More Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid growth and possibly planetesimal production via the streaming instability. We present evidence of dust trapping in the bright ring of the planet-hosting disk Elias 2-24, from the analysis of 1.3 mm and 3 mm ALMA observations at high spatial resolution (0.029 arcsec, 4.0 au). We leverage the high spatial resolution to demonstrate that larger grains are more efficiently trapped and place constraints on the local turbulence ($8 \times 10^{-4} < α_\mathrm{turb} < 0.03$) and the gas-to-dust ratio ($Σ_g / Σ_d < 30$) in the ring. Using a scattering-included marginal probability analysis we measure a total dust disk mass of $M_\mathrm{dust} = 13.8^{+0.7}_{-0.5} \times 10^{-4} \ M_\odot$. We also show that at the orbital radius of the proposed perturber, the gap is cleared of material down to a flux contrast of 10$^{-3}$ of the peak flux in the disk. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, accepted to ApJ

Evidence for Non-zero Turbulence in the Protoplanetary disc around IM Lup

Authors: Kevin Flaherty, A. Meredith Hughes, Jacob B. Simon, Alicia Smith Reina, Chunhua Qi, Xue-Ning Bai, Sean M. Andrews, David J. Wilner, Agnes Kospal

Abstract: The amount of turbulence in protoplanetary discs around young stars is critical for determining the efficiency, timeline, and outcomes of planet formation. It is also difficult to measure. Observations are still limited, but direct measurements of the non-thermal, turbulent gas motion are possible with the Atacama Large Millimeter/submillimeter Array (ALMA). Using CO(2-1)/$^{13}$CO(2-1)/C$^{18}$O(… ▽ More The amount of turbulence in protoplanetary discs around young stars is critical for determining the efficiency, timeline, and outcomes of planet formation. It is also difficult to measure. Observations are still limited, but direct measurements of the non-thermal, turbulent gas motion are possible with the Atacama Large Millimeter/submillimeter Array (ALMA). Using CO(2-1)/$^{13}$CO(2-1)/C$^{18}$O(2-1) ALMA observations of the disc around IM Lup at ~0.4" (~60 au) resolution we find evidence of significant turbulence, at the level of $δv_{\rm turb}=(0.18-0.30)$c$_s$. This result is robust against systematic uncertainties (e.g., amplitude flux calibration, midplane gas temperature, disc self-gravity). We find that gravito-turbulence as the source of the gas motion is unlikely based on the lack of an imprint on the rotation curve from a massive disc, while magneto-rotational instabilities and hydrodynamic instabilities are still possible, depending on the unknown magnetic field strength and the cooling timescale in the outer disc. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted by MNRAS, 17 pages, 12 figures

SMA 200-400 GHz Survey for Dust Properties in the Icy Class II Disks in the Taurus Molecular Cloud

Authors: Chia-Ying Chung, Sean M. Andrews, Mark A. Gurwell, Melvyn Wright, Feng Long, Wenrui Xu, Hauyu Baobab Liu

Abstract: We present a new SMA survey of 47 Class II sources in the Taurus-Auriga region. Our observations made 12 independent samples of flux densities over the 200-400 GHz frequency range. We tightly constrained the spectral indices of most sources to a narrow range of $2.0\pm0.2$; only a handful of spatially resolved (e.g., diameter $>$250 au) disks present larger spectral indices. The simplest interpret… ▽ More We present a new SMA survey of 47 Class II sources in the Taurus-Auriga region. Our observations made 12 independent samples of flux densities over the 200-400 GHz frequency range. We tightly constrained the spectral indices of most sources to a narrow range of $2.0\pm0.2$; only a handful of spatially resolved (e.g., diameter $>$250 au) disks present larger spectral indices. The simplest interpretation for this result is that the (sub)millimeter luminosities of all of the observed target sources are dominated by very optically thick (e.g., $τ\gtrsim$5) dust thermal emission. Some previous works that were based on the optically thin assumption thus might have underestimated optical depths by at least one order of magnitude. Assuming DSHARP dust opacities, this corresponds to underestimates of dust masses by a similar factor. Moreover, some population synthesis models show that to explain the observed, narrowly distributed spectral indices, the disks in our selected sample need to have very similar dust temperatures ($T_{\small{dust}}$). Given a specific assumption of median $T_{\small{dust}}$, the maximum grain sizes ($a_{\small{max}}$) can also be constrained, which is a few times smaller than 0.1 mm for $T_{\small{dust}}\sim$100 K and a few mm for $T_{\small{dust}}\sim$24 K. The results may indicate that dust grain growth outside the water snowline is limited by the bouncing/fragmentation barriers. In the Class II disks, the dust mass budget outside of the water snowline may be largely retained instead of being mostly consumed by planet formation. While Class II disks still possess sufficient dust masses to feed planet formation at a later time, it is unknown whether or not dust coagulation and planet formation can be efficient or natural outside of the water snowline. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: accepted by Astrophysical Journal Supplement

On Kinematic Measurements of Self-Gravity in Protoplanetary Disks

Authors: Sean M. Andrews, Richard Teague, Christopher P. Wirth, Jane Huang, Zhaohuan Zhu

Abstract: Using controlled injection and recovery experiments, we devised an analysis prescription to assess the quality of dynamical measurements of protoplanetary disk gas masses based on resolved (CO) spectral line data, given observational limitations (resolution, sampling, noise), measurement bias, and ambiguities in the geometry and physical conditions. With sufficient data quality, this approach perf… ▽ More Using controlled injection and recovery experiments, we devised an analysis prescription to assess the quality of dynamical measurements of protoplanetary disk gas masses based on resolved (CO) spectral line data, given observational limitations (resolution, sampling, noise), measurement bias, and ambiguities in the geometry and physical conditions. With sufficient data quality, this approach performed well for massive disks ($M_{\rm d}/M_\ast=0.1$): we inferred $M_{\rm d}$ posteriors that recovered the true values with little bias ($\lesssim$ 20%) and uncertainties within a factor of two (2$σ$). The gas surface density profiles for such cases are recovered with remarkable fidelity. Some experimentation indicates that this approach becomes insensitive when $M_{\rm d}/M_\ast\lesssim5$%, due primarily to degeneracies in the surface density profile parameters. Including multiple lines that probe different vertical layers, along with some improvements in the associated tools, might push that practical boundary down by another factor of $\sim$two in ideal scenarios. We also demonstrated this analysis approach using archival ALMA observations of the MWC 480 disk (Öberg et al. 2021): we measured $M_{\rm d}=0.13^{\: +0.04}_{\: -0.01} \: M_\odot$ (corresponding to $M_{\rm d}/M_\ast=7\pm1$%) and identified kinematic substructures consistent with surface density gaps around 65 and 135 au. Overall, this (and similar work) suggests that these dynamical measurements offer powerful new constraints with sufficient accuracy and precision to quantify gas masses and surface densities at the high end of the $M_{\rm d}/M_\ast$ distribution, and therefore can serve as key benchmarks for detailed thermo-chemical modeling. We address some prospects for improvements, and discuss various caveats and limitations to guide future work. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: ApJ, in press; 32 pages, 29 figures

Population Synthesis Models Indicate a Need for Early and Ubiquitous Disk Substructures

Authors: Luca Delussu, Tilman Birnstiel, Anna Miotello, Paola Pinilla, Giovanni Rosotti, Sean M. Andrews

Abstract: Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust… ▽ More Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport varying the initial conditions of the disk and substructure to find the best match to the observed distributions. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure, namely a perturbation strong enough to be able to trap particles, and that is formed early in the evolution of the disk, that is within 0.4Myr. Agreement is reached by relatively high initial disk masses ($10^{-2.3}M_{\star}\leqslant M_{disk}\leqslant10^{-0.5}M_{\star}$) and moderate levels of turbulence ($10^{-3.5}\leqslantα\leqslant 10^{-2.5}$). Other disk parameters play a weaker role. Only opacities with high absorption efficiency can reproduce the observed spectral indices. Our results extend to the whole population that substructure is likely ubiquitous, so far assessed only in individual disks and implies that most "smooth" disks hide unresolved substructure. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Planet formation regulated by galactic-scale interstellar turbulence

Authors: Andrew J. Winter, Myriam Benisty, Sean M. Andrews

Abstract: Planet formation occurs over a few Myr within protoplanetary discs of dust and gas, which are often assumed to evolve in isolation. However, extended gaseous structures have been uncovered around many protoplanetary discs, suggestive of late-stage in-fall from the interstellar medium (ISM). To quantify the prevalence of late-stage in-fall, we apply an excursion set formalism to track the local den… ▽ More Planet formation occurs over a few Myr within protoplanetary discs of dust and gas, which are often assumed to evolve in isolation. However, extended gaseous structures have been uncovered around many protoplanetary discs, suggestive of late-stage in-fall from the interstellar medium (ISM). To quantify the prevalence of late-stage in-fall, we apply an excursion set formalism to track the local density and relative velocity of the ISM over the disc lifetime. We then combine the theoretical Bondi-Hoyle-Lyttleton (BHL) accretion rate with a simple disc evolution model, anchoring stellar accretion time-scales to observational constraints. Disc lifetimes, as well as masses, stellar accretion rates and approximate gas radii as a function of stellar mass and age are remarkably well-reproduced by our simple model that includes only ISM accretion. If BHL accretion is the origin of the initial disc material, then $20{-}70$ percent of discs are mostly composed of material accreted in the most recent half of their lifetime, suggesting that disc properties are not a direct test of isolated evolution models. The ISM also deposits sufficient energy to drive observed turbulence in protoplanetary discs with viscous $α_\mathrm{SS} \sim 10^{-5}- 10^{-1}$, although we emphasise that angular momentum transport may still be driven internal processes. Overall, our findings represent a compelling case for regulation of planet formation by galactic-scale turbulence, with broad consequences for planet formation theory. Our simple approach can be easily applied to semi-analytic models, and our results represent urgent motivation for an unbiased survey to confirm or refute our findings. △ Less

Submitted 18 June, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: Submitted to ApJL

Sites of Planet Formation in Binary Systems. I. Evidence for Disk-Orbit Alignment in the Close Binary FO Tau

Authors: Benjamin M. Tofflemire, Lisa Prato, Adam L. Kraus, Dominique Segura-Cox, G. H. Schaefer, Rachel Akeson, Sean Andrews, Eric L. N. Jensen, Christopher M. Johns-Krull, J. J. Zanazzi, M. Simon

Abstract: Close binary systems present challenges to planet formation. As binary separations decrease, so too do the occurrence rates of protoplanetary disks in young systems and planets in mature systems. For systems that do retain disks, their disk masses and sizes are altered by the presence of the binary companion. Through the study of protoplanetary disks in binary systems with known orbital parameters… ▽ More Close binary systems present challenges to planet formation. As binary separations decrease, so too do the occurrence rates of protoplanetary disks in young systems and planets in mature systems. For systems that do retain disks, their disk masses and sizes are altered by the presence of the binary companion. Through the study of protoplanetary disks in binary systems with known orbital parameters, we seek to determine the properties that promote disk retention and, therefore, planet formation. In this work, we characterize the young binary-disk system, FO Tau. We determine the first full orbital solution for the system, finding masses of $0.35^{+0.06}_{-0.05}\ M_\odot$ and $0.34\pm0.05\ M_\odot$ for the stellar components, a semi-major axis of $22(^{+2}_{-1})$ AU, and an eccentricity of $0.21(^{+0.04}_{-0.03})$. With long-baseline ALMA interferometry, we detect 1.3mm continuum and $^{12}{\mathrm{CO}} \ (J=2-1)$ line emission toward each of the binary components; no circumbinary emission is detected. The protoplanetary disks are compact, consistent with being truncated by the binary orbit. The dust disks are unresolved in the image plane and the more extended gas disks are only marginally resolved. Fitting the continuum and CO visibilities, we determine the inclination of each disk, finding evidence for alignment of the disk and binary orbital planes. This study is the first of its kind linking the properties of circumstellar protoplanetary disks to a precisely known binary orbit. In the case of FO Tau, we find a dynamically placid environment (coplanar, low eccentricity), which may foster its potential for planet formation. △ Less

Submitted 11 April, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: AJ accepted, 29 pages, 14 figures

High-resolution Pan-STARRS and SMA observations of IRAS 23077+6707: A giant edge-on protoplanetary disk

Authors: Kristina Monsch, Joshua B. Lovell, Ciprian T. Berghea, Gordian Edenhofer, Garrett K. Keating, Sean M. Andrews, Ammar Bayyari, Jeremy J. Drake, David J. Wilner

Abstract: We present resolved images of IRAS 23077+6707 ("Dracula's Chivito") in 1.3 mm/225 GHz thermal dust and CO gas emission with the Submillimeter Array (SMA) and optical (0.5-$0.8\,μ\mathrm{m}$) scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observe… ▽ More We present resolved images of IRAS 23077+6707 ("Dracula's Chivito") in 1.3 mm/225 GHz thermal dust and CO gas emission with the Submillimeter Array (SMA) and optical (0.5-$0.8\,μ\mathrm{m}$) scattered light with the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). The Pan-STARRS data show a bipolar distribution of optically scattering dust that is characteristic for disks observed at high inclinations. Its scattered light emission spans ${\sim}14''$, with two highly asymmetric filaments extending along the upper bounds of each nebula by ${\sim}9''$. The SMA data measure 1.3 mm continuum dust as well as $^{12}$CO, $^{13}$CO and C$^{18}$O $J$=2$-$1 line emission over $12''$-$14''$ extents, with the gas presenting the typical morphology of a disk in Keplerian rotation, in both position-velocity space and in each CO line spectrum. IRAS 23077+6707 has no reported distance estimate, but if it is located in the Cepheus star-forming region (180-800 pc), it would have a radius spanning thousands of au. Taken together, we infer IRAS 23077+6707 to be a giant and gas-rich edge-on protoplanetary disk, which to our knowledge is the largest in extent so far discovered. △ Less

Submitted 14 May, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: published in The Astrophysical Journal Letters 967:L2 (9pp), 2024 May 20

Journal ref: ApJL 967:L2 (9pp), 2024

High-resolution ALMA observations of compact discs in the wide-binary system Sz 65 and Sz 66

Authors: J. M. Miley, J. Carpenter, R. Booth, J. Jennings, T. J. Haworth, M. Vioque, S. Andrews, D. Wilner, M. Benisty, J. Huang, L. Perez, V. Guzman, L. Ricci, A. Isella

Abstract: Substructures in disc density are ubiquitous in the bright extended discs that are observed with high resolution. These substructures are intimately linked to the physical mechanisms driving planet formation and disc evolution. Surveys of star-forming regions find that most discs are in fact compact, less luminous, and do not exhibit these same substructures. It remains unclear whether compact dis… ▽ More Substructures in disc density are ubiquitous in the bright extended discs that are observed with high resolution. These substructures are intimately linked to the physical mechanisms driving planet formation and disc evolution. Surveys of star-forming regions find that most discs are in fact compact, less luminous, and do not exhibit these same substructures. It remains unclear whether compact discs also have similar substructures or if they are featureless. This suggests that different planet formation and disc evolution mechanisms operate in these discs. We investigated evidence of substructure within two compact discs around the stars Sz 65 and Sz 66 using high angular resolution observations with ALMA at 1.3 mm. The two stars form a wide-binary system with 6.36 arcsec separation. The continuum observations achieve a synthesised beam major axis of 0.026 arcsec, equivalent to about 4.0 au, enabling a search for substructure on these spatial scales and a characterisation of the gas and dust disc sizes with high precision. We analysed the data in the image plane through an analysis of reconstructed images, as well as in the uv plane by modelling the visibilities and by an analysis of the 12CO emission line. Comparisons were made with high-resolution observations of compact discs and radially extended discs. We find evidence of substructure in the dust distribution of Sz 65, namely a shallow gap centred at approximately 20 au, with an emission ring exterior to it. Ninety percent of the measured continuum flux is found within 27 au, and the distance for 12CO is 142 au. The observations show that Sz 66 is very compact: 90 per cent of the continuum flux is contained within 16 au, and 48 au for the gas. While the overall prevalence and diversity of substructure in compact discs relative to larger discs is yet to be determined, we find evidence that substructures can exist in compact discs. △ Less

Submitted 13 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Journal ref: A&A, 682, A55 (2024)

SMA detection of an extreme millimeter flare from the young class III star HD 283572

Authors: Joshua Bennett Lovell, Garrett K. Keating, David J. Wilner, Sean M. Andrews, Meredith MacGregor, Ramisa Akther Rahman, Ramprasad Rao, Jonathan P. Williams

Abstract: We present evidence of variable 1.3 millimeter emission from the 1-3 Myr, SpT G2-G5 class III YSO, HD~283572. HD~283572 was observed on 8 dates with the Submillimeter Array between 2021 December and 2023 May, a total on-source time of 10.2 hours, probing a range of timescales down to 5.2 seconds. Averaging all data obtained on 2022 Jan 17 shows a 4.4 mJy ($8.8σ$) point source detection with a nega… ▽ More We present evidence of variable 1.3 millimeter emission from the 1-3 Myr, SpT G2-G5 class III YSO, HD~283572. HD~283572 was observed on 8 dates with the Submillimeter Array between 2021 December and 2023 May, a total on-source time of 10.2 hours, probing a range of timescales down to 5.2 seconds. Averaging all data obtained on 2022 Jan 17 shows a 4.4 mJy ($8.8σ$) point source detection with a negative spectral index ($α{=}{-2.7}{\pm}1.2$), with peak emission rising to 13.8 mJy in one 3 minute span, and 25 mJy in one 29.7 second integration ($L_ν=4.7\times10^{17}$ erg s$^{-1}$ Hz$^{-1}$). Combining our data for the other 7 dates shows no detection, with an rms noise of 0.24 mJy beam$^{-1}$. The stochastic millimeter enhancements on time frames of seconds--minutes--hours with negative spectral indices are most plausibly explained by synchrotron or gyro-synchrotron radiation from stellar activity. HD 283572's 1.3 mm light-curve has similarities with variable binaries, suggesting HD 283572's activity may have been triggered by interactions with an as-yet undetected companion. We additionally identify variability of HD 283572 at 10 cm, from VLASS data. This study highlights the challenges of interpreting faint mm emission from evolved YSOs that may host tenuous disks, and suggests that a more detailed temporal analysis of spatially unresolved data is generally warranted. The variability of class III stars may open up new ground for understanding the physics of flares in the context of terrestrial planet formation. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL. 12 pages, inc. 5 figures, 2 tables

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

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

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

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ

High-contrast JWST-MIRI spectroscopy of planet-forming disks for the JDISC Survey

Authors: Klaus M. Pontoppidan, Colette Salyk, Andrea Banzatti, Ke Zhang, Ilaria Pascucci, Karin I. Oberg, Feng Long, Carlos Munoz-Romero, John Carr, Joan Najita, Geoffrey A. Blake, Nicole Arulanantham, Sean Andrews, Nicholas P. Ballering, Edwin Bergin, Jenny Calahan, Douglas Cobb, Maria Jose Colmenares, Annie Dickson-Vandervelde, Anna Dignan, Joel Green, Phoebe Heretz, Greg Herczeg, Anusha Kalyaan, Sebastian Krijt , et al. (4 additional authors not shown)

Abstract: The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4… ▽ More The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4.9-28 micron spectral range. This is achieved using observations of Themis-family asteroids as precise empirical reference sources. High spectral contrast enables precise retrievals of physical parameters, searches for rare molecular species and isotopologues, and constraints on the inventories of carbon- and nitrogen-bearing species. JDISCS also offers significant improvements to the MRS wavelength and resolving power calibration. We describe the JDISCS calibrated data and demonstrate its quality using observations of the disk around the solar-mass young star FZ Tau. The FZ Tau MIRI spectrum is dominated by strong emission from warm water vapor. We show that the water and CO line emission originates from the disk surface and traces a range of gas temperatures of ~500-1500 K. We retrieve parameters for the observed CO and H2O lines, and show that they are consistent with a radial distribution represented by two temperature components. A high water abundance of n(H2O)~10^-4 fills the disk surface at least out to the 350 K isotherm at 1.5 au. We search the FZ Tau environs for extended emission detecting a large (radius of ~300 au) ring of emission from H2 gas surrounding FZ Tau, and discuss its origin. △ Less

Submitted 16 January, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the Astrophysical Journal

Millimeter emission in photoevaporating disks is determined by early substructures

Authors: Matías Gárate, Til Birnstiel, Paola Pinilla, Sean M. Andrews, Raphael Franz, Sebastian Markus Stammler, Giovanni Picogna, Barbara Ercolano, Anna Miotello, Nicolás T. Kurtovic

Abstract: [abridged]Photoevaporation and dust-trapping are individually considered to be important mechanisms in the evolution and morphology of protoplanetary disks. We studied how the presence of early substructures affects the evolution of the dust distribution and flux in the millimeter continuum of disks that are undergoing photoevaporative dispersal. We also tested if the predicted properties resemble… ▽ More [abridged]Photoevaporation and dust-trapping are individually considered to be important mechanisms in the evolution and morphology of protoplanetary disks. We studied how the presence of early substructures affects the evolution of the dust distribution and flux in the millimeter continuum of disks that are undergoing photoevaporative dispersal. We also tested if the predicted properties resemble those observed in the population of transition disks. We used the numerical code Dustpy to simulate disk evolution considering gas accretion, dust growth, dust-trapping at substructures, and mass loss due to X-ray and EUV (XEUV) photoevaporation and dust entrainment. Then, we compared how the dust mass and millimeter flux evolve for different disk models. We find that, during photoevaporative dispersal, disks with primordial substructures retain more dust and are brighter in the millimeter continuum than disks without early substructures, regardless of the photoevaporative cavity size. Once the photoevaporative cavity opens, the estimated fluxes for the disk models that are initially structured are comparable to those found in the bright transition disk population ($F_\textrm{mm} > 30\, \textrm{mJy}$), while the disk models that are initially smooth have fluxes comparable to the transition disks from the faint population ($F_\textrm{mm} < 30\, \textrm{mJy}$), suggesting a link between each model and population. Our models indicate that the efficiency of the dust trapping determines the millimeter flux of the disk, while the gas loss due to photoevaporation controls the formation and expansion of a cavity, decoupling the mechanisms responsible for each feature. In consequence, even a planet with a mass comparable to Saturn could trap enough dust to reproduce the millimeter emission of a bright transition disk, while its cavity size is independently driven by photoevaporative dispersal. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 679, A15 (2023)

TU Tau B: The Peculiar 'Eclipse' of a possible proto-Barium Giant

Authors: Richard O. Gray, Christopher J. Corbally, Michael M. Briley, Adam McKay, Forrest Sims, David Boyd, Christophe Boussin, Courtney E. McGahee, Robert Buchheim, Gary Walker, David Iadevaia, David Cejudo Fernandez, Damien Lemay, Jack Martin, Jim Grubb, Albert Stiewing, Joseph Daglen, Keith Shank, Sydney Andrews, Nick Barnhardt, Rebekah Clark, Hunter Corman, Sabina Gomes, Agastya Jonnalagadda, Theo McDaries , et al. (4 additional authors not shown)

Abstract: TU Tau (= HD 38218 = HIP 27135) is a binary system consisting of a C-N carbon star primary and an A-type secondary. We report on new photometry and spectroscopy which tracked the recent disappearance of the A-star secondary. The dimming of the A-star was gradual and irregular, with one or more brief brightenings, implying the presence of nonhomogeneities in the carbon star outflow. We also present… ▽ More TU Tau (= HD 38218 = HIP 27135) is a binary system consisting of a C-N carbon star primary and an A-type secondary. We report on new photometry and spectroscopy which tracked the recent disappearance of the A-star secondary. The dimming of the A-star was gradual and irregular, with one or more brief brightenings, implying the presence of nonhomogeneities in the carbon star outflow. We also present evidence that the A-star is actively accreting s-process enriched material from the carbon star and suggest that it will therefore eventually evolve into a Barium giant. This is an important system as well because the A-type star can serve as a probe of the outer atmosphere of the carbon star. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 9 pages, 9 figures, 4 tables, a number of amateur observatories made significant contributions to this research. Paper accepted for publication in The Astronomical Journal

Journal ref: Astronomical Journal 166, 161-170 (2023)

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

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

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

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

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

A Large Double-ring Disk around the Taurus M Dwarf J04124068+2438157

Authors: Feng Long, Bin B. Ren, Nicole L. Wallack, Daniel Harsono, Gregory J. Herczeg, Paola Pinilla, Dimitri Mawet, Michael C. Liu, Sean M. Andrews, Xue-Ning Bai, Sylvie Cabrit, Lucas A. Cieza, Doug Johnstone, Jarron M. Leisenring, Giuseppe Lodato, Yao Liu, Carlo F. Manara, Gijs D. Mulders, Enrico Ragusa, Steph Sallum, Yangfan Shi, Marco Tazzari, Taichi Uyama, Kevin Wagner, David J. Wilner , et al. (1 additional authors not shown)

Abstract: Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcse… ▽ More Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $σ$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3σ$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5σ$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST. △ Less

Submitted 25 March, 2023; originally announced March 2023.

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

An SMA Survey of Chemistry in Disks around Herbig AeBe Stars

Authors: Jamila Pegues, Karin I. Öberg, Chunhua Qi, Sean M. Andrews, Jane Huang, Charles J. Law, Romane Le Gal, Luca Matrà, David J. Wilner

Abstract: Protoplanetary disks around Herbig AeBe stars are exciting targets for studying the chemical environments where giant planets form. Save for a few disks, however, much of Herbig AeBe disk chemistry is an open frontier. We present a Submillimeter Array (SMA) $\sim$213-268 GHz pilot survey of mm continuum, CO isotopologues, and other small molecules in disks around five Herbig AeBe stars (HD 34282,… ▽ More Protoplanetary disks around Herbig AeBe stars are exciting targets for studying the chemical environments where giant planets form. Save for a few disks, however, much of Herbig AeBe disk chemistry is an open frontier. We present a Submillimeter Array (SMA) $\sim$213-268 GHz pilot survey of mm continuum, CO isotopologues, and other small molecules in disks around five Herbig AeBe stars (HD 34282, HD 36112, HD 38120, HD 142666, and HD 144432). We detect or tentatively detect $^{12}$CO 2--1 and $^{13}$CO 2--1 from four disks; C$^{18}$O 2--1 and HCO$^+$ 3--2 from three disks; HCN 3--2, CS 5--4, and DCO$^+$ 3--2 from two disks; and C$_2$H 3--2 and DCN 3--2 from one disk each. H$_2$CO 3--2 is undetected at the sensitivity of our observations. The mm continuum images of HD 34282 suggest a faint, unresolved source $\sim$5\farcs0 away, which could arise from a distant orbital companion or an extended spiral arm. We fold our sample into a compilation of T Tauri and Herbig AeBe/F disks from the literature. Altogether, most line fluxes generally increase with mm continuum flux. Line flux ratios between CO 2--1 isotopologues are nearest to unity for the Herbig AeBe/F disks. This may indicate emitting layers with relatively similar, warmer temperatures and more abundant CO relative to disk dust mass. Lower HCO$^+$ 3--2 flux ratios may reflect less ionization in Herbig AeBe/F disks. Smaller detection rates and flux ratios for DCO$^+$ 3--2, DCN 3--2, and H$_2$CO 3--2 suggest smaller regimes of cold chemistry around the luminous Herbig AeBe/F stars. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 43 pages (23 pages in the main document, 20 pages in the appendix), 16 figures, 16 tables. Accepted in ApJ (February 2023)

The earliest DT nuclear fusion discoveries

Authors: M. B. Chadwick, G. M. Hale, M. W. Paris, J. P. Lestone, C. Bates, J. B. Wilhelmy, S. A. Andrews, W. Tornow, S. W. Finch

Abstract: We describe the earliest measurements of the DT fusion cross section commissioned by the Manhattan Project, first at Purdue University in 1943 and then at Los Alamos 1945-6 and later, in 1951-2. The Los Alamos measurements led to the realization that a 3/2$^+$ resonance in the DT system enhances the fusion cross section by a factor of one hundred at energies relevant to applications. This was a tr… ▽ More We describe the earliest measurements of the DT fusion cross section commissioned by the Manhattan Project, first at Purdue University in 1943 and then at Los Alamos 1945-6 and later, in 1951-2. The Los Alamos measurements led to the realization that a 3/2$^+$ resonance in the DT system enhances the fusion cross section by a factor of one hundred at energies relevant to applications. This was a transformational discovery, making the quest for terrestrial fusion energy possible. The earliest measurements were reasonably accurate given the technology of the time and the scarcity of tritium, and were quickly improved to provide cross section data accurate to just a few percent. We provide a previously-unappreciated insight: that DT fusion was first reported in Ruhlig's 1938 University of Michigan experiment and likely influenced Konopinski in 1942 to suggest its usefulness for thermonuclear technologies. We report on preliminary work to repeat the 1938 measurement, and our simulations of that experiment. We also present some work by Fermi, from his 1945 Los Alamos lectures, showing that he used the S-factor concept about a decade before it was introduced by nuclear astrophysicists. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 10 pages, 8 figures

Report number: LA-UR-22-32324

Kinematic and thermal signatures of the directly imaged protoplanet candidate around Elias 2-24

Authors: Christophe Pinte, Iain Hammond, Daniel J. Price, Valentin Christiaens, Sean M. Andrews, Gaël Chauvin, Laura M. Pérez, Sebastián Jorquera, Himanshi Garg, Brodie J. Norfolk, Josh Calcino, Mickaël Bonnefoy

Abstract: We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require loc… ▽ More We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require local heating. While the foreground extinction hides the velocity channels associated with the planet, preventing a planet mass estimate, the level of gas heating implied by the CO emission indicates the presence of a warm, embedded giant planet. Comparison with models show this could either be a $\gtrsim 5$M$_\mathrm{Jup}$, or a lower mass ( $\gtrsim 2$M$_\mathrm{Jup}$) but accreting proto-planet. △ Less

Submitted 20 January, 2023; originally announced January 2023.

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

Molecular Mapping of DR Tau's Protoplanetary Disk, Envelope, Outflow, and Large-Scale Spiral Arm

Authors: Jane Huang, Edwin A. Bergin, Jaehan Bae, Myriam Benisty, Sean M. Andrews

Abstract: DR Tau has been noted for its unusually high variability in comparison with other T Tauri stars. Although it is one of the most extensively studied pre-main sequence stars, observations with millimeter interferometry have so far been relatively limited. We present NOEMA images of $^{12}$CO, $^{13}$CO, C$^{18}$O, SO, DCO$^+$, and H$_2$CO toward DR Tau at a resolution of $\sim0.5''$ ($\sim100$ au).… ▽ More DR Tau has been noted for its unusually high variability in comparison with other T Tauri stars. Although it is one of the most extensively studied pre-main sequence stars, observations with millimeter interferometry have so far been relatively limited. We present NOEMA images of $^{12}$CO, $^{13}$CO, C$^{18}$O, SO, DCO$^+$, and H$_2$CO toward DR Tau at a resolution of $\sim0.5''$ ($\sim100$ au). In addition to the protoplanetary disk, CO emission reveals an envelope, a faint asymmetric outflow, and a spiral arm with a clump. The $\sim1200$ au extent of the CO arm far exceeds that of the spiral arms previously detected in scattered light, which underlines the necessity of sensitive molecular imaging for contextualizing the disk environment. The kinematics and compact emission distribution of C$^{18}$O, SO, DCO$^+$, and H$_2$CO indicate that they originate primarily from within the Keplerian circumstellar disk. The SO emission, though, also exhibits an asymmetry that may be due to interaction with infalling material or unresolved substructure. The complex environment of DR Tau is reminiscent of those of outbursting FUor sources and some EXor sources, suggesting that DR Tau's extreme stellar activity could likewise be linked to disk instabilities promoted by large-scale infall. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 29 pages, accepted by ApJ

Mapping Protoplanetary Disk Vertical Structure with CO Isotopologue Line Emission

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

Abstract: High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we… ▽ More High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we can map the 2D (r, z) disk structure from the disk upper atmosphere, as traced by CO, to near the midplane, as probed by less abundant isotopologues. Here, we present high angular resolution (${\lesssim}$0."1 to ${\approx}$0."2; ${\approx}$15-30 au) observations of CO, $^{13}$CO, and C$^{18}$O in either or both J=2-1 and J=3-2 lines in the transition disks around DM Tau, Sz 91, LkCa 15, and HD 34282. We derived line emission surfaces in CO for all disks and in $^{13}$CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C$^{18}$O in any disk, which is instead consistent with C$^{18}$O emission originating from the midplane. Both the J=2-1 and J=3-2 lines show similar heights. Using the derived emission surfaces, we computed radial and vertical gas temperature distributions for each disk, including empirical temperature models for the DM Tau and LkCa 15 disks. After combining our sample with literature sources, we find that $^{13}$CO line emitting heights are also tentatively linked with source characteristics, e.g., stellar host mass, gas temperature, disk size, and show steeper trends than seen in CO emission surfaces. △ Less

Submitted 16 December, 2022; originally announced December 2022.

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

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

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

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

Submitted 11 December, 2022; originally announced December 2022.

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

Millimeter Dust Emission and Planetary Dynamics in the HD 106906 System

Authors: Anna Fehr, A. Meredith Hughes, Rebekah I. Dawson, Rachel E. Marino, Matan Ackelsberg, Jamar Kittling, Kevin M. Flaherty, Erika Nesvold, John Carpenter, Sean M. Andrews, Brenda Matthews, Katie Crotts, Paul Kalas

Abstract: Debris disks are dusty, optically thin structures around main sequence stars. HD 106906AB is a short-period stellar binary, host to a wide separation planet, HD 106906b, and a debris disk. Only a few known systems include a debris disk and a directly imaged planet, and HD 106906 is the only one in which the planet is exterior to the disk. The debris disk is edge-on and highly asymmetric in scatter… ▽ More Debris disks are dusty, optically thin structures around main sequence stars. HD 106906AB is a short-period stellar binary, host to a wide separation planet, HD 106906b, and a debris disk. Only a few known systems include a debris disk and a directly imaged planet, and HD 106906 is the only one in which the planet is exterior to the disk. The debris disk is edge-on and highly asymmetric in scattered light. Here we resolve the disk structure at a resolution of 0.38" (39 au) with the Atacama Large Millimeter/submillimeter Array (ALMA) at a wavelength of 1.3 mm. We model the disk with both a narrow and broad ring of material, and find that a radially broad, axisymmetric disk between radii of $\sim$50$-$100 au is able to capture the structure of the observations without evidence of any asymmetry or eccentricity, other than a tentative stellocentric offset. We place stringent upper limits on both the gas and dust content of a putative circumplanetary disk. We interpret the ALMA data in concert with scattered light observations of the inner ring and astrometric constraints on the planet's orbit, and find that the observations are consistent with a large-separation, low-eccentricity orbit for the planet. A dynamical analysis indicates that the central binary can efficiently stabilize planetesimal orbits interior to $\sim$100 au, which relaxes the constraints on eccentricity and semimajor axis somewhat. The observational constraints are consistent with in situ formation via gravitational instability, but cannot rule out a scattering event as the origin for HD 106906b's current orbit. △ Less

Submitted 23 September, 2022; originally announced September 2022.

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

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

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

Submitted 12 September, 2022; originally announced September 2022.

Comments: ApJL in press

Mapping the Complex Kinematic Substructure in the TW Hya Disk

Authors: Richard Teague, Jaehan Bae, Sean M. Andrews, Myriam Benisty, Edwin A. Bergin, Stefano Facchini, Jane Huang, Cristiano Longarini, David Wilner

Abstract: We present ALMA observations of CO $J = 2-1$ and CS $J = 5-4$ emission from the disk around TW~Hydrae. Both molecules trace a predominantly Keplerian velocity structure, although a slowing of the rotation velocity is detected at the outer edge of the disk beyond ${\approx}~140$~au in CO emission. This was attributed to the enhanced pressure support from the gas density taper near the outer edge of… ▽ More We present ALMA observations of CO $J = 2-1$ and CS $J = 5-4$ emission from the disk around TW~Hydrae. Both molecules trace a predominantly Keplerian velocity structure, although a slowing of the rotation velocity is detected at the outer edge of the disk beyond ${\approx}~140$~au in CO emission. This was attributed to the enhanced pressure support from the gas density taper near the outer edge of the disk. Subtraction of an azimuthally symmetric background velocity structure reveals localized deviations in the gas kinematics traced by each of the molecules. Both CO and CS exhibit a `Doppler flip' feature, centered nearly along the minor axis of the disk (${\rm PA} \sim 60\degr$) at a radius of $1\farcs35$, coinciding with the large gap observed in scattered light and mm~continuum. In addition, the CO emission, both through changes in intensity and its kinematics, traces a tightly wound spiral, previously seen with higher frequency CO $J = 3-2$ observations (Teague et al., 2019). Through comparison with linear models of the spiral wakes generated by embedded planets, we interpret these features in the context of interactions with a Saturn-mass planet within the gap at a position angle of ${\rm PA} = 60\degr$, consistent with the theoretical predictions of (Mentiplay et al. 2019). The lack of a corresponding spiral in the CS emission is attributed to the strong vertical dependence on the buoyancy spirals which are believed to only grow in the atmospheric of the disk, rather than those traced by CS emission. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted to ApJ

An ALMA 1.3 millimeter Search for Debris Disks around Solar-type Stars in the Pleiades

Authors: Devin Sullivan, David Wilner, Luca Matra, Mark C. Wyatt, Sean M. Andrews, Meredith A. MacGregor, Brenda Matthews

Abstract: Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust… ▽ More Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust at lower temperatures, corresponding to debris at radii comparable to the Kuiper Belt and beyond. The observations obtain a beam size of 1.5 arcsec (200 au) and a median rms noise of 54 mircoJy/beam, which corresponds to a fractional luminosity $L_{dust}/L_{star} \sim 10^{-4}$ for 40 K dust for a typical star in the sample. The ALMA images show no significant detections of the targeted stars. We interpret these limits in the context of a steady-state collisional cascade model for debris disk evolution that provides a good description of observations of the field population near the Sun but is not well-calibrated on younger populations.The ALMA non-detections of the Pleiades systems are compatible with the disk flux predictions of this model. We find no high fractional luminosity outliers from these ALMA data that could be associated with enhanced collisions resulting from activity not accounted for by steady-state evolution. However, we note that two systems (HII 1132 and HD 22680) show 24 micron excess much higher than the predictions of this model, perhaps due to unusually high dust production from dynamical events involving planets. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, accepted for publication in AJ

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/

Multiwavelength Vertical Structure in the AU Mic Debris Disk: Characterizing the Collisional Cascade

Authors: David Vizgan, A. Meredith Hughes, Evan S. Carter, Kevin M. Flaherty, Margaret Pan, Eugene Chiang, Hilke Schlichting, David J. Wilner, Sean M. Andrews, John M. Carpenter, Atilla Moór, Meredith A. MacGregor

Abstract: Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested… ▽ More Debris disks are scaled-up analogs of the Kuiper Belt in which dust is generated by collisions between planetesimals. In the "collisional cascade" model of debris disks, dust lost to radiation pressure and winds is constantly replenished by grinding collisions between planetesimals. The model assumes that collisions are destructive and involve large velocities; this assumption has not been tested beyond our Solar System. We present 0"25 ($\approx$2.4 au) resolution observations of the $λ$ = 450 $μ$m dust continuum emission from the debris disk around the nearby M dwarf AU Microscopii with the Atacama Large Millimeter/submillimeter Array. We use parametric models to describe the disk structure, and an MCMC algorithm to explore the posterior distributions of the model parameters; we fit the structure of the disk to both our data and archival $λ= 1.3$ mm data (Daley et al. 2019), from which we obtain two aspect ratio measurements at 1.3 mm ($h_{1300}$ = 0.025$^{+0.008}_{-0.002}$) and at 450 $μ$m ($h_{450}$ = 0.019$^{+0.006}_{-0.001}$), as well as the grain size distribution index $q =$ 3.03 $\pm$ 0.02. Contextualizing our aspect ratio measurements within the modeling framework laid out in Pan & Schlichting (2012), we derive a power law index of velocity dispersion as a function of grain size $p = 0.28 \pm 0.06$ for the AU Mic debris disk. This result implies that smaller bodies are more easily disrupted than larger bodies by collisions, which is inconsistent with the strength regime usually assumed for such small bodies. Possible explanations for this discrepancy are discussed. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 17 pages, 6 figures, 3 tables. Accepted for publication in ApJ

Distributions of gas and small and large grains in the LkH$α\,330$ disk trace a young planetary system

Authors: P. Pinilla, M. Benisty, N. T. Kurtovic, J. Bae, R. Dong, Z. Zhu, S. Andrews, J. Carpenter, C. Ginski, J. Huang, A. Isella, L. Pérez, L. Ricci, G. Rosotti, M. Villenave, D. Wilner

Abstract: [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''… ▽ More [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') mainly show an asymmetric ring located at 110au from the star. In addition to this asymmetry, there are two faint symmetric rings at 60au and 200au. The $^{12}$CO, $^{13}$CO, and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface-density variations, while the C$^{18}$O peaks slightly further away at around 60au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10$M_{\rm{Jup}}$ mass planet at 60au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$α\,330$ is an interesting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction. △ Less

Submitted 25 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A, minor changes after language edition

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

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

Authors: Evan A. Rich, John D. Monnier, Alicia Aarnio, Anna S. E. Laws, Benjamin R. Setterholm, David J. Wilner, Nuria Calvet, Tim Harries, Chris Miller, Claire L. Davies, Fred C. Adams, Sean M. Andrews, Jaehan Bae, Catherine Espaillat, Alexandra Z. Greenbaum, Sasha Hinkley, Stefan Kraus, Lee Hartmann, Andrea Isella, Melissa McClure, Rebecca Oppenheimer, Laura M. Pérez, Zhaohuan Zhu

Abstract: We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of… ▽ More We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 51 pages, 31 figures, 7 tables, accepted to AJ

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

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

Abstract: High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure o… ▽ More High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$. △ Less

Submitted 3 May, 2022; originally announced May 2022.

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

Gas and Dust Shadows in the TW Hydrae Disk

Authors: Richard Teague, Jaehan Bae, Myriam Benisty, Sean M. Andrews, Stefano Facchini, Jane Huang, David Wilner

Abstract: We present new observations of CO J=2-1 emission from the protoplanetary disk around TW Hya. Emission is detected out to 240 au (4") and found to exhibit azimuthal variations up to 20% beyond 180 au (3"), with the west side of the disk brighter than the east. This asymmetry is interpreted as tracing the shadow previously seen in scattered light. A reanalysis of the multi-epoch observations of the… ▽ More We present new observations of CO J=2-1 emission from the protoplanetary disk around TW Hya. Emission is detected out to 240 au (4") and found to exhibit azimuthal variations up to 20% beyond 180 au (3"), with the west side of the disk brighter than the east. This asymmetry is interpreted as tracing the shadow previously seen in scattered light. A reanalysis of the multi-epoch observations of the dust shadow in scattered light from Debes et al. (2017) suggests that an oscillatory motion would provide a better model of the temporal evolution of the dust shadow rather than orbital motion. Both models predict an angular offset between the dust shadow and the gas shadow of up to ~100 deg. We attribute this offset to the finite rate at which dust grains and gas molecules can exchange heat, dominated by the collisional rate between gas molecules and dust grains, $t_{\rm coll}$. The angular offsets derived are equivalent to collisional timescales that range from the near instantaneous up to $t_{\rm coll}$ ~ 10 years, depending on whether a straight or a curved dust shadow, as suggested by HST observations reported by Debes et al. (2017), is adopted. The inferred range of $t_{\rm coll}$ are consistent with those predictions based on representative gas densities, temperatures, gas-to-dust ratios and grain sizes. These results represent the first time empirical constraints can be placed on $t_{\rm coll}$. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: Accepted by ApJ

Directly Detecting the Envelopes of Low-mass Planets Embedded In Protoplanetary Discs and The Case For TW Hydrae

Authors: Zhaohuan Zhu, Avery Bailey, Enrique Macías, Takayuki Muto, Sean M. Andrews

Abstract: Despite many methods developed to find young massive planets in protoplanetary discs, it is challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1-D models and 3-D simulations to calculat… ▽ More Despite many methods developed to find young massive planets in protoplanetary discs, it is challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1-D models and 3-D simulations to calculate the envelopes around low-mass cores (several to tens of $M_{\oplus}$) with different luminosities, and derive their thermal fluxes at radio wavelengths. We find that, when the background disc is optically thin at radio wavelengths, radio observations can see through the disc and probe the denser envelope within the planet's Hill sphere. When the optically thin disc is observed with the resolution reaching one disc scale height, the radio thermal flux from the planetary envelope around a 10 M$_{\oplus}$ core is more than 10 % higher than the flux from the background disc. The emitting region can be extended and elongated. Finally, our model suggests that the au-scale clump at 52 au in the TW Hydrae disc revealed by ALMA is consistent with the envelope of an embedded 10-20 $M_{\oplus}$ planet, which can explain the detected flux, the spectral index dip, and the tentative spirals. The observation is also consistent with the planet undergoing pebble accretion. Future ALMA and ngVLA observations may directly reveal more such low-mass planets, enabling us to study core growth and even reconstruct the planet formation history using the embedded "protoplanet" population. △ Less

Submitted 16 September, 2022; v1 submitted 9 April, 2022; originally announced April 2022.

Comments: 20 pages, 14 figures, MNRAS in press

Gas Disk Sizes from CO Line Observations: A Test of Angular Momentum Evolution

Authors: Feng Long, Sean M. Andrews, Giovanni Rosotti, Daniel Harsono, Paola Pinilla, David J. Wilner, Karin I Öberg, Richard Teague, Leon Trapman, Benoît Tabone

Abstract: The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emis… ▽ More The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emission span from 50 to 1000\,au. This range could be explained by viscous evolution models with different $α$ values (mostly of $10^{-4}-10^{-3}$) and/or a spread of initial conditions. The CO sizes for most disks are also consistent with MHD wind models that directly remove disk angular momentum, but very large initial disk sizes would be required to account for the very extended CO disks in the sample. As no CO size evolution is observed across stellar ages of 0.5--20\,Myr in this sample, determining the dominant mechanism of disk evolution will require a more complete sample for both younger and more evolved systems. We find that the CO emission is universally more extended than the continuum emission by an average factor of $2.9\pm1.2$. The ratio of the CO to continuum sizes does not show any trend with stellar mass, mm continuum luminosity, or the properties of substructures. The GO Tau disk has the most extended CO emission in this sample, with an extreme CO to continuum size ratio of 7.6. Seven additional disks in the sample show high size ratios ($\gtrsim4$) that we interpret as clear signs of substantial radial drift. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

A large population study of protoplanetary disks: Explaining the millimeter size-luminosity relation with or without sub-structure

Authors: Apostolos Zormpas, Tilman Birnstiel, Giovanni P. Rosotti, Sean M. Andrews

Abstract: Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolut… ▽ More Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolution in a simple viscous accretion disk and radiative transfer calculations. We use a large grid of models ($10^{5}$ simulations) with and without planetary gaps, varying the initial conditions of the key parameters. We calculate the disk continuum emission and the effective radius for all models as a function of time. By selecting those simulations that continuously follow the SLR, we can derive constraints on the input parameters of the models. We confirm previous results that models of smooth disks in the radial drift regime are compatible with the observed SLR ($L_{mm}\propto r_{eff}^{2}$) but only smooth disks cannot be the reality. We show that the SLR is more widely populated if planets are present. However they tend to follow a different relation than smooth disks, potentially implying that a mixture of smooth and sub-structured disks are present in the observed sample. We derive a SLR ($L_{mm}\propto r_{eff}^{5/4}$) for disks with strong sub-structure. To be compatible with the SLR, models need to have an initially high disk mass ($\geq 2.5 \cdot 10^{-2}M_{\star}$) and low turbulence-parameter $α$ values ($\leq 10^{-3}$). Furthermore, we find that the grain composition and porosity drastically affects the evolution of disks on the size-luminosity diagram where relatively compact grains that include amorphous carbon are favoured. Moreover, a uniformly optically thick disk with high albedo ($0.9$) that follows the SLR cannot be formed from an evolutionary procedure. △ Less

Submitted 23 May, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: Replaced with the published version in A&A. 19+11 pages, 11+10 figures

Journal ref: A&A 661, A66 (2022)

Scattering and sublimation: a multi-scale view of $μ$m-sized dust in the inclined disc of HD 145718

Authors: Claire L. Davies, Evan A. Rich, Tim J. Harries, John D. Monnier, Anna S. E. Laws, Sean M. Andrews, Jaehan Bae, David J. Wilner, Narsireddy Anugu, Jacob Ennis, Tyler Gardner, Stefan Kraus, Aaron Labdon, Jean-Baptiste le Bouquin, Cyprien Lanthermann, Gail H. Schaefer, Benjamin R. Setterholm, Theo ten Brummelaar

Abstract: We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angl… ▽ More We present multi-instrument observations of the disc around the Herbig~Ae star, HD~145718, employing geometric and Monte Carlo radiative transfer models to explore the disc orientation, the vertical and radial extent of the near infrared (NIR) scattering surface, and the properties of the dust in the disc surface and sublimation rim. The disc appears inclined at $67-71^{\circ}$, with position angle, PA\,$=-1.0-0.6^{\circ}$, consistent with previous estimates. The NIR scattering surface extends out to $\sim75\,$au and we infer an aspect ratio, $h_{\rm{scat}}(r)/r\sim0.24$ in $J$-band; $\sim0.22$ in $H$-band. Our GPI images and VLTI+CHARA NIR interferometry suggest that the disc surface layers are populated by grains $\gtrsim λ/2π$ in size, indicating these grains are aerodynamically supported against settling and/or the density of smaller grains is relatively low. We demonstrate that our geometric analysis provides a reasonable assessment of the height of the NIR scattering surface at the outer edge of the disc and, if the inclination can be independently constrained, has the potential to probe the flaring exponent of the scattering surface in similarly inclined ($i\gtrsim70^{\circ}$) discs. In re-evaluating HD~145718's stellar properties, we found that the object's dimming events - previously characterised as UX~Or and dipper variability - are consistent with dust occultation by grains larger, on average, than found in the ISM. This occulting dust likely originates close to the inferred dust sublimation radius at $0.17\,$au. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 18 pages, 16 figures; Accepted for publication in MNRAS

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

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

Abstract: UV photochemistry in the surface layers of protoplanetary disks dramatically alters their composition relative to previous stages of star formation. The abundance ratio CN/HCN has long been proposed to trace the UV field in various astrophysical objects, however to date the relationship between CN, HCN, and the UV field in disks remains ambiguous. As part of the ALMA Large Program MAPS (Molecules… ▽ More UV photochemistry in the surface layers of protoplanetary disks dramatically alters their composition relative to previous stages of star formation. The abundance ratio CN/HCN has long been proposed to trace the UV field in various astrophysical objects, however to date the relationship between CN, HCN, and the UV field in disks remains ambiguous. As part of the ALMA Large Program MAPS (Molecules with ALMA at Planet-forming Scales), we present observations of CN N=1-0 transitions at 0.3'' resolution towards five disk systems. All disks show bright CN emission within $\sim$50-150 au, along with a diffuse emission shelf extending up to 600 au. In all sources we find that the CN/HCN column density ratio increases with disk radius from about unity to 100, likely tracing increased UV penetration that enhances selective HCN photodissociation in the outer disk. Additionally, multiple millimeter dust gaps and rings coincide with peaks and troughs, respectively, in the CN/HCN ratio, implying that some millimeter substructures are accompanied by changes to the UV penetration in more elevated disk layers. That the CN/HCN ratio is generally high (>1) points to a robust photochemistry shaping disk chemical compositions, and also means that CN is the dominant carrier of the prebiotically interesting nitrile group at most disk radii. We also find that the local column densities of CN and HCN are positively correlated despite emitting from vertically stratified disk regions, indicating that different disk layers are chemically linked. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

Abstract: During the main phase of evolution of a protoplanetary disk, accretion regulates the inner-disk properties, such as the temperature and mass distribution, and in turn, the physical conditions associated with planet formation. The driving mechanism behind accretion remains uncertain; however, one promising mechanism is the removal of a fraction of angular momentum via a magnetohydrodynamic (MHD) di… ▽ More During the main phase of evolution of a protoplanetary disk, accretion regulates the inner-disk properties, such as the temperature and mass distribution, and in turn, the physical conditions associated with planet formation. The driving mechanism behind accretion remains uncertain; however, one promising mechanism is the removal of a fraction of angular momentum via a magnetohydrodynamic (MHD) disk wind launched from the inner tens of astronomical units of the disk. This paper utilizes CO isotopologue emission to study the unique molecular outflow originating from the HD 163296 protoplanetary disk obtained with the Atacama Large Millimeter/submillimeter Array. HD~163296 is one of the most well-studied Class II disks and is proposed to host multiple gas-giant planets. We robustly detect the large-scale rotating outflow in the 12CO J=2-1 and the 13CO J=2-1 and J=1-0 transitions. We constrain the kinematics, the excitation temperature of the molecular gas, and the mass-loss rate. The high ratio of the rates of ejection to accretion (5 - 50), together with the rotation signatures of the flow, provides solid evidence for an MHD disk wind. We find that the angular momentum removal by the wind is sufficient to drive accretion through the inner region of the disk; therefore, accretion driven by turbulent viscosity is not required to explain HD~163296's accretion. The low temperature of the molecular wind and its overall kinematics suggest that the MHD disk wind could be perturbed and shocked by the previously observed high-velocity atomic jet. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

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

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

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

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

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

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

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

Abstract: Constraining dust properties of planet-forming disks via high angular resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distributi… ▽ More Constraining dust properties of planet-forming disks via high angular resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distribution of solid mass and maximum grain size in each disk, including or excluding dust scattering. In the nonscattering model, the dust surface density and maximum grain size profiles decrease from the inner disks to the outer disks, with local maxima at the bright ring locations, as expected from dust trapping models. The inferred maximum grain sizes from the inner to outer disks decrease from ~1 cm to 1 mm. For IM Lup, HD 163296, and MWC 480 in the scattering model, two solutions are compatible with their observed inner disk emission: one solution corresponding to a maximum grain size of a few millimeters (similar to the nonscattering model), and the other corresponding to a few hundred micrometer sizes. Based on the estimated Toomre parameter, only IM Lup -- which shows a prominent spiral morphology in millimeter dust -- is found to be gravitationally unstable. The estimated maximum Stokes number in all the disks lies between 0.01 and 0.3, and the estimated turbulence parameters in the rings of AS 209 and HD 163296 are close to the threshold where dust growth is limited by turbulent fragmentation. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

Abstract: We observed HCO$^+$ $J=1-0$ and H$^{13}$CO$^+$ $J=1-0$ emission towards the five protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 as part of the MAPS project. HCO$^+$ is detected and mapped at 0.3\arcsec\,resolution in all five disks, while H$^{13}$CO$^+$ is detected (SNR$>6 σ$) towards GM Aur and HD 163296 and tentatively detected (SNR$>3 σ$) towards the other disks by a… ▽ More We observed HCO$^+$ $J=1-0$ and H$^{13}$CO$^+$ $J=1-0$ emission towards the five protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 as part of the MAPS project. HCO$^+$ is detected and mapped at 0.3\arcsec\,resolution in all five disks, while H$^{13}$CO$^+$ is detected (SNR$>6 σ$) towards GM Aur and HD 163296 and tentatively detected (SNR$>3 σ$) towards the other disks by a matched filter analysis. Inside a radius of $R\sim 100$ au, the HCO$^+$ column density is flat or shows a central dip. At outer radii ($\gtrsim 100$ au), the HCO$^+$ column density decreases outwards, while the column density ratio of HCO$^+$/CO is mostly in the range of $\sim 10^{-5}-10^{-4}$. We derived the HCO$^+$ abundance in the warm CO-rich layer, where HCO$^+$ is expected to be the dominant molecular ion. At $R\gtrsim 100$ au, the HCO$^+$ abundance is $\sim 3 \times 10^{-11} - 3\times 10^{-10}$, which is consistent with a template disk model with X-ray ionization. At the smaller radii, the abundance decreases inwards, which indicates that the ionization degree is lower in denser gas, especially inside the CO snow line, where the CO-rich layer is in the midplane. Comparison of template disk models with the column densities of HCO$^+$, N$_2$H$^+$, and N$_2$D$^+$ indicates that the midplane ionization rate is $\gtrsim 10^{-18}$ s$^{-1}$ for the disks around IM Lup, AS 209, and HD 163296. We also find hints of an increased HCO$^+$ abundance around the location of dust continuum gaps in AS 209, HD 163296, and MWC 480. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

Comments: accepted to ApJS, 33 pages, 20 figures

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

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

Abstract: Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high angular resolution (10-50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find… ▽ More Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high angular resolution (10-50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find that the HCN emission arises in a temperate (20-30 K) layer in the disk, while C2H is present in relatively warmer (20-60 K) layers. In the case of HD 163296, we find a decrease in column density for HCN and C2H inside one of the dust gaps near 83 au, where a planet has been proposed to be located. We derive H2CO column density profiles assuming temperatures between 20 and 50 K, and find slightly higher column densities in the colder disks around T Tauri stars than around Herbig Ae stars. The H2CO column densities rise near the location of the CO snowline and/or millimeter dust edge, suggesting an efficient release of H2CO ices in the outer disk. Finally, we find that the inner 50 au of these disks are rich in organic species, with abundances relative to water that are similar to cometary values. Comets could therefore deliver water and key organics to future planets in these disks, similar to what might have happened here on Earth. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 September, 2021; originally announced September 2021.

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

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

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

Abstract: The precursors to larger, biologically-relevant molecules are detected throughout interstellar space, but determining the presence and properties of these molecules during planet formation requires observations of protoplanetary disks at high angular resolution and sensitivity. Here we present 0.3" observations of HC$_3$N, CH$_3$CN, and $c$-C$_3$H$_2$ in five protoplanetary disks observed as part… ▽ More The precursors to larger, biologically-relevant molecules are detected throughout interstellar space, but determining the presence and properties of these molecules during planet formation requires observations of protoplanetary disks at high angular resolution and sensitivity. Here we present 0.3" observations of HC$_3$N, CH$_3$CN, and $c$-C$_3$H$_2$ in five protoplanetary disks observed as part of the Molecules with ALMA at Planet-forming Scales (MAPS) Large Program. We robustly detect all molecules in four of the disks (GM Aur, AS 209, HD 163296 and MWC 480) with tentative detections of $c$-C$_3$H$_2$ and CH$_3$CN in IM Lup. We observe a range of morphologies -- central peaks, single or double rings -- with no clear correlation in morphology between molecule nor disk. Emission is generally compact and on scales comparable with the millimetre dust continuum. We perform both disk-integrated and radially-resolved rotational diagram analysis to derive column densities and rotational temperatures. The latter reveals 5-10 times more column density in the inner 50-100 au of the disks when compared with the disk-integrated analysis. We demonstrate that CH$_3$CN originates from lower relative heights in the disks when compared with HC$_3$N, in some cases directly tracing the disk midplane. Finally, we find good agreement between the ratio of small to large nitriles in the outer disks and comets. Our results indicate that the protoplanetary disks studied here are host to significant reservoirs of large organic molecules, and that this planet- and comet-building material can be chemically similar to that in our own Solar System. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series. △ Less

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

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

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

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

Abstract: Sulfur-bearing molecules play an important role in prebiotic chemistry and planet habitability. They are also proposed probes of chemical ages, elemental C/O ratio, and grain chemistry processing. Commonly detected in diverse astrophysical objects, including the Solar System, their distribution and chemistry remain, however, largely unknown in planet-forming disks. We present CS ($2-1$) observatio… ▽ More Sulfur-bearing molecules play an important role in prebiotic chemistry and planet habitability. They are also proposed probes of chemical ages, elemental C/O ratio, and grain chemistry processing. Commonly detected in diverse astrophysical objects, including the Solar System, their distribution and chemistry remain, however, largely unknown in planet-forming disks. We present CS ($2-1$) observations at $\sim0."3$ resolution performed within the ALMA-MAPS Large Program toward the five disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. CS is detected in all five disks, displaying a variety of radial intensity profiles and spatial distributions across the sample, including intriguing apparent azimuthal asymmetries. Transitions of C$_2$S and SO were also serendipitously covered but only upper limits are found. For MWC 480, we present complementary ALMA observations at $\sim0."5$, of CS, $^{13}$CS, C$^{34}$S, H$_2$CS, OCS, and SO$_2$. We find a column density ratio N(H$_{2}$CS)/N(CS)$\sim2/3$, suggesting that a substantial part of the sulfur reservoir in disks is in organic form (i.e., C$_x$H$_y$S$_z$). Using astrochemical disk modeling tuned to MWC 480, we demonstrate that $N$(CS)/$N$(SO) is a promising probe for the elemental C/O ratio. The comparison with the observations provides a super-solar C/O. We also find a depleted gas-phase S/H ratio, suggesting either that part of the sulfur reservoir is locked in solid phase or that it remains in an unidentified gas-phase reservoir. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

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

Abstract: Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission fro… ▽ More Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks - around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 - in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in 4 spectral setups, which together cover ~50 lines from over 20 different species. This paper introduces the ApJS MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical sub-structures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation. △ Less

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

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

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

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

Abstract: Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial resolution data from the Molecules with ALMA at Planet-forming Scales(MAPS) program. The observations of C$^{18}$O J=2-1 emission exhibit a… ▽ More Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial resolution data from the Molecules with ALMA at Planet-forming Scales(MAPS) program. The observations of C$^{18}$O J=2-1 emission exhibit a strong depression at 88 au overlapping with the positions of multiple gaps in millimeter dust continuum emission. We find that the observed CO column density is consistent with either gas surface-density perturbations or chemical processing, while C$_2$H column density traces changes in the C/O ratio rather than the H$_2$ gas surface density. However, the presence of a massive planet (> 0.2 M$_{Jup}$) would be required to account for this level of gas depression, which conflicts with constraints set by the dust emission and the pressure profile measured by gas kinematics. Based on our models, we infer that a local decrease of CO abundance is required to explain the observed structure in CO, dominating over a possible gap-carving planet present and its effect on the H$_2$ surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series. △ Less

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

Comments: 21 pages, 18 figures

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

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

Abstract: Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling f… ▽ More Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C$^{17}$O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H$_2$ abundance of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet-disk interactions, but some CO gaps are 5-10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

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

Abstract: Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk tempe… ▽ More Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk temperature structure. Using new Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program as well as archival ALMA observations, we construct a disk physical/chemical model of the protoplanetary disk GM Aur. Our model is in good agreement with the spatially resolved CO isotopolog emission from eleven rotational transitions with spatial resolution ranging from 0.15'' to 0.46'' (24-73 au at 159 pc) and the spatially unresolved HD J=1-0 detection from Herschel. Our best-fit model favors a cold protoplanetary disk with a total gas mass of approximately 0.2 solar masses, a factor of 10 reduction in CO gas inside roughly 100 au and a factor of 100 reduction outside of 100 au. Despite its large mass, the disk appears to be on the whole gravitationally stable based on the derived Toomre Q parameter. However, the region between 70 and 100 au, corresponding to one of the millimeter dust rings, is close to being unstable based on the calculated Toomre Q of <1.7. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 18 pages, 12 figures