EXISTENCE AND UNIQUENESS OF VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM FOR THE p th SUB-LAPLACIAN ASSOCIATED TO A SYSTEM OF HÖRMANDER VECTOR FIELDS arxiv:808.00849v [math.ap] 2 Aug 208 DUY-MINH NHIEU Abstract. We establish the existence and uniqueness of variational solution to the nonlinear Neumann boundary problem for the p th -Sub- Laplacian associated to a system of Hörmander vector fields. Introduction Two of the fundamental boundary value problems in PDE are the Dirchlet and Neumann problems. They have profound influences in the development of the theory of PDE. Aside from their theoretical applications, they describe many physical models. In the classical setting, these two problems have been advanced to a far extend. We would not attempt to give a bibliography here since we would inadvertently omit some of the important ones. Boundary value problems for Sub-elliptic operators were initiated by the pioneering works of Kohn and Nirenberg, Bony, Gaveau and complemented later by two important works of Jerison. These works established some fundamental aspects of subelliptic equations such as the Harnack inequality, smoothness of the Green s function up to the boundary, for boundaries that are nowhere characteristic with respect to the operator [KN],[B], explicit formula for the harmonic measure at the center of the Koranyi ball in the Heisenberg group H n [G] and the behavior of solutions to the Dirichlet problem near characteristic points [J, J2]. Of course, we only highlight those that are closely related to our development in this paper. For other early relevant results, the reader can see [OL, OL2, Ba, De]. These pioneering works subsequently drawed a significance of attention to boundary value problems in the subelliptic setting, see for instance the works [D],[D2], [HH], [Ci], [Da], [UL], [CG], [GV], [CGN], [CGN2], [GNg] and the references therein. However, none of the works cited above addressed the Neumann problem. Our purpose is to establish existence and uniqueness of variational solutions of the Neumann problem associated to the non-linear Sub-elliptic p th -Laplacian arising from a system of Hörmander vector fields. This is a continuation and extension of the works [N] and [DGN]. To set the stage we fix an open set U R n with diam(u) < and X = (X,...,X m )asystemofvector fieldswithsmoothcoefficients satisfying Hörmander s finite rank condition at every point x U: Supported by the Ministry of Science and Technology, R.O.C., Grant MOST 06-25- M-008-03.
2 DUY-MINH NHIEU ranklie[x,...,x m ] n Fix p <. Let L p be the p th -Sub-Laplacian associated to the system of vector fields X = (X,...,X m ) defined on smooth functions u by L p u = m Xj( Xu p 2 X j u) j= on R n. In the above the X j = n k= b(j) k (x) x k are smooth vector fields on R n (satisfying Hörmander s finiterankcondition though such an assumption is not necessary to introduce the operator) and Xj is the formal adjoint of X j given by X j g(x) = n Xu = k= x k (b (j) k (x)g(x)). We have also indicate by ( m ) j= (X ju) 2 2. For a smooth (at least C 2 ) domain, the classical Neumann problem for L p consists of seeking a function u that can be differentiated enough times and satisfying (.) { L p u = f in, m j= Xu p 2 (X j η)x j u = ν on where η is the unit outer (Euclidean) normal of, ν and f are continuous functions satisfy the compatibility condition (.2) < ν, > = f dx. When ν = g for some function g, then (.2) simply becomes g dµ = f dx. In the above, the notation <,> indicated the pairing between an element x Z andν Z for somebanach spacez andits dualz (weusemake use of a more general description of the left hand side of (.2)) in preparation for a more abstract statement later). For the sake of convenience, we would also indicate the sum m j= (X ju)(x j v) by the notation < Xu,Xv >, as long as it is clear from the context. In the case p = 2, the operator L p is known as the (real part of the) Kohn-Laplacian. Our results here extend the ones established for p = 2 in [N],[DGN]. In the classical setting, the method of layer potentials has met success in tackling such problems. However, even in the simplest prototype of the Hörmander type vector field, that is, the Heisenberg group, one faces significant difficulties in inverting the corresponding operators on the boundary of the domain (for an appropriate class of functions) due to the presence of characteristic points, see [J]. Besides, layer potentials will not help in dealing with the case p 2. In view of this obstacle, we seek a different method to establish the problem of existence of the non-linear Neumann problem in this broad context. Oneapproach to the problem (.) is to broaden the class and the concept of solutions in such a way that existence can be established. One would then
VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM 3 establish that solutions are regular in some context. It is our purpose to explore the first aspect in this paper and devote the second one to another occasion [DGMN]. We now introduce variational solutions. Definition.. Fix p <, q <, f L q (,dx), ν B q s q(,dµ) where /q +/q =. A function u L,q (,dx) is called a variational solution of (.) if for any φ L,q (,dx) (.3) Xu p 2 < Xu,Xφ > dx = < ν,tr(φ) > f φdx In the above, R n is a domain for which there is a continuous (trace) operator Tr : L,q (,dx) B q s q(,dµ). In the above, L,q (,dx) denotes the by now standard subelliptic Sobolev space associated to the system of vector fields X. For the definition and some properties of the Besov spaces B q s q(,dµ) we refer the reader to Section 2 below. The measure µ is an s-ahlfors measure (again see section 2 for the definition) supported in with 0 < s < min(q, n+q 2 ) Remark.2. In order for Definition. to make sense, it is required that the continous operator T r exists and this was established in [DGN, Theorem 0.6] under the assumptions. is an extension domain (i.e. X (ǫ,δ)-domain) and 2. µ satisfies the estimate (see definition 2.) (.4) µ(b(x,r)) M B(x,r) r s. In particular,. and 2. are fulfilled when = {ρ < } where ρ is the homogenous gauge in a Carnot group of step two and when dµ = Xρ dσ where dσ is the surface measure on, and s = in (.4), see e.g. [CG], [CGN] and [CGN2]. Therefore, in the sequel we will work under such hypothesis. Whenever an expression involves function values on the boundary, we mean the trace of such function. When the context is clear will omit the Tr notation. 3. By taking the test function φ(x), we arrive at the compatibility condition (.2), which is necessary for the existence of solutions to the above variational Neumann problem. The close subspace of L,q (,dx) defined by { L,q (,dx) = f L,q (,dx) f = is an appropriate space to treat the Neuamann problem. } f(x)dx = 0
4 DUY-MINH NHIEU Remark.3. Due to Poincaré inequality (see Theorem 2.9 below), an equivalent norm on L,q () is given by f L,q () = ( Xf q dx ) q. Our aim is to establish the existence and uniqueness of variational solution in such spaces. Our main result is the following Theorem.4. Let U be as in Theorem 2.8. Fix < p < and < q p. Let f L q (,dx), U R n be an X (ǫ,δ) domain with diam() < R o /2, > 0, µ an upper s-ahlfors measure and ν B q s/q (,dµ) satisfying the compatibility condition (.2). There exists a unique u L,p (,dx) that solves the Neumann problem for L p in the sense of Definition.. In the above, the parameter R o is from definition 2. below. Furthermore, the solution u satisfies the estimate (.5) Xu L p (,dx) C( ν B q + f s/q (,dµ) L q (,dx) ) p, where C depends on and various parameters such as R o, p,q and so on but does not depend on u,ν,f. To establish this main result, we show that the functional J p defined in (3.) has a unique minimum and it is the variational solution of (.). Conversely, variational solutions of (.) are minimizers of J p. The proof of this fact follows the classical approach of the direct method of Calculus of Variations. First, we show that the functional J p is sequentially lower semicontinuous in the weak topology of L,q (,dx). This is valid for all q <. We then establish the coercivity of J p in L,q (,dx). At the moment, this is limited to the range < q p. A weakness as it may seems in comparison to its Euclidean counterpart. However, in the Euclidean setting, for the case p = 2 on Lipschitz domains, similar result only holds for < q < 2 + ǫ where ǫ depends on the domain. The class of domains (the X (ǫ,δ) domains) that we treat here have characteristic points on the boundary. Such singular points act as corners (or perhaps as cusps) of the domain if seen from the Euclidean perspective. It is well-known that in the Euclidean case the inclusion of the classes of domains (.6) Lipschitz NTA (ǫ,δ) hold. The so call (ǫ,δ) domains were introduce by [Jo] and they are the largest class on which Sobolev functions can be extended. It is important to observe th at for a Carnot-Carathéodory space the last inclusion of (.6) continue to hold, see [CT]. A version of Jerison and Kenig s theory on the boundary behaviour of harmonic functions [JK] has been extended to the case of Carnot groups of step two [CG] on X NTA domains. Contrary to the Euclidean case where Lipschitz domains are probably the largest class of domains where a rich theory of boundary value problems can be developed, the analogue notion of Lipschitz domain in the subelliptic setting even in
VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM 5 the simplest case of the Heisenberg group is almost non-existence, see [CG]. Therefore part of our task is to identify a class of domains, as large as possible for which a theory of boundary value problems can be developed. The examples in [CT] and [CG] are an indication that the class of X (ǫ,δ) domains treated here is reasonably large. Due to the inclusion (.6), the domains considered in this paper, when confined to the Euclidean case, is larger than the Lipschitz domains. Hence, it is not supprising if a limitation on the parameter q is occurs. However, we have not been able to determine the sharp range at this moment. We start with some background materials in section 2 and section 3 is devoted in establishing our main result. Finally, we like to point out that, to the best of our knowledge, the results here are new even in the classical setting, as far as the class of domains treated is concerned. 2. Preliminaries We collect previously established results needed for subsequent developments. Some of these results are more general than what we require here. We present the version that is already adapted to our setting and omit the details that are relevant in a more broader context. Let d be the Carnot-Carathéodory distance associated with the system X. It is by now well known that if the system X satisfies Hörmander s finite rank condition [H], then d(x,y) < for any x,y R n ([NSW], [Chw], [Rs]) and the metric balls B(x,r) of d satisfy a doubling condition [NSW]. Denote the Borel measures on the metric space (R n,d) by B d. For a set E R n, E denotes the Lebesgue measure of E. Definition 2. ([DGN]). Given s 0, a measure µ B d will be called an uppers-ahlfors measure, if there exist M,R o > 0, such that for x R n, 0 < r R o, one has µ(b(x,r)) M B(x,r) r s. We will say that µ is a lower s-ahlfors measure, if for some M,R o > 0 one has instead for x and r as above B(x,r) µ(b(x,r)) M r s. The (dual of the) following subelliptic Besov space will serve as the space for the Neumann datum on the boundary. Definition 2.2 ([DGN]). Let µ B d having supp µ F, where F is a closed subset of R n. For p <, 0 < β <, we introduce the semi-norm N p β (f,f,dµ) = { F F ( ) f(x) f(y) p d(x,y) s } d(x,y) β B(x,d(x,y)) dµ(y)dµ(x) p. The Besov space on F, relative to the measure µ, is defined as
6 DUY-MINH NHIEU B p β (F,dµ) = {f Lp (F,dµ) N p β (f,f,dµ) < }. If f B p β (F,dµ), we define the Besov norm of f as f B p β (F,dµ) = f L p (F,dµ) + N p β (f,f,dµ). We note in passing the following result [DGN, Theorem.], which motivates the use of Besov spaces (since their dual are larger) instead of Lebesgue spaces on the boundary for the Neumann data. However, our main result here does not make use of Theorem 2.3 below. Note also that since we only needed µ to be an upper s-ahlfors measure, a comparison between Besov and Lebegues spaces as candidates for the boundary datum also require that we assume µ to be a lower s-ahlfors measure as well. In the subelliptic settings, ample supply of measures that satisfy both an upper and lower s-ahlfors condition are found in [DGN] and the references therein. Theorem 2.3 (Embedding a Besov space into a Lebesgue space). Given a bounded set U R n having characteristic local parameters C,R o, and local homogeneous dimension Q, let U be an open set with diam < R o /2. Let p, 0 < β <. Suppose µ is a lower s-ahlfors measure with 0 < s n+βp, s < Q βp, and such that supp µ = F. There exists a continuous embedding B p β (F,dµ) Lq (,dµ), and, in fact, for f B p β (F,dµ) one has f L q (,dµ) C where q = p Q s Q s βp, ) } {( + diam()β µ(f) β/(q s) N p β (f,f,dµ) + µ(f) β/(q s) f L p (,dµ), where C = C(,C,R o,p,β,s,m) > 0. Furthermore, ( f(x) f,µ q dµ(x) ( C F ) q where f,µ denotes the average µ() F f(x) f(y) p d(x,y)s β p f dµ. ) B(x,d(x,y)) dµ(y)dµ(x) p, A wide class of domains to which our results hold is the following: Definition 2.4 ([N, DGN]). An open set R n is called an X (ǫ,δ)- domain if there exist 0 < δ, 0 < ǫ, such that for any pair of points p,q, if d(p,q) δ, then one can find a continuous, rectifiable curve γ : [0,T], for which γ(0) = p, γ(t) = q, and
VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM 7 l(γ) ǫ d(p,q), d(z, ) ǫ min {d(p,z),d(z,q)} for all z {γ}. The following compact embedding theorem [GN, Theorem.28] plays an important role in our proof. Theorem 2.5 (Compact embedding). Let R n be an X (ǫ,δ) domain with diam() < R 0 2. Then, one has the following: (I) The embedding BV X (,dx) L q (,dx) is compact for any q < Q Q. (II) For any p < Q the embedding L,p (,dx) L q (,dx) is compact provided that q < Qp Q p. (III) For any Q p < and any q <, the embedding L,p (,dx) L q (,dx) is compact. As an easy consequence of the above Theorem we have Theorem 2.6. With the same assumptions in Theorem 2.5, for any p < the embedding L,p (,dx) L p (,dx) is compact. We can deduce from Theorem 2.6 easily the following Corollary 2.7. For any < p <, if {u h } h= is a sequence in L,p (,dx) such that u h u for some u L,p (,dx) then u h u in L p (,dx) hence u h u in L (,dx) also. (We note explicitly the notation means weak convergence whereas means convergence in norm). The proof of this fact is rather standard but we include it here for the sake of convenience of the reader. Proof. The weak convergence assumption together with the Banach-Steinhaus (uniform boundedness principle) theorem implies that {u h } h= is bounded in L,p (,dx). Now Theorem 2.6 implies that if {u hj } is any subsequence of {u h }, {u hj } has a subsequence we denote it by {v hj } such that v hj v L p (,dx) for some function v. Since u h u in L p (,dx), we must have u = v andtheoriginal sequence{u h }converge innormtouinl p (,dx). The traces of Sobolev functions on the boundary of a domain is a delicate matter. It was the purpose of [DGN] to develop such a theory in the setting of a Carnot-Carathéodory space. It is also an indispensible tool in dealing with the Neumann problem, even the non-linear version. We thus recall [DGN, Theorem 0.6]: Theorem 2.8 (Trace theorem on the boundary). Let U R n be a bounded set with characteristic local parameters C,R o, and let p >. There is σ = σ(x,u) > 0 such that, if U is a bounded X (ǫ,δ)-domain with rad() > 0,diam() < Ro 2σ, dist(, U) > R o, and µ is an upper s- Ahlfors measure for some 0 < s < p, having supp µ, then for every 0 < β s/p there exist a linear operator
8 DUY-MINH NHIEU T r : L,p (,dx) B p β (,dµ), and a constant C = C(U,p,s,M,β,ǫ,δ,rad()) > 0, such that (2.) Trf B p β (,dµ) C f L,p (,dx). Furthermore, if f C () L,p (,dx), then T rf = f on. Our final pillar is the following result [GN, Corollary.6] specialized to our setting: Theorem 2.9. Let be an X (ǫ,δ) domain, p < Q. There exists some constant C = C(ǫ,δ,p,n) > 0 such that for any k Q Q p and any u L,p () ( ) ( ) u u kp kp dx C diam() Xu p p dx 3. Existence and Uniqueness of solutions to the Neumann Problem Let be an X (ǫ,δ) domain. Consider the functional J p : L,q (,dx) R given by (3.) J p (u) = p Xu(x) p +f(x)u(x)dx < ν,tr(u) >. Note that in appriori, unless q = p or u C () L.q (,dx), the functional J p can take on the value of ± and the assumption that is an X (ǫ,δ) domain is necessary for the trace operator tr (as in Definition.) to be defined. Hence, this assumption is needed in the above definition and throughout the paper. Lemma 3.. Let be an X (ǫ,δ) domain. For any q the functional J p given by (3.) is convex on L,q (,dx). If in addition we assume that > 0 and p > then the functional J p is strictly convex in the subspace L,q (,dx). Proof. It suffices to establish the proposition for the non-linear part of J p, namely (3.2) I(u) = p Xu(x) p dx. It is an elementary fact that the function g : R m R given by g(z) = z p is convex for p < (see e.g., [Roc, Theorem 5.] and the remark there). This implies that any functional of the form F(u) = g(xu)(x)dx is convex. Hence, I is convex. We turn to the second part of the proposition. For p > the function g is strictly convex. Hence for any 0 < t < and any u,v L,q (,dx) we have g(t(xu) + ( t)(xv)) < tg(xu) ( t)g(xv)
VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM 9 unless Xu = Xv, that is, for all i =,..,m, X i (u v) = 0. To continue, observe that the system of Hörmander vector field X = (X,...,X m ) satisfy the property that i =,..m, for a.e.x : X i g(x) = 0 implies g(x) = constant for a.e.x. This imply u = v +c for some constant c. Since u,v L,p (,dx), taking into account that > 0, we see that c must be zero. Remark 3.2. Two remarks are in order. (i) The uniqueness of minimizer would follow from the strict convexity of J p, that is, minimizer would be unique in the spaces L,q (,dx). (ii) Due to the compatibility condition (.2), we have for any constant c R and any u L,q (,dx), J p (u+c) = J p (u). Hence, if a (unique) minimizer u o L,q (,dx) exists for J p, then u o remains to be a minimizer for J p in L,q (,dx) since for any u L,q (,dx), J p (u) = J p (u u ) J p (u o ) and that for any constant c R, u o +c is also a minimizer for J p in L,q (,dx). Next, weturntothe(sequential) lower semicontinuity of thefunctional J p. The following result (in many equivalent form) from Calculus of Variations is valuable to us and is available from many sources. For the sake of our purpose, we apply the following theorem from [Gi, Theorem 4.5]. Theorem 3.3. Let be an open set in R n, M a closed set in R N, and let F(x,u,z) be a function defined in M R l such that (i) F is a Caratheodory function, that is measurable in x for every (u,z) M R l and continuous in (u,z) for almost every x (ii) F(x,u,z) is convex in z for almost every x and for every u M. (iii) F 0. Let u h,u L (,M,dx),z h,z L (,R l ) and assume that u h u and z h z in L loc (,dx). Then, F(x,u,z)dx liminf F(x,u h,z h )dx h Corollary 3.4. Let be an X (ǫ,δ) domain with diam() < Ro 2 (and f L q (,dx)). The functional J p defined in (3.) is (sequentially) lower semicontinuous on the spaces L,q (,dx) for any q < with respect to the weak topology, that is, if {u h } h= is a sequence in L,q (,dx) such that u h u for some u L,q (,dx) then J p (u) liminf h J p(u h ). Proof. Byassumptionu L,q (,dx),f L q (,dx)andν B q s q(,dµ) the linear part of J is bounded by Hölder s inequality and the trace inequality (2.) therefore continuous. The sum of two lower semicontinuous function
0 DUY-MINH NHIEU is lower semicontinuous. Hence, it suffices to establish the lower semicontinuity of the non-linear part of J p, namely I defined in (3.2). Let u,{u h } h= be as in the hypothesis. By Corollary 2.7 we have u h u in L (,dx) and also X j u h X j u in L q loc (,dx) hence also in L (,dx) for j =,..m. Now we apply Theorem 3.3 with N =,M = R,l = m, z = Xu = (X u,...,x m u) and F(x,u,z) = z p to reach the conclusion. Lemma 3.5. Let q <, f L q (,dx), be X (ǫ,δ) domain. A function u L,q (,dx) minimizes the functional J p given by (3.) if and only if it is a variational solution to the Neumann problem (.) in the sense of Definition.. Proof. Let u L,q (,dx) be a minimizer of J p. We have for any ǫ > 0 and any φ L,q (,dx) the function ψ(ǫ) = J p (u+ǫφ) reaches a local minimum at ǫ = 0. Therefore ψ (0) = 0. On the other hand, differentiating under the integral sign yields ψ (ǫ) = Hence 0 = ψ (0) = Xu+ǫXφ p 2 < Xu+ǫXφ,Xφ > +f(x)φ(x)dx < ν,φ >. Xu p 2 < Xu,Xφ > + f(x)φ(x)dx < ν,φ >, which is (.3) in Definition.. We now prove the converse. Let u L,q (,dx) be a variational solution in the sense of Definition.. For any v L,q (,dx) we let φ = u v L,q (,dx) in (.3), adding and subtracting the terms p Xv p dx we obtain the identity J p (u) = J p (v) = J p (v) p Xu p dx, p Xv p + p Xu p Xu p 2 < Xu,Xv > dx p Xv p p Xu p Xu p 2 < Xu,Xv Xu > dx. To show that the last term in the above is non-negative, we recall that the convexity of the function g(z) = z p for p can be rephrase as: For every z,w R m, z 0 then (3.3) p w p p z p + z p 2 m z i (w i z i ). Letting z i = X i u and w i = X i v in the inequality (3.3) and integrating we have p Xv p p Xu p Xu p 2 < Xu,Xv Xu > dx 0. Hence, we have J p (u) J p (v) and this completes the proof. i=
VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM Our final result leading to the existence of minimizers of J p (hence, variational solutions to the Neumann problem in L,q (,dx)) is the following Lemma 3.6. Let be an X (ǫ,δ) domain. For any q p and f L q (,dx), the functional J p is coercive in L,q (,dx), that is, if {u h } is a sequence in L,q (,dx) suchthat lim h u h L,q (,dx) = lim h Xu h L q (,dx) = then lim h J p (u h ) =. Proof. Observe that due to the Trace theorem 2.8 and the fact that f L q (,dx), the linear part of J p is bounded. Therefore, it suffices to establish the Lemma for the non-linear part of J p, namely I given by (3.2). Let {u h } be a sequence in L,q (,dx) such that lim h Xu h L q (,dx) =. By Hölder s inequality, we have for q p that Xu h L q (,dx) q p Xu h L p (,dx) hence lim h I(u h ) = lim h Xu h L p (,dx) = as well. We now come to our main result which follows from a standard line of argument of the Calculus of Variations, Lemma 3.6 and Corollary 3.4. We include the proof for the sake of completeness. Theorem 3.7. Let < q p, be an X (ǫ,δ) domain with diam() < R o /2, > 0, f L q (,dx), ν B q s q(,dµ) satisfying the compatibility condition (.2). The functional J p given by (3.) has a unique minimizer in L,q (,dx). Proof. Let l = inf{j p (u) u L,q (,dx)} > since the linear part of J p is bounded and the non-linear part I given by (3.2) is non-negative. Let {u h } L,q (,dx) be a minimizing sequence, that is lim h J p (u h ) = l. Clearly, J p (u h ) is bounded and hence Lemma 3.6 implies that {u h } is a bounded sequence in L,q (,dx). Since for q >, the balls in L,q (,dx) are weakly compact, {u h } contains a subsequence (still denoted by {u h }) and u L,q (,dx) such that u h u o. Corollary 3.4 then imply l J p (u o ) liminf h J p(u h ) = lim h u h = l, hence J p (u o ) = l. Now from Remark 3.2, u o remains to be a minimizer of J p in L,q (,dx). Finally, we now come to the proof of our main result. Proof of Theorem.4. Theorem 3.7 and Lemma 3.5 yield a unique solution u L,q (,dx) in the sense of Definition.. To establish the estimate (.5) we take φ = u in Definition. and recalling Remark.3 we have
2 DUY-MINH NHIEU Xu p L p (,dx) = Xu p 2 < Xu,Xu > dx < ν,tr(u) > + f udx ν B q s/q (,dµ) tr(u) B q s/q (,dµ) + f L q (,dx) u L q (,dx) (by Theorems 2.8,2.9) C ν B q s/q (,dµ) u L,q (,dx) + f L q (,dx) Xu L q (,dx) (by Remark.3) C( ν B q s/q (,dµ) + f L q (,dx) ) Xu L q (,dx) (since q p) C( ν B q s/q (,dµ) + f L q (,dx) ) Xu L p (,dx) and therefore Xu L p (,dx) C( ν B q s/q (,dµ) + f L q (,dx) ) p. In view of Remark.3 again this shows u L,p (,dx) and the proof is now complete. Remark 3.8. Some intermediate results such as Lemmas 3., 3.5 and Corollary 3.4 hold for any q. The limitation of < q p comes in Lemma 3.6 and therefore also in Theorem.4. The case for q < p can also be infered from the case q = p due to the containment of the spaces B q s/q (,dµ) B p s/p (,dµ) and L,p (,dx) L,q (,dx). References [Ba] M. S. Baouendi, Sur une classe dopétrateurs elliptiques dégénérés, Bull. Soc. Math. Fr., 95 (967), 45-87. [B] J. M. Bony, Principe du maximum, inégalité de Harnack et unicité du problème de Cauchy pour les operateurs elliptique degeneres, Ann. Inst. Fourier, Grenoble () 9 (969), 277304. [CG] L. Capogna& N. Garofalo, Boundary behavior of non-negative solutions of subelliptic equations in NTA domains for Carnot-Carathéodory metrics, Journal of Fourier Anal. and Appl., 4-5, 4 (998), 403-432. [CGN] L. Capogna, N. Garofalo & D. M. Nhieu, Properties of harmonic measures in the Dirichlet problem for nilpotent Lie groups of Heisenberg type, Amer. J. Math., 24 (2002), no. 2, 273306. [CGN2] Mutual absolute continuity of harmonic and surface measures for Hrmander type operators, Perspectives in partial differential equations, harmonic analysis and applications, Proc. Sympos. Pure Math., 79 (2008), 4900. [CT] L. Capogna, P. Tang, Uniform domains and quasiconformal mappings on the Heisenberg group, Manuscripta Math., 86 (995), no. 3, 26728. [Ci] G. Citti, Wiener estimates at boundary points for Hrmanders operators, Boll. Un. Mat. Ital. 2-B, (988), 66768.
VARIATIONAL SOLUTION TO THE NEUMANN PROBLEM 3 [Chw] W. L. Chow, Über Systeme von linearen partiellen Differentialgleichungen erster Ordnung, Math. Ann., 7 (939), 98-05. [D] E. Damek, Harmonic functions on semidirect extensions of type H nilpotent groups, Trans. Amer. Math. Soc., 290 (985), 375384. [D2] A Poisson kernel on Heisenberg type nilpotent groups, Colloq. Math., 53 (987), 239247. [Da] D. Danielli, Regularity at the boundary for solutions of nonlinear subelliptic equations, Indiana Univ. Math. J., 44 (995), 269286. [DGMN] D. Danielli, N. Garofalo, I. Munive & D. M. Nhieu, The Neumann problem on the Heisenberg group, regularity and boundary behaviour of solutions, in preparation. [De] M. Derridj, Un probleme aux limites pour une classe dopérateurs hypeolliptiques du second ordre, Ann. Inst. Fourier (Grenoble) 2, 4 (97), 99-48. [GNg] N. Garofalo & C. P. Nguyen, Boundary behavior of p-harmonic functions in the Heisenberg group, Math. Ann. 35 (20), 587632. [GV] N. Garofalo & D. Vassilev, Regularity near the characteristic set in the non-linear Dirichlet problem and conformal geometry of sub-laplacians on Carnot groups, Math. Ann. 38 (2000), 453-56. [GN] N. Garofalo & D.M. Nhieu, Isoperimetric and Sobolev inequalities for Carnot- Carathéodory spaces and the existence of minimal surfaces, Comm. Pure Appl. Math., 49 (996), 08-44. [G] B. Gaveau, Principe de moindre action, propagation de la chaleur et estimees sous elliptiques sur certain groupes nilpotents, Acta Math. 39 (977), 9553. [DGN] D. Donatella, N. Garofalo and D.M. Nhieu, Non-doubling Ahlfors measures, perimeter measures, and the characterization of the trace spaces of Sobolev functions in Carnot-Carathéodory spaces, Mem. Amer. Math. Soc.,82, (2006), no. 857. (998). [Gi] E. Giusti, Direct methods in the calculus of variations. World Scientific Publishing Co., Inc., River Edge, NJ, 2003. ISBN: 98-238-043-4. [HH] W. Hansen & H. Hueber, The Dirichlet Problem for Sublaplacians on Nilpotent Lie Groups Geometric Criteria for Regularity, Math. Ann. 276 (987), 537-547. [H] L. Hörmander, Hypoelliptic second order differential equations, Acta Math., 9 (967), 47-7. [J] D. S. Jerison, The Dirichlet problem for the Kohn Laplacian on the Heisenberg group. I, J. Funct. Anal., 43 (98), no., 97-42. [J2] The Dirichlet problem for the Kohn Laplacian on the Heisenberg group. II, J. Funct. Anal., 43 (98), no. 2, 224-257. [JK] D. Jerison & C.Kenig, Boundary behavior of harmonic functions in non-tangentially accessible domains, Adv. in Math., 46 (982), 80-47. [Jo] P. Jones, Quasiconformal mappings and extendability of functions in Sobolev spaces, Acta Math., 47 (98), no. -2, 788. [KN] J. J. Kohn & L. Nirenberg, Non-coercive boundary value problems, Comm. Pure and Appl. Math. 8 (965), 443492. [NSW] A. Nagel, E. M. Stein & S. Wainger, Balls and metrics defined by vector fields I: basic properties, Acta Math., 55 (985), 03-47. [N] D. M. Nhieu, The Neumann problem for sub-laplacians on Carnot groups and the extension theorem for Sobolev spaces, Ann. Mat. Pura Appl. (4) 80 (200), no., 25. [OL] O. A. OLEINIK, On the smoothness of solutions of degenerate elliptic and parabolic equations, Dokl. Akad. Nauk SSSR, 63 (965), 577-580. [Russian] [OL2] On linear second order equations with non-negative characteristic form, Mat. Sb., 69 (966), I-40. [Russian] [Rs] P. K. Rashevsky, Any two points of a totally nonholonomic space may be connected by an admissible line, Uch. Zap. Ped. Inst. im. Liebknechta, Ser. Phys. Math., (Russian) 2 (938), 83-94. [Roc] R. T. Rockafellar, Convex analysis. Princeton Mathematical Series, No. 28 Princeton University Press, Princeton, N.J. 970. [UL] F. Uguzzoni, E. Lanconelli On the Poisson kernel for the Kohn Laplacian. Rend. Mat. Appl. (7) 7 (997), no. 4, 659677.
4 DUY-MINH NHIEU Department of Mathematics, National Central University, Zhongli District, Taoyuan City 3200,, Taiwan, R.O.C. E-mail address, Duy-Minh Nhieu: dmnhieu@math.ncu.edu.tw