By employing a pseudoorthonormal coordinate-free approach, the Dirac equation for particles in the Kerr–Newman spacetime is separated into its radial and angular parts. In the massless case to which a special attention is given, the general Heun-type equations turn into their confluent form. We show how one recovers some results previously obtained in literature, by other means.
1. Introduction
After Carter found that the scalar wave function is separable in the Kerr–Newman–de Sitter geometries [1], the solutions to the Teukolsky equations [2], for massless fields in the Kerr metrics, have been analytically expressed in the form of a series of various functions [3, 4].
Starting with the work of Chandrasekhar [5], general properties of a massive Dirac field equation in the Kerr background have been extensively studied.
The recent interest in the so-called quasinormal modes of a Dirac field in the Kerr background is motivated by the detection of gravitational waves [6–8], whose phase can be described in terms of the proper oscillation frequencies of the black hole.
In terms of techniques, after the Dirac equation in the Kerr–Newman background was separated [9, 10], using the Kinnersley tetrad [11], the Newman-Penrose formalism [12] has been considered a valuable tool for dealing with this subject [13]. This formalism as well as the Geroch–Held–Penrose variant has been used for the Teukolsky Master Equation describing any massless field of different spins, in the Kerr black hole and for an arbitrary vacuum spacetime [14, 15].
In [16, 17], it was shown that, for Kerr-de Sitter and Kerr–Newman–de Sitter geometries, both angular and radial equations for the Teukolsky equation, for massless fields, are transformed into Heun’s equation [18, 19] and analytic solutions can be derived in the form of a series of hypergeometric functions. More recently, solutions of the Dirac equation in the near horizon geometry of an extreme Kerr black hole have been found analytically in terms of the confluent Heun functions in [20].
The massive case was tackled within the WKB approach [21] or numerically, using the convergent Frobenius method [22]. Very recently, in [22], after a tedious calculation, using a generalised Kinnersley null tetrad in the Newman-Penrose formalism, the Dirac equation for a massive fermion has been separated in its radial and angular parts, the solutions being expressed in terms of generalised Heun functions.
Our work is proposing an alternative, free of coordinates, method based on Cartan’s formalism. Thus, we are computing all the geometrical essentials for dealing with the Dirac equation in its SO3,1×U1-gauge covariant formulation. Our approach is generalizing the theory developed in [23], where for massless fermions on the Kerr spacetime, the authors are switching between canonical and pseudoorthonormal bases and the solutions are derived using numerical techniques.
By imposing the necessary condition for a polynomial form of the Heun confluent functions [18, 19], we obtain the resonant frequencies, which are of crucial importance for getting information on the black holes interacting with different quantum fields [24]. In the last years, the Heun functions in either their general or confluent forms have been obtained by many authors, for example, [25–40] and the references therein.
The structure of this paper is as follows: in Section 2, we present all the necessary ingredients for writing down the massive Dirac equation in the Kerr–Newman background. We show that, by using an orthonormal tetrad adapted to the Kerr–Newman metric, one can separate the massive Dirac equation. For slowly rotating objects, the solutions of the radial equations can be expressed in terms of the confluent Heun functions. As an application, we compute the modal radial current vector. In Section 3, we turn our attention to the massless Dirac fermions and show that the Dirac equations can be solved exactly in the Kerr case and also in the extremal case, a result previously known in literature, obtained by other means. The final section is dedicated to conclusions.
2. The SO3,1×U1-Gauge Covariant Dirac Equation
Let us start with the four-dimensional Kerr–Newman metric in the usual Boyer–Lindquist coordinates:
(1)ds2=ρ2Δdr2+ρ2dθ2+sin2θρ2adt−r2+a2dφ2−Δρ2dt−asin2θdφ2,where Δ=r2−2Mr+a2+Q2, ρ2=r2+a2cos2θ, and M, Q, and a are the black hole’s mass, charge, and angular momentum per unit mass. The electromagnetic background of the black hole is given by the four-vector potential, in coordinate basis:
(2)Aidxi=Qrρ2dt−asin2θdφ.
Within a SO3,1-gauge covariant formulation, we introduce the pseudoorthonormal frame Eaa=1,4¯, whose corresponding dual base is
(3)Ω1=ρdθ,Ω2=sinθρr2+a2dφ−asinθρdt,Ω3=ρΔdr,Ω4=−aΔρsin2θdφ+Δρdt,leading to the expressions
(4)dθ=1ρΩ1,dφ=1ρsinθΩ2+aρΔΩ4,dr=ΔρΩ3,dt=aρsinθΩ2+r2+a2ρΔΩ4.
Thus, using the relations gikdxidxk=gikEaiEbkΩaΩb=δabΩaΩb, i.e., dxi=EaiΩa, one may write down the pseudoorthonormal frame:
(5)E1=1ρ∂θ,E2=1ρsinθ∂φ+aρsinθ∂t,E3=Δρ∂r,E4=aρΔ∂φ+r2+a2ρΔ∂t.
Using (3), the first Cartan’s equation
(6)dΩa=Γ.bcaΩb∧Ωc,with 1≤b<c≤4 and Γ.bca=Γ.bca−Γ.cba, can be explicitely worked out as
(7)dΩ1=−Δρ2ρ,3Ω1∧Ω3,dΩ2=1sinθsinθρ,1Ω1∧Ω2−rΔρ3Ω2∧Ω3+2arsinθρ3Ω3∧Ω4,dΩ3=ρ,1ρ2Ω1∧Ω3,dΩ4=−2aΔρ3cosθΩ1∧Ω2−a2ρ3sinθcosθΩ1∧Ω4+Δρ,3Ω3∧Ω4,where ·,1 and ·,3 are the derivatives with respect to θ and r, leading to the following complete list of nonzero connection coefficients in the Cartan frames Ωa,Eaa=1,4¯:
(8)Γ122=−Γ212=−1sinθsinθρ,1,Γ124=−Γ214=−Γ412=−aΔρ3cosθ,Γ131=−Γ311=Δρ2ρ,3,Γ133=−Γ313=−ρ,1ρ2,Γ232=−Γ322=rΔρ3,Γ234=−Γ324=arρ3sinθ,Γ142=−Γ412=−aΔρ3cosθ,Γ144=−Γ414=−a2ρ3sinθcosθ,Γ241=−Γ421=Γ412=aΔρ3cosθ,Γ243=−Γ423=−Γ234=−arρ3sinθ,Γ342=−Γ432=−Γ234=−arρ3sinθ,Γ344=−Γ434=Δρ,3.
Now, one has all the essentials to write down the SO3,1×U1-gauge covariant Dirac equation for the fermion of mass μ:
(9)γaΨ;a+μΨ=0,where “;” stands for the covariant derivative:
(10)Ψ;a=Ψ∣a+14ΓbcaγbγcΨ−iqAaΨ,with Ψ∣a=EaΨ.
In view of the relations (8), the term expressing the Ricci spin connection
(11)14Γbcaγaγbγc=12Γ212+Γ313−Γ414γ1+12Γ131+Γ232−Γ434γ3+i2Γ241γ3γ5+i2Γ234γ1γ5has the concrete expression
(12)14Γbcaγaγbγc=12cotθρ−a2ρ3sinθcosθγ1+12Δ,3ρ+rΔρ3γ3+iar2ρ3sinθγ1γ5+iaΔ2ρ3cosθγ3γ5,where γ5=−iγ1γ2γ3γ4, while the kinetic term reads
(13)γaΨ∣a=1ργ1Ψ,1+γ21ρsinθΨ,2+aρsinθΨ,4+Δργ3Ψ,3+γ4aρΔΨ,2+r2+a2ρΔΨ,4.
Putting everything together, the Dirac equation (9) has the explicit form:
(14)γ11ρΨ,1+cotθ2ρ−a22ρ3sinθcosθ+iar2ρ3sinθγ5Ψ+γ21ρsinθΨ,2+aρsinθΨ,4+γ3ΔρΨ,3+Δ,32ρ+rΔ2ρ3+iaΔ2ρ3cosθγ5Ψ+γ4aρΔΨ,2+r2+a2ρΔΨ,4−iqγ4A4Ψ+μΨ=0,where the proper four-potential component, coming from Aicdxi=A4Ω4, with Aic given in (2), reads
(15)A4=QrρΔ.
For ease of calculations, the choice for γa matrices is important and we are going to employ Weyl’s representation:
(16)γμ=−iβαμ,γ4=−iβ,with
(17)αμ=σμ00−σμ,β=0−I−I0,so that
(18)γ5=−iγ1γ2γ3γ4=I00−I.
Thus, for the bispinor written in terms of two-component spinors as
(19)Ψ=ζη,the general equation (14) leads to the following system of coupled equations for the spinors ζ and η:
(20)σ11ρζ,1+cotθ2ρ+iasinθ2ρ3ρ+ζ+σ21ρsinθζ,2+aρsinθζ,4+σ3Δρζ,3+Δ,32ρ+Δ2ρ3ρ+ζ+aρΔζ,2+r2+a2ρΔζ,4−iqA4ζ−iμη=0(21)σ11ρη,1+cotθ2ρ−iasinθ2ρ3ρ−η+σ21ρsinθη,2+aρsinθη,4+σ3Δρη,3+Δ,32ρ+Δ2ρ3ρ−η−aρΔη,2−r2+a2ρΔη,4+iqA4η+iμζ=0,where ρ±=r±iacosθ and ρ2=ρ+ρ−.
Due to the time independence and symmetry of the spacetime, we can assume that the wave function can be written as
(22)ζ=Δ−1/4ρ−−1/2eimφ−ωtXρ,θ,η=Δ−1/4ρ+−1/2eimφ−ωtYρ,θ,where the factors Δ−1/4ρ±−1/2 have been introduced in order to pull some terms out of equations (20) and (21).
With the new functions Xρ,θ and Yρ,θ, equations (20) and (21) can be put into the transparent form:
(23)σ1DθX+iσ2HX+σ3ΔX,3+iKX−iμρ−Y=0,σ1DθY+iσ2HY+σ3ΔY,3−iKY+iμρ+X=0,where we have introduced the operators:
(24)Dθ=∂∂θ+cotθ2,H=msinθ−ωasinθ,K=1Δma−ωr2+a2−qQr.
Finally, by applying the separation ansatz
(25)X1=R1rT1θ,X2=R2rT2θ,Y1=R2rT1θ,Y2=R1rT2θ,one gets the system:
(26)R1Dθ−HT1−T2Δ∂r−iKR2−iμρ−R1T2=0,R2Dθ+HT2+T1Δ∂r+iKR1−iμρ−R2T1=0,R2Dθ−HT1−T2Δ∂r+iKR1+iμρ+R2T2=0,R1Dθ+HT2+T1Δ∂r−iKR2+iμρ+R1T1=0,which leads to the radial and angular equations:
(27)Δ∂r+iKR1=λ+iμrR2,Δ∂r−iKR2=λ−iμrR1,Dθ−HT1=λ+μacosθT2,Dθ+HT2=−λ+μacosθT1,where λ is a separation constant.
The first-order angular equations may be combined to obtain the so-called Chandrasekhar-Page angular equation and have been discussed in detail in [41].
From the radial equations in (27), one gets the following second-order differential equation for the R1 component:
(28)ΔR1′′+r−M−iμΔλ+iμrR1′+iΔK′+μKΔλ+iμr+K2−λ2−μ2r2R1=0,and i→−i, for R2.
Similar relations have been obtained in [22], by a different approach, namely, using the Newman-Penrose formalism. In the generalised Kinnersley frame, the null tetrad has been constructed directly from the tangent vectors of the principal null geodesics. Even though the radial and angular equations coming from (27) have been reduced to generalised Heun differential equations [18, 19], the solutions are not physically transparent since they look quite complicated and there are many open questions especially related to their normalization or to the behavior around the singular points.
However, for large values of the coordinate r, equation (28), with K given in (24), reads
(29)r21−2MrR1′′+M−a2+Q2rR1′+−2iωr−iqQ+2iΩrr−Mr−2M+Ω2r3r−2M−λ2−μ2r2R1=0,with the notation
(30)Ω=ω+qQr+a2r2ω−ma,where one may identify the fermion’s quanta energy, ω, the standard Coulomb energy, qQ/r, and the internal centrifugal energy with the quantum resonant correction, i.e., ω−m/a.
To the first order in a, meaning a slowly rotating object, for which
(31)Ω≈ω+qQr−mar2,Ω2≈ω+qQr2−2ωmar2,the solutions of (29) are given in terms of the Heun confluent functions [18, 19] as
(32)R1~eiprr−2M14+γ2r34C1rβ2HeunCα,β,γ,δ,η,r2M+C2r−β2HeunCα,−β,γ,δ,η,r2M,with the parameters written in the physical transparent form as
(33)α=4ipM,β=94−2imaM≈32−2ima3M,γ=4iMΩ∗+i4M2+364M21/2≈4iMΩ∗,δ=8M2ωω+qQ2M−μ22,η=58−λ2,where p2=ω2−μ2 and Ω∗ is the energy computed on the Schwarzschild horizon, i.e.,
(34)Ω∗=ω+qQ2M−ma4M2.
The second component, R2, is given by the complex conjugated expression of (32).
One may notice that, for Δ≈rr−2M and ρ+≈ρ−≈r, the first component in ζ defined in (22) reads
(35)ζ1=eαx/2eimφe−iωtx−1γ/2x±β/2HeunCα,±β,γ,δ,η,x,with x=r/2M.
Moreover, since R12=R22, the modal radial current (of quantum origin), computed as
(36)jr=iΨ¯γ3Ψ=Ψ†α3Ψ=ζ†σ3ζ−η†σ3η,vanishes. The only nonvanishing component is the azimuthal one, which is given by the expression
(37)jφ=iΨ¯γ2Ψ=Imeαxx−1γxβHeunCα,β,γ,δ,η,x2T1T2=x3/2Imexp4ipMx+4iMΩ∗lnx−1−2ima3MlnxHeunC2T1T2.
The current has the generic representation given in Figure 1, for x>1, i.e., r>2M. One may notice the oscillating behavior, with both positive and negative regions, vanishing at infinity. Also, there is a dominant positive maximum, just after the (Schwarzschild) horizon r=2M of the slowly rotating black hole, where the Heun functions have a regular singularity.
The radial part of the current (37), for x>1.
For the asymptotic behavior in the neighborhood of the singular point at infinity, where the two solutions of the confluent Heun equation exist, one may use the formula [24]
(38)HeunCα,β,γ,δ,η,x≈D1x−β+γ+22+δα+D2e−αxx−β+γ+22−δα=e−αx2x−β+γ+22D1eαx2x−δα+D2e−αx2xδα=De−αx2x−β+γ+22sin−iαx2+iδαlnx+σ,so that the two independent solutions in (32) are given by the simple expression
(39)R1=Dsinpr+2Mpωω+qQ2M−μ22logr2M+σ,where σω is the phase shift, D=const, and p=ω2−μ2.
Thus, the first component of Ψ defined in (19), (22), and (25) has the following (physical) behavior for large r values:
(40)ζ1≈R1reimφ−ωtT1θ,and similarly for the other three spinor’s component built with (25).
Such analytical solutions of the radial part of the Dirac equation, computed far from the black hole, are useful to investigate the scattering of charged massive fermions.
3. The Massless Case3.1. The Kerr Metric
In the particular case of massless fermions, the Dirac equation can be solved exactly, with its solutions being given by the Heun confluent functions [20] (fermionic one-particle states in Kerr backgrounds have been considered in [42]).
In view of the analysis presented in the previous section, for μ=0, the system (27) gets the simplified form:
(41)Δ∂r+iK0R1=λR2,Δ∂r−iK0R2=λR1,Dθ−HT1=λT2,Dθ+HT2=−λT1,which firstly leads to the radial Teukolsky equations:
(42)ΔRA′′+r−MRA′+±iΔK0′+K02−λ2RA=0,where A=1,2, the prime denotes the derivative with respect to r, and K0 can be written from (24) putting q=0.
The corresponding solutions can be expressed in terms of Heun confluent functions [18, 19] as
(43)R1=Δ1/4eαz/2z−1γ/2×C1zβ/2HeunCα,β,γ,δ,η,z+C2z−β/2HeunCα,−β,γ,δ,η,z,of variable
(44)z=r−r−r+−r−,where r±=M±M2−a2 are the outer and inner horizons and parameters
(45)α=2iωr+−r−,β=12+2ir+−r−2ωMr−−ma,γ=−12+2ir+−r−2ωMr+−ma,δ=ω4Mω−ir+−r−,η=ω4Mω−ir−−2ω2a2−2a2M2−a2ωa−m22−λ2+38.
For the case under consideration with a<M, the two horizons are real, while for an overspinning Kerr spacetime with a>M, the quantities r+ and r− are complex. The solutions to Heun’s confluent equations are computed as power series expansions around the regular singular point z=0, i.e., r=r−. The series converges for z<1, where the second regular singularity is located. An analytic continuation of the HeunC function is obtained by expanding the solution around the regular singularity z=1 (i.e., r=r+) and overlapping the series.
For the polynomial form of the Heun functions, one has to impose the necessary condition [18, 19]:
(46)δα=−n+1+β+γ2,which gives us the resonant frequencies associated with the massless fermion.
In view of the parameters in (45), it turns out that only the component multiplied by C1 gets a polynomial expression, with the energy ω having the real and imaginary parts given by
(47)ωR=mar−r++r−,ωI=n+12r+−r−2r−r++r−,where m and n are the azimuthal and the principal quantum numbers.
To the first order in a2/M2, the above expressions become
(48)ωR≈ma,ωI≈n+12Ma2,and they depend only on the BH parameters.
Next, for a polynomial which truncates at the order n, once we set the n+1 coefficient in the series expansion to vanish, we get the separation constant λ expressed in terms of the black hole’s parameters.
For the asymptotic behavior at infinity, one may use formula (38) and expression (43) turns into the simplified form:
(49)R1≈Δ1/4rsinωr+2ωM−i2logr2M+σ≈Δ1/4rexpiωr+2ωMlogr2M+σ,where σω is the phase shift.
In order to study the radiation emitted by the black hole, one has to write down the wave function components near the exterior horizon, r→r+. Using (43), for z→1, the (radial) components of Ψ defined in (19), (22), and (25) can be written as
(50)Ψout~e−iωtr−r+ir+−r−2ωMr+−ma.
By definition, the component ψout near the event horizon should asymptotically have the form [24]
(51)Ψout~r−rhi/2κhω−ωh,and the scattering probability at the exterior event horizon surface is given by
(52)Γ=Ψoutr>r+Ψoutr<r+2=exp−2πκhω−ωh.
In our case, using the explicit expressions
(53)κh=r+−r−4Mr+=M2−a22MM+M2−a2,ωh=ma2Mr+,we get the Bose–Einstein distribution for the emitted particles:
(54)N=Γ1−Γ=1λeω/T−1,with
(55)T=κh2π=M2−a24πMM+M2−a256λ=exp−2πωhκh=exp−4πmar+−r−.
One may notice that the expression of the temperature (55) agrees with the one obtained following the usual thermodynamical procedure. Thus, by using the formula of the entropy S=πr+2+a2, with r+=M+M2−a2, and a=J/M, we express the mass in terms of the entropy as
(57)M=r+2+a22r+=S4π+πJ2S1/2and compute the temperature on the event horizon as the following derivative:
(58)T=∂M∂S=S2−4π2J24πS2SS2+4π2J21/2=r+2−a24πr+r+2+a2=M2−a24πMr+.
The corresponding heat capacity at constant angular momentum, i.e.,
(59)CJ=T∂S∂TJ=−2SS4−16π4J4S4−24π2J2S2−48π4J4=2πr+2−a2r+2+a223a4+6r+2a2−r+4,is positive for the following range of the parameter a/M:
(60)23−31/2<aM<1,for which the thermal system is stable on the event horizon.
For a slowly rotating black hole with a/M<23−31/2, the heat capacity becomes negative, corresponding to a thermodynamically unstable phase.
A particular value of a/M where the Kerr–Newman black hole undergoes a phase transition and the heat capacity has an infinite discontinuity was found many years ago by Davies [43].
Secondly, the angular equations coming from the system (27), i.e.,
(61)TA′′+cotθTA′+−cot2θ4−12∓H′−H2+λ2TA=0,where prime means the derivative with respect to θ, for ξ=cosθ, is the spheroidal Teukolsky equation. However, for y=cos2θ/2, the solutions are given by the Heun confluent functions as
(62)T1=eωacosθcosθ2γC1sinθ2βHeunCα,β,γ,δ,η,y+C2sinθ2−βHeunCα,−β,γ,δ,η,y,and similarly for T2, with the real parameters
(63)α=4ωa,β=m+12,γ=m−12,δ=−2ωa,η=1−2mωa−λ2+m22+38.
As expected, for the given parameters of the black hole M,a, the Dirac solutions are enumerated by the half-integer-positive multipole number m±1/2. Since β is not an integer, the two functions in (62) form linearly independent solutions of the confluent Heun differential equation.
Similar expressions have been obtained for the solutions of the Klein-Gordon equation describing a charged massive scalar field in the Kerr–Newman spacetime [24, 44].
Up to a normalization constant A, the first component of Ψ defined in (19), (22), and (25) has the following behavior for large r values:
(64)Ψ1≈Arexpiωr+2ωMlogr2M+σeimφ−ωtT1θ,while the other components can be easily built using the relations (25).
Let us notice that, by introducing the new coordinate r∗=r+2Mlogr/2M, the radial part of the above component has the form obtained by Starobinsky, for the Klein-Gordon equation in the Kerr metric [45], namely,
(65)R~1rAeiωr∗+Be−iωr∗,where A and B are for the incident and reflected wave coefficients, respectively.
3.2. The Extreme Kerr Metric
The extreme Kerr metric can be easily written from (1), by setting the Kerr parameter a equal to M, so that there is a single (degenerate) horizon at r=M with zero Hawking temperature and horizon angular velocity ΩH=1/2M. Thus, for the massless case, the radial equation (42) has the same form, but with Δ=r−M2 and
(66)K0=mM−ωr2+M2r−M.
The solutions are given by the Heun double confluent functions [18, 19] as being
(67)R1~C1exp−iωr−M−ikMr−MHeunDα,β,γ,δ,ζ+C2expiωr−M+ikMr−MHeunD−α,β,γ,δ,ζ,with
(68)k=2Mω−m2M,ζ=r−M+ikM/ωr−M−ikM/ω,and the parameters
(69)α=−8iωkM,β=−32ωMωkM+16ω2M2−4λ2,γ=2α,δ=−32ωMωkM−16ω2M2+4λ2.
Usually, the double confluent Heun functions are obtained from the confluent ones, through an additional confluence process [18, 19].
One may notice that, for ωm=m/2M and r=M, one has to deal with the irregular singularities, at ζ=±1. For ω<m/2M, the variable in (68) is real.
4. Conclusions
Since the pioneering works of Teukolsky [2] and Chandrasekhar [5], the study of the solutions of the massive Dirac equation in the background of an electrically charged black hole has a long history.
The method used in the present paper, while based on Cartan’s formalism with an orthonormal base, is an alternative to the Newman-Penrose (NP) formalism [12], which is usually employed for solving the Dirac equation describing fermions in the vicinity of different types of black holes.
The solutions to the radial Teukolsky equations (42), with two regular singularities at r=r± and an irregular singularity at r=∞, have been written in the form of a series of hypergeometric functions [4]. Similar expressions as the ones in (43) have been found for the exact solutions of the Teukolsky Master Equation for electromagnetic perturbations of the Kerr metric [46] and in the study of bosons in a Kerr-Sen black hole [47].
By imposing the necessary condition for a polynomial form of the Heun confluent functions [18, 19], we get the resonant frequencies, which are of crucial importance for getting information on the black holes interacting with different quantum fields [24].
By identifying the out modes near the r+ horizon, one is able to compute the scattering probability (52) and the Bose–Einstein distribution of the emitted particles. For a=0, we identify the expected Hawking black body radiation and the Hawking temperature Th=1/8πM. It is worth mentioning that, for computing the temperature on the event horizon, we have used the analytical solutions of the Dirac equation, expressed in terms of Heun functions, as an alternative method to the one usually employed in literature. Expression (55) agrees with the one obtained in other works devoted to the thermodynamics of the Kerr–Newman black hole, for example, in [48, 49].
As avenues for further work, it will be interesting to further investigate the connection between the asymptotic resonant frequencies found in this work in the context of the recent research on black hole quantization [39, 50]. In particular, it was shown in those works that the imaginary part of the asymptotic resonant frequencies could further confirm the Bekenstein conjecture [51, 52] of the quantization of the area and entropy of the black holes. We intend to address these issues in further work that will be published elsewhere.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
CarterB.Hamilton-Jacobi and Schrodinger separable solutions of Einstein’s equations196810428031010.1007/BF03399503TeukolskyS. A.Perturbations of a rotating black hole. I. Fundamental equations for gravitational, electromagnetic, and neutrino-field perturbations197318563510.1086/152444LeaverE. W.An analytic representation for the quasi-normal modes of Kerr black holes1985402182328529810.1098/rspa.1985.0119ManoS.SuzukiH.TakasugiE.Analytic solutions of the Teukolsky equation and their low frequency expansions19969561079109610.1143/PTP.95.10792-s2.0-0030499647ChandrasekharS.The solution of Dirac’s equation in Kerr geometry1976349165957157510.1098/rspa.1976.00902-s2.0-0016962863LIGO Scientific and Virgo CollaborationsObservation of gravitational waves from a binary black hole merger20161166, article 06110210.1103/PhysRevLett.116.0611022-s2.0-84958078327LIGO Scientific and VIRGO CollaborationsGW170104: observation of a 50-solar-mass binary black hole coalescence at redshift 0.2201711822, article 22110110.1103/PhysRevLett.118.2211012-s2.0-85020048154LIGO Scientific and Virgo CollaborationsGW151226: observation of gravitational waves from a 22-solar-mass binary black hole coalescence201611624, article 24110310.1103/PhysRevLett.116.2411032-s2.0-84975517093PageD. N.Dirac equation around a charged, rotating black hole19761461509151010.1103/PhysRevD.14.15092-s2.0-0000820433LeeC. H.Massive spin-wave around a Kerr-Newman black hole197768215215610.1016/0370-2693(77)90189-72-s2.0-0009567472KinnersleyW.Field of an arbitrarily accelerating point mass196918651335133610.1103/PhysRev.186.13352-s2.0-0000680538NewmanE.PenroseR.An approach to gravitational radiation by a method of spin coefficients19623356657810.1063/1.17242572-s2.0-36849117833HeX.JingJ.Late-time evolution of charged massive Dirac fields in the Kerr-Newman background20067551-331332810.1016/j.nuclphysb.2006.08.0152-s2.0-33749266918BiniD.CherubiniC.JantzenR. T.RuffiniR. J.Teukolsky Master Equation: de Rham wave equation for gravitational and electromagnetic fields in vacuum2002107596799210.1143/PTP.107.9672-s2.0-0036562884HuangY.LiuD. J.ZhaiX.LiX.Massive charged Dirac fields around Reissner-Nordström black holes: quasibound states and long-lived modes2017966, article 06500210.1103/PhysRevD.96.0650022-s2.0-85031718009SuzukiH.TakasugiE.UmetsuH.Perturbations of Kerr-de Sitter black holes and Heun's equations1998100349150510.1143/PTP.100.4912-s2.0-0032278190SuzukiH.TakasugiE.UmetsuH.Analytic solutions of the Teukolsky equation in Kerr-de Sitter and Kerr-Newman-de Sitter geometries1999102225327210.1143/PTP.102.2532-s2.0-0033243047RonveauxA.1995New York, NY, USAOxford University PressSlavyanovS. Y.LayW.2000Oxford Mathematical Monographs, Oxford University PressSakalliI.HalilsoyM.Solution of the Dirac equation in the near horizon geometry of an extreme Kerr black hole20046912, article 12401210.1103/PhysRevD.69.1240122-s2.0-3843094837ChoH. T.Dirac quasinormal modes in Schwarzschild black hole spacetimes2003682, article 02400310.1103/PhysRevD.68.0240032-s2.0-4944229117KraniotisG. V.The massive Dirac equation in the Kerr-Newman-de Sitter and Kerr-Newman black hole spacetimes201933, article 03502610.1088/2399-6528/ab1046CasalsM.DolanS. R.NolanB. C.OttewillA. C.WinstanleyE.Quantization of fermions on Kerr space-time2013876, article 06402710.1103/PhysRevD.87.0640272-s2.0-84875474515VieiraH. S.BezerraV. B.Confluent Heun functions and the physics of black holes: resonant frequencies, Hawking radiation and scattering of scalar waves2016373284210.1016/j.aop.2016.06.0162-s2.0-84978977367TarloyanA. S.IshkhanyanT. A.IshkhanyanA. M.Four five-parametric and five four-parametric independent confluent Heun potentials for the stationary Klein-Gordon equation20165283-426427110.1002/andp.2015003142-s2.0-84975769685BaticD.SchmidH.Heun equation, Teukolsky equation, and type-D metrics2007484, article 04250210.1063/1.27202772-s2.0-34247616195FizievP. P.Classes of exact solutions to the Teukolsky Master Equation20102713, article 13500110.1088/0264-9381/27/13/1350012-s2.0-77953111062DariescuC.DariescuM. A.SteleaC.The SO(3,1) × U(1)-gauge invariant approach to charged bosons in relativistic magnetars2017491215310.1007/s10714-017-2314-82-s2.0-85034791493DariescuM.DariescuC.SteleaC.Heun-type solutions of the Klein-Gordon and Dirac equations in the Garfinkle-Horowitz-Strominger dilaton black hole background201920198576956410.1155/2019/57695642-s2.0-85062856418DariescuM.-A.DariescuC.SteleaC.The SO(3,1) × U(1)-gauge covariant Dirac equation in relativistic magnetars2018501012610.1007/s10714-018-2449-22-s2.0-85053235430DariescuM.-A.DariescuC.SteleaC.Massless fermions on static general prolate metrics and their Heun solutions2020357, article 205003610.1142/S0217732320500364HortaçsuM.The radial Teukolsky equation for Kerr–Newman–de Sitter geometry: revisited202113611310.1140/epjp/s13360-020-01003-5BirkandanT.HortaçsuM.The zero mass limit of Kerr and Kerr-(anti-)de-Sitter space-times: exact solutions and wormholes20185032810.1007/s10714-018-2351-y2-s2.0-85042331512HortaçsuM.The zero mass limit of Kerr and Kerr-(anti) de Sitter space-times: revisited2020135327910.1140/epjp/s13360-020-00283-1BirkandanT.HortaçsuM.Heun-type solutions for Schwarzschild metric with electromagnetic fields20171192, article 2000210.1209/0295-5075/119/200022-s2.0-85032192461HortacsuM.Heun functions and some of their applications in physics20182018, article 862157311410.1155/2018/86215732-s2.0-85051140930Al-BadawiA.SakalliI.Solution of the Dirac equation in the rotating Bertotti–Robinson spacetime2008495, article 05250110.1063/1.29127252-s2.0-44649197711ÖvgünA.Sakallıİ.SaavedraJ.Quasinormal modes of a Schwarzschild black hole immersed in an electromagnetic universe20184210, article 10510210.1088/1674-1137/42/10/1051022-s2.0-85053230710SakalliI.Analytical solutions in rotating linear dilaton black holes: resonant frequencies, quantization, greybody factor, and Hawking radiation2016948, article 08404010.1103/PhysRevD.94.0840402-s2.0-84992710905Sakallıİ.JusufiK.ÖvgünA.Analytical solutions in a cosmic string Born–Infeld-dilaton black hole geometry: quasinormal modes and quantization2018501012510.1007/s10714-018-2455-42-s2.0-85053261098DolanS. R.GairJ. R.The massive Dirac field on a rotating black hole spacetime: angular solutions20092617, article 17502010.1088/0264-9381/26/17/1750202-s2.0-70349489883LoranF.Fermionic one-particle states in curved spacetimes2018201877110.1007/JHEP07(2018)0712-s2.0-85049884519DaviesP. C. W.The thermodynamic theory of black holes1977353167549952110.1098/rspa.1977.0047KraniotisG. V.The Klein–Gordon–Fock equation in the curved spacetime of the Kerr–Newman (anti) de Sitter black hole20163322, article 22501110.1088/0264-9381/33/22/2250112-s2.0-84994577142StarobinskyA. A.Amplification of waves during reflection from a rotating “black hole”19733712832StaicovaD. R.FizievP. P.The spectrum of electromagnetic jets from Kerr black holes and naked singularities in the Teukolsky perturbation theory2011332238540110.1007/s10509-010-0520-x2-s2.0-79952447115VieiraH. S.BezerraV. B.Charged scalar fields in a Kerr–Sen black hole: exact solutions, Hawking radiation, and resonant frequencies2019433, article 03510210.1088/1674-1137/43/3/0351022-s2.0-85063584516ZhangJ. Y.ZhaoZ.New coordinates for Kerr–Newman black hole radiation20056181-4142210.1016/j.physletb.2005.05.0242-s2.0-20444377288VieiraH. S.BezerraV. B.MunizC. R.Exact solutions of the Klein–Gordon equation in the Kerr–Newman background and Hawking radiation2014350142810.1016/j.aop.2014.07.0112-s2.0-84905250897HodS.Asymptotic quasinormal mode spectrum of rotating black holes2003http://arxiv.org/abs/gr-qc/0307060BekensteinJ. D.The quantum mass spectrum of the Kerr black hole197411946747010.1007/BF027627682-s2.0-51249187920BekensteinJ. D.Black holes and entropy1973782333234610.1103/PhysRevD.7.23332-s2.0-33750485765