Galactic foregrounds are the main obstacle to observations of the cosmic microwave background (CMB) $B$-mode polarization. In addition to obscuring the inflationary $B$-mode signal by several orders of magnitude, Galactic foregrounds have non-trivial spectral signatures that are partially unknown and distorted by averaging effects along the line-of-sight, within the pixel/beam window, and by various analysis choices (e.g., spherical harmonic transforms and filters). Statistical moment expansion methods provide a powerful tool for modeling the effective Galactic foreground emission resulting from these averaging effects in CMB observations, while blind component separation treatments can handle unknown foregrounds. In this work, we combine these two approaches to develop a new semi-blind component separation method at the intersection of parametric and blind methods, called constrained moment ILC (cMILC). This method adds several constraints to the standard ILC method to de-project the main statistical moments of the Galactic foreground emission. Applications to maps are performed in needlet space and when compared to the NILC method, this helps significantly reducing residual foreground contamination (bias, variance, and skewness) in the reconstructed CMB $B$-mode map, power spectrum, and tensor-to-scalar ratio. We consider sky-simulations for experimental settings similar to those of LiteBIRD and PICO, illustrating which trade-offs between residual foreground biases and degradation of the constraint on $r$ can be expected within the new cMILC framework. We also outline several directions that require more work in preparation for the coming analysis challenges.

Observations of the CMB today allow us to answer detailed questions about the properties of our Universe, targeting both standard and non-standard physics. In this paper, we study the effects of varying fundamental constants (i.e., the fine-structure constant, $\alpha_{\rm EM}$, and electron rest mass, $m_{\rm e}$) around last scattering using the recombination codes CosmoRec and Recfast++. We approach the problem in a pedagogical manner, illustrating the importance of various effects on the free electron fraction, Thomson visibility function and CMB power spectra, highlighting various degeneracies. We demonstrate that the simpler Recfast++ treatment (based on a three-level atom approach) can be used to accurately represent the full computation of CosmoRec. We also include explicit time-dependent variations using a phenomenological power-law description. We reproduce previous Planck 2013 results in our analysis. Assuming constant variations relative to the standard values, we find the improved constraints $\alpha_{\rm EM}/\alpha_{\rm EM,0}=0.9993\pm 0.0025$ (CMB only) and $m_{\rm e}/m_{\rm e,0}= 1.0039 \pm 0.0074$ (including BAO) using Planck 2015 data. For a redshift-dependent variation, $\alpha_{\rm EM}(z)=\alpha_{\rm EM}(z_0)\,[(1+z)/1100]^p$ with $\alpha_{\rm EM}(z_0)\equiv\alpha_{\rm EM,0}$ at $z_0=1100$, we obtain $p=0.0008\pm 0.0025$. Allowing simultaneous variations of $\alpha_{\rm EM}(z_0)$ and $p$ yields $\alpha_{\rm EM}(z_0)/\alpha_{\rm EM,0} = 0.9998\pm 0.0036$ and $p = 0.0006\pm 0.0036$. We also discuss combined limits on $\alpha_{\rm EM}$ and $m_{\rm e}$. Our analysis shows that existing data is not only sensitive to the value of the fundamental constants around recombination but also its first time derivative. This suggests that a wider class of varying fundamental constant models can be probed using the CMB.

The precision of recent experiments such as ${\it Planck}$ have allowed us to constrain standard and non-standard physics (e.g., due to dark matter annihilation or varying fundamental constants) during the recombination epoch. However, we can also probe this era of cosmic history using model-independent variations of the free electron fraction, $X_{\rm e}$, which in turn affects the temperature and polarization anisotropies of the cosmic microwave background. In this paper, we improve on the previous efforts to construct and constrain these generalised perturbations in the ionization history, deriving new optimized eigenmodes based on the full Planck 2015 likelihood data, introducing the new module 'FEARec++'. We develop a direct likelihood sampling method for attaining the numerical derivatives of the standard and non-standard parameters, and discuss complications arising from the stability of the likelihood code. We improve the amplitude constraints of the ${\it Planck}$ 2015 principal components constructed here, $\mu_1=-0.09\pm0.12$, $\mu_2=-0.17\pm0.20$ and $\mu_3=-0.30\pm0.35$, finding no indication for departures from the standard recombination scenario. The third mode error is reduced by $60\%$. We utilise an efficient eigen-analyser that keeps the cross-correlations of the first three eigenmodes to $\xi_{\mu,\mu'}<0.1\%$ after marginalisation for all the considered data combinations. We also propose a new projection method for estimating constraints on the parameters of non-standard recombination. Using our eigenmodes, this allows us to recreate the ${\it Planck}$ constraint on the two-photon decay rate, $A_{\rm 2s1s}=7.60\pm0.64$, giving an estimate to within $\simeq 0.05\sigma$ of the full MCMC result. The improvements on the analysis using the ${\it Planck}$ data will allow us to implement this new method for analysis with fundamental constant variations in the future.

The cosmological recombination radiation (CRR) is one of the inevitable $\Lambda$CDM spectral distortions of the cosmic microwave background (CMB). While it shows a rich spectral structure across dm-mm wavelengths, it is also one of the smallest signals to target. Here we carry out a detailed forecast for the expected sensitivity levels required to not only detect but also extract cosmological information from the CRR in the presence of foregrounds. We use ${\tt CosmoSpec}$ to compute the CRR including all important radiative transfer effects and modifications to the recombination dynamics. We confirm that detections of the overall CRR signal are possible with spectrometer concepts like ${\it SuperPIXIE}$. However, for real exploitation of the cosmological information, a $\simeq 50$ times more sensitive spectrometer is required. While extremely futuristic, this could provide independent constraints on the primordial helium abundance, $Y_p$, and probe the presence of extra relativistic degrees of freedom during BBN and recombination. Significantly improving the constraints on other cosmological parameters requires even higher sensitivity (another factor of $\simeq 5$) when considering a combination of a CMB spectrometer with existing CMB data. To a large part this is due to astrophysical foregrounds which interestingly do not degrade the constraints on $Y_p$ and $N_{\rm eff}$ as much. A future CMB spectrometer could thus open a novel way of probing non-standard BBN scenarios, dark radiation and sterile neutrinos. In addition, inflation physics could be indirectly probed using the CRR in combination with existing and forthcoming CMB anisotropy data.

The method of weighted addition of multi-frequency maps, more commonly referred to as {\it Internal Linear Combination} (ILC), has been extensively employed in the measurement of cosmic microwave background (CMB) anisotropies and its secondaries along with similar application in 21cm data analysis. Here we argue and demonstrate that ILC methods can also be applied to data from absolutely-calibrated CMB experiments to extract average-sky signals in addition to the conventional CMB anisotropies. The performance of the simple ILC method is, however, limited, but can be significantly improved by adding constraints informed by physics and existing empirical information. In recent work, a moment description has been introduced as a technique of carrying out high precision modeling of foregrounds in the presence of inevitable averaging effects. We combine these two approaches to construct a heavily constrained form of the ILC, dubbed \milc, which can be used to recover tiny monopolar spectral distortion signals in the presence of realistic foregrounds and instrumental noise. This is a first demonstration for measurements of the monopolar and anisotropic spectral distortion signals using ILC and extended moment methods. We also show that CMB anisotropy measurements can be improved, reducing foreground biases and signal uncertainties when using the \milc. While here we focus on CMB spectral distortions, the scope extends to the 21cm monopole signal and $B$-mode analysis. We briefly discuss augmentations that need further study to reach the full potential of the method.

We present updated constraints on the variation of the fine structure constant, $\alpha_{\rm EM}$, and effective electron rest mass, $m_{\rm e}$, during the cosmological recombination era. These two fundamental constants directly affect the ionization history at redshift $z\simeq 1100$ and thus modify the temperature and polarisation anisotropies of the cosmic microwave background (CMB) measured precisely with {\it Planck }. The constraints on $\alpha_{\rm EM}$ tighten slightly due to improved {\it Planck} 2018 polarisation data but otherwise remain similar to previous CMB analysis. However, a comparison with the 2015 constraints reveals a mildly discordant behaviour for $m_{\rm e}$, which from CMB data alone is found below its local value. Adding baryon acoustic oscillation data brings $m_{\rm e}$ back to the fiducial value, $m_{\rm e}=(1.0078\pm0.0067) m_{\rm e,0}$, and also drives the Hubble parameter to $H_0=69.1\pm 1.2$ [in units of ${\rm km \, s^{-1} \, Mpc^{-1} }$]. Further adding supernova data yields $m_{\rm e}=(1.0190\pm0.0055) m_{\rm e,0}$ with $H_0=71.24\pm0.96$. We perform several comparative analyses using the latest cosmological recombination calculations to further understand the various effects. Our results indicate that a single-parameter extension allowing for a slightly increased value of $m_{\rm e}$ ($\simeq 3.5\sigma$ above $m_{\rm e,0}$) could play a role in the Hubble tension.

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.

The cosmic microwave background (CMB) monopole temperature evolves with the inverse of the cosmological scale factor, independent of many cosmological assumptions. With sufficient sensitivity, real-time cosmological observations could thus be used to measure the local expansion rate of the Universe using the cooling of the CMB. We forecast how well a CMB spectrometer could determine the Hubble constant via this method. The primary challenge of such a mission lies in the separation of Galactic and extra-Galactic foreground signals from the CMB at extremely high precision. However, overcoming these obstacles could potentially provide an independent, highly robust method to shed light on the current low-/high-$z$ Hubble tension. An experiment with 3000 linearly spaced bins between 5 GHz and 3 THz with a sensitivity of 1 $\mathrm{mJy\sqrt{yr}~sr^{-1}}$ per bin, could measure $H_0$ to 3% over a 10 year mission, given current foreground complexity. This sensitivity would also enable high-precision measurements of the expected $\Lambda$CDM spectral distortions, but remains futuristic at this stage.

With increasing sensitivity, angular resolution, and frequency coverage, future cosmic microwave background (CMB) experiments like PICO will allow us to access new information about galaxy clusters through the relativistic thermal Sunyaev-Zeldovich (SZ) effect. We will be able to map the temperature of relativistic electrons across the entire sky, going well beyond a simple detection of the relativistic SZ effect by cluster stacking methods that currently define the state-of-the-art. Here, we propose a new map-based approach utilizing SZ-temperature moment expansion and constrained-ILC methods to extract electron gas temperature maps from foreground-obscured CMB data. This delivers a new independent map-based observable, the electron temperature power spectrum $T_{\rm e}^{yy}(\ell)$, which can be used to constrain cosmology in addition to the Compton-$y$ power spectrum $C_\ell^{yy}(\ell)$ . We find that PICO has the required sensitivity, resolution, and frequency coverage to accurately map the electron gas temperature of galaxy clusters across the full sky, covering a broad range of angular scales. Frequency-coverage at $\nu\gtrsim 300\,{\rm GHz}$ plays an important role for extracting the relativistic SZ effect in the presence of foregrounds. For Coma, PICO will allow us to directly reconstruct the electron temperature profile using the relativistic SZ effect. Coma's average electron temperature will be measured to $10\sigma$ significance after foreground removal using PICO. Low-angular resolution CMB experiment like LiteBIRD could achieve $2\sigma$ to $3\sigma$ measurement of the electron temperature of this largest cluster. Our analysis highlights a new spectroscopic window into the thermodynamic properties of galaxy clusters and the diffuse electron gas at large angular scales.

The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (tSZ) effect. The third measure determines the lowest order relativistic correction to the tSZ signal, which is seeing increased observational relevance. Our procedure allows us to accurately model the relativistic SZ (rSZ) contribution and we show that a $\gtrsim 10\%-40\%$ underestimation of this rSZ cluster temperature is expected when applying standard X-ray relations. The correction also exhibits significant mass and redshift evolution, as we demonstrate here. We present the mass dependence of each temperature measure alongside their profiles and a short analysis of the temperature dispersion as derived from the aforementioned simulations. We also discuss a new relation connecting the temperature and Compton-$y$ parameter, which can be directly used for rSZ modelling. Simple fits to the obtained scaling relations and profiles are provided. These should be useful for future studies of the rSZ effect and its relevance to cluster cosmology.

Spectral distortions and anisotropies of the CMB provide independent and complementary probes to study energy injection processes in the early universe. Here we discuss the synergy between these observables, and show the promising future of spectral distortion missions to constrain both exotic and non-exotic energy injections. We show that conventional probes such as Big Bang Nucleosynthesis and CMB anisotropies can benefit from and even be surpassed by future spectral distortion experiments. For this, we have implemented a unified framework within the Boltzmann code CLASS to consistently treat the thermal evolution of photons and baryons. Furthermore, we give an extensive and pedagogical introduction into the topic of spectral distortions and energy injections throughout the thermal history of the universe, highlighting some of their unique features and potential as a novel probe for cosmology and particle physics.

In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field.

The characterization of the spectral energy distribution (SED) of dust emission has become a critical issue in the quest for primordial B-modes. The dust SED is often approximated by a modified blackbody (MBB) emission law but to what extent is it accurate? This paper addresses this question expanding the dust SED, at the power spectrum level, in moments around the MBB law, related to derivatives with respect to the dust spectral index. We present the mathematical formalism and apply it to simulations and Planck total intensity data, from 143 to 857 GHz, because today no polarized data provide the required sensitivity to perform this analysis. With simulations, we demonstrate the ability of high order moments to account for spatial variations in MBB parameters. Neglecting these moment leads to poor fits and a bias in the recovered dust spectral index. We identify the main moments that are required to fit the Planck data. The comparison with simulations help us to disentangle the respective contributions from dust and the Cosmic Infrared Background to the high order moments, but the simulations give an insufficient description of the actual Planck data. Extending our model to CMB B-mode analyses within a simplified framework, we find that ignoring the dust SED distortions, or trying to model them with a single decorrelation parameter, could lead to biases larger than the targeted sensitivity for the next generation of CMB B-mode experiments.

Constraints on dark matter annihilation or decay offer unique insights into the nature of dark matter. We illustrate how surveys dedicated to detect the highly redshifted 21 cm signal from the dark ages will offer a new window into properties of particle dark matter. The 21 cm intensity mapping signal and its fluctuations are sensitive to energy injection from annihilation or decay of long-lived particles in a way that is complementary to other probes. We present forecasted constraints from forthcoming and next-generation radio surveys. We show that, while SKA might be capable of a detection for some cases, the most promising opportunity to detect dark matter in the 21 cm intensity mapping signal is with a futuristic radio array on the lunar far-side, with the potential to detect a signal many orders of magnitude weaker than current or maximal constraints from other probes.

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.

Compton scattering between electrons and photons plays a crucial role in astrophysical plasmas. Many important aspects of this process can be captured by using the so-called Compton scattering kernel. For isotropic media, exact analytic expressions (valid at all electron and photon energies) do exist but are hampered by numerical issues and often are presented in complicated ways. In this paper, we summarize, simplify and improve existing analytic expressions for the Compton scattering kernel, with an eye on clarity and physical understanding. We provide a detailed overview of important properties of the kernel covering a wide range of energies and highlighting aspects that have not been appreciated as much previously. We discuss analytic expressions for the moments of the kernel, comparing various approximations and demonstrating their precision. We also illustrate the properties of the scattering kernel for thermal electrons at various temperatures and photon energies, introducing new analytic approximations valid to high temperatures. The obtained improved formulae for the kernel and its moments should prove useful in many astrophysical computations, one of them being the evolution of spectral distortions of the cosmic microwave background in the early Universe. A novel code, CSpack, for efficient computations of the Compton scattering kernel and its properties (in the future also including anisotropies in the initial electron and photon distributions) is being developed in a series of papers and will be available within one month.

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)

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.

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.

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

