arxiv: v3 [math-ph] 23 Apr 2013

Transcription

1 O THE RIGOROUS DERIVATIO OF THE 3D CUBIC OLIEAR SCHRÖDIGER EQUATIO WITH A QUADRATIC TRAP XUWE CHE Dedicated to Xuqing. arxiv:14.15v3 [math-ph] 3 Apr 13 Abstract. We consider the dynamics of the 3D -body Schrödinger equation in the presence of a quadratic trap. We assume the pair interaction potential is 3β 1 V β x. We justify the mean-field approximation and offer a rigorous derivation of the 3D cubic LS with a quadratic trap. We establish the space-time bound conjectured by Klainerman and Machedon [3] for β,/7] by adapting and simplifying an argument in Chen and Pavlović [7] which solves the problem for β,1/4 in the absence of a trap. Contents 1. Introduction Comparison with Chen and Pavlović [7] The Anisotropic Version of the Main Theorem 5. Proof of the Main Theorem Theorem Lens Transform 7 4. Proof of Theorem The Uniqueness of Hierarchy 4. Proof of Theorem The Space-Time Bound of the BBGKY Hierarchy for β,/7] Proof of Theorem Conclusion 1 8. Acknowledgements 9. Appendix I: Proof of Corollary 1 1. Appendix II: Proof of Theorem 7 3 References 8 1. Introduction It is widely believed that the cubic nonlinear Schrödinger equation LS i t φ = Lφ+ φ φ, where L is the Laplacian or the Hermite operator +ω x, describes the physical phenomenon of Bose-Einstein condensation BEC. This belief is one of the main motivations for studying the cubic LS. BEC is the phenomenon that particles of integer spin Bosons occupy a macroscopic quantum state. This unusual state of matter was first predicted theoretically by Einstein for non-interacting particles. The first experimental observation of BEC in an interacting atomic gas did not occur until 1995 using laser cooling techniques [1, 14]. E. A. Cornell, W. Ketterle, and C. E. Wieman were awarded the 1 obel Prize in Physics for observing BEC. Many similar successful experiments were performed later on [4, 13, 7, 37]. Let t R be the time variable and x = x 1,x,...,x R 3 be the position vector of particles in R 3. Then BEC naively means that the -body wave function ψ t,x satisfies ψ t,x = φt,x j 1.1 up to a phase factor solely depending on t, for some one particle state φ. In other words, every particle is in the same quantum state. It is believed that the one particle wave function φ which models the condensate satisfies the cubic LS. Gross [5, 6] and Pitaevskii [34] proposed such a description. However, the cubic LS is a phenomenological mean field type equation and its validity needs to be established rigorously from Date: V3 for ARMA, 4/4/13. 1 Mathematics Subject Classification. Primary 35Q55, 35A, 81V7; Secondary 35A3, 35B45, 81Q5. Key words and phrases. Gross-Pitaevskii Hierarchy, Quadratic Trap, Collapsing Estimate, Klainerman-Machedon Board Game and Space-time Bound. 1

2 XUWE CHE the many body system which it is supposed to characterize. As a result, we investigate the procedure of laboratory experiments of BEC according to [1, 14]. Step A. Confine a large number of Bosons inside a trap e.g., the magnetic fields in [1, 14]. Cool it down so that the many body system reaches its ground state. It is expected that this ground state is a BEC state / factorized state. This step corresponds to the mathematical problem. Problem 1. Show that the ground state of the -body Hamiltonian 1 x j + ω x j + 1 3β V β x i x j is a factorized state. 1 i<j We use the quadratic potential x to represent the trap. This simplified yet reasonably general model is expected to capture the salient features of the actual trap: on the one hand the quadratic potential varies slowly, on the other hand it tends to as x. In the physics literature, Lieb, Seiringer and Yngvason remarked in [31] that the confining potential is typically x in the available experiments. Mathematically speaking, the strongest trap we can deal with in the usual regularity setting of LS is the quadratic trap since the work [39] by Yajima and Zhang points out that the ordinary Strichartz estimates start to fail as the trap exceeds quadratic. Step B. Switch the trap in order to enable measurement or direct observation. It is assumed that such a shift of the confining potential is instant and does not destroy the BEC obtained from Step A. To be more precise about the word switch : in [1, 14] the trap is removed, in [37] the initial magnetic trap is switched to an optical trap, in [4] the trap is enhanced, in [13] the trap is turned off in spatial directions to generate a D Bose gas. Hence we have a different trap after the switch, in other words, the trapping potential becomes ω x. The system is then time dependent unless ω = ω. Therefore, the factorized structure obtained in Step A must be preserved in time for the observation of BEC. Mathematically, this step stands for the following problem. Problem. Take the BEC state obtained in Step A. as initial datum, show that the solution to the many body Schrödinger equation i t ψ = 1 x j +ω x j ψ + 1 3β V β x i x j ψ 1. is a BEC state / factorized state. 1 i<j We first remark that neither of the problems listed above admits a factorized state solution. It is also unrealistic to solve the equations in Problems 1 and for large. Moreover, both problems are linear so that it is not clear how the cubic LS arises from either problem. Therefore, in order to justify the statement that the cubic LS depicts BEC, we have to show mathematically that, in an appropriate sense, ψ t,x φt,x j as for some one particle state φ which solves a cubic LS. However, when φ φ φt,x j φ t,x j as. L i.e. our desired limit the BEC state is not stable against small perturbations. One way to circumvent this difficulty is to use the concept of the k-particle marginal density γ k associated with ψ defined as γ k t,x k;x k = ψ t,x k,x k ψ t,x k,x k dx k, x k,x k R 3k 1.3 and show that γ k t,x k;x k k φt,x j φt,x j as. Penrose and Onsager[33] suggested such a formulation. Another approach is to add a second order correction to the mean field approximation. See [11, 3, 4]. For Problem 1, Lieb, Seiringer, Solovej and Yngvason showed that the ground state of the Hamiltonian exhibits complete BEC in [3], provided that the trapping potential V trap x satisfies inf x >R V trap x

3 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 3 for R and the interaction potential is spherically symmetric. To be more precise, let ψ, be the ground state, then γ 1, φ GP φ GP as, where γ 1, is the corresponding one particle marginal density defined via formula 1.3 and φ GP minimizer of the Gross-Pitaevskii energy functional with coupling constant 4πa 1, φ +V trap x φ +4πa φ 4 dx. is the So far, there has not been any work regarding the Hamiltonian evolution in Step B in the case when ω. Motivated by the above considerations, we aim to investigate the evolution of a many-body Boson system with a quadratic trap and study the dynamics after the switch of the trap. We derive rigorously the 3D cubic LS with a quadratic trap from the body linear equation 1.. To be specific, we establish the following theorem in this paper. Theorem 1. Main Theorem Let { γ k } be the family of marginal densities associated with ψ, the solution of the -body Schrödinger equation 1. for some β, 7]. Assume that the pair interaction V is a nonnegative L 1 R 3 H R 3 W, R 3 spherically symmetric function. Moreover, suppose the initial datum of equation 1. verifies the following the conditions: a the initial datum is normalized i.e. b the initial datum is asymptotically factorized ψ L = 1, for some one particle wave function φ H 1 R 3. c the initial datum has bounded energy per particle i.e. γ lim Tr 1,x 1;x 1 φ x 1φ x 1 =, 1.4 sup 1 ψ,h ψ <, where the Hamiltonian H is H = 1 x j +ω x j i<j 3β V β x i x j. Then t, k 1, we have the convergence in the trace norm that lim Tr γk t,x k;x k φt,x j φt,x j =, where φt,x is the solution to the 3D cubic LS with a quadratic trap i t φ = φ,x = φ x 1 x +ω x and the coupling constant b = R Vxdx. 3 { For the ω = case, the approach which uses the marginal densities φ+b φ φ in R γ k } for the dynamics problem has been proven to be successful in the fundamental papers [15, 16, 17, 18, 19,, 1] by Elgart, Erdös, Schlein, and Yau. As pointed out in [1], their work corresponds to the evolution after the removal of the traps. Motivated by a kinetic formulation of Spohn [36], their { program } consists of two principal parts: in one part, they prove that an appropriate limit of the sequence γ k as solves the Gross-Pitaevskii hierarchy i t + 1 x k 1 k x γ k = b k B j,k+1 γ k+1, k = 1,...,n, Here a is the limit of times the scattering length of 3β 1 V β. i.e. a is the scattering length of V when β = 1 and a should be V /8π when β,1. See [18].

