Type Ia Supernovae (SNeIa) used as standardizable candles have been instrumental in the discovery of cosmic acceleration, usually attributed to some form of dark energy (DE). Recent studies have raised the issue of whether intrinsic SNeIa luminosities might evolve with redshift. While the evidence for cosmic acceleration is robust to this possible systematic, the question remains of how much the latter can affect the inferred properties of the DE component responsible for cosmic acceleration. This is the question we address in this work. We use SNeIa distance moduli measurements from the Pantheon and JLA samples. We consider models where the DE equation of state is a free parameter, either constant or time-varying, as well as models where DE and dark matter interact, and finally a model-agnostic parametrization of effects due to modified gravity (MG). When SNeIa data are combined with Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements, we find strong degeneracies between parameters governing the SNeIa systematics, the DE parameters, and the Hubble constant $H_0$. These degeneracies significantly broaden the DE parameter uncertainties, in some cases leading to ${\cal O}(\sigma)$ shifts in the central values. However, including low-redshift Baryon Acoustic Oscillation and Cosmic Chronometer measurements, as well as CMB lensing measurements, considerably improves the previous constraints, and the only remaining effect of the examined systematic is a $\lesssim 40\%$ broadening of the uncertainties on the DE parameters. The constraints we derive on the MG parameters are instead basically unaffected by the systematic in question. We therefore confirm the overall soundness of dark energy properties.

Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of $H_0$ measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high $\ell$ polarisation data constrain the sterile neutrino mass to be less than approximately 1~eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in a narrow range around 1~eV remains consistent with both cosmology and SBL data.

We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Namely, we have included all T2K measurements as of December 2019, the most recent NO$\nu$A antineutrino statistics, and data collected by the Daya Bay and RENO reactor experiments. Concerning the atmospheric and solar sectors, we have also updated our analyses of DeepCore and SNO data, respectively. All in all, these new analyses result in more accurate measurements of $\theta_{13}$, $\theta_{12}$, $\Delta m_{21}^2$ and $|\Delta m_{31}^2|$. The best fit value for the atmospheric angle $\theta_{23}$ lies in the second octant, but first octant solutions remain allowed at $\sim2\sigma$. Regarding CP violation measurements, the preferred value of $\delta$ we obtain is 1.20$\pi$ (1.54$\pi$) for normal (inverted) neutrino mass ordering. These new results should be regarded as extremely robust due to the excellent agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a preference for the normal neutrino mass ordering at the $2.7\sigma$ level. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the significance to $3.1\sigma$ within a conservative approach. A more aggressive data set combination of cosmological observations leads to a stronger preference for normal with respect to inverted mass ordering, at the $3.3\sigma$ level. This cosmological data set provides $2\sigma$ upper limits on the total neutrino mass corresponding to $\sum m_\nu<0.13$ ($0.15$)~eV in the normal (inverted) neutrino mass ordering scenario. These bounds are among the most complete ones in the literature, as they include all currently available neutrino physics inputs.

One of the most criticized features of Bayesian statistics is the fact that credible intervals, especially when open likelihoods are involved, may strongly depend on the prior shape and range. Many analyses involving open likelihoods are affected by the eternal dilemma of choosing between linear and logarithmic prior, and in particular in the latter case the situation is worsened by the dependence on the prior range under consideration. In this letter, we revive a simple method to obtain constraints that depend neither on the prior shape nor range and, using the tools of Bayesian model comparison, extend it to overcome the possible dependence of the bounds on the choice of free parameters in the numerical analysis. An application to the case of cosmological bounds on the sum of the neutrino masses is discussed as an example.

