Next Article in Journal
Low-Order Moments of Velocity Gradient Tensors in Two-Dimensional Isotropic Turbulence
Previous Article in Journal
Approximate N5LO Higgs Boson Decay Width Γ(Hγγ)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unified Supersymmetric Description of Shape-Invariant Potentials within and beyond the Natanzon Class

1
Szeberényi Gusztáv Adolf Evangélikus Gimnázium, 5600 Békéscsaba, Hungary
2
Faculty of Science and Technology, University of Debrecen, 4032 Debrecen, Hungary
3
HUN-REN Institute for Nuclear Research (ATOMKI), P.O. Box 51, 4001 Debrecen, Hungary
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(2), 174; https://doi.org/10.3390/sym16020174
Submission received: 28 December 2023 / Revised: 25 January 2024 / Accepted: 30 January 2024 / Published: 1 February 2024
(This article belongs to the Section Physics)

Abstract

:
The transformations of supersymmetric quantum mechanics are discussed within a formalism that employs a six-parameter function, from which the superpotential and the supersymmetric partner potentials V ( x ) and V + ( x ) are constructed in a general form. By specific choice of the parameters, V ( x ) and V + ( x ) are matched with the general form of PI class potentials and their rationally extended versions. The choice of the parameters also determines which of the four possible SUSY transformations T i , i = 1 , 4 is in effect. After this general discussion, the formulae are specified to the three members of this potential class, the Scarf I, Scarf II and generalized Pöschl–Teller potentials. Due to the different domains of definition and their consequences on the boundary conditions, the results turn out to be rather diverse for the three potentials, while the mathematical formalism and the network of the potentials interconnected by the SUSYQM transformations still remains common to a large extent. The general framework allows a unified and consistent interpretation of earlier isolated findings. It also helps to connect the results to further potential classes and to place them into a more general context within the zoo of exactly solvable potentials.

1. Intoduction

Simple physical models have proven to be useful tools to understand the basic principles of the physical world. By increasing the complexity of these models, their range of applicability can be extended; however, these efforts are usually limited after some point by the increasing technical difficulties. Finding the point up to which the models improve without becoming overly complicated is a real challenge in general. An example for this scenario is studying exactly solvable quantum mechanical potentials. This field has been active since the introduction of modern quantum mechanics a century ago. The first simple models obtained by the exact solution of the Schrödinger equation were the harmonic oscillator and the Coulomb, Morse and Pöschl–Teller potentials, for example. Some of these problems have been incorporated into more developed models and have also become part of introductory quantum mechanical courses. The range of exactly solvable potentials still keeps extending, thanks to new innovative ideas and technical developments that appear from time to time (sometimes recycling or generalizing earlier ideas and developments). The technical difficulties one encounters in these efforts are usually of a mathematical nature. Fortunately, one can rely on the rich arsenal of mathematical physics, especially the theory of differential equations and special functions compiled well before the introduction of quantum mechanics. As the range of exactly solvable quantum mechanical potentials increased, the demand for finding principles along which they can be classified also appeared in a natural way.
The most straightforward and most widespread method of solving the Schrödinger equation (and also some further fundamental equations of physics) is mapping it to some second-order differential equation of mathematical physics by means of a well-chosen variable transformation. Many exactly solvable potentials have been found in this way. However, one cannot be sure whether these random and isolated results exhaust the set of all exactly solvable problems or not. This problem was handled by reversing the procedure: rather than setting out from a particular potential and arriving at a known differential equation, systematic procedures were formulated by which—in principle—all the exactly solvable problems with solutions containing a single special function of mathematical physics could be derived. These efforts resulted, for example, in the introduction of Natanzon potentials [1], for which the wave function is expressed in terms of a hypergeometric function. Natanzon confluent potentials [2] are related to the confluent hypergeometric function. The bound-state solutions of these problems are written in terms of Jacobi and generalized Laguerre polynomials, which are obtained from the named functions as special cases [3].
Applying this method also implied a natural classification of exactly solvable potentials in terms of the special function used in the particular case. However, this first level of classification leaves all the potentials within a potential class unidentified. As a second level of classification, the variable transformations applied in the mapping procedure proved to be a natural choice. Such an effort was proposed in 1962 by observing that a specific differential equation can be formulated for the function governing the variable transformation [4] by requiring that a constant term (corresponding to the energy eigenvalue) has to appear in the resulting transformed equation (the Schrödinger equation). The systematic identification of these secondary potential classes has been introduced [5] for potentials with bound-state wave functions containing Jacobi and generalized Laguerre polynomials (as well as Hermite polynomials, which can be considered special cases of the latter polynomials). The relation of these potential classes with other classification schemes of potentials has also been discussed in Ref. [5].
Another rather successful method of discussing exactly solvable potentials is supersymmetric quantum mechanics (SUSYQM) [6] (for a review, see Refs. [7,8], for example). In this approach, the focus is not on the structure of the solutions; rather, it is based on a symmetry, which implies a degeneracy of the energy levels of two potentials that are connected by a specific supersymmetric transformation that can be formulated by using a known particular (not necessarily physical) solution of one of the potentials. Supersymmetric quantum mechanics became a rather successful method of finding new exactly solvable potentials from an already known one. Not surprisingly, SUSYQM is also rooted in earlier theories, e.g., the Darboux transformation [9], going back well before the time of quantum mechanics.
SUSYQM can be used to generate new potentials both within and outside the range of exactly solvable potentials: it is equally applicable in numerical calculations. However, applying it to known exactly solvable potentials helped to gain further insight into their structure and classification. One such case was when the two potentials connected by SUSYQM turned out to have the same mathematical structure, with different parameters. This is how the concept of shape invariance was introduced [10]. The most well-known textbook examples for exactly solvable potentials (e.g., the harmonic oscillator, Coulomb, Morse, Pöschl–Teller, etc. potentials) were found to exhibit shape-invariance. In fact, these are the potentials that arise within the two-step classification procedure mentioned before as the most natural examples [5], and also in the factorization method [11] (introduced by E. Schrödinger [12,13]), which can be considered a precursor of the concept of shape-invariance.
Another notable case is when the two potentials connected by SUSYQM belong to two different classes in the variable transformation framework, i.e., their solutions are different special functions of mathematical physics. Examples of this scenario are the rationally extended versions of some of the shape-invariant potentials. The solutions of these potentials are exceptional orthogonal polynomials [14,15], in particular, the X 1 -type Jacobi and Laguerre polnomials. These differ from their classical counterparts in that the degree of the lowest polynomial in the sequence is one, rather than zero. The range of the potentials derived from these exceptional polynomials contain the rationally extended Scarf I [16], Pöschl–Teller [17] and Scarf II [18] potentials for the X 1 -type Jacobi and the rationally extended harmonic oscillator [16] for the X 1 -type Laguerre polynomials. As mentioned above, these new potentials can be obtained by SUSYQM transformation, with their classical counterparts as starting potentials. Furthermore, these rationally extended potentials are characterized by the shape-invariance property. Certainly, these potentials are outside the Natanzon class, because their solutions cannot be expressed in terms of the hypergeometric function, the confluent hypergeometric function or polynomials to which these functions can be reduced. However, the hypergeometric and the confluent hypergeometric functions can also be expressed as special cases of more general functions, e.g., the confluent Heun function [19]. It can be proven that this function can be specified in polynomial form in various ways, and both the X 1 -type Jacobi and the X 1 -type Laguerre polynomials can be obtained from the confluent Heun function. See, e.g., Refs. [20,21]. We also note that some of the potentials discussed here can also be obtained from the transformation of the general (rather than the confluent) Heun equation [22].
Here, we combine the two methods outlined above to discuss a number of potentials both within and outside the Natanzon class. We introduce a general formalism based on a rather general six-parameter function that satisfies the requirements of SUSYQM and introduce two potentials that are SUSYQM partners by construction. We then select the six parameters in such a way that the two potentials coincide with the Scarf I, Scarf II and generalized Pöschl–Teller potentials and/or their rational extension. The chosen parameters set the behavior of the original six-parameter function at the boundaries, which, in turn, defines the type of SUSYQM transformation and decides whether the transformation corresponds to broken or unbroken supersymmetry.
This project has been discussed already for the radial harmonic oscillator and its rational extension [23]. The connection between these potentials by different types of SUSYQM transformations has been analyzed, and their solutions, as well as their shape-invariance, were discussed in a unified and consistent way. In the present study, we implement this procedure to the Scarf I, Scarf II and the generalized Pöschl–Teller potentials and their rational extensions. Furthermore, the discussion is carried out on a higher level of abstraction in the sense that the formalism is developed in a general form valid for PI class potentials, to which all three mentioned potentials belong [5]. This potential class corresponds to type A potentials in the scheme of the factorization method [11], although it was not complete at the time of publication of the latter work. These potentials also belong to the same class in other classification schemes: (b) in Ref. [24] and A in Ref. [25]. The individual potentials are specified only at the stage where the domain of definition of the problems and the corresponding boundary conditions start to play a role: these differ in the case of the three potentials, while their general mathematical structure remains rather similar, as do also part of the conclusions. With this, a number of isolated earlier results appear in a general consistent framework. This general approach also helps to interpret the role of the exceptional polynomials; these are relatively new actors within the traditional formalisms that have been used previously to discuss exactly solvable potentials.
This paper is organized as follows. In Section 2, the basic properties of SUSYQM are reviewed, and the essential formulae for the general PI class potentials and their rational extension are collected for further reference. In Section 3, the χ ( r ) function generating the SUSYQM transformation is introduced, and the most general form of the two partner potentials is derived. The partner potentials are then matched with the general PI class potential and its rational extension. The general results are specified for the three particular potentials in Section 4. Finally, in Section 5, the results are summarized.

2. Preliminaries

Here, we discuss some general aspects of the solutions of the radial Schrödinger equation (using the units = 2 m = 1 )
d 2 ψ d x 2 + [ E V ( x ) ] ψ ( x ) = 0 ,
where x [ x , x + ] , a domain which will be specified later. Normalizable solutions, i.e., bound-state wave functions, have to satisfy the boundary conditions ψ ( x ) = ψ ( x + ) = 0 .

2.1. The Basics of Supersymmetric Quantum Mechanics (SUSYQM)

