Impact Factor 2.638 | CiteScore 2.3
More on impact ›

Original Research ARTICLE

Front. Phys., 30 July 2019 |

The Birman-Schwinger Operator for a Parabolic Quantum Well in a Zero-Thickness Layer in the Presence of a Two-Dimensional Attractive Gaussian Impurity

Sergio Albeverio1,2, Silvestro Fassari1,3*, Manuel Gadella4, Luis M. Nieto4 and Fabio Rinaldi1,3
  • 1CERFIM, Locarno, Switzerland
  • 2Institut für Angewandte Mathematik, HCM, Universität Bonn, Bonn, Germany
  • 3Dipartimento di Fisica Nucleare, Subnucleare e delle Radiazioni, Università degli Studi Guglielmo Marconi, Rome, Italy
  • 4Departamento de Física Teórica, Atómica y Óptica and IMUVA, Valladolid, Spain

In this note we consider a quantum mechanical particle moving inside an infinitesimally thin layer constrained by a parabolic well in the x-direction and, moreover, in the presence of an impurity modeled by an attractive Gaussian potential. We investigate the Birman-Schwinger operator associated to a model assuming the presence of a Gaussian impurity inside the layer and prove that such an integral operator is Hilbert-Schmidt, which allows the use of the modified Fredholm determinant in order to compute the bound states created by the impurity. Furthermore, we consider the case where the Gaussian potential degenerates to a δ-potential in the x-direction and a Gaussian potential in the y-direction. We construct the corresponding self-adjoint Hamiltonian and prove that it is the limit in the norm resolvent sense of a sequence of corresponding Hamiltonians with suitably scaled Gaussian potentials. Satisfactory bounds on the ground state energies of all Hamiltonians involved are exhibited.

1. Introduction

The study of point potentials in Quantum Physics has recently received a lot of attention for a wide range of interests. First of all, point potentials serve as solvable or quasi-solvable models that approximate the action of intense and very short range potentials [13]. They have been used to model several kinds of extra thin structures [4, 5], to mimic point defects in materials, or to study heterostructures [69]. In addition, point potentials play a role in modeling impurities in quantum field theory [1013]. Furthermore, they play an important role after a recent interpretation of the Casimir effect [14, 15]. The unexpected relations between contact potentials and group theory should also be noted [16]. They also play a role in modeling Kronig-Penney crystals in condensed matter physics in various dimensions [1, 1722].

More examples of physical applications of this kind of interactions are: Bose-Einstein condensation in a harmonic trap with a tight and deep “dimple” potential, modeled by a Dirac delta function [23]; non-perturbative study of the entanglement of two directed polymers subjected to repulsive interactions given by a Dirac δ-function potential [24]; a periodic array of Dirac delta interactions is useful to investigate the light propagation in a one-dimensional realistic dielectric superlattice, which has been investigated for the transverse electric and magnetic fields and for omnidirectional polarization modes [2527].

One-dimensional quantum models with contact interactions are also useful to study a wide range of quantum properties, including scattering, since these models are quite often solvable. They also serve to acquire experience in order to analyse systems with contact potentials in higher dimensions. However, an important difference is to be pointed out: while one-dimensional contact potentials are usually defined through matching conditions at isolated points, their proper definition in higher dimensions requires a process of regularization.

Quantum two-dimensional systems are particularly interesting for their physical applications. In this particular context, the graphene deserves a special mention because of its importance, although this is not the only one two-dimensional quantum system of interest in physics. From a theoretical point of view, quantum theory in two dimensions has not yet been developed to the same extent of its one-dimensional and three-dimensional analogs, in spite of its enormous interest. Although two-dimensional quantum systems look rather simple, due to the presence of logarithmic singularities in the resolvent kernel of their free Hamiltonian, their level of complexity and, hence, their difficulty of study is far higher than the one-dimensional case (and, to a great extent, the three-dimensional case).

An important contribution to the understanding of two-dimensional quantum mechanics was provided by Duclos' article [28] on the two-dimensional hydrogen atom perturbed by a point interaction, a model that had not been dealt with in Albeverio et al. [1]. One of the main results of that paper is that the free Hamiltonian of that model, namely that of the 2D hydrogen atom, was rigorously shown to be the norm resolvent limit of the Hamiltonian of the 3D Hamiltonian of the hydrogen atom confined to an infinite planar slab of width a > 0 as a → 0+.

