Measuring the sum of the three active neutrino masses, $M_\nu$, is one of the most important challenges in modern cosmology. Massive neutrinos imprint characteristic signatures on several cosmological observables in particular on the large-scale structure of the Universe. In order to maximize the information that can be retrieved from galaxy surveys, accurate theoretical predictions in the non-linear regime are needed. Currently, one way to achieve those predictions is by running cosmological numerical simulations. Unfortunately, producing those simulations requires high computational resources -- seven hundred CPU hours for each neutrino mass case. In this work, we propose a new method, based on a deep learning network (U-Net), to quickly generate simulations with massive neutrinos from standard $\Lambda$CDM simulations without neutrinos. We computed multiple relevant statistical measures of deep-learning generated simulations, and conclude that our method accurately reproduces the 3-dimensional spatial distribution of matter down to non-linear scales: $k < 0.7$ h/Mpc. Finally, our method allows us to generate massive neutrino simulations 10,000 times faster than the traditional methods.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

We study the cross-correlation between the Planck CMB lensing convergence map and the eBOSS quasar overdensity obtained from the Sloan Digital Sky Survey (SDSS) IV, in the redshift range $0.9 < z < 2.2$. We detect the CMB lensing convergence-quasar cross power spectrum at $5.4 \sigma$ significance. The cross power spectrum provides a quasar clustering bias measurement that is expected to be particularly robust against systematic effects. The redshift distribution of the quasar sample has a median redshift $z \approx 1.55$, and an effective redshift about $1.51$. The best fit bias of the quasar sample is $b_q = 2.43 \pm 0.45$, corresponding to a host halo mass of $\log_{10}\left( \frac{M}{h^{-1} M_\odot} \right) = 12.54^{+0.25}_{-0.36}$. This is broadly consistent with the previous literature on quasars with a similar redshift range and selection. Since our constraint on the bias comes from the cross-correlation between quasars and CMB lensing, we expect it to be robust to a wide range of possible systematic effects that may contaminate the auto correlation of quasars. We checked for a number of systematic effects from both CMB lensing and quasar overdensity, and found that all systematics are consistent with null within $2 \sigma$. The data is not sensitive to a possible scale dependence of the bias at present, but we expect that as the number of quasars increases (in future surveys such as DESI), it is likely that strong constraints on the scale dependence of the bias can be obtained.

The gravitational potential well of an $M>10^{13}$ $M_\odot$ galaxy will lead to a gravitational redshift differential of order 1 km/s between the galaxy core and its outskirts. Current surveys of massive galaxies with spatially resolved spectroscopy have reached a size which makes feasible attempts to detect gravitational redshifts within galaxies. We use spectra from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) experiment to attempt a measurement of the averaged stellar redshift profile of large elliptical galaxies. We find that systematic effects (possibly related to charge transfer or wavelength calibration accuracy) make the standard MaNGA data pipeline unsuitable for measuring the relevant sub km/s wavelength shifts. We therefore develop a cross-correlation technique to mitigate these effects, but caution that we are working beyond the design accuracy of the MaNGA experiment. With a sample of $272$ galaxies in halos with $\log (M/M_\odot)>13$, we attempt a measurement of the gravitational redshift profile, achieving $1 \sigma$ errors of size $\sim 0.5$ km/s, but are unable to make a significant detection of the signal. Even without a detection, our measurement can be used to limit the dark matter mass within the half light radius of elliptical galaxies to $1.2$ times the stellar mass, at the 68% confidence level. We also perform weighting scheme tests and split sample tests, and address target selection issues and other relativistic effects, including the transverse Doppler effect and relativistic beaming of stars. Future detections could lead to new constraints on the galaxy mass distribution that are different from kinematic and lensing determinations and open a window on galaxy properties and tests of gravity.

Gravitational redshifts and other relativistic effects are beginning to be studied in the context of galaxy clustering. Distortions consistent with those expected in General Relativity have been measured in galaxy cluster redshift profiles by Wojtak et al. and others and in the the cross-correlation function of galaxy populations by Alam et al. On scales below ~20 Mpc/$h$ simulations have shown that gravitational redshift dominates over other effects. However, this signal is related to the shape and depth of gravitational potentials, and therefore the matter density in galaxies and galaxy clusters that is responsible for them. We investigate the effects of baryonic physics on the gravitational redshift profiles of massive (group and cluster-sized) halos. We compare the profiles of different components in halos taken from the MassiveBlack-II cosmological hydrodynamic simulation and a dark matter-only version of the same simulation. We find that inclusion of baryons, cooling, star formation and feedback significantly alters the relevant inner density profiles. These baryonic effects lead to overall increases in both gravitational redshifts and the transverse relativistic Doppler effects by up to ~50% for group size halos. We show how modified Navarro Frenk White halo profiles can be used to parametrize these differences, and provide relevant halo profile fits.