We provide a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, $\beta$-decay and neutrinoless double-$\beta$ decay ($0\nu\beta\beta$) experiments. We directly constrain the full $3+1$ active-sterile mixing matrix elements $|U_{\alpha4}|^2$, with $\alpha \in ( e,\mu ,\tau )$, and the mass-squared splitting $\Delta m^2_{41} \equiv m_4^2-m_1^2$. We find that results for a $3+1$ case differ from previously studied $1+1$ scenarios where the sterile is only coupled to one of the neutrinos, which is largely explained by parameter space volume effects. Limits on the mass splitting and the mixing matrix elements are currently dominated by the cosmological data sets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below $|U_{\alpha4}|^2 \lesssim 10^{-3}$. Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos are in serious tension with these bounds, irrespective of prior assumptions. We also translate the bounds from the cosmological analysis into constraints on the parameters probed by laboratory searches, such as $m_\beta$ or $m_{\beta \beta}$, the effective mass parameters probed by $\beta$-decay and $0\nu\beta\beta$ searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of $m_\beta < 0.09$ eV and $m_{\beta \beta} < 0.07$ eV at 95\% C.L, more stringent than the limits from current laboratory experiments.

The Cosmic Neutrino Background is a prediction of the standard cosmological model, but it has been never observed directly. Experiments with the aim of detecting relic CNB neutrinos are under development. For such experiments, the expected event rate depends on the local number density of relic neutrinos. Since massive neutrinos can be attracted by the gravitational potential of our galaxy and cluster locally, a local overdensity of relic neutrinos should exist at Earth. We report the status of our knowledge on the clustering of neutrinos and the consequences for future direct detection experiments.

The standard cosmological model predicts the existence of a Cosmic Neutrino Background, which has not yet been observed directly. Some experiments aiming at its detection are currently under development, despite the tiny kinetic energy of the cosmological relic neutrinos, which makes this task incredibly challenging. Since massive neutrinos are attracted by the gravitational potential of our Galaxy, they can cluster locally. Neutrinos should be more abundant at the Earth position than at an average point in the Universe. This fact may enhance the expected event rate in any future experiment. Past calculations of the local neutrino clustering factor only considered a spherical distribution of matter in the Milky Way and neglected the influence of other nearby objects like the Virgo cluster, although recent $N$-body simulations suggest that the latter may actually be important. In this paper, we adopt a back-tracking technique, well established in the calculation of cosmic rays fluxes, to perform the first three-dimensional calculation of the number density of relic neutrinos at the Solar System, taking into account not only the matter composition of the Milky Way, but also the contribution of the Andromeda galaxy and the Virgo cluster. The effect of Virgo is indeed found to be relevant and to depend non-trivially on the value of the neutrino mass. Our results show that the local neutrino density is enhanced by 0.53% for a neutrino mass of 10 meV, 12% for 50 meV, 50% for 100 meV or 500% for 300 meV.

Light sterile neutrinos with a mass around 1 eV have been studied for many years as a possible explanation of the so called short-baseline neutrino oscillation anomalies. Recently, several neutrino oscillation experiments reported preferences for non-zero values of the mixing angles and squared mass differences for active-sterile mixing, which however are not always in agreement. I review our current knowledge of the light sterile neutrino in the 3+1 model, starting with a separate discussion on the status of the most relevant searches and then analyzing the problems that arise when combining different probes in a global fit. A short summary on the tension with cosmological observations is also provided.

The Magnificent CE$\nu$NS Workshop (2018) was held November 2 & 3 of 2018 on the University of Chicago campus and brought together theorists, phenomenologists, and experimentalists working in numerous areas but sharing a common interest in the process of coherent elastic neutrino-nucleus scattering (CE$\nu$NS). This is a collection of abstract-like summaries of the talks given at the meeting, including links to the slides presented. This document and the slides from the meeting provide an overview of the field and a snapshot of the robust CE$\nu$NS-related efforts both planned and underway.

The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.

In the framework of a 3+1 scheme with an additional inert state, we consider the thermalisation of sterile neutrinos in the early Universe taking into account the full $4\times4$ mixing matrix. The evolution of the neutrino energy distributions is found solving the momentum-dependent kinetic equations with full diagonal collision terms, as in previous analyses of flavour neutrino decoupling in the standard case. The degree of thermalisation of the sterile state is shown in terms of the effective number of neutrinos, $N_{\rm eff}$, and its dependence on the three additional mixing angles ($\theta_{14}$, $\theta_{24}$, $\theta_{34}$) and on the squared mass difference $\Delta m^2_{41}$ is discussed. Our results are relevant for fixing the contribution of a fourth light neutrino species to the cosmological energy density, whose value is very well constrained by the final Planck analysis. For the preferred region of active-sterile mixing parameters from short-baseline neutrino experiments, we find that the fourth state is fully thermalised ($N_{\rm eff}\simeq 4$).