Departures of the energy spectrum of the cosmic microwave background (CMB) from a perfect blackbody probe a fundamental property of the universe -- its thermal history. Current upper limits, dating back some 25 years, limit such spectral distortions to 50 parts per million and provide a foundation for the Hot Big Bang model of the early universe. Modern upgrades to the 1980's-era technology behind these limits enable three orders of magnitude or greater improvement in sensitivity. The standard cosmological model provides compelling targets at this sensitivity, spanning cosmic history from the decay of primordial density perturbations to the role of baryonic feedback in structure formation. Fully utilizing this sensitivity requires concurrent improvements in our understanding of competing astrophysical foregrounds. We outline a program using proven technologies capable of detecting the minimal predicted distortions even for worst-case foreground scenarios.

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 consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich (SZ) power spectrum and cluster counts analyses, highlighting subtle dependencies on the parameterisation and data combination. In $\Lambda$CDM, with an X-ray mass bias corresponding the expected hydrostatic mass bias, i.e., $b\simeq0.2$, our constraints from Planck SZ data are consistent with the latest results from SPT, DES-Y1 and KiDS+VIKING-450. In $\nu\Lambda$CDM, without prior information on $b$, our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint, as well as $(1-b)=0.64\pm0.04$~(68\%~CL). For forecasts, we find that a combination of a mock cosmic variance limited SZ power spectrum with primary CMB and BAO can improve the uncertainty on the total neutrino mass by 14\% with respect to CMB combined with BAO. This requires masking the heaviest clusters and probing the small-scale SZ power spectrum up to multipoles of $\ell_\mathrm{max}=10^4$. For a future competitive measurement of the total neutrino mass using CMB and the SZ power spectrum, but excluding BAO and lensing power spectrum, we find that a 1\% precision on the mass calibration is needed. Although SZ power spectrum-based measurements of the neutrino masses are challenging, we find that the SZ power spectrum can be used to tightly constrain intra-cluster medium properties.

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.

Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

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 investigate the impact of a stochastic background of Primordial Magnetic Fields (PMF) generated before recombination on the ionization history of the Universe and on the Cosmic Microwave Background radiation (CMB). Pre-recombination PMFs are dissipated during recombination and reionization via decaying MHD turbulence and ambipolar diffusion. This modifies the local matter and electron temperatures and thus affects the ionization history and Thomson visibility function. We use this effect to constrain PMFs described by a spectrum of power-law type, extending our previous study (based on a scale-invariant spectrum) to arbitrary spectral index. We derive upper bounds on the integrated amplitude of PMFs due to the separate effect of ambipolar diffusion and MHD decaying turbulence and their combination. We show that ambipolar diffusion is relevant for $n_{\rm B}>0$ whereas for $n_{\rm B}<0$ MHD turbulence is more important. The bound marginalized over the spectral index on the integrated amplitude of PMFs with a sharp cut-off is $\sqrt{\langle B^2 \rangle}<0.83$ nG. We discuss the quantitative relevance of the assumptions on the damping mechanism and the comparison with previous bounds.

We present a detailed analysis of the stacked frequency spectrum of a large sample of galaxy clusters using Planck data, together with auxiliary data from the AKARI and IRAS missions. Our primary goal is to search for the imprint of relativistic corrections to the thermal Sunyaev-Zeldovich effect (tSZ) spectrum, which allow to measure the temperature of the intracluster medium. We remove Galactic and extragalactic foregrounds with a matched filtering technique, which is validated using simulations with realistic mock data sets. The extracted spectra show the tSZ signal at high significance and reveal an additional far-infrared (FIR) excess, which we attribute to thermal emission from the galaxy clusters themselves. This excess FIR emission from clusters is accounted for in our spectral model. We are able to measure the tSZ relativistic corrections at $2.2\sigma$ by constraining the mean temperature of our cluster sample to $4.4^{+2.1}_{-2.0} \, \mathrm{keV}$. We repeat the same analysis on a subsample containing only the 100 hottest clusters, for which we measure the mean temperature to be $6.0^{+3.8}_{-2.9} \, \mathrm{keV}$, corresponding to $2.0\sigma$. The temperature of the emitting dust grains in our FIR model is constrained to $\simeq 20 \, \mathrm{K}$, consistent with previous studies. Control for systematic biases is done by fitting mock clusters, from which we also show that using the non-relativistic spectrum for SZ signal extraction will lead to a bias in the integrated Compton parameter $Y$, which can be up to 14% for the most massive clusters. We conclude by providing an outlook for the upcoming CCAT-prime telescope, which will improve upon Planck with lower noise and better spatial resolution.

Measurements of the thermal Sunyaev-Zel'dovich (tSZ) effect have long been recognized as a powerful cosmological probe. Here we assess the importance of relativistic temperature corrections to the tSZ signal on the power spectrum analysis of the Planck Compton-$y$ map, developing a novel formalism to account for the associated effects. The amplitude of the tSZ power spectrum is found to be sensitive to the effective electron temperature, $\bar{T}_e$, of the cluster sample. Omitting the corresponding modifications leads to an underestimation of the $yy$-power spectrum amplitude. Relativistic corrections thus add to the error budget of tSZ power spectrum observables such as $\sigma_8$. This could help alleviate the tension between various cosmological probes, with the correction scaling as $\Delta \sigma_8/\sigma_8 \simeq 0.019\,[k\bar{T}_e\,/\,5\,{\rm keV}]$ for Planck. At the current level of precision, this implies a systematic shift by $\simeq 1\sigma$, which can also be interpreted as an overestimation of the hydrostatic mass bias by $\Delta b \simeq 0.046\,(1-b)\,[k\bar{T}_e\,/\,5\,{\rm keV}]$, bringing it into better agreement with hydrodynamical simulations. It is thus time to consider relativistic temperature corrections in the processing of current and future tSZ data.

Since the measurements of COBE/FIRAS in the mid-90's we know that the energy spectrum of the cosmic microwave background (CMB) is extremely close to that of a perfect blackbody at an average temperature T0~2.726K. However, a number of early-universe processes are expected to create CMB spectral distortions - departures of the average CMB energy spectrum from a blackbody - at a level that is within reach of present-day technology. This provides strong motivation to study the physics of CMB spectral distortions and ask what these small signals might be able to tell us about the Universe we live in. In this lecture, I will give a broad-brush overview of recent theoretical and experimental developments, explaining why future spectroscopic measurements of the CMB will open an unexplored new window to early-universe and particle physics. I will give an introduction about the different types of distortions, how they evolve and thermalize and highlight some of the physical processes that can cause them. I hope to be able to convince you that CMB spectral distortions could open an exciting new path forward in CMB cosmology, which is complementary to planned and ongoing searches for primordial B-mode polarization signals. Spectral distortions should thus be considered very seriously as part of the activities in the next decades.

The origin of cosmic magnetic fields is an unsolved problem and magnetogenesis could have occurred in the early Universe. We study the evolution of such primordial magnetic fields across the cosmological recombination epoch via 3D magnetohydrodynamic numerical simulations. We compute the effective or net heating rate of baryons due to decaying magnetic fields and its dependence on the magnetic field strength and spectral index. In the drag-dominated regime ($z \gtrsim 1500$), prior to recombination, we find no real heating is produced. Our simulations allow us to smoothly trace a new transition regime ($600 \lesssim z \lesssim 1500$), where magnetic energy decays, at first, into the kinetic energy of baryons. A turbulent velocity field is built up until it saturates, as the net heating rate rises from a low value at recombination to its peak towards the end of the transition regime. This is followed by a turbulent decay regime ($z \lesssim 600$) where magnetic energy dissipates via turbulent decay of both magnetic and velocity fields while net heating remains appreciable and declines slowly. Both the peak of the net heating rate and the onset of turbulent decay are delayed significantly beyond recombination, by up to 0.5 Myr (until $z\simeq 600-700$), for scale-invariant magnetic fields. We provide analytic approximations and present numerical results for a range of field strengths and spectral indices, illustrating the redshift-dependence of dissipation and net heating rates. These can be used to study cosmic microwave background constraints on primordial magnetic fields.

Correlations between cosmic microwave background (CMB) temperature, polarization and spectral distortion anisotropies can be used as a probe of primordial non-Gaussianity. Here, we perform a reconstruction of $\mu$-distortion anisotropies in the presence of Galactic and extragalactic foregrounds, applying the so-called Constrained ILC component separation method to simulations of proposed CMB space missions (PIXIE, LiteBIRD, CORE, PICO). Our sky simulations include Galactic dust, Galactic synchrotron, Galactic free-free, thermal Sunyaev-Zeldovich effect, as well as primary CMB temperature and $\mu$-distortion anisotropies, the latter being added as correlated field. The Constrained ILC method allows us to null the CMB temperature anisotropies in the reconstructed $\mu$-map (and vice versa), in addition to mitigating the contaminations from astrophysical foregrounds and instrumental noise. We compute the cross-power spectrum between the reconstructed (CMB-free) $\mu$-distortion map and the ($\mu$-free) CMB temperature map, after foreground removal and component separation. Since the cross-power spectrum is proportional to the primordial non-Gaussianity parameter, $f_{\rm NL}$, on scales $k\simeq 740$ Mpc$^{-1}$, this allows us to derive $f_{\rm NL}$-detection limits for the aforementioned future CMB experiments. Our analysis shows that foregrounds degrade the theoretical detection limits (based mostly on instrumental noise) by more than one order of magnitude, with PICO standing the best chance at placing upper limits on scale-dependent non-Gaussianity. We also discuss the dependence of the constraints on the channel sensitivities and chosen bands. Like for $B$-mode polarization measurements, extended coverage at frequencies $\nu\lesssim 40\,{\rm GHz}$ and $\nu\gtrsim 400\,{\rm GHz}$ provides more leverage than increased channel sensitivity.

One of the major goals of future Cosmic Microwave Background measurements is the accurate determination of the effective number of neutrinos $N_{\rm eff}$. Reaching an experimental sensitivity of $\Delta N_{\rm eff} = 0.013$ could indeed falsify the presence of any non-standard relativistic particles at $95 \%$ c.l.. In this paper, we test how this future constraint can be affected by the removal of two common assumptions: a negligible running of the inflationary spectral index $n_{\rm run}$ and a precise determination of the neutron lifetime $\tau_n$. We first show that the constraints on $N_{\rm eff}$ could be significantly biased by the unaccounted presence of a running of the spectral index. Considering the Stage-IV experiment, a negative running of ${\rm d}n/{\rm d}\ln k= - 0.002$ could mimic a positive variation of $\Delta N_{\rm eff} = 0.03$. Moreover, given the current discrepancies between experimental measurements of the neutron lifetime $\tau_n$, we show that the assumption of a conservative error of $\Delta\tau_n \sim 10$s could bring to a systematic error of $\Delta N_{\rm eff} = 0.02$. Complementary cosmological constraints on the running of the spectral index and a solution to the neutron lifetime discrepancy are therefore needed for an accurate and reliable future CMB bound of $N_{\rm eff}$ at percent level.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

Future high-sensitivity measurements of the cosmic microwave background (CMB) anisotropies and energy spectrum will be limited by our understanding and modeling of foregrounds. Not only does more information need to be gathered and combined, but also novel approaches for the modeling of foregrounds, commensurate with the vast improvements in sensitivity, have to be explored. Here, we study the inevitable effects of spatial averaging on the spectral shapes of typical foreground components, introducing a moment approach, which naturally extends the list of foreground parameters that have to be determined through measurements or constrained by theoretical models. Foregrounds are thought of as a superposition of individual emitting volume elements along the line of sight and across the sky, which then are observed through an instrumental beam. The beam and line of sight averages are inevitable. Instead of assuming a specific model for the distributions of physical parameters, our method identifies natural new spectral shapes for each foreground component that can be used to extract parameter moments (e.g., mean, dispersion, cross-terms, etc.). The method is illustrated for the superposition of power-laws, free-free spectra, gray-body and modified blackbody spectra, but can be applied to more complicated fundamental spectral energy distributions. Here, we focus on intensity signals but the method can be extended to the case of polarized emission. The averaging process automatically produces scale-dependent spectral shapes and the moment method can be used to propagate the required information across scales in power spectrum estimates. The approach is not limited to applications to CMB foregrounds but could also be useful for the modeling of X-ray emission in clusters of galaxies.

