Next Article in Journal
Lie Symmetry Analysis, Particular Solutions and Conservation Laws of Benjiamin Ono Equation
Previous Article in Journal
Periodic Orbits of Nonlinear Ordinary Differential Equations Computed by a Boundary Shape Function Method
Previous Article in Special Issue
Dyadic Helmholtz Green’s Function for Electromagnetic Wave Transmission/Diffraction through a Subwavelength Nano-Hole in a 2D Quantum Plasmonic Layer: An Exact Solution Using “Contact Potential”-like Dirac Delta Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Supersymmetric Partners of the One-Dimensional Infinite Square Well Hamiltonian: Special Cases

by
Manuel Gadella
*,† and
Carlos San Millán
Departamento de Física Teórica, Atómica y Optica, Universidad de Valladolid, Paseo Belén 7, 47011 Valladolid, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2022, 14(7), 1314; https://doi.org/10.3390/sym14071314
Submission received: 20 May 2022 / Revised: 16 June 2022 / Accepted: 21 June 2022 / Published: 24 June 2022
(This article belongs to the Special Issue Contact Interactions and Solvable Models in Quantum Mechanics)

Abstract

:
In a previous paper, we used a classification of the self adjoint extensions, also called self-adjoint determinations, of the differential operator d 2 / d x 2 in order to obtain the whole list of Supersymmetric (SUSY) partners of those selfadjoint determinations for which the ground state has strictly positive energy. The existence of self adjoint determinations with a ground state of zero or even negative energy is a proved fact. In this paper, we analyze the possibility of constructing SUSY partners for those determinations. We also study those cases for which the ground state has a degeneracy, the study of their SUSY partners should be analyzed separately. So far, we have studied those determinations having an exactly solvable eigenvalue problem. On the present study, we also included some comments in relation to determinations not exactly solvable from this point of view. In addition, the use of self adjoint determinations for which the ground state wave function has nodes (zeroes) produces formal SUSY partners with a finite number of eigenvalues or even with a purely continuous spectrum. We give some worked examples of these situations.

1. Introduction

