Showing 1–50 of 169 results for author: Joachimi, B

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

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

The PAU Survey: Photometric Calibration of Narrow Band Images

Authors: F. J. Castander, S. Serrano, M. Eriksen, E. Gaztanaga, R. Casas, A. Alarcon, A. H. Bauer, E. Fernandez, D. Navarro-Girones, N. Tonello, L. Cabayol, J. Carretero, J. De Vicente, J. Garcia-Bellido, H. Hildebrandt, H. Hoekstra, B. Joachimi, R. Miquel, C. Padilla, P. Renard, E. Sanchez, I. Sevilla-Noarre, P. Tallada-Crespi

Abstract: The Physics of the Accelerating Universe (PAU) camera is an optical narrow band and broad band imaging instrument mounted at the prime focus of the William Herschel Telescope. We describe the image calibration procedure of the PAU Survey data. We rely on an external photometric catalogue to calibrate our narrow band data using stars that have been observed by both datasets. We fit stellar template… ▽ More The Physics of the Accelerating Universe (PAU) camera is an optical narrow band and broad band imaging instrument mounted at the prime focus of the William Herschel Telescope. We describe the image calibration procedure of the PAU Survey data. We rely on an external photometric catalogue to calibrate our narrow band data using stars that have been observed by both datasets. We fit stellar templates to the stellar broad band photometry of the Sloan Digital Sky Survey and synthesise narrow band photometry that we compare to the PAUS narrow band data to determine their calibration. Consequently, the PAUS data are in the AB system as inherited from its reference calibrator. We do several tests to check the performance of the calibration. We find it self-consistent when comparing repeated observations of the same objects, with a good overall accuracy to the AB system which we estimate to be at the 2\% precision level and no significant trends as a function of narrow band filter or wavelength. Repeated observations allow us to build a spatial map of the illumination pattern of the system. We also check the wavelength dependence of the calibration comparing to stellar spectra. We find that using only blue stars reduces the effects of variations in the stellar template fitting to broad-band colours, improving the overall precision of the calibration to around 1\% and its wavelength uniformity. The photometric redshift performance obtained with the PAUS data attests to the validity of our calibration to reach the PAUS science goals. △ Less

Submitted 10 June, 2024; originally announced June 2024.

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

Third-order intrinsic alignment of SDSS BOSS LOWZ galaxies

Authors: Laila Linke, Susan Pyne, Benjamin Joachimi, Christos Georgiou, Kai Hoffmann, Rachel Mandelbaum, Sukhdeep Singh

Abstract: Cosmic shear is a powerful probe of cosmology, but it is affected by the intrinsic alignment (IA) of galaxy shapes with the large-scale structure. Upcoming surveys like Euclid and Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require an accurate understanding of IA, particularly for higher-order cosmic shear statistics that are vital for extracting the most cosmological inform… ▽ More Cosmic shear is a powerful probe of cosmology, but it is affected by the intrinsic alignment (IA) of galaxy shapes with the large-scale structure. Upcoming surveys like Euclid and Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require an accurate understanding of IA, particularly for higher-order cosmic shear statistics that are vital for extracting the most cosmological information. In this paper, we report the first detection of third-order IA correlations using the LOWZ galaxy sample from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS). We compare our measurements with predictions from the MICE cosmological simulation and an analytical NLA-inspired model informed by second-order correlations. We also explore the dependence of the third-order correlation on the galaxies' luminosity. We find that the amplitude $A_\mathrm{IA}$ of the IA signal is non-zero at the $4.7σ$ ($7.6σ$) level for scales between $6 h^{-1} \mathrm{Mpc}$ ($1 h^{-1} \mathrm{Mpc}$) and $20 h^{-1} \mathrm{Mpc}$. For scales above $6 h^{-1}\mathrm{Mpc}$ the inferred AIA agrees both with the prediction from the simulation and estimates from second-order statistics within $1σ$ but deviations arise at smaller scales. Our results demonstrate the feasibility of measuring third-order IA correlations and using them for constraining IA models. The agreement between second- and third-order IA constraints also opens the opportunity for a consistent joint analysis and IA self-calibration, promising tighter parameter constraints for upcoming cosmological surveys. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 9 pages + appendix, 4 figures, submitted to Astronomy & Astrophysics

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. 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'

KiDS-1000 and DES-Y1 combined: Cosmology from peak count statistics

Authors: Joachim Harnois-Deraps, Sven Heydenreich, Benjamin Giblin, Nicolas Martinet, Tilman Troester, Marika Asgari, Pierre Burger, Tiago Castro, Klaus Dolag, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Angus H. Wright

Abstract: We analyse the fourth data release of the Kilo Degree Survey (KiDS-1000) and extract cosmological parameter constraints based on the cosmic shear peak count statistics. Peaks are identified in aperture mass maps in which the filter is maximally sensitive to angular scales in the range 2-4arcmin, probing deep into the non-linear regime of structure formation. We interpret our results with a simulat… ▽ More We analyse the fourth data release of the Kilo Degree Survey (KiDS-1000) and extract cosmological parameter constraints based on the cosmic shear peak count statistics. Peaks are identified in aperture mass maps in which the filter is maximally sensitive to angular scales in the range 2-4arcmin, probing deep into the non-linear regime of structure formation. We interpret our results with a simulation-based inference pipeline, sampling over a broad $w$CDM prior volume and marginalising over uncertainties on shape calibration, photometric redshift distribution, intrinsic alignment and baryonic feedback. Our measurements constrain the structure growth parameter and the amplitude of the non-linear intrinsic alignment model to $Σ_8 \equiv σ_8\left[Ω_{\rm m}/0.3\right]^{0.60}=0.765^{+0.030}_{-0.030}$ and $A_{\rm IA}= 0.71^{+0.42}_{-0.42}$, respectively, in agreement with previous KiDS-1000 results based on two-point shear statistics. These results are robust against modelling of the non-linear physics, different scale cuts and selections of tomographic bins. The posterior is also consistent with that from the Dark Energy Survey Year-1 peak count analysis presented in Harnois-Déraps et al (2021), and hence we jointly analyse both surveys. We obtain $Σ_8^{\rm joint} \equiv σ_8\left[Ω_{\rm m}/0.3\right]^{0.57}=0.759^{+0.020}_{-0.017}$, in agreement with the Planck $w$CDM results. The shear-CMB tension on this parameter increases to $3.1σ$ when forcing $w=-1.0$, and to $4.1σ$ if comparing instead with $S_{8,Λ{\rm CDM}}^{\rm joint} = 0.736^{+0.016}_{-0.018}$, one of the tightest constraints to date on this quantity. (abridged) △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 25 pages, 20 figures. comments welcome

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

KiDS-SBI: Simulation-Based Inference Analysis of KiDS-1000 Cosmic Shear

Authors: Maximilian von Wietersheim-Kramsta, Kiyam Lin, Nicolas Tessore, Benjamin Joachimi, Arthur Loureiro, Robert Reischke, Angus H. Wright

Abstract: We present a simulation-based inference (SBI) cosmological analysis of cosmic shear two-point statistics from the fourth weak gravitational lensing data release of the ESO Kilo-Degree Survey (KiDS-1000). KiDS-SBI efficiently performs non-Limber projection of the matter power spectrum via Levin's method, and constructs log-normal random matter fields on the curved sky for arbitrary cosmologies, inc… ▽ More We present a simulation-based inference (SBI) cosmological analysis of cosmic shear two-point statistics from the fourth weak gravitational lensing data release of the ESO Kilo-Degree Survey (KiDS-1000). KiDS-SBI efficiently performs non-Limber projection of the matter power spectrum via Levin's method, and constructs log-normal random matter fields on the curved sky for arbitrary cosmologies, including effective prescriptions for intrinsic alignments and baryonic feedback. The forward model samples realistic galaxy positions and shapes based on the observational characteristics, incorporating shear measurement and redshift calibration uncertainties, as well as angular anisotropies due to variations in depth and point-spread function. To enable direct comparison with standard inference, we limit our analysis to pseudo-angular power spectra. The SBI is based on sequential neural likelihood estimation to infer the posterior distribution of spatially-flat $Λ$CDM cosmological parameters from 18,000 realisations. We infer a mean marginal of the growth of structure parameter $S_{8} \equiv σ_8 (Ω_\mathrm{m} / 0.3)^{0.5} = 0.731\pm 0.033$ ($68 \%$). We present a measure of goodness-of-fit for SBI and determine that the forward model fits the data well with a probability-to-exceed of $0.42$. For fixed cosmology, the learnt likelihood is approximately Gaussian, while constraints widen compared to a Gaussian likelihood analysis due to cosmology dependence in the covariance. Neglecting variable depth and anisotropies in the point spread function in the model can cause $S_{8}$ to be overestimated by ${\sim}5\%$. Our results are in agreement with previous analysis of KiDS-1000 and reinforce a $2.9 σ$ tension with constraints from cosmic microwave background measurements. This work highlights the importance of forward-modelling systematic effects in upcoming galaxy surveys. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 44 pages, 30 figures. Submitted to Astronomy & Astrophysics

Euclid preparation. XLII. A unified catalogue-level reanalysis of weak lensing by galaxy clusters in five imaging surveys

Authors: Euclid Collaboration, M. Sereno, S. Farrens, L. Ingoglia, G. F. Lesci, L. Baumont, G. Covone, C. Giocoli, F. Marulli, S. Miranda La Hera, M. Vannier, A. Biviano, S. Maurogordato, L. Moscardini, N. Aghanim, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, F. Bellagamba, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann , et al. (199 additional authors not shown)