A remarkable feature of two-dimensional models with contact interactions, manifesting itself even in the simple case of the negative Laplacian perturbed by a point interaction, is represented by their peculiar dependence of the bound state energies on the coupling constant. As is well known, the one-dimensional model exhibits a single bound state only when the point interaction is attractive and the eigenvalue is a quadratic function of the strength of the interaction. The three-dimensional case also exhibits a single bound state only in the attractive case but the eigenvalue depends quadratically on the reciprocal of the renormalized coupling constant. In two dimensions the bound state keeps existing even if the contact interaction is repulsive and the dependence becomes exponential (see [1]). The latter behavior is confirmed even when a confinement potential is present in addition to the contact potential, which physically mimics the presence of impurities or thin barriers in the material inside which the quantum particle is moving [29, 30].

In this note, partly motivated by Duclos' paper, we intend to study a different two-dimensional model with the free Hamiltonian given by:

H0=(12d2dx2+x22)12d2dy2,    (1.1)

to which we add an attractive impurity assumed to be modeled by the isotropic Gaussian potential

W(x,y)=λV(x,y)=λe(x2+y2), λ>0,    (1.2)

so that the total Hamiltonian is

Hλ=H0+W(x,y)=H0λV(x,y)=H0λe(x2+y2),λ>0.    (1.3)

It is worth mentioning that the recent literature [3134] has shown a renewed interest in the spectral analysis of the one-dimensional Hamiltonian with a Gaussian potential, namely

H:=12d2dx2λex2/2,λ>0.    (1.4)

Therefore, (1.3) could also be regarded as a possible two-dimensional generalization of (1.4).

At this point, it is interesting to recall that a three-dimensional material with confinement in only one dimension is said to be a quantum well [35], while a 3D material with two-dimensional confinement is called a quantum wire. Therefore, in the limiting case of a quantum well inside a layer with zero thickness, it makes sense to consider the model in which the confining potential is parabolic. Due to the mathematical subtleties required, in this note we have chosen to omit the proof of the resolvent convergence of the Hamiltonian of a three-dimensional parabolic quantum well inside a thin layer to the 2D Hamiltonian (1.3) as the thickness of the layer vanishes.

Instead we start directly by writing the Green function of the two-dimensional Hamiltonian with a one-dimensional harmonic potential. Once a perturbation given by an attractive two-dimensional Gaussian potential is added, we study the properties of the corresponding Birman-Schwinger operator, that is to say the crucial part in the interaction term of the resolvent of the perturbed Hamiltonian. We remind the reader that the resolvent is the key to obtain the energy eigenvalues. We later consider the model in which the 2D Gaussian impurity potential gets replaced by one having a Dirac delta for the coordinate subjected to the harmonic confinement maintaining the Gaussian character for the other coordinate.

It may be worth pointing out that the Hamiltonian (1.1) has been used by Dell'Antonio and collaborators [36] as the free Hamiltonian in the model of a quantum system consisting of two one-dimensional particles, one of which is harmonically bound to its equilibrium position, mutually interacting by means of the contact interaction δ(xy). In other words, the interaction studied in [36] will be replaced by the one in (1.2).

Solving the eigenvalue problem for this kind of Hamiltonians is not, in general, an easy task and often requires rather sophisticated tools. One of the most widely used is the Birman-Schwinger operator, namely the integral operator

BE=(sgn W)|W|12(H0E)1|W|12,    (1.5)

and the related technique: as in most applications BE can be shown to be compact, the solutions of the eigenvalue problem for the Hamiltonian are given by those values of E for which BE has an eigenvalue equal to -1 (see [22, 37, 38] and references therein as well as [39], p. 99). Therefore, the detailed study of the properties of the Birman-Schwinger operator arising from our model is quite relevant. In the present note, we show that the Birman-Schwinger operator is Hilbert-Schmidt. We also show that Hλ is self-adjoint and bounded from below.

In addition, Hλ has a special relation with a kind of two-dimensional contact operator that will be studied in section 2.1. This is given by the Hamiltonian described heuristically by

Hλδ=H0λπδ(x)ey2,λ>0,    (1.6)

where δ(x) is the Dirac delta centered at the origin. We show that Hλδ is self-adjoint on a natural domain and can be obtained as the limit in the norm resolvent sense as n ↦ ∞ of the following sequence of Hamiltonians:

Hn,λ:=H0λn e(n2x2+y2),λ>0,    (1.7)

thus with Gaussian type potentials (as was the case for Hλ) which become increasingly more attractive and anisotropic as n goes to infinity.

Finally, it will be shown that the Hamiltonian Hλ (resp. Hλδ) is bounded from below and its lower bound can be obtained using a certain transcendental equation.

2. The Birman-Schwinger Operator for Our Model

Starting from the Hamiltonian H0 in (1.1), it is rather straightforward to infer that the associated Green function, namely the integral kernel of the resolvent operator, reads for any E<12:

(H0E)1(x,x,y,y)=n=0e2(n+12E)|yy|2(n+12E)ϕn(x)ϕn(x),    (2.1)

where ϕn(x) is the normalized n-th eigenfunction of the one-dimensional harmonic oscillator.

