Carleman estimates for third order operators of KdV and non KdV-type and applications

In this paper we study a class of variable coefficient third order partial differential operators on $\mathbb{R}^{n+1}$, containing, as a subclass, some variable coefficient operators of KdV-type in any space dimension. For such a class, as well as for the adjoint class, we obtain a Carleman estimate and the local solvability at any point of $\mathbb{R}^{n+1}$. A discussion of possible applications in the context of dispersive equations is provided.


Introduction
In this paper, we will continue the investigation of some variable coefficient PDOs (partial differential operators) built from a system of real smooth vector fields, initiated in the works [18], [19], [13] and [15] (see also [17], [12], [14]).Since when the celebrated works by Kolmogorov [33] first, and by H örmander [25] afterwards, about the hypoellipticity of operators written as sums of squares of vector fields were published, a lot of connected problems have been investigated.In particular, H örmander's hypoelliptic theorem in [25] opened up the study of sub-Laplacians on Lie groups (see [42]), of parametrices for such operators, of sharp estimates in appropriate functional spaces (see the pioneering works [23] and [42]), but also, among other questions, that of the local solvability of closely related models (see, for instance, [41], [37], [18], [19], [13], [15] and references therein).
In the series of works [18], [19], [13] and [15], the authors focus on the local solvability of some classes of degenerate second order PDOs built from a system of real smooth vector fields, plus a somewhere vanishing function.The use of the vanishing function is twofold: on one side it adds degeneracy to a model which is already degenerate by itself, on the other side it permits to include in the treatment some operators like the Kannai operator, that is operators having a changing sign principal symbol.The introduction of the vanishing function requires a price to pay, price which is given in terms of conditions on the lower order terms of the operators (see [17], [12], [14] for an overview of such results).Nevertheless, the classes studied in these works encompass different lower order terms, which makes it possible to include parabolic-type operator, Schr ödinger-type operators, a blend of the two, and some prototypes with non-smooth coefficients.
What we do in the present paper, is to take inspiration by these works and by the form of some operators relevant in physics, and build a class of operators of order three (instead of two) starting from a system of smooth real vector fields.Let us say right away that we will not involve any additional somewhere vanishing function, since the model under consideration is already highly degenerate and complicated without this adjustment.
Besides the connection with the aforementioned previous works, the interest for such kind of operators has several motivations (see [32] and references therein).One of them is represented by the applications to variable coefficient dispersive equations of KdV-type, where we refer to KdV-type operators as those of the form i∂ t + L(t, x, D x ), with L being a third order PDO with smooth coefficients.
In the last decades variable coefficient Schr ödinger equations have attracted lots of attention.Smoothing and Strichartz estimates have been proved under different hypotheses (see [2], [43], [36], [31], [21], [22], [20], [11], [39] and references therein), nonlinear problems have been solved (see [36], [31], [21], [22], [20], [11]), and uniqueness results have been proved (see [10], [8], [7] and [29] in the constant coefficient case, and [3], [16] and references therein for variable coefficient cases).For KdV-type equations the investigation has not been pushed that far, possibly because of the unknown real analogue of this equation in dimensions higher than two, and also because variable coefficient third order equations can be much more challenging to study.Here we wish to start a local analysis of some variable coefficient cases, and prove Carleman type estimates to be employed to reach uniqueness results.Let us underline that such estimates have a pivotal role in establishing uniqueness properties, not only at a local (in space) level, but also for global results.In the context of dispersive equations, results for KdV and ZK (Zakharov-Kuznetsov) equations can be found, for instance, in [4], [9] and [46], while for nonlocal operators we refer the interested reader to [29] and [40].Note that all these cases deal with constant coefficients operators, while here we are concerned with the situation where the coefficients are variable.As far as the author knows, there are no results in the space variable coefficients setting.Recently, some results on pevolution equations in one space dimension have appeared in [1].In that work the operators under investigation have time-dependent coefficients in the leading part, and the class only includes KdV operators (since the space domain is one dimensional) with time variable coefficients in the leading part, but not ZK operators.Here we are concerned with space-dependent coefficients in any space dimension.
An other motivation driving this work, is that of the local solvability problem for multiple characteristics PDOs.This problem is very hard to attack, and general results often require precise geometrical conditions on the characteristic set.We refer the interested reader to [44] for a nice overview of the classical local solvability problem, and to [5] for the resolution of the Nirenberg-Treves local solvability conjecture.
Let us give a description of the family of operators we will be considering in this paper.
Given a system of vector fields {iX j (x, D)} N j=0 on R n+1 -where N is a positive integer not necessarily equal to n or n + 1, and D = (D 1 , . . ., D n+1 ) = (−i∂ x 1 , . . ., −i∂ x n+1 ) -such that • the vector fields iX j (x, D), j = 0, . . ., N, are smooth and real; • the vector fields iX j (x, D); j = 1, . . ., N, form a global involutive distribution, that is, [X l , X m ](x, D) = −[iX l , iX m ] is uniquely determined by a linear combination with smooth coefficients of the vector fields {iX j (x, D)} j=1,...,N , for all l, m = 1, . . ., N, and for all x ∈ R n+1 ; • X 1 is nondegenerate, namely it is nowhere vanishing on R n+1 ; we define P 1 as the third order PDO of the form The latter is the precise class of operators we will be dealing with, which, depending on N and on the choice of the system of vector fields, includes (linear) KdV-type operators of different kind.
We have already mentioned that our goal is to prove Carleman estimates for such class of operators.To be precise, we will not only obtain Carleman estimates for the class of operators represented by P 1 , but also for the adjoint class, that is the family of operators P * 1 obtained by taking the adjoint of P 1 , namely which coincides with P 1 up to a differential operator of order two.
The main motivation for the Carleman estimate for P * 1 , is our interest in the local solvability of P 1 .We recall that the role of Carleman estimates -which is well-known to be determinant in the analysis of unique continuation problems -is also crucial to prove local solvability properties of partial differential equations with variable coefficients.An estimate of this sort for an operator P, yields a local solvability result for the adjoint operator P * .Hence, we can take advantage of the Carleman estimate for P * 1 to get the local solvability of P 1 .The connection between our operator P 1 and some important operators coming from physics, can be seen by the analyzing the following two examples: the (linear) KdV operator and its higher dimensional generalization, the ZK operator.The (linear) KdV operator is described by iP 1 with x , iX 0 = ∂ t .and iX 1 = −∂ x .However, by taking iX 1 (x) = a(x)∂ x , with a(x) 0 for all x ∈ R, then we get an operator, again of the form P 1 , that describes a variable coefficient KdV operator.If we consider x , iX 0 = ∂ t and iX j = −∂ x j , j = 1, 2, then iP 1 describes the ZK operator.Again, by taking nonvanishing space-variable coefficient real vector fields satisfying condition (H 1 ) in (2.1), we get a variable coefficient ZK model.
Note that both the classical KdV and the ZK operator (with constant coefficients) are built by taking an elliptic operator in R n x , i.e. the Laplacian in R and R 2 respectively, and then using the two vector fields iX 0 = ∂ t and iX 1 = ∂ x 1 to define a third order operator.Observe that P 1 is built exactly the same way, but employing variable coefficient real vector fields instead of constant coefficient ones.In fact, to define P 1 we take operators being sums of squares, that is N j=1 X * j X j (not necessarily elliptic in a subspace of dimension n of the n + 1-dimensional domain), and cook up a third order operator by using X 1 and X 0 .An example of variable coefficient KdV-type operator in a three-dimensional Euclidean space domain, can be constructed, for instance, by using the canonical basis of the Lie algebra of the Heisenberg group.Let be four real smooth vector fields, where iX 1 , iX 2 are the generators of the stratified Lie algebra is an elliptic operator on R 3 x , and is a variable coefficient KdV-type operator of the form P 1 for every choice of j 0 = 1, 2, 3, since X j 0 is nondegenerate for all j 0 = 1, 2, The examples above show that P 1 -type operators include KdV and ZK variable coefficient PDOs.However, our class is not limited to be just a generalization of KdV-type operators in any space dimension, since we are free to choose n ≥ 3, N ≤ n + 1, but also iX 0 ∂ t , giving rise to a very wide family of degenerate operators.
Let us finally conclude this introduction with the plan of the paper.Section 2 is devoted to the the main results of this work, specifically Theorem 2.8 about a Carleman estimate for P 1 and for P * 1 , and Theorem 2.10 giving a local solvability result for both the aforementioned operators.Section 3 is dedicated to explicit examples of operators of the form under study and of KdV-type.The last section, Section 4, contains final remarks and the discussion of possible applications in the context of dispersive equations.We recall once more that in this section we are concerned with the proofs of the main results of the article: a Carleman estimate for P 1 and P * 1 , and a local solvability result for these operators.The Carleman estimate is the most important result of the paper, in that it is the key tool to get the solvability property.Moreover, it can find applications in other problems of our interest, like in uniqueness/unique continuation problems for dispersive equations.On the other hand, the investigation of the local solvability of highly degenerate operators is quite complicated, thus our results represent a step ahead in the understanding of certain kind of operators.
Before stating and proving our theorems, let us make precise the objects we are working with.Recall that, on R n+1 x , x = (x 1 , . . ., x n , x n+1 ), n ≥ 1, we define the operator P 1 as where The condition we shall assume on the system of vector fields on R n+1 which make up P 1 is the following: (H 1 ) We say that a system of real smooth vector fields {iX l } l=1,...,N , 1 ≤ N ≤ n + 1, satisfies condition (H 1 ) if the distribution Γ({iX l } l=1,...,N ) is a global involutive distribution, that is, for all j, k = 1, . . ., N, there exist (uniquely determined) where [X, Y] := XY − YX denotes the commutator of the operators X and Y.
In other words (H 1 ) amounts to the global involutivity of the system of real smooth vector fields {iX j } N j=1 on R n+1 .Remark 2.1.Due to the global involutivity of the system, and the fact that N ≥ 1, there exists at least one vector field iX j 0 ∈ {iX j } N j=1 being nondegenerate on R n+1 , or, in other words, a nowhere vanishing vector field on R n+1 .By possibly renaming the vector fields, we can always assume j 0 = 1, so we will use this convention in the rest of the paper when assuming condition (H 1 ) on a system of vector fields.This assumption will be crucial to derive the Carleman estimate below, and thus the solvability result as well.Note also that, if N ≤ n + 1, and Γ({iX j }) is a global involutive distribution of rank = N, then X j is nondegenerate in R n+1 for all j = 1, . . ., N, hence we are free to take any vector field in the system as the one playing the role of X 1 .Finally, observe that the requirement N ≤ n + 1 is not restrictive.Due to the global involutivity, and the fact that the rank of the distribution is M ≤ n + 1, then if N > n + 1 we can rewrite P 1 in terms of M ≤ n + 1 < N vector fields in the second order part.Remark 2.2.Notice that formula (2.1) holds true independently of the rank of the distribution.If the system of vector fields satisfies (H 1 ), and the rank of where some coefficients will be identically zero.Since we will not care about the precise rank of the distribution, we will simply use (2.1).
Remark 2.3.An additional remark is that, by imposing only (H 1 ) on the system of vector fields, then P 1 generates a class of operators larger than the one of KdV-type operators, in that the latter is modeled by using a system of vector fields satisfying the additional conditions Since we will be giving a local solvability result for P 1 and P * 1 , for completeness we recall here the definition of H s − H s ′ locally solvable partial differential operator at a point x 0 ∈ Ω, where, adopting the usual notation, H s (Ω) stands for the Sobolev space of order s.Definition 2.4.Let P be a partial differential operator with smooth coefficients defined on an open subset Ω of R n , and let x 0 ∈ Ω.We say that P is H s − H s ′ locally solvable at x 0 if there exists a compact K containing x 0 in its interior U K , such that, for every f ∈ H −s loc (Ω) there exists u ∈ H −s ′ loc (Ω) for which Pu = f in U K .We say that P is H s − H s ′ locally solvable in Ω, if it is locally solvable at each point of Ω.
When P is H s − H s ′ locally solvable at x 0 with s = s ′ = 0, we will just say that the operator is L 2 − L 2 locally solvable at x 0 .
Notations.Below N will always be a positive integer, λ will be a real number, whereas A = {a ji (x)} j=0,...,N i=1,...,n+1 will denote the variable coefficient matrix containing the coefficients of the vector fields iX j , for j = 0, . . ., N. For a matrix A as above and α ∈ N n+1 {0} n+1 , we define Throughout the paper we will write • H s without specifying the set Ω where the L ∞ (Ω), L 2 (Ω), and the H s (Ω) norms are taken.Given that we will be working with compactly supported functions, the set Ω will be the fixed support of the function (or its interior), so we omit it for simplicity.However, note that we can take Finally, we shall use the notation Γ({X j } N j=1 ) for the distribution spanned by the system of vector fields {X j } N j=1 .
In order to prove our Carleman inequality for P 1 and for P * 1 , we will make use of the subsequent fundamental lemma.Lemma 2.5.Let f ∈ C ∞ (R n+1 ), N ≥ 1, and λ ≥ 1.Let also P 1 be as in (1.1) and {iX j } N j=1 satisfying (H 1 ).Then, for every compact set K of R n+1 , there exists a positive constant where with c j1 k , as in (2.1), and c Proof.We start the proof by computing the quantity where L 2 , L 1 , L 0 are the partial differential differential operators of order 2, 1, and 0, given, respectively, by so, for simplicity, we handle each term on the right-hand side separately.
Recall that, since the system of vector fields {iX j } j=1,...,N is globally involutive, that is satisfies (H 1 ), for all j = 1, . . ., N, there exist By the same token, there exist some Now, on putting v := e −λ f u, and using the previous observation, we have Adopting the notation , while C K denotes a constant depending on the L ∞ (K)-norm of some smooth functions, specifically functions depending on the coefficients of the vector fields iX j and on their derivatives.Abusing of the notation, below we shall write C K for any constant of this form, that is any constant depending on the L ∞ (K)-norms of f , of the coefficients of the vector fields iX j , j = 0, . . ., N, and of the derivatives of these functions.Sometimes C K will also depend on the L ∞ (K)-norm of c kl j , k, l = 1, . . ., N, j = . . ., n + 1.
As for the term containing L 2 , recalling that f is a real function and that X 1 f takes purely imaginary values, we obtain (2.4) By similar considerations and Finally, by (2.3), (2.4), (2.5), (2.6), and the fact that λ ≥ 1, we can find a new suitable positive constant C K such that with This completes the proof.
The same result as Lemma 2.5, with suitable small adjustments, is still valid for the operator P * 1 , the adjoint of P 1 .We state the precise result in Lemma 2.6 below.Lemma 2.6.Let f ∈ C ∞ (R n+1 ), N ≥ 1, and λ ≥ 1.Let also P * 1 be as in (1.2), and {iX j } N j=1 satisfying (H 1 ).Then, for every compact set K of R n+1 , there exists a positive constant C Im(e −λ f P * 1 u, e −λ f u) where , iR) for all j, k = 1, . . ., N, and Q j , c ′ j , as in (2.2) for all j = 1, . . ., N.
Proof.To prove this lemma we make use of Lemma 2.5.By standard arguments we have Next, as in the proof of Lemma 2.5, we compute where .
Using the same notations as in the proof of Lemma 2.5, we have where Q j is as in (2.2).
Next, we define g j (x) := N k=1 (c )), and have that, since g j and (X 1 d j ) are real valued functions, and since the c j1 k take purely imaginary values, then and where Inserting the last three inequalities, that is (2.10), (2.9) and (2.8), into (2.7), and using that λ ≥ 1, by rearranging the indices we conclude that is Q ′ j differs from Q j in the coefficients of the second term in (2.11), and in the appearance of the function c ′ j − c ′′ j instead of c ′ j in the last term in (2.11).This completes the proof.Remark 2.7.The estimates proved so far do not necessitate of the nondegeneracy of X 1 .The global nonvanishing requirement on X 1 will come into play in the subsequent results.
In the next theorem we will refer to the so called H örmander's condition of rank r on a system of real smooth vector fields.Let us recall that this means that the system of vector fields, along with their commutators of length r (being assumed that ]],... are counted as commutators of length 1, 2, 3,.. and so on), generate the whole tangent space at each point, namely, in our case, R n+1 .Theorem 2.8 (Carleman estimate for P 1 and P * 1 ).Let f ∈ C ∞ (R n+1 ), {iX i } N i=1 be a system of real smooth vector fields satisfying (H 1 ), and K a compact subset of R n+1 .Then, (ii) If the system of vector fields {iX i } n i=1 satisfies Hörmander's condition of rank r, and if there exists Proof.The proof relies on the use of Lemma 2.5 and of Lemma 2.6 when dealing with the statement for P 1 and for P * 1 respectively.Below we shall give a detailed proof of the result for P 1 , and we omit the one for P * 1 .This is done because the proofs of the estimates for P * 1 and for P 1 are exactly the same, up to the appearance of different constants depending on the L ∞ (K)-norm of some smooth functions (because the difference between the estimates in Lemma 2.5 and Lemma 2.6 is just in the appearance of possibly different smooth coefficients in Q ′ j and Q j ).Since there is no substantial change in the two cases, we focus on P 1 .By Lemma 2.5, for all λ ≥ 1, given a compact K of R n on which −iX 1 f ≥ C 0 > 0, we have Im(e −λ f P 1 u, e −λ f u) with where ic j1 k , ic ′ k , k = 1, . . ., N, are smooth and real valued functions.Therefore, once λ is big enough, the proof of (i) and of (ii) depends on the estimate we can prove for the second order term n j=1 Re(Q j e −λ f u, e −λ f u).Then, let us focus on this term and show how different hypotheses lead to different lower bounds.
In the sequel we will use again the notation v := e −λ f u.First observe that, by Cauchy-Schwarz inequality, ), and with a new positive constant C 1 := max j,k=1,...,N {|c j1 k |}, and where C K is a positive constant depending on all the L ∞ -norms of the functions X j d 1 and where C K is a different suitable positive constant depending on the L ∞ -norms of f, a jk , c ′ j and of their derivatives.Therefore, big enough, and δ 0 , δ ′ 0 , δ 1 sufficiently small, so that where C K is a new positive constant since we have fixed δ 0 , δ ′ 0 , δ 1 .The estimate above yields where, recall, C K is a constant depending on some fixed constants N, δ 0 , δ ′ 0 , δ 1 and on the L ∞ -norms on K of f, a jk , c ′ j , and of their derivatives.
Finally, by Cauchy-Schwarz inequality, which gives, taking δ small enough and not necessarily depending on K (since we can always assume C K 2 • λ > 1 and take, for instance, δ = 1/2 < λ 2 C K 4 ), we reach which proves part (i) of the theorem.
To prove part (ii) of the statement, we apply the Rothschield-Stein inequality for systems of vector fields satisfying H örmander's condition of rank r on the first term on the right hand side of (2.12), that is This completes the proof.
Remark 2.9 (The case N = n = 1).Note that when N = n = 1 (H 1 ) is trivially satisfied.For the same reason, result (ii) does not apply to this case, since X 1 alone cannot generate R 2 .
An immediate application of the Carleman estimates in Theorem 2.8 is the following local solvability result.Theorem 2.10.Let {iX j } N j=1 be a system of vector fields satisfying (H 1 ) such that X 1 (x, D) 0 for all x ∈ R n+1 .Let also P 1 and P * 1 be as in (1.1) and (1.2) respectively.Then the following properties hold.
(i) P 1 and P * 1 are L 2 − L 2 locally solvable at any point of R n+1 .(ii) If {iX j } N j=1 satisfies Hörmander's condition of rank r, P 1 and P * 1 are H −1/r − L 2 locally solvable at any point of R n+1 .
Proof.The proof is based on the use of Theorem 2.8 and of some standard functional analysis argument.More specifically, the H s − L 2 local solvability of P 1 at a point x ∈ R n+1 , is equivalent to the validity of the following solvability estimate for the adjoint P * 1 : there exists C > 0 and a compact K containing x in its interior U K , such that The solvability estimate, Hahn-Banach and Riesz representation theorems, alltogether give the result (see, for instance, Lemma 1.2.30 in [35]).Since X 1 is nondegenerate (see also Remark 2.1), given C 0 > 0 and x 0 ∈ R n+1 , it is always possible to find a smooth function f and a compact set K ⊂ R n+1 containing x 0 in its interior, such that iX 1 f (x) ≥ C 0 > 0 for all x ∈ K. Hence, by the Carleman estimate for P * 1 , for x 0 ∈ R n+1 , and for the compact K chosen as above and containing x 0 in its interior, there exists λ 0 (K) > 1 and C > 0 such that, for all λ > λ 0 , where s = 0 in case (i) and s = −1/r in case (ii).By estimating from above and from below e −λ f on K with positive constants, (2.13) easily implies the solvability estimate with s = 0 under the hypothesis in part (i), and with s = −1/r under the hypothesis in part (ii).From the solvability estimate, by the standard considerations mentioned at the beginning of the proof (see [35]) we get the L 2 − L 2 and H − 1 r − L 2 local solvability of P 1 at x 0 according to the hypotheses in (i) and (ii) respectively.Finally, by the arbitrariness of x 0 , we conclude the local solvability of P 1 at any point of R n+1 .
For the local solvability of P * 1 , one repeats the same steps reversing the roles of P 1 and P * 1 .Remark 2.11.When the vector field iX 1 is nondegenerate only on a compact K ⊂ R n+1 , using the previous strategy we can reach the local solvability of P 1 and of P * 1 at any point of the interior U K of K.

