Showing 1–50 of 151 results for author: Cataldi, G

Tensor Networks for Lattice Gauge Theories beyond one dimension: a Roadmap

Authors: Giuseppe Magnifico, Giovanni Cataldi, Marco Rigobello, Peter Majcen, Daniel Jaschke, Pietro Silvi, Simone Montangero

Abstract: Tensor network methods are a class of numerical tools and algorithms to study many-body quantum systems in and out of equilibrium, based on tailored variational wave functions. They have found significant applications in simulating lattice gauge theories approaching relevant problems in high-energy physics. Compared to Monte Carlo methods, they do not suffer from the sign problem, allowing them to… ▽ More Tensor network methods are a class of numerical tools and algorithms to study many-body quantum systems in and out of equilibrium, based on tailored variational wave functions. They have found significant applications in simulating lattice gauge theories approaching relevant problems in high-energy physics. Compared to Monte Carlo methods, they do not suffer from the sign problem, allowing them to explore challenging regimes such as finite chemical potentials and real-time dynamics. Further development is required to tackle fundamental challenges, such as accessing continuum limits or computations of large-scale quantum chromodynamics. In this work, we review the state-of-the-art of Tensor Network methods and discuss a possible roadmap for algorithmic development and strategies to enhance their capabilities and extend their applicability to open high-energy problems. We provide tailored estimates of the theoretical and computational resource scaling for attacking large-scale lattice gauge theories. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 14 pages, 6 figures

Search for photons above 10$^{18}$ eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, L. Andrade Dourado, S. Andringa, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato , et al. (342 additional authors not shown)

Abstract: The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above $10^{17}$ eV. It measures extensive air showers generated by ultra high energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the… ▽ More The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above $10^{17}$ eV. It measures extensive air showers generated by ultra high energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum ($X_{max}$) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV (1 EeV = $10^{18}$ eV) has been developed by combining the fluorescence detector-based measurement of $X_{max}$ with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and 0.0021 km$^{-2}$ sr$^{-1}$ yr$^{-1}$ above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultra-high energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 19 pages, 22 figures

Measurement of the Depth of Maximum of Air-Shower Profiles with energies between $\mathbf{10^{18.5}}$ and $\mathbf{10^{20}}$ eV using the Surface Detector of the Pierre Auger Observatory and Deep Learning

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, L. Andrade Dourado, S. Andringa, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato , et al. (342 additional authors not shown)

Abstract: We report an investigation of the mass composition of cosmic rays with energies from 3 to 100 EeV (1 EeV=$10^{18}$ eV) using the distributions of the depth of shower maximum $X_\mathrm{max}$. The analysis relies on ${\sim}50,000$ events recorded by the Surface Detector of the Pierre Auger Observatory and a deep-learning-based reconstruction algorithm. Above energies of 5 EeV, the data set offers a… ▽ More We report an investigation of the mass composition of cosmic rays with energies from 3 to 100 EeV (1 EeV=$10^{18}$ eV) using the distributions of the depth of shower maximum $X_\mathrm{max}$. The analysis relies on ${\sim}50,000$ events recorded by the Surface Detector of the Pierre Auger Observatory and a deep-learning-based reconstruction algorithm. Above energies of 5 EeV, the data set offers a 10-fold increase in statistics with respect to fluorescence measurements at the Observatory. After cross-calibration using the Fluorescence Detector, this enables the first measurement of the evolution of the mean and the standard deviation of the $X_\mathrm{max}$ distributions up to 100 EeV. Our findings are threefold: (1.) The evolution of the mean logarithmic mass towards a heavier composition with increasing energy can be confirmed and is extended to 100 EeV. (2.) The evolution of the fluctuations of $X_\mathrm{max}$ towards a heavier and purer composition with increasing energy can be confirmed with high statistics. We report a rather heavy composition and small fluctuations in $X_\mathrm{max}$ at the highest energies. (3.) We find indications for a characteristic structure beyond a constant change in the mean logarithmic mass, featuring three breaks that are observed in proximity to the ankle, instep, and suppression features in the energy spectrum. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: submitted to Phys. Rev. D, 28 pages, 18 figures, 5 tables

Inference of the Mass Composition of Cosmic Rays with energies from $\mathbf{10^{18.5}}$ to $\mathbf{10^{20}}$ eV using the Pierre Auger Observatory and Deep Learning

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, L. Andrade Dourado, S. Andringa, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato , et al. (342 additional authors not shown)

Abstract: We present measurements of the atmospheric depth of the shower maximum $X_\mathrm{max}$, inferred for the first time on an event-by-event level using the Surface Detector of the Pierre Auger Observatory. Using deep learning, we were able to extend measurements of the $X_\mathrm{max}$ distributions up to energies of 100 EeV ($10^{20}$ eV), not yet revealed by current measurements, providing new ins… ▽ More We present measurements of the atmospheric depth of the shower maximum $X_\mathrm{max}$, inferred for the first time on an event-by-event level using the Surface Detector of the Pierre Auger Observatory. Using deep learning, we were able to extend measurements of the $X_\mathrm{max}$ distributions up to energies of 100 EeV ($10^{20}$ eV), not yet revealed by current measurements, providing new insights into the mass composition of cosmic rays at extreme energies. Gaining a 10-fold increase in statistics compared to the Fluorescence Detector data, we find evidence that the rate of change of the average $X_\mathrm{max}$ with the logarithm of energy features three breaks at $6.5\pm0.6~(\mathrm{stat})\pm1~(\mathrm{sys})$ EeV, $11\pm 2~(\mathrm{stat})\pm1~(\mathrm{sys})$ EeV, and $31\pm5~(\mathrm{stat})\pm3~(\mathrm{sys})$ EeV, in the vicinity to the three prominent features (ankle, instep, suppression) of the cosmic-ray flux. The energy evolution of the mean and standard deviation of the measured $X_\mathrm{max}$ distributions indicates that the mass composition becomes increasingly heavier and purer, thus being incompatible with a large fraction of light nuclei between 50 EeV and 100 EeV. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: submitted to Phys. Rev. Lett., 10 pages, 3 figures, 1 table

