arxiv:math/ v1 [math.ap] 19 May 2004

Transcription

1 SEMICLASSICAL LIMIT OF THE HARTREE EQUATION WITH HARMONIC POTENTIAL RÉMI CARLES, NORBERT J. MAUSER, AND HANS PETER STIMMING arxiv:math/ v1 [math.ap] 19 May 004 Abstract. Nonlinear Schrödinger Equations NLS of the Hartree type occur in the modeling of quantum semiconductor devices. Their semiclassical limit of vanishing scaled Planck constant is both a mathematical challenge and practically relevant when coupling quantum models to classical models. With the aim of describing the semi-classical limit of the 3D Schrödinger Poisson system with an additional harmonic potential, we study some semi-classical limits of the Hartree equation with harmonic potential in space dimension n. The harmonic potential is confining, and causes focusing periodically in time. We prove asymptotics in several cases, showing different possible nonlinear phenomena according to the interplay of the size of the initial data and the power of the Hartree potential. In the case of the 3D Schrödinger Poisson system with harmonic potential, we can only give a formal computation since the need of modified scattering operators for this long range scattering case goes beyond current theory. We also deal with the case of an additional local nonlinearity given by a power of the local density - a model that is relevant when incorporating the Pauli principle in the simplest model given by the Schrödinger-Poisson-Xα equation. Further we discuss the connection of our WKB based analysis to the Wigner function approach to semiclassical limits. 1. Introduction Nonlinear Schrödinger Equations NLS are important both for many different applications as well as a source of rich mathematical theory with several hard challenges still open. The NLS in the most common meaning contains a local nonlinearity given by a power of the local density, in particular the defocusing cubic NLS which arises e.g. in nonlinear optics or for Bose Einstein condensates. In 1-d this NLS is an integrable system and the semi-classical limit high wave number limit can be performed by methods of inverse scatteringsee e.g. [0] and [] for results on the defocusing and focusing case. A class of NLS with a nonlocal nonlinearity that we call Hartree type occur in the modeling of quantum semiconductor devices. Their semi-classical limit of vanishing scaled Planck constant is both a mathematical challenge and practically relevant when coupling quantum models to classical models. Incorporating the Pauli principle for fermions in the simplest possible model yields the case of a Hartree equation with an additional local nonlinearity given by a power of the local density, the Schrödinger-Poisson-Xα equation see [5]. In this paper we deal with the semi-classical limit of nonlinear Schrödinger equations of Hartree type, with a harmonic potential and a weak nonlinearity which is a convolution of the density with a more or less singular potential. 000 Mathematics Subject Classification. 35B33, 35B40, 35C0, 35Q40, 81Q0, 81S30. The authors acknowledge support by the Austrian START award project FWF, contract No. Y-137-TEC of N.J.M. and by the Wissenschaftskolleg doctoral school Differential Equations FWF, contract No. W8 as well as the European network HYKE, funded by the EC as contract HPRN-CT

2 R. CARLES, N. MAUSER, AND H. P. STIMMING In three space dimensions, for the case where we convolute with the Newtonian potential 1/ x, the Hartree equation is the Schrödinger Poisson system with harmonic potential : iε t u ε + 1 ε u ε = x uε +Vxu ε, 1.1 V = u ε, u ε t=0 = uε 0, with x R 3. This equation arises typically if we consider the quantum mechanical time evolution of electrons in the mean field approximation of the many body effects, modeled by the Poisson equation, with a confinement modeled by the quadratic potential of the harmonic oscillator. The limit ε 0 in such a quantum model corresponds to a classical limit of vanishing Planck constant = ε 0. We adopt the terminology semi-classical limit for what should properly be called classical limit see the discussion in [31], the term semi-classical being actually more appropriate for the situation of the homogenization limit from a Schrödinger equation with periodic potential see e.g. []. The problem of the mathematically rigorous classical limit of the Schrödinger- Poissonsystemis highly nontrivial. First resultsofweaklimits ε 0to the Vlasov- Poisson system where given in [3] and [4] using Wigner transform techniques for the mixed state case, where additional strong assumptions on the initial data can be imposed which are necessary to guarantee a uniform bound on the Wigner function. In [31] this assumption could be removed for the 1-d case and the classical limit for the pure state case could be performed, where the notorious problem of non-uniqueness of the Vlasov-Poisson system with measure valued initial data reappears. For an overview of this kind of semi-classical limits of Hartree equations see [6]. For an introduction to Wigner transforms and their comparison to WKB methods for the linear case see [11] and [9]. Up to a constant, 1.1 is equivalent to the Hartree equation 1. iε t u ε + 1 ε u ε = x uε + x 1 u ε u ε ; u ε t=0 = uε 0. We restrict our attention to small data cases with u ε 0 = εα/ f, where f is independent of ε and α 1. Notice that we can allow for more general data with initial plane oscillations, 1.3 u ε t=0 = εα/ fxe i x.ξ 0 ε for ξ 0 R 3, since the change of variables given in [6] 1.4 u ε t,x = u ε t,x ξ 0 sinte ix ξ 0 sint.ξ 0 cost/ε, yields the solution of 1.. This change of variable could also be used in Equation 1.6 below and hence our results also hold for the more general ε-dependent class of data 1.3. Note that small data can be equivalently written as small nonlinearity, since with the change of the unknown u ε = ε α/ u ε, 1. becomes 1.5 iε t u ε + 1 ε u ε = x uε +ε α x 1 u ε u ε ; u ε t=0 = f. We will consider the more general semi-classical Hartree equation 1.6 iε t u ε + 1 ε u ε = x uε +ε α x γ u ε u ε ; u ε t=0 = f,

