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.

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

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]

The detection of gravitational waves (GWs) propagating through cosmic structures can provide invaluable information on the geometry and content of our Universe, as well as on the fundamental theory of gravity. In order to test possible departures from General Relativity, it is essential to analyse, in a modified gravity setting, how GWs propagate through a perturbed cosmological space-time. Working within the framework of geometrical optics, we develop tools to address this topic for a broad class of scalar-tensor theories, including scenarios with non-minimal, derivative couplings between scalar and tensor modes. We determine the corresponding evolution equations for the GW amplitude and polarization tensor. The former satisfies a generalised evolution equation that includes possible effects due to a variation of the effective Planck scale; the latter can fail to be parallely transported along GW geodesics unless certain conditions are satisfied. We apply our general formulas to specific scalar-tensor theories with unit tensor speed, and then focus on GW propagation on a perturbed space-time. We determine corrections to standard formulas for the GW luminosity distance and for the evolution of the polarization tensor, which depend both on modified gravity and on the effects of cosmological perturbations. Our results can constitute a starting point to disentangle among degeneracies from different sectors that can influence GW propagation through cosmological space-times.

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.

Cosmic voids offer an extraordinary opportunity to study the effects of massive neutrinos on cosmological scales. Because they are freely streaming, neutrinos can penetrate the interior of voids more easily than cold dark matter or baryons, which makes their relative contribution to the mass budget in voids much higher than elsewhere in the Universe. In simulations it has recently been shown how various characteristics of voids in the matter distribution are affected by neutrinos, such as their abundance, density profiles, dynamics, and clustering properties. However, the tracers used to identify voids in observations (e.g., galaxies or halos) are affected by neutrinos as well, and isolating the unique neutrino signatures inherent to voids becomes more difficult. In this paper we make use of the DEMNUni suite of simulations to investigate the clustering bias of voids in Fourier space as a function of their core density and compensation. We find a clear dependence on the sum of neutrino masses that remains significant even for void statistics extracted from halos. In particular, we observe that the amplitude of the linear void bias increases with neutrino mass for voids defined in dark matter, whereas this trend gets reversed and slightly attenuated when measuring the relative void-halo bias using voids identified in the halo distribution. Finally, we argue how the original behaviour can be restored when considering observations of the total matter distribution (e.g. via weak lensing), and comment on scale-dependent effects in the void bias that may provide additional information on neutrinos in the future.

We test a theoretical description of the void size distribution function against direct estimates from halo catalogues of the DEMNUni suite of large cosmological simulations. Besides standard $\Lambda$CDM, we consider deviations of the dark energy equation of state from $w=-1$, corresponding to four combinations in the popular Chevallier-Polarski-Linder parametrisation: $w_0=-0.9;-1.1$, $w_a=-0.3;0.3$. The theoretical void size function model, relying on the Sheth & van de Weygaert double barrier excursion set formalism, provides an accurate description of the simulation measurements for the different dark energy models considered, within the statistical errors. The model remains accurate for any value of the threshold for void formation $\delta_\mathrm{v}$. Its robust consistency with simulations demonstrates that the theoretical void size function can be applied to real data as a sensitive tool to constrain dark energy.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

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.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects depend on the observed void tracers. As the neutrino mass increases, the number of small voids traced by cold dark matter particles increases and the number of large voids decreases. Surprisingly, when massive, highly biased, halos are used as tracers, we find the opposite effect. The scale at which voids cluster, as well as the void correlation, is similarly sensitive to the sum of neutrino masses and the tracers. This scale dependent trend is not due to simulation volume or halo density. The interplay of these signatures in the void abundance and clustering leaves a distinct fingerprint that could be detected with observations and potentially help break degeneracies between different cosmological parameters. This paper paves the way to exploit cosmic voids in future surveys to constrain the mass of neutrinos.

We estimate the velocity field in a large set of $N$-body simulations including massive neutrino particles, and measure the auto-power spectrum of the velocity divergence field as well as the cross-power spectrum between the cold dark matter density and the velocity divergence. We perform these measurements at four different redshifts and within four different cosmological scenarios, covering a wide range in neutrino masses. We find that the nonlinear correction to the velocity power spectra largely depend on the degree of nonlinear evolution with no specific dependence on the value of neutrino mass. We provide a fitting formula, based on the value of the r.m.s. of the matter fluctuations in spheres of $8h^{-1}$Mpc, describing the nonlinear corrections with 3\% accuracy on scales below $k=0.7\; h$ Mpc$^{-1}$.