One dimensional quantum models are often useful for a better understanding of quantum systems. The infinite one dimensional square well (or particle of mass m in a one-dimensional box) is characterized by an infinite value of the potential outside a finite interval, in which this potential is zero. Elementary books in Quantum Mechanics give an exact solution to the eigenvalue problem (Schrödinger equation with boundary conditions) provided by this model, showing a countably infinite number of energy levels given by the formula E n = ( n π ) 2 / 2 m L 2 , n = 1 , 2 , 3 , , where m has the dimensions of a mass and L is the length of the interval in which the potential is equal to zero. These energy levels has been determined using the boundary condition that the solutions of the Schrödinger equation 2 2 m d 2 ϕ ( x ) d x 2 = E ϕ ( x ) vanishes at the edges of the interval with zero potential. We call this solution the textbook solution.
From the mathematical point of view, this is the same that solving the eigenvalue problem H ψ = E ψ with H = d 2 / d x 2 as an operator on the Hilbert space L 2 [ a , a ] , where we have chosen L = 2 a and 2 / 2 m = 1 for simplicity. So far, we always write D 2 to denote the operator d 2 / d x 2 on the Hilbert space L 2 [ a , a ] , unless otherwise stated.
The choice of the interval is really irrelevant under the condition that it be finite. If we demand that the solutions vanish inside the region where the potential is infinite and assume that they should be continuous at all points, which implies that they should be zero at the edges (in our case at x = ± a ), we obtain the textbook solution for the energy levels. This is all well known.
However, this is not the whole story. Recall that a Hamiltonian should be a self adjoint operator (we discard other considerations such that pseudo-hermiticity or PT-symmetric Hamiltonians). In addition, although we always should require to the wave functions to be zero inside the region of infinite potential, there is no clear reason why they should be continuous at all points. They have to be square integrable and this allows for finite jumps of the wave function at some points. Thus, we may consider wave functions which do not vanish at the edges x = ± a .
The operator D 2 : = d 2 / d x 2 on the Hilbert space L 2 [ a , a ] is an unbounded operator, which in particular means that it is not defined on all functions of the Hilbert space, but instead on all functions belonging to the dense subspace D of all absolutely continuous square integrable functions on [ a , a ] with absolutely continuous first derivative and whose second derivative, which is well defined except possibly on a set of zero Lebesgue measure, is also square integrable on the same interval. Let us also consider the following scalar product of functions f ( x ) , g ( x ) L 2 [ a , a ] , D : = d / d x , and use integration by parts:
f | ( D 2 ) g = a a f * ( x ) ( D 2 g ) ( x ) d x = { f * ( a ) g ( a ) f * ( a ) g ( a ) } + { f * ( a ) g ( a ) f * ( a ) g ( a ) } + ( D 2 ) f | g .
Now, the objective is to characterize all self adjoints determinations (also called extensions [1], both expressions will be used indistinctly along the paper) of the differential operator D 2 on L 2 [ a , a ] . To do it, we need to characterize their respective domains, which are dense subspaces of L 2 [ a , a ] contained in D . The procedure is well know, since functions f ( x ) belong to one of these domains if and only if they fulfil the following relation:
{ f * ( a ) g ( a ) f * ( a ) g ( a ) } { f * ( a ) g ( a ) f * ( a ) g ( a ) } = 0 .
The terminology of self adjoint determinations or self adjoint extensions of a symmetric operator, in particular the differential operator D 2 , is explained in Appendix C.
As proven in [2], this condition may be written in matrix form as f ( x ) L 2 [ a , a ] :
2 a f ( a ) i f ( a ) 2 a f ( a ) + i f ( a ) = U 2 a f ( a ) + i f ( a ) 2 a f ( a ) i f ( a ) ,
where U is an arbitrary 2 × 2 unitary matrix. In fact, according to a Theorem due to von Neumann [3], each unitary 2 × 2 matrix gives different choices of domains for D 2 in which D 2 is self adjoint. Each of the 2 × 2 unitary matrices is characterized by the particular choice of a set of four independent real parameters. The most general form of these matrices is [2]:
U = e i ψ m 0 i m 3 m 2 i m 1 m 2 i m 1 m 0 + i m 3 ,
with ψ [ 0 , π ] , m 0 2 + m 1 2 + m 2 2 + m 3 2 = 1 . Note that there are four independent parameters, ψ and three out of the m i , where the dependent parameter may be conveniently chosen. Each of the choices of these parameters gives relations between the boundary values of functions and their derivatives and these relations determine the domains of each self adjoint determination of D 2 . A pedagogical presentation of the self-adjoint determinations of D 2 is given in [2]. All them have some features in common such as having a purely discrete spectrum with an infinite number of eigenvalues and eigenfunctions, as happens for the particular case studied in textbooks. This fact is shown in [4], vol 2, p. 90.
The study of supersymmetric (SUSY) partners of a given Hamiltonian has been an object of study during the past few decades. As a result, many articles have been published. For the benefit of the reader, we have selected to list in the bibliography just a few [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. This kind of supersymmetry intends the construction of a sequence of Hamiltonians, derived from an original Hamiltonian, with a closely related energy spectrum, which could be obtained from the spectrum of this original Hamiltonian.
In a previous paper [20], we have obtained sequences of supersymmetric partners from a wide class of these self adjoint determinations of D 2 . In the studied cases, being given a self adjoint determination of D 2 on L 2 [ a , a ] , we obtain a sequence of partners with the property that each of the Hamiltonians of the sequence loses the ground state of the previous one; all the rest of the spectrum remains the same. All potentials in one and the same sequence have the form D 2 + V , where V is a potential obtained with the ground state of the previous Hamiltonian using a well known procedure [17] and that we describe in Appendix A. In all cases studied in [20], the energy of the ground state of the self-adjoint determination of D 2 used is positive, no degenerate and its corresponding wave function has no nodes, i.e., no zeroes.
In spite of its appearance, not all self adjoint determinations of D 2 have ground state with positive energy. Some determinations have ground state with zero energy and, most surprising, with negative energy, so that these determinations are not positive definite operators. In both cases, a result due to Naimark [4] shows that we may have either one (possible with double degeneracy) or two negative energy states, or a ground state with negative energy and a first excited state with zero energy, or a zero energy ground state, possibly with double degeneracy. No other possibility exists. Since the study of SUSY partners for these determinations is somehow different than for the rest, we call exceptional cases to this situations. The category of exceptional cases also include those determinations for which the ground state wave function has a node (zero). Their first formal partner may have a completely different properties than the original Hamiltonian, for instance it may have a finite number of eigenvalues or even a purely continuous spectrum. By formal partner, we mean a first SUSY partner constructed following the general rules as explained in Appendix A. To close the category of exceptional cases, we include those which are not exactly solvable in the sense that we cannot give for the first partner either the exact values of all eigenvalues or their respective wave functions.
The objective of the present article is the determination and study of the first partner in these exceptional cases. We should note that this determination means the construction of the potential of the Hamiltonian partner.
The present article is organized as follows: in Section 2, we summarize the main results given in [20] and makes the present article self contained. This summary will help us to understand the problem under consideration. We say in passing that we have corrected here an imprecision in [20]. In Section 3, we make a first attempt to classify the self adjoint extensions with ground state with zero or negative energy including degeneracies, such that the wave functions of these states still have symmetry properties. In Section 4, we construct their first or second SUSY partner. Note that for degenerate ground states a different technique should be used. Finally, on a brief Section 5, we give an example on how the absence of the mentioned symmetries for the wave function of the ground state may give a non-Hermitian first SUSY partner. The paper closes with three Appendices including the construction of first and second SUSY partners, with or without degeneracy on the ground state and one more on extensions of symmetric operators. This is known material that we add for the benefit of the reader.

2. Regular Cases

By regular cases, we mean the study of the SUSY partners of those self adjoint determinations of the operator D 2 on L 2 [ a , a ] (original Hamiltonian) which are exactly solvable and have strictly positive eigenvalues only. In addition to positivity, we add the extra conditions m 2 = m 3 = 0 in order to assure the solvability of both the eigenvalue problem associated to the original Hamiltonian and the analytic determination of all SUSY partners at all orders. Along this section, by positivity we mean strict positivity, in the sense that all eigenvalues are strictly positive.
In the introduction of Formula (4), we have given an operator U that characterize all self adjoint determinations of D 2 on L 2 [ a , a ] . This operator gives the boundary conditions that give the different domains of all self-adjoint determinations. Now, it is convenient to provide of a notation for the domain defined by U as D U , so that D U is the domain of the self adjoint determination of D 2 fixed by U.
Why is it possible to have determinations of D 2 with negative energy? Assume that ϕ ( x ) D ( U ) . Then,
ϕ | D 2 ϕ = i D ϕ | i D ϕ + ϕ * ( a ) ϕ ( a ) ϕ * ( a ) ϕ ( a ) .
Observe that the right hand side of (5) may be negative for some ϕ ( x ) D U . We give some examples in Section 3.
Regular cases have been studied in [20] and it is interesting to give a summary of the obtained results here. First of all, we have to characterize those self adjoint determinations of D 2 on L 2 [ a , a ] corresponding to the regular cases. Since all eigenvalues are positive, their general form should be E = s 2 / ( 2 a ) 2 > 0 , where s is a real discrete parameter. Determinations of D 2 satisfying this condition are strictly positive. Recall that all extensions have an infinite number of eigenvalues. These eigenvalues must satisfy the time independent Schrödinger equation D 2 ϕ ( x ) = E ϕ ( x ) . The general form of the eigenfunctions for all self adjoint extensions should be:
ϕ ( x ) = A cos s x 2 a + B sin s x 2 a ,
where the coefficients A and B are determined by the boundary conditions. Now, let us consider strictly positive determinations. Using (6) in (4), we obtain the following system of linear equations on the parameters A and B:
L ( s ) U M ( s ) A B = N ( s ) A B = 0 ,
where N ( s ) = ( L ( s ) U M ( s ) ) , with
L ( s ) : = s sin s 2 i cos s 2 s cos s 2 + i sin s 2 s sin s 2 + i cos s 2 s cos s 2 + i sin s 2 ,
and
M ( s ) = s sin s 2 + i cos s 2 s cos s 2 i sin s 2 s sin s 2 i cos s 2 s cos s 2 + i sin s 2 .
Since det N ( s ) = 0 , otherwise the solutions of (7) would have been A = B = 0 , one of the eigenvalues of N ( s ) must be zero and the other the trace of the matrix N ( s ) , Tr [ N ( s ) ] . The determinant of N ( s ) is easy to calculate and, after some algebra, the expression det N ( s ) = 0 yields:
s sin s = m 1 cos s sin ψ m 0 cos ψ ± m 1 cos s sin ψ m 0 cos ψ + m 0 + cos ψ m 0 cos ψ sin 2 s .
These are two transcendental equations, one with plus sign and the other with minus sign, which give the values of s and, therefore, the values of the energy in terms of the parameters ψ , m 0 and m 1 . This expression is valid for all positive determinations of D 2 including those with m 2 0 and m 3 0 . However, the energy values do not depend on m 2 and m 3 .
We also want to obtain the values for A ( s ) and B ( s ) in order to obtain the wave function which corresponds to the eigenvalue E given by s, through E = s 2 / ( 2 a ) 2 . Although this coefficients have been obtained in [20], we recalculate them here as the result given in [20] was not fully correct. One interesting point is that, no matter if the energy levels do not depend on m 2 0 and m 3 0 , the wave functions do depend on these parameters in general. To begin with, let us write the general form for the eigenvalues of the matrix N as:
λ ± = Tr N 2 ± Tr N 2 2 .
Since N is a complex matrix, its zero eigenvalue is not necessarily λ as claimed in [20]. Instead, the zero eigenvector is either λ + or λ depending whether the argument of Tr N lies on the left, { z C ; Real z < 0 } , or the right, { z C ; Real z > 0 } , complex half-planes, respectively. Then, if we write the matrix N in simplified form as:
N = a b c d ,
the eigenvectors ( A ( s ) , B ( s ) ) of N ( s ) with zero eigenvalue should have the form,
( A ( s ) , B ( s ) ) = ( d , c ) , π 2 < arg ( Tr N ) < π 2 , ( a , c ) , π 2 < arg ( Tr N ) < 3 π 2 .
This gives the following result, which has been corrected respect the result given in [20]:
A ( s ) = s cos s 2 1 + e i ψ ( m 0 i m 1 + m 2 + i m 3 ) i sin s 2 1 + e i ψ ( m 0 i m 1 + m 2 + i m 3 ) π 2 < arg Tr N < π 2 , i cos s 2 1 + e i ψ ( m 0 + i m 1 + m 2 i m 3 ) + s sin s 2 1 e i ψ ( m 0 + i m 1 + m 2 i m 3 ) , π 2 < arg Tr N < 3 π 2
and
B ( s ) = s 1 + e i ψ ( m 0 + i m 1 + m 2 i m 3 ) sin s 2 + i + e i ψ ( i m 0 m 1 + i m 2 + m 3 ) cos s 2 .
From these results, one may study the behaviour of the self adjoint extensions preserving positivity, i.e., those with a ground state with positive energy. Partial results may be found in [2,20]. The form for these coefficients with m 2 = m 3 = 0 is obvious.
In addition to the angular parameter ψ , the relation i = 0 3 m i 2 = 1 suggests the possibility of writing all other parameters in terms of some angular variables θ i , i = 0 , 1 , 2 as follows:
m 0 = cos θ 1 cos θ 0 , m 1 = cos θ 1 sin θ 0 , m 2 = sin θ 1 cos θ 2 , m 3 = sin θ 1 sin θ 2 .
We have shown that the energy levels of those (strictly) positive extensions of D 2 are solely determined by the values of three parameters: ( ψ , m 0 , m 1 ) . Alternatively, they depend on three angular variables only, ( ψ , θ 0 , θ 1 ) , as we may infer from (16). However, due to the complicated form of the transcendental equation (10), it is not possible in general to obtain an explicit expression of the form s = s ( ψ , θ 0 , θ 1 ) from (10) and (16). This explicit expression would have been very useful in order to obtain the spectrum of all the considered self adjoint extensions. Nevertheless, in some particular cases, simple transcendental equations can be given, which have been studied in [20]. For instance, if the wave equation for the ground state, ϕ ( x ) , is parity invariant, which means | ϕ ( x ) | 2 = | ϕ ( x ) | 2 . This happens if and only if at least one of the following three relations hold:
m 3 = 0 ,
sin s = 0 ,
s sin s = m 2 + cos s cos ψ m 3 sin ψ ± m 2 + cos s cos ϕ m 3 sin ψ + m 3 + sin ψ m 3 sin ψ sin 2 s .
Before going on, it is important to say that there exists a group of extensions for which the solution, ψ ( x , t ) , of the time dependent Schrödinger equation is time reversal invariant in the sense that ψ ( x , + t ) = ψ ( x , t ) . These are characterized by m 2 = 0 .
Then, the three relations (19) shows that these tractable cases may be grouped into three clases:
  • Extensions with parity and time reversal invariance. These are those for which m 2 = m 3 = 0 and, therefore, m 0 2 + m 1 2 = 1 , so that m 0 = cos θ 0 and m 1 = sin θ 0 . Now, the transcendental Equation (10) becomes much simpler and the sign depends on the possible values of s. Thus, the sign is positive if s lies on the intervals ( 0 , π ) , ( 2 π , 3 π ) , and is negative if s belongs to the intervals ( π , 2 π ) , ( 3 π , 4 π ) , . In both cases, the transcendental equations given the values of the energy are, respectively,
    s tan s 2 = cot ψ + θ 0 2 , plus sign in ( 10 )
    s cot s 2 = cot ψ θ 0 2 , minus sign in ( 10 ) .
    Note that ψ + θ 0 2 and ψ θ 0 2 are two independent angles varying from 0 to π / 2 . The solutions in s, s n , n = 0 , 1 , 2 , of (20) and (21), give the energy levels, E = s n 2 / ( 2 a ) 2 . The respective eigenfunctions are:
    ϕ n ( x ) = cos s n x 2 a , n = 0 , 2 , 4 , , ϕ n ( x ) = sin s n x 2 a , n = 1 , 3 , 5 , .
  • Parity preserving extensions fulfilling sin s = 0 . The energy levels are E n = ( n 2 π 2 ) / ( 2 a ) 2 , n = 1 , 2 , , which correspond to the energy levels obtained for the extension described in textbooks.
  • Parity preserving extensions fulfilling (19).
In [20], we have obtained sequences of SUSY partners for these extensions. It is time to drop the positivity condition and consider more general situations and this is precisely the goal of the present article, to be discussed next.

3. Transcendental Equation for Other Cases

In addition to positive extensions, which are those with ground state with strictly positive energy, two other situations are possible" either the ground state has zero energy or negative energy. In the latter case, as mentioned in the Introduction, there are one or two negative energy levels, so that E = s 2 / ( 2 a ) 2 implies s = i r with r real for these levels. Zero energy means one solution with s = 0 .
Then, let us go back to the matrix N in (12). This matrix has different forms depending on the three cases considered, which are: (i) Positive ground state energy E > 0 , N E > 0 ( s ) ; (ii) Negative ground state energy, E < 0 , N E < 0 ( r ) ; (iii) Zero ground state energy, E = 0 , N E = 0 . These matrices are rather complicated and their explicit form is not very interesting. The condition that their respective determinants are equal to zero gives the following respective transcendental equations:
(i) For E = s 2 > 0 , det N E > 0 ( s ) = 0 gives:
( m 0 + cos ψ ) sin s + 2 ( m 1 cos s sin ψ ) s + sin s ( cos ψ m 0 ) s 2 = 0 ;
(ii) For E = r 2 < 0 , det N E < 0 ( r ) = 0 gives,
( m 0 + cos ψ ) sinh r + 2 ( m 1 cosh r sin ψ ) r + sinh r ( m 0 cos ψ ) r 2 = 0 ;
(iii) For E = 0 , det N E = 0 = 0 gives:
m 0 2 m 1 + cos ψ 2 sin ψ = 0 .
While (23) gives all possible positive values of the energy in terms of s, Equation (24) gives the possible values with negative energy (two maximum). Equation (25) gives the relation that the parameters should satisfy for the existence of a ground state with zero energy.
So far, we have concentrated ourselves in positive extensions or determinations of D 2 on L 2 [ a , a ] and the calculation of their SUSY partners. After we drop the positivity condition, we have to select solvable extensions in the sense that we may, at least, find the explicit for for the wave function of the ground state. As mentioned earlier, this wave function is essential in order to obtain the first partner.
The “solvable” cases object of our next study are those preserving time reversal invariance and parity inversion, with the previously discussed senses. This implies that at the same time our determinations of D 2 verify two conditions: (i) m 2 = 0 , for the time reversal invariance and (ii) any or several of (17)–(19).
In terms of the new independent angular parameters ( θ 0 , θ 1 , θ 2 , ψ ) as defined in (16), we are considering extensions with θ 1 = θ 2 = 0 , so that these extensions are labelled by ( θ 0 , 0 , 0 , ψ ) .
Now, transcendental equations given the energy values are quadratic on the square root of the energy, which is given by s for the positive energy solutions and by r for the negative energies. Transcendental equations may be factorized into two equations on the respective variables s and r. They have even or odd ground state wave functions, respectively.
In the next table, we we list these transcendental equation with their respective ground state wave function.
In Table 1, we have obtained the transcendental equation given the negative energy eigenvalues (two maximum) by making the change s i r , as previously mentioned. It is convenient the use of graphics in order to describe the behaviour of the solutions of the above transcendental equations.
In Figure 1, we observe the zone of negative energies from the plane s = 0 . We have divided the figure into several regions. Orange and blue color represent the values of ( θ 0 , ψ ) for which there are negative eigenvalues of even or odd eigenfunctions, respectively. In the absence of color (white) the values of ( θ 0 , ψ ) determine a ground state with positive energy. In the regions (a) and (b), we have only one state with negative energy, which is the ground state. Contrarily, regions (d) and (c) display two eigenvalues with negative energy. In the case of (d) the orange color superposes the blue color. This means that the ground state is odd and the first excited state is even. The opposite is true in the region (c). Dashed lines correspond to ground states with zero energy and the red line shows negative energy doubly degenerated.
In Figure 2, we include a three-dimensional picture, where the third coordinate corresponds to the (negative) energy. It is clear here that the degenerate levels are in the intersection of the two latter surfaces on Table 1. Each of these surfaces are written on the form F ( ψ , θ 0 , r ) = 0 .
We give a summary of these results on Table 2. In the first column, a label denotes the corresponding regions in Figure 1. The last column gives the signs of the energy of the ground state and the first few excited states corresponding to the classes of self adjoint extensions described on columns 2 and 3.
Note that the parity of the first excited state having positive energy may be either odd or even, independently of the parity of the ground state. Meanwhile, if the first excited state has negative or zero energy, the parity of its wave function is the opposite to the parity of the ground state wave function.
In Figure 3, we depict the regions in the space of parameters ( θ 0 , θ 1 , ψ ) for which there exists a ground state with negative energy. The empty regions give the values of the parameters for those extensions with ground state with positive energy. The black net covers those surfaces for which there exists zero energy levels, characterized by m 0 + 2 m 1 + cos ψ 2 sin ψ = 0 .
The black dashed lines in Figure 1 mean the intersection of the plane E = 0 with the two last surfaces in the mid column of Table 1. The intersection of these two surfaces give the degeneration of the negative energy levels. The equations on the plane E = 0 of these dashed lines are as follows:
(i) The intersection of E = 0 with the forth surface in Table 1 is given by
ψ = θ 0 + 2 arccot 2 .
(ii) The intersection of E = 0 with the third surface in Table 1 is given by
ψ = θ 0 + π .
(iii) There is only one situation in which there exists one degenerate level with zero energy. This is given by the following value of the parameters:
( θ 0 , ψ ) = π 2 arccot 2 , π 2 + arccot 2 .
The corresponding even and odd eigenfunctions are given as the third and forth entry in the third column of Table 1, respectively.
When the ground state is doubly degenerate (28), we may also obtain the remaining values of the spectrum, which are all positive. For even and odd wave functions, the values of the energies are respectively given by the following transcendental equations:
s sin s 2 = 0 , s cos s 2 2 sin s 2 = 0 .
Obviously, both equations have a solution for s = 0 , hence the double degeneracy of the ground level. The solutions of the first equation in (29) are obvious and of the second are given by the intersections of cot ( s / 2 ) with 2 / s .
In Figure 4, we show the values of the parameters ψ and θ 0 for which the first positive excited state has wave functions with parity odd (blue) or even (orange). The point at which the dotted lines cross give the degenerate level with zero energy. Negative energy levels have been depicted in Figure 1.
In Figure 5, we represent the surfaces f ( ψ , θ 0 , s ) = 0 given the states with positive energy with wave equations with parities even (orange) or odd (blue). Contrary to the first energy levels classified in Table 2, there are plenty subcategories, so that the equivalent of Table 2 in this case would have been extremely complicated. Note that Figure 4 is nothing else than the projection into the plane E = 0 of Figure 5.
Concerning solutions with ground state with positive energy, we arrive to the following conclusions for spectra and eigenfunctions of extensions in terms of the parameters ( θ 0 , ψ ) :
  • The energy levels and the eigenfunctions for the extensions characterized by the values ( θ 0 , ψ ) = ( 0 , 0 ) , ( θ 0 , ψ ) = ( 2 π , 0 ) and ( θ 0 , ψ ) = ( π , π ) coincide with those given in textbooks:
    E n : = n π 2 a 2 , ϕ n ( x ) = cos n π x 2 a , n = 1 , 3 , 5 , sin n π x 2 a , n = 2 , 4 , 6 ,
    with = 2 m = 1 . We call this result the textbook solution.
  • There are two extensions with a double degenerate positive spectrum, which are given by the values θ 0 , ψ = π 2 , π 2 and θ 0 , ψ = 3 π 2 , π 2 . The former has a ground state of odd parity and its energy levels are:
    E n = n 2 π 2 a 2 , n = 0 , 1 , 2 , ,
    with respective wave functions:
    ϕ 0 ( x ) = 3 2 a 3 x , n = 0 , ϕ n ( x ) = cos n π a x and sin n π a x , n = 1 , 2 , .
    From (32), we observe that excited states are doubly degenerate. The extension given by θ 0 , ψ = 3 π 2 , π 2 has a ground state with even parity with E 0 = 0 . The remaining energy values are
    ϕ ( x ) = 1 2 a , n = 0 , ϕ n ( x ) = cos ( 2 n 1 ) π 2 a x and sin ( 2 n 1 ) π 2 a x ,
    with n = 0 , 1 , 2 , . Again we observe a double degeneracy for all states save for the ground state.
  • In general, it is not possible the determination of all the values of the energy, since they come after the solution of a transcendental equation. Using Mathematica or another numerical method, we may obtain a finite number of eigenvalues for each extensions, but no more. On the other hand, once we have the value of the energy, we have the corresponding eigenfunction. Nevertheless the higher the value of n the closer the energy levels are to the levels of the textbook solution given in (30).

4. Supersymmetric Partners

Along the present Section, we classify the SUSY partners of those self adjoint extensions having an exactly solvable eigenvalue problem with zero or negative energy ground state, either degenerate or not. The general procedure for the construction of SUSY partners of a given self adjoint determination of D 2 is discussed in Appendix A. However, we have to consider the case in which this ground state has a double degeneracy. The construction of the first SUSY partner in this situation is somehow different as this construction works for the non degenerate ground state, as described in Appendix B. Ground states for exactly solvable determinations of D 2 have a fixed parity.
In order to comply with the usual notation, let us use H ( 0 ) : = D 2 from now on, in order to denote an arbitrary self adjoint extension.
For the construction of SUSY partners of the self adjoint extensions that are the object of our attention here, we begin with a self adjoint extension with doubly degenerate ground state. The general formalism should be adapted to this new situation as discussed in the sequel.

4.1. Second Order Confluent SUSY

In Appendix A, we have sketched the construction of the first SUSY partner of a given self adjoint operator. However, if the condition that the ground state of the operator H ( 0 ) is doubly degenerate applies, we have to use a modification of the mentioned construction. As customary, we start with a self adjoint extension H ( 0 ) with spectrum { E 0 ( 0 ) , E 1 ( 0 ) , E 2 ( 0 ) , } and respective eigenfunctions { ϕ 0 ( 0 ) ( x ) , ϕ 1 ( 0 ) ( x ) , ϕ 2 ( 0 ) ( x ) , } . Then, we assume that E 0 ( 0 ) = E 1 ( 0 ) , so that the ground state is doubly degenerate. Thus, the functions ϕ 0 ( 0 ) ( x ) and ϕ 1 ( 0 ) ( x ) are different eigenfunctions of E 0 ( 0 ) , which may be taken as mutually orthogonal. Then, we use the second order confluent SUSY formalism introduced in [21]. The crucial point lies in the use of a generalized wave function defined with the help of any of the orthogonal eigenfunctions of E 0 ( 0 ) :
v j ( x ) = u j ( x ) w 0 u j 2 ( x ) d x u j 2 ( x ) d x .
In (34), the function u j ( x ) is either ϕ 0 ( 0 ) ( x ) or ϕ 1 ( 0 ) ( x ) , we may choose one of them, so that j = 0 , 1 .
The constant w 0 is chosen so as to avoid the presence of nodes in the function v j ( x ) [21]. As seen in Appendix A, first SUSY partners are well defined for those extensions for which ground states have no zeros, also called nodes. In Appendix B, we have discussed the difficulties for a proper consideration of the first order SUSY partner of those extensions with ground state with degeneracy, so that we have to go directly to the second order partner and skip the first one.
Thus, the second order SUSY partner of H ( 0 ) is given by [21,22]:
H ( 2 , 2 ) = x 2 + V ( 2 , 2 ) , V ( 2 , 2 ) = V ( 0 ) 2 x 2 log W ( v j ( x ) , ϕ j ( 0 ) ( x ) ) , j = 0 , 1 ,
with x : = d / d x , and W denotes the Wronskian. The eigenfunctions are:
ϕ n ( 2 , 2 ) ( x ) = W ( v j ( x ) , ϕ j ( 0 ) ( x ) , ϕ n ( 0 ) ) ( x ) W ( v j ( x ) , ϕ j ( 0 ) ( x ) ) , j = 0 , 1 , n 2 .
Then, we list the second order SUSY partners of the extensions with negative or zero ground state degenerate or not, which are exactly solvable. We classify these partners using the values of the parameters that fix the extension, with the additional criterium whether the ground state is either degenerate or not. Here, we do not go beyond the second partner, since the study and classification of higher order partners goes as previously studied in [20]. All the extensions depend on the parameters ( cos θ 0 , sin θ 0 , 0 , 0 , ψ ) .

4.1.1. Degenerate Ground State with Zero Energy

The only exactly solvable extension having a degenerate ground state with zero energy is characterized by the following values of the parameters:
( θ 0 , ψ ) = π 2 + arccot 2 , π 2 arccot 2 ,
being the remainder parameters equal to zero. The two dimensional eigenspace of functions with energy zero is spanned by one odd and one even eigenfunctions. The explicit values of the energy for excited states depend on whether they have even parity, with energy levels E 2 n , n > 1 , or odd parity, with energy levels E 2 n + 1 . For even levels, E 2 n is proportional to ( 2 n ) 2 π 2 , while for odd levels their energy is given by a transcendental equation. For high values of n, these levels go to E 2 n + 1 ( 2 n + 1 ) 2 π 2 . The list of the eigenfunctions is:
ϕ 0 ( 0 ) ( x ) , ϕ 1 ( 0 ) ( x ) , ϕ 2 ( 0 ) ( x ) , ϕ 3 ( 0 ) ( x ) , = 1 2 a , x , cos s 2 x 2 a , sin s 3 x 2 a , ,
where we have omitted the normalization constants. Their corresponding values of the energy are:
E 0 ( 0 ) , E 1 ( 0 ) , E 2 ( 0 ) , E 3 ( 0 ) , = 0 , 0 , s 2 2 a 2 , s 3 2 a 2 , ,
with s 2 n = 2 n π . The eigenfunctions are displayed at Figure 6, and the deviation of the odd labeled levers from the even are given by Figure 7.
We recall that in Appendix B, we have justify the reasons for going directly to the second order SUSY partner, avoiding the consideration of the first order partner. There are two second order partners that are constructed with the aid of the functions v i ( x ) , i = 1 , 2 , defined in (34). The function v 0 ( x ) comes after (34) and the even eigenfunction of the ground state, ϕ 0 ( 0 ) ( x ) , while v 1 ( x ) is the result of using in (34) the other eigenfunction of the ground state, i.e., ϕ 1 ( 0 ) ( x ) .
Let us find the second order partner for the even eigenfunction ϕ 0 ( 0 ) ( x ) . After (34) and (38), the function v 0 ( x ) is:
v 0 ( x ) = 2 a w 0 x x 2 2 .
Following the procedure explained in Appendix B, see also [21,22], we obtain the following second order eigenfunctions:
ϕ n ( 2 , 2 ) ( x ) W v 0 ( x ) , ϕ 0 ( 0 ) ( x ) , ϕ n ( 0 ) W v 0 ( x ) , ϕ 0 ( 0 ) ( x ) = x ϕ n ( 0 ) 2 a w 0 x E n ϕ n ( 0 ) .
The second order potential is:
V ( 2 , 2 ) = 2 d d x 2 W v 0 , ϕ 0 ( 0 ) = 2 ( x 2 a w 0 ) 2 ,
where we choose w 0 to avoid the singularity in the denominators of (41) and (42). Using the change on the variable given by x = ( t + 1 ) 2 a w 0 , the Schrödinger equation for the partner H ( 2 , 2 ) is transformed so as to become:
d 2 d t 2 ϕ n ( 2 , 2 ) ( t ) + 2 t 2 ϕ n ( 2 , 2 ) ( t ) = ( 2 a w 0 ) 2 E n ϕ n ( 2 , 2 ) ( t ) .
This is a Bessel equation with index 3 / 2 , for which the general solution after undoing the change of variables is:
ϕ n ( 2 , 2 ) ( x ) = C 1 x 2 a w 0 1 J 3 2 s n x 2 a w 0 + C 2 x 2 a w 0 1 J 3 2 s n x 2 a w 0 ,
where C 1 and C 2 are arbitrary constants, which should be fixed using (41). The spectrum of H ( 2 , 2 ) coincides with the spectrum of H ( 0 ) , although it has lost the degeneracy in the first level: 0 , E 1 , E 2 , E 3 , .
If we begin with the odd ground eigenfunction, ϕ 1 ( 0 ) ( x ) , we have to use the following function as given in (34):
v 1 ( x ) = 4 a 3 w 0 + x 3 ,
while the potential is
V ( 2 , 2 ) ( x ) = 2 d d x 2 W { v 1 , ϕ 1 ( 0 ) } = 3 4 a 3 w 0 x + x 4 x 3 2 a 3 w 0 2 ,
where | w 0 | > 1 / 2 in order to avoid singularities. The eigenfunctions of H ( 2 , 2 ) are proportional to:
ϕ n ( 2 , 2 ) ( x ) 1 E n W { v 1 , ϕ 1 ( 0 ) , ϕ n ( 0 ) ( x ) } W { v 1 , ϕ 1 ( 0 ) } = 3 x 2 x ϕ n ( 0 ) ( x ) 3 x ϕ n ( 0 ) ( x ) x 3 2 a 3 w 0 E n ϕ n ( 0 ) ( x ) ,
with energy levels 0 , E 1 , E 2 , E 3 , . The potential V ( 2 , 2 ) is a quotient of polynomials and the resulting Schrödinger equation has three singular regular points located at the cubic roots of 2 w 0 a 3 and a singular irregular point at the infinity, resulting of a confluence of two points. Thus, the Schrödinger equation is a Fuchsian equation of five points, taking into account the confluence of two of these points at the infinity.
We change the variable x into the variable z as x = ( 2 w 0 ) 1 / 3 a z and the function solution, with a new one y ( z ) given by:
ϕ n ( 2 , 2 ) ( x ) = e τ x / 2 ( z z 1 ) γ 1 ( z z 2 ) γ 2 ( z z 3 ) γ 3 y ( z ) .
After these changes, we arrive at the following Fuchs differential equation:
d 2 y ( z ) d z 2 + Γ + γ 1 z z 1 + γ 2 z z 2 + γ 3 z z 3 d y ( z ) d z + q 1 z z 1 + q 2 z z 2 + q 3 z z 3 y ( z ) = 0 ,
where the constants are given by:
γ 1 = 2 , γ 2 = 4 , γ 3 = 4 , q 1 = 4 2 i ϵ , q 2 = 2 ( 2 + ( 3 + i ) ϵ ) 1 3 1 , q 3 = 4 i ( ( 3 3 i ) ϵ + i 3 3 ) 3 3 i , Γ = 2 i ϵ , ϵ = 2 3 / 2 a 2 w 0 2 / 3 E n ( 2 ) , z 1 = 1 , z 2 = ( 1 ) 1 / 3 , z 3 = ( 1 ) 2 / 3 .
In Figure 8, we depict the first five energy levels for H ( 2 , 2 ) when the ground state of H ( 0 ) is doubly degenerate and we use the confluent SUSY using ϕ 1 ( 0 ) as initial eigenfunction. In Figure 9, we depict the same energy levels under similar conditions, this time using ϕ 0 ( 0 ) as initial eigenfunction.

4.1.2. Degenerate Ground State with Negative Energy

Now, we intend to reproduce the previous analysis when the extension of H ( 0 ) has a doubly degenerate ground state with negative energy equal to E 0 ( 0 ) = r 0 2 / ( 4 a 2 ) , r 0 > 0 . There are several extensions with this property characterized by some values of the parameters, which are given by the red curves that appear in Figure 1. The eigenfunctions and the corresponding eigenvalues for these extensions are:
ϕ 0 ( 0 ) ( x ) , ϕ 1 ( 0 ) ( x ) , ϕ 2 ( 0 ) ( x ) , ϕ 3 ( 0 ) ( x ) , = cosh r 0 x 2 a , sinh r 0 x 2 a , cos s 2 x 2 a , sin s 3 x 2 a , ,
and
E 0 ( 0 ) , E 1 ( 0 ) , E 2 ( 0 ) , E 3 ( 0 ) , = r 0 2 a 2 , r 0 2 a 2 , s 2 2 a 2 , s 3 2 a 2 , ,
where the values of r 0 and s n , n > 2 , are given by transcendental equations as explained in Section 3.
Since the ground state is doubly degenerate, we obtain directly the second order partner using the technique described in the precedent subsection. To implement this technique, we begin with the even eigenfunction for the ground state ϕ 0 ( 0 ) = cosh r 0 x 2 a and determine its associated function following the recipe given in (34). This associate function is:
v 0 ( x ) = ( 2 w 0 x ) sinh r 0 x 2 a .
Then, the potential for the second order partner Hamiltonian H ( 0 , 0 ) is given by
V ( 2 , 2 ) ( x ) = 2 x 2 log W v 0 , ϕ 0 ( 0 ) = 2 r 0 2 r 0 ( 2 w 0 x ) sinh r 0 x a + 2 a cosh r 0 x a + 2 a a a sinh r 0 x a + r 0 ( x 2 w 0 ) 2 .
As always, the choice of w 0 intends to avoid singularities.
We may have began instead with the odd eigenfunction for the ground state, ϕ 1 ( 0 ) = sinh r 0 x 2 a . Then we would have the following associate function:
v 1 ( x ) = ( 2 w 0 + x ) cosh r 0 x 2 a .
This gives the following partner potential:
V ( 2 , 2 ) ( x ) = 2 x 2 log W v 1 , ϕ 1 ( 0 ) = 2 r 0 2 r 0 ( 2 w 0 + x ) sinh r 0 x a 2 a cosh r 0 x a + 2 a a r 0 ( 2 w 0 + x ) a sinh r 0 x a 2 ,
where the choice w 0 is arbitrarily determined in order to avoid singularities.

4.1.3. Other Solvable Cases Up to Second Order

First, we look for the second SUSY partner directly for two cases in which the ground state has non degenerate negative energy and zero energy for the first excited state, respectively. As customary, we just obtain the second partner since the technique to obtain successive partners is repetitive. We look for those extensions which we have named as “exactly solvable”, which are given by some determination of the parameters.
When the parameters θ 0 and ψ fulfil either one of the following relations:
ψ = θ 0 + 2 arccot ( 2 ) 2 π and ψ = θ 0 + 2 arccot ( 2 ) ,
then the self adjoint determination H ( 0 ) have a negative energy ground state with even eigenfunction and a first excited state with zero energy and odd eigenfunction, respectively. The eigenfunctions of H ( 0 ) are
ϕ 0 ( 0 ) ( x ) , ϕ 1 ( 0 ) ( x ) , ϕ 2 ( 0 ) ( x ) , ϕ 3 ( 0 ) ( x ) , = cosh r 0 x 2 a , x , cos s 2 x 2 a , sin s 3 x 2 a , ,
with respective eigenvalues
E 0 ( 0 ) , E 1 ( 0 ) , E 2 ( 0 ) , E 3 ( 0 ) , = r 0 2 a 2 , 0 , s 2 2 a 2 , s 3 2 a 2 , ,
where r 0 and s n , n > 2 depend on the parameters θ 0 and ψ through a transcendental equation.
The second partner potential is:
V ( 2 , 2 ) ( x ) = r 0 2 2 a 2 2 a 2 cosh r 0 x a + 2 a 2 + r 0 2 x 2 2 a 2 r 0 x sinh r 0 x 2 a 2 a cosh r 0 x 2 a 2 .
The regular formalism on the construction of SUSY partners asserts that the eigenvalues of H ( 2 , 2 ) are well determined, as are those of the original Hamiltonian, marked as H ( 0 ) , except for the two first in the list. This list in (59). However, the Schrödinger equation with potential (60) is rather complicated and so is the explicit for of its eigenfunctions. In any case, this explicit form is not especially relevant for our discussion. For the same reason, we cannot obtain H ( 3 , 2 ) and successive partners.
When the relation between the parameters is ψ = θ 0 + π , then the self adjoint extensions of the given H ( 0 ) have a non degenerate negative energy ground state with odd eigenfunction and a zero energy first excited state. Their eigenfunctions have the form:
ϕ 0 ( 0 ) ( x ) , ϕ 1 ( 0 ) ( x ) , ϕ 2 ( 0 ) ( x ) , ϕ 3 ( 0 ) ( x ) , = sinh r 0 x 2 a , 1 2 a , cos s 2 x 2 a , sin s 3 x 2 a , ,
with respective eigenvalues,
E 0 ( 0 ) , E 1 ( 0 ) , E 2 ( 0 ) , E 3 ( 0 ) , = r 0 2 a 2 , 0 , s 2 2 a 2 , s 3 2 a 2 , ,
which, again, depend on the values of the parameters θ 0 and ψ through a transcendental equation. Observe that the ground state of all these extensions has a node at the origin x = 0 . This means that a formal construction of the first SUSY partner will give a potential that either does not have eigenvalues or has a finite number thereof. These are the singular potential partners of which we provide some examples at the end of this Section. Nevertheless, the second order SUSY potential partner, V ( 2 , 2 ) , is regular and gives:
V ( 2 , 2 ) ( x ) = 2 x log W ( ϕ 0 ( 0 ) , ϕ 1 ( 0 ) ) = r 2 2 a 2 1 cosh r 0 x 2 a .
This potential is exactly solvable, as is a hyperbolic Pöschl-Teller potential. These eigenfunctions are:
ϕ n ( 2 , 2 ) ( x ) = C 1 P 1 i s n tanh r 0 x 2 a + C 2 Q 1 i s n tanh r 0 x 2 a = s n cos s n x 2 a + r sin s n x 2 a tanh r 0 x 2 a if n = 2 , 4 , 6 , s n sin s n x 2 a + r cos s n x 2 a tanh r 0 x 2 a if n = 3 , 5 , 7 , .
Next, let us consider the extension given by the values { θ 0 , ψ } = { 3 π / 2 , π / 2 } . This extension has a ground state with zero energy and excited states with double degeneracy. The eigenfunctions and the respective eigenvalues for this extension are respectively given by:
Eigenfunctions : 1 2 a cos ( 2 n 1 ) π 2 a x , sin ( 2 n 1 ) π 2 a x n = 1 , Spectrum : { 0 } ( 2 n 1 ) π 2 a 2 n = 1 .
We see that the eigenfunction for the ground state is a constant and, therefore, an even function. The remaining energy levels, labelled with the index n, are double degenerate as becomes clear from (69).
We may construct the first Hamiltonian partner, H ( 1 , 2 ) = d 2 / d x 2 + V ( 1 , 2 ) , for this extension, where the potential V ( 1 , 2 ) is:
V ( 1 , 2 ) = E 0 ( 0 ) 2 d 2 d x 2 log ϕ 0 ( 0 ) ( x ) 0 .
The rule giving this potential is presented in Appendix A and obeys to Formula (A3). Here, ϕ 0 ( 0 ) ( x ) and E 0 ( 0 ) are the eigenfunction of the ground state of H ( 0 ) and its energy value, respectively. Since ϕ 0 ( 0 ) ( x ) is a constant and E 0 ( 0 ) = 0 , this potential vanishes. The eigenfunctions and spectrum of H ( 1 , 2 ) are the same as in (65) with the exclusion of the first one, both for eigenfunctions ( 1 / 2 a ) as for eigenvalues (zero energy).
We observe that the ground state corresponds to n = 1 and is twice degenerate. Due to this double degeneracy, we may use the confluent SUSY technique here taken one of the eigenfunctions of the ground state. Thus,
u 0 ( x ) = ϕ 1 ( 1 , 2 ) ( x ) = cos π x 2 a , v 0 ( x ) = ( 2 w 0 x ) sin π x 2 a ,
V ( 3 , 2 ) = 2 π 2 π ( x 2 w 0 ) sin π x a + 2 a cos π x a + 2 a a a sin π x a + π ( x 2 w 0 ) 2 .
The potential V ( 3 , 2 ) produces a Schrödinger equation, which may be explicitly written. However, this equation is too complicated and does not match with any studied differential equation. We depict this potential in Figure 10. As usual, the implicit form of the wave functions is given by:
ϕ n ( 3 , 2 ) ( x ) = W v 0 , ϕ 1 ( 1 , 2 ) , ϕ n ( 1 , 2 ) W v 0 , ϕ 1 ( 1 , 2 ) .
Their final form is a bit too long and not quite interesting, so we omit it and just ass that the term w ( x ) : = W ( v 0 , ϕ 0 ( 0 ) ) always appears in the denominator.
Similarly, the extension given by the { θ 0 , ψ } = { π / 2 , π / 2 } has zero energy ground state and doubly degenerate excited states. Its eigenfunctions and respective eigenvalues are:
Eigenfunctions : 3 2 a 3 x cos n π a x , sin n π a x , Spectrum : { 0 } n π a 2 n = 1 .
Note that the bound state wave function is ϕ 0 ( 0 ) x , which is odd and has a node at the origin. This has some consequences. First of all, the first order SUSY Hamiltonian is singular at the origin:
H ( 1 , 2 ) = d 2 d x 2 2 x 2 log x = d 2 d x 2 + 2 x 2 .
The potential of this hamiltonian is displayed at Figure 11. The eigenfunctions of H ( 1 , 2 ) are the following:
ϕ n ( 1 , 2 ) ( x ) = W x , ϕ n ( 0 ) ( x ) x = x J 3 / 2 π n x a π n sin π n x a a + cos π n x a x , x J 3 / 2 π n x a π n cos π n x a a sin π n x a x .
These solutions including the Bessel functions J 3 / 2 and J 3 / 2 come from even and odd eigenfunctions of H ( 0 ) , respectively. Analogously, these solutions are even and odd, so that the parity is preserved after the transformation. However, both solutions are not square integrable, which is absolutely not a surprise, since H ( 1 , 2 ) in (71) cannot have bound states and only have continuous spectrum. This is a consequence of the presence of a node in the wave function of the odd ground state.
The use of the technique described in 4.1 allows for a direct construction of V ( 2 , 2 ) . This is given by:
V ( 2 , 2 ) = 2 x 2 log W { x , cos π x 2 a = π 2 2 a 2 π 2 x 2 2 cos π x a 2 π x sin π x 2 a + 2 cos π x 2 a 2 .
This potential is depicted in Figure 12. The energy eigenvalues for this potential are:
E n ( 2 ) = E n ( 1 ) E 1 ( 1 ) = E n ( 1 ) E 1 ( 0 ) = n ( n 1 ) π 2 a 2 , n = 2 , 3 , .
Observe that E 1 ( 2 ) is not degenerate, although all other energy levels are degenerate.

5. A Comment and Example on General Cases

So far, we have studied the exactly solvable cases here and in [20], all others admit numerical solutions only. In addition, we have considered situations in which the ground state wave function is real. What happens if this wave function were complex? Then, the first SUSY Hamiltonian H ( 1 , 2 ) is not even Hermitian. Let us see an example thereof.
We determine a self adjoint extension were the wave function of the ground state have not a defined parity, so that m 2 0 , and is not “time reversal invariant” in the sense given in [2], which means that m 3 0 . We may use either one of the two kinds of parameters, with equivalence given in (16), in order to fix a generic self adjoint extension without this two mentioned symmetries. We choose to fix the angles θ i , i = 0 , 1 , 2 , and ψ , so that:
θ 0 = θ 1 = θ 2 = ψ = π 10 { m 0 , m 1 , m 2 , m 3 , ψ } { 0.90 , 0.29 , 0.29 , 0.095 , 0.31 } .
The determination of the wave functions for the ground level and the first excited levels is possible. One sees that these wave functions are complex, so that we may depict their probability density, which we do in Figure 13. In Figure 14, we depict the logarithmic deviation of the square root of the energy levels with respect to the standard extension discussed in the textbooks.
The wave function of the ground state is:
ϕ 0 ( 0 ) ( x ) = ( 0.22 0.669 i ) cosh 0.111 x a ( 0.507 + 0.373 i ) sinh 0.111 x a .
Observe that this ground state wave function is complex, contrarily to the cases previously studied. The ground state energy is now,
E 0 ( 0 ) 0.012 1 a 2 ,
which is negative. In the present example, the only level having a negative energy is the ground level.
With these ingredients, we may obtain the Hamiltonian for the first SUSY partner, which is not Hermitian due to the fact that ϕ 0 ( 0 ) ( x ) is not real. The potential for H ( 1 , 2 ) is:
V ( 1 , 2 ) ( x ) = E 0 ( 0 ) 2 x 2 log ϕ 0 ( 0 ) ( x ) = 0.012 1 a 2 1 a 2 0.041 + 0.012 i ( 1 . cosh 0.111 x a + ( 0.276 0.845 i ) sinh 0.111 x a ) 2 .
In Figure 15 and Figure 16, we represent the real and the imaginary parts of the potential in (78) and the graphics for the probability density corresponding to the first six eigenfunctions of H ( 1 , 2 ) .
SUSY partners with complex potentials have previously been obtained. For instance, it typically appear when, instead a bound state, we use a decaying state or Gamow state function for a resonance to construct the partner instead the ground state of a self adjoint operator with discrete spectrum [23].

6. Discussion and Conclusions

Supersymmetric (SUSY) is a technique that may serve to construct Hamiltonians with discrete spectrum, similar to the spectrum of a given Hamiltonian. One of the most celebrated models in non-relativistic quantum mechanics is the one dimensional infinite square well, which is mathematically modelled by the operator D 2 : = d 2 / d x 2 on L 2 [ a , a ] , a > 0 being finite. However, there is not a unique self adjoint operator defined by D 2 on L 2 [ a , a ] , but instead a family of independent self adjoint determination of this operator depending on the values of four real parameters.
The objective of the present research is the investigation of the SUSY partners for these self adjoint determinations (also called extensions). For a large number of these extensions, under certain conditions such as a ground state with positive energy, no zeroes (also called nodes) and some symmetry conditions, all these SUSY partners have been classified in a previous paper [20]. For each of these determinations, we have obtained an infinite chain of partners, each one differing from the previous one by an eigenvalue. The shape of the potential and the eigenfunctions may differ from one partner to another.
The objective of the present article is the investigation of the SUSY partners of those extensions that do not satisfy the previous conditions. Some of them may have a ground state with negative energy that may have double degeneracy. Others have a ground state with zero energy that may also be doubly degenerate. Still keeping symmetry conditions for the wave function of the ground state, such as having a definite parity, we may find the exact form of the first SUSY partner Hamiltonian, and, hence, following a technique explained in [20], of all partners. When a double degeneracy on the ground state appears, we need to use a different technique that goes directly to the second partner. We show that extensions with a ground state with nodes in the wave function may lead to formal partners which do not satisfy the condition of having a similar discrete spectrum than of the original Hamiltonian and that may even have no eigenvalue.
Partners of extensions without symmetry conditions are difficult to classify. Here, the ground state wave function may even be complex and this gives rise to a complex partner potential, as to a non-Hermitian partner Hamiltonian. It would be of some interest to search for the possibility of finding PT-symmetric Hamiltonians among these non-Hermitian Hamiltonians obtained by this method.

Author Contributions

Conceptualization, M.G. and C.S.M.; methodology, M.G. and C.S.M.; software, M.G. and C.S.M.; validation, M.G. and C.S.M.; formal analysis, M.G. and C.S.M.; investigation, M.G. and C.S.M.; resources, M.G. and C.S.M.; data curation, M.G. and C.S.M.; writing—original draft preparation, M.G. and C.S.M.; writing—review and editing, M.G. and C.S.M.; visualization, M.G. and C.S.M.; supervision, M.G. and C.S.M.; project administration, M.G. and C.S.M.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by MCIN with funding from the European Union NextGenerationEU (PRTRC17.I1) and the Consejería de Educación through QCAYLE project, as well as the PID2020-113406GB-I0 project by MCIN of Spain.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Sara Cruz y Cruz and David J. Fernández C. for their support and interest in the present research.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. SUSY Partners

In this Appendix, we briefly review the foundations of the construction of the Supersymmetric (SUSY) partners of a given self adjoint operator with discrete spectrum. For a more complete account, see the bibliography on the subject listed at the end of the paper and particularly [10]. We use the notation introduced in [20].
To begin with, let us consider an initial one dimensional Hamiltonian H ( 0 ) = d 2 / d x 2 + V ( 0 ) ( x ) , where for simplicity we always assume that = 1 and m = 1 / 2 . Assume that H ( 0 ) has a non-degenerate discrete spectrum E i ( 0 ) i = 0 n with respective eigenfunctions ϕ i ( 0 ) i = 0 n , where n is either finite or infinite.
According to the standard formalism [10], the operator H ( 0 ) has two SUSY partners:
H ( 1 , 1 ) = d 2 d x 2 + V ( 1 , 1 ) ( x ) , H ( 1 , 2 ) = d 2 d x 2 + V ( 1 , 2 ) ( x ) ,
The superscripts of a given operator O ( p , q ) mean: p is the SUSY order, to be clarified later, and q = 1 , 2 . Here the indices ( 1 , 1 ) and ( 1 , 2 ) means: first order SUSY partner and partner one or partner two, respectively. Here, we need to define the partner potentials V ( 1 , i ) ( x ) , i = 1 , 2 , in (A1). To this end, let us take the ground state, ϕ 0 ( 0 ) ( x ) , of H ( 0 ) . Then,
V ( 1 , 1 ) ( x ) : = V 0 ( x ) + x 2 ϕ 0 ( 0 ) ( x ) ϕ 0 ( 0 ) ( x ) = V 0 ( x ) E 0 ( 0 ) ,
where x and x 2 stand for the first and second derivative with respect to the variable x. The second identity in (A2) is obvious, since ϕ 0 ( 0 ) ( x ) satisfies the Schrödinger equation for V 0 ( x ) and E 0 ( 0 ) . The second partner potential is
V ( 1 , 2 ) ( x ) = V ( 1 , 1 ) ( x ) 2 d 2 d x 2 log ϕ 0 ( 0 ) ( x ) .
It is interesting to consider the following two operators on L 2 ( R ) :
A ( 1 ) : = x x ϕ 0 ( 0 ) ( x ) ϕ 0 ( 0 ) ( x ) , ( A ( 1 ) ) : = x x ϕ 0 ( 0 ) ( x ) ϕ 0 ( 0 ) ( x ) ,
These operators are well defined provided that ϕ 0 ( 0 ) ( x ) does not have zeroes, that are usually called nodes. They intertwine the Hamiltonian (A1) in the following sense:
H ( 1 , 1 ) A ( 1 ) = ( A ( 1 ) ) H ( 1 , 2 ) , H ( 1 , 2 ) ( A ( 1 ) ) = A ( 1 ) H ( 1 , 1 ) .
The Hamiltonians H ( 1 , 1 ) and H ( 1 , 2 ) have the following properties:
  • Both have an infinite number of eigenvalues.
  • The eigenfunctions of H ( 1 , 1 ) , { ϕ i ( 1 , 1 ) ( x ) } i = 1 , coincide with the eigenfunctions of the original Hamiltonian H ( 0 ) , { ϕ i ( 0 ) ( x ) } i = 1 , although the spectrum of H ( 0 ) is of the form { ϵ i ( 0 ) : = E i ( 0 ) E 0 ( 0 ) } i = 0 , where E i ( 0 ) are the energy levels of H ( 0 ) and E 0 ( 0 ) its ground state. Thus, zero is an eigenvalue of H ( 1 , 1 ) .
  • The eigenfunctions of H ( 1 , 2 ) , { ϕ i ( 1 , 2 ) ( x ) } i = 1 are different from those of H ( 0 ) , while its spectrum coincides with the spectrum of H ( 1 , 1 ) with the ground level removed: { ϵ i ( 0 ) : = E i ( 0 ) E 0 ( 0 ) } i = 1 . Note that ϵ i ( 0 ) > 0 for i = 1 , 2 , .
  • Since the eigenfunctions of H ( 1 , 1 ) are known for being those of H ( 0 ) , we need a receipt so as to obtain { ϕ i ( 1 , 2 ) ( x ) } i = 1 , the eigenfunctions of H ( 1 , 2 ) . This is
    ϕ n ( 1 , 2 ) ( x ) = ϵ n 1 / 2 A ( 1 ) ϕ n ( 1 , 1 ) ( x ) , n = 1 , 2 , .
    In addition, we have the following inverse relation:
    ϕ n ( 1 , 1 ) ( x ) = ϵ n 1 / 2 ( A ( 1 ) ) ϕ n ( 1 , 2 ) ( x ) , n = 1 , 2 , .
    From (A4), (A6) and (A7), we understand the importance of the condition that ϕ 0 ( 0 ) ( x ) has no nodes.
This is the origin of a chain of SUSY partners of the given H ( 0 ) . Take H ( 1 , 2 ) , and consider a new pair of potentials, H ( 2 , 1 ) and H ( 2 , 2 ) , for which their potentials have the following form:
V ( 2 , 1 ) ( x ) : = V ( 1 , 2 ) ( x ) + x 2 ϕ 1 ( 1 , 2 ) ( x ) ϕ 1 ( 1 , 2 ) ( x ) ,
and
V ( 2 , 2 ) ( x ) : = V ( 2 , 1 ) ( x ) 2 d 2 d x 2 log ϕ 1 ( 1 , 2 ) ( x ) .
The operators (A4) are here replaced by
A ( 2 ) : = x x ϕ 1 ( 1 , 2 ) ( x ) ϕ 1 ( 1 , 2 ) ( x ) , ( A ( 2 ) ) : = x x ϕ 1 ( 1 , 2 ) ( x ) ϕ 1 ( 1 , 2 ) ( x ) .
We recall that x : = d / d x .
Thus, the Hamiltonian H ( 2 , 2 ) has a spectrum equal to ϵ n ϵ 1 = E n ( 0 ) E 0 ( 0 ) E 1 ( 0 ) + E 0 ( 0 ) = E n ( 0 ) E 1 ( 0 ) , n = 2 , 3 , .
Thus, we have started the construction of a chain of Hamiltonian SUSY partners of H ( 0 ) , where the first two Hamiltonians are H ( 1 , 2 ) and H ( 2 , 2 ) , first and second SUSY partners, and whose respective spectrum can be easily obtain form H ( 0 ) . We can proceed with this chain of SUSY partners of H ( 0 ) up to infinite. At the -th, we have another pair of Hamiltonians:
H ( , 1 ) : = d 2 d x 2 + V ( , 1 ) ( x ) , H ( , 2 ) : = d 2 d x 2 + V ( , 2 ) ( x ) ,
with
V ( , 1 ) ( x ) = V ( 0 ) ( x ) E 1 ( 0 ) + x 2 log W { ϕ 0 ( 0 ) , ϕ 1 ( 0 ) , , ϕ 2 ( 0 ) } ,
V ( , 2 ) ( x ) = V ( 0 ) ( x ) E 1 ( 0 ) x 2 log W { ϕ 0 ( 0 ) , ϕ 1 ( 0 ) , , ϕ 1 ( 0 ) } ,
where W denotes the Wronskian of the functions under brackets.
The eigenfunctions of the -th SUSY partner H ( , 2 ) are
ϕ n ( , 2 ) ( x ) = C n ( , 2 ) W { ϕ 0 ( 0 ) , ϕ 1 ( 0 ) , , ϕ 1 ( 0 ) , ϕ n ( 0 ) } W { ϕ 0 ( 0 ) , ϕ 1 ( 0 ) , , ϕ 1 ( 0 ) } ,
where C n ( , 2 ) stand for normalization constants. This general formula allows us to have the eigenfunctions of all SUSY partner Hamiltonians of the chain without having to solve the Schrödinger equation.
This is the procedure that we have to apply in Section 4.

Appendix B. Second Order SUSY Partners: Some Peculiarities for the Confluent Case

The purpose of this Appendix is the review on the construction of the second order partner directly from the original Hamiltonian. Although we essentially follow the presentation in [21] for the general case, we consider that we need to justify the reasons why we jump from the original Hamiltonian to the second order partner in the confluent case mentioned in Section 4.1. At the same time, we intend a pedagogical presentation for the benefit of the reader.
We begin with the original self adjoint Hamiltonian H ( 0 ) with its eigenfunctions { ϕ n ( 0 ) ( x ) } and respective eigenvalues { E n ( 0 ) } . The first partner potential H ( 1 , 1 ) has potential V ( 1 , 1 ) = V ( 0 ) E 0 ( 0 ) , same eigenfunctions as H ( 0 ) and eigenvalues those of H ( 0 ) sifted by E 0 ( 0 ) as described in the previous Appendix.
The usual procedure to obtain directly the second partner uses the construction of a second order differential operator B such that
B H ( 1 , 1 ) = H ( 2 , 2 ) B .
Compare (A15) with (A5). In the latter, the operator A ( 1 ) and its adjoint are first order operators as we remain in the first order. In (A15), we jump from first order to second order so that we need add one more order of differentiation to the intertwining operator. Thus, the operator B should have the form [21]:
B ( x ) : = d 2 d x + η ( x ) d d x + γ ( x ) ,
where η ( x ) and γ ( x ) are functions to be determined, by using (A16) into (A15). We obtain the following relations (recall that x : = d / d x and the prime also means derivation with respect to x):
x 2 V ( 1 , 1 ) ( x ) + η ( x ) x V ( 1 , 1 ) ( x ) + γ ( x ) V ( 1 , 1 ) ( x ) + γ ( x ) = γ ( x ) V ( 2 , 2 ) ( x ) ,
2 x V ( 1 , 1 ) ( x ) + η ( x ) V ( 1 , 1 ) ( x ) + 2 γ ( x ) + η ( x ) = η ( x ) V ( 2 , 2 ) ( x ) ,
V ( 2 , 2 ) ( x ) = V ( 1 , 1 ) ( x ) + 2 η ( x ) .
Equation (A18) permits to write γ ( x ) in terms of η ( x ) as
γ ( x ) = γ 0 V ( 1 , 1 ) ( x ) η ( x ) 2 + η ( x ) 2 2 ,
where γ 0 is an integration constant. Then, we need to obtain η ( x ) . Let us multiply (A18) by η ( x ) and integrate. The result is
4 η ( x ) 2 γ 0 V ( 1 , 1 ) ( x ) η ( x ) + 2 η ( x ) η ( x ) η ( x ) 2 + η ( x ) 4 + η 0 = 0 ,
where η 0 is a new integration constant. Equation (A21) is second order and non-linear and looks intractable at first sight.
Nevertheless, there is a way out. The operator B should annihilate the ground state, ϕ 0 ( 1 , 1 ) , of H ( 1 , 1 ) . This gives the following differential equation to be satisfied by η ( x ) :
B ( ϕ 0 ( 1 , 1 ) ) = 0 η ( x ) = 2 γ 0 + 2 η ( x ) x ϕ 0 ( 1 , 1 ) ( x ) ϕ 0 ( 1 , 1 ) ( x ) + η ( x ) 2 .
Then, if we use (A22) in (A21), we obtain a relation between both integration constants, which is η 0 = 4 γ 0 2 . If we choose the particular solution η 0 = 0 , we obtain the following solution [21,22]:
η ( x ) = d d x log ( w ( x ) ) , w ( x ) = w 0 x ω 2 ( y ) d y , ω ( x ) = ϕ 0 ( 1 , 1 ) ( x ) ,
where x f ( y ) d y denotes the primitive of the function f ( x ) .
Since B is a second order operator, it should exists a second function v ( x ) linearly independent with ϕ 0 ( 1 , 1 ) ( x ) and such that is annihilated by B: B ( v ( x ) ) = 0 . This implies that [21,22]:
η ( x ) = x log W ( v ( x ) , ϕ 0 ( 1 , 1 ) ( x ) ) .
If we compare (A24) to the first relation in (A23), we have that
W ( v ( x ) , ϕ 0 ( 1 , 1 ) ( x ) ) = w ( x ) .
This is a first order differential equation on v ( x ) . Choosing the value of the arbitrary constant equal to one, we have the following solution:
v ( x ) = ϕ 0 ( 1 , 1 ) ( x ) x w ( y ) ( ϕ 0 ( 1 , 1 ) ( y ) ) 2 d y ,
which is well defined in absence of nodes for ϕ 0 ( 1 , 1 ) ( x ) . The function v ( x ) is not an eigenfunction of H ( 1 , 1 ) . Instead,
H ( 1 , 1 ) ( v ( x ) ) = ϕ 0 ( 1 , 1 ) ( x ) [ H ( 1 , 1 ) ] 2 v ( x ) = H ( 1 , 1 ) ( ϕ 0 ( 1 , 1 ) ( x ) ) = 0 .
Thus, v ( x ) is an eigenfunction of the square of the Hamiltonian H ( 1 , 1 ) .
In Section 4.1, we have gone directly from the original Hamiltonian H ( 0 ) to the second partner H ( 2 , 2 ) without considering the first partner H ( 1 , 2 ) . This has a justification which is not given in [21,22]. To begin with this, we recall that H ( 1 , 1 ) and H ( 1 , 2 ) should be interconnected via a first order differential operator A 1 . Another first order differential operator, A 2 interconnects H ( 2 , 1 ) and H ( 2 , 2 ) . We have the following relations:
A 1 H ( 1 , 1 ) = H ( 1 , 2 ) A 1 , A 2 H ( 2 , 1 ) = H ( 2 , 2 ) A 2 , B H ( 1 , 1 ) = H ( 2 , 2 ) B ,
and H ( 1 , 2 ) = H ( 2 , 1 ) + E 1 ( 1 ) . From (A28), it comes that
B = A 2 A 1 .
In summary, operators A 1 , A 2 and B interconnect the mentioned Hamiltonian as follows:
H ( 1 , 1 ) A 1 H ( 1 , 2 ) H ( 2 , 1 ) A 2 H ( 2 , 2 ) , H ( 1 , 1 ) B H ( 2 , 2 ) .
So far, the treatment is quite general. Now, let us focus in the confluent case, where the ground state is doubly degenerate. With no degeneracy, B annihilates both ϕ 0 ( 1 , 1 ) and ϕ 1 ( 1 , 1 ) . With double degeneracy for the ground state, B annihilates ϕ 0 ( 1 , 1 ) and v ( x ) . The former is an eigenvalue of H ( 1 , 1 ) , although this is not the case for v ( x ) . It is not difficult to finds the pair of first order differential operators A 1 and A 2 verifying (A28) and (A29). In fact, two solutions are possible for this pair of operators. The first solution is
A 1 = x x ϕ 0 ( 1 , 1 ) ( x ) ϕ 0 ( 1 , 1 ) ( x ) , A 2 = x x W ( ϕ 0 ( 1 , 1 ) , v ) ( x ) ϕ 0 ( 1 , 1 ) ( x ) ,
and the second solution is
A 1 = x v ( x ) v ( x ) , x x W ( v , ϕ 0 ( 1 , 1 ) ) ( x ) v ( x ) .
The first solution does not add anything new. For this solution, A 1 ( ϕ 0 ( 1 , 1 ) ) = 0 and then, H ( 1 , 2 ) is the first order partner without the ground state of H ( 0 ) . In addition, the state A 1 ( v ) is not an eigenfunction of H ( 1 , 2 ) , since v ( x ) is not an eigenfunction of H ( 1 , 1 ) .
Take the second solution, which is the solution that we should take into account for the confluent case. Here A 1 ( v ) = 0 , although A 1 does not annihilate an eigenstate as in the previous case, since v ( x ) is not ein eigenfunction of H ( 1 , 1 ) . According to the general theory, the eigenfunctions of H ( 1 , 2 ) should have the form:
ϕ n ( 1 , 2 ) = A 1 ( ϕ n ( 1 , 1 ) ) = A 1 ( ϕ n ( 0 ) ) = W ( { v ( x ) , ϕ n ( 0 ) ( x ) } ) v ( x ) .
Let us see whether the form of these functions is compatible with the properties they should have as eigenfunctions of H ( 1 , 2 ) . In fact,
H ( 1 , 2 ) ϕ n ( 1 , 2 ) = H ( 1 , 2 ) A 1 ϕ n ( 1 , 1 ) = A 1 H ( 1 , 1 ) ϕ n ( 1 , 1 ) = E n ( 0 ) A 1 ( ϕ n ( 1 , 1 ) ) = E n ( 1 ) A 1 ( ϕ n ( 0 ) ) .
According to (A32) and (A33), the eigenfunctions ϕ n ( 1 , 1 ) ( x ) , n = 1 , 2 , , should depend on v ( x ) , which is not possible, since
V ( 1 , 2 ) = V ( 1 , 1 ) 1 w 0 ϕ 0 ( 1 , 1 ) ( x ) 2 d x ϕ 0 ( 1 , 1 ) ( x ) 2 d x ,
which does not depend on v ( x ) . Thus, this second solution is spurious. Therefore, the first partner H ( 1 , 2 ) is ill defined. Thus is the motivation that we should go directly to the second partner when the bound state shows degeneracy.

Appendix C. Self Adjoint Extensions of Symmetric Operators

The terminology self adjoint determinations or self adjoint extensions comes from the following: The domain, D A , of a (linear) operator A is the subspace of the Hilbert space, H , of all vectors on which A acts. Bounded operators have the whole Hilbert space as domain. However, this is not true for the unbounded operators such as the differential operator D 2 . Most of operators used in quantum mechanics, position, momentum, angular momentum, Hamiltonians, are unbounded operators. An unbounded operator has two parts, the prescription on how it acts on the vector of its domain, for instance ( D 2 ) f ( x ) = f ( x ) , and the domain itself. Two operators with analogous action on their vectors, for instance take the second derivative, having different domains are different, so that an unbounded operator is a pair ( A , D A ) . Two operators ( A , D A ) and ( B , D B ) are equal if and only if their domains are equal, D A D B , and being given an arbitrary f D A D B , A f = B f .
Let ( A , D A ) be an unbounded operator such that the domain D A is dense in the Hilbert space H . This means that for any g H and for any neighbourhood V of g, there exists at least one f D A such that f V . The adjoint A of A is an operator with domain
D A = { g H such that h H ; h | f = g | A f f D A } .
In this case,
A g = h , g D A .
We call A the adjoint operator of A. It is well defined if and only if D A is dense in H .
An operator ( B , D B ) extends another operator ( A , D A ) , A B , if and only if D A D B and for any g D A , A g = B g .
An operator ( A , D A ) is symmetric (also called Hermitian) if for all f , g D A , one has that
f | A g = A f | g .
After (A38), we note that for all symmetric operator A, A A . A symmetric operator A is self adjoint if A = A . This means in particular that D A D A .
Let I be the identity operator, i.e., I f = f , f H and N ( A ) the nucleus (kernel) of the operator A. The deficiency indices of an unbounded symmetric operator A are n ± : = dim N ( A ± i I ) , where dim V means the dimension of the subspace V.
A well know theorem [3] establishes that A has self adjoint extensions if and only if n + = n 0 . If n + = n = 0 , then A admits one and only one self adjoint extension, that coincides with A if this is self adjoint.
In the case of the operator D 2 on L 2 [ a , a ] , we take as its domain, D , the functions f ( x ) L 2 [ a , a ] satisfying the following properties:
(i). The function f ( x ) is absolutely continuous on the interval [ a , a ] , having a square integrable second derivative, on the same interval, defined almost everywhere.
(ii). For any f ( x ) D , one has that
f ( a ) = f ( a ) = f ( a ) = f ( a ) = 0 .
This operator is symmetric with deficiency indices n ± = 2 [1]. It has an infinite number of self adjoint extensions (also called determinations along the present paper) depending on four real parameters as described in the Introduction, see (4) and comment thereafter.

References

  1. Kurasov, P. Distribution Theory for Discontinuous Test Functions and Differential Operators with Generalized Coefficients. J. Math. Ann. Appl. 1996, 201, 297–323. [Google Scholar] [CrossRef] [Green Version]
  2. Bonneau, G.; Faraut, J.; Valent, G. Self adjoint extensions of operators and the teaching of quantum mechanics. Am. J. Phys. 2001, 69, 322–331. [Google Scholar] [CrossRef] [Green Version]
  3. Reed, M.; Simon, B. Fourier Analysis. Self Adjointness; Academic Press: New York, NY, USA, 1975. [Google Scholar]
  4. Naimark, M.A. Linear Differential Operators; Dover: New York, NY, USA, 2014. [Google Scholar]
  5. Fernández, D.J. Supersymmetric quantum mechanics. AIP Conf. Proc. 2010, 1287, 3–36. [Google Scholar]
  6. Infeld, L.; Hull, T.E. The factorization method. Rev. Mod. Phys. 1951, 23, 21–68. [Google Scholar] [CrossRef]
  7. Lahiri, A.; Roy, P.K.; Bagchi, B. Supersymmetry in quantum mechanics. Int. J. Mod. Phys. A 1990, 5, 1383–1456. [Google Scholar] [CrossRef] [Green Version]
  8. Roy, B.; Roy, P.; Roychoudhury, R. On the solution of quantum eigenvalue problems. A supersymmetric point of view. Fortschr. Phys. 1991, 39, 211–258. [Google Scholar] [CrossRef]
  9. Bagchi, B. Supersymmetry in Quantum and Classical Mechanics; Chapman and Hall: Boca Raton, FL, USA, 2001. [Google Scholar]
  10. Cooper, F.; Khare, A.; Sukhatme, U. Supersymmetry and Quantum Mechanics. Phys. Rep. 1995, 251, 267–385. [Google Scholar] [CrossRef] [Green Version]
  11. Díaz, J.I.; Negro, J.; Nieto, L.M.; Rosas-Ortiz, O. The supersymmetric modified Pöschl-Teller and delta well potential. J. Phys. A Math. Gen. 1999, 32, 8447–8460. [Google Scholar] [CrossRef]
  12. Mielnik, B.; Nieto, L.M.; Rosas-Ortiz, O. The finite difference algorithm for higher order supersymmetry. Phys. Lett. A 2000, 269, 70–78. [Google Scholar] [CrossRef] [Green Version]
  13. Fernández, D.J.; Negro, N.; Nieto, L.M. Regularized scarf potentials: Energy band structure and supersymmetry. J. Phys. A Math. Gen. 2004, 37, 6987–7001. [Google Scholar]
  14. Ioffe, M.V.; Negro, J.; Nieto, L.M.; Nishiniandze, D. New two dimensional integrable quantum models from SUSY intertwining. J. Phys. A Math. Gen. 2006, 39, 9297–9308. [Google Scholar] [CrossRef]
  15. Correa, F.; Jakubsk, V.; Nieto, L.M.; Plyushchay, M.S. Self-isospectrality, special supersymmetry, and their effect in band structure. Phys. Rev. Lett. 2008, 101, 030403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Fernández, D.J.; Gadella, M.; Nieto, L.M. Supersymmetry transformations for delta potentials. SIGMA 2001, 7, 029. [Google Scholar] [CrossRef]
  17. Fernández, D.J.; Morales-Salgado, V.S. Supersymmetric partners of the harmonic oscillator with an infinite potential barrier. J. Phys. A Math. Theor. 2014, 47, 035304. [Google Scholar] [CrossRef] [Green Version]
  18. Cruz y Cruz, S.; Santiago-Cruz, C. Position dependent mass Scarf Hamiltonians generated via the Riccati equation. Math. Meth. Appl. Sci. 2019, 42, 4909–4924. [Google Scholar] [CrossRef] [Green Version]
  19. Cruz y Cruz, S.; Romero-Osnaya, A.; Rosas-Ortiz, O. Balanced Gain-and-Loss Optical Waveguides: Exact Solutions for Guided Modes in Susy-QM. Symmetry 2021, 13, 1583. [Google Scholar] [CrossRef]
  20. Gadella, M.; Hernández-Muñoz, J.; Nieto, L.M.; San Millán, C. Supersymmetric partners of the one-dimensional infinite square well hamiltonian. Symmetry 2021, 13, 350. [Google Scholar] [CrossRef]
  21. Fernández, D.J.; Salinas-Hernández, E. Wronskian formula for confluent second order supersymmetric quantum mechanics. Phys. Lett. A 2005, 338, 13–18. [Google Scholar] [CrossRef] [Green Version]
  22. Fernández, D.J.; Salinas-Hernández, E. The confluent algorithm in second-order supersymmetric quantum mechanics. J. Phys. A Math. Gen. 2003, 36, 2537–2543. [Google Scholar]
  23. Çevik, D.; Gadella, M.; Kuru, Ş.; Negro, J. Resonances and antibound states for the Pöschl-Teller potential: Ladder operators and SUSY partners. Phys. Lett. A 2016, 380, 1600–1609. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Regions in the parameter space ( θ 0 , ψ ) for which we have one or two negative eigenvalue with negative energy, as explained in the text. White zones correspond to solutions with positive energy eigenvalues only.
Figure 1. Regions in the parameter space ( θ 0 , ψ ) for which we have one or two negative eigenvalue with negative energy, as explained in the text. White zones correspond to solutions with positive energy eigenvalues only.
Symmetry 14 01314 g001
Figure 2. The two last surfaces on Table 1, corresponding to even (orange) and odd (blue) wave function. The parameter space is here ( θ 0 , ψ , E ) with E = r 2 / ( 2 a ) .
Figure 2. The two last surfaces on Table 1, corresponding to even (orange) and odd (blue) wave function. The parameter space is here ( θ 0 , ψ , E ) with E = r 2 / ( 2 a ) .
Symmetry 14 01314 g002
Figure 3. Values of ( θ 0 , θ 1 , ψ ) giving extensions with ground state of negative energy, the black net represents where a state with null energy is located. Note that these regions are three dimensional.
Figure 3. Values of ( θ 0 , θ 1 , ψ ) giving extensions with ground state of negative energy, the black net represents where a state with null energy is located. Note that these regions are three dimensional.
Symmetry 14 01314 g003
Figure 4. Regions in which the first excited state with positive energy is even, orange or odd, in blue.
Figure 4. Regions in which the first excited state with positive energy is even, orange or odd, in blue.
Symmetry 14 01314 g004
Figure 5. Surfaces corresponding the positive eigenvalues. Orange and blue colors represent even or odd states, respectively.
Figure 5. Surfaces corresponding the positive eigenvalues. Orange and blue colors represent even or odd states, respectively.
Symmetry 14 01314 g005
Figure 6. Eigenfunctions of the self adjoint extension of H ( 0 ) corresponding to the values of the parameters θ 0 , ψ = π 2 + arccot ( 2 ) , π 2 arccot ( 2 ) . Here, the gound state has zero energy and is doubly degenerate. We represent the two zero energy wave functions in black. The black continuous line represents the even wave function and the dashed line the odd wave function. Then, we depict the first four positive energy levels.
Figure 6. Eigenfunctions of the self adjoint extension of H ( 0 ) corresponding to the values of the parameters θ 0 , ψ = π 2 + arccot ( 2 ) , π 2 arccot ( 2 ) . Here, the gound state has zero energy and is doubly degenerate. We represent the two zero energy wave functions in black. The black continuous line represents the even wave function and the dashed line the odd wave function. Then, we depict the first four positive energy levels.
Symmetry 14 01314 g006
Figure 7. Deviations of the energy levels with odd parity corresponding to the extension labelled with θ 0 , ψ = π 2 + arccot ( 2 ) , π 2 arccot ( 2 ) , with respect to the odd parity energy levels corresponding to the extension with θ 0 , ψ = π 2 , π 2 .
Figure 7. Deviations of the energy levels with odd parity corresponding to the extension labelled with θ 0 , ψ = π 2 + arccot ( 2 ) , π 2 arccot ( 2 ) , with respect to the odd parity energy levels corresponding to the extension with θ 0 , ψ = π 2 , π 2 .
Symmetry 14 01314 g007
Figure 8. Eigenfunctions of H ( 2 , 2 ) when the ground state of H ( 0 ) is doubly degenerate with zero energy. We use ϕ 1 ( 0 ) as the point of departure for the confluent SUSY. The parameter w 0 has been chosen so that the singularity in the potential lies at x = 1.2 a .
Figure 8. Eigenfunctions of H ( 2 , 2 ) when the ground state of H ( 0 ) is doubly degenerate with zero energy. We use ϕ 1 ( 0 ) as the point of departure for the confluent SUSY. The parameter w 0 has been chosen so that the singularity in the potential lies at x = 1.2 a .
Symmetry 14 01314 g008
Figure 9. Same as in Figure 8, starting with ϕ 0 ( 0 ) .
Figure 9. Same as in Figure 8, starting with ϕ 0 ( 0 ) .
Symmetry 14 01314 g009
Figure 10. Potential for the confluent first order SUSY corresponding to the self adjoint extension of H ( 0 ) determined by the parameter values { θ 0 , ψ } = { 3 π / 2 , π / 2 } .
Figure 10. Potential for the confluent first order SUSY corresponding to the self adjoint extension of H ( 0 ) determined by the parameter values { θ 0 , ψ } = { 3 π / 2 , π / 2 } .
Symmetry 14 01314 g010
Figure 11. Potential V ( 1 , 2 ) for the self adjoint extension of H ( 0 ) fixed by { θ 0 , ψ } = { π / 2 , π / 2 } . This potential has continuous spectrum only. This anomaly is due to the presence a node at the origin in the odd wave function of the ground state.
Figure 11. Potential V ( 1 , 2 ) for the self adjoint extension of H ( 0 ) fixed by { θ 0 , ψ } = { π / 2 , π / 2 } . This potential has continuous spectrum only. This anomaly is due to the presence a node at the origin in the odd wave function of the ground state.
Symmetry 14 01314 g011
Figure 12. Potential V ( 2 , 2 ) corresponding to the self adjoint extension of H ( 0 ) with { θ 0 , ψ } = { π / 2 , π / 2 } .
Figure 12. Potential V ( 2 , 2 ) corresponding to the self adjoint extension of H ( 0 ) with { θ 0 , ψ } = { π / 2 , π / 2 } .
Symmetry 14 01314 g012
Figure 13. Probability densities of the first six bound states corresponding to the self adjoint extension under our study.
Figure 13. Probability densities of the first six bound states corresponding to the self adjoint extension under our study.
Symmetry 14 01314 g013
Figure 14. Logarithmic deviation of the square root of the energy levels with respect to the standard textbook solution. Blue and red colours represent odd and even labels, respectively.
Figure 14. Logarithmic deviation of the square root of the energy levels with respect to the standard textbook solution. Blue and red colours represent odd and even labels, respectively.
Symmetry 14 01314 g014
Figure 15. Real and imaginary parts of V ( 1 , 2 ) .
Figure 15. Real and imaginary parts of V ( 1 , 2 ) .
Symmetry 14 01314 g015
Figure 16. Probability densities of the first five bound states corresponding to the first SUSY partner constructed with the negative energy ground state (76).
Figure 16. Probability densities of the first five bound states corresponding to the first SUSY partner constructed with the negative energy ground state (76).
Symmetry 14 01314 g016
Table 1. Transcendental equations and ground state functions for the studied extensions. Note that we two solutions for each case, depending whether the wave function is even or odd. Under the label “kind”, we give the sign of the energy of the ground state and the parity of this wave function.
Table 1. Transcendental equations and ground state functions for the studied extensions. Note that we two solutions for each case, depending whether the wave function is even or odd. Under the label “kind”, we give the sign of the energy of the ground state and the parity of this wave function.
KindTranscendental EquationGround State Eigenfunction
Positive and even s sin s 2 sin θ 0 + ψ 2 + cos s 2 cos θ 0 + ψ 2 = 0 s a ( s + sin s ) cos s x 2 a ,
Positive and odd s cos s 2 sin θ 0 ψ 2 + sin s 2 cos θ 0 ψ 2 = 0 s a ( s sin s ) sin s x 2 a ,
Zero and even m 0 + 2 m 1 + cos ψ 2 sin ψ = 0 1 2 a ,
Zero and odd m 0 + 2 m 1 + cos ψ 2 sin ψ = 0 3 2 a 3 x ,
Negative and even r sinh r 2 sin θ 0 + ψ 2 cosh r 2 cos θ 0 + ψ 2 = 0 r a ( sinh r + r ) cosh r x 2 a ,
Negative and odd r cosh r 2 sin θ 0 ψ 2 + sinh r 2 cos θ 0 ψ 2 = 0 r a ( sinh r r ) sinh r x 2 a .
Table 2. In the first column, we list the regions in accordance to Figure 1. In the second column, we give the sign of the energy of the first excited state and the parity of its wave function, provided that the ground state have zero or negative energy. Same for the first excited state in the third. Finally, the last column gives the signs of the forth first labels. All other levels have positive energy. The red line shows that the ground state has negative energy and is double degenerate. The red dot in Figure 1 shows a double degenerate ground state with zero energy, as indicated in the last row.
Table 2. In the first column, we list the regions in accordance to Figure 1. In the second column, we give the sign of the energy of the first excited state and the parity of its wave function, provided that the ground state have zero or negative energy. Same for the first excited state in the third. Finally, the last column gives the signs of the forth first labels. All other levels have positive energy. The red line shows that the ground state has negative energy and is double degenerate. The red dot in Figure 1 shows a double degenerate ground state with zero energy, as indicated in the last row.
Region LabelGround StateFirst Excited StateSign of the Energy for the First Levels
(a)Negative and evenPositive(−,+,+,+,…)
(b)Negative and oddPositive(−,+,+,+,…)
(c)Negative and evenNegative and odd(−,−,+,+,…)
(d)Negative and oddNegative and even(−,−,+,+,…)
(e)PositivePositive(+,+,+,+,…)
(b) ∩ (d)Negative and oddZero and even(−,0,+,+,…)
(a) ∩ (c)Negative and evenZero and odd(-,0,+,+,…)
(e) ∩ (b)Zero oddPositive(0,+,+,+,…)
(e) ∩ (a)Zero evenPositive(0,+,+,+,…)
Red lineNegative and doubly degenerateSame(−,−,+,+,…)
Red dotZero and doubly degenerateSame(0,0,+,+,…)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gadella, M.; San Millán, C. Supersymmetric Partners of the One-Dimensional Infinite Square Well Hamiltonian: Special Cases. Symmetry 2022, 14, 1314. https://doi.org/10.3390/sym14071314

AMA Style

Gadella M, San Millán C. Supersymmetric Partners of the One-Dimensional Infinite Square Well Hamiltonian: Special Cases. Symmetry. 2022; 14(7):1314. https://doi.org/10.3390/sym14071314

Chicago/Turabian Style

Gadella, Manuel, and Carlos San Millán. 2022. "Supersymmetric Partners of the One-Dimensional Infinite Square Well Hamiltonian: Special Cases" Symmetry 14, no. 7: 1314. https://doi.org/10.3390/sym14071314

APA Style

Gadella, M., & San Millán, C. (2022). Supersymmetric Partners of the One-Dimensional Infinite Square Well Hamiltonian: Special Cases. Symmetry, 14(7), 1314. https://doi.org/10.3390/sym14071314

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