Showing 1–50 of 100 results for author: Schafer, C

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

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

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

Submitted 9 July, 2024; originally announced July 2024.

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

The flux of ultra-high-energy cosmic rays along the supergalactic plane 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, 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: Ultra-high-energy cosmic rays are known to be mainly of extragalactic origin, and their propagation is limited by energy losses, so their arrival directions are expected to correlate with the large-scale structure of the local Universe. In this work, we investigate the possible presence of intermediate-scale excesses in the flux of the most energetic cosmic rays from the direction of the supergala… ▽ More Ultra-high-energy cosmic rays are known to be mainly of extragalactic origin, and their propagation is limited by energy losses, so their arrival directions are expected to correlate with the large-scale structure of the local Universe. In this work, we investigate the possible presence of intermediate-scale excesses in the flux of the most energetic cosmic rays from the direction of the supergalactic plane region using events with energies above 20 EeV recorded with the surface detector array of the Pierre Auger Observatory up to 31 December 2022, with a total exposure of 135,000 km^2 sr yr. The strongest indication for an excess that we find, with a post-trial significance of 3.1σ, is in the Centaurus region, as in our previous reports, and it extends down to lower energies than previously studied. We do not find any strong hints of excesses from any other region of the supergalactic plane at the same angular scale. In particular, our results do not confirm the reports by the Telescope Array collaboration of excesses from two regions in the Northern Hemisphere at the edge of the field of view of the Pierre Auger Observatory. With a comparable exposure, our results in those regions are in good agreement with the expectations from an isotropic distribution. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: submitted to ApJ

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

Overstability of the 2:1 mean motion resonance: Exploring disc parameters with hydrodynamic simulations

Authors: Zahra Afkanpour, Sareh Ataiee, Alexandros Ziampras, Anna B. T. Penzlin, Rafael Sfair, Christoph Schäfer, Wilhelm Kley, Hilke Schlichting

Abstract: Resonant planetary migration in protoplanetary discs can lead to an interplay between the resonant interaction of planets and their disc torques called overstability. While theoretical predictions and N-body simulations hinted at its existence, there was no conclusive evidence until hydrodynamical simulations were performed. Our primary purpose is to find a hydrodynamic setup that induces overst… ▽ More Resonant planetary migration in protoplanetary discs can lead to an interplay between the resonant interaction of planets and their disc torques called overstability. While theoretical predictions and N-body simulations hinted at its existence, there was no conclusive evidence until hydrodynamical simulations were performed. Our primary purpose is to find a hydrodynamic setup that induces overstability in a planetary system with two moderate-mass planets in a first-order 2:1 mean motion resonance. We also aim to analyse the impact of key disc parameters, namely the viscosity, surface density, and aspect ratio, on the occurrence of overstability in this planetary system when the masses of the planets are kept constant. We performed 2D locally isothermal hydrodynamical simulations of two planets, with masses of 5 and 10 $M_{\oplus}$, in a 2:1 resonance. Upon identifying the fiducial model in which the system exhibits overstability, we performed simulations with different disc parameters to explore the effects of the disc on the overstability of the system. We observe an overstable planetary system in our hydrodynamic simulations. In the parameter study, we note that overstability occurs in discs characterised by low surface density and low viscosity. Increasing the surface density reduces the probability of overstability within the system. A limit cycle was observed in a specific viscous model with $α_ν = 10^{-3}$. In almost all our models, planets create partial gaps in the disc, which affects both the migration timescale and structure of the planetary system. We demonstrate the existence of overstability using hydrodynamic simulations but find deviations from the analytic approximation and show that the main contribution to this deviation can be attributed to dynamic gap opening. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 16 pages (including appendices), 27 figures, 3 tables, accepted for publication by A&A

Journal ref: A&A 686, A277 (2024)

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.

Multiplicative and additive systematics in galaxy density fluctuations and clustering measurements

Authors: Federico Berlfein, Rachel Mandelbaum, Scott Dodelson, Chad Schafer

Abstract: Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical… ▽ More Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical uncertainties get smaller, correcting for these systematic contaminants becomes increasingly important for unbiased cosmological analysis. We present and validate a new method for understanding and mitigating these systematics in galaxy clustering measurements (two-point function) by identifying and characterizing contaminants in the galaxy overdensity field (one-point function) using a maximum-likelihood estimator (MLE). We test this methodology with KiDS-like mock galaxy catalogs and synthetic systematic template maps. We estimate the cosmological impact of such mitigation by quantifying uncertainties and possible biases in the inferred relationship between the observed and the true galaxy clustering signal. Our method robustly corrects the clustering signal to the sub-percent level and reduces numerous additive and multiplicative systematics from 1.5$σ$ to less than 0.1$σ$ for the scenarios we tested. In addition, we provide an empirical approach to identifying the functional form (additive, multiplicative, or other) by which specific systematics contaminate the galaxy number density. Even though this approach is tested and geared towards systematics contaminating the galaxy number density, the methods can be extended to systematics mitigation for other two-point correlation measurements. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 19 pages, 10 figures; submitted to MNRAS

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)

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