Measurements of cosmic microwave background spectral distortions have profound implications for our understanding of physical processes taking place over a vast window in cosmological history. Foreground contamination is unavoidable in such measurements and detailed signal-foreground separation will be necessary to extract cosmological science. We present MCMC-based spectral distortion detection forecasts in the presence of Galactic and extragalactic foregrounds for a range of possible experimental configurations, focusing on the Primordial Inflation Explorer (PIXIE) as a fiducial concept. We consider modifications to the baseline PIXIE mission (operating 12 months in distortion mode), searching for optimal configurations using a Fisher approach. Using only spectral information, we forecast an extended PIXIE mission to detect the expected average non-relativistic and relativistic thermal Sunyaev-Zeldovich distortions at high significance (194$\sigma$ and 11$\sigma$, respectively), even in the presence of foregrounds. The $\Lambda$CDM Silk damping $\mu$-type distortion is not detected without additional modifications of the instrument or external data. Galactic synchrotron radiation is the most problematic source of contamination in this respect, an issue that could be mitigated by combining PIXIE data with future ground-based observations at low frequencies ($\nu < 15-30$GHz). Assuming moderate external information on the synchrotron spectrum, we project an upper limit of $|\mu| < 3.6\times 10^{-7}$ (95\% c.l.), slightly more than one order of magnitude above the fiducial $\Lambda$CDM signal from the damping of small-scale primordial fluctuations, but a factor of $\simeq 250$ improvement over the current upper limit from COBE/FIRAS. This limit could be further reduced to $|\mu| < 9.4\times 10^{-8}$ (95\% c.l.) with more optimistic assumptions about low-frequency information. (Abridged)

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.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

The small-scale crisis, discrepancies between observations and N-body simulations, may imply suppressed matter fluctuations on subgalactic distance scales. Such a suppression could be caused by some early-universe mechanism (e.g., broken scale invariance during inflation), leading to a modification of the primordial power spectrum at the onset of the radiation-domination era. Alternatively, it may be due to nontrivial dark-matter properties (e.g., new dark-matter interactions or warm dark matter) that affect the matter power spectrum at late times, during radiation domination, after the perturbations re-enter the horizon. We show that early- and late-time suppression mechanisms can be distinguished by measurement of the $\mu$ distortion to the frequency spectrum of the cosmic microwave background. This is because the $\mu$ distortion is suppressed, if the power suppression is primordial, relative to the value expected from the dissipation of standard nearly scale-invariant fluctuations. We emphasize that the standard prediction of the $\mu$ distortion remains unchanged in late-time scenarios even if the dark-matter effects occur before or during the era (redshifts $5\times 10^4 \lesssim z \lesssim 2\times 10^6$) at which $\mu$ distortions are generated.

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.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

Cosmic baryon evolution during the Cosmic Dawn and Reionization results in redshifted 21-cm spectral distortions in the cosmic microwave background (CMB). These encode information about the nature and timing of first sources over redshifts 30-6 and appear at meter wavelengths as a tiny CMB distortion along with the Galactic and extragalactic radio sky, which is orders of magnitude brighter. Therefore, detection requires precise methods to model foregrounds. We present a method of foreground fitting using maximally smooth (MS) functions. We demonstrate the usefulness of MS functions over traditionally used polynomials to separate foregrounds from the Epoch of Reionization (EoR) signal. We also examine the level of spectral complexity in plausible foregrounds using GMOSS, a physically motivated model of the radio sky, and find that they are indeed smooth and can be modeled by MS functions to levels sufficient to discern the vanilla model of the EoR signal. We show that MS functions are loss resistant and robustly preserve EoR signal strength and turning points in the residuals. Finally, we demonstrate that in using a well-calibrated spectral radiometer and modeling foregrounds with MS functions, the global EoR signal can be detected with a Bayesian approach with 90\% confidence in 10 minutes' integration.

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

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 examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

We present Global MOdel for the radio Sky Spectrum (GMOSS) -- a novel, physically motivated model of the low-frequency radio sky from 22 MHz to 23 GHz. GMOSS invokes different physical components and associated radiative processes to describe the sky spectrum over 3072 pixels of $5^{\circ}$ resolution. The spectra are allowed to be convex, concave or of more complex form with contributions from synchrotron emission, thermal emission and free-free absorption included. Physical parameters that describe the model are optimized to best fit four all-sky maps at 150 MHz, 408 MHz, 1420 MHz and 23 GHz and two maps at 22 MHz and 45 MHz generated using the Global Sky Model of de Oliveira-Costa et al. (2008). The fractional deviation of model to data has a median value of $6\%$ and is less than $17\%$ for $99\%$ of the pixels. Though aimed at modeling of foregrounds for the global signal arising from the redshifted 21-cm line of Hydrogen during Cosmic Dawn and Epoch of Reionization (EoR) - over redshifts $150\lesssim z \lesssim 6$, GMOSS is well suited for any application that requires simulating spectra of the low-frequency radio sky as would be observed by the beam of any instrument. The complexity in spectral structure that naturally arises from the underlying physics of the model provides a useful expectation for departures from smoothness in EoR foreground spectra and hence may guide the development of algorithms for EoR signal detection. This aspect is further explored in a subsequent paper.

Anisotropies in distortions to the frequency spectrum of the cosmic microwave background (CMB) can be created through spatially varying heating processes in the early Universe. For instance, the dissipation of small-scale acoustic modes does create distortion anisotropies, in particular for non-Gaussian primordial perturbations. In this work, we derive approximations that allow describing the associated distortion field. We provide a systematic formulation of the problem using Fourier-space window functions, clarifying and generalizing previous approximations. Our expressions highlight the fact that the amplitudes of the spectral-distortion fluctuations induced by non-Gaussianity depend also on the homogeneous value of those distortions. Absolute measurements are thus required to obtain model-independent distortion constraints on primordial non-Gaussianity. We also include a simple description for the evolution of distortions through photon diffusion, showing that these corrections can usually be neglected. Our formulation provides a systematic framework for computing higher order correlation functions of distortions with CMB temperature anisotropies and can be extended to describe correlations with polarization anisotropies.

We present an extension to the general relativistic radiation magnetohydrodynamic code HARMRAD to account for emission and absorption by thermal cyclo-synchrotron, double Compton, bremsstrahlung, low-temperature OPAL opacities as well as Thomson and Compton scattering. We approximate the radiation field as a Bose-Einstein distribution and evolve it using the radiation number-energy-momentum conservation equations in order to track photon hardening. We perform various simulations to study how these extensions affect the radiative properties of magnetically-arrested disks accreting at Eddington to super-Eddington rates. We find that double Compton dominates bremsstrahlung in the disk within a radius of $r\sim 15r_g$ (gravitational radii) at a hundred times the Eddington accretion rate, and within smaller radii at lower accretion rates. Double Compton and cyclo-synchrotron regulate radiation and gas temperatures in the corona, while cyclo-synchrotron regulates temperatures in the jet. Interestingly, as the accretion rate drops to Eddington, an optically thin corona develops whose gas temperature of $T\sim 10^9$K is $\sim 100$ times higher than the disk's black body temperature. Our results show the importance of double Compton and synchrotron in super-Eddington disks, magnetized coronae, and jets.

We present direct estimates of the mean sky brightness temperature in observing bands around 99GHz and 242GHz due to line emission from distant galaxies. These values are calculated from the summed line emission observed in a blind, deep survey for specrtal line emission from high redshift galaxies using ALMA (the 'ASPECS' survey). In the 99 GHz band, the mean brightness will be dominated by rotational transitions of CO from intermediate and high redshift galaxies. In the 242GHz band, the emission could be a combination of higher order CO lines, and possibly [CII] 158$\mu$m line emission from very high redshift galaxies ($z \sim 6$ to 7). The mean line surface brightness is a quantity that is relevant to measurements of spectral distortions of the cosmic microwave background, and as a potential tool for studying large-scale structures in the early Universe using intensity mapping. While the cosmic volume and the number of detections are admittedly small, this pilot survey provides a direct measure of the mean line surface brightness, independent of conversion factors, excitation, or other galaxy formation model assumptions. The mean surface brightness in the 99GHZ band is: $T_B = 0.94\pm 0.09$ $\mu$K. In the 242GHz band, the mean brightness is: $T_B = 0.55\pm 0.033$ $\mu$K. These should be interpreted as lower limits on the average sky signal, since we only include lines detected individually in the blind survey, while in a low resolution intensity mapping experiment, there will also be the summed contribution from lower luminosity galaxies that cannot be detected individually in the current blind survey.

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

Ever refined cosmological measurements have established the $\Lambda$CDM concordance model, with the key cosmological parameters being determined to percent-level precision today. This allows us to make explicit predictions for the spectral distortions of the cosmic microwave background (CMB) created by various processes occurring in the early Universe. Here, we summarize all guaranteed CMB distortions and assess their total uncertainty within $\Lambda$CDM. We also compare simple methods for approximating them, highlighting some of the subtle aspects when it comes to interpreting future distortion measurements. Under simplified assumptions, we briefly study how well a PIXIE-like experiment may measure the main distortion parameters (i.e., $\mu$ and $y$). Next generation CMB spectrometers are expected to detect the distortion caused by reionization and structure formation at extremely high significance. They will also be able to constrain the small-scale power spectrum through the associated $\mu$-distortion, improving limits on running of the spectral index. Distortions from the recombination era, adiabatic cooling of matter relative to the CMB and dark matter annihilation require a higher sensitivity than PIXIE in its current design. The crucial next step is an improved modeling of foregrounds and instrumental aspects, as we briefly discuss here.

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

If local supersymmetry is the correct extension of the standard model of particle physics, then following Inflation the early universe would have been populated by gravitinos produced from scatterings in the hot plasma during reheating. Their abundance is directly related to the magnitude of the reheating temperature. The gravitino lifetime is fixed as a function of its mass, and for gravitinos with lifetimes longer than the age of the universe at redshift $z\simeq 2\times 10^{6}$ (or roughly $6\times 10^6{\rm s}$), decay products can produce spectral distortion of the cosmic microwave background. Currently available COBE/FIRAS limits on spectral distortion can, in certain cases, already be competitive with respect to cosmological constraints from primordial nucleosynthesis for some gravitino decay scenarios. We show how the sensitivity limits on $\mu$ and $y$ distortions that can be reached with current technology would improve constraints and possibly rule out a significant portion of the parameter space for gravitino masses and Inflation reheating temperatures.

The cosmic microwave background (CMB) energy spectrum is a near-perfect blackbody. The standard model of cosmology predicts small spectral distortions to this form, but no such distortion of the sky-averaged CMB spectrum has yet been measured. We calculate the largest expected distortion, which arises from the inverse Compton scattering of CMB photons off hot, free electrons, known as the thermal Sunyaev-Zel'dovich (tSZ) effect. We show that the predicted signal is roughly one order of magnitude below the current bound from the COBE-FIRAS experiment, but can be detected at enormous significance ($\gtrsim 1000\sigma$) by the proposed Primordial Inflation Explorer (PIXIE). Although cosmic variance reduces the effective signal-to-noise to $230\sigma$, this measurement will still yield a sub-percent constraint on the total thermal energy of electrons in the observable universe. Furthermore, we show that PIXIE can detect subtle relativistic effects in the sky-averaged tSZ signal at $30\sigma$, which directly probe moments of the optical depth-weighted intracluster medium electron temperature distribution. These effects break the degeneracy between the electron density and temperature in the mean tSZ signal, allowing a direct inference of the mean baryon density at low redshift. Future spectral distortion probes will thus determine the global thermodynamic properties of ionized gas in the universe with unprecedented precision. These measurements will impose a fundamental "integral constraint" on models of galaxy formation and the injection of feedback energy over cosmic time.

The energy spectrum of the cosmic microwave background (CMB) provides a powerful tool for constraining standard and non-standard physics in the primordial Universe. Previous studies mainly highlight spectral distortions (mu-, y- and r-type) created by episodes of early energy release; however, several processes also lead to copious photon production, which requires a different treatment. Here, we carry out a first detailed study for the evolution of distortions caused by photon injection at different energies in the CMB bands. We provide detailed analytical and numerical calculations illustrating the rich phenomenology of the associated distortion signals. We show that photon injection at very high and very low frequencies creates distortions that are similar to those from pure energy release. In the mu-era (z>3x10^5), a positive or negative chemical potential can be formed, depending on the balance between added photon energy and number. At lower redshifts (z<3x10^5), partial information about the photon injection process (i.e., injection time and energy) can still be recovered, with the distortion being found in a partially comptonized state. We briefly discuss current and future constraints on scenarios with photon production. We also argue that more detailed calculations for different scenarios with photon injection may be required to assess in which regimes these can be distinguished from pure energy release processes.

We present the first fast and detailed computation of the cosmological recombination radiation released during the hydrogen (redshift z ~ 1300) and helium (z ~ 2500 and z ~ 6000) recombination epochs, introducing the code CosmoSpec. Our computations include important radiative transfer effects, 500-shell bound-bound and free-bound emission for all three species, the effects of electron scattering and free-free absorption as well as interspecies (HeII --> HeI --> HI) photon feedback. The latter effect modifies the shape and amplitude of the recombination radiation and CosmoSpec improves significantly over previous treatments of it. Utilizing effective multilevel atom and conductance approaches, one calculation takes only ~ 15 seconds on a standard laptop as opposed to days for previous computations. This is an important step towards detailed forecasts and feasibility studies considering the detection of the cosmological recombination lines and what one may hope to learn from the ~ 6.1 photons emitted per hydrogen atom in the three recombination eras. We briefly illustrate some of the parameter dependencies and discuss remaining uncertainties in particular related to collisional processes and the neutral helium atom model.

