arxiv: v2 [math-ph] 17 Sep 2010

Transcription

1 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS arxiv: v [math-h] 17 Se 010 ROBERT SEIRIGER Abstract. We investigate the low energy excitation sectrum of a Bose gas with weak, long range reulsive interactions. In articular, we rove that the Bogoliubov sectrum of elementary excitations with linear disersion relation for small momentum becomes exact in the mean-field limit. 1. Introduction and Main Results Bogoliubov s seminal 1947 aer [1] contains several imortant results concerning the low energy behavior of bosonic systems. Among its striking redictions is the fact that the excitation sectrum is made u of elementary excitations whose energy is linear in the momentum for small momentum. Bogoliubov s method is based on various aroximations and crucially uses a formalism on Fock sace that does not conserve article number. Mathematically, the validity of his method has so far only been established for the ground state energy of certain systems, see [1, 13, 14, 4, 8, 16]. In articular, there are no rigorous results on the low energy excitation sectrum of interacting Bose gases, with the notable excetion of exactly solvable models in one dimension [7, 9, 10,, 15]. In this article, we shall rove the validity of Bogoliubov s aroximation scheme for a Bose gas in arbitrary dimension in the mean-field Hartree limit, where the interaction strength is roortional to the inverse article number, and its range extends over the whole system. In articular, we verify that the low energy excitation sectrum for such a system equals the sum of elementary excitations, as redicted by Bogoliubov. As a corollary, we observe that the lowest energy in the sector of total momentum P deends linearly on P, a roerty that is crucial for the suerfluid behavior of the system. The mean-field limit has served as a convenient and instructive toy model for several asects of bosonic systems over the years. We refer to [5, 6] for a review and further references. We consider a homogeneous system of bosons on the flat unit torus T d, d 1. The bosons interact with a weak two-body interaction which we write for convenience as 1 vx. We assume that v is ositive, bounded, eriodic with eriod one, and v x = vx. We alsoassume that v is of ositive tye, i.e., it has only non-negative Fourier coefficients. With denoting the usual Lalacian on T d, the Hamiltonian equals H = i + 1 vx i x j i=1 i<j Date: Set. 17, 010. c 010 by the author. This work may be reroduced, in its entirety, for non-commercial uroses. 1

2 ROBERT SEIRIGER in suitable units. It acts on L sym Td, the ermutation-symmetric square integrable functions of variables x i T d. Let E 0 denote the ground state energy of H. The Bogoliubov aroximation [1, 11] redicts that E 0 is close to 1 v0+ebog, where E Bog = 1 + v 4 + v. 1 The sum runs over πz d, and v = vxe ix dx T d are the Fourier coefficients of v. ote that the sum above converges, since the summands behave like v / for large. More imortantly, the Bogoliubov aroximation redicts that the excitation sectrum of H is made u of elementary excitations of momentum with corresonding energy e = 4 + v. One noteworthy feature of is that it is linear in for small, in contrast to the case when interactions are absent. Our main results can be summarized as follows. THEOREM 1. The ground state energy E 0 of H equals E 0 = v0+ebog +O/, 3 with E Bog defined in 1. Moreover, the sectrum of H E 0 below an energy ξ is equal to finite sums of the form e n +O ξ 3/ /, 4 πz d \{0} where e is given in and n {0,1,,...} for all 0. The error term O/ in 3 refers to an exression that is bounded, in absolute value, by a constant times / for large, where the constant deends only on the interaction otential v; likewise for the error term Oξ 3/ / in 4. The deendence on v is rather comlicated but can be deduced from our roof, which gives exlicit bounds. Keeing track of this deendence allows to draw conclusions about the sectrum for large even if v deends on. In the case of a fixed, -indeendent v, Theorem 1 imlies that the Bogoliubov aroximation becomes exact in the mean-field Hartree limit. As long as ξ 1/3, each individual excitation energy ξ is of the form e n with error o1. Moreover, as long as ξ, it is of the form e n 1 +o1, i.e., the error is small relative to the magnitude of the excitation energy. In other words, in the mean field limit the whole excitation sectrum with energy ξ is given in terms of sums of elementary excitations. The condition ξ can be exected to be otimal, since only under this condition a large fraction of the articles are guaranteed to occuy the zero momentum mode, one of the key assumtions in the Bogoliubov aroximation. For excitation energies of the order and larger, the sectrum will not be comosed of elementary excitations anymore but has a more comlicated structure.

