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.

We investigate the possibility of measuring the primordial gravitational wave (GW) signal across 23 decades in frequencies, using the cosmic microwave background (CMB), pulsar timing arrays (PTA), and direct detection with laser and atomic interferometers. For the CMB and PTA experiments we consider the LiteBIRD mission and the Square Kilometer Array (SKA), respectively. For the interferometers we consider space mission proposals including the Laser Interferometer Space Antenna (LISA), the Big Bang Observer (BBO), the Deci-hertz Interferometer Gravitational wave Observatory (DECIGO), the $\mu$Ares experiment, the Decihertz Observatory (DO), and the Atomic Experiment for Dark Matter and Gravity Exploration in Space (AEDGE), as well as the ground-based Einstein Telescope (ET) and Cosmic Explorer (CE) proposals. We implement the mathematics needed to compute sensitivities for both CMB and interferometers, and derive the response functions for the latter from the first principles. We also evaluate the effect of the astrophysical foreground contamination in each experiment. We present binned sensitivity curves and error bars on the energy density parameter, $\Omega_{GW}h^2$, as a function of frequency for two representative classes of models for the stochastic background of primordial GW: the quantum vacuum fluctuation in the metric from single-field slow-roll inflation, and the source-induced tensor perturbation from the spectator axion-SU(2) inflation models. We find excellent prospects for joint measurements of the GW spectrum by CMB and space-borne direct detection mission proposals.

In this work we present a Neural Network (NN) algorithm for the identification of the appropriate parametrization of diffuse polarized Galactic emissions in the context of Cosmic Microwave Background (CMB) $B$-mode multi-frequency observations. In particular, we have focused our analysis on low frequency foregrounds relevant for polarization observation: namely Galactic Synchrotron and Anomalous Microwave Emission (AME). We have implemented and tested our approach on a set of simulated maps corresponding to the frequency coverage and sensitivity represented by future satellite and low frequency ground based probes. The NN efficiency in recognizing the right parametrization of foreground emission in different sky regions reaches an accuracy of about $90\%$. We have compared this performance with the $\chi^{2}$ information following parametric foreground estimation using multi-frequency fitting, and quantify the gain provided by a NN approach. Our results show the relevance of model recognition in CMB $B$-mode observations, and highlight the exploitation of dedicated procedures to this purpose.

We report a measurement of the E-mode polarization power spectrum of the cosmic microwave background (CMB) using 150 GHz data taken from July 2014 to December 2016 with the POLARBEAR experiment. We reach an effective polarization map noise level of $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ across an observation area of 670 square degrees. We measure the EE power spectrum over the angular multipole range $500 \leq \ell <3000$, tracing the third to seventh acoustic peaks with high sensitivity. The statistical uncertainty on E-mode bandpowers is $\sim 2.3 \mu {\rm K}^2$ at $\ell \sim 1000$ with a systematic uncertainty of 0.5$\mu {\rm K}^2$. The data are consistent with the standard $\Lambda$CDM cosmological model with a probability-to-exceed of 0.38. We combine recent CMB E-mode measurements and make inferences about cosmological parameters in $\Lambda$CDM as well as in extensions to $\Lambda$CDM. Adding the ground-based CMB polarization measurements to the Planck dataset reduces the uncertainty on the Hubble constant by a factor of 1.2 to $H_0 = 67.20 \pm 0.57 {\rm km\,s^{-1} \,Mpc^{-1}}$. When allowing the number of relativistic species ($N_{eff}$) to vary, we find $N_{eff} = 2.94 \pm 0.16$, which is in good agreement with the standard value of 3.046. Instead allowing the primordial helium abundance ($Y_{He}$) to vary, the data favor $Y_{He} = 0.248 \pm 0.012$. This is very close to the expectation of 0.2467 from Big Bang Nucleosynthesis. When varying both $Y_{He}$ and $N_{eff}$, we find $N_{eff} = 2.70 \pm 0.26$ and $Y_{He} = 0.262 \pm 0.015$.

