Showing papers by "Nabila Aghanim published in 2023"

Euclid preparation. XXIX. Forecasts for ten different higher-order weak lensing statistics

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 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 ($\Omega_{\rm m}$, $\sigma_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.

Can decaying dark matter models be a solution to the S 8 tension from the thermal Sunyaev-Zel’dovich effect?

Abstract: ,

Euclid preparation. XXIX. Water ice in spacecraft. I. The physics of ice formation and contamination

Abstract: Molecular contamination is a well-known problem in space flight. Water is the most common contaminant and alters numerous properties of a cryogenic optical system. Too much ice means that Euclid's calibration requirements and science goals cannot be met. Euclid must then be thermally decontaminated, a long and risky process. We need to understand how iced optics affect the data and when a decontamination is required. This is essential to build adequate calibration and survey plans, yet a comprehensive analysis in the context of an astrophysical space survey has not been done before. In this paper we look at other spacecraft with well-documented outgassing records, and we review the formation of thin ice films. A mix of amorphous and crystalline ices is expected for Euclid. Their surface topography depends on the competing energetic needs of the substrate-water and the water-water interfaces, and is hard to predict with current theories. We illustrate that with scanning-tunnelling and atomic-force microscope images. Industrial tools exist to estimate contamination, and we must understand their uncertainties. We find considerable knowledge errors on the diffusion and sublimation coefficients, limiting the accuracy of these tools. We developed a water transport model to compute contamination rates in Euclid, and find general agreement with industry estimates. Tests of the Euclid flight hardware in space simulators did not pick up contamination signals; our in-flight calibrations observations will be much more sensitive. We must understand the link between the amount of ice on the optics and its effect on Euclid's data. Little research is available about this link, possibly because other spacecraft can decontaminate easily, quenching the need for a deeper understanding. In our second paper we quantify the various effects of iced optics on spectrophotometric data.

Systematic error mitigation for the PIXIE Fourier transform spectrometer

Abstract: The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission concept to measure the spectrum and polarization of the cosmic microwave background. Cosmological signals are small compared to the instantaneous instrument noise, requiring strict control of instrumental signals. The instrument design provides multiple levels of null operation, signal modulation, and signal differences, with only few-percent systematic error suppression required at each level. Jackknife tests based on discrete instrument symmetries provide an independent means to identify, model, and remove remaining instrumental signals. We use detailed time-ordered simulations, including realistic performance and tolerance parameters, to evaluate the instrument response to broad classes of systematic errors for both spectral distortions and polarization. The largest systematic errors contribute additional white noise at the few-percent level compared to the dominant photon noise. Coherent instrumental effects which do not integrate down are smaller still, and remain several orders of magnitude below the targeted cosmological signals.

Euclid preparation. XXVII. A UV-NIR spectral atlas of compact planetary nebulae for wavelength calibration

Abstract: The Euclid mission will conduct an extragalactic survey over 15000 deg$^2$ of the extragalactic sky. The spectroscopic channel of the Near-Infrared Spectrometer and Photometer (NISP) has a resolution of $R\sim450$ for its blue and red grisms that collectively cover the $0.93$--$1.89 $\micron;range. NISP will obtain spectroscopic redshifts for $3\times10^7$ galaxies for the experiments on galaxy clustering, baryonic acoustic oscillations, and redshift space distortion. The wavelength calibration must be accurate within $5$\AA to avoid systematics in the redshifts and downstream cosmological parameters. The NISP pre-flight dispersion laws for the grisms were obtained on the ground using a Fabry-Perot etalon. Launch vibrations, zero gravity conditions, and thermal stabilisation may alter these dispersion laws, requiring an in-flight recalibration. To this end, we use the emission lines in the spectra of compact planetary nebulae (PNe), which were selected from a PN data base. To ensure completeness of the PN sample, we developed a novel technique to identify compact and strong line emitters in Gaia spectroscopic data using the Gaia spectra shape coefficients. We obtained VLT/X-SHOOTER spectra from $0.3$ to $2.5$ \micron;for 19 PNe in excellent seeing conditions and a wide slit, mimicking Euclid's slitless spectroscopy mode but with 10 times higher spectral resolution. Additional observations of one northern PN were obtained in the $0.80$--$1.90$ \micron range with the GMOS and GNIRS instruments at the Gemini North observatory. The collected spectra were combined into an atlas of heliocentric vacuum wavelengths with a joint statistical and systematic accuracy of 0.1 \AA in the optical and 0.3 \AA in the near-infrared. The wavelength atlas and the related 1D and 2D spectra are made publicly available.

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

Abstract: The future Euclid space satellite mission will offer an invaluable opportunity to constrain modifications to general relativity at cosmic scales. We focus on modified gravity models characterised, at linear scales, by a scale-independent growth of perturbations while featuring different testable types of derivative screening mechanisms at smaller nonlinear scales. We consider 3 specific models, namely Jordan-Brans-Dicke (JBD), the normal branch of Dvali-Gabadadze-Porrati (nDGP) gravity and $k$-mouflage (KM) gravity. We provide forecasts from spectroscopic and photometric primary probes by Euclid on the cosmological parameters and the extra parameters of the models, respectively, $\omega_{\rm BD}$, $\Omega_{\rm rc}$ and $\epsilon_{2,0}$, which quantify the deviations from general relativity. This analysis will improve our knowledge of the cosmology of these modified gravity models. The forecasts analysis employs the Fisher matrix method applied to weak lensing (WL); photometric galaxy clustering (GC$_{ph}$); spectroscopic galaxy clustering (GC$_{sp}$) and the cross-correlation (XC) between GC$_{ph}$ and WL. For the Euclid survey specifications we define three scenarios, characterised by different cuts in $\ell$ and $k$, to assess the constraining power of nonlinear scales. For each model we consider two fiducial values for the corresponding model parameter. In an optimistic setting at 68.3\% confidence interval, with Euclid alone we find the following percentage relative errors: for $\log_{10}{\omega_{\rm BD}}$, with a fiducial value of $\omega_{\rm BD}=800$, 35% using GC$_{sp}$ alone, 3.6% using GC$_{ph}$+WL+XC and 3.3% using GC$_{ph}$+WL+XC+GC$_{sp}$; for $\log_{10}{\Omega_{\rm rc}}$, with a fiducial value of $\Omega_{\rm rc}=0.25$, we find respectively 90%, 20% and 17%; finally, for $\epsilon_{2,0}=-0.04$ respectively 5%, 0.15% and 0.14%. (abridged)

Simulating the LOcal Web (SLOW): I. Anomalies in the local density field

Abstract: Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Methods (abridged): We used a 500 Mpc/h large constrained simulation of the LU with initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalogue and follow galaxy formation physics directly in the hydro-dynamical simulations to base the comparison on stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 Mpc/h large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that haloes and galaxies in our constrained simulation trace the local DM density field very differently. Thereby, this simulation reproduces the observed 50% under-density of galaxy clusters and groups within the sphere of ~100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ~1.5$\sigma$ feature. At the same time, the simulation reproduces the observed over-density of massive galaxy clusters within the same sphere, which on its own also corresponds to a ~1.5$\sigma$ feature. Interestingly, we find that only 44 out of 15635 random realizations (i.e. 0.28%) are matching both anomalies, making the LU to be a ~3$\sigma$ environment. We finally compared a mock galaxy catalogue with the observed distribution of galaxies in the LU, finding also a match to the observed factor of two over-density at ~16 Mpc as well as the observed 15% under-density at ~40 Mpc distance.

A VLT/VIMOS view of two Planck multiple-cluster systems: Structure and galaxy properties

Abstract: We analysed spectroscopic data obtained with VLT-VIMOS for two multiple-cluster systems, PLCKG$214.6+36.9$ and PLCKG$334.8-38.0$, discovered via their thermal Sunyaev-Zel'dovich signal by $Planck$. Combining the Optical spectroscopy, for the redshift determination, and photometric data from galaxy surveys (SDSS, WISE, DESI), we were able to study the structure of the two multiple-cluster systems, to determine their nature and the properties of their member galaxies. We found that the two systems are populated mainly with passive galaxies and that PLCKG$214.6+36.9$ consists of a pair of clusters at redshift $z = 0.445$ and a background isolated cluster at $z = 0.498$, whereas the system PLCKG$334.8-38.0$ is a chance association of three independent clusters at redshifts $z = 0.367$, $z =0.292$, and $z = 0.33$. We also find evidence for remaining star formation activity in the highest-redshift cluster of PLCKG$214.6+36.9$, at $z = 0.498$.

Euclid: Validation of the MontePython forecasting tools

Abstract: The Euclid mission of the European Space Agency will perform a survey of weak lensing cosmic shear and galaxy clustering in order to constrain cosmological models and fundamental physics. We expand and adjust the mock Euclid likelihoods of the MontePython software in order to match the exact recipes used in previous Euclid Fisher matrix forecasts for several probes: weak lensing cosmic shear, photometric galaxy clustering, the cross-correlation between the latter observables, and spectroscopic galaxy clustering. We also establish which precision settings are required when running the Einstein-Boltzmann solvers CLASS and CAMB in the context of Euclid. For the minimal cosmological model, extended to include dynamical dark energy, we perform Fisher matrix forecasts based directly on a numerical evaluation of second derivatives of the likelihood with respect to model parameters. We compare our results with those of other forecasting methods and tools. We show that such MontePython forecasts agree very well with previous Fisher forecasts published by the Euclid Collaboration, and also, with new forecasts produced by the CosmicFish code, now interfaced directly with the two Einstein-Boltzmann solvers CAMB and CLASS. Moreover, to establish the validity of the Gaussian approximation, we show that the Fisher matrix marginal error contours coincide with the credible regions obtained when running Monte Carlo Markov Chains with MontePython while using the exact same mock likelihoods. The new Euclid forecast pipelines presented here are ready for use with additional cosmological parameters, in order to explore extended cosmological models.

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

Abstract: ,

Testing decaying dark matter models as a solution to the S8 tension with the thermal Sunyaev-Zel'dovich effect

Abstract: Considering possible solutions to the $S_8$ tension between the Planck cosmic microwave background (CMB) measurement and low-redshift probes, we extended the standard $\Lambda$CDM cosmological model by including decay of dark matter (DDM). We first tested the DDM model in which dark matter decays into a form of noninteracting dark radiation. Under this DDM model, we investigated the impacts of DDM on the Sunyaev Zel'dovich (SZ) effect by varying the decay lifetime, $\Gamma^{-1}$, including the background evolution in cosmology and the nonlinear prescription in the halo mass function. We performed a cosmological analysis under the assumption of this extended cosmological model by combining the latest high-redshift Planck CMB measurement and low-redshift measurements of the SZ power spectrum as well as the baryonic acoustic oscillations (BAO) and luminosity distances to type Ia supernovae (SNIa). Our result shows a preference for $\Gamma^{-1} \sim 220$ Gyr with a lower bound on the decay lifetime of $\sim$ 38 Gyr at 95\% confidence level. Additionally, we tested the other DDM model in which dark matter decays into warm dark matter and dark radiation. This model supports $\Gamma^{-1} \sim 137$ Gyr to resolve the $S_8$ tension with a lower bound on the decay lifetime of $\sim$ 24 Gyr at 95\% confidence level. Comparing these two models, we find that the second leads to slightly better reconciliation of the $S_8$ tension.

The cosmic web around the Coma cluster from constrained cosmological simulations

Abstract: Galaxy clusters in the Universe occupy the important position of nodes of the cosmic web. They are connected among them by filaments, elongated structures composed of dark matter, galaxies, and gas. The connection of galaxy clusters to filaments is important, as it is related to the process of matter accretion onto the former. For this reason, investigating the connections to the cosmic web of massive clusters, especially well-known ones for which a lot of information is available, is a hot topic in astrophysics. In a previous work, we performed an analysis of the filament connections of the Coma cluster of galaxies, as detected from the observed galaxy distribution. In this work we resort to a numerical simulation whose initial conditions are constrained to reproduce the local Universe, including the region of the Coma cluster to interpret our observations in an evolutionary context. We detect the filaments connected to the simulated Coma cluster and perform an accurate comparison with the cosmic web configuration we detect in observations. We perform an analysis of the halos’ spatial and velocity distributions close to the filaments in the cluster outskirts. We conclude that, although not significantly larger than the average, the flux of accreting matter on the simulated Coma cluster is significantly more collimated close to the filaments with respect to the general isotropic accretion flux. This paper is the first example of such a result and the first installment in a series of publications which will explore the build-up of the Coma cluster system in connection to the filaments of the cosmic web as a function of redshift.

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

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 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\sigma$) cosmological biases in one or more parameters, and accordingly must be accounted for. In total, over all effects, we find biases in $\Omega_{\rm m}$, $\Omega_{\rm b}$, $h$, and $\sigma_{8}$ of $0.73\sigma$, $0.28\sigma$, $0.25\sigma$, and $-0.79\sigma$, respectively, for flat $\Lambda$CDM. For the $w_0w_a$CDM case, we find biases in $\Omega_{\rm m}$, $\Omega_{\rm b}$, $h$, $n_{\rm s}$, $\sigma_{8}$, and $w_a$ of $1.49\sigma$, $0.35\sigma$, $-1.36\sigma$, $1.31\sigma$, $-0.84\sigma$, and $-0.35\sigma$, respectively; which are increased relative to the $\Lambda$CDM due to additional degeneracies as a function of redshift and scale.

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

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

Euclid preparation. XXX. Performance assessment of the NISP red grism through spectroscopic simulations for the wide and deep surveys

Abstract: This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$\sigma$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

Velocity waves in the Hubble diagram: signature of local galaxy clusters

Abstract: The Universe expansion rate is modulated around local inhomogeneities due to their gravitational potential. Velocity waves are then observed around galaxy clusters in the Hubble diagram. This paper studies them in a ∼ 738 Mpc wide, with 2048 3 particles, cosmological simulation of our cosmic environment (a.k.a. CLONE: Constrained LOcal & Nesting Environment Simulation). For the first time, the simulation shows that velocity waves that arise in the lines-of-sight of the most massive dark matter halos agree with those observed in local galaxy velocity catalogs in the lines-of-sight of Coma and several other local (Abell) clusters. For the best-constrained clusters such as Virgo and Centaurus, i.e