One of the most powerful cosmological datasets when it comes to constraining neutrino masses is represented by galaxy power spectrum measurements, $P_{gg}(k)$. The constraining power of $P_{gg}(k)$ is however severely limited by uncertainties in the modeling of the scale-dependent galaxy bias $b(k)$. In this Letter we present a new method to constrain $b(k)$ by using the cross-correlation between the Cosmic Microwave Background (CMB) lensing signal and galaxy maps ($C_\ell^{\rm \kappa g}$) using a simple but theoretically well-motivated parametrization for $b(k)$. We apply the method using $C_\ell^{\rm \kappa g}$ measured by cross-correlating Planck lensing maps and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 (DR11) CMASS galaxy sample, and $P_{gg}(k)$ measured from the BOSS DR12 CMASS sample. We detect a non-zero scale-dependence at moderate significance, which suggests that a proper modeling of $b(k)$ is necessary in order to reduce the impact of non-linearities and minimize the corresponding systematics. The accomplished increase in constraining power of $P_{gg}(k)$ is demonstrated by determining a 95% C.L. upper bound on the sum of the three active neutrino masses $M_{\nu}$ of $M_{\nu}<0.19\, {\rm eV}$. This limit represents a significant improvement over previous bounds with comparable datasets. Our method will prove especially powerful and important as future large-scale structure surveys will overlap more significantly with the CMB lensing kernel providing a large cross-correlation signal.

Using some of the latest cosmological datasets publicly available, we derive the strongest bounds in the literature on the sum of the three active neutrino masses, $M_\nu$, within the assumption of a background flat $\Lambda$CDM cosmology. In the most conservative scheme, combining Planck cosmic microwave background (CMB) temperature anisotropies and baryon acoustic oscillations (BAO) data, as well as the up-to-date constraint on the optical depth to reionization ($\tau$), the tightest $95\%$ confidence level (C.L.) upper bound we find is $M_\nu<0.151$~eV. The addition of Planck high-$\ell$ polarization data, which however might still be contaminated by systematics, further tightens the bound to $M_\nu<0.118$~eV. A proper model comparison treatment shows that the two aforementioned combinations disfavor the IH at $\sim 64\%$~C.L. and $\sim 71\%$~C.L. respectively. In addition, we compare the constraining power of measurements of the full-shape galaxy power spectrum versus the BAO signature, from the BOSS survey. Even though the latest BOSS full shape measurements cover a larger volume and benefit from smaller error bars compared to previous similar measurements, the analysis method commonly adopted results in their constraining power still being less powerful than that of the extracted BAO signal. Our work uses only cosmological data; imposing the constraint $M_\nu>0.06\,{\rm eV}$ from oscillations data would raise the quoted upper bounds by ${\cal O}(0.1\sigma)$ and would not affect our conclusions.

In a galaxy redshift survey the objects to be targeted for spectra are selected from a photometrically observed sample. The observed magnitudes and colours of galaxies in this parent sample will be affected by their peculiar velocities, through relativistic Doppler and relativistic beaming effects. In this paper we compute the resulting expected changes in galaxy photometry. The magnitudes of the relativistic effects are a function of redshift, stellar mass, galaxy velocity and velocity direction. We focus on the CMASS sample from the Sloan Digital Sky Survey (SDSS), Baryon Oscillation Spectroscopic Survey (BOSS), which is selected on the basis of colour and magnitude. We find that 0.10\% of the sample ($\sim 585$ galaxies) has been scattered into the targeted region of colour-magnitude space by relativistic effects, and conversely 0.09\% of the sample ($\sim 532$ galaxies) has been scattered out. Observational consequences of these effects include an asymmetry in clustering statistics, which we explore in a companion paper. Here we compute a set of weights which can be used to remove the effect of modulations introduced into the density field inferred from a galaxy sample. We conclude by investigating the possible effects of these relativistic modulation on large scale clustering of the galaxy sample.

Large redshift surveys of galaxies and clusters are providing the first opportunities to search for distortions in the observed pattern of large-scale structure due to such effects as gravitational redshift. We focus on non-linear scales and apply a quasi-Newtonian approach using N-body simulations to predict the small asymmetries in the cross-correlation function of two galaxy different populations. Following recent work by Bonvin et al., Zhao and Peacock and Kaiser on galaxy clusters, we include effects which enter at the same order as gravitational redshift: the transverse Doppler effect, light-cone effects, relativistic beaming, luminosity distance perturbation and wide-angle effects. We find that all these effects cause asymmetries in the cross-correlation functions. Quantifying these asymmetries, we find that the total effect is dominated by the gravitational redshift and luminosity distance perturbation at small and large scales, respectively. By adding additional subresolution modelling of galaxy structure to the large-scale structure information, we find that the signal is significantly increased, indicating that structure on the smallest scales is important and should be included. We report on comparison of our simulation results with measurements from the SDSS/BOSS galaxy redshift survey in a companion paper.