Once the above attractive Gaussian perturbation (1.2) is added, the total Hamiltonian is Hλ in (1.3). Therefore, its associated Birman-Schwinger integral kernel [22, 31, 32, 37, 38, 40, 41] given by (1.5) is:

BE=λB˜E(x,x1,y,y1)=λ|V|12(H0E1)|V|12(x,x1,y,y1)=      =λe(x2+y2)/2[n=0e2(n+12E)|yy1|2(n+12E)ϕn(x)ϕn(x1)]e(x12+y12)/2,    (2.2)

The main goal of this brief note is to rigorously prove that such an integral operator is Hilbert-Schmidt, that is to say tr(B~E2)<, given the evident positivity of the operator B~E (and our choice λ > 0). As the kernel of the positive operator B~E2 is clearly

B˜E2(x,x2,y,y2)=e(x2+y2)/2{[m=0e2(m+12E)|yy|2(m+12E)ϕm(x)ϕm(x)]                                   ×e(x2+y2)[n=0e2(n+12E)|yy2|2(n+12E)ϕn(x)ϕn(x2)]dxdy}e(x22+y22)/2,    (2.3)

its trace reads

tr(B˜E2)=B˜E2(x,x,y,y)dxdy              =dxdydxdye(x2+y2)e(x2+y2)                           ×m=0n=0e(2(m+12E)+2(n+12E))|yy|2(m+12E)(n+12E)                                      ϕm(x)ϕm(x)ϕn(x)ϕn(x).    (2.4)

The latter multiple integral can be rewritten as:

m=0n=0[ey2e(2(m+12E)+2(n+12E))|yy|2(m+12E)(n+12E)ey2dydy]ϕm,e(·)2ϕn2    (2.5)

where 〈f, g〉 denotes the standard scalar product of the two functions.

Let us consider the double integral in (2.5). With the notation,

f(yy):=e(2(m+12E)+2(n+12E))|yy|,  g(y)=ey2,    (2.6)

the second integral becomes the convolution (f * g)(y), so that the double integral may be written as

12(m+12E)(n+12E)ey2[(f*g)(y)]dy.    (2.7)

Using the Schwarz inequality (2.7) is smaller than or equal to

12(m+12E)(n+12E)e(·)22f*g2,    (2.8)

where ||·||p denotes the norm in Lp(ℝ). Young's inequality [42] shows that

f*grfpgq,  with  1p+1q=1r+1.    (2.9)

Therefore, with p = r = 2 and q = 1, it follows that (2.8) is smaller than or equal to

12(m+12E)(n+12E)e(·)222f1.    (2.10)

The two norms in (2.9) yield two integrals which can be easily computed, so as to obtain

π232(m+12E)(n+12E)(2(m+12E)+2(n+12E))π4(m+12E)34(n+12E)34.    (2.11)

Hence, the trace (2.4) is bounded by:

tr(B˜E2)π4m=0n=0ϕm,e(·)2ϕn2(m+12E)34(n+12E)34            =π324m=0n=0ϕmϕ0,ϕ0ϕn2(m+12E)34(n+12E)34.    (2.12)

The scalar products inside the double series can be estimated using Wang's results on integrals of products of eigenfunctions of the harmonic oscillator [43]. While the scalar product clearly vanishes if m + n = 2s + 1, when both indices are either even or odd we get:

ϕ2mϕ0,ϕ0ϕ2n2=12π[[2(m+n)]!(m+n)!]2124(m+n)(2m)!(2n)!                                            12π[[2(m+n)]!22(m+n)[(m+n)!]2]2                                            =ϕ2(m+n)4(0)2,    (2.13)
ϕ2m+1ϕ0,ϕ0ϕ2n+12=12π[[2(m+n+1)]!(m+n+1)!]2                                                                      124(m+n+1)(2m+1)!(2n+1)!                                                          12π[[2(m+n+1)]!22(m+n+1)[(m+n+1)!]2]2                                                          =ϕ2(m+n+1)4(0)2,    (2.14)

the final equalities in (2.13) and (2.14) resulting from Fassari and Inglese [44] and Mityagin and Siegl [45]. Therefore, the r.h.s. of (2.12) is bounded from above by:

tr(B˜E2)π328[m,n=0ϕ2(m+n)4(0)(2m+12E)34(2n+12E)34            +m,n=0ϕ2(m+n+1)4(0)(2m+32E)34(2n+32E)34].    (2.15)

As can be gathered from Mityagin and Siegl [45] using Stirling's formula,

                                        ϕ2n4(0)1π2n,  n1,               ϕ2(m+n)4(0)1π2(m+n),  m,n1,ϕ2(m+n+1)4(0)1π2(m+n+1),   m,n0,    (2.16)

which implies that (2.15) is bounded by

tr(B˜E2)18π12[1(12E)32+2(12E)34n=11n(2n+12E)34]                   +18π12m=1n=11m12(2m+12E)34n12(2n+12E)34                   +18π12m=0n=01(m+12)12(2m+32E)34(n+12)12(2n+32E)34             18π12[1(12E)32+2(12E)34n=11n12(2n+12E)34]                   +18π12[n=11n12(2n+12E)34]2                   +18π12[n=01(n+12)12(2n+32E)34]2<,    (2.17)

since both series involved in the final expression are clearly absolutely convergent given that the summands are positive sequences decaying like n-54. Hence, the trace of the square of the Birman-Schwinger operator, i.e., its Hilbert-Schmidt norm, is finite for any E<12.

Our result is not surprising at all since the norm could have been bounded by that of the Birman-Schwinger operator with the same impurity but with the resolvent of our H0 replaced by that of -Δ2 in two dimensions, which is known to be finite [1]. However, it provides us with a far more accurate estimate of the norm, which in turn leads to a more precise determination of the spectral lower bound resulting from the use of the Hilbert-Schmidt norm of the Birman-Schwinger operator in the KLMN theorem [46]. In fact, the latter bound is what we wish to achieve by further estimating the bottom lines of (2.17).

The series in (2.17) can be bounded from above by their respective improper integrals as follows:

S1=n=11n12(2n+12E)34=n=12(2n)12(2n+12E)34      <02 dx(2x)12(2x+12E)34     =3240s12 ds(s+12E)743240ds(s+12E)54        =32(12E)14.    (2.18)
S2=n=01(n+12)12(2n+32E)34=2(12E)34      +n=12(2n+1)12(2n+32E)34      <2(12E)34+02 dx(2x+1)12(2x+32E)34     =2(12E)342(12E)34         +3240(s+1)12 ds(s+32E)743240ds(s+32E)54     =32(32E)14.    (2.19)

Therefore, the bottom lines of (2.17) are bounded by:

tr(B˜E2)18π12[1(12E)32+6212E+18(12E)12]             +18π1218(32E)12             =18π12(12E)12[32+1(12E)12]2             +14π129(32E)12.    (2.20)

Hence, our estimate of the Hilbert-Schmidt norm of the Birman-Schwinger operator is:

tr(B˜E2)=||ex2+y22(H0E)1ex2+y22||22            18π12(12E)12[32+1(12E)12]2+14π129(32E)12.

As is well known [22, 31, 32, 37, 38, 40, 41], the operator


is isospectral to the Birman-Schwinger operator so that their Hilbert-Schmidt norms are identical. Hence, what has been achieved so far can be summarized by means of the following claim.

Theorem 2.1 The integral operators


are Hilbert-Schmidt and their Hilbert-Schmidt norms satisfy

      ||(H0E)1/2e(x2+y2)(H0E)1/2||22                         =||ex2+y22(H0E)1ex2+y22||2218π12(12E)12[32+1(12E)12]2+14π129(32E)12.    (2.21)

As an immediate consequence of the above theorem we get:

Corollary 1 The Hamiltonian


defined in the sense of quadratic forms, is self-adjoint and bounded from below by E(λ), the solution of the equation:

12(12E)12[32+1(12E)12]2+9(32E)12=4π12λ2.    (2.22)

Proof. For any E < 0 and ψQ(H0)=D(H012) (the form domain of H0):

λψ,e(x2+y2)ψ=λ(H0E)1/2ψ,[(H0E)1/2e(x2+y2)(H0E)1/2](H0E)1/2ψ    λ||(H0E)1/2e(x2+y2)(H0E)1/2||2||(H0E)1/2ψ||22   λ[18π12(12E)12[32+1(12E)12]2+14π129(32E)12]12                 [(ψ,H0ψ)E||ψ||22].    (2.23)

By taking E sufficiently negative, the first factor in the bottom line of (2.23) can be made arbitrarily small, which ensures that the Gaussian perturbation is infinitesimally small with respect to H0 in the sense of quadratic forms. Hence, we need only invoke the KLMN theorem (see [42]) to infer that Hλ is self-adjoint and bounded from below by the quantity

λ2π14[12(12E)12[32+1(12E)12]2+9(32E)12]12E,E<0,    (2.24)

so that the supremum of such lower bounds is attained for that particular value of E solving (2.22).

In the following subsections we first consider a Hamiltonian with a point interaction all along the x-direction in place of the Gaussian potential and then we investigate in detail the solution of (2.22), that is to say the lower bound of the spectrum of Hλ.

2.1. Hamiltonian With a Point Interaction Along the x-Direction

Let us consider now the Hamiltonian


