Showing 1–50 of 115 results for author: Guzman, V

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

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

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

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

Submitted 1 May, 2024; originally announced May 2024.

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

Radial and vertical constraints on the icy origin of H$_{2}$CO in the HD 163296 Protoplanetary Disk

Authors: Claudio Hernández-Vera, Viviana V. Guzmán, Elizabeth Artur de la Villarmois, Karin I. Öberg, L. Ilsedore Cleeves, Michiel R. Hogerheijde, Chunhua Qi, John Carpenter, Edith C. Fayolle

Abstract: H$_2$CO is a small organic molecule widely detected in protoplanetary disks. As a precursor to grain-surface formation of CH$_3$OH, H$_2$CO is considered an important precursor of O-bearing organic molecules that are locked in ices. Still, since gas-phase reactions can also form H$_2$CO, there remains an open question on the channels by which organics form in disks, and how much the grain versus t… ▽ More H$_2$CO is a small organic molecule widely detected in protoplanetary disks. As a precursor to grain-surface formation of CH$_3$OH, H$_2$CO is considered an important precursor of O-bearing organic molecules that are locked in ices. Still, since gas-phase reactions can also form H$_2$CO, there remains an open question on the channels by which organics form in disks, and how much the grain versus the gas pathways impact the overall organic reservoir. We present spectrally and spatially resolved Atacama Large Millimeter/submillimeter Array observations of several ortho- and para-H$_2$CO transitions toward the bright protoplanetary disk around the Herbig Ae star HD 163296. We derive column density, excitation temperature, and ortho-to-para ratio (OPR) radial profiles for H$_2$CO, as well as disk-averaged values of $N_{\mathrm{T}}\sim4\times 10^{12}$ cm$^{-2}$, $T_{\mathrm{ex}}\sim20$ K, and $\mathrm{OPR}\sim2.7$, respectively. We empirically determine the vertical structure of the emission, finding vertical heights of $z/r\sim0.1$. From the profiles, we find a relatively constant $\mathrm{OPR}\sim2.7$ with radius, but still consistent with $3.0$ among the uncertainties, a secondary increase of $N_{\mathrm{T}}$ in the outer disk, and low $T_{\mathrm{ex}}$ values that decrease with disk radius. Our resulting radial, vertical, and OPR constraints suggest an increased UV penetration beyond the dust millimeter edge, consistent with an icy origin but also with cold gas-phase chemistry. This Herbig disk contrasts previous results for the T Tauri disk, TW Hya, which had a larger contribution from cold gas-phase chemistry. More observations of other sources are needed to disentangle the dominant formation pathway of H$_2$CO in protoplanetary disks. △ Less

Submitted 24 May, 2024; v1 submitted 9 April, 2024; originally announced April 2024.

Comments: 21 pages, 4 tables, 10 figures; V2: minor changes, matches the published version in ApJ

Journal ref: ApJ 967(1), 68 (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)

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

One loop reduced QED for massive fermions within an innovative formalism

Authors: Victor Miguel Banda Guzmán, Adnan Bashir, Luis Albino, Dania Rodríguez-Tzintzun

Abstract: We carry out a detailed study of the three-point fermion-photon interaction vertex at one loop order for massive fermions in reduced quantum electrodynamics. This calculation is carried out in arbitrary covariant gauges and space-time dimensions within a recently proposed innovative approach based upon an efficient combination of the first and second order formalisms of quantum electrodynamics. Th… ▽ More We carry out a detailed study of the three-point fermion-photon interaction vertex at one loop order for massive fermions in reduced quantum electrodynamics. This calculation is carried out in arbitrary covariant gauges and space-time dimensions within a recently proposed innovative approach based upon an efficient combination of the first and second order formalisms of quantum electrodynamics. This procedure provides a natural decomposition of the vertex into its components which are longitudinal and transverse to the photon momentum. It also separates the spin and scalar degrees of freedom of a fermion interacting electromagnetically, allowing us to readily establish the gauge-independence of the Pauli form factor and compute it in an expeditious manner. All incoming and outgoing momenta are taken off-shell at the outset. However, we present results for cases of particular kinematic interest whenever required. For the sake of completeness, we also provide expressions for the massive fermion self energy and photon vacuum polarization, verifying known expressions for massless reduced quantum electrodynamics and computing the renormalization constants ${\cal Z}_1$, ${\cal Z}_2$ and ${\cal Z}_3$. As we provide general expressions for the computed Green functions, we readily reproduce and confirm the results for standard quantum electrodynamics. Comparing the two cases, we infer that the Pauli form factor for reduced quantum electrodynamics is $8/3$ times that for the standard QED in four dimensions, implying a higher Landé $g$-factor. We expect our perturbative calculation of the fermion-photon vertex to serve as a guide for any non-perturbative construction of this Green function, invariably required in the Schwinger-Dyson equation studies of the subject. △ Less

Submitted 14 October, 2023; originally announced October 2023.

Comments: 14 pages, 5 figures

The Perseus ALMA Chemical Survey (PEACHES). III. Sulfur-bearing species tracing accretion and ejection processes in young protostars

Authors: E. Artur de la Villarmois, V. V. Guzman, Y. -L. Yang, Y. Zhang, N. Sakai

Abstract: (Abridged) Sulfur chemistry is poorly understood in the process of low-mass star and planet formation, where the main carriers of sulfur are still unknown. Despite the fact that simple S-bearing molecules are usually detected toward embedded sources, large surveys of S-bearing molecules with high angular resolution and sensitive observations are currently lacking. The goal of this work is to prese… ▽ More (Abridged) Sulfur chemistry is poorly understood in the process of low-mass star and planet formation, where the main carriers of sulfur are still unknown. Despite the fact that simple S-bearing molecules are usually detected toward embedded sources, large surveys of S-bearing molecules with high angular resolution and sensitive observations are currently lacking. The goal of this work is to present an unbiased survey of simple sulfur-bearing species in protostars and provide new statistics. In addition, we investigate the role of S-bearing molecules in accretion processes and the connection between (non-)detection of complex organic molecules (COMs) and S-related species. We present the observations of sulfur-bearing species that are part of the Perseus ALMA Chemical Survey (PEACHES). We analyzed a total of 50 Class 0/I sources with an average angular resolution of about 0.6" (~180 au) in ALMA band 6. We present detection rates for CS, SO, 34SO, and SO2. The SO/34SO ratio is lower than the canonical value of 22 and the lowest values are found for those sources rich in COMs. This ratio, therefore, seems to be a good tracer of the inner high-density envelope. The detection of multiple COMs seems to be related to the presence of collimated outflows and SO2 emission seems to trace the warm gas in those sources where CH3OH is also detected. The SO2 abundances toward the PEACHES sample are, on average, two orders of magnitude lower than values from the Ophiuchus star-forming region and comparable with sources in Taurus, suggesting that the sulfur depletion in the gas-phase could depend on the external UV radiation. Finally, the SO2 emission detected in different evolutionary stages seems to arise from different physical mechanisms: high column density of warm material in Class 0 sources, shocks in Class I/II, and exposure to UV radiation from the protostar in more evolved Class II disks. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 31 pages, 17 Figures (including appendix)

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