We consider the distortions of the cosmic microwave background (CMB) dipole anisotropy related to the primordial recombination radiation (PRR) and primordial $y$- and $\mu$-distortions. The signals arise due to our motion relative to the CMB restframe and appear as a frequency-dependent distortion of the CMB temperature dipole. To leading order, the expected relative distortion of the CMB dipole does not depend on the particular observation directions and reaches the level of $10^{-6}$ for the PRR- and $\mu$-distortions and $10^{-5}$ for the $y$-distortion in the frequency range 1 -- 700 GHz. The temperature differences arising from the dipole anisotropy of the relic CMB distortions depend on the observation directions. For mutually opposite directions, collinear to the CMB dipole axis, the temperature differences because of the PRR- and $\mu$-dipole anisotropy attain values $\Delta T\simeq 10\,$nK in the considered range. The temperature difference arising from the $y$-dipole anisotropy may reach values of up to $1\,\mu$K. The key features of the considered effect are as follow: (i) an observation of the effect does not require absolute calibration; (ii) patches of sky with minimal foreground contamination can be chosen. Future measurements of the CMB dipole distortion thus will provide an alternative method for direct detection of the PRR-, $y$-, and $\mu$-distortions. The $y$-distortion dipole may be detectable with PIXIE at a few standard deviations.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

We propose a new method to constrain elastic scattering between dark matter (DM) and standard model particles in the early Universe. Direct or indirect thermal coupling of non-relativistic DM with photons leads to a heat sink for the latter. This results in spectral distortions of the cosmic microwave background (CMB), the amplitude of which can be as large as a few times the DM-to-photon number ratio. We compute CMB spectral distortions due to DM-proton, DM-electron and DM-photon scattering for generic energy-dependent cross sections and DM mass m_DM >~ 1 keV. Using FIRAS measurements we set constraints on the cross sections for m_DM <~ 0.1 MeV. In particular, for energy-independent scattering we obtain sigma[DM-proton] <~ 10^(-24) cm^2 (keV/m_DM)^(1/2), sigma[DM-electron] <~ 10^(-27) cm^2 (keV/m_DM)^(1/2) and sigma[DM-photon] <~ 10^(-39) cm^2 (m_DM/keV). An experiment with the characteristics of PIXIE would extend the regime of sensitivity up to masses m_DM ~ 1 GeV.

All cosmological observations to date are consistent with adiabatic, Gaussian and nearly scale invariant initial conditions. These findings provide strong evidence for a particular symmetry breaking pattern in the very early universe (with a close to vanishing order parameter, $\epsilon$), widely accepted as conforming to the predictions of the simplest realizations of the inflationary paradigm. However, given that our observations are only privy to perturbations, in inferring something about the background that gave rise to them, it should be clear that many different underlying constructions project onto the same set of cosmological observables. Features in the primordial correlation functions, if present, would offer a unique and discriminating window onto the parent theory in which the mechanism that generated the initial conditions is embedded. In certain contexts, simple linear response theory allows us to infer new characteristic scales from the presence of features that can break the aforementioned degeneracies among different background models, and in some cases can even offer a limited spectroscopy of the heavier degrees of freedom that couple to the inflaton. In this review, we offer a pedagogical survey of the diverse, theoretically well grounded mechanisms which can imprint features into primordial correlation functions in addition to reviewing the techniques one can employ to probe observations. These observations include cosmic microwave background anisotropies and spectral distortions as well as the matter two and three point functions as inferred from large-scale structure and potentially, 21 cm surveys.

Many inflation models predict that primordial density perturbations have a nonzero three-point correlation function, or bispectrum in Fourier space. Of the several possibilities for this bispectrum, the most commmon is the local-model bispectrum, which can be described as a spatial modulation of the small-scale (large-wavenumber) power spectrum by long-wavelength density fluctuations. While the local model predicts this spatial modulation to be scale-independent, many variants have some scale-dependence. Here we note that this scale dependence can be probed with measurements of frequency-spectrum distortions in the cosmic microwave background (CMB), in particular highlighting Compton-$y$ distortions. Dissipation of primordial perturbations with wavenumbers $50\,{\rm Mpc}^{-1} \lesssim k \lesssim 10^4\,{\rm Mpc}^{-1}$ give rise to chemical-potential ($\mu$) distortions, while those with wavenumbers $1\,{\rm Mpc}^{-1} \lesssim k \lesssim 50\,{\rm Mpc}^{-1}$ give rise to Compton-$y$ distortions. With local-model non-Gaussianity, the distortions induced by this dissipation can be distinguished from those due to other sources via their cross-correlation with the CMB temperature $T$. We show that the relative strengths of the $\mu T$ and $yT$ correlations thus probe the scale-dependence of non-Gaussianity and estimate the magnitude of possible signals relative to sensitivities of future experiments. We discuss the complementarity of these measurements with other probes of squeezed-limit non-Gaussianity.

Spectral distortions of the CMB have recently experienced an increased interest. One of the inevitable distortion signals of our cosmological concordance model is created by the cosmological recombination process, just a little before photons last scatter at redshift $z\simeq 1100$. These cosmological recombination lines, emitted by the hydrogen and helium plasma, should still be observable as tiny deviation from the CMB blackbody spectrum in the cm--dm spectral bands. In this paper, we present a forecast for the detectability of the recombination signal with future satellite experiments. We argue that serious consideration for future CMB experiments in space should be given to probing spectral distortions and, in particular, the recombination line signals. The cosmological recombination radiation not only allows determination of standard cosmological parameters, but also provides a direct observational confirmation for one of the key ingredients of our cosmological model: the cosmological recombination history. We show that, with present technology, such experiments are futuristic but feasible. The potential rewards won by opening this new window to the very early universe could be considerable.

Primordial magnetic fields (PMF) damp at scales smaller than the photon diffusion and free-streaming scale. This leads to heating of ordinary matter (electrons and baryons), which affects both the thermal and ionization history of our Universe. Here, we study the effect of heating due to ambipolar diffusion and decaying magnetic turbulence. We find that changes to the ionization history computed with recfast are significantly overestimated when compared with CosmoRec. The main physical reason for the difference is that the photoionization coefficient has to be evaluated using the radiation temperature rather than the matter temperature. A good agreement with CosmoRec is found after changing this aspect. Using Planck 2013 data and considering only the effect of PMF-induced heating, we find an upper limit on the r.m.s. magnetic field amplitude of B0 < 1.1 nG (95% c.l.) for a stochastic background of PMF with a nearly scale-invariant power spectrum. We also discuss uncertainties related to the approximations for the heating rates and differences with respect to previous studies. Our results are important for the derivation of constraints on the PMF power spectrum obtained from measurements of the cosmic microwave background anisotropies with full-mission Planck data. They may also change some of the calculations of PMF-induced effects on the primordial chemistry and 21cm signals.

Photons emitted during the epochs of Hydrogen ($500 \lesssim z \lesssim 1600$) and Helium recombination ($1600 \lesssim z \lesssim 3500$ for HeII $\rightarrow$ HeI, $5000 \lesssim z \lesssim 8000$ for HeIII $\rightarrow$ HeII) are predicted to appear as broad, weak spectral distortions of the Cosmic Microwave Background. We present a feasibility study for a ground-based experimental detection of these recombination lines, which would provide an observational constraint on the thermal ionization history of the Universe, uniquely probing astrophysical cosmology beyond the last scattering surface. We find that an octave band in the 2--6 GHz window is optimal for such an experiment, both maximizing signal-to-noise ratio and including sufficient line spectral structure. At these frequencies the predicted signal appears as an additive quasi-sinusoidal component with amplitude about $8$ nK that is embedded in a sky spectrum some nine orders of magnitude brighter. We discuss an algorithm to detect these tiny spectral fluctuations in the sky spectrum by foreground modeling. We introduce a \textit{Maximally Smooth} function capable of describing the foreground spectrum and distinguishing the signal of interest. With Bayesian statistical tests and mock data we estimate that a detection of the predicted distortions is possible with 90\% confidence by observing for 255 days with an array of 128 radiometers using cryogenically cooled state-of-the-art receivers. We conclude that detection is in principle feasible in realistic observing times; we propose APSERa---Array of Precision Spectrometers for the Epoch of Recombination---a dedicated radio telescope to detect these recombination lines.

Spectral distortions of the cosmic microwave background (CMB) may become a powerful probe of primordial perturbations at small scales. Existing studies of spectral distortions focus almost exclusively on primordial scalar metric perturbations. Similarly, vector and tensor perturbations should source CMB spectral distortions. In this paper, we give general expressions for the effective heating rate caused by these types of perturbations, including previously neglected contributions from polarization states and higher multipoles. We then focus our discussion on the dissipation of tensors, showing that for nearly scale invariant tensor power spectra, the overall distortion is some six orders of magnitudes smaller than from the damping of adiabatic scalar modes. We find simple analytic expressions describing the effective heating rate from tensors using a quasi-tight coupling approximation. In contrast to adiabatic modes, tensors cause heating without additional photon diffusion and thus over a wider range of scales, as recently pointed out by Ota et. al 2014. Our results are in broad agreement with their conclusions, but we find that small-scale modes beyond k< 2x10^4 Mpc^{-1} cannot be neglected, leading to a larger distortion, especially for very blue tensor power spectra. At small scales, also the effect of neutrino damping on the tensor amplitude needs to be included.

