arxiv: v1 [math-ph] 26 Oct 2017

Transcription

1 ORM APPROXIMATIO FOR MAY-BODY QUATUM DYAMICS: FOCUSIG CASE I LOW DIMESIOS PHA THÀH AM AD MARCI APIÓRKOWSKI arxiv: v [math-ph] 26 Oct 207 Abstract. We study the norm approximation to the Schrödinger dynamics of bosons in R d (d =,2) with an interaction potential of the form dβ w( β (x y)). Here we are interested in the focusing case w 0. Assuming that there is complete Bose-Einstein condensation in the initial state, we show that in the large limit, the evolution of the condensate is effectively described by a nonlinear Schrödinger equation and the evolution of the fluctuations around the condensate is governed by a quadratic Hamiltonian, resulting from Bogoliubov approximation. Our result holds true for all β > 0 when d = and for all 0 < β < when d = 2.. Introduction Bose-Einstein condensation is a remarkable phenomenon of dilute Bose gases at very low temperatures, when many particles occupy a common single quantum state. This was predicted in 924 by Bose and Einstein [7, 20] and produced experimentally in 995 by Cornell, Wieman and Ketterle [2, 25]. Since then, there have remained fundamental questions in the rigorous understanding of the condensation and fluctuations around the condensate. The latter is essential for the emergence of superfluidity and other interesting quantum effects. From first principles of quantum mechanics, a system of identical (spinless) bosons in R d is described by a wave function in the bosonic Hilbert space H = L 2 (R d ). The evolution of the system is governed by the Schrödinger equation sym i t Ψ (t) = H Ψ (t) () where H is the Hamiltonian of the system. We will focus on the typical Hamiltonian of the form H = xj + w (x j x k ). (2) j= j<k In order to model a short-range interaction, we will take w (x) = dβ w( β x). (3) where β > 0 and w : R d R are fixed. We put the factor ( ) in front of the interaction potential in order to ensure that the interaction energy and the kinetic energy have the same order of magnitude. We can think of the initial state Ψ (0) as a ground state of a trapped system described by the Hamiltonian H V = H + V(x i ) j= Date: October 27, 207.

2 2 P.T. AM AD M. APIÓRKOWSKI with V an external trapping potential. When the trapping potential is turned off, Ψ (0) is no longer a ground state of H, and its nontrivial evolution is described by Ψ (t) in (). Although the Schrödinger equation () is linear, its complexity increases dramatically when becomes large. Therefore, for computational purposes, it is important to derive effective descriptions for collective behaviors of the quantum system. To the leading order, complete Bose-Einstein condensation means Ψ (0) u(0) in an appropriate sense. If webelievethatcondensationisstableundertheschrödingerflow, namelyψ (t) u(t), then, by formally substituting the ansatz u(t) into the Schrödinger equation (), we obtain the Hartree equation i t u(t,x) = ( x +(w u(t,.) 2 )(x) µ (t) ) u(t,x), (4) where µ (t) is a phase factor. In the large limit, the Hartree evolution can be further approximated by the (-independent) nonlinear Schrödinger equation (LS) i t ϕ(t,x) = ( x +a ϕ(t,x) 2 µ(t))ϕ(t,x). (5) ote that the meaning of the approximation Ψ (t) u(t) has to be understood properly as a convergence of one-body reduced density matrices. In particular, this approximation does not hold true in the norm topology of H, except for w = 0 (non-interacting case). The rigorous derivation of the Hartree and nonlinear Schrödinger equation has been the subject of a vast literature, initiated by Hepp [3], Ginibre and Velo [2] and Spohn [46]. The results for 0 < β (dilute regime) has been studied by Erdös, Schlein and Yau [7, 8, 9] for d = 3, Adami, Golse and Teta [] for d = and Kirkpatrick, Schlein and Staffilani [26] for d = 2 (see also [24]). All these works deal with the defocusing case w 0. In the focusing case w 0, the LS is only well-posed globally for d 2, and it is natural to restrict to these low dimensions. The derivation of the focusing LS has been achieved very recently by Chen and Hölmer [4, 3] by means of the BBGKY approach, and then by Jeblick and Pickl [23] by another approach. In these works, it is crucial to consider the -body Hamiltonian H V with a trapping potential like V(x) = x s, and restrict to 0 < β < when d = [4] and 0 < β < (s+)/(s +2) when d = 2 [3, 23]. The presence of the trap and the restriction on β allow to use the stability of the second kind H V C by Lewin, am and Rougerie [35, 36]. This stability is important to control the (negative) interaction potential by the kinetic operator. In the present paper, we are interested in the norm approximation, which is much more precise than the convergence of density matrices. It requires to understand not only the condensate but also the fluctuations around the condensate. Motivated by rigorous results on ground states of trapped systems [45, 22, 34, 6, 43, 5], we will assume that the initial datum satisfies the norm approximation Ψ (0) u(0) ( n) s ψ n (0). n=0 Bogoliubov s approximation [6] suggests that Ψ (t) u(t) ( n) s ψ n (t) (6) n=0 where (ψ n (t)) n=0 is governed by an effective Hamiltonian in Fock space which is quadratic in creation and annihilation operators.

3 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 3 The norm approximation (6) has been established by Lewin, am and Schlein [33] for β = 0 (see [39] for a similar result). See [40, 4, 8] for extensions to 0 < β < with w 0 in d = 3 dimensions. The goal of the present paper is to extend the norm approximation (6) to lower dimensions, in both defocusing and focusing cases. The focusing case is more interesting as we have to rule out the instability of the system. Our method can treat all β > 0 for d = and 0 < β < for d = 2. So in particular, we can recover and improve the leading order results in [4, 3, 23]. More precisely, we can deal with a larger (and much more natural) range of β and do not need to add a trapping potential (although our method works with the presence of an external potential as well). To achieve our result we will use a new localization technique, which allows us to go beyond the stability regime established in [35, 36]. Moreover, our main result not only deals with the condensate, but also provides a detailed description for the fluctuations. An analogue of (6) related to the fluctuations around coherent states in Fock space has attracted many studies [29, 27, 32, 4, 28]. All these works concentrate on the defocusing case in 3D. It is straightforward to apply our method to the Fock space setting to treat the focusing case in one- and two-dimensions. Finally, let us mention that in 3D, the focusing LS may blow-up at finite time and it is interesting to derive effective equations locally in time. We refer to [44, 40, 5] for relevant results in this direction. Our method should be also useful for this problem. We hope to come back this issue in the future. The precise statement of our result is given in the next section. Acknowledgments. We thank Christian Brennecke, Mathieu Lewin, icolas Rougerie, Benjamin Schlein and Robert Seiringer for helpful discussions. The support of the ational Science Centre (C) project r. 206/2/D/ST/02430 is gratefully acknowledged (M). 2. Main result First, we start with our assumption on the interaction potential. Recall from (3) that we are choosing w (x) = dβ w( β x). We will always assume that w L (R d ), w(x) = w( x) R. (7) When d =, there is no further assumption is needed (indeed, our proof can be even extended to a delta interaction ±δ 0 ). In d = 2, we need an additional assumption that w L (R d ), w < a, w = max( w,0). (8) R 2 Here a is the optimal constant in the Gagliardo irenberg inequality ( )( ) f 2 f 2 a f 4, f H (R 2 ). (9) R 2 R d 2 R 2 Indeed, it is well-known [47, 48, 30, 38] that a := Q 2 L 2 where Q H (R 2 ) is the unique positive radial solution to Q+Q Q 3 = 0. (0) The condition (8) is essential for the stability of the 2D system; see [35] for a detailed discussion. This condition has been used in the derivation of the nonlinear Schrödinger equation for ground states of trapped system [35, 36], as well as for the dynamics [3, 23].

