The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

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.

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian streaming model to calculate the quasi non-linear redshift space matter power spectrum multipoles from a non-linear prescription for the matter power spectrum, which remains agnostic about gravity or dark energy. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model. When comparing the redshift space multipoles for halos, we find that a linear bias implementation is insufficient to be competitive with the TNS-based model used in the BOSS clustering analysis, which however has additional degrees of freedom. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.

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.

This paper is a pedagogical introduction to models of gravity and how to constrain them through cosmological observations. We focus on the Horndeski scalar-tensor theory and on the quantities that can be measured with a minimum of assumptions. Alternatives or extensions of General Relativity have been proposed ever since its early years. Because of Lovelock theorem, modifying gravity in four dimensions typically means adding new degrees of freedom. The simplest way is to include a scalar field coupled to the curvature tensor terms. The most general way of doing so without incurring in the Ostrogradski instability is the Horndeski Lagrangian and its extensions. Testing gravity means therefore, in its simplest term, testing the Horndeski Lagrangian. Since local gravity experiments can always be evaded by assuming some screening mechanism or that baryons are decoupled, or even that the effects of modified gravity are visible only at early times, we need to test gravity with cosmological observations in the late universe (large-scale structure) and in the early universe (cosmic microwave background). In this work we review the basic tools to test gravity at cosmological scales, focusing on model-independent measurements.

We discuss the cosmological phenomenology of biscalar-tensor models displaying a maximally symmetric Einstein-frame kinetic sector and constructed on the basis of scale symmetry and volume-preserving diffeomorphisms. These theories contain a single dimensionful parameter $\Lambda_0$-associated with the invariance under the aforementioned restricted coordinate transformations-and a massless dilaton field. At large field values these scenarios lead to inflation with no generation of isocurvature perturbations. The corresponding predictions depend only on two dimensionless parameters, which characterize the curvature of the field--manifold and the leading order behavior of the inflationary potential. For $\Lambda_0=0$ the scale symmetry is unbroken and the dilaton admits only derivative couplings to matter, evading all fifth force constraints. For $\Lambda_0\neq 0$ the field acquires a run-away potential that can support a dark energy dominated era at late times. We confront a minimalistic realization of this appealing framework with present data sets. The impact of possible consistency relations among the early and late Universe dynamics that could appear within this setting is also discussed.

The gravitational slip parameter is an important discriminator between large classes of gravity theories at cosmological and astrophysical scales. In this work we use a combination of simulated information of galaxy cluster mass profiles, inferred by Strong+Weak lensing analyses and by the study of the dynamics of the cluster member galaxies, to reconstruct the gravitational slip parameter $\eta$ and predict the accuracy with which it can be constrained with current and future galaxy cluster surveys. Performing a full-likelihood statistical analysis, we show that galaxy cluster observations can constrain $\eta$ down to the percent level already with a few tens of clusters. We discuss the significance of possible systematics, and show that the cluster masses and numbers of galaxy members used to reconstruct the dynamics mass profile have a mild effect on the predicted constraints.

In the context of the effective field theory of dark energy (EFT) we perform agnostic explorations of Horndeski gravity. We choose two parametrizations for the free EFT functions, namely a power law and a dark energy density-like behaviour on a non trivial Chevallier-Polarski-Linder background. We restrict our analysis to those EFT functions which do not modify the speed of propagation of gravitational waves. Among those, we prove that one specific function cannot be constrained by data, since its contribution to the observables is below the cosmic variance, although we show it has a relevant role in defining the viable parameter space. We place constraints on the parameters of these models combining measurements from present day cosmological datasets and we prove that the next generation galaxy surveys can improve such constraints by one order of magnitude. We then verify the validity of the quasi-static limit within the sound horizon of the dark field, by looking at the phenomenological functions $\mu$ and $\Sigma$, associated respectively to clustering and lensing potentials. Furthermore, we notice up to $5\%$ deviations in $\mu, \Sigma$ with respect to General Relativity at scales smaller than the Compton one. For the chosen parametrizations and in the quasi-static limit, future constraints on $\mu$ and $\Sigma$ can reach the $1\%$ level and will allow us to discriminate between certain models at more than $3\sigma$, provided the present best-fit values remain.