Galaxy redshift surveys are one of the pillars of the current standard cosmological model and remain a key tool in the experimental effort to understand the origin of cosmic acceleration. To this end, the next generation of surveys aim at achieving sub-percent precision in the measurement of the equation of state of dark energy $w(z)$ and the growth rate of structure $f(z)$. This however requires comparable control over systematic errors, stressing the need for improved modelling methods. In this contribution we review at the introductory level some highlights of the work done in this direction by the {\it Darklight} project. Supported by an ERC Advanced Grant, {\it Darklight} developed novel techniques for clustering analysis, which were tested through numerical simulations before being finally applied to galaxy data as in particular those of the recently completed VIPERS redshift survey. We focus in particular on: (a) advances on estimating the growth rate of structure from redshift-space distortions; (b) parameter estimation through global Bayesian reconstruction of the density field from survey data; (c) impact of massive neutrinos on large-scale structure measurements. Overall, {\it Darklight} has contributed to paving the way for forthcoming high-precision experiments, such as {\it Euclid}, the next ESA cosmological mission.

We perform a complete study of the gravitational lensing effect beyond the Born approximation on the Cosmic Microwave Background (CMB) anisotropies using a multiple-lens raytracing technique through cosmological N-body simulations of the DEMNUni suite. The impact of second order effects accounting for the non-linear evolution of large-scale structures is evaluated propagating for the first time the full CMB lensing jacobian together with the light rays trajectory. We carefully investigate the robustness of our approach against several numerical effects in the raytracing procedure and in the N-body simulation itself, and find no evidence of large contaminations. We discuss the impact of beyond-Born corrections on lensed CMB observables, and compare our results with recent analytical predictions that appeared in the literature, finding a good agreement, and extend these results to smaller angular scales. We measure the gravitationally-induced CMB polarization rotation that appears at second order, and compare this result with the latest analytical predictions. We then present the detection prospect of beyond-Born effects with the future CMB-S4 experiment. We show that corrections to the temperature power spectrum can be measured only if a good control of the extragalactic foregrounds is achieved. Conversely, the beyond-Born corrections on E and B-modes power spectra will be much more difficult to detect.

The main effect of massive neutrinos on the large-scale structure consists in a few percent suppression of matter perturbations on all scales below their free-streaming scale. Such effect is of particular importance as it allows to constraint the value of the sum of neutrino masses from measurements of the galaxy power spectrum. In this work, we present the first measurements of the next higher-order correlation function, the bispectrum, from N-body simulations that include massive neutrinos as particles. This is the simplest statistics characterising the non-Gaussian properties of the matter and dark matter halos distributions. We investigate, in the first place, the suppression due to massive neutrinos on the matter bispectrum, comparing our measurements with the simplest perturbation theory predictions, finding the approximation of neutrinos contributing at quadratic order in perturbation theory to provide a good fit to the measurements in the simulations. On the other hand, as expected, a linear approximation for neutrino perturbations would lead to O($f_{\nu}$) errors on the total matter bispectrum at large scales. We then attempt an extension of previous results on the universality of linear halo bias in neutrino cosmologies, to non-linear and non-local corrections finding consistent results with the power spectrum analysis.

The clustering ratio is defined as the ratio between the correlation function and the variance of the smoothed overdensity field. In LCDM cosmologies not accounting for massive neutrinos, it has already been proved to be independent from bias and redshift space distortions on a range of linear scales. It therefore allows for a direct comparison of measurements (from galaxies in redshift space) to predictions (for matter in real space). In this paper we first extend the applicability of such properties of the clustering ratio to cosmologies that include massive neutrinos, by performing tests against simulated data. We then investigate the constraining power of the clustering ratio when cosmological parameters such as the total neutrino mass and the equation of state of dark energy are left free. We analyse the joint posterior distribution of the parameters that must satisfy, at the same time, the measurements of the galaxy clustering ratio in the SDSS DR12, and the angular power spectrum of temperature and polarization anisotropies of the CMB measured by the Planck satellite. We find the clustering ratio to be very sensitive to the CDM density parameter, but not very much so to the total neutrino mass. Lastly, we forecast the constraining power the clustering ratio will achieve with forthcoming surveys, predicting the amplitude of its errors in a Euclid-like galaxy survey. In this case, we find it is expected to improve the constraint at 95% level on the CDM density by 40% and on the total neutrino mass by 14%.