In the expanding Universe, the average temperature of the cosmic microwave background (CMB) is expected to depend like TCMB~(1+z) on redshift z. Adiabatic photon production (or destruction) or deviations from isotropy and homogeneity could modify this scaling and several observational tests have been carried out in response. Here, we explain why `adiabatic' conditions are extremely difficult to establish in the redshift range targeted by these tests. Thus, instead of leading to a simple rescaling of the CMB temperature, a spectral distortion should be produced, which can be constrained using COBE/FIRAS. For scenarios with late photon production, tests of the temperature-redshift relation (TRR) should therefore be reinterpreted as weak spectral distortion limits, directly probing the energy dependence of the photon production process. For inhomogeneous cosmologies, an average y-type distortion is produced, but this type of distortion can be created in several other ways. Here, we briefly discuss possible effects that may help disentangling different contributions to the distortion signal, finding this to be very challenging. We furthermore argue that tests of the TRR, using the SZ effect have limited applicability and that for non-gravitational changes to the TRR the CMB anisotropy spectrum should exhibit an additional y-type dependence.

Aberration kernels describe how harmonic-space multipole coefficients of cosmic microwave background (CMB) observables transform under Lorentz boosts of the reference frame. For spin-weighted CMB observables, transforming like the CMB temperature (i.e. Doppler weight $d = 1$), we show that the aberration kernels are the matrix elements of a unitary boost operator in harmonic space. Algebraic properties of the rotation and boost generators then give simple, exact recursion relations that allow us to raise or lower the multipole quantum numbers $\ell$ and $m$, and the spin weight $s$. Further recursion relations express kernels of other Doppler weights $d \neq 1$ in terms of the $d = 1$ kernels. From those we show that on the full sky, to all orders in $\beta=v/c$, $E$- and $B$-mode polarization observables do not mix under aberration if and only if $d = 1$. The new relations, fully non-linear in the boost velocity $\beta$, form the basis of a practical recursive algorithm to accurately compute any aberration kernel. In addition, we develop a second, fast algorithm in which aberration kernels are obtained using a set of ordinary differential equations. This system can also be approximately solved at small scales, providing simple asymptotic formulae for the aberration kernels. The results of this work will be useful for further studying the effect of aberration on future CMB temperature and polarization analysis, and might provide a basis for relativistic radiative transfer schemes.

The measurements of COBE/FIRAS have shown that the CMB spectrum is extremely close to a perfect blackbody. There are, however, a number of processes in the early Universe that should create spectral distortions at a level which is within reach of present day technology. In this talk, I will give a brief overview of recent theoretical and experimental developments, explaining why future measurements of the CMB spectrum will open up an unexplored window to early-universe and particle physics with possible non-standard surprises but also several guaranteed signals awaiting us.

For the simplest inflation models, the recent detection of a large primordial B-mode polarization signal by the BICEP2 experiment indicates a slight tension with the upper limit on the tensor-to-scalar ratio, $r$, from the Planck satellite. Here, we discuss spatially varying $r$ as a possible explanation for this discrepancy. This idea seems attractive since it may also explain part of the hemispherical temperature power asymmetry seen by WMAP and Planck at large angular scales. If these two aspects are indeed connected, the model suggests that in the Northern hemisphere $r$ should be much smaller, a hypothesis that could be confirmed with future B-mode experiments, providing a test for the stationarity of primordial tensor contributions across the sky. The BICEP2 measurement furthermore rules out that a simple dipolar modulation of $r$ alone can be responsible for the full hemispherical power asymmetry.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

The physical processes affecting the thermalization of cosmic microwave background spectral distortions are very simple and well understood. This allows us to make precise predictions for the distortions signals caused by various energy release scenarios, where the theoretical uncertainty is largely dominated by the physical ingredients that are used for the calculation. Here, we revisit various approximations for the distortion visibility function -- defined using the fraction of the released energy that does not thermalize -- and early $\mu$-type distortions. Our approach is based on a perturbative expansion, which allows us to identify and clarify the origin of different improvements over earlier approximations. It provides a better than ~0.1%-1% description of our numerical results over a wide range of parameters. In particular, we are able to capture the high-frequency part of the mu-distortion, which directly depends on the time derivative of the electron temperature. We also include lowest order double Compton and Compton scattering relativistic corrections, finding that because of cancelation they increase the thermalization efficiency in the tail of the distortion visibility function by only ~10% (at $z\lesssim 6\times10^6$), although individually their effect can reach ~20%-40%.

We study the dissipation of small-scale adiabatic perturbations at early times when the Universe is hotter than T ~ 0.5 keV. When the wavelength falls below the damping scale 1/kD, the acoustic modes diffuse and thermalize, causing entropy production. Before neutrino decoupling, kD is primarily set by the neutrino shear viscosity, and we study the effect of acoustic damping on the relic neutrino number, primordial nucleosynthesis, dark-matter freeze-out, and baryogenesis. This sets a new limit on the amplitude of primordial fluctuations of DeltaR^2 < 0.007 at 10^4/Mpc< k < 10^5/Mpc and a model dependent limit of DeltaR^2 < 0.3 at k < 10^{20-25}/Mpc.

Future high-resolution, high-sensitivity Sunyaev-Zeldovich (SZ) observations of individual clusters will allow us to study the dynamical state of the intra cluster medium (ICM). To exploit the full potential of this new observational window, it is crucial to understand the different contributions to the total SZ signal. Here, we investigate the signature caused by multiple scatterings at lowest order of the electron temperature. Previous analytic discussions of this problem used the isotropic scattering approximation (ISA), which even for the simplest cluster geometries is rather rough. We take a step forward and consistently treat the anisotropy of the ambient radiation field caused by the first scattering. We show that the multiple scattering SZ signal directly probes line of sight anisotropies of the ICM, thereby delivering a new set of SZ observables which could be used for 3D cluster-profile reconstruction. The multiple scattering signal should furthermore correlate spatially with the cluster's X-ray and SZ polarization signals, an effect that could allow enhancing the detectability of this contribution.

Departures of the Cosmic Microwave Background (CMB) frequency spectrum from a blackbody - commonly referred to as spectral distortions - encode information about the thermal history of the early Universe (redshift z < few x 10^6). While the signal is usually characterized as mu- and y-type distortion, a smaller residual (non-y/non-mu) distortion can also be created at intermediate redshifts 10^4 < z < 3 x 10^5. Here, we construct a new set of observables, mu_k, that describes the principal components of this residual distortion. The principal components are orthogonal to temperature shift, y- and mu-type distortion, and ranked by their detectability, thereby delivering a compression of all valuable information offered by the CMB spectrum. This method provides an efficient way of analyzing the spectral distortion for given experimental settings, and can be applied to a wide range of energy-release scenarios. As an illustration, we discuss the analysis of the spectral distortion signatures caused by dissipation of small-scale acoustic waves and decaying/annihilating particles for a PIXIE-type experiments. We provide forecasts for the expected measurement uncertainties of model parameters and detections limits in each case. We furthermore show that a PIXIE-type experiments can in principle distinguish dissipative energy release from particle decays for a nearly scale-invariant primordial power spectrum with small running. Future CMB spectroscopy thus offers a unique probe of physical processes in the primordial Universe.

We study the multiple scattering Sunyaev-Zeldovich (SZ) signature, extending our previous analysis to high-temperature clusters. We consistently treat the anisotropy of the ambient radiation field caused by the first scattering and also consider lowest order kinematic terms. We show that due to temperature corrections monopole through octupole anisotropy of the singly scattered SZ signal attain different spectra in the second scattering. The difference becomes more pronounced at high temperature, and thus could be used to constrain individual line of sight moments of the electron density and temperature profiles. While very challenging from the observational point of view, this further extends the list of possible SZ observables that will be important for 3D cluster-profile reconstruction, possibly helping to break geometric degeneracies caused by projection effects. We also briefly discuss the scattering of primordial CMB anisotropies by SZ clusters.

Deviations of the cosmic microwave background (CMB) frequency spectrum from a pure blackbody tell an exciting story about the thermal history of our Universe. In this paper we illustrate how well future CMB measurements might decipher this tale, envisioning a PIXIE-like spectrometer, which could improve the distortion constraints obtained with COBE/FIRAS some 20 years ago by at least three orders of magnitude. This opens a large discovery space, offering deep insights to particle and early-universe physics, opportunities that no longer should be left unexplored. Specifically, we consider scenarios with annihilating and decaying relic particles, as well as signatures from the dissipation of primordial small-scale power. PIXIE can potentially rule out different early-universe scenarios, and moreover will allow unambiguous detections in many of the considered cases, as we demonstrate here. We also discuss slightly more futuristic experiments, with several times improved sensitivities, to highlight the large potential of this new window to the pre-recombination universe.

Our motion relative to the cosmic-microwave-background (CMB) rest frame deflects light rays giving rise to shifts as large as L -> L(1+beta), where beta=0.00123 is our velocity (in units of the speed of light) on measurements CMB fluctuations. Here we present a novel harmonic-space approach to this CMB aberration that improves upon prior work by allowing us to (i) go to higher orders in beta, thus extending the validity of the analysis to measurements at L > 1/beta ~ 800; and (ii) treat the effects of window functions and pixelization in a more accurate and computationally efficient manner. We calculate precisely the magnitude of the systematic bias in the power spectrum inferred from the partial sky, and show that aberration shifts the multipole moment by Delta L/L ~ beta<cos(theta)>, with <cos(theta)> averaged over the survey footprint. Such a shift, if ignored, would bias the measurement of the sound-horizon size theta_* at the 0.01%-level, which is comparable to the measurement uncertainties of Planck. The bias can then propagate into cosmological parameters such as the angular-diameter distance, Hubble parameter and dark-energy equation of state. We study the effect of aberration for current Planck, South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) data and show that the bias cannot be neglected. We suggest that the small tension between Planck, ACT, and SPT may be due partially to aberration. An Appendix shows how the near constancy of the full-sky power spectrum under aberration follows from unitarity of the aberration kernel.

The damping of primordial perturbations at small scales gives rise to distortions of the cosmic microwave background (CMB). Here, the dependence of the distortion on the different types of cosmological initial conditions is explored, covering adiabatic, baryon/cold dark matter isocurvature, neutrino density/velocity isocurvature modes and some mixtures. The radiation transfer functions for each mode are determined and then used to compute the dissipative heating rates and spectral distortion signatures, utilizing both analytic estimates and numerical results from the thermalization code CosmoTherm. Along the way, the early-time super-horizon behavior for the resulting fluid modes is derived in conformal Newtonian gauge, and tight-coupling transfer function approximations are given. CMB spectral distortions caused by different perturbation modes can be estimated using simple k-space window functions which are provided here. Neutrinos carry away some fraction of the primordial perturbation power, introducing an overall efficiency factor that depends on the perturbation type. It is shown that future measurements of the CMB frequency spectrum have the potential to probe different perturbation modes at very small scales (corresponding to wavenumbers 1 Mpc^{-1} < k < few x 10^4 Mpc^{-1}). These constraints are complementary to those obtained at large scales and hence provide an exciting new window to early-universe physics.

Energy release in the early Universe leads to spectral distortions of the cosmic microwave background (CMB) which in the future might allow probing different physical processes in the pre-recombination (z>~10^3) epoch. Depending on the energy injection history, the associated distortion partially thermalizes due to the combined action of Compton scattering, double Compton scattering and Bremsstrahlung emission, a problem that in general is hard to solve. Various analytic approximations describing the resulting distortion exist, however, for small distortions and fixed background cosmology the Green's function of the problem can be pre-computed numerically. Here we show that this approach gives very accurate results for a wide range of thermal histories, allowing fast and quasi-exact computation of the spectral distortion given the energy release rate. Our method is thus useful for forecasts of possible constraints on early-universe physics obtained from future measurements of the CMB spectrum.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

The cosmological recombination of H and He at z \sim 1000 and the formation of H2 during the dark ages produce a non-thermal photon excess in the Wien tail of the cosmic microwave background (CMB) blackbody spectrum. Here we compute the effect of these photons on the H- photodetachment and H2+ photodissociation processes. We discuss the implications for the chemical evolution of the Universe in the post-recombination epoch, emphasizing how important a detailed account of the full vibrational manifold of H2 and H2+ in the chemical network is. We find that the final abundances of H2, H2+, H3+ and HD are significantly smaller than in previous calculations that neglected the effect of non-thermal photons. The suppression is mainly caused by extra hydrogen recombination photons and could affect the formation rate of first stars. We provide simple analytical approximations for the relevant rate coefficients and briefly discuss the additional effect of dark matter annihilation on the considered reaction rates.

Physical models for the hemispherical power asymmetry in the cosmic microwave background (CMB) reported by the Planck Collaboration must satisfy CMB constraints to the homogeneity of the Universe and quasar constraints to power asymmetries. We survey a variety of models for the power asymmetry and show that consistent models include a modulated scale-dependent isocurvature contribution to the matter power spectrum or a modulation of the reionization optical depth, gravitational-wave amplitude, or scalar spectral index. We propose further tests to distinguish between the different scenarios.

Recent data from the WMAP, ACT and SPT experiments provide precise measurements of the cosmic microwave background temperature power spectrum over a wide range of angular scales. The combination of these observations is well fit by the standard, spatially flat LCDM cosmological model, constraining six free parameters to within a few percent. The scalar spectral index, n_s = 0.9690 +/- 0.0089, is less than unity at the 3.6 sigma level, consistent with simple models of inflation. The damping tail of the power spectrum at high resolution, combined with the amplitude of gravitational lensing measured by ACT and SPT, constrains the effective number of relativistic species to be N_eff = 3.28 +/- 0.40, in agreement with the standard model's three species of light neutrinos.

Future high resolution, high sensitivity Sunyaev-Zeldovich (SZ) observations of individual clusters will provide an exciting opportunity to answer specific questions about the dynamical state of the intra-cluster medium (ICM). In this paper we develop a new method that clearly shows the connection of the SZ signal with the underlying cluster model. We include relativistic temperature and kinematic corrections in the single-scattering approximation, allowing studies of hot clusters. In our approach, particular moments of the temperature and velocity field along the line-of-sight determine the precise spectral shape and morphology of the SZ signal. We illustrate how to apply our method to different cluster models, highlighting parameter degeneracies and instrumental effects that are important for interpreting future high-resolution SZ data. Our analysis shows that line-of-sight temperature variations can introduce significant biases in the derived SZ temperature and peculiar velocity. We furthermore discuss how the position of the SZ null is affected by the cluster's temperature and velocity structure. Our computations indicate that the SZ signal around the null alone is rather insensitive to different cluster models and that high frequency channels add a large leverage in this respect. We also apply our method to recent high sensitivity SZ data of the Bullet cluster, showing how the results can be linked to line-of-sight variations in the electron temperature. The tools developed here as part of SZpack should be useful for analyzing high-resolution SZ data and computing SZ maps from simulated clusters.

The primordial CMB at small angular scales is sensitive to the ionization and expansion history of the universe around the time of recombination. This dependence has been exploited to constrain the helium abundance and the effective number of relativistic species. Here we focus on allowed ionization fraction trajectories, $\Xe (z)$, by constraining low-order principal components of perturbations to the standard recombination scenario ($\Xe$-eigenmodes) in the circa 2011 SPT, ACT and WMAP7 data. Although the trajectories are statistically consistent with the standard recombination, we find that there is a tension similar to that found by varying the helium fraction. As this paper was in press, final SPT and ACT datasets were released and we applied our framework to them: we find the tension continues, with slightly higher significance, in the new 2012 SPT data, but find no tension with the standard model of recombination in the new 2012 ACT data. We find that the prior probabilities on the eigenamplitudes are substantially influenced by the requirement that $\Xe$ trajectories conserve electron number. We propose requiring a sufficient entropy decrease between posterior and prior marginalized distributions be used as an $\Xe$-mode selection criterion. We find that in the case of the 2011 SPT/ACT+WMAP7 data only two modes are constrainable, but upcoming ACTPol, Planck and SPTPol data will be able to test more modes and more precisely address the current tension.

New generation ground and space-based CMB experiments have ushered in discoveries of massive galaxy clusters via the Sunyaev-Zeldovich (SZ) effect, providing a new window for studying cluster astrophysics and cosmology. Many of the newly discovered, SZ-selected clusters contain hot intracluster plasma (kTe > 10 keV) and exhibit disturbed morphology, indicative of frequent mergers with large peculiar velocity (v > 1000 km s^{-1}). It is well-known that for the interpretation of the SZ signal from hot, moving galaxy clusters, relativistic corrections must be taken into account, and in this work, we present a fast and accurate method for computing these effects. Our approach is based on an alternative derivation of the Boltzmann collision term which provides new physical insight into the sources of different kinematic corrections in the scattering problem. This allows us to obtain a clean separation of kinematic and scattering terms which differs from previous works. We briefly mention additional complications connected with kinematic effects that should be considered when interpreting future SZ data for individual clusters. One of the main outcomes of this work is SZpack, a numerical library which allows very fast and precise (<~0.001% at frequencies h nu <~ 20kT_g) computation of the SZ signals up to high electron temperature (kT_e ~ 25 keV) and large peculiar velocity (v/c ~ 0.01). The accuracy is well beyond the current and future precision of SZ observations and practically eliminates uncertainties related to more expensive numerical evaluation of the Boltzmann collision term. Our new approach should therefore be useful for analyzing future high-resolution, multi-frequency SZ observations as well as computing the predicted SZ effect signals from numerical simulations.

Mixing of blackbodies with different temperatures creates a spectral distortion which, at lowest order, is a y-type distortion, indistinguishable from the thermal y-type distortion produced by the scattering of CMB photons by hot electrons residing in clusters of galaxies. This process occurs in the radiation-pressure dominated early Universe, when the primordial perturbations excite standing sound waves on entering the sound horizon. Photons from different phases of the sound waves, having different temperatures, diffuse through the electron-baryon plasma and mix together. This diffusion, with the length defined by Thomson scattering, dissipates sound waves and creates spectral distortions in the CMB. Of the total dissipated energy, 2/3 raises the average temperature of the blackbody part of spectrum, while 1/3 creates a distortion of y-type. It is well known that at redshifts 10^5< z< 2x10^6, comptonization rapidly transforms y-distortions into a Bose-Einstein spectrum. The chemical potential of the Bose-Einstein spectrum is again 1/3 the value we would get if all the dissipated energy was injected into a blackbody spectrum but no extra photons were added. We study the mixing of blackbody spectra, emphasizing the thermodynamic point of view, and identifying spectral distortions with entropy creation. This allows us to obtain the main results connected with the dissipation of sound waves in the early Universe in a very simple way. We also show that mixing of blackbodies in general, and dissipation of sound waves in particular, leads to creation of entropy.

Silk damping of primordial small-scale perturbations in the photon-baryon fluid due to diffusion of photons inevitably creates spectral distortions in the CMB. With the proposed CMB experiment PIXIE it might become possible to measure these distortions and thereby constrain the primordial power spectrum at comoving wavenumbers 50 Mpc^{-1} < k < 10^4 Mpc^{-1}. Since primordial fluctuations in the CMB on these scales are completely erased by Silk damping, these distortions may provide the only way to shed light on otherwise unobservable aspects of inflationary physics. A consistent treatment of the primordial dissipation problem requires going to second order in perturbation theory, while thermalization of these distortions necessitates consideration of second order in Compton scattering energy transfer. Here we give a full 2x2 treatment for the creation and evolution of spectral distortions due to the acoustic dissipation process, consistently including the effect of polarization and photon mixing in the free streaming regime. We show that 1/3 of the total energy (9/4 larger than previous estimates) stored in small-scale temperature perturbations imprints observable spectral distortions, while the remaining 2/3 only raises the average CMB temperature, an effect that is unobservable. At high redshift dissipation is mainly mediated through the quadrupole anisotropies, while after recombination peculiar motions are most important. During recombination the damping of the higher multipoles is also significant. We compute the average distortion for several examples using CosmoTherm, analyzing their dependence on parameters of the primordial power spectrum. For one of the best fit WMAP7 cosmologies, with n_S=1.027 and n_run=-0.034, the cooling of baryonic matter practically compensates the heating from acoustic dissipation in the mu-era. (abridged)

In this paper we refine the calculation of primordial helium recombination, accounting for several additional effects that were neglected or treated more approximately in previous studies. These include consideration of (i) time-dependent radiative transfer interaction between the 2^1 P_1 - 1^1 S_0 and 2^3 P_1 - 1^1 S_0 resonances; (ii) time-dependent radiative transfer for the partially overlapping n^1 P_1 - 1^1 S_0, n^1 D_2 - 1^1 S_0 and n^3 P_1 - 1^1 S_0 series with 3 <= n <= 10; (iii) electron scattering within a kernel approach. We also briefly discuss the effect of electron scattering and HI quadrupole lines on the recombination of hydrogen. Although the physics of all the considered processes is interesting and subtle, for the standard cosmology, with Y_p \sim 0.24, the overall correction to the ionization history during helium recombination with respect to the previous implementation of CosmoRec remains smaller than |DeltaNe/Ne| \sim 0.05%. The dominant improvement is caused by consistent inclusion of resonance scattering for the 2^1 P_1 - 1^1 S_0 resonance. For cosmologies with a large helium fraction, Y_p \sim 0.4, the difference reaches |D N_e/N_e|\sim 0.22% at z\sim1800, however, the overall correction to the CMB power spectra is small, exceeding |D C_l/C_l| \sim 0.05 % only at l>\sim 3000. In comparison to the current version of Recfast the difference reaches |D C_l/C_l|\sim0.4% at l\sim 3000 for Y_p \sim 0.4, and also for the standard value Y_p \sim 0.24 we find differences |DC_l/C_l|>\sim 0.1% at l\sim 2500. The new processes are now included by the cosmological recombination code CosmoRec and can be activated as needed for most settings without affecting its runtime significantly.

The difference in the adiabatic indices of photons and non-relativistic baryonic matter in the early Universe causes the electron temperature to be slightly lower than the radiation temperature. Thermalization of photons with a colder plasma results in the accumulation of photons in the Rayleigh-Jeans tail, aided by stimulated recoil, while the higher frequency spectrum tries to approach Planck spectrum at the electron temperature $T_{\gamma}^{final}=\Te<T_{\gamma}^{initial}$; i.e., Bose-Einstein condensation of photons occurs. We find new solutions of the Kompaneets equation describing this effect. No actual condensate is, in reality, possible since the process is very slow and photons drifting to low frequencies are efficiently absorbed by bremsstrahlung and double Compton processes. The spectral distortions created by Bose-Einstein condensation of photons are within an order of magnitude (for the present range of allowed cosmological parameters), with exactly the same spectrum but opposite in sign, of those created by diffusion damping of the acoustic waves on small scales corresponding to comoving wavenumbers $45< k< 10^4\, Mpc^{-1}$. The initial perturbations on these scales are completely unobservable today due to their being erased completely by Silk damping. There is partial cancellation of these two distortions, leading to suppression of $\mu$ distortions expected in the standard model of cosmology. The net distortion depends on the scalar power index $n_S$ and its running $d n_S/d\ln k$, and may vanish for special values of parameters, for example, for a running spectrum with, $n_S=1,d n_S/d\ln k=-0.038$. We arrive at an intriguing conclusion: even a null result, non-detection of $\mu$-type distortion at a sensitivity of $10^{-9}$, gives a quantitative measure of the primordial small-scale power spectrum.

In the early Universe, energy stored in small-scale density perturbations is quickly dissipated by Silk-damping, a process that inevitably generates mu- and y-type spectral distortions of the cosmic microwave background (CMB). These spectral distortions depend on the shape and amplitude of the primordial power spectrum at wavenumbers k < 10^4 Mpc^{-1}. Here we study constraints on the primordial power spectrum derived from COBE/FIRAS and forecasted for PIXIE. We show that measurements of mu and y impose strong bounds on the integrated small-scale power, and we demonstrate how to compute these constraints using k-space window functions that account for the effects of thermalization and dissipation physics. We show that COBE/FIRAS places a robust upper limit on the amplitude of the small-scale power spectrum. This limit is about three orders of magnitude stronger than the one derived from primordial black holes in the same scale range. Furthermore, this limit could be improved by another three orders of magnitude with PIXIE, potentially opening up a new window to early Universe physics. To illustrate the power of these constraints, we consider several generic models for the small-scale power spectrum predicted by different inflation scenarios, including running-mass inflation models and inflation scenarios with episodes of particle production. PIXIE could place very tight constraints on these scenarios, potentially even ruling out running-mass inflation models if no distortion is detected. We also show that inflation models with sub-Planckian field excursion that generate detectable tensor perturbations should simultaneously produce a large CMB spectral distortion, a link that could potentially be established by PIXIE.

The results for the total multi-photon decay rates of the 3p and 4s levels of hydrogen, presented by D. Solovyev and L. Labzowsky within the cascade approximation, are revisited. The corrected results for certain decay channels differ from original ones of those authors sometimes by order of magnitude. Some aspects with respect to the cosmological recombination process are clarified.

Cosmological parameter measurements from CMB experiments such as Planck, ACTpol, SPTpol and other high resolution follow-ons fundamentally rely on the accuracy of the assumed recombination model, or one with well prescribed uncertainties. Deviations from the standard recombination history might suggest new particle physics or modified atomic physics. Here we treat possible perturbative fluctuations in the free electron fraction, $\Xe(z)$, by a semi-blind expansion in densely-packed modes in redshift. From these we construct parameter eigenmodes, which we rank order so that the lowest modes provide the most power to probe the $\Xe(z)$ with CMB measurements. Since the eigenmodes are effectively weighed by the fiducial $\Xe$ history, they are localized around the differential visibility peak, allowing for an excellent probe of hydrogen recombination, but a weaker probe of the higher redshift helium recombination and the lower redshift highly neutral freeze-out tail. We use an information-based criterion to truncate the mode hierarchy, and show that with even a few modes the method goes a long way towards morphing a fiducial older {\sc Recfast} $X_{\rm e,i} (z)$ into the new and improved {\sc CosmoRec} and {\sc HyRec} $X_{\rm e,f} (z)$ in the hydrogen recombination regime, though not well in the helium regime. Without such a correction, the derived cosmic parameters are biased. We discuss an iterative approach for updating the eigenmodes to further hone in on $X_{\rm e,f} (z)$ if large deviations are indeed found. We also introduce control parameters that downweight the attention on the visibility peak structure, e.g., focusing the eigenmode probes more strongly on the $\Xe (z)$ freeze-out tail, as would be appropriate when looking for the $\Xe$ signature of annihilating or decaying elementary particles.

The energy spectrum of the cosmic microwave background (CMB) allows constraining episodes of energy release in the early Universe. In this paper we revisit and refine the computations of the cosmological thermalization problem. For this purpose a new code, called CosmoTherm, was developed that allows solving the coupled photon-electron Boltzmann equation in the expanding, isotropic Universe for small spectral distortion in the CMB. We explicitly compute the shape of the spectral distortions caused by energy release due to (i) annihilating dark matter; (ii) decaying relict particles; (iii) dissipation of acoustic waves; and (iv) quasi-instantaneous heating. We also demonstrate that (v) the continuous interaction of CMB photons with adiabatically cooling non-relativistic electrons and baryons causes a negative mu-type CMB spectral distortion of DI_nu/I_nu ~ 10^{-8} in the GHz spectral band. We solve the thermalization problem including improved approximations for the double Compton and Bremsstrahlung emissivities, as well as the latest treatment of the cosmological recombination process. At redshifts z <~ 10^3 the matter starts to cool significantly below the temperature of the CMB so that at very low frequencies free-free absorption alters the shape of primordial distortions significantly. In addition, the cooling electrons down-scatter CMB photons introducing a small late negative y-type distortion at high frequencies. We also discuss our results in the light of the recently proposed CMB experiment Pixie, for which CosmoTherm should allow detailed forecasting. Our current computations show that for energy injection because of (ii) and (iv) Pixie should allow to improve existing limits, while the CMB distortions caused by the other processes seem to remain unobservable with the currently proposed sensitivities and spectral bands of Pixie.

(Context) We calculate constraints from current and future cosmic microwave background (CMB) measurements on annihilating dark matter (DM) with masses below the electroweak scale: m_{DM}=5-100 GeV. In particular, we focus our attention on the lower end of this mass range, as DM particles with masses m_{DM} ~ 10 GeV have been recently claimed to be consistent with the CoGeNT and DAMA/LIBRA results, while also providing viable DM candidates to explain the measurements of Fermi and WMAP haze. (Aims) We study the model (in)dependence of CMB spectrum on particle physics DM models, large scale structure formation and cosmological uncertainties. We attempt to find a simple and practical recipe for estimating current and future CMB bounds on a broad class of DM annihilation models. (Results) We show that in the studied DM mass range the CMB signal of DM annihilations is independent of the details of large scale structure formation, distribution and profile of DM halos and other cosmological uncertainties. All particle physics models of DM annihilation can be described with only one parameter, the fraction of energy carried away by neutrinos in DM annihilation. As the main result we provide a simple and rather generic fitting formula for calculating CMB constraints on the annihilation cross section of light WIMPs. We show that thermal relic DM in the CoGeNT, DAMA favoured mass range is in a serious conflict with present CMB data for the annihilation channels with few neutrinos, and will definitely be tested by the Planck mission for all possible DM annihilation channels. Also, our findings strongly disfavor the claim that thermal relic DM annihilations with m_{DM} ~ 10 GeV and $<sigma_av> ~ 9x10^{-25} cm^{3}s^{-1} could be a cause for Fermi and WMAP haze.

It is well known that our motion with respect to the cosmic microwave background (CMB) rest frame introduces a large dipolar CMB anisotropy, with an amplitude ~beta=v/c~10^{-3}. In addition it should lead to a small breaking of statistical isotropy which becomes most notable at higher multipoles. In principle this could be used to determine our velocity with respect to the CMB rest frame using high angular resolution data from Planck, without directly relying on the amplitude and direction of the CMB dipole, allowing us to constrain cosmological models in which the cosmic dipole arises partly from large-scale isocurvature perturbations instead of being fully motion-induced. Here we derive simple recursion relations that allow precise computation of the motion-induced coupling between different spherical harmonic coefficients. Although the lowest order approximations for the coupling kernel can be deficient by factors of 2-5 at multipoles l~1000-3000, using our results for the aberration kernel we explicitly confirm that for a statistical detection of the aberration effect only first order terms in beta matter. However, the expressions given here are not restricted to beta~10^{-3}, but can be used at much higher velocities. We demonstrate the robustness of these formulae, illustrating the dependence of the kernel on beta, as well as the spherical harmonic indices l and m.

We use the new cosmological recombination code, CosmoRec, for parameter estimation in the context of (future) precise measurements of the CMB temperature and polarization anisotropies. We address the question of how previously neglected physical processes in the recombination model of Recfast affect the determination of key cosmological parameters, for the first time performing a model-by-model computation of the recombination problem. In particular we ask how the biases depend on different combinations of parameters, e.g. when varying the helium abundance or the effective number of neutrino species in addition to the standard six parameters. We also forecast how important the recombination corrections are for a combined Planck, ACTPol and SPTpol data analysis. Furthermore, we ask which recombination corrections are really crucial for CMB parameter estimation, and whether an approach based on a redshift-dependent correction function to Recfast is sufficient in this context.

A new approach to the cosmological recombination problem is presented, which completes our previous analysis on the effects of two-photon processes during the epoch of cosmological hydrogen recombination, accounting for ns-1s and nd-1s Raman events and two-photon transitions from levels with n>=2. The recombination problem for hydrogen is described using an effective 400-shell multi-level approach, to which we subsequently add all important recombination corrections discussed in the literature thus far. We explicitly solve the radiative transfer equation of the Lyman-series photon field to obtain the required modifications to the rate equations of the resolved levels. In agreement with earlier computations we find that 2s-1s Raman scattering leads to a delay in recombination by DN_e/N_e~0.9% at z~920. Two-photon decay and Raman scattering from higher levels (n>3) result in a small additional modifications, and precise results can be obtained when including their effect for the first 3-5 shells. This work is a major step towards a new cosmological recombination code (CosmoRec) that supersedes the physical model included in Recfast, and which, owing to its short runtime, can be used in the analysis of future CMB data from the Planck Surveyor.

In this paper we discuss the effect of recombinations to highly excited states (n > 100) in hydrogen during the cosmological recombination epoch. For this purpose, we developed a new ODE solver for the recombination problem, based on an implicit Gear's method. This solver allows us to include up to 350 l-resolved shells or ~61 000 separate levels in the hydrogen model and to solve the recombination problem for one cosmology in ~27 hours. This is a huge improvement in performance over our previous recombination code, for which a 100-shell computation (5050 separate states) already required ~150 hours on a single processor. We show that for 350 shells down to redshift z ~200 the results for the free electron fraction have practically converged. The final modification in the free electron fraction at z ~200 decreases from about \DeltaNe/Ne ~2.8% for 100 shells to \DeltaNe/Ne ~1.6% for 350 shells. However, the associated changes in the CMB power spectra at large multipoles l are rather small, so that for accurate computations in connection with the analysis of Planck data already ~100 shells are expected to be sufficient. Nevertheless, the total value of \tau could still be affected at a significant level. We also briefly investigate the effect of collisions on the recombination dynamics. With our current estimates for the collisional rates we find a correction of \DeltaNe/Ne ~ -0.088% at z ~ 700, which is mainly caused by l-changing collisions with protons. Furthermore, we present results on the cosmological recombination spectrum, showing that at low frequencies collisional processes are important. However, the current accuracy of collisional rates is insufficient for precise computations of templates for the recombination spectrum at \nu<~1 GHz, and also the effect of collisions on the recombination dynamics suffers from the uncertainty in these rates.

[Abridged] We use our most recent training set for the RICO code to estimate the impact of recombination uncertainties on the posterior probability distributions which will be obtained from future CMB experiments, and in particular the Planck satellite. Using a MCMC analysis to sample the posterior distribution of the cosmological parameters, we find that Planck will have biases of -0.7, -0.3 and -0.4 sigmas for n_S, Omega_b h2 and log(As), respectively, in the minimal 6-parameter LCDM model, if the description of the recombination history given by RICO is not used. The remaining parameters are not significantly affected. We also show, that the cosmology dependence of the corrections to the recombination history modeled with RICO has a negligible impact on the posterior distributions obtained for the case of the Planck satellite. In practice, this implies that the inclusion of additional corrections to existing recombination codes can be achieved using simple cosmology-independent `fudge functions'. Finally, we also investigated the impact of some recent improvements in the treatment of hydrogen recombination which are still not included in the current version of our training set for Rico, by assuming that the cosmology dependence of those corrections can be neglected. In summary, with our current understanding of the complete recombination process, the expected biases in the cosmological parameters inferred from Planck might be as large as -2.3, -1.7 and -1 sigmas for n_S, Omega_b h2 and log(As) respectively, if all those corrections are not taken into account. We note that although the list of physical processes that could be of importance for Planck seems to be nearly complete, still some effort has to be put in the validation of the results obtained by the different groups.

