Abstract

Most of the theoretical physics known today is described by using a small number of differential equations. For linear systems, different forms of the hypergeometric or the confluent hypergeometric equations often suffice to describe the system studied. These equations have power series solutions with simple relations between consecutive coefficients and/or can be represented in terms of simple integral transforms. If the problem is nonlinear, one often uses one form of the Painlevé equations. There are important examples, however, where one has to use higher order equations. Heun equation is one of these examples, which recently is often encountered in problems in general relativity and astrophysics. Its special and confluent forms take names as Mathieu, Lamé, and Coulomb spheroidal equations. For these equations whenever a power series solution is written, instead of a two-way recursion relation between the coefficients in the series, we find one between three or four different ones. An integral transform solution using simpler functions also is not obtainable. The use of this equation in physics and mathematical literature exploded in the later years, more than doubling the number of papers with these solutions in the last decade, compared to time period since this equation was introduced in 1889 up to 2008. We use SCI data to conclude this statement, which is not precise, but in the correct ballpark. Here this equation will be introduced and examples for its use, especially in general relativity literature, will be given.

1. Introduction

Most of the theoretical physics known today is described by using a small number of differential equations. If we study only linear systems, different forms of the hypergeometric or the confluent hypergeometric equations often suffice to describe the system studied. These equations have power series solutions with simple relations between consecutive coefficients and/or can be represented in terms of simple integral transforms. If the problem is described in terms of nonlinear differential equations, then one often uses one form of the Painlevé equations.

There are important examples, however, where one has to use higher order equations. Such an equation was proposed by Karl Heun in 1889 [1]. This equation and its confluent forms become indispensable in general relativity if one studies exact solutions of wave equations in the background of certain metrics. A well-known example is the Kerr metric [2]. Although it is possible to solve the wave equations in the background of some metrics in terms of hypergeometric functions or its confluent forms, this is not possible for the much studied Kerr metric. If we also study even the trivially extended forms of some metrics by adding a flat dimension to the existing metric, we may have to solve the Heun equation to obtain the exact solution.

Here we will introduce the Heun equation and its confluent forms and mention some of the properties of the Heun equation. Then we will give some examples in physics, mainly in gravitational physics, where one can find many recent papers. This part is meant to be a survey of the work done in the field of General Relativity and Quantum Gravity concentrating on the last decades. In another section we will give an example where the Heun equation emerges from a trivial extension of a wave equation in the background of the Eguchi-Hanson instanton metric [3]. We will end with some concluding remarks.

2. Heun Equation

Let us review some well-known facts about second-order differential equations. Differential equations are classified according to their singularity structure [4, 5]. If a differential equation has no singularities over the full complex plane, it can only be a constant. Singularities are classified as regular singular and irregular singular points. If the coefficient of the first derivative has at most single poles and the coefficient of the term without a derivative has at most double poles when the coefficient of the second derivative is unity, this second-order differential equation has regular singularities, which gives us one regular solution while expanding around this singular point. In general the second solution has a pole or a branch point singularity. If the poles of these coefficients are higher, we have irregular singularities and the general solution has an essential singularity [6].

As stated in Morse and Feshbach [4] an example of a second-order differential equation with one regular singular point is This equation has one solution which is constant. The second solution blows up at infinity. The differential equation has one irregular singularity at infinity which gives an essential singularity at this point. The equation has two regular singular points, at zero and at infinity.

In physics an often used equation is the hypergeometric equation This equation has three regular singular points, at zero, one, and infinity. Jacobi, Legendre, Gegenbauer, and Tchebycheff equations are special forms of this equation. When the singular points at z=1 and z which equal infinity are “coalesced” at infinity, we get the confluent hypergeometric equation with an essential singularity at infinity and a regular singularity at zero. Bessel, Laguerre, and Hermite equations can be reduced to this form.

An important property of all these equations is that they allow infinite series solutions about one of their regular singular points where a recursion relation can be found between two consecutive coefficients. This fact allows one to have an idea about the general properties of the solution, as the asymptotic behaviour at distant points, the radius of convergence of the series, etc.

A new equation was introduced in 1889 by Karl M. W. L. Heun [1]. This is an equation with four regular singular points at zero and one and an arbitrary point f between zero and one and infinity. This equation is discussed in the book edited by Ronveaux [7]. Most of the general information we give below is taken from this book. As discussed there, any equation with four regular singular points can be transformed to the equation given below: There is a relation between the constants given as . This relation is not related to the regularity of the singularity at infinity. It just gives the exponents of the term multiplying the series solution around infinity in terms of as .

If we try to obtain a solution in terms of a power series, one cannot get a recursion relation between two consecutive coefficients. We have a relation at least between three coefficients.