Supersymmetric theories handle bosonic and fermionic quantities within a joint formalism and contain operators that satisfy a set of commutation and anticommutation relations. This is also the case in supersymmetric quantum mechanics, which was introduced as a toy model to study the breakdown of supersymmetry [6], but later evolved into a rather successful and robust theory in quantum mechanics. It was found that a handful of operators satisfying a few simple equations can lead to a rich variety of results concerning quantum mechanical potential problems [7,8]. This is also the case for N = 2 SUSYQM, which is based on the superalgebra
{ Q , Q } = H , Q 2 = ( Q ) 2 = 0 , [ Q , H ] = [ Q , H ] = 0 .
Here H is the supersymmetric Hamiltonian, while Q and Q are called the supersymmetric charge operators. A simple realization of this superalgebra makes use of 2 × 2 matrices with operator entries
Q = 0 0 A 0 , Q = 0 A 0 0 ,
H = A A 0 0 A A H 0 0 H + .
H and H + with eigenfunctions ψ ( ) and ψ ( + ) , respectively, are traditionally called the bosonic and fermionic Hamiltonians, respectively, and are interpreted as supersymmetric partners. The supersymmetric charge operators connect the bosonic and fermionic “sectors” via
Q ψ ( ) 0 = 0 A ψ ( ) , Q 0 ψ ( + ) = A ψ ( + ) 0 .
A straightforward way to obtain the one-dimensional Schrödinger equation as a second-order linear differential equation is its factorization in terms of two first-order differential operators. The formalism of SUSYQM is suitable for this if one defines the A and A operators as
A = d d x + W ( x ) and A = d d x + W ( x ) .
Then, the bosonic and fermionic Hamiltonians take the form
H ± ψ ( ± ) ( x ) = d 2 d x 2 + V ± ( x ) ψ ( ± ) ( x ) = E ( ± ) ψ ( ± ) ( x ) .
The V ( x ) and V + ( x ) potentials are called supersymmetric partner potentials and are obtained from the superpotential W ( x ) as
V ± ( x ) = W 2 ( x ) ± d d x W ( x ) .
Straightforward calculations reveal that H and H + have common energy eigenvalues. In particular, if ψ ( ) ( x ) ( ψ ( + ) ( x ) ) is the eigenfunction of H ( H + ), then A ψ ( ) ( x ) ( A ψ ( + ) ( x ) ) is the eigenfunction of H + ( H ).
It is customary to set the energy scale such that the ground state of H is at E = 0 , i.e., H ψ 0 ( ) ( x ) = 0 . This result also follows from the relation
A ψ 0 ( ) = 0 .
It can be shown that there will be no eigenstate of H + at E = 0 , while its other energy eigenvalues coincide with those of H , i.e.,
E n + 1 ( ) = E n ( + ) ( n = 0 , 1 , 2 , ) with E 0 ( ) = 0 .
Furthermore, Equation (9) allows expressing the superpotential in terms of the ground-state wave function of H as
W ( x ) = d d x ln ψ 0 ( ) ( x ) .
In summary, in the knowledge of the bound-state wave functions of a potential V ( x ) , one can construct its supersymmetric partner potential V + ( x ) with (almost) the same energy spactrum, and its bound-state wave functions can also be easily obtained from those of V ( x ) .
Supersymmetric partner potentials can also be generated from further, not necessarily physical, solutions of H [26]. Taking the solution χ ( x ) with eigenvalue ϵ called the factorization energy [27],
H χ ( x ) A A χ ( x ) = ϵ χ ( x ) .
Equations (7) and (8) lead to the potential
V ( x ) = χ ( x ) χ ( x ) + ϵ = χ ( x ) χ ( x ) 2 χ ( x ) χ ( x ) + ϵ ,
which will be a continuous function of x, provided that χ ( x ) is nodeless. The analogues of Equations (8) and (11) are then
V ± ( x ) = W ˜ 2 ( x ) ± d d x W ˜ ( x ) + ϵ
and
W ˜ ( x ) = d d x ln χ ( x ) .
Note the presence of W ^ ( x ) = χ ( x ) / χ ( x ) in Equation (13). It is often possible to factorize χ ( x ) as
χ ( x ) = ψ 0 ( ) ( x ) ξ ( x ) ,
which leads to the relation connecting W ( x ) in (11) and W ˜ ( x ) in (15):
W ˜ ( x ) = W ( x ) d d x ln ξ ( x ) .
The ϵ = 0 choice recovers the special case χ ( x ) = ψ 0 ( ) ( x ) , which is interpreted as unbroken supersymmetry, while for ϵ 0 supersymmetry is said to be broken. For broken SUSY, H + will also possess a level at E = E 0 ( ) .
The degeneracy of the energy levels of H and H + ultimately originates from the superalgebra (2), which contains relations implying that the supersymmetric Hamiltonian H commutes with the supersymmetric charge operators. A Hamiltonian commuting with another operator indicates a symmetry, which generally implies a degeneracy. The origin of the usage of broken and unbroken supersymmetry is also related to these operators: the two-component ground state is the eigenfunction of H , but it can also be the eigenfunction of Q if A ψ ( ) = 0 holds, i.e., for ϵ = 0 . For A ψ ( ) 0 , this does not occur, so supersymmetry is spontaneously broken.
The behavior of the possible nodeless χ ( x ) solutions (12) at the boundaries, i.e., their normalizability or non-normalizability, is also related to the question of whether SUSY is broken or not in the actual situation. In fact, depending on the boundary conditions satisfied by χ ( x ) , four different SUSYQM transformations can be defined, as shown in Table 1 [26,28,29,30]. These have characteristic effects on the energy spectrum. The results are generally not affected by the domain of definition of the potential [ x , x + ] , i.e., whether it is defined on a finite domain ( < x , x + < ), on the full real x axis ( x , x + < ) or on the positive real semi-axis ( x = 0 , x + < ). The last case contains radial potentials, where it is customary to denote x with r. When x or x + is finite, then it is advisable to inspect the possible singularities of the problem at these points. Let us assume that x = 0 and the solution behaves like χ ( x ) x σ there. It is straightforward to demonstrate by taking χ ( x ) that the potential behaves like V ( x ) σ ( σ 1 ) x 2 there. This is the case, for example, for radial potentials, when the centrifugal term is l ( l + 1 ) r 2 , while σ = l + 1 for regular solutions and σ = l for irregular ones, where l is the orbital angular momentum with non-negative integer values. (Note that this applies to radial problem in three spatial dimensions.) The T 1 and T 3 transformations increase σ by one unit, while the T 2 and T 4 transformations decrease it by one unit. In the latter two cases, one has to consider σ > 1 in order to avoid attractive singularities in the SUSY partner potential V + ( x ) at x = 0 .
The same prescription applies to potentials corresponding to non-integer values of σ . Furthermore, the singularities have to be handled in a similar way at both boundaries of problems defined on a finite domain (i.e., with finite x and x + ). In these cases, x x and x + x take the role of x. For the subtleties of singular potentials, see, e.g., Ref. [31].
When x and/or x + is infinite, the boundary conditions are handled in a different way. If the potentials are assumed to decay exponentially in the asymptotic limit(s), then the solutions also follow an exponential tail: χ ( x ) exp ( κ x ) . Then, the sign of κ determines whether χ ( x ) converges or diverges asymptotically.
It has to be mentioned that for the transformations T 2 , T 3 and T 4 ϵ , the factorization energy can take on any value below E 0 ( ) , the ground-state energy of H . In these cases, it also appears in the spectrum of H + , and it can also be considered a new parameter of V + ( x ) . One finds that some ϵ values play a special role for a certain class of potentials: in these cases, the mathematical form of V + ( x ) and V ( x ) turns out to be similar, and only the parameters differ. This is the case for shape-invariant potentials mentioned previously.

2.2. The General PI Class Potentials and Their Rational Extensions

Here, we collect the essential formulae for the potentials discussed in the present work, both in their general form as PI class potentials [5,32] and in their form specified for the three concrete cases, i.e., the generalized Pöschl–Teller, Scarf I and Scarf II potentials and their rationally extended version. The procedure of obtaining the general form of these potentials by transforming the differential equation of the Jacobi polynomials and that of the X 1 -type exceptional Jacobi polynomials is presented in Appendix A.
The expression for the general PI type potentials [5,32] is
V ( α , β , x ) = C 1 z 2 ( x ) α + β 2 2 + α β 2 2 1 4 + 2 C z ( x ) 1 z 2 ( x ) α + β 2 α β 2
while the bound-state energy wave functions and the corresponding energy eigenvalues are written, respectively, as
ψ n ( α , β ; x ) = C n ( α , β ) ( 1 z ( x ) ) α 2 + 1 4 ( 1 + z ( x ) ) β 2 + 1 4 P n ( α , β ) ( z ( x ) )
and
E n = C n + α + β + 1 2 2 .
Here P n ( α , β ) ( z ) is a Jacobi polynomial [3]. Conditions for normalizability and regularity at the boundaries follow from the boundary conditions, which can be defined only after specifying the z ( x ) function and its domain of definition.
As discussed in Appendix A, there are three particular solutions belonging to the PI potential class, each with a different domain of definition. The generalized Pöschl–Teller potential is obtained for z ( x ) = cosh ( a x ) with C = a 2 , and it is defined on the positive real axis x [ 0 , ) corresponding to z [ 1 , ) . The Scarf I potential is the trigonometric version of the generalized Pöschl–Teller potential defined on a finite domain: z ( x ) = cos ( a x ) , C = a 2 , x [ 0 , π / a ] , z [ 1 , 1 ] . This potential is often derived from the z ( x ) = sin ( a x ) function defined on the x [ π / ( 2 a ) , π / ( 2 a ) ] domain, but this is merely a shifted version of the original potential. Finally, the Scarf II potential can be obtained by the substitutions z ( x ) = i sinh ( a x ) , C = a 2 , x ( , ) and z running along the whole imaginary axis. Note that although z ( x ) is an imaginary function, potential (18) will be real for an appropriate choice of α and β , as will also the bound-state wave functions (19), apart from an unimportant phase. We note that there are two more potentials, named Pöschl–Teller I and II, that are traditionally considered distinct entities; however, these can be derived from the Scarf I and the generalized Pöschl–Teller potentials by a trivial mapping x x / 2 [7].
According to the results derived in Appendix A, the rationally extended versions of the same potentials are obtained for the same three choices of z ( x ) and C, so their domain of definition will be the same as their classical correspondents. The general formulae corresponding to Equations (18)–(20) are as follows:
V ^ ( α , β , x ) = V ( α , β ; x ) + 2 C ( α + β ) ( α β ) z ( x ) + α + β + 2 C ( α β ) 2 ( α + β ) 2 ( α β ) z ( x ) + α + β 2 ,
ψ ^ n ( α , β ; x ) = C ^ n ( α , β ) ( 1 z ( x ) ) α 2 + 1 4 ( 1 + z ( x ) ) β 2 + 1 4 [ ( α β ) z ( x ) + α + β ] 1 P ^ n ( α , β ) ( z ( x ) ) ,
and
E ^ n = C n 1 + α + β + 1 2 2 .
Here, P ^ n ( α , β ) ( z ) are the X 1 -type exceptional Jacobi polynomials [14,15] discussed in Appendix A. Note that, of the three (five) linearly independent terms in Equation (A17) (Equation (A28)), only two (four) remain in Equation (19) (Equation (22)), because the third (fifth) constant term was incorporated into the energy E.
Before closing this section, we note that the P n ( α , β ) ( z ) Jacobi and P ^ n ( α , β ) ( z ) X 1 -type exceptional Jacobi polynomials form an infinite orthonormal polynomial system only when z [ 1 , 1 ] . For z outside this domain, they form a finite orthonormal system. See, e.g., Ref. [33].

