# Expansions of the Solutions of the General Heun Equation Governed by Two-Term Recurrence Relations for Coefficients.

1. IntroductionThe general Heun equation [1-3], which is the most general second-order linear ordinary differential equation having four regular singular points, is currently widely encountered in physics and mathematics research (see, e.g., [1-14] and references therein). However, this equation is much less studied than its immediate predecessor, the Gauss hypergeometric equation, which is the most general equation having three regular singular points. A reason for the slow progress in the development of the theory is that the solutions of the Heun equation (as well as its four confluent reductions) in general are not expressed in terms of definite or contour integrals involving simpler functions [2,3]. Furthermore, the convergence regions of power-series expansions near different singularities are rather restricted and several complications arise in studying the relevant connection problems [2,3, 15]. Another general problematic point is that the power-series solutions of the Heun equation are governed by three-term recurrence relations between successive coefficients of expansions [1-3], instead of two-term ones appearing in the hypergeometric case [16-18]. As a result, in general the coefficients are not determined explicitly.

In the present paper, we show that there exist some particular choices of the involved parameters for which the recurrence relations governing the power-series expansions become two-term. In these cases the solution of the general Heun equation can be written either as a linear combination of a finite number of the Gauss hypergeometric functions or in terms of a single generalized hypergeometric function. This is a main result of the present paper.

Another major result we report here is that in the case of the expansions of the solutions of the Heun equation in terms of hypergeometric functions there also exist particular choices of the involved parameters for which the governing three-term recurrence relations for expansion coefficients become two-term. In these cases the coefficients are explicitly written in terms of the gamma functions.

Expansions of the solutions of the Heun equation in terms of the Gauss hypergeometric functions [sub.2][F.sub.1], initiated by Svartholm [19], suggest a notable extension of the series technique. This is a useful approach applicable to many differential equations including those of more general type whose nature, outside a certain region of the extended complex plane containing only two regular singular points, is not necessary to be specified exactly. Expansions involving functions other than powers have been applied to the general and confluent Heun equations by many authors. The ordinary hypergeometric [19-25], confluent hypergeometric [26-29], Coulomb wave functions [30,31], Bessel and Hankel functions [32], incomplete Beta and Gamma functions [33-35], Hermite functions [36,37], Goursat and Appell generalized hypergeometric functions of two variables of the first kind [38,39], and other known special functions have been used as expansion functions. A useful property suggested by these expansions is the possibility of deriving finite-sum solutions by means of termination of the series.

As far as the expansions of the general Heun equation in terms of the hypergeometric functions are concerned, in the early papers by Svartholm [19], Erdelyi [20,21], and Schmidt [22], the intuitive intention was to apply hypergeometric functions with parameters so chosen as to match the Heun equation as closely as possible. For this reason they used functions of the form [sub.2][F.sub.1] ([lambda] + n, [mu] - n; [gamma]; z), which have matching behavior in two singular points, z = 0 and 1. These functions have the following Riemann P-symbol representation:

[mathematical expression not reproducible]. (1)

Here [lambda] and [mu] may adopt several values provided 1 + [lambda] + [mu] = [gamma] + [delta] (see [25]). It is clear that the functions have matching characteristic exponents at z = 0 and z = 1, and their behavior does not match that of the Heun function at the third singular point of the hypergeometric equation z = [infinity].

However, it has been shown that one can also use functions that have matching behavior at only one singular point [38]. Exploring this idea, in the present paper we discuss the hypergeometric expansions of the solutions of the Heun equation in terms of functions of the form [sub.2][F.sub.1] ([alpha], [beta]; [[gamma].sub.0] [+ or -] n; z) which have matching behavior (i.e., characteristic exponents) only at the singular point z = [infinity]. These functions are presented by the Riemann P-symbol

[mathematical expression not reproducible], (2)

where 1 + [alpha] + [beta] = [[gamma].sub.0] + [[delta].sub.0], and the parameters [[gamma].sub.0], [[delta].sub.0] are chosen so that the Fuchsian condition for the general Heun equation (3), i.e., 1 + [alpha] + [beta] = [gamma] + [delta] + [epsilon], is fulfilled: [[gamma].sub.0] + [[delta].sub.0] = [gamma] + [delta] + [epsilon]. Note that these functions differ also from the Jacobi-polynomials used by Kalnins and Miller whose functions can be written in terms of hypergeometric functions of the form [sub.2][F.sub.1]([lambda] + n, [mu] - n; v + 2n; z) [23].

In the present paper we discuss several expansions in terms of the mentioned hypergeometric functions. In general, the coefficients of the expansions obey three-term recurrence relations similar to those known from previous developments [19-25]. However, tor certain choices ot the involved parameters the recurrence relations reduce to two-term ones. In these exceptional cases the coefficients of the expansions are explicitly calculated. The result is expressed in terms of the gamma functions.

Discussing the conditions for deriving finite-sum solutions by means of termination of the presented series, we show that the termination is possible if a singularity of the Heun equation is an apparent one. Furthermore, we show that any finite-sum solution of the general Heun equation derived in this way has a representation through a single generalized hypergeometric function [sub.p][F.sub.q]. The general conclusion is then that in any such case the power-series expansion of the Heun function is governed by a two-term recurrence relation (obviously, this is the relation obeyed by the corresponding power-series for [sub.p][F.sub.q]).

2. Hypergeometric Expansions

The general Heun equation written in its canonical form is [1]

[d.sup.2]u/d[z.sup.2] + ([gamma]/z + [delta]/z - 1 + [epsilon]/z - a) du/dz + [alpha]pbeta]z - q/z(z - 1)(z - a)u = 0, (3)

where the parameters satisfy the Fuchsian relation 1 + [alpha] + [beta] = [gamma] + [delta] + [epsilon]. We introduce an expansion of this equation's solution of the form

u = [[infinity].summation over (n=0)] [c.sub.n] x [sub.2][F.sub.1]([alpha], [beta]; [[gamma].sub.0] + n; z) (4)

with the involved Gauss hypergeometric functions obeying the equation

[d.sup.2][u.sub.n]/d[z.sup.2] + ([[gamma].sub.0] + n/z + [delta] + [epsilon] + [gamma] - [[gamma].sub.0] - n/z - 1) d[u.sub.n]/dz + [alpha][beta]/z(z - 1)[u.sub.n] = 0. (5)

Substitution of (4)-(5) into (3) gives

[summation over (n)][c.sub.n][([gamma] - [[gamma].sub.0] - n/z - [epsilon] + [gamma] - [[gamma].sub.0] - n/z - 1 + [epsilon]/z - a)d[u.sub.n]/dz + [alpha][beta]a - q/z(z - 1)(z - a)[u.sub.n]] = 0 (6)

or [summation over (n)] [c.sub.n]{[(a - a)([epsilon] + [gamma] - [[gamma].sub.0] - n)z - a([gamma] - [[gamma].sub.0] - n)(z - 1)]d[u.sub.n]/dz + ([alpha][beta]a - q) x [u.sub.n]} = 0. (7)

Now, using the following relations between the involved hypergeometric functions [16-18]:

z [d[u.sub.n]/dz] = [[gamma].sub.n-1][[u.sub.n-1] - [u.sub.n]], (8)

(z - 1) d[u.sub.n]/dz = -[[delta].sub.n+1][u.sub.n] + ([[delta].sub.n+1] - [alpha][beta]/[[gamma].sub.n])[u.sub.n+1], (9)

where [[gamma].sub.n] = [[gamma].sub.0] + n and [[delta].sub.n] = [delta] + [epsilon] + [gamma] - [[gamma].sub.n], this equation is rewritten as

[mathematical expression not reproducible], (10)

from which we get a three-term recurrence relation for the coefficients of the expansion (4)

[R.sub.n][c.sub.n] + [Q.sub.n-1][c.sub.n-1] + [P.sub.n-2][c.sub.n-2] = 0 (11)

with

[R.sub.n] = (a - 1)([epsilon] + [gamma] - [[gamma].sub.n])([[gamma].sub.n] - 1), (12)

[Q.sub.n] = -[R.sub.n] + a([gamma] - [[gamma].sub.n])([alpha] + [beta] - [[gamma].sub.n]) + ([alpha][beta]a - q), (13)

[P.sub.n] = a/[[gamma].sub.n]([gamma] - [[gamma].sub.n])([alpha] - [[gamma].sub.n])([beta] - [[gamma].sub.n]). (14)

From the initial conditions [c.sub.0] = 1 and [c.sub.-1] = [c.sub.-2] = 0 we get ([epsilon] + [gamma] - [[gamma].sub.0])([[gamma].sub.0] - 1) = 0. Since [[gamma].sub.0] = 1 is forbidden (it causes division by zero at n = 1 in [P.sub.-1]) we obtain that the only possibility is [[gamma].sub.0] = [epsilon] + [gamma]. Hence, the expansion finally reads

u = [[infinity].summation over (n=0)][c.sub.n] x [sub.a][F.sub.1]([alpha], [beta]; [gamma] + [epsilon] + n; z), (15)

and the coefficients of the three-term recurrence relation (11) are

[R.sub.n] = (1 - a) n ([epsilon] + [gamma] + n - 1), (16)

[Q.sub.n] = -[R.sub.n] + a(1 + n - [delta])(n + [epsilon]) + (a[alpha][beta] - q), (17)

[P.sub.n] = -a/n + [epsilon] + [gamma](n + [epsilon])(n + [epsilon] + [gamma] - [alpha])(n + [epsilon] + [gamma] - [beta]). (18)

The expansion applies if [alpha], [beta], and [gamma] + [epsilon] are not zero or negative integers. The restrictions on [alpha] and [beta] assure that the hypergeometric functions are not polynomials of fixed degree.

The derived expansion terminates if two successive coefficients vanish. If [c.sub.N] is the last nonzero coefficient and [c.sub.N+1] = [c.sub.N+2] = 0 for some N = 0,1,2, ..., we obtain from (11) that it should be [P.sub.N] = 0 so that the termination is possible if

[epsilon] = -N or [epsilon] + [gamma] - [alpha] = -N or [epsilon] + [gamma] - [beta] = -N. (19)

Note that the equation [c.sub.N+1] = 0 results in a polynomial equation of degree N + 1 for the accessory parameter q. This equation is convenient for rewriting the recurrence relation (11) in the following matrix form:

[mathematical expression not reproducible]. (20)

The vanishing of the determinant of the above matrix gives the polynomial equation for q defining in general N + 1 values for which the termination occurs.

One may consider a mirror expansion

u = [[infinity].summation over (n=0)] [c.sub.n] x [sub.2][F.sub.1]([alpha], [beta]; [[gamma].sub.0] - n; z), (21)

which differs from expansion (4) only by the sign of n in the lower parameter of the involved hypergeometric functions. This change of the sign leads to a three-term recurrence relation (11) with the coefficients

[R.sub.n] = a/[[gamma].sub.0] - n([gamma] - [[gamma].sub.0] + n)([alpha] - [[gamma].sub.0] + n)([beta] - [[gamma].sub.0] + n), (22)

[Q.sub.n] = -[P.sub.n] + a([gamma] - [[gamma].sub.0] + n)([alpha] + [beta] - [[gamma].sub.0] + n) + [alpha][beta]a - q, (23)