4 4 XUWE CHE where B j,k+1 = Bj,k+1 1 B j,k+1 and Bj,k+1 1 γ k+1 t,x k ;x k = γ k+1 t,x k ;x k = B j,k+1 δx j x k+1 δx k+1 x k+1γ k+1 t,x k+1 ;x k+1dx k+1 dx k+1, δx j x k+1δx k+1 x k+1 γk+1 t,x k+1 ;x k+1 dx k+1dx k+1 ; in another part, they show that hierarchy 1.6 has a unique solution which is therefore a completely factorized state. However, as remarked by Terence Tao, the uniqueness theory for hierarchy 1.6 is surprisingly delicate due to the fact that it is a system of infinitely many coupled equations over an unbounded number of { variables. In [3], by assuming a space-time bound on the limit of γ k }, Klainerman and Machedon gave another proof of the uniqueness in [18] through a collapsing estimate originated from the ordinary multilinear Strichartz estimates in their null form paper [9] and a board game argument inspired by the Feynman graph argument in [18]. Later, the method in Klainerman and Machedon [3] was taken up by Kirkpatrick, Schlein, and Staffilani [8], who studied the corresponding problem in D; by Chen and Pavlović [5, 6], who considered the 1D and D 3-body interaction problem and the general existence theory of hierarchy 1.6; and by the author [1], who investigated the trapping problem in D. In [8, 9], Chen, Pavlović and Tzirakis worked out the virial and Morawetz identities for hierarchy 1.6. Recently, for the 3D case without traps, Chen and Pavlović [7] proved that, for β,1/4, the limit of { γ k } actually satisfies the space-time bound assumed by Klainerman and Machedon [3] as. This has been a well-known open problem in the field. Moreover, they showed that the solution to the BBGKY hierarchy converges strongly to the solution to hierarchy 1.6 in H 1 without assuming asymptotically factorized initial datum. In this paper, we adapt and simplify their argument in establishing the Klainerman-Machedon space-time bound. We also extend the range of β from, 1/4 in Chen and Pavlović[7] to, /7]. Through simple functional analysis, we obtain a convergence result without assuming asymptotically factorized initial datum as well. See Corollary 1 But we are not claiming it as a main result in this paper. We compare our result and the one in Chen and Pavlović [7] briefly in Section 1.1 below Comparison with Chen and Pavlović [7]. For comparison purpose, we transform Theorem 4 in Section, which implies Theorem 1, into the general convergence result below without the assumption of asymptotically factorized initial data since the result in Chen and Pavlović [7] is under the same regularity setting condition 1.7 as Theorem 4. Corollary 1. For comparison purpose only, not a main result. Let γ k t,x k;x k be the k-marginal density associated with the solution of the -body Schrödinger equation 1. for some β,/7] and γ k t,x k ;x k be the solution to the Gross-Pitaevskii hierarchy with a quadratic trap hierarchy 3.4. Assume the initial data of γ k and γk satisfy the following conditions a γ lim Tr k,x k;x k γ k,x k ;x k =, b ψ,hψ k C k k. 1.7 Then we have the convergence of the evolution in trace norm γ lim Tr k t,x k;x k γk t,x k ;x k =. Proof. This result for the non-trap case should be credited to Erdös, Schlein and Yau though they did not state it as one of their main theorems. They mentioned it on page 97 of their paper [1]. Once we have established Theorem 4, its proof together with some simple functional analysis proves this corollary. We include the proof in Appendix I Section 9 for completeness. Briefly, the main result Theorem 3.1 in Chen and Pavlović [7] is the following. Theorem. [7] Let ω =. Suppose γ k t,x k;x k is the k-marginal density associated with the solution of the -body Schrödinger equation 1. for some β,1/4 and γ k t,x k ;x k is the solution to the ordinary Gross-Pitaevskii hierarchy without traps. If the initial data of γ k and γk satisfy the following conditions a lim 1 xj x γ k j,x k;x k γk,x k ;x k =, L dx k dx k

5 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 5 b ψ,h k ψ C k k, then we have the convergence of the evolution in H 1 norm 1 1 lim 1 xj 1 x γ k j t,x k;x k γk t,x k ;x k =. L dx k dx k One can then easily tell the following: on the one hand, the result in this work allows a quadratic trap and a larger range of β in the analysis these are also the main novelty and the main technical improvement of this paper; on the other hand, the result in [7] yields a stronger convergence H 1 convergence when the initial data admits a stronger convergence. The main purpose of this paper is to justify the mean-field approximation and offer a rigorous derivation of the 3D cubic LS with a quadratic trap Theorem 1. Thus we establish the convergence of probability densities which is the trace norm convergence in this context. The Chen-Pavlović result is crucial for their work on the Cauchy problem of Gross-Pitaevskii hierarchies. Before we delve into the proofs, we remark that taking the coupling level k to be ln in Section 6..3 is exactly the place where we follow the original idea of Chen and Pavlović [7]. 1.. The Anisotropic Version of the Main Theorem. It is of interest to use anisotropic traps in laboratory experiments. See, for example, [13]. Our proof for Theorem 1 also applies to the case with anisotropic traps. In fact, we have the following theorem. { Theorem 3. Define a 3*3 diagonal matrix Q whose diagonal entries are ω 1,ω,ω 3. Let family of marginal densities associated with ψ, the solution of the -body Schrödinger equation i t ψ = 1 x j + 1 xt j Qx j ψ + 1 3β V β x i x j ψ 1 i<j γ k } be the for some β, 7]. Assume that the pair interaction V is a nonnegative L 1 R 3 H R 3 W, R 3 spherically symmetric function. Moreover, suppose the initial datum of equation 1. verifies Conditions a, b, and c in the assumption of Theorem 1. Then t, k 1, we have the convergence in the trace norm that lim Tr t,x k;x k φt,x j φt,x j =, γk where φt,x is the solution to the 3D cubic LS with anisotropic quadratic traps i t φ = 1 x + 1 xt Qx φ+b φ φ in R 3+1 and b = R 3 Vxdx. φ,x = φ x The anisotropic version of the main theorem stated above yields to the same techniques as Theorem 1 but is technically more complicated. Therefore we prove only the latter in detail, merely suggesting during the course of the proof the appropriate modifications needed to obtain the more general theorem.. Proof of the Main Theorem Theorem 1 We establish Theorem 1 with a smooth approximation argument and the following theorem. { } Theorem 4. Let be the family of marginal densities associated with ψ, the solution of the - γ k body Schrödinger equation 1. for some β, 7]. Assume that the pair interaction V is a nonnegative L 1 R 3 H R 3 W, R 3 spherically symmetric function. Moreover, suppose the initial datum of equation 1. is normalized and asymptotically factorized in the sense of a and b of Theorem 1, and verifies the condition that there is a C independent of k or such that ψ,h k ψ C k k Then t and k 1, we have the convergence in the trace norm that lim Tr k t,x k;x k φt,x j φt,x j =, γk