Reinterpreting Fundamental Plane Correlations with Machine Learning

Authors: Chad Schafer, Sukhdeep Singh, Yesukhei Jagvaral

Abstract: This work explores the relationships between galaxy sizes and related observable galaxy properties in a large volume cosmological hydrodynamical simulation. The objectives of this work are to both develop a better understanding of the correlations between galaxy properties and the influence of environment on galaxy physics in order to build an improved model for the galaxy sizes, building off of t… ▽ More This work explores the relationships between galaxy sizes and related observable galaxy properties in a large volume cosmological hydrodynamical simulation. The objectives of this work are to both develop a better understanding of the correlations between galaxy properties and the influence of environment on galaxy physics in order to build an improved model for the galaxy sizes, building off of the {\it fundamental plane}. With an accurate intrinsic galaxy size predictor, the residuals in the observed galaxy sizes can potentially be used for multiple cosmological applications, including making measurements of galaxy velocities in spectroscopic samples, estimating the rate of cosmic expansion, and constraining the uncertainties in the photometric redshifts of galaxies. Using projection pursuit regression, the model accurately predicts intrinsic galaxy sizes and have residuals which have limited correlation with galaxy properties. The model decreases the spatial correlation of galaxy size residuals by a factor of $\sim$ 5 at small scales compared to the baseline correlation when the mean size is used as a predictor. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 16 pages, 12 figures, MNRAS

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

Stability of coorbital planets around binaries

Authors: Stefan Adelbert, Anna B. T. Penzlin, Christoph M. Schäfer, Wilhelm Kley, Billy Quarles, Rafael Sfair

Abstract: In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the para… ▽ More In previous hydrodynamical simulations, we found a mechanism for nearly circular binary stars, like Kepler-413, to trap two planets in a stable 1:1 resonance. Therefore, the stability of coorbital configurations becomes a relevant question for planet formation around binary stars. Here, we investigate the coorbital planet stability using a Kepler-413 analogue as example and then expanding the parameters to study general n-body stability of planet pairs in eccentric horseshoe orbits around binaries. The stability is tested by evolving the planet orbits for $10^5$ binary periods with varying initial semi-major axes and planet eccentricities. The unstable region of a single circumbinary planet is used as a comparison to the investigated coorbital configurations in this work. We confirm previous findings on the stability of single planets and find a first order linear relation between orbit eccentricity and pericentre to identify stable orbits for various binary configurations. Such a linear relation is also found for the stability of 1:1 resonant planets around binaries. Stable orbits for eccentric horseshoe configurations exist with a pericentre closer than seven binary separations and, in the case of Kepler-413, the pericentre of the first stable orbit can be approximated by $r_{c,peri} = (2.88 e_p + 2.46) a_{bin}$. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 10 figures, accepted for publication by A&A

Journal ref: A&A 680, A29 (2023)

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.

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

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

Observability of Photoevaporation Signatures in the Dust Continuum Emission of Transition Discs

Authors: Giovanni Picogna, Carolina Schäfer, Barbara Ercolano, Christian Rab, Rafael Franz, Matías Gárate

Abstract: Photoevaporative disc winds play a key role in our understanding of circumstellar disc evolution, especially in the final stages, and they might affect the planet formation process and the final location of planets. The study of transition discs (i.e. discs with a central dust cavity) is central for our understanding of the photoevaporation process and disc dispersal. However, we need to distingui… ▽ More Photoevaporative disc winds play a key role in our understanding of circumstellar disc evolution, especially in the final stages, and they might affect the planet formation process and the final location of planets. The study of transition discs (i.e. discs with a central dust cavity) is central for our understanding of the photoevaporation process and disc dispersal. However, we need to distinguish cavities created by photoevaporation from those created by giant planets. Theoretical models are necessary to identify possible observational signatures of the two different processes, and models to find the differences between the two processes are still lacking. In this paper we study a sample of transition discs obtained from radiation-hydrodynamic simulations of internally photoevaporated discs, and focus on the dust dynamics relevant for current ALMA observations. We then compared our results with gaps opened by super Earths/giant planets, finding that the photoevaporated cavity steepness depends mildly on gap size, and it is similar to that of a 1 Jupiter mass planet. However, the dust density drops less rapidly inside the photoevaporated cavity compared to the planetary case due to the less efficient dust filtering. This effect is visible in the resulting spectral index, which shows a larger spectral index at the cavity edge and a shallower increase inside it with respect to the planetary case. The combination of cavity steepness and spectral index might reveal the true nature of transition discs. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 10 pages, 7 figures, accepted for publication in MNRAS

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

