<ph f="cmr"> </ph><ph f="cmbx">The General Quasilinear Ultrahyperbolic Schrödinger Equation</ph>

C. E. Kenig, G. Ponce, C. Rolvung, and L. Vega

1 Introduction

In this article we consider nonlinear Schrödinger equations of the form
∂ t u = − i ∂ x j ( a j k ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ∂ x k u ) + b ⃗ 1 ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ⋅ ∇ u + b ⃗ 2 ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ⋅ ∇ u ¯ + c 1 ( x , t , u , u ¯ ) u + c 2 ( x , t , u , u ¯ ) u ¯ + f ( x , t ) , (1.1)
where x ∈ R n   , t > 0   , and A = ( a j k ( ⋅ ) ) j , k = 1 , . . , n   is a real, symmetric matrix.
Our aim is to study the existence, uniqueness and regularity of local solutions to the initial value problem (IVP) associated to the equation (1.1).
In the case where A = ( a j k ( ⋅ ) ) j , k = 1 , . . , n   is assumed to be elliptic the local solvability of the IVP associated to (1.1) was recently established in [20] . Hence, in this work we should be concerned with the case where ( a j k ( ⋅ ) ) j , k = 1 , . . , n   is just a non-degenerate matrix.
Equations of the form described in (1.1) with A = ( a j k ( ⋅ ) ) j , k = 1 , . . , n   merely invertible arise in water wave problems, and as higher dimensions completely integrable models, see for example [1] , [7] , [8] , [15] , [27] , and [30] .
There are significant differences in the arguments required for the local solvability in the case where A   is a non-degenerate matrix in comparison with the elliptic case treated in [20] . To illustrate them as well as to review some of the previous related results we consider first the semi-linear equation with constant coefficients (for more details and further references and comments see [19] , [20] , [21] , and references therein)
∂ t u = − i ( ∂ x 1 2 + . . + ∂ x k 2 − ∂ x k + 1 2 − . . − ∂ x n 2 ) u + P ( u , u ¯ , ∇ x u , ∇ x u ¯ ) , (1.2)
where P ( ⋅ )   is the non-linearity (by simplicity a polynomial in its variables without constant or linear terms).
In [18] based on the smoothing effects (homogeneous and inhomogeneous, see [16] , [22] , [4] , [25] , [29] , [17] ) associated to the group { e − i t ( ∂ x 1 2 + . . + ∂ x k 2 − ∂ x k + 1 2 − . . − ∂ x n 2 ) : t ∈ R }   the local wellposedness for “small” data for the IVP associated to  1.2 was deduced. In [12] for the one dimensional case ( n = 1   ), Hayashi and Ozawa eliminated the size restriction on the data in [18] . Their argument was based on a change of variable which transforms the equation into a new system without the term ∂ x u   , so that the standard energy estimate yields the desired local result. In [5] Chihara, for the elliptic case (i.e. k = n   in  1.2 ), removed the size restriction on the data in any dimension n   . Roughly speaking, the argument there first uses the ellipticity to diagonalize the system for ( u , u ¯ )   , and then introduced an operator K   so that the commutator i [ K ; Δ ]   “controls” the term K b ⃗ ( x ) ⋅ ∇ x   . This is achieved by combining some result of Doi [9] concerning the local smoothing effects in the solution with the sharp version of Garding inequality.
If instead of “controlling” it one asks for the operator K   to verify that
− i [ K ; Δ ] + K b ⃗ ( x ) ⋅ ∇ x = 0 + order zero , (1.3)
one finds that K   has symbol
k ( x , ξ ) = e x p ( − ∫ 0 ∞ b ⃗ ( x + 2 s ξ ) ⋅ ξ d s ) , (1.4)
which is in the non-standard class studied by Craig-Kappeler-Strauss [6] . In particular, it satisfies that | ∂ x α ∂ ξ β k ( x , ξ ) | ≤ c α β 〈 x 〉 | β | 〈 ξ 〉 − | β | , α , β ∈ ( Z + ) n .   However, in the non-elliptic case with coefficients depending just on the space variable x   , the geometric assumptions in [6] (Chapter 3.1, section 3.1) does not hold for the relevant symbols.
The local wellposedness of the IVP associated to the equation in  1.2 ( 1 < k ≤ n   ) was established in [19] . The method of proof there, among other arguments, utilizes the symbol class S 0 , 0 0   of Calderón-Vaillancourt [3] . However, this approach does not seem to extend to the variable coefficients case.
Next, we consider the linear IVP
∂ t u = − i ∂ x k a j k ( x ) ∂ x j u + b ⃗ 1 ( x ) ⋅ ∇ u + b ⃗ 2 ( x ) ⋅ ∇ u ¯ + f ( x , t ) , u ( x , 0 ) = u 0 ( x ) . (1.5)
We recall the notion of the bicharacteristic flow associated to the symbol of the principal part of the operator − ∂ x k a j k ( x ) ∂ x j   , i.e. h ( x , ξ ) = a j k ( x ) ξ k ξ j   .
Let ( X ( s ; x 0 , ξ 0 ) , Ξ ( s ; x 0 , ξ 0 ) )   denote the solution of the Hamiltonian system
d d s X j ( s ; x 0 , ξ 0 ) = 2 a j k ( X ( s ; x 0 , ξ 0 ) ) Ξ k ( s ; x 0 , ξ 0 ) = ∂ ξ j h , d d s Ξ j ( s ; x 0 , ξ 0 ) = − ∂ x j a l k ( X ( s ; x 0 , ξ 0 ) ) Ξ k Ξ l ( s ; x 0 , ξ 0 ) = − ∂ x j h , (1.6)
for j = 1 , . . , n   , with data ( X ( 0 ; x 0 , ξ 0 ) , Ξ ( 0 ; x 0 , ξ 0 ) ) = ( x 0 , ξ 0 ) .   Under mild regularity assumptions on the coefficients a j k ( x )   's the bicharacteristic flow ( X ( s ; x 0 , ξ 0 ) , Ξ ( s ; x 0 , ξ 0 ) )   is defined in the interval s ∈ ( − δ , δ )   with δ = δ ( x 0 , ξ 0 ) > 0   depending continuously on ( x 0 , ξ 0 )   . If the operator − a j k ( x ) ∂ x k x j 2   is elliptic, i.e. ( a j k ( x ) )   is positive definite, using that the flow preserves h   , i.e.
H h h = ∑ j = 1 n ( ∂ ξ j h ∂ x j p − ∂ x j h ∂ ξ j p ) = { h , p } = 0 ,   one has that ν − 2 | ξ 0 | 2 ≤ | Ξ ( s ; x 0 , ξ 0 ) | 2 ≤ ν 2 | ξ 0 | 2 .   Hence δ = ∞   , i.e. the bicharacteristic flow is globally defined. In the non-elliptic case one needs to prove it.
In [14] Ichinose established that for the L 2   -local wellposedness of IVP  1.5 , with − ∂ x k a j k ( x ) ∂ x j   elliptic and b 2 ( x ) ≡ 0   , it is necessary that
sup x ∈ R n , ω ∈ S n − 1 , R > 0 | I m ∫ 0 R b ⃗ 1 ( X ( s ; x , ω ) ) ⋅ Ξ ( s ; x , ω ) d s | < ∞ . (1.7)
This extends previous results of Mizohata [23] and Takeuchi [28] for the constant coefficient case, where ( X ( s ; x 0 , ξ 0 ) , Ξ ( s ; x 0 , ξ 0 ) ) = ( x 0 + 2 s ξ 0 , ξ 0 )   . Notice that in this variable coefficient case the “integrating factor” in  1.4 reads
k ( x , ξ ) = e x p ( − ∫ 0 ∞ b ⃗ 1 ( X ( s ; x , ξ ) ) ⋅ Ξ ( s ; x , ξ ) d s ) . (1.8)
The condition  1.7 justifies the following “non-trapping” assumption. The bicharacteristic flow  1.6 is non-trapped if the set { X ( s ; x 0 , ξ 0 ) : s ∈ R + }   is unbounded in R n   for each ( x 0 , ξ 0 ) ∈ R n × R n − { 0 }   .
As it was already mentioned the IVP for the equation  1.1 with ( a j k ( ⋅ ) ) j , k = 1 , . . , n   elliptic was studied in [20] . There the local solvability was obtained under regularity and decay assumptions on the coefficients. Also a non-trapping character of the bicharacteristic flow associated to the principal symbol of the elliptic operator − ∂ x j ( a j k ( x , t , u , u ¯ , ∇ x u , ∇ x u ¯ ) ∂ x k   when evaluated at the data u ( x , 0 ) = u 0 ( x )   and at a time t = 0   was assumed, i.e. the bicharacteristic flow associated to the symbol
h ( u 0 ) = h u 0 ( x , ξ ) = a j k ( x , 0 , u 0 , u ¯ 0 , ∇ x u 0 , ∇ x u ¯ 0 ) ξ j ξ k . (1.9)
The proofs of the semi-linear results in [18] , [19] follow a fixed point theorem (via contraction principle). This approach does not extend to the quasi-linear case, thus in [20] the proof used the so called “artificial viscosity” method.
We recall that one of the advantages of the contraction principle approach is that it also provides the continuity (in fact, the analyticity) of the solution upon the data.
Returning to the non-degenerate case, in [21] the following semi-linear IVP was studied
∂ t u = − i ∂ x k a j k ( x ) ∂ x j u + b ⃗ 1 ( x ) ⋅ ∇ u + b ⃗ 2 ( x ) ⋅ ∇ u ¯ + c 1 ( x ) u + c 2 ( x ) u ¯ + P ( u , ∇ u , u ¯ , ∇ u ¯ ) , u ( x , 0 ) = u 0 ( x ) , (1.10)
where the non-linearity P   is given by a polynomial without linear or constant terms.
Under the following assumptions:
(a) (Non-degeneracy) There exists ν ∈ ( 0 , 1 )   ν | ξ | ≤ | A ( x ) ξ | ≤ ν − 1 | ξ | , x , ξ ∈ R n ,   where A ( x ) = ( a j k ( x ) ) j , k = 1 , . . , n   .
(b) (Asymptotic Flatness) There exist c 0 > 0   and N , M ~ ∈ Z +   large enough such that
| ∂ x α ( a j k ( x ) − a j k 0 ) | ≤ c 0 〈 x 〉 N , | α | ≤ M ~ , j , k = 1 , . . , n , (1.11)
where
A h = ( a j k 0 ) j , k = 1 , . . , n = ( I k × k 0 0 − I ( n − k ) × ( n − k ) ) . (1.12)
(c) (Non-trapping condition) The initial data u 0   satisfies that the bicharacteristic flow associated to  1.9 is non-trapping.
(d) (Growth condition of the first order coefficients) There exist c 0 > 0   and N , M ~ ∈ Z +   large enough such that
| ∂ x α b j k ( x ) | ≤ c 0 〈 x 〉 N , | α | ≤ M ~ , j = 1 , 2 , k = 1 , . . , n . (1.13)
(e) Regularity of the coefficients For J ∈ Z +   sufficiently large a j k , b m j , c m ∈ C b J ( R n : C ) , j , k = 1 , . . , n , m = 1 , 2 .   the following result was established in [21] .
Theorem 1.1. [21] There exist s , N > 0   , s > N   , and N ~ > 0   , depending only on n   , such that for u 0 ∈ H s ( R n ) ∩ L 2 ( 〈 x 〉 N d x )   there exists T = T ( ∥ u 0 ∥ s , 2 , ∥ 〈 x 〉 N / 2 u 0 ∥ 2 ) > 0   such that the IVP  1.10 has a unique solution u   defined in the time interval [ 0 , T ]   satisfying u ∈ C ( [ 0 , T ] : H s ( R n ) ∩ L 2 ( 〈 x 〉 N d x ) ) ≡ X T s , N .   Moreover, ∫ 0 T ∫ | J s + 1 / 2 u ( x , t ) | 2 〈 x 〉 − N ~ d x d t < ∞ ,   and for every u 0 ∈ H s ( R n ) ∩ L 2 ( 〈 x 〉 N d x )   there exists a neighborhood U   of u 0   and a T ′ > 0   such that the map data →   solution of is continuous from U   into X T ′ s , N   .
Here ∥ ⋅ ∥ s , 2   denotes the norm in the Sobolev spaces H s ( R n )   .
As in the previous semi-linear cases the proof in [21] was based on the contraction principle.
In this non-degenerate case the operator describing the “integrating factor” in  1.8 has not been shown to be an L 2   -bounded operator. However, thanks to the local smoothing effects it suffices to solve  1.3 up to “small” first order term. This is achieved in [21] by introducing a new class of symbols.
Definition 1.2. (i) It will be said that a ∈ S ( R n : S 1 , 0 m )   (where S 1 , 0 m   the class of ψ   .d.o's of classical symbols of order m   ) if a ∈ C ∞ ( R n × R n × R n )   and a   satisfies
| 〈 z 〉 μ ∂ s α ∂ x β ∂ ξ γ a ( s ; x , ξ ) | ≤ c μ α β γ 〈 ξ 〉 m − | γ | , ∀ z , x , ξ ∈ ( Z + ) n , ∀ μ , α , β , γ ∈ ( Z + ) n . (1.14)
(ii) For a ∈ S ( R n : S 1 , 0 m )   let
b ( x , ξ ) = χ ( | ξ | ) a ( P ( x , A h ξ ) ; x , ξ ) , (1.15)
where P ( y , z ) = y − ( y ⋅ z ) z / | z | 2   for y , z ∈ R n , z ≠ 0   , is the projection of y   onto the hyperplane perpendicular to z   , A h   as in  1.12 , and χ ∈ C ∞ ( R n )   with χ ( l ) = 0   for | l | < 1   and χ ( l ) = 1   for | l | > 2   .
In fact, we showed in [21] that it suffices to have  1.14 for sufficiently large μ , α , β , γ ∈ ( Z ) n   .
We observe that if a ∈ S ( R n : S 1 , 0 m )   , then ∂ ξ α ( ξ β a ( ⋅ ) ) ∈ S ( R n : S 1 , 0 k ) , k = m + | β | − | α | ,   i.e. this class is closed with respect to differentiation and multiplication in the ξ   -variable.
We shall show below that this class is also closed with respect to differentiation in the x   -variable.
Roughly speaking, for r ∈ Z +   large enough 〈 x 〉 − r χ ( | ξ | ) a ( P ( x , A h ξ ) ; x , ξ ) = 〈 x 〉 − r b ( x , ξ ) ,   is a symbol in the class S 1 , 0 m   .
In [21] , we deduce several properties of operators with symbol in our class. These include their continuity from H m ( R n )   to L 2   and their composition rules with classical differential operators P ( x , ∂ x )   with decaying coefficients, i.e. with P ( x , ∂ x ) = φ α ( x ) ∂ x α   with | x | l φ α   bounded for l ∈ Z   sufficiently large.
The proof of the nonlinear results concerning the IVP  1.10 relies on two key linear estimates. The first one is concerned with the smoothing effect described for solutions of the IVP  1.5 with ( a j k ( x ) )   being just an invertible matrix
∫ 0 T ∫ R n | J s + 1 / 2 u ( x , t ) | 2 〈 x 〉 − N ~ d x d t ≤ c ( 1 + T ) sup 0 ≤ t ≤ T ∥ u ( t ) ∥ s , 2 2 + c ∫ 0 T ∫ R n | f ( x , t ) | 2 d x d t , (1.16)
or
∫ 0 T ∫ R n | J s + 1 / 2 u ( x , t ) | 2 〈 x 〉 − N ~ d x d t ≤ c ( 1 + T ) sup 0 ≤ t ≤ T ∥ u ( t ) ∥ s , 2 2 + ∫ 0 T ∫ R n | J s − 1 / 2 f ( x , t ) | 2 〈 x 〉 N ~ d x d t (1.17)
for N ~ > 1   .
The second is related with the local wellposedness in L 2   (and in H s   ) of the IVP  1.5 . To establish it we followed an indirect approach. First, we truncated at infinity the operator ℒ ( x ) = − ∂ x k a j k ( x ) ∂ x j   using θ ∈ C 0 ∞ ( R n )   with θ ( x ) = 1 , | x | ≤ 1   , and θ ( x ) = 0 , | x | ≥ 2   .
For R > 0   we define
ℒ R ( x ) = θ ( x / R ) ℒ ( x ) + ( 1 − θ ( x / R ) ) ℒ 0 , (1.18)
where ℒ 0 = − a j k 0 ∂ x j x k 2   , A h = ( a j k 0 ) j , k = 1 , . . , n   as in  1.12 , with the decay assumption a j k ( x ) − a j k 0 ∈ S ( R n )   , j , k = 1 , . . , n   , (the same proof worked if one just assumed that the corresponding estimate held for a sufficiently large number of semi-norms of S ( R n )   ).
Thus, ℒ ( x ) = ℒ R ( x ) + ℰ R ( x ) .   For R   large enough we considered the bicharacteristic flow ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) )   associated to the operator ℒ R ( x )   and the corresponding integrating factor K R   , i.e. the operator with symbol as in  1.8 but evaluated in the bicharacteristic flow ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) )   .
We deduced several estimates concerning the operator K R   . In particular, we showed in [21] that there exists N 0 ∈ Z +   (depending only on the dimension n   ) such that for any M ∈ Z +   there exist N 1 = N 1 ( M ) ∈ Z +   and R 0   sufficiently large such that for R ≥ R 0   it follows that
sup 0 ≤ t ≤ T ∥ K R u ( t ) ∥ 2 2 ≤ c 0 R N 0 ∥ u ( 0 ) ∥ 2 2 + R − M ∫ 0 T ∫ R n | J 1 / 2 u | 2 〈 x 〉 − N ~ d x d t + c 0 T R N 0 + N 1 ( M ) sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 , (1.19)
for N ~   large enough depending only on the dimension n   . In [21] the coefficients in  1.10 were taken in the Schwartz class S ( R n )   . However, it is clear from the proof there that the same argument works if one just assumes that a fixed large number of semi-norms of S ( R n )   of the coefficients are bounded. In this case, M   will be chosen depending only on the decay of the coefficients. More precisely, one chooses M = M ( N )   , M ( N ) ↑ ∞   as N ↑ ∞   , with N   as in  1.11 ,  1.13 .
To complete the estimate one needs to consider the operator E R = I − K ~ R ( K R ) *   ,where the symbol of K ~ R   differs from that of K R   only in the sign of the exponent, and ( K R ) *   is the adjoint of K R   . It was established that E R u ( t )   satisfies an estimate similar to that in  1.16 . Combining these results we get that
sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 ≤ c 0 R N 0 ∥ u ( 0 ) ∥ 2 2 + R − M ∫ 0 T ∫ R n | J 1 / 2 u | 2 〈 x 〉 − N ~ d x d t + c 0 T R N 0 + N 1 ( M ) sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 + c 0 R N 0 ∫ 0 T ∫ R n | f ( x , t ) | 2 d x d t . (1.20)
From  1.16 and  1.20 fixing T   small enough one gets that
sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 ≤ c 0 R N 0 ∥ u ( 0 ) ∥ 2 2 + c 0 R N 0 ∫ 0 T ∫ R n | f ( x , t ) | 2 d x d t . (1.21)
This allows to use the contraction principle to obtain in [21] the desired nonlinear result.
Returning to the IVP  1.1 we shall assume that the coefficients satisfy the following hypotheses:
(H1) Non-degeneracy.Given r 0 > 0   there exists γ r 0 ∈ ( 0 , 1 )   such that for any
( x , t , z ⃗ ) ∈ R n × R × ( ( B r 0 ( 0 ) ) ¯ ) 2 n + 2 ≡ D r 0 , γ r 0 | ξ | ≤ | a j k ( x , t , z ⃗ ) ξ | ≤ γ r 0 − 1 | ξ | , ∀ ξ ∈ R n , (1.22)
where B r 0 ( 0 ) ¯ = { z ∈ C : | z | ≤ r 0 }   , and with
A ( x , 0 , 0 ⃗ ) = ( a j k ( x , 0 , 0 ⃗ ) ) j , k = 1 , . . , n = A 0 ( x ) + A h = ( a 0 , j k ( x ) ) j , k = 1 , . . , n + ( a j k 0 ) j , k = 1 , . . , n (1.23)
where for some N , M ~ ∈ Z +   large enough
| ∂ x α a 0 , j k ( x ) | ≤ c 0 〈 x 〉 N , | α | ≤ M ~ j , k = 1 , . . , n , (1.24)
and A h   as in  1.12 .
(H2) Asymptotic flatness.There exists c > 0   such that for any ( x , t ) ∈ R n × R   and 0 < | α | ≤ M ~   , 0 ≤ | α ′ | ≤ M ~   it follows that
| ∂ x α a j k ( x , t , 0 ⃗ ) | + | ∂ t ∂ x α ′ a j k ( x , t , 0 ⃗ ) | ≤ c 〈 x 〉 N , (1.25)
for k , j = 1 , . . , n   with N , M ~   as in  1.24 .
(H3) Growth of the first order coefficients.There exist c , c 0 > 0   such that for any x ∈ R n   and any ( x , t ) ∈ R n × R  
| ∂ x α b m j ( x , 0 , 0 ⃗ ) | ≤ c 0 〈 x 〉 N , | α | ≤ M ~ , m = 1 , 2 , j = 1 , . . , n , | ∂ t ∂ x α b m j ( x , t , 0 ⃗ ) | ≤ c 〈 x 〉 N , | α | ≤ M ~ , m = 1 , 2 , j = 1 , . . , n , (1.26)
where b ⃗ l = ( b l 1 , . . , b l n ) , l = 1 , 2   with N , M ~   as in  1.24 .
(H4) Regularity.For any J ∈ Z +   and r 0 > 0   the coefficients
a j k , b 1 j , b 2 j ∈ C b J ( R n × R × ( ( B r 0 ( 0 ) ) ¯ ) 2 n + 2 ) , (1.27)
and
c 1 , c 2 ∈ C b J ( R n × R × ( ( B r 0 ( 0 ) ) ¯ ) 2 ) . (1.28)
Our main result is the following.
Theorem 1.3. Under the hypotheses (H1)-(H4) there exists N = N ( n ) ∈ Z +   such that given any
u 0 ∈ H s ( R n ) with 〈 x 〉 N ∂ x α u 0 ∈ L 2 ( R n ) , | α | ≤ s 1 , (1.29)
and
f ∈ L 1 ( R + : H s ( R n ) ) with 〈 x 〉 N ∂ x α f ∈ L 1 ( R + : L 2 ( R n ) ) , | α | ≤ s 1 , (1.30)
where s , s 1 ∈ Z +   , sufficiently large with s > s 1 + 4   , for which the Hamiltonian flow H h ( u 0 )   associated to the symbol
h ( u 0 ) = h u 0 ( x , ξ ) = ∑ j , k = 1 n a j k ( x , 0 , u 0 , u ¯ 0 , ∇ u 0 , ∇ u ¯ 0 ) ξ j ξ k , (1.31)
is non-trapping, there exist T 0 > 0   , depending on
λ = ∥ u 0 ∥ s , 2 + ∑ | α 1 | ≤ s 1 ∥ 〈 x 〉 N ∂ x α u 0 ∥ 2
+ ∫ 0 ∞ ∥ f ( t ) ∥ s , 2 d t + ∑ | α 1 | ≤ s 1 ∫ 0 ∞ ∥ 〈 x 〉 N ∂ x α f ( t ) ∥ 2 d t ,
the constants in (H1)-(H4) and on the non-trapping condition (H5), and a unique solution u = u ( x , t )   of the equation in  1.1 with initial data u ( x , 0 ) = u 0 ( x )   on the time interval [ 0 , T 0 ]   satisfying
u ∈ C ( [ 0 , T 0 ] : H s − 1 ) ∩ L ∞ ( [ 0 , T 0 ] : H s ) ∩ C 1 ( ( 0 , T 0 ) : H s − 3 ) , 〈 x 〉 N ∂ x α u ∈ C ( [ 0 , T 0 ] : L 2 ) , | α | ≤ s 1 . (1.32)
Moreover, if ( u 0 , f ) ∈ H s ′ ( R n ) × L 1 ( R + : H s ′ ( R n ) )   with s ′ > s   then  1.32 holds with s ′   instead of s   in the same time interval [ 0 , T 0 ]   .
RemarksHere, we are not concerned with the problem of estimating the optimal values of s , s 1   or N   in Theorem  1.3 .
Also we shall not attempt to obtain the sharp form of the persistence property of the solution, (i.e. u ∈ C ( [ 0 , T 0 ] : H s ( R n ) )   ) as well as the continuous dependence of the solution upon the data. These can be established by combining the argument in [2] with our key a priori estimates in Lemma 2.1.
Similarly, from our arguments it is easy to deduce that the local solution possesses the local smoothing effect, i.e. if u 0 ∈ H s ( R n )   , then J s + 1 / 2 u ∈ L 2 ( R n × [ 0 , T 0 ] : 〈 x 〉 − N ~ d x d t )   .
The use of the weights in Theorem  1.3 comes from two sources. First, in order to convert L 1   conditions, such as  1.7 , into L 2   conditions. Secondly, one needs them in order to maintain the asymptotic flatness condition (H2), when 0 ⃗   is replaced by ( u , u ¯ , ∇ x u , ∇ x u ¯ )   , for a solution u   . As it was mentioned quasilinear results as those in Theorem  1.3 cannot be obtained by using just a fixed point argument. Instead, as in [20] , we shall rely on the artificial viscosity method. First, we consider the linear problem
∂ t u = − ε Δ 2 u − i ∂ x j ( a j k ( x , t ) ∂ x k u ) + b ⃗ 1 ( x , t ) ⋅ ∇ u + b ⃗ 2 ( x , t ) ⋅ ∇ u ¯ + c 1 ( x , t ) u + c 2 ( x , t ) u ¯ + f ( x , t ) ≡ − ε Δ 2 u + L ( x , t ) u + f ( x , t ) , u ( x , 0 ) = u 0 ( x ) , (1.33)
with ε ∈ ( 0 , 1 )   . The main step is to obtain the following a priori estimate for solutions of the linear IVP  1.33 : there exits T > 0   such that
sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 ≤ c T ( ∥ u 0 ∥ 2 2 + ∫ 0 T ∥ f ( t ) ∥ 2 2 d t ) , (1.34)
with T , c   independent of ε   and depending on an appropriate manner on the coefficients in  1.33 .
The inequality  1.34 will be proved under general hypotheses on the coefficients in  1.33 .
This will allow us to find a class of functions such that when the coefficients of the equation in  1.1 are evaluated in an element in this class they satisfy these general hypotheses.
We denote by u ε   the solution of the following nonlinear IVP associated to equation in  1.1 
∂ t u = − ε Δ 2 u − i ∂ x j ( a j k ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ∂ x j k u ) + b ⃗ 1 ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ⋅ ∇ u + b ⃗ 2 ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ⋅ ∇ u ¯ + c 1 ( x , t , u , u ¯ ) u + c 2 ( x , t , u , u ¯ ) u ¯ + f ( x , t ) , u ( x , 0 ) = u 0 ( x ) . (1.35)
The viscosity method provides the solution u ε   , in the appropriate class, in the time interval [ 0 , T ε ]   with T ε = O ( ε )   . Evaluating the coefficients in  1.35 in this solution we get a linear problem as that in  1.33 for which the estimate  1.34 holds. This a priori estimate allows us to extend the solution u ε   , in the same class, to a time interval [ 0 , T 0 ]   with T 0   independent of ε ∈ ( 0 , 1 )   .
Once the estimate  1.34 is available the proof of Theorem  1.3 follows an argument quite similar to that explained in detail in [20] . So we shall concentrate in the proof of the inequality  1.1 . This will be given in section 2, Lemma  2.1 .
Finally, we point out another difference between the elliptic quasi-linear case and the non-degenerate one.
In [20] the following general class of quasilinear equation was considered
∂ t u = − i ∂ x j ( a j k ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ∂ x k u + ∂ x j ( b j k ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ∂ x k u + b ⃗ 1 ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ⋅ ∇ u + b ⃗ 2 ( x , t , u , u ¯ , ∇ u , ∇ u ¯ ) ⋅ ∇ u ¯ + c 1 ( x , t , u , u ¯ ) u + c 2 ( x , t , u , u ¯ ) u ¯ + f ( x , t ) , (1.36)
where ( b j k )   is a symmetric complex valued matrix.
Under the ellipticity assumption : given r 0 > 0   there exists γ r 0 > 0   such that for any ( x , t , z ⃗ ) ∈ R n × R × ( ( B r 0 ( 0 ) ) ¯ 2 n + 2 ,   it follows that 〈 a j k ( x , t , z ⃗ ) ξ , ξ 〉 − | 〈 b j k ( x , t , z ⃗ ) ξ , ξ 〉 | ≥ γ r 0 | ξ | 2 , ∀ ξ ∈ R n ,   with B r 0 ( 0 )   as in  1.22 and similar assumptions on the asymptotic flatness, the growth of first order coefficients, the regularity of the coefficients, and a non-trapping hypothesis it was established in [20] that the IVP associated to the equation  1.36 is locally well posed.
The non-trapping hypothesis in [20] was the following : Given u 0 ∈ H r ( R n ) , r > n / 2 + 2   , define ϖ 2 ( x , 0 , ξ ) = − a j k ( x , 0 , u 0 , u ¯ 0 , ∇ u 0 , ∇ u ¯ 0 ) ξ k ξ j ,   ϖ 3 ( x , 0 , ξ ) = − b j k ( x , 0 , u 0 , u ¯ 0 , ∇ u 0 , ∇ u ¯ 0 ) ξ k ξ j ,   and κ ( x , ξ ) = ϖ 2 2 ( x , 0 , ξ ) − | ϖ 3 ( x , 0 , ξ ) | 2 .   It was said that u 0   satisfies the non-trapping hypothesis if there exists 0 < η < 1   and functions a ( x , ξ ) , a 1 ( x , ξ )   such that κ ( x , ξ ) = a ( x , ξ ) + η a 1 ( x , ξ ) ,   (with a ( x , ξ )   real, homogeneous of degree 2   , with | ∂ x β a ( x , ξ ) | ∈ C 1 , 1 ( R n × R n )   , | β | ≤ N ( n )   and θ ( ξ ) a ( x , ξ ) ∈ C N ( n ) ( R n × R n )   , where θ ≡ 1   for | ξ | > 1   , θ ≡ 0   for | ξ | < 1 / 2   , θ ∈ C ∞ ( R n )   , with a 1   verifying similar estimates), the Hamiltonian flow H a   associated to the symbol a   is non-trapping, for more details see [20] .
In the non-degenerate setting the equation in  1.36 does not allow for a more general one that that considered in  1.1 . This is due to the following linear algebra result (whose proof follows by induction).
Lemma 1.1. Let A , B   be two n × n   matrices, A   be real symmetric non-positive and non-degenerate, and B   a symmetric complex valued one such that 〈 A ξ , ξ 〉 = 0 implies 〈 B ξ , ξ 〉 = 0 .   Then there exists λ ∈ C   such that B = λ A   .

2 The linear problem

In this section we shall consider the linear IVP
∂ t u = − i ∂ x j ( a j k ( x , t ) ∂ x k u ) + b ⃗ 1 ( x , t ) ⋅ ∇ u + b ⃗ 2 ( x , t ) ⋅ ∇ u ¯ + c 1 ( x , t ) u + c 2 ( x , t ) u ¯ + f ( x , t ) ≡ L ( x , t ) u + f ( x , t ) , u ( x , 0 ) = u 0 ( x ) , (2.1)
where x ∈ R n , n > 1   , t ∈ [ 0 , T ]   , with T > 0   , and its associated ε   -viscosity version
∂ t u = − ε Δ 2 u + L ( x , t ) u + f ( x , t ) , ε ∈ ( 0 , 1 ) , u ( x , 0 ) = u 0 ( x ) , (2.2)
under the following assumptions:
(H l   1) Non-degeneracy: A ( x , t ) = ( a j k ( x , t ) ) j , k = 1 , . . , n   is a real symmetric matrix and there exist γ , γ 0 ∈ ( 0 , 1 )   such that for any ξ ∈ R n , ( x , t ) ∈ R n × [ 0 , T ]  
γ | ξ | ≤ | A ( x , t ) ξ | ≤ γ − 1 | ξ | , γ 0 | ξ | ≤ | A ( x , 0 ) ξ | ≤ γ 0 − 1 | ξ | , (2.3)
with
A ( x , 0 ) = ( a j k ( x , 0 ) ) j , k = 1 , . . , n = A 0 ( x ) + A h = ( a 0 , j k ( x ) ) j , k = 1 , . . , n + ( a j k 0 ) j , k = 1 , . . , n (2.4)
where for some N , M ~ ∈ Z +   large enough
| ∂ x α a 0 , j k ( x ) | ≤ c 0 〈 x 〉 N , | α | ≤ M ~ , j , k = 1 , . . , n , (2.5)
and A h   as in  1.12 .
(H l   2) Asymptotic Flatness: There exists c > 0   such that for all ξ ∈ R n , ( x , t ) ∈ R n × [ 0 , T ]   , j , k = 1 , . . , n   , and 0 < | α | ≤ M ~ , 0 ≤ | α ′ | ≤ M ~  
| ∂ x α a j k ( x , t ) | + | ∂ t ∂ x α ′ a j k ( x , t ) | ≤ c 〈 x 〉 N , (2.6)
with N , M ~   as in H l 1   .
(H l   3) Growth Condition of the First Order Coefficients: There exist c , c 0 > 0   such that for all ξ ∈ R n , ( x , t ) ∈ R n × [ 0 , T ]  
| ∂ x α b m j ( x , 0 ) | ≤ c 0 〈 x 〉 N , | α | ≤ M ~ , m = 1 , 2 , j = 1 , . . , n , (2.7)
and
| ∂ x α b m j ( x , t ) | ≤ c 〈 x 〉 N , | α | ≤ M ~ , m = 1 , 2 , j = 1 , . . , n , (2.8)
with N , M ~   as in H l 1   . (H l   4) Regularity of the Coefficients:
a j k , b m j , c m ∈ C b J ( R n × [ 0 , T ] ) , j , k = 1 , . . , n , m = 1 , 2 , (2.9)
with J = J ( n ) ∈ Z +   sufficiently large such that the proofs below involving ψ   .d.o's can be carried out, with norm
∥ ∂ x α ∂ t r d ∥ L ∞ ( R n × [ 0 , T ] ) ≤ c J , | α | + r ≤ J . (2.10)
with d = a j k , b m j   or c m   , j , k = 1 , . . , n , m = 1 , 2   .
(H l   5) Non-trapping Condition: The bicharacteristic flow associated to the symbol of A ( x , 0 )   (see  2.4 ), i.e.
h 2 0 ( x , ξ ) = ∑ j , k = 1 n a j k ( x , 0 ) ξ j ξ k (2.11)
is non-trapping.
We shall use c 0   to denote a generic constant which only depends on the coefficients evaluated at time t = 0   .
The main ingredient in the proof of our main result Theorem  1.3 is the following estimate for the solution of the linear IVP  2.1 .
Lemma 2.1. There exist c ~ 0 = c ~ 0 ( γ 0 , c 0 , n ) > 0   with γ 0   as in  2.3 , c 0   as in  2.5 and  2.7 , and the constant c 0   in Lemma  2.4 below, N ~ = N ~ ( n )   , and K 0 > 0 , T 0 ∈ ( 0 , T ]   depending on  2.3 - 2.10 such that for T ~ ∈ ( 0 , T 0 )   the solution of the IVP  2.2  u ε   satisfies
sup 0 ≤ t ≤ T ~ ∥ u ε ( t ) ∥ 2 + ∫ 0 T ~ ∫ | J 1 / 2 u ε | 2 〈 x 〉 − N ~ d x d t ≤ c ~ 0 e K 0 T ~ ( ∥ u 0 ∥ 2 + ∫ 0 T ~ ∥ f ( ⋅ , t ) ∥ 2 d t ) . (2.12)
Moreover,  2.12 still holds if we replace its last term by ( ∫ 0 T ~ ∫ | J − 1 / 2 f ( x , t ) | 2 〈 x 〉 − N ~ d x d t ) 1 / 2 .  
We recall the class of ψ   .d.o's introduced in [21] :
Ψ a f ( x ) = ∫ e i x ⋅ ξ a ( x , ξ ) f ^ ( ξ ) d ξ ,   with symbol a ( x , ξ ) = χ ( | ξ | ) a ( P ( x , A h ξ ) ; x , ξ ) ,   with χ   as in  1.15 , P ( x , A h ξ ) = x − ( x ⋅ A h ξ | ξ | 2 ) A h ξ ,   i.e. P ( x , A h ξ )   is the projection of x   into the hyperplane perpendicular to A h ξ   , and a = a ( s ; x , ξ ) ∈ S ( R n : S 1 , 0 m ) ,   S ( ⋅ )   denoting the Schwartz class, and S 1 , 0 m   the class of classical symbol of ψ   .d.o's of order m   .
In [21] (Theorem 3.2.1) we prove that
∥ Ψ a f ∥ 2 ≤ c ∥ f ∥ m , 2 . (2.13)
We shall use that our class is closed under differentiation of the symbol in the x   -variable.
Lemma 2.2. Let a α ( x , ξ ) = ∂ x α ( a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) ) ,   with a ( s ; x , ξ )   as above. Then a α ( x , ξ )   defines a symbol in our class. Moreover,
∥ Ψ a α f ∥ 2 ≤ c α ∥ f ∥ m , 2 . (2.14)
  • Proof. First we consider the case α = ( 1 , 0 , . . , 0 )   . So
    ∂ x 1 a ( P ( x , A h ξ ) ; x , ξ )
    = ∂ a ∂ x 1 ( P ( x , A h ξ ) ; x , ξ ) + ∑ j = 1 n ∂ a ∂ z j ( P ( x , A h ξ ) ; x , ξ ) ∂ P ∂ x 1 ( x , A h ξ ) .
    Since ∂ x 1 ( P ( x , A h ξ ) ) j = ∂ x 1 ( x j − x ⋅ A h ξ | ξ | 2 ( A h ξ ) j ) = δ 1 j − ( A h ξ ) 1 ( A h ξ ) j | ξ | 2 ,   it follows that (see remark 3.1.3(b) in [21] ) ∂ x 1 a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) = b 1 ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) ) ,   with b 1 = b 1 ( z ; x , ξ ) ∈ S ( R n : S 1 , 0 m )   which yields the result.
    The proof of the general case combines the above argument and induction in | α |   . □
