Next Article in Journal
Faber Polynomial Coefficient Estimates for Bi-Close-to-Convex Functions Defined by the q-Fractional Derivative
Previous Article in Journal
Fault Detection and Identification with Kernel Principal Component Analysis and Long Short-Term Memory Artificial Neural Network Combined Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Overlapping of Lévai’s and Milson’s e-Tangent-Polynomial Potentials along Symmetric Curves

AI-Solutions Inc., Lanham, MD 20706, USA
Axioms 2023, 12(6), 584; https://doi.org/10.3390/axioms12060584
Submission received: 2 May 2023 / Revised: 5 June 2023 / Accepted: 6 June 2023 / Published: 12 June 2023
(This article belongs to the Section Mathematical Physics)

Abstract

:
The paper examines common elements between Lévai’s and Milson’s potentials obtained by Liouville transformations of two rational canonical Sturm–Liouville equations (RCSLEs) with even density functions which are exactly solvable via Jacobi polynomials in a real or accordingly imaginary argument. We refer to the polynomial numerators of the given rational density function as ‘tangent polynomial’ (TP) and thereby term the aforementioned potentials as ‘e-TP’. Special attention is given to the overlap between the two potentials along symmetric curves which represent two different rational forms of the Ginocchio potential exactly quantized via Gegenbauer and Masjed-Jamei polynomials, respectively. Our analysis reveals that the actual interconnection between Lévai’s parameters for these two rational realizations of the Ginocchio potential is much more complicated than one could expect based on the striking resemblance between two quartic equations derived by Lévai for ‘averaged’ Jacobi indexes.

1. Introduction

Nearly 60 years ago, Bose [1,2] made a revolutionary discovery that the most solvable potentials known at that time can be obtained by the ‘Liouville transformation’ [3,4] of some special cases of a rational Sturm–Liouville equation [4] written in its canonical (no first derivative) form and simply referred to below as ‘RCSLE’. Several years later, the author [5] made use of Bose’s discovery to construct the family of rational potentials exactly solvable in terms of superpositions of two hypergeometric functions of a real variable as well as its confluent counterpart exactly solvable in terms of superpositions of confluent hypergeometric series and Whittaker function (see, e.g., Ch. 5 and Ch. 6 in [3]). It was shown that the corresponding eigenfunctions are expressible in terms of classical Jacobi and classical Laguerre polynomials with degree-dependent indexes in general. Thereby, we refer to these two families of solvable rational potentials as ‘Jacobi-reference’ ( J Ref) and ‘Laguerre-reference’ ( L Ref) potentials.
The Bose approach was further advanced by Milson [6], who demonstrated that there is another family of solvable potentials obtained by the Liouville transformation of the real Fuschian SLE with two poles at complex-conjugated points. It was proven by us later [7] under the influence of Alvarez-Castillo and Kirchbach’s breakthrough paper [8] that the eigenfunctions of this RCSLE are expressible in terms of Romanovski–Routh (R–Routh) polynomials with degree-dependent indexes. (For the reasons explained in [9,10], we prefer to refer to pseudo-Jacobi polynomials [11,12,13] as Routh polynomials [14], so the epithet ‘R–Routh polynomials’ is fully consistent with the term ‘Romanovski/pseudo-Jacobi polynomials’ in Leski’s [15,16] classification scheme of Romanovski polynomials [17,18]).
While making some inspirational remarks on the original draft of [19], Kirchbach drew the author’s attention to the reference to Milson’s paper [6] in the brand-new review article [20]. It became clear that Milson has already addressed this problem to a large extent. However, our study revealed some new important elements not discussed in [6]. First, it was shown that hypergeometric polynomials constructed in [6] are nothing but R–Routh polynomials. Secondly, the author took advantage of Stevenson’s idea [21] (also pointed to the author by Kirchbach) to express an analytically continued solution in terms of hypergeometric polynomials in a complex argument. It was just confirmed that the latter formally complex polynomials can be converted to real R–Routh polynomials in the real argument and that the Sturm–Liouville problem in question is indeed exactly solvable.
Our current analysis explicitly takes advantage of the interconnection [8] between the real R–Routh polynomials and Jacobi polynomials with complex-conjugated indexes, provided that the argument of the complex polynomials is restricted to the imaginary axis [22,23]. To our knowledge, it was Bagrov et al. [24] who brought the latter polynomials into the quantum mechanics using them as polynomial components of the eigenfunctions for the very unique trigonometric limit of the potentials solvable in terms of the R–Routh polynomials—the so-called ‘Rosen–Morse II’ potential in the Cooper–Khare–Sukhatme (CKS) [25,26] classification scheme of solvable rational potentials. In following [27], we refer to this trigonometric modification of the Rosen–Morse potential [28] simply as ‘t-RM’. The cited textbook [24] made no mathematical references and did not provide any arguments in support of the conjecture that the constructed eigenfunctions are real.
A few years later, Dabrowska et al. [29] ran into these polynomials while examining the eigenfunctions for the Gendenshtein potential [30] (the Scarf II potential in the CKS classification scheme [25,26]). The cited authors explicitly computed the first three eigenfunctions and did confirm that they are real and mutually orthogonal. Under the evident influence of these authors, Lévai [31], in his search for ‘shape-invariant’ [30] rational potentials, extended to complex Jacobi polynomials the systematic method suggested by Bhattacharjie and Sudarshan [32] for constructing potentials exactly solvable via classical Jacobi polynomials. He found that the list of shape-invariant potentials compiled in [29] missed the aforementioned t-RM potential.
Though the fact that the eigenvalues of the J Ref CSLE coincide with one of real roots of a quartic equation was originally recognized by Grosche [33], it was Lévai [34] who proved that the average of the Jacobi indexes of the polynomial forming a ‘quasi-rational’ [35] eigenfunction is determined by a root of a similar quartic equation in the particular case of an even density function. It was originally overlooked by the author [7,10] that Lévai (apparently unaware of Milson’s paper [6]) implicitly discussed both J Ref and Routh-reference (‘ R Ref’) CSLEs with even density functions. (We use the epithet ‘ R Ref’ to stress that the CSLE in question have ‘quasi-rational’ [35] solutions (q-RSs) expressible in terms of Routh polynomials.) In particular, the quartic equation derived by us [7] for complex-conjugated indexes of the Jacobi polynomials forming eigenfunctions of the R Ref CSLE happened to be another representation for the equation obtained by Lévai [34] in a more general context. In following the terminology suggested in our papers [7,36,37], we refer to the numerator of the rational density function in the RCSLEs of our interest as the ‘tangent polynomial’ (TP), so we term the Liouville potential associated with the even density function as the ‘e-TP’ potential for briefness.
In [38], Lévai presented a more thorough analysis comparing his approach [31,34,39] and our original technique [5] used for constructing the general 1D Schrödinger equation exactly solvable via a superposition of hypergeometric functions. An examination of his arguments illuminated in the next section revealed that, while starting from the complex Jacobi equation, Lévai implicitly converts it to the CSLE initially introduced by us [5] under the real field, with the only (though very important for his analysis) difference being that the new CSLE was expressed in terms of the variable 2z(x) − 1 instead of the variable z(x) used in [5] to convert the 1D Schrödinger equation to the hypergeometric equation in z. Obviously, the conversion of the CSLE by a linear transformation of the variable may not affect the results, so there is no surprise that introducing the solved-by-polynomials equation [31] and then converting it to its canonical form leads to exactly the same results as the Bose technique [1,2] directly starting from the resultant CSLE. Contrary to Lévai’s assertion in [40], it is the new variable (not the difference in the approach) ‘suits better the formulation of solvable P T -symmetric potentials than that in [5]’. The author has already taken advantage of this variable introduced in Lévai’s earlier works [31,39] while discussing the exactly solvable J Ref and R Ref problems in parallel in [19]. Regrettably, this re-examination of our original study [5] (as well as its extension by Milson [6]) was disregarded in Lévai’s papers [38,40,41] (see Section 2 for more details).
Though the primary focus of this paper is the real rational Sturm–Liouville problems solvable by polynomials, we (under the influence of Lévai’s cited works [31,39]) start our analysis in Section 2 from a complex (non-self-adjoint) RCSLE such that its two real-field reductions result in the self-adjoint RCSLEs with J Ref and R Ref Liouville potentials. It is worth stressing in this context that our persistent references to the Schrödinger equation with solvable rational potentials are made solely because the RCSLEs represented in their Liouville form have broad applications in quantum mechanics and are much better known to scientists. Also, note that Everitt’s [4] catalogue of Sturm–Liouville differential equations includes some renowned examples of the solvable J Ref and L Ref Liouville potentials, while disregarding their generalizations presented in our paper [5] as well as the extension of our technique to the Fuschian RCSLE with two poles on the imaginary axis [6].
In Section 3, we present a unified approach to the Liouville potentials of the J Ref and R Ref CSLEs with even density functions by treating them as two real branches of the generally complex potential of the ‘Lévai class’. We term these two branches as Lévai’s and Milson’s e-TP potentials. The new element of our analysis of Lévai’s e-TP potential, compared with [34], is the proof [36] that the (n + 1)-th eigenfunction (specified by the label c in our notation) for ‘nearly symmetric’ potentials is accompanied by a triplet of q-RSs composed of Jacobi polynomials of the same degree n which belong to three different types a, b, and d. Since any solution vanishing at one of the quantization ends (types a and b) and lying below the ground-energy level are necessarily nodeless [42], the q-RSs of these two types can be used as the factorization functions (FFs) for the so-called [43] ‘rational Darboux transformations’ (RDTs) giving rise to new exactly solvable rational potentials [37]. We also point to some important details in both Milson’s [6] and our [7,10] papers which are absent in Lévai’s sketch of this little-known problem of exactly solvable in terms of R–Routh polynomials with degree-dependent indexes.
As has been already asserted by Lévai [34], the symmetric reduction of the R Ref potential is nothing but an alternative representation for the Ginocchio potential on the line [44]. Though the assertion itself happened to be correct [45], the interrelationship between the two alternative parametrizations of the Ginocchio potential examined in Section 4 turned out to be much more complicated than one would expect based on Lévai’s observation that the quartic equations for the averaged Jacobi indexes in both quantization schemes may be rewritten in the unified fashion.
The most important consequence of the presented proof is that, in addition to the renowned quantization scheme by classical Gegenbauer polynomials [44], the cited symmetric potential can be alternatively quantized by R–Routh polynomials of a definite parity termed by us [7,45] as ‘Masjed-Jamei polynomials’ to give credit to the scrupulous analysis of these polynomials in [46]. In other words, the RCSLEs with even J Ref [44] and even R Ref [7] Bose invariants constitute the same spectral problem unambiguously defined by the common Liouville form of its two rational realizations—the main result of this study proven in Section 4.