It is known that [8] any second-order differential equation with n regular singular points has a family of local solutions, which splits into 2n sets of equivalent expressions, each set defining one of the two Frobenius solutions in the neighborhood of a singular point. The n! factor comes from permuting the n singular points and the factor from negating exponent differences. Maier [8] gave the list of 192 local solutions for the Heun equation.

The set of transformations that can be applied to the Fuchian equation with n singular points to generate alternative expressions for this equation has order and acts on the parameter space of the equation. This group of transformations is isomorphic to the Coxeter group . These transformations generate solutions. For the Heun case n=4, and this group is isomorphic to , a group of order 192. These transformations will be the combination of Mobius transformations and transformations which multiply the desired solution by powers.

It turns out that the Mobius group PGL(2,C), which takes x to , for nonvanishing AD-BC, can be used where takes values from the different singular points. For Heun equation with four regular singular points, this transformation takes each singular point to five other points, which have zeroes at the same value. These points are given below:,,,

Any one of these transformations maps three of the four points, 0,1,, and infinity, into 0,1, and infinity but generally changes the value of , which takes one of the six possible values: , and . Each value is taken four times.

Just recall the Heun equation: written in terms of the real variable . One writes the solution to the Heun equation in the form This changes the form of the differential equation. For (i) or , (ii) or , and (iii) or , however, the resulting equation has the Heun form. The values given above are the exponents at the singularities [9, 10].

Of course, the parameters of the equations change. For each such combination, say for , there are four possible values s and t can take, namely, both equal to zero; , and , and , and Thus we get three more solutions for each solution. Another factor of six comes from the six different possible values can take. In total for expansions around a single regular singular point, we have twenty-four equivalent solutions, obtained by simply transforming the original equation.

The presence of two different indices for expansion around each singular point doubles the number of equivalent solutions, resulting in 48 solutions for expansions around each singular point. Four singular points multiply this number by four giving the total of 192 local solutions.

It turns out that, for infinite set of values of the parameter q, there are solutions which are analytic at 0 and at 1. These are called Heun functions, whereas those which are analytic only at one point are called local Heun functions [11].

For integer values of one of and for special finite values of , solutions analytic at three singularities exist, the so-called Heun polynomials. A special case is for and , where are eigenvalues of a tridiagonal matrix, we get the solution as a polynomial of degree n, which is analytic at three singular points, 0,1, and [12].

“No example has been given of a solution of Heun’s equation expressed in the form of a definite integral or contour integral involving only functions which are, in some sense, simpler” [13]. This statement does not exclude the possibility of having an infinite series of integrals with “simpler” integrands.

One can obtain different confluent forms of this equation. When we “coalesce” two regular singular points, we get the confluent Heun equation: the standard form of the confluent form equation is given as [14] with solution Another version of this equation can be written as Special forms of this equation are obtained in problems with two Coulombic centers, whose special form, when , is the spheroidal equation Another form is the algebraic form of the Mathieu equation is If we coalesce two regular singular points pairwise, we obtain the double confluent form Here . We can reduce the new equation to the Mathieu equation, an equation with two irregular singularities at zero and at infinity if we reduce this equation to the form Another form is the biconfluent form, where three regular singularities are coalesced. The result is an equation with a regular singularity at zero and an irregular singularity at infinity of higher order: The anharmonic equation in three dimensions can be reduced to this equation: In the triconfluent case, all regular singular points are “coalesced” at infinity which gives the equation below: These different forms are used in different problems in physics.

3. Some Examples of the Heun Equation in Physical Applications

In SCI we found about one hundred thirty papers when Heun functions were searched in the summer of 2010. Now, at the end of April 2018, the number exceeded 330. The number of published articles in SCI more than doubled in the last eight years. More than three fourths of these papers were published in the last ten years. The rest of the papers were published between 1990 and 2005, except a single paper in 1986 [15]. These numbers may differ depending on the institution where one uses the SCI, since different universities in Turkey start their search from different dates. We think we are still in the correct ball park. This shows that although the Heun equation was found in 1889, it was largely neglected in the physics literature until recently. Earlier papers on this topic are mostly articles in mathematics journals. If one looks for books on this topic, published before the year 2000, one finds out the list of books is not very long. There is a book edited by A. Ronveaux, which is a collection of papers presented in the “Centennial Workshop on Heun’s Equations: Theory and Application. Sept. 3–8 1989, Schloss Ringberg.” It was published by the Oxford University Press in 1995 by the title Heun’s Differential Equations [7]. There are two books on functions which are special cases of the Heun Equation: Mathieusche Funktionen und Sphaeroidfunktionen mit anwendungen auf physikalische und technische Probleme by Joseph Meixner and Friedrich Wilhelm Schaefke, published by Springer Verlag in 1954 [16] and a Dover reprint of a book first published in 1946, Theory and Applications of Mathieu Functions by N. W. McLachlan in 1963 [17]. Classical mathematical physics books, such as Morse and Feshbach [4], Whittaker and Watson [18], or the Batemann Manuscript [19], have sections or chapters on the special forms of the Heun equation like Mathieu, Lamé, or spheroidal functions. Some papers on different mathematical properties of these functions can be found in [8, 2025].