We propose here a \emph{novel} method which singles out the \emph{a priori} unavoidable dependence on the underlying cosmological model when extracting parameter constraints, providing robust limits which only depend on the considered dataset. Interestingly, when dealing with several possible cosmologies and interpreting the Bayesian preference in terms of the Gaussian statistical evidence, the preferred model is much less favored than when only two cases are compared. As a working example, we apply our approach to the cosmological neutrino mass bounds, which play a fundamental role not only in establishing the contribution of relic neutrinos to the dark matter of the Universe, but also in the planning of future experimental searches of the neutrino character and of the neutrino mass ordering.

In this talk I will review the current status of the constraints on the neutrino properties from cosmological measurements, with a particular focus on their mass and effective number. I will also discuss the existing tensions within the context of the \lcdm\ model, including the discrepancies on the Hubble parameter and on the matter fluctuations at small scales, and how neutrinos could help to alleviate the aforementioned problems.

We present a detailed description of the electromagnetic filter for the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Starting with an initial estimate for the orbital magnetic moment, the higher-order drift process of ExB is configured to balance the gradient-B drift motion of the electron in such a way as to guide the trajectory into the standing voltage potential along the mid-plane of the filter. As a function of drift distance along the length of the filter, the filter zooms in with exponentially increasing precision on the transverse velocity component of the electron kinetic energy. This yields a linear dimension for the total filter length that is exceptionally compact compared to previous techniques for electromagnetic filtering. The parallel velocity component of the electron kinetic energy oscillates in an electrostatic harmonic trap as the electron drifts along the length of the filter. An analysis of the phase-space volume conservation validates the expected behavior of the filter from the adiabatic invariance of the orbital magnetic moment and energy conservation following Liouville's theorem for Hamiltonian systems.

The ordering of the neutrino masses is a crucial input for a deep understanding of flavor physics, and its determination may provide the key to establish the relationship among the lepton masses and mixings and their analogous properties in the quark sector. The extraction of the neutrino mass ordering is a data-driven field expected to evolve very rapidly in the next decade. In this review, we both analyze the present status and describe the physics of subsequent prospects. Firstly, the different current available tools to measure the neutrino mass ordering are described. Namely, reactor, long-baseline (accelerator and atmospheric) neutrino beams, laboratory searches for beta and neutrinoless double beta decays and observations of the cosmic background radiation and the large scale structure of the universe are carefully reviewed. Secondly, the results from an up-to-date comprehensive global fit are reported: the Bayesian analysis to the 2018 publicly available oscillation and cosmological data sets provides \emph{strong} evidence for the normal neutrino mass ordering versus the inverted scenario, with a significance of 3.5 standard deviations. This preference for the normal neutrino mass ordering is mostly due to neutrino oscillation measurements. Finally, we shall also emphasize the future perspectives for unveiling the neutrino mass ordering. In this regard, apart from describing the expectations from the aforementioned probes, we also focus on those arising from alternative and novel methods, as 21~cm cosmology, core-collapse supernova neutrinos and the direct detection of relic neutrinos.

We propose to achieve the proof-of-principle of the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Each of the technological challenges described in [1,2] will be targeted and hopefully solved by the use of the latest experimental developments and profiting from the low background environment provided by the LNGS underground site. The first phase will focus on the graphene technology for a tritium target and the demonstration of TES microcalorimetry with an energy resolution of better than 0.05 eV for low energy electrons. These technologies will be evaluated using the PTOLEMY prototype, proposed for underground installation, using precision HV controls to step down the kinematic energy of endpoint electrons to match the calorimeter dynamic range and rate capabilities. The second phase will produce a novel implementation of the EM filter that is scalable to the full target size and which demonstrates intrinsic triggering capability for selecting endpoint electrons. Concurrent with the CNB program, we plan to exploit and develop the unique properties of graphene to implement an intermediate program for direct directional detection of MeV dark matter [3,4]. This program will evaluate the radio-purity and scalability of the graphene fabrication process with the goal of using recently identified ultra-high radio-purity CO2 sources. The direct detection of the CNB is a snapshot of early universe dynamics recorded by the thermal relic neutrino yield taken at a time that predates the epochs of Big Bang Nucleosynthesis, the Cosmic Microwave Background and the recession of galaxies (Hubble Expansion). Big Bang neutrinos are believed to have a central role in the evolution of the Universe and a direct measurement with PTOLEMY will unequivocally establish the extent to which these predictions match present-day neutrino densities.