General relativistic effects have long been predicted to subtly influence the observed large-scale structure of the universe. The current generation of galaxy redshift surveys have reached a size where detection of such effects is becoming feasible. In this paper, we report the first detection of the redshift asymmetry from the cross-correlation function of two galaxy populations which is consistent with relativistic effects. The dataset is taken from the Sloan Digital Sky Survey DR12 CMASS galaxy sample, and we detect the asymmetry at the $2.7\sigma$ level by applying a shell-averaged estimator to the cross-correlation function. Our measurement dominates at scales around $10$ h$^{-1}$Mpc, larger than those over which the gravitational redshift profile has been recently measured in galaxy clusters, but smaller than scales for which linear perturbation theory is likely to be accurate. The detection significance varies by 0.5$\sigma$ with the details of our measurement and tests for systematic effects. We have also devised two null tests to check for various survey systematics and show that both results are consistent with the null hypothesis. We measure the dipole moment of the cross-correlation function, and from this the asymmetry is also detected, at the $2.8 \sigma$ level. The amplitude and scale-dependence of the clustering asymmetries are approximately consistent with the expectations of General Relativity and a biased galaxy population, within large uncertainties. We explore theoretical predictions using numerical simulations in a companion paper.

We study the asymmetry in the two-point cross-correlation function of two populations of galaxies focusing in particular on the relativistic effects that include the gravitational redshift. We derive the cross-correlation function on small and large scales using two different approaches: General Relativistic and Newtonian perturbation theory. Following recent work by Bonvin et al., Gaztanaga et al. and Croft, we calculate the dipole and the shell estimator with the two procedures and we compare our results. We find that while General Relativistic Perturbation Theory (GRPT) is able to make predictions of relativistic effects on very large, obviously linear scales (r > 50 Mpc/h), the presence of non-linearities physically occurring on much smaller scales (down to those describing galactic potential wells) can strongly affect the asymmetry estimators. These can lead to cancellations of the relativistic terms, and sign changes in the estimators on scales up to r ~ 50 Mpc/h. On the other hand, with an appropriate non-linear gravitational potential, the results obtained using Newtonian theory can successfully describe the asymmetry on smaller, non-linear scales (r < 20 Mpc/h) where gravitational redshift is the dominant term. On larger scales the asymmetry is much smaller in magnitude, and measurement is not within reach of current observations. This is in agreement with the observational results obtained by Gaztnaga et al. and the first detection of relativistic effects (on (r < 20 Mpc/h) scales) by Alam et al.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

In the preprint arxiv:1703.03425 "strong evidence" for the normal neutrino mass ordering is claimed. The authors obtain Bayesian odds of 42:1 in favour of the normal ordering. Their conclusion is based on adopting a flat logarithmic prior for the three neutrino masses. Such an assumption favours a hierarchical spectrum for the masses, which is much easier to accommodate for the normal mass ordering, and hence their prior assumption makes the inverted ordering much less likely a priori. We argue that the claimed "evidence" for normal ordering is almost entirely driven by the adopted prior and not due to the data itself.

