arxiv:math-ph/ v1 9 Mar 2006

Transcription

1 Effective models for excitons in carbon nanotubes nd of February, 6 arxiv:math-ph/636v 9 Mar 6 Horia D. Cornean, Pierre Duclos, Benjamin Ricaud 3 Abstract We analyse the low lying spectrum of a model of excitons in carbon nanotubes. Consider two particles with a Coulomb self-interaction, placed on an infinitely long cylinder. If the cylinder radius becomes small, the low lying spectrum is well described by a one-dimensional effective Hamiltonian which is exactly solvable. Contents Introduction The mathematical model 3. Center of mass separation in the longitudinal direction The self-adjointness problem An effective operator for the low lying spectrum Spectral analysis of H r eff 3. A solvable comparison operator H C The eigenvalues and eigenvectors of H C Approximation of H r eff by H C Reduction of H r to Heff r 6 4. V offdiag is H diag -form bounded Comparison between H r and H diag The main theorem and some conclusions 3 Dept. Math., Aalborg University, Fredrik Bajers Vej 7G, 9 Aalborg, Denmark; cornean@math.aau.dk Centre de Physique Théorique UMR 67 - Unité Mixte de Recherche du CNRS et des Universités Aix-Marseille I, Aix-Marseille II et de l université du Sud Toulon-Var - Laboratoire affilié à la FRUMAM, Luminy Case 97, F-388 Marseille Cedex 9 France; duclos@univtln.fr 3 Centre de Physique Théorique UMR 67 - Unité Mixte de Recherche du CNRS et des Universités Aix-Marseille I, Aix-Marseille II et de l université du Sud Toulon-Var - Laboratoire affilié à la FRUMAM, Luminy Case 97, F-388 Marseille Cedex 9 France; ricaud@cpt.univmrs.fr