One of the not-yet determined properties of neutrinos is the ordering of their mass eigenstates. We combine the available data from neutrino oscillations, neutrinoless double beta decay and Cosmic Microwave Background observations to derive robust constraints on the mass ordering in a Bayesian context. Based on arxiv:1801.04946.

We examine the results from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES), which has recently claimed the detection of a strong absorption in the 21 cm hyperfine transition line of neutral hydrogen, at redshifts demarcating the early stages of star formation. More concretely, we study the compatibility of the shape of the EDGES absorption profile, centered at a redshift of $z \sim 17.2$, with measurements of the reionization optical depth, the Gunn-Peterson optical depth, and Lyman-$\alpha$ emission from star-forming galaxies, for a variety of possible reionization models within the standard $\Lambda$CDM framework (that is, a Universe with a cosmological constant and cold dark matter). When, conservatively, we only try to accommodate the location of the absorption dip, we identify a region in the parameter space of the astrophysical parameters that successfully explains all of the aforementioned observations. However, one of the most abnormal features of the EDGES measurement is the absorption amplitude, which is roughly a factor of two larger than the maximum allowed value in the $\Lambda$CDM framework. We point out that the simple considered astrophysical models that produce the largest absorption amplitudes are unable to explain the depth of the dip and of reproducing the observed shape of the absorption profile.

A deep understanding of the Epoch of Reionization is still missing in our knowledge of the universe. While future probes will allow us to test the precise evolution of the free electron fraction from redshifts between $z\simeq 6$ and $z\simeq 20$, at present one could ask what kind of reionization processes are allowed by present Cosmic Microwave Background temperature and polarization measurements. An early contribution to reionization could imply a departure from the standard picture where star formation determines the reionization onset. BBy considering a broad class of possible reionization parameterizations, we find that current data do not require an early reionization component in our universe and that only one marginal class of models, based on a particular realization of reionization, may point to that. In addition, the frequentist Akaike Information Criterion (AIC) provides strong evidence against alternative reionization histories, favoring the most simple reionization scenario, which describes reionization by means of only one (constant) reionization optical depth $\tau$.

We present a Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and CMB observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and c) find that there is a weak-to-moderate preference for logarithmic priors. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass. We also show that the limits on the total neutrino mass $\sum m_\nu$ can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.

We consider the phenomenological implications of the violation of the Pauli exclusion principle for neutrinos, focusing on cosmological observables such as the spectrum of Cosmic Microwave Background anisotropies, Baryon Acoustic Oscillations and the primordial abundances of light elements. Neutrinos that behave (at least partly) as bosonic particles have a modified equilibrium distribution function that implies a different influence on the evolution of the Universe that, in the case of massive neutrinos, can not be simply parametrized by a change in the effective number of neutrinos. Our results show that, despite the precision of the available cosmological data, only very weak bounds can be obtained on neutrino statistics, disfavouring a more bosonic behaviour at less than $2\sigma$.

The Cosmic Neutrino Background is a prediction of the standard cosmological model, but it has been never observed directly. In the experiments with the aim of detecting relic CNB neutrinos, currently under development, the expected event rate depends on the local density of relic neutrinos. Since massive neutrinos can be attracted by the gravitational potential of our galaxy and cluster around it, a local overdensity of cosmic neutrinos should exist. Considering the minimal masses guaranteed by neutrino oscillations, we review the computation of the local density of relic neutrinos and we present realistic prospects for a PTOLEMY-like experiment.

