In recent years forecasting activities have become a very important tool for designing and optimising large scale structure surveys. To predict the performance of such surveys, the Fisher matrix formalism is frequently used as a fast and easy way to compute constraints on cosmological parameters. Among them lies the study of the properties of dark energy which is one of the main goals in modern cosmology. As so, a metric for the power of a survey to constrain dark energy is provided by the Figure of merit (FoM). This is defined as the inverse of the surface contour given by the joint variance of the dark energy equation of state parameters $\{w_0,w_a\}$ in the Chevallier-Polarski-Linder parameterisation, which can be evaluated from the covariance matrix of the parameters. This covariance matrix is obtained as the inverse of the Fisher matrix. Inversion of an ill-conditioned matrix can result in large errors on the covariance coefficients if the elements of the Fisher matrix have been estimated with insufficient precision. The conditioning number is a metric providing a mathematical lower limit to the required precision for a reliable inversion, but it is often too stringent in practice for Fisher matrices with size larger than $2\times2$. In this paper we propose a general numerical method to guarantee a certain precision on the inferred constraints, like the FoM. It consists on randomly vibrating (perturbing) the Fisher matrix elements with Gaussian perturbations of a given amplitude, and then evaluating the maximum amplitude that keeps the FoM within the chosen precision. The steps used in the numerical derivatives and integrals involved in the calculation of the Fisher matrix elements can then be chosen accordingly in order to keep the precision of the Fisher matrix elements below this maximum amplitude...

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The influence of considering a generalized dark matter (GDM) model, which allows for a non-pressure-less dark matter and a non-vanishing sound speed in the non-linear spherical collapse model is discussed for the Einstein-de Sitter-like (EdSGDM) and $\Lambda$GDM models. By assuming that the vacuum component responsible for the accelerated expansion of the Universe is not clustering and therefore behaving similarly to the cosmological constant $\Lambda$, we show how the change in the GDM characteristic parameters affects the linear density threshold for collapse of the non-relativistic component ($\delta_{\rm c}$) and its virial overdensity ($\Delta_{\rm V}$). We found that the generalized dark matter equation of state parameter $w_{\rm gdm}$ is responsible for lower values of the linear overdensity parameter as compared to the standard spherical collapse model and that this effect is much stronger than the one induced by a change in the generalized dark matter sound speed $c^2_{\rm s, gdm}$. We also found that the virial overdensity is only slightly affected and mostly sensitive to the generalized dark matter equation of state parameter $w_{\rm gdm}$. These effects could be relatively enhanced for lower values of the matter density. Finally, we found that the effects of the additional physics on $\delta_{\rm c}$ and $\Delta_{\rm V}$, when translated to non-linear observables such as the halo mass function, induce an overall deviation of about 40\% with respect to the standard $\Lambda$CDM model at late times for high mass objects. However, within the current linear constraints for $c^2_{\rm s, gdm}$ and $w_{\rm gdm}$, we found that these changes are the consequence of properly taking into account the correct linear matter power spectrum for the GDM model while the effects coming from modifications in the spherical collapse model remain negligible.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

The $\Lambda$CDM concordance model is very successful at describing our Universe with high accuracy and few parameters. Despite its successes, a few tensions persist; most notably, the best-fit $\Lambda$CDM model, as derived from the Planck CMB data, largely overpredicts the abundance of SZ clusters when using their standard mass calibration. Whether this is a sign of an incorrect calibration or the need for new physics remains a matter of debate. Here we examined two simple extensions of the standard model and their ability to release this tension: massive neutrinos and a simple modified gravity model via a non-standard growth index $\gamma$. We used both the Planck CMB and SZ cluster counts as datasets, with or without local X-ray clusters. In the case of massive neutrinos, the SZ calibration $(1-b)$ is constrained to $0.59^{+0.03}_{-0.04}$ (68\%), more than 5$\sigma$ away from its standard value $\sim0.8$. We found little correlation between $\sum m_\nu$ and $(1-b)$, corroborating previous conclusions derived from X-ray clusters; massive neutrinos do not alleviate the cluster-CMB tension. With our simple $\gamma$ model, we found a large correlation between calibration and growth index but contrary to local X-ray clusters, SZ clusters are able to break the degeneracy between the two thanks to their extended $z$ range. The calibration $(1-b)$ was then constrained to $0.60^{+0.05}_{-0.07}$, leading to an interesting constraint on $\gamma=0.60\pm 0.13$. When both massive neutrinos and modified gravity were allowed, preferred values remained centred on standard $\Lambda$CDM values, but $(1-b)\sim0.8$ was allowed (though only at the $2\sigma$ level) provided $\sum m_\nu\sim0.34 $ eV and $\gamma\sim0.8$. We conclude that massive neutrinos do not relieve the cluster-CMB tension and that a calibration close to the standard value $0.8$ would call for new physics in the gravitational sector.