In this paper we explore the potential effects of DM annihilations on the cosmological recombination spectrum. With this example we want to demonstrate that the cosmological recombination spectrum in principle is sensitive to details related to possible extra energy release during recombination. We restrict ourselves to DM models which produce a negligible primordial distortion of the CMB energy spectrum. However, since during the epoch of cosmological recombination a large fraction of the deposited energy can directly go into ionizations and excitations of neutral atoms, both the cosmological recombination spectrum and ionization history can still be affected significantly. We compute the modifications to the cosmological recombination spectrum using our multi-level HI and HeI recombination code, showing that additional photons are created due to uncompensated loops of transitions which are induced by DM annihilations. As we illustrate here, the results depend on the detailed branching of the deposited energy into heating, ionizations and excitations. This dependence in principle should allow us to shed light on the nature of the underlying annihilating DM model (or more generally speaking, the mechanism leading to energy injection) when measuring the cosmological recombination spectrum. However, for current upper limits on the potential DM annihilation rate during recombination the cosmological recombination spectrum is only affected at the level of a few percent. Nevertheless, we argue here that the cosmological recombination spectrum would provide another independent and very direct way of checking for the presence of sources of extra ionizing or exciting photons at high redshifts. This would open an new window to possible (non-standard) processes occurring (abridged)