Quantum Many-Body Scarring in a Non-Abelian Lattice Gauge Theory

Authors: Giuseppe Calajò, Giovanni Cataldi, Marco Rigobello, Darvin Wanisch, Giuseppe Magnifico, Pietro Silvi, Simone Montangero, Jad C. Halimeh

Abstract: Quantum many-body scarring (QMBS) is an intriguing mechanism of ergodicity breaking that has recently spurred significant attention. Particularly prominent in Abelian lattice gauge theories (LGTs), an open question is whether QMBS nontrivially arises in non-Abelian LGTs. Here, we present evidence of robust QMBS in a non-Abelian SU(2) LGT with dynamical matter. Starting in product states that requi… ▽ More Quantum many-body scarring (QMBS) is an intriguing mechanism of ergodicity breaking that has recently spurred significant attention. Particularly prominent in Abelian lattice gauge theories (LGTs), an open question is whether QMBS nontrivially arises in non-Abelian LGTs. Here, we present evidence of robust QMBS in a non-Abelian SU(2) LGT with dynamical matter. Starting in product states that require little experimental overhead, we show that prominent QMBS arises for certain quenches, facilitated through meson and baryon-antibaryon excitations, highlighting its non-Abelian nature. The uncovered scarred dynamics manifests as long-lived coherent oscillations in experimentally accessible local observables as well as prominent revivals in the state fidelity. Our findings bring QMBS to the realm of non-Abelian LGTs, highlighting the intimate connection between scarring and gauge symmetry, and are amenable for observation in a recently proposed trapped-ion qudit quantum computer. △ Less

Submitted 24 May, 2024; v1 submitted 21 May, 2024; originally announced May 2024.

Comments: 13 pages, 8 figures

Impact of the Magnetic Horizon on the Interpretation of the Pierre Auger Observatory Spectrum and Composition Data

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato, A. Bartz Mocellin , et al. (342 additional authors not shown)

Abstract: The flux of ultra-high energy cosmic rays reaching Earth above the ankle energy (5 EeV) can be described as a mixture of nuclei injected by extragalactic sources with very hard spectra and a low rigidity cutoff. Extragalactic magnetic fields existing between the Earth and the closest sources can affect the observed CR spectrum by reducing the flux of low-rigidity particles reaching Earth. We perfo… ▽ More The flux of ultra-high energy cosmic rays reaching Earth above the ankle energy (5 EeV) can be described as a mixture of nuclei injected by extragalactic sources with very hard spectra and a low rigidity cutoff. Extragalactic magnetic fields existing between the Earth and the closest sources can affect the observed CR spectrum by reducing the flux of low-rigidity particles reaching Earth. We perform a combined fit of the spectrum and distributions of depth of shower maximum measured with the Pierre Auger Observatory including the effect of this magnetic horizon in the propagation of UHECRs in the intergalactic space. We find that, within a specific range of the various experimental and phenomenological systematics, the magnetic horizon effect can be relevant for turbulent magnetic field strengths in the local neighbourhood of order $B_{\rm rms}\simeq (50-100)\,{\rm nG}\,(20\rm{Mpc}/{d_{\rm s})( 100\,\rm{kpc}/L_{\rm coh}})^{1/2}$, with $d_{\rm s}$ the typical intersource separation and $L_{\rm coh}$ the magnetic field coherence length. When this is the case, the inferred slope of the source spectrum becomes softer and can be closer to the expectations of diffusive shock acceleration, i.e., $\propto E^{-2}$. An additional cosmic-ray population with higher source density and softer spectra, presumably also extragalactic and dominating the cosmic-ray flux at EeV energies, is also required to reproduce the overall spectrum and composition results for all energies down to 0.6~EeV. △ Less

Submitted 4 April, 2024; originally announced April 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

Testing Hadronic-Model Predictions of Depth of Maximum of Air-Shower Profiles and Ground-Particle Signals using Hybrid Data of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato, A. Bartz Mocellin , et al. (346 additional authors not shown)

Abstract: We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, $X_{max}$, and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, $S(1000)$, using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists in fitting the measured two-dimensional ($S(1000)$, $X_{max}$) distri… ▽ More We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, $X_{max}$, and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, $S(1000)$, using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists in fitting the measured two-dimensional ($S(1000)$, $X_{max}$) distributions using templates for simulated air showers produced with hadronic interaction models EPOS-LHC, QGSJet II-04, Sibyll 2.3d and leaving the scales of predicted $X_{max}$ and the signals from hadronic component at ground as free fit parameters. The method relies on the assumption that the mass composition remains the same at all zenith angles, while the longitudinal shower development and attenuation of ground signal depend on the mass composition in a correlated way. The analysis was applied to 2239 events detected by both the fluorescence and surface detectors of the Pierre Auger Observatory with energies between $10^{18.5}$ to $10^{19.0}$ eV and zenith angles below $60^\circ$. We found, that within the assumptions of the method, the best description of the data is achieved if the predictions of the hadronic interaction models are shifted to deeper $X_{max}$ values and larger hadronic signals at all zenith angles. Given the magnitude of the shifts and the data sample size, the statistical significance of the improvement of data description using the modifications considered in the paper is larger than $5σ$ even for any linear combination of experimental systematic uncertainties. △ Less