We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the $n_s/r$ plane. We study the following neutrino properties: (i) the total neutrino mass $ M_\nu =\sum_i m_i$; (ii) the number of relativistic degrees of freedom $N_{eff}$; and (iii) the neutrino hierarchy: whereas previous literature assumed 3 degenerate neutrino masses or two massless neutrino species (that do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce $< 1 \sigma$ shift of the probability contours in the $n_s/r$ plane with both current or upcoming data. We find that the choice of neutrino hierarchy has a negligible impact. However, the minimal cutoff on the total neutrino mass $M_{\nu,{min}}=0 $ that accompanies previous works using the degenerate hierarchy does introduce biases in the $n_s/r$ plane and should be replaced by $M_{\nu,min}= 0.059$ eV as required by oscillation data. Using current CMB data from Planck and Bicep/Keck, marginalizing over $ M_\nu$ and over $r$ can lead to a shift in the mean value of $n_s$ of $\sim0.3\sigma$ towards lower values. However, once BAO measurements are included, the standard contours in the $n_s/r$ plane are basically reproduced. Larger shifts of the contours in the $n_s/r$ plane (up to 0.8$\sigma$) arise for nonstandard values of $N_{eff}$. We also provide forecasts for the future CMB experiments COrE and Stage-IV and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over $M_\nu$, could induce a shift of $\sim0.4\sigma$ towards lower values in the determination of $n_s$ (or a $\sim 0.8\sigma$ shift if one marginalizes over $N_{eff}$). Comparison to specific inflationary models is shown.

Recent Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the Planck mission have significantly improved previous constraints on the neutrino masses as well as the bounds on extended models with massless or massive sterile neutrino states. However, due to parameter degeneracies, additional low redshift priors are mandatory in order to sharpen the CMB neutrino bounds. We explore here the role of different priors on low redshift quantities, such as the Hubble constant, the cluster mass bias, and the reionization optical depth $\tau$. Concerning current priors on the Hubble constant and the cluster mass bias, the bounds on the neutrino parameters may differ appreciably depending on the choices adopted in the analyses. With regard to future improvements in the priors on the reionization optical depth, a value of $\tau=0.05\pm 0.01$, motivated by astrophysical estimates of the reionization redshift, would lead to $\sum m_\nu<0.0926$~eV at $90\%$~CL, when combining the full \textit{Planck} measurements, Baryon Acoustic Oscillation and Planck clusters data, thereby opening the window to unravel the neutrino mass hierarchy with existing cosmological probes.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

The simplest inflationary models predict a primordial power spectrum (PPS) of the curvature fluctuations that can be described by a power-law function that is nearly scale-invariant. It has been shown, however, that the low-multipole spectrum of the CMB anisotropies may hint the presence of some features in the shape of the scalar PPS, which could deviate from its canonical power-law form. We study the possible degeneracies of this non-standard PPS with the neutrino anisotropies, the neutrino masses, the effective number of relativistic species and a sterile neutrino or a thermal axion mass. The limits on these additional parameters are less constraining in a model with a non-standard PPS when only including the temperature auto-correlation spectrum measurements in the data analyses. The inclusion of the polarization spectra noticeably helps in reducing the degeneracies, leading to results that typically show no deviation from the $\Lambda$CDM model with a standard power-law PPS.

The most recent measurements of the temperature and low-multipole polarization anisotropies of the Cosmic Microwave Background (CMB) from the Planck satellite, when combined with galaxy clustering data from the Baryon Oscillation Spectroscopic Survey (BOSS) in the form of the full shape of the power spectrum, and with Baryon Acoustic Oscillation measurements, provide a $95\%$ confidence level (CL) upper bound on the sum of the three active neutrinos $\sum m _\nu< 0.183$ eV, among the tightest neutrino mass bounds in the literature, to date, when the same datasets are taken into account. This very same data combination is able to set, at $\sim70\%$ CL, an upper limit on $\sum m _\nu$ of $0.0968$ eV, a value that approximately corresponds to the minimal mass expected in the inverted neutrino mass hierarchy scenario. If high-multipole polarization data from Planck is also considered, the $95\%$ CL upper bound is tightened to $\sum m _\nu< 0.176$ eV. Further improvements are obtained by considering recent measurements of the Hubble parameter. These limits are obtained assuming a specific non-degenerate neutrino mass spectrum; they slightly worsen when considering other degenerate neutrino mass schemes. Current cosmological data, therefore, start to be mildly sensitive to the neutrino mass ordering. Low-redshift quantities, such as the Hubble constant or the reionization optical depth, play a very important role when setting the neutrino mass constraints. We also comment on the eventual shifts in the cosmological bounds on $\sum m_\nu$ when possible variations in the former two quantities are addressed.

Recent cosmic microwave background data in temperature and polarization have reached high precision in estimating all the parameters that describe the current so-called standard cosmological model. Recent results about the integrated Sachs-Wolfe effect from cosmic microwave background anisotropies, galaxy surveys, and their cross-correlations are presented. Looking at fine signatures in the cosmic microwave background, such as the lack of power at low multipoles, the primordial power spectrum and the bounds on non-Gaussianities, complemented by galaxy surveys, we discuss inflationary physics and the generation of primordial perturbations in the early Universe. Three important topics in particle physics, the bounds on neutrinos masses and parameters, on thermal axion mass and on the neutron lifetime derived from cosmological data are reviewed, with attention to the comparison with laboratory experiment results. Recent results from cosmic polarization rotation analyses aimed at testing the Einstein equivalence principle are presented. Finally, we discuss the perspectives of next radio facilities for the improvement of the analysis of future cosmic microwave background spectral distortion experiments.

We present new, tight, constraints on the cosmological background of gravitational waves (GWs) using the latest measurements of CMB temperature and polarization anisotropies provided by the Planck, BICEP2 and Keck Array experiments. These constraints are further improved when the GW contribution $N^{\rm GW}_{\rm eff}$ to the effective number of relativistic degrees of freedom $N_{\rm eff}$ is also considered. Parametrizing the tensor spectrum as a power law with tensor-to-scalar ratio $r$, tilt $n_\mathrm{t}$ and pivot $0.01\,\mathrm{Mpc}^{-1}$, and assuming a minimum value of $r=0.001$, we find $r < 0.089$, $n_\mathrm{t} = 1.7^{+2.1}_{-2.0}$ ($95\%\,\mathrm{CL}$, no $N^{\rm GW}_{\rm eff}$) and $r < 0.082$, $n_\mathrm{t} = -0.05^{+0.58}_{-0.87}$ ($95\%\,\mathrm{CL}$, with $N^{\rm GW}_{\rm eff}$). When the recently released $95\,\mathrm{GHz}$ data from Keck Array are added to the analysis, the constraints on $r$ are improved to $r < 0.067$ ($95\%\,\mathrm{CL}$, no $N^{\rm GW}_{\rm eff}$), $r < 0.061$ ($95\%\,\mathrm{CL}$, with $N^{\rm GW}_{\rm eff}$). We discuss the limits coming from direct detection experiments such as LIGO-Virgo, pulsar timing (European Pulsar Timing Array) and CMB spectral distortions (FIRAS). Finally, we show future constraints achievable from a COrE-like mission: if the tensor-to-scalar ratio is of order $10^{-2}$ and the inflationary consistency relation $n_\mathrm{t} = -r/8$ holds, COrE will be able to constrain $n_\mathrm{t}$ with an error of $0.16$ at $95\%\,\mathrm{CL}$. In the case that lensing $B$-modes can be subtracted to $10\%$ of their power, a feasible goal for COrE, these limits will be improved to $0.11$ at $95\%\,\mathrm{CL}$.

The Integrated Sachs-Wolfe (ISW) effect predicts additional anisotropies in the Cosmic Microwave Background due to time variation of the gravitational potential when the expansion of the universe is not matter dominated. The ISW effect is therefore expected in the early universe, due to the presence of relativistic particles at recombination, and in the late universe, when dark energy starts to dominate the expansion. Deviations from the standard picture can be parameterized by $A_{e\text{ISW}}$ and $A_{l\text{ISW}}$, which rescale the overall amplitude of the early and late ISW effects. Analyzing the most recent CMB temperature spectra from the Planck 2015 release, we detect the presence of the early ISW at high significance with $A_{e\mathrm{ISW}} = 1.06\pm0.04$ at 68% CL and an upper limit for the late ISW of $A_{l\mathrm{ISW}} < 1.1$ at 95% CL. The inclusion of the recent polarization data from the Planck experiment results in $A_{e\mathrm{ISW}} = 0.999\pm0.028$ at 68% CL, in better agreement with the value $A_{e\mathrm{ISW}} = 1$ of a standard cosmology. When considering the recent detections of the late ISW coming from correlations between CMB temperature anisotropies and weak lensing, a value of $A_{l\mathrm{ISW}}=0.85\pm0.21$ is predicted at 68% CL, showing a $4\sigma$ evidence. We discuss the stability of our result in the case of an extra relativistic energy component parametrized by the effective neutrino number $N_\mathrm{eff}$ and of a CMB lensing amplitude $A_\mathrm{L}$.

The value of the tensor-to-scalar ratio $r$ in the region allowed by the latest $Planck$ 2015 measurements can be associated to a large variety of inflationary models. We discuss here the potential of future Cosmic Microwave Background cosmological observations in disentangling among the possible theoretical scenarios allowed by our analyses of current $Planck$ temperature and polarization data. Rather than focusing only on $r$, we focus as well on the running of the primordial power spectrum, $\alpha_s$ and the running of thereof, $\beta_s$. Our Fisher matrix method benefits from a detailed and realistic appraisal of the expected foregrounds. Future cosmological probes, as the COrE mission, may be able to reach an unprecedented accuracy in the extraction of $\beta_s$ and rule out the most favoured inflationary models.

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.

Axions currently provide the most compelling solution to the strong CP problem. These particles may be copiously produced in the early universe, including via thermal processes. Therefore, relic axions constitute a hot dark matter component and their masses are strongly degenerate with those of the three active neutrinos, as they leave identical signatures in the different cosmological observables. In addition, thermal axions, while still relativistic states, also contribute to the relativistic degrees of freedom, parameterised via $N_{eff}$. We present the cosmological bounds on the relic axion and neutrino masses, exploiting the full Planck mission data, which include polarization measurements. In the mixed hot dark matter scenario explored here, we find the tightest and more robust constraint to date on the sum of the three active neutrino masses, $\sum m_\nu <0.136$ eV at $95\%$ CL, obtained in the well-known linear perturbation regime. The Planck Sunyaev-Zeldovich cluster number count data further tightens this bound, providing a $95\%$ CL upper limit of $\sum m_\nu <0.126$ eV in this very same mixed hot dark matter model, a value which is very close to the expectations in the inverted hierarchical neutrino mass scenario. Using this same combination of data sets we find the most stringent bound to date on the thermal axion mass, $m_a<0.529$ eV at $95\%$ CL.

We present cosmological bounds on the thermal axion mass in an extended cosmological scenario in which the primordial power spectrum of scalar perturbations differs from the usual power-law shape predicted by the simplest inflationary models. The power spectrum is instead modeled by means of a "piecewise cubic Hermite interpolating polynomial" (PCHIP). When using Cosmic Microwave Background measurements combined with other cosmological data sets, the thermal axion mass constraints are degraded only slightly. The addition of the measurements of $\sigma_8$ and $\Omega_m$ from the 2013 Planck cluster catalogue on galaxy number counts relaxes the bounds on the thermal axion mass, mildly favouring a $\sim 1$~eV axion mass, regardless of the model adopted for the primordial power spectrum.However, in general, such a preference disappears if the sum of the three active neutrino masses is also considered as a free parameter in our numerical analyses, due to the strong correlation between the masses of these two hot thermal relics.

We examine the impact of a non-minimal coupling of the inflaton to the Ricci scalar, $\frac12 \xi R\phi^2$, on the inflationary predictions. Such a non-minimal coupling is expected to be present in the inflaton Lagrangian on fairly general grounds. As a case study, we focus on the simplest inflationary model governed by the potential $V\propto \phi^2$, using the latest combined 2015 analysis of Planck and BICEP2/Keck Array. We find that the presence of a coupling $\xi$ is favoured at a significance of $99\%$ CL, assuming that nature has chosen the potential $V\propto \phi^2$ to generate the primordial perturbations and a number of e-foldings $N=60$. Within the context of the same scenario, we find that the value of $\xi$ is different from zero at the $2\sigma$ level. When considering the cross-correlation polarization spectra from BICEP2/Keck Array and Planck, a value of $r=0.038_{-0.030}^{+0.039}$ is predicted in this particular non-minimally coupled scenario. Future cosmological observations may therefore test these values of $r$ and verify or falsify the non-minimally coupled model explored here.

In this work, we analyze two possible alternative and model-independent approaches to describe the inflationary period. The first one assumes a general equation of state during inflation due to Mukhanov, while the second one is based on the slow-roll hierarchy suggested by Hoffman and Turner. We find that, remarkably, the two approaches are equivalent from the observational viewpoint, as they single out the same areas in the parameter space, and agree with the inflationary attractors where successful inflation occurs. Rephrased in terms of the familiar picture of a slowly rolling, canonically normalized scalar field, the resulting inflaton excursions in these two approaches are almost identical. Furthermore, once the Galactic dust polarization data from Planck are included in the numerical fits, inflaton excursions can safely take sub-Planckian values.

We revisit the cosmological viability of the Hu $\&$ Sawicki modified gravity scenario. The impact of such a modification on the different cosmological observables, including gravitational waves, is carefully described. The most recent cosmological data, as well as constraints on the relationship between the clustering parameter $\sigma_8$ and the current matter mass-energy density $\Omega_m$ from cluster number counts and weak lensing tomography, are considered in our numerical calculations. The strongest bound we find is $|f_{R0}| < 3.7 \times 10^{-6}$ at $95\%$~CL. Forthcoming cluster surveys covering 10,000 deg$^2$ in the sky, with galaxy surface densities of $\mathcal{O}(10)$~arcmin$^{-2}$ could improve the precision in the $\sigma_8$-$\Omega_m$ relationship, tightening the above constraint.

We investigate the axion dark matter scenario (ADM), in which axions account for all of the dark matter in the Universe, in light of the most recent cosmological data. In particular, we use the Planck temperature data, complemented by WMAP E-polarization measurements, as well as the recent BICEP2 observations of B-modes. Baryon Acoustic Oscillation data, including those from the Baryon Oscillation Spectroscopic Survey, are also considered in the numerical analyses. We find that, in the minimal ADM scenario, the full dataset implies that the axion mass m_a = 82.2 pm 1.1 {\mu}eV (corresponding to the Peccei-Quinn symmetry being broken at a scale f_a = (7.54 pm 0.10)*10^10 GeV), or m_a = 76.6 pm 2.6 {\mu}eV (f_a = (8.08 pm 0.27)*10^10 GeV) when we allow for a non- standard effective number of relativistic species Neff . We also find a 2{\sigma} preference for Neff > 3.046. The limit on the sum of neutrino masses is Sum m_{\nu} < 0.25eV at 95% CL for Neff = 3.046, or Sum m_{\nu} < 0.47 eV when Neff is a free parameter. Considering extended scenarios where either the dark energy equation-of-state parameter w, the tensor spectral index nt or the running of the scalar index dns/dlnk are allowed to vary does not change significantly the axion mass-energy density constraints. However, in the case of the full dataset exploited here, there is a preference for a non-zero tensor index or scalar running, driven by the different tensor amplitudes implied by the Planck and BICEP2 observations.

We present up to date cosmological bounds on the sum of active neutrino masses as well as on extended cosmological scenarios with additional thermal relics, as thermal axions or sterile neutrino species. Our analyses consider all the current available cosmological data in the beginning of year 2014, including the very recent and most precise Baryon Acoustic Oscillation (BAO) measurements from the Baryon Oscillation Spectroscopic Survey. In the minimal three active neutrino scenario, we find Sum m_nu < 0.22 eV at 95% CL from the combination of CMB, BAO and Hubble Space Telescope measurements of the Hubble constant. A non zero value for the sum of the three active neutrino masses of about 0.3 eV is significantly favoured at more than 3 standard deviations when adding the constraints on sigma_8 and Omega_m from the Planck Cluster catalog on galaxy number counts. This preference for non zero thermal relic masses disappears almost completely in both the thermal axion and massive sterile neutrino schemes. Extra light species contribute to the effective number of relativistic degrees of freedom, parameterised via Neff. We found that when the recent detection of B mode polarization from the BICEP2 experiment is considered, an analysis of the combined CMB data in the framework of LCDM+r models gives Neff=4.00pm0.41, suggesting the presence of an extra relativistic relic at more than 95 % c.l. from CMB-only data.

New measurements of the cosmic microwave background (CMB) by the Planck mission have greatly increased our knowledge about the Universe. Dark radiation, a weakly interacting component of radiation, is one of the important ingredients in our cosmological model which is testable by Planck and other observational probes. At the moment the possible existence of dark radiation is an unsolved question. For instance, the discrepancy between the value of the Hubble constant, H_0, inferred from the Planck data and local measurements of H_0 can to some extent be alleviated by enlarging the minimal LambdaCDM model to include additional relativistic degrees of freedom. From a fundamental physics point of view dark radiation is no less interesting. Indeed, it could well be one of the most accessible windows to physics beyond the standard model. An example of this is that sterile neutrinos, hinted at in terrestrial oscillation experiments, might also be a source of dark radiation, and cosmological observations can therefore be used to test specific particle physics models. Here we review the most recent cosmological results including a complete investigation of the dark radiation sector in order to provide an overview of models that are still compatible with new cosmological observations. Furthermore we update the cosmological constraints on neutrino physics and dark radiation properties focussing on tensions between data sets and degeneracies among parameters that can degrade our information or mimic the existence of extra species.

We present here bounds on neutrino masses from the combination of recent Planck Cosmic Microwave Background measurements and galaxy clustering information from the Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey-III. We use the full shape of either the photometric angular clustering (Data Release 8) or the 3D spectroscopic clustering (Data Release 9) power spectrum in different cosmological scenarios. In the Lambda$CDM scenario, spectroscopic galaxy clustering measurements improve significantly the existing neutrino mass bounds from Planck data. We find sum m_nu< 0.39 eV at 95% confidence level for the combination of the 3D power spectrum with Planck CMB data (with lensing included) and Wilkinson Microwave Anisoptropy Probe 9-year polarization measurements. Therefore, robust neutrino mass constraints can be obtained without the addition of the prior on the Hubble constant from HST. In extended cosmological scenarios with a dark energy fluid or with non flat geometries, galaxy clustering measurements are essential to pin down the neutrino mass bounds, providing in the majority of cases better results than those obtained from the associated measurement of the Baryon Acoustic Oscillation scale only. In the presence of a freely varying (constant) dark energy equation of state, we find sum m_nu<0.49 eV at 95% confidence level for the combination of the 3D power spectrum with Planck CMB data (with lensing included) and Wilkinson Microwave Anisoptropy Probe 9-year polarization measurements. This same data combination in non flat geometries provides the neutrino mass bound sum m_nu<0.35 eV at 95% confidence level.

Recent Cosmic Microwave Background measurements at high multipoles from the South Pole Telescope and from the Atacama Cosmology Telescope seem to disagree in their conclusions for the neutrino and dark radiation properties. In this paper we set new bounds on the dark radiation and neutrino properties in different cosmological scenarios combining the ACT and SPT data with the nine year data release of the Wilkinson Microwave Anisotropy Probe (WMAP9), Baryon Acoustic Oscillation data, Hubble Telescope measurements of the Hubble constant, and Supernovae Ia luminosity distance data. In the standard three massive neutrino case, the two high multipole probes give similar results if Baryon Acoustic Oscillation data are removed from the analyses and Hubble Telescope measurements are also exploited. A similar result is obtained within a standard cosmology with Neff massless neutrinos, although in this case the agreement between these two measurements is also improved when considering simultaneously Baryon Acoustic Oscillation data and Hubble Space Telescope measurements. In the Neff massive neutrino case the two high multipole probes give very different results regardless of the external data sets used in the combined analyses. In the case in which a dark radiation background with unknown clustering properties is also considered, SPT data seem to exclude the standard value for the dark radiation viscosity cvis=1/3 at the 2 sigma CL, finding evidence for massive neutrinos only when combining SPT data with BAO measurements.

An extra dark radiation component can be present in the universe in the form of sterile neutrinos, axions or other very light degrees of freedom which may interact with the dark matter sector. We derive here the cosmological constraints on the dark radiation abundance, on its effective velocity and on its viscosity parameter from current data in dark radiation-dark matter coupled models. The cosmological bounds on the number of extra dark radiation species do not change significantly when considering interacting schemes. We also find that the constraints on the dark radiation effective velocity are degraded by an order of magnitude while the errors on the viscosity parameter are a factor of two larger when considering interacting scenarios. If future Cosmic Microwave Background data are analysed assuming a non interacting model but the dark radiation and the dark matter sectors interact in nature, the reconstructed values for the effective velocity and for the viscosity parameter will be shifted from their standard 1/3 expectation, namely ceff=0.34 (+0.006 -0.003) and cvis=0.29 (+0.002 -0.001) at 95% CL for the future COrE mission data.

Cosmological constraints on the sum of neutrino masses and on the effective number of neutrino species in standard and non-standard scenarios are computed using the most recent available cosmological data. Our cosmological data sets include the measurement of the Baryonic Acoustic Oscillation (BAO) feature in the Data Release 9 CMASS sample of the Baryon Oscillation Spectroscopic Survey (BOSS). We study in detail the different degeneracies among the parameters, as well as the impact of the different data sets used in the analyses. When considering bounds on the sum of the three active neutrino masses, the information in the BAO signal from galaxy clustering measurements is approximately equally powerful as the shape information from the matter power spectrum. The most stringent bound we find is sum m_nu<0.32 eV at 95 % CL. When non-standard neutrino scenarios with neff massless or massive neutrino species are examined, power spectrum shape measurements provide slightly better bounds than the BAO signal only, due to the breaking of parameter degeneracies. Recent BOSS data combined with CMB and Hubble Space Telescope measurements give neff=3.66^{+0.20 +0.73}_{-0.21 -0.69} in the massless neutrino scenario, and similar results are obtained in the massive case. The evidence for extra radiation neff>3 often claimed in the literature therefore remains at the 2 sigma level when considering up-to-date cosmological data sets. Measurements from the Wilkinson Microwave Anisotropy Probe combined with a prior on the Hubble parameter from the Hubble Space Telescope are very powerful in constraining either the sum of the three active neutrino masses or the number of massless neutrino species. If the former two parameters are allowed to freely vary, however, the bounds from the combination of these two cosmological probes get worse by an order of magnitude.

Recent cosmological data have provided evidence for a "dark" relativistic background at high statistical significance. Parameterized in terms of the number of relativistic degrees of freedom Neff, however, the current data seems to indicate a higher value than the one expected in the standard scenario based on three active neutrinos. This dark radiation component can be characterized not only by its abundance but also by its clustering properties, as its effective sound speed and its viscosity parameter. It is therefore crucial to study the correlations among the dark radiation properties and key cosmological parameters, as the dark energy equation of state or the running of the scalar spectral index, with current and future CMB data. We find that dark radiation with viscosity parameters different from their standard values may be misinterpreted as an evolving dark energy component or as a running spectral index in the power spectrum of primordial fluctuations.

We present neutrino mass bounds using 900,000 luminous galaxies with photometric redshifts measured from Sloan Digital Sky Survey III Data Release Eight (SDSS DR8). The galaxies have photometric redshifts between $z = 0.45$ and $z = 0.65$, and cover 10,000 square degrees and thus probe a volume of 3$h^{-3}$Gpc$^3$, enabling tight constraints to be derived on the amount of dark matter in the form of massive neutrinos. A new bound on the sum of neutrino masses $\sum m_\nu < 0.26$ eV, at 95% confidence level (CL), is obtained after combining our sample of galaxies, which we call "CMASS", with WMAP 7 year Cosmic Microwave Background (CMB) data and the most recent measurement of the Hubble parameter from the Hubble Space Telescope (HST). This constraint is obtained with a conservative multipole range choice of $30 < \ell < 200$ in order to minimize non-linearities, and a free bias parameter in each of the four redshift bins. We study the impact of assuming this linear galaxy bias model using mock catalogs, and find that this model causes a small ($\sim 1-1.5 \sigma$) bias in $\Omega_{\rm DM} h^2$. For this reason, we also quote neutrino bounds based on a conservative galaxy bias model containing additional, shot noise-like free parameters. In this conservative case, the bounds are significantly weakened, e.g. $\sum m_\nu < 0.36$ eV (95% confidence level) for WMAP+HST+CMASS ($\ell_{\rm max}=200$). We also study the dependence of the neutrino bound on multipole range ($\ell_{\rm max}=150$ vs $\ell_{\rm max}=200$) and on which combination of data sets is included as a prior. The addition of supernova and/or Baryon Acoustic Oscillation data does not significantly improve the neutrino mass bound once the HST prior is included. [abridged]

Sterile massive neutrinos are a natural extension of the Standard Model of elementary particles. The energy density of the extra sterile massive states affects cosmological measurements in an analogous way to that of active neutrino species. We perform here an analysis of current cosmological data and derive bounds on the masses of the active and the sterile neutrino states as well as on the number of sterile states. The so-called (3+2) models with three sub-eV active massive neutrinos plus two sub-eV massive sterile species is well within the 95% CL allowed regions when considering cosmological data only. If the two extra sterile states have thermal abundances at decoupling, Big Bang Nucleosynthesis bounds compromise the viability of (3+2) models. Forecasts from future cosmological data on the active and sterile neutrino parameters are also presented. Independent measurements of the neutrino mass from tritium beta decay experiments and of the Hubble constant could shed light on sub-eV massive sterile neutrino scenarios.

Cosmological measurements are affected by the energy density of both active and sterile massive neutrinos. We extend here a recent analysis of current cosmological data to non minimal cosmologies. Several possible scenarios are examined: a constant w \neq -1 dark energy equation of state, a non flat universe, a time varying dark energy component and coupled dark matter dark energy universes or modified gravity scenarios. When considering cosmological data only, (3+2) massive neutrino models with ~0.5 eV sterile species are allowed at 95% CL. This scenario has been shown to reconcile reactor, LSND and MiniBooNE positive signals with null results from other searches. Big Bang Nucleosynthesis bounds could compromise the viability of (3+2) models if the two sterile species are fully thermalized states at decoupling.

Determination of whether the Harrison--Zel'dovich spectrum for primordial scalar perturbations is consistent with observations is sensitive to assumptions about the reionization scenario. In light of this result, we revisit constraints on inflationary models using more general reionization scenarios. While the bounds on the tensor-to-scalar ratio are largely unmodified, when different reionization schemes are addressed, hybrid models are back into the inflationary game. In the general reionization picture, we reconstruct both the shape and amplitude of the inflaton potential. We find a broader spectrum of potential shapes when relaxing the simple reionization restriction. An upper limit of $10^{16}$ GeV to the amplitude of the potential is found, regardless of the assumptions on the reionization history.

Inflation predicts primordial scalar perturbations with a nearly scale-invariant spectrum and a spectral index approximately unity (the Harrison--Zel'dovich (HZ) spectrum). The first important step for inflationary cosmology is to check the consistency of the HZ primordial spectrum with current observations. Recent analyses have claimed that a HZ primordial spectrum is excluded at more than 99% c.l.. Here we show that the HZ spectrum is only marginally disfavored if one considers a more general reionization scenario. Data from the Planck mission will settle the issue.