Spectral distortions of the cosmic microwave background (CMB) provide a unique tool for learning about the early phases of cosmic history, reaching deep into the primordial Universe. At redshifts $z<10^6$, thermalization processes become inefficient and existing limits from COBE/FIRAS imply that no more than $\Delta \rho/\rho<6\times 10^{-5}$ (95% c.l.) of energy could have been injected into the CMB. However, at higher redshifts, when thermalization is efficient, the constraint weakens and $\Delta \rho/\rho \simeq 0.01-0.1$ could in principle have occurred. Existing computations for the evolution of distortions commonly assume $\Delta \rho/\rho \ll 1$ and thus become inaccurate in this case. Similarly, relativistic temperature corrections become relevant for large energy release, but have previously not been modeled as carefully. Here we study the evolution of distortions and the thermalization process after single large energy release at $z>10^5$. We show that for large distortions the thermalization efficiency is significantly reduced and that the distortion visibility is sizeable to much earlier times. This tightens spectral distortions constraints on low-mass primordial black holes with masses $M_{\rm PBH} < 6\times 10^{11}$ g. Similarly, distortion limits on the amplitude of the small-scale curvature power spectrum at wavenumbers $k>10^4\,{\rm Mpc}^{-1}$ and short-lived decaying particles with lifetimes $t_X< 10^7$ s are tightened, however, these require a more detailed time-dependent treatment. We also briefly discuss the constraints from measurements of the effective number of relativistic degrees of freedom and light element abundances and how these complement spectral distortion limits.

We study the double Compton (DC) process for a wide range of particle energies, extending previous treatments well beyond the soft photon limit, employing both numerical and analytical methods. This allows us to investigate the physics of the DC process up to the highly relativistic regime relevant to electromagnetic particle cascades in the early Universe and photon-dominated astrophysical plasmas. Generalized exact analytic expressions for the DC emissivity in the soft photon limit are obtained. These are compared to existing approximations, for the first time studying the ultra-relativistic regime. We also numerically integrate the full DC collision term calculating the DC emissivity at general particle energies. A careful treatment of DC infrared divergences inside astrophysical plasmas, including subtle effects related to the presence of stimulated DC emission, is discussed. The obtained results can be efficiently represented using the code DCpack, which also allows one to compute average emissivities for general incoming electron and photon distributions. This puts the modelling of the DC process inside astrophysical plasmas on a solid footing and should find applications in particular for computations of the cosmological thermalization problem in the early Universe.

We study the position-dependent power spectrum and the integrated bispectrum statistic for 2D cosmological fields on the sphere (integrated angular bispectrum). First, we derive a useful, $m$-independent, formula for the full-sky integrated angular bispectrum, based on the construction of azimuthally symmetric patches. We then implement a pipeline for integrated angular bispectrum estimation, including a mean-field correction to account for spurious isotropy-breaking effects in realistic conditions (e.g., inhomogenous noise, sky masking). Finally, we show examples of applications of this estimator to CMB analysis, both using simulations and actual Planck data. Such examples include $f_\mathrm{NL}$ estimation, analyses of non-Gaussianity from secondary anisotropies (ISW-lensing and ISW-tSZ-tSZ bispectra) and studies of non-Gaussian signatures from foreground contamination.

Electron-ion Bremsstrahlung (free-free) emission and absorption occur in many astrophysical plasmas for a wide range of physical conditions. This classical problem has been studied multiple times, and many analytical and numerical approximations exist. However, accurate calculations of the transition from the non-relativistic to the relativistic regime remain sparse. Here we provide a comprehensive study of the free-free Gaunt factors for ions with low charge (Z<=10). We compute the Gaunt factor using the expressions for the differential cross section given by Elwert & Haug (EH) and compare to various limiting cases. We develop a new software package, BRpack, for direct numerical applications. This package uses a combination of pre-computed tables and analytical approximations to efficiently cover a wide range of electron and photon energies, providing a representation of the EH Gaunt factor to better than 0.03% precision for Z<=2. Our results are compared to those of previous studies highlighting the improvements achieved here. BRpack should be useful in computations of spectral distortions of the cosmic microwave background, radiative transfer problems during reionization or inside galaxy clusters, and the modeling of galactic free-free foregrounds. The developed computational methods can furthermore be extended to higher energies and ion charge.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