3. Unified Discussion

In this section, we consider a particular form of the solution χ ( x ) in Equation (12). It is constructed to be nodeless and singularity-free within [ x , x + ] , in accordance with the requirements outlined in Section 2.1. In the next step, we use Equation (15) to generate the corresponding superpotential W ˜ ( x ) (from which we drop the tilde from now on, to ease the notation). Then, we express the V ( x ) and V + ( x ) partner potentials using Equation (14). Next, we analyze the mathematical form of these potentials and attempt to set the parameters such that V ( x ) and V + ( x ) take the form of either the general PI type potential (18) or its rationally extended version (21). These parameters also determine the behavior of χ ( x ) at the boundaries, also identifying the type of SUSYQM transformation in Table 1. We also calculate the factorization energy ϵ . This can be conducted up to a constant energy shift, unless one requires that the potentials be free from additive constants.

3.1. The General Results

The discussion starts with selecting a sufficiently general χ ( x ) function that reflects the requirements related to the boundary conditions specific to the possible T i transformations. Let us consider the χ ( x ) function depending on six parameters
χ ( x ) = ( 1 z ( x ) ) t ( 1 + z ( x ) ) s ( p + z ( x ) ) k ( q + z ( x ) ) j .
t and s (and especially their sign) determine the behavior of χ ( x ) at the boundaries z = 1 and z = 1 , while t + s + k j sets its asymptotic properties for z ± (whichever is applicable). This behavior also determines which of the four types of SUSY transformations χ ( x ) corresponds to from Table 1. Note that sometimes it is more natural to apply ( z ( x ) 1 ) rather than ( 1 z ( x ) ) ; however, this change in sign does not have a significant effect on the results. We also choose k and j to be non-negative integers and assume that p q : the p = q case is included in the formula as the j = 0 and k k j choice. The 1 / χ ( r ) operation that reverses the behavior of the solution at the boundaries and changes W ( r ) to W ( r ) can also be obtained by the ( t , s , k , j , p , q ) ( t , s , j , k , q , p ) transformation. The natural requirement that χ ( x ) has to be nodeless and singularity-free leads to the prescriptions z ( x ) p and z ( x ) q .
According to Equation (15), the superpotential is then given by
W ( x ) = C 1 / 2 ( 1 z 2 ( x ) ) 1 / 2 t 1 z ( x ) + s 1 + z ( x ) + k p + z ( x ) j q + z ( x ) ,
where the pre-factor originates from z ( x ) in Equation (A15). From this result and Equation (14), the partner potentials are obtained as
V ( x ) ϵ = C ( t + s + k j ) 2 + C 1 z 2 ( x ) ( 2 t 2 + 2 s 2 t s ) + C z ( x ) 1 z 2 ( x ) ( 2 t 2 2 s 2 t + s ) + C k ( k 1 ) ( 1 p 2 ) ( p + z ( x ) ) 2 + C j ( j + 1 ) ( 1 q 2 ) ( q + z ( x ) ) 2 + C k p + z ( x ) p ( 2 t + 2 s + 2 k 1 ) 2 t + 2 s 2 j 1 p 2 q p + C j q + z ( x ) q ( 2 t + 2 s 2 j 1 ) + 2 t 2 s + 2 k 1 q 2 q p
and
V + ( x ) ϵ = C ( t + s + k j ) 2 + C 1 z 2 ( x ) ( 2 t 2 + 2 s 2 + t + s ) + C z ( x ) 1 z 2 ( x ) ( 2 t 2 2 s 2 + t s ) + C k ( k + 1 ) ( 1 p 2 ) ( p + z ( x ) ) 2 + C j ( j 1 ) ( 1 q 2 ) ( q + z ( x ) ) 2 + C k p + z ( x ) p ( 2 t + 2 s + 2 k + 1 ) 2 t + 2 s 2 j 1 p 2 q p + C j q + z ( x ) q ( 2 t + 2 s 2 j + 1 ) + 2 t 2 s + 2 k 1 q 2 q p .
Note that the six terms of (26) and (27) contain terms with structure similar to the two and four terms appearing in (18) and (21), respectively, which allows us to link the parameters t, s, k, j, p and q with the α and β parameters of the Jacobi and the X 1 -type exceptional Jacobi polynomials.
Note also that if we choose V ( x ) and V + ( x ) to be free from an additive constant, e.g., in Equations (A17) and (A28), then the factorization energy is found to be
ϵ = C ( t + s + k j ) 2 .
In what follows, we focus on special cases that recover exactly solvable potentials.

3.2. Special Case: j = 0 , k = 0

For k = j = 0 , only the first two terms remain in Equations (26) and (27). Matching the coefficients of (26) with those of (18), one finds that
α + β 2 2 + α β 2 2 1 4 = 2 t 2 + 2 s 2 t s 2 α + β 2 α β 2 = 2 t 2 2 s 2 t + s .
Subtracting and adding the two equations, one obtains
α + 1 2 α 1 2 = 2 t ( 2 t 1 ) β + 1 2 β 1 2 = 2 s ( 2 s 1 ) .
There are four different solutions to this system of equations, corresponding to four distinct SUSY transformations. The actual values of t and s also determine the SUSY partner potential V + ( x ) . The results are summarized in Table 2 for j = k = 0 , where the results of the four different transformations are displayed. The table also contains the quantity E 0 ( ) ϵ , where
ϵ = C ( t + s ) 2
holds in this case. When this difference is zero, SUSY is unbroken, and we have the T 1 transformation from Table 1, while when it is negative, SUSY is broken and the T 2 , T 3 and T 4 transformations may emerge, depending on the behavior of χ ( x ) at the boundaries. These latter findings can be analyzed only when the allowed z ( x ) functions, and thus the potentials, are specified.

3.3. Special Case: j = 0 , k 0

Now, two more terms appear both in V ( x ) in (26) and V + ( x ) in (27). One of these can be cancelled in the former with k = 1 , while to cancel the second in order to recover the general PI class potential (18), a specific choice of p has to be made:
p = 2 t 2 s 2 t + 2 s + 1 .
Then one obtains
V ( x ) = V ( α , β , x )
under the same conditions as in Section 3.2, i.e., for t = ± α / 2 + 1 / 4 and s = ± β / 2 + 1 / 4 . However, further conditions follow from matching V + ( x ) to V ^ ( α , β , x ) (in particular, to its last two terms) in Equation (21). These conditions rule out two possibilities in which α and β have the same sign in t and s. The two remaining transformations are displayed in Table 2 for j = 0 , k = 1 . The factorization energy is now
ϵ = C ( t + s + 1 ) 2

3.4. Special Case: j 0 , k = 0

This case is the inverse of that discussed in Section 3.3 in the sense that now V + ( x ) in (27) will play the role of V ( α , β , x ) in (18), while V ( x ) has to be matched with V ^ ( α , β , x ) in (21). For this, the choice
q = 2 t 2 s 2 t + 2 s 1
has to be made, while one finds that again, the possibilities with the same sign of α and β in t and s have to be ruled out. With the actual value of the factorization energy
ϵ = C ( t + s 1 ) 2 ,
the results are summarized in Table 2 for j = 1 , k = 0 .

3.5. Special Case: j 0 , k 0

Now the six terms in Equations (26) and (27) do not vanish in general. However, for j = 1 and k = 1 , they can be both reduced to five terms, while with the appropriate choice of p and q, both potentials can be reduced to four terms, and it becomes possible to match them with the rationally extended version of the general PI class potential (21). Again, all four possibilities t = ± α / 2 + 1 / 4 , s = ± β / 2 + 1 / 4 follow from matching the first two terms, corresponding to the general PI class potential, while the elimination of one term each by the choice of p and q rules out the possibilities with different signs of α and β in t and s, respectively. The supersymmetric partner of V ( x ) = V ^ ( α , β , x ) can be determined in both cases from the actual values of t and s. The results are summarized in Table 2 for j = k = 1 , together with the consequences of the actual value of the factorization energy
ϵ = C ( t + s ) 2 .

4. Specializing to the Three PI Class Potentials

In this section, we substitute the three z ( x ) functions into the general formulae and inspect the prescriptions imposed by the boundary conditions on the parameters α and β . Further conditions follow for the rationally extended potentials V ^ ( α , β , x ) from the requirement that the denominator of the last two potential terms not have zeros within the domain of definition of z ( x ) , i.e.,
( α β ) z ( x ) + ( α + β ) 0 for x ( x , x + ) .
All these prescriptions determine t and s, which, in turn, decide which of the four transformations in Table 1 the given parameter combinations correspond to. A summary of these conditions is presented in Table 3.
In order to simplify the notation, in what follows we introduce the expressions
A ( α , β ) α + β 2 2 + α β 2 2 1 4 , B ( α , β ) α + β 2 α β 2 .

4.1. z ( x ) = cosh ( a x ) : The Generalized Pöschl–Teller Potential and Its Rational Extension