2. Quartic Equation for the Average of Indexes of Jacobi Polynomial Forming a q-RS of Complex Fuschian CSLE with Three Singular Points

Let us start our analysis from the complex (non-self-adjoint) Fuschian CSLE with three singular points
    d 2     d η 2 + I [ η ; Λ ¯ ; T K ; ε ]   Φ [ η ; Λ ¯ ; T K ; ε ] =   0 ,
representing its Bose invariant [1,2,6] as
I [ η ; Λ ¯ ; T K ; ε ] = I o [ η ; Λ ¯ ] + ρ   [ η ; T K ]   ε ,
with both reference polynomial fraction (RefPF)
I o [ η ; Λ ¯ ] =   h o ; 4 ( η + 1 ) 2   h o ; + 4 ( 1 η ) 2   O 0 o 4 ( 1 η 2 )
and the density function
ρ   [ η ; T K ] = T K [ η ] ( 1 η 2 ) 2                     ( K 2 )
having second-order poles at η = ±1. Note that, in contrast with J RefPF (3), Bose invariant (2) also depends on the coefficients of the TP of degree K ≤ 2:
T K [ η ; a , c ± ] = 1 2 [ c + ( 1 + η ) + c ( 1 η ) ] a ( 1 η 2 )
= 1 4 [ c + ( 1 + η ) 2 + c ( 1 η ) 2 d ( 1 η 2 ) ] ,
where
d 4 a c + c
The real-field self-adjoint reduction of J Ref CSLE (1) was initially introduced by us in [19] to treat J Ref and R Ref CSLEs in parallel, to a large extent under the influence of Lévai’s renowned papers [31,39]. Examination of (4) in [38] shows that
I [ η ; Λ ¯ ; T K ; ε ] = R [ η ; Λ ¯ ; T K ; ε ] 1 2 Q [ η ; Λ ¯ ; T K ; ε ] 1 4 Q [ η ; Λ ¯ ; T K ; ε ] ,
while Lévai’s second-order differential Equation (2) can be rewritten as
( 1 η 2 ) d 2       d η 2 + τ [ η ; Λ ¯ ; T K ; ε ] d     d η + ε   F [ η ; Λ ¯ ; T K ; ε ] = 0
with
τ [ η ; Λ ¯ ; T K ; ε ] ( 1 η 2 ) Q [ η ; Λ ¯ ; T K ; ε ]
and
( 1 η 2 ) R [ η ; Λ ¯ ; T K ; ε ] = ε .
Note that the parameters p I ,   p I I , and p I I I in (9) in [38] are related to the coefficients a, c ± of TP (5) in the trivial fashion:
p I a ,     p I I 1 2 ( c + + c ) ,     p I I I 1 2 ( c + c ) .
Lévai’s starting Formula (13) in [38], with
ϕ [ η ] T K [ η ] = ( 1 η 2 ) 2 ρ   [ η ; a ,   c ± ]              
in our notation, is nothing but the conventional representation of the given Liouville potential rewritten as
V [ η ; h o ; ± , O 0 o ;   a ,   c ± ] + 1 2 { η , x } = ρ 1   [ η ; a ,   c ± ]     I o [ η ; h o ; ± , O 0 o ;   a ,   c ± ] ,
where the Schwarzian derivative { η , x } is expressed in terms of η (cf. (2.6) in [19]), provided that the change of variable η ( x ; T K ) satisfies the ordinary differential equation (ODE)
d x d η = 1 η 2 T K [ η ] ,
with prime denoting the derivative with respect to x. The parameters s I ,     s I I , and s I I I in the mentioned formula (or similarly in (1.29) in [41]) are thus related to the parameters Λ ¯ of J RefPF (3) above as follows:
s I 1 4   ( O 0 o h o ; h o ; + ) , s I I 1 2 ( h o ; + h o ; ) , s I I I 1 2 ( h o ; + h o ; + ) .
Disregarding the enhancements suggested by us in [19], Lévai simply compared his approach with our initial scheme [5] using the variable z(x) satisfying the ODE
d x d z = 2 z ( 1 z ) T 1 K [ z ] .
This transformation converts the Schrödinger equation into the J Ref CSLE with the Bose invariant
I 1 [ z ; Λ ¯ ; T 1 K ; ε ] 4   I [ 2 z 1 ; Λ ¯ ; T K ; ε ] ,
where T 1 K are the coefficients of the TP
T 1 K [ z ] T K [ 2 z 1 ] .
Though our initial technique [5] dealing with hypergeometric functions of the variable z ∈ (0, 1) is necessary to prove that the given potential is indeed exactly solvable (as well as to derive close-form expressions for the scattering amplitudes [47,48]), the use of the variable η = 2z − 1 ∈ (−1, 1) allows one to treat J Ref and R Ref potentials in a symmetric fashion and also makes it easier to examine P T -symmetric reductions of the complex J RefPF potential as it has been done by Lévai [34,38,40,41]. Obviously, the conversion of the CSLE by a linear transformation of the variable may not affect the results, so there is no surprise that starting from the solved-by-polynomials equation [31] and then converting it to the canonical form leads to exactly the same results as the Bose technique [1,2] directly starting from the resultant CSLE. Lévai’s assertion in [40] that the discussion of this problem in [38] ‘revealed that his approach suits better the formulation of solvable P T -symmetric potentials than that in [5]’ is not precisely accurate—it is the change of variable
η ( x ; T K ) = 2 z ( x ; T 1 K ) 1
that suits better his analysis [38] compared with the variable z ( x ; T K ) . Namely, if the reflection η → − η keeps the TP unchanged—the case of our current interest—then one can choose η ( x ; T K ) to be an odd function of x. If all the TP coefficients are required to be real, then the Schwarzian derivative { η , x } becomes an even real function of x and the P T -transformation is equivalent to the complex conjugation of RefPF (3) followed by the reflection of the argument η.
As another novel development inspired by Lévai’s works [31,34,38,39], the author can point to our recent idea [9,49,50] to introduce complex (non-self-adjoint) CSLE (1) and then examine its so-called [36,51,52] ‘almost-everywhere holomorphic’ (AEH) solutions
  ϕ k , m [ η ; Λ ¯ ; T K ] =       =   ± ( 1 + η )     ρ ; k , m Π m [ η ; η ¯ ( m ) ]
which exist at some energies ε k , m inside the vertical band |Re η| < 1 in the complex plane. It will be proven below that the monomial product
Π m [ η ; η ¯ ( m ) ] l = 1 m ( η η l ( m ) )
coincides with the monic Jacobi polynomial with generally complex indexes
Π m [ η ; η ¯ ( m ) ] = P ^ m ( λ + ; k , m , λ ; k , m ) ( η ) .
(It would be more accurate to use the notation ρ ± ; k , m ( Λ ¯ ; T K ) , λ ± ; k , m ( Λ ¯ ; T K ) , and ε k , m ( Λ ¯ ; T K ) , but we disregard the dependence of the latter quantities on the parameters Λ ¯ and T K for briefness.) Note that we use the term ‘AEH solutions’, instead of the equivalent epithet ‘q-RSs’ appearing in the section title, to stress that we deal with complex functions analytically continued from the real axis into the complex plane.
We then took advantage of the fact that the so-called [53,54] ‘differential polynomial system’ (DPS) composed of complex Jacobi polynomials [55] allows the second real-field reduction formed by Routh polynomials [14], in addition to the one formed by conventional (real) Jacobi polynomials.
To derive the necessary and sufficient conditions for CSLE (1) to have an AEH solution (21), first note that characteristic exponents (ChExps) of these solutions for the poles of CSLE (1) at ±1 satisfy the indicial equations
( 2 ρ ± ; k , m 1 ) 2 = h o ; ± + 1 c ± ε k , m .
Introducing the complex exponent differences (ExpDiffs)
λ ± ; k , m 2 ρ ± ; k , m 1 ,
we come to the equation
λ ± ; k , m 2 = h o ; ± + 1 c ± ε k , m .
Representing RefPF (3) as
I o [ η ; h o ; ± , μ o ] =   h o ; + h o ; + + ( h o ; + h o ; ) η 2 ( 1 η 2 ) 2   + μ o 2 1 4 ( 1 η 2 ) ,
where the parameter
μ o h o ; + + h o ; O 0 o + 1
is chosen to coincide with the ExpDiff for the pole of the J Ref CSLE at infinity, and examining asymptotic behavior of AEH solutions near this pole, we can write the closing equation for the given system of algebraic equations in λ ± ; k , m and ε k , m as follows:
a   ε k , m = 1 4   μ o 2 ρ ; k , m ( ρ ; k , m + 1 ) ,
where
ρ ; k , m = ω k , m m 1
is the ChExp for the pole of CSLE (1) at infinity and
ω k , m 1 2 ( λ ; k , m + λ + ; k , m ) .
It can be shown that the derived system of algebraic Equations (26), (29)–(31) is simply another representation for coupled Equations (14)–(16) in [34].
Squaring (31) gives
2 λ + ; k , m   λ ; k , m = 4 ω k , m 2 λ + ; k , m 2 λ ; k , m 2 ,
while substituting (30) into (29) brings us to the following energy dispersion formula:
a   ε k , m = 1 4   ( μ o 2 + 1 ) ( ω k , m + m + 1 2 ) 2 ( K = 2 ,     a 0 ) .
Rewriting (26) as
λ ± ; k , m 2 = h ˜ o ; ± + 1 + c ± ( ω k , m + m + 1 2 ) 2 / a ,
where
h ˜ o ; ± h o ; ± 1 4 c ± μ o 2 / a ,
squaring (32) and making use of (34), we come to the following quartic equation
1 4 a 2 [ 4 ω k , m 2 h ˜ o ; + h ˜ o ; 2 ( c + + c ) ( ω k , m + m + 1 2 ) 2 / a ] 2 [ a (   h ˜ o ; + + 1 ) c + ( ω k , m + m + 1 2 ) 2 ] × [ a   ( h ˜ o ; + 1 ) c ( ω k , m + m + 1 2 ) 2 ] = 0
in ω k , m . Note that the leading coefficient of quartic Equation (36) coincides with the TP discriminant
Δ T ( a , c ± ) = 1 4 ( c + c ) 2 2 a ( c + + c ) + 4 a 2 .
For each of four generally complex roots of this equation (or for each of three roots if the TP discriminant vanishes [56]), the corresponding pair of the ExpDiffs λ ± ; k , m coincides with two roots of the quadratic equation
λ ± ; k , m 2 2 ω k , m λ ± ; k , m + 1 2 ω k , m 2 1 2 ( h ˜ o ; + h ˜ o ; + ) 1 1 2 ( c + c + )   ( ω k , m + m + 1 2 ) 2 / a = 0 .
Expressing again λ ± ; k , m 2 in terms of ω k , m 2 via (34), one finds
  λ ; k , m = h ˜ o ; h ˜ o ; ± ( c ± c )   ( ω k , m + m + 1 2 ) 2 / a 2 ω k , m ω k , m .
