Heun Functions and their uses in Physics

3 downloads 0 Views 181KB Size Report
Jan 3, 2011 - tion methods. If one studies certain problems in astronomy or general ... These equations have power series solutions with simple relations ...
arXiv:1101.0471v1 [math-ph] 3 Jan 2011

Heun Functions and their uses in Physics M. Horta¸csu Istanbul Technical University, Department of Physics, Istanbul, Turkey January 4, 2011 Abstract 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 this problem. These equations have power series solutions with simple relations between consecutive coefficients and can be generally represented in terms of simple integral transforms. If the problem is nonlinear, one often uses one form of the Painlev´e equation. There are important examples, however, where one has to use more complicated equations. An example often encountered in quantum mechanics is the hydrogen atom in an external electric field, the Stark effect. One often bypasses this difficulty by studying this problem using perturbation methods. If one studies certain problems in astronomy or general relativity, encounter with Heun equation is inevitable. This is a general equation whose special forms take names as Mathieu, Lam´e and Coulomb spheroidal equations. Here the coefficients in a power series expansions do not have two way recursion relations. We have a relation at least between three or four different coefficients. A simple integral transform solution also is not obtainable. Here I will try to introduce this equation whose popularity increased recently, mostly among theoretical physicists,and give some examples where the result can be expressed in terms of solutions of this equation.

1

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 this problem. These equations have power series solutions with simple relations between consecutive coefficients and can be generally represented in terms of simple integral transforms. If the problem is nonlinear, one often uses one form of the Painlev´e equation. Let us review some well known facts about linear second order differential equations. Differential equations are classified according to their singularity structure [1, 2]. If a differential equation has no singularities over the full complex plane, it can only be a constant. Singularities are classified as regular or 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. Then we have one regular solution while expanding around these singular points. The second solution around a regular singular point has a branch cut. If the poles of these coefficients are higher, we have irregular singularities and the general solution has essential singularities around these points [3]. As stated in Morse and Feshbach [4], an example of a second order differential equation with one regular singular point is d2 w = 0. (1) dz 2 This equation has one solution which is constant. The second solution blows up at infinity. The differential equation d2 w + k 2 w = 0. (2) dz 2 has one irregular singularity at infinity which gives an essential singularity at this point. The equation dw d2 w + (1 + a) = 0. (3) 2 dz dz has two regular singular points, at zero and at infinity. In physics an often used equation is the hypergeometric equation z

dw d2 w + [c − (1 + a + b)z] − abw = 0. (4) 2 dz dz This equation has three regular singular points, at zero, one and infinity. Jacobi, Legendre, Gegenbauer, Tchebycheff equations are special forms of this equation. When the singular points at z=1 and z equals infinity are ” coalesced” [5] at infinity, we get the confluent hypergeometric equation z(1 − z)

z

d2 w dw + (c − z) − aw = 0, dz 2 dz 2

(5)

with an essential singularity at infinity and a regular singulariy at zero. Bessel, Laguerre, 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 and a recursion relation can be found between two consecutive coefficients of the series. This fact allows us having an idea about the general properties of the solution, as 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 [6]. This is an equation with four regular singular points. This equation is discussed in the book edited by Ronveaux [7]. All the general information we give below is taken from this book. As discussed there [8], any equation with four regular singular points can be transformed to the equation given below: d2 w c d e dw abz − q +[ + + ] − w = 0. 2 dz z z − 1 z − f dz z(z − 1)(z − f )

(6)

There is a relation between the constants given as a + b + 1 = c + d + e . Then this equation has regular singularities at zero, one, f and infinity. If we try to obtain a solution in terms of a power series, one can not get a recursion relation between two consecutive coefficients. Such a relation exists between at least three coefficients. A simple solution as an integral transform also can not be found [1]. One can obtain different confluent forms of this equation. When we ”coalesce” two regular singular points at zero and infinity, we get the confluent Heun equation [9]: d dw m2 + s2 + 2msz ((z 2 − 1) ) + [−p2 (z 2 − 1) + 2pβz − λ − ]w = 0. dz dz (z 2 − 1)

(7)