6 6 XUWE CHE where φt,x is the solution to the 3D cubic LS with a quadratic trap and b = R 3 Vxdx. i t φ = φ,x = φ x 1 x +ω x φ+b φ φ in R 3+1 Before presenting the proof of Theorem 4, we discuss how to deduce Theorem 1 from Theorem 4. It is a well-known smooth approximation argument. We include it for completeness. For technical details, we refer the readers to Erdös-Schlein-Yau [, 1] and Kirkpatrick-Schlein-Staffilani [8]. Write the spaces of compact operators and trace class operators of L R 3k as K k and L 1 k. Then K k = { L 1 k. Via } the fact that K k is separable, we select a dense countable subset of the unit ball of K k and call it J k i K k. We have J k i 1 where is the operator norm. We set up the following metric on i 1 L 1 k, for γk, γ k L 1 k : define d k γ k, γ k = i TrJ k i γ k γ k. i=1 Then a uniformly bounded sequence γ k L1 k converges to γk L 1 k with respect to the weak* topology if and only if lim d kγ k,γk =. Fix κ > and χ Cc R, with χ 1, χs = 1, for s 1, and χs =, for s. We regularize the initial data of the -body Schrödinger equation 1. with ψ = χκh /ψ χκh /ψ, and we denote { ψ } t the solution of the -body Schrödinger equation 1. subject to this regularized initial data and γ k t the family of marginal densities associated with ψ t. k=1 With these notations, if κ > small enough, on the one hand, we have ψ,h ψ k C k k, and lim TrJk γ k k φ φ =, for every J k K k. Thus γ k,x k;x k k φ φ as in the weak*topology. Since k φ φ is an orthogonal projection, the convergence in the weak* topology is equivalent to the convergence in the trace norm. Consequently, the conditions of Theorem 4 are verified and it implies that t [,T ] and k 1, lim Tr γk t,x k;x k φt,x j φt,x j =. On the other hand, there is a constant C independent of and κ such that k J TrJ k γ t γk t C k κ 1 for every J k K k. Therefore, k TrJk γ k t,x k;x k φt,x j φt,x j TrJ k γ k + t γk t k TrJk γ k t φt,x j φt,x j CJ k κ 1 + TrJ k k γ k t φt,x j φt,x j.

7 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 7 As the above inequality holds for all κ > small enough, we know γ k t,x k;x k k φt,x jφt,x j as in the weak* topology. This convergence is again equivalent to the convergence in the trace norm because k φt,x jφt,x j is an orthogonal projection as well. Whence we have established our main theorem Theorem 1 through Theorem 4. It remains to show Theorem 4. We prove Theorem 4 with the help of the lens transform. We define the lens transform and show its related properties in the next section then we establish Theorem 4 in Section Lens Transform In this section, we first define the lens transform and review its relevant properties, then we prove an energy estimate Proposition which relates the energy on the two sides of the lens transform. It aids in the proof of Theorem 4 in the sense that it links the analysis of x +ω x to the analysis of y which is a better understood operator. We denote t,x the space-time on the Hermite side and,y the space-time on the Laplacian side. We now define the lens transform we need. Definition 1. Let x,y R 3. We define the lens transform for L functions M : L dy L dx and its inverse by x M u t,x = e iωtanωt u cosωt 3 tanωt ω, x cosωt ω M 1 ψ y,y = ei 1+ω 1+ω 3 4 M is unitary by definition and the variables are related by ψ arctanω, ω y 1+ω. = tanωt ω, y = x k cosωt. Definition. Let x k,x k,y k,y k R3k. We define the lens transform for Hilbert-Schmidt kernels T k : L dy k dy k L dx k dx k and its inverse by e iωtanωt xk T k u k t,x k ;x k = cosωt 3k x k ω y k y k T 1 k γk,y k ;y k = ei 1+ω 1+ω 3k u k tanωt ω, x k cosωt ; x k cosωt γ k arctanω, ω T k is unitary by definition as well and the variables are again related by y k 1+ω ; y k ω. = tanωt ω, y k = x k cosωt and y k = x k cosωt. Before we characterize the exact effect of the lens transform, we clarify the motivation of such definitions by a lemma. Lemma 1. [3, 1] Define α and β via the system αt+ηtαt =,α =, α = 1, βt+ηtβt =,β = 1, β =. If β is nonzero in the time interval [,T], then the solution of the 1D Schrödinger equation with a time dependent quadratic trap i t ψ = 1 x + 1 ηtx ψ in R 1+1 ψ,x = fx L R in [,T] is given by βt x ψt,x = ei βt u αt βt 1 βt, x βt, if u,y solves the 1D free Schrödinger equation subject to the same initial data. i u = 1 yu in R 1+1

8 8 XUWE CHE Proof. See [3] for a proof by direct computation and [1] for an algebraic proof using the metaplectic representation. When ηt = ω, such a transformation has a long history, we refer the readers to [3] and the references within. To make formulas shorter, let us write formula 3.1 as T 1 k γk,y k ;y k = γ k t,x k ;x kh k,y k ;y k where h k,y k ;y k represents the factor away from γk in formula 3.1, then precisely, the lens transform has the following effect. Proposition 1. i + 1 y k 1 y k T 1 k γk,y k ;y k = hk,y k ;y k 1+ω i t 1 x k +ω x k + Proof. This is a direct computation. 1 x k +ω x k γ k t,x k ;x k Via Proposition 1, we know how the lens transform acts on the Bogoliubov Born Green Kirkwood Yvon BBGKY hierarchy and the Gross-Pitaevskii hierarchy. Lemma. BBGKY hierarchy under the lens transform Write V x = 3β V β x { }. γ k solves the 3D BBGKY hierarchy with a quadratic trap i t γ k 1 x k +ω x k γ k + 1 x +ω x k γ k k 3. = 1 V x i x j V x i x j γ k 1 i<j k + k k V x j x k+1 V x j x k+1 γ k+1 t,x k,x k+1 ;x k,x k+1 dx k+1, { } in [ T,T ] if and only if u k = T 1 k γk solves the hierarchy i + 1 y k 1 y k u k y i y j y i = 1+ω V V 1 i<j k 1+ω 1 y j u k 1+ω 1 + k 1 k 1+ω V y j y k+1 V 1+ω 1 y j y k+1 1+ω 1 in [ tanωt ω ], tanωt ω. u k+1,y k,y k+1 ;y k,y k+1dy k+1, Lemma 3. Gross-Pitaevskii hierarchy under the lens transform { γ k} solves the 3D Gross-Pitaevskii hierarchy with a quadratic trap i t γ k 1 x k +ω x k γ k + 1 k x +ω x k γ k = b k B j,k+1 γ k+1, 3.4 in [ T,T ] if and only if { u k = T 1 k γk} solves the hierarchy i + 1 y k 1 y k u k = k 1+ω 1 b B j,k+1 u k+1, 3.5 in [ tanωt ω ], tanωt ω. The lens transform for the Hilbert-Schmidt kernels is not only by definition a unitary transform on L R 6k, it is also an isometry on the space of the self-adjoint trace class operator kernels.