The three-point correlation function (3PCF) is a powerful probe to investigate the clustering of matter in the Universe in a complementary way with respect to lower-order statistics, providing additional information with respect to the two-point correlation function and allowing us to shed light on biasing, nonlinear processes, and deviations from Gaussian statistics. In this paper, we analyse the first data release of the VIMOS Public Extragalactic Redshift Survey (VIPERS), determining the dependence of the three-point correlation function on luminosity and stellar mass at $z=[0.5,1.1]$. We exploit the VIPERS Public Data Release 1, consisting of more than 50,000 galaxies with B-band magnitudes in the range $-21.6\lesssim M_{\rm B}-5\log(h)\lesssim-19.9$ and stellar masses in the range $9.8\lesssim\log(M_\star[h^{-2}\,M_\odot])\lesssim 10.7$. We measure both the connected 3PCF and the reduced 3PCF in redshift space, probing different configurations and scales, in the range $2.5<r\,$[Mpc/h]$<20$. We find a significant dependence of the reduced 3PCF on scales and triangle shapes, with stronger anisotropy at larger scales ($r\sim10$ Mpc/h) and an almost flat trend at smaller scales, $r\sim2.5$ Mpc/h. Massive and luminous galaxies present a larger connected 3PCF, while the reduced 3PCF is remarkably insensitive to magnitude and stellar masses in the range we explored. These trends, already observed at low redshifts, are confirmed for the first time to be still valid up to $z=1.1$, providing support to the hierarchical scenario for which massive and bright systems are expected to be more clustered. The possibility of using the measured 3PCF to provide independent constraints on the linear galaxy bias $b$ has also been explored, showing promising results in agreement with other probes.

The set-up of the initial conditions in cosmological N-body simulations is usually implemented by rescaling the desired low-redshift linear power spectrum to the required starting redshift consistently with the Newtonian evolution of the simulation. The implementation of this practical solution requires more care in the context of massive neutrino cosmologies, mainly because of the non-trivial scale-dependence of the linear growth that characterises these models. In this work we consider a simple two-fluid, Newtonian approximation for cold dark matter and massive neutrinos perturbations that can reproduce the cold matter linear evolution predicted by Boltzmann codes such as CAMB or CLASS with a 0.1% accuracy or below for all redshift relevant to nonlinear structure formation. We use this description, in the first place, to quantify the systematic errors induced by several approximations often assumed in numerical simulations, including the typical set-up of the initial conditions for massive neutrino cosmologies adopted in previous works. We then take advantage of the flexibility of this approach to rescale the late-time linear power spectra to the simulation initial redshift, in order to be as consistent as possible with the dynamics of the N-body code and the approximations it assumes. We implement our method in a public code providing the initial displacements and velocities for cold dark matter and neutrino particles that will allow accurate, i.e. one-percent level, numerical simulations for this cosmological scenario.

We present, for the first time in the literature, a full reconstruction of the total (linear and non-linear) ISW/Rees-Sciama effect in the presence of massive neutrinos, together with its cross-correlations with CMB-lensing and weak-lensing signals. The present analyses make use of all-sky maps extracted via ray-tracing across the gravitational potential distribution provided by the "Dark Energy and Massive Neutrino Universe" (DEMNUni) project, a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as separate collisionless particles. We correctly recover, at $1-2\%$ accuracy, the linear predictions from CAMB. Concerning the CMB-lensing and weak-lensing signals, we also recover, with similar accuracy, the signal predicted by Boltzmann codes, once non-linear neutrino corrections to Halofit are accounted for. Interestingly, in the ISW/Rees-Sciama signal, and its cross correlation with lensing, we find an excess of power with respect to the massless case, due to free streaming neutrinos, roughly at the transition scale between the linear and non-linear regimes. The excess is $\sim 5-10\%$ at $l\sim 100$ for the ISW/Rees-Sciama auto power spectrum, depending on the total neutrino mass $M_\nu$, and becomes a factor of $\sim 4$ for $M_\nu=0.3$ eV, at $l\sim 600$, for the ISW/Rees-Sciama cross power with CMB-lensing. This effect should be taken into account for the correct estimation of the CMB temperature bispectrum in the presence of massive neutrinos.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

(abridged) We analyse the clustering features of Large Scale Structures (LSS) in the presence of massive neutrinos, employing a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as a separate collisionless fluid. The volume of 8$\cGpc$, combined with a resolution of about $8\times 10^{10}\Ms$ for the cold dark matter (CDM) component, represents a significant improvement over previous N-body simulations in massive neutrino cosmologies. We show that most of the nonlinear evolution is generated exclusively by the CDM component. We find that accounting only for the nonlinear evolution of the CDM power spectrum allows to recover the total matter power spectrum with the same accuracy as the massless case. Indeed, we show that, the most recent version of the \halofit\ formula calibrated on $\Lambda$CDM simulations can be applied directly to the linear CDM power spectrum without requiring additional fitting parameters in the massive case. As a second step, we study the abundance and clustering properties of CDM halos, confirming that, in massive neutrino cosmologies, the proper definition of the halo bias should be made with respect to the {\em cold} rather than the {\em total} matter distribution, as recently shown in the literature. Here we extend these results to the redshift space, finding that, when accounting for massive neutrinos, an improper definition of the linear bias can lead to a systematic error of about 1-$2 \%$ in the determination of the linear growth rate from anisotropic clustering. This result is quite important if we consider that future spectroscopic galaxy surveys, as \eg\ Euclid, are expected to measure the linear growth-rate with statistical errors less than about $3 \%$ at $z\lesssim1$.

