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 address key points for an efficient implementation of likelihood codes for modern weak lensing large-scale structure surveys. Specifically, we focus on the joint weak lensing convergence power spectrum-bispectrum probe and we tackle the numerical challenges required by a realistic analysis. Under the assumption of (multivariate) Gaussian likelihoods, we have developed a high performance code that allows highly parallelised prediction of the binned tomographic observables and of their joint non-Gaussian covariance matrix accounting for terms up to the 6-point correlation function and super-sample effects. This performance allows us to qualitatively address several interesting scientific questions. We find that the bispectrum provides an improvement in terms of signal-to-noise ratio (S/N) of about 10% on top of the power spectrum, making it a non-negligible source of information for future surveys. Furthermore, we are capable to test the impact of theoretical uncertainties in the halo model used to build our observables; with presently allowed variations we conclude that the impact is negligible on the S/N. Finally, we consider data compression possibilities to optimise future analyses of the weak lensing bispectrum. We find that, ignoring systematics, 5 equipopulated redshift bins are enough to recover the information content of a Euclid-like survey, with negligible improvement when increasing to 10 bins. We also explore principal component analysis and dependence on the triangle shapes as ways to reduce the numerical complexity of the problem.

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.

Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to $k_{\rm max} = 4 h/{\rm Mpc}$, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, $f(R)$ gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as \emph{Euclid}. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.

