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.

In this work, we investigate cosmologies where the gravitational constant varies in time, with the aim of explaining the accelerated expansion without a cosmological constant. We achieve this by considering a phenomenological extension to general relativity, modifying Einstein's field equations such that $G$ is a function of time, $G(t)$, and we preserve the geometrical consistency (Bianchi identity) together with the usual conservation of energy by introducing a new tensor field to the equations. In order to have concrete expressions to compare with cosmological data, we posit additional properties to this tensor field, in a way that it can be interpreted as a response of spacetime to a variation of $G$. Namely, we require that the energy this tensor represents is nonzero only when there is a time variation of $G$, and its energy depends on the scale factor only because of its coupling to $G$ and the matter and radiation energy densities. Focusing on the accelerated expansion period, we use type Ia supernovae and baryon acoustic oscillation data to determine the best fit of the cosmological parameters as well as the required variation in the gravitational constant. As a result, we find that it is possible to explain the accelerated expansion of the Universe with a variation of $G$ and no cosmological constant. The obtained variation of $G$ stays under 10 \% of its current value in the investigated redshift range and it is consistent with the local observations of $\dot{G}/G$.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

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.

With the recent increase in precision of our cosmological datasets, measurements of $\Lambda$CDM model parameter provided by high- and low-redshift observations started to be in tension, i.e., the obtained values of such parameters were shown to be significantly different in a statistical sense. In~this work we tackle the tension on the value of the Hubble parameter, $H_0$, and the weighted amplitude of matter fluctuations, $S_8$, obtained from local or low-redshift measurements and from cosmic microwave background (CMB) observations. We combine the main approaches previously used in the literature by extending the cosmological model and accounting for extra systematic uncertainties. With such analysis we aim at exploring non standard cosmological models, implying deviation from a cosmological constant driven acceleration of the Universe expansion, in the presence of additional uncertainties in measurements. In more detail, we reconstruct the Dark Energy equation of state as a function of redshift, while we study the impact of type-Ia supernovae (SNIa) redshift-dependent astrophysical systematic effects on these tensions. We consider a SNIa intrinsic luminosity dependence on redshift due to the star formation rate in its environment, or the metallicity of the progenitor. We find that the $H_0$ and $S_8$ tensions can be significantly alleviated, or~even removed, if we account for varying Dark Energy for SNIa and CMB data. However, the tensions remain when we add baryon acoustic oscillations (BAO) data into the analysis, even after the addition of extra SNIa systematic uncertainties. This points towards the need of either new physics beyond late-time Dark Energy, or other unaccounted systematic effects (particulary in BAO measurements), to fully solve the present tensions.

In this work we investigate cosmologies with a time-varying gravitational constant, $G(t)$, as an alternative to having a cosmological constant. Since, it is known in the literature that, in general relativity, having a variable gravitational constant breaks the conservation of energy unless further dynamical components are added, we consider a Newtonian approach, in which the energy conservation is ensured. We investigate whether these cosmologies are able to reproduce the cosmological acceleration without a cosmological constant, nor any other sort of dark energy fluid. After constructing the Friedmann-Lema\^itre equation for a late-stage, matter dominated universe by starting with Newton's second law, where $G$ is taken as a function of time, we create specific models to test against cosmological observations. In this Friedmann-Lema\^itre equation we obtain a second gravitational constant, $G^*$, related to the original $G$ in the acceleration equation. For the tests we focus on the acceleration period and use low-redshift probes: type Ia supernovae (SNIa), baryon acoustic oscillations, and cosmic chronometers, taking also into account a possible change in the supernova intrinsic luminosity with redshift. As a result, we obtain several models without a cosmological constant with similar or even slightly better $\chi^2$ values than the standard $\Lambda$CDM cosmology. When SNIa luminosity dependence on redshift is added, a model with a $G$ exponentially decreasing to zero today can explain the observations. In the cases where SNIa are assumed to be standard candles, the observations favour a negative $G$ at large scales to produce cosmic acceleration. We conclude that these models offer a viable interpretation of the late-time universe observations.

The concordance model in cosmology, $\Lambda$CDM, is able to fit the main cosmological observations with a high level of accuracy. However, around 95% of the energy content of the Universe within this framework remains still unknown. In this work we focus on the dark matter component and we investigate the generalized dark matter (GDM) model, which allows for non-pressure-less dark matter and a non-vanishing sound speed and viscosity. We first focus on current observations, showing that GDM could alleviate the tension between cosmic microwave background and weak lensing observations. We then investigate the ability of the photometric Euclid survey (photometric galaxy clustering, weak lensing, and their cross-correlations) to constrain the nature of dark matter. We conclude that Euclid will provide us with very good constraints on GDM, enabling us to better understand the nature of this fluid, but a non-linear recipe adapted to GDM is clearly needed in order to correct for non-linearities and get reliable results down to small scales.