Abstract: Precise and accurate mass calibration is required to exploit galaxy clusters as astrophysical and cosmological probes in the Euclid era. Systematic errors in lensing signals by galaxy clusters can be empirically estimated by comparing different surveys with independent and uncorrelated systematics. To assess the robustness of the lensing results to systematic errors, we carried out end-to-end test… ▽ More Precise and accurate mass calibration is required to exploit galaxy clusters as astrophysical and cosmological probes in the Euclid era. Systematic errors in lensing signals by galaxy clusters can be empirically estimated by comparing different surveys with independent and uncorrelated systematics. To assess the robustness of the lensing results to systematic errors, we carried out end-to-end tests across different data sets. We performed a unified analysis at the catalogue level by leveraging the Euclid combined cluster and weak-lensing pipeline (COMB-CL). COMB-CL will measure weak lensing cluster masses for the Euclid Survey. Heterogeneous data sets from five independent, recent, lensing surveys (CHFTLenS, DES~SV1, HSC-SSP~S16a, KiDS~DR4, and RCSLenS), which exploited different shear and photometric redshift estimation algorithms, were analysed with a consistent pipeline under the same model assumptions. We performed a comparison of the amplitude of the reduced excess surface density and of the mass estimates using lenses from the Planck PSZ2 and SDSS redMaPPer cluster samples. Mass estimates agree with literature results collected in the LC2 catalogues. Mass accuracy was further investigated considering the AMICO detected clusters in the HSC-SSP XXL North field. The consistency of the data sets was tested using our unified analysis framework. We found agreement between independent surveys, at the level of systematic noise in Stage-III surveys or precursors. This indicates successful control over systematics. If such control continues in Stage-IV, Euclid will be able to measure the weak lensing masses of around 13000 (considering shot noise only) or 3000 (noise from shape and large-scale-structure) massive clusters with a signal-to-noise ratio greater than 3. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 25 pages; in press on A&A

The Wide-field Spectroscopic Telescope (WST) Science White Paper

Authors: Vincenzo Mainieri, Richard I. Anderson, Jarle Brinchmann, Andrea Cimatti, Richard S. Ellis, Vanessa Hill, Jean-Paul Kneib, Anna F. McLeod, Cyrielle Opitom, Martin M. Roth, Paula Sanchez-Saez, Rodolfo Smiljanic, Eline Tolstoy, Roland Bacon, Sofia Randich, Angela Adamo, Francesca Annibali, Patricia Arevalo, Marc Audard, Stefania Barsanti, Giuseppina Battaglia, Amelia M. Bayo Aran, Francesco Belfiore, Michele Bellazzini, Emilio Bellini , et al. (192 additional authors not shown)

Abstract: The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integ… ▽ More The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate △ Less

Submitted 12 April, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

Euclid: Improving the efficiency of weak lensing shear bias calibration. Pixel noise cancellation and the response method on trial

Authors: H. Jansen, M. Tewes, T. Schrabback, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, V. F. Cardone, J. Carretero, S. Casas, M. Castellano, S. Cavuoti, A. Cimatti, G. Congedo, L. Conversi, Y. Copin, L. Corcione, F. Courbin , et al. (96 additional authors not shown)

Abstract: To obtain an accurate cosmological inference from upcoming weak lensing surveys such as the one conducted by Euclid, the shear measurement requires calibration using galaxy image simulations. We study the efficiency of different noise cancellation methods that aim at reducing the simulation volume required to reach a given precision in the shear measurement. Explicitly, we compared fit methods wit… ▽ More To obtain an accurate cosmological inference from upcoming weak lensing surveys such as the one conducted by Euclid, the shear measurement requires calibration using galaxy image simulations. We study the efficiency of different noise cancellation methods that aim at reducing the simulation volume required to reach a given precision in the shear measurement. Explicitly, we compared fit methods with different noise cancellations and a method based on responses. We used GalSim to simulate galaxies both on a grid and at random positions in larger scenes. Placing the galaxies at random positions requires their detection, which we performed with SExtractor. On the grid, we neglected the detection step and, therefore, the potential detection bias arising from it. The shear of the simulated images was measured with the fast moment-based method KSB, for which we note deviations from purely linear shear measurement biases. For the estimation of uncertainties, we used bootstrapping as an empirical method. We find that each method we studied on top of shape noise cancellation can further increase the efficiency of calibration simulations. The improvement depends on the considered shear amplitude range and the type of simulations (grid-based or random positions). The response method on a grid for small shears provides the biggest improvement. In the more realistic case of randomly positioned galaxies, we still find an improvement factor of 70 for small shears using the response method. Alternatively, the runtime can be lowered by a factor of 7 already using pixel noise cancellation on top of shape noise cancellation. Furthermore, we demonstrate that the efficiency of shape noise cancellation can be enhanced in the presence of blending if entire scenes are rotated instead of individual galaxies. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted by A&A. 24 pages, 15+7 figures, Euclid standard paper

Improved Weak Lensing Photometric Redshift Calibration via StratLearn and Hierarchical Modeling

Authors: Maximilian Autenrieth, Angus H. Wright, Roberto Trotta, David A. van Dyk, David C. Stenning, Benjamin Joachimi

Abstract: Discrepancies between cosmological parameter estimates from cosmic shear surveys and from recent Planck cosmic microwave background measurements challenge the ability of the highly successful $Λ$CDM model to describe the nature of the Universe. To rule out systematic biases in cosmic shear survey analyses, accurate redshift calibration within tomographic bins is key. In this paper, we improve phot… ▽ More Discrepancies between cosmological parameter estimates from cosmic shear surveys and from recent Planck cosmic microwave background measurements challenge the ability of the highly successful $Λ$CDM model to describe the nature of the Universe. To rule out systematic biases in cosmic shear survey analyses, accurate redshift calibration within tomographic bins is key. In this paper, we improve photo-$z$ calibration via Bayesian hierarchical modeling of full galaxy photo-$z$ conditional densities, by employing $\textit{StratLearn}$, a recently developed statistical methodology, which accounts for systematic differences in the distribution of the spectroscopic training/source set and the photometric target set. Using realistic simulations that were designed to resemble the KiDS+VIKING-450 dataset, we show that $\textit{StratLearn}$-estimated conditional densities improve the galaxy tomographic bin assignment, and that our $\textit{StratLearn}$-Bayesian framework leads to nearly unbiased estimates of the target population means. This leads to a factor of $\sim 2$ improvement upon the previously best photo-$z$ calibration method. Our approach delivers a maximum bias per tomographic bin of $Δ\langle z \rangle = 0.0095 \pm 0.0089$, with an average absolute bias of $0.0052 \pm 0.0067$ across the five tomographic bins. △ Less

Submitted 12 March, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 24 pages, 20 figures, 3 appendices

Euclid preparation: TBD. The pre-launch Science Ground Segment simulation framework

Authors: Euclid Collaboration, S. Serrano, P. Hudelot, G. Seidel, J. E. Pollack, E. Jullo, F. Torradeflot, D. Benielli, R. Fahed, T. Auphan, J. Carretero, H. Aussel, P. Casenove, F. J. Castander, J. E. Davies, N. Fourmanoit, S. Huot, A. Kara, E. Keihanen, S. Kermiche, K. Okumura, J. Zoubian, A. Ealet, A. Boucaud, H. Bretonniere , et al. (251 additional authors not shown)

Abstract: The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previous… ▽ More The European Space Agency's Euclid mission is one of the upcoming generation of large-scale cosmology surveys, which will map the large-scale structure in the Universe with unprecedented precision. The development and validation of the SGS pipeline requires state-of-the-art simulations with a high level of complexity and accuracy that include subtle instrumental features not accounted for previously as well as faster algorithms for the large-scale production of the expected Euclid data products. In this paper, we present the Euclid SGS simulation framework as applied in a large-scale end-to-end simulation exercise named Science Challenge 8. Our simulation pipeline enables the swift production of detailed image simulations for the construction and validation of the Euclid mission during its qualification phase and will serve as a reference throughout operations. Our end-to-end simulation framework starts with the production of a large cosmological N-body & mock galaxy catalogue simulation. We perform a selection of galaxies down to I_E=26 and 28 mag, respectively, for a Euclid Wide Survey spanning 165 deg^2 and a 1 deg^2 Euclid Deep Survey. We build realistic stellar density catalogues containing Milky Way-like stars down to H<26. Using the latest instrumental models for both the Euclid instruments and spacecraft as well as Euclid-like observing sequences, we emulate with high fidelity Euclid satellite imaging throughout the mission's lifetime. We present the SC8 data set consisting of overlapping visible and near-infrared Euclid Wide Survey and Euclid Deep Survey imaging and low-resolution spectroscopy along with ground-based. This extensive data set enables end-to-end testing of the entire ground segment data reduction and science analysis pipeline as well as the Euclid mission infrastructure, paving the way to future scientific and technical developments and enhancements. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 38 pages, 25 figures, A&A submitted

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

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. Forecast impact of super-sample covariance on 3x2pt analysis with Euclid

Authors: Euclid Collaboration, D. Sciotti, S. Gouyou Beauchamps, V. F. Cardone, S. Camera, I. Tutusaus, F. Lacasa, A. Barreira, A. Gorce, M. Aubert, P. Baratta, R. E. Upham, M. Bonici, C. Carbone, S. Casas, S. Ilić, M. Martinelli, Z. Sakr, A. Schneider, R. Maoli, R. Scaramella, S. Escoffier, W. Gillard, N. Aghanim, A. Amara , et al. (199 additional authors not shown)