that is to say the energy operator given by the same H0 as before but with the interaction term having a point interaction in place of the Gaussian along the x-direction. Our goal is to prove that such an operator is self-adjoint and that it is the limit in the norm resolvent sense of the sequence H0 − λVn(x, y) as n → ∞, with Vn(x,y)=ne-(n2x2+y2). As is to be expected, our approximating sequence is quite different from the one used in Albeverio et al. [1] to get the Laplacian perturbed by a point interaction in two dimensions. Before stating and proving the main result of this section, we wish to provide the reader with the visualization of the approximating potentials in Figure 1.


Figure 1. Approximating potentials Vn(x,y)=ne-(n2x2+y2) for n = 15 (pink), 30 (magenta), and 70 (cyan).

Corollary 2 The Hamiltonian Hλδ=H0-λπδ(x)e-y2, defined in the sense of quadratic forms, is self-adjoint and is the norm resolvent limit of the sequence of Hamiltonians


Furthermore, Hλδ is bounded from below by Eδ(λ), the solution of the equation:

1(12E)12[1(12E)12+32]2=4π12λ2.    (2.25)

Proof. First of all, it is quite straightforward to show that the integral operator


is Hilbert-Schmidt with the square of the Hilbert-Schmidt norm given by:


As the double integral involving the convolution has already been estimated in (2.11), the latter expression is bounded by:

π3/24m=0n=0ϕ2m2(0)ϕ2n2(0)(2m+12E)34(2n+12E)34=π3/24[n=0ϕ2n2(0)(2n+12E)34]2,    (2.26)

which, using (2.17), is bounded by:

14π1/2[1(12E)34+n=11n12(2n+12E)34]2          14π1/2[1(12E)34+32(12E)14]2                                     =14π12(12E)12[1(12E)12+32]2,    (2.27)

having taken advantage of (2.18). As the right hand side of (2.27) can be made arbitrarily small by taking E < 0 large in absolute value, the KLMN theorem ensures, as was done previously in the case of Hλ, the self-adjointness of Hλδ as well as the existence of the spectral lower bound Eδ(λ) given by the solution of (2.25).

As to the convergence of Hn, λ to Hλδ, we start by noting that, for any E < 0, the operator (H0-E)-1/2ne-(n2x2+y2)(H0-E)-1/2 converges weakly to (H0-E)-1/2πδ(x)e-y2(H0-E)-1/2 as n → ∞. Furthermore,

||(H0E)1/2ne(n2x2+y2)(H0E)1/2||22=||nen2x2+y22(H0E)1en2x2+y22||22=l,m=0[ey2e(2(l+12E)+2(m+12E))|yy|2(l+12E)(m+12E)ey2dydy]        ϕl,nen2(·)2ϕm2​​.    (2.28)


ϕl,nen2(·)2ϕm=nen2x2ϕl(x)ϕm(x)dx                                       =ex2ϕl(x/n)ϕm(x/n)dxπϕl(0)ϕm(0),

as n → ∞, the right hand side of (2.28) converges to

π2l=0m=0ϕ2l2(0)ϕ2m2(0)    [ey2e(2(2l+12E)+2(2m+12E))|yy|(2l+12E)(2m+12E)ey2dydy]    =π||(δ(x)ey2)12(H0E)1(δ(x)ey2)12||22   =π||(H0E)1/2δ(x)ey2(H0E)1/2||22.

Hence, the Hilbert-Schmidt norm of (H0-E)-1/2ne-(n2x2+y2)(H0-E)-1/2 converges to the Hilbert-Schmidt norm of (H0-E)-1/2πδ(x)e-y2(H0-E)-1/2 as n → ∞. Due to Theorem 2.21 in Simon [47], this fact and the previous weak convergence imply that the convergence actually takes place in the Hilbert-Schmidt norm. Then, the norm convergence of these integral operators ensures the norm resolvent convergence of Hn, λ to Hλδ, as guaranteed by Theorem VIII.25 in Reed and Simon [46], which completes our proof of Corollary 2.

2.2. The Lower Bound of σ(H0-λe-(x2+y2))

As anticipated earlier, the lower bound of the spectrum of Hλ in (1.3) is the function E(λ) given implicitly by the equation (2.22). From this expression, some approximate results can be easily obtained in two different regimes. For example, it is possible to prove that the asymptotic behavior of (2.22) for large values of both variables E(λ), λ is

E(λ)=814π2λ4.    (2.29)

On the other hand, for small values of λ we can prove that (2.22) behaves approximately as follows

E(λ)=1214π2/3λ4/3.    (2.30)

A plot of the λ-dependence (λ being the strength of the potential of Hλ) of the lower bound of the energy E(λ), resulting from the solution of (2.22), as well as those of the two approximations given by (2.29) and (2.30), are given in Figure 2.


