Advances in High Energy Physics
Volume 2021 (2021), Article ID 5512735, 10 pages
https://doi.org/10.1155/2021/5512735

Research Article

Dirac Equation on the Kerr–Newman Spacetime and Heun Functions

Ciprian Dariescu,¹ Marina-Aura Dariescu,¹ and Cristian Stelea²

¹Faculty of Physics, “Alexandru Ioan Cuza” University of Iasi, 11 Bd. Carol I, Iasi 700506, Romania
²Science Research Department, Institute of Interdisciplinary Research, “Alexandru Ioan Cuza” University of Iasi, 11 Bd. Carol I, Iasi 700506, Romania

Correspondence should be addressed to Cristian Stelea; cristian.stelea@uaic.ro

Received 12 February 2021; Revised 3 March 2021; Accepted 2 April 2021; Published 14 April 2021

Academic Editor: Hasan El moumni

Copyright © 2021 Ciprian Dariescu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP³.

Abstract.

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 -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 -Gauge Covariant Dirac Equation

Let us start with the four-dimensional Kerr–Newman metric in the usual Boyer–Lindquist coordinates:

where , , and , , and 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:

Within a -gauge covariant formulation, we introduce the pseudoorthonormal frame , whose corresponding dual base is

leading to the expressions

Thus, using the relations , i.e., , one may write down the pseudoorthonormal frame:

Using (3), the first Cartan’s equation

with and , can be explicitely worked out as

where and are the derivatives with respect to and , leading to the following complete list of nonzero connection coefficients in the Cartan frames :

Now, one has all the essentials to write down the -gauge covariant Dirac equation for the fermion of mass :

where “” stands for the covariant derivative:

with .

In view of the relations (8), the term expressing the Ricci spin connection

has the concrete expression

where , while the kinetic term reads

Putting everything together, the Dirac equation (9) has the explicit form:

where the proper four-potential component, coming from , with given in (2), reads

For ease of calculations, the choice for matrices is important and we are going to employ Weyl’s representation:

with

so that

Thus, for the bispinor written in terms of two-component spinors as

the general equation (14) leads to the following system of coupled equations for the spinors and :

where and .

Due to the time independence and symmetry of the spacetime, we can assume that the wave function can be written as

where the factors have been introduced in order to pull some terms out of equations (20) and (21).

With the new functions and , equations (20) and (21) can be put into the transparent form:

where we have introduced the operators:

Finally, by applying the separation ansatz

one gets the system:

which leads to the radial and angular equations:

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 component:

and , for .

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 , equation (28), with given in (24), reads

with the notation

where one may identify the fermion’s quanta energy, , the standard Coulomb energy, , and the internal centrifugal energy with the quantum resonant correction, i.e., .

To the first order in , meaning a slowly rotating object, for which

the solutions of (29) are given in terms of the Heun confluent functions [18, 19] as

with the parameters written in the physical transparent form as

where and is the energy computed on the Schwarzschild horizon, i.e.,

The second component, , is given by the complex conjugated expression of (32).

One may notice that, for and , the first component in defined in (22) reads

with .

Moreover, since , the modal radial current (of quantum origin), computed as

vanishes. The only nonvanishing component is the azimuthal one, which is given by the expression

The current has the generic representation given in Figure 1, for , i.e., . 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 of the slowly rotating black hole, where the Heun functions have a regular singularity.

Figure 1: The radial part of the current (37), for .

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]

so that the two independent solutions in (32) are given by the simple expression

where is the phase shift, , and .

Thus, the first component of defined in (19), (22), and (25) has the following (physical) behavior for large values:

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 Case

3.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 , the system (27) gets the simplified form:

which firstly leads to the radial Teukolsky equations:

where , the prime denotes the derivative with respect to , and can be written from (24) putting .

The corresponding solutions can be expressed in terms of Heun confluent functions [18, 19] as

of variable

where are the outer and inner horizons and parameters

For the case under consideration with , the two horizons are real, while for an overspinning Kerr spacetime with , the quantities and are complex. The solutions to Heun’s confluent equations are computed as power series expansions around the regular singular point , i.e., . The series converges for , where the second regular singularity is located. An analytic continuation of the HeunC function is obtained by expanding the solution around the regular singularity (i.e., ) and overlapping the series.

