We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. The net effect of the improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is ($3366.6 \pm 2.7$)$\mu$K, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of $\tau = 0.051 \pm 0.006$, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first (and significant) detection of this signal in the tSZ maps, and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

The symmetron is a typical example of screened modified gravity, wherein the symmetron force is dynamically suppressed in dense environments. This allows it to hide in traditional tests of gravity. However, the past decade has seen great experimental progress towards measuring screened forces in the laboratory or in space. Screening relies on nonlinearities in the equation of motion, which significantly complicates the theoretical analysis of such forces. Here, we present a calculation of the symmetron force between a dense plate and sphere surrounded by vacuum. This is done via semi-analytical approaches in two limiting cases, based on the size of the sphere: large spheres are analyzed via the proximity force approximation, whilst small spheres are treated as screened test particles. In the intermediate regime we solve the problem numerically. Our results allow us to make contact with Casimir force experiments, which often employ a plate and sphere configuration for practical reasons, and may therefore be used to constrain symmetrons. We use our results to forecast constraints on the symmetron's parameters for a hypothetical Casimir experiment that is based on the current state of the art. The forecasts compare favorably to other leading laboratory tests of gravity, particularly atom interferometry and bouncing neutrons. We thus conclude that near-future Casimir experiments will be capable of placing tight new bounds on symmetrons. Our results for the symmetron force are derived in a scale-invariant way, such that although we here focus on Casimir experiments, they may be applied to any other plate-sphere system, ranging from microscopic to astrophysical scales.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $\Lambda$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $\Lambda$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales.

We study the first cosmological implications of the mimetic theory of massive gravity recently proposed by Chamseddine and Mukhanov. This is a novel theory of ghost-free massive gravity which additionally contains a mimetic dark matter component. In an echo of other modified gravity theories, there are self-accelerating solutions which contain a ghost instability. In the ghost-free region of parameter space, the effect of the graviton mass on the cosmic expansion history amounts to an effective negative cosmological constant, a radiation component, and a negative curvature term. This allows us to place constraints on the model parameters---the graviton mass and the St\"uckelberg vacuum expectation value---by insisting that the effective radiation and curvature terms be within observational bounds. The late-time acceleration must be accounted for by a separate positive cosmological constant or other dark energy sector. We impose further constraints at the level of perturbations by demanding linear stability. We comment on the possibility of distinguishing this theory from $\Lambda$CDM with current and future large-scale structure surveys.

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged).

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $\beta_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We review the advanced version of the KKLT construction and pure $d=4$ de Sitter supergravity, involving a nilpotent multiplet, with regard to various conjectures that de Sitter state cannot exist in string theory. We explain why we consider these conjectures problematic and not well motivated, and why the recently proposed alternative string theory models of dark energy, ignoring vacuum stabilization, are ruled out by cosmological observations at least at the $3\sigma$ level, i.e. with more than $99.7\%$ confidence.

This paper presents the Planck Multi-frequency Catalogue of Non-thermal (i.e. synchrotron-dominated) Sources (PCNT) observed between 30 and 857 GHz by the ESA Planck mission. This catalogue was constructed by selecting objects detected in the full mission all-sky temperature maps at 30 and 143 GHz, with a signal-to-noise ratio (S/N)>3 in at least one of the two channels after filtering with a particular Mexican hat wavelet. As a result, 29400 source candidates were selected. Then, a multi-frequency analysis was performed using the Matrix Filters methodology at the position of these objects, and flux densities and errors were calculated for all of them in the nine Planck channels. The present catalogue is the first unbiased, full-sky catalogue of synchrotron-dominated sources published at millimetre and submillimetre wavelengths and constitutes a powerful database for statistical studies of non-thermal extragalactic sources, whose emission is dominated by the central active galactic nucleus. Together with the full multi-frequency catalogue, we also define the Bright Planck Multi-frequency Catalogue of Non-thermal Sources PCNTb, where only those objects with a S/N>4 at both 30 and 143 GHz were selected. In this catalogue 1146 compact sources are detected outside the adopted Planck GAL070 mask; thus, these sources constitute a highly reliable sample of extragalactic radio sources. We also flag the high-significance subsample PCNThs, a subset of 151 sources that are detected with S/N>4 in all nine Planck channels, 75 of which are found outside the Planck mask adopted here. The remaining 76 sources inside the Galactic mask are very likely Galactic objects.