The $\Lambda$CDM model is the current standard model in cosmology thanks to its ability to reproduce the observations. Its first observational evidence appeared from the type Ia supernovae (SNIa) Hubble diagram. However, there has been some debate in the literature concerning the statistical treatment of SNIa. In this paper we relax the standard assumption that SNIa intrinsic luminosity is independent of the redshift, and we examine whether it may have an impact on the accelerated nature of the expansion of the Universe. In order to be as general as possible, we reconstruct the expansion rate of the Universe through a cubic spline interpolation fitting observations of different probes: SNIa, baryon acoustic oscillations (BAO), and the high-redshift information from the cosmic microwave background (CMB). We show that when SNIa intrinsic luminosity is not allowed to vary as a function of the redshift, cosmic acceleration is definitely proven in a model-independent approach. However, allowing for a redshift dependence, a non-accelerated reconstruction of the expansion rate is able to fit, as well as $\Lambda$CDM, the combination of SNIa and BAO data, both treating the BAO standard ruler $r_d$ as a free parameter, or adding the recently published prior from CMB observations. We further extend the analysis by including the CMB data, and we show that a non-accelerated reconstruction is able to nicely fit this combination of low and high-redshift data. In this work we present a model-independent reconstruction of a non-accelerated expansion rate of the Universe that is able to nicely fit all the main background cosmological probes. However, the predicted value of $H_0$ is in tension with recent direct measurements. Our analysis points out that a final, reliable, and consensual value for $H_0$ would be critical to definitively prove the cosmic acceleration in a model-independent way. [Abridged]

Context: The cosmological concordance model ($\Lambda$CDM) matches the cosmological observations exceedingly well. This model has become the standard cosmological model with the evidence for an accelerated expansion provided by the type Ia supernovae (SNIa) Hubble diagram. However, the robustness of this evidence has been addressed recently with somewhat diverging conclusions. Aims: The purpose of this paper is to assess the robustness of the conclusion that the Universe is indeed accelerating if we rely only on low-redshift (z$\lesssim$2) observations, that is to say with SNIa, baryonic acoustic oscillations, measurements of the Hubble parameter at different redshifts, and measurements of the growth of matter perturbations. Methods: We used the standard statistical procedure of minimizing the $\chi^2$ function for the different probes to quantify the goodness of fit of a model for both $\Lambda$CDM and a simple nonaccelerated low-redshift power law model. In this analysis, we do not assume that supernovae intrinsic luminosity is independent of the redshift, which has been a fundamental assumption in most previous studies that cannot be tested. Results: We have found that, when SNIa intrinsic luminosity is not assumed to be redshift independent, a nonaccelerated low-redshift power law model is able to fit the low-redshift background data as well as, or even slightly better, than $\Lambda$CDM. When measurements of the growth of structures are added, a nonaccelerated low-redshift power law model still provides an excellent fit to the data for all the luminosity evolution models considered. Conclusions: Without the standard assumption that supernovae intrinsic luminosity is independent of the redshift, low-redshift probes are consistent with a nonaccelerated universe.

In this paper we study the consequences of relaxing the hypothesis of the pressureless nature of the dark matter component when determining constraints on dark energy. To this aim we consider simple generalized dark matter models with constant equation of state parameter. We find that present-day low-redshift probes (type-Ia supernovae and baryonic acoustic oscillations) lead to a complete degeneracy between the dark energy and the dark matter sectors. However, adding the cosmic microwave background (CMB) high-redshift probe restores constraints similar to those on the standard $\Lambda$CDM model. We then examine the anticipated constraints from the galaxy clustering probe of the future Euclid survey on the same class of models, using a Fisher forecast estimation. We show that the Euclid survey allows us to break the degeneracy between the dark sectors, although the constraints on dark energy are much weaker than with standard dark matter. The use of CMB in combination allows us to restore the high precision on the dark energy sector constraints.

In this work we study the consequences of allowing non pressureless dark matter when determining dark energy constraints. We show that present-day dark energy constraints from low-redshift probes are extremely degraded when allowing this dark matter variation. However, adding the cosmic microwave background (CMB) we can recover the $w_{DM} = 0$ case constraints. We also show that with the future Euclid redshift survey we expect to largely improve these constraints; but, without the complementary information of the CMB, there is still a strong degeneracy between dark energy and dark matter equation of state parameters.

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.