3 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS 3 Our roof shows that for each value of the {n } there exists exactly one eigenvalue of the form 4. Moreover, the eigenfunction corresonding to an eigenvalue with given {n } has total momentum n. Given this fact we readily deduce the following corollary from Theorem 1. Corollary 1. Let E P denote the ground state energy of H in the sector of total momentum P. We have E P E 0 = min e n +O P 3/ /. 5 {n }, n=p In articular, E P E 0 P min v+ +O P 3/ /. Thelineardeendence ofe P on P isofcrucialimortanceforthesuerfluid behavior of the Bose gas; see, e.g., the detailed discussion in [3]. It is a result of the interactions among the articles. The exression 5 differs markedly from the corresonding result for an ideal, non-interacting gas, esecially if v0 is large. Finally, we note that under the unitary transformation U = ex iq j=1 x j, q πz d, the Hamiltonian H transforms as U H U = H + q qp, where P = i j=1 j denotes the total momentum oerator. Hence our results aly equally also to the arts of the sectrum of H with excitation energies close to q, corresondingto collective excitations where the articles move uniformly with momentum q. The remainder of this aer is devoted to the roof of Theorem 1. The main strategy is to comare H with Bogoliubov s aroximatehamiltonian. The latter has to be suitably modified to take article number conservation into account.. Preliminaries We denote by P the rojection onto the constant function in L T d, and Q = 1 P. The oerator that counts the number of articles outside the zero momentum mode will be denoted by >, i.e., > = Q i. i=1 We shall also use the symbol T for the kinetic energy T = i=1 i. Lemma 1 below gives simle uer and lower bounds on the ground state energy E 0 of H, as well as an uer bound on the exectation value of T in a low energy state. Lemma 1. The ground state energy of H satisfies the bounds 0 E 0 v0 v0 v0. 6 Moreover, in any -article state Ψ with Ψ H Ψ v0+µ we have π Ψ > Ψ Ψ T Ψ v0 v0+µ. 7

4 4 ROBERT SEIRIGER Proof. The uer bound to the ground state energy follows from using a constant trial function. Since v 0, we have v e ixj 0. πz d \{0} This inequality can be rewritten as vx i x j v0 v0 8 and thus The rest follows easily. 1 i<j H v0+t j=1 v0 v0. In the following, we shall also need a bound on the exectation value of higher owers of >. More recisely, we shall use the following lemma. Lemma. Let Ψ be an -article wave function in the sectral subsace of H corresonding to energy E E 0 +µ. Then π Ψ > T Ψ v0+ µ + µ+3v0+ v0µ+v0 v0. In articular, Ψ > T Ψ is bounded above by an exression deending only on µ and v but not on. Proof. Since Ψ is ermutation symmetric, Ψ > T Ψ = Ψ Q1 S Ψ + Ψ > H E 0 1 µ Ψ where S = E µ 1 i<j vx i x j. Using Schwarz s inequality, the last term can be bounded as Ψ > H E 0 1 µ µ Ψ Ψ > 1/ Ψ. We slit S into two arts, S = S a +S b, with S a = E 0 + µ 1 vx i x j and S b = 1 i<j vx 1 x j. j= ote that S a does not deend on x 1. Using the ositivity of v as in 8, but with relaced by 1, as well as the uer bound on E 0 in 6, we see that In articular, this imlies that S a 1 µ+ v0+v0. Ψ Q 1 S a Ψ 1 µ+ v0+v0 Ψ > Ψ.

5 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS 5 To bound the contribution of S b, we use Ψ Q 1 S b Ψ = Ψ Q 1 vx 1 x Ψ = Ψ Q 1 Q vx 1 x Ψ + Ψ Q 1 P vx 1 x P Ψ + Ψ Q 1 P vx 1 x Q Ψ. The second term on the right side is ositive. For the first and the third, we use Schwarz s inequality and v = v0 to conclude Since Ψ Q 1 S b Ψ v0 Ψ Q 1 Q Ψ 1/ +v0 Ψ Q 1 Ψ. Ψ Q 1 Q Ψ Ψ > Ψ we have thus shown that Ψ > T Ψ 1 µ+ v0+3v0 Ψ > Ψ Using > π T this yields + v0+ 1 µ Ψ > Ψ 1/. Ψ > T Ψ v0+ 1 π µ The result then follows from Lemma 1. +µ+3v0+ v0 Ψ > Ψ. 3. The Bogoliubov Hamiltonian The main strategy in the roof of Theorem 1 is to comare the Hamiltonian H with the Bogoliubov Hamiltonian. We will use a slightly modified version of it, which is, in articular, article number conserving. Let a and a denote the usual creationand annihilation oeratorson Focksace for a article with momentum πz d, satisfying the canonical commutation relations [a,a q ] = δ q. For 0, let b = a a 0 9 and define the Bogoliubov Hamiltonian H Bog = [ b b + 1 v ] b b +b b +b b. 10 ote that H Bog conserves article number. We are interested in H Bog in the sector of exactly articles. ote also that v = v for all πz d. Let and A = + v B = v.