9 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 9 Lemma 4. t [ T,T ], Ky k,y k the kernel of a self-adjoint trace class operator on L R 3k. If Ky k,y k fy k dy k = λfy k, then T k Kx k,x k M fx k dx k = λm fx k. In other words, the eigenvectors of the kernel T k Kx k,x k are exactly the lens transform Lemma 1 of the eigenvectors of the kernel Ky k,y k with the same eigenvalues. In particular, we have Tr T k K = Tr K. Proof. This is a straight forward computation. We remark that we have defined the generalized lens transform for a function and a kernel separately via Definitions 1 and. Once we have proved the following proposition which relates the energy of the two sides of the lens transform, we can start the proof of Theorem 4. Proposition. Let ψ t,x be the solution to equation 1. for some β,3/5 subject to initial ψ which satisfies the energy condition ψ,h k ψ C k k. If u,y = M 1 ψ, then there is a C, for all k, k such that u, 1 yj u C k, for all and all [ tanωt ω, tanωt ω ] provided that T < π ω. The rest of this section is the proof of Proposition. We prove it for ω > through Lemmas 5 and 6 since the case ω = has already been studied in [15, 18]. Lemma 5. For β,3/5, there is a C, for all k, k such that ϕ,h k ϕ C k k ϕ, xj +ω x j ϕ, for all and all ϕ L s R 3. Proof. TheproofbasicallyfollowsProposition1in[15]stepbystepifonereplaces 1 xj by xj +ω x j and notices that for ω >, 1 α f L C α +ω x α when one uses Sobolev. There are some extra error terms which can be easily handled. We illustrate the control of the extra error terms through the following example. Write S i = xi +ω x i, consider 3β ϕ,s 1...S n+1 V β x 1 x ϕ +c.c. f L [38] = 3β S 3...S n+1 ϕ,s 1 S V β x 1 x S 3...S n+1 ϕ +c.c., where c.c. denotes complex conjugates. eglecting S 3...S n+1, we have 3β ϕ,s1 S V β x 1 x ϕ +c.c. = 3β ϕ, x1 +ω x 1 x +ω x V β x 1 x ϕ = 3β ϕ, x1 x V β x 1 x ϕ +c.c. + 3β ϕ, x1 ω x V β x 1 x ϕ +c.c. + 3β ϕ, ω x 1 x V β x 1 x ϕ +c.c. + 3β ϕ, ω x 1 ω x V β x 1 x ϕ +c.c. = I +II +III +IV +c.c.

10 1 XUWE CHE Compared to [15], the extra error terms are II,III, and IV. Since IV is positive, we only look at II 4β x1 ω x ϕ, V β x 1 x ω x ϕ C 4β α x1 ω x ϕ, V β x 1 x x1 ω x ϕ C 4β α 1 ω x ϕ, V β x 1 x ω x ϕ C 4β α ϕ,s 1 S ϕ +α 1 β ω x ϕ,s 1 ω x ϕ Sobolev at the second term C 3β ϕ,s 1 S ϕ α = β C 3β ϕ,s 1 S ϕ As long as β, 3, we can absorb the extra error terms into the main term ϕ,s 1 S ϕ. The Sobolev we used is V β x 1 x ϕ dx 1 dx β V L 3 ϕ 6 dx dx C β I V 3 x1 L 1 ϕ dx 1 dx C β V 3 L x1 +ω x 1 1 ϕ Lemma 6. Let P x t = i x cosωt ωxsinωt. If u,y = M 1 ψ, then u, yj u = P xj tψ t,p xj tψ t = ψ t,px j tψ t. dx 1 dx. Proof. We provide a proof through direct computation here. We remark that P x t is in fact the evolution of momentum. See [1]. Without lose of generality, we may assume = 1, then tan ωt iω e x P x tψ 1 t = P x t cosωt 3 u 1 tanωt ω, x cos ωt tan ωt = e iω x tanωt i cosωt 3 x cosωtu 1 ω, x cos ωt = e iω tan ωt x cosωt 3 i y u 1,y. Thus u 1, y u 1 = y u 1,y y u 1,y dy = cosωt 3 dx = P x tψ 1 t,p x tψ 1 t. Proof of Proposition. We first notice that u, 1 yj u = ψ t, Since f, 1+Px j t f 1+Px j t ψ t Lemma 6. = f,f +cos ωt f, f +ω sin ωt xf,xf ωsinωtcosωtim f,xf C f, x +ω x f, we know ψ t, 1+Px j t ψ t C k ψ t, xj +ω x j ψ t Ck k ψ t,h k ψ t. Lemma 5

11 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 11 From the energy condition on the initial datum ψ, we then deduce Proposition that is u, 1 yj u C k. 4. Proof of Theorem 4 We devote this section to establishing Theorem 4. The main idea is to first prove that, in the time period [ ] { }, tanωt ω with T < π ω, as, u k = T 1 k γk, the lens transform of the solution to the BBGKY hierarchy 3., converges to { u k = T 1 k γk}, the lens transform of the solution to the Gross- Pitaevskii hierarchy 3.4 in the trace norm, then use Lemma 4 to conclude the convergence γ k γk in the trace norm as in the time period t [,T ]. Then the time translation invariance of equation 1. proves Theorem 4 for all time. Step I. In Proposition, we have already established the energy estimate for { u = M 1 ψ }, sup u, 1 yj u C k. [, tan ωt ω ] which becomes { } u k = T 1 k γk sup Tr [, tan ωt ω ] 1 yj u k Ck, 4.1 for. Therefore we can utilize the proof in Erdös-Schlein-Yau [18] or Kirkpatrick- { } Schlein-Staffilani [8] to show that the sequence u k = T 1 k γk is compact with respect to the weak* topology on the trace class operators and every limit point { u k} solves hierarchy 3.5 which, we recall, is i + 1 y k 1 y k u k = k 1+ω b B j,k+1 u k This is a fixed time argument. We omit the details here. otation 1. To make formulas shorter, let us write g = 1+ω 1 from here on. The only property of g we are going to need is that < c g C < in any finite time period. Step II. In this step, we use the a-priori estimate 4.1 to provide a space-time bound of u k so that we can employ Theorem 6 in Step III. We transform estimate 4.1 into the following theorem. Theorem 5. Main Auxiliary Theorem Assume V is a nonnegative L 1 R 3 H R 3 function. Let Ṽ, y = 3β Ṽ β y = 3β g 3 Vg β y be the interaction potential with the interaction parameter β ], k 7. Suppose that u solves hierarchy 3.3 in [,T] [ ], tanωt ω, which, written in the integral form, is u k = U k u k, i 1 i<j k i k subject to the condition that k U k s sup [,T] U k s Tr Ṽ,s y i y j Ṽ,sy i y j B,j,k+1,s gs gs u k+1 ds, u k s,y k;y k ds yj u k Ck, 4.4