Nonzero neutrino masses are required by the existence of flavour oscillations, with values of the order of at least 50 meV. We consider the gravitational clustering of relic neutrinos within the Milky Way, and used the $N$-one-body simulation technique to compute their density enhancement factor in the neighbourhood of the Earth with respect to the average cosmic density. Compared to previous similar studies, we pushed the simulation down to smaller neutrino masses, and included an improved treatment of the baryonic and dark matter distributions in the Milky Way. Our results are important for future experiments aiming at detecting the cosmic neutrino background, such as the Princeton Tritium Observatory for Light, Early-universe, Massive-neutrino Yield (PTOLEMY) proposal. We calculate the impact of neutrino clustering in the Milky Way on the expected event rate for a PTOLEMY-like experiment. We find that the effect of clustering remains negligible for the minimal normal hierarchy scenario, while it enhances the event rate by 10 to 20% (resp. a factor 1.7 to 2.5) for the minimal inverted hierarchy scenario (resp. a degenerate scenario with 150 meV masses). Finally we compute the impact on the event rate of a possible fourth sterile neutrino with a mass of 1.3 eV.

Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold thermal relic. Considering the cases where the non-cold dark matter species could be either a fermion or a boson, we derive consistent upper limits on the non-cold dark relic energy density for a very large range of velocity dispersions, covering the entire range from dark radiation to cold dark matter. To this end, we employ the latest Planck Cosmic Microwave Background data, the recent BOSS DR11 and other Baryon Acoustic Oscillation measurements, and also constraints on the number of Milky Way satellites, the latter of which provides a measure of the suppression of the matter power spectrum at the smallest scales due to the free-streaming of the non-cold dark matter component. We present the results on the fraction $f_{\rm ncdm}$ of non-cold dark matter with respect to the total dark matter for different ranges of the non-cold dark matter masses. We find that the 2$\sigma$ limits for non-cold dark matter particles with masses in the range 1--10 keV are $f_{\rm ncdm}\leq0.29$ (0.23) for fermions (bosons), and for masses in the 10--100 keV range they are $f_{\rm ncdm}\leq0.43$ (0.45), respectively.

Recent analysis of the WMAP and Planck data have shown the presence of a dip and a bump in the spectrum of primordial perturbations at the scales $k=0.002$ Mpc${}^{-1}$ and $k=0.0035$ Mpc${}^{-1}$ respectively. We analyze for the first time the effects a local feature in the inflaton potential to explain the observed deviations from scale invariance in the primordial spectrum. We perform a best fit analysis of the cosmic microwave background (CMB) radiation temperature and polarization data. The effects of the features can improve the agreement with observational data respect to the featureless model. The best fit local feature affects the primordial curvature spectrum mainly in the region of the bump, leaving the spectrum unaffected on other scales.

We explore an extended cosmological scenario where the dark matter is an admixture of cold and additional non-cold species. The mass and temperature of the non-cold dark matter particles are extracted from a number of cosmological measurements. Among others, we consider tomographic weak lensing data and Milky Way dwarf satellite galaxy counts. We also study the potential of these scenarios in alleviating the existing tensions between local measurements and Cosmic Microwave Background (CMB) estimates of the $S_8$ parameter, with $S_8=\sigma_8\sqrt{\Omega_m}$, and of the Hubble constant $H_0$. In principle, a sub-dominant, non-cold dark matter particle with a mass $m_X\sim$~keV, could achieve the goals above. However, the preferred ranges for its temperature and its mass are different when extracted from weak lensing observations and from Milky Way dwarf satellite galaxy counts, since these two measurements require suppressions of the matter power spectrum at different scales. Therefore, solving simultaneously the CMB-weak lensing tensions and the small scale crisis in the standard cold dark matter picture via only one non-cold dark matter component seems to be challenging.

Current measurements of the temperature and polarization anisotropy power spectra of the Cosmic Microwave Background (CMB) seem to indicate that the naive expectation for the slow-roll hierarchy within the most simple inflationary paradigm may not be respected in nature. We show that a primordial power spectra with localized features could in principle give rise to the observed slow-roll anarchy when fitted to a featureless power spectrum. Future CMB missions have the key to disentangle among the two possible paradigms and firmly establish the slow-roll mechanism as the responsible one for the inflationary period in the early universe. From a model comparison perspective, and assuming that nature has chosen a featureless primordial power spectrum, we find that, while with mock Planck data there is only weak evidence against a model with localized features, upcoming CMB measurements may provide strong evidence against such a non-standard primordial power spectrum.