For the polynomial form of the Heun functions, one has to impose the necessary condition [18, 19]:

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 gets a polynomial expression, with the energy having the real and imaginary parts given by

where and are the azimuthal and the principal quantum numbers.

To the first order in , the above expressions become

and they depend only on the BH parameters.

Next, for a polynomial which truncates at the order , once we set the 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:

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, . Using (43), for , the (radial) components of defined in (19), (22), and (25) can be written as

By definition, the component near the event horizon should asymptotically have the form [24]

and the scattering probability at the exterior event horizon surface is given by

In our case, using the explicit expressions

we get the Bose–Einstein distribution for the emitted particles:

with

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 , with , and , we express the mass in terms of the entropy as

and compute the temperature on the event horizon as the following derivative:

The corresponding heat capacity at constant angular momentum, i.e.,

is positive for the following range of the parameter :

for which the thermal system is stable on the event horizon.

For a slowly rotating black hole with , the heat capacity becomes negative, corresponding to a thermodynamically unstable phase.

A particular value of 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.,

where prime means the derivative with respect to , for , is the spheroidal Teukolsky equation. However, for , the solutions are given by the Heun confluent functions as

and similarly for , with the real parameters

As expected, for the given parameters of the black hole , the Dirac solutions are enumerated by the half-integer-positive multipole number . 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 , the first component of defined in (19), (22), and (25) has the following behavior for large values:

while the other components can be easily built using the relations (25).

Let us notice that, by introducing the new coordinate , the radial part of the above component has the form obtained by Starobinsky, for the Klein-Gordon equation in the Kerr metric [45], namely,

where and 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 equal to , so that there is a single (degenerate) horizon at with zero Hawking temperature and horizon angular velocity . Thus, for the massless case, the radial equation (42) has the same form, but with and

The solutions are given by the Heun double confluent functions [18, 19] as being

with

and the parameters

Usually, the double confluent Heun functions are obtained from the confluent ones, through an additional confluence process [18, 19].

One may notice that, for and , one has to deal with the irregular singularities, at . For , 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 and an irregular singularity at , 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 horizon, one is able to compute the scattering probability (52) and the Bose–Einstein distribution of the emitted particles. For , we identify the expected Hawking black body radiation and the Hawking temperature . 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.

References

B. Carter, “Hamilton-Jacobi and Schrodinger separable solutions of Einstein’s equations,” Communications in Mathematical Physics, vol. 10, no. 4, pp. 280–310, 1968.

S. A. Teukolsky, “Perturbations of a rotating black hole. I. Fundamental equations for gravitational, electromagnetic, and neutrino-field perturbations,” The Astrophysical Journal, vol. 185, p. 635, 1973.

E. W. Leaver, “An analytic representation for the quasi-normal modes of Kerr black holes,” Proceeding of the Royal Society of London, vol. 402, no. 1823, pp. 285–298, 1985.

S. Mano, H. Suzuki, and E. Takasugi, “Analytic solutions of the Teukolsky equation and their low frequency expansions,” Progress of Theoretical Physics, vol. 95, no. 6, pp. 1079–1096, 1996.

S. Chandrasekhar, “The solution of Dirac’s equation in Kerr geometry,” Proceeding of the Royal Society of London, vol. 349, no. 1659, pp. 571–575, 1976.

LIGO Scientific and Virgo Collaborations, “Observation of gravitational waves from a binary black hole merger,” Physical Review Letters, vol. 116, no. 6, article 061102, 2016.

LIGO Scientific and VIRGO Collaborations, “GW170104: observation of a 50-solar-mass binary black hole coalescence at redshift 0.2,” Physical Review Letters, vol. 118, no. 22, article 221101, 2017.

LIGO Scientific and Virgo Collaborations, “GW151226: observation of gravitational waves from a 22-solar-mass binary black hole coalescence,” Physical Review Letters, vol. 116, no. 24, article 241103, 2016.

D. N. Page, “Dirac equation around a charged, rotating black hole,” Physical Review D, vol. 14, no. 6, pp. 1509-1510, 1976.

C. H. Lee, “Massive spin-wave around a Kerr-Newman black hole,” Physics Letters, vol. 68, no. 2, pp. 152–156, 1977.

W. Kinnersley, “Field of an arbitrarily accelerating point mass,” Physics Review, vol. 186, no. 5, pp. 1335-1336, 1969.