Abstract: Deviations from Gaussianity in the distribution of the fields probed by large-scale structure surveys generate additional terms in the data covariance matrix, increasing the uncertainties in the measurement of the cosmological parameters. Super-sample covariance (SSC) is among the largest of these non-Gaussian contributions, with the potential to significantly degrade constraints on some of the pa… ▽ More Deviations from Gaussianity in the distribution of the fields probed by large-scale structure surveys generate additional terms in the data covariance matrix, increasing the uncertainties in the measurement of the cosmological parameters. Super-sample covariance (SSC) is among the largest of these non-Gaussian contributions, with the potential to significantly degrade constraints on some of the parameters of the cosmological model under study -- especially for weak lensing cosmic shear. We compute and validate the impact of SSC on the forecast uncertainties on the cosmological parameters for the Euclid photometric survey, obtained with a Fisher matrix analysis, both considering the Gaussian covariance alone and adding the SSC term -- computed through the public code PySSC. The photometric probes are considered in isolation and combined in the `3$\times$2pt' analysis. We find the SSC impact to be non-negligible -- halving the Figure of Merit of the dark energy parameters ($w_0$, $w_a$) in the 3$\times$2pt case and substantially increasing the uncertainties on $Ω_{{\rm m},0}, w_0$, and $σ_8$ for cosmic shear; photometric galaxy clustering, on the other hand, is less affected due to the lower probe response. The relative impact of SSC does not show significant changes under variations of the redshift binning scheme, while it is smaller for weak lensing when marginalising over the multiplicative shear bias nuisance parameters, which also leads to poorer constraints on the cosmological parameters. Finally, we explore how the use of prior information on the shear and galaxy bias changes the SSC impact. Improving shear bias priors does not have a significant impact, while galaxy bias must be calibrated to sub-percent level to increase the Figure of Merit by the large amount needed to achieve the value when SSC is not included. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 22 pages, 13 figures

Euclid preparation. XXXI. The effect of the variations in photometric passbands on photometric-redshift accuracy

Authors: Euclid Collaboration, Stéphane Paltani, J. Coupon, W. G. Hartley, A. Alvarez-Ayllon, F. Dubath, J. J. Mohr, M. Schirmer, J. -C. Cuillandre, G. Desprez, O. Ilbert, K. Kuijken, N. Aghanim, B. Altieri, A. Amara, N. Auricchio, M. Baldi, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco , et al. (192 additional authors not shown)

Abstract: The technique of photometric redshifts has become essential for the exploitation of multi-band extragalactic surveys. While the requirements on photo-zs for the study of galaxy evolution mostly pertain to the precision and to the fraction of outliers, the most stringent requirement in their use in cosmology is on the accuracy, with a level of bias at the sub-percent level for the Euclid cosmology… ▽ More The technique of photometric redshifts has become essential for the exploitation of multi-band extragalactic surveys. While the requirements on photo-zs for the study of galaxy evolution mostly pertain to the precision and to the fraction of outliers, the most stringent requirement in their use in cosmology is on the accuracy, with a level of bias at the sub-percent level for the Euclid cosmology mission. A separate, and challenging, calibration process is needed to control the bias at this level of accuracy. The bias in photo-zs has several distinct origins that may not always be easily overcome. We identify here one source of bias linked to the spatial or time variability of the passbands used to determine the photometric colours of galaxies. We first quantified the effect as observed on several well-known photometric cameras, and found in particular that, due to the properties of optical filters, the redshifts of off-axis sources are usually overestimated. We show using simple simulations that the detailed and complex changes in the shape can be mostly ignored and that it is sufficient to know the mean wavelength of the passbands of each photometric observation to correct almost exactly for this bias; the key point is that this mean wavelength is independent of the spectral energy distribution of the source}. We use this property to propose a correction that can be computationally efficiently implemented in some photo-z algorithms, in particular template-fitting. We verified that our algorithm, implemented in the new photo-z code Phosphoros, can effectively reduce the bias in photo-zs on real data using the CFHTLS T007 survey, with an average measured bias Delta z over the redshift range 0.4<z<0.7 decreasing by about 0.02, specifically from Delta z~0.04 to Delta z~0.02 around z=0.5. Our algorithm is also able to produce corrected photometry for other applications. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 19 pages, 13 figures; Accepted for publication in A&A

Journal ref: A&A 681, A66 (2024)

Optimal data compression for Lyman-$α$ forest cosmology

Authors: Francesca Gerardi, Andrei Cuceu, Benjamin Joachimi, Seshadri Nadathur, Andreu Font-Ribera

Abstract: The Lyman-$α$ (Ly$α$) three-dimensional correlation functions have been widely used to perform cosmological inference using the baryon acoustic oscillation (BAO) scale. While the traditional inference approach employs a data vector with several thousand data points, we apply near-maximal score compression down to tens of compressed data elements. We show that carefully constructed additional data… ▽ More The Lyman-$α$ (Ly$α$) three-dimensional correlation functions have been widely used to perform cosmological inference using the baryon acoustic oscillation (BAO) scale. While the traditional inference approach employs a data vector with several thousand data points, we apply near-maximal score compression down to tens of compressed data elements. We show that carefully constructed additional data beyond those linked to each inferred model parameter are required to preserve meaningful goodness-of-fit tests that guard against unknown systematics, and to avoid information loss due to non-linear parameter dependencies. We demonstrate, on suites of realistic mocks and DR16 data from the Extended Baryon Oscillation Spectroscopic Survey, that our compression approach is lossless and unbiased, yielding a posterior that is indistinguishable from that of the traditional analysis. As an early application, we investigate the impact of a covariance matrix estimated from a limited number of mocks, which is only well-conditioned in compressed space. △ Less

Submitted 13 February, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: Matches version published in MNRAS

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

Euclid Preparation XXXIII. Characterization of convolutional neural networks for the identification of galaxy-galaxy strong lensing events

Authors: Euclid Collaboration, L. Leuzzi, M. Meneghetti, G. Angora, R. B. Metcalf, L. Moscardini, P. Rosati, P. Bergamini, F. Calura, B. Clément, R. Gavazzi, F. Gentile, M. Lochner, C. Grillo, G. Vernardos, N. Aghanim, A. Amara, L. Amendola, S. Andreon, N. Auricchio, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia , et al. (194 additional authors not shown)

Abstract: Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (… ▽ More Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (CNNs) in particular can efficiently process large volumes of images. We assess and compare the performance of three network architectures in the classification of strong lensing systems on the basis of their morphological characteristics. We train and test our models on different subsamples of a data set of forty thousand mock images, having characteristics similar to those expected in the wide survey planned with the ESA mission \Euclid, gradually including larger fractions of faint lenses. We also evaluate the importance of adding information about the colour difference between the lens and source galaxies by repeating the same training on single-band and multi-band images. Our models find samples of clear lenses with $\gtrsim 90\%$ precision and completeness, without significant differences in the performance of the three architectures. Nevertheless, when including lenses with fainter arcs in the training set, the three models' performance deteriorates with accuracy values of $\sim 0.87$ to $\sim 0.75$ depending on the model. Our analysis confirms the potential of the application of CNNs to the identification of galaxy-scale strong lenses. We suggest that specific training with separate classes of lenses might be needed for detecting the faint lenses since the addition of the colour information does not yield a significant improvement in the current analysis, with the accuracy ranging from $\sim 0.89$ to $\sim 0.78$ for the different models. △ Less

Submitted 26 January, 2024; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 23 pages,12 figures

Journal ref: Astronomy & Astrophysics,2024, 681, A68

KiDS-1000: Cosmology with improved cosmic shear measurements

Authors: Shun-Sheng Li, Henk Hoekstra, Konrad Kuijken, Marika Asgari, Maciej Bilicki, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Lance Miller, Jan Luca van den Busch, Angus H. Wright, Arun Kannawadi, Robert Reischke, HuanYuan Shan

Abstract: We present refined cosmological parameter constraints derived from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000). Our main improvements include enhanced galaxy shape measurements made possible by an updated version of the lensfit code and improved shear calibration achieved with a newly developed suite of multi-band image simulations. Additionally, we inc… ▽ More We present refined cosmological parameter constraints derived from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000). Our main improvements include enhanced galaxy shape measurements made possible by an updated version of the lensfit code and improved shear calibration achieved with a newly developed suite of multi-band image simulations. Additionally, we incorporated recent advancements in cosmological inference from the joint Dark Energy Survey Year 3 and KiDS-1000 cosmic shear analysis. Assuming a spatially flat standard cosmological model, we constrain $S_8\equivσ_8(Ω_{\rm m}/0.3)^{0.5} = 0.776_{-0.027-0.003}^{+0.029+0.002}$, where the second set of uncertainties accounts for the systematic uncertainties within the shear calibration. These systematic uncertainties stem from minor deviations from realism in the image simulations and the sensitivity of the shear measurement algorithm to the morphology of the galaxy sample. Despite these changes, our results align with previous KiDS studies and other weak lensing surveys, and we find a ${\sim}2.3σ$ level of tension with the Planck cosmic microwave background constraints on $S_8$. △ Less

Submitted 2 November, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: 20 pages, 13 figures, 4 tables, minor revisions to match the final accepted version

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

DES Y3 + KiDS-1000: Consistent cosmology combining cosmic shear surveys

Authors: Dark Energy Survey, Kilo-Degree Survey Collaboration, :, T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, M. Asgari, S. Avila, D. Bacon, K. Bechtol, M. R. Becker, G. M. Bernstein, E. Bertin, M. Bilicki, J. Blazek, S. Bocquet, D. Brooks, P. Burger, D. L. Burke, H. Camacho, A. Campos, A. Carnero Rosell , et al. (138 additional authors not shown)

Abstract: We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of… ▽ More We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of $0.790^{+0.018}_{-0.014}$. The mean marginal is lower than the maximum a posteriori estimate, $S_8=0.801$, owing to skewness in the marginal distribution and projection effects in the multi-dimensional parameter space. Our results are consistent with $S_8$ constraints from observations of the cosmic microwave background by Planck, with agreement at the $1.7σ$ level. We use a Hybrid analysis pipeline, defined from a mock survey study quantifying the impact of the different analysis choices originally adopted by each survey team. We review intrinsic alignment models, baryon feedback mitigation strategies, priors, samplers and models of the non-linear matter power spectrum. △ Less

Submitted 19 October, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 21 figures, 15 tables, accepted Open Journal of Astrophysics. Download the chains from https://des.ncsa.illinois.edu/releases/y3a2/Y3key-joint-des-kids or create your own chains with CosmoSIS using https://github.com/joezuntz/cosmosis-standard-library/blob/main/examples/des-y3_and_kids-1000.ini Watch the core team discuss this analysis at https://cosmologytalks.com/2023/05/26/des-kids

Report number: FERMILAB-PUB-23-267-PPD

GLASS: Generator for Large Scale Structure

Authors: Nicolas Tessore, Arthur Loureiro, Benjamin Joachimi, Maximilian von Wietersheim-Kramsta, Niall Jeffrey

Abstract: We present GLASS, the Generator for Large Scale Structure, a new code for the simulation of galaxy surveys for cosmology, which iteratively builds a light cone with matter, galaxies, and weak gravitational lensing signals as a sequence of nested shells. This allows us to create deep and realistic simulations of galaxy surveys at high angular resolution on standard computer hardware and with low re… ▽ More We present GLASS, the Generator for Large Scale Structure, a new code for the simulation of galaxy surveys for cosmology, which iteratively builds a light cone with matter, galaxies, and weak gravitational lensing signals as a sequence of nested shells. This allows us to create deep and realistic simulations of galaxy surveys at high angular resolution on standard computer hardware and with low resource consumption. GLASS also introduces a new technique to generate transformations of Gaussian random fields (including lognormal) to essentially arbitrary precision, an iterative line-of-sight integration over matter shells to obtain weak lensing fields, and flexible modelling of the galaxies sector. We demonstrate that GLASS readily produces simulated data sets with per cent-level accurate two-point statistics of galaxy clustering and weak lensing, thus enabling simulation-based validation and inference that is limited only by our current knowledge of the input matter and galaxy properties. △ Less

Submitted 22 March, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 23 pages, 15 figures; v2 accepted by OJAp with minor changes; code available at https://github.com/glass-dev/glass

Euclid preparation. XXXII. Evaluating the weak lensing cluster mass biases using the Three Hundred Project hydrodynamical simulations

Authors: Euclid Collaboration, C. Giocoli, M. Meneghetti, E. Rasia, S. Borgani, G. Despali, G. F. Lesci, F. Marulli, L. Moscardini, M. Sereno, W. Cui, A. Knebe, G. Yepes, T. Castro, P. -S. Corasaniti, S. Pires, G. Castignani, L. Ingoglia, T. Schrabback, G. W. Pratt, A. M. C. Le Brun, N. Aghanim, L. Amendola, N. Auricchio, M. Baldi , et al. (191 additional authors not shown)

Abstract: The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consiste… ▽ More The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consistent with the implementation within the Euclid Consortium of the dedicated processing function and find that, when jointly modelling mass and the concentration parameter of the Navarro-Frenk-White halo profile, the weak lensing masses tend to be, on average, biased low by 5-10% with respect to the true mass, up to z=0.5. Using a fixed value for the concentration $c_{200} = 3$, the mass bias is diminished below 5%, up to z=0.7, along with its relative uncertainty. Simulating the weak lensing signal by projecting along the directions of the axes of the moment of inertia tensor ellipsoid, we find that orientation matters: when clusters are oriented along the major axis, the lensing signal is boosted, and the recovered weak lensing mass is correspondingly overestimated. Typically, the weak lensing mass bias of individual clusters is modulated by the weak lensing signal-to-noise ratio, related to the redshift evolution of the number of galaxies used for weak lensing measurements: the negative mass bias tends to be larger toward higher redshifts. However, when we use a fixed value of the concentration parameter, the redshift evolution trend is reduced. These results provide a solid basis for the weak-lensing mass calibration required by the cosmological application of future cluster surveys from Euclid and Rubin. △ Less

Submitted 18 October, 2023; v1 submitted 1 February, 2023; originally announced February 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 681, A67 (2024)

KiDS-1000: cross-correlation with Planck cosmic microwave background lensing and intrinsic alignment removal with self-calibration

Authors: Ji Yao, Huanyuan Shan, Pengjie Zhang, Xiangkun Liu, Catherine Heymans, Benjamin Joachimi, Marika Asgari, Maciej Bilicki, Hendrik Hildebrandt, Konrad Kuijken, Tilman Tröster, Jan Luca van den Busch, Angus Wright, Ziang Yan

Abstract: Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared… ▽ More Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared to the standard Planck $Λ$CDM (cosmological constant and cold dark matter) cosmology prediction. In this work, we investigate how IA affects the Kilo-Degree Survey (KiDS) galaxy lensing shear - Planck CMB lensing convergence cross-correlation and compare it to previous treatments with or without IA taken into consideration. More specifically, we compare marginalization over IA parameters and the IA self-calibration (SC) method (with additional observables defined only from the source galaxies) and prove that SC can efficiently break the degeneracy between the CMB lensing amplitude $A_{\rm lens}$ and the IA amplitude $A_{\rm IA}$. We further investigate how different systematics affect the resulting $A_{\rm IA}$ and $A_{\rm lens}$, and validate our results with the MICE2 simulation. We find that by including the SC method to constrain IA, the information loss due to the degeneracy between CMB lensing and IA is strongly reduced. The best-fit values are $A_{\rm lens}=0.84^{+0.22}_{-0.22}$ and $A_{\rm IA}=0.60^{+1.03}_{-1.03}$, while different angular scale cuts can affect $A_{\rm lens}$ by $\sim10\%$. We show that appropriate treatment of the boost factor, cosmic magnification, and photometric redshift modeling is important for obtaining the correct IA and cosmological results. △ Less

Submitted 10 April, 2023; v1 submitted 31 January, 2023; originally announced January 2023.

Comments: match version accepted by A&A

Journal ref: A&A 673, A111 (2023)

Euclid Preparation. XXVIII. Forecasts for ten different higher-order weak lensing statistics

Authors: Euclid Collaboration, V. Ajani, M. Baldi, A. Barthelemy, A. Boyle, P. Burger, V. F. Cardone, S. Cheng, S. Codis, C. Giocoli, J. Harnois-Déraps, S. Heydenreich, V. Kansal, M. Kilbinger, L. Linke, C. Llinares, N. Martinet, C. Parroni, A. Peel, S. Pires, L. Porth, I. Tereno, C. Uhlemann, M. Vicinanza, S. Vinciguerra , et al. (189 additional authors not shown)

Abstract: Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten… ▽ More Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten different HOS on a common set of $Euclid$-like mocks, derived from N-body simulations. In this first paper of the HOWLS series, we computed the nontomographic ($Ω_{\rm m}$, $σ_8$) Fisher information for the one-point probability distribution function, peak counts, Minkowski functionals, Betti numbers, persistent homology Betti numbers and heatmap, and scattering transform coefficients, and we compare them to the shear and convergence two-point correlation functions in the absence of any systematic bias. We also include forecasts for three implementations of higher-order moments, but these cannot be robustly interpreted as the Gaussian likelihood assumption breaks down for these statistics. Taken individually, we find that each HOS outperforms the two-point statistics by a factor of around two in the precision of the forecasts with some variations across statistics and cosmological parameters. When combining all the HOS, this increases to a $4.5$ times improvement, highlighting the immense potential of HOS for cosmic shear cosmological analyses with $Euclid$. The data used in this analysis are publicly released with the paper. △ Less

Submitted 10 July, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 33 pages, 24 figures, main results in Fig. 19 & Table 5, version published in A&A

Journal ref: A&A 675, A120 (2023)

A simulation-based inference pipeline for cosmic shear with the Kilo-Degree Survey

Authors: Kiyam Lin, Maximilian von Wietersheim-Kramsta, Benjamin Joachimi, Stephen Feeney

Abstract: The standard approach to inference from cosmic large-scale structure data employs summary statistics that are compared to analytic models in a Gaussian likelihood with pre-computed covariance. To overcome the idealising assumptions about the form of the likelihood and the complexity of the data inherent to the standard approach, we investigate simulation-based inference (SBI), which learns the lik… ▽ More The standard approach to inference from cosmic large-scale structure data employs summary statistics that are compared to analytic models in a Gaussian likelihood with pre-computed covariance. To overcome the idealising assumptions about the form of the likelihood and the complexity of the data inherent to the standard approach, we investigate simulation-based inference (SBI), which learns the likelihood as a probability density parameterised by a neural network. We construct suites of simulated, exactly Gaussian-distributed data vectors for the most recent Kilo-Degree Survey (KiDS) weak gravitational lensing analysis and demonstrate that SBI recovers the full 12-dimensional KiDS posterior distribution with just under $10^4$ simulations. We optimise the simulation strategy by initially covering the parameter space by a hypercube, followed by batches of actively learnt additional points. The data compression in our SBI implementation is robust to suboptimal choices of fiducial parameter values and of data covariance. Together with a fast simulator, SBI is therefore a competitive and more versatile alternative to standard inference. △ Less

Submitted 23 August, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 14 pages, 8 figures; updated to match version accepted by MNRAS

Report number: MN-22-4878-MJ

Journal ref: MNRAS 524 (2023) 6167-6180

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

KiDS-1000: Combined halo-model cosmology constraints from galaxy abundance, galaxy clustering and galaxy-galaxy lensing

Authors: Andrej Dvornik, Catherine Heymans, Marika Asgari, Constance Mahony, Benjamin Joachimi, Maciej Bilicki, Elisa Chisari, Hendrik Hildebrandt, Henk Hoekstra, Harry Johnston, Konrad Kuijken, Alexander Mead, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke, Sandra Unruh, Angus H. Wright

Abstract: We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matt… ▽ More We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matter field. Our magnitude-limited bright galaxy sample combines 9-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter $S_8=σ_8\sqrt{Ω_{\mathrm{m}}/0.3}=0.773^{+0.028}_{-0.030}$, and the matter density parameter $Ω_{\mathrm{m}}=0.290^{+0.021}_{-0.017}$. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint $3\times2{\mathrm{pt}}$ clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often-excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias and the impact of baryon feedback, this work demonstrates the potential and a way forward to include non-linear scales in cosmological analyses. Varying the width of the satellite galaxy distribution with an additional parameter yields a strong preference for sub-Poissonian variance, improving the goodness of fit by 0.18 in reduced $χ^{2}$ value compared to a fixed Poisson distribution. △ Less

Submitted 13 June, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 25 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 675, A189 (2023)

Constraints on the cosmic expansion rate at redshift 2.3 from the Lyman-$α$ forest

Authors: Andrei Cuceu, Andreu Font-Ribera, Seshadri Nadathur, Benjamin Joachimi, Paul Martini

Abstract: We determine the product of the expansion rate and angular-diameter distance at redshift $z=2.3$ from the anisotropy of Lyman-$α$ (Ly$α$) forest correlations measured by the Sloan Digital Sky Survey (SDSS). Our result is the most precise from large-scale structure at $z>1$. In flat $Λ$CDM we determine the matter density to be $Ω_\mathrm{m}=0.36^{+0.03}_{-0.04}$ from Ly$α$ alone. This is a factor o… ▽ More We determine the product of the expansion rate and angular-diameter distance at redshift $z=2.3$ from the anisotropy of Lyman-$α$ (Ly$α$) forest correlations measured by the Sloan Digital Sky Survey (SDSS). Our result is the most precise from large-scale structure at $z>1$. In flat $Λ$CDM we determine the matter density to be $Ω_\mathrm{m}=0.36^{+0.03}_{-0.04}$ from Ly$α$ alone. This is a factor of two tighter than baryon acoustic oscillation results from the same data due to our use of a wide range of scales ($25<r<180$ $h^{-1}\text{Mpc}$). Using a nucleosynthesis prior, we measure the Hubble constant to be $H_0=63.2\pm2.5$ km/s/Mpc. In combination with other SDSS tracers, we find $H_0=67.2\pm0.9$ km/s/Mpc and measure the dark energy equation-of-state parameter to be $w=-0.90\pm0.12$. Our work opens a new avenue for constraining cosmology at high redshift. △ Less

Submitted 19 June, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

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

Euclid preparation: XXII. Selection of Quiescent Galaxies from Mock Photometry using Machine Learning

Authors: Euclid Collaboration, A. Humphrey, L. Bisigello, P. A. C. Cunha, M. Bolzonella, S. Fotopoulou, K. Caputi, C. Tortora, G. Zamorani, P. Papaderos, D. Vergani, J. Brinchmann, M. Moresco, A. Amara, N. Auricchio, M. Baldi, R. Bender, D. Bonino, E. Branchini, M. Brescia, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander , et al. (184 additional authors not shown)

Abstract: The Euclid Space Telescope will provide deep imaging at optical and near-infrared wavelengths, along with slitless near-infrared spectroscopy, across ~15,000 sq deg of the sky. Euclid is expected to detect ~12 billion astronomical sources, facilitating new insights into cosmology, galaxy evolution, and various other topics. To optimally exploit the expected very large data set, there is the need t… ▽ More The Euclid Space Telescope will provide deep imaging at optical and near-infrared wavelengths, along with slitless near-infrared spectroscopy, across ~15,000 sq deg of the sky. Euclid is expected to detect ~12 billion astronomical sources, facilitating new insights into cosmology, galaxy evolution, and various other topics. To optimally exploit the expected very large data set, there is the need to develop appropriate methods and software. Here we present a novel machine-learning based methodology for selection of quiescent galaxies using broad-band Euclid I_E, Y_E, J_E, H_E photometry, in combination with multiwavelength photometry from other surveys. The ARIADNE pipeline uses meta-learning to fuse decision-tree ensembles, nearest-neighbours, and deep-learning methods into a single classifier that yields significantly higher accuracy than any of the individual learning methods separately. The pipeline has `sparsity-awareness', so that missing photometry values are still informative for the classification. Our pipeline derives photometric redshifts for galaxies selected as quiescent, aided by the `pseudo-labelling' semi-supervised method. After application of the outlier filter, our pipeline achieves a normalized mean absolute deviation of ~< 0.03 and a fraction of catastrophic outliers of ~< 0.02 when measured against the COSMOS2015 photometric redshifts. We apply our classification pipeline to mock galaxy photometry catalogues corresponding to three main scenarios: (i) Euclid Deep Survey with ancillary ugriz, WISE, and radio data; (ii) Euclid Wide Survey with ancillary ugriz, WISE, and radio data; (iii) Euclid Wide Survey only. Our classification pipeline outperforms UVJ selection, in addition to the Euclid I_E-Y_E, J_E-H_E and u-I_E,I_E-J_E colour-colour methods, with improvements in completeness and the F1-score of up to a factor of 2. (Abridged) △ Less

Submitted 5 December, 2022; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: 37 pages (including appendices), 26 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 671, A99 (2023)

The Alcock-Paczyński effect from Lyman-$α$ forest correlations: Analysis validation with synthetic data

Authors: Andrei Cuceu, Andreu Font-Ribera, Paul Martini, Benjamin Joachimi, Seshadri Nadathur, James Rich, Alma X. González-Morales, Hélion du Mas des Bourboux, James Farr

Abstract: The three-dimensional distribution of the Ly$α$ forest has been extensively used to constrain cosmology through measurements of the baryon acoustic oscillations (BAO) scale. However, more cosmological information could be extracted from the full shapes of the Ly$α$ forest correlations through the Alcock-Paczyński (AP) effect. In this work, we prepare for a cosmological analysis of the full shape o… ▽ More The three-dimensional distribution of the Ly$α$ forest has been extensively used to constrain cosmology through measurements of the baryon acoustic oscillations (BAO) scale. However, more cosmological information could be extracted from the full shapes of the Ly$α$ forest correlations through the Alcock-Paczyński (AP) effect. In this work, we prepare for a cosmological analysis of the full shape of the Ly$α$ forest correlations by studying synthetic data of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We use a set of one hundred eBOSS synthetic data sets in order to validate such an analysis. These mocks undergo the same analysis process as the real data. We perform a full-shape analysis on the mean of the correlation functions measured from the one hundred eBOSS realizations, and find that our model of the Ly$α$ correlations performs well on current data sets. We show that we are able to obtain an unbiased full-shape measurement of $D_M/D_H(z_\mathrm{eff})$, where $D_M$ is the transverse comoving distance, $D_H$ is the Hubble distance, and $z_\mathrm{eff}$ is the effective redshift of the measurement. We test the fit over a range of scales, and decide to use a minimum separation of $r_\mathrm{min}=25\ h^{-1}\text{Mpc}$. We also study and discuss the impact of the main contaminants affecting Ly$α$ forest correlations, and give recommendations on how to perform such analysis with real data. While the final eBOSS Ly$α$ BAO analysis measured $D_M/D_H(z_\mathrm{eff}=2.33)$ with $4\%$ statistical precision, a full-shape fit of the same correlations could provide a $\sim2\%$ measurement. △ Less

