Variable coefficient Schrödinger flows for ultrahyperbolic operators

C.E. KenigPartially supported by NSF and IBERDROLA program of Profesores Visitantes University of Chicago, Chicago Il., USA 60637, cek@math.uchicago.edu

G. PoncePartially supported by NSF and IBERDROLA program of Profesores Visitantes University of California, Santa Barbara, Ca, USA ponce@math.ucsb.edu

C. Rolvung Nykredit Markets & Asset Management, Kalvevod Bridge 1–3, DK-1780 Copenhagen V, Denmark, cro@nykredit.dk

L. VegaPartially supported by a MECyT grant Universidad del Pais Vasco, Apdo. 644, 48080, Bilbao, Spain, mtpvegol@lg.ehu.es

Contents

1 INTRODUCTION

In this paper we shall consider nonlinear Schrödinger equations of the form
t u = i ( x ) u + b 1 ( x ) x u + b 2 ( x ) x u ¯ + c 1 ( x ) u + c 2 ( x ) u ¯ + P ( u , u ¯ , x u , x u ¯ ) , (1.1)
where x R n   , t > 0   , ( x ) = j , k = 1 n x j ( a j k ( x ) x k )   , A ( x ) = ( a j k ( x ) ) j , k = 1 , . . , n   is a real, symmetric and nondegenerate variable coefficient matrix, and P   is a polynomial with no linear or constant terms.
Equations of the form described in ( 1.1 ) with A ( x )   merely invertible as opposed to positive definite arise in connection with water wave problems, and in higher dimensions as completely integrable models, see [5, [14, and [32.
In this work we shall study the existence, uniqueness and regularity of the local solutions to the initial value problem (IVP) associated to the equation ( 1.1 ). The class of equations is rather general and appropriate assumptions have to be imposed on the smoothness and decay of the coefficients a j k , b 1 , b 2 , c 1   and c 2   and on the initial data u ( x , 0 ) = u 0 ( x )   as well as on the asymptotic behavior of a j k ( x )   as | x |   . Also it will be necessary to measure the regularity of solutions in weighted Sobolev spaces of high indexes. The main result we obtain in this direction is Theorem  6.2.1 -also see Remark  6.2.2 , in Section  6 .
One of the main difficulties in ( 1.1 ) is that the nonlinear terms incur in the so called “loss of derivatives”. This can be avoided if P   is assumed to have a special symmetric form and b 1   is real valued. In this case, the standard energy method gives local well-posedness of the corresponding IVP in H s ( R n )   for s > n / 2 + 1   independently of the dispersive nature of ( 1.1 ), see [15, [23. Another approach used to overcome this loss of derivatives is to restrict oneself to working with = Δ   , b 1 = b 2 = 0   in suitable analytic function spaces, see [10and references therein.
In [18, C. E. Kenig, G. Ponce and L. Vega used linear dispersive smoothing effects of the associated linear equation to show that the IVP for the equation ( 1.1 ) with = Δ   , b 1 = b 2 = 0   , c 1 = c 2 = 0   and general P   is locally well posed in (possibly weighted) Sobolev spaces with high index for small initial data. For the case n = 1   , N. Hayashi and T. Ozawa [11removed the smallness condition by using an integrating factor which reduces the problem to a system where the energy method applies. H. Chihara [2removed the smallness assumptions in weighted Sobolev spaces [18in any dimension n   by considering systems of two equations which he diagonalized to essentially eliminate the conjugate first order terms. The remaining first order terms are treated by a method similar to the one used by S. Mizohata [26and S. Doi [7to solve linear Schrödinger equations with lower order terms. It consists in applying a pseudo-differential operator K   to the equation. The commutator i [ K Δ Δ K ]   basically absorbs the first orderterm to overcome the loss of derivatives in a way related to the method of integrating factors. In Chihara's approach the ellipticity of = Δ   is key in the diagonalization argument. C. E. Kenig, G. Ponce and L. Vega [22obtained local well-posedness for the IVP for ( 1.1 ) in the non-elliptic constant coefficient case = j = 1 k x j 2 j = k + 1 n x j 2 , k = 1 , . . , n 1 ,   using the pseudo-differential operators of [1in the linear problem to avoid the diagonalization process. Furthermore, their results are valid in Sobolev spaces with no weights if P   has no quadratic terms.
It may be gathered from this short summary of background literature that a thorough understanding of linear Schrödinger equations is important in the attempt to solve the nonlinear problem for ( 1.1 ).
Our approach in this work will be illustrated with the special case of ( 1.1 )
{ t u = i ( x ) u + b 1 ( x ) x u + u x 1 u , u ( x , 0 ) = u 0 ( x ) , (1.2)
or equivalently
{ t u = i ( x ) u + [ b 1 ( x ) x u + u 0 ( x ) x 1 u ] + ( u u 0 ) x 1 u , u ( x , 0 ) = u 0 ( x ) . (1.3)
The nonlinear part ( u u 0 ) x 1 u   of ( 1.3 ) should be small for small t   because of the factor ( u u 0 )   , but the factor x 1 u   still incurs loss of one derivative. The linear part of ( 1.3 ) has a modified first order coefficient.
It is therefore useful to study linear Schrödinger equations of the form
{ t u = i ( x ) u + b 1 ( x ) x u + b 2 ( x ) x u ¯ + c 1 ( x ) u + c 2 ( x ) u ¯ + f ( x , t ) , u ( x , 0 ) = u 0 ( x ) . (1.4)
Solutions u   of ( 1.4 ) gain one derivative compared to f   and 1 / 2   derivative compared to u 0   , on the average in time and modulo spatial weights, under suitable assumptions. More precisely, for s Z +   and N Z +   , N > 1   , the solution of ( 1.4 ) satisfies
0 T R n | J s + 1 / 2 u ( x , t ) | x N d x d t c ( ( 1 + T ) u 0 H s 2 + 0 T R n | J s 1 / 2 f ( x , t ) | x N d x d t ) (1.5)
For the proof of these results for the constant coefficient case ( x ) = Δ   , b 1 = b 2 = 0   see [16, [24, [3, [29, [31, [17, [18). The estimate ( 1.5 ) allows to overcome the loss of one derivative introduced by the nonlinear part of ( 1.3 ) and, more generally, to solve ( 1.1 ).
Consider ( 1.4 ) with ( x ) = Δ   , b 1 = ( i , 0 , . . , 0 )   , b 2 = c 1 = c 2 = f = 0   . The solution of this constant coefficient equation is given via Fourier transform by u ^ ( ξ , t ) = exp ( t ( i | ξ | 2 + ξ 1 ) ) u ^ 0 ( ξ )   .
The multiplier exp ( t ( i | ξ | 2 + ξ 1 ) )   is unbounded for t 0   , so ( 1.4 ) is not wellposed in L 2   in this case. In fact, the following condition, deduced by S. Mizohata [26, has been proven to be necessary for the well-posedness in L 2   of ( 1.4 ) with ( x ) = Δ   , b 2 = 0  
sup x R n , ω S n 1 | 0 Im b 1 ( x + r ω ) ω d r | < . (1.6)
So the decay assumptions on Im b 1   are natural. An application of J s = ( I Δ ) s / 2   to ( 1.4 ) gives a new equation with x a j k ( x )   appearing in the first order coefficient. Well posedness of ( 1.4 ) is of interest for any s   , so decay assumptions on x a j k ( x )   seem also natural.
To justify the decay assumptions on b 2 ( x )   , we use the result in [20. Consider the IVP
{ t u = i x 1 x 2 2 u + b 2 ( x ) x u ¯ , x R 2 , t > 0 , u ( x , 0 ) = u 0 ( x ) (1.7)
with b 2 = ( i , 0 )   . It was shown in [20that the IVP ( 1.7 ) is wellposed, however its solutions gain only 1 / 4   of derivative compared with u 0   instead of the expected 1 / 2   , i.e. ( 1.5 ) holds with f = 0   and 1 / 4   instead of 1 / 2   . Moreover, if we replace x 1 x 2 2   by x 1 2 + x 2 2   the expected gain of 1 / 2   derivatives is obtained. So one has that in the non-elliptic case, decay assumptions on b 2   are also necessary to obtain ( 1.5 ). The main result we obtain concerning ( 1.4 ) is Theorem  5.1.1 in Subsection  5.1 .
When considering the variable coefficient equation in ( 1.1 ) one should study of the bicharacteristic flow, i.e. solutions ( X ( s ; x , ξ ) , Ξ ( s ; x , ξ ) )   of the system
{ d d s X j ( s ; x , ξ ) = 2 k = 1 n a j k ( X ( s ; x , ξ ) ) Ξ k ( s ; x , ξ ) , d d s Ξ j ( s ; x , ξ ) = k , l = 1 n x j a k l ( X ( s ; x , ξ ) ) Ξ k ( s ; x , ξ ) Ξ l ( s ; x , ξ ) , ( X ( 0 ; x , ξ ) , Ξ ( 0 ; x , ξ ) ) = ( x , ξ ) . (1.8)
In the constant coefficient case, = a j k 0 x j x k   , one has that the bicharacteristic flow is ( X ( s ; x , ξ ) , Ξ ( s ; x , ξ ) ) = ( x + 2 s A 0 ξ , ξ ) , A 0 = ( a j k 0 ) j , k = 1 , . . , n .   For = Δ   one gets ( X ( s ; x , ξ ) , Ξ ( s ; x , ξ ) ) = ( x + 2 s ξ , ξ )   and the condition in ( 1.6 ) can be seen as an integrability one along the bicharacteristics. As we will see in this work (Sections 4-5), roughly speaking the operator K   whose symbol is
k ( x , ξ ) = exp ( 0 b ( X ( s ; x , ξ ) , Ξ ( s ; x , ξ ) ) d s ) , with b ( x , ξ ) = Im b 1 ( x ) ξ , (1.9)
will play the role of the “integrating factor” introduced in [11. Such constructions were also previously used in the works [4and [8. The commutator term i [ K K ]   , used to cancel the term b 1 ( x )   , corresponds to differentiation of K   along the bicharacteristic flow. Unfortunately the symbol in ( 1.9 ) is not in a standard class. It satisfies
| x α ξ β k ( x , ξ ) | c α β x | β | ξ | β | . (1.10)
We observe that in the particular constant coefficient elliptic case, = Δ   we have
k ( x , ξ ) = exp ( 0 Im b 1 ( x + s ξ ) ξ d s )
= exp ( 0 Im b 1 ( x + s ξ ^ ) ξ ^ d s ) , ξ ^ = ξ / | ξ | ,
which is related to Mizohata's condition in ( 1.6 ).
In the case where ( x )   is elliptic, the class described in ( 1.10 ) was introduced and studied by W. Craig, T. Kappeler and W. Strauss in [4. However, it should be pointed out that in the non-elliptic case, i.e. ( x )   is just nondegenerate, the geometric assumption (4.2) in [4is not satisfied by symbols of interest -see subsection  3.1 in Section  3 . Moreover, we observe that the Hamiltonian h 2 ( x , ξ ) = j , k a j k ( x ) ξ k ξ j   is preserved under the flow. Then one of the main differences of the flows considered here with those associated to elliptic operators is that ellipticity gives the a priori estimate ν 2 | ξ 0 | 2 | Ξ ( s ; x 0 , ξ 0 ) | 2 ν 2 | ξ 0 | 2 ,   which guarantees that the solutions of the system ( 1.8 ) are globally defined.
We will assume that the bicharacteristic flow is non-trapping, i.e. for each ( x 0 , ξ 0 ) R n × ( R n { 0 } )   and for each μ > 0   there exists s 0 > 0   such that
| X ( s 0 ; x 0 , ξ 0 ) | μ . (1.11)
The non-trapping condition appears naturally, since Ichinose [13showed that a necessary condition for the well-posedness in L 2   of ( 1.1 ) with   elliptic, b 2 0   , c 1 0 ,   c 2 0 ,   and P 0 ,   is that the analog of ( 1.6 ) must hold in this case, with the integration taking place along the bicharacterisitics.
The non-trapping condition also is essential in the works of [4and [8. In fact even when also b 1 0   and f 0   , Doi showed in [9that it is necessary for ( 1.5 ) to hold.
In the ultra-hyperbolic case (i.e. with a merely non-degenerate matrix A   ), under appropriate decay assumptions and asymptotic behavior as | x |   on the coefficients a j k ( x )   , we shall prove that the bicharacteristic flow is globally defined and “uniformly non-trapping”. Moreover in order to keep the structure of the conjugate first order terms, so that after applying the operator K   we can obtain energy estimates, we need the symbol of K   to be even -see Definition  5.2.1 (iv) in Section  5 . Therefore we have to study carefully the bicharacteristic for backward and forward time. In particular when looking at the forward bicharacteristic the more delicate part is when it is not outgoing -see Theorem  4.1.1 of Section  4 . In that region we prove that outside a bounded ball, in the x   variable, it behaves in dyadic annuli as the free flow -see Theorem  4.1.1 in Section  4 for a precise statement.
As in [18, [22, the proof of the nonlinear results relies on two kinds of linear estimates. The first one is concerned with the smoothing effect described in ( 1.5 ) for solutions of the IVP ( 1.4 )
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 ) H s 2 + c 0 T R n | f ( x , t ) | 2 d x d t , (1.12)
and
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 ) H s 2 + 0 T R n | J s 1 / 2 f ( x , t ) | 2 x N d x d t
for N > 1   .
The second kind is related with the local well-posedness in L 2   (and in H s   ) of the IVP ( 1.4 ). To establish this result we follow an indirect approach. First we truncate at infinity the operator ( x )   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 ,   where 0 = a j k 0 x j x k   , A 0 = ( a j k 0 ) j , k = 1 , . . , n   is a (constant) matrix, with the decay assumption a j k ( x ) a j k 0 S ( R n )   , j , k = 1 , . . , n   , (although we will work in the S ( R n )   class, it will be clear from our proofs that the same results hold if we just assume that the corresponding estimate holds for a finite number of seminorms in ( 2.2 ) Section 2). Thus, ( x ) = R ( x ) + R ( x ) .   For R   large enough we consider 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.9 ) but evaluated in the bicharacteristic flow ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) )   . To obtain the L 2   local well-posedness of the IVP ( 1.4 ) we show that there exists N 0   depending only on the dimension such that for any M Z +   there exists R 0 = R 0 ( M )   such that for R R 0  
sup 0 t T K R u ( t ) L 2 2 c R N 0 u ( 0 ) L 2 2
+ R M 0 T R n | J 1 / 2 u | 2 x N d x d t + c T R N 0 sup 0 t T u ( t ) L 2 2 .
Next, we deduce several estimates concerning the operator K R   . In particular, for 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   , which allows us to treat E R u ( t )   as an error term -see Lemma  5.2.6 in Section  5 . Collecting these results we get that
sup 0 t T u ( t L 2 2 R N 0 u ( 0 ) L 2 2
+ R M 0 T R n | J 1 / 2 u | 2 x N d x d t + c T sup 0 t T u ( t ) L 2 2 ,
which combined with ( 1.12 ) yields the desired estimate, i.e. the local well-posedness in L 2   of the IVP ( 1.4 ) for T   sufficiently small sup 0 t T u ( t ) L 2 2 c ( T ) ( u 0 L 2 2 + 0 T f ( t ) L 2 2 d t ) .   The smoothing effect and local well-posedness in H s   of ( 1.4 ) in the case where   has elliptic variable coefficients will be proven in Section  2 -see Lemma  2.2.2 and Theorem  2.3.1 of that section. This builds on S. Doi's pseudo-differential method in [7, [8and on H. Chihara's diagonalization method for systems [2and uses only classical pseudo-differential operators.
The diagonalization method cannot be used when   is ultrahyperbolic. When   has constant coefficients it is possible to cancel the loss of derivatives using a pseudo-differential transformation which falls under the scope of Calderón-Vaillancourt's theorem -see [22, but this does not seem to extend to the variable coefficient case. Instead one is led to study a new class of symbols and this is done in Section  3 . As we already mentioned the corresponding operators in the elliptic case were studied in [4.
A typical example of the symbols that need to be considered (take n = 2   for simplicity) is a ( x ( ξ 2 , ξ 1 ) / | ξ | ) χ ( ξ ) ,   where χ C   , χ 1   for | ξ | 2   , χ 0   for | ξ | < 1   , and a C 0 ( R )   , in the elliptic case; and a ( x ( ξ 2 , ξ 1 ) / | ξ | ) χ ( ξ ) ,   in the ultrahyperbolic one with χ   and a   as before.
As is explained in [21, the operators of the elliptic case are easily reduced to classical pseudo-differential operators by expressing a   in terms of its Fourier transform and, given τ R ,   using the invertiblechange of variable ξ ξ + τ ( ξ 2 , ξ 1 ) / | ξ |   . In the ultrahyperbolic setting this approach fails since the corresponding mappings are not invertible, and hence the theory, in particular the L 2   boundedness, is more delicate -see Theorem  3.2.1 in Section  3 . The proofs of the rest of the results concerning the calculus of the operators arising from these symbols-Theorem  3.3.1 , Theorem  3.3.2 and Theorem  3.3.3 of the same section, are reduced after some manipulations to the L 2   boundedness.
In Section  4 we study the bicharacteristic flow in the ultrahyperbolic case for ( x )   and its truncated version R ( x )   . There we shall deduce several estimates to be used in establishing the smoothing effect and the local well-posedness of ( 1.4 ) with ( x )   non-elliptic which will be given in Section  5 .
This also relies on the calculus of Section  3 .
Finally, the smoothing effect in ( 1.4 ) is used to solve ( 1.1 ) in Section  6 . Solutions of ( 1.1 ) are fixed points of an integral mapping which is a contraction on a suitable function space in a small time interval, so Banach's contraction mapping principle applies.
The results in the elliptic case, i.e. those in Section  2 , are due to C. Rolvung, and appear in his PhD dissertation [27. The results in the ultrahyperbolic case, for   a C 0   perturbation of a constant coefficient operator 0   also appear in [27.

2 THE LINEAR ELLIPTIC EQUATION

The local well posedness and smoothing effect for linear elliptic equations are considered in this section. This builds on S. Doi`s method involving classical pseudo-differential operators and the sharp Gårding inequality [8as well as on a diagonalization as in [2.

2.1 Pseudo-differential Operators

First we will recall some results from the theory of pseudo-differential operators.
The class S m = S 1 , 0 m   of classical symbols of order m R   is defined by
S m = { p ( x , ξ ) C ( R n × R n ) : | p | S m ( j ) < , j N } (2.1)
where
| p | S m ( j ) = sup { ξ m + | α | ξ α x β p ( , ) L ( R n × R n ) : | α + β | j } (2.2)
and ξ = ( 1 + | ξ | 2 ) 1 / 2   .
The pseudo-differential operators Ψ p   associated to the symbol p S m   is defined by
Ψ p f ( x ) = R n e i x ξ p ( x , ξ ) f ^ ( ξ ) d ξ ( 2 π ) n / 2 , f S ( R n ) . (2.3)
For example, a partial differential operator P = | α | N a α ( x ) x α ,   with a α C b ( R n )   is a pseudo-differential operator P = Ψ p   with symbol p ( x , ξ ) = | α | N a α ( x ) ( i ξ ) α .   The fractional differentiation operator J s = Ψ ξ s   is also a pseudo-differential operator. The collection of symbol classes S m , m R   , is in some cases closed under the division and square root operations. This is not the case for polynomials in ξ   and sometimes allows one to construct approximate inverses and square roots of pseudo-differential operators.
The following facts will be used throughout this work and the proofs can be found for example in [25.
Theorem 2.1.1 (Sobolev boundedness) Let m R   , p S m   and s R   . Then Ψ p   extends to a bounded linear operator from H m + s ( R n )   to H s ( R n )   . Moreover, there exist j = j ( n ; m ; s ) N   and c = c ( n ; m ; s )   such that
Ψ p f H s c | p | S m ( j ) f H m + s . (2.4)
Finally for p S 0   and λ ( | x | ) = ( 1 + | x | 2 ) N / 2 = x N   , N > 1   there exists j = j ( n , N )   such that
R n | Ψ p f ( x ) | 2 λ ( | x | ) d x c | p | S 0 ( j ) R n | f ( x ) | 2 λ ( | x | ) d x . (2.5)
The proof of ( 2.5 ) can be seen for example in [22.
Theorem 2.1.2 (Symbolic calculus) Let m 1 , m 2 R   , p 1 S m 1   , p 2 S m 2   . Then there exist p 3 S m 1 + m 2 1   , p 4 S m 1 + m 2 2   , and p 5 S m 1 1   such that
Ψ p 1 Ψ p 2 = Ψ p 1 p 2 + Ψ p 3 , Ψ p 1 Ψ p 2 Ψ p 2 Ψ p 1 = Ψ i { p 1 , p 2 } + Ψ p 4 , ( Ψ p 1 ) * = Ψ p ¯ 1 + Ψ p 5 (2.6)
where { p 1 , p 2 }   denotes the Poisson bracket, i.e.
{ p 1 , p 2 } = j = 1 n ( ξ j p 1 x j p 2 x j p 1 ξ j p 2 ) ,   and such that for any j N   there exist j N   and c 1 = c 1 ( n ; m 1 ; m 2 ; j )   , c 2 = c 2 ( n ; m 1 ; j )   such that
| p 3 | S m 1 + m 2 1 ( j ) + | p 4 | S m 1 + m 2 2 ( j ) c 1 | p 1 | S m 1 ( j ) | p 2 | S m 2 ( j ) | p 5 | S m 1 1 ( j ) c 2 | p 1 | S m 1 ( j ) . (2.7)
Theorem 2.1.3 (Sharp Gårding inequality) Let p S 1   and suppose that there exists R 0   such that Re p ( x , ξ ) 0   for ξ R   . Then there exist j = j ( n )   and c = c ( n ; R )   such that
Re Ψ p f ; f L 2 c | p | S 1 ( j ) f L 2 2 , f S ( R n ) . (2.8)
This result is due to L. Hörmander [12.

2.2 The Bicharacteristic Flow

The basic idea is to apply a pseudo-differential operator K   to the equation ( 1.4 ) of Section 1 in such a way that the commutator i [ K K ]   cancels b 1 ( x ) x   . It turns out that i [ K K ]   corresponds to differentiation along the bicharacteristic flow which will now be introduced.
Let A ( x ) = ( a j k ( x ) )   be a real, and symmetric n × n   matrix of functions a j k C b   . We will assume that
| a j k ( x ) | = o ( | x | 1 ) as | x | , j , k = 1 , . . , n , (2.9)
and that A ( x )   is positive definite, i.e.
ν > 0 x , ξ R n ν 1 | ξ | 2 | j , k = 1 n a j k ( x ) ξ j ξ k | ν | ξ | 2 . (2.10)
Let h 2   be the principal symbol of = x j a j k ( x ) x k   , i.e.
h 2 ( x , ξ ) = j , k = 1 n a j k ( x ) ξ j ξ k . (2.11)
The bicharacteristic flow is the flow of the Hamiltonian vector field
H h 2 = j = 1 n [ ξ j h 2 x j x j h 2 ξ j ] (2.12)
and is denoted by ( X ( s ; x 0 , ξ 0 ) , Ξ ( s ; x 0 , ξ 0 ) )   , i.e.
{ d d s X j ( s ; x 0 , ξ 0 ) = 2 k = 1 n a j k ( X ( s ; x 0 , ξ 0 ) ) Ξ k ( s ; x 0 , ξ 0 ) , d d s Ξ j ( s ; x 0 , ξ 0 ) = k , l = 1 n x j a l k ( X ( s ; x 0 , ξ 0 ) ) Ξ k ( s ; x 0 , ξ 0 ) Ξ l ( s ; x 0 , ξ 0 ) (2.13)
for j = 1 , . . , n   , with
( X ( 0 ; x 0 , ξ 0 ) , Ξ ( 0 ; x 0 , ξ 0 ) ) = ( x 0 , ξ 0 ) . (2.14)
The bicharacteristic flow exists in the time interval s ( δ , δ )   with δ = δ ( x 0 , ξ 0 )   , and δ ( )   depending continuously on ( x 0 , ξ 0 )   .
The bicharacteristic flow preserves h 2   , so ellipticity gives
ν 2 | ξ 0 | 2 | Ξ ( s ; x 0 , ξ 0 | 2 ν 2 | ξ 0 | 2 , (2.15)
and hence δ =   .
It will be assumed that the bicharacteristic flow is non-trapped which means that the set { X ( s ; x 0 , ξ 0 ) : s R }   is unbounded in R n   for each ( x 0 , ξ 0 ) R n × R n { 0 }   .
Note that h 2   is homogeneous of degree 2   in ξ   so that
{ X ( s ; x , t ξ ) = X ( t s ; x , ξ ) , Ξ ( s ; x , t ξ ) = t Ξ ( t s ; x , ξ ) . (2.16)
The next result shows that the Hamiltonian vector field is differentiation along the bicharacteristics.
Lemma 2.2.1 Let φ C ( R n × R n )   . Then
( H h 2 φ ) ( x , ξ ) = s [ φ ( X ( s ; x , ξ ) , Ξ ( s ; x , ξ ) ) ] | s = 0 . (2.17)
The following key lemma is due to S. Doi [8(Lemmas 2.3-2.5).
Lemma 2.2.2 Let A ( x )   and its bicharacteristic flow satisfy the assumptions above. Suppose λ L 1 ( [ 0 , ) ) C ( [ 0 , ) )   is strictly positive and nonincreasing. Then there exist c > 0   and a real symbol p S 0   , both depending on h 2   and λ   , such that
H h 2 p = { h 2 , p } ( x , ξ ) λ ( | x | ) | ξ | c , ( x , ξ ) R n × R n . (2.18)
An extension of this result to the case of invertible A ( x )   will be given in Section  5 , Lemma  5.1.1 .

2.3 Linear Elliptic Smoothing Effects

In this subsection we consider the IVP associated to the linear Schrödinger equation
{ t u = i j , k = 1 n x j ( a j k ( x ) x k u ) + b 1 ( x ) x u + b 2 ( x ) x u ¯ + c 1 ( x ) u + c 2 ( x ) u ¯ + f ( x , t ) , u ( x , 0 ) = u 0 ( x ) , (2.19)
where A ( x ) = ( a j k ) j , k = 1 , . . , n   satisfies ( 2.9 )-( 2.10 ) and its bicharacteristic flow satisfies the assumptions in the previous subsection, b l = ( b l 1 , . . , b l n ) ( C b ) n   , l = 1 , 2   and c 1 , c 2 C b   .
Combining the equation in ( 2.9 ) and its complex conjugate we obtain a system in w = ( u , u ¯ ) T  
{ t w = ( i H + B + C ) w + f , w ( x , 0 ) = w 0 ( x ) , (2.20)
where H =
( 0 0 ) , C = ( c 11 c 12 c 21 c 22 ) ,
B = ( B 11 B 12 B 21 B 22 ) ,   with = j , k = 1 n x j ( a j k ( x ) x k ) ,   B l m u ( x , t ) = j = 1 n b l m j ( x ) x j u ( x , t ) , b l m j C b l , m = 1 , 2 , j = 1 , . . , n ,   and c l m = c l m ( x ) C b , l , m = 1 , 2 , f ( x , t ) = ( f ( x , t ) , f ¯ ( x , t ) ) T .   The following well-posedness and smoothing results contain three parts depending on the regularity and the decay of the external force f ( x , t )   .
Theorem 2.3.1 Let w 0 = ( u 0 , u ¯ 0 ) T ( H s ( R n ) ) 2   , s R   . Assume that there exist N > 1   and a constant c 0   such that if λ ( | x | ) = x N   then | Im b l l j ( x ) | c 0 λ ( | x | ) , l , m = 1 , 2 , j = 1 , . . , n , x R n .   Then (a) If f ( L 1 ( [ 0 , T ] : ( H s ( R n ) ) ) 2   then the IVP ( 2.20 ) has a unique solution w C ( [ 0 , T ] : ( H s ( R n ) ) 2 )   satisfying
sup 0 t T w ( t ) ( H s ) 2 c 1 e c 2 T ( w 0 ( H s ) 2 + 0 T f ( t ) ( H s ) 2 d t ) . (2.21)
(b) If f ( L 2 ( [ 0 , T ] : ( H s ( R n ) ) ) 2   then the IVP ( 2.20 ) has a unique solution w C ( [ 0 , T ] : ( H s ( R n ) ) 2 )   satisfying
sup 0 t T w ( t ) ( H s ) 2 2 + 0 T R n | J s + 1 / 2 w ( x , t ) | 2 λ ( | x | ) d x d t c 1 e c 2 T ( w 0 ( H s ) 2 2 + 0 T f ( t ) ( H s ) 2 2 d t ) . (2.22)
(c) If J s 1 / 2 f ( L 2 ( R n × [ 0 , T ] : ( λ ( | x | ) 1 d x d t ) ) 2   then the IVP ( 2.20 ) has a unique solution w C ( [ 0 , T ] : ( H s ( R n ) ) 2 )   satisfying
sup 0 t T w ( t ) ( H s ) 2 2 + 0 T R n | J s + 1 / 2 w ( x , t ) | 2 λ ( | x | ) d x d t c 1 e c 2 T ( w 0 ( H s ) 2 2 + 0 T R n | J s 1 / 2 f ( x , t ) | 2 ( λ ( | x | ) ) 1 d x d t ) . (2.23)
Here c 1   depends on n , s , ν , ( a j k ) j , k = 1 , . . , n , c 0 , ( b l m j ) l , m = 1 , 2 ; j = 1 , . . , n   , and c 2   depends in addition on ( c l m ) l , m = 1 , 2   .
Corollary 2.3.1 Let s R   and u 0 H s ( R n   ) and suppose λ   satisfies the assumptions of Theorem  2.3.1 . Then (a) If f L 1 ( [ 0 , T ] : H s ( R n ) )   , then there exists a unique solution u C ( [ 0 , T ] : H s ( R n ) )   of ( 2.19 ) satisfying sup 0 t T u ( , t ) H s c 1 e c 2 T ( u 0 H s + 0 T f ( , t ) H s d t )   (b) If f L 2 ( [ 0 , T ] : H s ( R n ) )   , then there exists a unique solution u C ( [ 0 , T ] : H s ( R n ) )   of ( 2.19 ) satisfying
sup 0 t T u ( , t ) H s 2 + 0 T R n | J s + 1 / 2 u ( x , t ) | 2 λ ( | x | ) d x d t
c 1 e c 2 T ( u 0 H s 2 + 0 T f ( , t ) H s 2 d t ) ;
(c) If J s 1 / 2 f L 2 ( R n × [ 0 , T ] : ( λ ( | x | ) ) 1 d x d t )   , then there exists a unique solution u C ( [ 0 , T ] : H s ( R n ) )   of ( 2.19 ) satisfying
sup 0 t T u ( , t ) H s 2 + 0 T R n | J s + 1 / 2 u ( x , t ) | 2 λ ( | x | ) d x d t
c 1 e c 2 T ( u 0 H s 2 + 0 T R n | J s 1 / 2 f ( x , t ) | 2 ( λ ( | x | ) ) 1 d x d t ) .
Here c 1   depends on n , s , ν , ( a j k ) j , k = 1 , . . , n , c 0 , ( b lmj ) l , m = 1 , 2 ; j = 1 , . . , n   , and c 2   depends in addition on ( c l m ) l , m = 1 , 2   .
The following a priori estimate is needed for the proof of Theorem  2.3.1 .
Lemma 2.3.1 Let s R   and suppose λ   satisfies the assumptions of Theorem  2.3.1 . Then there exist c 1   depending on n , s   , ν , ( a j k ) j , k = 1 , . . , n   , c 0 ,   and finitely many derivatives of ( b l m j ) l , m = 1 , 2 ; j = 1 , . . , n   , and c 2   depending in addition on finitely many derivatives of ( c l m ) l , m = 1 , 2   such that for all w ( C ( [ 0 , T ] : H s + 2 ( R n ) ) C 1 ( [ 0 , T ] : H s + 2 ( R n ) ) ) 2   the following four estimates hold :
( i ) sup 0 t T w ( t ) ( H s ) 2 c 1 e c 2 T ( w 0 ( H s ) 2
+ 0 T ( t ( i H + B + C ) ) w ( , t ) ( H s ) 2 d t ) ,
( i i ) sup 0 t T w ( t ) ( H s ) 2 c 1 e c 2 T ( w ( , T ) ( H s ) 2
+ 0 T ( t + ( i H + B + C ) * ) w ( , t ) ( H s ) 2 d t ) ,
( i i i ) sup 0 t T w ( t ) ( H s ) 2 2 + 0 T R n | J s + 1 / 2 w ( x , t ) | 2 λ ( | x | ) d x d t
c 1 e c 2 T ( w 0 ( H s ) 2 2 + 0 T ( t ( i H + B + C ) ) w ( , t ) ( H s ) 2 2 d t ) ,
( i v ) sup 0 t T w ( t ) ( H s ) 2 2 + 0 T R n | J s + 1 / 2 w ( x , t ) | 2 λ ( | x | ) d x d t
c 1 e c 2 T w 0 ( H s ) 2 2
+ c 1 e c 2 T 0 T R n | J s 1 / 2 ( t ( i H + B + C ) ) w ( , t ) | 2 ( λ ( | x | ) ) 1 d x d t .
We observe that (ii) follows from (i) by applying (i), with i H + B + C   replaced by ( i H + B + C ) *   , to w ( , T t )   , so it suffices to prove (i), (iii), and (iv) of Lemma  2.3.1 . The idea of the proof is to apply transformations Λ   and K   to the system. Λ   will diagonalize B   and essentially transforms the system into two single equations where the pseudo-differential calculus applies. The idea of this diagonalization came from the work of H. Chihara [2. K   will eliminate the loss of derivatives of the first order terms. This idea is due to S. Doi [7, [8and S. Mizohata [26.
Proof of Lemma  2.3.1 Let w ( C ( [ 0 , T ] : H s + 2 ( R n ) ) C 1 ( [ 0 , T ] : H s ( R n ) ) ) 2   . Set h 2 ( x , ξ ) = j , k = 1 n a j k ( x ) ξ j ξ k , h 1 ( x , ξ ) = i j , k = 1 n x j a j k ( x ) ξ k ,   so that = Ψ h 1 + h 2   . Since ( a j k ( x ) )   is positive definite there exist c = c ( a j k ; x j a j k )   and R = R ( a j k ; x j a j k )   such that | h 1 ( x , ξ ) + h 2 ( x , ξ ) | c | ξ | 2 , | ξ | R .   Choosing φ C 0 ( R n )   with φ ( y ) = 1   if | y | R   and φ ( y ) = 0   if | y | 2 R   we define h ~ ( x , ξ ) = ( h 1 ( x , ξ ) + h 2 ( x , ξ ) ) 1 ( 1 φ ( ξ ) ) and ~ = Ψ h ~ .   So h ~ S 2   and ~ = I + Ψ r 1 ,   where r 1 = r 1 ( ( a j k ) j , k = 1 , . . , n ) S 1   by the symbolic calculus in Theorem  2.1.2 in the sense that for any l N   , | r 1 | S 1 l   depends on the ellipticity constant μ   and on finitely many derivatives of the a j k   's.
We define B d i a g   , B a n t i   with symbols in ( S 1 ) 2 × 2 , S 12 , S 21   with symbols in S 1   and S   with symbol in ( S 1 ) 2 × 2   by
B d i a g = ( B 11 0 0 B 22 ) , B a n t i = ( 0 B 12 B 21 0 ) ,
S 12 = 1 2 i B 12 ~ , S 21 = 1 2 i B 21 ~ ,   S = ( 0 S 12 S 21 0 ) ,   and the diagonalizing transform Λ   with symbol in ( S 0 ) 2 × 2   by Λ = I S .   Then S = S ( n ; ν ; ( a j k ) ; ( b l m j ) )   in the sense above. Letting f = ( t ( i H + B + C ) ) w ,   and applying Λ   , one obtains
t Λ w = i Λ H w + Λ B w + Λ C w + Λ f . (2.24)
We shall show that the system ( 2.24 ) is diagonalized modulo operators with symbols in S 0   . So wewrite
i Λ H + Λ B = i H Λ + ( i Λ H i H Λ ) + B S B = i H Λ + ( i Λ H i H Λ ) + B d i a g + B a n t i S B = ( i H Λ + B d i a g Λ ) + ( B a n t i + i Λ H i H Λ ) + ( B d i a g S S B ) . (2.25)
Since the operator in the first parenthesis is diagonalized and the operator in the last parenthesis has order 0   , it suffices to consider only the operator in the second parenthesis in ( 2.25 ). Thus,
B a n t i + i Λ H i H Λ = B a n t i + i H S i S H =
= ( 0 B 12 B 21 0 ) + ( 0 0 ) ( 0 1 2 B 12 ~ 1 2 B 21 ~ 0 )
( 0 1 2 B 12 ~ 1 2 B 21 ~ 0 ) ( 0 0 )
= ( 0 B 12 1 2 B 12 ~ 1 2 B 12 ~ B 21 1 2 B 21 ~ 1 2 B 21 ~ 0 ) .
We observe that B 12 ~ = B 12 ~ + Ψ r 2   , where r 2 S 0   and B 12 1 2 B 12 ~ 1 2 B 12 ~ = B 12 B 12 ~ 1 2 Ψ r 2 = B 12 Ψ r 1 1 2 Ψ r 2 .   A similar calculus argument handles the term involving B 21   . Therefore, we have that B a n t i + i Λ H i H Λ   has order zero, which allows to conclude that
t Λ w = i H Λ w + B d i a g Λ w + Ψ r 3 w + Λ f , (2.26)
where r 3 ( S 0 ) 2 × 2   and r 3 = r 3 ( n ; ν ; ( a j k ) ; ( b l m j ) ; ( c l m ) ) ,   in the sense that for any j 0 Z +   , | r 3 | ( S 0 ) 2 × 2 ( j 0 )   depends on n , ν   and finitely many of the derivatives of a j k   , b l m j   and c l m   .
By Lemma  2.2.2 there exists a real-valued p S 0   and C > 0   , both depending on ( a j k )   and c 0   , such that
{ h 2 ( x , ξ ) ; p ( x , ξ ) } C c 0 λ ( | x | ) | ξ | C , x , ξ R n , (2.27)
with C = C ( n )   to be determined. Let k ( x , ξ ) = ( e p ( x , ξ ) ( 1 + ( h 2 ( x , ξ ) ) 2 ) s / 4 0 0 e p ( x , ξ ) ( 1 + ( h 2 ( x , ξ ) ) 2 ) s / 4 )   and K = Ψ k   . Note that e ± p ( x , ξ ) ( 1 + ( h 2 ( x , ξ ) ) 2 ) s / 4 S s   and is elliptic, since p S 0   is real and ( 1 + ν 2 | ξ | 4 ) 1 / 4 ( 1 + ( h 2 ( x , ξ ) ) 2 ) 1 / 4 ( 1 + ν 2 | ξ | 4 ) 1 / 4   where ν   is the ellipticity constant of ( a j k )   .
The norm N   on ( H s ) 2   is defined by
( N ( v ) ) 2 = K Λ v L 2 2 + v H s 1 2 . (2.28)
It will be shown that N   is equivalent to the standard H s   -norm.
Let k ~ ( x , ξ ) = ( e p ( x , ξ ) ( 1 + ( h 2 ( x , ξ ) ) 2 ) s / 4 0 0 e p ( x , ξ ) ( 1 + ( h 2 ( x , ξ ) ) 2 ) s / 4 )   and K ~ = Ψ k ~   . Then k ~ ( S s ) 2 × 2   and by the symbolic calculus in Theorem  2.1.2  K ~ K = I + Ψ r 4 ,   for some r 4 ( S 1 ) 2 × 2   , where | r 4 | ( S 1 ) 2 × 2 ( j 0 )   depends on ν , ( a j k )   and c 0   for each j 0 N   . Therefore K ~ K Λ = ( I + Ψ r 4 ) ( I S ) = I ( S Ψ r 4 + Ψ r 4 S ) ,   where S Ψ r 4 + Ψ r 4 S   has order 1   . By the Sobolev boundedness (Theorem  2.1.1 )
v H s 2 2 K ~ K Λ v H s 2 + 2 ( S Ψ r 4 + Ψ r 4 S ) v H s 2
c ( K Λ v L 2 2 + v H s 1 2 ) c v H s 2 ,
for a sufficient large constant c = c ( n , s , ν , ( a j k ) , c 0 , ( b l m j ) )   independent of v   . This shows the equivalence of the norms.
Next, we shall estimate the norm N   to establish the inequalities (i), (iii) and (iv) in Lemma  2.3.1 In the following c j   will denote a constant depending on n , s , ν , ( a j k )   , c 0   and finitely many derivatives of ( b l m j )   and ( c l m )   .
To estimate the second term of ( N ( v ) ) 2   in ( 2.28 ) we write
t w H s 1 2 = t J s 1 w , J s 1 w L 2 = 2 Re J s 1 t w , J s 1 w L 2
= 2 Re J s 1 ( i H w + B w + C w + f ) , J s 1 w L 2
= 2 Re i H J s 1 w , J s 1 w L 2
+ 2 Re ( i [ J s 1 H H J s 1 ] + J s 1 B + J s 1 C ) w , J s 1 w L 2
+ 2 Re J s 1 f , J s 1 w L 2
c 1 ( N ( w ) ) 2 + c 2 min { N ( f ) N ( w ) ; ( λ ( | x | ) ) 1 J s 1 / 2 f , J s 1 / 2 f L 2 } .
Above we have used that H   is self-adjoint and diagonal, that N   and H s   -norm are equivalent, and that λ   is bounded above.
For the first term of ( N ( w ) ) 2   in ( 2.28 ) we write
t K Λ w L 2 2 = t K Λ w , K Λ w L 2 = 2 Re t K Λ w , K Λ w L 2 = 2 Re K t Λ w , K Λ w L 2 = 2 Re K ( i H + B d i a g ) Λ w + Ψ r 3 w + Λ f , K Λ w L 2 = 2 Re K ( i H + B d i a g ) Λ w , K Λ w L 2 + 2 Re K Ψ r 3 w , K Λ w L 2 + 2 Re K Λ f , K Λ w L 2 2 Re K ( i H + B d i a g ) Λ w , K Λ w L 2 + c ( N ( w ) ) 2 + 2 Re K Λ f , K Λ w L 2 = I + c ( N ( w ) ) 2 + I I I , (2.29)
since K Ψ r 3 w L 2 c w H s c ( N ( w ) )   by Sobolev boundedness and norm equivalence.
We should consider the terms I   and I I I   separately. First we have
I = 2 Re K ( i H + B d i a g ) Λ w , K Λ w L 2
= 2 Re ( i H + B d i a g ) K Λ w , K Λ w L 2 + 2 Re i [ K H H K ] Λ w , K Λ w L 2
+ 2 Re [ K B d i a g B d i a g K ] Λ w , K Λ w L 2
2 Re ( B d i a g K + i [ K H H K ] ) Λ w , K Λ w L 2 + c ( N ( w ) ) 2 ,
since H   is self-adjoint and [ K B d i a g B d i a g K ]   is a commutator of diagonal matrices and therefore has order s   . Using the commutator formula of the symbolic calculus on the diagonal matrices K   and H   it follows that i [ K H H K ] = Ψ q + Ψ r 5 ,   where r 5 S s   and q S s + 1   is given by q = ( { e p ( 1 + h 2 2 ) s / 4 ; h 2 + h 1 } 0 0 { e p ( 1 + h 2 2 ) s / 4 ; h 2 h 1 } ) .   Since { e p ( 1 + h 2 2 ) s / 4 ; h 2 + h 1 } = { h 2 ; p } e p ( 1 + h 2 2 ) s / 4 + { e p ( 1 + h 2 2 ) s / 4 ; h 1 } ,   where the last term is in S s   , it follows that q = ( { h 2 ; p } 0 0 { h 2 ; p } ) k + r 6 ,   for some r 6 S s   , and thus i [ K H H K ] = Ψ { h 2 ; p } K + Ψ r 7 ,   with r 7 S s   . Hence,
I = 2 Re K ( i H + B d i a g ) Λ w , K Λ w L 2
2 Re ( B d i a g Ψ { h 2 ; p } ) K Λ w , K Λ w L 2 + c ( N ( w ) ) 2 .
Next we apply the sharp Gårding inequality (Theorem  2.1.3 ) to the diagonal matrix B d i a g Ψ { h 2 ; p }   .
For l = 1 , 2   ,
Re ( i j = 1 n b l l j ( x ) ξ j { h 2 ; p } ) = j = 1 n Im b l l j ( x ) ξ j { h 2 ; p }
c 0 λ ( | x | ) j = 1 n | ξ j | C c 0 λ ( | x | ) | ξ | + c ( n C ) c 0 λ ( | x | ) | ξ | + C .
Choosing C = 1 + n   and using that ( 1 + | ξ | 2 ) 1 / 2 1 + | ξ |   , we obtain Re ( i j = 1 n b l l j ( x ) ξ j { h 2 ; p } ) c 0 λ ( | x | ) ( 1 + | ξ | 2 ) 1 / 2 + c 0 + C ,   and the sharp Gårding inequality yields
2 Re ( B d i a g Ψ { h 2 ; p } ) K Λ w , K Λ w L 2 + c ( N ( w ) ) 2
2 Re c 0 λ ( | x | ) J 1 K Λ w , K Λ w L 2 + c ( N ( w ) ) 2 .
Since ( λ ( | x | ) ) 1 / 2 C b   , one has λ J 1 = Ψ λ ξ = Ψ ( λ ) 1 / 2 ( ξ ) 1 / 2 Ψ ( λ ) 1 / 2 ( ξ ) 1 / 2 + Ψ r 8 ,   with r 8 S 0   and ξ = ( 1 + | ξ | 2 ) 1 / 2   . But Ψ ( λ ) 1 / 2 ( ξ ) 1 / 2 = ( Ψ ( λ ) 1 / 2 ( ξ ) 1 / 2 ) * + Ψ r 9 ,   for some r 9 S 1 / 2   , so 2 Re λ J 1 K Λ w , K Λ w L 2 2 ( λ ) 1 / 2 J 1 / 2 K Λ w L 2 2 + c ( N ( w ) ) 2 .   Recalling that I = K ~ K Λ + Ψ r 1   , r 1 S 1   and using the symbolic calculus on diagonal matrices we obtain that
λ 1 / 2 J s + 1 / 2 w L 2 ( λ 1 / 2 J 1 / 2 ) ( J s K ~ ) ( K Λ ) w L 2 + c N ( w )
( J s K ~ ) ( λ 1 / 2 J 1 / 2 ) ( K Λ ) w L 2 + c N ( w ) c ( ( λ 1 / 2 J 1 / 2 ) ( K Λ ) w L 2 + N ( w ) ) .
So the following estimate for the term I   in ( 2.29 ) is therefore obtained
I = 2 Re K ( i H + B d i a g ) Λ w , K Λ w L 2
c λ ( | x | ) J s + 1 / 2 w , J s + 1 / 2 w L 2 + c ( N ( w ) ) 2 .
The estimate for the term I I I   in ( 2.29 ) will depend on which inequality (i), (ii), (iv) in Lemma  2.3.1 is considered.
To obtain (i) we write, I I I = 2 Re K Λ f , K Λ w L 2 c N ( f ) N ( w ) .   Adding the estimates for the two terms of t ( N ( w ) ) 2   we get t ( N ( w ) ) 2 c ( N ( w ) ) 2 + c N ( f ) N ( w ) ,   so t N ( w ) c N ( w ) + c N ( f ) .   Hence, N ( w ( t ) ) e c t ( N ( w ( 0 ) ) + c 0 t N ( f ( τ ) ) d τ ) ,   which proves part (i) in Lemma  2.3.1 .
To obtain (iii), we use again that I I I c N ( f ) N ( w ) .   Therefore, adding estimates one finds that t ( N ( w ) ) 2 + c λ J s + 1 / 2 w , J s + 1 / 2 w L 2 c ( N ( w ) ) 2 + c N ( f ) N ( w ) .   Integration from 0   to t   yields (iii).
Finally to obtain (iv) we estimate the term I I I   as follows
I I I = 2 Re K Λ f , K Λ w L 2 = 2 Re J 1 / 2 J 1 / 2 K Λ f , K Λ w L 2
= 2 Re J 1 / 2 K Λ f , J 1 / 2 K Λ w L 2 = 2 Re λ 1 / 2 J 1 / 2 K Λ f , λ 1 / 2 J 1 / 2 K Λ w L 2
2 λ 1 / 2 J 1 / 2 K Λ f L 2 λ 1 / 2 J 1 / 2 K Λ w L 2
= 2 λ 1 / 2 ( J 1 / 2 K Λ J s + 1 / 2 ) ( J s 1 / 2 f ) L 2
λ 1 / 2 ( J 1 / 2 K Λ J s 1 / 2 ) ( J s + 1 / 2 w ) L 2
1 ε R n | ( J 1 / 2 K Λ J s + 1 / 2 ) ( J s 1 / 2 f ) | 2 λ 1 d x
+ ε R n | ( J 1 / 2 K Λ J s 1 / 2 ) ( J s + 1 / 2 w ) | 2 λ d x
c ε R n | J s 1 / 2 f | 2 λ 1 d x + c ε R n | J s + 1 / 2 w ) | 2 λ d x ,
where the last estimate follows from Theorem  2.1.1 since J 1 / 2 K Λ J s + 1 / 2   and J 1 / 2 K Λ J s 1 / 2   have order zero. ε   is a small constant to be chosen. Adding estimates we obtain
t ( N ( w ) ) 2 + ( c 1 c 2 ε ) λ J s + 1 / 2 w , J s + 1 / 2 w L 2
c ( N ( w ) ) 2 + ( c 3 ε + c 4 ) λ 1 J s 1 / 2 f , J s 1 / 2 f L 2 .
Then, choosing ε   sufficiently small, integrating from 0   to t   we obtain part (iv) of Lemma  2.3.1 .
Proof of Theorem  2.3.1 Uniqueness :
Assume w ( C ( [ 0 , T ] : H s ( R n ) ) ) 2   is a solution of the system in ( 2.20 ) with f = 0   and w 0 = 0   . Then w ( C ( [ 0 , T ] : H s ( R n ) ) C 1 ( [ 0 , T ] : H s 2 ( R n ) ) ) 2   . From part (i) in Lemma  2.3.1 one concludes that w = 0   .
Existence :Case 1 : f ( S ( R n + 1 ) ) 2   and w 0 ( S ( R n ) ) 2   :
The conjugate linear functional l *   is defined in the linear subspace ( t + ( i H + B + C ) * ) ( C 0 ( R n × [ 0 , T ) ) ) 2 ( L 1 [ 0 , T ] : H s ( R n ) ) ) 2 ,   by l * ( η ) = 0 T f , φ L 2 × L 2 d t + w 0 , φ ( 0 ) L 2 × L 2   for φ ( C 0 ( R n × [ 0 , T ) ) ) 2   and η = ( t + ( i H + B + C ) * ) φ   . This is well defined by the uniqueness part above.
By part (ii) of Lemma  2.3.1 with s   replaced by s   it follows that
| l * ( η ) | 0 T f ( H s ) 2 φ ( H s ) 2 d t + w 0 ( H s ) 2 φ ( 0 ) ( H s ) 2
c 1 e c 2 T ( f ( L 1 ( [ 0 , T ] : H s ) ) 2 + w 0 ( H s ) 2 ) η ( L 1 ( [ 0 , T ] : H s ) ) 2 .
Using the Hahn-Banach theorem to extend l *   , there exists w ( L ( [ 0 , T ] : H s ( R n ) ) ) 2   such that
0 T w , ( t + ( i H + B + C ) * ) φ ( H s ) 2 × ( H s ) 2 = 0 T f , φ L 2 × L 2 d t + w 0 , φ ( 0 ) L 2 × L 2 , φ ( C 0 ( R n × [ 0 , T ) ) ) 2 . (2.30)
Thus, ( t ( i H + B + C ) ) w = f   as distributions for 0 < t < T   . From this equation one has that t w ( L ( [ 0 , T ] : H s 2 ( R n ) ) ) 2   since f ( S ( R n + 1 ) ) 2   , so w ( C ( [ 0 , T ] : H s 2 ( R n ) ) ) 2   . Using the equation once more, w ( C 1 ( [ 0 , T ] : H s 4 ( R n ) ) ) 2   , and w ( 0 ) = w 0   by ( 2.30 ).
Since w 0 ( S ( R n ) ) 2   , s   can be replaced by s + 4   in the previous argument and there is a solution w   of ( 2.20 ) to which Lemma  2.3.1 parts (i)-(iv) hold.
Case 2 : w 0 ( H s ( R n ) ) 2   :
Choose a sequence ( v j )   in S ( R n ) ) 2   such that v j w 0   in ( H s ( R n ) ) 2   .
(A) If f ( L 1 ( [ 0 , T ] : H s ) ) 2   , choose a sequence ( f j )   in ( S ( R n + 1 ) ) 2   such that f j f   in ( L 1 ( [ 0 , T ] : H s ( R n ) ) ) 2   .
By case 1 there is a solution w j ( C ( [ 0 , T ] : H s + 2 ( R n ) ) ) 2   of ( 2.20 ) with f   and w 0   replaced by f j   and v j   respectively. Using Lemma  2.3.1 , part (i), it follows that ( w j )   is a Cauchy sequence in ( C ( [ 0 , T ] : H s ( R n ) ) 2   and that the limit w   is a solution of ( 2.20 ) satisfying part (i) in Lemma  2.3.1 .
(B) If f ( L 2 ( [ 0 , T ] : H s ) ) 2   , choose a sequence ( f j )   in ( S ( R n + 1 ) ) 2   such that f j f   in ( L 2 ( [ 0 , T ] : H s ( R n ) ) ) 2   . Procceding as in (A), there is a solution w ( C ( [ 0 , T ] : H s ( R n ) ) ) 2   of ( 2.20 ) satisfying (iii) of Lemma  2.3.1 .
(C) Let J s 1 / 2 f ( L 2 ( R n × [ 0 , T ] : d x d t λ ( | x | ) ) ) 2 .   By Theorem  2.1.1 there exists a sequence ( g j )   in ( S ( R n + 1 ) ) 2   , such that g j J s 1 / 2 f   in ( L 2 ( R n × [ 0 , T ] : d x d t λ ( | x | ) ) ) 2   .
Procceding as in (A) with f j   replaced by J s 1 / 2 g j ( S ( R n + 1 ) ) 2   , there is a solution w ( C ( [ 0 , T ] : H s ( R n ) ) ) 2   of ( 2.20 ) satisfying (iv) of Lemma  2.3.1 .
This completes the proof of Theorem  2.3.1 .
Remark 2.3.1 Suppose that the differential operators B 11   and B 22   in the entries of B   in ( 2.20 ) are replaced by pseudo-differential operators Ψ b 11   and Ψ b 22   of order 1   and suppose that C   in ( 2.20 ) is replaced by a 2 × 2   matrix of pseudo-differential operators of order 0   . Then the conclusion of Theorem  2.3.1 still holds if | R e b l l ( x , ξ ) | c 0 λ ( | x | ) ξ , l = 1 , 2 , x , ξ R n .   The reason is that the application of the sharp Gårding inequality in the proof of Lemma  2.3.1 goes through in exactly the same way.
Consequently, if b 1 ( x ) x   in ( 2.20 ) is replaced by Ψ b   with b S 1   and c 1 , c 2   in ( 2.20 ) are replaced by pseudo-differential operators of order 0   , then the conclusion of Corollary  2.3.1 holds if | R e b ( x , ξ ) | c 0 λ ( | x | ) ξ , x , ξ R n .  
This will be useful later.

3 A NEW CLASS OF SYMBOLS

As it has been shown in [18, 22to obtain local well posedness for nonlinear Schrödinger equations one relies on certain smoothing effects for the associated linear equation with lower order terms (order zero and one). In the previous section we have established these smoothing effects in equations with variable second order elliptic coefficients by using known properties of classical pseudo-differential operators. In an attempt to prove these smoothing effects for the non-elliptic case, one is led (see [4) to the study of certain operators with non-standard symbols. Our goal in this section is to study results concerning the L 2   -boundedness and composition of operators in this class by using geometric arguments. The elliptic case of our results were proved by Craig, Kappeler and Strauss in [4whose statements we follow. The differences between the elliptic and non-elliptic settings are highlighted in Proposition  3.1.2 below, -see also [21.

3.1 Symbol Properties

To begin with, the symbols of interest will be compared to the classical ones defined in ( 2.2 ) Section  2 .
We recall the following spaces
S ( R n ) = { u C ( R n ) : sup x R n x k | x α u ( x ) | < , k N , α N n } , (3.1)
with seminorms
| u | S , m = max k + | α | m x k x α u ( x ) L ( R n ) , (3.2)
and
C b ( R n ) = { u C ( R n ) : sup x R n | x α u ( x ) | < , α N n } , (3.3)
with seminorms
| u | C b , m = max | α | m x α u ( x ) L ( R n ) . (3.4)
The symbol a = a ( x , ξ ) C ( R n × R n )   will satisfy certain estimates and the operator Ψ a   associated with the symbol a   will be defined as Ψ a u ( x ) = R n e i x ξ a ( x , ξ ) u ^ ( ξ ) d ξ ( 2 π ) n / 2 , u S ( R n ) .  
Proposition 3.1.1 (i) Suppose a   is a classical symbol of order m R   , a S 1 , 0 m   , i.e. for α , β N n  
| x α ξ β a ( x , ξ ) | c α , β ξ m | β | , x , ξ R n . (3.5)
Then Ψ a   is a continuous map from S ( R n )   into S ( R n )   .
(ii) Suppose m R   and a   satisfies that for α N n  
| x α a ( x , ξ ) | c α ξ m , x , ξ R n . (3.6)
Then Ψ a   is a continuous map from S ( R n )   into C b ( R n )   .
Proof of Proposition  3.1.1  See [25.
The symbols of interest in this section satisfy estimates of the type
| x α ξ β a ( x , ξ ) | c α β x | β | ξ m | β | . (3.7)
This is better than ( 3.6 ) in part (ii) of Proposition  3.1.1 , but not as good as (i). In particular, we will see (Proposition  3.1.2 , part (iv)) that there exist a   satisfying ( 3.7 ) with m = 0   and v S ( R n )   such that Ψ a v / S ( R n )   .
Let A ( x )   be a real, symmetric and invertible n × n   matrix. Using a coordinate change (a rotation and dilations) in the x   -variable, there are essentially only the elliptic and ultrahyperbolic cases A e = I n and A h = ( I k 0 0 I n k ) , k { 1 , . . , n 1 } .   where I j   is the j × j   unit matrix.
Take χ C ( R n )   with χ ( t ) = 0   for | t | 1   and χ ( t ) = 1   for | t | 2   .
Definition 3.1.1 (i) It will be said that a S ( R n : S 1 , 0 m )   if a C ( R n × R n × R n )   and a   satisfies
| s μ s α x β ξ γ a ( s ; x , ξ ) | c μ α β γ ξ m | γ | , s , x , ξ R n , μ , α , β , γ N n . (3.8)
(ii) For a S ( R n : S 1 , 0 m )   , let
{ ( 1 ) b e ( x , ξ ) = χ ( | ξ | ) a ( P ( x , A e ξ ) ; x , ξ ) , ( 2 ) b h ( x , ξ ) = χ ( | ξ | ) a ( P ( x , A h ξ ) ; x , ξ ) , (3.9)
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   , (notice that P ( y , z )   is homogeneous of degree 0   in z   ).
Remark 3.1.1 (a) Although we shall work in the class S ( R n : S 1 , 0 m )   , it will be clear that all the results deduced for this class still hold when just a finite number of semi-norms in S   and S 1 , 0 m   are assumed to be finite, i.e. ( 3.8 ) with | μ | + | α | + | β | + | γ | N   for N   large enough.
(b) We observe that if a S ( R n : S 1 , 0 m )   , then ξ α ( ξ β a ( ) ) S ( R n : S 1 , 0 k ) , k = m + | β | | α | ,   and for M N   large enough x M χ ( | ξ | ) a ( P ( x , A l ξ ) ; x , ξ ) = x M b l ( x , ξ ) , with l = e , or l = h ,   is “roughly speaking” a symbol in the class S 1 , 0 m   (when only finitely many ξ   derivatives are taken into account, which is always the case in the sequel).
(c) Finally notice that if a S ( R n : S 1 , 0 m )   and b l   is defined as in ( 3.9 ) then b β   given by b l β = x | β | ξ β b l ( x , ξ ) , with l = e , or l = h ,   is a symbol of the same type and the corresponding bounds in ( 3.8 ) are controlled by those of a   .
(d) The symbols described here will be our basic building blocks in the study of variable coefficients Schrödinger operators, see Section  5 .
The symbols defined in ( 3.9 ) satisfy an estimate of the type given in Proposition  3.1.1 . More precisely, if u S ( R n )   then Ψ b e u S ( R n )   and Ψ b h u C b ( R n )   and is rapidly decreasing away from the characteristic directions, i.e. A x x = 0   . In the characteristics directions, Ψ b h u ( x )   decays as | x | 1 n   as | x |   .
Proposition 3.1.2 Let u S ( R n )   , and a   , b e   and b h   as in definition  3.1.1 (i) If | α + β | k   , then for all x   in R n   | x α x β Ψ b e u ( x ) | c k ( sup | α | k ; | β | k ; | γ | k s k | s α x β ξ γ a ( s ; x , ξ ) | ) | u ^ | S , 2 k + m + n + 1 .   (ii) If c [ 0 , 1 ]   and | α + β | k   , then if | A h x x | | A h x | | x | c   c k | x α x β Ψ b h u ( x ) | c k ( sup | α | k ; | β | k ; | γ | k s k | s α x β ξ γ a ( s ; x , ξ ) | ) | u ^ | S , 2 k + m + n + 1 .   (iii) If a C 0 ( R n : S 1 , 0 m )   with a ( s ; x , ξ ) = 0   if | s | > 1   and | β | k   , then | x β Ψ b h u ( x ) | c k ( sup | α | k ; | β | k ; | γ | k | s α x β ξ γ a ( s ; x , ξ ) | ) | u ^ | S , 2 k + m + n + 1 x 1 n , x R n .   (iv) There exist a S ( R n : S 1 , 0 m )   , in fact, a C 0 ( R n : S 1 , 0 m )   , v S ( R n )   , and c > 0   such that | Ψ b h v ( x ) | c | x | 1 n , with A h x x = 0 , | x | 10  
In the proof of Proposition  3.1.2 we will use the following results.
Proposition 3.1.3 (i) Let α N n   . Then ( ξ ξ ) ξ α = ξ α ( ξ ξ ) | α | ξ α .   (ii) Let T t   denote the transpose of the operator T   , i.e. T u v = u T t v   . Then [ 1 i x ξ ( ξ ξ ) ] t = 1 i x ξ ( ( ξ ξ ) ( n 1 ) I ) ,   and ( ξ ξ ) [ 1 i x ξ ] = 1 i x ξ .   (iii) Let φ   be a differentiable and homogeneous function of degree 0   on R n { 0 }   . Then ( ξ ξ ) φ ( ξ ) = 0 , ξ 0 .   (iv) Let α , β N n   , k N   , b = b e   or b = b h   . Then | ( ξ ξ ) k x α ξ β b ( x , ξ ) | c k ( sup | γ | k | s α x β ξ γ a ( s ; x , ξ ) | ) x | β | ξ m | β | .  
Proof of Proposition  3.1.3 (i) and (ii) are easily verified. As for (iii) we have
( ξ ξ ) [ φ ( ξ | ξ | ) ] = j = 1 n ξ j ξ j [ φ ( ξ | ξ | ) ]
= j = 1 n ξ j k = 1 n ( ξ k φ ) ( ξ | ξ | ) ξ j ( ξ k | ξ | ) = j , k ( ξ k φ ) ( ξ | ξ | ) ξ j ( δ j k | ξ | ξ j ξ k | ξ | 3 )
= k ( ξ k φ ) ( ξ | ξ | ) ξ k | ξ | k ( ξ k φ ) ( ξ | ξ | ) ξ k j ξ j 2 | ξ | 3 = 0 .
(iv) follows from (i), (iii) and Definition  3.1.1 , since P ( x , A ξ )   is homogeneous of degree 0   in ξ   .
Proof of Proposition  3.1.2 For simplicity of the exposition we shall assume that the constants c μ α β γ   which appear in ( 3.8 ) are all smaller than unity for | μ | k   , | α | k   , | β | k   and | γ | k   . Also we shall drop the powers of 2 π   which appear in the definition of Ψ b   .
(i) Let | α + β | k   . We consider three cases:Case 1 : | x | 2   . Here | x α x β Ψ b e u ( x ) | c k | u ^ | S , m + k + n + 1 .   Case 2 : | x | 2   and | P ( x , ξ ) | | x | / 2   . From ( 3.7 ), it suffices to estimate terms of the type
I = e i x ξ ξ α 1 [ ξ β 1 ] x β 2 ξ α 2 [ b e ( x , ξ ) ] ξ α 3 u ^ ( ξ ) d ξ , (3.10)
where α 1 + α 2 + α 3 = α   , and β 1 + β 2 = β   .
Pythagoras' theorem gives | x | 2 = | x ξ | 2 | ξ | 2 + | P ( x , ξ ) | 2 ,   so by hypothesis
| x ξ | = | x | | ξ | 1 | P ( x , ξ ) | 2 | x | 2 | x | | ξ | ( 1 | P ( x , ξ ) | 2 | x | 2 )
| x | | ξ | ( 1 1 4 ) 3 4 | x | | ξ | .
Above we have used that 1 t 1 t   for t [ 0 , 1 ]   . Now we use the identity e i x ξ = ( 1 i x ξ ( ξ ξ ) ) k e i x ξ .   By Proposition  3.1.3 (ii), if k 1 + k 2 + k 3 k   it suffices to estimate terms of the type
I I = e i x ξ ( i x ξ ) k ( ξ ξ ) k 1 [ ξ α 1 [ ξ β 1 ] ] ( ξ ξ ) k 2 x β 2 ξ α 2 b e ( x , ξ ) ( ξ ξ ) k 3 ξ α 3 u ^ ( ξ ) d ξ .
Using Proposition  3.1.3 (iv) it follows that
| I I | c k | ξ | 1 1 | x | k | ξ | k ξ | β 1 α 1 | + m | α 2 | x | α 2 | | u ^ | S , 2 k + m + n + 1 c k | u ^ | S , 2 k + m + n + 1 .
Case 3 : | x | 2   and | x | < 2 | P ( x , ξ ) |   . From ( 3.8 ) and ( 3.10 ) it follows that | I | c k | u ^ | S , m + k + n + 1 ,   which proves (i).
(ii) Let c [ 0 , 1 ]   , | α + β | k   . Suppose that | x | 2 | A h x x | c   with c > 0   , (otherwise the statement is trivial). We consider three cases:
Case 1 : | x | 4 / c   . Hence c k | x α x β Ψ b h u ( x ) | 4 k c k | u ^ | S , k + m + n + 1 .   Case 2 : | x | 4 / c   and | P ( x , A h ξ ) | c | x | / 2   .
Since x = x A h ξ | ξ | 2 A h ξ + P ( x , A h ξ ) ,   by Pythagoras' theorem, | x | 2 | ξ | 2 = | x A h ξ | 2 + | ξ | 2 | P ( x , A h ξ ) | 2 .   The matrix A h   is symmetric, so
| ξ A h x | = | x A h ξ | = | x | | ξ | 1 | P ( x , A h ξ ) | 2 | x | 2
| x | | ξ | ( 1 c 2 4 ) 3 2 | x | | ξ | .
Next we write ξ = ξ A h x | x | 2 A h x + P ( ξ , A h x ) .   By Pythagoras' theorem | P ( ξ , A h x ) | = | ξ | 1 | ξ A h x | 2 | ξ | 2 | x | 2 | ξ | 1 ( 1 c 2 4 ) = c 2 | ξ | .   Using that x ξ = | x | 2 ( ξ A h x ) ( A h x x ) + x P ( ξ , A h x ) ,   from the above estimates and our hypothesis we conclude that
| x ξ | c 3 2 | x | | ξ | | x | | P ( x , A h ξ ) | c 3 2 | x | | ξ | c 2 | x | ξ | c 4 | x | | ξ | (3.11)
Integrating by parts as in case 2 of part (i), (ii) follows.
Case 3 : | x | 4 / c   and | P ( x , A h ξ ) | c | x | / 2   .
From ( 3.9 ) and ( 3.10 ) as in case 3 of the elliptic case it follows that | x β Ψ b h u ( x ) | c k c 0 | u ^ | S , m + k + n + 1 .   (iii) It suffices to show the statement for | x | 10   . Let | β 1 + β 2 | k   and consider terms of the type I = e i x ξ ( i ξ ) β 1 x β 2 b h ( x , ξ ) u ^ ( ξ ) d ξ .   Since | P ( x , A h ξ ) | 1   proceeding as before x = x A h ξ | ξ | 2 A h ξ + P ( x , A h ξ ) ,   and by Pythagoras' theorem | x | 2 | ξ | 2 = | x A h ξ | 2 + | ξ | 2 | P ( x , A h ξ ) | 2 .   The matrix A h   is symmetric, so | ξ A h x | = | x A h ξ | = | x | | ξ | 1 | P ( x , A h ξ ) | 2 | x | 2 | x | | ξ | ( 1 | x | 2 ) .   Next we write ξ = ξ A h x | x | 2 A h x + P ( ξ , A h x ) .   By Pythagoras' theorem | P ( ξ , A h x ) | = | ξ | 1 | ξ A h x | 2 | ξ | 2 | x | 2 | ξ | 1 ( 1 | x | 2 ) 2 | ξ | | x | 1 ,   i.e.
| P ( ξ , A h x ) | | ξ | 2 | x | 1 . (3.12)
Hence ξ   is in a cone Γ x   with vertex at the origin, axis given by A h x   and opening angle θ   where s i n ( θ ) = 2 | x | 1   . In particular, θ 2 | x | 1   , because | x | 10   .
Therefore I = Γ x e i x ξ ( i ξ ) β 1 x β 2 b h ( x , ξ ) u ^ ( ξ ) d ξ   and
| I | c | u ^ | S , k + m + n + 1 Γ x ξ n 1 d ξ
c | u ^ | S , k + m + n + 1 ( | x | 1 ) n 1 R n ξ n 1 d ξ
= c | u ^ | S , k + m + n + 1 | x | n + 1 .
(iv) Choose φ C 0 ( R n )   with φ 0   , φ ( s ) = 1   if | s | 1 / 4   and φ ( s ) = 0   if | s | 1 / 2   . Let a ( s ; x , ξ ) = φ ( s ) ξ m   . Choose v S ( R n )   such that v ^ ( ξ ) = χ ( | ξ | ) ( 1 χ ( | ξ | ) ) 2   with χ ( )   given in Definition  3.1.1 . Then
Ψ b h v ( x ) = e i x ξ φ ( P ( x , A h ξ ) ) ξ m ( χ ( | ξ | ) ) 2 ( 1 χ ( | ξ | ) ) 2 d ξ . (3.13)
Assume that | x | 10   . An argument similar to that used to deduce ( 3.12 ) shows that | P ( ξ , A h x ) | ( 2 ) 1 | ξ | | x | 1 .   Also one has that x ξ = x P ( ξ , A h x )   when A h x x = 0   . Hence | x ξ | ( 2 ) 1 | ξ | 2   on the ξ   support of the integrand in ( 3.13 ). Thus R e Ψ b h v ( x ) cos 2 φ ( P ( x , A h ξ ) ) ξ m χ ( | ξ | ) 2 ( 1 χ ( | ξ | ) ) 2 d ξ .   Let Γ x   be the cone defined by Γ x = { ξ R n : | ξ | 1 | P ( ξ , A h x ) | ( 4 2 | x | ) 1 } .   Suppose ξ Γ x   . Since ξ = ξ A h x | x | 2 A h x + P ( ξ , A h x ) ,   it follows from Pythagoras' theorem that | x A h ξ | = | ξ A h x | = | ξ | | x | 1 | P ( ξ , A h x ) | 2 | x | 2 | ξ | | x | ( 1 1 32 | x | 2 ) .   Now write x = x A h ξ | ξ | 2 A h ξ + P ( x , A h ξ ) ,   and use Pythagoras' theorem to get | P ( x , A h ξ ) | 2 = | x | 2 | x A h ξ | 2 | ξ | 2 | x | 2 | x | 2 ( 1 1 32 | x | 2 ) 2 1 16 .   Therefore
R e Ψ b h v ( x cos ( 2 ) Γ x ξ m ( χ ( | ξ | ) ) 2 ( 1 χ ( | ξ | ) ) 2 d ξ
= cos ( 2 ) c n | x | 1 n R n ξ m ( χ ( | ξ | ) ) 2 ( 1 χ ( | ξ | ) ) 2 d ξ = c c n cos ( 2 ) | x | 1 n ,
where c = R n ξ m ( χ ( | ξ | ) ) 2 ( 1 χ ( | ξ | ) ) 2 d ξ > 0 .   This completes the proof of Proposition  3.1.2 .

3.2 L 2   -boundedness

Now we are ready to establish the L 2   -boundedness of our ultrahyperbolic operators.
Theorem 3.2.1 Suppose a S ( R n : S 1 , 0 m )   (see Definition  3.1.1 , ( 3.8 )-( 3.9 )). Then there exists c = c ( m )   and N = N ( n )   such that Ψ b h u L 2 c max μ + | α + β + γ | N s μ ξ m + | γ | s α x β ξ γ a L u H m , u S ( R n ) .  
The proof of this theorem in the 2-dimensional case ( n = 2   ) is more involved than in the higher dimensional case. As it was pointed out in Proposition  3.1.2 (iv) one can not expect enough decay in cones around the characteristic directions (i.e. A h x x = 0   ). In fact, for n = 2   the estimate | x | 1   is critical and after some decomposition in frequency and space we need to use Cotlar-Stein lemma to glue the pieces together.
Proof of Theorem  3.2.1 If a S ( R n : S 1 , 0 m )   , then ξ m a S ( R n : S 1 , 0 0 )   with
max μ + | α + β + γ | N s μ ξ | γ | s α x β ξ γ [ ξ m a ] L
c m max μ + | α + β + γ | N s μ ξ m s α x β ξ γ a L .
Furthermore, J m   is an isometry of H m   onto L 2   , so it can be assumed that m = 0   .
Next, we use an argument similar to that in [21Section 3 (which in the case of the elliptic operators gives a straightforward proof of the L 2   -continuity). This is based on a simple change of variables -see ( 3.14 ) below. In order to do it we make the following decomposition. Let { φ j } j N   be a smooth partition of unity on R   subordinate to the covering R = ( 1 , 1 ) j = 1 { t R : 2 j 2 < | t | < 2 j }   i.e. satisfying 0 φ j ( x ) 1 , 1 = j = 0 φ j   , supp φ 0 ( 1 , 1 )   and supp φ j { t R : 2 j 2 < | t | < 2 j }   .
Since the intervals in the covering have length at least 1   , it can be assumed that | d k d t k φ j ( t ) | c k , t R , k N .   Let φ C ( R )   with φ ( t ) = 1   for t 1   , φ ( t ) = 0   for t 2   and 0 φ 1   .
Denote by a ^   the Fourier transform of a   in the s   -variable. Thus, σ N 1 σ β 1 a ^ ( σ ; x , ξ ) L c max N 2 + | β 2 | 2 N 1 + n + 1 s N 2 s β 2 a L .   Now using that σ P ( x , A h ξ ) = x P ( σ , A h ξ )   we have
Ψ b h u ( x ) = e i x ξ a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) u ^ ( ξ ) d ξ
= ( e i x ( ξ + P ( σ , A h ξ ) ) a ^ ( σ ; x , ξ ) χ ( | ξ | ) u ^ ( ξ ) d ξ ) d σ
= j = 0 ( e i x ( ξ + P ( σ , A h ξ ) ) φ j ( | σ | ) a ^ ( σ ; x , ξ ) χ ( | ξ | ) ( 1 φ ( δ | ξ | 2 j + 1 ) ) u ^ ( ξ ) d ξ ) d σ
+ j = 0 ( e i x ( ξ + P ( σ , A h ξ ) ) φ j ( | σ | ) a ^ ( σ ; x , ξ ) χ ( | ξ | ) φ ( δ | ξ | 2 j + 1 ) u ^ ( ξ ) d ξ ) d σ
= I + I I ,
where δ > 0   is a small constant to be chosen. In I   we make the change of variables
η = ξ + P ( σ , A h ξ ) = ξ + ( σ σ A h ξ | ξ | 2 A h ξ ) . (3.14)
Then ( η ξ ) = I k = 1 n σ k | ξ | M k ( ξ ) ,   where M k   is a matrix whose entries are homogeneous of degree 0   in ξ   . The determinant function is continuous on R n 2   , so 1 2 | det ( η ξ ) | 2 ,   by fixing δ > 0   sufficiently small, we recall that | σ | 2 j , δ | ξ | 2 j + 1   . This gives
I = j = 0 φ j ( | σ | ) e i x η a ^ ( σ ; x , ξ ( η ) ) χ ( | ξ ( η ) | )
( 1 φ ( δ 2 j 1 | ξ ( η ) | ) ) u ^ ( ξ ( η ) ) | det ( η ξ ) | 1 d η d σ .
Combining Minkowski's integral inequality and the L 2   -boundedness of S 1 , 0 0   pseudo-differential operators we get I L x 2 c max μ + | α + β + γ | N s μ ξ | γ | s α x β ξ γ a L u L 2 j = 0 φ j ( | σ | ) σ n + 1 d σ ,   where j = 0 φ j ( | σ | ) σ n + 1 d σ < .   As for I I   , let a j k ( s ; x , ξ ) = φ k ( | s | ) e i σ s φ j ( | σ | ) a ^ ( σ ; x , ξ ) d σ , j , k N .   Then
I I = j , k = 0 e i x ξ a j k ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) φ ( δ 2 j 1 | ξ | ) u ^ ( ξ ) d ξ = j , k = 0 Ψ j k u ( x ) , (3.15)
where Ψ j k   is a pseudo-differential operator with symbol a j k ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) φ ( δ 2 j 1 | ξ | ) .   First, by Cauchy-Schwarz, | Ψ j k u ( x ) | c 2 j n / 2 a j k L u L 2 ,   so
χ B 2 k + 10 ( 0 ) Ψ j k u L 2 c 2 ( j + k ) n / 2 a j k L u L 2 . (3.16)
Next, we shall estimate ( 1 χ B 2 k + 10 ( 0 ) ) Ψ j k u L 2   by using Cotlar-Stein lemma which can be stated as follows.
Cotlar-Stein lemma.
Let { Ψ l } l = 0   be a sequence of bounded operators on L 2   and let { γ ( l ) } l = 0   be a sequence of positive numbers with l = 0 γ ( l ) <   . Suppose Ψ l 1 * Ψ l 2 , Ψ l 1 Ψ l 2 * ( γ ( l 2 l 1 ) ) 2 , l 1 , l 2 N , l 1 l 2 .   Then l = 0 Ψ l u   converges in L 2   for any u L 2   and l = 0 Ψ l u L 2 l = 0 γ ( l ) u L 2 .   For a proof of this lemma we refer to [30.
For fixed j , k 0   , let A l = { x R n : 2 k + l + 10 < | x | 2 k + l + 11 } , l N ,   and define Ψ l = Ψ l j k   by Ψ l u ( x ) = χ A l ( x ) Ψ j k u ( x ) , u S ( R n ) , l N .   Then l = 0 Ψ l u   converges to ( 1 χ B 2 k + 10 ( 0 ) ) Ψ j k u   uniformly on compact subsets of R n   . In particular in the distribution sense. By Cauchy-Schwarz Ψ l u L 2 c 2 ( j + k + l ) n / 2 a j k L u L 2 , u S ( R n ) ,   so Ψ l   is bounded on L 2   for l N   .
In order to apply Cotlar-Stein lemma, it suffices to consider Ψ l 1 Ψ l 2 *   for l 1 l 2   , because ( Ψ l 2 Ψ l 1 * ) * = Ψ l 1 Ψ l 2 *   and Ψ l 1 * Ψ l 2 = 0   for l 1 l 2   .
Let u S ( R n )   . Since
Ψ l * u ( x ) = e i ( x y ) ξ a ¯ j k ( P ( y , A h ξ ) ; y , ξ ) χ ( | ξ | ) φ ( δ 2 j 1 | ξ | ) χ A l ( y ) u ( y ) d y d ξ
= e i x ξ ( e i y ξ a ¯ j k ( P ( y , A h ξ ) ; y , ξ ) χ ( | ξ | ) φ ( δ 2 j 1 | ξ | ) χ A l ( y ) u ( y ) d y ) d ξ
it follows that
Ψ l 1 * v , Ψ l 2 * u = e i y ξ a ¯ j k ( P ( y , A h ξ ) ; y , ξ ) χ ( | ξ | ) φ ( δ 2 j 1 | ξ | ) χ A l 1 ( y ) v ( y )
e i z ξ a j k ( P ( z , A h ξ ) ; z , ξ ) χ ( | ξ | ) φ ( δ 2 j 1 | ξ | ) χ A l 2 ( z ) u ¯ ( z ) d ξ d z d y ,
for all v , u S ( R n )   . So Ψ l 1 Ψ l 2 * u ( y ) = K l 1 l 2 ( y , z ) u ( z ) d z ,   where
K l 1 l 2 ( y , z ) = χ A l 1 ( y ) χ A l 2 ( z ) e i ( y z ) ξ a j k ( P ( y , A h ξ ; y , ξ )
a ¯ j k ( P ( z , A h ξ ) ; z , ξ ) χ 2 ( | ξ | ) φ 2 ( δ 2 j 1 | ξ | ) d ξ .
It is convenient to fix the following terminology.
Definition 3.2.1 Γ   is the cone of angle θ ( 0 , 1 / 2 )   if there exists e 0 R n { 0 }   such that Γ = { x R n : | x | 1 | P ( x , e 0 ) | θ } .   Notice that | Γ A l | c θ n 1 2 n ( k + l )   .
Choose a collection { Γ m 1 }   of cones of angle at most 2 l 1 / 16   such that R n = m 1 Γ m 1 .   Let Γ m 1 * = { z R n : | z | 1 | P ( z , x ) | 2 l 1 16 for some x Γ m 1 } .   Then Γ m 1 *   is a cone of angle at most 2 l 1 / 4   , and it can be assumed that the collection { Γ m 1 * }   is locally finite on R n { 0 }   by choosing { Γ m 1 }   appropriately. Let S m 1 = Γ m 1 A l 1 , S m 1 * = Γ m 1 * A l 2 .  
Claim 3.2.1 :If x S m 1   , z A l 2   and a j k ( P ( x , A h ξ ) ; x , ξ ) a j k ( P ( z , A h ξ ) ; z , ξ ) 0   , then z S m 1 *   .
Proof of claim  3.2.1 : The goal is to show that z S m 1 *   . By assumption, | P ( x , A h ξ ) | , | P ( z , A h ξ ) | 2 k .   Since x = x A h ξ | ξ | 2 A h ξ + P ( x , A h ξ ) ,   by Pythagoras' theorem | x A h ξ | | x | | ξ | 1 | P ( x , A h ξ ) | 2 | x | 2 | x | | ξ | ( 1 2 2 l 1 20 ) .   Similarly, | z A h ξ | | z | | ξ | ( 1 2 2 l 2 20 ) .   Now we use the identity x z = | ξ | 2 ( x A h ξ ) ( z A h ξ ) + P ( x , A h ξ ) P ( z , A h ξ ) ,   to obtain | x z | | x | | z | ( 1 2 2 l 1 20 ) ( 1 2 2 l 2 20 ) 2 2 k .   But z = z x | x | 2 x + P ( z , x ) ,   so
| P ( z , x ) | 2 = | z | 2 ( z x ) 2 | x | 2 | z | 2 | z | 2 ( 1 2 2 l 1 19 ) + | x | 2 2 2 k + 1
2 2 l 1 19 | z | 2 + 2 2 l 1 19 2 2 l 1 16 | z | 2 ,
and therefore z Γ m 1 *   . This proves claim  3.2.1 .
Claim 3.2.2 : sup x , z R n | K l 1 l 2 ( x , z ) | c a j k L 2 2 n j 2 l 2 ( n 1 ) .  
Proof of claim  3.2.2 We recall that
K l 1 l 2 ( x , z ) = χ A l 1 ( x ) χ A l 2 ( z ) e i ( x z ) ξ a j k ( P ( x , A h ξ ) ; x , ξ )
a ¯ j k ( P ( z , A h ξ ) ; z , ξ ) χ 2 ( | ξ | ) φ 2 ( δ 2 j 1 | ξ | ) d ξ .
Suppose z A l 2   (otherwise K l 1 l 2 ( x , z ) = 0   ) and a ¯ j k ( P ( z , A h ξ ) ; z , ξ ) 0   so that | P ( z , A h ξ ) | 2 k   .
Then | ξ A h z | = | z A h ξ | | z | | ξ | ( 1 2 2 l 2 20 ) ,   as in the proof of claim  3.2.1 . Since ξ = ξ A h z | z | 2 A h z + P ( ξ , A h z ) ,   one has | P ( ξ , A h z ) | 2 = | ξ | 2 | ξ A h z | 2 | z | 2 | ξ | 2 | ξ | 2 ( 1 2 2 l 2 20 ) 2 | ξ | 2 2 2 l 2 19 ,   which shows that ξ   belongs to a cone Γ z   of angle 2 l 2 9   . Hence | K l 1 l 2 ( x , z ) | a j k L 2 | Γ z { δ | ξ | 2 j + 2 } | c a j k L 2 2 n j 2 l 2 ( n 1 ) .   This yields claim  3.2.2 .
Now to estimate Ψ l 1 Ψ l 2 * u L 2   we consider two separate cases. Case 1: l 1 l 2 l 1 + 10   .
Then
Ψ l 1 Ψ l 2 * u L 2 2 = | K l 1 l 2 ( x , z ) u ( z ) d z | 2 d x m 1 S m 1 | S m 1 * K l 1 l 2 ( x , z ) u ( z ) d z | 2 d x
m 1 | S m 1 | | S m 1 * | sup x , z | K l 1 l 2 ( x , z ) | 2 χ S m 1 * u L 2 2
c 2 l 1 ( n 1 ) 2 n ( k + l 1 ) 2 l 1 ( n 1 ) 2 n ( k + l 2 ) a j k L 4 2 2 n j 2 l 2 ( n 1 ) m 1 χ S m 1 * u L 2 2
c 2 2 n ( k + j ) 2 ( 2 n ) ( l 1 + l 2 ) a j k L 4 u L 2 2 .
The first inequality above follows from claim  3.2.1 , the second inequality from Cauchy-Schwarz, the third inequality from claim  3.2.2 , and the fourth inequality from the local finiteness of { S m 1 * }   .
Therefore,
Ψ l 1 Ψ l 2 * c 2 n ( k + j ) 2 ( 2 n ) ( l 1 + l 2 ) / 2 a j k L 2 , for l 1 l 2 l 1 + 10 . (3.17)
Before turning to the remaining case l 2 l 1 + 11   , it is useful to split A l 2   in sectors S m 2   where m 2   roughly speaking measures how noncharacteristic the directions in S m 2   are. More precisely, let
S m 2 = { z A l 2 : ( m 2 1 ) 2 l 2 | A h z z | | z | 2 m 2 2 l 2 } , m 2 = 1 , 2 , . . , 2 l 2 1 ,
S l 2 * = { z A l 2 : 1 2 | A h z z | | z | 2 1 } .
Then A l 2 = S l 2 * 2 l 2 1 m 2 = 1 S m 2   . The next result will be used to estimate the volume of S m 2   . Recall that A h = ( I k 0 0 I n k ) , k { 1 , . . , n 1 } ,   where I j   is the j × j   unit matrix.
Proposition 3.2.1 Suppose that 0 < ε a 1 / 2   , ε 1 / 16   . Let S a , ε = { z R n : | z | 1 and a ε | A h z z | | z | 2 a } .   Then | S a , ε | c n , k ε ,   where c n , k   is independent of a   .
Proof of Proposition  3.2.1 Let m j   denote the Lebesgue measure in R j   . Write S a , ε = S a , ε + S a , ε ,   where S a , ε ± = { z R n : | z | 1 and a ε ± | A h z z | | z | 2 a } .   Let z + = ( z 1 , . . , z k )   , z = ( z k + 1 , . . , z n )   . By Fubini's theorem m n ( S a , ε + ) = | z | 1 m k ( S a , ε + ( z ) ) d m n k ( z ) ,   where
S a , ε + ( z ) = { z + R k : ( z + , z ) S a , ε + }
= { z + R k : | z + | 2 + | z | 2 1 , a ε | z + | 2 | z | 2 | z + | 2 + | z | 2 a } .
By straightforward calculation, | z + | 2 | z | 2 | z + | 2 + | z | 2 [ a ε , a ] | z + | 2 [ 1 + ( a ε ) 1 ( a ε ) | z | 2 , 1 + a 1 a | z | 2 ] .   Let f ( t ) = 1 + t 1 t , t ( , 1 ) .   Then f ( t ) = 2 ( 1 t ) 2 ,   so sup [ 0 , 1 / 2 ] | f ( t ) | = 8   . By the fundamental theorem of calculus, | 1 + ( a ε ) 1 ( a ε ) 1 + a 1 a | 8 ε ,   so 1 + ( a ε ) 1 ( a ε ) 1 + a 1 a 8 ε .   Also,
1 + a 1 a 8 ε = 1 + a 1 a 1 8 ε 1 a 1 + a
1 + a 1 a ( 1 8 ε 1 a 1 + a ) 1 + a 1 a ( 1 8 ε ) .
It follows that S a , ε + ( z ) { z + R k : 1 + a 1 a ( 1 8 ε ) | z + | | z | 1 + a 1 a } ,   and hence m k ( S a , ε + ( z ) ) c k ( 1 + a 1 a ) k / 2 | z | k ε c k 3 k / 2 | z | k ε .   Therefore,
m n ( S a , ε + ) = | z | 1 m k ( S a , ε + ( z ) ) d m n k ( z )
c k ε | z | 1 | z | k d m n k ( z ) c n , k ε .
Replacing k   by n k   in the above argument for S a , ε +   , one gets m n ( S a , ε ) c n , k ε .   This proves Proposition  3.2.1 .
Claim 3.2.3 :Let m 2 { 1 , 2 , . . , 2 l 2 1 }   . Then | S m 2 | c 2 ( n 1 ) l 2 + n k .  
Proof of claim  3.2.3  This follows by Proposition  3.2.1 and homogeneity.
Claim 3.2.4 :Suppose l 2 l 1 + 11   , x A l 1   , z A l 2   , 10 m 2 2 l 2   , | A h z z | | z | 2 ( m 2 1 ) 2 l 2 ,   and a j k ( P ( x , A h ξ ) ; x , ξ ) a j k ( P ( z , A h ξ ) ; z , ξ ) 0 .   Then | ( z x ) ξ | c m 2 2 k | ξ |   .
Proof of claim  3.2.4  The identity z = z A h ξ | ξ | A h ξ + P ( z , A h ξ )   and Pythagoras' theorem give | ξ A h z | = | z A h ξ | = | z | | ξ | 1 P ( z , A h ξ ) | z | 2 | z | | ξ | ( 1 2 2 l 2 20 ) .   Also
ξ = ξ A h z | z | 2 A h z + P ( ξ , A h z ) , (3.18)
so | P ( ξ , A h z ) | 2 = | ξ | 2 | ξ A h z | 2 | z | 2 | ξ | 2 | ξ | 2 ( 1 2 2 l 2 20 ) 2 | ξ | 2 2 2 l 2 19 .   Using again ( 3.18 ) z ξ = | z | 2 ( ξ A h z ) ( A h z z ) + z P ( ξ , A h z ) ,   and therefore
| z ξ | | z | | ξ | ( 1 2 2 l 2 20 ) ( m 2 1 ) 2 l 2 | z | | ξ | 2 l 2 9
2 k + 10 | ξ | ( 1 2 2 l 2 20 ) ( m 2 1 ) 2 k + 2 | ξ | 2 k + 9 | ξ | m 2 .
Another application of ( 3.18 ) gives A h ξ ξ = | z | 2 ( ξ A h z ) ( A h ξ A h z ) + A h ξ P ( ξ , A h z ) .   Here A h ξ A h z = ξ z   , so
| A h ξ ξ | | z | 2 | z | | ξ | ( 1 2 2 l 2 20 ) 2 k + 9 | ξ | m 2 | ξ | 2 2 l 2 9
| ξ | 2 ( 1 2 2 l 2 20 ) 2 l 2 2 m 2 | ξ | 2 2 l 2 9 2 l 2 3 m 2 | ξ | 2 .
Finally, ( z x ) ξ = ( z x ) A h ξ | ξ | 2 ( A h ξ ξ ) + ξ P ( z , A h ξ ) ξ P ( x , A h ξ ) .   Here | ( z x ) A h ξ | | z A h ξ | | x | | ξ | | z | | ξ | ( 1 2 2 l 2 20 ) | x | | ξ | 1 2 | z | | ξ | ,   so | ( z x ) ξ | 1 2 | z | | ξ | m 2 2 l 2 3 2 k + 1 | ξ | m 2 2 k + 6 | ξ | 2 k + 1 | ξ | m 2 2 k + 5 | ξ | .   This proves claim  3.2.4 .
Claim 3.2.5 : Suppose l 2 l 1 + 11   , 10 m 2 2 l 2   and | A h z z | | z | 2 ( m 2 1 ) 2 l 2 .   Then | K l 1 l 2 ( x , z ) | c m 2 2 2 ( n 2 ) j 2 k ( n 1 ) l 2 max N 2 ( ξ ξ ) N a j k L 2 , x R n .  
Proof of claim  3.2.5  For x A l 1   and z A l 2   (otherwise K l 1 l 2 ( x , z ) = 0   ), it follows from claim  3.2.4 and parts (ii) and (iii) of Proposition  3.1.3 that
| K l 1 l 2 ( x , z ) | = | e i ( x z ) ξ ( ( 1 i ( x z ) ξ ξ ξ ) 2 ) t |
[ a j k ( P ( x , A h ξ ) ; x , ξ ) a ¯ j k ( P ( z , A h ξ ) ; z , ξ ) χ 2 ( | ξ | ) φ 2 ( δ 2 j 1 | ξ | ) ] d ξ
c m 2 2 2 ( n 2 ) j 2 k ( n 1 ) l 2 max N 2 ( ξ ξ ) N a j k L 2 ,
since the integrand is 0   for ξ   not in the cone with angle 2 l 2 9   (see the proof of claim  3.2.1 ). This proves claim  3.2.5 .
With these results Ψ l 1 Ψ l 2 * u L 2   can be estimated in the remaining case. Case 2: l 2 l 1 + 11   .
Ψ l 1 Ψ l 2 * u L 2 2 m 1 S m 1 | A l 2 K l 1 l 2 ( x , z ) u ( z ) d z | 2 d x .   Writing A l 2 = S l 2 * + m 2 = 1 9 S m 2 + m 2 = 10 2 l 2 1 S m 2 ,   and then Ψ l 1 Ψ l 2 * u L 2 11 ( I 1 + I 2 + I 3 ) ,   where
I 1 = m 1 S m 1 ( S l 2 * | K l 1 l 2 ( x , z ) u ( z ) | d z ) 2 d x ,
I 2 = m 1 S m 1 m 2 = 1 9 ( S m 2 | K l 1 l 2 ( x , z ) u ( z ) | d z ) 2 d x ,
I 3 = m 1 S m 1 | m 2 = 10 2 l 2 1 S m 2 | K l 1 l 2 ( x , z ) u ( z ) | d z | 2 d x .
Using claims  3.2.1 and  3.2.5 with m 2 = 2 l 2 1   (integration by parts),
| I 1 | c 2 n k + l 1 2 n k + n l 2 ( n 1 ) l 1 2 2 ( n 2 ) j 4 k 2 l 2 ( n 1 ) 4 l 2
max N 2 ( ξ ξ ) N a j k L 4 m 1 χ S m 1 * u L 2 2
c 2 2 ( n 2 ) j 2 k ( 2 n 4 ) 2 ( 2 n ) l 1 2 ( 2 n ) l 2 max N 2 ( ξ ξ ) N a j k L 4 u L 2 2 .
Using claims  3.2.2 and  3.2.3 (size)
| I 2 | c 2 n k + l 1 2 n k + ( n 1 ) l 2 2 2 n j 2 l 2 ( n 1 ) a j k L 4 m 1 χ S m 1 * u L 2 2
c 2 2 n ( k + j ) 2 l 1 l 2 ( n 1 ) a j k L 4 u L 2 2 .
Using claim  3.2.3 (size) and  3.2.5 with m 2 { 10 , 11 , . . . 2 l 2 1 }   (integration by parts), and the Cauchy-Schwarz inequality on the sum in m 2   ,
| I 3 | c 2 n k + l 1 ( m 2 = 10 2 l 2 1 1 m 2 2 ) ( m 2 = 10 2 l 2 1 m 2 2 2 n k + l 2 ( n 1 ) m 2 4 )
2 2 ( n 2 ) j 4 k 2 l 2 ( n 1 ) max N 2 ( ξ ξ ) N a j k L 4 m 1 χ S m 1 * u L 2 2
c 2 2 ( n 2 ) j 2 k ( 2 n 4 ) 2 l 1 l 2 ( n 1 ) max N 2 ( ξ ξ ) N a j k L 4 u L 2 2 .
Combining the above estimates for I 1 , I 2   and I 3   , one obtains for l 1 + 11 l 2  
Ψ l 1 Ψ l 2 * L 2 c 2 n ( k + j ) max N 2 ( ξ ξ ) N a j k L 2 2 ( l 1 l 2 ) / 2 2 l 2 n 2 2 . (3.19)
By ( 3.17 ), ( 3.19 ) is valid for all l 1 l 2   . Notice that from ( 3.17 ) and ( 3.19 ) we conclude that the estimates are much better for n > 2   than for n = 2   . In fact it is just in dimension n = 2   where Cotlar-Stein lemma is necessary.
Taking γ ( l ) = 2 l / 4   it follows from Cotlar-Stein lemma that l = 0 Ψ l u   converges in L 2   for any u S ( R n )   . As it was previously noted, l = 0 Ψ l u   converges to ( 1 χ B 2 k + 10 ) Ψ j k u   in the distribution sense, so ( 1 χ B 2 k + 10 ) Ψ j k u L 2 c 2 n ( k + j ) / 2 max N 2 ( ξ ξ ) N a j k L 2 u L 2 , u S ( R n ) .   Combining with ( 3.16 ) Ψ j k u c 2 n ( k + j ) / 2 max N 2 ( ξ ξ ) N a j k L 2 u L 2 .   Finally, one can use the following fact which is easily verified using standard Fourier transform arguments.
Lemma 3.2.1 For any N 1 , N 2 N   there exists N 3 N   such that 2 N 1 ( j + k ) ( ξ ξ ) N 2 a j k L c N 1 max k + | α | N 3 s k ( ξ ξ ) N 2 s α a L .  
Using Lemma  3.2.1 with N 1 n / 2 + 1   one obtains I I L x 2 max k + | α | N 1 ; N 2 s k ( ξ ξ ) 1 N s α a L u L x 2 ,   where I I   was defined in ( 3.15 ), and the proof of Theorem  3.2.1 is completed.

3.3 Composition results

First a few facts concerning oscillatory integrals will be listed. For details and further results, see [25, section 1.6.
Definition 3.3.1 Let m R   , τ 0   . Then A τ m   is the class of functions φ C ( R y n × R ξ n )   satisfying | y α ξ β φ ( y , ξ ) | c α , β ξ m y τ .   The class A   of amplitude functions is defined by A = m R τ 0 A τ m .  
Definition 3.3.2 Let a A   . Then O s e i y ξ a ( y , ξ ) d y d ξ = lim ε 0 e i y ξ a ( y , ξ ) χ ~ ( ε y , ε ξ ) d y d ξ ,   if χ ~ S ( R y n × R ξ n )   and χ ~ ( 0 , 0 ) = 1   .
The oscillatory integral in Definition  3.3.2 is well-defined because of the following lemma which allows one to integrate by parts and use Lebesgue's dominated convergence theorem.
Lemma 3.3.1 Let χ ~ S ( R y n )   with χ ~ ( 0 ) = 1   . Then
( i ) χ ~ ( ε y ) 1 in R n uniformly on compact sets ,
( i i ) y α [ χ ~ ( ε y ) ] 0 in R n uniformly for α N n { 0 } ,
( i i i ) α N n c α > 0 s.t.
| y α [ χ ~ ( ε y ) ] | c α ε σ y σ | α | , y R n , 0 σ | α | , 0 < ε < 1 .
Theorem 3.3.1 Let a ( s ; x , ξ ) S ( R n ; S 1 , 0 m )   , α N n   and φ S ( R n )   . Suppose N m + | α |   , N N   .
Let
c 1 ( x , ξ ) = | β | < N i | β | β ! φ ( x ) ξ β [ ( i ξ ) α ] x β b h ( x , ξ ) , (3.20)
and
c 2 ( x , ξ ) = | β | < N i | β | β ! ξ β b h ( x , ξ ) x β φ ( x ) ( i ξ ) α . (3.21)
Let
E 1 = φ x α Ψ b h Ψ c 1 , (3.22)
and
E 2 = Ψ b h φ x α Ψ c 2 . (3.23)
Then there exist N 1 N   and c > 0   such that for any u S ( R n )   and for j = 1 , 2  
E j u L 2 c max | α | + | β | N 1 x α x β φ L max | γ 1 + γ 2 + γ 3 + γ 4 | N 1 s γ 4 ξ m + | γ 3 | s γ 1 x γ 2 ξ γ 3 a L u L 2 . (3.24)
Proof of Theorem  3.3.1 For simplicity of the exposition we shall drop all the powers of 2 π   which appear in the definition of the operators Ψ b   .
Let u S ( R n )   . Then
E 1 u ( x ) = φ ( x ) x α e i x ξ b h ( x , ξ ) u ^ ( ξ ) d ξ
φ ( x ) e i x ξ | β | < N i | β | β ! ξ β [ ( i ξ ) α ] x β b h ( x , ξ ) u ^ ( ξ ) d ξ .
Notice that x α [ e i x ξ b h ( x , ξ ) ] = β α ( α β ) i | α β | ξ α β e i x ξ x β b h ( x , ξ ) ,   and that if β α   i | α | β ! ξ β [ ( i ξ ) α ] = i | α β | β ! α ! ( α β ) ! ξ α β = i | α β | ( α β ) ξ α β ,   then E 1 u ( x ) = φ ( x ) e i x ξ N | β | , β α ( α β ) i | α β | ξ α β x β b h ( x , ξ ) u ^ ( ξ ) d ξ .   If m > 0   , then E 1 = 0   .
If m 0   , N | β |   , and β α   , then φ ( x ) ξ α β x β b h ( x , ξ ) S 1 , 0 0 .   Hence, there exists c > 0   and N 1 N   such that ( 3.24 ) holds for j = 1   .Let us go back to E 2   . We shall assume that φ C 0 ( R n )   and vanish outside B R ( 0 )   . Then we shall obtain bounds of the form of powers of R   . Thus by introducing a partition of unity with respect to dyadic x   -annuli and summing the corresponding operators one obtains the case φ S ( R n )   .
Write E 2 u ( x ) = ( 1 χ ( ( 2 R ) 1 | x | ) ) E 2 u ( x ) + χ ( ( 2 R ) 1 | x | ) E 2 u ( x ) = I + I I .   To estimate I   it suffices to consider classical symbols. Indeed, integrating by parts with respect to y   in the second term
I = ( 1 χ ( ( 2 R ) 1 | x | ) ) O s ( e i ( x y ) ξ b h ( x , ξ ) φ ( y ) y α u ( y ) e i ( x y ) ξ | β | < N i | β | β ! ξ β b h ( x , ξ ) x β φ ( x ) y α u ( y ) ) d y d ξ . (3.25)
By Taylor's formula
φ ( y ) = | β | < N 1 β ! x β φ ( x ) ( y x ) β
+ N | β | = N 1 β ! ( y x ) β 0 1 x β φ ( θ y + ( 1 θ ) x ) ( 1 θ ) N 1 d θ .
Using ( y x ) β e i ( x y ) ξ = i | β | ξ β [ e i ( x y ) ξ ] ,   and integrating by parts with respect to ξ   in the first term of I   in  3.25 
I = N ( 1 χ ( ( 2 R ) 1 | x | ) ) 0 1 O s e i ( x y ) ξ | β | = N i | β | β ! ξ β b h ( x , ξ )
x β φ ( θ y + ( 1 θ ) x ) ( 1 θ ) N 1 y α u ( y ) d y d ξ d θ .
For each θ [ 0 , 1 ]   , the multiple symbol N ( 1 χ ( ( 2 R ) 1 | x | ) ) | β | = N i | β | β ! ξ β b h ( x , ξ ) x β φ ( θ y + ( 1 θ ) x ) ( 1 θ ) N 1 ,   is in S 1 , 0 m N   uniformly in θ   . Since m N + | α | 0   , there exist c > 0   and N 1 N   such that ( 3.24 ) holds with E j u   replaced by I. In II, note that χ ( ( 2 R ) 1 | x | ) x β φ ( x ) = 0   , so that
I I = ( 1 ) | α | χ ( ( 2 R ) 1 | x | )
O s e i ( x y ) ξ β α ( α β ) ( i ξ ) β b h ( x , ξ ) y α β φ ( y ) u ( y ) d y d ξ ,
by integration by parts with respect to y   . Therefore, it suffices to consider the pseudo-differential operator Ψ I I   given by
Ψ I I u ( x ) = χ ( ( 2 R ) 1 | x | )
O s e i ( x y ) ξ ξ β a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) φ ( y ) u ( y ) d y d ξ ,
for β α   . Choose M N   such that 2 M m | α | n + 1   . Integrating by parts with respect to ξ   and using | x y | | x | / 2 R > 0   ,
Ψ I I u ( x ) = χ ( ( 2 R ) 1 | x | ) e i ( x y ) ξ ( 1 ) M | x y | 2 M Δ ξ M [ ξ β a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) ] φ ( y ) u ( y ) d y d ξ , (3.26)
where the integral converges absolutely. Indeed,
| Δ ξ M [ ξ β a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) ] | c x 2 M ξ n 1 max | γ 1 + γ 2 | 2 M ξ m + | γ 2 | s γ 1 ξ γ 2 a L .
Therefore, letting | a | ( 2 M ) = max | γ 1 + γ 2 | 2 M ξ m + | γ 2 | s γ 1 ξ γ 2 a L ,   one has | Ψ I I u ( x ) | c R n / 2 φ L | a | ( 2 M ) u L 2 .   Define χ R ( | x | ) = χ ( | x | / R )   . Then by Cauchy-Schwarz ( 1 χ c 0 R ) Ψ I I u L 2 c R n φ L | a | ( 2 M ) u L 2 ,   where the constant c 0   will be fixed below (see ( 3.29 )). Now we turn to the estimate of
χ c 0 R Ψ I I u L 2 . (3.27)
We shall use that
1 | x y | 2 M = 1 | | x | 2 2 ( x y ) + | y | 2 | M = 1 | x | 2 M 1 | 1 2 ( x | x | y | x | ) + | y | 2 | x | 2 | M = 1 | x | 2 M ( 1 j = 1 ( 2 n + 2 ) M 1 P j ( x | x | ) Q j ( y | x | ) ) 1 , (3.28)
where P j , Q j   are monomials of degree no bigger than 2 M   with deg Q j 0   for j = 1 , 2 , . . . , ( 2 n + 2 ) M 1   . Since | x | > c 0 R   and | y | < R   it follows that | P j ( x | x | ) | a j , | Q j ( y | x | ) | b j ( | y | | x | ) d e g Q j b j c 0 d e g Q j ,   so | P j ( x | x | ) Q j ( y | x | ) | a j b j c 0 d e g Q j .   We take
c 0 2 ( 2 n + 2 ) M max { | a j b j | : j = 1 , . . , ( 2 n + 2 ) M 1 } (3.29)
and rewrite ( 3.28 ) as
1 | x | 2 M 1 1 j P j ( x | x | ) Q j ( y | x | ) = 1 | x | 2 M ( k = 0 ( j = 1 ( 2 n + 2 ) M 1 P j ( x | x | ) Q j ( y | x | ) ) k ) = 1 | x | 2 M k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 P j i ( x | x | ) Q j i ( y | x | ) | x | deg Q j i | x | deg Q j i . (3.30)
We observe that | x | deg Q j i Q j i ( y | x | )   depends just on y   and if | y | < R   then | | x | deg Q j i Q j i ( y | x | ) | b j i R deg Q j i .   So returning to the estimate ( 3.27 ) and using the argument in ( 3.26 ) it follows
χ c 0 R Ψ I I u L 2 = χ c 0 R k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 ( 1 ) M P j i ( x | x | ) | x | deg Q j i
e i ( x y ) ξ 1 | x | 2 M Δ ξ M ( ξ β a ( ; , ) χ ( | ξ | ) ) Q j i ( y | x | ) | x | deg Q j i φ ( y ) u ( y ) d y d ξ L 2 .
Using that b ~ ( x , ξ ) = χ c 0 R ( | x | ) 1 | x | 2 M Δ ξ M ( ξ β a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) )   is a symbol which falls under the scope of Theorem  3.2.1 -see Remark  3.1.1 (c), it follows that
e i x ξ b ~ ( x , ξ ) ( e i y ξ Q j i ( y ) φ ( y ) u ( y ) d y ) d ξ L x 2
c e i y ξ Q j i ( y ) φ ( y ) u ( y ) d y L ξ 2
c b j i R deg Q j i u L 2 .
Also
χ c 0 R k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 ( 1 ) M P j i ( x | x | ) | x | deg Q j i L k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 a j i ( c 0 R ) deg Q j i .
Combining the last two estimates we conclude that
χ c 0 R Ψ I I u L 2
c k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 a j i b j i R deg Q j i ( c 0 R ) deg Q j i u L 2 c k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 a j i b j i c 0 deg Q j i u L 2
c k = 0 ( j = 1 ( 2 n + 2 ) M 1 a j i b j i c 0 ) k u L 2 k = 0 1 2 k u L 2 c u L 2
and the proof of Theorem  3.3.1 is completed.
Theorem 3.3.2 Let a S ( R n : S 1 , 0 m )   and let b h   be defined as ( 3.9 ), N N   , N m   . Let c ( x , ξ ) = | α | < N i | α | α ! x α ξ α b ¯ h ( x , ξ )   and E = Ψ b h * Ψ c   . Let φ S ( R n )   . Then there exists N 1 = N 1 ( n ) N   such that for u S ( R n )  
φ E u L 2 c max | β | N 1 x N 1 x α φ L max β 4 + | β 1 + β 2 + β 3 | N 1 s β 4 ξ m + | β 3 | s β 1 x β 2 ξ β 3 a L u L 2 (3.31)
Proof of Theorem  3.3.2 We shall prove the Theorem for φ C 0 ( R n )   with support in B R ( 0 )   . Introducing a partition of unity and summing the corresponding terms we get the desired result.
Thus,
φ ( x ) E u ( x ) = O s e i ( x y ) ξ φ ( x ) { b ¯ h ( y , ξ ) | α | < N i | α | α ! x α ξ α b ¯ h ( x , ξ ) } u ( y ) d y d ξ
= O s e i ( x y ) ξ φ ( x ) ( 1 χ ( ( 2 R ) 1 | y | ) )
{ b ¯ h ( y , ξ ) | α | < N i | α | α ! x α ξ α b ¯ h ( x , ξ ) } u ( y ) d y d ξ
+ O s e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | )
{ b ¯ h ( y , ξ ) | α | < N i | α | α ! x α ξ α b ¯ h ( x , ξ ) } u ( y ) d y d ξ
= I + I I .
First consider I   . By a Taylor expansion of order N   ,
b h ( y , ξ ) = | α | < N 1 α ! ( y x ) α x α b h ( x , ξ )
+ | α | = N N α ! ( y x ) α 0 1 ( 1 θ ) N 1 x α b h ( θ y + ( 1 θ ) x , ξ ) d θ .
Since ( y x ) α e i ( x y ) ξ = i | α | ξ α e i ( x y ) ξ   , integration by parts with respect to ξ   gives
I = O s e i ( x y ) ξ φ ( x ) ( 1 χ ( ( 2 R ) 1 | y | ) ) N i N
| α | = N 1 α ! 0 1 ( 1 θ ) N 1 x α ξ α b ¯ h ( θ y + ( 1 θ ) x , ξ ) d θ u ( y ) d y d ξ .
Because of the compact support in x   and y   , the multiple symbol φ ( x ) ( 1 χ ( ( 2 R ) 1 | y | ) ) N i N | α | = N 1 α ! ( 1 θ ) N 1 x α ξ α b ¯ h ( θ y + ( 1 θ ) x , ξ )   is in S 1 , 0 m N S 1 , 0 0   uniformly in θ   , so an estimate of the type ( 3.31 ) holds with φ E u   replaced by I   . A factor R N   appears due to the differentiation with respect to ξ   of the quantity P ( θ y + ( 1 θ ) x , A h ξ )   .
Next we consider I I   . Choose M N   such that m 2 M n 1   . If φ ( x ) χ ( ( 2 R ) 1 | y | ) 0 ,   then | x y | R   so by integration by parts with respect to ξ  
I I = e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | ) ( 1 ) M | x y | 2 M Δ ξ M b ¯ h ( y , ξ ) u ( y ) d y d ξ
e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | ) ( 1 ) M | x y | 2 M
| α | < N i | α | α ! x α Δ ξ M ξ α b ¯ h ( x , ξ ) u ( y ) d y d ξ ,
where the integrals converge absolutely by the choice on M   . The multiple symbol in the second term above is in S 1 , 0 m 2 M S 1 , 0 0   because of the factor φ ( x )   , so the corresponding pseudo-differential operator is L 2   -bounded as in ( 3.31 ).
Therefore, it remains to show the L 2   -boundedness of the operator
Ψ ¯ I I u ( x ) = e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | )
( 1 ) M | x y | 2 M Δ ξ M [ a ¯ ( P ( y , A h ξ ) ; y , ξ ) χ ( | ξ | ) ] u ( y ) d y d ξ .
But Ψ ¯ I I   is the adjoint of the operator Ψ I I   defined in ( 3.26 ) in the proof of Theorem  3.3.1 (with β = 0   and φ   replaced by φ ¯   ). Since Ψ I I   was there proved to be bounded in L 2   , so is Ψ ¯ I I   , and the operator norms are equal.
This proves Theorem  3.3.2 .
Theorem 3.3.3 Let a 1 S ( R n : S 1 , 0 m 1 )   , a 2 S ( R n : S 1 , 0 m 2 )   , and let b 1 , b 2   be the corresponding symbols given in ( 3.9 ), Definition  3.1.1 , with A = A h   . Suppose N N   , N m 1 + m 2   and let c ( x , ξ ) = | α | < N i | α | α ! ξ α [ b 1 ( x , ξ ) x α b ¯ 2 ( x , ξ ) ] .   Let φ S ( R n )   . Then there exist c = c ( N )   and N 1 N   such that for any u S ( R n )  
φ ( Ψ b 1 Ψ b 2 * Ψ c ) u L 2 max | β | N 1 x N 1 x β φ L max β 4 + | β 1 + β 2 + β 3 | N 1 s β 4 ξ m 1 + | β 3 | s β 1 x β 2 ξ β 3 a 1 L max β 4 + | β 1 + β 2 + β 3 | N 1 s β 4 ξ m 2 + | β 3 | s β 1 x β 2 ξ β 3 a 2 L u L 2 . (3.32)
Proof of Theorem  3.3.3 As in the proof of Theorem  3.3.2 we shall assume that φ C 0 ( R n )   with support in B R ( 0 )   .
Introducing a partition of unity and summing the corresponding terms we get the desired result.
We have,
φ ( Ψ b 1 Ψ b 2 * Ψ c ) u ( x ) = O s e i ( x y ) ξ φ ( x ) [ b 1 ( x , ξ ) b ¯ 2 ( y , ξ ) c ( x , ξ ) ] u ( y ) d y d ξ
= O s e i ( x y ) ξ φ ( x ) ( 1 χ ( ( 2 R ) 1 | y | ) ) [ b 1 ( x , ξ ) b ¯ 2 ( y , ξ ) c ( x , ξ ) ] u ( y ) d y d ξ
+ O s e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | ) [ b 1 ( x , ξ ) b ¯ 2 ( y , ξ ) c ( x , ξ ) ] u ( y ) d y d ξ
= I + I I .
First we consider I   . By a Taylor expansion of order N   ,
b 2 ( y , ξ ) = | α | < N 1 α ! ( y x ) α x α b 2 ( x , ξ )
+ | α | = N N α ! ( y x ) α 0 1 ( 1 θ ) N 1 x β b 2 ( θ y + ( 1 θ ) x , ξ ) d θ .
Since ( y x ) α e i ( x y ) ξ = i | α | ξ α e i ( x y ) ξ   , integration by parts with respect to ξ   gives
I = e i ( x y ) ξ φ ( x ) ( 1 χ ( ( 2 R ) 1 | y | ) )
| α | = N N α ! i | α | 0 1 ( 1 θ ) N 1 ξ α [ b 1 ( x , ξ ) x α b ¯ 2 ( θ y + ( 1 θ ) x , ξ ) ] d θ u ( y ) d y d ξ .
Because of the compact support in x   and y   , the multiple symbol φ ( x ) ( 1 χ ( ( 2 R ) 1 | y | ) ) | α | = N N α ! i | α | ξ α [ b 1 ( x , ξ ) x α b ¯ 2 ( θ y + ( 1 θ ) x , ξ ) ]   is in S 1 , 0 m 1 + m 2 N S 1 , 0 0   with seminorms uniformly bounded in θ [ 0 , 1 ]   , so the estimate ( 3.32 ) holds with the left hand side replaced by I L x 2   . A factor R N 1   appears from differentiation with respect to ξ   of the quantity P ( θ y + ( 1 θ ) x , A h ξ )   , θ [ 0 , 1 ]   .
Next consider I I   . Note that if φ ( x ) χ ( ( 2 R ) 1 | y | ) 0   , then | x y | R   . Therefore, choosing M N   such that m 1 + m 2 2 M n 1   and integrating by parts with respect to ξ  
I I = e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | ) ( 1 ) M | x y | 2 M Δ ξ M [ b 1 ( x , ξ ) b ¯ 2 ( y , ξ ) ] u ( y ) d y d ξ e i ( x y ) ξ φ ( x ) χ ( ( 2 R ) 1 | y | ) ( 1 ) M | x y | 2 M Δ ξ M c ( x , ξ ) u ( y ) d y d ξ , (3.33)
where the integrals converge absolutely.
The multiple symbol in the second term in ( 3.33 ) is in S 1 , 0 n 1 S 1 , 0 0   , so an estimate of the type( 3.32 ) holds.
The adjoint of the operator corresponding to the first term in ( 3.33 ) has multiple symbol φ ¯ ( y ) χ ( ( 2 R ) 1 | x | ) ( 1 ) M | x y | 2 M Δ ξ M [ b ¯ 1 ( y , ξ ) b 2 ( x , ξ ) ] .   Replacing φ ¯   by φ   and b ¯ 1 ( y , ξ )   by ξ β   one obtains a symbol similar to that of the operator Ψ I I   in ( 3.26 ) in the proof of Theorem  3.3.1 , so we just need to sketch the proof.
Thus, we define
Ψ u ( x ) = e i ( x y ) ξ φ ¯ ( y ) χ ( ( 2 R ) 1 | x | )
( 1 ) M | x y | 2 M ( | γ 1 | + | γ 2 | = 2 M c γ 1 γ 2 ξ γ 1 b ¯ 1 ( y , ξ ) ξ γ 2 b 2 ( x , ξ ) ) u ( y ) d y d ξ .
Next consider φ ~ C 0   which vanishes outside the ball of radius 3 R / 2   and such that φ ~ φ = φ   . Then define b ~ γ 1 ( y , ξ ) = φ ~ ( y ) ξ | γ 1 | m 1 ξ γ 1 b ¯ 1 ( y , ξ ) ,   and b ~ γ 2 ( x , ξ ) = χ ( ( 2 R ) 1 | x | ) | x | 2 M ξ m 1 | γ 1 | ξ γ 2 b 2 ( x , ξ ) .   Using the notation introduced in the proof of Theorem  3.3.1 , see ( 3.28 )-( 3.30 ), we have 1 | x y | 2 M = 1 | x | 2 M k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 k i = 1 P j i ( x | x | ) Q j i ( y | x | ) .   Hence as in ( 3.28 )-( 3.30 ) we can reduce ourselves to study χ ( ( c 0 R ) 1 | x | ) Ψ u L 2   for c 0   large enough.
Then we have
χ ( ( c 0 R ) 1 | x | ) Ψ u L 2
c k = 0 j 1 . . j k = 1 ( 2 n + 2 ) M 1 | γ 1 | + | γ 2 | = 2 M χ ( ( c 0 R ) 1 | x | ) k i = 1 P j i ( x | x | ) 1 | x | deg Q j i
e i ( x y ) ξ φ ( y ) Q j i ( y ) b ~ γ 1 ( y , ξ ) b ~ γ 2 ( x , ξ ) u ( y ) d y d ξ L x 2 .
Then from Remark  3.1.1 (b)
e i ( x y ) ξ φ ( y ) Q j i ( y ) b ~ γ 1 ( y , ξ ) b ~ γ 2 ( x , ξ ) u ( y ) d y d ξ L x 2
= e i x ξ b ~ γ 2 ( x , ξ ) ( F ( y , ξ ) d y ) d ξ L x 2
F ( y , ξ ) d y L ξ 2 ,
with F ( y , ξ ) = e i y ξ φ ( y ) Q j i ( y ) b ~ γ 1 ( y , ξ ) u ( y ) .   Now observe that if a ( z , η ) = b ~ γ 1 ( η , z ) ,   then | z α η β a ( z , η ) | c α β η | β | R | α | + | β | .   Hence a S 1 , 0 0   and therefore F ( y , ξ ) d y L ξ 2 c R deg Q j i + N u L 2   for some large N   which just depends on the dimension.Gathering the above information with the argument used in the proof of Theorem  3.3.1 one completes the proof.

4 THE BICHARACTERISTIC FLOW

In this section we study the bicharacteristic flows associated to the second order ultrahyperbolic variable coefficient operator ( x )   and its truncated version R ( x )   , see ( 4.12 )-( 4.13 ) below. Assuming a non-trapping condition (basic assumption) for the flow generated by ( x )   we prove in subsection  4.1 that the bicharacteristic flow is “uniformly non-trapping” with respect to the parameter R   . The analysis of the flow is moredelicate when the trajectory is not outgoing 1   , -see Theorem  4.1.1 below. In that case we prove that outside a bounded ball (the same for all R   large enough) it behaves in the x   variable as the free flow, but just in dyadic annuli see Theorem  4.1.1 ( v) for a precise statement. On the other hand in the outgoing case, the trajectories are in fact perturbations of the free ones, as is proved by Craig, Kappeler, and Strauss in [4whose arguments we follow. The end of subsection  4.1 is devoted to prove that the non-trapping condition is stable under small perturbation in the coefficients.
In subsection  4.2 we deduce several estimates for the continuous dependence upon the initial data of the flows associated to the operators R ( x )   (with respect the parameter R   ). The arguments in [4are again very helpful. The estimates in subsections  4.1 and  4.2 will be used in the next section to deduce several properties of the integrating factor K R   .
One of the main differences of the flow in the non-elliptic setting with respect to the elliptic one studied in Section  2 is that the Hamiltonian in the elliptic case h 2 ( x , ξ ) = j , k a j k ( x ) ξ k ξ j   is preserved under the flow. So ellipticity gives the a priori estimate -see ( 2.15 ) in Section  2 , ν 2 | ξ 0 | 2 | Ξ ( s ; x 0 , ξ 0 ) | 2 ν 2 | ξ 0 | 2 .   In particular, one has that the bicharacteristic flow is globally defined. This is also true in the ultra-hyperbolic case under the ”asymptotic flatness” assumption, but does not follow immediately.

1   Using the homogeneity property ( 2.16 ) in Section  2 , it is enough to consider the flow for s 0   and | ξ | = 1   .

4.1 Uniformly non-trapping flows

Let ( x ) = x j a j k ( x ) x k ,   where A ( x ) = ( a j k ( x ) )   is a real, symmetric, and non-degenerated matrix, i.e.
ν > 0 ξ R n ν 1 | ξ | | A ( x ) ξ | ν | ξ | . (4.1)
We will assume that there exists a constant coefficient operator
0 = a j k 0 x j x k 2 , (4.2)
with
a j k ( x ) a j k 0 S ( R n ) , j , k = 1 , . . , n . (4.3)
After a change of variable we can assume, without loss of generality, that
A h = ( a j k 0 ) = ( I k × k 0 0 I ( n k ) × ( n k ) ) . (4.4)
We shall use the notation ξ ~ = A h ξ .   The bicharacteristic flow is defined by
{ d d s X j ( s ; x 0 , ξ 0 ) = 2 k = 1 n a j k ( X ( s ; x 0 , ξ 0 ) ) Ξ k ( s ; x 0 , ξ 0 ) , d d s Ξ j ( s ; x 0 , ξ 0 ) = k , l = 1 n x j a k l ( X ( s ; x 0 , ξ 0 ) ) Ξ k ( s ; x 0 , ξ 0 ) Ξ l ( s ; x 0 , ξ 0 ) ( X ( 0 ; x 0 , ξ 0 ) , Ξ ( 0 ; x 0 , ξ 0 ) ) = ( x 0 , ξ 0 ) . (4.5)
From the classical result of ode's we know that the bicharacteristic flow exists in the time interval s ( δ , δ )   with δ = δ ( x 0 , ξ 0 )   , and δ ( )   depending continuously on ( x 0 , ξ 0 )   .
By homogeneity of the symbol of ( x )   (see ( 2.16 ) in Section  2 ) for any t R   one has that X ( s ; , x , t ξ ) = X ( t s ; x , ξ ) , Ξ ( s ; x , t ξ ) = t Ξ ( t s : x , ξ ) .   BASIC ASSUMPTION We shall assume that ( x )   is non-trapping, i.e. for each ( x 0 , ξ 0 ) R n × ( R n { 0 } )   and for each μ > 0   there exists s 0   with 0 < s 0 < δ   such that 2   | X ( s 0 ; x 0 , ξ 0 ) | μ .   Our first result is the following.
Proposition 4.1.1
d d s | X ( s ; x 0 , ξ 0 ) | 2 = 4 X ( s ; x 0 , ξ 0 ) ; A ( X ( s ; x 0 , ξ 0 ) ) Ξ ( s ; x 0 , ξ 0 ) (4.6)
and
d 2 d s 2 | X ( s ; x 0 , ξ 0 ) | 2 = 8 | A ( X ( s ; x 0 , ξ 0 ) ) Ξ ( s ; x 0 , ξ 0 ) | 2 + 8 X ( s ; x 0 , ξ 0 ) ; A ( X ( s ; x 0 , ξ 0 ) ) A ( X ( s ; x 0 , ξ 0 ) ) Ξ ( s ; x 0 , ξ 0 ) Ξ ( s ; x 0 , ξ 0 ) 4 X ( s ; x 0 , ξ 0 ) ; A ( X ( s ; x 0 , ξ 0 ) ) A ( X ( s ; x 0 , ξ 0 ) ) Ξ ( s ; x 0 , ξ 0 ) Ξ ( s ; x 0 , ξ 0 ) . (4.7)
The proof of Proposition  4.1.1 follows directly from ( 4.5 ).
Lemma 4.1.1 There exists M > 0   which depends only on a finite number of S ( R n )   -seminorms of A ( x ) A h   and on ν   in ( 4.1 ) such that if | X ( s ; x 0 , ξ 0 ) | > M   , then
d 2 d s 2 | X ( s ; x 0 , ξ 0 ) | 2 4 | A ( X ( s ; x 0 , ξ 0 ) ) Ξ ( s ; x 0 , ξ 0 ) | 2 0 . (4.8)
The proof of Lemma  4.1.1 follows by combining ( 4.7 ) with ( 4.1 ) and ( 4.3 ).
Lemma 4.1.2 For any M ~ M   , M   as in Lemma  4.1.1 , there exist s M ~   , b   , b ~   , and c 1   depending only on M ~   and A   , such that for | x 0 | M ~   , and 1 / 2 | ξ 0 | 7 / 4   , (i) s M ~ ( 0 , δ ( x 0 , ξ 0 ) )   ; (ii) d d s | X ( s ; x 0 , ξ 0 ) | 2 | s = s M ~ b > 0 ;   (iii) d 2 d s 2 | X ( s ; x 0 , ξ 0 ) | 2 | s = s M ~ b ~ > 0 ;   (iv) for all s   with | s | s M ~   , 0 < c 1 1 | Ξ ( s ; x 0 , ξ 0 ) | c 1 < + ;   (v) for s ( s M ~ , δ )   , (ii) holds.
Proof of Lemma  4.1.2 For convenience we introduce the notation
N ( s ; x 0 , ξ 0 ) = | X ( s ; x 0 , ξ 0 ) | 2 . (4.9)
Because of the homogeneity property of the flow we can assume s > 0   and s 0 > 0   . We fix ( x 0 , ξ 0 ) R n × ( R n { 0 } )   , and let M 1 = max { M ; | x 0 | } ,   where M   is as in Lemma  4.1.1 . By the non-trapping assumption s ¯ = s ¯ ( x 0 , ξ 0 ) s . t . N ( s ¯ ; x 0 , ξ 0 ) ( M 1 + 1 ) 2 .   Also N ( 0 ; x 0 , ξ 0 ) = | x 0 | 2 M 1 2   . Define s 1 = sup { s [ 0 , s ¯ ] : N ( s ; x 0 , ξ 0 ) M 1 2 } .   Thus, s 1 [ 0 , s ¯ ) , N ( s 1 ; x 0 , ξ 0 ) = M 1 2   , and for s ( s 1 , s ¯ ]   one has N ( s ; x 0 , ξ 0 ) > M 1 2   . By the mean value theorem there exists s 2 ( s 1 , s ¯ )   such that N ( s 2 ; x 0 , ξ 0 ) > 0   . Moreover, N ( s 2 ; x 0 , ξ 0 ) > M 1 2 M 2   . By continuity there exists a neighborhood U   of ( x 0 , ξ 0 )   such that
N ( s 2 ; x ¯ , ξ ¯ ) N ( s 2 ; x 0 , ξ 0 ) / 2 > 0 , N ( s 2 ; x ¯ 0 , ξ ¯ 0 ) > M 2 , ( x ¯ , ξ ¯ ) U . (4.10)
Claim 4.1.1 :
N ( s ; x ¯ , ξ ¯ ) N ( s 2 ; x ¯ , ξ ¯ ) , s s 2 ( x ¯ , ξ ¯ ) U . (4.11)
Suppose not, so there exists s 3 > s 2 , ( x , ξ ) U   , and β > 0   such that N ( s 3 ; x , ξ ) < β < N ( s 2 ; x , ξ ) .   Let s 4 = inf { s [ s 2 , s 3 ] : N ( s ; x , ξ ) < β } ,   then s 4 ( s 2 , s 3 ]   , and N ( s 4 ; x , ξ ) = β   . By the mean value theorem, there exists s 5 ( s 2 , s 4 )   such that N ( s 5 ; x , ξ ) < 0 .   But for s [ s 2 , s 4 ) , N ( s ; x , ξ ) β   , and so N ( s 5 ; x , ξ ) > N ( s 2 ; x , ξ ) M 2   , which contradicts Lemma  4.1.1 , so we have established Claim  4.1.1 .
Next we cover K   by finitely many neighborhoods U r , r = 1 , . . , N   , and let a = max r = 1 , . . , N s 2 ( U r ) , ( see (eq4.1.104.10)-(eq4.1.114.11)) ,   and b = min r = 1 , . . , N { N ( s 2 ( U r ) ; x ¯ , ξ ¯ ) : ( x ¯ , ξ ¯ ) U r } .   From ( 4.10 ) it follows that b > 0   and from Claim  4.1.1 one has for s a   that N ( s ; x , ξ ) b , ( x , ξ ) K .   Hence only (iv) remains to be proved. The upper bound follows by compactness. For the lower bound it suffices to see that Ξ ( s ; x 0 , ξ 0 ) 0   for 0 < s < s 0 .   But this is a consequence of the uniqueness of the flow ( 4.5 ).
The proof of Lemma  4.1.2 is completed.
Let θ C 0 ( R n ) , 0 θ ( x ) 1   , s u p p θ { x : | x | 2 }   and θ ( x ) = 1   on { x : | x | 1 }   . For each 1 < R   define
A R ( x ) = θ ( x / R ) A ( x ) + ( 1 θ ( x / R ) ) A h = A h + θ ( x / R ) ( A ( x ) A h ) = a j k R ( x ) , if R > 1 , A = A ; (4.12)
R ( x ) = x j a j k R ( x ) x k if R > 1 , = . (4.13)
Notice that there exists R > 1   such that for R > R   one has
ν 1 | ξ | / 2 | A R ( x ) ξ | 2 ν | ξ | and A R ( x ) A h S ( R n ) (4.14)
uniformly in R, so that Lemma  4.1.1 will apply uniformly in R > R   .
Also
x l a j k R ( x ) = θ ( x / R ) x l a j k ( x ) + x l θ ( x / R ) R ( a j k ( x ) a j k 0 ) , j , k = 1 , . . , n , (4.15)
so for any τ > 0   there exists c τ   such that
| x l a j k R ( x ) | c τ ( 1 + | x | ) τ , j , k = 1 , . . , n . (4.16)
We will always consider R R * = max { R ; 4 M ; 4 }   with M   as in Lemma  4.1.1 .
Denote by ( X R ( s ; x 0 , ξ 0 ) , Ξ R ( s ; x 0 , ξ 0 ) )   the bicharacteristic flow associated to the operators R   .
Theorem 4.1.1 There exist c 1 , c 2 , c 3 > 0   and M 1 > 0   (sufficiently large) such that for any ( x , ξ ) R n × S n 1   and R > 2 M 1   one has (i) the bicharacteristic flow
( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) ) exists for any s R , (4.17)
(ii)
0 < c 1 1 | Ξ R ( s ; x , ξ ) | c 1 , s R , (4.18)
(iii) there exists s ~ 0 R   such that 3  
d d t | X R ( s ; x , ξ ) | 2 | s = s ~ 0 0 and | X R ( s ~ 0 ; x , ξ ) | M 1 , (4.19)
(iv) define  4.19  s 0 = inf { s ~ 0 0 : ( ) holds } ,   then
| X R ( s ; x , ξ ) | 2 c 2 ( s s 0 ) 2 + M 1 2 , s s 0 , (4.20)
(v) for 0 < s s 0   and k { 1 , 2 , . . }   the set
I k R = { s [ 0 , s 0 ] : 2 k | X R ( s ; x , ξ ) | 2 k + 1 } , (4.21)
verifies
| I k R | c 3 2 k . (4.22)
Proof of Theorem  4.1.1 Recall ( 4.5 ), ( 4.14 ), and ( 4.16 ). Then given any τ > 1   we can always assume, possibly by considering another τ ~   with 1 < τ ~ < τ   , that there is M 1 = 2 J   , J N   , 2 J ( 1 τ ) < 1 / 4   such that if | X R ( s ~ 0 ; x , ξ ) | > M 1 = 2 J   then
d 2 d s 2 | X R ( s ~ 0 ; x , ξ ) | 2 4 ν 2 | Ξ R ( s ~ 0 ; x , ξ ) | 2 , (4.23)
and
| d d s | Ξ R ( s ~ 0 ; x , ξ ) | | ν 1 100 | X R ( s ~ 0 ; x , ξ ) | τ | Ξ R ( s ~ 0 ; x , ξ ) | 2 . (4.24)
Assume as a first step that | x | M 1   and 1 / 2 | ξ | 7 / 4   . Since R M 1   then X R ( s ; x , ξ ) = X ( s ; x , ξ )   , Ξ R ( s ; x , ξ ) = Ξ ( s ; x , ξ )   as far as | X ( s ; x , ξ ) | R   . Then from Lemma  4.1.2 there exist c 0 > 0   and s M 1   such that
| X R ( s M 1 ; x , ξ ) | M 1 , (4.25)
c 0 1 | Ξ R ( s ; x , ξ ) | c 0 , if 0 s s M 1 , (4.26)
and
d d s | X R ( s ; x , ξ ) | 2 > b > 0 , if s > s M 1 , (4.27)
as far as | X ( s ; x , ξ | R .   We will also use the inequalities -see ( 4.6 ),
| d d s | X R ( s ; x , ξ ) | 2 | 4 ν | X R ( s ; x , ξ ) | | Ξ R ( s ; x , ξ ) | , (4.28)
and
1 w < 1 1 + w < 1 1 w < 1 + 3 w / 2 , if w ( 0 , 1 / 4 ) . (4.29)
Assume now for j Z   that 2 j + J | x | 2 j + 1 + J   and d d s | X R ( s ; x , ξ ) | 2 | s = 0 < 0 .   If j < 0   we have already proved parts (ii), (iii) and (v) of the theorem. Thus, consider j = 0 , 1 , 2 , . .   From ( 4.24 ) and as long as | X R ( t ; x , ξ ) | > 2 j + J   for t [ 0 , s ]   one has
| ξ | 1 + 2 ( j + J ) τ | ξ | ( s 0 ) / 100 ν | Ξ R ( s ; x , ξ ) | | ξ | 1 2 ( j + J ) τ | ξ | ( s 0 ) / 100 ν . (4.30)
Thus, combining ( 4.29 ) and ( 4.30 ) and if s < 2 ( j + J ) τ 25 ν   ,
| ξ | ( 1 2 ( j + J ) τ | ξ | s 100 ν ) | Ξ R ( s ; x , ξ ) | | ξ | ( 1 + 3 2 2 ( j + J ) τ | ξ | s 100 ν ) . (4.31)
So, in particular,
1 2 | ξ | | Ξ R ( s ; x , ξ ) | 7 4 | ξ | , (4.32)
as long as | X R ( t ; x , ξ ) | > 2 j + J   for t [ 0 , s ]   .
From ( 4.23 ), ( 4.28 ), and ( 4.32 ) we get (recall ξ S n 1   ) that
d d s | X R ( s ; x , ξ ) | 2 7 ν 2 j + J + 1 + 2 ν 1 s , (4.33)
as long as | X R ( t ; x , ξ ) | > 2 j + J   for t [ 0 , s ]   . Then either there is s ~ j + 1   , 0 s ~ j + 1 7 ν 2 2 j + J + 1   such that | X R ( s ~ j + 1 ; x , ξ | 2 j + J   or for 0 s 7 ν 2 2 j + J + 1   , | X R ( s ; x , ξ | 2 j + J   . In the second case and from ( 4.33 ) there exists s ~ 0   such that ( 4.19 ) holds with 0 s ~ 0 7 ν 2 2 j + J + 1   . Therefore s 0   exists and from ( 4.28 ) and ( 4.32 ) we have that | X ( s 0 ; x , ξ ) | ( 7 ν 3 + 1 ) 2 j + J   . Therefore (v) holds.
Assume now the first case. Define s j + 1 = inf { s : 0 s 7 ν 2 2 j + J + 1 such that | X R ( s ; x , ξ ) | 2 j + J } .   Hence 2 j + J | X R ( s ; x , ξ ) | ( 7 ν 3 + 1 ) 2 j + J + 1   for 0 s s j + 1   . Then from ( 4.31 ) we get | ξ | ( 1 2 ( j + J ) ( 1 τ ) ) | Ξ R ( s j + 1 ; x , ξ ) | | ξ | ( 1 + 2 ( j + J ) ( 1 τ ) ) .   We repeat the process changing j   by j 1   and ξ   by Ξ R ( s j + 1 ; x , ξ ) = : ξ j   , and taking s > s j + 1   . Then as before either the trajectory becomes outgoing at some s ~ 0 > s j + 1 + 7 ν 2 j + J   and we obtain (v), or there is s j   such that | X R ( s j ; x , ξ ) | 2 j 1 + J and s j s j + 1 7 ν 2 2 j + J ,   and | ξ j | ( 1 2 ( j 1 + J ) ( 1 τ ) ) | Ξ R ( s ; x , ξ ) | | ξ j | ( 1 + 2 ( j 1 + J ) ( 1 τ ) ) ,   for s j + 1 < s < s j   . We keep doing this until either the trajectory becomes outgoing or there are l = j , j 1 , . . . . , 1 ,   with | X R ( s l ; x , ξ ) | 2 l 1 + J and s l s l + 1 7 ν 2 2 l + J ,   and | ξ l | ( 1 2 ( l 1 + J ) ( 1 τ ) ) | Ξ R ( s ; x , ξ ) | | ξ l | ( 1 + 2 ( l 1 + J ) ( 1 τ ) ) .   Here ξ l = Ξ R ( s l + 1 ; x , ξ )   and s l + 1 < s < s l   .Assume this second possibility holds. From the condition 2 J ( 1 τ ) < 1 / 4   it follows that
j l = 1 ( 1 ± 2 ( l + J ) ( 1 τ ) ) = exp ( l = 1 j ln ( 1 ± 2 ( l + J ) ( 1 τ ) ) )
Hence, we get 1 / 2 | Ξ R ( s j ; x , ξ ) | 7 / 4 .   Then we can apply the first step and obtain ( 4.25 )-( 4.27 ). In particular from ( 4.25 ) we know that s 0 s 1 + s M 1   , and from ( 4.26 ) c 0 1 | Ξ R ( s 0 ; x , ξ ) | c 0 .   At this point we have proved (iii). By considering the cases k + 1 J   and k + 1 > J   (v) also holds assuming either d d s | X R ( s , x , ξ ) | 2 | s = 0 < 0   or | x | M 1   . Otherwise s 0 = 0   and (v) is void.
It remains to prove (i), (ii), and (iv) for s s 0   . From ( 4.23 ) we know that | X R ( s , x , ξ ) |   is nondecreasing for s > s 0   . Assume that | X R ( s 0 , x , ξ ) | 2 L + J ,   for some L = 0 , 1 , . . .   . Define ξ 0 = Ξ R ( s 0 ; x , ξ ) .   We know that min { c 0 1 , 1 / 2 } | ξ 0 | max { c 0 , 7 / 4 } .   Following the same argument as in ( 4.31 ) we know that ( X R ( s ; x , ξ ) ; Ξ R ( s ; x , ξ ) )   is defined for s s 0 + 25 ν | ξ 0 | 2 J + L = s 1   and from ( 4.23 ) | X R ( s 0 ; x , ξ ) | 2 L + J + 1 ;   | ξ 0 | ( 1 2 ( L + J ) τ | ξ 0 | s s 0 100 ν ) | Ξ R ( s ; x , ξ ) | | ξ 0 | ( 1 + 2 ( L + J ) τ | ξ 0 | s s 0 100 ν ) ,   and therefore for s 0 s s 1   | ξ 0 | 2 | Ξ R ( s ; x , ξ ) | 7 4 | ξ 0 | .   Then we repeat the process and construct { s k }   , k = 1 , 2 , . . .   such that s k = s k 1 + 25 ν | ξ 0 | 2 J + L + k 1 .   Hence ( X R ( s ; x , ξ ) ; Ξ R ( s ; x , ξ ) )   is defined for 0 < s < s k   and | X R ( s k ; x , ξ ) | 2 L + J + k ;   | ξ k | ( 1 2 ( k + L + J ) ( 1 τ ) ) | Ξ R ( s ; x , ξ ) | | ξ k | ( 1 + 2 ( k + L + J ) ( 1 τ ) ) .   Then reasoning as in the case s < s 0   we get that 1 / 2 | ξ 0 | | Ξ R ( s ; x , ξ ) | 7 / 4 | ξ 0 | ,   for s 0 < s < s k .   This proves (i) and (ii). Finally (iv) follows from (i) and ( 4.23 ). The proof of Theorem  4.1.1 is complete.
Lemma 4.1.3 For any ψ   such that | ψ ( x ) | c τ x τ , τ > 1   there exists c ~   such that sup s , x , | ξ | = 1 0 s ψ ( X R ( r ; x , ξ ) ) d r < c ~ ,   where the constant c ~   depends just on τ   and c τ   .
Proof of Lemma  4.1.3 It follows from Theorem  4.1.1 .
Corollary 4.1.1 Let A 0 ( x )   be as in ( 4.1 )-( 4.4 ). Let B 1 ( x )   be an n × n   real matrix with entries in S ( R n )   , and define A 1 ( x ) = A 0 ( x ) + ε B 1 ( x )   , where ε   is chosen so that for all ξ R n   ( 2 ν ) 1 | A 1 ( x ) ξ | 2 ν | ξ |   .
Then, there exists ε 0 > 0 ,   such that A 1   verifies the basic assumption (i.e. A 1   is non-trapping).
Here ε 0   depends on M 1   , s M 1   , ν   , B 1   and on a finite number of seminorms of the difference ( 4.3 ).
In order to establish Corollary  4.1.1 , we need the following elementary o.d.e. lemma.
Lemma 4.1.4 Let y 0 ( s )   verify { d d s y 0 ( s ) = f 0 ( y 0 ( s ) ) for 0 < s < T y 0 ( 0 ) = z 0 ( s ) .   Let M ~ = sup 0 < s < T | y 0 ( s ) |   , and suppose that f 1   is given, with K = sup | y | M ~ | y f 1 ( y ) | ,   and let ε 0 > 0   be given. Then there exists ε ~ = ε ~ ( ε 0 , T , K ) ,   such that, if
sup 0 t T | f 0 ( y 0 ( t ) ) f 1 ( y 0 ( t ) ) | ε ~ , (4.34)
then there exists a unique solution y 1   to { d d s y 1 ( s ) = f 1 ( y 1 ( s ) ) for 0 < s < T y 1 ( 0 ) = z 0 ( s ) .   Moreover sup 0 s T | y 0 ( s ) y 1 ( s ) | ε 0 ,  
The proof is elementary and will be ommitted.
Proof of Corollary  4.1.1 Assume that ξ S n 1 ,   | x | 2 M 1   . By Theorem  4.1.1 , there exist s 2 M 1   and c 0 ,   such that | X 0 ( s 2 M 1 ; x , ξ | 4 M 1 ,   and such that c 0 1 | Ξ 0 ( s ; x , ξ ) | c 0   for all s   . Let now M ~ = sup | ( X 0 , Ξ 0 ) | ,   where the supremum is taken over all 0 s s 2 M 1 ,   | ξ | = 1 ,   and | x | 2 M 1 .   Let K   now be defined as in Lemma  4.1.4 for the system ( 4.5 ) associated to A 1 .   For ε 0   to be chosen, T = s 2 M 1 ,   let ε ~   be as in Lemma  4.1.4 . Choose now ε   so small that, in the terminology of Lemma  4.1.4 , for f 0 ,   f 1 ,   the system ( 4.5 ) associated to A 0 ,   A 1 ,   ( 4.34 ) holds. Then, if ε 0   in Lemma  4.1.4 is chosen small enough, we can conclude that there exists 0 < s ¯ < T   such that | X 1 ( s 2 M 1 ; x , ξ ) | 2 9 M 1 2 ;   ( 2 c 0 ) 1 | Ξ 1 ( s ¯ ; x , ξ ) | 2 c 0 ;   d d s | X 1 ( s ; x , ξ ) | 2 | s = s ¯ > 0 ;   and d 2 d s 2 | X 1 ( s ; x , ξ ) | 2 | s = s ¯ 2 ν 1 | A 1 ( X 1 ) Ξ 1 | ,   as long as | X 1 | M 1 .   Inserting this information in the proof of Theorem  4.1.1 we obtain the conclusions of Thoerem  4.1.1 for A 1 .  

2   Using the homogeneity property of the hamiltonian flow together with the methods given below in Theorem  4.1.1 and simple compactness and connectivity arguments, one can show that the seemingly weaker assumption that there exists s 0 ( δ , δ )   such that | X ( s 0 ; x 0 , ξ 0 ) | μ   implies our basic assumption.

4.2 The Continuous Dependence

In this subsection we shall deduce estimates concerning the continuous dependence of the flow associated to the truncated operator R ( x )   with respect to the initial value. These estimates will be given in Theorem  4.2.1 as a function of the parameter R   .
To simplify the notation we shall omit the sub-indexes denoting component or coordinates. So instead of
{ d d s X j R ( s ; x , ξ ) = 2 k = 1 n a j k R ( X R ( s ; x , ξ ) ) Ξ k R ( s ; x , ξ ) , d d s Ξ j R ( s ; x , ξ ) = k , l = 1 n x j a k l R ( X R ( s ; x , ξ ) ) Ξ k R ( s ; x , ξ ) Ξ l R ( s ; x , ξ ) , (4.35)
with a j k R ( )   defined in ( 4.12 )-( 4.13 ), we shall write
{ d d s X R ( s ; x , ξ ) = 2 a R ( X R ( s ; x , ξ ) ) Ξ R ( s ; x , ξ ) , d d s Ξ R ( s ; x , ξ ) = a R ( X R ( s ; x , ξ ) ) Ξ R ( s ; x , ξ ) Ξ R ( s ; x , ξ ) . (4.36)
Also, we shall omit the variables.
Hence ξ X R ( ) = ξ X R ( s ; x , ξ )   , ξ Ξ R ( ) = ξ Ξ R ( s ; x , ξ )   satisfy the system
{ d d s ( ξ X R ( ) ) = 2 a R ( X R ( ) ) ( ξ X R ( ) ) Ξ R ( ) + 2 a R ( X R ( ) ) ( ξ Ξ R ( ) ) , d d s ( ξ Ξ R ( ) ) = 2 a R ( X R ( ) ) ( ξ X R ( ) ) Ξ R ( ) Ξ R ( ) 2 a R ( X R ( ) ) Ξ R ( ) ( ξ Ξ R ( ) ) . (4.37)
We observe that x X R ( ) = x X R ( s ; x , ξ )   , x Ξ R ( ) = x Ξ R ( s ; x , ξ )   satisfy the system obtained from ( 4.37 ) by substituting ξ   by x   everywhere.
As a final simplification c   will denote in what follows a generic constant independent of R   which can change from line to line.
By homogeneity of the symbol of R ( x )   (see ( 2.16 ) in Section  2 ) for any R > 0   ,
X R ( s ; , x , t ξ ) = X R ( t s ; x , ξ ) , Ξ R ( s ; x , t ξ ) = t Ξ R ( t s : x , ξ ) , (4.38)
and consequently
( x α ξ β X R ) ( s ; x , t ξ ) = t | β | ( x α ξ β X R ) ( t s ; x , ξ ) , ( x α ξ β Ξ R ) ( s ; x , t ξ ) = t 1 | β | ( x α ξ β Ξ R ) ( t s ; x , ξ ) . (4.39)
So we can take | ξ | = 1   and consider s 0   . Therefore combining ( 4.18 ) of Theorem  4.1.1 and ( 4.37 ) it can be deduced that { d d s | ξ X R ( ) | c | a R ( X R ( ) ) | | ξ X R ( ) | + 2 | a R ( X R ( ) ) | | ξ Ξ R ( ) | , d d s | ξ Ξ R ( ) | c | 2 a R ( X R ( ) ) | | ξ X R ( ) | + c | a R ( X R ( ) ) | | ξ Ξ R ( ) | .   From our hypothesis on the decay of a ( )   and Lemma  4.1.3 we can define
{ f ( s ) = | ξ X R ( s ; x , ξ ) | exp ( c 0 s | a R ( X R ( s ; x , ξ ) ) | d s ) | ξ X R ( s ; x , ξ ) | , g ( s ) = | ξ Ξ R ( s ; x , ξ ) | exp ( c 0 s | a R ( X R ( s ; x , ξ ) ) | d s ) | ξ Ξ R ( s ; x , ξ ) | . (4.40)
From ( 4.39 ) it follows that
{ f ( s ) c | a R ( X R ( s ; x , ξ ) ) | g ( s ) c g ( s ) , g ( s ) c | 2 a R ( X R ( s ; x , ξ ) ) | f ( s ) c | 2 a R ( X R ( s ; ) ) | f ( s ) . (4.41)
We observe that if | x | > R   then ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) ) = ( x + 2 s ξ ~ , ξ )   as long as | x + 2 s ξ ~ | > R   , where A h ξ = ξ ~   .
Case 1Assume | x | > R   with | ξ | = 1   .
We shall assume that there exists a first s 1 > 0   such that | X R ( s 1 ; x , ξ ) | = R   , otherwise ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) ) = ( x + 2 s ξ ~ , ξ )   , for all s > 0   . One has that s 1 | x | + R   .
Since for s ( 0 , s 1 ) ,   ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) ) = ( x + 2 s ξ ~ , ξ ) ,   then | ξ X R ( s ; x , ξ ) | 2 s   and | ξ Ξ R ( s ; x , ξ ) | 1 .   The case s s 1   reduces to our next step.
Case 2 ( s < s 0 )If | x | M 1   , M 1   as in Theorem  4.1.1 , the analysis is given in case 3. Assume M 1 | x | R ,   1 / 2 < | ξ | 3 / 2   and
d d s | X R ( s ; x , ξ ) | | s = 0 0 . (4.42)
We also assume that ( 4.42 ) holds for s ( 0 , s 0 )   with s 0   defined in Theorem  4.1.1 and therefore comparable in size to R   because of (v) of that theorem. Otherwise we would have reached the outgoing situation (case 4) in an intermediate step.
Consider the majorized system ( 4.41 ) with data ( f ( 0 ) , g ( 0 ) ) = ( a , b ) , a , b > 0   .
Integrating we have g ( s ) b + c 0 s | 2 a ( X ( θ ; x , ξ ) ) | f ( θ ) d θ ,   so
f ( s ) a + b s + c 0 s ( 0 l | 2 a ( X ( θ ; x , ξ ) ) | f ( θ ) d θ ) d l a + b s + c 0 s ( s θ ) | 2 a ( X ( θ ; x , ξ ) ) | f ( θ ) d θ . (4.43)
Let s 2 = min { s ( 0 , c 0 R ) : f ( s ) = 2 ( a + b s ) }   .
Claim | X R ( s 2 ; x , ξ ) | R / 2   . Otherwise, by ( 4.43 ) one should have a + b s 2 c 0 s 2 ( s 2 θ ) | 2 a ( X ( θ ; x , ξ ) ) | ( a + b θ ) d θ ,   so c s 2 2 2 τ + 2 R τ + 2 1 2 ,   which is a contradiction for R R 0   with R 0   sufficiently large.
We repeat the argument assuming that ( 4.42 ) holds in the interval ( 0 , s j )   defining s j + 1 = min { s ( s j , c 0 R ) : f ( s ) = 2 j ( a + b s ) } .   Claim | X R ( s j ; x , ξ ) | R / 2 j   . Otherwise, we would have 2 c ( s j + 1 s j ) 2 2 ( τ + 2 ) j R ( τ + 2 ) 1 ,   and consequently, from ( 4.22 ) in Theorem  4.1.1 , c R 2 j s j + 1 s j R τ + 2 2 c 2 ( τ + 2 ) j ,   which is a contradiction if R 10 c τ 2 j   .So we can repeat the argument k   -times until R 2 k   , with f ( s k ) = 2 k ( a + b s k ) R ( a + b R )   since by Theorem  4.1.1 (v) we have that s k R .   Similarly, one gets that g ( s k ) R ( a + b R )   and | X R ( s k ; x , ξ ) | M ~   , with M ~   independent of R   .
Restarting the variable s   we are led to the following case.
Case 3Assume | X R ( 0 ; x , ξ ) | = | X ( 0 ; x , ξ ) | M ~   .
We consider the majorized system ( 4.41 ) with data ( a 1 , b 1 )   . From Theorem  4.1.1 there exists s * > 0   (independently of R   ) such that | X R ( s * ; x , ξ ) | = M ~ + 1 , with d d s | X R ( s ; x , ξ ) | | s = s * > 0 .   Integrating the system ( 4.41 ) we find that f ( s ) a 1 e c s , g ( s ) b 1 e c s for any s [ 0 , s * ) .   After restarting s   we are reduced to the following case.
Case 4(outgoing, i.e. ( 4.20 ) holds and d d s | X R ( s ; x , ξ ) | > 0   ).
Assume | X R ( 0 ; x , ξ ) | = | x | M ~ + 1 , with d d s | X R ( s ; x , ξ ) | | s = 0 > 0 .   We consider the majorized system ( 4.41 ) with data ( a 2 , b 2 )   .
Define h ( s ) = f ( s ) ( 1 + s 2 ) 2 + g ( s ) 1 + s 2 .   Thus,
h ( s ) d 0 1 + s 2 g ( s ) 1 + s 2 + d 1 1 + s 2 ( 1 + s 2 ) 2 | 2 a ( X R ( s ; x , ξ ) ) | f ( s ) ( 1 + s 2 ) 2
c 1 + s 2 h ( s ) .
Hence, from Theorem  4.1.1 one has
| ξ X R ( s ; x , ξ ) | f ( s ) c ( 1 + s 2 ) 2 ( a 2 + b 2 ) c ( 1 + R 2 ) 2 ( a 2 + b 2 ) ,
| ξ Ξ R ( s ; x , ξ ) | g ( s ) c ( 1 + s 2 ) 2 ( a 2 + b 2 ) c ( 1 + R 2 ) 2 ( a 2 + b 2 ) ,
as far as | X R ( s ; x , ξ ) | R   .
Case 5Assume | X R ( 0 ; x , ξ ) | = | x | R   with d d s | X R ( s ; x , ξ ) | | s = 0 > 0   .
Also we have the initial values ( | ξ X R ( 0 ; x , ξ ) | , | ξ Ξ R ( 0 ; x , ξ ) | ) = ( a 3 , b 3 )   .
So we have the solution | ξ X R ( s ; x , ξ ) | a 3 + c s , | ξ Ξ R ( s ; x , ξ ) | b 3 .   Collecting the information above we get
| ξ X R ( s ; x , ξ ) | , | ξ Ξ R ( s ; x , ξ ) | c ( 1 + | s | ) R 6 , (4.44)
for R R 0   .
To estimate ( x X R ( s ; x , ξ ) , x Ξ R ( s ; x , ξ ) )   as above we observe that it satisfies the system obtained from ( 4.37 ) by substituting x   by ξ   everywhere. As in the previous case, we majorize ( | x X R ( s ; x , ξ ) | , | x Ξ R ( s ; x , ξ ) | )   by the system in ( 4.41 ). So the same argument shows that | x X R ( s ; x , ξ ) | , | x Ξ R ( s ; x , ξ ) | c ( 1 + | s | ) R 6 ,   for R R 0   .
To estimate the higher order derivatives we first observe that if ( h 1 R α β ( s ) , h 2 R α β ( s ) ) = ( x α ξ β X R ( s ; x , ξ ) , x α ξ β Ξ R ( s ; x , ξ ) ) ,   then h 1 R α β ( s ) , h 2 R α β ( s ) )   satisfies the system { d d s h 1 R α β ( s ) = 2 a R ( X R ( ) ) ( h 1 R α β ( s ) ) Ξ R ( ) + 2 a R ( X R ( ) ) ( h 2 R α β ( s ) ) + Q 1 R α β , d d s ( h 2 R α β ( s ) ) = 2 a R ( X R ( ) ) ( h 1 R α β ( s ) ) Ξ R ( ) Ξ R ( ) 2 a R ( X R ( ) ) Ξ R ( ) ( h 2 R α β ( s ) ) + Q 2 R α β ,   where
Q j R α β = Q j R α β ( ( γ a ) 1 | γ | | α | + | β | + 1 ; ( x ν ξ μ X R ) | ν | + | μ | < | α | + | β | ; ( x ν ξ μ Ξ R ) | ν | + | μ | < | α | + | β | ) ,
j = 1 , 2   is a polynomial in its variables. In others words, ( h 1 R α β ( s ) , h 2 R α β ( s ) )   satisfies a system similar to that in ( 4.37 ) with external forces depending on the previous steps. We observe that each term in Q j R α β   , j = 1 , 2   has a factor of the form γ a ( X R ( s ; x , ξ ) )   with | γ | 1   . This guarantees that γ a ( X R ( s ; x , ξ ) ) = 0   if | X R ( s ; x , ξ ) | 2 R   and if | X R ( s ; x , ξ ) | 2 R   then | s | c ( | x | + R )   , (see Theorem  4.1.1 ). So when integrating the Q j R α β   's one can substitute in ( 4.44 ) the factor ( 1 + | s | )   by c ( | x | + R )   to get a bound independent of s   until the trajecory gets the free regime (i.e.
| X R ( s ; x , ξ | > R   ) which provides a linear in s   global bound. Using a recursive argument we can obtain the following estimates which, although no sharp, suffices for our purpose here.
Theorem 4.2.1 For any ( x , ξ ) R n × S n 1   and for any R R 0   with R 0   sufficiently large, the derivatives of the bicharacteristic flow satisfy : given α , β Z n   with | α | + | β | 1   there exists μ α , β Z +   with μ α , β | α | + | β |   such that
| x α ξ β X R ( s ; x , ξ ) | c ( | s | + ( | x | + R ) μ α , β ) ,
| x α ξ β Ξ R ( s ; x , ξ ) | c ( | s | + ( | x | + R ) μ α , β ) ,
for any s R   .
Combining the results in Theorem  4.2.1 with the identities in ( 4.38 )-( 4.39 ) we get the result for ( x , ξ ) R n × R n { 0 }   .

5 LINEAR ULTRAHYPERBOLIC EQUATIONS

In this section we shall deduce the linear estimates to be used in Section  6 in the proof of our nonlinear results. The first result in subsection  5.1 is the ultrahyperbolic version of Doi`s lemma [8, (see Lemma  2.2.2 in Section  2 ). This will allow us to establish the smoothing effects for linear ultrahyperbolic equations with variable second order coefficients.
Subsection  5.2 is concerned with the L 2   (and H s   ) well posedness of the associated linear problem ( 1.4 ) in the Introduction. To establish this result we follow an indirect approach. As we did in Section  4 -see ( 4.12 ) and ( 4.13 )-, we consider the truncation at infinity R ( x )   of the operator ( x )   , R ( x ) = θ ( x / R ) ( x ) + ( 1 θ ( x / R ) ) 0 .   For R   large enough we consider the bicharacteristic flow ( X R ( s ; x , ξ ) , Ξ R ( s ; x , ξ ) )   (studied in Section  4 ) associated to the operator R ( x )   and the corresponding integrating factor K R   . To obtain the L 2   local well posedness of the linear problem we combine several estimates for the operator R ( x )   and its associated “errors”, as function of R   , with the local smoothing effect obtained in subsection  5.1 .

5.1 Linear Ultrahyperbolic Smoothing

We shall begin this subsection by proving the ultrahyperbolic version of Doi's lemma [8, (see Lemma  2.2.2 in Section  2 ).
Lemma 5.1.1 Assume that the bicharacteristic flow is non-trapped -see basic assumption in subsection  4.1 of Section 4, and that a j k ( x ) = o ( | x | 1 )   as | x |   for all j , k = 1 , . . , n   . Suppose that λ L 1 ( [ 0 , ) ) C ( [ 0 , ) is positive and nonincreasing .   Then there exist p S 1 , 0 0   and c > 0   such that ( H h 2 p ) ( x , ξ ) λ ( | x | ) | ξ | c x , ξ R n .  
Proof of Lemma  5.1.1 Let M > 0   be a constant to be chosen. Let ψ C ( R )   with ψ ( t ) = 0   for t M 2   , ψ ( t ) = 1   for t ( M + 1 ) 2   and ψ ( t ) 0   for t R   . Let p 1 ( x , ξ ) = ξ 1 ψ ( | x | 2 ) H h 2 ( | x | 2 ) = 4 ξ 1 ψ ( | x | 2 ) ( x A h ( x ) ξ ) .   By straightforward calculation -see Proposition  4.1.1 of Section  4 
( H h 2 p 1 ) ( x , ξ ) = ξ 1 ψ ( | x | 2 ) ( H h 2 ( | x | 2 ) ) 2
+ ψ ( | x | 2 ) ξ 1 [ 8 | A ( x ) ξ | 2 + 8 j , k , l , m x l x j a l m ( x ) a j k ( x ) ξ k ξ m
4 j , k , l , m , p x j a j k ( x ) x k a l m ( x ) ξ l ξ m ] .
With the assumptions on A ( x )   and c 1 = 4 ν 2   one can fix M   sufficiently large such that ( H p 1 ) ( x , ξ ) c 1 ψ ( | x | 2 ) ξ 1 | ξ | 2 , x , ξ R n .   Now choose φ 1 C 0 ( R n )   with φ 1 ( x ) = 1   for | x | M + 1   . For ξ 0   , let p 2 ( x , ξ ) = 0 φ 1 ( X ( s ; x , ξ ) ) Ξ ( s ; x , ξ ) d s .   By theorem  4.1.1 in Section  4 , for each ( x 0 , ξ 0 ) R n × R n { 0 }   there is a neighborhood U   of ( x 0 , ξ 0 )   such that the integral defining p 2 ( x , ξ )   is over a fixed compact interval for all ( x , ξ ) U   so p 2 ( , )   is smooth. Furthermore, by homogeneity of the bicharacteristic flow (see ( 2.16 ) in Section  2 ) and a change of variable p 2 ( x , ξ ) = | ξ | 1 0 φ 1 ( X ( s ; x , ξ | ξ | ) ) | ξ | Ξ ( s ; x , ξ | ξ | ) d s .   Choose φ 2 C ( R n )   with φ 2 ( ξ ) = 0   for | ξ | 1   and φ 2 ( ξ ) = 1   for | ξ | 2   . Let p 3 ( x , ξ ) = φ 1 ( x ) φ 2 ( ξ ) p 2 ( x , ξ ) , x , ξ R n .   Then p 3 S 1 , 0 0   by the support properties of φ 1   and φ 2   , and
( H h 2 p 3 ) ( x , ξ ) = [ 2 j k a j k ( x ) ξ k x j φ 1 ( x ) ] φ 2 ( ξ ) p 2 ( x , ξ )
+ φ 1 ( x ) H h 2 φ 2 ( ξ ) p 2 ( x , ξ ) + ( φ 1 ( x ) ) 2 φ 2 ( ξ ) ξ .
Now let p 4 ( x , ξ ) = c 2 p 1 ( x , ξ ) + p 3 ( x , ξ ) ,   with c 2 > 0   sufficiently large, p 4   then satisfies | x α ξ β p 4 ( x , ξ ) | c α β x ξ | β | , α , β N n ,   and ( H h 2 p 4 ) ( x , ξ ) c 3 | ξ | c 4 , x , ξ R n ,   where c α β , c 3 , c 4 > 0   are constants.
To complete the construction of p ( x , ξ )   we observe that the proof of Lemma 2.3 in [8applies verbatim with q ( x , ξ ) = p 4 ( x , ξ )   , since the proof does not depend on the ellipticity of A ( x )   assumed in [8.
With this, the proof of Lemma  5.1.1 is completed.
Consider now systems of the form
{ t w = i H w + B w + C w + f , ( x , t ) R n × ( 0 , T ) , w ( x , 0 ) = w 0 ( x ) . (5.1)
Here w   and f   are C 2   -value functions on R n × ( 0 , T )   , H = ( 0 0 ) ,   where = j , k x j ( a j k ( x ) x k )   and A ( x ) = ( a j k ( x ) )   satisfies the assumptions in subsection  4.1 of Section  4 , B = ( Ψ b 1 b 2 ( x ) b ¯ 2 ( x ) Ψ b ¯ 1 )   where b 1 S 1 , 0 1   is odd in ξ   and b 2 C b ( R n : R n )   , and C = ( Ψ c 11 Ψ c 12 Ψ c 21 Ψ c 22 )   where c l m S 1 , 0 1   for l , m = 1 , 2   .
Next we consider the ultrahyperbolic linear Schrödinger (scalar) equation
{ t u = i u + Ψ b 1 u + b 2 ( x ) u ¯ + Ψ c 1 u + Ψ c 2 u ¯ + f , u ( x , 0 ) = u 0 ( x ) . (5.2)
Taking f = ( f , f ¯ ) T   , w = ( u , u ¯ ) T   and suitable c   's, equation ( 5.2 ) is reduced to a system as ( 5.1 ).
Theorem 5.1.1 Let s R   . Then there exists N = N ( n ) N   such that if | x α ξ β b 1 ( x , ξ ) | c α β x N ξ 1 | β | , x , ξ R n ,   and | x α b 2 ( x ) | c α β x N , x R n ,   then there exists T > 0   so that thefollowing holds: Let u 0 H s   .
(A) If f L 1 ( [ 0 , T ] : H s )   , then there is a unique solution u C ( [ 0 , T ] : H s )   of the IVP ( 5.2 ) satisfying
sup 0 t T u ( t ) H s + ( 0 T R n | J s + 1 / 2 u ( x , t ) | 2 x N d x d t ) 1 / 2
c u 0 H s + c 0 T f ( t ) H s d t .
(B) If f L 2 ( [ 0 , T ] : H s )   , then there is a unique solution u C ( [ 0 , T ] : H s )   of the IVP ( 5.2 ) satisfying
sup 0 t T u ( t ) H s 2 + 0 T R n | J s + 1 / 2 u ( x , t ) | 2 x N d x d t
c u 0 H s 2 + c T 0 T f ( t ) H s 2 d t .
(C) If J s 1 / 2 f L 2 ( R n × [ 0 , T ] : x N d x d t )   then there is a unique solution u C ( [ 0 , T ] : H s )   of the IVP ( 5.2 ) satisfying
sup 0 t T u ( t ) H s 2 + 0 T R n | J s + 1 / 2 u ( x , t ) | 2 x N d x d t
c u 0 H s 2 + c T 0 T R n | J s 1 / 2 f ( x , t ) | 2 x N d x d t .
Here c   depends on s ,   ν ,   A ( x ) ,   b 1 ,   b 2 ,   and c l m .  
The goal is to prove Theorem  5.1.1 . The a priori estimates needed to prove Theorem  5.1.1 are essentially reduced to the case s = 0   by the following commutator result.
Lemma 5.1.2 J s ( i H + B + C ) = ( i H + B ~ + C ~ ) J s ,   where
B ~ = i s j , k , l x j a k l ( x ) x j x k x l 3 J 2 ( 1 0 0 1 ) + B , (5.3)
C ~ = ( Ψ c ~ 11 Ψ c ~ 11 Ψ c ~ 11 Ψ c ~ 11 )   and c ~ l m S 1 , 0 0   for l , m = 1 , 2   .
A similar result applies to the scalar equation ( 5.2 ).
Proof of Lemma  5.1.2 The lemma follows from the classical S 1 , 0 0   pseudo-differential calculus of subsection  2.1 in Section  2 .
The next step is to obtain an a priori estimate for Theorem  5.1.1 in the case s = 0   . The last step will be the reduction to the case s = 0   (via Lemma  5.1.2 ). From this point, c   will denote a constant (not necessarily the same at each appearance) depending on s , ν , A ( x ) , b 1 , b 2 , c l m   and the hypothesis of Theorem  5.1.1 is assumed to be satisfied. In order to simplify notation, ( H s ) 2 = H s × H s   will be denoted merely by H s   .
Lemma 5.1.3 Let T > 0   . For all w C ( [ 0 , T ] : H 2 ) C 1 ( [ 0 , T ] : L 2 )   0 T R n | J 1 / 2 w ( x , t ) | 2 x N d x d t c ( 1 + T ) sup 0 t T w ( t ) L 2 2 + 0 T f ( t ) L 2 2 d t   and
0 T R n | J 1 / 2 w ( x , t ) | 2 x N d x d t
c ( ( 1 + T ) sup 0 t T w ( t ) L 2 2 + 0 T R n | J 1 / 2 f | 2 x N d x d t ) ,
where f ( x , t ) = t w ( i H + B ~ + C ~ ) w   and B ~ , C ~   were given in Lemma  5.1.2 .
Proof of Lemma  5.1.3 By Lemma  5.1.1 there exists a real-valued p S 1 , 0 0   and c > 0   such that ( H h 2 p ) ( x , ξ ) c x N | ξ | c , x , ξ R n ,   where c = c ( s )   is to be determined. Let k ( x , ξ ) = ( exp ( p ( x , ξ ) ) 0 0 exp ( p ( x , ξ ) ) )   and K = Ψ k   . Then K   is a diagonal 2 × 2   matrix of S 1 , 0 0   pseudo-differential operators. One can now calculate as follows
t K w , w L 2 × L 2 = K t w , w L 2 × L 2 + K w , t w L 2 × L 2
= ( i [ K H H K ] + K B ~ + B ~ * K ) w , w L 2 × L 2
+ [ K C ~ + C ~ * K ] w , w L 2 × L 2 + ( K f , w L 2 × L 2 + K w , f L 2 × L 2 )
= I + I I + I I I .
Disregarding symbols of order 0   , the first order symbol of i [ K H H K ] + K B ~ + B ~ * K   is
e p ( ξ h 2 x p x h 2 ξ p 0 0 ξ h 2 x p x h 2 ξ p )
2 s e p ( j k l x j a k l ( x ) ξ j ξ k ξ l ξ 2 0 0 j k l x j a k l ( x ) ξ j ξ k ξ l ξ 2 )
+ e p ( 2 Re b 1 2 i b 2 ( x ) ξ 2 i b ¯ 2 ( x ) ξ 2 Re b 1 ) .
Choosing c ( s )   large enough and using the matrix version of the sharp Gårding inequality in [12, it follows that
| I | c ( λ J 1 0 0 λ J 1 ) w , w + c w L 2 2
c R n | J 1 / 2 w ( x , t ) | 2 x N d x + c w L 2 2 .
Using that both K   and C   are of order 0   | I I | c w L 2 2 .   The estimate of I I I   is split into two cases according to the desired norm of f   . This is identical to the elliptic case given in Section  2 .
Combining estimates, integrating in t   and using the L 2   -boundedness of K   , Lemma  5.1.3 follows.
Remark 5.1.1 We will use Lemma  5.1.3 for ω = ( J s u , J s u ¯ )   , where u   solves the scalar equation ( 5.2 ).

5.2 Linear Ultrahyperbolic L 2   -Well posedness

Our goal in this subsection is to established the following result,
Lemma 5.2.1 Let T > 0   . For all u C ( [ 0 , T ] ; H 2 ) C 1 ( [ 0 , T ] ; H 2 )   the following two estimates hold.
  • (A) sup 0 t T u L 2 c u ( 0 ) L 2 + c T sup 0 t T u L 2 + c 0 T f L 2 d t   ,
  • (B) sup 0 t T u L 2 2 c u ( 0 ) L 2 2 + c T sup 0 t T u L 2 2 + c 0 T f L 2 2 d t .  
Here f = t u ( i u + Ψ b 1 u + b 2 u ¯ + Ψ c 1 u + Ψ c 2 u ¯ )   analogously to Lemma  5.1.3 .
The proof of Lemma  5.2.1 involves several steps. The first one is to cut 0   at infinity. Therefore for R   large enough and to be fixed later on we define a j k R ( x ) = θ ( x R ) a j k ( x ) + ( 1 θ ( x R ) ) a j k 0 = a j k 0 + θ ( x R ) ( a j k ( x ) a j k 0 ) ,   with θ   a smooth cut off function such that θ ( x ) = 1   if | x | < 1   and θ ( x ) = 0   if | x | > 2   . Define R ( x ) = j , k x j ( a j k R ( x ) x k ) ,   h 2 R   its corresponding hamiltonian and R ( x ) = ( x ) R ( x ) .   Therefore R ( x ) = j , k x j ( e j k ( x ) x k )   with e j k ( x ) = ( 1 θ ( x R ) ) ( a j k ( x ) a j k 0 ) .   The main point is to apply a suitable pseudo-differential operator K R   to the corresponding system in order to cancel the first order terms. We begin with the definition and some properties of K R   and its symbol. We recall that χ C ( R )   , χ ( t ) = 0   for | t | 1   and χ ( t ) = 1   for | t | 2   .
Definition 5.2.1 Let R = 2 j 0   for j 0 N .   Recall the definition of B ~   in Lemma  5.1.2 . We define
  • (i) b R ( x , ξ ) = s j , k , l x j a j k R ( x ) ξ j ξ k ξ l ξ 2 Re b 1 ( x , ξ ) ,  
  • (ii) p R ( x , ξ ) = χ ( 1 2 | ξ | ) 0 b R ( X R ( σ ; x , ξ ) , Ξ ( σ ; x , ξ ) ) d σ ,  
  • (iii) p e R ( x , ξ ) = 1 2 ( p R ( x , ξ ) + p R ( x , ξ ) ) ,  
  • (iv) k R ( x , ξ ) = exp ( p e R ( x , ξ ) ) ,  
  • (v) k ~ R ( x , ξ ) = exp ( p e R ( x , ξ ) ) ,  
  • (vi) B R = Ψ b R , P R = Ψ p R , P e R = Ψ p e R , K R = Ψ k R , K ~ R = Ψ k ~ R ,  
  • (vii) b R = b 0 R + j = j 0 b j R   , with b j R ( x ) = ( θ ( x 2 j + 1 ) θ ( x 2 j ) ) b R ( x )   , and analogously p j R   , p e j R   , k j R   , and k ~ j R .  
Some comments about the above definition are in order.
We recall that applying the operator J s   to the equation in ( 5.2 ) the symbol of the first order term for J s u   is given by (see ( 5.3 )) s j , k , l x j a j k ( x ) ξ j ξ k ξ l ξ 2 Re b 1 ( x , ξ ) .   So b R ( x , ξ )   in (i) is the approximation (due to our truncation of the second order term) of this symbol.
The symbol of the integrating factor needed to “cancel” Ψ b R   is given in (iv). The reason to use the even function p e R ( x , ξ )   in (iii) instead of that in (ii) is to preserve the symmetry needed in the integration by parts use to handle the first order term in u ¯   .
Finally, K ~ R ( x , ξ )   is the symbol of a cuasi-inverse of Ψ K R   .
Lemma 5.2.2
  • (i) p e R , k R , k ~ R   are all even in ξ   , and real.
  • (ii) ( ξ h 2 R x p R x h 2 R ξ p R ) ( x , ξ ) = b R ( x , ξ ) + r 1 ( x , ξ )   where r 1 S 1 , 0 .  
  • (iii) ( ξ h 2 R x p e R x h 2 R ξ p e R ) ( x , ξ ) = b R ( x , ξ ) + r 2 ( x , ξ )   where r 2 S 1 , 0 .  
  • (iv) ( ξ h 2 R x k R x h 2 R ξ k R ) ( x , ξ ) = k R ( x , ξ ) b R ( x , ξ ) + r 3 ( x , ξ )   where r 3 S 1 , 0 .  
  • (v) ( ξ h 2 R x k ~ R x h 2 R ξ k ~ R ) ( x , ξ ) = k ~ R ( x , ξ ) b R ( x , ξ ) + r 4 ( x , ξ )   where r 4 S 1 , 0 .  
  • (vi) Let ψ C ( R )   , ψ ( t ) = 1   if t 1 2   and ψ ( t ) = 1   if t 1 2 , 1 ψ 1   . Define
    a j R ( z ; x , ξ ) = χ ( 1 102 j | x | ) p e j R ( z + 2 j ψ ( x A h ξ | x | | ξ | ) A h ξ | ξ | , ξ ) θ ( z 2 j + 1 ) , (5.4)
    j = 0 , 1 , . . .   . Then a j R C 0 ( B 2 j + 1 ( 0 ) ; S 1 , 0 0 )   and
    p e j R ( x , ξ ) = p e j R ( x , ξ ) ( 1 χ ( 1 102 j | x | ) ) + a j R ( P ( x , A h ξ ) ; x , ξ ) . (5.5)
  • (vii) Set q ( x , ξ ) = j p e j R ( 1 χ ( 1 102 j | x | ) )   , and a R ( z ; x , ξ ) = exp ( q ) ( exp ( j a j R ( z ; x , ξ ) ) 1 ) .   Then q S 1 , 0 0   , q ( , ξ ) S ( R n )   uniformly in ξ   , a R S ( R n ; S 1 , 0 0 )   , and
    k R ( x , ξ ) = a R ( P ( x , A h ξ ) ; x , ξ ) + exp ( q ( x , ξ ) ) . (5.6)
  • (viii) The seminorms of the remainders r k   , k = 1 , 2 , 3 , 4   , and of a j R   and q   grow as R N 0   with N 0 Z +   depending just on the dimension.
Parts (i)–(iii) are preliminary results needed in the proof of (iv). The crucial point (iv) shows at the level of the symbols that the commutator of K R   and R ( x )   “cancels” K R B R   . In (vi) and (vii) we prove that the symbol k R ( x , ξ )   is in the class introduced in Section  3 .
Proof of Lemma  5.2.2 For simplicity of the exposition we shall drop the index R   .
The proof of (i) is clear.
As for (ii), it is sufficient to work with each b j   and the corresponding p j   . We get
( ξ h 2 R x p j x h 2 R ξ p j ) ( x , ξ )
= 0 b j ( X ( σ ; x , ξ ) , Ξ ( σ ; x , ξ ) ) d σ ( x h 2 R ξ [ χ ( 1 2 | ξ | ) ] )
+ χ ( 1 2 | ξ | ) ( ξ h 2 R x x h 2 R x ) 0 b j ( X ( σ ; x , ξ ) , Ξ ( σ ; x , ξ ) ) d σ .
For the first factor of the first term we use Theorem  4.1.1 in Section  4 to see that is bounded by c N 2 N j   , and therefore independently of R   . Then by Theorem  4.2.1 in Section  4 we prove that the derivatives x α ξ β   inside the integral are bounded by c α β ( | σ | + ( | x | + R ) μ α β )   for α , β N n   . If | σ | > s 0   where s 0   is given in Theorem  4.1.1 of Section  4 , then | σ | < c | X ( σ ; x , ξ ) | < c 2 j   . Otherwise | σ | < c | x | < c R   , where the last inequality follows from the support properties of x h 2 R   . Therefore from the decay of b j   and A   we conclude that belongs to S 1 , 0   and the corresponding seminorms grow like powers of R   . By Lemma  2.2.1 in Section  2 , the second term equals χ ( 1 2 | ξ | ) b j ( x , ξ ) = b j ( x , ξ ) + ( χ ( 1 2 | ξ | ) 1 ) b j ( x , ξ ) .   Notice that the last term above is compactly supported in ξ   .
(iii) h 2 R   is homogeneous of degree 2 in ξ   , so
( ξ h 2 R x x h 2 R ξ ) [ p ( x , ξ ) ]
= ( ξ h 2 R ) ( x , ξ ) ( x p ) ( x , ξ ) + ( x h 2 R ) ( x , ξ ) ( ξ p ) ( x , ξ )
= ( ξ h 2 R x p + x h 2 R ξ p ) ( x , ξ ) .
By (ii), this equals b ( x , ξ ) r 1 ( x , ξ ) = b ( x , ξ ) r 1 ( x , ξ ) ,   since b   is odd in ξ   . This proves (iii).
(iv) By the chain rule -see also the proof of (ii),
[ ξ h 2 R x k x h 2 R ξ k ] ( x , ξ )
= k ( x , ξ ) ( ξ h 2 R x p e x h 2 R ξ p e ) ( x , ξ )
= k ( x , ξ ) ( b ( x , ξ ) + r 2 ( x , ξ ) ) ,
where r 2 ( x , ξ )   is compactly supported in ξ   .
(v) is similar to (iv).
(vi) From the definition of a j   we have supp a j ( ; x , ξ ) B 2 j + 1 ( 0 )   .
By homogeneity, b j ( X ( σ ; x , ξ ) , Ξ ( σ ; x , ξ ) ) = b j ( X ( | ξ | σ ; x , ξ | ξ | ) , | ξ | Ξ ( | ξ | σ ; x , ξ | ξ | ) ) ,   so it follows by a change of variable that
p e j ( x , ξ ) = 1 2 χ ( 1 2 | ξ | ) 1 | ξ | ( 0 b j ( X ( σ ; x , ξ | ξ | ) , | ξ | Ξ ( σ ; x , ξ | ξ | ) ) d σ 0 b j ( X ( σ ; x , ξ | ξ | ) , | ξ | Ξ ( σ ; x , ξ | ξ | ) ) d σ ) . (5.7)
Thus combining ( 5.7 ) and Theorem  4.2.1 in Section  4 and proceeding as in part (ii) we obtain that a j ( z ; , ) S 1 , 0 0   uniformly in z   , and so are in the derivatives with respect to z   .
Let us prove ( 5.5 ). We have to see that
χ ( 1 102 j | x | ) p e j ( x , ξ ) = χ ( 1 102 j | x | ) p e j ( P ( x , A h ξ ) + 2 j ψ ( x A h ξ | x | | ξ | ) A h ξ | ξ | , ξ ) θ ( P ( x , A h ξ ) 2 j + 1 ) . (5.8)
Consider different cases.
1 o ̲   | x | 102 j   : Here χ ( | | x | 102 j | ) = 0   so RHS=LHS in ( 5.8 ).
2 o ̲   | ξ | 2   : R H S = L H S = 0 .   3 o ̲   | P ( x , A h ξ ) | 2 j   and | x | 102 j   : First note that P ( P ( x , A h ξ ) + 2 j ψ ( x A h ξ | x | | ξ | ) A h ξ | ξ | , A h ξ ) = P ( x , A h ξ ) .   Next, if | P ( y , A h ξ ) | 2 j R   for some y R n   , then b j ( X ( σ ; y , ξ ) , Ξ ( σ ; y , ξ ) ) = 0   for all σ R   because b j   has x   -support in B 2 j ( 0 ) ,   and the bicharacteristics are lines. It follows that RHS=LHS=0 in ( 5.8 ).
4 o ̲   | x | 102 j , | ξ | 2 , | P ( x , A h ξ ) | 2 j   : Here χ ( | ξ | / 2 ) = 1   and θ ( P ( x , A h ξ ) 2 j + 1 ) = 1 .   | x A h ξ | = | x | | ξ | 1 | P ( x , A h ξ ) | 2 | x | 2 99 100 | x | | ξ | .   Now split in two subcases according to the sign of x A h ξ   .
4 o ̲ ( a )   x A h ξ 99 100 | x | | ξ |   : Here ψ ( x A h ξ | x | | ξ | ) = 1   . Suppose σ 1 2 ( x A h ξ | ξ | 2 j + 1 )   . Then ( x 2 σ A h ξ | ξ | ) A h ξ | ξ | = x A h ξ | ξ | 2 σ 2 j + 1 .   Hence X ( σ ; x , ξ | ξ | ) = x 2 σ A h ξ | ξ |   and b j ( X ( σ ; x , ξ | ξ | ) , | ξ | Ξ ( σ ; x , ξ | ξ | ) ) = 0 .   By ( 5.7 ) and since 1 2 ( x A h ξ | ξ | 2 j + 1 ) > 0   ,
p e j ( x , ξ ) =
1 2 χ ( 1 2 | ξ | ) 1 | ξ | 1 / 2 ( x A h ξ / | ξ | 2 j + 1 ) b j ( X ( σ ; x , ξ | ξ | ) , | ξ | Ξ ( σ ; x , ξ | ξ | ) ) d σ .
Doing the change of variable τ = σ 1 2 ( x A h ξ | ξ | 2 j + 1 )   ,
p e j ( x , ξ ) =
1 2 χ ( 1 2 | ξ | ) 1 | ξ | 0 b j ( X ( τ + 1 2 ( x A h ξ | ξ | 2 j + 1 ) ; x , ξ | ξ | ) ,
, | ξ | Ξ ( τ + 1 2 ( x A h ξ | ξ | 2 j + 1 ) ; x , ξ | ξ | ) ) d τ .
Now,
X ( τ + 1 2 ( x A h ξ | ξ | 2 j + 1 ) ; x , ξ | ξ | )
= X ( τ ; X ( 1 2 ( x A h ξ | ξ | 2 j + 1 ) ; x , ξ | ξ | ) , Ξ ( 1 2 ( x A h ξ | ξ | 2 j + 1 ) ; x , ξ | ξ | ) )
= X ( τ ; P ( x , A h ξ ) + 2 j + 1 A h ξ | ξ | , ξ | ξ | ) ,
and similarly, Ξ ( τ + 1 2 ( x A h ξ | ξ | 2 j + 1 ) ; x , ξ | ξ | ) = Ξ ( τ ; P ( x , A h ξ ) + 2 j + 1 A h ξ | ξ | , ξ | ξ | )   If τ 0   then b j ( X ( τ ; P ( x , A h ξ ) + 2 j + 1 A h ξ | ξ | , ξ | ξ | ) , | ξ | Ξ ( τ ; P ( x , A h ξ ) + 2 j + 1 A h ξ | ξ | , ξ | ξ | ) ) = 0 .   Hence by ( 5.7 ) p e j ( x , ξ ) = p e j ( P ( x , A h ξ ) + 2 j + 1 ψ ( x A h ξ | x | | ξ | ) A h ξ | ξ | , ξ ) .   Therefore, ( 5.8 ) holds. This finishes case 4 o ̲ ( a )   .
4 o ̲ ( b )   x A h ξ 99 100 | x | | ξ |   : Here ψ ( x A h ξ | x | | ξ | ) = 1   .
Suppose σ 1 2 ( x A h ξ | ξ | + 2 j + 1 )   . Then ( x 2 σ A h ξ | ξ | ) A h ξ | ξ | = x A h ξ | ξ | 2 σ 2 j + 1 .   Hence X ( σ ; x , ξ | ξ | ) = x 2 σ A h ξ | ξ |   and b j ( X ( σ ; x , ξ | ξ | ) , | ξ | Ξ ( σ ; x , ξ | ξ | ) ) = 0 .   Therefore since 1 2 ( 2 j + 1 + x A h ξ | ξ | ) < 0   ,
p e j ( x , ξ )
= 1 2 χ ( 1 2 | ξ | ) 1 | ξ | 1 / 2 ( 2 j + 1 + x A h ξ / | ξ | ) b j ( X ( σ ; x , ξ | ξ | ) , | ξ | Ξ ( σ ; x , ξ | ξ | ) ) d σ .
Let τ = σ 1 2 ( 2 j + 1 + x A h ξ | ξ | )   . Then
p e j ( x , ξ ) =
1 2 χ ( 1 2 | ξ | ) 1 | ξ | 0 b j ( X ( τ + 1 2 ( 2 j + 1 + x A h ξ | ξ | ) ; x , ξ | ξ | ) ,
, | ξ | Ξ ( τ + 1 2 ( 2 j + 1 + x A h ξ | ξ | ) ; x , ξ | ξ | ) ) d τ .
But X ( τ + 1 2 ( 2 j + 1 + x A h ξ | ξ | ) ; x , ξ | ξ | ) = X ( τ ; P ( x , A h ξ ) 2 j + 1 A h ξ | ξ | , ξ | ξ | )   and Ξ ( τ + 1 2 ( 2 j + 1 + x A h ξ | ξ | ) ; x , ξ | ξ | ) = Ξ ( τ ; P ( x , A h ξ ) 2 j + 1 A h ξ | ξ | , ξ | ξ | ) .   Since | X ( τ ; P ( x , A h ξ ) 2 j + 1 A h ξ | ξ | , ξ | ξ | ) | 2 j + 1 for τ 0 ,   it follows from ( 5.7 ) that p e j ( x , ξ ) = p e j ( P ( x , A h ξ ) + 2 j + 1 ψ ( x A h ξ | x | | ξ | ) A h ξ | ξ | , ξ )   and ( 5.8 ) holds.
For the proof of (vii) notice first that p e j ( 1 χ ( 1 102 j | x | ) )   has compact support in x   and therefore belongs to S 1 , 0 0   . Thanks to the decay properties of b j   we can summ in j   to obtain the same property for q ( x , ξ )   . The other conclusions follow by inspection from (vi) and the definition of k   .
Finally notice that (viii) follows from the previous steps. This completes the proof of Lemma  5.2.2 .
Definition 5.2.2 The equivalence relation =   is given by A 1 = A 2 ( A 1 A 2 ) is L 2 -bounded with norm C ( R ) ,   where C ( R )   grows at most polynomially in R .  
Lemma 5.2.3 Let q S 1 , 0 1   with q ( , ξ ) S ( R n )   uniformly in ξ   . Let a S ( R n ; S 1 , 0 0 )   and d ( x , ξ ) = a ( P ( x , A h ξ ) ; x , ξ ) χ ( | ξ | ) .   Then
  • (i) Ψ q d = Ψ q Ψ d = Ψ d Ψ q ,  
  • (ii) Ψ q k R = Ψ q K R = K R Ψ q = Ψ q ( K R ) * = ( K R ) * Ψ q ,  
  • (iii) Ψ q k ~ R = Ψ q K ~ R = K ~ R Ψ q = Ψ q ( K ~ R ) * = ( K ~ R ) * Ψ q ,  
  • (iv) i [ ( K R ) * R R ( K R ) * ] = Ψ k R b R ,  
  • (v) i [ ( K ~ R ) * R R ( K ~ R ) * ] = Ψ k ~ R b R .  
The constants for the above inequalities are bounded by R N   for some fixed power N   , and by some fixed number of seminorms of q   and a .  
Proof of Lemma  5.2.3  The proof is based on the calculus developed in Section  3 .(i) Let φ ( x ) = ( 1 + | x | 2 ) N = x 2 N   , and 1 φ q = q ~   . Then for q = φ q ~   one has
Ψ q Ψ d = Ψ q ~ φ Ψ d = Ψ q ~ φ Ψ d = Ψ q ~ d = Ψ q d .
Here it was used that φ d   behaves as a classical symbol because for our purposes just finitely many derivatives of the symbol in ξ   are needed. This number of derivatives determines the choice of N   .
Similarly,
Ψ d Ψ q = Ψ d φ ( I Δ ) J 2 Ψ q ~
= Ψ d φ ξ 2 Ψ ξ ) 2 q ~
= Ψ d φ ξ 2 ξ 2 q ~ = Ψ q ~ d .
In this case we used that φ ( I Δ )   is a partial differential operator with decay in the coefficients and that φ ( x ) ξ 2   behaves as a symbol in S 1 , 0 2   for all N   .
(ii) Ψ q k R = Ψ q K R = K R Ψ q   follows from the decomposition in Lemma  5.2.2 (vii) of k R   into a sum of an S 1 , 0 0   symbol and a symbol of the type d   in (i). For the remainder of (ii),
Ψ q ( K R ) * = ( K R Ψ q * ) * = ( K R Ψ q ¯ ) *
= ( Ψ k R q ¯ ) * = Ψ k ¯ R q = Ψ k R q
because k R q ¯ S 1 , 0 1   and k R   is real–valued. Similarly,
( K R ) * Ψ q = ( Ψ q * K R ) * = ( Ψ q ¯ K R ) *
( Ψ q ¯ k R ) * = Ψ q k R .
(iii) is similar to (ii).
(iv) k R b R S 1 , 0 1   and is real–valued so, taking adjoints, it suffices to show that i [ K R R R K R ] = Ψ k R b R .   Now write R = [ R ( A h ) ] + ( A h )   to see that R   is a compactly supported perturbation of a constant coefficient operator. By Theorem  3.3.1 in Section  3 and the decomposition in Lemma  5.2.2 (vii), it follows that i [ K R R R K R ] = Ψ ξ h 2 R x k R x h 2 R ξ k R = Ψ b R k R   using Lemma  5.2.2 (iv) in the last equivalence.
(v) is similar to (iv) using Lemma  5.2.2 (v) instead of Lemma  5.2.2 (iv).
This completes the proof of Lemma  5.2.3 .
In order to prove Lemma  5.2.1 we still need some technical results. Recall that = R + R   and that K R   and b R   were given in Definition  5.2.1 .
Lemma 5.2.4 There exists N 0   large enough such that
  • (i) ( K R ) * u L 2 = O ( R N 0 u L 2 ) ,  
  • (ii) i [ R , ( K R ) * ] u + ( K R ) * Ψ b R u L 2 c R N 0 u L 2 ,  
  • (iii) ( K R ) * b 2 ( x ) u ¯ b 2 ( x ) ( K R ) * u ¯ L 2 c R N 0 u L 2 .  
  • (iv) ( K R ) * Ψ I m b 1 u Ψ I m b 1 ( K R ) * u L 2 c R N 0 u L 2 .  
Proof of Lemma  5.2.4 Part (i) follows from the descomposition of Lemma  5.2.2 (vii) and Theorem  3.2.1 in Section  3 . Part (ii) follows from Lemma  5.2.2 (iv), Theorems  3.3.1  3.3.3 in Section  3 and that k R   is real. As for (iii) follows from Theorems  3.3.1 - 3.3.3 in Section  3 and that k R   is even in ξ   . Finally (iv) follows from Lemma  5.2.3 (ii).
Lemma 5.2.5 Given M   there are R   large enough and N ( M )   such that | i [ R , ( K R ) * ] u , ( K R ) * u | R M J 1 / 2 ( u ) x N L 2 2 + O ( R N ( M ) = N 0 u L 2 2 ) ,   with N 0   as in Lemma  5.2.4 .
Remark 5.2.1 Here M   can be taken arbitrary large since the coefficients are in S   . However it suffices to assume ( 4.16 ) in Section  4 for some τ   sufficiently large.
Proof of Lemma  5.2.5 We have i [ R , ( K R ) * ] = i ( R ( K R ) * ( K R ) * R ) .   Take adjoints to get
i ( K R ( R ) * ( R ) * K R ) . (5.9)
Recall that R = ( R ) = j , k x j ( e j k x k ) ,   with e j k = ( 1 θ ( x R ) ) ( a j k ( x ) a j k 0 ) .   Thus we can see R   as a second order differential operator with coefficients of the form e j k ( x ) = 1 R M e ~ j k R ( x )   , and e ~ j k R   with decay uniform in R   . Then ( R ) * = e j k R ( x ) 2 x j x k + b 3 ,   for some b 3   with the right decay. Therefore this term gives bounds of the type O ( R N 0 u L 2 2 )   in ( 5.9 ), just using Lemma  5.2.2 (ii).
Now we use Theorem  3.3.1 in Section  3 and Lemma  5.2.2 (iv) to get K R ( R ) * ( R ) * K R = Ψ α R + zero order terms ,   with α R ( x , ξ ) = j , k ( x e j k ξ j ξ k ξ k R e j k R ξ ( ξ j ξ k ) x k R ) ,   and the zero order terms have bounds in L 2   which grow as R N 0 .   Notice that α R ( x , ξ ) S 1 , 0 1   , with uniform O ( R M )   decay in x   . We need to study Ψ α R * u , ( K R ) * u .   We observe that Ψ α R * = Ψ α R ¯ + zero order   . Therefore we will work with Ψ α R   , being analogous the calculations for Ψ α R ¯   . Define φ ( x ) = x N   . Then
R M Ψ α R u , ( K R ) * u = R M φ φ 2 Ψ α R J 1 / 2 J 1 J 1 / 2 u , φ ( K R ) * u = φ J 1 / 2 φ 2 Ψ α R J 1 J 1 / 2 u , φ ( K R ) * u = φ J 1 / 2 φ φ 3 Ψ α R J 1 J 1 / 2 u , φ ( K R ) * u = φ J 1 / 2 φ 3 Ψ α R J 1 φ J 1 / 2 u , φ ( K R ) * u = J 1 / 2 φ φ 3 Ψ α R J 1 φ J 1 / 2 u , φ ( K R ) * u = φ 3 Ψ α R J 1 φ J 1 / 2 u , φ J 1 / 2 φ ( K R ) * u φ 3 Ψ α R J 1 φ J 1 / 2 u L 2 φ J 1 / 2 φ ( K R ) * u L 2 C φ J 1 / 2 u L 2 φ J 1 / 2 φ ( K R ) * u L 2 (5.10)
as desired. Lemma  5.2.5 is proved.
Next, with K ~ R   as in Definition  5.2.1 , we define E R = I K ~ R ( K R ) * .  
Lemma 5.2.6 There exists N 0   such that
  • (i) E R u L 2 C R N 0 u L 2 ,  
  • (ii) Let q S 1 , 0 1   with q ( . , ξ ) S ( R n )   uniformly in ξ   . Then Ψ q E R u L 2 + E R Ψ q u L 2 C R N 0 u L 2 .   As a consequence E R Ψ b 1 u L 2 + E R b 2 . u ¯ L 2 C R N 0 u L 2 ,  
  • (iii) | [ R , E R ] u , E R u | C R N 0 u L 2 2 ,  
  • (iv) | i [ , E R ] u , E R u | C R N 0 u L 2 2 .  
Proof of Lemma  5.2.6 Part (i) follows from the descomposition of Lemma  5.2.2 (vii) and Theorem  3.2.1 in Section  3 .
As for (ii) Ψ q E R = Ψ q Ψ q K ~ R ( K R ) * = Ψ q Ψ q k ~ ( K R ) * = Ψ q Ψ q k ~ k = 0 ,   by using Lemma  5.2.2 (iii) and (ii). Similarly E R Ψ q = Ψ q K ~ R ( K R ) * Ψ q = Ψ q K ~ R Ψ q k = Ψ q Ψ q k ~ k = 0 .   For (iii) it is enough to prove that [ R , E R ] u   is L 2   -bounded. But
i [ R , E R ] = i ( R K ~ R ( K R ) * K ~ R ( K R ) * R )
= i ( ( R K ~ R K ~ R R ) ( K R ) * + i K ~ R ( R ( K R ) * ( K R ) * R ) )
= Ψ ( k ~ R b R ) ( K R ) * K ~ R Ψ ( k R b R )
= Ψ ( k ~ R b R k R ) Ψ ( k R b R k ~ R ) = 0 .
Finally let us prove (iv). We have = R + R ,   and E R = I K ~ R ( K R ) * .   Part (iii) gives that [ R , E R ] u , E R u   has the right bound. Hence we need to understand [ R , K ~ R ( K R ) * ] .   We have
R K ~ R ( K R ) * K ~ R ( K R ) * R =
( R K ~ R K ~ R R ) ( K R ) * + K ~ R ( R ( K R ) * ( K R ) * R ) .
For both ( R K ~ R K ~ R R )   and ( R ( K R ) * ( K R ) * R )   we can use Theorems  3.3.1  3.3.3 in Section  3 so that they can be written as Ψ β j + zero order terms   with Ψ β j , j = 1 , 2   classical first order pseudo-differential operators with the right decay in x   as we did in the proof of Lemma  5.2.5 .
But from Lemma  5.2.3  Ψ β 1 ( K R ) * = Ψ β 1 k R , K ~ R Ψ β 2 = Ψ β 2 k ~ R ,   and we can apply part (ii).
Proof of Lemma  5.2.1 Notice first that part A follows making f = 0   in part B and using Duhamel's principle. We study the problem { t u = i u + Ψ b 1 u + b 2 ( x ) u ¯ + Ψ c 1 u + Ψ c 2 u ¯ + f u ( x , 0 ) = u 0 .   We also have = R + R .   Applying the operator ( K R ) *   and observing that if s = 0   in Definition  5.2.1 then ψ b R = ψ Re b 1 )   , we get from Lemma  5.2.4 
t ( K R ) * u = i [ R , ( K R ) * ] u + ( K R ) * Ψ b 1 u + i R ( K R ) * u
+ i R ( K R ) * u + i [ R , ( K R ) * ] u + ( K R ) * b 2 ( x ) u ¯
+ ( K R ) * f + zero order terms
= i ( K R ) * u + b 2 ( x ) ( K R ) * u ¯ + i [ E R , ( K R ) * ] u
+ i Ψ I m b 1 ( K R ) * u + ( K R ) * f + zero order terms .
Define v R = ( K R ) * u   . Then we have
t v R , v R = i v R , v R + i [ R , ( K R ) * ] u , v R + b 2 ( x ) v R ¯ , v R + i Ψ I m b 1 v R , v R
+ O ( R N 0 u L 2 v R L 2 ) + O ( R N 0 u L 2 2 ) + O ( f L 2 2 )
= i v R , v R + i [ R , ( K R ) * ] u , v R + i Ψ I m b 1 v R , v R
+ O ( R N 0 u L 2 v R L 2 ) + O ( R N 0 u L 2 2 ) + O ( f L 2 2 ) ,
where the last step follows by integration by parts. Taking the real part of both sides and using Lemma  5.2.5 and Garding's inequality for Re i Ψ I m b 1 v R , v R   we get after integration in the temporal variable
sup 0 t T v R ( t ) L 2 2 C R N 0 u ( 0 ) L 2 2
+ R M 0 T J 1 / 2 ( u ) x N L 2 2 d t + C R N ( M ) + N 0 T sup 0 < t < T u ( t ) L 2 2 + C 0 T f L 2 2
for T   small.
Similarly t E R u = i E R u + i [ , E R ] u + E R Ψ b 1 u + E R b 2 u ¯ + E R f .   Then from Lemma  5.2.6 we have sup 0 < t < T E R u L 2 2 C R N 0 u ( 0 ) L 2 2 + C R N 0 T sup 0 < t < T u ( t ) L 2 2 + C 0 T f L 2 2 .   Now using I = E R + K ~ R ( K R ) *   and the previous estimates we get
sup 0 < t < T u ( t ) L 2 2 C R N 0 sup v L 2 2 + C sup E R u L 2 2
C R N 0 { C R N 0 u ( 0 ) L 2 2 + R M 0 T J 1 / 2 ( u ) x N L 2 2 d t +
+ R T N ( M ) sup 0 < t < T u ( t ) L 2 2 } + C R N 0 T sup 0 < t < T u ( t ) L 2 2 + C 0 T f L 2 2 .
Since N 0   is fixed and M   arbitrary, this combined with Lemma  5.2.6 finishes the proof.
Proof of Theorem  5.1.1 Parts A and B follow by the a priori estimates in Lemma  5.1.3 and Lemma  5.2.1 just as we did in Theorem 4.1 of Section  2 for the elliptic case.
So just part C remains to be proved. It is enough to prove the case s=0 by Lemma  5.1.2 . It can be assumed that u 0 = 0   since (B) solves the case f = 0   and that f S ( R n + 1 )   . Let u   be the solution of ( 5.2 ) given by (B) and let φ S ( R n )   . The family of problems ( 5.2 ) is invariant under adjoints and time reversal, so (B) yields a solution v   of { t v = ( i + Ψ b 1 + b 2 ( x ) ( ¯ ) + Ψ c 1 + Ψ c 2 ( ¯ ) ) * v , v ( x , 0 ) = v 0 ( x ) .   Now
t u ( T ) , v ( T ) = u ( 0 ) , v ( 0 ) + 0 T ( t u ( t ) , v ( t ) + u ( t ) , t v ( t ) ) d t
i u + Ψ b 1 u + b 2 ( x ) u ¯ + Ψ c 1 u + Ψ c 2 u ¯ + f , v
+ 0 T u , ( i + Ψ b 1 + b 2 ( x ) ( ¯ ) + Ψ c 1 + Ψ c 2 ( ¯ ) ) * v d t
= 0 T f , v d t = 0 T x N / 2 J 1 / 2 f , x N / 2 J 1 / 2 v d t
( 0 T R n | J 1 / 2 f | 2 x N d x d t ) 1 / 2 ( 0 T R n | J 1 / 2 v | 2 x N d x d t ) 1 / 2
C ( 0 T R n | J 1 / 2 f | 2 x N d x d t ) 1 / 2 φ L 2 ,
where the last inequality follows from (B). Hence sup [ 0 , T ] u L 2 2 C 0 T R n | J 1 / 2 f | 2 x N d x d t .   Finally use Lemma  5.1.3 , and the proof of Theorem  5.1.1 is complete.

6 NONLINEAR EQUATIONS

The local smoothing results for linear Schrödinger equations will now be applied to obtain local well-posedness in weighted Sobolev spaces with high Sobolev index for a quite general class of nonlinear Schrödinger equations with initial data in the Schwartz class S ( R n )   . This follows the contraction mapping scheme as in [18and [22.

6.1 Linear solutions and weights

Suppose u   is a solution of the linear Schrödinger equation
{ t u = i u + b 1 u + b 2 u ¯ + c 1 u + c 2 u ¯ + f , u ( x , 0 ) = u 0 ( x ) . (6.1)
By Lemma  5.1.2 and lemma  5.1.3 in Section  5 , w = J s u   is a solution of the equation
{ t w = i w + Ψ b 3 w + b 2 w ¯ + Ψ r 1 w + Ψ r 2 w ¯ + J s f w | t = 0 = J s u 0 . (6.2)
where Ψ b 3 = b 1 i s j , k , l x j a k l x j x k x l 3 J 2 , r 1 , r 2 S 1 , 0 0 .   Let W j   denote the solution operator of (6. j   ), j = 1 , 2   , with f = 0   . Then one has J s W 1 u 0 = W 2 J s u 0   . A similar result is needed when J s   is replaced by the weight ( 1 + | x | 2 ) N   . It is useful to obtain a few results concerning commutators of weights and classical pseudo-differential operators.
Lemma 6.1.1 Let p S 1 , 0 m , α N 0 n   . Then x α Ψ p f Ψ p [ x α f ] = 0 < β α ( α β ) Ψ i β ξ β p [ x α β f ] , f S .  
Proof of Lemma  6.1.1 Using integration by parts and Leibniz' rule,
x α Ψ p f ( x ) = x α e i x ξ p ( x , ξ ) f ^ ( ξ ) d ξ
= ( i ξ ) α [ e i x ξ ] p ( x , ξ ) f ^ ( ξ ) d ξ = e i x ξ ( i ξ ) α [ p ( x , ξ ) f ^ ( ξ ) ] d ξ
= e i x ξ β α ( α β ) ( i ξ ) β p ( x , ξ ) ( i ξ ) α β f ^ ( ξ ) d ξ
= e i x ξ β α ( α β ) ( i ξ ) β p ( x , ξ ) ( x α β f ) . ^ ( ξ ) d ξ
= Ψ p [ x α f ] ( x ) + 0 < β α ( α β ) Ψ i β ξ β p [ x α β f ] ( x ) .
This proves Lemma  6.1.1 .
Lemma 6.1.2 Let p S 1 , 0 m , N N   . Then
( 1 + | x | 2 ) N Ψ p f = Ψ p [ ( 1 + | x | 2 ) N f ] + 2 N j Ψ i ξ j p [ x j ( 1 + | x | 2 ) N 1 f ]
+ | α + β | 2 N , | α | 2 , | β | 2 N 2 c α β Ψ ξ α p [ x β f ] , f S ( R n ) .
Proof of Lemma  6.1.2 Lemma  6.1.1 yields the identity ( 1 + | x | 2 ) Ψ p Ψ p ( 1 + | x | 2 ) = 2 j Ψ i ξ j p x j Ψ Δ ξ p .   Lemma  6.1.2 follows by induction of N   and further applications of Lemma  6.1.1 .
Now we study weighted Sobolev norms of solutions of ( 6.1 ) and ( 6.2 ). Let u ( t ) = W 1 ( t ) u 0   be the solution of the linear equation { t u = i u + Ψ b 1 u + b 2 u ¯ + Ψ c 1 u + Ψ c 2 u ¯ u ( x , 0 ) = u 0 ( x ) .  
Lemma 6.1.3 Let N N , s R   . Suppose x 2 N u 0 H s + 2 N .   Then
sup 0 t T x 2 N W 1 ( t ) u 0 H s 2 j = 0 2 N c j T j x 2 N j u 0 H s + j 2 , (6.3)
and
sup 0 t T x 2 N W 1 ( t ) u 0 H s 2 c ( 1 + T 2 N ) x 2 N u 0 H s + 2 N 2 . (6.4)
Proof of Lemma  6.1.3 Let u = W 1 u 0   . Then
t [ x 2 N u ] = x 2 N t u
= x 2 N { i u + Ψ b 1 u + b 2 u ¯ + Ψ c 1 u + Ψ c 2 u ¯ } .
Using Lemma  6.1.2 , this equals to
i x 2 N u + Ψ b 1 x 2 N u + b 2 x 2 N u ¯ + Ψ c ~ 1 x 2 N u +
+ Ψ c ~ 2 x 2 N u ¯ + + 2 N j Ψ ξ j h 2 x j x 2 N 2 u .
Hence x 2 N W 1 u 0   satisfies a linear equation with initial data x 2 N u 0   and forcing term f = 2 N j Ψ ξ j h 2 x j x 2 N 2 W 1 u 0 .   Hence,
sup 0 t T x 2 N W 1 ( t ) u 0 H s 2 c x 2 N u 0 H s 2
+ c 0 T j x j x 2 N 2 W 1 ( t ) u 0 H s + 1 2 d t
c { x 2 N u 0 H s 2 + T j sup 0 t T x j x 2 N 2 W 1 ( t ) u 0 H s + 1 2 } .
Using Lemma  6.1.1 instead of Lemma  6.1.2 and arguing as above, it follows that x j x 2 N 2 W 1 u 0   satisfies a linear equation with initial data x j x 2 N 2 u 0   and forcing term f = Ψ q x 2 N 2 W 1 u 0 ,   where q S 1 , 0 1   . Therefore,
sup 0 t T x j x 2 N 2 W 1 ( t ) u 0 H s + 1 2
c x j x 2 N 1 u 0 H s + 1 2 + c 0 T Ψ q x 2 N 2 W 1 ( t ) u 0 H s + 1 2 d t
c x 2 N 1 u 0 H s + 1 2 + c T sup 0 t T x 2 N 2 W 1 ( t ) u 0 H s + 2 2 ,
and
sup 0 t T x 2 N W 1 ( t ) u 0 H s 2 c x 2 N u 0 H s 2
+ c 1 T x 2 N 1 u 0 H s + 1 2 + c 2 T 2 x 2 N 2 W 1 ( t ) u 0 H s + 2 2 .
Now apply this result N 1   times with N   replaced by N 1 , N 2 , . . . , 1 ,   and the proof of ( 6.3 ) is complete.
For ( 6.4 ), it suffices to note that
x 2 N j u 0 H s + j = x j x 2 N u 0 H s + j
c x 2 N u 0 H s + j c x 2 N u 0 H s + 2 N
if j { 0 , 1 , . . . , 2 N }   since x j S 1 , 0 0   . This proves Lemma  6.1.3 .

6.2 The nonlinear Cauchy problem

Consider the initial value problem
{ t u = i u + b 1 u + b 2 u ¯ + c 1 u + c 2 u ¯ + P ( u , u , u ¯ , u ¯ ) , u ( x , 0 ) = u 0 ( x ) , (6.5)
where u ( x ) = j , k x j ( a j k ( x ) x k u ( x ) ) , A ( x ) = ( a j k ( x ) ) j , k = 1 , . . . , n   is a real, symmetric n × n   matrix, and P   is any polynomial with no linear or constant terms.
Concerning the variable coefficients, it will be assumed that a j k , b 1 , b 2 , c 1 , c 2 C b .   Assume further that the matrix A ( x ) = ( a j k ( x ) ) ) j , k = 1 , . . . , n   is positive definite or invertible. There will be additional hypotheses on a j k , b 1   and b 2   in each of the two cases. More precisely, Elliptic case Suppose ν 1 | ξ | 2 | A ( x ) ξ ξ | ν | ξ | 2 , x , ξ R n .   Then assume in addition the following:
  • (a) A ( x )   generates a bicharacteristic flow with non-trapped bicharacteristics.
  • (b) There exist N > 1   and a constant C   such that if λ ( | x | ) = x N   | a j k ( x ) | , | Im b 1 ( x ) | C λ ( | x | ) , x R n .  
Ultrahyperbolic case Suppose ν 1 | ξ | | A ( x ) ξ | ν | ξ | , x , ξ R n .   Then assume in addition the following:
  • (a) a j k ( x ) a j k 0 S ( R n )   for j , k = 1 , . . , n   where A 0 = ( a j k 0 )   is a real symmetric n × n   constant matrix.
  • (b) The bicharacteristics are non-trapped.
  • (c) b 1 , b 2 S ( R n : C n ) .  
Remark 6.2.1 It will be clear from our previous results and the proof below that assumptions in (a) and (c) in the hyperbolic case can be relaxed so that one only needs a finite number of seminorms in ( 3.2 ) in Section  3 .
Under the above assumptions, the following result holds.
Theorem 6.2.1 Let u 0 S ( R n )   and s n + 4 N + 13   . Then there exists T = T ( u 0 H s , x 2 N J s 3 / 2 u 0 L 2 )   such that ( 6.5 ) 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 ) ) .   Let X T s = { u C ( [ 0 , T ] ; H s ( R n ) ) : max j = 1 , 2 , 3 λ j ( w ) < }   where
λ 1 ( w ) = sup 0 t T w H s
λ 2 = 0 T R n | J s + 1 / 2 w ( x , t ) | 2 x 2 N d x d t
λ 3 ( w ) = sup 0 t T x 2 N t w ( t ) H s / 2 + n / 2 + 3 .
Then for every u 0 S ( R n )   there exists a neighborhood U   of u 0   in S   and a T > 0   such that the map data   solution of ( 6.5 ) is continuous from U   into X T s   .
Remark 6.2.2 The classical pseudo-differential theory in subsection 2.1 of Section  2 and the new operator calculus in Section  3 both basically rely on Taylor expansions of finite order and finitely many integrations by parts. Consequently, by Sobolev's theorem, the assumption u 0 S   in Theorem  6.2.1 can be relaxed to x 2 N J s 1 u 0 L 2   for some large s 1 = s 1 ( n , s )   . The solution u   of ( 6.5 ) is then in C ( [ 0 , T ] : H s ( R n ) ) L 2 ( x N d x ) )   . It is an interesting problem to determine the optimal regularity and decay of u 0   needed in specific examples of equation ( 6.5 ).
Proof of Theorem  6.2.1 Let s 0 2 N   with s 0 + 1 2 n + 4 N + 13   . Let v = J s 0 u   . Then u   solves ( 6.5 ) if and only if v   solves { t v = i v + Ψ b ~ 1 v + b 2 v ¯ + Ψ c ~ 1 v + Ψ c ~ 2 v ¯ + J s 0 [ P ( J s 0 v , J s 0 v , J s 0 v ¯ , J s 0 v ¯ ) ] v ( x , 0 ) = J s 0 u 0 ( x ) .   Using Leibniz' rule,
J s 0 [ P ( J s 0 v , J s 0 v , J s 0 v ¯ , J s 0 v ¯ ) ]
= j Q 1 , j ( J s 0 v , J s 0 v , J s 0 v ¯ , J s 0 v ¯ ) x j v
+ j Q 2 , j ( J s 0 v , J s 0 v , J s 0 v ¯ , J s 0 v ¯ ) x j v ¯
+ R 1 ( { x α J s 0 v , x α J s 0 v ¯ } | α | s 0 / 2 + 1 ) Ψ p 1 v
+ R 2 ( { x α J s 0 v , x α J s 0 v ¯ } | α | s 0 / 2 + 1 ) Ψ p 2 v ¯ ,
where Q 1 , j , Q 2 , j , j = 1 , . . . , n   , and R 1 , R 2   are polynomials with no constant terms and p 1 , p 2 S 1 , 0 0   . The right hand linear factors in the 2 n + 2   terms above arise from the highest order derivative in each term of the Leibniz sum. Now let
Q ~ 1 ( v ) =
j [ Q 1 , j ( J s 0 v , J s 0 v , J s 0 v ¯ , J s 0 v ¯ ) Q 1 , j ( u 0 , u 0 , u 0 ¯ , u 0 ¯ ) ] x j v ,
Q ~ 2 ( v ) =
j [ Q 2 , j ( J s 0 v , J s 0 v , J s 0 v ¯ , J s 0 v ¯ ) Q 2 , j ( u 0 , u 0 , u 0 ¯ , u 0 ¯ ) ] x j v ,
R ~ 1 ( v ) =
[ R 1 ( { x α J s 0 v , x α J s 0 v ¯ } | α | s 0 / 2 + 1 ) R 1 ( { x α u 0 , x α u 0 ¯ } | α | s 0 / 2 + 1 ) ] Ψ p 1 v ,
R ~ 2 ( v ) =
[ R 2 ( { x α J s 0 v , x α J s 0 v ¯ } | α | s 0 / 2 + 1 ) R 2 ( { x α u 0 , x α u 0 ¯ } | α | s 0 / 2 + 1 ) ] Ψ p 2 v ¯ ,
and
b ~ 3 ( x , ξ ) = b ~ 1 ( x , ξ ) + j Q 1 , j ( u 0 , u 0 , u 0 ¯ , u 0 ¯ ) ( x ) i ξ j ,
b ~ 4 ( x ) = b 2 ( x ) + Q 2 ( u 0 , u 0 , u 0 ¯ , u 0 ¯ ) ( x ) ,
c ~ 3 = c ~ 1 + R 1 ( { x α u 0 , x α u 0 ¯ } | α | s 0 / 2 + 1 ) p 1 ,
c ~ 4 = c ~ 2 + R 2 ( { x α u 0 , x α u 0 ¯ } | α | s 0 / 2 + 1 ) p 2 .
Then it suffices to solve the following nonlinear equation for v   , { t v = i v + Ψ b ~ 3 v + b ~ 4 v ¯ + Ψ c ~ 3 v + Ψ c ~ 4 v ¯ + Q ~ 1 ( v ) + Q ~ 2 ( v ) + R ~ 1 ( v ) + R ~ 2 ( v ) v ( x , 0 ) = J s 0 u 0 ( x ) .   This corresponds to solving the integral equation v ( t ) = W 1 ( t ) J s 0 u 0 + 0 t W 1 ( t t ) [ Q ~ 1 ( v ) + Q ~ 2 ( v ) + R ~ 1 ( v ) + R ~ 2 ( v ) ] ( t ) d t ,   where w ( t ) = W 1 ( t ) w 0   is the solution of the linear homogeneous equation { t w = i w + Ψ b ~ 3 w + b ~ 4 w ¯ + Ψ c ~ 3 w + Ψ c ~ 4 w ¯ , w ( x , 0 ) = w 0 ( x ) .   The solution of the integral equation is a fixed point of the following map which will turn out to be a contraction on a suitable function space. Let
[ Φ u 0 ( w ) ] ( t ) = W 1 ( t ) J s 0 u 0
+ 0 t W 1 ( t t ) [ Q ~ 1 ( v ) + Q ~ 2 ( v ) + R ~ 1 ( v ) + R ~ 2 ( v ) ] ( t ) d t ,
λ 1 ( w ) = sup 0 t T w ( t ) H 1 / 2 ,
λ 2 ( w ) = ( 0 T R n | J 1 w ( x , t ) | 2 x 2 N d x d t ) 1 / 2 ,
λ 3 ( w ) = sup 0 t T x 2 N t w ( t ) H s 0 / 2 + n / 2 + 7 / 2 ,
Λ ( w ) = max { λ j ( w ) : j = 1 , 2 , 3 } ,
λ 4 ( w ) = sup 0 t T x 2 N w ( t ) H s 0 / 2 + n / 2 + 7 / 2 ,
X T a = { w : R n × [ 0 , T ] C : Λ ( w ) a } .
For suitable a   and sufficiently small T   it will be shown that Φ u 0   maps the complete metric space X T a   into X T a   and is a contraction.
Let w X T a   . The first goal is to show that Φ u 0 ( w ) X T a   . Notice that v = Φ u 0 ( w )   solves the linear equation { t v = i v + Ψ b ~ 3 v + b ~ 4 v ¯ + Ψ c ~ 3 v + Ψ c ~ 4 v ¯ + Q ~ 1 ( w ) + Q ~ 2 ( w ) + R ~ 1 ( w ) + R ~ 2 ( w ) v ( x , 0 ) = J s 0 u 0 ( x ) .   By the linear smoothing effect,
λ 1 2 ( Φ u 0 ( w ) ) + λ 2 2 ( Φ u 0 ( w ) ) c { u 0 H s 0 + 1 / 2 2
+ 0 T R n | Q ~ 1 ( w ) | 2 x 2 N d x d t + 0 T R n | Q ~ 2 ( w ) | 2 x 2 N d x d t
+ 0 T R n | R ~ 1 ( w ) | 2 x 2 N d x d t + 0 T R n | R ~ 2 ( w ) | 2 x 2 N }
= c { u 0 H s 0 + 1 / 2 2 + I 1 + I 2 + I I 1 + I I 2 } .
Concerning the first term,
I 1 = 0 T R n | Q 1 ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ )
Q 1 ( u 0 , u 0 , u 0 ¯ , u 0 ¯ ) ] w | 2 x 2 N d x d t
c sup 0 t T x 2 N [ Q 1 ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ )
Q 1 ( u 0 , u 0 , u 0 ¯ , u 0 ¯ ) ] L x 2 0 T R n | w | 2 x 2 N d x d t
= c I 1 a I 1 b .
For j { 1 , . . . , n }   one has x j = ( x j J 1 ) J 1   . Observe that x j J 1   is a 0th order classical pseudo-differential operator and is hence bounded on L 2 ( x 2 N d x )   by Lemma  2.3.1 in Section  2 .
Therefore, I 1 b c 0 T R n | J 1 w | 2 x 2 N d x d t = c λ 2 ( w ) .   In order to estimate I 1 a   , let j { 1 , . . . , n }   and t [ 0 , T ]   . Then
Q 1 ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ ) Q 1 ( u 0 , u 0 , u 0 ¯ , u 0 ¯ )
= 0 t t [ Q 1 ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ ) ( t ) ] d t .
By the product rule for t   -differentiation, Lemma  6.1.2 and Sobolev's theorem, I 1 a c T λ 3 2 ( w ) S 1 ( λ 1 2 ( w ) )   where S 1   is some polynomial of one variable. In the rest of the proof, S j , j = 2 , 3 . . .   , will denote other such polynomials. Combining estimates for I 1 a   and I 1 b   , I 1 c T λ 3 2 ( w ) λ 2 2 ( w ) S 1 ( λ 1 2 ( w ) ) .   The estimate for I 2   is similar. Concerning the last two terms,
I I 1 =
0 T R n | [ R 1 ( { x α J s 0 w , x α J s 0 w ¯ } | α | s 0 / 2 + 1 ) R 1 ( { x α u 0 , x α u 0 ¯ } | α | s 0 / 2 + 1 ) ]
Ψ p 1 w | 2 x 2 N d x d t
c T λ 1 2 ( w )
x 2 N [ R 1 ( { x α J s 0 w , x α J s 0 w ¯ } | α | s 0 / 2 + 1 ) R 1 ( { x α u 0 , x α u 0 ¯ } ) | α | s 0 / 2 + 1 ) ] L x 2
= c T λ 1 2 ( w ) I I 1 a .
To estimate I I 1 a   let t [ 0 , T ]   . Then
R 1 ( { x α J s 0 w , x α J s 0 w ¯ } | α | s 0 / 2 + 1 ) ( t ) R 1 ( { x α u 0 , x α u 0 ¯ } | α | s 0 / 2 + 1 )
= 0 t t [ R 1 ( { x α J s 0 w , x α J s 0 w ¯ } | α | s 0 / 2 + 1 ) ( t ) ] d t .
By the product rule for t   -differentiation, Lemma  6.1.2 and Sobolev's theorem, I I 1 a C T λ 3 2 ( w ) S 2 ( λ 1 2 ( w ) ) .   Hence I I 1 C T 2 λ 3 2 ( w ) λ 1 2 S 2 ( λ 1 2 ( w ) ) .   The estimate for I I 2   is similar. Combining estimates for I 1 , I 2 , I I 1   and I I 2   , one gets that ( λ 1 2 + λ 2 2 ) ( Φ u 0 ( w ) ) c u 0 H s 0 + 1 / 2 2 + c T ( 1 + T ) Λ 4 ( w ) S 3 ( Λ 2 ( w ) ) .   Next λ 4   will be estimated and then used in the estimate of λ 3   . By Lemma  6.1.3 (b),
λ 4 ( Φ u 0 ( w ) ) sup 0 t T x 2 N W 1 ( t ) J s 0 u 0 H s 0 / 2 + n / 2 + 7 / 2
+ T sup 0 t t T x 2 N W 1 ( t t ) [ Q ~ 1 ( w ) + Q ~ 2 ( w )
+ R ~ 1 ( w ) + R ~ 2 ( w ) ] ( t ) H s 0 / 2 + n / 2 + 7 / 2
c ( 1 + T N ) x 2 N J s 0 u 0 H 1 + c T ( 1 + T N )
{ sup 0 t T x 2 N Q ~ 1 ( w ( t ) ) H 1 + sup 0 t T x 2 N Q ~ 2 ( w ( t ) ) H 1
+ sup 0 t T x 2 N R ~ 1 ( w ( t ) ) H 1 + sup 0 t T x 2 N R ~ 2 ( w ( t ) ) H 1 }
= c ( 1 + T N ) x 2 N J s 0 u 0 H 1 + c T ( 1 + T N ) { I 1 + I 2 + I I 1 + I I 2 } .
By the product rule for x j   -differentiation and the fundamental theorem of calculus in the t   -variable,
Q ~ 1 ( w ( t ) ) =
j x j [ w ( t ) 0 t t [ Q 1 , j ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ ) ( t ) ] d t ]
+ w ( t ) j 0 t x j t [ Q 1 , j ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ ) ( t ) ] d t
so that, by Sobolev's theorem,
I 1 c j sup 0 t T x 2 N w ( t )
0 t t [ Q 1 , j ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ ) ( t ) ] d t L 2
+ c j sup 0 t T x 2 N w ( t )
0 t x j t [ Q 1 , j ( J s 0 w , J s 0 w , J s 0 w ¯ , J s 0 w ¯ ) ( t ) ] d t L 2
c T λ 3 ( w ) λ 1 ( w ) S 4 ( λ 1 ( w ) ) .
Similarly, I 2 C T λ 3 ( w ) λ 1 ( w ) S 5 ( λ 1 ( w ) ) .   Concerning I I 1   , one has R ~ 1 ( w ( t ) ) = Ψ p 1 w ( t ) 0 t t [ R 1 ( { x α J s 0 w , x α J s 0 w ¯ } | α | s 0 / 2 + 1 ) ( t ) ] d t .   By Sobolev's theorem and the L 2   boundedness of Ψ p 1   ,
I I 1 = sup 0 t T x 2 N Ψ p 1 w ( t ) 0 t t [ R 1 ( { x α J s 0 w , x α J s 0 w ¯ } | α | s 0 / 2 + 1 ) ( t ) ] d t L 2
c T λ 3 ( w ) λ 1 ( w ) S 6 ( λ 1 ( w ) ) .
Similarly I I 2 c T λ 3 ( w ) λ 1 ( w ) S 7 ( λ 1 ( w ) ) .   Combining estimates for I 1 , I 2 , I I 1   and I I 2   , it follows that λ 4 ( Φ u 0 ( w ) ) c ( 1 + T N ) x 2 N J s 0 u 0 H 1 2 + c T 2 ( 1 + T N ) Λ 2 ( w ) S 8 ( Λ ( w ) ) .   This estimate of λ 4   will be used in that of λ 3   . Let v = Φ u 0 ( w )   and note that t v = ( i + Ψ b ~ 3 + Ψ c ~ 3 ) v + ( b ~ 4 + Ψ c ~ 4 ) v ¯ + Q ~ 1 ( w ) + Q ~ 2 ( w ) + R ~ 1 ( w ) + R ~ 2 ( w )   where i + Ψ b ~ 3 + Ψ c ~ 3   and b ~ 4 + Ψ c ~ 4   are classical pseudo-differential operators of order 2 and 1, respectively. By Lemma  6.1.2 ,
λ 3 ( v ) = λ 4 ( t v ) c sup 0 t T x 2 N v H s 0 / 2 + n / 2 + 11 / 2
+ c sup 0 t T x 2 N [ Q ~ 1 ( w ) + Q ~ 2 ( w ) + R ~ 1 ( w ) + R ~ 2 ( w ) ] H s 0 / 2 + n / 2 + 7 / 2 .
The first term can be estimated as λ 4 ( v )   since s 0   is sufficiently large. The second term is dominated by I 1 + I 2 + I I 1 + I I 2   in the estimate of λ 4 ( v )   . Hence λ ( Φ u 0 ( w ) ) c ( 1 + T N ) x 2 N J s 0 u 0 H 1 2 + c T ( 1 + T N + 1 ) Λ 2 ( w ) S 9 ( Λ ( w ) ) .   T   will later be chosen small, so it can be assumed that T 1   . Combining estimates for λ 1 , λ 2   and λ 3   ,
Λ ( Φ u 0 ( w ) ) c ( u 0 H s 0 + 1 / 2 + x 2 N J s 0 u 0 H 1 )
+ c T 1 / 2 Λ 2 ( w ) ( 1 + Λ ρ ( w ) )
where ρ > 0   . First fix a > 2 C ( u 0 H s 0 + 1 / 2 + x 2 N J s 0 u 0 H 1 ) .   Next choose T = min { ( 2 c a ( 1 + a ρ ) ) 2 , 1 } .   Then Λ ( Φ u 0 ( w ) ) a   so Φ u 0   maps X T a   into X T a   . To see that Φ u 0   is a contraction on X T a   for sufficiently small T   , notice that
( Φ u 0 ( w 1 ) Φ u 0 ( w 2 ) ) ( t ) = 0 t W 1 ( t t ) [ ( Q ~ 1 ( w 1 ) Q ~ 1 ( w 2 ) )
+ ( Q ~ 2 ( w 1 ) Q ~ 2 ( w 2 ) ) + ( R ~ 1 ( w 1 ) R ~ 1 ( w 2 ) ) + ( R ~ 2 ( w 1 ) R ~ 2 ( w 2 ) ) ] ( t ) d t .
The estimates used in showing Φ u 0 : X T a X T a   therefore give Λ ( Φ u 0 ( w 1 ) Φ u 0 ( w 2 ) ) C T 1 / 2 Λ ( w 1 w 2 ) a ( 1 + a ρ ) .   Now choose T = min { 1 , 2 C a ( 1 + a ρ ) ) 2 }   . Then Λ ( Φ u 0 ( w 1 ) Φ u 0 ( w 2 ) ) 1 2 Λ ( w 1 w 2 ) ,   so Φ u 0   is a contraction. By Banach's contraction mapping principle there is a unique fixed point of Φ u 0   which solves the nonlinear equation. Now let u   and v   be solutions of ( 6.5 ) with initial values u 0   and v 0   respectively. Then the estimates above give Λ ( u v ) C ( u 0 v 0 H s + x 2 N J s 3 / 2 ( u 0 v 0 ) L 2 ) + C T 1 / 2 Λ ( u v )   if u   and v   are both in an open ball in X T s   . Here C   depends on the radius of the ball. Now choose T > 0   so small that C ( T ) 1 / 2 1 2   . The proof of Theorem  6.2.1 is complete since H s + x 2 N J s 3 / 2 L 2   is dominated by finitely many seminorms in S   (see ( 3.2 )).
References

  1. Calderón, A.P., and Vaillancourt,R., On the boundedness of pseudo-differential operators, J. Math. Soc. Japan 23 (1971), 374–378.
  2. Chihara, H., Local existence for semilinear Schrödinger equations, Math. Japan 42 (1995), 35–51.
  3. Constantin, P., and Saut, J. C., Local smoothing properties of dispersive equations J. Amer. Math. Soc. 1 (1989), 413–446.
  4. Craig, W., Kappeler, T., and Strauss, W., Microlocal dispersive smoothing for the Schrödinger equation, Comm. Pure Appl. Math. 48 (1995), 769–860.
  5. Davey, A., and Stewartson, K., On three–dimensional packets of surface waves Proc. Roy. Soc. London Ser. 338 (1974), 101–110.
  6. Djordjevic, V. D., and Redekopp, L. G., On two-dimensional packets of capillary-gravity waves, J. Fluid Mech. 79 (1977), 703-714.
  7. Doi, S., On the Cauchy problem for Schrödinger type equations and the regularity of solutions, J. Math. Kyoto Univ. 34 (1994), 319–328.
  8. Doi, S., Remarks on the Cauchy problem for Schrödinger–type equations, Comm. Partial Differential Equations 21 (1996), 163–178.
  9. Doi, S. , Smoothing effects for Schrodinger evolution equation and global behavior of geodesic flow Math.-Ann. 318 (2000), 355–389.
  10. Hayashi, N., Global existence of small analytic solutions to nonlinear Schrödinger equations, Duke Math. J. 62 (1991), 575–592.
  11. Hayashi, N., and Ozawa, T., Remarks on nonlinear Schrödinger equations in one space dimension, Differential Integral Equations 7 (1994), 453–461.
  12. Hörmander, L., Pseudo-differential operators and non–elliptic boundary problems, Ann. of Math. (2) 83 (1966), 129–209.
  13. Ichinose, W., On L 2   well-posedness of the Cauchy problem for Schrödinger type equations on a Riemannian manifold and Maslov theory, Duke Math. J. 56 (1988), 549–588.
  14. Ishimori, Y., Multi vortex solutions of a two dimensional nonlinear wave equation, Progr. Theor. Phys. 72 (1984), 33–37.
  15. Kato, T., Quasi–linear equations of evolution, with applications to partial differential equations, Lecture Notes in Math. 448 (1975), Springer–Verlag, 27–70.
  16. Kato, T., 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. Kenig, C. E., Ponce, G., and Vega, L., Oscillatory integrals and regularity of dispersive equations, Indiana University Math. J. 40 (1991), 33–69.
  18. Kenig, C. E., Ponce, G., and Vega, L., Small solutions to nonlinear Schrödinger equations, Ann. Inst. H. Poincaré Anal. Non Linéaire 10 (1993), 255–288.
  19. Kenig, C. E., Ponce, G., and Vega, L., On the Zakharov and Zakharov–Schulman systems, J. Funct. Anal. 127 (1995), 204–234.
  20. Kenig, C. E., Ponce, G., and Vega, L., On the smoothing properties of some dispersive hyperbolic systems, Nonlinear waves(Sapporo, 1995), 221-229: GAKUTO Internat. ser. Math. Sci. Appl.,10, Tokyo (1997).
  21. Kenig, C. E., Ponce, G., and Vega, L., On the Cauchy problem for linear Schrödinger systems with variable coefficient lower order terms, Harmonic analysis and number theory (Montreal, PQ, 1996), 205–22: CMS Conf. Proc., 21, Amer. Math. Soc., Providence, RI (1997).
  22. Kenig, C. E., Ponce, G., and Vega, L., Smoothing effects and local existence theory for the generalized nonlinear Schrödinger equations, Invent. Math. 134 (1998), 489-545.
  23. Klainerman, S., Long time behavior of solutions to nonlinear evolution equations, Arch. Ration. Mech. and Analysis 78 (1981), 73–98.
  24. Kruzhkov, D. J., and Faminskii, A. V., Generalized solutions for the Cauchy problem for the Korteweg-de Vries equation, Math. USSR Sb 48 (1990), 93–138.
  25. Kumano–Go, H., Pseudo-differential operators, MIT Press, Cambridge (1981).
  26. Mizohata, S., On the Cauchy Problem, Notes and Reports in Mathematics in Science and Engineering,3, Science press and Academic Press (1985).
  27. Rolvung, C., Non-isotropic Schrödinger equations PhD. dissertation, University of Chicago (1998).
  28. Simon, J., and Taflin, E., Wave operators and analytic solutions for systems of nonlinear Klein-Gordon equations and of non-linear Schrödinger equations Comm. Math. Phys. 99 (1985), 541–562.
  29. Sjölin, P., Regularity of solutions to the Schrödinger equations, Duke Math. J. 55 (1987), 699–715.
  30. Stein, E. M., Harmonic analysis: real–variable methods, orthogonality, and oscillatory integrals, Princeton University Press (1993).
  31. Vega, L. , The Schrödinger equation: pointwise convergence to the initial date, Proc. Amer. Math. Soc. 102 (1988), 874–878.
  32. Zakharov, V. E., and Schulman, E. I., Small solutions to nonlinear Schrödinger equations, Physica 1D (1980), 185–250.