In recent years, the amplitude of matter fluctuations inferred from low-redshift probes has been found to be generally lower than the value derived from CMB observations in the $\Lambda$CDM model. This tension has been exemplified by Sunyaev-Zel'dovich and X-ray cluster counts which, when using their Planck standard cluster mass calibration, yield a value of $\sigma_8$ , appreciably lower than estimations based on the latest Planck CMB measurements. In this work we examine whether non-minimal neutrino masses can alleviate this tension substantially. We used the cluster X-ray temperature distribution function derived from a flux-limited sample of local X-ray clusters, combined with Planck CMB measurements. These datasets were compared to $\Lambda$CDM predictions based on recent mass function, adapted to account for the effects of massive neutrinos. Treating the clusters mass calibration as a free parameter, we examined whether the data favours neutrino masses appreciably higher than the minimal 0.06 eV value. Using Markov chain Monte Carlo methods, we found no significant correlation between the mass calibration of clusters and the sum of neutrino masses, meaning that massive neutrinos do not noticeably alleviate the above-mentioned Planck CMB--clusters tension. The addition of other datasets (BAO and Ly-$\alpha$) reinforces those conclusions. As an alternative possible solution to the tension, we introduced a simple, phenomenological modification of gravity by letting the growth index $\gamma$ vary as an additional free parameter. We find that the cluster mass calibration is robustly correlated with the $\gamma$ parameter, insensitively to the presence of massive neutrinos or/and additional data used. We conclude that the standard Planck mass calibration of clusters, if consolidated, would represent evidence for new physics beyond $\Lambda$CDM with massive neutrinos.

Evolution and abundance of the large-scale structures we observe today, such as clusters of galaxies, is sensitive to the statistical properties of dark matter primordial density fluctuations, which is assumed to follow a Gaussian probability distribution function. Within this assumption, a significant disagreement have been found between clusters counts made by Planck and their prediction when calibrated by CMB angular power spectrum. The purpose of this work is to relax the Gaussianty assumption and test if Non-Gaussianity in dark matter primordial density fluctuations, could alleviate the tension.

The abundance of clusters is a classical cosmological probe sensitive to both the geometrical aspects and the growth rate of structures. The abundance of clusters of galaxies measured by Planck has been found to be in tension with the prediction of the LCDM models normalized to Planck CMB fluctuations power spectra. The same tension appears with X-ray cluster local abundance. Massive neutrinos and modified gravity are two possible solutions to fix this tension. Alternatively, others options include a bias in the selection procedure or in the mass calibration of clusters. We present a study, based on our recent work, updating the present situation on this topic and discuss the likelihood of the various options.

In this work we examine what are the cosmological implications of allowing the geometrical curvature density to behave independently from the energy density contents. Using the full data extracted by Planck mission from CMB, combined with BAO and SNIa measurements, we derive, in the light of this approach, new constraints on the cosmological parameters. In particular we determine the behavior of the curvature dark energy degeneracy when allowing a varying equation of state for the latter. We also examine whether this approach could bridge the gap recently found between the Hubble parameter value determined from CMB and that from the local universe measurements

Despite the ability of the cosmological concordance model ($\Lambda$CDM) to describe the cosmological observations exceedingly well, power law expansion of the Universe scale radius, $R(t)\propto t^n$, has been proposed as an alternative framework. We examine here these models, analyzing their ability to fit cosmological data using robust model comparison criteria. Type Ia supernovae (SNIa), baryonic acoustic oscillations (BAO) and acoustic scale information from the cosmic microwave background (CMB) have been used. We find that SNIa data either alone or combined with BAO can be well reproduced by both $\Lambda$CDM and power law expansion models with $n\sim 1.5$, while the constant expansion rate model $(n=1)$ is clearly disfavored. Allowing for some redshift evolution in the SNIa luminosity essentially removes any clear preference for a specific model. The CMB data are well known to provide the most stringent constraints on standard cosmological models, in particular, through the position of the first peak of the temperature angular power spectrum, corresponding to the sound horizon at recombination, a scale physically related to the BAO scale. Models with $n\geq 1$ lead to a divergence of the sound horizon and do not naturally provide the relevant scales for the BAO and the CMB. We retain an empirical footing to overcome this issue: we let the data choose the preferred values for these scales, while we recompute the ionization history in power law models, to obtain the distance to the CMB. In doing so, we find that the scale coming from the BAO data is not consistent with the observed position of the first peak of the CMB temperature angular power spectrum for any power law cosmology. Therefore, we conclude that when the three standard probes are combined, the $\Lambda$CDM model is very strongly favored over any of these alternative models, which are then essentially ruled out.

Halo Occupation Distribution (HOD) is a model giving the average number of galaxies in a dark matter halo, function of its mass and other intrinsic properties, like distance from halo center, luminosity and redshift of its constituting galaxies. It is believed that these parameters could also be related to the galaxy history of formation. We want to investigate more this relation in order to test and better refine this model. To do that, we extract HOD indicators from EUCLID mock catalogs for different luminosity cuts and for redshifts ranges going from 0.1 < z < 3.0. We study and interpret the trends of indicators function of these variations and tried to retrace galaxy formation history following the idea that galaxy evolution is the combination rather than the conflict of the two main proposed ideas nowadays: the older hierarchical mass merger driven paradigm and the recent downsizing star formation driven approach.