For this choice, the actual PI class potential is the generalized Pöschl–Teller potential, which is
V G ( α , β , x ) = a 2 sinh 2 ( a x ) A ( α , β ) + 2 a 2 cosh ( a x ) sinh 2 ( x ) B ( α , β )
with bound-state eigenfunctions
ψ n ( α , β ; x ) = C n ( α , β ) ( cosh ( a x ) 1 ) α 2 + 1 4 ( cosh ( a x ) + 1 ) β 2 + 1 4 P n ( α , β ) ( cosh ( x ) )
and energy eigenvalues
E n = a 2 n + α + β + 1 2 2 .
For the rationally extended generalized Pöschl–Teller potential, the corresponding formulae are as follows:
V ^ G ( α , β , x ) = V G ( α , β , x ) 2 a 2 ( α + β ) ( α β ) cosh ( a x ) + α + β 2 a 2 ( α + β ) 2 ( α β ) 2 ( α β ) cosh ( a x ) + α + β 2
with bound-state eigenfunctions
ψ ^ n ( α , β ; x ) = C ^ n ( α , β ) ( cosh ( a x ) 1 ) α 2 + 1 4 ( cosh ( a x ) + 1 ) β 2 + 1 4 [ ( α β ) cosh ( a x ) + α + β ] 1 P ^ n ( α , β ) ( cosh ( x ) )
and energy eigenvalues
E ^ n = a 2 n 1 + α + β + 1 2 2 .
Here, n runs from n = 1 . These expressions recover the findings of Ref. [17] with α = A B 1 / 2 and β = A + B 1 / 2 , noting also that the roles of V ( x ) and V + ( x ) are reversed in that publication.
The behavior of both potentials at x = 0 is of the type x 2 , while that of the bound-state wave functions is x α + 1 / 2 . For this reason, the prescription α > 1 / 2 holds in general, while for transformations decreasing α by one unit, α > 1 / 2 has to hold in order to guarantee the proper behavior of the bound-state solutions of the partner potential. Asymptotically decaying solutions are obtained for α + β + 1 < 2 n for the generalized Pöschl–Teller potential and α + β + 1 < 2 ( n + 1 ) for its rationally extended version. In both cases, the condition for having at least one bound state is α + β + 1 < 0 . In summary, the conditions valid for the α and β parameters for the potential V G ( α , β , x ) are
± 1 2 < α
α + β + 1 < 0
These conditions also hold for the potential V ^ ( α , β , x ) , but there is also an extra condition to guarantee its non-singularity for finite positive x values:
1 < α + β α β .
Elementary calculations show that combining conditions (46) and (47) leads to 0 < ( α + β + 1 ) + 1 ± 1 < α β , which, in turn, leads to
1 < α
when combined with Equation (48). These conditions have to be considered for the generalized Pöschl–Teller potential V G ( α , β , x ) and its rational extension V ^ G ( α , β , x ) in the admissible transformation schemes, applying them to the actual parameter values.
The SUSY transformations summarized for the general PI class potential V G ( α , β , x ) and its rational extension V ^ G ( α , β , x ) are specified for the case of the generalized Pöschl–Teller potential in Table 4.
The first column indicates the relevant values of j and k, which set the actual form of the solution χ ( x ) in Equation (24) that generates the SUSY transformation. The second column displays V + ( x ) , the supersymmetric partner of V ( x ) , the role of which is played by V G ( α , β , x ) in (40) or V ^ G ( α , β , x ) (43), depending on j and k. In the third column, the lower boundaries of α are presented, originating from Equations (46) and (49), whenever applicable. When both equations are relevant, both lower values are presented. The fourth and fifth columns contain t and s, while the sixth one is simply t + s + k j , setting the asymptotic behavior of the χ ( x ) solution. The seventh and eighth columns contain the ratios that determine whether the potentials playing the role of V ( x ) and V + ( x ) satisfy the condition of being non-singular at finite x values. These appear only when these potentials are of the type V ^ G ( α , β , x ) . In the ninth column, the difference in the ground state energy of V ( x ) and the factorization energy ϵ is displayed. When this difference is zero, SUSY is unbroken, and the relevant transformation is of the type T 1 , while when it is positive, SUSY is broken and the transformation can be of the type T 2 , T 3 or T 4 . Finally, the last column presents the actual type of transformation determined by the sign of t and t + s + k j in the fourth and sixth column, respectively, and confirmed by the ninth column. The values appearing in the inequalities in the third and in the sixth to ninth columns are calculated based on the conditions formulated in Equations (46)–(49).
It is seen that the SUSY transformation systematically changes α and β by one unit and connects the same type of potentials (for j = k ) or different types (for j k ). The T 1 - and T 2 -type transformations confirm the shape-invariance of the generalized Pöschl–Teller potential V G ( α , β , x ) (for j = k = 0 ) and of its rational extension V ^ G ( α , β , x ) (for j = k = 1 ).

4.2. z ( x ) = cos ( a x ) : The Scarf I Potential and Its Rational Extension

Now, the essential formulae for the Scarf I potential are the following:
V I ( α , β , x ) = a 2 sin 2 ( a x ) A ( α , β ) + 2 a 2 cos ( a x ) sin 2 ( x ) B ( α , β ) ,
ψ n ( α , β ; x ) = C n ( α , β ) ( 1 cos ( a x ) ) α 2 + 1 4 ( 1 + cos ( a x ) ) β 2 + 1 4 P n ( α , β ) ( cos ( x ) ) ,
E n = a 2 n + α + β + 1 2 2 .
Note that now the energy egenvalues are in the positive domain.
The corresponding results for the rationally extended Scarf I potential are
V ^ I ( α , β , x ) = V I ( α , β , x ) + 2 a 2 ( α + β ) ( α β ) cos ( a x ) + α + β 2 a 2 ( α + β ) 2 ( α β ) 2 ( α β ) cos ( a x ) + α + β 2 ,
ψ ^ n ( α , β ; x ) = C ^ n ( α , β ) ( 1 cos ( a x ) 1 ) α 2 + 1 4 ( 1 + cos ( a x ) ) β 2 + 1 4 [ α β ) cos ( a x ) + α + β ] 1 P ^ n ( α , β ) ( cos ( x ) ) ,
E ^ n = a 2 n 1 + α + β + 1 2 2 .
Again, n runs from n = 1 in accordance with the findings of Ref. [16], which correspond to the choice α = A B 1 / 2 and β = A + B 1 / 2 .
In spite of the mathematical similarities to the generalized Pöschl–Teller potential and its rational extension, the boundary conditions differ significantly. Now both boundaries occur at finite real values, where the singularity of the potential is of the type ξ 2 , where ξ = x and π / a x . The strength of these singularities are set independently by α and β . The SUSY transformations summarized in Table 2 for the general case either increase or decrease α and β by one unit, which set the minimal values of these parameters to 1 / 2 or 1 / 2 , respectively, for the given transformation. The actual requirements for α and β set by the boundary conditions of V I ( α , β , x ) (and V ^ I ( α , β , x ) ) are now
± 1 2 < α
± 1 2 < β
depending of whether the actual parameter value is increased (−) or decreased (+) by one unit. Besides these prescriptions set by the boundary conditions, the non-singularity requirement of V ^ I ( α , β , x ) introduces further conditions. These are now
α + β α β < 1
or
1 < α + β α β
depending on whether α < β or β < α holds. Note that, since α + β > 0 holds in most cases, the sign of the quantities in Equations (58) and (59) is determined by that of α β .
Table 5 summarizes the SUSY transformations specified for the Scarf I potential V I ( α , β , x ) and its rational extension V ^ I ( α , β , x ) . The notation in Table 5 is essentially the same as in Table 4, with the difference in the sixth and seventh columns, where the quantity setting the behavior of χ ( x ) at the upper boundary ( x + = π / a here and x + in Table 4) is displayed.
The general pattern of the SUSY transformations is similar to that observed for the generalized Pöschl–Teller potential and its rational extension. The shape-invariance of the Scarf I potential V I ( α , β , x ) and that of its rational extension V ^ I ( α , β , x ) is confirmed by the T 1 - and T 2 -type transformations with j = k = 0 and j = k = 1 , respectively.

4.3. z ( x ) = i sinh ( a x ) : The Scarf II Potential—The Real Version and Its Rational Extension

The general form of the Scarf II potential is obtained from an imaginary z ( x ) function; nevertheless, real potentials with real energy spectra follow from the formulae in a straightforward way. The potential obtained from the general PI class expressions is
V I I ( α , β , x ) = a 2 cosh 2 ( a x ) A ( α , β ) 2 i a 2 sinh ( a x ) cosh 2 ( x ) B ( α , β ) .
The bound-state eigenfunctions are
ψ n ( α , β ; x ) = C n ( α , β ) ( 1 i sinh ( a x ) ) α 2 + 1 4 ( 1 + i sinh ( a x ) ) β 2 + 1 4 P n ( α , β ) ( i sinh ( x ) ) ,
and the energy eigenvalues are
E n = a 2 n + α + β + 1 2 2 ,
which are formally the same as those obtained for the generalized Pöschl–Teller potential (42).
This potential is real for α * = β . (Note that in this setting, the wave functions (61) can also be expressed in terms of Romanowski polynomials, which can be obtained from Jacobi polynomials with imaginary arguments [34]. See also the more recent Ref. [33].) In this case, ( α + β ) / 2 = α R and ( α β ) / 2 = i α I . As another similarity to the generalized Pöschl–Teller potential, the regularity of the wave functions for x follows for
α + β + 1 < 0 ,
i.e.,
α R < 1 2 .
However, now this is also the condition for the regularity at x , and no further conditions arise for the parameters α and β .
Another novelty is that the rational extension of the real Scarf II potential cannot be defined, because the ( α β ) i sinh ( a x ) + α + β = 2 α I sinh ( a x ) + 2 α R function will take on a zero at a finite x for sinh ( a x ) = α R / α I . For this reason, we discuss only the real Scarf II potential here, so Table 6 contains only the case with j = k = 0 .
The first line in Table 6 describes a standard T 1 -type SUSYQM transformation associated with unbroken SUSY: it eliminates the state of V + ( x ) corresponding to the ground state of V ( x ) . The T 2 -type transformation in the second line is its inverse: it introduces a new ground state of V + ( x ) at the energy corresponding to the factorization energy ϵ . The remaining two transformations are, in principle, T 2 -type, because the χ ( x ) function that generates the SUSYQM transformation is unbound for x ± . However, the energy eigenvalues are complex (actually, the two energies are complex conjugates of each other), so they cannot be interpreted as valid transformations. These energy eigenvalues fit into the sequence of resonances that correspond to the irregular solutions of the real Scarf II potential, i.e., its resonances. Note that these can be obtained by the α α replacement, which leave potential (60) invariant, but converts the energy eigenvalues (62) into complex.
Before closing this section, we note that a much richer network of SUSYQM transformations can be generated for the PT -symmetric Scarf II potential. (For a recent review of PT -symmetric quantum mechanics, see the monograph [35]). PT -symmetric potentials are special complex potentials that are invariant under the simultaneous space and time reflection, the latter represented by complex conjugation. In practical terms, for one-dimensional potentials PT , symmetry corresponds to the requirement V * ( x ) = V ( x ) . The PT -symmetric version of the Scarf II potential corresponds to real values of α and β for unbroken PT symmetry. With this choice, the denominators appearing in the formulae will be defined on the complex plane and avoid zero, cancelling any singularities along the real x axis. As another specific feature, the discrete spectrum of these potentials is real (for unbroken PT symmetry), which is secured by the real values of α and β . Furthermore, in many PT -symmetric potentials, the bound-state (normalizable) solutions belong to two series discriminated by the quasi-parity quantum number Q, which takes on the values + 1 and 1 or + and − in short. In the case of the PT -symmetric Scarf II potential, this dual system is implemented by replacing α by Q α . This means that the PT -symmetric Scarf II potential has two sets of solutions, including two ground states, and consequently, two SUSYQM partner potentials discriminated by Q [36]. This result can be reproduced by adapting the four transformations for j = k = 0 in Table 2. Further findings will be discussed in detail elsewhere, leaving the present work focusing on real potentials only. For the general complex Scarf II potential, see, e.g., the recent study [37].