We provide constraints on the accuracy with which the neutrino mass fraction, $f_{\nu}$, can be estimated when exploiting measurements of redshift-space distortions, describing in particular how the error on neutrino mass depends on three fundamental parameters of a characteristic galaxy redshift survey: density, halo bias and volume. In doing this, we make use of a series of dark matter halo catalogues extracted from the BASICC simulation. The mock data are analysed via a Markov Chain Monte Carlo likelihood analysis. We find a fitting function that well describes the dependence of the error on bias, density and volume, showing a decrease in the error as the bias and volume increase, and a decrease with density down to an almost constant value for high density values. This fitting formula allows us to produce forecasts on the precision achievable with future surveys on measurements of the neutrino mass fraction. For example, a Euclid-like spectroscopic survey should be able to measure the neutrino mass fraction with an accuracy of $\delta f_{\nu} \approx 6.7\times10^{-4}$, using redshift-space clustering once all the other cosmological parameters are kept fixed to the $\Lambda$CDM case.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We use the "Dark Energy and Massive Neutrino Universe" (DEMNUni) simulations to compare the constraining power of "sufficient statistics" with the standard matter power spectrum on the sum of neutrino masses, $M_\nu \equiv \sum m_\nu$. In general, the power spectrum, even supplemented with higher moments of the distribution, captures only a fraction of the available cosmological information due to correlations between the Fourier modes. In contrast, the non-linear transform of sufficient statistics, approximated by a logarithmic mapping A=ln(1+\delta), was designed to capture all the available cosmological information contained in the matter clustering; in this sense it is an optimal observable. Our analysis takes advantage of the recent analytical model developed by Carron et al. 2014 to estimate both the matter power spectrum and the A-power spectrum covariance matrices. Using a Fisher information approach, we find that using sufficient statistics increases up to 8 times the available information on the total neutrino mass at z=0, thus tightening the constraints by almost a factor of 3 compared to the matter power spectrum.

We study the gravitational lensing effect on the Cosmic Microwave Background (CMB) anisotropies performing a ray-tracing of the primordial CMB photons through intervening large-scale structures (LSS) distribution predicted by N-Body numerical simulations with a particular focus on the precise recovery of the lens-induced polarized counterpart of the source plane. We apply both a multiple plane ray-tracing and an effective deflection approach based on the Born approximation to deflect the CMB photons trajectories through the simulated lightcone. We discuss the results obtained with both these methods together with the impact of LSS non-linear evolution on the CMB temperature and polarization power spectra. We compare our results with semi-analytical approximations implemented in Boltzmann codes like, e.g., CAMB. We show that, with our current N-body setup, the predicted lensing power is recovered with good accuracy in a wide range of multipoles while excess power with respect to semi-analytic prescriptions is observed in the lensing potential on scales $\ell \gtrsim 3000$. We quantify the impact of the numerical effects connected to the resolution in the N-Body simulation together with the resolution and band-limit chosen to synthesise the CMB source plane. We found these quantities to be particularly important for the simulation of B-mode polarization power spectrum.

Future galaxy redshift surveys aim at probing the clustering of the cosmic large-scale structure with unprecedented accuracy, thus complementing cosmic microwave background experiments in the quest to deliver the most precise and accurate picture ever of our Universe. Analyses of such measurements are usually performed within the context of the so-called vanilla LCDM model - the six-parameter phenomenological model which, for instance, emerges from best fits against the recent data obtained by the Planck satellite. Here, we show that such an approach is prone to subtle systematics when the Gaussianity of primordial fluctuations is concerned. In particular, we demonstrate that, if we neglect even a tiny amount of primordial non-Gaussianity - fully consistent with current limits - we shall introduce spurious biases in the reconstruction of cosmological parameters. This is a serious issue that must be properly accounted for in view of accurate (as well as precise) cosmology.

Massive neutrinos are expected to influence the formation of the large-scale structure of the Universe, depending on the value of their total mass, $\Sigma m_\nu$. In particular Planck data indicate that a non-zero $\Sigma m_\nu$ may help to reconcile CMB data with Sunyaev-Zel'dovich (SZ) cluster surveys. In order to study the impact of neutrinos on the SZ and X-ray cluster properties we run a set of six very large cosmological simulations (8$h^{-3}$ Gpc$^3$ comoving volume) that include a massive neutrino particle component: we consider the values of $\Sigma m_\nu$ = (0, 0.17, 0.34) eV in two cosmological scenarios to test possible degeneracies. Using the halo catalogues extracted from their outputs we produce 50 mock light-cones and, assuming suitable scaling relations, we determine how massive neutrinos affect SZ and X-ray cluster counts, the $y$-parameter and its power spectrum. We provide forecasts for the South Pole Telescope (SPT) and eROSITA cluster surveys, showing that the number of expected detections is reduced by 40 per cent when assuming $\Sigma m_\nu$ =0.34 eV with respect to a model with massless neutrinos. However the degeneracy with $\sigma_8$ and $\Omega_m$ is strong, in particular for X-ray data, requiring the use of additional probes to break it. The $y$-parameter properties are also highly influenced by the neutrino mass fraction, $f_\nu$, with $<\,y\,>\propto(1-f_\nu)^{20}$, considering the cluster component only, and the normalization of the SZ power spectrum is proportional to $(1-f_\nu)^{25-30}$. Comparing our findings with SPT and Atacama Cosmology Telescope measurements at $\ell$ = 3000 indicates that, when Planck cosmological parameters are assumed, a value of $\Sigma m_\nu\simeq0.34$ eV is required to fit with the data.