6 6 ROBERT SEIRIGER A simle comutation comare with [1, Thm. 6.3] shows that A b b +b b +B b b +b b b +α b b +α b b b = A B +α b +α b 1 α + 1 α 1 A A B [b,b ]+[b,b ], where α = 1 A A B B if B > 0, α = 0 if B = 0. ote that 0 α v/ + v. In articular, su α < 1, and α v/ for large. Define further c = b +α b 1 α for 0. What the above calculation shows is that H Bog = 1 [b,b A A B ]+[b,b ] + c e c 11 with e defined in. The commutators equal [b,b ] = a 0 a 0 a a Proof of Theorem 1: Lower Bound Recall that P denotes the rojection onto the constant function in L T d, and Q = 1 P. Denote the two-article multilication oerator vx 1 x by v for short. Using translation invariance and the Schwarz inequality P Q+Q PvQ Q+Q QvP Q+Q P εp Q+Q PvP Q+Q P ε 1 Q QvQ Q which follows from ositivity of v we conclude that v P PvP P +P PvQ Q+Q QvP P +1 εp Q+Q PvP Q+Q P ε 1 Q QvQ Q 13 for any ε > 0. The last term can be bounded by v0q Q. In second quantized language, this means that H is bounded from below by the restriction of a a + v0 ε > > > > 1 + v 1 εb b +b b +b b > > 1v0 ε

7 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS 7 to the -article sector. Here, we use again the definition 9 of b. From now on, we shall work with oerators on Fock sace, but it is always understood that we are only concerned with the sector of articles. ext, we observe that a a b b. Moreover, vb b v0 >, and hence H v0+hbog E ε where the Bogoliubov Hamiltonian H Bog was defined in 10 and E ε = 1 T + > > 1 v0+ v0 +ε v0 ε >. 14 In articular, using 11 and 1, H +E ε v0+ EBog + e c c, 15 where E Bog is the Bogoliubov energy defined in 1. The last term on the right side of 15 is ositive and can be droed for a lower bound on the ground state energy of H. For the choice ε = O/, the exected value of E ε in the ground state of H is bounded above by O/, as the bounds in Lemma 1 and show. This roves the desired lower bound on E 0. To obtain lower bounds on excited eigenvalues, it remains to investigate the ositive last term in 15. We do this via a unitary transformation. Let U = e X, where X = β b b b b with β 0 determined by tanhβ = α. ote that X is anti-hermitian and hence U is unitary. If a 0 and a 0 were relaced by, U would be the usual Bogoliubov transformation. Our modified U has the advantage of being article number conserving, however. The rice to ay for this modification is that U a q U can not be calculated anymore so easily. A second order Taylor exansion yields e tx a q e tx = a q t[x,a q ]+ for any t > 0. We comute t 0 [X,a q ] = β qa q a 0 t se sx [X,[X,a q ]]e sx ds