5. Summary and Conclusions

We presented a formalism in which the network of SUSYQM transformations between some shape-invariant potentials belonging to the PI class and their rationally extended versions can be discussed in a unified framework. We specified this formalism for the three members of this class, the generalized Pöschl–Teller, the Scarf I and the Scarf II potentials. The common element behind these potentials that allowed their unified treatment is that the variable transformation that maps the Schrödinger equation into the differential equation of the Jacobi polynomials and that of the X 1 -type Jacobi polynomials satisfies the same differential equation.
The formalism was based on a general six-parameter function χ ( x ) , from which the W ( x ) superpotential and the supersymmetric partner potentials V ( x ) and V + ( x ) could be constructed. With a careful choice of parameters, V ( x ) could be brought to the form of the general PI class potential or its rational extension. The choice of parameters also determined the behavior of χ ( x ) at the boundaries, selecting one of the four possible types of SUSYQM transformations, corresponding to unbroken or broken supersymmetry. This construction led to a network of SUSYQM transformations within the PI class, between the PI class potential and its rationally extended version, as well as within the latter potential class.
We specified the general results for the members of the PI class and established the same network of SUSY transformations for the generalized Pöschl–Teller, Scarf I and Scarf II potentials, as well as for their rational extensions. At this point, the difference between the dissimilar domain of definition of the three potentials and the boundary conditions implied by them resulted in differences in the results, too. It was found that, although the generalized Pöschl–Teller and the Scarf I potential are defined on the positive real axis and on a finite domain, respectively, due to the fact that their mathematical formalism is based on related functions (hyperbolic and trigonometric), the general patterns of the basic results are rather similar. In contrast to these findings, the real Scarf II potential differed significantly from the other two: its rational extension could not be constructed, because no parameter sets could be chosen to guarantee that the potential is singularity-free on the real coordinate axis. In this case, SUSYQM transformations are restricted to the Scarf II potentials differing in their potential parameters. However, without presenting the details, we noted that the PT -symmetric version of the Scarf II potential and its rational extension are associated with a rather rich network of SUSYQM transformations.
The present formalism was suitable to place the rationally extended PI class potentials into a unified framework of exactly solvable potentials. Furthermore, the fact that they fulfill the requirement of shape-invariance could be proven in a straightforward way. Their place beyond the Natanzon class could also be verified: the X 1 -type Jacobi polynomials, which appear in the bound-state solutions, represent a special polynomial solution of the confluent Heun differential equation, functions that also contain the hypergeometric function, and thus, the Jacobi polynomials as special cases.
The procedure presented here is analogous to a recent one devised to discuss the relation of the radial harmonic oscillator and its rationally extended version [23]. There, the exceptional polynomial was the X 1 -type Laguerre polynomial. A similar network of SUSYQM transformations between the named potentials was identified. The difference between that study and the present one is that the PI class potentials could be discussed on an abstract unified level, as there are several (three) members of a potential class, while the radial harmonic oscillator stands alone in its potential class, named LI [5].
The formalism presented here and in Ref. [23] also explains why some potentials have rational extensions, in contrast to others. The explanation lies in the key system of Equations (A18)–(A20) for the Jacobi polynomials and (A29)–(A33) for the X 1 -type Jacobi polynomials (as well as the analogous set of equations in Ref. [23] for the X 1 -type Laguerre polynomials), which have to be solved in order to obtain exactly solvable potentials. In all these cases, only one equation depends on n, the degree of the polynomial, which discriminates the individual energy levels. Setting only the corresponding p i parameter to a non-zero value implies that the energy eigenvalues E n are directly related to n, while the potential parameter s i , as well as α and β , will be free from any dependence on n. In any other case, α and β necessarily depend on n, which complicates the formalism beyond the point that would allow n-independent (i.e., state-independent) potentials.
It has to be added that some of these results have been known in the literature in some form. The findings were usually isolated and typically covered only certain aspects of the problem. They also focused on different contexts and often used different formalisms. Reference [17] is an example where all three potentials were discussed in some detail, and some arguments exposing analogies between them were also formulated. We hope that the unified treatment of a class of potentials in a general context will help to arrange the fragments into a systematic and coherent framework.
The generalization of these results can also be envisaged. Here, the χ ( x ) function contained expressions ( p + z ( x ) ) corresponding to a first-order polynomial of z. This construction is based on potentials solved in terms of X 1 -type Laguerre polynomials [14,15]. It would be a straightforward choice to consider χ ( x ) functions with higher (nodeless) polynomial forms. One can expect that this could be possible considering potentials generated from X m - rather than X 1 -type Jacobi and Laguerre polynomials [38].

Author Contributions