Figure 2. A plot of the lower bound of the energy E as a function of λ resulting from the solution of Equation (2.22) (blue curve), the approximate expression valid for large values of E and λ obtained in (2.29) (yellow curve) and the approximation for small values of λ as in (2.30) (green curve). In the inset we have enlarged the region where λ and E are small. While the similarity between the solution of (2.22) and the funtion (2.29) is quite acceptable for a wide range of the parameters, the solution of (2.22) is satisfactorily approximated by (2.30) only for very small values of λ.

3. Final Remarks

In this note we have analyzed in detail the Birman-Schwinger operator of the two-dimensional Hamiltonian Hλ=H0-λe-(x2+y2), namely the integral operator -λe-x2+y22(H0-E)-1e-x2+y22 where H0=(-12d2dx2+x22)-12d2dy2. In particular, we have rigorously shown that the operator is Hilbert-Schmidt and have estimated its Hilbert-Schmidt norm. This fact has enabled us to use the KLMN theorem to determine a lower bound for the spectrum of Hλ, that is to say E(λ), the implicit function representing the solution of an equation involving the energy parameter and the coupling constant. Furthermore, we have investigated the Hamiltonian Hλδ=H0-λπδ(x)e-y2, having the Gaussian impurity in the direction subjected to the harmonic confinement replaced by a point impurity.

As anticipated in the introduction, the proof of the resolvent convergence, as the thickness of the layer vanishes, of the Hamiltonian of a three-dimensional parabolic quantum well inside a thin layer to the 2D Hamiltonian (1.3) has been put off as it may deserve a separate paper.

The results of this article will enable us to study the lowest bound states created by the Gaussian impurity potential of the aforementioned Hamiltonian by means of the modified Fredholm determinant det2[1-λex2+y22(H0-E)-1e-x2+y22], the regularized determinant used to handle Hilbert-Schmidt operators. Work in this direction is in progress.

Data Availability

All datasets generated for this study are included in the manuscript and/or the supplementary material.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.


This work was partially supported by the Spanish MINECO (MTM2014-57129-C2-1-P), Junta de Castilla y León and FEDER projects (BU229P18, VA057U16, and VA137G18).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The authors are indebted to both referees for their suggestions and corrections resulting in the significant improvement of the paper.


1. Albeverio S, Gesztesy F, Høegh-Krohn R, Holden H. Solvable Models in Quantum Mechanics, 2nd ed. Providence, RI: AMS Chelsea Publishing (2004). doi: 10.1090/chel/350

CrossRef Full Text | Google Scholar

2. Seba P. Some remarks on the δ-interaction in one dimension Czech. J Phys. (1986) 36:667–73. doi: 10.1007/BF01597402

CrossRef Full Text | Google Scholar

3. Bräunlich G, Hainzl C, Seiringer R. On contact interactions as limits of short-range potentials Met. Funct Anal Top. (2013) 19:364–75.

Google Scholar

4. Zolotaryuk AV, Zolotaryuk Y A zero-thickness limit of multilayer structures: a resonant-tunnelling δ′-potential. J Phys A. (2015) 48:035302. doi: 10.1088/1751-8113/48/3/035302

CrossRef Full Text | Google Scholar

5. Zolotaryuk AV, Zolotaryuk Y. Controllable resonant tunnelling through single-point potentials: a point triode. Phys Lett A. (2015) 379:511–7. doi: 10.1016/j.physleta.2014.12.016

CrossRef Full Text | Google Scholar

6. Gadella M, Heras FJH, Negro J, Nieto LM. A delta well with a mass jump. J Phys A. (2009) 42:465207. doi: 10.1088/1751-8113/42/46/465207

CrossRef Full Text | Google Scholar

7. Kulinskii VL, Panchenko DY. Physical structure of point-like interactions for one-dimensional Schrödinger and the gauge symmetry. Phys B Cond Matt. (2015) 472:78–83. doi: 10.1016/j.physb.2015.05.011

CrossRef Full Text | Google Scholar

8. Kulinskii VL, Panchenko DY. Mass-jump and mass-bump boundary conditions for singular self-adjoint extensions of the Schrödinger operator in one dimension. Ann Phys. (2019) 404:47–56. doi: 10.1016/j.aop.2019.03.001

CrossRef Full Text | Google Scholar

9. Kulinskii VL, Panchenko DY. Point-like rashba interactions as singular self-adjoint extensions of the Schrödinger operator in one dimension. Front Phys. (2019) 7:44. doi: 10.3389/fphy.2019.00044

CrossRef Full Text | Google Scholar

10. Muñoz-Castañeda JM, Mateos-Guilarte J, Mosquera AM. Quantum vacuum energies and Casimir forces between partially transparent δ-function plates. Phys Rev D (2013) 87:105020. doi: 10.1103/PhysRevD.87.105020

CrossRef Full Text | Google Scholar