3 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 3 with γ > 0, α 1 and x R n, where the space dimension n may be different from 3. The first point to notice is that in the linear case, the harmonic potential causes focusing at the origin resp. at 1 k ξ 0 in the case 1.4 at times t = π/+kπ, for any k N. The solution u ε free of the linear equation 1.7 iε t u ε free + 1 ε u ε free = x uε free ; u ε free t=0 = f, is initially of size Ø1. At time t = π/, the solution focuses at the origin and is of order Øε n/ ; it is of order Ø1 for t = π, and so on for a more precise analysis, see [6]. This phenomenon is easy to read from Mehler s formula see e.g. [10, 18]: for 0 < t < π, we have 1.8 u ε free t,x = e inπ 4 πεsint n/ e i εsint R n x +y cos t x y fydy. Essentially, one can apply a stationary phase formula for t ]0,π/[ ]π/,π[ u ε free is Ø1, while it is not possible at t = π/ u ε free is Øε n/. Following the same approach as in [3], we get the following distinctions: α > γ α = γ α > 1 Linear WKB, Linear WKB, linear focus nonlinear focus α = 1 Nonlinear WKB, Nonlinear WKB, linear focus nonlinear focus The expression linear WKB means that the nonlinear Hartree interaction term is negligible away from the focus when the WKB approximation is valid; linear focus means that the nonlinearity is negligible near the focus; the WKB régime resp. the focus is nonlinear when the Hartree term has a leading order influence away from resp. in the neighborhood of the focus, in the limit ε 0. This terminology follows [19]. We did not obtain a rigorous description of the case α = γ = 1, which corresponds to the Schrödinger Poisson system 1.1 when n = 3. This problem seems out of reach for the methods currently available in this field. On the other hand, we study rigorously the three other cases in an exhaustive way: In Section 3, we prove that the Hartree term has no influence at leading order when α > γ = 1. Back to 1., this shows that initial data of size ε α/ with α > 1 yield a linearizable solution. The expected critical size is ε; this heuristic is reinforced by the next three sections. In Section 4, we study the case α = 1 > γ. We prove that the nonlinear term must be taken into account to describe the solution u ε. It is so through a slowly oscillating phase term. On the other hand, no nonlinear effect occurs at leading order near the focus. In Section 5, we show that when α = γ > 1, nonlinear effects occur at leading order at the focuses, while they are negligible elsewhere. This phenomenon is the same as in [6] for the nonlinear Schrödinger equation; each focus crossing is described in terms of the scattering operator associated to the Hartree equation 1.9 i t ψ + 1 ψ = x γ ψ ψ. In Section 6, we perform a formal computation suggested by the results of Sections 4 and 5. This can be seen as a further evidence that nonlinear effects are always relevant in the case α = γ = 1, along with a precise idea of the nature of these nonlinear effects, which we expect to be true. We add a brief discussion of

4 4 R. CARLES, N. MAUSER, AND H. P. STIMMING the case of an additional local nonlinearity in the equation and some remarks on the Wigner measures in view of the ill-posedness results of [5]. This program is very similar to the one achieved in [3]. We want to underscore at least two important differences. First, we have to adapt the notion of oscillatory integral to incorporate the presence of the harmonic potential see Section 4.1. Second, the power-like nonlinearity treated in [3] is replaced by a Hartree-type nonlinearity. This yields different and less technical proofswe do not use Strichartz estimates in Sections 3 and 4, and makes a more complete description of the above table possible; the case nonlinear WKB, linear focus was treated very partially in [3], due to the lack of regularity of the map z z σ z for small σ > 0. This technical difficulty does not occur in the present case, and the main result of Section 4 Proposition 4.1 is proved with no restriction. The content of this article is as explained above, plus a paragraph dedicated to a quick review of the facts we will need about the Cauchy problem 1.6 Section. We will use the following notation throughout this paper. Notation. If a ε ε ]0,1] and b ε ε ]0,1] are two families of numbers, we write a ε b ε if there exists C independent of ε ]0,1] such that for any ε ]0,1], a ε Cb ε.. The Cauchy problem Before studying semi-classical limits, we recall some known facts about the initial value problem 1.6. We will always assume that the initial datum f is in the space Σ defined by Σ := { φ H 1 R n ; φ Σ := φ + xφ + φ < + }. This space is natural in the case of Schrödinger equations with harmonic potential, since Σ is the domain of + x see for instance [7]. Local existence results for 1.6 follow for instance from Strichartz inequalities one can do without these inequalities, see [7]. Global existence results then stem from conservation laws see.3 below. From Mehler s formula 1.8, Strichartz type estimates are available for e i t ε ε +x =: U ε t. Definition. Let n. A pair q,r is admissible if r < n n resp. r < if n = and 1 q = δr n 1. r Following [6], we have the following scaled Strichartz inequalities: Proposition.1. Let I be a finite time interval. 1 For any admissible pair q,r, there exists C r I such that.1 ε 1 q U ε tφ L q I;L r C ri φ. For any admissible pairs q 1,r 1 and q,r, there exists C r1,r I such that. ε 1 q q U ε t sfsds C r1,r I F L q I;L r. L q 1I;L r 1 I {s t} The above constants are independent of ε. The main result of this section follows from [7, 1]. Denote YI = {φ CI,Σ ; φ, x φ, x φ L q loc I,Lr x, q,r admissible}.

