Numerical study of the SWKB condition of novel classes of exactly solvable systems Yuta Nasuda∗ and Nobuyuki Sawado† Department of Physics, Graduate School of Science and Technology, Tokyo University of Science, Noda, Chiba 278-8510, Japan (Dated: January 1, 2021) The supersymmetric WKB (SWKB) condition is supposed to be exact for all known exactly solvable quantum mechanical systems with the shape invariance. Recently, it was claimed that the SWKB condition was not exact for the extended radial oscillator, whose eigenfunctions consisted of the the exceptional orthogonal polynomial, even the system possesses the shape invariance. In this paper, we examine the SWKB condition for the two novel classes of exactly solvable systems: one has the multi-indexed Laguerre and Jacobi polynomials as the main parts of the eigenfunctions, and the other has the Krein–Adler Hermite, Laguerre and Jacobi polynomials. For all of them, one can always remove the ~-dependency from the condition, and it is satisfied with a certain degree of accuracy. I. INTRODUCTION and also a potential which is isospectral to the 1D har- monic oscillator, both of which are exactly solvable but The well-known WKB quantization condition is given are not shape invariant. They arrived at the conjecture by that the shape invariance could be a necessary condition for the exactness of the SWKB condition. xR 1 As was mentioned in Ref. [4], any qualitative argument n V (x) dx = n + π~ (n Z>0) , Z E − 2 ∈ or proof of the exactness of the SWKB condition is still xL p (1) absent at present, and then the conjecture of Ref. [7] in which xL and xR are the “turning points”; V (xL) = has not been proved yet. In order to come closer to the V (xR) = n. In the context of supersymmetric quan- answer of this question, it is worthwhile to examine the tum mechanicsE (SUSY QM) [1–3], the potential V (x) is SWKB condition for new examples which are all exactly formally given by the groundstate eigenfunction of the solvable but are not shape-invariant. Khare and Varshni system φ0(x) as [7] reported the non-exactness of the SWKB condition for the above-mentioned potentials, which are not shape 2 2 2 V (x)= ~ (∂x ln φ (x) ) + ∂ ln φ (x) . (2) invariant. Around the same time, it was also pointed | 0 | x | 0 | out that the SWKB condition is neither exact nor never It should be emphasized that the above potential V (x) worse than WKB for the Abraham–Moses systems [8], corresponds to the vanishing groundstate energy, 0 = 0. which are solvable, by DeLaney et al. [9]. On the other hand, a WKB-like condition inE SUSY Recently, Bougie et al. [10] showed that the SWKB QM (the SWKB condition) proposed by Comtet et al. [4] condition is not always exact for shape-invariant poten- reads tials. Their claim that the shape invariance is not suffi- cient condition for the SWKB exactness is suggestive. b ~ 2 ~ Z However, their analysis is wrong because it lacks the n ( ∂x ln φ0(x) ) dx = nπ (n >0) , ~ Za qE − | | ∈ proper treatment of . We now realize that there still (3) is a possibility that the shape invariance does not guar- 2 where a,b are the roots of (~∂x ln φ0(x) ) = n. This antee the exactness of the SWKB condition. Section II condition is exact for the groundstate| by| construction.E is devoted to a further numerical analysis on the novel For all conventional shape-invariant potentials [5], it has class of shape-invariant potentials. been demonstrated that the SWKB condition (3) repro- In Refs. [11, 12], the authors tried to understand the duces the exact bound-state spectra by Dutt et al. [6]. It non-exactness of the SWKB condition for the systems is well known that the shape invariance is a sufficient (but without shape invariance in a context of the semi-classical arXiv:2004.04927v2 [math-ph] 29 Dec 2020 not a necessary) condition for the exact solvability of the regime, where the SWKB condition is obtained from the Schr¨odinger equation. A natural question that arises here standard WKB formula. Although most of the litera- is whether or not the shape invariance of the potential is ture employs the unit of ~ = 1, to simplify the analy- necessary in order that the SWKB condition reproduces ses, we retain ~ in this paper for rigorous discussions. the exact bound-state spectrum. This was already dis- In fact, as we shall show, ~ is removed completely from cussed by Khare et al. [7] with the Ginocchio potential the SWKB condition (3). It apparently means that the SWKB exactness should not be discussed in the context of the semi-classical regime of the quantum system. Here, a fundamental question arises: what does the ∗ [email protected] SWKB condition imply? The exactness of the condition † [email protected] might indicate that the potential belongs to the class 2 of conventional shape-invariant ones. Breaking of the The authors employed the extended radial oscillator, SWKB condition somehow reflects the discrepancy be- which was equivalent to the exceptional Laguerre or tween the potential of our concern and the corresponding the type II X1-Laguerre polynomial [16–22], to show shape-invariant system, which might be depicted by an that the SWKB condition was not exact for all additive unknown factor. For the moment, we have only a vague shape-invariant potentials. The extended radial oscilla- notion to realize it. Thus, it is worth performing the nu- tor can be seen as a special case of the multi-indexed merical study of several novel classes of solvable systems systems. We note that they alleged that the additive to get an insight for the factor. shape invariance was realized for the parameters ai such In this paper, we study the systems with the multi- that ai+1 = ai + ~ (also in Refs.[23–26]) and expanded indexed Laguerre and Jacobi polynomials [13] as the the superpotential W in power of ~ where the conven- main parts of the eigenfunctions and the Krein–Adler sys- tional shape-invariant superpotentials are the lowest or- tems [14, 15]. They are obtained by deforming the har- der (~0). Note stress that the SWKB condition has no monic oscillator (H), the radial oscillator, or the pseudo- dependencyO on ~ and no room for their expansion param- harmonic oscillator, (L) and the P¨oschl–Teller potential eter ~. The analysis is thus quite questionable. (Another (J) (we set 2m = 1): drawback of their analysis is discussed in Appendix A.) We carry out a further numerical analysis not only (∗) ~2 2 (∗) = ∂x + V (x) , (4) on the extended radial oscillator but also on the multi- H − ω2x2 ~ω = H indexed systems, which are more general cases of the − ∗ ~2g(g 1) novel shape-invariant potentials. ω2x2 + − ~ω(2g + 1) x2 − (∗) V (x)= =L ~2g(g 1) ~2h(h 1) ∗ − + − ~2(g + h)2 A. Multi-indexed Laguerre/Jacobi systems sin2 x cos2 x − =J ∗ The multi-indexed Laguerre and Jacobi polynomi- with the Schr¨odinger equations als [13] are obtained through deformations of two of the three exactly solvable systems (4) via the virtual-state (∗) (∗) (∗) (∗) Z φ (x)= φ (x) (n ≥0), (5) wavefunctions ϕn(x) , which are defined as H n En n ∈ { } 2n~ω = H z g (g− 1 ) (∗) ~ ∗ (L),I 2 2 2 n = 4n ω =L , ϕn (x) :=e z Ln ( z) , (6) E 2 ∗ − 4~ n(n + g + h) =J z 1−g 1 (L),II − ( 2 −g) ∗ 2 2 2 ϕn (x) :=e z Ln (z) , (7) − ωx e 2~ , = H g 1−h (∗) 2 ∗ 2 2 1 1 − ωx g (J),I 1 y 1+ y (g− 2 , 2 −h) φ0 (x)= e 2~ x , =L . ϕ (x) := − Pn (y) , n 2 2 (sin x)g (cos x)h, ∗ =J ∗ (8) 1−g h The shape invariance is achieved by g g + 1 and 2 2 1 1 → (J),II 1 y 1+ y ( 2 −g,h− 2 ) h h+1. We note that the significant parts of the eigen- ϕn (x) := − Pn (y) , → (∗) 2 2 functions φn (x) are described by Hermite polynomials { } (9) Hn for the harmonic oscillator (H), Laguerre polynomials (α) Ln for the radial oscillator (L) and Jacobi polynomials where (α,β) Pn for the P¨oschl–Teller potential (J), respectively. We shall see that for both of the multi-indexed and the ω ξ x, z ξ2 , y cos2x , (10) Krein–Adler systems, the SWKB condition is exact only ≡ r ~ ≡ ≡ for the groundstates n = 0, which is obvious by construc- tion, but they could be a good approximation especially and the parameters g,h must satisfy for sufficiently higher excited states. We shall also discuss the mathematical implications of the results in Sec. IV. 3 1 L: g > max N + , dII + , (11) 2 j 2 1 II. MULTI-INDEXED SYSTEMS AND THE II J: g > max N +2, dj + , SWKB CONDITION 2 1 h> max M +2, dI + . (12) In this section, we examine the SWKB condition for j 2 the novel class of shape-invariant potentials, i.e., the multi-indexed systems. We employ the virtual-state wavefunctions with n = I II I I II II I ∈ D Bougie et al. [10] discuss non-exactness of the SWKB = d1,...,dM d1 ,...,dN and d1 < < DI ∪ DII { II }∪{ } ··· condition for the additive shape-invariant potentials. d , d < < d Z> as the seed solutions, and M 1 ··· N ∈ 0 3 deform φn(x) through the multiple Darboux transfor- For the case of the radial oscillator (L), with the mation{ to construct} the multi-indexed systems.