We reconstruct shear maps and angular power spectra from simulated weakly lensed total intensity (TT) and polarised (EB) maps of the Cosmic Microwave Background (CMB) anisotropies, obtained using Born approximated ray-tracing through the N-body simulated Cold Dark Matter (CDM) structures in the Millennium Simulations (MS). We compare the recovered signal with the $\Lambda$CDM prediction, on the whole interval of angular scales which is allowed by the finite box size, extending from the degree scale to the arcminute, by applying a quadratic estimator in the flat sky limit; we consider PRISM-like instrumental specification for future generation CMB satellites, corresponding to arcminute angular resolution of 3.2' and sensitivity of 2.43 ${\mu}$K-arcmin. -arcmin. The noise bias in the simulations closely follows the estimator prediction, becoming dominated by limits in the angular resolution for the EB signal, at l${\simeq}$1500. The de-biased signal shows no visible departure from predictions of the weak lensing power within uncertainties, when considering TT and EB data singularly. In particular, the reconstruction precision reaches the level of a few percent in bins with $\Delta$l ${\simeq}$100 in the angular multiple interval 1000<l<2000 for TT, and about 10$%$ for EB. Within the adopted specifications, polarisation data do represent a significant contribution to the lensing shear, which appear to faithfully trace the underlying N-body structure down to the smallest angular scales achievable with the present setup, validating at the same time the latter with respect to semi-analytical predictions from ${\Lambda}$CDM cosmology at the level of CMB lensing statistics. This work demonstrates the feasibility of CMB lensing studies based on large scale simulations of cosmological structure formation in the context of the current and future high resolution and sensitivity CMB experiment.

We produce lensing potential and deflection-angle maps in order to simulate CMB weak-lensing via ray-tracing through the COupled Dark Energy Cosmological Simulations (CoDECS). The constructed maps reflect the N-body cosmic structures on a range of scales going from the arcminute to the degree scale. We investigate the variation of the lensing pattern due to the DE dynamics, characterised by different background and perturbation behaviours as a consequence of the interaction between the DE field and CDM. We study the results from three models differing in the background and perturbations evolution, with the purpose to isolate their imprints in the lensing observables. The scenarios investigated include a reference LCDM cosmology, a standard coupled DE (cDE) scenario, and a "bouncing" cDE scenario. For the standard cDE scenario, we find that differences in the lensing potential result from two effects: the enhanced growth of linear CDM density fluctuations with respect to the LCDM case, and the modified nonlinear dynamics of collapsed structures induced by the DE-CDM interaction. As a consequence, CMB lensing highlights the DE impact in the cosmological expansion, even in the degenerate case where the amplitude of the linear matter density perturbations, parametrised through sigma_8, is the same in both the standard cDE and LCDM cosmologies. For the bouncing scenario, we find that the two opposite behaviours of the lens density contrast and of the matter abundance lead to a counter-intuitive effect, making the power of the lensing signal lower by 10% than in the LCDM scenario. Moreover, we compare the behaviour of CDM and baryons in CoDECS separately, in order to isolate effects coming from the coupling with the DE component. We find that, in the bouncing scenario, baryons show an opposite trend with respect to CDM, due to the coupling of the latter with the DE component. [abridged]

We discuss the ability of the planned Euclid mission to detect deviations from General Relativity using its extensive redshift survey of more than 50 Million galaxies. Constraints on the gravity theory are placed measuring the growth rate of structure within 14 redshift bins between z=0.7 and z=2. The growth rate is measured from redshift-space distortions, i.e. the anisotropy of the clustering pattern induced by coherent peculiar motions. This is performed in the overall context of the Euclid spectroscopic survey, which will simultaneously measure the expansion history of the universe, using the power spectrum and its baryonic features as a standard ruler, accounting for the relative degeneracies of expansion and growth parameters. The resulting expected errors on the growth rate in the different redshift bins, expressed through the quantity f\sigma_8, range between 1.3% and 4.4%. We discuss the optimisation of the survey configuration and investigate the important dependence on the growth parameterisation and the assumed cosmological model. We show how a specific parameterisation could actually drive the design towards artificially restricted regions of the parameter space. Finally, in the framework of the popular "\gamma -parameterisation", we show that the Euclid spectroscopic survey alone will already be able to provide substantial evidence (in Bayesian terms) if the growth index differs from the GR value \gamma=0.55 by at least \sim 0.13. This will combine with the comparable inference power provided by the Euclid weak lensing survey, resulting in Euclid's unique ability to provide a decisive test of modified gravity.