5 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 5 Proposition.. Fix ε ]0,1] and let f Σ. Then 1.6 has a unique solution u ε YR. Moreover, the following quantities are independent of time:.3 Mass: u ε t, Energy: 1 ε xu ε t + 1 xuε t +ε R n x γ u ε u ε t,x dx It was noticed in [6] that this result can be retrieved very simply thanks to the following lemma, which we will use to prove asymptotics. Lemma.3 [6]. Define the operators.4 J ε t = x ε sint icost x ; H ε t = xcost+iεsint x. J ε and H ε satisfy the following properties. They are Heisenberg observables:.5 J ε t = iu ε t x U ε t ; H ε t = U ε txu ε t. The commutation relation: ].6 [J ε t,iε t + ε x = ] [H ε t,iε t + ε x = 0. Denote M ε t = e ix ε tant, and Q ε t = e ix ε cott, then.7 J ε t = icostm ε t x M ε t ; H ε t = iεsintq ε t x Q ε t. The modified Sobolev inequalities. Let r n exists C r independent of ε such that, for any φ Σ,.8 n φ L r C r cost δr φ 1 δr L J ε tφ δr L, φ L r C r εsint δr φ 1 δr L H ε tφ δr L. Action on nonlinear Hartree term: for φ = φt,x, r < if n = ; there.9 J ε t x γ φ φ = x γ φ J ε tφ+re x γ φj ε tφ φ. The same holds for H ε t. Remark. Property.6 follows from.5, which is the way J ε and H ε appear in the linear theory see e.g. [30, p. 108]. Property.8 is a consequence of Gagliardo Nirenberg inequalities and.7. Finally,.9 stems from Very weak nonlinearity case In this section, we study the semi-classical limit of u ε when γ = 1 and α > 1 which is equivalent to very small data in our context cf. 1.. This case includes the 3D Schrödinger Poisson equation with very small data. We prove that the Hartree term plays no role at leading order. Proposition 3.1. Let f Σ, n, and assume α > γ = 1. Then for any T > 0, u ε u ε free L [0,T]; = Ø ε α 1 ln 1 as ε 0, ε and for any δ > 0 α 1 δ > 0, A ε tu ε u ε free L [0,T]; = Ø ε α 1 δ as ε 0, where A ε is either of the operators J ε or H ε, and u ε free is the solution of 1.7. Remark. Using modified Sobolev inequalities.8, we can deduce L p estimates for u ε u ε free for p n/n p < if n = from the above result.

6 6 R. CARLES, N. MAUSER, AND H. P. STIMMING Remark. We could probably get the logarithmic estimate for the second part of the statement as well, using Strichartz estimates. The proof given below is not technically involved, and suffices for our purpose: we do not seek sharp results. Proof. Denote w ε = u ε u ε free. It solves the initial value problem iε t w ε + 1 ε w ε = x wε +ε α x 1 u ε u ε ; w ε t=0 = 0. Standard energy estimates for Schrödinger equations yield 3.1 ε t w ε t ε α x 1 u ε u ε. From Hölder s inequality, we have 3. x 1 u ε u ε x 1 u ε L r u ε L k, for 1 r + 1 k = 1. From the Hardy Littlewood Sobolev inequality, 3.3 x 1 u ε t L r u ε t L, for 1 < r, p p < and 1+ 1 r = p + 1 n. Therefore, 3.1 yields 3.4 ε t w ε t ε α u ε t L p uε t L k, where p and k satisfy the properties stated in 3. and 3.3. For k =, r = and p = n/n 1, the algebraic identities stated in 3. and 3.3 are satisfied. Now since the conditions 1 < r < and 1 < p/ < are open, a continuity argument shows that we can find p and k satisfying all the properties stated in 3. and 3.3. Notice that they imply the relation δp +δk = 1, hence δp, δk < 1; this allows us to use weighted Gagliardo Nirenberg inequalities. We have w ε t=0 = 0, and from Proposition., wε CR + ;Σ. Therefore, there exists t ε > 0 such that 3.5 J ε tw ε 1, for 0 t t ε. The argument of the proof then follows [8] see also [6]. Recall that from Lemma.3, J ε tu ε free = f. Because of.6, J ε u ε free solves the linear Schrödinger equation with harmonic potential, and J ε tu ε free f L. So long as 3.5 holds, we have, from.8, u ε C 0 t L p ; u ε C 0 t cost δp L k cost, δk for some C 0 independent of ε and t. Then 3.4 yields ε t w ε t ε α εα = cost δp+δk cost. Integration in time on [0,t] yields, so long as 3.5 holds, w ε L [0,t]; ε α 1 t For t < π/, we get, so long as 3.5 holds, From.6, J ε tw ε solves 0 dτ cosτ, w ε L [0,t]; ε α 1 ln π t. iε t J ε w ε + 1 ε J ε w ε = x Jε w ε +ε α J ε x 1 u ε u ε ; J ε w ε t=0 = 0.

7 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 7 Using.9, energy estimate for J ε w ε yields ε t J ε tw ε ε α x 1 u ε J ε tu ε + x 1 u ε J ε u ε u ε ε α x 1 u ε L J ε tu ε + x 1 u ε J ε u ε u ε For the first term of the right hand side, use the easy estimate x 1 f f L n + f L n+ where n res. n + stands for n η resp. n+η for any small η > 0. We have x 1 u ε L u ε t L κ + uε t n Lκ+, with κ = n 1 It is at this stage that we lose the logarithmic ratewe cannot use Hardy Littlewood Sobolev inequality when an exponent is infinite: using Strichartz estimates see Section 5, we believe that we could recover that rate, with a more technically involved proof. For the second term, we proceed as in the beginning of the proof. From Hölder s inequality, 3.6 x 1 u ε J ε u ε u ε x 1 u ε J ε u ε L r u ε L σ, with 1 r + 1 σ = 1 From the Hardy Littlewood Sobolev inequality, this is estimated, up to a constant, by 3.7 u ε J ε u ε L p u ε L σ, with 1+ 1 r = 1 p + 1 n for 1 < r,p <. Using Hölder s inequality again yields an estimate by 3.8 u ε L k J ε u ε u ε L σ, with 1 p = k Take r = n, σ = n/n, k = and p = 1: the algebraic identities from 3.6, 3.7and3.8aresatisfied,butnottheboundp > 1. Decreasingslightlyσ increases p take σ large but finite when n =, so we can find indices satisfying 3.6, 3.7 and 3.8 by a continuity argument. Note that they satisfy δk + δσ = 1, and each term is positive. Gathering all these estimates together yields the energy estimate ε t J ε tw ε ε α u ε t L κ + uε t L κ+ + uε L k u ε L σ J ε tu ε So long as 3.5 holds, we deduce from.8, ε t J ε tw ε ε α 1 cost + 1 δκ cost + 1 δκ+ ε α 1 cost + 1 ε α. δκ+ cost cost 1+ Integrate this, so long as 3.5 holds: J ε w ε L [0,t]; ε α 1 π 0 t Fix δ,λ > 0. So long as 3.5 holds, we infer, for t π/ Λε, J ε w ε L [0,t]; ε α 1 Λε δ. cost δk+δσ Therefore, there exists ε Λ > 0 such that, for 0 < ε ε Λ, 3.5 holds up to time π/ Λε, with the estimates 3.9 w ε L [0,π/ Λε]; ε α 1 ln 1 ε ; J ε w ε L [0,π/ Λε]; ε α 1 δ.

8 8 R. CARLES, N. MAUSER, AND H. P. STIMMING An estimate similar to that of J ε w ε then follows for H ε w ε, since from.3, H ε tu ε f Σ. Denote IΛ ε = [π/ Λε,π/+Λε]. Mimicking the above computations, we have w ε L IΛ ε;l w ε π Λε L +ε α 1 u ε τ L p uε τ L kdτ, IΛ ε where p and k satisfy 3. and 3.3. Recall that they satisfy δp + δk = 1. Using the conservations of mass and energy.3, along with Gagliardo Nirenberg inequalities, we have, for any t, We deduce w ε L I ε Λ ;L u ε t L p ε δp ; u ε t L k ε δk. w ε π Λε L +ε α 1 ε δp δk IΛ ε ε α 1 ln 1 ε +Λεα 1. The same method yields, since.3 shows that H ε tu ε f Σ : To treat the case of J ε w ε, introduce H ε w ε L I ε Λ ;L ε α 1 δ, for any δ > 0. z ε t = sup J ε τw ε π. Λε τ t Proceeding as above, we have z ε t J ε π t w Λε ε L +ε α 1 π 3.10 Λε t ε α 1+ +ε α 1 z ε τ+ J ε τu ε free dτ. ε 1 + π Λε J ε τ x 1 u ε u ε dτ We can then apply the Gronwall lemma recall that J ε τu ε free f : z ε t ε α 1+. Gathering these informations we get, for any δ > 0: w ε L I ε Λ ;L ε α 1 ln 1 ε, J ε w ε L I ε Λ ;L + H ε w ε L I ε Λ ;L ε α 1 δ. For t [π/+ε,π], we can use the same proof as for t [0,π/ ε], to obtain: w ε L [0,π]; ε α 1 ln 1 ε, J ε w ε L [0,π]; + H ε w ε L [0,π]; ε α 1 δ. Repeating the same argument a finite number of times covers any given time interval [0,T] and completes the proof of Proposition Nonlinear propagation and linear focus In this paragraph, we assume α = 1 and γ < 1. We define 4.1 gt,x = x γ f x t 0 dτ cosτ γ. This function is well defined for any t, since γ < 1. We will see later on how this function appears.

9 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 9 Proposition 4.1. Let n, f Σ, and assume γ < α = 1. Let A ε be either of the operators Id, J ε or H ε. For 0 t < π/, the following asymptotics holds: sup 0 τ t Aε τ For π/ < t π, sup Aε τ t τ π u ε τ,x 1 cosτ n/f x cosτ u ε τ,x e inπ cosτ n/f x cosτ e ix ε tanτ+igτ, x cosτ e ix ε tanτ+igτ, x cosτ x x 0. ε 0 0. ε 0 For t = π/, Bε u ε π 1 ε n/f fe igπ ε 0, ε 0 where B ε is either of the operators Id, x ε or ε x, and the Fourier transform is defined by 4. Fφξ = φξ = 1 e ix ξ φxdx. π n/ R n Remark. We can also prove estimates for arbitrarily large time intervals, with the same proof as below. Remark. The difference between the asymptotics before and after the focus is measured only by the Maslov index, through the phase shift e inπ/ : no nonlinear phenomenon occurs at leading order near the focus. On the other hand, nonlinear effects are relevant outside the focus, as shown by the presence of g Oscillatory integrals. The main tool for proving Proposition 4.1 is the same as in linear cases [9], see also [1, 3] for applications in nonlinear settings: we represent the solution u ε as an oscillatory integral. Recall that u ε CR;Σ and that e i t ε ε +x = U ε t is a unitary group on. Define a ε by 4.3 a ε t,x = U ε tu ε t,x. We first seek a limit as ε 0 for a ε before the focus. This is suggested by a formal computation as in [4], and the following lemma: Lemma 4.. For t [0,π/[ ]π/,π], define V ε by 1 x cost n/φ cost 4.4 V ε tφx = e inπ/ cost n/φ x cost e ix ε tant if 0 t < π/, e ix ε tant if π/ < t π. For any φ H 1 R n, any θ ]0,1/], and any t [0,π/[ ]π/,π], U ε tφ V ε tφ εtant θ φ H 1. Proof. Notice that from Mehler s formula 1.8, we can write, for 0 < t < π, U ε t = M ε t Dε t FMε t where M ε t x = e i x εtan t, D ε t φx = 1 x iεsint n/φ sint and the Fourier transform is defined by 4.. We infer U ε tφ V ε tφ = 1 e i x y L iπtant n/ ε tan t fydy fx From Parseval formula, 1 e i x y 1 iπtant n/ εtan t fydy = e iεtantξ π n/ +ix ξ Ffξdξ,,