Conceptualization, G.L.; Methodology, T.S.; Investigation, T.S. and L.F.P.; Writing—original draft, G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research, Development and Innovation Office (NKFIH, Hungary), Grant No. K128729.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Here, we adapt a well-known transformation method [4,5,32] to the Jacobi and X 1 -type exceptional Jacobi polynomials to generate the general PI class potentials and their rational extension. This transformation method maps the second-order differential equation of some F ( z ) special function
d 2 F d z 2 + Q ( z ) d F d z + R ( z ) F ( z ) = 0 ,
into the one-dimesional Schrödinger equation
d 2 ψ d x 2 + ( E V ( x ) ) ψ ( x ) = 0 ,
by applying the variable transformation z ( x ) . (For the sake of simplicity, the units 2 m = 1 and = 1 are used.) The domain of definition of F ( z ) is not specified at this stage, while the Schrödinger equation is assumed to be defined on the positive real semi-axis ( x [ 0 , ) , for radial potentials), the full real axis ( x ( , ) ) or a finite domain ( x [ x , x + ] ) of it. The procedure is applicable to each of these situations. We assume that the solution of the Schrödinger equation is written as
ψ ( x ) = f ( x ) F ( z ( x ) ) ,
Substituting this into Equation (A2) and comparing the appropriate terms, one finds that
E V ( x ) = z ( x ) 2 z ( x ) 3 4 z ( x ) z ( x ) 2 + ( z z ( x ) ) 2 R ( z ( x ) ) 1 2 d Q d z 1 4 Q 2 ( z ( x ) ) .
Furthermore, the solutions of the Schrödinger equation can be expressed as
ψ ( x ) ( z ( x ) ) 1 2 exp 1 2 z ( x ) Q ( z ) d z F ( z ( x ) ) .
Up to this point, Q ( z ) and R ( z ) are defined by the choice of the special function F ( z ) , but the z ( x ) function is unspecified. z ( x ) defines the variable transformation, so one would expect that it satisfies certain criteria, e.g., that it is monontonous and single-valued. However, not any z ( x ) function leads to an exactly solvable potential V ( x ) . This is because nothing guarantees that one finds the solutions of the same  V ( x ) potential for all energy eigenvalues E. This problem was handled by the proposal outlined in Ref. [4] after noticing that the constant E on the left-hand side of Equation (A4) should have a correspondent on the right-hand side among the various terms appearing there. Equating one (or more) term with a constant leads to a differential equation for z ( x ) . Furthermore, it is advisable to chose from the terms originating from the Q ( z ) and R ( z ) functions, because in this case, the differential equation will contain only z ( x ) , the first derivative of z. Furthermore, these terms contain the parameters appearing in the special function F ( z ) . The two terms containing higher derivatives depend only on the parameters that appear in the z ( x ) function, so they are less suitable to generate the E energy eigenvalues. These two terms together are called the Schwarzian derivative
{ z , x } z ( x ) z ( x ) 3 2 z ( x ) z ( x ) 2 .
Let us prescribe the differential equation
d z d x 2 Φ ( z ) = C .
Its solution can be found by direct integration
Φ 1 / 2 ( z ) d z = C 1 / 2 x + x 0 ,
although one obtains the x ( z ) function this way, which is the inverse of z ( x ) . Nothing guarantees that one can also obtain z ( x ) in closed form from x ( z ) . Nevertheless, both V ( x ) and the solutions (A5) can be generated without any complication. In this case, the potential is called implicit. For the most well-known potentials, however, it is possible to obtain z ( x ) in an explicit form. Note that an integration constant, x 0 , appears in Equation (A8). It is a simple coordinate shift that usually does not play an important role. It has no effect on the energy spectrum, but it can be used to set z ( 0 ) = 0 , for example. However, in PT -symmetric quantum mechanics, it can be chosen to be imaginary, which allows one to avoid singularities appearing on the real x line.
Combining Equations (A4) and (A7), z ( x ) can be replaced in the former expression, leaving V ( x ) and E as the function of Q ( z ) , R ( z ) , Φ ( z ) , z ( x ) and C:
E V ( x ) = 1 2 { z , x } + C Φ ( z ( x ) ) R ( z ( x ) ) 1 2 d Q d z 1 4 Q 2 ( z ( x ) ) .
when F ( z ) is chosen to be the hypergeometric function F 1 2 ( a , b ; c ; z ) , one obtains the Natanzon-class potentials [1]. The bound-state solutions of these potentials are written in terms of Jacobi polynomials, which are obtained from the hypergeometric function for a = n or b = n . (Similarly, Natanzon confluent potentials [2] are obtained in an analogous way from the confluent hypergeometric function F 1 1 ( a ; c ; z ) , which reduces to generalized Laguerre polynomials for a = n .) Considering F ( z ) functions that are different from or more general than the (confluent) hypergeometric function naturally leads to potentials beyond the Natanzon class.
Specifying the procedure for the Jacobi polynomial [32] F ( z ) = P n ( α , β ) ( z ) corresponds to the choice [3]
Q ( z ) = ( β + α + 2 ) z + ( α β ) 1 z 2 ,
R ( z ) = n ( n + α + β + 1 ) 1 z 2 .
With these choices and introducing for simplicity the definitions
ω = α + β 2 , ρ = α β 2 ,
the actual form of Equation (A4) becomes
E V ( x ) = 1 2 { z , x } + ( z ( x ) ) 2 ( 1 z 2 ( x ) ) 2 n + ω + 1 2 2 1 4 ( 1 z 2 ( x ) ) + ( 1 ω 2 ρ 2 ) 2 ω ρ z ( x ) ,
while Equation (A7) turns into
d z d x 2 Φ ( z ( x ) ) d z d x 2 φ ( z ( x ) ) ( z 2 ( x ) 1 ) 2 = C ,
i.e.,
d z d x = C 1 / 2 ( z 2 ( x ) 1 ) [ φ ( z ( x ) ) ] 1 / 2 ,
where
φ ( z ( x ) ) = p I ( 1 z 2 ( x ) ) + p II + p III z ( x ) .
The x ( z ) function, and, whenever possible, its inverse, z ( x ) is then determined from the integration of Equation (A15).
From the structure of (A13), one can conclude that the most general form of the derived potential is
V ( x ) = 1 2 { z , x } + C φ ( z ( x ) ) s I ( 1 z 2 ( x ) ) + s II + s III z ( x ) .
Note that φ ( z ( x ) ) in Equation (A16) and the expressions in square brackets in Equations (A13) and (A17) contain the same linearly independent terms.
Substituting (A17) and ( z ( x ) ) 2 from (A14) into Equation (A13) and comparing the corresponding terms, one finds that the following three equations have to be satisfied simultaneously:
n + ω + 1 2 2 1 4 + s I p I E C = 0 ,
1 ω 2 ρ 2 + s II p II E C = 0 ,
2 ω ρ + s III p III E C = 0 ,
This set of equations points out the intimate and subtle relation that connects the parameters p i appearing in the z ( x ) function, the coupling coefficients ( s i ) of the potential (A17), the energy eigenvalue E and the parameters of the Jacobi polynomials (see Equation (A1) with Equations (A10) and (A11)).
The PI class potentials are obtained by setting p I = 1 and p II = p III = 0 , corresponding to φ ( z ) = 1 z 2 . Note that this choice secures that E and n appear only in Equation (A18), so the dependence of E n on n directly follows from this equation:
E n = C n + ω + 1 2 2 = C n + α + β + 1 2 2
if we set the energy scale such that s I = 1 4 . The remaining two equations directly set the coupling coefficients as s II = ω 2 + ρ 2 1 and s III = 2 ω ρ . The general form of the PI class potential can be calculated after evaluating the contribution of the Schwarzian derivative (A6), which is found to be 3 C / ( 1 z 2 ( x ) ) C / 4 :
V ( α , β , x ) = C 1 z 2 ( x ) α + β 2 2 + α β 2 2 1 4 + 2 C z ( x ) 1 z 2 ( x ) + α + β 2 α β 2 .
The bound-state wave functions are
ψ n ( α , β ; x ) = C n ( α , β ) ( 1 z ( x ) ) α 2 + 1 4 ( 1 + z ( x ) ) β 2 + 1 4 P n ( α , β ) ( z ( x ) ) .
Other choices of φ ( z ) lead to further classes of potentials. Taking only p II = 1 to be non-zero, the PII potential class is obtained [5], which contains three potentials for three different z ( x ) functions: the Rosen–Morse I and II and the Eckart potentials, while from p III = 1 , the PIII potential follows [39], which is an implicit potential, meaning that x ( z ) cannot be inverted into an explicit z ( x ) function. Nevertheless, all the relevant quantities can be calculated exactly. Further potentials can be obtained by choosing two or three different p i parameters to be non-zero. Examples are the Ginocchio [40], generalized Ginocchio [41] and DKV [42] potentials and a potential from which both the PI and PII class potentials can be obtained as special cases [43]. For a complete list of these potentials, see Ref. [32]. Alternative approaches to potentials [40,43] can be found in Ref. [44], while the relation of these studies is discussed in Ref. [45].
The analogous procedure for the X 1 -type exceptional Jacobi polynomials F ( z ) = P ^ n ( α , β ) ( z ) follows from the choice [14,15]
Q ( z ) = ( β + α ) z ( β α ) z 2 1 ( β α ) ( β α ) z ( β + α ) ,
R ( z ) = ( β α ) z ( n 1 ) ( n + α + β ) z 2 1 ( β α ) 2 ( β α ) z ( β + α ) .
With the notation in Equation (A12), Equation (A4) becomes
E V ( x ) = z ( x ) 2 z ( x ) 3 4 z ( x ) z ( x ) 2 + ( z ( x ) ) 2 ( 1 z 2 ( x ) ) 2 n + ω 1 2 2 1 4 ( 1 z 2 ( x ) ) + ( 1 ω 2 ρ 2 ) 2 ω ρ z ( x ) + 2 ω z 2 ( x ) 1 ρ z ( x ) + ω + 2 ( ρ 2 ω 2 ) z 2 ( x ) 1 ( ρ z ( x ) + ω ) 2 .
Now, the φ ( z ) function contains two more linearly independent terms compared to that in Equation (A16):
φ ( z ( x ) ) = p I ( 1 z 2 ( x ) ) + p II + p III z ( x ) + p IV z 2 ( x ) 1 ρ z ( x ) + ω + p V ( z 2 ( x ) 1 ) 2 ( ρ z ( x ) + ω ) 2 .
Two new terms also appear in the potential
V ^ ( x ) = 1 2 { z , x } + C φ ( z ( x ) ) s I ( 1 z 2 ( x ) ) + s II + s III z ( x ) + s IV ( z 2 ( x ) 1 ) ρ z ( x ) + ω + s V ( z 2 ( x ) 1 ) ( ρ z ( x ) + ω ) 2 .
Substituting (A28) and ( z ( x ) ) 2 into Equation (A26), the set of key equations is enlarged:
n 1 + ω + 1 2 2 1 4 + s I p I E C = 0 ,
1 ω 2 ρ 2 + s II p II E C = 0 ,
2 ω ρ + s III p III E C = 0 ,
2 ω + s IV p IV E C = 0 ,
2 ( ρ 2 ω 2 ) + s V p V E C = 0 .
Note that the first three of these equations coincide with Equations (A18)–(A20), with the difference that in (A29) n 1 appears rather than n. This corresponds to the fact that the X 1 -type exceptional Jacobi polynomials have n 1 zeros within the domain of integration. According to this, the energy eigenvalues of this potential class are written as
E ^ n = C n 1 + ω + 1 2 2 = C n 1 + α + β + 1 2 2 ,
where n = 1 , 2, …. The ground-state energy is E ^ 1 . The corresponding equations for the potential and the bound-state wave functions are, respectively,
V ^ ( α , β , x ) = C 1 z 2 ( x ) α + β 2 2 + α β 2 2 1 4 + 2 C z ( x ) 1 z 2 ( x ) α + β 2 α β 2 + 2 C ( α + β ) ( α β ) z ( x ) + α + β + 2 C ( α β ) 2 ( α + β ) 2 ( α β ) z ( x ) + α + β 2 ,
ψ n ( α , β ; x ) = C ^ n ( α , β ) ( 1 z ( x ) ) α 2 + 1 4 ( 1 + z ( x ) ) β 2 + 1 4 ( α β ) z ( x ) + α + β ) 1 P ^ n ( α , β ) ( z ( x ) ) .
As discussed above, the P ^ n ( α , β ) ( z ) X 1 -type exceptional Jacobi polynomials [14,15], similarly to the X 1 -type exceptional Laguerre polynomials, bear most characteristics of their classical counterparts except that the degree of the lowest polynomial is larger than zero: in the case of X 1 -type exceptional orthogonal polynomials, the sequence starts with 1. This is because one of their zeros falls outside the domain of integration. It is thus customary to use ν = n 1 as their index: in this case, ν = 0 , 1 , and ν indicates the degree of the polynomial within the domain of integration. The X 1 -type exceptional Jacobi polynomials can be expressed in terms of two classical Jacobi polynomials as
P ^ n ( α , β ) ( z ) = ( α β ) z + α + β 2 ( α β ) α + β ( α β ) ( α + β + 2 n 2 ) P n 1 ( α , β ) ( z ) + 1 α + β + 2 n 2 P n 2 ( α , β ) ( z )
The X 1 -type exceptional orthogonal polynomials have an even wider class, the X m -type [46,47,48,49], which contains the X 1 -type polynomials as a special case. Based on these, a more general rationally extended version of PI class potentials can be constructed; however, we do not consider this option here.

