1. Introduction
The successful hot big bang model after inflation predicts that neutrinos produced in the early universe still exist in the current universe. After the temperature of the universe dropped below
, weak interactions became ineffective and neutrinos would have decoupled from thermal plasma. Analogous to photons that make up the Cosmic Microwave Background (CMB), these decoupled neutrinos are called the Cosmic Neutrino Background (C
B). The existence of these relic neutrinos is confirmed indirectly by the observations of primordial abundances of light elements from the Big Bang Nucleosynthesis (BBN), the anisotropies of the CMB and the distribution of Large Scale Structure (LSS) of the universe. In particular, observations from the Planck satellite impose the severe constraint on the effective number of relativistic species
, which describes the total neutrino energy in the Standard Model (SM), and the sum of the neutrino masses at
CL as [
1]
where
and
are the energy densities of photons and radiation, which is composed of photons and neutrinos in the SM, respectively.
Future observations of the C
B will be developed both indirectly and directly. In fact, CMB-S4 observations are expected to determine
with a very good precision of ∼0.03 at 68% C.L. [
2]. Thus, an estimation of
in the SM with
precision will be important towards the future CMB-S4 observation. In addition, although it is still very difficult to observe the C
B in a direct way at present, it is inconceivable that the C
B will never be directly observed. Among the various discussions on the direct observations, the most promising method of direct detection of the C
B is neutrino capture on
-decaying nuclei [
3,
4],
, where there is no threshold energy for relic cosmic neutrinos. In both cases, the theoretical prediction of the relic neutrino spectrum is a crucial ingredient since the radiation energy density in
and the direct detection rates depend on the spectrum, and their deviations from the SM suggest physics beyond the SM.
Soon after the decoupling of neutrinos, -pairs start to annihilate and heat photons when the temperature of the universe is . If neutrinos decoupled instantaneously and all electrons and positrons annihilated into photons, the ratio for the temperatures of cosmic photons and neutrinos would be , due to entropy conservation of the universe. However, the temperatures of neutrino decoupling and -annihilations are so close that -pairs slightly annihilate into neutrinos, which leads to non-thermal distortions in neutrino spectra and a less increase in the photon temperature. These modifications are also parametrized by an increase in from 3.
The non-thermal distortions of relic neutrino spectra and the precise value of
have long been studied by solving kinetic equations for neutrinos, which are the Boltzmann equations and the continuity equation. First, several studies solved the Boltzmann equations for neutrino distribution functions [
5,
6,
7,
8,
9,
10,
11,
12]. Then the kinetic equations were solved with including finite temperature radiative corrections at leading order
[
13,
14,
15,
16,
17,
18], and then including three-flavor neutrino oscillations the Boltzmann equations for a neutrino density matrix formalism were solved [
19,
20,
21]. A fast and precise method to calculate effective neutrino temperature for all neutrino species and
was also proposed [
22,
23]. Recently, the authors in ref. [
24] pointed out that the finite temperature corrections to electromagnetic plasma at the next-to-leading order
are expected to decrease
by
. After that, the present authors found a precise value of
[
25] by solving the Boltzmann equations for the neutrino density matrix including the corrections to electron mass and electromagnetic plasma up to
but neglecting off-diagonal parts derived from self-interactions of neutrinos. Later, the authors in refs. [
26,
27] estimate
and
, respectively, including off-diagonal parts of the collision term derived by neutrino self-interactions. However, QED corrections to weak interaction rates at the order
and forward scattering of neutrinos via their self-interactions have not been precisely taken into account in the above references so far. Recent studies [
23,
28] suggest that these omissions might still induce uncertainties of
in
.
If we observe the C
B in a direct way in addition to its indirect observations, we might see neutrino decoupling directly. In the current universe, since the average momentum of the C
B is
, two massive neutrinos at least are non-relativistic. Under such a situation, it is quite nontrivial to quantize neutrinos in the flavor basis. To reveal the contribution of
-annihilation in neutrino decoupling to the spectrum of the C
B, we calculated the spectra, number densities and energy densities for relic neutrinos in the mass-diagonal basis in the current homogeneous and isotropic universe [
25,
29].
In this article, we present a review of the distorted spectra of relic cosmic neutrinos from neutrino decoupling to the current universe based on refs. [
25,
29]. First, in
Section 2, we describe the kinetic equations for cosmic neutrinos. In
Section 3, we present our results of relic neutrino spectra and
. Here we also discuss the uncertainties in
. In
Section 4, we calculate the number density and energy density of the C
B in the present universe. In
Section 5, the impact of the distortions of the spectra in neutrino decoupling on neutrino capture experiments is also discussed. One of such experiments, which is called the PTOLEMY-type experiment [
30,
31], uses 100 g of tritium [
29,
32,
33,
34,
35] as a target through the reaction,
. Tritium is an appropriate candidate for the target due to its availability, high neutrino capture cross section, low Q-value and long half lifetime of
years. Here we also include the effects of gravitational clustering of the C
B by our Galaxy and nearby galaxies based on the results in ref. [
36]. Finally, conclusions and discussion are given in
Section 6.
2. Kinetic Equations for Neutrinos in Their Decoupling
To follow relic neutrino spectra from neutrino decoupling to the current homogeneous and isotropic universe, we first discuss the field operators and the density matrix for relativistic and non-relativistic neutrinos. Then we introduce the kinetic equations for neutrinos, which are the Boltzmann equations for the evolution of the neutrino density matrix known as the quantum kinetic equations. The continuity equations for the evolution of the total energy density are also introduced.
2.1. Field Operators and Density Matrix
We consider field operators of neutrinos and their density matrices in a homogeneous and isotropic system. With neutrino masses, we cannot define annihilation and creation operators for neutrinos in flavor basis due to their off-diagonal masses in the conventional way, where we interpret these operators as operators that annihilate and create a state with eigenvalues of energy and momentum. On the other hand, in the mass-diagonal basis, we can define such annihilation and creation operators, including neutrino masses. We also compare relic cosmic neutrino spectra obtained in the two bases and confirm their match.
In the ultra-relativistic limit, the field operators for left-handed flavor neutrinos in terms of 4-component spinors, which are composed of only active states for Majorana neutrinos and both active and sterile states for Dirac neutrinos, are expanded in terms of plane wave solutions as
where
and
are annihilation operators for negative-helicity neutrinos and positive-helicity anti-neutrinos, respectively, and
H is the Hamiltonian.
and
are a flavor index and a three dimensional momentum with
, respectively.
denotes the Dirac spinor for a massless negative-helicity particle (positive-helicity anti-particle), which is normalized to be
. The annihilation and creation operators satisfy the anti-commutation relations,
For freely evolving massless neutrinos without any interactions,
and
and the Dirac spinors satisfy free Dirac equations,
. On the other hand, for free massive neutrinos in the flavor basis,
and
cannot be expanded in terms of a plane wave with an eigenvalue of their energy due to off-diagonal neutrino masses. Then we cannot interpret
and
as annihilation operators except in the ultra-relativistic case.
The density matrices for neutrinos and anti-neutrinos in the flavor basis are defined through the following expectation values of these operators concerning the initial states,
where
. Due to the reversed order of flavor indices in
, both density matrices transform in the same way under a unitary transformation of flavor space. Here the diagonal parts are the usual distribution functions of flavor neutrinos and the off-diagonal parts represent non-zero in the presence of flavor mixing.
On the other hand, in the mass-diagonal basis, the field operators for the negative helicity neutrinos
1 can be expanded as, including neutrino masses,
where
denotes a mass eigenstate,
,
and
is the neutrino mass in the mass basis.
denotes the Dirac spinor for negative-helicity particles (positive-helicity anti-particles), which is also normalized to be
. For freely evolving neutrinos,
and the Dirac spinors satisfy
and
As in the flavor basis, the commutation relations for
and
, and the density matrix are defined in the same way except for the exchange of the subscripts,
.
The diagonalization of the mass matrix for left-handed neutrinos in the flavor basis is achieved through the transformations,
where
represents a component of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix
. Due to Equation (6), in the ultra-relativistic limit, the relation of the density matrices in the flavor and the mass bases is described as
In addition, after neutrino decoupling, the off-diagonal parts of the density matrix in the mass basis are zero,
, since all neutrino interactions are ineffective and the oscillations do not occur after neutrino decoupling. In this case, the relations of distribution function in the two bases are simply
2 Note that Equation (8) is only valid when neutrinos are relativistic and decoupled with thermal plasma. Our numerical calculations also confirm Equation (8).
2.2. Boltzmann Equations
In this section, we derive the Boltzmann equations for the neutrino density matrix, known as quantum kinetic equations, including neutrino oscillations in vacuum, forward scattering with
-background, corresponding to neutrino oscillations in matter, and the collision process at tree level. The resulting Boltzmann equations for neutrinos are summarized in
Section 2.5, where we will also discuss the approximations we used in our numerical calculations.
2.2.1. Boltzmann Equations in a Homogeneous and Isotropic System
The Boltzmann equations for neutrinos, including flavor conversion effects, are derived from the Heisenberg equations for the neutrino density operator,
where
represents the commutator of matrices with a flavor (or mass) index and
is the neutrino density operator,
H is the full Hamiltonian in a system, which can be separated into
where
is the free Hamiltonian and
is the interaction Hamiltonian. We assume interactions are enough small that collisions occur individually. Then any fields can be regarded as free ones except during interactions. When the interaction Hamiltonian can be treated perturbatively, the density operator evolves at the first order of
,
where
is the initial time and
is the interaction Hamiltonian as a function of freely evolving fields, which are solutions of free Dirac equations, and
is the free density operator evolved as
The first order solution (12) includes only neutrino oscillation in vacuum and forward (momentum conserving) scattering with a medium in the system.
To take into account momentum changing collisions, we consider the evolution equation for the density operator at second order of
, substituting Equation (12) into Equation (9),
and an analogous equation for anti-neutrinos [
37],
, which is not solved in this article since we assume no lepton asymmetry. Here
is also the free Hamiltonian as a function of freely evolving fields, where we neglect would-be tiny corrections in the presence of interactions. We also ignore the tiny modification of oscillation and forward scattering,
compared with
and
. Note that the differential Equation (14) is not closed for both
and
.
To close and simplify the differential Equation (14), we impose additional approximations. We may set
and
in the integral range since the time step of the change of
,
t, may be chosen to be small enough compared to the timescale of the evolution of the universe and large enough compared to the timescale of one collision,
. In addition, at
, the free density operator coincides with the full one,
. Then Equation (14) can be rewritten as
Thus, the time evolution of the expectation value of
concerning the initial state,
, is given by
Equation (16) will be valid at all times, even at , if in two or more collisions, the correlation of the particles in each collision is independent. This assumption is called molecular chaos in the derivation of the Boltzmann equation. In general, n-point correlation functions are produced by both forward and non-forward collisions. Under the assumption of molecular chaos, n-point correlation functions are reduced to combinations of two-point correlation functions as in ordinary scattering theory. Here two-point correlation functions correspond to distribution functions and neutrino density matrix.
The first term in the right hand side (RHS) represents neutrino oscillations in vacuum and the second term represents forward scattering of neutrinos with background in the system, which is called refractive effects and corresponds to neutrino oscillations in matter. These two terms do not change neutrino momenta but induce flavor conversions. The third term represents scattering and annihilation including both momentum conserving and changing processes, usually rewritten as
where
is called the collision term. In the following sections, we calculate the formulae of these three terms. The resulting Boltzmann equations for the neutrino density matrix are summarized in
Section 2.5.
2.2.2. Neutrino Oscillation in Vacuum
The calculation of the first term in the RHS of Equation (16) is well established in the mass basis. The free Hamiltonian of neutrinos in the mass basis is given by
where
are the gamma matrices. After substituting the free operators for left-handed neutrinos, the free Hamiltonian becomes
The first term in the RHS of Equation (16) in the mass basis is written as
where
and
denote mass-eigenstates. In the flavor basis, as in discussed in
Section 2.1, it is quite nontrivial to quantize neutrinos in the flavor basis with non-zero masses. When we calculate the first term in the RHS of Equation (16) in the flavor basis directly, we replace the free annihilation operators
and
with
and
as in [
37], where
. Then we also obtain the first term of Equation (16), following the similar procedure in the mass basis,
where
is the neutrino mass matrix in the flavor basis. For anti-neutrinos, the corresponding term is obtained by adding a minus sign for the reverse indices in the anti-neutrino density matrix (4),
.
2.2.3. Forward Scattering with -Background
In the following of
Section 2, we consider the flavor basis of neutrinos. Forward scattering of neutrinos with background in the system called refractive effects modifies neutrino oscillations through the one-loop thermal interaction as given in
Figure 1. Since the temperature in thermal plasma is ∼2
in neutrino decoupling, particles except for photons, electrons, neutrinos and their anti-particles are already annihilated due to their heavy masses. Then we consider only
-background. The interaction Hamiltonian is described as
where
and
are the full propagator of
boson and
boson,
Here
are the electroweak coupling constant, the
boson mass and the
boson mass, respectively. The neutral current and the charged current are given by
where
with
Here is the weak mixing angle, is the field operator for electron and positron and is the field operator for neutrinos and anti-neutrinos with a flavor .
The interaction Hamiltonian is divided into the two parts corresponding to the neutral current interaction,
, and to the charged current interaction,
. For the charged current interactions, the second term in the RHS of Equation (16), which represents forward scattering of neutrinos with
-background, is given by [
37,
38]
where
is the Fermi coupling constant and
and
are the number density, energy density and pressure for
-background, respectively, which are described in the flavor basis as
where
. In the temperature of MeV scale in neutrino decoupling, the densities for muons and tauons are enough suppressed by their heavy masses. We neglect forward scattering of neutrinos with muons and tauons, which corresponds to the second and third diagonal components in Equation (28).
For the neutral current interactions, the second term in the RHS of Equation (16), which represents forward scattering of neutrinos via neutrino self-interactions, is given by [
37,
38]
where
and
are the number and energy densities for the density matrices of
-background, respectively, which are described in the flavor basis as
where we neglect neutrino masses since neutrinos are relativistic in neutrino decoupling.
For anti-neutrinos, the corresponding terms,
and
, are obtained by adding an overall minus sign for the reverse indices in the anti-neutrino density matrix (4) and replacing
and
for an opposite evolution of anti-neutrinos due to the lepton asymmetry in Equations (27) and (29) [
37].
If there is a large lepton asymmetry, the terms proportional to and will be important. Note that even if there is no lepton asymmetry, the off-diagonal parts of have non-zero contribution since the density matrices for neutrinos and anti-neutrinos follow the same evolution, , in the case of no lepton asymmetry.
2.2.4. Collision Term
Finally we discuss the third term in the RHS of Equation (16) called the collision term. The temperature of ∼2
in neutrino decoupling is much lower than the electroweak scale of ∼
. After integrating out
and
bosons in the instantaneous interaction limit, the interaction Hamiltonian in neutrino decoupling can be written as
The interaction Hamiltonian can be divided into the part including both neutrinos and electrons (and their anti-particles), and the one only including neutrinos and anti-neutrinos,
, while we ignore the part including only electrons and positrons,
with
Here we have used the following Fierz transformation in the charged currents,
Due to this contribution of the charged current, only electron-type neutrinos and anti-neutrinos interact with electrons and positrons via different magnitudes of interactions, compared to other flavor neutrinos with
.
The Hamiltonian of Equation (32) can be further divided as
where
is the term including operators of (anti-)particles,
and
d. In the following, we neglect
since this Hamiltonian does not contribute the evolution of neutrinos. In addition, we only consider contributions proportional to the following terms as a function of freely evolving fields in the collision term in Equation (16),
The other terms also denote forward scattering, which would give tiny modifications of Equations (27) and (29). The first term in Equation (36) denotes the annihilation of neutrinos and anti-neutrinos into -pairs, which mainly contribute to the distortion of neutrino spectrum in their decoupling. The second term denotes the scattering between neutrinos and electrons (positrons). The third term represents the scattering process including only neutrinos while the fourth term denotes the annihilation and scattering processes of neutrinos and anti-neutrinos.
In a schematic manner, the collision term for two-body reactions
at tree level takes the following expressions,
where
denote the neutrino density matrix, not the energy density, and
for
and
while
for
.
is a matrix depending on
,
,
and/or
.
is a part of
, where
S is the symmetric factor and
is the squared matrix element summed over spins of all particles except for the first one. The formulae of
for the relevant reaction in neutrino decoupling are shown in
Table 1. Nine integrals in the collision term in Equation (37) can be reduced analytically to two integrals as in
Appendix B.
In the following, we rewrite the collision terms
including Equation (36) with neutrino density matrices and the distribution functions of electrons and positrons. The formulae of the collision terms for neutrino density matrix are originally given in refs. [
37,
39], and for numerical calculations of neutrino spectra, these formulae are developed in refs. [
20,
26].
- (i)
The collision term for the annihilation process including
,
, comes from the term proportional to
. We can calculate the corresponding collision terms, which are denoted as
,
where
Here is the distribution function for electrons and positrons, respectively.
- (ii)
The collision term for the scatterings including
,
, comes from the term proportional to
. We can similarly calculate the corresponding collision term, which is denoted as
and
, respectively,
and
where
- (iii)
and
The collision terms for the scatterings including only neutrinos and anti-neutrinos,
and
, come from the term proportional to
and
, respectively. The corresponding collision terms, which are denoted as
and
, respectively, are calculated as
where
,
and
denote contributions from scatterings for
, scatterings and annihilations for
, respectively,
where
.
Finally, we obtain the collision term in Equation (16),
, combining Equations (38), (40), (41), (43) and (44),
The collision terms for anti-neutrinos can be obtained by appropriately replacing the density matrices and momenta,
and
[
37,
40]. Changing the collision term for
to
corresponds replacing
and
in this collision term while changing that for
to
corresponds
and
. One may consider the transpose in the collision terms is necessary for the reverse indices in the anti-neutrino density matrix (4), but this is not necessary since the collision terms are invariant under the transpose.
2.3. Continuity Equation
In addition to the Boltzmann equations for the neutrino density matrix, the energy conservation law must be satisfied,
where
and
P are the total energy density and pressure of
around MeV-scale temperature, respectively. The continuity equation corresponds to the evolution of the photon temperature
.
Though we will discuss finite temperature corrections from QED to
and
in the next section, in the ideal gas limit, they are given as follows, which are denoted by
and
, respectively,
The Hubble parameter in Equation (49) is calculated using the usual relation, with being the Planck mass, where we ignore the curvature term and the cosmological constant because they are negligible in the radiation dominated epoch.
2.4. Finite Temperature QED Corrections to , and P up to
QED interactions at finite temperature modify the energy density and pressure of electromagnetic plasma from the ideal gas limit. In addition, their interactions change the electron mass (and produce an effective photon mass). These corrections affect the kinetic equations for neutrinos discussed in the former sections. The corrections to the electron mass modify the weak interaction rates and the distribution function for
. Through the direct modifications of
and
P, the expansion rate
H is also changed. Note that QED interactions also modify weak interaction rates in the collision term
and the Hamiltonian for the forward scattering (60) at order
directly. In our numerical calculations, we consider corrections to weak interaction rates only due to the change of
. We will discuss other QED corrections to weak interaction rates and their uncertainties in
in
Section 3.3.1.
The corrections to the grand canonical partition function
Z by interactions at finite temperature are well established perturbatively and can be calculated by the similar procedure of the functional integrals of Quantum Field Theory (QFT) at zero temperature after changing
.
P and
are described by
Z as
where
T and
V are the temperature and volume in the system, respectively. Then we can expand
in powers of the QED coupling constant
e as
, where
. In the isotropic and lepton symmetric universe, the corresponding corrections to
P and
at
,
, are [
41]
where
and
is the sum of the distribution functions for
,
The next-to-leading order of thermal corrections to
is
, not
. These non-trivial corrections come from the resummation of ring diagrams in the photon propagator at all orders. The thermal corrections to
at
,
, are [
24,
41],
where
Finally, we read the total energy density and the total pressure of electromagnetic plasma up to
corrections as
The thermal corrections to the
mass at
is given by, through modifications of the
self energy [
42],
The last logarithmic terms in Equations (52) and (57) give less than
corrections to these equations around the decoupling temperature and the average momentum of electrons [
43]. These terms also give contributions less than
to
[
24,
27]. In the following, we neglect the logarithmic corrections. Note that thermal corrections to
at
do not appear because
corrections stem from ring diagrams in the photon propagator.
2.5. Summary and Approximations
In this section, we summarize the closed system of the resulting Boltzmann equations for the neutrino density matrix and the continuity equation in neutrino decoupling. We also discuss the approximations we used in our numerical calculations. The following Equations (58)–(61) have already been presented in the previous sections.
The closed system of the equations of motion for the neutrino density matrix and the continuity equation, which reads the equation of the evolution for the photon temperature, in the expanding universe are [
37,
39]
and analogous Boltzmann equations for anti-neutrinos [
37,
40], which is not solved in this article since we assume no lepton asymmetry. Here
is the Hubble parameter,
is the Hamiltonian which governs the neutrino oscillation in vacuum and the forward scattering of neutrinos in the
-background,
is the collision term describing the momentum changing scatterings and annihilations, and
represents the commutator of matrices with a flavor (or mass) index.
and
P in Equation (59) are the total energy density and the pressure for
, respectively. Including QED finite temperature corrections up to
,
and
P are given by Equation (56) (see also Equations (50), (52) and (54) for the detail of Equation (56)).
The effective Hamiltonian for the neutrino oscillations in vacuum and the forward scattering of neutrinos in the
-background is given by
3 where
GF is the Fermi coupling constant and
mW,
mZ are the W and Z boson masses, respectively.
The first term in the RHS of Equation (60) denotes neutrino oscillations in vacuum and
is the mass-squared matrix. In the flavor basis, we can write
, where
. The other terms describe the forward scattering of neutrinos in the background of thermal plasma which comes from one-loop thermal contributions to neutrino self energy.
are defined in the flavor basis around the temperature of MeV scale as
where
and
is the distribution function of
.
is the QED finite temperature correction to
, which is given by Equation (57) up to
. Here we neglect the contributions of
and
since the densities of these charged particles are significantly suppressed.
In the following, we assume that electrons and positrons are always in thermal equilibrium and follow the Fermi–Dirac distributions since electrons, positrons and photons interact with each other through rapid electromagnetic interactions. In addition we neglect lepton asymmetry since neutrino oscillations leading to flavor equilibrium before the BBN imposes a stringent constraint on this asymmetry [
44,
45,
46,
47,
48,
49,
50]. The standard baryogenesis scenarios via the sphaleron process in leptogenesis models predict that the lepton asymmetry is of the order of the current baryon asymmetry,
, which is much smaller than the above constraint. We also neglect any CP-violating phase in the PMNS matrix for simplicity. Note that from the recent global analysis of neutrino oscillation experiments [
51,
52], the CP-conserving PMNS matrix is excluded at approximately
confidence level. Strictly speaking, ignoring the CP-violating phase is inconsistent with the experimental results, but we adopt this assumption to save computational time. In fact, since effects of CP-violating phase on neutrino oscillations are sub-dominant, this ignorance will not affect the resultant neutrino spectra and
significantly. Under these assumptions, neutrinos and anti-neutrinos satisfy the same density matrices and the same evolutions in the Universe,
, and electrons and positrons follow the same Fermi–Dirac distributions with
and no chemical potential.
Note that without lepton asymmetry,
due to
. However, in the following, we neglect it for reducing computational time. We will discuss this uncertainty in
Section 3.3. In addition, as in refs. [
19,
20,
21,
25], we replace
as
for simplicity. Strictly, this replacement is valid only in the ultra-relativistic limit [
38]. However, since in the non-relativistic region
is suppressed by the Boltzmann factor, these difference would be quite small. Ref. [
27] reported this difference in
is no more than
.
The final term in the RHS of Equation (58) represents both the momentum conserving and changing collisions of neutrinos with neutrinos, electrons and their anti-particles. In this term, collisions are dominated by two-body reactions
, i.e.,
, where
is the Fermi coupling constant. The detailed formula for
is given by Equation (48) (see also Equations (38), (40), (41), (43) and (44) in this review and refs. [
20,
26]). Nine integrals in the collision term (37) can be reduced analytically to two integrals as in
Appendix B. We deal with both diagonal and off-diagonal collision terms in Equations (38), (40) and (41) for the processes which involve electrons and positrons,
and
. On the other hand, we do not treat the off-diagonal terms in Equations (43) and (44) for the self-interactions of neutrinos,
and
, since the annihilations of electrons and positrons are important for the heating process of neutrinos while the self-interactions of neutrinos less contribute to this heating process. We treat this collision term from neutrino self-interaction in Equation (A13) of
Appendix A. In refs. [
26,
27], the authors solve kinetic equations for neutrinos including the full collision term at tree level and reported almost the same results with very small difference in
,
[
27]. Here, we take into account finite temperature corrections to
up to
in the collision term as
. However, we neglect other sub-leading contributions to the collision term, i.e., other QED corrections to weak interaction rates. We also discuss these uncertainties in
Section 3.3.
2.6. Computational Method, Initial Conditions and Values of Neutrino Masses and Mixing
We solve kinetic equations for neutrinos of Equations (58) and (59) with the following comoving variables instead of the cosmic time
t, the momentum
p, and the photon temperature
,
where we choose an arbitrary mass scale in
x to be the electron mass
and
a is the scale factor of the universe, normalized as
in high temperature limit. The resultant kinetic equations for neutrinos in the comoving variables are described in
Appendix A.
Since the Boltzmann Equations (58) are integro-differential equations due to integrations in the collision terms, their equations were solved by a discretization in a momentum grid
in refs. [
8,
9,
10,
16,
19,
20,
21], by an expansion of the distortions of neutrinos from the Fermi–Dirac distribution in refs. [
11,
14,
15], or by a hybrid method combining the previous two methods in ref. [
18]. In this study, we adopt the discretization method we mentioned first and take 100 grid points for
, equally spaced in the region
with the Simpson method. We have used MATLAB ODE solver, in particular, ode15s with an absolute and relative tolerance of
. In these tolerances, we confirm that numerical errors for relic neutrino spectra and
are typically
or less.
We have numerically estimated the evolution of the density matrix for neutrinos and the photon temperature in
. We have set
as an initial time. Since neutrinos are kept in thermal equilibrium with the electromagnetic plasma at
, the initial values of density matrix
are regarded as
The initial dimensionless photon temperature at
,
, slightly deviates from 1 because a tiny amount of
-pairs have already been annihilated at
. Due to the entropy conservation of electromagnetic plasma, neutrinos and anti-neutrinos,
is estimated as in [
10],
We take
as a final time, when the neutrino density matrix and
z can be regarded as frozen.
Finally we comment on values of neutrino masses and mixing we use in our numerical simulation. We use the best-fit values in the global analysis in 2019 [
53], but assume CP-symmetry,
. We note that in 2020 their best-fit values are updated [
51,
52] though their differences are very small. Their parameters include small uncertainties of about
at
confidence level. Effects of their uncertainties on
is investigated in ref. [
27] and slightly change
by
. In our numerical simulation, we confirmed that relic neutrino spectra and the value of
with
precision are the same for both neutrino mass ordering. In the following, we show the results in the normal mass ordering,
, not in the inverted ordering,
, because the results do not change significantly.
3. Effective Number of Neutrino Species
To describe the process of neutrino decoupling, we first numerically solve a set of Equations (58) and (59) and show relic neutrino spectra in the flavor basis. Then we present a precise value of the effective number of neutrino species, , and discuss effects of neutrino oscillations and finite temperature corrections to and P up to on . We also comment on uncertainties of ingredients we ignored in estimating .
3.1. Relic Neutrino Spectra in the Flavor Basis
In the left panel of
Figure 2, we show the distortions of the flavor neutrino spectra for a comoving momentum
, where we plot the neutrino spectra
as a function of the normalized cosmic scale factor
x.
is the neutrino distribution function if neutrinos decoupled instantaneously and all
-pairs annihilated into photons,
At high temperature with , the temperature differences between photons and neutrinos are negligible and neutrinos are in thermal equilibrium with electrons and positrons. In the intermediate regime with , weak interactions gradually become ineffective with shifting from small to large momenta. In this period, the neutrino spectra are distorted since the energies of electrons and positrons partially convert into those of neutrinos coupled with electromagnetic plasma. Finally, at low temperature with , the collision term becomes ineffective and the distortions are frozen.
The difference between the spectrum and the spectrum without flavor mixing arises from the fact that only electron-type neutrinos interact with electrons and positrons through the weak charged currents. On the other hand, in the cases with neutrino mixing, neutrino oscillations mix the distortions of the flavor neutrinos too.
In the right panel of
Figure 2, we show the frozen values of the flavor neutrino spectra
as a function of a comoving momentum
y for both cases with and without neutrino mixing. This figure shows the fact that neutrinos with higher energies interact with electrons and positrons until a later epoch. In addition, we see neutrino oscillations tend to equilibrate the flavor neutrino distortions. Although the neutrino spectra
with low energies are very slightly less than unity, these extractions of low energy neutrinos stem from an energy boost through the scattering by electrons, positrons, (and neutrinos) with sufficiently high energies, which are not yet annihilated and hence still effective at neutrino decoupling process.
3.2. Value of the Effective Number of Neutrino Species
The effective number of neutrinos
can be rewritten,
where
and
. In
Table 2 and
Table 3, we present final values (at
) of the dimensionless photon temperature
, the difference of energy densities and number densities of flavor neutrinos from those where neutrinos decoupled instantaneously denoted by
and
, and the effective number of neutrinos
.
By comparing values of
in the cases without QED corrections and with QED corrections to
and
P up to
and
in
Table 2, we find that the QED corrections at
and
shift
by
and
, respectively, which is very close to the value estimated in the instantaneous decoupling limit [
24].
In the cases with neutrino mixing,
Table 3 shows that the energy densities of
-type neutrinos increase more while those of electron-type neutrinos increase less, compared to the cases without neutrino mixing. This modification leads to the enhancement of the total energy density for neutrinos with final values of
with QED corrections to
and
P up to
. Since the blocking factor for electron neutrinos,
is decreased by neutrino mixing, the annihilation of electrons and positrons into electron neutrinos increases. Although the annihilation into the other neutrinos decreases, electron neutrinos contribute to the neutrino heating most efficiently, and neutrino oscillations enhance the annihilation of electrons and positrons into neutrinos. From these processes, we conclude that neutrino oscillations slightly promote neutrino heating and the difference of
is
, which agrees with the results of previous works [
12,
20,
23].
To conclude, our numerical calculation with neutrino oscillations and QED finite temperature corrections to
and
P up to
finds
. This value is in excellent agreement with later independent works [
26,
27].
3.3. Discussions of Uncertainties in
We comment on possible errors of the results for relic neutrino spectra and due to approximations in Equations (58) and (59) and the choice of physical parameters. Our numerical calculations converge very well since we have directly computed in the mass basis as will be done in the next section and obtained .
First we neglect the off-diagonal parts for neutrino self-interactions in the collision term,
and
. Later, in refs. [
26,
27], the authors solve kinetic equations for neutrinos including their off-diagonal parts in the collision term and report the difference in
is
[
27]. We also neglect the
logarithmic terms and terms above
in QED finite temperature corrections to
,
and
P. Their corrections to
and
P are reported to contribute
to
in refs. [
24,
27]. Though their corrections to
are not taken into account, the corrections to
even at
contribute
to
[
27] and we have also confirmed it.
The neutrino masses and mixing parameters contain 10–20% uncertainties at
confidence level. Since in our estimations, neutrino oscillations contribute
to
, their uncertainties are expected to be quite small. In ref. [
27], the authors report that their uncertainties are
. We also neglect the CP-violating phase
in the PMNS matrix. No one has yet computed precise neutrino evolution in the decoupling including three-flavor oscillations with CP violating phase. However, since effect of the CP-violating phase on neutrino oscillations is sub-dominant, we expect neutrino and anti-neutrino spectra might not change significantly. In addition, the total energy density, i.e.,
would change much less than
since the changes for the energy densities of neutrinos and anti-neutrinos would be canceled out. See also discussion in Appendix F of ref. [
26] and ref. [
54]. Other physical parameters for electroweak interaction are measured very precisely and will not affect neutrino spectra and
.
However, QED corrections to weak interaction rates at order and forward scattering of neutrinos via their self-interactions have not been precisely taken into account in the whole literature so far.
3.3.1. QED Corrections to Weak Interaction Rates at Order
QED interactions also modify the weak interaction rates in the collision term
and the Hamiltonian for the forward scattering of neutrinos (60) at order
in addition to the modification of the energy density and pressure for electromagnetic plasma,
and
P. These corrections are partially taken into account by considering thermal QED corrections on
so far. See also Section 3.1.2 in ref. [
27].
QED corrections to the weak interaction rates (see also the diagrams in
Figure 3) are categorized as (i) additional photon emission and absorption, (ii) corrections to the dispersion relation for external
, (iii) vertex corrections, and (iv) corrections mediated by photon propagator. The interference among the weak interaction at leading order
and corrections (i)–(iv) produce modifications to the weak interaction rates at the next-to-leading order
.
The correction (i) might be the most dominant contribution to since the photon emission processes, e.g., , would not be suppressed by the distribution function of photons in the Boltzmann equations. The photon emission processes reduce . However, there are many processes in the categories (ii), (iii) and (iv). In total, these contributions to might be as large as that from the correction (i).
For category (ii), corrections to the dispersion relation for
produce a thermal electron mass as Equation (57). One can incorporate corrections (i) in the weak interaction rates by shifting
, but it is numerically difficult to take into account the momentum-dependent part of
, which corresponds to the logarithmic
corrections to
. These logarithmic
corrections to
are less than
of corrections at
to
around neutrino decoupling [
43], and corrections even at leading
to the weak interaction rates (i.e.,
) contributes
to
[
27] and we confirmed it. Thus, we would properly be able to incorporate corrections (i) to
with
precision. However, we should carefully derive these corrections to the weak interaction rates and consider effects of the logarithmic
corrections and other sub-dominant neglected contributions in the collision term in the future.
For categories (i), (iii) and (iv), corrections to the weak interaction rates are typically momentum-dependent. It would be quite difficult to solve the Boltzmann equation, which is the integro-differential equation, including such momentum-dependent corrections. In ref. [
55], the authors consider energy loss rate of a stellar plasma, including corrections on
at order
and found such corrections modify the energy loss rate of a stellar plasma by a few percent. In ref. [
23], the author suggests
due to correction (i) by roughly extrapolating the results in ref. [
55] and using a precise and simple evaluation method of
proposed in ref. [
23]. The contributions of (i), (iii) and (iv) to
should be evaluated in the future in a more precise way.
3.3.2. Forward Scattering of Neutrinos via Their Self-Interactions
In the Hamiltonian (60) in the Boltzmann Equations (58), the forward scattering terms of neutrinos via their self-interactions correspond to
Even in the case without lepton asymmetry,
due to
in general, where
is
Though
might be small, forward scattering via neutrino self-interactions could be more dominant than neutrino oscillation in vacuum, with a typical dimensional analysis,
In ref. [
28], the authors suggest forward scattering of neutrinos via their self-interactions contributes
to
by solving a simplified kinetic equations for neutrinos. In the future, relic neutrino spectra and
should be estimated including the above forward scattering of neutrinos more precisely.
Though recent estimations might contain uncertainties of in , would still be one of very good reference values in .
4. Relic Cosmic Neutrino Spectra in the Current Homogeneous and Isotropic Universe
In the current universe, two neutrino species at least are non-relativistic. Then relic neutrino spectra in the mass basis will be important observable to detect the C
B in a direct way as discussed in
Section 2.1. In this section, we present the spectrum (as a function of comoving momenta), number density and energy density of the C
B in the current homogeneous and isotropic universe, including non-thermal distortions due to
-annihilation during neutrino decoupling.
4.1. Relic Neutrino Spectra in the Mass Basis
We present relic neutrino spectra in the mass basis by solving a set of Equations (58) and (59) in the mass basis directly. We can also obtain the same result by transforming relic neutrino spectra in the flavor basis through Equation (8).
In the mass basis, the neutral and charged currents including left-handed neutrino fields in Equation (24) are given by, using
as in Equation (6),
Then, using the relations of Equation (70) and (34), we obtain the 4-point interaction Hamiltonian (32) in the mass basis
with
Then we obtain the Boltzmann equation for the neutrino density matrix in the mass basis after replacements of and analogous to in Equation (58) for the flavor basis.
In the left panel of
Figure 4, we show the evolution of the neutrino spectra,
, for a comoving momentum
as a function of the normalized scale factor
x. In the right panel of
Figure 4, we show the asymptotic values of the neutrino spectra
4 as a function of
y. The differences of distortions for each neutrino species arise from the charged current interactions between neutrinos and electrons weighted by the PMNS matrix with mass species
i,
, as in Equation (71). Note that neutral currents between neutrinos in the mass basis are the same as that in the flavor basis except for the subscript,
. Then the scattering and annihilation among neutrinos and electrons and their anti-particles induce the spectral distortions in
Figure 4.
Finally we comment on . After we directly solve a set of Equations (58) and (59) in the mass basis, including vacuum three-flavor neutrino oscillations, forward scatterings in -background, and QED corrections to , and P up to , we find , which is an excellent agreement with our calculation in the flavor basis. The tiny difference from in the flavor basis may come from ignoring the off-diagonal parts for self-interaction processes in the Boltzmann equations and/or numerical errors.
4.2. Neutrino Number Density and Energy Density in the Current Homogeneous and Isotropic Universe
In
Table 4, we show the final values of the dimensionless photon temperature
, the relativistic energy densities
and number densities
of neutrinos in the mass basis after neutrino decoupling. Note that the expression of energy density for a relativistic particle is not applicable to the first and second heaviest neutrinos today because they are non-relativistic in the current universe.
After neutrino decoupling, the neutrino momentum distribution in the homogeneous and isotropic universe can be parametrized as
is the effective neutrino temperature, which is
and normalized as
in high temperature limit. Under this definition of
, neutrino spectral distortions,
, can be rewritten as
given in the right panel of
Figure 4. At
in the current universe,
satisfies
where
is the effective photon temperature in the current universe [
56]. Then the effective neutrino temperature in the current universe is
Neutrino number density and energy density per one degree of freedom in the current universe are also parametrized as
where
and
are given by
Then
and
are given in
Table 4. The values of neutrino number density in the current universe are listed in
Table 5.
In the current universe, two species of cosmic relic neutrinos at least are non-relativistic because of
. On the other hand, the lightest neutrinos might be relativistic in the current universe because the lightest neutrino mass is not yet determined. In
Table 6, we show energy density for the lightest neutrinos in the case of
. Here we consider both the normal mass ordering,
, and the inverted mass ordering,
.
To estimate the effects of
-annihilation into neutrinos during neutrino decoupling on neutrino number density and energy density, it is useful to compare the neutrino number density and relativistic energy density per one degree of freedom in the case when all
-pairs annihilate into photons,
and
, respectively,
where
. We show the deviation of neutrino number density from the case when all
-pairs annihilate into photons,
, in
Table 7. The number densities for all neutrino species are enhanced by about
due to
-annihilations to neutrinos during neutrino decoupling and the number density for
is most efficiently enhanced.
4.3. Helicity of Relic Majorana Neutrinos vs. Dirac Neutrinos
The weak interaction is chiral, which is manifest in the Lagrangian. Due to its chirality, the left-chiral states for SM fermions interact with the weak bosons while the right-chiral states do not. In the early universe, only left-chiral neutrinos and right-chiral anti-neutrinos, i.e., left-handed neutrinos and right-handed anti-neutrinos are produced via the weak interaction. Note that chirality is different from helicity in general, which is defined as the projection of the spin vector onto the momentum vector.
During free streaming of relic neutrinos after their decoupling, the chirality for non-relativistic neutrinos is not conserved since the chiral symmetry in the free neutrino Lagrangian is broken due to their masses. On the other hand, the helicity for relic neutrinos is conserved in the homogeneous and isotropic universe. Thus, we should estimate the spectrum for each helicity state of relic cosmic neutrinos in the current universe.
In the early universe, both chirality and helicity for relic neutrinos are conserved and then neutrino helicity and chirality have one-to-one correspondence since neutrinos are approximately massless in the early universe. We define left (right) helical neutrinos with helicity
such that they correspond to left (right) handed neutrinos in the early universe. Then the spectra for the left-handed neutrinos (right-handed anti-neutrinos) produced in the early universe are translated into the left-helical neutrinos (right-helical anti-neutrinos) [
34],
where
is given by Equation (73) and
if we neglect lepton asymmetry. Here right-helical neutrinos,
with
, (left-helical anti-neutrinos,
with
,) corresponds to right-handed neutrinos (left-handed anti-neutrinos), which are sterile states. We assume sterile neutrinos are not produced in the early universe due to very weak interactions with the SM particles or have already decayed if sterile neutrinos are right-handed heavy Majorana particles as required for the see-saw mechanism.
For Majorana neutrinos, right-handed active anti-neutrinos are regarded as right-handed active neutrinos due to the lepton number violation. Then
for
are given by
where
denotes a sterile state of neutrino. Note that even in the case of Majorana neutrinos lepton asymmetry can be interpreted as chiral asymmetry between left-handed and right-handed neutrinos. Then
and
are different strictly speaking but almost the same approximately.
For Dirac neutrinos, since right-handed neutrinos and left-handed anti-neutrinos are sterile,
for
are given by
where
denotes a sterile state of anti-neutrino.
From Equations (81) and (82), the magnitude of relic neutrino spectra summed over helicity for Majorana and Dirac neutrinos differ by a factor of two, which is first pointed out in ref. [
34],
Then number density and energy density summed over helicity for Majorana and Dirac neutrinos also differ by a factor of two,
5. Implications for the Capture Rates on Cosmic Neutrino Capture on Tritium
Finally we discuss how neutrino spectral distortions from
-annihilations during neutrino decoupling affect direct detection of the C
B on tritium target, with emphasis on the PTOLEMY-type experiment [
30,
31], where cosmic neutrinos can be captured on tritium by the inverse beta decay process without threshold energy for neutrinos,
. Tritium is one of appropriate candidates for the target because of its availability, high capture rate for neutrinos, low Q-value and long half lifetime of
years. Here we take 100 g of tritium as the target. We take into account gravitational clustering for cosmic neutrinos in our Galaxy and nearby galaxies because we would observe the C
B directly inside our Galaxy. We also comment on gravitational helicity flipping and annual modulation for the C
B. Then we discuss the potential of direct measurements of such cosmological effects although it would be still extremely difficult to observe such effects directly. In particular, we compute the capture rates of cosmic relic neutrinos on tritium, including such cosmological effects.
5.1. Gravitational Effects for the CB
5.1.1. Clustering for the CB by Our Galaxy and Nearby Galaxies
Near the Earth, non-relativistic relic neutrinos cluster locally in the gravitational potential of our Galaxy and nearby galaxies. Then the local distribution function is distorted and the local number density is enhanced compared with the global distribution function and number density. The local number density for relic neutrinos in the current universe is described as
where
is an enhancement factor by the gravitational attraction by galaxies, which is estimated in refs. [
36,
57,
58,
59,
60,
61]. For reference, we display some of these values, estimated in a recent numerical study [
36], in
Table 8, where the authors consider the gravitational potential in the Milky Way, Virgo cluster, and Andromeda galaxy. Note that so far, when evaluating values of
, effects of
-annihilations into
during neutrino decoupling have not been taken into account simultaneously. For
, spectral distortions to the momentum distributions for relic cosmic neutrinos by the gravitational clustering have not also been explicitly estimated (see ref. [
58] for spectral distortions by gravitational clustering for relic neutrinos with
).
In the following, we discuss only the case where
and the lightest neutrino mass is quite small because the Planck satellite suggests
. Then the local number density for relic neutrino can be parametrized as, using linear approximation,
where
is the enhancement factor by
-annihilations into
and
during neutrino decoupling given in
Table 7.
5.1.2. Helicity Flipping and Annual Modulation for the CB
We shortly comment on gravitational helicity flipping and annual modulations for relic neutrinos. Gravitational clustering for massive neutrinos may induce mixing of relic neutrino helicity [
34,
35,
62] since the direction of neutrino momentum would change in the gravitational potential for our Galaxy whereas its spin does not. Although the quantitative calculations have not yet been achieved, the capture rates on tritium would not change since their capture rates depend on neutrino number density summed over helicities at leading order as we will see in the next section. In addition, an annual modulation for relic neutrinos might occur in a direct detection experiment for the C
B since their velocity relative to the Earth could be anisotropic due to neutrino clustering and the gravitational focusing for the C
B by the Sun could also occur. The former effect is negligible since the capture rates on tritium target are independent of neutrino velocity as we will see in the next section. The latter effect is expected to change the capture rates by much less than 1% for
[
63]. In the following, we neglect helicity flipping and annual modulation for relic neutrinos.
5.2. Precise Capture Rates on Tritium including Sub-Dominant Cosmological Effects
In
Table 6, non-thermal distortions during neutrino decoupling enhance the number density of the C
B by about
. To properly incorporate such effects into the capture rates of the C
B on tritium, we discuss the formula of their capture rate with
precision.
Cosmic relic neutrinos can be captured on tritium by the following inverse beta decay process,
The total capture rate for the C
B in this process,
, can be written
where
is the number of (mass) species of neutrinos.
is the capture rate for a given mass-eigenstate of neutrino
, given by
where
is the number of tritium,
is the total tritium mass in the experimental setup, and
is the atomic mass of tritium.
and
are helicity, velocity and the total cross section in the inverse beta decay on tritium, respectively.
is the local momentum distribution for relic cosmic neutrinos around the Earth, which satisfies
.
In cosmic neutrino capture on tritium, the spins of the outgoing electron and nucleus would not be measured. In addition, the spin of the initial nucleus would not be identified either. On the other hand, the helicity state for cosmic neutrinos in the Dirac case is polarized as in
Section 4.3. Then we compute the spin-polarized cross section for
. After averaging over the spin of
and summing over the spin of outgoing
and
, the formulae of
with
precision reduces to (see
Appendix D for detail calculations)
where
is a component of the Cabibbo-Kobayashi-Maskawa (CKM) matrix,
and
are the nuclear masses of
and
,
and
are the axial and vector coupling constant, and
and
are the reduced matrix elements of the Fermi and Gamow-Teller (GT) operators, respectively. The Fermi function
is an enhancement factor by the Coulombic attraction of the outgoing electron and proton, which is approximately given by [
64]
where
is the fine structure constant.
Z is the atomic number of the daughter nucleus and
for
. The energy and momentum for an emitted electron
and
depend on the neutrino masses and momenta strictly because of momentum conservation in the inverse
-decay process. However, since the contributions of the neutrino masses and momenta to
and
are very small,
and
are approximately given by (see
Appendix C for details)
where
is the beta decay endpoint kinetic energy for massless neutrinos given by
is so small compared to
and
that we can safely neglect
in Equation (92).
Then we obtain
with
precision substituting Equation (90) into Equation (89),
where
is the (unnormalized) average magnitude of velocity for
given by
Typically, contributes more than to . If , due to , we can drop in the formula of Equation (94) with precision. Here is the average momentum of the CB in the current universe. We also comment on whether we can use further approximations with precision to write Equation (94) into a simpler form. For massless neutrinos, due to , the (unnormalized) velocity is written as . For non-relativistic neutrinos , due to , is approximately written as , where and . We note that gravitational helicity flipping for massive neutrinos by neutrino clustering would be negligible since the helicity-dependent part in is already suppressed by .
5.2.1. Majorana vs. Dirac Neutrinos
For non-relativistic neutrinos, i.e.,
, if we set
in Equation (94),
is porportional to
and left-helical and right-helical components for relic neutrinos interact with tritium with the same magnitude via the weak interaction. Then the capture rate on tritium for Majorana neutrinos
is twice that for Dirac neutrinos [
34],
On the other hand, for relativistic neutrinos, i.e.,
, only the left-helical neutrinos interact with tritium via the weak interaction since helicity coincides with chirality in the relativistic limit. Then in both Majorana and Dirac cases, the capture rates are the same [
35],
Note again that the approximations in Equations (96) and (97) might not be valid for the capture rates with precision. To estimate the capture rates with precision, the term that depends on in Equation (94) should be included precisely.
5.2.2. Values of the Capture Rates on Tritium with
For references, we show values of the capture rates including cosmological effects discussed in
Section 4.2 and
Section 5.1 in the case of
. We choose other neutrino masses and their ordering to satisfy the observed values of neutrino squared-mass differences from neutrino oscillation experiments [
51,
52],
In both neutrino mass ordering we take the following values of the PMNS matrix,
Note that neutrino squared-mass differences and neutrino mixing parameters currently include a few percent (about ) uncertainties even at () confidence level.
In
Table 9, we show values of the capture rates on 100 g of tritium in both the cases of NO and IO for Majorana and Dirac neutrinos with
.
denotes the differences between the cases with and without effects of
-annihilation during neutrino decoupling and
denotes the differences with and without gravitational clustering for relic neutrinos in nearby galaxies.
For Majorana neutrinos, the capture rates for the first and second heaviest neutrinos are slightly less than twice those for Dirac neutrinos because of . On the other hand, the capture rates for massless (or almost massless) neutrinos in the cases of Majorana and Dirac neutrinos are the same because of .
5.2.3. Discussions on Exposure and Uncertainties in the Capture Rates
In this section, we discuss the required amount of tritium to observe the sub-leading cosmological effects themselves, , and the estimated error of the capture rates for relic neutrinos on tritium in more detail.
To observe
, we need a large number of events to satisfy typically
where
is the exposure time and
is a background rate. Even if the background is successfully removed, we need 10
–10
events of the C
B signal (
10
–10
) because of
for
. This requirement corresponds to the need for 10–10
kg yr of exposure of tritium. Currently, it is extremely difficult to obtain such amount of the exposure. In the next
Section 5.3, we comment on
-decay background, which is one of main background in cosmic neutrino capture on tritium.
The estimated error of the neutrino capture rates mainly comes from the uncertainties of the neutrino mixing parameter,
, and the undetermined value of the lightest neutrino mass,
. The current errors of PMNS matrix are about a few percent (about
) at
confidence level [
51,
52]. The current upper bound of
is ≲0.8
[
65]. Thus, unfortunately, it is still difficult to incorporate cosmological sub-dominant contributions into the value of
precisely. However,
for
is correctly estimated since uncertainties of
are canceled out in
. Future neutrino oscillation experiments will reduce uncertainties of PMNS matrix (see, e.g., [
66,
67,
68]). In addition, measurement of large
-decay background in the PTOLEMY-type experiment might determine the value of
very precisely [
31].
We also note that the theoretical calculation of
still includes the uncertainty of a few %, although the estimation of
through the observation of the tritium half-life and the value of the Fermi operator,
, only involves uncertainty of
[
69].
For a large value of , gravitational clustering effects of relic neutrinos are typically more dominant than effects of -annihilation during neutrino decoupling. Although the CB itself with a large value of would be easier to observe due to a large gravitational clustering, it is also a very difficult task to distinguish the effects of -annihilation during neutrino decoupling from gravitational clustering effect of relic neutrinos.
Based on the evaluation in this section, it is still extremely difficult to observe -annihilation during neutrino decoupling in the PTOLEMY-type experiment. However, the precise capture rates including cosmological sub-dominant contributions might be useful to distinguish the SM from physics beyond the SM properly in the future.
5.3. -Decay Background and the Energy Resolution of the Detector to Distinguish the CB Signal from It
Finally we comment on -background and the required energy resolution of the detector to distinguish the CB signal from this background, which is one of main difficulties to observe the CB directly in the inverse -decay process.
The main background comes from tritium
-decay process,
The
-decay spectrum and the capture rate for the
-decay process are given by [
70] (see also
Appendix D)
where
is the maximal energy of the emitted electron for
, where the electron is emitted in opposite direction to both
and
(see also
Appendix C),
Then the maximal energy for the emitted electron in the
-decay process called the energy at
-decay endpoint is
where
is the lightest neutrino mass. We can see that the
-decay spectrum
vanishes for
. Then the total tritium
-decay rate is obtained as
Since the event number of -decay background is extremely larger than that of the CB signal, we must distinguish the two signals clearly.
To distinguish the C
B signal and
-decay background, we need a tiny energy resolution of the detector
. The energy resolution of a detector characterizes the smallest separation where two signals can be distinguished. The
-decay background closest to the C
B signal is the electron signal with the maximal energy
. To distinguish the C
B signal for a mass species
from
-decay background near the endpoint, the required energy resolution
is expected to be (see
Appendix C for details)
where
is the emitted electron energy from the C
B signal,
given by Equation (92).
To take into account the energy resolution of the detector
in the spectrum and the number of events for the C
B signal and the
-decay background, we model the would-be observed spectrum of the emitted electron as a Gaussian-smeared version of the actual spectrum. This is achieved by convolving both the C
B signal and the
-decay background with a Gaussian of full width at half maximum (FWHM) equal to
, where
is the Gaussian standard deviation,
Substituting Equation (88) into Equation (108), the smeared spectrum of the emitted electron from the C
B signal can be written as
where
Equation (110) is a Fredholm integral equation of the first kind and is a would-be observed quantity. After solving Equation (110) inversely, the spectrum of the CB, , can be in principle reconstructed though we might need a significantly large number of observations for the CB events. We leave the detailed study for the reconstruction of the CB spectrum on tritium as future work.
In
Figure 5, we show the expected spectra for the emitted electrons from the C
B signals (solid lines) and the
-decay background (dashed lines) with
and 100 g of tritium, the energy resolution
(left panel) and
(right panel) considering the case of Dirac neutrinos and both the normal (fine red) and inverted (bold blue) mass hierarchies. In these figures, we neglect spectral distortions for the C
B from
-annihilation during their neutrino decoupling and the gravitational clustering for simplicity. We can see that the C
B signal is distinguished from the
-decay background if
. It is easier to distinguish the C
B signal from the
-decay background in the inverted mass ordering than the normal ordering. This is because we can obtain a larger number of events for the heaviest neutrinos in the inverted case due to the large value of
. In addition,
-decay spectrum near the endpoint is smaller in the inverted case because in the inverted case the
-decay spectrum near the endpoint is composed of
with small
while in the normal ordering that is composed of
with large
.
6. Conclusions
In the near future, CMB-S4 will determine with a very good precision of ∼0.03 at C.L., and consequently confirm neutrino decoupling process in the SM and/or impose severe constraints on many scenarios in physics beyond the SM. In addition, in the future, a direct observation of the CB might bring us more information about the early universe and neutrino physics. In both observations, the CB spectrum is one of crucial ingredients to estimate and a direct detection rate.
In this article, we review the formula of kinetic equations for neutrinos in the early universe, which are the quantum Boltzmann equations for neutrinos and the continuity equation and the possible spectral distortions due to -annihilation in neutrino decoupling. We also discuss the impact of the distortion of the CB spectrum in neutrino decoupling on direct observation of the CB on tritium, with emphasis on the PTOLEMY-type experiment.
We find
[
25,
26,
27] by solving the kinetic equations for neutrino density matrix in the early universe, including vacuum three-flavor oscillations, oscillations in
-background, finite temperature corrections to
,
and
P up to the next-to-leading order
(see also ref. [
24] for the first suggestion on the importance of this contribution), and the collision term where we consider full diagonal parts and off-diagonal parts derived from charged current interactions but neglect off-diagonal parts derived from neutral current interactions. Later, the authors in refs. [
26,
27] also find
and
, respectively, including off-diagonal parts in the collision term derived from neutrino neutral current interactions. Effects of their off-diagonal parts, and the choice of neutrino mass and mixing parameters on
are quite small,
[
27]. In refs. [
25,
26,
27], the Dirac CP-violating phase in neutrino mixing parameters is neglected. This contribution to
is expected to be also quite small since increases and decreases for the energy densities of neutrinos and anti-neutrinos due to the Dirac CP-violating phase would be canceled out (see also ref. [
54]). However, QED corrections to weak interaction rates at order
and forward scattering of neutrinos via their self-interactions have not been precisely taken into account. Recent studies [
23,
28] suggest that these neglects might still induce uncertainties of
in
. Although we should consider their contributions to
in the future,
is still a very good reference value.
We have revealed the spectrum, number and energy density of the C
B in the current homogeneous and isotropic universe, including the spectral distortions in neutrino decoupling, as in the right panel of
Figure 4 and
Table 4 and
Table 5. Then we have discussed the capture rates of the C
B on tritium with
precision to observe effects of
enhancement of the number density of the C
B by the spectral distortions due to
-annihilation during neutrino decoupling. Unfortunately, it is extremely difficult to observe such sub-dominant effects since we will need more than 10 kg of tritium. The precise capture rates of the C
B on tritium will be also useful to distinguish the SM from physics beyond the SM properly.
If observations and theoretical estimations of the C
B spectrum are improved significantly, we will obtain much richer information about neutrino physics and the early universe. Through direct observations of the C
B, one can impose significant constraints on neutrino decays and lifetimes in the region of the age of the universe,
[
34,
71]. The C
B spectrum would also have fluctuations imprinted by inflationary perturbations. Towards a precise estimation of anisotropy of the C
B as the CMB, one would need to solve kinetic equations for neutrinos in an anisotropic background, develop a detection method of the anisotropy, and reduce uncertainties of physical constants such as neutrino mass and mixing parameters, and Newton constant.