Keeping in mind that AEH solutions (21) obey the ODE
  ϕ k , m [ η ; h o ; ± , μ o ; T 2 ] +   I   [   η ; h o ; ± ; T 2 ; ε k , m ]   ϕ k , m [ η ; h o ; ± , μ o ; T 2 ] =   0
with the Bose invariant
I [ η ; h o ; ± , μ o ; T 2 ; ε k , m ] =   h o ; c   ε k , m 4 ( η + 1 ) 2   h o ; + c +   ε k , m 4 ( 1 η ) 2       + h o ; + + h o ; + 1 μ o 2 + d   ε k , m 4 ( η 2 1 )
and substituting (21) and (23) into (40) then brings us to the ODE
( 1 η 2 ) P   m ( λ + ; k , m , λ ; k , m ) ( η ) 2 P 1 ( λ + ; k , m , λ ; k , m ) ( η ) P   m ( λ + ; k , m , λ ; k , m ) ( η ) + m ( m + λ + ; k , m + λ ; k , m + 1 ) P m ( λ + ; k , m , λ ; k , m ) ( η ) = 0
if we set
2 ρ + ; k , m   ρ ; k , m + 1 4 ( h o ; + + h o ; μ o 2 + 1 + d ε k , m ) = m ( 2 ω k , m + m + 1 )
One can easily verify the latter relation by substituting (7), (25), (26), and (33) into the left-hand side of (43).
We thus proved that q-RSs (21) have the form
  ϕ k , m [ η ; Λ ¯ ; a , c ± ]       =   ± ( 1 + η )     ρ ; k , m P m ( λ + ; k , m , λ ; k , m ) ( η )
and thereby arrived at the starting point of Lévai’s approach [34,38,40,41] while moving in the opposite direction.
Setting
δ 1 2 ( c + + c ) / a ,
Σ 1 2 ( h ˜ o ; + + h ˜ o ; + 2 ) δ 1 4 ,
and also taking into account that, as a direct consequence of (34) and (45),
1 2 ( λ + ; k , m 2 + λ ; k , m 2 ) = 1 2 ( h ˜ o ; + + h ˜ o ; + 2 ) δ ( ω k , m + m + 1 2 ) 2 ,
we come to (10) in [34] with the Jacobi indexes α and β dependent on the polynomial degree; namely,
δ ( ω k , m + m + 1 2 ) 2 δ + 1 2 ( λ + ; k , m 2 + λ ; k , m 2 ) 1 4 = Σ
in our notation.
The simplification utilized in Lévai’s aforementioned papers takes place in the particular case:
c + = c = δ / C
when the parameter
Λ k , m 1 4 ( λ ; k , m 2 λ + ; k , m 2 )
becomes independent of the polynomial degree:
Λ k , m Λ = 1 4 ( h o ; h o ; + ) .
Subtracting one of Equation (34) from another, we come to the crucial relation [34]
λ ; k , m 2 λ + ; k , m 2 = h o ; h o ; + = 4 Λ
which allows one to combine two separate quartic equations in λ ± ; k , m into a single quartic equation in ω k , m . Namely, making use of (31), one finds
λ ; k , m λ + ; k , m = 2 Λ / ω k , m ,
so [34]
λ ± ; k , n = ω k , n   Λ / ω k , n
and
1 2 ( λ + ; k , m 2 + λ ; k , m 2 ) = ω k , m 2 + Λ 2 / ω k , m 2 .
We can thus rewrite (47) as
ω k , m 2 + Λ 2 / ω k , m 2 = 1 2 ( h ˜ o ; + + h ˜ o ; + 2 ) δ ( ω k , m + m + 1 2 ) 2 .
Substituting (46) into the right-hand side of this equation gives
δ ( ω k , m + m + 1 2 ) 2 + ω k , m 2 Σ δ 1 4 + Λ 2 / ω k , m 2 = 0 ,
which brings us to quartic Equation (13) in [34]:
( 1 + δ ) ω k , m 4   + δ ( 2 m + 1 ) ω k , m 3 + [ δ ( m + 1 2 ) 2 δ Ξ 1 4 ] ω k , m 2 + Λ 2 = 0
The trivial real-field reduction of CSLE (1) obtained by choosing all six parameters Λ ¯ , T K to be real leads to the very special representative of the family of the J Ref potentials [5,36] referred to by us as ‘Lévai’s e-TP potential’. However, as it has been noticed by Lévai himself [34], there is another family of real CSLEs discovered by Milson [6] and thereby referred to by us [45] as ‘Milson’s e-TP potential’. It is obtained by requiring the parameters h o ; + and h o ; to be complex-conjugated while keeping real the parameter O 0 o . Both RefPF and density function thus become real when expressed in terms of the new variable η I = i η . Lévai’s formalism outlined above made it possible to treat both potentials in the unified matter outlined in Section 3.

3. Two Real-Field Reductions of Complex J Ref SLE

The important common feature of the J Ref and R Ref potentials on the line is that the corresponding density functions
ρ   [ η ; T 2 ] = T 2 [ η ] ( 1 η 2 ) 2                 ( T 2 [ ± 1 ] 0 )
and
ρ i   [ η I ; T i 2 ] = T i 2 [ η I ] ( 1 + η I 2 ) 2        
have second-order poles at the endpoints ±1 and ±∞ accordingly. As a direct consequence of this observation, the algebraic Schrödinger equation
  d     d ξ ρ ι   1 2 [ ξ ; T ι 2 ] d       d ξ V ι   [ ξ ; Λ ¯ ι ; T ι 2 ] ε ρ ι   1 2 [ ξ ; T ι 2 ]   Ψ ι [ ξ ; Λ ¯ ι ; T ι 2 ; ε ]  
(with ξ η or η I for ι = ‘ ’ or i, respectively) has the so-called [37] ‘prime’ form; namely, the sum of two ChExps equal to 0, and as a result, the principal Frobenius solution is unambiguously determined by the Dirichlet boundary conditions (DBCs). In [37], we have proved the cited below theorem which allows one to formulate the given spectral problem as the Dirichlet problem for the prime SLE.
Theorem 1. 
If the density function  ρ ι   [ ξ ; T ι 2 ]  of the given RCSLE has second-order poles at both endpoints, then the corresponding Dirichlet problem formulated for the prime SLE is fully equivalent to the requirement that the eigenfunctions of RCSLE are square-integrable with the weight  ρ ι   [ ξ ; T ι 2 ] .
One of the advantages of formulating the given spectral problem as the Dirichlet problem for the prime SLE is that we can automatically adopt the rigorous theorems proven in [42] for SLEs solved under the DBCs. As explained in [37], the requirement for the density function to have second-order poles at the endpoints automatically assures that the corresponding Liouville transformation leads to the 1D Schrödinger equation on the line. Again Gesztesy, Simon, and Teschl’s meticulous proofs made it possible to extend the conventional results of the regular Sturm–Liouville theory [57] to singular Sturm–Liouville problems including the 1D Schrödinger equation on the line.
It is worth stressing in this context that the theorem does not hold if the density function has a simple pole at one of the endpoints and thereby is invalid for the radial Schrödinger equation [37,58].
In this paper, we focus solely on the J Ref and R Ref CSLEs with the numerators of density functions (59) and (60) formed by the even second-degree polynomials
T e [ η ; δ , C ] = C 1 ( 1 + δ η 2 ) ( δ / C > 0 )
and
T i e [ η I ; a i , κ + ] = a i ( η I 2 + κ + )   ( a i ,     κ + 1 + δ + >   0 ) .

3.1. Lévai’s Real e-TP Potential

As initially noticed in [59], the variable η ; s ( x ; δ , C ) obtained by solving the ODE
η ; s ( x ; δ , C )   =   ρ   L 1 2   [ η ; s ( x ; δ , C ) ; δ , C ] ,
under the boundary condition
η ; s ( 0 ; δ , C ) = 0
with
  ρ L   [ η ; δ , C ] 1 + δ η 2 C ( 1 η 2 ) 2 ,
coincides with variable (60) in [44] if one sets
δ = ( λ 2 1 ) 1 , C = δ λ 4
in Ginocchio’s notation. Namely,
η ; s [ y ; λ ] = λ   y Y [ y 2 ; λ 2 ] ,
where the function y(x ; λ ) is a solution of the ODE [60,61]
y ( x , λ )     = [ 1 y 2 ( x , λ )   ]   Y [ y 2 ( x , λ )   ; λ 2 ]   for 1   <   y   < + 1
and
Y [ y 2 ; λ 2 ] 1 + ( λ 2 1 )   y 2 .
To prove Lévai’s assertion [59]. one can follow Wu’s pioneering arguments [60,61] in support of the claim that the radial Ginocchio potential [62] is nothing but another representation of the implicit radial potential exactly solvable in terms of hypergeometric functions [47,48]. First multiplying the square of (68) by its reverse
y 2 = η ; s 2 N [ η ; s 2 ; λ 2 ] ,
where
N [ η 2 ; λ 2 ] λ 2 + ( 1 λ 2 ) η 2
gives
N [ η 2 ; λ 2 ] = λ 2   Y [ y 2 ; λ 2 ] .
Setting
η ; s ( x ; λ ) η ; s [ y ( x , λ ) ; λ ] = η ; s ( x ; δ , δ λ 4 )
with δ = ( λ 2 1 ) 1 and also taking into account that
N [ η 2 ; 1 + 1 / δ ] = ( 1 + δ η 2 ) / δ ,
we can rewrite (64) as
η ; s ( x ; λ )   = λ 2   [ 1 η ; s 2 ( x ; λ ) ] N [ η ; s 2 ( x ; λ ) ; λ 2 ]
or, which is equivalent,
η ; s ( x ; λ )   = [ 1   y 2 ( x ; λ ) ]   N [ η ; s 2 ( x ; λ ) ; λ 2 ] .
On the other hand, we come to the same relation multiplying (70) by the derivative
d η ; s d y           = λ 2 N 3 2 [ η ; s 2 ; λ 2 ]
and making use of (73).
Algebraic Schrödinger Equation (61) thus takes the form
  d     d η ρ   L 1 2   [ η ; δ , C ] d d η ρ   L 1 2   [ η ; δ , C ] V L   [ η ; Λ ¯ ; δ , C ] ε   ×   Ψ L [ η ; Λ ¯ ; δ , C ; ε ] = 0 ,