Submitted 3 May, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: Published version

Journal ref: Phys. Rev. D 109 (2024) 102001

Constraints on metastable superheavy dark matter coupled to sterile neutrinos with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato, A. Bartz Mocellin , et al. (346 additional authors not shown)

Abstract: Dark matter particles could be superheavy, provided their lifetime is much longer than the age of the universe. Using the sensitivity of the Pierre Auger Observatory to ultra-high energy neutrinos and photons, we constrain a specific extension of the Standard Model of particle physics that meets the lifetime requirement for a superheavy particle by coupling it to a sector of ultra-light sterile ne… ▽ More Dark matter particles could be superheavy, provided their lifetime is much longer than the age of the universe. Using the sensitivity of the Pierre Auger Observatory to ultra-high energy neutrinos and photons, we constrain a specific extension of the Standard Model of particle physics that meets the lifetime requirement for a superheavy particle by coupling it to a sector of ultra-light sterile neutrinos. Our results show that, for a typical dark coupling constant of 0.1, the mixing angle $θ_m$ between active and sterile neutrinos must satisfy, roughly, $θ_m \lesssim 1.5\times 10^{-6}(M_X/10^9~\mathrm{GeV})^{-2}$ for a mass $M_X$ of the dark-matter particle between $10^8$ and $10^{11}~$GeV. △ Less

Submitted 14 March, 2024; v1 submitted 24 November, 2023; originally announced November 2023.

Comments: Version accepted for publication in PRD Letters, including Supplemental Material

Journal ref: Phys. Rev. D 109 L081101 (2024)

A Primordial Origin for the Gas-Rich Debris Disks Around Intermediate-Mass Stars

Authors: Riouhei Nakatani, Neal J. Turner, Yasuhiro Hasegawa, Gianni Cataldi, Yuri Aikawa, Sebastián Marino, Hiroshi Kobayashi

Abstract: While most debris disks consist of dust with little or no gas, a fraction has significant amounts of gas detected via emission lines of CO, ionized carbon, and/or atomic oxygen. Almost all such gaseous debris disks known are around A-type stars with ages up to 50 Myr. We show, using semi-analytic disk evolution modeling, that this can be understood if the gaseous debris disks are remnant protoplan… ▽ More While most debris disks consist of dust with little or no gas, a fraction has significant amounts of gas detected via emission lines of CO, ionized carbon, and/or atomic oxygen. Almost all such gaseous debris disks known are around A-type stars with ages up to 50 Myr. We show, using semi-analytic disk evolution modeling, that this can be understood if the gaseous debris disks are remnant protoplanetary disks that have become depleted of small grains compared to the interstellar medium. Photoelectric heating by the A-stars' FUV radiation is then inefficient, while the stars' EUV and X-ray emissions are weak owing to a lack of surface convective zones capable of driving magnetic activity. In this picture, stars outside the range of spectral types from A through early B are relatively hard to have such long-lived gas disks. Less-massive stars have stronger magnetic activity in the chromosphere, transition region, and corona with resulting EUV and X-ray emission, while more-massive stars have photospheres hot enough to produce strong EUV radiation. In both cases, primordial disk gas is likely to photoevaporate well before 50 Myr. These results come from 0D disk evolution models where we incorporate internal accretion stresses, MHD winds, and photoevaporation by EUV and X-ray photons with luminosities that are functions of the stellar mass and age. A key issue this work leaves open is how some disks become depleted in small dust so that FUV photoevaporation slows. Candidates include grains' growth, settling, radial drift, radiation force, and incorporation into planetary systems. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: Submitted to ApJL

Radio Measurements of the Depth of Air-Shower Maximum at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, Anukriti, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato , et al. (350 additional authors not shown)

Abstract: The Auger Engineering Radio Array (AERA), part of the Pierre Auger Observatory, is currently the largest array of radio antenna stations deployed for the detection of cosmic rays, spanning an area of $17$ km$^2$ with 153 radio stations. It detects the radio emission of extensive air showers produced by cosmic rays in the $30-80$ MHz band. Here, we report the AERA measurements of the depth of the s… ▽ More The Auger Engineering Radio Array (AERA), part of the Pierre Auger Observatory, is currently the largest array of radio antenna stations deployed for the detection of cosmic rays, spanning an area of $17$ km$^2$ with 153 radio stations. It detects the radio emission of extensive air showers produced by cosmic rays in the $30-80$ MHz band. Here, we report the AERA measurements of the depth of the shower maximum ($X_\text{max}$), a probe for mass composition, at cosmic-ray energies between $10^{17.5}$ to $10^{18.8}$ eV, which show agreement with earlier measurements with the fluorescence technique at the Pierre Auger Observatory. We show advancements in the method for radio $X_\text{max}$ reconstruction by comparison to dedicated sets of CORSIKA/CoREAS air-shower simulations, including steps of reconstruction-bias identification and correction, which is of particular importance for irregular or sparse radio arrays. Using the largest set of radio air-shower measurements to date, we show the radio $X_\text{max}$ resolution as a function of energy, reaching a resolution better than $15$ g cm$^{-2}$ at the highest energies, demonstrating that radio $X_\text{max}$ measurements are competitive with the established high-precision fluorescence technique. In addition, we developed a procedure for performing an extensive data-driven study of systematic uncertainties, including the effects of acceptance bias, reconstruction bias, and the investigation of possible residual biases. These results have been cross-checked with air showers measured independently with both the radio and fluorescence techniques, a setup unique to the Pierre Auger Observatory. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: Submitted to Phys. Rev. D

Demonstrating Agreement between Radio and Fluorescence Measurements of the Depth of Maximum of Extensive Air Showers at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, Anukriti, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, F. Barbato , et al. (350 additional authors not shown)

