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.

We present a novel method for determining the total matter surface density of the Galactic disk by analysing the kinematics of a dynamically cold stellar stream that passes through or close to the Galactic plane. The method relies on the fact that the vertical component of energy for such stream stars is approximately constant, such that their vertical positions and vertical velocities are interrelated via the matter density of the Galactic disk. By testing our method on mock data stellar streams, with realistic phase-space dispersions and Gaia uncertainties, we demonstrate that it is applicable to small streams out to a distance of a few kilo-parsec, and that the surface density of the disk can be determined to a precision of 6 %. This method is complementary to other mass measurements. In particular, it does not rely on any equilibrium assumption for stars in the Galactic disk, and also makes it possible to measure the surface density to good precision at large distances from the Sun. Such measurements would inform us of the matter composition of the Galactic disk and its spatial variation, place stronger constraints on dark disk sub-structure, and even diagnose possible non-equilibrium effects that bias other types of dynamical mass measurements.

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 Standard Cosmological Model predicts the existence of relic neutrinos, which are indirectly probed through the effective number of relativistic species in the early Universe. In addition, from neutrino flavour oscillations we know that at least two of the neutrino mass states have a non-zero mass. Since the expansion of the Universe has diluted the energy of relic neutrinos, those that are massive are also non relativistic today. This means that they can be trapped in strong gravitational potentials, such as the one of the Milky Way. We review the calculation of the local overdensity of relic neutrinos produced from the gravitational attraction of our Galaxy, Andromeda and the Virgo cluster, commenting on the implications for an experiment aiming at relic neutrino detection, such as the PTOLEMY project.

The local density of dark matter is an important quantity. On the one hand, its value is needed for dark matter direct detection searches. On the other hand, a precise and robust determination of the local dark matter density would help us learn about the shape of the dark matter halo of our Galaxy, which plays an important role in dark matter indirect detection searches, as well as in many studies in astrophysics and cosmology. There are different methods available to determine the local dark matter density. Among them, it is common to study either the vertical kinematics of a selected group of tracers or the rotation curve of the Milky Way. Recent estimates of the local dark matter density have used the precise observations conducted by the ESA/Gaia mission. However, in spite of the quality of the data released by Gaia's observations, different analyses of the local dark matter density produce dissimilar results. After a brief review of the most common methods to estimate the local density of dark matter, here we argue about different explanations for the discrepancies in the results of recent analyses. We finish by presenting new approaches that have been proposed in the literature and could help us improve our knowledge of this important quantity.

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.

The rotation curve of the Milky Way is commonly used to estimate the local dark matter density $\rho_{{\rm DM},\odot}$. However, the estimates are subject to the choice of the distribution of baryons needed in this type of studies. In this work we explore several Galactic mass models that differ in the distribution of baryons and dark matter, in order to determine $\rho_{{\rm DM},\odot}$. For this purpose we analyze the precise circular velocity curve measurement of the Milky Way up to $\sim 25$ kpc from the Galactic centre obtained from Gaia DR2 [1]. We find that the estimated value of $\rho_{{\rm DM},\odot}$ stays robust to reasonable changes in the spherical dark matter halo. However, we show that $\rho_{{\rm DM},\odot}$ is affected by the choice of the model for the underlying baryonic components. In particular, we find that $\rho_{{\rm DM},\odot}$ is mostly sensitive to uncertainties in the disk components of the Galaxy. We also show that, when choosing one particular baryonic model, the estimate of $\rho_{{\rm DM},\odot}$ has an uncertainty of only about $10\%$ of its best-fit value, but this uncertainty gets much bigger when we also consider the variation of the baryonic model. In particular, the rotation curve method does not allow to exclude the presence of an additional very thin component, that can increase $\rho_{{\rm DM},\odot}$ by more than a factor of 8 (the thin disk could even be made of dark matter). Therefore, we conclude that exclusively using the rotation curve of the Galaxy is not enough to provide a robust estimate of $\rho_{{\rm DM},\odot}$. For all the models that we study without the presence of an additional thin component, our resulting estimates of the local dark matter density take values in the range $\rho_{{\rm DM},\odot} \simeq \text{0.3--0.4}\,\mathrm{GeV/cm^3}$, consistent with many of the estimates in the literature.

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

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

We describe the main features of a new and updated version of the program PArthENoPE, which computes the abundances of light elements produced during Big Bang Nucleosynthesis. As the previous first release in 2008, the new one, PArthENoPE 2.0, will be soon publicly available and distributed from the code site, http://parthenope.na.infn.it. Apart from minor changes, which will be also detailed, the main improvements are as follows. The powerful, but not freely accessible, NAG routines have been substituted by ODEPACK libraries, without any significant loss in precision. Moreover, we have developed a Graphical User Interface (GUI) which allows a friendly use of the code and a simpler implementation of running for grids of input parameters. Finally, we report the results of PArthENoPE 2.0 for a minimal BBN scenario with free radiation energy density.

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.

Neutrino oscillations are a widely observed and well established phenomenon. It is also well known that deviations with respect to flavor conversion probabilities in vacuum arise due to neutrino interactions with matter. In this work, we analyze the impact of new interactions between neutrinos and the dark matter present in the Milky Way on the neutrino oscillation pattern. The dark matter-neutrino interaction is modeled by using an effective coupling proportional to the Fermi constant $G_F$ with no further restrictions on its flavor structure. For the galactic dark matter profile we consider an homogeneous distribution as well as several density profiles, estimating in all cases the size of the interaction required to get an observable effect at different neutrino energies. Our discussion is mainly focused in the PeV neutrino energy range, to be explored in observatories like IceCube and KM3NeT. The obtained results may be interpreted in terms of a light $\mathcal{O}$(sub-eV--keV) or WIMP-like dark matter particle or as a new interaction with a mediator of $\mathcal{O}$(sub-eV--keV) mass.

We study the decoupling process of neutrinos in the early universe in the presence of three-flavour oscillations. The evolution of the neutrino spectra is found by solving the corresponding momentum-dependent kinetic equations for the neutrino density matrix, including for the first time the proper collision integrals for both diagonal and off-diagonal elements. This improved calculation modifies the evolution of the off-diagonal elements of the neutrino density matrix and changes the deviation from equilibrium of the frozen neutrino spectra. However, it does not vary the contribution of neutrinos to the cosmological energy density in the form of radiation, usually expressed in terms of the effective number of neutrinos, N_eff. We find a value of N_eff=3.045, in agreement with previous theoretical calculations and consistent with the latest analysis of Planck data. This result does not depend on the ordering of neutrino masses. We also consider the effect of non-standard neutrino-electron interactions (NSI), predicted in many theoretical models where neutrinos acquire mass. For two sets of NSI parameters allowed by present data, we find that N_eff can be reduced down to 3.040 or enhanced up to 3.059.

We consider the case of very low reheating scenarios ($T_{\rm RH}\sim\mathcal{O}({\rm MeV})$) with a better calculation of the production of the relic neutrino background (with three-flavor oscillations). At 95% confidence level, a lower bound on the reheating temperature $T_{\rm RH}>4.1$ MeV is obtained from Big Bang Nucleosynthesis, while $T_{\rm RH}>4.3$ MeV from Planck data for very light ($\sum m_i = 0.06$ eV) neutrinos. If neutrino masses are allowed to vary, Planck data yield $T_{\rm RH}>4.7$ MeV, the most stringent bound on the reheating temperature to date. Neutrino masses as large as 1 eV are possible for very low reheating temperatures.