4 4 P.T. AM AD M. APIÓRKOWSKI Hartree equation. In our paper, the condensate u(t) is governed by the Hartree equation (4) with the phase µ (t) = u(t,x) 2 w (x y) u(t,y) 2 dxdy 2 R d R d which is chosen to ensure an energy compatibility (see [33] for further explanations). Here we have written u(t) = u(t,.) and ignored the -dependencein the notation for simplicity. The well-posedness of the Hartree equation will be discussed in Section 3. It is easy to see that when, u(t) convergence to the solution of the cubic nonlinear Schrödinger equation (5). Bogoliubov equation. To describe the fluctuations around the condensate, it is convenient to turn to the grand-canonical setting of Fock space n F(H) = H n, H n = H, H = L 2 (R d ) n=0 where the number of particles can vary (and indeed, as we will soon see, the number of excitations will not be fixed). On Fock space, we define the creation and annihilation operators a (f), a(f), with f H, by n+ (a (f)ψ)(x,...,x n+ ) = f(x j )Ψ(x,...,x j,x j+,...,x n+ ), n+ (a(f)ψ)(x,...,x n ) = n j= sym f(x n )Ψ(x,...,x n )dx n, Ψ H n, n. These operators satisfy the canonical commutation relations (CCR) [a(f),a(g)] = [a (f),a (g)] = 0, [a(f),a (g)] = f,g, f,g H. We can also define the operator-valued distributions a x and a x, with x R d, by a (f) = f(x)a xdx, a(f) = f(x)a x dx, f H, R d R d which satisfy [a x,a y ] = [a x,a y ] = 0, [a x,a y ] = δ(x y), x,y Rd. The Hamiltonian H can be extended to Fock space as H = a x ( x)a x dx+ w (x y)a x 2() a y a xa y dxdy. For every one-body operator h, we will use the short hand notation n dγ(h) := a x h xa x dx = 0 h j. In particular, = dγ() is called the number operator. We look for the norm approximation of the form Ψ (t) u(t) ( n) s ψ n (t) := n=0 n=0 n=0 j= (a (u(t))) n ( n)! ψ n (t). Here the particles outside of the condensate are described by a unit vector Φ(t) = (ψ n (t)) n=0