Special forms of this equation are obtained in problems with two Coulombic centers [10], dw m2 d ((z 2 − 1) ) + [−p2 (z 2 − 1) + 2pβz − λ − 2 ]w = 0, dz dz (z − 1)

(8)

whose special form is the spheroidal equation [10], d dw m2 ((z 2 − 1) ) + [−p2 (z 2 − 1) − λ − 2 ]w = 0. dz dz (z − 1)

(9)

dw 1 d ((z 2 − 1) ) + [−p2 (z 2 − 1) − λ − ]w = 0. dz dz 4(z 2 − 1)

(10)

Another form is the algebraic form of the Mathieu equation [10]:

If we coalesce two regular singular points pairwise, we obtain the double confluent form [11]: D2 w + (α1 z +

α−1 α1 α1 α−1 α−1 1 )Dw + [(B1 + )z + (B0 + ) + (B−1 − ) )]w = 0. z 2 2 2 z (11) 3

d Here D = z dz . If we equate the coefficients αi , i = −1, 1 and set Bi , i = −1, 1 to zero, we can reduce the new equation to the Mathieu equation [12], an equation with two irregular singularities at zero and at infinity. 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 [13]:

z2

dw d2 w +z w + (A0 + A1 z + A2 z 2 + A3 z 3 − z 4 )w = 0. 2 dz dz

(12)

The anharmonic equation in three dimensions can be reduced to this equation [14]: d2 w ν + (E − 2 − µr2 − λr4 − ηr6 )w = 0. (13) 2 dr r In the triconfluent case, all regular singular points are ”coalesced” at infinity which gives the equation below [15]: d2 w 9 + (A0 + A1 z + A2 z 2 − z 4 )w = 0. 2 dz 4

(14)

These different forms are used in different physics problems. In SCI I found 94 papers when I searched for Heun functions. More than two thirds of these papers were published in the last ten years. The rest of the papers were published between 1990 and 1999, except a single paper in 1986 on a mathematical problem [16] . 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. The list of books on this topic also is not very long. There is a book edited by A.Ronveau, 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 is an extensive bibliography at the end of this book. I also found frequent references to two other books: Mathieusche Funktionen und Sphaeroidfunktionen mit Anwendungen auf physikalische und technische Probleme by Joseph Meixner and Friedrich Wilhelm Schaefke, published by Springer Verlag in 1954 [17] and a 1963 Dover reprint of a book first published in 1946, Theory and Applications of Mathieu Functions by N.W. McLachlan [18] . These are on functions which are special cases of the Heun equation. Some well known papers on different mathematical properties of these functions can be found in references [19] - [23]. 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 especially in ” the theoretical theory” section of particle physics anymore and people doing research on this field 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 4

the Eguchi-Hanson case as an example. The wave equation in the background of the Eguchi-Hanson metric [24] in four dimensions has hypergeometric functions as solutions [25] whereas the Nutku helicoid [26, 27] metric, another instanton metric which is little bit more complicated than the Eguchi-Hanson one, gives us Mathieu function solutions [28], a member of the Heun function set. We also find that the same equation in the background of the Eguchi-Hanson metric, trivially extended to five dimensions, gives Heun type solutions [29]. Note that the problem need not be very complicated to end up with these equations. We encounter Mathieu functions if we consider two dimensional problems with elliptical boundaries[30]. Let us use x = 12 a cosh µ cos θ, y = 1 2 a sinh µ sin θ, where a is the distance from the origin to the focal point. Then the Helmholtz equation can be written as 1 ∂µµ ψ + ∂θθ ψ + a2 k 2 [cosh2 µ − cos2 θ]ψ = 0 4

(15)

which separates into two equations d2 H + (b − h2 cos2 θ)H = 0, dθ2

(16)

d2 M + (b − h2 cosh2 µ)M = 0. dµ2

(17)



a, k, h are constants given in [30]. 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 [31]. Solution of the Schr¨odinger equation to symmetric double Morse 2 potentials, like V (x) = B4 sinh2x − (s + 12 )Bcoshx , where s = (0, 1/2, 1, ) [32], also needs these functions. Similar problems are treated in references [33, 34, 35]. o In atomic physics further problems such as separated double wells, Stark effect, hydrogen molecule ion use these functions. Physics problems which end up in these equations are given in the book by S.Y. Slavyanov and S. Lay [36]. 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 [37], also treated by S.Yu Slavyanov [38]), when all the relevant constants, namely Planck constant over 2π, electron mass and electron charge are set to unity, the Schrodinger equation for the hydrogen atom in a constant electric field of magnitude F in the z direction is given by 

1  ∆ + 2[E − (F z − )] Ψ = 0. r

Here ∆ is the laplacian operator. 5

(18)

Using parabolic coordinates, √ where the√cartesian ones are given in terms of the new coordinates by x = ξηcosφ, y = ξηsinφ, z = ξ−η 2 and writing the wave function in the product form p (19) Ψ = ξηV (ξ)U (η)exp(imφ],

we get two separated equations:

E d2 V β1 F 1 − m2  V (ξ) = 0, + + + ξ + dξ 2 2 ξ 4 4ξ 2

(20)

β2 F 1 − m2  d2 U  E U (η) = 0. + + + η + dη 2 2 η 4 4η 2

(21)

Here β1 and β2 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 reference [36], (original reference is [39]). 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 +r2 +r2 coordinates ξ = r12c , η = r12c are used. Here c is the distance between the two centers. Assuming p ψ = ξηV (ξ)U (η)exp(imφ), (22) we achieve separation into two equations:

 dV   2 2 m2 d (1 − ξ 2 ) + λ ξ − κξ − + µ V = 0, dξ dξ 1 − ξ2

(23)

 dU   2 2 m2 d + µ U = 0. (1 − η 2 ) + λ η − 2 dη dη 1−η

(24)

Both equations are of the confluent Heun type. If we mention some recent papers on atomic physics with Heun type solutions we find three relatively recent papers which treat atoms in magnetic fields: o Exact low-lying states of two interacting equally charged particles in a magnetic field are studied in by Truong and Bazzali [40]. o The energy spectrum of a charged particle on a sphere under a magnetic field and Coulomb force is studied by Ralko and Truong [41]. Both papers reduce the problem to a biconfluent Heun equation. o B.S. Kandemir presented an analytical analysis of the two-dimensional Schr¨odinger equation for two interacting electrons subjected to a homogeneous magnetic field and confined by a two-dimensional external parabolic potential. Here also a biconfluent Heun equation is used [42]. o Dislocation movement in crystalline materials, quantum diffusion of kinks along dislocations are some solid state applications of this equation. The book by S.Y. Slavyanov and S. Lay [36] is a general reference on problems solved before 2000. 6

o In a relatively recent work P. Dorey, J. Suzuki, R. Tateo [43] show that equations in finite lattice systems also reduce to Heun equations. In the rest of this paper I will comment only on papers on particle physics and general relativity. o In a relatively early work, Teukolsky studied the perturbations of the Kerr metric and found out that they were described by two coupled singly confluent Heun equations [44]. o Quasi-normal modes of rotational gravitational singularities were also studied by by E.W. Leaver [45] with Heun type equations. In recent applications in general relativity, these equations become indispensible when one studies phenomena in higher dimensions, or in different geometries. Some references are: o D. Batic, H. Schmid, M. Winklmeier studied the Dirac equation in the Kerr-Newman metric and static perturbations of the non-extremal ReisnerNordstrom solution and ended up with Heun type solutions [46]. D. Batic and H. Schmid also studied the Dirac equation for the Kerr-Newman metric and computed its propagator [47]. They found that the equation satisfied is a form of a generalized Heun equation described in Reference [46]. In later work Batic, with collaborators, continued studying Heun equations and their generalizations [48, 49]. Prof. P.P. Fiziev studied problems whose solutions are Heun equations extensively: o In a paper published in gr-qc/0603003, he studied the exact solutions of the Regge-Wheeler equation in the Schwarschild black hole interior [50]. In another paper he obtained exact solutions of the same equation and quasi-normal modes of compact objects [51]. o He presented a novel derivation of the Teukolsky-Starobinsky identities, based on properties of the confluent Heun functions [52]. These functions define analytically all exact solutions to the Teukolsky master equation, as well as to the Regge-Wheeler and Zerilli ones. o 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 [53]. o All classes of exact solutions to the Teukolsky master equation were described in terms of confluent Heun functions in Reference [54]. o In reference [55] 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 subclasses of confluent Heun’s functions are introduced and studied. A new alternative derivation of confluent Heun’s polynomials is presented.

7

o In another paper [56] he, with a collaborator, noted that weak gravitational, electromagnetic, neutrino and scalar fields, considered as perturbations on Kerr background satisfied the Teukolsky Master Equation. 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. Among other papers on this subject one may cite the following papers: o In reference [57], R. Manvelyan, H.J.W. M¨ uller Kirsten, J.Q. Liang and Y. Zhang, calculated the absorption rate of a scalar field by a D3 brane in ten dimensions in terms of modified Mathieu functions, and obtained the S-matrix . o In reference [58], 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. o In reference [59], 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 reduced to a Heun equation. o A. Al-Badawi and I. Sakallı studied the Dirac equation in the rotating Bertotti-Robinson spacetime [60], ending up with a confluent Heun equation for the angular part of the equation for a massless particle. o Mirjam Cvetic, 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 [61]. When they studied the similar problem for the rotating black hole with four U (1) charges, they again obtained a confluent Heun equation for the radial component of the Klein Gordon equation, which they reduced to the hypergeometric form by making approximations [62]. These two papers are partly repeated in [63]. Same equations were obtained which were reduced to approximate forms which gave solutions in the hypergeometric form. Other relevant references I could find are listed as references [64] - [71]. I first ”encountered” this type of equation when we tried to solve the scalar wave equation in the background of the Nutku helicoid instanton [26]. In this case one gets the Mathieu equation which is a special case of the Heun class [28]. o The helicoid instanton is a double-centered solution. As remarked above, for the simpler instanton solution of Eguchi-Hanson [24], hypergeometric solutions are sufficient [25]. Here one must remark that another paper using the Eguchi-Hanson metric ends up with the confluent Heun equation [72]. These two papers show that sometimes judicious choice of the coordinate system and separation ansatz matters. ¨ o Sucu and Unal obtained closed solutions for the spinor particle written in the background of the Nutku helicoid instanton metric [25]. One can show that 8

these solutions 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 [73]. In the subsequent sections I will summarize some work Tolga Birkandan and I have done on this topic [29, 74, 75].

2

Dirac equation in the background of the Nutku helicoid metric [74, 75]

The Nutku helicoid metric is given as ds2

=

1 q 1+

  a2 [dr2 + (r2 + a2 )dθ2 + 1 + 2 sin2 θ dy 2 r a2 r2

2



  a a2 2 sin 2θdydz+ 1 + cos θ dz 2 ]. r2 r2

(25)

where 0 < r < ∞, 0 ≤ θ ≤ 2π, y and z are along the Killing directions and will be taken to be periodic coordinates on a 2-torus [28]. If we make the following transformation r = a sinh x,

(26)

the metric is written as ds2

=

a2 sinh 2x(dx2 + dθ2 ) 2 2 + [(sinh2 x + sin2 θ)dy 2 sinh 2x − sin 2θdydz + (sinh2 x + cos2 θ)dz 2 ].

(27)

We write the system in the form Lψ = Λψ, where L is the Dirac operator, Λ is the eigenvalue, and try to obtain the solutions for the different components. We use the Newman-Penrose formalism [76, 77] to write the Dirac equation, with four components, in this metric. We use a solution of the type ψi = ei(ky y+kz z) Ψi (x, θ). We also use ky = kcosφ, kz = ksinφ. The transformation 1 f1,2 is used for the upper components to have similar equations Ψ1,2 = √sinh2x for all components. Then these equations read: a (∂x + i∂θ )Ψ3 + iak[cos(θ − φ + ix)]Ψ4 = Λ √ f1 , 2

(28)

a (∂x − i∂θ )Ψ4 − iak[cos(θ − φ − ix)]Ψ3 = Λ √ f2 , 2

(29)

9

(−∂x + i∂θ )f1 + iak[cos(θ − φ + ix)]f2 = −Λ

asinh2x √ Ψ3 , 2

(30)

(−∂x − i∂θ )f2 − iak[cos(θ − φ − ix)]f1 = −Λ

asinh2x √ Ψ4 . 2

