Generalized-hypergeometric solutions of the general Fuchsian linear ODE having five regular singularities
Axioms , Article
Generalized-Hypergeometric Solutions of the General Fuchsian Linear ODE Having Five Regular Singularities
Artur Ishkhanyan and Clemente Cesarano * Department of General Physics, Russian-Armenian University, Yerevan 0051, Armenia Matter Wave Physics Department, Institute for Physical Research, Ashtarak0203, Armenia Section of Mathematics-International Telematic University Uninettuno, C.so Vittorio Emanuele II, 39, 00186 Roma, Italy * Correspondence: [email protected] Received: 23 June 2019; Accepted: 27 August 2019; Published: 2 September 2019
Abstract:
We show that a Fuchsian differential equation having five regular singular points admits solutions in terms of a single generalized hypergeometric function for infinitely many particular choices of equation parameters. Each solution assumes four restrictions imposed on the parameters: two of the singularities should have non-zero integer characteristic exponents and the accessory parameters should obey polynomial equations.
Keywords:
Fuchsian equation; generalized hypergeometric function; recurrence relation
1. Introduction
During the last decades, the five Heun equations and their generalizations became a subject of intensive investigations in the context of numerous advanced problems of contemporary fundamental and applied research. Being natural generalizations of the equations of the hypergeometric class, these equations frequently appear, for instance, in contemporary classical and non-classical physics. Some representative examples include nuclear reactions, neutrino oscillations, quantum mechanics and quantum information theory, nanophysics, ultralow-temperature physics (Bose–Einstein condensation), fluid mechanics, statistical physics, astrophysics, cosmology, etc. (see, e.g., [1,2] and references therein)
The general Heun equation, from which the other equations originate, is the most general
Fuchsian second-order linear differential equation having four regular singular points [3,4]. Here we consider the linear ordinary differential equation:
22 1 012 1 1
01 ( 1)( )( ) z zd u du uz z z a z a dz z z z a z adz αβ θ θεγ δ ε − −+ + + + + = − − − − − − (1) which presents the most general Fuchsian equation having five regular singular points [3,4]. The singularities are located at z a a = and at z = ∞ . The regularity of the singularity at infinity is ensured by the Fuchsian condition α β γ δ ε ε+ + = + + + (2) If ε = and / a a θ αβ θ= − , Equation (1) becomes the general Heun equation [5] xioms , , 102 2 of 8
01 ( 1)( ) d u du z q uz z z a dz z z z adz γ δ ε αβ − + + + + = − − − − (3) with / q a θ= − , which is the most general Fuchsian equation having four regular singular points [6–8]. The parameter q which is not related to the characteristic exponents of the singularities, is referred to as the accessory parameter. For Equation (1), there are two accessory parameters − θ and θ . There are many other choices of parameters that reduce Equation (1) to the general Heun Equation (3). For instance, an obvious choice is a ε θ= = = . This and the previous reduction cases can be made symmetric if the parameters θ and θ are expressed through new parameters q and q such that q q θ = + and qq θ = . Then, for the triads ( , , ) (0, 0, 0) a q ε = and ( , , ) (0, 0, 0) a q ε = Equation (1) becomes the general Heun equation with the third finite singularity being located at z a = and z a = , respectively. For this reason, in this case, one may conventionally think of the parameters q and q as accessory parameters associated with the singularities a and a , respectively. To find solutions of Equation (1) in terms of the generalized hypergeometric functions [9,10], we apply the approach developed in [11,12] for Heun-type equations which originate, via coalescence of singularities, from the general Heun equation [6–8]. A main result we report in the present paper is that there exist infinitely many particular choices of parameters for which Equation (1) having five irreducible regular singularities admits solutions in terms of a single generalized-hypergeometric function. We should keep in the mind the possible reductions of Equation (1) to the Heun equation in order to, in due course, identify the cases when the constructed solutions in fact apply to the general Heun equation, not to Equation (1) with irreducible singularities. These reducible cases are not taken into account.
2. The Approach
To construct the solutions, we examine the cases, when a power-series Frobenius solution [13] of Equation (1) reduces to the generalized hypergeometric series [9,10] ( ) ,..., ; ,..., ; nr s r s nn cF a a b b z z ∞= = (4) for which the coefficients obey the two-term recurrence relation ( )( )
11 1
11 1 k kkrn sn k c n a nnc b =− = − +− += ∏∏ (5) The Frobenius solution of Equation (1) for the vicinity of the regular singularity located at z = is given as ( ) nnn u z z zc μ ∞= = (6) where the characteristic exponent μ may adopt two values: μ= or γ− , and the expansion coefficients obey a four -term recurrence relation [3,4]: n n n n n n n n R c Q c P c S c − − − − − − + + + = (7) For the exponent μ= the coefficients are given as ( 1) n R aa n n γ= + − (8) xioms , , 102 3 of 8 ( ) ( 1)( ) n Q nn a a aa a a aa θ γ ε ε δ− + − + + += + + (9) ( ) ( )(1 ) n n a a a a nP θ α β δ ε ε+ + + + − − −= + (10) ( )( ) n S n n α β= − + + (11) We note that the series (6) with such coefficients may terminate only if α or β is a non-positive integer (besides, θ and θ should obey certain polynomial equations). Having this in the mind and based on the experience gained in treating the general and single-confluent Heun equations [11,12], we try the ansatz (5) for r N = + , s N = + , and ,...., , , , ,1 ,...,1 N N N N a a a a e e α β + + = + + (12) ,...., , , ,...,
N N N b b b e e γ + = (13) This means that we look for a solution of Equation (1) as ( , ,1 ,...,1 ; , ,..., ; ) N N N N u F e e e e z α β γ σ + + = + + (14) where σ is a scaling factor to be defined later. The series representation for this function is governed by a two term recurrence relation for coefficients explicitly written as ( )( )( )
11 101
N Nk n kk n kkn kn c nsn n e nn S e nn e ec n R α βγσ − −= = − + − + + +=− + − + − += ∏ ∏ (15) where s aa σ= − . With this, the recurrence relation (7) is rewritten as ( ) ( )( ) ( )
12 0.3
N Nn n k n kk kN Nn n k n n knk k s S S e n s S e ns R P e n R R e nQ − − −= =− − −− −= = + − +− + + +− ++ = ∏ ∏∏ ∏ (16) This is a polynomial equation in summation index n of the degree utmost N n + . To ensure the fulfillment of this equation for all n , we demand the coefficients to identically vanish. For the coefficient of the highest order term we have ( )( )( ) a s a s aa s + + + = , hence, , , s a a aa = − − − (17) The equation for the term proportional to N n + reads ( ) ( )( ) a a a a N ε ε− − + + = , hence, since , 0,1 a a ≠ , N ε ε+ = − (18) Calculating the sum of all coefficients, we derive another relation: ( )
10 0 1 N kk N kk es e θ αβ = = += ∏∏ (19) Finally, one more relation is derived if one checks the value adopted by the product ( ) Nk k ε = + ∏ . Notably, the product turns to be zero. In the force of Equation (18), this means that ε and ε are non-positive integers such that N ε ε+ = − . With these observations, Equation (16) is reduced to a polynomial equation for n of the degree N : mN mm A n = = , where m A are functions of the parameters of the starting Equation (1) and the parameters θ and ,..., N e e . Resolving then the equations m A = , m N = , we arrive at a polynomial equation of the degree xioms , , 102 4 of 8 N + for θ and the corresponding parameters ,..., N e e for the generalized-hypergeometric solution (14) of Equation (1).
3. Explicit Solutions of Lowest Orders N = Consider the choice s aa = − , that is σ= . The case N = is almost trivial: N = : ε ε= = , ( ) a a θ αβ= + , aa θ αβ= − , ( , ; ; ) u F z α β γ= (20) In this case the parameters θ and θ are such that the singularities at z a = and z a = both disappear so that Equation (1) is reduced to the ordinary Gauss hypergeometric equation. For N = we have two choices: ( , ) ( 1, 0) ε ε = − and ( , ) (0, 1) ε ε = − , both reproducing the result for the general Heun equation [14–18,12]. Here is the explicit solution for the first choice: N = : ε = − , ε = : ( , ,1 ; , ; ) u F e e z α β γ= + (21)
01 1 1 / a a αβ θθ −= ,
10 1 1 eaa e θ αβ += − (22) ( )( ) ( ) a e e e e α β γ− − + − =− (23) Because of the particular value of θ , in this case the singularity at z a = disappears so that Equation (1) is reduced to the general Heun equation. Resolving the second Equation (22) with respect to e and substituting the result into (23), we arrive at a quadratic equation for θ , which provides two values of this accessory parameter for which the solution (21) applies. This equation reads ( ) ( )
20 02 11 a a aaa θ θα β αβ γ αβ α β γ+ + + + − + + + − = (24) For N = there are three possibilities: ( , ) ( 2, 0) ε ε = − , ( , ) (0, 2) ε ε = − , and ( , ) ( 1, 1) ε ε = − − . For the first two cases a singularity disappears and the result reproduces that for the general Heun equation. For the first choice the solution reads [12,18] N = : ε= − , ε = : ( , ,1 ,1 ; , , ; ) u F e e e e z α β γ= + + (25) with parameters , e e being defined by equations e eaa e e θ αβ + += − , ( )( )( )
22 121 k k kk k k kk e ea e e γα β == = = − += − − ∏∏ (26) This solution applies if
01 1 1 / a a αβ θθ −= and θ obeys the cubic equation ( ) ( ) ( ) ( ) ( ) ( ) ( )( )( ) a a aa a a ω ω γ α β αβω α β γ ω α β− − + − − + − ×+ + + + − + − + + = (27) where / a a ω αβ θ= + . The third choice ( , ) ( 1, 1) ε ε = − − suggests a non-trivial solution again in terms of the functions F : N = : ε= − , ε = − : ( , ,1 ,1 ; , , ; ) u F e e e e z α β γ= + + (28) xioms , , 102 5 of 8 ( ) ( ) ( ) ( )( )
21 11 1 21 k kkk kk ea a a a e θ α γ αα β δ α ==== + − + − + + − − − + − = ∏∏ (29) e eaa e e θ αβ + += − (30) ( )( ) ( ) ( )( ) ( ) ( )( ) k kk kk kk kk kk k e ea a e e α α γ α β β γ βα β δ α β = == == == = + − + − + − + −− − − = −− − ∏ ∏∏ ∏ (31) ( )( )( )
211 21 k kk kkkkk e eaa e e γα β ==== − +−= − ∏∏ (32) where Equations (29) and (30) parameterize θ , θ through the parameters e and e , while the last two equations define the latter parameters. Eliminating , e e from Equations (29) and (30), it can be shown that θ can be presented as a rational function of θ : ( ) ( )( )( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) a aaa a a a a a a aaa a a a aa a a a a θ γ α βα α β β γ β γ α βα β α β γ θα β α β γγ β α β θ θ− + + + − + + + + + − + + + + + + + + − + + = + + + − + + − + + + + + + + (33) while θ obeys a (rather cumbersome) fourth -order polynomial equation the coefficients of which are functions of the parameters , , , , a a α γ δ : A A A A θ θ θ θ+ + + + = (34) with ( ) ( )( ) ( )
13 1
A aa a a αβ α β γ+ + += + − − (35) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( )
A a aa a a a a aa a a a α α β α α β α α βγ α β β α β β γγ αβ α β αβ α β γ γ+ + + + + + + + −+ + + + + + + + + − +− + + + + + + + − += + (36) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ( ) ( ) ( )( ) ( ) ( ) ( ) ( ) ( ( (2 1 1 2 2 1 4 7 12 3 5 3 2 5 7 3 3 12 1 1 1 2 )2 3 5 3 2 5 7 3 8 1 3 1 8 6 26 4 4 2 3 2 1 1 2 )1 2 aA aa a a aa aa a aa a aa a α β α α β α α β α α βγ α α β α α β α αα β αβ α β γ γ αβ α βγ α β β α β β β β βγ αβ α β α β αβ α β γγ γ αβ α β γ+ + + + + + + + + −+ + + + + + + ++ + + + + − + + −+ + + + + + + + + + + ++= + + + + + + + + + − +− + + − ( ) + (37) xioms , , 102 6 of 8 ( )( ) ( ) ( ) ( )( ) ( ) ( ) ( ) ( )( ) ( )( ) ( ) ( )( ) ( ) ( 1 1 1( 3 4 4 4 1 1 2 3 3 4 21 1 2 1 1 2 )). aA a a a aa a a aa a αβ γ α β γγ γ αβ α β α β αβ α β γ γα β α β γ α β γ+ + − + ++ + + − + + + + + + − + +− + + + + − + −= + (38) This result has been checked by Mathematica.
4. Discussion
Thus, we have shown that there exist infinitely many particular choices of parameters for which a Fuchsian differential equation having five regular singular points admits solutions in terms of a single generalized-hypergeometric function. This is the main result of this paper. It should be noted that there have been several studies towards the solution of Fuchsian equations with more than four singularities (these are more complicated equations than the general Heun equation) in terms of simpler mathematical functions of the hypergeometric class (see, e.g., [18–24]). In particular, based on Scheffé ′ s lemma [19], several such solutions have been derived in [20,21]. A computer algebra solver to find solutions having rational function arguments has been presented in [22] (note that the solver works with equations for which the coefficients do not involve continuously variable parameters). The results we have reported here, that is the closed-form explicit solutions for equations having five regular singular points, are in line with similar results derived for the Heun-type equations [11,12]. The results are in agreement with the conjecture by Takemura [18] which suggests that generalized-hypergeometric solutions exist for any Fuchsian differential equation having M + regular singularities of which M are apparent [25]. As we have seen, in our case each generalized-hypergeometric solution assumes four restrictions imposed on the equation parameters, namely, two of the singularities should have non-zero integer characteristic exponents and two accessory parameters should obey certain polynomial equations. It can be shown that the latter restrictions warrant that the two singularities with non-zero integer characteristic exponents are indeed apparent. Fuchsian linear ordinary differential equations having five regular singularities appear in many branches of contemporary physics and mathematics research (see, e.g., [1,26,27]). Based on the experience gained in application of such solutions in the case of the Heun-type equations (see, e.g., [28,29]), one may envisage many applications of these solutions in the future. In particular, one may apply the above solution for N = to generate a new exactly solvable potential for the one-dimensional stationary Schrödinger equation. Other possibilities may appear if the technique used to derive the above-presented results is extended to involve the multivariable generalizations of the hypergeometric functions [30–33]. Finally, taking into account the relationship of the Fuchsian equations having five regular singularities with the Painlevé I-VI equations, applications to nonlinear problems are also possible [34,35]. Author Contributions:
The authors contributed equally to this work. The authors read and approved the final version of the manuscript.
Funding:
This research received no external funding.
Conflicts of Interest:
The authors declare no conflict of interest.
References The
Heun
Project:
Heun
Functions,
Their
Generalizations and
Applications.
Available online : https://theheunproject.org/bibliography.html (accessed on 30 August 2019). 2. Hortaçsu, M. Heun functions and some of their applications in physics. Adv. High Energy Phys. , , Fuchs, L. Zur
Theorie der linearen
Differentialgleichungen mit veränderlichen
Coefficienten.
J. Reine Angew. Math. , , xioms , , 102 7 of 8 Fuchs, L. Zur
Theorie der linearen
Differentialgleichungen mit veränderlichen
Coefficienten. (Ergänzungen zu der im Bande dieses
Journals enthaltenen
Abhandlung).
J. Reine Angew. Math. , , Heun, K. Zur
Theorie der
Riemann’schen
Functionen
Zweiter
Ordnung mit
Verzweigungspunkten.
Math. Ann. , , Ronveaux, A. Heun’s Differential Equations ; Oxford
University
Press:
London,
UK,
Slavyanov,
S.Y.;
Lay, W. Special Functions ; Oxford
University
Press:
Oxford,
UK,
Olver,
F.W.J.;
Lozier,
D.W.;
Boisvert,
R.F.;
Clark,
C.W. (Eds.)
NIST Handbook of Mathematical Functions ; Cambridge
University
Press:
New
York,
NY, USA,
Slater,
L.J.
Generalized Hypergeometric Functions ; Cambridge
University
Press:
Cambridge,
UK,
Bailey,
W.N.
Generalized Hypergeometric Series ; Stechert-Hafner
Service
Agency:
London,
UK,
Ishkhanyan,
T.A.;
Ishkhanyan,
A.M.
Generalized confluent hypergeometric solutions of the Heun confluent equation.
Appl. Math. Comput. , , Ishkhanyan,
A.M.
Generalized hypergeometric solutions of the Heun equation . arXiv , arXiv:1802.04263. 13.
Frobenius,
F.G.
Ueber die
Integration der linearen
Differentialgleichungen durch
Reihen.
J. Reine Angew. Math. , , Letessier, J. Co-recursive associated
Jacobi polynomials.
J. Comput. Appl. Math. , , Letessier,
J.;
Valent,
G.;
Wimp, J. Some differential equations satisfied by hypergeometric functions. Intern. Ser. Numer. Math. , , Maier,
R.S.
P-symbols, Heun Identities, and 3F2 Identities.
Contemp. Math. , , Ishkhanyan,
T.A.;
Shahverdyan,
T.A.;
Ishkhanyan,
A.M.
Expansions of the solutions of the general Heun equation governed by two-term recurrence relations for coefficients. Adv. High Energy Phys. , , Takemura, K. Heun’s equation, generalized hypergeometric function and exceptional
Jacobi polynomial.
J. Phys. A , , Scheffé, H. Linear differential equations with two-term recurrence formulas.
J. Math. Phys. , , Crowson,
H.L. An analysis of a second order linear ordinary differential equation with five regular singular points. J. Math. Phys. , , Crowson,
H.L.
Hypergeometric solutions of a second-order linear ordinary differential equation with n-regular singular points. J. Math. Phys. , , Van
Hoeij,
M.,
Kunwar,
V.J.
Finding type solutions of differential equations with singularities. ACM Commun. Comput. Algebra , , Kunwar,
V.J.
Hypergeometric
Solutions of Linear
Differential
Equations with
Rational
Function
Coefficients.
Ph.D.
Thesis,
Florida
State
University, Tallahassee, FL, USA,
Kruglov,
V.E.
Solution of the linear differential equation of n th-order with four singular points. Ann. Univ. Sci. Bp. Sect. Comp. , , Ince,
E.L.
Ordinary Differential Equations ; Dover:
New
York,
NY, USA,
Redkov,
V.M.,
Ovsiyuk,
E.M.
Quantum Mechanics in Spaces of Constant Curvature ; Nova
Science
Publishers:
New
York,
NY,
USA,
Marin, M. Cesaro means in thermoelasticity of dipolar bodies. Acta Mech. , , Ishkhanyan,
A.M.
The third exactly solvable hypergeometric quantum-mechanical potential.
EPL , , Ishkhanyan,
A.M.
Schrödinger potentials solvable in terms of the general Heun functions.
Ann. Phys. , , Exton, H. Multiple Hypergeometric Functions and Applications ; Halsted
Press:
Chichester,
UK,
Dattoli,
G.;
Cesarano, C. On a new family of Hermite polynomials associated to parabolic cylinder functions. Appl. Math. Comput. , , Cesarano, C. Generalized special functions in the description of fractional diffusive equations. Commun. Appl. Ind. Math. , , Ishkhanyan,
A.M.
Appell hypergeometric expansions of the solutions of the general Heun equation.
Constr. Approx. , , xioms , , 102 8 of 8 Slavyanov,
S.Y.
Relations between linear equations and
Painlevé’s equations.
Constr. Approx. , , Iwasaki,
K.,
Kimura,
H.,
Shimomura,
S.,
Yoshida, M. From Gauss to Painlevé: A Modern Theory of Special Functions, Aspects of Mathematics ; Vieweg:
Braunschweig,
Germany,
16. © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).16. © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).