HCN emission from translucent gas and UV-illuminated cloud edges revealed by wide-field IRAM 30m maps of Orion B GMC: Revisiting its role as tracer of the dense gas reservoir for star formation

Authors: M. G. Santa-Maria, J. R. Goicoechea, J. Pety, M. Gerin, J. H. Orkisz, F. Le Petit, L. Einig, P. Palud, V. de Souza Magalhaes, I. Bešlić, L. Segal, S. Bardeau, E. Bron, P. Chainais, J. Chanussot, P. Gratier, V. V. Guzmán, A. Hughes, D. Languignon, F. Levrier, D. C. Lis, H. S. Liszt, J. Le Bourlot, Y. Oya, K. Öberg , et al. (6 additional authors not shown)

Abstract: We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V <… ▽ More We present 5 deg^2 (~250 pc^2) HCN, HNC, HCO+, and CO J=1-0 maps of the Orion B GMC, complemented with existing wide-field [CI] 492 GHz maps, as well as new pointed observations of rotationally excited HCN, HNC, H13CN, and HN13C lines. We detect anomalous HCN J=1-0 hyperfine structure line emission almost everywhere in the cloud. About 70% of the total HCN J=1-0 luminosity arises from gas at A_V < 8 mag. The HCN/CO J=1-0 line intensity ratio shows a bimodal behavior with an inflection point at A_V < 3 mag typical of translucent gas and UV-illuminated cloud edges. We find that most of the HCN J=1-0 emission arises from extended gas with n(H2) ~< 10^4 cm^-3, even lower density gas if the ionization fraction is > 10^-5 and electron excitation dominates. This result explains the low-A_V branch of the HCN/CO J=1-0 intensity ratio distribution. Indeed, the highest HCN/CO ratios (~0.1) at A_V < 3 mag correspond to regions of high [CI] 492 GHz/CO J=1-0 intensity ratios (>1) characteristic of low-density PDRs. Enhanced FUV radiation favors the formation and excitation of HCN on large scales, not only in dense star-forming clumps. The low surface brightness HCN and HCO+ J=1-0 emission scale with I_FIR (a proxy of the stellar FUV radiation field) in a similar way. Together with CO J=1-0, these lines respond to increasing I_FIR up to G0~20. On the other hand, the bright HCN J=1-0 emission from dense gas in star-forming clumps weakly responds to I_FIR once the FUV radiation field becomes too intense (G0>1500). The different power law scalings (produced by different chemistries, densities, and line excitation regimes) in a single but spatially resolved GMC resemble the variety of Kennicutt-Schmidt law indexes found in galaxy averages. As a corollary for extragalactic studies, we conclude that high HCN/CO J=1-0 line intensity ratios do not always imply the presence of dense gas. △ Less

Submitted 18 September, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: accepted for publication in A&A. 24 pages, 18 figures, plus Appendix. Abridged Abstract

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

Neural network-based emulation of interstellar medium models

Authors: Pierre Palud, Lucas Einig, Franck Le Petit, Emeric Bron, Pierre Chainais, Jocelyn Chanussot, Jérôme Pety, Pierre-Antoine Thouvenin, David Languignon, Ivana Bešlić, Miriam G. Santa-Maria, Jan H. Orkisz, Léontine E. Ségal, Antoine Zakardjian, Sébastien Bardeau, Maryvonne Gerin, Javier R. Goicoechea, Pierre Gratier, Viviana V. Guzman, Annie Hughes, François Levrier, Harvey S. Liszt, Jacques Le Bourlot, Antoine Roueff, Albrecht Sievers

Abstract: The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolatio… ▽ More The interpretation of observations of atomic and molecular tracers in the galactic and extragalactic interstellar medium (ISM) requires comparisons with state-of-the-art astrophysical models to infer some physical conditions. Usually, ISM models are too time-consuming for such inference procedures, as they call for numerous model evaluations. As a result, they are often replaced by an interpolation of a grid of precomputed models. We propose a new general method to derive faster, lighter, and more accurate approximations of the model from a grid of precomputed models. These emulators are defined with artificial neural networks (ANNs) designed and trained to address the specificities inherent in ISM models. Indeed, such models often predict many observables (e.g., line intensities) from just a few input physical parameters and can yield outliers due to numerical instabilities or physical bistabilities. We propose applying five strategies to address these characteristics: 1) an outlier removal procedure; 2) a clustering method that yields homogeneous subsets of lines that are simpler to predict with different ANNs; 3) a dimension reduction technique that enables to adequately size the network architecture; 4) the physical inputs are augmented with a polynomial transform to ease the learning of nonlinearities; and 5) a dense architecture to ease the learning of simple relations. We compare the proposed ANNs with standard classes of interpolation methods to emulate the Meudon PDR code, a representative ISM numerical model. Combinations of the proposed strategies outperform all interpolation methods by a factor of 2 on the average error, reaching 4.5% on the Meudon PDR code. These networks are also 1000 times faster than accurate interpolation methods and require ten to forty times less memory. This work will enable efficient inferences on wide-field multiline observations of the ISM. △ Less

Submitted 4 September, 2023; originally announced September 2023.

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

The extremely sharp transition between molecular and ionized gas in the Horsehead nebula

Authors: C. Hernández-Vera, V. V. Guzmán, J. R. Goicoechea, V. Maillard, J. Pety, F. Le Petit, M. Gerin, E. Bron, E. Roueff, A. Abergel, T. Schirmer, J. Carpenter, P. Gratier, K. Gordon, K. Misselt

