Showing 1–50 of 317 results for author: Ferreira, P

Euclid preparation. Exploring the properties of proto-clusters in the Simulated Euclid Wide Survey

Authors: Euclid Collaboration, H. Böhringer, G. Chon, O. Cucciati, H. Dannerbauer, M. Bolzonella, G. De Lucia, A. Cappi, L. Moscardini, C. Giocoli, G. Castignani, N. A. Hatch, S. Andreon, E. Bañados, S. Ettori, F. Fontanot, H. Gully, M. Hirschmann, M. Maturi, S. Mei, L. Pozzetti, T. Schlenker, M. Spinelli, N. Aghanim, B. Altieri , et al. (241 additional authors not shown)

Abstract: Galaxy proto-clusters are receiving an increased interest since most of the processes shaping the structure of clusters of galaxies and their galaxy population are happening at early stages of their formation. The Euclid Survey will provide a unique opportunity to discover a large number of proto-clusters over a large fraction of the sky (14 500 square degrees). In this paper, we explore the expec… ▽ More Galaxy proto-clusters are receiving an increased interest since most of the processes shaping the structure of clusters of galaxies and their galaxy population are happening at early stages of their formation. The Euclid Survey will provide a unique opportunity to discover a large number of proto-clusters over a large fraction of the sky (14 500 square degrees). In this paper, we explore the expected observational properties of proto-clusters in the Euclid Wide Survey by means of theoretical models and simulations. We provide an overview of the predicted proto-cluster extent, galaxy density profiles, mass-richness relations, abundance, and sky-filling as a function of redshift. Useful analytical approximations for the functions of these properties are provided. The focus is on the redshift range z= 1.5 to 4. We discuss in particular the density contrast with which proto-clusters can be observed against the background in the galaxy distribution if photometric galaxy redshifts are used as supplied by the ESA Euclid mission together with the ground-based photometric surveys. We show that the obtainable detection significance is sufficient to find large numbers of interesting proto-cluster candidates. For quantitative studies, additional spectroscopic follow-up is required to confirm the proto-clusters and establish their richness. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: Submitted to Astronomy and Astrophysics, 24 pages, 28 figures

Euclid preparation. Forecasting the recovery of galaxy physical properties and their relations with template-fitting and machine-learning methods

Authors: Euclid Collaboration, A. Enia, M. Bolzonella, L. Pozzetti, A. Humphrey, P. A. C. Cunha, W. G. Hartley, F. Dubath, S. Paltani, X. Lopez Lopez, S. Quai, S. Bardelli, L. Bisigello, S. Cavuoti, G. De Lucia, M. Ginolfi, A. Grazian, M. Siudek, C. Tortora, G. Zamorani, N. Aghanim, B. Altieri, A. Amara, S. Andreon, N. Auricchio , et al. (238 additional authors not shown)

Abstract: Euclid will collect an enormous amount of data during the mission's lifetime, observing billions of galaxies in the extragalactic sky. Along with traditional template-fitting methods, numerous Machine Learning algorithms have been presented for computing their photometric redshifts and physical parameters (PP), requiring significantly less computing effort while producing equivalent performance me… ▽ More Euclid will collect an enormous amount of data during the mission's lifetime, observing billions of galaxies in the extragalactic sky. Along with traditional template-fitting methods, numerous Machine Learning algorithms have been presented for computing their photometric redshifts and physical parameters (PP), requiring significantly less computing effort while producing equivalent performance measures. However, their performance is limited by the quality and amount of input information, to the point where the recovery of some well-established physical relationships between parameters might not be guaranteed. To forecast the reliability of Euclid photo-$z$s and PPs calculations, we produced two mock catalogs simulating Euclid photometry. We simulated the Euclid Wide Survey (EWS) and Euclid Deep Fields (EDF). We tested the performance of a template-fitting algorithm (Phosphoros) and four ML methods in recovering photo-$z$s, stellar masses, star-formation rates, and the SFMS. To mimic the Euclid processing as closely as possible, the models were trained with Phosphoros-recovered labels. For the EWS, we found that the best results are achieved with a Mixed Labels approach, training the models with Wide survey features and labels from the Phosphoros results on deeper photometry, i.e., with the best possible set of labels for a given photometry. This imposes a prior, helping the models to better discern cases in degenerate regions of feature space, i.e., when galaxies have similar magnitudes and colors but different redshifts and PPs, with performance metrics even better than those found with Phosphoros. We found no more than $3\%$ performance degradation using a COSMOS-like reference sample or removing $u$ band data, which will not be available until after data release DR1. The best results are obtained for the EDF, with appropriate recovery of photo-$z$, PPs, and the SFMS. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 26 pages, 13 figures. Submitted to A&A

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

Tomographic constraints on the production rate of gravitational waves from astrophysical sources

Authors: David Alonso, Mehraveh Nikjoo, Arianna I. Renzini, Emilio Bellini, Pedro G. Ferreira

Abstract: Using an optimal quadratic estimator, we measure the large-scale cross-correlation between maps of the stochastic gravitational-wave intensity, constructed from the first three LIGO-Virgo observing runs, and a suite of tomographic samples of galaxies covering the redshift range $z\lesssim 2$. We do not detect any statistically significant cross-correlation, but the tomographic nature of the data a… ▽ More Using an optimal quadratic estimator, we measure the large-scale cross-correlation between maps of the stochastic gravitational-wave intensity, constructed from the first three LIGO-Virgo observing runs, and a suite of tomographic samples of galaxies covering the redshift range $z\lesssim 2$. We do not detect any statistically significant cross-correlation, but the tomographic nature of the data allows us to place constraints on the (bias-weighted) production rate density of gravitational waves by astrophysical sources as a function of cosmic time. Our constraints range from $\langle b\dotΩ_{\rm GW}\rangle<3.0\times10^{-9}\,{\rm Gyr}^{-1}$ at $z\sim0.06$ to $\langle b\dotΩ_{\rm GW}\rangle<2.7\times10^{-7}\,{\rm Gyr}^{-1}$ at $z\sim1.5$ (95\% C.L.), assuming a frequency spectrum of the form $f^{2/3}$ (corresponding to an astrophysical background of binary mergers), and a reference frequency $f_{\rm ref}=25\,{\rm Hz}$. Although these constraints are $\sim2$ orders of magnitude higher than the expected signal, we show that a detection may be possible with future experiments. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 14 pages, 11 figures

Euclid preparation. Sensitivity to non-standard particle dark matter model

Authors: Euclid Collaboration, J. Lesgourgues, J. Schwagereit, J. Bucko, G. Parimbelli, S. K. Giri, F. Hervas-Peters, A. Schneider, M. Archidiacono, F. Pace, Z. Sakr, A. Amara, L. Amendola, S. Andreon, N. Auricchio, H. Aussel, C. Baccigalupi, M. Baldi, S. Bardelli, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann , et al. (227 additional authors not shown)

Abstract: The Euclid mission of the European Space Agency will provide weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and its extensions, with an opportunity to test the properties of dark matter beyond the minimal cold dark matter paradigm. We present forecasts from the combination of these surveys on the parameters describing four int… ▽ More The Euclid mission of the European Space Agency will provide weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and its extensions, with an opportunity to test the properties of dark matter beyond the minimal cold dark matter paradigm. We present forecasts from the combination of these surveys on the parameters describing four interesting and representative non-minimal dark matter models: a mixture of cold and warm dark matter relics; unstable dark matter decaying either into massless or massive relics; and dark matter experiencing feeble interactions with relativistic relics. We model these scenarios at the level of the non-linear matter power spectrum using emulators trained on dedicated N-body simulations. We use a mock Euclid likelihood to fit mock data and infer error bars on dark matter parameters marginalised over other parameters. We find that the Euclid photometric probe (alone or in combination with CMB data from the Planck satellite) will be sensitive to the effect of each of the four dark matter models considered here. The improvement will be particularly spectacular for decaying and interacting dark matter models. With Euclid, the bounds on some dark matter parameters can improve by up to two orders of magnitude compared to current limits. We discuss the dependence of predicted uncertainties on different assumptions: inclusion of photometric galaxy clustering data, minimum angular scale taken into account, modelling of baryonic feedback effects. We conclude that the Euclid mission will be able to measure quantities related to the dark sector of particle physics with unprecedented sensitivity. This will provide important information for model building in high-energy physics. Any hint of a deviation from the minimal cold dark matter paradigm would have profound implications for cosmology and particle physics. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 31 pages, 21 figures

Report number: TTK-24-26

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

Euclid. V. The Flagship galaxy mock catalogue: a comprehensive simulation for the Euclid mission

Authors: Euclid Collaboration, F. J. Castander, P. Fosalba, J. Stadel, D. Potter, J. Carretero, P. Tallada-Crespí, L. Pozzetti, M. Bolzonella, G. A. Mamon, L. Blot, K. Hoffmann, M. Huertas-Company, P. Monaco, E. J. Gonzalez, G. De Lucia, C. Scarlata, M. -A. Breton, L. Linke, C. Viglione, S. -S. Li, Z. Zhai, Z. Baghkhani, K. Pardede, C. Neissner , et al. (344 additional authors not shown)

Abstract: We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from… ▽ More We present the Flagship galaxy mock, a simulated catalogue of billions of galaxies designed to support the scientific exploitation of the Euclid mission. Euclid is a medium-class mission of the European Space Agency optimised to determine the properties of dark matter and dark energy on the largest scales of the Universe. It probes structure formation over more than 10 billion years primarily from the combination of weak gravitational lensing and galaxy clustering data. The breath of Euclid's data will also foster a wide variety of scientific analyses. The Flagship simulation was developed to provide a realistic approximation to the galaxies that will be observed by Euclid and used in its scientific analyses. We ran a state-of-the-art N-body simulation with four trillion particles, producing a lightcone on the fly. From the dark matter particles, we produced a catalogue of 16 billion haloes in one octant of the sky in the lightcone up to redshift z=3. We then populated these haloes with mock galaxies using a halo occupation distribution and abundance matching approach, calibrating the free parameters of the galaxy mock against observed correlations and other basic galaxy properties. Modelled galaxy properties include luminosity and flux in several bands, redshifts, positions and velocities, spectral energy distributions, shapes and sizes, stellar masses, star formation rates, metallicities, emission line fluxes, and lensing properties. We selected a final sample of 3.4 billion galaxies with a magnitude cut of H_E<26, where we are complete. We have performed a comprehensive set of validation tests to check the similarity to observational data and theoretical models. In particular, our catalogue is able to closely reproduce the main characteristics of the weak lensing and galaxy clustering samples to be used in the mission's main cosmological analysis. (abridged) △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. IV. The NISP Calibration Unit

Authors: Euclid Collaboration, F. Hormuth, K. Jahnke, M. Schirmer, C. G. -Y. Lee, T. Scott, R. Barbier, S. Ferriol, W. Gillard, F. Grupp, R. Holmes, W. Holmes, B. Kubik, J. Macias-Perez, M. Laurent, J. Marpaud, M. Marton, E. Medinaceli, G. Morgante, R. Toledo-Moreo, M. Trifoglio, Hans-Walter Rix, A. Secroun, M. Seiffert, P. Stassi , et al. (310 additional authors not shown)

Abstract: The near-infrared calibration unit (NI-CU) on board Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid is a mission in ESA's Cosmic Vision 2015-2025 framework, to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and da… ▽ More The near-infrared calibration unit (NI-CU) on board Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid is a mission in ESA's Cosmic Vision 2015-2025 framework, to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and dark energy. Calibrating photometric and spectrometric measurements of galaxies to better than 1.5% accuracy in a survey homogeneously mapping ~14000 deg^2 of extragalactic sky requires a very detailed characterisation of near-infrared (NIR) detector properties, as well their constant monitoring in flight. To cover two of the main contributions - relative pixel-to-pixel sensitivity and non-linearity characteristics - as well as support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous (<12% variations) and temporally stable illumination (0.1%-0.2% over 1200s) over the NISP detector plane, with minimal power consumption and energy dissipation. NI-CU is covers the spectral range ~[900,1900] nm - at cryo-operating temperature - at 5 fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of >=100 from ~15 ph s^-1 pixel^-1 to >1500 ph s^-1 pixel^-1. For this functionality, NI-CU is based on LEDs. We describe the rationale behind the decision and design process, describe the challenges in sourcing the right LEDs, as well as the qualification process and lessons learned. We also provide a description of the completed NI-CU, its capabilities and performance as well as its limits. NI-CU has been integrated into NISP and the Euclid satellite, and since Euclid's launch in July 2023 has started supporting survey operations. △ Less

Submitted 10 July, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: Paper accepted for publication in A&A as part of the special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. III. The NISP Instrument

Authors: Euclid Collaboration, K. Jahnke, W. Gillard, M. Schirmer, A. Ealet, T. Maciaszek, E. Prieto, R. Barbier, C. Bonoli, L. Corcione, S. Dusini, F. Grupp, F. Hormuth, S. Ligori, L. Martin, G. Morgante, C. Padilla, R. Toledo-Moreo, M. Trifoglio, L. Valenziano, R. Bender, F. J. Castander, B. Garilli, P. B. Lilje, H. -W. Rix , et al. (412 additional authors not shown)

Abstract: The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the proc… ▽ More The Near-Infrared Spectrometer and Photometer (NISP) on board the Euclid satellite provides multiband photometry and R>=450 slitless grism spectroscopy in the 950-2020nm wavelength range. In this reference article we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated, and its technical potentials and limitations. Links to articles providing more details and technical background are included. NISP's 16 HAWAII-2RG (H2RG) detectors with a plate scale of 0.3" pix^-1 deliver a field-of-view of 0.57deg^2. In photo mode, NISP reaches a limiting magnitude of ~24.5AB mag in three photometric exposures of about 100s exposure time, for point sources and with a signal-to-noise ratio (SNR) of 5. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux ~2x10^-16erg/s/cm^2 integrated over two resolution elements of 13.4A, in 3x560s grism exposures at 1.6 mu (redshifted Ha). Our calibration includes on-ground and in-flight characterisation and monitoring of detector baseline, dark current, non-linearity, and sensitivity, to guarantee a relative photometric accuracy of better than 1.5%, and relative spectrophotometry to better than 0.7%. The wavelength calibration must be better than 5A. NISP is the state-of-the-art instrument in the NIR for all science beyond small areas available from HST and JWST - and an enormous advance due to its combination of field size and high throughput of telescope and instrument. During Euclid's 6-year survey covering 14000 deg^2 of extragalactic sky, NISP will be the backbone for determining distances of more than a billion galaxies. Its NIR data will become a rich reference imaging and spectroscopy data set for the coming decades. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue 'Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. II. The VIS Instrument

Authors: Euclid Collaboration, M. Cropper, A. Al-Bahlawan, J. Amiaux, S. Awan, R. Azzollini, K. Benson, M. Berthe, J. Boucher, E. Bozzo, C. Brockley-Blatt, G. P. Candini, C. Cara, R. A. Chaudery, R. E. Cole, P. Danto, J. Denniston, A. M. Di Giorgio, B. Dryer, J. Endicott, J. -P. Dubois, M. Farina, E. Galli, L. Genolet, J. P. D. Gow , et al. (403 additional authors not shown)

Abstract: This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift ran… ▽ More This paper presents the specification, design, and development of the Visible Camera (VIS) on the ESA Euclid mission. VIS is a large optical-band imager with a field of view of 0.54 deg^2 sampled at 0.1" with an array of 609 Megapixels and spatial resolution of 0.18". It will be used to survey approximately 14,000 deg^2 of extragalactic sky to measure the distortion of galaxies in the redshift range z=0.1-1.5 resulting from weak gravitational lensing, one of the two principal cosmology probes of Euclid. With photometric redshifts, the distribution of dark matter can be mapped in three dimensions, and, from how this has changed with look-back time, the nature of dark energy and theories of gravity can be constrained. The entire VIS focal plane will be transmitted to provide the largest images of the Universe from space to date, reaching m_AB>24.5 with S/N >10 in a single broad I_E~(r+i+z) band over a six year survey. The particularly challenging aspects of the instrument are the control and calibration of observational biases, which lead to stringent performance requirements and calibration regimes. With its combination of spatial resolution, calibration knowledge, depth, and area covering most of the extra-Galactic sky, VIS will also provide a legacy data set for many other fields. This paper discusses the rationale behind the VIS concept and describes the instrument design and development before reporting the pre-launch performance derived from ground calibrations and brief results from the in-orbit commissioning. VIS should reach fainter than m_AB=25 with S/N>10 for galaxies of full-width half-maximum of 0.3" in a 1.3" diameter aperture over the Wide Survey, and m_AB>26.4 for a Deep Survey that will cover more than 50 deg^2. The paper also describes how VIS works with the other Euclid components of survey, telescope, and science data processing to extract the cosmological information. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Euclid preparation. Sensitivity to neutrino parameters

Authors: Euclid Collaboration, M. Archidiacono, J. Lesgourgues, S. Casas, S. Pamuk, N. Schöneberg, Z. Sakr, G. Parimbelli, A. Schneider, F. Hervas Peters, F. Pace, V. M. Sabarish, M. Costanzi, S. Camera, C. Carbone, S. Clesse, N. Frusciante, A. Fumagalli, P. Monaco, D. Scott, M. Viel, A. Amara, S. Andreon, N. Auricchio, M. Baldi , et al. (224 additional authors not shown)

Abstract: The Euclid mission of the European Space Agency will deliver weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and extensions thereof. We present forecasts from the combination of these surveys on the sensitivity to cosmological parameters including the summed neutrino mass $M_ν$ and the effective number of relativistic species… ▽ More The Euclid mission of the European Space Agency will deliver weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and extensions thereof. We present forecasts from the combination of these surveys on the sensitivity to cosmological parameters including the summed neutrino mass $M_ν$ and the effective number of relativistic species $N_{\rm eff}$ in the standard $Λ$CDM scenario and in a scenario with dynamical dark energy ($w_0 w_a$CDM). We compare the accuracy of different algorithms predicting the nonlinear matter power spectrum for such models. We then validate several pipelines for Fisher matrix and MCMC forecasts, using different theory codes, algorithms for numerical derivatives, and assumptions concerning the non-linear cut-off scale. The Euclid primary probes alone will reach a sensitivity of $σ(M_ν)=$56meV in the $Λ$CDM+$M_ν$ model, whereas the combination with CMB data from Planck is expected to achieve $σ(M_ν)=$23meV and raise the evidence for a non-zero neutrino mass to at least the $2.6σ$ level. This can be pushed to a $4σ$ detection if future CMB data from LiteBIRD and CMB Stage-IV are included. In combination with Planck, Euclid will also deliver tight constraints on $ΔN_{\rm eff}< 0.144$ (95%CL) in the $Λ$CDM+$M_ν$+$N_{\rm eff}$ model, or $ΔN_{\rm eff}< 0.063$ when future CMB data are included. When floating $(w_0, w_a)$, we find that the sensitivity to $N_{\rm eff}$ remains stable, while that to $M_ν$ degrades at most by a factor 2. This work illustrates the complementarity between the Euclid spectroscopic and imaging/photometric surveys and between Euclid and CMB constraints. Euclid will have a great potential for measuring the neutrino mass and excluding well-motivated scenarios with additional relativistic particles. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 48 pages, 33 figures

Report number: TTK-24-17

Euclid preparation. Improving cosmological constraints using a new multi-tracer method with the spectroscopic and photometric samples

Authors: Euclid Collaboration, F. Dournac, A. Blanchard, S. Ilić, B. Lamine, I. Tutusaus, A. Amara, S. Andreon, N. Auricchio, H. Aussel, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, S. Brau-Nogue, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, J. Carretero, S. Casas, M. Castellano, S. Cavuoti, A. Cimatti , et al. (218 additional authors not shown)

Abstract: Future data provided by the Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of Euclid, namely the three-dimensional clustering of the spectroscopic galaxy sample… ▽ More Future data provided by the Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of Euclid, namely the three-dimensional clustering of the spectroscopic galaxy sample, and the so-called 3x2pt signal from the photometric sample (i.e., the weak lensing signal, the galaxy clustering, and their cross-correlation). So far, these two probes have been treated as independent. In this paper, we introduce a new observable given by the ratio of the (angular) two-point correlation function of galaxies from the two surveys. For identical (normalised) selection functions, this observable is unaffected by sampling noise, and its variance is solely controlled by Poisson noise. We present forecasts for Euclid where this multi-tracer method is applied and is particularly relevant because the two surveys will cover the same area of the sky. This method allows for the exploitation of the combination of the spectroscopic and photometric samples. When the correlation between this new observable and the other probes is not taken into account, a significant gain is obtained in the FoM, as well as in the constraints on other cosmological parameters. The benefit is more pronounced for a commonly investigated modified gravity model, namely the $γ$ parametrisation of the growth factor. However, the correlation between the different probes is found to be significant and hence the actual gain is uncertain. We present various strategies for circumventing this issue and still extract useful information from the new observable. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 16 pages, 7 figures, submitted to A&A

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.

Symmetry restoration and vacuum decay from accretion around black holes

Authors: James Marsden, Josu C. Aurrekoetxea, Katy Clough, Pedro G. Ferreira

Abstract: Vacuum decay and symmetry breaking play an important role in the fundamental structure of the matter and the evolution of the universe. In this work we study how the purely classical effect of accretion of fundamental fields onto black holes can lead to shells of symmetry restoration in the midst of a symmetry broken phase. We also show how it can catalyze vacuum decay, forming a bubble that expan… ▽ More Vacuum decay and symmetry breaking play an important role in the fundamental structure of the matter and the evolution of the universe. In this work we study how the purely classical effect of accretion of fundamental fields onto black holes can lead to shells of symmetry restoration in the midst of a symmetry broken phase. We also show how it can catalyze vacuum decay, forming a bubble that expands asymptotically at the speed of light. These effects offer an alternative, purely classical mechanism to quantum tunnelling for seeding phase transitions in the universe. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 4 pages, 4 figures, 1 movie: https://youtu.be/c2vcNoX8fVs Comments welcome!

Tomographic redshift dipole: Testing the cosmological principle

Authors: Pedro da Silveira Ferreira, Valerio Marra

Abstract: The cosmological principle posits that the universe is statistically homogeneous and isotropic on large scales, implying all matter shares the same rest frame. This principle suggests that velocity estimates of our motion from various sources should agree with the cosmic microwave background (CMB) dipole's inferred velocity of 370 km/s. Yet, for over two decades, analyses of radio galaxy and quasa… ▽ More The cosmological principle posits that the universe is statistically homogeneous and isotropic on large scales, implying all matter shares the same rest frame. This principle suggests that velocity estimates of our motion from various sources should agree with the cosmic microwave background (CMB) dipole's inferred velocity of 370 km/s. Yet, for over two decades, analyses of radio galaxy and quasar catalogs have found velocities at odds with the CMB dipole, with tensions up to 5$σ$. In a blind analysis of BOSS and eBOSS spectroscopic data from galaxies and quasars across $0.2<z<2.2$, we applied a novel dipole estimator for a tomographic approach, robustly correcting biases and quantifying uncertainties with realistic mock catalogs. Our findings with eBOSS data ($0.6<z<2.2$), indicating a velocity of $195^{+99}_{-81}$ km/s, demonstrate a $2σ$ agreement with the CMB dipole when considering the full 3D vector distribution and a 3-to-6$σ$ tension with previous number count studies. This result supports the cosmological principle, emphasizing the consistency of our motion with the CMB across vast cosmic distances. On the other hand, the BOSS data revealed potential unmodeled systematics; the estimator could not be minimized using the LOWZ set ($0.2<z<0.4$), and the CMASS set ($0.4<z<0.6$) presented results that pointed towards the southern hemisphere, conflicting with the CMB dipole. Addressing the disparities with earlier number count analyses and understanding possible systematics in spectroscopic measurements will be essential to further validate the cosmological principle. △ Less

Submitted 5 June, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 19 pages, 11 figures; v2: added direct comparison with previous number count studies; v3: revised analysis and results following corrections to the data processing pipeline

syren-halofit: A fast, interpretable, high-precision formula for the $Λ$CDM nonlinear matter power spectrum

Authors: Deaglan J. Bartlett, Benjamin D. Wandelt, Matteo Zennaro, Pedro G. Ferreira, Harry Desmond

Abstract: Rapid and accurate evaluation of the nonlinear matter power spectrum, $P(k)$, as a function of cosmological parameters and redshift is of fundamental importance in cosmology. Analytic approximations provide an interpretable solution, yet current approximations are neither fast nor accurate relative to numerical emulators. We use symbolic regression to obtain simple analytic approximations to the n… ▽ More Rapid and accurate evaluation of the nonlinear matter power spectrum, $P(k)$, as a function of cosmological parameters and redshift is of fundamental importance in cosmology. Analytic approximations provide an interpretable solution, yet current approximations are neither fast nor accurate relative to numerical emulators. We use symbolic regression to obtain simple analytic approximations to the nonlinear scale, $k_σ$, the effective spectral index, $n_{\rm eff}$, and the curvature, $C$, which are required for the halofit model. We then re-optimise the coefficients of halofit to fit a wide range of cosmologies and redshifts. We explore the space of analytic expressions to fit the residuals between $P(k)$ and the optimised predictions of halofit. Our results are designed to match the predictions of EuclidEmulator2, but are validated against $N$-body simulations. Our symbolic expressions for $k_σ$, $n_{\rm eff}$ and $C$ have root mean squared fractional errors of 0.8%, 0.2% and 0.3%, respectively, for redshifts below 3 and a wide range of cosmologies. The re-optimised halofit parameters reduce the root mean squared fractional error (compared to EuclidEmulator2) from 3% to below 2% for wavenumbers $k=9\times10^{-3}-9 \, h{\rm Mpc^{-1}}$. We introduce syren-halofit (symbolic-regression-enhanced halofit), an extension to halofit containing a short symbolic correction which improves this error to 1%. Our method is 2350 and 3170 times faster than current halofit and hmcode implementations, respectively, and 2680 and 64 times faster than EuclidEmulator2 (which requires running class) and the BACCO emulator. We obtain comparable accuracy to EuclidEmulator2 and BACCO when tested on $N$-body simulations. Our work greatly increases the speed and accuracy of symbolic approximations to $P(k)$, making them significantly faster than their numerical counterparts without loss of accuracy. △ Less

Submitted 15 April, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

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

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

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

The VISCACHA survey -- IX. The SMC Southern Bridge in 8D

Authors: M. C. Parisi, R. A. P. Oliveira, M. Angelo, B. Dias, F. Maia, S. Saroon, C. Feinstein, J. F. C. Santos Jr., E. Bica, B. Pereira Lima Ferreira, J. G. Fernández-Trincado, P. Westera, D. Minniti, E. R. Garro, O. J. Katime Santrich, B. De Bortoli, S. Souza, L. Kerber, A. Pérez-Villegas

Abstract: The structure of the Small Magellanic Cloud (SMC) outside of its main body is characterised by tidal branches resulting from its interactions mainly with the Large Magellanic Cloud (LMC). Characterising the stellar populations in these tidal components helps to understand the dynamical history of this galaxy and of the Magellanic system in general. We provide full phase-space vector information fo… ▽ More The structure of the Small Magellanic Cloud (SMC) outside of its main body is characterised by tidal branches resulting from its interactions mainly with the Large Magellanic Cloud (LMC). Characterising the stellar populations in these tidal components helps to understand the dynamical history of this galaxy and of the Magellanic system in general. We provide full phase-space vector information for Southern Bridge clusters. We performed a photometric and spectroscopic analysis of twelve SMC clusters, doubling the number of SMC clusters with full phase-space vector information known to date. We reclassify the sample considering 3D distances and 3D velocities. We found that some of the clusters classified as Southern Bridge objects according to the projected 2D classification actually belong to the Main Body and Counter-Bridge in the background. The comparison of the kinematics of the genuine foreground Bridge clusters with those previously analysed in the same way reveals that Southern Bridge clusters are moving towards the LMC and share the kinematics of the Northern Bridge. Adding to our sample clusters from the literature with CaT metallicity determinations we compare the age-metallicity relation of the Southern Bridge with the one of the Northern Bridge. We reinforce the idea that both regions do not seem to have experienced the same chemical enrichment history and that there is a clear absence of clusters in the Northern Bridge older than 3Gyr and more metal-poor than -1.1, which would not seem to be due to a selection effect. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 18 pages, 9 figures, accepted for publication in MNRAS

Euclid preparation. TBD. Galaxy power spectrum modelling in real space

Authors: Euclid Collaboration, A. Pezzotta, C. Moretti, M. Zennaro, A. Moradinezhad Dizgah, M. Crocce, E. Sefusatti, I. Ferrero, K. Pardede, A. Eggemeier, A. Barreira, R. E. Angulo, M. Marinucci, B. Camacho Quevedo, S. de la Torre, D. Alkhanishvili, M. Biagetti, M. -A. Breton, E. Castorina, G. D'Amico, V. Desjacques, M. Guidi, M. Kärcher, A. Oddo, M. Pellejero Ibanez , et al. (224 additional authors not shown)

Abstract: We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resoluti… ▽ More We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resolution simulations. These models are benchmarked against comoving snapshots of the Flagship I N-body simulation at $z=(0.9,1.2,1.5,1.8)$, which have been populated with H$α$ galaxies leading to catalogues of millions of objects within a volume of about $58\,h^{-3}\,{\rm Gpc}^3$. Our analysis suggests that both models can be used to provide a robust inference of the parameters $(h, ω_{\rm c})$ in the redshift range under consideration, with comparable constraining power. We additionally determine the range of validity of the EFTofLSS model in terms of scale cuts and model degrees of freedom. From these tests, it emerges that the standard third-order Eulerian bias expansion can accurately describe the full shape of the real-space galaxy power spectrum up to the maximum wavenumber $k_{\rm max}=0.45\,h\,{\rm Mpc}^{-1}$, even with a measurement precision well below the percent level. In particular, this is true for a configuration with six free nuisance parameters, including local and non-local bias parameters, a matter counterterm, and a correction to the shot-noise contribution. Fixing either tidal bias parameters to physically-motivated relations still leads to unbiased cosmological constraints. We finally repeat our analysis assuming a volume that matches the expected footprint of Euclid, but without considering observational effects, as purity and completeness, showing that we can get consistent cosmological constraints over this range of scales and redshifts. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 38 pages, 19 figures

Journal ref: A&A 687, A216 (2024)

The effect of wave dark matter on equal mass black hole mergers

Authors: Josu C. Aurrekoetxea, Katy Clough, Jamie Bamber, Pedro G. Ferreira

Abstract: For dark matter to be detectable with gravitational waves from binary black holes, it must reach higher than average densities in their vicinity. In the case of light (wave-like) dark matter, the density of dark matter between the binary can be significantly enhanced by accretion from the surrounding environment. Here we show that the resulting dephasing effect on the last ten orbits of an equal m… ▽ More For dark matter to be detectable with gravitational waves from binary black holes, it must reach higher than average densities in their vicinity. In the case of light (wave-like) dark matter, the density of dark matter between the binary can be significantly enhanced by accretion from the surrounding environment. Here we show that the resulting dephasing effect on the last ten orbits of an equal mass binary is maximized when the Compton wavelength of the scalar particle is comparable to the orbital separation, $2π/μ\sim d$. The phenomenology of the effect is different from the channels that are usually discussed, where dynamical friction (along the orbital path) and radiation of energy and angular momentum drive the dephasing, and is rather dominated by the radial force (the spacetime curvature in the radial direction) towards the overdensity between the black holes. Whilst our numerical studies limit us to scales of the same order, this effect may persist at larger separations and/or particle masses, playing a significant role in the merger history of binaries. △ Less

Submitted 28 May, 2024; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: 4 pages, 4 figures, 1 appendix, 1 movie: https://youtu.be/2VJIfqCp7D8. Matches version published in PRL

A precise symbolic emulator of the linear matter power spectrum

Authors: Deaglan J. Bartlett, Lukas Kammerer, Gabriel Kronberger, Harry Desmond, Pedro G. Ferreira, Benjamin D. Wandelt, Bogdan Burlacu, David Alonso, Matteo Zennaro

Abstract: Computing the matter power spectrum, $P(k)$, as a function of cosmological parameters can be prohibitively slow in cosmological analyses, hence emulating this calculation is desirable. Previous analytic approximations are insufficiently accurate for modern applications, so black-box, uninterpretable emulators are often used. We utilise an efficient genetic programming based symbolic regression fra… ▽ More Computing the matter power spectrum, $P(k)$, as a function of cosmological parameters can be prohibitively slow in cosmological analyses, hence emulating this calculation is desirable. Previous analytic approximations are insufficiently accurate for modern applications, so black-box, uninterpretable emulators are often used. We utilise an efficient genetic programming based symbolic regression framework to explore the space of potential mathematical expressions which can approximate the power spectrum and $σ_8$. We learn the ratio between an existing low-accuracy fitting function for $P(k)$ and that obtained by solving the Boltzmann equations and thus still incorporate the physics which motivated this earlier approximation. We obtain an analytic approximation to the linear power spectrum with a root mean squared fractional error of 0.2% between $k = 9\times10^{-3} - 9 \, h{\rm \, Mpc^{-1}}$ and across a wide range of cosmological parameters, and we provide physical interpretations for various terms in the expression. Our analytic approximation is 950 times faster to evaluate than camb and 36 times faster than the neural network based matter power spectrum emulator BACCO. We also provide a simple analytic approximation for $σ_8$ with a similar accuracy, with a root mean squared fractional error of just 0.1% when evaluated across the same range of cosmologies. This function is easily invertible to obtain $A_{\rm s}$ as a function of $σ_8$ and the other cosmological parameters, if preferred. It is possible to obtain symbolic approximations to a seemingly complex function at a precision required for current and future cosmological analyses without resorting to deep-learning techniques, thus avoiding their black-box nature and large number of parameters. Our emulator will be usable long after the codes on which numerical approximations are built become outdated. △ Less

Submitted 15 April, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 9 pages, 5 figures. Accepted for publication in A&A

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

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)

Euclid preparation. Modelling spectroscopic clustering on mildly nonlinear scales in beyond-$Λ$CDM models

Authors: Euclid Collaboration, B. Bose, P. Carrilho, M. Marinucci, C. Moretti, M. Pietroni, E. Carella, L. Piga, B. S. Wright, F. Vernizzi, C. Carbone, S. Casas, G. D'Amico, N. Frusciante, K. Koyama, F. Pace, A. Pourtsidou, M. Baldi, L. F. de la Bella, B. Fiorini, C. Giocoli, L. Lombriser, N. Aghanim, A. Amara, S. Andreon , et al. (207 additional authors not shown)

Abstract: We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$Λ$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calcu… ▽ More We investigate the approximations needed to efficiently predict the large-scale clustering of matter and dark matter halos in beyond-$Λ$CDM scenarios. We examine the normal branch of the Dvali-Gabadadze-Porrati model, the Hu-Sawicki $f(R)$ model, a slowly evolving dark energy, an interacting dark energy model and massive neutrinos. For each, we test approximations for the perturbative kernel calculations, including the omission of screening terms and the use of perturbative kernels based on the Einstein-de Sitter universe; we explore different infrared-resummation schemes, tracer bias models and a linear treatment of massive neutrinos; we employ two models for redshift space distortions, the Taruya-Nishimishi-Saito prescription and the Effective Field Theory of Large-Scale Structure. This work further provides a preliminary validation of the codes being considered by Euclid for the spectroscopic clustering probe in beyond-$Λ$CDM scenarios. We calculate and compare the $χ^2$ statistic to assess the different modelling choices. This is done by fitting the spectroscopic clustering predictions to measurements from numerical simulations and perturbation theory-based mock data. We compare the behaviour of this statistic in the beyond-$Λ$CDM cases, as a function of the maximum scale included in the fit, to the baseline $Λ$CDM case. We find that the Einstein-de Sitter approximation without screening is surprisingly accurate for all cases when comparing to the halo clustering monopole and quadrupole obtained from simulations. Our results suggest that the inclusion of multiple redshift bins, higher-order multipoles, higher-order clustering statistics (such as the bispectrum) and photometric probes such as weak lensing, will be essential to extract information on massive neutrinos, modified gravity and dark energy. △ Less

Submitted 11 July, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 29 pages, 10 figures, 4 tables, 3 appendices. Journal accepted version

Euclid Preparation. TBD. Impact of magnification on spectroscopic galaxy clustering

Authors: Euclid Collaboration, G. Jelic-Cizmek, F. Sorrenti, F. Lepori, C. Bonvin, S. Camera, F. J. Castander, R. Durrer, P. Fosalba, M. Kunz, L. Lombriser, I. Tutusaus, C. Viglione, Z. Sakr, N. Aghanim, A. Amara, S. Andreon, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, V. Capobianco , et al. (204 additional authors not shown)

Abstract: In this paper we investigate the impact of lensing magnification on the analysis of Euclid's spectroscopic survey, using the multipoles of the 2-point correlation function for galaxy clustering. We determine the impact of lensing magnification on cosmological constraints, and the expected shift in the best-fit parameters if magnification is ignored. We consider two cosmological analyses: i) a full… ▽ More In this paper we investigate the impact of lensing magnification on the analysis of Euclid's spectroscopic survey, using the multipoles of the 2-point correlation function for galaxy clustering. We determine the impact of lensing magnification on cosmological constraints, and the expected shift in the best-fit parameters if magnification is ignored. We consider two cosmological analyses: i) a full-shape analysis based on the $Λ$CDM model and its extension $w_0w_a$CDM and ii) a model-independent analysis that measures the growth rate of structure in each redshift bin. We adopt two complementary approaches in our forecast: the Fisher matrix formalism and the Markov chain Monte Carlo method. The fiducial values of the local count slope (or magnification bias), which regulates the amplitude of the lensing magnification, have been estimated from the Euclid Flagship simulations. We use linear perturbation theory and model the 2-point correlation function with the public code coffe. For a $Λ$CDM model, we find that the estimation of cosmological parameters is biased at the level of 0.4-0.7 standard deviations, while for a $w_0w_a$CDM dynamical dark energy model, lensing magnification has a somewhat smaller impact, with shifts below 0.5 standard deviations. In a model-independent analysis aiming to measure the growth rate of structure, we find that the estimation of the growth rate is biased by up to $1.2$ standard deviations in the highest redshift bin. As a result, lensing magnification cannot be neglected in the spectroscopic survey, especially if we want to determine the growth factor, one of the most promising ways to test general relativity with Euclid. We also find that, by including lensing magnification with a simple template, this shift can be almost entirely eliminated with minimal computational overhead. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 23 pages, 14 figures

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

Optimal Inflationary Potentials

Authors: Tomás Sousa, Deaglan J. Bartlett, Harry Desmond, Pedro G. Ferreira

Abstract: Inflation is a highly favoured theory for the early Universe. It is compatible with current observations of the cosmic microwave background and large scale structure and is a driver in the quest to detect primordial gravitational waves. It is also, given the current quality of the data, highly under-determined with a large number of candidate implementations. We use a new method in symbolic regres… ▽ More Inflation is a highly favoured theory for the early Universe. It is compatible with current observations of the cosmic microwave background and large scale structure and is a driver in the quest to detect primordial gravitational waves. It is also, given the current quality of the data, highly under-determined with a large number of candidate implementations. We use a new method in symbolic regression to generate all possible simple scalar field potentials for one of two possible basis sets of operators. Treating these as single-field, slow-roll inflationary models we then score them with an information-theoretic metric ("minimum description length") that quantifies their efficiency in compressing the information in current data. We explore two possible priors on the parameter space of potentials, one related to the functions' structural complexity and one that uses a Katz back-off language model to prefer functions that may be theoretically motivated. This enables us to identify the inflaton potentials that optimally balance simplicity with accuracy at explaining current data, which may subsequently find theoretical motivation. Our exploratory study opens the door to extraction of fundamental physics directly from data, and may be augmented with more refined theoretical priors in the quest for a complete understanding of the early Universe. △ Less

Submitted 15 April, 2024; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: 13+4 pages, 4 figures; Accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 109, 083524 (2024)

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