5 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 5 in the excited Fock space n F + (t) = H n +, Hn + = H + (t), H + (t) = {u(t)} L 2 (R d ). n=0 sym As explained in [33], Bogoliubov approximation suggests that the vector Φ(t) solves the equation i t Φ(t) = H(t)Φ(t) () where H(t) is a (-dependent) quadratic Hamiltonian on the Fock space F of the form H(t) = dγ(h(t))+ ) (K 2 (t,x,y)a 2 xa y +K 2 (t,x,y)a x a y dx dy. Here R d R d h(t) = + u(t, ) 2 w µ (t)+q(t) K (t)q(t), K 2 (t) = Q(t) Q(t) K 2 (t) where the kernel of the operator K (t) and the 2-body function K 2 (t) H 2 are K (t,x,y) = u(t,x)w (x y)u(t,y), K2 (t,x,y) = u(t,x)w (x y)u(t,y). The well-posedness of () will be revisited in Section 4. ote that Φ(t) belongs to F + (t), which is not obviously seen from the equation (). Main result. ow we are ready to state our main result. Theorem (orm approximation). Let β > 0 when d = and 0 < β < when d = 2. Assume that the interaction potential w satisfies (7) (8). Let u(t) satisfy the Hartree equation (4) with u(0) H d+2 (R d ) C. Let Φ(t) = (ϕ n (t)) n=0 satisfy the Bogoliubov equation () with Φ(0),dΓ( )Φ(0) C. Consider the -body Schrödinger evolution Ψ (t) in () with the initial state Ψ (0) = u(0) ( n) s ψ n (0). n=0 Take α = /2 when d = and 0 < α < ( β)/3 arbitrarily when d = 2. Then for all t > 0, there exists a constant C t > 0 independent of such that for all large, Ψ (t) u(t) ( n) s ψ n (t) 2 C t α. (2) H n=0 Remark 2. We have some immediate remarks concerning the result in Theorem. () The initial state Ψ (0) is not normalized, but its norm converges to very fast when (as we will explain in the proof). Hence, we ignore a trivial normalization in the statement of Theorem. (2) Our approach is quantitative and our result applies equally well for trapped systems. Moreover, the initial values u(0) and Φ(0) can be chosen -dependently, as soon as u(0) H d+2 (R d ) and Φ(0),dΓ( )Φ(0) grow slowly enough. We do not include these extensions to simplify the representation. (3) When d = 2, our result holds true if (8) is replaced by the weaker condition from [35] u(x) 2 u(y) 2 w(x y)dxdy inf R 2 R 2 u H (R 2 ) 2 u 2 L 2 (R 2 ) u 2 >. L 2 (R 2 ) The latter condition is enough to ensure the well-posedness of the Hartree equation.

6 6 P.T. AM AD M. APIÓRKOWSKI (4) A simplified formulation for Bogoliubov equation () can be given in terms of density matrices. Recall that for a vector Φ in Fock space, its one-body density matrices γ Ψ : H H and α Ψ : H H H can be defined by f,γ Ψ g = Ψ,a (g)a(f)ψ, f,α Ψ g = Ψ,a(g)a(f)Ψ, f,g H. As explained in [40], if Φ(t) solves (), then (γ(t),α(t)) = (γ Φ(t),α Φ(t) ) is the unique solution to the system i t γ = hγ γh+k 2 α α K2, i t α = hα+αh T +K 2 +K 2 γ T +γk 2, (3) γ(0) = γ Φ(0), α(0) = α Φ(0). ote that (3) is similar (but not identical) to the equations studied in the Fock space setting [27, 32, 3]. Reversely, if Φ(0) is a quasi-free state, then Φ(t) is a quasi-free state for all t > 0, and in this case the equation (3) is equivalent to the Bogoliubov equation (). As we have mentioned in the introduction, the norm approximation is much more precise than the convergence of density matrices. Recall that the one-body density matrix of a -body wave function Ψ is obtained by taking the partial trace γ () Ψ = Tr 2 Ψ Ψ. Equivalently, γ () Ψ is a trace class operator on L 2 (R d ) with kernel γ () Ψ (x,y) = Ψ (x,x 2,...,x )Ψ (y,x 2,...,x )dx 2...dx. The following result is a direct consequence of Theorem. Corollary 3 (Convergence of reduced density). Under the same conditions in Theorem, we have the convergence in trace class lim γ() Ψ (t) = ϕ(t) ϕ(t), t > 0, where ϕ(t) is the solution to the following cubic nonlinear Schrödinger equation, i t ϕ(t,x) = ( x +a ϕ(t,x) 2 )ϕ(t,x), ϕ(0,x) = u(0,x), a = w. (4) ote that (4) is different from the equation (5) in the introduction because there is no phase µ(t). However, the phase plays no role when we consider the projection ϕ(t) ϕ(t). Corollary 3 recovers and improves the existing results on the leading order: [26, 2] for defocusing 2D with 0 < β < ; [4] for focusing D with 0 < β < ; [3, 23] for focusing 2D with 0 < β < (s + )/(s + 2) with a trapping potential like x s. Here our range of β is larger, i.e. β > 0 for d = and 0 < β < for d = 2, which is optimal to some extent. Moreover, we do not have to restrict to trapped systems, which is consistent with our interpretation that Ψ (t) is the evolution of a ground state of a trapped system when the trap is turned off (but our proof works equally well with trapped systems). To improve the range of β and to remove the trap, we will not rely on the stability of the second kind H C in [35, 36]. Indeed, thanks to a new localization method on the number of particles, we will only need a weaker version of the stability on the sector of very few particles, which is much easier to achieve. This weaker stability is enough for our purpose because the fluctuations around the condensate involve only very few particles (most of particles are already in the condensate mode). Organization of the paper. The paper is organized as follows. We will always focus on the more difficult case d = 2, and only explain the necessary changes for d = at the end. We will revise the well-posedness of the Hartree equation (4) in Section 3 and the

7 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 7 Bogoliubov equation () Section 4. In Section 5, we reformulate the Bogoliubov approximation using a unitary transformation from H to a truncated Fock space, following ideas in [34, 33]. Then we provide several estimates which are useful to implement Bogoliubov s approximation. In Section 6 we explain the localization method. Then we prove Theorem is presented in Section 7, for d = 2. All the changes needed to prove Theorem for d = are explained in Section Hartree dynamics In this section, we discuss the well-posedness and provide various estimates for Hartree equation (4). Our proof will require bounds on the solutions of the Hartree equation (4). Under our assumptions on the nonlinearity, it is well-known that the equation is globally well-posed in H (R d ) (see, e.g., [, Cor. 6..2]). However, since the potential w depends on, it is not clear if the norm u(t,.) H is bounded uniformly in. The same question applies to other norms we will use in the proof. We are going to prove that it is indeed the case. Here we will consider the case d = 2 in detail. Remarks about the corresponding results for d = will be given in Section 8. Lemma 4. Let d = 2. Assume w L (R 2 ) L (R 2 ) and w < a. For every u 0 H 4 (R 2 ) with u 0 L 2 =, equation (5) has a unique solution u(t, ) in H 4 (R 2 ) and we have for all t > 0 the bounds u(t, ) H (R 2 ) C, u(t, ) H 2 (R 2 ) Cexp(Cexp(Ct)), u(t, ) L (R 2 ) Cexp(Ct), t u(t, ) L (R 2 ) Cexp(exp(Cexp(Ct))). for a constant C > 0 independent of t and. Proof. For convenience, we will work with the equation i t u = u+(w u 2 )u, u(0,x) = u 0 (x), x R d. (5) To go from equation (5) to equation (4) it is enough consider a gauge transformation u e i t 0 µ (s)ds u with µ (t) = u(t,x) 2 w (x y) u(t,y) 2 dxdy. 2 R d R d The -dependence in the desired bounds is not affected by this change. Step. First, the bound in H follows from the energy conservation u(t,.) 2 L + u(t,x) 2 w 2 (x y) u(t,y) 2 dxdy = C, 2 the simple estimate u(t,x) 2 w (x y) u(t,y) 2 dxdy u(t,x) 2 [w (x y)] u(t,y) 2 dxdy u(t,x) 4 + u(t,y) 4 [w (x y)] dxdy 2 = u(t,.) 4 L w 4 and the bounds (8) and (9).

8 8 P.T. AM AD M. APIÓRKOWSKI with Step 2. ext, to prove the H 2 bound, we use Duhamel s formula It follows that We compute [ ] G = (w u 2 )u u(t,x) = e it u 0 (x)+ t 0 e i (t s) G(s,x)ds G(s, ) = (w u(s, ) 2 )u(s, ). u(t, ) L 2 u 0 L 2 + = (w u 2 ) u+ For the first term we easily find 2 i= t 0 G(s, ) L 2 ds. [ ] [ ] 2 w ( xi ( u 2 )) xi u+ w ( ( u 2 )) u. (w u 2 ) u L 2 C w L u(s, ) 2 L u L 2. Since xi ( u 2 ) = ū xi u+u xi ū in the second term of (6) we need to bound w (y)ū( y) xi ū( y)dy xi u( ) L 2 C u L w (y) ( xi ū( y) 2 + xi u( ) 2) L 2 C u L w L u 2 L 2. To bound the last term we use the the two dimensional Sobolev inequality u 2 4 C( u 2 H + u 2 L 2 ). Treating other terms in (6) in a similar way we obtain Thus we deduce that G(s, ) L 2 C(+ u(s, ) L 2)(+ u(s, ) 2 L ). u(t, ) H 2 C +C t 0 u(s, ) H 2(+ u(s, ) 2 L )ds. We can now use the two dimensional Brezis Gallouet Wainger (or log-sobolev) inequality [9, 0] (6) v(s, ) 2 L C(+log(+ v(s, ) H2)) (7) which holds true for functions v with v H =. By the assumption u 0 L 2 = and mass preservation it follows that u(t, ) H. Using this and u(t, ) H C from step we deduce from the Brezis Gallouet Waigner inequality for v = u/ u H that Denote Then u(t, ) H 2 C +C F(t) = C +C t 0 t 0 u(s, ) H 2(+log(+ u(s, ) H 2))ds. u(s, ) H 2(+log(+ u(s, ) H 2))ds. F (t) = C u(t, ) H 2(+log(+ u(t, ) H 2)) CF(t)[+log(+F(t))] which implies that d log(+log(+f(t))) C. dt

9 This allows us to conclude that FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 9 u(t,.) H 2 exp(cexp(ct)) which using (7) immediately implies the bound on u(t, ) L. Step 3. Let us now prove the last bound. We will again use use the Brezis Gallouet Waigner inequality, this time for t u. It reads where t u(t, ) 2 L C 0(t)(+log(+ t u(t, ) H 2)) (8) C 0 (t) = max{ t u(t, ) 2 H,}. Since v H v H 2 it is clear that we need to obtain a bound on t u(t, ) H 2. From the Hartree (5) equation we get The norm equivalence thus implies t u(t, ) H 2 u(t, ) H 2 + t 0 G(s, ) H 2. c ( f 2 L 2 + f 2 L 2 ) /2 f H 2 C ( f 2 L 2 + f 2 L 2 ) /2 t u(t, ) H 2 C ( u L ) t u L 2 + C( G(s, ) L 2 + G(s, ) L 2)ds. Bounds obtained in the previous steps thus lead to 0 t u(t, ) H 2 C(+t)exp(Cexp(Ct))+C u(t, ) H 4. (9) By Duhamel s formula, for any integer k, we have u(t, ) H k u 0 ( ) H k + t 0 G(s, ) H k ds. (20) We shall first get a bound on u(t, ) H 3 and therefore we first look at G(t, ) H 3. To this end we notice that by the multi-index Leibniz formula we have G(t, ) H 3 = D l G L 2 (2) l 3 k l m k l 3 ( 3 k )( ) k m w (y) m u(x y) k m ū(x y)dy l k u(x) L 2. We will derive a Grönwall type inequality. Since our goal is to obtain a bound that is independent of, in our bounds we will need to extract the L norm of w ; otherwise the -dependence will appear. Thus, to do the dy integration we need to use the L bounds on the derivatives of u. By (7) we have where i u(t, ) 2 L C 0(t)(+log(+ i u(t, ) H 2)) C 0 (t)(+log(+ u(t, ) H 2+ i )) (22) C 0 (t) = Cmax{ i u(t, ) 2 H,} Cmax{ u(t, ) 2 H + i,}. Let us now look at different terms in the sum in (2). The ingredients for the case l = 0 can be trivially bounded by For l = we obtain the bound l = 0 term C u 2 L w u L 2 Cexp(Ct). (23) l = term C u 2 L w u H Cexp(Ct). (24)

10 0 P.T. AM AD M. APIÓRKOWSKI For the l = 2 term we have two possibilities. Either both derivatives hit one function or they distribute among two functions. This leads to l = 2 term C u 2 L w u H 2 +C u L i u L w u H C(+exp(Cexp(Ct)))(+log(+ u H 3)) (25) where we used (22). It remains to bound the term corresponding to l = 3. We again notice that if all derivatives hit one function then we get a bound Cexp(Ct) u H 3. In the case when two derivatives hit one function and and one derivative hits one of the others we get a bound as in (25). The same bound works for the case when one derivative hits each function. Altogether we arrive at the bound for the l = 3 term of the form l = 3 term C(+exp(Cexp(Ct)))(+ u H 3). (26) Collecting the bounds (23)-(26) and inserting them into (2) yields G(t, ) H 3 C(+exp(Cexp(Ct)))(+ u H 3). Inserting this into (20) implies by Grönwall inequality that u(t, ) H 3 Cexp(Cexp(Cexp(Ct))). (27) With this result we are finally able to proceed to the last estimate, that is the bound on u(t, ) H 4. In this case the analysis proceeds in the same way. We again write out the H 4 counterpart of (2) with l now being not bigger than four. The analysis above shows that all terms in the sum which correspond to l 3 can be bounded in terms of Cexp(Cexp(Cexp(Ct))). Let us now look at different terms corresponding to l = 4. If all derivatives hit one function then obviously we bound it by Cexp(Ct) u H 4. If two derivatives hit one function and two another one, then using (22) with i = 2 we obtain a bound of the form Cexp(Cexp(Cexp(Ct)))( + log( + u H 4)). Other terms can be bounded independently of u H 4 by Cexp(Cexp(Cexp(Ct))). Again, using Grönwall, we arrive at he bound u(t, ) H 4 C(exp(Cexp(Cexp(Cexp(Ct)))). (28) Inserting (28) into (9) and this into (8), we see that finally t u(t, ) 2 L Cexp(exp(Cexp(Ct))). This ends the proof. From now on, we will often omit the explicit time-dependence on the constant and replace it by a general notation C t, for simplicity. Also, we will focus on d = 2. The main result of this section is 4. Bogoliubov dynamics Lemma 5. Assume that Φ(0) satisfies Φ(0), dγ( )Φ(0) C. Then the Bogoliubov equation () has a unique global solution Φ C([0, ),F(H)) L loc ((0, ),Q(dΓ( ))). Moreover, Φ(t),dΓ( )Φ(t) Ct,ε ε. The well-posedness in Lemma 5 follows from [33, Theorem 7]. The new result is the kinetic estimate. This follows from the following quadratic form estimates on F(H).

11 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS Lemma 6. For every ε > 0 and η > 0, ( ) ± H(t)+dΓ( ) ηdγ( )+C t,ε (+η )( + ε ), ± t H(t) ηdγ( )+C t,ε (+η )( + ε ), ±i[h(t),] ηdγ( )+C t,ε (+η )( + ε ), Let us assume Lemma 6 for the moment and provide Proof of Lemma 5. We will use Gronwall s argument. By Lemma 6, we have A(t) := H(t)+C t,ε ( + ε ) 2 dγ( ), for C t,ε sufficiently large and, by the equation for Φ(t), d Φ(t),A(t)Φ(t) = Φ(t), t A(t)Φ(t) + Φ(t),i[H(t),A(t)]Φ(t) dt = Φ(t), t (H(t)+ t C t,ε ( + ε ))Φ(t) + Φ(t),i[H(t),]Φ(t) Thus C t,ε Φ(t),A(t)Φ(t). Φ(t),A(t)Φ(t) e Ct,ε Φ(0),A(0)Φ(0). The left side is bounded from below by 2 Φ(t),dΓ( )Φ(t). The right side can be bounded by using A(0) C ε (dγ( )+ ε ). We thus conclude that Φ(t),dΓ( )Φ(t) Cε e Ct,ε( Φ(0),dΓ( )Φ(0) ) + ε). This ends the proof. ow we turn to the proof of Lemma 6. We will need two preliminary results. The first is a lower bound on general quadratic Hamiltonians, taken from [42, Lemma 9]. Lemma 7. Let H > 0 be a self-adjoint operator on H. Let K : H H H be an operator with kernel K(x,y) H 2. Assume that KH K H and that H /2 K is Hilbert-Schmidt. Then dγ(h)+ ( K(x,y)a 2 xa y +K(x,y)a x a y )dxdy 2 H /2 K 2 HS. The second is the following kernel estimate. Lemma 8. For all ε > 0 we have ( x ) /2 K 2 (t,, ) 2 L 2 + ( x ) /2 t K 2 (t,, ) 2 L 2 C t,ε ε. Proof of Lemma 8. We will focus on t K 2 (t) as K 2 (t) can be treated similarly. By the definition we have K 2 (t,, ) = Q(t) Q(t) K 2 (t,, ), K2 (t,x,y) = u(t,x)w (x y)u(t,y). t K 2 (t) = t Q(t) Q(t) K 2 (t)+q(t) t Q(t) K 2 (t)+q(t) Q(t) t K2 (t). First, letuscompare t K 2 (t)with t K2 (t). Since t Q(t) = t u(t) u(t) u(t) t u(t), we have t Q(t) Q(t) K 2 (t,, ) L 2 ( t Q(t) ) K 2 (t,, ) L 2 ( t u(t) u(t) ) K 2 (t,, ) L 2 + ( u(t) t u(t) ) K 2 (t,, ) L 2

12 2 P.T. AM AD M. APIÓRKOWSKI For the first term, it is straightforward to see that ( t u u ) K 2 (t,, ) 2 L 2 = u(t,z)u(t,z)w (z y)u(t,y)dz u(t, ) 4 L w 2 L u(t, ) 2 L 2 t u(t, ) 2 L 2 C t. The second term is also bounded by the same way. Thus we find that By the same argument, we also obtain and In summary, Since ( x ) /2 on L 2, we obtain t Q(t) Q(t) K 2 (t,, ) L 2 C t. Q(t) t Q(t) K 2 (t,, ) L 2 C t ( Q(t) Q(t)) t K2 (t,, ) L 2 C t. t K 2 (t,, ) t K2 (t,, ) L 2 C t. 2 t u(t,x) 2 dxdy (29) ( x ) /2 t K 2 (t,, ) ( x ) /2 t K2 (t,, ) L 2 C t. It remains to bound ( x ) /2 t K2 (t,, ). We have t K2 (t,, ) = t u(t,x)w (x y)u(t,y)+u(t,x)w (x y) t u(t,y) and it suffices to consider f(t,x,y) = t u(t,x)w (x y)u(t,y) (the other term is similar). Inspired by the idea in [27], we compute the Fourier transform: f(t,p,q) = u(t,x)w (x y)( t u)(t,y)e 2πi(p x+q y) dxdy = u(t,y +z)w (z)( t u)(t,y)e 2πi(p (y+z)+q y) dzdy = w (z)(u z t u)(t,p+q)e 2πip z dz where u z (t, ) := u(t,z + ). By the Cauchy-Schwarz inequality, f(t,p,q) 2 w L w (z) (u z t u)(t,p+q) 2 dz and by Lemma 4, (u z t u)(t,p+q) 2 dq = (u z t u)(t, ) 2 L C 2 t. Thus by Plancherel s Theorem, for all ε > 0 we have ( x ) /2 ε f(t,, ) 2 L = (+ 2πp 2 ) 2ε 2 f(t,p,q) dpdq 2 w L (+ 2πp 2 ) 2ε w (z) (u z t u)(t,p+q) 2 dpdqdz C t,ε. Here we have used the facts that w L = w L and (+ 2πp 2 ) 2ε dp C ε <. Moreover, by Lemma 4 we have the simple estimate f(t,, ) 2 L 2 C u(t, ) 2 L tu(t, ) 2 L 2 w 2 L 2 C t 2β.

13 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 3 By interpolation, we thus obtain ( x ) /2 ε f(t,, ) 2 L C 2 t,ε ε, ε > 0. This ends the proof. ow we are ready to give Proof of Lemma 6. We have H(t)+dΓ( ) = dγ(h+ )+ 2 ( K 2 (t,x,y)a xa y +K 2 (t,x,y)a x a y )dxdy. By Lemma 4, it is straightforward to see that h+ = u(t, ) 2 w µ (t)+q(t) K (t)q(t) C t, and hence ±dγ(h+ ) C t. ow we consider the paring term. ote that K 2 C t. We apply Lemma 7 with H = η( )+η K 2 2, η > 0, K = ±K 2 and then use Lemma 8. We get ηdγ( )+η ± ( K 2 (t,x,y)a x 2 a y +K 2(t,x,y)a x a y )dxdy Thus 2 η ( ) /2 K 2 2 HS C t,ε η ε, ε > 0. ( ) ± H(t)+dΓ( ) ηdγ( )+C t (+η ) +C t,ε η ε for all η > 0,ε > 0. The bound for t H(t) is obtained by the same way and we omit the details. Moreover, it is straightforward to see that ( i[h(t),] = ik 2 (t,x,y)a xa y +ik 2 (t,x,y)a x a y )dxdy and the bound for i[h(t),] also follows from the above argument. This completes the proof. 5. Bogoliubov s approximation As in [34, Sec. 2.3], for every normalized vector Ψ H we can write uniquely as Ψ = u(t) ( n) s ψ n = n=0 n=0 with ψ n H + (t) n. This gives rise the unitary operator Thus (a (u(t))) n ( n)! ψ n U (t) : H F+ (t) := H + (t) n n=0 Ψ ψ 0 ψ ψ. Φ (t) := U (t)ψ (t) describes the fluctuations around the condensate u(t). As proved in [33], Φ (t) belongs to F+ (t) and solves the equation { i t Φ (t) = G (t)φ (t), Φ (0) = Φ(0). (30)

14 4 P.T. AM AD M. APIÓRKOWSKI Here m =( m) is the projection onto the truncated Fock space and with F m = C H H m G (t) = [ ] H(t)+E (t) E (t) = 2 4 (R j +Rj), (3) R 0 = R 0 = dγ(q(t)[w u(t) 2 + K (t) µ (t)]q(t)), R = 2 a(q(t)[w u(t) 2 ]u(t)), ( ) ( )( ) R 2 = K 2 (t,x,y)a xa ydxdy, R 3 = ( Q(t)w Q(t) Q(t))(x,y;x,y ) R 4 = R4 = 2() j=0 u(t,x)a ya x a y dxdydx dy, (Q(t) Q(t)w Q(t) Q(t))(x,y;x,y ) a xa ya x a y dxdydx dy. Here, in R 0 and R we write w for the function w (x), while in R 3 and R 4 we write w for the two-body multiplication operator w (x y). Bogoliubov s approximation suggests that the error term E (t) should be small. In the following we will justify this on the sectors of few particles. Lemma 9 (Bogoliubov approximation). For every m and ε > 0 we have the quadratic form estimate m ± m E (t) m C t,ε ε dγ( ) (32) m ± m t E (t) m C t,ε ε dγ( ) (33) m ± m i[e (t),] m C t,ε ε dγ( ) (34) Proof. We will use the decomposition (3). We will denote by Φ an arbitrary normalized vector in F+ m. ote that (32) and (33) follow from m Φ,R j Φ C t,ε ε Φ,dΓ( )Φ, 0 j 4 (35) and m Φ, t R j Φ C t,ε ε Φ,dΓ( )Φ, 0 j 4 (36) respectively. Moreover, since [R 0,] = [R 4,] = 0, [R,] = R, [R 2,] = 2R 2, [R 3,] = R 3, the commutator bound (34) also follows from (35). ow we prove (35)-(36) term by term.

15 j = 0 Recall that From the operator bounds FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 5 ( R 0 = dγ Q(t)[w u(t) 2 + K ) (t) µ (t)]q(t). we have Consequently, Similarly, Q(t)[w u(t) 2 + K (t) µ (t)]q(t) C t ( +) ±R 0 C t. ± m R 0 m C tm. ± t R 0 C t( +), ± m t R 0 m C tm. j = For every Φ F m + Φ,R Φ = 2 (t), by the Cauchy-Schwarz inequality we have Φ, ( ) a Q(t)[w u(t) 2 ]u(t) Φ 2 ( ) Φ a Q(t)[w u(t) 2 ]u(t) Φ. Using the elementary inequality a (v)a(v) v 2 L and 2 Q(t)[w u(t) 2 ]u(t) [w u(t) 2 ]u(t) L 2 L 2 C t we obtain 2 ( ) Φ a Q(t)[w u(t) 2 δ ]u(t) C () 2 Φ 2 +C t δ /2 Φ 2 C mδ Φ,Φ +C tδ Φ,Φ which choosing δ = /m yields the result m Φ,R Φ C t Φ,Φ. (37) Similarly, using t (Q(t)[w u(t) 2 ]u(t)) L 2 C t we get j = 2 We have Φ,R 2 Φ = = m Φ, t R Φ C t Φ,Φ. (38) ( K 2 (t,x,y) Φ,a ( )( ) x a y w (x y)u(t,x)u(t,y) Φ,a xa y ( ( )( ) ) Φ dxdy ) Φ dxdy.

16 6 P.T. AM AD M. APIÓRKOWSKI Here we have replaced K 2 (t) = Q(t) Q(t) K 2 (t) by K 2 (t), namely ignored the projection projection Q(t), because Φ belongs to the excited Fock space F + (t). By the Cauchy- Schwarz inequality, Φ,R 2 Φ w (x y). u(t,x). u(t,y). ( +) /2 a x a y Φ dxdy ( ( +)/2( ( ( )( ) ) /2 w (x y) u(t,x) 2 u(t,y) 2 dxdy ) Φ ) /2 w (x y) ( +) /2 a x a y Φ 2 dxdy ( )( ) ) Φ. ( +)/2( On the truncated Fock space F+ m, we have the operator inequality ( ( +) /2( ( )( ) ) ) 2 2m ( +). Moreover, from Sobolev s embedding in R 2 we have the operator estimate on L 2 ((R 2 ) 2 ) Consequently, w (x y) C ε ε ( x ), ε > 0. (39) w (x y) a x a y ( +) a x a y dxdy = ( +3) a xa y w (x y) a x a y dxdy C ε ε ( +3) a xa y( ) x a x a y dxdy Thus, since we obtain C ε ε dγ( ). ( w (x y) u(t,x) 2 u(t,y) 2 dxdy) /2 C t Φ,R 2 Φ C t δc ε ε Φ,dΓ( )Φ +C t δ m Φ,( +)Φ which choosing δ = m/ finally leads to ext, we consider Φ, t R 2 Φ = Φ,R 2 Φ C ε ε m Φ,dΓ( )Φ. [ ( t Q(t) + t Q(t)) K 2 (t,x,y)+ t K2 (t,x,y)] ( Φ,a xa ( )( ) ) y Φ dxdy. Here we have used the decomposition t K 2 (t) = t Q(t) Q(t) K 2 (t)+q(t) t Q(t) K 2 (t)+q(t) Q(t) t K2 (t)

17 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 7 and omitted the projection Q(t) again using Φ F+ m (t). The term involving t K 2 (t,x,y) can be bounded as above. The term involving ( t Q(t) ) K 2 (t,x,y) is bounded as ( t Q(t) ) K2 ( (t,x,y) Φ,a ( )( ) ) x a y Φ dxdy ( ( t Q(t) ) K2 ) (t,x,y) 2 /2 ( ) /2 dxdy ( +) /2 a x a y Φ 2 dxdy ( ) ( )( ) ( +)/2 Φ C tm Φ,( +)Φ. Here we used (29). The term involving t Q(t) can be bounded similarly. Thus Φ, t R 2 Φ C t,ε ε m Φ,dΓ( )Φ. Since m the desired estimate follows. j = 3 By the Cauchy-Schwarz inequality and Sobolev s inequality (39) we have Φ,R 3 Φ = w (x y)u(t,x) Φ, a ya y a x Φ dxdy w (x y) u(t,x) a y Φ ay a x Φ dxdy u(t, ) L ( ( ) /2 w (x y) a x a y Φ 2 dxdy ) /2 w (x y) a y Φ 2 dxdy ext, C t Φ,C ε ε dγ( )Φ /2 Φ,Φ /2 C t,ε ε m /2 /2 Φ,dΓ( )Φ. (40) Φ, t R 3 Φ = [( Q(t)w Q(t) Q(t) ) (x,y;x,y ) t u(t,x) ( ( + t Q(t)w Q(t) Q(t) )) ] (x,y;x,y )u(t,x) Φ, a ya x a y Φ dxdydx dy. The term involving t u(t,x) can be estimated similarly to (40). For the other term, we decompose t ( Q(t)w Q(t) Q(t) ) = t Q(t)w Q(t) Q(t)+ Q(t)w t Q(t) Q(t) We will use the kernel estimate where + Q(t)w Q(t) t Q(t) ( t Q(t))(z;z ) = t u(t,z)u(t,z )+u(t,z) t u(t,z ) q(z)q(z ) (4) q(t,z) := u(t,z) + t u(t,z), q(t, ) L 2 C t, q(t, ) L C t.

18 8 P.T. AM AD M. APIÓRKOWSKI For the first term involving t Q(t)w Q(t) Q(t), we can estimate ( t Q(t)w Q(t) Q(t))(x,y;x,y )u(t,x) Φ, a ya x a y Φ dxdydx dy = ( t Q(t))(y;y )w (x y )δ(x x )u(t,x) Φ, a ya x a y Φ dxdydx dy q(t,y) q(t,y ) w (x y ) q(t,x) a y Φ ax a y Φ dxdydy q(t, ) ( ) /2 ( L q(t,y) 2 /2 dy a y Φ dy) 2 ( ) /2 ( ) /2 w (x y ) q(t,x) 2 dxdy w (x y ) a x a y Φ 2 dxdy C t Φ,Φ /2 Φ,C ε ε dγ( )Φ /2 C t,ε ε m /2 Φ,dΓ( )Φ. /2 For the term involving Q(t)w t Q(t) Q(t), we have ( Q(t)w t Q(t) Q(t))(x,y;x,y )u(t,x) Φ, a ya x a y Φ dxdydx dy = w (x y)( t Q(t))(x,x )δ(y y )u(t,x) Φ, a y a x a y Φ dxdydx dy w (x y) q(t,x) q(t,x ) a y Φ ax a y Φ dxdydx ( /2 q(t, ) 2 L w L a x a y Φ 2 dx dy) ( ) q(t,x ) 2 /2 a y Φ 2 dx dy C Φ,2 Φ /2 Φ,Φ /2 Cm/2 Φ,Φ. /2 The term involving Q(t)w Q(t) t Q(t) is bounded similarly as above. Thus j = 4 By Sobolev s inequality (39) Φ, t R 3 Φ C t,ε ε m /2 Φ,dΓ( )Φ. /2 R 4 C ε ε dγ( ).

19 Therefore, ext, we consider Φ, t R 4 Φ = FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 9 ± m R 4 m C ε ε m dγ( ). 2() R ) t (Q(t) Q(t)w Q(t) Q(t) (x,y;x,y ) Φ,a xa ya x a y Φ dxdydx dy. ( Let us decompose t Q(t) Q(t)w Q(t) Q(t) ) into four terms, and consider for example t Q(t) Q(t)w Q(t) Q(t). Using (4) again and Sobolev s inequality (39), we have ( ) t Q(t) Q(t)w Q(t) Q(t) (x,y;x,y ) Φ,a xa ya x a y Φ dxdydx dy = ( t Q(t))(x,x )w (x y)δ(y y ) Φ,a x a y a x a y Φ dxdydx dy q(t,x) q(t,x ) w (x y) a x a y Φ a x a y Φ dxdydx q(t, ) ( ) /2 L w (x y) a x a y Φ 2 dxdydx ( ) /2 w (x y) q(t,x) 2 a x a y Φ 2 dxdydx Thus C t Φ,2 Φ /2 C ε ε Φ,dΓ( )Φ /2 C t,ε ε m Φ,dΓ( )Φ. Φ, t R 4 Φ C t,ε ε m Φ,dΓ( )Φ. This completes the proof. As a simple consequence of the above estimates, we have an a-priori upper bound for the kinetic energy of Φ. Lemma 0. For all β > 0 we have Φ (t),dγ( )Φ (t) C t,ε ( + 2β ). Proof. Since Ψ (t) is the Schrödinger dynamics, it preserves energy: Ψ (t),h Ψ (t) = Ψ (0),H Ψ (0). From [36, Lemma 2] we have the lower bound H ε ( i ) C ε 2β. i= with fixed ε > 0 small enough. We conclude that Ψ, ( i )Ψ +C ε 2β +C ε Ψ,H Ψ C ε ( + 2β ) i=

20 20 P.T. AM AD M. APIÓRKOWSKI which follows from Ψ (0),H Ψ (0) C. The last inequality follows from the assumption on Ψ (0) and bounds derived in Lemmas 5 and 9. Finally, we have Ψ (t), ( i )Ψ (t) = Φ (t),u (t)dγ( )U (t)φ (t) i= = Φ (t),dγ( )Φ (t) + u(t) 2 dxφ (t),( )Φ (t) [ + Φ (t), ] a(q(t)( )u(t))φ (t) +h.c.. Since u(t) 2 dxφ (t),( )Φ (t) C t and Φ (t), a(q(t)( )u(t))φ (t) Φ a(q(t)( )u(t))φ (t) + Q(t)( )u(t) 2 2 C t we conclude that Φ (t),dγ( )Φ (t) C t,ε ( + 2β ). 6. Truncated dynamics As we explained, the main goal is to compare Φ and Φ. Instead of doing this directly, we introduce an intermediate dynamics living in the sector of very few particles. Related ideas have been used in our works [4, 8] on defocusing 3D systems. Take and let Φ,M (t) be the solution to M = δ, δ (0,) i t Φ,M (t) = M G (t) M Φ,M (t), Φ,M (0) = M Φ (0). (42) Putting differently, Φ,M (t) is the dynamics in the truncated excited Fock space F+ M, generated by the quadratic form G (t) restricted on F+ M. The existence and uniqueness of Φ,M (t) follows from [33, Theorem 7]. Let us briefly explain why Φ,M (t) indeed belongs to F+ M. The identity d dt >M Φ,M (t) 2 = Φ,M (t),i[ M G (t) M, >M ]Φ,M (t) = 0 together with the initial condition >M Φ,M (0) = 0 imply that >M Φ,M (t) = 0 and thus indeed Φ,M (t) belongs to the truncated space with no more than M particles. The fact that Φ,M (t) belongs to the excited space follows from the same argument applied to a(u(t))φ,m (t) 2 (see also proof of [33, Thm. 7]. The kinetic energy of this truncated dynamics is the subject of the following Lemma. When M = δ with a constant δ (0,) we have Φ,M (t),dγ( )Φ,M (t) C t,ε ε, ε > 0.

21 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 2 Proof. By using Lemmas 6, 9 and the assumption M δ we have the following quadratic form estimate on F M : ±(G (t)+dγ( )) 2 dγ( )+C t,ε( + ε ), ε > 0, ± t G (t) C t dγ( ), ±i[g (t),] C t dγ( ). The desired estimate follows by applying Gronwall s argument as in the proof of Lemma 5, with A(t) replaced by M G (t) M +C t,ε ( + ε ). The main result of this section is the following comparison. Lemma 2. When M = δ with a constant δ (0,), Proof. ote that Φ (t) Φ,M (t) 2 C t,ε ε( β M + M /2 ), ε > 0. Φ (t) = Φ (0) = Φ(0), Φ,M (t) = Φ,M (0) = M Φ(0). Therefore, ( ) Φ (t) Φ,M (t) 2 2 RΦ (t),φ,m (t). Take a parameter M/2 m M 3 and write Φ (t),φ,m (t) = Φ (t), m Φ,M (t) +Φ (t), >m Φ,M (t). For the many-particle sectors, using the Cauchy-Schwarz inequality and Lemma we can estimate Φ (t), >m Φ,M (t) Φ (t). >m Φ,M (t) Φ,M (t),(/m)φ,m (t) /2 C t,ε ε M /2. (43) For the few-particle sectors, we use the equations of Φ (t) and Φ,M (t). We have d ( dt Φ (t), m Φ,M (t) = i Φ (t), G (t) m m M G (t) M) Φ,M (t) = Φ (t),i[g (t), m ]Φ,M (t) Here in the last equality we have used m M G (t) M = m G (t) (44) which follows from the choice m M 3 and the fact that G (t) contains at most 2 creation operators and at most 2 annihilation operators. ow we average over m [M/2,M 3]. Lemma 3. When M = δ with a constant δ (0,), then on the excited Fock space we have ± i[g (t), m ] C t,ε ε M/2 2 M dγ( ). Proof. We decompose [G (t), m ] = m G (t) >m >m G (t) m.

22 22 P.T. AM AD M. APIÓRKOWSKI We will focus on >m G (t) m. The other term can be treated similarly. Let us denote A = (Q(t) Q(t)w Q(t) )(x,y;x,y )u(t,x )a x 2 a y a y dxdydx dy, a (Q(t)[w u(t) 2 ]u(t)) =: A 3 +A, A 2 = K 2 (t,x,y)a 2 xa ydxdy. We then have >m G (t) m = >m( ( )( ) ) A +A 2 m ( )( ) = A ( = m)+a 2 (m m). Here we used the fact that A contains one creation operator, while A 2 contains two creation operators. It follows that >m G (t) m = A (M/2 M 3) ( )( ) +A 2 [(M/2 < M 3)+(M/2 < M 3)]. Let us now bound the first term. For any X F +, using Cauchy-Schwarz we have X,A 3 (M/2 M 3)X w (x y) u(t,x) a x a y X 2 a y (M/2 M 3)X dxdy C w (x y) a x a y X 2 dxdy + C w (x y) u(t,x) 2 a y (M/2 M 3)X 2 dxdy C ε ε X,dΓ( )X +C tx,x C t,ε ε dγ( ) where we used (39) and u(t) C t. The terms involving A and A 2 can be bounded in the same way. This ends the proof. Remark 4. From the proof, we also obtain ± M/2 2 i[h(t), m ] C t,ε ε M dγ( ). ow we come back to the proof of Lemma 2. Using Lemma 3, the Cauchy-Schwarz inequality and the kinetic estimates in Lemmas 0,, we can estimate M/2 2 = M/2 2 d dt Φ (t), m Φ,M (t) Φ (t),i[g (t), m ]Φ,M (t) C t,ε ε /2 Φ (t),dγ( )Φ (t) Φ,M (t),dγ( )Φ,M (t) M C t,ε 3 2 ε β, ε > 0. M /2

23 Consequently, FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 23 R M/2 2 R M/2 2 Φ (t), m Φ,M (t) Φ (0), m Φ,M (0) C t,ε 3 2 εβ. M Moreover, since Φ (0) = Φ(0), Φ,M (0) = M Φ(0) and Φ(0),Φ(0) C we have Thus Φ (0), m Φ,M (0) = Φ(0), m Φ(0) R M/2 2 Recall that from (43) we have immediately ±R M/2 2 = Φ(0), >m Φ(0) Φ(0),(/m)Φ(0) C M. Φ (t), m Φ,M (t) C M C t,ε ε β M. Combining the latter two estimates, we arrive at Φ (t), >m Φ,M (t) C t,ε ε M /2. RΦ (t),φ,m (t) C t,ε ε( β M + ) M /2. Therefore, we conclude that Φ (t) Φ,M (t) 2 2( RΦ (t),φ,m (t) ) C t,ε ε( β M + ) M /2, ε > Proof of Theorem (d = 2) In this section, we prove Theorem in case d = 2. We have compare Φ and Φ,M in the previous section. The last main step is the following comparison Lemma 5. When M = δ with a constant δ (0,), Φ,M (t) Φ(t) 2 C t,ε ε( M + M ), ε > 0. Proof. We follow the same strategy as in the proof of Lemma 2. Again, we use and write, with M/2 m M 3, Φ,M (t) Φ(t) 2 2( RΦ,M (t),φ(t) ) Φ,M (t),φ(t) = Φ,M (t), m Φ(t) +Φ,M (t), >m Φ(t). For the many-particle sectors, by using the bounds in Lemmas 5, we have Φ,M (t), >m Φ(t) >m Φ,M (t). >m Φ(t) Φ,M (t),(/m)φ,m (t) /2 Φ(t),(/m)Φ(t) /2 C t,ε ε M. (45)

24 24 P.T. AM AD M. APIÓRKOWSKI For the few-particle sectors, by using the equations of Φ,M (t) and Φ(t) and the identity (cf. (44)) M G (t) M m = G (t) m we can write d dt Φ,M(t), m Φ(t) ( ) = i Φ,M (t), M G (t) M m m H Φ(t) ( ) = i Φ,M (t), (G (t) H) m +[H, m ]) Φ(t). For the term involving (G (t) H), by Lemma 9 we have m ± m+2 (G (t) H) m+2 C t,ε ε dγ( ), ε > 0. Therefore, by the Cauchy-Schwarz inequality and Lemmas 5,, we can estimate Φ,M (t),(g (t) H) Φ(t) m = Φ,M (t), m+2 (G (t) H) m+2 Φ(t) m m /2 /2 C t,ε ε Φ,M (t),dγ( )Φ,M (t) m Φ(t),dΓ( ) m Φ(t) M C t,ε 2ε, ε > 0. For the part involving i[h, m ], we will take the average over m [M/2,]. Recall the remark after Lemma 3: ± M/2 2 i[h(t), m ] C t,ε ε M dγ( ). By the Cauchy-Schwarz inequality and Lemmas 5, again, we can thus bound M/2 2 Φ,M (t),i[h, m ]Φ(t) C t,ε ε /2 Φ,M (t),dγ( )Φ,M (t) Φ(t), dγ( )Φ(t) M C t,ε 2ε, ε > 0. M Thus we have proved that M/2 2 Consequently, R M/2 2 R M/2 2 d dt Φ,M(t), m Φ(t) C t,ε ε( M + M Φ,M (t), m Φ(t) Φ,M (0), m Φ(0) C t,ε ε( M + M /2 ), ε > 0. ).

25 Moreover, Thus R M/2 2 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 25 Φ,M (0), m Φ(0) Φ(0), m Φ(0) C M. Φ,M (t), m Φ(t) C M C t,ε ε( M + M Averaging (45) over m [M/2,M 3] and combining with the latter estimate, we arrive at RΦ,M (t),φ(t) C t,ε ε( M + ), ε > 0, M and hence Φ,M (t) Φ(t) 2 2( RΦ,M (t),φ(t) ) C t,ε ε( M + ), ε > 0. M ow we are ready to conclude Proof of Theorem in case d = 2. Since U (t) is a unitary operator from H to F +, we have the identity Ψ (t) U (t) Φ(t) = U (t)ψ (t) Φ(t) Φ (t) Φ(t). Using the triangle inequality Φ (t) Φ(t) Φ (t) Φ,M (t) + Φ,M (t) Φ(t), where Φ,M is introduced in (42) with M = δ, and using the comparison results in Lemmas 2, 5, we find that Φ (t) Φ(t) 2 C t,ε ε( β M + M ) M /2 +, ε > 0. By choosing M = 2β+ 3 we obtain This ends the proof. Φ (t) Φ(t) 2 C t,ε ε β 3, ε > 0. Proof of Corollary 3. It is well-known that the trace norm controls the partial trace norm, and hence from the norm convergence we can deduce easily that γ lim Tr () Ψ (t) u(t) u(t) = 0 (see [33, Corollary 2] for a detailed explanation). Moreover, since w aδ 0 weakly with a = w, it is straightforward to check that the Hartree solution converges strongly in L 2 (R d ) to the nonlinear Schödinger evolution in (5) with the initial value ϕ(0,x) = u(0,x) and the phase factor µ(t) = lim µ (t) = a u(t,x) 4 dx. 2 R d Thus we conclude that γ lim Tr () Ψ (t) ϕ(t) ϕ(t) = 0 To go from equation (5) to equation (4), we only need to use a gauge transformation ϕ(t) e i t 0 µ(s)ds ϕ(t). ).

26 26 P.T. AM AD M. APIÓRKOWSKI This transformation, however, does not change the projection ϕ(t) ϕ(t) and the desired conclusion follows. 8. Proof of Theorem (d = ) In this section, we prove Theorem in case d =. The scheme of the proof remains the same as in the case of d = 2, but the proof for d = is somewhat simpler, thanks to the stronger Sobolev s inequality in one-dimension. Below we will sketch the main steps again and point out where the differences lie. First, weconsiderthehartreeequation (4). Whenw (x) = β w( β x)withw L (R) andβ > 0, we can show that for every u(0) H 3 (R) theequation (4) has a uniquesolution u(t,.) in H 3 (R) and for all t > 0 we have the bounds u(t, ) u(t, ) H (R) C, u(t, ) H 2 (R) Ce Ct, t u(t, ) Ce Ct where the constant C > 0 is independent of t and (it depends only on w L and u(0) H 3 (R)). This result is obtained by following the same strategy as in the proof of Lemma 4, plus Sobolev s inequality in D (see e.g. [37]) f 2 L (R) f L 2 (R) x f L 2 (R), f H (R). ext, we consider Lemma 8. By following the proof of Lemma 8 and using R d (+ 2πp 2 ) dp < when d =, we obtain the improved version of Lemma 8 when d = : ( x ) /2 K 2 (t,, ) 2 L 2 + ( x ) /2 t K 2 (t,, ) 2 L 2 C t. This leads to an improved estimate in Lemma 6, namely for d = and any η > 0 we have the bounds ( ) ± H(t)+dΓ( ) ηdγ( )+C t (+η ), ± t H(t) ηdγ( )+C t (+η ), ±i[h(t),] ηdγ( )+C t (+η ). Consequently, we have immediately an improved kinetic estimate for the Bogoliubov dynamics Φ(t): Φ(t),dΓ( )Φ(t) Ct. (47) ow we consider the operator inequality (39). When d =, Sobolev s inequality implies that w (x y) C( x ). (48) As a consequence we are able to bound the error term E in Lemma 9 in terms of the kinetic energy on the whole truncated Fock space (not only on the sector with less than m < particles as for d = 2). The one dimensional counterpart of Lemma 9 reads, for any η > 0, ± E (t) ηdγ( )+C t (+η ) ± t E (t) ηdγ( )+C t (+η ) ± i[e (t),] ηdγ( )+C t (+η ). The bounds (46) and (49) imply that we can apply the Grönwall argument to the quantity Φ (t),(g (t) + C t )Φ (t) and derive a better, β-independent bound on the kinetic energy of the full many-body dynamics Φ (t), i.e. (46) (49) Φ (t),dγ( )Φ (t) C t. (50)

27 FOCUSIG DYAMICS FOR LOW DIMESIOAL BOSOS 27 The bounds (47) and (50) imply in particular that now, when d =, the bounds in Lemmas 2 and 5 become independent of β and ε, namely Φ (t) Φ,M (t) 2 C t M /2 and ( Φ,M (t) Φ(t) 2 M C t + ). M Then we can proceed exactly as in the proof for d = 2 in Section 7 and arrive at the desired final estimate. References [] R. Adami, F. Golse, and A. Teta, Rigorous Derivation of the Cubic LS in Dimension One, J. Stat. Phys., 27 (6) (2007), pp [2] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, E.A. Cornell, Observation of Bose Einstein Condensation in a Dilute Atomic Vapor, Science 269 (995), pp [3] V. Bach, S. Breteaux, T. Chen, J. Fröhlich and I.M. Sigal, The time-dependent Hartree-Fock- Bogoliubov equations for Bosons, e-print arxiv: (205). [4] C. Boccato, S. Cenatiempo, and B. Schlein, Quantum many-body fluctuations around nonlinear Schrödinger dynamics, Annales Henri Poincaré, 8 (207), pp [5] C. Boccato, C. Brennecke, S. Cenatiempo and B. Schlein, The excitation spectrum of Bose gases interacting through singular potentials, Preprint arxiv: (207). [6].. Bogoliubov, About the theory of superfluidity, Izv. Akad. auk SSSR, (947), p. 77. [7] S.. Bose, Plancks Gesetz und Lichtquantenhypothese, Zeitschrift für Physik 26 (924), pp [8] C. Brennecke, P.T. am. M. apiórkowski and B. Schlein, Fluctuations of -particle quantum dynamics around the nonlinear Schrödinger equation, Preprint (207). [9] H. Brezis and T. Gallouet, onlinear Schrödinger evolution equations, onlinear Analysis: Theory, Methods and Applications 4 (980), pp [0] H. Brezis and S. Wainger, A note on limiting cases of Sobolev embeddings, Comm. Partial Diff. Equations 5 (980), pp [] T. Cazenave, Semilinear Schrödinger Equations, Courant Lecture otes vol. 0, American Mathematical Society, [2] X. Chen, Collapsing Estimates and the Rigorous Derivation of the 2d Cubic onlinear Schrödinger Equation with Anisotropic Switchable Quadratic Traps, J. Math. Pure. Appl. 98 (202), no. 5, [3] X. Chen and J. Holmer, The Rigorous Derivation of the 2D Cubic Focusing LS from Quantum Many Body Evolution, International Mathematics Research otices, Vol. 206, o. 00, pp. 44 doi: 0.093/imrn/rnw3 [4] X. Chen and J. Holmer, Focusing Quantum Many-body Dynamics: The Rigorous Derivation of the D Focusing Cubic onlinear Schrödinger Equation, Arch. Rational Mech. Anal. 22 (206), pp doi:0.007/s [5] J.J.W Chong, Dynamics of Large Boson Systems with Attractive Interaction and A Derivation of the Cubic Focusing LS in R 3, Preprint: arxiv: [6] J. Dereziński and M. apiórkowski, Excitation spectrum of interacting bosons in the mean-field infinite-volume limit, Annales Henri Poincaré, 5 (204), pp Erratum: Annales Henri Poincaré 6 (205), pp [7] L. Erdös, B. Schlein, and H.-T. Yau, Derivation of the cubic non-linear Schrödinger equation from quantum dynamics of many-body systems, Invent. Math., 67 (2007), pp [8], Rigorous derivation of the Gross-Pitaevskii equation with a large interaction potential, J. Amer. Math. Soc., 22 (2009), pp [9], Derivation of the Gross-Pitaevskii equation for the dynamics of Bose-Einstein condensate, Ann. of Math. (2), 72 (200), pp [20] A. Einstein, Quantentheorie des einatomigen idealen Gases, Sitzungsberichte der Preussischen Akademie der Wissenschaften (925), pp. 3. [2] J. Ginibre and G. Velo, The classical field limit of scattering theory for nonrelativistic many-boson systems. I, Commun. Math. Phys., 66 (979), pp [22] P. Grech and R. Seiringer, The excitation spectrum for weakly interacting bosons in a trap, Comm. Math. Phys., 322 (203), pp [23] M. Jeblick and P. Pickl, Derivation of the Time Dependent Two Dimensional Focusing LS Equation, Preprint arxiv: (207).

28 28 P.T. AM AD M. APIÓRKOWSKI [24] M. Jeblick,. Leopold and P. Pickl, Derivation of the Time Dependent Gross-Pitaevskii Equation in Two Dimensions, Preprint arxiv: (206). [25] K.B. Davis, M.-O. Mewes, M.R. Andrews,.J. van Druten, D.S. Durfee, D.M. Kurn, W. Ketterle, Bose Einstein condensation in a gas of sodium atoms, Phys. Rev. Lett. 75 (995), pp [26] K. Kirkpatrick, B. Schlein and G. Staffilani, Derivation of the two-dimensional nonlinear Schrödinger equation from many body quantum dynamics, American Journal of Mathematics, 33(), 9-30 (20) [27] M. Grillakis and M. Machedon, Pair excitations and the mean field approximation of interacting Bosons I, Commun. Math. Phys., 324 (203), pp [28] M. Grillakis and M. Machedon, Pair excitations and the mean field approximation of interacting Bosons II, Communications in Partial Differential Equations 42 (207), [29], Second-order corrections to mean field evolution of weakly interacting bosons. II, Adv. Math., 228 (20), pp [30] Y. Guo and R. Seiringer. On the Mass Concentration for BoseEinstein Condensates with Attractive Interactions, Lett. Math. Phys. 04 (204), [3] K. Hepp, The classical limit for quantum mechanical correlation functions, Comm. Math. Phys., 35 (974), pp [32] E. Kuz, Exact Evolution versus Mean Field with Second-order correction for Bosons Interacting via Short-range Two-body Potential, Differential and Integral Equations 30 (207), [33] M. Lewin, P. T. am, and B. Schlein, Fluctuations around Hartree states in the mean-field regime, Amer. J. Math., 37 (205), pp [34] M. Lewin, P. T. am, S. Serfaty, and J. P. Solovej, Bogoliubov spectrum of interacting Bose gases, Comm. Pure Appl. Math., 68 (205), pp [35] M. Lewin, P.T. am, and. Rougerie, The mean-field approximation and the nonlinear Schrödinger functional for trapped Bose gases, Trans. Amer. Math. Soc. 369 (206), [36] M. Lewin, P.T. am, and. Rougerie, A note on 2D focusing many-boson systems, Proc. Amer. Math. Soc. 45 (207), [37] E. H. Lieb and M. Loss, Analysis, Graduate Studies in Mathematics, American Mathematical Society, 200. [38] M. Maeda, On the symmetry of the ground states of nonlinear Schrdinger equation with potential. Adv. onlinear Stud. 0 (200), [39] D. Mitrouskas, S. Petrat and P. Pickl, Bogoliubov corrections and trace norm convergence for the Hartree dynamics, e-print arxiv: (206) [40] P. T. am and M. apiórkowski, Bogoliubov correction to the mean-field dynamics of interacting bosons, Adv. Theor. Math. Phys. 2 (207), [4] P. T. am and M. apiórkowski, A note on the validity of Bogoliubov correction to mean-field dynamics, J. Math. Pure. Appl. 08 (207), [42] P. T. am, M. apiórkowski and J. P. Solovej, Diagonalization of bosonic quadratic Hamiltonians by Bogoliubov transformations, J. Funct. Anal., 270 () (206), pp [43] P. T. am and R. Seiringer, Collective excitations of Bose gases in the mean-field regime, Arch. Rational Mech. Anal., 25 (205), pp [44] P. Pickl, Derivation of the Time Dependent Gross Pitaevskii Equation Without Positivity Condition on the Interaction, J. Stat. Phys., 40 (200), pp [45] R. Seiringer, The excitation spectrum for weakly interacting bosons, Commun. Math. Phys., 306 (20), pp [46] H. Spohn, Kinetic equations from Hamiltonian dynamics: Markovian limits, Rev. Modern Phys., 52 (980), pp [47] M. I. Weinstein, onlinear Schrödinger equations and sharp interpolation estimates, Commun. Math. Phys. 87 (983), pp [48] J. Zhang, Stability of Attractive Bose Einstein Condensates, J. Stat. Phys. 0 (2000), Department of Mathematics, LMU Munich, Theresienstrasse 39, Munich, Germany address: nam@math.lmu.de Department of Mathematical Methods in Physics, Faculty of Physics, University of Warsaw, Pasteura 5, Warszawa, Poland address: marcin.napiorkowski@fuw.edu.pl