Partial wave analysis of the Dirac fermions scattered from Schwarzschild black holes

Asymptotic analytic solutions of the Dirac equation, giving the scattering modes (of the continuous energy spectrum, E>mc2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E>mc^2$$\end{document}) in Schwarzschild’s chart and Cartesian gauge, are used for building the partial wave analysis of Dirac fermions scattered by black holes. In this framework, the analytic expressions of the differential cross section and induced polarization degree are derived in terms of scattering angle, mass of the black hole, and energy and mass of the fermion. Moreover, the closed form of the absorption cross section due to the scattering modes is derived showing that in the high-energy limit this tends to the event horizon area regardless of the fermion mass (including zero). A graphical study presents the differential cross section analyzing the forward/backward scattering (known also as glory scattering) and the polarization degree as functions of scattering angle. The graphical analysis shows the presence of oscillations in scattering intensity around forward/backward directions, phenomena known as spiral scattering. The energy dependence of the differential cross section is also established by using analytical and graphical methods.


Introduction
The complex problem of the quantum particles scattered from black holes was studied extensively, mainly considering massless scalar [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] or other massless boson fields [17][18][19], since the equations of massive fields cannot be solved analytically in the Schwarzschild geometry. This is one of the reasons why the scattering of massive Dirac fermions by black holes was studied either in particular cases [20,21] or by using combined analytical and numerical methods [22][23][24][25][26]. Thus in Refs. [25,26] such methods were applied for investigating the absorption cross section, the scattering a e-mail: cota@physics.uvt.ro b e-mail: crucean@physics.uvt.ro c e-mail: ciprian.sporea89@e-uvt.ro intensity (pointing out the glory and spiral scattering) and the polarization degree of these scattering processes. However, in the case of the interactions of quantum particles with black holes the numerical methods may represent a difficult task since one must combine quantities at quantum scale (mass, energy of the fermion) with quantities at galactic scale (as the black hole mass). Therefore, it is obvious that any new analytic study may improve this investigation, helping us to understand the quantum mechanisms governing the scattering process.
For this reason we would like to propose in this paper the partial wave analysis of the Dirac fermions scattered from Schwarzschild black holes constructed applying exclusively analytical methods. In order to do this, we exploit the analytical properties of the (approximative) asymptotic solutions of the Dirac equation we have found some time ago [27]. These solutions were obtained in the chart with Schwarzschild coordinates where we considered the Cartesian gauge that preserves the global central symmetry of the field equations [28], allowing the separation of spherical variables just as in the central problems of special relativity [29]. We note that thanks to this gauge one of us (IIC) succeeded to solve analytically the Dirac equation on the central charts of the de Sitter and anti-de Sitter spacetimes [30][31][32][33]. In the case of the Schwarzschild chart and Cartesian gauge, after the separation of the angular variables, we remain with a pair of simple radial equations depending only on the gauge field components, which can be approximatively solved by using the Novikov radial coordinate [34,35]. We obtained thus the asymptotic radial solutions that are either of scattering type or even possible bound states [27]. However, these last mentioned states are less studied and seem to be unstable [36], so that we restrict ourselves to considering only the scattering solutions (with E > mc 2 ) for performing the partial wave analysis of the Dirac fermions scattered from black holes.
In general, when one considers the analytic expressions of the scattering quantum modes, one must impose suitable boundary conditions (in the origin or at the event horizon) in order to fix the integration constants determining the asymptotic behavior and implicitly the phase shifts of the partial wave analysis. Having here only the asymptotic form of the Dirac spinors in the black hole field we must replace the boundary conditions with suitable asymptotic conditions playing the same role in determining the integration constants of our asymptotic solutions. Fortunately, this can be done assuming that for large values of the angular momentum the collision becomes elastic approaching the Newtonian limit [20,21]. This hypothesis is enough for determining completely the integration constants and deriving the analytical forms of the phase shifts, scattering amplitudes, differential cross section, and polarization degree. The nice surprise is that this approach emphasizes the absorption of the fermions by black holes in a natural manner, laying out the absorption cross section in an analytical closed form with some associated selection rules.
The paper is organized as follows. In Sect. 2 we briefly present our method of separating the variables of the Dirac equation on central backgrounds and Cartesian gauge focusing on the mentioned approximative solutions of the Dirac equation in Schwarzschild's charts, corresponding to the continuous energy spectrum [27]. Section 3 is devoted to our partial wave analysis based on the asymptotic condition (discussed in Appendix C), which fixes the integration constants giving the analytical form of the phase shifts encapsulating both cases studied here, the elastic scattering of fermions and their absorption by a black hole. We verify that our scattering amplitudes have the correct Newtonian limit and we derive the scattering intensity and polarization degree. Special attention is paid to the absorption cross section, for which we give the analytic expression of the partial cross sections and the selection rules indicating in which partial wave we can find absorption. Section 4 is devoted to the graphical analysis and discussion of the physical consequences of our results. Here we show that our analytical results concerning the elastic scattering are very similar to those obtained by using analytical-numerical methods [25,26], but there are some differences as concerns the absorption. Other conclusions are summarized in the last section.
In the following we use natural units with c =h = G = 1.