(31)

We solve our equations in terms of f1,2 and substitute these expressions in equations, given above. This substitution gives us second order, but uncoupled equations for the lower components:   a2 2 2 ∂xx + ∂θθ + [k (−cos[2(θ − φ)] − cosh2x) − Λ sinh2x] Ψ3,4 = 0. (32) 2 We can separate this equation into two ordinary differential equations by the ansatz Ψ3,4 = R(x)S(θ − φ). Using this ansatz gives us two ordinary differential equations. The equation for S reads   2 a 2 k cos(2Θ) − n S(Θ) = 0, (33) ∂ΘΘ S(Θ) − 2 where Θ = (θ − φ). This equation is of the Mathieu type and the solution can be written immediately.

S(θ) = C1 Se(n,

a2 k 2 a2 k 2 , θ − φ) + C2 So(n, , θ − φ), 4 4

(34)

where Se, So are the even and odd Mathieu functions. The solutions should be periodic in the angular variable Θ. This fact forces n, the separation constant, to take discrete values. It is known that the angular Mathieu functions satisfy an orthogonality relation such that functions with different n values are perpendicular to each other. Here, we integrate the angular variable from zero to 2π. The equation for R(x) reads   a2 2 2 ∂xx − [ (k cosh2x) + Λ sinh2x] + n R(x) = 0. (35) 2 This solution is of the double confluent form which can be reduced to the form R(x) = D1 Se(n, A6 , i(x + b)) + D2 So(n, A6 , i(x + b)).

(36)

by several transformations. Here C1 , C2 , D1 , D2 are arbitrary constants. A6 , b are given in the original reference [74]. As a result of this analysis we see that the solutions of the Dirac equation, written in the background of the Nutku helicoid metric, can be expressed in 10

terms of Mathieu functions, which is a special form of Heun function [74, 75]. Mathieu function is a related but much more studied function with singularity structure sme as the double confluent Heun equation. o One can also show [74, 75] that one can use the similar metric in five dimensions and obtain, in general, double confluent Heun functions which can be reduced to Mathieu functions.

3

Scalar field in the background of the extended Eguchi-Hanson solution [29]

To go to five dimensions, we can add a time component to the Eguchi-Hanson metric [24] so that we have ds2 = −dt2 +

1 a4 2 2 2 2 2 2 dr + r (σ + σ ) + r (1 − )σ 4 x y r4 z 1 − ar4

(37)

where

1 (− cos ξdθ − sin θ sin ξdφ) 2 1 σy = (sin ξdθ − sin θ cos ξdφ) 2 1 σz = (−dξ − cos θdφ). 2 This is a vacuum solution. If we take the solution to the wave equation Φ as σx =

1

Φ = eikt einφ ei(m+ 2 )ξ ϕ(r, θ),

(38) (39) (40)

(41)

we can write the reduced wave equation as

∆ϕ(r, θ)

r 4 − a4 3r4 + a4 4a4 m2 2 2 ∂ + ∂ + k r + + rr r r2 r3 a4 − r 4 8mn cos θ − 4(m2 + n2 ) 4∂θθ + 4 cot θ∂θ + )ϕ(r, θ) = 0. (42) sin2 θ

=