11. Asorey M, García-Alvarez D, Muñoz-Castañeda JM. Casimir effect and global theory of boundary conditions. J Phys A Math Gen. (2006) 39:6127–36. doi: 10.1088/0305-4470/39/21/S03

CrossRef Full Text | Google Scholar

12. Muñoz-Castañeda JM, Kirsten K, Bordag M. QFT over the finite line. Heat kernel coefficients, spectral zeta functions and selfadjoint extensions. Lett Math Phys. (2015) 105:523–49. doi: 10.1007/s11005-015-0750-5

CrossRef Full Text | Google Scholar

13. Muñoz-Castañeda JM, Mateos-Guilarte J. δ − δ′ generalized Robin boundary conditions and quantum vacuum fluctuations. J Phys Rev D. (2015) 91:025028. doi: 10.1103/PhysRevD.91.025028

CrossRef Full Text | Google Scholar

14. Bordag M. Monoatomically thin polarizable sheets. Phys Rev D. (2014) 89:125015. doi: 10.1103/PhysRevD.89.125015

CrossRef Full Text | Google Scholar

15. Asorey M, Muñoz-Castañeda JM. Attractive and repulsive Casimir vacuum energy with general boundary conditions. Nucl Phys B. (2013) 874:852–76. doi: 10.1016/j.nuclphysb.2013.06.014

CrossRef Full Text | Google Scholar

16. Gadella M, Mateos-Guilarte J, Muñoz-Castañeda JM, Nieto LM. Two-point one-dimensional δ − δ′ interactions: non-abelian addition law and decoupling limit. J Phys A Math. Theor. (2016) 49:015204. doi: 10.1088/1751-8113/49/1/015204

CrossRef Full Text | Google Scholar

17. Fairbairn WM, Glasser ML, Steslicka M. Relativistic theory of surface states. Surf Sci. (1973) 36:462–71. doi: 10.1016/0039-6028(73)90394-4

CrossRef Full Text | Google Scholar

18. Fassari S. On the Schrödinger operator with periodic point interactions in the three-dimensional case. J Math Phys. (1984) 25:2910–7. doi: 10.1063/1.526038

CrossRef Full Text | Google Scholar

19. Fassari S. Spectral properties of the Kronig-Penney Hamiltonian with a localized impurity. J Math Phys. (1989) 30:1385–92. doi: 10.1063/1.528320

CrossRef Full Text | Google Scholar

20. Fassari S. On the bound states of non-relativistic Kronig-Penney Hamiltonians with short range impurities. Helv Phys Acta. (1990) 63:849–83.

Google Scholar

21. Fassari S. On the bound states of relativistic Kronig-Penney Hamiltonians with short range impurities. Helv Phys Acta. (1990) 63:884–907.

Google Scholar

22. Fassari S, Klaus M. Coupling constant thresholds of perturbed periodic Hamiltonians. J Math Phys. (1998) 39:4369–16. doi: 10.1063/1.532516

CrossRef Full Text | Google Scholar

23. Uncu H, Tarhan D, Demiralp E, Müstecaplıoğlu ÖE. Bose-Einstein condensate in a harmonic trap decorated with Dirac δ-functions. Phys Rev A. (2007) 76:013618. doi: 10.1103/PhysRevA.76.013618

CrossRef Full Text | Google Scholar

24. Ferrari F, Rostiashvilii VG, Vilgis TA. Directed polymers with constrained winding angle. Phys Rev E. (2005) 71:061802. doi: 10.1103/PhysRevE.71.061802

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Alvarado-Rodríguez I, Halevi P, Sánchez-Mondragón JJ. Density of states for a dielectric superlattice: TE polarization. Phys Rev E (1999) 59:3624–30. doi: 10.1103/PhysRevE.59.3624

CrossRef Full Text | Google Scholar

26. Zurita-Sánchez JR, Halevi P. Density of states for a dielectric superlattice. II. TM polarization. Phys Rev E. (2000) 61:5802–7. doi: 10.1103/PhysRevE.61.5802

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lin MC, Jao RF. Quantitative analysis of photon density of states for a realistic superlattice with omnidirectional light propagation. Phys Rev E. (2006) 74:046613. doi: 10.1103/PhysRevE.74.046613

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Duclos P, Štovícek P, Tušek M. On the two-dimensional Coulomb-like potential with a central point interaction. J Phys A Math Theor. (2010) 43:474020. doi: 10.1088/1751-8113/43/47/474020

CrossRef Full Text | Google Scholar

29. Fassari S, Gadella M, Glasser ML, Nieto LM, Rinaldi F. Level crossings of eigenvalues of the Schrödinger Hamiltonian of the isotropic harmonic oscillator perturbed by a central point interaction in different dimensions. Nanosyst Phys Chem Math. (2018) 9:179–86. doi: 10.17586/2220-8054-2018-9-2-179-186