The Dirac equation in central charts and Cartesian gauge
In a given frame {x; e}, the Dirac equation of a free spinor field ψ of mass m has the form where γα are the point-independent Dirac matrices that satisfy {γα, γβ } = 2ηαβ and define the generators of the spinor representation of the SL(2, C) group, Sαβ = i 4 [γα, γβ ], which give the spin connections of the covariant derivatives, depending on the connection components in local frameŝ σ μν = e α μ e β ν (êσ γ γ αβ −êσ β,α ) where the notation γ αβ stands for the usual Christoffel symbols.
In this approach the Dirac equation (1) takes the explicit form where g = det(g μν ). Moreover, from the conservation of the electric charge, one deduces that the time-independent relativistic scalar product of two spinors [30] ( is given by the integral over the space domain D of the local chart under consideration. In general, a manifold with central symmetry has a static central chart with spherical coordinates (t, r, θ, φ), associated to the Cartesian ones (t, x), with r = | x|, covering the space domain D = D r × S 2 , i.e. r ∈ D r while θ and φ cover the sphere S 2 . The form of the Dirac equation in any chart is strongly dependent on the choice of the tetrad gauge. For this reason, our approach is based on the virtues of the mentioned Cartesian gauge which is defined by the 1-forms [28,30], expressed in terms of three arbitrary functions of r , denoted by u, v, and w, which allow us to write the general line element We have shown that in this gauge the last term of Eq. (3) does not contribute and, moreover, there is a simple transformation, ψ → vw − 3 2 ψ, able to eliminate the terms containing the derivatives of the functions u, v, and w, leading thus to a simpler reduced Dirac equation [30]. The advantage of this equation is that its spherical variables can be separated just as in the case of the central problems in Minkowski spacetime [29]. Consequently, the Dirac field can be written as a linear combination of particular solutions of given energy, E. Those of positive frequency, are particle-like energy eigenspinors expressed in terms of radial wave functions, f ± E,κ , and usual four-component angular spinors ± m j ,κ [29]. It is well known that these spinors are orthogonal to each other, being labeled by the angular quantum numbers m j and which encapsulates the information as regards the quantum numbers l and j = l ± 1 2 as defined in Refs. [29,37] (while in Ref. [26] κ is of opposite sign). The spherical spinors are normalized to unity with respect to their own angular scalar product. We note that the antiparticle-like energy eigenspinors can be obtained directly using the charge conjugation as in the flat case [39].
Thus the problem of the angular motion is completely solved for any central background. We are left with a pair of radial wave functions, f ± (denoted from now on without indices), which satisfy two radial equations that can be written in compact form as the eigenvalue problem H r F = EF of the radial Hamiltonian [30], in the space of two-component vectors (or doublets), F = ( f + , f − ) T , equipped with the radial scalar product [30] ( resulting from the general formula (4) where we have to take √ −g(x) = w(r ) 4 u(r )v(r ) 2 r 2 sin θ , the spinors U and U of the form (10) and e 0 μ γμ = 1 w(r ) γ 0 . This scalar product has to select the 'good' radial wave functions, i.e. tempered distributions, which enter in the structure of the particle-like energy eigenspinors.
In the central charts each particular solution (10) gives rise to a partial radial current defined as where we take into account that in our Cartesian gauge we have Exploiting then the property γ 0 γ x ± = ±i ∓ and the orthogonality of the spherical spinors [29] we obtain the final result Now it is obvious that these currents are conserved and independent on r , since ∂ r J rad = 0 whenever the functions f ± satisfy the above radial equations. These will represent useful significant constants of motion in the scattering process.

Analytic asymptotic solutions in Scwarzschild's charts with Cartesian gauge
Let us assume that a Dirac particle of mass m is moving freely (as a perturbation) in the central gravitational field of a black hole of mass M with the Schwarzschild line element defined on the radial domain D r = (r 0 , ∞) where r 0 = 2M. Hereby we identify the functions which give the radial Hamiltonian (12). The resulting radial problem cannot be solved analytically as it stays forcing one to resort to numerical methods [26] or to some approximations. In Ref. [27] we proposed an effective method of approximating this radial problem expanding the radial equations in terms of the Novikov dimensionless coordinate [34,35], Using this new variable and introducing the notations we rewrite the exact radial problem as ⎛ Moreover, from Eq. (13) we find that the radial scalar product now takes the form We note that the singular form above, due to our special parametrization and use of the Novikov variable, is merely apparent, since this comes from the usual regular scalar product (4). Moreover, we have shown that there exist square integrable spinors with respect to this scalar product [27].
For very large values of x, we can use the Taylor expansion with respect to 1 x of the operator (21) neglecting the terms of the order O(1/x 2 ). We obtain thus the asymptotic radial problem [27], which can be rewritten as after reversing between themselves the lines of the matrix operator. This is necessary for diagonalizing simultaneously the term containing derivatives and the one proportional to x, as in the Dirac-Coulomb case [37]. This can be done by using the matrix for transforming the radial doublet as where ν = 2 − μ 2 . These equations can be solved analytically for any values of . In Ref. [27] we derived the spinors corresponding to the continuous spectrum ∈ [μ, ∞), but without investigating scattering effects.
Here we would like to complete this study by calculating the fermion scattering amplitudes starting with the asymptotic solutions of Ref. [27] we briefly present below. The radial equations (25) can be analytically solved for the continuous energy spectrum, > μ, in terms of the Whittaker functions as [27] f where we denote The integration constants must satisfy [27] C − We observe that these solutions are similar to those of the relativistic Dirac-Coulomb problem. The functions

Partial wave analysis
We consider now the scattering of Dirac fermions on a black hole. This is described by the energy eigenspinor U whose asymptotic form, for r → ∞ (where the gravitational field vanishes) is given by the plane wave spinor of momentum p and the free spherical spinors of the flat case behaving as Here we fix the geometry such that p = p e 3 , while the direction of the scattered fermion is given by the scattering angles θ and φ, which are just the spheric angles of the unit vector n. Then the scattering amplitude depends on two scalar amplitudes, f (θ ) and g(θ ), which can be studied by using partial wave analysis.

Asymptotic conditions and phase shifts
The partial wave analysis exploits the asymptotic form of the exact analytic solutions which satisfy suitable boundary conditions that in our case might be fixed at the event horizon (where x = 0). Unfortunately, we have here only the asymptotic solutions (26) and (27) whose integration constants cannot be related to those of the approximative solutions near event horizon [27] without resorting to numerical methods [20][21][22][23][24][25][26]. Therefore, as long as we restrict ourselves only to an analytical study, we cannot match boundary conditions near the event horizon, being forced to find suitable asymptotic conditions for determining the integration constants. Then the behavior near horizon can be discussed, a posteriori, focusing on the partial radial currents and especially on the absorption cross section, which can easily be derived in our approach. When the Dirac fermions are scattered from the central potentials, the particular spinors of the parameters ( p, κ) represent partial waves having the asymptotic form where the phase shifts δ κ can be, in general, complex numbers which allow one to write down the amplitudes in terms of the Lagrange polynomials of cos θ . Each partial wave gives rise to the radial current (16), which can be put now in the form The scattering is elastic when the phase shifts δ κ are real numbers such that the partial radial currents vanish. Whenever there are inelastic processes the phase shifts become complex numbers, whose imaginary parts are related to the inelastic behavior, producing the non-vanishing radial currents indicating absorption.
In Appendix C we show that in our approach it is necessary to adopt the general asymptotic condition C + 2 = C − 2 = 0 in order to have elastic collisions with a correct Newtonian limit for large angular momentum. It is remarkable that this asymptotic condition selects the asymptotic spinors that are regular in x = 0; but this cannot be interpreted as a boundary condition since the exact solutions have a different structure near the event horizon [27].
Under such circumstances we have to use the asymptotic forms of the radial functions for large x resulting from Eq. (A2). We observe that the first term from (A2) is dominant forf + , while the second term in (A2) is dominant forf − , such that for x → ∞ the radial functions behave aŝ Hereby we derive the asymptotic form (33) of the doublet F = TF observing that the argument of the trigonometric functions can be deduced as 1 2 arg Then by using Eq. (29a) and taking into account that Eqs. (19) and (31b) allow us to replace νx 2 = p(r − r 0 ), we obtain the definitive asymptotic form of the radial functions for the scattering by black holes, whose point-independent phase shifts δ κ give the quantities Notice that the values of κ and l are related as in Eq. (11), i.e. l = |κ| − 1 2 (1 − sign κ). The remaining radially dependent phase, which does not depend on angular quantum numbers, may be ignored as in the Dirac-Coulomb case [26,37].
We arrived thus at the final result (38) depending on the parameters introduced above, which can be expressed in terms of physical quantities by using Eq. (31b) as The parameters k, q, λ ∈ R + are positively defined and satisfy the identity In the particular case of the massless fermions (m = 0) we are left with the unique parameter k = q = 2 pM, since λ = 0.
The parameter s has a special position since this can take either real values or pure imaginary ones. Let us briefly discuss these two cases.
The case of elastic scattering arises for the values of κ (at given p) that satisfy the condition where n = floor(k) is the largest integer less than k. Then s = √ κ 2 − k 2 ∈ R and the identity (43) guarantees that the phase shifts of Eq. (38) are real numbers such that |S κ | = 1.
The case of absorption is present in the partial waves for which we have Here we meet a branch point in s = 0 and two solutions s = ±i|s| = ±i √ k 2 − κ 2 ; among them we must choose s = −i|s|, since only in this manner we select the physical case of |S κ | < 1. More specific, by substituting s = −i|s| in Eq. (38) we obtain the simple closed form showing that 0 < |S κ | < 1, since |s| < q for any ( p, κ) obeying the condition (45). Moreover, we can verify that in the limit of the large momentum (or energy) the absorption tends to become maximal, since regardless of the fermion mass.
Finally, we note that for any values of the above parameters the phase shifts δ l and δ −l are related as which means that for massless Dirac fermions we have δ l = δ −l . This is the basic framework of the relativistic partial wave analysis of the Dirac fermions scattered by Schwarzschild black holes in which we consider exclusively the contribution of the scattering modes. Our results are in accordance with the Newtonian limit since in the large-l limits and for very small momentum we can take s ∼ |κ| ∼ l and λ ∼ q such that our phase shifts (38) become just the Newtonian ones, Eq. (B1) [20,21].
On the other hand, we observe that our phase shifts are given by a formula which has the same form as that deduced for the Dirac-Coulomb scattering [37]. Indeed, by replacing in Eq. (38) the parameters we recover the phase shifts of the Dirac-Coulomb scattering. However, there are significant differences between these parameterizations showing that these two systems are of different natures.

Partial amplitudes and cross sections
Now we have all the elements for calculating the analytic expressions of the amplitudes and cross sections in terms of the phase shifts defined by Eq. (38). In the relativistic theory, the scalar amplitudes of Eq. (32), depend on the following partial amplitudes [26,37], In our problem the quantities S κ are given by the analytic expression (38) such that we can apply analytic methods for studying the above amplitudes, which give the elastic scattering intensity or the differential cross section, and the polarization degree, This last quantity is interesting for the scattering of massive fermions representing the induced polarization for an unpolarized initial beam. Let us verify first that our approach recovers the correct Newtonian limit in weak gravitational fields (with small M) where the partial amplitudes can be expanded as l + · · · l = 1, 2, . . . .
Our algebraic codes indicate that the cases l = 0 and l > 0 must be studied separately since the expansion of f l does not commute with its limit for l → 0. However, this is not surprising since a similar phenomenon can be met in the Dirac-Coulomb problem. Then, according to Eqs. (52), (53), and (38), we obtain first the expansion for l = 0, where γ is the Euler constant. Furthermore, for any l > 0 we find the terms of first order, depending on the digamma function ψ, and the more complicated ones of second order, which show the dependence on l and p. Then, from Eq. (B3) we observe that f are of Newtonian form, increasing with l because of the digamma function. The other terms decrease when l is increasing so that we verify again that our partial wave analysis has a correct Newtonian limit for large l, as happens in the well-known case of the scalar particles [20,21].
In other respects, we observe that the partial amplitudes of the massless fermions are regular in p = 0, while those of the massive particles diverge in any order. Thus the problem of removing the infrared catastrophe in the massive case seems to remain as a serious challenge.
Important global quantities, independent on θ , are the total cross sections. The elastic cross section, is obtained by integrating the scattering intensity (54) over θ and φ, according to the normalization integral The general expression (64) gives the elastic cross section even in the presence of absorption when |S κ | < 1. For the genuine elastic scattering (without absorption and |S κ | = 1), the elastic cross section σ e represents just the total cross section σ t such that we can write and one deduces that the absorption cross section, σ a = σ t − σ e , reads [26] σ a = π since for s = −i|s| we have |S −κ | = |S κ | as in Eq. (46). This last cross section deserves special attention, since this can be calculated at any time as a finite sum indicating how the fermions can be absorbed by black holes. This is a function of the momentum p since all our parameters, including n, depend on it. For this reason it is convenient to represent the absorption cross section as in terms of the partial cross sections whose definitive closed form, is derived according to Eqs. (46) and (67), while the condition (45) introduces the Heaviside step function θ(k − l). Hereby we understand that the absorption arises in the partial wave l for the values of p satisfying the condition k > l. This means that for large values of l there may appear non-vanishing thresholds, indicating that the fermions with |κ| = l can be absorbed by a black hole only if p ≥ p l . Obviously, when √ 3m M < 1 we have p l = 0 for any l such that we meet a non-vanishing threshold in every partial wave. This happens for massless fermions too when we have the equidistant thresholds p l = l 2M for any l ≥ 1.
The existence of these thresholds is important, since these keep under control the effect of the singularities in p = 0. Thus for any partial wave with p l > 0 the Heaviside function prevents the partial cross section to be singular, but when p l = 0 then p reaches the singularity point p = 0, where the partial cross section diverges. Obviously, for √ 3m M < 1 all the partial sections are finite.
Finally, we must specify that in the high-energy limit all these absorption cross sections tend to the event horizon (apparent) area, indifferent to the fermion mass m ≥ 0, as results from Eq.
since n(n + 1) ∼ k 2 ∼ 4M 2 p 2 . This asymptotic value is less than the geometrical optics value of 27π M 2 [38], which is the high-energy limit of the absorption cross sections obtained applying analytical-numerical methods [25,26]. The explanation could be that here we neglected the effects of possible unstable bound states [27,36] that may give rise to a resonant scattering.

Numerical examples
Our purpose now is to use the graphical analysis for understanding the physical consequences of our analytical results encapsulated in quite complicated formulas and the infinite series (50) and (51), which are poorly convergent or even divergent since their coefficients a l and b l are increasing with l faster than l ln l, respectively, l −1 . This is an unwanted effect of the singularity at θ = 0 but which can be attenuated adopting the method of Ref. [41], which consists in replacing the series (50) and (51) by the mth reduced ones, The recurrence relations satisfied by the Legendre polynomials P l (x) , P 1 l (x) lead to the iterative rules giving the reduced coefficients in any order, if we start with a (0) l = a l and b (0) l = b l as defined by Eqs. (52), (53) and (38). Note that in this last equation we replace s → s −i| s| in order to cover automatically the two cases of interest here, elastic collision and absorption. Then we will see that this method is very effective, ensuring the convergence of the reduced series for any value of θ apart from the singularity in θ = 0. Here we present the numerical results obtained by using the second iteration for f (m 1 = 2) and the first one for g (m 2 = 1), which seems to be satisfactory without distorting the analytical results.
In this approach the elastic and total cross sections cannot be calculated since their series remain divergent under any conditions because of the mentioned singularity. Even if we replace the amplitudes f and g with their reduced series we cannot ensure the convergence since then we introduce additional factors of the type (1 − cos θ) −n , n ≥ 2, leading to divergent integrals over θ when we calculate σ e as in Eq. (64). Therefore, we are able to study only the absorption cross section (69), which is given by Fig. 1 The differential cross section as a function of θ for different values of the fermion speed a finite sum, followed then by focusing on the scattering amplitudes, the scattering intensity, and the polarization degree.
All these quantities depend on three free parameters: the fermion mass m and momentum p, and the black hole mass M. Since we work in the asymptotic zone where the fermion energy is E = m 2 + p 2 , we can use the fermion velocity v = p/E as an auxiliary parameter. With these parameters one can construct two relevant dimensionless quantities, which in the usual units read G M E/(hc 3 ) and mG M/(hc). The last one can be seen as (proportional to) the ratio of the black hole horizon to the fermion Compton wavelength. In our natural units (with c =h = G = 1) these quantities appear as M E, respectively, m M, these being used for labeling our graphs.
We should mention that our graphical analysis is performed in the following only for the case of small or micro black holes, for which ME and m M take relatively small values, since in this manner we can compare our results with those obtained by using analytical-numerical methods [22][23][24][25][26]. However, our analytical results presented above are valid for any values of these parameters.

Forward and backward scattering
We begin the graphical analysis by plotting the differential cross section (54), as a function of the angle θ for different numerical values of m M and M E. Since E = m 2 + p 2 , the condition M E ≥ m M must be always satisfied. In addition, multiplying by M the expression of energy we obtain m M = M E √ 1 − v 2 , which gives the connection between the pair of parameters that define our analytical formulas. The graphs in Figs. 1 and 2 show how the scattering intensity depends on the scattering angle for small/large fermion velocities. In order to observe the oscillations in the scattering intensity around θ = π , corresponding to backward scattering, we shorten the axis of θ because the cross section is divergent in θ = 0.
In Figs. 1 and 2 we observe the presence of a maximum in scattering intensity in the backward direction. This is known as glory scattering [42], while the oscillations in the scattering intensity at intermediate scattering angles, around θ = π , are known as orbiting or spiral scattering [42]. Another important observation that emerges from these graphs is that the scattering intensity has large values only for small fermion Fig. 2 The differential cross section as a function of θ for different values of the fermion speed velocities, while at large ones the scattering intensity is observably smaller. The conclusion is that the glory and orbiting scattering are significant only for non-relativistic fermions.
At θ = 0 the scattering intensity becomes divergent for both small or large values of the fermion speed. For this reason, in order to obtain the behavior of scattering intensity at small scattering angles for different fermion velocities, we restrict ourselves only to values of θ close to 0.
Our graphs in Fig. 3 show that the scattering intensity becomes divergent as θ → 0 (forward scattering) and they also show the presence of oscillations around small values of θ , i.e. the orbiting scattering [42]. We also observe that the scattering intensity in the forward direction increases with the parameter M E. The forward scattering is in fact a diffraction on the black hole horizon. This could be possible only in the case when the wavelength of the incident particle is comparable with the size of the event horizon. From Fig. 3 we observe that the oscillatory behavior of scattering intensity in the forward direction is more pronounced in the case of small fermion velocities comparatively with the case of relativistic ones. It is also worth to mention that the oscillatory behavior of the scattering intensity is increasing with the parameter M E (see Fig.  4). The conclusion is that the orbiting scattering is negligible for relativistic particles, while in the case of particles with small velocities this phenomenon becomes important.
Furthermore, we address the problem of variation of the cross section with the black hole mass given in Fig. 4. We observe that the scattering intensity in the backward direction is increasing with the black hole mass but without changing its general profile. This suggests that the positions of the relative extremes are independent on the black hole mass.
Let us comment now on the results of the scattering intensities obtained in graphs Figs. 1, 2, 3, and 4. One knows that the differential cross section represents the area that the incident particle must cross in the target zone in order to be detected in the solid angle d . In the case of a scattering process between two quantum objects this quantity is very small. So in our case it is not surprising that the differential cross section becomes very large, since the target is of the size of the black hole event horizon. This can be better understood if we recall the result from classical physics according to which Fig. 3 The differential cross at small θ, for different values of the fermion speed the scattering intensity for a classical particle moving on a spiral trajectory is larger that the scattering intensity for a particle moving on a straight line. A classical particle which has nonzero angular momentum will always cross a larger area in the target zone. The situation is the same in the case of a quantum particle scattered by a black hole but with the observation that in this case the notion of trajectory is not well defined. Taking into consideration that the minimum area of the target is of the size of the event horizon we see that the area crossed by the fermion to be detected in a solid angle could be very large. As a final remark we can observe that the scattering intensity in the forward/backward direction increases with the black hole mass. This result is expected since the area of the event horizon also increases with the black hole mass. Our results are compatible with those obtained in the literature using analytical-numerical methods [26].

Dependence on energy
We study now the behavior of the differential cross section in terms of the energy by plotting Eq. (54) as a function of the ratio E/m for different scattering angles. Since for θ = 0 the differential cross section is divergent, our analysis is done for θ = π/3, π/4. Another observation that emerges from Figs. 5 and 6 is that the energy dependence of scattering intensity has a more pronounced oscillatory behavior for small scattering angles and large values of m M. On the other hand, we also observe from Figs. 5 and 6 that the scattering intensity is decreasing when the energy increases and becomes divergent in the limit of small energies. The shapes of these graphs are the result of the fact that our scattering intensity is proportional with the usual factor 1/E 2 and the oscillatory effect is given by the more complicated dependence of the energy on the phase shift given by Eq. (38).

Polarization degree
If we consider that the fermions from the incident beam are not polarized, then after the interaction with the black hole, the scattered beam could become partially polarized. It is interesting to study this effect by plotting the degree of polarization (55) as a function of scattering angle for given values of M E and different fermion velocities. Plotting the polar-  As we may observe from Figs. 7 and 8, the polarization is a very oscillatory function of the scattering angle. These oscillations are the result of the forward and backward scat-tering as well as the orbiting scattering. These three types of scattering induce the oscillatory behavior of polarization, since the scattering intensity also oscillates with the scattering angle. Another result that is worth mentioning is the oscillatory behavior of the polarization, Figs. 7 and 8, which  is modified as we change the parameter M E. If the fermion energy is fixed, then we can draw the conclusion that the oscillatory behavior of the polarization depends on the black hole mass M, becoming more pronounced as we increase the black hole mass.
To see how the spin of the fermion is aligned with a given direction after scattering off a black hole, we present the polar plots for the degree of polarization in Figs. 9, 10. These graphs are plotted for different values of M E and different velocities.
From our Figs. 9 and 10 we observe that the scattered wave can be partially polarized in the direction orthogonal to the scattering plane. This phenomenon is similar to the Mott polarization [43], which appears in the electromagnetic scattering. This conclusion was also underlined in Ref. [26].

Absorption cross section
The above examples show that our analytical approach reproduces correctly all the results concerning the elastic scattering obtained by using analytical-numerical methods [25,26].
However, there are significant differences in what concerns the partial absorption cross sections (69) or the total one (68).
We specify first that in our case the partial waves |S −κ | = |S κ | give the same contribution to the partial absorption cross sections. As observed in Sect. III.B., when √ 3m M < 1 all these partial sections remain finite in E = m as we can see in the left side of Fig. 11. However, if this quantity becomes larger than 1, then the partial sections σ l a ( p) with 1 < l < √ 3m M reach the singularity point E = m where these are divergent, as in the right side of Fig. 11, where only σ 1 a ( p) becomes singular. Thus the profile of σ a ( p) is strongly dependent on the parameter m M, as shown in Fig. 12. In the first panel (Fig. 12) we plot the cross section σ a ( p) = σ 1 a ( p) given by |S −1 | = |S 1 |, which are the only contributions in the energy range under consideration. The next panels show the profiles of the absorption cross sections for different values of m M, pointing out the mentioned effect of the singularity in E = m and the asymptotic behavior resulting from Eq. (71).
We note that the analytical-numerical results obtained so far [25,26] are somewhat different from the analytical ones we present here. The first difference is that one obtains numerically |S −κ | = |S κ | such that the quantity is maximal for |κ| = 1-when σ −1 a (with l = 0) diverges for E → m while σ 1 a ( p) remains finitevanishing then rapidly with increasing |κ| [25]. The second difference is the numerical asymptote that corresponds to the geometrical photon value 27π M 2 , while in our case this asymptote (71) is much smaller. On the other hand, despite these two discrepancies, there are notable similarities as the general profile of the cross sections and the values of the thresholds of the partial waves become equidistant for m = 0.
Under such circumstances it seems that something is missing in our analytic approach, namely a possible resonant scattering that may be dominant in the s wave (with l = 0) and increases the absorption cross section significantly; this could explain satisfactorily the differences discussed above.

Concluding remarks
In this paper we performed the partial wave analysis of Dirac fermions scattered by black holes by using the approximative analytical scattering modes of the Dirac equation in Schwarzschild's charts with Cartesian gauge. Selecting a suitable asymptotic condition we obtained the analytical expression of the phase shifts that allow us to derive the scattering amplitudes, differential cross sections, and degrees of polarization. These depend on a parameter (s) that takes real values in elastic collisions, becoming pure imaginary when the fermion is absorbed by a black hole. Thus by using the same formalism we can study the elastic scattering and absorption deriving the closed form of the absorption cross section produced by the scattering modes.
From our analytical and graphical results we established that the spinor wave can be scattered in forward and backward directions and that the spiral scattering (orbital scattering) is present in both these cases. The scattering intensity in both forward and backward directions is increasing with the black hole mass. Also the oscillations in the scattering intensity around θ = π become more pronounced as the black hole mass is increasing. The polarization degree has an oscillatory behavior in terms of the scattering angle and depends on the black hole mass too. Our polar plots for the polarizations show what the directions are in which the spin is aligned after the interaction. We note that thanks to our analytical formulas the graphical analysis can be performed for any values of the parameters m M, M E, and v, leading to similar conclusions.
The general conclusion is that our formalism is suitable for describing the elastic fermion-black hole collisions for which we recover the entire phenomenology pointed out by analytical-numerical methods [22][23][24][25][26]. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 . , with arbitrary integration constants, which holds for any value of the parameter s, which can take either real values, s = |s|, or pure imaginary ones, s = ±i|s|. We observe that the phase shifts depend now only on the two real integration constants C and the relative phase θ = θ 1 − θ 2 . The elastic scattering can arise only when we have e 2iδ κ = 1.
There are two cases. In the first one, when s = |s|, Eq. (C3) has two real solutions, C = 0 and The second case is of pure imaginary s = ±i|s| when the above equation has no real solutions. Furthermore, we observe that the phase shifts have the correct Newtonian limits (B1) for large l only if we choose the asymptotic condition C = 0 (i.e. C + 2 = C − 2 = 0) when the phase shifts are completely determined, being given by Eq. (38). Otherwise, if we consider the solution (C4) we obtain non-determinate phase shifts, which are still depending on the arbitrary phase θ . Obviously, in this case we cannot speak about the Newtonian limit. The conclusion is that we must consider the asymptotic condition C = 0 giving the correct phase shifts in elastic collisions for s = |s|. However, it is natural to keep the same condition for s = ±i|s| when the collision is no longer elastic because of the absorption of the fermions by a black hole.