In this work, we use recent data on the Hubble expansion rate $H(z)$, the quantity $f\sigma_8(z)$ from redshift space distortions and the statistic $E_g$ from clustering and lensing observables to constrain in a model-independent way the linear anisotropic stress parameter $\eta$. This estimate is free of assumptions about initial conditions, bias, the abundance of dark matter and the background expansion. We denote this observable estimator as $\eta_{{\rm obs}}$. If $\eta_{{\rm obs}}$ turns out to be different from unity, it would imply either a modification of gravity or a non-perfect fluid form of dark energy clustering at sub-horizon scales. Using three different methods to reconstruct the underlying model from data, we report the value of $\eta_{{\rm obs}}$ at three redshift values, $z=0.29, 0.58, 0.86$. Using the method of polynomial regression, we find $\eta_{{\rm obs}}=0.57\pm1.05$, $\eta_{{\rm obs}}=0.48\pm0.96$, and $\eta_{{\rm obs}}=-0.11\pm3.21$, respectively. Assuming a constant $\eta_{{\rm obs}}$ in this range, we find $\eta_{{\rm obs}}=0.49\pm0.69$. We consider this method as our fiducial result, for reasons clarified in the text. The other two methods give for a constant anisotropic stress $\eta_{{\rm obs}}=0.15\pm0.27$ (binning) and $\eta_{{\rm obs}}=0.53 \pm 0.19$ (Gaussian Process). We find that all three estimates are compatible with each other within their $1\sigma$ error bars. While the polynomial regression method is compatible with standard gravity, the other two methods are in tension with it.

The Higgs-Dilaton model is a scale-invariant extension of the Standard Model non-minimally coupled to gravity and containing just one additional degree of freedom on top of the Standard Model particle content. This minimalistic scenario predicts a set of measurable consistency relations between the inflationary observables and the dark-energy equation-of-state parameter. We present an alternative derivation of these consistency relations that highlights the connections and differences with the $\alpha$-attractor scenario. We study in how far these constraints allow to distinguish the Higgs-Dilaton model from $\Lambda$CDM and $w$CDM cosmologies. To this end we first analyze existing data sets using a Markov Chain Monte Carlo approach. Second, we perform forecasts for future galaxy surveys using a Fisher matrix approach, both for galaxy clustering and weak lensing probes. Assuming that the best fit values in the different models remain comparable to the present ones, we show that both Euclid- and SKA2-like missions will be able to discriminate a Higgs-Dilaton cosmology from $\Lambda$CDM and $w$CDM.

Modified Gravity theories generally affect the Poisson equation and the gravitational slip (effective anisotropic stress) in an observable way, that can be parameterized by two generic functions ($\eta$ and $\mu$) of time and space. We bin the time dependence of these functions in redshift and present forecasts on each bin for future surveys like Euclid. We consider both Galaxy Clustering and Weak Lensing surveys, showing the impact of the non-linear regime, treated with two different semi-analytical approximations. In addition to these future observables, we use a prior covariance matrix derived from the Planck observations of the Cosmic Microwave Background. Our results show that $\eta$ and $\mu$ in different redshift bins are significantly correlated, but including non-linear scales reduces or even eliminates the correlation, breaking the degeneracy between Modified Gravity parameters and the overall amplitude of the matter power spectrum. We further decorrelate parameters with a Zero-phase Component Analysis and identify which combinations of the Modified Gravity parameter amplitudes, in different redshift bins, are best constrained by future surveys. We also extend the analysis to two particular parameterizations of the time evolution of $\mu$ and $\eta$ and consider, in addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future surveys will be able to constrain the current values of $\eta$ and $\mu$ at the $2-5\%$ level when using only linear scales (wavevector k < 0.15 h/Mpc), depending on the specific time parameterization; sensitivity improves to about $1\%$ when non-linearities are included.

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range $z=0-1.6$ and wave modes below $k=10 \text{h/Mpc}$. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, $\beta^{2}$, with precision smaller than 4\% and all other cosmological parameters better than 1\%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications.

We investigate the formation and dissipation of large scale neutrino structures in cosmologies where the time evolution of dynamical dark energy is stopped by a growing neutrino mass. In models where the coupling between neutrinos and dark energy grows with the value of the scalar cosmon field, the evolution of neutrino lumps depends on the neutrino mass. For small masses the lumps form and dissolve periodically, leaving only a small backreaction of the neutrino structures on the cosmic evolution. This process heats the neutrinos to temperatures much above the photon temperature such that neutrinos acquire again an almost relativistic equation of state. The present equation of state of the combined cosmon-neutrino fluid is very close to -1. In contrast, for larger neutrino masses the lumps become stable. The highly concentrated neutrino structures entail a large backreaction similar to the case of a constant neutrino-cosmon coupling. A present average neutrino mass of around 0.5 eV seems so far compatible with observation. For masses lower than this value, neutrino induced gravitational potentials remain small, making the lumps difficult to detect.