Abstract: We show, for the first time, radio measurements of the depth of shower maximum ($X_\text{max}$) of air showers induced by cosmic rays that are compared to measurements of the established fluorescence method at the same location. Using measurements at the Pierre Auger Observatory we show full compatibility between our radio and the previously published fluorescence data set, and between a subset of… ▽ More We show, for the first time, radio measurements of the depth of shower maximum ($X_\text{max}$) of air showers induced by cosmic rays that are compared to measurements of the established fluorescence method at the same location. Using measurements at the Pierre Auger Observatory we show full compatibility between our radio and the previously published fluorescence data set, and between a subset of air showers observed simultaneously with both radio and fluorescence techniques, a measurement setup unique to the Pierre Auger Observatory. Furthermore, we show radio $X_\text{max}$ resolution as a function of energy and demonstrate the ability to make competitive high-resolution $X_\text{max}$ measurements with even a sparse radio array. With this, we show that the radio technique is capable of cosmic-ray mass composition studies, both at Auger and at other experiments. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: Submitted to Phys. Rev. Lett

Scalable digital quantum simulation of lattice fermion theories with local encoding

Authors: Marco Ballarin, Giovanni Cataldi, Giuseppe Magnifico, Daniel Jaschke, Marco Di Liberto, Ilaria Siloi, Simone Montangero, Pietro Silvi

Abstract: We numerically analyze the feasibility of a platform-neutral, general strategy to perform quantum simulations of fermionic lattice field theories under open boundary conditions. The digital quantum simulator requires solely one- and two-qubit gates and is scalable since integrating each Hamiltonian term requires a finite (non-scaling) cost. The exact local fermion encoding we adopt relies on auxil… ▽ More We numerically analyze the feasibility of a platform-neutral, general strategy to perform quantum simulations of fermionic lattice field theories under open boundary conditions. The digital quantum simulator requires solely one- and two-qubit gates and is scalable since integrating each Hamiltonian term requires a finite (non-scaling) cost. The exact local fermion encoding we adopt relies on auxiliary $\mathbb{Z}_2$ lattice gauge fields by adding a pure gauge Hamiltonian term akin to the Toric Code. By numerically emulating the quantum simulator real-time dynamics, we observe a timescale separation for spin- and charge-excitations in a spin-$\frac{1}{2}$ Hubbard ladder in the $t-J$ model limit. Additionally, we show that such local fermion encoding is also useful in tensor network numerical simulations. △ Less

Submitted 5 December, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 17 pages, 10 figures

Ground observations of a space laser for the assessment of its in-orbit performance

Authors: The Pierre Auger Collaboration, O. Lux, I. Krisch, O. Reitebuch, D. Huber, D. Wernham, T. Parrinello, :, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, Anukriti, L. Apollonio, C. Aramo, P. R. Araújo Ferreira , et al. (358 additional authors not shown)

Abstract: The wind mission Aeolus of the European Space Agency was a groundbreaking achievement for Earth observation. Between 2018 and 2023, the space-borne lidar instrument ALADIN onboard the Aeolus satellite measured atmospheric wind profiles with global coverage which contributed to improving the accuracy of numerical weather prediction. The precision of the wind observations, however, declined over the… ▽ More The wind mission Aeolus of the European Space Agency was a groundbreaking achievement for Earth observation. Between 2018 and 2023, the space-borne lidar instrument ALADIN onboard the Aeolus satellite measured atmospheric wind profiles with global coverage which contributed to improving the accuracy of numerical weather prediction. The precision of the wind observations, however, declined over the course of the mission due to a progressive loss of the atmospheric backscatter signal. The analysis of the root cause was supported by the Pierre Auger Observatory in Argentina whose fluorescence detector registered the ultraviolet laser pulses emitted from the instrument in space, thereby offering an estimation of the laser energy at the exit of the instrument for several days in 2019, 2020 and 2021. The reconstruction of the laser beam not only allowed for an independent assessment of the Aeolus performance, but also helped to improve the accuracy in the determination of the laser beam's ground track on single pulse level. The results presented in this paper set a precedent for the monitoring of space lasers by ground-based telescopes and open new possibilities for the calibration of cosmic-ray observatories. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 10 pages, 10 figures

The Pierre Auger Observatory Open Data

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, Anukriti, L. Apollonio, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova , et al. (352 additional authors not shown)

Abstract: The Pierre Auger Collaboration has embraced the concept of open access to their research data since its foundation, with the aim of giving access to the widest possible community. A gradual process of release began as early as 2007 when 1% of the cosmic-ray data was made public, along with 100% of the space-weather information. In February 2021, a portal was released containing 10% of cosmic-ray d… ▽ More The Pierre Auger Collaboration has embraced the concept of open access to their research data since its foundation, with the aim of giving access to the widest possible community. A gradual process of release began as early as 2007 when 1% of the cosmic-ray data was made public, along with 100% of the space-weather information. In February 2021, a portal was released containing 10% of cosmic-ray data collected from 2004 to 2018, during Phase I of the Observatory. The Portal included detailed documentation about the detection and reconstruction procedures, analysis codes that can be easily used and modified and, additionally, visualization tools. Since then the Portal has been updated and extended. In 2023, a catalog of the 100 highest-energy cosmic-ray events examined in depth has been included. A specific section dedicated to educational use has been developed with the expectation that these data will be explored by a wide and diverse community including professional and citizen-scientists, and used for educational and outreach initiatives. This paper describes the context, the spirit and the technical implementation of the release of data by the largest cosmic-ray detector ever built, and anticipates its future developments. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: 19 pages, 8 figures

AugerPrime Surface Detector Electronics

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, Anukriti, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, F. Barbato , et al. (346 additional authors not shown)