2 Introduction In order to understand the quantum mechanics of carbon nanotubes, one has to reconsider many classical problems in which the systems are restricted to low dimensional configuration spaces. The effects induced by these special shapes are significant. For example, optical properties and electrical conductivity in nanotubes and nanowires are highly influenced by their geometry. In a periodic structure, bands of allowed and forbidden energies are characteristic for non-interacting electrons. When applying an external perturbation, such as light, electrons can only absorb the amount of energy which allows them to jump from an occupied level of energy to a free one. In the particular case of a semi-conductor, at low temperatures, the band of energies are either full (valence bands) or empty (conduction bands). So the electron must absorb a fairly large amount of energy to jump to the conduction band. When the self-interaction is also considered, the mathematical problem of the optical response becomes very difficult, and there are not many rigorous results in this direction. Here is a sketchy description of what physicists generally do (see for example the book of Fetter and Walecka [FW]:. Work in the grand-canonical ensemble, at zero temperature, and the Fermi energy E F is in the middle of an energy gap;. Switch to an electron-hole representation, via a Bogoliubov unitary transformation. The new non-interacting ground state is the tensor product of two vacua. If before an excited state meant that an electron was promoted from an occupied energy state from below E F to an empty state above E F, in the electron-hole representation it simply means that an electron-hole pair was created; 3. Try to diagonalize in a way or the other the true, self-interacting many-body Hamiltonian by restricting it to a certain subspace of physically relevant excited states; this usually amounts to formulate a Hartree-Fock problem; 4. Try to obtain an effective one-body Hamiltonian, whose spectrum approximates in some way the original problem in a neighborhood of E F ; 5. The bound states of this effective one-body operator are called excitons. They describe virtual, not real particles; 6. Use the exciton energies to correct the optical response formulas derived in the non-interacting case. Now this one-body effective Hamiltonian is a complicated object in general. If one makes a number of further assumptions like:. There is only one conduction band above E F and only one valence band below E F ;

3 . The dispersion law of these two energy bands is replaced with paraboloids (effective mass approximation), then this one body effective Hamiltonian is nothing but the one describing the relative motion of a positively charged particle (a hole) and a negatively charged particle (electron), interacting through an eventually screened Coulomb potential. Let us stress that this procedure is generally accepted as physically sound in the case when the crystal is periodic in all three dimensions. If a special geometry is imposed (i.e. the electrons are confined on long and thin cylinders made out of carbon atoms) then the above procedure has to be completely reconsidered. The problem is even more complicated, because two dimensions are on a torus and the band structure only arises from the longitudinal variable. Even the position of the Fermi level moves when the radius of the cylinder varies. It has been argued by physicists [Ped, Ped] that one can still write down a Hartree-Fock type eigenvalue problem which describes the excited states near the Fermi level. This operator is a two-body one, which does not in general allow a complete reduction of its mass center. A mathematically sound formulation of this Hartree-Fock problem would be of certain interest, but it is not what we do in this paper. We rather perform the spectral analysis of an operator which has been conjectured by physicists as being the relevant one. The main point in investigating these low dimensional structures, is that the interaction between electrons is enhanced and gives much stronger exciton effects than in three dimensions. This means that some new energy states appear deep inside the forbidden energy band. The smaller the radius, the more important these new energies are. That is because they allow photons with much smaller energy that the band gap to be absorbed into the material. We will therefore consider two charged particles living on a cylinder and interacting through a Coulomb potential. As we have already pointed out above, this operator models an effective Hamiltonian for excitons in carbon nanotubes, according to [Ped] and [Ped]. Remember that it has nothing to do with real particles living on the nanotube, the exciton being just a mathematical artifact describing virtual particles. We also hope that our current results could also describe excitons living in nanowires [Aki], or two particles in a strong magnetic field as in [BD]. Let us mention that our paper is an improvement and a continuation of a previous work done in [CDP]. The mathematical model Our configuration space is a cylinder of radius r and infinite length, space denoted by C r = R rs, S being the unit circle. The coordinates on the cylinder are (x, y) (R rs ) where x is the variable along the tube axis and y is the transverse coordinate. 3

4 The two virtual particles live in the Hilbert space L (C r C r ). We formally consider the Hamiltonian ( H r = x + x + y + ) y V r (x x, y y ), (.) m m m m where e e V r (x, y) := ε x + 4r sin ( ) (.) y r (x i, y i ) are the coordinates on the cylinder of the two charged particles, m i their masses, and e i their charges. Here ε is the electric permittivity of the material. The potential V r is the three dimensional Coulomb potential simply restricted to the cylinder. We justify the expression of V r by Pythagora s theorem. The cylinder is embedded in R 3. The distance ρ from one particle to the other in R 3 is: ( ) ρ = (x x ) + 4r sin y y r where r sin ( y y ) r is the length of the chord joining two points of coordinate y and y on the circle. Now consider the space D = { f C (C r C r ) : α,, β,, γ, N (.3) } x α xα Dβ x Dx β Dy γ Dy γ f(x, y, x, y ) C α,β,γ of Schwartz functions along x and smooth and πr-periodic along y. Clearly D is dense in the Sobolev space H (Cr ). Let us define on D the kinetic quadratic form t [ψ] = ( x ψ + y ψ ) + ( x ψ + y ψ ) (.4) m m and the quadratic form associated with the Coulomb potential on the cylinder: t V [ψ] = V r ψ, (.5) { D ψ L (Cr ), } V r (x x, y y ) ψ(x, x, y, y ) dx dx dy dx <. C r Finally, define the form t H := t t V on the domain D. (.6) The sesquilinear form induced by t is densely defined, closable, symmetric, nonnegative, and its corresponding selfadjoint operator H is with periodic boundary conditions in the y variables. Its form domain is H (Cr ), and is essentially self-adjoint on D. 4

5 . Center of mass separation in the longitudinal direction We introduce the total mass M := m + m and the reduced mass µ := mm m +m. Denote as usual with { X = m x +m x m +m, x = x x, Y = y, y = y y then, { x = X + m M x, x = X m M x, y = y + Y, y = Y Unfortunately, for the y, y variables we cannot use Jacobi coordinates because the transformation does not leave invariant the domain of the Laplacian, because of the periodic boundary conditions. That is why we use atomic coordinates y and Y instead. In these new coordinates, the total Hilbert space splits in a tensor product L (R) L [R (rs ) ]. More precisely, if we denote by U : L (R) L [R (rs ) ] L (C r ), [U f](x, y, x, y ) = f(x, x, Y, y) then it is quite standard to show that after this variable change we can separate away the X variable and for f, g U[D ] we get: t H (U f, U g) = M Xf, X g + m Y f, Y g (.7) + µ xf, x g + µ yf, y g y f, Y g V r (x, y)f, V r (x, y)g. m Note that the subset U[D ] has the same properties as D but in the new variables. Therefore we can concentrate on the reduced form t h (f, g) = m Y f, Y g (.8) + µ xf, x g + µ yf, y g y f, Y g V r (x, y)f, V r (x, y)g, m densely defined on smooth enough functions in L [R (rs ) ], decaying along the x variable. Consider the decomposition L [R (rs ) ] = k Z L (R rs ) implemented by the Fourier series f(x, y, Y ) = k Z e iy k r πr ˆfk,r (x, y). 5

6 where ˆf k,r (x, y) = πr f(x, y, Y )e iy k r dy πr Then for our form t h we get: t h =: t hk, (.9) k Z where t hk is defined on a domain E k, t hk (f, g) = µ xf, x g + µ yf, y g + k f, g (.) m r ik m r f, yg V (x, y)f, V (x, y)g. E k = { f C (R rt) x α D β x Dγ y f(x, y) C αβγ, α, β, γ N }. Now remember that we are only interested in the low lying spectrum of our original operator. We will now show that for small r, only t h contributes to the bottom of the spectrum. Indeed, let us concentrate on the operator µ y + k m r i k m r y defined on rs with periodic boundary conditions. Via the discrete Fourier transformation it is unitarily equivalent to: ( p µ r + k m r kp ) m r (.) p Z where p µ r + k m r kp m r = ( µ m (p, k) r m m )( p k). (.) A simple calculation shows that both eigenvalues of the above matrix are positive; denote with λ the smaller one. Then the operator in (.) obeys: p Z ( p µ r + k m r Using this in (.) we obtain the inequality kp ) m r λ y + λ k r. (.3) 6

7 t hk (f, f) min {, µλ } µ [ xf, x f + y f, y f)] V r (x, y)f, V r (x, y)f + λ k r f = t h (f, f) + λ k r f (.4) where t h is obviously defined by the previous line. Now one of the results obtained in this paper will be that the specrum of the self-adjoint operator associated to a form like t h is bounded from below by a numerical constant times (ln(r)). Hence if k and r is small enough, all t hk will be positive and only t h will contribute to the negative part of the spectrum.. The self-adjointness problem Due to (.7), (.8), (.9) and (.4), it is clear that it is enough to concentrate on t h. If we can prove that it is bounded from below, then all other forms with k will also have this property, and the total Hamiltonian will be a direct sum of Friederichs extensions. Because we can anyway scale the masses and charges away, and in order to simplify the notation, let us consider the sesquilinear form: on t H (f, g) := [ xf, x g + y f, y g ] V r (x, y)f, V r (x, y)g =: t (f, g) t V (f, g) (.5) E = { f C (R rt) x α D β x Dγ y f(x, y) C αβ, α, β, γ N }, now where V r is as in (.) but with e e =. We will now construct a self-adjoint operator out of this form. Proposition.. The form t V is a relatively compact perturbation of t, with relative bound. Thus the form t H is bounded from below, closable, and the closure defines a self adjoint operator H r whose form domain is the Sobolev space H (C r ). Moreover, σ ess ( H r ) = [, ). Proof. We identify the cylinder C r with R [ rπ, rπ]. For every α > we define H α (C r ) to be the set of all functions which (at least formally) can be expressed as: f(x, y) = π e ikx+imy/r ˆfm (p) dp, r m Z R ( + p α + m α ) ˆf m (p) dp <. (.6) m Z R 7

8 Let χ be the characteristic function of the interval [ r/, r/]. Since near the boundary of the strip V r (x, y) ( χ(y)) is bounded, we only have to look at Ṽ (x, y) := V (x, y) χ(y). Then we can find a constant C such that everywhere in C r to have Ṽ (x, y) C x + y. Denote with ρ := x + y. Choose a function χ C (R) with support in ( 3r/, 3r/), such that χχ = χ. Then the operator of multiplication by χ is bounded from H α (C r ) to H α (R ) and vice versa, because it does not touch the boundary (the proof of this fact is more or less standard). Moreover, if is the operator associated to t, then we have H α (C r ) = ( + ) α/ L (C r ). Note that E is dense in any H α (C r ). Moreover, for every ψ E we have ψ, Ṽ ψ L (C r) C χ ψ, ρ χ ψ L (R ). (.7) We have that χ ψ S(R ). Then we can write χ ψ, (/ρ)χ ψ L (R ) = π χ (ρ sin(θ))ψ(ρ, θ) ψ(ρ, θ)dρdθ and after integration by parts in the radial integral we obtain χ ψ, (/ρ)χ ψ L (R ) = π { ρ [χ (ρ sin(θ))ψ(ρ, θ) χ (ρ sin(θ))ψ(ρ, θ)]}ρdρdθ. (.8) Then using the estimate ρ (χ ψ) (χ ψ), and with the Cauchy-Schwarz inequality: χ ψ, (/ρ)χ ψ L (R ) C χ ψ L (R ) (χ ψ) L (R ) const ψ L (C r) ψ H (C r). (.9) Now for an arbitrarily small ǫ > we have ψ, Ṽ ψ L (C r) (C /ǫ) ψ L (C r) + C ǫ t (ψ, ψ), where C is just a numerical constant. The density of E in H (C r ) finishes the proof of relative boundedness, and we can define H r as the Frederichs extension. 8

9 Until now we have shown in an ellementary way that V r ( + ) / is bounded, but one can do much better than that. In [Bo] it has been proven a two dimensional version of an inequality of Kato, which states the following: ψ, x ψ L (R ) Γ(/4)4 4π ψ, ψ L (R ). (.) This inequality immediately implies that V r : H /4 (C r ) L (C r ) is bounded. Now let us show that the operator V r H / (C r ) L (C r ) is compact. We will in fact prove the sufficient condition that the operator T := x / ( + ) / defined on L (R ) is compact. Indeed, let us denote by χ n the characteristic function of the ball of radius n >, centred at the origin in R. Then we can write: T = χ n (x)tχ n ( ) + [( χ n )(x)]t + χ n (x)t[( χ n )( )]. First, the operator χ n (x)tχ n ( ) is Hilbert-Schmidt (its integral kernel is an L (R 4 ) function), thus compact. Second, the sequence of operators [( χ n )(x)]t converges in norm to zero. Third, the sequence χ n (x)t[( χ n )( )] can be expressed in the following way: χ n (x)t[ χ n ( )] = {χ n (x) x / ( + ) /4 }( + ) /4 [( χ n )( )], where the first factor is uniformly bounded in n, while the second one converges in norm to zero. Thus T can be approximated in operator norm with a sequence of compact operators, hence it is compact. Therefore V r is a relatively compact form perturbation to, hence the essential spectrum is stable, and the proof is over..3 An effective operator for the low lying spectrum We will show in this section that at small r, the negative spectrum of H r can be determined by studying a one dimensional effective operator Heff r. It is natural to expect that the high transverse modes do not contribute much to the low region of the spectrum. First, we separate H r into different parts taking advantage of the cylindrical geometry, that is to say, we represent H r as a sum of orthogonal transverse modes using the periodic boundary conditions along the circumference of the cylinder. Second, we analyse which part is relevant when the radius tends to zero. We recall that H r is formally given by H r = x y V r in the space L (C r ) L (R) L (rs ). The domain contains all ψ H (C r ) with the property that in distribution sense we have ( x ) y V r ψ L (C r ). (.) 9

10 This does not mean that the domain is H (C r ) because V r is too singular at the origin. Our problem has two degrees of freedom. We consider the orthonormal basis of eigenvectors of y with domain H per((, πr)) H (rs ). Here, the Sobolev space Hper(( πr, πr)) denote functions which are πr-periodic with first and second derivatives in the distribution sense in L. We can write y = n= E r nπ r n where the one dimensional projectors Π r n are defined by Π r n = χ r n, χ r n, χ r n(y) = πr e in y r We now introduce a family of orthogonal projectors and E r n = n r, n Z. Π r n := Π r n, (.) which project from L (C r ) into what we call the n th transverse mode. The operator H r can be split as follows: H r = n,m Π r n H r Π r m =: n H r n ( χ n χ n ) + n m H r n,m ( χ n χ m ), (.3) where the sum is a direct sum, since the projectors are orthogonal. By a natural unitary identification, we can work in a new Hilbert space: H = l [Z; L (R)], H ψ = {ψ n } n Z, ψ n L (R). (.4) Therefore our original operator is an infinite matrix now, {H n,m } n,m Z whose elements are operators in L (R). If n m, the only contribution comes from V r, and the corresponding operator is a multiplication operator given by (x ): V r n,m(x) := πr πr πr V (x, y y i(m n) )e r dy, x. (.5) If n = m, then the corresponding diagonal element is given by the operator: H r n = x V r eff + n r, (.6) where Veff r is deduced from V n,m r when m = n and is given by V r eff (x) = πr πr πr V (x, y)dy. (.7)

11 Finally, let us introduce a last notation for what will be our effective one-dimensional comparison operator: and note that d Heff r := dx V eff(x) r (.8) Hn r = Heff r + n r. (.9) One can see that for n, the diagonal entries of our infinite operator valued matrix are pushed up by a term proportional with /r. Thus a natural candidate for a comparison operator for the negative spectrum of H r is Heff r. In the next section we will perform a careful study of this operator. 3 Spectral analysis of H r eff We now want to study the spectrum of the operator Heff r when r becomes small. We recall that: where d Heff r = dx V eff(x) r Veff r (x) = πr πr πr x + 4r sin y r with form domain Q(Heff r ) = H (R). We are going to use perturbation theory around r =, which will turn out to be quite a singular limit. The strategy is to approximate the form associated to the potential Veff r (x) around r = by another quadratic form which provides a solvable approximation. Let us define the sesquilinear form on S(R) (later on we will show that it is bounded on H (R)): C (f, g) := ln(x) [f (x)g(x) + f(x)g (x)]dx (3.) + ln( x) [f (x)g(x) + f(x)g (x)]dx. (3.) Note that up to an integration by parts, and for functions supported away from zero, we have C (f, g) = f, x g. But this is of course not true for functions which are not zero at the origin. The main result os this subsection is contained in the following proposition: Proposition 3.. For r < and for every f H (R) we have the estimate ( r 4 9 f, V r efff = ln(r) f() + C (f, f) + O dy ) f H (R). (3.3)

12 Proof. The argument is a bit long, and we split it in several lemmas. Let us start by listing some of the properties of Veff r. First note that it scales like a delta function, i.e. it is homogeneous of order : Veff r (x) = r V eff (x ), (3.4) r The next observation is that due to the integral with respect to y it is much less singular than V r : Lemma 3.. The behaviour of Veff r (x) is logarithmic at. Proof. There exists a constant c > large enough such that for every y π and x R we have c (x + y ) x 4 + sin (y) c(x + y ). Therefore, at small x, Veff behaves like: π Veff (x) ( ) x + y dy = ln π π x ln( x ) x ln( x ) + O(). (3.5) We now define on R a comparison function Y r (x) := x +4r. and we also denote by Y r the associated quadratic form defined on H (R). For the following, let us recall an important Sobolev estimate in one dimension, f L (R) const f H (R) (3.6) which comes from the fact that f(x) (p + i) (p + i) ˆf(p) dp const( f + f ). (3.7) π R Lemma 3.3. We have the following properties: (i) V eff Y, V eff (x) Y (x) = O( x 5 ) for x, and V eff Y L (R) = ln(4); (3.8) (ii) f, (Veff r Y r)f = ln(4) f() + O(r 4 9) f H (R) r.

13 Proof. (i). To show that Veff Y, one uses sin( ). The second estimate for x follows from: V eff(x) Y (x) = x = x π ( π π π ( π + 4 sin y x π ) dy ( + 4x ) sin y x dy + O(x 4 ) + x + O(x 4 ) = O( x 5 ), x. (3.9) Before computing the L norm of (3.8), let us notice that none of the terms is in L. We first integrate over x and get: π Veff Y L (R) = R π π x + 4 sin y π dydx x + 4 = 4 π π then we do the integral with respect to y: π ln(siny)dy = π ln + π = π ln() + π 4 ln(sin y)dy, (3.) ln(sin y )dy + π = π ln() + π 4 ln[siny]dy + ln[siny]dy + ln(cos y )dy π 4 π π 4 ln[sin(π/ y)]dy ln[sin y]dy = π ln(), (3.) ), so (3.8) is proved. Let us now prove (ii). We have due to the scaling properties: f, (Veff r Y r )f = (Veff Y )(x/r) f(x) dx r R = (Veff Y )(x) f(rx) dx. (3.) R Then, we substract the term V eff Y L f() = ln(4) f() which gives f, (Veff r Y r)f ln(4) f() = { } = (Veff Y )(x) [f(rx) f()]f(rx) + [f(rx) f()]f() dx. (3.3) R 3

14 Let α (, ) a real number. We split the above integral in two regions: x r α and x r α. We have, using (3.6): (Veff Y )(x)[f(rx) f()]f(rx)dx (3.4) x r α f H (V eff Y )(x) dx x r α f H x r x 5 dx if r α α O(r 4α ) f H. (3.5) For the region x r α (and similarly for the other one), we can write: (Veff Y )(x) f(rx) f() f(rx) dx (3.6) x r α r α rx (Veff Y )(x) f (t)dt f(rx) dx V eff Y L f L and the Cauchy-Schwarz inequality yields: r α r α f (t) dt f (t) dt r α f H. Then we set α as the solution of ( α)/ = 4α which gives α = 9. We now concentrate ourselves on Y r when r is small. For the next two lemmas, we need to introduce the the following caracteristic function: { if x χ(x) = otherwise. (3.7) Then we have the following lemma: Lemma 3.4. Consider the self-adjoint operator of multiplication by ln( )χ defined on its natural domain in L (R). This operator is relatively bounded to p x := id/dx, with relative bound zero. Proof. Indeed, ln( x )χ( x ) (p x + iλ), λ > is Hilbert-Schmidt since we have, from [RS3, XI.3]: ln( )χ (p x + iλ) HS const ln( )χ L ( + iλ) L const λ. (3.8) 4

15 Note that by a similar argument as the one in (3.6) we get the estimate: (p x + iλ) L L Then a standard argument finishes the proof. const λ. (3.9) We can now characterize the form C introduced in (3.): Lemma 3.5. The quadratic form induced by C admits a continuous extension to H (R). Moreover, C is infinitesimally form bounded with respect to the form associated to p x = d /dx. Proof. Fix some ǫ (, ). Then for every f S(R) we can write: C (f, f) = ln(x) (d x f )(x)dx + ln( x) (d x f )(x)dx First we have Then using (3.9) we have = ln( x) (d x f )(x)dx ε + ln(ε) ( f(ε) + f( ε) ) + R\[ ε,ε] ε R\[ ε,ε] x f(x) dx ε f. ln(x) (d x f )(x)dx x f(x) dx. (3.) sup f(t) const (p x + iλ)f, (3.) t R λ which takes care of the terms containing f(±ε). The remaining two integrals can be treated with the following argument: ε ε ln(x) (d x f )(x)dx ln(x) f(x) f (x) dx χ ln( )f f H const λ (p x + iλ)f f H, (3.) where in the second inequality we used the Cauchy inequality, while in the third inequality we used (3.8). These estimates allow us to find two constants A and B (growing when λ grows, such that: C (f, f) const λ f + A f f + B f. 5

16 But we use the inequality f f λa f + λa f, which finally allows us to say that for any < a < we can find b > such that and the proof is over. C (f, f) a p x f + b f, (3.3) The final ingredient in proving Proposition 3. is contained in the following lemma: Lemma 3.6. Recall that Y r (x) = (x + 4r ) /. Then for every r <, and for every f S(R), we have the estimate: ( ) f, Y r f = ln(r) f() + C (f, f) + O r 4 9 f H (R). Proof. Integrating by parts we obtain: and: f, Y r f = ln(r) f() + ln(x + x + 4r ) [f (x)f(x) + f(x)f (x)]dx ln( x + x + 4r ) [f (x)f(x) + f(x)f (x)]dx (3.4) f, Y r f C (f, f) + ln(r) f() = (3.5) + [ln(x + x + 4r ) ln(x)] (d x f(x) )dx [ln( x + x + 4r ) ln( x)] (d x f(x) )dx. (3.6) The idea is to show that the last two integrals are small when r is small. We only consider the first integral, since the argument is completely analoguous for the second one. Fix some < α < (its optimal value will be chosen later), and assume that r is small enough such that r α /. We split the domain of integration into two regions: one in which x > r α, and the other one where x r α. For the first region we have: ln(x + x + 4r ) ln(x) = ln r x = O(r /x ). (3.7) 6

17 Then by integration, and using (3.6) together with the Cauchy inequality, we get: [ln(x + x + 4r ) ln(x)]f (x)f(x)dx r α [ln(x + x + 4r ) ln(x)]f (x)dx f H r α ( [ln(x + ) x + 4r ) ln(x)] dx f H r α O(r 3 α ) f H. (3.8) For the region where x r α, we use the monotonicity of the logarithm and write: ln(x + x + 4r ) ln(x) ln(x + x + 4r ) + ln(x) Then we can write r α [ln(x + x + 4r ) ln(x)]f (x)f(x)dx ( r α const f H [ln(x) + ln(r)] dx ln(r) + ln(x). (3.9) ) = O(r α lnr ) f H. (3.3) Compairing (3.8) and (3.3), we see that we can take α arbitrarily close to. In particular, we can find some α such that 3α/ > 4/9 and α/ > 4/9 and we are done. We can now conclude the proof of Proposition 3. by putting together the estimates from Lemma 3.3 and Lemma A solvable comparison operator H C Remember that we are interested in the negative spectrum of Heff r, operator given in (.8). Lemmas 3. and 3.4 tell us that its operator domain is H (R), while the form domain is H (R). Proposition 3. indicates that a good approximation for Heff r at small r would be the operator formally defined as H C := p x + ln(r)δ x. (3.3) 7

18 Of course, as it is written above H C makes no sense. It has to be understood in the following way: consider the sesquilinear form on S(R) given by t C (f, g) := f, g + ln(r) f()g() C (f, g). (3.3) A standard consequence of (3.) and (3.3) is that the quadratic form associated with t C is closable, bounded from below, and the domain of its closure is H (R). Then H C is the self-adjoint operator generated by t C, and its operator domain D C is characterised by: D C := {ψ H (R) : t C (φ, ψ) const φ, φ H (R)}. (3.33) Moreover, if ψ D C, then we have the equality: t C (φ, ψ) = φ, H C ψ, φ H (R)}. (3.34) Another representation for ψ D C is that there exists f ψ = H C ψ L (R) such that the distribution ψ is a regular distribution on R \ {} and is given by: ψ (x) = x ψ(x) f ψ(x). (3.35) One important consequence is that ψ H (R \ {}), and ψ is continuous on R \ {}. Let us now introduce the parity operators P + and P P ± : H (R) H (R), f(x) (P ± f)(x) = f(x) ± f( x). (3.36) We have that P + + P =. We will call RanP + the even sector and RanP the odd sector. The following lemma is an easy application of definitions, and we give it without proof: Lemma 3.7. We have that t C (P ± f, P g) = for all f, g H (R). Then the domain D C of H C is left invariant by P ± ; moreover, H C commutes with P ±. A standard consequence of the elliptic regularity (see also (3.35)) is the following lemma, given again without proof: Lemma 3.8. The eigenvectors of H C belong to C (R\{}). A less obvious result is the following characterisation of D C : Lemma 3.9. Every ψ in D C obeys the following boundary condition at : [ ψ ( ε) ψ ] (ε) lim + ln(r)ψ() ln(ε)ψ() =. (3.37) ε 8

19 Proof. Now remember that ψ is continuous outside the origin if ψ D C, so ψ (±ε) makes sense. Moreover, for every φ H (R) we have: φ (x)ψ (x)dx + ln(r)φ()ψ() C (φ, ψ) = φ, H C ψ. (3.38) R We can write for x > (see (3.35)): ψ (x) = ψ () + x and then for x close to we have: ψ (y)dy = ψ () x ( ψ(y) y ) + f ψ (y) dy ψ (x) const + lnx ψ. (3.39) The same estimate is true for negative x near, and together with the estimate (3.6) it follows that ψ diverges at not faster than a logarithm. Now we can integrate by parts outside the origin and write: φ (x)ψ (x)dx (3.4) R = φ( ε)ψ ( ε) ε + φ (x)ψ (x)dx ε ε φ(x)ψ (x)dx φ(ε)ψ (ε) where the last integral will converge to zero with ε. After a similar integration by parts we obtain: C (φ, ψ) = ln( x) (d x (φψ))(x)dx ε ( ) + ln(ε) φ(ε)ψ(ε) + φ( ε)ψ( ε) + ε R\[ ε,ε] Following the reasoning in (3.), one can prove that : and thus: ε C (φ, ψ) ε = ln(ε) ε φ(x)ψ (x)dx ln(x) (d x (φψ))(x)dx φ(x)ψ(x)dx. (3.4) x ln(x) [d x (φψ)](x)dx ε = O(ε lnε ) φ H ψ H ( ) φ(ε)ψ(ε) + φ( ε)ψ( ε) + R\[ ε,ε] x φ(x)ψ(x)dx. (3.4) 9

20 Putting (3.4) and (3.4) in (3.38), and using (3.35), we eventually get: [ ( ) lim φ( ε)ψ ( ε) φ(ε)ψ (ε) + ln(r)φ()ψ() ε ( )] ln(ε) φ(ε)ψ(ε) + φ( ε)ψ( ε) =. (3.43) The last ingredient is the embedding H (R) C / δ (R), and the estimate (3.39) which allows us to simplify the limit: [ ψ lim φ() ( ε) ψ ] (ε) + ln(r)ψ() ln(ε)ψ() = (3.44) ε for all φ H (R). The lemma is proved. 3. The eigenvalues and eigenvectors of H C In this subsection we give analytic expressions for eigenvalues and eigenvectors of H C corresponding to the negative, discrete spectrum; much of the information we need about special functions is borrowed from [L]. We want to have the same formal expression for our eigenvalue problem outside z = as in that paper, namely d dz ψ 4 ψ + α z ψ = where ψ will be an eigenfunction with an associated eigenvalue E = α. Let us now do this in a rigorous manner. We want to implement the change of variables x = αz, α >, which amounts to defining a unitary operator on L (R): (U W ψ)(z) = α ψ(αz/), (U W ψ)(z) = α Now assume φ is a normalized eigenvector for H C satisfying ψ(z/α). (3.45) H C φ = Eφ, E <. (3.46) Instead of solving the above equation, we will reformulate it in terms of ψ = U W φ, and φ = U W ψ. To do that, we need to fulfill two conditions. The first one is: [U W HU W ψ](z) = α [ ψ (z) α ψ(z)] = Eψ(z), z, or z ψ (z) 4 ψ(z) + α ψ(z) =, z, E = z α. (3.47) The second thing is to see what condition at z = should ψ obey in order to be sure that U W ψ is in the domain of H C. If we replace ψ in (3.37) by φ = U W ψ, then we get the modified condition:

21 [ ψ ( α lim ε) ψ ( α ε) ] + ln(r)ψ() ln(ε)ψ() = (3.48) ε α or [ ψ ( ε) ψ ] (ε) lim + α(ln r ln(αε))ψ() =. (3.49) ε Therefore we reduced the problem of finding the eigenfunctions and eigenvalues of H C to solving the ordinary differential equations in (3.47), with the boundary condition given in (3.49). We will see that L solutions can be constructed only if α obeys some conditions. A priori, α can be any positive real number. First assume that α is not a positive integer. Then if we solve (3.47) for z >, we see that the only square integrable solution at + is the one given by a Whittaker function: W α, (z) = ze z U( α,, z), (3.5) where U is the confluent hypergeometric function or Kummer function, see [AS]. If α is a positive integer, the solution is obtained as the limit of W α, (z) when α tends to a positive integer N and get: lim W α, α N (z) = e z z N L N (z) (3.5) where L N is an associated Laguerre polynomial. We denote with Γ(z) and ψ(z) = Γ (z)/γ(z) the usual gamma and digamma functions. We have the following first result: Proposition 3.. (i). All negative eigenvalues of H C are nondegenerate. The eigenfunctions of H C are also eigenfunctions of P ±. (ii). There exists an infinite number of odd eigenfunctions φ odd,k, k Z +, corresponding to every α {,,... }. (iii). There also exists an infinite number of even eigenfunctions φ even,k, k Z +, each corresponding to a certain α k (k, k) for every k Z +. Proof. (i). Choose any eigenfunction φ of H C corresponding to E <. Make the change ψ = U W φ, and then look at the associated differential equation: ψ (z) 4 ψ(z) + α ψ(z) =, z, E = z α, (3.5) [ ψ ( ε) ψ ] (ε) lim + α(ln r ln(αε))ψ() =. (3.53) ε First assume α = N Z +. The theory of ordinary differential equations insure the existence of two constants C and C such that ψ(z) = C e z z N L N (z), z >, ψ(z) = C e z z N L N ( z), z <. (3.54)

22 By inspection (and by continuity) we get that ψ() =. If we put ψ in (3.53), and using the explicit form of the Laguerre polynomials, we get that the boundary condition is fulfilled only if C = C which amounts to ψ(z) = ψ( z), i.e. there is one and only one solution which is also odd. Now assume α Z +. The theory of ordinary differential equations insure the existence of two constants C 3 and C 4 such that ψ(z) = C 3 W α, (z), z >, ψ(z) = C 4 W α, ( z), z <. (3.55) By inspection we see that W α, () = Γ( α), hence by continuity at zero we must have C 3 = C 4, hence ψ(z) = ψ( z) and the eigenfunction must be even. (ii). The proof is already contained in (i), since the boundary condition is trivially fulfilled for odd functions. There is exactly one eigenfunction, an odd one, coresponding to every α Z +. (iii). We saw in (i) that if there are eigenfunctions corresponding to α Z +, then they must be even. In order to get all possible α s which are compatible with the boundary condition (3.53), we compute (note that ψ is odd): [ ] lim W ε α, (ε) + α(ln r ln(αε)) = (3.56) Γ( α) and using the explicit expression of these special functions we obtain the condition on α: f(α, r) := ψ( α) + γ + lnα + lnr =. (3.57) α where ψ here means the digamma function and γ is Euler s constant. Since the digamma function is stricly increasing from to + on each interval of the form ( m, m + ), m Z +, one can easily see that f(, r) is stricly decreasing from + to when α varies in an interval of the form (k, k) for every k Z +. Therefore we have a unique solution α k (k, k) of the equation f(α k, r) for every k Z +. The proposition is proved. The previous proposition stated that only the eigenvalues from the even sector can vary with r. Let us now further investigate this dependence. Corollary 3.. (i). The excited states with even parity tend to those with odd parity when r is small. More precisely, for k, we have that lim r α k = k ; (ii). For k =, we have the following asymptotic behavior of the ground state: α (r) = ln(r) { + o r()}, E (r) = [ln(r)] { + o r ()}. (3.58)

23 Proof. (i). The limit follows easily from the properties of the digamma function. (ii). We apply the implicit function theorem. Define the function F(α, y) := α + α[γ + ψ( α) ln(α)] y, for (α, y) in a small disk around the origin in R. This function is C near (, ), ( α F)(, ) =, and F(, ) =. Then for every y > small enough there exists α(y) > such that F(α(y), y) =. Now put y = / ln(r) and we are done because (3.57) is also satisfied with this α. 3.3 Approximation of H r eff by H C We will now show that the negative spectrum of Heff r converges in a certain sense to the one of H C. This is made precise in the next proposition, but before we need a definition. For a given subset S of R, and for any ǫ > we define S ǫ := x S B ǫ (x). (3.59) If S is a discrete, finite set, then S ǫ is a finite union of intervals of length ǫ, centered at the points of S. Proposition 3.. The following three statements hold true: (i). Fix a <, and denote by A := σ(h C ) (, a] and B := σ(h r eff ) (, a]. Then for every ǫ >, there exists r ǫ > such that for every r < r ǫ we have A B ǫ, B A ǫ. (3.6) (ii). The groundstate of Heff r is nondegenerate, has even parity, and diverges to when r. Moreover: lim inf σ(h C) inf σ(heff r ) =. (3.6) r (iii). Fix a compact interval [a, b] (, ) and suppose that H C has exactly one eigenvalue of a given parity E C in [a, b], for all r < r. Then Heff r has exactly one eigenvalue of the same parity E eff in this interval and lim E eff E C =. r Proof. Let us introduce the resolvents R eff (z) = (H r eff z) for all z ρ(h r eff ) and R C (z) = (H C z) for all z ρ(h C ). The first ingredient in the proof is contained by the following lemma: 3

24 Lemma 3.3. There exists a constant K > sufficiently large, and r small enough, such that for every r < r we have that the form defined on L (R) L (R) (see also (3.3)) ) V C (f, g) := t C ([p x / + λ r] / f, [p x / + λ r] / g + λ r f, [p x / + λ r] g f, g = [p x / + λ r] / {H C + λ r }[p x / + λ r] / Id, λ r := K ln (r). (3.6) generates a bounded operator on L (R) denoted in the same way. Moreover, sup <r<r V C /. And we have {H C + λ r } = [p x / + λ r] / {Id + V C } [p x / + λ r] /. (3.63) Proof. The key estimate is contained in (p x/ + λ) / L L const, (3.64) λ/4 obtained with an argument as in (3.9). Then if we have ln(r) / λ r small enough, then the delta function part of t C will be small uniformly in r < r. Using also (3.3), and the definition (3.3), then one can show that (3.6) is a bounded sesquilinear form on L (R), with a norm which can be made arbitrarily small if K is chosen large enough. Now the equality (3.63) is easy, and note that this is also compatible with (3.58). Introduce the notation: Ṽ eff := (p x/ + λ r ) / V r eff(p x/ + λ r ) /. (3.65) The second ingredient in the proof of the above proposition is the following estimate, which is an easy consequence of Proposition 3.: V C Ṽeff = O(r 4/9 ), r < r. (3.66) We also have that Ṽeff /3 if r is small enough, uniformly in r < r, and then {H r eff + λ r} = [p x / + λ r] / {Id + Ṽeff} [p x / + λ r] /. (3.67) It is clear that a similar identity would hold for any other λ λ r, and this already tells us that the spectrum of H r eff is contained in an interval of the type ( λ r, ), thus justifying the discussion after (.4). From (3.67), (3.66) and (3.63), we get the crucial estimate: R eff ( λ r ) R C ( λ r ) const r4/9 λ r, r < r. (3.68) 4

25 This estimate allows us to prove (i). Introduce the notation d C (z) := dist(z, σ(h C )). (3.69) Choose z ρ(h C ) (thus d C (z) > ). From the identity: (H C z)r C ( λ r ) = Id (z + λ r )R C ( λ r ) (3.7) we get that the right hand side is invertible and: {Id (z + λ)r C ( λ r )} = (H C + λ r )R C (z) = Id + (z + λ r )R C (z). (3.7) The first equality implies that R C (z) = R C ( λ r ) {Id (z + λ r )R C ( λ r )}, (3.7) while the second one gives the norm estimate: {Id (z + λ r )R C ( λ r )} [ + ( z + λr )/d C (z)]. (3.73) Note the important fact that (3.7) is just another form of the resolvent identity, valid for any self-adjoint operator. If we could replace R C ( λ r ) by R eff ( λ r ), then the right hand side would immediately imply that z ρ(h eff ). We can restrict ourselves to those z s which obey z λ r. Then using (3.68) and (3.73), we get that for r < r and d C (z) r /3, the operator Id (z + λ r )R eff ( λ r ) is invertible and we get the estimate: {Id (z + λ r )R eff ( λ r )} const d C (z) λ r [ const λ r r 4/9 /d C (z)]. (3.74) Therefore we have proved that for every z which obeys z λ r and d C (z) > r /3, the operator R eff ( λ r ) {Id (z + λ r )R eff ( λ r )} (3.75) exists and defines R eff (z). It means that the spectrum of Heff r is close to that of H C, and the distance between them is going to zero at least like r /3. Let us now prove (ii). We know that the groundstate of H C diverges like ln (r) for small r, and it is isolated from the rest of the spectrum. Choose a circular contour Γ of radius around this groundstate. It means that d C (z) = for z Γ. Then (3.75), (3.74) and (3.68) imply the estimate sup R eff (z) R C (z) const r 4/9, r < r. (3.76) z Γ 5

26 Now we can employ the regular perturbation theory, see [K], by using Riesz projections defined as complex integrals of the resolvents on contours like Γ. Then the estimate (3.6) is straightforward. Finally, let us prove (iii). We know that for small r, the excited states of H C tend to cluster in pairs. The eigenvalues from the odd parity sector are independent of r, while those from the even parity sector will converge from above to the odd ones (see Proposition 3.). Consider such a pair of eigenvalues, which will always remain separated from the rest of the spectrum if r < r and r is small enough. Then we can find a contour Γ which contains them and inf z Γ d C (z) is bounded from below uniformly in r < r. Then we can again write an estimate like (3.76), and then apply the regular perturbation theory. The proof is over. 4 Reduction of H r to H r eff We are now ready to go back to (.4), and argue why only the diagonal entries of the infinite operator-valued matrix {H m,n } m,n Z are important for the low lying spectrum of H r. Let us formally write H r as: H r = H diag + V offdiag, where H diag = n Z (Hr eff + n r ), and V offdiag contains all the non-diagonal entries of the form Vm,n r, m n, (see (.5)), and zero on the diagonal. We will prove in this section that V offdiag is relatively form bounded with respect to H diag, and moreover, it is a small perturbation when r is small. The main result is very similar to Proposition 3., where we only have to change H C by H diag and Heff r by H r. Parity here only refers to the x variable. Therefore we will start comparing the two operators. Before that, let us note that the negative spectrum of H diag is given by Heff r if r is small enough. 4. V offdiag is H diag -form bounded Let λ r = K ln (r) with K large enough and r < r, as in the previous section. We know that λ r ρ(h diag ), and denote by R diag ( λ r ) the resolvent (H diag +λ r ). Then the main technical result of this subsection will be the following estimate: there exists δ > and r (δ) such that R / diag ( λ r)v offdiag R / diag ( λ r) = B (l (Z;L (R))) O(rδ λr / ), r < r. (4.) Let us first notice that we can replace R / diag by a simpler operator, namely n Z (ǫp x + + n r ) /, where ǫ is a small enough positive number. Indeed, we can write f, [H r eff + λ r + n /(r )]f f, [ǫp x + λ r/ + n /(r )]f, 6

27 where we used that for ǫ small enough we can show that: (/ ǫ)p x V r eff + λ r /, r < r. This means that [ǫp x + λ r/ + n /(r )] / [Heff r + λ r + n /(r )] /. (4.) B (L (R)) Define the bounded operators in L (R) (see (.5)): Ṽ r m,n := [ǫp x + λ r/ + m /(r )] / V r m,n [ǫp x + λ r/ + n /(r )] /, m n, Ṽ r m,m :=, m Z. (4.3) Then (4.) would be implied by the following, stronger estimate: } r {Ṽ m,n = O(r δ λr / ), r < r. (4.4) m,n Z B (l (Z;L (R))) By an easy application of the Schur-Holmgren lemma, one can prove the estimate: } r {Ṽ m,n sup Ṽ m,n r B (L (R)). (4.5) B (l (Z;L (R))) m,n Z m Z n Z We now concentrate on the norms Ṽ r m,n B (L (R)) and study their behavior in r, m, and n. Remember that only the case m n is of interest, since the diagonal terms are zero. Before anything else, let us do a unitary rescaling of L (R) by (Uf)(x) := r / f(rx) and (U f)(x) := r / f(x/r). Then due to various homogeneity properties we get: UṼ r m,n U = r [ǫp x + r λ r / + m /] / V m,n [ǫp x + r λ r / + n /] /. (4.6) We first give an important estimate for V m,n, stated in the next lemma: Lemma 4.. Let < α < and m n. Fix any < ǫ <. Then there exists a constant C = C(α, ǫ) such that we have the following estimate: { } Vm,n (x) C n m α x α +, x r ǫ, (4.7) n m and r 3ǫ Vm,n (x) const n m, x r ǫ. (4.8) 7

28 Proof. Due to symmetry properties we can write V m,n(x) = π Integrating by parts we get: V m,n(x) = π π π(n m) cos[(n m)y] [x + 4 sin dy. (4.9) (y/)] / π π sin[(n m)y] sin(y) [x + 4 sin (y/)] 3/dy. (4.) This equality immediately proves (4.8). So we now focus on x r ǫ. We can split the integral in two: one in which y π/, and where the integrand has no singularities when x is small, and the second where y π/. In that region we can use the same idea as in Lemma 3. of replacing sin (y/) by y. We hence get: ( ) Vm,n const π/ sin[(n m)y] sin(y) (x) + dy. (4.) n m π/ [x + y ] 3/ Now we employ the inequalities (here < α < is arbitrary): sin[(n m)y] n m α y α, sin(y) y, then we make the change of variables s = y/ x and write: Vm,n (x) const n m α s ( α ) + n m x α ds. (4.) [ + s ] 3/ Thus the lemma is proved. Now let us go back to (4.6), and estimate the various norms. If we write V m,n = V m,n χ( r ǫ ) + V m,n χ( > r ǫ ), then we have two different types of estimates. When we keep V m,n χ( > r ǫ ), which is bounded, then for the two resolvents we can use the usual B (L ) norm, which together with (4.8) gives a contribution: const r 3ǫ r λ r + n, n m. (4.3) r λ r + m m n When we keep V m,n χ( r ǫ ), the estimate from (4.7) gives us that [ V m,n χ( )]/ is an L function, hence the operator V m,n χ( r ǫ )[ǫp x + r λ r / + n /] / 8

29 is Hilbert-Schmidt. Thus we have a product of two Hilbert-Schmidt operators, and we can give an upper bound for the B norm of their product of the form: const(α) r ǫ (r λ r + n ) /4 (r λ r + m ) /4 m n α, n m. (4.4) Therefore we obtained an upper bound for the norm of the operator in (4.6) of the form: Ṽ m,n r r const r 3ǫ r λ r + n r λ r + m m n r const r ǫ + (r λ r + n ) /4 (r λ r + m ) /4 m n α, m n. (4.5) Remember that one is interested in the right hand side of (4.5). Now choose / < α <. We have to investigate several cases:. When m = and n. Then the first term in (4.5) will behave like λ / r r 3ǫ n. The second term will behave like r / ǫ λr /4 n / α. Both contributions are summable with respect to n. Note that if ǫ is small enough, both exponents of r are positive. Denote by δ the smaller one.. Fix m, and consider all n m. When n =, we get similar terms as above. If n, then we remain with the problem of summing up something like sup m / m n,n m n / n m α, / < α. We can either use Hölder s inequality, or we can split the above sum in the following way: n / n m α n,n m = + n / n m α n,n m, n n m n,n m, n > n m ( ) n + const(α). (4.6) α+/ n m /+α n,n m We therefore consider (4.) as proved. 9

30 4. Comparison between H r and H diag If r is small enough, we have the identity: ( H { } r + λ r ) = R / diag ( λ r) Id + R / diag ( λ r)v offdiag R / diag ( λ / r) R diag ( λ r). Moreover, this implies: ( H r + λ r ) R diag ( λ r ) const λ r R / diag ( λ r)v offdiag R / diag ( λ r) const rδ λr /. (4.7) λ r This is the same as what we had in (3.68), but with λ / r r δ instead of r 4/9. Therefore we can repeat the arguments of Proposition 3. and prove a similar kind of spectrum stability for H r and H diag. 5 The main theorem and some conclusions It is particularly hard to state a concentrated main result of our paper. Let us now go back to the very first Hamiltonian which was declared to be relevant for the exciton problem. This is H r, written in (.). Because of the Coulomb singularity, the best way to look at the spectral problem is to consider its form t H, given by (.4), (.5) and (.6). We then managed to separate the mass center motion in the longitudinal direction, and we got a simpler form t h in (.8). The center of the mass cannot be separated in the transverse direction because of the cylindrical geometry, but at least we can write t h as a direct sum of k Z t h k. A crucial observation has been stated in (.4), which says that only t h is responsible for the lowest lying spectrum of the original form. This gave us the possibility of renaming t h with t H in (.5), and declare it as the central object of study. Then in Proposition. we constructed its associated self-adjoint operator H r, where we had to take care of the Coulomb-type singularity in two dimensions. Then after a unitary transformation induced by the discrete Fourier transform with respect to the y variable, we can see H r as an infinite operator valued matrix acting on the Hilbert space l (Z; L (R)). We then decomposed H r as the sum of a diagonal operator H diag and an off-diagonal part V offdiag. Eventually we proved in Section 4 that the low lying spectrum of H r is only slightly influenced by the off-diagonal part for small r, and therefore the relevant object remains H diag. But this diagonal part has the nice feature that each of its entry is of the form Heff r + n r, n Z, where H r eff is given in (.8) and (.7). Then in Section 3, more precisely in Proposition 3. we prove that the low lying spectrum of Heff r is well approximated by the spectrum of a solvable operator, H C, which we discussed in Proposition 3.. 3

31 We are now ready to collect all these results in the main theorem of our paper: Theorem 5.. The following three statements hold true: (i). Fix a <, and denote by A := σ(h C ) (, a] and B := σ( H r ) (, a]. Then for every ǫ >, there exists r ǫ > such that for every r < r ǫ we have A B ǫ, B A ǫ. (5.) (ii). The groundstate of H r is nondegenerate, and diverges to when r. The corresponding eigenfunction has even parity with respect to both variables. Moreover: lim inf σ(h C) inf σ( H r ) =. (5.) r (iii). Fix a compact interval [a, b] (, ) and suppose that H C has exactly one eigenvalue E C in [a, b], of parity p = ±, for all r < r. Then if r is small enough, H r has exactly one eigenvalue Ẽ in this interval and lim Ẽ E C =. r Moreover, the corresponding eigenfunction has even parity with respect to y and parity p with respect to x. Another important aspect of this problem is to determine how fast these limits are assumed. We have not touched this issue here, but we will study the numerical and physical implications of our results in a consequent paper. Acknowledgments. The authors thank T.G. Pedersen for many fruitful discussions. H.C. acknowledges support from the Danish F.N.U. grant Mathematical Physics and Partial Differential Equations, and partial support through the European Union s IHP network Analysis & Quantum HPRN-CT--77. References [AS] [Aki] [Bo] Abramowitz M., Stegun I. A., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover, New York, 964. Akiyama H.: One dimensional excitons in GaAs quantum nanowires, J. Phys.: Condens. Matter, (998). Bouzouina, A.: Stability of the two-dimensional Brown-Ravenhall operator. Proc. Roy. Soc. Edinburgh Sect. A 3 no. 5, (). 3

32 [BD] [CDP] [FW] [K] [KCM] Brummelhuis R., Duclos P.: Effective Hamiltonians for atoms in very strong magnetic fields. Few-Body Systems 3, -6 (). Cornean, H.D., Duclos, P., Pedersen, T.G.: One dimensional models of excitons in carbon nanotubes, Few-Body Systems 34 no. -3, 55-6 (3). Fetter, A.L., Walecka, J.D.: Quantum theory of many-particle systems, Dover, 3. Kato T.: Perturbation theory for linear operators, Springer-Verlag, 966. Kostov M. K., Cole M. W., Mahan G. D.: Variational approach to the coulomb problem on a cylinder, Phys. Rev. B 66, 7547 (). [L] Loudon R.: One-dimensional hydrogen atom, Amer. J. Phys. 7, (959). [Ped] Pedersen, T.G.: Variational approach to excitons in carbon nanotubes, Phys. Rev. B 67, 734 (3). [Ped] Pedersen, T.G.: Exciton effects in carbon nanotubes, Carbon 4, 7 (4). [RS3] [SS66] Reed, M., Simon, B.: Methods of modern mathematical physics. III: Scattering theory, Academic press, 97. Shinada M., Sugano, S.: Interband optical transitions in extremely anisotropic semiconductors. I. Bound and unbound exciton absorption, J. Phys. Soc. Jpn. no.7, 936 (966). 3