Showing 1–50 of 135 results for author: Joudaki, S

Euclid preparation. First investigation of the impact of cross-contamination on spectroscopic redshift measurements with pixel-level simulations

Authors: Euclid Collaboration, F. Passalacqua, S. Anselmi, P. Monaco, C. Sirignano, S. Dusini, N. Fourmanoit, M. Fumana, E. Lecrivain, K. S. McCarthy, M. Moresco, F. Oppizzi, A. Renzi, M. Scodeggio, L. Stanco, A. Troja, S. Bruton, C. Carbone, S. de la Torre, B. R. Granett, G. Lavaux, S. Lee, K. Markovic, W. J. Percival, I. Risso , et al. (277 additional authors not shown)

Abstract: We present a study on simulated data focused on understanding the performance of the spectroscopic redshift measurements with the Near-Infrared Spectrometer and Photometer (NISP) instrument on Euclid. Simulations include scenarios with different levels of cross-contamination arising from overlapping spectra of nearby sources, which represents one of the main drawbacks of slitless spectroscopy. We… ▽ More We present a study on simulated data focused on understanding the performance of the spectroscopic redshift measurements with the Near-Infrared Spectrometer and Photometer (NISP) instrument on Euclid. Simulations include scenarios with different levels of cross-contamination arising from overlapping spectra of nearby sources, which represents one of the main drawbacks of slitless spectroscopy. We present a new analysis based on pixel-level simulations of the NISP images, with the data processed using the Euclid spectroscopic pipeline. We first consider an idealised case with non-overlapping spectra to assess the accuracy and reliability of the redshift measurement as a function of the flux of the H-alpha emission line and galaxy size. We then introduce more realistic contamination scenarios, distinguishing between two contributions: contamination from H-alpha emitters, which are the Euclid targets for cosmological analyses, and contamination from all other galaxies. In the second case, we analyse the impact of cross-contamination with an increasing number of contaminants, from the brighter to the fainter galaxies. Given that our results show no clear evidence that sources fainter than magnitude 20 degrade redshift measurements, we conservatively restrict our analysis to galaxies with magnitudes up to 24. In particular, we provide a preliminary estimate that contamination from galaxies within the same redshift range as the target sample contributes to about 4% of the total degradation due to cross-contamination from all galaxies. △ Less

Submitted 23 June, 2026; originally announced June 2026.

Comments: 18 pages, 15 pages

KiDS-Legacy: Joint analysis of second- and third-order cosmic shear

Authors: L. Linke, L. Porth, P. Burger, J. Harnois-Déraps, S. Heydenreich, P. Schneider, M. Asgari, M. Bilicki, C. Georgiou, C. Heymans, H. Hildebrandt, H. Hoekstra, P. Jalan, B. Joachimi, S. Joudaki, K. Kuijken, S. Li, L. Moscardini, M. Radovich, R. Reischke, B. Stölzner, A. H. Wright, Z. Yan, Y. -H. Zhang

Abstract: Weak lensing by large-scale structure is a powerful cosmological probe. While most analyses rely on second-order correlations, these are primarily sensitive to the parameter combination $S_8 = σ_8 (Ω_m/0.3)^{0.5}$, limiting their ability to constrain $Ω_m$ and other cosmological parameters independently. Higher-order statistics capture non-Gaussian features of the density field and can therefore b… ▽ More Weak lensing by large-scale structure is a powerful cosmological probe. While most analyses rely on second-order correlations, these are primarily sensitive to the parameter combination $S_8 = σ_8 (Ω_m/0.3)^{0.5}$, limiting their ability to constrain $Ω_m$ and other cosmological parameters independently. Higher-order statistics capture non-Gaussian features of the density field and can therefore break parameter degeneracies and extract more cosmological information from weak lensing surveys. We present a joint analysis of second- and third-order cosmic shear in the final data release of the Kilo-Degree Survey (KiDS-Legacy). We combine COSEBIs (Complete Orthogonal Sets of E-/B-mode Integrals) at scales between 2' and 300' with third-order aperture mass moments at scales between 4' and 32' to perform a joint analysis of second- and third-order statistics. Compared to previous KiDS analyses, we implement several methodological advances: an intrinsic alignment model with redshift and mass dependence, a baryon correction model validated on multiple hydrodynamical simulations, and corrections for reduced shear and source clustering. Combining COSEBIs with third-order aperture mass statistics in KiDS-Legacy yields $Ω_m = 0.297^{+0.056}_{-0.040}$ and $S_8 = 0.806^{+0.025}_{-0.023}$, significantly tightening the $Ω_m$ constraints and more than doubling the figure of merit in the $Ω_m$--$S_8$ plane compared to the two-point analysis alone. The third-order measurements pass stringent internal consistency tests, are fully compatible with the KiDS-Legacy 2-point constraints, other 2+3-point lensing results and with Planck CMB measurements within $1σ$, providing no evidence for an $S_8$ tension and demonstrating the maturity of 3-point cosmic shear as a key probe for forthcoming surveys. △ Less

Submitted 10 June, 2026; originally announced June 2026.

Comments: 14 pages plus appendix, 15 figures, submitted to A&A

Euclid preparation. Three-dimensional galaxy clustering in configuration space: Three-point correlation function estimation

Authors: Euclid Collaboration, A. Veropalumbo, M. Moresco, F. Marulli, E. Branchini, M. Guidi, A. Farina, A. Pugno, E. Sefusatti, D. Tavagnacco, F. Rizzo, E. Romelli, S. de la Torre, A. Eggemeier, E. Sihvola, M. Viel, N. Aghanim, B. Altieri, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, S. Bardelli, P. Battaglia, A. Biviano , et al. (281 additional authors not shown)

Abstract: Higher-order correlation functions are firmly established as a fundamental tool for the statistical analysis of clustering in modern galaxy surveys. It was demonstrated that they greatly enrich the information content extracted by two-point statistics, allowing us to break the degeneracies between model parameters and constrain departures from Gaussianity. This paper presents the statistical estim… ▽ More Higher-order correlation functions are firmly established as a fundamental tool for the statistical analysis of clustering in modern galaxy surveys. It was demonstrated that they greatly enrich the information content extracted by two-point statistics, allowing us to break the degeneracies between model parameters and constrain departures from Gaussianity. This paper presents the statistical estimators adopted to evaluate the galaxy three-point correlation function and its numerical implementation within the data analysis pipeline of the Euclid Science Ground Segment. Two different algorithms are adopted to count triplets: a direct and exact counting method capable of providing a robust three-point correlation function measurement for any triangular configuration, and a more efficient method based on spherical harmonic decomposition, designed to address the computational challenges of measuring the three-point statistics for data sets as large as those of the final Euclid survey. The spherical harmonic decomposition estimates the Legendre coefficients of the three-point correlation function up to a finite expansion order. Despite being an approximation, the three-point function measured with this approach satisfies the scientific requirements of the mission. We also introduce, implement, and validate the random split technique, which reduces the computational cost of counting triplets in the reference random sample by a factor of 10, without significantly compromising numerical accuracy. We evaluated the robustness, precision, and accuracy of the numerical estimates through an extensive campaign of validation tests, the results of which are presented. Finally, we quantify the computational requirements and their scaling with the expected size of Euclid data set, showing that a complete three-point analysis of the final Euclid survey is within computational reach. △ Less

Submitted 4 May, 2026; originally announced May 2026.

Comments: 20 pages, 9 Figures

Euclid preparation. Impact of redshift distribution uncertainties on the joint analysis of photometric galaxy clustering and weak gravitational lensing

Authors: Euclid Collaboration, K. A. Bertmann, A. Porredon, V. Duret, J. Fonseca, H. Hildebrandt, I. Tutusaus, S. Camera, S. Escoffier, N. Aghanim, B. Altieri, A. Amara, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, S. Bardelli, P. Battaglia, A. Biviano, E. Branchini, M. Brescia, G. Cañas-Herrera, V. Capobianco, C. Carbone, V. F. Cardone , et al. (268 additional authors not shown)

Abstract: One of the $\textit{Euclid}$ mission's key projects is the so-called 3$\times$2pt analysis, that is, the combination of cosmic shear, photometric galaxy clustering, and galaxy-galaxy lensing. Although $\textit{Euclid}$ has established quality requirements for the photo-$z$ accuracy needed for the weak lensing galaxy sample, no such requirements have been set for the photometric clustering sample.… ▽ More One of the $\textit{Euclid}$ mission's key projects is the so-called 3$\times$2pt analysis, that is, the combination of cosmic shear, photometric galaxy clustering, and galaxy-galaxy lensing. Although $\textit{Euclid}$ has established quality requirements for the photo-$z$ accuracy needed for the weak lensing galaxy sample, no such requirements have been set for the photometric clustering sample. In this paper, we investigate the impact of redshift uncertainties on $\textit{Euclid}$'s photometric galaxy clustering analysis and its combination with weak gravitational lensing, focusing on data release 1 (DR1). In particular, we study whether having precise knowledge of the mean of the redshift distributions per bin is sufficient to avoid biases in the resulting cosmological constraints or whether accuracy in the higher-order moments of the distribution is required. We evaluate the results based on their constraining power on $w_{\mathrm{0}}$ and $w_{a}$ and define thresholds for the precision and accuracy of $\textit{Euclid}$'s redshift distribution of the photometric clustering sample. We find that the redshift distributions of the photometric clustering sample must be known at an accuracy of 0.004(1+$z$) in the mean in order to recover 80$\%$ of the constraining power in $\textit{Euclid}$'s DR1 $w_{\mathrm{0}}w_{a}$CDM 3$\times$2pt analysis. The impact of the uncertainty on the width is negligible, provided the mean redshift is constrained with sufficient accuracy. For most sources of redshift distribution error, attaining the requirement on the mean will also reduce uncertainty in the width well below the required level. △ Less

Submitted 1 April, 2026; originally announced April 2026.

Comments: 13+6 pages, 5+9 figures, 5+5 tables

Euclid preparation. Cosmology Likelihood for Observables in Euclid (CLOE). 2. Code implementation

Authors: Euclid Collaboration, S. Joudaki, V. Pettorino, L. Blot, M. Bonici, S. Camera, G. Cañas-Herrera, V. F. Cardone, P. Carrilho, S. Casas, S. Davini, S. Di Domizio, S. Farrens, L. W. K. Goh, S. Gouyou Beauchamps, S. Ilić, F. Keil, A. M. C. Le Brun, M. Martinelli, C. Moretti, A. Pezzotta, Z. Sakr, A. G. Sánchez, D. Sciotti, K. Tanidis , et al. (323 additional authors not shown)

Abstract: We provide a description of the code implementation and structure of Cosmology Likelihood for Observables in Euclid (CLOE), developed by members of the Euclid Consortium. CLOE is a modular Python code for computing the theoretical predictions of cosmological observables and evaluating them against state-of-the-art data from galaxy surveys such as Euclid in a unified likelihood. This primarily incl… ▽ More We provide a description of the code implementation and structure of Cosmology Likelihood for Observables in Euclid (CLOE), developed by members of the Euclid Consortium. CLOE is a modular Python code for computing the theoretical predictions of cosmological observables and evaluating them against state-of-the-art data from galaxy surveys such as Euclid in a unified likelihood. This primarily includes the core observables of weak gravitational lensing, photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering, but also extended probes such as the clusters of galaxies and cross-correlations of galaxy positions and shapes with the cosmic microwave background. While CLOE has been developed to serve as the unified framework for the parameter inferences in Euclid, it has general capabilities that can serve the broader cosmological community. It is different from other comparable cosmological tools in that it is written entirely in Python, performs the full likelihood calculation, and includes both photometric and spectroscopic observables. We will focus on the primary probes of Euclid and will describe the overall code structure, rigorous code development practices, extensive documentation, unique features, speed optimization, and future development plans. CLOE is publicly available at https://github.com/cloe-org/cloe. △ Less

Submitted 23 March, 2026; originally announced March 2026.

Comments: Second in a series of six papers presenting CLOE, the Euclid likelihood code; 43 pages, 11 figures, accepted for publication in A&A (with a table of contents added)

Euclid preparation. 3D reconstruction of the cosmic web with simulated Euclid Deep spectroscopic samples

Authors: Euclid Collaboration, K. Kraljic, C. Laigle, M. Balogh, P. Jablonka, U. Kuchner, N. Malavasi, F. Sarron, C. Pichon, G. De Lucia, M. Bethermin, F. Durret, M. Fumagalli, C. Gouin, M. Magliocchetti, J. G. Sorce, O. Cucciati, F. Fontanot, M. Hirschmann, Y. Kang, M. Spinelli, N. Aghanim, A. Amara, S. Andreon, N. Auricchio , et al. (299 additional authors not shown)

Abstract: The ongoing Euclid mission aims to measure spectroscopic redshifts for approximately two million galaxies using the H $α$ line emission detected in near-infrared slitless spectroscopic data from the Euclid Deep Fields (EDFs). These measurements will reach a flux limit of $5\times 10^{-17}\,{\rm erg}\,{\rm cm}^{-2}\,{\rm s}^{-1}$ in the redshift range $0.4<z<1.8$, opening the door to numerous inves… ▽ More The ongoing Euclid mission aims to measure spectroscopic redshifts for approximately two million galaxies using the H $α$ line emission detected in near-infrared slitless spectroscopic data from the Euclid Deep Fields (EDFs). These measurements will reach a flux limit of $5\times 10^{-17}\,{\rm erg}\,{\rm cm}^{-2}\,{\rm s}^{-1}$ in the redshift range $0.4<z<1.8$, opening the door to numerous investigations involving galaxy evolution, extending well beyond the mission's core objectives. The achieved H $α$ luminosity depth will lead to a sufficiently high sampling, enabling the reconstruction of the large-scale galaxy environment. We assess the quality of the reconstruction of the galaxy cosmic web environment with the expected spectroscopic dataset in EDFs. The analysis is carried out on the Flagship and GAEA galaxy mock catalogues. The quality of the reconstruction is first evaluated using geometrical and topological statistics measured on the cosmic web, namely the length of filaments, the area of walls, the volume of voids, and its connectivity and multiplicity. We then quantify how accurately gradients in galaxy properties with distance from filaments can be recovered. As expected, the small-scale redshift-space distortions, have a strong impact on filament lengths and connectivity, but can be mitigated by compressing galaxy groups before skeleton extraction. The cosmic web reconstruction is biased when relying solely on H $α$ emitters. This limitation can be mitigated by applying stellar mass weighting during the reconstruction. However, this approach introduces non-trivial biases that need to be accounted for when comparing to theoretical predictions. Redshift uncertainties pose the greatest challenge in recovering the expected dependence of galaxy properties, though the well-established stellar mass transverse gradients towards filaments can still be observed. △ Less