Submitted 19 June, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Journal ref: MNRAS, Volume 523, Issue 3, August 2023, Pages 3773-3790

Euclid preparation XXVI. The Euclid Morphology Challenge. Towards structural parameters for billions of galaxies

Authors: Euclid Collaboration, H. Bretonnière, U. Kuchner, M. Huertas-Company, E. Merlin, M. Castellano, D. Tuccillo, F. Buitrago, C. J. Conselice, A. Boucaud, B. Häußler, M. Kümmel, W. G. Hartley, A. Alvarez Ayllon, E. Bertin, F. Ferrari, L. Ferreira, R. Gavazzi, D. Hernández-Lang, G. Lucatelli, A. S. G. Robotham, M. Schefer, L. Wang, R. Cabanac, H. Domínguez Sánchez , et al. (193 additional authors not shown)

Abstract: The various Euclid imaging surveys will become a reference for studies of galaxy morphology by delivering imaging over an unprecedented area of 15 000 square degrees with high spatial resolution. In order to understand the capabilities of measuring morphologies from Euclid-detected galaxies and to help implement measurements in the pipeline, we have conducted the Euclid Morphology Challenge, which… ▽ More The various Euclid imaging surveys will become a reference for studies of galaxy morphology by delivering imaging over an unprecedented area of 15 000 square degrees with high spatial resolution. In order to understand the capabilities of measuring morphologies from Euclid-detected galaxies and to help implement measurements in the pipeline, we have conducted the Euclid Morphology Challenge, which we present in two papers. While the companion paper by Merlin et al. focuses on the analysis of photometry, this paper assesses the accuracy of the parametric galaxy morphology measurements in imaging predicted from within the Euclid Wide Survey. We evaluate the performance of five state-of-the-art surface-brightness-fitting codes DeepLeGATo, Galapagos-2, Morfometryka, Profit and SourceXtractor++ on a sample of about 1.5 million simulated galaxies resembling reduced observations with the Euclid VIS and NIR instruments. The simulations include analytic Sérsic profiles with one and two components, as well as more realistic galaxies generated with neural networks. We find that, despite some code-specific differences, all methods tend to achieve reliable structural measurements (10% scatter on ideal Sérsic simulations) down to an apparent magnitude of about 23 in one component and 21 in two components, which correspond to a signal-to-noise ratio of approximately 1 and 5 respectively. We also show that when tested on non-analytic profiles, the results are typically degraded by a factor of 3, driven by systematics. We conclude that the Euclid official Data Releases will deliver robust structural parameters for at least 400 million galaxies in the Euclid Wide Survey by the end of the mission. We find that a key factor for explaining the different behaviour of the codes at the faint end is the set of adopted priors for the various structural parameters. △ Less

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