A reason why more physicists are interested in the Heun equation recently may be, perhaps, a demonstration of the fact that we do not have simple problems in theoretical physics anymore. Mathematical physicists have to tackle more difficult problems, either with more difficult metrics or in higher dimensions. Both of these extensions may necessitate the use of the Heun functions among the solutions. We can give the Eguchi-Hanson case as an example. The wave equation for the scalar particle in the background of the Eguchi-Hanson metric [3] in four dimensions has hypergeometric functions as solutions [26], whereas the Nutku helicoid [27, 28] metric, the next higher one, gives us Mathieu functions [29], a member of the Heun function set, if the method of separation of variables is used to get a solution. We also find that the scalar particle, in the background of the Eguchi-Hanson metric, trivially extended to five dimensions gives Heun type solutions. [30].

Note that the problem does not need to be very complicated to work with these equations. We encounter Mathieu functions if we consider two-dimensional problems with elliptical shapes [31]. Let us use and , where is the distance from the origin to the focal point. Then the Helmholtz equation can be written as which separates into two equations The solutions to these two equations can be represented as Mathieu and modified Mathieu functions.

If we combine different inverse powers of r, starting from first up to the fourth, or if we combine the quadratic potentials with inverse even powers of two, four and six, we see that the solution of the Schrodinger equation involves Heun functions [23]. Solution to symmetric double Morse potentials also needs these functions, like , where [23]. Similar problems are treated in [3234]

In atomic physics further problems such as separated double wells, Stark effect, and hydrogen molecule ion use these functions. Physics problems which end up with these equations are given in the book by S. Y. Slavyanov and S. Lay [35]. Here we see that even the Stark effect, hydrogen atom in the presence of an external electric field, gives rise to this equation. As described in page 166 of Slavyanov’s book, cited above (original reference is Epstein [36], also treated by S. Yu Slavyanov [37]), when all the relevant constants, namely, Planck constant over , electron mass, and electron charge, are set to unity, the Schrodinger equation for the hydrogen atom in a constant electric field of magnitude in the direction is given by Here is the Laplacian operator. Using parabolic coordinates, where the Cartesian ones are given in terms of the new coordinates by , and and writing the wave function in the product form we get two separated equations: Here and are separation constants that must add to one. We note that these equations are of the biconfluent Heun form.

The hydrogen molecule also is treated in [38]. When the hydrogen-molecule ion is studied in the Born-Oppenheimer approximation, where the ratio of the electron mass to the proton mass is very small, one gets two singly confluent Heun equations if the prolate spheroidal coordinates and are used. Here is the distance between the two centers. Assuming we get two confluent Heun equations: Some additional physics papers with Heun type solutions include the following:

Three relatively recent papers which treat atoms in magnetic fields are as follows.

Exact low-lying states of two interacting equally charged particles in a magnetic field are studied by Truong and Bazzali [39].

The energy spectrum of a charged particle on a sphere under a magnetic field and Coulomb force are studied by Ralko and Truong [40].

B. S. Kandemir presented an analytical analysis of the two-dimensional Schrodinger equation for two interacting electrons subjected to a homogeneous magnetic field and confined by a two-dimensional external parabolic potential. Here a biconfluent Heun (BHE) equation is used [41].

Arda and Sever in one instance with Aydodu studied Schrodinger equation with different potentials and in two cases found Heun and confluent Heun solutions [42, 43].

In two papers Hammann et al. [44, 45] solved the one-dimensional Schrodinger equation for position-dependent masses and obtained Heun solutions. The importance of these papers is the derivation and use of relations between Heun functions which are functions of and , which can be used for obtaining the reflection and transition amplitudes for scattering problems for waves described in terms of Heun functions.

Recently Ishkhanyan showed that the solution of the Schrodinger equation for can be given as a derivative of a triconfluent Heun function [46]. In another paper, solution for the same potential is given [47] as a linear combination of two confluent hypergeometric functions. For another potential which is an inverse square root near the origin and vanishes exponentially at infinity, solution is given in terms of linear combination of Gauss hypergeometric functions [48]. These potentials belong to the Heun class.

Downing showed that the solution to the one-dimensional Schrodinger equation with a hyperbolic double well potential is obtained by a transformation of the confluent Heun equation [49].