In the very last years, cosmological models where the properties of the dark components of the Universe - dark matter and dark energy - are accounted for by a single "dark fluid" have drawn increasing attention and interest. Amongst many proposals, Unified Dark Matter (UDM) cosmologies are promising candidates as effective theories. In these models, a scalar field with a non-canonical kinetic term in its Lagrangian mimics both the accelerated expansion of the Universe at late times and the clustering properties of the large-scale structure of the cosmos. However, UDM models also present peculiar behaviours, the most interesting one being the fact that the perturbations in the dark-matter component of the scalar field do have a non-negligible speed of sound. This gives rise to an effective Jeans scale for the Newtonian potential, below which the dark fluid does not cluster any more. This implies a growth of structures fairly different from that of the concordance LCDM model. In this paper, we demonstrate that forthcoming large-scale surveys will be able to discriminate between viable UDM models and LCDM to a good degree of accuracy. To this purpose, the planned Euclid satellite will be a powerful tool, since it will provide very accurate data on galaxy clustering and the weak lensing effect of cosmic shear. Finally, we also exploit the constraining power of the ongoing CMB Planck experiment. Although our approach is the most conservative, with the inclusion of only well-understood, linear dynamics, in the end we also show what could be done if some amount of non-linear information were included.

We present forecast errors on a wide range of cosmological parameters obtained from a photometric cluster catalogue of a future wide-field Euclid-like survey. We focus in particular on the total neutrino mass as constrained by a combination of the galaxy cluster number counts and correlation function. For the latter we consider only the shape information and the Baryon Acoustic Oscillations (BAO), while marginalising over the spectral amplitude and the redshift space distortions. In addition to the cosmological parameters of the standard LCDM+nu model we also consider a non-vanishing curvature, and two parameters describing a redshift evolution for the dark energy equation of state. For completeness, we also marginalise over a set of "nuisance" parameters, representing the uncertainties on the cluster mass determination. We find that combining cluster counts with power spectrum information greatly improves the constraining power of each probe taken individually, with errors on cosmological parameters being reduced by up to an order of magnitude. In particular, the best improvements are for the parameters defining the dynamical evolution of dark energy, where cluster counts break degeneracies. Moreover, the resulting error on neutrino mass is at the level of \sigma(M_\nu)\sim 0.9 eV, comparable with that derived from present Ly-alpha forest measurements and Cosmic Microwave background (CMB) data in the framework of a non-flat Universe. Further adopting Planck priors and reducing the number of free parameters to a LCDM+nu cosmology allows to place constraints on the total neutrino mass of \sigma(M_\nu) \sim 0.08 eV, close to the lower bound enforced by neutrino oscillation experiments. [abridged]

Cosmological neutrinos strongly affect the evolution of the largest structures in the Universe, i.e. galaxies and galaxy clusters. We use large box-size full hydrodynamic simulations to investigate the non-linear effects that massive neutrinos have on the spatial properties of cold dark matter (CDM) haloes. We quantify the difference with respect to the concordance LambdaCDM model of the halo mass function and of the halo two-point correlation function. We model the redshift-space distortions and compute the errors on the linear distortion parameter beta introduced if cosmological neutrinos are assumed to be massless. We find that, if not taken correctly into account and depending on the total neutrino mass, these effects could lead to a potentially fake signature of modified gravity. Future nearly all-sky spectroscopic galaxy surveys will be able to constrain the neutrino mass if it is larger than 0.6 eV, using beta measurements alone and independently of the value of the matter power spectrum normalisation. In combination with other cosmological probes, this will strengthen neutrino mass constraints and help breaking parameter degeneracies.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We consider cosmological parameters estimation in the presence of a non-zero isocurvature contribution in the primordial perturbations. A previous analysis showed that even a tiny amount of isocurvature perturbation, if not accounted for, could affect standard rulers calibration from Cosmic Microwave Background observations such as those provided by the Planck mission, affect Baryon Acoustic Oscillations interpretation and introduce biases in the recovered dark energy properties that are larger than forecasted statistical errors from future surveys. Extending on this work, here we adopt a general fiducial cosmology which includes a varying dark energy equation of state parameter and curvature. Beside Baryon Acoustic Oscillations measurements, we include the information from the shape of the galaxy power spectrum and consider a joint analysis of a Planck-like Cosmic Microwave Background probe and a future, space-based, Large Scale Structure probe not too dissimilar from recently proposed surveys. We find that this allows one to break the degeneracies that affect the Cosmic Microwave Background and Baryon Acoustic Oscillations combination. As a result, most of the cosmological parameter systematic biases arising from an incorrect assumption on the isocurvature fraction parameter $f_{iso}$, become negligible with respect to the statistical errors. We find that the Cosmic Microwave Background and Large Scale Structure combination gives a statistical error $\sigma(f_{iso}) \sim 0.008$, even when curvature and a varying dark energy equation of state are included, which is smaller than the error obtained from Cosmic Microwave Background alone when flatness and cosmological constant are assumed. These results confirm the synergy and complementarity between Cosmic Microwave Background and Large Scale Structure, and the great potential of future and planned galaxy surveys.