Comments: Accepted by A&A. 30 pages, 23+6 figures, Euclid pre-launch key paper. Companion paper: Euclid Collaboration XXV: Merlin et al. 2022 Minor corrections after journal review

Journal ref: A&A 671, A102 (2023)

Euclid preparation. XXV. The Euclid Morphology Challenge -- Towards model-fitting photometry for billions of galaxies

Authors: Euclid Collaboration, E. Merlin, M. Castellano, H. Bretonnière, M. Huertas-Company, U. Kuchner, D. Tuccillo, F. Buitrago, J. R. Peterson, C. J. Conselice, F. Caro, P. Dimauro, L. Nemani, A. Fontana, M. Kümmel, B. Häußler, W. G. Hartley, A. Alvarez Ayllon, E. Bertin, P. Dubath, F. Ferrari, L. Ferreira, R. Gavazzi, D. Hernández-Lang, G. Lucatelli , et al. (196 additional authors not shown)

Abstract: The ESA Euclid mission will provide high-quality imaging for about 1.5 billion galaxies. A software pipeline to automatically process and analyse such a huge amount of data in real time is being developed by the Science Ground Segment of the Euclid Consortium; this pipeline will include a model-fitting algorithm, which will provide photometric and morphological estimates of paramount importance fo… ▽ More The ESA Euclid mission will provide high-quality imaging for about 1.5 billion galaxies. A software pipeline to automatically process and analyse such a huge amount of data in real time is being developed by the Science Ground Segment of the Euclid Consortium; this pipeline will include a model-fitting algorithm, which will provide photometric and morphological estimates of paramount importance for the core science goals of the mission and for legacy science. The Euclid Morphology Challenge is a comparative investigation of the performance of five model-fitting software packages on simulated Euclid data, aimed at providing the baseline to identify the best suited algorithm to be implemented in the pipeline. In this paper we describe the simulated data set, and we discuss the photometry results. A companion paper (Euclid Collaboration: Bretonnière et al. 2022) is focused on the structural and morphological estimates. We created mock Euclid images simulating five fields of view of 0.48 deg2 each in the $I_E$ band of the VIS instrument, each with three realisations of galaxy profiles (single and double Sérsic, and 'realistic' profiles obtained with a neural network); for one of the fields in the double Sérsic realisation, we also simulated images for the three near-infrared $Y_E$, $J_E$ and $H_E$ bands of the NISP-P instrument, and five Rubin/LSST optical complementary bands ($u$, $g$, $r$, $i$, and $z$). To analyse the results we created diagnostic plots and defined ad-hoc metrics. Five model-fitting software packages (DeepLeGATo, Galapagos-2, Morfometryka, ProFit, and SourceXtractor++) were compared, all typically providing good results. (cut) △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: 29 pages, 33 figures. Euclid pre-launch key paper. Companion paper: Bretonniere et al. 2022

Journal ref: A&A 671, A101 (2023)

Direct cosmological inference from three-dimensional correlations of the Lyman-$α$ forest

Authors: Francesca Gerardi, Andrei Cuceu, Andreu Font-Ribera, Benjamin Joachimi, Pablo Lemos