Residual Neural Networks for the Prediction of Planetary Collision Outcomes

Authors: Philip M. Winter, Christoph Burger, Sebastian Lehner, Johannes Kofler, Thomas I. Maindl, Christoph M. Schäfer

Abstract: Fast and accurate treatment of collisions in the context of modern N-body planet formation simulations remains a challenging task due to inherently complex collision processes. We aim to tackle this problem with machine learning (ML), in particular via residual neural networks. Our model is motivated by the underlying physical processes of the data-generating process and allows for flexible predic… ▽ More Fast and accurate treatment of collisions in the context of modern N-body planet formation simulations remains a challenging task due to inherently complex collision processes. We aim to tackle this problem with machine learning (ML), in particular via residual neural networks. Our model is motivated by the underlying physical processes of the data-generating process and allows for flexible prediction of post-collision states. We demonstrate that our model outperforms commonly used collision handling methods such as perfect inelastic merging and feed-forward neural networks in both prediction accuracy and out-of-distribution generalization. Our model outperforms the current state of the art in 20/24 experiments. We provide a dataset that consists of 10164 Smooth Particle Hydrodynamics (SPH) simulations of pairwise planetary collisions. The dataset is specifically suited for ML research to improve computational aspects for collision treatment and for studying planetary collisions in general. We formulate the ML task as a multi-task regression problem, allowing simple, yet efficient training of ML models for collision treatment in an end-to-end manner. Our models can be easily integrated into existing N-body frameworks and can be used within our chosen parameter space of initial conditions, i.e. where similar-sized collisions during late-stage terrestrial planet formation typically occur. △ Less

Submitted 9 October, 2022; originally announced October 2022.

Comments: 13 pages, 7 figures, 7 tables

MSC Class: 70F16 ACM Class: E.1; I.6.6; I.2.0

Analytical models for the design of photophoretically levitating macroscopic sensors in the stratosphere

Authors: Benjamin C. Schafer, Jong-hyoung Kim, Joost J. Vlassak, David W. Keith

Abstract: Photophoretic forces could levitate thin 10 centimeter-scale structures in Earth's stratosphere indefinitely. We develop analytical models of the thermal transpiration force on a bilayer sandwich structure in the stratosphere. Lofting is maximized when the layers are separated by an air gap equal to the mean free path (MFP), when about half of the layers' surfaces consist of holes with radii < MFP… ▽ More Photophoretic forces could levitate thin 10 centimeter-scale structures in Earth's stratosphere indefinitely. We develop analytical models of the thermal transpiration force on a bilayer sandwich structure in the stratosphere. Lofting is maximized when the layers are separated by an air gap equal to the mean free path (MFP), when about half of the layers' surfaces consist of holes with radii < MFP, and when the top layer is solar-transmissive and infrared-emissive while the bottom layer is solar-absorptive and infrared-transmissive. We use the models to design a 10 cm diameter device with sufficient strength to withstand forces that might be encountered in transport, deployment, and flight. The device has a payload of about 285 mg at an altitude of 25 km; enough to support bidirectional radio communication at over 10 Mb/s and limited navigation. Such devices could be useful for atmospheric science or telecommunications on Earth and Mars. Structures a few times larger might have payloads of a few grams. △ Less

Submitted 26 November, 2022; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: Main: 17 pages, 5 figures. Supplementary information: 8 pages, 8 figures

Effects of impact and target parameters on the results of a kinetic impactor: predictions for the Double Asteroid Redirection Test (DART) mission

Authors: Angela M. Stickle, Mallory E. DeCoster, Christoph Burger, Wendy K. Caldwell, Dawn Graninger, Kathryn M. Kumamoto, Robert Luther, Jens Ormö, Sabina Raducan, Emma Rainey, Christoph M. Schäfer, James D. Walker, Yun Zhang, Patrick Michel, J. Michael Owen, Olivier Barnouin, Andy F. Cheng, Sidney Cochron, Gareth S. Collins, Thomas M. Davison, Elisabetta Dotto, Fabio Ferrari, M. Isabel Herreros, Stavro L. Ivanovski, Martin Jutzi , et al. (8 additional authors not shown)