10 10 R. CARLES, N. MAUSER, AND H. P. STIMMING therefore U ε tφ V ε 1 tφ = π n/ e iεtantξ 1 e ix ξ Ffξdξ = e iεtantξ 1 Ffξ, fromplancherelformula. Thelemmathenfollowsfromtheestimate e is 1 s θ, for 0 θ 1/. From Duhamel s principle, we have u ε t = U ε tf i Using 4.3, we deduce t 0 U ε t s x γ u ε u ε sds. 4.5 t a ε t = iu ε t x γ u ε u ε t. Now the formal computation begins. Assume a ε a as ε 0, in some suitable sense. Then u ε t U ε tat, and from Lemma 4., u ε 1 x t,x ε 0 cost n/a t, e ix ε tant for 0 t < π/. cost Plugging this into 4.5 and using Lemma 4. again with U ε t instead of U ε t, the result still holds yields t at,x = i cost γ x γ a at,x. Recall that a t=0 = u ε t=0 = f, and notice that from the above ordinary differential equation, t a = 0: we have at,x = fxe igt,x, where t gt,x = 1 cost γ x γ f x ; g t=0 = 0. Integrating this equation yields the definition of gt, x given in 4.1. Proposition 4.1 stems from the more precise following proposition, Lemma 4. and a density argument. In view of a rigorous justification, denote 4.6 b ε t,x = a ε t,xe igt,x = e igt,x U ε tu ε t,x. Proposition 4.3. Let f Σ H R n. Fix δ > 0. There exists C δ such that sup b ε t f Σ 0 t π π 0 t b ε t Σ dt C δ ε 1 γ δ. The first inequality is trivial. We prove the second one in three steps: i On [0,π/ Λε] for any Λ > 0, with a constant depending on δ and Λ. ii On [π/ Λε,π/+Λε], with a constant depending on δ and Λ. iii On [π/+λε,π], with a constant depending on δ and Λ. As in Section 3, the parameter Λ > 0 is arbitrary, while it has to be large in the case α = γ > 1 see Section 5 and [6]: this situation is typical from a case where the focus is linear.

11 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL Asymptotics before the focus. Fix Λ, δ > 0. We prove that there exists C Λ,δ such that 4.7 Denote π Λε 0 t b ε t Σ dt C Λ,δ ε 1 γ δ. y ε t = t 0 t b ε τ H 1dτ. From 4.5 and the definition 4.6, t b ε t = Uε t x γ u ε u ε 1 t x γ cost γ f a ε t = x γ u ε u ε 1 t cost γuε t x γ f a ε 4.8 t. Lemma 4. suggests that we can replace U ε with V ε in the last expression, up to a controllable error. Before going further into details, we prove two lemmas which will be of constant use in the proof of Proposition 4.3. Lemma 4.4. Assume γ < 1, and let 0 < δ < 1 γ. There exist p and q with δp = γ + δ < 1, p < n γ ; δq = γ +1 + δ 4 and such that there exists C such that for any φ C c Rn, x γ φ L C φ L 1 + φ L p, x γ φ L C φ L 1 + φ L q. < 1, q < n γ +1, Proof. We have δp = γ when p = n/γ, and δq = γ+1 when q = n/γ +1. Therefore p < n/γ and q < n/γ +1 if δp = γ +δ/ and δq = γ+1 + δ 4. Let χ Cc R +,[0,1] with χ 1 on [0,1]. We have x γ φ L χ x γ φ L + 1 χ x γ φ L χ x γ L p φ L p + 1 χ x γ L φ L 1 C φ L p + φ L 1, where we have used x x γ L p loc Rn because p < n/γ. The other estimate is similar, since x γ = Ø x γ 1. Lemma 4.5. Let γ < 1 and f Σ H R n. Recall that g is defined by 4.1. We have: x γ f W, ; g L loc R;W, ; fe ig, x γ f fe ig L loc R;H. Proof. From Lemma 4.4 and Sobolev embeddings, x γ f L f + f p f H 1, x γ f L f + f q f H 1, x γ f L x γ f L f f + f q f q f H. Since t cost γ L 1 loc R, we infer that g L loc R;W,. The last two properties follow easily.