12 1 XUWE CHE where B,j,k+1, = B,j,k+1, 1 B,j,k+1, with B1,j,k+1, u k+1,y k ;y k = B,j,k+1, u k+1,y k ;y k = and U k is the solution operator to the free equation, that is Ṽ, y j y k+1 u k+1,y k,y k+1 ;y k,y k+1 dy k+1, Ṽ, y j y k+1 u k+1,y k,y k+1 ;y k,y k+1 dy k+1. U k = e i y k e i y k. Then there is a C independent of j, k and such that T R k B,j,k+1, u k+1 where R k = L yj y. j d C k. Proof. We prove our main auxiliary theorem in Section 6. This theorem establishes the Klainerman- Machedon space-time bound for β, 7]. Via the above theorem, we infer that every limit point { u k} { } of satisfies the space time bound tan ωt ω for some C > and all 1 j k. R k B j,k+1 u k+1, ; u k d L R 3k R 3k Ck, Step III. Regarding the solution to the infinite hierarchy 4., we have the following uniqueness theorem. Theorem 6. Let { u k} be a solution of the infinite hierarchy 4. in [s,t] [ ], tanωt ω subject to zero initial data that is u k s,y k ;y k =, k, and the space time bound T s R k B j,k+1 u k+1, ; d L R 3k R 3k Ck, 4.5 for some C > and all 1 j k. Then k, [s,t], we have R k u k, ; =. L R 3k R 3k Proof. See Section 5. Since we have shown the space-time bound 4.5 in Step III, we apply the above uniqueness theorem and deduce that u k,y k ;y k = φ,y j φ,y j, 4.6 where φ,y solves the 3D LS 1 b φ i φ = φ y φ+ in R g φ,y = φ. { } Hence the compact sequence has only one limit point. So u k u k k φ,y j φ,y j as in the weak* topology. Since u k is an orthogonal projection, the convergence in the weak* topology is then equivalent to the convergence in the trace norm.

13 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 13 Example 1. At the suggestion of Professor Walter Strauss, we give a brief explanation on why a factorized state like formula 4.6 is a solution to the Gross-Pitaevskii hierarchy. Consider k = 1, then plugging φ,y 1 φ,y 1 into the infinite hierarchy yields φ,y1 φ,y 1 i + 1 y 1 1 y 1 φ φ,y 1 φ,y 1 = b b φ,y1 g = b B 1 1, φ φ g φ,y1 φ,y φ,y 1 φ,y b B1, g,y 1 which is i + 1 y 1 1 y 1 u 1 = 1 g b B 1, u. Step IV. In Step III, we have concluded the convergence lim Tr,y k;y k φ,y j φ,y j =, uk otice that u k = T 1 k γk and the lens transform of uk is where φt, x solves i t φ = γ k t,x k ;x k = φ,x = φ x φt,x j φt,x j, 1 x +ω x φ,y1 φ,y φ,y 1 φ,y φ+b φ φ in R 3+1 [, tanωt ]. ω i.e. equation 1.5. Thence we conclude that lim Tr k γk t,x k;x k φt,x j φt,x j =, t [,T ], as a result of the fact that the lens transform preserves the trace norm Lemma 4. Since equation 1. is time translation invariant, we have established Theorem 4 and consequently Theorem 1. The purpose of the rest of this paper is to prove Theorems 6 and 5. Remark 1. Since φ = M 1 1 φ, the global H1 well-posedness of equation 4.7 is implied by the global wellposedness in the scattering space of equation 1.5 which comes from the Strichartz estimates. Through the lens transform, we always have a L solution to equation 4.7. Lemma 6 then shows φ L. 5. The Uniqueness of Hierarchy 4. Proof of Theorem 6 In this section, we produce Theorem 6 with Lemmas 7 and 8. For convenience, we set the coupling constant b in the infinite hierarchy 4. to be 1. Lemma 7. [3] Assume u k+1 verifies i + 1 y k+1 1 y k+1 u k+1 =, then there is a C >, independent of j, k, and u k+1 s.t. R k B j,k+1 u k+1,y k ;y k L R R 3k R 3k CR k+1 u k+1,y k+1 ;y k+1 L. R 3k+1 R 3k+1, Proof. This is Theorem 1.3 of [3]. For some other estimates of this type, see [1, 1,, 8].

14 14 XUWE CHE Lemma 8. Assuming zero initial data i.e. u k s,y k ;y k =, k, then one can express u 1 1, ; in the Gross-Pitaevskii hierarchy 4. as a sum of at most 4 n terms of the form J n+1,µ m u n+1 n+1 d n+1, D or in other words, u 1 1, ; = J n+1,µ m u n+1 n+1 d n m D Here n+1 =, 3,..., n+1, D [s, 1 ] n, µ m are a set of maps from {,...,n+1} to {1,...,n} satisfying µ m = 1 and µ m j < j for all j, and 1 n J n+1,µ m u n+1 n+1 = g j+1 U 1 1 B 1, U 3 B µm 3,... U n n n+1 B µm n+1,n+1u n+1 n+1, ;. Proof. The RHS of formula 5.1 is in fact an application of Duhamel s principle involving only the inhomogeneous terms since we have zero initial data. The parameter n is the coupling level we take. This lemma follows from the proof of Theorem 3.4 in [3]. One needs only notice that factors depending solely on, e.g. commutes with U k and B j,k+1 j,k. 1 g j+1 With Lemmas 7 and 8, we prove Theorem 6. Let D = { 3,..., n+1, 3,..., n+1 D} where D is as in Lemma 8. Given that we have already checked that R 1 u 1 s, =, L R 3 R 3 applying Lemma 8 to [s, 1 ] [s,t] [, tanωt ω ], we have R 1 u 1 1, m m m = m m R1 [s, 1] n D inf [, tan ωt ω ] C n L R 3 R 3 L J n+1,µ m u n+1 n+1 d n+1 R 3 R 3 R 1 J n+1,µ m u n+1 L n+1 d R 3 R 3 n+1 [s, 1] n C n 1 s 1 g C n C 1 s 1 m Iterate n times... m n [s, 1] n R 1 B 1, U 3 B µm 3,... [s, 1] n 1 C n C 1 s n 1 CC 1 s n 1. [s, 1] n 1 1 s R 1 U 1 1 B 1, U 3 B µm 3,... d L R 3 R 3 n+1 d L R 3 R 3 n+1 R 1 B 1,U 3 B µm3,... 1 d d3...d n+1 L R 3 R 3 R B µm 3,U d n+1 Lemma 7 L R 3 R d 3 R n B µm n+1,n+1u n+1 n+1, ; d n+1 L R 3 R 3 Let 1 s be sufficiently small, and n, we infer that R 1 u 1 1, = in [s, 1 ]. L R 3 R 3