Using the notation
ℒ = ℒ ( x , t ) = − ∂ x j ( a j k ( x , t ) ∂ x k ⋅ ) , (2.15)
and taking complex conjugate in the equation  2.2 we obtain the system
∂ t w ⃗ = − ε Δ 2 I w ⃗ + i H w ⃗ + B w ⃗ + C w ⃗ + F ⃗ , w ⃗ ( x , 0 ) = w ⃗ 0 ( x ) , (2.16)
where
w ⃗ = ( u u ¯ ) , F ⃗ = ( f f ¯ ) , (2.17)
H = ( ℒ 0 0 − ℒ ) , B = ( Ψ b 1 b ⃗ 2 ( x , t ) ⋅ ∇ b ⃗ 2 ( x , t ) ¯ ⋅ ∇ Ψ b ¯ 1 ) , (2.18)
with b 1 ∈ S 1 , 0 1   with odd symbol in ξ   , and C   a 2 × 2   matrix of ψ   .d.o's of order zero.
At some point we will take derivative of the equation in  1.1 , so the new coefficients of ∇ u   will be a combination of the original ones and some derivatives of the a j k   's. In fact, these coefficients depend on the order of the derivative just as a multiplicative constant.
For this reason we consider a general b 1 ∈ S 1 , 0 1   with odd symbol in ξ   in  2.15 ,  2.18 .
Lemma 2.3. There exists N ∈ Z +   such that for any T 0 ∈ ( 0 , T ]   and any ε ∈ ( 0 , 1 )   the solution w ⃗ ε = w ⃗   of the IVP  2.16 satisfies
∫ 0 T 0 ∫ | J 1 / 2 w ⃗ | 2 〈 x 〉 − N ~ d x d t ≤ ( c 0 + c T 0 ) sup 0 ≤ t ≤ T 0 ∥ w ⃗ ( t ) ∥ 2 2 + c 0 ε ∫ 0 T 0 ∥ Δ w ⃗ ( t ) ∥ 2 2 d t + c 0 ∫ 0 T 0 ∥ F ⃗ ( t ) ∥ 2 2 d t , (2.19)
where c 0   depends only on the coefficients evaluated at t = 0   and on Lemma 2.1 and c   depends on the estimates in  2.3 - 2.10 .
To prove Lemma  2.3 we will follow the argument in [21] (Lemma 5.1.3). First, we recall Lemma 5.1.1 in [21] .
Lemma 2.4. Let A ( x , 0 )   be as in  2.4 . Assume that the bicharacteristic flow is non-trapped, i.e.
{ X ( s ; x 0 , ξ 0 ) : s ∈ R + }   is unbounded for any ( x 0 , ξ 0 ) ∈ R n × R n − { 0 }   . Suppose λ ∈ L 1 ( [ 0 , ∞ ) ) ∩ C ( [ 0 , ∞ ) )   is strictly positive and non-increasing. Then there exist c 0 > 0   and a real symbol p ∈ S 1 , 0 0   , both depending on h 2 0 ( x , ξ )   in  2.11 and λ   such that H h 2 0 p = { h 2 0 , p } ( x , ξ ) ≥ λ ( | x | ) | ξ | − c 0 , ∀ ( x , ξ ) ∈ R n × R n .  
We fix N   as in  2.5 ,  2.6 ,  2.7 , and  2.8 , and choose λ ( ρ ) = 1 / ( 1 + ρ 2 ) N ~ / 2 = 〈 ρ 〉 N ~ .   to obtain the following time dependent version of Lemma  2.4 .
Lemma 2.5. These exists T 0 > 0   depending only on H l   2 such that for any ( x , ξ ) ∈ R n × R n  
H h 2 p = { h 2 ; p } ( x , t , ξ ) ≥ | ξ | 2 〈 x 〉 N ~ − 2 c 0 , (2.20)
where
h 2 = h 2 ( x , t , ξ ) = a j k ( x , t ) ξ k ξ j , (2.21)
and p ∈ S 1 , 0 0   is the operator in Lemma  2.4 .
  • Proof. We recall the notation h 2 0 ( x , ξ ) = a j k ( x , 0 ) ξ k ξ j .   By definition
    H h 2 p = ∂ ξ j h 2 ∂ x j p − ∂ x j h 2 ∂ ξ j p
    = 2 a j k ( x , t ) ξ k ∂ x j p − ∂ x j a j k ( x , t ) ξ k ξ j ∂ ξ j p .
    From our hypothesis H l   2 it follows that for any t ∈ ( 0 , T 0 ]  
    | ∂ ξ j h 2 ( x , t , ξ ) − ∂ ξ j h 2 0 ( x , ξ ) |
    = 2 | ( a j k ( x , t ) − a j k ( x , 0 ) ) ξ k | ≤ c T 0 〈 x 〉 N | ξ | ,
    and
    | ∂ x j h 2 ( x , t , ξ ) − ∂ x j h 2 0 ( x , ξ ) |
    = | ( ∂ x j a j k ( x , t ) − ∂ x j a j k ( x , 0 ) ) ξ k ξ j | ≤ c T 0 〈 x 〉 N | ξ | 2 .
    Since p ∈ S 1 , 0 0   choosing N ≫ N ~   one has that | H h 2 p − H h 2 0 p | ≤ c T 0 〈 x 〉 N ~ | ξ | .   From Lemma  2.4 we know that H h 2 0 p ( x , ξ ) ≥ | ξ | 〈 x 〉 N ~ − c 0 ,   so taking T 0   sufficiently small we obtain the desired result  2.20 . □