We examine whether future, nearly all-sky galaxy redshift surveys, in combination with CMB priors, will be able to detect the signature of the cosmic neutrino background and determine the absolute neutrino mass scale. We also consider what constraints can be imposed on the effective number of neutrino species. In particular we consider two spectroscopic strategies in the near-IR, the so-called "slitless" and "multi-slit" approaches, whose examples are given by future space-based galaxy surveys, as EUCLID for the slitless case, or SPACE, JEDI, and possibly WFIRST in the future, for the multi-slit case. We find that, in combination with Planck, these galaxy probes will be able to detect at better than 3--sigma level and measure the mass of cosmic neutrinos: a) in a cosmology-independent way, if the sum of neutrino masses is above 0.1 eV; b) assuming spatial flatness and that dark energy is a cosmological constant, otherwise. We find that the sensitivity of such surveys is well suited to span the entire range of neutrino masses allowed by neutrino oscillation experiments, and to yield a clear detection of non-zero neutrino mass. The detection of the cosmic relic neutrino background with cosmological experiments will be a spectacular confirmation of our model for the early Universe and a window into one of the oldest relic components of our Universe.

We investigate the constraints on primordial non-Gaussianity with varied bispectrum shapes that can be derived from the power spectrum of galaxies and clusters of galaxies detected in future wide field optical/near-infrared surveys. Having in mind the proposed ESA space mission \emph{Euclid} as a specific example, we combine the spatial distribution of spectroscopically selected galaxies with that of weak lensing selected clusters. We use the physically motivated halo model in order to represent the correlation function of arbitrary tracers of the Large Scale Structure in the Universe. As naively expected, we find that galaxies are much more effective in jointly constrain the level of primordial non-Gaussianity $f_\mathrm{NL}$ and the amplitude of the matter power spectrum $\sigma_8$ than clusters of galaxies, due to the much lower abundance of the latter that is not adequately compensated by the larger effect on the power spectrum. Nevertheless, combination of the galaxy power spectrum with the cluster-galaxy cross spectrum can decrease the error on the determination of $f_\mathrm{NL}$ by up to a factor of $\sim 2$. This decrement is particularly evident for the less studied non-Gaussian bispectrum shapes, the so-called enfolded and the orthogonal ones. Setting constraints on these models can shed new light on various aspects of the physics of the early Universe, and it is hence of extreme importance. By combining the power spectra of clusters and galaxies with the cluster-galaxy cross spectrum we find constraints on primordial non-Gaussianity of the order $\Delta f_\mathrm{NL} \sim $ a few, competitive and possibly superior to future CMB experiments.

We study the impact of the cosmological parameters uncertainties on the measurements of primordial non-Gaussianity through the large-scale non-Gaussian halo bias effect. While this is not expected to be an issue for the standard LCDM model, it may not be the case for more general models that modify the large-scale shape of the power spectrum. We consider the so-called local non-Gaussianity model and forecasts from planned surveys, alone and combined with a Planck CMB prior. In particular, we consider EUCLID- and LSST-like surveys and forecast the correlations among $f_{\rm NL}$ and the running of the spectral index $\alpha_s$, the dark energy equation of state $w$, the effective sound speed of dark energy perturbations $c^2_s$, the total mass of massive neutrinos $M_\nu=\sum m_\nu$, and the number of extra relativistic degrees of freedom $N_\nu^{rel}$. Neglecting CMB information on $f_{\rm NL}$ and scales $k > 0.03 h$/Mpc, we find that, if $N_\nu^{\rm rel}$ is assumed to be known, the uncertainty on cosmological parameters increases the error on $f_{\rm NL}$ by 10 to 30% depending on the survey. Thus the $f_{\rm NL}$ constraint is remarkable robust to cosmological model uncertainties. On the other hand, if $N_\nu^{\rm rel}$ is simultaneously constrained from the data, the $f_{\rm NL}$ error increases by $\sim 80%$. Finally, future surveys which provide a large sample of galaxies or galaxy clusters over a volume comparable to the Hubble volume can measure primordial non-Gaussianity of the local form with a marginalized 1--$\sigma$ error of the order $\Delta f_{\rm NL} \sim 2-5$, after combination with CMB priors for the remaining cosmological parameters. These results are competitive with CMB bispectrum constraints achievable with an ideal CMB experiment.

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