In this paper we consider the re-processing of high frequency photons emitted by HeII and HeI during the epoch of cosmological recombination by HeI and HI. We demonstrate that, in comparison to computations which neglect all feedback processes, the number of cosmological recombination photons that are related to the presence of helium in the early Universe could be increased by ~40%-70%. Our computations imply that per helium nucleus ~3-6 additional photons could be produced. Therefore, a total of ~12-14 helium-related photons are emitted during cosmological recombination. This is an important addition to cosmological recombination spectrum which in the future may render it slightly easier to determine the primordial abundance of helium using differential measurements of the CMB energy spectrum. Also, since these photons are the only witnesses of the feedback process at high redshift, observing them in principle offers a way to check our understanding of the recombination physics. Here most interestingly, the feedback of HeII photons on HeI leads to the appearance of several additional, rather narrow spectral features in the HeI recombination spectrum at low frequencies. Consequently, the signatures of helium-related features in the CMB spectral distortion due to cosmological recombination at some given frequency can exceed the average level of ~17% several times. We find that in particular the bands around nu ~10GHz, ~35GHz, ~80GHz, and ~200GHz seem to be affected strongly. In addition, we computed the changes in the cosmological ionization history, finding that only the feedback of primary HeI photons on the dynamics of HeII-->HeI recombination has an effect, producing a change of DN_e/N_e ~+ 0.17% at z~2300. This result seems to be ~2-3 times smaller than the one obtained in earlier computations for this process (abridged).

The physical ingredients to describe the epoch of cosmological recombination are amazingly simple and well-understood. This fact allows us to take into account a very large variety of physical processes, still finding potentially measurable consequences for the energy spectrum and temperature anisotropies of the Cosmic Microwave Background (CMB). In this contribution we provide a short historical overview in connection with the cosmological recombination epoch and its connection to the CMB. Also we highlight some of the detailed physics that were studied over the past few years in the context of the cosmological recombination of hydrogen and helium. (abrigded)

In this paper we consider the effects of resonance and electron scattering on the escape of Lyman alpha photons during cosmological hydrogen recombination. We pay particular attention to the influence of atomic recoil, Doppler boosting and Doppler broadening using a Fokker-Planck approximation of the redistribution function describing the scattering of photons on the Lyman alpha resonance of moving hydrogen atoms. We extend the computations of our recent paper on the influence of the 3d/3s-1s two-photon channels on the dynamics of hydrogen recombination, simultaneously including the full time-dependence of the problem, the thermodynamic corrections factor, leading to a frequency-dependent asymmetry between the emission and absorption profile, and the quantum-mechanical corrections related to the two-photon nature of the 3d/3s-1s emission and absorption process on the exact shape of the Lyman alpha emission profile. We show here that due to the redistribution of photons over frequency hydrogen recombination is sped up by DN_e/N_e~-0.6% at z~900. For the CMB temperature and polarization power spectra this results in |DC_l/C_l|~0.5%-1% at l >~ 1500, and therefore will be important for the analysis of future CMB data in the context of the PLANCK Surveyor, SPT and ACT. The main contribution to this correction is coming from the atomic recoil effect (DN_e/N_e~-1.2% at z~900), while Doppler boosting and Doppler broadening partially cancel this correction, again slowing hydrogen recombination down by DN_e/N_e~0.6% at z~900. The influence of electron scattering close to the maximum of the Thomson visibility function at z~1100 can be neglected. (abridged)

We give a formulation of the radiative transfer equation for Lyman alpha photons which allows us to include the two-photon corrections for the 3s-1s and 3d-1s decay channels during cosmological hydrogen recombination. We use this equation to compute the corrections to the Sobolev escape probability for Lyman alpha photons during hydrogen recombination, which then allow us to calculate the changes in the free electron fraction and CMB temperature and polarization power spectra. We show that the effective escape probability changes by DP/P ~+ 11% at z~1400 in comparison with the one obtained using the Sobolev approximation. This speeds up of hydrogen recombination by DN_e/N_e ~- 1.6% at z~1190, implying |DC_l/C_l| ~1%-3% at l >~ 1500 with shifts in the positions of the maxima and minima in the CMB power spectra. These corrections will be important for the analysis of future CMB data. The total correction is the result of the superposition of three independent processes, related to (i) time-dependent aspects of the problem, (ii) corrections due to quantum mechanical deviations in the shape of the emission and absorption profiles in the vicinity of the Lyman alpha line from the normal Lorentzian, and (iii) a thermodynamic correction factor, which occurs to be very important. All these corrections are neglected in the Sobolev-approximation, but they are important in the context of future CMB observations. All three can be naturally obtained in the two-photon formulation of the Lyman alpha absorption process. However, the corrections (i) and (iii) can also be deduced in the normal '1+1' photon language, without necessarily going to the two-photon picture. Therefore only (ii) is really related to the quantum mechanical aspects of the two-photon process (abridged)

We consider the effects connected with the detailed radiative transfer during the epoch of cosmological recombination on the ionization history of our Universe. We focus on the escape of photons from the hydrogen Lyman-alpha resonance at redshifts 600<~ z <~ 2000, one of two key mechanisms defining the rate of cosmological recombination. We approach this problem within the standard formulation, and corrections due to two-photon interactions are deferred to another paper. As a main result we show here that within a non-stationary approach to the escape problem, the resulting correction in the free electron fraction, N_e, is about ~1.6-1.8% in the redshift range 800<~z<~1200. Therefore the discussed process results in one of the largest modifications to the ionization history close to the maximum of Thomson-visibility function at z~1100 considered so far. We prove our results both numerically and analytically, deriving the escape probability, and considering both Lyman-alpha line emission and line absorption in a way different from the Sobolev approximation. In particular, we give a detailed derivation of the Sobolev escape probability during hydrogen recombination, and explain the underlying assumptions. We then discuss the escape of photons for the case of coherent scattering in the lab frame, solving this problem analytically in the quasi-stationary approximation and also in the time-dependent case. We show here that during hydrogen recombination the Sobolev approximation for the escape probability is not valid at the level of DP/P~5-10%. This is because during recombination the ionization degree changes significantly over a characteristic time Dz/z~10%, so that at percent level accuracy the photon distribution is not evolving (abridged)

We present Rico, a code designed to compute the ionization fraction of the Universe during the epoch of hydrogen and helium recombination with an unprecedented combination of speed and accuracy. This is accomplished by training the machine learning code Pico on the calculations of a multi-level cosmological recombination code which self-consistently includes several physical processes that were neglected previously. After training, Rico is used to fit the free electron fraction as a function of the cosmological parameters. While, for example at low redshifts (z<~900), much of the net change in the ionization fraction can be captured by lowering the hydrogen fudge factor in Recfast by about 3%, Rico provides a means of effectively using the accurate ionization history of the full recombination code in the standard cosmological parameter estimation framework without the need to add new or refined fudge factors or functions to a simple recombination model. Within the new approach presented here it is easy to update Rico whenever a more accurate full recombination code becomes available. Once trained, Rico computes the cosmological ionization history with negligible fitting error in ~10 milliseconds, a speed-up of at least 10^6 over the full recombination code that was used here. Also Rico is able to reproduce the ionization history of the full code to a level well below 0.1%, thereby ensuring that the theoretical power spectra of CMB fluctuations can be computed to sufficient accuracy and speed for analysis from upcoming CMB experiments like Planck. Furthermore it will enable cross-checking different recombination codes across cosmological parameter space, a comparison that will be very important in order to assure the accurate interpretation of future cosmic microwave background data.

We discuss a new solution of the Kompaneets-equation for physical situations in which low frequency photons, forming relatively narrow spectral details, are Compton scattered in an isotropic, infinite medium with an intense ambient blackbody field that is very close to full thermodynamic equilibrium with the free electrons. In this situation the background-induced stimulated Compton scattering slows down the motion of photons toward higher frequencies by a factor of 3 in comparison with the solution that only takes into account Doppler-broadening and boosting. This new solution is important for detailed computations of cosmic microwave background spectral distortions arising due to uncompensated atomic transitions of hydrogen and helium in the early Universe. In addition we derive another analytic solution that only includes the background-induced stimulated Compton scattering and is valid for power-law ambient radiation fields. This solution might have interesting applications for radio lines arising inside of bright extra-galactic radio source, where according to our estimates line shifts because of background-induced stimulated scattering could be amplified and even exceed the line broadening due to the Doppler-effect.

Energy release in the early Universe (z<~ 2x10^6) should lead to some broad spectral distortion of the cosmic microwave background (CMB) radiation field, which can be characterized as y-type distortion when the injection process started at redshifts z<~ 5x10^4. Here we demonstrate that if energy was released before the beginning of cosmological hydrogen recombination (z~1400), closed loops of bound-bound and free-bound transitions in HI and HeII lead to the appearance of (i) characteristic multiple narrow spectral features at dm and cm wavelengths, and (ii) a prominent sub-millimeter feature consisting of absorption and emission parts in the far Wien tail of CMB spectrum. The additional spectral features are generated in the pre-recombinational epoch of HI (z>~1800) and HeII (z>~7000), and therefore differ from those arising due to normal cosmological recombination in the undisturbed CMB blackbody radiation field. We present the results of numerical computations including 25 atomic shells for both HI and HeII, and discuss the contributions of several individual transitions in detail. As examples, we consider the case of instantaneous energy release (e.g. due to phase transitions) and exponential energy release because of long-lived decaying particles. Our computations show that due to possible pre-recombinational atomic transitions the variability of the CMB spectral distortion increases when comparing with the distortions arising in the normal recombination epoch. The existence of these narrow spectral features would open an unique way to separate y-distortions due to pre-recombinational ($1400<~ z <~5x10^4) energy release from those arising in the post-recombinational era at redshifts z<~800. (abridged)