References

  1. Natanzon, G.A. General properties of potentials for which the Schrödinger equation can be solved by means of hypergeometric functions. Teor. Mat. Fiz. 1979, 38, 146–153. [Google Scholar] [CrossRef]
  2. Cordero, P.; Salamó, S. Algebraic solution for the Natanzon confluent potentials. J. Phys. A Math. Gen. 1991, 24, 5299–5305. [Google Scholar] [CrossRef]
  3. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions; Dover: New York, NY, USA, 1970. [Google Scholar]
  4. Bhattacharjie, A.; Sudarshan, E.C.G. A class of solvable potentials. Nuovo Cim. 1962, 25, 864–879. [Google Scholar] [CrossRef]
  5. Lévai, G. A search for shape-invariant solvable potentials. J. Phys. A Math. Gen. 1989, 22, 689–702. [Google Scholar] [CrossRef]
  6. Witten, E. Dynamical Breaking of Supersymmetry. Nucl. Phys. B 1981, 188, 513. [Google Scholar] [CrossRef]
  7. Cooper, F.; Khare, A.; Sukhatme, U. Supersymmetry and Quantum Mechanics. Phys. Rep. 1995, 251, 267. [Google Scholar] [CrossRef]
  8. Junker, G. Supersymmetric Methods in Quantum and Statistical Physics; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  9. Darboux, G. Sur une proposition relative aux équations linéaires. C. R. Acad. Sci. 1982, 94, 1456. [Google Scholar]
  10. Gendenshtein, L.E. Derivation of exact spectra of the Schrödinger equation by means of supersymmetry. JETP Lett. 1983, 38, 356–359. [Google Scholar]
  11. Infeld, L.; Hull, T.D. The Factorization Method. Rev. Mod. Phys. 1951, 23, 21. [Google Scholar] [CrossRef]
  12. Schrödinger, E. A Method of Determining Quantum-Mechanical Eigenvalues and Eigenfunctions. Proc. R. Ir. Acad. A 1940, 46, 9. [Google Scholar]
  13. Schrödinger, E. The factorization of the hypergeometric equation. Proc. R. Ir. Acad. A 1941/1942, 47, 53–54. [Google Scholar]
  14. Gómez-Ullate, D.; Kamran, N.; Milson, R. An extended class of orthogonal polynomials defined by a Sturm-Liouville problem. J. Math. Anal. Appl. 2009, 359, 352. [Google Scholar] [CrossRef]
  15. Gómez-Ullate, D.; Kamran, N.; Milson, R. An extension of Bochner’s problem: Exceptional invariant subspaces. J. Approx. Theor. 2010, 162, 987. [Google Scholar] [CrossRef]
  16. Quesne, C. Exceptional orthogonal polynomials, exactly solvable potentials and supersymmetry. J. Phys. A Math. Theor. 2008, 41, 392001. [Google Scholar] [CrossRef]
  17. Bagchi, B.; Quesne, C.; Roychoudhury, R. Isospectrality of conventional and new extended potentials, second-order supersymmetry and role of symmetry. Pramana J. Phys. 2009, 73, 337. [Google Scholar] [CrossRef]
  18. Midya, B.; Roy, B. Infinite families of (non)-Hermitian Hamiltonians associated with exceptional Xm Jacobi polynomials. J. Phys. A Math. Theor. 2013, 46, 175201. [Google Scholar] [CrossRef]
  19. Ronveaux, A. (Ed.) Heun’s Differential Equations; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  20. Takemura, K. Heun’s equation, generalized hypergeometric function and exceptional Jacobi polynomial. J. Phys. A Math. Theor. 2012, 45, 085211. [Google Scholar] [CrossRef]
  21. Lévai, G. Potentials from the Polynomial Solutions of the Confluent Heun Equation. Symmetry 2023, 15, 461. [Google Scholar] [CrossRef]
  22. Ishkhanyan, A.M. Schrödinger potentials solvable in terms of the general Heun functions. Ann. Phys. 2018, 388, 456–471. [Google Scholar] [CrossRef]
  23. Lévai, G.; Soltész, T. Unified supersymmetric transformations for the harmonic oscillator and its rational extension. Eur. J. Phys. 2020, 41, 025403. [Google Scholar] [CrossRef]
  24. Cooper, F.; Ginocchio, J.N.; Khare, A. Relationship between supersymmetry and solvable potentials. Phys. Rev. D 1987, 36, 2458. [Google Scholar] [CrossRef]
  25. Odake, S.; Sasaki, R. Krein-Adler transformations for shape-invariant potentials and pseudo virtual states. J. Phys. A 2013, 46, 245201. [Google Scholar] [CrossRef]
  26. Sukumar, C.V. Supersymmetric quantum mechanics and the inverse scattering method. J. Phys. A Math. Gen. 1985, 18, 2937. [Google Scholar] [CrossRef]
  27. Sukumar, C.V. Supersymmetric quantum mechanics of one dimensional systems. J. Phys. A Math. Gen. 1985, 18, 2917. [Google Scholar] [CrossRef]
  28. Ancarani, L.U.; Baye, D. Iterative supersymmetric construction of phase-equivalent potentials. Phys. Rev. A 1992, 46, 206. [Google Scholar] [CrossRef]
  29. Lévai, G.; Baye, D.; Sparenberg, J.-M. Analytical results on generating phase-equivalent potentials by supersymmetry: Removal and addition of bound states. Lect. Not. Phys. 1997, 488, 363. [Google Scholar]
  30. Baye, D.; Sparenberg, J.-M. Supersymmetric transformations of real potentials on the line. J. Phys. A Math. Gen. 1995, 28, 5079. [Google Scholar]
  31. Frank, W.M.; Land, D.J.; Spector, R.M. Singular potentials. Rev. Mod. Phys. 1971, 43, 36. [Google Scholar] [CrossRef]
  32. Lévai, G. PT Symmetry in Natanzon-class Potentials. Int. J. Theor. Phys. 2015, 54, 2724. [Google Scholar] [CrossRef]
  33. Quesne, C. Extending Romanowski polynomials in quantum mechanics. J. Math. Phys. 2013, 54, 122103. [Google Scholar] [CrossRef]
  34. Raposo, A.P.; Weber, H.J.; Alvarez-Castillo, D.E.; Kirchbach, M. Romanowski polynomials in selected physics models. Cent. Eur. J. Phys. 2007, 5, 253. [Google Scholar]
  35. Bender, C.M.; Dorey, P.E.; Dunning, C.; Fring, A.; Hook, D.W.; Jones, H.F.; Kuzhel, S.; Lévai, G.; Tateo, R. PT Symmetry in Quantum and Classical Physics; World Scientific Publishing Europe Ltd.: London, UK, 2019. [Google Scholar]
  36. Lévai, G.; Znojil, M. The interplay of supersymmetry and PT symmetry in quantum mechanics: A case study for the Scarf II potential. J. Phys. A Math. Gen. 2002, 35, 8793. [Google Scholar] [CrossRef]
  37. Acar, Y.C.; Acevedo, L.; Kuru, S. Unusual isospectral factorizations of shape invariant Hamiltonians with Scarf II potential. Phys. Scr. 2023, 98, 125229. [Google Scholar] [CrossRef]
  38. Yadav, R.K.; Kumari, N.; Khare, A.; Mandal, B.P. Group theoretic approach to rationally extended shape invariant potentials. Ann. Phys. 2015, 359, 46–54. [Google Scholar] [CrossRef]
  39. Lévai, G. A class of exactly solvable potentials related to the Jacobi polynomials. J. Phys. A Math. Gen. 1991, 24, 131. [Google Scholar] [CrossRef]
  40. Ginocchio, J.N. A class of exactly solvable potentials. I. One-dimensional Schrödinger equation. Ann. Phys. 1984, 152, 203. [Google Scholar] [CrossRef]
  41. Ginocchio, J.N. A class of exactly solvable potentials II. The three-dimensional Schrödinger equation. Ann. Phys. 1985, 159, 467. [Google Scholar] [CrossRef]
  42. Dutt, R.; Khare, A.; Varshni, Y.P. New class of conditionally exactly solvable potentials in quantum mechanics. J. Phys. A Math. Gen. 1995, 28, L107. [Google Scholar] [CrossRef]
  43. Lévai, G. Gradual spontaneous breakdown of PT symmetry in a solvable potential. J. Phys. A Math. Theor. 2012, 45, 444020. [Google Scholar] [CrossRef]
  44. Milson, R. Liouville transformation and exactly solvable Schrödinger equations. Int. J. Theor. Phys. 1998, 37, 1735. [Google Scholar] [CrossRef]
  45. Natanson, G. Overlapping of Lévai’s and Milson’s e-tangent-polynomial potentials. along symmetric curves. Axioms 2023, 12, 584. [Google Scholar] [CrossRef]
  46. Odake, S.; Sasaki, R. Infinitely many shape invariant potentials and new orthogonal polynomials. Phys. Lett. B 2009, 679, 414. [Google Scholar] [CrossRef]
  47. Odake, S.; Sasaki, R. Another set of infinitely many exceptional (Xl) Laguerre polynomials. Phys. Lett. B 2010, 684, 173. [Google Scholar] [CrossRef]
  48. Gómez-Ullate, D.; Kamran, N.; Milson, R. Exceptional orthogonal polynomials and the Darboux transformation. J. Phys. A Math. Theor. 2010, 43, 434016. [Google Scholar] [CrossRef]
  49. Gómez-Ullate, D.; Kamran, N.; Milson, R. On orthogonal polynomials spanning a non-standard flag. Contemp. Math. 2012, 563, 51. [Google Scholar]