where V L   [ η ; s ; Λ ¯ ; C , δ ] is the Liouville potential converted back to the variable η ; s ( x ; δ , C ) , i.e.,
V L [ η ; s ( x ; λ ) ; Λ ¯ ; δ , δ λ 4 ] =   [ η ; s ( x ; λ ) ] 2 I o [ η ; s ( x ; λ ) ; Λ ¯ ] 1 2 { η ; s ( x ; λ ) , x } ,
where again δ = ( λ 2 1 ) 1 . It has been shown in [45] that the Schwarzian derivative expressed in terms of the variable η ; s can be written as
{ η ; s , x } = 2 C   1 + δ η ; s 2 + C   ( 1 η ; s 2 ) 2 2 ( 1 + δ η ; s 2 ) 2   1 4 1 η ; s 2 + 5 ( δ + 1 )   1 + δ η ; s 2 .
Substituting (3), (64), (66), and (81) into potential function (80) expressed in terms of the variable η ; s , one finds
V L [ η ; Λ ¯ ; δ , C ] / C = h o ; ( 1 η ) 2 + h o ; +   ( 1 + η ) 2       + O 0 o   ( 1 η 2 ) + 4 4 ( 1 + δ η 2 ) + 1 η 2 ( 1 + δ η 2 ) 2   ( 1 η 2 ) 2 4 ( 1 + δ η 2 ) 2   1 + 5 ( δ + 1 )   1 + δ η 2
or, alternatively,
V L [ η ; Λ ¯ ; δ , C ] / C = 1 4 ( μ o 2 2 ) + Ξ 1 + δ η 2 + 3 δ ( 3 δ + 2 ) 4 ( 1 + δ η 2 ) 2   5 δ 2 ( δ + 1 ) 4 ( 1 + δ η 2 ) 3   + ( h o ; + h o ; ) η 2 ( 1 + δ η 2 ) ,
where
Σ 1 2 ( h o ; + + h o ; + 2 ) + 1 4 [ δ ( μ o 2 4 ) 1 ]
is simply another representation for parameter (46) taking into account (28) and (35).
It is essential that solutions of prime SLE (79) are related to solutions of the real J Ref CSLE with density function (66) via the simple gauge transformation
Ψ L [ η ; Λ ¯ ; δ , C ; ε ]     ρ   L 1 4 [ η ; δ , C ]   Φ L [ η ; Λ ¯ ; δ , C ; ε ]   ,
and therefore the eigenfunctions of the prime SLE and the eigenfunctions of the J Ref CSLE,
  ϕ c n [ η ; Ξ , Λ ; δ , C ]   =       =   ± ( 1 + η )   1 2 λ ; c n P n ( λ + ; c n , λ ; c n ) ( η ) ,
are interrelated via the elementary formula
  ψ c n [ η ; Ξ , Λ ; δ , C ]     ρ   L 1 4 [ η ; C , δ ]     ϕ c n [ η ; Ξ , Λ ; δ , C ] .
Examination of the asymptotic behavior of eigenfunctions (86) near poles of the J Ref CSLE at ±1 reveals that solutions (87) obey the DBCs
l i m η ± 1   ψ c n [ η ; Ξ , Λ ; δ , C ; ε ] = 0
if
λ ± ; c n > 0 .
One can then directly verify that constraints (89) are the necessary and sufficient conditions for q-RSs (86) to be square-integrable with nonnegative weight (66):
1 + 1 ϕ c n 2 [ η ; Ξ , Λ ; δ , C ] ρ     L [ η ; δ , C ]   d η < ,
as prescribed by Theorem 1. It directly follows from definition (31) of the sought-for roots, coupled with constraints (89), that the root ω c n associated with the (n + 1)-th eigenfunction must be positive.
Examination of asymptotic behavior of potential (82) near the singular endpoints shows that
V L [ ± 1 ; Λ ¯ ; δ , C ] / C = ( h o ; ± + 1 ) / δ .
We choose the energy reference point from the requirement that the potential vanishes at η = +1 setting h o ; + = 1 . We also require that the potential takes a nonnegative value at the lower endpoint by choosing
h o ; + 1 λ o 2 4 Λ 0 ,
where the parameter λ o coincides with the zero-energy ExpDiff for the pole of J Ref CSLE (1) at η = −1. Taking into account (28), we can thus represent J Ref PF (3) as
I o [ η ; λ o , μ o ] =   1 λ o 2 4 ( η + 1 ) 2   + 1 4 ( 1 η ) 2   + μ o 2 + 1 λ o 2 4 ( 1 η 2 ) .
With such a choice of the energy reference point, two positive roots λ ± ; c n of quadratic Equation (38) specifying the (n + 1)-th eigenfunction satisfy the inequality
λ ; c n λ + ; c n > 0 ,
so the sought-for roots take the form
λ ± ; c n = ω c n Λ / ω c n .
Comparing (95) with (12) in [34], we conclude that Lévai’s definition of the e-TP potential differs from ours by a reflection of the argument, i.e., z = η ; s and α = λ ; c n , β = λ + ; c n in our terms.
In the limiting case of the Ginocchio potential on the line ( λ o = Λ = 0 ) , quartic Equation (58) takes the form:
ω n 2 [ δ ( ω n + n + 1 2 ) 2 + ω n 2 Σ δ 1 4 ] = 0
with the double-zero root
ω n = 0
and two other roots determined by the quadratic equation
ω + n 2   + δ ( ω + n + n + 1 2 ) 2 δ Ξ 1 4 = 0 .
The average of the indexes of the Jacobi polynomials in the right-hand side of (86) thus coincides with the positive root ω c n of quadratic Equation (98). Since the TP may not have zeros between −1 and +1, the parameter δ must be either positive or smaller than −1, i.e., the linear coefficient of the quadratic equation
ω + n 2   + δ δ + 1 ( n + 1 2 ) ω + n + [ δ ( n + 1 2 ) 2 δ Ξ 1 4 ] / ( δ + 1 ) = 0
is necessarily positive, and therefore the quadratic equation has a positive root if its free term is negative, i.e., if
( n + 1 2 ) 2 < ( Ξ + 1 4 ) / δ + 1
in agreement with (23) in [34]. Setting Λ = 0 in (95) confirms that
λ + ; c n ( Ξ ; δ ) =   λ ; c n ( Ξ ; δ ) = ω c n ,
as expected. The negative root ω d n of quadratic Equation (99) specifies the q-RS composed of the Jacobi polynomial with two equal negative indexes:
λ + ; d n ( Ξ ; δ ) =   λ ; d n ( Ξ ; δ ) = ω d n < 0 .
Examination of quartic Equation (58) reveals that
lim Λ 0 ( Λ / ω n ) = δ + Ξ + 1 4 δ ( n + 1 2 ) 2
and therefore
λ ± ; n ( Ξ ; δ ) lim Λ 0 λ ± ; n ( Ξ , Λ ; δ ) = δ + Ξ + 1 4 δ ( n + 1 2 ) 2 .
So,
λ + ; n ( Ξ ; δ )   λ ; n ( Ξ ; δ ) < 0
as far as condition (100) holds, and therefore both Jacobi indexes are real. We thus conclude that
=     a     o r     b .
We thus confirmed that the (n + 1)-th eigenfunction is indeed accompanied by a triplet of q-RSs an, bn, and dn, in agreement with our more general results [36] for asymptotically levelled (AL) potentials ( λ o = 0 ,     c +     c ) .
As the asymmetry parameter Λ > 0 increases, the e-TP J Ref potential has at least n + 1 bound energy levels as far as the positive root of quartic Equation (58) remains larger than Λ , and thereby both Jacobi indexes stay positive, as prescribed by (89). Indeed, the latter equation may not have a positive double root larger than Λ , because this would imply the existence of two eigenfunctions with exactly the same number of nodes inside the quantization interval. Therefore, two positive roots may not merge giving rise to a complex-conjugated pair as far as ω c n 2 > Λ . The (n + 1)-th bound energy level may disappear if Λ reaches the positive root of the quadratic equation
( δ + 2 ) Λ + ( 2 n + 1 ) Λ + ( n + 1 2 ) 2 ( Ξ + 1 4 ) / δ 1 = 0
in Λ .

3.2. Milson’s e-TP Potential

As stressed above, choosing all the parameters of RefPF (3) to be real and making the Liouville transformation on the real axis [1,2] brings us directly to the conventional J Ref potential [5]. However, as mentioned in Section 2, there is another less obvious option to formulate the real Sturm–Liouville problem utilized by Milson [6]. Namely, keeping real only O 0 o and a 2 while allowing the parameters in each pair h o ; ± and c ± to be complex and also restricting CSLE (1) to the imaginary axis η = i η I , we come to the real CSLE [6]:
  d 2     d η I 2 + I [ η I ; h o , O i 0 o ; T i K ] +   ρ i   [ η I ; T i K ] ε ˜   Φ [ η I ; h o , O i 0 o ; T i K ; ε ˜ ] =   0
with the R Ref PF
I i o [ η I ; h o , O i 0 o ] =   h o 4 ( 1 + i η I ) 2   + h o 4 ( 1 i η I ) 2   +   O i 0 o 4 ( 1 + η I 2 )
and the positive density function
ρ i M   [ η I ; T i K ] = T i K [ η I ] ( 1 + η I 2 ) 2                     ( K = 0     o r     2 )
(assuming that the TP has a positive leading coefficient and a negative discriminant Δ i T for K = 2).
The conventional Liouville transformation of CSLE (108) then leads to the real R Ref potential in our classification scheme. Though the full credit for the discovery of the latter potential should be certainly given to Milson [6], two particular cases of CSLE (108) with a constant TP (K = 0) and the second-degree TP
T i 2 [ η I ] =   a i ( 1 + η I 2 )
associated with the translationally shape-invariant (TSI) R Ref potentials (t-RM and Gendenshtein potentials, respectively) have already implicitly appeared in Lévai’s pioneering paper [31], which originally drew the author’s attention to this extraneous class of rational potentials solvable in terms of Jacobi polynomials with complex-conjugated indexes and imaginary argument.
Let us now discuss the energy spectrum of R Ref CSLE (108). Since the corresponding Liouville potential is asymptotically levelled, we can choose the energy reference point via the requirement that the potential vanishes at both quantization endpoints ±∞. This can be achieved by setting [7,10]
O i 0 o = 2 h o ; R + 1 ,
where h o ; R is the real part of the complex parameter
h o ; h o h o ; R + i h o ; I .
In other words, we make parameter (28) equal to zero so h ˜ o ; ± h o ; ± and R Ref CSLE (108) takes the form [7,10]
d 2     d η + ; s 2 + I i o   [ η + ; s ; h o ]   + ρ i M [ η + ; s ; κ + ]   ε ˜     Φ i [ η + ; s ; h o ; κ ; ε ˜ ] =   0 ,
where
    I i o [ η + ; s ; h o ; R + i   h o ; I ]     = h o ; R + h o ; I η + ; s ( η + ; s 2 + 1 ) 2 + 1 4 ( η + ; s 2 + 1 ) ,
