The Epoch of Reionization (EoR) depends on the complex astrophysics governing the birth and evolution of the first galaxies and structures in the intergalactic medium. EoR models rely on cosmic microwave background (CMB) observations, and in particular the large-scale E-mode polarization power spectra (EE PS), to help constrain their highly uncertain parameters. However, rather than directly forward-modelling the EE PS, most EoR models are constrained using a summary statistic -- the Thompson scattering optical depth, $\tau_e$. Compressing CMB observations to $\tau_e$ requires adopting a basis set for the EoR history. The common choice is the unphysical, redshift-symmetric hyperbolic tangent (Tanh) function, which differs in shape from physical EoR models based on hierarchical structure formation. Combining public EoR and CMB codes, 21cmFAST and CLASS, here we quantify how inference using the $\tau_e$ summary statistic impacts the resulting constraints on galaxy properties and EoR histories. Using the last Planck 2018 data release, we show that the marginalized constraints on the EoR history are more sensitive to the choice of the basis set (Tanh vs physical model) than to the CMB likelihood statistic ($\tau_e$ vs PS). For example, EoR histories implied by the growth of structure show a small tail of partial reionization extending to higher redshifts: the lower limit on the volume-averaged neutral hydrogen fraction at $z=10$ changes from $\bar{x}_{\rm HI} \ge 0.925 (1 \sigma)$ using Tanh to $\bar{x}_{\rm HI} \ge 0.849 (1 \sigma)$ using the physical model. However, biases in inference using $\tau_e$ are negligible for the Planck 2018 data.

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests.

Dark photon dark matter will resonantly convert into visible photons when the dark photon mass is equal to the plasma frequency of the ambient medium. In cosmological contexts, this transition leads to an extremely efficient, albeit short-lived, heating of the surrounding gas. Existing work in this field has been predominantly focused on understanding the implications of these resonant transitions in the limit that the plasma frequency of the Universe can be treated as being perfectly homogeneous, \ie neglecting inhomogeneities in the electron number density. In this work we focus on the implications of heating from dark photon dark matter in the presence of inhomogeneous structure (which is particularly relevant for dark photons with masses in the range $10^{-15} \, {\rm eV} \, \lesssim m_{A^\prime} \lesssim 10^{-12}$ eV), emphasizing both the importance of inhomogeneous energy injection, as well as the sensitivity of cosmological observations to the inhomogeneities themselves. More specifically, we derive modified constraints on dark photon dark matter from the Ly-$\alpha$ forest, and show that the presence of inhomogeneities allows one to extend constraints to masses outside of the range that would be obtainable in the homogeneous limit, while only slightly relaxing their strength. We then project sensitivity for near-future cosmological surveys that are hoping to measure the 21cm transition in neutral hydrogen prior to reionization, and demonstrate that these experiments will be extremely useful in improving sensitivity to masses near $\sim 10^{-14}$ eV, potentially by several orders of magnitude. Finally, we discuss implications for reionization, early star formation, and late-time $y$-type spectral distortions, and show that probes which are inherently sensitive to the inhomogeneous state of the Universe could resolve signatures unique to the light dark photon...

Even if massive ($10\,M_\odot \lesssim M \lesssim 10^4 M_\odot$) primordial black holes (PBHs) can only account for a small fraction of the dark matter (DM) in the universe, they may still be responsible for a sizable fraction of the coalescence events measured by LIGO/Virgo, and/or act as progenitors of the supermassive black holes (SMBHs) observed already at high redshift ($z\gtrsim 6$). In presence of a dominant, non-PBH DM component, the bounds set by CMB via an altered ionization history are modified. We revisit the cosmological accretion of a DM halo around PBHs via toy models and dedicated numerical simulations, deriving updated CMB bounds which also take into account the last Planck data release. We prove that these constraints dominate over other constraints available in the literature at masses $M\gtrsim 20-50\,M_\odot$ (depending on uncertainty in accretion physics), reaching the level $f_{\rm PBH}<3\times 10^{-9}$ around $M\sim 10^{4}\,M_\odot$. These tight bounds are nonetheless consistent with the hypothesis of a primordial origin of the SMBH massive seeds.

