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 Laser Interferometer Space Antenna will be the first Gravitational Wave observatory in space. It is scheduled to fly in the early 2030's. LISA design predicts sensitivity levels that enable the detection a Stochastic Gravitational Wave Background signal. This stochastic type of signal is a superposition of signatures from sources that cannot be resolved individually and which are of various types, each one contributing with a different spectral shape. In this work we present a fast methodology to assess the detectability of a stationary, Gaussian, and isotropic stochastic signal in a set of frequency bins, combining information from the available data channels. We derive an analytic expression of the Bayes Factor between the instrumental noise-only and the signal plus instrumental noise models, that allows us to compute the detectability bounds of a given signal, as a function of frequency and prior knowledge on the instrumental noise spectrum.

The Laser Interferometer Space Antenna (LISA) is a European Space Agency mission that aims to measure gravitational waves in the millihertz range. Laser frequency noise enters the interferometric measurements and dominates the expected gravitational signals by many orders of magnitude. Time-delay interferometry (TDI) is a technique that reduces this laser noise by synthesizing virtual equal-arm interferometric measurements. Laboratory experiments and numerical simulations have confirmed that this reduction is sufficient to meet the scientific goals of the mission in proof-of-concept setups. In this paper, we show that the on-board antialiasing filters play an important role in TDI's performance when the flexing of the constellation is accounted for. This coupling was neglected in previous studies. To reach an optimal reduction level, filters with vanishing group delays must be used on board or synthesized off-line. We propose a theoretical model of the residual laser noise including this flexing-filtering coupling. We also use two independent simulators to produce realistic measurement signals and compute the corresponding TDI Michelson variables. We show that our theoretical model agrees with the simulated data with exquisite precision. Using these two complementary approaches, we confirm TDI's ability to reduce laser frequency noise in a more realistic mission setup. The theoretical model provides insight on filter design and implementation.

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)

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)