12 1 R. CARLES, N. MAUSER, AND H. P. STIMMING We can now replace U ε with V ε in 4.8, up to the following error. From Lemmas 4., 4.4 and 4.5, U ε t V ε t x γ f a ε t εtant θ x γ f a ε t H 1 εtant θ x γ f W 1, a ε t H 1 εtant θ a ε t + x a ε t εtant θ f + x b ε e ig εtant θ 1+ x b ε t εtant θ 1+ x b ε t f t εtant 1+ θ t b ε τ H 1 dτ, for 0 < θ 1/ to be fixed later. Plugging this estimate into 4.8 yields t b ε t x γ u ε u ε 1 t cost γvε t x γ f a ε t εtant θ cost γ 1+y εt. We check that cost γvε t x γ f φ = x γ V ε tf V ε tφ. Since we expect V ε ta ε t to be close to U ε ta ε t = u ε t as ε 0, we estimate the difference x γ V ε tf V ε ta ε t U ε ta ε t x γ V ε tf L V ε t U ε tb ε e ig V ε tf + Vε tf L εtant θ b ε te igt H p 1 1+ cost δp εtant θ b ε t f H 1 + f H 1, 0 using the modified Sobolev inequality.8. Since δp = n p 4.9 that > γ, we infer from 4.11 t b ε t εtant θ cost δp 1+y εt + x γ u ε t V ε tf u ε t. From Lemma 4.4, the last term is estimated, up to a constant, by u ε t V ε tf L 1 + u ε t V ε tf L p 4.1 u ε t V ε t fe igt L u ε t + V ε t fe igt L + u ε t V ε t fe igt u ε V t Lp p + ε t fe igt. Lp For the first term of the right hand side, we have, since U ε is unitary on, U ε t b ε e ig V ε t fe ig b ε t f + U ε t V ε t fe igt L y ε t+ εtant θ fe igt. H 1

13 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 13 In addition, notice that u ε t = V ε tf = f. The second term is estimated thanks to the modified Gagliardo Nirenberg inequality.8, u ε t V ε t fe igt L p cost δp u ε t V ε t fe igt 1 δp L J ε t u ε t V ε t fe igt δp. The first norm was estimated just above. For the second one, notice that J ε tu ε t = iu ε t x ; J ε tv ε t = iv ε t x, therefore: J ε t u ε t V ε t fe igt L U ε t b ε te igt V ε t fe igt L b ε te igt fe igt L + U ε t V ε t fe igt L y ε t+ εtant θ, where we have used Lemmas 4. and 4.5. We infer that u ε t V ε t fe igt L y p cost δp ε t+ εtant θ. We have explicitly V ε t fe igt L = cost δp f p p cost δp. Proceeding as above, we have u ε t p cost δp u ε 1 δp J ε tu ε δp, with J ε tu ε b ε t f H 1 + f H 1. These estimates will eventually lead to an inequality of the form y ε t aty εt+bty ε t κ +ct, for some κ > 1. To avoid that situation, we proceed as in Section 3; there exists t ε > 0 such that 4.13 b ε t H 1 f H 1, for t [0,t ε ]. So long as 4.13 holds, we have from the above estimates 4.14 t b ε t cost δp y ε t+ εtant θ. To prove that 4.13 holds up to time π/ Λε for 0 < ε ε Λ along with the error estimate 4.7, we estimate the norm of x t b ε. From 4.5 and 4.6, x t b ε t = i x gt t b ε t ie igt x U ε t x γ u ε u ε t 1 x γ cost γ f a ε t. The first term is controlled thanks to Lemma 4.5 and For the other term, we notice that since U ε is unitary on, from.5 its norm is equal to: Jε t x γ u ε u ε i t+ cost γuε t x x γ f a ε t. We proceed as before: we first replace U ε with V ε in the last term, up to an error of cost γ times: U ε t V ε t x x γ f a ε L U ε t V ε t x x γ f b ε fe ig + U ε t V ε t x x γ f fe ig.

14 14 R. CARLES, N. MAUSER, AND H. P. STIMMING For the first term, we do not use Lemma 4., but roughly the fact that U ε and V ε are unitary on. It is not larger than x x γ f b ε fe ig b ε t f H 1, from Lemma 4.5. The second term is controlled thanks to Lemmas 4. and 4.5, U ε t V ε t x x γ f fe ig εtant θ. We now have, so long as 4.13 holds, 4.15 t b ε t H 1 cost δp y ε t+ εtant θ + J ε t x γ u ε u ε i + cost γvε t x x γ f a ε t L. Using the identity J ε tv ε t = iv ε t x and 4.10, we have to estimate x t Jε γ u ε u ε 1 cost γvε t x γ f a ε t 4.16 = J ε t x γ u ε u ε x γ V ε tf V ε ta ε x γ u ε J ε tu ε x γ V ε tf J ε tv ε ta ε + cost x x γ u ε u ε x x γ V ε tf V ε ta ε. We replace V ε with U ε in the first term of the right hand side, up to the error x γ V ε tf J ε tv ε t U ε ta ε V ε tf + Vε tf p V ε t U ε t x a ε cost δp V ε t U ε t x b ε fe ig + cost δp V ε t U ε t x fe ig cost δp b ε t f H 1 + εtant θ, from the above computation. Therefore, the first term of the right hand side of 4.16 is estimated by cost δp y ε t+ εtant θ + x γ u ε V ε tf J ε tu ε. So long as 4.13 holds, J ε tu ε 1, and the last term is estimated by x γ u ε V ε tf L, which already appeared above and was estimated in 4.1. We are left with the second term of the right hand side of Using Lemma 4.4 with q instead of p now, x x γ V ε tf V ε t U ε ta ε V ε tf + Vε tf q V ε t U ε ta ε cost δq y ε t+ εtant θ. The final term to estimate is x x γ u ε V ε tf u ε u ε V ε tf L 1 + u ε V ε tf L q. The right hand side was already estimated in 4.1 with p instead of q. We finally have, so long as 4.13 holds, y t cost δp + cost 1 δp yε t+ εtant θ.

15 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 15 Now recall that given δ > 0, δp and δq are explicit, hence y t cost γ δ yε t+ εtant θ. It is now time to fix θ. In view of 4.7, it is natural to take θ = 1 γ δ. This yields, so long as 4.13 holds, 4.17 y ε t cost γ δ yε t+ εtant 1 γ δ cost γ δ yε t+ ε1 γ δ. cost 1 δ The maps t cost γ δ and t cost 1+δ are locally integrable we can assume γ + δ/ < 1 δ/, otherwise 4.7 is of no interest. From the Gronwall lemma, so long as 4.13 holds, we infer 4.18 y ε t ε 1 γ δ. Therefore, there exists ε Λ > 0 such that for 0 < ε ε Λ, 4.13 holds up to time π/ Λε, with The estimate for x t b ε then is easy, we leave out this part; this proves 4.7. Remark. One might believe that we could deduce Proposition 4.3 in one shot from 4.17, and wonder why we split the proof into three steps. The reason is that we cannot apply Lemma 4. which was used to get 4.17 near t = π/. On the other hand, we will see below that computations near t = π/ are far simpler Near the focus and beyond. Keep Λ, δ > 0 fixed. We prove that there exists C Λ,δ such that 4.19 π +Λε t b ε t Σ dt C Λ,δ ε 1 γ δ. π Λε A rough estimate in 4.8 yields 4.0 t b ε t x γ u ε u ε t + u ε t + uε t p u ε t + 1 cost γ x γ f a ε t 1 cost γ uε t. The conservation of mass yields u ε t = f. The conservations of mass and energy.3 yield, along with Gagliardo Nirenberg inequalities, u ε t p ε δp. Using this estimate which is sharp near the focus, and only near the focus and integrating 4.0, we get π +Λε t b ε t dt Λε 1 δp π +Λε dt + π Λε π Λε cost γ ε 1 γ δ +ε 1 γ. The term x t b ε t is estimated the same way, since the conservation of energy yields an a priori bound for H ε u ε. For x t b ε t, we proceed as in Section 3, 3.10 to get an estimate from Gronwall lemma; the details are left to the reader. Finally, one can prove that there exists C Λ,δ such that π π +Λε t b ε t Σ dt C Λ,δ ε 1 γ δ by mimicking the computations performed in Section 4., and the proof of Proposition 4.3 is complete.

