-
Euclid preparation. Sensitivity to non-standard particle dark matter model
Authors:
Euclid Collaboration,
J. Lesgourgues,
J. Schwagereit,
J. Bucko,
G. Parimbelli,
S. K. Giri,
F. Hervas-Peters,
A. Schneider,
M. Archidiacono,
F. Pace,
Z. Sakr,
A. Amara,
L. Amendola,
S. Andreon,
N. Auricchio,
H. Aussel,
C. Baccigalupi,
M. Baldi,
S. Bardelli,
R. Bender,
C. Bodendorf,
D. Bonino,
E. Branchini,
M. Brescia,
J. Brinchmann
, et al. (227 additional authors not shown)
Abstract:
The Euclid mission of the European Space Agency will provide weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and its extensions, with an opportunity to test the properties of dark matter beyond the minimal cold dark matter paradigm. We present forecasts from the combination of these surveys on the parameters describing four int…
▽ More
The Euclid mission of the European Space Agency will provide weak gravitational lensing and galaxy clustering surveys that can be used to constrain the standard cosmological model and its extensions, with an opportunity to test the properties of dark matter beyond the minimal cold dark matter paradigm. We present forecasts from the combination of these surveys on the parameters describing four interesting and representative non-minimal dark matter models: a mixture of cold and warm dark matter relics; unstable dark matter decaying either into massless or massive relics; and dark matter experiencing feeble interactions with relativistic relics. We model these scenarios at the level of the non-linear matter power spectrum using emulators trained on dedicated N-body simulations. We use a mock Euclid likelihood to fit mock data and infer error bars on dark matter parameters marginalised over other parameters. We find that the Euclid photometric probe (alone or in combination with CMB data from the Planck satellite) will be sensitive to the effect of each of the four dark matter models considered here. The improvement will be particularly spectacular for decaying and interacting dark matter models. With Euclid, the bounds on some dark matter parameters can improve by up to two orders of magnitude compared to current limits. We discuss the dependence of predicted uncertainties on different assumptions: inclusion of photometric galaxy clustering data, minimum angular scale taken into account, modelling of baryonic feedback effects. We conclude that the Euclid mission will be able to measure quantities related to the dark sector of particle physics with unprecedented sensitivity. This will provide important information for model building in high-energy physics. Any hint of a deviation from the minimal cold dark matter paradigm would have profound implications for cosmology and particle physics.
△ Less
Submitted 26 June, 2024;
originally announced June 2024.
-
Euclid preparation. Observational expectations for redshift z<7 active galactic nuclei in the Euclid Wide and Deep surveys
Authors:
Euclid Collaboration,
M. Selwood,
S. Fotopoulou,
M. N. Bremer,
L. Bisigello,
H. Landt,
E. Bañados,
G. Zamorani,
F. Shankar,
D. Stern,
E. Lusso,
L. Spinoglio,
V. Allevato,
F. Ricci,
A. Feltre,
F. Mannucci,
M. Salvato,
R. A. A. Bowler,
M. Mignoli,
D. Vergani,
F. La Franca,
A. Amara,
S. Andreon,
N. Auricchio,
M. Baldi
, et al. (238 additional authors not shown)
Abstract:
We forecast the expected population of active galactic nuclei (AGN) observable in the Euclid Wide Survey (EWS) and Euclid Deep Survey (EDS). Starting from an X-ray luminosity function (XLF) we generate volume-limited samples of the AGN expected in the survey footprints. Each AGN is assigned an SED appropriate for its X-ray luminosity and redshift, with perturbations sampled from empirical distribu…
▽ More
We forecast the expected population of active galactic nuclei (AGN) observable in the Euclid Wide Survey (EWS) and Euclid Deep Survey (EDS). Starting from an X-ray luminosity function (XLF) we generate volume-limited samples of the AGN expected in the survey footprints. Each AGN is assigned an SED appropriate for its X-ray luminosity and redshift, with perturbations sampled from empirical distributions. The photometric detectability of each AGN is assessed via mock observation of the assigned SED. We estimate 40 million AGN will be detectable in at least one band in the EWS and 0.24 million in the EDS, corresponding to surface densities of 2.8$\times$10$^{3}$ deg$^{-2}$ and 4.7$\times$10$^{3}$ deg$^{-2}$. Employing colour selection criteria on our simulated data we select a sample of 4.8$\times$10$^{6}$ (331 deg$^{-2}$) AGN in the EWS and 1.7$\times$10$^{4}$ (346 deg$^{-2}$) in the EDS, amounting to 10% and 8% of the AGN detectable in the EWS and EDS. Including ancillary Rubin/LSST bands improves the completeness and purity of AGN selection. These data roughly double the total number of selected AGN to comprise 21% and 15% of the detectable AGN in the EWS and EDS. The total expected sample of colour-selected AGN contains 6.0$\times$10$^{6}$ (74%) unobscured AGN and 2.1$\times$10$^{6}$ (26%) obscured AGN, covering $0.02 \leq z \lesssim 5.2$ and $43 \leq \log_{10} (L_{bol} / erg s^{-1}) \leq 47$. With this simple colour selection, expected surface densities are already comparable to the yield of modern X-ray and mid-infrared surveys of similar area. The relative uncertainty on our expectation for detectable AGN is 6.7% for the EWS and 12.5% for the EDS, driven by the uncertainty of the XLF.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
Euclid preparation. Detecting globular clusters in the Euclid survey
Authors:
Euclid Collaboration,
K. Voggel,
A. Lançon,
T. Saifollahi,
S. S. Larsen,
M. Cantiello,
M. Rejkuba,
J. -C. Cuillandre,
P. Hudelot,
A. A. Nucita,
M. Urbano,
E. Romelli,
M. A. Raj,
M. Schirmer,
C. Tortora,
Abdurro'uf,
F. Annibali,
M. Baes,
P. Boldrini,
R. Cabanac,
D. Carollo,
C. J. Conselice,
P. -A. Duc,
A. M. N. Ferguson,
L. K. Hunt
, et al. (247 additional authors not shown)
Abstract:
Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid…
▽ More
Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics.
△ Less
Submitted 29 May, 2024; v1 submitted 22 May, 2024;
originally announced May 2024.
-
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.
-
Euclid. IV. The NISP Calibration Unit
Authors:
Euclid Collaboration,
F. Hormuth,
K. Jahnke,
M. Schirmer,
C. G. -Y. Lee,
T. Scott,
R. Barbier,
S. Ferriol,
W. Gillard,
F. Grupp,
R. Holmes,
W. Holmes,
B. Kubik,
J. Macias-Perez,
M. Laurent,
J. Marpaud,
M. Marton,
E. Medinaceli,
G. Morgante,
R. Toledo-Moreo,
M. Trifoglio,
Hans-Walter Rix,
A. Secroun,
M. Seiffert,
P. Stassi
, et al. (310 additional authors not shown)
Abstract:
The near-infrared calibration unit (NI-CU) onboard Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid is a mission in ESA's 'Cosmic Vision 2015-2025' framework, to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and d…
▽ More
The near-infrared calibration unit (NI-CU) onboard Euclid's Near-Infrared Spectrometer and Photometer (NISP) is the first astronomical calibration lamp based on light-emitting diodes (LEDs) to be operated in space. Euclid is a mission in ESA's 'Cosmic Vision 2015-2025' framework, to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and dark energy. Calibrating photometric and spectrometric measurements of galaxies to better than 1.5% accuracy in a survey homogeneously mapping ~14000 deg^2 of extragalactic sky requires a very detailed characterisation of near-infrared (NIR) detector properties, as well their constant monitoring in flight. To cover two of the main contributions - relative pixel-to-pixel sensitivity and non-linearity characteristics - as well as support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous (<12% variations) and temporally stable illumination (0.1%-0.2% over 1200s) over the NISP detector plane, with minimal power consumption and energy dissipation. NI-CU is covers the spectral range ~[900,1900] nm - at cryo-operating temperature - at 5 fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of >=100 from ~15 ph s^-1 pixel^-1 to >1500 ph s^-1 pixel^-1. For this functionality, NI-CU is based on LEDs. We describe the rationale behind the decision and design process, describe the challenges in sourcing the right LEDs, as well as the qualification process and lessons learned. We also provide a description of the completed NI-CU, its capabilities and performance as well as its limits. NI-CU has been integrated into NISP and the Euclid satellite, and since Euclid's launch in July 2023 has started supporting survey operations.
△ Less
Submitted 22 May, 2024;
originally announced May 2024.
-
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.
-
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.
-
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.
-
Euclid preparation. Improving cosmological constraints using a new multi-tracer method with the spectroscopic and photometric samples
Authors:
Euclid Collaboration,
F. Dournac,
A. Blanchard,
S. Ilić,
B. Lamine,
I. Tutusaus,
A. Amara,
S. Andreon,
N. Auricchio,
H. Aussel,
M. Baldi,
S. Bardelli,
C. Bodendorf,
D. Bonino,
E. Branchini,
S. Brau-Nogue,
M. Brescia,
J. Brinchmann,
S. Camera,
V. Capobianco,
J. Carretero,
S. Casas,
M. Castellano,
S. Cavuoti,
A. Cimatti
, et al. (218 additional authors not shown)
Abstract:
Future data provided by the \Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of \Euclid, namely the three-dimensional clustering of the spectroscopic galaxy samp…
▽ More
Future data provided by the \Euclid mission will allow us to better understand the cosmic history of the Universe. A metric of its performance is the figure-of-merit (FoM) of dark energy, usually estimated with Fisher forecasts. The expected FoM has previously been estimated taking into account the two main probes of \Euclid, namely the three-dimensional clustering of the spectroscopic galaxy sample, and the so-called 3$\times$2\,pt signal from the photometric sample (i.e., the weak lensing signal, the galaxy clustering, and their cross-correlation). So far, these two probes have been treated as independent. In this paper, we introduce a new observable given by the ratio of the (angular) two-point correlation function of galaxies from the two surveys. For identical (normalised) selection functions, this observable is unaffected by sampling noise, and its variance is solely controlled by Poisson noise. We present forecasts for \Euclid where this multi-tracer method is applied and is particularly relevant because the two surveys will cover the same area of the sky. This method allows for the exploitation of the combination of the spectroscopic and photometric samples. When the correlation between this new observable and the other probes is not taken into account, a significant gain is obtained in the FoM, as well as in the constraints on other cosmological parameters. The benefit is more pronounced for a commonly investigated modified gravity model, namely the $γ$ parametrisation of the growth factor. However, the correlation between the different probes is found to be significant and hence the actual gain is uncertain. We present various strategies for circumventing this issue and still extract useful information from the new observable.
△ Less
Submitted 18 April, 2024;
originally announced April 2024.
-
Non-linear Eulerian Hydrodynamics of Dark Energy: Riemann problem and Finite Volume Schemes
Authors:
Linda Blot,
Pier Stefano Corasaniti,
Fabian Schmidt
Abstract:
Upcoming large-scale-structure surveys can shed new light on the properties of dark energy. In particular, if dark energy is a dynamical component, it must have spatial perturbations. Their behaviour is regulated by the speed of sound parameter, which is currently unconstrained. In this work we present the numerical methods that will allow to perform cosmological simulations of inhomogeneous dark…
▽ More
Upcoming large-scale-structure surveys can shed new light on the properties of dark energy. In particular, if dark energy is a dynamical component, it must have spatial perturbations. Their behaviour is regulated by the speed of sound parameter, which is currently unconstrained. In this work we present the numerical methods that will allow to perform cosmological simulations of inhomogeneous dark energy scenarios where the speed of sound is small and non-vanishing. We treat the dark energy component as an effective fluid and build upon established numerical methods for hydrodynamics to construct a numerical solution of the effective continuity and Euler equations. In particular, we develop conservative finite volume schemes that rely on the solution of the Riemann problem, which we provide here in both exact and approximate forms for the case of a dark energy fluid.
△ Less
Submitted 21 April, 2023; v1 submitted 10 October, 2022;
originally announced October 2022.
-
Euclid: Fast two-point correlation function covariance through linear construction
Authors:
E. Keihanen,
V. Lindholm,
P. Monaco,
L. Blot,
C. Carbone,
K. Kiiveri,
A. G. Sánchez,
A. Viitanen,
J. Valiviita,
A. Amara,
N. Auricchio,
M. Baldi,
D. Bonino,
E. Branchini,
M. Brescia,
J. Brinchmann,
S. Camera,
V. Capobianco,
J. Carretero,
M. Castellano,
S. Cavuoti,
A. Cimatti,
R. Cledassou,
G. Congedo,
L. Conversi
, et al. (87 additional authors not shown)
Abstract:
We present a method for fast evaluation of the covariance matrix for a two-point galaxy correlation function (2PCF) measured with the Landy-Szalay estimator. The standard way of evaluating the covariance matrix consists in running the estimator on a large number of mock catalogs, and evaluating their sample covariance. With large random catalog sizes (data-to-random objects ratio M>>1) the computa…
▽ More
We present a method for fast evaluation of the covariance matrix for a two-point galaxy correlation function (2PCF) measured with the Landy-Szalay estimator. The standard way of evaluating the covariance matrix consists in running the estimator on a large number of mock catalogs, and evaluating their sample covariance. With large random catalog sizes (data-to-random objects ratio M>>1) the computational cost of the standard method is dominated by that of counting the data-random and random-random pairs, while the uncertainty of the estimate is dominated by that of data-data pairs. We present a method called Linear Construction (LC), where the covariance is estimated for small random catalogs of size M = 1 and M = 2, and the covariance for arbitrary M is constructed as a linear combination of these. We validate the method with PINOCCHIO simulations in range r = 20-200 Mpc/h, and show that the covariance estimate is unbiased. With M = 50 and with 2 Mpc/h bins, the theoretical speed-up of the method is a factor of 14. We discuss the impact on the precision matrix and parameter estimation, and derive a formula for the covariance of covariance.
△ Less
Submitted 24 May, 2022;
originally announced May 2022.
-
Dark Energy Survey Year 3 Results: Galaxy mock catalogs for BAO analysis
Authors:
I. Ferrero,
M. Crocce,
I. Tutusaus,
A. Porredon,
L. Blot,
P. Fosalba,
A. Carnero Rosell,
S. Avila,
A. Izard,
J. Elvin-Poole,
K. C. Chan,
H. Camacho,
R. Rosenfeld,
E. Sanchez,
P. Tallada-Crespí,
J. Carretero,
I. Sevilla-Noarbe,
E. Gaztanaga,
F. Andrade-Oliveira,
J. De Vicente,
J. Mena-Fernández,
A. J. Ross,
D. Sanchez Cid,
A. Ferté,
A. Brandao-Souza
, et al. (61 additional authors not shown)
Abstract:
The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the D…
▽ More
The calibration and validation of scientific analysis in simulations is a fundamental tool to ensure unbiased and robust results in observational cosmology. In particular, mock galaxy catalogs are a crucial resource to achieve these goals in the measurement of baryon acoustic oscillation (BAO) in the clustering of galaxies. Here we present a set of 1952 galaxy mock catalogs designed to mimic the Dark Energy Survey (DES) Year 3 BAO sample over its full photometric redshift range $0.6<z_{\rm photo}<1.1$. The mocks are based upon 488 ICE-COLA fast $N$-body simulations of full-sky light cones and were created by populating halos with galaxies, using a hybrid halo occupation distribution - halo abundance matching model. This model has ten free parameters, which were determined, for the first time, using an automatic likelihood minimization procedure. We also introduced a novel technique to assign photometric redshift for simulated galaxies, following a two-dimensional probability distribution with VIMOS Public Extragalactic Redshift Survey (VIPERS) data. The calibration was designed to match the observed abundance of galaxies as a function of photometric redshift, the distribution of photometric redshift errors, and the clustering amplitude on scales smaller than those used for BAO measurements. An exhaustive analysis was done to ensure that the mocks reproduce the input properties. Finally, mocks were tested by comparing the angular correlation function $w(θ)$, angular power spectrum $C_\ell$, and projected clustering $ξ_p(r_\perp)$ to theoretical predictions and data. The impact of volume replication in the estimate of the covariance is also investigated. The success in accurately reproducing the photometric redshift uncertainties and the galaxy clustering as a function of redshift render this mock creation pipeline as a benchmark for future analyses of photometric galaxy surveys.
△ Less
Submitted 13 December, 2021; v1 submitted 9 July, 2021;
originally announced July 2021.
-
Cosmological Model Parameter Dependence of the Matter Power Spectrum Covariance from the DEUS-PUR $Cosmo$ Simulations
Authors:
Linda Blot,
Pier-Stefano Corasaniti,
Yann Rasera,
Shankar Agarwal
Abstract:
Future galaxy surveys will provide accurate measurements of the matter power spectrum across an unprecedented range of scales and redshifts. The analysis of these data will require one to accurately model the imprint of non-linearities of the matter density field. In particular, these induce a non-Gaussian contribution to the data covariance that needs to be properly taken into account to realise…
▽ More
Future galaxy surveys will provide accurate measurements of the matter power spectrum across an unprecedented range of scales and redshifts. The analysis of these data will require one to accurately model the imprint of non-linearities of the matter density field. In particular, these induce a non-Gaussian contribution to the data covariance that needs to be properly taken into account to realise unbiased cosmological parameter inference analyses. Here, we study the cosmological dependence of the matter power spectrum covariance using a dedicated suite of N-body simulations, the Dark Energy Universe Simulation - Parallel Universe Runs (DEUS-PUR) {\it Cosmo}. These consist of 512 realizations for 10 different cosmologies where we vary the matter density $Ω_m$, the amplitude of density fluctuations $σ_8$, the reduced Hubble parameter $h$ and a constant dark energy equation of state $w$ by approximately $10\%$. We use these data to evaluate the first and second derivatives of the power spectrum covariance with respect to a fiducial $Λ$CDM cosmology. We find that the variations can be as large as $150\%$ depending on the scale, redshift and model parameter considered. By performing a Fisher matrix analysis we explore the impact of different choices in modelling the cosmological dependence of the covariance. Our results suggest that fixing the covariance to a fiducial cosmology can significantly affect the recovered parameter errors and that modelling the cosmological dependence of the variance while keeping the correlation coefficient fixed can alleviate the impact of this effect.
△ Less
Submitted 2 November, 2020; v1 submitted 29 July, 2020;
originally announced July 2020.
-
Euclid preparation: II. The EuclidEmulator -- A tool to compute the cosmology dependence of the nonlinear matter power spectrum
Authors:
Euclid Collaboration,
Mischa Knabenhans,
Joachim Stadel,
Stefano Marelli,
Doug Potter,
Romain Teyssier,
Laurent Legrand,
Aurel Schneider,
Bruno Sudret,
Linda Blot,
Saeeda Awan,
Carlo Burigana,
Carla Sofia Carvalho,
Hannu Kurki-Suonio,
Gabriele Sirri
Abstract:
We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum. It is based on a spectral decomposition method called polynomial chaos expansion. All steps in the construction of the emulator have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a si…
▽ More
We present a new power spectrum emulator named EuclidEmulator that estimates the nonlinear correction to the linear dark matter power spectrum. It is based on a spectral decomposition method called polynomial chaos expansion. All steps in the construction of the emulator have been tested and optimized: the large high-resolution N-body simulations carried out with PKDGRAV3 were validated using a simulation from the Euclid Flagship campaign and demonstrated to have converged up to wavenumbers $k\approx 5\,h\,{\rm Mpc}^{-1}$ for redshifts $z\leq 5$. The emulator is constructed using the uncertainty quantification software UQLab and it has been optimized first by creating mock emulators based on Takahashi's HALOFIT. We show that it is possible to successfully predict the performance of the final emulator in this way prior to performing any N-body simulations. We provide a C-code to calculate the nonlinear correction at a relative accuracy of $\sim0.3\%$ with respect to N-body simulations within 50 ms. The absolute accuracy of the final nonlinear power spectrum is comparable to one obtained with N-body simulations, i.e. $\sim 1\%$ for $k\lesssim 1\,h\,{\rm Mpc}^{-1}$ and $z\lesssim 3.5$. This enables efficient forward modeling in the nonlinear regime allowing for maximum likelihood estimation of cosmological parameters. EuclidEmulator has been compared to HALOFIT and CosmicEmu, an alternative emulator based on the Mira-Titan Universe, and shown to be more accurate than these other approaches. This work paves a new way for optimal construction of future emulators that also consider other cosmological observables, use higher resolution input simulations and investigate higher dimensional cosmological parameter spaces.
△ Less
Submitted 20 February, 2019; v1 submitted 12 September, 2018;
originally announced September 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices III: Bispectrum
Authors:
Manuel Colavincenzo,
Emiliano Sefusatti,
Pierluigi Monaco,
Linda Blot,
Martin Crocce,
Martha Lippich,
Ariel G. Sánchez,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in…
▽ More
We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.
△ Less
Submitted 8 October, 2018; v1 submitted 25 June, 2018;
originally announced June 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices II: Power spectrum multipoles
Authors:
Linda Blot,
Martin Crocce,
Emiliano Sefusatti,
Martha Lippich,
Ariel G. Sánchez,
Manuel Colavincenzo,
Pierluigi Monaco,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo densi…
▽ More
We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.
△ Less
Submitted 18 February, 2019; v1 submitted 25 June, 2018;
originally announced June 2018.
-
Comparing approximate methods for mock catalogues and covariance matrices I: correlation function
Authors:
Martha Lippich,
Ariel G. Sánchez,
Manuel Colavincenzo,
Emiliano Sefusatti,
Pierluigi Monaco,
Linda Blot,
Martin Crocce,
Marcelo A. Alvarez,
Aniket Agrawal,
Santiago Avila,
Andrés Balaguera-Antolínez,
Richard Bond,
Sandrine Codis,
Claudio Dalla Vecchia,
Antonio Dorta,
Pablo Fosalba,
Albert Izard,
Francisco-Shu Kitaura,
Marcos Pellejero-Ibanez,
George Stein,
Mohammadjavad Vakili,
Gustavo Yepes
Abstract:
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that fol…
▽ More
This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.
△ Less
Submitted 13 May, 2019; v1 submitted 25 June, 2018;
originally announced June 2018.
-
The modified gravity lightcone simulation project I: Statistics of matter and halo distributions
Authors:
Christian Arnold,
Pablo Fosalba,
Volker Springel,
Ewald Puchwein,
Linda Blot
Abstract:
We introduce a set of four very high resolution cosmological simulations for exploring $f(R)$-gravity, with $2048^3$ particles in $768\,h^{-1}\textrm{Mpc}$ and $1536\,h^{-1}\textrm{Mpc}$ simulation boxes, both for a $|\overline{f_{R0}}| = 10^{-5}$ model and a $Λ$CDM comparison universe, making the set the largest simulations of $f(R)$-gravity to date. In order to mimic real observations, the simul…
▽ More
We introduce a set of four very high resolution cosmological simulations for exploring $f(R)$-gravity, with $2048^3$ particles in $768\,h^{-1}\textrm{Mpc}$ and $1536\,h^{-1}\textrm{Mpc}$ simulation boxes, both for a $|\overline{f_{R0}}| = 10^{-5}$ model and a $Λ$CDM comparison universe, making the set the largest simulations of $f(R)$-gravity to date. In order to mimic real observations, the simulations include a continuous 2D and 3D lightcone output which is dedicated to study lensing and clustering statistics in modified gravity. In this work, we present a detailed analysis and resolution study for the matter power spectrum in $f(R)$-gravity over a wide range of scales. We also analyse the angular matter power spectrum and lensing convergence on the lightcone. In addition, we investigate the impact of modified gravity on the halo mass function, matter and halo auto-correlation functions, linear halo bias and the concentration-mass relation. We find that the impact of $f(R)$-gravity is generally larger on smaller scales and increases with decreasing redshift. Comparing our simulations to state-of-the-art hydrodynamical simulations we confirm a degeneracy between $f(R)$-gravity and baryonic feedback in the matter power spectrum on small scales, but also find that scales around $k = 1\, h\, {\rm Mpc}^{-1}$ are promising to distinguish both effects. The lensing convergence power spectrum is increased in $f(R)$-gravity. Interestingly available numerical fits are in good agreement overall with our simulations for both standard and modified gravity, but tend to overestimate their relative difference on non-linear scales by few percent. We also find that the halo bias is lower in $f(R)$-gravity compared to general relativity, whereas halo concentrations are increased for unscreened halos.
△ Less
Submitted 24 May, 2018;
originally announced May 2018.
-
Assessment of the Information Content of the Power Spectrum and Bispectrum
Authors:
Kwan Chuen Chan,
Linda Blot
Abstract:
The covariance matrix of the matter and halo power spectrum and bispectrum are studied. Using a large suite of simulations, we find that the non-Gaussianity in the covariance is significant already at mildly nonlinear scales. We compute the leading disconnected non-Gaussian correction to the matter bispectrum covariance using perturbation theory, and find that the corrections result in good agreem…
▽ More
The covariance matrix of the matter and halo power spectrum and bispectrum are studied. Using a large suite of simulations, we find that the non-Gaussianity in the covariance is significant already at mildly nonlinear scales. We compute the leading disconnected non-Gaussian correction to the matter bispectrum covariance using perturbation theory, and find that the corrections result in good agreement in the mildly nonlinear regime. The shot noise contribution to the halo power spectrum and bispectrum covariance is computed using the Poisson model, and the model yields decent agreement with simulation results. However, when the shot noise is estimated from the individual realization, which is usually done in reality, we find that the halo covariance is substantially reduced and gets close to the Gaussian covariance. This is because most of the non-Gaussianity in the covariance arises from the fluctuations in the Poisson shot noise. We use the measured non-Gaussian covariance to access the information content of the power spectrum and bispectrum. The signal-to-noise ratio, S/N, of the matter and halo power spectrum levels off in the mildly nonlinear regime, $k \sim 0.1 - 0.2 \,\mathrm{Mpc} h^{-1}$. In the nonlinear regime the S/N of the matter and halo bispectrum increases but much slower than the Gaussian results suggest. We find that both the S/N for power spectrum and bispectrum are overestimated by the Gaussian covariances, but the problem being much more serious for the bispectrum. Because the bispectrum is affected strongly by nonlinearity and shot noise, inclusion of the bispectrum only adds modest amount of S/N compared to that of the power spectrum.
△ Less
Submitted 25 July, 2017; v1 submitted 20 October, 2016;
originally announced October 2016.
-
Non-Linear Matter Power Spectrum Covariance Matrix Errors and Cosmological Parameter Uncertainties
Authors:
Linda Blot,
Pier Stefano Corasaniti,
Luca Amendola,
Thomas D. Kitching
Abstract:
The covariance matrix of the matter power spectrum is a key element of the statistical analysis of galaxy clustering data. Independent realisations of observational measurements can be used to sample the covariance, nevertheless statistical sampling errors will propagate into the cosmological parameter inference potentially limiting the capabilities of the upcoming generation of galaxy surveys. Th…
▽ More
The covariance matrix of the matter power spectrum is a key element of the statistical analysis of galaxy clustering data. Independent realisations of observational measurements can be used to sample the covariance, nevertheless statistical sampling errors will propagate into the cosmological parameter inference potentially limiting the capabilities of the upcoming generation of galaxy surveys. The impact of these errors as function of the number of independent realisations has been previously evaluated for Gaussian distributed data. However, non-linearities in the late time clustering of matter cause departures from Gaussian statistics. Here, we address the impact of non-Gaussian errors on the sample covariance and precision matrix errors using a large ensemble of numerical N-body simulations. In the range of modes where finite volume effects are negligible ($0.1\lesssim k\,[h\,{\rm Mpc^{-1}}]\lesssim 1.2$) we find deviations of the estimated variance of the sample covariance with respect to Gaussian predictions above $\sim 10\%$ level. These reduce to about $\sim 5\%$ in the case of the precision matrix. Finally, we perform a Fisher analysis to estimate the effect of covariance errors on the cosmological parameter constraints. In particular, assuming Euclid-like survey characteristics we find that a number of independent realisation larger than $\gtrsim 5000$ is necessary to reduce the contribution of sample covariance errors to the cosmological parameter uncertainties at sub-percent level. We also show that restricting the analysis to large scales $k\lesssim0.2\,h\,{\rm Mpc^{-1}}$ results in a considerable loss in constraining power, while using the linear covariance to include smaller scales leads to an underestimation of the errors on the cosmological parameters.
△ Less
Submitted 13 April, 2016; v1 submitted 16 December, 2015;
originally announced December 2015.
-
Matter Power Spectrum Covariance Matrix from the DEUS-PUR ΛCDM simulations: Mass Resolution and non-Gaussian Errors
Authors:
Linda Blot,
Pier Stefano Corasaniti,
Jean-Michel Alimi,
Vincent Reverdy,
Yann Rasera
Abstract:
The upcoming generation of galaxy surveys will probe the distribution of matter in the universe with unprecedented accuracy. Measurements of the matter power spectrum at different scales and redshifts will provide stringent constraints on the cosmological parameters. However, on non-linear scales this will require an accurate evaluation of the covariance matrix. Here, we compute the covariance mat…
▽ More
The upcoming generation of galaxy surveys will probe the distribution of matter in the universe with unprecedented accuracy. Measurements of the matter power spectrum at different scales and redshifts will provide stringent constraints on the cosmological parameters. However, on non-linear scales this will require an accurate evaluation of the covariance matrix. Here, we compute the covariance matrix of the 3D matter density power spectrum for the concordance $Λ$CDM cosmology from an ensemble of N-body simulations of the Dark Energy Universe Simulation - Parallel Universe Runs (DEUS-PUR). This consists of 12288 realisations of a $(656\,h^{-1}\,\textrm{Mpc})^3$ simulation box with $256^3$ particles. We combine this set with an auxiliary sample of 96 simulations of the same volume with $1024^3$ particles. We find N-body mass resolution effect to be an important source of systematic errors on the covariance at high redshift and small intermediate scales. We correct for this effect by introducing an empirical statistical method which provide an accurate determination of the covariance matrix over a wide range of scales including the Baryon Oscillations interval. Contrary to previous studies that used smaller N-body ensembles, we find the power spectrum distribution to significantly deviate from expectations of a Gaussian random density field at $k\gtrsim 0.25\,h\,\textrm{Mpc}^{-1}$ and $z<0.5$. This suggests that in the case of finite volume surveys an unbiased estimate of the ensemble averaged band power at these scales and redshifts may require a careful assessment of non-Gaussian errors more than previously considered.
△ Less
Submitted 27 October, 2014; v1 submitted 10 June, 2014;
originally announced June 2014.