Abstract: (Abridged) Massive stars can determine the evolution of molecular clouds with their strong ultraviolet (UV) radiation fields. Moreover, UV radiation is relevant in setting the thermal gas pressure in star-forming clouds, whose influence can extend from the rims of molecular clouds to entire star-forming galaxies. Probing the fundamental structure of nearby molecular clouds is therefore crucial to… ▽ More (Abridged) Massive stars can determine the evolution of molecular clouds with their strong ultraviolet (UV) radiation fields. Moreover, UV radiation is relevant in setting the thermal gas pressure in star-forming clouds, whose influence can extend from the rims of molecular clouds to entire star-forming galaxies. Probing the fundamental structure of nearby molecular clouds is therefore crucial to understand how massive stars shape their surrounding medium and how fast molecular clouds are destroyed, specifically at their UV-illuminated edges, where models predict an intermediate zone of neutral atomic gas between the molecular cloud and the surrounding ionized gas whose size is directly related to the exposed physical conditions. We present the highest angular resolution (~$0.5$", corresponding to $207$ au) and velocity-resolved images of the molecular gas emission in the Horsehead nebula, using CO J=3-2 and HCO$^+$ J=4-3 observations with ALMA. We find that CO and HCO$^+$ are present at the edge of the cloud, very close to the ionization (H$^+$/H) and dissociation fronts (H/H$_2$), suggesting a very thin layer of neutral atomic gas (<$650$ au) and a small amount of CO-dark gas ($A_V=0.006-0.26$ mag) for stellar UV illumination conditions typical of molecular clouds in the Milky Way. The new ALMA observations reveal a web of molecular gas filaments with an estimated thermal gas pressure of $P_{\mathrm{th}} = (2.3 - 4.0) \times 10^6$ K cm$^{-3}$, and the presence of a steep density gradient at the cloud edge that can be well explained by stationary isobaric PDR models with pressures consistent with our estimations. However, in the HII region and PDR interface, we find $P_{\mathrm{th,PDR}} > P_{\mathrm{th,HII}}$, suggesting the gas is slightly compressed. Therefore, dynamical effects cannot be completely ruled out and even higher angular observations will be needed to unveil their role. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 15 pages, 1 table, 9 figures; Accepted for publication in A&A

Journal ref: A&A 677, A152 (2023)

Baryon-meson scattering amplitude in the $1/N_c$ expansion

Authors: Victor Miguel Banda Guzmán, Ruben Flores-Mendieta, Johann Hernandez

Abstract: The baryon-meson scattering amplitude is computed within the $1/N_c$ expansion of QCD, where $N_c$ is the number of color charges. The most general expression is obtained by accounting for explicitly the effects of the decuplet-octet baryon mass difference and perturbative flavor $SU(3)$ symmetry breaking. Although the resultant expression is general enough that it can be applied to any incoming a… ▽ More The baryon-meson scattering amplitude is computed within the $1/N_c$ expansion of QCD, where $N_c$ is the number of color charges. The most general expression is obtained by accounting for explicitly the effects of the decuplet-octet baryon mass difference and perturbative flavor $SU(3)$ symmetry breaking. Although the resultant expression is general enough that it can be applied to any incoming and outgoing baryons and mesons, provided that the Gell-Mann--Nishijima scheme is respected, results for $Nπ$ scattering processes are explicitly dealt with. With these, some isospin relations are verified to be valid at the physical value $N_c=3$. The expressions obtained here represent a first effort toward understanding scattering processes in the context of the $1/N_c$ expansion. △ Less

Submitted 20 September, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: 87 pages, 1 figure

One-loop fermion-photon vertex in arbitrary gauge and dimensions: a novel approach

Authors: Victor Miguel Banda Guzmán, Adnan Bashir

Abstract: We compute one-loop electron-photon vertex with fully off-shell external momenta in an arbitrary covariant gauge and space-time dimension. There exist numerous efforts in literature where one-loop off-shell vertex is calculated by employing the standard first order Feynman rules in different covariant gauges and space-time dimensions of interest. The tensor structure which decomposes this three-po… ▽ More We compute one-loop electron-photon vertex with fully off-shell external momenta in an arbitrary covariant gauge and space-time dimension. There exist numerous efforts in literature where one-loop off-shell vertex is calculated by employing the standard first order Feynman rules in different covariant gauges and space-time dimensions of interest. The tensor structure which decomposes this three-point vertex into the components transverse and longitudinal to the photon momentum gets intertwined in this first order formalism. The Ward-Takahashi identity is explicitly invoked to untangle the two pieces and the results are expressed in a preferred basis of twelve spin-amplitudes. We propose a novel approach based upon an efficient combination of the first and second order formalisms of quantum electrodynamics to compute this one-loop vertex. Among some conspicuous advantages is the fact that this less known second order formalism separates the spin and scalar degrees of freedom of an electron interacting electromagnetically. More noticeably, the longitudinal and transverse contributions naturally disentangle from the onset in our approach. Moreover, this decomposition leads to identities between one-loop scalar Feynman integrals with higher powers in the propagators and shifted space-time dimensions that can be used to prove the Ward-Takahashi identity at one-loop order without the need to evaluate any Feynman integral. Additionally, this natural decomposition allows us to establish the gauge-independence of the Pauli form factor through explicit cancellations of scalar Feynman integrals that depend on the gauge parameter. These cancellations naturally lead to a compact expression for the Pauli form factor in arbitrary dimensions. Wherever necessary and insightful, we make comparisons with earlier works. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 14 pages, 3 Figures

Report number: JLAB-THY-23-3790

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 gap-sharing planet pair shaping the crescent in HD 163296: a disk sculpted by a resonant chain

Authors: Juan Garrido-Deutelmoser, Cristobal Petrovich, Carolina Charalambous, Viviana V. Guzmán, Ke Zhang

Abstract: ALMA observations of the disk around HD 163296 have resolved a crescent-shape substructure at around 55 au, inside and off-center from a gap in the dust that extends from 38 au to 62 au. In this work we propose that both the crescent and the dust rings are caused by a compact pair (period ratio $\simeq 4:3$) of sub-Saturn-mass planets inside the gap, with the crescent corresponding to dust trapped… ▽ More ALMA observations of the disk around HD 163296 have resolved a crescent-shape substructure at around 55 au, inside and off-center from a gap in the dust that extends from 38 au to 62 au. In this work we propose that both the crescent and the dust rings are caused by a compact pair (period ratio $\simeq 4:3$) of sub-Saturn-mass planets inside the gap, with the crescent corresponding to dust trapped at the $L_5$ Lagrange point of the outer planet. This interpretation also reproduces well the gap in the gas recently measured from the CO observations, which is shallower than what is expected in a model where the gap is carved by a single planet. Building on previous works arguing for outer planets at $\approx 86$ and $\approx 137$ au, we provide with a global model of the disk that best reproduces the data and show that all four planets may fall into a long resonant chain, with the outer three planets in a 1:2:4 Laplace resonance. We show that this configuration is not only an expected outcome from disk-planet interaction in this system, but it can also help constraining the radial and angular position of the planet candidates using three-body resonances. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 11 pages, 8 figures, re-submitted to ApJL, comments welcome

Accurate boundary-integral formulations for the calculation of electrostatic forces with an implicit-solvent model

Authors: Ian Addison-Smith, Horacio V. Guzmán, Christopher D. Cooper

Abstract: An accurate force calculation with the Poisson-Boltzmann equation is challenging, as it requires the electric field on the molecular surface. Here, we present a calculation of the electric field on the solute-solvent interface that is exact for piece-wise linear variations of the potential and analyze four different alternatives to compute the force using a boundary element method. We performed a… ▽ More An accurate force calculation with the Poisson-Boltzmann equation is challenging, as it requires the electric field on the molecular surface. Here, we present a calculation of the electric field on the solute-solvent interface that is exact for piece-wise linear variations of the potential and analyze four different alternatives to compute the force using a boundary element method. We performed a verification exercise for two cases: the isolated and two interacting molecules. Our results suggest that the boundary element method outperforms the finite difference method, as the latter needs a much finer mesh than in solvation energy calculations to get acceptable accuracy in the force, whereas the same surface mesh than a standard energy calculation is appropriate for the boundary element method. Among the four evaluated alternatives of force calculation, we saw that the most accurate one is based on the Maxwell stress tensor. However, for a realistic application, like the barnase-barstar complex, the approach based on variations of the energy functional, which is less accurate, gives equivalent results. This analysis is useful towards using the Poisson-Boltzmann equation for force calculations in applications where high accuracy is key, for example, to feed molecular dynamics models or to enable the study of the interaction between large molecular structures, like viruses adsorbed onto substrates. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 33 pages

ACM Class: J.2

Cold Deuterium Fractionation in the Nearest Planet-Forming Disk

Authors: Carlos E. Muñoz-Romero, Karin I. Öberg, Charles J. Law, Richard Teague, Yuri Aikawa, Jennifer B. Bergner, David J. Wilner, Jane Huang, Viviana V. Guzmán, L. Ilsedore Cleeves

Abstract: Deuterium fractionation provides a window to the thermal history of volatiles in the solar system and protoplanetary disks. While evidence of active molecular deuteration has been observed towards a handful of disks, it remains unclear whether this chemistry affects the composition of forming planetesimals due to limited observational constraints on the radial and vertical distribution of deuterat… ▽ More Deuterium fractionation provides a window to the thermal history of volatiles in the solar system and protoplanetary disks. While evidence of active molecular deuteration has been observed towards a handful of disks, it remains unclear whether this chemistry affects the composition of forming planetesimals due to limited observational constraints on the radial and vertical distribution of deuterated molecules. To shed light on this question, we introduce new ALMA observations of DCO$^+$ and DCN $J=2-1$ at an angular resolution of $0.5"$ (30 au) and combine them with archival data of higher energy transitions towards the protoplanetary disk around TW Hya. We carry out a radial excitation analysis assuming both LTE and non-LTE to localize the physical conditions traced by DCO$^+$ and DCN emission in the disk, thus assessing deuterium fractionation efficiencies and pathways at different disk locations. We find similar disk-averaged column densities of $1.9\times10^{12}$ and $9.8\times10^{11}$ cm$^{-2}$ for DCO$^{+}$ and DCN, with typical kinetic temperatures for both molecules of 20-30K, indicating a common origin near the comet- and planet-forming midplane. The observed DCO$^+$/DCN abundance ratio, combined with recent modeling results, provide tentative evidence of a gas phase C/O enhancement within $<40$ au. Observations of DCO$^+$ and DCN in other disks, as well as HCN and HCO$^+$, will be necessary to place the trends exhibited by TW Hya in context, and fully constrain the main deuteration mechanisms in disks. △ Less

Submitted 5 January, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: Fixed author information

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

Gas kinematics around filamentary structures in the Orion B cloud

Authors: Mathilde Gaudel, Jan H. Orkisz, Maryvonne Gerin, Jérôme Pety, Antoine Roueff, Antoine Marchal, François Levrier, Marc-Antoine Miville-Deschênes, Javier R. Goicoechea, Evelyne Roueff, Franck Le Petit, Victor de Souza Magalhaes, Pierre Palud, Miriam G. Santa-Maria, Maxime Vono, Sébastien Bardeau, Emeric Bron, Pierre Chainais, Jocelyn Chanussot, Pierre Gratier, Viviana Guzman, Annie Hughes, Jouni Kainulainen, David Languignon, Jacques Le Bourlot , et al. (5 additional authors not shown)

Abstract: Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used t… ▽ More Understanding the initial properties of star-forming material and how they affect the star formation process is key. From an observational point of view, the feedback from young high-mass stars on future star formation properties is still poorly constrained. In the framework of the IRAM 30m ORION-B large program, we obtained observations of the translucent and moderately dense gas, which we used to analyze the kinematics over a field of 5 deg^2 around the filamentary structures. We used the ROHSA algorithm to decompose and de-noise the C18O(1-0) and 13CO(1-0) signals by taking the spatial coherence of the emission into account. We produced gas column density and mean velocity maps to estimate the relative orientation of their spatial gradients. We identified three cloud velocity layers at different systemic velocities and extracted the filaments in each velocity layer. The filaments are preferentially located in regions of low centroid velocity gradients. By comparing the relative orientation between the column density and velocity gradients of each layer from the ORION-B observations and synthetic observations from 3D kinematic toy models, we distinguish two types of behavior in the dynamics around filaments: (i) radial flows perpendicular to the filament axis that can be either inflows (increasing the filament mass) or outflows and (ii) longitudinal flows along the filament axis. The former case is seen in the Orion B data, while the latter is not identified. We have also identified asymmetrical flow patterns, usually associated with filaments located at the edge of an HII region. This is the first observational study to highlight feedback from HII regions on filament formation and, thus, on star formation in the Orion B cloud. This simple statistical method can be used for any molecular cloud to obtain coherent information on the kinematics. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: 45 pages. Abridged abstract. Accepted by Astronomy and Astrophysics

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

Physical properties of accretion shocks toward the Class I protostellar system Oph-IRS 44

Authors: E. Artur de la Villarmois, V. V. Guzmán, J. K. Jørgensen, L. E. Kristensen, E. A. Bergin, D. Harsono, N. Sakai, E. F. van Dishoeck, S. Yamamoto

Abstract: (Abridged) Physical processes such as accretion shocks are thought to be common in the protostellar phase, where the envelope component is still present, and they can release molecules from the dust to the gas phase, altering the original chemical composition of the disk. Consequently, the study of accretion shocks is essential for a better understanding of the physical processes at disk scales an… ▽ More (Abridged) Physical processes such as accretion shocks are thought to be common in the protostellar phase, where the envelope component is still present, and they can release molecules from the dust to the gas phase, altering the original chemical composition of the disk. Consequently, the study of accretion shocks is essential for a better understanding of the physical processes at disk scales and their chemical output. The purpose of this work is to assess the characteristics of accretion shocks traced by sulfur-related species. We present ALMA high angular resolution observations (0.1") of the Class I protostar Oph-IRS 44. The continuum emission at 0.87 mm is observed, together with sulfur-related species such as SO, SO$_{2}$, and $^{34}$SO$_{2}$. Six lines of SO$_{2}$, two lines of $^{34}$SO$_{2}$, and one line of SO are detected toward IRS 44. The emission of all the detected lines peaks at ~0.1" (~14 au) from the continuum peak and we find infalling-rotating motions inside 30 au. However, only redshifted emission is seen between 50 and 30 au. Colder and more quiescent material is seen toward an offset region located at a distance of ~400 au from the protostar, and we do not find evidence of a Keplerian profile in these data. Accretion shocks are the most plausible explanation for the high temperatures, high densities, and velocities found for the SO$_{2}$ emission. When material enters the disk--envelope system, it generates accretion shocks that increase the dust temperature and desorb SO$_{2}$ molecules from dust grains. High-energy SO$_{2}$ transitions (~200 K) seem to be the best tracers of accretion shocks that can be followed up by future higher angular resolution ALMA observations and compared to other species to assess their importance in releasing molecules from the dust to the gas phase. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 16 pages, 16 figures

Journal ref: A&A 667, A20 (2022)

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

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

Abstract: We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is… ▽ More We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion. △ Less

Submitted 12 July, 2022; originally announced July 2022.

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

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

The Isotopic Links from Planet Forming Regions to the Solar System

Authors: H. Nomura, K. Furuya, M. A. Cordiner, S. B. Charnley, C. M. O'D. Alexander, C. A. Nixon, V. V. Guzman, H. Yurimoto, T. Tsukagoshi, T. Iino

Abstract: Isotopic ratios provide a powerful tool for understanding the origins of materials, including the volatile and refractory matter within solar system bodies. Recent high sensitivity observations of molecular isotopologues, in particular with ALMA, have brought us new information on isotopic ratios of hydrogen, carbon, nitrogen and oxygen in star and planet forming regions as well as the solar syste… ▽ More Isotopic ratios provide a powerful tool for understanding the origins of materials, including the volatile and refractory matter within solar system bodies. Recent high sensitivity observations of molecular isotopologues, in particular with ALMA, have brought us new information on isotopic ratios of hydrogen, carbon, nitrogen and oxygen in star and planet forming regions as well as the solar system objects. Solar system exploration missions, such as Rosetta and Cassini, have given us further new insights. Meanwhile, the recent development of sophisticated models for isotope chemistry including detailed gas-phase and grain surface reaction network has made it possible to discuss how isotope fractionation in star and planet forming regions is imprinted into the icy mantles of dust grains, preserving a record of the initial isotopic state of solar system materials. This chapter reviews recent progress in observations of molecular isotopologues in extra-solar planet forming regions, prestellar/protostellar cores and protoplanetary disks, as well as objects in our solar system -- comets, meteorites, and planetary/satellite atmospheres -- and discusses their connection by means of isotope chemical models. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: To appear in Protostars and Planets VII, University of Arizona Press, eds. Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Obtaining fully polarised amplitudes in gauge invariant form

Authors: Naser Ahmadiniaz, Victor Miguel Banda Guzman, Christian Schubert, Fiorenzo Bastianelli, Olindo Corradini, James P. Edwards

Abstract: We describe progress applying the \textit{Worldline Formalism} of quantum field theory to the fermion propagator dressed by $N$-photons to study multi-linear Compton scattering processes, explaining how this approach -- whose calculational advantages are well-known at multi-loop order -- yields compact and manifestly gauge invariant scattering amplitudes. We describe progress applying the \textit{Worldline Formalism} of quantum field theory to the fermion propagator dressed by $N$-photons to study multi-linear Compton scattering processes, explaining how this approach -- whose calculational advantages are well-known at multi-loop order -- yields compact and manifestly gauge invariant scattering amplitudes. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: 3 pages, 1 figure. Prepared for the Proceedings of the 20th Lomonosov Conference on Elementary Particle Physics, held in August 2021, based on a talk given by James P. Edwards on the results of [arXiv:2004.01391 [hep-th]] and [arXiv:2107.00199 [hep-th]

Thermal desorption of interstellar ices. A review on the controlling parameters and their implications fromsnowlines to chemical complexity

Authors: Marco Minissale, Yuri Aikawa, Edwin Bergin, M. Bertin, Wendy A. Brown, Stephanie Cazaux, Steven B. Charnley, Audrey Coutens, Herma M. Cuppen, Victoria Guzman, Harold Linnartz, Martin R. S. McCoustra, Albert Rimola, Johanna G. M. Schrauwen, Celine Toubin, Piero Ugliengo, Naoki Watanabe, Valentine Wakelam, Francois Dulieu

Abstract: The evolution of star-forming regions and their thermal balance are strongly influenced by their chemical composition, that, in turn, is determined by the physico-chemical processes that govern the transition between the gas phase and the solid state, specifically icy dust grains (e.g., particles adsorption and desorption). Gas-grain and grain-gas transitions as well as formation and sublimation o… ▽ More The evolution of star-forming regions and their thermal balance are strongly influenced by their chemical composition, that, in turn, is determined by the physico-chemical processes that govern the transition between the gas phase and the solid state, specifically icy dust grains (e.g., particles adsorption and desorption). Gas-grain and grain-gas transitions as well as formation and sublimation of interstellar ices are thus essential elements of understanding astrophysical observations of cold environments (e.g., pre-stellar cores) where unexpected amounts of a large variety of chemical species have been observed in the gas phase. Adsorbed atoms and molecules also undergo chemical reactions which are not efficient in the gas phase. Therefore, the parameterization of the physical properties of atoms and molecules interacting with dust grain particles is clearly a key aspect to interpret astronomical observations and to build realistic and predictive astrochemical models. In this consensus evaluation, we focus on parameters controlling the thermal desorption of ices and how these determine pathways towards molecular complexity and define the location of snowlines, which ultimately influence the planet formation process. We review different crucial aspects of desorption parameters both from a theoretical and experimental point of view. We critically assess the desorption parameters commonly used in the astrochemical community for astrophysical relevant species and provide tables with recommended values. In addition, we show that a non-trivial determination of the pre-exponential factor nu using the Transition State Theory can affect the binding energy value. Finally, we conclude this work by discussing the limitations of theoretical and experimental approaches currently used to determine the desorption properties with suggestions for future improvements. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 115 pages, 17 figures

Journal ref: ACS Earth and Space Science, 2022

Tuning contact angles of aqueous droplets on hydrophilic and hydrophobic surfaces by surfactants

Authors: Fabio Staniscia, Horacio V. Guzman, Matej Kanduč

Abstract: Adsorption of small amphiphilic molecules occurs in various biological and technological processes, sometimes desired, the other times unwanted (e.g., contamination). Surface-active molecules preferentially bind to interfaces and affect their wetting properties. We study the adsorption of short-chained alcohols (simple surfactants) to the water-vapor interface and solid surfaces of various polarit… ▽ More Adsorption of small amphiphilic molecules occurs in various biological and technological processes, sometimes desired, the other times unwanted (e.g., contamination). Surface-active molecules preferentially bind to interfaces and affect their wetting properties. We study the adsorption of short-chained alcohols (simple surfactants) to the water-vapor interface and solid surfaces of various polarities using molecular dynamics simulations. The analysis enables us to establish a theoretical expression for the adsorption coefficient, which exponentially scales with the molecular surface area and the surface wetting coefficient, and which is in good agreement with the simulation results. The competition of the adsorptions to both interfaces of a sessile droplet alters its contact angle in a nontrivial way. The influence of surfactants is strongest on very hydrophilic and very hydrophobic surfaces, whereas droplets on surfaces of moderate hydrophilicity are much less affected. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Journal ref: J. Phys. Chem. B 2022, 126, 17, 3374-3384

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) 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 (MAPS) XIX. Spiral Arms, a Tail, and Diffuse Structures Traced by CO around the GM Aur Disk

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Abstract: We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Kepleria… ▽ More We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Keplerian rotation in the projected velocity fields of both disks ($\lesssim\!5\%$ of the local Keplerian velocity), suggestive of large-scale (10s of au in size), coherent flows. By accounting for the azimuthal dependence on the projection of the velocity field, the velocity fields were decomposed into azimuthally averaged orthogonal components, $v_φ$, $v_r$ and $v_z$. Using the optically thick $^{12}$CO emission as a probe of the gas temperature, local variations of $\approx\! 3$ K ($\approx\! 5 \%$ relative changes) were observed and found to be associated with the kinematic substructures. The MWC 480 disk hosts a suite of tightly wound spiral arms. The spirals arms, in conjunction with the highly localized perturbations in the gas velocity structure (kinematic planetary signatures), indicate a giant planet, $\sim\! 1$ $M_{\rm Jup}$, at a radius of $\approx 245$ au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. (2018a) and Teague et al. (2018a) advocating for multiple $\sim\! 1$ $M_{\rm Jup}$ planets embedded in the disk. These results demonstrate that molecular line observations that characterize the dynamical structure of disks can be used to search for the signatures of embedded planets. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Abstract: Understanding the temperature structure of protoplanetary disks is key to interpreting observations, predicting the physical and chemical evolution of the disk, and modeling planet formation processes. In this study, we constrain the two-dimensional thermal structure of the disk around Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spa… ▽ More Understanding the temperature structure of protoplanetary disks is key to interpreting observations, predicting the physical and chemical evolution of the disk, and modeling planet formation processes. In this study, we constrain the two-dimensional thermal structure of the disk around Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spatially resolved ALMA observations (~0.12 arcsec (13 au) - 0.25 arcsec (26 au)) of CO J = 2-1, 13CO J = 1-0, 2-1, C18O J = 1-0, 2-1, and C17O J = 1-0, the HD J = 1-0 flux upper limit, the spectral energy distribution (SED), and continuum morphology. The final model incorporates both a radial depletion of CO motivated by a time scale shorter than typical CO gas-phase chemistry (0.01 Myr) and an enhanced temperature near the surface layer of the the inner disk (z/r <= 0.21). This model agrees with the majority of the empirically derived temperatures and observed emitting surfaces derived from the J = 2-1 CO observations. We find an upper limit for the disk mass of 0.35 Msun, using the upper limit of the HD J = 1-0 and J = 2-1 flux. With our final thermal structure, we explore the impact that gaps have on the temperature structure constrained by observations of the resolved gaps. Adding a large gap in the gas and small dust additionally increases gas temperature in the gap by only 5-10%. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

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

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

Abstract: The Molecules with ALMA at Planet-forming Scales large program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high angular resolution (0.15" and 0.30" beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3 - 1.3 mJy/beam and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our mu… ▽ More The Molecules with ALMA at Planet-forming Scales large program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high angular resolution (0.15" and 0.30" beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3 - 1.3 mJy/beam and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our multi-stage workflow -- built around the CASA tclean image deconvolution procedure -- that we used to generate the core data product of the MAPS LP: the position-position-velocity image cubes for each spectral line. Owing to the expansive nature of the survey, we encountered a range of imaging challenges; some are familiar to the sub-mm protoplanetary disk community, like the benefits of using an accurate CLEAN mask, and others less well-known, like the incorrect default flux scaling of the CLEAN residual map first described in Jorsater & van Moorsel 1995 (the "JvM effect"). We distill lessons learned into recommended workflows for synthesizing image cubes of molecular emission. In particular, we describe how to produce image cubes with accurate fluxes via the "JvM correction," a procedure that is generally applicable to any image synthesized via CLEAN deconvolution but is especially critical for low S/N emission. We further explain how we used visibility tapering to promote a common, fiducial beam size and contextualize the interpretation of signal to noise ratio when detecting molecular emission from protoplanetary disks. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. △ Less

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

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

Worldline master formulas for the dressed electron propagator, part 2: On-shell amplitudes

Authors: Naser Ahmadiniaz, Victor Miguel Banda Guzman, Fiorenzo Bastianelli, Olindo Corradini, James P. Edwards, Christian Schubert

Abstract: In the first part of this series, we employed the second-order formalism and the ``symbol'' map to construct a particle path-integral representation of the electron propagator in a background electromagnetic field, suitable for open fermion-line calculations. Its main advantages are the avoidance of long products of Dirac matrices, and its ability to unify whole sets of Feynman diagrams related by… ▽ More In the first part of this series, we employed the second-order formalism and the ``symbol'' map to construct a particle path-integral representation of the electron propagator in a background electromagnetic field, suitable for open fermion-line calculations. Its main advantages are the avoidance of long products of Dirac matrices, and its ability to unify whole sets of Feynman diagrams related by permutation of photon legs along the fermion lines. We obtained a Bern-Kosower type master formula for the fermion propagator, dressed with $N$ photons, in terms of the ``$N$-photon kernel,'' where this kernel appears also in ``subleading'' terms involving only $N$-$1$ of the $N$ photons. In this sequel, we focus on the application of the formalism to the calculation of on-shell amplitudes and cross sections. Universal formulas are obtained for the fully polarised matrix elements of the fermion propagator dressed with an arbitrary number of photons, as well as for the corresponding spin-averaged cross sections. A major simplification of the on-shell case is that the subleading terms drop out, but we also pinpoint other, less obvious simplifications. We use integration by parts to achieve manifest transversality of these amplitudes at the integrand level and exploit this property using the spinor helicity technique. We give a simple proof of the vanishing of the matrix element for ``all $+$'' photon helicities in the massless case, and find a novel relation between the scalar and spinor spin-averaged cross sections in the massive case. Testing the formalism on the standard linear Compton scattering process, we find that it reproduces the known results with remarkable efficiency. Further applications and generalisations are pointed out. △ Less

Submitted 30 June, 2021; originally announced July 2021.

Comments: 40 pages, 3 figures. Part 2 of a series started by arXiv:2004.01391 [hep-th]

An ALMA Survey of Chemistry in Disks around M4-M5 Stars

Authors: Jamila Pegues, Karin I. Oberg, Jennifer B. Bergner, Jane Huang, Ilaria Pascucci, Richard Teague, Sean M. Andrews, Edwin A. Bergin, L. Ilsedore Cleeves, Viviana V. Guzman, Feng Long, Chunhua Qi, David J. Wilner

Abstract: M-stars are the most common hosts of planetary systems in the Galaxy. Protoplanetary disks around M-stars thus offer a prime opportunity to study the chemistry of planet-forming environments. We present an ALMA survey of molecular line emission toward a sample of five protoplanetary disks around M4-M5 stars (FP Tau, J0432+1827, J1100-7619, J1545-3417, and Sz 69). These observations can resolve che… ▽ More M-stars are the most common hosts of planetary systems in the Galaxy. Protoplanetary disks around M-stars thus offer a prime opportunity to study the chemistry of planet-forming environments. We present an ALMA survey of molecular line emission toward a sample of five protoplanetary disks around M4-M5 stars (FP Tau, J0432+1827, J1100-7619, J1545-3417, and Sz 69). These observations can resolve chemical structures down to tens of AU. Molecular lines of $^{12}$CO, $^{13}$CO, C$^{18}$O, C$_2$H, and HCN are detected toward all five disks. Lines of H$_2$CO and DCN are detected toward 2/5 and 1/5 disks, respectively. For disks with resolved C$^{18}$O, C$_2$H, HCN, and H$_2$CO emission, we observe substructures similar to those previously found in disks around solar-type stars (e.g., rings, holes, and plateaus). C$_2$H and HCN excitation conditions estimated interior to the pebble disk edge for the bright disk J1100-7619 are consistent with previous measurements around solar-type stars. The correlation previously found between C$_2$H and HCN fluxes for solar-type disks extends to our M4-M5 disk sample, but the typical C$_2$H/HCN ratio is higher for the M4-M5 disk sample. This latter finding is reminiscent of the hydrocarbon enhancements found by previous observational infrared surveys in the innermost ($<$10AU) regions of M-star disks, which is intriguing since our disk-averaged fluxes are heavily influenced by flux levels in the outermost disk, exterior to the pebble disk edge. Overall, most of the observable chemistry at 10-100AU appears similar for solar-type and M4-M5 disks, but hydrocarbons may be more abundant around the cooler stars. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 46 pages (20 pages in the main document, 26 pages in the appendix), 48 figures, 7 tables. Accepted in ApJ (February 2021)

Journal ref: The Astrophysical Journal, 911(2), 150 (2021)

Contribution of $SU(3)$ quadratic Casimir squared to rotational bands in the interacting boson model

Authors: Victor Miguel Banda Guzman, Ruben Flores-Mendieta, Johann Hernandez

Abstract: Rotational bands are commonly used in the analysis of the spectra of atomic nuclei. The early version of the interacting boson model of Arima and Iachello has been foundational to the description of rotations in nuclei. The model is based on a unitary spectrum generating algebra $U(6)$ and an orthogonal (angular momentum) symmetry algebra $SO(3)$. A solvable limit of the model contains $SU(3)$ in… ▽ More Rotational bands are commonly used in the analysis of the spectra of atomic nuclei. The early version of the interacting boson model of Arima and Iachello has been foundational to the description of rotations in nuclei. The model is based on a unitary spectrum generating algebra $U(6)$ and an orthogonal (angular momentum) symmetry algebra $SO(3)$. A solvable limit of the model contains $SU(3)$ in its dynamical symmetry chain. The corresponding Hamiltonian is written as a linear combination of linear and quadratic Casimir invariants of all the algebras in the chain. Prompted by these facts, a Hamiltonian containing the $SU(3)$ quadratic Casimir squared is proposed to evaluate its effects on rotational bands. The additional term yields three undetermined parameters into the theory, which need be obtained from experiment. The lack of data does not allow one to perform a detailed numerical analysis, but a rather restricted one in terms of a one-parameter fit. Nevertheless, these additional terms provide a good description of rotational bands of two nuclei of interest, $^{156}\mathrm{Ga}$ and $^{234}\mathrm{U}$. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 13 pages, 2 figures, 1 table

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

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

Abstract: Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project… ▽ More Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at $15-20$ au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, $E_u=29-97$ K. We find the ortho-to-para ratio of $\rm C_3H_2$ is consistent with 3 throughout extent of the emission, and the total abundance of both $\rm C_3H_2$ isomers is $(7.5-10)\times10^{-11}$ per H atom, or $1-10$% of the previously published $\rm C_2H$ abundance in the same source. We find $\rm C_3H_2$ comes from a layer near the surface that extends no deeper than $z/r=0.25$. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside $\sim30$ au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridonial flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 15 pages, 8 figures, Accepted in ApJ

Exploring HNC and HCN line emission as probes of the protoplanetary disk temperature

Authors: Feng Long, Arthur D. Bosman, Paolo Cazzoletti, Ewine F. van Dishoeck, Karin I. Oberg, Stefano Facchini, Marco Tazzari, Viviana V. Guzman, Leonardo Testi

Abstract: The distributions and abundances of molecules in protoplanetary disks are powerful tracers of the physical and chemical disk structures. The abundance ratios of HCN and its isomer HNC are known to be sensitive to gas temperature. Their line ratios might therefore offer a unique opportunity to probe the properties of the emitting gas. Using the 2D thermochemical code DALI, we ran a set of disk mode… ▽ More The distributions and abundances of molecules in protoplanetary disks are powerful tracers of the physical and chemical disk structures. The abundance ratios of HCN and its isomer HNC are known to be sensitive to gas temperature. Their line ratios might therefore offer a unique opportunity to probe the properties of the emitting gas. Using the 2D thermochemical code DALI, we ran a set of disk models accounting for different stellar and disk properties, with an updated chemical network for the nitrogen chemistry. These modeling results were then compared with observations, including new ALMA observations of HNC $J=3-2$ for the TW Hya disk and HNC $J=1-0$ for 29 disks in Lupus. Similar to CN, HCN and HNC have brighter line emission in models with larger disk flaring angles and higher UV fluxes. HNC and HCN are predicted to be abundant in the warm surface layer and outer midplane region, which results in ring-shaped emission patterns. However, the precise emitting regions and emission morphology depend on the probed transition, as well as on other parameters such as C and O abundances. The modeled HNC-to-HCN line intensity ratio increases from $<0.1$ in the inner disk to up to 0.8 in the outer disk regions, which can be explained by efficient HNC destruction at high temperatures. Disk-integrated HNC line fluxes from current scarce observations and its radial distribution in the TW Hya disk are broadly consistent with our model predictions. The HNC-to-HCN flux ratio robustly increases with radius (decreasing temperature), but its use as a chemical thermometer in disks is affected by other factors, including UV flux and C and O abundances. High-spatial resolution ALMA disk observations of HNC and HCN that can locate the emitting layers would have the great potential to constrain both the disk thermal and UV radiation structures, and also to verify our understanding of the nitrogen chemistry. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: accepted for publication in A&A

Journal ref: A&A 647, A118 (2021)

Dynamical Masses and Stellar Evolutionary Model Predictions of M-Stars

Authors: Jamila Pegues, Ian Czekala, Sean M. Andrews, Karin I. Öberg, Gregory J. Herczeg, Jennifer B. Bergner, L. Ilsedore Cleeves, Viviana V. Guzmán, Jane Huang, Feng Long, Richard Teague, David J. Wilner

Abstract: In this era of Gaia and ALMA, dynamical stellar mass measurements provide benchmarks that are independent of observations of stellar characteristics and their uncertainties. These benchmarks can then be used to validate and improve stellar evolutionary models, which can lead to both imprecise and inaccurate mass predictions for pre-main-sequence, low-mass stars. We present the dynamical stellar ma… ▽ More In this era of Gaia and ALMA, dynamical stellar mass measurements provide benchmarks that are independent of observations of stellar characteristics and their uncertainties. These benchmarks can then be used to validate and improve stellar evolutionary models, which can lead to both imprecise and inaccurate mass predictions for pre-main-sequence, low-mass stars. We present the dynamical stellar masses derived from disks around three M-stars (FP Tau, J0432+1827, and J1100-7619) using ALMA observations of $^{12}$CO (J=2--1) and $^{13}$CO (J=2--1) emission. These are the first dynamical stellar mass measurements for J0432+1827 and J1100-7619 and the most precise measurement for FP Tau. Fiducial stellar evolutionary model tracks, which do not include any treatment of magnetic activity, agree with the dynamical measurement of J0432+1827 but underpredict the mass by $\sim$60\% for FP Tau and $\sim$80\% for J1100-7619. Possible explanations for the underpredictions include inaccurate assumptions of stellar effective temperature, undetected binarity for J1100-7619, and that fiducial stellar evolutionary models are not complex enough to represent these stars. In the former case, the stellar effective temperatures would need to be increased by $\sim$40K to $\sim$340K to reconcile the fiducial model predictions with the dynamically-measured masses. In the latter case, we show that the dynamical masses can be reproduced using results from stellar evolutionary models with starspots, which incorporate fractional starspot coverage to represent the manifestation of magnetic activity. Folding in low-mass M-stars from the literature and assuming that the stellar effective temperatures are imprecise but accurate, we find tentative evidence of a relationship between fractional starspot coverage and observed effective temperature for these young, cool stars. △ Less

Submitted 14 January, 2021; originally announced January 2021.

Comments: 37 pages (14 pages in the main document, 23 pages in the appendix), 35 figures, 5 tables. Accepted in ApJ (December 2020)