Submitted 15 January, 2026; originally announced January 2026.

Comments: 30 pages, 23 figures. Accepted for publication in A&A

Journal ref: A&A 708, A164 (2026)

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Authors: Euclid Collaboration, K. Hoffmann, R. Paviot, B. Joachimi, N. Tessore, P. Tallada-Crespí, N. E. Chisari, E. J. Gonzalez, A. Loureiro, P. Fosalba, J. Blazek, C. Laigle, Y. Dubois, C. Pichon, B. Altieri, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, S. Bardelli, F. Bernardeau, A. Biviano, E. Branchini, M. Brescia, S. Camera , et al. (276 additional authors not shown)

Abstract: Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes… ▽ More Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments. △ Less

Submitted 12 January, 2026; originally announced January 2026.

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

Authors: Euclid Collaboration, R. Paviot, B. Joachimi, K. Hoffmann, S. Codis, I. Tutusaus, D. Navarro-Gironés, J. Blazek, F. Hervas-Peters, B. Altieri, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, S. Bardelli, A. Biviano, E. Branchini, M. Brescia, S. Camera, G. Cañas-Herrera, V. Capobianco, C. Carbone, V. F. Cardone, J. Carretero, S. Casas , et al. (271 additional authors not shown)

Abstract: We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Li… ▽ More We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples. △ Less

Submitted 12 January, 2026; originally announced January 2026.

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

Authors: Euclid Collaboration, M. Kärcher, M. -A. Breton, S. de la Torre, A. Veropalumbo, A. Eggemeier, M. Crocce, E. Sefusatti, E. Sarpa, R. E. Angulo, B. Camacho Quevedo, L. Castiblanco, E. Castorina, A. Chudaykin, V. Desjacques, A. Farina, G. Gambardella, M. Guidi, D. Linde, F. Marulli, A. Moradinezhad Dizgah, M. Moresco, C. Moretti, K. Pardede, A. Pezzotta , et al. (287 additional authors not shown)

Abstract: The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-… ▽ More The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space. △ Less

Submitted 8 January, 2026; originally announced January 2026.

Comments: 29 pages, 13 figures, abstract abridged for arxiv submission

The Power of DESI for Photometric Redshift Calibration: A Case Study with KiDS-1000

Authors: Diana Blanco, Alexie Leauthaud, Johannes Ulf Lange, Angus H. Wright, Hendrik Hildebrandt, Sven Heydenreich, Darshika Ravulapalli, Joshua Ratajczak, Kyle S. Dawson, Jamie McCullough, Biprateep Dey, Jessica N. Aguilar, Steven Ahlen, Abhijeet Anand, Davide Bianchi, Chris Blake, David Brooks, Francisco J. Castander, Todd Claybaugh, Andrei Cuceu, Axel de la Macorra, John Della Costa, Arjun Dey, Ann Elliott, Ni Putu Audita Placida Emas , et al. (42 additional authors not shown)

Abstract: Accurate redshift estimates are a critical requirement for weak lensing surveys and one of the main uncertainties in constraints on dark energy and large-scale cosmic structure. In this paper, we study the potential to calibrate photometric redshift (photo-z) distributions for gravitational lensing using the Dark Energy Spectroscopic Instrument (DESI). Since beginning its science operations in 202… ▽ More Accurate redshift estimates are a critical requirement for weak lensing surveys and one of the main uncertainties in constraints on dark energy and large-scale cosmic structure. In this paper, we study the potential to calibrate photometric redshift (photo-z) distributions for gravitational lensing using the Dark Energy Spectroscopic Instrument (DESI). Since beginning its science operations in 2021, DESI has collected more than 50 million redshifts, adding about one million monthly. In addition to its large-scale structure samples, DESI has also acquired over 256k high-quality spectroscopic redshifts (spec-zs) in the COSMOS and XMM and VVDS fields. This is already a factor of 3 larger than previous spec-z calibration compilations in these two regions. Here, we explore calibrating photo-zs for the subset of KiDS-1000 galaxies that fall into joint self-organizing map (SOM) cells overlapping the DESI COSMOS footprint using the DESI COSMOS observations. Estimating the redshift distribution in KiDS-1000 with the new DESI data, we find broad consistency with previously published results while also detecting differences in the mean redshift in some tomographic bins with an average shifts of Delta Mean(z) = -0.028 in the mean and Delta Median(z) = +0.011 in the median across tomographic bins. However, we also find that incompleteness per SOM cell, i.e., groups of galaxies with similar colors and magnitudes, can modify n(z) distributions. Finally, we comment on the fact that larger photometric catalogs, aligned with the DESI COSMOS and DESI XMM and VVDS footprints, would be needed to fully exploit the DESI dataset and would extend the coverage to nearly eight times the area of existing 9-band photometry. △ Less

Submitted 21 December, 2025; v1 submitted 17 December, 2025; originally announced December 2025.

Comments: 21 pages, 14 figures

Clustering redshift distribution calibration of weak lensing surveys using the DESI-DR1 spectroscopic dataset

Authors: R. Ruggeri, C. Blake, J. Aguilar, S. Ahlen, D. Bianchi, D. Brooks, F. J. Castander, T. Claybaugh, A. Cuceu, K. S. Dawson, A. de la Macorra, B. Dey, P. Doel, A. Elliott, N. Emas, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, C. Garcia-Quintero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez, J. Guy, B. Hadzhiyska, H. K. Herrera-Alcantar , et al. (41 additional authors not shown)

Abstract: We estimate the source redshift distribution of current weak lensing surveys by applying the clustering-based redshift calibration technique, using the galaxy redshift sample provided by the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1). We cross-correlate the Bright Galaxy Survey (BGS), Luminous Red Galaxies (LRGs) and Emission Line Galaxies (ELGs) from DESI, within the redshift… ▽ More We estimate the source redshift distribution of current weak lensing surveys by applying the clustering-based redshift calibration technique, using the galaxy redshift sample provided by the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1). We cross-correlate the Bright Galaxy Survey (BGS), Luminous Red Galaxies (LRGs) and Emission Line Galaxies (ELGs) from DESI, within the redshift range $0.1 < z < 1.6$, with overlapping tomographic source samples from the Dark Energy Survey (DES), Kilo-Degree Survey (KiDS), and Hyper Suprime-Cam (HSC) survey. Using realistic mock catalogues, we test the stability of the clustering-redshift signal to fitting scale, reference-sample choice, and the evolution of source galaxy bias, and we explicitly model and marginalise over magnification contributions, which become non-negligible at $z \gtrsim 1$ due to the depth of the DESI ELG sample. We then compare the resulting bias-weighted redshift distributions to those calibrated using self-organising map (SOM) techniques, finding agreement within uncertainties for all surveys and tomographic bins. Our results demonstrate that clustering redshifts enabled by DESI's unprecedented spectroscopic sample provides a robust, complementary, and independent constraint capable of reducing one of the dominant systematic uncertainties in weak lensing cosmology. △ Less

Submitted 17 December, 2025; originally announced December 2025.

Comments: 20 pages, 10 figures

Cosmological Constraints from Full-Scale Clustering and Galaxy-Galaxy Lensing with DESI DR1

Authors: Johannes U. Lange, Alexandra Wells, Andrew Hearin, Gillian Beltz-Mohrmann, Alexie Leauthaud, Sven Heydenreich, Chris Blake, Jessica Nicole Aguilar, Steven Ahlen, Abhijeet Anand, Davide Bianchi, David Brooks, Francisco Javier Castander, Todd Claybaugh, Shaun Cole, Andrei Cuceu, Kyle Dawson, Axel de la Macorra, Biprateep Dey, Peter Doel, Ann Elliott, Ni Putu Audita Placida Emas, Simone Ferraro, Andreu Font-Ribera, Jaime E. Forero-Romero , et al. (54 additional authors not shown)

Abstract: We present constraints on cosmic structure growth from the analysis of galaxy clustering and galaxy-galaxy lensing with galaxies from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1. We analyze four samples drawn from the Bright Galaxy Survey (BGS) and the Luminous Red Galaxy (LRG) target classes. Projected galaxy clustering measurements from DESI are supplemented with lensing measu… ▽ More We present constraints on cosmic structure growth from the analysis of galaxy clustering and galaxy-galaxy lensing with galaxies from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1. We analyze four samples drawn from the Bright Galaxy Survey (BGS) and the Luminous Red Galaxy (LRG) target classes. Projected galaxy clustering measurements from DESI are supplemented with lensing measurements from the Dark Energy Survey (DES), the Kilo-Degree Survey (KiDS), and the Hyper Suprime-Cam (HSC) survey around the same targets. Our method relies on a simulation-based modeling framework using the AbacusSummit simulations and a complex halo occupation distribution model that incorporates assembly bias. We analyze scales down to $0.4 \, h^{-1} \, \mathrm{Mpc}$ for clustering and $2.5 \, h^{-1} \, \mathrm{Mpc}$ for lensing, leading to stringent constraints on $S_8 = σ_8 \sqrt{Ω_\mathrm{m} / 0.3}$ and $Ω_\mathrm{m}$ when fixing other cosmological parameters to those preferred by the CMB. We find $S_8 = 0.794 \pm 0.023$ and $Ω_\mathrm{m} = 0.295 \pm 0.012$ when using lensing measurements from DES and KiDS. Similarly, for HSC, we find $S_8 = 0.793 \pm 0.017$ and $Ω_\mathrm{m} = 0.303 \pm 0.010$ when assuming the best-fit photometric redshift offset suggested by the HSC collaboration. Overall, our results are in good agreement with other results in the literature while continuing to highlight the constraining power of non-linear scales. △ Less

Submitted 17 December, 2025; originally announced December 2025.

Comments: 16 pages, 7 figures

Joint cosmological fits to DESI-DR1 full-shape clustering and weak gravitational lensing in configuration space

Authors: A. Semenaite, C. Blake, A. Porredon, J. Aguilar, S. Ahlen, D. Bianchi, D. Brooks, F. J. Castander, T. Claybaugh, A. Cuceu, K. S. Dawson, A. de la Macorra, Biprateep Dey, P. Doel, A. Eggemeier, A. Elliott, N. Emas, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, C. Garcia-Quintero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez, J. Guy , et al. (47 additional authors not shown)

Abstract: We present a joint $3\times2$-pt cosmological analysis of auto- and cross-correlations between the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1) Bright Galaxy Survey (BGS) and Luminous Red Galaxy (LRG) samples and overlapping shear measurements from the KiDS-1000, DES-Y3 and HSC-Y3 weak lensing surveys. We perform our analysis in configuration space and, in addition to the cosmic… ▽ More We present a joint $3\times2$-pt cosmological analysis of auto- and cross-correlations between the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1) Bright Galaxy Survey (BGS) and Luminous Red Galaxy (LRG) samples and overlapping shear measurements from the KiDS-1000, DES-Y3 and HSC-Y3 weak lensing surveys. We perform our analysis in configuration space and, in addition to the cosmic shear correlation functions for each weak lensing dataset, we fit the tangential shear of the weak lensing source galaxies around DESI lens galaxies. Finally, we make use of the anisotropic BGS and LRG clustering information by fitting the full shape of the two-point correlation function multipoles measured over the full DESI-DR1 footprint, presenting the first full-shape analysis of DESI measurements in configuration space. We find that the addition of weak lensing information serves to improve, with respect to the clustering-only case, the measurements of the power spectrum amplitude parameters $\ln(10^{10}A_{\rm{s}})$ and $σ_{12}$ by $15\%$ and $36\%$, respectively. It also improves measurements of the linear bias of the lens galaxies by $15-20\%$, depending on the tracer. Our results show excellent consistency, regardless of the weak lensing survey considered, and are furthermore consistent with a companion analysis that fits $3\times2$-pt correlations including DESI projected clustering measurements, as well as the results published by the weak lensing collaborations themselves. Our measured values for weak lensing amplitude are $S_{8}^{\mathrm{DESI\times HSC}}=0.787\pm0.020$, $S_{8}^{\mathrm{DESI\times DES}}=0.791\pm0.016$, $S_{8}^{\mathrm{DESI\times KiDS}}=0.771\pm0.017$, which are $1.9σ-2.9σ$ below the $S_8$ value preferred by Planck. Finally, our clustering-only results are in good agreement with the Fourier space full-shape analysis of all DESI tracers. △ Less

Submitted 15 May, 2026; v1 submitted 17 December, 2025; originally announced December 2025.

Comments: 27 pages, 12 Figures; Published in the Open Journal of Astrophysics

DESI-DR1 $3 \times 2$-pt analysis: consistent cosmology across weak lensing surveys

Authors: A. Porredon, C. Blake, J. U. Lange, N. Emas, J. Aguilar, S. Ahlen, A. Bera, D. Bianchi, D. Brooks, F. J. Castander, T. Claybaugh, J. Coloma Nadal, A. Cuceu, K. S. Dawson, A. de la Macorra, Biprateep Dey, P. Doel, A. Elliott, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, C. Garcia-Quintero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez , et al. (48 additional authors not shown)

Abstract: We present a joint cosmological analysis of projected galaxy clustering observations from the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1), and overlapping weak gravitational lensing observations from three datasets: the Kilo-Degree Survey (KiDS-1000), the Dark Energy Survey (DES-Y3), and the Hyper-Suprime-Cam Survey (HSC-Y3). This combination of large-scale structure probes allo… ▽ More We present a joint cosmological analysis of projected galaxy clustering observations from the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1), and overlapping weak gravitational lensing observations from three datasets: the Kilo-Degree Survey (KiDS-1000), the Dark Energy Survey (DES-Y3), and the Hyper-Suprime-Cam Survey (HSC-Y3). This combination of large-scale structure probes allows us to measure a set of $3 \times 2$-pt correlation functions, breaking the degeneracies between parameters in cosmological fits to individual observables. We obtain mutually-consistent constraints on the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3} = 0.786^{+0.022}_{-0.019}$ from the combination of DESI-DR1 and DES-Y3, $S_8 = 0.760^{+0.020}_{-0.018}$ from KiDS-1000, and $S_8 = 0.771^{+0.026}_{-0.027}$ from HSC-Y3. These parameter determinations are consistent with fits to the Planck Cosmic Microwave Background dataset, albeit with $1.5-2σ$ lower values in the $S_8-Ω_{\rm m}$ plane. We perform our analysis with a unified pipeline tailored to the requirements of each cosmic shear survey, which self-consistently determines cosmological and astrophysical parameters. We generate an analytical covariance matrix for the correlation data including all cross-covariances between probes, and we design a new blinding procedure to safeguard our analysis against confirmation bias, whilst leaving goodness-of-fit statistics unchanged. Our study is part of a suite of papers that present joint cosmological analyses of DESI-DR1 and weak gravitational lensing datasets. △ Less