The Short BaseLine (SBL) neutrino oscillations anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such neutrino is highly incompatible with the cosmological data, in particular from the Cosmic Microwave Background (CMB), if no new physics is assumed. An interesting possibility for reconciling the 1 eV sterile neutrino presence in cosmology is related to the existence of a new pseudoscalar interaction. If the sterile neutrinos experience such a pseudoscalar interaction, the cosmological analyses of the full CMB data prefer a sterile neutrino mass that is fully compatible with the SBL determinations. The additional interaction allows to obtain also an improved compatibility of the cosmological predictions with the local measurements of the Hubble parameter.

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

In this paper, we extend an alternative, phenomenological approach to inflation by means of an equation of state and a sound speed, both of them functions of the number of $e$-folds and three phenomenological parameters. This approach captures a number of possible inflationary models, including those with non-canonical kinetic terms or scale-dependent non-gaussianities. We perform Markov Chain Monte Carlo analyses using the latest cosmological publicly available measurements, which include Cosmic Microwave Background (CMB) data from the Planck satellite. In the context of the phenomenological description studied here, the sound speed of the primordial curvature perturbation is found to be $c_s> 0.2$ at $95\%$~CL for a number of $e$-folds $N=60$. Within this parametrization, we discard scale invariance with a significance of about $10\sigma$, and the running of the spectral index is constrained as $\alpha_s=-0.62\,^{+0.08}_{-0.09} \times 10^{-3}$. The limit on the tensor-to-scalar ratio is $r<0.003$ at $95\%$~CL from CMB data alone. The maximum amplitude of the equilateral non-gaussianity that we obtain, $|f^{\text{equil}}_{\text{NL}}|\sim 2$, is much smaller than the current Planck mission errors. Future high-redshift, all-sky surveys could reach the required accuracy on equilateral non-gaussianities to allow for additional and independent tests of the parametrization explored here.

The simplest inflationary models predict a primordial power spectrum (PPS) of the curvature fluctuations that can be described by a power-law function that is nearly scale-invariant. It has been shown, however, that the low-multipole spectrum of the CMB anisotropies may hint the presence of some features in the shape of the scalar PPS, which could deviate from its canonical power-law form. We study the possible degeneracies of this non-standard PPS with the neutrino anisotropies, the neutrino masses, the effective number of relativistic species and a sterile neutrino or a thermal axion mass. The limits on these additional parameters are less constraining in a model with a non-standard PPS when only including the temperature auto-correlation spectrum measurements in the data analyses. The inclusion of the polarization spectra noticeably helps in reducing the degeneracies, leading to results that typically show no deviation from the $\Lambda$CDM model with a standard power-law PPS.

We investigate a phenomenological non-gravitational coupling between dark energy and dark matter, where the interaction in the dark sector is parameterized as an energy transfer either from dark matter to dark energy or the opposite. The models are constrained by a whole host of updated cosmological data: cosmic microwave background temperature anisotropies and polarization, high-redshift supernovae, baryon acoustic oscillations, redshift space distortions and gravitational lensing. Both models are found to be compatible with all cosmological observables, but in the case where dark matter decays into dark energy, the tension with the independent determinations of $H_0$ and $\sigma_8$, already present for standard cosmology, increases: this model in fact predicts lower $H_0$ and higher $\sigma_8$, mostly as a consequence of the higher amount of dark matter at early times, leading to a stronger clustering during the evolution. Instead, when dark matter is fed by dark energy, the reconstructed values of $H_0$ and $\sigma_8$ nicely agree with their local determinations, with a full reconciliation between high- and low-redshift observations. A non-zero coupling between dark energy and dark matter, with an energy flow from the former to the latter, appears therefore to be in better agreement with cosmological data.

