Showing 1–50 of 90 results for author: Pourtsidou, A

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

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

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

Submitted 22 May, 2024; originally announced May 2024.

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

Euclid. III. The NISP Instrument

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

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

Submitted 22 May, 2024; originally announced May 2024.

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

Euclid. II. The VIS Instrument

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

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

Submitted 22 May, 2024; originally announced May 2024.

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

Euclid. I. Overview of the Euclid mission

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

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

Submitted 22 May, 2024; originally announced May 2024.

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

Power spectrum multipoles and clustering wedges during the Epoch of Reionization

Authors: Zhaoting Chen, Alkistis Pourtsidou

Abstract: We study the viability of using power spectrum clustering wedges as summary statistics of 21cm surveys during the Epoch of Reionization (EoR). For observations along a large lightcone $z\sim 7-9$, the power spectrum is subject to large anisotropic effects due to the evolution along the light-of-sight. Information on the physics of reionization can be extracted from the anisotropy using the power s… ▽ More We study the viability of using power spectrum clustering wedges as summary statistics of 21cm surveys during the Epoch of Reionization (EoR). For observations along a large lightcone $z\sim 7-9$, the power spectrum is subject to large anisotropic effects due to the evolution along the light-of-sight. Information on the physics of reionization can be extracted from the anisotropy using the power spectrum multipoles. Signals of the power spectrum monopole are highly correlated at scales smaller than the typical ionization bubble, which can be disentangled by including higher-order multipoles. By simulating observations of the low frequency part of the SKA Observatory, we find that the sampling of the cylindrical wavenumber $k$-space is highly non-uniform due to the baseline distribution. Measurements in clustering wedges can be used for isolating parts of $k$-space with relatively uniform sampling, allowing for more precise parameter inference. Using Fisher Matrix forecasts, we find that the reionization model can be inferred with per-cent level precision with $\sim 120$hrs of integration time using SKA-Low. Comparing to using only the power spectrum monopole above the foreground wedge, model inference using multipole power spectra in clustering wedges yields a factor of $\sim 3$ improvement on parameter constraints. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 14 pages, 11 figures. Submitted to MNRAS

Radio Frequency Interference from Radio Navigation Satellite Systems: simulations and comparison to MeerKAT single-dish data

Authors: Brandon Engelbrecht, Mario G. Santos, José Fonseca, Yichao Li, Jingying Wang, Melis O. Irfan, Stuart E. Harper, Keith Grainge, Philip Bull, Isabella P. Carucci, Steven Cunnington, Alkistis Pourtsidou, Marta Spinelli, Laura Wolz

Abstract: Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital source… ▽ More Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital sources in the Radio Navigational Satellite System within the 1100-1350 MHz frequency. This simulation can be split into two parts: (I) satellite positioning, emission power, and beam response on the telescope and (II) fitting of the satellite signal to data in order to improve the original model. We use previously observed single dish MeerKAT L-band data which needs to be specially calibrated to include data contaminated by satellite-based RFI. We find that due to non-linearity effects, it becomes non-trivial to fit the satellite power. However, when masking regions where this non-linearity is problematic, we can recreate the satellite contamination with high accuracy around its peak frequencies. The simulation can predict satellite movements and signal for past and future observations, which can help in RFI avoidance and testing novel cleaning methods. The predicted signal from simulations sits below the noise in the target cosmology window for the L-band (970 - 1015 MHz) making it difficult to confirm any out-of-band emission from satellites. However, a power spectrum analysis shows that such signal can still contaminate the 21cm power spectrum at these frequencies. In our simulations, this contamination overwhelms the auto-power spectrum but still allows for a clean detection of the signal in cross-correlations with mild foreground cleaning. Whether such contamination does exist one will require further characterization of the satellite signals far away from their peak frequencies. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: 20 pages, 27 figures, 4 tables, Comments welcome

Stage-IV Cosmic Shear with Modified Gravity and Model-independent Screening

Authors: Maria Tsedrik, Benjamin Bose, Pedro Carrilho, Alkistis Pourtsidou, Sefa Pamuk, Santiago Casas, Julien Lesgourgues

Abstract: We forecast constraints on minimal model-independent parametrisations of several Modified Gravity theories using mock Stage-IV cosmic shear data. We include nonlinear effects and screening, which ensures recovery of General Relativity on small scales. We introduce a power spectrum emulator to accelerate our analysis and evaluate the robustness of the growth index parametrisation with respect to tw… ▽ More We forecast constraints on minimal model-independent parametrisations of several Modified Gravity theories using mock Stage-IV cosmic shear data. We include nonlinear effects and screening, which ensures recovery of General Relativity on small scales. We introduce a power spectrum emulator to accelerate our analysis and evaluate the robustness of the growth index parametrisation with respect to two cosmologies: $Λ$CDM and the normal branch of the DGP model. We forecast the uncertainties on the growth index $γ$ to be of the order $\sim 10\%$. We find that our halo-model based screening approach demonstrates excellent performance, meeting the precision requirements of Stage-IV surveys. However, neglecting the screening transition results in biased predictions for cosmological parameters. We find that the screening transition shows significant degeneracy with baryonic feedback, requiring a much better understanding of baryonic physics for its detection. Massive neutrinos effects are less prominent and challenging to detect solely with cosmic shear data. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 22 pages, 16 figures, 2 tables plus Appendices. Our power spectrum boost emulators for extended cosmologies are publicly available at https://github.com/nebblu/ReACT-emus/tree/main/emulators

Beyond-$Λ$CDM Cosmologies with Bayesian Neural Networks II: Massive Neutrinos, Baryonic Feedback, and the Theoretical Error

Authors: L. Thummel, B. Bose, A. Pourtsidou, L. Lombriser

Abstract: We study the capacity of Bayesian Neural Networks (BNNs) to detect new physics in the dark matter power spectrum. As in previous studies, the Bayesian Cosmological Network (BaCoN) classifies spectra into one of 5 classes: $Λ$CDM, $f(R)$, $w$CDM, Dvali-Gabadaze-Porrati (DGP) gravity and a 'random' class, with this work extending it to include the effects of massive neutrinos and baryonic feedback.… ▽ More We study the capacity of Bayesian Neural Networks (BNNs) to detect new physics in the dark matter power spectrum. As in previous studies, the Bayesian Cosmological Network (BaCoN) classifies spectra into one of 5 classes: $Λ$CDM, $f(R)$, $w$CDM, Dvali-Gabadaze-Porrati (DGP) gravity and a 'random' class, with this work extending it to include the effects of massive neutrinos and baryonic feedback. We further develop the treatment of theoretical errors in BaCoN-II, investigating several approaches and identifying the one that best allows the trained network to generalise to other power spectrum modelling prescriptions. In particular, we compare power spectra data produced by EuclidEmulator2, HMcode and halofit, all supplemented with the halo model reaction to model beyond-$Λ$CDM physics. We investigate BNN-classifiers trained on these sets of spectra, adding in Stage-IV survey noise and various theoretical error models. Using our optimal theoretical error model, our fiducial classifier achieves a total classification accuracy of $\sim$ 95% when it is trained on EuclidEmulator2-based spectra with modification parameters drawn from a Gaussian distribution centred around $Λ$CDM ($f(R)$: $σ_{fR0} = 10^{-5.5}$, DGP: $σ_{r\mathrm{c}} = 0.173$, $w$CDM: $σ_{w0} = 0.097$, $σ_{wa}=0.32$). This strengthens the promise of this method to glean the maximal amount of unbiased gravitational and cosmological information from forthcoming Stage-IV galaxy surveys. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Dark Scattering: accelerated constraints from KiDS-1000 with $\tt{ReACT}$ and $\tt{CosmoPower}$

Authors: Karim Carrion, Pedro Carrilho, Alessio Spurio Mancini, Alkistis Pourtsidou, Juan Carlos Hidalgo

Abstract: We present constraints on the Dark Scattering model through cosmic shear measurements from the Kilo Degree Survey (KiDS-1000), using an accelerated pipeline with novel emulators produced with $\tt{CosmoPower}$. Our main emulator, for the Dark Scattering non-linear matter power spectrum, is trained on predictions from the halo model reaction framework, previously validated against simulations. Addi… ▽ More We present constraints on the Dark Scattering model through cosmic shear measurements from the Kilo Degree Survey (KiDS-1000), using an accelerated pipeline with novel emulators produced with $\tt{CosmoPower}$. Our main emulator, for the Dark Scattering non-linear matter power spectrum, is trained on predictions from the halo model reaction framework, previously validated against simulations. Additionally, we include the effects of baryonic feedback from $\tt{HMcode2016}$, whose contribution is also emulated. We analyse the complete set of statistics of KiDS-1000, namely Band Powers, COSEBIs and Correlation Functions, for Dark Scattering in two distinct cases. In the first case, taking into account only KiDS cosmic shear data, we constrain the amplitude of the dark energy - dark matter interaction to be $\vert A_{\rm ds} \vert \lesssim 20$ $\rm b/GeV$ at 68% C.L. Furthermore, we add information from the cosmic microwave background (CMB) from Planck, along with baryon acoustic oscillations (BAO) from 6dFGS, SDSS and BOSS, approximating a combined weak lensing + CMB + BAO analysis. From this combination, we constrain $A_{\rm ds} = 10.6^{+4.5}_{-7.3}$ $\rm b/GeV$ at 68% C.L. We confirm that with this estimated value of $A_{\rm ds}$ the interacting model considered in this work offers a promising alternative to solve the $S_8$ tension. △ Less

Submitted 11 July, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 8 pages, 8 figures, and 1 appendix. Version 2: Including additional explanations and minor corrections. Accepted for publication by MNRAS. Emulators available in GitHub repository https://github.com/karimpsi22/DS-emulators

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

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

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

Submitted 1 December, 2023; originally announced December 2023.

Comments: 38 pages, 19 figures

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

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

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

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

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

Euclid preparation. TBD. Forecast impact of super-sample covariance on 3x2pt analysis with Euclid

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

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

Submitted 24 October, 2023; originally announced October 2023.

Comments: 22 pages, 13 figures

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

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

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

Submitted 23 October, 2023; originally announced October 2023.

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

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

Euclid: Constraints on f(R) cosmologies from the spectroscopic and photometric primary probes

Authors: S. Casas, V. F. Cardone, D. Sapone, N. Frusciante, F. Pace, G. Parimbelli, M. Archidiacono, K. Koyama, I. Tutusaus, S. Camera, M. Martinelli, V. Pettorino, Z. Sakr, L. Lombriser, A. Silvestri, M. Pietroni, F. Vernizzi, M. Kunz, T. Kitching, A. Pourtsidou, F. Lacasa, C. Carbone, J. Garcia-Bellido, N. Aghanim, B. Altieri , et al. (104 additional authors not shown)

Abstract: $\textit{Euclid}$ will provide a powerful compilation of data including spectroscopic redshifts, the angular clustering of galaxies, weak lensing cosmic shear, and the cross-correlation of these last two photometric observables. In this study we extend recently presented $\textit{Euclid}$ forecasts into the Hu-Sawicki $f(R)… ▽ More $\textit{Euclid}$ will provide a powerful compilation of data including spectroscopic redshifts, the angular clustering of galaxies, weak lensing cosmic shear, and the cross-correlation of these last two photometric observables. In this study we extend recently presented $\textit{Euclid}$ forecasts into the Hu-Sawicki $f(R)$ cosmological model, a popular extension of the Hilbert-Einstein action that introduces an universal modified gravity force in a scale-dependent way. Our aim is to estimate how well future $\textit{Euclid}$ data will be able to constrain the extra parameter of the theory, $f_{R0}$, for the range in which this parameter is still allowed by current observations. For the spectroscopic probe, we use a phenomenological approach for the scale dependence of the growth of perturbations in the terms related to baryon acoustic oscillations and redshift-space distortions. For the photometric observables, we use a fitting formula that captures the modifications in the non-linear matter power spectrum caused by the $f(R)$ model. We show that, in an optimistic setting, and for a fiducial value of $f_{R0} = 5 \times 10^{-6}$, $\textit{Euclid}$ alone will be able to constrain the additional parameter $\log f_{R0}$ at the $3\%$ level, using spectroscopic galaxy clustering alone; at the $1.4\%$ level, using the combination of photometric probes on their own; and at the $1\%$ level, using the combination of spectroscopic and photometric observations. This last constraint corresponds to an error of the order of $6 \times 10^{-7}$ at the $1σ$ level on the model parameter $f_{R0} = 5 \times 10^{-6}$. We report also forecasted constraints for $f_{R0} = 5 \times 10^{-5}$ and $f_{R0} = 5 \times 10^{-7}$ and show that in the optimistic scenario, $\textit{Euclid}$ will be able to distinguish these models from $Λ\mathrm{CDM}$ at more than 3$σ$. (abridged) △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 10 Figures, 3 Tables, 1 Appendix

Modified gravity and massive neutrinos: constraints from the full shape analysis of BOSS galaxies and forecasts for Stage IV surveys

Authors: Chiara Moretti, Maria Tsedrik, Pedro Carrilho, Alkistis Pourtsidou

Abstract: We constrain the growth index $γ$ by performing a full-shape analysis of the power spectrum multipoles measured from the BOSS DR12 data. We adopt a theoretical model based on the Effective Field theory of the Large Scale Structure (EFTofLSS) and focus on two different cosmologies: $γ$CDM and $γν$CDM, where we also vary the total neutrino mass. We explore different choices for the priors on the pri… ▽ More We constrain the growth index $γ$ by performing a full-shape analysis of the power spectrum multipoles measured from the BOSS DR12 data. We adopt a theoretical model based on the Effective Field theory of the Large Scale Structure (EFTofLSS) and focus on two different cosmologies: $γ$CDM and $γν$CDM, where we also vary the total neutrino mass. We explore different choices for the priors on the primordial amplitude $A_s$ and spectral index $n_s$, finding that informative priors are necessary to alleviate degeneracies between the parameters and avoid strong projection effects in the posterior distributions. Our tightest constraints are obtained with 3$σ$ Planck priors on $A_s$ and $n_s$: we obtain $γ= 0.647 \pm 0.085$ for $γ$CDM and $γ= 0.612^{+0.075}_{-0.090}$, $M_ν< 0.30$ for $γν$CDM at 68\% c.l., in both cases $\sim 1σ$ consistent with the $Λ$CDM prediction $γ\simeq 0.55$. Additionally, we produce forecasts for a Stage-IV spectroscopic galaxy survey, focusing on a DESI-like sample. We fit synthetic data-vectors for three different galaxy samples generated at three different redshift bins, both individually and jointly. Focusing on the constraining power of the Large Scale Structure alone, we find that forthcoming data can give an improvement of up to $\sim 85\%$ in the measurement of $γ$ with respect to the BOSS dataset when no CMB priors are imposed. On the other hand, we find the neutrino mass constraints to be only marginally better than the current ones, with future data able to put an upper limit of $M_ν< 0.27~{\rm eV}$. This result can be improved with the inclusion of Planck priors on the primordial parameters, which yield $M_ν< 0.18~{\rm eV}$. △ Less

Submitted 12 December, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: updated with accepted version

Journal ref: JCAP12(2023)025

The foreground transfer function for HI intensity mapping signal reconstruction: MeerKLASS and precision cosmology applications

Authors: Steven Cunnington, Laura Wolz, Philip Bull, Isabella P. Carucci, Keith Grainge, Melis O. Irfan, Yichao Li, Alkistis Pourtsidou, Mario G. Santos, Marta Spinelli, Jingying Wang

Abstract: Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensi… ▽ More Blind cleaning methods are currently the preferred strategy for handling foreground contamination in single-dish HI intensity mapping surveys. Despite the increasing sophistication of blind techniques, some signal loss will be inevitable across all scales. Constructing a corrective transfer function using mock signal injection into the contaminated data has been a practice relied on for HI intensity mapping experiments. However, assessing whether this approach is viable for future intensity mapping surveys where precision cosmology is the aim, remains unexplored. In this work, using simulations, we validate for the first time the use of a foreground transfer function to reconstruct power spectra of foreground-cleaned low-redshift intensity maps and look to expose any limitations. We reveal that even when aggressive foreground cleaning is required, which causes ${>}\,50\%$ negative bias on the largest scales, the power spectrum can be reconstructed using a transfer function to within sub-percent accuracy. We specifically outline the recipe for constructing an unbiased transfer function, highlighting the pitfalls if one deviates from this recipe, and also correctly identify how a transfer function should be applied in an auto-correlation power spectrum. We validate a method that utilises the transfer function variance for error estimation in foreground-cleaned power spectra. Finally, we demonstrate how incorrect fiducial parameter assumptions (up to ${\pm}100\%$ bias) in the generation of mocks, used in the construction of the transfer function, do not significantly bias signal reconstruction or parameter inference (inducing ${<}\,5\%$ bias in recovered values). △ Less

Submitted 23 May, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 25 pages, 20 figures. See Figure 4 for the main demonstration of the transfer function's performance for reconstructing signal loss from foreground cleaning. Accepted for publication in MNRAS

Euclid preparation: XXVIII. Modelling of the weak lensing angular power spectrum

Authors: Euclid Collaboration, A. C. Deshpande, T. Kitching, A. Hall, M. L. Brown, N. Aghanim, L. Amendola, N. Auricchio, M. Baldi, R. Bender, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, G. P. Candini, V. Capobianco, C. Carbone, V. F. Cardone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou , et al. (178 additional authors not shown)

Abstract: This work considers which higher-order effects in modelling the cosmic shear angular power spectra must be taken into account for Euclid. We identify which terms are of concern, and quantify their individual and cumulative impact on cosmological parameter inference from Euclid. We compute the values of these higher-order effects using analytic expressions, and calculate the impact on cosmological… ▽ More This work considers which higher-order effects in modelling the cosmic shear angular power spectra must be taken into account for Euclid. We identify which terms are of concern, and quantify their individual and cumulative impact on cosmological parameter inference from Euclid. We compute the values of these higher-order effects using analytic expressions, and calculate the impact on cosmological parameter estimation using the Fisher matrix formalism. We review 24 effects and find the following potentially need to be accounted for: the reduced shear approximation, magnification bias, source-lens clustering, source obscuration, local Universe effects, and the flat Universe assumption. Upon computing these explicitly, and calculating their cosmological parameter biases, using a maximum multipole of $\ell=5000$, we find that the magnification bias, source-lens clustering, source obscuration, and local Universe terms individually produce significant ($\,>0.25σ$) cosmological biases in one or more parameters, and accordingly must be accounted for. In total, over all effects, we find biases in $Ω_{\rm m}$, $Ω_{\rm b}$, $h$, and $σ_{8}$ of $0.73σ$, $0.28σ$, $0.25σ$, and $-0.79σ$, respectively, for flat $Λ$CDM. For the $w_0w_a$CDM case, we find biases in $Ω_{\rm m}$, $Ω_{\rm b}$, $h$, $n_{\rm s}$, $σ_{8}$, and $w_a$ of $1.49σ$, $0.35σ$, $-1.36σ$, $1.31σ$, $-0.84σ$, and $-0.35σ$, respectively; which are increased relative to the $Λ$CDM due to additional degeneracies as a function of redshift and scale. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 20 pages, submitted to A&A

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

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

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

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

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

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

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

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

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

Submitted 23 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures

MSC Class: 85A40

Fast and accurate predictions of the nonlinear matter power spectrum for general models of Dark Energy and Modified Gravity

Authors: B. Bose, M. Tsedrik, J. Kennedy, L. Lombriser, A. Pourtsidou, A. Taylor

Abstract: We embed linear and nonlinear parametrisations of beyond standard cosmological physics in the halo model reaction framework, providing a model-independent prescription for the nonlinear matter power spectrum. As an application, we focus on Horndeski theories, using the Effective Field Theory of Dark Energy (EFTofDE) to parameterise linear and quasi-nonlinear perturbations. In the nonlinear regime… ▽ More We embed linear and nonlinear parametrisations of beyond standard cosmological physics in the halo model reaction framework, providing a model-independent prescription for the nonlinear matter power spectrum. As an application, we focus on Horndeski theories, using the Effective Field Theory of Dark Energy (EFTofDE) to parameterise linear and quasi-nonlinear perturbations. In the nonlinear regime we investigate both a nonlinear parameterised-post Friedmannian (nPPF) approach as well as a physically motivated and approximate phenomenological model based on the error function (Erf). We compare the parameterised approaches' predictions of the nonlinear matter power spectrum to the exact solutions, as well as state-of-the-art emulators, in an evolving dark energy scenario and two well studied modified gravity models, finding sub-percent agreement in the reaction using the Erf model at $z\leq1$ and $k\leq 5~h/{\rm Mpc}$. This suggests only an additional 3 free constants, above the background and linear theory parameters, are sufficient to model nonlinear, non-standard cosmology in the matter power spectrum at scales down to $k \leq 3h~/{\rm Mpc}$ within $2\%$ accuracy. We implement the parametrisations into ver.2.0 of the ReACT code. △ Less

Submitted 6 January, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: Journal accepted version. 21 pages, 10 figures. Download ReACT v2: https://github.com/nebblu/ACTio-ReACTio

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

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

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

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

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

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

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

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

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

Submitted 26 September, 2022; originally announced September 2022.

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

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

Euclid preparation. XXIV. Calibration of the halo mass function in $Λ(ν)$CDM cosmologies

Authors: Euclid Collaboration, T. Castro, A. Fumagalli, R. E. Angulo, S. Bocquet, S. Borgani, C. Carbone, J. Dakin, K. Dolag, C. Giocoli, P. Monaco, A. Ragagnin, A. Saro, E. Sefusatti, M. Costanzi, A. M. C. Le Brun, P. -S. Corasaniti, A. Amara, L. Amendola, M. Baldi, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, S. Camera , et al. (157 additional authors not shown)

Abstract: Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be… ▽ More Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be subdominant with respect to other sources of uncertainty in recovering cosmological parameters from Euclid cluster counts. Our model is calibrated against a suite of N-body simulations using a Bayesian approach taking into account systematic errors arising from numerical effects in the simulation. First, we test the convergence of HMF predictions from different N-body codes, by using initial conditions generated with different orders of Lagrangian Perturbation theory, and adopting different simulation box sizes and mass resolution. Then, we quantify the effect of using different halo-finder algorithms, and how the resulting differences propagate to the cosmological constraints. In order to trace the violation of universality in the HMF, we also analyse simulations based on initial conditions characterised by scale-free power spectra with different spectral indexes, assuming both Einstein--de Sitter and standard $Λ$CDM expansion histories. Based on these results, we construct a fitting function for the HMF that we demonstrate to be sub-percent accurate in reproducing results from 9 different variants of the $Λ$CDM model including massive neutrinos cosmologies. The calibration systematic uncertainty is largely sub-dominant with respect to the expected precision of future mass-observation relations; with the only notable exception of the effect due to the halo finder, that could lead to biased cosmological inference. △ Less

Submitted 16 March, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 25 pages, 21 figures, 5 tables, 3 appendixes; v2 matches published version

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

Cosmology with the EFTofLSS and BOSS: dark energy constraints and a note on priors

Authors: Pedro Carrilho, Chiara Moretti, Alkistis Pourtsidou

Abstract: We analyse the BOSS DR12 galaxy power spectrum data jointly with BAO data for three models of dark energy. We use recent measurements using a windowless estimator, and an independent and fast pipeline based on EFTofLSS implemented via the FAST-PT algorithm to compute the redshift-space loop corrections. We accelerate our analysis by using the BACCO linear emulator instead of a Boltzmann solver. We… ▽ More We analyse the BOSS DR12 galaxy power spectrum data jointly with BAO data for three models of dark energy. We use recent measurements using a windowless estimator, and an independent and fast pipeline based on EFTofLSS implemented via the FAST-PT algorithm to compute the redshift-space loop corrections. We accelerate our analysis by using the BACCO linear emulator instead of a Boltzmann solver. We perform two sets of analyses: one with $3σ$ Planck priors on $A_s$ and $n_s$, and another that is CMB-free, without such priors. Firstly, we study $Λ$CDM, reproducing previous results obtained with the same estimator. We find a low value of $A_s$ in the CMB-free case, in agreement with many previous full-shape analyses of the BOSS data. We then study $w$CDM, finding a lower value of the amplitude in the CMB-free run, coupled with a preference for phantom dark energy with $w=-1.17^{+0.12}_{-0.11}$, again in broad agreement with previous results. Finally, we investigate the dark scattering model, which we label $wA$CDM. In the CMB-free analysis, we find a large degeneracy between the interaction strength $A$ and the amplitude $A_s$, hampering measurements of those parameters. On the contrary, in our run with a CMB prior, we are able to constrain the dark energy parameters to be $w=-0.972^{+0.036}_{-0.029}$ and $A = 3.9^{+3.2}_{-3.7}$, which show a 1$σ$ hint of interacting dark energy. This is the first measurement of this parameter and demonstrates the ability of this model to alleviate the $σ_8$ tension. Our analysis can be used as a guide for any model with scale-independent growth. Finally, we study the dependence of the results on the priors of the nuisance parameters and find these priors to be informative, with their broadening creating shifts in the contours. We argue for an in depth study of this issue, which can affect current and forthcoming analyses of LSS data. △ Less

Submitted 16 January, 2023; v1 submitted 29 July, 2022; originally announced July 2022.

Comments: 34+17 pages, 18 figures. Baseline analyses and results in Section 3, discussion on priors in appendix A. V2: More detailed comparisons with previous work and additional references. Matches version accepted for publication in JCAP

Interacting dark energy from the joint analysis of the power spectrum and bispectrum multipoles with the EFTofLSS

Authors: Maria Tsedrik, Chiara Moretti, Pedro Carrilho, Federico Rizzo, Alkistis Pourtsidou

Abstract: Interacting dark energy models have been suggested as alternatives to the standard cosmological model, $Λ$CDM. We focus on a phenomenologically interesting class of dark scattering models that is characterised by pure momentum exchange between dark energy and dark matter. This model extends the parameter space with respect to $Λ$CDM by two parameters, $w$ and $A$, which define the dark energy equa… ▽ More Interacting dark energy models have been suggested as alternatives to the standard cosmological model, $Λ$CDM. We focus on a phenomenologically interesting class of dark scattering models that is characterised by pure momentum exchange between dark energy and dark matter. This model extends the parameter space with respect to $Λ$CDM by two parameters, $w$ and $A$, which define the dark energy equation of state and the strength of the coupling between dark energy and dark matter, respectively. In order to test non-standard cosmologies with Stage-IV galaxy clustering surveys, it is crucial to model mildly nonlinear scales and perform precision vs accuracy tests. We use the Effective Field Theory of Large-Scale Structure, and we perform validation tests by means of an MCMC analysis using a large set of N-body simulations. We find that adding the bispectrum monopole to the power spectrum multipoles improves the constraints on the dark energy parameters by $\sim 30 \%$ for $k_{\mathrm{max}, B}^{l=0} = 0.11$ $h$ Mpc$^{-1}$ without introducing biases in the parameter estimation. We also find that the same improvement can be achieved with more moderate scale cuts and the use of bias relations, or with the addition of the bispectrum quadrupole. Finally, we study degeneracies between the dark energy parameters and the scalar amplitude $A_\mathrm{s}$ and discuss the corresponding projection effects, as well as degeneracies with other cosmological parameters. △ Less

Submitted 23 January, 2023; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 22 pages, 13 figures; V2: Matches version accepted for publication in MNRAS, includes a new analysis of degeneracies with cosmological parameters

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Authors: Euclid Collaboration, L. Bisigello, C. J. Conselice, M. Baes, M. Bolzonella, M. Brescia, S. Cavuoti, O. Cucciati, A. Humphrey, L. K. Hunt, C. Maraston, L. Pozzetti, C. Tortora, S. E. van Mierlo, N. Aghanim, N. Auricchio, M. Baldi, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone , et al. (174 additional authors not shown)

Abstract: Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshift… ▽ More Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning. △ Less

Submitted 4 January, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: accepted for publication in MNRAS, 21 pages, 22 figures, 6 tables

Interferometric HI intensity mapping: perturbation theory predictions and foreground removal effects

Authors: Alkistis Pourtsidou

Abstract: We provide perturbation theory predictions for the HI intensity mapping power spectrum multipoles using the Effective Field Theory of Large Scale Structure (EFTofLSS), which should allow us to constrain cosmological parameters exploiting mildly nonlinear scales. Assuming survey specifications typical of proposed interferometric HI intensity mapping experiments like CHORD and PUMA, and realistic ra… ▽ More We provide perturbation theory predictions for the HI intensity mapping power spectrum multipoles using the Effective Field Theory of Large Scale Structure (EFTofLSS), which should allow us to constrain cosmological parameters exploiting mildly nonlinear scales. Assuming survey specifications typical of proposed interferometric HI intensity mapping experiments like CHORD and PUMA, and realistic ranges of validity for the perturbation theory modelling, we run mock full shape MCMC analyses at a redshift bin centred at $z=0.5$, and compare with Stage-IV optical galaxy surveys. We include the impact of 21cm foreground removal using simulations-based prescriptions, and quantify the effects on the precision and accuracy of the parameter estimation. We vary 11 parameters in total: 3 cosmological parameters, 7 bias and counterterms parameters, and the HI brightness temperature. Amongst them, the 4 parameters of interest are: the cold dark matter density, $ω_{\rm c}$, the Hubble parameter, $h$, the primordial amplitude of the power spectrum, $A_{\rm s}$, and the linear HI bias, $b_1$. For the best case scenario, we obtain unbiased constraints on all parameters with $<3\%$ errors at $68\%$ confidence level. When we include the foreground removal effects, the parameter estimation becomes strongly biased for $ω_{\rm c}, h$, and $b_1$, while $A_{\rm s}$ is less biased ($< 2σ$). We find that scale cuts $k_{\rm min} \geq 0.03 \, h/{\mathrm{Mpc}}$ are required to return accurate estimates for $ω_{\rm c}$ and $h$, at the price of a decrease in the precision, while $b_1$ remains strongly biased. We comment on the implications of these results for real data analyses. △ Less

Submitted 9 January, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 11 pages, 1 table, 10 figures; matches version accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 519 (2023) 4, 6246-6256

Euclid: Forecasts from the void-lensing cross-correlation

Authors: M. Bonici, C. Carbone, S. Davini, P. Vielzeuf, L. Paganin, V. Cardone, N. Hamaus, A. Pisani, A. J. Hawken, A. Kovacs, S. Nadathur, S. Contarini, G. Verza, I. Tutusaus, F. Marulli, L. Moscardini, M. Aubert, C. Giocoli, A. Pourtsidou, S. Camera, S. Escoffier, A. Caminata, M. Martinelli, M. Pallavicini, V. Pettorino , et al. (107 additional authors not shown)

Abstract: The Euclid space telescope will survey a large dataset of cosmic voids traced by dense samples of galaxies. In this work we estimate its expected performance when exploiting angular photometric void clustering, galaxy weak lensing and their cross-correlation. To this aim, we implement a Fisher matrix approach tailored for voids from the Euclid photometric dataset and present the first forecasts on… ▽ More The Euclid space telescope will survey a large dataset of cosmic voids traced by dense samples of galaxies. In this work we estimate its expected performance when exploiting angular photometric void clustering, galaxy weak lensing and their cross-correlation. To this aim, we implement a Fisher matrix approach tailored for voids from the Euclid photometric dataset and present the first forecasts on cosmological parameters that include the void-lensing correlation. We examine two different probe settings, pessimistic and optimistic, both for void clustering and galaxy lensing. We carry out forecast analyses in four model cosmologies, accounting for a varying total neutrino mass, $M_ν$, and a dynamical dark energy (DE) equation of state, $w(z)$, described by the CPL parametrisation. We find that void clustering constraints on $h$ and $Ω_b$ are competitive with galaxy lensing alone, while errors on $n_s$ decrease thanks to the orthogonality of the two probes in the 2D-projected parameter space. We also note that, as a whole, the inclusion of the void-lensing cross-correlation signal improves parameter constraints by $10-15\%$, and enhances the joint void clustering and galaxy lensing Figure of Merit (FoM) by $10\%$ and $25\%$, in the pessimistic and optimistic scenarios, respectively. Finally, when further combining with the spectroscopic galaxy clustering, assumed as an independent probe, we find that, in the most competitive case, the FoM increases by a factor of 4 with respect to the combination of weak lensing and spectroscopic galaxy clustering taken as independent probes. The forecasts presented in this work show that photometric void-clustering and its cross-correlation with galaxy lensing deserve to be exploited in the data analysis of the Euclid galaxy survey and promise to improve its constraining power, especially on $h$, $Ω_b$, the neutrino mass, and the DE evolution. △ Less

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

Comments: 22 pages, 9 figures - published in A&A

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

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

Authors: Euclid Collaboration, R. Saglia, S. De Nicola, M. Fabricius, V. Guglielmo, J. Snigula, R. Zöller, R. Bender, J. Heidt, D. Masters, D. Stern, S. Paltani, A. Amara, N. Auricchio, M. Baldi, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, M. Castellano , et al. (161 additional authors not shown)

Abstract: The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak l… ▽ More The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration. △ Less

Submitted 7 September, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: A&A language edited version

Journal ref: A&A 664, A196 (2022)

HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies

Authors: Steven Cunnington, Yichao Li, Mario G. Santos, Jingying Wang, Isabella P. Carucci, Melis O. Irfan, Alkistis Pourtsidou, Marta Spinelli, Laura Wolz, Paula S. Soares, Chris Blake, Philip Bull, Brandon Engelbrecht, José Fonseca, Keith Grainge, Yin-Zhe Ma