Hartmann and Portnoi calculated the bound modes of two-dimensional massless Dirac fermions confined within a hyperbolic secant potential [50]

Portnoi et al. continued studying the two-dimensional Dirac particles in two papers, first confined in nonuniform magnetic fields and second in Poschl-Teller waveguide [51, 52] in terms of confluent Heun functions.

In a relatively recent work P. Dorey [53] showed that equations in finite lattice systems also reduce to Heun equations.

Dislocation movement in crystalline materials and quantum diffusion of kinks along dislocations are some solid state applications of this equation. The book by S. Y. Slavyanov and S. Lay [35] is a general reference on problems solved before 2000.

We also cite a recent mathematical application by A. M. Ishkhanyan et al. where “total fifteen potentials for which the stationary Klein-Gordon equation is solvable in terms of the confluent Heun functions are presented.. Only nine of the potentials are independent due to the transposition symmetry of regular singular points of the equation. Four of these equations can be reduced to the hypergeometric form. The remaining five independent Heun potentials are four-parametric and have solutions only in terms of irreducible confluent Heun functions [54]. Prof Ishkhanyan expands the Heun solution in terms of hypergeometric functions and shows that the sum has only finite number of terms in his cases. Prof. Ishkhanyan wrote additional papers after this one using the same method for other potentials. We will not comment on them, however, since from this point on, we will confine ourselves only to papers on general relativity and cosmology.

Among the papers in general relativity, we also will not be able to comment on all the works of some experts like Prof Fiziev on this field, who wrote scores of papers on Heun equations. We will give only the earlier papers and leave the reader to investigate the later ones in the ArXiv.

In general relativity, in a relatively early work, Teukolsky studied the perturbations of the Kerr metric [55]. If we take for the scalar particle we get two equations. Here is the separation constant, .

Teukolsky just stated these equations [55]. Later these equations were found to be two coupled singly confluent Heun equations [56].

Quasi-normal modes of rotational gravitational singularities were also studied by solving these equations by E.W. Leaver [57].

In recent applications in general relativity, Heun type equations become indispensable when one studies phenomena in higher dimensions, or in different geometries. We must note that even the simplest black hole metric, the Schwarzschild, has solutions in the Heun form [58, 59].

Some other references for general relativity applications are:

D. Batic, H. Schmid, M. Winklmeier where the Dirac equation in the Kerr-Newman metric and static perturbations of the non-extremal Reisner-Nordstrom solution are studied [60]. D. Batic and H. Schmid also studied the Dirac equation for the Kerr-Newman metric and looked for its propagator [61]. They found that the equation satisfied is a form of a general Heun equation described in Reference [60]. In later work Batic, with collaborators continued studying Heun equations and their generalizations [62]. In his most recent paper Batic, with collaborators studied Semi commuting and commuting operators for the Heun family [63].

Prof. P.P. Fiziev studied problems whose solutions are Heun equations extensively.

In a paper published in gr-qc/0603003, he studied the exact solutions of the Regge-Wheeler equation in the Schwarschild black hole interior [58].

He presented a novel derivation of the Teukolsky-Starobinsky identities, based on properties of the confluent Heun functions [64]. These functions define analytically all exact solutions to the Teukolsky master equation, as well as to the Regge-Wheeler and Zerilli ones.

In a talk given at 29th Spanish Relativity Meeting (ERE 2006), he depicted in more detail the exact solutions of Regge-Wheeler equation,which described the axial perturbations of Schwarzschild metric in linear approximation, in the Schwarzschild black hole interior and on Kruskal-Szekeres manifold in terms of the confluent Heun functions [65].

All classes of exact solutions to the Teukolsky master equation were described in terms of confluent Heun functions in Reference [66, 67].

In reference [68] he reveals important properties of the confluent Heun’s functions by deriving a set of novel relations for confluent Heun’s functions and their derivatives of arbitrary order. Specific new sub classes of confluent Heun’s functions are introduced and studied. A new alternative derivation of confluent Heun’s polynomials is presented.

In another paper [69] he, with a collaborator, noted that weak gravitational, electromagnetic, neutrino and scalar fields, considered as perturbations on Kerr background satisfied Teukolsky Master Equation. The two non-trivial equations were obtained after separating the variables, one equation only with the polar angle and another using only the radial variable. These were solved by transforming each one into the form of a confluent Heun equation.

Fiziev is an expert in this topic. Two further articles by him and his collaborator are Solving systems of transcendental equations involving the Heun functions, [70] and Application of the confluent Heun functions for finding the quasinormal modes of non rotating black holes [71].