ε ˜ ε ,     a n d
ρ   M [ η + ; s ; κ + , a + ] =     a + ( η + ; s 2 + κ + ) ( η + ; s 2 + 1 ) 2       ( a + > 0 ) .
(The density function with a constant TP, giving rise to the translationally form-invariant (TFI) CSLE of Group B [63] with the trigonometric Liouville potential [27], requires a special consideration.) In Lévai’s notation [34],
κ + 1 = δ + > 1 ,     a + = C + 1 > 0 ,
and
h o ; R δ + + 3 4 Σ + ,
while parameter (51) becomes imaginary:
Λ + 1 4 ( λ ; k , m 2 λ + ; k , m 2 ) = 1 2 h i o ; I   i .
The change of variable defined by the ODE
η + ; s ( x ; κ + , a + )   =   ρ   M 1 2   [ η + ; s ( x ; κ + , a + ) ; κ + , a + ]
solved under the boundary condition
η + ; s ( 0 ; κ + , a + ) = 0
converts CSLE (114) into the Schrödinger equation with Milson’s [6] potential
V M   ( x ; h o ; κ + , a + ) V M   [ η + ; s ( x ; κ + , a + ) ; h o ; κ + , a + ] = ( η + ; s ) 2 I i o [ η + ; s ( x ; κ + , a + ) ; h o ] 1 2 { η + ; s , x } ,
where [45]
  {   η + ; s , x } = 2   a +   ( η + ; s 2 + κ + ) + 2 ( 1 + η + ; s 2 ) a + ( η + ; s 2 + κ + ) 2 ( 1 + η + ; s 2 ) 2   2 a +   ( η + ; s 2 + κ + ) 2 1 + 5 κ + η + ; s 2 + κ + .
Setting
ω i k , n R e   λ ± ; k , n
and again taking into account that parameter (119) is independent of the polynomial degree, we come to the quartic equation
( 1 + δ + ) ω i k , n 4   + δ + ( 2 n + 1 ) ω i k , n 3 + [ δ + ( n + 1 2 ) 2 δ + Ξ + 1 4 ]   ω i k , n 2 + Λ + 2 = 0 ,
while combining (54) with (119) gives
λ ± ; k , n = ω i k , n i h o ; I 2 ω i k , n .
Though quartic Equation (125) seems identical to (58), there is a catch: the roots specifying the eigenvalues of the corresponding real CSLEs with density functions (66) and (110) have opposite signs. Namely, while ω c n is positive, its counterpart ω i c n specifying the eigenvalue of the Schrödinger equation with Milson’s e-TP potential is restricted by the negative upper bound specified by (139) below. Though Lévai’s assertion that the symmetric R Ref potential coincides with the Ginocchio potential on the line [44] is formally correct [45], its proof relies on rather sophisticated arguments thoroughly illuminated in Section 4.
The quartic equation similar to (125) was independently derived by us in [19] under the influence of Alvarez-Castillo and Kirchbach’s breakthrough study [8] on the quantization of the Gendenshtein (‘Scarf II) potential [30] in terms of Romanovski polynomials [18]. The author realized that the corresponding Schrödinger equation can be obtained by the Liouville transformation of R Ref CSLE (114) with the density function
ρ i M [ η I ; 1 ] =     1 η I 2 + 1 .
By searching for eigenfunctions of the given CSLE in the ‘quasi-rational’ [35] form
ϕ i c n [ η I ; h o ; κ + ]     ( η I 2 + 1 ) ω i c n + 1   1   i   η I 1 + i   η I   λ c n ; I   Π n [ η I ; h o ; κ + ]
and examining the ExpDiff for the pole of CSLE (114) at infinity, we found that n-degree polynomials in the right-hand side of (128) satisfy the Routh equation [14]
( η I 2 + 1 ) d 2     d η I 2 + 2 [ ( ω i c n + 1 ) η I λ c n ; I ] d         d η I n ( n + 2 ω c n 1 ) × n ( ω i c n + i   λ c n ; I ) [ η I ] = 0 ,
such that its polynomial solutions n ( λ i c n ) [ η I ] after being rewritten in the monic form coincide with the pseudo-Jacobi polynomials
^ n ( λ i c n ) [ η ] =   P n ( η ; λ c n ; I , ω i c n 1 )
defined via (9.9.1) in [11]. It was found that the eigenvalues of R Ref CSLE (114) are unambiguously determined by the real part ω i c n of the complex index λ i c n of the Routh polynomial:
ε i c n ( h o ; κ + ) = ρ i ; c n 2 < 0 ,
with
ρ i ; c n =   ω i c n n 1   2 > 0 ,
while examination of the ExpDiffs for the poles of CSLE (114) brought us to the algebraic equation
λ i c n 2 = h o ; R + 1 + i h o ; I + ( κ + 1 )   ε i c n ( h o ; κ + ) ,
similar to (26) except that the coefficients c + and c are real and coincide with each other. It then directly follows from (133) that
I m   λ i c n 2 = 2 ω i c n   λ c n ; I = h o ; I ,
bringing us back to Lévai’s formula (126).
Substituting (131) and (132) into the right-hand side of (133) and eliminating λ c n ; I from the resultant formula
ω i c n 2 λ c n ; I 2 = h o ; R + 1 + ( 1 κ + )   ( ω i c n + n + 1   2 ) 2 ,
we [19] came to the quartic equation
ω i c n 2 [ ω i c n 2 h o ; R 1 + ( κ + 1 )   ( ω i c n + n + 1   2 ) 2 ] 1   4 h o ; I 2 = 0 ,
which is nothing but an alternative form of (125) with κ + = 1 + δ + and Lévai’s parameters Ξ + and Λ + defined via (118) and (119) accordingly.
By requiring eigenfunctions (128) to be square-integrable with weight (116):
+ ϕ i c n 2 [ η I ; h o ; κ + ] ρ i   M [ η I ; κ + ]   d η I < ,
we come to the constraint
ρ i ; c n =   ω i c n n 1   2 >     0 ,
which implies that the real part of the Routh index   λ i c n must be restricted by the upper bound
ω i c n <   n 1   2 .
We define R–Routh polynomials (Romanovski/pseudo-Jacobi polynomials in Leski’s terms [15,16]) via (3.5) and (3.6) in [64]:
R n ( 2 α I , α R + 1 ) ( x ) = ( i ) n   P n ( α R + i α I , α R i α I ) ( i x ) n ( α R + i α I ) ( x )
for
n N 0 α R 1   2 .
(Remember that we use the term ‘pseudo-Jacobi polynomials’ as a synonym for Routh polynomials and do not require the parameter N 0 [12] to be a nonnegative integer, so our use of this term is similar to the definition of these polynomials in [13].)
The requirement for quartic equation
( ρ i ; c n + n + 1   2 ) 2 [ ( ρ i ; c n + n + 1   2 ) 2 h o ; R 1 + ( κ + 1 ) ρ i ; c n 2 ] 1   4 h o ; I 2 = 0
to have an odd number of negative roots [7] results in Milson’s formula for the number n o of discrete energy levels [6]
n o =     1 2 ( h o ; R + 1 + | h o + 1 |   ) + 1 2 .
Finally, let us point to the fact that each eigenfunction is accompanied by another q-RS associated with the second real root of the quartic equation. If the Routh polynomial forming the second q-RS has no real zeros, then the q-RS in question can be used as the FF for the RDT, giving rise to a new solvable rational potential [10] (as originally suggested by Quesne [64] for the Scarf II potential).
It is also worth mentioning that our argumentation in [7] was restricted solely to the real (self-adjoint) CSLE, and within this approach, we were unable to prove that the real part of sought-for complex ExpDiff λ i c n coincides with a negative real root of a quartic equation in the general case of the density function [9]
ρ i   [ η ; κ + + i   κ I ] =     η 2 + κ I   η + κ + ( η 2 + 1 ) 2
with | κ I |   > 0 , assuming that the TP has a negative discriminant
Δ T =     1 4 κ I 2 κ + < 0 .
To cover the general case [9], we (under the influence of Lévai’s works [31,34,38,39]) had to start from complex (non-self-adjoint) CSLE (1) and then consider the second real-field reduction of Bochner’s [55] complex Jacobi DPS to the DPS formed by Routh polynomials [14].

4. Ginocchio Potential as an Overlap of Lévai’s and Milson’s e-TP Potentials