15 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 15 Such a choice of 1 s works for all of [s,t]. Accordingly, we have R k u k, L R 3 R 3 =, k, [s,t] by iterating the above argument. Thence we have attained Theorem The Space-Time Bound of the BBGKY Hierarchy for β,/7] Proof of Theorem 5 We establish Theorem 5 in this section. This section also serves as a simplification and an extension of Chen-Pavlovic [7]. Without loss of generality, we may assume k = 1 that is T R 1 B,1,, u d C. 6.1 L We are going to prove estimate 6.1 for a sufficiently small time T determined by the controlling constant in condition 4.4 and independent of, then the bootstrapping argument in Section 5 Proof of Theorem 6 and condition 4.4 provide the bound for every finite time T [, tanωt ω ]. Since we work with L norms here, we transform condition 4.4 into the H 1 energy bound: yj 1 y u k j dy k dy k Tr 1 yj u k. To obtain the above estimate, one notices 1 y 1 1 y 1 φy,rφy,rdr dydy = 1 y 1 φy,r1 y 1 φy,rdr dydy 1 y 1 φy,r1 y 1 φy,rdr 1 y 1 φy,r1 y 1 φy,rdr dydy = φy,r1 y φy,rdydr, the energybound then follows fromthe definition ofu k. The analysisoftheorem5alsoinvolves B,j,k+1, which approximates B j,k+1 for every, so we generalize Lemma 7 to the following collapsing estimate. Theorem 7. Suppose u,y 1,y,y solves the Schrödinger equation iu + 1 y 1 u+ 1 y u 1 y u = in R u,y 1,y,y = fy 1,y,y, and g c >, then there is a C independent of and u such that g y 3 V gy y δy y u,y,y ;y dy dy dyd R 3+1 Cb y1 y y f. where b = V dy. Proof. Theorem 7 follows from a slightly modified version of the proof of Theorem of [1]. We include it in Appendix II for completeness. We now present the proof of estimate 6.1. In order to more conveniently apply the Klainerman-Machedon board game, let us start by rewriting hierarchy 4.3 as where u k k = U k k u k, + k + k Ṽ k, uk,y k;y k = 1 B k+1, u k+1 = k k 1 i<j k Ṽ k U k, k k+1 k+1 u k g k+1 U k k k+1 k+1 B k+1, k+1 u k+1 k+1 g k+1 d k+1 Ṽ, y i y j Ṽ,y i y j u k,y k;y k B,j,k+1, u k+1 = k d k B1,j,k+1, B,j,k+1, u k+1.

16 16 XUWE CHE We omit the imaginary unit in front of the potential term and the interaction term so that we do not need to keep track of its exact power. Iterate Duhamel s principle equation 6.3 k times, we have u = U u, + = U u, Ṽ U, 3 3 u g 3 U 3 g 3 3 d 3 + B 3 U, 3 3 U 3 3 u 3, d 3 + g 3 B 3, 3 3 Ṽ 3 U 3, u 3 4 g 4 B 3 3, 3 U 3 g 3 = FreePart k +PotentialPart k +InteractionPart k where FreePart k k j = U u, + j=3 j 1 j 1... l= U j j u j, l=3 +1 l B 3 U, 3 3 u 3 3 d 3 g 3 d 4 d 3 B 4 U 3, u 4 4 d 4 d 3 g 4 Ṽ U, 3 3 u 3 d 3 g 3 1 g l+1 U 3 B 3, 3...U j 1 j 1 j B j, j d 3...d j, = PotentialPart k Ṽ U, 3 3 u 3 d 3 + g 3 j j 1 j 1... l= k j j=3 l=3 +1 l 1 g l+1 U 3 B 3, 3...U j 1 j 1 j B j, j Ṽ j U j, j j+1 j+1 u j g j+1 j+1 d j+1 d 3...d j, = InteractionPart k k+1 +1 l l=3 k k... g l+1 l= u k+1 k+1 d 3...d k+1 1 U 3 B 3, 3...U k k k+1 B k+1, k+1 From here on out, the k s in the formulas are the number of Duhamel iterations we take to prove estimate 6.1. We will call it the coupling level for short. It is distinct from the k in the statement of Theorem 5. We are going to argue T R 1 B,1,, FreePart k d C 6.4 T T R 1 B,1,, PotentialPart k R 1 B,1,, InteractionPart k L L L d C 6.5 d C 6.6

17 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 17 for some C and a sufficiently small T determined by the controlling constant in condition 4.4 and independent j of. We observe that B, j has j terms inside so that each summand of u contains factorially many terms k!. So we use the Klainerman-Machedon board game to reduce the number of terms. Define j 1 J j+1 f = g l+1 U 3 B 3, 3...U j j j+1 B j+1, j+1 f, l= where j+1 means 3,..., j+1, then the Klainerman-Machedon board game implies the lemma. Lemma 9. [3] One can express j... J j+1 fd j+1 as a sum of at most 4 j 1 terms of the form J j+1,µ m fd j+1, or in other words, D j... J j+1 fd j+1 = m D J j+1,µ m fd j+1. Here D [, ] j 1, µ m are a set of maps from {3,...,j+1} to {,...,j} satisfying µ m 3 = and µ m l < l for all l, and j 1 J j+1,µ m f = g l+1 U 3 B,,3,3 U B,µm 4,4,... 4 l= U j j j+1 B,µm j+1,j+1, j+1 f. Remark. There is no difference between Lemma 9 and the one we used for the uniqueness of hierarchy 4. Lemma 8. We have restated it to remind the reader of its exact form since we start from u here. With the above lemma and the collapsing estimate Theorem 7, we have the following relation, which is essentially part of the proof of Theorem 6, to help establishing estimates 6.4, 6.5, and 6.6. T L R1 B,1,, j+1,µ m fd j+1 d 6.7 DJ T j 1 = g l+1 R 1 B,1,, U 3 B,,3,3...d 3...d j+1 d D l= L j+1 inf g R 1 B,1,, U 3 B,,3,3... d d 3...d j+1 [, tan ωt ω ] [,T] j L C j 1 T 1 R 1 B,1,,U 3 B,,3, d d3...d j+1 L [,T] j 1 Cauchy-Schwarz C j 1 CT 1 [,T] j 1 Iterate j times... CT 1 j 1 T R B,,3,3 U d 3...d j+1 Theorem 7 R j B,µm j+1,j+1, j+1 f L d j+1. We show estimate 6.4 in Section Assuming Proposition 3, whose proof is postponed to Section 6..4, we derive estimate 6.5 in Section 6... Finally, by taking the coupling level k to be ln, we check estimate 6.6 in Section We remark that the proof of estimates 6.4 and 6.5 is independent of the choice of the coupling level k. Taking the coupling level k to be ln in the estimate of the interaction part is exactly the place where we follow the original idea of Chen and Pavlović[7]. Moreover, the proof of estimate 6.5 Section 6.. is the only place which relies on β, 7] in this paper.

18 18 XUWE CHE Estimate of the Free Part of u. Applying Lemma 9 and relation 6.7 to the free part of u, we obtain T R 1 B,1,, FreePart k d L CT 1 CT 1 CT 1 R u, R u, R u,, k + L + L j=3 m k j=3 +C L m T R1 B,1,, T CCT 1 j D k 4 j CT 1 j 1 R j u j j=3 J j,µ m U j j u j R j 1 B,µm j,j, j U j j u j, L CT 1 R u +C CT 1 j 1 C j Condition 4.4 L j=3 C < for T small enough. Whence, we have shown estimate Estimate of the Potential Part of u. We have T R 1 B,1,, PotentialPart k d L T + Ṽ R 1 B,1,, U, 3 3 u 3 d 3 d g 3 L k T j=3 m R1 B,1,, thus same procedure in Section 6..1 deduces, T CT 1 R Ṽ, 3 u 3 d 3 L k T + CCT 1 j CT 1 j=3 m T R Ṽ, 3 u 3 D j Ṽ j J j,µ m U j, t j t j+1 j+1 u j g j+1 Rj 1 B,µm j,j, j L d 3 +C j=3 j j+1 Ṽ j U j, j j+1 j+1 u j g j+1 k T 4 j CT 1 j 1 R j Ṽ j, j+1 u j Assume for the moment that we have the estimate T R k Ṽ k, uk L d C C k T where C and C are independent of T, k and, then T R 1 B,1,, PotentialPart k d C <. L, d j L, L d d j. d j+1 d j L j+1 d j+1 L j+1 L d, d j d j+1 for a sufficiently small T and a C independent of k and. As a result, we complete the proof of estimate 6.5 with the following proposition. Proposition 3. Assume β, 7] and V is a nonnegative L 1 R 3 H R 3 function. Then, given T [ ], tanωt ω, there are C and C independent of T, k and such that T R k Ṽ k, uk L d C C k T. Proof. The proof is elementary and we relegate it to Section Remark 3. This proposition is exactly the reason we restrict β, 7] in this paper..