Abstract: The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material pro… ▽ More The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material properties. Currently, very little is known about Dimorphos and its material properties that introduces uncertainty in the results of the deflection efficiency observables, including crater formation, ejecta distribution, and Beta. The DART Impact Modeling Working Group (IWG) is responsible for using impact simulations to better understand the results of the DART impact. Pre-impact simulation studies also provide considerable insight into how different properties and impact scenarios affect momentum enhancement following a kinetic impact. This insight provides a basis for predicting the effects of the DART impact and the first understanding of how to interpret results following the encounter. Following the DART impact, the knowledge gained from these studies will inform the initial simulations that will recreate the impact conditions, including providing estimates for potential material properties of Dimorphos and Beta resulting from DARTs impact. This paper summarizes, at a high level, what has been learned from the IWG simulations and experiments in preparation for the DART impact. While unknown, estimates for reasonable potential material properties of Dimorphos provide predictions for Beta of 1-5, depending on end-member cases in the strength regime. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted to PSJ Didymos-DART Focus Issue

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

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)

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

Binary orbital evolution driven by a circumbinary disc

Authors: Anna B. T. Penzlin, Wilhelm Kley, Hugo Audiffren, Christoph M. Schäfer

Abstract: The question whether the interaction of a circumbinary disc with the central binary system leads to shrinking or expansion of the binary orbit has attracted considerable interest as it impacts the evolution of binary black holes and stellar binary stars in their formation phase. We performed two-dimensional hydrodynamical simulations of circumbinary discs for a large parameter set of disc viscosit… ▽ More The question whether the interaction of a circumbinary disc with the central binary system leads to shrinking or expansion of the binary orbit has attracted considerable interest as it impacts the evolution of binary black holes and stellar binary stars in their formation phase. We performed two-dimensional hydrodynamical simulations of circumbinary discs for a large parameter set of disc viscosities and thicknesses and two different binary mass ratios for binaries on circular orbits. For those we measured carefully the net angular momentum and mass transfer between disc and binary system, and evaluate the normalised specific angular momentum accretion, $j_\mathrm{s}$ . This is compared to the theoretical, critical specific angular momentum change $j_\mathrm{s,crit}$ that separates contracting from expanding cases which depends on the the binary's mass ratio and the relative accretion onto the two stars. Using finite and infinite disc models we show that the inferred binary evolution is very similar for both setups and confirm that $j_\mathrm{s}$ can be measured accurately with cylindrical simulations that do not include the central binary. However, to obtain the relative accretion onto the stars for non-equal mass binaries, simulations that cover the whole domain including the binary are required. We find that for thick discs with aspect ratio $h = 0.1$ the binaries expand for all viscosities, while discs with $h = 0.05$ lead to an expansion only for larger viscosities with $α$ exceeding $\sim 0.005$. Overall, the regime of binary expansion extends to a much wider parameter space than previously anticipated, but for thin, low viscosity discs the orbits shrink. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 16 pages, 23 figures, accepted by A&A

Journal ref: A&A 660, A101 (2022)

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

Protoplanet Collisions: Statistical Properties of Ejecta

Authors: Samuele Crespi, Ian Dobbs-Dixon, Nikolaos Georgakarakos, Nader Haghighipour, Thomas I. Maindl, Christoph M. Schäfer, Philip Matthias Winter

Abstract: The last phase of the formation of rocky planets is dominated by collisions among Moon- to Mars-sized planetary embryos. Simulations of this phase need to handle the difficulty of including the post-impact material without saturating the numerical integrator. A common approach is to include the collision-generated material by clustering it into few bodies with the same mass and uniformly scatterin… ▽ More The last phase of the formation of rocky planets is dominated by collisions among Moon- to Mars-sized planetary embryos. Simulations of this phase need to handle the difficulty of including the post-impact material without saturating the numerical integrator. A common approach is to include the collision-generated material by clustering it into few bodies with the same mass and uniformly scattering them around the collision point. However, this approach oversimplifies the properties of the collision material by neglecting features that can play important roles in the final structure and composition of the system. In this study, we present a statistical analysis of the orbital architecture, mass, and size distributions of the material generated through embryo-embryo collisions and show how they can be used to develop a model that can be directly incorporated into the numerical integrations. For instance, results of our analysis indicate that the masses of the fragments follow an exponential distribution with an exponent of $-2.21\pm0.17$ over the range of $10^{-7}$ to $2\times 10^{-2}$ Earth-masses. The distribution of the post-impact velocities show that a large number of fragments are scattered toward the central star. The latter is a new finding that may be quite relevant to the delivery of material from the outer regions of the asteroid belt to the accretion zones of terrestrial planets. Finally, we present an analytical model for the 2D distribution of fragments that can be directly incorporated into numerical integrations. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 10 pages, 13 figures, accepted for publication in MNRAS

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

Design and implementation of the AMIGA embedded system for data acquisition