8 8 ROBERT SEIRIGER and [X,[X,a q ]] = 4βq a a 0 a 0 q 4β q β a a a a 0 a 0 +1 q For any t > 0 we thus have =: 4β q a q +J q. 16 e tx a q e tx = a q + t β qa q a 0 + Iterating this identity leads to with t 0 t se sx 4β q a q +J q e sx ds. U a q U = coshβ q a q +sinhβ q a q +K q =: d q +K q K q = In articular, we see that 1 for any λ > 0. We further have d d = c c a a 1 α 0 a 0 e sx J q e sxsinhβ q1 s ds. β q U a qa q U = d qd q +K qk q +d qk q +K qd q 1+λd qd q +1+λ 1 K qk q 17 a 0 a 0 1 α a a 1 α a 0 a 0 a 0 a 0 c c + a a > 1 α. 18 Using Schwarz s inequality, we can bound K q K q as K q K q coshβ q 1 β q 1 0 e sx J q J qe sxsinhβ q1 s β q ds. 19 To get an uer bound on J qj q, we write J q as the sum of two terms, J q 1 +J q, where J q is the second term on the right side in the first line of 16, and J q 1 = 4βq a q a 0 a 0 1. Using J q J q J q 1 J q 1 +J q J q as well as β a a q a q a q q 0β β > > 1, we obtain the bound J q J q C 1 β q > +1, C 1 = 64 su q q 0β + β. 0 To roceed, we need an uer bound on e sx > +1 e sx for 0 s 1. For this urose, let us comute [X, > ] = β q b q b q +b qb q. q 0

9 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS 9 We have [X, > ] 8 β b b q b q b q +8 q 0β b b β β q b q b q q 0 8 [ β > > 1 +1 > + > + > +1 > + > > ] 16 β With the aid of Schwarz s inequality, we obtain > [X, > +1 ] = > +1[X, > ]+[X, > ] > +1 for any η > 0. In articular, with We conclude that η > +1 η 1 [X, > ] [X, > +1 ] C > +1 C = 8 q 0 β q. t e tx > +1 e tx > +1 +C e sx > +1 e sx ds for any t > 0. Iterating this bound gives 0 e tx > +1 e tx e tc > +1. In combination, 19, 0 and yield the bound K q K q > +1 C 1 e Ccoshβ q 1 16βq. 3 By combining 17 with 18 and 3, we obtain e c c 1 1+λ U e a a U > T 1+λ 1 > +1 C 1 e C q su e 1 α e q coshβ q 1 16β q. 4 ote that q βq vq / q for large q, hence the last sum is finite. Alying this bound to 15, we have thus shown that H +Ẽε,λ v0+ EBog λ U e a a U

10 10 ROBERT SEIRIGER where Ẽ ε,λ = E ε + > T su e 1 α +1+λ 1 > +1 C 1 e C q coshβ q 1 e q 16βq, with E ε defined in 14. The desired lower bound now follows easily from the min-max rincile, and the fact that the sectrum of e a a equals e n, with n {0,1,,...} for all πz d. In fact, for any function Ψ in the sectral subsace of H corresonding to energy E E 0 +ξ, we have Ψ Ẽε,λ Ψ O ε+ ξ +ξ ε 1 +λ 1 according to Lemmas 1 and. The choice ε = O ξ/ = λ then leads to the conclusionthat the sectrum H belowan energye 0 +ξ is bounded from below by the corresonding sectrum of v0+ebog + e a a O This comletes the roof of the lower bound. ξ 3/ /. 5. Proof of Theorem 1: Uer Bound We roceed in essentially the same way as in the lower bound. In analogy to 13, we have v P PvP P +P PvQ Q+Q QvP P +1+εP Q+Q PvP Q+Q P+1+ε 1 Q QvQ Q for any ε > 0. Together with this imlies the uer bound where b b a a 1 > H v0+hbog +F ε F ε = > 1 T +ε v0 > +1+ε 1 > > 1v0. 5 From a lower bound to the commutator 1, namely we get [b,b ]+[b,b ] = a 0 a 0 a a a a 3>, H Bog E Bog 1 3> + c e c. 6

11 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS 11 Finally, to investigate the last term on the right side of this exression, we roceed as in 17 4, with the obvious modifications to get an uer bound instead of a lower bound. In relacement of 18 we use The result is d d c c 1 a a 1 α e c c 1 1 λ U + T for any λ > 0. Altogether, this shows that su α > > +1 1 α. e a a U + > > +1 e α 1 α e 1 α +λ 1> +1 C 1 e C q coshβ q 1 e q 16βq H v0+ebog λ U e a a U + F ε,λ, with F 3 ε,λ given by the sum of F ε in 5, > 1 1 E Bog from 6 and the last three terms in 7. To comlete the uer bound, we need a bound on U F ε,λ U. For this urose, we find it convenient to first bound F ε,λ by ε F ε,λ C 1 3 +λ 1 T +1 + ε+ T for an aroriate constant C 3 > 0. What remains to be shown is that 7 UT +1 U e C4 T +1 8 for some constant C 4 > 0. Given 8, we obtain UH U v0+ebog + 1 e a 1 λ a ε +C 3 e C4 1 +λ 1 T +1 + ε+ T. 9 The sectrum of the oerator on the right side of this inequality has exactly the desired form. Given an eigenvalue of e a a with value ξ, we choose ε = O ξ/ = λ to obtain v0 + EBog + ξ + Oξ 3/ / for the right side of 9. This gives the desired uer bound. It remains to rove 8. This can be done in essentially the same way as in the roof of. In fact, [X,T] = β q q b qb q +b q b q q and hence, similarly to 1, [X,T] 16 q 4 βq > +1. q