19 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP Estimate of the Interaction Part of u. We proceed like Sections 6..1 and 6... T R 1 B,1,, InteractionPart k d L T R1 B,1,, J k+1,µ m u k+1 L k+1 d k+1 d m D T CCT 1 k 1 R k B,µm k+1,k+1, k+1 u k+1 k+1 d k+1. m L Then the next step is to investigate T R k B,µm k+1,k+1, k+1 u k+1 k+1 L d k+1. Without loss of generality, set µ m k +1 = 1 and look at B,1,k+1, 1 k+1, we have R k B1,1,k+1, k+1 u k+1 k+1 dy k dy k = R k Ṽ,k+1 y 1 y k+1 u k+1 k+1,y k,y k+1 ;y k,y k+1 dy k+1 dy k dy k C y1 Ṽ,k+1 y k+1 yj y j u k+1 k+1,y k,y k+1 ;y k,y k+1dy k+1 dy k dy k j= +C Ṽ,k+1 y 1 y k+1 R k u k+1 k+1,y k,y k+1 ;y k,y k+1dy k+1 dy k dy k = CI +II. oticing that Ṽ H C V H uniformly for [ ], tanωt ω, we can then estimate I y1 = Ṽ,k+1 y k+1 yj y j u k+1 k+1,y k,y k+1 ;y k,y k+1dy k+1 dy k dy k j= y1 dy k dy k y Ṽ,k+1 k+1 dy k+1 yj y j u k+1 k+1,y k,y k+1 ;y k,y k+1 dy k+1 Cauchy-Schwarz j= C 5β V L yj y j u k+1 k+1,y k,y k+1 ;y k,y k+1 dy k+1 dy k dy k j= C 5β V L dy k dy k yk+11 y yj k+1 y j u k+1 k+1,y k,y k+1 ;y k,y k+1 dy k+1 dy k+1 Trace Theorem C 5β V L Ck+1 Condition 4.4 j= and II = Ṽ,k+1 y 1 y k+1 R k u k+1 k+1,y k,y k+1 ;y k,y k+1dy k+1 dy k dy k C 3β V L Ck+1 Same method as I.

20 XUWE CHE Accordingly, R k B1,1,k+1, k+1 u k+1 k+1 dy k dy k C 5β V H Ck+1. Remark 4. The estimates of I and II may not be optimal. But they are good enough for proving estimate 6.6 for arbitrary β >. Thence T R 1 B,1,. InteractionPart k L d T CCT 1 k 1 R k B,µm k+1,k+1, k+1 u k+1 k+1 m 4 k 1 CCT 1 k 1 T C 5β V H C k+1 L d k+1. C V H T 1 k+ 5β C k. Take the coupling level k = ln, we have T Selecting T such that ensures that and thence R 1 B,1,. InteractionPart k d C V L H T 1 +ln 5β c. T T e 5β+C T 1 ln 5β c 1 R 1 B,1,. InteractionPart k L d C where C is independent of. We remind the reader that, at this point, we have obtained estimates 6.4, 6.5, and 6.6 for a sufficiently small T determined by the controlling constant in condition 4.4 and independent of. Thus one can repeat the argument to acquire the estimates for any finite time T [ ], tanωt ω through bootstrapping and condition 4.4. Whence, we have earned estimate 6.1 and established Theorem 5. The rest of this section is the proof of Proposition Proof of Proposition 3. We will utilize the lemma. Lemma 1. [18] Vx 1 x fx 1,x dx 1 dx C V L 1 1 x1 1 1 x 1 fx 1,x dx 1 dx. In particular, Vx 1 x fx 1,x dx 1 dx C V 1 x1 L 1 1 x 1 fx 1,x dx 1 dx. Proof. This is Lemma A.3 in [18]. Without loss of generality, we show Proposition 3 for k = which corresponds to T y1 y y y Ṽ, y 1 y 1 u,y ;y d C C T 6.8 L and T y1 y y y 1 Ṽ, y 1 y u,y ;y d C C T. L By similarity we only prove estimate 6.8. For a general k, there are k terms in Ṽ k,, whence Rk Ṽ k, γk has about 8k terms by Leibniz s rule. Since 8k can be absorbed into C k, the method here applies.

21 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 1 First of all, β, 7] implies the following properties of Ṽ, / : Ṽ,/ = 3β 1 L Ṽ < 1 7 Ṽ, L L / Ṽ, = 4β 1 3β p Ṽ = 6 7p pε 3 L Ṽ, decays up to p = 6, L p L p p / L Ṽ, 7 Ṽ ε. L Ṽ C V H H uniformly for [ ], tanωt ω. On the other hand we assume V is a nonnegative L 1 R 3 H R 3 function. Thus we know Ṽ L and Ṽ L 6. We compute y1 y y y 1 1 Ṽ, y 1 y u,y ;y y = y1 y 1 Ṽ, y 1 y 1 y u,y ;y + y1 1 Ṽ, y 1 y y y y 1 u,y ;y + y 1 Ṽ, y 1 y y1 y y 1 u,y ;y + 1 Ṽ, y 1 y y1 y y y 1 u,y ;y. where ε is 7 β. On the one hand, But and sup = sup sup y y1 y 1 Ṽ, y 1 y y 1 u,y ;y 1 y1 Ṽ, y y 1 y u,y ;y dy dy 1 y Ṽ, dy 1 y1 1 1 y 1 Lemma 1 C 7ε V L sup 1 y1 1 1 y 1 7ε V L C Condition 4.4 sup = sup sup L y y y 1 1 y u 1 Ṽ, y 1 y y1 y y y 1 u,y ;y 1 Ṽ,y 1 y y1 y y y 1 1 Ṽ, 1 y1 1 1 y 1 L 7 V L sup 1 y1 1 1 y 1 L u u,y ;y dy dy y 1 1 y u,y ;y dy dy,y ;y dy dy,y ;y dy dy y 1 1 y u,y ;y dy dy C 7 V H C Sobolev and Condition 4.4. Since the same method applies to the middle terms, we have obtained estimate 6.8 and hence Proposition Conclusion In this paper, we have rigorously derived the 3D cubic nonlinear Schrödinger equation with a quadratic trap from the -body linear Schrödinger equation. The main novelty is that we allow a quadratic trap in our analysis and the main technical improvements are the simplified proof of the Klainerman-Machedon space-time bound as compared to non-trap case in Chen and Pavlović [7], and the extension of the range of β from,1/4 in Chen and Pavlović [7] to,/7]. Compared to the D work [1] which is also by the author, the 3D problem in this paper is of critical regularity. To explain what we mean by critical: in D one easily obtains the 1 -space-time bound needed for the uniqueness theorem by a trace theorem; in 3D the only way to obtain the space-time bound 4.5 is through smoothing estimates since one does not have

22 XUWE CHE enough regularity to apply a trace theorem. Thence the key arguments in 3D are more involved and totally different from the D case which is a subcritical problem. Moreover, we have established the trace norm convergence in the main theorem which is a stronger result than the Hilbert-Schmidt norm convergence. 8. Acknowledgements The author would like to first thank his advisors, Professor Matei Machedon and Professor Manoussos G. Grillakis, for the discussion related to this work and their guidance on the author s PhD dissertation which consists of [11, 1]. It is due to their guidance, the author became keenly interested in the subject and wrote this paper. The author is greatly indebted to Professor Thomas Chen and Professor ataša Pavlović for telling the author about their important work [7] and giving the author many helpful suggestions on finding his first job during his visit to Austin, and for their very detailed and helpful comments on this paper. The author s thanks go to Professor Sergiu Klainerman for the discussion related to this work, and to Professor Clément Mouhot for sharing the physical meaning of the trace norm in this setting with the author. The author also would like to thank the anonymous referee for many insightful comments and helpful suggestions. 9. Appendix I: Proof of Corollary 1 For the purpose of this Appendix I, we may assume ω =. Or in other words, we skip Steps I and IV of the proof of Theorem 4 here. When the desired limit is an orthogonal projection, one does not need this proof. This is a functional analysis argument and all operators mentioned in this Appendix I acts on L R 3k. We prove Corollary 1 by verifying the hypothesis of the following lemma. Lemma 11. [35] Assume the operator sequence {A n } satisfies that, as bounded operators, A n A, A n A and A n A in the weak sense. If then lim Tr A n = Tr A, lim Tr A n A =. Proof. This is Theorem. in [35]. It implies the Grümm s convergence theorem Theorem.19 of [35] used in [1]. We first observe that condition 1.7 implies the a-priori estimate Tr 1 xj γ k Ck. sup t [,T] Thus we have the compactness argument and the uniqueness argument to conclude that, as trace class operator kernels, γ k t,x k;x k γ k t,x k ;x k weak*. 9.1 Remark 5. Because we assume ω =, both of the Erdös-Schlein-Yau uniqueness theorem [18] and the Klainerman-Machedon uniqueness theorem [3] apply here. For the general case, one has to use the main argument in this paper. Let H k be the Hilbert-Schmidt operators on L R 3k. Recall that the test functions for weak* convergence in L 1 k come from K k and the test functions for weak* convergence in H k come from H k. Thus the weak* convergence 9.1 as trace class operator kernels infers that as Hilbert-Schmidt kernels, γ k t,x k;x k γk t,x k ;x k weak*, because H k K k i.e. there are fewer test functions. Since H k is reflexive, the weak* convergence is no different from the weak convergence. Thus we know that as Hilbert-Schmidt kernels and hence as bounded operator kernels, γ k t,x k;x k γ k t,x k ;x k weak. At this point, we have verified that, as bounded operators, A n A and A n A in the weak sense since γ k and γk are self adjoint. We now check A n A.

23 SPACE-TIME BOUD DERIVATIO OF 3D LS WITH A QUADRATIC TRAP 3 To check A n A, one first notices that γ k and γk has only real eigenvalues since they are self adjoint. Moreover, γ k t,x k;x k has no negative eigenvalues by definition 1.3, in fact, φx,z φy,zdz fydy fxdx = dz φx,z fxdx. Since fx k fx k in the above estimate is also a Hilbert-Schmidt kernel, the definition of weak convergence in L dx k dx k gives γ k t,x k ;x k fx k fx k dx kdx k = lim γ k t,x k;x k fx k fx k dx kdx k, as γ k t,x k;x k γk t,x k ;x k weakly in L dx k dx k. So we have checked A n A because A n = A n and A = A. To prove Corollary 1, by Lemma 11, it remains to argue otice that we have the conservation of trace γ k t,x k;x k dx k = γ k t,x k ;x k dx k = γ lim Tr k t,x k;x γ k = Tr k t,x k ;x. k γ k,x k;x k dx k By definition 1.3 γ k,x k ;x k dx k [8]. and we have shown that γ k t,x k;x k and γk t,x k ;x k have no negative eigenvalues, hence Trγ k t,x k;x k = γ k t,x k;x k dx k = γ k,x k;x k dx k, Trγ k t,x k ;x k = γ k t,x k ;x k dx k = γ k,x k ;x k dx k. On the one hand, γ k,x k ;x k dx k = Trγ k,x k ;x γ k = lim Tr k,x k;x γ k = lim Tr k t,x k;x k, on the other hand, That is Trγ k t,x k ;x k = γ k t,x k ;x k dx k = γ k,x k ;x k dx k. Trγ k t,x k ;x γ k = lim Tr k t,x k;x k. Whence we conclude the proof of Corollary 1 by Lemma Appendix II: Proof of Theorem 7 In this Appendix II, we prove Theorem 7. We will make use of the lemma. Lemma 1. [3] Let ξ R 3 and P be a d plane or sphere in R 3 with the usual induced surface measure ds. 1 Suppose < a, b <, a+b >, then dsη P ξ η a η b C ξ a+b. Suppose ε = 1 1, then P dsη C ξ ξ η η ε ε η Both constants in the above estimates are independent of P. ξ 3 ε. Proof. See pages of [3].

24 4 XUWE CHE In this appendix II, we write which has the property that,ξ, V x = g 3 V gx sup V ξ Vx dx = b. By duality, to gain Theorem 7, it suffices to prove that h,y y Vy y δy y u,y,y ;y dy dy dyd R 3+1 Cb h y1 y y f. From equation 6., we compute the spatial Fourier transform of y V y y δy y u,y,y ;y dy dy to be ξ 1 V ξ +ξ e i ξ 1 ξ ξ + ξ ξ ˆfξ1 ξ ξ,ξ,ξ dξ dξ, thus we have h,y y V y y δy y u,y,y ;y dy dy dyd R 3+1 = ξ 1 V ξ +ξ e i ξ 1 ξ ξ + ξ ξ ˆfξ1 ξ ξ,ξ,ξ ĥ,ξ 1ddξ 1 dξ dξ = spatial Fourier transform on h V ξ +ξ ξ 1 e i ξ 1 ξ ξ + ξ ξ ĥ,ξ 1 d ˆfξ 1 ξ ξ,ξ,ξ dξ 1 dξ dξ Ih y1 y y f L Cauchy-Schwarz where ξ1 V ξ +ξ e i ξ 1 ξ ξ + ξ ξ ĥ,ξ 1 d Ih = ξ1 ξ ξ ξ ξ dξ 1 dξ dξ. So our purpose in the remainder of this Appendix II is to show that Ih Cb h L. oticing that, away from the factor V, the integral Ih is symmetric in ξ1 ξ ξ and ξ, it suffices that we deal with the region: ξ 1 ξ ξ > ξ only since our proof treats V as a harmless factor. We separate this region into two parts, Cases I and II. Away from the region ξ1 ξ ξ > ξ, there are other restrictions on the integration regions in Cases I and II. We state the restrictions in the beginning of both Cases I and II. Due to the limited space near, we omit the actual region. The alert reader should bear this mind Case I: Ih restricted to the region ξ < ξ with integration order dξ prior to dξ. Write the phase function of the d integral inside Ih as ξ1 ξ ξ + ξ ξ1 ξ ξ = + ξ ξ 1 ξ ξ. The change of variable ξ,new = ξ,old ξ 1 ξ leads to the expression ξ 1 ξ1 V ξ + ξ 1 +ξ e i ξ + ξ ξ ĥ,ξ 1 d Ih = ξ ξ 1 ξ ξ + ξ 1 ξ dξ 1 dξ dξ. ξ 1.1