Spin oscillations of neutrinos scattered off a rotating black hole

Spin oscillations of neutrinos, gravitationally scattered off a black hole (BH), are studied. The cases of nonrotating and rotating BHs are analyzed. We derive the analytic expressions for the transition and survival probabilities of spin oscillations when neutrinos interact with these gravitational backgrounds. The obtained transition probabilities depend on the impact parameter, as well as the neutrino energy and the particle mass. We find that there is a possibility of spin oscillations of ultrarelativistic neutrinos scattering off a rotating BH. Then, considering the neutrino scattering off BH surrounded by background matter, we derive the effective Schrodinger equation for spin oscillations. The numerical solution of this equation is obtained in the case of a supermassive BH with a realistic accretion disk. Spin effects turn out to be negligible in the neutrino scattering in the Schwarzschild metric. In the Kerr metric, we find that the observed neutrino fluxes can be reduced almost 10% because of spin oscillations when ultrarelativistic neutrinos experience gravitational scattering. The neutrino interaction with an accretion disk results in the additional modification of the intensities of outgoing neutrino fluxes. We consider the applications of the obtained results for the neutrino astronomy.


Introduction
The recent successful studies of oscillations of accelerator [1] and atmospheric [2] neutrinos demonstrate that the masses of these particles are nonzero and there is a mixing between different neutrino flavors. This fact is the direct indication to the physics beyond the standard model. One can study transitions between neutrinos belonging to different flavors, or neutrino flavor oscillations. This type of oscillations is the most plausible solution of the solar neutrino problem [3].
However, in the present work, we shall concentrate on neutrino spin oscillations, which were first proposed in Ref. [4]. Neutrinos are left polarized particles in the standard model. If such a neutrino changes its polarization under the action of some external backgrounds, it cannot be observed since right neutrinos are sterile particles. This process leads to the suppresion of the emitted flux of left neutrinos. The study of neutrino spin and spin-flavor oscillations in various external fields can provide a valuable information about neutrino magnetic moments [5].
The neutrino interactions with external fields, e.g., the electroweak interaction with background matter [6], are known to influence the process of neutrino oscillations. The neutrino interaction with gravitational fields of astrophysical objects, despite of its weakness, can also change the dynamics of neutrino oscillations. Earlier, we examined quasiclassically neutrino spin oscillations under the influence of external fields in curved spacetime in Refs. [7][8][9], where both static metrics and time dependent backgrounds, like a gravitational wave, were studied. The quantum description of the fermion spin evolution in curved spacetime was developed in Ref. [10].
We should make a remark on the studies of neutrino spin oscillations in gravitational fields in Refs. [7,8,11,12]. In these works, neutrinos were supposed to be gravitationally captured by a massive object, like a black hole (BH). In this case, the effect of spin oscillations is not observable since a neutrino detector is typically outside the BH region. Thus one has to consider neutrino spin oscillations in the neutrino gravitational scattering, when one can potentially measure the helicities of both incoming and outgoing particles. Note that flavor oscillations of neutrinos gravitationally scattered off BH were studied in Refs. [13,14].
Before the recent observation of the event horizon silhouette of a supermassive BH (SMBH) [15], which is the main motivation for the present work, all the experimental manifestations of the general relativity, including the direct detection of gravitational waves [16], corresponded to the weak field limit. The bright halo around this silhouette is formed by photons emitted by the accretion disk around SMBH. The size of the event horizon silhouette, computed in Ref. [17], turns out to be in the agreement with the prediction of the general relativity.
However, besides photons, a significant flux of neutrinos is expected to be emitted by dense and hot matter of an accretion disk. These neutrinos were found in Ref. [18] to modify the r-process nucleosynthesis in the vicinity of BH surrounded by an accretion disk. The spin of emitted neutrinos can precess in the gravitational field of a central BH. In this work, we shall examine how a strong gravitational field of BH and the neutrino interaction with an accretion disk can cause their helicity change.
The neutrino gravitational scattering was studied recently [19], in order to compute the size and the shape of the BH shadow formed by these particles [20]. In our work, we shall study neutrino spin oscillations in the gravitational scattering. This process is expected to suppress the flux of particles measured with a terrestrial neutrino telescope.
It was found in Ref. [21] that photons, which form a bright halo around the event horizon silhouette of BH, interact with both its gravitational field and plasma which surrounds BH. This interaction with plasma leads to the modification of the size and the shape of the BH shadow. In the present work, the role of the neutrino interaction with background matter, e.g., with an accretion disk, for the detected flux of gravitationally scattered neutrinos is examined.
In this our work, we continue our studies of neutrino spin oscillations in Refs. [7][8][9]. We start in Sec. 2 with the analysis of the neutrino spin evolution when a particle gravitationally scatters off BH. We find the general expressions for the transition and survival probabilities. Then, we apply our results for the description of spin oscillations in the neutrino scattering off a nonrotating BH in Sec. 2.1 and a rotating one in Sec. 2.2. In Sec. 3, we formulate the effective Schrödinger equation for neutrino spin oscillations in the scattering off BH surrounded by background matter. We study astrophysical applications in Sec. 4. In particular, we consider the effect of spin oscillations on the measured neutrino fluxes when particles scatter off SMBH with a realistic accretion disk. The situations of nonrotating and rotating BHs are studied. Finally, in Sec. 5, we discuss our results. We remind how a scalar particle moves in the Schwarzschild metric in Appendix A and in Kerr metric in Appendix B.