Submitted 27 April, 2026; v1 submitted 17 December, 2025; originally announced December 2025.

Comments: 21 pages, 15 figures; Published in the Open Journal of Astrophysics

Journal ref: The Open Journal of Astrophysics Vol. 9, 2026

KiDS-Legacy: Constraining dark energy, neutrino mass, and curvature

Authors: Robert Reischke, Benjamin Stölzner, Benjamin Joachimi, Angus H. Wright, Marika Asgari, Maciej Bilicki, Nora Elisa Chisari, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Joachim Harnois-Déraps, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Shahab Joudaki, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Arthur Loureiro, Constance Mahony, Lauro Moscardini, Lucas Porth, Mario Radovich, Tilman Tröster, Maximilian von Wietersheim-Kramsta , et al. (3 additional authors not shown)

Abstract: We constrained minimally extended cosmological models with the cosmic shear analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy) in combination with external probes. Due to the consistency of the KiDS-Legacy analysis with the cosmic microwave background (CMB), we could combine these datasets reliably for the first time. Additionally, we used CMB lensing, galaxy redshift-spa… ▽ More We constrained minimally extended cosmological models with the cosmic shear analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy) in combination with external probes. Due to the consistency of the KiDS-Legacy analysis with the cosmic microwave background (CMB), we could combine these datasets reliably for the first time. Additionally, we used CMB lensing, galaxy redshift-space distortions, and baryon acoustic oscillations. We assessed, in turn, the effects of spatial curvature, varying neutrino masses, and an evolving dark energy component on cosmological constraints from KiDS-Legacy alone and from KiDS-Legacy combined with external probes. We find KiDS-Legacy to be consistent with the fiducial flat $Λ$-cold dark matter ($Λ$CDM) analysis with $c^2 \sum m_ν\leq 1.5\,$eV, $w_0 = -1.0\pm 0.7$, and $w_a = -1.3^{+1.9}_{-2.0}$ while $Ω_K = 0.08^{+0.16}_{-0.17}$ (1$σ$ bounds) with an almost equal goodness of fit. The $w_0w_a$CDM model is not a significant improvement over $Λ$CDM when cosmic shear and CMB lensing are combined, yielding a Bayes factor $B = 0.07$. If all probes are combined, however, $B$ increases to 2.73, corresponding to a $2.6σ$ suspiciousness tension. The constraint on $S_8 = σ_8\sqrt{Ω_\mathrm{m}/0.3}$ is robust to opening up the parameter space for cosmic shear. Adding all external datasets to KiDS-Legacy, we find $S_8 = 0.816 \pm 0.006$ in $Λ$CDM and $S_8 = 0.837 \pm 0.008$ in $w_0 w_a$CDM for all probes combined. △ Less

Submitted 10 May, 2026; v1 submitted 11 December, 2025; originally announced December 2025.

Comments: 14 pages, 8 figures, published in A&A

Journal ref: A&A 709, A82 (2026)

KiDS-Legacy: Constraints on Horndeski gravity from weak lensing combined with galaxy clustering and cosmic microwave background anisotropies

Authors: Benjamin Stölzner, Robert Reischke, Matteo Grasso, Matteo Cataneo, Benjamin Joachimi, Arthur Loureiro, Alessio Spurio Mancini, Angus H. Wright, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Shahab Joudaki, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Constance Mahony, Lauro Moscardini, Lucas Porth, Mario Radovich, Tilman Tröster, Maximilian von Wietersheim-Kramsta , et al. (3 additional authors not shown)

Abstract: We present constraints on modified gravity from a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy) in combination with DESI measurements of baryon acoustic oscillations, eBOSS observations of redshift space distortions, and cosmic microwave background anisotropies from Planck. We study the Horndeski class of modified gravity models in an effective field theor… ▽ More We present constraints on modified gravity from a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy) in combination with DESI measurements of baryon acoustic oscillations, eBOSS observations of redshift space distortions, and cosmic microwave background anisotropies from Planck. We study the Horndeski class of modified gravity models in an effective field theory framework employing a parameterisation that satisfies stability conditions by construction and, for the first time, present a cosmological analysis in this inherently stable parameter basis. Cosmic shear constrains the Horndeski parameter space significantly, matching or surpassing the CMB contribution. Adopting the de-mixed kinetic term of the scalar field perturbation, $D_{\rm kin}$, and the deviation of the Planck mass from its fiducial value, $ΔM_*^2\equiv M_*^2-1$, as model parameters, we constrain their present values to be $Δ\hat{M}_*^2=0.32^{+0.07}_{-0.21}$ and $\hat{D}_{\rm kin} = 3.74^{+0.69}_{-1.92}$, which deviate from general relativity at $1.5σ$ and $1.9σ$, respectively. We derive constraints on the structure growth parameter $S_8=0.813^{+0.008}_{-0.011}$, which is compatible with the $Λ$CDM constraint at $0.54σ$. We obtain the deviation of the effective Newtonian coupling from the GR value as $Δμ_{\infty,{\rm eff}}=0.066\pm0.023$, corresponding to a $2.9σ$ significance. Although modified gravity provides a slightly better fit to the data, a model comparison shows only a weak preference for modified gravity at the $1.4σ$ level. When adopting a dynamical dark energy model of the background cosmology, the inferred modified gravity parameter constraints are stable with respect to a $Λ$CDM background, while a mild preference at $1.57σ$ for dynamical dark energy remains. △ Less

Submitted 24 March, 2026; v1 submitted 11 December, 2025; originally announced December 2025.

Comments: 21 pages, 13 figures, 5 tables, published in A&A

Journal ref: A&A 707, A323 (2026)

Euclid preparation. XCII. Controlling angular systematics in the Euclid spectroscopic galaxy sample

Authors: Euclid Collaboration, P. Monaco, M. Y. Elkhashab, B. R. Granett, J. Salvalaggio, E. Sefusatti, C. Scarlata, B. Zabelle, M. Bethermin, S. Bruton, C. Carbone, S. de la Torre, S. Dusini, A. Eggemeier, L. Guzzo, G. Lavaux, S. Lee, K. Markovic, K. S. McCarthy, M. Moresco, F. Passalacqua, W. J. Percival, I. Risso, A. G. Sánchez, D. Scott , et al. (276 additional authors not shown)

