Here we review the physics of solar neutrinos, with neutrino production, neutrino flavor mixing which is modified by matter interactions and thermal motion of nuclei, and finally the neutrino energy spectra of different nuclear reactions. A simple analysis of supernova neutrino spectra shows that neutrino spectra for supernova are not well understood.
Neutrinos and antineutrinos are particles produced in many nuclear reactions such as beta decays,
In such reactions, charged current weak interaction plays a role which takes a down quark in neutron to a up quark while releasing electrons and antielectron neutrino,
More generally, positron emission and electron capture are also neutrino related nuclear reactions which is explained in table Neutrino related nuclear or leptonic reactions.. In the context of astrophysics, (anti)neutrinos also participate in nuclear reaction chains in stars, synthesis of heavy and rare elements and more.
LaTeX Code for The Table
\begin{table}[ht]
\centering
\begin{tabular}{c  c  c}
\hline
Reaction & Equation & Boson \\ [0.5ex]
\hline
Electron emission & :math:`{}^A_Z X \to {}^A_{Z+1}X + e^ +\bar \nu_e` & :math:`W` \\
Positron emission & :math:`{}^A_Z X \to {}^A_{Z1}X + e^+ + \nu_e` & :math:`W` \\
Electron capture & :math:`{}^A_Z X + e^ \to {}^A_{Z1}X + \nu_e` & :math:`W` \\
Positron capture & :math:`{}^A_Z X + e^+ \to {}^A_{Z+1}X + \bar\nu_e` & :math:`W` \\
[0.5ex]
\hline
Electron annihilation & :math:`e^ + e^+ \to \nu_e + \bar\nu_e ` & :math:`W` \\
Electron annihilation & :math:`e^ + e^+ \to \nu + \bar\nu ` & :math:`Z` \\
[0.5ex]
\hline
Neutrino capture & :math:`{}^A_{Z}X + \overset{()}{\nu_e} \to {}^A_{Z\mp 1}X + e^\pm ` & W\\
[1ex]
\hline
:math:`e^\nu` scattering & :math:`e^ + \overset{()}{\nu_e} \to e^ + \overset{()}{\nu_e} ` & :math:`W` \\
:math:`e^\nu` scattering & :math:`e^{\pm} + \overset{()}{\nu_e} \to e^{\pm} + \overset{()}{\nu_e} ` & :math:`Z` \\
Neutrino scattering & :math:` {}^A_Z X + \overset{()}{\nu} \to {}^A_Z X + \overset{()}{\nu} ` & Z\\
[0.5ex]
\hline
\end{tabular}
\caption{Neutrino related nuclear or leptonic reactions}
\end{table}
Astrophysical neutrino sources such as cores of stars, AGN, Gammaray bursts, and supernovae, reveals a lot of information about the sources. Though neutrinos are are weakly interacting particles thus rendering them hard to detect, it is an important subject to theoretically inspect neutrino productions in astrophysical processes.
One of the neutrino factories is the stellar core. Numerous nuclear reactions as well thermal neutrino production produce high number density neutrinos emissions. In this section, the nuclear reactions in stars are reviewed as well as the neutrino oscillation in the stars.
pp chain and CNO cycle are the most important nuclear reactions in a star. Through different stages of its life, a star experience different nuclear reactions. Fig. 14.3 shows the dominant energy source of solar mass stars. In order to calculate the neutrino spectrum we need the neutrino production rate in each reaction and the branching ratios. The processes that emits neutrinos are the red reactions in Fig. 14.3, which explicitly are
LaTeX Code for The Diagram
\begin {figure*}%[!hbtp]
\centering
\begin{adjustbox}{width=\textwidth}
\begin{tikzpicture}[sibling distance=15em,
every node/.style = {shape=rectangle,
draw, align=center}
edge from parent/.style = {draw, latex},]]
\node {\color{red}:math:`\mathrm{p+p\to {}^2H + e^+ +\nu_e}` }
child { node {:math:`\mathrm{p+{}^2H \to {}^3He + \gamma}`}
child { node {:math:`\mathrm{{}^3He+{}^3He \to {}^4 He + 2p }`}
edge from parent node [left] {83.30\% (ppI) } }
child { node {:math:`\mathrm{{}^3He+{}^4He \to {}^7 Be + \gamma }`}
child { node {
\color{red}:math:`\mathrm{{}^7Be + e^ \to {}^7Li + \nu_e}`
}
child { node { :math:`\mathrm{{}^7Li + p \to 2{}^4He }`} }
edge from parent node [left] {99.88\% (ppII) } }
child { node { :math:`\mathrm{{}^7 Be + p \to {}^8 B + \gamma}`}
child { node { \color{red}:math:`\mathrm{{}^8B \to {}^8Be^* +e^+ +\nu_e}` }
child { node { :math:`\mathrm{{}^8Be^* \to 2 {}^4He }` } }}
edge from parent node [right] {0.12\% (ppIII) } }
edge from parent node [right] {16.70\% } }};
\end{tikzpicture}
\end{adjustbox}
\caption{ pp chain with branching ratio. From Altmann 2001. }
\end{figure*}
Solar neutrinos are mostly produced in pp reaction, Be electron capture and B decay, which are called pp neutrinos, Be neutrinos and B neutrinos. Even without knowledge of the detailed reactions, the conservation of lepton numbers will lead to the overall neutrino production
where it is important to notice that two neutrinos are produced in each reaction. This overall reaction can either be pp chain or CNO cycle.
Using this simple relation, we can estimate the neutrino flux emitted by our sun. Given the kinetic energy produced in each reaction is the difference between the initial masses and the final masses, \(Q_{pp}=4m_p+2m_em_{He4}=26.7\mathrm{MeV}\) where the mass of neutrinos are neglected since they are small compared with every other particle. On average, each neutrino carries away 0.2MeV energy and the rest will be mostly in the form of thermal energy \(Q_\gamma=26.3\mathrm{MeV}\). Number flux of thermal photons near Earth can be calculated using the solar constant \(S_0\),
Since we know each reaction produces 2 neutrinos while producing \(Q_\gamma\), which means that the number flux of neutrinos near Earth is roughly twice of the number flux of photons, i.e.,
For such a large flux, understanding the role in stellar nuclear reaction and spectra is important.
Inside our Sun, two additional reactions also produce neutrinos which are called pep and hep neutrinos.
pep neutrinos are produced in
which is only has a branching ratio 0.4% instead of the 99.6% of pp reaction.
hep neutrinos are produced in
which has a branching ratio of \(2\times 10^{5}\%\). As a comparison, the \(\mathrm{{}^3He + {}^3He}\) has a branching ratio 85% and \(\mathrm{{}^3He + {}^4He}\) has a branching ratio 15%.
Neutrinos are special particles that their flavor eigenstates are not the propagation eigenstates, which leads to neutrino flavor oscillations. Since neutrinos with different flavor interact with matter with different cross section, we need to investigate the neutrino flavor carefully. Even though only electron flavor neutrinos are produced, what we detect on Earth is different in flavor, which depends on two phenomena, neutrino vacuum oscillation and Mikheyev–Smirnov–Wolfenstein effect.
To understand the neutrino vacuum oscillation phenomenon, we use two flavor neutrino as an example. In vacuum, propagation states are mass eigenstates, which is different from flavor eigenstates. The wave function in flavor eigenstates basis is related to wave function in mass eigenstates through an unitary matrix \(\mathbf U\),
where \(\Psi_f\) is wave function in flavor basis and \(\Psi_v\) is the wave function in vacuum mass eigenstate basis. The rotation matrix is
In vacuum basis, the Hamiltonian is free propagation, which is given by
where
We assume the neutrinos have almost the same momentum, which is true since their mass is small, i.e., \(p_i \approx E\). To first order, the Hamiltonian becomes
where the identity matrices only give us an overall phase so we drop them. With the definition that \(\Delta m^2 = m_2^2  m_1^2\) The vacuum Hamiltonian in vacuum basis is simplify
which leads to the simple solution for the wave function in vacuum basis
where the initial condition is
In flavor basis, the wave function at anytime is
As seen in the nuclear reactions in the solar core, electron neutrinos are most abundant flavor. Initial condition is assumed to be electron flavor in the calculation which leads to the survival probability of electron flavor
Since neutrinos travel with velocity approximately the speed of light, we use \(L = t\) where L is the distance travelled. The survival probability is
The important parameter is the oscillation length of the neutrino flavor conversion. Here we have the oscillation frequency \(\omega = \frac{\Delta m^2}{2E}\).
Oscillation Length
Need a figure about the oscillation lengh here.
The nature of neutrino oscillation means that flavor conversion occurs as long as their propagation eigenstates are not flavor eigenstates. We expect neutrino propagation eigenstates in matter are different from flavor states in general. [wolf78] Using the fact that neutral current interactions between different flavor neutrinos and matter is independent of flavor, we only include the charged current, which will produce a effective potential for electron flavor. In flavor basis, the effective potential is
where \(G_F\) is Fermi constant, \(n_e\) is number density of electrons. We also removed the identity in this matrix since it doesn’t change our survival probability. The Hamiltonian with matter effect is the combination of vacuum oscillation and matter effect, which is, in flavor basis, explicitly,
where we used the result of flavor basis vacuum oscillation Hamiltonian
Applying Pauli matrices and \(\lambda = \frac{\sqrt{2}G_F n_e}{2}\) to this total Hamiltonian, it is rewritten as
Due to the offdiagonal terms in the Hamiltonian, the system will experience oscillations in flavor. A resonance, i.e., maximum mixing, dominates the system when the diagonal terms becomes zero,
which gives us the MSW resonance condition.
The importance of matter effect to our understanding of solar neutrinos is that it modifies the oscillation, which depends on the matter profile. For a solar mass star, we have almost adiabatic evolution of the neutrinos, which means that the instantaneous eigenstates and eigenvectors of Hamiltonian is good enough for the time dependent Schrodinger equation.
For simplicity, we define the vacuum frequency and the hatted quantities
The eigenstates, derived by diagonalizing the Hamiltonian, are
Fig. 14.5 shows the two energy levels. For very high matter density, which interact with electron neutrinos more through charged current, electron flavor is composed almost with heavy propagation eigenstate. However, as the matter density becomes lower, the heavy propagation state will be gradually transformed to the other flavors since electron flavor in vacuum is composed mostly the light mass state. The resonance, which is the closest point of energy levels, happens at density \(n_e = 2\omega \cos(2\theta_v)/\sqrt{2}G_F\) which depends on \(\omega = \Delta m^2/2E\). Neutrinos with different energies have different resonance, which will significantly reshape the neutrino spectra for different flavors. More explicitly, the conversion of flavor is shown in Fig. 14.6 which is taken from Smirnov. [Smirnov2003] Since we are discussion adiabatic evolution, the probability of each energy eigenstates doesn’t change, as the boxes of each energy eigenstates shown in Fig. 14.6 doesn’t change in size. The extreme dense case shows that matter converts a lot of electron flavor to muon flavor.
In summary, even though only electron flavor neutrinos are produced in the core of a solar mass star, the neutrino flavor conversion to the other flavors is enhanced by matter interaction, in addition to the vacuum oscillation. However, the actual neutrino flavor conversion is much more complicated than just MSW effect and almost impossible to calculate without knowing the very exact matter profile of the Sun even with the timedependent small perturbations to the density. As an approximation, MSW transition is good enough for the solar neutrinos. [Lopes2013a] In the case of the sun, the flavor conversion is calculated by Ilídio Lopes and shown in Fig. 14.7. [Lopes2013]
The simplest model for the total neutrino flux generated due to solar nuclear reaction is simple linear superposition of solar neutrino flux from each reactions, since solar neutrinos are not dense enough to have a significant selfinteraction. For the total neutrino flux F,
where \({}_{\alpha}\) stands for flavor and E is the energy of neutrinos. To generate the final solar neutrino spectra, the neutrino energy for each reaction should be examined.
For each reaction, the spectral shape of solar neutrinos are mostly spectrum inferred from lab experiments, with modifications of thermal motion and slight modifications from relativistic effect. [Bahcall1991]
In the lab experiments, we only have low energy nuclei and less dens environment, while in the solar core, the temperature is high and the density is large. The first question to ask is whether the neutrino production is modified by a thermal equilibrium environment or it is too fast that no equilibrium can not be maintained by the electromagnetic scattering. The answer is that the neutrino production is in equilibrium which ensures the equilibrium statistics of the neutrino spectra. [Bahcall1991] The two quantities that is related to this problem is the characteristic time of neutrino production and the characteristic electromagnetic scattering. The time scale of electromagnetic scattering is [Bahcall1991]
where E is the energy of ions which is in a bath of p plasma and \(\rho\) is the density of plasma. For the sun, this is of order \(10^{12}\mathrm{s}\). Meanwhile, the time scale of neutrino production \(\tau_{\nu}\) is shown in table Neutrino production characteristic time. Reproduced from J. N. Bahcall, “Shapes of solarneutrino spectra: Unconventional tests of the standard electroweak model,” (1991)..
LaTeX Code for The Table
\begin{table}[ht]
\centering
\begin{tabular}{cc}
\hline
Nuclear Reaction & :math:`\tau_{\nu}` \\
\hline
pp & :math:`10^{10}` years \\
:math:`\mathrm{ {}^7Be }` & :math:`10^{12}` years \\
:math:`\mathrm{ {}^8B }` & :math:`1` second \\
:math:`\mathrm{ {}^{13}N }` & :math:`10^3` seconds \\
:math:`\mathrm{ {}^{15}O }` & :math:`10^2` seconds \\
:math:`\mathrm{ {}^{17}F }` & :math:`10^2` seconds\\
\hline
\end{tabular}
\caption{Neutrino production characteristic time.}
\end{table}
This huge difference between \(\tau_{EM}\) and \(\tau_{neutrino}\) keeps the thermal equilibrium of the ions that produces neutrinos.
Therefore we consider the neutrino spectra of all the reactions with thermal motions of the nuclei in a thermal bath. John Bahcall explained thermal corrections to neutrino spectra of beta decay using a simple argument that [Bahcall1991]
where \(v_z\) is the velocity that causes the spectrum redshift, i.e., recoil velocity, and the integral should be over \(v_z\) not \(q'\). \(F_{lab}\) is the spectrum in lab experiments. \(f(v_z)\) is the corrections due to thermal motion of the nuclei. The relation between \(dq\) and \(dq'\) is set up using relativistic velocity transformation
The final result for the spectrum is
The upper limit of the integral is given by the cut off of momentum \(q'_m\),
where \(q'_m:math:\) is given by the Q values of the nuclear reactions. The maximum momentum is limited by the energy released in the reaction.
Since the core temperature is not high enough to make the heavy nuclei relativistic, we expand the spectrum using Taylor expansion and keep only first order,
Thus the correction depends on the energy of neutrinos. The correction only plays a role when the energy of neutrino is close to endpoint energy of the beta decay. The energy of the neutrinos will be blueshifted beyond the beta decay spectrum endpoint in lab experiments. However, the correction is as small as \(10^{6}\) of the lab spectrum in the case of solar core. [Bahcall1991]
The other correction is due to the gravitational redshift. Combining the thermal correction and gravitational redshift, we get the neutrino spectrum that should be observed on the Earth. As an example, John Bahcall calculated the pp reaction spectrum, which is given in Fig. 14.9.
This overall spectra are the summation of neutrinos produced at different radius where the temperature and density are different. One of them is the fact that nuclear reactions happen at different rates when the temperature changes. Adelberger et al calculated the stellar energy contribution of pp chain and CN cycle at different temperatures. [Adelberger2011]
Standard solar model produces a complete set of neutrino spectra from different reactions, which is shown in Fig. 14.11. The total neutrino spectrum we expect from the sun is the superposition of all the neutrino spectra from different reactions. From Fig. 14.10 we expect the pp chain neutrino flux is \(2\sim 3\) orders of magnitude larger than CNO cycle neutrino flux, which is checked in Fig. 14.11. Meanwhile, another interesting question to look into is the neutrino production at different radius. Since nuclear reaction rates depend on temperature and density, neutrino flux generated at different radius inside the Sun is very different. A calculation done by Ilídio Lopes shows that most neutrinos are produced at 0.05 radius of the Sun, which is in Fig. 14.12. [Lopes2013]
The solar neutrino spectra is not very different from lab experiments since solar neutrino flux is not high enough to interact with solar medium significantly. However, in a supernova explosion, \(10^{58}\) neutrinos are released from the protoneutron star, which is of radius \(10\mathrm{km}\), in a few seconds. The huge number density of neutrinos and large density of matter have a huge effect on the neutrino spectra, especially for different flavors since matter has a huge effect on flavor conversion as we have already seen in MSW effect. The matter effect will be much more than MSW since the violent matter motion. In addition, neutrino neutrino interaction will be efficient because of the high neutrino number density.
Apart from the emission of neutrinos from nuclear reactions of electron capture and positron emission in the solar interior, supernova environment also gives rise to Bremsstrahlung pair neutrino production, electronpositron neutrino pair production, which brings all three flavors and also antineutrinos into the spectra. However, even the with the presence of intensive interaction between neutrinos and the leptons and hadrons, which thermalize the neutrinos in the supernova core, the neutrino spectrum escaping from the supernova core is not completely FermiDirac spectrum. Nonetheless, it is possible to parametrize it using nominal FermiDirac distribution, [raffelt2004]
However, numerical results show that there is a deviation from this FermiDirac distribution. [Totani1998] [Keil2003] Keil Mathias and Georg Raffelt showed that it is good enough to approximate the neutrino spectrum from supernova in Monte Carlo simulations using the so called “alpha fit”,
where \(\langle E\rangle\) is the average energy, or the first moment of energy. The values from Monte Carlo simulations falls into the range \(\alpha = 2.5\sim 5\), which clearly shows the spectra are pinched as explained in Fig. 14.13. Detection of deviation from nominal FermiDirac distribution will show evidence of corecollapse information.
In this review, we presented the solar neutrino production, thermal modification, gravitational effect and flavor conversion, which leads to the theoretical predicted solar neutrino spectra for each reaction. Even though we understand solar neutrino well, overall the neutrino spectra of supernova are not so to our complete knowledge. Phenomena such as spectral split due to neutrinoneutrino interaction and matter effect reshape the spectra significantly. Future supernova neutrino observation data is needed to a better understanding of the supernova physics.
[Smirnov2003] 

[Adelberger2011] 

[wolf78] 

[Lopes2013a] 

[Lopes2013]  (1, 2, 3)

[raffelt2004] 

[Bahcall1991]  (1, 2, 3, 4, 5)

[Totani1998] 

[Keil2003]  0M. T. Keil, Supernova Neutrino Spectra and Applications to Flavor Oscillations, Ph.D. thesis (2003), arXiv:0308228 [astroph]. 
© 2021, Lei Ma  Created with Sphinx and .  On GitHub  Physics Notebook Statistical Mechanics Notebook Index  Page Source