Madelung, Gross-Pitaevskii Korteweg Rémi Carles

Raphaël Danchin
email: danchin@univ-paris12.fr

Jean-Claude Saut
email: jean-claude.saut@math.u-psud.fr

Rémi Carles
email: remi.carles@math.cnrs.fr

Émi Madelung

RAPHA ËL DANCHIN Jean-Claude Carles

Saut

la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

In his seminal work [START_REF] Madelung | Quanten theorie in Hydrodynamischer Form[END_REF] (see also [START_REF] Landau | Tome III: Mécanique quantique. Théorie non relativiste[END_REF]), E. Madelung introduced the so-called Madelung transform in order to relate the linear Schrödinger equation to a hydrodynamic type system. This system takes (slightly) different forms according to the context: linear or nonlinear equation with various nonlinearities. For the semiclassical nonlinear Schrödinger equation (shortened in NLS in what follows):

iε∂ t ψ ε + ε 2 2 ∆ψ ε = f |ψ ε | 2 ψ ε ,
the Madelung transform amounts to setting ψ ε (t, x) = ρ(t, x)e iφ(t,x)/ε , so as to get the following system for ρ and v = ∇φ:

(1.1)

     ∂ t v + v • ∇v + ∇f (ρ) = ε 2 2 ∇ ∆ √ ρ √ ρ , ∂ t ρ + div (ρv) = 0.
This system is referred to as the hydrodynamic form of NLS because of its similarity with the compressible Euler equation (which corresponds to ε=0). The additional term on the right-hand side is the so-called quantum pressure.

Madelung transform is crucial to investigate qualitative properties of the nonlinear Schrödinger equation with nonzero boundary conditions at infinity whenever the solution is not expected to vanish "too often". Of particular interest is the so-called Gross-Pitaevskii equation

iε∂ t ψ ε + ε 2 2 ∆ψ ε = (|ψ ε | 2 -1)ψ ε ,
which corresponds to f (r) = r -1 or, more generally, the case where f is a smooth function vanishing at some r 0 and such that f ′ (r 0 ) < 0. This covers in particular the "cubic-quintic" NLS (f (r) = -α 1 + α 3 r -α 5 r 2 with α 1 , α 3 , α 5 > 0). Using Madelung transform is also rather popular to study the semi-classical limit of NLS. This requires ρ to be nonvanishing, though. We shall see below to what extent the presence of vacuum (that is of points where ρ vanishes) is merely a technical problem due to the approach related to Madelung transform.

The hydrodynamic form of NLS may be seen as a particular case of the system of quantum fluids (QHD) with a suitable choice of the pressure law, which, in turn, enters in the class of Korteweg (or capillary) fluids. We aim here at further investigating the link between these three a priori disjoint domains, through various uses of the Madelung transform. In passing, we will also review some more or less known facts pertaining to the theories of quantum fluids and Korteweg fluids. The most original part of this work is a simple proof of a recent result of Antonelli and Marcati [START_REF] Antonelli | The quantum hydrodynamics system in two dimensions[END_REF][START_REF]On the finite energy weak solutions to a system in quantum fluid dynamics[END_REF], on the global existence of weak solutions to a quantum fluids system.

1.1. Organization of the paper. The paper is organized as follows. The second section is a review of the connexions of the Madelung transform with the semiclassical limit of the nonlinear Schrödinger equation. In passing we briefly present the state-of-the-art for the Cauchy problem for the Gross-Pitaevskii equation. In Section 3, we solve the quantum hydrodynamical system (1.1) by a direct method based on the use of an extended formulation, and explain how it may be adapted to tackle general Korteweg fluids. In Section 4, we review the use of the Madelung transform and of the hydrodynamical form of the Gross-Pitaevskii equation to study the existence and properties of its traveling wave solutions and of its transonic limit, both in the steady and unsteady cases. We give in Section 5 a simple proof of the aforementioned result of Antonelli and Marcati ( [START_REF] Antonelli | The quantum hydrodynamics system in two dimensions[END_REF][START_REF]On the finite energy weak solutions to a system in quantum fluid dynamics[END_REF]). Lastly we list in Appendix the basic conservation laws for the Schrödinger, the QHD and the compressible Euler equations, and show that some of these laws naturally carry over to general Korteweg fluids.

Notations.

• We denote by |•| p (1 p ∞) the standard norm of the Lebesgue spaces L p (R d ).

• The standard H s (R d ) Sobolev norm will be denoted • s .

• We use the Fourier multiplier notation: f (D)u is defined as F (f (D)u)(ξ) = f (ξ) u(ξ), where F and • stand for the Fourier transform. • The operator Λ = (1 -∆) 1/2 is equivalently defined using the Fourier multiplier notation to be Λ = (1 + |D| 2 ) 1/2 . • The partial derivatives will be denoted with a subscript, e.g. u t , u x , . . . or ∂ t u, ∂ x u, . . . or even ∂ j (to designate ∂ xj ). • C will denote various nonnegative absolute constants, the meaning of which will be clear from the context.

Madelung transform and the semi-classical limit of NLS

As already mentioned in the Introduction, we consider the equation (2.1)

iε∂ t ψ ε + ε 2 2 ∆ψ ε = f |ψ ε | 2 ψ ε ; ψ ε (0, x) = ρ 0 (x)e iφ0(x)/ε .
Here we have in mind the limit ε → 0. The space variable x belongs to R d in this paragraph. The periodic case x ∈ T d being of particular interest for numerical simulations however (see e.g. [START_REF] Bao | Numerical study of time-splitting spectral discretizations of nonlinear Schrödinger equations in the semiclassical regimes[END_REF][START_REF] Degond | An asymptotic preserving scheme for the Schrödinger equation in the semiclassical limit[END_REF] and references therein), we will state some analogous results in that case, which turn out to be a little simpler than in the case x ∈ R d . In this section, we shall focus on two types of nonlinearity:

• Cubic, defocusing nonlinearity1 :

f |ψ ε | 2 ψ ε = |ψ ε | 2 ψ ε .
The Hamiltonian associated to (2.1) then reads

H ε NLS (ψ ε ) = ε∇ψ ε 2 L 2 + ψ ε 4 L 4 . It is well defined on the Sobolev space H 1 (R d ) in dimension d 4.
• Gross-Pitaevskii equation:

f |ψ ε | 2 ψ ε = |ψ ε | 2 -1 ψ ε .
In this case, the natural Hamiltonian associated to (2.1) is

H ε GP (ψ ε ) = ε∇ψ ε 2 L 2 + |ψ ε | 2 -1 2 L 2 .
In both cases, the Hamiltonian defines an energy space, in which existence and uniqueness for the Cauchy problem (2.1) have been established.

The case of a defocusing cubic nonlinearity is now well understood. In dimension d 4 the corresponding NLS equation is globally well-posed in H 1 (R d ), and the additional H s (R d ) regularity (s 1) is propagated (see the textbooks [START_REF] Cazenave | Semilinear Schrödinger equations[END_REF][START_REF] Linares | Introduction to nonlinear dispersive equations[END_REF][START_REF] Tao | Nonlinear dispersive equations[END_REF] and the references therein).

The situation is more complicated for the Gross-Pitaevskii equation, where the finite energy solutions ψ ε cannot be expected to be in

L 2 (R d ), since |ψ ε | 2 -1 ∈ L 2 (R d ).
As noticed in [START_REF] Zhidkov | The Cauchy problem for a nonlinear Schrödinger equation[END_REF][START_REF] Vries | nonlinear Schrödinger equations: qualitative theory[END_REF], and extended in [START_REF] Gallo | Schrödinger group on Zhidkov spaces[END_REF], a convenient space to study the Gross-Pitaevskii equation is the Zhidkov space:

X s (R d ) = {ψ ∈ L ∞ (R d ) ; ∇ψ ∈ H s-1 (R d )} with s > d/2.
Remark 2.1. In the case x ∈ T d , the spaces H s (T d ) and X s (T d ) (with obvious definitions) are the same.