The main goal of this work is to calculate the contributions to the cosmological recombination spectrum due to bound-bound transitions of helium. We show that due to the presence of helium in the early Universe unique features appear in the total cosmological recombination spectrum. These may provide a unique observational possibility to determine the relative abundance of primordial helium, well before the formation of first stars. We include the effect of the tiny fraction of neutral hydrogen atoms on the dynamics of HeII -> HeI recombination at redshifts $z\sim 2500$. As discussed recently, this process significantly accelerates HeII -> HeI recombination, resulting in rather narrow and distinct features in the associated recombination spectrum. In addition this process induces some emission within the hydrogen Lyman-$\alpha$ line, before the actual epoch of hydrogen recombination round $z\sim 1100-1500$. We also show that some of the fine structure transitions of neutral helium appear in absorption, again leaving unique traces in the Cosmic Microwave Background blackbody spectrum, which may allow to confirm our understanding of the early Universe and detailed atomic physics.

We study the two-photon process for the transitions ns --> 1s and nd --> 1s in hydrogen up to large n. For n<=20 we provide simple analytic fitting formulae to describe the non-resonant part of the two-photon emission profiles. Combining these with the analytic form of the cascade-term yields a simple and accurate description of the full two-photon decay spectrum, which only involves a sum over a few intermediate states. We demonstrate that the cascade term naturally leads to a nearly Lorentzian shape of the two-photon profiles in the vicinity of the resonances. However, due to quantum-electrodynamical corrections, the two-photon emission spectra deviate significantly from the Lorentzian shape in the very distant wings of the resonances. We investigate up to which distance the two-photon profiles are close to a Lorentzian and discuss the role of the interference term. We then analyze how the deviation of the two-photon profiles from the Lorentzian shape affects the dynamics of cosmological hydrogen recombination. Since in this context the escape of photons from the Lyman-alpha resonance plays a crucial role, we concentrate on the two-photon corrections in the vicinity of the Lyman-alpha line. Our computations show that the changes in the ionization history due to the additional two-photon process from high shell (n>2) likely do not reach the percent-level. For conservative assumptions we find a correction DN_e/N_e~-0.4% at redshift z~1160. This is numerically similar to the result of another recent study; however, the physics leading to this conclusion is rather different. In particular, our calculations of the effective two-photon decay rates yield significantly different values, where the destructive interference of the resonant and non-resonant terms plays a crucial role in this context (abridged)

The physical ingredients to describe the epoch of cosmological recombination are amazingly simple and well-understood. This fact allows us to take into account a very large variety of processes, still finding potentially measurable consequences. In this contribution we highlight some of the detailed physics that were recently studied in connection with cosmological hydrogen and helium recombination. The impact of these considerations is two-fold: (i) the associated release of photons during this epoch leads to interesting and unique deviations of the Cosmic Microwave Background (CMB) energy spectrum from a perfect blackbody, which, in particular at decimeter wavelength, may become observable in the near future. Despite the fact that the abundance of helium is rather small, it also contributes a sizeable amount of photons to the full recombination spectrum, which, because of differences in the dynamics of the helium recombinations and the non-trivial superposition of all components, lead to additional distinct spectral features. Observing the spectral distortions from the epochs of hydrogen and helium recombination, in principle would provide an additional way to determine some of the key parameters of the Universe (e.g. the specific entropy, the CMB monopole temperature and the pre-stellar abundance of helium), not suffering from limitations set by cosmic variance. Also it permits us to confront our detailed understanding of the recombination process with direct observational evidence. (ii) with the advent of high precision CMB data, e.g. as will be available using the Planck Surveyor or CMBpol, a very accurate theoretical understanding of the ionization history of the Universe becomes necessary for the interpretation of the CMB temperature and polarization anisotropies. (abridged)

The physical ingredients to describe the epoch of cosmological recombination are amazingly simple and well-understood. This fact allows us to take into account a very large variety of processes, still finding potentially measurable consequences. In this contribution we highlight some of the detailed physics that were recently studied in connection with cosmological hydrogen recombination. The impact of these considerations is two-fold: (i) the associated release of photons during this epoch leads to interesting and unique deviations of the Cosmic Microwave Background (CMB) energy spectrum from a perfect blackbody, which, in particular at decimeter wavelength, may become observable in the near future. Observing these distortions, in principle would provide an additional way to determine some of the key parameters of the Universe (e.g. the specific entropy, the CMB monopole temperature and the pre-stellar abundance of helium), not suffering from limitations set by cosmic variance. Also it permits us to confront our detailed understanding of the recombination process with direct observational evidence. In this contribution we illustrate how the theoretical spectral template for the cosmological recombination spectrum may be utilized for this purpose. (ii) with the advent of high precision CMB data, e.g. as will be available using the Planck Surveyor or Cmbpol, a very accurate theoretical understanding of the ionization history of the Universe becomes necessary for the interpretation of the CMB temperature and polarization anisotropies. (abridged)

The recombination history of the Universe depends exponentially on the temperature, T_0, of the cosmic microwave background. Therefore tiny changes of T_0 are expected to lead to significant changes in the free electron fraction. Here we show that even the current 1sigma-uncertainty in the value of T_0 results in more than half a percent ambiguity in the ionization history, and more than 0.1% uncertainty in the TT and EE power spectra at small angular scales. We discuss how the value of T_0 affects the highly redshifted cosmological hydrogen recombination spectrum and demonstrate that T_0 could, in principle, be measured by looking at the low frequency distortions of the cosmic microwave background spectrum. For this no absolute measurements are necessary, but sensitivities on the level of ~30nK are required to extract the quasi-periodic frequency-dependent signal with typical Delta nu/nu~0.1 coming from cosmological recombination. We also briefly mention the possibility of obtaining additional information on the specific entropy of the Universe, and other cosmological parameters.

We compute the corrections to the cosmological hydrogen recombination history due to delayed feedback of Lyman-series photons and the escape in the Lyman-continuum. The former process is expected to slightly delay recombination, while the latter should allow the medium to recombine a bit faster. It is shown that the subsequent feedback of released Lyman-n photons on the lower lying Lyman-(n-1) transitions yields a maximal correction of DN_e/N_e 0.22% at z~ 1050. Including only Lyman-\beta feedback onto the Lyman-\alpha transition, accounts for most of the effect. We find corrections to the cosmic microwave background TT and EE power spectra \change{with typical peak to peak amplitude |DC^{TT}_l/C^{TT}_l|~0.15% and |\Delta C^{EE}_l/C^{EE}_l|~0.36% at l<~3000. The escape in the Lyman-continuum and feedback of Lyman-\alpha photons on the photoionization rate of the second shell lead to modifications of the ionization history which are very small (less than |DN_e/N_e|~few x 10^{-6}).

We provide simple and accurate analytic approximations for the low frequency double Compton emission coefficient that are applicable in a broad range of physical situations up to mildly relativistic temperatures. These approximations may be useful for checking in which circumstances the double Compton process is important. We perform series expansions of the double Compton emission integral for low energies of the incident photon and electron and compare the derived analytic expressions with the results obtained by numerical integrations of the full double Compton cross section. We explicitly derived analytic approximations to the low frequency double Compton emission coefficient for initial monochromatic photons and Wien spectra. We show that combining the analytic approximations given in this paper, an accuracy of better than 5% over a very broad range of temperatures and under various physical conditions can be achieved. The double Compton emissivity strongly depends on the ratio of the incoming photon and electron energies: for hard photons and cold electrons the emission is strongly suppressed compared to the case of similar photon and electron energy, whereas in the opposite situation, i.e. hot electrons and soft initial photons, the emission is enhanced. For photons and electrons close to thermodynamic equilibrium the double Compton emissivity increases less rapidly with temperature than in the Lightman-Thorne approximation and the corrections exceed ~10% for temperatures above 4 keV.

We present results for the spectral distortions of the Cosmic Microwave Background (CMB) arising due to bound-bound transitions during the epoch of cosmological hydrogen recombination at frequencies down to nu~100MHz. We extend our previous treatment of the recombination problem now including the main collisional processes and following the evolution of all the hydrogen angular momentum sub-states for up to 100 shells. We show that, due to the low baryon density of the Universe, even within the highest considered shell full statistical equilibrium (SE) is not reached and that at low frequencies the recombination spectrum is significantly different when assuming full SE for n>2. We also directly compare our results for the ionization history to the output of the Recfast code, showing that especially at low redshifts rather big differences arise. In the vicinity of the Thomson visibility function the electron fraction differs by roughly -0.6% which affects the temperature and polarization power spectra by <~1%. Furthermore we shortly discuss the influence of free-free absorption and line broadening due to electron scattering on the bound-bound recombination spectrum and the generation of CMB angular fluctuations due to scattering of photons within the high shells.

In this paper we compute the emission coming from the direct recombination of free electrons to a given shell (n>=2) during the epoch of cosmological hydrogen recombination. This contribution leads to a total of one photon per recombined hydrogen atom and therefore a ~30-88% increase of the recombination spectrum within the frequency range 1 GHz<=nu<=100 GHz. In particular the Balmer-continuum emission increases the distortion at nu~ 690 GHz by ~92%. With our 100 shell calculations for the hydrogen atom we find that a total of ~5 photons per hydrogen atom are emitted when including all the bound-bound transitions, the 2s two-photon decay channel and the optically thin free-bound transitions. Since the direct recombination continuum at high n is very broad only a few n-series continuua are distinguishable and most of this additional emission below nu<~30 GHz is completely featureless.

We compute the spectral distortions of the Cosmic Microwave Background (CMB) arising during the epoch of cosmological hydrogen recombination within the standard cosmological (concordance) model for frequencies in the range 1 GHz-3500 GHz. We follow the evolution of the populations of the hydrogen levels including states up to principle quantum number $n=30$ in the redshift range $500\leq z\leq 3500$. All angular momentum sub-states are treated individually, resulting in a total number of 465 hydrogen levels. The evolution of the matter temperature and the fraction of electrons coming from HeII are also included. We present a detailed discussion of the distortions arising from the main dipolar transitions, e.g. Lyman and Balmer series, as well as the emission due to the two-photon decay of the hydrogen 2s level. Furthermore, we investigate the robusteness of the results against changes in the number of shells considered. The resulting spectral distortions have a characteristic oscillatory behaviour, which might allow experimentalists to separate them from other backgrounds. The relative distortion of the spectrum exceeds a value of $10^{-7}$ at wavelengths longer than 21cm. Our results also show the importance of detailed follow-up of the angular momentum sub-states, and their effect on the amplitude of the lines. The effect on the residual electron fraction is only moderate, and mainly occurs at low redshifts. The CMB angular power spectrum is changed by less than 1%. Finally, our computations show that if the primordial radiation field is described by a pure blackbody, then there is no significant emission from any hydrogen transition at redshifts greater than $z \sim 2000$. This is in contrast to some earlier works, where the existence of a `pre-recombination' peak was claimed.

Ultraluminous X-ray sources fed by Wolf-Rayet star winds and X-ray bursters in ultracompact binaries with He or C white dwarfs have accretion disks, whose properties may significantly differ from those of pure H alpha-accretion disks. We have therefore included the dependence on charge number Z and mean molecular weights mu_{e/I} into the Shakura and Sunyaev (1973) scaling relations for the key parameters of the disk. Furthermore, we also consider the case of the pseudo-Newtonian potential of Paczynsky and Wiita (1980). These scaling relations might become useful, e.g., when estimating the illumination efficiency of the external parts of the disk. We also address the changes in the structure of the boundary (spreading) layer on the surface of neutron stars, occurring in the case of H depleted accretion disks.

Induced emission due to the presence of soft CMB photons slightly increases the two-photon decay rate of the 2s level of hydrogen defining the rate of cosmological recombination. This correspondingly changes the degree of ionization, the visibility function and the resulting primordial temperature anisotropies and polarization of the CMB on the percent level. These changes exceed the precision of the widely used CMBFAST and CAMB codes by more than one order of magnitude and can be easily taken into account.

In this work we consider the changes of the SZ cluster brightness, flux and number counts induced by the motion of the Solar System with respect to the frame defined by the cosmic microwave background (CMB). These changes are connected with the Doppler effect and aberration and exhibit a strong spectral and spatial dependence. The correction to the SZ cluster brightness and flux has an amplitude and spectral dependence, which is similar to the first order cluster peculiar velocity correction to the thermal SZ effect. Due to the change in the received cluster CMB flux the motion of the Solar System induces a dipolar asymmetry in the observed number of clusters above a given flux level. Similar effects were discussed for $\gamma$-ray bursts and radio galaxies, but here, due to the very peculiar frequency-dependence of the thermal SZ effect, the number of observed clusters in one direction of the sky can be both, decreased or increased depending on the frequency band. A detection of this asymmetry should be possible using future full sky CMB experiments with mJy sensitivities.

The CMB angular temperature fluctuations observed by COBE and WMAP enable us to place a lower limit on the spectral distortions of the CMB at any angular scale. These distortions are connected with the simple fact, that the superposition of blackbodies with different temperatures in general is not a blackbody. We show, that in the limit of small temperature fluctuations, the superposition of blackbodies leads to a y-type spectral distortion. It is known, that the CMB dipole induces a y-type spectral distortion with quadrupole and monopole angular distribution leading to a corresponding whole sky y-parameter of y~2.6\cdot 10^{-7}.We show here that taking the difference of the CMB signal in the direction of the maximum and minimum of the CMB dipole due to the superposition of two blackbodies leads to a spectral distortion with y~3.1\cdot 10^{-6}. The amplitude of this distortion can be calculated to the same precision as the CMB dipole, i.e. 0.3% today. Therefore it may be used as a source with brightness of several or tens of \muK to cross calibrate and calibrate different frequency channels of CMB surveys with precision of a few tens or hundreds of nK. Furthermore, we show in this work that primordial anisotropies for multipoles 2\leq l \leq 1000 also lead to spectral distortions but with much smaller y-parameter, i.e. y\sim 10^{-11}-10^{-9}.

We show how the temperature and the polarisation of the cosmic microwave background are affected by bulk rotation of clusters of galaxies owing to the kinetic Sunyaev-Zeldovich effect. The main effects of rotation are (i) a shift of the position of the peak of the temperature fluctuation relative to the center of the cluster by a few percent of the core radius and (ii) a tilt of the direction of the plane of linear polarisation by several degrees.