12 1 ROBERT SEIRIGER Recall that β q vq / q 4 for large q, which imlies the finiteness of the sum since v is bounded by assumtion and hence, in articular, square integrable on T d. By Schwarz, for any η > 0. In articular, for some C 4 > 0. We conclude that [X,T +1 ] ηt +1 +η 1 [X,T] [X,T +1 ] C 4 T +1 t e tx T +1 e tx T +1 +C 4 e sx T +1 e sx ds for any t > 0. Iterating this bound yields 8. This comletes the roof of the uer bound, and hence the roof of Theorem 1. Acknowledgments. It is a leasure to thank J. Fröhlich for insiring discussions and for drawing my attention towards studying the mean-field limit. Financial suort by the U.S. ational Science Foundation grant PHY is gratefully acknowledged. References [1].. Bogoliubov, On the theory of suerfluidity, J. Phys. U.S.S.R. 11, [] F. Calogero, Ground State of a One-Dimensional -Body System, J. Math. Phys. 10, Solution of the One-Dimensional -Body Problems with Quadratic and/or Inversely Quadratic Pair Potentials, J. Math. Phys. 1, [3] H.D. Cornean, J. Dereziński, P. Ziń, On the infimum of the energy-momentum sectrum of a homogeneous Bose gas, J. Math. Phys. 50, [4] L. Erdős, B. Schlein, H.-T. Yau, Ground-state energy of a low-density Bose gas: A secondorder uer bound, Phys. Rev. A 78, [5] J. Fröhlich, E. Lenzmann, Mean-Field Limit of Quantum Bose Gases and onlinear Hartree equation, Séminaire É. D. P. XVIII, [6] J. Fröhlich, A. Knowles, S. Schwarz, On the Mean-Field Limit of Bosons with Coulomb Two-Body Interaction, Commun. Math. Phys. 88, [7] M. Girardeau, Relationshi between Systems of Imenetrable Bosons and Fermions in One Dimension, J. Math. Phys. 1, [8] A. Giuliani, R. Seiringer, The Ground State Energy of the Weakly Interacting Bose Gas at High Density, J. Stat. Phys. 135, [9] E.H. Lieb, W. Liniger, Exact Analysis of an Interacting Bose Gas. I. The General Solution and the Ground State, Phys. Rev. 130, [10] E.H. Lieb, Exact Analysis of an Interacting Bose Gas. II. The Excitation Sectrum, Phys. Rev. 130, [11] E.H. Lieb, R. Seiringer, J.P. Solovej, J. Yngvason, The Mathematics of the Bose Gas and its Condensation, Oberwolfach Seminars, Vol. 34, Birkhäuser 005. Also available at arxiv:cond-mat/ [1] E.H. Lieb, J.P. Solovej, Ground State Energy of the One-Comonent Charged Bose Gas, Commun. Math. Phys. 17, Errata 5, [13] E.H. Lieb, J.P. Solovej, Ground State Energy of the Two-Comonent Charged Bose Gas, Commun. Math. Phys. 5, [14] J.P. Solovej, Uer Bounds to the Ground State Energies of the One- and Two-Comonent Charged Bose Gases, Commun. Math. Phys. 66, [15] B. Sutherland, Quantum Many-Body Problem in One Dimension: Ground State, J. Math. Phys. 1, Quantum Many-Body Problem in One Dimension: Thermodynamics, J. Math. Phys. 1, [16] H.-T. Yau, J. Yin, The Second Order Uer Bound for the Ground Energy of a Bose Gas, J. Stat. Phys. 136,

13 THE EXCITATIO SPECTRUM FOR WEAKLY ITERACTIG BOSOS 13 Deartment of Mathematics and Statistics, McGill University, Burnside Hall, 805 Sherbrooke Street West, Montreal, Quebec H3A K6, Canada address: