Generalized-hypergeometric solutions of the general Fuchsian linear ODE having five regular singularities

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.


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: which presents the most general Fuchsian equation having five regular singular points [3,4].The singularities are located at  ( 1 ) ( ) , which is the most general Fuchsian equation having four regular singular points [6][7][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 − 1 θ and 0 θ .
There are many other choices of parameters that reduce Equation (1) to the general Heun Equation (3).For instance, an obvious choice is 0 0 a ε θ = = = .This and the previous reduction cases can be made symmetric if the parameters 1 θ and 0 θ are expressed through new parameters q and 1 q such that 1 1 q q θ = + and 0 1 qq θ = .Then, for the triads ( , , ) (0, 0, 0) a q ε = and 1 1 1 ( , , ) (0, 0, 0) a q ε = Equation (1) becomes the general Heun equation with the third finite singularity being located at 1 z a = and z a = , respectively.For this reason, in this case, one may conventionally think of the parameters q and 1 q as accessory parameters associated with the singularities a and 1 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][7][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.

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] ( ) ,..., ; ,..., ; for which the coefficients obey the two-term recurrence relation ( ) The Frobenius solution of Equation ( 1) for the vicinity of the regular singularity located at where the characteristic exponent μ may adopt two values: 0 μ = or 1 γ − , and the expansion coefficients obey a four-term recurrence relation [3,4]: For the exponent 0 μ = the coefficients are given as ( ) ( ) We note that the series ( 6) with such coefficients may terminate only if α or β is a non- positive integer (besides, 0 θ and 1 θ 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 2 r N = + , 1 s N = + , and ,...., , , ,..., This means that we look for a solution of Equation ( 1) as 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 ) where . With this, the recurrence relation ( 7) is rewritten as This is a polynomial equation in summation index n of the degree utmost 4 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 ( )( )( ) The equation for the term proportional to Calculating the sum of all coefficients, we derive another relation: ( ) Finally, one more relation is derived if one checks the value adopted by the product . Notably, the product turns to be zero.In the force of Equation (18), this means that ε and 1 ε are non-positive integers such that

Explicit Solutions of Lowest Orders
The case 0 N = is almost trivial: In this case the parameters 1 θ and 0 θ are such that the singularities at z a = and 1 z a = both disappear so that Equation ( 1) is reduced to the ordinary Gauss hypergeometric equation.For 1 N = we have two choices: ε ε = − , both reproducing the result for the general Heun equation [14][15][16][17][18]12].Here is the explicit solution for the first choice: Because of the particular value of 1 θ , in this case the singularity at Equation ( 1) is reduced to the general Heun equation.Resolving the second Equation ( 22) with respect to 1 e and substituting the result into (23), we arrive at a quadratic equation for 0 θ , which provides two values of this accessory parameter for which the solution ( 21) applies.This equation reads ( ) ( ) This solution applies if where The third choice F : where Equations ( 29) and ( 30) parameterize 0 θ , 1 θ through the parameters 1 e and 2 e , while the last two equations define the latter parameters.Eliminating 1 2 , e e from Equations ( 29) and ( 30), it can be shown that 1 θ can be presented as a rational function of 0 θ : while 0 θ obeys a (rather cumbersome) fourth-order polynomial equation the coefficients of which are functions of the parameters 1 , , , , a a α γ δ : This result has been checked by Mathematica.

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][19][20][21][22][23][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 3 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 2 N = to generate a new exactly solvable potential for the onedimensional 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][31][32][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].

1 N
ε ε + = − .With these observations, Equation (16) is reduced to a polynomial equation for n of the degree N : of the parameters of the starting Equation (1) and the parameters 1 θ and 1 ,..., N e e .Resolving then the equations0 m A = , 0,1,.., m N =, we arrive at a polynomial equation of the degree 1 N + for 1 θ and the corresponding parameters 1 ,..., N e e for the generalized-hypergeometric solution(14) of Equation (1).
two cases a singularity disappears and the result reproduces that for the general Heun equation.For the first choice the solution reads[12,18] − − suggests a non-trivial solution again in terms of the functions 4 3