At present, the strongest upper limit on $\sum m_{\nu}$, the sum of neutrino masses, is from cosmological measurements. However, this bound assumes that the neutrinos are stable on cosmological timescales, and is not valid if the neutrino lifetime is less than the age of the universe. In this paper, we explore the cosmological signals of theories in which the neutrinos decay into invisible dark radiation on timescales of order the age of the universe, and determine the bound on the sum of neutrino masses in this scenario. We focus on the case in which the neutrinos decay after becoming non-relativistic. We derive the Boltzmann equations that govern the cosmological evolution of density perturbations in the case of unstable neutrinos, and solve them numerically to determine the effects on the matter power spectrum and lensing of the cosmic microwave background. We find that the results admit a simple analytic understanding. We then use these results to perform a Monte Carlo analysis based on the current data to determine the limit on the sum of neutrino masses as a function of the neutrino lifetime. We show that in the case of decaying neutrinos, values of $\sum m_{\nu}$ as large as 0.9 eV are still allowed by the data. Our results have important implications for laboratory experiments that have been designed to detect neutrino masses, such as KATRIN and KamLAND-ZEN.

We explore the cosmological signals of theories in which the neutrinos decay into invisible dark radiation after becoming non-relativistic. We show that in this scenario, near-future large scale structure measurements from the Euclid satellite, when combined with cosmic microwave background data from Planck, may allow an independent determination of both the lifetime of the neutrinos and the sum of their masses. These parameters can be independently determined because the Euclid data will cover a range of redshifts, allowing the growth of structure over time to be tracked. If neutrinos are stable on cosmological timescales, these observations can improve the lower limit on the neutrino lifetime by seven orders of magnitude, from $\mathcal{O}(10)$ years to $2\times 10^8$ years ($95\%$ C.L.), without significantly affecting the measurement of neutrino mass. On the other hand, if neutrinos decay after becoming non-relativistic but on timescales less than $\mathcal{O}(100)$ million years, these observations may allow, not just the first measurement of the sum of neutrino masses, but also the determination of the neutrino lifetime from cosmology.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

We present a detailed investigation of a sub-dominant oscillating scalar field ('early dark energy', EDE) in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until ${\rm log}_{10}(z_c)\sim3.5$, reaching $\rho_{\rm EDE}(z_c)/\rho_{\rm tot}\sim10$%, and diluting faster than matter afterwards can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES collaboration into agreement. A scalar field potential which scales as $V(\phi) \propto \phi^{2n}$ with $2\lesssim n\lesssim 3.4$ around the minimum is preferred at the 68% confidence level, and the {\em Planck} polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefers a potential which flattens at large field displacements. An MCMC analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the $E$-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $n\approx 2$. This instability and ensuing potentially nonlinear, spatially inhomogenoues, dynamics may provide unique signatures of this scenario.

Early dark energy (EDE) that behaves like a cosmological constant at early times (redshifts $z\gtrsim3000$) and then dilutes away like radiation or faster at later times can solve the Hubble tension. In these models, the sound horizon at decoupling is reduced resulting in a larger value of the Hubble parameter $H_0$ inferred from the cosmic microwave background (CMB). We consider two physical models for this EDE, one involving an oscillating scalar field and another a slowly-rolling field. We perform a detailed calculation of the evolution of perturbations in these models. A Markov Chain Monte Carlo search of the parameter space for the EDE parameters, in conjunction with the standard cosmological parameters, identifies regions in which $H_0$ inferred from {\it Planck} CMB data agrees with the SH0ES local measurement. In these cosmologies, current baryon acoustic oscillation and supernova data are described as successfully as in \LCDM, while the fit to {\it Planck} data is slightly improved. Future CMB and large-scale-structure surveys will further probe this scenario.

