FRINK Linked Data Fragments server

Matches in SemOpenAlex for { <https://semopenalex.org/work/W2078778880> ?p ?o ?g. }

• W2078778880 endingPage "426" @default.
• W2078778880 startingPage "412" @default.
• W2078778880 abstract "first_page settings Order Article Reprints Font Type: Arial Georgia Verdana Font Size: Aa Aa Aa Line Spacing:    Column Width:    Background: Open AccessArticle Confluent Supersymmetric Partners of Quantum Systems Emerging from the Spheroidal Equation by Axel Schulze-Halberg 1,* and Jie Wang 2 1 Department of Mathematics and Actuarial Science and Department of Physics, Indiana University Northwest, 3400 Broadway, Gary, IN 46408, USA 2 Department of Computer Information Systems, Indiana University Northwest, 3400 Broadway, Gary, IN 46408, USA * Author to whom correspondence should be addressed. Symmetry 2015, 7(2), 412-426; https://doi.org/10.3390/sym7020412 Received: 8 March 2015 / Revised: 14 April 2015 / Accepted: 20 April 2015 / Published: 22 April 2015 (This article belongs to the Special Issue Supersymmetry 2014) Download Download PDF Download PDF with Cover Download XML Download Epub Browse Figures Versions Notes Abstract: We construct confluent supersymmetric partners of quantum systems that emerge from the spheroidal equation. Properties of the systems and of their transformed counterparts are discussed. Keywords: confluent supersymmetry algorithm; spheroidal equation; series solutions 1. IntroductionThe supersymmetry (SUSY) formalism is a method for generating solvable models in quantum mechanics. Mathematically based on the Darboux transformation [1–3], it converts a known quantum system into a transformed counterpart, the so-called SUSY partner. As this formalism is compatible with basically any exactly-solvable quantum model governed by the Schrödinger equation and many systems described by the Dirac equation, a large amount of applications can be found in the literature. An overview and details on the general topic can be found in one of the reviews [4–6] and references therein. The formalism of SUSY can be split into two different versions, commonly referred to as standard and confluent algorithm, respectively. While the vast majority of applications uses the standard algorithm, its confluent counterpart has received considerably less attention. As a consequence, many of its properties are not yet well understood. Recent works on the confluent algorithm include its application to the inverted oscillator [7] and to the Dirac equation for pseudoscalar potentials [8]. Further results concern its equivalence with the zero-mode supersymmetry scheme [9], a representation using a differential formula [10,11] and higher-order confluent SUSY transformations [12]. The purpose of the present note is to apply the confluent SUSY algorithm to models that are governed by the spheroidal equation [13], a particular case of the confluent Heun class [14]. The spheroidal equation has many applications in physics, including a sextic oscillator system with centrifugal barrier [15], a quantum model featuring two Coulombian centers [16], the dynamical theory of zonal tides in a global ocean [17], just to mention a few. While solutions of the spheroidal equation cannot be given in closed form, for certain parameter settings one can construct them through a series expansion in terms of Gegenbauer polynomials [18]. In this work we will construct such series solutions for two systems that are associated with the spheroidal equation. The first system is a generalization of the Mathews-Lakshmanan oscillator [19–21], while the second features a potential of Pöschl-Teller type that contains an aditional cosine term. After obtaining the solutions, we will generate confluent SUSY partners of the two systems. The remainder of this paper is organized as follows. In Section 2 we give a summary of the confluent SUSY algorithm. Section 3 is devoted to the spheroidal equation, the algorithm for the construction of its solutions and a discussion of their properties. In Section 4 the findings are applied to the two particular quantum systems mentioned above, resulting in the generation of their confluent SUSY partner systems.Implementation notes: All results on spectral values and solutions that are presented subsequently were obtained by implementing the algorithm described in Section 2.2 using MATLAB (The MathWorks, Inc., Natick, MA, USA). Calculations were done at 40-digit precision. All digits of numerical values stated throughout this work are significant. 2. The Confluent SUSY AlgorithmWe will now give a brief review of the confluent SUSY formalism. For further details, particularly on spectral design, we refer the reader to [11,22,23] and references therein. We start out by considering an initial equation that has the general form f Ψ ″ + g Ψ ′ + ( E − V 0 ) Ψ = 0 , (1) where the coefficient functions f, g are required to be sufficiently smooth and the potential V0 must be continuous. Next, for a positive integer N ≥ 2, assume that N solutions u0, vj, j = 0,…, N − 2, of the equation f u ″ + g u ′ + ( λ − V 0 ) u = 0 , (2) are given, introducing a constant λ, commonly referred to as factorization energy. Note that the solutions u0, vj, j = 0,…, N − 2 do not have to be normalizable in the L2-sense, that is, they are not required to represent bound states. We will return to this point at the end of this paragraph. Now we define functions uj, j = 1,…, N − 1, as follows u j = ∑ k = 0 j − 1 ∂ k v k ∂ λ k + 1 j ! ∂ j u 0 ∂ λ j , j = 1 , … , N − 1. (3) As an example, let us mention that in the simplest case N = 2 we obtain from Equation (3) a single function u1, given by u 1 = v 0 + ∂ u 0 ∂ λ , (4) where v0 and u0 are solutions of Equation (2). We will now use the functions uj, j = 0,…, N − 1, called transformation functions, to set up a SUSY transformation of order N, applied to the solution Ψ of Equation (1). It reads Φ = f N 2 W u 0 , … , u N − 1 , Ψ W u 0 , … , u N − 1 . (5) Here, the symbol W stands for the Wronskian of the functions in its index. The Equation (5) is a solution of the supersymmetric partner equation to Equation (1), given by f Φ ″ + g Φ ′ + ( E − V 1 ) Φ = 0. (6) The new potential V1 is linked to its initial counterpart V0 in Equation (1) by means of the following relation [24] V 1 = V 0 + N ( N − 2 ) [ ( f ′ 2 f ) 2 − f ″ 2 ] + N g f ′ 2 f − N g ′ − f ′ W u 0 , u 1 , … , u N − 1 ′ W u 0 , u 1 , … , u N − 1 − 2 f d d x W u 0 , u 1 , … , u N − 1 ′ W u 0 , u 1 , … , u N − 1 . (7) The functions V0 and V1 are usually called supersymmetric partner potentials. There are two principal differences between the confluent SUSY algorithm and its standard counterpart, the first of which contains the transformation order. The confluent algorithm uses a SUSY transformation of order two or higher, while in the conventional case first-order transformations are possible. In other words, there is no first-order confluent SUSY transformation. In applications, such as the recent works in optics [25,26], one must consider SUSY transformations of at least order two, before the confluent SUSY algorithm becomes applicable. The second principal difference between the confluent SUSY algorithm and its standard counterpart lies in the definition of the transformation functions uj, j = 0,…, N − 1, affecting both transformed solution (5) and potential (7). In the standard framework of SUSY each transformation function is a solution of the initial Equation (1), associated to its own factorization energy, while the confluent algorithm requires a single factorization energy and defines the transformation functions through the differential formula Equation (3). Note that the latter functions can also be expressed through an integral representation. Since this representation is equivalent to Equation (3) and for our purposes less convenient, we omit to state it here. Before we conclude this paragraph, let us comment once more on the transformation functions uj, j = 0,…, N − 1. The main purpose of these functions is to render the transformed potential (7) real-valued and free of singularities inside its domain. This condition can be fulfilled if the Wronskian that appears in the denominator of the terms in Equation (7), is free of zeros, which essentially relates to linear independence of the transformation functions. However, in the present models that will be considered in the subsequent sections, we only have access to solutions representing bound states. As a consequence, we do not need to consider transformation functions that are of more general type. 3. The Spheroidal EquationThe spheroidal equation, equipped with boundary conditions of Dirichlet type, can be written in the following way ( 1 − y 2 ) χ ″ − 2 y χ ′ + [ A − h 2 y 2 + M 2 y 2 − 1 ] χ = 0 , y ∈ ( − 1 , 1 ) (8) χ ( − 1 ) = χ ( 1 ) = 0 , (9) where A, h and M are constants. Comparison with Equation (1) shows that Equation (8) is a special case. While Equation (8) cannot be solved in closed form, under certain conditions it admits series solutions in terms of Gegenbauer polynomials. For the sake of convenience, we will now briefly summarize how to construct such series solutions. For further details, the reader may refer to [18]. It is known that there exists an infinite, discrete and ascending sequence of characteristic values A = Aj, j ∈ ℕ ∪ {0}, for each of which the problem Equation (8), Equation (9) admits a solution that has exactly j zeros inside (−1, 1). The solutions form an orthogonal set with respect to the usual inner product of the Hilbert space L2(−1, 1). All solutions have parity, but do not feature a closed-form representation. Instead, they can be written as a series expansion in terms of Gegenbauer polynomials T [27] via one of the formulas χ 2 j = ( 1 − y 2 ) M 2 ∑ n = 0 ∞ ( d 2 n ) | A = A 2 j T 2 n M ( y ) (10) χ 2 j + 1 = ( 1 − y 2 ) M 2 ∑ n = 0 ∞ ( d 2 n + 1 ) | A = A 2 j + 1 T 2 n + 1 M ( y ) , (11) for j ∈ ℕ ∪ {0}, where Equation (10) and Equation (11) represent even and odd parity functions, respectively. The expansion coefficients in the series are determined through a three-term recursion relation that reads a n d n − 2 + b n d n + c n d n + 2 = 0 , n ∈ ℕ ∪ { 0 } , d − 2 = d − 1 = 0 , (12) where the following abbreviations are involved: a n = n ( n − 1 ) h 2 ( 2 n + 2 M − 1 ) ( 2 n + 2 M − 3 ) (13) b n = h 2 { 2 [ ( n + M ) ( n + M + 1 ) − M 2 ] − 1 } ( 2 n + 2 M + 3 ) ( 2 n + 2 M − 1 ) ) + ( n + M ) ( n + M + 1 ) − A (14) c n = h 2 ( n + 2 M + 1 ) ( n + 2 M + 2 ) ( 2 n + 2 M + 3 ) ( 2 n + 2 M + 5 ) . (15) It is important to point out that the series Equations (10) and (11) are in general divergent. Convergence is guaranteed only if the parameter A in Equation (14) takes a characteristic value A = Aj, j ∈ ℕ ∪ {0}. The algorithm for determining these values, as well as the expansion coefficients in the series Equations (10) and (11), consists of two steps.Step 1. We solve the recursion (12) with respect to d2/d0, first substituting n = 0 and the second time n = 2. This gives the following two relations d 2 d 0 = − b 0 c 0 d 2 d 0 = − a 2 b 2 + c 2 d 4 d 2 . (16) Combining these two relations, we obtain − b 0 c 0 = − a 2 b 2 + c 2 d 4 d 2 . (17) Now we will modify the right side of this equation. We shift all indices n in Equation (12) up by two, set n = 4, solve for d4/d2, and substitute the result into the right-hand side of Equation (16). An iteration of this process gives rise to the following continued fraction expansion: − b 0 c 0 = − a 2 b 2 + c 2 − a 4 b 4 + c 4 − a 6 b 6 + c 6 … . (18) We stop this iteration after N > 0 steps by reasoning that dN+2/dN is very small, so we set dN+2/dN = 0. This turns Equation (18) into an algebraic equation that can be solved for the parameter A, yielding approximations for the characteristic values belonging to solutions of even parity. Observe that Equation (18) generally has several solutions, each of which is associated with a characteristic value. Also, the number of solutions increases with the number N of steps that are done. In a similar way we obtain a continued fraction expansion for odd indices: − b 1 c 1 = − a 3 b 3 + c 3 − a 5 b 5 + c 5 − a 7 b 7 + c 7 … , (19) that leads to the characteristic values associated with the odd solutions.Step 2. We set d0 = 1. Relation (16) then gives d2, which can be plugged into the recursion relation (12) for n = 2 in order to obtain d4. The continuation of this process yields the following expansion coefficients d6, d8, d10,…. The coefficients for odd indices are constructed in a similar manner, starting out with d1 = 1. The recursion (12) for n = 1 generates d3, which can be substituted into (12) for n = 3 in order to give d5. This iteration produces the subsequent expansion coefficients d7, d9, d11,…. The solutions of the boundary-value problem (8), (9) can now be approximated by plugging the expansion coefficients and the characteristic values into the series representations (10), (11). 4. ApplicationsWe will now use point transformations for converting the spheroidal equation, such that it describes different physically interesting models. In particular, we focus on two systems, featuring a double oscillator potential of Mathews-Lakshmanan type and a generalized trigonometric Pöschl-Teller interaction, respectively. 4.1. A Double Oscillator SystemThe first model that will be considered here can be generated by means of the following point transformation. χ [ y ( x ) ] = ( μ x 2 − 1 ) − 1 4 Ψ ( x ) , y ( x ) = − μ x , (20) where µ > 0 is a real constant. In addition, we make the following definitions in Equation (8) A = 2 E μ + 1 μ 2 − 1 μ − 1 4 h = − q μ M = 1 2 μ μ + 4 μ − 4 , (21) for a real constant q. Note that although h can become imaginary, it enters in the spheroidal equation and in the recursion Equations (12)–(15) for the series coefficients merely as the real-valued h2, therefore preventing the spectral values from becoming complex. After substituting both Equation (20) and the settings (21) into the boundary-value problem (8), (9), we obtain after some elementary manipulations ( 1 − μ x 2 ) Ψ ″ − μ x Ψ ′ + [ 2 E − q x 2 − ( 1 − μ ) x 2 1 − μ x 2 ] Ψ = 0 , x ∈ ( − 1 μ , 1 μ ) (22) Ψ ( − 1 μ ) = Ψ ( 1 μ ) = 0. (23) Note that we applied the point transformation (20) to the initial boundary conditions (9) as well, yielding the result Equation (23). The model described by Equations (22) and (23) can be interpreted as the quantum counterpart of a classical system that is governed by the following Lagrangian ℒ = x ˙ 2 2 − 2 μ x 2 − q 2 x 2 − ( 1 − μ ) x 2 2 − 2 μ x 2 , (24) The specific form of the kinetic term in Equation (24) can be explained by introducing a position-dependent mass function. This function m and the potential V featured in Equation (24) have the explicit form m = 1 1 − μ x 2 V 0 = q 2 x 2 + ( 1 − μ ) x 2 2 − 2 μ x 2 . (25) The interaction V represents a nonlinear, two-parameter double oscillator, containing a harmonic term and a contribution of Mathews-Lakshmanan type. By virtue of its two parameters µ and q, the potential (25) has several well-known limiting cases, a few of which are listed in Table 1.Next, we will show that the point transformation (20) preserves orthogonality and normalizability of the bound-state solutions stemming from the initial system (8), (9), as long as we require the solutions of Equations (22) and (23) to be located in the weighted Hilbert space L w 2 ( − 1 / μ , 1 / μ ) weight function w = 1 1 − μ x 2 . (26) Let k, l be two nonnegative integers ands let N be a constant. We obtain by combining orthogonality of the Equations (10) and (11) with the point transformation (20) and (26) the following result N δ k l = ∫ − 1 1 χ k * χ l d y = ∫ − 1 μ 1 μ Ψ k * Ψ l μ 1 1 − μ x 2 d x = μ ∫ − 1 μ 1 μ Ψ k * Ψ l 1 1 − μ x 2 d x . (27) Combining the very left and the very right term in this chain of equations, we obtain ∫ − 1 μ 1 μ Ψ k * Ψ l 1 1 − μ x 2 d x = N μ δ k l , (28) implying orthogonality of the set (Ψj) with respect to the weighted Hilbert space L w 2 ( − 1 / μ , 1 / μ ). Normalizability in the same Hilbert space follows then from ‖ Ψ k ‖ = ∫ − 1 μ 1 μ Ψ k * Ψ k 1 1 − μ x 2 d x = N μ ∫ − 1 1 χ k * χ k d y = N μ ‖ χ k ‖ . (29) Now that we have established orthogonality and normalizability of the solutions, we are able to construct an explicit form for them. First, the point transformation (20) gives Ψ ( x ) = ( μ x 2 − 1 ) 1 4 χ ( − μ x ) (30) In the next step we substitute the series form (10), (11) for the solutions of the spheroidal Equation (8) into Equation (30). Note that due to the settings (21), the characteristic values of A will translate into discrete spectral values of E, denoted by Ej, j a nonnegative integer. After adding an index for taking into account the bound-state character of the solutions, we obtain Ψ 2 j = ( 1 − μ x 2 ) 1 4 + 1 4 μ μ + 4 μ − 4 ∑ n = 0 ∞ ( d 2 n ) | E = E 2 j T 2 n 1 2 μ μ + 4 μ − 4 ( − μ x ) (31) Ψ 2 j + 1 = ( 1 − μ x 2 ) 1 4 + 1 4 μ μ + 4 μ − 4 ∑ n = 0 ∞ ( d 2 n + 1 ) | E = E 2 j + 1 T 2 n + 1 1 2 μ μ + 4 μ − 4 ( − μ x ) , (32) We are now ready to compute confluent SUSY partners of our model (22), (23), the governing equation of which is clearly a special case of Equation (8) for the settings f = 1 − μ x 2 g = − μ x . (33) For the sake of simplicity, here we will restrict ourselves to transformations of second order, given by the relation (5) for N = 2. The principal reason for this restriction lies in the computational effort that is required for the calculation of the transformed potential (7), in particular for the Wronskians and their derivatives. Each function uj, j = 0,…, N − 1, is constructed by means of the algorithm shown in Section 4, such that it is represented by a series of the form (31) or (32). Consequently, Wronskians for values of N ≥ 3 contain at least three of such series including at least their second derivatives. This means a large amount of computational work, especially since we must maintain our prescribed numerical precision. For these reasons, here we prefer to use only the case N = 2 (second-order transformations), but we will get back to this issue below. Now we need two transformation functions u0 and u1, the first of which must be a solution of Equation (22). After determining u0 we can generate u1 by means of Equation (4). While in general we can choose a transformation function u0 for any factorization energy λ, in the present case only bound-state solutions of the series form (31), (32) are available. Next we focus on the function u1 as given in Equation (4). Since v0 must be a solution of Equation (8) for the same energy λ as u0, we can only pick u0 itself or a multiple of it. Note that applying the reduction-of-order method for generating another solution of Equation (8) would lead to integrals of rational functions that feature high-degree polynomials in the denominator. Such integrals can in general not be resolved in closed form. Now we observe that the transformation functions u0 and u1 enter in the potential (7) only through their Wronskian. Consequently, if v0 in Equation (4) and u0 were not linearly independent, we would obtain the same transformed potential as for v0 = 0. For a nonnegative integer j we make the choice u 1 = ∂ u 0 ∂ λ | λ = E j , (34) while u0 will be one of the solutions given in Equations (31) and (32), pertaining to the spectral value Ej. We can now perform a confluent SUSY transformation of second order by substituting the settings Equation (34) into Equations (5) and (7). The representation (31) for u0 and u1 requires truncation of the series. Despite its convergence being very fast [18], we still have to take into account a large number of terms. For this reason we do not show the explicit calculations that lead to the subsequently shown results.In the Figures 1–4 we display examples of potentials V1 that result from a confluent SUSY transformation of second order. These potentials are obtained from Equation (7) by substituting N = 2 and the Wronskians of the transformation functions u0, u1. Each of the four plots contained in the latter figures shows the initial potential and two of its confluent SUSY partners at the factorization energies λ = E0 and λ = E1, respectively. We chose a variety of parameter settings for µ and q in order to start out from four different shapes of the initial potential. 4.2. A generalized trigonometric Pöschl-Teller systemWe will now generate another quantum model from the spheroidal Equation (8). Since the process of construction follows similar steps as in the previous section, we will not give details for all our calculations. This time we employ the following point transformation χ [ y ( x ) ] = 1 cos ( μ 2 x ) Ψ ( x ) , y ( x ) = − sin ( μ 2 x ) , (35) while we keep the parameter settings (21). The transformation (35) converts the boundary-value problem (8), (9) into the new form Ψ ″ + ( E − V 0 ) Ψ = 0 , x ∈ ( − π 2 π , π 2 π ) (36) Ψ ( − π 2 π ) = Ψ ( π 2 π ) = 0 (37) where the potential V0 is given explicitly by V 0 = 1 − μ 2 μ cos 2 ( μ 2 x ) − q 4 μ cos ( 2 μ x ) + q 4 μ − 1 2 μ + 1 2 (38) We observe that in contrast to the prior system (22), (23), the present governing Equation (36) has standard Schrödinger form. By combining the point transformation (35) with the series (10) and (11), we obtain solutions Ψj, j a nonnegative integer, of bound-state type to our boundary-value problem (36), (37) in the form Ψ 2 j = cos 1 2 + 1 2 μ μ + 4 μ − 4 ( μ 2 z ) ∑ n = 0 ∞ ( d 2 n ) | E = E 2 j T 2 n 1 2 μ μ + 4 μ − 4 [ − sin ( μ 2 z ) ] (39) Ψ 2 j + 1 = cos 1 2 + 1 2 μ μ + 4 μ − 4 ( μ 2 z ) ∑ n = 0 ∞ ( d 2 n + 1 ) | E = E 2 j + 1 T 2 n + 1 1 2 μ μ + 4 μ − 4 [ − sin ( μ 2 z ) ] (40) Orthogonality and normalizabiliy of these solutions with respect to the Hilbert space L 2 ( − π / 2 μ , π / 2 μ ) can be established similarly to Equations (27)–(29), such that we omit to give details. In order to perform our confluent SUSY transformations, we need to determine factorization functions. This time we will compute second and third-order transformations, such that we need three transformation functions u0, u1 and u2. Note that the computational effort for determining the Wronskians and their derivatives in Equation (7) is significantly lower then in the previous application, because our governing Equation (36) takes a much simpler form when compared to its counterpart (22). Now, for a nonnegative integer j we pick u 1 = ∂ u 0 ∂ λ | λ = E j u 2 = 1 2 ∂ 2 u 0 ∂ λ 2 | λ = E j , (41) while u0 will be taken from the solution set given by Equations (39) and (40), pertaining to the spectral value Ej. Figures 5–8 show potentials V1 that result from confluent SUSY transformations of second and third order, respectively. In order to start out from different shapes of the initial potential, we picked a variety of parameter settings for q and µ. Due to the length of the expressions involved in the calculations of the SUSY partner potentials, we omit to give details. 5. ConclusionsIn this work we have constructed confluent SUSY partners of two quantum systems that are related to the spheroidal equation by means of a point transformation. Although our models do not admit closed-form solutions, we can generate SUSY partners through a series representation for them. Since these representations are only available for bound-state type solutions, we are restricted in the choice of transformation functions and factorization energies. Our approach allows for an extension to more general models by varying the point transformation that connects the spheroidal equation and the targeted system. AcknowledgmentsJie Wang acknowledges Indiana University Northwest Research Support Grant.PACS classifications: 03.65.Ge, 03.65.PmAuthor ContributionsJie Wang implemented and tested the algorithm, which was used by Axel Schulze-Halberg to perform calculations. The final version of the manuscript was written in collaboration.Conflicts of InterestThe authors declare no conflict of interest.ReferencesBagrov, V.G.; Samsonov, B.F. Darboux transformation of the Schrödinger equation. Phys. Part. Nucl. 1997, 28, 374–397. [Google Scholar]Darboux, G. Sur une proposition relative aux équations lińeaires. Comptes Rendus Acad. Sci. Paris 1882, 94, 1456–1459. [Google Scholar]Matveev, V.B.; Salle, M.A. Darboux Transformations and Solitons; Springer: Berlin, Germany, 1991. [Google Scholar]Cooper, F.; Khare, A.; Sukhatme, U. Supersymmetry and Quantum Mechanics. Phys. Rep. 1995, 251, 267–388. [Google Scholar]Fernandez, C.; D.J. Supersymmetric quantum mechanics. AIP Conf. Proc. 2010, 1287, 3–36. [Google Scholar]Junker, G. Supersymmetric Methods in Quantum and Statistical Physics. Texts and Monographs in Physics; Springer: Heidelberg, Germany, 1996. [Google Scholar]Bermudez, D.; Fernandez, C.; D.J. Factorization method and new potentials from the inverted oscillator. Ann. Phys. 2013, 333, 290–306. [Google Scholar]Contreras-Astorga, A.; Schulze-Halberg, A. The confluent supersymmetry algorithm for Dirac equations with pseudoscalar potentials. J. Math. Phys. 2014, 55. [Google Scholar] [CrossRef]Contreras-Astorga, A.; Schulze-Halberg, A. The generalized zero-mode supersymmetry scheme and the confluent algorithm. Ann. Phys. 2015, 354, 353–364. [Google Scholar]Bermudez, D.; Fernandez, C.; D.J.; Fernandez-Garcia, N. Wronskian differential formula for confluent supersymmetric quantum mechanics. Phys. Lett. A 2012, 376, 692–696. [Google Scholar]Fernandez, C.; D.J.; Salinas-Hernandez, E. Hyperconfluent third-order supersymmetric quantum mechanics. J. Phys. A 2011, 44. [Google Scholar] [CrossRef]Schulze-Halberg, A. Wronskian representation for confluent supersymmetric transformation chains of arbitrary order. Eur. Phys. J. Plus 2013, 128, 68–85. [Google Scholar]Slavianov, S.Y.; Lay, W. Special Functions: A Uunified Theory Based on Singularities; Oxford University Press: New York, NY, USA, 2000. [Google Scholar]Ronveaux, A. Heun’s Differential Equations; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]Raduta, A.A.; Buganu, P. Application of the sextic oscillator with a centrifugal barrier and the spheroidal equation for some X(5) candidate nuclei. J. Phys. G 2013, 40. [Google Scholar] [CrossRef]Figueiredo, B.D.B. On some solutions to generalized spheroidal wave equations and applications. J. Phys. A 2002, 35. [Google Scholar] [CrossRef]O’Connor, W.P. On the Application of the Spheroidal Wave Equation to the Dynamical Theory of the Long-Period Zonal Tides in a Global Ocean. Proc. Royal Soc. Lond. A 1992, 439, 189–196. [Google Scholar]Morse, P.M.; Feshbach, H. Methods of Theoretical Physics; Mc Graw-Hill: New York, NY, USA, 1953; Volume 2. [Google Scholar]Mathews, P.M.; Lakshmanan, M. On a unique nonlinear oscillator. Quart. Appl. Math. 1974, 32, 215–218. [Google Scholar]Midya, B.; Roy, B. A generalized quantum nonlinear oscillator. J. Phys. A 2009, 42. [Google Scholar] [CrossRef]Schulze-Halberg, A.; Morris, J.R. Special function solutions of a spectral problem for a nonlinear quantum oscillator. J. Phys. A 2012, 45. [Google Scholar] [CrossRef]Bermudez, D. Algebras de Heisenberg Polinomiales y Ecuaciones de Painleve. Ph.D. Thesis, Physics Department, Cinvestav, Mexico, 2013. [Google Scholar]Fernandez, C.; D.J.; Salinas-Hernandez, E. Wronskian formula for confluent second-order supersymmetric quantum mechanics. Phys. Lett. A 2005, 338, 13–18. [Google Scholar]Schulze-Halberg, A.; Pozdeeva, E.; Suzko, A. Explicit Darboux transformations of arbitrary order for generalized time-dependent Schrödinger equations. J. Phys. A 2009, 42. [Google Scholar] [CrossRef]Miri, M.-A.; Heinrich, M.; El-Ganainy, R.; Christodoulides, D.N. Supersymmetric optical structures. Phys. Rev. Lett. 2013, 110. [Google Scholar] [CrossRef]Miri, M.-A.; Heinrich, M.; Christodoulides, D.N. SUSY-inspired one-dimensional transformation optics. Optica 2014, 1, 89–95. [Google Scholar]Abramowitz, M.; Stegun, I. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Dover Publications: New York, NY, USA, 1964. [Google Scholar] Figure 1. Graph of the initial potential (light-gray curve) and two of its second-order, confluent supersymmetry (SUSY) partners. We set q = 2, µ = 0.2, λ = E0 ≈ 2.072601182582 (gray curve) and λ = E1 ≈ 3.5831316 (black curve). Figure 1. Graph of the initial potential (light-gray curve) and two of its second-order, confluent supersymmetry (SUSY) partners. We set q = 2, µ = 0.2, λ = E0 ≈ 2.072601182582 (gray curve) and λ = E1 ≈ 3.5831316 (black curve). Figure 2. Graph of the initial potential (light-gray curve) and two of its second-order, confluent supersymmetry (SUSY) partners. We set q = 2, µ = 1.1, λ = E0 ≈ 0.5308008 (gray curve) and λ = E1 ≈ 2.5414273 (black curve). Figure 2. Graph of the initial potential (light-gray curve) and two of its second-order, confluent supersymmetry (SUSY) partners. We set q = 2, µ = 1.1, λ = E0 ≈ 0.5308008 (gray curve) and λ = E1 ≈ 2.5414273 (black curve). Figure 3. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = −2, µ = 0.6, λ = E0 ≈ −0.0926509 (gray curve) and λ = E1 ≈ 1.058759001 (black curve). Figure 3. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = −2, µ = 0.6, λ = E0 ≈ −0.0926509 (gray curve) and λ = E1 ≈ 1.058759001 (black curve). Figure 4. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = −2, µ = 1.1, λ = E0 ≈ 0.4957628 (gray curve) and λ = E1 ≈ 1.5751028 (black curve). Figure 4. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = −2, µ = 1.1, λ = E0 ≈ 0.4957628 (gray curve) and λ = E1 ≈ 1.5751028 (black curve). Figure 5. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = 2, µ = 1.05, λ = E0 ≈ 0.23187263 (gray curve) and λ = E1 ≈ 2.53550068 (black curve). Figure 5. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = 2, µ = 1.05, λ = E0 ≈ 0.23187263 (gray curve) and λ = E1 ≈ 2.53550068 (black curve). Figure 6. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = −2, µ = 1.05, λ = E0 ≈ 0.7978696 (gray curve) and λ = E1 ≈ 1.53303989 (black curve). Figure 6. Graph of the initial potential (light-gray curve) and two of its second-order, confluent SUSY partners. We set q = −2, µ = 1.05, λ = E0 ≈ 0.7978696 (gray curve) and λ = E1 ≈ 1.53303989 (black curve). Figure 7. Graph of the initial potential (light-gray curve) and two of its third-order, confluent SUSY partners. We set q = 2, µ = 0.9, λ = E0 ≈ 1.225181834 (gray curve) and λ = E1 ≈ 2.5291772 (black curve). Figure 7. Graph of the initial potential (light-gray curve) and two of its third-order, confluent SUSY partners. We set q = 2, µ = 0.9, λ = E0 ≈ 1.225181834 (gray curve) and λ = E1 ≈ 2.5291772 (black curve). Figure 8. Graph of the initial potential (light-gray curve) and two of its third-order, confluent SUSY partners. We set q = −2, µ = 0.9, λ = E0 ≈ −0.1871158 (gray curve) and λ = E1 ≈ 1.3976904 (black curve). Figure 8. Graph of the initial potential (light-gray curve) and two of its third-order, confluent SUSY partners. We set q = −2, µ = 0.9, λ = E0 ≈ −0.1871158 (gray curve) and λ = E1 ≈ 1.3976904 (black curve). Table 1. Limiting cases of the potential V in Equation (25). Table 1. Limiting cases of the potential V in Equation (25). Parameter SettingPotentialName of Special Caseq = −1, µ → 0V0 = 0Free particleq = 1, µ → 0V0 = x2Harmonic oscillatorq = −3, µ → 0V0 = − x2Inverted harmonic oscillatorq = 0 V 0 = ( 1 − μ ) x 2 2 − 2 μ x 2Mathews-Lakshmanan oscillator © 2015 by the authors; licensee MDPI, Basel, Switzerland This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/). Share and Cite MDPI and ACS Style Schulze-Halberg, A.; Wang, J. Confluent Supersymmetric Partners of Quantum Systems Emerging from the Spheroidal Equation. Symmetry 2015, 7, 412-426. https://doi.org/10.3390/sym7020412 AMA Style Schulze-Halberg A, Wang J. Confluent Supersymmetric Partners of Quantum Systems Emerging from the Spheroidal Equation. Symmetry. 2015; 7(2):412-426. https://doi.org/10.3390/sym7020412 Chicago/Turabian Style Schulze-Halberg, Axel, and Jie Wang. 2015. Confluent Supersymmetric Partners of Quantum Systems Emerging from the Spheroidal Equation Symmetry 7, no. 2: 412-426. https://doi.org/10.3390/sym7020412 Find Other Styles Article Metrics No No Article Access Statistics For more information on the journal statistics, click here. Multiple requests from the same IP address are counted as one view." @default.
• W2078778880 created "2016-06-24" @default.
• W2078778880 creator A5017867619 @default.
• W2078778880 creator A5075877965 @default.
• W2078778880 date "2015-04-22" @default.
• W2078778880 modified "2023-09-28" @default.
• W2078778880 title "Confluent Supersymmetric Partners of Quantum Systems Emerging from the Spheroidal Equation" @default.
• W2078778880 cites W1971135854 @default.
• W2078778880 cites W1978315564 @default.
• W2078778880 cites W2019806305 @default.
• W2078778880 cites W2020350214 @default.
• W2078778880 cites W2029318039 @default.
• W2078778880 cites W2031673422 @default.
• W2078778880 cites W2042822903 @default.
• W2078778880 cites W2062033781 @default.
• W2078778880 cites W2067482542 @default.
• W2078778880 cites W2079163912 @default.
• W2078778880 cites W2102245148 @default.
• W2078778880 cites W2171340515 @default.
• W2078778880 cites W2562397760 @default.
• W2078778880 cites W3098311644 @default.
• W2078778880 cites W3099062828 @default.
• W2078778880 cites W3100172235 @default.
• W2078778880 cites W3101245455 @default.
• W2078778880 cites W3105536310 @default.
• W2078778880 doi "https://doi.org/10.3390/sym7020412" @default.
• W2078778880 hasPublicationYear "2015" @default.
• W2078778880 type Work @default.
• W2078778880 sameAs 2078778880 @default.
• W2078778880 citedByCount "7" @default.
• W2078778880 countsByYear W20787788802015 @default.
• W2078778880 countsByYear W20787788802016 @default.
• W2078778880 countsByYear W20787788802017 @default.
• W2078778880 countsByYear W20787788802018 @default.
• W2078778880 crossrefType "journal-article" @default.
• W2078778880 hasAuthorship W2078778880A5017867619 @default.
• W2078778880 hasAuthorship W2078778880A5075877965 @default.
• W2078778880 hasBestOaLocation W20787788801 @default.
• W2078778880 hasConcept C121332964 @default.
• W2078778880 hasConcept C33332235 @default.
• W2078778880 hasConcept C37914503 @default.
• W2078778880 hasConcept C41008148 @default.
• W2078778880 hasConcept C62520636 @default.
• W2078778880 hasConcept C74650414 @default.
• W2078778880 hasConcept C84114770 @default.
• W2078778880 hasConceptScore W2078778880C121332964 @default.
• W2078778880 hasConceptScore W2078778880C33332235 @default.
• W2078778880 hasConceptScore W2078778880C37914503 @default.
• W2078778880 hasConceptScore W2078778880C41008148 @default.
• W2078778880 hasConceptScore W2078778880C62520636 @default.
• W2078778880 hasConceptScore W2078778880C74650414 @default.
• W2078778880 hasConceptScore W2078778880C84114770 @default.
• W2078778880 hasIssue "2" @default.
• W2078778880 hasLocation W20787788801 @default.
• W2078778880 hasLocation W20787788802 @default.
• W2078778880 hasOpenAccess W2078778880 @default.
• W2078778880 hasPrimaryLocation W20787788801 @default.
• W2078778880 hasRelatedWork W1972665047 @default.
• W2078778880 hasRelatedWork W1975516339 @default.
• W2078778880 hasRelatedWork W1983117717 @default.
• W2078778880 hasRelatedWork W2066832011 @default.
• W2078778880 hasRelatedWork W2086733001 @default.
• W2078778880 hasRelatedWork W2146691474 @default.
• W2078778880 hasRelatedWork W2782882015 @default.
• W2078778880 hasRelatedWork W3102427995 @default.
• W2078778880 hasRelatedWork W3103234735 @default.
• W2078778880 hasRelatedWork W3130965325 @default.
• W2078778880 hasVolume "7" @default.
• W2078778880 isParatext "false" @default.
• W2078778880 isRetracted "false" @default.
• W2078778880 magId "2078778880" @default.
• W2078778880 workType "article" @default.