We also cite one of the last papers of Fiziev on the mathematical properties of this subject which can have applications in physics. In [72], the author “introduces and studies a novel type of solutions to the general Heun equation.” His approach is based “on the symmetric form of the Heun differential equation yielded by development of the Papperitz-Klein symmetric form of the Fuchsian equations with an arbitrary number of regular singular points greater than 4. The symmetry group of these equations turns to be a proper extension of the Mobius group.” He also introduces and studies “new series solutions and derives solutions for the four singular point case which treats simultaneously and on an equal footing all singular points.”

Among other papers on this subject one may cite the following papers.

R. Manvelyan, H. J. W. Muller Kirsten, J. Q. Liang, and Y. Zhang calculated the absorption rate of a scalar by a D3 brane in ten dimensions in terms of modified Mathieu functions and obtained the S-matrix in [73].

T.Oota and Y.Yasui studied the scalar Laplacian on a wide class of five-dimensional toric Sasaki-Einstein manifolds, ending in two Heun’s differential equations in [74].

S. Musiri and G. Siopsis found out that the wave equation, obtained in calculating the asymptotic form of the quasi-normal frequencies for large AdS black holes in five dimensions, reduces to a Heun equation, in [75].

A. Al-Badawi and I. Sakalli studied the Dirac equation in the rotating Bertotti-Robinson space-time [76] ending up with a Heun type equation.

I first encountered this type of equation when we tried to solve the scalar wave in the background of the Nutku helicoid instanton [29]. In this case, for a scalar particle in this background metric, one gets the Mathieu equation which is a special case of the Heun equation. In the same paper, the solutions in four dimensions involve the product of two exponentials and two Heun functions. These solutions can be summed to give Green’s function for this problem in a closed form. We could not succeed in obtaining a closed form solution for the Green function when the similar problem is studied in five dimensions [77, 78].

The helicoid instanton is a double-centered solution. As remarked above, for the simpler instanton solution of Eguchi-Hanson [3] hypergeometric solutions are sufficient [26]. Here one must remark that another paper using the Eguchi-Hanson metric ends up with the confluent Heun equation [79]. These two papers show that sometimes judicious choice of the coordinate system and separation ansatz matters.

Sucu and Ünal also obtained closed solutions for the spinor particle written in the background of the Nutku helicoid instanton [26], whereas using the separation of variables method gives us an infinite series of product of two Mathieu functions [77]

One can show that the solutions of Sucu and Ünal can be expanded in terms of Mathieu functions if one attempts to use the separation of variables method, as described by L.Chaos-Cador and E. Ley-Koo [80].

Tolga Birkandan and I also found an extension of the Heun equation with five singular points [30] and calculated the solution of a scalar field in the background of the Eguchi-Hanson equation trivially extended to five dimensions [30]. Then the solution for the radial component turned out to be given in terms of the confluent Heun equation.

Mirjam Cvetič and Finn Larsen studied grey body factors and event horizons for rotating black holes with two rotation parameters and five charges in five dimensions. When the Klein-Gordon equation for a scalar particle in this background is written, one gets a confluent Heun equation. In the asymptotic region this equation turns into the hypergeometric form [81]. When they studied the similar problem for the rotating black hole with four charges, they again obtained a confluent Heun equation for the radial component of the Klein-Gordon equation, which they reduce to the hypergeometric form by making approximations [82]. These two papers are partly repeated in [83]. The same equations were obtained which were reduced to approximate forms which gave solutions in the hypergeometric form.

M. Cvetič encounters this function in several of her publications and reduces them to the hypergeometric form by giving physical arguments to drop certain terms in the equation. The hypergeometric solution points to the presence of conformal symmetry in the reduced model [84, 85]. The method is going to the extreme and near-extreme (Kerr/CFT correspondence) limits, going to the boundary and in some cases using a “subtracted metric" using a warp factor which preserves all the near-horizon properties of the black hole such as the entropy and the thermodynamic potentials, and if necessary dropping certain terms which are negligible in these limits [8688].

“In general, conformal symmetry does not exist in the non-extremal cases. The solutions often turn out to be of the Heun form. In the extremal case two horizons overlap. In the near extremal case they are very close to each other. In these two cases and in the near horizon limit, we find conformal symmetry, resulting in solutions which are hypergeometric functions, or one of its confluent forms. If we want conformal symmetry without going to the extremal or the near horizon limit, we have to change the ‘warp factor’. When the warp factor is changed, the rest of the metric preserves its initial form. The thermodynamic potentials and entropy do not change. You have to drop some terms resulting in solutions in the hypergeometric form. This is equivalent to putting the black hole into a conic box. If you go to the asymptotic or to the scaling limit, this is seen clearly. In these limits the Einstein equations are not satisfied unless the energy-momentum tensor, on the right side of the Einstein equations are also changed, to account for putting the system into the conic box.” [89]