Examples and applications: KdV and non KdV-type operators with variable coefficients
In this section we would like to list a few examples of operators of P 1 -type.We pointed out more than once that we are particularly interested in dispersive models, therefore we will mainly focus on such cases here, even though the class is more general and we could exhibit different type of operators.

Variable coefficient KdV operators. Let
where a ∈ C ∞ (R x ; R), a(x) 0 for all x ∈ R. Then which is a KdV operator with space-variable coefficients.Our solvability theorem applies to this model, so we can assert that this operator, as well as iP 1 of course, is L 2 − L 2 locally solvable at any point of R 2 .More generally, one can consider the coefficient a as a function of both (t, x), hence a(t, x), or even just a time dependent function, provided that tit is everywhere nonvanishing.In this last case one can also describe some of the operators treated in [1].

Variable coefficient ZK operators. Let
x , and X 0 (t, x, D t , D x ) = D t , and X j (t, x, D t , D x ) = X j (x, D x ), j = 1, 2, with X 1 , X 2 such that X * 1 X 1 (x, D x ) + X 2 * X 2 (x, D) is an elliptic operator on R 2 x .In other words, we are assuming the vector fields iX j (x, D) to be nondegenerate and linearly independent at any point x ∈ R 2 x , so they form a global involutive distribution both in R 2 and in R 3 .Under these assumptions the operator ) is a prototype of a variable coefficient ZK-operator to which our solvability theorem applies.Hence, we have the L 2 − L 2 local solvability property for P 1 in R n+1 by Theorem 2.10.
Notice that the "ellipticity in space" of X * 1 X 1 (x, D x ) + X * 2 X 2 (x, D) implies global involutivity, thus any prototype of this sort fits in our class.On the other hand, we assumed the ellipticity in space of X * 1 X 1 (x, D x ) + X * 2 X 2 (x, D) just in order to represent a variable coefficient analogue of the ZK operator, and not because P 1 -type operators need to satisfy this requirement.In fact, we can take more general cases, provided that the global involutivity is satisfied.In general, to describe a variable coefficient P 1 -type ZK operator, one can take X 1 , X 2 of the form x , with a ij ∈ C ∞ (R 2 ; R) for all i, j = 1, 2, where the global involutivity is given by the third condition.In such a case our solvability theorem still applies.
Non KdV-type operators.To give a more complete picture of the objects of this paper, we provide immediate examples of non KdV-type operators of P 1 -type.
Let N ≤ n and R n+1 = R t ×R x .Let also {iX j } N j=1 be a global involutive system of real smooth vector fields in R n x (depending only on the space variables and on space derivatives), X 1 nondegenerate on R n+1 , and X 0 ≡ 0. Then is of P 1 -type both on R n and on R n+1 but not of KdV-type on R n+1 , since there is no time-evolution here.Note that we could also take N ≤ n + 1 and define an operator whose leading part has variable coefficients depending both on space and time.
More degenerate variable coefficient KdV-type operators.Still in the case n = 2, one can consider models where N = 1, iX 1 (x, D t , D x ) = X 1 (x, D x ) is a nonvanishing real vector field with smooth coefficients, and X 0 = D t .In such a case we have no ellipticity in space for X * 1 X 1 , since X * 1 X 1 cannot be elliptic under our assumption that iX 1 has real coefficients.Nevertheless, all the hypotheses, namely (H1) and the nondegeneracy of X 1 , are satisfied, hence our results apply -the Carleman and the solvability theorem -and the operator is locally solvable at any point of R 3 .
This example can be generalized in any dimension n ≥ 2. By taking {iX j (x, D x )} N j=1 , N ≤ n, being a globally involutive system of real smooth vector fields (indexed in such a way that X 1 is nondegenerate), and taking iX 0 = ∂ t , the corresponding operator P 1 is built from an operator L := N j=1 X * j X j which is not necessarily elliptic in space.
In the rest of this section we will give concrete examples of KdV and non KdV P 1 -type operators constructed via some Lie algebras, specifically stratified Lie algebras.These cases are taken into account due to their global involutive structure.One can, of course, design many other different examples.
KdV-type operators built via the Heisenberg Lie algebra.Let us start with the example mentioned in the introduction and related to the Heisenberg group.We restrict ourselves to N = n = 3 -hence P 1 is defined on R 4 t,x -which can be immediately generalized to the case N = n = 2k + 1, for any positive integer k.
We take, as before, X 0 = D t and iX j , j = 1, 2, 3, as the vector fields giving the canonical basis of the Lie algebra of H 1 .Recall that the Lie algebra h 1 is stratified of step 2, and that X 1 , X 2 are the so-called "generators" of the stratified Lie algebra.These vector fields, together with their commutator, generate the whole Lie algebra h 3 .This guarantees that the vector fields X j , j = 1, 2, 3, form a global involutive structure.For completeness, we recall the expression of the vector fields below Note that X j = X * j , thus P 1 = P * 1 in this case.Moreover, X j (x, D) is nondegenerate for all j = 0, . . ., 3, implying that the operators P 1, j 0 (t, x, D t , D x ) = X j 0 (x, D x ) 3 j=1 X 2 j (x, D x ) + D t , j 0 = 1, 2, 3, (t, x) ∈ R 4 , are locally solvable at any point of R 4 t,x by Theorem 2.10.To define the same kind of operators in higher dimension it suffices to use the Heisenberg Lie algebra h k of dimension 2k + 1, and take N = 2k + 1 = n.
KdV-type operators built via the Heisenberg Lie algebra in higher dimensional spaces.By taking the same vector fields as above generating h 1 , one can cook up more singular operators in higher dimensional spaces.Take N = 3 < n, R n+1 = R t × R 3 x × R n−3 y , and P 1, j 0 as before, but now defined on R n+1 , that is P 1, j 0 (t, x, y, D t , D x , D y ) = X j 0 (x, D x ) 3 j=1 X 2 j (x, D x ) + D t , j 0 = 1, 2, 3, (t, x, y) ∈ R t × R 3 x × R n−3 y .
Then we have that the system of vector fields {iX j (x, D x )} j=1,2,3 is still globally involutive in R n+1 t,x,y , and the hypotheses of Theorem 2.10 are satisfied.Hence, we get local solvability for these operators at any point of R n+1 t,x,y .Similar operators can be written by using the Heisenberg Lie algebra h k , for every k ≥ 1, and N = 2k + 1 < n.
KdV-type operators built via stratified Lie algebras: general construction.Following the previous constructions, we can formalize the procedure to build operators as above with any stratified That said, it should be possible to reach global (in all but at most one direction) Carleman estimates for P 1 and P * 1 from our local estimates in Theorem 2.8.These global estimates can be applied to study uniqueness problems for dispersive operators belonging to our class.For what we just explained above, by introducing suitable conditions on the coefficients of the operator and on the weight function f , one should be able to produce the desired global inequality.However, even if the road map to get a global inequality is in the proof of Theorem 2.8, its derivation and the application to uniqueness problems is nontrivial.
Let us also stress that, if one is interested in the study of local unique continuation properties, then our estimates are already strong enough to pursue this goal, and one should focus more on the study of the geometric conditions leading to the result, that is, in other words, on the choice of the suitable weight function f .Other KdV-type models and relative questions.In our discussion, we have motivated our analysis with the connection of the class P 1 with KdV-type operators.However, it must be said that we have a rigorous derivation of the ZK equation as a sort of generalization of the KdV one in (space) dimension 2, but for the higher dimensional case the argument is much harder and not established to the author knowledge.On the other hand, it is known that third order equations play a central role in water waves, nonlinear optic and related fields.A detailed interesting description of third order dispersive equations and of their physical role is contained in [32], where a class of third order operators with constant coefficients is studied.Let us remark that a subclass of the class in [32] of constant coefficient operators is included in our general variable coefficient class P 1 .That said, it would be interesting to investigate Carleman estimates, local solvability, and connected problems, for operators in the very general form Note that this class, which generalizes and contains P 1 , clearly still contains KdV-type operators, both with constant and with variable coefficients, as well as the whole class studied in [32].Due to these considerations, it would be interesting to understand weather P 2 is the right generalization of KdV-type operators, or if P 1 provides the best description, being clear that P 2 is a class wider than P 1 .In any case, the validity of Carleman estimates in this generality is still an open problem.

2 .
Carleman estimate and local solvability for P 1 and P * 1