Abstract: We present the strategy used to identify and mitigate potential sources of angular systematics in the \textit{Euclid} spectroscopic galaxy survey, and we quantify their impact on galaxy clustering measurements and cosmological parameter estimation. We first surveyed the \textit{Euclid} processing pipeline to identify all evident, potential sources of systematics, and classified them into two broad… ▽ More We present the strategy used to identify and mitigate potential sources of angular systematics in the \textit{Euclid} spectroscopic galaxy survey, and we quantify their impact on galaxy clustering measurements and cosmological parameter estimation. We first surveyed the \textit{Euclid} processing pipeline to identify all evident, potential sources of systematics, and classified them into two broad classes: angular systematics, which modulate the galaxy number density across the sky, and catastrophic redshift errors, which lead to interlopers in the galaxy sample. We then used simulated spectroscopic surveys to test our ability to mitigate angular systematics by constructing a random catalogue that represents the `visibility mask' of the survey; this is a dense set of intrinsically unclustered objects, subject to the same selection effects as the data catalogue. The construction of this random catalogue relies on a detection model, which gives the probability of reliably measuring the galaxy redshift as a function of the signal-to-noise ratio (S/N) of its emission lines. We demonstrate that, in the ideal case of a perfect knowledge of the visibility mask, the galaxy power spectrum in the presence of systematics is recovered, to within sub-per cent accuracy, by convolving a theory power spectrum with a window function obtained from the random catalogue itself. In the case of only approximate knowledge of the visibility mask, we test the stability of power spectrum measurements and cosmological parameter posteriors by using perturbed versions of the random catalogue. We find that significant effects are limited to very large scales, and parameter estimation remains robust; the most impacting effects are connected to the calibration of the detection model. △ Less

Submitted 27 March, 2026; v1 submitted 25 November, 2025; originally announced November 2025.

Comments: A&A in press; 36 pages, 18 figures, 4 appendices

A Unified Photometric Redshift Calibration for Weak Lensing Surveys using the Dark Energy Spectroscopic Instrument

Authors: Johannes U. Lange, Diana Blanco, Alexie Leauthaud, Angus Wright, Abigail Fisher, Joshua Ratajczak, Jessica Nicole Aguilar, Steven Ahlen, Stephen Bailey, Davide Bianchi, Chris Blake, David Brooks, Todd Claybaugh, Andrei Cuceu, Kyle Dawson, Axel de la Macorra, Joseph DeRose, Arjun Dey, Peter Doel, Ni Putu Audita Placida Emas, Simone Ferraro, Andreu Font-Ribera, Jaime E. Forero-Romero, Cristhian Garcia-Quintero, Enrique Gaztañaga , et al. (39 additional authors not shown)

Abstract: The effective redshift distribution $n(z)$ of galaxies is a critical component in the study of weak gravitational lensing. Here, we introduce a new method for determining $n(z)$ for weak lensing surveys based on high-quality redshifts and neural network-based importance weights. Additionally, we present the first unified photometric redshift calibration of the three leading stage-III weak lensing… ▽ More The effective redshift distribution $n(z)$ of galaxies is a critical component in the study of weak gravitational lensing. Here, we introduce a new method for determining $n(z)$ for weak lensing surveys based on high-quality redshifts and neural network-based importance weights. Additionally, we present the first unified photometric redshift calibration of the three leading stage-III weak lensing surveys, the Dark Energy Survey (DES), the Hyper Suprime-Cam (HSC) survey and the Kilo-Degree Survey (KiDS), with state-of-the-art spectroscopic data from the Dark Energy Spectroscopic Instrument (DESI). We verify our method using a new, data-driven approach and obtain $n(z)$ constraints with statistical uncertainties of order $σ_{\bar z} \sim 0.01$ and smaller. Our analysis is largely independent of previous photometric redshift calibrations and, thus, provides an important cross-check in light of recent cosmological tensions. Overall, we find excellent agreement with previously published results on the DES Y3 and HSC Y1 data sets while there are some differences on the mean redshift with respect to the previously published KiDS-1000 results. We attribute the latter to mismatches in photometric noise properties in the COSMOS field compared to the wider KiDS SOM-gold catalog. At the same time, the new $n(z)$ estimates for KiDS do not significantly change estimates of cosmic structure growth from cosmic shear. Finally, we discuss how our method can be applied to future weak lensing calibrations with DESI data. △ Less

Submitted 26 December, 2025; v1 submitted 29 October, 2025; originally announced October 2025.

Comments: 23 pages, 9 figures, accepted for publication in ApJ

Euclid preparation. Cosmology Likelihood for Observables in Euclid (CLOE). 6: Impact of systematic uncertainties on the cosmological analysis

Authors: Euclid Collaboration, L. Blot, K. Tanidis, G. Cañas-Herrera, P. Carrilho, M. Bonici, S. Camera, V. F. Cardone, S. Casas, S. Davini, S. Di Domizio, S. Farrens, L. W. K. Goh, S. Gouyou Beauchamps, S. Ilić, S. Joudaki, F. Keil, A. M. C. Le Brun, M. Martinelli, C. Moretti, V. Pettorino, A. Pezzotta, Z. Sakr, A. G. Sánchez, D. Sciotti , et al. (287 additional authors not shown)

Abstract: Extracting cosmological information from the Euclid galaxy survey will require modelling numerous systematic effects during the inference process. This implies varying a large number of nuisance parameters, which have to be marginalised over before reporting the constraints on the cosmological parameters. This is a delicate process, especially with such a large parameter space, which could result… ▽ More Extracting cosmological information from the Euclid galaxy survey will require modelling numerous systematic effects during the inference process. This implies varying a large number of nuisance parameters, which have to be marginalised over before reporting the constraints on the cosmological parameters. This is a delicate process, especially with such a large parameter space, which could result in biased cosmological results. In this work, we study the impact of different choices for modelling systematic effects and prior distribution of nuisance parameters for the final Euclid Data Release, focusing on the 3$\times$2pt analysis for photometric probes and the galaxy power spectrum multipoles for the spectroscopic probes. We explore the effect of intrinsic alignments, linear galaxy bias, magnification bias, multiplicative cosmic shear bias and shifts in the redshift distribution for the photometric probes, as well as the purity of the spectroscopic sample. We find that intrinsic alignment modelling has the most severe impact with a bias up to $6\,σ$ on the Hubble constant $H_0$ if neglected, followed by mis-modelling of the redshift evolution of galaxy bias, yielding up to $1.5\,σ$ on the parameter $S_8\equivσ_8\sqrt{Ω_{\rm m} /0.3}$. Choosing a too optimistic prior for multiplicative bias can also result in biases of the order of $0.7\,σ$ on $S_8$. We also find that the precision on the estimate of the purity of the spectroscopic sample will be an important driver for the constraining power of the galaxy clustering full-shape analysis. These results will help prioritise efforts to improve the modelling and calibration of systematic effects in Euclid. △ Less

Submitted 11 October, 2025; originally announced October 2025.

Comments: 22 pages, 16 figures. Submitted to A&A

Euclid preparation. XCVI. Cosmology Likelihood for Observables in Euclid (CLOE). 3. Inference and Forecasts

Authors: Euclid Collaboration, G. Cañas-Herrera, L. W. K. Goh, L. Blot, M. Bonici, S. Camera, V. F. Cardone, P. Carrilho, S. Casas, S. Davini, S. Di Domizio, S. Farrens, S. Gouyou Beauchamps, S. Ilić, S. Joudaki, F. Keil, A. M. C. Le Brun, M. Martinelli, C. Moretti, V. Pettorino, A. Pezzotta, Z. Sakr, A. G. Sánchez, D. Sciotti, K. Tanidis , et al. (315 additional authors not shown)

Abstract: The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scal… ▽ More The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scale structure observables improves. In this study, we present the first scientific forecasts using CLOE (Cosmology Likelihood for Observables in Euclid), a dedicated Euclid cosmological pipeline developed to support this endeavour. Using advanced Bayesian inference techniques applied to synthetic Euclid-like data, we sample the posterior distribution of cosmological and nuisance parameters across a variety of cosmological models and Euclid primary probes: cosmic shear, angular photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering. We validate the capability of CLOE to produce reliable cosmological forecasts, showcasing Euclid's potential to achieve a figure of merit for the dark energy parameters $w_0$ and $w_a$ exceeding 400 when combining all primary probes. Furthermore, we illustrate the behaviour of the posterior probability distribution of the parameters of interest given different priors and scale cuts. Finally, we emphasise the importance of addressing computational challenges, proposing further exploration of innovative data science techniques to efficiently navigate the Euclid high-dimensional parameter space in upcoming cosmological data releases. △ Less

Submitted 6 May, 2026; v1 submitted 10 October, 2025; originally announced October 2025.

Comments: Third in a series of six papers presenting CLOE, the Euclid likelihood code; 39 pages, 21 figures, submitted to A&A

Euclid preparation. XCVIII. Cosmology Likelihood for Observables in Euclid (CLOE). 5: Extensions beyond the standard modelling of theoretical probes and systematic effects

Authors: Euclid Collaboration, L. W. K. Goh, A. Nouri-Zonoz, S. Pamuk, M. Ballardini, B. Bose, G. Cañas-Herrera, S. Casas, G. Franco-Abellán, S. Ilić, F. Keil, M. Kunz, A. M. C. Le Brun, F. Lepori, M. Martinelli, Z. Sakr, F. Sorrenti, E. M. Teixeira, I. Tutusaus, L. Blot, M. Bonici, C. Bonvin, S. Camera, V. F. Cardone, P. Carrilho , et al. (279 additional authors not shown)

Abstract: Euclid is expected to establish new state-of-the-art constraints on extensions beyond the standard LCDM cosmological model by measuring the positions and shapes of billions of galaxies. Specifically, its goal is to shed light on the nature of dark matter and dark energy. Achieving this requires developing and validating advanced statistical tools and theoretical prediction software capable of test… ▽ More Euclid is expected to establish new state-of-the-art constraints on extensions beyond the standard LCDM cosmological model by measuring the positions and shapes of billions of galaxies. Specifically, its goal is to shed light on the nature of dark matter and dark energy. Achieving this requires developing and validating advanced statistical tools and theoretical prediction software capable of testing extensions of the LCDM model. In this work, we describe how the Euclid likelihood pipeline, Cosmology Likelihood for Observables in Euclid (CLOE), has been extended to accommodate alternative cosmological models and to refine the theoretical modelling of Euclid primary probes. In particular, we detail modifications made to CLOE to incorporate the magnification bias term into the spectroscopic two-point correlation function of galaxy clustering. Additionally, we explain the adaptations made to CLOE's implementation of Euclid primary photometric probes to account for massive neutrinos and modified gravity extensions. Finally, we present the validation of these CLOE modifications through dedicated forecasts on synthetic Euclid-like data by sampling the full posterior distribution and comparing with the results of previous literature. In conclusion, we have identified in this work several functionalities with regards to beyond-LCDM modelling that could be further improved within CLOE, and outline potential research directions to enhance pipeline efficiency and flexibility through novel inference and machine learning techniques. △ Less

Submitted 6 May, 2026; v1 submitted 10 October, 2025; originally announced October 2025.

Comments: Accepted for publication in A&A

Euclid preparation. Cosmology Likelihood for Observables in Euclid (CLOE). 4: Validation and Performance

Authors: Euclid Collaboration, M. Martinelli, A. Pezzotta, D. Sciotti, L. Blot, M. Bonici, S. Camera, G. Cañas-Herrera, V. F. Cardone, P. Carrilho, S. Casas, S. Davini, S. Di Domizio, S. Farrens, L. W. K. Goh, S. Gouyou Beauchamps, S. Ilić, S. Joudaki, F. Keil, A. M. C. Le Brun, C. Moretti, V. Pettorino, A. G. Sánchez, Z. Sakr, K. Tanidis , et al. (312 additional authors not shown)

Abstract: The Euclid satellite will provide data on the clustering of galaxies and on the distortion of their measured shapes, which can be used to constrain and test the cosmological model. However, the increase in precision places strong requirements on the accuracy of the theoretical modelling for the observables and of the full analysis pipeline. In this paper, we investigate the accuracy of the calcula… ▽ More The Euclid satellite will provide data on the clustering of galaxies and on the distortion of their measured shapes, which can be used to constrain and test the cosmological model. However, the increase in precision places strong requirements on the accuracy of the theoretical modelling for the observables and of the full analysis pipeline. In this paper, we investigate the accuracy of the calculations performed by the Cosmology Likelihood for Observables in Euclid (CLOE), a software able to handle both the modelling of observables and their fit against observational data for both the photometric and spectroscopic surveys of Euclid, by comparing the output of CLOE with external codes used as benchmark. We perform such a comparison on the quantities entering the calculations of the observables, as well as on the final outputs of these calculations. Our results highlight the high accuracy of CLOE when comparing its calculation against external codes for Euclid observables on an extended range of operative cases. In particular, all the summary statistics of interest always differ less than $0.1\,σ$ from the chosen benchmark, and CLOE predictions are statistically compatible with simulated data obtained from benchmark codes. The same holds for the comparison of correlation function in configuration space for spectroscopic and photometric observables. △ Less

Submitted 10 October, 2025; originally announced October 2025.

Comments: Fourth in a series of six papers presenting CLOE, the Euclid likelihood code; 19 pages, 13 figures, submitted to A&A

Cosmology Likelihood for Observables in \Euclid (CLOE). 1. Theoretical recipe

Authors: Euclid Collaboration, V. F. Cardone, S. Joudaki, L. Blot, M. Bonici, S. Camera, G. Cañas-Herrera, P. Carrilho, S. Casas, S. Davini, S. Di Domizio, S. Farrens, L. W. K. Goh, S. Gouyou Beauchamps, S. Ilić, F. Keil, A. M. C. Le Brun, M. Martinelli, C. Moretti, V. Pettorino, A. Pezzotta, A. G. Sánchez, Z. Sakr, D. Sciotti, K. Tanidis , et al. (301 additional authors not shown)

Abstract: As the statistical precision of cosmological measurements increases, the accuracy of the theoretical description of these measurements needs to increase correspondingly in order to infer the underlying cosmology that governs the Universe. To this end, we have created the Cosmology Likelihood for Observables in Euclid (CLOE), which is a novel cosmological parameter inference pipeline developed with… ▽ More As the statistical precision of cosmological measurements increases, the accuracy of the theoretical description of these measurements needs to increase correspondingly in order to infer the underlying cosmology that governs the Universe. To this end, we have created the Cosmology Likelihood for Observables in Euclid (CLOE), which is a novel cosmological parameter inference pipeline developed within the Euclid Consortium to translate measurements and covariances into cosmological parameter constraints. In this first in a series of six papers, we describe the theoretical recipe of this code for the Euclid primary probes. These probes are composed of the photometric 3x2pt observables of cosmic shear, galaxy-galaxy lensing, and galaxy clustering, along with spectroscopic galaxy clustering. We provide this description in both Fourier and configuration space for standard and extended summary statistics, including the wide range of systematic uncertainties that affect them. This includes systematic uncertainties such as intrinsic galaxy alignments, baryonic feedback, photometric and spectroscopic redshift uncertainties, shear calibration uncertainties, sample impurities, photometric and spectroscopic galaxy biases, as well as magnification bias. The theoretical descriptions are further able to accommodate both Gaussian and non-Gaussian likelihoods and extended cosmologies with non-zero curvature, massive neutrinos, evolving dark energy, and simple forms of modified gravity. These theoretical descriptions that underpin CLOE will form a crucial component in revealing the true nature of the Universe with next-generation cosmological surveys such as Euclid. △ Less

Submitted 10 October, 2025; originally announced October 2025.

Comments: First in a series of six papers presenting CLOE, the Euclid likelihood code; 27 pages, 10 figures, A&A submitted

Validation of the DESI-DR1 3x2-pt analysis: scale cut and shear ratio tests

Authors: N. Emas, A. Porredon, C. Blake, J. DeRose, J. Aguilar, S. Ahlen, D. Bianchi, D. Brooks, F. J. Castander, T. Claybaugh, A. Cuceu, A. de la Macorra, A. Dey, B. Dey, P. Doel, S. Ferraro, J. E. Forero-Romero, C. Garcia-Quintero, E. Gaztañaga, S. Gontcho A Gontcho, G. Gutierrez, S. Heydenreich, K. Honscheid, D. Huterer, M. Ishak , et al. (32 additional authors not shown)

Abstract: Combined survey analyses of galaxy clustering and weak gravitational lensing (3x2-pt studies) will allow new and accurate tests of the standard cosmological model. However, careful validation is necessary to ensure that these cosmological constraints are not biased by uncertainties associated with the modelling of astrophysical or systematic effects. In this study we validate the combined 3x2-pt a… ▽ More Combined survey analyses of galaxy clustering and weak gravitational lensing (3x2-pt studies) will allow new and accurate tests of the standard cosmological model. However, careful validation is necessary to ensure that these cosmological constraints are not biased by uncertainties associated with the modelling of astrophysical or systematic effects. In this study we validate the combined 3x2-pt analysis of the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1) spectroscopic galaxy clustering and overlapping weak lensing datasets from the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES), and the Hyper-Suprime-Cam Survey (HSC). By propagating the modelling uncertainties associated with the non-linear matter power spectrum, non-linear galaxy bias and baryon feedback, we design scale cuts to ensure that measurements of the matter density and the amplitude of the matter power spectrum are biased by less than 30% of the statistical error. We also test the internal consistency of the data and weak lensing systematics by performing new measurements of the lensing shear ratio. We demonstrate that the DESI-DR1 shear ratios can be successfully fit by the same model used to describe cosmic shear correlations, and analyse the additional information that can be extracted about the source redshift distributions and intrinsic alignment parameters. This study serves as crucial preparation for the upcoming cosmological parameter analysis of these datasets. △ Less

Submitted 29 March, 2026; v1 submitted 6 October, 2025; originally announced October 2025.

Comments: 21 pages, 11 figures, Published in the Open Journal of Astrophysics

Taking the Weight Off: Mitigating Parameter Bias from Catastrophic Outliers in 3$\times$2pt Analysis

Authors: Carolyn McDonald Mill, C. Danielle Leonard, Markus Michael Rau, Cora Uhlemann, Shahab Joudaki

Abstract: Stage IV cosmological surveys will map the universe with unprecedented precision, reducing statistical uncertainties to levels where unmodelled systematics can significantly bias inference. In particular, photometric redshift (photo-z) errors and intrinsic alignments (IA) must be robustly accounted for to ensure accurate inference of cosmological parameters. The increasing depth of Stage IV survey… ▽ More Stage IV cosmological surveys will map the universe with unprecedented precision, reducing statistical uncertainties to levels where unmodelled systematics can significantly bias inference. In particular, photometric redshift (photo-z) errors and intrinsic alignments (IA) must be robustly accounted for to ensure accurate inference of cosmological parameters. The increasing depth of Stage IV surveys exacerbates these challenges by producing low signal-to-noise galaxy populations prone to inaccurate photo-z measurements. Catastrophically misidentified redshifts are especially problematic for 3$\times$2pt inferences that combine weak lensing and galaxy clustering information. We demonstrate that even modest outlier fractions (e.g. 5%) can lead to substantial biases in cosmological parameter estimates: up to 1.8$σ$ in $Ω_M$ and $σ_8$, and over 8$σ$ in the IA redshift evolution parameter $η$. To address this, we introduce a flexible weighting scheme at the likelihood level that down-weights the most contamination-sensitive elements of the data vector during inference. This method mitigates biases without inflating the parameter space, reducing cosmological parameter biases to below 1$σ$ without substantially degrading constraining power. Our approach offers a practical solution for future analyses, enabling robust cosmological inference in the presence of catastrophic redshift errors. △ Less

Submitted 10 December, 2025; v1 submitted 9 September, 2025; originally announced September 2025.

Comments: 25 pages, 10 figures

Testing gravitational physics by combining DESI DR1 and weak lensing datasets using the E_G estimator

Authors: S. J. Rauhut, C. Blake, U. Andrade, H. E. Noriega, J. Aguilar, S. Ahlen, S. BenZvi, D. Bianchi, D. Brooks, T. Claybaugh, A. Cuceu, A. de la Macorra, J. DeRose, P. Doel, N. Emas, S. Ferraro, J. E. Forero-Romero, C. Garcia-Quintero, E. Gaztañaga, G. Gutierrez, S. Heydenreich, K. Honscheid, C. Howlett, D. Huterer, M. Ishak , et al. (39 additional authors not shown)

Abstract: The action of gravitational physics across space-time creates observable signatures in the behaviour of light and matter. We perform combined-probe studies using data from the Baryon Oscillation Spectroscopic Survey (BOSS) and Dark Energy Spectroscopic Instrument survey Data Release 1 (DESI-DR1), in combination with three existing weak lensing surveys, the Kilo-Degree Survey (KiDS), the Dark Energ… ▽ More The action of gravitational physics across space-time creates observable signatures in the behaviour of light and matter. We perform combined-probe studies using data from the Baryon Oscillation Spectroscopic Survey (BOSS) and Dark Energy Spectroscopic Instrument survey Data Release 1 (DESI-DR1), in combination with three existing weak lensing surveys, the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES), and the Hyper Suprime-Cam Survey (HSC), to test and constrain General Relativity (GR) in the context of the standard model of cosmology (LCDM). We focus on measuring the gravitational estimator statistic, E_G, which describes the relative amplitudes of weak gravitational lensing and galaxy velocities induced by a common set of overdensities. By comparing our amplitude measurements with their predicted scale- and redshift-dependence within the GR+LCDM model, we demonstrate that our results are consistent with the predictions of the Planck cosmology. The redshift span of the DESI dataset allows us to perform these E_G measurements at the highest redshifts achieved to date, z ~ 1. △ Less

Submitted 2 October, 2025; v1 submitted 21 July, 2025; originally announced July 2025.

Comments: 20 pages, 9 figures, published in the Open Journal of Astrophysics

Euclid preparation: Expected constraints on initial conditions

Authors: Euclid Collaboration, F. Finelli, Y. Akrami, A. Andrews, M. Ballardini, S. Casas, D. Karagiannis, Z. Sakr, J. Valiviita, G. Alestas, N. Bartolo, J. R. Bermejo-Climent, S. Nesseris, D. Paoletti, D. Sapone, I. Tutusaus, A. Achúcarro, G. Cañas-Herrera, J. Jasche, G. Lavaux, N. Aghanim, B. Altieri, A. Amara, L. Amendola, S. Andreon , et al. (285 additional authors not shown)

Abstract: The Euclid mission of the European Space Agency will deliver galaxy and cosmic shear surveys, which will be used to constrain initial conditions and statistics of primordial fluctuations. We present highlights for the Euclid scientific capability to test initial conditions beyond LCDM with the main probes, i.e. 3D galaxy clustering from the spectroscopic survey, the tomographic approach to 3x2pt s… ▽ More The Euclid mission of the European Space Agency will deliver galaxy and cosmic shear surveys, which will be used to constrain initial conditions and statistics of primordial fluctuations. We present highlights for the Euclid scientific capability to test initial conditions beyond LCDM with the main probes, i.e. 3D galaxy clustering from the spectroscopic survey, the tomographic approach to 3x2pt statistics from photometric galaxy survey, and their combination. We provide Fisher forecasts from the combination of Euclid spectroscopic and photometric surveys for spatial curvature, running of the spectral index of the power spectrum of curvature perturbations, isocurvature perturbations, and primordial features. For the parameters of these models we also provide the combination of Euclid forecasts (pessimistic and optimistic) with current and future measurements of the cosmic microwave background (CMB) anisotropies., i.e. Planck, the Simons Observatory (SO), and CMB-S4. We provide Fisher forecasts for how the power spectrum and bispectrum from the Euclid spectroscopic survey will constrain the local, equilateral, and orthogonal shapes of primordial non-Gaussianity. We also review how Bayesian field-level inference of primordial non-Gaussianity can constrain local primordial non-Gaussianity. We show how Euclid, with its unique combination of the main probes, will provide the tightest constraints on low redshift to date. By targeting a markedly different range in redshift and scale, Euclid's expected uncertainties are complementary to those obtained by CMB primary anisotropy, returning the tightest combined constraints on the physics of the early Universe. △ Less

Submitted 21 July, 2025; originally announced July 2025.

Comments: Abstract abridged, 25 pages, 6 tables, 11 figures

AMICO galaxy clusters in KiDS-1000: cosmological sample

Authors: M. Maturi, M. Radovich, L. Moscardini, G. F. Lesci, G. Castignani, F. Marulli, E. A. Puddu, M. Romanello, M. Sereno, C. Giocoli, L. Ingoglia, S. Bardelli, B. Giblin, H. Hildebrandt, S. Joudaki

Abstract: Context. Galaxy clusters provide key insights into cosmic structure formation, galaxy formation and are essential for cosmological studies. Aims. We present a catalog of galaxy clusters detected in the Kilo-Degree Survey (KiDS-DR4) optimized for cosmological analyses and investigations of cluster properties. Each detection includes probabilistic membership assignments for the KiDS-DR4 galaxies wit… ▽ More Context. Galaxy clusters provide key insights into cosmic structure formation, galaxy formation and are essential for cosmological studies. Aims. We present a catalog of galaxy clusters detected in the Kilo-Degree Survey (KiDS-DR4) optimized for cosmological analyses and investigations of cluster properties. Each detection includes probabilistic membership assignments for the KiDS-DR4 galaxies within the magnitude range $15<r'<24$. Methods. Using the Adaptive Matched Identifier of Clustered Objects (AMICO) algorithm, we identified 23965 clusters over an effective area of about 839 deg2 in the redshift range $0.1\le z \le0.9$, with a signal-to-noise ratio $S/N>3.5$. The sample is highly homogeneous across the entire survey thanks to the restrictive galaxy selection criteria we adopted. Spectroscopic data from the GAMA survey were used to calibrate the clusters photometric redshift and assess their uncertainties. We introduced algorithmic enhancements to AMICO to mitigate border effects among neighbor tiles. Quality flags are also provided for each cluster detection. The sample purity and completeness assessments have been estimated using the SinFoniA data driven approach, thus avoiding strong assumptions embedded in numerical simulations. We introduced a blinding scheme of the selection function meant to support the cosmological analyses. Results. Our cluster sample includes 321 cross-matches with the X-ray eRASS1 "primary" sample and 235 matches with the ACT-DR5 cluster sample. We derived a mass-proxy scaling relation based on intrinsic richness, $λ_*$, using masses from the eRASS1 catalog. Conclusions. The KiDS-DR4 cluster catalog provides a valuable data set for investigating galaxy cluster properties and contributes to cosmological studies by offering a large, well-characterized cluster sample. △ Less

Submitted 23 July, 2025; v1 submitted 18 July, 2025; originally announced July 2025.

Comments: 23 pages, 24 figures, submitted to A&A

AMICO galaxy clusters in KiDS-1000: cosmological constraints and mass calibration from counts and weak lensing

Authors: G. F. Lesci, F. Marulli, L. Moscardini, M. Maturi, M. Sereno, M. Radovich, M. Romanello, C. Giocoli, A. H. Wright, S. Bardelli, M. Bilicki, G. Castignani, H. Hildebrandt, L. Ingoglia, S. Joudaki, A. Kannawadi, E. Puddu

Abstract: We present the joint modelling of weak-lensing and count measurements of the galaxy clusters detected with the AMICO code, in the fourth data release of the Kilo Degree Survey (KiDS-1000). The analysed sample comprises about 8000 clusters, covering an effective area of 839 deg$^{2}$ and extending up to a redshift of $z = 0.8$. Stacked cluster weak-lensing and count measurements have been derived i… ▽ More We present the joint modelling of weak-lensing and count measurements of the galaxy clusters detected with the AMICO code, in the fourth data release of the Kilo Degree Survey (KiDS-1000). The analysed sample comprises about 8000 clusters, covering an effective area of 839 deg$^{2}$ and extending up to a redshift of $z = 0.8$. Stacked cluster weak-lensing and count measurements have been derived in bins of redshift and intrinsic richness, $λ^*$. Based on self-organising maps, we reconstructed the true redshift distributions of the background galaxy samples. We accounted for the systematic uncertainties arising from impurities in the background and cluster samples, biases in the cluster $z$ and $λ^*$, projection effects, halo orientation and miscentring, truncation of cluster halo mass distributions, matter correlated with cluster haloes, multiplicative shear bias, baryonic matter, geometric distortions in the lensing profiles, uncertainties in the theoretical halo mass function, and super-sample covariance. We also employed a blinding strategy based on perturbing the cluster sample completeness. The improved statistics and photometry compared to the previous KiDS data release, KiDS-DR3, have led to a halving of the uncertainties on $Ω_{\rm m}$ and $σ_8$, as we obtained $Ω_{\rm m}=0.22\pm0.02$ and $σ_8=0.86\pm0.03$. The constraint on $S_8 \equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$, $S_8=0.74\pm0.03$, is in excellent agreement with recent cluster count and KiDS-1000 cosmic shear analyses, while it shows a $2.8σ$ tension with Planck cosmic microwave background results. The constraints on the $\logλ^*-\log M_{200}$ relation imply a mass precision of 8%, on average. In addition, the result on the intrinsic scatter of the $\logλ^*-\log M_{200}$ relation, $σ_{\rm intr}=0.05\pm0.02$, confirms that $λ^*$ is an excellent mass proxy. △ Less

Submitted 25 September, 2025; v1 submitted 18 July, 2025; originally announced July 2025.

Comments: 26 pages, 19 figures. Submitted to A&A. We addressed the first referee report

Journal ref: A&A 703, A25 (2025)

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

Authors: Euclid Collaboration, P. Monaco, G. Parimbelli, M. Y. Elkhashab, J. Salvalaggio, T. Castro, M. D. Lepinzan, E. Sarpa, E. Sefusatti, L. Stanco, L. Tornatore, G. E. Addison, S. Bruton, C. Carbone, F. J. Castander, J. Carretero, S. de la Torre, P. Fosalba, G. Lavaux, S. Lee, K. Markovic, K. S. McCarthy, F. Passalacqua, W. J. Percival, I. Risso , et al. (281 additional authors not shown)

Abstract: We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating… ▽ More We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1). △ Less

Submitted 26 September, 2025; v1 submitted 16 July, 2025; originally announced July 2025.

Comments: 19 pages, A&A in press

Journal ref: A&A 704, A306 (2025)

Intrinsic alignment demographics for next-generation lensing: Revealing galaxy property trends with DESI Y1 direct measurements

Authors: J. Siegel, J. McCullough, A. Amon, C. Lamman, N. Jeffrey, B. Joachimi, H. Hoekstra, S. Heydenreich, A. J. Ross, J. Aguilar, S. Ahlen, D. Bianchi, C. Blake, D. Brooks, F. J. Castander, T. Claybaugh, A. de la Macorra, J. DeRose, P. Doel, N. Emas, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho , et al. (38 additional authors not shown)

Abstract: We present direct measurements of the intrinsic alignments (IA) of over 2 million spectroscopic galaxies using DESI Data Release 1 and imaging from four lensing surveys: DES, HSC, KiDS, and SDSS. In this uniquely data-rich regime, we take initial steps towards a more tailored IA modelling approach by building a library of IA measurements across colour, luminosity, stellar mass, and redshift. We ma… ▽ More We present direct measurements of the intrinsic alignments (IA) of over 2 million spectroscopic galaxies using DESI Data Release 1 and imaging from four lensing surveys: DES, HSC, KiDS, and SDSS. In this uniquely data-rich regime, we take initial steps towards a more tailored IA modelling approach by building a library of IA measurements across colour, luminosity, stellar mass, and redshift. We map the dependence between galaxy type -- in terms of rest-frame colour, strength of the 4000 Angstrom break, and specific star formation rate -- and IA amplitude; the bluest galaxies have an alignment consistent with zero, across low ($0.05<z<0.5$) and high ($0.8<z<1.55$) redshifts. In order to construct cosmic shear samples that are minimally impacted by IA but maintain maximum sample size and statistical power, we map the dependence of alignment with colour purity. Red, quenched galaxies are strongly aligned and the amplitude of the signal increases with luminosity, which is tightly correlated with stellar mass in our catalogues. For DESI galaxies between $0<z<1.5$, trends in luminosity and colour alone are sufficient to explain the alignments we measure -- with no need for an explicit redshift dependence. In a companion paper (Jeffrey et al., in prep), we perform detailed modelling of the IA signals with significant detections, including model comparison. Finally, to direct efforts for future IA measurements, we juxtapose the colour-magnitude-redshift coverage of existing IA measurements against modern and future lensing surveys. △ Less

Submitted 16 July, 2025; v1 submitted 15 July, 2025; originally announced July 2025.

Comments: 22 pages, submitted for MNRAS

Euclid preparation. Full-shape modelling of 2-point and 3-point correlation functions in real space

Authors: Euclid Collaboration, M. Guidi, A. Veropalumbo, A. Pugno, M. Moresco, E. Sefusatti, C. Porciani, E. Branchini, M. -A. Breton, B. Camacho Quevedo, M. Crocce, S. de la Torre, V. Desjacques, A. Eggemeier, A. Farina, M. Kärcher, D. Linde, M. Marinucci, A. Moradinezhad Dizgah, C. Moretti, K. Pardede, A. Pezzotta, E. Sarpa, A. Amara, S. Andreon , et al. (286 additional authors not shown)

Abstract: We investigate the accuracy and range of validity of the perturbative model for the 2-point (2PCF) and 3-point (3PCF) correlation functions in real space in view of the forthcoming analysis of the Euclid mission spectroscopic sample. We take advantage of clustering measurements from four snapshots of the Flagship I N-body simulations at z = {0.9, 1.2, 1.5, 1.8}, which mimic the expected galaxy pop… ▽ More We investigate the accuracy and range of validity of the perturbative model for the 2-point (2PCF) and 3-point (3PCF) correlation functions in real space in view of the forthcoming analysis of the Euclid mission spectroscopic sample. We take advantage of clustering measurements from four snapshots of the Flagship I N-body simulations at z = {0.9, 1.2, 1.5, 1.8}, which mimic the expected galaxy population in the ideal case of absence of observational effects such as purity and completeness. For the 3PCF we consider all available triangle configurations given a minimal separation. First, we assess the model performance by fixing the cosmological parameters and evaluating the goodness-of-fit provided by the perturbative bias expansion in the joint analysis of the two statistics, finding overall agreement with the data down to separations of 20 Mpc/h. Subsequently, we build on the state-of-the-art and extend the analysis to include the dependence on three cosmological parameters: the amplitude of scalar perturbations As, the matter density ωcdm and the Hubble parameter h. To achieve this goal, we develop an emulator capable of generating fast and robust modelling predictions for the two summary statistics, allowing efficient sampling of the joint likelihood function. We therefore present the first joint full-shape analysis of the real-space 2PCF and 3PCF, testing the consistency and constraining power of the perturbative model across both probes, and assessing its performance in a combined likelihood framework. We explore possible systematic uncertainties induced by the perturbative model at small scales finding an optimal scale cut of rmin = 30 Mpc/h for the 3PCF, when imposing an additional limitation on nearly isosceles triangular configurations included in the data vector. This work is part of a Euclid Preparation series validating theoretical models for galaxy clustering. △ Less

Submitted 27 June, 2025; originally announced June 2025.

Comments: 21 pages, 11 figures

Journal ref: A&A 707, A228 (2026)

Lensing Without Borders: Measurements of galaxy-galaxy lensing and projected galaxy clustering in DESI DR1

Authors: S. Heydenreich, A. Leauthaud, C. Blake, Z. Sun, J. U. Lange, T. Zhang, M. DeMartino, A. J. Ross, J. Aguilar, S. Ahlen, D. Bianchi, D. Brooks, F. J. Castander, T. Claybaugh, A. Cuceu, A. de la Macorra, J. DeRose, Arjun Dey, Biprateep Dey, P. Doel, N. Emas, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, C. Garcia-Quintero , et al. (42 additional authors not shown)

Abstract: We present Galaxy-Galaxy Lensing measurements obtained by cross-correlating spectroscopically observed galaxies from the first data release of the Dark Energy Spectroscopic Instrument (DESI) with source galaxies from the Hyper Suprime-Cam Subaru Strategic Survey, the Kilo-Degree Survey, the Sloan Digital Sky Survey, and the Dark Energy Survey. Specifically, we measure the excess surface mass densi… ▽ More We present Galaxy-Galaxy Lensing measurements obtained by cross-correlating spectroscopically observed galaxies from the first data release of the Dark Energy Spectroscopic Instrument (DESI) with source galaxies from the Hyper Suprime-Cam Subaru Strategic Survey, the Kilo-Degree Survey, the Sloan Digital Sky Survey, and the Dark Energy Survey. Specifically, we measure the excess surface mass density $ΔΣ$ and tangential shear $γ_\mathrm{t}$ for the Bright Galaxy Sample and Luminous Red Galaxies measured within the first year of observations with DESI. To ensure robustness, we test the measurements for systematic biases, finding no significant trends related to the properties of the \acrshort{desi} lens galaxies. We identify a significant trend with the average redshift of source galaxies, however, this trend vanishes once we apply shifts to the Hyper Suprime-Cam Subaru Strategic Survey redshift distributions that are also favored by their fiducial cosmology analysis. Additionally, we compare the observed scatter in the measurements with the theoretical covariance and find excess scatter, driven primarily by small-scale measurements of $r\leq 1 \, \mathrm{Mpc}/h$; measurements on larger scales are consistent at the $2\,σ$ level. We further present the projected clustering measurements $w_p$ of the galaxy samples in the the first data release of DESI. These measurements, which will be made publicly available, serve as a foundation for forthcoming cosmological analyses. △ Less

Submitted 26 June, 2025; originally announced June 2025.

Euclid preparation. LXXXIX. Accurate and precise data-driven angular power spectrum covariances

Authors: Euclid Collaboration, K. Naidoo, J. Ruiz-Zapatero, N. Tessore, B. Joachimi, A. Loureiro, N. Aghanim, B. Altieri, A. Amara, L. Amendola, S. Andreon, N. Auricchio, C. Baccigalupi, D. Bagot, M. Baldi, S. Bardelli, P. Battaglia, A. Biviano, E. Branchini, M. Brescia, S. Camera, V. Capobianco, C. Carbone, V. F. Cardone, J. Carretero , et al. (258 additional authors not shown)

Abstract: We develop techniques for generating accurate and precise internal covariances for measurements of clustering and weak-lensing angular power spectra. These methods have been designed to produce non-singular and unbiased covariances for Euclid's large anticipated data vector and will be critical for validation against observational systematic effects. We constructed jackknife segments that are equa… ▽ More We develop techniques for generating accurate and precise internal covariances for measurements of clustering and weak-lensing angular power spectra. These methods have been designed to produce non-singular and unbiased covariances for Euclid's large anticipated data vector and will be critical for validation against observational systematic effects. We constructed jackknife segments that are equal in area to a high precision by adapting the binary space partition algorithm to work on arbitrarily shaped regions on the unit sphere. Jackknife estimates of the covariances are internally derived and require no assumptions about cosmology or galaxy population and bias. Our covariance estimation, called DICES (Debiased Internal Covariance Estimation with Shrinkage), first estimated a noisy covariance through conventional delete-1 jackknife resampling. This was followed by linear shrinkage of the empirical correlation matrix towards the Gaussian prediction, rather than linear shrinkage of the covariance matrix. Shrinkage ensures the covariance is non-singular and therefore invertible, which is critical for the estimation of likelihoods and validation. We then applied a delete-2 jackknife bias correction to the diagonal components of the jackknife covariance that removed the general tendency for jackknife error estimates to be biased high. We validated internally derived covariances, which used the jackknife resampling technique, on synthetic Euclid-like lognormal catalogues. We demonstrate that DICES produces accurate, non-singular covariance estimates, with the relative error improving by 33% for the covariance and 48% for the correlation structure in comparison to jackknife estimates. These estimates can be used for highly accurate regression and inference. △ Less

Submitted 13 April, 2026; v1 submitted 10 June, 2025; originally announced June 2025.

Comments: 20 pages, 21 figures, matches published version in A&A. The skysegmentor package developed in this analysis is available here: see http://pypi.org/project/skysegmentor/

Journal ref: A&A 708, A167 (2026)

Euclid preparation. The impact of redshift interlopers on the two-point correlation function analysis

Authors: Euclid Collaboration, I. Risso, A. Veropalumbo, E. Branchini, E. Maragliano, S. de la Torre, E. Sarpa, P. Monaco, B. R. Granett, S. Lee, G. E. Addison, S. Bruton, C. Carbone, G. Lavaux, K. Markovic, K. McCarthy, G. Parimbelli, F. Passalacqua, W. J. Percival, C. Scarlata, E. Sefusatti, Y. Wang, M. Bonici, F. Oppizzi, N. Aghanim , et al. (295 additional authors not shown)

Abstract: The Euclid survey aims to measure the spectroscopic redshift of emission-line galaxies by identifying the H$\,α$ line in their slitless spectra. This method is sensitive to the signal-to-noise ratio of the line, as noise fluctuations or other strong emission lines can be misidentified as H$\,α$, depending on redshift. These effects lead to catastrophic redshift errors and the inclusion of interlop… ▽ More The Euclid survey aims to measure the spectroscopic redshift of emission-line galaxies by identifying the H$\,α$ line in their slitless spectra. This method is sensitive to the signal-to-noise ratio of the line, as noise fluctuations or other strong emission lines can be misidentified as H$\,α$, depending on redshift. These effects lead to catastrophic redshift errors and the inclusion of interlopers in the sample. We forecast the impact of such redshift errors on galaxy clustering measurements. In particular, we study the effect of interloper contamination on the two-point correlation function (2PCF), the growth rate of structures, and the Alcock-Paczynski (AP) parameters. We analyze 1000 synthetic spectroscopic catalogues, the EuclidLargeMocks, designed to match the area and selection function of the Data Release 1 (DR1) sample. We estimate the 2PCF of the contaminated catalogues, isolating contributions from correctly identified galaxies and from interlopers. We explore different models with increasing complexity to describe the measured 2PCF at fixed cosmology. Finally, we perform a cosmological inference and evaluate the systematic error on the inferred $fσ_8$, $α_{\parallel}$ and $α_{\perp}$ values associated with different models. Our results demonstrate that a minimal modelling approach, which only accounts for an attenuation of the clustering signal regardless of the type of contaminants, is sufficient to recover the correct values of $fσ_8$, $α_{\parallel}$, and $α_{\perp}$ at DR1. The accuracy and precision of the estimated AP parameters are largely insensitive to the presence of interlopers. The adoption of a minimal model induces a 1%-3% systematic error on the growth rate of structure estimation, depending on the redshift. However, this error remains smaller than the statistical error expected for the Euclid DR1 analysis. △ Less

Submitted 27 March, 2026; v1 submitted 7 May, 2025; originally announced May 2025.

Comments: 28 pages, published on A&A

Journal ref: A&A 707, A233 (2026)

Redshift Assessment Infrastructure Layers (RAIL): Rubin-era photometric redshift stress-testing and at-scale production

Authors: The RAIL Team, Jan Luca van den Busch, Eric Charles, Johann Cohen-Tanugi, Alice Crafford, John Franklin Crenshaw, Sylvie Dagoret, Josue De-Santiago, Juan De Vicente, Qianjun Hang, Benjamin Joachimi, Shahab Joudaki, J. Bryce Kalmbach, Arun Kannawadi, Shuang Liang, Olivia Lynn, Alex I. Malz, Rachel Mandelbaum, Grant Merz, Irene Moskowitz, Drew Oldag, Jaime Ruiz-Zapatero, Mubdi Rahman, Markus M. Rau, Samuel J. Schmidt , et al. (8 additional authors not shown)

Abstract: Virtually all extragalactic use cases of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require the use of galaxy redshift information, yet the vast majority of its sample of tens of billions of galaxies will lack high-fidelity spectroscopic measurements thereof, instead relying on photometric redshifts (photo-$z$) subject to systematic imprecision and inaccuracy best encap… ▽ More Virtually all extragalactic use cases of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require the use of galaxy redshift information, yet the vast majority of its sample of tens of billions of galaxies will lack high-fidelity spectroscopic measurements thereof, instead relying on photometric redshifts (photo-$z$) subject to systematic imprecision and inaccuracy best encapsulated by photo-$z$ probability density functions (PDFs). We present the version 1 release of Redshift Assessment Infrastructure Layers (RAIL), an open source Python library for at-scale probabilistic photo-$z$ estimation, initiated by the LSST Dark Energy Science Collaboration (DESC) with contributions from the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) Frameworks team. RAIL's three subpackages provide modular tools for end-to-end stress-testing, including a forward modeling suite to generate realistically complex photometry, a unified API for estimating per-galaxy and ensemble redshift PDFs by an extensible set of algorithms, and built-in metrics of both photo-$z$ PDFs and point estimates. RAIL serves as a flexible toolkit enabling the derivation and optimization of photo-$z$ data products at scale for a variety of science goals and is not specific to LSST data. We thus describe to the extragalactic science community, including and beyond Rubin the design and functionality of the RAIL software library so that any researcher may have access to its wide array of photo-$z$ characterization and assessment tools. △ Less

Submitted 12 March, 2026; v1 submitted 5 May, 2025; originally announced May 2025.

Comments: Accepted by OJAp, 21 pages, 6 figures, 5 tables

The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics

Authors: Eleonora Di Valentino, Jackson Levi Said, Adam Riess, Agnieszka Pollo, Vivian Poulin, Adrià Gómez-Valent, Amanda Weltman, Antonella Palmese, Caroline D. Huang, Carsten van de Bruck, Chandra Shekhar Saraf, Cheng-Yu Kuo, Cora Uhlemann, Daniela Grandón, Dante Paz, Dominique Eckert, Elsa M. Teixeira, Emmanuel N. Saridakis, Eoin Ó Colgáin, Florian Beutler, Florian Niedermann, Francesco Bajardi, Gabriela Barenboim, Giulia Gubitosi, Ilaria Musella , et al. (516 additional authors not shown)

Abstract: The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-t… ▽ More The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-time inflationary physics, and underpinned by gravitation through general relativity. There have always been open questions about the soundness of the foundations of the standard model. However, recent years have shown that there may also be questions from the observational sector with the emergence of differences between certain cosmological probes. In this White Paper, we identify the key objectives that need to be addressed over the coming decade together with the core science projects that aim to meet these challenges. These discordances primarily rest on the divergence in the measurement of core cosmological parameters with varying levels of statistical confidence. These possible statistical tensions may be partially accounted for by systematics in various measurements or cosmological probes but there is also a growing indication of potential new physics beyond the standard model. After reviewing the principal probes used in the measurement of cosmological parameters, as well as potential systematics, we discuss the most promising array of potential new physics that may be observable in upcoming surveys. We also discuss the growing set of novel data analysis approaches that go beyond traditional methods to test physical models. [Abridged] △ Less

Submitted 4 August, 2025; v1 submitted 2 April, 2025; originally announced April 2025.

Comments: 416 pages, 81 figures

Journal ref: Phys. Dark Univ. 49 (2025) 101965

KiDS-Legacy: Consistency of cosmic shear measurements and joint cosmological constraints with external probes

Authors: Benjamin Stölzner, Angus H. Wright, Marika Asgari, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Benjamin Joachimi, Konrad Kuijken, Shun-Sheng Li, Constance Mahony, Robert Reischke, Mijin Yoon, Maciej Bilicki, Pierre Burger, Nora Elisa Chisari, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Joachim Harnois-Déraps, Priyanka Jalan, Anjitha John William, Shahab Joudaki, Giorgio Francesco Lesci, Laila Linke, Arthur Loureiro , et al. (11 additional authors not shown)

Abstract: We present a cosmic shear consistency analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy). By adopting three tiers of consistency metrics, we compare cosmological constraints between subsets of the KiDS-Legacy dataset split by redshift, angular scale, galaxy colour and spatial region. We also review a range of two-point cosmic shear statistics. With the data passing all ou… ▽ More We present a cosmic shear consistency analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy). By adopting three tiers of consistency metrics, we compare cosmological constraints between subsets of the KiDS-Legacy dataset split by redshift, angular scale, galaxy colour and spatial region. We also review a range of two-point cosmic shear statistics. With the data passing all our consistency metric tests, we demonstrate that KiDS-Legacy is the most internally consistent KiDS catalogue to date. In a joint cosmological analysis of KiDS-Legacy and DES Y3 cosmic shear, combined with data from the Pantheon+ Type Ia supernovae compilation and baryon acoustic oscillations from DESI Y1, we find constraints consistent with Planck measurements of the cosmic microwave background with $S_8\equiv σ_8\sqrt{Ω_{\rm m}/0.3} = 0.814^{+0.011}_{-0.012}$ and $σ_8 = 0.802^{+0.022}_{-0.018}$. △ Less

Submitted 20 October, 2025; v1 submitted 25 March, 2025; originally announced March 2025.

Comments: 25 pages, 16 figures, 7 tables, published in A&A

Journal ref: A&A 702, A169 (2025)

KiDS-Legacy: Cosmological constraints from cosmic shear with the complete Kilo-Degree Survey

Authors: Angus H. Wright, Benjamin Stölzner, Marika Asgari, Maciej Bilicki, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Benjamin Joachimi, Konrad Kuijken, Shun-Sheng Li, Robert Reischke, Maximilian von Wietersheim-Kramsta, Mijin Yoon, Pierre Burger, Nora Elisa Chisari, Jelte de Jong, Andrej Dvornik, Christos Georgiou, Joachim Harnois-Déraps, Priyanka Jalan, Anjitha John William, Shahab Joudaki, Giorgio Francesco Lesci, Laila Linke , et al. (13 additional authors not shown)

Abstract: We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-in… ▽ More We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-infrared, along with an extra $23$ square degrees of KiDS-like calibration observations of deep spectroscopic surveys. Improvements in our redshift distribution estimation methodology, combined with our enhanced calibration data and multi-band image simulations, allow us to extend our lensed sample out to a photometric redshift of $z_{\rm B}\leq2.0$. Compared to previous KiDS analyses, the increased survey area and redshift depth results in a $\sim32\%$ improvement in constraining power in terms of $Σ_8\equivσ_8\left(Ω_{\rm m}/0.3\right)^α= 0.821^{+0.014}_{-0.016}$, where $α= 0.58$ has been optimised to match the revised degeneracy direction of $σ_8$ and $Ω_{\rm m}$. We adopt a new physically motivated intrinsic alignment model that depends jointly on the galaxy sample's halo mass and spectral type distributions, and that is informed by previous direct alignment measurements. We also marginalise over our uncertainty on the impact of baryon feedback on the non-linear matter power spectrum. Comparing to previous KiDS analyses, we conclude that the increase seen in $S_8$ primarily results from our improved redshift distribution estimation and calibration, as well as new survey area and improved image reduction. Our companion paper Stölzner et al. (submitted) presents a full suite of internal and external consistency tests, finding the KiDS-Legacy data set to be the most internally robust sample produced by KiDS to date. △ Less

Submitted 21 October, 2025; v1 submitted 25 March, 2025; originally announced March 2025.

Comments: 38 pages, 24 figures, 7 tables, accepted for publication in A&A

Journal ref: A&A 703, A158 (2025)

KiDS-Legacy: Redshift distributions and their calibration

Authors: Angus H. Wright, Hendrik Hildebrandt, Jan Luca van den Busch, Maciej Bilicki, Catherine Heymans, Benjamin Joachimi, Constance Mahony, Robert Reischke, Benjamin Stölzner, Anna Wittje, Marika Asgari, Nora Elisa Chisari, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Henk Hoekstra, Priyanka Jalan, Anjitha John William, Shahab Joudaki, Konrad Kuijken, Giorgio Francesco Lesci, Shun-Sheng Li, Laila Linke, Arthur Loureiro, Matteo Maturi , et al. (8 additional authors not shown)

Abstract: We present the redshift calibration methodology and bias estimates for the cosmic shear analysis of the fifth and final data release (DR5) of the Kilo-Degree Survey (KiDS). KiDS-DR5 includes a greatly expanded compilation of calibrating spectra, drawn from $27$ square degrees of dedicated optical and near-IR imaging taken over deep spectroscopic fields. The redshift distribution calibration levera… ▽ More We present the redshift calibration methodology and bias estimates for the cosmic shear analysis of the fifth and final data release (DR5) of the Kilo-Degree Survey (KiDS). KiDS-DR5 includes a greatly expanded compilation of calibrating spectra, drawn from $27$ square degrees of dedicated optical and near-IR imaging taken over deep spectroscopic fields. The redshift distribution calibration leverages a range of new methods and updated simulations to produce the most precise $N(z)$ bias estimates used by KiDS to date. Improvements to our colour-based redshift distribution measurement method (SOM) mean that we are able to use many more sources per tomographic bin for our cosmological analyses, and better estimate the representation of our source sample given the available spec-$z$. We validate our colour-based redshift distribution estimates with spectroscopic cross-correlations (CC). We find that improvements to our cross-correlation redshift distribution measurement methods mean that redshift distribution biases estimated between the SOM and CC methods are fully consistent on simulations, and the data calibration is consistent to better than $2σ$ in all tomographic bins. △ Less

Submitted 21 October, 2025; v1 submitted 25 March, 2025; originally announced March 2025.

Comments: 26 pages, 18 figures, 8 tables, accepted for publication in A&A

Journal ref: A&A 703, A144 (2025)

The fifth data release of the Kilo Degree Survey: Multi-epoch optical/NIR imaging covering wide and legacy-calibration fields

Authors: Angus H. Wright, Konrad Kuijken, Hendrik Hildebrandt, Mario Radovich, Maciej Bilicki, Andrej Dvornik, Fedor Getman, Catherine Heymans, Henk Hoekstra, Shun-Sheng Li, Lance Miller, Nicola R. Napolitano, Qianli Xia, Marika Asgari, Massimo Brescia, Hugo Buddelmeijer, Pierre Burger, Gianluca Castignani, Stefano Cavuoti, Jelte de Jong, Alastair Edge, Benjamin Giblin, Carlo Giocoli, Joachim Harnois-Déraps, Priyanka Jalan , et al. (29 additional authors not shown)

Abstract: We present the final data release of the Kilo-Degree Survey (KiDS-DR5), a public European Southern Observatory (ESO) wide-field imaging survey optimised for weak gravitational lensing studies. We combined matched-depth multi-wavelength observations from the VLT Survey Telescope and the VISTA Kilo-degree INfrared Galaxy (VIKING) survey to create a nine-band optical-to-near-infrared survey spanning… ▽ More We present the final data release of the Kilo-Degree Survey (KiDS-DR5), a public European Southern Observatory (ESO) wide-field imaging survey optimised for weak gravitational lensing studies. We combined matched-depth multi-wavelength observations from the VLT Survey Telescope and the VISTA Kilo-degree INfrared Galaxy (VIKING) survey to create a nine-band optical-to-near-infrared survey spanning $1347$ deg$^2$. The median $r$-band $5σ$ limiting magnitude is 24.8 with median seeing $0.7^{\prime\prime}$. The main survey footprint includes $4$ deg$^2$ of overlap with existing deep spectroscopic surveys. We complemented these data in DR5 with a targeted campaign to secure an additional $23$ deg$^2$ of KiDS- and VIKING-like imaging over a range of additional deep spectroscopic survey fields. From these fields, we extracted a catalogue of $126\,085$ sources with both spectroscopic and photometric redshift information, which enables the robust calibration of photometric redshifts across the full survey footprint. In comparison to previous releases, DR5 represents a $34\%$ areal extension and includes an $i$-band re-observation of the full footprint, thereby increasing the effective $i$-band depth by $0.4$ magnitudes and enabling multi-epoch science. Our processed nine-band imaging, single- and multi-band catalogues with masks, and homogenised photometry and photometric redshifts can be accessed through the ESO Archive Science Portal. △ Less

Submitted 25 March, 2025; originally announced March 2025.

Comments: 60 pages, 43 figures, 20 tables, published in A&A

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

Euclid Quick Data Release (Q1). Galaxy shapes and alignments in the cosmic web

Authors: Euclid Collaboration, C. Laigle, C. Gouin, F. Sarron, L. Quilley, C. Pichon, K. Kraljic, F. Durret, N. E. Chisari, U. Kuchner, N. Malavasi, M. Magliocchetti, H. J. McCracken, J. G. Sorce, Y. Kang, C. J. R. McPartland, S. Toft, N. Aghanim, B. Altieri, A. Amara, S. Andreon, N. Auricchio, H. Aussel, C. Baccigalupi, M. Baldi , et al. (319 additional authors not shown)

Abstract: Galaxy morphologies and shape orientations are expected to correlate with their large-scale environment, since they grow by accreting matter from the cosmic web and are subject to interactions with other galaxies. Cosmic filaments are extracted in projection from the Euclid Quick Data Release 1 (covering 63.1 $\mathrm{deg}^2$) at $0.5<z<0.9$ in tomographic slices of 170 comoving… ▽ More Galaxy morphologies and shape orientations are expected to correlate with their large-scale environment, since they grow by accreting matter from the cosmic web and are subject to interactions with other galaxies. Cosmic filaments are extracted in projection from the Euclid Quick Data Release 1 (covering 63.1 $\mathrm{deg}^2$) at $0.5<z<0.9$ in tomographic slices of 170 comoving $h^{-1}\mathrm{Mpc}$ using photometric redshifts. Galaxy morphologies are accurately retrieved thanks to the excellent resolution of VIS data. The distribution of massive galaxies ($M_* > 10^{10} M_\odot$) in the projected cosmic web is analysed as a function of morphology measured from VIS data. Specifically, the 2D alignment of galaxy shapes with large-scale filaments is quantified as a function of Sérsic indices and masses. We find the known trend that more massive galaxies are closer to filament spines. At fixed stellar masses, morphologies correlate both with densities and distances to large-scale filaments. In addition, the large volume of this data set allows us to detect a signal indicating that there is a preferential alignment of the major axis of massive early-type galaxies along projected cosmic filaments. Overall, these results demonstrate our capabilities to carry out detailed studies of galaxy environments with Euclid, which will be extended to higher redshift and lower stellar masses with the future Euclid Deep Survey. △ Less

Submitted 19 March, 2025; originally announced March 2025.

Comments: Paper submitted as part of the A&A Special Issue `Euclid Quick Data Release (Q1)', 22 pages, 19 figures

Euclid Quick Data Release (Q1) -- Data release overview

Authors: Euclid Collaboration, H. Aussel, I. Tereno, M. Schirmer, G. Alguero, B. Altieri, E. Balbinot, T. de Boer, P. Casenove, P. Corcho-Caballero, H. Furusawa, J. Furusawa, M. J. Hudson, K. Jahnke, G. Libet, J. Macias-Perez, N. Masoumzadeh, J. J. Mohr, J. Odier, D. Scott, T. Vassallo, G. Verdoes Kleijn, A. Zacchei, N. Aghanim, A. Amara , et al. (385 additional authors not shown)

Abstract: The first Euclid Quick Data Release, Q1, comprises 63.1 sq deg of the Euclid Deep Fields (EDFs) to nominal wide-survey depth. It encompasses visible and near-infrared space-based imaging and spectroscopic data, ground-based photometry in the u, g, r, i and z bands, as well as corresponding masks. Overall, Q1 contains about 30 million objects in three areas near the ecliptic poles around the EDF-No… ▽ More The first Euclid Quick Data Release, Q1, comprises 63.1 sq deg of the Euclid Deep Fields (EDFs) to nominal wide-survey depth. It encompasses visible and near-infrared space-based imaging and spectroscopic data, ground-based photometry in the u, g, r, i and z bands, as well as corresponding masks. Overall, Q1 contains about 30 million objects in three areas near the ecliptic poles around the EDF-North and EDF-South, as well as the EDF-Fornax field in the constellation of the same name. The purpose of this data release -- and its associated technical papers -- is twofold. First, it is meant to inform the community of the enormous potential of the Euclid survey data, to describe what is contained in these data, and to help prepare expectations for the forthcoming first major data release DR1. Second, it enables a wide range of initial scientific projects with wide-survey Euclid data, ranging from the early Universe to the Solar System. The Q1 data were processed with early versions of the processing pipelines, which already demonstrate good performance, with numerous improvements in implementation compared to pre-launch development. In this paper, we describe the sky areas released in Q1, the observations, a top-level view of the data processing of Euclid and associated external data, the Q1 photometric masks, and how to access the data. We also give an overview of initial scientific results obtained using the Q1 data set by Euclid Consortium scientists, and conclude with important caveats when using the data. As a complementary product, Q1 also contains observations of a star-forming area in Lynd's Dark Nebula 1641 in the Orion~A Cloud, observed for technical purposes during Euclid's performance-verification phase. This is a unique target, of a type not commonly found in Euclid's nominal sky survey. △ Less

Submitted 19 March, 2025; originally announced March 2025.

Comments: 27 pages, 12 figures, data release at https://www.cosmos.esa.int/en/web/euclid/euclid-q1-data-release paper submitted to the special A&A issue

Data Release 1 of the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, M. Abdul Karim, A. G. Adame, D. Aguado, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, L. Allen, C. Allende Prieto, O. Alves, A. Anand, U. Andrade, E. Armengaud, S. Avila, A. Aviles, H. Awan, S. Bailey, A. Baleato Lizancos, O. Ballester, A. Bault, J. Bautista, R. Bean , et al. (285 additional authors not shown)

Abstract: In 2021 May the Dark Energy Spectroscopic Instrument (DESI) collaboration began a 5-year spectroscopic redshift survey to produce a detailed map of the evolving three-dimensional structure of the universe between $z=0$ and $z\approx4$. DESI's principle scientific objectives are to place precise constraints on the equation of state of dark energy, the gravitationally driven growth of large-scale st… ▽ More In 2021 May the Dark Energy Spectroscopic Instrument (DESI) collaboration began a 5-year spectroscopic redshift survey to produce a detailed map of the evolving three-dimensional structure of the universe between $z=0$ and $z\approx4$. DESI's principle scientific objectives are to place precise constraints on the equation of state of dark energy, the gravitationally driven growth of large-scale structure, and the sum of the neutrino masses, and to explore the observational signatures of primordial inflation. We present DESI Data Release 1 (DR1), which consists of all data acquired during the first 13 months of the DESI main survey, as well as a uniform reprocessing of the DESI Survey Validation data which was previously made public in the DESI Early Data Release. The DR1 main survey includes high-confidence redshifts for 18.7M objects, of which 13.1M are spectroscopically classified as galaxies, 1.6M as quasars, and 4M as stars, making DR1 the largest sample of extragalactic redshifts ever assembled. We summarize the DR1 observations, the spectroscopic data-reduction pipeline and data products, large-scale structure catalogs, value-added catalogs, and describe how to access and interact with the data. In addition to fulfilling its core cosmological objectives with unprecedented precision, we expect DR1 to enable a wide range of transformational astrophysical studies and discoveries. △ Less

Submitted 4 March, 2026; v1 submitted 18 March, 2025; originally announced March 2025.

Comments: 64 pages, 7 figures, 15 tables, accepted to The Astronomical Journal

Euclid preparation. 3-dimensional galaxy clustering in configuration space. Part I. 2-point correlation function estimation

Authors: Euclid Collaboration, S. de la Torre, F. Marulli, E. Keihänen, A. Viitanen, M. Viel, A. Veropalumbo, E. Branchini, D. Tavagnacco, F. Rizzo, J. Valiviita, V. Lindholm, V. Allevato, G. Parimbelli, E. Sarpa, Z. Ghaffari, A. Amara, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, S. Bardelli, A. Basset, D. Bonino, M. Brescia , et al. (275 additional authors not shown)

Abstract: The 2-point correlation function of the galaxy spatial distribution is a major cosmological observable that enables constraints on the dynamics and geometry of the Universe. The Euclid mission aims at performing an extensive spectroscopic survey of approximately 20--30 million H$α$-emitting galaxies up to about redshift two. This ambitious project seeks to elucidate the nature of dark energy by ma… ▽ More The 2-point correlation function of the galaxy spatial distribution is a major cosmological observable that enables constraints on the dynamics and geometry of the Universe. The Euclid mission aims at performing an extensive spectroscopic survey of approximately 20--30 million H$α$-emitting galaxies up to about redshift two. This ambitious project seeks to elucidate the nature of dark energy by mapping the 3-dimensional clustering of galaxies over a significant portion of the sky. This paper presents the methodology and software developed for estimating the 3-dimensional 2-point correlation function within the Euclid Science Ground Segment. The software is designed to overcome the significant challenges posed by the large and complex Euclid data set, which involves millions of galaxies. Key challenges include efficient pair counting, managing computational resources, and ensuring the accuracy of the correlation function estimation. The software leverages advanced algorithms, including kd-tree, octree, and linked-list data partitioning strategies, to optimise the pair-counting process. The implementation also includes parallel processing capabilities using shared-memory open multi-processing to further enhance performance and reduce computation times. Extensive validation and performance testing of the software are presented. The results indicate that the software is robust and can reliably estimate the 2-point correlation function, which is essential for deriving cosmological parameters with high precision. Furthermore, the paper discusses the expected performance of the software during different stages of the Euclid Wide Survey observations and forecasts how the precision of the correlation function measurements will improve over the mission's timeline, highlighting the software's capability to handle large data sets efficiently. △ Less

Submitted 27 January, 2025; originally announced January 2025.

Comments: 17 pages, 13 figures, submitted to A&A

Journal ref: A&A 700, A78 (2025)

The DESI-Lensing Mock Challenge: large-scale cosmological analysis of 3x2-pt statistics

Authors: C. Blake, C. Garcia-Quintero, S. Ahlen, D. Bianchi, D. Brooks, T. Claybaugh, A. de la Macorra, J. DeRose, A. Dey, P. Doel, N. Emas, S. Ferraro, J. E. Forero-Romero, G. Gutierrez, S. Heydenreich, K. Honscheid, C. Howlett, M. Ishak, S. Joudaki, E. Jullo, R. Kehoe, D. Kirkby, A. Kremin, A. Krolewski, M. Landriau , et al. (19 additional authors not shown)

Abstract: The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study… ▽ More The current generation of large galaxy surveys will test the cosmological model by combining multiple types of observational probes. Realising the statistical promise of these new datasets requires rigorous attention to all aspects of analysis including cosmological measurements, modelling, covariance and parameter likelihood. In this paper we present the results of an end-to-end simulation study designed to test the analysis pipeline for the combination of the Dark Energy Spectroscopic Instrument (DESI) Year 1 galaxy redshift dataset and separate weak gravitational lensing information from the Kilo-Degree Survey, Dark Energy Survey and Hyper-Suprime-Cam Survey. Our analysis employs the 3x2-pt correlation functions including cosmic shear and galaxy-galaxy lensing, together with the projected correlation function of the spectroscopic DESI lenses. We build realistic simulations of these datasets including galaxy halo occupation distributions, photometric redshift errors, weights, multiplicative shear calibration biases and magnification. We calculate the analytical covariance of these correlation functions including the Gaussian, noise and super-sample contributions, and show that our covariance determination agrees with estimates based on the ensemble of simulations. We use a Bayesian inference platform to demonstrate that we can recover the fiducial cosmological parameters of the simulation within the statistical error margin of the experiment, investigating the sensitivity to scale cuts. This study is the first in a sequence of papers in which we present and validate the large-scale 3x2-pt cosmological analysis of DESI-Y1. △ Less

Submitted 6 March, 2025; v1 submitted 17 December, 2024; originally announced December 2024.

Comments: 29 pages, 17 figures, accepted for publication in The Open Journal of Astrophysics

Super sample covariance and the volume scaling of galaxy survey covariance matrices

Authors: Greg Schreiner, Alex Krolewski, Shahab Joudaki, Will J. Percival

Abstract: Super sample covariance (SSC) is important when estimating covariance matrices using a set of mock catalogues for galaxy surveys. If the underlying cosmological simulations do not include the variation in background parameters appropriate for the simulation sizes, then the scatter between mocks will be missing the SSC component. The coupling between large and small modes due to non-linear structur… ▽ More Super sample covariance (SSC) is important when estimating covariance matrices using a set of mock catalogues for galaxy surveys. If the underlying cosmological simulations do not include the variation in background parameters appropriate for the simulation sizes, then the scatter between mocks will be missing the SSC component. The coupling between large and small modes due to non-linear structure growth makes this pernicious on small scales. We compare different methods for generating ensembles of mocks with SSC built in to the covariance, and contrast against methods where the SSC component is computed and added to the covariance separately. We find that several perturbative expansions, developed to derive background fluctuations, give similar results. We then consider scaling covariance matrices calculated for simulations of different volumes to improve the accuracy of the covariance matrix for a given computational time. On large scales, we find that the primary limitation is from the discrete number of modes contributing to the measured power spectrum, and we propose a new method for correcting this effect. Correct implementation of SSC and the effect of discrete mode numbers allows covariance matrices created from mocks to be scaled between volumes, potentially leading to a significant saving on computational resources when producing covariance matrices. We argue that a sub-percent match is difficult to achieve because of the effects of modes on scales between the box sizes, which cannot be easily included. Even so, when working in real space and cubic boxes, we show that a 3% match in the dark matter power spectrum covariance is achievable on scales of interest for current surveys scaling the simulation volume by 512x, costing a small fraction of the computational time of running full-sized simulations. △ Less

Submitted 3 March, 2025; v1 submitted 25 November, 2024; originally announced November 2024.

Comments: 26 pages, 9 figures. Published in JCAP

Journal ref: JCAP02(2025)022

KiDS-Legacy: Angular galaxy clustering from deep surveys with complex selection effects

Authors: Ziang Yan, Angus H. Wright, Nora Elisa Chisari, Christos Georgiou, Shahab Joudaki, Arthur Loureiro, Robert Reischke, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Catherine Heymans, Hendrik Hildebrandt, Priyanka Jalan, Benjamin Joachimi, Giorgio Francesco Lesci, Shun-Sheng Li, Laila Linke, Constance Mahony, Lauro Moscardini, Nicola R. Napolitano, Benjamin Stoelzner, Maximilian Von Wietersheim-Kramsta, Mijin Yoon

Abstract: Photometric galaxy surveys, despite their limited resolution along the line of sight, encode rich information about the large-scale structure (LSS) of the Universe thanks to the high number density and extensive depth of the data. However, the complicated selection effects in wide and deep surveys can potentially cause significant bias in the angular two-point correlation function (2PCF) measured… ▽ More Photometric galaxy surveys, despite their limited resolution along the line of sight, encode rich information about the large-scale structure (LSS) of the Universe thanks to the high number density and extensive depth of the data. However, the complicated selection effects in wide and deep surveys can potentially cause significant bias in the angular two-point correlation function (2PCF) measured from those surveys. In this paper, we measure the 2PCF from the newly published KiDS-Legacy sample. Given an $r$-band $5σ$ magnitude limit of $24.8$ and survey footprint of $1347$ deg$^2$, it achieves an excellent combination of sky coverage and depth for such a measurement. We find that complex selection effects, primarily induced by varying seeing, introduce over-estimation of the 2PCF by approximately an order of magnitude. To correct for such effects, we apply a machine learning-based method to recover an organised random (OR) that presents the same selection pattern as the galaxy sample. The basic idea is to find the selection-induced clustering of galaxies using a combination of self-organising maps (SOMs) and hierarchical clustering (HC). This unsupervised machine learning method is able to recover complicated selection effects without specifying their functional forms. We validate this SOM+HC method on mock deep galaxy samples with realistic systematics and selections derived from the KiDS-Legacy catalogue. Using mock data, we demonstrate that the OR delivers unbiased 2PCF cosmological parameter constraints, removing the $27σ$ offset in the galaxy bias parameter that is recovered when adopting uniform randoms. Blinded measurements on the real KiDS-Legacy data show that the corrected 2PCF is robust to the SOM+HC configuration near the optimal set-up suggested by the mock tests. △ Less

Submitted 3 February, 2025; v1 submitted 30 October, 2024; originally announced October 2024.

Comments: 29 pages, 27 figures, 4 tables; Accepted for publication on Astronomy & Astrophysics; The code used for this work is published on https://github.com/yanzastro/tiaogeng

Journal ref: A&A 694, A259 (2025)

KiDS-Legacy: Covariance validation and the unified OneCovariance framework for projected large-scale structure observables

Authors: Robert Reischke, Sandra Unruh, Marika Asgari, Andrej Dvornik, Hendrik Hildebrandt, Benjamin Joachimi, Lucas Porth, Maximilian von Wietersheim-Kramsta, Jan Luca van den Busch, Benjamin Stölzner, Angus H. Wright, Ziang Yan, Maciej Bilicki, Pierre Burger, Nora Elisa Chisari, Joachim Harnois-Deraps, Christos Georgiou, Catherine Heymans, Priyanka Jalan, Shahab Joudaki, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Constance Mahony, Davide Sciotti , et al. (2 additional authors not shown)

Abstract: We introduce OneCovariance, an open-source software designed to accurately compute covariance matrices for an arbitrary set of two-point summary statistics across a variety of large-scale structure tracers. Utilising the halo model, we estimated the statistical properties of matter and biased tracer fields, incorporating all Gaussian, non-Gaussian, and super-sample covariance terms. The flexible c… ▽ More We introduce OneCovariance, an open-source software designed to accurately compute covariance matrices for an arbitrary set of two-point summary statistics across a variety of large-scale structure tracers. Utilising the halo model, we estimated the statistical properties of matter and biased tracer fields, incorporating all Gaussian, non-Gaussian, and super-sample covariance terms. The flexible configuration permits user-specific parameters, such as the complexity of survey geometry, the halo occupation distribution employed to define each galaxy sample, or the form of the real-space and/or Fourier space statistics to be analysed. We illustrate the capabilities of OneCovariance within the context of a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy). Upon comparing our estimated covariance with measurements from mock data and calculations from independent software, we ascertain that OneCovariance achieves accuracy at the per cent level. When assessing the impact of ignoring complex survey geometry in the cosmic shear covariance computation, we discover misestimations at approximately the $10\%$ level for cosmic variance terms. Nonetheless, these discrepancies do not significantly affect the KiDS-Legacy recovery of cosmological parameters. We derive the cross-covariance between real-space correlation functions, bandpowers, and COSEBIs, facilitating future consistency tests among these three cosmic shear statistics. Additionally, we calculate the covariance matrix of photometric-spectroscopic galaxy clustering measurements, validating the jackknife covariance estimates for calibrating KiDS-Legacy redshift distributions. The OneCovariance can be found on GitHub, together with comprehensive documentation and examples. △ Less

Submitted 16 September, 2025; v1 submitted 9 October, 2024; originally announced October 2024.

Comments: 37 pages, 11 figures, published in A&A, code available at https://github.com/rreischke/OneCovariance

Journal ref: A&A 699, A124 (2025)

6x2pt: Forecasting gains from joint weak lensing and galaxy clustering analyses with spectroscopic-photometric galaxy cross-correlations

Authors: Harry Johnston, Nora Elisa Chisari, Shahab Joudaki, Robert Reischke, Benjamin Stölzner, Arthur Loureiro, Constance Mahony, Sandra Unruh, Angus H. Wright, Marika Asgari, Maciej Bilicki, Pierre Burger, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Lucas Porth, HuanYuan Shan, Tilman Tröster, Jan Luca van den Busch , et al. (3 additional authors not shown)

Abstract: We explore the enhanced self-calibration of photometric galaxy redshift distributions, $n(z)$, through the combination of up to six two-point functions. Our $\rm 3\times2pt$ configuration is comprised of photometric shear, spectroscopic galaxy clustering, and spectroscopic-photometric galaxy-galaxy lensing (GGL). We further include spectroscopic-photometric cross-clustering; photometric GGL; and p… ▽ More We explore the enhanced self-calibration of photometric galaxy redshift distributions, $n(z)$, through the combination of up to six two-point functions. Our $\rm 3\times2pt$ configuration is comprised of photometric shear, spectroscopic galaxy clustering, and spectroscopic-photometric galaxy-galaxy lensing (GGL). We further include spectroscopic-photometric cross-clustering; photometric GGL; and photometric auto-clustering, using the photometric shear sample as density tracer. We perform simulated likelihood forecasts of the cosmological and nuisance parameter constraints for Stage-III- and Stage-IV-like surveys. For the Stage-III-like case, we employ realistic but perturbed redshift distributions, and distinguish between "coherent" shifting in one direction, versus more internal scattering and full-shape errors. For perfectly known $n(z)$, a $\rm 6\times2pt$ analysis gains $\sim40\%$ in Figure of Merit (FoM) in the $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3}$ and $Ω_{\rm m}$ plane relative to the $\rm 3\times2pt$ analysis. If untreated, coherent and incoherent redshift errors lead to inaccurate inferences of $S_8$ and $Ω_{\rm m}$, respectively. Employing bin-wise scalar shifts $δ{z}_i$ in the tomographic mean redshifts reduces cosmological parameter biases, with a $\rm 6x2pt$ analysis constraining the shift parameters with $2-4$ times the precision of a photometric $\rm 3^{ph}\times2pt$ analysis. For the Stage-IV-like survey, a $\rm 6\times2pt$ analysis doubles the FoM($σ_8{-}Ω_{\rm m}$) compared to any $\rm 3\times2pt$ or $\rm 3^{ph}\times2pt$ analysis, and is only $8\%$ less constraining than if the $n(z)$ were perfectly known. A Gaussian mixture model for the $n(z)$ reduces mean-redshift errors and preserves the $n(z)$ shape. It also yields the most accurate and precise cosmological constraints for any $N\rm\times2pt$ configuration given $n(z)$ biases. △ Less

Submitted 25 September, 2024; originally announced September 2024.

Comments: 38 pages, 20 figures, to be submitted to A&A

Journal ref: A&A 699, A127 (2025)