CC BY
31
NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2016, 7 (2), P. 268-289
Spectral properties of a symmetric three-dimensional quantum dot with a pair of identical attractive 5-impurities symmetrically situated around the origin
S. Albeverio1,2’3’4, S. Fassari2,5’6 and F. Rinaldi2,6,7
4Institut fur Angewandte Mathematik, HCM, IZKS, BiBoS, Universitat Bonn, Endenicheralee 60, D-53115 Bonn, Germany 2CERFIM, PO Box 1132, CH-6601 Locarno, Switzerland 3Centre Interfacultaire Bernoulli, EPFL, CH-1015 Lausanne, Switzerland 4Chair Professorship, Department of Mathematics and Statistics,
King Fahd University of Petroleum and Minerals, Dhahran, KSA 5ISR, Aeulistr.10, CH-9470 Buchs, Switzerland 6Universita degli Studi Guglielmo Marconi, Via Plinio 44, I-00193 Rome, Italy 7BIS Group of Institutions, Gagra-Moga, Punjab under Punjab Technical University,
Punjab, India
silvestro.fassari@isr.ch
PACS 02.30.Gp, 02.30.Hq, 02.30.Hq, 02.30.Lt, 02.30.Sa, 02.30.Tb, 03.65.Db,
03.65.Ge, 68.65.Hb
DOI 10.17586/2220-8054-2016-7-2-268-289
In this presentation, we wish to provide an overview of the spectral features for the self-adjoint Hamiltonian of the three-dimensional isotropic harmonic oscillator perturbed by either a single attractive ^-interaction centered at the origin or by a pair of identical attractive ^-interactions symmetrically situated with respect to the origin. Given that such Hamiltonians represent the mathematical model for quantum dots with sharply localized impurities, we cannot help having the renowned article by Burning, Geyler and Lobanov [1] as our key reference. We shall also compare the spectral features of the aforementioned three-dimensional models with those of the self-adjoint Hamiltonian of the harmonic oscillator perturbed by an attractive ^'-interaction in one dimension, fully investigated in [2,3], given the existence in both models of the remarkable spectral phenomenon called ”level crossing”. The rigorous definition of the self-adjoint Hamiltonian for the singular double well model will be provided through the explicit formula for its resolvent (Green’s function). Furthermore, by studying in detail the equation determining the ground state energy for the double well model, it will be shown that the concept of “positional disorder”, introduced in [1] in the case of a quantum dot with a single ^-impurity, can also be extended to the model with the twin impurities in the sense that the greater the distance between the two impurities is, the less localized the ground state will be. Another noteworthy spectral phenomenon will also be determined; for each value of the distance between the two centers below a certain threshold value, there exists a range of values of the strength of the twin point interactions for which the first excited symmetric bound state is more tightly bound than the lowest antisymmetric bound state. Furthermore, it will be shown that, as the distance between the two impurities shrinks to zero, the 3D-Hamiltonian with the singular double well (requiring renormalization to be defined) does not converge to the one with a single ^-interaction centered at the origin having twice the strength, in contrast to its one-dimensional analog for which no renormalization is required. It is worth stressing that this phenomenon has also been recently observed in the case of another model requiring the renormalization of the coupling constant, namely the one-dimensional Salpeter Hamiltonian perturbed by two twin attractive ^-interactions symmetrically situated at the same distance from the origin.
Spectral properties of a symmetric three-dimensional quantum dot ...
269
Keywords: level crossing, degeneracy, point interactions, renormalisation, Schrodinger operators, quantum dots, perturbed quantum oscillators.
Received: 7 November 2015
1. Introduction
In this work, we wish to further extend our previous research on various types of point perturbations of Schrodinger Hamiltonians with or without harmonic confinement (see, e.g., [2-11]) to the physical models utilized to describe three-dimensional quantum dots.
This presentation will try to be an expansion of the detailed spectral analysis carried out in [1], in the sense that, after briefly reviewing in Section 2 the findings of [1,9] for the case of a three-dimensional quantum dot with a single impurity consisting of a point interaction centered at the origin, we are going to investigate in Section 3 the model with two twin attractive point impurities symmetrically situated with respect to the bottom of the harmonic confining potential. In Section 2, by analyzing the behavior of the new energy levels created by a point perturbation as functions of the parameter labeling the self-adjoint extensions, physically characterized by being proportional to the inverse scattering length (see [1,4]), we are going to stress the analogy of the presence of the so-called “level crossings” between the newly created eigenenergies of bound states with a given symmetry and those of the unperturbed harmonic oscillator unaffected by the singular perturbation pertaining to bound states having the opposite symmetry. We will also underline the difference between this three-dimensional model and the one-dimensional one, studied in [2,3], in which the harmonic oscillator is perturbed by a so-called ^'-interaction, whose quadratic form is given by the expression |d') (d'|.
The first step in the investigation outlined in Section 3 will necessarily be the rigorous definition of the self-adjoint Hamiltonian representing the energy operator for our model. This will be achieved by first introducing the usual ultraviolet energy cut-off and then making the coupling constant dependent on the cut-off itself in such a way that the typical cancellation of divergences will take place in the norm-resolvent limit once the cut-off is removed.
Once the explicit expression of the resolvent (Green’s function) is obtained, its poles will provide us with the eigenvalues (energy levels) for the three-dimensional harmonic oscillator perturbed by singular interactions (3D quantum dot with point impurities). The noticeable difference between the bound state equation pertaining to the model with a single impurity centered at the origin (described in Section 2) and the one for the singular double well is the inevitable absence of the simple expression containing the ratio of Gamma functions enabling us to study all the eigenvalues on an equal footing. Although it is still possible, at least conceptually, to carry out the analysis of the entire discrete spectrum, the operational task is made quite challenging by the fact that the equation for an eigenvalue pertaining to a bound state with a given symmetry may not be used for the eigenvalue pertaining to the next upper bound state with the same symmetry without the modifications needed to cancel the inevitable appearance of divergences, a procedure outlined in [2,8-11] in the case of the other aforementioned models.
As a consequence, we have necessarily restricted our task to a rather detailed description of the behavior of the lowest lying eigenvalues (energy levels) of the spectrum with respect to the variations of the two key parameters of the model, one labelling the self-adjoint extensions of the operator (or its reciprocal) and the other one given by the distance between either impurity and the origin.
Given that the findings of our analysis in Section 3 are unavoidably linked to some possible avenues of future research, we refrain from anticipating them in an abbreviated fashion at this stage, as they will be summarized in detail in Section 4, the final section in which
270 S. Albeverio, S. Fassari and F. Rinaldi
the most important conclusions will be drawn and some perspectives for future work will be sketched.
2. The three-dimensional isotropic harmonic oscillator perturbed by a single attractive point perturbation centered at the origin
We wish to start our analysis by revisiting the model used in [1] to describe a 3D-quantum dot, namely the three-dimensional isotropic harmonic oscillator perturbed by an attractive point perturbation. The spectroscopy of such a Schrodinger Hamiltonian was also fully investigated in [9] under the restriction of having the impurity, that is to say the point perturbation, situated at the origin. Hence, the Hamiltonian was given by the self-adjoint operator making sense of the heuristic expression (see [9]):
(2.1)
1 3 / d2 \ 3 3
He = H0 - в5(x), H0 = 2 ( -dX2 + x2) = a+ai + 2
i=1 ' * ' i=1
(ai and its adjoint being the usual annihilation and creation operators associated with H0).
The eigenfunctions of H0, the Hamiltonian of the three-dimensional harmonic oscillator, can clearly be expressed in terms of those of their one-dimensional counterparts, namely:
Фп(х) = Фт (х)фП2 {у)фпз (z).
Although, in this case, the radial symmetry of the model could obviously be exploited, as was done in [1], we adhere to the Cartesian framework used in [9] in view of our prospective investigation of the model in which H0 is perturbed by two attractive deltas. Before moving forward, it is crucial to notice that the function:
(Ho - E)-1(0,x)
E
|n|=0
Ф2т(0) (|2n| + § -
E
^2n(X)
(2.2)
is square integrable. In fact, by considering without loss of generality only negative values of E, by means of Parseval’s identity, we obtain the following estimate:
E
Ф|л(0)
|n|
=0 (|2n| + 2 - E)
,2 —
ф2й(0)
|n|=0 (2n1 + I) / (2n2 + D / (2n3 + 2)
П2/3 2 /
^ ф2га(0)
,П=0 (2n + 1/2)2/3
< to,
(2.3)
taking into account the well-known pointwise decay of the harmonic oscillator eigenfunctions (see, e.g., [12,13]), precisely that lim п1/4фт (x) < to, for any fixed x.
As was shown in [9], using the coupling constant renormalization procedure, H is rigorously defined (for any в different from zero) by means of its resolvent, namely:
(He - E)-1
(H0 - E)-1 +
| ( H0 -
в-1 -
E)-1(-, 0)><(H0 - E)-1(0,-)|
E
Vе
2^h
ф2я(0)
'|n|=° (|2n|+§)(|2n| + §-E)
(2.4)
for any negative E below E0 (в), the ground state energy (lowest eigenvalue) of Hp, that is to say the first root of the equation:
Spectral properties of a symmetric three-dimensional quantum dot ...
271
в-1
E
*2n(o)
1=0 d2nl + D d2ni +1 - E
(2.5)
as was done in [9] (and in similar contexts in [2,8,10,11]), the series on the right hand side of (2.5) can be recast as an integral. This can be achieved by taking advantage of two crucial properties that will be exploited throughout this work: the well-known integral relationship between the resolvent and the semigroup of any semibounded operator for any E below the lowest point in the spectrum (see, e.g., [14] page 204), as well as the fact that the integral kernel of the semigroup of the three-dimensional harmonic oscillator is perfectly separable (see [1] page 1278), such that (2.5) reads for any E < 3/2:
or equivalently,
в-1
С© о
1 f e2* (eEt - !) dt
V^J (eit - 1) 2 ’
0
в-1
e 2 t (eEt — l)
3
(1 — e-1t) 2
dt.
(2.5a)
(2.5b)
By means of a simple change of variable inside the latter integral, (2.5b) can be rewritten as an integral over a bounded interval (involving Mehler's kernel):
By setting a = l/в and
в-1
1
0
e1 (e-E -1) (i - e1)2
de.
a0 = в0-1
1
0
1 - e2 (1 - e1)2
de
Г(3/4)
n ■ Г(1/4)
^ 0.107585,
(2.5c) can be further transformed into:
(2.5c)
a - a0
1
0
e2-E -1 (1 - e2)3/2
de.
(2.5d)
In perfect analogy with what was done in the aforementioned articles, the integral on the right hand side can be converted into a ratio of Gamma functions as follows:
a0 - a
Г(4 - f) пГ( 1 - f).
(2.5e)
It is crucial to realize that (2.5e) holds also to the right of E = 3/2, enabling us to determine all the eigenvalues and not only the ground state energy. Of course, our analysis in this work is entirely consistent with what was observed in [9] (see the final remark in [2] as well). The graph of the ground state energy as a function of the extension parameter a = 1/в is shown below (Fig. 1). As is evident, E0(a) ^ 3/2 as a ^ +ro.
272
S. Albeverio, S. Fassari and F. Rinaldi
Fig. 1. The ground state energy of the 3D-isotropic harmonic oscillator perturbed by an attractive point interaction situated at the origin as a function of the extension parameter a = 1/в
At this stage, it might be worth recalling that the corresponding bound state equation for the one-dimensional Hamiltonian:
h(a) = hi/a = 1 ^-d-_ + x^ - a\8'(x)){8'(x)\,
defined rigorously by means of its resolvent in two different ways in [2, 3], is instead:
a0D — a
2-
r {%- E
a0D
2
Г(3/4)
0.675978.
Г (4 — E)’ ^0 “Г(1/4)
Conversely, the corresponding bound state equation for the one-dimensional Hamiltonian is:
(2-6)
h(X) = 2 ^—d— + — Ah(x) = 1 ^—— A|h(x))(h(x)1,
taking account of the obvious identity in the sense of quadratic forms between the ^-interaction, namely the quadratic form |b) (b|, and the so-called b-potential (whose quadratic form is instead (■, b-)), reads:
A
2r (4 — f) r( 4 — f),
(2-7)
(see [2,8,15]).
The fact that (2.7), unlike (2.5e) and (2.6), involves only the coupling constant is a clear manifestation of the fact that h(A) can be easily defined by means of the KLMN theorem (see [16]) without the renormalization required for the other two models.
The higher energy levels are even more interesting, as they exhibit the phenomenon called “level crossing”. As was pointed out in [2] (see also a similar remark in [9]): “as a
consequence of the presence of the point perturbation acting only on states with zero angular
7 (2l + 1) (2l + 2) . , ^ o; 3
momentum, the ------------------degeneracy of the eigenvalue E2t = 2l + - gets lowered by
one due to the emergence of the simple eigenvalue generated by the perturbation. Such a simple eigenvalue, regarded as a function of the extension parameter a, does cross the next
3
lower unperturbed eigenvalue E2l-i = (2l — 1) + - ”.
Spectral properties of a symmetric three-dimensional quantum dot ...
273
As can be seen in Fig. 2, it is quite noteworthy that all these level crossings take place exactly at a = a0 (approximately equal to 0.107585), adopting the renormalization used in [9]. As was stressed in [2]:
“by using instead the alternative renormalization (3.6), the location of all the level crossings would be exactly a(3^=0, leading to the graph shown in Fig. 4(a) of the aforementioned paper by Bruning, Geyler and Lobanov ”.
E
: . :
: —- ^ : / ■
: /
3-2-10 1
Fig. 2. The ground state energy and the next two symmetric eigenenergies of the 3D-isotropic harmonic oscillator with an attractive 5-perturbation situated at the origin as functions of the extension parameter a = Iff
The very same phenomenon observed in [2] for the one-dimensional Hamiltonian h(a), for which all the level crossings take place exactly at a = a[,1D) (approximately equal to 0.675978) using the renormalization (2.4a) of that paper, is to be regarded as being even more remarkable. Citing again our considerations in [2], it can be stated that:
“from the point of view of the structure of the spectral curves representing the eigenvalues (energy levels) as functions of the extension parameter and neglecting the degeneracy of the three-dimensional eigenvalues, the harmonic oscillator perturbed by the point interaction considered here (|5'(x))(5'(x)|) seems to be a more legitimate one-dimensional counterpart than the delta distribution as the latter bears no resemblance of the 3D-level crossing involving eigenstates of different symmetry. Having stated that analogy, an important difference must also be pointed out: whilst in the case of the perturbed isotropic oscillator the eigenenergy of the simple eigenvalue created by the point interaction (emerging out of a degenerate level with an even value of the total angular momentum) can cross the next lower unperturbed level (having an odd value of the total angular momentum) beyond a certain threshold of the key parameter of that model, something of an opposite nature occurs in the one-dimensional model being studied here: each perturbed odd eigenvalue can fall below the next lower unperturbed even eigenvalue beyond a certain threshold. As a consequence, the symmetry of the ground state wave function can change in the case of a sufficiently strong 5'-interaction (the ground state wave function being given by an odd function discontinuous at the origin)”.
The spectral features of the Hamiltonian for our perturbed isotropic oscillator described in the previous remarks can be visualized in the above graph (Fig. 2) displaying the new ground state energy E0(a) and those of the two symmetric bound states E2(a) and E4(a) which emerge from the 7f 2 and 11/2 levels as a result of the perturbation and are confined inside the intervals (3f2, 7f2) and (7f2,11f2), as functions of the extension parameter, as well as their level crossings with the unperturbed antisymmetric states E = 5f2 and E = 9f2. Both of these
274
S. Albeverio, S. Fassari and F. Rinaldi
occur, as anticipated above, at a = ao (approximately equal to 0.107585). It is worth reminding the reader that the degeneracy of E = 7/2 has been lowered to five and that of E =11/2 to fourteen.
3. The three-dimensional isotropic harmonic oscillator perturbed by two twin attractive point perturbations symmetrically situated with respect to the origin
We wish to continue our analysis of the three-dimensional model by considering the heuristic Hamiltonian:
H{e,x0} = Ho - в [£(x - Xo) + S(x + Xo)] , xo = (xo, 0, 0), x0 > 0. (3.1)
As pointed out in the final remarks in [17] (see also the references therein), the model with a singular double well has been increasingly utilized in the recent condensed matter theory literature “since, in combination with a cubic nonlinearity, it leads to a particular type of the nonlinear Schrodinger equation for Bose-Einstein condensates (Gross-Pitaevskii equation)”.
By essentially mimicking what was done in [2, 9], we can start by performing the standard ultraviolet cut-off and making the coupling p = p(t,f) > 0 dependent upon the related upper limit 0 < t < for the energy, which leads to the introduction of a perfectly defined self-adjoint Hamiltonian:
t
H{£j3,x0} = Ho - p(t,e) ^ |ФГ) [ФГ(xo) Фт (xo) + ФГ(-xo) Фm (-xo)] (Фл|, (3.2)
|zl=o,|m |=o
which implies that for any sufficiently negative E we can write:
Hft,s,*0} - E =(Ho - E)1/2{l - p(t,e)[|(4t - E)-1/2(-, -x0))((4t - E)-1/2(-f0, -)| + \(H‘ - E)-1/2(-,xo))((Ht - E)"1/2(fo, -)|]}(Ho - E)1/2.
(3.3)
Therefore, in order to explicitly write the resolvent of the left hand side of (3.2), it is crucial to determine the inverse of the operator inside the curly brackets on the right hand side of (3.3). We wish to remind the reader that the expression of the resolvent given by the inverse of (3.3) is known as Tiktopoulos’ formula, commonly used in dealing with perturbations of semibounded operators (see, e.g., Theor. VIII.25 in [18]). It is worth noting that in the latter expansion of the resolvent the crucial role is played by an operator of the type:
(A + I)-1/2(An - A)(A + I)-1/2, A > 0,
which is isospectral to the more widely investigated Birman-Schwinger operator (see, e.g., [19, 20] in the case of perturbations of the negative Laplacian in one dimension). In perfect analogy with what was done in [7,11,17], it is possible to rewrite the operator of rank two inside the square brackets on the right hand side of (3.3) as follows:
|(H - E)-1/2(-, -xo)) ((Ht - E)-1/2(-xo, -)| + |(Ht - E)"1/2(-,fo)) ((Ht - E)-1/2(xo, -)| =
2[|(Н - E);1/2(-,xo)) ((H - E)-1/2(xo, -)| + |(Ht - E)-,1/2(.,fo)) ((H - E)-,1/2(fo, -)|],
(3.4)
where
Ф2n(Xo )
(Ht - E);1/2(fo,x):= 5]
|n|
=o (|2n| + 2 - E)
1/2
Ф2П (x),
(Ho - E)-1/2(xo,x) :=
Spectral properties of a symmetric three-dimensional quantum dot ...
275
£
f2n1 + l(Xo)f2n2 (0)^2ns (0)
T/2 ^2ni+1 (x)f2 n2 (У)^2пз (z) ,
|n|=0 (2n1 + 2n2 + 2n3 + § — E)
which implies their mutual orthogonality.
By essentially mimicking what was done in the aforementioned papers, it is rather straightforward to obtain the algebra for the powers of the rank-two operator on the right hand side of (3.4):
m+1
2 [|(H - £)-1/2МоЖ(Я0 - E)-1/2(xo,-)| +
|(H> - E)-,1/2(.,fo)X(H0 - E)-,1/2(fo, ■)
2m+1 { [(Я0 - E)-1 (fo,Xo)]m |(H0 - E)-1/2(-,fo)> ((H‘ - E):1/2(lo, 0| +
[(H - E);J1(xo,?o)]m|(H0 - E)-,1/2(.,xo)>{(Ho' - E)-1/2(xo,-)|
given that:
(3.5)
2
(H0 - E)-1/2(xo, ■) = £
^§n (xo)
|n|=0
|2n| + § - E
(H0 - E)g 1 (Xo,Xo)
1 /2
and similarly for the square of the norm of (HO - E)as (X0, ■). Hence, the inverse of the operator inside the curly brackets of (3.3) can be written as a new rank-two operator:
b) - (H0 - E)/ (x0,x0)
2^m
i
(H0 - E)-1/2 (•,.*,))((H0 - E)-1/2 (X0, ■)
(H0 - E )-s1/2 (',.x,))((H0 - E )-1/2 (■*. ■)
+
2^(Г,в) (H0 E)as (x0,x0)
Therefore, the inverse of the operator on the left hand side of (3.3) is:
1
(3.6)
(H{Wo} - E) 1 = (H0 - E)-1 +
i
2М(Г„3)
- (h0 - e)e 1 (X0,x0)
i
+
1
2М(Г,в)
- (H0 - E)J (x0,x0)
(H0: - E)-1 (.,f0^((H0 - E)-1 (X0, ■)
(h0: - e )-,1 (•,.x=))((H0: - e )-; (.as ■)
(3.7)
At this stage, we need to thoroughly investigate what occurs when the ultraviolet cut-off gets removed, that is to say the behavior of (3.7) as t ^ +<ro.
As a consequence of the simple estimate:
£
^2n(x0)
,2 —
£
^2n(x0)
|n| =
=0 (|2n| + § - E) |n|=0 (|2n| + § - E)
2—
18 IT о 1 18' 3 I
_n=0 (2n + 1/2)2/3_ _n=0 (2n + 1/2)2/3_
< to, (3.8)
i
1
2
276
S. Albeverio, S. Fassari and F. Rinaldi
it is immediate to realize that — E) s (x°, ■) will converge in the Hilbert space norm to the
following:
(H° — E)s 1 (x0,x) := Y »2+(3X°— EФ2n(x).
(3.9)
|Й|=° 1 1 2
The latter convergence implies in a straightforward manner that the rank-one operator:
(я0 — e);1 (•,.х°)}((я0 — e);1 (.x°,■),
will converge in the trace class norm to the following:
|(H° — e);1 (•,.x°)X(H° — e);1 (,x°, .)|.
Similarly, it follows that:
(h° — e);,1 (;x„))((h° — e);,1 (x°,-)| ^ |(h° — e);1 (;x°)X(h° — e);,1 (x°,-)|
with
(H° — EL1 (x°,x) := Y ( Eni+1(xo)En2 (0)^2пз (°)^2rai+1(x)^2n2 (У)ф2пз (z). (3.10)
|й|=0 (2П1 + 2П2 + 2Пз + § — E)'
Therefore, the only problem caused by the removal of the ultraviolet cut-off is the divergence of the series inside each denominator in the second and third term on the right hand side of (3.7). However, we still have the freedom of choosing the £-dependence of у(£,в) in such a way that the typical “cancellation of infinities” may take place. Let us set:
^_ = 1 , v ф2п(°)
М£,в) в n=0 |2n| + 3 ’
(3.11)
or equivalently:
M£,e) = в
i+в Y, *2n(0)
|П|=°
|2n| + 2
-1
(3.11b)
In perfect accordance with the use of the term “attractive” in [2, 3, 9, 17], it is clear that у(£,в) > ° for the large values of £ involved in the limit, regardless of the sign of в. This makes the singular interaction attractive because of the presence of the negative sign in the second term in (3.2). Hence, for any E < 3/2:
1 — (h° — E )s 1 (йА) = Ф +1Y
ф|п(°)
£
Ф2п(х°)
2^(£,в)
2в' 2 |П=°|2n| + 2 |П|=°|2n| +3 — E
(3.12)
and
1
2^(£,в)
(H° — E)sl (x° ,x°) =
1 1 1 ^n(0)
2в + 2 n^° |2n| + 3
|n|=^ 2
£
E
|n|=°
ф(х°)2га1+1ф2п2 (0)ф2пз (0)
2n1 + 2n2 + 2пз + I — E
(3.12a)
(as will be seen shortly, considered separately, the latter expression is well defined for any E < 5/2).
By taking advantage of the integral representation for the resolvent of the three-dimensional harmonic oscillator for any E < 3/2, that is to say the aforementioned Mehler's kernel, as well
Spectral properties of a symmetric three-dimensional quantum dot ...
277
as of the explicit expressions of the symmetric and antisymmetric components of the integral kernel of the semigroup for the one-dimensional harmonic oscillator established in [11], the limit of the right hand side of (3.12), as t ^ +<ro can be written as:
1 1
— +
2в 2n3/2
£2
(1 - £2)3/2
d£ -
1 E
г £2 E
' _T2i-i _x2i+?
e x0 1+g + g x0 1-g
(1 - £2)3/2
d£
1 1 r £1/2 1 - £-E (g-x01+« + g-x01-
+ 2n3/2 J
..2 i-g _21+;
(1 - £2)3/2
d£ < to.
(3.13)
Similarly, for any E < 5/2, the limit of the other denominator is given by:
1 1 r £1/2 1 - £-E (g-x21- - g-x21-
— +
2f 2n3/2
(1 - £2)3/2
d£ < to,
(3.13a)
1
taking into account that, as pointed out in [11], we have:
g
-x
2 1-g 0 1+g
g
-x
2 1+g
0 1-g
—> 4x
2
0,
as £ ^ 0+, which implies that £1/2 E as long as E < 5/2.
g
-x
2 1-g 0 1+g
g
-x
2 1+g
0 1-g
has an integrable singularity at 0+
Hence, for any x0 > 0 and any E < 3/2 (resp. E < 5/2), the difference of the two divergent series on the right hand side of (3.12) (resp. (3.12a)) yields a finite limit. Then, it is not difficult to prove that, as a consequence of the first resolvent identity, the same holds at any other point away from a(H0).
by:
Therefore, the norm resolvent limit of (3.7) after removing the ultraviolet cut-off is given
(Ho - E)-1 +
1
2в
(Ho -
+ lim
l^ + <^
E)- (-,Хо))((Яо - E)- (xo,-)|
1 ^2n(0) *2»(X0)
2 |Пт=0|2n| +3 |Пт=0 l2n| +3 - E
(H - E)-sr (•■Xq^((Hq - E)- (X0, •)
1
2в
+ lim
l^ + <^
1 Ф2га(0)
2 |i=0 l2nl + 1
l
E
|n|=0
^jn1+1(x0)^in2 (0)^2ra3 (0) 2nr + 2n2 + 2n3 + I - E
(3.14)
278
S. Albeverio, S. Fassari and F. Rinaldi
In particular, for any E < 3/2, we have the following integral representation:
(Ho - E)-1 +
(Ho - E)- (;Xo)){(Ho - E)- (xo, •)
1 1 ^ C1/2 [l - (e-x2l+ + e-x§1-t
2e + 2П3/2,/
dC
(1 - C2)3/2
I (Ho - E)- Mo))((Ho - E)- (Xo,-)|
1 1 1 C1/2 [l - C-E {e-x2°^ - e-x2i-
2в + 2n3/2 J
(3.14a)
(1 - C2)3/2
dC
By essentially mimicking the proofs used in some of the above-mentioned papers [2,9, 17], it can be rigorously shown that the latter operator is indeed the resolvent of a semibounded self-adjoint operator for any E below Eo(5), the lowest zero of the denominator in the second term in (3.14) and (3.14a). Furthermore, the explicit form of the resolvent clearly implies its analyticity in norm in close proximity to 0 (since we have an analytic family in the sense of Kato, see [14,21]) as a function of в.
The results obtained so far can thus be summarized in the following theorem.
Theorem 3.1. The rigorous Hamiltonian of the three-dimensional isotropic oscillator, perturbed by two identical attractive point interactions situated symmetrically with respect to the origin at the points ±Xo = (±xo, 0, 0), xo > 0, making sense of the merely formal expression (3.1), is the self-adjoint operator H{e,Xo} whose resolvent is given by the bounded operator (3.14). The latter is the limit of the resolvents (3.7) (of the Hamiltonians with the energy cut-off defined by (3.2)) in the norm topology of bounded operators on L2(R3) once the energy cut-off is removed. Furthermore, H{e,Xo}, regarded as a function of в, is an analytic family in the sense of Kato.
Of course, unlike the Hamiltonian analyzed in Section 2, we can no longer hope to have the bound state equation expressed in terms of a ratio of Gamma functions.
Nevertheless, since the limit of the difference of the divergent series appearing in the first (resp. second) denominator of (3.14) gives rise to a smooth multibranch function of the energy parameter, with its vertical asymptotes given by the symmetric (resp. antisymmetric) eigenvalues of the unperturbed Hamiltonian, we can rely on the findings of [8-11] to have it expressed in terms of an integral that will have to be suitably modified for each single level in order to avoid the unpleasant appearance of divergences.
Remark. In perfect accordance with what was observed in [10] for the analogous one-dimensional model with a single point impurity centered away from the origin, and in [11] for the model having a pair of twin impurities symmetrically situated with respect to the origin, it is worth pointing out that whenever the location of the impurities coincides with a node of an eigenfunction for the one-dimensional harmonic oscillator (here we are obviously excluding the trivial case xo = 0), then the corresponding eigenenvalue will not be affected at all by the singular perturbation. In other words, if
^2mi (xo) = 0, xo > 0 (respectively ^mi+^xo) = 0),
then
^2m (x), m =(ш1,Ш2,Шз) (resp. f2mi+1(x)^2m2 (y)f2m3 (z )t)
Spectral properties of a symmetric three-dimensional quantum dot ...
279
belongs to the kernel of the first (resp. second) rank-one operator in (3.14) and the corresponding eigenvalue
, 3
\2m \ + ^
resp. 2mi + 2m2 + 2m3 +
5
2
will stay in the spectrum with its original multiplicity.
Hence, for any E < 3/2, the equation determining the ground state energy is given by:
1 1 j f1/2 f-E (e-x21+g + e-x21—- 1
в = n/2 J (1 - f2)3/2
df.
(3.15)
By setting a = 1/в and
a+(xo)
1
i ! —x
2 1—g
x21+g
e x0 1+g + e xo 1+g
) - f1/2
n3/2
(1 - f2)3/2
df < to,
(3.15a)
the equation can be exactly recast in a manner that is completely analogous to (2.5d), that is to say:
a - a+ (x0)
1
(f1/2
E
n3/2
i
1
' x2 1+g
e Xo 1+« + e
-x
2 1+g
0 1—g
(1 - f2)3/2
df.
(3.15b)
The graph for the solution of the above equation, i.e. the ground state energy as a function of the extension parameter, is shown below in Fig. 3 for x0 = 0.2.
E
r
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4. /' • y’ . / • /' • /
0.0 0.5 1.0 1.5
Fig. 3. The ground state energy for the 3D-isotropic harmonic oscillator perturbed by a pair of identical attractive point interactions symmetrically situated at a distance x0 = 0.2 from the origin as a function of the extension parameter
a = 1/в
Similarly, the equation determining the energy of the lowest antisymmetric bound state for any E < 5/2, namely:
1
в
1
п3/2
f1/2 f
_ rp t _x21+i _x21±g
--E I e xo 1+g e xo 1—g
(1 - f2)3/2
i
1
Adf,
(3.16)
0
280
S. Albeverio, S. Fassari and F. Rinaldi
after setting
can be recast as:
ao Ы
1
..2 1 —
x2 1+i
1 1 e xoi+1 — e X0 1 —
) - i1/2
n3/2
(i - i 2)3/2
di < ж,
(3.16a)
a — a0 (x0)
1
(i1/2-E
— D(
_х2 1— x21+;
e X0 1+5 _ e X0 1-5
П3/2
(i — i2)3/2
di.
1
(3.16b)
0
The graph for the solution of the above equation, representing the energy of the lowest antisymmetric bound state as a function of the extension parameter, is shown below in Fig. 4 for
x0 = 0.2.
FIG. 4. The energy for the lowest antisymmetric bound state of the 3D-isotropic harmonic oscillator perturbed by a pair of identical attractive point interactions symmetrically situated at a distance x0 = 0.2 from the origin as a function of the extension parameter a = 1/в
By further exploiting the technique developed in the aforementioned papers [8-11] to obtain cancellation of the divergences as we move to the right of E = 3/2, we can also study the equation determining the second symmetric bound state energy.
Given that:
lim
1 v ф2п(°)
2 n0 |2n| + §
|n|=0 2
1Ф0 (0)
2 3/2
i
E
|n|=0
^L(xq)
|2n| + 3 — E
+ lim 1 T, T^2^
3/2 — E i^+~ 2 |2n| + 2
|n|=^ 2
Ф2п(Х0)
№=1 |2nl + I — E
and
фр(0)
3/2
1
n3/2
1
f i 1/2di,
Ф0Ы
3/2 — E
x2 e Xo
n3/2
1
i1/2-E di
(3.17)
0
0
Spectral properties of a symmetric three-dimensional quantum dot ...
the equation determining the second symmetric bound state energy is given by:
1 £1/2 [(i - a3/2 -1]
281
2g-x°
+
в (3/2 — E) n3/2 3n3/2 n3/2
(1 — £2)3/2
_d£+
1
£1/2-E e-xg(1-f)/(1+f) + g-x°(1+?)/(1-?) — 2e-xg (i — £2)3/2
П3/2
d£. (3-18)
o (1 — £2)3/2 Since the quantity
g-*§(1-0/(1+0 + g-x°(1+?)/(1-?) — 2e-x° ^1 — £2)3/2
behaves like c£2 + O(£4) in the region to the right of the origin, it is quite evident that the numerator inside the second integral on the right hand side of (3.18) has an integrable singularity at the lower limit of integration for any E < 7/2, that is to say the third eigenvalue of the unperturbed 3D-isotropic harmonic oscillator. One can immediately notice that both integrals on the right hand side of (3.18) diverge as £ ^ 1-. However, by rewriting their sum as:
1
£1/2
£2)3/2 — 1 — £ E |~2e-x° (1 — £2)3/2 — g-x§(1-£)/(1+£) — g-*°(1+0/(1-0
П3/2
(1 — £2)3/2
u
it is a bit tedious but straightforward to check that:
£1/2 j (1 — £2)3/2 — 1 — £-E 2g-x2 (1 — £2)3/2 — g-x°(1-?)/(1+?) — g-x°(1+?)/(1-?) j
d£,
—У Le < oc,
1 — £2
as £ ^ 1-, which implies the integrable nature of the singularity at the upper limit of integration and, therefore, the cancellation of the two divergent quantities on the right hand side of (3.18). Hence, for any 3/2 < E < 7/2, the equation determining the second symmetric bound state energy can be recast as:
2g-x°
21
+ ^7 x
1 _
в (3/2 — E) n3/2 3n3/2 ' n3/2
1 £1/2 j (1 — £2)3/2 — £-E 2g-x° (1 — £2)3/2 — g-xg(1-?)/(1+?) — g-xg(1+£)/(1-£)
J (1 — £ 2)3/2
0 V 7
The latter equation can be further transformed into:
d£. (3.18a)
a — a+(x0)
2g-x°
21
+ x
(3/2 — E) n3/2 3n3/2 n3/2
£ 1/2 j (1 — £2)3/2 — £-E 2g-x§ (1 — £2)3/2 — g-x°(1-0/(1+0 — g-x°(1+?)/(1-?) |
(1 — £ 2)3/2
— g-x°(1-?)/(1+?) — g-x§(1+S)/(1-S) (1 — £—)3/2
d£+
(3.18b)
Figure 5 shows the plot for the counterpart of Fig. 2 for x0 = 0.2 with the five lowest eigenvalues, namely E0(a), E1(a), E2(a) (clearly omitting their dependence on x0 to make the
1
1
1
1
282
S. Albeverio, S. Fassari and F. Rinaldi
notation less cumbersome), the new eigenenergies created by the pair of twin point interactions, in addition to 5/2 and 7/2, whose degeneracy has now been lowered to two and five respectively.
Fig. 5. The ground state energy and the next two eigenenergies created by two identical attractive ^-perturbations situated symmetrically at the same distance from the origin, x0 = 0.2, as functions of the extension parameter a = 1/в
It is worth stressing that, in addition to level crossings of the same type as those encountered in the aforementioned articles [1,9], which occur in the case of a single 4-impurity, involving the new eigenvalues created by the point perturbations and the eigenvalues of the unperturbed harmonic oscillator preserved under the perturbation due to their degeneracy, two further level crossings of a new type involving Ei(a) and E2(a) (the two new energy levels created by the singular double-well perturbation above the ground state one), can be observed at the values a2 (approximately equal to -0.126478) and a3 (approximately equal to 0.309201). It seems rather remarkable to us that, for such a small value of the distance x0, there is a range of values for the extension parameter a, namely the interval (a2, a3), inside which E2 (a) < E1(a). In other words, when the coupling is sufficiently strong (taking account of the fact that a = 1/в) and the distance between the two centers is sufficiently small (this point will be made more precise shortly), the second symmetric bound state created by the singular perturbation can become more tightly bound than the lowest antisymmetric bound state. Furthermore, the energy level created by the perturbation at a2 (resp. a3) is doubly degenerate, since E1(a2) = E2(a2) (resp. E1(a3) = E2(a3)).
In addition to the value a4 (approximately equal to 0.356236), where E2(a) = 5/2, it is also worth including the value a1 (approximately equal to -0.198151) where E1(a) = 3/2, even though E = 3/2 is no longer an eigenvalue. Then, the spectral features displayed in Fig. 5 can be summarized in Table 1.
Next, we will show how an increase in x0 affects the energy levels considered before. As clearly shown by the comparison between the graph of the ground state energy as a function of a in Fig. 6 for x0 = 0.4 and its analogue in Fig. 3 for x0 = 0.2, we have E0(a; 0.2) < E0(a; 0.4) for any a.
Although we cannot strictly use the definition of “positional disorder ” given in [1] for the 3D-isotropic quantum dot with a single impurity because of the presence of the vector -x0 in our model, we, nevertheless, believe that the concept is still valid in the sense that the greater the distance between the two impurities is, the less localized the ground state will be. It might be worth noting that this is in full agreement with our findings in some of the aforementioned papers [7,11,17] on one-dimensional singular double wells (see also [10]).
Spectral properties of a symmetric three-dimensional quantum dot ...
283
Table 1. The ordering of the five lowest eigenvalues over the corresponding five intervals of the extension parameter for x0 = 0.2
a < ai E0(a) < Ei(a) < 3/2 < E2(a) < 5/2
ai < a < a2 E0 (a) < 3/2 < Ei(a) < E2 (a) < 5/2
a2 < a < a3 E0 (a) < 3/2 < E2(a) < Ei (a) < 5/2
a3 < a < a4 E0 (a) < 3/2 < Ei(a) < E2 (a) < 5/2
a4 < a E0(a) < 3/2 < Ei(a) < 5/2 < E2(a) < 7/2
E
- _______—:
" / : / " / / . / ■ /
•0.2 0.0 0.2 0.4 0.6 0.S 1.0
Fig. 6. The ground state energy of the 3D-isotropic harmonic oscillator perturbed by a pair of identical attractive point interactions symmetrically situated at a distance x0 = 0.4 from the origin as a function of the extension parameter
a = 1/в
The graph of the lowest antisymmetric eigenvalue created by the singular double well when the separation distance becomes x0 = 0.4 is also provided below (Fig.7).
We now wish to provide the analog of Fig. 5 for x0 = 0.3 (Fig. 8).
As can be clearly seen in the plot, a remarkable change in the energy spectrum has now occurred: for any value of the extension parameter we have Efa) < E2(a), that is to say, the new energy levels created by the double well perturbation completely avoid each other, borrowing the terminology used in articles such as [22] on the spectral phenomenon called “level repulsion” or “avoided level crossing”. Of course, we still have the level crossing between E2(a) and 5/2, the antisymmetric eigenvalue of the 3D-isotropic harmonic oscillator whose degeneracy has now been lowered to two by the perturbation.
The threshold value of x0, denoted by Xt in the following, can be determined by observing that for that particular value, the two spectral curves Ei(a,Xt) and E2(a,Xt) must have the same value as well as the same partial derivative with respect to the extension parameter a. Hence, Xt is the second coordinate of the point representing the solution of the following
284
S. Albeverio, S. Fassari and F. Rinaldi
Fig. 7. The energy for the lowest antisymmetric bound state of the 3D-isotropic harmonic oscillator perturbed by a pair of identical attractive point interactions symmetrically situated at a distance x0 = 0.4 from the origin as a function of the extension parameter a = 1/в
Fig. 8. The ground state energy and the next two eigenenergies created by two attractive b-perturbations situated symmetrically at the same distance from the origin, x0 = 0.3, as functions of the extension parameter a = 1/в
system:
Ei(a, xo) = E2(a,xo);
д д
—Ei (a, xo) = —E2(a,xo).
(3.19)
For any fixed value of the distance x0, the function El(a, x0) with domain given by the entire real axis and range (-го, 5/2) has its inverse given by al(E,x0) with domain (-го, 5/2) and range equal to the entire real axis.
Similarly, if we restrict E2(a,x0) to the domain (-го, a*] and range (-го, 5/2], with E2(a*,x0) = 5/2, its inverse function is a2(E,x0) with domain (-го, 5/2] and range (-го, a*]. Both functions can actually be written in terms of explicit integrals and functions:
ai(E,xo)
1
1 £1/2 E-e L-x2ok+l - e-x0T-|
1
n3/2
(1 - £2)3/2
d£,
Spectral properties of a symmetric three-dimensional quantum dot ...
285
a2(E,xo) r t1/2
2
2e-x°
+
(3/2 — E) n3/2 3n3/2 n3/2
_ t2)3/2 — 1 — t-E 2e-x° (1 — 12)3/2 — e-x2(1-?)/(1+?) — e-x°(1+?)/(1-?)
(i —12)3/2
dt
2
1
(of course, we could alternatively use (3.16b) and (3.18b)).
Consequently, also taking into account the well-known theorem on the derivative of the inverse function, instead of obtaining Xt by solving the system (3.19), it is easier from a computational point of view to derive its value by solving the system:
aq(E,Xo) = «2 (E,x o);
d_
dE
a1(E,Xo)
d_
dE
a2(E,Xo).
(3.19a)
The numerical solution of (3.19a) is the point with coordinates E = Et (approximately equal to 2.29949), x0 = Xt (approximately equal to 0.2339644), where both functions a1 and a2 attain the value at (approximately equal to 0.0829). The plot of the tangential contact between the two spectral curves for xo = Xt is provided below in Fig. 9.
Fig. 9. The curves of the two eigenenergies E1(a,Xt) and E2(a,Xt) (Xt being approximately equal to 0.2339644) intersecting each other tangentially at a = at (approximately equal to 0.0829)
Hence, what has been seen regarding the three different spectral configurations for the two energy levels E1(a,x0) and E2(a,x0) can be summarized by means of the following Table 2, which clearly shows the role played by the threshold value x0 = Xt.
The table thus makes our earlier statement written below Fig. 5 more precise: for any x0 < Xt, there exists an interval (a2(x0), a3(x0)) (using the same indices used above) of values for the extension parameter a (or equivalently of its reciprocal в) for which E2(a,x0) < E1(a,x0). In other words, as the values of |a| in that range are rather small (a2(0.2) approximately equal to —0.126478 and a3(0.2) approximately equal to 0.309201 in our previous example) and the corresponding values of |в | are large, such a singular configuration of those two eigenvalues (energy levels), characterized by having the second excited symmetric
286
S. Albeverio, S. Fassari and F. Rinaldi
Table 2. The three possible configurations involving the two energy levels Ei(a,x0) and E2(a,x0) as the distance parameter x0 varies
x0 < Xt = 0.2339644 two distinct intersections between Ei(a,x0) and E2(a,x0)
x0 = Xt = 0.2339644, at = 0.0829 one tangential intersection between Ei(a,x0) and E2(a,x0)
x0 > Xt ^ 0.2339644 no intersection since Ei(a,x0) < E2(a,x0)
bound state more tightly bound than the lowest antisymmetric one, can only occur by shrinking the distance between the two impurities below the threshold and raising their strength. Moreover, it can be seen that the interval expands as x0 shrinks to zero.
Therefore, it might be fair to say that a wider separation of the two zero-range impurities makes the structure of the spectrum more similar to that of the unperturbed isotropic harmonic oscillator. As a matter of fact, it is not difficult to convince oneself that, for any fixed value of the coupling в (or, equivalently, of the extension parameter a), the denominators of both rank-one operators in (3.14) or (3.14a) are bound to diverge as x0 ^ +rc>. Hence, the resolvent of our Hamiltonian will approach that of the 3D-isotropic harmonic oscillator in the limit x0 ^ +rcx As a result, that is in perfect agreement with its one-dimensional counterpart established in Theorem 2.2(a) in [11].
However, a sharp contrast with the one-dimensional model occurs for the opposite limit analyzed in Theorem 2.2(b) in [11]. In fact, by noting that the right hand sides of both (3.15) and (3.16) become divergent as x0 ^ 0+, it is clear that the lowest eigenvalue of our Hamiltonian will decrease without any lower bound as the distance shrinks to zero. Therefore, in sharp contrast with the one-dimensional case, we cannot expect in the limit x0 ^ 0+ to obtain the operator studied in Section 2 with the strength of the point interaction doubled. It might also be worth reminding the reader that, as is well known to Quantum Chemistry students, three-dimensional interactions with a nonzero range do not manifest this singular behavior as the distance between the two centers shrinks to zero, as the classical textbook example of H+ smoothly approaching He+ in the limit R ^ 0+ clearly shows (see [23-25]).
The same singular spectral phenomenon has been recently observed in [17] where another model requiring a renormalization procedure, that is to say the one-dimensional Salpeter Hamiltonian with a double well of identical attractive point interactions, was investigated. Our spectral finding in this work confirms the general conjecture formulated in the above-mentioned paper that, when we deal with point interactions requiring the renormalization of the coupling constant (either because of the higher dimension or because of the different kinetic term), we must inevitably give up “the intuitive idea that in the limit (x0 ^ 0+) the eigenvalue should approach the ground state energy pertaining to a single delta with double strength ”.
In analogy with what was proved in [17], it will be shown in a forthcoming paper that the only possible way of obtaining the expected behavior in the limit also in the case of the operator investigated in this work is to make the coupling constant gif, в) used in the renormalization procedure suitably dependent on the separation distance x0.
Spectral properties of a symmetric three-dimensional quantum dot ... 287
4. Conclusions and perspectives
After reviewing the spectral features of a three-dimensional isotropic quantum dot with an attractive point impurity situated at the origin (the bottom of the confining harmonic potential), a model investigated in papers such as [1,9,26], the model in which a pair of twin attractive point impurities gets inserted into the 3D-isotropic quantum dot has been studied by first defining in a rigorous manner the self-adjoint Hamiltonian representing the energy operator for such a model. Since the energy eigenvalues are given by the poles of the resolvent, its explicit expression has enabled us to start the spectral analysis of the model. Although the bound state equation no longer has the simple expression involving a ratio of Gamma functions, as was the case in the presence of a single impurity centered at the origin, it is still possible to carry out, at least conceptually, the required spectral analysis.
Here, we have restricted our task to a rather detailed description of the behavior of the lowest lying spectrum eigenvalues (energy levels) with respect to variations of the two key parameters of the model, one labeling the self-adjoint extensions of the operator (or its reciprocal) and the other one given by the distance between either impurity and the origin.
The important spectral phenomenon called ’’positional disorder” singled out in [1] (see also the references cited therein as well as [10]) in the case of a single point impurity located away from the origin retains its validity in the sense that the more separated the twin impurities are the less localized the ground state is.
In addition to the level crossings involving the spectral curve of each new singular perturbation-created eigenenergy and the next lower eigenvalue of the unperturbed harmonic oscillator (still present in the spectrum because of its degeneracy), another seemingly more noteworthy spectral phenomenon occurs. There exists a threshold for the separation distance between the twin impurities, above which, the spectral curve E2(a,x0), representing the energy of the second symmetric bound state, lies entirely above Ei(a,x0), the curve pertaining to the energy for the lowest antisymmetric bound state. Below that threshold, there exists a range of values of the extension parameter, which is physically characterized by being proportional to the inverse scattering length, over which E2(a,x0) < E1(a,x0), implying the existence of two new types of level crossings. From the point of view of what could be called “spectral engineering”, the latter implies that the second symmetric bound state can become more tightly bound than the lowest antisymmetric bound state provided the separation distance is below the threshold and the extension parameter is sufficiently small. At the threshold value we have shown the existence of a single tangential intersection between the two spectral curves. As a consequence, the narrower the separation between the impurities is, the wider the aforementioned range is.
The current analysis could be further extended to the next higher eigenvalues to explore the existence of other distance thresholds leading to level crossings of the new type.
It has also been pointed out that, in full accordance with the findings of [17] for another singular double well model requiring renormalization of the coupling constant, and in total contrast to those of [11] for the one-dimensional counterpart of the model studied in this note, not requiring the renormalization procedure, the limit of our Hamiltonian as the separation distance shrinks to zero is not given by the Hamiltonian of the 3D-isotropic quantum dot with a single point impurity centered at the origin having twice the strength. It is our intention to prove in a forthcoming paper that the expected behavior in the limit can be achieved by making the coupling constant used in the renormalization procedure also dependent on the separation distance.
Finally, we also intend to extend our investigation to the two-dimensional counterpart of the model given the growing relevance of that dimension in Nanophysics, as attested by some
288
S. Albeverio, S. Fassari and F. Rinaldi
contributions to the relevant literature, such as [27-30]. Of course, as can be seen in the case of the 2D-isotropic quantum dot with a single impurity, fully studied in [31], some logarithmic divergences, which are typical of two-dimensional Quantum Mechanics, are inevitably bound to arise.
Acknowledgments
Silvestro Fassari and Fabio Rinaldi wish to express their heartfelt gratitude to Prof. I. Popov and Prof. I. Lobanov for their kind invitation to present the main contents of this note during one of the sessions of the conference “Mathematical Challenge of Quantum Transport in Nanosystems”, held at ITMO University, St. Petersburg, Russian Federation (9-11 September 2015).
Sergio Albeverio is very grateful to Prof. Nicolas Monod, Director of “Centre Interfac-ultaire Bernoulli - EPFL”, Lausanne, for the warm hospitality during his stay there as organiser jointly with Prof. A. B. Cruzeiro and Prof. D. Holm of the semester programme ’’Geometric Mechanics, Variational and Stochastic Methods” (January - June 2015).
S. Fassari and F. Rinaldi also wish to thank Prof. A. B. Cruzeiro and Prof. J. C. Zambrini for their warm hospitality on the occasion of their visit to CIB - EPFL.
References
[1] Bruning J., Geyler V., Lobanov I. Spectral properties of a short-range impurity in a quantum dot. Journal of Mathematical Physics, 2004, 45, P. 1267-1290.
[2] 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(38), 385305.
[3] 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(2), 025302.
[4] Albeverio S., Gesztesy F., H0egh-Krohn R., Holden H. Solvable models in Quantum Mechanics. AMS (Chelsea Series) second edition (with an appendix by P. Exner), 2004.
[5] Albeverio S., Kurasov P. Singular Perturbations of Differential Operators: Solvable Type Operators. Cambridge University Press, 2000.
[6] Albeverio S., Dabrowski L., Kurasov P. Symmetries of Schrodinger operators with point interactions. Letters in Mathematical Physics, 1998, 45, P. 33-47.
[7] Fassari S., Rinaldi F. On the spectrum of the Schrodinger Hamiltonian with a particular configuration of three point interactions. Reports on Mathematical Physics, 2009, 64(3), P. 367-393.
[8] Fassari S., Inglese G. On the spectrum of the harmonic oscillator with a Vtype perturbation. Helvetica Physica Acta,1994, 67, P. 650-659.
[9] Fassari S., Inglese G. Spectroscopy of a three-dimensional isotropic harmonic oscillator with a Vtype perturbation. Helvetica Physica Acta, 1996, 69, P. 130-140.
[10] Fassari S., Inglese G. On the spectrum of the harmonic oscillator with a Vtype perturbation II. Helvetica Physica Acta, 1997, 70, P. 858-865.
[11] Fassari S., Rinaldi F. On the spectrum of the Schrodinger Hamiltonian of the one-dimensional harmonic oscillator perturbed by two identical attractive point interactions. Reports on Mathematical Physics, 2012, 69(3), P. 353-370.
[12] Plancherel M., Rotach W. Sur les valeurs asymptotiques des polynomes d’Hermite. Commentarii Mathematici Helvetici, 1929, 1(1), P. 227-254.
[13] Mityagin B., Siegl P. Root system of singular perturbations of the harmonic oscillator type operators. ArXiv: 1307.6245, 2013.
[14] Reed M., Simon B. Analysis of Operators - Methods of Modern Mathematical Physics IV, Academic Press NY, 1978.
[15] Avakian M.P., Pogosyan G.S., Sissakian A.N., Ter-Antonyan V.M. Spectroscopy of a singular linear oscillator. Physics Letters A, 1987, 124( 4-5), P. 233-236.
Spectral properties of a symmetric three-dimensional quantum dot ...
289
[16] Reed M., Simon B. Fourier Analysis, Self-adjointness - Methods of Modern Mathematical Physics II, Academic Press NY, 1975.
[17] Albeverio S., Fassari S., Rinaldi F. The discrete spectrum of the spinless Salpeter Hamiltonian perturbed by ^-interactions. Journal of Physics A: Mathematical and Theoretical, 2015, 48(18), 185301.
[18] Reed M., Simon B. Functional Analysis - Methods of Modern Mathematical Physics I, Academic Press NY, 1972.
[19] Klaus M. A remark about weakly coupled one-dimensional Schrodinger operators. Helvetica Physica Acta, 1979, 52, P. 223-229.
[20] Fassari S. An estimate regarding one-dimensional point interactions. Helvetica Physica Acta, 1995, 68, P 121-125.
[21] Kato T. Perturbation Theory for Linear Operators, Springer Verlag Heidelberg-NY, 1995.
[22] Demkov Yu.N., Kurasov P.B. Von Neumann-Wigner Theorem: level repulsion and degenerate eigenvalues. Theoretical and Mathematical Physics, 2007, 153(1), P. 1407-1422.
[23] Schmidtke H.-H. Quantenchemie Weinheim: VCH, 1987 (in German).
[24] Byers Brown W., Steiner E. On the electronic energy of a one-electron diatomic molecule near the united atom. Journal of Chemical Physics, 1966, 44, P. 3934.
[25] Klaus M. On H+ for small internuclear separation. Journal of Physics A: Mathematical and General, 1983, 16, P. 2709-2720.
[26] Gonzalez-Santander C., Dominguez-Adame F. Non-local separable solutions of two interacting particles in a harmonic trap. Physics Letters A, 2011, 375, P. 314-317.
[27] Gonzalez-Santander C., Dominguez-Adame F. Exciton states and optical absorption in quantum wires under laser radiation. Physics Letters A, 2010, 374, P. 2259-2261.
[28] Gonzalez-Santander C., Dominguez-Adame F. Modelling of Coulomb interaction in parabolic quantum wires. Physica E, 2009, 41, P. 1645-1647.
[29] Lima R.P.A., Amado M., Dominguez-Adame F. A solvable model of hydrogenic impurities in quantum dots. Nanotechnology, 2008, 19, P. 135402-5406.
[30] Lopez S., Dominguez-Adame F. Non-local potential approach to the ground state of confined excitons in quantum dots. Semiconductor Science Technology, 2002, 17, P. 227-229.
[31] Tusek M. Analysis of Two-dimensional Quantum Models with Singular Potentials. Diploma Thesis Czech Technical University in Prague, Faculty of Nuclear Sciences and Physical Engineering, 2006.