In the case x ∈ R d , as a consequence of the Hardy-Littlewood-Sobolev inequality (see e.g. [START_REF] Hörmander | The analysis of linear partial differential operators[END_REF]Th. 4.5.9] or [28, Lemma 7]), one may show that if d 2 and ψ ∈ D ′ (R d ) is such that ∇ψ ∈ L p (R d ) for some p ∈]1, d[, then there exists a constant γ such that ψ -γ ∈ L q (R d ), with 1/p = 1/q + 1/d. Morally, γ is the limit of ψ at infinity. If d 3 then we can take p = 2, so every function in X s (R d ) satisfies the above property; for d = 2, the above assumption requires a little more decay on ∇ψ than general functions in X s (R2 ).

The well-posedness issue in the natural energy space

E(R d ) = {ψ ∈ H 1 loc (R d ), ∇ψ ∈ L 2 (R d ), |ψ| 2 -1 ∈ L 2 (R d
)} associated to the Gross-Pitaevskii equation has been investigated only very recently by C. Gallo in [START_REF]The Cauchy problem for defocusing nonlinear Schrödinger equations with nonvanishing initial data at infinity[END_REF] and P. Gérard in [START_REF]The Cauchy problem for the Gross-Pitaevskii equation[END_REF][START_REF]The Gross-Pitaevskii equation in the energy space, Stationary and time dependent Gross-Pitaevskii equations[END_REF], in dimension d 4.

Let us emphasize that in the R d case, the energy space is no longer a linear space (contrary to the case of zero boundary condition a infinity where it is

H 1 (R n )), hence solving the Gross-Pitaevskii equation in E(R d ) is more complicated. How- ever, if d = 3, 4, one may show that E(R d ) coincides with {ψ = c(1 + v), c ∈ S 1 , ψ ∈ H 1 (R d ), 2 Re ψ + |ψ| 2 ∈ L 2 (R d )},
which allows to endow it with a structure of a metric space. The case d = 2 is slightly more technical.

2.1. Some issues related to the use of the Madelung transform. In the semi-classical context, Madelung transform consists in seeking

(2.2) ψ ε (t, x) = ρ(t, x)e iφ(t,x)/ε
for some ρ 0 and real-valued function φ. Plugging (2.2) into (2.1) and separating real and imaginary values yields:

(2.3)

       √ ρ ∂ t φ + 1 2 |∇φ| 2 + f (ρ) = ε 2 2 ∆ √ ρ ; φ |t=0 = φ 0 , ∂ t √ ρ + ∇φ • ∇ √ ρ + 1 2 √ ρ∆φ = 0 ; ρ |t=0 = ρ 0 .
Two comments are in order at this stage: the first equation shows that φ depends on ε and the second equation shows that so does ρ in general. We shall underscore this fact by using the notation (φ ε , ρ ε ). Second, the equation for φ ε can be simplified, provided that ρ ε has no zero. Introducing the velocity v ε = ∇φ ε , (2.3) yields the system of quantum hydrodynamics (1.1) presented in the introduction.

To study the limit ε → 0 (the Euler limit ), it is natural to consider the following compressible Euler equation:

(2.4) ∂ t v + v • ∇v + ∇f (ρ) = 0 ; v |t=0 = ∇φ 0 , ∂ t ρ + div (ρv) = 0 ; ρ |t=0 = ρ 0 .
Note that the solution to (2.4) must not be expected to remain smooth for all time, even if the initial data are smooth. In [START_REF] Makino | Sur la solution à support compact de l'équation d'Euler compressible[END_REF] (see also [START_REF] Xin | Blowup of smooth solutions of the compressible Navier-Stokes equation with compact density[END_REF]), it is shown that compactly supported initial data lead to the formation of singularities in finite time. In [START_REF] Sideris | Formation of Singularities in Three-Dimensional Compressible Fluids[END_REF], the author constructs a solution developing singularities in finite time, in the absence of vacuum. In [START_REF] Alinhac | Temps de vie des solutions régulières des équations d'Euler compressibles axisymétriques en dimension deux[END_REF], it is shown that for rotationally invariant two-dimensional data that are perturbation of size ε of a rest state, blow-up occurs at time

T ε ∼ τ /ε 2 .
Remark 2.2. In the framework of this section, we have f ′ = 1 > 0. As pointed out above, the case f (r) = r σ , σ ∈ N, could be considered as well. On the other hand, if f ′ < 0 (corresponding to the semiclassical limit for a focusing nonlinearity), (2.4) becomes an elliptic system which may be solved locally-in-time for analytic data (see [START_REF] Gérard | Remarques sur l'analyse semi-classique de l'équation de Schrödinger non linéaire[END_REF][START_REF] Thomann | Instabilities for supercritical Schrödinger equations in analytic manifolds[END_REF]). At the same time, in the case d = 1 and f ′ = -1, it has been shown in [START_REF] Métivier | Remarks on the well-posedness of the nonlinear Cauchy problem, Geometric analysis of PDE and several complex variables[END_REF] that there are smooth initial data for which the Cauchy problem (2.4) has no solution. In short, working with analytic data in this context is not only convenient, it is mandatory.

Let us emphasize that the hydrodynamical formulations for the cubic NLS and the Gross-Pitaevskii equations are exactly the same as in both cases we have ∇f (ρ) = ∇ρ. From this point of view, studying either of the equations is mainly a matter of boundary conditions at infinity: H s (R d ) is the appropriate space for the cubic NLS equation whereas X s (R d ) is adapted to the Gross-Pitaevskii equation.

In the sequel, Z s denotes either H s (R d ) or X s (R d ). In addition, we set

Z ∞ = s>d/2 Z s . Theorem 2.3. Let ρ 0 , φ 0 ∈ C ∞ (R d ) with √ ρ 0 , ∇φ 0 ∈ Z s for some s > d/2 + 1.
There exists a unique maximal solution (v, ρ) ∈ C([0, T max ); Z s ) to (2.4). In addition, T max is independent of s > d/2 + 1 and

T max < +∞ =⇒ Tmax 0 (v, √ ρ)(t) W 1,∞ dt = +∞.
Finally, if ∇φ 0 and ρ 0 are smooth, nonzero and compactly supported then T max is finite.

We investigate the following natural questions:

Question 1. Assume that ρ 0 (x) > 0 for all x ∈ R d . Can we say that ρ ε (t, x) > 0 for all x ∈ R d and t ∈ [0, T max )? If not, what is the maximal interval allowed for t?

We will also recall that despite the appearance, the presence of vacuum (existence of zeroes of ψ ε ) is merely a technical problem: the Madelung transform ceases to make sense, but a rigorous WKB analysis is available, regardless of the presence of vacuum. See §2.3.

Question 2. Let ρ 0 , φ 0 ∈ C ∞ (R d ) with ρ 0 , ∇φ 0 ∈ H s (R d ) for some s > d/2 + 1. Suppose that the solution (v, ρ) to (2.4) satisfies ρ(t, x) > 0 for (t, x) ∈ [0, τ ] × R d .
Can we construct a solution to (1.1) in C([0, τ * ]; H s ) (possibly with 0 < τ * < τ )?

We will see that in general, the answer for this question is no. Even though from the answer to the first question, (1.1) makes sense formally, the analytical properties associated to (1.1) are not as favorable as for (2.4). Typically, the right-hand side of the equation for the quantum velocity need not belong to L 2 (R d ).

Question 3. Let ρ 0 , φ 0 ∈ C ∞ (R d ) with ρ 0 , ∇φ 0 ∈ X s (R d ) for some s > d/2 + 1.
Suppose that the solution (v, ρ) to (2.4) satisfies ρ(t, x) > 0 for (t, x) ∈ [0, τ ] × R d . Can we construct a solution to (1.1) in C([0, τ * ]; X s ) (possibly with 0 < τ * < τ )? 2.2. Proof of Theorem 2.3. We shall simply give the main ideas of the proof of Theorem 2.3. Complete proofs can be found in [START_REF] Makino | Sur la solution à support compact de l'équation d'Euler compressible[END_REF] for the cubic defocusing NLS equation in Sobolev spaces, and in [START_REF] Alazard | WKB analysis for the Gross-Pitaevskii equation with non-trivial boundary conditions at infinity[END_REF] for the Gross-Pitaevskii equation in Zhidkov spaces.

In the framework of this paper, we have

∇f (ρ) = ∇ρ.
Introduce formally the auxiliary function a = √ ρ. This nonlinear change of variable makes (2.4) hyperbolic symmetric:

(2.5)    ∂ t v + v • ∇v + 2a∇a = 0 ; v |t=0 = ∇φ 0 , ∂ t a + v • ∇a + 1 2 a∇ • v = 0 ; a |t=0 = √ ρ 0 .
This system is of the form

∂ t u + d j=1 A j (u)∂ j u = 0, where u =      v 1 . . . v d a      ∈ R d+1 ,
and the matrices A j are symmetrized by the constant multiplier

S = I d 0 0 4 .
Standard analysis (see e.g. [START_REF] Taylor | Partial differential equations[END_REF]) shows that (2.5) has a unique maximal solution (v, a) in C([0, T max ); Z s ), provided that s > d/2 + 1 and that, in addition,

T max < +∞ =⇒ Tmax 0 (v, a)(t) W 1,∞ dt = +∞.
We can then define ρ by the linear equation

∂ t ρ + div (ρv) = 0 ; ρ |t=0 = ρ 0 .
By uniqueness for this linear equation, ρ = a 2 (a is real-valued, so ρ is non-negative), and (v, ρ) solves (2.4).

We now briefly explain why compactly supported initial data lead to the formation of singularities in finite time. The first remark is that in this case, the solution to (2.4) has a finite speed of propagation, which turns out to be zero: so long as (v, ρ) is smooth, it remains supported in the same compact as its initial data. To see this, consider the auxiliary system (2.5): the first equation is a Burgers' equation with source term 2a∇a; the second equation is an ordinary differential equation along the trajectories of the particles. Define the trajectory by

(2.6) d dt x(t, y) = v (t, x(t, y)) ; x(0, y) = y.
For 0 t < T max , this is a global diffeomorphism of R d , as shown by the equation d dt ∇ y x(t, y) = ∇v (t, x(t, y)) ∇ y x(t, y) ; ∇ y x(0, y) = Id, and Gronwall lemma. Therefore, for a smooth function f ,

(∂ t f + v • ∇f ) (t, x(t, y)) = ∂ t (f (t, x(t, y))) ,
and (2.5) can be viewed as a system of ordinary differential equations.

Once the non-propagation of the support of smooth solutions is established, the end of the proof relies on a virial computation (like in [START_REF] Zakharov | Collapse of Langmuir waves[END_REF][START_REF] Glassey | On the blowing up of solutions to the Cauchy problem for nonlinear Schrödinger equations[END_REF], see also [START_REF] Cazenave | Semilinear Schrödinger equations[END_REF]). This computation shows that global in time smooth solutions to (2.5) are dispersive (see also [START_REF] Serre | Solutions classiques globales des équations d'Euler pour un fluide parfait compressible[END_REF]). This is incompatible with the zero propagation speed of smooth compactly supported solutions. Therefore, singularities have to appear in finite time.

2.3.

A review of WKB analysis associated to (2.1). We consider initial data which are a little more general than in (2.1), namely (2.7) ψ ε (0, x) = a ε 0 (x)e iϕ0(x)/ε , where the initial amplitude a ε 0 is assumed to be smooth, complex-valued, and possibly depending on ε. Typically, we assume that there exist a 0 , a 1 ∈ Z ∞ independent of ε such that

(2.8) a ε 0 = a 0 + εa 1 + O(ε 2 ) in Z s , ∀s > d/2.

First order approximation.

Introduce the solution to the quasilinear system (2.9)

     ∂ t ϕ + 1 2 |∇ϕ| 2 + f |a| 2 = 0 ; ϕ |t=0 = ϕ 0 , ∂ t a + ∇ϕ • ∇a + 1 2 a∆ϕ = 0 ; a |t=0 = a 0 .
Theorem 2.3 shows that (2.9) has a unique, smooth solution with a, ∇ϕ ∈ Z ∞ . The main remark consists in noticing that (2.9) implies that (∇ϕ, |a| 2 ) has to solve (2.4) (a 0 may be complex-valued): Theorem 2.3 yields v, ρ ∈ Z ∞ . We can then define a as the solution to the linear transport equation

∂ t a + v • ∇a + 1 2 a∇ • v = 0 ; a |t=0 = a 0 .
Now |a| 2 and ρ solve the same linear transport equation, with the same initial data, hence ρ = |a| 2 . Using this information in the equation for the velocity, define

ϕ(t) = ϕ 0 - t 0 1 2 |v(τ )| 2 + f |a(τ )| 2 dτ.
We easily check that ∂ t (∇ϕ -v) = ∇∂ t ϕ -∂ t v = 0, and that (ϕ, a) solves (2.9). Introduce the solution to the linearization of (2.9), with an extra source term:

     ∂ t ϕ (1) + ∇ϕ • ∇ϕ (1) + 2 Re aa (1) f ′ |a| 2 = 0 ; ϕ (1) |t=0 = 0, ∂ t a (1) + ∇ϕ • ∇a (1) + ∇ϕ (1) • ∇a + 1 2 a (1) ∆ϕ + 1 2 a∆ϕ (1) = i 2 ∆a ; a (1) |t=0 = a 1 .
We also check that it has a unique smooth solution, with a (1) , ∇ϕ (1) ∈ Z ∞ . The main result that we will invoke is the following:

Proposition 2.4. Let a ε 0 , a 0 , a 1 , ϕ 0 be smooth, with a ε 0 , a 0 , a 1 , ∇ϕ 0 ∈ Z ∞ . Assume that (2.8) holds. Then (2.10) ψ ε = ae iϕ (1) + O (ε) e iϕ/ε in L ∞ ([0, τ ]; Z s ), ∀τ < T max and ∀s 0. This result was established in [START_REF] Grenier | Semiclassical limit of the nonlinear Schrödinger equation in small time[END_REF] when Z s = H s (R d ), and in [START_REF] Alazard | WKB analysis for the Gross-Pitaevskii equation with non-trivial boundary conditions at infinity[END_REF] when Z s = X s (R d ). Note the shift between the order of the approximation between the initial data (known up to O(ε 2 )) and the approximation (of order O(ε) only): this is due to the fact that we consider a regime which is super-critical as far as WKB analysis is concerned (see e.g. [START_REF] Carles | Semi-classical analysis for nonlinear Schrödinger equations[END_REF]). In particular, the phase modulation ϕ (1) is a function of ϕ 0 , a 0 and a 1 . It is non-trivial in general, and since we are interested here in real-valued a 0 , we shall merely mention two cases (see [16, pp. 69-70]):

• If a 1 = 0 is real-valued, then ϕ (1) = 0 in general.

• If a 1 = 0 (or more generally if a 1 ∈ iR), then ϕ (1) = 0. To see the first point, it suffices to notice that the equation for ϕ (1) gives (recall that f ′ = 1):

∂ t ϕ (1)
|t=0 = -2a 0 a 1 . This shows that for (2.2)-(1.1) to yield a relevant description of the solution to (2.1), we have to assume a 0 = √ ρ 0 and a 1 = 0. Otherwise, a phase modulation is necessary to describe ψ ε at leading order, by (2.10), which is incompatible with the form (2.2), unless the Madelung phase φ ε admits a corrector of order ε. But formal asymptotics in (1.1) give v ε = v +O(ε 2 ) = ∇ϕ+O(ε 2 ). Hence the Madelung transform has a chance to give a relevant result only if a 1 = 0.

To check the second point of the above assertion, we set α = Re aa (1) . Direct computations show that (ϕ (1) , α) solves, as soon as a 0 ∈ R and a 1 ∈ iR:

   ∂ t ϕ (1) + ∇ϕ • ∇ϕ (1) + 2α = 0 ; ϕ (1) |t=0 = 0, ∂ t α + ∇ϕ • ∇α = - 1 2 div |a| 2 ∇ϕ (1) -α∆ϕ ; α |t=0 = 0.
This is a linear, homogeneous system, with zero initial data, so its solution is identically zero.

To conclude this paragraph, we briefly outline the proof of Proposition 2.4. The approach in [START_REF] Alazard | WKB analysis for the Gross-Pitaevskii equation with non-trivial boundary conditions at infinity[END_REF] is the same as in [START_REF] Grenier | Semiclassical limit of the nonlinear Schrödinger equation in small time[END_REF], with slightly different estimates. For simplicity, and in view of the above discussion, we assume a ε 0 = a 0 independent of ε. We write the solution ψ ε as ψ ε = a ε e iϕ ε /ε (exact formula), where we impose

(2.11)      ∂ t ϕ ε + 1 2 |∇ϕ ε | 2 + f |a ε | 2 = 0 ; ϕ ε |t=0 = ϕ 0 , ∂ t a ε + ∇ϕ ε • ∇a ε + 1 2 a ε ∆ϕ ε = i ε 2 ∆a ε ; a ε |t=0 = a 0 .
Note that both ϕ ε and a ε depend on ε, because the right-hand side of the equation for a ε depends on ε, and because of the coupling between the two equations. Note also that with this approach, one abandons the possibility of considering a realvalued amplitude a ε . It is not hard to construct a solution to (2.11) in Z s , for s > d/2 + 2, and then check that asymptotic expansions are available in Z s :

a ε = a + εa (1) + O ε 2 ; ϕ ε = ϕ + εϕ (1) + O ε 2 .
Back to ψ ε , this yields Proposition 2.4. We see that the general loss in the precision (from O(ε 2 ) in the initial data to O(ε) in the approximation for t > 0) is due to the division of ϕ ε by ε. Note finally that even though a 1 = 0, one has a (1) = 0: the corrector a (1) solves a linear equation, with a purely imaginary (non trivial) source term, and so a (1) ∈ iR is not trivial, while ϕ (1) = 0 since Re aa (1) = 0.

Higher order approximation and formal link with quantum hydrodynamics.

One can actually consider an asymptotic expansion to arbitrary order,

a ε = a + εa (1) + . . . + ε N a (N ) + O ε N +1 , ϕ ε = ϕ + εϕ (1) + . . . + ε N ϕ (N ) + O ε N +1 , ∀N ∈ N.
For j 1, the coefficient a (j) is given by a linear system for (ϕ (j) , a (j) ), with source terms involving (ϕ (k) , a (k) ) 0 k j-1 . In the case a 0 ∈ R and a 1 = 0, we know that ϕ (1) = 0, and a (2) is given by

∂ t a (2) + ∇ϕ • ∇a (2) + 1 2 a (2) ∆ϕ + ∇ϕ (2) • ∇a + 1 2 a∆ϕ (2) = i 2 ∆a (1) ; a (2) |t=0 = 0.
We check that in the case a ε 0 = a 0 ∈ R (which includes the case where Madelung transform is used, a 0 0), all the profiles a and a (2j) , j 1, are real-valued, while a (2j+1) , j 0, are purely imaginary. Moreover, ϕ (2j+1) = 0 for all j ∈ N. This is formally in agreement with (1.1): indeed, (1.1) suggests that φ ε and ρ ε have asymptotic expansions of the form (2.12)

φ ε ≈ φ+ ε 2 φ (2) + . . .+ ε 2j φ (2j) + . . . ; ρ ε ≈ ρ+ ε 2 ρ (2) + . . .+ ε 2j ρ (2j) + . . .
On the other hand, we have

ρ ε = |a ε | 2 ≈ a + εa (1) + . . . 2 ≈ a + . . . + ε 2j a (2j) + . . . 2 -εa (1) + . . . + ε 2j+1 a (2j+1) + . . . 2 ,
since the a (2j+1) 's are purely imaginary. This is in agreement with the second formal asymptotics in (2.12). We can check similarly that (2.12) is in agreement with the higher order generalization of (2.10), in view of the special properties of the ϕ (j) 's and a (j) 's pointed out above.

2.4. Absence of vacuum before shocks. Lemma 2.5. In Theorem 2.3, assume that

ρ 0 (x) > 0 for all x ∈ R d (absence of vacuum). Then ρ > 0 on [0, T max ) × R d .
Proof. As in Section 2.2, we use the fact that on [0, T max ), the equation for the density is just an ordinary differential equation. Introduce the Jacobi determinant J t (y) = det∇ y x(t, y), where x(t, y) is given by (2.6). We have seen that

J t (y) > 0 for (t, y) ∈ [0, T max ) × R d . Change the unknown ρ to r, with r(t, y) = ρ (t, x(t, y)) J t (y).
Then for 0 t < T max , the continuity equation is equivalent to: ∂ t r = 0. Therefore,

ρ(t, x) = 1 J t (y(t, x)) ρ 0 (y(t, x)) ,
where x → y(t, x) denotes the inverse mapping of y → x(t, y).

We infer:

Proposition 2.6. Under the assumptions of Theorem 2.3, assume that ρ 0 (x) > 0 for all x ∈ R d (absence of vacuum). Let 0 < T < T max , and K be a compact set in R d . There exists ε(T, K) > 0 such that for 0 < ε ε(T, K),

|ψ ε | > 0 on [0, T ] × K.
Proof. Proposition 2.4 shows that

|ψ ε | = |a| + O(ε) in L ∞ [0, T ] × R d .
Note that the constant involved in this O(ε) depends on T in general. Recalling that a = √ ρ, Lemma 2.5 shows that min

(t,x)∈[0,T ]×K a(t, x) = c(T, K) > 0. Now for 0 < ε ε(T, K) ≪ 1, ||ψ ε | -|a|| 1 2 c(T, K) in L ∞ [0, T ] × R d ,
and the result follows.

In the case x ∈ T d , this shows that before the formation of shocks in the Euler equation, and provided that ε is sufficiently small, the amplitude remains positive: the right-hand side of (1.1) makes sense. This point was remarked initially in [START_REF] Gérard | Remarques sur l'analyse semi-classique de l'équation de Schrödinger non linéaire[END_REF]. Note that the result of [START_REF] Gamba | On the blowing up of solutions to quantum hydrodynamic models on bounded domains[END_REF] in the one-dimensional case x ∈ [0, 1] shows that suitable boundary conditions lead to the existence of finite time blow-up for (1.1). Therefore, the above result is qualitatively sharp (qualitatively only, for it might happen that the solution to (1.1) remains smooth longer than the solution to (2.4)).

Finally, we show that the compactness assumption in Proposition 2.6 can be removed in the case of the Gross-Pitaevskii equation. As regards the nonlinear Schrödinger equation on R d , this issue seems much more delicate and will not be addressed in this paper. Assume that the Gross-Pitaevskii equation is associated with the boundary condition at infinity

|ψ ε (t, x)| -→ |x|→∞ 1.
Such a condition is used frequently in physics, possibly with a stronger one, of the form (see e.g. [START_REF] Lin | Semiclassical limit of the Gross-Pitaevskii equation in an exterior domain[END_REF] and references therein)

ψ ε (t, x) -e iv ∞ •x/ε -→ |x|→∞ 0,
for some fixed asymptotic "velocity" v ∞ ∈ R d . Now, as regards the Gross-Pitaevskii equation, putting together the continuity of ψ ε over [0, T ] × R d , the compactness of the time interval [0, T ] and Proposition 2.6, we get Corollary 2.7. Under the assumptions of Theorem 2.3, assume that ρ 0 (x) > 0 for all x ∈ R d (absence of vacuum). Assume moreover that the Gross-Pitaevskii equation is associated with the boundary condition at infinity

2 |ψ ε (t, x)| -→ |x|→∞ 1.
Let 0 < T < T max . There exist ε(T ), c(T ) > 0 such that

|ψ ε | c(T ) on [0, T ] × R d , for all 0 < ε ε(T ).
2.5. Functional spaces associated to the Madelung transform. It is rather easy to see that the answer to Question 2 is no, in general. Consider for ρ 0 the function in the Schwartz class

ρ 0 (x) = e -|x| 2k , k ≥ 1.
At time t = 0, the quantum pressure (right-hand side of (1.1)) grows like |x| 2k-1 hence the velocity v ε has no chance to belong to H s for general initial data in H ∞ . Thus, working in Sobolev spaces for general initial data does not make sense for (1.1), while the results in [START_REF] Grenier | Semiclassical limit of the nonlinear Schrödinger equation in small time[END_REF] show that it is a fairly reasonable framework to study the semi-classical limit of (2.1).

On the other hand, like for the absence of vacuum, the answer to Question 3 is positive, at least if we consider some special boundary conditions at infinity.

Solving the QHD system by a direct approach

In the present section, we describe an efficient method to solve directly the hydrodynamic form of (2.1) given by (1.1), once performed the Madelung transform. This method enables us to study the corresponding initial value problem for (1.1) with data (v 0 , ρ 0 ) such that (v 0 , ∇ρ 0 ) has a high order Sobolev regularity and ρ 0 is positive and bounded away from zero. In addition to local-in-time well-posedness results, we get (see Theorems 3.1 and 3.3 below) nontrivial lower bounds on the first appearance of a zero for the solution, which are of particular interest for the study of long-wavelength asymptotics if the data are a perturbation of a constant state of modulus one.

We here closely follow the approach that has been initiated in [START_REF] Béthuel | On the linear wave regime of the Gross-Pitaevskii equation[END_REF]. To help the reader to compare the present results with those of the previous section however, we keep on using the semi-classical scaling given by (2.1) (whereas ε = 1 in [START_REF] Béthuel | On the linear wave regime of the Gross-Pitaevskii equation[END_REF]).

The use of a suitable extended formulation for (1.1) and of weighted Sobolev estimates will be the key to our approach. Let us stress that, recently, similar extended formulations have proved to be efficient in other contexts for both numerical (see [START_REF] Coquel | Extended formulations for Van der Waals models[END_REF]) and theoretical purposes. As a matter of fact, in the last paragraph of this section, we shall briefly explain how our approach based on such an extended formulation carries over to the more complicated case of Korteweg fluids.

3.1. Solving the QHD system by means of an extended formulation. The "improved" WKB method that has been described in the previous section amounts to writing the sought solution ψ ε as:

ψ ε = a ε e iϕ ε /ε
for some complex valued function a ε and real valued function ϕ ε . In this section, we rather start from the Madelung transform

ψ ε = √ ρ ε e iφ ε where √ ρ ε = |ψ ε |,
then write

ψ ε = e iΦ ε /ε with Φ ε = φ ε -i ε 2 log ρ ε ,
and consider the redundant system that is satisfied by both ρ ε and z ε = ∇Φ ε = v ε + iw ε with v ε = ∇φ ε and w ε = -ε 2ρ ε ∇ρ ε . In order to obtain the system for z ε , we first differentiate the density equation in (1.1). This yields

∂ t w ε + ∇(v ε • w ε ) = ε 2 ∇ div v ε .
Next, we notice that

ε 2 2 ∆ √ ρ ε √ ρ ε = - ε 2 div w ε + 1 2 |w ε | 2 .
In consequence, the equation for v ε rewrites

∂ t v ε + v ε • ∇v ε - 1 2 ∇|w ε | 2 + ε 2 ∇ div w ε + ∇f (ρ ε ) = 0.
Of course, z ε is a potential vector-field, hence ∇ div z ε = ∆z ε so that we eventually get the following "extended" system for (ρ ε , z ε ):

(3.1)

∂ t z ε + 1 2 ∇(z ε • z ε ) + ∇f (ρ ε ) = i ε 2 ∆z ε , ∂ t ρ ε + div(ρ ε v ε ) = 0,
where we agree that a • b := Let us now explain how Sobolev estimates may be derived for (ρ ε , z ε ) in the case where ρ ε = 1 + b ε for some b ε going to 0 at infinity, and f (ρ ε ) = ρ ε (to simplify). The following computations are borrowed from [START_REF] Béthuel | On the linear wave regime of the Gross-Pitaevskii equation[END_REF]. For notational simplicity, we omit the superscripts ε.

In order to get the basic energy estimate, we compute3 :

d dt R d (1 + b)|z| 2 + b 2 = 2 R d (1 + b) z, ∂ t z I1 + R d b ∂ t b I2 + R d ∂ t b |z| 2 I3
where the notation a, b

= d j=1
Re a j Re b j + Im a j Im b j has been used in I 1 .

Further computations yield I 1 = I 1,1 + I 1,2 + I 1,3 + I 1,4 + I 1,5 and I 2 = I 2,1 + I 2,2 with

I 1,1 = - z, ∇b , I 2,1 = -b div v, I 1,2 = -b z, ∇b , I 2,2 = -b div(bv). I 1,3 = z, i ε 2 ∆z , I 1,4 = b z, i ε 2 ∆z , I 1,5 = - 1 2 ρ z, ∇(z • z) ,
Using obvious integrations by parts we readily get

I 1,1 + I 2,1 = 0, I 1,2 + I 2,2 = 0 and I 1,3 = 0.
Therefore, integrating by parts in I 1,4 also, we get

d dt (1 + b)|z| 2 + b 2 = -ε z, i∇z • ∇b -2 ρ z, ∇z • z + ∂ t b |z| 2 .
For "general" functions b and z, the appearance of the terms ∇b and ∇z would preclude any attempt to "close" the estimates. In our case however, as the algebraic relation -ε∇b = 2ρw holds true, one may avoid this loss of one derivative for one may write

-ε z, i∇z • ∇b -2ρ z, ∇z • z = -2ρ z, ∇z • v .
Now, integrating by parts an ultimate time, we conclude that

d dt ρ|z| 2 + b 2 = -ρv • ∇|z| 2 + ∂ t b |z| 2 = 0.
Hence, ρ|z| 2 + b 2 is a conserved quantity.

The same algebraic cancellations may be used for getting higher order Sobolev (or Besov) estimates. Indeed consider an "abstract" pseudo-differential operator A(D) (for instance a differential, a fractional derivatives or a spectral localization operator). Then one may write

d dt (1 + b)|A(D)z| 2 + (A(D)b) 2 = 2 (1 + b) A(D)z, ∂ t A(D)z I1 + A(D)b ∂ t A(D)b I2 + ∂ t b |A(D)z| 2

I3

.

We notice that I 1 = I 1,1 + I 1,2 + I 1,3 + I 1,4 + I 1,5 and I 2 = I 2,1 + I 2,2 with

I 1,1 = - A(D)z, ∇A(D)b , I 2,1 = -A(D)b div A(D)v, I 1,2 = -b A(D)z, ∇A(D)b , I 2,2 = -A(D)b A(D) div(bv). I 1,3 = A(D)z, i ε 2 ∆A(D)z , I 1,4 = b A(D)z, i ε 2 ∆A(D)z , I 1,5 = - 1 2 ρ A(D)z, A(D)∇(z • z) ,
As above, obvious integrations by parts ensure that I 1,1 + I 2,1 = 0 and I 1,3 = 0. Next, using again integrations by parts, we notice that

I 1,2 + I 2,2 = A(D)b div[b, A(D)]v.
Finally, integrating by parts in I 1,4 and using the fact that -ε∇b = 2ρw yields

I 1,4 = A(D)z, i(∇A(D)z) • (ρw) ,
and we have

I 1,5 = -A(D)z, A(D)∇z • (ρv + iρw) + ρ A(D)z, ∇A(D)z • z -A(D)(∇z • z) .
Therefore, using the fact that ∂ t ρ + div(ρv) = 0, we conclude that 2(I 1,4 + I 1,5

) + I 3 = 2 ρ A(D)z, ∇A(D)z • z -A(D)(∇z • z) .
Putting all the above equalities together, we thus get 1 2

d dt ρ|A(D)z| 2 + (A(D)b) 2 = A(D)b div[b, A(D)]v + ρ A(D)z, ∇A(D)z • z -A(D)(∇z • z) .
If, say, A(D) is a fractional derivatives operator, then one may show by means of classical commutator estimates that the right-hand side may be bounded by

C ρ L ∞ (Db, Dz) L ∞ (A(D)b, A(D)z) 2 L 2 . Therefore, d dt ρ|A(D)z| 2 + (A(D)b) 2 C ρ L ∞ (Db, Dz) L ∞ (A(D)b, A(D)z) 2 L 2 .
Denoting by E 2 A (t) the left-hand side and resorting to Gronwall lemma, we thus get

(3.2) E A (t) E A (0) exp C t 0 ρ L ∞ ρ -1 L ∞ (Db, Dz) L ∞ dτ .
It is now clear that whenever Db and Dz are bounded in L 1 ([0, T ]; L ∞ ) and ρ is bounded from below and from above then we get a control of A(D)b and

A(D)z in L ∞ ([0, T ]; L 2 ).
Taking A(D) = Λ s and performing a time integration, (3.2) implies that

(3.3) (b, z)(t) H s ≤ C (b 0 , z 0 ) H s + T 0 (Db, Dz) L ∞ (b, z)(t) H s dt for some constant C = C(s, d, ρ ±1 L ∞ ).
So assuming that s > 1 + d/2 and using Sobolev embedding and Gronwall's inequality, one may conclude by elementary methods to the following statement.

Theorem 3.1. Let s > 1 + d/2. Assume that ρ 0 = 1 + b 0 for some b 0 ∈ H s+1 (R d ) such that 1 + b 0 > 0, and that v 0 ∈ H s (R d ).
Then there exists a time

T T 0 := C b 0 H s + ε ∇b 0 H s + v 0 H s+1 with C = C(s, d, ρ ±1 0 L ∞ ) such that (1.1) has a unique solution (v ε , ρ ε ) on [0, T ] × R d with ρ ε = 1 + b ε bounded away from 0 and (v ε , b ε ) ∈ C([0, T ]; H s × H s+1 ) ∩ C 1 ([0, T ]; H s-2 × H s-1 ).
Remark 3.2. Combining basic energy estimates for the wave equation with the above result, one may control the discrepancy between (b, v) and the solution to the acoustic wave equation

(3.4) ∂ t v + ∇ ḃ = 0 ; v|t=0 = v 0 , ∂ t ḃ + div v = 0 ; ḃ|t=0 = b 0 .
We have, up to time T 0 ,

(v ε -v, b ε -ḃ)(t) H s-2 ≤ C t (b 0 , u 0 ) 2 H s+1 ×H s + ε 2 t (b 0 , u 0 ) H s+1 ×H s . Note also that (3.
3) provides a blow-up criterion involving the W 2,∞ norm of b and the Lipschitz norm of v. In particular, this implies that for given data in H s (s > 1 + d/2), the lifespan in H s is the same as the lifespan in H s ′ , for any 1 + d/2 < s ′ < s.

Dispersive properties and improved lower bounds for the lifespan.

The system for (b, v) reads

(3.5)    ∂ t v + v • ∇v + ∇b = ε 2 4 ∇ 1 ρ ∆b - 1 2ρ 2 |∇b| 2 , ∂ t b + div v = -div(bv).
Therefore the linearized system about (0, 0) is not (3.4) but rather

(3.6)    ∂ t v + ∇ ḃ = ε 2 4 ∇∆ ḃ, ∂ t ḃ + div v = 0.
A straightforward spectral analysis (based on the Fourier transform) shows that the above linear system behaves as the wave equation with speed 1 for frequencies small with respect to 1/ε, and as the Schrödinger equation with coefficient ε/2 in the high frequency regime. In fact, in dimension d 2, it is possible to prove Strichartz inequalities (related to the wave and Schrödinger equations for low and high frequencies, respectively) for (3.6). In the case of small data (b 0 , v 0 ), these inequalities allow to improve the lower bound for the lifespan (see also [START_REF] Danchin | Zero Mach number limit in critical spaces for compressible Navier-Stokes equations[END_REF] where a similar idea has been used in the context of the incompressible limit for compressible flows). For the sake of simplicity, let us just state the result in dimension d 4 (the reader is referred to [START_REF] Béthuel | On the linear wave regime of the Gross-Pitaevskii equation[END_REF] for the case d = 2, 3 and for more details concerning the approximation of the solution by (3.6)): Theorem 3.3. Under the assumptions of Theorem 3.1 with s > 2 + d/2, then the lifespan T may be bounded from below by

T 1 := C ( b 0 H s + ε ∇b 0 H s + v 0 H s+1 ) 2 ,
and the discrepancy between (b, v) and the solution ( ḃ, v) to (3.6) with the same data may be bounded in terms of t and of the data, up to time T 1 .

3.3.

Extended formulation for Korteweg fluids. Compared to the "improved" WKB method, the main drawback of the direct approach based on an extended formulation for solving (1.1) is that vanishing solutions cannot be handled.

On the other hand, the direct method is robust enough so as to be used to solve locally more complicated models such as the following system governing the evolution of inviscid capillary fluids:

(3.7) ∂ t v + v • ∇v + ∇f (ρ) = ∇ κ(ρ)∆ρ + 1 2 κ ′ (ρ)|∇ρ| 2 , ∂ t ρ + div(ρv) = 0.
Physically, the function κ correspond to the capillary coefficient. Obviously, System (1.1) is included in (3.7) (take κ(ρ) = ε 2 /(4ρ)). In the general case, introducing

a(ρ) = ρκ(ρ), w = - κ(ρ) ρ ∇ρ and z = v + iw,
we get the following extended formulation for (3.7):

(3.8)

∂ t z + v • ∇z + i∇z • w + ∇f (ρ) = i∇(a(ρ) div z), ∂ t ρ + div(ρv) = 0.
Note that in the potential case (namely curl z = 0) then v • ∇z + i∇z • w = 1 2 ∇(z • z). Note also that in the general case, the second order term ∇(a(ρ) div z) is degenerate.

The case of System (1.1) is particularly simple inasmuch as a is the constant function ε/2 and ∇ div z = ∆z.

For general capillarity coefficients, one may prove a local well-posedness result, similar to that of Theorem 3.1. This has been done in [START_REF] Benzoni-Gavage | On the well-posedness for the Euler-Korteweg model in several space dimensions[END_REF]. The proof relies on the use of weighted Sobolev estimates, with a weight depending both on ρ and on the order of differentiation.

As for the QHD system, we expect the potential part of the solution to System (3.7) to have dispersive properties in dimension d 2. The general situation is much more complicated however, because those properties are related to those of the quasilinear Schrödinger equation. To our knowledge, this aspect has been investigated only very recently in a work by C. Audiard [START_REF] Audiard | Dispersive smoothing for the Euler-Korteweg model[END_REF] that concerns potential flows.

Asymptotics for the Gross-Pitaevskii equation

This section is concerned with the existence and asymptotics of traveling wave solutions for the Gross-Pitaevskii equation (4.1)

iψ t + ∆ψ + (1 -|ψ| 2 )ψ = 0,
which may be obtained from (2.1) (with f (r) = r -1), up to the factor 1 2 , after performing the change of unknown:

ψ ε (t, x) = ψ(ε -1 t, ε -1 x).
Equation (4.1) is associated to the Ginzburg-Landau energy (or Hamiltonian):

(4.2) H(ψ) = 1 2 R d |∇ψ| 2 + 1 4 R d (1 -|ψ| 2 ) 2 .
As a consequence, in contrast with the cubic NLS equation, the natural energy space for (4.1) is not H 1 (R d ) but rather

E(R d ) = ψ ∈ H 1 loc (R d ), s.t. H(ψ)
< +∞ • As pointed out before, for H(ψ) to be finite, |ψ| must, in some sense, tend to 1 at infinity. This "nontrivial" boundary condition provides (4.1) with a richer dynamics than in the case of null condition at infinity which, for a defocusing NLS type equation, is essentially governed by dispersion and scattering. For instance, in nonlinear optics, the "dark solitons" are localized nonlinear waves (or "holes") which exist on a stable continuous wave background. The boundary condition |ψ(t, x)| → 1 at infinity is due to this nonzero background. In the context of superfluids, 1 is the density of the fluid at infinity.

Similarly to the energy, the momentum

P(ψ) = 1 2 R d i∇ψ , ψ ,
is formally conserved. This quantity is well defined for ψ ∈ H 1 (R d ) but not for solutions with a finite Ginzburg-Landau energy. A major difficulty in the theory of the Gross-Pitaevskii equation is to find an appropriate definition of the momentum which leads to a conserved quantity. A natural definition would be

P(ψ) = 1 2 R d i∇ψ , ψ -1 , but this would require for instance that ψ -1 ∈ L 2 (R d ).
In any case we will denote by p the (scalar) first component of P which will play an important role in this section.

Recall that the use of the Madelung transform4 ψ = ̺e iφ , leads to the following hydrodynamic form of the equation for ̺ and v = 2∇φ :

(4.3)      ∂ t v + v • ∇v + 2∇̺ 2 = 2∇ ∆̺ ̺ , ∂ t ̺ 2 + div (̺ 2 v) = 0.
As pointed out before, if we discard the right-hand side of the first equation and look at ̺ 2 as the density of a fluid with velocity v, then the above system coincides with the Euler equations for a compressible fluid with pressure law P (ρ) = ρ 2 . In particular, the speed of sound waves near the constant solution v = 1 is given by

c s = √ 2.
As we will see below, this sound speed (the value and relevance of which is not so obvious if looking at the original Gross-Pitaevskii equation (4.1)), plays an important role in various aspects of the dynamics of (4.1). The value of c s may be also found by neglecting the quantum pressure term and linearizing for a perturbation ψ = (1 + ̺) exp(i φ). This leads to the wave equation:

∂ 2 t ̺ -2∆̺ = 0.
Note that if the quantum pressure is included (as in (3.6)) then the linearization reads:

∂ 2 t ̺ -2∆̺ -∆ 2 ̺ = 0, which, roughly, is the factorization of two linear Schrödinger operators.

The Madelung transform and the hydrodynamic form of (4.1) turn out to be of great interest to study the Gross-Pitaevskii equation with finite Ginzburg-Landau energy since the solution is expected to have very few "vortices" (or cancellations). Even in the "Euler limit" that has been presented in Section 2, one can use it outside the vortices ( [START_REF] Béthuel | Vortex rings for the Gross-Pitaevskii equation[END_REF][START_REF] Béthuel | Travelling waves for the Gross-Pitaevskii equation. I[END_REF]) to study the traveling waves of sufficiently small velocities. Let us also stress that the hydrodynamic form of the (one dimensional) Gross-Pitaevskii equation is needed in order to define a generalized momentum in the context of the orbital stability of the black solitons (such solitary waves have zeroes. . . ); see [START_REF] Béthuel | Orbital stability of the black soliton for the Gross-Pitaevskii equation[END_REF].

In the present section, we shall concentrate on the transonic limit of solutions to the Gross-Pitaevskii equation. We shall first present a result pertaining to the asymptotics of traveling waves with speed c tending to the sound speed c s , in the case d = 2, in connexion with the (KP I) equation (see below). Next, for the onedimensional case, we give an accurate description of the transonic long wave limit of (4.1) in terms of solutions to the KdV equation.

4.1. The transonic limit of finite energy traveling waves. Finite energy traveling wave solutions of (4.1) are solutions of the form ψ(x, t) = θ(x 1 -ct, x ⊥ ) where H(θ) < +∞ and x ⊥ denotes the transverse variables x 2 , • • • , x d . The profile θ satisfies the following equation:

(4.4) ic∂ 1 θ + ∆θ + θ(1 -|θ| 2 ) = 0.
A suitable functional setting for the study of such traveling waves is the space:

W (R d ) = {1} + V (R d ), with V (R d ) = {ψ : R d → C, (∇ψ, Re ψ) ∈ L 2 (R d ) 2 , Im ψ ∈ L 4 (R d ), ∇ Re ψ ∈ L 4/3 (R d )}• Indeed, given that W (R d
) is a subset of the energy space E(R d ), for any data in W (R d ), Equation (4.1) admits a unique solution. In addition, one may show that this solution stays in W (R d ) (see [START_REF]The Cauchy problem for the Gross-Pitaevskii equation[END_REF][START_REF]The Gross-Pitaevskii equation in the energy space, Stationary and time dependent Gross-Pitaevskii equations[END_REF]). Furthermore, the quantity i∂ 1 ψ, ψ-1 is integrable whenever ψ ∈ W (R d ), so that the scalar momentum p(ψ) is well-defined. This is a consequence of the identity

(4.5) i∂ 1 ψ, ψ -1 = ∂ 1 (Re ψ) Im ψ -∂ 1 (Im ψ)(Re ψ -1),
and various Hölder's inequalities. So finally, H and p are continuous on W (R d ) and all finite energy subsonic solutions to (4.4) have to belong to W (R d ). Moreover, if ψ ∈ W (R d ) may be lifted as ψ = ̺ exp(iφ) then

(4.6) p(ψ) = 1 2 R d i∂ 1 ψ , ψ -1 = 1 2 R d (1 -̺ 2 )∂ 1 φ.
Notice that for maps which may be lifted, with ̺ 1 2 , the last integral makes sense, even if we assume that ψ only belongs to the energy space E(R d ).

To simplify the presentation, we shall focus on the simpler case d = 2. We shall denote x = x 1 and y = x ⊥ = x 2 .

It is proven in [START_REF] Béthuel | Existence and properties of travelling waves for the Gross-Pitaevskii equation, Stationary and time dependent Gross-Pitaevskii equations[END_REF] that, for any p > 0, the following minimization problem

(4.7) H min (p) = inf{H(ψ), ψ ∈ W (R 2 ), p(ψ) = p},
has a solution u p which is a nontrivial traveling wave. We call it a ground state.

In the rest of this subsection, we focus on the asymptotics p going to 0 for u p , in connexion with the Kadomtsev-Petviashvili I (KP I) equation

u t + uu x + u xxx -∂ -1
x u yy = 0, where the antiderivative is defined in Fourier variables by ∂ -1

x f (ξ) = 1 iξ f (ξ). The following proposition states that the corresponding speed c(u p ) tends to c s and gives the first term in the asymptotic expansion for both H min (p) and c(u p ). Proposition 4.1. There exist positive constants p 1 , K 0 and K 1 such that we have the asymptotic behaviors

(4.8) 48 √ 2 S 2 KP p 3 -K 0 p 4 √ 2p -H min (p) K 1 p 3 , ∀0 p p 1 ,
where S KP stands for the action of the KP I ground state N of velocity 1, that is

(4.9) S KP = 1 2 R 2 (∂ x N ) 2 + 1 2 R 2 (∂ -1 x (∂ y N )) 2 - 1 6 R 2 N 3 + 1 2 R 2 N 2 .
Moreover, the map u p has no zeroes and there exist some positive constants p 2 , K 2 and K 3 such that

(4.10) K 2 p 2 √ 2 -c(u p ) K 3 p 2 , ∀0 p < p 2 .
Actually, it was established in [START_REF] Béthuel | Existence and properties of travelling waves for the Gross-Pitaevskii equation, Stationary and time dependent Gross-Pitaevskii equations[END_REF] that if ψ is a finite energy traveling wave of sufficiently small energy, then

1 2 < |ψ| < 1.
This implies that small energy traveling waves have no zeroes and thus can be lifted according to the Madelung transformation. This is in particular the case of minimizers u p corresponding to small enough values of p.

The transonic limit to KP I for traveling waves is obtained through the following change of scales:

x = ǫ(ψ)x, ỹ = ǫ(ψ) 2 √ 2 y.
We then set

η = 1 -|ψ| 2 , N ψ (x, y) = 6 ǫ(ψ) 2 η x ǫ(ψ) , √ 2y ǫ(ψ) 2 .
It turns out that N ψ converges to a traveling wave solution of the (differentiated) KP I equation 5 as ǫ(ψ) → 0, that is it approximately solves the equation

Θ p (x, y) = 6 √ 2 ǫ p φ p x ǫ p , √ 2y ǫ 2 p ,
we have (see [START_REF]On the KP I transonic limit of two-dimensional Gross-Pitaevskii travelling waves[END_REF]):

Theorem 4.2. There exist a subsequence (p n ) n∈N tending to 0, as n → +∞, and a ground state w of the KP I equation such that both N pn and Θ pn tend to w in W k,q (R 2 ) (for any k ∈ N and q ∈ (1, +∞]) as n goes to +∞.

Remark 4.3. Those results on the transonic limit of solitary waves have been recently extended in [START_REF] Chiron | Rarefaction pulses for the nonlinear Schrödinger equation in the transonic limit[END_REF] to the three-dimensional case with also general nonlinearities (see also [START_REF] Mariş | Travelling waves for nonlinear Schrödinger equations with non zero boundary conditions at infinity[END_REF] for other existence results). The additional serious difficulty is that the Gross-Pitaevskii ground states solutions are no longer global minimizers of the energy with fixed momentum and thus not expected to be stable.

4.2. The unsteady transonic long wave limit of the Gross-Pitaevskii equation. The Madelung transform is also crucial to derive and justify the transonic (weak amplitude, long wave) limit of the Gross-Pitaevskii equation.

With the scaling which is used in this section and for data which are perturbations of order ε 2 of a constant state with modulus 1, Theorem 3.1 and the remark that follows ensure that the linear wave equation gives a good approximation of the solution for t = o(ε -3 ). In the present paragraph, we describe what happens at next order. We shall see in particular that, up to times of order ε -3 , the Korteweg-de Vries (KdV) equation u t + uu x + u xxx = 0, gives an accurate approximation of the one-dimensional Gross-Pitaevskii equation (4.12)

i∂ t ψ + ∂ xx ψ + (1 -|ψ| 2 )ψ = 0
with data which are small long-wave perturbations of the constant one, namely ψ = ̺e iφ with (4.13) 5 Actually, to a ground state solution of the KP I equation, that is a minimizer of the Hamiltonian with fixed L 2 norm.

̺ 0 = 1 - ε 2 6 N 0 ε (εx) 1/2 , φ 0 = ε 6 √ 2 Θ 0 ε (εx),
for 0 < ε ≪ 1, and N 0 ε and W 0 ε = ∂ x Θ 0 ε are uniformly bounded in some Sobolev space H k (R) for sufficiently large k.

We refer to [START_REF]On the Korteweg-de Vries long-wave approximation of the Gross-Pitaevskii equation. I[END_REF][START_REF]On the Korteweg-de Vries long-wave approximation of the Gross-Pitaevskii equation II[END_REF] for a detailed analysis and will only summarize the limit to (long) waves propagating in two directions and following a coupled system of KdV equations.

Recall that the one-dimensional Gross-Pitaevskii equation (4.12) is globally well posed in the Zhidkov type spaces

Y k , Y k (R) = ψ ∈ L 1 loc (R; C), 1 -|ψ| 2 ∈ L 2 (R), ∂ x ψ ∈ H k-1 (R) , for any integer k 1.
Moreover the Ginzburg-Landau energy H(ψ(t)) is conserved by the flow and, provided H(ψ 0 ) < 2 √ 2

3 , the corresponding solution ψ(t) does not vanish so that one may write ψ = ̺ exp(iφ), for some continuous function φ.

We consider data as in (4.13) with small enough ε and assume in addition that

N 0 ε M(R) + ∂ x Θ 0 ε M(R) < +∞. Here, • M(R) denotes the norm defined on L 1 loc (R) by f M(R) = sup (a,b)∈R 2 b a f (x)dx .
We next introduce the slow coordinates

x -= ε(x + √ 2t), x + = ε(x - √ 2t), and τ = ε 3 2 √ 2 t.
The definition of the coordinates x -and x + corresponds to reference frames traveling to the left and to the right, respectively, with speed √ 2 in the original coordinates (t, x). We define accordingly the rescaled functions N ± ε and Θ ± ε as follows:

(4.14)

           N ± ε (τ, x ± ) = 6 ε 2 η(t, x) = 6 ε 2 η 2 √ 2τ ε 3 , x ± ε ± 4τ ε 3 with η = 1 -̺ 2 , Θ ± ε (τ, x ± ) = 6 √ 2 ε φ(t, x) = 6 √ 2 ε φ 2 √ 2τ ε 3 , x ± ε ± 4τ ε 3 . Setting (4.15)      U - ε (τ, x -) = 1 2 N - ε (τ, x -) + ∂ x -Θ -(τ, x -) , U + ε (τ, x + ) = 1 2 N + ε (τ, x + ) -∂ x + Θ + (τ, x + ) ,
the main result is (see [START_REF]On the Korteweg-de Vries long-wave approximation of the Gross-Pitaevskii equation II[END_REF] for details):

Theorem 4.4. Let k 0 and ε > 0 be given. Assume that the initial data ψ 0 belongs to Y k+6 (R) and satisfies the assumption

N 0 ε M(R) + ∂ x Θ 0 ε M(R) + N 0 ε H k+5 (R) + ε ∂ k+6 x N 0 ε L 2 (R) + ∂ x Θ 0 ε H k+5 (R) K 0 .
Let U -and U + denote the solutions to the Korteweg-de Vries equations

∂ τ U -+ ∂ 3 x -U -+ U -∂ x -U -= 0, and ∂ τ U + -∂ 3 x + U + -U + ∂ x + U + = 0,
with the same initial value as U - ε and U + ε , respectively. Then, there exist positive constants ε 1 and K 1 , depending only on k and K 0 , such that

U - ε (τ, •) -U -(•, τ ) H k (R) + U + ε (τ, •) -U + (•, τ ) H k (R) K 1 ε 2 exp (K 1 |τ |) , for any τ ∈ R provided ε ε 1 .
We now turn to the two (or higher) dimensional case, which is studied in [START_REF] Chiron | The KdV/KP-I limit of the nonlinear Schrödinger equation[END_REF] for a general nonlinear Schrödinger equation of the form similar to (4.1)

i∂ t ψ + ∆ψ = f (|ψ| 2 )ψ,
where f is smooth and satisfies f (1) = 0, f ′ (1) > 0.

One also uses a "weakly transverse transonic" scaling, namely

T = cε 3 t, X 1 = ε(x 1 -ct), X j = ε 2 x j , j = 2, • • • , d.
After performing the ansatz

ψ ε (t, X) = 1 + ε 2 A ε (t, X) exp(iεφ ε (t, X)),
the hydrodynamic reformulation of the Gross-Pitaevskii equation is used to recast the problem as a singular limit for an hyperbolic system in the spirit of [START_REF] Grenier | Semiclassical limit of the nonlinear Schrödinger equation in small time[END_REF]. Then smooth H s solutions are proven to exist on an interval independent of the small parameter ε. Passing to the limit by a compactness argument yields the convergence of the solutions to that of the KP-I equation. Note however that this method does not provide a convergence rate with respect to ε, contrary to the KdV case considered above.

In comparison, for such data, Theorem 3.3 would ensure that the linear system (3.6) gives a good description of the solution only for times that are o(ε -3 ) (see the introduction of [START_REF] Béthuel | On the linear wave regime of the Gross-Pitaevskii equation[END_REF] for more details).

Global existence of weak solutions to a quantum fluids system

We aim at providing an elementary proof of the result by P. Antonelli and P. Marcati in [START_REF] Antonelli | The quantum hydrodynamics system in two dimensions[END_REF][START_REF]On the finite energy weak solutions to a system in quantum fluid dynamics[END_REF] concerning global finite energy weak solutions to the QHD system. Here is the statement: Theorem 5.1. Let the initial data (ρ 0 , Λ 0 ) ∈ W 1,1 × L 2 be "well-prepared" in the sense that there exists some wave function ψ 0 ∈ H 1 such that

ρ 0 = |ψ 0 | 2 and J 0 := √ ρ 0 Λ 0 = Im( ψ0 ∇ψ 0 ). Assume that f (r) = r σ for some integer σ such that W 1,1 (R d ) ֒→ L σ+1 (R d ).
There exist some vector-field Λ ∈ L ∞ (R; L 2 ) and some nonnegative function

ρ ∈ L ∞ (R; L 1 ∩ L σ+1 ) with ∇ √ ρ ∈ L ∞ (R; L 2
) such that the following system holds true in the distributional sense6 :

(5.1)

                 ∂ t ρ + div J = 0, ∂ t J + div (Λ ⊗ Λ) + ∇(P (ρ)) = 1 4 ∆∇ρ -div (∇ √ ρ ⊗ ∇ √ ρ), ∂ j J k -∂ k J j = 2Λ k ∂ j √ ρ -2Λ j ∂ k √ ρ for all (j, k) ∈ {1, • • • , d} 2 , (ρ, J) |t=0 = (ρ 0 , J 0 ),
that of the velocity equation multiplied by ρ in (1.1) with ε = 1, and the third equation just means that there exists some function φ such that J = ρ∇φ.

with J := √ ρΛ and P (ρ

) := ρf (ρ) -F (ρ) with F (ρ) = ρ 0 f (ρ ′ ) dρ ′ . In addition, the energy R d 1 2 |Λ| 2 + 1 2 |∇ √ ρ| 2 + F (ρ)
is conserved for all time.

Proof. Let us first prove the statement in the smooth case, namely we assume that the data ψ 0 is in H s for some large enough s. It is well known that (5.2)

i∂ t ψ + 1 2 ∆ψ = f |ψ| 2 ψ ; ψ |t=0 = ψ 0
has a unique solution ψ in C(R; H s ) whenever s 1 (see e.g. [START_REF] Cazenave | Semilinear Schrödinger equations[END_REF][START_REF] Linares | Introduction to nonlinear dispersive equations[END_REF][START_REF] Tao | Nonlinear dispersive equations[END_REF]) and that

(5.3) ∀t ∈ R, R d 1 2 |∇ψ(t)| 2 + F |ψ(t)| 2 = R d 1 2 |∇ψ 0 | 2 + F |ψ 0 | 2 .
Let us set ρ := |ψ| 2 and

(5.4)

φ(x) := |ψ(x)| -1 ψ(x) if ψ(x) = 0, 0 if ψ(x) = 0.
We claim that The first term in the right-hand side converges pointwise to Re( φ∇ψ) hence in L 1 loc owing to Lebesgue's theorem as it is bounded by |∇ψ|.

As for the last term, we have

Re(ψ∇ φε ) = ε 2 (ε 2 + |ψ| 2 ) 3/2 Re( ψ∇ψ).
It is clear that the right-hand side converges pointwise to 0 and is bounded by |∇ψ|. Hence it also converges to 0 in L 1 loc . So finally, putting these two results together with (5.6) and (5.7), one may conclude to (5.5).

Let Λ := Im( φ∇ψ) and J := √ ρΛ. We claim that (ρ, Λ, J) satisfies (5.1). Indeed, from (5.2), we see that

∂ t (|ψ| 2 ) = 2 Re( ψ∂ t ψ) = -Im( ψ∆ψ) = -Im div ( ψ∇ψ), hence (5.8) ∂ t ρ + div J = 0.
Next, we compute

∂ j J k -∂ k J j = Im ∂ j ( ψ∂ k ψ) -∂ k ( ψ∂ j ψ) = 2 Im ∂ j ψ∂ k ψ .
Recall (see e.g. Theorem 6.19 in [START_REF] Lieb | Analysis[END_REF]) that (5.9) ∇ψ = 0 a. e. on ψ -1 ({0}) whenever ψ is locally in W 1,1 .

Therefore, given that |φ| = 1 on ψ -1 (C \ {0}), one may write a. e.

Im ∂ j ψ∂ k ψ = Im φ∂ j ψ φ∂ k ψ = Re(φ∂ j ψ) Im( φ∂ k ψ) + Re( φ∂ k ψ) Im(φ∂ j ψ).

So we eventually get (5.10)

∂ j J k -∂ k J j = 2∂ j √ ρ Λ k -2∂ k √ ρ Λ j .
Next, we compute

∂ t J = Im ∂ t ψ ∇ψ + ψ∇∂ t ψ .
Hence, using the equation satisfied by ψ, we get

∂ t J = 1 2 Re ψ∇∆ψ -∆ ψ∇ψ A + Re f (|ψ| 2 ) ψ∇ψ -ψ∇(f (|ψ| 2 )ψ) B .
If ψ is C 1 and the function f has a derivative at every point of R + , then straightforward computations show that (5.11) B = -ρf ′ (ρ)∇ρ = -∇(P (ρ)).

Next, we see that (still in the smooth case)

∇∆|ψ| 2 = 4 Re ∇ 2 ψ : ∇ ψ + 2 Re ∇ψ∆ ψ + 2 Re ψ∇∆ψ . Hence A = 1 2 ∇∆|ψ| 2 -2 Re ∆ ψ∇ψ -2 Re ∇ 2 ψ : ∇ ψ .
So we get (5.12)

1 2 A = 1 4 ∇∆|ψ| 2 -div Re(∇ ψ ⊗ ∇ψ).
Now, using again (5.9), one may write at almost every point of R d , Re(∇ ψ ⊗ ∇ψ) = Re(φ∇ ψ ⊗ φ∇ψ), = Re( φ∇ψ) ⊗ Re( φ∇ψ) + Im( φ∇ψ) ⊗ Im( φ∇ψ).

Therefore, we have

(5.13) Re(∇ ψ ⊗ ∇ψ) = ∇ √ ρ ⊗ ∇ √ ρ + Λ ⊗ Λ.
Putting this together with (5.11) and (5.12), one may conclude that (5.14)

∂ t J + div (Λ ⊗ Λ) + ∇(P (ρ)) = 1 4 ∆∇ρ -div (∇ √ ρ ⊗ ∇ √ ρ).
Of course, as owing to (5.9)

(5.15)

|∇ψ| 2 = (Re( φ∇ψ)) 2 + (Im( φ∇ψ)) 2 = |∇ √ ρ| 2 + |Λ| 2 a. e.,
the energy equality for ψ recasts in (5.16)

R d 1 2 |Λ(t)| 2 + 1 2 |∇ ρ(t)| 2 +F (ρ(t)) = R d 1 2 |Λ 0 | 2 + 1 2 |∇ √ ρ 0 | 2 +F (ρ 0 ) .
This completes the proof in the smooth case.

Let us now treat the rough case where ψ 0 belongs only to H 1 . Then we fix some sequence (ψ 0,n ) n∈N of functions in H s (with s large) converging to ψ 0 in H 1 . Let us denote by ψ n the solution of (5.2) in C(R; H s ) corresponding to the data ψ 0,n , and by ψ ∈ C(R; H 1 ) the solution to (5.2) with data ψ 0 .

From the first part of the proof, we know that there exists some sequence (φ n ) n∈N of functions with modulus at most 1 such that if we set Λ n := Im( φn ∇ψ n ),

ρ n := |ψ n | 2 and J n := √ ρ n Λ n then ∇ √ ρ n = Re( φn ∇ψ n ) in L 2 ,
and (ρ n , Λ n , J n ) satisfies (5.1) with data (ρ 0,n , J 0,n ), together with the energy equality

R d 1 2 |Λ n (t)| 2 + 1 2 |∇ ρ n (t)| 2 + F (ρ n (t)) (5.17) = R d 1 2 |Λ 0,n | 2 + 1 2 |∇ √ ρ 0,n | 2 + F (ρ 0,n ) .
We now have to prove the convergence of (ρ n , Λ n , J n ) to some solution (ρ, Λ, J) of (5.1) satisfying the energy equality. On the one hand, standard stability estimates (based on Strichartz inequalities) guarantee that (5.18) ψ n -→ ψ in L ∞ loc (R; H 1 ). On the other hand, because (φ n ) n∈N is bounded by 1, it converges (up to extraction) in L ∞ weak * to some function φ such that φ L ∞ 1. As ∇ψ n converges strongly to ∇ψ in L 2 , this implies that (5.19) φn ∇ψ n ⇀ φ∇ψ in L 2 .

In turn, as obviously √ ρ n → √ ρ in L 2 and as ∇ √ ρ n = Re( φn ∇ψ n ), we deduce that ∇ √ ρ = Re( φ∇ψ) and ∇ √ ρ n ⇀ ∇ √ ρ in L 2 . Given that Λ n = Im( φn ∇ψ n ), (5.19) also ensures that Λ n ⇀ Λ := Im( φ∇ψ) in L 2 .

In order to establish that strong convergence in L 2 holds true, it suffices to show that (up to an omitted extraction)

(5.20) ∇ √ ρ n L 2 → ∇ √ ρ L 2 and Λ n L 2 → Λ L 2 .
On the one hand, the weak convergence ensures that

∇ √ ρ L 2 lim inf ∇ √ ρ n L 2 and Λ L 2 lim inf Λ n L 2 ;
on the other hand, given that (∇ψ n ) n∈N converges strongly to ∇ψ in L 2 and that (5.15) holds true for ψ and ψ n , we may write

∇ √ ρ 2 L 2 + Λ 2 L 2 = ∇ψ 2 L 2 , = lim n→+∞ ∇ψ n 2 L 2 , = lim n→+∞ ∇ √ ρ n 2 L 2 + Λ n 2 L 2 , lim inf ∇ √ ρ n L 2 2 + lim inf Λ n L 2 2 .
Therefore (5.20) is satisfied. As a conclusion, we thus have established that

√ ρ n → √ ρ in H 1 and Λ n → Λ in L 2 .
Of course, this implies that J n → J := √ ρΛ in L 1 so it is easy to pass to the limit in (5.1) and in the energy equality (5.16). The details are left to the reader.

Momentum of inertia or virial.

A direct computation shows that the momentum of inertia (or virial) (A.9) This equality can be proven more simply by using the fact that the operators iε∂ t + ε 2 2 ∆ and x + iεt∇ commute. The conservation law for Z is not independent from the preceding ones: by expanding the square of the modulus, one observes that Z(t) = t 2 H -tF + I. To prove (A.6), one writes (using the Einstein summation convention),

I

Consequently

A i = ε ijk x j ρv k
with ε ijk = 0 if two of the indices are equal, and equal to signature of (i j k) otherwise.

Using the conservative form of the equation for ρv and integrating by parts, one thus deduces

d dt A i = ε ijk x j ∂ ℓ K ℓk dx -ε ijk x j ∂ k P -ε ijk x j ∂ ℓ (ρv k v ℓ ) dx, = -ε ijk K jk dx + 0 -ε ijk ρv j v k dx.
The tensors K and ρv ⊗ v being symmetric, the right-hand member of the previous equality vanishes.

The simplest way to prove the conservation law on F , is to use the second equation in (A.12). Integrating by parts the capillary term and treating the other terms as in the compressible Euler equation one obtains

dF dt = x • div K dx -x • ∇P + div (ρv ⊗ v) dx, = x j ∂ i K ij dx + (dP + ρ|v| 2 ) dx.
Using the expression of the tensor K and integrating by parts, one gets

x j ∂ i K ij dx = -tr K dx, = -d 1 2 (κ + ρκ ′ )|∇ρ| 2 + ρκ∆ρ dx + κ|∇ρ| 2 , = d 2 (ρκ) ′ |∇ρ| 2 dx + κ|∇ρ| 2 dx,
from which the claimed equality results.

Remark A.2. The case of QHD corresponds to (ρκ) ′ = 0. For capillary fluids, the presence of extra terms in the conservation law for F seems relatively harmless if ρ → ρκ(ρ) is an increasing function. One should be able to prove without difficulty that with a pressure law such that dP -2F 0, one has

I(t) ∼ t→+∞ Ht 2 .
In the case where ρ → ρκ is decreasing, we expect the conservation law for F to be the key to proving finite time blow-up results similar to those of [START_REF] Sideris | Formation of Singularities in Three-Dimensional Compressible Fluids[END_REF].

a

  j b j for a and b in C d .

(4. 11 )

 11 -w xx -w yy + w xxxx + (ww x ) x = 0. More precisely, denoting N p = 1 -|u p | 2 where u p = |u p |e iφp is a minimizer of the energy with fixed momentum p, and

5

 5 Re( φ∇ψ) a. e. Indeed, for ε > 0, let us set φ ε := ψ/ |ψ| 2 + ε 2 . Then (φ ε ) converges pointwise to φ and an easy computation shows that (φε ψ) = Re( φε ∇ψ) + Re(ψ∇ φε ).

  Pseudo-conformal invariance. Let ϕ the pseudo-conformal transform of ψ defined by + iεt∇)ψ|2 + t 2 F (|ψ| 2 ) dx, |∇ √ ρ| 2 + t 2 F (|ψ| 2 ) dx.

					:=	1 2		|x| 2 |ψ| 2 dx =	1 2	|x| 2 ρ dx
	satisfies	dI dt	= F .
					ϕ(t, x) =	e i t i|x| 2 2εt ε d/2	ψ ε 2 t	,	x t	•
	One notices that
			iε∂ t ϕ +	ε 2 2	∆ϕ (t, x) =	ε 2 t 2	e (i t i|x| 2 2εt ε ) n 2		iε∂ t ψ +	ε 2 2	∆ψ	ε 2 t	,	x t	.
	Thus
				iε∂ t ϕ +	ε 2 2	∆ϕ =	ε 2 t 2 f			t ε
	For (A.1), this yields
	(v, ρ)(t, x) solution ⇐⇒	x t	-	v t	ε 2 t	,	x t	,	ε t	d	ρ	ε 2 t	,	x t	solution.
	Using the conservation of energy for ϕ, one deduces after a lengthy computation
	that
	(A.10)				dZ dt	+ t (dP -2F )(|ψ| 2 ) dx = 0
	with
			Z(t) := |(x = 1 2 1 2 ρ|x -tv| 2 +	1 2	ε 2 t 2

d |ϕ| 2 ϕ.

For the L 2 -critical power nonlinearity f (r) = r 2/d , one checks that ψ(t, x) solution of (A.2) ⇐⇒ ϕ(t, x) solution of (A.2).

We could consider more general defocusing nonlinearities such as f |ψ ε |

ψ ε = |ψ ε | 2σ ψ ε , σ ∈ N.We consider the exactly cubic case for the simplicity of the exposition only.

Note that this implies that there exists c > 0 such that ρ 0 (x) c for all x ∈ R d .

The method may seem uselessly complicated. However, the algebraic cancellations that are going to be used remain the same when estimating higher order Sobolev norms.

In this section, we use the normalization of[START_REF] Béthuel | Existence and properties of travelling waves for the Gross-Pitaevskii equation, Stationary and time dependent Gross-Pitaevskii equations[END_REF][START_REF]On the KP I transonic limit of two-dimensional Gross-Pitaevskii travelling waves[END_REF], instead of ψ = √ ρe iφ .

In the smooth non-vanishing case, the right-hand side of the second equation coincides with

This work was supported by the French ANR projects Equa-disp (ANR-07-BLAN-0250-01) and R.A.S. (ANR-08-JCJC-0124-01).

Appendix A. Conservation laws

In this Appendix we review conservations laws for the nonlinear Schrödinger, QHD and compressible Euler equations. Even though most of the results are classical (as concerns the Schrödinger equation, they may be found in the textbooks [START_REF] Sulem | The nonlinear Schrödinger equation, self-focusing and wave collapse[END_REF][START_REF] Tao | Nonlinear dispersive equations[END_REF] for instance; links between Schrödinger and Euler conservation laws may be found in [START_REF] Carles | Semi-classical analysis for nonlinear Schrödinger equations[END_REF]), we believe the relationships between the aforementioned equations to be of interest. In addition, those conservation laws are still meaningful for the less classical framework of general Korteweg fluids.

A.1. The case of Schrödinger, QHD and compressible Euler equations. For the time being, we consider the following system (A.1)

which is the QHD system if ε > 0, and the compressible Euler equation if ε = 0, and the nonlinear Schrödinger equation:

Recall that for ε > 0, one may pass formally from (A.2) to (A.1) by setting

In what follows, the function f is assumed to be continuous on R + and, say, C 1 on (0, +∞[, standard cases being f (r) = r σ and f (r) = r -1. We denote by F the anti-derivative of f which vanishes at 0, and set P (ρ) = ρf (ρ) -F (ρ). As pointed out before, from a physical viewpoint, P is the pressure.

The first part of the appendix aims at listing (and deriving formally) the classical conservation laws for (A.2) and (A.1).

Phase invariance. For every α ∈ R, one has ψ solution of (A.2) ⇐⇒ e iα ψ solution of (A.2).

The phase invariance is not seen at the level of (A.1) (this amount to changing φ into φ + εα).

By Noether's theorem or by an easy computation, this leads to the conservation of mass:

Time translation invariance. For every τ ∈ R, one has ψ(t, x) solution of (A.2) ⇐⇒ ψ(t + τ, x) solution of (A.2).

By the same time translation, this is expressed at the level of (A.1) and leads to the conservation of the energy (or Hamiltonian):

Space translation invariance. For every x 0 ∈ R d , one has

By the same space translation, this is expressed at the level of (A.1) and leads to the conservation of momentum:

Invariance by spatial rotation. Let R be a spatial rotation. Then ψ(t, x) solution of (A.2) ⇐⇒ ψ(t, Rx) solution of (A.2).

By the same spatial rotation on a solution of (A.1) this leads to the conservation of angular momentum, which we write in the case of R 3 for the sake of simplicity:

Galilean invariance. For every ξ 0 ∈ R d , one has ψ(t, x) solution of (A.2) ⇐⇒ e -iξ0•x e -i tε 2 |ξ0| 2 ψ(t, x + εξ 0 t) solution of (A.2). For (A.1), this implies (v, ρ)(t, x) solution ⇐⇒ v(t, x + εξ 0 t) -εξ 0 , ρ(t, x + εξ 0 t) solution, and leads to (A.7) dX dt = P with X := x|ψ| 2 dx = xρ dx.

Scale invariance. When f (r) = r σ , one has for all λ > 0, ψ(t, x) solution of (A.2) ⇐⇒ λ 1/σ ψ(λ 2 t, λx) solution of (A.2), which implies for (A.1)

The associated conservation law is (recall that d stands for the space dimension):

Equality (A.8) remains formally true if the nonlinearity is not a pure power.

The Carles-Nakamura conservation law. Formally, the quantity

is constant. This can be checked by an indirect way (see [START_REF] Carles | Nonlinear Schrödinger equations with Stark potential[END_REF]) or by a direct computation. The interpretation in terms of (A.1) is

This conservation law also results from the conservation of momentum and from the conservation law for X since in both cases, (A.2) or (A.1), one has

A.2. Conservation laws for general capillary fluids. We here study to what extent the conservation laws listed in the previous subsection are relevant for general inviscid capillary fluids. We recall that such fluids are governed by System (3.7). Throughout, we assume the capillarity κ to be a differentiable function on R + .

For smooth solutions with a non vanishing density, System (3.7) recasts in the following conservative form 7 : (A.12)

Theorem A.1. If (ρ, v) is a sufficiently smooth solution of (3.7) which decays at infinity, with ρ non vanishing, then the equalities (A.3), (A.5), (A.6), (A.7) and (A.9) are still valid, and the energy

Furthermore, the quantity F defined in (A.8) satisfies

Proof. It is very likely that most of those conservation laws could be derived from Noether theorem. They can also been obtained by pedestrian computations. For (A.3), (A.7) and (A.9), there is no difference with the QHD since the only equation on ρ is concerned.

For (A.5) the proof is obvious from the conservative form (A.12). In order to prove the conservation of energy, the simplest method is to take the L 2 scalar product with ρv the first equation of (3.7). The treatment of the terms on the left-hand side of the equation for v is the same as in the compressible Euler 7 With the convention (div K) j := i ∂ i K ij