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...

We show that, in the context of patchy reionisation, an accurate description of the angular power spectrum of the kinetic Sunyaev-Zel'dovich (kSZ) effect is not possible with simple scaling relations between the amplitude of the spectrum and global parameters, such as the reionisation midpoint and its duration. We introduce a new parameterisation of this spectrum, based on a novel description of the power spectrum of the free electrons density contrast Pee (k,z) in terms of the reionisation global history and morphology. We directly relate features of the spectrum to the typical ionised bubble size at different stages in the process and, subsequently, to the angular scale at which the patchy kSZ power spectrum reaches its maximum. We successfully calibrated our results on a custom set of advanced radiative hydrodynamical simulations and later found our parameterisation to be a valid description of a wide range of other simulations and, therefore, reionisation physics. In the end, and as long as the global reionisation history is known, two parameters are sufficient to derive the angular power spectrum. Such an innovative framework applied to cosmic microwave background data and combined with 21cm intensity mapping will allow a first consistent detection of the amplitude and shape of the patchy kSZ signal, giving in turn access to the physics of early light sources.

The key probes of the growth of large-scale structure are its rate $f$ and amplitude $\sigma_8$. Redshift space distortions in the galaxy power spectrum allow us to measure only the combination $f\sigma_8$, which can be used to constrain the standard cosmological model or alternatives. By using measurements of the galaxy-galaxy lensing cross-correlation spectrum or of the galaxy bispectrum, it is possible to break the $f\sigma_8$ degeneracy and obtain separate estimates of $f$ and $\sigma_8$ from the same galaxy sample. Currently there are only a handful of such separate measurements, but even this allows for improved constraints on cosmological models. We consider what can be achieved by a future nominal sample that delivers a $\sim 1\%$ constraint on $f$ and $\sigma_8$ separately, compared to the case with a similar precision on the combination $f\sigma_8$. For the six cosmological parameters of $\Lambda$CDM, we find improvements of $\sim\! 5$--$50\%$ on their constraints. For modified gravity models in the Horndeski class, the improvements on these standard parameters are $\sim\! 0$--$15\%$. However, the precision on the sum of neutrino masses improves by 65\% and there is a significant increase in the precision on the background and perturbation Horndeski parameters.

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]

We perform the first test of Dark Matter (DM) stress-energy evolution through cosmic history, using Cosmic Microwave Background measurements supplemented with Baryon Acoustic Oscillation data and the Hubble Space Telescope key project data. We constrain the DM equation of state (EoS) in 8 redshift bins, and its sound speed and (shear) viscosity in 9 redshift bins, finding no convincing evidence for non-$\Lambda$CDM values in any of the redshift bins. Despite this enlarged parameter space, the sound speed and viscosity are constrained relatively well at late times (due to the inclusion of CMB lensing), whereas the EoS is most strongly constrained around recombination. These results constrain for the first time the level of ''coldness'' required of DM across various cosmological epochs at both the background and perturbative levels. We show that simultaneously allowing time dependence for both the EoS and sound speed parameters shifts the posterior of the DM abundance before recombination to a higher value, whilst keeping the present day DM abundance similar to the $\Lambda$CDM value. This shifts the posterior for the present day Hubble constant compared to $\Lambda$CDM, suggesting that DM with time-dependent parameters is well-suited to explore possible solutions to persistent tensions within the $\Lambda$CDM 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.

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.

We explore the use of the cosmic infrared background as a tracer of the LSS for cross-correlating with the CMB and exploit the ISW. We use the improved linear CIB model of Maniyar et al (2018) and derive the theoretical CIBxISW cross-correlation for Planck HFI frequencies 217, 353, 545 and 857 GHz and IRAS 3000 GHz. We predict a positive cross-correlation between the CIB and CMB whose amplitude decreases rapidly at small scales. We perform a signal-to-noise ratio (SNR) analysis on this cross-correlation. In the ideal case of the cross-correlation obtained over 70% (40%) of the sky with no residual contaminants (e.g. galactic dust) in maps, the SNR ranges from 4.2-5.6 (3.2-4.3) with the highest for 857 GHz. A Fisher matrix analysis shows that an ISW signal detected with such high SNR on the 40% sky can improve the constraints on the cosmological parameters considerably; constraints on the equation of state of the dark energy are improved by 80%. We then perform a more realistic analysis with the effect of residual galactic dust contamination in CIB maps. We calculate the dust power spectra for several frequencies and sky fractions which dominate over CIB at lower multipoles we are interested in. Considering conservative 10% residual level of galactic dust in the CIB power spectra, we find that the SNR drops drastically making ISW detection difficult. To check the capability of current maps to detect ISW via this method, we measure the cross-correlation of the CIB and the CMB Planck maps on so called GASS field covering an area of 11% in the southern hemisphere. We find that with such a small sky fraction and dust residuals present in CIB maps, we do not detect ISW signal and the measured signal is consistent with zero. In order not to degrade the SNR for the ISW measurement by more than 10% on the 40% sky, we find that the dust needs to be cleaned up to 0.01% level on the power spectrum.

Cold Dark Matter (CDM) is a crucial constituent of the current concordance cosmological model. Having a vanishing equation of state (EoS), its energy density scales with the inverse cosmic volume and is thus uniquely described by a single number, its present abundance. We test the inverse cosmic volume law for Dark Matter (DM) by allowing its EoS to vary independently in eight redshift bins in the range $z=10^5$ and $z=0$. We use the latest measurements of the Cosmic Microwave Background radiation from the Planck satellite and supplement them with Baryon Acoustic Oscillation (BAO) data from the 6dF and SDSS-III BOSS surveys, and with the Hubble Space Telescope (HST) key project data. We find no evidence for nonzero EoS in any of the eight redshift bins. With Planck data alone, the DM abundance is most strongly constrained around matter-radiation equality $\omega^{\rm eq}_g = 0.1193^{+0.0036}_{-0.0035}$ (95% c.l.), whereas its present day value is more weakly constrained $\omega^{(0)}_g = 0.16^{+0.12}_{-0.10}$ (95% c.l.). Adding BAO or HST data does not significantly change the $\omega^{\rm eq}_g$ constraint, while $\omega^{(0)}_g$ tightens to $0.160^{+0.069}_{-0.065} $ (95% c.l.) and $0.124^{+0.081}_{-0.067}$ (95% c.l.) respectively. Our results constrain for the first time the level of "coldness" required of the DM across various cosmological epochs and show that the DM abundance is strictly positive at all times.

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.

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.

Motivated by the current constraints on the epoch of reionisation from recent cosmic microwave background observations, ionising background measurements of star-forming galaxies, and low redshifts line-of-sight probes, we propose a new data-motivated parameterisation of the history of the average ionisation fraction. This parameterisation describes a flexible redshift-asymmetric reionisation process in two regimes that is capable of fitting all the current constraints.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.