×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

Neutrino oscillations in the Non-Kerr black hole with quantum phenomenon

Get Citation
Husan Alibekov, Farruh Atamurotov, Ahmadjon Abdujabbarov and Vokhid Khamidov. Neutrino oscillations in the Non-Kerr black hole with quantum phenomenon[J]. Chinese Physics C. doi: 10.1088/1674-1137/ad3c2c
Husan Alibekov, Farruh Atamurotov, Ahmadjon Abdujabbarov and Vokhid Khamidov. Neutrino oscillations in the Non-Kerr black hole with quantum phenomenon[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ad3c2c shu
Milestone
Received: 2024-01-07
Article Metric

Article Views(783)
PDF Downloads(14)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Neutrino oscillations in the Non-Kerr black hole with quantum phenomenon

    Corresponding author: Husan Alibekov, alibekov@astrin.uz
  • 1. Ulugh Beg Astronomical Institute, Astronomy St. 33, Tashkent 100052, Uzbekistan
  • 2. University of Tashkent for Applied Sciences, Str. Gavhar 1, Tashkent 100149, Uzbekistan
  • 3. New Uzbekistan University, Movarounnahr street 1, Tashkent 100000, Uzbekistan
  • 4. Tashkent State Technical University, Tashkent 100095, Uzbekistan
  • 5. Institute of Theoretical Physics, National University of Uzbekistan, Tashkent 100174, Uzbekistan
  • 6. Institute of Fundamental and Applied Research, National Research University TIIAME, Kori Niyoziy 39, Tashkent 100000, Uzbekistan
  • 7. Tashkent University of Information Technologies named after Muhammad al Khwarizmi, Amir Temur 108, Tashkent 100014, Uzbekistan

Abstract: In this study, we have investigated the mathematical components of the Dirac equation in curved spacetime and how they can be applied to the analysis of neutrino oscillations. More specifically, we have developed a method for calculating the phase shift in flavor neutrino oscillations by utilizing a Taylor series expansion of the action that takes into account $ \Delta m^4 $ orders. In addition, we have used this method to assess how the phase difference in neutrino mass eigenstates changes according to the gravitational field described by the Johannsen spacetime.

    HTML

    I.   INTRODUCTION
    • In the field of neutrino physics, a phenomenon called neutrino oscillation has been reported for which neutrinos change from one flavor to another while moving through space. This oscillation results from the interaction between the three recognized neutrino flavors: electron neutrino $ \nu_e $, muon neutrino $ \nu_\mu $, and tau neutrino $ \nu_\tau $. Pontecorvo [1] first proposed the concept of neutrino oscillation. This author also suggested that neutrinos possess a mass that was previously thought to be nonexistent. The mathematical description of neutrino oscillation involves the use of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [2], which connects the flavor states of neutrinos to their mass eigenstates. This matrix includes four parameters: three mixing angles and one phase. The mixing angles determine the likelihood of a neutrino transitioning from one flavor to another, while the phase influences the relative probabilities of oscillation between different flavors. Neutrino oscillation has significant implications in astrophysics, particle physics, and cosmology, facilitating our understanding of neutrino properties and their role in the universe [3, 4], including their contribution to dark matter. It is a captivating phenomenon that has opened up new avenues of research in particle physics and astrophysics [5]. The discovery of neutrino oscillation has challenged our knowledge on neutrinos and provided valuable insights into the nature of the universe The significant breakthrough of this discovery has received recognition through various awards, notably the 2015 Nobel Prize in Physics. This prestigious honor was bestowed upon Takaaki Kajita and Arthur B. McDonald for their remarkable contributions to the Super-Kamiokande [6] and SNO experiments.

      Sudbury Neutrino Observatory (SNO) [7], The Super-Kamiokande [8], and MINOS experiments [9, 10] played a vital role in enabling the discovery of neutrino oscillation. These experiments observed the phenomenon of different types of neutrinos disappearing and reappearing as they passed through the Earth's atmosphere or matter.

      Different research studies aim to understand the phenomenon of neutrino lensing caused by gravitational sources. They are revealing an intriguing connection between the probability of neutrino oscillation and the individual masses of neutrinos. This connection is explained through an analysis of the impact of weak lensing induced by a Schwarzschild mass [11]. This analysis explores the implications of gravitationally modified neutrino oscillations in realistic scenarios involving two or three flavors, such as the influence of the gravitational field of a a supernova on the travel of emitted neutrinos. This influence could have observable effects on the neutrino signal [12]. Furthermore, the propagation of neutrinos in a strong gravitational field regime has been studied by considering electromagnetic interactions using the WKB approximation [13]. The behavior of neutrino oscillations in the Schwarzschild spacetime has also been investigated taking into account spin precession in the presence of a magnetic field [14, 15]. Notably, both radial and nonradial propagations of neutrinos in the Schwarzschild spacetime have been examined [16]. Finally, the effects of universe expansion and torsion on neutrino oscillations have also been studied [17]. The mass hierarchy of neutrinos refers to how the three types of neutrinos are arranged in terms of their relative sizes. There are electron, muon, and tau neutrinos, each with their own antineutrinos. Despite being light compared to other particles, neutrinos do have small but nonzero masses according to current knowledge on neutrino physics. The mass hierarchy of neutrinos can be classified as normal or inverted. In the normal hierarchy, the masses are ordered as $ m_1<m_2<m_3 $, with $ m_1 $ corresponding to the lightest, $ m_2 $ corresponding to the second lightest, and $ m_3 $ corresponding to the heaviest neutrino [18]. Conversely, in the inverted hierarchy, the masses are arranged as $ m_3<m_1<m_2 $. Determining the neutrino mass hierarchy is an important topic in neutrino physics because it affects various astrophysical and cosmological phenomena. Current evidence for the neutrino mass hierarchy comes from the observation of neutrino oscillations, which refer to the ability of neutrinos to change their type as they travel in space owing to quantum mechanical mixing between the three types of neutrinos. The probability of oscillation is influenced by the differences in the squared masses of the three neutrino types and the mixing angles between them. Experiments such as Super-Kamiokande and Daya Bay have provided valuable information about the neutrino mass differences and mixing angles. According to this information, it is highly probable that the neutrino mass hierarchy is normal [19]. However, future experiments such as the Deep Underground Neutrino Experiment (DUNE) will provide more accurate measurements, allowing for a definitive determination of the mass hierarchy.

      The quantum field theory of neutrinos coupled to gravity serves as the theoretical framework for studying neutrino oscillation in curved spacetime [14]. In this framework, the probability of oscillation depends on various factors such as neutrino energy, mass-squared differences, and curvature of spacetime. The metric tensor describes the curvature of spacetime, which is influenced by the gravitational field as well as the distribution of matter and energy. Numerous studies have explored the effects of curved spacetime on neutrino oscillation probability, considering aspects such as the gravitational redshift and curvature-induced potential. These studies have demonstrated that the gravitational field can modify the oscillation probability, leading to potentially observable consequences. Examination of neutrino oscillation in curved spacetime is an active area of research with significant implications for astrophysics and cosmology [17, 20]. The development of theoretical models that describe the quantum field theory of neutrinos coupled to gravity and the impact of curved spacetime on the oscillation probability have been addressed in several studies (see, for example, [21]).

      The rotation of spacetime under weak gravity conditions has been extensively investigated in relation to neutrino oscillations, particularly when neutrinos travel along the equatorial plane. By using the asymptotic form of the Kerr metric, it has been demonstrated that the rotation of a gravitational source significantly changes the phase of neutrinos. Specifically, when neutrinos are generated near a black hole with angular momentum and detected on the same side without the influence of gravitational lensing, the probability of oscillation differs greatly to that observed in the Schwarzschild spacetime [18]. The effects of gravitational lensing on neutrino oscillations within the framework of the γ-spacetime have also been studied employing a quantum-mechanical approach for relativistic neutrinos [22]. This study examined both radial and non-radial propagations, taking into account the phase of neutrino oscillations within this specific spacetime. Additionally, the presence of massive objects in the universe can impact the probability of neutrino oscillation, which has implications in the prediction of the cosmic neutrino background.

    II.   DIRAC EQUATION IN CURVED SPACETIME
    • It is crucial to first understand the fundamental methods used for solving the Dirac equation in curved spacetime before delving into the properties of neutrinos in the presence of gravity. The Dirac equation, which governs the behavior of a massive spinor field on a torsion-free pseudo-Riemannian manifold, can be easily extended as [23, 24]

      $ \begin{array}{*{20}{l}} {\rm i} \hbar \gamma^\mu \ {\cal D}_\mu \ \psi(x) = m c\ \psi(x) \ , \end{array} $

      (1)

      and using the relation $ {\cal D}_\mu = \left(\partial_\mu + \Gamma_\mu\right) $, which is the covariant derivative for a spinor field, Eq. (1) can be expressed as [25, 26]

      $ \begin{array}{*{20}{l}} \left[{\rm i} \hbar \gamma^\mu \ \left(\partial_\mu + \Gamma_\mu\right)- m c\right] \psi(x) = 0 \ . \end{array} $

      (2)

      Here, $ \Gamma_\mu $ represents the spin connection and $ \gamma^\mu $ is associated with the covariant Dirac matrices, which are linked to spacetime through the following relations [27, 28]:

      $ \begin{array}{*{20}{l}} \{ \gamma^\mu, \gamma^\nu \} = 2 g^{\mu\nu} \ . \end{array} $

      (3)

      In this context, the anti-commutation operation is denoted by curly brackets, and spin connection $ \Gamma_\mu $ is determined by the following condition [23]:

      $\frac{\partial \gamma_\nu}{\partial x^\mu} - \Gamma_{\nu\mu}^\lambda \gamma_\lambda - \Gamma_\mu \gamma_\nu + \gamma_\nu \Gamma_\mu = 0. $

      (4)

      Let us now define constant Dirac matrices $ \gamma^{(a)} $ as

      $ \begin{array}{*{20}{l}} \gamma^{(a)} = e_\mu^{(a)} \gamma^\mu. \end{array} $

      (5)

      Here, $ e_\mu^{(a)} $ is the orthogonal tetrad that fulfills the following relationship:

      $ \begin{array}{*{20}{l}} g_{\mu\nu} = e_\mu^{(a)} e_\nu^{(b)} \eta_{ab}. \end{array} $

      (6)

      In the convention of flat metric, where $\eta_{ab} ={\rm diag}(-c^2, 1, 1, 1)$, the expression for the spin connection can be expressed using these constant Dirac matrices [25]:

      $ \Gamma_\mu = \frac{\{\gamma^{(a)},\gamma^{(b)}\}}{8} \ g_{\nu\lambda} e_{(a)}^\nu \ \nabla_\mu e_{(b)}^\lambda. $

      (7)

      The action corresponding to Eq. (1) is

      $ {\cal S} = \int {\rm d}^4 x \ \sqrt{g} \ {\cal L}_D. $

      (8)

      Here, g is defined as $ g=g^{\mu\nu} g_{\mu\nu} $, and we can express the Lagrangian as [29]

      $ {\cal L}_D = \frac{\rm i}{2} \left[\bar{\psi}\gamma^\mu \ {\cal D}_\mu \ \psi - \left({\cal D}_\mu \bar{\psi}\right)\gamma^\mu \psi\right]- m \bar{\psi} \psi. $

      (9)
    III.   WKB APPROXIMATION FOR DIRAC EQUATION
    • We are searching for an approximate solution to Eq. (1) by employing the Wentzel–Kramers–Brillouin (WKB) approximation method [30]. Different forms of the WKB approximation have been proposed. Investigation on complex spinor $ \psi(x) $ revealed that it can be decomposed into two components: amplitude $ \xi = \xi(x) $ and semi-classical phase $ {\cal S} = {\cal S}(x) $ as [31]

      $ \begin{array}{*{20}{l}} \psi (x) = {\rm e}^{-\frac{\rm i}{\hbar} {\cal S}(x)} \xi(x) \ . \end{array} $

      (10)

      The wave function in the presence of spin connection can be expressed as [32, 33, 29]

      $ \psi (x) = {\rm e}^{-\frac{\rm i}{\hbar} {\cal S}(x)} {\rm e}^{-\Gamma_\mu x^\mu} \sum\limits_{n=0}^{\infty} \left(\frac{\hbar}{i}\right)^n \xi_n(x) \ . $

      (11)

      The Dirac matrix product in the spin connection term is expressed as [34, 27]

      $ \begin{array}{*{20}{l}} \gamma^a \{\gamma^b, \gamma^c\} = 2 \eta^{ab} \gamma^c-2 \eta^{ac} \gamma^b-2{\rm i}\epsilon^{dabc}\gamma_5 \gamma_d \end{array}. $

      (12)

      Here, $ \eta^{ab} $ represents the metric in a flat space, while $ \epsilon^{abcd} $ denotes the totally antisymmetric tensor in the same flat space. The spin connection can be expressed in terms of the matrix that violates parity as

      $ \Gamma_\mu =\frac{\gamma_5}{2i} \sqrt{-g} \ {\cal A}_\mu, $

      (13)

      where

      $ {\cal A}_\mu=\frac{\sqrt{-g}}{4} e_a^\mu \epsilon^{abcd} \left(e_{b\nu,\sigma}-e_{b\sigma,\nu}\right)e_c^\nu e_d^\sigma. $

      (14)

      One can infer that the additional phase factor in Eq. (11) effectively confirms the interaction between the metric and the spin orientation of the spinor.

      By substituting Eq. (10) into Eq. (1) and equating terms with equal powers of $ \hbar $, we obtain a series of recursive equations for amplitudes $ \xi_n $:

      $ \begin{array}{*{20}{l}} &\left[\gamma^\nu \partial_\nu {\cal S}(x) + m c\right] \xi_0(x)=0 \ , \end{array} $

      (15)

      $ \begin{array}{*{20}{l}} & \left[\gamma^\nu \partial_\nu {\cal S}(x) + m c\right] \xi_n(x)=\left[\gamma^\nu \partial_\nu {\cal S}(x) + m c \right]\xi_{n-1}(x) \end{array}. $

      (16)

      As a result of multiplying by $ \left[\gamma^\nu \partial_\nu {\cal S}(x) - m \right] $ the left-hand side of Eq. (15), the Hamilton-Jacobi equation for a massive particle in a curved spacetime can be expressed as [35, 36, 37]

      $ \begin{array}{*{20}{l}} g^{\mu\nu} \partial_\mu {\cal S}(x) \partial_\nu {\cal S}(x)-m^2c^2=0. \end{array} $

      (17)

      As long as the four-momentum of the particle is known, we can link this expression to the classical action of a particle with mass m on a torsion-free pseudo-Riemannian manifold, allowing us to associate phase $ {\cal S}(x) $ with this action:

      $ p_\mu = m g_{\mu\nu} \frac{{\rm d}x^\nu}{{\rm d}\tau}. $

      (18)

      If we recognize Eq. (15) and establish this identification, it becomes equivalent to the mass-shell condition,

      $ \begin{array}{*{20}{l}} p_\mu = \partial_\mu {\cal S}(x) \end{array}, $

      (19)

      and the solution of Eq. (15) can be expressed as

      $ \begin{array}{*{20}{l}} {\cal S}(x)= \int p_\mu {\rm d}x^\mu \end{array}. $

      (20)

      The Lagrangian expression, which characterizes the geodesic motion, can be formulated as $ {\cal L}=\sqrt{g_{\mu\nu}\dot{x}^\mu \dot{x}^\nu} $, and the geodesic equation can be expressed as

      $ \begin{array}{*{20}{l}} \ddot{x}^\mu + \Gamma_{\mu \sigma}^\nu \dot{x}^\mu \dot{x}^\sigma = 0 \ . \end{array} $

      (21)
    • A.   Dynamics of neutrino spin in external fields within a curved spacetime

    • The study of spin oscillations of massive Dirac neutrinos in the presence of background matter, electromagnetic fields, and gravitational fields is a complex and ongoing research area in theoretical physics. Neutrinos are fundamental particles with extremely weak interaction with matter, making their study challenging and intriguing. Next, we delve deeper into some of the key aspects of studying spin oscillations of massive Dirac neutrinos in background matter, electromagnetic fields, and gravitational fields:

      ● Matter Effects and Neutrino Oscillations: When neutrinos propagate through a medium, such as the dense matter found in the core of a star or the early universe, their interactions with the medium can modify their oscillation behavior. This is known as the matter effect or MSW effect, named after the physicists who first studied it (Mikheyev, Smirnov, and Wolfenstein).

      This matter effect arises from the presence of charged particles in the medium. Neutrinos can experience forward scattering interactions with these charged particles, leading to an effective potential that depends on the neutrino flavor. As a result, the flavor oscillation probability of neutrinos can be significantly altered with respect to their behavior in a vacuum. The matter effect can induce resonances, which maximally modifies the oscillation probability, leading to interesting phenomena in neutrino oscillation experiments.

      ● Electromagnetic Fields and Spin Precession: Neutrinos, being electrically neutral particles, do not directly interact with electromagnetic fields. However, they possess a magnetic dipole moment, which allows for an indirect interaction with magnetic fields. When neutrinos propagate through regions with magnetic fields, such as in astrophysical environments or laboratory experiments, they can experience spin precession.

      Spin precession refers to the rotation of the neutrino's spin around the direction of the magnetic field. This precession can modify the flavor oscillation probability of neutrinos and introduce new effects that depend on the relative orientation between the neutrino's momentum, magnetic field, and direction of propagation. The study of spin precession in neutrinos requires a careful treatment of the neutrino's magnetic properties and their interactions with magnetic fields.

      ● Gravitational Fields and General Relativity: Neutrinos, like all particles, are influenced by gravitational fields according to the principles of general relativity. In the presence of a gravitational field, the curvature of spacetime affects the propagation of neutrinos. This can lead to modifications in their oscillation behavior and introduce additional complexities.

      The gravitational interaction can cause the trajectory of neutrinos to deviate and induce effects such as gravitational redshift and time dilation. These gravitational effects can impact the neutrino oscillation probability and potentially generate spin oscillations as well. The study of neutrino oscillations in the context of general relativity requires a combination of quantum field theory, general relativity, and development of suitable theoretical frameworks.

      ● Experimental Probes and Future Directions: Experimental efforts play a crucial role in studying the spin oscillations of massive Dirac neutrinos in various physical environments. Neutrino oscillation experiments, conducted at particle accelerators, underground laboratories, or using astrophysical neutrino sources, provide valuable data for testing theoretical predictions and exploring the properties of neutrinos.

      Future experiments, such as the Deep Underground Neutrino Experiment (DUNE) and Jiangmen Underground Neutrino Observatory (JUNO), aim to study neutrino oscillations with higher precision and investigate matter effects, electromagnetic field interactions, and the impact of gravitational fields on neutrino's behavior.

      In [37], neutrino spin oscillations within external fields in curved spacetime were investigated. The contributions of this study have been highly valuable to the field. The researchers analyzed the evolution of neutrino's spin in the presence of background matter and an external electromagnetic field within a curved spacetime. The primary motivation behind this study was to provide evidence supporting the validity of the quasiclassical equation governing neutrino spin evolution. They successfully derived a covariant equation for this purpose, starting from the Dirac equation that describes the interaction between a massive neutrino and external fields in a curved spacetime as

      $ \left[{\rm i} \hbar \gamma^\mu \ {\cal D}_\mu - \frac{\mu}{2} F_{\mu\nu} \sigma^{\mu\nu} - \frac{V^\mu}{2} \gamma_\mu \left(1-\gamma^5\right)\right] \ \psi(x) = m c\ \psi(x) \ , $

      (22)

      where $ F_{\mu\nu} $ is the Faraday tensor defined as

      $ \begin{array}{*{20}{l}} F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu \ , \end{array} $

      (23)

      and $ A_\mu $ is the four-vector potential of the electromagnetic field. In Eq. (15), $ {\cal D}_\mu = \left(\partial_\mu + \Gamma_\mu\right) $ represents the covariant derivative, where $ \Gamma_\mu $ denotes the spin connection. The symbols $ \gamma^\mu = \gamma^\mu (x) $, $\sigma^{\mu\nu}= \dfrac{\rm i}{2}\left(\gamma_\mu \gamma_\nu -\gamma_\nu \gamma_\mu\right)$, and $\gamma^5= - \dfrac{\rm i}{4!} E^{\mu\nu\alpha\beta }\gamma_\mu \gamma_\nu \gamma_\alpha \gamma_\beta$ refer to the coordinate-dependent Dirac matrices. Here, $ E^{\mu\nu\alpha\beta } = \varepsilon^{\mu\nu\alpha\beta }/\sqrt{-g} $ represents the covariant antisymmetric tensor in the context of curved spacetime, where $ g={\rm det}(g_{\mu\nu}) $ and $ g_{\mu\nu} $ denote the metric tensor. Symbol μ represents the magnetic moment of a neutrino. The expression for $ V^\mu = (V^0, \mathbf{V}) $, which serves as the effective potential governing the interaction of a neutrino with arbitrarily polarized and moving matter, can be found in [38].

      The relationship between Eqs. (15) and (1) becomes apparent when we observe that neglecting the impact of background matter and electromagnetic fields results in considering only gravitational effects. Consequently, Eq. (15) can be expressed equivalently to Eq. (1) as

      $ \begin{array}{*{20}{l}} \left[{\rm i} \hbar \gamma^\mu \ \left(\partial_\mu + \Gamma_\mu\right)- m c\right] \psi(x) = 0 \ . \end{array} $

      (24)

      The expression for the Dirac equation in a locally Minkowskian frame was derived in [38]. Recognizing the work of these authors in this area, we consider unnecessary to duplicate their efforts. Instead, we focus on obtaining the Dirac equation that specifically characterizes the interaction between a massive neutrino and external fields within a curved spacetime, tailored to our unique case:

      $ \begin{aligned}[b] & \Bigg[{\rm i} \hbar \gamma^\mu \ \partial_\mu + \frac{ \hbar }{2} \gamma^\mu \gamma_5\sqrt{-g} \ {\cal A}_\mu\\&\quad- \frac{\mu}{2} F_{\mu\nu} \sigma^{\mu\nu} - \frac{V^\mu}{2} \gamma_\mu \left(1-\gamma^5\right) - m c\Bigg] \ \psi(x) = 0 \ , \end{aligned} $

      (25)

      and, according to Ref. [39], the covariant equation governing the quasiclassical evolution of the neutrino spin, denoted as $ S^\mu $ in the presence of general external fields, can be derived. This derivation relies on the Heisenberg equation applied to the corresponding spin operator, taking into account the influence of the external fields. Subsequently, the equation is subjected to an averaging process over the neutrino wave packet. By employing Eq. (25), we obtain the Lorentz invariant expression for the evolution equation of neutrino spin $ S^\mu $, which accounts for the general interactions with external fields as

      $ \frac{{\rm d} S^\mu}{{\rm d} \tau} = 2 \mu \left(F^{\mu\nu} S_\nu - u^\mu F^{\lambda\rho} u_\lambda S_\rho \right) + \sqrt{2} G_F {\cal K}^{\mu\nu} S_\nu + {\cal G}^{\mu\nu} S_\nu, $

      (26)

      where $G_F = 1.17 \times 10^{-5}\; {\rm GeV}^{-2}$ is the Fermi constant. Additionally, the expressions for tensors $ {\cal K}^{\mu\nu} $ and $ {\cal G}^{\mu\nu} $ are

      $ \begin{aligned}[b] &{\cal K}^{\mu\nu} = \varepsilon^{\mu\nu\alpha\beta} V_\alpha u_\beta \ , \\ &{\cal G}^{\mu\nu} = \left( \gamma^{\mu\nu\lambda}+\gamma^{\lambda\mu\nu}+\gamma^{\nu\lambda\mu}\right) u_\lambda \ . \end{aligned} $

      (27)
    IV.   NEUTRINOS
    • The production and detection of neutrinos occur in various flavor eigenstates represented by $ |\nu_\alpha\rangle $. These flavor eigenstates are combinations of mass eigenstates represented by $ |\nu_i\rangle $. Therefore, a flavor eigenstate can be expressed in terms of mass eigenstates as discussed in [16, 40],

      $ |\nu_\alpha\rangle=\sum\limits_{i}U_{\alpha i}^* |\nu_i\rangle\ , $

      (28)

      for a set of three neutrino flavors $ \alpha=\{e, \mu, \tau\} $ and a set of three generations $ i=\{1,2,3\} $. Maki-Nakagawa-Sakata-Pontecorvo (MNSP) matrix U, also known as the neutrino flavor mixing unitary matrix, plays a role similar to the Cabibbo-Kobayashi-Maskawa matrix in governing quark mixing. In the case of three generations of neutrinos, the MNSP matrix is characterized by three mixing angles $ \theta_i $, phase δ that describes CP-violation, and two additional phases $ \alpha_1 $ and $ \alpha_2 $, which can only be non-zero if neutrinos are Majorana particles. If neutrinos are Dirac particles, then $ \alpha_1 = \alpha_2 = 0 $.

      In a process of propagation, the neutrino moves from source S to detector D positioned at $ x_S $ and $ x_D $, respectively. The amplitude associated with the event of detecting a neutrino of flavor α at position $ x_S $ and observing it as a neutrino of flavor β at position $ x_D $ is described by

      $ {\cal A}_{\beta\alpha} = |\langle\nu_\beta(x_D)|\nu_\alpha(x_S)\rangle| = \sum\limits_{i} U_{\alpha i}^* U_{\beta i} |\langle\nu_i(x_D)|\nu_i(x_S)\rangle|. $

      (29)

      To estimate spinor $ \nu_i $, the WKB approximation in Eq. (11) incorporates both action $ {\cal S}_i(x) $ for the i-th mass eigenstate and spin connection $ \Gamma_\mu $. Previous studies [41, 42, 43] have explored three distinct scenarios that give rise to neutrino oscillation: (a) a flat spacetime, (b) a curved spacetime within a non-rotating frame, and (c) a curved spacetime within a rotating frame. In scenarios (a) and (b), the phase difference of the neutrinos relies solely on $ {\cal S}(x) $ and does not involve the spin connection, resulting in the following phase difference [44]:

      $ {\cal S}\left(m_i, x_D-x_S\right) \approx {\cal S}_i(x_D) - {\cal S}_i(x_S) = \int_{x_S}^{x_D} p_\mu {\rm d}x^\mu. $

      (30)

      The presence of parity-violating matrix $ \gamma_5 $ in the representation of the spin connection [14] in Eq. (13) reveals that in scenario (c), there is an additional contribution to the phase shift when there exist differences in spin orientation between the two mixing eigenstates $ \nu_i $. Moving forward, we will exclusively examine massive neutrinos with identical spin orientations, disregarding any contributions stemming from different spin orientations. This decision is based on the fact that our research is primarily centered around neutrino flavor oscillations that occur under the influence of action $ {\cal S}(x) $. Consequently, we will solely consider the same spin orientation for massive neutrinos, dismissing any interaction between the neutrino spin and the metric that arises from $ \Gamma_\mu $.

      The expression for the amplitude of the neutrino flavor transition can be formulated as follows:

      $ {\cal A}_{\beta\alpha} = \sum\limits_{i} U_{\alpha i}^* U_{\beta i} {\rm e}^{-{\rm i} {\cal S} \left(m_i, x_D-x_S\right)}. $

      (31)

      The difference in phase between two mass eigenstates can be expressed in the following form [37]:

      $ \begin{array}{*{20}{l}} \Phi_{ij} = {\cal S} \left(m_i, x_D-x_S\right) - {\cal S} \left(m_j, x_D-x_S\right) \end{array}. $

      (32)

      The probability of a transition in neutrino flavor from the initially produced α flavor to the detection points can be calculated as [15, 45]

      $ {\cal P}_{\beta\alpha}= \big|{\cal A}_{\beta\alpha} \big|^2=\sum\limits_{i,j}U_{\beta i}U_{\beta j}^*U_{\alpha j}U_{\alpha i}^*{\rm e}^{-{\rm i}\Phi_{ij}}\ . $

      (33)

      It is important to note that the action assumes the following form when m is small:

      $ {\cal S} \left(m_i, x_D-x_S\right) = \sum\limits_{n=0}^{\infty} \frac{(m_i^2)^n}{n!} {\cal S}^{(n)} \left( x_D-x_S\right), $

      (34)

      where

      $ {\cal S}^{(n)} \left( x_D-x_S\right) = \frac{\partial^{(n)} {\cal S} \left(m_i, x_D-x_S\right)}{\partial^{(n)} m_i^2}. $

      (35)

      As a result, using the Taylor series, one can express the phase difference as

      $ \Phi_{ij} = \Delta m_{ij}^2 \ {\cal S}^{(1)} \left( x_D-x_S\right) + \frac{\Delta m_{ij}^4}{2} {\cal S}^{(2)} \left( x_D-x_S\right) + ..., $

      (36)

      where

      $ \begin{array}{*{20}{l}} \Delta m_{ij}^2 = m_i^2-m_j^2 \ , \qquad \Delta m_{ij}^4 = \left(m_i^2+m_j^2\right) \Delta m_{ij}^2 \end{array}. $

      (37)
    V.   TWO-FLAVOUR NEUTRINO OSCILLATIONS
    • The result expressed by Eq. (33) is applicable to all numbers of neutrino generations and neutrino energies. However, the likelihood of conversion into a specific neutrino flavor is reduced in certain scenarios, including solar neutrino mixing. In these situations, the MNSP matrix [46, 47] simplifies to a member of the $ S O(2)$ group and can be represented by a single mixing angle, Θ, as follows:

      $ \begin{array}{*{20}{l}} U = \begin{pmatrix} \cos\Theta & \sin\Theta \\ -\sin\Theta& \cos\Theta \end{pmatrix} \end{array}. $

      (38)

      In the case of two flavors of neutrino oscillations, there is a single difference in mass denoted as $ \Delta m_{12}^2 = m_1^2 - m_2^2 $ and another difference $ \Delta m_{12}^4 = \left(m_1^2 + m_2^2\right) \Delta m_{12}^2 $. By simplifying Eq. (33), the final form of the probability of neutrino oscillation can be expressed as [48, 49]

      $ \begin{array}{*{20}{l}} {\cal P}_{\beta\alpha}&= \left\{ \begin{matrix} \sin^2\Theta\sin^2\Phi_{12} \ , & \beta \neq \alpha\\ 1-\sin^2\Theta\sin^2\Phi_{12} \ ,& \beta = \alpha \end{matrix}\right. \end{array}. $

      (39)

      Here, phase shift $ \Phi_{12} $ up to second order can be determined using Eq. (36) as

      $ \Phi_{12} = \Delta m_{12}^2 \ {\cal S}^{(1)} \left( r_D-r_S\right) + \frac{\Delta m_{12}^4}{2} {\cal S}^{(2)} \left( r_D-r_S\right). $

      (40)

      Therefore, the probability expression takes the following form:

      $ {\cal P}_{\beta\alpha}= \left\{ \begin{split}& \sin^2\Theta\sin^2\Big[ {\cal S}^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2 \\& \quad+ \dfrac{1}{2} {\cal S}^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ ,\quad \beta \neq \alpha\\ & 1-\sin^2\Theta\sin^2\Big[ {\cal S}^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2 \\& \quad+ \dfrac{1}{2} {\cal S}^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ ,\quad \beta = \alpha \end{split}\right. $

      (41)

      and here, we can insert the expressions for $ {\cal S}^{(1)} \left(r_D-r_S\right) $ and $ {\cal S}^{(2)} \left(r_D-r_S\right) $ given by the Eqs. (53) and (54), respectively. As a result, the Eq. (41) becomes a general formula to calculate the probability of a neutrino undergoing flavor transition as it propagates from its source to a receiver.

    VI.   JOHANNSEN SPACETIME
    • We next focus on evaluating the Johannsen spacetime, which is a more general version of the Kerr spacetime, and it can be defined using the following metric [50]:

      $ \begin{aligned}[b] {\rm d}s^2=\;&-\frac{\tilde\Sigma\left(\Delta-a^2A_2^2(r)\sin^2\theta\right)}{B^2}{\rm d}t^2+\frac{\tilde\Sigma}{A_5(r)\Delta }{\rm d}r^2 +\frac{\tilde\Sigma}{A_6(\theta)} {\rm d}\theta^2- \\& - \frac{2a\left[\left(r^2+a^2\right)A_1(r)A_2(r)-\Delta\right]\tilde\Sigma\sin^2\theta}{B^2}{\rm d}t{\rm d}\phi + \\& +\frac{\left[\left(r^2+a^2\right)^2A_1^2(r)-a^2\Delta\sin^2\theta\right]\tilde\Sigma\sin^2\theta}{B^2}{\rm d}\phi^2 \ , \end{aligned} $

      (42)

      where

      $ \begin{aligned}[b] &B=A_1(r)A_3(\theta)(r^2+a^2)-A_2(r)A_4(\theta)a^2\sin^2\theta \\ & \tilde\Sigma=\Sigma+f(r)+g(\theta) \ , \quad \Delta=r^2-2Mr+a^2 \ , \\ & \Sigma=r^2+a^2\cos^2\theta. \end{aligned} $

      (43)

      Generally speaking, it is not possible to locate the stationary points of the effective potential of the function within the given background spacetime expressed by Eq. (42). Nevertheless, if we specifically select the profile functions as follows, it may become feasible:

      $ \begin{aligned}[b] A_1(r)=1+\sum\limits_{k=3}^{\infty}\alpha_{1k}\left(\frac{M}{r}\right)^k\ ,\quad A_2(r)=1+\sum\limits_{k=2}^{\infty}\alpha_{3k}\left(\frac{M}{r}\right)^k\ , \\ \end{aligned} $

      $ \begin{aligned}[b] &A_5(r)=1+\sum\limits_{k=2}^{\infty}\alpha_{5k}\left(\frac{M}{r}\right)^k\ ,\quad f(r)=r^2\sum\limits_{k=3}^\infty\epsilon_k\left(\frac{M}{r}\right)^k\ , \\ &A_3(\theta)=A_4(\theta)=A_6(\theta)=1\ , \qquad g(\theta)=0\ , \end{aligned} $

      (44)

      in which the equatorial plane at $ \theta_0=\pi/2 $ is where the stationary points of the function effective potential can be found. It is crucial to highlight that the Johannsen spacetime is distinguished by a set of parameters, namely, $ \alpha_{1k} $, $ \alpha_{3k} $, $ \alpha_{5k} $, and $ \epsilon_k $, which vary as a function of the mass and spin of the black hole. It is also important to emphasize the following key features of this parametrization: (a) the metric retains its smoothness at all points both inside and outside the event horizon, and (b) it has been convincingly demonstrated that certain black hole solutions can be accurately reproduced in alternative theories of gravity by appropriately choosing the deformation parameters [51]. In addition, the spacetime is derived by enforcing the separability of the Hamilton-Jacobi equations, despite lacking a theoretical basis for doing so. It is worth noting that alternative theories of gravity exist in which non-Kerr black hole solutions do not meet this condition. However, maintaining separability may aid in certain calculations. This spacetime model finds application in conducting phenomenological calculations in the field of black hole astrophysics [52, 53].

      The equation for massive neutrinos, known as the Hamilton-Jacobi equation, can be expressed as [54]

      $ -2\frac{\partial {\cal S}}{d\tau} = g^{\mu\nu}\frac{\partial{\cal S}}{\partial x^\mu}\frac{\partial {\cal S}}{\partial x^\nu}. $

      (45)

      Using the metric geometry expressed by Eq. (42), Eq. (45) takes the form

      $ \begin{aligned}[b] -2\frac{\partial {\cal S}}{d\tau} =\;&g^{tt}\left(\frac{\partial{\cal S}_t}{\partial t}\right)^2 + 2g^{t\phi} \frac{\partial{\cal S}_t}{\partial t} \frac{\partial{\cal S}_\phi}{\partial \phi} + g^{\phi\phi}\left(\frac{\partial{\cal S}_\phi}{\partial \phi}\right)^2 \\& +g^{rr}\left(\frac{\partial{\cal S}_r}{\partial r}\right)^2+g^{\theta\theta}\left(\frac{\partial{\cal S}_\theta}{\partial\theta}\right)^2 \end{aligned}. $

      (46)

      The given form of the Hamilton-Jacobi function is [50]

      $ {\cal S}=\frac{1}{2}m^2\tau-Et+L\phi+{\cal S}_r+{\cal S}_\theta. $

      (47)

      After performing certain algebraic calculations, the HJ equation can be separated as follows:

      $ \begin{aligned}[b] \frac{1}{\Delta} \left[(r^2+a^2)A_1(r)E-aA_2(r)L\right]^2 \end{aligned} $

      $ \begin{aligned}[b]&-m^2\left[r^2+f(r)\right]-A_5(r) \Delta \left(\frac{\partial{\cal S}_r}{\partial r}\right)^2 \\ =\;& \left(\frac{\partial{\cal S}_\theta}{\partial\theta}\right)^2 +\left(\frac{A_3(\theta)L}{\sin\theta}-a A_4(\theta)E\sin\theta\right)^2 +m^2a^2\cos^2\theta \ , \end{aligned} $

      (48)

      and the same constant K, which is known as the Carter constant, can be assigned to both sides of the equation:

      $ \begin{aligned}[b] &\left(\frac{\partial{\cal S}_\theta}{\partial\theta}\right)^2 = K -\left(\frac{A_3(\theta)L}{\sin\theta}-a A_4(\theta)E\sin\theta\right)^2 -m^2a^2\cos^2\theta \ , \\ &\left(\frac{\partial{\cal S}_r}{\partial r}\right)^2 = \frac{\left[(r^2 + a^2)A_1(r)E - aA_2(r)L\right]^2 - \Delta \left[K - m^2\left(r^2 + f(r)\right)\right]}{A_5(r) \Delta^2} \end{aligned}. $

      (49)

      It is possible to achieve motion on a plane with fixed angle $ \theta = \theta_0 $ by selecting the indicated angle:

      $ K = \left(\frac{A_3(\theta_0)L}{\sin\theta_0}-a A_4(\theta_0)E\sin\theta_0\right)^2 +m^2a^2\cos^2\theta_0. $

      (50)

      If this is the case, the action is interpreted as

      $ {\cal S}=-Et+L\phi+ \int \frac{\sqrt{R(r)}}{\Delta} {\rm d}r, $

      (51)

      where

      $ R(r) = \frac{\left[(r^2+a^2)A_1(r)E-aA_2(r)L\right]^2-\Delta \left[K-m^2\left(r^2+f(r)\right)\right]}{A_5(r)}. $

      (52)

      As mentioned in Eq. (34), we can consider the action up to the second order in $ m^2 $ and expand it from $ r_S $ (source distance) to $ r_D $ (detector distance) as

      $ {\cal S}^{(1)}\left( r_D-r_S\right) = \int_{r_S}^{r_D} \frac{\left(r^2+f(r)\right) \sqrt{A_5(r)}}{2\sqrt{\left((r^2+a^2)A_1(r)E-aA_2(r)L\right)^2-\Delta K}}{\rm d}r, $

      (53)

      and

      $ {\cal S}^{(2)}\left( r_D-r_S\right) = \int_{r_S}^{r_D} \frac{\left(r^2+f(r)\right)^2 A_5^{3/2}(r)}{4\left[\left((r^2 + a^2)A_1(r)E - aA_2(r)L\right)^2 - \Delta K\right]^{3/2}}{\rm d}r. $

      (54)

      Thus, the Johannsen spacetime reduces to the Kerr spacetime when the parameterizations in Eq. (44) are not taken into account. The expression for $ {\cal S}^{(1)} $ becomes the same as the one obtained in [36], where the first-order result for neutrino oscillations in a Kerr metric was derived. Next, we explore the process of determining the probability for the two-flavor case of neutrino oscillations.

    • A.   Radial propagation

    • The reason why purely radial motion is not possible in the Kerr spacetime is the influence of the black hole's rotation, which introduces a phenomenon known as frame-dragging. This phenomenon is a consequence of the spacetime curvature caused by the rotating black hole. In general relativity, the presence of mass or energy curves the surrounding spacetime, affecting the motion of objects within it. In the case of a rotating black hole described by the Kerr spacetime, the rotation creates a twisting or dragging effect on the nearby spacetime. When a test particle moves radially, aiming directly toward or away from the black hole, it is still subject to the curvature of spacetime caused by the rotating black hole. This curvature, combined with the rotation, causes the test particle to experience an additional angular momentum component. As a result, the trajectory of the particle deviates from a purely radial path, leading to a combination of radial and angular motion. This can be viewed as if the rotating black hole "drags" the nearby spacetime around with it, causing objects to be influenced by the rotation even if they are initially moving radially. This effect prevents purely radial motion in the Kerr spacetime and results in a spiraling or helical trajectory for the test particle. Therefore, in the Kerr spacetime, a test particle cannot move along purely radial trajectories. Owing to the curvature of spacetime caused by the rotating black hole, the particle will experience a combination of radial and angular motion, even if it initially moves radially.

      The Johannsen spacetime, on the other hand, is a modification of the Kerr spacetime that introduces additional parameters to account for potential deviations from the standard Kerr geometry. The specific effects of these deviations on the motion of test particles depend on the particular form of the Johannsen metric and the values of the parameters. Therefore, it is possible that in certain cases or parameter regimes, the Johannsen spacetime may allow for different types of motion, including potentially radial motion. The specific effects of these deviations on the motion of test particles depend on the particular form of the Johannsen metric and the values of the parameters. Without knowing the specific form of the Johannsen metric and the parameter values, it is challenging to provide a definitive answer regarding the existence or nature of purely radial motion in the Johannsen spacetime. The additional parameters introduced in the Johannsen metric could potentially alter the gravitational field in such a way that purely radial motion becomes possible. However, it is also possible that the additional forces or curvature effects introduced by the deviations from the Kerr geometry would still prevent purely radial motion.

      The specific form of the Johannsen metric, which incorporates these additional parameters, determines the exact nature of the deviations from the Kerr spacetime and how they affect the motion of test particles. Depending on the values of these parameters, it is conceivable that the additional forces or curvature effects introduced in the Johannsen spacetime could allow for purely radial motion in certain cases. For example, if the additional parameters in the Johannsen metric introduce modifications to the gravitational field that counteract or weaken the frame-dragging effect caused by the black hole's rotation, purely radial motion could be potentially allowed. This could occur if the additional parameters change the geometry of spacetime in such a way that the angular momentum component induced by the rotation becomes negligible or is canceled out. The characterization of $ \Phi^m_{ij} $ will be influenced by the geodesic parameters, including energy and angular momentum. Next, we consider the scenario of radial propagation with $ L = K = 0 $, where

      $ {\cal S}=-Et+ \int \frac{\sqrt{R_*(r)}}{\Delta} {\rm d}r, $

      (55)

      and

      $ R_*(r) = \frac{(r^2+a^2)^2 A_1^2(r)E^2+ m^2 \Delta \left(r^2+f(r)\right)}{A_5(r)}. $

      (56)

      By expanding the action to its second order in $ m^2 $, we obtain

      $ {\cal S}_*^{(1)}\left( r_D-r_S\right) = \frac{1}{2 E} \int_{r_S}^{r_D} \frac{\left(r^2+f(r)\right) \sqrt{A_5(r)}}{ A_1(r)(r^2+a^2)}{\rm d}r, $

      (57)

      and

      $ {\cal S}_*^{(2)}\left( r_D-r_S\right) = \frac{1}{4 E^3} \int_{r_S}^{r_D} \frac{\Delta \left(r^2+f(r)\right)^2 A_5^{3/2}(r)}{A^3_1(r)(r^2+a^2)^{3}}{\rm d}r. $

      (58)

      From Eq. (41), the probability expression for radial propagation adopts the following form

      $ {\cal P}_{\beta\alpha}= \left\{ \begin{split}& \sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2 \\ &\quad+ \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ , \quad \beta \neq \alpha\\ &1-\sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2 \\&\quad+ \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ ,\quad \beta = \alpha \end{split}\right. $

      (59)

      Consequently, from Eq. (59), it becomes evident that in the case of radial propagation with $ L = K = 0 $, the contribution of ϕ to the phase of oscillations is not discernible. This implies that the phase of oscillations remains independent of ϕ. However, it is important to stress that the specific implications of the Johannsen metric on test particle dynamics and the conditions under which purely radial motion might be allowed depend on the precise form of the metric and the values of the additional parameters. The exploration of these effects and conditions requires detailed analysis and numerical simulations specific to the Johannsen metric [52]. Numerical simulations are used to study the behavior of the Johannsen metric and its effect on various physical phenomena. By simulating the motion of test particles, electromagnetic fields, or other relevant quantities in the Johannsen spacetime, the results from observations or theoretical predictions can be compared [53]. These simulations can help constrain the values of the additional parameters that best match the observed data or desired physical behavior.

      Subsequently, conducting a comparative analysis of the outcomes obtained from the Schwarzschild, Kerr, and Johannsen spacetimes is of great interest. The Johannsen metric is determined by the mass (M) and spin (a) of the black hole, along with four independent functions that account for potential deviations from the Kerr solution. When $ \alpha_{1k} = \alpha_{3k} = \alpha_{5k} = \epsilon_k = 0 $, the metric simplifies to the Kerr solution [50]. In this paper, for the sake of simplicity, we specifically examine two scenarios: one where only $ \alpha_{52} $ is non-zero, and another where only $ \epsilon_3 $ is set to zero. Consequently, we can derive the expressions for $ {\cal S}_*^{(1)}\left(r_D-r_S\right) $ and $ {\cal S}_*^{(2)}\left(r_D-r_S\right) $ as follows:

      $ \begin{aligned}[b] \\[-12pt]{\cal S}_*^{(1)}\left( r_D-r_S\right) = &\frac{1}{2 E} \cdot \frac{1}{ a^3 \cdot M^2 \cdot r_D^2} \cdot \Bigg[a \cdot M^2 \cdot \sqrt{\alpha_{52} \cdot M^2+r_D^2} \cdot \left(a^2 \cdot r_D-\epsilon_3 \cdot M^3\right) + \ \epsilon_3 \cdot M^5 \cdot r_D \cdot \sqrt{\alpha_{52} M^2-a^2} \cdot \cot ^{-1}\\&\left(\frac{a \cdot \sqrt{\alpha_{52} \cdot M^2-a^2}}{a^2-r_D \cdot \sqrt{\alpha_{52} \cdot M^2+r_D^2}+r_D^2}\right) \Bigg] -\frac{1}{2 E} \cdot \frac{1}{ a^3 \cdot M^2 \cdot r_S^2} \cdot \Bigg[a \cdot M^2 \cdot \sqrt{\alpha_{52} \cdot M^2+r_S^2} \cdot \left(a^2 \cdot r_S-\epsilon_3 \cdot M^3\right) \\& + \ \epsilon_3 \cdot M^5 \cdot r_S \cdot \sqrt{\alpha_{52} M^2-a^2} \cdot \cot ^{-1}\left(\frac{a \cdot \sqrt{\alpha_{52} \cdot M^2-a^2}}{a^2-r_S \cdot \sqrt{\alpha_{52} \cdot M^2+r_S^2}+r_S^2}\right) \Bigg] \ . \end{aligned} $

      (60)

      and

      $ \begin{aligned}[b] {\cal S}_*^{(2)}\left( r_D-r_S\right) =\; &\frac{1}{4 E^3} \cdot \frac{r_D^2 \sqrt{\left(\alpha_{52} \cdot M^2+r_D^2\right)}}{8 \left(\epsilon_3 \cdot M^3+r_D^3\right)^2} \cdot \left(\frac{\epsilon_3 \cdot M^3}{r_D}+r_D^2\right)^2 \cdot \Bigg[8 +\frac{12 \left( a^{2} +3 M r_D - \epsilon_3 M^2 ) \right)}{3 \cdot \left(a^2+r_D^2\right) \cdot} \\& -\frac{ M^3 \cdot (\alpha_{52} \cdot r_D+2 \epsilon_3 \cdot (M-2 r_D))}{3 \cdot a^2 \cdot \left(a^2+r_D^2\right) } -\frac{2 \cdot \epsilon_3 \cdot M^5 \cdot (2 \alpha_{52} \cdot (2 M+r_D)+\epsilon_3 \cdot M) }{3 \cdot a^4 \cdot \left(a^2+r_D^2\right) } \\& +\frac{2 a^2 \epsilon_3^2 M^7 (\alpha_{52} M+3 r_D)-11 \alpha_{52} \epsilon_3^2 M^9 r_D}{3 a^8 \left(a^2+r_D^2\right)} -\frac{16 M \log \left(\sqrt{\alpha_{52} M^2+r_D^2}+r_D\right)}{\left(\alpha_{52} M^2+r_D^2\right)^{1/2}} +\frac{1}{a^9 \cdot \sqrt{a^2-\alpha_{52} \cdot M^2} \cdot \left(\alpha_{52} \cdot M^2+r_D^2\right)^{1/2}} \\& +\frac{\epsilon_3 \cdot M^5 \cdot \log (r_D) \cdot \left(3 \epsilon_3 a^4 -8 \alpha_{52} a^2 M^2 (4 \alpha_{52}+3 \epsilon_3)+24 \alpha_{52}^2 \epsilon_3 M^4\right)}{a^8 \cdot \sqrt{\alpha_{52}} \cdot \left(a_{52} \alpha_{52}M^2+r_D^2\right)^{1/2}} \\& +\frac{\epsilon_3 \cdot M^5 \cdot \left( 8 a^2 \alpha_{52} M^2 (4 \alpha_{52}+3 \epsilon_3) -3 \epsilon_3 a^4 -24 \alpha_{52}^2 \epsilon_3 M^4\right) \cdot \log \left(\sqrt{\alpha_{52}} \cdot \sqrt{\alpha_{52} M^2+r_D^2}+a_{52} \cdot M\right)}{a^8 \sqrt{\alpha_{52}} \left(\alpha_{52} M^2+r_D^2\right)^{1/2}} \\& -\frac{12 M \left(a^2-\alpha_{52} M^2\right) \left(a^6 r_D-2 a^4 \epsilon_3 M^3-\epsilon_3^2 M^6 r_D\right) }{\left(\alpha_{52} M^2+r_D^2\right)^{1/2}} -\frac{16 \epsilon_3 M^5 \left(3 a^4 \alpha_{52}-4 a^2 \epsilon_3 M^2+9 \alpha_{52} \epsilon_3 M^4\right)}{a^8 r_D \left(\alpha_{52} M^2+r_D^2\right)^{1/2}} \\& +\frac{3 \epsilon_3^2 M^6 \left(8 \alpha_{52} M^2-5 a^2\right)}{a^6 r_D^2 \left(\alpha_{52} M^2+r_D^2\right)^{1/2}} -\frac{6 \alpha_{52} \epsilon_3^2 M^8}{a^4 r_D^4 \left(\alpha_{52} M^2+r_D^2\right)^{1/2}} +\frac{16 \alpha_{52} \epsilon_3^2 M^9}{a^6 r_D^3 \left(\alpha_{52} M^2+r_D^2\right)^{1/2}} \Bigg] \\& -\frac{1}{4 E^3} \cdot \frac{r_S^2 \sqrt{\left(\alpha_{52} \cdot M^2+r_S^2\right)}}{8 \left(\epsilon_3 \cdot M^3+r_S^3\right)^2} \cdot \cdot \left(\frac{\epsilon_3 \cdot M^3}{r_S}+r_S^2\right)^2 \cdot \Bigg[8 +\frac{12 \left( a^{2} +3 M r_S - \epsilon_3 M^2 ) \right)}{3 \cdot \left(a^2+r_S^2\right) \cdot} \\&+\frac{2 a^2 \epsilon_3^2 M^7 (\alpha_{52} M+3 r_S)-11 \alpha_{52} \epsilon_3^2 M^9 r_S}{3 a^8 \left(a^2+r_S^2\right)} -\frac{ M^3 \cdot (\alpha_{52} \cdot r_S+2 \epsilon_3 \cdot (M-2 r_S))}{3 \cdot a^2 \cdot \left(a^2+r_S^2\right) } \\&-\frac{2 \cdot \epsilon_3 \cdot M^5 \cdot (2 \alpha_{52} \cdot (2 M+r_S)+\epsilon_3 \cdot M) }{3 \cdot a^4 \cdot \left(a^2+r_S^2\right) } -\frac{16 M \log \left(\sqrt{\alpha_{52} M^2+r_S^2}+r_S\right)}{\left(\alpha_{52} M^2+r_S^2\right)^{1/2}} + \frac{1}{a^9 \cdot \sqrt{a^2-\alpha_{52} \cdot M^2} \cdot \left(\alpha_{52} \cdot M^2+r_S^2\right)^{1/2}} \\&+\frac{\epsilon_3 \cdot M^5 \cdot \log (r_S) \cdot \left(3 \epsilon_3 a^4 -8 \alpha_{52} a^2 M^2 (4 \alpha_{52}+3 \epsilon_3)+24 \alpha_{52}^2 \epsilon_3 M^4\right)}{a^8 \cdot \sqrt{\alpha_{52}} \cdot \left(\alpha_{52} M^2+r_S^2\right)^{1/2}} \\& +\frac{\epsilon_3 \cdot M^5 \cdot \left( 8 a^2 \alpha_{52} M^2 (4 \alpha_{52}+3 \epsilon_3) -3 \epsilon_3 a^4 -24 \alpha_{52}^2 \epsilon_3 M^4\right) \cdot \log \left(\sqrt{\alpha_{52}} \cdot \sqrt{\alpha_{52} M^2+r_S^2}+\alpha_{52} \cdot M\right)}{a^8 \sqrt{\alpha_{52}} \left(\alpha_{52} M^2+r_S^2\right)^{1/2}} \\& -\frac{12 M \left(a^2-\alpha_{52} M^2\right) \left(a^6 r_S-2 a^4 \epsilon_3 M^3-\epsilon_3^2 M^6 r_S\right) }{\left(\alpha_{52} M^2+r_S^2\right)^{1/2}} -\frac{16 \epsilon_3 M^5 \left(3 a^4 a_{52}-4 a^2 \epsilon_3 M^2+9 \alpha_{52} \epsilon_3 M^4\right)}{a^8 r_S \left(\alpha_{52} M^2+r_S^2\right)^{1/2}} \\& +\frac{3 \epsilon_3^2 M^6 \left(8 \alpha_{52} M^2-5 a^2\right)}{a^6 r^2 \left(\alpha_{52} M^2+r_S^2\right)^{1/2}} -\frac{6 \alpha_{52} \epsilon_3^2 M^8}{a^4 r_S^4 \left(\alpha_{52} M^2+r_S^2\right)^{1/2}} +\frac{16 \alpha_{52} \epsilon_3^2 M^9}{a^6 r_S^3 \left(\alpha_{52} M^2+r_S^2\right)^{1/2}} \Bigg]. \end{aligned} $

      (61)

      As a result, the probability expression for radial propagation in Johannsen spacetime adopts the following form:

      $ {\cal P}_{\beta\alpha}^{\rm Johannsen}= \left\{ \begin{split}& \sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2\\&\quad + \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ , \quad \beta \neq \alpha\\& 1-\sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2\\&\quad + \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ ,\quad \beta = \alpha \end{split}\right. $

      (62)

      When $ a_{52} = \epsilon_3 = 0 $, the metric expressed by Eq. (42) is simplified and reduces to the Kerr solution. Consequently, the expressions for $ {\cal S}*^{(1)}\left(r_D-r_S\right) $ and $ {\cal S}*^{(2)}\left(r_D-r_S\right) $ also reduce in a manner similar to that observed in Kerr spacetime [36, 18, 29]:

      $ {\cal S}_*^{(1)}\left( r_D-r_S\right) = \frac{1}{2 E} \cdot \left[r_D -r_S-a\left( \arctan\frac{r_D}{a} -\arctan\frac{r_S}{a}\right)\right] \ . $

      (63)

      and

      $ \begin{aligned}[b] {\cal S}_*^{(2)}\left( r_D-r_S\right) = \frac{1}{8 E^3} \cdot \Bigg[ 2 \left(r_D-r_S\right)-3 a\left( \arctan\frac{r_D}{a} -\arctan\frac{r_S}{a}\right) \end{aligned} $

      $ \begin{aligned}[b]&+\frac{a^2 (r_D-4 M)}{a^2+r_D^2}-\frac{a^2 (r_S-4 M)}{a^2+r_S^2} \\& -2 M \ln \left(\frac{a^2+r_D^2}{a^2+r_S^2}\right)+\frac{a^4 M}{\left(a^2+r_D^2\right)^2}-\frac{a^4 M}{\left(a^2+r_S^2\right)^2}\Bigg] \ . \end{aligned} $

      (64)

      As a result, the probability expression for radial propagation in Kerr spacetime adopts the following form:

      $ {\cal P}_{\beta\alpha}^{Kerr}= \left\{ \begin{split}& \sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2\\&\quad + \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ , \quad \beta \neq \alpha\\ &1-\sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2 \\&\quad+ \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ ,\quad \beta = \alpha \end{split}\right. $

      (65)

      If we set $ a_{52} $, $ \epsilon_3 $, and a to zero, the metric given by Eq. (42) is simplified and we obtain the Schwarzschild solution. Consequently, the expressions for $ {\cal S}_*^{(1)}\left(r_D-r_S\right) $ and $ {\cal S}_*^{(2)}\left(r_D-r_S\right) $ also simplify, yielding the same result as in Schwarzschild spacetime [16, 11]:

      $ {\cal S}_*^{(1)}\left( r_D-r_S\right) = \frac{1}{2 E} \cdot \left(r_D -r_S\right) \ . $

      (66)

      and

      $ {\cal S}_*^{(2)}\left( r_D-r_S\right) = \frac{1}{4 E^3} \cdot \Bigg[ r_D-r_S -2 M \cdot \ln \frac{r_D}{r_S} \Bigg] \ . $

      (67)

      As a result, the probability expression for radial propagation in Schwarzschild spacetime adopts the following form:

      $ {\cal P}_{\beta\alpha}^{\rm Schwrzd}= \left\{ \begin{split} &\sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2 \\&\quad+ \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ , \quad \beta \neq \alpha\\& 1-\sin^2\Theta\sin^2\Big[ {\cal S}_*^{(1)} \left( r_D-r_S\right) \Delta m_{12}^2\\&\quad + \frac{1}{2} {\cal S}_*^{(2)} \left( r_D-r_S\right)\Delta m_{12}^4\Big] \ ,\quad \beta = \alpha \end{split}\right. $

      (68)
    VII.   PROPER DICTANCE
    • The neutrino propagates across its proper distance, while $ {\rm d}r $ is simply a coordinate. The proper distance can be expressed as [14]

      $ {\rm d}L_p=\sqrt{\left(\frac{g_{0\nu}g_{0\mu}}{g_{00}}-g_{\nu\mu}\right){\rm d}x^{\nu}{\rm d}x^{\mu}}. $

      (69)

      In the context of Johannsen spacetime, there exists

      $ {\rm d}L_p=\sqrt{-g_{rr} {\rm d}r^2+\left(\frac{g_{t\phi}^2}{g_{tt}}-g_{\phi\phi}\right){\rm d}\phi^2}. $

      (70)

      The following expression is obtained by multiplying by ${\rm d}r^2$

      $ L_p=\int_{r_S}^{r_D}\sqrt{-g_{rr} +\left(\frac{g_{t\phi}^2}{g_{tt}}-g_{\phi\phi}\right)\frac{\dot\phi^2}{\dot r^2}}{\rm d}r. $

      (71)

      With the use of four-velocity normalization, $ g_{\mu\nu}{\dot x}^\mu{\dot x}^\nu=-1 $, and considering the aforementioned expressions, $ \dot{r} $ and $ \dot{\phi} $ can be expressed as

      $ \begin{aligned}[b] \dot{r} =\;& \frac{A_5(r)}{m \tilde{\Sigma}}\sqrt{R(r)} = \frac{\sqrt{A_5(r)}}{m \tilde{\Sigma}}\\& \times \sqrt{\left[(r^2+a^2)A_1(r)E-aA_2(r)L\right]^2-\Delta \left[K-m^2\left(r^2+f(r)\right)\right]} \ , \\ \dot{\phi} =\;& \frac{L}{\tilde\Sigma\sin^2\theta}\frac{\left[A_1(r)(r^2+a^2)-A_2(r)a^2\sin^2\theta \right]^2 }{\left[\left(r^2+a^2\right)^2A_1^2(r)-a^2\Delta\sin^2\theta\right]} \end{aligned}. $

      (72)

      It is apparent that performing an analytical calculation of Eq. (71) is more complex.

    VIII.   FINDINGS AND PROSPECTS FOR THE FUTURE
    • In this study, we delved into the mathematical aspects of the Dirac equation within a curved spacetime and investigated its application in analyzing neutrino oscillations. To achieve this objective, we used the WKB approximation. In particular, we devised a technique for determining the phase shift in flavor neutrino oscillations by employing a Taylor series expansion of the action, considering contributions up to fourth order in $ \Delta m^4 $. In Sec. IV, we examined the intricate dynamics of transition probabilities within our framework, revealing their intricate nature despite fluctuations in mass representation. Furthermore, this method has been employed to evaluate the variation in the phase difference of neutrino mass eigenstates caused by the gravitational field described by the Johannsen spacetime.

      It is well known that the phenomenon of neutrino oscillation, which is a quantum phenomenon, takes place in both flat and curved spacetime metrics. It is anticipated that the presence of massive objects such as stars and black holes, with their gravitational fields, can impact the propagation of neutrinos and alter their oscillation patterns. Research on neutrino oscillation in curved spacetime is being actively pursued, and the resulting findings hold significant implications for the fields of astrophysics and cosmology. Furthermore, it is widely acknowledged that the application of gravitational lensing methods offers compelling proof of the presence of Dark Matter. According to the principles of general relativity, the trajectory of the light can be bent when encountering massive objects or gravitational fields. This deflection is closely linked to the mass of the object and can be likened to the focusing effect of a lens. In upcoming research, we will investigate the gravitational lensing effect on neutrino oscillations to provide evidence of the existence of dark matter.

      We would like to emphasize the significant progress made in a previous study[38] regarding the intricacies of neutrino spin oscillations within a curved spacetime, particularly in the presence of background matter and an external electromagnetic field. The authors successfully derived the Dirac equation that governs neutrino oscillations under these external field conditions. They also conducted numerical solutions to explore the behavior of these oscillations. Numerical simulations and analysis of experimental data are crucial in elucidating the behavior of neutrinos in realistic scenarios. Neutrino oscillation experiments, such as those conducted at particle accelerators or involving atmospheric and solar neutrinos, provide invaluable data for testing theoretical predictions and gaining insights into the spin oscillations of massive Dirac neutrinos in diverse physical environments. This research requires a significant amount of time and extensive numerical calculations owing to the complexity of studying spin oscillations of massive Dirac neutrinos in the presence of background matter, electromagnetic fields, and gravitational fields. In future research, our objective is to explore and uncover the potential for spin precession in the presence of background matter and an external electromagnetic field within a curved spacetime. The derivation of the neutrino spin evolution equation provided in this context relies on the general spin evolution equation in the Heisenberg representation. Adopting this approach will allow for careful analysis of the contributions of various external fields previously mentioned to the evolution of neutrino spin.

Reference (54)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return