Now we shall prove Lemma  2.3 .
  • Proof. Let
    k ( x , ξ ) = ( e x p ( p ( x , ξ ) ) 0 0 − e x p ( p ( x , ξ ) ) ) , (2.22)
    where p ∈ S 1 , 0 0   is the symbol in Lemma  2.5 so K = Ψ k   is a diagonal matrix of ψ   .d.o's of order zero. We calculate
    ∂ t 〈 K w ⃗ , w ⃗ 〉 = 〈 K ∂ t w ⃗ , w ⃗ 〉 + 〈 K w ⃗ , ∂ t w ⃗ 〉 = − ε ( 〈 K Δ 2 I w ⃗ , w ⃗ 〉 + 〈 K w ⃗ , Δ 2 I w ⃗ 〉 ) + 〈 ( i [ K H − H K ] + K B + B * K ) w ⃗ , w ⃗ 〉 〈 ( K C + C * K ) w ⃗ , w ⃗ 〉 + 〈 K F ⃗ , w ⃗ 〉 + 〈 K w ⃗ , F ⃗ 〉 . (2.23)
    We disregard the symbols of order zero and consider first the symbol of i [ K H − H K ] + K B + B * K   , i.e.
    σ ( i [ K H − H K ] + K B + B * K ) ,   which is equal up to a symbol of order zero to
    − e p ( H h 2 p 0 0 H h 2 p ) + e p ( b 1 ( x , t , ξ ) 2 i b ⃗ 2 ( x , t ) ⋅ ξ − 2 i b ¯ ⃗ 2 ( x , t ) ⋅ ξ b 1 ( x , t , ξ ) ¯ ) (2.24)
    We write b ⃗ 2 ( x , t ) = b ⃗ 2 ( x , 0 ) + t b ⃗ 2 ( x , t ) − b ⃗ 2 ( x , 0 ) t ,   with a similar identity for b 1 ( x , t , ξ )   . Also since K = Ψ k   has order zero it is easy to see that | − ε ( 〈 K Δ 2 I w ⃗ , w ⃗ 〉 + 〈 K w ⃗ , Δ 2 I w ⃗ 〉 ) | ≤ c 0 ε ∥ w ⃗ ( t ) ∥ 2 , 2 2 .   This combined with the matrix version of the sharp Garding inequality, the hypothesis H l   3, and Lemma  2.5 yields, after integration by part, the desired result, i.e. inequality  2.19 (for details we refer to the proof of Lemma 5.1.3 in [21] ).
    □