Cvetič also studied black holes in supergravity with Birkandan. Heun solutions also exist for the Wu Black Hole which is the most general solution of maximally supersymmetric gauged supergravity in D=5 [90]. Here they did not study the limiting cases. For the massless Klein-Gordon equation in the background of the most general black hole in four dimensions and N=2 gauge supersymmetry with gauge symmetry (Chow-Compere solution [91]), the angular equation gives Heun type solutions. The radial equation has five regular singularities, which reduce to hypergeometric functions in the near-horizon extremal limit [92].

We should also mention two papers by H. R. Christiansen and M. S. Cunha with Heun type solutions. These are Confluent Heun Functions in Gauge Theories on Thick Braneworlds [93] and Kalb-Ramond Excitations in a Thick-Brane Scnario with Dilaton [94]. In the first paper, the propagation modes of gauge fields in an infinite Randall-Sundrum scenario are investigated. Here a sine-Gordon soliton represents the thick four-dimensional braneworld while an exponentially coupled scalar field acts for the dilaton. For the gauge field motion a differential equation is found which can be transformed into a confluent Heun equation. In the second paper a similar scenario is used. Here a bulk Kalb-Ramond field is coupled to a dilaton, in a warped space-time in the presence of a brane field in five dimensions. Full spectrum and eigenstates are studied. In the general case, the solutions to the field equations are given in terms of the confluent Heun function, which reduces to the confluent hypergeometric function for special values of the parameters.

Other relevant references I could find are listed as [95102].

The more recent papers on this subject include The Quantum Treatment of the 5D-Warped Friedman-Robertson-Walker Universe in Schrodinger Picture [103]. Here the time-evolving Schrodinger version of the Wheeler-De Witt equation, written for the five-dimensional warped k=0-FRW Universe, is studied. For small values of the cosmological scale factor, , the wave function of the Universe is expressed in terms of the Heun Double Confluent functions, whereas for large values of this parameter the solution becomes the Hermite associated functions. Two papers by the same authors using Heun type functions are Fermions in Magnestar’s Crust in terms of Heun Double Confluent Functions [104] and The Approximative Analytic Study of Fermions in Magnetar’s Crust; Ultra-Relativistic Plane Waves, Heun and Mathieu Solutions and Beyond [105].

In Fermi Surfaces and Analytic Green’s Functions from Conformal Gravity [106], T2-symmetric charged AdS black holes are constructed in conformal gravity. The most general solution up to an overall conformal factor contains three nontrivial parameters: the mass, electric charge, and a quantity that can be identified as the massive spin-2 hair. The Dirac equation for the charged massless spinor in this background can be solved in terms of the general Heun function for generic frequency and wave number . This allows us to obtain the analytic Green function for both extremal and nonextremal black holes. For some special choice of black hole parameters, the Green function reduces to simpler hypergeometric or confluent hypergeometric functions.

Two of the authors of the paper quoted above had calculated Green’s functions in terms of the Heun function in an earlier paper, Exact Green’s Functions from Conformal Gravity [107].

Another paper is Quantized Black Hole and Heun Function by D. Momeni, K. Yerzhanov, and R. Myrzakulov [108] where a black hole is quantized using the Bohr method. The solution turns to be of the Heun type.

In On an Approach to Constructing Static Ball Models in General Relativity by A. M. Baranov, some solutions of the Einstein equation were described by Heun functions [109].

In an paper on analytic solutions of wave equations in regular coordinate systems on Schwarzschild background Dennis Philipp and Volker Perlick claim that they find “the wave equation for the propagation of (massless) scalar, electromagnetic and gravitational waves on fixed Schwarzschild background spacetime, which is described by the general time-dependent Regge-Wheeler equation, can be transformed to usual Schwarzschild, Eddington-Finkelstein, Painleve Gullstrand and Kruskal-Szekeres coordinates. In the first three cases, but not in the last one, it is possible to separate a harmonic time-dependence. Then the resulting radial equations belong to the class of confluent Heun equations” [59].

Among additional papers we can also cite the article of Bezerra et al., Exact Solutions of the Klein-Gordon Equation in the Kerr-Newman Background and Hawking Radiation, where both the radial and angular solutions are given in terms of confluent Heun functions [110]. In the particular case corresponding to an extreme Kerr-Newman black hole, the solution is given by the double confluent Heun functions [111]. Biconfluent Heun functions were obtained for the exact solution of the Schrodinger equation for a particle (galaxy) moving in a Newtonian universe with a cosmological constant [112].