Authors: The Pierre Auger Collaboration, A. Aab, 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 , et al. (361 additional authors not shown)

Abstract: The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km$^2$ large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down t… ▽ More The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km$^2$ large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down to about 10$^{17}$ eV. At the depth of 2.3 m the electromagnetic component of cosmic ray showers is almost entirely absorbed so that the buried scintillators provide an independent and direct measurement of the air showers muon content. This work describes the design and implementation of the AMIGA embedded system, which provides centralized control, data acquisition and environment monitoring to its detectors. The presented system was firstly tested in the engineering array phase ended in 2017, and lately selected as the final design to be installed in all new detectors of the production phase. The system was proven to be robust and reliable and has worked in a stable manner since its first deployment. △ Less

Submitted 20 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: Accepted for publication at JINST. Published version, 34 pages, 15 figures, 4 tables

Report number: FERMILAB-PUB-21-032-AD-AE-E-TD

Journal ref: JINST 16 T07008 (2021)

The FRAM robotic telescope for atmospheric monitoring at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, 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 , et al. (355 additional authors not shown)

Abstract: FRAM (F/Photometric Robotic Atmospheric Monitor) is a robotic telescope operated at the Pierre Auger Observatory in Argentina for the purposes of atmospheric monitoring using stellar photometry. As a passive system which does not produce any light that could interfere with the observations of the fluorescence telescopes of the observatory, it complements the active monitoring systems that use lase… ▽ More FRAM (F/Photometric Robotic Atmospheric Monitor) is a robotic telescope operated at the Pierre Auger Observatory in Argentina for the purposes of atmospheric monitoring using stellar photometry. As a passive system which does not produce any light that could interfere with the observations of the fluorescence telescopes of the observatory, it complements the active monitoring systems that use lasers. We discuss the applications of stellar photometry for atmospheric monitoring at optical observatories in general and the particular modes of operation employed by the Auger FRAM. We describe in detail the technical aspects of FRAM, the hardware and software requirements for a successful operation of a robotic telescope for such a purpose and their implementation within the FRAM system. △ Less

Submitted 26 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: Published version, 27 pages, 8 figures

Journal ref: JINST 16 P06027 (2021)

Deep-Learning based Reconstruction of the Shower Maximum $X_{\mathrm{max}}$ using the Water-Cherenkov Detectors of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, 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 , et al. (348 additional authors not shown)

Abstract: The atmospheric depth of the air shower maximum $X_{\mathrm{max}}$ is an observable commonly used for the determination of the nuclear mass composition of ultra-high energy cosmic rays. Direct measurements of $X_{\mathrm{max}}$ are performed using observations of the longitudinal shower development with fluorescence telescopes. At the same time, several methods have been proposed for an indirect e… ▽ More The atmospheric depth of the air shower maximum $X_{\mathrm{max}}$ is an observable commonly used for the determination of the nuclear mass composition of ultra-high energy cosmic rays. Direct measurements of $X_{\mathrm{max}}$ are performed using observations of the longitudinal shower development with fluorescence telescopes. At the same time, several methods have been proposed for an indirect estimation of $X_{\mathrm{max}}$ from the characteristics of the shower particles registered with surface detector arrays. In this paper, we present a deep neural network (DNN) for the estimation of $X_{\mathrm{max}}$. The reconstruction relies on the signals induced by shower particles in the ground based water-Cherenkov detectors of the Pierre Auger Observatory. The network architecture features recurrent long short-term memory layers to process the temporal structure of signals and hexagonal convolutions to exploit the symmetry of the surface detector array. We evaluate the performance of the network using air showers simulated with three different hadronic interaction models. Thereafter, we account for long-term detector effects and calibrate the reconstructed $X_{\mathrm{max}}$ using fluorescence measurements. Finally, we show that the event-by-event resolution in the reconstruction of the shower maximum improves with increasing shower energy and reaches less than $25~\mathrm{g/cm^{2}}$ at energies above $2\times 10^{19}~\mathrm{eV}$. △ Less

Submitted 27 July, 2021; v1 submitted 8 January, 2021; originally announced January 2021.

Comments: Published version, 29 pages, 12 figures

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

Journal ref: JINST 16 P07019 (2021)

Calibration of the underground muon detector of the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, 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 , et al. (348 additional authors not shown)