(

If we take ϕ(r, θ) = f (r)g(θ), the solution of the radial part is expressed in terms of confluent Heun (HC ) functions.   1 2 2 a2 + r 2 1 2 2 1 2 1 f (r) = −a + r HC 0, m, m, k a , m − λ − k a , 2 2 4 4 2a2   − 1 m 2 1 m 1 1 1 a2 + r 2 1 . + a2 + r 2 2 r − a2 2 HC 0, −m, m, k 2 a2 , m2 − λ − k 2 a2 , 2 2 4 4 2a2 (43) 4

4

 21 m

11

√ If the variable transformation r = a cosh x is made, the solution can be expressed as

f (x) =

  1 1 1 1 1 1 1 {(sinh (x))m+1 HC 0, m, m, k 2 a2 , m2 − λ − k 2 a2 , cosh (x) + sinh x 2 2 4 4 2 2 1

1

1

1

+ (2 cosh (x) + 2)− 2 m+ 2 (2 cosh (x) − 2) 2 m+ 2

(44)

  1 1 1 1 1 1 ×HC 0, −m, m, k 2 a2 , m2 − λ − k 2 a2 , cosh (x) + }. 2 2 4 4 2 2 2

2

+r to We tried to express the equation for the radial part in terms of u = a 2a 2 see the singularity structure more clearly. Then the radial differential operator reads     d2 m2 1 1 d 1 1 4 2 +4 + 2 . (45) + + k 2 a2 + du u − 1 u du u−1 u u (1 − u)2

This 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 angular solution is in terms of hypergeometric functions.

g (θ) =

1 p { 2 − 2 cos (θ) sin θ



1 1 cos (θ) − 2 2 1

 12 m 

1 1 cos (θ) − 2 2

− 21 n

.

− 1 n− 1 m

[(2 cos (θ) + 2) 2 2 2 1 1 1√ 1√ 1 1 λ + 1 + , −n − λ + 1 + ], [1 − n − m], cos (θ) + ) × 2 F1 ([−n + 2 2 2 2 2 2 (46) 1 + 1 n+ 1 m + (2 cos (θ) + 2) 2 2 2 1√ 1 1 1√ 1 1 ×2 F1 ([m + λ + 1+ ,m− λ + 1 + ], [1 + n + m], cos (θ) + )]}. 2 2 2 2 2 2

4

Derivation of a new generalized Heun equation [29]

The Dirac equation, with Λ = 0, written in the background of the Nutku helicoids metric, reads as follows: (∂x + i∂θ )Ψ3 + iak[cos(θ − φ + ix)]Ψ4 = 0, 12

(47)

(∂x − i∂θ )Ψ4 − iak[cos(θ − φ − ix)]Ψ3 = 0,

(48)

(−∂x − i∂θ )f2 − iak[cos(θ − φ − ix)]f1 = 0.

(50)

Θ = θ − φ − ix0

(55)

(−∂x + i∂θ )f1 + iak[cos(θ − φ + ix)]f2 = 0,

(49)

These equations have simple solutions [25] which can also be expanded in terms of products of radial and angular Mathieu functions [73], [74]. Problem arises when these solutions are restricted to boundary [29, 75]. To impose these boundary conditions we need to write the little Dirac equation, the Dirac equation restricted to the boundary, where the variable x takes a fixed value x0 . We choose to write the equations in the form, √ 2 d {i Ψ3 + ikacos(θ − φ + ix0 )Ψ4 } = λf1 , (51) a dθ √ d 2 {−i Ψ4 − iakcos(θ − φ − ix0 )Ψ3 } = λf2 , (52) a dθ √ d 2 {−i f1 − iakcos(θ − φ + ix0 )f2 } = λΨ3 , (53) a dθ √ 2 d {i f2 + iakcos(θ − φ − ix0 )f1 } = λΨ4 . (54) a dθ Here λ is the eigenvalue of the little Dirac equation. We take λ = 0 as the simplest case. The transformation

can be used. Then we solve f1 in the latter two equations in terms of f2 : d (ak)2 d2 f −tan Θ f + [cos(2Θ) cosh(2x0 )−i sin(2Θ) sinh(2x0 )+cosh(2x0 )]f2 = 0. 2 2 dΘ2 dΘ 2 (56) When we make the transformation



u = e2iΘ ,

(57)

the equation reads, d 1 d (ak)2 d2 + ]−2iu(u−1) + (u+1)[ue−2x0 + e2x0 +cosh(2x0 )]}f2 = 0. 2 du du du 2 u (58) This equation has irregular singularities at u = 0 and ∞ and a regular singularity at u = −1. This is still another generalized Heun equation, different from the one given by by reference [46]. {4(u+1)u[u

Acknowledgement: I thank the organizers of 13th Regional Conference on Mathematical Physics for inviting me to this meeting and dedicate this work to the memory of my dear friend Faheem Hussain. I thank Tolga Birkandan ¨ ITAK, ˙ for collaboration. This work is supported by TUB the Scientific and ¨ Technological Council of Turkey and by TUBA, the Academy of Sciences of Turkey. 13

References [1] Philip M. Morse and Herman Feshbach, Methods of Theoretical Physics, McGraw Hill Company, New York (1953) [2] E.L. Ince, Ordinary Differential Equations, Dover Publications (1926, 1956) [3] Philip M. Morse and Herman Feshbach, Methods of Theoretical Physics, McGraw Hill Company, New York (1953), p.532 [4] Philip M. Morse and Herman Feshbach, Methods of Theoretical Physics, McGraw Hill Company, New York (1953), pp. 667-673 [5] P. Dennery, A. Krzywicki, Mathematics for Physicists , Dover Publications (1967) [6] K. Heun, Mathematische Annalen 33, 161 (1889) [7] A. Ronveaux (ed.), Heun’s Differential Equations, Oxford University Press (1995) [8] F. M. Arscott in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995),p.7 [9] S.Yu. Slavyanov in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p.92 [10] S.Yu. Slavyanov in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p.94 [11] D.Schmidt and G. Wolf in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p.135 [12] D.Schmidt and G. Wolf in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p. 140-141 [13] P. Maroni in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p.193 [14] P. Maroni in the book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p.199 [15] A.Duval in tyhe book edited by A. Ronveaux, Heun’s Differential Equations, Oxford University Press (1995), p. 253 [16] G. Valent, Siam Journal on Mathematical Analysis 17, 688-703 (1986) [17] J. Meixner and F. W. Schaefke, Mathieusche Funktionen und Sphaeroidfunktionen mit Anwendungen auf physikalische und technische Probleme, Springer Verlag (1954)

14

[18] N.W. McLachlan, Theory and Applications of Mathieu Functions, Dover reprint from 1946 edition (1963). [19] R. Schafke, D. Schmidt, SIAM Journal of Mathematical Analysis 11, 848 (1980) [20] R.S. Maier, math.CA/0408317 [21] R.S. Maier, math.CA/0203264, also Journal of Differential Equations 213, 171-203 (2005) [22] N. Gurappa, P.K. Panigrahi, Journal of Physics A: Mathematical and General 37, L605-L608 (2004) [23] K. Kuiken, SIAM Journal of Mathematical Analysis 10, 655-567 (1979) [24] T. Eguchi, A.J. Hanson, Physics Letters 74B, 249 (1978) ¨ [25] Y. Sucu, N. Unal, Classical and Quantum Gravity 21, 1443 (2004) [26] Y. Nutku, Physical Review Letters 77, 4702 (1996) [27] D. Lorenz-Petzold, Journal of Mathematical Physics 24, 2632 (1983) [28] A.N. Aliev, M. Horta¸csu, J. Kalaycı, Y. Nutku, Classical and Quantum Gravity 16, 631 (1999) [29] T. Birkandan, M. Horta¸csu, Journal of Mathematical Physics 49, 0054101 (2008) [30] Philip M. Morse and Herman Feshbach, Methods of Theoretical Physics McGraw Hill Company, New York (1953), p.1407 [31] B.D.B. Figueiredo, Journal of Mathematical Physics 46, 113503 (2005), math-phys/0509013 [32] B.D.B. Figueiredo, Journal of Mathematical Physics 48, 013503 (2007) [33] B.D.B. Figueiredo, math-ph 10402071 [34] L.J. El-Jaick, B.D.B. Figueiredo, Journal of Mathematical Physics 49, 083508 (2008) [35] L.J. El-Jaick, B.D.B. Figueiredo, Journal of Mathematical Physics 50,123511 (2008) [36] S.Y. Slavyanov, W. Lay, Special Functions, A Unified Theory Based on Singularities, Oxford University Press (2000) [37] P.S. Epstein, Physical Review 2, 695 (1926)

15

[38] S.Y. Slavyanov, Asymptotic Solutions of the One-dimensional Schrodinger Equation (Leningrad University Press) (in Russian) (1991), Translation into English: S.Y. Slavyanov, Asymptotic Solutions of the One-dimensional Schrodinger Equation (Amer. Math. Soc. Trans. of Math. Monographs) 151 (1996) [39] A.H. Wilson, Proceedings of Royal Society A118, 617-647 (1928), also reference [36], p.167 [40] T.T. Truong, D. Bazzali, Physics Letters A 269, 186-193 (2000) [41] A. Ralko, T.T. Truong, Journal of Physics A-Mathematical and General 35, 9573-9583 (2002) [42] B.S. Kandemir, Journal of Mathematical Physics 46, 032110 (2005) [43] P.Dorey, J. Suzuki and R. Tateo, Journal of Physics A: Mathematical and General 37, 2047-2062 (2004) [44] S.A. Teukolsky, Physical Review Letters 29, 1114 (1972), also reference [36], p.170 [45] E.W. Leaver, Proceedings of Royal Society London A 402, 285 (1985), also reference [36], p.171 [46] D. Batic, H. Schmid, M. Winklmeier, Journal of Physics A: Mathematical and General 39, 12559 (2006), gr-qc/0607017 [47] D. Batic, H. Schmid, Progress in Theoretical Physics 116, 517-544 (2006) [48] D. Batic, H. Schmid, Journal of Mathematical Physics 48, 042502(2007) [49] D. Batic, M. Sandoval, e-Print: arXiv:0805.4399 [gr-qc] [50] P.P. Fiziev, e-print gr-qc/0603003 [51] P.P. Fiziev, Classical and Quantum Gravity 23, 2447 (2006), gr-qc/0509123 [52] P.P. Fiziev, Physical Review D 80, 124001 (2009) [53] P.P. Fiziev, Journal of Physics Conference Series 66, 012016, (2007), gr-qc/0702014 [54] P.P. Fiziev, Classical and Quantum Gravity 27, 135001 2010 [55] P.P. Fiziev, Journal of Physics-Mathematical and Theoretical 43, 035203 (2010) [56] Roumen S. Borissov, Plamen P. Fiziev, e-print: arXiv:0903.3617 [gr-qc] [57] R. Manvelyan, H.J.W. M¨ uller-Kirsten, J.Q. Liang and Y. Zhang, Nuclear Physics B579, 177-208 (2000), hep-th/0001179 16

[58] T. Oota, Y. Yasui, Nuclear Physics B 742, 275-294 (2006) [59] S. Musiri, G. Siopsis, hep-th/0308196,

Physics Letters B 576,

309-313 (2003),

[60] A. Al-Badawi, I. Sakallı, Journal of Mathematical Physics 49, 052501 (2008) [61] M. Cvetic and F. Larsen, Physical Review D 56, 4994-5007 (1997), hep-th/9705192 [62] M. Cvetic and F. Larsen, Nuclear Physics B 506, 107-120 (1997), hep-th/9706071 [63] M.Cvetic and F. Larsen, Journal of High Energy Physics, 0909:088 (2009), arXiv:0908.1136 [hep-th] [64] G. Siopsis, Nuclear Physics B 715, 483-498 (2005), hep-th/0407157 [65] L. Anguelva, P. Langfelder, Journal of High Energy Physics, Issue: 3 Article 057, March 2003 [66] S.R. Lau, Classical and Quantum Gravity 21, 4147-4192 (2004), also Journal of Computational Physics 199, 376-422 (2004) [67] H. Suzuki, E. Takasugi and H. Umetsu, Progress of Theoretical Physics 100, 491-505 (1998) [68] A. Zecca, Nuovo Cimento B 125, 191-199 (2010) [69] A. Ensico and N. Kamran, Communications in Mathematical Physics 290, 105-127 (2009) [70] G. Esposito and R. Roychowdhury, General Relativity and Gravitation 42, 1221-1238 (2010) [71] S. Yoshida, N. Uchikata and T. Fumatase, Physical Review D 81, 044005 (2010) [72] A. Malmendier, Journal of Mathematical Physics 44, 4308-4343 (2003) [73] L. Chaos-Cador, E. Ley-Koo, Revista Mexicana de Fisica 48, 67 (2002) [74] T. Birkandan, M. Horta¸csu, Journal of Physics A: Mathematical and Theoretical 40, 1105-1116 (2007) [75] T. Birkandan, M. Horta¸csu, Journal of Mathematical Physics 48, 092301 (2007) [76] E.T. Newman and R. Penrose, Journal of Mathematical Physics 3, 566 (1962) [77] E.T. Newman and R. Penrose, Journal of Mathematical Physics 4, 998 ((1962) 17