CC BY
76
Level crossings of eigenvalues of the Schrodinger Hamiltonian of the isotropic harmonic oscillator perturbed by a central point interaction in different dimensions
S. Fassari1'2'3, M. Gadella1, M.L. Glasser1'4, L.M. Nieto1, F. Rinaldi2'3
departamento de Física Teórica, Atómica y (Optica, U. Valladolid, 47011 Valladolid, Spain 2CERFIM, PO Box 1132, Via F. Rusca 1, CH-6601 Locarno, Switzerland
3Dipartimento di Fisica Nucleare, Subnucleare e delle Radiazioni, Univ. degli Studi Guglielmo Marconi, Via Plinio 44, I-00193 Rome, Italy 4Department of Physics, Clarkson University, Potsdam, NY, 13699 USA [email protected], [email protected], [email protected], [email protected], [email protected]
PACS 02.30.Gp, 02.30.Lt, 02.30.Sa, 02.30.Tb, 03.65.Db, 03.65.Ge, 68.65.Hb DOI 10.17586/2220-8054-2018-9-2-179-186
In this brief presentation, some striking differences between level crossings of eigenvalues in one dimension (harmonic or conic oscillator with a central nonlocal ¿'-interaction) or three dimensions (isotropic harmonic oscillator with a three-dimensional delta located at the origin) and those occurring in the two-dimensional analogue of these models will be highlighted.
Keywords: harmonic oscillator, point interactions, renormalisation, eigenvalues, level crossing.
Received: 8 January 2018 Revised: 24 January 2018
1. Introduction
The aim of this note is to compare the discrete spectrum of the isotropic harmonic oscillator perturbed by a point interaction centered at the origin, the bottom of a parabolic well, across various dimensions. The physical motivation behind this detailed study is the fact that this Hamiltonian, at least in the three-dimensional case, has been used over the last twenty years to model two-dimensional quantum wires and three-dimensional quantum dots.
It is important to point out and is well known that, while in two or three dimensions there is only one type of point interaction, namely the 2/3D delta (requiring a renormalization procedure if it is to be a perturbation of the negative Laplacian or any other free Hamiltonian), in one-dimensional Quantum Mechanics we have four different types of point interactions, out of which, the only one requiring such renormalization is the nonlocal S'—interaction, defined by its quadratic form (•, |S'} (S'|, •) (see [1-3]). A self-adjoint Hamiltonian with the nonlocal S'—interaction can also be defined as the norm resolvent limit of a Hamiltonian involving three Dirac distributions (see [4-6]). We remind the reader that this interaction, being an element of H-2\H-1, may not be combined with the Dirac distribution, obviously an element of H-1 (the reader is referred to [2] for a thorough description of the Sobolev spaces labelled by a negative index).
While the spectral features of the three-dimensional isotropic oscillator perturbed by a single point interaction centered at the origin have been extensively investigated (see [7-10]), less attention has been paid to those of its two-dimensional analog. The remarkable spectral feature of the three-dimensional case is the existence of level crossings: each eigenvalue created by the point perturbation, clearly affecting only the s-states of the model, regarded as a function of either the parameter a labelling the possible self-adjoint extensions or its reciprocal ¡ (acting as the coupling constant), crosses the next lower unperturbed energy level (pertaining to an antisymmetric state with an odd value of the total angular momentum). Even more remarkable, all the level crossings take place for the same coupling parameter value.
As was shown in [11] (respectively [12,13]), level crossings also occur in the entirely discrete spectrum of the one-dimensional harmonic (respectively conic) oscillator perturbed by a nonlocal S'-interaction. This one-dimensional model shares the property of having all the crossings situated at the same point when either the extension parameter or the coupling parameter is chosen. However, the nature of such crossings is different from the aforementioned in three dimensions, as the eigenenergy of each antisymmetric bound state, regarded as a function of either the extension parameter a, labeling the self-adjoint extensions, or the coupling ¡ , its reciprocal, crosses the next lower unperturbed energy level (pertaining to a symmetric state). With reference to the terminology
introduced in [14], one could describe the one-dimensional model as Fermionic and the three-dimensional model as Bosonic.
Therefore, we believe it may be worth investigating the spectral properties of the two-dimensional analog in depth given the increasing importance of this particular dimension in modern Nanophysics. As expected, level crossings also occur in the two-dimensional case. However, it will be seen that the location of these crossings is no longer at a unique parameter value. Furthermore, while in odd dimensions the interaction must be attractive in order to produce the crossings, a weaker repulsive interaction will suffice in the two-dimensional model.
2. Level crossings in the discrete spectrum of the one-dimensional harmonic or conic oscillator perturbed by a central nonlocal ¿'-interaction
In this section, we shall review the findings of [11] (respectively [12,13]) showing the occurrence of level crossings in the discrete spectrum of the harmonic (resp. conic) oscillator perturbed by a nonlocal S'—interaction centred at the minimum of the parabolic (resp. V-shaped) potential.
It is well known that the Schrodinger Hamiltonian perturbed by a Dirac distribution need not be defined by means of renormalization, since such a perturbation is infinitesimally small in the sense of quadratic forms (see [1,2]). Therefore, the Dirac distribution can be regarded as an element of H-1.
As anticipated earlier, the nonlocal S'-perturbation of the one-dimensional negative Laplacian, being an element of H-2\H-1, requires instead either the self-adjoint extension method or coupling constant renormalisation to produce a self-adjoint Hamiltonian (see [1,2]).
The spectral effects of such a singular perturbation on the harmonic oscillator were studied in [11] starting
from the explicit expression for the resolvent of the self-adjoint operator Hp acting like H0 = 1
the functions satisfying the condition 0(0+) — 0(0-) = —fi0'(0):
2
dx2
on
[H Ej- -m Ej- , 2 l*(E))ME )l (1)
[H3 - Ej - [H° - Ej + 1 1 2r(3/4-E/2) , (1)
3 3o + r(1/4-E/2)
r(1/4)
with ft - 2^3/4)) ^ 1.47934 and:
E) ^ (2n +1)1/V2n(0) , () (2)
*(x> E) — X, (2n + 3/2 - E) ^2"+l(x)' (2)
0n(x) being the n-th normalized eigenfunction of the harmonic oscillator. The latter function, having a jump
discontinuity satisfying the condition 0(0+) - 0(O-) — -,00'(O), is square summable since the sequence 02n(0) behaves like n-1/4 as n ^ œ. Hence, the equation determining the bound state energies of H3 is:
1 1 _ 2r(3/4 - E/2) (3)
fi0 fi r(l/4 — E/2) •
Clearly, given that the perturbation does not affect the symmetric subspace, the even eigenvalues E2n = 2n + 1/2 of the unperturbed harmonic oscillator are still in the discrete spectrum for any value of fi. Each new odd energy level created by the perturbation is a decreasing function of the coupling parameter fi (as can be easily seen by means of implicit differentiation), namely E2n+1(fi). For fi = 0, the operator trivially coincides with H0, so the odd eigenvalues are also those of the harmonic oscillator E2n+1 = 2n + 3/2. Something remarkable occurs for fi = fi0: as the left hand side of (3) vanishes, the only values of the energy that can satisfy (3) are those that make the Gamma function in the denominator diverge, which implies that 1/4 — E/2 = —n, or equivalently E2n+1(fi0) = 2n +1/2. Therefore, the discrete spectrum of Hp0 consists only of the doubly degenerate eigenvalues E2n(fi0) = E2„+1(fi0) = 2n +1/2. As E2„+1(fi) < E2„(fi) = 2n + 1/2 for any fi > fi0, we have infinitely many level crossings to the right of the critical value fi0 of the coupling parameter. The plot shown in Fig. 1 depicts the three lowest level crossings.
It is worth pointing out from a physical point of view that, since the symmetry of the ground state eigenfunc-tion changes from symmetric to antisymmetric for any fi > fi0, the model exhibits a quantum phase transition. Furthermore, the effect produced by the singular S'—perturbation is far stronger than the Zeldovich effect produced by the S—interaction due to the fact that the horizontal asymptotes of the curves are situated below the horizontal lines pertaining to the unperturbed symmetric eigenvalues (see [15-17] for a detailed description of the spectrum of the 1D harmonic oscillator perturbed by Dirac distributions).
-2
-4
-4
-2
Fig. 1. The curves of the four lowest eigenenergies of Hp as functions of the coupling param-
r(1/4)
eter fi, and the level crossings occurring at fi = fi0 = 2p(3/4) ~ 1.47934. The horizontal
lines are the even energy levels of the unperturbed harmonic oscillator: E0 = 1/2 (yellow), E2 = 5/2 (green), E4 = 9/2 (red)
The same spectral pattern has recently been encountered in the discrete spectrum of the self-adjoint Hamiltonian
1 T d2 1
of the conic oscillator HC acting like HC = - — + |x| on functions satisfying the condition —(0+) ——(0 ) = p 2 L dx2 J
—ß—'(0), whose resolvent is explicitly given by (see [12,13] for the procedure requiring renormalization needed to define HC):
H — E ] 1 = [H0 — E]-1 +
|*C(E ))(*C(E )|
with:
Ai(0)Ai'(-2E)-Ai'(0)Ai(-2E) Ai(0)Ai(-2E)
—2n(x)
21/2£^ (E2n — E)'
(4)
(5)
where ^2n(x) is the 2n-th normalized eigenfunction (pertaining to an antisymmetric bound state) of H0 and E2n the corresponding eigenvalue, fully investigated in [18]. The function E) is also a square summable function since the sequence (E2n - E)-1 behaves like n-2/3 as n ^ to and has a jump discontinuity at the origin satisfying the condition ^(0+) - ^(0-) = -fi^'(0).
As a consequence, the bound state equation reads:
Ai(0)Ai'(-2E) - Ai'(0)Ai(-2E)
or equivalently:
1
ß
ß'-
Ai(0)Ai(—2E) Ai(0)Ai(—2E)
(6)
(7)
Ai(0)Ai'(-2E) - Ai'(0)Ai(-2E)"
Given that Ai'(-2E2n-1) = 0 for any n > 1, it follows that, if the energy is exactly equal to any odd eigenvalue (pertaining to a symmetric bound state), the right hand side of (7) is equal to:
ßc = Ai(°)
ß0 = —
Ai'(0)
(8)
which implies that such eigenvalues cause the denominator of the rank one operator on the right hand side of (4)
for the special value of the coupling parameter fi = fiC = - ^ffl ~ 1.37172 to vanish. This leads to their double
Ai (0)
degeneracy as eigenvalues of the Hamiltonian H^ and the ensures level crossings at that particular value of the coupling.
E
4
2
0
ß
0
2
4
The blue curves in Fig. 2 depict E0(ft),E2(ft) and E4(ft), the three lowest lying energy levels pertaining to antisymmetric bound states generated by the singular perturbation as functions of the coupling parameter. The yellow horizontal line represents E' and the green one E3, the energies of the two lowest symmetric bound states. Quite obviously, the spacing between two consecutive energy levels of HC, the Hamiltonian of the unperturbed conic oscillator, is not constant in contrast with what happens for the harmonic oscillator. Nevertheless, the spectral features are qualitatively identical to those observed in the latter case: all the level crossings take place for the same value of the coupling parameter which is located to the right of the origin. Identically to what has already been seen in the case of the harmonic oscillator, the effect produced by this perturbation is far stronger than the one named after Zeldovich occurring for an attractive ¿-interaction.
E
-4 0 1.37172 4
2.0
1.5
1.0
0.5
Fig. 2. The curves of the three lowest eigenenergies as functions of the coupling parameter ft, and the level crossings occurring at ft = ftC _ Ai(0)
~ 1.37172. The horizontal lines are the
Ai' (0)
energy levels of the unperturbed conic oscillator: E' (yellow) and E3 (green)
|8
-4
0 1.37172
4
3. Level crossings in the discrete spectrum of the three-dimensional isotropic harmonic oscillator perturbed by a central ¿-interaction
It is well known that the three-dimensional negative Laplacian perturbed by the Dirac distribution cannot be defined by means of the KLMN theorem since such a perturbation is no longer infinitesimally small in the sense of quadratic forms (see [1,2]). The same obviously holds in the presence of an harmonic confining potential.
The self-adjoint operator Hp, representing the Hamiltonian of he three-dimensional isotropic harmonic oscillator perturbed by a central ¿-interaction, defined either by means of the self-adjoint extension method or renormalization, has been investigated in depth in [7-10]. The entirely discrete spectrum of the operator, in particular the lowest lying eigenvalues, has also been studied in detail in those papers. The operator is rigorously defined by means of the explicit formula for its resolvent, namely:
H - E]-' = [Ho - E]-' + 2|*<E»WE^
l _ X + r(3/4-E/2) ' P Po + nr(l/4-E/2)
*T(1/4)
with ft0 = „A ,./ ~ 9.29495 and:
r(3/4)
TO TO TO
if,r T?\ ST ST ST ^2ni (o)02n2 (O)02na(0) , , s, , s, , s
ni=0 n2=0 n3 = 0
whose norm, for any E < 0, is bounded by:
2n' + 2n2 + 2n3 + 3 - E
=0 (2n + l)2/3
,n=0
3/2
< TO,
(9)
(10)
(11)
since the sequence ^2n(0) behaves like n 1/4 as n ^ to.
3
As is well known, the eigenenergy 2ni + 2n2 + 2n3 + - in the discrete spectrum of the three-dimensional
isotropic harmonic oscillator has degeneracy ^(2n1 + 2n2 + 2n3 + l)(2n1 + 2n2 + 2n3 + 2). Out of all the eigenstates sharing this energy, only the one having radial symmetry will be affected by the 3D point perturbation and give rise to a new simple eigenvalue, thus lowering the original degeneracy of the level to - (2n1 + 2n2 + 2n3 + 1)(2n1 + 2n2 + 2n3 + 2) - l. This new simple eigenvalue is determined as a function of the coupling parameter by the equation:
1 1 r(3/4 - E/2)
£c £ nr(1/4 - E/2)' It is worth noting that, apart from the different constant multiplying the ratio of Gamma functions on the right
hand side, the latter equation is of the same type of (3). However, the zeros of the entire function
5
*T(1/4 - E/2)
are now given by E = 2N + -, the eigenenergies pertaining to the antisymmetric states, so that for the particular
value £ = £c the degeneracy of the next lower antisymmetric state increases by one. Hence, we get level crossings
nr(1/4)
occurring for the unique value of the coupling parameter, namely £c —
9.29495. Therefore, the 3D
r(3/4)
point perturbation must be strongly attractive in order to exhibit the level crossing phenomenon.
In Fig. 3, the lowest level crossings between the new eigenvalues created by the point perturbation and the
antisymmetric levels E = -, -, —, clearly unaffected by the singular perturbation, are depicted. Therefore, also in the three-dimensional case, the effect is stronger than the Zeldovich effect that would occur in the presence of any short range interaction.
Fig. 3. The three lowest level crossings between the new eigenvalues created by the point perturbation and the antisymmetric levels E = -, -, —, clearly unaffected by the singular perturbation
E
8
6
4
2
0
-2
-4
10
-5
0
5
10
4. Level crossings in the discrete spectrum of the two-dimensional isotropic harmonic oscillator perturbed by a central ¿-interaction
In perfect analogy with the three-dimensional case, the two-dimensional self-adjoint Hamiltonian Hp representing the Hamiltonian of the two-dimensional isotropic harmonic oscillator perturbed by a central ¿-interaction, can be defined by means of the self-adjoint extension theory or by renormalization. Its discrete spectrum, in particular the lowest lying eigenvalues, has been investigated far less than the three-dimensional analog (see [19,20]).
As the intermediate steps would be identical to those for the three-dimensional case (see [7]), we omit them and write the limit as N ^ to in the norm topology of the resolvents, that is to say the counterpart of (4) in [7]:
(Hß - E)-1 = (Ho - E)-1 +
l*(E )><*(E)|
1 _ E ^ ^ Ki (°)^2n2 (0)
ß + 2n2 + 1)(2ni + 2n2 + 1 - E)
ni=0 n2=0 v ' K '
with:
oo oo
*(x,y; E) = ]T ]T
ni=0 n2 = 0
whose norm, for any E < 0, is bounded by:
02ni (O)02n2 (0) 2n1 + 2n2 + 1 - E
02 ni
-" <
=0 (2n + 2)
(13)
(14)
(15)
since the sequence 02n(0) behaves like n-1/4 as n ^ to.
Given that the eigenvalues of the operator Hp are exactly the poles of its resolvent, in order to achieve an accurate description of the spectrum (obviously being exclusively discrete) of the operator, we need only look for the solutions of the equation:
1
1 = eV
ß
CO CO
02ni (O)02n2 (0)
ß
(2ni + 2n2 + 1)(2ni + 2n2 + 1 - E)'
ni =0 n2 =0
(16)
The striking difference between the latter equation and its three-dimensional counterpart is that the right hand side cannot be recast as a ratio of Gamma functions; this can be seen as follows. By exploiting the well-known integral relationship between the resolvent and the semigroup for any semibounded operator for any E below the lowest point in the spectrum (see, e.g., [21] page 204), as well as the fact that the integral kernel of the semigroup of the two-dimensional harmonic oscillator is perfectly separable, we get that the right hand side of the above equation reads for any E < 1:
1 = 1, eH^it W e- - i) d =1
ß n J e2t - 1 n J 1 - e2t n ,
e-E -1 1 - e2
de.
(17)
It is not difficult to show that the last integral on the right hand side can be recast as a limit of an expression involving Gamma functions so that the bound state equation becomes:
a = — lim r(e)
r(1/2 - E/2)
r(1/2)
r(1/2 - E/2 + e) r(1/2 + e)
(18)
By exploiting the well-known identity:
r(e)r(1 - e) = ——-, sin(ne)
and the limits:
-in(ne) lim -= n,
e
equation (18) can be transformed into:
1 1
a = — lim -2n £
lim r(1 - e) = 1,
r(1/2 - E/2)
r(1/2)
_r(1/2 - E/2 + e) r(1/2 + e)_ ' (19)
The limit on the right hand side exists, which leads to an expression involving the difference of two different values of the logarithmic derivative of the Gamma function, namely:
1
r'(1/2) r'(1/2 - E/2)
r(1/2) r(1/2 - E/2)
1
(lnr)'(1/2) - (lnr)'(1/2 - E/2)
By noting that:
(lnr)'(x) = F(x - 1) = -y -
C
V ,
' \ n + x — 1
n=1
(20)
(21)
where y denotes the Euler-Mascheroni constant and F(x — 1) is the Digamma function, the latter expression can be recast as:
F(-1/2) - F(-1/2 - E/2)
E v
n=1
1
4^(n - 1) (n - W
(22)
1
1
n
Hence, the equation determining the perturbed eigenvalues of our Hamiltonian can be written either in terms of the difference of two different values of the Digamma function or in terms of a simple series (reducing to a telescoping one when E is any integer with the exclusion of all odd positive ones):
1
1
1
F(-1/2) - F(-1/2 - E/2)
E
1
4n n=o(n + i) (n + W
(23)
In particular, for any E given by a positive even number, the total energy of an antisymmetric state, the series on the right hand side becomes a simple telescoping series so that it is rather straightforward to find the location of the points where level crossings occur. For the lowest antisymmetric level with E = 2 we find = —n, whereas
3
for the upper one E = 4 we have £2 = — ^ n.
Figure 4 depicts the two lowest level crossings between the new eigenvalues of the perturbed Hamiltonian as functions of the strength of the point interaction and the two lowest antisymmetric levels of the 2D harmonic oscillator clearly belonging to the spectrum of perturbed Hamiltonian. As anticipated in the introduction, the first striking difference between Fig. 4 and Fig. 3, its three-dimensional counterpart, is that in the plot of the two-dimensional model the level crossings take place to the left of the origin, which implies that in two dimensions the point perturbation need only be weakly repulsive to produce the level crossing. Furthermore, each level crossing occurs at a different value of the coupling parameter.
-4
-2
Fig. 4. Plot of the two lowest level crossings between the eigenvalues created by the point perturbation and the two lowest antisymmetric levels of the harmonic oscillator which also belong to the spectrum of the perturbed Hamiltonian
E
4
2
0
P
2
4
0
2
4
5. Final remarks
The spectral phenomenon known as "level crossing of eigenvalues" has been investigated by considering the entirely discrete spectrum of four different energy operators: the Hamiltonian of the one-dimensional harmonic oscillator perturbed by a nonlocal ¿'-interaction, the Hamiltonian of the one-dimensional conic oscillator perturbed by the same singular interaction, the Hamiltonian of the two-dimensional harmonic oscillator perturbed by a ¿-interaction and its three-dimensional analog. All the point perturbations considered are highly singular so that the KLMN theorem used to define the ¿-perturbation in one dimension may not be invoked.
Although the level crossings in the one-dimensional models share with the three-dimensional model the property of occurring for a unique positive value of the coupling parameter (representing an attractive interaction), their nature is fundamentally different since in one dimension the eigenvalues pertaining to the antisymmetric bound states created by the singular perturbation cross the unperturbed even energy levels, exactly the opposite of what happens in three dimensions.
The two-dimensional model exhibits entirely different features; the crossings take place for different negative values of the coupling parameter so that the singular interaction need only be weakly repulsive to produce the crossings. From a mathematical point of view, the difference is due to the absence of a ratio of Gamma functions in the bound state equation and the appearance of a difference of Digamma functions as a consequence of the
logarithmic divergences of two-dimensional quantum mechanics. Further work to achieve a deeper understanding of the two-dimensional model is ongoing.
Acknowledgements
S. Fassari wishes to acknowledge Prof. I. Popov's kind invitation to present the main points of this article in one of the sessions of the conference "Mathematical Challenge of Quantum Transport in Nanosystems, Pierre Duclos Workshop", held at ITMO University, Saint Petersburg, Russia (September 26-27, 2017). Stimulating discussions with Prof. M.J. Lopez and Prof. J. Negro are kindly acknowledged. Partial financial support is acknowledged to the Spanish Junta de Castilla y Leon and FEDER (Project VA057U16) and MINECO (Project MTM2014-57129-C2-1-P). S. Fassari and M. L. Glasser also wish to thank the entire staff at Departamento de Física Teórica, Atómica y (Optica, Universidad de Valladolid, for their warm hospitality throughout their stays.
References
[1] Albeverio S., Gesztesy F., H0egh-Krohn R., Holden H. Solvable Models in Quantum Mechanics (AMS Chelsea Series), 2nd edition, American Mathematical Society, Providence, RI, 2004.
[2] Albeverio S., Kurasov P. Singular Perturbations of Differential Operators. Publishing House, Cambridge UK, 2000.
[3] Albeverio S., Dabrowski L., Kurasov P. Symmetries of Schrodinger operators with point interactions. Lett. Math. Phys., 1998, 45, P. 33-47.
[4] Exner P., Neidhardt H., Zagrebnov V.A. Potential approximations to ¿': an inverse Klauder phenomenon with norm-resolvent convergence. Commun. Math. Phys., 2001, 224, P. 593-612.
[5] Fassari S., Rinaldi F. On the spectrum of the Schrodinger Hamiltonian with a particular configuration of three point interactions. Rep. Math. Phys., 2009, 64, P. 367-393.
[6] Albeverio S., Fassari S., Rinaldi F. The Hamiltonian of the harmonic oscillator with an attractive ¿'-interaction centred at the origin as approximated by the one with a triple of attractive ¿-interactions. Journal of Physics A: Mathematical and Theoretical, 2016, 49, 025302.
[7] Fassari S., Inglese G. Spectroscopy of a three-dimensional isotropic harmonic oscillator with a delta-type perturbation. Helv. Phys. Acta, 1996, 69, P. 130-140.
[8] Brüning J., Geyler V., Lobanov I. Spectral properties of a short-range impurity in a quantum dot. J. Math. Phys., 2004, 45, P. 1267-1290.
[9] Albeverio S., Fassari S., Rinaldi F. Spectral properties of a symmetric three-dimensional quantum dot with a pair of identical attractive ¿-impurities symmetrically situated around the origin II. Nanosystems: Physics, Chemistry, Mathematics, 2016, 7 (5), P. 803-815.
[10] Albeverio S., Fassari S., Rinaldi F. Spectral properties of a symmetric three-dimensional quantum dot with a pair of identical attractive ¿-impurities symmetrically situated around the origin. Nanosystems: Physics, Chemistry, Mathematics, 2016, 7 (2), P. 268-289.
[11] Albeverio S., Fassari S., Rinaldi F. A remarkable spectral feature of the Schrodinger Hamiltonian of the harmonic oscillator perturbed by an attractive ¿'-interaction centred at the origin: double degeneracy and level crossing. Journal of Physics A: Mathematical and Theoretical, 2013, 46, 385305.
[12] Fassari S., Gadella M., Glasser M.L., Nieto L.M. On the spectrum of the Schrodinger Hamiltonian of the one-dimensional conic oscillator perturbed by a point interaction. ArXiv preprint 1706.04916, 2017.
[13] Fassari S., Gadella M., Glasser M.L., Nieto L.M. Spectroscopy of a one-dimensional V-shaped quantum well with a point impurity. Annals of Physics, 2018, 389, P. 48-62.
[14] Cheon T., Shigehara T. Fermion-Boson duality of one-dimensional quantum particles with generalized contact interactions. Phys. Rev. Lett., 1999, 82, 2536.
[15] Fassari S., Inglese G. On the spectrum of the harmonic oscillator with a delta-type perturbation. Helv. Phys. Acta, 1994, 67, P. 650-659.
[16] Fassari S., Inglese G. On the spectrum of the harmonic oscillator with a delta-type perturbation. II. Helv Phys. Acta, 1997, 70, P. 858-865.
[17] Fassari S., Rinaldi F. On the spectrum of the Schrodinger Hamiltonian of the one-dimensional harmonic oscillator perturbed by two identical attractive point interactions. Rep. Math. Phys., 2012, 69, P. 353-370.
[18] Lorinczi J., Ma lecki J. Spectral properties of the massless relativistic harmonic oscillator. J. Diff. Eq., 2012, 253, P. 2846-2871.
[19] Tusek M. Analysis of Two-Dimensional Quantum Models with Singular Potentials, Diploma Thesis, Czech Technical University, Prague, 2006.
[20] Tunali S. Point Interactions in Quantum Mechanics, MS Thesis, Izmir Institute of Technology, Izmir, 2014.
[21] Reed M., Simon B. Methods in Modern Mathematical Physics IV - Analysis of Operators, Academic Press, NY, 1978.