Abstract: When performing cosmological inference, standard analyses of the Lyman-$α$ (Ly$α$) three-dimensional correlation functions only consider the information carried by the distinct peak produced by baryon acoustic oscillations (BAO). In this work, we address whether this compression is sufficient to capture all the relevant cosmological information carried by these functions. We do this by performing… ▽ More When performing cosmological inference, standard analyses of the Lyman-$α$ (Ly$α$) three-dimensional correlation functions only consider the information carried by the distinct peak produced by baryon acoustic oscillations (BAO). In this work, we address whether this compression is sufficient to capture all the relevant cosmological information carried by these functions. We do this by performing a direct fit to the full shape, including all physical scales without compression, of synthetic Ly$α$ auto-correlation functions and cross-correlations with quasars at effective redshift $z_{\rm{eff}}=2.3$, assuming a DESI-like survey, and providing a comparison to the classic method applied to the same dataset. Our approach leads to a $3.5\%$ constraint on the matter density $Ω_{\rm{M}}$, which is about three to four times better than what BAO alone can probe. The growth term $f σ_{8} (z_{\rm{eff}})$ is constrained to the $10\%$ level, and the spectral index $n_{\rm{s}}$ to $\sim 3-4\%$. We demonstrate that the extra information resulting from our `direct fit' approach, except for the $n_{\rm{s}}$ constraint, can be traced back to the Alcock-Paczyński effect and redshift space distortion information. △ Less

Submitted 13 December, 2022; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: Matches version published in MNRAS

The PAU Survey & Euclid: Improving broad-band photometric redshifts with multi-task learning

Authors: L. Cabayol, M. Eriksen, J. Carretero, R. Casas, F. J. Castander, E. Fernández, J. Garcia-Bellido, E. Gaztanaga, H. Hildebrandt, H. Hoekstra, B. Joachimi, R. Miquel, C. Padilla, A. Pocino, E. Sanchez, S. Serrano, I. Sevilla, M. Siudek, P. Tallada-Crespí, N. Aghanim, A. Amara, N. Auricchio, M. Baldi, R. Bender, D. Bonino , et al. (101 additional authors not shown)

Abstract: Current and future imaging surveys require photometric redshifts (photo-zs) to be estimated for millions of galaxies. Improving the photo-z quality is a major challenge but is needed to advance our understanding of cosmology. In this paper we explore how the synergies between narrow-band photometric data and large imaging surveys can be exploited to improve broadband photometric redshifts. We used… ▽ More Current and future imaging surveys require photometric redshifts (photo-zs) to be estimated for millions of galaxies. Improving the photo-z quality is a major challenge but is needed to advance our understanding of cosmology. In this paper we explore how the synergies between narrow-band photometric data and large imaging surveys can be exploited to improve broadband photometric redshifts. We used a multi-task learning (MTL) network to improve broadband photo-z estimates by simultaneously predicting the broadband photo-z and the narrow-band photometry from the broadband photometry. The narrow-band photometry is only required in the training field, which also enables better photo-z predictions for the galaxies without narrow-band photometry in the wide field. This technique was tested with data from the Physics of the Accelerating Universe Survey (PAUS) in the COSMOS field. We find that the method predicts photo-zs that are 13% more precise down to magnitude i_{AB} < 23; the outlier rate is also 40% lower when compared to the baseline network. Furthermore, MTL reduces the photo-z bias for high-redshift galaxies, improving the redshift distributions for tomographic bins with z>1. Applying this technique to deeper samples is crucial for future surveys such as \Euclid or LSST. For simulated data, training on a sample with i_{AB} <23, the method reduces the photo-z scatter by 16% for all galaxies with i_{AB}<25. We also studied the effects of extending the training sample with photometric galaxies using PAUS high-precision photo-zs, which reduces the photo-z scatter by 20% in the COSMOS field. △ Less

Submitted 23 January, 2023; v1 submitted 21 September, 2022; originally announced September 2022.

Comments: 20 pages, 16 figures

Journal ref: A&A 671, A153 (2023)

Euclid: Calibrating photometric redshifts with spectroscopic cross-correlations

Authors: K. Naidoo, H. Johnston, B. Joachimi, J. L. van den Busch, H. Hildebrandt, O. Ilbert, O. Lahav, N. Aghanim, B. Altieri, A. Amara, M. Baldi, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou , et al. (104 additional authors not shown)

Abstract: Cosmological constraints from key probes of the Euclid imaging survey rely critically on the accurate determination of the true redshift distributions, $n(z)$, of tomographic redshift bins. We determine whether the mean redshift, $<z>$, of ten Euclid tomographic redshift bins can be calibrated to the Euclid target uncertainties of $σ(<z>)<0.002\,(1+z)$ via cross-correlation, with spectroscopic sam… ▽ More Cosmological constraints from key probes of the Euclid imaging survey rely critically on the accurate determination of the true redshift distributions, $n(z)$, of tomographic redshift bins. We determine whether the mean redshift, $<z>$, of ten Euclid tomographic redshift bins can be calibrated to the Euclid target uncertainties of $σ(<z>)<0.002\,(1+z)$ via cross-correlation, with spectroscopic samples akin to those from the Baryon Oscillation Spectroscopic Survey (BOSS), Dark Energy Spectroscopic Instrument (DESI), and Euclid's NISP spectroscopic survey. We construct mock Euclid and spectroscopic galaxy samples from the Flagship simulation and measure small-scale clustering redshifts up to redshift $z<1.8$ with an algorithm that performs well on current galaxy survey data. The clustering measurements are then fitted to two $n(z)$ models: one is the true $n(z)$ with a free mean; the other a Gaussian Process modified to be restricted to non-negative values. We show that $<z>$ is measured in each tomographic redshift bin to an accuracy of order 0.01 or better. By measuring the clustering redshifts on subsets of the full Flagship area, we construct scaling relations that allow us to extrapolate the method performance to larger sky areas than are currently available in the mock. For the full expected Euclid, BOSS, and DESI overlap region of approximately 6000 deg$^{2}$, the uncertainties attainable by clustering redshifts exceeds the Euclid requirement by at least a factor of three for both $n(z)$ models considered, although systematic biases limit the accuracy. Clustering redshifts are an extremely effective method for redshift calibration for Euclid if the sources of systematic biases can be determined and removed, or calibrated-out with sufficiently realistic simulations. We outline possible future work, in particular an extension to higher redshifts with quasar reference samples. △ Less

Submitted 24 February, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: 14 pages, 8 figures, accepted for publication in A&A

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

KiDS-1000 Cosmology: Constraints from density split statistics

Authors: Pierre A. Burger, Oliver Friedrich, Joachim Harnois-Déraps, Peter Schneider, Marika Asgari, Maciej Bilicki, Hendrik Hildebrandt, Angus H. Wright, Tiago Castro, Klaus Dolag, Catherine Heymans, Benjamin Joachimi, Konrad Kuijken, Nicolas Martinet, HuanYuan Shan, Tilman Tröster

Abstract: Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based o… ▽ More Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = σ_8 \sqrt{Ω_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $Ω_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$. △ Less

Submitted 10 January, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 18 pages, 23 Figures

Report number: A69

Journal ref: A&A 669, A69 (2023)

The PAU Survey: Narrow-band image photometry

Authors: S. Serrano, E. Gaztañaga, F. J. Castander, M. Eriksen, R. Casas, A. Alarcon, A. Bauer, L. Cabayol, J. Carretero, E. Fernandez, D. Navarro-Gironés, C. Neissner, P. Renard, P. Tallada-Crespí, N. Tonello, I. Sevilla-Noarbe, M. Crocce, J. García-Bellido, H. Hildebrandt, H. Hoekstra, B. Joachimi, R. Miquel, C. Padilla, E. Sanchez, J. de Vicente

Abstract: PAUCam is an innovative optical narrow-band imager mounted at the William Herschel Telescope built for the Physics of the Accelerating Universe Survey (PAUS). Its set of 40 filters results in images that are complex to calibrate, with specific instrumental signatures that cannot be processed with traditional data reduction techniques. In this paper we present two pipelines developed by the PAUS da… ▽ More PAUCam is an innovative optical narrow-band imager mounted at the William Herschel Telescope built for the Physics of the Accelerating Universe Survey (PAUS). Its set of 40 filters results in images that are complex to calibrate, with specific instrumental signatures that cannot be processed with traditional data reduction techniques. In this paper we present two pipelines developed by the PAUS data management team with the objective of producing science-ready catalogues from the uncalibrated raw images. The Nightly pipeline takes care of all image processing, with bespoke algorithms for photometric calibration and scatter-light correction. The Multi-Epoch and Multi-Band Analysis (MEMBA) pipeline performs forced photometry over a reference catalogue to optimize the photometric redshift performance. We verify against spectroscopic observations that the current approach delivers an inter-band photometric calibration of 0.8% across the 40 narrow-band set. The large volume of data produced every night and the rapid survey strategy feedback constraints require operating both pipelines in the Port d'Informació Cientifica data centre with intense parallelization. While alternative algorithms for further improvements in photo-z performance are under investigation, the image calibration and photometry presented in this work already enable state-of-the-art photometric redshifts down to iAB=23.0. △ Less

Submitted 27 June, 2023; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 32 pages, 26 figures, MNRAS in press

Optimising the shape of photometric redshift distributions with clustering cross-correlations

Authors: Benjamin Stölzner, Benjamin Joachimi, Andreas Korn, the LSST Dark Energy Science Collaboration

Abstract: We present an optimisation method for the assignment of photometric galaxies into a chosen set of redshift bins. This is achieved by combining simulated annealing, an optimisation algorithm inspired by solid-state physics, with an unsupervised machine learning method, a self-organising map (SOM) of the observed colours of galaxies. Starting with a sample of galaxies that is divided into redshift b… ▽ More We present an optimisation method for the assignment of photometric galaxies into a chosen set of redshift bins. This is achieved by combining simulated annealing, an optimisation algorithm inspired by solid-state physics, with an unsupervised machine learning method, a self-organising map (SOM) of the observed colours of galaxies. Starting with a sample of galaxies that is divided into redshift bins based on a photometric redshift point estimate, the simulated annealing algorithm repeatedly reassigns SOM-selected subsamples of galaxies, which are close in colour, to alternative redshift bins. We optimise the clustering cross-correlation signal between photometric galaxies and a reference sample of galaxies with well-calibrated redshifts. Depending on the effect on the clustering signal, the reassignment is either accepted or rejected. By dynamically increasing the resolution of the SOM, the algorithm eventually converges to a solution that minimises the number of mismatched galaxies in each tomographic redshift bin and thus improves the compactness of their corresponding redshift distribution. This method is demonstrated on the synthetic LSST cosmoDC2 catalogue. We find a significant decrease in the fraction of catastrophic outliers in the redshift distribution in all tomographic bins, most notably in the highest redshift bin with a decrease in the outlier fraction from 57 per cent to 16 per cent. △ Less

Submitted 7 February, 2023; v1 submitted 26 May, 2022; originally announced May 2022.

Comments: 13 pages, 9 figures, version accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 519, Issue 2, February 2023, Pages 2438-2450

Fast and realistic large-scale structure from machine-learning-augmented random field simulations

Authors: Davide Piras, Benjamin Joachimi, Francisco Villaescusa-Navarro

Abstract: Producing thousands of simulations of the dark matter distribution in the Universe with increasing precision is a challenging but critical task to facilitate the exploitation of current and forthcoming cosmological surveys. Many inexpensive substitutes to full $N$-body simulations have been proposed, even though they often fail to reproduce the statistics of the smaller, non-linear scales. Among t… ▽ More Producing thousands of simulations of the dark matter distribution in the Universe with increasing precision is a challenging but critical task to facilitate the exploitation of current and forthcoming cosmological surveys. Many inexpensive substitutes to full $N$-body simulations have been proposed, even though they often fail to reproduce the statistics of the smaller, non-linear scales. Among these alternatives, a common approximation is represented by the lognormal distribution, which comes with its own limitations as well, while being extremely fast to compute even for high-resolution density fields. In this work, we train a generative deep learning model, mainly made of convolutional layers, to transform projected lognormal dark matter density fields to more realistic dark matter maps, as obtained from full $N$-body simulations. We detail the procedure that we follow to generate highly correlated pairs of lognormal and simulated maps, which we use as our training data, exploiting the information of the Fourier phases. We demonstrate the performance of our model comparing various statistical tests with different field resolutions, redshifts and cosmological parameters, proving its robustness and explaining its current limitations. When evaluated on 100 test maps, the augmented lognormal random fields reproduce the power spectrum up to wavenumbers of $1 \ h \ \rm{Mpc}^{-1}$, and the bispectrum within 10%, and always within the error bars, of the fiducial target simulations. Finally, we describe how we plan to integrate our proposed model with existing tools to yield more accurate spherical random fields for weak lensing analysis. △ Less

Submitted 1 February, 2023; v1 submitted 16 May, 2022; originally announced May 2022.

Comments: 16 pages, 11 figures. Matches MNRAS published version, which includes more tests with e.g. varying cosmological parameters

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 520, Issue 1, March 2023, Pages 668-683

Three-point intrinsic alignments of dark matter halos in the IllustrisTNG simulation

Authors: Susan Pyne, Ananth Tenneti, Benjamin Joachimi

Abstract: We use the IllustrisTNG suite of cosmological simulations to measure intrinsic alignment (IA) bispectra of dark matter subhalos between redshifts 0 and 1. We decompose the intrinsic shear field into E- and B-modes and find that the bispectra $B_{δδ\mathrm{E}}$ and $B_{δ\mathrm{EE}}$, between the matter overdensity field, $δ$, and the E-mode field, are detected with high significance. We also model… ▽ More We use the IllustrisTNG suite of cosmological simulations to measure intrinsic alignment (IA) bispectra of dark matter subhalos between redshifts 0 and 1. We decompose the intrinsic shear field into E- and B-modes and find that the bispectra $B_{δδ\mathrm{E}}$ and $B_{δ\mathrm{EE}}$, between the matter overdensity field, $δ$, and the E-mode field, are detected with high significance. We also model the IA bispectra analytically using a method consistent with the two-point non-linear alignment model. We use this model and the simulation measurements to infer the intrinsic alignment amplitude $A_\mathrm{IA}$ and find that values of $A_\mathrm{IA}$ obtained from IA power spectra and bispectra agree well at scales up to $k_\mathrm{max}= 2 \, h \mathrm{Mpc}^{-1}$, for example at $z=1$ $A_\mathrm{IA} = 2.13 \pm$ 0.02 from the cross power spectrum between the matter overdensity and E-mode fields and $A_\mathrm{IA} =2.11 \pm$ 0.03 from $B_{δδ\mathrm{E}}$. This demonstrates that a single physically motivated model can jointly model two-point and three-point statistics of intrinsic alignments, thus enabling a cleaner separation between intrinsic alignments and cosmological weak lensing signals. △ Less

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

Comments: Consistent with version published in MNRAS

Consistent lensing and clustering in a low-$S_8$ Universe with BOSS, DES Year 3, HSC Year 1 and KiDS-1000

Authors: A. Amon, N. C. Robertson, H. Miyatake, C. Heymans, M. White, J. DeRose, S. Yuan, R. H. Wechsler, T. N. Varga, S. Bocquet, A. Dvornik, S. More, A. J. Ross, H. Hoekstra, A. Alarcon, M. Asgari, J. Blazek, A. Campos, R. Chen, A. Choi, M. Crocce, H. T. Diehl, C. Doux, K. Eckert, J. Elvin-Poole , et al. (83 additional authors not shown)

Abstract: We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data a… ▽ More We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck, with $S_8=0.83$, and a Lensing cosmology with $S_8=0.76$. For a joint analysis limited to scales with $R>5.25h^{-1}$Mpc, we find that both cosmologies provide an acceptable fit to the data. Full utilisation of the small-scale clustering and lensing measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic scaling parameter for each redshift bin, $A$, that decouples the lensing and clustering to capture any inconsistency. When a wide range of scales ($0.15<R<60h^{-1}$Mpc) are incorporated, we find different results for the consistency of clustering and lensing between the two cosmologies. Limiting the analysis to the bins for which the impact of the selection of the lens sample is expected to be minimal, for the low-$S_8$ Lensing cosmology, the measurements are consistent with $A$=1; $A=0.91\pm0.04$ using DES+KiDS and $A=0.97\pm0.06$ using HSC. For the Planck cosmology case, we find a discrepancy: $A=0.79\pm0.03$ using DES+KiDS and $A=0.84\pm0.05$ using HSC. We demonstrate that a kSZ-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of these small-scale measurements, but also indicates that caution is still warranted given current uncertainties in modelling baryonic effects, assembly bias, and selection effects in the foreground sample. △ Less

Submitted 13 October, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 28 pages, 11 figures

The PAU Survey: Measurements of the 4000 Å spectral break with narrow-band photometry

Authors: Pablo Renard, Małgorzata Siudek, Martin B. Eriksen, Laura Cabayol, Zheng Cai, Jorge Carretero, Ricard Casas, Francisco J. Castander, Enrique Fernandez, Juan García-Bellido, Enrique Gaztanaga, Henk Hoekstra, Benjamin Joachimi, Ramon Miquel, David Navarro-Girones, Cristóbal Padilla, Eusebio Sanchez, Santiago Serrano, Pau Tallada-Crespí, Juan De Vicente, Anna Wittje, Angus H. Wright

Abstract: The D4000 spectral break index is one of the most important features in the visible spectrum, as it is a proxy for stellar ages and is also used in galaxy classification. However, its direct measurement has always been reserved to spectroscopy. Here, we present a general method to directly measure the D4000 with narrow-band (NB) photometry; it has been validated using realistic simulations, and th… ▽ More The D4000 spectral break index is one of the most important features in the visible spectrum, as it is a proxy for stellar ages and is also used in galaxy classification. However, its direct measurement has always been reserved to spectroscopy. Here, we present a general method to directly measure the D4000 with narrow-band (NB) photometry; it has been validated using realistic simulations, and then evaluated with PAUS NBs, cross-matched with VIPERS spectra ($i_{\rm AB} < 22.5$, $0.562 < z < 0.967$). We also reconstruct the D4000 with the SED-fitting code CIGALE; the use of PAUS NBs instead of broad bands significantly improves the SED fitting results. For D4000$_{\rm n}$, the direct measurement has $\rm \langle SNR \rangle \sim 4$, but we find that for $i_{\rm AB}<21$ all direct D4000 measurements have $\rm SNR>3$. The CIGALE D4000$_{\rm n}$ has $\rm \langle SNR \rangle \sim 20$, but underestimates the error by $>$50\%. Furthermore, the direct method recreates well the D4000-SFR relation, as well as the D4000-mass relation for blue galaxies (for red galaxies, selection effects impact the results). On the other hand, CIGALE accurately classifies galaxies into red and blue populations. We conclude that the direct measurement of D4000 with narrow-band photometry is a promising tool to determine average properties of galaxy samples, with results compatible with spectroscopy. △ Less

Submitted 22 July, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 21 pages, 18 figures. Accepted in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 515, Issue 1, September 2022, Pages 146-166

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

Euclid: Covariance of weak lensing pseudo-$C_\ell$ estimates. Calculation, comparison to simulations, and dependence on survey geometry

Authors: R. E. Upham, M. L. Brown, L. Whittaker, A. Amara, N. Auricchio, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, V. Capobianco, C. Carbone, J. Carretero, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou, G. Congedo, L. Conversi, Y. Copin, L. Corcione, M. Cropper, A. Da Silva, H. Degaudenzi, M. Douspis, F. Dubath , et al. (80 additional authors not shown)

Abstract: An accurate covariance matrix is essential for obtaining reliable cosmological results when using a Gaussian likelihood. In this paper we study the covariance of pseudo-$C_\ell$ estimates of tomographic cosmic shear power spectra. Using two existing publicly available codes in combination, we calculate the full covariance matrix, including mode-coupling contributions arising from both partial sky… ▽ More An accurate covariance matrix is essential for obtaining reliable cosmological results when using a Gaussian likelihood. In this paper we study the covariance of pseudo-$C_\ell$ estimates of tomographic cosmic shear power spectra. Using two existing publicly available codes in combination, we calculate the full covariance matrix, including mode-coupling contributions arising from both partial sky coverage and non-linear structure growth. For three different sky masks, we compare the theoretical covariance matrix to that estimated from publicly available N-body weak lensing simulations, finding good agreement. We find that as a more extreme sky cut is applied, a corresponding increase in both Gaussian off-diagonal covariance and non-Gaussian super-sample covariance is observed in both theory and simulations, in accordance with expectations. Studying the different contributions to the covariance in detail, we find that the Gaussian covariance dominates along the main diagonal and the closest off-diagonals, but further away from the main diagonal the super-sample covariance is dominant. Forming mock constraints in parameters describing matter clustering and dark energy, we find that neglecting non-Gaussian contributions to the covariance can lead to underestimating the true size of confidence regions by up to 70 per cent. The dominant non-Gaussian covariance component is the super-sample covariance, but neglecting the smaller connected non-Gaussian covariance can still lead to the underestimation of uncertainties by 10--20 per cent. A real cosmological analysis will require marginalisation over many nuisance parameters, which will decrease the relative importance of all cosmological contributions to the covariance, so these values should be taken as upper limits on the importance of each component. △ Less

Submitted 17 February, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: 15 pages, 8 figures; matches version accepted by A&A; code available at https://github.com/robinupham/shear_pcl_cov

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

Forecasting the potential of weak lensing magnification to enhance LSST large-scale structure analyses

Authors: Constance Mahony, Maria Cristina Fortuna, Benjamin Joachimi, Andreas Korn, Henk Hoekstra, Samuel J. Schmidt, The LSST Dark Energy Science Collaboration

Abstract: Recent works have shown that weak lensing magnification must be included in upcoming large-scale structure analyses, such as for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), to avoid biasing the cosmological results. In this work we investigate whether including magnification has a positive impact on the precision of the cosmological constraints, as well as being necessary… ▽ More Recent works have shown that weak lensing magnification must be included in upcoming large-scale structure analyses, such as for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), to avoid biasing the cosmological results. In this work we investigate whether including magnification has a positive impact on the precision of the cosmological constraints, as well as being necessary to avoid bias. We forecast this using an LSST mock catalog and a halo model to calculate the galaxy power spectra. We find that including magnification has little effect on the precision of the cosmological parameter constraints for an LSST galaxy clustering analysis, where the halo model parameters are additionally constrained by the galaxy luminosity function. In particular, we find that for the LSST gold sample ($i < 25.3$) including weak lensing magnification only improves the galaxy clustering constraint on $Ω_{\rm{m}}$ by a factor of 1.03, and when using a very deep LSST mock sample ($i<26.5$) by a factor of 1.3. Since magnification predominantly contributes to the clustering measurement and provides similar information to that of cosmic shear, this improvement would be reduced for a combined galaxy clustering and shear analysis. We also confirm that not modelling weak lensing magnification will catastrophically bias the cosmological results from LSST. Magnification must therefore be included in LSST large-scale structure analyses even though it does not significantly enhance the precision of the cosmological constraints. △ Less

Submitted 28 April, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: v2: Version accepted in MNRAS. Minor clarifications, conclusions unchanged