Abstract: To obtain direct measurements of the muon content of extensive air showers with energy above $10^{16.5}$ eV, the Pierre Auger Observatory is currently being equipped with an underground muon detector (UMD), consisting of 219 10 $\mathrm{m^2}$-modules, each segmented into 64 scintillators coupled to silicon photomultipliers (SiPMs). Direct access to the shower muon content allows for the study of b… ▽ More To obtain direct measurements of the muon content of extensive air showers with energy above $10^{16.5}$ eV, the Pierre Auger Observatory is currently being equipped with an underground muon detector (UMD), consisting of 219 10 $\mathrm{m^2}$-modules, each segmented into 64 scintillators coupled to silicon photomultipliers (SiPMs). Direct access to the shower muon content allows for the study of both of the composition of primary cosmic rays and of high-energy hadronic interactions in the forward direction. As the muon density can vary between tens of muons per m$^2$ close to the intersection of the shower axis with the ground to much less than one per m$^2$ when far away, the necessary broad dynamic range is achieved by the simultaneous implementation of two acquisition modes in the read-out electronics: the binary mode, tuned to count single muons, and the ADC mode, suited to measure a high number of them. In this work, we present the end-to-end calibration of the muon detector modules: first, the SiPMs are calibrated by means of the binary channel, and then, the ADC channel is calibrated using atmospheric muons, detected in parallel to the shower data acquisition. The laboratory and field measurements performed to develop the implementation of the full calibration chain of both binary and ADC channels are presented and discussed. The calibration procedure is reliable to work with the high amount of channels in the UMD, which will be operated continuously, in changing environmental conditions, for several years. △ Less

Submitted 14 April, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: Published version, 25 pages, 12 figures, 0 tables

Report number: FERMILAB-PUB-20-678-AD-E-TD

Journal ref: JINST (2021) 16 P04003

Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, 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, 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. (335 additional authors not shown)

Abstract: AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each det… ▽ More AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m$^2$ detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m$^2$. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector. △ Less

Submitted 25 January, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 40 pages, 33 figures

Report number: FERMILAB-PUB-20-605-AD-E-TD

Journal ref: JINST 16 (2021) P01026

A search for ultra high energy neutrinos from TXS 0506+056 using the Pierre Auger Observatory

Authors: The Pierre Auger Collaboration, A. Aab, 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, 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. (342 additional authors not shown)

Abstract: Results of a search for ultra-high-energy neutrinos with the Pierre Auger Observatory from the direction of the blazar TXS 0506+056 are presented. They were obtained as part of the follow-up that stemmed from the detection of high-energy neutrinos and gamma rays with IceCube, \textit{Fermi}-LAT, MAGIC, and other detectors of electromagnetic radiation in several bands. The Pierre Auger Observatory… ▽ More Results of a search for ultra-high-energy neutrinos with the Pierre Auger Observatory from the direction of the blazar TXS 0506+056 are presented. They were obtained as part of the follow-up that stemmed from the detection of high-energy neutrinos and gamma rays with IceCube, \textit{Fermi}-LAT, MAGIC, and other detectors of electromagnetic radiation in several bands. The Pierre Auger Observatory is sensitive to neutrinos in the energy range from 100 PeV to 100 EeV and in the zenith angle range from $θ=60^\circ$ to $θ=95^\circ$, where the zenith angle is measured from the vertical direction. No neutrinos from the direction of TXS 0506+056 have been found. The results were analyzed in three periods: One of 6 months around the detection of IceCube-170922A, coinciding with a flare period of TXS 0506+056, a second one of 110 days during which the IceCube collaboration found an excess of 13 neutrinos from a direction compatible with TXS 0506+056, and a third one from 2004 January 1 up to 2018 August 31, over which the Pierre Auger Observatory has been taking data. The sensitivity of the Observatory is addressed for different spectral indices by considering the fluxes that would induce a single expected event during the observation period. For indices compatible with those measured by the IceCube collaboration the expected number of neutrinos at the Observatory is well-below one. Spectral indices as hard as 1.5 would have to apply in this energy range to expect a single event to have been detected. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: Published version, 4 figures

Journal ref: Published in The Astrophysical Journal, 902:105, 2020

Features of the energy spectrum of cosmic rays above $2.5{\times} 10^{18}$ eV using the Pierre Auger Observatory

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