Other papers on general relativity written in 2015 also include New Results for Electromagnetic Quasinormal and Quasibound Modes of Kerr Black Holes, by D. Staicova and P. Fiziev [113], where the authors solve Teukolsky equations with confluent Heun solutions numerically. In Heun Functions Describing Fermions Evolving in Parallel and Magnetic Fields, by C. Dariescu and M. A. Dariescu [114], the solutions are in terms of double confluent Heun functions. The same authors also published Quantum Analysis of k=-1 Robert-Walker Universe, where they solved the Wheeler-DeWitt equation [115]. The solutions turned out to be Heun functions. M. C. E. Cedeno and J. C. N. Araujo show that, for Master equation solutions in the linear regime of characteristic formulation of general relativity, the solution is in terms of confluent Heun’s functions for radiative case in the Schwarzschild’s background [116]. In Massless Dirac Particles in the Vacuum C-Metric, D. Bini et al. show that the Dirac equation, written in the background of the C-metric, can be reduced to a radial and an angular equation, both of which can be solved in terms of general Heun functions [117]. Vieira et al. [118] show that for Charged Massive Scalar Fields are Considered in the Gravitational and Electromagnetic Field Produced by a Dyonic Black Hole with a Cosmic String along its Axis of Symmetry “exact solutions of both angular and radial parts of the covariant Klein-Gordon equation in this background can be obtained, and are given in terms of the confluent Heun functions.” In [119], Kofron separates test fields equations on the nonrotating C-metric background. He finds that the resulting equations are of the Heun or confluent Heun form for the general case. These equations, however, can be reduced to hypergeometric functions in the static, axial symmetric, and the extremal case where the inner and outer horizons coalesce. In another paper [120], the same author studies the similar phenomena on the background of the rotating C-metric. For the general case, the radial equation has five regular singularities. In the extremal, static, and axial symmetric cases, one obtains a polynomial solution.

Some other papers published in 2016 in the field of general relativity where solutions to field equations are in the background of different metrics are as follows.

Valtancoli [121] found Heun solutions for the radial part of the Klein-Gordon equation when the scalar field is conformally coupled to a charged BTZ black hole.

Vieira and Bezerra [122] study “resonant frequencies, Hawking radiation and scattering of scalar waves” and find confluent Heun solutions. They also study [123] the class of solutions of the Wheeler-DeWitt equation in the Friedmann-Robertson-Walker universe. In still another paper [124], these authors find confluent Heun solutions for the massless Klein-Gordon equation in the background metric of the three-dimensional rotating and four-dimensional canonical acoustic black holes.

Sakalli [125] finds analytical solutions in rotating linear dilaton black holes.

Kraniotis [126] studies the Klein-Gordon equation in the background metric of the Kerr-Newman (anti-)de Sitter black hole. He first reduces the radial and angular equations to the Heun form and writes the solution in terms of local Heun and confluent Heun functions. In my opinion this paper should be also praised for the introduction of the “false singular point” concept, which reduces the solution to hypergeometric functions for certain values of the physical parameters in the equation.

Since we updated this paper in February 2017, we find close to thirty new publications if one searches for the word Heun Functions in the index Web of Knowledge in April 2018. Many of these papers are on the mathematical aspects of the equation and solving Schrodinger equations for different new potentials in terms of Heun or linear combinations of Heun functions. There are also solutions in terms of Heun functions for equations used in different branches of physics. Here we will attempt to review only the papers for applications in physics related to general relativity.

In [127], Arda et al. solve the energy relations obtained with the help of the quantization rule for the Klein-Gordon equation with a linear plus an inverse-linear potential in terms of biconfluent Heun equations. Vieira wrote two papers [128, 129] where he first studied Resonant Frequencies of a Hydrodynamic Vortex. The radial equation has solutions in terms of double confluent Heun functions. In the second paper, analytic solutions for sound perturbations in the presence of a rotating acoustic black hole which is an analogue of the conical Kerr metric were studied. In the massless case, the radial equation has Heun type solutions. Vieira also wrote another paper with coauthors [130], where massive scalar fields are considered in the gravitational field produced by a Schwarzschild black hole with a global monopole in f(R) gravity. The exact solution of the radial part of the Klein-Gordon equation in this background is given in terms of the general Heun functions. The properties of the general Heun functions are applied to study the Hawking radiation and the resonant frequencies of scalar particles.

Ciprian Dariescu wrote two papers with collaborators [131, 132]. In the first paper, using a perturbative method, Klein-Gordon equation for a charged massive field in the background of a magnetar is solved both in the interior solution and outside the star. Equations can be separated with general and confluent Heun function solutions. With special conditions on parameters, polynomial solutions can be found and first-order transition amplitudes are computed [131]. In the second paper, for the spatially open Friedmann-Robertson-Walker (FRW) Universe with stiff matter and radiation as noninteracting matter sources, the scale function coming from the integration of the Friedmann equation is expressed in terms of elliptic integrals. For a negative cosmological constant, the allowed ranges for the models parameters are identified. Within the quantum analysis, the Wheeler–DeWitt (WDW) equation turns into a modified Morse equation whose solutions are Mathieu and Heun functions. [132].