The Large Scale Polarization Explorer (LSPE) is a cosmology program for the measurement of large scale curl-like features (B-modes) in the polarization of the Cosmic Microwave Background. Its goal is to constrain the background of inflationary gravity waves traveling through the universe at the time of matter-radiation decoupling. The two instruments of LSPE are meant to synergically operate by covering a large portion of the northern microwave sky. LSPE/STRIP is a coherent array of receivers planned to be operated from the Teide Observatory in Tenerife, for the control and characterization of the low-frequency polarized signals of galactic origin; LSPE/SWIPE is a balloon-borne bolometric polarimeter based on 330 large throughput multi-moded detectors, designed to measure the CMB polarization at 150 GHz and to monitor the polarized emission by galactic dust above 200 GHz. The combined performance and the expected level of systematics mitigation will allow LSPE to constrain primordial B-modes down to a tensor/scalar ratio of $10^{-2}$. We here report the status of the STRIP pre-commissioning phase and the progress in the characterization of the key subsystems of the SWIPE payload (namely the cryogenic polarization modulation unit and the multi-moded TES pixels) prior to receiver integration.

The search for primordial gravitational waves through the $B$-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary $B$-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper, we revise estimates of the reionization-induced $B$-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-$\alpha$ forest. To derive these $B$-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Ly$\alpha$ data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles $\ell\approx 100$, reionization induces $B$-mode power with $\ell(\ell+1)C_\ell^{BB}/2\pi\approx 4\times 10^{-6}\,\mu$K$^2$. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio $r<10^{-4}$, and can increase by a factor of $\sim 50$ if reionization is sourced by highly clustered sources residing in haloes with mass of $\sim 10^{11}$ M$_\odot$. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments.

Using only cosmic microwave background polarization data from the POLARBEAR experiment, we measure $B$-mode polarization delensing on subdegree scales at more than $5\sigma$ significance. We achieve a 14% $B$-mode power variance reduction, the highest to date for internal delensing, and improve this result to 2% by applying for the first time an iterative maximum a posteriori delensing method. Our analysis demonstrates the capability of internal delensing as a means of improving constraints on inflationary models, paving the way for the optimal analysis of next-generation primordial $B$-mode experiments.

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.

The cross-correlation between fluctuations in the electron scattering optical depth $\tau_{\rm es}$ as probed by future Cosmic Microwave Background (CMB) experiments, and fluctuations in the 21cm differential brightness temperature $\Delta T_{\rm 21cm}$ as probed by ground-based radio interferometers, will trace the reionization history of the Universe. In particular, the $\tau_{\rm es}-$21cm cross-correlation should yield a determination of the characteristic bubble size distribution and ionization fraction as a function of redshift. When assuming that the cross-correlation signal is limited by instrumental noise rather than by foregrounds, we estimate its potential detectability by upcoming experiments. Specifically, the combination of HERA and Simons Observatory, CMB-S4 and PICO should yield a signal-to-noise ratio around 3 - 6, while and the exploitation of the SKA should increase it to 10-20. Finally, we have discussed how such levels of detectability can be affected when (simply modeled) 21cm foregrounds are present. For the most promising PICO$\times$SKA configuration, an efficiency of foreground removal to a level of $7\times 10^{-4}$ is needed to achieve a $5\sigma$ detection of the cross-correlation signal; in addition, safe avoidance of foreground contamination in the line-of-sight Fourier modes above $0.03 \,h\,\rm Mpc^{-1}$ would guarantee a detection significance around $3\sigma$.

We present a measurement of the gravitational lensing deflection power spectrum reconstructed with two seasons cosmic microwave background polarization data from the POLARBEAR experiment. Observations were taken at 150 GHz from 2012 to 2014 which survey three patches of sky totaling 30 square degrees. We test the consistency of the lensing spectrum with a Cold Dark Matter (CDM) cosmology and reject the no-lensing hypothesis at a confidence of 10.9 sigma including statistical and systematic uncertainties. We observe a value of A_L = 1.33 +/- 0.32 (statistical) +/- 0.02 (systematic) +/- 0.07 (foreground) using all polarization lensing estimators, which corresponds to a 24% accurate measurement of the lensing amplitude. Compared to the analysis of the first year data, we have improved the breadth of both the suite of null tests and the error terms included in the estimation of systematic contamination.