16 16 R. CARLES, N. MAUSER, AND H. P. STIMMING 5. Linear propagation and nonlinear focus We now consider the case where α = γ > 1 in 1.6. Our results are similar to those of [6]. Before stating the main result, we recall some points of the scattering theory for 1.9. Proposition 5.1 [13, 17]. Assume ψ Σ and 1 < γ < min4,n. If γ > 4/3 or if ψ Σ is sufficiently small, then: There exists a unique ψ CR t,σ solution of 1.9, such that lim t ψ U tψt Σ = 0, where Ut = e i t. There exists a unique ψ + Σ such that The scattering operator is S : ψ ψ +. Our main result in this section is: lim ψ + U tψt Σ = 0. t + Proposition 5.. Suppose n. Let f Σ, 1 < γ = α < min4,n, and k N. Assume either γ > 4/3 or f Σ is sufficiently small. Then the asymptotics of u ε for π/+k 1π < a b < π/+kπ is given by sup Aε t u ε t,x e ink π F S k F 1 x L f e ix a t b cost n/ ε tant 0, cost ε 0 x where A ε is either of the operators Id, J ε or H ε, and S k denotes the k-th iterate of S which is well defined under our assumptions on f. At the focuses: Bε u ε π +kπ e ink π F S k ε f 0, ε n/ ε 0 where B ε is either of the operators Id, x ε or ε x. With Lemma 4. in mind, this shows that nonlinear effects are negligible away from focuses, while they have an influence at leading order near the focuses: each caustic crossing is described in average by the nonlinear scattering operator S the phase shift e ink π is the Maslov index, present in the linear case [9]. The proof of Proposition 5. is very similar to the one in [6], which relies on scaled Strichartz estimates. We will refrain from repeating everything in detail and limit ourselves to prove the main technical proposition and present an outline for the rest of the proof. One main difference to the problem in [6] is the action of the operators J ε t, H ε t on the Hartree nonlinearity as described by.9. We start by reformulating Equation 1.6 by the Duhamel formula 5.1 u ε t = U ε t t 0 u ε 0 iεγ 1 t t 0 U ε t sf ε u ε sds iε 1 t t 0 U ε t sh ε sds. This equation generalizes Eq. 1.6 to the case of an additional source term and a general nonlinear term F ε. The main technical result which is used throughout the proof of Proposition 5. is:

17 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 17 Proposition 5.3. Let t 1 > t 0, with t 1 t 0 π. Let q,r,s,k [1, ] be such that: 1 a r = 1 r + s + γ n 1 and s < n n γ, 1 b 5. q = 1 q + k, c q,r is an admissible pair, d 0 < 1 < δs < 1. k Assume that there exists a constant C independent of t and ε such that for t 0 t t 1, 5.3 F ε u ε t L r C t x cost +ε δs uε L r x, and define /k t1 A ε dt t 0,t 1 :=. t 0 cost +ε kδs Then there exists C independent of ε, t 0 and t 1 such that for any admissible pair ρ,σ, 5.4 u ε Lq t 0,t 1;L r C ε 1/q u ε 0 +C q,ρε 1 1 q 1 ρ h ε L ρ t 0,t 1;L σ Mostly the following corollary is applied: +C ε δs 1 k A ε t 0,t 1 u ε Lq t 0,t 1;L r, Corollary 5.4. Suppose the assumptions of Prop. 5.3 are satisfied. Assume moreover that C ε δs k 1 A ε t 0,t 1 1/, which holds in either of the two cases, 0 t 0 t 1 π Λε, with Λ Λ 0 sufficiently large, t 0,t 1 [ π Λε, π t1 t0 +Λε], with ε η sufficiently small. Then 5.5 u ε L t 0,t 1; C u ε 0 +C q,ρ ε 1 1 ρ h ε L ρ t 0,t 1;L σ. To prove Proposition 5.3, we first prove the following algebraic lemma: Lemma 5.5. Let n, and assume 1 < γ < min4,n. Then there exist q,r,s,k [1, ] satisfying the conditions 5.. Proof. Note that a is equivalent to demanding γ/ = δr+δs and γ/ > δs. Case γ : Suppose γ/ = δs. Then, by the first half of a δr = 0 and q,r =, by c. With k =, b and d are satisfied. Now choose s such that 1/ < δs < γ/, but close enough to γ/ for 5. still to be valid by continuity for example δs = γ 1. Then 5. is satisfied. Case γ > : In this case take s such that δs = 1, e.g. s = n n. Up to a continuity argument as in the previous case, δs < 1 and 5. is satisfied. Proof of Proposition 5.3. Application of thescaled Strichartz estimatesprop..1 to equation 5.1 yields u ε L q t 0,t 1;L r Cε 1/q u ε 0 +C q,ρ ε 1 1 q 1 ρ h ε L ρ t 0,t 1;L σ +Cε γ 1 q F ε u ε L q t 0,t 1;L r. Then by the assumptions on F ε u ε, after an application of Hölder inequality in time, the statement follows.