Abstract: Operating since 2004, the Pierre Auger Observatory has led to major advances in our understanding of the ultra-high-energy cosmic rays. The latest findings have revealed new insights that led to the upgrade of the Observatory, with the primary goal of obtaining information on the primary mass of the most energetic cosmic rays on a shower-by-shower basis. In the framework of the upgrade, called Aug… ▽ More Operating since 2004, the Pierre Auger Observatory has led to major advances in our understanding of the ultra-high-energy cosmic rays. The latest findings have revealed new insights that led to the upgrade of the Observatory, with the primary goal of obtaining information on the primary mass of the most energetic cosmic rays on a shower-by-shower basis. In the framework of the upgrade, called AugerPrime, the 1660 water-Cherenkov detectors of the surface array are equipped with plastic scintillators and radio antennas, allowing us to enhance the composition sensitivity. To accommodate new detectors and to increase experimental capabilities, the electronics is also upgraded. This includes better timing with up-to-date GPS receivers, higher sampling frequency, increased dynamic range, and more powerful local processing of the data. In this paper, the design characteristics of the new electronics and the enhanced dynamic range will be described. The manufacturing and test processes will be outlined and the test results will be discussed. The calibration of the SD detector and various performance parameters obtained from the analysis of the first commissioning data will also be presented. △ Less

Submitted 8 October, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

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

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

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

Submitted 22 August, 2023; originally announced August 2023.

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

Search for UHE Photons from Gravitational Wave Sources with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato , et al. (346 additional authors not shown)

