Measurements of the cosmic microwave background (CMB) spectral distortions (SDs) will open a new window on the very early universe, providing new information complementary to that gathered from CMB temperature and polarization anisotropies. In this paper, we study their synergy as a function of the characteristics of the considered experiments. In particular, we examine a wide range of sensitivities for possible SD measurements, spanning from FIRAS up to noise levels 1000 times better than PIXIE, and study their constraining power when combined with current or future CMB anisotropy experiments such as Planck or LiteBIRD plus CMB-S4. We consider a number of different cosmological models such as the $\Lambda$CDM, as well as its extensions with the running of the scalar spectral index, the decay or the annihilation of dark matter (DM) particles. While upcoming CMB anisotropy experiments will be able to decrease the uncertainties on inflationary parameters such as $A_s$ and $n_s$ by about a factor 2 in the $\Lambda$CDM case, we find that an SD experiment 10 times more sensitive than PIXIE (comparable to the proposed PRISM satellite) could potentially further contribute to constrain these parameters. This is even more significant in the case of the running of the scalar spectral index. Furthermore, as expected, constraints on DM particles decaying at redshifts probed by SDs will improve by orders of magnitude even with an experiment 10 times worse than PIXIE as compared to CMB anisotropies or Big Bang Nucleosynthesis bounds. On the contrary, DM annihilation constraints will not significantly improve over CMB anisotropy measurements. For some of the cases considered, we study the impact of marginalizing over the contribution from reionization and structure formation. Finally, we forecast the constraints obtainable with sensitivities achievable either from the ground or from a balloon.

The Ly-$\alpha$ forest 1D flux power spectrum is a powerful probe of several cosmological parameters. Assuming a $\Lambda$CDM cosmology including massive neutrinos, we find that the latest SDSS DR14 BOSS and eBOSS Ly-$\alpha$ forest data is in very good agreement with current weak lensing constraints on $(\Omega_m, \sigma_8)$ and has the same small level of tension with Planck. We did not identify a systematic effect in the data analysis that could explain this small tension, but we show that it can be reduced in extended cosmological models where the spectral index is not the same on the very different times and scales probed by CMB and Ly-$\alpha$ data. A particular case is that of a $\Lambda$CDM model including a running of the spectral index on top of massive neutrinos. With combined Ly-$\alpha$ and Planck data, we find a slight (3$\sigma$) preference for negative running, $\alpha_s= -0.010 \pm 0.004$ (68% CL). Neutrino mass bounds are found to be robust against different assumptions. In the $\Lambda$CDM model with running, we find $\sum m_\nu <0.11$ eV at the 95% confidence level for combined Ly-$\alpha$ and Planck (temperature and polarisation) data, or $\sum m_\nu < 0.09$ eV when adding CMB lensing and BAO data. We further provide strong and nearly model-independent bounds on the mass of thermal warm dark matter. For a conservative configuration consisting of SDSS data restricted to $z<4.5$ combined with XQ-100 \lya data, we find $m_X > 5.3\;\mathrm{keV}$ (95\%CL).

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.

Many attempts to solve the Hubble tension with extended cosmological models combine an enhanced relic radiation density, acting at the level of background cosmology, with new physical ingredients affecting the evolution of cosmological perturbations. Several authors have pointed out the ability of combined Baryon Acoustic Oscillation (BAO) and Big Bang Nucleosynthesis (BBN) data to probe the background cosmological history independently of both CMB maps and supernovae data. Using state-of-the-art assumptions on BBN, we confirm that combined BAO, deuterium, and helium data are in tension with the SH0ES measurements under the $\Lambda$CDM assumption at the 3.2$\sigma$ level, while being in close agreement with the CMB value. We subsequently show that floating the radiation density parameter $N_\mathrm{eff}$ only reduces the tension down to the 2.6$\sigma$ level. This conclusion, totally independent of any CMB data, shows that a high $N_\mathrm{eff}$ accounting for extra relics (either free-streaming or self-interacting) does not provide an obvious solution to the crisis, not even at the level of background cosmology. To circumvent this strong bound, (i) the extra radiation has to be generated after BBN to avoid helium bounds, and (ii) additional ingredients have to be invoked at the level of perturbations to reconcile this extra radiation with CMB and LSS data.

Einstein-Boltzmann Solvers (EBSs) are run on a massive scale by the cosmology community when fitting cosmological models to data. We present a new concept for speeding up such codes with neural networks. The originality of our approach stems from not substituting the whole EBS by a machine learning algorithm, but only its most problematic and least parallelizable step: the integration of perturbation equations over time. This approach offers two significant advantages: the task depends only on a subset of cosmological parameters, and it is blind to the characteristics of the experiment for which the output must be computed (for instance, redshift bins). These allow us to construct a fast and highly re-usable network. In this proof-of-concept paper, we focus on the prediction of CMB source functions, and design our networks according to physical considerations and analytical approximations. This allows us to reduce the depth and training time of the networks compared to a brute-force approach. Indeed, the calculation of the source functions using the networks is fast enough so that it is not a bottleneck in the EBS anymore. Finally, we show that their accuracy is more than sufficient for accurate MCMC parameter inference from Planck data. This paves the way for a new project, CosmicNet, aimed at gradually extending the use and the range of validity of neural networks within EBSs, and saving massive computational time in the context of cosmological parameter extraction.

We present a new efficient method to compute the angular power spectra of large-scale structure observables that circumvents the numerical integration over Bessel functions, expanding on a recently proposed algorithm based on FFTlog. This new approach has better convergence properties. The method is explicitly implemented in the CLASS code for the case of number count $C_\ell$'s (including redshift-space distortions, weak lensing, and all other relativistic corrections) and cosmic shear $C_\ell$'s. In both cases our approach speeds up the calculation of the exact $C_\ell$'s (without the Limber approximation) by a factor of order 400 at a fixed precision target of 0.1%.