Next, we shall recall some notations and definitions used in [21] . First, we define
a j k R ( x , 0 ) = θ ( x R ) a j k ( x , 0 ) + ( 1 − θ ( x R ) ) a j k 0 , (2.25)
where θ ∈ C 0 ∞ ( R n )   , θ ( x ) = 1   for | x | < 1   , and θ ( x ) = 0   for | x | > 2   ,
ℒ R ( x ) = ∂ x j ( a j k R ( x , 0 ) ∂ x k ) , and ℒ ( x , t ) = ∂ x j ( a j k ( x , t ) ∂ x k ) . (2.26)
We apply the operator ( K R ) *   , which is independent of t   , to the equation in  2.2 and use that ℒ ( x , t ) = ℒ ( x , t ) − ℒ ( x , 0 ) + ℒ R ( x ) + ℰ R ( x ) ,   with ℰ R ( x ) = ∂ x j ( ( 1 − θ ( x R ) ) ( a j k ( x , 0 ) − a j k 0 ) ∂ x k ) ,   and b ⃗ j ( x , t ) ⋅ ∇ = ( b ⃗ j ( x , t ) − b ⃗ j ( x , 0 ) ) ⋅ ∇ + b ⃗ j ( x , 0 ) ⋅ ∇ , j = 1 , 2 .   We shall need the following symbols used in [21] ,
b 1 ( x , ξ ) = − b 1 ( x , 0 ) i ξ , p R ( x , ξ ) = − χ ( | ξ | ) ∫ − ∞ 0 b 1 ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) ) d s , p e R ( x , ξ ) = 1 2 ( p R ( x , ξ ) + p R ( x , − ξ ) ) , k R ( x , ξ ) = exp ( p e R ( x , ξ ) ) , k R ( x , ξ ) = exp ( − p e R ( x , ξ ) ) . (2.27)
where ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) ) ,   denotes the bicharacteristic flow associated to the symbol of the truncated operator ℒ R ( x )   defined in  2.26 and χ   as in  1.15 .
We remark that after taking s   -derivatives in the equation  1.1 and linearizing the resulting equation the new term b 1 ( x , ξ )   obtained has the form s ∑ j , k , l ∂ x j a j k R ( x , 0 ) ξ j ξ k ξ l 〈 ξ 〉 − 2 − b 1 ( x , 0 ) i ξ .   From our hypotheses it is clear that this new b R ( x , ξ )   satisfies similar estimates than that in  2.27 . Thus, from  2.2 we get
∂ t ( K R ) * u = − ε ( K R ) * Δ 2 u − i [ ( K R ) * ; ℒ R ] u + i ℒ R ( K R ) * u + i [ ( K R ) * ; ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ] u + i ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ( K R ) * u + i ℰ R ( K R ) * u + i [ ℰ R ; ( K R ) * ] u + ( K R ) * b ⃗ 1 ( x , 0 ) ⋅ ∇ u + ( K R ) * ( b ⃗ 1 ( x , t ) − b ⃗ 1 ( x , 0 ) ) ⋅ ∇ u + ( K R ) * b ⃗ 2 ( x , 0 ) ⋅ ∇ u ¯ + ( K R ) * ( b ⃗ 2 ( x , t ) − b ⃗ 2 ( x , 0 ) ) ⋅ ∇ u ¯ + ( K R ) * f + terms of order zero in u . (2.28)
Next, we shall estimate
∂ t 〈 ( K R ) * u , ( K R ) * u 〉 = d d t ∫ | ( K R ) * u | 2 ( x , t ) d x . (2.29)
After using the equation  2.28 to estimate  2.29 we separate the terms obtained into four groups I K ; j , j = 1 , . . , 4   .
In I K ; 1   we set the terms with coefficients independent of t   , i.e. those involving the operators ℒ R ( x ) , b ⃗ 1 ( x , 0 ) ⋅ ∇ , b ⃗ 2 ( x , 0 ) ⋅ ∇ , ℰ R ( x ) ,   in the equation.
In I K ; 2   one has the terms involving the difference of the coefficients at the time t   and time 0   , i.e. the terms containing the operators ( ℒ ( x , t ) − ℒ ( x , 0 ) ) , ( b ⃗ 1 ( x , t ) − b ⃗ 1 ( x , 0 ) ) ⋅ ∇ , ( b ⃗ 2 ( x , t ) − b ⃗ 2 ( x , 0 ) ) ⋅ ∇ .   I K ; 3   contains the terms coming from the ε   viscosity part of the equation, i.e.
− ε 〈 ( K R ) * Δ 2 u , ( K R ) * u 〉 − ε 〈 ( K R ) * , ( K R ) * Δ 2 u 〉 .   In I K ; 4   we collect the terms of order zero and those coming from the inhomogeneous term f ( x , t )   . From the L 2   -continuity of ( K R ) *   they are bounded by
c R 2 N 0 ∥ u ( t ) ∥ 2 2 + c R 2 N 0 ∥ u ( t ) ∥ 2 ∥ f ∥ 2 , (2.30)
with N 0   depending only on the dimension n   . Similarly, we shall estimate
∂ t 〈 E R u , E R u 〉 = d d t ∫ | E R u | 2 ( x , t ) d x . (2.31)
where I = E R + K ~ R ( K R ) * ,   see Lemma 5.2.6 in [21] .
Thus,
∂ t E R u = − ε E R Δ 2 u + i E R ℒ ( x , t ) u + E R b ⃗ 1 ( x , t ) ⋅ ∇ u + E R b ⃗ 2 ( x , t ) ⋅ ∇ u ¯ + E R c 1 ( x , t ) u + E R c 2 ( x , t ) u ¯ + E R f . (2.32)
Inserting the identity
E R ℒ ( x , t ) = E R ℒ ( x , 0 ) + E R ( ℒ ( x , t ) − ℒ ( x , 0 ) ) = ℒ ( x , 0 ) E R + [ E R ; ℒ ( x , 0 ) ] + ( ℒ ( x , t ) − ℒ ( x , 0 ) ) E R + [ E R ; ℒ ( x , t ) − ℒ ( x , 0 ) ] , (2.33)
into the equation  2.32 and  2.31 we split the terms obtained into four groups I E ; j , j = 1 , . . , 4   .
In I E ; 1   we place the terms with coefficients independent of t   coming from the expressions ℒ R ( x ) , b ⃗ 1 ( x , 0 ) ⋅ ∇ , b ⃗ 2 ( x , 0 ) ⋅ ∇ , ℰ R ( x ) ,   in the equation.
I E ; 2   contains the terms involving the operator ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ,   (in this case the terms involving the operators E R Ψ b 1 , E R b ⃗ 2 ( x , t ) ⋅ ∇   can be directly bounded by using Lemma 5.2.6 in [21] ).
I E ; 3   contains the terms coming from the ε   viscosity part of the equation, i.e.
− ε 〈 E R Δ 2 u , E R u 〉 − ε 〈 E R , E R Δ 2 u 〉 .   In I E ; 4   we collect the remainder terms of order zero which using Lemma 5.2.6 in [21] are bounded by
c R 2 N 0 ∥ u ( t ) ∥ 2 2 + c R 2 N 0 ∥ u ( t ) ∥ 2 ∥ f ∥ 2 , (2.34)
with N 0   as in  2.30 .
The terms in I K ; 1   and I E ; 1   (all involving time independent coefficients) were already considered in Section 5 of [21] , (i.e. their contribution to the equations of  2.29 and  2.31 ).
Thus, to bound them we shall use the following inequalities proved in Lemmas 4.2.4-5.2.6 in [21] .
Proposition 2.1. There exist M = M ( N ) ∈ Z +   , N   denoting the decay of the coefficients in  2.5 - 2.8 , with M ( N ) ↑ ∞   as N ↑ ∞   , and R 0   large enough such that for R ≥ R 0   one has
ℛ e i 〈 ℒ R ( x ) ( K R ) * u , ( K R ) * u 〉 ≡ 0 , (2.35)
ℛ e i 〈 ℒ ( x , 0 ) E R u , E R u 〉 ≡ 0 , (2.36)
| 〈 i [ ℒ R ; ( K R ) * ] u + ( K R ) * b ⃗ 1 ( x , 0 ) ⋅ ∇ u , ( K R ) * u 〉 | + | 〈 ℰ ℛ ( K R ) * u , ( K R ) * u 〉 | + | 〈 [ ℰ R ; ( K R ) * ] u , ( K R ) * u 〉 | + | 〈 ( K R ) * b ⃗ 2 ( x , 0 ) ⋅ ∇ u ¯ , ( K R ) * u 〉 | ≤ c 0 R − M ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c 0 R N 0 + N 1 ( M ) ∥ u ( t ) ∥ 2 2 + c 0 ∥ f ∥ 2 2 . (2.37)
and
| 〈 [ E R ; ℒ ( x , 0 ) ] u , E R u 〉 | + | 〈 E R b ⃗ 1 ( x , 0 ) ⋅ ∇ u , E R u 〉 | + | 〈 E R b ⃗ 2 ( x , 0 ) ⋅ ∇ u ¯ , E R u 〉 | ≤ c 0 R N 0 ∥ u ( t ) ∥ 2 2 + c 0 ∥ f ∥ 2 2 , (2.38)
with N 0   as in  2.30 , N ~   depending only on n   , with N ≫ N ~   , and N 1   depending only on M   .
To handle the contribution coming from the terms in I K ; 2   we shall establish the following inequalities.
Proposition 2.2. Take N 0   as in  2.30 , then
| 〈 [ ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ; ( K R ) * ] u , ( K R ) * u 〉 | ≤ c T R N 0 ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c R N 0 ∥ u ( t ) ∥ 2 2 , (2.39)
ℛ e i 〈 ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ( K R ) * u , ( K R ) * u 〉 ≡ 0 , (2.40)
| 〈 ( K R ) * ( b ⃗ 1 ( x , t ) − b ⃗ 1 ( x , 0 ) ) ⋅ ∇ u , ( K R ) * u 〉 | ≤ c T R N 0 ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c R N 0 ∥ u ( t ) ∥ 2 2 , (2.41)
and
| 〈 ( K R ) * ( b ⃗ 2 ( x , t ) − b ⃗ 2 ( x , 0 ) ) ⋅ ∇ u ¯ , ( K R ) * u ¯ 〉 | ≤ c T R N 0 ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c R N 0 ∥ u ( t ) ∥ 2 2 , (2.42)
with N ~   depending only on n   and N ≫ N ~   .
  • Proof. We write
    ( [ ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ; ( K R ) * ] ) *
    = ( ℒ ( x , t ) − ℒ ( x , 0 ) ) * K R − K R ( ℒ ( x , t ) − ℒ ( x , 0 ) ) *
    where
    ( ℒ ( x , t ) − ℒ ( x , 0 ) ) *
    = ( a j k ( x , t ) − a j k ( x , 0 ) ) ∂ j k 2 + ∂ j ( a j k ( x , t ) − a j k ( x , 0 ) ) ∂ k
    = β j k ( x , t ) ∂ j k 2 + β ~ k ( x , t ) ∂ k = Θ 1 + Θ 2 ,
    with the notation ∂ j   instead of ∂ x j   . Notice that the coefficients β j k ( x , t )   and β ~ k ( x , t )   and their derivatives up to order p   (large enough) have strong decay at infinity uniformly in t ∈ [ 0 , T ]   . In fact, a sufficient number of semi-norms in the Schwartz class are bounded by c T   for all t ∈ [ 0 , T ]   .
    The term Θ 2   yields 〈 ( β ~ k ( x , t ) ∂ k K R − K R β ~ k ( x , t ) ∂ k ) u , ( K R ) * u 〉   which using Theorem 3.3.1 in [21] and the L 2   -continuity of K R   can be bounded by T R N 0 ∥ u ( t ) ∥ 2 2 , ∀ t ∈ ( 0 , T ) .   To handle Θ 1   we use again Theorem 3.3.1 in [21] to write that
    β j k ( x , t ) ∂ j k 2 K R − K R β j k ( x , t ) ∂ j k 2 = Ψ d R + zero order terms , (2.43)
    with d R ( x , t , ξ ) = β j k ( x , t ) ∂ ξ l ( ξ j ξ k ) ∂ x l K R − ∂ x l β j k ( x , t ) ξ j ξ k ∂ ξ l K R .   We recall that a symbol in our class a ( x , ξ )   when multiplied by φ ( x )   , a fast decaying function in the x   variable, becomes a classical symbol of the same order.
    To each term in the symbol of d R ( x , t , ξ )   we can apply the following argument.
    Claim: Let φ ( x , t )   be an smooth function with strong decay at infinity in the x   -variable uniformly in t ∈ [ 0 , T ]   , with a sufficient number of semi-norms in the Schwartz class bounded by c T   . Then
    | 〈 φ ( x , t ) ∂ j K R u , ( K R ) * u 〉 | ≤ c T R 2 N 0 ∥ ( J 1 / 2 u ( t ) ) 〈 x 〉 − N ~ ∥ 2 2 + ∥ u ( t ) ∥ 2 2 (2.44)
    NotationWe shall use the notation ≃   to denote terms whose difference can be bounded by a multiple of ∥ u ∥ 2 2   or by an operator of order zero.
    • Proof. We have
      〈 φ ( x , t ) ∂ j K R u , ( K R ) * u 〉
      = 〈 〈 x 〉 − N ~ 〈 x 〉 2 N ~ φ ( x , t ) ∂ j K R J 1 / 2 J − 1 J 1 / 2 u , 〈 x 〉 − N ~ ( K R ) * u 〉
      ≃ 〈 〈 x 〉 − N ~ J 1 / 2 〈 x 〉 2 N ~ φ ( x , t ) ∂ j K R J − 1 J 1 / 2 u , 〈 x 〉 − N ~ ( K R ) * u 〉
      ≃ 〈 〈 x 〉 − N ~ J 1 / 2 〈 x 〉 3 N ~ φ ( x , t ) ∂ j K R J − 1 〈 x 〉 − N ~ J 1 / 2 u , 〈 x 〉 − N ~ ( K R ) * u 〉
      ≃ 〈 J 1 / 2 〈 x 〉 − N ~ 〈 x 〉 3 N ~ φ ( x , t ) ∂ j K R J − 1 〈 x 〉 − N ~ J 1 / 2 u , 〈 x 〉 − N ~ ( K R ) * u 〉
      ≃ 〈 〈 x 〉 3 N ~ φ ( x , t ) ∂ j K R J − 1 〈 x 〉 − N ~ J 1 / 2 u , 〈 x 〉 − N ~ J 1 / 2 〈 x 〉 − N ~ ( K R ) * u 〉
      = Λ 3 .
      Since 〈 x 〉 3 N ~ φ ( x , t ) ∂ j K R J − 1 ∈ S 1 , 0 0   , i.e. a classical ψ   .d.o. of order zero, it follows that | Λ 3 | ≤ c T R N 0 ∥ 〈 x 〉 − N ~ J 1 / 2 u ∥ 2 ∥ 〈 x 〉 − N ~ J 1 / 2 〈 x 〉 − N ~ ( K R ) * u ∥ 2 .   Finally, using that ( 〈 x 〉 − N ~ J 1 / 2 〈 x 〉 − N ~ ( K R ) * ) * = K R 〈 x 〉 − N ~ J 1 / 2 〈 x 〉 − N ~ ,   we obtain the desired inequality  2.44 .
      Returning to the operator Ψ d R   whose symbol is d R ( x , ξ )   and applying our claim one gets that | 〈 Ψ d R u , ( K R ) * u 〉 | ≤ c T R N 0 ∥ 〈 x 〉 − N ~ J 1 / 2 u ( t ) ∥ 2 2 + ∥ u ( t ) ∥ 2 2 ,   which proves  2.39 .
      Integration by parts yields  2.40 .
      To prove  2.41 we write ( ( K R ) * ( b 1 ( x , t ) − b 1 ( x , 0 ) ) ⋅ ∇ ) ) * ≃ ( ( b 1 ( x , t ) − b 1 ( x , 0 ) ) ⋅ ∇ ) K R = G R * ∈ S 1 , 0 1 ,   and Ψ G R * ≃ Ψ G ¯ R   , so using an argument similar to that given in the proof of the claim above we get  2.41 . Similarly for  2.42 .
      □
    To bound the terms in I E ; 2   we have the following estimates.
    Proposition 2.3. There exists N 0 ∈ Z +   depending only on the dimension n   such that for any t ∈ ( 0 , T 0 ]  
    | 〈 [ ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ; E R ] u , E R u 〉 | ≤ c T R N 0 ∥ u ( t ) ∥ 2 2 , (2.45)
    and
    ℛ e i 〈 ( ℒ ( x , t ) − ℒ ( x , 0 ) ) E R u , E R u 〉 = 0 , (2.46)
    • Proof. With the same notation that in the previous proof we have ( ℒ ( x , t ) − ℒ ( x , 0 ) ) u = ∂ j ( β j k ( x , t ) ∂ k u ) .   So using that E R = I − K ~ R ( K R ) *   , one sees that
      [ ( ℒ ( x , t ) − ℒ ( x , 0 ) ) ; E R ] u
      = K ~ R [ ( K R ) * ; ∂ j β j k ( x , t ) ∂ k ] + [ K ~ R ; ∂ j β j k ( x , t ) ∂ k ] ( K R ) * .
      First we consider the term involving [ ( K R ) * ; ∂ j β j k ( x , t ) ∂ k ]   . Taking the adjoint it follows that
      ( [ ( K R ) * ; ∂ j β j k ( x , t ) ∂ k ] ) * = − K R ∂ j β j k ( x , t ) ∂ k + ∂ j β j k ( x , t ) ∂ k K R
      = − K R ∂ j ( β j k ( x , t ) ) ∂ k + ∂ j ( β j k ( x , t ) ) ∂ k K R
      − K R β j k ( x , t ) ∂ j k 2 + β j k ( x , t ) ∂ j k 2 K R
      ≃ − K R β j k ( x , t ) ∂ j k 2 + β j k ( x , t ) ∂ j k 2 K R ,
      since − K R ∂ j ( β j k ) ∂ k + ∂ j ( β j k ) ∂ k K R   is an operator of order zero, see Theorem 3.3.1 in [21] .
      Up to symbols of order zero (bounded operator in L 2   ) the symbol of β j k ( x , t ) ∂ j k 2 K R − K R β j k ( x , t ) ∂ j k 2   is equal to η ( x , t , ξ ) = β j k ∂ ξ l ( ξ k ξ j ) ∂ x l k R ( x , ξ ) − ∂ x l β j k ξ k ξ j ∂ ξ l k R ( x , ξ ) .   So η ( x , t , ξ ) ∈ S 1 , 0 1   uniformly in t ∈ [ 0 , T ]   and η ( x , t , ξ ) ≃ Ψ d β ~ l ( x , t ) ∂ l   with d = d ( x , ξ )   in our class. Similarly for η * ( x , t , ξ )   , then η * ( x , t , ξ ) ≃ β ~ l ( x , t ) ∂ l Ψ d 1   . Inserting this in  2.45 we get
      〈 K ~ R [ ( K R ) * ; ∂ j β j k ( x , t ) ∂ k ] u , E R u 〉 ≃ 〈 K ~ R β ~ l ( x , t ) ∂ l Ψ d 1 ( x , ξ ) u , E R u 〉
      ≃ 〈 Ψ b Ψ d 1 ( x , ξ ) u , E R u 〉 ≃ 〈 Ψ d 1 ( x , ξ ) u , ( Ψ b ) * E R u 〉 ,
      which can be bounded by ∥ u ( t ) ∥ 2 2   , by using that Ψ b   is a classical ψ   .d.o. of order 1 and the continuities properties of E R   deduce in Lemma 5.2.6 in [21] .
      A similar argument provides the bound for the term [ K ~ R , ∂ j β j k ( x , t ) ∂ k ] ( K R ) *   . Collecting this information we get  2.45 .
      The proof of  2.46 follows by integration by parts.
      □
    The terms coming from the artificial viscosity term ε Δ 2   , i.e. the terms in I K ; 3   and I E ; 3   will be handled by using the following inequalities.
    Proposition 2.4.
    〈 ε ( K R ) * Δ 2 u , ( K R ) * u 〉 = ε 〈 Δ ( K R ) * u , Δ ( K R ) * u 〉 + Λ 1 , (2.47)
    and
    〈 ε E R Δ 2 u , E R u 〉 = ε 〈 Δ E R u , Δ E R u 〉 + Λ 2 , (2.48)
    where | Λ j | ≤ ε R 2 N 0 ∥ Δ u ( t ) ∥ 2 ( ∥ u ( t ) ∥ 2 + ∥ ∇ u ( t ) ∥ 2 ) , j = 1 , 2 .  
    • Proof. To obtain  2.47 we write
      〈 ( K R ) * Δ 2 u , ( K R ) * u 〉 = 〈 [ ( K R ) * ; Δ ] Δ u , ( K R ) * u 〉 + 〈 Δ ( K R ) * Δ u , ( K R ) * u 〉 = 〈 [ ( K R ) * ; Δ ] Δ u , ( K R ) * u 〉 + 〈 ( K R ) * Δ u , Δ ( K R ) * u 〉 = 〈 [ ( K R ) * ; Δ ] Δ u , ( K R ) * u 〉 + 〈 ( K R ) * Δ u , [ Δ ; ( K R ) * ] u 〉 + 〈 ( K R ) * Δ u , ( K R ) * Δ u 〉 ≡ Ω 1 + Ω 2 + Ω 3 . (2.49)
      Thus, we consider [ Δ ; ( K R ) * ] = Δ ( K R ) * − ( K R ) * Δ = T ,   and its adjoint T * = − ( Δ K R − K R Δ ) .   Let k R ( x , ξ )   denote the symbol of K R   such that K R f ( x ) = ∫ e i x ⋅ ξ k R ( x , ξ ) f ^ ( ξ ) d ξ ,   so Δ K R f ( x ) = ∫ e i x ⋅ ξ { − | ξ | 2 k R ( x , ξ ) + 2 i ξ j ∂ x j k R ( x , ξ ) + Δ x k R ( x , ξ ) } f ^ ( ξ ) d ξ .   Therefore, T * f ( x ) = ∫ e i x ⋅ ξ { 2 i ξ j ∂ x j k R ( x , ξ ) + Δ x k R ( x , ξ ) } f ^ ( ξ ) d ξ .   From Lemma  2.2 above and Theorem 3.2.1 in [21] one has that the operators with symbols ∂ x j k R ( x , ξ )   and Δ x k R ( x , ξ )   are bounded in L 2   . Hence, we can write that
      T * = C j ∂ x j + C 0 , and T = C 0 * + ∂ x j ( − C j * ) , (2.50)
      with C j , j = 0 , 1 , . . , n   denoting L 2   -bounded operators. Also we have that
      ∫ e i x ⋅ ξ i ξ j ∂ x j k R ( x , ξ ) f ^ ( ξ ) d ξ
      = ∂ x j ( ∫ e i x ⋅ ξ ∂ x j k R ( x , ξ ) f ^ ( ξ ) d ξ ) − ∫ e i x ⋅ ξ ∂ x j x j 2 k R ( x , ξ ) f ^ ( ξ ) d ξ ,
      so that
      T * = ∂ x j C ~ j + C ~ 0 , and T = C ~ 0 * − C ~ j * ∂ x j , (2.51)
      with C ~ j , j = 0 , 1 , . . , n   denoting L 2   bounded operators.
      To bound Ω 1   we see that
      〈 [ ( K R ) * ; Δ ] Δ u , ( K R ) * u 〉 = − 〈 T Δ u , ( K R ) * u 〉
      = − 〈 ( C 0 * + ∂ x j ( − C j * ) ) Δ u , ( K R ) * u 〉
      = − 〈 C 0 * Δ u , ( K R ) * u 〉 − 〈 ( − C j * ) Δ u , ∂ x j ( K R ) * u 〉 .
      Since an explicit computation shows that ∂ x j ( K R ) * ≃ ( K R ) * ∂ x j ,   i.e. their difference is an L 2   -bounded operator, it follows that
      | Ω 1 | ≤ c 0 ∥ Δ u ∥ 2 ( ∥ u ∥ 2 + ∥ ( K R ) * ∇ u ∥ 2 )
      ≤ c 0 ∥ Δ u ∥ 2 ( ∥ u ∥ 2 + R N 0 ∥ ∇ u ∥ 2 ) .
      Using again  2.50 - 2.51 we have that
      〈 ( K R ) * Δ u , [ Δ ; ( K R ) * ] u 〉 = 〈 ( K R ) * Δ u , T u 〉 = 〈 ( K R ) * Δ u , ( C ~ 0 * − C ~ j * ∂ x j ) u 〉 ,
      so Ω 2   can be bounded as
      | Ω 2 | ≤ c 0 ∥ ( K R ) * Δ u ∥ 2 ( ∥ u ∥ 2 + ∥ ∇ u ∥ 2 )
      ≤ c 0 R N 0 ∥ Δ u ∥ 2 ( ∥ u ∥ 2 + ∥ ∇ u ∥ 2 ) .
      Inserting this information in  2.49 we obtain  2.47 . To obtain  2.48 we recall that E R = I − K ~ R ( K R ) *   so that E R Δ 2 = Δ E R Δ + [ E R ; Δ ] Δ ,   with
      [ E R ; Δ ] = − [ K ~ R ( K R ) * ; Δ ] = − ( K ~ R ( K R ) * Δ − Δ K ~ R ( K R ) * )
      = − K ~ R [ ( K R ) * ; Δ ] + [ K ~ R ; Δ ] ( K R ) * ≡ Γ 1 + Γ 2 .
      Using  2.50 - 2.51 and an explicit computation it follows that Γ 1 = K ~ R T = K ~ R ( C 0 * − ∂ x j C j * ) = − ∂ x j K ~ R C j * + Q 1 ,   where Q j , j = 1 , . . , 4   will denote L 2   bounded operators. So
      | 〈 Γ 1 Δ u , E R u 〉 | = | 〈 ( ∂ x j K ~ R C j * + Q 1 ) Δ u , E R u 〉 |
      ≤ c 0 ∥ Δ u ∥ 2 ( ∥ ∂ x j E R u ∥ 2 + ∥ E R u ∥ 2 ) .
      Combining ∂ x j E R = ∂ x j − ∂ x j K ~ R ( K R ) *   and an explicit computation one gets that
      ∥ ∂ x j E R u ∥ 2 = ∥ ( ∂ x j − ∂ x j K ~ R ( K R ) * ) u ∥ 2
      ≤ ∥ u ∥ 2 + ∥ ∂ x j K ~ R ( K R ) * u ∥ 2 ≃ ∥ u ∥ 2 + ∥ K ~ R ∂ x j ( K R ) * u ∥ 2
      ≤ ∥ u ∥ 2 + R N 0 ∥ ∂ x j ( K R ) * u ∥ 2 = ∥ u ∥ 2 + R N 0 ∥ K R ∂ x j u ∥ 2
      ≤ ∥ u ∥ 2 + R 2 N 0 ∥ ∂ x j u ∥ 2 ,
      and consequently | 〈 Γ 1 Δ u , E R u 〉 | ≤ c 0 R 2 N 0 ∥ Δ u ∥ 2 ( ∥ u ∥ 2 + ∥ ∇ u ∥ 2 ) .   For Γ 2   we reproduce the argument in  2.50 - 2.51 for K ~ R   instead of K R   . So using the same notation we have Γ 2 = [ K ~ R ; Δ ] ( K R ) * = ( C 0 + ∂ x j C j ) ( K R ) * ,   so as before
      | 〈 Γ 2 Δ u , E R u 〉 | = | 〈 ( C 0 + ∂ x j C j ) ( K R ) * Δ u , E R u 〉 |
      ≤ | 〈 C 0 ( K R ) * Δ u , E R u 〉 | + | 〈 C j ( K R ) * Δ u , ∂ x j E R u 〉 |
      ≤ ∥ Δ u ∥ 2 ( ∥ E R u ∥ 2 + ∥ ∂ x j E R u ∥ 2 ) ≤ c 0 R 2 N 0 ∥ Δ u ∥ 2 ( ∥ u ∥ 2 + ∥ ∇ u ∥ 2 ) .
      Collecting the results in Propositions  2.1 - 2.4 we get that
      d d t ∥ ( K R ) * u ( t ) ∥ 2 2 + ε ∥ ( K R ) * Δ u ( t ) ∥ 2 2 ≤ ( c 0 R − M + c T R N 0 ) ∥ ( J 1 / 2 u ) ( t ) 〈 x 〉 − N ~ ∥ 2 2 + c ( R N 0 + N 1 ( M ) + T R N 0 ) ∥ u ( t ) ∥ 2 2 + c 0 ∥ f ∥ 2 2 + c 0 ε R 2 N 0 ∥ Δ u ( t ) ∥ 2 ( ∥ u ( t ) ∥ 2 + ∥ ∇ u ( t ) ∥ 2 ) , (2.52)
      and
      d d t ∥ E R u ∥ 2 2 + ε ∥ E R Δ u ∥ 2 2 ≤ c ( R N 0 + T R N 0 ) ∥ u ( t ) ∥ 2 2 + c ∥ f ∥ 2 2 + c 0 ε R 2 N 0 ∥ Δ u ( t ) ∥ 2 ( ∥ u ( t ) ∥ 2 + ∥ ∇ u ( t ) ∥ 2 ) . (2.53)
      We will use that E R = I − K ~ R ( K R ) *   so that
      ∥ v ∥ 2 = ∥ ( E R + K ~ R ( K R ) * ) v ∥ 2 ≤ ∥ E R v ∥ 2 + R N 0 ∥ ( K R ) * v ∥ 2 . (2.54)
      Also, we need the following interpolation estimates : for any l ∈ Z +  
      R l ∥ v ∥ 2 ∥ Δ v ( t ) ∥ 2 ≤ R 2 l ∥ v ∥ 2 2 + ∥ Δ v ∥ 2 2 , (2.55)
      and
      R l ∥ ∇ v ∥ 2 ∥ Δ v ∥ 2 ≤ c R l ∥ v ∥ 2 1 / 2 ∥ Δ v ∥ 2 3 / 2 ≤ c R 4 l ∥ v ∥ 2 2 + ∥ Δ v ∥ 2 2 . (2.56)
      So combining  2.52 - 2.53 with  2.54 - 2.56 we find that
      d d t ( ∥ ( K R ) * u ( t ) ∥ 2 2 + ∥ E R u ( t ) ∥ 2 2 ) + ε R − 2 N 0 ∥ Δ u ( t ) ∥ 2 2 ≤ ( c T R N 0 + c 0 R − M ) ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c ( R N 0 + N 1 ( M ) + T R N 0 ) ∥ u ( t ) ∥ 2 2 + c 0 ∥ f ( t ) ∥ 2 2 + c 0 ε R 2 N 0 ∥ Δ u ( t ) ∥ 2 ( ∥ u ( t ) ∥ 2 + ∥ ∇ u ( t ) ∥ 2 ) ≤ ( c T R N 0 + c 0 R − M ) ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c ( R N 0 + N 1 ( M ) + T R N 0 ) ∥ u ( t ) ∥ 2 2 + c 0 ∥ f ( t ) ∥ 2 2 + ε 2 R − 2 N 0 ∥ Δ u ( t ) ∥ 2 2 + c 0 ε R 14 N 0 ∥ u ( t ) ∥ 2 2 . (2.57)
      Thus,
      d d t ( ∥ ( K R ) * u ( t ) ∥ 2 2 + ∥ E R u ( t ) ∥ 2 2 ) + ε 2 R − 2 N 0 ∥ Δ u ( t ) ∥ 2 2 ≤ ( c T R N 0 + c 0 R − M ) ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 + c ( R N 0 + N 1 ( M ) + T R N 0 ) ∥ u ( t ) ∥ 2 2 + c 0 ∥ f ( t ) ∥ 2 2 + c 0 ε R 14 N 0 ∥ u ( t ) ∥ 2 2 . (2.58)
      Integrating the above inequality in the time interval ( 0 , T )   and inserting in the result the estimate  2.19 we get
      sup 0 ≤ t ≤ T ( ∥ ( K R ) * u ( t ) ∥ 2 2 + ∥ E R u ( t ) ∥ 2 2 ) + ε 2 R − 2 N 0 ∫ 0 T ∥ Δ u ( t ) ∥ 2 2 d t ≤ ∥ ( K R ) * u ( 0 ) ∥ 2 2 + ∥ E R u ( 0 ) ∥ 2 2 + ( c T R N 0 + c 0 R − M ) ∫ 0 T ∥ ( J 1 / 2 u ) 〈 x 〉 − N ~ ∥ 2 2 d t + ( c R N 0 + N 1 ( M ) + T R N 0 + c 0 ε R 14 N 0 ) ∫ 0 T ∥ u ( t ) ∥ 2 2 d t + c 0 ∫ 0 T ∥ f ( t ) ∥ 2 2 d t ≤ ∥ ( K R ) * u ( 0 ) ∥ 2 2 + ∥ E R u ( 0 ) ∥ 2 2 + ( c T R N 0 + c 0 R − M ) ( ( c 0 + c T ) sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 + c 0 ε ∫ 0 T ∥ Δ u ( t ) ∥ 2 2 d t + c 0 ∫ 0 T ∥ f ( t ) ∥ 2 2 d t ) + c ( R N 0 + N 1 ( M ) + T R N 0 + c 0 ε R 14 N 0 ) T sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 + c 0 ∫ 0 T ∥ f ( t ) ∥ 2 2 d t . (2.59)
      Since M ( N ) ↑ ∞   as N ↑ ∞   , we take N   in our hypotheses large enough such that M = 100 N 0   (we recall that N 0   depends only on the dimension n   ). Next, we fix R ≥ R 0   sufficiently large and then T = T ( N 0 , M , R ) > 0   small enough such that the following inequalities holds ( c T R N 0 + c 0 R − M ) ( c 0 + c T ) ≤ R − 2 N 0 4 ,   ( c R N 0 + N 1 ( M ) + c T R N 0 + c 0 R 14 N 0 ) T ≤ R − 2 N 0 4   and c 0 ( c T R N 0 + c 0 R − M ) ≤ R − N 0 / 4 .   Combining these inequalities,  2.59 , and  2.54 we get the estimate sup 0 ≤ t ≤ T ∥ u ( t ) ∥ 2 2 + ε 4 R − 2 N 0 ∫ 0 T ∥ Δ u ( t ) ∥ 2 2 d t ≤ c 0 R 2 N 0 ∥ u ( 0 ) ∥ 2 2 + c 0 ∫ 0 T ∥ f ( t ) ∥ 2 2 d t .   which proves  2.12 .
      □
    □