E. Newman and R. Penrose, “An approach to gravitational radiation by a method of spin coefficients,” Journal of Mathematical Physics, vol. 3, no. 3, pp. 566–578, 1962.

X. He and J. Jing, “Late-time evolution of charged massive Dirac fields in the Kerr-Newman background,” Nuclear Physics B, vol. 755, no. 1-3, pp. 313–328, 2006.

D. Bini, C. Cherubini, R. T. Jantzen, and R. J. Ruffini, “Teukolsky Master Equation: de Rham wave equation for gravitational and electromagnetic fields in vacuum,” Progress of Theoretical Physics, vol. 107, no. 5, pp. 967–992, 2002.

Y. Huang, D. J. Liu, X. Zhai, and X. Li, “Massive charged Dirac fields around Reissner-Nordström black holes: quasibound states and long-lived modes,” Physical Review D, vol. 96, no. 6, article 065002, 2017.

H. Suzuki, E. Takasugi, and H. Umetsu, “Perturbations of Kerr-de Sitter black holes and Heun's equations,” Progress of Theoretical Physics, vol. 100, no. 3, pp. 491–505, 1998.

H. Suzuki, E. Takasugi, and H. Umetsu, “Analytic solutions of the Teukolsky equation in Kerr-de Sitter and Kerr-Newman-de Sitter geometries,” Progress of Theoretical Physics, vol. 102, no. 2, pp. 253–272, 1999.

A. Ronveaux, Ed., Heun’s Differential Equations, Oxford University Press, New York, NY, USA, 1995.

S. Y. Slavyanov and W. Lay, Special Functions, A Unified Theory Based on Singularities, Oxford Mathematical Monographs, Oxford University Press, 2000.

I. Sakalli and M. Halilsoy, “Solution of the Dirac equation in the near horizon geometry of an extreme Kerr black hole,” Physical Review D, vol. 69, no. 12, article 124012, 2004.

H. T. Cho, “Dirac quasinormal modes in Schwarzschild black hole spacetimes,” Physical Review D, vol. 68, no. 2, article 024003, 2003.

G. V. Kraniotis, “The massive Dirac equation in the Kerr-Newman-de Sitter and Kerr-Newman black hole spacetimes,” Journal of Physics Communications, vol. 3, no. 3, article 035026, 2019.

M. Casals, S. R. Dolan, B. C. Nolan, A. C. Ottewill, and E. Winstanley, “Quantization of fermions on Kerr space-time,” Physical Review D, vol. 87, no. 6, article 064027, 2013.

H. S. Vieira and V. B. Bezerra, “Confluent Heun functions and the physics of black holes: resonant frequencies, Hawking radiation and scattering of scalar waves,” Annals of Physics, vol. 373, pp. 28–42, 2016.

A. S. Tarloyan, T. A. Ishkhanyan, and A. M. Ishkhanyan, “Four five-parametric and five four-parametric independent confluent Heun potentials for the stationary Klein-Gordon equation,” Annalen der Physik, vol. 528, no. 3-4, pp. 264–271, 2016.

D. Batic and H. Schmid, “Heun equation, Teukolsky equation, and type-D metrics,” Journal of Mathematical Physics, vol. 48, no. 4, article 042502, 2007.

P. P. Fiziev, “Classes of exact solutions to the Teukolsky Master Equation,” Classical and Quantum Gravity, vol. 27, no. 13, article 135001, 2010.

C. Dariescu, M. A. Dariescu, and C. Stelea, “The SO(3,1) × U(1)-gauge invariant approach to charged bosons in relativistic magnetars,” General Relativity and Gravitation, vol. 49, no. 12, p. 153, 2017.

M. Dariescu, C. Dariescu, and C. Stelea, “Heun-type solutions of the Klein-Gordon and Dirac equations in the Garfinkle-Horowitz-Strominger dilaton black hole background,” Advances in High Energy Physics, vol. 2019, Article ID 5769564, 8 pages, 2019.

M.-A. Dariescu, C. Dariescu, and C. Stelea, “The SO(3,1) × U(1)-gauge covariant Dirac equation in relativistic magnetars,” General Relativity and Gravitation, vol. 50, no. 10, p. 126, 2018.

M.-A. Dariescu, C. Dariescu, and C. Stelea, “Massless fermions on static general prolate metrics and their Heun solutions,” Modern Physics Letters A, vol. 35, no. 7, article 2050036, 2020.