The main purpose of this section is to prove that Lévai’s [34] and Milson’s [6] e-TP reductions of the J Ref and R Ref potentials defined via (80) and (122) overlap along the symmetric potential curves
V L [ η ; s ( x ; δ , C ) ; 0 , μ o ; δ , C ] = V M [ η + ; s ( x ; κ + , a + ) ; h o ; R ; κ + , a + ]
obtained by setting corresponding asymmetry parameters (92) and (119) to zero:
Λ = Λ + = 0
(assuming that the J RefPF in the right-hand side of (80) is parametrized by the parameters λ o and μ o according to (93).) As explained below, this assertion can be considered as the corollary of the following theorem:
Theorem 2. 
The RCSLEs with even J Ref and even R Ref Bose invariants constitute two rational realizations of the same Sturm–Liouville problem interrelated via an algebraic change of variable.
Proof of Theorem 2. 
Let us first prove that the solutions of ODEs (64) and (120) under the boundary conditions (65) and (121), respectively, are interrelated via the following algebraic formulas:
η + ; s = η ; s 1 η ; s 2
or
η ; s = η + ; s η + ; s 2 + 1
if we choose
1 κ + = ( κ + 1 ) 1
and
a + = a ( κ + 1 )
or, in Lévai’s notation,
δ + δ = 1
and
C = δ C + ,
so
δ ± + 1 ± κ ±
and
a ± = ± C ± 1
with δ δ ,     C C in the notation of Section 3.1.
If we, following [45], introduce two nonnegative variables
z ± η ± ; s 2
satisfying the ODEs
z ± ( x ; κ ± , a ± ) = 2 z ± ( 1 ± z ± )   a ±   ( z ± + κ ± ) ,
then algebraic relation (148) is equivalent to the linear fractional transformation
z +   = z 1 z ,
assuming that the parameters κ ± ,   a ± are interrelated as specified above.
To prove (158), let us introduce the auxiliary PF
ζ   [ z ] z 1 z
and demonstrate that the function
ζ ( x ; κ , a ) ζ   [ z ( x ; κ , a ) ]
satisfies ODE (157) for the variable z + ( x ; κ + , a + ) with
ζ ( 0 ; κ , a ) = z ± ( 0 ; κ + , a + ) = 0 .
Indeed, differentiating (160) with respect to x and making use of (157) with the lower subscript, one finds
ζ ( x ; κ , a ) =     2 [ ζ ( x ; κ , a ) + 1 ]   z ( x ; κ , a ) a   [ z ( x ; κ , a ) + κ ] ,
where we also took into account that
ζ   [ z ] + 1 = 1 1 z .
Under constraints (150) and (151), the denominator of the PF in the right-hand side of (162) can be rewritten as follows:
a [ z ( x ; κ , a ) + κ ] = a + [ 1 z ( x ; κ , a ) ] × [ ζ ( x ; κ , a ) + κ + ]
which gives
ζ ( x ; κ , a ) =     2   [ ζ ( x ; κ , a ) + 1 ]   ζ ( x ; κ , a ) a + [ ζ ( x ; κ , a ) + κ + ] .
This confirms that
z + ( x ; κ + , a + ) ζ ( x ; κ , a )
as asserted.
Our next step is to prove that the parameters μ o and h o ; R can be interrelated in such a way that
I o [ η ; s ; 0 , μ o ] = d η + ; s d η ; s 2 I i o [ η + ; s ; h o ; R ] + 1 2 { η + ; s , η ; s } ,
provided that the coefficients of e-TPs (62) and (63) obey constraints (150)–(153) and therefore, as proven above, the variables z ± are related via linear fractional transformation (158), so [65]
{ z + , z } = 0 .
This implies that interrelation Formula (167) is equivalent to the requirement that the RefPFs
I 1 o [ z ; μ o ] = d η ; s d     z   2 I o [ z ; 0 , μ o ] 1 2 { η ; s 2 , η ; s }
and
I 1 + o [ z + ; h o ; R ] d η + ; s d     z +   2 I i o [ z + ; h o ; R ]   1 2 { η + ; s 2 , η + ; s } ,
where
{ η 2 , η } = 3 2 η 2 ,
are interrelated as follows:
I 1 o [ z ; μ o ] d   z + d   z   2 I 1 + o [ z + ; h o ; R ] .
Substituting (93) with λ o = 0 and (115) with h o ; I = 0 , into the right-hand sides of (169) and (170), accordingly, we can represent the RefPFs in question as
I 1 o [ z ; μ o ] = 3 16 z 2 + 1 4 z ( 1 z ) 2 + μ o 2 1 16 z ( 1 z )
and
I 1 + o [ z + ; h o ; R ] = 3 16 z + 2 + h o ; R 4 z + ( z + + 1 ) 2 + 1 16 z + ( z + + 1 ) .
Expressing the right-hand side of (173) in terms of the variable z , we then rewrite (172) as
( 1 z ) 4 I 1 o [ z ; μ o ] = 1 16 z + 2 ( z + + 1 ) 2 3 + z + ( μ o 2 + 3 ) + 4 z + 2
and thereby confirm that the right-hand side of the latter formula turns into (174) if we choose [45]
h o ; R + 1 = 1 4   μ o 2 ,
which concludes the proof. □
It seems beneficiary to also reformulate Theorem 2 as the following proposition:
Corollary 1. 
The Liouville transformations of the RCSLEs with even J Ref and even R Ref Bose invariants on the intervals (−1, +1) and (−∞, +∞) accordingly results in the Schrödinger equation with exactly the same potential symmetric under the reflection of its argument x.
This constitutes the main result of this paper.
As illuminated in detail in [45], the Liouville transformation of the CSLE
    d 2     d z 2 + I 1 o [ z ; μ o ]     ( z δ 1 )   ε 4 C z ( 1 z ) 2   Φ 1 [ z ; μ o ; δ , C ; ε ] =   0
on the finite interval (0, 1) converts it to the Schrödinger equation with the non-singular radial J Ref potential. By reflecting the latter potential around the origin, the author [5] originally constructed the symmetric potential later rediscovered by Ginocchio [44] in the form
V G [ y ] = ( 1 y 2 ) λ 2 ν ( ν + 1 ) 1 4 ( 1 λ 2 ) [ 5 ( 1 λ 2 ) y 4 ( 7 λ 2 ) y 2 + 2 ] ,
where the variable y is related to the variable z as follows:
z = λ 2   y 2 1 + ( λ 2 1 )   y 2 .
It was Wu [60] who recognized the equivalence of two representations as soon as Ginocchio presented his results.
Later, Lévai [59] pointed to the fact that the Ginocchio potential on the line turns into the symmetric J Ref potential associated with e-TP (62) if we choose
μ o = 2 ν + 1 , δ = ( λ 2 1 ) 1 , C = δ λ 4
and then express (178) in terms of variable (68). Setting h o ; + = h o ; = 1 in (84) and also making use of constraint (49) to replace λ ± ; k , m in (26) for
λ ; c n =   λ + ; c n ω c n
thus gives
Σ + δ + 1 4 = 1 4 δ μ o 2
and
  ε c n =   λ 4   ω c n 2 ,
respectively. As expected, energy dispersion formula (183) and quadratic equation
( λ 2 1 )   ω c n 2 + (   ω c n + n + 1 2 ) 2 1 4 μ o 2 = 0
match, respectively, (3.9) and (3.10) in [44], with μ standing for ω c n in our notation. Taking into account that both leading and linear coefficients of the quadratic equation are positive, the latter may have a positive root
  ω c n = λ 2 1 4 λ 2 μ o 2 ( λ 2 1 ) ( n + 1 2 ) 2 n 1 2 > 0
iff its free term is negative, i.e., iff
0 n < 1 2 ( μ o 1 ) .
On the other hand, substituting (176) into (118) shows that
Σ + + δ + + 1 4 = 1 4   μ o 2 ,
so the quadratic polynomial in the brackets in quartic Equation (136) takes the form
ω i c n 2   + ( λ 2 1 ) ( ω i c n + n + 1 2 ) 2 1 4   μ o 2 = 0 ( ω i c n < n 1   2 ) .
Setting side by side quadratic Equations (184) and (188), we conclude that the positive root of the former equation specifies ChExp (132) of the (n + 1)-th eigenfunction for the pole of CSLE (114) at infinity:
ρ i ; c n = ω c n =   ω i c n n 1   2 > 0 ,  
provided that
0 n < h o + 1 1 2 ,
as prescribed by (143) with h o = h o ; R .
Finally, comparing (187) with (182), we find that Lévai’s parameters (84) and (118) are related via a nontrivial formula
Σ + = δ + Σ 3 4 + 3 ,
which confirms our assertion that, despite the formal similarity between quartic Equations (58) and (125), the two equations are related in a rather complicated fashion in the symmetric limit Λ = Λ ± = 0 .
As already pointed to in the Introduction, the fact that the symmetric reduction of the R Ref potential leads to a certain subclass of the J Ref potentials has been already recognized by Milson [6], who cited in this connection the quadratic transformation of the hypergeometric function using the substitution
1 2 ( 1 i   η + ; s ) η + ; s 2 + 1 ,  
which is reminiscent of (156) for the variable z + . However, the author was unable to figure out all the details necessary for the suggested (possibly alternative) representation of the Schrödinger equation with the given symmetric potential.
As originally discovered by Ginocchio [44], the eigenfunctions of the Schrödinger equation with potential (178) can be expressed in terms of classical Gegenbauer polynomials with degree-dependent indexes after being converted to the variable η , s :
ϕ c n [ η , s ; μ o ; λ ] =   ( 1 η ; s 2 ) 1 2 ( ω c n + 1 )   C n ( ω c n + 1 2 ) ( η , s ) .
On the other hand, as demonstrated in [7,45], the eigenfunctions of R Ref CSLE (114) with real h o are expressible in terms of orthogonal Masjed-Jamei polynomials [46]
ϕ i c n [ η + ; s ; μ o ; λ ]       ( η + ; s 2 + 1 ) 1 2 ω i c n + 1 2 I n ( p n ) ( η + ; s )
with degree-dependent indexes p n = 1 2 ω i c n larger than n + 1.

5. Example: Sech-Squared Potential

Lévai’s [34] and Milson’s [6] e-TP potentials for a = 0       ( a κ = 1 ) and κ + =     a + = 1 turn into the Rosen–Morse [28] and, respectively, Gendenshtein [30] TSI potentials:
V R M [ t a n h ( 2 x ) ; λ o , μ o ] = 1 4 ( μ o 2 1 ) s e c h 2 ( 2 x ) + 1 2 λ o 2 [ 1 t a n h   ( 2 x ) ]
and
V G   [ s i n h   ( 2 x ) ; h o ; R + i   h o ; I ] =   s e c h 2 ( 2 x ) [ h o ; R + 3 4 + h o ; I   s i n h   ( 2 x ) ] ,
where
V R M [ η , s ; λ o , μ o ] = 1 2 ( 1 η , s ) [ 1 2 ( μ o 2 1 ) ( 1 + η , s ) λ o 2 ]
and
V G   [ η + , s ; h o ; R + i   h o ; I ] =   4   h o ; R + 3 + 4 h o ; I   η + , s 4 ( η + , s 2 + 1 ) .
Our choice of the TP coefficients a = 0   ,     κ = + , and κ + =     a + = 1 assures that conditions (150) and (151) automatically hold and, as a result, the variables
z ( x )   = t a n h 2 ( 2 x )
and
z + ( x )   = s i n h 2 ( 2 x )
are related via linear fractional Formula (158).
The corresponding prime SLEs take the form
  ( 1 η 2 ) d     d η ( 1 η 2 ) d       d η V R M [ η ; λ o , μ o ] ε   Ψ [ η ; λ o , μ o ; ε ] = 0
and
  ( 1 + η 2 ) d     d η 1 + η 2 d       d η V G [ η ; h o ] ε ˜   Ψ i [ η ; h o ; ε ˜ ] = 0 ,