Voyage 2050 White Paper highlighting the unique science opportunities using spectral distortions of the cosmic microwave background (CMB). CMB spectral distortions probe many processes throughout the history of the Universe. Precision spectroscopy, possible with existing technology, would provide key tests for processes expected within the cosmological standard model and open an enormous discovery space to new physics. This offers unique scientific opportunities for furthering our understanding of inflation, recombination, reionization and structure formation as well as dark matter and particle physics. A dedicated experimental approach could open this new window to the early Universe in the decades to come, allowing us to turn the long-standing upper distortion limits obtained with COBE/FIRAS some 25 years ago into clear detections of the expected standard distortion signals.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

We develop a linear perturbation theory for the spectral $y$-distortions of the cosmic microwave background (CMB). The $y$-distortions generated during the recombination epoch are usually negligible because the energy transfer due to the Compton scattering is strongly suppressed at that time, but they can be significant if there is a considerable amount of compensated isocurvature perturbation (CIP), which is not tightly constrained from the present CMB observations. The linear $y$-distortions explicitly depend on the baryon density fluctuations, therefore $y$ anisotropies can completely resolve the degeneracy between the baryon isocurvature perturbations and the cold dark matter ones. This novel method is free from lensing contaminations that can affect the previous approach to the CIPs based on the nonlinear modulation of the CMB anisotropies. We compute the cross correlation functions of the $y$-distortions with the CMB temperature and the $E$ mode polarization anisotropies. They are sensitive to the correlated CIPs parameterized by $f'\equiv\mathcal P_{\rm CIP\zeta}/\mathcal P_{\zeta \zeta}$ with $\mathcal P_{\zeta \zeta}$ and $\mathcal P_{\rm CIP\zeta}$ being the auto correlation of the adiabatic perturbations and the cross correlation between them and the CIPs. We investigate how well the $y$ anisotropies will constrain $f'$ in future observations such as those provided by a PIXIE-like and a PRISM-like survey, LiteBIRD and a cosmic variance limited (CVL) survey, taking into account the degradation in constraining power due to the presence of Sunyaev Zel'dovich effect from galaxy clusters. For example, our forecasts show that it is possible to achieve an upper limit of $f'< 2 \times 10^{5}$ at 68% C.L. with LiteBIRD, and $f'<2\times 10^{4}$ with CVL observations.

Cross-correlations between Cosmic Microwave Background (CMB) temperature and $y$-spectral distortions anisotropies have been previously proposed as a way to measure the local bispectrum parameter $f_{\rm NL}^{\rm loc.}$ in a range of scales inaccessible to either CMB ($T$, $E$) bispectra or $T$-$\mu$ correlations. This is useful e.g. to test scale dependence of primordial non-Gaussianity. Unfortunately, the primordial $y$-T signal is strongly contaminated by the late-time correlation between the Integrated Sachs Wolfe and Sunyaev-Zel'dovich (SZ) effects. Moreover, SZ itself generates a large noise contribution in the $y$-parameter map. We consider two original ways to address these issues. In order to remove the bias due to the SZ-CMB temperature coupling, while also adding new signal, we include in the analysis the cross-correlation between $y$-distortions and CMB {\em polarization}. In order to reduce the noise, we propose to clean the $y$-map by subtracting a SZ template, reconstructed via cross-correlation with external tracers (CMB and galaxy-lensing signals). We combine this SZ template subtraction with the previously adopted solution of directly masking detected clusters. Our final forecasts show that, using $y$-distortions, a PRISM-like survey can achieve $1\sigma(f_{\rm NL}^\text{loc.}) = 300$, while an ideal experiment will achieve $1\sigma(f_{\rm NL}^\text{loc.}) = 130$, with improvements of a factor $\sim 3$ from adding the $y$-$E$ signal, and a further $20-30 \%$ from template cleaning. These forecasts are much worse than current $f_{\rm NL}^\text{loc.}$ boundaries from {\em Planck}, but we stress again that they refer to completely different scales.

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

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We present a minimally parametric, model independent reconstruction of the shape of the primordial power spectrum. Our smoothing spline technique is well-suited to search for smooth features such as deviations from scale invariance, and deviations from a power law such as running of the spectral index or small-scale power suppression. We use a comprehensive set of the state-of the art cosmological data: {\it Planck} observations of the temperature and polarisation anisotropies of the cosmic microwave background, WiggleZ and Sloan Digital Sky Survey Data Release 7 galaxy power spectra and the Canada-France-Hawaii Lensing Survey correlation function. This reconstruction strongly supports the evidence for a power law primordial power spectrum with a red tilt and disfavours deviations from a power law power spectrum including small-scale power suppression such as that induced by significantly massive neutrinos. This offers a powerful confirmation of the inflationary paradigm, justifying the adoption of the inflationary prior in cosmological analyses.