18 18 R. CARLES, N. MAUSER, AND H. P. STIMMING Proof of Corollary 5.4. The additional assumption implies that the last term in 5.4 can be absorbed by the left hand side, and we get 5.6 u ε Lq t 0,t 1;L r Cε 1/q u ε 0 +Cε 1 1 q 1 ρ h ε L ρ t 0,t 1;L σ. Another application of Strichartz estimates to equation 5.1, with indices, on the left and ρ,σ respectively q,r on the right, yields u ε L t 0,t 1; C u ε 0 +Cε 1 1 ρ h ε L ρ t 0,t 1;L σ As before, +Cε γ 1 1 q F ε u ε L q t 0,t 1;L r. ε γ 1 1 q F ε u ε L q t 0,t 1;L r Cε1 q ε δs 1 k A ε t 0,t 1 u ε L q t 0,t 1;L r and the statement now follows from 5.6. Cε 1 q u ε L q t 0,t 1;L r, The proof of Proposition 5. consists of three parts: the propagation before the focus, the matching between the two regimes, and proof that near the focus, the harmonic potential is negligible. In all parts the main tool to derive the major statements is Prop Since the proof is very similar to the one in [6], we do not repeat everything in detail but give a detailed proof only for the first part to show how the methods of [6] are applied. We now show the proof for the propagation before the focus, that is the approximation of u ε t by u ε free t for 0 t π Λε, in the limit Λ +. We prove: lim sup ε 0 A sup ε tu ε t,x u ε 0 t π Λε free t,x L x 0, Λ + with A ε t being either of the operators Id, J ε or H ε. Define the remainder w ε = u ε u ε free. It solves iε t w ε + 1 ε w ε = Vxw ε +ε γ x γ u ε u ε, w t=0 ε = 0. From Duhamel s principle, this can be written as t 5.7 w ε t = U ε tr ε iε γ 1 U ε t s x γ u ε u ε sds. 0 Since u ε free solves the linear equation 1.7, so does Jε tu ε free from.6, and u ε free t = f ; Jε tu ε free = f. From the Sobolev inequality.8, u ε free t L s C cost δs f 1 δs L f δs for any s [, n n [. Therefore there exists C 0 such that 5.8 u ε freet L s C 0 cost δs From Prop.., for fixed ε > 0, u ε CR,Σ, and the same obviously holds for u ε free. Therefore, wε CR,Σ, and there exists t ε > 0 such that 5.9 w ε t L s C 0 cost δs,

19 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIAL 19 for any t [0,t ε ]. So long as 5.9 holds, we have u ε t L s C 0 cost δs, and we can apply Prop Takeh ε = ε γ x γ u ε u ε free andfε w ε = x γ u ε w ε andletq,k,r,s [1, ] satisfy the assumptions of Prop Now by Hölder s inequality, F ε w ε t L r x γ u ε t L β w ε t L r with β such that 1 r = 1 r + 1 β. By the Hardy Littlewood Sobolev inequality and the above estimate, F ε w ε t L r uε t L s wε t L r C 0 cost δs wε t L r. Note that the second statement of 5.a ensures that s,β 1, so the Hardy- Littlewood-Sobolevinequalityisapplicablehere. Assume5.9holdsfor0 t T ε. If 0 t T ε π Λε, then ε cost, and the above estimate shows that Fε satisfies assumption 5.3. From Corollary 5.4, if Λ is sufficiently large, we get for 0 t T ε π Λε: w ε L 0,T; C σ ε γ 1 1 ρ x γ u ε u ε free L ρ 0,T;L σ for any admissible ρ,σ. Now take ρ,σ = q,r and proceed as above in space, and apply Hölder inequality in time: x γ u ε u ε free L q 0,T;L r C γ,n u ε L k 0,T;L s uε free L q 0,T;L r. The first term of the right-hand side is estimated through 5.8 and 5.9: u ε L k 0,T;L s C π T δs 1/k, the last term is estimated the same way, for 5.8 still holds when replacing s with r: u ε free L q 0,T;L r C. π T δr 1/q We infer: x γ u ε u ε free L q 0,T;L r C, π T γ 1 1 q thus ε 5.10 w ε L 0,T; C π T γ 1 1 q. Now apply the operator J ε to 5.7. Since J ε and U ε commute, it yields, t J ε tw ε = U ε tj ε 0r ε iε γ 1 U ε t sj ε s x γ u ε u ε sds. 0 The action of J ε on the nonlinear term is described by.9. In order to apply Prop. 5.3 like before, we take now h ε = ε γ 1 x γ u ε J ε tu ε free +εγ 1 Re x γ u ε J ε tu ε free u ε and 5.11 F ε w ε = x γ u ε J ε tw ε + Re x γ u ε J ε tw ε u ε.

20 0 R. CARLES, N. MAUSER, AND H. P. STIMMING The first term on the r.h.s of 5.11 leads to an equation which is very similar to 5.7, with w ε replaced by J ε w ε and is treated by the same computations as above. For the second term, we estimate by Hölder, by the Hardy-Littlewood-Sobolev inequality and then again by Hölder: Re x γ u ε J ε tw ε u ε L r x γ u ε J ε tw ε L β 1 u ε t L s u ε t L s J ε tw ε t L r u ε t L s with r,s as stated in 5. and 1 r = 1 β s. Here the condition to use the Hardy- Littlewood-Sobolev inequality is γ > δr + δs, which is always satisfied by 5.. By applying now 5.9 we continue to estimate C 0 cost δs Jε tw ε t L r. Then we apply, like before, Prop. 5.3 and estimate the term for h ε as above, to obtain: γ 1 1 ε 5.1 J ε w ε L 0,T; C r ε q +C π. T Combining 5.10 and 5.1 yields, along with.8, t [0,T], w ε C ε t L s r ε cost δs H 1 + π t γ 1 1 q. Therefore, choosing ε sufficiently small and Λ sufficiently large, we deduce that we can take T = π Λε. With the result of Lemma 4. on the limit of uε free, this yields Prop. 5., away from the focus, for A ε = Id and J ε. The case A ε = H ε on this time interval is now straightforward. The remaining parts of the prooffor Prop. 5. are done as in [6] with the method changes as in the part shown above. It remains to show that the approximations in the two different regimes match at t = π Λε, and that the influence of the harmonic potential is small near the focus so that the propagation there is given by 5.13 v ε t,x = 1 t π ε ψ, x, n/ ε ε where ψ is the solution of 1.9 subject to the following initial condition at t = lim sup ε 0 π Λε t π +Λε U tψt t= = e inπ 4 f. This solution exists according to Proposition 5.1. Then the following asymptotic is proven: sup A ε tu ε t,x v ε t,x x 0, Λ + with A ε t being one of the operators Id, J ε or H ε. Since these parts are quite similar to the treatment in [6], we do not repeat them. After the crossing of the first focus, the solution is again propagated linearly and at subsequent focusing points this process is iterated.

21 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIA1 6. Formal Computations and Discussions 6.1. The case α = γ = 1 in 3-d: Schrödinger-Poisson. We saw in Section 5 that when α = γ > 1, the nonlinear term in 1.6 has a leading order influence near the focuses, and only in these regions. On the other hand, if α = 1 and γ < 1, Section 4 shows that the Hartree term cannot be neglected away from the focuses. These two cases suggest that when α = γ = 1, the nonlinearinfluence is everywhere relevant. The aim of this final section is to give convincing arguments that this is the case. For the influence near the focuses, we need the scattering theory for 1.9 at γ = 1. In this long range scattering case, modified scattering operators are needed instead of the ones described in Prop Hayashi and Naumkin [16] obtained an asymptotic completeness result for n with smoothness assumptions which are applicable to our situation. On the other hand, they could not obtain wave operators. Ginibre and Velo [14, 15] obtained modified wave operators for 1.9 with γ = 1 using Gevrey spaces by a technically involved method. A drawback of both these results is that they include a loss in regularity. To show how the long range scattering theory fits into our framework we report a particular case of the result of Hayashi and Naumkin [16]. Proposition 6.1 [16]. Assume n = 3, ϕ Σ, and δ = ϕ Σ is sufficiently small. Let ψ CR,Σ be the solution of 1.9 with ψ t=0 = ϕ. Then there exists a unique function ψ + H σ,0 H 0,σ, 1 < σ < 1, such that ψt exp i x 1 ψ + x L log t Utψ + t 0, t + where H α,β = {φ S 1+ x β/ 1+ α/ φ < }. To summarize very roughly, the results in [14, 15] consist in showing that given some ψ + or ψ for an asymptotic behavior for t, one can find ψ solving 1.9 such that the above asymptotics holds. Analogoulsy to the treatment of long-range scattering in [4], one can now define g ε t,x := x 1 f xlog cost ε comparewith 4.1 andadd the phasegε t=0 to the initial data in 1.6. This yields: u ε t=0 = fxe i x 1 f xlogε. Using the modified scattering operators from the results of [14, 15] we get, at least formally, for 0 t < π/, u ε 1 x t,x e cost 3/f ix ε tant+igε x t, cost as ε 0. cost This asymptotic also stems from the same computations as those performed in Section 4.1. Notice that the matching for t π = Øε is similar to the one in [6], except that we now have to take the presence of g ε into account. This is where changing the integration from 0 to t in 4.1 into the above definition of g ε makes the matching possible. Indeed, for t π = Øε, we compare uε with the function v ε given by 5.13, where ψ is now the solution given by the long range wave operators constructed in [14, 15]. To make this statement more precise and the link between 4.1 and the definition of g ε more explicit, notice that we have, as t π : g ε x t, x 1 f x π π cost t log t phase shift for v ε ε x 1 f x t dτ compare with 4.1. cost cosτ arccos ε

22 R. CARLES, N. MAUSER, AND H. P. STIMMING The effects of the nonlinearity show up in g ε. Using the scaling 5.13 we can then formally continue with Prop. 6.1: for π/ < t < 3π/, u ε t,x e i3π F cost S F 1 x f e ix 3/ ε tant+ihε x t, cost as ε 0, cost where S is the map S : ψ ψ +, where ψ is the asymptotic state of the result of [14], which yields some solution ψ to 1.9, and ψ + is provided by Prop h ε is given by h ε t,x := x 1 F S F 1 f cost xlog. ε TheactionofF S F 1 onf accountsfornonlineareffectstakingplaceatthefocus, and h ε for nonlinear effects after the focus. So the influence of the nonlinearity will be relevant at all times. Theimpossibilityto define ascatteringoperatorforthiscaseisoneofthe reasons why this argument is only formal. Remark. A rigorous result could be obtained with the same approach as in [3]. It would consist in studying the system of linear equations with a nonlinear coupling, iε t u ε + 1 ε u ε = x uε, iε t u ε + 1 ε u ε = x uε +ε x 1 u ε u ε. The first equation is solved explicitly thanks to Mehler s formula, and the second one is a linear Schrödinger equation with a harmonic potential and a time-dependent perturbation. With the oscillatory integral used in Section 4, and adapting the results of [8], one could prove similar asymptotics to those stated above. 6.. The case of an additional local strong nonlinearity. We now consider equation 1.6 with an additional nonlinear term that is a multiplication operator with a power of the density u ε. Such equations arise in the modeling of effective one particle Schrödinger equations where exchange terms like in the Hartree-Fock equation are simplified to functionals of the local densities, i.e. time dependent density functional theory, with the Schrödinger- Poisson-Xα equation as the simplest of such models see [5] and [1] for a heuristic derivation and numerical simulations. Note that the additional local term has the opposite sign than the Hartree term corresponding to the physical fact that the exchange-correlation hole weakens the direct Coulomb interaction. We will hence consider the following class of semi-classical Hartree equations 6.1 iε t u ε + 1 ε u ε = x uε +ε α x γ u ε u ε ε β u ε σ u ε, with α 1, β 1, γ > 0 for x R n, and with a σ that is sub-critical with respect to finite time blow-up, i.e. n > σ > 0. We can now discern the influence of the two nonlinear terms in the classical limit in terms of: The size of the scaling exponents α, β with respect to the critical value. The relation between the scaling and the strength of the nonlinearities determined respectively by γ and σ. If we take α > 1 and β > 1, by [6] and Section 5 we find that the classical limit is given by the linear propagation as long as no focusing occurs. At the focus, the relevant discrimination is σ = β/n or < β/n for the power nonlinearity and γ = α

23 SEMICLASSICAL HARTREE EQUATION WITH HARMONIC POTENTIA3 or < α for the Hartree term. If σ = β/n and γ < α, the crossing of the focus will be described by the scattering operator for NLS when it is defined, if on the other hand σ < β/n and γ = α and the assumptions of Prop. 5.1 are satisfied, focus crossing will be determined by the scattering operator of Prop If both nonlinearities are at the critical strength σ = β/n and γ = α, then both will have an influence in crossingthe caustic. If, on the other hand, both σ < β/n and γ < α, the nonlinear influence will be negligible everywhere. If at least one of the scaling exponents α and β is equal to 1 and, at the same time, both σ < β/n and γ < α, the corresponding nonlinear term will be relevant in the WKB propagation before the focusing. At the focus, the nonlinear terms will not be relevant and the crossing of the focus will be as in Prop If σ = β/n and γ = α then there will be a nonlinear influence everywhere and long range scattering for NLS and/or Hartree has to be taken into account. The influence of the nonlinear action for the single power NLS and the Hartree equation is summed up in two tables, for Hartree the table is given in the introduction, for single power nonlinear Schrödinger equation, it is stated in [3]. The behavior of 6.1 can be described by independently superposing these two tables. The following table is an extract from that superposition: α > γ and β > σn α = γ or β = σn α > 1 and β > 1 Linear WKB, Linear WKB, linear focus nonlinear focus α = 1 or β = 1 Nonlinear WKB, Nonlinear WKB, linear focus nonlinear focus Nonlinear WKB respectively nonlinear focus here stands for an influence from at least one of the nonlinear terms away from the focus or close to the focus Wigner measures. We already mentioned in the introduction the work of Zhang, Zheng and Mauser [31] where the semiclassical limit of the Schrödinger- Poisson equation with no smallness assumption on the initial data or the nonlinearity is studied by means of Wigner measures. Wigner measures have proven to be efficient tools for linear semi-classical problems and for homogenization limits; see [6] for an overview on Wigner measure limits of Hartree equations. Wigner measureshavethe merit that in phase space the caustics ofphysical spaceare somewhat unfolded and that generally, results globally in time are possible. In [5], the Wigner measure of the nonlinear Schrödinger equation with power-like nonlinearity studied in [3] is investigated. It is shown that the Wigner measure leads to an ill-posed problem whenever nonlinear effects at the focal points come into play. In other words, the Wigner measure can only be valid as long as no caustic appears. We briefly discuss the Wigner measures of 1.6 in view of these results. The Wigner measure of the family u ε t 0<ε 1, which is bounded in, is the weak limit under ε 0 up to an extraction of its Wigner transform, W ε u ε t,x,ξ = 1 π n u ε t,x vε u ε t,x+ vε e iξ v dv. This limit is a positive radon measure µ and is in general not a unique limit. linear case: Case α > γ, α > 1: By the result of Prop. 3.1 and the asymptotics of u free in Lemma 4., the Wigner