Under certain conditions the collision and intercommutation of two cosmic strings can result in the formation of a third string, with the three strings then remaining connected at Y-junctions. The kinematics and dynamics of collisions of this type have been the subject of analytical and numerical analyses in the special case in which the strings are Nambu-Goto. Cosmic strings, however, may well carry currents, in which case their dynamics is not given by the Nambu-Goto action. Our aim is to extend the kinematic analysis to more general kinds of string model. We focus in particular on the collision of strings described by conservative elastic string models, characteristic of current carrying strings, and which are expected to form in a cosmological context. As opposed to Nambu-Goto strings collisions, we show that in this case the collision cannot lead to the formation of a third elastic string: if dynamically such a string forms then the joining string must be described by a more general equation of state. This process will be studied numerically in a forthcoming publication.

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]

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20{\deg}$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $\beta = 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $\beta = 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{\deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $\Lambda$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5{\sigma} from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($\mu$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $\Lambda$CDM parameters consistent with those from temperature alone.

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

Although the inflationary paradigm is the most widely accepted explanation for the current cosmological observations, it does not necessarily correspond to what actually happened in the early stages of our Universe. To decide on this issue, two paths can be followed: first, all the possible predictions it makes must be derived thoroughly and compared with available data, and second, all imaginable alternatives must be ruled out. Leaving the first task to all other contributors of this volume, we concentrate here on the second option, focusing on the bouncing alternatives and their consequences.

The simplest possible classical model leading to a cosmological bounce is examined in the light of the non-Gaussianities it can generate. Concentrating solely on the transition between contraction and expansion, and assuming initially purely Gaussian perturbations at the end of the contracting phase, we find that the bounce acts as a source such that the resulting value for the post-bounce $f_{\mathrm{NL}}$ may largely exceed all current limits, to the point of potentially casting doubts on the validity of the perturbative expansion. We conjecture that if one can assume that the non-Gaussianity production depends only on the bouncing behavior of the scale factor and not on the specifics of the model examined, then many realistic models in which a nonsingular classical bounce takes place could exhibit a generic non-Gaussianity excess problem that would need to be addressed for each case.

We compute the level of non-gaussianities produced by a cosmological bouncing phase in the minimal non-singular setup that lies within the context of General Relativity when the matter content consists of a simple scalar field with a standard kinetic term. Such a bouncing phase is obtained by requiring that the spatial sections of the background spacetime be positively curved. We restrict attention to the close vicinity of the bounce by Taylor expanding the scale factor, the scalar field and its potential in powers of the conformal time around the bounce. We find that possibly large non-gaussianities are generically produced at the bounce itself and also discuss which shapes of non-gaussianities are mostly likely to be produced.

We study a cosmological scenario in which inflation is preceded by a bounce. In this scenario, the primordial singularity, one of the major shortcomings of inflation, is replaced by a non-singular bounce, prior to which the universe undergoes a phase of contraction. Our starting point is the bouncing cosmology investigated in Falciano et al. (2008), which we complete by a detailed study of the transfer of cosmological perturbations through the bounce and a discussion of possible observational effects of bouncing cosmologies. We focus on a symmetric bounce and compute the evolution of cosmological perturbations during the contracting, bouncing and inflationary phases. We derive an expression for the Mukhanov-Sasaki perturbation variable at the onset of the inflationary phase that follows the bounce. Rather than being in the Bunch-Davies vacuum, it is found to be in an excited state that depends on the time scale of the bounce. We then show that this induces oscillations superimposed on the nearly scale-invariant primordial spectra for scalar and tensor perturbations. We discuss the effects of these oscillations in the cosmic microwave background and in the matter power spectrum. We propose a new way to indirectly measure the spatial curvature energy density parameter in the context of this model.

A nonabelian generalization of the neutral Witten current-carrying string model is discussed in which the bosonic current-carrier belongs to a two dimensional representation of SU(2). We find that the current-carrying solutions can be of three different kinds: either the current spans a U(1) subgroup, and in which case one is left with an abelian current-carrying string, or the three currents are all lightlike, travelling in the same direction (only left or right movers). The third, genuinely nonabelian situation, cannot be handled within a cylindrically symmetric framework, but can be shown to depend on all possible string Lorentz invariant quantities that can be constructed out of the phase gradients.

We first examine the microstructure of a cosmic string endowed with two simple Abelian currents. This microstructure depends on two state parameters. We then provide the macroscopic description of such a string and show that it depends on an additional Lorentz-invariant state parameter that relates the two currents. We find that in most of the parameter space, the two-current string is essentially equivalent to the single current-carrying string, i.e., only one field condenses onto the defect. In the regions where two currents are present, we find that as far as stability is concerned, one can approximate the dynamics with good accuracy using an analytic model based on either a logarithmic (on the electric side, i.e., for timelike currents) or a rational (on the magnetic side, i.e., for spacelike currents) worldsheet Lagrangian.

The aim of these lectures is to give an introduction to several topics which lie at the intersection of string theory, gravity theory and gravity phenomenology. One successively reviews: (i) the "membrane" approach to the dissipative dynamics of classical black holes, (ii) the current experimental tests of gravity, and their theoretical interpretation, (iii) some aspects of the string-inspired phenomenology of the gravitational sector, and (iv) some possibilities for observing string-related signals in cosmology (including a discussion of gravitational wave signals from cosmic superstrings).

Gravitational waves generated by the final merger of double neutron star (DNS) binary systems are a key target for the gravitational wave (GW) interferometric detectors, such as LIGO, and the next generation detectors, Advanced LIGO. The cumulative GW signal from DNS mergers in interferometric data will manifest as "geometrical noise": a non-continuous stochastic background with a unique statistical signature dominated by the spatial and temporal distribution of the sources. Because geometrical noise is highly non-Gaussian, it could potentially be used to identify the presence of a stochastic GW background from DNS mergers. We demonstrate this by fitting to a simulated distribution of transients using a model for the DNS merger rate and idealized Gaussian detector noise. Using the cosmological "probability event horizon" concept and recent bounds for the Galactic DNS merger rate, we calculate the evolution of the detectability of DNS mergers with observation time. For Advanced LIGO sensitivities and a detection threshold assuming optimal filtering, there is a 95% probability that a minimum of one DNS merger signal will be detectable from the ensemble of events comprising the stochastic background during 12-211 days of observation. For initial LIGO sensitivities, we identify an interesting regime where there is a 95% probability that at least one DNS merger with signal-to-noise ratio > unity will occur during 4-68 days of observation. We propose that there exists an intermediate detection regime with pre-filtered signal-noise-ratio less than unity, where the DNS merger rate is high enough that the geometrical signature could be identified in interferometer data.

The locations of the peaks of the CMB spectrum are sensitive indicators of cosmological parameters, yet there is no known analytic formula which accurately describes their dependence on them. We parametrize the location of the peaks as l_m = l_A(m - \phi_m), where l_A is the analytically calculable acoustic scale and m labels the peak number. Fitting formulae for the phase shifts \phi_m for the first three peaks and the first trough are given. It is shown that in a wide range of parameter space, the acoustic scale l_A can be retrieved from actual CMB measurements of the first three peaks within one percent accuracy. This can be used to speed up likelihood analysis. We describe how the peak shifts can be used to distinguish between different models of dark energy.

The CMB data recently released by BOOMERANG and MAXIMA suggest that the anisotropy spectrum has a third peak in the range 800<l_3<900. A combination of this result with constraints from large-scale structure permit us to differentiate between different quintessence models. In particular, we find that inverse power law models with power \alpha >1 are disfavoured. Models with more than 5% quintessence before last scattering require a spectral index greater than 1. These constraints are compared with supernovae observations. We also show that the CMB alone now provides strong evidence for an accelerating universe.

A hypothetical dark energy component may have an equation of state that is different from a cosmological constant and possibly even changing in time. The spacing of the cosmic microwave background peaks is sensitive to the ratio of the horizon sizes today and at last scattering and therefore to the details of the cosmological evolution. Together with independent measurements of today's cosmological parameters, this restricts quintessence models in the epoch before and during decoupling. For fixed \Omega^b_0, \Omega^{cdm}_0 and h, we give an analytic estimate of the spacing that depends on \Omega^{\phi}_{ls}, \Omega^{\phi}_0 and an effective equation of state of the dark energy component.

We propose that it should be possible to use the CMB to discriminate between dark energy models with different equations of state, including distinguishing a cosmological constant from many models of quintessence. The separation of peaks in the CMB anisotropies can be parametrised by three quantities: the amount of quintessence today, the amount at last scattering, and the averaged equation of state of quintessence. In particular, we show that the CMB peaks can be used to measure the amount of dark energy present before last scattering.

In cosmological first-order phase transitions, the progress of true-vacuum bubbles is expected to be significantly retarded by the interaction between the bubble wall and the hot plasma. We examine the evolution and collision of slow-moving true-vacuum bubbles. Our lattice simulations indicate that phase oscillations, predicted and observed in systems with a local symmetry and with a global symmetry where the bubbles move at speeds less than the speed of light, do not occur inside collisions of slow-moving local-symmetry bubbles. We observe almost instantaneous phase equilibration which would lead to a decrease in the expected initial defect density, or possibly prevent defects from forming at all. We illustrate our findings with an example of defect formation suppressed in slow-moving bubbles. Slow-moving bubble walls also prevent the formation of `extra defects', and in the presence of plasma conductivity may lead to an increase in the magnitude of any primordial magnetic field formed.

We point out that the lower bound on the primordial magnetic field required to seed the galactic dynamo is significantly relaxed in an open universe or in a universe with a positive cosmological constant. It is shown that, for reasonable cosmological parameters, primordial seed fields of strength 10^{-30} Gauss or less at the time of galaxy formation could explain observed galactic magnetic fields. As a consequence, mechanisms of primordial magnetic seed-field generation that have previously been ruled out could well be viable. We also comment on the implications of the observation of micro-Gauss magnetic fields in galaxies at high redshift.