Neutrino spin evolution in gravitational scattering
In this section, we study how the spin of a neutrino evolves when a particle scatters off BH. First, we derive the general expressions for the transition and survival probabilities. Then, we apply this result for spin oscillations in the Schwarzschild and Kerr metrics. In case of the Kerr metric, we study the neutrino motion in the equatorial plane only. In these situations, one can analytically solve the spin evolution equation and obtain the probabilities for spin oscillations in quadratures.
In Refs. [7,8], we found that the neutrino invariant spin ζ, defined in a locally Minkowskian frame, evolves in an external gravitational field as where t is the time in world coordinates and Ω g is the vector accounting for the gravity contribution. If a neutrino interacts with a Schwarzschild or Kerr BH, Ω g in Eq.
(1) has only one nonzero component [7,8], Ω g " p0, Ω 2 , 0q. We are interested in neutrino spin oscillations, i.e. in the change of the neutrino helicity, h " pζuq{|u|, where u is the spatial part of the neutrino four velocity in the locally Minkowskian frame. Therefore, besides the study of the neutrino spin in Eq. (1), we should account for the evolution of u.
In principle, we can avoid simultaneous tracking of both ζ and u. However, we should fix the initial and final neutrino polarizations. As we shall see in Secs. 2.1 and 2.2, at r Ñ 8, upt Ñ˘8q " u˘8 "´˘"E 2´m2 ‰ 1{2 {m, 0, 0¯, where r is the distance between the BH center and a neutrino, E is the neutrino energy, which is the integral of motion in the considered metrics, and m is the neutrino mass. Thus the asymptotic neutrino motion happens along the first axis in the locally Minkowskian frame. In this frame, an incoming neutrino (t Ñ´8) propagates oppositely the first axis. An outgoing particle (t Ñ`8) moves along this axis.
Since only Ω 2 ‰ 0, the nonzero neutrino spin components are ζ 1,3 ‰ 0, and ζ 2 " 0. It is convenient to represent where α is the angle between ζ and the positive direction of the first axis in the locally Minkowskian frame. Now we have to specify the initial condition for Eq. (1). We suppose that, initially, at r Ñ 8, an incoming neutrino is left polarized, i.e. the helicity is negative, h´8 " pζ´8u´8q{|u´8| "´1. Accounting for the expression for u´8 above, we get that ζ´8 1 " 1 and ζ´8 3 " 0, or α´8 " 0 in Eq. (2). The helicity of an outgoing neutrino has the form, h`8 " pζ`8u`8q{|u`8|, where ζ`8 " pcos α`8, 0, sin α`8q and u`8 is given above. Using Eq. (2), we get that h`8 " cos α`8. The transition P LR and survival P LL probabilities for neutrino spin oscillations are where the upper sign stays for P LR and the lower one for P LL .

Schwarzschild metric
First, we study the neutrino motion in the field of a nonrotating BH. Using the spherical coordinates pr, θ, φq, the interval in this case has the form [22, p. 284], where r g " 2M is the gravitational radius, and M is the BH mass. Since the Schwarzschild metric in Eq. (4) is spherically symmetric, we can take that a neutrino moves in the equatorial plane with θ " π{2, i.e. dθ " 0. The nonzero component of Ω g in the Schwarzschild metric has the form [7], where L is the conserved angular momentum of a neutrino. In Eq. (5), U t " dt{dτ " Ep1´r g {rq´1{m is the component of the four velocity U µ in world coordinates.
The expression for u in the Schwarzschild metric was also obtained in Ref. [7], where the signs˘stay for outgoing and incoming neutrinos respectively [see Eq. (33)]. Using Eq. (6) at r Ñ 8, we obtain the expression for u˘8 proposed above. Then we use Eqs. (1), (5), and (33) to determine the evolution of α in Eq. (2). It obeys the equation, where the signs˘stay for outgoing and incoming neutrinos. Now, accounting for the initial condition α´8 " 0 and the fact that α`8 is twice the angle corresponding to the minimal distance between a neutrino and BH, we get that α`8 reads, where y " b{r g , b " L{E a 1´γ´2 is the impact parameter, γ " E{m in the Lorentz factor at the infinity, and x m is the maximal root of the equation Note that y ą y 0 , where for a neutrino not to fall to BH (see Appendix A). For ultrarelativistic neutrinos with γ " 1, we obtain the well known result, y 0 " 3 ? 3{2.
3{2 ď y ă 8. Hence, using Eq. (3), we get that P LR " p1`cos α`8q{2 " 0 at E " m for the arrbirtary impact parameter. It means that there is no spin flip of ultrarelativistic neutrinos when they scatter off a Schwarzschild BH. This fact is in agreement with the results of Ref. [23]. Now we can correct the result of Ref. [24], where the nonzero transition probability for spin oscillations of ultrarelativistic neutrinos scattering off a Schwarzschild BH was obtained. That incorrect result is a consequence of the extra factor 1{U t in Eq. (1) used in Ref. [24]. The vector Ω g in Eq. (1) already accounts for the change of the proper time τ to the world time t.
At the end of this section, we return to the case of the finite neutrino energy. Considering the situation y " y 0 , the roots of Eq. (9) have the form, When these x 0,1,2 in Eq. (12) are used in Eq. (8), the value of α`8 turns out to be infinite. Thus one cannot express α`8 as a power series of γ´1 at y " y 0 . It means that the neutrino spin makes infinite revolutions with respect to the neutrino velocity when a particle asymptotically approaches a nonrotating BH. In Fig. 1, we show the transition probabilities of spin oscillations of massive neutrinos scattered off a nonrotating BH versus the impact parameter. The maximal transition probability is reached when y Ñ y 0 . One can see that the greater γ " E{m is the smaller P LR is. This behavior of P LR confirms our result that there are no spin oscillations in scattering of massless neutrinos in the Schwarzschild metric.

Kerr metric
Now we study spin evolution of neutrinos scattered off a rotating BH. The space-time in this case is described by the Kerr metric. In Boyer-Lindquist coordinates pt, r, θ, φq, this metric has the form [25], where ∆ " rpr´r g q`a 2 , Σ " r 2`a2 cos 2 θ, Ξ "`r 2`a2˘Σ`r r g a 2 sin 2 θ.
The parameter a in Eqs. (13) and (14) can be in the range, 0 ă a ă r g {2. The angular momentum of BH is J " M a.  Figure 1: The transition probability P LR of spin oscillations versus the dimensioless impact parameter y for neutrinos scattering off a Schwarzschild BH for different neutrino energies.
We study a neutrino moving in the equatorial plane of a Kerr BH with θ " π{2 and dθ " 0. As in the Schwarzschild case, we also have two integrals of motion: the energy E and the neutrino angular momentum L. However, we should consider two cases L ą 0 and L ă 0 in the Kerr metric.
The vector Ω g has the nonzero component [8], where b 2 "´a rpr´r g q`a 2 2r 2 a rrr 3`a2 pr`r g qs rp2r 3´r g a 2 qU φ`r g aU t s, e 1 " r g 2r 2 rap3r 2`a2 qU φ´p r 2`a2 qU t s a rrr 3`a2 pr`r g qs , are the components the gravi-magnetic and gravi-electric fields, rpr´r g q`a 2 , u 3 " pr 3`r a 2`r g a 2 qU φ´a r g U t a rrr 3`a2 pr`r g qs , are the components of the neutrino four velocity u a " pu 0 , u 1 , 0, u 3 q in the locally Minkowskian frame, and U µ " pU t , U r , 0, U φ q " dx µ {dτ is the neutrino four velocity in the world coordi-nates. For example, U t " rr 3`a2 pr`r g qsE´aLr g mrrrpr´r g q`a 2 s .
The expressions for U r and U φ can be found on the basis of Eqs. (18) and (35).
Equations (2) and (3) are valid for the motion in the equatorial plane of a Kerr BH. Using Eqs. (16)- (18) and applying the similar technique as in Sec. 2.1, we get the expression for the angle α`8, which determines the probabilities of spin oscillations, in the form, where z " a{r g and x m is the maximal root of the equation, One can check that, at z Ñ 0, F K pxq Ñ F S pxq and R K pxq Ñ R S pxq, i.e., using Eqs. (20) and (21), we reproduce the spin evolution in the Schwarzschild metric, studied in Sec. 2.1.
It should be noted that, while studying the scattering off a rotating BH, we should distinguish the cases of the direct scattering with L ą 0 and the retrograde one L ă 0 (see Fig. 9 in Appendix B). We can formally suggest that the impact parameter is positive in the retrograde scattering, but consider the expressions similar to Eqs. (20) and (21), where all terms with odd powers of y have opposite signs.
As in Sec. 2.1, it is interesting to study the case of ultrareletivistic particles. Using Eqs. (20) and (21) in the limit γ " 1, we get that and The simultaneous solution of the equations R K pxq " 0 and R 1 K pxq " 0 gives one the critical impact parameter, where the upper signs stay for the direct scattering and the lower ones for the retrograde scattering. Note that, we made that both y 0 in Eq. (24) are positive. If y ă y 0 , a neutrino asymptotically falls to BH. The transition probability of spin oscillations P LR " p1`cos α`8q{2, computed on the basis of Eqs. (19), (22), and (23) turns out to be nonzero. It means that there is a possibility of transitions between left and right polarized ultrarelarivistic neutrinos interacting with a rotating BH. The corresponding transition probabilities are shown in Fig. 2.
Putting z " 0 in Eqs. (19), (22), and (23), one reproduces Eq. (11). Thus, using the results of Sec. 2.1, one gets that the transition probability of spin oscillations is vanishing at a small angular momentum of BH. The same feature results from the comparison of Figs. 2(a) and 2(b) with Figs. 2(c) and 2(d).

Neutrino gravitational scattering accounting for the matter interaction
In this section, we formulate the neutrino spin evolution equation in background matter under the influence of a gravitational field when a neutrino scatters off BH. Then, we derive the effective Schrödinger equation for spin oscillations of scattered neutrinos. Using the forward scattering approximation, one gets that the neutrino interaction with background matter is described by the following effective Lagrangian in the Minkowski spacetime [26]: where ν is the neutrino bispinor, γ µ and γ 5 are the Dirac matrices, and G F " 1.171 0´5 GeV´2 is the Fermi constant. The four vector G µ is the linear combination of the hydrodynamic currents and polarizations of background fermions. It depends on the chemical composition of matter and the type of a neutrino. The explicit form of G µ can be found in Ref. [27]. Basing on Eq. (25), the influence of the neutrino interaction with background matter on its spin evolution in curved spacetime was studied in Refs. [8,9]. It results in the appearance of the additional components of the vector Ω g in Eq. (5): Ω g Ñ Ω " Ω g`Ωmatt , where the vector Ω matt has the form, where g a " pg 0 , gq " e a µ G µ is the four vector of the effective neutrino interaction with background matter in the locally Minkowskian frame with coordinates x a , e a µ " Bx a {Bx µ are the vierbein vectors, and G µ is the analogue of the four vector in Eq. (25) given in curved spacetime.
If we study the neutrino interaction with nonrelativistic unpolarized matter, only G 0 " n eff U t f ‰ 0, where U t f is the time component of the four velocity of plasma. For spin oscillations of electron neutrinos in the electrically neutral hydrogen plasma one has n eff " n e , where n e is the electron number density. The expressions for n eff for other neutrino oscillations channels and various types of background fermions can be found in Ref. [27].
Instead of dealing with Eq. (1) for the spin precession, it is convenient to study the neutrino polarization density matrix, ρ " 1 2 r1`pσζqs, which obeys the equation, i 9 ρ " rH, ρs, where H "´pσΩq and Ω includes both the gravity and matter contributions in Eqs. (5) or (15), and (26). Here σ are the Pauli matrices.
Since the Liouville-von Neumann equation for the density matrix is rather complicated for the analysis, we can use the Schrödinger equation, i 9 ψ " Hψ. As we mentioned in Sec. 2, neutrinos move along the first axis in the locally Minkowskian frame at r Ñ 8. Hence, it is convenient to use this axis for the spin quantization. It mean that we should replace the Hamiltonian H Ñ U 2 HU : 2 , where U 2 " exppiπσ 2 {4q. This procedure brings the meaning to the effective wave function ψ. As in Secs. 2.1 and 2.2, it is convenient to rewrite the Schrödinger equation using the normalized radial coordinate x " r{r g , where Ω x " r g Ωdt{dr. Now we are ready to write down the effective Hamiltonian H x for both the Schwarzschild and Kerr metrics. One has that U t f " p1´r g {rq´1 and U f " 0 for nonrelativistic plasma near a nonrotating BH. We found in Ref. [7] that e 0 µ " p a 1´r g {r, 0, 0, 0q is the vierbein vector in the Schwarzschild metric. Then, in Eq. (26), g 0 " n eff p1´r g {rq´1 {2 and g " 0.
Using Eqs. (6) and (33), we obtain all the components of the vector Ω x in the form, where V " G F n eff r g { ? 2 is the dimensionless effective potential, the functions F S pxq and R S pxq are given in Eqs. (8) and (9), the upper and the lower signs stay for outgoing and incoming neutrinos.
In the case of a rotating BH, the four velocity of nonrelativistic matter has the components, and U f " 0. The four vector g a was obtained in Ref. [8] in the form, g a " n eff U t f ? r a r 3`a2 pr`r g qˆb rpr´r g q`a 2 , 0, 0,´a r g r˙, where we take again that matter is nonrelativistic. Basing on Eqs. (17), (30), and (35), we obtain the components of Ω x in the Kerr metric as where the functions F K pxq and R K pxq are given in Eqs. (20) and (21). We revealed in Sec. 2.2 that there is a spin conversion of ultrarelativistic neutrinos scattered off a Kerr BH. Thus, it is reasonable to rewrite Eq. (31) for such particles. We get when γ " 1. The expression for Ω x2 straightforwardly results from Eq. (22). In Eq. (32), the function R K pxq is given by Eq. (23). Equation (27) should be supplied with the initial condition ψ T 8 " p1, 0q, which means that all incoming neutrinos are left polarized. Since the neutrino velocity u changes the direction at t Ñ`8, the transition probability reads P LR " |ψ p1q 8 | 2 , and, correspondingly, the survival probability is P LL " |ψ The solution of Eq. (27), with Ω x in Eqs. (28) or (31), can be found only numerically because of the nontrivial dependence of Ω x on x. Moreover, in Sec. 4, we discuss the situation when n eff " n eff prq, or V " V pxq, which makes the analysis more complicated.
We also mention, that we cannot integrate Eqs. (27) to the turn point x m and then automatically reconstruct ψ`8, as we made in Secs. 2.1 and 2.2 to find α`8. In the presence of the background matter, the neutrino spin precesses around the axis with nonconstant direction. Moreover, the components Ω x2,3 in Eqs. (28) and (31)

Astrophysical applications
In this section, we present the numerical solutions of Eq. (27) for the neutrino scattering off SMBH surrounded by an accretion disk. We discuss both rotating and nonrotating SMBH, as well as the cases of ultrarelativistic neutrinos and neutrinos having a finite energy. The measurable neutrino fluxes are obtained.
First we notice, that standard model neutrinos are produced as left polarized particles. If they gravitationally interact with BH, some incoming left neutrinos become right polarized after scattering. A neutrino detector can observe only left neutrinos. Hence, the observed flux of neutrinos is F ν " P LL F 0 , where F 0 is the flux of scalar particles. The value of F 0 is proportional to the differential cross section, F 0 " dσ{dΩ, which is studied in Appendices A and B.
We assume that the neutrino beam scatters off a SMBH surrounded by an accretion disk. For example, we can suppose that such a SMBH is in the center of a Seyfert galaxy. We take that the plasma density in the disk scales as n e 9r´β. The value of β is very model dependent. For example, β « 0.5 in an advection dominated accretion disk studied in Ref. [28]. If we take that the mass of SMBH in question is M " 10 8 M d , the plasma density in the vicinity of SMBH can be up to n e " 10 18 cm´3 [29]. Thus, the dimensionless effective potential V prq " G F n e prqr g { ? 2, reads V pxq " V max x´β, where x " r{r g . The plasma motion in an accretion disk is driven not only by the gravitational interaction with a central BH. It also depends on the interaction between charged particles. Thus, the plasma angular velocity can be a complicated function of r, which is model dependent. We omit this additional factor in the description of the neutrino spin evolution and assume that the plasma motion in an accretion disk is nonrelativistic. First, in Fig. 3, we show the ratio of the fluxes of massive neutrinos, F pgq ν , scattered off a nonrotating BH and scalar particles F 0 of the same mass, with spin oscillations of neutrinos in the gravitational field being accounted for. These fluxes are proportional to the differential cross section, F " dσ{dΩ, where dΩ " 2π sin χdχ. The calculation of the cross section for scalar particles is presented in Appendix A and the result is shown in Fig. 8. We do not take into account the neutrino interaction with matter in Fig. 3.
The neutrino energy is E " 10m in Fig. 3. We have seen in Sec. 2.1 that neutrino spin oscillations in the scattering in the Schwarzschild metric are vanishing if E " m. Thus, the further enhancement of the neutrino energy will result in F pgq ν practically coinciding with F 0 . Consideration of neutrino energies smaller than in Fig. 3 is inexpedient from the point of view of possible astrophysical applications accounting for the current upper bound on neutrino masses in Ref. [30].
One can see in Fig. 3 that F pgq ν ă F 0 . This fact is owing to the survival probability P LL ă 1, which F pgq ν is proportional to, for the gravitational scattering of massive neutrinos. The maximal difference between the fluxes is for the backward neutrino scattering at χ " π. However, the maximal deviation of F pgq ν from F 0 is less than 1%. Accounting for the rather small neutrino energy, it makes difficult to observe the effect of spin oscillations in the neutrino scattering off a nonrotating BH. Now we study the influence of the interaction with an accretion disk matter on the flux of neutrinos scattered off a nonrotating SMBH with M " 10 8 M d . In Fig. 4, we show the result of the numerical solution of Eqs. (27) and (28). We take that n pmaxq e " 2ˆ10 18 cm´3 at r " r g and β " 0.2. The flux of scattered neutrinos with E " 10m accounting for the matter interaction F pdiskq ν , normalized by the flux of scalar particles, is shown in Fig. 4(a). The ratio of F pgq ν and F pdiskq ν is presented in Fig. 4(b). One can see in Figs. 3 and 4 that F pgq ν ă F pdiskq ν ă F 0 . The former inequality results from the fact that the matter interaction makes neutrino spin oscillations to be out of the resonance. Hence P pgq LL ă P pdiskq LL . The latter inequality is a consequence of P pdiskq LL ă 1. As in Fig. 3, the  Fig. 4 are at χ " π, Nevertheless, taking into account the magnitude of these deviations, one concludes that neutrino spin oscillations are negligible in the gravitational scattering off a nonrotating BH. Now we turn to the consideration of spin oscillations of neutrinos scattered off a rotating BH. We have revealed in Sec. 2.2 that ultrarelativistic neutrinos can change their polarization in this situation. Moreover, using Eqs. (28) and (31), one gets that the effective Hamiltonian for spin oscillations in the Kerr metric coincides with that for the Schwarzschild metric at z Ñ 0. Therefore, to highlight the manifestation of spin effects in the gravitational scattering off a Kerr BH, we consider ultrarelativistic neutrinos and a maximally rotating BH with z " 1{2 or a " M .
In Fig. 5, we show the ratios of the fluxes of gravitationally scattered ultrarelativistic neutrinos F pgq νd,r and scalar particles F 0d,r in the case of a rotating BH. We present the case of the direct scattering in Fig. 5(a) and the retrograde one in Fig. 5(b). As for a Schwarzschild BH, F pgq νd,r " F 0d,r at χ " 0. However, F pgq νd,r is about 10% less than F 0d,r for backwardly scattered neutrinos with χ " π. Now we take into account the neutrino interaction with an accretion disk. We consider a situation of a maximally rotating SMBH with M " 10 8 M d surrounded by an accretion disk with n pmaxq e " 10 18 cm´3 at r " r g and β " 0.5. In Fig. 6, we depict the fluxes of ultrarelativistic neutrinos F pdiskq νd,r , accounting for the interaction with background matter normalized by F 0d,r , as well as the ratios of F pgq νd,r and F pdiskq νd,r . We present the cases of the direct scattering in Figs. 6(a) and 6(b), and the retrograde one in Figs. 6(c) and 6(d).
We can see in Fig. 6 that the fluxes F pdiskq νd,r and F 0d,r , as well as F pgq νd,r and F pdiskq νd,r , differ by about 5%. The maximal difference is at χ " π. The inequality F pgq ν ă F pdiskq ν ă F 0 , established for the Schwarzschild metric above, remains valid for a Kerr BH as well.
Now it is interesting to compare the direct and retrograde scatterings. We show the ratios of the corresponding fluxes in Fig. 7. One can see that that the difference between the fluxes for the retrograde and direct scatterings can be about 20%. This asymmetry remains valid for both the only gravitational scattering, shown in Fig. 7(a), and when the matter contribution is accounted for, which is depicted in Fig. 7(b).

Discussion
In the present work, we have studied spin oscillations in the neutrino scattering off BHs. Both nonrotating and rotating BHs have been discussed. The neutrino spin evolution in curved spacetime has been accounted for quasiclassically basing on the approach developed in Refs. [7,8]. As an application of the obtained results, we have examined the neutrino scattering off SMBH with a realistic accretion disk. In Sec. 2, we have derived the general expressions for the transition and survival probabilities for neutrino spin oscillations. Then, in Sec. 2.1, we have applied these results for the neutrino scattering off a nonrotating BH. The expression for the spin rotation angle, which is valid for the arbitrary neutrino energy and the impact parameter, has been presented in Eqs. (8) and (9). Using this result, we have obtained that an ultrarelativistic neutrino preserves its helicity while scattering in the Schwarzschild metric. This feature is valid for any impact parameter y ě y 0 . This our finding is in agreement with the results of Ref. [23]. The transition probabilities P LR of spin oscillations of massive neutrinos in their gravitational scattering turn out to be nonzero, and are shown in Fig. 1 for several neutrino energies. However, P LR is rather small for reasonable Lorentz factors E{m to be observed.
Then, in Sec. 2.2, we have considered spin oscillations of neutrinos scattered off a rotating BH. As in Sec. 2.1, we have derived the general formula for the spin rotation angle in Eqs. (19)- (21). The obtained expressions reproduce the results of Sec. 2.1 in the limit z Ñ 0. The transition probabilities of spin oscillations are nonzero even for ultrarelativistic neutrinos. In Fig. 2, we have shown P LR for the direct and retrograde scatterings of ultrarelativistic neutrinos for different angular momenta of BH.
Note that the fact that the helicity of ultrarelativistic, or even massless, fermions can be changed under the influence of a gravitational field was noticed earlier in Refs. [31,32]. .
We have found that spin oscillations of ultrarelativistic neutrinos happen only when particles scatter off a rotating BH. The polarization change of ultrarelativistic particles (photons) in their gravitational scattering off a rotating star was studied previously in Ref. [33]. Then, in Sec. 3, we have derived the effective Schrödinger equation for a neutrino scattering off BH surrounded by background matter with a nonuniform density. In the case of only the gravitational scattering, studied in Secs. 2.1 and 2.2, it was possible to obtain the analytical transition and survival probabilities for some impact parameters. If, besides gravity, a neutrino interacts with a background matter, the probabilities can be derived only in the numerical solution of Eq. (27).
The effective Hamiltonians for neutrino spin oscillations for both the Schwarzschild and Kerr metrics have been obtained in Sec. 3 in the approximation of a slowly rotating accretion disk. The effective Hamiltonian for ultrarelativistic neutrinos scattering off a rotating BH surrounded by an accretion disk has been also derived in Sec. 3.
In Sec. 4, we have considered the astrophysical applications of our results. In particular, we have studied the effect of spin oscillations on the neutrino scattering off SMBH surrounded by an accretion disk. We have taken the parameters of the accretion disk, such as the maximal number density and the profile of the mass distribution, close to the values resulting from observations and hydrodynamics simulations.
First, we have studied the case of a nonrotating BH. Using the numerical solution of Eqs. (27) and (28), we have found the observed fluxes of outgoing neutrinos for only the gravitational scattering and when the neutrino interaction with the accretion disk is accounted for. The contribution of spin oscillations to the neutrino fluxes, shown in Figs. 3 and 4, appears to be negligible for reasonable neutrino energies and the current upper bound on the neutrino mass.
Greater spin effects have been revealed for neutrinos gravitationally scattered off a rotating BH surrounded by an accretion disk. In this situation, we have studied the ultrarelativistic neutrinos scattering off a maximally rotating SMBH to highlight the effect of spin oscillations. We have shown in Fig. 5 that the observed fluxes of gravitationally scattered neutrinos can be reduced by almost 10%. The contribution of the neutrino interaction with matter changes the fluxes by about 5%, see Fig. 6.
As one can see in Figs. 3-6, there is no deviation of the fluxes for the forward neutrino scattering at χ " 0 if one compares them with the fluxes of scalar particles. It means that neutrino spin oscillations do not affect the size of a BH shadow. The major effect of spin oscillations is for the backward neutrino scattering at χ " π. Thus the intensity of the glory flux for ultrarelativistic neutrinos is almost 10% less than for scalar particles in case a rotating BH.
The influence of the plasma interaction on the gravitational scattering of photons was thoroughly studied previously [21]. The photons propagation in plasma surrounding a nonrotating BH was examined in Ref. [34], where it was found that the form of the BH shadow is not changed. However its size can be enlarged. The shape of the shadow of a rotating BH can be deformed [34]. Using Fig. 4(b), we conclude that there is an asymmetry in the observed neutrino fluxes depending on the orientation of the neutrino trajectory with respect to a slim accretion disk. This asymmetry is maximal for the backward neutrino scattering. Unfortunately, this effect is rather small in the Schwarzschild metric. The asymmetry in the outgoing neutrino fluxes may well exist for the Kerr BH surrounded by an accretion disk. Basing on our results, we cannot quantitatively describe this effect since we rely on the equatorial neutrino motion only.
We have found that spin oscillations in the gravitational scattering off a rotating BH are sizable for ultrarelativistic neutrinos. Thus our results are of interest for the neutrino astronomy [35], which is a rapidly developing area of the cosmic rays physics. It is known that neutrinos with energies in the PeV range were detected [36]. Moreover, several sources of ultrahigh energy neutrinos can be identified with with some astronomical objects such as active galactic nuclei [37]. We have demonstrated in our work that, if the incoming flux of cosmic neutrinos experience the gravitational lensing, in some cases, the observed flux can be reduced by down to 10%, compared to its initial value, because of neutrino spin oscillations.
There is another possible application of the obtained results. The r-process nucleosynthesis in the vicinity of BH, surrounded by an accretion disk, was found in Ref. [18] to be affected by the neutrino radiation of the accretion disk. To influence the weak nuclear reactions, which the nucleosynthesis is based on, emitted neutrinos should be active. In the present work, we predict the significant conversion of left active neutrinos to right sterile particles in the neutrino scattering off a rotating BH. Thus the nucleosynthesis near such BHs is further modified because of the neutrino gravitational interaction. Note that the effect of active to sterile neutrinos oscillations on the dynamics of the supernova explosion was studied in Ref. [38]. The differential cross section of the gravitational scattering of scalar particles off a nonrotating BH, normalized by r 2 g , versus χ{π. The particle energy is E " 10m.
While computing the differential cross section, dσ{dΩ, where dΩ " 2π sin χdχ, we should take into account that a particle, before being scattered off, can make multiple revolutions around BH, both clockwisely and anticlockwisely. One should account for this fact in the determination of the angle χ, fixing the position of a detector, which is in the range 0 ă χ ă π.
In Fig. 8, we present the result of the numerical computation of the cross section. While building this plot, we take that y 0 ă y ă 30y 0 and account for up to two revolutions of a particle around BH in both directions. This our result is used in Sec. 4 when we study the neutrino scattering off a realistic BH surrounded by an accretion disk.

B Particle motion in the Kerr metric
In this Appendix, we describe the motion of a scalar particle in the gravitational field of a rotating BH. The trajectory of a particle in the Kerr metric, given in Eqs. (13) and (14), is rather complicated in the general case. However, if we consider the particle motion in the equatorial plane, the trajectory and the equation of motion can be expressed in the relatively simple form (see, e.g., Ref. [25]), dt dr "˘? r rpr´r g q`a 2 rr 3`a2 pr g`r qsE´aLr g a rr 3`a2 pr g`r qsE 2´2 aLEr g´p r´r g qL 2´m2 rrrpr´r g q`a 2 s , dφ dr "˘? r rpr´r g q`a 2 pr´r g qL`aEr g a rr 3`a2 pr g`r qsE 2´2 aLEr g´p r´r g qL 2´m2 rrrpr´r g q`a 2 s , where one has two integrals of motion: the particle angular momentum L and its energy E. Analogously to Eq. (34), we get the expression for φ m , where R K pxq is given in Eq. (21), and x m is the maximal root of Eq. (21). There is an important difference between the particle scattering in the Schwarzschild and Kerr metrics. We can take that L ą 0 in the Schwarzschild metric. In the Kerr metric, the cases L ą 0 and L ă 0 are different. It can be illustrated in Fig. 9. One can call the situation, when the position of a detector is bent towards the BH rotation direction, as the direct scattering. It is shown in Fig. 9(a). The opposite situation, depicted in Fig. 9(b) can be called the retrograde scattering. The corresponding quantities for the retrograde scattering can be obtained by replacing the signs in all terms with odd powers of y.
Of course, the total flux of, e.g., directly scattered particles contains contributions with L ă 0. Indeed, despite 0 ă χ ă π in Fig. 9(a), particles, which make one full revolution around BH, move oppositely the BH rotation, i.e. they have L ă 0. This fact should be accounted for while computing the total flux of scattered particles.
The calculation of the differential cross section of the gravitational scattering of scalar particles off a Kerr BH is analogous to the Schwarzschild metric studied in Appendix A. Thus we omit the details.