Abstract: We report a measurement of the energy spectrum of cosmic rays above $2.5{\times} 10^{18}$ eV based on $215,030$ events. New results are presented: at about $1.3{\times} 10^{19}$ eV, the spectral index changes from $2.51 \pm 0.03 \textrm{ (stat.)} \pm 0.05 \textrm{ (sys.)}$ to $3.05 \pm 0.05 \textrm{ (stat.)}\pm 0.10\textrm{ (sys.)}$, evolving to… ▽ More We report a measurement of the energy spectrum of cosmic rays above $2.5{\times} 10^{18}$ eV based on $215,030$ events. New results are presented: at about $1.3{\times} 10^{19}$ eV, the spectral index changes from $2.51 \pm 0.03 \textrm{ (stat.)} \pm 0.05 \textrm{ (sys.)}$ to $3.05 \pm 0.05 \textrm{ (stat.)}\pm 0.10\textrm{ (sys.)}$, evolving to $5.1\pm0.3\textrm{ (stat.)} \pm 0.1\textrm{ (sys.)}$ beyond $5{\times} 10^{19}$ eV, while no significant dependence of spectral features on the declination is seen in the accessible range. These features of the spectrum can be reproduced in models with energy-dependent mass composition. The energy density in cosmic rays above $5{\times} 10^{18}$ eV is $(5.66 \pm 0.03 \textrm{ (stat.)} \pm 1.40 \textrm{ (sys.)} ) {\times} 10^{53}~$erg Mpc$^{-3}$. △ Less

Submitted 6 October, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Published version; 2 figures, 1 table

Report number: FERMILAB-PUB-20-431-E-TD

Journal ref: Phys. Rev. Lett. 125, 121106 (2020)

Measurement of the cosmic-ray energy spectrum above $2.5{\times} 10^{18}$ eV using the Pierre Auger Observatory

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

Abstract: We report a measurement of the energy spectrum of cosmic rays for energies above $2.5 {\times} 10^{18}~$eV based on 215,030 events recorded with zenith angles below $60^\circ$. A key feature of the work is that the estimates of the energies are independent of assumptions about the unknown hadronic physics or of the primary mass composition. The measurement is the most precise made hitherto with th… ▽ More We report a measurement of the energy spectrum of cosmic rays for energies above $2.5 {\times} 10^{18}~$eV based on 215,030 events recorded with zenith angles below $60^\circ$. A key feature of the work is that the estimates of the energies are independent of assumptions about the unknown hadronic physics or of the primary mass composition. The measurement is the most precise made hitherto with the accumulated exposure being so large that the measurements of the flux are dominated by systematic uncertainties except at energies above $5 {\times} 10^{19}~$eV. The principal conclusions are: (1) The flattening of the spectrum near $5 {\times} 10^{18}~$eV, the so-called "ankle", is confirmed. (2) The steepening of the spectrum at around $5 {\times} 10^{19}~$eV is confirmed. (3) A new feature has been identified in the spectrum: in the region above the ankle the spectral index $γ$ of the particle flux ($\propto E^{-γ}$) changes from $2.51 \pm 0.03~{\rm (stat.)} \pm 0.05~{\rm (sys.)}$ to $3.05 \pm 0.05~{\rm (stat.)} \pm 0.10~{\rm (sys.)}$ before changing sharply to $5.1 \pm 0.3~{\rm (stat.)} \pm 0.1~{\rm (sys.)}$ above $5 {\times} 10^{19}~$eV. (4) No evidence for any dependence of the spectrum on declination has been found other than a mild excess from the Southern Hemisphere that is consistent with the anisotropy observed above $8 {\times} 10^{18}~$eV. △ Less

Submitted 6 October, 2020; v1 submitted 14 August, 2020; originally announced August 2020.

Comments: Published version; 16 figures, 6 table

Report number: FERMILAB-PUB-20-432-E-TD

Journal ref: Phys. Rev. D 102, 062005 (2020)

Reconstruction of Events Recorded with the Surface Detector of the Pierre Auger Observatory

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

Abstract: Cosmic rays arriving at Earth collide with the upper parts of the atmosphere, thereby inducing extensive air showers. When secondary particles from the cascade arrive at the ground, they are measured by surface detector arrays. We describe the methods applied to the measurements of the surface detector of the Pierre Auger Observatory to reconstruct events with zenith angles less than $60^\circ$ us… ▽ More Cosmic rays arriving at Earth collide with the upper parts of the atmosphere, thereby inducing extensive air showers. When secondary particles from the cascade arrive at the ground, they are measured by surface detector arrays. We describe the methods applied to the measurements of the surface detector of the Pierre Auger Observatory to reconstruct events with zenith angles less than $60^\circ$ using the timing and signal information recorded using the water-Cherenkov detector stations. In addition, we assess the accuracy of these methods in reconstructing the arrival directions of the primary cosmic ray particles and the sizes of the induced showers. △ Less

Submitted 5 November, 2020; v1 submitted 17 July, 2020; originally announced July 2020.

Comments: Published version, 36 pages, 14 figures

Report number: FERMILAB-PUB-20-487-AD-E-TD

Journal ref: JINST 15 (2020) P10021

Studies on the response of a water-Cherenkov detector of the Pierre Auger Observatory to atmospheric muons using an RPC hodoscope

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

Abstract: Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowled… ▽ More Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure. △ Less

Submitted 9 September, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: Published version, 25 pages, 9 figures