LimberJack.jl: auto-differentiable methods for angular power spectra analyses

Authors: J. Ruiz-Zapatero, D. Alonso, C. García-García, A. Nicola, A. Mootoovaloo, J. M. Sullivan, M. Bonici, P. G. Ferreira

Abstract: We present LimberJack.jl, a fully auto-differentiable code for cosmological analyses of 2 point auto- and cross-correlation measurements from galaxy clustering, CMB lensing and weak lensing data written in Julia. Using Julia's auto-differentiation ecosystem, LimberJack.jl can obtain gradients for its outputs up to an order of magnitude faster than traditional finite difference methods. This makes… ▽ More We present LimberJack.jl, a fully auto-differentiable code for cosmological analyses of 2 point auto- and cross-correlation measurements from galaxy clustering, CMB lensing and weak lensing data written in Julia. Using Julia's auto-differentiation ecosystem, LimberJack.jl can obtain gradients for its outputs up to an order of magnitude faster than traditional finite difference methods. This makes LimberJack.jl greatly synergistic with gradient-based sampling methods, such as Hamiltonian Monte Carlo, capable of efficiently exploring parameter spaces with hundreds of dimensions. We first prove LimberJack.jl's reliability by reanalysing the DES Y1 3$\times$2-point data. We then showcase its capabilities by using a O(100) parameters Gaussian Process to reconstruct the cosmic growth from a combination of DES Y1 galaxy clustering and weak lensing data, eBOSS QSO's, CMB lensing and redshift-space distortions. Our Gaussian process reconstruction of the growth factor is statistically consistent with the $Λ$CDM Planck 2018 prediction at all redshifts. Moreover, we show that the addition of RSD data is extremely beneficial to this type of analysis, reducing the uncertainty in the reconstructed growth factor by $20\%$ on average across redshift. LimberJack.jl is a fully open-source project available on Julia's general repository of packages and GitHub. △ Less

Submitted 15 March, 2024; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: Accepted to OJA, corrected bug displaying wrong Fig. 9

Underdetermination of dark energy

Authors: William J. Wolf, Pedro G. Ferreira

Abstract: There is compelling evidence that the Universe is undergoing a late phase of accelerated expansion. One of the simplest explanations for this behaviour is the presence of dark energy. A plethora of microphysical models for dark energy have been proposed. The hope is that, with the ever increasing precision of cosmological surveys, it will be possible to precisely pin down the model. We show that t… ▽ More There is compelling evidence that the Universe is undergoing a late phase of accelerated expansion. One of the simplest explanations for this behaviour is the presence of dark energy. A plethora of microphysical models for dark energy have been proposed. The hope is that, with the ever increasing precision of cosmological surveys, it will be possible to precisely pin down the model. We show that this is unlikely and that, at best, we will have a phenomenological description for the microphysics of dark energy. Furthermore, we argue that the current phenomenological prescriptions are ill-equipped for shedding light on the fundamental theory of dark energy. △ Less

Submitted 15 November, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Updated to match published version and fixed minor typos. 14 pages, 8 figures