We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters.

Using the ${\it Planck}$ full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the ${\it Planck}$ foreground-cleaned `${\tt 2D-ILC\,}$' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized ${\it rms}$ dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of $1526$ clusters is $\langle \left(\Delta T/T \right)^{2} \rangle = (1.64 \pm 0.48) \times 10^{-11}$. However, comparison with analytic calculations and simulations suggest that around $0.7\,\sigma$ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be $\langle \left(\Delta T/T \right)^{2} \rangle = (1.35 \pm 0.48) \times 10^{-11}$, which gives a detection at the $2.8\,\sigma$ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is $\langle v^{2} \rangle =(123\,000 \pm 71\,000)\,({\rm km}\,{\rm s}^{-1})^{2}$, which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of $600\,h^{-1}{\rm Mpc}$. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35\mu$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level.

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects.

We study the implications of the recent detection of gravitational waves emitted by a pair of merging neutron stars and their electromagnetic counterpart, events GW170817 and GRB170817A, on the viability of the doubly coupled bimetric models of cosmic evolution, where the two metrics couple directly to matter through a composite, effective metric. We demonstrate that the bounds on the speed of gravitational waves place strong constraints on the doubly coupled models, forcing either the two metrics to be proportional at the background level or the models to become singly coupled. Proportional backgrounds are particularly interesting as they provide stable cosmological solutions with phenomenologies equivalent to that of $\Lambda$CDM at the background level as well as for linear perturbations, while nonlinearities are expected to show deviations from the standard model.

Over the last few years, a large family of cosmological attractor models has been discovered, which can successfully match the latest inflation-related observational data. Many of these models can also describe a small cosmological constant $\Lambda$, which provides the most natural description of the present stage of the cosmological acceleration. In this paper, we study $\alpha$-attractor models with dynamical dark energy, including the cosmological constant $\Lambda$ as a free parameter. Predominantly, the models with $\Lambda > 0$ converge to the asymptotic regime with the equation of state $w=-1$. However, there are some models with $w\neq -1$, which are compatible with the current observations. In the simplest models with $\Lambda = 0$, one has the tensor to scalar ratio $r=\frac{12\alpha}{N^2}$ and the asymptotic equation of state $w=-1+\frac{2}{9\alpha}$ (which in general differs from its present value). For example, in the seven disk M-theory related model with $\alpha = 7/3$ one finds $r \sim 10^{-2}$ and the asymptotic equation of state is $w \sim -0.9$. Future observations, including large-scale structure surveys as well as B-mode detectors will test these, as well as more general models presented here. We also discuss gravitational reheating in models of quintessential inflation and argue that its investigation may be interesting from the point of view of inflationary cosmology. Such models require a much greater number of $e$-folds, and therefore predict a spectral index $n_{s}$ that can exceed the value in more conventional models by about $0.006$. This suggests a way to distinguish the conventional inflationary models from the models of quintessential inflation, even if they predict $w = -1$.

We propose a simple, nonlocal modification to general relativity (GR) on large scales, which provides a model of late-time cosmic acceleration in the absence of the cosmological constant and with the same number of free parameters as in standard cosmology. The model is motivated by adding to the gravity sector an extra spin-2 field interacting nonlocally with the physical metric coupled to matter. The form of the nonlocal interaction is inspired by the simplest form of the Deser-Woodard (DW) model, $\alpha R\frac{1}{\Box}R$, with one of the Ricci scalars being replaced by a constant $m^{2}$, and gravity is therefore modified in the infrared by adding a simple term of the form $m^2\frac{1}{\Box}R$ to the Einstein-Hilbert term. We study cosmic expansion histories, and demonstrate that the new model can provide background expansions consistent with observations if $m$ is of the order of the Hubble expansion rate today, in contrast to the simple DW model with no viable cosmology. The model is best fit by $w_0\sim-1.075$ and $w_a\sim0.045$. We also compare the cosmology of the model to that of Maggiore and Mancarella (MM), $m^2R\frac{1}{\Box^2}R$, and demonstrate that the viable cosmic histories follow the standard-model evolution more closely compared to the MM model. We further demonstrate that the proposed model possesses the same number of physical degrees of freedom as in GR. Finally, we discuss the appearance of ghosts in the local formulation of the model, and argue that they are unphysical and harmless to the theory, keeping the physical degrees of freedom healthy.

The existence (and abundance) of primordial black holes (PBHs) is governed by the power spectrum of primordial perturbations generated during inflation. So far no PBHs have been observed, and instead, increasingly stringent bounds on their existence at different scales have been obtained. Up until recently, this has been exploited in attempts to constrain parts of the inflationary power spectrum that are unconstrained by cosmological observations. We first point out that the simple translation of the PBH non-observation bounds into constraints on the primordial power spectrum is inaccurate as it fails to include realistic aspects of PBH formation and evolution. We then demonstrate, by studying two examples of uncertainties from the effects of critical and non-spherical collapse, that even though they may seem small, they have important implications for the usefulness of the constraints. In particular, we point out that the uncertainty induced by non-spherical collapse may be much larger than the difference between particular bounds from PBH non-observations and the general maximum cap stemming from the condition $\Omega \leq 1$ on the dark-matter density in the form of PBHs. We therefore make the cautious suggestion of applying only the overall maximum dark-matter constraint to models of early Universe, as this requirement seems to currently provide a more reliable constraint, which better reflects our current lack of detailed knowledge of PBH formation. These, and other effects, such as merging, clustering and accretion, may also loosen constraints from non-observations of other primordial compact objects such as ultra-compact minihalos of dark matter.

The six parameters of the standard $\Lambda$CDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We investigate these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium $\tau$, the baryon density $\omega_{\rm b}$, the matter density $\omega_{\rm m}$, the angular size of the sound horizon $\theta_*$, the spectral index of the primordial power spectrum, $n_{\rm s}$, and $A_{\rm s}e^{-2\tau}$ (where $A_{\rm s}$ is the amplitude of the primordial power spectrum), we examine the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment $\ell<800$ in the Planck temperature power spectrum) and an all angular-scale data set ($\ell<2500$ Planck temperature power spectrum), each with a prior on $\tau$ of $0.07\pm0.02$. We find that the shifts, in units of the 1$\sigma$ expected dispersion for each parameter, are $\{\Delta \tau, \Delta A_{\rm s} e^{-2\tau}, \Delta n_{\rm s}, \Delta \omega_{\rm m}, \Delta \omega_{\rm b}, \Delta \theta_*\} = \{-1.7, -2.2, 1.2, -2.0, 1.1, 0.9\}$, with a $\chi^2$ value of 8.0. We find that this $\chi^2$ value is exceeded in 15% of our simulated data sets, and that a parameter deviates by more than 2.2$\sigma$ in 9% of simulated data sets, meaning that the shifts are not unusually large. Comparing $\ell<800$ instead to $\ell>800$, or splitting at a different multipole, yields similar results. We examine the $\ell<800$ model residuals in the $\ell>800$ power spectrum data and find that the features there... [abridged]

We discuss the cosmological implications of nonlocal modifications of general relativity containing tensorial structures. Assuming the presence of standard radiation- and matter-dominated eras, we show that, except in very particular cases, the nonlocal terms contribute a rapidly growing energy density. These models therefore generically do not have a stable cosmological evolution.

Measurements of flux density are described for five planets, Mars, Jupiter, Saturn, Uranus, and Neptune, across the six Planck High Frequency Instrument frequency bands (100-857 GHz) and these are then compared with models and existing data. In our analysis, we have also included estimates of the brightness of Jupiter and Saturn at the three frequencies of the Planck Low Frequency Instrument (30, 44, and 70 GHz). The results provide constraints on the intrinsic brightness and the brightness time-variability of these planets. The majority of the planet flux density estimates are limited by systematic errors, but still yield better than 1% measurements in many cases. Applying data from Planck HFI, the Wilkinson Microwave Anisotropy Probe (WMAP), and the Atacama Cosmology Telescope (ACT) to a model that incorporates contributions from Saturn's rings to the planet's total flux density suggests a best fit value for the spectral index of Saturn's ring system of $\beta _\mathrm{ring} = 2.30\pm0.03$ over the 30-1000 GHz frequency range. The average ratio between the Planck-HFI measurements and the adopted model predictions for all five planets (excluding Jupiter observations for 353 GHz) is 0.997, 0.997, 1.018, and 1.032 for 100, 143, 217, and 353 GHz, respectively. Model predictions for planet thermodynamic temperatures are therefore consistent with the absolute calibration of Planck-HFI detectors at about the three-percent-level. We compare our measurements with published results from recent cosmic microwave background experiments. In particular, we observe that the flux densities measured by Planck HFI and WMAP agree to within 2%. These results allow experiments operating in the mm-wavelength range to cross-calibrate against Planck and improve models of radiative transport used in planetary science.

Many theories of modified gravity with higher order derivatives are usually ignored because of serious problems that appear due to an additional ghost degree of freedom. Most dangerously, it causes an immediate decay of the vacuum. However, breaking Lorentz invariance can cure such abominable behavior. By analyzing a model that describes a massive graviton together with a remaining Boulware-Deser ghost mode we show that even ghostly theories of modified gravity can yield models that are viable at both classical and quantum levels and, therefore, they should not generally be ruled out. Furthermore, we identify the most dangerous quantum scattering process that has the main impact on the decay time and find differences to simple theories that only describe an ordinary scalar field and a ghost. Additionally, constraints on the parameters of the theory including some upper bounds on the Lorentz-breaking cutoff scale are presented. In particular, for a simple theory of massive gravity we find that a breaking of Lorentz invariance is allowed to happen even at scales above the Planck mass. Finally, we discuss the relevance to other theories of modified gravity.

Theories of massive gravity with one or two dynamical metrics generically lack stable and observationally-viable cosmological solutions that are distinguishable from $\Lambda$CDM. We consider an extension to trimetric gravity, with three interacting spin-2 fields which are not plagued by the Boulware-Deser ghost. We systematically explore every combination with two free parameters in search of background cosmologies that are competitive with $\Lambda$CDM. For each case we determine whether the expansion history satisfies viability criteria, and whether or not it contains beyond-$\Lambda$CDM phenomenology. Among the many models we consider, there are only three cases that seem to be both viable and distinguishable from standard cosmology. One of the models has only one free parameter and displays a crossing from above to below the phantom divide. The other two provide scaling behavior, although they contain future singularities that need to be studied in more detail. These models possess interesting features that make them compelling targets for a full comparison to observations of both cosmological expansion history and structure formation.

We discuss the cosmological implications of the $R~\Box^{-2}R$ nonlocal modification to standard gravity. We relax the assumption of special initial conditions in the local formulation of the theory, perform a full phase-space analysis of the system, and show that the late-time cosmology of the model exhibits two distinct evolution paths, on which a large range of values for the present equation of state can be reached. We then compare the general solutions to supernovae data and place constraints on the parameters of the model. In particular, we find that the mass parameter of the theory should be smaller than 1.2 in Hubble units.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

Massive bigravity models are interesting alternatives to standard cosmology. In most cases, however, these models have been studied for a simplified scenario in which both metrics take homogeneous and isotropic forms [Friedmann-Lema\^itre-Robertson-Walker (FLRW)] with the same spatial curvatures. The interest to consider more general geometries arises, in particular, in view of the difficulty so far encountered in building stable cosmological solutions with homogeneous and isotropic metrics. Here we consider a number of cases in which the two metrics take more general forms, namely FLRW with different spatial curvatures---Lema\^itre, Lema\^itre-Tolman-Bondi (LTB), and Bianchi I---as well as cases where only one metric is linearly perturbed. We discuss possible consistent combinations and find that only some special cases of FLRW--Lema\^itre, LTB--LTB, and FLRW--Bianchi I combinations give consistent, nontrivial solutions.

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.

Bimetric theory describes gravitational interactions in the presence of an extra spin-2 field. Previous work has suggested that its cosmological solutions are generically plagued by instabilities. We show that by taking the Planck mass for the second metric, $M_f$, to be small, these instabilities can be pushed back to unobservably early times. In this limit, the theory approaches general relativity with an effective cosmological constant which is, remarkably, determined by the spin-2 interaction scale. This provides a late-time expansion history which is extremely close to $\Lambda$CDM, but with a technically-natural value for the cosmological constant. We find $M_f$ should be no larger than the electroweak scale in order for cosmological perturbations to be stable by big-bang nucleosynthesis. We further show that in this limit the helicity-0 mode is no longer strongly-coupled at low energy scales.

We study the integrated Sachs-Wolfe (ISW) effect in ghost-free, massive bigravity. We focus on the infinite-branch bigravity (IBB) model which exhibits viable cosmic expansion histories and stable linear perturbations, while the cosmological constant is set to zero and the late-time accelerated expansion of the Universe is due solely to the gravitational interaction terms. The ISW contribution to the CMB auto-correlation power spectrum is predicted, as well as the cross-correlation between the CMB temperature anisotropies and large-scale structure. We use ISW amplitudes as observed by the WMAP 9-year temperature data together with galaxy and AGN data provided by the WISE mission in order to compare the theoretical predictions to the observations. The ISW amplitudes in IBB are found to be larger than the corresponding ones in the standard LCDM model by roughly a factor of 1.5, but are still consistent with the observations.

There is a no-go theorem forbidding flat and closed FLRW solutions in massive gravity on a flat reference metric, while open solutions are unstable. Recently it was shown that this no-go theorem can be overcome if at least some matter couples to a hybrid metric composed of both the dynamical and the fixed reference metric. We show that this is not compatible with the standard description of cosmological sources in terms of effective perfect fluids, and the predictions of the theory become sensitive either to the detailed field-theoretical modelling of the matter content or to the presence of additional dark degrees of freedom. This is a serious practical complication. Furthermore, we demonstrate that viable cosmological background evolution with a perfect fluid appears to require the presence of fields with highly contrived properties. This could be improved if the equivalence principle is broken by coupling only some of the fields to the composite metric, but viable self-accelerating solutions due only to the massive graviton are difficult to obtain.

We study the cosmic expansion history of massive bigravity with a viable matter coupling which treats both metrics on equal footing. We derive the Friedmann equation for the effective metric through which matter couples to the two metrics, and study its solutions. For certain parameter choices, the background cosmology is identical to that of LCDM. More general parameters yield dynamical dark energy, which can still be in agreement with observations of the expansion history. We study specific parameter choices of interest, including minimal models, maximally-symmetric models, and a candidate partially-massless theory.

Nonlinear, ghost-free massive gravity has two tensor fields; when both are dynamical, the mass of the graviton can lead to cosmic acceleration that agrees with background data, even in the absence of a cosmological constant. Here the question of the stability of linear perturbations in this bimetric theory is examined. Instabilities are presented for several classes of models, and simple criteria for the cosmological stability of massive bigravity are derived. In this way, we identify a particular self-accelerating bigravity model, infinite-branch bigravity (IBB), which exhibits both viable background evolution and stable linear perturbations. We discuss the modified gravity parameters for IBB, which do not reduce to the standard $\Lambda$CDM result at early times, and compute the combined likelihood from measured growth data and type Ia supernovae. IBB predicts a present matter density $\Omega_{m0}=0.18$ and an equation of state $w(z)=-0.79+0.21z/(1+z)$. The growth rate of structure is well-approximated at late times by $f(z)\approx\Omega_{m}^{0.47}[1+0.21z/(1+z)]$. The implications of the linear instability for other bigravity models are discussed: the instability does not necessarily rule these models out, but rather presents interesting questions about how to extract observables from them when linear perturbation theory does not hold.

We study the asymmetric primordial fluctuations in a model of inflation in which translational invariance is broken by a domain wall. We calculate the corrections to the power spectrum of curvature perturbations; they are anisotropic and contain dipole, quadrupole, and higher multipoles with non-trivial scale-dependent amplitudes. Inspired by observations of these multipole asymmetries in terms of two-point correlations and variance in real space, we demonstrate that this model can explain the observed anomalous power asymmetry of the cosmic microwave background (CMB) sky, including its characteristic feature that the dipole dominates over higher multipoles. We test the viability of the model and place approximate constraints on its parameters by using observational values of dipole, quadrupole, and octopole amplitudes of the asymmetry measured by a local-variance estimator. We find that a configuration of the model in which the CMB sphere does not intersect the domain wall during inflation provides a good fit to the data. We further derive analytic expressions for the corrections to the CMB temperature covariance matrix, or angular power spectra, which can be used in future statistical analysis of the model in spherical harmonic space.

The ghost-free theory of massive gravity with two dynamical metrics has been shown to produce viable cosmological expansion, where the late-time acceleration of the Universe is due to the finite range of the gravitational interaction rather than a nonzero cosmological constant. Here the cosmological perturbations are studied in this theory. The full perturbation equations are presented in a general gauge and analyzed, focusing on subhorizon scales in the quasistatic limit during the matter-dominated era. An evolution equation for the matter inhomogeneities and the parameters quantifying the deviations from general relativistic structure formation are expressed in terms of five functions whose forms are determined directly by the coupling parameters in the theory. The evolution equation has a similar structure to Horndeski-type scalar-tensor theories, exhibiting a modified growth rate and scale-dependence at intermediate wavenumbers. Predictions of the theory are confronted with observational data on both background expansion and large-scale structure, although care must be taken to ensure a model is stable. It is found that while the stable models fit the data well, they feature deviations from the standard cosmology that could be detected or ruled out by near-future experiments.

The possibility of matter coupling to two metrics at once is considered. This appears natural in the most general ghost-free, bimetric theory of gravity, where it unlocks an additional symmetry with respect to the exchange of the metrics. This double coupling, however, raises the problem of identifying the observables of the theory. It is shown that if the two metrics couple minimally to matter, then there is no physical metric to which all matter would universally couple, and that moreover such an effective metric generically does not exist even for an individual matter species. By studying point particle dynamics, a resolution is suggested in the context of Finsler geometry.

We revisit the question of hemispherical power asymmetry in the WMAP and Planck temperature sky maps by measuring the local variance over the sky and on disks of various sizes. For the 2013 Planck sky map we find that none of the 1000 available isotropic Planck "Full Focal Plane" simulations have a larger variance asymmetry than that estimated from the data, suggesting the presence of an anisotropic signature formally significant at least at the $3.3\sigma$ level. For the WMAP 9-year data we find that 5 out of 1000 simulations have a larger asymmetry. The preferred direction for the asymmetry from the Planck data is $(l,b)=(212^\circ,-13^\circ)$, in good agreement with previous reports of the same hemispherical power asymmetry.

A ghost-free theory of gravity with two dynamical metrics both coupled to matter is shown to be consistent and viable. Its cosmological implications are studied, and the models, in particular in the context of partially massless gravity, are found to explain the cosmic acceleration without resorting to dark energy.

We study the background cosmology of the ghost-free, bimetric theory of gravity. We perform an extensive statistical analysis of the model using both frequentist and Bayesian frameworks and employ the constraints on the expansion history of the Universe from the observations of supernovae, the cosmic microwave background and the large scale structure to estimate the model's parameters and test the goodness of the fits. We explore the parameter space of the model with nested sampling to find the best-fit chi-square, obtain the Bayesian evidence, and compute the marginalized posteriors and mean likelihoods. We mainly focus on a class of sub-models with no explicit cosmological constant (or vacuum energy) term to assess the ability of the theory to dynamically cause a late-time accelerated expansion. The model behaves as standard gravity without a cosmological constant at early times, with an emergent extra contribution to the energy density that converges to a cosmological constant in the far future. The model can in most cases yield very good fits and is in perfect agreement with the data. This is because many points in the parameter space of the model exist that give rise to time-evolution equations that are effectively very similar to those of the $\Lambda$CDM. This similarity makes the model compatible with observations as in the $\Lambda$CDM case, at least at the background level. Even though our results indicate a slightly better fit for the $\Lambda$CDM concordance model in terms of the $p$-value and evidence, none of the models is statistically preferred to the other. However, the parameters of the bigravity model are in general degenerate. A similar but perturbative analysis of the model as well as more data will be required to break the degeneracies and constrain the parameters, in case the model will still be viable compared to the $\Lambda$CDM.

The recent discovery of a ghost-free, non-linear extension of the Fierz-Pauli theory of massive gravity, and its bigravity formulation, introduced new possibilities of interpreting cosmological observations, in particular, the apparent late-time accelerated expansion of the Universe. Here we discuss such possibilities by studying the background cosmology of the model and comparing its predictions to different cosmological measurements. We place constraints on the model parameters through an extensive statistical analysis of the model, and compare its viability to that of the standard model of cosmology. We demonstrate that the model can yield perfect fits to the data and is capable of explaining the cosmic acceleration in the absence of an explicit cosmological constant or dark energy, but there are a few caveats that must be taken into account when interpreting the results.

For a Gaussian spectrum of primordial density fluctuations, ultracompact minihalos (UCMHs) of dark matter are expected to be produced in much greater abundance than, e.g., primordial black holes. Forming shortly after matter-radiation equality, these objects would develop very dense and spiky dark matter profiles. In the standard scenario where dark matter consists of thermally-produced, weakly-interacting massive particles, UCMHs could thus appear as highly luminous gamma-ray sources, or leave an imprint in the cosmic microwave background by changing the reionisation history of the Universe. We derive corresponding limits on the cosmic abundance of UCMHs at different epochs, and translate them into constraints on the primordial power spectrum. We find the resulting constraints to be quite severe, especially at length scales much smaller than what can be directly probed by the cosmic microwave background or large-scale structure observations. We use our results to provide an updated compilation of the best available constraints on the power of density fluctuations on all scales, ranging from the present-day horizon to scales more than 20 orders of magnitude smaller.

Events like inflation or phase transitions can produce large density perturbations on very small scales in the early Universe. Probes of small scales are therefore useful for e.g. discriminating between inflationary models. Until recently, the only such constraint came from non-observation of primordial black holes (PBHs), associated with the largest perturbations. Moderate-amplitude perturbations can collapse shortly after matter-radiation equality to form ultracompact minihalos (UCMHs) of dark matter, in far greater abundance than PBHs. If dark matter self-annihilates, UCMHs become excellent targets for indirect detection. Here we discuss the gamma-ray fluxes expected from UCMHs, the prospects of observing them with gamma-ray telescopes, and limits upon the primordial power spectrum derived from their non-observation by the Fermi Large Area Space Telescope.

Weak-scale supersymmetry is one of the most favoured theories beyond the Standard Model of particle physics that elegantly solves various theoretical and observational problems in both particle physics and cosmology. In this thesis, I describe the theoretical foundations of supersymmetry, issues that it can address and concrete supersymmetric models that are widely used in phenomenological studies. I discuss how the predictions of supersymmetric models may be compared with observational data from both colliders and cosmology. I show why constraints on supersymmetric parameters by direct and indirect searches of particle dark matter are of particular interest in this respect. Gamma-ray observations of astrophysical sources, in particular dwarf spheroidal galaxies, by the Fermi satellite, and recording nuclear recoil events and energies by future ton-scale direct detection experiments are shown to provide powerful tools in searches for supersymmetric dark matter and estimating supersymmetric parameters. I discuss some major statistical issues in supersymmetric global fits to experimental data. In particular, I further demonstrate that existing advanced scanning techniques may fail in correctly mapping the statistical properties of the parameter spaces even for the simplest supersymmetric models. Complementary scanning methods based on Genetic Algorithms are proposed.

Models of weak-scale supersymmetry offer viable dark matter (DM) candidates. Their parameter spaces are however rather large and complex, such that pinning down the actual parameter values from experimental data can depend strongly on the employed statistical framework and scanning algorithm. In frequentist parameter estimation, a central requirement for properly constructed confidence intervals is that they cover true parameter values, preferably at exactly the stated confidence level when experiments are repeated infinitely many times. Since most widely-used scanning techniques are optimised for Bayesian statistics, one needs to assess their abilities in providing correct confidence intervals in terms of the statistical coverage. Here we investigate this for the Constrained Minimal Supersymmetric Standard Model (CMSSM) when only constrained by data from direct searches for dark matter. We construct confidence intervals from one-dimensional profile likelihoods and study the coverage by generating several pseudo-experiments for a few benchmark sets of pseudo-true parameters. We use nested sampling to scan the parameter space and evaluate the coverage for the benchmarks when either flat or logarithmic priors are imposed on gaugino and scalar mass parameters. The sampling algorithm has been used in the configuration usually adopted for exploration of the Bayesian posterior. We observe both under- and over-coverage, which in some cases vary quite dramatically when benchmarks or priors are modified. We show how most of the variation can be explained as the impact of explicit priors as well as sampling effects, where the latter are indirectly imposed by physicality conditions. For comparison, we also evaluate the coverage for Bayesian credible intervals, and observe significant under-coverage in those cases.

Weakly interacting massive particles (WIMPs) are amongst the most interesting dark matter (DM) candidates. Many DM candidates naturally arise in theories beyond the standard model (SM) of particle physics, like weak-scale supersymmetry (SUSY). Experiments aim to detect WIMPs by scattering, annihilation or direct production, and thereby determine the underlying theory to which they belong, along with its parameters. Here we examine the prospects for further constraining the Constrained Minimal Supersymmetric Standard Model (CMSSM) with future ton-scale direct detection experiments. We consider ton-scale extrapolations of three current experiments: CDMS, XENON and COUPP, with 1000 kg-years of raw exposure each. We assume energy resolutions, energy ranges and efficiencies similar to the current versions of the experiments, and include backgrounds at target levels. Our analysis is based on full likelihood constructions for the experiments. We also take into account present uncertainties on hadronic matrix elements for neutralino-quark couplings, and on halo model parameters. We generate synthetic data based on four benchmark points and scan over the CMSSM parameter space using nested sampling. We construct both Bayesian posterior PDFs and frequentist profile likelihoods for the model parameters, as well as the mass and various cross-sections of the lightest neutralino. Future ton-scale experiments will help substantially in constraining supersymmetry, especially when results of experiments primarily targeting spin-dependent nuclear scattering are combined with those directed more toward spin-independent interactions.

The Constrained Minimal Supersymmetric Standard Model (CMSSM) is one of the simplest and most widely-studied supersymmetric extensions to the standard model of particle physics. Nevertheless, current data do not sufficiently constrain the model parameters in a way completely independent of priors, statistical measures and scanning techniques. We present a new technique for scanning supersymmetric parameter spaces, optimised for frequentist profile likelihood analyses and based on Genetic Algorithms. We apply this technique to the CMSSM, taking into account existing collider and cosmological data in our global fit. We compare our method to the MultiNest algorithm, an efficient Bayesian technique, paying particular attention to the best-fit points and implications for particle masses at the LHC and dark matter searches. Our global best-fit point lies in the focus point region. We find many high-likelihood points in both the stau co-annihilation and focus point regions, including a previously neglected section of the co-annihilation region at large m_0. We show that there are many high-likelihood points in the CMSSM parameter space commonly missed by existing scanning techniques, especially at high masses. This has a significant influence on the derived confidence regions for parameters and observables, and can dramatically change the entire statistical inference of such scans.

The dwarf galaxy Segue 1 is one of the most promising targets for the indirect detection of dark matter. Here we examine what constraints 9 months of Fermi-LAT gamma-ray observations of Segue 1 place upon the Constrained Minimal Supersymmetric Standard Model (CMSSM), with the lightest neutralino as the dark matter particle. We use nested sampling to explore the CMSSM parameter space, simultaneously fitting other relevant constraints from accelerator bounds, the relic density, electroweak precision observables, the anomalous magnetic moment of the muon and B-physics. We include spectral and spatial fits to the Fermi observations, a full treatment of the instrumental response and its related uncertainty, and detailed background models. We also perform an extrapolation to 5 years of observations, assuming no signal is observed from Segue 1 in that time. Results marginally disfavour models with low neutralino masses and high annihilation cross-sections. Virtually all of these models are however already disfavoured by existing experimental or relic density constraints.

We present a new model universe based on the junction of FRW to flat Lemaitre-Tolman-Bondi (LTB) solutions of Einstein equations along our past light cone, bringing structures within the FRW models. The model is assumed globally to be homogeneous, i.e. the cosmological principle is valid. Local inhomogeneities within the past light cone are modeled as a flat LTB, whereas those outside the light cone are assumed to be smoothed out and represented by a FRW model. The model is singularity free, always FRW far from the observer along the past light cone, gives way to a different luminosity distance relation as for the CDM/FRW models, a negative deceleration parameter near the observer, and correct linear and non-linear density contrast. As a whole, the model behaves like a FRW model on the past light cone with a special behavior of the scale factor, Hubble and deceleration parameter, mimicking dark energy.