CrossRef Full Text | Google Scholar

30. Fassari S, Gadella M, Glasser ML, Nieto LM, Rinaldi F. Spectral properties of the 2D Schrödinger Hamiltonian with various solvable confinements in the presence of a central point perturbation. Phys Scr. (2019) 94:055202. doi: 10.1088/1402-4896/ab0589

CrossRef Full Text | Google Scholar

31. Muchatibaya G, Fassari S, Rinaldi, F, Mushanyu J. A note on the discrete spectrum of Gaussian wells (I): the ground state energy in one dimension. Adv Math Phys. (2016) 2016:2125769. doi: 10.1155/2016/2125769

CrossRef Full Text | Google Scholar

32. Fassari S, Gadella M, Nieto L M, Rinaldi F. On the spectrum of the 1D Schrödinger Hamiltonian perturbed by an attractive Gaussian potential. Acta Polytech. (2017) 57:385–90. doi: 10.14311/AP.2017.57.0385

CrossRef Full Text | Google Scholar

33. Fernández FM. Quantum Gaussian wells and barriers. Am J Phys. (2011) 79:752–4. doi: 10.1119/1.3574505

CrossRef Full Text | Google Scholar

34. Nandi S. The quantum Gaussian well. Am J Phys. (2010) 78:1341–5. doi: 10.1119/1.3474665

CrossRef Full Text | Google Scholar

35. Harrison P. Quantum wells, wires, dots. In: Theoretical and Computational Physics of Semiconductor Nanostructures. 3rd ed. Chichester: Wiley (2009).

Google Scholar

36. Correggi M, Dell'Antonio G, Finco D. Spectral analysis of a two-body problem with zero-range perturbation. J Funct Anal. (2008) 255:502–531. doi: 10.1016/j.jfa.2008.04.005

CrossRef Full Text | Google Scholar

37. Klaus M. A remark about weakly coupled one-dimensional Schrödinger operators. Helv Phys Acta. (1979) 52:223.

Google Scholar

38. Klaus M. Some applications of the Birman-Schwinger principle. Helv Phys Acta. 55, 49–68.

Google Scholar

39. Reed M, Simon B. Analysis of Operators. Methods in modern mathematical physics. New York, NY: Academic Press (1978).

Google Scholar

40. Fassari S. An estimate regarding one-dimensional point interactions. Helv Phys Acta. (1995) 68:121–5.

Google Scholar

41. Fassari S, Rinaldi F. On the spectrum of the Schrödinger Hamiltonian with a particular configuration of three point interactions. Rep Math Phys. (2009) 64:367–93. doi: 10.1016/S0034-4877(10)00004-2

CrossRef Full Text | Google Scholar

42. Reed M, Simon B. Fourier Analysis, Self-adjointness. Methods in modern mathematical physics. New York, NY: Academic Press (1975).

Google Scholar

43. Wang WM. Pure point spectrum of the floquet Hamiltonian for the quantum harmonic oscillator under time quasi-periodic perturbations commun. Math Phys. (2008) 277:459–96. doi: 10.1007/s00220-007-0379-z

CrossRef Full Text | Google Scholar

44. Fassari S, Inglese G. On the spectrum of the harmonic oscillator with a δ-type perturbation. Helv Phys Acta. (1994) 67:650–9.

Google Scholar

45. Mityagin BS, Siegl P. Root system of singular perturbations of the harmonic oscillator type operators. Lett Math Phys. (2016) 106:147–67. doi: 10.1007/s11005-015-0805-7

CrossRef Full Text | Google Scholar

46. Reed M, Simon B. Functional Analysis. Methods in Modern Mathematical Physics. New York, NY: Academic Press (1972).

Google Scholar

47. Simon B. Trace Ideals and Their Applications. Cambridge: Cambridge University Press (1979).

Google Scholar

Keywords: Gaussian potential, Birman-Schwinger operator, Hilbert-Schmidt operator, contact interaction, quantum well

Citation: Albeverio S, Fassari S, Gadella M, Nieto LM and Rinaldi F (2019) The Birman-Schwinger Operator for a Parabolic Quantum Well in a Zero-Thickness Layer in the Presence of a Two-Dimensional Attractive Gaussian Impurity. Front. Phys. 7:102. doi: 10.3389/fphy.2019.00102

Received: 30 March 2019; Accepted: 02 July 2019;
Published: 30 July 2019.

Edited by:

Manuel Asorey, University of Zaragoza, Spain

Reviewed by:

Yuriy Golovaty, Lviv University, Ukraine
Q. H. Liu, Hunan University, China

Copyright © 2019 Albeverio, Fassari, Gadella, Nieto and Rinaldi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Silvestro Fassari,