This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($\gtrsim10$GV), i.e. independent of the propagation scenario: at $1\sigma$, the diffusion slope $\delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $\Delta_h\sim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($\sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.

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.

The axion is a hypothetical, well-motivated dark-matter particle whose existence would explain the lack of charge-parity violation in the strong interaction. In addition to this original motivation, an `axiverse' of ultra-light axions (ULAs) with masses $10^{-33}\,{\rm eV}\lesssim m_{\rm a}\lesssim 10^{-10}\,{\rm eV}$ also emerges from string theory. Depending on the mass, such a ULA contributes to the dark-matter density, or alternatively, behaves like dark energy. At these masses, ULAs' classical wave-like properties are astronomically manifested, potentially mitigating observational tensions within the $\Lambda$CDM paradigm on local-group scales. ULAs also provide signatures on small scales such as suppression of structure, interference patterns and solitons to distinguish them from heavier dark matter candidates. Through their gravitational imprint, ULAs in the presently allowed parameter space furnish a host of observational tests to target in the next decade, altering standard predictions for microwave background anisotropies, galaxy clustering, Lyman-$\alpha$ absorption by neutral hydrogen along quasar sightlines, pulsar timing, and the black-hole mass spectrum.

We discuss the origin of the anti-helium-3 and -4 events possibly detected by AMS-02. Using up-to-date semi-analytical tools, we show that spallation from primary hydrogen and helium nuclei onto the ISM predicts a $\overline{{}^3{\rm He}}$ flux typically one to two orders of magnitude below the sensitivity of AMS-02 after 5 years, and a $\overline{{}^4{\rm He}}$ flux roughly 5 orders of magnitude below the AMS-02 sensitivity. We argue that dark matter annihilations face similar difficulties in explaining this event. We then entertain the possibility that these events originate from anti-matter-dominated regions in the form of anti-clouds or anti-stars. In the case of anti-clouds, we show how the isotopic ratio of anti-helium nuclei might suggest that BBN has happened in an inhomogeneous manner, resulting in anti-regions with a anti-baryon-to-photon ratio $\bar{\eta}\simeq10^{-3}\eta$. We discuss properties of these regions, as well as relevant constraints on the presence of anti-clouds in our Galaxy. We present constraints from the survival of anti-clouds in the Milky-Way and in the early Universe, as well as from CMB, gamma-ray and cosmic-ray observations. In particular, these require the anti-clouds to be almost free of normal matter. We also discuss an alternative where anti-domains are dominated by surviving anti-stars. We suggest that part of the unindentified sources in the 3FGL catalog can originate from anti-clouds or anti-stars. AMS-02 and GAPS data could further probe this scenario.

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

We perform a comprehensive analysis of the most common early- and late-Universe solutions to the $H_0$, Ly-$\alpha$, and $S_8$ discrepancies. When considered on their own, massive neutrinos provide a natural solution to the $S_8$ discrepancy at the expense of increasing the $H_0$ tension. If all extensions are considered simultaneously, the best-fit solution has a neutrino mass sum of $\sim 0.4$ eV, a dark energy equation of state close to that of a cosmological constant, and no additional relativistic degrees of freedom. However, the $H_0$ tension, while weakened, remains unresolved. Motivated by this result, we perform a non-parametric reconstruction of the evolution of the dark energy fluid density (allowing for negative energy densities), together with massive neutrinos. When all datasets are included, there exists a residual $\sim1.9\sigma$ tension with $H_0$. If this residual tension remains in the future, it will indicate that it is not possible to solve the $H_0$ tension solely with a modification of the late-Universe dynamics within standard general relativity. However, we do find that it is possible to resolve the tension if either galaxy BAO or JLA supernovae data are omitted. We find that \textit{negative} dark energy densities are favored near redshift $z\sim2.35$ when including the Ly-$\alpha$ BAO measurement (at $\sim 2\sigma$). This behavior may point to a negative curvature, but it is most likely indicative of systematics or at least an underestimated covariance matrix. Quite remarkably, we find that in the extended cosmologies considered in this work, the neutrino mass sum is always close to $0.4$ eV regardless of the choice of external datasets, as long as the $H_0$ tension is solved or significantly decreased.

We investigate the hypothesis that Coulomb-type interactions between dark matter (DM) and baryons explain the anomalously low 21cm brightness-temperature minimum at redshift z ~ 17 that was recently measured by the EDGES experiment. In particular, we reassess the validity of the scenario where a small fraction of the total DM is millicharged, focusing on newly derived constraints from Planck 2015 cosmic microwave background (CMB) data. Crucially, the CMB power spectrum is sensitive to DM-baryon scattering if the fraction of interacting DM is larger than (or comparable to) the fractional uncertainty in the baryon energy density. Meanwhile, there is a mass-dependent lower limit on the fraction for which the required interaction to cool the baryons sufficiently is so strong that it drives the interacting-DM temperature to the baryon temperature prior to their decoupling from the CMB. If this occurs as early as recombination, the cooling saturates. We precisely determine the viable parameter space for millicharged DM, and find that only a fraction (m_chi/MeV) 0.0115% <~ f <~ 0.4% of the entire DM content, and only for DM-particle masses between 0.5 MeV - 35 MeV, can be charged at the level needed to marginally explain the anomaly, without violating limits from SLAC, CMB, Big-Bang nucleosynthesis (BBN), or stellar and SN1987A cooling. In reality, though, we demonstrate that at least moderate fine tuning is required to both agree with the measured absorption profile and overcome various astrophysical sources of heating. Finally, we point out that a ~0.4% millicharged DM component which is tightly coupled to the baryons at recombination may resolve the current 2-sigma tension between the BBN and CMB determinations of the baryon energy density. Future CMB-S4 measurements will be able to probe this scenario directly.

We perform an improved cosmic microwave background (CMB) analysis to search for dark matter-proton scattering with a momentum-transfer cross section of the form $\sigma_0 v^n$ for $n=-2$ and $n=-4$. In particular, we present a new and robust prescription for incorporating the relative bulk velocity between the dark matter and baryon fluids into the standard linear Boltzmann calculation. Using an iterative procedure, we self-consistently include the effects of the bulk velocities in a cosmology in which dark matter interacts with baryons. With this prescription, we derive CMB bounds on the cross section, excluding $\sigma_0 > 2.3 \times 10^{-33}~\mathrm{cm}^2$ for $n=-2$ and $\sigma_0 > 1.7 \times 10^{-41}~\mathrm{cm}^2$ for $n=-4$ at $95\%$ confidence, for dark matter masses below 10 MeV. Furthermore, we investigate how these constraints change when only a subcomponent of dark matter is interacting. We show that Planck limits vanish if $\lesssim 0.4\%$ of dark matter is tightly coupled to baryons. We discuss the implications of our results for present and future cosmological observations.

Cosmological observations are used to test for imprints of an ultra-light axion-like field (ULA), with a range of potentials $V(\phi)\propto[1-\cos(\phi/f)]^n$ set by the axion-field value $\phi$ and decay constant $f$. Scalar field dynamics dictate that the field is initially frozen and then begins to oscillate around its minimum when the Hubble parameter drops below some critical value. For $n\!=\!1$, once dynamical, the axion energy density dilutes as matter; for $n\!=\!2$ it dilutes as radiation and for $n\!=\!3$ it dilutes faster than radiation. Both the homogeneous evolution of the ULA and the dynamics of its linear perturbations are included, using an effective fluid approximation generalized from the usual $n=1$ case. ULA models are parameterized by the redshift $z_c$ when the field becomes dynamical, the fractional energy density $f_{z_c} \equiv \Omega_a(z_c)/\Omega_{\rm tot}(z_c)$ in the axion field at $z_c$, and the effective sound speed $c_s^2$. Using Planck, BAO and JLA data, constraints on $f_{z_c}$ are obtained. ULAs are degenerate with dark energy for all three potentials if $1+z_c \lesssim 10$. When $3\times10^4 \gtrsim 1+z_c \gtrsim 10 $, $f_{z_c}$ is constrained to be $ \lesssim 0.004 $ for $n=1$ and $f_{z_c} \lesssim 0.02 $ for the other two potentials. The constraints then relax with increasing $z_c$. These results strongly constrain ULAs as a resolution to cosmological tensions, such as discrepant measurements of the Hubble constant, or the EDGES measurement of the global 21 cm signal.

Using cosmic-ray boron to carbon ratio (B/C) data recently released by the AMS-02 experiment, we find tantalizing indications ({\it decisive evidence}, in Bayesian terms) in favor of a diffusive origin for the broken power-law spectra found in protons ($p$) and helium nuclei (He). The result is robust with respect to currently estimated uncertainties in the cross sections, and in the presence of a small component of primary boron, expected because of spallation at the acceleration site. Reduced errors at high energy as well as further cosmic ray nuclei data (as absolute spectra of C, N, O, Li, Be) may definitively confirm this scenario.

We devise a new user-friendly tool interfaced with the Boltzmann code CLASS to deal with any kind of exotic electromagnetic energy injection in the universe and its impact on anisotropies of the Cosmic Microwave Background. It makes use of the results from standard electromagnetic cascade calculations develop in the context of WIMP annihilation, generalized to incorporate any injection history. We first validate it on a specific WIMP scenario, the Higgs Portal model, confirming that the standard effective on-the-spot treatment is accurate enough. We then analyze the more involved example of evaporating Primordial Black Holes (PBHs) with masses in the range $[3\times10^{13},5\times10^{16}]$g, for which the standard approximations break down. We derive robust CMB bounds on the relic density of evaporating PBHs, ruling out the possibility for PBHs with a monochromatic distribution of masses in the range $[3\times10^{13},2.5\times10^{16}]$g to represent all of the Dark Matter in our Universe. Remarkably, we confirm with an accurate study that the CMB bounds are several orders of magnitude stronger than those from the galactic gamma-ray background in the range $[3\times10^{13},3\times10^{14}]$g. A future CMB experiment like CORE+, or an experiment attempting at measuring the 21 cm signal from the Dark Ages could greatly improve the sensitivity to these models.

Stellar-mass Primordial Black Holes (PBH) have been recently reconsidered as a Dark Matter (DM) candidate, after the aLIGO discovery of several binary BH mergers with masses of tens of $M_\odot$. Matter accretion on such massive objects leads to the emission of high-energy photons, capable of altering the ionization and thermal history of the universe. This in turn affects the statistical properties of the cosmic microwave background (CMB) anisotropies. Previous analyses have assumed spherical accretion. We argue that this approximation likely breaks down and that an accretion disk should form in the dark ages. Using the most up-to-date tools to compute the energy deposition in the medium, we derive constraints on the fraction of DM in PBH. Provided that disks form early on, even under conservative assumptions for accretion, these constraints exclude a monochromatic distribution of PBH with masses above $\sim 2\, M_\odot$ as the dominant form of DM. The bound on the median PBH mass gets more stringent if a broad, log-normal mass function is considered. A deepened understanding of non-linear clustering properties and BH accretion disk physics would permit an improved treatment and possibly lead to more stringent constraints.

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.

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]

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.

Two years ago, the AMS collaboration released the most precise measurement of the cosmic ray positron flux. It confirms that pure secondary predictions fall below the data above 10 GeV, suggesting the presence of a primary component, e.g. annihilations of WIMPs dark matter. Most analyses have focused on the high-energy part of the spectrum, disregarding the GeV energy region where cosmic ray transport is harder to model and solar modulation comes into play. Given the high quality of AMS measurements, it is timely to re-examine the positron anomaly over the entire energy range, taking into account transport processes so far neglected, e.g. convection or diffusive re-acceleration. We devise a new semi-analytical method to take into account transport processes so far neglected, but important below a few GeV. It is based on the pinching of inverse Compton and synchrotron energy losses inside the Galactic disc. It allows to carry out extensive scans over the cosmic ray propagation parameters, which we strongly constrain by requiring that the secondary component does not overshoot the AMS measurements. Only models with large diffusion coefficients survive this test. The positron flux is a powerful and independent probe of cosmic ray propagation, complementary to the boron-to-carbon ratio. We then scan over WIMP mass to fit the annihilation cross section and branching ratios, exploring both direct annihilations into standard model particles or through light mediators. In the former case, the best fit yields a $p$-value of 0.4% for a mass of 264 GeV, a value that does not allow to reproduce the highest energy data points. Worse quality fits are found in the latter case. The interpretation of the positron excess in terms of single DM species annihilations is strongly disfavored. This conclusion is based solely on the positron data, and no other observation needs to be invoked.

The interplay between cosmology and earth based experiments is crucial in order to pin down neutrino physics. Indeed cosmology can provide very tight, yet model dependent, constraints on some neutrino properties. Here we focus on the neutrino mass sum, reviewing the up to date current bounds and showing the results of our forecast of the sensitivity of future experiments. Finally, we discuss the case for sterile neutrinos, explaining how non standard sterile neutrino self-interactions can reconcile the oscillation anomalies with cosmology.

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 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 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 compute cosmic microwave background (CMB) anisotropy constraints on exotic forms of energy injection in electromagnetic (e.m.) channels over a large range of timescales. These constraints are very powerful around or just after recombination, although CMB keeps some sensitivity e.g. to decaying species with lifetimes as long as $10^{25}\,$s. We review here complementary with CMB spectral distortions and primordial nucleosynthesis bounds, which dominate at earlier timescales. For the first time, we describe the effects of the e.m. energy injection on the CMB power spectra as a function of the injection epoch, using the lifetime of a decaying particle as proxy. We identify a suitable on-the-spot approximation. Our results are of interest not only for early universe relics constituting (a fraction of) the dark matter, but also for other exotic injection of e.m. radiation. For illustration, we apply our formalism to: i) Primordial black holes of mass $\in [10^{13.5},10^{16.8}]$ g, showing that the constraints are comparable to the ones obtained from gamma-ray background studies and even dominate below $\sim 10^{14}$g. ii) To a peculiar mass-mixing range in the sterile neutrino parameter space, complementary to other astrophysical and laboratory probes. iii) Finally, we provide a first estimate of the room for improvement left for forthcoming 21 cm experiments, comparing it with the reach of proposed CMB spectral distortion (PiXiE) and CMB angular power spectrum (CORE) missions. We show that the best and most realistic opportunity to look for this signal (or to improve over current constraints) in the 21 cm probe is to focus on the Cosmic Dawn epoch, $15\lesssim z\lesssim30$, where the qualitatively unambiguous signature of a spectrum in emission can be expected for models that evade all current constraints.

Future Cosmic Microwave Background experiments together with upcoming galaxy and 21-cm surveys will provide extremely accurate measurements of different cosmological observables located at different epochs of the cosmic history. The new data will be able to constrain the neutrino mass sum with the best precision ever. In order to exploit the complementarity of the different redshift probes, a deep understanding of the physical effects driving the impact of massive neutrinos on CMB and large scale structures is required. The goal of this work is to describe these effects, assuming a summed neutrino mass close to its minimum allowed value. We find that parameter degeneracies can be removed by appropriate combinations, leading to robust and model independent constraints. A joint forecast of the sensitivity of Euclid and DESI surveys together with a CORE-like CMB experiment leads to a $1\sigma$ uncertainty of $14$~meV on the summed neutrino mass. However this particular combination gives rise to a peculiar degeneracy between $M_\nu$ and the optical depth at reionization. Independent constraints from 21-cm surveys can break this degeneracy and decrease the $1\sigma$ uncertainty down to $12$~meV.

We consider a cosmological model in which a fraction $f$ of the Dark Matter (DM) is allowed to decay in an invisible relativistic component, and compute the resulting constraints on both the decay width (or inverse lifetime) $\Gamma$ and $f$ from purely gravitational arguments. We report a full derivation of the Boltzmann hierarchy, correcting a mistake in previous literature, and compute the impact of the decay --as a function of the lifetime-- on the CMB and matter power spectra. From CMB only, we obtain that no more than 3.8 % of the DM could have decayed in the time between recombination and today (all bounds quoted at 95 % CL). We also comment on the important application of this bound to the case where primordial black holes constitute DM, a scenario notoriously difficult to constrain. For lifetimes longer than the age of the Universe, the bounds can be cast as $f\Gamma < 6.3\times10^{-3}$ Gyr$^{-1}$. For the first time, we also checked that degeneracies with massive neutrinos are broken when information from the large scale structure is used. Even secondary effects like CMB lensing suffice to this purpose. Decaying DM models have been invoked to solve a possible tension between low redshift astronomical measurements of $\sigma_8$ and $\Omega_{\rm m}$ and the ones inferred by Planck. We reassess this claim finding that with the most recent BAO, HST and $\sigma_8$ data extracted from the CFHT survey, the tension is only slightly reduced despite the two additional free parameters, loosening the bound to $f\Gamma < 15.9\times10^{-3}$ Gyr$^{-1}$. The bound however improves to $f\Gamma < 5.9\times10^{-3}$ Gyr$^{-1}$ if only data consistent with the CMB are included. This highlights the importance of establishing whether the tension is due to real physical effects or unaccounted systematics, for settling the reach of achievable constraints on decaying DM.

It is well known that annihilations in the homogeneous fluid of dark matter (DM) can leave imprints in the cosmic microwave background (CMB) anisotropy power spectrum. However, the relevance of DM annihilations in halos for cosmological observables is still subject to debate, with previous works reaching different conclusions on this point. Also, all previous studies used a single type of parameterization for the astrophysical reionization, and included no astrophysical source for the heating of the intergalactic medium. In this work, we revisit these problems. When standard approaches are adopted, we find that the ionization fraction does exhibit a very particular (and potentially constraining) pattern, but the currently measurable optical depth to reionization is left almost unchanged: In agreement with the most of the previous literature, for plausible halo models we find that the modification of the signal with respect to the one coming from annihilations in the smooth background is tiny, below cosmic variance within currently allowed parameter space. However, if different and probably more realistic treatments of the astrophysical sources of reionization and heating are adopted, a more pronounced effect of the DM annihilation in halos is possible. We thus conclude that within currently adopted baseline models the impact of the virialised DM structures cannot be uncovered by CMB power spectra measurements, but a larger impact is possible if peculiar models are invoked for the redshift evolution of the DM annihilation signal or different assumptions are made for the astrophysical contributions. A better understanding (both theoretical and observational) of the reionization and temperature history of the universe, notably via the 21 cm signal, seems the most promising way for using halo formation as a tool in DM searches, improving over the sensitivity of current cosmological probes.

Using the updated proton and helium fluxes just released by the AMS-02 experiment we reevaluate the secondary astrophysical antiproton to proton ratio and its uncertainties, and compare it with the ratio preliminarly reported by AMS-02. We find no unambiguous evidence for a significant excess with respect to expectations. Yet, some preference for a flatter energy dependence of the diffusion coefficient starts to emerge. Also, we provide a first assessment of the room left for exotic components such as Galactic Dark Matter annihilation or decay, deriving new stringent constraints.

In Poulin and Serpico [Phys. Rev. Lett. 114, 091101 (2015)] we have recently argued that when the energy of a photon injected in the primordial plasma falls below the pair-production threshold, the universality of the non-thermal photon spectrum from the standard theory of electromagnetic cascades onto a photon background breaks down. We showed that this could reopen or widen the parameter space for an exotic solution to the 'lithium problem'. Here we discuss another application, namely the impact that this has on non-thermal big bang nucleosynthesis constraints from 4He, 3He and 2H, using the parametric example of monochromatic photon injection of different energies. Typically, we find tighter bounds than those existing in the literature, up to more than one order of magnitude. As a consequence of the non-universality of the spectrum, the energy-dependence of the photodissociation cross-sections is important. We also compare the constraints obtained with current level and future reach of cosmic microwave background spectral distortion bounds.

The standard theory of electromagnetic cascades onto a photon background predicts a quasi-universal shape for the resulting non-thermal photon spectrum. This has been applied to very disparate fields, including non-thermal big bang nucleosynthesis (BBN). However, once the energy of the injected photons falls below the pair-production threshold the spectral shape is very different, a fact that has been overlooked in past literature. This loophole may have important phenomenological consequences, since it generically alters the BBN bounds on non-thermal relics: for instance it allows to re-open the possibility of purely electromagnetic solutions to the so-called "cosmological lithium problem", which were thought to be excluded by other cosmological constraints. We show this with a proof-of-principle example and a simple particle physics model, compared with previous literature.

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

The positron fraction in cosmic rays was found to be a steadily increasing in function of energy, above $\sim$ 10 GeV. This behaviour contradicts standard astrophysical mechanisms, in which positrons are secondary particles, produced in the interactions of primary cosmic rays during the propagation in the interstellar medium. The observed anomaly in the positron fraction triggered a lot of excitement, as it could be interpreted as an indirect signature of the presence of dark matter species in the Galaxy. Alternatively, it could be produced by nearby astrophysical sources, such as pulsars. Both hypotheses are probed in this work in light of the latest AMS-02 positron fraction measurements. The transport of the primary and secondary positrons in the Galaxy is described using a semi-analytic two-zone model. MicrOMEGAs is used to model the positron flux generated by dark matter species. The description of the positron fraction from astrophysical sources is based on the pulsar observations included in the ATNF catalogue. We find that the mass of the favoured dark matter candidates is always larger than 500 GeV. The only dark matter species that fulfils the numerous gamma ray and cosmic microwave background bounds is a particle annihilating into four leptons through a light scalar or vector mediator, with a mixture of tau (75%) and electron (25%) channels, and a mass between 0.5 and 1 TeV. The positron anomaly can also be explained by a single astrophysical source and a list of five pulsars from the ATNF catalogue is given. Those results are obtained with the cosmic ray transport parameters that best fit the B/C ratio. Uncertainties in the propagation parameters turn out to be very significant. In the WIMP annihilation cross section to mass plane for instance, they overshadow the error contours derived from the positron data.