Report number: FERMILAB-PUB-20-285-PPD-TD

Journal ref: JINST 15 (2020) P09002

Galaxy cluster cores as seen with VLT/MUSE: new strong-lensing analyses of RX J2129.4+0009, MS 0451.6-0305 & MACSJ2129.4-0741

Authors: Mathilde Jauzac, Baptiste Klein, Jean-Paul Kneib, Johan Richard, Markus Rexroth, Christoph Schäfer, Aurélien Verdier

Abstract: We present strong-lensing analyses of three galaxy clusters, RXJ2129.4+0009 (z=0.235), MS0451.6-0305 (z=0.55), and MACSJ2129.4-0741 (z=0.589), using the powerful combination of Hubble Space Telescope (HST) multi-band observations, and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. In RXJ2129, we newly spectroscopically confirm 15 cluster members. Our resulting mass model uses 8 multiple im… ▽ More We present strong-lensing analyses of three galaxy clusters, RXJ2129.4+0009 (z=0.235), MS0451.6-0305 (z=0.55), and MACSJ2129.4-0741 (z=0.589), using the powerful combination of Hubble Space Telescope (HST) multi-band observations, and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. In RXJ2129, we newly spectroscopically confirm 15 cluster members. Our resulting mass model uses 8 multiple image systems as we include a galaxy-galaxy lensing system North-East of the cluster, and is composed of 71 halos including one dark matter cluster-scale halo and 2 galaxy-scale halos optimized individually. For MS0451, we report the spectroscopic identification of 2 new systems of multiple images in the Northern region, and 112 cluster members. Our mass model uses 16 multiple image systems, and 146 halos, including 2 large-scale halos, and 7 galaxy-scale halos independently optimized. For MACSJ2129, we report the spectroscopic identification of one new multiple image system at z=4.41, and newly measure spectroscopic redshifts for 4 cluster members. Our mass model uses 14 multiple image systems, and is composed of 151 halos, including 2 large-scale halos and 4 galaxy-scale halos independently optimized. Our best models have rms of 0.29'', 0.6'', 0.74'' in the image plane for RXJ2129, MS0451, and MACSJ2129 respectively. This analysis presents a detailed comparison with the existing literature showing excellent agreements, and discuss specific studies of lensed galaxies, e.g. a group of submilimeter galaxies at z=2.9 in MS0451, and a bright z=2.1472 red singly imaged galaxy in MACSJ2129. △ Less

Submitted 3 August, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 25 pages, 10 figures, 13 tables, accepted in MNRAS

Search for magnetically-induced signatures in the arrival directions of ultra-high-energy cosmic rays measured at the Pierre Auger Observatory

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

Abstract: We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same… ▽ More We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory. △ Less

Submitted 27 July, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

Comments: Published version; correction of a minor typo

Report number: FERMILAB-PUB-20-177-E-TD

Journal ref: Published in JCAP06(2020)017

Lenstool-HPC: A High Performance Computing based mass modelling tool for cluster-scale gravitational lenses

Authors: Christoph Schäfer, Gilles Fourestey, Jean-Paul Kneib

Abstract: With the upcoming generation of telescopes, cluster scale strong gravitational lenses will act as an increasingly relevant probe of cosmology and dark matter. The better resolved data produced by current and future facilities requires faster and more efficient lens modeling software. Consequently, we present Lenstool-HPC, a strong gravitational lens modeling and map generation tool based on High… ▽ More With the upcoming generation of telescopes, cluster scale strong gravitational lenses will act as an increasingly relevant probe of cosmology and dark matter. The better resolved data produced by current and future facilities requires faster and more efficient lens modeling software. Consequently, we present Lenstool-HPC, a strong gravitational lens modeling and map generation tool based on High Performance Computing (HPC) techniques and the renowned Lenstool software. We also showcase the HPC concepts needed for astronomers to increase computation speed through massively parallel execution on supercomputers. Lenstool-HPC was developed using lens modelling algorithms with high amounts of parallelism. Each algorithm was implemented as a highly optimised CPU, GPU and Hybrid CPU-GPU version. The software was deployed and tested on the Piz Daint cluster of the Swiss National Supercomputing Centre (CSCS). Lenstool-HPC perfectly parallel lens map generation and derivative computation achieves a factor 30 speed-up using only 1 GPUs compared to Lenstool. Lenstool-HPC hybrid Lens-model fit generation tested at Hubble Space Telescope precision is scalable up to 200 CPU-GPU nodes and is faster than Lenstool using only 4 CPU-GPU nodes. △ Less

Submitted 14 April, 2020; originally announced April 2020.

Comments: 9 pages + 4 appendix, published