M. Hortaçsu, “The radial Teukolsky equation for Kerr–Newman–de Sitter geometry: revisited,” The European Physical Journal Plus, vol. 136, no. 1, p. 13, 2021.

T. Birkandan and M. Hortaçsu, “The zero mass limit of Kerr and Kerr-(anti-)de-Sitter space-times: exact solutions and wormholes,” General Relativity and Gravitation, vol. 50, no. 3, p. 28, 2018.

M. Hortaçsu, “The zero mass limit of Kerr and Kerr-(anti) de Sitter space-times: revisited,” The European Physical Journal Plus, vol. 135, no. 3, p. 279, 2020.

T. Birkandan and M. Hortaçsu, “Heun-type solutions for Schwarzschild metric with electromagnetic fields,” EPL, vol. 119, no. 2, article 20002, 2017.

M. Hortacsu, “Heun functions and some of their applications in physics,” The European Physical Journal Plus, vol. 2018, article 8621573, pp. 1–14, 2018.

A. Al-Badawi and I. Sakalli, “Solution of the Dirac equation in the rotating Bertotti–Robinson spacetime,” Journal of Mathematical Physics, vol. 49, no. 5, article 052501, 2008.

A. Övgün, İ. Sakallı, and J. Saavedra, “Quasinormal modes of a Schwarzschild black hole immersed in an electromagnetic universe,” Chinese Physics C, vol. 42, no. 10, article 105102, 2018.

I. Sakalli, “Analytical solutions in rotating linear dilaton black holes: resonant frequencies, quantization, greybody factor, and Hawking radiation,” Physical Review D, vol. 94, no. 8, article 084040, 2016.

İ. Sakallı, K. Jusufi, and A. Övgün, “Analytical solutions in a cosmic string Born–Infeld-dilaton black hole geometry: quasinormal modes and quantization,” General Relativity and Gravitation, vol. 50, no. 10, p. 125, 2018.

S. R. Dolan and J. R. Gair, “The massive Dirac field on a rotating black hole spacetime: angular solutions,” Classical and Quantum Gravity, vol. 26, no. 17, article 175020, 2009.

F. Loran, “Fermionic one-particle states in curved spacetimes,” Journal of High Energy Physics, vol. 2018, no. 7, p. 71, 2018.

P. C. W. Davies, “The thermodynamic theory of black holes,” Proceeding of the Royal Society of London, vol. 353, no. 1675, pp. 499–521, 1977.

G. V. Kraniotis, “The Klein–Gordon–Fock equation in the curved spacetime of the Kerr–Newman (anti) de Sitter black hole,” Classical and Quantum Gravity, vol. 33, no. 22, article 225011, 2016.

A. A. Starobinsky, “Amplification of waves during reflection from a rotating “black hole”,” Journal of Experimental and Theoretical Physics, vol. 37, no. 1, pp. 28–32, 1973.

D. R. Staicova and P. P. Fiziev, “The spectrum of electromagnetic jets from Kerr black holes and naked singularities in the Teukolsky perturbation theory,” Astrophysics and Space Science, vol. 332, no. 2, pp. 385–401, 2011.

H. S. Vieira and V. B. Bezerra, “Charged scalar fields in a Kerr–Sen black hole: exact solutions, Hawking radiation, and resonant frequencies,” Chinese Physics C, vol. 43, no. 3, article 035102, 2019.

J. Y. Zhang and Z. Zhao, “New coordinates for Kerr–Newman black hole radiation,” Physics Letters B, vol. 618, no. 1-4, pp. 14–22, 2005.

H. S. Vieira, V. B. Bezerra, and C. R. Muniz, “Exact solutions of the Klein–Gordon equation in the Kerr–Newman background and Hawking radiation,” Annals of Physics, vol. 350, pp. 14–28, 2014.

S. Hod, “Asymptotic quasinormal mode spectrum of rotating black holes,” 2003, http://arxiv.org/abs/gr-qc/0307060.

J. D. Bekenstein, “The quantum mass spectrum of the Kerr black hole,” Lettere Al Nuovo Cimento Series 2, vol. 11, no. 9, pp. 467–470, 1974.

J. D. Bekenstein, “Black holes and entropy,” Physical Review D, vol. 7, no. 8, pp. 2333–2346, 1973.