which are solved under the DBCs at ±1 and ±∞, respectively. The representation of the RCSLEs of our interest in the prime form [37] makes the ChExps of two Frobenius solutions differ only by sign, and therefore the DBC automatically selects the principal solution.
Setting κ + = 1 in (136) brings us to the following quadratic equation
ω i m 2 ( ω i m 2 h o ; R ) 1   4 h o ; I 2 = 0
in ω i m 2 , where we changed c n for m to indicate that the given equation holds for the averaged indexes of the Jacobi polynomials forming not only the eigenfunctions of CSLE (114), but also any other q-RSs. The crucial feature of this quadratic equation is that its coefficients are independent of the polynomial degree, and therefore, as a direct consequence of (126), this is also true for the complex-conjugated pairs of the Jacobi indexes in question. So, in contrast with RM potential (195), Gendenshtein potential (196) belongs to Group A (not B!) in Odake and Sasaki’s [66] classification scheme of the rational TSI potentials. This is the direct consequence of our observation [51] that the ExpDiffs for the poles of SLE (202) at ±∞ are energy-independent, so the TFI CSLE under consideration belongs to Group A [63]. Since all the q-RSs in this case are specified by a single series of Maya diagrams [63], all the solvable rational Darboux–Crum transforms of SLE (202) with a complex parameter h o can be obtained using admissible Wronskians of Routh polynomials [67] with the common complex index. It is convenient to choose the sequence which starts from a finite orthogonal set of R–Routh polynomials forming n o eigenfunctions of SLE (202). As a result [66,67], the corresponding Wronskian transforms of R–Routh polynomials form finite sequences of exceptional orthogonal polynomials (EOSs) in Quesne’s terms [64].
The symmetric reductions λ o = 0 and h o ; I = 0 of potentials (195) and (196) represent exactly the same sech-squared potential if the parameters μ o and h o ; R are related via (176). This potential can be thus quantized either via classical Gegenbauer polynomials with degree-dependent indexes [44] or via Masjed-Jamei polynomials [46] with degree-independent indexes [45]. As pointed to in [63], the A or B grouping suggested by Odake and Sasaki [66] for rational TSI potentials is actually an attribute of the particular rational realization of the given potential rather than the potential itself. The h-PT and Morse potentials represent two other instances of such an atypical dualism [63,68].
The latter quantization scheme using R–Routh polynomials of a definite parity is obviously preferable for constructing the symmetric Darboux–Crum transforms of the sech-squared potential [66]. Since it is the symmetric limit of the Gendenshtein potential, the q-RSs of a definite parity in this limiting case are also specified by a single series of Maya diagrams, and therefore all the rational Darboux–Crum transforms (RDC T ) of SLE (202) with a real parameter h o can be obtained using admissible Wronskians of Routh polynomials of a definite parity with the common real index such that the first n o polynomials in the given infinite polynomial sequence constitute a finite orthogonal set of Masjed-Jamei polynomials. As a result, the corresponding Wronskian transforms of Masjed-Jamei polynomials must form finite sequences of EOPs. A detailed study on this remarkable family of finite EOP sequences is currently under way.

6. Discussion

The presented analysis scrutinizes the striking resemblance between quartic Equations (58) and (125) suggested by Lévai [34] for computing the energy spectrum of real Fuschian CSLEs (1) and (114) with two poles located symmetrically around the origin on either the real or imaginary axis. In both cases, even density functions (66) and, accordingly, (110) remain positive within the corresponding quantization intervals, so the change of variable η ; s ( x ; δ , C ) or subsequently η + ; s ( x ; δ + , C + ) converting the given CSLE into the Schrödinger equation is an odd real function of x. It has been shown that the two variables are interrelated via simple algebraic Formula (148) or its reverse (149).
It was also proven that the corresponding Liouville potentials referred to in the paper as Lévai’s [34] and Milson’s [6] e-TP potentials overlap along the symmetric curves, which are nothing but two alternative rational representations of Ginocchio’s [44] potential function (178).
Since the Liouville transformation converts R Ref CSLE (114) with real h o into the Schrödinger equation with the Ginocchio potential on the line, the eigenfunctions of the latter equations can be expressed in terms of Masjed-Jamei polynomials [46] with degree-dependent index [45]. In the TSI limit represented by sech-squared potential, the index of the Masjed-Jamei polynomials forming the corresponding eigenfunctions becomes degree-independent and, as a result, the Wronskians of the eigenfunctions turn into weighted Wronskians of the orthogonal Masjed-Jamei polynomials with the common index. This finite polynomial sequence starts an infinite sequence of Routh polynomials of a definite parity with the same real index. If the Wronskian formed by a subset of the latter polynomials does not have real roots, then the corresponding Wronskian transforms of the Masjed-Jamei polynomials form a finite EOP sequence. The EOP sequences constructed in such a way represent the symmetric reduction of the EOP sequences generated by us in [67] for Gendenshtein potential (196) using the admissible Wronskians of Routh polynomials with a common complex index.

7. Conclusions

Though the main purpose of the paper was to reveal some remarkable overlapping features of two real-field reductions of the complex (non-self-adjoint) J Ref CSLE with an even density function, the starting point for our analysis has a much broader area of application. Namely, the density function does not have to be even—we only need to require that it remains positive inside the finite real interval or on the imaginary axis. The Liouville transformation on these intervals results in two families of complex Liouville potentials referred to by us as complex potentials of the Lévai class.
As suggested by Lévai [34], one can significantly simplify the analysis of the resultant Schrödinger equation by choosing the even density function. Such a choice of the density function assures that the variable η(x) used to convert the CSLE to the Schrödinger equation in x is an odd function of x. The direct consequence of this constraint is that the corresponding Liouville potential is P T -symmetric if this is true for the J RefPF in question.
The important advantage of the suggested complexification of the J Ref CSLE comes from the fact that it has a quartet of AEH solutions which can be used as FFs for RDTs whether the Liouville transformation is done on the real or imaginary axis. In particular, it will be shown in a separate publication that Lévai’s P T -symmetric potential [34,38,40] obtained via the Liouville transformation of the P T -symmetric J Ref CSLE on the real axis has AEH solutions composed of R–Routh polynomials in an imaginary argument i η I . Since these functions may not have zeros on the real axis < η I ( x ) < + , each of them can be used as the FF to construct a new solvable P T -symmetric potential. The sketched complexification of the J Ref CSLE thus provides a new mechanism for constructing P T -symmetric potentials with real energy spectra specified by positive roots of quartic Equation (58).
The presented analysis of the AEH solutions of the complex J Ref CSLE with the even density function thus opens new horizons for Lévai’s unified approach [34] to two families of generally complex Liouville potentials obtained via the Liouville transformations of the mentioned CSLE, either inside the interval (−1, +1) or on the imaginary axis.

Funding

This research received no external funding.

Data Availability Statement

No new data were created.

Acknowledgments