We consider an Early Dark Energy (EDE) cosmological model, and perform an analysis which takes into account both background and perturbation effects via the parameters $c^{2}_{\rm eff}$ and $c^{2}_{\rm vis}$, representing effective sound speed and viscosity, respectively. By using the latest available data we derive constraints on the amount of dark energy at early times and the present value of the equation of state. Our focus is on the effect that early dark energy has on the Cosmic Microwave Background (CMB) data, including polarization and lensing, in a generalized parameter space including a varying total neutrino mass, and tensor to scalar ratio, besides the 6 standard parameters of the minimal cosmological model. We find that the inclusion of Baryonic Acoustic Oscillations (BAO) data and CMB lensing significantly improves the constraints on the EDE parameters, while other high redshift data like the Quasar Hubble diagram and the Lyman-$\alpha$ forest BAO have instead a negligible impact. We find $\Omega_{\rm eDE} < 0.00391 $ and $w_{0} < -0.9500 $ at the $95 \%$ C.L. for EDE accounting for its clustering through the inclusion of perturbation dynamics. This limit becomes weaker $\Omega_{\rm eDE} < 0.00337$ if perturbations are neglected. The constraints on the EDE parameters are remarkably stable even when $\Sigma m_{\nu}$, and $r$ parameters are varied, with weak degeneracies between $\Omega_{\rm eDE}$ and $r$ or $\Sigma m_{\nu}$. Perturbation parameters are not constrained with current data sets, and tensions between the CMB derived $H_0$ and $\sigma_8$ values and those measured with local probes are not eased. This work demonstrates the capability of CMB probes to constrain the total amount of EDE well below the percent level.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

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 report a 4.8$\sigma$ measurement of the cross-correlation signal between the cosmic microwave background (CMB) lensing convergence reconstructed from measurements of the CMB polarization made by the POLARBEAR experiment and the infrared-selected galaxies of the Herschel-ATLAS survey. This is the first measurement of its kind. We infer a best-fit galaxy bias of $b = 5.76 \pm 1.25$, corresponding to a host halo mass of $\log_{10}(M_h/M_\odot) =13.5^{+0.2}_{-0.3}$ at an effective redshift of $z \sim 2$ from the cross-correlation power spectrum. Residual uncertainties in the redshift distribution of the sub-mm galaxies are subdominant with respect to the statistical precision. We perform a suite of systematic tests, finding that instrumental and astrophysical contaminations are small compared to the statistical error. This cross-correlation measurement only relies on CMB polarization information that, differently from CMB temperature maps, is less contaminated by galactic and extra-galactic foregrounds, providing a clearer view of the projected matter distribution. This result demonstrates the feasibility and robustness of this approach for future high-sensitivity CMB polarization experiments.

Foreground components in the Cosmic Microwave Background (CMB) are sparse in a needlet representation, due to their specific morphological features (anisotropy, non-Gaussianity). This leads to the possibility of applying needlet thresholding procedures as a component separation tool. In this work, we develop algorithms based on different needlet-thresholding schemes and use them as extensions of existing, well-known component separation techniques, namely ILC and template-fitting. We test soft- and hard-thresholding schemes, using different procedures to set the optimal threshold level. We find that thresholding can be useful as a denoising tool for internal templates in experiments with few frequency channels, in conditions of low signal-to-noise. We also compare our method with other denoising techniques, showing that thresholding achieves the best performance in terms of reconstruction accuracy and data compression while preserving the map resolution. The best results in our tests are in particular obtained when considering template-fitting in an LSPE like experiment, especially for B-mode spectra.

We study how the shape of the spectrum of primordial gravitational waves can be constrained by future experiments looking at the B-mode of the Cosmic Microwave Background (CMB) polarization. We implement a Principal Component Analysis (PCA) including the effects of diffuse foreground residuals, following component separation, in the uncertainty of CMB angular power spectra, and taking into account the gravitational lensing by Large Scale Structure. We perform our study by considering the capabilities of future B-mode CMB experiments such as LiteBIRD, the Simons Observatory (SO) and Stage-IV (CMB-S4), in particular in detecting deviations of the primordial tensor spectrum from the scale-invariant behavior. We find that diffuse foreground residuals impact substantially both the derivation of the PCA basis and the corresponding constraining power, in all cases. In particular, depending on which experimental specifications and which value $r$ of tensor-to-scalar ratio for cosmological perturbations are considered, adding foregrounds residuals can determine an increase as large as a factor $\sim 4$ both on the uncertainty on $r$ and on the recovery of the PCA modes. We study the limitations of the methodology, including the effect of physicality priors on the PCA, which we quantify via a Monte Carlo Markov chain (MCMC) analysis of the combined cosmological and PCA power spectrum parameter space.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

We present a measurement of the $B$-mode polarization power spectrum of the cosmic microwave background (CMB) using taken from July 2014 to December 2016 with the POLARBEAR experiment. The CMB power spectra are measured using observations at 150 GHz with an instantaneous array sensitivity of $\mathrm{NET}_\mathrm{array}=23\, \mu \mathrm{K} \sqrt{\mathrm{s}}$ on a 670 square degree patch of sky centered at (RA, Dec)=($+0^\mathrm{h}12^\mathrm{m}0^\mathrm{s},-59^\circ18^\prime$). A continuously rotating half-wave plate is used to modulate polarization and to suppress low-frequency noise. We achieve $32\,\mu\mathrm{K}$-$\mathrm{arcmin}$ effective polarization map noise with a knee in sensitivity of $\ell = 90$, where the inflationary gravitational wave signal is expected to peak. The measured $B$-mode power spectrum is consistent with a $\Lambda$CDM lensing and single dust component foreground model over a range of multipoles $50 \leq \ell \leq 600$. The data disfavor zero $C_\ell^{BB}$ at $2.2\sigma$ using this $\ell$ range of POLARBEAR data alone. We cross-correlate our data with Planck high frequency maps and find the low-$\ell$ $B$-mode power in the combined dataset to be consistent with thermal dust emission. We place an upper limit on the tensor-to-scalar ratio $r < 0.90$ at 95% confidence level after marginalizing over foregrounds.

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.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

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).

We investigate the merging rates of compact binaries in galaxies, and the related detection rate of gravitational wave (GW) events with AdvLIGO/Virgo and with the Einstein Telescope. To this purpose, we rely on three basic ingredients: (i) the redshift-dependent galaxy statistics provided by the latest determination of the star formation rate functions from UV+far-IR/(sub)millimeter/radio data; (ii) star formation and chemical enrichment histories for individual galaxies, modeled on the basis of observations; (iii) compact remnant mass distribution and prescriptions for merging of compact binaries from stellar evolution simulations. We present results for the intrinsic birthrate of compact remnants, the merging rates of compact binaries, GW detection rates and GW counts, attempting to differentiate the outcomes among BH-BH, NS-NS, and BH-NS mergers, and to estimate their occurrence in disk and spheroidal host galaxies. We compare our approach with the one based on cosmic SFR density and cosmic metallicity, exploited by many literature studies; the merging rates from the two approaches are in agreement within the overall astrophysical uncertainties. We also investigate the effects of galaxy-scale strong gravitational lensing of GW in enhancing the rate of detectable events toward high-redshift. Finally, we discuss the contribution of undetected GW emission from compact binary mergers to the stochastic background.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

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 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).

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

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 Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

We present the mean spectral energy distribution (SED) of a sample of optically selected quasars (QSOs) at redshifts of $1 \le z \le 5$. To derive it, we exploit photometric information from SDSS, UKIDSS, and WISE surveys in combination with a stacking analysis of \textit{Herschel}, \textit{AKARI}, and \textit{Planck} maps at the location of the QSOs. The near-UV and optical parts of the reconstructed mean rest-frame SED are similar to those found in other studies. However, the SED shows an excess at 1-2 $\mu$m (when compared to the aforementioned SEDs normalized in the near-UV) and a prominent bump around 4-6 $\mu$m, followed by a decrease out to $\sim 20 \,\mu$m and a subsequent far-IR bump. From the fitted SEDs we estimate the average active galactic nuclei (AGN) luminosity $L_{\rm AGN}$ and star formation rate (SFR) as function of cosmic time, finding typical $L_{\rm AGN} \sim 10^{46} - 10^{47}$ erg/s and SFR $\sim 50 - 1000\, M_{\odot}/$yr. We develop mid-IR based criteria to split the QSO sample, finding that these allow us to move along the average relationship in the SFR vs. $L_{\rm AGN}$ diagram toward increasing AGN luminosities. When interpreted in the context of the in-situ coevolution scenario presented by Lapi et al. 2014, our results suggest that the detection in the far-IR band is an effective criterion to select objects where the star formation is on the verge of being affected by energy/momentum feedback from the central AGN.

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.

The polarization of the atmosphere has been a long-standing concern for ground-based experiments targeting cosmic microwave background (CMB) polarization. Ice crystals in upper tropospheric clouds scatter thermal radiation from the ground and produce a horizontally-polarized signal. We report the detailed analysis of the cloud signal using a ground-based CMB experiment, POLARBEAR, located at the Atacama desert in Chile and observing at 150 GHz. We observe horizontally-polarized temporal increases of low-frequency fluctuations ("polarized bursts," hereafter) of $\lesssim$0.1 K when clouds appear in a webcam monitoring the telescope and the sky. The hypothesis of no correlation between polarized bursts and clouds is rejected with $>$24$\sigma$ statistical significance using three years of data. We consider many other possibilities including instrumental and environmental effects, and find no other reasons other than clouds that can explain the data better. We also discuss the impact of the cloud polarization on future ground-based CMB polarization experiments.

We analyze the Southern Sky emission in linear polarization at 2.3 GHz as observed by the S-band Polarization All Sky Survey S-PASS. Our purpose is to study the properties of the diffuse Galactic polarized synchrotron as a contaminant to CMB B-mode observations. We study the angular distribution of the S-PASS signal at intermediate and high Galactic latitudes by means of angular power spectra. Power spectra, show a decay of the spectral amplitude as a function of multipole for \ell<200, typical of the diffuse emission. Spectra can be approximated by a power law C_{\ell}\propto\ell^{alpha}, with alpha~-3, and characterized by a B-to-E ratio of ~0.5. We study the synchrotron SED in polarization by computing power spectra of the low frequency WMAP and Planck maps. Results show that the SED, in the frequency range 2.3-33 GHz, is compatible with a power law with beta_s=-3.22\pm0.08. Combining S-PASS polarization maps with those coming from WMAP and Planck we derived a map of the synchrotron spectral index at angular resolution of 2{\deg} on about 30% of the sky. The recovered distribution peaks at the value around -3.2. We also measure a significant spatial correlation between synchrotron and thermal dust signals, as traced by the Planck 353 GHz channel. This correlation reaches about 40% on the larger angular scales, decaying considerably at the degree scales. Finally, we use the S-PASS maps to assess the polarized synchrotron contamination to CMB observations of the B-modes. Moreover, by combining S-PASS data with Planck 353 GHz observations, we recover a map of the minimum level of total polarized foreground contamination to B-modes, finding that there is no region of the sky, at any frequency, where this contamination lies below equivalent tenor-to-scalar ratio ~10^-3. This result confirms the importance of observing both high and low frequency foregrounds in CMB B-mode measurements.

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.

The Simons Observatory (SO) will make precise temperature and polarization measurements of the cosmic microwave background (CMB) using a set of telescopes which will cover angular scales between 1 arcminute and tens of degrees, contain over 60,000 detectors, and observe at frequencies between 27 and 270 GHz. SO will consist of a 6 m aperture telescope coupled to over 30,000 transition-edge sensor bolometers along with three 42 cm aperture refractive telescopes, coupled to an additional 30,000+ detectors, all of which will be located in the Atacama Desert at an altitude of 5190 m. The powerful combination of large and small apertures in a CMB observatory will allow us to sample a wide range of angular scales over a common survey area. SO will measure fundamental cosmological parameters of our universe, constrain primordial fluctuations, find high redshift clusters via the Sunyaev-Zel`dovich effect, constrain properties of neutrinos, and trace the density and velocity of the matter in the universe over cosmic time. The complex set of technical and science requirements for this experiment has led to innovative instrumentation solutions which we will discuss. The large aperture telescope will couple to a cryogenic receiver that is 2.4 m in diameter and nearly 3 m long, creating a number of technical challenges. Concurrently, we are designing the array of cryogenic receivers housing the 42 cm aperture telescopes. We will discuss the sensor technology SO will use and we will give an overview of the drivers for and designs of the SO telescopes and receivers, with their cold optical components and detector arrays.

EBEX was a long-duration balloon-borne experiment to measure the polarization of the cosmic microwave background. The experiment had three frequency bands centered at 150, 250, and 410 GHz and was the first to use a kilo-pixel array of transition edge sensor (TES) bolometers aboard a balloon platform; shortly after reaching float we operated 504, 342, and 109 TESs at each of the bands, respectively. We describe the design and characterization of the array and the readout system. We give the distributions of measured thermal conductances, normal resistances, and transition temperatures. With the exception of the thermal conductance at 150 GHz. We measured median low-loop-gain time constants $\tau_{0}=$ 88, 46, and 57 ms and compare them to predictions. Two measurements of bolometer absorption efficiency show high ($\sim$0.9) efficiency at 150 GHz and medium ($\sim$0.35, and $\sim$0.25) at the two higher bands, respectively. We measure a median total optical load of 3.6, 5.3 and 5.0 pW absorbed at the three bands, respectively. EBEX pioneered the use of the digital version of the frequency domain multiplexing (FDM) system which multiplexed the bias and readout of 16 bolometers onto two wires. We present accounting of the measured noise equivalent power. The median per-detector noise equivalent temperatures referred to a black body with a temperature of 2.725 K are 400, 920, and 14500 $\mu$K$\sqrt{s}$ for the three bands, respectively. We compare these values to our pre-flight predictions and to a previous balloon payload, discuss the sources of excess noise, and the path for a future payload to make full use of the balloon environment.

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.

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.

An estimation of the sky signal from streams of Time Ordered Data (TOD) acquired by Cosmic Microwave Background (\cmb) experiments is one of the most important steps in the context of \cmb data analysis referred to as the map-making problem. The continuously growing \cmb data sets render the \cmb map-making problem more challenging in terms of computational cost and memory in particular in the context of ground based experiments. In this context, we study a novel class of the Preconditioned Conjugate Gradient (PCG) solvers which invoke two-level preconditioners. We compare them against PCG solvers commonly used in the map-making context considering their precision and time-to-solution. We compare these new methods on realistic, simulated data sets reflecting the characteristics of current and forthcoming \cmb ground-based experiment. We develop an embarrassingly parallel implementation of the approach where each processor performs a sequential map-making for a subset of the TOD. We find that considering the map level residuals the new class of solvers permits achieving tolerance of up to 3 orders of magnitude better than the standard approach, where the residual level often saturates before convergence is reached. This corresponds to an important improvement in the precision of recovered power spectra in particular on the largest angular scales. The new method also typically requires fewer iterations to reach a required precision and thus shorter runtimes for a single map-making solution. However, the construction of an appropriate two-level preconditioner can be as costly as a single standard map-making run. Nevertheless, if the same problem needs to be solved multiple times, e.g., as in Monte Carlo simulations, this cost has to be incurred only once, and the method should be competitive not only as far as its precision but also its performance is concerned.

We study the signal from patchy reionization in view of the future high accuracy polarization measurements of the Cosmic Microwave Background (CMB). We implement an extraction procedure of the patchy reionization signal analogous to CMB lensing. We evaluate the signal to noise ratio (SNR) for the future Stage IV (S4) CMB experiment. The signal has a broad peak centered on the degree angular scales, with a long tail at higher multipoles. The CMB S4 experiment can effectively constrain the properties of reionization by measuring the signal on degree scales. The signal amplitude depends on the properties of the structure determining the reionization morphology. We describe bubbles having radii distributed log-normally. The expected S/N is sensitive to the mean bubble radius: $\bar{R}=5$ Mpc implies $S/N \approx 4$, $\bar{R}=10$ Mpc implies $S/N \approx 20$. The spread of the radii distribution strongly affects the integrated SNR, that changes by a factor of $10^2$ when $\sigma_{lnr}$ goes from $\ln 2$ to $\ln3$. Future CMB experiments will thus place important constraints on the physics of reionization.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

We combine the latest datasets obtained with different surveys to study the frequency dependence of polarized emission coming from Extragalactic Radio Sources (ERS). We consider data over a very wide frequency range starting from $1.4$ GHz up to $217$ GHz. This range is particularly interesting since it overlaps the frequencies of the current and forthcoming Cosmic Microwave Background (\cmb) experiments. Current data suggest that at high radio frequencies, ($ \nu \geq 20$ GHz) the fractional polarization of ERS does not depend on the total flux density. Conversely, recent datasets indicate a moderate increase of polarization fraction as a function of frequency, physically motivated by the fact that Faraday depolarization is expected to be less relevant at high radio-frequencies. We compute ERS number counts using updated models based on recent data, and we forecast the contribution of unresolved ERS in CMB polarization spectra. Given the expected sensitivities and the observational patch sizes of forthcoming \cmb\ experiments about $\sim 200 $ ( up to $\sim 2000 $ ) polarized ERS are expected to be detected. Finally, we assess that polarized ERS can contaminate the cosmological B-mode polarization if the tensor-to-scalar ratio is $r< 0.05$ and they have to be robustly controlled to de-lens \cmb\ B-modes at the arcminute angular scales.

We discuss a systematic effect associated with measuring polarization with a continuously rotating half-wave plate. The effect was identified with the data from the E and B Experiment (EBEX), which was a balloon-borne instrument designed to measure the polarization of the CMB as well as that from Galactic dust. The data show polarization fraction larger than 10\% while less than 3\% were expected from instrumental polarization. We give evidence that the excess polarization is due to detector non-linearity in the presence of a continuously rotating HWP. The non-linearity couples intensity signals into polarization. We develop a map-based method to remove the excess polarization. Applying this method for the 150 (250) GHz bands data we find that 81\% (92\%) of the excess polarization was removed. Characterization and mitigation of this effect is important for future experiments aiming to measure the CMB B-modes with a continuously rotating HWP.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

We report an improved measurement of the cosmic microwave background (CMB) $B$-mode polarization power spectrum with the POLARBEAR experiment at 150 GHz. By adding new data collected during the second season of observations (2013-2014) to re-analyzed data from the first season (2012-2013), we have reduced twofold the band-power uncertainties. The band powers are reported over angular multipoles $500 \leq \ell \leq 2100$, where the dominant $B$-mode signal is expected to be due to the gravitational lensing of $E$-modes. We reject the null hypothesis of no $B$-mode polarization at a confidence of 3.1$\sigma$ including both statistical and systematic uncertainties. We test the consistency of the measured $B$-modes with the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) framework by fitting for a single lensing amplitude parameter $A_L$ relative to the Planck best-fit model prediction. We obtain $A_L = 0.60 ^{+0.26} _{-0.24} ({\rm stat}) ^{+0.00} _{-0.04}({\rm inst}) \pm 0.14 ({\rm foreground}) \pm 0.04 ({\rm multi})$, where $A_{L}=1$ is the fiducial $\Lambda$CDM value, and the details of the reported uncertainties are explained later in the manuscript.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

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

The E and B Experiment (EBEX) was a long-duration balloon-borne cosmic microwave background polarimeter that flew over Antarctica in 2013. We describe the experiment's optical system, receiver, and polarimetric approach, and report on their in-flight performance. EBEX had three frequency bands centered on 150, 250, and 410 GHz. To make efficient use of limited mass and space we designed a 115 cm$^{2}$sr high throughput optical system that had two ambient temperature mirrors and four anti-reflection coated polyethylene lenses per focal plane. All frequency bands shared the same optical train. Polarimetry was achieved with a continuously rotating achromatic half-wave plate (AHWP) that was levitated with a superconducting magnetic bearing (SMB). Rotation stability was 0.45 % over a period of 10 hours, and angular position accuracy was 0.01 degrees. This is the first use of a SMB in astrophysics. The measured modulation efficiency was above 90 % for all bands. To our knowledge the 109 % fractional bandwidth of the AHWP was the broadest implemented to date. The receiver that contained one lens and the AHWP at a temperature of 4 K, the polarizing grid and other lenses at 1 K, and the two focal planes at 0.25 K performed according to specifications giving focal plane temperature stability with fluctuation power spectrum that had $1/f$ knee at 2 mHz. EBEX was the first balloon-borne instrument to implement technologies characteristic of modern CMB polarimeters including high throughput optical systems, and large arrays of transition edge sensor bolometric detectors with mutiplexed readouts.

A continuously rotating half-wave plate (CRHWP) is a promising tool to improve the sensitivity to large angular scales in cosmic microwave background (CMB) polarization measurements. With a CRHWP, single detectors can measure three of the Stokes parameters, $I$, $Q$ and $U$, thereby avoiding the set of systematic errors that can be introduced by mismatches in the properties of orthogonal detector pairs. We focus on the implementation of CRHWPs in large aperture telescopes (i.e. the primary mirror is larger than the current maximum half-wave plate diameter of $\sim$0.5 m), where the CRHWP can be placed between the primary mirror and focal plane. In this configuration, one needs to address the intensity to polarization ($I{\rightarrow}P$) leakage of the optics, which becomes a source of 1/f noise and also causes differential gain systematics that arise from CMB temperature fluctuations. In this paper, we present the performance of a CRHWP installed in the POLARBEAR experiment, which employs a Gregorian telescope with a 2.5 m primary illumination pattern. The CRHWP is placed near the prime focus between the primary and secondary mirrors. We find that the $I{\rightarrow}P$ leakage is larger than the expectation from the physical properties of our primary mirror, resulting in a 1/f knee of 100 mHz. The excess leakage could be due to imperfections in the detector system, i.e. detector non-linearity in the responsivity and time-constant. We demonstrate, however, that by subtracting the leakage correlated with the intensity signal, the 1/f noise knee frequency is reduced to 32 mHz ($\ell \sim$39 for our scan strategy), which is very promising to probe the primordial B-mode signal. We also discuss methods for further noise subtraction in future projects where the precise temperature control of instrumental components and the leakage reduction will play a key role.

We present a model for simulating Carbon Monoxide (CO) rotational line emission in molecular clouds, taking account of their 3D spatial distribution in galaxies with different geometrical properties. The model implemented is based on recent results in the literature and has been designed for performing Monte-Carlo simulations of this emission. We compare the simulations produced with this model and calibrate them, both on the map level and on the power spectrum level, using the second release of data from the Planck satellite for the Galactic plane, where the signal-to-noise ratio is highest. We use the calibrated model to extrapolate the CO power spectrum at low Galactic latitudes where no high sensitivity observations are available yet. We then forecast the level of unresolved polarized emission from CO molecular clouds which could contaminate the power spectrum of Cosmic Microwave Background (CMB) polarization B-modes away from the Galactic plane. Assuming realistic levels of the polarization fraction, we show that the level of contamination is equivalent to a cosmological signal with $r \lesssim 0.02$. The Monte-Carlo MOlecular Line Emission (MCMole3D) Python package, which implements this model, is being made publicly available.

The E and B Experiment (EBEX) was a long-duration balloon-borne instrument designed to measure the polarization of the cosmic microwave background (CMB) radiation. EBEX was the first balloon-borne instrument to implement a kilo-pixel array of transition edge sensor (TES) bolometric detectors and the first CMB experiment to use the digital version of the frequency domain multiplexing system for readout of the TES array. The scan strategy relied on 40 s peak-to-peak constant velocity azimuthal scans. We discuss the unique demands on the design and operation of the payload that resulted from these new technologies and the scan strategy. We describe the solutions implemented including the development of a power system designed to provide a total of at least 2.3 kW, a cooling system to dissipate 590 W consumed by the detectors' readout system, software to manage and handle the data of the kilo-pixel array, and specialized attitude reconstruction software. We present flight performance data showing faultless management of the TES array, adequate powering and cooling of the readout electronics, and constraint of attitude reconstruction errors such that the spurious B-modes they induced were less than 10% of CMB B-mode power spectrum with $r=0.05$.

We investigate the Alcock Paczynski (AP) test applied to the Baryon Acoustic Oscillation (BAO) feature in the galaxy correlation function. By using a general formalism that includes relativistic effects, we quantify the importance of the linear redshift space distortions and gravitational lensing corrections to the galaxy number density fluctuation. We show that redshift space distortions significantly affect the shape of the correlation function, both in radial and transverse directions, causing different values of galaxy bias to induce offsets up to 1% in the AP test. On the other hand, we find that the lensing correction around the BAO scale modifies the amplitude but not the shape of the correlation function and therefore does not introduce any systematic effect. Furthermore, we investigate in details how the AP test is sensitive to redshift binning: a window function in transverse direction suppresses correlations and shifts the peak position toward smaller angular scales. We determine the correction that should be applied in order to account for this effect, when performing the test with data from three future planned galaxy redshift surveys: Euclid, the Dark Energy Spectroscopic Instrument (DESI) and the Square Kilometer Array (SKA).

Analysis of cosmic microwave background (CMB) datasets typically requires some filtering of the raw time-ordered data. Filtering is frequently used to minimize the impact of low frequency noise, atmospheric contributions and/or scan synchronous signals on the resulting maps. In this work we explicitly construct a general filtering operator, which can unambiguously remove any set of unwanted modes in the data, and then amend the map-making procedure in order to incorporate and correct for it. We show that such an approach is mathematically equivalent to the solution of a problem in which the sky signal and unwanted modes are estimated simultaneously and the latter are marginalized over. We investigate the conditions under which this amended map-making procedure can render an unbiased estimate of the sky signal in realistic circumstances. We then study the effects of time-domain filtering on the noise correlation structure in the map domain, as well as impact it may have on the performance of the popular pseudo-spectrum estimators. We conclude that although maps produced by the proposed estimators arguably provide the most faithful representation of the sky possible given the data, they may not straightforwardly lead to the best constraints on the power spectra of the underlying sky signal and special care may need to be taken to ensure this is the case. By contrast, simplified map-makers which do not explicitly correct for time-domain filtering, but leave it to subsequent steps in the data analysis, may perform equally well and be easier and faster to implement. We focus on polarization-sensitive measurements targeting the B-mode component of the CMB signal and apply the proposed methods to realistic simulations based on characteristics of an actual CMB polarization experiment, POLARBEAR.

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.

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than $3\,\sigma$. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

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.

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.