We examine the expansions of the solutions of the general Heun equation in terms of the Gauss hypergeometric functions. We present several expansions using functions, the forms of which differ from those applied before. In general, the coefficients of the expansions obey three-term recurrence relations. However, there exist certain choices of the parameters for which the recurrence relations become two-term. The coefficients of the expansions are then explicitly expressed in terms of the gamma functions. Discussing the termination of the presented series, we show that the finite-sum solutions of the general Heun equation in terms of generally irreducible hypergeometric functions have a representation through a single generalized hypergeometric function. Consequently, the power-series expansion of the Heun function for any such case is governed by a two-term recurrence relation.
Ministry of Education and Science of the Russian FederationArmenian State Committee of Science18RF-139Armenian National Science and Education FundPS-4986Leading Russian Research UniversitiesFTI_24_2016SPIE1. Introduction
The 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 F12, 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], Erdélyi [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 F12(λ+n,μ-n;γ;z), which have matching behavior in two singular points, z=0 and 1. These functions have the following Riemann P-symbol representation:(1)P01∞00λ+nz1-γ1-δμ-n.Here λ and μ may adopt several values provided 1+λ+μ=γ+δ (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=∞.
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 F12(α,β;γ0±n;z) which have matching behavior (i.e., characteristic exponents) only at the singular point z=∞. These functions are presented by the Riemann P-symbol(2)P01∞00αz1-γ0±n1-δ0∓nβ,where 1+α+β=γ0+δ0, and the parameters γ0, δ0 are chosen so that the Fuchsian condition for the general Heun equation (3), i.e., 1+α+β=γ+δ+ε, is fulfilled: γ0+δ0=γ+δ+ε. 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 F12(λ+n,μ-n;ν+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, for certain choices of 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 pFq. 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 pFq).
2. Hypergeometric Expansions
The general Heun equation written in its canonical form is [1](3)d2udz2+γz+δz-1+εz-adudz+αβz-qzz-1z-au=0,where the parameters satisfy the Fuchsian relation 1+α+β=γ+δ+ε. We introduce an expansion of this equation’s solution of the form(4)u=∑n=0∞cn·F12α,β;γ0+n;zwith the involved Gauss hypergeometric functions obeying the equation(5)d2undz2+γ0+nz+δ+ε+γ-γ0-nz-1dundz+αβzz-1un=0.Substitution of (4)-(5) into (3) gives(6)∑ncnγ-γ0-nz-ε+γ-γ0-nz-1+εz-adundz+αβa-qzz-1z-aun=0(7)or∑ncna-1ε+γ-γ0-nz-aγ-γ0-nz-1dundz+αβa-qun=0.Now, using the following relations between the involved hypergeometric functions [16–18]:(8)zdundz=γn-1un-1-un,(9)z-1dundz=-δn+1un+δn+1-αβγnun+1,where γn=γ0+n and δn=δ+ε+γ-γn, this equation is rewritten as(10)∑ncna-1ε+γ-γnγn-1un-1-un+aγ-γnδn-1un-δn-1-αβγnun+1+αβa-qun=0,from which we get a three-term recurrence relation for the coefficients of the expansion (4)(11)Rncn+Qn-1cn-1+Pn-2cn-2=0with(12)Rn=a-1ε+γ-γnγn-1,(13)Qn=-Rn+aγ-γnα+β-γn+αβa-q,(14)Pn=aγnγ-γnα-γnβ-γn.From the initial conditions c0=1 and c-1=c-2=0 we get (ε+γ-γ0)(γ0-1)=0. Since γ0=1 is forbidden (it causes division by zero at n=1 in P-1) we obtain that the only possibility is γ0=ε+γ. Hence, the expansion finally reads(15)u=∑n=0∞cn·F21α,β;γ+ε+n;z,and the coefficients of the three-term recurrence relation (11) are(16)Rn=1-anε+γ+n-1,(17)Qn=-Rn+a1+n-δn+ε+aαβ-q,(18)Pn=-an+ε+γn+εn+ε+γ-αn+ε+γ-β.The expansion applies if α,β, and γ+ε are not zero or negative integers. The restrictions on α and β assure that the hypergeometric functions are not polynomials of fixed degree.
The derived expansion terminates if two successive coefficients vanish. If cN is the last nonzero coefficient and cN+1=cN+2=0 for some N=0,1,2,…, we obtain from (11) that it should be PN=0 so that the termination is possible if(19)ε=-Norε+γ-α=-Norε+γ-β=-N.Note that the equation cN+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:(20)Q0R10P0Q1R200P1Q2R3⋱⋱⋱0PN-1QNc0c1c2⋮cN=000⋮0.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(21)u=∑n=0∞cn·F12α,β;γ0-n;z,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(22)Rn=aγ0-nγ-γ0+nα-γ0+nβ-γ0+n,(23)Qn=-Pn+aγ-γ0+nα+β-γ0+n+αβa-q,(24)Pn=a-1ε+γ-γ0+nγ0-n-1,where(25)γ0=γ or α orβ.This expansion applies if α and β are not zero or negative integers and γ is not an integer.
In order for the series to terminate at some n=N we put PN=0 so that this time we derive(26)ε,ε+γ-αor ε+γ-β=-Nfor the expansions with γ0=γ or α or β, respectively. Then the equation cN+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 ε=-N. The involved hypergeometric functions have the form un=F21(α,β,γ-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 γ0=γ:{F21(α,β,γ;z), F21(α,β,γ-1;z),…,F12(α,β,γ-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: ε+γ-α=-N and ε+γ-β=-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 α↔β the finite-sum solutions produced by the choices ε+γ-α=-N and ε+γ-β=-N are of the same form. Furthermore, by applying the formula [16](27)F12α,β;α+k;z=1-zk-βF12k,α-β+k;α+k;zto the involved hypergeometric functions F12(α,β;α-N+n;z) or F12(α,β;β-N+n;z) we see that the sum is a quasi-polynomial, namely, a product of (1-z)1-δ and a polynomial in z. Here are the first two of the solutions:(28)ε+γ-α=0,q=aγδ-1,u=1-z1-δ,(29)ε+γ-α=-1,q2+α-1-aδ-2+γ2δ-3q-aγδ-2α-a1+γδ-1=0,(30)u=1-z1-δ1-α+1-δα-1z+q-aαβ+ε-δε1-aα-11-z.
Note that, since 1-δ is a characteristic exponent of the Heun equation, the transformation u=(1-z)1-δw(z) results in another Heun equation for w(z). Hence, the derived finite-sum solutions corresponding to the choices ε+γ-α=-N and ε+γ-β=-N are generated from the polynomial solutions of the equation for w(z).
More interesting is the case ε=-N, when the finite-sum solutions involve N+1 hypergeometric functions irreducible, in general, to simpler functions. The case ε=0 produces the trivial result q=aαβ, when the Heun equation is degenerated into the hypergeometric equation with the solution u=F12(α,β;γ;z). The solution for the first nontrivial case ε=-1 reads(31)u=F12α,β;γ-1;z+q-aαβ+a1-δ1-aγ-1·F12α,β;γ;z,where q is a root of the equation(32)q-aαβ+a1-δq-aαβ+a-11-γ-a1-a1+α-γ1+β-γ=0.Note that the second term in (31) vanishes if q-aαβ+a(1-δ)=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, γ+δ=2 and α or β equals γ-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 ε=-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 F23 with an upper parameter exceeding a lower one by unity [40–42]:(33)uu0=F23α,β,e+1;γ,e;z,where the parameter e is given as(34)e=aαβq-aαβ.Note that using this parameter e the solution of (32) is parameterized as [41](35)q=aαβ1+ee,a=ee-γ+1e-αe-β.We will now show that a similar generalized hypergeometric representation holds also for ε=-2 and for all ε∈Z, ε≠1.
4. The Case ε≤-2, ε∈Z
For ε=-2 the termination equation cN+1=0 for the accessory parameter q is written as(36)q-aαβ2+q-aαβ4a-2-3+α+βa+γ+2aa-1αβ×q-aαβ-21+α+βa-2+2γ+q-aαβ2aa-1αβ+1+α+β=0.The solution of the Heun equation for a root of this equation is given as(37)u=F12α,β;γ-2;z+B1·F12α,β;γ-1;z+B2·F12α,β;γ;z,with(38)B1=q-aαβ+a1-δ1-aγ-2,B2=a1+α-γ1+β-γq-aαβ+21-aγ-1γ-1B1.It is now checked that this solution is presented by the hypergeometric function F43 as [42](39)uu0=F43α,β,e+1,r+1;γ,e,r;z,where u0=1+B1+B2 and the parameters e, r solve the equations (compare with (35))(40)q=aαβ1+e1+rer,(41)a=er2er+e+r+12-γαβ1+e2+1+r2-1+er2er-4-e+r+3α+β-1.It is further checked that this system of equations admits a unique solution e, r (up to the transposition e↔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 z0, z1, and z2 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 ε=-2. This is straightforwardly verified by checking the power-series solution u=∑n=0∞cn(z-a)n with c0≠0 for the neighborhood of the point z=a. In calculating c3 a division by zero will occur, unless q satisfies (36), in which case the equation for c3 will be identically satisfied.
It can be checked that generalized hypergeometric representations are achieved also for ε=-3,-4,-5 [42]. The conjecture is that for any negative integer ε=-N, N=1,2,3,… there exists a generalized hypergeometric solution of the Heun equation given by the ansatz(42)u=FN+21+Nα,β,e1+1,…,eN+1;γ,e1,…,eN;zprovided 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 ε=0, for which the Heun function degenerates to the Gauss hypergeometric function u=F21(α,β;γ;z) provided q=0.
Finally, we note that by the elementary power change u=(z-a)1-εw a Heun equation with a positive exponent parameter ε>2 is transformed into the one with a negative parameter 2-ε<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 ε=N, N=2,3,…. Thus, the only exception is the case ε=1.
It is a basic knowledge that the generalized hypergeometric function pFq 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 pFq, 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 pFq).
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(43)a=1/2,γ+δ=2and q=aαβ+a1-δε,when the coefficient Qn 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(44)u=∑k=0∞ε/2kγ+ε-α/2kγ+ε-β/2kk!γ+ε/2k1+γ+ε/2k·F21α,β;γ+ε+2k;z,where (…)k denotes the Pochhammer symbol. The values u(0), u′(0) and u(1), and u′(1) can then be written in terms of generalized hypergeometric series [17, 18]. For instance,(45)u0=F32γ+ε-α2,γ+ε-β2,ε2;γ+ε2,1+γ+ε2;1.
However, as it was already mentioned above, the case (43) is a rather simple one because the transformation(46)uz=z1-γ1-za1-εw4z1-zreduces 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(47)uz=z1-γ1-za1-εF211-α+δ2,1-β+δ2;δ;41-zz.
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):(48)cn=1n∏k=1N+2ak-1+n∏k=1N+1bk-1+ncn-1,where, having in mind the coefficients Rn, Qn, and Pn given by (16)-(18), we put(49)a1,…,aN,aN+1,aN+2=1+e1,…,1+eN,γ+ε-α,γ+ε-β,(50)b1,…,bN,bN+1=e1,…,eN,γ+εwith parameters e1,…,eN to be defined later. Note that this ansatz implies that e1,…,eN are not zero or negative integers.
The ratio cn/cn-1 is explicitly written as(51)cncn-1=γ+ε-α-1+nγ+ε-α-1+nγ+ε-1+nn∏k=1Nek+nek-1+n.With this, the recurrence relation (11) is rewritten as(52)Rnγ+ε-α-1+nγ+ε-β-1+nγ+ε-1+nn∏k=1Nek+nek-1+n+Qn-1+Pn-2γ+ε-2+nn-1γ+ε-α-2+nγ+ε-β-2+n∏k=1Nek-2+nek-1+n=0.Substituting Rn and Pn-2 from (16) and (18) and cancelling the common denominator, this equation becomes(53)1-aγ+ε-α-1+nγ+ε-β-1+n∏k=1Nek+n+Qn-1∏k=1Nek-1+n-aε+n-2n-1∏k=1Nek-2+n=0.This is a polynomial equation in n. Notably, it is of degree N+1, not N+2, because the highest-degree term proportional to nN+2 identically vanishes. Hence, we have an equation of the form(54)∑m=0N+1Ama,q;α,β,γ,δ,ε;e1,…,eNnm=0.Then, equating to zero the coefficients Am warrants the satisfaction of the three-term recurrence relation (11) for all n. We thus have N+2 equations Am=0, m=0,1,…,N+1, of which N equations serve for determination of the parameters e1,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 AN+1 of the term proportional to nN+1 which is readily shown to be 2+N-δ. Hence,(55)δ=2+N.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+α+β=γ+δ+ε, we have(56)γ+ε-α-1=β-δ=β-2-N,(57)γ+ε-β-1=α-δ=α-2-N,(58)γ+ε-1=α+β-δ=α+β-2-N,so that the two-term recurrence relation (51) can be rewritten as (c0=1)(59)cn=α-2-N+nβ-2-N+nα+β-2-N+nn∏k=1Nek+nek-1+ncn-1,n≥1.Note that it follows from this relation, since α, β, and e1,…,eN are not zero or negative integers, that cn may vanish only if α is a positive integer such that 0<α<2+N or β is a positive integer such that 0<β<2+N.
Resolving the recurrence (51), the coefficients of expansion (15)-(18) are explicitly written in terms of the gamma functions as(60)cn=Γγ+εΓn+γ+ε-αΓn+γ+ε-βn!Γγ+ε-αΓγ+ε-βΓn+γ+ε∏k=1Nek+nek,n≥1.
Here are the explicit solutions of the recurrence relation (11) for N=0 and N=1.(61)N=0:δ=2,(62)q=aγ+α-1β-1,(63)cn=Γγ+εΓn+γ+ε-αΓn+γ+ε-βn!Γγ+ε-αΓγ+ε-βΓn+γ+ε.(64)N=1:δ=3,(65)q2-q4+a-3α-3β+2αβ+3aγ+α-2α-1β-2β-1+a4+2a-4α-4β+3αβγ+2a2γ2=0,(66)e1=-q+a1+γ-1+α-1β-1,(67)cn=Γγ+εΓn+γ+ε-αΓn+γ+ε-βn!Γγ+ε-αΓγ+ε-βΓn+γ+εe1+ne1.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 F21(α,β;γ0+n;z) with γ0=ε+γ and expansions in terms of functions F21(α,β;γ0-n;z) with γ0=γ,α,β. For any set of parameters of the Heun equation provided that γ+ε, γ,α,β are not all simultaneously integers at least one of these expansions can be applied. Obviously, the expansions are meaningless if αβ=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 F11(α0+n;γ0+n;s0z), F11(α0+n;γ0;s0z), and F11(α0;γ0+n;s0z) 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 Schrödinger equation is listed in [37] and references therein.
Regarding the closed-form solutions produced by the presented expansions, this happens in three cases: ε=-N, ε+γ-α=-N, ε+γ-β=-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 ε result in quasi-polynomial solutions, while, in the first case, when ε 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 pFq.
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 Pöschl-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.
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.
HeunK.Zur Theorie der Riemann'schen Functionen zweiter Ordnung mit vier Verzweigungspunkten188833216117910.1007/BF01443849MR1510535RonveauxA.1995Oxford, UKOxford University Press10.1016/j.amc.2006.11.170Zbl1126.34307SlavyanovS. Y.LayW.2000Oxford, UKOxford University PressOxford Mathematical MonographsMR1858237RenziE.SammarcoP.The hydrodynamics of landslide tsunamis: Current analytical models and future research directions2016136136913772-s2.0-8500384866710.1007/s10346-016-0680-zRenardyM.On the eigenfunctions for Hookean and FENE dumbbell models2013575131113242-s2.0-8488376551910.1122/1.4816631AfonsoM. M.VincenziD.Nonlinear elastic polymers in random flow2005540991082-s2.0-2694448146210.1017/S0022112005005951Zbl1082.76006FonsecaI. C.BakkeK.Quantum Effects on an Atom with a Magnetic Quadrupole Moment in a Region with a Time-Dependent Magnetic Field20175812-s2.0-85000730411XieQ.ZhongH.BatchelorM. T.LeeC.The quantum Rabi model: solution and dynamics20175011113001, 4010.1088/1751-8121/aa5a65MR3622572Zbl1362.81105DowningC. A.PortnoiM. E.Massless Dirac fermions in two dimensions: Confinement in nonuniform magnetic fields2016941610.1103/PhysRevB.94.165407FizievP.StaicovaD.Application of the confluent Heun functions for finding the quasinormal modes of nonrotating black holes201184122-s2.0-84855327131BaticD.Mills-HowellD.NowakowskiM.Potentials of the Heun class: the triconfluent case2015565052106, 1710.1063/1.4921344MR3390991Zbl1319.81049VieiraH. S.BezerraV. B.Confluent Heun functions and the physics of black holes: resonant frequencies, HAWking radiation and scattering of scalar waves2016373284210.1016/j.aop.2016.06.016MR3548887Zbl1380.831672-s2.0-84978977367HallR. L.SaadN.Exact and approximate solutions of Schrödinger’s equation with hyperbolic double-well potentials2016131810.1140/epjp/i2016-16277-1HortacsuM.2013World Scientific, SingaporeSchäfkeR.SchmidtD.The connection problem for general linear ordinary differential equations at two regular singular points with applications in the theory of special functions198011584886210.1137/0511076MR586913ErdélyiA.MagnusW.OberhettingerF.TricomiF. G.19553New York, NY, USAMcGraw-HillSlaterL. J.1966Cambridge, UKCambridge University PressMR0201688BaileyW. N.1964Stechert-Hafner Service AgencySvartholmN.Die Lösung der Fuchs'schen Differentialgleichung zweiter Ordnung durch Hypergeometrische Polynome1939116141342110.1007/BF01597365MR1513235ErdélyiA.The Fuchsian equation of second order with four singularities19429485810.1215/S0012-7094-42-00906-2MR0005981ErdélyiA.Certain expansions of solutions of the Heun equation1944156269MR0011506SchmidtD.Die Lösung der linearen Differentialgleichung 2. Ordnung um zwei einfache Singularitäten durch Reihen nach hypergeometrischen Funktionen197930912714810.1515/crll.1979.309.127MR542043KalninsE. G.MillerJ.Hypergeometric expansions of Heun polynomials19912251450145910.1137/0522093MR1112519ManoS.SuzukiH.TakasugiE.Analytic solutions of the Teukolsky equation and their low frequency expansions19969561079109610.1143/PTP.95.1079MR14019742-s2.0-0030499647SleemanB. D.KuznetsovV. B.Heun Functions: Expansions in Series of Hypergeometric Functionshttp://dlmf.nist.gov/31.11KurthT.SchmidtD.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 functions19861751086110310.1137/0517077MR853518ManoS.TakasugiE.Analytic solutions of the Teukolsky equation and their properties199797221323210.1143/PTP.97.213MR14471612-s2.0-0031506340El-JaickL. J.FigueiredoB. D.Solutions for confluent and double-confluent Heun equations2008498083508, 2810.1063/1.2970150MR2440706Zbl1152.81421IshkhanyanT. A.IshkhanyanA. M.Expansions of the solutions to the confluent Heun equation in terms of the Kummer confluent hypergeometric functions20144087132LeaverE. W.Solutions to a generalized spheroidal wave equation: Teukolsky's equations in general relativity, and the two-center problem in molecular quantum mechanics19862751238126510.1063/1.527130MR836023Zbl0612.340412-s2.0-0001528256El-JaickL. J.FigueiredoB. D.Confluent Heun equations: convergence of solutions in series of Coulomb wavefunctions2013468085203, 2910.1088/1751-8113/46/8/085203MR3030189Zbl1264.81175FigueiredoB. D. B.Generalized spheroidal wave equation and limiting cases200748013503IshkhanyanA.Incomplete beta-function expansions of the solutions to the confluent Heun equation20053828L491L49810.1088/0305-4470/38/28/L02MR2166618Zbl1080.34571Cheb-TerrabE. S.Solutions for the general, confluent and biconfluent Heun equations and their connection with Abel equations200437429923994910.1088/0305-4470/37/42/007MR2098526Zbl1072.34002IshkhanyanT. A.Pashayan-LeroyY.GevorgynM. R.LeroyC.IshkhanyanA. M.Expansions of the solutions of the biconfluent Heun equation in terms of incomplete Beta and Gamma functions201651229236HautotA.Sur des combinaisons linéaires d'un nombre fini de fonctions transcendantes comme solutions d'équations différentielles du second ordre1971401323MR0299853IshkhanyanT. A.IshkhanyanA. M.Solutions of the bi-confluent Heun equation in terms of the Hermite functions2017383799110.1016/j.aop.2017.04.015MR3682014IshkhanyanA.SuominenK.-A.New solutions of Heun's general equation2003365L81L8510.1088/0305-4470/36/5/101MR1960079Zbl1066.34088LeroyC.IshkhanyanA. M.Expansions of the solutions of the confluent Heun equation in terms of the incomplete Beta and the Appell generalized hypergeometric functions201526645145910.1080/10652469.2015.1019490MR3327473Zbl1326.33032LetessierJ.ValentG.WimpJ.Some differential equations satisfied by hypergeometric functions1994119Boston, Mass, USABirkhauser Boston371381Internat. Ser. Numer. Math.MR1333630MaierR. S.P-symbols, Heun identities, and 3 F2 identities2008471Amer. Math. Soc., Providence, RI139159Contemp. Math.10.1090/conm/471/09211MR2367417TakemuraK.Heun's equation, generalized hypergeometric function and exceptional Jacobi polynomial2012458085211, 1410.1088/1751-8113/45/8/085211MR2897019Zbl1247.34132ShaninA. V.CrasterR. V.Removing false singular points as a method of solving ordinary differential equations200213661763910.1017/S0956792502004916MR1949726Zbl1022.340782-s2.0-0036904148IshkhanyanA.The third exactly solvable hypergeometric quantum-mechanical potential201611522-s2.0-84984874847IshkhanyanA. M.GrigoryanA. E.Fifteen classes of solutions of the quantum two-state problem in terms of the confluent Heun function20144746465205, 2210.1088/1751-8113/47/46/465205MR3279984Zbl1307.81081CarlitzL.Some orthogonal polynomials related to elliptic functions196027443459MR012303010.1215/S0012-7094-60-02742-3Zbl0096.269032-s2.0-84972526614KuikenK.Heun's equation and the hypergeometric equation197910365565710.1137/0510060MR529081Zbl0415.34038GinocchioJ. N.A class of exactly solvable potentials. I. One-dimensional Schrödinger equation1984152120321910.1016/0003-4916(84)90084-8MR734212Zbl0545.340182-s2.0-0000682647ValentG.An integral transform involving Heun functions and a related eigenvalue problem198617368870310.1137/0517049MR838248Zbl0604.44004JoyceG. S.On the cubic lattice Green functions1994445192446347710.1098/rspa.1994.0072MR1276913Zbl0808.33015JoyceG. S.DelvesR. T.Exact product forms for the simple cubic lattice Green function. I200437113645367110.1088/0305-4470/37/11/008MR2042148JoyceG. S.DelvesR. T.Exact product forms for the simple cubic lattice Green function. II200437205417544710.1088/0305-4470/37/20/012MR2065679MaierR. S.On reducing the Heun equation to the hypergeometric equation2005213117120310.1016/j.jde.2004.07.020MR2139342Zbl1085.340352-s2.0-18144366618van HoeijM.VidunasR.Belyi functions for hyperbolic hypergeometric-to-Heun transformations201544160965910.1016/j.jalgebra.2015.06.008MR3391940VidunasR.FilipukG.Parametric transformations between the Heun and Gauss hypergeometric functions201356227132110.1619/fesi.56.271MR3114824Zbl1279.33030VidunasR.FilipukG.A classification of coverings yielding Heun-to-hypergeometric reductions2014514867903MR3273870Zbl1309.330242-s2.0-84908314190