Table 1. SUSYQM transformations of a potential problem defined on the domain x [ x , x + ] . These belong to four different types of solutions χ ( x ) .
Table 1. SUSYQM transformations of a potential problem defined on the domain x [ x , x + ] . These belong to four different types of solutions χ ( x ) .
Transformation T 1 T 2 T 3 T 4
ϵ ϵ = E 0 ( ) ϵ < E 0 ( ) ϵ < E 0 ( ) ϵ < E 0 ( )
lim x x χ convergentdivergentconvergentdivergent
lim x x + χ convergentdivergentdivergentconvergent
Spectrum modificationdeletes g.s.adds new g.s.nonenone
Table 2. SUSY transformations for the general form of PI class potentials. V ( x ) and V + ( x ) are expressed in terms of potentials (18) or (21), respectively.
Table 2. SUSY transformations for the general form of PI class potentials. V ( x ) and V + ( x ) are expressed in terms of potentials (18) or (21), respectively.
j, k V ( x ) tspq V + ( x ) ϵ E 0 ( ) ϵ
0, 0 V ( α , β , x )
α 2 + 1 4 β 2 + 1 4 V ( α + 1 , β + 1 , x ) C α + β + 1 2 2 0
α 2 + 1 4 β 2 + 1 4 V ( α 1 , β 1 , x ) C α β + 1 2 2 C ( α + β )
α 2 + 1 4 β 2 + 1 4 V ( α + 1 , β 1 , x ) C α β + 1 2 2 C ( α + 1 ) β
α 2 + 1 4 β 2 + 1 4 V ( α 1 , β + 1 , x ) C α + β + 1 2 2 C α ( β + 1 )
0, 1 V ( α , β , x )
α 2 + 1 4 β 2 + 1 4 α + β α β + 2 V ^ ( α + 1 , β 1 , x ) C α β + 3 2 2 C ( α + 2 ) ( β 1 )
α 2 + 1 4 β 2 + 1 4 α + β α β 2 V ^ ( α 1 , β + 1 , x ) C α + β + 3 2 2 C ( α 1 ) ( β + 2 )
1, 0 V ^ ( α , β , x )
α 2 + 1 4 β 2 + 1 4 α + β α β V ( α + 1 , β 1 , x ) C α β 1 2 2 C α ( β + 1 )
α 2 + 1 4 β 2 + 1 4 α + β α β V ( α 1 , β + 1 , x ) C α + β 1 2 2 C β ( α + 1 )
1, 1 V ^ ( α , β , x )
α 2 + 1 4 β 2 + 1 4 α + β + 2 α β α + β α β V ^ ( α + 1 , β + 1 , x ) C α + β + 1 2 2 0
α 2 + 1 4 β 2 + 1 4 α + β 2 α β α + β α β V ^ ( α 1 , β 1 , x ) C α β + 1 2 2 C ( α + β )
Table 3. The key values for the three potentials. The last column displays conditions for p and q.
Table 3. The key values for the three potentials. The last column displays conditions for p and q.
z ( x ) C x x + z z + p, q
cosh ( x ) a 2 01 1 < p , q
cos ( x ) a 2 0 π a 1 1 p , q < 1 , 1 < p , q
i sinh ( x ) a 2 n. a.
Table 4. SUSY transformations for the generalized Pöschl–Teller potential V G ( α , β , x ) (40) and its rational extension V ^ G ( α , β , x ) (43). These potentials play the role of V + ( x ) , depending on the parameters t, s, j and k, while V ( x ) = V G ( α , β , x ) for j = 1 and V ( x ) = V ^ G ( α , β , x ) for j = 1 . The type of transformation ( T i in Table 1 is determined by these parameters and the conditions prescribed for the potential parameters α and β in this subsection.)
Table 4. SUSY transformations for the generalized Pöschl–Teller potential V G ( α , β , x ) (40) and its rational extension V ^ G ( α , β , x ) (43). These potentials play the role of V + ( x ) , depending on the parameters t, s, j and k, while V ( x ) = V G ( α , β , x ) for j = 1 and V ( x ) = V ^ G ( α , β , x ) for j = 1 . The type of transformation ( T i in Table 1 is determined by these parameters and the conditions prescribed for the potential parameters α and β in this subsection.)
j, k V + ( x ) α
Condition
ts t + s + k j α + β α β
in V ( x )
α + β α β
in V + ( x )
E 0 ( ) ϵ T i
0, 0
V G ( α + 1 , β + 1 , x ) > 1 2 α 2 + 1 4 β 2 + 1 4 α + β + 1 2 < 0 0 T 1
V G ( α 1 , β 1 , x ) > 1 2 α 2 + 1 4 β 2 + 1 4 1 α + β + 1 2 > 0 a 2 ( α + β ) > 0 T 2
V G ( α + 1 , β 1 , x ) > 1 2 α 2 + 1 4 β 2 + 1 4 α + 1 α + β + 1 2 > 0 a 2 ( α + 1 ) ( β ) > 0 T 3
V G ( α 1 , β + 1 , x ) > 1 2 α 2 + 1 4 β 2 + 1 4 α + β + 1 2 α < 0 a 2 ( β + 1 ) ( α ) > 0 T 4
0, 1
V ^ G ( α + 1 , β 1 , x ) > 1 2 α 2 + 1 4 β 2 + 1 4 α + 1 α + β + 1 2 > 0 α + β α β + 2 > 1 a 2 ( α + 2 ) ( 1 β ) T 3
V ^ G ( α 1 , β + 1 , x ) > 1 > 1 2 α 2 + 1 4 β 2 + 1 4 α + β + 1 2 α < 0 α + β α β 2 > 1 a 2 ( 1 α ) ( β + 2 ) > 0 T 4
1, 0
V G ( α + 1 , β 1 , x ) > 0 > 1 2 α 2 + 1 4 β 2 + 1 4 α α + β + 1 2 > 0 α + β α β > 1 a 2 ( α ) ( β + 1 ) T 3
V ^ G ( α 1 , β + 1 , x ) > 1 > 1 2 α 2 + 1 4 β 2 + 1 4 α + β + 1 2 α < 0 α + β α β > 1 a 2 ( α + 1 ) ( β ) > 0 T 4
1, 1
V G ( α + 1 , β + 1 , x ) > 0 > 1 2 α 2 + 1 4 β 2 + 1 4 α + β + 1 2 < 0 α + β α β > 1 α + β + 2 α β > 1 0 T 1
V ^ G ( α 1 , β 1 , x ) > 1 > 1 2 α 2 + 1 4 β 2 + 1 4 1 α + β + 1 2 > 0 α + β α β > 1 α + β 2 α β > 1 a 2 ( α + β ) > 0 T 2
Table 5. The same as Table 4 for the Scarf I potential V I ( α , β , x ) (50) and its rational extension I ^ I ( α , β , x ) (53). In the seventh and eighth columns, 1 < ( A + B ) / ( A B ) < 1 indicates that the first inequality holds for A > B , and the second one for A < B .
Table 5. The same as Table 4 for the Scarf I potential V I ( α , β , x ) (50) and its rational extension I ^ I ( α , β , x ) (53). In the seventh and eighth columns, 1 < ( A + B ) / ( A B ) < 1 indicates that the first inequality holds for A > B , and the second one for A < B .
j, k V + ( x ) α
Condition
β
Condition
ts α + β α β
in V ( x )
α + β α β
in V + ( x )
E 0 ( ) ϵ T i
0, 0
V I ( α + 1 , β + 1 , x ) > 1 2 1 2 α 2 + 1 4 β 2 + 1 4 0 T 1
V I ( α 1 , β 1 , x ) > 1 2 > 1 2 α 2 + 1 4 β 2 + 1 4 a 2 ( α + β ) > 0 T 2
V I ( α + 1 , β 1 , x ) > 1 2 > 1 2 α 2 + 1 4 β 2 + 1 4 a 2 ( α + 1 ) ( β ) > 0 T 3
V I ( α 1 , β + 1 , x ) > 1 2 > 1 2 α 2 + 1 4 β 2 + 1 4 a 2 ( β + 1 ) ( α ) > 0 T 4
0, 1
V ^ I ( α + 1 , β 1 , x ) > 1 2 > 1 > 1 2 α 2 + 1 4 β 2 + 1 4 1 < α + β α β + 2 < 1 a 2 ( α + 2 ) ( β 1 ) > 0 T 3
V ^ I ( α 1 , β + 1 , x ) > 1 > 1 2 > 1 2 α 2 + 1 4 β 2 + 1 4 1 < α + β α β 2 < 1 a 2 ( α 1 ) ( β + 2 ) > 0 T 4
1, 0
V I ( α + 1 , β 1 , x ) > 0 > 1 2 > 1 > 1 2 α 2 + 1 4 β 2 + 1 4 1 < α + β α β < 1 a 2 α ( β + 1 ) > 0 T 3
V I ( α 1 , β + 1 , x ) > 1 > 1 2 > 0 > 1 2 α 2 + 1 4 β 2 + 1 4 1 < α + β α β < 1 a 2 ( α + 1 ) β > 0 T 4
1, 1
V ^ I ( α + 1 , β + 1 , x ) > 0 > 1 2 > 0 > 1 2 α 2 + 1 4 β 2 + 1 4 1 < α + β α β < 1 1 < α + β + 2 α β < 1 0 T 1
V ^ I ( α 1 , β 1 , x ) > 1 > 1 2 > 1 > 1 2 α 2 + 1 4 β 2 + 1 4 1 < α + β α β < 1 1 < α + β 2 α β < 1 a 2 ( α + β ) > 0 T 2
Table 6. The same as Table 4 for the real Scarf II potential V I I ( α , β , x ) .
Table 6. The same as Table 4 for the real Scarf II potential V I I ( α , β , x ) .
j, k V + ( x ) α R
Condition
ts Re ( t + s ) ϵ T i
0, 0
V I I ( α + 1 , β + 1 , x ) < 1 2 α 2 + 1 4 β 2 + 1 4 α R + 1 2 < 0 a 2 α R + 1 2 2 = E 0 ( ) T 1
V I I ( α 1 , β 1 , x ) < 1 2 α 2 + 1 4 β 2 + 1 4 α R + 1 2 > 1 a 2 α R 1 2 2 < E 0 ( ) T 2
V I I ( α + 1 , β 1 , x ) < 1 2 α 2 + 1 4 β 2 + 1 4 1 2 > 0 a 2 i α I + 1 2 2 complex ( T 2 )
V I I ( α 1 , β + 1 , x ) < 1 2 α 2 + 1 4 β 2 + 1 4 1 2 > 0 a 2 i α I + 1 2 2 complex ( T 2 )
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

Soltész, T.; Pethő, L.F.; Lévai, G. Unified Supersymmetric Description of Shape-Invariant Potentials within and beyond the Natanzon Class. Symmetry 2024, 16, 174. https://doi.org/10.3390/sym16020174

AMA Style

Soltész T, Pethő LF, Lévai G. Unified Supersymmetric Description of Shape-Invariant Potentials within and beyond the Natanzon Class. Symmetry. 2024; 16(2):174. https://doi.org/10.3390/sym16020174

Chicago/Turabian Style

Soltész, Tibor, Levente Ferenc Pethő, and Géza Lévai. 2024. "Unified Supersymmetric Description of Shape-Invariant Potentials within and beyond the Natanzon Class" Symmetry 16, no. 2: 174. https://doi.org/10.3390/sym16020174

APA Style

Soltész, T., Pethő, L. F., & Lévai, G. (2024). Unified Supersymmetric Description of Shape-Invariant Potentials within and beyond the Natanzon Class. Symmetry, 16(2), 174. https://doi.org/10.3390/sym16020174

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