Abstract: We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7σ$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($Ω_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constrai… ▽ More We present a detection of correlated clustering between MeerKAT radio intensity maps and galaxies from the WiggleZ Dark Energy Survey. We find a $7.7σ$ detection of the cross-correlation power spectrum, the amplitude of which is proportional to the product of the HI density fraction ($Ω_{\rm HI}$), HI bias ($b_{\rm HI}$) and the cross-correlation coefficient ($r$). We therefore obtain the constraint $Ω_{\rm HI} b_{\rm HI} r\,{=}\,[0.86\,{\pm}\,0.10\,({\rm stat})\,{\pm}\,0.12\,({\rm sys})]\,{\times}\,10^{-3}$, at an effective scale of $k_{\rm eff}\,{\sim}\,0.13\,h\,{\rm Mpc}^{-1}$. The intensity maps were obtained from a pilot survey with the MeerKAT telescope, a 64-dish pathfinder array to the SKA Observatory (SKAO). The data were collected from 10.5 hours of observations using MeerKAT's L-band receivers over six nights covering the 11hr field of WiggleZ, in the frequency range $1015-973\,{\rm MHz}$ (0.400$\,{<}\,z\,{<}\,$0.459 in redshift). This detection is the first practical demonstration of the multi-dish auto-correlation intensity mapping technique for cosmology. This marks an important milestone in the roadmap for the cosmology science case with the full SKAO. △ Less

Submitted 25 October, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 11 pages, 7 figures. See Fig.5 for the main result, showing the cross-power spectrum between MeerKAT and WiggleZ, revealing the 7.7σdetection. v2 is the version accepted for publication by MNRAS. The main changes from v1 were in introducing the radio intensity mapping data

Unveiling the Universe with Emerging Cosmological Probes

Authors: Michele Moresco, Lorenzo Amati, Luca Amendola, Simon Birrer, John P. Blakeslee, Michele Cantiello, Andrea Cimatti, Jeremy Darling, Massimo Della Valle, Maya Fishbach, Claudio Grillo, Nico Hamaus, Daniel Holz, Luca Izzo, Raul Jimenez, Elisabeta Lusso, Massimo Meneghetti, Ester Piedipalumbo, Alice Pisani, Alkistis Pourtsidou, Lucia Pozzetti, Miguel Quartin, Guido Risaliti, Piero Rosati, Licia Verde

Abstract: The detection of the accelerated expansion of the Universe has been one of the major breakthroughs in modern cosmology. Several cosmological probes (CMB, SNe Ia, BAO) have been studied in depth to better understand the nature of the mechanism driving this acceleration, and they are being currently pushed to their limits, obtaining remarkable constraints that allowed us to shape the standard cosmol… ▽ More The detection of the accelerated expansion of the Universe has been one of the major breakthroughs in modern cosmology. Several cosmological probes (CMB, SNe Ia, BAO) have been studied in depth to better understand the nature of the mechanism driving this acceleration, and they are being currently pushed to their limits, obtaining remarkable constraints that allowed us to shape the standard cosmological model. In parallel to that, however, the percent precision achieved has recently revealed apparent tensions between measurements obtained from different methods. These are either indicating some unaccounted systematic effects, or are pointing toward new physics. Following the development of CMB, SNe, and BAO cosmology, it is critical to extend our selection of cosmological probes. Novel probes can be exploited to validate results, control or mitigate systematic effects, and, most importantly, to increase the accuracy and robustness of our results. This review is meant to provide a state-of-art benchmark of the latest advances in emerging beyond-standard cosmological probes. We present how several different methods can become a key resource for observational cosmology. In particular, we review cosmic chronometers, quasars, gamma-ray bursts, standard sirens, lensing time-delay with galaxies and clusters, cosmic voids, neutral hydrogen intensity mapping, surface brightness fluctuations, stellar ages of the oldest objects, secular redshift drift, and clustering of standard candles. The review describes the method, systematics, and results of each probe in a homogeneous way, giving the reader a clear picture of the available innovative methods that have been introduced in recent years and how to apply them. The review also discusses the potential synergies and complementarities between the various probes, exploring how they will contribute to the future of modern cosmology. △ Less

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

Comments: Invited review article for Living Reviews in Relativity. Accepted version. 191 pages, 54 figures, 11 tables

On the road to percent accuracy VI: the nonlinear power spectrum for interacting dark energy with baryonic feedback and massive neutrinos

Authors: Pedro Carrilho, Karim Carrion, Benjamin Bose, Alkistis Pourtsidou, Juan Carlos Hidalgo, Lucas Lombriser, Marco Baldi

Abstract: Understanding nonlinear structure formation is crucial for fully exploring the data generated by stage IV surveys, requiring accurate modelling of the power spectrum. This is challenging for deviations from $Λ$CDM, but we must ensure that alternatives are well tested, to avoid false detections. We present an extension of the halo model reaction framework for interacting dark energy. We modify the… ▽ More Understanding nonlinear structure formation is crucial for fully exploring the data generated by stage IV surveys, requiring accurate modelling of the power spectrum. This is challenging for deviations from $Λ$CDM, but we must ensure that alternatives are well tested, to avoid false detections. We present an extension of the halo model reaction framework for interacting dark energy. We modify the halo model including the additional force present in the Dark Scattering model and implement it into ReACT. The reaction is combined with a pseudo spectrum from EuclidEmulator2 and compared to N-body simulations. Using standard mass function and concentration-mass relation, we find predictions to be 1 % accurate at $z=0$ up to $k=0.8~h/{\rm Mpc}$ for the largest interaction strength tested ($ξ=50$ b/GeV), improving to $2~h/{\rm Mpc}$ at $z=1$. For smaller interaction strength ($10$ b/GeV), we find 1 % agreement at $z=1$ up to scales above $3.5~h/{\rm Mpc}$, being close to $1~h/{\rm Mpc}$ at $z=0$. Finally, we improve our predictions with the inclusion of baryonic feedback and massive neutrinos and study degeneracies between the effects of these contributions and those of the interaction. Limiting the scales to where our modelling is 1 % accurate, we find a degeneracy between the interaction and feedback, but not with massive neutrinos. We expect the degeneracy with feedback to be resolvable by including smaller scales. This work represents the first analytical tool for calculating the nonlinear spectrum for interacting dark energy models. △ Less

Submitted 7 March, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 11 pages, 4 figures. Version 2: version accepted for publication in MNRAS, including additional explanations and minor corrections. The version of ReACT used in this work is available at https://github.com/PedroCarrilho/ReACT/tree/react_with_interact_baryons

KiDS-1000 Cosmology: machine learning -- accelerated constraints on Interacting Dark Energy with COSMOPOWER

Authors: A. Spurio Mancini, A. Pourtsidou

Abstract: We derive constraints on a coupled quintessence model with pure momentum exchange from the public $\sim$1000 deg$^2$ cosmic shear measurements from the Kilo-Degree Survey and the $\it{Planck}$ 2018 Cosmic Microwave Background data. We compare this model with $Λ$CDM and find similar $χ^2$ and log-evidence values. We accelerate parameter estimation by sourcing cosmological power spectra from the neu… ▽ More We derive constraints on a coupled quintessence model with pure momentum exchange from the public $\sim$1000 deg$^2$ cosmic shear measurements from the Kilo-Degree Survey and the $\it{Planck}$ 2018 Cosmic Microwave Background data. We compare this model with $Λ$CDM and find similar $χ^2$ and log-evidence values. We accelerate parameter estimation by sourcing cosmological power spectra from the neural network emulator COSMOPOWER. We highlight the necessity of such emulator-based approaches to reduce the computational runtime of future similar analyses, particularly from Stage IV surveys. As an example, we present MCMC forecasts on the same coupled quintessence model for a $\it{Euclid}$-like survey, revealing degeneracies between the coupled quintessence parameters and the baryonic feedback and intrinsic alignment parameters, but also highlighting the large increase in constraining power Stage IV surveys will achieve. The contours are obtained in a few hours with COSMOPOWER, as opposed to the few months required with a Boltzmann code. △ Less

Submitted 24 February, 2022; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: 5 pages, 2 min. summary video available at https://youtu.be/c2x8hzApAgE. Emulators available in the COSMOPOWER GitHub repository, https://github.com/alessiospuriomancini/cosmopower. Matches version published in MNRAS Letters

Journal ref: Monthly Notices of the Royal Astronomical Society: Letters, Volume 512, Issue 1, May 2022, Pages L44-L48

Contributions from primordial non-Gaussianity and General Relativity to the galaxy power spectrum

Authors: Rebeca Martinez-Carrillo, Juan Carlos Hidalgo, Karim A. Malik, Alkistis Pourtsidou

Abstract: We compute the real space galaxy power spectrum, including the leading order effects of General Relativity and primordial non-Gaussianity from the $f_{\mathrm{NL}}$ and $g_{\mathrm{NL}}$ parameters. Such contributions come from the one-loop matter power spectrum terms dominant at large scales, and from the factors of the non-linear bias parameter $b_{\mathrm{NL}}$ (akin to the Newtonian $b_φ$). We… ▽ More We compute the real space galaxy power spectrum, including the leading order effects of General Relativity and primordial non-Gaussianity from the $f_{\mathrm{NL}}$ and $g_{\mathrm{NL}}$ parameters. Such contributions come from the one-loop matter power spectrum terms dominant at large scales, and from the factors of the non-linear bias parameter $b_{\mathrm{NL}}$ (akin to the Newtonian $b_φ$). We assess the detectability of these contributions in Stage-IV surveys. In particular, we note that specific values of the bias parameter may erase the primordial and relativistic contributions to the configuration space power spectrum. △ Less

Submitted 30 November, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: Version accepted for publication in JCAP

Journal ref: JCAP12(2021)025

SKAO HI Intensity Mapping: Blind Foreground Subtraction Challenge

Authors: Marta Spinelli, Isabella P. Carucci, Steven Cunnington, Stuart E. Harper, Melis O. Irfan, José Fonseca, Alkistis Pourtsidou, Laura Wolz

Abstract: Neutral Hydrogen Intensity Mapping (HI IM) surveys will be a powerful new probe of cosmology. However, strong astrophysical foregrounds contaminate the signal and their coupling with instrumental systematics further increases the data cleaning complexity. In this work, we simulate a realistic single-dish HI IM survey of a $5000$~deg$^2$ patch in the $950 - 1400$ MHz range, with both the MID telesc… ▽ More Neutral Hydrogen Intensity Mapping (HI IM) surveys will be a powerful new probe of cosmology. However, strong astrophysical foregrounds contaminate the signal and their coupling with instrumental systematics further increases the data cleaning complexity. In this work, we simulate a realistic single-dish HI IM survey of a $5000$~deg$^2$ patch in the $950 - 1400$ MHz range, with both the MID telescope of the SKA Observatory (SKAO) and MeerKAT, its precursor. We include a state-of-the-art HI simulations and explore different foreground models and instrumental effects such as non-homogeneous thermal noise and beam side-lobes. We perform the first Blind Foreground Subtraction Challenge for HI IM on these synthetic data-cubes, aiming to characterise the performance of available foreground cleaning methods with no prior knowledge of the sky components and noise level. Nine foreground cleaning pipelines joined the Challenge, based on statistical source separation algorithms, blind polynomial fitting, and an astrophysical-informed parametric fit to foregrounds. We devise metrics to compare the pipeline performances quantitatively. In general, they can recover the input maps' 2-point statistics within 20 per cent in the range of scales least affected by the telescope beam. However, spurious artefacts appear in the cleaned maps due to interactions between the foreground structure and the beam side-lobes. We conclude that it is fundamental to develop accurate beam deconvolution algorithms and test data post-processing steps carefully before cleaning. This study was performed as part of SKAO preparatory work by the HI IM Focus Group of the SKA Cosmology Science Working Group. △ Less

Submitted 20 October, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted version for publication in MNRAS

Momentum transfer models of interacting dark energy

Authors: Mark S. Linton, Robert Crittenden, Alkistis Pourtsidou

Abstract: We consider two models of interacting dark energy, both of which interact only through momentum exchange. One is a phenomenological one-parameter extension to $w$CDM, and the other is a coupled quintessence model described by a Lagrangian formalism. Using a variety of high and low redshift data sets, we perform a global fitting of cosmological parameters and compare to $Λ$CDM, uncoupled quintessen… ▽ More We consider two models of interacting dark energy, both of which interact only through momentum exchange. One is a phenomenological one-parameter extension to $w$CDM, and the other is a coupled quintessence model described by a Lagrangian formalism. Using a variety of high and low redshift data sets, we perform a global fitting of cosmological parameters and compare to $Λ$CDM, uncoupled quintessence, and $w$CDM. We find that the models are competitive with $Λ$CDM, even obtaining a better fit when certain data sets are included. △ Less

Submitted 3 August, 2022; v1 submitted 7 July, 2021; originally announced July 2021.

Comments: 22 pages, 10 figures, 3 tables; matches published version in JCAP

Interacting dark energy from redshift-space galaxy clustering

Authors: Pedro Carrilho, Chiara Moretti, Benjamin Bose, Katarina Markovič, Alkistis Pourtsidou

Abstract: Interacting dark energy models have been proposed as attractive alternatives to $Λ$CDM. Forthcoming Stage-IV galaxy clustering surveys will constrain these models, but they require accurate modelling of the galaxy power spectrum multipoles on mildly non-linear scales. In this work we consider a dark scattering model with a simple 1-parameter extension to $w$CDM - adding only $A$, which describes a… ▽ More Interacting dark energy models have been proposed as attractive alternatives to $Λ$CDM. Forthcoming Stage-IV galaxy clustering surveys will constrain these models, but they require accurate modelling of the galaxy power spectrum multipoles on mildly non-linear scales. In this work we consider a dark scattering model with a simple 1-parameter extension to $w$CDM - adding only $A$, which describes a pure momentum exchange between dark energy and dark matter. We then provide a comprehensive comparison of three approaches of modeling non-linearities, while including the effects of this dark sector coupling. We base our modeling of non-linearities on the two most popular perturbation theory approaches: TNS and EFTofLSS. To test the validity and precision of the modelling, we perform an MCMC analysis using simulated data corresponding to a $Λ$CDM fiducial cosmology and Stage-IV surveys specifications in two redshift bins, $z=0.5$ and $z=1$. We find the most complex EFTofLSS-based model studied to be better suited at both, describing the mock data up to smaller scales, and extracting the most information. Using this model, we forecast uncertainties on the dark energy equation of state, $w$, and on the interaction parameter, $A$, finding $σ_w=0.06$ and $σ_A=1.1$ b/GeV for the analysis at $z=0.5$ and $σ_w=0.06$ and $σ_A=2.0$ b/GeV for the analysis at $z=1$. In addition, we show that a false detection of exotic dark energy up to 3$σ$ would occur should the non-linear modelling be incorrect, demonstrating the importance of the validation stage for accurate interpretation of measurements. △ Less

Submitted 8 September, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 30 pages, 10 figures. Version accepted for publication in JCAP

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

Authors: Euclid Collaboration, S. Ilić, N. Aghanim, C. Baccigalupi, J. R. Bermejo-Climent, G. Fabbian, L. Legrand, D. Paoletti, M. Ballardini, M. Archidiacono, M. Douspis, F. Finelli, K. Ganga, C. Hernández-Monteagudo, M. Lattanzi, D. Marinucci, M. Migliaccio, C. Carbone, S. Casas, M. Martinelli, I. Tutusaus, P. Natoli, P. Ntelis, L. Pagano, L. Wenzl , et al. (185 additional authors not shown)

Abstract: The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological para… ▽ More The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data. △ Less

Submitted 10 September, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 20 pages, 8 figures, 5 tables, 1 appendix; updated to match version accepted by journal

Journal ref: A&A 657, A91 (2022)

The RSD Sorting Hat: Unmixing Radial Scales in Projection

Authors: Peter L. Taylor, Katarina Markovič, Alkistis Pourtsidou, Eric Huff

Abstract: Future data sets will enable cross-correlations between redshift space distortions (RSD) and weak lensing (WL). While photometric lensing and clustering cross-correlations have provided some of the tightest cosmological constraints to date, it is not well understood how to optimally perform similar RSD/WL joint analyses in a lossless way. RSD is typically measured in $3D$ redshift space, but WL is… ▽ More Future data sets will enable cross-correlations between redshift space distortions (RSD) and weak lensing (WL). While photometric lensing and clustering cross-correlations have provided some of the tightest cosmological constraints to date, it is not well understood how to optimally perform similar RSD/WL joint analyses in a lossless way. RSD is typically measured in $3D$ redshift space, but WL is inherently a projected signal, making angular statistics a natural choice for the combined analysis. Thus, we determine the amount of RSD information that can be extracted using projected statistics. Specifically we perform a Fisher analysis to forecast constraints and model bias comparing two different Fingers-of-God (FoG) models using both, the $3D$ power spectrum, $P(k, μ)$, and tomographic $C(\ell)$. We find that because naïve tomographic projection mixes large scales with poorly modelled nonlinear radial modes, it does not provide competitive constraints to the $3D$ RSD power spectrum without the model bias becoming unacceptably large. This is true even in the limit of narrow tomographic bins. In light of this we propose a new radial weighting scheme which unmixes radial RSD scales in projection yielding competitive constraints to the $3D$ RSD power spectrum, while keeping the model bias small. This work lays the groundwork for optimal joint analyses of RSD and cosmic shear. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 17 pages. 11 Figures. PRD submitted

Gaussian Process Regression for foreground removal in HI intensity mapping experiments

Authors: Paula S. Soares, Catherine A. Watkinson, Steven Cunnington, Alkistis Pourtsidou

Abstract: We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source Python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use G… ▽ More We apply for the first time Gaussian Process Regression (GPR) as a foreground removal technique in the context of single-dish, low redshift HI intensity mapping, and present an open-source Python toolkit for doing so. We use MeerKAT and SKA1-MID-like simulations of 21cm foregrounds (including polarisation leakage), HI cosmological signal and instrumental noise. We find that it is possible to use GPR as a foreground removal technique in this context, and that it is better suited in some cases to recover the HI power spectrum than Principal Component Analysis (PCA), especially on small scales. GPR is especially good at recovering the radial power spectrum, outperforming PCA when considering the full bandwidth of our data. Both methods are worse at recovering the transverse power spectrum, since they rely on frequency-only covariance information. When halving our data along frequency, we find that GPR performs better in the low frequency range, where foregrounds are brighter. It performs worse than PCA when frequency channels are missing, to emulate RFI flagging. We conclude that GPR is an excellent foreground removal option for the case of single-dish, low redshift HI intensity mapping in the absence of missing frequency channels. Our Python toolkit gpr4im and the data used in this analysis are publicly available on GitHub. △ Less

Submitted 4 October, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Comments: 19 pages, 1 table, 13 figures, main results in figures 6 and 9, code and data available at https://github.com/paulassoares/gpr4im. Version accepted by MNRAS

On the road to percent accuracy V: the non-linear power spectrum beyond $Λ$CDM with massive neutrinos and baryonic feedback

Authors: Benjamin Bose, Bill S. Wright, Matteo Cataneo, Alkistis Pourtsidou, Carlo Giocoli, Lucas Lombriser, Ian G. McCarthy, Marco Baldi, Simon Pfeifer, Qianli Xia

Abstract: In the context of forthcoming galaxy surveys, to ensure unbiased constraints on cosmology and gravity when using non-linear structure information, percent-level accuracy is required when modelling the power spectrum. This calls for frameworks that can accurately capture the relevant physical effects, while allowing for deviations from $Λ$CDM. Massive neutrino and baryonic physics are two of the mo… ▽ More In the context of forthcoming galaxy surveys, to ensure unbiased constraints on cosmology and gravity when using non-linear structure information, percent-level accuracy is required when modelling the power spectrum. This calls for frameworks that can accurately capture the relevant physical effects, while allowing for deviations from $Λ$CDM. Massive neutrino and baryonic physics are two of the most relevant such effects. We present an integration of the halo model reaction frameworks for massive neutrinos and beyond-$Λ$CDM cosmologies. The integrated halo model reaction, combined with a pseudo power spectrum modelled by HMCode2020 is then compared against $N$-body simulations that include both massive neutrinos and an $f(R)$ modification to gravity. We find that the framework is 4% accurate down to at least $k\approx 3 \, h/{\rm Mpc}$ for a modification to gravity of $|f_{\rm R0}|\leq 10^{-5}$ and for the total neutrino mass $M_ν\equiv \sum m_ν\leq 0.15$ eV. We also find that the framework is 4% consistent with EuclidEmulator2 as well as the Bacco emulator for \B{most of the considered} $νw$CDM cosmologies down to at least $k \approx 3 \, h$/Mpc. Finally, we compare against hydrodynamical simulations employing HMCode2020's baryonic feedback modelling on top of the halo model reaction. For $νΛ$CDM cosmologies we find 2% accuracy for $M_ν\leq 0.48$eV down to at least $k\approx 5h$/Mpc. Similar accuracy is found when comparing to $νw$CDM hydrodynamical simulations with $M_ν= 0.06$eV. This offers the first non-linear, theoretically general means of accurately including massive neutrinos for beyond-$Λ$CDM cosmologies, and further suggests that baryonic, massive neutrino and dark energy physics can be reliably modelled independently. △ Less

Submitted 25 September, 2021; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 13 pages, 11 figures. MNRAS accepted version. Download ReACT: https://github.com/nebblu/ReACT/tree/react_with_neutrinos

Euclid preparation: XII. Optimizing the photometric sample of the Euclid survey for galaxy clustering and galaxy-galaxy lensing analyses

Authors: Euclid Collaboration, A. Pocino, I. Tutusaus, F. J. Castander, P. Fosalba, M. Crocce, A. Porredon, S. Camera, V. Cardone, S. Casas, T. Kitching, F. Lacasa, M. Martinelli, A. Pourtsidou, Z. Sakr, S. Andreon, N. Auricchio, C. Baccigalupi, A. Balaguera-Antolínez, M. Baldi, A. Balestra, S. Bardelli, R. Bender, A. Biviano, C. Bodendorf , et al. (135 additional authors not shown)

Abstract: The accuracy of photometric redshifts (photo-zs) particularly affects the results of the analyses of galaxy clustering with photometrically-selected galaxies (GCph) and weak lensing. In the next decade, space missions like Euclid will collect photometric measurements for millions of galaxies. These data should be complemented with upcoming ground-based observations to derive precise and accurate p… ▽ More The accuracy of photometric redshifts (photo-zs) particularly affects the results of the analyses of galaxy clustering with photometrically-selected galaxies (GCph) and weak lensing. In the next decade, space missions like Euclid will collect photometric measurements for millions of galaxies. These data should be complemented with upcoming ground-based observations to derive precise and accurate photo-zs. In this paper, we explore how the tomographic redshift binning and depth of ground-based observations will affect the cosmological constraints expected from Euclid. We focus on GCph and extend the study to include galaxy-galaxy lensing (GGL). We add a layer of complexity to the analysis by simulating several realistic photo-z distributions based on the Euclid Consortium Flagship simulation and using a machine learning photo-z algorithm. We use the Fisher matrix formalism and these galaxy samples to study the cosmological constraining power as a function of redshift binning, survey depth, and photo-z accuracy. We find that bins with equal width in redshift provide a higher Figure of Merit (FoM) than equipopulated bins and that increasing the number of redshift bins from 10 to 13 improves the FoM by 35% and 15% for GCph and its combination with GGL, respectively. For GCph, an increase of the survey depth provides a higher FoM. But the addition of faint galaxies beyond the limit of the spectroscopic training data decreases the FoM due to the spurious photo-zs. When combining both probes, the number density of the sample, which is set by the survey depth, is the main factor driving the variations in the FoM. We conclude that there is more information that can be extracted beyond the nominal 10 tomographic redshift bins of Euclid and that we should be cautious when adding faint galaxies into our sample, since they can degrade the cosmological constraints. △ Less

Submitted 12 April, 2021; originally announced April 2021.

Comments: 21 pages, 13 figures. Abstract abridged

Journal ref: A&A 655, A44 (2021)

The HI intensity mapping bispectrum including observational effects

Authors: Steven Cunnington, Catherine Watkinson, Alkistis Pourtsidou

Abstract: The bispectrum is a 3-point statistic with the potential to provide additional information beyond power spectra analyses of survey datasets. Radio telescopes which broadly survey the 21cm emission from neutral hydrogen (HI) are a promising way to probe LSS and in this work we present an investigation into the HI intensity mapping (IM) bispectrum using simulations. We present a model of the redshif… ▽ More The bispectrum is a 3-point statistic with the potential to provide additional information beyond power spectra analyses of survey datasets. Radio telescopes which broadly survey the 21cm emission from neutral hydrogen (HI) are a promising way to probe LSS and in this work we present an investigation into the HI intensity mapping (IM) bispectrum using simulations. We present a model of the redshift space HI IM bispectrum including observational effects from the radio telescope beam and 21cm foreground contamination. We validate our modelling prescriptions with measurements from robust IM simulations, inclusive of these observational effects. Our foreground simulations include polarisation leakage, on which we use a Principal Component Analysis cleaning method. We also investigate the effects from a non-Gaussian beam including side-lobes. For a MeerKAT-like single-dish IM survey at $z=0.39$, we find that foreground removal causes a 8% reduction in the equilateral bispectrum's signal-to-noise ratio $S/N$, whereas the beam reduces it by 62%. We find our models perform well, generally providing $χ^2_\text{dof}\sim 1$, indicating a good fit to the data. Whilst our focus is on post-reionisation, single-dish IM, our modelling of observational effects, especially foreground removal, can also be relevant to interferometers and reionisation studies. △ Less

Submitted 28 July, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: 19 pages, 14 figures. The main results for telescope beam and foreground effects on the bispectrum can be found in Sec 3.2 and 3.3. Minor changes to text in Version 2 in response to referee report and an amendment to Fig 12 to include more information. Version 2 accepted for publication in MNRAS

HI constraints from the cross-correlation of eBOSS galaxies and Green Bank Telescope intensity maps

Authors: Laura Wolz, Alkistis Pourtsidou, Kiyoshi W. Masui, Tzu-Ching Chang, Julian E. Bautista, Eva-Maria Mueller, Santiago Avila, David Bacon, Will J. Percival, Steven Cunnington, Chris Anderson, Xuelei Chen, Jean-Paul Kneib, Yi-Chao Li, Yu-Wei Liao, Ue-Li Pen, Jeffrey B. Peterson, Graziano Rossi, Donald P. Schneider, Jaswant Yadav, Gong-Bo Zhao

Abstract: We present the joint analysis of Neutral Hydrogen (HI) Intensity Mapping observations with three galaxy samples: the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples from the eBOSS survey, and the WiggleZ Dark Energy Survey sample. The HI intensity maps are Green Bank Telescope observations of the redshifted 21cm emission on 100deg2 covering the redshift range $0.6<z<1.0$. We proce… ▽ More We present the joint analysis of Neutral Hydrogen (HI) Intensity Mapping observations with three galaxy samples: the Luminous Red Galaxy (LRG) and Emission Line Galaxy (ELG) samples from the eBOSS survey, and the WiggleZ Dark Energy Survey sample. The HI intensity maps are Green Bank Telescope observations of the redshifted 21cm emission on 100deg2 covering the redshift range $0.6<z<1.0$. We process the data by separating and removing the foregrounds with FastICA, and construct a transfer function to correct for the effects of foreground removal on the HI signal. We cross-correlate the cleaned HI data with the galaxy samples and study the overall amplitude as well as the scale-dependence of the power spectrum. We also qualitatively compare our findings with the predictions by a semi-analytic galaxy evolution simulation. The cross-correlations constrain the quantity $Ω_{HI} b_{HI} r_{{HI},{opt}}$ at an effective scale $k_{eff}$, where $Ω_{HI}$ is the HI density fraction, $b_{HI}$ is the HI bias, and $r_{{HI},{opt}}$ the galaxy-hydrogen correlation coefficient, which is dependent on the HI content of the optical galaxy sample. At $k_{eff}=0.31 \, h/{Mpc}$ we find $Ω_{HI} b_{HI} r_{{HI},{Wig}} = [0.58 \pm 0.09 \, {(stat) \pm 0.05 \, {(sys)}}] \times 10^{-3}$ for GBT-WiggleZ, $Ω_{HI} b_{HI} r_{HI,{ELG}} = [0.40 \pm 0.09 \, {(stat) \pm 0.04 \, {(sys)}}] \times 10^{-3}$ for GBT-ELG, and $Ω_{HI} b_{HI} r_{{HI},{LRG}} = [0.35 \pm 0.08 \, {(stat) \pm 0.03 \, {(sys)}}] \times 10^{-3}$ for GBT-LRG, at $z\simeq 0.8$. We also report results at $k_{eff}=0.24 \, h/{Mpc}$ and $k_{eff}=0.48 \, h/{Mpc}$. With little information on HI parameters beyond our local Universe, these are amongst the most precise constraints on neutral hydrogen density fluctuations in an underexplored redshift range. △ Less

Submitted 9 December, 2021; v1 submitted 9 February, 2021; originally announced February 2021.

Comments: 18 pages, 19 figures, 1 table; accepted by MNRAS

Euclid: Forecasts for $k$-cut $3 \times 2$ Point Statistics

Authors: P. L. Taylor, T. Kitching, V. F. Cardone, A. Ferté, E. M. Huff, F. Bernardeau, J. Rhodes, A. C. Deshpande, I. Tutusaus, A. Pourtsidou, S. Camera, C. Carbone, S. Casas, M. Martinelli, V. Pettorino, Z. Sakr, D. Sapone, V. Yankelevich, N. Auricchio, A. Balestra, C. Bodendorf, D. Bonino, A. Boucaud, E. Branchini, M. Brescia , et al. (70 additional authors not shown)

Abstract: Modelling uncertainties at small scales, i.e. high $k$ in the power spectrum $P(k)$, due to baryonic feedback, nonlinear structure growth and the fact that galaxies are biased tracers poses a significant obstacle to fully leverage the constraining power of the {\it Euclid} wide-field survey. $k$-cut cosmic shear has recently been proposed as a method to optimally remove sensitivity to these scales… ▽ More Modelling uncertainties at small scales, i.e. high $k$ in the power spectrum $P(k)$, due to baryonic feedback, nonlinear structure growth and the fact that galaxies are biased tracers poses a significant obstacle to fully leverage the constraining power of the {\it Euclid} wide-field survey. $k$-cut cosmic shear has recently been proposed as a method to optimally remove sensitivity to these scales while preserving usable information. In this paper we generalise the $k$-cut cosmic shear formalism to $3 \times 2$ point statistics and estimate the loss of information for different $k$-cuts in a $3 \times 2$ point analysis of the {\it Euclid} data. Extending the Fisher matrix analysis of~\citet{blanchard2019euclid}, we assess the degradation in constraining power for different $k$-cuts. We work in the idealised case and assume the galaxy bias is linear, the covariance is Gaussian, while neglecting uncertainties due to photo-z errors and baryonic feedback. We find that taking a $k$-cut at $2.6 \ h \ {\rm Mpc} ^{-1}$ yields a dark energy Figure of Merit (FOM) of 1018. This is comparable to taking a weak lensing cut at $\ell = 5000$ and a galaxy clustering and galaxy-galaxy lensing cut at $\ell = 3000$ in a traditional $3 \times 2$ point analysis. We also find that the fraction of the observed galaxies used in the photometric clustering part of the analysis is one of the main drivers of the FOM. Removing $50 \% \ (90 \%)$ of the clustering galaxies decreases the FOM by $19 \% \ (62 \%)$. Given that the FOM depends so heavily on the fraction of galaxies used in the clustering analysis, extensive efforts should be made to handle the real-world systematics present when extending the analysis beyond the luminous red galaxy (LRG) sample. △ Less

Submitted 20 July, 2021; v1 submitted 8 December, 2020; originally announced December 2020.

Comments: 10 pages, 5 figures. Accepted by the Open Journal of Astrophysics

HI intensity mapping with MeerKAT: Calibration pipeline for multi-dish autocorrelation observations

Authors: Jingying Wang, Mario G. Santos, Philip Bull, Keith Grainge, Steven Cunnington, Jose Fonseca, Melis O. Irfan, Yichao Li, Alkistis Pourtsidou, Paula S. Soares, Marta Spinelli, Gianni Bernardi, Brandon Engelbrecht

Abstract: While most purpose-built 21cm intensity mapping experiments are close-packed interferometer arrays, general-purpose dish arrays should also be capable of measuring the cosmological 21cm signal. This can be achieved most efficiently if the array is used as a collection of scanning autocorrelation dishes rather than as an interferometer. As a first step towards demonstrating the feasibility of this… ▽ More While most purpose-built 21cm intensity mapping experiments are close-packed interferometer arrays, general-purpose dish arrays should also be capable of measuring the cosmological 21cm signal. This can be achieved most efficiently if the array is used as a collection of scanning autocorrelation dishes rather than as an interferometer. As a first step towards demonstrating the feasibility of this observing strategy, we show that we are able to successfully calibrate dual-polarisation autocorrelation data from 64 MeerKAT dishes in the L-band (856-1712 MHz, 4096 channels), with 10.5 hours of data retained from six nights of observing. We describe our calibration pipeline, which is based on multi-level RFI flagging, periodic noise diode injection to stabilise gain drifts and an absolute calibration based on a multi-component sky model. We show that it is sufficiently accurate to recover maps of diffuse celestial emission and point sources over a 10 deg x 30 deg patch of the sky overlapping with the WiggleZ 11hr field. The reconstructed maps have a good level of consistency between per-dish maps and external datasets, with the estimated thermal noise limited to 1.4 x the theoretical noise level (~ 2 mK). The residual maps have rms amplitudes below 0.1 K, corresponding to <1% of the model temperature. The reconstructed Galactic HI intensity map shows excellent agreement with the Effelsberg-Bonn HI Survey, and the flux of the radio galaxy 4C+03.18 is recovered to within 3.6%, which demonstrates that the autocorrelation can be successfully calibrated to give the zero-spacing flux and potentially help in the imaging of MeerKAT interferometric data. Our results provide a positive indication towards the feasibility of using MeerKAT and the future SKA to measure the HI intensity mapping signal and probe cosmology on degree scales and above. △ Less

Submitted 13 May, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 24 pages, 34 figures, accepted for publication in MNRAS

Euclid: impact of nonlinear prescriptions on cosmological parameter estimation from weak lensing cosmic shear

Authors: M. Martinelli, I. Tutusaus, M. Archidiacono, S. Camera, V. F. Cardone, S. Clesse, S. Casas, L. Casarini, D. F. Mota, H. Hoekstra, C. Carbone, S. Ilić, T. D. Kitching, V. Pettorino, A. Pourtsidou, Z. Sakr, D. Sapone, N. Auricchio, A. Balestra, A. Boucaud, E. Branchini, M. Brescia, V. Capobianco, J. Carretero, M. Castellano , et al. (69 additional authors not shown)

Abstract: Upcoming surveys will map the growth of large-scale structure with unprecented precision, improving our understanding of the dark sector of the Universe. Unfortunately, much of the cosmological information is encoded by the small scales, where the clustering of dark matter and the effects of astrophysical feedback processes are not fully understood. This can bias the estimates of cosmological para… ▽ More Upcoming surveys will map the growth of large-scale structure with unprecented precision, improving our understanding of the dark sector of the Universe. Unfortunately, much of the cosmological information is encoded by the small scales, where the clustering of dark matter and the effects of astrophysical feedback processes are not fully understood. This can bias the estimates of cosmological parameters, which we study here for a joint analysis of mock Euclid cosmic shear and Planck cosmic microwave background data. We use different implementations for the modelling of the signal on small scales and find that they result in significantly different predictions. Moreover, the different nonlinear corrections lead to biased parameter estimates, especially when the analysis is extended into the highly nonlinear regime, with both the Hubble constant, $H_0$, and the clustering amplitude, $σ_8$, affected the most. Improvements in the modelling of nonlinear scales will therefore be needed if we are to resolve the current tension with more and better data. For a given prescription for the nonlinear power spectrum, using different corrections for baryon physics does not significantly impact the precision of Euclid, but neglecting these correction does lead to large biases in the cosmological parameters. In order to extract precise and unbiased constraints on cosmological parameters from Euclid cosmic shear data, it is therefore essential to improve the accuracy of the recipes that account for nonlinear structure formation, as well as the modelling of the impact of astrophysical processes that redistribute the baryons. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: 18 pages, 6 figures. Prepared for submission to Astronomy & Astrophysics

Report number: IFT-UAM/CSIC-20-146

21cm foregrounds and polarization leakage: a user's guide on cleaning and mitigation strategies

Authors: Steven Cunnington, Melis O. Irfan, Isabella P. Carucci, Alkistis Pourtsidou, Jérôme Bobin

Abstract: The success of HI intensity mapping is largely dependent on how well 21cm foreground contamination can be controlled. In order to progress our understanding further, we present a range of simulated foreground data from four different $\sim3000$\,deg$^2$ sky regions, with and without effects from polarization leakage. Combining these with underlying cosmological HI simulations creates a range of si… ▽ More The success of HI intensity mapping is largely dependent on how well 21cm foreground contamination can be controlled. In order to progress our understanding further, we present a range of simulated foreground data from four different $\sim3000$\,deg$^2$ sky regions, with and without effects from polarization leakage. Combining these with underlying cosmological HI simulations creates a range of single-dish intensity mapping test cases that require different foreground treatments. This allows us to conduct the most generalized study to date into 21cm foregrounds and their cleaning techniques for the post-reionization era. We first provide a pedagogical review of the most commonly used blind foreground removal techniques (PCA/SVD, FASTICA, GMCA). We also trial a non-blind parametric fitting technique and discuss potential hybridization of methods. We highlight the similarities and differences in these techniques finding that the blind methods produce near equivalent results, and we explain the fundamental reasons for this. The simulations allow an exact decomposition of the resulting cleaned data and we analyse the contribution from foreground residuals. Our results demonstrate that polarized foreground residuals should be generally subdominant to HI on small scales ($k\gtrsim0.1\,h\,\text{Mpc}^{-1}$). However, on larger scales, results are more region dependent. In some cases, aggressive cleans severely damp HI power but still leave dominant foreground residuals. We also demonstrate the gain from cross-correlations with optical galaxy surveys, where extreme levels of residual foregrounds can be circumvented. However, these residuals still contribute to errors and we discuss the optimal balance between over- and under-cleaning. △ Less

Submitted 24 March, 2021; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 21 pages, 16 figures. Changes made in response to reviewer report. Noted changes from V1: analysis now only performed on two foreground regions (excluding Galactic Plane and Barnard's loop since these were very high emission regions). General results and conclusions unchanged though. Version 2 accepted for publication in MNRAS

Power spectrum multipole expansion for HI intensity mapping experiments: unbiased parameter estimation

Authors: Paula S. Soares, Steven Cunnington, Alkistis Pourtsidou, Chris Blake

Abstract: We assess the performance of the multipole expansion formalism in the case of single-dish HI intensity mapping, including instrumental and foreground removal effects. This formalism is used to provide MCMC forecasts for a range of HI and cosmological parameters, including redshift space distortions and the Alcock-Paczynski effect. We first determine the range of validity of our power spectrum mode… ▽ More We assess the performance of the multipole expansion formalism in the case of single-dish HI intensity mapping, including instrumental and foreground removal effects. This formalism is used to provide MCMC forecasts for a range of HI and cosmological parameters, including redshift space distortions and the Alcock-Paczynski effect. We first determine the range of validity of our power spectrum modelling by fitting to simulation data, concentrating on the monopole, quadrupole, and hexadecapole contributions. We then show that foreground subtraction effects can lead to severe biases in the determination of cosmological parameters, in particular the parameters relating to the transverse BAO rescaling, the growth rate and the HI bias ($α_\perp$, $\overline{T}_\text{HI} fσ_8$, and $\overline{T}_\text{HI} b_\text{HI} σ_8$, respectively). We attempt to account for these biases by constructing a 2-parameter foreground modelling prescription, and find that our prescription leads to unbiased parameter estimation at the expense of increasing the estimated uncertainties on cosmological parameters. In addition, we confirm that instrumental and foreground removal effects significantly impact the theoretical covariance matrix, and cause the covariance between different multipoles to become non-negligible. Finally, we show the effect of including higher-order multipoles in our analysis, and how these can be used to investigate the presence of instrumental and systematic effects in HI intensity mapping data. △ Less

Submitted 5 January, 2021; v1 submitted 27 August, 2020; originally announced August 2020.

Comments: 18 pages, 12 figures, 5 tables, main concepts in figures 2 and 3, main results in figures 5 and 10

The degeneracy between primordial non-Gaussianity and foregrounds in 21cm intensity mapping experiments

Authors: Steven Cunnington, Stefano Camera, Alkistis Pourtsidou

Abstract: Potential evidence for primordial non-Gaussianity (PNG) is expected to lie in the largest scales mapped by cosmological surveys. Forthcoming 21cm intensity mapping experiments will aim to probe these scales by surveying neutral hydrogen (HI) within galaxies. However, foreground signals dominate the faint 21cm emission, meaning foreground cleaning is required to recover the cosmological signal. The… ▽ More Potential evidence for primordial non-Gaussianity (PNG) is expected to lie in the largest scales mapped by cosmological surveys. Forthcoming 21cm intensity mapping experiments will aim to probe these scales by surveying neutral hydrogen (HI) within galaxies. However, foreground signals dominate the faint 21cm emission, meaning foreground cleaning is required to recover the cosmological signal. The effect this has is to damp the HI power spectrum on the largest scales, especially along the line-of-sight. Whilst there is agreement that this contamination is potentially problematic for probing PNG, it is yet to be fully explored and quantified. In this work we carry out the first forecasts on $f_\text{NL}$ that incorporate simulated foreground maps that are removed using techniques employed in real data. Using an MCMC analysis, we demonstrate that foreground cleaned data recovers hugely biased values ($f_\text{NL} = -102.1_{-7.96}^{+8.39}$ [68% CL]) on our $f_\text{NL}=0$ fiducial input. Introducing a model with fixed parameters for the foreground contamination allows us to recover unbiased results ($f_\text{NL} = -2.94_{-11.9}^{+11.4}$). However, it is not clear that we will have sufficient understanding of foreground contamination to allow for such rigid models. Treating the main parameter $k_\parallel^\text{FG}$ in our foreground model as a nuisance parameter and marginalizing over it, still recovers unbiased results but at the expense of much larger errors ($f_\text{NL} = 0.75^{+40.2}_{-44.5}$), that can only be reduced by imposing the Planck 2018 prior. Our results show that significant progress on understanding and controlling foreground removal effects is necessary in order to study PNG with HI intensity mapping. △ Less

Submitted 19 October, 2020; v1 submitted 23 July, 2020; originally announced July 2020.

Comments: 15 pages, 8 figures. The busy reader should study Fig 8 which demonstrates the degeneracy between f_NL and k_FG the main parameter in our model to account for foreground removal effects. Accepted by MNRAS