Abstract: A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localizati… ▽ More A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at $\sim$7 MeV cm$^{-2}$ (time period 1000~s) and $\sim$35 MeV cm$^{-2}$ (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total gravitational wave energy. These are the first limits on UHE photons from GW sources. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 15 pages

Journal ref: ApJ 952 (2023) 91

(2+1)D SU(2) Yang-Mills Lattice Gauge Theory at finite density via tensor networks

Authors: Giovanni Cataldi, Giuseppe Magnifico, Pietro Silvi, Simone Montangero

Abstract: We numerically simulate a non-Abelian lattice gauge theory in two spatial dimensions, with Tensor Networks (TN). We focus on the SU(2) Yang-Mills model in Hamiltonian formulation, with dynamical matter and minimally truncated gauge field (hardcore gluon). Thanks to the TN sign-problem-free approach, we characterize the phase diagram of the model at zero and finite baryon number as a function of th… ▽ More We numerically simulate a non-Abelian lattice gauge theory in two spatial dimensions, with Tensor Networks (TN). We focus on the SU(2) Yang-Mills model in Hamiltonian formulation, with dynamical matter and minimally truncated gauge field (hardcore gluon). Thanks to the TN sign-problem-free approach, we characterize the phase diagram of the model at zero and finite baryon number as a function of the quark bare mass and color charge. Already at intermediate system sizes, we distinctly detect a liquid phase of quark-pair bound-state quasi-particles (baryons), whose mass is finite towards the continuum limit. Interesting phenomena arise at the transition boundary where color-electric and color-magnetic terms are maximally frustrated: for low quark masses, we see traces of potential deconfinement, while for high masses, signatures of a possible topological order. △ Less

Submitted 8 January, 2024; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 24 pages, 16 figures. Added references; extended Sec. 1, Sec. 2, Appendix A, and Appendix B with generalized and scalable formalism for arbitrary gauge field truncations

Constraining models for the origin of ultra-high-energy cosmic rays with a novel combined analysis of arrival directions, spectrum, and composition data measured at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu , et al. (349 additional authors not shown)

Abstract: The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearb… ▽ More The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around $20^\circ$ for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of $4.5σ$ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the $γ$-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions. △ Less

Submitted 14 January, 2024; v1 submitted 26 May, 2023; originally announced May 2023.

Journal ref: JCAP01(2024)022

Primordial or Secondary? Testing models of debris disk gas with ALMA

Authors: Gianni Cataldi, Yuri Aikawa, Kazunari Iwasaki, Sebastian Marino, Alexis Brandeker, Antonio Hales, Thomas Henning, Aya E. Higuchi, A. Meredith Hughes, Markus Janson, Quentin Kral, Luca Matrà, Attila Moór, Göran Olofsson, Seth Redfield, Aki Roberge

Abstract: The origin and evolution of gas in debris disks is still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only few C observations available. We create an overview of the C and CO content of debris disk gas and use it test state-of-the-art models. We use new and archival ALMA obs… ▽ More The origin and evolution of gas in debris disks is still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only few C observations available. We create an overview of the C and CO content of debris disk gas and use it test state-of-the-art models. We use new and archival ALMA observations of CO and CI emission, complemented by CII data from Herschel, for a sample of 14 debris disks. This expands the number of disks with ALMA measurements of both CO and CI by ten disks. We present new detections of CI emission towards three disks: HD 21997, HD 121191 and HD 121617. We use a simple disk model to derive gas masses and column densities. We find that current state-of-the-art models of secondary gas production overpredict the neutral carbon content of debris disk gas. This does not rule out a secondary origin, but might indicate that the models require an additional C removal process. Alternatively, the gas might be produced in transient events rather than a steady-state collisional cascade. We also test a primordial gas origin by comparing our results to a simplified thermo-chemical model. This yields promising results, but more detailed work is required before a conclusion can be reached. Our work demonstrates that the combination of C and CO data is a powerful tool to advance our understanding of debris disk gas. △ Less

Submitted 5 June, 2023; v1 submitted 20 May, 2023; originally announced May 2023.

Comments: 90 pages, 60 figures. Accepted for publication in ApJ. version 2: additional acknowledgement. versions 3, 4: minor edits

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 Catalog of the Highest-Energy Cosmic Rays Recorded During Phase I of Operation of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, P. Allison, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, M. Ave, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova , et al. (354 additional authors not shown)

Abstract: A catalog containing details of the highest-energy cosmic rays recorded through the detection of extensive air-showers at the Pierre Auger Observatory is presented with the aim of opening the data to detailed examination. Descriptions of the 100 showers created by the highest-energy particles recorded between 1 January 2004 and 31 December 2020 are given for cosmic rays that have energies in the r… ▽ More A catalog containing details of the highest-energy cosmic rays recorded through the detection of extensive air-showers at the Pierre Auger Observatory is presented with the aim of opening the data to detailed examination. Descriptions of the 100 showers created by the highest-energy particles recorded between 1 January 2004 and 31 December 2020 are given for cosmic rays that have energies in the range 78 EeV to 166 EeV. Details are also given of a further nine very-energetic events that have been used in the calibration procedure adopted to determine the energy of each primary. A sky plot of the arrival directions of the most energetic particles is shown. No interpretations of the data are offered. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Article accepted for publication in ApJS

Constraining the sources of ultra-high-energy cosmic rays across and above the ankle with the spectrum and composition data measured at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato , et al. (343 additional authors not shown)

Abstract: In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above $6 \times 10^{17}$ eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around $5\times 10^{18}$ eV (the so-called "ankle" feature) as a transition between these two components. We fin… ▽ More In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above $6 \times 10^{17}$ eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around $5\times 10^{18}$ eV (the so-called "ankle" feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources. To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs). △ Less

Submitted 17 April, 2023; v1 submitted 5 November, 2022; originally announced November 2022.

Comments: Submitted to JCAP

Journal ref: JCAP05(2023)024

Searches for Ultra-High-Energy Photons at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, J. A. Bellido , et al. (340 additional authors not shown)

Abstract: The Pierre Auger Observatory, being the largest air-shower experiment in the world, offers an unprecedented exposure to neutral particles at the highest energies. Since the start of data taking more than 18 years ago, various searches for ultra-high-energy (UHE, $E\gtrsim10^{17}\,\text{eV}$) photons have been performed: either for a diffuse flux of UHE photons, for point sources of UHE photons or… ▽ More The Pierre Auger Observatory, being the largest air-shower experiment in the world, offers an unprecedented exposure to neutral particles at the highest energies. Since the start of data taking more than 18 years ago, various searches for ultra-high-energy (UHE, $E\gtrsim10^{17}\,\text{eV}$) photons have been performed: either for a diffuse flux of UHE photons, for point sources of UHE photons or for UHE photons associated with transient events like gravitational wave events. In the present paper, we summarize these searches and review the current results obtained using the wealth of data collected by the Pierre Auger Observatory. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: Review article accepted for publication in Universe (special issue on ultra-high energy photons)

Report number: FERMILAB-PUB-22-810-AD-PPD-SCD-TD

Journal ref: Universe 8 (2022) 579

Mid-infrared time-domain study of recent dust production events in the extreme debris disc of TYC 4209-1322-1

Authors: A. Moór, P. Ábrahám, Á. Kóspál, K. Y. L. Su, G. H. Rieke, G. Cataldi, A. Bódi, Zs. Bognár, B. Cseh, G. Csörnyei, N. Egei, A. Farkas, O. Hanyecz, B. Ignácz, Cs. Kalup, R. Könyves-Tóth, L. Kriskovics, L. Mészáros, A. Pál, A. Ordasi, K. Sárneczky, B. Seli, Á. Sódor, R. Szakáts, J. Vinkó , et al. (1 additional authors not shown)

Abstract: Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescop… ▽ More Extreme debris discs are characterized by unusually strong mid-infrared excess emission, which often proves to be variable. The warm dust in these discs is of transient nature and is likely related to a recent giant collision occurring close to the star in the terrestrial region. Here we present the results of a 877 days long, gap-free photometric monitoring performed by the Spitzer Space Telescope of the recently discovered extreme debris disc around TYC 4209-1322-1. By combining these observations with other time-domain optical and mid-infrared data, we explore the disc variability of the last four decades with particular emphasis on the last 12 yr. During the latter interval the disc showed substantial changes, the most significant was the brightening and subsequent fading between 2014 and 2018 as outlined in WISE data. The Spitzer light curves outline the fading phase and a subsequent new brightening of the disc after 2018, revealing an additional flux modulation with a period of ~39 days on top of the long-term trend. We found that all these variations can be interpreted as the outcome of a giant collision that happened at an orbital radius of ~0.3 au sometime in 2014. Our analysis implies that a collision on a similar scale could have taken place around 2010, too. The fact that the disc was already peculiarly dust rich 40 yr ago, as implied by IRAS data, suggests that these dust production events belong to a chain of large impacts triggered by an earlier even more catastrophic collision. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 12 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 516, Issue 4, November 2022, Pages 5684-5701

Search for photons above 10$^{19}$ eV with the surface detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, J. A. Bellido , et al. (343 additional authors not shown)

Abstract: We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above $10^{19}$ eV. Photons in the zenith angle range from 30$^\circ$ to 60$^\circ$ can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the arr… ▽ More We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above $10^{19}$ eV. Photons in the zenith angle range from 30$^\circ$ to 60$^\circ$ can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the array and the steeper lateral distribution of shower particles reaching ground. Applying the search method to data collected between January 2004 and June 2020, upper limits at 95\% CL are set to an $E^{-2}$ diffuse flux of ultra-high energy photons above $10^{19}$ eV, $2{\times}10^{19}$ eV and $4{\times}10^{19}$ eV amounting to $2.11{\times}10^{-3}$, $3.12{\times}10^{-4}$ and $1.72{\times}10^{-4}$ km$^{-2}$ sr$^{-1}$ yr$^{-1}$, respectively. While the sensitivity of the present search around $2 \times 10^{19}$ eV approaches expectations of cosmogenic photon fluxes in the case of a pure-proton composition, it is one order of magnitude above those from more realistic mixed-composition models. The inferred limits have also implications for the search of super-heavy dark matter that are discussed and illustrated. △ Less

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

Comments: 21 pages, 8 figures

Cosmological implications of photon-flux upper limits at ultra-high energies in scenarios of Planckian-interacting massive particles for dark matter

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu , et al. (352 additional authors not shown)

Abstract: Using the data of the Pierre Auger Observatory, we report on a search for signatures that would be suggestive of super-heavy particles decaying in the Galactic halo. From the lack of signal, we present upper limits for different energy thresholds above ${\gtrsim}10^8$\,GeV on the secondary by-product fluxes expected from the decay of the particles. Assuming that the energy density of these super-h… ▽ More Using the data of the Pierre Auger Observatory, we report on a search for signatures that would be suggestive of super-heavy particles decaying in the Galactic halo. From the lack of signal, we present upper limits for different energy thresholds above ${\gtrsim}10^8$\,GeV on the secondary by-product fluxes expected from the decay of the particles. Assuming that the energy density of these super-heavy particles matches that of dark matter observed today, we translate the upper bounds on the particle fluxes into tight constraints on the couplings governing the decay process as a function of the particle mass. Instantons, which are non-perturbative solutions to Yang-Mills equations, can give rise to decay channels otherwise forbidden and transform stable particles into meta-stable ones. Assuming such instanton-induced decay processes, we derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $α_X \lesssim 0.09$, for $10^{9} \lesssim M_X/\text{GeV} < 10^{19}$. Conversely, we obtain that, for instance, a reduced coupling constant $α_X = 0.09$ excludes masses $M_X \gtrsim 3\times 10^{13}~$GeV. In the context of dark matter production from gravitational interactions alone during the reheating epoch, we derive constraints on the parameter space that involves, in addition to $M_X$ and $α_X$, the Hubble rate at the end of inflation, the reheating efficiency, and the non-minimal coupling of the Higgs with curvature. △ Less

Submitted 15 December, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 15 pages, 8 figures, Accompanying paper of arXiv:2203.08854, accepted for publication in Phys. Rev. D

Journal ref: Phys.Rev.D 107 (2023) 4, 042002

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/

Arrival Directions of Cosmic Rays above 32 EeV from Phase One of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu , et al. (350 additional authors not shown)

Abstract: A promising energy range to look for angular correlation between cosmic rays of extragalactic origin and their sources is at the highest energies, above few tens of EeV ($1\:{\rm EeV}\equiv 10^{18}\:$eV). Despite the flux of these particles being extremely low, the area of ${\sim}\:3{,}000 \: \text{km}^2$ covered at the Pierre Auger Observatory, and the 17-year data-taking period of the Phase 1 of… ▽ More A promising energy range to look for angular correlation between cosmic rays of extragalactic origin and their sources is at the highest energies, above few tens of EeV ($1\:{\rm EeV}\equiv 10^{18}\:$eV). Despite the flux of these particles being extremely low, the area of ${\sim}\:3{,}000 \: \text{km}^2$ covered at the Pierre Auger Observatory, and the 17-year data-taking period of the Phase 1 of its operations, have enabled us to measure the arrival directions of more than 2,600 ultra-high energy cosmic rays above $32\:\text{EeV}$. We publish this data set, the largest available at such energies from an integrated exposure of $122{,}000 \: \text{km}^2\:\text{sr}\:\text{yr}$, and search it for anisotropies over the $3.4π$ steradians covered with the Observatory. Evidence for a deviation in excess of isotropy at intermediate angular scale, with ${\sim}\:15^\circ$ Gaussian spread or ${\sim}\:25^\circ$ top-hat radius, is obtained at the $4\:σ$ significance level for cosmic-ray energies above ${\sim}\:40\:\text{EeV}$. △ Less

Submitted 5 September, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: Published version, 33 pages, 10 figures. Dataset and analysis software available at https://doi.org/10.5281/zenodo.6504276

Report number: FERMILAB-PUB-22-491-AD-PPD-SCD-TD

Journal ref: ApJ 935 170 (2022)

Vertical evolution of exocometary gas: I. How vertical diffusion shortens the CO lifetime

Authors: S. Marino, G. Cataldi, M. R. Jankovic, L. Matrà, M. C. Wyatt

Abstract: Bright debris discs can contain large amounts of CO gas. This gas was thought to be a protoplanetary remnant until it was recently shown that it could be released in collisions of volatile-rich solids. As CO is released, interstellar UV radiation photodissociates CO producing CI, which can shield CO allowing a large CO mass to accumulate. However, this picture was challenged because CI is ineffici… ▽ More Bright debris discs can contain large amounts of CO gas. This gas was thought to be a protoplanetary remnant until it was recently shown that it could be released in collisions of volatile-rich solids. As CO is released, interstellar UV radiation photodissociates CO producing CI, which can shield CO allowing a large CO mass to accumulate. However, this picture was challenged because CI is inefficient at shielding if CO and CI are vertically mixed. Here, we study for the first time the vertical evolution of gas to determine how vertical mixing affects the efficiency of shielding by CI. We present a 1D model that accounts for gas release, photodissociation, ionisation, viscous evolution, and vertical mixing due to turbulent diffusion. We find that if the gas surface density is high and the vertical diffusion weak ($α_{\rm v}/α<[H/r]^2$) CO photodissociates high above the midplane, forming an optically thick CI layer that shields the CO underneath. Conversely, if diffusion is strong ($α_{\rm v}/α>[H/r]^2$) CI and CO become well mixed, shortening the CO lifetime. Moreover, diffusion could also limit the amount of dust settling. High-resolution ALMA observations could resolve the vertical distribution of CO and CI, and thus constrain vertical mixing and the efficiency of CI shielding. We also find that the CO and CI scale heights may not be good probes of the mean molecular weight, and thus composition, of the gas. Finally, we show that if mixing is strong the CO lifetime might not be long enough for CO to spread interior to the planetesimal belt where gas is produced. △ Less

Submitted 27 June, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS (18 pages, 14 figures)

Investigating Hadronic Interactions at Ultra-High Energies with the Pierre Auger Observatory

Authors: Isabel Goos, :, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova , et al. (352 additional authors not shown)

Abstract: The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid de… ▽ More The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles that arrive at the ground. This way, observables that are sensitive to hadronic interactions at ultra-high energies can be obtained. While the hadronic interaction cross-section can be assessed from the longitudinal profiles, the number of muons and their fluctuations measured with the ground detectors are linked to other physical properties. In addition to these direct studies, we discuss here how measurements of the atmospheric depth of the maximum of air-shower profiles and the characteristics of the muon signal at the ground can be used to test the self-consistency of the post-LHC hadronic models. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: Proceedings of VHEPU 56th Rencontres de Moriond 2022

A search for photons with energies above $2{\times}10^{17}$ eV using hybrid data from the low-energy extensions of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu , et al. (351 additional authors not shown)

Abstract: Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding… ▽ More Ultra-high-energy photons with energies exceeding $10^{17}$ eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the $10^{15}$ eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding $2{\times}10^{17}$ eV using about 5.5 years of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between $10^{17}$ and $10^{18}$ eV. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: Accepted for publication in ApJ

Report number: FERMILAB-PUB-22-441-AD-PPD-SCD-TD

Journal ref: ApJ 933 (2022) 125

Limits to gauge coupling in the dark sector set by the non-observation of instanton-induced decay of Super-Heavy Dark Matter in the Pierre Auger Observatory data

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, J. Ammerman Yebra, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, E. Avocone, A. M. Badescu, A. Bakalova, A. Balaceanu , et al. (352 additional authors not shown)

Abstract: Instantons, which are non-perturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decayi… ▽ More Instantons, which are non-perturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decaying in the Galactic halo. These particles could have been produced during the post-inflationary epoch and match the relic abundance of dark matter inferred today. The non-observation of the signatures searched for allows us to derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $α_X \lesssim 0.09$, for $10^{9} \lesssim M_X/{\rm GeV} < 10^{19}$. Conversely, we obtain that, for instance, a reduced coupling constant $α_X = 0.09$ excludes masses $M_X \gtrsim 3\times 10^{13}~$GeV. In the context of dark matter production from gravitational interactions alone, we illustrate how these bounds are complementary to those obtained on the Hubble rate at the end of inflation from the non-observation of tensor modes in the cosmological microwave background. △ Less

Submitted 15 December, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 9 pages, 3 figures, accepted for publication in Phys. Rev. Lett

Journal ref: Phys.Rev.Lett. 130 (2023) 6, 061001

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

Authors: The ANTARES collaboration, A. Albert, S. Alves, M. André, M. Anghinolfi, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, B. Belhorma, M. Bendahman, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (1025 additional authors not shown)

Abstract: For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for corre… ▽ More For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs. △ Less

Submitted 23 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 39 pages, 7 figures, 4 tables; updated source files including xml authorlist

Report number: FERMILAB-PUB-22-033-AD-PPD-SCD-TD

Journal ref: ApJ 934 164 (2022)

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, K. Almeida Cheminant, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, E. Arnone, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, J. A. Bellido , et al. (352 additional authors not shown)

Abstract: Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constrai… ▽ More Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+δ_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $δ_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $δ_{γ,0} > -10^{-21}$, $δ_{γ,1} > -10^{-40}$ eV$^{-1}$ and $δ_{γ,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $δ_{\mathrm{had},0} < 10^{-19}$, $δ_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $δ_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$σ$ CL. △ Less

Submitted 19 January, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Comments: Published version, 22 pages, 8 figures, 2 tables

Report number: FERMILAB-PUB-21-716-AD-PPD-SCD-TD

Journal ref: JCAP 01 (2022) 023

The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, J. M. Albury, I. Allekotte, A. Almela, J. Alvarez-Muñiz, R. Alves Batista, G. A. Anastasi, L. Anchordoqui, B. Andrada, S. Andringa, C. Aramo, P. R. Araújo Ferreira, J. C. Arteaga Velázquez, H. Asorey, P. Assis, G. Avila, A. M. Badescu, A. Bakalova, A. Balaceanu, F. Barbato, R. J. Barreira Luz, K. H. Becker, J. A. Bellido , et al. (352 additional authors not shown)

Abstract: We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this sp… ▽ More We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays. △ Less

Submitted 20 April, 2022; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: Published version, 25 pages, 20 figures, 11 tables Updated April 2022 to fix typo in Table 10

Report number: FERMILAB-PUB-21-474-AD-AE-SCD-TD

Journal ref: Eur. Phys. J. C (2021) 81:966

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Comments: accepted to ApJS, 33 pages, 20 figures

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

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

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

Submitted 13 September, 2021; originally announced September 2021.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Comments: 21 pages, 18 figures

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

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

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

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

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

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

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

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

Submitted 13 September, 2021; originally announced September 2021.

Comments: 18 pages, 12 figures

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

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

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

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

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