As galaxy surveys become more precise and push to smaller scales, the need for accurate covariances beyond the classical Gaussian formula becomes more acute. Here, I investigate the analytical implementation and impact of non-Gaussian covariance terms that I previously derived for galaxy clustering. Braiding covariance is such a class of terms and it gets contribution both from in-survey and super-survey modes. I present an approximation for braiding covariance which speeds up the numerical computation. I show that including braiding covariance is a necessary condition for including other non-Gaussian terms: the in-survey 2-, 3- and 4-halo covariance, which yield covariance matrices with negative eigenvalues if considered on their own. I then quantify the impact on parameter constraints, with forecasts for a Euclid-like survey. Compared to the Gaussian case, braiding and in-survey covariances significantly increase the error bars on cosmological parameters, in particular by 50% for w. The Halo Occupation Distribution (HOD) error bars are also affected between 12% and 39%. Accounting for super-sample covariance (SSC) also increases parameter errors, by 90% for w and between 7% and 64% for HOD. In total, non-Gaussianity increases the error bar on w by 120% (between 15% and 80% for other cosmological parameters), and the error bars on HOD parameters between 17% and 85%. Accounting for the 1-halo trispectrum term on top of SSC is not sufficient for capturing the full non-Gaussian impact: braiding and the rest of in-survey covariance have to be accounted for. Finally, I discuss why the inclusion of non-Gaussianity generally eases up parameter degeneracies, making cosmological constraints more robust to astrophysical uncertainties. The data and a Python notebook reproducing the results and plots of the article are available at \url{https://github.com/fabienlacasa/BraidingArticle}. [Abridged]

Interest rises to exploit the full shape information of the galaxy power spectrum, as well as pushing analyses to smaller non-linear scales. Here I use the halo model to quantify the information content in the tomographic angular power spectrum of galaxies, for future high resolution surveys : Euclid and SKA2. I study how this information varies as a function of the scale cut applied, either with angular cut $\ell_{max}$ or physical cut kmax. For this, I use analytical covariances with the most complete census of non-Gaussian terms, which proves critical. I find that the Fisher information on most cosmological and astrophysical parameters follows a striking behaviour. Beyond the perturbative regime we first get decreasing returns : the information keeps rising but the slope slows down until reaching a saturation. The location of this plateau is a bit beyond the reach of current modeling methods : k $\sim$ 2 Mpc$^{-1}$ and slightly depends on the parameter and redshift bin considered. I explain the origin of this plateau, which is due to non-linear effects both on the power spectrum, and more importantly on non-Gaussian covariance terms. Then, pushing further on I find that information rises again in the highly non-linear regime. I find that the cosmological information in this small scale miracle can indeed be disentangled from astrophysical information and yield large improvements. Pushing SKA2 analysis from kmax=1 Mpc$^{-1}$ to kmax=10 Mpc$^{-1}$ can improve the error bar on $\sigma_8$ by a factor 9 and the error bar on the Dark Energy equation of state $w_0$ by a factor 5. Finally I show that high order statistics beyond the power spectrum should yield further significant improvements in this regime, with the improvements increasing when pushing kmax. Data and notebooks reproducing all plots and results will be made available at \url{https://github.com/fabienlacasa/SmallScaleMiracle}

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.

We propose a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$\times$galaxy $-$ galaxy$\times$IM; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As a pure cross-correlation estimator, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. For cosmic variance dominated surveys the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of DES cross-correlated with HIRAX, the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

We present a numerically cheap approximation to super-sample covariance (SSC) of large scale structure cosmological probes, first in the case of angular power spectra. It necessitates no new elements besides those used for the prediction of the considered probes, thus relieving analysis pipelines from having to develop a full SSC modeling, and reducing the computational load. The approximation is asymptotically exact for fine redshift bins $\Delta z \rightarrow 0$. We furthermore show how it can be implemented at the level of a Gaussian likelihood or a Fisher matrix forecast, as a fast correction to the Gaussian case without needing to build large covariance matrices. Numerical application to a Euclid-like survey show that, compared to a full SSC computation, the approximation recovers nicely the signal-to-noise ratio as well as Fisher forecasts on cosmological parameters of the $w$CDM cosmological model. Moreover it allows for a fast prediction of which parameters are going to be the most affected by SSC and at which level. In the case of photometric galaxy clustering with Euclid-like specifications, we find that $\sigma_8$, $n_s$ and the dark energy equation of state $w$ are particularly heavily affected. We finally show how to generalize the approximation for probes other than angular spectra (correlation functions, number counts and bispectra), and at the likelihood level, allowing for the latter to be non-Gaussian if needs be. We release publicly a Python module allowing to implement the SSC approximation, as well as a notebook reproducing the plots of the article, at https://github.com/fabienlacasa/PySSC

We present angular diameter distance measurements obtained by locating the BAO scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1318 deg$^2$ with $0.6 < z_{\rm photo} < 1$ and a typical redshift uncertainty of $0.03(1+z)$. This sample was selected, as fully described in a companion paper, using a color/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the co-moving transverse separation, and spherical harmonics. Further, we compare results obtained from template based and machine learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, $D_A$, at the effective redshift of our sample divided by the true physical scale of the BAO feature, $r_{\rm d}$. We obtain close to a 4 per cent distance measurement of $D_A(z_{\rm eff}=0.81)/r_{\rm d} = 10.75\pm 0.43 $. These results are consistent with the flat $\Lambda$CDM concordance cosmological model supported by numerous other recent experimental results. All data products are publicly available here: https://des.ncsa.illinois.edu/releases/y1a1/bao

Magnetized plasmas within halos of galaxies leave their footprint on the polarized anisotropies of the cosmic microwave background. The two dominant effects for astrophysical halos are Faraday rotation generating rotation of the plane of linear polarization, and Faraday conversion inducing a leakage from linear polarization to circular polarization. We revisit these sources of secondary anisotropies by computing the angular power spectra of the Faraday rotation angle and of the Faraday conversion rate by the large scale structures. To this end, we use the halo model and we pay special attention to the impact of magnetic field projections. Assuming magnetic fields of halos to be uncorrelated, we found a vanishing 2-halo term, and angular power spectra peaking at multipoles $\ell\sim10^4$. The Faraday rotation angle is dominated by the contribution of thermal electrons. For the Faraday conversion rate, we found that both thermal electrons and relativistic, non-thermal electrons contribute equally in the most optimistic case for the density and Lorentz factor of relativistic electrons, while in more pessimistic cases the thermal electrons give the dominant contribution. Assuming the magnetic field to be independent of the halo mass, the angular power spectra for both effects roughly scale with the amplitude of matter perturbations as $\sim\sigma_8^3$, and with a very mild dependence with the density of cold dark matter. Introducing a dependence of the magnetic field strength with the halo mass leads to an increase of the scaling with the amplitude of matter fluctuations, up to $\sim\sigma_8^{9.5}$ for Faraday rotation and $\sim\sigma_8^{15}$ for Faraday conversion for a magnetic field strength scaling linearly with the halo mass.

We explore the impact of an update to the typical approximation for the shape noise term in the analytic covariance matrix for cosmic shear experiments that assumes the absence of survey boundary and mask effects. We present an exact expression for the number of galaxy pairs in this term based on the the survey mask, which leads to more than a factor of three increase in the shape noise on the largest measured scales for the Kilo-Degree Survey (KIDS-450) real-space cosmic shear data. We compare the result of this analytic expression to several alternative methods for measuring the shape noise from the data and find excellent agreement. This update to the covariance resolves any internal model tension evidenced by the previously large cosmological best-fit $\chi^2$ for the KiDS-450 cosmic shear data. The best-fit $\chi^2$ is reduced from 161 to 121 for 118 degrees of freedom. We also apply a correction to how the multiplicative shear calibration uncertainty is included in the covariance. This change, along with a previously known update to the reported effective angular values of the data vector, jointly shift the inferred amplitude of the correlation function to higher values. We find that this improves agreement of the KiDS-450 cosmic shear results with Dark Energy Survey Year 1 and Planck results.

We use data from the first-year (Y1) observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature using a template-fitting method. We test our methodology in a sample of 1800 DES Y1-like mock catalogs. The APS is measured with the pseudo-$C_\ell$ method, using pixelized maps constructed from the mock catalogs and the DES mask. The covariance matrix of the $C_\ell$'s in these tests are also obtained from the mock catalogs. We use templates to model the measured spectra and estimate template parameters firstly from the $C_\ell$'s of the mocks using two different methods, a maximum likelihood estimator and a MCMC, finding consistent results with a good reduced $\chi^2$. Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. After these tests on mocks, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalog comprises galaxies within an effective area of 1318 deg$^2$ and $0.6<z<1.0$. We fit the observed spectra with our optimized templates, considering models with and without BAO features. We find that the DES Y1 data favors a model with BAO at the $2.6\,\sigma$ C.L. with a best-fit shift parameter of $\alpha=1.023\pm 0.047$. However, the goodness-of-fit is somewhat poor, with $\chi^2/$(dof) = 1.49. We identify a possible cause of this issue and show that using a theoretical covariance matrix obtained from $C_\ell$'s that are better adjusted to data results in an improved value of $\chi^2/$(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Our results correspond to a distance measurement of $D_A(z_{\rm eff} = 0.81)/r_d = 10.65 \pm 0.49$, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.

As the determination of density fluctuations becomes more precise with larger surveys, it becomes more important to account for the increased covariance due to the non-linearity of the field. Here I have focussed on the galaxy density, with analytical prediction of the non-Gaussianity using the halo model coupled with standard perturbation theory in real space. I carried out an exact and exhaustive derivation of all tree-level terms of the non-Gaussian covariance of the galaxy $C_\ell$, with the computation developed up to the third order in perturbation theory and local halo bias, including the non-local tidal tensor effect. A diagrammatic method was used to derive the involved galaxy 3D trispectra, including shot-noise contributions. The projection to the angular covariance was derived in all trispectra cases with and without Limber's approximation, with the formulae being of potential interest for other observables than galaxies. The effect of subtracting shot-noise from the measured spectrum is also discussed, and does simplify the covariance, though some non-Gaussian shot-noise terms still remain. I make the link between this complete derivation and partial terms which have been used previously in the literature, including super-sample covariance (SSC). I uncover a wealth of additional terms which were not previously considered, including a whole new class which I dub braiding terms as it contains multipole-mixing kernels. The importance of all these new terms is discussed with analytical arguments. I find that they become comparable to, if not bigger than, SSC if the survey is large or deep enough to probe scales comparable with the matter-radiation equality $k_\mathrm{eq}$. A short self-contained summary of the equations is provided in Section 9 for the busy reader, ready to be implemented numerically for analysis of current and future galaxy surveys.

Super-sample covariance (SSC) is the dominant source of statistical error on large scale structure (LSS) observables for both current and future galaxy surveys. In this work, we concentrate on the SSC of cluster counts, also known as sample variance, which is particularly useful for the self-calibration of the cluster observable-mass relation; our approach can similarly be applied to other observables, such as galaxy clustering and lensing shear. We first examined the accuracy of two analytical approximations proposed in the literature for the flat sky limit, finding that they are accurate at the 15% and 30-35% level, respectively, for covariances of counts in the same redshift bin. We then developed a harmonic expansion formalism that allows for the prediction of SSC in an arbitrary survey mask geometry, such as large sky areas of current and future surveys. We show analytically and numerically that this formalism recovers the full sky and flat sky limits present in the literature. We then present an efficient numerical implementation of the formalism, which allows fast and easy runs of covariance predictions when the survey mask is modified. We applied our method to a mask that is broadly similar to the Dark Energy Survey footprint, finding a non-negligible negative cross-z covariance, i.e. redshift bins are anti-correlated. We also examined the case of data removal from holes due to, for example bright stars, quality cuts, or systematic removals, and find that this does not have noticeable effects on the structure of the SSC matrix, only rescaling its amplitude by the effective survey area. These advances enable analytical covariances of LSS observables to be computed for current and future galaxy surveys, which cover large areas of the sky where the flat sky approximation fails.

We give an analytical interpretation of how subsample-based internal covariance estimators lead to biased estimates of the covariance, due to underestimating the super-sample covariance (SSC). This includes the jackknife and bootstrap methods as estimators for the full survey area, and subsampling as an estimator of the covariance of subsamples. The limitations of the jackknife covariance have been previously presented in the literature because it is effectively a rescaling of the covariance of the subsample area. However we point out that subsampling is also biased, but for a different reason: the subsamples are not independent, and the corresponding lack of power results in SSC underprediction. We develop the formalism in the case of cluster counts that allows the bias of each covariance estimator to be exactly predicted. We find significant effects for a small-scale area or when a low number of subsamples is used, with auto-redshift biases ranging from 0.4% to 15% for subsampling and from 5% to 75% for jackknife covariance estimates. The cross-redshift covariance is even more affected; biases range from 8% to 25% for subsampling and from 50% to 90% for jackknife. Owing to the redshift evolution of the probe, the covariances cannot be debiased by a simple rescaling factor, and an exact debiasing has the same requirements as the full SSC prediction. These results thus disfavour the use of internal covariance estimators on data itself or a single simulation, leaving analytical prediction and simulations suites as possible SSC predictors.

We present the methodology for and detail the implementation of the Dark Energy Survey (DES) 3x2pt DES Year 1 (Y1) analysis, which combines configuration-space two-point statistics from three different cosmological probes: cosmic shear, galaxy-galaxy lensing, and galaxy clustering, using data from the first year of DES observations. We have developed two independent modeling pipelines and describe the code validation process. We derive expressions for analytical real-space multi-probe covariances, and describe their validation with numerical simulations. We stress-test the inference pipelines in simulated likelihood analyses that vary 6-7 cosmology parameters plus 20 nuisance parameters and precisely resemble the analysis to be presented in the DES 3x2pt analysis paper, using a variety of simulated input data vectors with varying assumptions. We find that any disagreement between pipelines leads to changes in assigned likelihood $\Delta \chi^2 \le 0.045$ with respect to the statistical error of the DES Y1 data vector. We also find that angular binning and survey mask do not impact our analytic covariance at a significant level. We determine lower bounds on scales used for analysis of galaxy clustering (8 Mpc$~h^{-1}$) and galaxy-galaxy lensing (12 Mpc$~h^{-1}$) such that the impact of modeling uncertainties in the non-linear regime is well below statistical errors, and show that our analysis choices are robust against a variety of systematics. These tests demonstrate that we have a robust analysis pipeline that yields unbiased cosmological parameter inferences for the flagship 3x2pt DES Y1 analysis. We emphasize that the level of independent code development and subsequent code comparison as demonstrated in this paper is necessary to produce credible constraints from increasingly complex multi-probe analyses of current data.

The Sunyaev-Zel'dovich (SZ) effect is a powerful probe of the evolution of structures in the universe, and is thus highly sensitive to cosmological parameters $\sigma_8$ and $\Omega_m$, though its power is hampered by the current uncertainties on the cluster mass calibration. In this analysis we revisit constraints on these cosmological parameters as well as the hydrostatic mass bias, by performing (i) a robust estimation of the tSZ power-spectrum, (ii) a complete modeling and analysis of the tSZ bispectrum, and (iii) a combined analysis of galaxy clusters number count, tSZ power spectrum, and tSZ bispectrum. From this analysis, we derive as final constraints $\sigma_8 = 0.79 \pm 0.02$, $\Omega_{\rm m} = 0.29 \pm 0.02$, and $(1-b) = 0.71 \pm 0.07$. These results favour a high value for the hydrostatic mass bias compared to numerical simulations and weak-lensing based estimations. They are furthermore consistent with both previous tSZ analyses, CMB derived cosmological parameters, and ancillary estimations of the hydrostatic mass bias.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

The abundance of clusters and the clustering of galaxies are two of the important cosmological probes for current and future large scale surveys of galaxies, such as the Dark Energy Survey. In order to combine them one has to account for the fact that they are not independent quantities, since they probe the same density field. It is important to develop a good understanding of their correlation in order to extract parameter constraints. We present a detailed modelling of the joint covariance matrix between cluster number counts and the galaxy angular power spectrum. We employ the framework of the halo model complemented by a Halo Occupation Distribution model (HOD). We demonstrate the importance of accounting for non-Gaussianity to produce accurate covariance predictions. Indeed, we show that the non-Gaussian covariance becomes dominant at small scales, low redshifts or high cluster masses. We discuss in particular the case of the super-sample covariance (SSC), including the effects of galaxy shot-noise, halo second order bias and non-local bias. We demonstrate that the SSC obeys mathematical inequalities and positivity. Using the joint covariance matrix and a Fisher matrix methodology, we examine the prospects of combining these two probes to constrain cosmological and HOD parameters. We find that the combination indeed results in noticeably better constraints, with improvements of order 20\% on cosmological parameters compared to the best single probe, and even greater improvement on HOD parameters, with reduction of error bars by a factor 1.4-4.8. This happens in particular because the cross-covariance introduces a synergy between the probes on small scales. We conclude that accounting for non-Gaussian effects is required for the joint analysis of these observables in galaxy surveys.

In the Halo Model, galaxies are hosted by dark matter halos, while the halos themselves are biased tracers of the underlying matter distribution. Measurements of galaxy correlation functions include contributions both from galaxies in different halos, and from galaxies in the same halo (the so-called 1-halo terms). We show that, for highly biased tracers, the 1-halo term of the power spectrum obeys a steep scaling relation in terms of bias. We also show that the 1-halo term of the trispectrum has a steep scaling with bias. The steepness of these scaling relations is such that the 1-halo terms can become key contributions to the $n$-point correlation functions, even at large scales. We interpret these results through analytical arguments and semi-analytical calculations in terms of the statistical properties of halos.

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.

We have constructed all-sky y-maps of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 30 to 857 GHz frequency channel maps from the Planck satellite survey. These reconstructed y-maps are delivered as part of the Planck 2015 release. The y-maps are characterised in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales and CIB and extragalactic point sources at small angular scales. Specific masks are defined to minimize foreground residuals and systematics. Using these masks we compute the y-map angular power spectrum and higher order statistics. From these we conclude that the y-map is dominated by tSZ signal in the multipole range, 20-600. We compare the measured tSZ power spectrum and higher order statistics to various physically motivated models and discuss the implications of our results in terms of cluster physics and cosmology.

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Present and future large scale surveys offer promising probes of cosmology. For example the Dark Energy Survey (DES) is forecast to detect ~300 millions galaxies and thousands clusters up to redshift ~1.3. I here show ongoing work to combine two probes of large scale structure : cluster number counts and galaxy 2-point function (in real or harmonic space). The halo model (coupled to a Halo Occupation Distribution) can be used to model the cross-covariance between these probes, and I introduce a diagrammatic method to compute easily the different terms involved. Furthermore, I compute the joint non-Gaussian likelihood, using the Gram-Charlier series. Then I show how to extend the methods of Bayesian hyperparameters to Poissonian distributions, in a first step to include them in this joint likelihood.

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.

[Abridged] This PhD thesis was defended in September 2013 at Universit\'e Paris-Sud, France. It studies the non-Gaussianity (NG) of extragalactic foregrounds to the Cosmic Microwave Background (CMB). This is particularly important for the analysis of CMB data from the Planck satellite, in which I was involved. I studied in particular the extragalactic foregrounds tracing the Large Scale Structure (LSS) of the universe : radio and infrared point-sources and the thermal Sunyaev-Zel'dovich effect (tSZ). The motivations and summary of this thesis cannot be shortened into the length constraints of this abstract and are thus described in an inner abstract and introduction. In a broad summary, this thesis is organised in the following parts : statistics of random variables and fields, spherical analysis in particular of polyspectra, some cosmology (inflation and its NG, CMB, LSS and halo model), extragalactic foregrounds (tSZ, point-sources), and measuring the NG of these foregrounds with Planck. More specifically and in a disorganised manner, the following were most of my contributions in the field : developments on the effect of masking on the angular bispectrum, introduction of an efficient bispectrum representation, derivation of the bispectrum covariance, generation of non-gaussian simulations, projection of 3D Fourier polyspectra onto 2D spherical ones, treatment of shot-noise terms with the halo model, introduction of a diagrammatic formalism for elegant and simple derivation of galaxy polyspectra, halo model description of the bispectrum of the tSZ effect and of the Cosmic Infrared Background (CIB), phenomenological prescription of the NG of clustered sources and corresponding KSW-like estimator, effect of point-sources NG for Planck, measurement of the CIB bispectrum with Planck data, measurement of the tSZ bispectrum with Planck data and cosmological constraints with it.

Using a full analytical computation of the bispectrum based on the halo model together with the halo occupation number, we derive the bispectrum of the cos- mic infrared background (CIB) anisotropies that trace the clustering of dusty-star- forming galaxies. We focus our analysis on wavelengths in the far-infrared and the sub-millimeter typical of the Planck/HFI and Herschel/SPIRE instruments, 350, 550, 850, and 1380 um. We explore the bispectrum behaviour as a function of several models of evolution of galaxies and show that it is strongly sensitive to that ingredient. Contrary to the power spectrum, the bispectrum, at the four wavelengths, seems dominated by low redshift galaxies. Such a contribution can be hardly limited by applying low flux cuts. We also discuss the contributions of halo mass as a function of the redshift and the wavelength, recovering that each term is sensitive to a different mass range. Furthermore, we show that the CIB bispectrum is a strong contaminant of the Cosmic Microwave Background bispectrum at 850 um and higher. Finally, a Fisher analysis of the power spectrum, bispectrum alone and of the combination of both shows that degeneracies on the HOD parameters are broken by including the bispectrum information, leading to tight constraints even when including foreground residuals.

We present the first halo model based description of the Cosmic Infrared Background (CIB) non-Gaussianity (NG) that is fully parametric. To this end, we introduce, for the first time, a diagrammatic method to compute high order polyspectra of the 3D galaxy density field. It allows an easy derivation and visualisation of the different terms of the polyspectrum. We apply this framework to the power spectrum and bispectrum, and we show how to project them on the celestial sphere in the purpose of the application to the CIB angular anisotropies. Furthermore, we show how to take into account the particular case of the shot noise terms in that framework. Eventually, we compute the CIB angular bispectrum at 857 GHz and study its scale and configuration dependencies, as well as its variations with the halo occupation distribution parameters. Compared to a previously proposed empirical prescription, such physically motivated model is required to describe fully the CIB anisotropies bispectrum. Finally, we compare the CIB bispectrum with the bispectra of other signals potentially present at microwave frequencies, which hints that detection of CIB NG should be possible above 220 GHz.

The Planck nominal mission cosmic microwave background (CMB) maps yield unprecedented constraints on primordial non-Gaussianity (NG). Using three optimal bispectrum estimators, separable template-fitting (KSW), binned, and modal, we obtain consistent values for the primordial local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result fNL^local= 2.7+/-5.8, fNL^equil= -42+/-75, and fNL^ortho= -25+\-39 (68% CL statistical). NG is detected in the data; using skew-C_l statistics we find a nonzero bispectrum from residual point sources, and the ISW-lensing bispectrum at a level expected in the LambdaCDM scenario. The results are based on comprehensive cross-validation of these estimators on Gaussian and non-Gaussian simulations, are stable across component separation techniques, pass an extensive suite of tests, and are confirmed by skew-C_l, wavelet bispectrum and Minkowski functional estimators. Beyond estimates of individual shape amplitudes, we present model-independent, 3-dimensional reconstructions of the Planck CMB bispectrum and thus derive constraints on early-Universe scenarios that generate primordial NG, including general single-field models of inflation, excited initial states (non-Bunch-Davies vacua), and directionally-dependent vector models. We provide an initial survey of scale-dependent feature and resonance models. These results bound both general single-field and multi-field model parameter ranges, such as the speed of sound, c_s \geq 0.02 (95% CL), in an effective field theory parametrization, and the curvaton decay fraction r_D \geq 0.15 (95% CL). The Planck data significantly limit the viable parameter space of the ekpyrotic/cyclic scenarios. The amplitude of the 4-point function in the local model tauNL < 2800 (95% CL). These constraints represent the highest precision tests to date of physical mechanisms for the origin of cosmic structure.

We present new measurements of CIB anisotropies using Planck. Combining HFI data with IRAS, the angular auto- and cross frequency power spectrum is measured from 143 to 3000 GHz, and the auto-bispectrum from 217 to 545 GHz. The total areas used to compute the CIB power spectrum and bispectrum are about 2240 and 4400 deg^2, respectively. After careful removal of the contaminants, and a complete study of systematics, the CIB power spectrum and bispectrum are measured with unprecedented signal to noise ratio from angular multipoles ell~150 to 2500, and ell~130 to 1100, respectively. Two approaches are developed for modelling CIB power spectrum anisotropies. The first approach takes advantage of the unique measurements by Planck at large angular scales, and models only the linear part of the power spectrum, with a mean bias of dark matter halos hosting dusty galaxies at a given redshift weighted by their contribution to the emissivities. The second approach is based on a model that associates star-forming galaxies with dark matter halos and their subhalos, using a parametrized relation between the dust-processed infrared luminosity and (sub-)halo mass. The two approaches simultaneously fit all auto- and cross- power spectra very well. We find that the star formation history is well constrained up to z~2. However, at higher redshift, the accuracy of the star formation history measurement is strongly degraded by the uncertainty in the spectral energy distribution of CIB galaxies. We also find that CIB galaxies have warmer temperatures as redshift increases. The CIB bispectrum is steeper than that expected from the power spectrum, although well fitted by a power law; this gives some information about the contribution of massive halos to the CIB bispectrum.

The multi-frequency capability of the Planck satellite provides information both on the integrated history of star formation (via the cosmic infrared background, or CIB) and on the distribution of dark matter (via the lensing effect on the cosmic microwave background, or CMB). The conjunction of these two unique probes allows us to measure directly the connection between dark and luminous matter in the high redshift (1 < z <3) Universe. We use a three-point statistic optimized to detect the correlation between these two tracers. Following a thorough discussion of possible contaminants and a suite of consistency tests, using lens reconstructions at 100, 143 and 217 GHz and CIB measurements at 100-857 GHz, we report the first detection of the correlation between the CIB and CMB lensing. The well matched redshift distribution of these two signals leads to a detection significance with a peak value of 42 \sigma at 545 GHz and a correlation as high as 80% across these two tracers. Our full set of multi-frequency measurements (both CIB auto- and CIB-lensing cross-spectra) are consistent with a simple halo-based model, with a characteristic mass scale for the halos hosting CIB sources of log_{10}(M/M_sun) = 10.5 \pm 0.6. Leveraging the frequency dependence of our signal, we isolate the high redshift contribution to the CIB, and constrain the star formation rate (SFR) density at z>1. We measure directly the SFR density with around 2 sigma significance for three redshift bins between z=1 and 7, thus opening a new window into the study of the formation of stars at early times.

We devise a fast and optimal estimator for the amplitude of the bispectrum of clustered Infrared (IR) point-sources. We show how this estimator can account for the cases of partial sky coverage and inhomogeneous noise. Expected detection significance are presented in terms of signal-to-noise, finding that the IR bispectrum will realistically be undetectable below 220 GHz with a Planck-like experiment; on the contrary detection may be achieved at, or above, 220 GHz especially if the CMB is removed. We also show how this estimator can be combined with estimators of radio and CMB non-Gaussianity to build up joint robust constraints. On the one hand, we find that, for a Planck-like experiment, CMB non-Gaussianity estimation can be decoupled from point-source contributions, unless few sources are masked. On the other hand, we find that the estimation of radio and IR non-Gaussianity are strongly coupled, which diminishes their separate detection significance.

The population of compact extragalactic sources contribute to the non-Gaussianity at Cosmic Microwave Background frequencies. We study their non-Gaussianity using publicly available full-sky simulations. We introduce a parametrisation to visualise efficiently the bispectrum and we describe the scale and frequency dependences of the bispectrum of radio and IR point-sources. We show that the bispectrum is well fitted by an analytical prescription. We find that the clustering of IR sources enhances their non-Gaussianity by several orders of magnitude, and that their bispectrum peaks in the squeezed triangles. Examining the impact of these sources on primordial non-Gaussianity estimation, we find that radio sources yield an important positive bias to local fNL at low frequencies but this bias is efficiently reduced by masking detectable sources. IR sources produce a negative bias at high frequencies, which is not dimmed by the masking, as their clustering is dominated by faint sources.

This study, using publicly available simulations, focuses on the characterisation of the non-Gaussianity produced by radio point sources and by infrared (IR) sources in the frequency range of the cosmic microwave background from 30 to 350 GHz. We propose a simple prescription to infer the angular bispectrum from the power spectrum of point sources considering independent populations of sources, with or without clustering. We test the accuracy of our prediction using publicly available all-sky simulations of radio and IR sources and find very good agreement. We further characterise the configuration dependence and the frequency behaviour of the IR and radio bispectra. We show that the IR angular bispectrum peaks for squeezed triangles and that the clustering of IR sources enhances the bispectrum values by several orders of magnitude at scales l \sim 100. At 150 GHz the bispectrum of IR sources starts to dominate that of radio sources on large angular scales, and it dominates over the whole multipole range at 350 GHz. Finally, we compute the bias on f_NL induced by radio and IR sources. We show that the positive bias induced by radio sources is significantly reduced by masking the sources. We also show, for the first time, that the form of the IR bispectrum mimics a primordial 'local' bispectrum f_NL. The IR sources produce a negative bias which becomes important for Planck-like resolution and at high frequencies (Delta f_NL ~ -6 at 277 GHz and Delta f_NL \sim -60-70 at 350 GHz). Most of the signal being due to the clustering of faint IR sources, the bias Delta f_NL^IR is not reduced by masking sources above a flux limit and may, in some cases, even be increased due to the reduction of the shot-noise term.