Sobhani et al. [133] wrote a paper where the thermodynamic properties of the anharmonic oscillator cosmic string framework are studied. The Schrodinger equation is written in the cosmic string framework and anharmonic oscillations are investigated. The wave function and energy spectrum are derived using confluent Heun functions.

Birkandan was also active in this period. He wrote four papers. In the first paper, with Bouaziz, he showed that the deformed Schrodinger equation for a singular inverse square potential in coordinate space with a minimal length is solved in terms of Heun functions [134]. In his second paper with a collaborator, confluent Heun solutions to the radial equations of two Halilsoy-Badawi metrics are found. For the first metric, the radial part of the massless Dirac equation and, for the second case, the radial part of the massless Klein-Gordon equation are studied [135], both with Heun type solutions. In the third paper, he and his collaborator showed that Heun-type exact solutions emerged for both the radial and the angular equations for the case of a scalar particle coupled to the zero-mass limit of both the Kerr and Kerr-(anti-)de-Sitter space-times. Since any type D metric has Heun-type solutions, it is interesting that this property is retained when the black hole has a zero-mass limit. This work further refuted the claims that mass of the black hole, going to zero limit of the Kerr metric, was both locally and globally the same as the Minkowski metric [136]. We comment on the fourth paper in the Conclusion.

A comprehensive bibliography can be found at the bibliography section of http://tcpa.uni-sofia.bg/heun/home.html, compiled by Professors Plamen Fiziev and Denitsa Staicova.

Just to give an example of how the Heun function emerges in a simple problem, in the next section, our work in [30] for the scalar particle in the background metric of the extended Eguchi-Hanson solution will be sketched.

4. Scalar Field in the Background of the Extended Eguchi-Hanson Solution

To go to five dimensions, we can add a time component to the Eguchi-Hanson [3] metric so that we havewhereThis is a vacuum solution. If we takewe find the scalar equation asIf we take , the solution of the radial part is expressed in terms of confluent Heun () functions. If the variable transformation is made, one solution can be expressed asWe tried to express the equation for the radial part in terms of to see the singularity structure more clearly. Then the radial differential operator readsThis operator has two regular singularities at zero and one and an irregular singularity at infinity, the singularity structure of the confluent Heun equation. This is different from the hypergeometric equation, which has regular singularities at zero, one and infinity.

The solution of the angular equation which is regular at for greater than is given below in terms of hypergeometric functions.

5. Conclusion

In this paper, first the Heun function is introduced; then some of its uses in physics, especially in the field of general relativity and gravitation, are demonstrated. We have to note that most of the physicists that bluntly state their solution is in terms of Heun functions are mainly from the third world. We see physicists from Bulgaria, Romania, Brazil, Armenia, and even Turkey in this group. There are mathematicians from the western world, though, who are experts in this field. Batic, a mathematician, although he now works in UAE may be considered from the western world. Ronveaux from Belgium and many other mathematicians are from the western world.

They are not really many exceptions to this observation. Cvetič and Larsen demonstrate what the physicists from the western world do. They try to express their solutions in terms of hypergeometric functions, by going to the asymptotics, to the extremal, or to the near-extremal limit, or putting the solution into a conic box, by changing the energy momentum term if necessary, but keeping the thermodynamic potentials the same. A long endeavor was necessary to label Teukolsky Master Equations as belonging to the Heun class [56]. Only recently was the equation given by ’t Hooft [137] shown to belong to the Heun class if it were not modified [138]. When modified the solution is the manageable hypergeometric function. We agree that this impression may be wrong, but it is just an observation.

The first version of this paper was submitted to the 13th Regional Conference on Mathematical Physics, which was held in Antalya, Turkey, on 27–31 October 2010 and printed in [139].

Disclosure

This paper is a revised and many-times-updated version of the conference talk by the same author, published in Proceedings of the 13th Regional Conference on Mathematical Physics, Antalya, Turkey, 27–31 October 2010, edited by Uğur Camcı and Ibrahim Semiz, pp. 23–39, World Scientific (2013).

Conflicts of Interest

The author declares that they have no conflicts of interest.

Acknowledgments

The author is grateful to Professors Cemsinan Deliduman and Kayhan for providing the author with a shelter at Mimar Sinan Fine Arts University during the author’s days in retirement. The author is indebted to Tolga Birkandan for collaboration and technical assistance and is grateful to Professor Dr. André Ronveaux for informing the author of a slight error in [9]. The author thanks Science Academy, Turkey, for support.