Journal ref: Phys. Rev. D 108 (2023) 103519

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, 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. (336 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 26 July, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 18 pages, 8 figures

Angular Correlation Function from sample covariance with BOSS and eBOSS LRG

Authors: Paula S. Ferreira, Ribamar R. R. Reis

Abstract: The Baryon Acoustic Oscillations (BAO) are one of the most used probes to understand the accelerated expansion of the Universe. Traditional methods rely on fiducial model information within their statistical analysis, which may be a problem when constraining different families of models. This work aims to provide a method that constrains $θ_{BAO}$ through a model-independent approach using the cov… ▽ More The Baryon Acoustic Oscillations (BAO) are one of the most used probes to understand the accelerated expansion of the Universe. Traditional methods rely on fiducial model information within their statistical analysis, which may be a problem when constraining different families of models. This work aims to provide a method that constrains $θ_{BAO}$ through a model-independent approach using the covariance matrix from the galaxy sample from thin redshift bins, later validated with a mock sample covariance matrix. We used widths of $δz = 0.002$ separation for all samples as the basis for a sample covariance matrix weighted by the statistical importance of the redshift bin. Each sample belongs to the Sloan Digital Sky Survey: BOSS1, BOSS2, and eBOSS, with effective redshift $z_{eff}$: 0.35, 0.51, 0.71, and different numbers of bins with 50, 100, and 200. To get $θ_{BAO}$, we correct the angular separation from the polynomial fit ($θ_{fit}$) by comparing each bin correlation function with the correlation function of the whole set, a parameter named $\tildeα$. We also tested such correction by choosing the bin at $z_{eff}$ and found that for eBOSS $θ_{BAO}$ is in $1 σ$ agreement with the Planck 18 model. Finally, we found that the sample covariances are noisy compared to the mocks for lower $z$ samples, something expected due to nonlinear effects. Such noise impact can be seen in the parameter constraints but does not affect the eBOSS covariance sample. It is shown that mocks' results do tend to its chosen fiducial cosmology $θ_{BAO}$. BOSS1 and BOSS2 showed agreement with Planck 18 and an agreement with Pantheon + S$H_0$ES when $\tildeα$ is based on the bin $z=z_{eff}$. △ Less

Submitted 6 May, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 18 pages Changed figures and included comparison between old and new method with the mocks

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.

Euclid preparation. XXXIV. The effect of linear redshift-space distortions in photometric galaxy clustering and its cross-correlation with cosmic shear

Authors: Euclid Collaboration, K. Tanidis, V. F. Cardone, M. Martinelli, I. Tutusaus, S. Camera, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, E. Branchini, M. Brescia, J. Brinchmann, V. Capobianco, C. Carbone, J. Carretero, S. Casas, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou, G. Congedo, L. Conversi , et al. (185 additional authors not shown)

Abstract: The cosmological surveys that are planned for the current decade will provide us with unparalleled observations of the distribution of galaxies on cosmic scales, by means of which we can probe the underlying large-scale structure (LSS) of the Universe. This will allow us to test the concordance cosmological model and its extensions. However, precision pushes us to high levels of accuracy in the th… ▽ More The cosmological surveys that are planned for the current decade will provide us with unparalleled observations of the distribution of galaxies on cosmic scales, by means of which we can probe the underlying large-scale structure (LSS) of the Universe. This will allow us to test the concordance cosmological model and its extensions. However, precision pushes us to high levels of accuracy in the theoretical modelling of the LSS observables, so that no biases are introduced into the estimation of the cosmological parameters. In particular, effects such as redshift-space distortions (RSD) can become relevant in the computation of harmonic-space power spectra even for the clustering of the photometrically selected galaxies, as has previously been shown in literature. In this work, we investigate the contribution of linear RSD, as formulated in the Limber approximation by a previous work, in forecast cosmological analyses with the photometric galaxy sample of the Euclid survey. We aim to assess their impact and to quantify the bias on the measurement of cosmological parameters that would be caused if this effect were neglected. We performed this task by producing mock power spectra for photometric galaxy clustering and weak lensing, as is expected to be obtained from the Euclid survey. We then used a Markov chain Monte Carlo approach to obtain the posterior distributions of cosmological parameters from these simulated observations. When the linear RSD is neglected, significant biases are caused when galaxy correlations are used alone and when they are combined with cosmic shear in the so-called 3$\times$2pt approach. These biases can be equivalent to as much as $5\,σ$ when an underlying $Λ$CDM cosmology is assumed. When the cosmological model is extended to include the equation-of-state parameters of dark energy, the extension parameters can be shifted by more than $1\,σ$. △ Less

Submitted 22 April, 2024; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: 15 pages, 6 figures. Version matching publication at journal

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

Euclid: Constraining linearly scale-independent modifications of gravity with the spectroscopic and photometric primary probes

Authors: N. Frusciante, F. Pace, V. F. Cardone, S. Casas, I. Tutusaus, M. Ballardini, E. Bellini, G. Benevento, B. Bose, P. Valageas, N. Bartolo, P. Brax, P. G. Ferreira, F. Finelli, K. Koyama, L. Legrand, L. Lombriser, D. Paoletti, M. Pietroni, A. Rozas-Fernández, Z. Sakr, A. Silvestri, F. Vernizzi, H. A. Winther, N. Aghanim , et al. (105 additional authors not shown)

Abstract: The future Euclid space satellite mission will offer an invaluable opportunity to constrain modifications to Einstein's general relativity at cosmic scales. We focus on modified gravity models characterised, at linear scales, by a scale-independent growth of perturbations while featuring different testable types of derivative screening mechanisms at smaller non-linear scales. We considered three s… ▽ More The future Euclid space satellite mission will offer an invaluable opportunity to constrain modifications to Einstein's general relativity at cosmic scales. We focus on modified gravity models characterised, at linear scales, by a scale-independent growth of perturbations while featuring different testable types of derivative screening mechanisms at smaller non-linear scales. We considered three specific models, namely JBD, a scalar-tensor theory with a flat potential, the nDGP gravity, a braneworld model in which our Universe is a four-dimensional brane embedded in a five-dimensional Minkowski space-time, and $k$-mouflage (KM) gravity, an extension of $k$-essence scenarios with a universal coupling of the scalar field to matter. In preparation for real data, we provide forecasts from spectroscopic and photometric primary probes by Euclid on the cosmological parameters and the additional parameters of the models, respectively, $ω_{\rm BD}$, $Ω_{\rm rc}$ and $ε_{2,0}$. The forecast analysis employs the Fisher matrix method applied to weak lensing (WL); photometric galaxy clustering (GCph), spectroscopic galaxy clustering (GCsp) and the cross-correlation (XC) between GCph and WL. In an optimistic setting at 68.3\% confidence interval, we find the following percentage relative errors with Euclid alone: for $\log_{10}{ω_{\rm BD}}$, with a fiducial value of $ω_{\rm BD}=800$, 27.1\% using GCsp alone, 3.6\% using GCph+WL+XC and 3.2\% using GCph+WL+XC+GCsp; for $\log_{10}{Ω_{\rm rc}}$, with a fiducial value of $Ω_{\rm rc}=0.25$, we find 93.4\%, 20\% and 15\% respectively; and finally, for $ε_{2,0}=-0.04$, we find 3.4\%, 0.15\%, and 0.14\%. (abridged) △ Less

Submitted 11 July, 2024; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 23 pages, 8 figures, 4 tables, 1 appendix, matches Journal version

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

Relativistic drag forces on black holes from scalar dark matter clouds of all sizes

Authors: Dina Traykova, Rodrigo Vicente, Katy Clough, Thomas Helfer, Emanuele Berti, Pedro G. Ferreira, Lam Hui

Abstract: We use numerical simulations of scalar field dark matter evolving on a moving black hole background to confirm the regime of validity of (semi-)analytic expressions derived from first principles for both dynamical friction and momentum accretion in the relativistic regime. We cover both small and large clouds (relative to the de Broglie wavelength of the scalars), and light and heavy particle mass… ▽ More We use numerical simulations of scalar field dark matter evolving on a moving black hole background to confirm the regime of validity of (semi-)analytic expressions derived from first principles for both dynamical friction and momentum accretion in the relativistic regime. We cover both small and large clouds (relative to the de Broglie wavelength of the scalars), and light and heavy particle masses (relative to the BH size). In the case of a small dark matter cloud, the effect of accretion is a non-negligible contribution to the total force on the black hole, even for small scalar masses. We confirm that this momentum accretion transitions between two regimes (wave- and particle-like) and we identify the mass of the scalar at which the transition between regimes occurs. △ Less

Submitted 19 February, 2024; v1 submitted 17 May, 2023; originally announced May 2023.

Comments: 11 pages, 5 figures. Minor corrections and references added to match published version

Journal ref: Phys.Rev.D 108 (2023) 12, L121502

Priors for symbolic regression

Authors: Deaglan J. Bartlett, Harry Desmond, Pedro G. Ferreira

Abstract: When choosing between competing symbolic models for a data set, a human will naturally prefer the "simpler" expression or the one which more closely resembles equations previously seen in a similar context. This suggests a non-uniform prior on functions, which is, however, rarely considered within a symbolic regression (SR) framework. In this paper we develop methods to incorporate detailed prior… ▽ More When choosing between competing symbolic models for a data set, a human will naturally prefer the "simpler" expression or the one which more closely resembles equations previously seen in a similar context. This suggests a non-uniform prior on functions, which is, however, rarely considered within a symbolic regression (SR) framework. In this paper we develop methods to incorporate detailed prior information on both functions and their parameters into SR. Our prior on the structure of a function is based on a $n$-gram language model, which is sensitive to the arrangement of operators relative to one another in addition to the frequency of occurrence of each operator. We also develop a formalism based on the Fractional Bayes Factor to treat numerical parameter priors in such a way that models may be fairly compared though the Bayesian evidence, and explicitly compare Bayesian, Minimum Description Length and heuristic methods for model selection. We demonstrate the performance of our priors relative to literature standards on benchmarks and a real-world dataset from the field of cosmology. △ Less

Submitted 2 June, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 8+2 pages, 2 figures. Accepted for The Genetic and Evolutionary Computation Conference (GECCO) 2023 Workshop on Symbolic Regression

Analytical marginalisation over photometric redshift uncertainties in cosmic shear analyses

Authors: Jaime Ruiz-Zapatero, Boryana Hadzhiyska, David Alonso, Pedro G. Ferreira, Carlos García-García, Arrykrishna Mootoovaloo

Abstract: As the statistical power of imaging surveys grows, it is crucial to account for all systematic uncertainties. This is normally done by constructing a model of these uncertainties and then marginalizing over the additional model parameters. The resulting high dimensionality of the total parameter spaces makes inferring the cosmological parameters significantly more costly using traditional Monte-Ca… ▽ More As the statistical power of imaging surveys grows, it is crucial to account for all systematic uncertainties. This is normally done by constructing a model of these uncertainties and then marginalizing over the additional model parameters. The resulting high dimensionality of the total parameter spaces makes inferring the cosmological parameters significantly more costly using traditional Monte-Carlo sampling methods. A particularly relevant example is the redshift distribution, $p(z)$, of the source samples, which may require tens of parameters to describe fully. However, relatively tight priors can be usually placed on these parameters through calibration of the associated systematics. In this paper we show, quantitatively, that a linearisation of the theoretical prediction with respect to these calibratable systematic parameters allows us to analytically marginalise over these extra parameters, leading to a factor $\sim30$ reduction in the time needed for parameter inference, while accurately recovering the same posterior distributions for the cosmological parameters that would be obtained through a full numerical marginalisation over 160 $p(z)$ parameters. We demonstrate that this is feasible not only with current data and current achievable calibration priors but also for future Stage-IV datasets. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 11 pages, 8 figures, prepared for submission to MNRAS, comments welcome

On the functional form of the radial acceleration relation

Authors: Harry Desmond, Deaglan J. Bartlett, Pedro G. Ferreira

Abstract: We apply a new method for learning equations from data -- Exhaustive Symbolic Regression (ESR) -- to late-type galaxy dynamics as encapsulated in the radial acceleration relation (RAR). Relating the centripetal acceleration due to baryons, $g_\text{bar}$, to the total dynamical acceleration, $g_\text{obs}$, the RAR has been claimed to manifest a new law of nature due to its regularity and tightnes… ▽ More We apply a new method for learning equations from data -- Exhaustive Symbolic Regression (ESR) -- to late-type galaxy dynamics as encapsulated in the radial acceleration relation (RAR). Relating the centripetal acceleration due to baryons, $g_\text{bar}$, to the total dynamical acceleration, $g_\text{obs}$, the RAR has been claimed to manifest a new law of nature due to its regularity and tightness, in agreement with Modified Newtonian Dynamics (MOND). Fits to this relation have been restricted by prior expectations to particular functional forms, while ESR affords an exhaustive and nearly prior-free search through functional parameter space to identify the equations optimally trading accuracy with simplicity. Working with the SPARC data, we find the best functions typically satisfy $g_\text{obs} \propto g_\text{bar}$ at high $g_\text{bar}$, although the coefficient of proportionality is not clearly unity and the deep-MOND limit $g_\text{obs} \propto \sqrt{g_\text{bar}}$ as $g_\text{bar} \to 0$ is little evident at all. By generating mock data according to MOND with or without the external field effect, we find that symbolic regression would not be expected to identify the generating function or reconstruct successfully the asymptotic slopes. We conclude that the limited dynamical range and significant uncertainties of the SPARC RAR preclude a definitive statement of its functional form, and hence that this data alone can neither demonstrate nor rule out law-like gravitational behaviour. △ Less

Submitted 1 March, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 12+4 pages, 4 figures, 3 tables; minor revision to match MNRAS published version

Journal ref: MNRAS 2023 521 (2): 1817-1831

The CMB cold spot under the lens: ruling out a supervoid interpretation

Authors: Stephen Owusu, Pedro da Silveira Ferreira, Alessio Notari, Miguel Quartin

Abstract: The Cosmic Microwave Background (CMB) anisotropies are thought to be statistically isotropic and Gaussian. However, several anomalies are observed, including the CMB Cold Spot, an unexpected cold $\sim 10^{\circ}$ region with $p$-value $\lesssim 0.01$ in standard $Λ$CDM. One of the proposed origins of the Cold Spot is an unusually large void on the line of sight, that would generate a cold region… ▽ More The Cosmic Microwave Background (CMB) anisotropies are thought to be statistically isotropic and Gaussian. However, several anomalies are observed, including the CMB Cold Spot, an unexpected cold $\sim 10^{\circ}$ region with $p$-value $\lesssim 0.01$ in standard $Λ$CDM. One of the proposed origins of the Cold Spot is an unusually large void on the line of sight, that would generate a cold region through the combination of integrated Sachs-Wolfe and Rees-Sciama effects. In the past decade extensive searches were conducted in large scale structure surveys, both in optical and infrared, in the same area for $z \lesssim 1$ and did find evidence of large voids, but of depth and size able to account for only a fraction of the anomaly. Here we analyze the lensing signal in the Planck CMB data and rule out the hypothesis that the Cold Spot could be due to a large void located anywhere between us and the surface of last scattering. In particular, computing the evidence ratio we find that a model with a large void is disfavored compared to $Λ$CDM, with odds 1 : 13 (1 : 20) for SMICA (NILC) maps, compared to the original odds 56 : 1 (21 : 1) using temperature data alone. △ Less

Submitted 14 June, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: v2: typos corrected in some equations; extra analysis in Section 5. Accepted for publication in JCAP

Journal ref: JCAP06(2023)040

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

Euclid preparation. XXVII. Covariance model validation for the 2-point correlation function of galaxy clusters

Authors: Euclid Collaboration, A. Fumagalli, A. Saro, S. Borgani, T. Castro, M. Costanzi, P. Monaco, E. Munari, E. Sefusatti, N. Aghanim, N. Auricchio, M. Baldi, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, R. Cledassou , et al. (169 additional authors not shown)

Abstract: Aims. We validate a semi-analytical model for the covariance of real-space 2-point correlation function of galaxy clusters. Methods. Using 1000 PINOCCHIO light cones mimicking the expected Euclid sample of galaxy clusters, we calibrate a simple model to accurately describe the clustering covariance. Then, we use such a model to quantify the likelihood analysis response to variations of the covaria… ▽ More Aims. We validate a semi-analytical model for the covariance of real-space 2-point correlation function of galaxy clusters. Methods. Using 1000 PINOCCHIO light cones mimicking the expected Euclid sample of galaxy clusters, we calibrate a simple model to accurately describe the clustering covariance. Then, we use such a model to quantify the likelihood analysis response to variations of the covariance, and investigate the impact of a cosmology-dependent matrix at the level of statistics expected for the Euclid survey of galaxy clusters. Results. We find that a Gaussian model with Poissonian shot-noise does not correctly predict the covariance of the 2-point correlation function of galaxy clusters. By introducing few additional parameters fitted from simulations, the proposed model reproduces the numerical covariance with 10 per cent accuracy, with differences of about 5 per cent on the figure of merit of the cosmological parameters $Ω_{\rm m}$ and $σ_8$. Also, we find that the cosmology-dependence of the covariance adds valuable information that is not contained in the mean value, significantly improving the constraining power of cluster clustering. Finally, we find that the cosmological figure of merit can be further improved by taking mass binning into account. Our results have significant implications for the derivation of cosmological constraints from the 2-point clustering statistics of the Euclid survey of galaxy clusters. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures

MSC Class: 85A40

Exhaustive Symbolic Regression

Authors: Deaglan J. Bartlett, Harry Desmond, Pedro G. Ferreira

Abstract: Symbolic Regression (SR) algorithms attempt to learn analytic expressions which fit data accurately and in a highly interpretable manner. Conventional SR suffers from two fundamental issues which we address here. First, these methods search the space stochastically (typically using genetic programming) and hence do not necessarily find the best function. Second, the criteria used to select the equ… ▽ More Symbolic Regression (SR) algorithms attempt to learn analytic expressions which fit data accurately and in a highly interpretable manner. Conventional SR suffers from two fundamental issues which we address here. First, these methods search the space stochastically (typically using genetic programming) and hence do not necessarily find the best function. Second, the criteria used to select the equation optimally balancing accuracy with simplicity have been variable and subjective. To address these issues we introduce Exhaustive Symbolic Regression (ESR), which systematically and efficiently considers all possible equations -- made with a given basis set of operators and up to a specified maximum complexity -- and is therefore guaranteed to find the true optimum (if parameters are perfectly optimised) and a complete function ranking subject to these constraints. We implement the minimum description length principle as a rigorous method for combining these preferences into a single objective. To illustrate the power of ESR we apply it to a catalogue of cosmic chronometers and the Pantheon+ sample of supernovae to learn the Hubble rate as a function of redshift, finding $\sim$40 functions (out of 5.2 million trial functions) that fit the data more economically than the Friedmann equation. These low-redshift data therefore do not uniquely prefer the expansion history of the standard model of cosmology. We make our code and full equation sets publicly available. △ Less

Submitted 29 May, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 15 pages, 7 figures, 2 tables. Accepted for publication in the IEEE Transactions on Evolutionary Computation

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

Combining cosmic shear data with correlated photo-$z$ uncertainties: constraints from DESY1 and HSC-DR1

Authors: Carlos García-García, David Alonso, Pedro G. Ferreira, Boryana Hadzhiyska, Andrina Nicola, Carles Sánchez, Anže Slosar

Abstract: An accurate calibration of the source redshift distribution $p(z)$ is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often u… ▽ More An accurate calibration of the source redshift distribution $p(z)$ is a key aspect in the analysis of cosmic shear data. This, one way or another, requires the use of spectroscopic or high-quality photometric samples. However, the difficulty to obtain colour-complete spectroscopic samples matching the depth of weak lensing catalogs means that the analyses of different cosmic shear datasets often use the same samples for redshift calibration. This introduces a source of statistical and systematic uncertainty that is highly correlated across different weak lensing datasets, and which must be accurately characterised and propagated in order to obtain robust cosmological constraints from their combination. In this paper we introduce a method to quantify and propagate the uncertainties on the source redshift distribution in two different surveys sharing the same calibrating sample. The method is based on an approximate analytical marginalisation of the $p(z)$ statistical uncertainties and the correlated marginalisation of residual systematics. We apply this method to the combined analysis of cosmic shear data from the DESY1 data release and the HSC-DR1 data, using the COSMOS 30-band catalog as a common redshift calibration sample. We find that, although there is significant correlation in the uncertainties on the redshift distributions of both samples, this does not change the final constraints on cosmological parameters significantly. The same is true also for the impact of residual systematic uncertainties from the errors in the COSMOS 30-band photometric redshifts. Additionally, we show that these effects will still be negligible in Stage-IV datasets. Finally, the combination of DESY1 and HSC-DR1 allows us to constrain the ``clumpiness'' parameter to $S_8 = 0.768^{+0.021}_{-0.017}$. This corresponds to a $\sim\sqrt{2}$ improvement in uncertainties with respect to either DES or HSC alone. △ Less

Submitted 21 December, 2022; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 32 pages, 12 figures, comments welcome; v2: minor changes, accepted version in JCAP