Acknowledgments C. E. K. and G. P. were supported by NSF grants. L. V. was supported by a DGICYT grant. This work was carried out at the IAS which support we acknowledge. C. E. K. was supported at IAS by The Von Newmann Fund, The Weyl Fund, The Oswald Veblem Fund and The Bell Companies Fellowship. References

  1. M. J. Ablowitz, R. Haberman, Nonlinear evolution equations in two and three dimensions, Phys. Rev. Lett. 35 (1975), 1185–1188.
  2. J. L. Bona, R. Smith, The initial value problem for the Korteweg-de Vries equation, Roy. Soc. London Ser A 278 (1978), 555–601.
  3. A. P. Calderón, R. Vaillancourt, On the boundedness of pseudo-differential operators, J. Math. Soc. Japan 23 (1971), 374–378.
  4. P. Constantin, J-C. Saut, Local smoothing properties of dispersive equations, Journal A.M.S. 1 (1989), 413–446.
  5. H. Chihara, Local existence for semi-linear Schrödinger equations, Math. Japan 42 (1995), 35–51.
  6. W. Craig, T. Kappeler, W. Strauss, Microlocal dispersive smoothing for the Schrödinger equation, Comm. Pure Appl. Math. 48 (1995), 769–860.
  7. A. Davey, K. Stewartson, On three–dimensional packets of surface waves, Proc. Roy. Soc. London Ser. 338 (19740), 101–110.
  8. V. D. Djordjevic, L. G. Redekopp, On two-dimensional packets of capillary-gravity waves J. Fluid Mech. 79 (1977), 703–714.
  9. S. Doi, On the Cauchy problem for Schrödinger type equations and the regularity of solution J. Math. Kyoto Univ. 34 (1994), 319–328.
  10. S. Doi, Remarks on the Cauchy problem for Schrödinger–type equations Comm. Partial Differential Equations 21 (1996), 163–178.
  11. S. Doi, Smoothing effects for Schrodinger evolution equation and global behavior of geodesic flow Math. Ann. 318 (2000), 355–389.
  12. N. Hayashi, T. Ozawa, Remarks on nonlinear Schrödinger equations in one space dimension Differential Integral Equations 7 (1994), 453–461.
  13. L. Hörmander, Pseudo–differential operators and non–elliptic boundary problems, Ann. of Math. 2 (1966), 129–209.
  14. W. Ichinose, On L 2   wellposedness of the Cauchy problem for Schrödinger type equations on a Riemannian manifold and Maslov theory, Duke Math. J 56 (1988), 549–588.
  15. Y. Ishimori, Multi vortex solutions of a two dimensional nonlinear wave equation, Progr. Theor. Phys 72 (1984), 33–37.
  16. T. Kato, On the Cauchy problem for the (generalized) Korteweg-de Vries equation, Advances in Math. Supp. Studies, Studies in Applied Math. 8 (1983), 93–128.
  17. C. E. Kenig, G. Ponce, L. Vega, Oscillatory integrals and regularity of dispersive equations, Indiana University Math. J. 40 (1991), 33–69.
  18. C. E. Kenig, G. Ponce, L. Vega, Small solutions to nonlinear Schrödinger equations, Ann. Inst. H. Poincaré Anal. Non Linéaire 10 (1993), 255–288.
  19. C. E. Kenig, G. Ponce, L. Vega, Smoothing effects and local existence theory for the generalized nonlinear Schrödinger equations, Invent. Math. 134 (1998), 489-545.
  20. C. E. Kenig, G. Ponce, L. Vega, The Cauchy problem for the quasi-linear Schrödinger equations, Invent. Math. 158 (2004), 343-388.
  21. C. E. Kenig, G. Ponce, C. Rolvung, L. Vega, Variable coefficient Schrödinger flows for ultrahyperbolic operators to appear in Advances in Math.
  22. D. J. Kruzhkov, A. V. Faminskii, Generalized solutions for the Cauchy problem for the Korteweg-de Vries equation, Math. USSR Sb, 48 (1990), 93–138.
  23. S. Mizohata, On the Cauchy Problem Notes and Reports in Mathematics in Science and Engineering,3, Science Press and Academic Press (1985).
  24. C. Rolvung, Non-isotropic Schrödinger equations PhD. dissertation, University of Chicago, (1998).
  25. P. Sjölin, Regularity of solutions to the Schrödinger equation, Duke Math. J. 55 (1987), 699–715.
  26. E. M. Stein, Harmonic analysis: real–variable methods, orthogonality, and oscillatory integrals, Princeton University Press, (1993).
  27. C. Sulem, P-L. Sulem, The nonlinear Schrödinger Equation. Self-focusing and collapse, Applied Mathematical Sciences 139. Springer (1999).
  28. J. Takeuchi, Le problème de Cauchy pour certaines equations aux derivées partielles du type de Schrödinger symétrisations indépendendantes du temps, C. R. Acad. Sci. Paris, t315 serie 1 (1992), 1055-1058.
  29. L. Vega, Schrödinger equations: pointwise convergence to the initial data, Proc A.M.S. 102 (1988), 874–878.
  30. V. E. Zakharov, E. I. Schulman, Degenerated dispersion laws, motion invariant and kinetic equations, Physica 1D (1980), 185-250.

Carlos E. Kenig Department of Mathematics University of Chicago Chicago, Il. 60637 USA E-mail: cek@math.uchicago.edu Gustavo Ponce Department of Mathematics University of California Santa Barbara, CA 93106 USA E-mail: ponce@math.ucsb.edu Christian Rolvung Nykredit Markets & Asset Management Kalvevod Bridge 1-3 DK-1780 Copenhagen V Denmark E-mail: cro@nykredit.dk Luis Vega Departamento de Matematicas Universidad del Pais Vasco Apartado 644 48080 Bilbao Spain E-mail: mtpvegol@lg.ehu.es