Loading [MathJax]/jax/output/HTML-CSS/jax.js

Fermion scattering by a class of Bardeen black holes

  • In this study, the scattering of fermions by a class of Bardeen black holes is investigated. After obtaining the scattering modes by solving the Dirac equation in this geometry, we use the partial wave method to derive an analytical expression for the phase shifts that enter into the definitions of partial amplitudes that define the scattering cross sections and induced polarization. It is shown that, similar to Schwarzschild and Reissner-Nordström black holes, the phenomena of glory and spiral scattering are present.
  • The existence of black holes is theoretically well motivated in general relativity (GR) and other modified theories of gravity. Thus far, black holes have not been directly detected or observed. However, plenty of indirect evidence indicate that such objects exist in nature. As shown by Penrose and Hawking [1], the occurrence of singularities is inevitable in GR. This raises many issues such as the black hole information paradox [2-4] or the more recently discovered paradox of black hole's "firewalls" [5-9]. These paradoxes emerge because of the incompatibility between the quantum theory and GR. It is widely believed that in the quantum theory of gravity, the singularities contained in black holes will be removed. However, even in the early stages of the investigations on singularities [10, 11] in GR, several black hole models, which could avoid the occurrence of a singularity, have been proposed. These black holes are called "regular," because they are singularity-free.

    The first proposal of a regular black hole solution was made by Bardeen in Ref. [12], and since then many other models of spherically symmetric regular black holes have been presented in literature [13-21]. In Ref. [22] the authors showed that the Bardeen black hole model can be physically interpreted as the gravitational field produced by a nonlinear magnetic monopole. Later, this interpretation was extended to also include nonlinear electric charges so that regular black holes models can have a nonlinear electromagnetic field as a source. More recently, in Ref. [23] the authors extended the Bardeen solution to an entire class of Bardeen-like black holes that can or cannot be regular.

    In this study, the scattering of fermions (spin 1/2) by a Bardeen regular black hole and by a Bardeen-class of black holes (as constructed in [23]) is investigated. We use the partial wave method to obtain analytical expressions for the phase shifts that enter into the definition of partial amplitudes, defining the scattering cross sections and the induced polarization. To the best of our knowledge, this is the first study that reports analytical phase shifts for 1/2-spin wave scattering by regular black holes. The absorption of fermions by Bardeen black holes was numerically investigated in [24], while in Refs. [25-28], the case of massless scalar scattering by regular black holes was also treated numerically. Studies dedicated to fermion scattering by other types of spherically symmetric black holes can be found as examples in Refs. [29-43].

    Both Bardeen regular black holes and the Bardeen-class type of black holes possess nonlinear magnetic (monopole) charges [22, 23]. This implies the existence of an electromagnetic potential of the form A=Qmcosθdϕ , with Qm the total magnetic charge. In this work, we neglect the interaction between this potential and the charge of the fermion and only focus on the scattering due to the "pure" gravitational interaction between the fermion and the black hole. However, even in this approximation, the black hole magnetic monopole charge will still influence the scattering patterns through the presence of Qm into the metric function and into the resulting scattering modes of the radial Dirac equation.

    The remainder of the paper is structured as follows. In Section 2, the Bardeen-class of black holes is presented very briefly. Section 3 begins with a very short review of the Dirac equation in spherically symmetric black hole geometries, and then it discusses the search for scattering modes in the Bardeen spacetime. Finally, the main result of the paper is presented, namely the form of the analytical phase shifts resultant from applying the partial wave method on the scattering modes derived earlier. Section 4 presents a graphical analysis of the induced polarization and the scattering cross sections in which the presence of a backward "glory" and "spiral scattering" (orbiting) oscillations are shown to be present. The main conclusions and some final remarks are presented in Section 5.

    In Ref. [23], a class of spherically symmetric and asymptotically flat Bardenn-like black holes depending on two-parameters was constructed with the following line element

    ds2=h(r)dt2dr2h(r)r2(dθ2+sin2θdϕ2),h(r)=12Msr2Mrγ1(r2+Q2)γ.

    (1)

    These black holes have a singularity if the parameter Ms , hereinafter, referred to as the "Schwarzschild mass," has a nonzero value. Otherwise, if Ms=0 , regular black holes are obtained and the Bardeen black hole corresponds to the particular choice Ms=0 and γ=3 . The term M can be interpreted as the mass of the nonlinear magnetic monopole. Moreover, the sum

    MADM=Ms+M

    (2)

    constitutes the ADM mass of the black hole obtained from the asymptotic form of the metric function h(r). The parameter Q in Eq. (1) is related to the magnetic monopole charge Qm ; Q=2Q2m/M .

    As shown in Ref. [23], the Lagrangian density, for which Eq. (1) is a solution of the coupled Einstein-Maxwell filed equations, is given by

    L=4γα(αFμνFμν)5/4(1+αFμνFμν)1+γ/2,

    (3)

    where α has the dimensions of length squared, and γ is a dimensionless constant. In the weak field limit, a vector field, which is slightly stronger when compared with a Maxwell field, is obtained.

    The Dirac equation:

    iγaDaψmψ=0

    (4)

    The following explicit form can be obtained from the Dirac equation [44]

    (iγaeμaμm)ψ+i21gμ(geμa)γaψ14{γa,Sbc}ωcabψ=0,

    (5)

    where g=det(gμν) and the covariant derivative is defined by Da=a+i2Sbcωcab . Sab=i4[γa,γb] the generators of the SL(2,C) group, γa are the point-independent Dirac matrices obeying {γa,γb}=2ηab and ωcab is a spin-connection

    ωcab=eμaeνb(ˆecλΓλμνˆecν,μ)

    (6)

    Γλμν represents the usual Christoffel symbols. The tetrad fields ea(x) and ˆea(x) are point-dependent defining (nonholonomic) local frames and co-frames. The following relations hold:

    ea=eμaμ,ˆea=ˆeaμdxμea(x)eb(x)=ηab,ˆea(x)ˆeb(x)=ηabds2=ηabˆeaμdxμˆebνdxν=gμν(x)dxμdxν.

    (7)

    The so-called Cartesian gauge [44-46] for a spherically symmetric line element of the form of Eq. (1) is defined by the following tetrad fields:

    ˆe0=hdtˆe1=1hsinθcosϕdr+rcosθcosϕdθrsinθsinϕdϕˆe2=1hsinθsinϕdr+rcosθsinϕdθ+rsinθcosϕdϕˆe3=1hcosθdrrsinθdθ,

    (8)

    The Dirac equation (5) can be reduced to only a radial equation. The angular part of the Dirac equation is the same as in the Dirac theory from flat spacetime, and its solutions are the usual four-component angular spinors Φ±m,κ(θ,ϕ) [47, 48]. This is because, in the Cartesian gauge (8), the Dirac equation is manifestly covariant under rotations [46]. Using this gauge, it was possible to find complete analytical solutions to the Dirac equation on de Sitter/anti-de Sitter spacetime [49-52] and approximate analytical solutions in black hole geometries [31-33, 53, 54], which were later used to study different aspects of the scattering problem in those spacetimes [31-33, 55-58].

    The remaining unsolved radial part of the Dirac equation was obtained by assuming the following type of particle-like solutions with a given energy E

    ψ(x)=ψE,j,m,κ(t,r,θ,ϕ)=eiEtrh(r)1/4{f+E,κ(r)Φ+m,κ(θ,ϕ)+fE,κ(r)Φm,κ(θ,ϕ)}

    (9)

    where f±E,κ(r) are two unknown radial wave functions. As in Refs. [31, 53] the radial Dirac equation can be expressed in a matrix form

    (mh(r)h(r)ddr+κrh(r)h(r)ddr+κrh(r)mh(r))×(f+E,κ(r)fE,κ(r))=E(f+E,κ(r)fE,κ(r)).

    (10)

    Substituting the line element (1) in Eq. (10), we obtain a system of two differential equations that has no analytical solutions owing to the complex form of the line element. However, because we are only interested in finding the scattering modes, we can approximate Eq. (10) in the asymptotic region of the Bardeen black hole specified by the line element in Eq. (1) and find approximative analytical solutions. A partial wave method was used on the obtained solutions to compute the (elastic) scattering cross section and the induced polarization that results after the interaction of a fermion beam with the black hole.

    Let us start by introducing the new variable

    x=zr+1,z=r2+Q2,

    (11)

    where r+ is the radius of the black hole horizon. The function h(r), in terms of x, is

    h=12Msr+1(1+x2)2δ2Mr+11+x2[1δ(1+x2)2]γ12

    (12)

    where the notation δ=(Q/r+)2 was introduced. The radial Dirac equation (10) is a system of two differential equations for the radial wave functions f±(x) that, after multiplying each equation by x/[r+(1+x2)] and expressing all the terms as a function of x given by Eq. (11), is equivalent to Eq. (13)

    [11+x21δ(1+x2)2{...}12ddx+κ(1+x2)2δ×{...}12]f+E,κ(x)[μ{...}12+ϵ(x+1x)]fE,κ(x)=0[11+x21δ(1+x2)2{...}12ddxκ(1+x2)2δ×{...}12]fE,κ(x)[μ{...}12ϵ(x+1x)]f+E,κ(x)=0where{...}=(1+x2)2x22Msr+(1+x2)2x2(1+x2)2δ2Mr+1+x2x2[1δ(1+x2)2]γ12,μ=r+m,ϵ=r+E.

    (13)

    The system of differential equations obtained in Eq. (13) cannot directly be solved analytically. However, an analytical solution can be obtained if a domain is restricted far away from the black hole event horizon. A simpler system of differential equations that has analytical solutions can be obtained. By Taylor expansion with respect to 1/x and discarding the terms of the order O(1/x2) and higher, the system of equations, which is valid in the asymptotic region of the black hole for the two radial wave function, reduces to

    (12ddx+κx)f+E,κ(x)x(ϵ+μ)fE,κ(x)1x[ϵ+μ(1MADMr+)]fE,κ(x)=0(12ddxκx)fE,κ(x)+x(ϵμ)f+E,κ(x)+1x[ϵμ(1MADMr+)]f+E,κ(x)=0.

    (14)

    We observe that if M0 then the ratio MADMr+Msrs=12 and μ(1MADMr+)12μ , recovering the Schwarzschild case discussed in Ref. [53].

    To find the scattering modes for the case ϵ>μ , it is useful to introduce new radial wave functions ˆf± , present in the following combination of the old wave functions f±

    ˆf+E,κ=i2f+E,κϵ+μ+12fE,κϵμˆfE,κ=i2f+E,κϵ+μ+12fE,κϵμ.

    (15)

    Using the above relation we can obtain the following equations, which are satisfied by the functions ˆf+E,κ and ˆfE,κ :

    [νxddx+2i(μ2(1MADMr+)ϵ2ν2x2)]ˆf+E,κ(2κνiϵμMADMr+)ˆfE,κ=0[νxddx2i(μ2(1MADMr+)ϵ2ν2x2)]ˆfE,κ(2κν+iϵμMADMr+)ˆf+E,κ=0,

    (16)

    where ν=ϵ2μ2 . Eq. (16) can now be solved using Maple or Mathematica and the analytical solutions can be written as a combination of Whittaker functions

    ˆf+E,κ=C11xMρ+,s(z)+C21xWρ+,s(z)ˆfE,κ=1κ2+λ2[(siα)(κ+iλ)C11xMρ,s(z)(κ+iλ)C21xWρ,s(z)],

    (17)

    where z=2iνx2 and the following parameters were introduced

    s=[κ2+μ2(1MADMr+)2ϵ2]12,λ=ϵμνMADMr+ρ±=12iα,α=1ν[ϵ2μ2(1MADMr+)].

    (18)

    3.2.1 Pure Bardeen black hole case

    As already mentioned, by choosing Ms=0 and γ=3 in Eq. (1) the original Bardeen black hole solution is recovered. This solution is a black hole with two distinct horizons only if |Q|<4/27 , it has degenerate horisons if Q=4/27 and for Q>4/27 there are no horizons present. Because the function h(r) has now a more simpler form,

    h(r)=12Mr2(r2+Q2)3,

    (19)

    one can find and write an analytical expression for the location of the black hole outer horizon, denoted from now on by r+ . Moreover, one can easily show that r+=Mf(Q/M) , where the function f depends only on the ratio Q/M. The existence of horizons implies the constraint |Q|/M4/27 [14].

    Following the same steps as in section 3.2, one finds the same scattering modes as given by Eq. (17), but with the new parameters

    s=[κ2+μ2(1Mr+)2ϵ2]12,λ=ϵμνMr+ρ±=12iα,α=1ν[ϵ2μ2(1Mr+)].

    (20)

    We will see in section 4 that this is enough to produce a noticeable difference in the scattering patterns.

    By choosing Ms=0 and γ=3 in Eq. (1), the original Bardeen black hole solution can be recovered. This solution is a black hole with two distinct horizons only if |Q|<4/27 ; it has degenerate horizons if Q=4/27 , and for Q>4/27 there are no horizons present. Because the function h(r) can now be expressed in a simpler form:

    In spinor wave scattering theory [48, 59] the differential scattering cross section for an unpolarized incident beam is the sum of the squares of two scalar functions

    dσdΩ=|f(θ)|2+|g(θ)|2,

    (21)

    which only depend on the scattering angle θ :

    f(θ)=l=012ip[(l+1)(e2iδl11)+l(e2iδl1)]P0l(cosθ),g(θ)=l=112ip[e2iδl1e2iδl]P1l(cosθ),

    (22)

    where p is the incident momentum, and P0l(cosθ),P1l(cosθ) are the Legendre and associated Legendre polynomials, which are special cases of Legendre functions, respectively [60].

    The phase shifts δl can be computed as in Refs. [31, 32] by applying the partial wave method on the scattering modes (17). The asymptotic form of the radial wave functions f± can be written as [31]

    (f+E,κfE,κ)E+msinEmcos(prπl2+δκ+ϑ(r))

    (23)

    where ϑ(r)=pr++αln[2p(rr+)] represents a radially dependent phase that is independent of any angular quantum numbers; therefore, it does not contribute to the scattering cross sections and may be neglected as in the Dirac-Coulomb case [30, 48].

    The resultant final form for the point-independent phase shifts δκ is given by the following expression (see also Appendix A)

    e2iδκ=κiλsiαΓ(1+siα)Γ(1+s+iα)eiπ(ls),

    (24)

    where the sign convention for κ is the same as in [48], such that κ=±(j+1/2) and l=|κ|(1signκ)/2 .

    Series (22) are poorly convergent as a direct consequence of the singularity present at θ=0 ; an infinite number of Legendre polynomials are necessary to describe it. To make the series more convergent, the mth reduced series can be defined

    (1cosθ)m1f(θ)=l0a(m1)lPl(cosθ),(1cosθ)m2g(θ)=l1b(m2)lP1l(cosθ),

    (25)

    as first proposed in [61] and more recently used in [30-32]. The coefficients a(i)l and b(i)l are computed using the recurrence relations

    a(i+1)l=a(i)ll+12l+3a(i)l+1l2l1a(i)l1,b(i+1)l=b(i)ll+22l+3b(i)l+1l12l1b(i)l1,

    (26)

    where a(0)l=[(l+1)(e2iδl11)+l(e2iδl1)]/2ip and b(0)l=[e2iδl1e2iδl]/2ip from Eq. (22). We found that using only two iterations m1=2 for the function f(θ) and one iteration m2=1 for g(θ) , it is sufficient to make the series converge without distorting the analytical results.

    In this section, we present and discuss the main features of fermion scattering by Bardeen regular black holes and also by a Bardeen-class of black holes. The analysis will focus on scattering by small or micro black holes (with MADM10151022 Kg) because in this case the glory and orbiting scattering phenomena are shown to be significant.

    The following parameters are used to label the figures: v=p/E is the speed of the incident fermions; EM can be considered a dimensionless measure of the gravitational coupling because (restoring the units) it forms the dimensionless quantity: ε=GEMc3=πrsλv , with rs the Schwazschild radius and λ=h/p the associated quantum particle wavelength; and the ratios q=Q/Ms and g=M/Ms that appear when writing the black hole horizon radius as r+=Msf(MMs,QMs) . Moreover, because in the asymptotic zone E=m2+p2 , one can easily show that mM=EM1v2 such that the condition EMmM is always satisfied. In the following analysis we take γ=3 in all the plots (except those in Fig. 7) because the same conclusions are also obtained for the cases with γ3 .

    Figure 7

    Figure 7.  (color online) Comparison of the scattering cross section for a Bardeen-class black hole with different values of the parameter γ . Note that the scattering profile keeps the same shape for different values of γ .

    In Fig. 1, the differential scattering cross section is presented as a function of the scattering angle for fixed values of the ratios: Q/M (left panels corresponding to pure Bardeen case); q=Q/Ms and g=M/Ms ; a fixed value of the speed (v) of incoming fermions, while the parameter EM (or EMs for Bardenn-class) takes different values. From Fig. 1, it can be seen that the scattering pattern takes a simple form for small values of EM. However, as the value of EM increases, more complex scattering patterns occur. This includes the presence of a maximum in the backward direction ( θ=π ), which is also known as "glory" scattering [62, 63] and the presence of oscillations in the scattering intensity that give rise to orbiting or "spiral scattering" [64, 65] (that may occur when the particle's "classical" orbit passes the scattering center multiple times). As the speed of the incoming fermion increases, the peak in the π -direction moves toward the left, and the maximum of the scattering intensity that occurs at θ=π for nonrelativistic fermions transforms into a minimum if the fermions are massless (v= 1). From the bottom panels of Fig. 1, it can be seen that the magnitude of the spiral scattering oscillations and their angular frequency increase with the black hole mass.

    Figure 1

    Figure 1.  (color online) Plot of the differential scattering cross section as a function of the scattering angle for a regular Bardeen black hole, Eq. (19) (left panels) and for a Bardeen-class black bole, Eq. (1) (right panels). Phenomena of glory (scattering in the backward direction) and spiral scattering (oscillations in the scattering intensity) are present for both types of black holes. The value q/Mem=4/27 corresponds to a degenerate regular Bardeen black hole (it has only one horizon). The parameter q gives the ratio q=Q/Ms and g=M/Ms is the ratio between the magnetic monopole mass and the "Schwarzschild mass".

    The scattering cross section dσ/dΩ for the Bardeen-class black holes (1) is plotted in Fig. 2. The left panels show the scattering patterns for Bardeen-class black holes that have the "Schwarzschild mass" greater than the mass of the magnetic monopoles, while the right side panels show those with a lower mass. By comparing our analysis results with the Schwarzschild black hole, the scattering intensity in the backward direction ( θ=π ) is higher for the scattering by a Bardeen-class black hole. From Fig. 2, one can also observe that in all the plots, the value of g increases as the number of oscillations of the scattering intensity increases. Furthermore, if we assume Ms is fixed, then the ADM-mass of the black hole increases (which is equivalent to increasing the mass M, Eq. (2), if Ms remains constant), and as a consequence the angular frequency of the oscillations present in the scattering intensity also increase. This feature was also true for fermion scattering by Schwarzschild and Reissner-Nordström black holes [30-33]. Thus, the increase in the frequency of the oscillations of the spiral scattering with the black hole mass for any type of (spherically symmetric) black holes is of interest.

    Figure 2

    Figure 2.  (color online) Variation of the scattering cross section with the parameter g=M/Ms for given values of EMs and v. In the top panels, the Schwarzschild scattering (for which g= 0) is compared. The value of g increases as the spiral scattering becomes more pronounced.

    Comparing the scattering intensity, in Fig. 3, of a Bardeen regular black hole (blue and red curves) with that of a Schwarzschild black hole (the black dotted curve), we observe that the glory peak is higher for the scattering by a Schwarzschild black hole. As the ratio Q/M increases, the maximum in the backward direction becomes lower, and the frequency of the oscillations in the spiral scattering slightly decrease simultaneously. The curve with Q/M=4/27 corresponds to the Bardeen degenerate case when only one horizon is present. The type of black hole can be determined by analyzing the patterns in the scattering intensity when the scattering of a beam of fermions by a black hole can be observed and measured. If it is a Bardeen black hole, with a known mass, then the value of Q can also be found.

    Figure 3

    Figure 3.  (color online) Scattering cross section for a regular Bardeen black hole for different values of Q/M, while the parameters EM and v are kept fixed.

    In Fig. 4, we plotted the differential scattering cross section dσdΩ as a function of the scattering angle θ for a Schwarzschild black hole (blue solid lines), regular Bardeen black hole (black dotted lines) and Bardeen-class black hole (red dash-dotted lines).

    Figure 4

    Figure 4.  (color online) Comparison between the fermion scattering by a Schwarzschild black hole, regular Bardeen black hole, and Bardeen-class black hole. For the left panel Q/M= 0.2, q= 0.2, and g= 0.5 were used, while for the right panel Q/M= 0.4, q= 0.1, and g= 0.2.

    An incident unpolarized beam of massive fermions becomes partially polarized after it gets scattered by the black hole. The induced polarization degree can be computed using the following formula [30, 59]

    P=ifgfg|f|2+|g|2n

    (27)

    where n is a unit vector in a direction orthogonal to the plane of scattering.

    In Fig. 5, the dependence of the polarization on the scattering angle is plotted for given values of the parameters EM, v, g, and Q/M. It has a pronounced oscillatory behavior. From the top panels in Fig. 5, one can observe that if the parameter EM increases, then the frequency of the oscillations present in the polarization increases as well. Now, if we assume that the energy of the fermion is fixed, then, the oscillations present in the polarization are more pronounced for black holes with higher masses. Compared with the Schwarzschild polarization (see bottom panels in Fig. 5) the Bardeen and Bardeen-class black hole polarizations exhibit a slightly less oscillatory behavior. The oscillations present in the polarization can be observed as a consequence of the oscillations present in the glory and spiral scattering.

    Figure 5

    Figure 5.  (color online) Partial polarization P as a function of the scattering angle θ . Top left panel: typical Bardeen black hole (19) with Q/M= 0.4 for incoming fermions with speed v= 0.8 in units of c; Top right panel: Bardeen-class black hole (1) with g=M/Ms=0.2 for incoming fermions with speed v=0.4 in units of c; Bottom left panel: comparing Schwarzschild and pure Bardeen black hole polarizations for fixed EM= 1 and v= 0.3; Bottom right panel: comparing Schwarzschild and Bardeen-class black hole polarizations for fixed EMs=1.5 , q=Q/Ms=0.2 and v= 0.5.

    The alignment of the scattered fermions with the forward on-axis direction can be visualized using the polar plots representations of the polarization degree as shown in Fig. 6. One can observe the Mott polarization in the direction orthogonal to the scattering plane, a phenomena that has been reported before in literature for Schwarzschild [30, 31] and Reissner-Nordström [32] black holes.

    Figure 6

    Figure 6.  (color online) Polar plots of pure Bardeen (left panel) and Bardeen-class (right panel) polarization P(θ) for 0θ2π showing the alignment of the scattered fermion's spin with a given direction. The Mott polarization (polarization in the direction orthogonal to the scattering plane) can also be observed.

    In Fig. 7, the differential scattering cross section is plotted for a Bardeen-class black hole with different values of the parameter γ . Note that the scattering pattern maintains the same profile shape, and the difference between the values of dσ/dΩ (computed using different values of γ ) is very small. Moreover, from our analysis results, it can be seen that, for most cases, this difference is even smaller than for the examples presented in Fig. 7. This can be attributed to the difference between the black hole horizon radii r+ obtained for different values of γ (while all the other parameters are kept constant) is very small.

    We studied the scattering of fermions by a class of Bardeen black holes that also include the original Bardeen regular black hole solution. A partial wave method was used on a set of scattering modes that were obtained by solving the Dirac equation in the asymptotic region of these black hole geometries. TO the best of our knowledge, for the first time, time analytical phase shifts for fermion scattering by Bardeen regular black holes were obtained. The phenomena of glory (scattering in the backward direction) and spiral scattering (oscillations in the scattering intensity) were also discussed. We also saw that an incident unpolarized beam could become partially polarized after interaction with a black hole.

    In Figs. 1-6, other than the parameters EM (that can be associated with a measure of the gravitational coupling) and v (speed of the fermion) we also used the ratios Q/M, q=Q/Ms , and g=M/Ms to label the figures. The departure of the scattering pattern from the Schwarzschild case was significant as Q/M and g increased. For the original Bardeen regular black hole, Q/M=4/27 was the maximum allowed value, and it corresponds to the degenerate case when the two horizons coincide. If the magnetic charge Q0 , then M0 and g0 , and the scattering by a Schwarzschild black hole can be recovered [31].

    The glory and spiral scattering become significant for values of the parameter EM~1 or greater (in geometrical units with G==c=1 ) because the associated wavelength for these values of the incident fermions is of the same order of magnitude as the black hole horizon radius; thus, diffraction patterns occur. Another feature, which was also shown to be present for fermion scattering by Schwarzschild [30, 31] and Reissner-Nordström black holes [32], is that as the total mass of the black hole increases, the oscillations present in the scattering intensity become more frequent, that is, spiral scattering is significantly enhanced with an increase in black hole mass.

    The aim of this Appendix is to briefly present how the phase shifts (24) are obtained.

    The asymptotic representation of the Whittaker function M for large values of |z| is given by the following equation [60]

    Mκ,μ(z)Γ(1+2μ)Γ(12+μκ)e12zzκ(1+O(z1))+ei(12+μκ)πΓ(1+2μ)Γ(12+μ+κ)e12zzκ(1+O(z1))

    valid for 12π<phz<32π .

    To obtain the phase shifts (24), the condition C2=0 must be imposed on Eq. (17) in order to obtain the correct Newtonian phase shifts for large values of angular momenta. As an observation, this condition also selects the spinors that are regular in x= 0 owing to the regularity of the function Mρ±,s(2iνx2)=(2iνx2)s+12[1+O(x2)] . More arguments for letting C2=0 can be found in Appendix C of our previous paper [31].

    Using Eq. (28) on the ˆf±(x) functions and by observing that for ˆf+ the dominant term is the first one in (28), while for ˆf the second term is dominant, then one obtains the following asymptotic expressions

    ˆf+(x)C1e12παΓ(2s+1)Γ(1+s+iα)ei[νx2+αln(2νx2)]ˆf(x)(siα)(κ+iλ)κ2+λ2C1e12παΓ(2s+1)Γ(1+siα)×ei[νx2+πs+αln(2νx2)].

    The last step is to compute the argument of the trigonometric functions (23) as 12arg(ˆf+ˆf) from which the expression for e2iδκ will result.

    [1] S.W. Hawking and G. F. R. Ellis, The Large Scale Structure of Spacetime (Cambridge: Cambridge University Press, 1973)
    [2] Stephen Hawking, Phys. Rev. D, 14: 2460 (1976) doi: 10.1103/PhysRevD.14.2460
    [3] G. T. Horowitz and J. M. Maldacena, JHEP, 0402: 008 (2004)
    [4] S. D. Mathur, Lect. Notes Phys., 769: 3 (2009) doi: 10.1007/978-3-540-88460-6
    [5] H. Zhang, Sci. Rep., 7: 4000 (2017) doi: 10.1038/s41598-017-03854-y
    [6] A. Almheiri, D. Marolf, J. Polchinski, and J. Sully, JHEP, 1302: 062 (2013)
    [7] A. Almheiri, D. Marolf, J. Polchinski et al, JHEP, 1309: 018 (2013)
    [8] P. Chen, Y. Chin Ong, D. N. Page et al, Phys. Rev. Lett., 116: 161304 (2016) doi: 10.1103/PhysRevLett.116.161304
    [9] J. Maldacena and L. Susskind, Fortsch. Phys., 61: 781 (2013) doi: 10.1002/prop.v61.9
    [10] A. D. Sakharov, Zh. Eksp. Teor. Fiz. 49, 345 (1966) [Sov.Phys. JETP 22, 241 (1966)].
    [11] E. B. Gliner, Sov. Phys. JETP, 22: 378 (1966)
    [12] J. M. Bardeen, in: Conference Proceedings of GR5, Tbilisi, USSR, 1968, p. 174
    [13] A. Borde, Phys. Rev. D, 50: 3692 (1994) doi: 10.1103/PhysRevD.50.3692
    [14] A. Borde, Phys. Rev. D, 55: 7615 (1997)
    [15] C. Barrabes and V. P. Frolov, Phys. Rev. D, 53: 3215 (1996)
    [16] A. Cabo and E. Ayon-Beato, Int. J. Geom. Methods Mod. Phys. A, 14: 2013 (1999) doi: 10.1142/S0217751X99001019
    [17] S. A. Hayward, Phys. Rev. Lett., 96: 031103 (2006) doi: 10.1103/PhysRevLett.96.031103
    [18] E. Ayon-Beato and A. Garcia, Phys. Rev. Lett., 80: 5056 (1998) doi: 10.1103/PhysRevLett.80.5056
    [19] E. Ayon-Beato and A. Garcia, Phys. Lett. B, 464: 25 (1999) doi: 10.1016/S0370-2693(99)01038-2
    [20] E. Ayon-Beato and A. Garcia, Gen. Relativ. Gravit., 37: 635 (2005) doi: 10.1007/s10714-005-0050-y
    [21] M. S. Ma, Ann. Phys. (Amsterdam), 362: 529 (2015) doi: 10.1016/j.aop.2015.08.028
    [22] E. Ayon-Beato and A. Garcia, Phys. Lett B, 493: 149 (2000) doi: 10.1016/S0370-2693(00)01125-4
    [23] Z.-Y. Fan and X. Wang, Phys. Rev. D, 94: 124027 (2016) doi: 10.1103/PhysRevD.94.124027
    [24] H. Huang, M. Jiang, J. Chen, and Y. Wang, Gen. Relativ. Gravit., 47: 8 (2015) doi: 10.1007/s10714-014-1850-8
    [25] H. Huang, P. Liao, J. Chen, and Y. Wang, J. Grav. 2014,231727 (2014)
    [26] C.F.B. Macedo and L.C.B. Crispino, Phys. Rev. D, 90: 064001 (2014) doi: 10.1103/PhysRevD.90.064001
    [27] C.F.B. Macedo, L.C.B. Crispino, and E.S. de Oliveira, Int. J. Mod. Phys. D, 25: 1641008 (2016)
    [28] C.F.B. Macedo, E.S. de Oliveira, and L.C.B. Crispino, Phys. Rev. D, 92: 024012 (2015) doi: 10.1103/PhysRevD.92.024012
    [29] W. G. Unruh, Phys. Rev. D, 14: 3251 (1976)
    [30] S. Dolan, C. Doran, and A. Lasenby, Phys. Rev. D, 74: 064005 (2006) doi: 10.1103/PhysRevD.74.064005
    [31] I. I. Cotăescu, C. Crucean, and C. A. Sporea, Eur. Phys. J. C, 76: 102 (2016) doi: 10.1140/epjc/s10052-016-3936-9
    [32] I. I. Cotăescu, C. Crucean, and C. A. Sporea, Eur. Phys. J. C, 76: 413 (2016) doi: 10.1140/epjc/s10052-016-4260-0
    [33] C. A. Sporea, Chinese Physics C, 41(12): 123101 (2017) doi: 10.1088/1674-1137/41/12/123101
    [34] S. R. Das, G. Gibbons, and S. D. Mathur, Phys. Rev. Lett., 78: 417 (1997) doi: 10.1103/PhysRevLett.78.417
    [35] W. M. Jin, Classical and Quantum Gravity, 15(10): 3163 (1998) doi: 10.1088/0264-9381/15/10/018
    [36] C. Doran and A. Lasenby, Phys. Rev, D, 66: 024006 (2002) doi: 10.1103/PhysRevD.66.024006
    [37] E. Jung, S. H. Kim, and D. K. Park, JHEP, 09: 005 (2004)
    [38] H. Cho and Y. Lin, Classical and Quantum Gravity, 22(5): 775 (2005) doi: 10.1088/0264-9381/22/5/001
    [39] A. B. Gaina and G. A. Chizhov, Moscow Univ. Phys. Bull., 38N2: 1-7 (1983)
    [40] A. B. Gaina, Moscow VINITI, No. 1970-80 Dep. , 20 pp. (1980)
    [41] M. Rogatko and A. Szyplowska, Phys. Rev. D, 79: 104005 (2009) doi: 10.1103/PhysRevD.79.104005
    [42] H. Liao, J-H. Chen, P. Liao, and Y-J. Wang, Commun. Theor. Phys., 62: 227-234 (2014)
    [43] A. Ghosh and P. Mitra, Phys. Rev. D, 50: 7389-7393 (1994)
    [44] I.I. Cotăescu, Mod. Phys. Lett. A, 13: 2923 (1998) doi: 10.1142/S0217732398003107
    [45] V. M. Villalba, Mod. Phys. Lett. A, 8: 2351-2364 (1993)
    [46] I. I. Cotăescu, J. Phys. A Math. Gen., 33: 1977 (2000)
    [47] B. Thaller. The Dirac Equation. (Springer Verlag, Berlin Heidelberg, 1992)
    [48] V. B. Berestetski, E. M. Lifshitz, and L. P. Pitaevski, Quantum Electrodynamics, (Oxford: Pergamon Press, 1982)
    [49] I. I. Cotăescu, Phys. Rev. D, 60: 124006 (1999) doi: 10.1103/PhysRevD.60.124006
    [50] I. I. Cotăescu, Phys. Rev. D, 65: 084008 (2002) doi: 10.1103/PhysRevD.65.084008
    [51] I. I. Cotaescu and C. Crucean, Mod. Phys. Lett. A, 22: 3707-3720 (2008)
    [52] I. I. Cotaescu, R. Racoceanu, and C. Crucean, Mod. Phys. Lett. A, 21: 1313-1318 (2006) doi: 10.1142/S0217732306020044
    [53] I.I. Cotăescu, Mod. Phys. Lett. A, 22: 2493 (2007)
    [54] C.A. Sporea, Mod. Phys. Lett. A, 30: 1550145 (2015)
    [55] C. Crucean, Mod. Phys. Lett. A, 22: 2573-2585 (2007) doi: 10.1142/S0217732307025480
    [56] C. Crucean, Mod. Phys. Lett. A, 25: 1679-1687 (2010) doi: 10.1142/S0217732310033414
    [57] C. Crucean, R. Racoceanu, and A. Pop, Phys. Lett. B, 665: 409-411 (2008) doi: 10.1016/j.physletb.2008.06.034
    [58] I. I. Cotaescu and C. Crucean, Int. J. Mod. Phys. A, 23: 1351-1359 (2008) doi: 10.1142/S0217751X08039487
    [59] M. Rose. Relativistic Electron Theory (New York: John Wiley & Sons, 1961)
    [60] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, NIST Handbook of Mathematical Functions (Cambridge: Cambridge University Press, 2010)
    [61] D. R. Yennie, D. G. Ravenhall, and R. N. Wilson, Phys. Rev., 95: 500 (1954) doi: 10.1103/PhysRev.95.500
    [62] K. W. Ford and J. A. Wheeler, Ann. Phys. (N.Y.), 7: 259 (1959)
    [63] K. W. Ford and J. A. Wheeler, Ann. Phys. (N.Y.), 7: 287 (1959)
    [64] P. Anninos, C. DeWitt-Morette, R. A. Matzner, P. Yioutas, and T. R. Zhang, Phys. Rev. D, 46: 4477 (1992) doi: 10.1103/PhysRevD.46.4477
    [65] R. A. Matzner, C. DeWitte-Morette, B. Nelson, and T.-R. Zhang, Phys. Rev. D, 31: 1869 (1985)
  • [1] S.W. Hawking and G. F. R. Ellis, The Large Scale Structure of Spacetime (Cambridge: Cambridge University Press, 1973)
    [2] Stephen Hawking, Phys. Rev. D, 14: 2460 (1976) doi: 10.1103/PhysRevD.14.2460
    [3] G. T. Horowitz and J. M. Maldacena, JHEP, 0402: 008 (2004)
    [4] S. D. Mathur, Lect. Notes Phys., 769: 3 (2009) doi: 10.1007/978-3-540-88460-6
    [5] H. Zhang, Sci. Rep., 7: 4000 (2017) doi: 10.1038/s41598-017-03854-y
    [6] A. Almheiri, D. Marolf, J. Polchinski, and J. Sully, JHEP, 1302: 062 (2013)
    [7] A. Almheiri, D. Marolf, J. Polchinski et al, JHEP, 1309: 018 (2013)
    [8] P. Chen, Y. Chin Ong, D. N. Page et al, Phys. Rev. Lett., 116: 161304 (2016) doi: 10.1103/PhysRevLett.116.161304
    [9] J. Maldacena and L. Susskind, Fortsch. Phys., 61: 781 (2013) doi: 10.1002/prop.v61.9
    [10] A. D. Sakharov, Zh. Eksp. Teor. Fiz. 49, 345 (1966) [Sov.Phys. JETP 22, 241 (1966)].
    [11] E. B. Gliner, Sov. Phys. JETP, 22: 378 (1966)
    [12] J. M. Bardeen, in: Conference Proceedings of GR5, Tbilisi, USSR, 1968, p. 174
    [13] A. Borde, Phys. Rev. D, 50: 3692 (1994) doi: 10.1103/PhysRevD.50.3692
    [14] A. Borde, Phys. Rev. D, 55: 7615 (1997)
    [15] C. Barrabes and V. P. Frolov, Phys. Rev. D, 53: 3215 (1996)
    [16] A. Cabo and E. Ayon-Beato, Int. J. Geom. Methods Mod. Phys. A, 14: 2013 (1999) doi: 10.1142/S0217751X99001019
    [17] S. A. Hayward, Phys. Rev. Lett., 96: 031103 (2006) doi: 10.1103/PhysRevLett.96.031103
    [18] E. Ayon-Beato and A. Garcia, Phys. Rev. Lett., 80: 5056 (1998) doi: 10.1103/PhysRevLett.80.5056
    [19] E. Ayon-Beato and A. Garcia, Phys. Lett. B, 464: 25 (1999) doi: 10.1016/S0370-2693(99)01038-2
    [20] E. Ayon-Beato and A. Garcia, Gen. Relativ. Gravit., 37: 635 (2005) doi: 10.1007/s10714-005-0050-y
    [21] M. S. Ma, Ann. Phys. (Amsterdam), 362: 529 (2015) doi: 10.1016/j.aop.2015.08.028
    [22] E. Ayon-Beato and A. Garcia, Phys. Lett B, 493: 149 (2000) doi: 10.1016/S0370-2693(00)01125-4
    [23] Z.-Y. Fan and X. Wang, Phys. Rev. D, 94: 124027 (2016) doi: 10.1103/PhysRevD.94.124027
    [24] H. Huang, M. Jiang, J. Chen, and Y. Wang, Gen. Relativ. Gravit., 47: 8 (2015) doi: 10.1007/s10714-014-1850-8
    [25] H. Huang, P. Liao, J. Chen, and Y. Wang, J. Grav. 2014,231727 (2014)
    [26] C.F.B. Macedo and L.C.B. Crispino, Phys. Rev. D, 90: 064001 (2014) doi: 10.1103/PhysRevD.90.064001
    [27] C.F.B. Macedo, L.C.B. Crispino, and E.S. de Oliveira, Int. J. Mod. Phys. D, 25: 1641008 (2016)
    [28] C.F.B. Macedo, E.S. de Oliveira, and L.C.B. Crispino, Phys. Rev. D, 92: 024012 (2015) doi: 10.1103/PhysRevD.92.024012
    [29] W. G. Unruh, Phys. Rev. D, 14: 3251 (1976)
    [30] S. Dolan, C. Doran, and A. Lasenby, Phys. Rev. D, 74: 064005 (2006) doi: 10.1103/PhysRevD.74.064005
    [31] I. I. Cotăescu, C. Crucean, and C. A. Sporea, Eur. Phys. J. C, 76: 102 (2016) doi: 10.1140/epjc/s10052-016-3936-9
    [32] I. I. Cotăescu, C. Crucean, and C. A. Sporea, Eur. Phys. J. C, 76: 413 (2016) doi: 10.1140/epjc/s10052-016-4260-0
    [33] C. A. Sporea, Chinese Physics C, 41(12): 123101 (2017) doi: 10.1088/1674-1137/41/12/123101
    [34] S. R. Das, G. Gibbons, and S. D. Mathur, Phys. Rev. Lett., 78: 417 (1997) doi: 10.1103/PhysRevLett.78.417
    [35] W. M. Jin, Classical and Quantum Gravity, 15(10): 3163 (1998) doi: 10.1088/0264-9381/15/10/018
    [36] C. Doran and A. Lasenby, Phys. Rev, D, 66: 024006 (2002) doi: 10.1103/PhysRevD.66.024006
    [37] E. Jung, S. H. Kim, and D. K. Park, JHEP, 09: 005 (2004)
    [38] H. Cho and Y. Lin, Classical and Quantum Gravity, 22(5): 775 (2005) doi: 10.1088/0264-9381/22/5/001
    [39] A. B. Gaina and G. A. Chizhov, Moscow Univ. Phys. Bull., 38N2: 1-7 (1983)
    [40] A. B. Gaina, Moscow VINITI, No. 1970-80 Dep. , 20 pp. (1980)
    [41] M. Rogatko and A. Szyplowska, Phys. Rev. D, 79: 104005 (2009) doi: 10.1103/PhysRevD.79.104005
    [42] H. Liao, J-H. Chen, P. Liao, and Y-J. Wang, Commun. Theor. Phys., 62: 227-234 (2014)
    [43] A. Ghosh and P. Mitra, Phys. Rev. D, 50: 7389-7393 (1994)
    [44] I.I. Cotăescu, Mod. Phys. Lett. A, 13: 2923 (1998) doi: 10.1142/S0217732398003107
    [45] V. M. Villalba, Mod. Phys. Lett. A, 8: 2351-2364 (1993)
    [46] I. I. Cotăescu, J. Phys. A Math. Gen., 33: 1977 (2000)
    [47] B. Thaller. The Dirac Equation. (Springer Verlag, Berlin Heidelberg, 1992)
    [48] V. B. Berestetski, E. M. Lifshitz, and L. P. Pitaevski, Quantum Electrodynamics, (Oxford: Pergamon Press, 1982)
    [49] I. I. Cotăescu, Phys. Rev. D, 60: 124006 (1999) doi: 10.1103/PhysRevD.60.124006
    [50] I. I. Cotăescu, Phys. Rev. D, 65: 084008 (2002) doi: 10.1103/PhysRevD.65.084008
    [51] I. I. Cotaescu and C. Crucean, Mod. Phys. Lett. A, 22: 3707-3720 (2008)
    [52] I. I. Cotaescu, R. Racoceanu, and C. Crucean, Mod. Phys. Lett. A, 21: 1313-1318 (2006) doi: 10.1142/S0217732306020044
    [53] I.I. Cotăescu, Mod. Phys. Lett. A, 22: 2493 (2007)
    [54] C.A. Sporea, Mod. Phys. Lett. A, 30: 1550145 (2015)
    [55] C. Crucean, Mod. Phys. Lett. A, 22: 2573-2585 (2007) doi: 10.1142/S0217732307025480
    [56] C. Crucean, Mod. Phys. Lett. A, 25: 1679-1687 (2010) doi: 10.1142/S0217732310033414
    [57] C. Crucean, R. Racoceanu, and A. Pop, Phys. Lett. B, 665: 409-411 (2008) doi: 10.1016/j.physletb.2008.06.034
    [58] I. I. Cotaescu and C. Crucean, Int. J. Mod. Phys. A, 23: 1351-1359 (2008) doi: 10.1142/S0217751X08039487
    [59] M. Rose. Relativistic Electron Theory (New York: John Wiley & Sons, 1961)
    [60] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, NIST Handbook of Mathematical Functions (Cambridge: Cambridge University Press, 2010)
    [61] D. R. Yennie, D. G. Ravenhall, and R. N. Wilson, Phys. Rev., 95: 500 (1954) doi: 10.1103/PhysRev.95.500
    [62] K. W. Ford and J. A. Wheeler, Ann. Phys. (N.Y.), 7: 259 (1959)
    [63] K. W. Ford and J. A. Wheeler, Ann. Phys. (N.Y.), 7: 287 (1959)
    [64] P. Anninos, C. DeWitt-Morette, R. A. Matzner, P. Yioutas, and T. R. Zhang, Phys. Rev. D, 46: 4477 (1992) doi: 10.1103/PhysRevD.46.4477
    [65] R. A. Matzner, C. DeWitte-Morette, B. Nelson, and T.-R. Zhang, Phys. Rev. D, 31: 1869 (1985)
  • 加载中

Cited by

1. Sporea, C.A.. Scattering of quantum fields by a MOG black hole[J]. Modern Physics Letters A, 2020, 35(15): 2050121. doi: 10.1142/S0217732320501217
2. Sporea, C.A.. Scattering of spin half quantum fields by a Kiselev black hole[J]. AIP Conference Proceedings, 2020. doi: 10.1063/5.0001038
3. Sporea, C.A.. Quasibound states of the Dirac field in Schwarzschild and Reissner-Nordström black hole backgrounds[J]. Modern Physics Letters A, 2019, 34(39): 1950323. doi: 10.1142/S0217732319503231
4. Rodrigues, M.E., Silva, M.V.D.S. Regular multihorizon black holes in f (G) gravity with nonlinear electrodynamics[J]. Physical Review D, 2019, 99(12): 124010. doi: 10.1103/PhysRevD.99.124010

Figures(7)

Get Citation
Ciprian A. Sporea. Fermion scattering by a class of Bardeen black holes[J]. Chinese Physics C. doi: 10.1088/1674-1137/43/3/035104
Ciprian A. Sporea. Fermion scattering by a class of Bardeen black holes[J]. Chinese Physics C.  doi: 10.1088/1674-1137/43/3/035104 shu
Milestone
Received: 2018-10-03
Revised: 2018-12-29
Article Metric

Article Views(3138)
PDF Downloads(29)
Cited by(4)
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:

Fermion scattering by a class of Bardeen black holes

    Corresponding author: Ciprian A. Sporea, ciprian.sporea@e-uvt.ro
  • West University of Timişoara, V. Pârvan Ave. 4, RO-300223 Timişoara, Romania

Abstract: In this study, the scattering of fermions by a class of Bardeen black holes is investigated. After obtaining the scattering modes by solving the Dirac equation in this geometry, we use the partial wave method to derive an analytical expression for the phase shifts that enter into the definitions of partial amplitudes that define the scattering cross sections and induced polarization. It is shown that, similar to Schwarzschild and Reissner-Nordström black holes, the phenomena of glory and spiral scattering are present.

    HTML

    1.   Introduction
    • The existence of black holes is theoretically well motivated in general relativity (GR) and other modified theories of gravity. Thus far, black holes have not been directly detected or observed. However, plenty of indirect evidence indicate that such objects exist in nature. As shown by Penrose and Hawking [1], the occurrence of singularities is inevitable in GR. This raises many issues such as the black hole information paradox [2-4] or the more recently discovered paradox of black hole's "firewalls" [5-9]. These paradoxes emerge because of the incompatibility between the quantum theory and GR. It is widely believed that in the quantum theory of gravity, the singularities contained in black holes will be removed. However, even in the early stages of the investigations on singularities [10, 11] in GR, several black hole models, which could avoid the occurrence of a singularity, have been proposed. These black holes are called "regular," because they are singularity-free.

      The first proposal of a regular black hole solution was made by Bardeen in Ref. [12], and since then many other models of spherically symmetric regular black holes have been presented in literature [13-21]. In Ref. [22] the authors showed that the Bardeen black hole model can be physically interpreted as the gravitational field produced by a nonlinear magnetic monopole. Later, this interpretation was extended to also include nonlinear electric charges so that regular black holes models can have a nonlinear electromagnetic field as a source. More recently, in Ref. [23] the authors extended the Bardeen solution to an entire class of Bardeen-like black holes that can or cannot be regular.

      In this study, the scattering of fermions (spin 1/2) by a Bardeen regular black hole and by a Bardeen-class of black holes (as constructed in [23]) is investigated. We use the partial wave method to obtain analytical expressions for the phase shifts that enter into the definition of partial amplitudes, defining the scattering cross sections and the induced polarization. To the best of our knowledge, this is the first study that reports analytical phase shifts for 1/2-spin wave scattering by regular black holes. The absorption of fermions by Bardeen black holes was numerically investigated in [24], while in Refs. [25-28], the case of massless scalar scattering by regular black holes was also treated numerically. Studies dedicated to fermion scattering by other types of spherically symmetric black holes can be found as examples in Refs. [29-43].

      Both Bardeen regular black holes and the Bardeen-class type of black holes possess nonlinear magnetic (monopole) charges [22, 23]. This implies the existence of an electromagnetic potential of the form A=Qmcosθdϕ , with Qm the total magnetic charge. In this work, we neglect the interaction between this potential and the charge of the fermion and only focus on the scattering due to the "pure" gravitational interaction between the fermion and the black hole. However, even in this approximation, the black hole magnetic monopole charge will still influence the scattering patterns through the presence of Qm into the metric function and into the resulting scattering modes of the radial Dirac equation.

      The remainder of the paper is structured as follows. In Section 2, the Bardeen-class of black holes is presented very briefly. Section 3 begins with a very short review of the Dirac equation in spherically symmetric black hole geometries, and then it discusses the search for scattering modes in the Bardeen spacetime. Finally, the main result of the paper is presented, namely the form of the analytical phase shifts resultant from applying the partial wave method on the scattering modes derived earlier. Section 4 presents a graphical analysis of the induced polarization and the scattering cross sections in which the presence of a backward "glory" and "spiral scattering" (orbiting) oscillations are shown to be present. The main conclusions and some final remarks are presented in Section 5.

    2.   A class of Bardeen-like black holes
    • In Ref. [23], a class of spherically symmetric and asymptotically flat Bardenn-like black holes depending on two-parameters was constructed with the following line element

      ds2=h(r)dt2dr2h(r)r2(dθ2+sin2θdϕ2),h(r)=12Msr2Mrγ1(r2+Q2)γ.

      (1)

      These black holes have a singularity if the parameter Ms , hereinafter, referred to as the "Schwarzschild mass," has a nonzero value. Otherwise, if Ms=0 , regular black holes are obtained and the Bardeen black hole corresponds to the particular choice Ms=0 and γ=3 . The term M can be interpreted as the mass of the nonlinear magnetic monopole. Moreover, the sum

      MADM=Ms+M

      (2)

      constitutes the ADM mass of the black hole obtained from the asymptotic form of the metric function h(r). The parameter Q in Eq. (1) is related to the magnetic monopole charge Qm ; Q=2Q2m/M .

      As shown in Ref. [23], the Lagrangian density, for which Eq. (1) is a solution of the coupled Einstein-Maxwell filed equations, is given by

      L=4γα(αFμνFμν)5/4(1+αFμνFμν)1+γ/2,

      (3)

      where α has the dimensions of length squared, and γ is a dimensionless constant. In the weak field limit, a vector field, which is slightly stronger when compared with a Maxwell field, is obtained.

    3.   Dirac fermions and scattering cross sections

      3.1.   Dirac equation. Preliminaries

    • The Dirac equation:

      iγaDaψmψ=0

      (4)

      The following explicit form can be obtained from the Dirac equation [44]

      (iγaeμaμm)ψ+i21gμ(geμa)γaψ14{γa,Sbc}ωcabψ=0,

      (5)

      where g=det(gμν) and the covariant derivative is defined by Da=a+i2Sbcωcab . Sab=i4[γa,γb] the generators of the SL(2,C) group, γa are the point-independent Dirac matrices obeying {γa,γb}=2ηab and ωcab is a spin-connection

      ωcab=eμaeνb(ˆecλΓλμνˆecν,μ)

      (6)

      Γλμν represents the usual Christoffel symbols. The tetrad fields ea(x) and ˆea(x) are point-dependent defining (nonholonomic) local frames and co-frames. The following relations hold:

      ea=eμaμ,ˆea=ˆeaμdxμea(x)eb(x)=ηab,ˆea(x)ˆeb(x)=ηabds2=ηabˆeaμdxμˆebνdxν=gμν(x)dxμdxν.

      (7)

      The so-called Cartesian gauge [44-46] for a spherically symmetric line element of the form of Eq. (1) is defined by the following tetrad fields:

      ˆe0=hdtˆe1=1hsinθcosϕdr+rcosθcosϕdθrsinθsinϕdϕˆe2=1hsinθsinϕdr+rcosθsinϕdθ+rsinθcosϕdϕˆe3=1hcosθdrrsinθdθ,

      (8)

      The Dirac equation (5) can be reduced to only a radial equation. The angular part of the Dirac equation is the same as in the Dirac theory from flat spacetime, and its solutions are the usual four-component angular spinors Φ±m,κ(θ,ϕ) [47, 48]. This is because, in the Cartesian gauge (8), the Dirac equation is manifestly covariant under rotations [46]. Using this gauge, it was possible to find complete analytical solutions to the Dirac equation on de Sitter/anti-de Sitter spacetime [49-52] and approximate analytical solutions in black hole geometries [31-33, 53, 54], which were later used to study different aspects of the scattering problem in those spacetimes [31-33, 55-58].

      The remaining unsolved radial part of the Dirac equation was obtained by assuming the following type of particle-like solutions with a given energy E

      ψ(x)=ψE,j,m,κ(t,r,θ,ϕ)=eiEtrh(r)1/4{f+E,κ(r)Φ+m,κ(θ,ϕ)+fE,κ(r)Φm,κ(θ,ϕ)}

      (9)

      where f±E,κ(r) are two unknown radial wave functions. As in Refs. [31, 53] the radial Dirac equation can be expressed in a matrix form

      (mh(r)h(r)ddr+κrh(r)h(r)ddr+κrh(r)mh(r))×(f+E,κ(r)fE,κ(r))=E(f+E,κ(r)fE,κ(r)).

      (10)
    • 3.2.   Scattering modes

    • Substituting the line element (1) in Eq. (10), we obtain a system of two differential equations that has no analytical solutions owing to the complex form of the line element. However, because we are only interested in finding the scattering modes, we can approximate Eq. (10) in the asymptotic region of the Bardeen black hole specified by the line element in Eq. (1) and find approximative analytical solutions. A partial wave method was used on the obtained solutions to compute the (elastic) scattering cross section and the induced polarization that results after the interaction of a fermion beam with the black hole.

      Let us start by introducing the new variable

      x=zr+1,z=r2+Q2,

      (11)

      where r+ is the radius of the black hole horizon. The function h(r), in terms of x, is

      h=12Msr+1(1+x2)2δ2Mr+11+x2[1δ(1+x2)2]γ12

      (12)

      where the notation δ=(Q/r+)2 was introduced. The radial Dirac equation (10) is a system of two differential equations for the radial wave functions f±(x) that, after multiplying each equation by x/[r+(1+x2)] and expressing all the terms as a function of x given by Eq. (11), is equivalent to Eq. (13)

      [11+x21δ(1+x2)2{...}12ddx+κ(1+x2)2δ×{...}12]f+E,κ(x)[μ{...}12+ϵ(x+1x)]fE,κ(x)=0[11+x21δ(1+x2)2{...}12ddxκ(1+x2)2δ×{...}12]fE,κ(x)[μ{...}12ϵ(x+1x)]f+E,κ(x)=0where{...}=(1+x2)2x22Msr+(1+x2)2x2(1+x2)2δ2Mr+1+x2x2[1δ(1+x2)2]γ12,μ=r+m,ϵ=r+E.

      (13)

      The system of differential equations obtained in Eq. (13) cannot directly be solved analytically. However, an analytical solution can be obtained if a domain is restricted far away from the black hole event horizon. A simpler system of differential equations that has analytical solutions can be obtained. By Taylor expansion with respect to 1/x and discarding the terms of the order O(1/x2) and higher, the system of equations, which is valid in the asymptotic region of the black hole for the two radial wave function, reduces to

      (12ddx+κx)f+E,κ(x)x(ϵ+μ)fE,κ(x)1x[ϵ+μ(1MADMr+)]fE,κ(x)=0(12ddxκx)fE,κ(x)+x(ϵμ)f+E,κ(x)+1x[ϵμ(1MADMr+)]f+E,κ(x)=0.

      (14)

      We observe that if M0 then the ratio MADMr+Msrs=12 and μ(1MADMr+)12μ , recovering the Schwarzschild case discussed in Ref. [53].

      To find the scattering modes for the case ϵ>μ , it is useful to introduce new radial wave functions ˆf± , present in the following combination of the old wave functions f±

      ˆf+E,κ=i2f+E,κϵ+μ+12fE,κϵμˆfE,κ=i2f+E,κϵ+μ+12fE,κϵμ.

      (15)

      Using the above relation we can obtain the following equations, which are satisfied by the functions ˆf+E,κ and ˆfE,κ :

      [νxddx+2i(μ2(1MADMr+)ϵ2ν2x2)]ˆf+E,κ(2κνiϵμMADMr+)ˆfE,κ=0[νxddx2i(μ2(1MADMr+)ϵ2ν2x2)]ˆfE,κ(2κν+iϵμMADMr+)ˆf+E,κ=0,

      (16)

      where ν=ϵ2μ2 . Eq. (16) can now be solved using Maple or Mathematica and the analytical solutions can be written as a combination of Whittaker functions

      ˆf+E,κ=C11xMρ+,s(z)+C21xWρ+,s(z)ˆfE,κ=1κ2+λ2[(siα)(κ+iλ)C11xMρ,s(z)(κ+iλ)C21xWρ,s(z)],

      (17)

      where z=2iνx2 and the following parameters were introduced

      s=[κ2+μ2(1MADMr+)2ϵ2]12,λ=ϵμνMADMr+ρ±=12iα,α=1ν[ϵ2μ2(1MADMr+)].

      (18)

      3.2.1 Pure Bardeen black hole case

      As already mentioned, by choosing Ms=0 and γ=3 in Eq. (1) the original Bardeen black hole solution is recovered. This solution is a black hole with two distinct horizons only if |Q|<4/27 , it has degenerate horisons if Q=4/27 and for Q>4/27 there are no horizons present. Because the function h(r) has now a more simpler form,

      h(r)=12Mr2(r2+Q2)3,

      (19)

      one can find and write an analytical expression for the location of the black hole outer horizon, denoted from now on by r+ . Moreover, one can easily show that r+=Mf(Q/M) , where the function f depends only on the ratio Q/M. The existence of horizons implies the constraint |Q|/M4/27 [14].

      Following the same steps as in section 3.2, one finds the same scattering modes as given by Eq. (17), but with the new parameters

      s=[κ2+μ2(1Mr+)2ϵ2]12,λ=ϵμνMr+ρ±=12iα,α=1ν[ϵ2μ2(1Mr+)].

      (20)

      We will see in section 4 that this is enough to produce a noticeable difference in the scattering patterns.

      By choosing Ms=0 and γ=3 in Eq. (1), the original Bardeen black hole solution can be recovered. This solution is a black hole with two distinct horizons only if |Q|<4/27 ; it has degenerate horizons if Q=4/27 , and for Q>4/27 there are no horizons present. Because the function h(r) can now be expressed in a simpler form:

    • 3.3.   Analytical phase shifts and scattering cross sections

    • In spinor wave scattering theory [48, 59] the differential scattering cross section for an unpolarized incident beam is the sum of the squares of two scalar functions

      dσdΩ=|f(θ)|2+|g(θ)|2,

      (21)

      which only depend on the scattering angle θ :

      f(θ)=l=012ip[(l+1)(e2iδl11)+l(e2iδl1)]P0l(cosθ),g(θ)=l=112ip[e2iδl1e2iδl]P1l(cosθ),

      (22)

      where p is the incident momentum, and P0l(cosθ),P1l(cosθ) are the Legendre and associated Legendre polynomials, which are special cases of Legendre functions, respectively [60].

      The phase shifts δl can be computed as in Refs. [31, 32] by applying the partial wave method on the scattering modes (17). The asymptotic form of the radial wave functions f± can be written as [31]

      (f+E,κfE,κ)E+msinEmcos(prπl2+δκ+ϑ(r))

      (23)

      where ϑ(r)=pr++αln[2p(rr+)] represents a radially dependent phase that is independent of any angular quantum numbers; therefore, it does not contribute to the scattering cross sections and may be neglected as in the Dirac-Coulomb case [30, 48].

      The resultant final form for the point-independent phase shifts δκ is given by the following expression (see also Appendix A)

      e2iδκ=κiλsiαΓ(1+siα)Γ(1+s+iα)eiπ(ls),

      (24)

      where the sign convention for κ is the same as in [48], such that κ=±(j+1/2) and l=|κ|(1signκ)/2 .

      Series (22) are poorly convergent as a direct consequence of the singularity present at θ=0 ; an infinite number of Legendre polynomials are necessary to describe it. To make the series more convergent, the mth reduced series can be defined

      (1cosθ)m1f(θ)=l0a(m1)lPl(cosθ),(1cosθ)m2g(θ)=l1b(m2)lP1l(cosθ),

      (25)

      as first proposed in [61] and more recently used in [30-32]. The coefficients a(i)l and b(i)l are computed using the recurrence relations

      a(i+1)l=a(i)ll+12l+3a(i)l+1l2l1a(i)l1,b(i+1)l=b(i)ll+22l+3b(i)l+1l12l1b(i)l1,

      (26)

      where a(0)l=[(l+1)(e2iδl11)+l(e2iδl1)]/2ip and b(0)l=[e2iδl1e2iδl]/2ip from Eq. (22). We found that using only two iterations m1=2 for the function f(θ) and one iteration m2=1 for g(θ) , it is sufficient to make the series converge without distorting the analytical results.

    4.   Results and discussion
    • In this section, we present and discuss the main features of fermion scattering by Bardeen regular black holes and also by a Bardeen-class of black holes. The analysis will focus on scattering by small or micro black holes (with MADM10151022 Kg) because in this case the glory and orbiting scattering phenomena are shown to be significant.

      The following parameters are used to label the figures: v=p/E is the speed of the incident fermions; EM can be considered a dimensionless measure of the gravitational coupling because (restoring the units) it forms the dimensionless quantity: ε=GEMc3=πrsλv , with rs the Schwazschild radius and λ=h/p the associated quantum particle wavelength; and the ratios q=Q/Ms and g=M/Ms that appear when writing the black hole horizon radius as r+=Msf(MMs,QMs) . Moreover, because in the asymptotic zone E=m2+p2 , one can easily show that mM=EM1v2 such that the condition EMmM is always satisfied. In the following analysis we take γ=3 in all the plots (except those in Fig. 7) because the same conclusions are also obtained for the cases with γ3 .

      Figure 7.  (color online) Comparison of the scattering cross section for a Bardeen-class black hole with different values of the parameter γ . Note that the scattering profile keeps the same shape for different values of γ .

      In Fig. 1, the differential scattering cross section is presented as a function of the scattering angle for fixed values of the ratios: Q/M (left panels corresponding to pure Bardeen case); q=Q/Ms and g=M/Ms ; a fixed value of the speed (v) of incoming fermions, while the parameter EM (or EMs for Bardenn-class) takes different values. From Fig. 1, it can be seen that the scattering pattern takes a simple form for small values of EM. However, as the value of EM increases, more complex scattering patterns occur. This includes the presence of a maximum in the backward direction ( θ=π ), which is also known as "glory" scattering [62, 63] and the presence of oscillations in the scattering intensity that give rise to orbiting or "spiral scattering" [64, 65] (that may occur when the particle's "classical" orbit passes the scattering center multiple times). As the speed of the incoming fermion increases, the peak in the π -direction moves toward the left, and the maximum of the scattering intensity that occurs at θ=π for nonrelativistic fermions transforms into a minimum if the fermions are massless (v= 1). From the bottom panels of Fig. 1, it can be seen that the magnitude of the spiral scattering oscillations and their angular frequency increase with the black hole mass.

      Figure 1.  (color online) Plot of the differential scattering cross section as a function of the scattering angle for a regular Bardeen black hole, Eq. (19) (left panels) and for a Bardeen-class black bole, Eq. (1) (right panels). Phenomena of glory (scattering in the backward direction) and spiral scattering (oscillations in the scattering intensity) are present for both types of black holes. The value q/Mem=4/27 corresponds to a degenerate regular Bardeen black hole (it has only one horizon). The parameter q gives the ratio q=Q/Ms and g=M/Ms is the ratio between the magnetic monopole mass and the "Schwarzschild mass".

      The scattering cross section dσ/dΩ for the Bardeen-class black holes (1) is plotted in Fig. 2. The left panels show the scattering patterns for Bardeen-class black holes that have the "Schwarzschild mass" greater than the mass of the magnetic monopoles, while the right side panels show those with a lower mass. By comparing our analysis results with the Schwarzschild black hole, the scattering intensity in the backward direction ( θ=π ) is higher for the scattering by a Bardeen-class black hole. From Fig. 2, one can also observe that in all the plots, the value of g increases as the number of oscillations of the scattering intensity increases. Furthermore, if we assume Ms is fixed, then the ADM-mass of the black hole increases (which is equivalent to increasing the mass M, Eq. (2), if Ms remains constant), and as a consequence the angular frequency of the oscillations present in the scattering intensity also increase. This feature was also true for fermion scattering by Schwarzschild and Reissner-Nordström black holes [30-33]. Thus, the increase in the frequency of the oscillations of the spiral scattering with the black hole mass for any type of (spherically symmetric) black holes is of interest.

      Figure 2.  (color online) Variation of the scattering cross section with the parameter g=M/Ms for given values of EMs and v. In the top panels, the Schwarzschild scattering (for which g= 0) is compared. The value of g increases as the spiral scattering becomes more pronounced.

      Comparing the scattering intensity, in Fig. 3, of a Bardeen regular black hole (blue and red curves) with that of a Schwarzschild black hole (the black dotted curve), we observe that the glory peak is higher for the scattering by a Schwarzschild black hole. As the ratio Q/M increases, the maximum in the backward direction becomes lower, and the frequency of the oscillations in the spiral scattering slightly decrease simultaneously. The curve with Q/M=4/27 corresponds to the Bardeen degenerate case when only one horizon is present. The type of black hole can be determined by analyzing the patterns in the scattering intensity when the scattering of a beam of fermions by a black hole can be observed and measured. If it is a Bardeen black hole, with a known mass, then the value of Q can also be found.

      Figure 3.  (color online) Scattering cross section for a regular Bardeen black hole for different values of Q/M, while the parameters EM and v are kept fixed.

      In Fig. 4, we plotted the differential scattering cross section dσdΩ as a function of the scattering angle θ for a Schwarzschild black hole (blue solid lines), regular Bardeen black hole (black dotted lines) and Bardeen-class black hole (red dash-dotted lines).

      Figure 4.  (color online) Comparison between the fermion scattering by a Schwarzschild black hole, regular Bardeen black hole, and Bardeen-class black hole. For the left panel Q/M= 0.2, q= 0.2, and g= 0.5 were used, while for the right panel Q/M= 0.4, q= 0.1, and g= 0.2.

      An incident unpolarized beam of massive fermions becomes partially polarized after it gets scattered by the black hole. The induced polarization degree can be computed using the following formula [30, 59]

      P=ifgfg|f|2+|g|2n

      (27)

      where n is a unit vector in a direction orthogonal to the plane of scattering.

      In Fig. 5, the dependence of the polarization on the scattering angle is plotted for given values of the parameters EM, v, g, and Q/M. It has a pronounced oscillatory behavior. From the top panels in Fig. 5, one can observe that if the parameter EM increases, then the frequency of the oscillations present in the polarization increases as well. Now, if we assume that the energy of the fermion is fixed, then, the oscillations present in the polarization are more pronounced for black holes with higher masses. Compared with the Schwarzschild polarization (see bottom panels in Fig. 5) the Bardeen and Bardeen-class black hole polarizations exhibit a slightly less oscillatory behavior. The oscillations present in the polarization can be observed as a consequence of the oscillations present in the glory and spiral scattering.

      Figure 5.  (color online) Partial polarization P as a function of the scattering angle θ . Top left panel: typical Bardeen black hole (19) with Q/M= 0.4 for incoming fermions with speed v= 0.8 in units of c; Top right panel: Bardeen-class black hole (1) with g=M/Ms=0.2 for incoming fermions with speed v=0.4 in units of c; Bottom left panel: comparing Schwarzschild and pure Bardeen black hole polarizations for fixed EM= 1 and v= 0.3; Bottom right panel: comparing Schwarzschild and Bardeen-class black hole polarizations for fixed EMs=1.5 , q=Q/Ms=0.2 and v= 0.5.

      The alignment of the scattered fermions with the forward on-axis direction can be visualized using the polar plots representations of the polarization degree as shown in Fig. 6. One can observe the Mott polarization in the direction orthogonal to the scattering plane, a phenomena that has been reported before in literature for Schwarzschild [30, 31] and Reissner-Nordström [32] black holes.

      Figure 6.  (color online) Polar plots of pure Bardeen (left panel) and Bardeen-class (right panel) polarization P(θ) for 0θ2π showing the alignment of the scattered fermion's spin with a given direction. The Mott polarization (polarization in the direction orthogonal to the scattering plane) can also be observed.

      In Fig. 7, the differential scattering cross section is plotted for a Bardeen-class black hole with different values of the parameter γ . Note that the scattering pattern maintains the same profile shape, and the difference between the values of dσ/dΩ (computed using different values of γ ) is very small. Moreover, from our analysis results, it can be seen that, for most cases, this difference is even smaller than for the examples presented in Fig. 7. This can be attributed to the difference between the black hole horizon radii r+ obtained for different values of γ (while all the other parameters are kept constant) is very small.

    5.   Conclusions
    • We studied the scattering of fermions by a class of Bardeen black holes that also include the original Bardeen regular black hole solution. A partial wave method was used on a set of scattering modes that were obtained by solving the Dirac equation in the asymptotic region of these black hole geometries. TO the best of our knowledge, for the first time, time analytical phase shifts for fermion scattering by Bardeen regular black holes were obtained. The phenomena of glory (scattering in the backward direction) and spiral scattering (oscillations in the scattering intensity) were also discussed. We also saw that an incident unpolarized beam could become partially polarized after interaction with a black hole.

      In Figs. 1-6, other than the parameters EM (that can be associated with a measure of the gravitational coupling) and v (speed of the fermion) we also used the ratios Q/M, q=Q/Ms , and g=M/Ms to label the figures. The departure of the scattering pattern from the Schwarzschild case was significant as Q/M and g increased. For the original Bardeen regular black hole, Q/M=4/27 was the maximum allowed value, and it corresponds to the degenerate case when the two horizons coincide. If the magnetic charge Q0 , then M0 and g0 , and the scattering by a Schwarzschild black hole can be recovered [31].

      The glory and spiral scattering become significant for values of the parameter EM~1 or greater (in geometrical units with G==c=1 ) because the associated wavelength for these values of the incident fermions is of the same order of magnitude as the black hole horizon radius; thus, diffraction patterns occur. Another feature, which was also shown to be present for fermion scattering by Schwarzschild [30, 31] and Reissner-Nordström black holes [32], is that as the total mass of the black hole increases, the oscillations present in the scattering intensity become more frequent, that is, spiral scattering is significantly enhanced with an increase in black hole mass.

    Appendix A
    • The aim of this Appendix is to briefly present how the phase shifts (24) are obtained.

      The asymptotic representation of the Whittaker function M for large values of |z| is given by the following equation [60]

      Mκ,μ(z)Γ(1+2μ)Γ(12+μκ)e12zzκ(1+O(z1))+ei(12+μκ)πΓ(1+2μ)Γ(12+μ+κ)e12zzκ(1+O(z1))

      valid for 12π<phz<32π .

      To obtain the phase shifts (24), the condition C2=0 must be imposed on Eq. (17) in order to obtain the correct Newtonian phase shifts for large values of angular momenta. As an observation, this condition also selects the spinors that are regular in x= 0 owing to the regularity of the function Mρ±,s(2iνx2)=(2iνx2)s+12[1+O(x2)] . More arguments for letting C2=0 can be found in Appendix C of our previous paper [31].

      Using Eq. (28) on the ˆf±(x) functions and by observing that for ˆf+ the dominant term is the first one in (28), while for ˆf the second term is dominant, then one obtains the following asymptotic expressions

      ˆf+(x)C1e12παΓ(2s+1)Γ(1+s+iα)ei[νx2+αln(2νx2)]ˆf(x)(siα)(κ+iλ)κ2+λ2C1e12παΓ(2s+1)Γ(1+siα)×ei[νx2+πs+αln(2νx2)].

      The last step is to compute the argument of the trigonometric functions (23) as 12arg(ˆf+ˆf) from which the expression for e2iδκ will result.

Reference (65)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return