The description of the abundance and clustering of halos for non-Gaussian initial conditions has recently received renewed interest, motivated by the forthcoming large galaxy and cluster surveys, which can potentially yield constraints of order unity on the non-Gaussianity parameter f_{NL}. We present tests on N-body simulations of analytical formulae describing the halo abundance and clustering for non-Gaussian initial conditions. We calibrate the analytic non-Gaussian mass function of Matarrese et al.(2000) and LoVerde et al.(2008) and the analytic description of clustering of halos for non-Gaussian initial conditions on N-body simulations. We find excellent agreement between the simulations and the analytic predictions if we make the corrections delta_c --> delta_c X sqrt{q} and delta_c --> \delta_c X q where q ~ 0.75, in the density threshold for gravitational collapse and in the non-Gaussian fractional correction to the halo bias, respectively. We discuss the implications of this correction on present and forecasted primordial non-Gaussianity constraints. We confirm that the non-Gaussian halo bias offers a robust and highly competitive test of primordial non-Gaussianity.

We have constructed the first all-sky CMB temperature and polarization lensed maps based on a high-resolution cosmological N-body simulation, the Millennium Simulation (MS). We have exploited the lensing potential map obtained using a map-making procedure (Carbone et al. 2008) which integrates along the line-of-sight the MS dark matter distribution by stacking and randomizing the simulation boxes up to $z = 127$, and which semi-analytically supplies the large-scale power in the angular lensing potential that is not correctly sampled by the N-body simulation. The lensed sky has been obtained by properly modifying the latest version of the LensPix code (Lewis 2005) to account for the MS structures. We have also produced all-sky lensed maps of the so-called $\psi_E$ and $\psi_B$ potentials, which are directly related to the electric and magnetic types of polarization. The angular power spectra of the simulated lensed temperature and polarization maps agree well with semi-analytic estimates up to $l < 2500$, while on smaller scales we find a slight excess of power which we interpret as being due to non-linear clustering in the MS. We also observe how non-linear lensing power in the polarised CMB is transferred to large angular scales by suitably misaligned modes in the CMB and the lensing potential. This work is relevant in view of the future CMB probes, as a way to analyse the lensed sky and disentangle the contribution from primordial gravitational waves.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the deflection angle, aiming at gravitational lensing of the CMB, with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray-tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to $z\sim 11$, producing continuous all-sky maps with arcminute angular resolution. A randomization scheme avoids repetition of structures along the line of sight and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcminutes, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analyzing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees-Sciama effect.

We forecast constraints on primordial non-Gaussianity achievable from forthcoming surveys by exploiting the scale-dependent halo bias introduced on large scales by non-Gaussian initial conditions. We explore the performance of exploiting both the shape of the galaxy power-spectrum on large scales and the cross-correlation of galaxies with Cosmic Microwave Background maps through the Integrated Sachs Wolfe effect. We find that future surveys can detect primordial non-Gaussianity of the local form with a non-Gaussianity parameter $|f_{\rm NL}|$ of order unity. This is particularly exciting because, while the simplest single-field slow-roll models of inflation predict a primordial $|f_{\rm NL}|\ll 1$, this signal sources extra contributions to the effective $f_{\rm NL}$ of large-scale structures that are expected to be above our predicted detection threshold.

The current knowledge of cosmological structure formation suggests that Cold Dark Matter (CDM) halos possess a nonspherical density profile, implying that cosmic structures can be potential sources of gravitational waves via power transfer from scalar perturbations to tensor metric modes in the nonlinear regime. By means of a previously developed mathematical formalism and a triaxial collapse model, we numerically estimate the stochastic gravitational-wave background generated by CDM halos during the fully nonlinear stage of their evolution. Our results suggest that the energy density associated with this background is comparable to that produced by primordial tensor modes at frequencies $\nu\approx10^{-18}-10^{-17}$Hz if the energy scale of inflation is $V^{1/4}\approx 1-2\times 10^{15}$ GeV, and that these gravitational waves could give rise to several cosmological effects, including secondary CMB anisotropy and polarization.

We study the evolution of cosmological perturbations, using a hybrid approximation scheme which upgrades the weak-field limit of Einstein's field equations to account for post-Newtonian scalar and vector metric perturbations and for leading-order source terms of gravitational waves, while including also the first and second-order perturbative approximations. Our equations, which are derived in the Poisson gauge, provide a unified description of matter inhomogeneities in a Universe filled with a pressureless and irrotational fluid and a cosmological constant, ranging from the linear to the highly non-linear regime. The derived expressions for scalar, vector and tensor modes can have a wide range of cosmological applications, ranging from secondary CMB anisotropy and polarization effects, cosmographic relations in a inhomogeneous Universe, gravitational lensing effects and the stochastic gravitational-wave backgrounds generated by non-linear cosmic structures.