[P.sub.n] = (a - 1)([epsilon] + [gamma] - [[gamma].sub.0] + n)([[gamma].sub.0] - n - 1), (24)

where

[[gamma].sub.0] = [gamma] or [alpha] or [beta]. (25)

This expansion applies if [alpha] and [beta] are not zero or negative integers and [gamma] is not an integer.

In order for the series to terminate at some n = N we put [P.sub.N] = 0 so that this time we derive

[epsilon], [epsilon] + [gamma] - [alpha] or [epsilon] + [gamma] - [beta] = -N (26)

for the expansions with [[gamma].sub.0] = [gamma] or [alpha] or [beta], respectively. Then the equation [c.sub.N+1] = 0 again gives a (N + 1)-degree polynomial equation for those values of the accessory parameter q for which the termination occurs.

3. Finite-Sum Hypergeometric Solutions

It is readily shown that the finite-sum solutions derived from the above two types of expansions by the described termination procedure coincide, as can be expected because of apparent symmetry. For example, consider the expansion (15) in the case [epsilon] = -N. The involved hypergeometric functions have the form [u.sub.n] = [sub.2][F.sub.1]([alpha], [beta], [gamma] - N + n; z). Since n = 0 ,1,..., N, we see that the set of the involved hypergeometric functions is exactly the same as in the case of the second type expansion (21) with [[gamma].sub.0] = [gamma] : {[sub.2][F.sub.1]([alpha], [beta], [gamma]; z), [sub.2][F.sub.1]([alpha], [beta], [gamma] - 1; z),..., [sub.2][F.sub.1]([alpha], [beta], [gamma] - N; z)}. Furthermore, examination of (20) shows that the equation for the accessory parameter q and the expansion coefficients are also the same for the two expansions. The same happens to other two cases: [epsilon] + [gamma] - [alpha] = -N and [epsilon] + [gamma] - [beta] = -N. Thus, while different in general, the expansions (15)-(18) and (21)-(25) lead to the same finite-sum closed-form solutions.

Consider the explicit forms of these solutions examining, for definiteness, the expansion (15)-(18). An immediate observation is that because of the symmetry of the Heun equation with respect to the interchange [alpha] [left and right arrow] [beta] the finite-sum solutions produced by the choices [epsilon] + [gamma] - [alpha] = -N and [epsilon] + [gamma] - [beta] = -N are of the same form. Furthermore, by applying the formula [16]

[sub.2][F.sub.1]([alpha], [beta]; [alpha] + k; z) = [(1 - z).sup.k-[beta]] [sub.2][F.sub.1](k, [alpha] - [beta] + k; [alpha] + k; z) (27)

to the involved hypergeometric functions [sub.2][F.sub.1]([alpha], [beta]; [alpha] - N + n; z) or [sub.2][F.sub.1]([alpha], [beta]; [beta] - N + n; z) we see that the sum is a quasipolynomial, namely, a product of [(1 - z).sup.1-[delta]] and a polynomial in z. Here are the first two of the solutions:

[epsilon] + [gamma] - [alpha] = 0, q = a[gamma]([delta] - 1), u = [(1 - z).sup.1-[delta]], [epsilon] + [gamma] - [alpha] = -1, (28)

[q.sup.2] + [[alpha] - 1 - a([delta] - 2 + [gamma](2[delta] - 3))]q - a[gamma]([delta] - 2) x ([alpha] - a(1 + [gamma])([delta] - 1)) = 0, (29)

u = [(1 - z).sup.1-[delta]] x (1 - [alpha] + 1 - [delta]/[alpha] - 1 z + q - a([alpha][beta] + [epsilon] - [delta][epsilon])/(1 - a) ([alpha] - 1)(1 - z)). (30)

Note that, since 1 - [delta] is a characteristic exponent of the Heun equation, the transformation u = [(1 - z).sup.1-[delta]] w(z) results in another Heun equation for w(z). Hence, the derived finite-sum solutions corresponding to the choices [epsilon] + [gamma] - [alpha] = -N and [epsilon] + [gamma] - [beta] = -N are generated from the polynomial solutions of the equation for w(z).

More interesting is the case [epsilon] = -N, when the finite-sum solutions involve N + 1 hypergeometric functions irreducible, in general, to simpler functions. The case [epsilon] = 0 produces the trivial result q = a[alpha][beta], when the Heun equation is degenerated into the hypergeometric equation with the solution u = [sub.2][F.sub.1]([alpha], [beta]; [gamma]; z). The solution for the first nontrivial case [epsilon] = -1 reads

u = [sub.2][F.sub.1]([alpha], [beta]; [gamma] - 1; z) + q - a[alpha] + a(1 - [delta])/(1 - a)([gamma] - 1) x [sub.2][F.sub.1]([alpha], [beta]; [gamma]; z), (31)

where q is a root of the equation

(q - a[alpha][beta] + a(1 - [delta]))(q - a[alpha][beta] + (a - 1)(1 - [gamma])) - a(1 - a)(1 + [alpha] - [gamma])(1 + [beta] - [gamma]) = 0. (32)

Note that the second term in (31) vanishes if q - a[alpha][beta] + a(1 - [delta]) = 0 so that in this degenerate case the solution involves one, not 2 = N + 1, terms. It is seen from (32) that this situation is necessarily the case if a = 1/2, [gamma] + [delta] = 2 and [alpha] or [beta] equals [gamma] - 1. We will see that the solution in this case is a member of a family of specific solutions for which the expansion is governed by two-term recurrence relations for the coefficients.

The solutions (31) and (32) have been noticed on several occasions [40-44]. It has been shown that, for [epsilon] = -1, when the characteristic exponents of z = a are 0,2 so that they differ by an integer, (32) provides the condition for the singularity z = a to be apparent (or "simple"); that is, no logarithmic terms are involved in the local Frobenius series expansion [40-42]. In fact, the Frobenius solution in this case degenerates to a Taylor series. It has further been observed that the solution (31) can be expressed in terms of the Clausen generalized hypergeometric function [sub.3][F.sub.2] with an upper parameter exceeding a lower one by unity [40-42]:

u/u(0) = [sub.3][F.sub.2]([alpha], [beta], e + 1; [gamma], e; z), (33)

where the parameter e is given as

e = a[alpha][beta]/q - a[alpha][beta]. (34)

Note that using this parameter e the solution of (32) is parameterized as [41]

q = a[alpha][beta] 1 + e/e, a = e(e - [gamma] + 1)/(e - [alpha](e - [beta]). (35)

We will now show that a similar generalized hypergeometric representation holds also for [epsilon] = -2 and for all [epsilon] [member of] Z, [epsilon] [not equal to] 1.

4. The Case [epsilon] [less than or equal to] -2, [epsilon] [member of] Z

For [epsilon] = -2 the termination equation [c.sub.N+1] = 0 for the accessory parameter q is written as

[mathematical expression not reproducible]. (36)

The solution of the Heun equation for a root of this equation is given as

u = [sub.2][F.sub.1]([alpha], [beta]; [gamma] - 2; z) + [B.sub.1] x [sub.2][F.sub.1]([alpha], [beta]; [gamma] - 1; z) + [B.sub.2] x [sub.2][F.sub.1]([alpha], [beta]; [gamma]; z), (37)

with

[mathematical expression not reproducible]. (38)

It is now checked that this solution is presented by the hypergeometric function [sub.4][F.sub.3] as [42]

u/u(0) = [sub.4][F.sub.3]([alpha], [beta], e + 1, r + 1; [gamma], e, r; z), (39)

where u(0) = 1 + [B.sub.1] + [B.sub.2] and the parameters e, r solve the equations (compare with (35))

q = a[alpha][beta][(1 + e)(1 + r)/er], (40)

[mathematical expression not reproducible]. (41)

It is further checked that this system of equations admits a unique solution e, r (up to the transposition e [left and right arrow] r).

The presented result is derived in a simple way by substituting the ansatz (39) into the general Heun equation and expanding the result in powers of z. The equations resulting in cancelling the first three terms proportional to [z.sup.0], [z.sup.1], and [z.sup.2] are that given by (40), (41), and (36), respectively. It is then shown that these three equations are enough for the Heun equation to be satisfied identically.

A further remark is that (36) presents the condition for the singularity z = a to be apparent for [epsilon] = -2. This is straightforwardly verified by checking the power-series solution u = [[summation].sup.[infinity].sub.n=0] [c.sub.n][(z - a).sup.n] with [c.sub.0] [not equal to] 0 for the neighborhood of the point z = a. In calculating [c.sub.3] a division by zero will occur, unless q satisfies (36), in which case the equation for [c.sub.3] will be identically satisfied.

It can be checked that generalized hypergeometric representations are achieved also for [epsilon] = -3,-4,-5 [42]. The conjecture is that for any negative integer [epsilon] = -N, N = 1,2,3,... there exists a generalized hypergeometric solution of the Heun equation given by the ansatz

u = [sub.N+2][F.sub.1+N]([alpha], [beta], [e.sub.1] + 1,...,[e.sub.N] + 1; [gamma], [e.sub.1],...,[e.sub.N]; z) (42)

provided the singularity at z = a is an apparent one. For the latter condition to be the case, the accessory parameter q should satisfy a (N + 1)-degree polynomial equation which forces the above expansions (15)-(18) and (21)-(25) to terminate at Nth term. Note that (42) applies also for N = 0, that is, for [epsilon] = 0, for which the Heun function degenerates to the Gauss hypergeometric function u = [sub.2][F.sub.1]([alpha], [beta]; [gamma]; z) provided q = 0.

Finally, we note that by the elementary power change u = [(z - a).sup.1-[epsilon]] a Heun equation with a positive exponent parameter e >2 is transformed into the one with a negative parameter 2 - [epsilon] < 0. Hence, it is understood that a similar generalized hypergeometric representation of the solution of the Heun equation can also be constructed for positive integer [epsilon] = N, N = 2,3, Thus, the only exception is the case [epsilon] = 1.

It is a basic knowledge that the generalized hypergeometric function [sub.p][F.sub.q] is given by a power-series with coefficients obeying a two-term recurrence relation. Since any finite-sum solution of the general Heun equation derived via termination of a hypergeometric series expansion has a representation through a single generalized hypergeometric function [sub.p][F.sub.q], the general conclusion is that in each such case the power-series expansion of the Heun function is governed by a two-term recurrence relation (obviously, by the relation obeyed by the corresponding power-series for [sub.p][F.sub.q]).

5. Hypergeometric Expansions with Two-Term Recurrence Relations for the Coefficients

In this section we explore if the three-term recurrence relations governing the above-presented hypergeometric expansions can be reduced to two-term ones. We will see that the answer is positive. Two-term reductions are achieved for an infinite set of particular choices of the involved parameters.

First, we mention a straightforward case which actually turns to be rather simple because in this case the Heun equation is transformed into the Gauss hypergeometric equation by a variable change. This is the case if

a = 1/2, [gamma] + [delta] = 2 and q = a[alpha][beta] + a(1 - [delta])[epsilon], (43)

when the coefficient [Q.sub.n] in (11) identically vanishes so that the recurrence relation between the expansion coefficients straightforwardly becomes two-term for both expansions (15)-(18) and (21)-(25). The coefficients of the expansions are then explicitly calculated. For instance, expansion (15) is written as

[mathematical expression not reproducible], (44)

where [(...).sub.k] denotes the Pochhammer symbol. The values m(0), u'(0) and u(1), and u'(1) can then be written in terms of generalized hypergeometric series [17,18]. For instance,

u(0) = [sub.3][F.sub.2]([gamma] + [epsilon] - [alpha]/2, [gamma] + [epsilon] - [beta]/2, [epsilon]/2; [gamma] + [epsilon]/2, 1 + [gamma] + [epsilon]/2; 1). (45)

However, as it was already mentioned above, the case (43) is a rather simple one because the transformation

u(z) = [z.sup.1-[gamma]][(1 - z/a).sup.1-[epsilon]] w(4z (1 - z)) (46)

reduces the Heun equation to the Gauss hypergeometric equation for the new function w. The solution of the general Heun equation is then explicitly written as

u(z) = [z.sup.1-[gamma]][(1 - z/a).sup.1-[epsilon]] x [sub.2][F.sub.1](1 - [alpha] + [delta]/2, 1 - [beta] + [delta]/2; [delta]; 4(1 - z)z). (47)

Now, we will show that there exist nontrivial cases of two-term reductions of the three-term recurrence (11) with (16)-(18) or (22)-(24). These reductions are achieved by the following ansatz guessed by examination of the structure of solutions (33) and (39):

[c.sub.n] = (1/n [[PI].sup.N+2.sub.k+1]([a.sub.k] - 1 + n)/[[PI].sup.N+1.sub.k=1]([b.sub.k] - 1 + n)[c.sub.n-1], (48)

where, having in mind the coefficients [R.sub.n], [Q.sub.n], and [P.sub.n] given by (16)-(18), we put

[a.sub.1],...[a.sub.N], [a.sub.N+1], [a.sub.N+2] = 1 + [e.sub.1],...,1 + [e.sub.N], [gamma] + [epsilon] - [alpha], [gamma] + [epsilon] - [beta], (49)

[b.sub.1],...,[b.sub.N], [b.sub.N+1] = [e.sub.1],...,[e.sub.N], [gamma] + [epsilon] (50)

with parameters [e.sub.1], ..., [e.sub.N] to be defined later. Note that this ansatz implies that [e.sub.1], ..., [e.sub.N] are not zero or negative integers.

The ratio [c.sub.n]/[c.sub.n-1] is explicitly written as

[mathematical expression not reproducible]. (51)

With this, the recurrence relation (11) is rewritten as

[mathematical expression not reproducible]. (52)

Substituting [R.sub.n] and [P.sub.n-2] from (16) and (18) and cancelling the common denominator, this equation becomes

[mathematical expression not reproducible]. (53)

This is a polynomial equation in n. Notably, it is of degree N + 1, not N + 2, because the highest-degree term proportional to [n.sup.N+2] identically vanishes. Hence, we have an equation of the form

[N+1.summation over (m=0)][A.sub.m](a, q; [alpha], [beta], [gamma], [delta], [epsilon]; [e.sub.1], ..., [e.sub.N])[n.sup.m] = 0. (54)

Then, equating to zero the coefficients [A.sub.m] warrants the satisfaction of the three-term recurrence relation (11) for all n. We thus have N + 2 equations [A.sub.m] = 0, m = 0,1,...,N +1,of which N equations serve for determination of the parameters [e.sub.1,2,...,N] and the remaining two impose restrictions on the parameters of the Heun equation.

One of these restrictions is derived by calculating the coefficient [A.sub.N+1] of the term proportional to [n.sup.N+1] which is readily shown to be 2 + N - [delta]. Hence,

[delta] = 2 + N. (55)

The second restriction imposed on the parameters of the Heun equation is checked to be a polynomial equation of the degree N + 1 for the accessory parameter q.

With the help of the Fuchsian condition 1 + [alpha] + [beta] = [gamma] + [delta] + [epsilon], we have

[gamma] + [epsilon] - [alpha] - 1 = [beta] - [delta] = [beta] - 2 - N, (56)

[gamma] + [epsilon] - [beta] - 1 = [alpha] - [delta] = [alpha] - 2 - N, (57)

[gamma] + [epsilon] - 1 = [alpha] + [beta] - [delta] = [alpha] + [beta] - 2 - N, (58)

so that the two-term recurrence relation (51) can be rewritten as ([c.sub.0] = 1)

[mathematical expression not reproducible]. (59)

Note that it follows from this relation, since [alpha], [beta], and [e.sub.1],...,[e.sub.N] are not zero or negative integers, that [c.sub.n] may vanish only if [alpha] is a positive integer such that 0 < [alpha] < 2 + N or [beta] is a positive integer such that 0 < [beta] < 2 + N.

Resolving the recurrence (51), the coefficients of expansion (15)-(18) are explicitly written in terms of the gamma functions as

[mathematical expression not reproducible]. (60)

Here are the explicit solutions of the recurrence relation (11) for N = 0 and N = 1.

N = 0: [delta] = 2, (61)

q = a[hamma] + ([alpha] - 1)([beta] - 1), (62)

[mathematical expression not reproducible]. (63)

N = 1: [delta] = 3, (64)

[mathematical expression not reproducible], (65)

[e.sub.1] = -q + a(1 + [gamma]) - 1 + ([alpha] - 1)([beta] - 1), (66)

[mathematical expression not reproducible]. (67)

These results are readily checked by direct verification of the recurrence relation (11) with coefficients (16)-(18). We conclude by noting that similar explicit solutions can be straightforwardly derived for the expansion (21)-(25) as well.

6. Discussion

Thus, we have presented an expansion of the solutions of the Heun equation in terms of hypergeometric functions having the form [sub.2][F.sub.1]([alpha], [beta]; [[gamma].sub.0] + n; z) with [[gamma].sub.0] = [epsilon] + [gamma] and expansions in terms of functions [sub.2][F.sub.1]([alpha], [beta]; [[gamma].sub.0] - n; z) with [[gamma].sub.0] = [gamma], [alpha], [beta]. For any set of parameters of the Heun equation provided that [gamma] + [epsilon], [gamma], [alpha], [beta] are not all simultaneously integers at least one of these expansions can be applied. Obviously, the expansions are meaningless if [alpha][beta] = 0 since then the involved hypergeometric functions are mere constants and for the solution the summation produces the trivial result u = 0.

The applied technique is readily extended to the four confluent Heun equations. For instance, the solutions of the single- and double-confluent Heun equations using the Kummer confluent hypergeometric functions of the forms [sub.1][F.sub.1]([[alpha].sub.0]+ n; [[gamma].sub.0] + n; [s.sub.0]z), [sub.1][F.sub.1]([[alpha].sub.0] + n; [[gamma].sub.0]; [s.sub.0]z), and [sub.1][F.sub.1]([[alpha].sub.0]; [[gamma].sub.0] + n; [s.sub.0]z) are straightforward. By slight modification, equations of more general type, e.g., of the type discussed by Schmidt [22], can also be considered. In all these cases the termination of the series results in closed-form solutions appreciated in many applications. A representative example is the determination of the exact complete return spectrum of a quantum two-state system excited by a laser pulse of Lorentzian shape and of a double level-crossing frequency detuning [45]. A large set of recent applications of the finite-sum expansions of the biconfluent Heun equation in terms of the Hermite functions to the Schrodinger equation is listed in [37] and references therein.

Regarding the closed-form solutions produced by the presented expansions, this happens in three cases: [epsilon] = -N, [epsilon] + [gamma] - [alpha] = -N, [epsilon] + [gamma] - [beta] = -N, N = 0,1,2,3,.... In each case the general Heun equation admits finite-sum solutions in general at N + 1 choices of the accessory parameter q defined by a polynomial equation of the order of N + 1. The last two choices for [epsilon] result in quasi-polynomial solutions, while, in the first case, when [epsilon] is a negative integer, the solutions involve N + 1 hypergeometric functions generally irreducible to simpler functions. Discussing the termination of this series, we have shown that this is possible if a singularity of the Heun equation is an apparent one. We have further shown that the corresponding finite-sum solution of the general Heun equation has a representation through a single generalized hypergeometric function. The general conclusion suggested by this result is that in any such case the power-series expansion of the Heun function is governed by the two-term recurrence relation obeyed by the power-series for the corresponding generalized hypergeometric function [sub.p][F.sub.q].

There are many examples of application of finite-sum solutions of the Heun equation both in physics and mathematics [46-56], for instance, the solution of a class of free boundary problems occurring in groundwater flow in liquid mechanics and the removal of false singular points of Fuchsian ordinary differential equations in applied mathematics [43]. Another example is the derivation of the third independent exactly solvable hypergeometric potential, after the Eckart and the Poschl-Teller potentials, which is proportional to an energy-independent parameter and has a shape that is independent of this parameter [44]. Some other recent examples can be found in references listed in [14].

Finally, we have shown that there exist infinitely many choices of the involved parameters for which the three-term recurrence relations governing the hypergeometric expansions of the solutions of the general Heun equation are reduced to two-term ones. The coefficients of the expansions are then explicitly expressed in terms of the gamma functions. We have explicitly presented two such cases.

https://doi.org/10.1155/2018/4263678

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research has been conducted within the scope of the International Associated Laboratory IRMAS (CNRS-France & SCS-Armenia). The work has been supported by the Russian-Armenian (Slavonic) University at the expense of the Ministry of Education and Science of the Russian Federation, the Armenian State Committee of Science (SCS Grant no. 18RF-139), the Armenian National Science and Education Fund (ANSEF Grant no. PS-4986), and the project "Leading Russian Research Universities" (Grant no. FTI_24_2016 of the Tomsk Polytechnic University). T. A. Ishkhanyan acknowledges the support from SPIE through a 2017 Optics and Photonics Education Scholarship and thanks the French Embassy in Armenia for a doctoral grant as well as the Agence Universitaire de la Francophonie and Armenian State Committee of Science for a Scientific Mobility grant.

References

[1] K. Heun, "Zur Theorie der Riemann'schen Functionen zweiter Ordnung mit vier Verzweigungspunkten," Mathematische Annalen, vol. 33, no. 2, pp. 161-179, 1888.

[2] A. Ronveaux, Ed., Heuns Differential Equations, Oxford University Press, Oxford, UK, 1995.

[3] S. Y. Slavyanov and W. Lay, Special functions, Oxford Mathematical Monographs, Oxford University Press, Oxford, UK, 2000.

[4] E. Renzi and P. Sammarco, "The hydrodynamics of landslide tsunamis: Current analytical models and future research directions," Landslides , vol. 13, no. 6, pp. 1369-1377, 2016.

[5] M. Renardy, "On the eigenfunctions for Hookean and FENE dumbbell models," Journal of Rheology, vol. 57, no. 5, pp. 1311-1324, 2013.

[6] M. M. Afonso and D. Vincenzi, "Nonlinear elastic polymers in random flow," Journal of Fluid Mechanics, vol. 540, pp. 99-108, 2005.

[7] I. C. Fonseca and K. Bakke, "Quantum Effects on an Atom with a Magnetic Quadrupole Moment in a Region with a Time-Dependent Magnetic Field," Few-Body Systems, vol. 58, no. 1, 2017.

[8] Q. Xie, H. Zhong, M. T. Batchelor, and C. Lee, "The quantum Rabi model: solution and dynamics," Journal of Physics A: Mathematical and General, vol. 50, no. 11, 113001, 40 pages, 2017.

[9] C. A. Downing and M. E. Portnoi, "Massless Dirac fermions in two dimensions: Confinement in nonuniform magnetic fields," Physical Review B: Condensed Matter and Materials Physics, vol. 94, no. 16, 2016.

[10] P. Fiziev and D. Staicova, "Application of the confluent Heun functions for finding the quasinormal modes of nonrotating black holes," Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 84, no. 12, 2011.

[11] D. Batic, D. Mills-Howell, and M. Nowakowski, "Potentials of the Heun class: the triconfluent case," Journal of Mathematical Physics, vol. 56, no. 5, 052106, 17 pages, 2015.

[12] 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.

[13] R. L. Hall and N. Saad, "Exact and approximate solutions of Schrodinger's equation with hyperbolic double-well potentials," The European Physical Journal Plus, vol. 131, no. 8, 2016.

[14] M. Hortacsu, Proc. 13th Regional Conference on Mathematical Physics, World Scientific, Singapore, 2013.

[15] R. Schafke and D. Schmidt, "The connection problem for general linear ordinary differential equations at two regular singular points with applications in the theory of special functions," SIAM Journal on Mathematical Analysis, vol. 11, no. 5, pp. 848-862, 1980.

[16] A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Higher Transcendental Functions, vol. 3, McGraw-Hill, New York, NY, USA, 1955.

[17] L. J. Slater, Generalized Hypergeometric Functions, Cambridge University Press, Cambridge, UK, 1966.

[18] W. N. Bailey, Generalized Hypergeometric Series, Stechert-Hafner Service Agency, 1964.

[19] N. Svartholm, "Die Loosung der Fuchs'schen Differentialgleichung zweiter Ordnung durch Hypergeometrische Polynome," Mathematische Annalen, vol. 116, no. 1, pp. 413-421, 1939.

[20] A. Erdelyi, "The Fuchsian equation of second order with four singularities," Duke Mathematical Journal, vol. 9, pp. 48-58, 1942.

[21] A. Erdelyi, "Certain expansions of solutions of the Heun equation," Quarterly Journal of Mathematics, vol. 15, pp. 62-69, 1944.

[22] D. Schmidt, "Die Losung der linearen Differentialgleichung 2. Ordnung um zwei einfache Singularitaoten durch Reihen nach hypergeometrischen Funktionen," Journal fur die Reine und Angewandte Mathematik, vol. 309, pp. 127-148, 1979.

[23] E. G. Kalnins and J. Miller, "Hypergeometric expansions of Heun polynomials," SIAM Journal on Mathematical Analysis, vol. 22, no. 5, pp. 1450-1459, 1991.

[24] S. Mano, H. Suzuki, and E. Takasugi, "Analytic solutions of the Teukolsky equation and their low frequency expansions," Progress of Theoretical and Experimental Physics, vol. 95, no. 6, pp. 1079-1096, 1996.

[25] B. D. Sleeman and V. B. Kuznetsov, "Heun Functions: Expansions in Series of Hypergeometric Functions," http://dlmf.nist .gov/31.11.

[26] T. Kurth and D. Schmidt, "On the global representation of the solutions of second-order linear differential equations having an irregular singularity of rank one in by series in terms of confluent hypergeometric functions," SIAM Journal on Mathematical Analysis, vol. 17, no. 5, pp. 1086-1103, 1986.

[27] S. Mano and E. Takasugi, "Analytic solutions of the Teukolsky equation and their properties," Progress of Theoretical and Experimental Physics, vol. 97, no. 2, pp. 213-232, 1997.

[28] L. J. El-Jaick and B. D. Figueiredo, "Solutions for confluent and double-confluent Heun equations," Journal of Mathematical Physics, vol. 49, no. 8, 083508, 28 pages, 2008.

[29] T. A. Ishkhanyan and A. M. Ishkhanyan, "Expansions of the solutions to the confluent Heun equation in terms of the Kummer confluent hypergeometric functions," AIP Advances, vol. 4, Article ID 087132, 2014.

[30] E. W. Leaver, "Solutions to a generalized spheroidal wave equation: Teukolsky's equations in general relativity, and the two-center problem in molecular quantum mechanics," Journal of Mathematical Physics, vol. 27, no. 5, pp. 1238-1265, 1986.

[31] L. J. El-Jaick and B. D. Figueiredo, "Confluent Heun equations: convergence of solutions in series of Coulomb wavefunctions," Journal of Physics A: Mathematical and General, vol. 46, no. 8, 085203, 29 pages, 2013.

[32] B. D. B. Figueiredo, "Generalized spheroidal wave equation and limiting cases," Journal of Mathematical Physics, vol. 48, Article ID 013503, 2007.

[33] A. Ishkhanyan, "Incomplete beta-function expansions of the solutions to the confluent Heun equation," Journal of Physics A: Mathematical and General, vol. 38, no. 28, pp. L491-L498, 2005.

[34] E. S. Cheb-Terrab, "Solutions for the general, confluent and biconfluent Heun equations and their connection with Abel equations," Journal of Physics A: Mathematical and General, vol. 37, no. 42, pp. 9923-9949, 2004.

[35] T. A. Ishkhanyan, Y. Pashayan-Leroy, M. R. Gevorgyn, C. Leroy, and A. M. Ishkhanyan, "Expansions of the solutions of the biconfluent Heun equation in terms of incomplete Beta and Gamma functions," Journal of Contemporary Physics (Armenian Academy of Sciences), vol. 51, pp. 229-236, 2016.

[36] A. Hautot, "Sur des combinaisons lineaires d'un nombre fini de fonctions transcendantes comme solutions d'equations differentielles du second ordre," Bulletin de la Societe Royale des Sciences de Liege, vol. 40, pp. 13-23, 1971.

[37] T. A. Ishkhanyan and A. M. Ishkhanyan, "Solutions of the biconfluent Heun equation in terms of the Hermite functions," Annals of Physics, vol. 383, pp. 79-91, 2017.

[38] A. Ishkhanyan and K.-A. Suominen, "New solutions of Heun's general equation," Journal of Physics A: Mathematical and General, vol. 36, no. 5, pp. L81-L85, 2003.

[39] C. Leroy and A. M. Ishkhanyan, "Expansions of the solutions of the confluent Heun equation in terms of the incomplete Beta and the Appell generalized hypergeometric functions," Integral Transforms and Special Functions, vol. 26, no. 6, pp. 451-459, 2015.

[40] J. Letessier, G. Valent, and J. Wimp, "Some differential equations satisfied by hypergeometric functions," in Approximation and computation (West Lafayette, IN, 1993), vol. 119 of Internat. Ser. Numer. Math., pp. 371-381, Birkhauser Boston, Boston, Mass, USA, 1994.

[41] R. S. Maier, "P-symbols, Heun identities, and 3 F2 identities," in Special Functions and Orthogonal Polynomials, vol. 471 of Contemp. Math., pp. 139-159, Amer. Math. Soc., Providence, RI, 2008.

[42] K. Takemura, "Heun's equation, generalized hypergeometric function and exceptional Jacobi polynomial," Journal of Physics A: Mathematical and General, vol. 45, no. 8, 085211, 14 pages, 2012.

[43] A. V. Shanin and R. V. Craster, "Removing false singular points as a method of solving ordinary differential equations," European Journal of Applied Mathematics, vol. 13, no. 6, pp. 617-639, 2002.

[44] A. Ishkhanyan, "The third exactly solvable hypergeometric quantum-mechanical potential," EPL (Europhysics Letters), vol. 115, no. 2, 2016.

[45] A. M. Ishkhanyan and A. E. Grigoryan, "Fifteen classes of solutions of the quantum two-state problem in terms of the confluent Heun function," Journal of Physics A: Mathematical and General, vol. 47, no. 46, 465205, 22 pages, 2014.

[46] L. Carlitz, "Some orthogonal polynomials related to elliptic functions," Duke Mathematical Journal, vol. 27, pp. 443-459, 1960.

[47] K. Kuiken, "Heun's equation and the hypergeometric equation," SIAM Journal on Mathematical Analysis, vol. 10, no. 3, pp. 655-657, 1979.

[48] J. N. Ginocchio, "A class of exactly solvable potentials. I. One-dimensional Schrodinger equation," Annals of Physics, vol. 152, no. 1, pp. 203-219, 1984.

[49] G. Valent, "An integral transform involving Heun functions and a related eigenvalue problem," SIAM Journal on Mathematical Analysis, vol. 17, no. 3, pp. 688-703, 1986.

[50] G. S. Joyce, "On the cubic lattice Green functions," Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 445, no. 1924, pp. 463-477, 1994.

[51] G. S. Joyce and R. T. Delves, "Exact product forms for the simple cubic lattice Green function. I," Journal of Physics A: Mathematical and General, vol. 37, no. 11, pp. 3645-3671, 2004.

[52] G. S. Joyce and R. T. Delves, "Exact product forms for the simple cubic lattice Green function. II," Journal of Physics A: Mathematical and General, vol. 37, no. 20, pp. 5417-5447, 2004.

[53] R. S. Maier, "On reducing the Heun equation to the hypergeometric equation," Journal of Differential Equations, vol. 213, no. 1, pp. 171-203, 2005.

[54] M. van Hoeij and R. Vidunas, "Belyi functions for hyperbolic hypergeometric-to-Heun transformations," Journal of Algebra, vol. 441, pp. 609-659, 2015.

[55] R. Vidunas and G. Filipuk, "Parametric transformations between the Heun and Gauss hypergeometric functions," Funkcialaj Ekvacioj. Serio Internacia, vol. 56, no. 2, pp. 271-321, 2013.

[56] R. Vidunas and G. Filipuk, "A classification of coverings yielding Heun-to-hypergeometric reductions," Osaka Journal of Mathematics, vol. 51, no. 4, pp. 867-903, 2014.

T. A. Ishkhanyan, (1,2,3) T. A. Shahverdyan, (3) and A. M. Ishkhanyan (iD) (1,4)

(1) Russian-Armenian University, H. Emin 123, 0051 Yerevan, Armenia

(2) Moscow Institute of Physics and Technology, Dolgoprudny, Moscow Region 141700, Russia

(3) Institute for Physical Research, NAS of Armenia, Ashtarak 0203, Armenia

(4) Institute of Physics and Technology, National Research Tomsk Polytechnic University, Tomsk 634050, Russia

Correspondence should be addressed to A. M. Ishkhanyan; aishkhanyan@gmail.com

Received 27 January 2018; Accepted 4 May 2018; Published 4 July 2018

Academic Editor: Saber Zarrinkamar

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Ishkhanyan, T.A.; Shahverdyan, T.A.; Ishkhanyan, A.M. |

Publication: | Advances in High Energy Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 6930 |

Previous Article: | Dark Matter in the Standard Model Extension with Singlet Quark. |

Next Article: | Apparent Horizon and Gravitational Thermodynamics of Universe in the Eddington-Born-Infeld Theory. |

Topics: |