In this Thesis I discuss several recent results obtained using the CMB spectra measured by Planck and several other cosmological probes. Extensions of the $\Lambda$CDM model are studied, including the presence of an additional sterile neutrino (motivated by the short-baseline oscillation anomalies) and of a thermal axion. The degeneracies of the cosmological effects of these particles with the power spectrum of primordial perturbations are tested. We also show that the power spectrum of initial scalar perturbations can be degenerate with the presence of primordial non-Gaussianities, thus affecting the constraints on the non-Gaussianity parameter $f_{NL}$. Finally, an effective interaction between dark energy and dark matter is studied.

We perform a cosmological analysis in which we allow the primordial power spectrum of scalar perturbations to assume a shape that is different with respect to the usual power-law, arising from the simplest models of cosmological inflation. We parametrize the primordial power spectrum with a piecewise monotone cubic Hermite function and we use it to investigate how the constraints on the various cosmological parameters change: we find that the obtained limits are relaxed with respect to the power-law case, if CMB polarization data are not included. Moreover, the cosmological analyses provide us some indications about the shape of the reconstructed primordial power spectrum, where we notice possible features around $k\simeq0.002\,\mathrm{Mpc}^{-1}$ and $k\simeq0.0035\,\mathrm{Mpc}^{-1}$. If confirmed in future analyses involving enhanced experimental data, these features suggests that the simplest cosmological inflation models may be incomplete.

We briefly discuss the main effects of the presence of a Light Sterile Neutrino (LS$\nu$) in Cosmology and how its properties can be constrained by Cosmic Microwave Background (CMB) and other cosmological measurements.

The theory and phenomenology of light sterile neutrinos at the eV mass scale is reviewed. The reactor, Gallium and LSND anomalies are briefly described and interpreted as indications of the existence of short-baseline oscillations which require the existence of light sterile neutrinos. The global fits of short-baseline oscillation data in 3+1 and 3+2 schemes are discussed, together with the implications for beta-decay and neutrinoless double-beta decay. The cosmological effects of light sterile neutrinos are briefly reviewed and the implications of existing cosmological data are discussed. The review concludes with a summary of future perspectives.

A relaxed primordial power spectrum (PPS) of scalar perturbations arising from inflation can impact the dark radiation constraints obtained from Cosmic Microwave Background and other cosmological measurements. If inflation produces a non-standard PPS for the initial fluctuations, a fully thermalized light sterile neutrino can be favoured by CMB observations, instead of being strongly disfavoured. In the case of a thermal axion, the constraints on the axion mass are relaxed when the PPS is different from the standard power law.

The simplest models of inflation predict small non-gaussianities and a featureless power spectrum. However, there exist a large number of well-motivated theoretical scenarios in which large non-gaussianties could be generated. In general, in these scenarios the primordial power spectrum will deviate from its standard power law shape. We study, in a model-independent manner, the constraints from future large scale structure surveys on the local non-gaussianity parameter $f_{\rm NL}$ when the standard power law assumption for the primordial power spectrum is relaxed. If the analyses are restricted to the large scale-dependent bias induced in the linear matter power spectrum by non-gaussianites, the errors on the $f_{\rm NL}$ parameter could be increased by $60\%$ when exploiting data from the future DESI survey, if dealing with only one possible dark matter tracer. In the same context, a nontrivial bias $|\delta f_{\rm NL}| \sim 2.5$ could be induced if future data are fitted to the wrong primordial power spectrum. Combining all the possible DESI objects slightly ameliorates the problem, as the forecasted errors on $f_{\rm NL}$ would be degraded by $40\%$ when relaxing the assumptions concerning the primordial power spectrum shape. Also the shift on the non-gaussianity parameter is reduced in this case, $|\delta f_{\rm NL}| \sim 1.6$. The addition of Cosmic Microwave Background priors ensure robust future $f_{\rm NL}$ bounds, as the forecasted errors obtained including these measurements are almost independent on the primordial power spectrum features, and $|\delta f_{\rm NL}| \sim 0.2$, close to the standard single-field slow-roll paradigm prediction.