I wish to use this opportunity to thank Mariana Kirchbach for her inspiring advice, which greatly influenced the early stage of my studies on the quantization of the R Ref potential by Romanovski–Routh polynomials. The author is also grateful to the referees for their stimulating critical suggestions to tune the paper results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bose, A.K. Solvable potentials. Phys. Lett. 1963, 7, 245–246. [Google Scholar] [CrossRef]
  2. Bose, A.K. A class of solvable potentials. Nuovo Cim. 1964, 32, 679–688. [Google Scholar] [CrossRef]
  3. Erdelyi, A.; Bateman, H. Transcendental Functions; McGraw Hill: New York, NY, USA, 1953; Volume 1. [Google Scholar]
  4. Everitt, W.N. A Catalogue of Sturm-Liouville Differential Equations. In Sturm-Liouville Theory, Past and Present; Amrein, W.O., Hinz, A.M., Pearson, D.B., Eds.; Birkhäuser Verlag: Basel, Switzerland, 2005; pp. 271–331. [Google Scholar]
  5. Natanzon, G.A. Study of the one-dimensional Schrödinger equation generated from the hypergeometric equation. Vestn. Leningr. Univ. 1971, 10, 22–28. Available online: arxiv.org/PS_cache/physics/pdf/9907/9907032v1.pdf (accessed on 20 July 1999).
  6. Milson, R. Liouville transformation and exactly solvable Schrödinger equations. Int. J. Theor. Phys. 1998, 37, 1735–1752. [Google Scholar] [CrossRef]
  7. Natanson, G. Exact quantization of the Milson potential via Romanovski-Routh polynomials. arXiv 2015, arXiv:1310.0796v3. [Google Scholar]
  8. Avarez-Castillo, D.E.; Kirchbach, M. Exact spectrum and wave functions of the hyperbolic Scarf potential in terms of finite Romanovski polynomials. Rev. Mex. Fis. E 2007, 53, 143–154. [Google Scholar]
  9. Natanson, G. Routh Polynomials: Hundred Years in Obscurity. 2022. pp. 1–61. Available online: researchgate.net/publication/326522529 (accessed on 1 November 2022).
  10. Natanson, G. Rediscovery of Routh Polynomials after Hundred Years in Obscurity. In Recent Research in Polynomials; Özger, F., Ed.; IntechOpen: London, UK, 2023; 27p. [Google Scholar]
  11. Koekoek, R.; Lesky, P.A.; Swarttouw, R.F. Hypergeometric Orthogonal Polynomials and Their q-Analogues; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  12. Koornwinder, T. Additions to the Formula Lists in ‘Hypergeometric Orthogonal Polynomials and Their q-Analogues’ by Koekoek, Lesky and Swarttouw. 2022. pp. 1–45. Available online: Staff.fnwi.uva.nl/t.h.koornwinder/art/informal/KLSadd.pdf (accessed on 4 February 2022).
  13. Jordaan, K.; Toókos, F. Orthogonality and asymptotics of Pseudo-Jacobi polynomials for non-classical parameters. J. Approx. Theory 2014, 178, 1–12. [Google Scholar] [CrossRef]
  14. Routh, E.J. On some properties of certain solutions of a differential equation of second order. Proc. Lond. Math. Soc. 1884, 16, 245–261. [Google Scholar] [CrossRef] [Green Version]
  15. Lesky, P.A. Vervollstandigung der klassischen Orthogonalpolynome durch Erganzungen zum Askey-Schema der hypergeometrischen orthogonalen Polynome. Ost. Akad. Wiss. 1995, 204, 151–166. [Google Scholar]
  16. Lesky, P.A. Endliche und unendliche Systeme von kontinuierlichen klassichen Othogonalpolynomen. Z. Angew. Math. Mech. 1996, 76, 181–184. [Google Scholar] [CrossRef]
  17. Romanovsky, V. Sur la généralisation des courbes de Pearson. Atti Congresso Intern. Mat. 1928, 6, 107–110. [Google Scholar]
  18. Romanovski, V.I. Sur quelques classes nouvelles de polynomes orthogonaux. C. R. Acad. Sci. Paris 1929, 188, 1023–1025. [Google Scholar]
  19. Natanson, G. Exact quantization of the Milson potential via Romanovski polynomials. arXiv 2013, arXiv:1310.0796v1. [Google Scholar]
  20. Dereziński, J.; Wrochna, M. Exactly solvable Schrödinger Operators. Ann. Henri Poincaré 2011, 12, 397–418. [Google Scholar] [CrossRef] [Green Version]
  21. Stevenson, A.F. Note on the ‘Kepler problem’ in a spherical space, and the factorization method of solving eigenvalue problems. Phys. Rev. 1941, 59, 842–843. [Google Scholar] [CrossRef]
  22. Cryer, C.W. Rodrigues’ formulas and the classical orthogonal polynomials. Boll. Unione Mat. Ital. 1970, 25, 1–11. [Google Scholar]
  23. Askey, R. An integral of Ramanujan and orthogonal polynomials. J. Indian Math. Soc. 1987, 51, 27–36. [Google Scholar]
  24. Bagrov, V.G.; Gitman, D.M.; Ternov, I.M.; Khalilov, V.R.; Shapovalov, V.N. Tochnye Reshenya Relativiskikh Volnovykh Uravnenii; Nauka: Novosibirsk, Russia, 1982. (In Russian) [Google Scholar]
  25. Cooper, F.; Khare, A.; Sukhatme, U.P. Supersymmetry and quantum mechanics. Phys. Rep. 1995, 251, 267–385. [Google Scholar] [CrossRef] [Green Version]
  26. Cooper, F.; Khare, A.; Sukhatme, U.P. Supersymmetry in Quantum Mechanics; World Scientific: Denver, CO, USA, 2001. [Google Scholar]
  27. Compean, C.B.; Kirchbach, M. The trigonometric Rosen–Morse potential in the supersymmetric quantum mechanics and its exact solutions. J. Phys. A 2006, 39, 547–557. [Google Scholar] [CrossRef] [Green Version]
  28. Rosen, N.; Morse, P.M. On the vibrations of polyatomic molecules. Phys. Rev. 1932, 42, 210–217. [Google Scholar] [CrossRef]
  29. Dabrowska, J.W.; Khare, A.; Sukhatme, U.P. Explicit wavefunctions for shape-invariant potentials by operator techniques. J. Phys. A 1988, 21, L195–L200. [Google Scholar] [CrossRef]
  30. Gendenshtein, L.E. Derivation of exact spectra of the Schrödinger equation by means of supersymmetry. JETP Lett. 1983, 38, 356–359. [Google Scholar]
  31. Lévai, G. A search for shape-invariant solvable potentials. J. Phys. A 1989, 22, 689–702. [Google Scholar] [CrossRef]
  32. Bhattacharjie, A.; Sudarshan, E.C.G. A class of solvable potentials. Nuovo Cim. 1962, 25, 864–879. [Google Scholar] [CrossRef]
  33. Grosche, C. The general Besselian and Legendrian path integrals. J. Phys. A 1996, 29, L183–L189. [Google Scholar] [CrossRef]
  34. Lévai, G. Gradual spontaneous breakdown of 𝒫𝒯 symmetry in a solvable potential. J. Phys. A 2012, 45, 444020. [Google Scholar] [CrossRef]
  35. Gibbons, J.; Veselov, A.P. On the rational monodromy-free potentials with sextic growth. J. Math. Phys. 2009, 50, 013513. [Google Scholar] [CrossRef]
  36. Natanson, G. Survey of nodeless regular almost-everywhere holomorphic solutions for exactly solvable Gauss-reference Liouville potentials on the line I. Subsets of nodeless Jacobi-seed solutions co-existent with discrete energy spectrum. arXiv 2016, arXiv:1606.08758. [Google Scholar]
  37. Natanson, G. Darboux-Crum Nets of Sturm-Liouville Problems Solvable by Quasi-Rational Functions I. General Theory. 2018. pp. 1–99. Available online: researchgate.net/publication/323831953 (accessed on 1 March 2018).
  38. Lévai, G. PT Symmetry in Natanzon-class Potentials. Int. J. Theor. Phys. 2015, 54, 2724–2736. [Google Scholar] [CrossRef] [Green Version]
  39. Lévai, G. A class of exactly solvable potentials related to the Jacobi polynomials. J. Phys. A 1991, 24, 131–146. [Google Scholar] [CrossRef]
  40. Lévai, G. Accidental crossing of energy eigenvalues in PT-symmetric Natanzon-class potentials. Ann. Phys. 2017, 38, 1–11. [Google Scholar] [CrossRef]
  41. Lévai, G. Exactly Solvable PT-Symmetric Models. In PT Symmetry in Classical and Quantum Physics; Bender, C., Ed.; World Scientific: London, UK, 2019; pp. 221–260. [Google Scholar]
  42. Gesztesy, F.; Simon, B.; Teschl, G. Zeros of the Wronskian and renormalize oscillation theory. Am. J. Math. 1996, 118, 571–594. [Google Scholar]
  43. Garcia-Ferrero, M.; Gȯmez-Ullate, D.; Milson, R. A Bochner type classification theorem for exceptional orthogonal polynomials. J. Math. Anal. Appl. 2019, 472, 584–626. [Google Scholar] [CrossRef]
  44. Ginocchio, J. A class of exactly solvable potentials: I. One-dimensional Schrödinger equation. Ann. Phys. 1984, 152, 203–219. [Google Scholar] [CrossRef]
  45. Natanson, G. Quantization of One-Dimensional Ginocchio Potential by Masjed-Jamei Polynomials with Degree-Dependent Indexes. 2022. pp. 1–52. Available online: researchgate.net/publication/360950671 (accessed on 1 May 2022).
  46. Masjedjamei, M. Three finite classes of hypergeometric orthogonal polynomials and their application in functions approximation. Integr. Transf. Spec. Funct. 2002, 13, 169–190. [Google Scholar] [CrossRef]
  47. Natanzon, G.A. Construction of the Jost function and the S-matrix for a general potential allowing solution of the Schrödinger equation in terms of hypergeometric functions. Sov. Phys. J. 1978, 21, 855–859. [Google Scholar] [CrossRef]
  48. Natanzon, G.A. General properties of potentials for which the Schrödinger equation can be solved by means of hypergeometric functions. Theor. Math. Phys. 1979, 38, 146–153. [Google Scholar] [CrossRef]
  49. Natanson, G. Single-source nets of Fuschian rational canonical Sturm-Liouville equations with common simple-poles density functions and related sequences of multi-indexed orthogonal Heine eigenpolynomials. In Proceedings of the 14th International Symposium on Orthogonal Polynomials, Special Functions and Applications, Canterbury, UK, 3–7 July 2017; Available online: https://www.researchgate.net/publication/317643178 (accessed on 1 July 2017).
  50. Natanson, G. Use of Wronskians of Jacobi Polynomials with Common Complex Indexes for Constructing X-DPSs and Their Infinite and Finite Orthogonal Subsets. 2019. pp. 1–57. Available online: researchgate.net/publication/331638063 (accessed on 1 March 2019).
  51. Natanson, G. Gauss-seed nets of Sturm-Liouville problems with energy-independent characteristic exponents and related sequences of exceptional orthogonal polynomials I. Canonical Darboux transformations using almost-everywhere holomorphic factorization functions. arXiv 2013, arXiv:1305.7453v1. [Google Scholar]
  52. Natanson, G. Single-source nets of algebraically-quantized reflective Liouville potentials on the line I. Almost-everywhere holomorphic solutions of rational canonical Sturm-Liouville equations with second-order poles. arXiv 2015, arXiv:1503.04798v2. [Google Scholar]
  53. Everitt, W.N.; Littlejohn, L.L. Orthogonal polynomials and spectral theory: A survey. In Orthogonal Polynomials and Their Applications; IMACS Annals on Computing and Applied, Mathematics; Brezinski, C., Gori, L., Ronveaux, A., Eds.; J. C. Baltzer AG Publishers: Tulsa, OK, USA, 1991; Volume 9, pp. 21–55. [Google Scholar]
  54. Everitt, W.N.; Kwon, K.H.; Littlejohn, L.L.; Wellman, R. Orthogonal polynomial solutions of linear ordinary differential equations. J. Comp. Appl. Math. 2001, 133, 85–109. [Google Scholar] [CrossRef] [Green Version]
  55. Bochner, S. Über Sturm-Liouvillesche Polynomsysteme. Math. Z. 1929, 29, 730–736. [Google Scholar] [CrossRef]
  56. Natanson, G. Dutt-Khare-Varshni Potential and Its Quantized-by-Heun-Polynomial SUSY Partners as Nontrivial Examples of Solvable Potentials Explicitly Expressible in Terms of Elementary Functions. arXiv 2015, arXiv:1508.04738v1. [Google Scholar]
  57. Courant, R.; Hilbert, D. Methods of Mathematical Physics; Interscience: New York, NY, USA, 1953; Volume 1. [Google Scholar]
  58. Natanson, G. Breakup of SUSY quantum mechanics in the Limit-Circle region of the reflective Kratzer oscillator. arXiv 2014, arXiv:1405.2059v1. [Google Scholar]
  59. Lévai, G. Non-compact groups and solvable potentials. In Proceedings of the XXI International Colloquium on Group Theoretical Methods in Physics, Goslar, Germany, 15–20 July 1996; Doebner, H.-D., Nattermann, P., Scherer, W., Eds.; World Scientific: Singapore, 1997; Volume 1, pp. 461–466. [Google Scholar]
  60. Wu, J. Group Theory Approach to Scattering. Ph.D. Thesis, Yale University, New Haven, CT, USA, 1985. [Google Scholar]
  61. Wu, J.; Alhassid, Y.; Gùrsey, F. Group theory approach to scattering. IV. Solvable potentials associated with SO(2,2). Ann. Phys. 1989, 196, 163–181. [Google Scholar] [CrossRef]
  62. Ginocchio, J. A Class of exactly solvable potentials II. The three-dimensional Schrodinger equation. Ann. Phys. 1985, 159, 467–480. [Google Scholar] [CrossRef]
  63. Natanson, G. Equivalence Relations for Darboux-Crum Transforms of Translationally Form-Invariant Sturm-Liouville Equations. 2021. pp. 1–58. Available online: researchgate.net/publication/353131294 (accessed on 1 August 2021).
  64. Quesne, C. Extending Romanovski polynomials in quantum mechanics. J. Math. Phys. 2013, 54, 122103. [Google Scholar] [CrossRef] [Green Version]
  65. McIntosh, H.V. Complex Analysis; Universidad Autónoma de Puebla: Puebla, Mexico, 2001. [Google Scholar]
  66. Odake, S.; Sasaki, R. Krein–Adler transformations for shape-invariant potentials and pseudo virtual states. J. Phys. A 2013, 46, 245201. [Google Scholar] [CrossRef] [Green Version]
  67. Natanson, G. On Finite Exceptional Orthogonal Polynomial Sequences Composed of Wronskian Transforms of Romanovski-Routh Polynomials. 2022. pp. 1–38. Available online: https://www.researchgate.net/publication/364350392 (accessed on 1 October 2022).
  68. Natanson, G. Quantization of rationally deformed Morse potentials by Wronskian transforms of Romanovski-Bessel polynomials. Acta Polytec. 2022, 62, 100–117. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Natanson, G. Overlapping of Lévai’s and Milson’s e-Tangent-Polynomial Potentials along Symmetric Curves. Axioms 2023, 12, 584. https://doi.org/10.3390/axioms12060584

AMA Style

Natanson G. Overlapping of Lévai’s and Milson’s e-Tangent-Polynomial Potentials along Symmetric Curves. Axioms. 2023; 12(6):584. https://doi.org/10.3390/axioms12060584

Chicago/Turabian Style

Natanson, Gregory. 2023. "Overlapping of Lévai’s and Milson’s e-Tangent-Polynomial Potentials along Symmetric Curves" Axioms 12, no. 6: 584. https://doi.org/10.3390/axioms12060584

APA Style

Natanson, G. (2023). Overlapping of Lévai’s and Milson’s e-Tangent-Polynomial Potentials along Symmetric Curves. Axioms, 12(6), 584. https://doi.org/10.3390/axioms12060584

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop