Phenomenology of (3+1) Four Flavour Neutrino Oscillations in Dense Matter

Ashish Shrivastava1, Bipin Singh Koranga2,* and Seema Rawat3

1Department of Applied Physics, Delhi Technological University, Delhi-110042, India
2Department of Physics, Kirori Mal College, University of Delhi, Delhi-110007, India
3Department of Physics, Zakir Husain Delhi College, University of Delhi, Delhi-110007, India
E-mail: bskoranga@kmc.du.ac.in
*Corresponding Author

Received 11 August 2024; Accepted 03 November 2024

Abstract

The mass square difference varies as neutrinos interact with the matter, assuming four flavours of neutrino oscillations. Three regions – vacuum-dominated, resonance-dominated, and matter-dominated – are where the neutrinos oscillate. This research goes on to examine how each of these three locations affects the likelihood of neutrino oscillations and the mass square difference in distinct ways. There are three alternative analytical formulas for the mass square difference. In this study, we have studied the behaviour of neutrino oscillations for four flavours in vacuum-dominated, resonance-dominated, and matter-dominated regions.

Keywords: Four flavor neutrino oscillations, sterile neutrino, dense matter, mass square difference.

1 Introduction

With the help of the findings of numerous solar, atmospheric, accelerator, and reactor experiments [13], neutrino oscillation has been very well defined. Since the discovery of neutrino oscillations, the behaviour of the particles and their propensities to oscillate between various flavours have been examined. Since the mixing angle and the mass square difference affect likelihood. It was discovered that these are impacted by the matter effect, which occurs when neutrinos contact with matter. Although the precise mass of neutrinos is still unknown, multiple medium and long baseline experiments [4, 5] have provided the values of the mass square difference of various mass eigenstates and mixing angles.

The solar neutrino mass square difference measured by the LSND experiment is 0.2<Δ21<100eV2 [6], whereas the solar mass square difference measured by the MiniBooNe experiment is 0.01<Δ21<1.0eV2 [7]. This experimental finding indicated the potential for sterile (3+1) neutrinos [4, 5]. Four flavour neutrino oscillation is conceivable if sterile neutrinos are predicted to exist. When neutrinos interact with matter through neutral and charged current interactions, a rise in potential is seen, which affects the oscillation probability and mass square difference [8]. The outcomes of these interactions are determined by the neutrino energy as well as the matter density, or the concentration of electrons and neutrons in the matter. The three regions that make up the matter effect are the resonance region (also known as the MSW effect), the matter dominated region (also known as AΔ21,|Δ31| [12], and the vacuum dominated region (also known as AΔ21,|Δ31|,Δ41.

2 Framework for Four-Flavor Neutrino Oscillation in Theory

Four-flavor mixing can be expressed as [9]: in terms of mass eigenstates.

(vevμvτvs)=U(v1v2v3v4)

where

U=R34(θ34,δ34)R24(θ34)R14(θ14,δ14)
R23(θ23)R13(θ13,δ13)R12(θ12)
Rij(θij,δ)=(cosθijsinθije-iδ-sinθijeiδcosθij).

Here, δ represents dirac CP violating phase.

The following can be used to represent a Hamiltonian in four flavours in a vacuum [10]

Hv=12EU(00000Δm2120000Δm3120000Δm412)U.

Thus the likelihood of neutrino oscillation can be expressed as [11]:

Pvαvβ=δαβ-4i>j(UβiUαi*Uβj*Uαj)+2i>j(UβiUαi*Uβj*Uαj)

We already know that the hamiltonian changes as neutrinos go through matter because their potential increases. This increase is the result of neutral current scattering and coherent forward elastic weak charged current. Charged current scattering is caused by the interaction of an electron and neutrino pair. Consequently, the updated Hamiltonian may be written as [10]:

Hm=Hv+A =12E(U(00000Δm2120000Δm3120000Δm412)U
+(ACC00000000000000ANC))

where

ACC=7.63×10-5ρEandANC=-3.815×10-5ρE

The matter terms for current charge and neutral charge, respectively, are ACCand ANC. We will utilise the standard hierarchical approach to tackle the issue for |Δ31|.

2.1 Vacuum Dominated Region (ACCΔ21,|Δ31|,Δ41)

As stated in the section’s headline, the matter effect in vacuum-dominated regions is small, meaning that interactions between neutrinos are insignificant and do not raise the effective potential of moving neutrinos. In light of this circumstance, we may therefore deduce the approximate values of the matter terms ACC0 and ANC0.

Effective Hamiltonian can be expressed as follows from these values:

Hm=Hv=12EU(00000Δm2120000Δm3120000Δm412)U

This, in vacuum, is equivalent to the Hamiltonian.

ij={0ijΔi1i=j

There is no such shift in mass square disparities in vacuum-dominated regions.

2.2 Resonance Region (AΔ21,|Δ31|)

The neutral and charge current in the resonance zone is sufficient to increase the neutrinos’ potential, which will modify their mass square difference and alter their chance of oscillations as well. The following is an example of an effective Hamiltonian for this area:

Hm=Hv+A =12E(U(00000Δm2120000Δm3120000Δm412)U
+(ACC00000000000000ANC))

This, as a result of coherent forward elastic weak charged current and neutral current scattering, closely resembles the effective hamiltonian. You can alternatively write the above form as:

UHmU=(H11H12H13H14H21H22H23H24H31H32H33H34H41H42H43H44)=

where:

ij={ACCUei*Uej+ANCUsi*UsjijΔi1+ACC|Uei|2+ANC|Usi|2i=j

Using Jacobi-Like Method [13] after the first rotation of the matrix :

Following the initial rotation of the matrix H, one can use the Jacobi-Like Method [13] to:

λ-=H11+H22tan2ω1-2A1tanω11+tan2ω1

and

λ+=H22+H11tan2ω1-2A1tanω11+tan2ω1.

When the expressions’ defined parameters can be expressed as:

tanω1=2A1(H22-H11)+(H22-H11)2+4A12,
A1=+|H12|

and

Φ1=Arg(sign(A1)H12*)

After the second rotation:

λ-=λ++H33tan2ω2-2A2tanω21+tan2ω2

and

λ+=H33+λ+tan2ω1-2A2tanω21+tan2ω2

When the parameters in the expressions have the following representations:

tanω2=2A2(H33-λ+)+(H33-λ+)2+4A22
A2=+|H23|,
Φ2=Arg(sign(A2)H23*)

and

H23=H13tanω1eiΦ1+H231+tan2ω1.

So now we can say that:

Δ21m=λ--λ-,
Δ31m=λ+-λ-

and

Δ41m=H44-λ-

images

Figure 1 Variation of relative percentage change in Mass Square difference in resonance region where energy varies from 0GeV to 0.5GeV where mixing angles are θ12=340θ13=10,0θ14=340,θ23=3.60,θ24=40 and θ34=18.50, mass square difference in vacuum are Δ21=8×10-5eV2, Δ31=2×10-3eV2 and Δ41=1.7eV2 and constant density 3g/cc.

images

Figure 2 Variation of relative percentage change in Mass Square difference in resonance region where energy varies from 0.5 GeV to 10 GeV where mixing angles are θ12=340θ13=100, θ14=340, θ23=3.60, θ24=40 and θ34=18.50, mass square difference in vacuum are Δ21=8×10-5eV2, Δ31=2×10-3eV2 and Δ41=1.7eV2 and constant density 3g/cc.

2.3 Matter Dominated Region(ACCΔ21,|Δ31)

In matter dominated region, neutrinos must be travelling in a dense matter. The condition mentioned in the heading of this section AΔ21,|Δ31| can be approximated as AΔi10 where i=2,3. Using this condition in the expression of resonance region, the effective hamiltonian can be written as:

Neutrinos must be passing through a dense matter in an area where matter predominates. The condition AΔ21,|Δ31| given in the section’s headline can be roughly expressed as AΔi10 where i=2,3. The effective Hamiltonian can be stated as follows using this condition in the formulation of resonance region:

UHmU=(H11H12H13H14H21H22H23H24H31H32H33H34H41H42H43H44)=

where:

ij={ACCUei*Uej+ANCUsi*UsjijACC(|Uei|2+ANCACC|Usi|2+Δi1Acci=j

Using a method akin to Jacobi [13]. The formulas obtained for the resonance zone will be the same for this as well, with the exception that the matrix’s parameters will change once the condition for various values of ij is applied. The graph below shows the variation in ACC for dense matter, i.e. AΔ21,|Δ31|, from less dense to very dense matter. Due to the fact that Δ21m=λ--λ-, Δ31m=λ+-λ- and Δ41m=H44-λ-. The connection H44-λ-+Δ21m can be obtained from the expression Δ41m. λ- and Δ21m do not affect the value of Δ41m at resonance region, but at high density, the values are comparable, therefore the change in Δ41m is also visible in the derived graph. This change may become saturated at larger densities, just like Δ21m and Δ31m.

images

Figure 3 Variation of relative percentage change in Mass Square difference in matter dominated region where density varies from 0 g/cc to 1500 g/cc where mixing angles are θ12=340θ13=100, θ14=340, θ23=3.60, θ24=40 and θ34=18.50, mass square difference in vacuum are Δ21=8×10-5eV2, Δ31=2×10-3eV2 and Δ41=1.7eV2 and constant energy 3GeV.

3 Conclusions

When neutrinos interact with matter, the probabilities of neutrino oscillation and the mass square difference are affected. Only in the case of neutrinos is the normal hierarchy computation completed. The probability in vacuum-dominated regions are nearly identical to those in vacuum, and the percentage change in the relative mass square difference with regard to matter is essentially identical to those in vacuum. The effective mass square difference is growing in the resonance zone as the neutrino energy rises and the oscillation probability changes. While there are no such impacts on Δ41m in this location, the percentage change in mass square difference is found to be in the range of 0 GeV to 10 GeV. However, there is a significant change in the resonance region from region 0 GeV to 0.5 GeV, followed by a harsh change in Δ21m, a significant change in Δ31m, as seen in the graph, and another drastic change after 3 GeV. The relative percentage change in the mass square difference with respect to the matter varies depending on whether the region is vacuum or resonance in a region where matter predominates.

Both the Δ21m and Δ31m alter and become saturated at about 99%. But even though the sterile neutrino only interacts with gravity, the change in Δ41m has begun to rise. According to the expressions derived from the relationship between Δ21m and Δ41m, the extreme shift in in a dense substance that caused the effective change observed has an impact on Δ41m%.

References

[1] John N. Bahcall et al., Phys. Rev. Lett. 40. pp. 1351–1354, 1978.

[2] Patrick Huber, Thomas Schwetz, Phys. Lett. B 669. pp. 294–300, 2008.

[3] Zeno Dixon Greenwood et al., Phys. Rev. D 53. pp. 6054–6064, 1996.

[4] G.J. Feldman et al., Advances in High Energy Physics, Vol. 2013, pp. 1–30, 2013-01, 2013.

[5] Chiaki Yanagisawa et al., Phys. Rev. D 83, pp. 072002, 2011.

[6] G. Conforto et al., Phys. Lett. B 447, pp. 122–126, 1999.

[7] MiniBoone Collaboration, A.A. Aguilar-Arevalo et al., Phys. Rev D 103, pp. 1–5, 2021.

[8] A. Yu. Smirnov, Phys. Scripta T 121, pp. 57–64, 2005.

[9] B. S. Koranga et al., Int. J. Theor. Phys. 60, 9 pp. 3548–3565, 2021.

[10] Ashish Shrivastava, B. S. Koranga et al., arxiv:2106.15945 (hep-ph).

[11] S. Dev et al., Nuclear. Phys B 491, pp. 401, 2019.

[12] S. Luo et al., Moscow Univ. Phys. Bull. 77 pp. 462–463, 2022.

[13] Bao-Biao Yue, Jia-Jie Ling et al., Chin. Phys. C 44, 10 pp. 103001, 2020.

Biographies

Ashish Shrivastava is a post graduate student in the Department of Applied Physics, Delhi Technological University. He is in the final year of M.Sc. Physics. He received his bachelor’s degree from Kirori Mal College, University of Delhi in Physical Science. He completed his dissertation in Phenomenology of Neutrino Oscillations during his under graduation. His work has been published in Journal of Graphic Era along with the conference. He is currently working on a project in which he is studying Quantum Entanglement in Neutrino Oscillations.

Bipin Singh Koranga is Professor in the Department of Physics, Kirorimal College, University of Delhi. He joined the Theoretical Physics Group at the Indian Institute of Technology Bombay in 2006 and earned his Ph.D. in physics (neutrino masses and mixings) there in 2007. For the past 17 years, he has been instructing graduate-level students in foundational physics and mathematical physics courses. The genesis of the cosmos, physics beyond the standard model, theoretical nuclear physics, quantum mechanical neutrino oscillation, and a few astrological issues are among his areas of interest. He has more than 60 scientific publications to his credit in different international journals. His current areas of focus in study are linked phenomenology and neutrino mass theories.

Seema Rawat is an Associate Professor of physics at Zakir Husain College, University of Delhi, where she has taught Quantum mechanics and Theortical Physics to undergraduate students. He received his master’s degree in physics from Kumaun University Nainital and doctoral degree in physics from Kumaun University Nainital. He has more than 30 scientific publications to his credit in different international journals. His current areas of focus in high energy particle physics.

Abstract

1 Introduction

2 Framework for Four-Flavor Neutrino Oscillation in Theory

2.1 Vacuum Dominated Region (ACCΔ21,|Δ31|,Δ41)

2.2 Resonance Region (AΔ21,|Δ31|)

images

images

2.3 Matter Dominated Region(ACCΔ21,|Δ31)

images

3 Conclusions

References

Biographies