We present cosmological bounds on the thermal axion mass in an extended cosmological scenario in which the primordial power spectrum of scalar perturbations differs from the usual power-law shape predicted by the simplest inflationary models. The power spectrum is instead modeled by means of a "piecewise cubic Hermite interpolating polynomial" (PCHIP). When using Cosmic Microwave Background measurements combined with other cosmological data sets, the thermal axion mass constraints are degraded only slightly. The addition of the measurements of $\sigma_8$ and $\Omega_m$ from the 2013 Planck cluster catalogue on galaxy number counts relaxes the bounds on the thermal axion mass, mildly favouring a $\sim 1$~eV axion mass, regardless of the model adopted for the primordial power spectrum.However, in general, such a preference disappears if the sum of the three active neutrino masses is also considered as a free parameter in our numerical analyses, due to the strong correlation between the masses of these two hot thermal relics.

We perform a cosmological analysis in which we allow the primordial power spectrum of scalar perturbations to assume a shape that is different from the usual power-law predicted by the simplest models of cosmological inflation. We parameterize the free primordial power spectrum with a "piecewise cubic Hermite interpolating polynomial" (PCHIP). We consider a 3+1 neutrino mixing model with a sterile neutrino having a mass at the eV scale, which can explain the anomalies observed in short-baseline neutrino oscillation experiments. We find that the freedom of the primordial power spectrum allows to reconcile the cosmological data with a fully thermalized sterile neutrino in the early Universe. Moreover, the cosmological analysis gives us some information on the shape of the primordial power spectrum, which presents a feature around the wavenumber $k=0.002\,\text{Mpc}^{-1}$.

The recent discovery of B-modes in the polarization pattern of the Cosmic Microwave Background by the BICEP2 experiment has important implications for neutrino physics. We revisit cosmological bounds on light sterile neutrinos and show that they are compatible with all current cosmological data provided that the mass is relatively low. Using CMB data, including BICEP-2, we find an upper bound of $m_s < 0.85$ eV ($2\sigma$ Confidence Level). This bound is strengthened to 0.48 eV when HST measurements of $H_0$ are included. However, the inclusion of SZ cluster data from the Planck mission and weak gravitational measurements from the CFHTLenS project favours a non-zero sterile neutrino mass of $0.44^{+0.11}_{-0.16}$ eV. Short baseline neutrino oscillations, on the other hand, indicate a new mass state around 1.2 eV. This mass is highly incompatible with cosmological data if the sterile neutrino is fully thermalised ($\Delta \chi^2>10$). However, if the sterile neutrino only partly thermalises it can be compatible with all current data, both cosmological and terrestrial.

We introduce a cosmological invisible decay of the sterile neutrino with the eV-scale mass indicated by short-baseline neutrino oscillation experiments in order to allow its full thermalization in the early Universe. We show that the fit of the cosmological data is practically as good as the fit obtained with a stable sterile neutrino without mass constraints, which has been recently considered by several authors for the explanation of the observed suppression of small-scale matter density fluctuations and for a solution of the tension between the Planck and BICEP2 measurements of the tensor to scalar ratio of large-scale fluctuations. Moreover, the extra relativistic degree of freedom corresponding to a fully thermalized sterile neutrino is correlated with a larger value of the Hubble constant, which is in agreement with local measurements.

We analyze the most recent cosmological data, including Planck, taking into account the possible existence of a sterile neutrino with a mass at the eV scale indicated by short-baseline neutrino oscillations data in the 3+1 framework. We show that the contribution of local measurements of the Hubble constant induces an increase of the value of the effective number of relativistic degrees of freedom above the Standard Model value, giving an indication in favor of the existence of sterile neutrinos and their contribution to dark radiation. Furthermore, the measurements of the local galaxy cluster mass distribution favor the existence of sterile neutrinos with eV-scale masses, in agreement with short-baseline neutrino oscillations data. In this case there is no tension between cosmological and short-baseline neutrino oscillations data, but the contribution of the sterile neutrino to the effective number of relativistic degrees of freedom is likely to be smaller than one. Considering the Dodelson-Widrow and thermal models for the statistical cosmological distribution of sterile neutrinos, we found that in the Dodelson-Widrow model there is a slightly better compatibility between cosmological and short-baseline neutrino oscillations data and the required suppression of the production of sterile neutrinos in the early Universe is slightly smaller.