$\text{1}$

## On Hessian measures for non-commuting vector fields.

### November 27, 2006

Abstract
Previous results on Hessian measures by Trudinger and Wang are extended to the subelliptic case. Specifically we prove the weak continuity of the 2-Hessian operator, with respect to local ${L}^{1}$  convergence, for a system of $m$  vector fields of step 2 and derive gradient estimates for the corresponding $k$  -convex functions, $1\le k\le m$  .

1 Introduction

In the paper [14, we introduced the notion of $k$  -convexity, $k=1\dots n$  , for functions $u$  defined on domains $\Omega$  in Euclidean space, ${\mathbb{R}}^{n}$  . Namely, for $u\in {C}^{2}\left(\Omega \right)$  , we call $u$  $k$  -convex in $\Omega$  if
 $\begin{array}{c}{F}_{j}\left[u\right]:={F}_{j}\left({D}^{2}u\right):={S}_{j}\left(\lambda \right)\ge 0,\end{array}$ (1.1)
for $j=1\dots k$  , where $\lambda =\left({\lambda }_{1},\dots ,{\lambda }_{n}\right)$  are the eigenvalues of the Hessian matrix ${D}^{2}u$  of second derivatives of $u$  and ${S}_{j}$  denotes the $j$  th elementary symmetric function, that is
 $\begin{array}{c}{S}_{j}\left(\lambda \right)={\sum }_{{i}_{1}<\dots <{i}_{j}}{\lambda }_{{i}_{1}}\dots {\lambda }_{{i}_{j}},j=1\dots n.\end{array}$ (1.2)
When there is no confusion we use the same notation ${F}_{j}$  for both the operator and the function on ${\mathbb{R}}^{n}×{\mathbb{R}}^{n}$  . Equivalently, $u$  is $k$  -convex in $\Omega$  if $u$  is subharmonic with respect to the operator ${F}_{k}$  and this is the basis for our definition of $k$  -convexity for non-smooth functions in our sequel papers [15, [16, [17. The core result in our paper [15, is that the mapping $u↦{F}_{k}\left[u\right]$  is weakly continuous as a mapping from ${L}_{loc}^{1}\left(\Omega \right)$  to ${M}_{loc}\left(\Omega \right)$  , the space of locally finite measures in $\Omega$  , that is for any subdomain ${\Omega }^{\prime }\subset \subset \Omega$  , $\eta \in {C}_{0}^{0}\left({\Omega }^{\prime }\right)$  and positive constant $\epsilon$  , there exists a constant $\delta$  such that
 $\begin{array}{c}|{\int }_{\Omega }\eta \left({F}_{k}\left(u\right)-{F}_{k}\left(v\right)\right)|<\epsilon \end{array}$ (1.3)
whenever ${\int }_{{\Omega }^{\prime }}|u-v|<\delta ,{\int }_{{\Omega }^{\prime }}|u+v|<1,$  for arbitrary $k$  -convex $u$  and $v$  . This result enables us to define for any locally integrable $k$  -convex function $u$  , the Hessian measure, ${\mu }_{k}\left[u\right]$  , as an extension of ${F}_{k}\left[u\right]$  . In our first paper [14, we only proved the continuity of ${\mu }_{k}$  from ${C}^{0}\left(\Omega \right)$  to ${M}_{loc}\left(\Omega \right)$  but this was enough for the cases $k>\frac{n}{2}$  , which included the Monge-Ampére measure when $k=n$  . The weak continuity of the Monge-Ampére measure is a fundamental result of Aleksandrov (see eg.
[12).
In this paper we extend our results in [15to the case of non-commuting vector fields but only prove the corresponding weak continuity for the case $k=2$  . Our approach follows [15, with some help from [16and [17, and is inspired by the recent paper [8on the special case of the Heisenberg group ${\mathbb{H}}^{1}$  , by Gutiérrez and Montanari, where the more restrictive approach in [14was adequate. To formulate the main theorem, we let $X={X}_{1}\dots {X}_{m}$  denote a system of vector fields in ${\mathbb{R}}^{n}$  , that is first order differential operators of the form
 $\begin{array}{c}{X}_{i}={\sum }_{j=1}^{n}{b}^{ij}{D}_{j}\end{array}$ (1.4)
with coefficients ${b}^{ij}\in {C}^{\infty }\left(\overline{\Omega }\right)$  (although weaker regularity will suffice).
Then, for $k=1\dots m$  , we call a function $u\in {C}^{2}\left(\Omega \right)$  , $k$  -convex, with respect to $X$  if
 $\begin{array}{c}{F}_{j}\left[u\right]:={F}_{j}\left({X}_{s}^{2}u\right):={S}_{j}\left(\lambda \right)\ge 0,j=1\dots k\end{array}$ (1.5)
where now $\lambda =\left({\lambda }_{1}\dots {\lambda }_{m}\right)$  denote the eigenvalues of the symmetric Hessian,
 $\begin{array}{c}{X}_{s}^{2}u={\left[\frac{1}{2}\left({X}_{i}{X}_{j}+{X}_{j}{X}_{i}\right)u\right]}_{i,j=1\dots m}.\end{array}$ (1.6)
Our hypotheses on the vector fields ${X}_{1}\dots {X}_{m}$  are that:
• (i) they are anti-self adjoint, namely
${X}_{i}^{*}=-{X}_{i},i=1\dots m;$
• (ii) they satisfy the Hormander condition, namely the Lie algebra generated by them spans ${\mathbb{R}}^{n}$  and
• (iii) the second commutators formed from any two vector fields vanish.
These conditions will be automatically satisfied by the vector fields generating an homogeneous group of Heisenberg type. We can now state the main theorem.
Theorem 1.1. The mappings
 $\begin{array}{c}u↦{F}_{2}\left[u\right]+\alpha {\sum }_{i (1.7)
for $u$  2-convex in $\Omega$  , are weakly continuous from ${L}_{loc}^{1}\left(\Omega \right)$  to ${M}_{loc}\left(\Omega \right)$  , for any constant $\alpha$  .
As mentioned above, the special case of the Heisenberg group ${\mathbb{H}}^{1}$  , given by
 $\begin{array}{ccc}{X}_{1}u& =& {D}_{1}u-\frac{1}{2}{x}_{2}{D}_{3}u,\end{array}$
 $\begin{array}{ccc}{X}_{2}u& =& {D}_{2}u+\frac{1}{2}{x}_{1}{D}_{3}u,\end{array}$ (1.8)
 $\begin{array}{ccc}\left[{X}_{1},{X}_{2}\right]u& =& {D}_{3}u,\end{array}$
is proved in [8. Here ${L}_{loc}^{1}\left(\Omega \right)$  convergence is equivalent to local uniform convergence in $\Omega$  and the proof is much simpler.
Theorem  1.1 enables us to assign a Borel measure ${\mu }_{2}\left[u\right]$  to any ${L}_{loc}^{1}\left(\Omega \right)$  limit of smooth $2$  -convex functions, which extends ${F}_{2}\left[u\right]$  and is weakly continuous.
Letting ${\Phi }^{2}\left(\Omega \right)$  denote the space of such functions we also see that the commutators $\left[{X}_{i},{X}_{j}\right]u\in {L}_{loc}^{2}\left(\Omega \right)$  for $u\in {\Phi }^{2}\left(\Omega \right)$  .
This paper is arranged as follows. In the next section we generalize the basic divergence identity of Gutiérrez and Montanari [8, [9on Heisenberg groups to vector fields satisfying conditions  (i) and  (iii) . In Section 3, we employ our approach in [15to obtain integral estimates for the subelliptic gradient $Xu$  under conditiions  (i) and  (ii) , Theorem  3.1 . As we cannot extend all of our argument in [15, we have to rely strongly on the subelliptic potential estimates in [17. In Section 4, we carry out the arguments, again adapting [15to the non-commutative case, to conclude the local boundedness and weak continuity of the functionals in Theorem  1.1 , thereby completing the proof. Finally in Section 5, we extend our previous results to the classes ${\phi }^{k}\left(\Omega \right)$  of ${L}_{loc}^{1}$  limits of $k$  -convex functions. Additional remarks at the ends of Sections 4 and 5 treat the removal of condition  (iii) and more general definitions of $k$  -convexity.
We are grateful for useful comments and discussions with N. Chaudhuri, X-J Wang, A. Montanari and T. Nguyen.

2 Divergence structure and monotonicity

It is well known that if the vector fields ${X}_{1}\dots {X}_{m}$  commute, then
 $\begin{array}{c}{X}_{i}{F}_{k}^{ij}\left({X}^{2}u\right)=0,j=1\dots m,\end{array}$ (2.1)
where
 $\begin{array}{c}{X}^{2}u={X}_{s}^{2}u=\left[{X}_{i}{X}_{j}u{\right]}_{i,j=1\dots m}\end{array}$ (2.2)
and
 $\begin{array}{c}{F}_{k}^{ij}\left(r\right):=\frac{\partial }{\partial {r}_{ij}}{F}_{k}\left(r\right)\end{array}$ (2.3)
The identity ( 2.1 ), which means that the columns of the linearized coefficient matrix ( 2.3 ) are divergence free, was the basis for our approach in [14, [15.
Now suppose, more generally, that the second commutators formed from any two vector fields vanish, that is for any $i,j=1\dots m$  ,
 $\begin{array}{ccc}0& =& \left[{X}_{i},\left[{X}_{i},{X}_{j}\right]\right]\end{array}$
 $\begin{array}{ccc}& =& {X}_{i}\left[{X}_{i},{X}_{j}\right]-\left[{X}_{i},{X}_{j}\right]{X}_{i}\end{array}$
 $\begin{array}{ccc}& =& {X}_{i}{X}_{j}{X}_{j}-2{X}_{i}{X}_{j}{X}_{i}+{X}_{j}{X}_{i}{X}_{i}\end{array}$
 $\begin{array}{ccc}& =& {X}_{j}\left({X}_{i}{X}_{i}\right)+{X}_{i}\left({X}_{i}{X}_{j}-2{X}_{j}{X}_{i}\right).\end{array}$
Then, defining for any real matrix $r\in {\mathbb{R}}^{n}×{\mathbb{R}}^{n}$  ,
 $\begin{array}{ccc}{\mathfrak{F}}_{2}\left(r\right)& :=& \frac{1}{2}\left\{\left({r}_{ii}{\right)}^{2}-{r}_{ij}{r}_{ji}+\frac{1}{2}\left({r}_{ij}-{r}_{ji}{\right)}^{2}\right\}\end{array}$ (2.4)
 $\begin{array}{ccc}& =& {F}_{2}\left[r\right]+\frac{3}{4}{\sum }_{i
we have the identity
 $\begin{array}{c}{X}_{i}{\mathfrak{F}}_{2}^{ij}\left({X}^{2}u\right)=0,\end{array}$ (2.5)
where
 $\begin{array}{ccc}{\mathfrak{F}}_{2}^{ij}\left(r\right)& =& \frac{\partial {\mathfrak{F}}_{2}}{\partial {r}_{ij}}\left(r\right)\end{array}$ (2.6)
 $\begin{array}{ccc}& =& \text{(trace}r\text{)}{\delta }_{ij}+{r}_{ij}-2{r}_{ji},\end{array}$
which extends ( 2.1 ) in the case $k=2$  . The identity ( 2.5 ) was discovered by Gutiérrez and Montanari [8, [9for the Heisenberg groups ${\mathbb{H}}^{n}$  , (see also [5, [7).
From ( 2.5 ), we infer the monotonicity formula for the operator ${\mathfrak{F}}_{2}$  defined by
 $\begin{array}{c}{\mathfrak{F}}_{2}\left[u\right]={\mathfrak{F}}_{2}\left({X}^{2}u\right),\end{array}$ (2.7)
extending Lemma 2.1 in [14for $k=2$  .
Lemma 2.1. Let $u,v\in {C}^{2}\left(\Omega \right)\cap {C}^{0}\left(\overline{\Omega }\right)$  satisfy $u\le v$  in $\Omega$  , $u=v$  on $\partial \Omega$  with the operator ${F}_{2}$  degenerate elliptic with respect to their sum $u+v$  , that is
 $\begin{array}{ccc}{\mathfrak{F}}_{2}^{ij}\left({X}^{2}\left(u+v\right)\right){\xi }_{i}{\xi }_{j}& =& {F}_{2}^{ij}\left({X}^{2}\left(u+v\right)\right){\xi }_{i}{\xi }_{j}\end{array}$ (2.8)
 $\begin{array}{ccc}& \ge & 0\end{array}$
for all $\xi \in {\mathbb{R}}^{m}$  . Then, if the vector fields ${X}_{1}\dots {X}_{m}$  satisfy conditions  (i) and  (iii) , we have
 $\begin{array}{c}{\int }_{\Omega }{\mathfrak{F}}_{2}\left[v\right]\le {\int }_{\Omega }{\mathfrak{F}}_{2}\left[u\right]\end{array}$ (2.9)
• Proof. By integration by parts and the identity ( 2.5 ), we have, for $u,v\in {C}^{2}\left(\overline{\Omega }\right)$  ,  $\begin{array}{ccc}{\int }_{\Omega }\left({\mathfrak{F}}_{2}\left[u\right]-{\mathfrak{F}}_{2}\left[v\right]\right)& =& {\int }_{0}^{1}dt{\int }_{\Omega }{\mathfrak{F}}_{2}^{ij}\left[{X}^{2}\left(tu+\left(1-t\right)v\right)\right]{X}_{i}{X}_{j}\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& {\int }_{0}^{1}dt{\int }_{\partial \Omega }{\mathfrak{F}}_{2}^{ij}\left({X}_{i}.\gamma \right){X}_{j}\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& {\int }_{0}^{1}dt{\int }_{\partial \Omega }{\mathfrak{F}}_{2}^{ij}\left({X}_{i}.\gamma \right)\left({X}_{j}.\gamma \right)|D\left(u-v\right)|\end{array}$
 $\begin{array}{ccc}& \ge & 0\end{array}$
Here $\gamma$  denotes the outer unit normal to $\partial \Omega$  and
${X}_{i}.\gamma ={b}^{ij}{\gamma }_{j}.$
The general case $u,v\in {C}^{0}\left(\overline{\Omega }\right)\cap {C}^{2}\left(\Omega \right)$  follows by approximation.
For weak continuity with respect to ${C}^{0}\left(\Omega \right)$  and for groups of Heisenberg type we may proceed exactly as in [14. In the next section we present the basic gradient estimates for $k$  -convex functions needed to handle the general case.

In this section we provide the necessary gradient estimates for our proof of weak continuity. For these we do not have to restrict to the case $k=2$  and moreover we only need to assume the vector fields ${X}_{1}\dots {X}_{m}$  satisfy conditions  (i) and  (ii) . First we note that since $k$  -convexity implies 1-convexity, $k$  -convex functions $u$  are subharmonic with respect to the sub-Laplacian associated with ${X}_{1}\dots {X}_{n}$  , that is
 $\begin{array}{c}{\Delta }_{X}u:={X}_{i}{X}_{i}u\ge 0\end{array}$ (3.1)
in $\Omega$  . From ( 3.1 ) we infer immediately a bound from above, namely, for any ${\Omega }^{\prime }\subset \subset \Omega$  ,
 $\begin{array}{c}{sup}_{{\Omega }^{\prime }}u\le C{\int }_{\Omega }|u|\end{array}$ (3.2)
where the constant $C$  depends on ${X}_{1}\dots {X}_{m}$  and $dist\left({\Omega }^{\prime },\partial \Omega \right)$  . As we are only dealing with local estimates in this paper, we will always assume, without loss of generality, that $u\in {L}^{1}\left(\Omega \right)$  . Following [15, our treatment of gradient estimates depends on the relation between $k$  -convexity and the subelliptic p-Laplacian operators ${\Delta }_{p}$  defined by
 $\begin{array}{c}{\Delta }_{p}u={X}_{i}\left({|Xu|}^{p-2}{X}_{i}u\right)\end{array}$ (3.3)
for $p>1$  .
Lemma 3.1. Let $u$  be $k$  -convex in $\Omega$  . Then $u$  is subharmonic with respect to ${\Delta }_{p}$  for $p-1\le k\left(m-1\right)/\left(m-k\right)$  .
• Proof. Although this is just the special case $l=1$  in Lemma 4.2 of [15, we include it for completeness as it is simpler that the cases $l>1$  . We use the notation  $\begin{array}{c}{S}_{k,i}\left(\lambda \right)={S}_{k}\left(\lambda \right){|}_{{\lambda }_{i}=0}\end{array}$ (3.4)
so that  $\begin{array}{c}{S}_{j,i}\left(\lambda \right)\ge 0\end{array}$ (3.5)
for all $j\le k-1$  , if ${S}_{j}\left(\lambda \right)\ge 0$  , for all $j=1,\dots k$  . It follows then that  $\begin{array}{c}0\le {S}_{k}\left(\lambda \right)={S}_{k,i}\left(\lambda \right)+{S}_{k-1,i}\left(\lambda \right){\lambda }_{i},i=1,\dots m,\end{array}$ (3.6)
whence  $\begin{array}{ccc}-{\lambda }_{i}& \le & \frac{{S}_{k,i}}{{S}_{k-1,i}}\left(\lambda \right)\end{array}$ (3.7)
 $\begin{array}{ccc}& \le & \frac{\left(m-k\right)}{k\left(m-1\right)}{S}_{1,i}\left(\lambda \right)\end{array}$
by MacLaurins' inequality for ratios of elementary symmetric functions.
Consequently, if $u$  is $k$  -convex,  $\begin{array}{ccc}{\Delta }_{p}u& =& {X}_{i}\left({|Xu|}^{p-2}{X}_{i}u\right)\end{array}$ (3.8)
 $\begin{array}{ccc}& =& {|Xu|}^{p-2}\left\{\Delta u+\left(p-2\right)\frac{{X}_{i}u{X}_{j}u}{{|Xu|}^{2}}{X}_{i}{X}_{j}u\right\}\end{array}$
 $\begin{array}{ccc}& \ge & 0,\end{array}$
for $p-1\le k\left(m-1\right)/\left(m-k\right)$  , by taking ${\lambda }_{1}\dots {\lambda }_{m}$  to be the eigenvalues of ${X}_{s}^{2}u$  .
Note that Lemma  3.1 also includes the case $k=m,p=\infty$  , when ${\Delta }_{p}$  is the subelliptic $\infty$  -Laplacian
 $\begin{array}{c}{\Delta }_{\infty }u={X}_{i}u{X}_{j}u{X}_{ij}u\end{array}$ (3.9)
Our gradient estimates now follow immediately from [17but to express them we need the concept of homogeneous dimension. For our purposes here, we define the C-C (Carnot-Caratheodory) metric induced from the vector fields ${X}_{1}\dots {X}_{m}$  by $d\left(x,y\right)=inf\left\{T>0|\exists \text{a sub-unitary}\gamma :\left[0,T\right]\to {\mathbb{R}}^{n}$
 $\begin{array}{c}\text{with}\gamma \left(0\right)=x,\gamma \left(T\right)=y\right\},\end{array}$ (3.10)
where a piecewise ${C}^{1}$  curve $\gamma :\left[0,T\right]\to {\mathbb{R}}^{n}$  is said to be sub-unitary, with respect to ${X}_{1}\dots {X}_{m}$  , if for every $\xi \in {\mathbb{R}}^{n}$  and $t\in \left(0,T\right)$  ,
 $\begin{array}{c}{|{\gamma }^{\prime }\left(t\right).\xi |}^{2}\le {\sum }_{i=1}^{m}\left({X}_{i}\left(\gamma \left(t\right)\right).\xi {\right)}^{2}.\end{array}$ (3.11)
Let ${B}_{R}\left(x\right)$  denote the C-C ball $\left\{y\in {\mathbb{R}}^{n}|d\left(x,y\right)  , and let $\Omega$  be a bounded domain in ${\mathbb{R}}^{n}$  . The the fundamental result of Nagel, Stein and Wainger [11asserts that there exist positive constants $C,{R}_{0}$  and positive integer $Q$  , depending on $X$  and $\Omega$  such that
 $\begin{array}{c}|{B}_{tR}\left(x\right)|\ge C{t}^{Q}|{B}_{R}\left(x\right)|\end{array}$ (3.12)
for any $x\in \Omega$  , $t\in \left(0,1\right)$  and $R<{R}_{0}$  , where $||$  denotes the Lebesgue volume.
The number $Q$  ( $\ge n$  ), is chosen as the least integer for which ( 3.12 ) holds and is called the homogeneous dimension of $X$  in $\Omega$  . For ( 3.12 ) we only need the Hormander condition  (ii) and we could replace it more generally in this paper by simply the validity of ( 3.12 ).
Theorem 3.1. For any $k$  -convex function $u$  in $\Omega$  , and subdomain ${\Omega }^{\prime }\subset \subset \Omega$  , we have the estimates
 $\begin{array}{c}||Xu|{|}_{{L}^{q}\left({\Omega }^{\prime }\right)}\le C\left({\int }_{\Omega }|u|\right)\end{array}$ (3.13)
 $\begin{array}{c}{\int }_{{\Omega }^{\prime }}{|Xu|}^{r}{\Delta }_{X}u\le C{\left({\int }_{\Omega }|u|\right)}^{1+r}\end{array}$ (3.14)
for $1\le q  , $0\le r  where $C$  depends on $\Omega$  , ${\Omega }^{\prime }$  , ${X}_{1}\dots {X}_{m}$  and $q$  or $r$  as appropriate.
• Proof. The estimate ( 3.13 ) follows from Lemma  3.1 and [17. For ( 3.14 ), we have from Lemma  3.1 ,  $\begin{array}{c}{|Xu|}^{r}{\Delta }_{X}u\le \frac{m\left(k-1\right)}{m\left(k-1\right)-r\left(m-k\right)}{\Delta }_{p}u\end{array}$ (3.15)
and since  $\begin{array}{c}{\int }_{\Omega }\eta {\Delta }_{p}u\le {\int }_{\Omega }|X\eta |{|Xu|}^{p-1},\end{array}$ (3.16)
for any $\eta >0,\in {C}_{0}^{1}\left(\Omega \right)$  , we infer ( 3.14 ) from ( 3.13 ).
By using the subelliptic Sobolev inequality [2, we obtain corresponding ${L}^{p}$  estimates, namely
 $\begin{array}{c}||u|{|}_{{L}^{p}\left({\Omega }^{\prime }\right)}\le C||u|{|}_{{L}^{1}\left(\Omega \right)}\end{array}$ (3.17)
where
$1\le p<\frac{Qk\left(m-1\right)}{\left(Q-1\right)m-\left(Q+m-2\right)k}$
for $\left(Q-1\right)m\ge \left(Q+m-2\right)k$  and $p=\infty$  if $\left(Q-1\right)m<\left(Q+m-2\right)k$  . For this last case we have a Hölder estimate [3, [17,
 $\begin{array}{c}{sup}_{{\Omega }^{\prime }}\frac{|u\left(x\right)-u\left(y\right)|}{{|d\left(x,y\right)|}^{\alpha }}\le C||u|{|}_{{L}^{1}\left(\Omega \right)}\end{array}$ (3.18)
where
$\alpha =\frac{k\left(Q+m-2\right)-m\left(Q-1\right)}{k\left(m-1\right)},$
if $k  ( $\alpha <$  1 if $k=m$  ).

4 Weak Continuity

In this section, we complete the proof of Theorem  1.1 . First we prove a local bound for ${\mathfrak{F}}_{2}$  . For convenience we use the notation
 $\begin{array}{c}{\mathfrak{E}}_{2}\left[u\right]={\sum }_{i (4.1)
so that
 $\begin{array}{c}{\mathfrak{F}}_{2}\left[u\right]={F}_{2}\left[u\right]+\frac{3}{4}{\mathfrak{E}}_{2}\left[u\right]\end{array}$ (4.2)
Lemma 4.1. Let $u\in {C}^{2}\left(\Omega \right)$  be 2-convex in $\Omega$  with respect to ${X}_{1},\dots {X}_{m}$  , satisfying hypotheses  (i) to  (iii) . Then, for any subdomain ${\Omega }^{\prime }\subset \subset \Omega$  , we have
 $\begin{array}{c}{\int }_{{\Omega }^{\prime }}{\mathfrak{F}}_{2}\left[u\right]\le C{\left({\int }_{\Omega }|u|\right)}^{2},\end{array}$ (4.3)
where $C$  depends on $dist\left({\Omega }^{\prime },\partial \Omega \right)$  , ${X}_{1},\dots {X}_{m}$  .
In particular, Lemma  4.1 provides a local ${L}^{2}$  estimate for the commutators $\left[{X}_{i},{X}_{j}\right]u$  .
• Proof. Letting $\eta >0\in {C}_{0}^{1}\left(\Omega \right)$  be a cut-off function, we have  $\begin{array}{ccc}{\int }_{\Omega }{\eta }^{2}{\mathfrak{F}}_{2}\left[u\right]& =& \frac{1}{2}{\int }_{\Omega }{\eta }^{2}{\mathfrak{F}}_{2}^{ij}{X}_{i}{X}_{j}u\end{array}$ (4.4)
 $\begin{array}{ccc}& =& -\frac{1}{2}{\int }_{\Omega }{\mathfrak{F}}_{2}^{ij}{X}_{i}{\eta }^{2}{X}_{j}u\end{array}$
 $\begin{array}{ccc}& =& -\frac{1}{2}{\int }_{\Omega }\left\{{F}_{2}^{ij}{X}_{i}{\eta }^{2}{X}_{j}u+\frac{3}{4}\left[{X}_{i},{X}_{j}\right]u{X}_{i}{\eta }^{2}{X}_{j}u\right\}\end{array}$
 $\begin{array}{ccc}& \le & C{\int }_{\Omega }\left\{|X{\eta }^{2}||Xu|{\Delta }_{X}u+\eta {\mathfrak{E}}_{2}^{\frac{1}{2}}\left[u\right]|X\eta ||Xu|\right\}.\end{array}$
Consequently,  $\begin{array}{c}{\int }_{\Omega }{\eta }^{2}{\mathfrak{F}}_{2}\left[u\right]\le C{\int }_{\Omega }\left\{|X{\eta }^{2}||Xu|{\Delta }_{X}u+{|X\eta |}^{2}{|Xu|}^{2}\right\}\end{array}$ (4.5)
and ( 4.3 ) follows from Theorem  3.1 .
• Proof of Theorem  1.1 . Letting $u$  and $v$  be 2-convex in $\Omega$  , ${u}_{t}=tu+\left(1-t\right)v,0\le t\le 1$  and $\eta \ge 0,\in {C}_{0}^{2}\left(\Omega \right)$  , we now have  $\begin{array}{ccc}& & {\int }_{\Omega }\eta \left({\mathfrak{F}}_{2}\left[u\right]-{\mathfrak{F}}_{2}\left[v\right]\right)\end{array}$ (4.6)
 $\begin{array}{ccc}& =& {\int }_{0}^{1}dt{\int }_{\Omega }\eta {\mathfrak{F}}_{2}^{ij}\left[{u}_{t}\right]{X}_{i}{X}_{j}\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}{\int }_{\Omega }\eta {\mathfrak{F}}_{2}^{ij}\left[u+v\right]{X}_{i}{X}_{j}\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& -\frac{1}{2}{\int }_{\Omega }{\mathfrak{F}}_{2}^{ij}{X}_{i}\eta {X}_{j}\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& -\frac{1}{2}{\int }_{\Omega }\left\{{\mathfrak{F}}_{2}^{ji}+\frac{3}{2}\left[{X}_{i},{X}_{j}\right]\left(u+v\right)\right\}{X}_{i}\eta {X}_{j}\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}{\int }_{\Omega }\left\{{\mathfrak{F}}_{2}^{ij}{X}_{i}{X}_{j}\eta \left(u-v\right)-\frac{3}{2}\left[{X}_{i},{X}_{j}\right]\left(u+v\right){X}_{i}\eta {X}_{j}\left(u-v\right)\right\}\end{array}$
 $\begin{array}{ccc}& =& \frac{1}{2}{\int }_{\Omega }{F}_{2}^{ij}\left[u+v\right]\left({X}_{i}{X}_{j}\eta \right)\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& & +\frac{3}{8}{\int }_{\Omega }\left[{X}_{i},{X}_{j}\right]\left(u+v\right)\left\{\left({X}_{i}{X}_{j}\eta \right)\left(u-v\right)-2{X}_{i}\eta {X}_{j}\left(u-v\right)\right\}\end{array}$
 $\begin{array}{ccc}& :=& \frac{1}{2}{I}_{1}+\frac{3}{8}{I}_{2}\end{array}$
The estimation of the integral ${I}_{1}$  is similar to the corresponding term in the commuting case [15. Namely  $\begin{array}{ccc}|{I}_{1}|& =& |{\int }_{\Omega }{F}_{2}^{ij}\left[u+v\right]\left({X}_{i}{X}_{j}\eta \right)\left(u-v\right)|\end{array}$ (4.7)
 $\begin{array}{ccc}& \le & {\int }_{\Omega }{\Delta }_{X}\left(u+v\right)|{X}^{2}\eta ||u-v|.\end{array}$
Denoting  $\begin{array}{c}\delta ={\int }_{\Omega }|u-v|,K={\int }_{\Omega }|u+v|,\end{array}$ (4.8)
we have, for any $\epsilon >0$  ,
$|u-v|<\epsilon$
except on a set ${A}_{\epsilon }$  of measure $|{A}_{\epsilon }|\le \delta /\epsilon$  . We then estimate for a further cut-off function $\stackrel{~}{\eta }\in {C}_{0}^{1}\left(\Omega \right)$  , $0\le \stackrel{~}{\eta }\le 1$  ,  $\begin{array}{ccc}& & {\int }_{\Omega }\stackrel{~}{\eta }{\Delta }_{X}\left(u+v\right)\left(u-v{\right)}^{+}\end{array}$ (4.9)
 $\begin{array}{ccc}& \le & \epsilon {\int }_{\Omega }\stackrel{~}{\eta }{\Delta }_{X}\left(u+v\right)+{\int }_{\Omega }\stackrel{~}{\eta }\left(u-v-\epsilon {\right)}^{+}{\Delta }_{X}\left(u+v\right)\end{array}$
 $\begin{array}{ccc}& =& \epsilon {\int }_{\Omega }\stackrel{~}{\eta }{\Delta }_{X}\left(u+v\right)-{\int }_{{A}_{\epsilon }}{X}_{i}\left(u+v\right){X}_{i}\left\{\stackrel{~}{\eta }\left(u-v-\epsilon {\right)}^{+}\right\}\end{array}$
 $\begin{array}{ccc}& \le & C\left(\epsilon +{|{A}_{\epsilon }|}^{1-\frac{2}{q}}\right)\end{array}$
 $\begin{array}{ccc}& \le & C\left\{\epsilon +{\left(\frac{\delta }{\epsilon }\right)}^{1-\frac{2}{q}}\right\}\end{array}$
by Theorem  3.1 , where $q$  is chosen so that
$2
and $C$  depends on $\stackrel{~}{\eta }$  , $K$  , $\delta$  , ${X}_{1},\dots {X}_{m}$  . It then follows that  $\begin{array}{c}|{I}_{1}|\le C\left\{\epsilon +{\left(\frac{\delta }{\epsilon }\right)}^{1-\frac{2}{q}}\right\}\end{array}$ (4.10)
where $C$  depends on $\eta$  , $K$  , $\delta$  , ${X}_{1},\dots {X}_{m}$  . To estimate ${I}_{2}$  , we first use the bound for ${\mathfrak{E}}_{2}$  in Lemma  4.1 , to obtain  $\begin{array}{c}{I}_{2}\le C{\left\{{\int }_{supp\eta }{|X\left(u-v\right)|}^{2}+{|u-v|}^{2}\right\}}^{\frac{1}{2}},\end{array}$ (4.11)
where $C$  depends on $\eta$  , $K$  , $\delta$  , ${X}_{1},\dots {X}_{m}$  . The first part of the above integral may be estimated similarly to ${I}_{1}$  , since  $\begin{array}{ccc}{\int }_{\Omega }\stackrel{~}{\eta }{|X\left(u-v\right)|}^{2}& =& -{\int }_{\Omega }{X}_{i}\left(\stackrel{~}{\eta }{X}_{i}\left(u-v\right)\right)\left(u-v\right)\end{array}$ (4.12)
 $\begin{array}{ccc}& \le & {\int }_{\Omega }\left\{|X\stackrel{~}{\eta }||X\left(u-v\right)|+\stackrel{~}{\eta }{\Delta }_{X}\left(u-v\right)\right\}|u-v|\end{array}$
 $\begin{array}{ccc}& \le & {\int }_{\Omega }\stackrel{~}{\eta }{\Delta }_{X}\left(u-v\right)|u-v|+C{\left\{{\int }_{supp\stackrel{~}{\eta }}{|u-v|}^{2}\right\}}^{\frac{1}{2}},\end{array}$
while the second parts of ( 4.11 ) and ( 4.12 ) are handled readily by the estimate ( 3.17 ). It follows that  $\begin{array}{c}{I}_{2}\le C{\left\{\epsilon +{\left(\frac{\delta }{\epsilon }\right)}^{1-\frac{2}{q}}\right\}}^{\frac{1}{2}}.\end{array}$ (4.13)
With appropriate choice of $\delta$  , we then conclude Theorem  1.1 from the estimates ( 4.10 ) and ( 4.13 ).
Remark. By inspection of the above proofs, we see that condition  (iii) may be weakened to only requiring that the vector fields ${Y}_{j}:={\sum }_{i=1}^{m}\left[{X}_{i},\left[{X}_{i},{X}_{j}\right]\right],j=1\dots m,$  lie in the span of ${X}_{i},\left[{X}_{i},{X}_{j}\right]$  , $i,j=1\dots m.$  Without condition  (iii) additional terms
 $\begin{array}{c}-\frac{1}{2}{\int }_{\Omega }{\eta }^{2}Xu.Yu\end{array}$ (4.14)
and
 $\begin{array}{c}-\frac{1}{2}{\int }_{\Omega }\eta X\left(u-v\right).Y\left(u+v\right)+\frac{1}{2}{\int }_{\Omega }\left(u-v\right)X\eta .Y\left(u+v\right)\end{array}$ (4.15)
will arise in the right hand sides of ( 4.4 ) and ( 4.6 ) respectively and these are then automatically controlled. More general hypotheses are clearly possible.
Note also that if $Y$  commutes with $X$  , then the quantity
 $\begin{array}{c}{\mathfrak{F}}_{2}^{*}\left[u\right]:={\mathfrak{F}}_{2}\left[u\right]+\frac{1}{2}Xu.Yu\end{array}$ (4.16)
also satisfies the monotonicity property ( 2.9 ). For the Engel group, this was also observed in [7.
More generally if we integrate by parts in ( 4.15 ), we find
 $\begin{array}{ccc}& & {\int }_{\Omega }\eta X\left(u-v\right).Y\left(u+v\right)\end{array}$ (4.17)
 $\begin{array}{ccc}& =& -{\int }_{\Omega }\eta \left(u+v\right)Y.X\left(u-v\right)+\left(u+v\right)Y\eta .X\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& {\int }_{\Omega }\eta \left\{X\left(u+v\right).Y\left(u-v\right)+\left(u+v\right)Z\left(u-v\right)\right\}\end{array}$
 $\begin{array}{ccc}& & -{\int }_{\Omega }\left(u+v\right)\left\{Y\eta .X\left(u-v\right)-X\eta .Y\left(u-v\right)\right\}\end{array}$
where $Z$  is the vector field given by
 $\begin{array}{c}Z=X.Y-Y.X.\end{array}$ (4.18)
Inserting ( 4.16 ) into ( 4.6 ) and ( 4.15 ) we may infer weak continuity results, in the absence of condition  (iii) , with respect to stronger topologies. For example ${F}_{2}$  will be weakly continuous with respect to ${L}_{loc}^{1}\left(\Omega \right)$  on 2-convex functions, with uniformly bounded $X_$  gradients and $Y,Z_$  gradients in ${L}_{loc}^{\infty }\left(\Omega \right)$  and ${L}_{loc}^{1}\left(\Omega \right)$  respectively.
Finally without hypothesis  (iii) in Lemma 2.1, we obtain, from the proof,
 $\begin{array}{ccc}{\int }_{\Omega }\left({\mathfrak{F}}_{2}\left[u\right]-{\mathfrak{F}}_{2}\left[v\right]\right)& \ge & -\frac{1}{2}{\int }_{\Omega }Y\left(u+v\right).X\left(u-v\right)\end{array}$
 $\begin{array}{ccc}& =& -\frac{1}{4}{\int }_{\Omega }\left\{Y\left(u+v\right).X\left(u-v\right)+Y\left(u-v\right)X\left(u+v\right)\end{array}$
 $\begin{array}{ccc}& & -Z\left(u+v\right)\left(u-v\right)\right\}\end{array}$
so that
 $\begin{array}{ccc}{\int }_{\Omega }\left({\mathfrak{F}}_{2}^{*}\left[u\right]-{\mathfrak{F}}_{2}^{*}\left[v\right]\right)& \ge & \frac{1}{4}{\int }_{\Omega }\left(u-v\right)Z\left(u+v\right)\end{array}$ (4.19)
 $\begin{array}{ccc}& \ge & 0\end{array}$
if $Z\left(u+v\right)\le 0$  . In a similar fashon, we may combine ( 4.15 ) and ( 4.17 ) in ( 4.6 ) to conclude,
 $\begin{array}{ccc}{\int }_{\Omega }\eta \left({\mathfrak{F}}_{2}^{*}\left[u\right]-{\mathfrak{F}}_{2}^{*}\left[v\right]\right)& =& \frac{1}{2}{I}_{1}+\frac{3}{8}{I}_{2}+\frac{3}{4}{\int }_{\Omega }\left(u-v\right)X\eta .Y\left(u+v\right)\end{array}$ (4.20)
 $\begin{array}{ccc}& +& \frac{1}{4}{\int }_{\Omega }\left({u}^{2}-{v}^{2}\right)X.Y\eta +\frac{1}{4}{\int }_{\Omega }\eta \left(u-v\right)Z\left(u+v\right)\end{array}$

5 General $k$  -convex functions

For our purposes here, we define a function $u\in {L}_{loc}^{1}\left(\Omega \right)$  to be $k$  -convex in $\Omega$  , with respect to the system of vector fields $X=\left({X}_{1}\dots {X}_{m}\right)$  if, for any ${\Omega }^{\prime }\subset \subset \Omega$  , there exists a sequence of $k$  -convex functions in ${C}^{2}\left({\Omega }^{\prime }\right)$  converging to $u$  in ${L}_{loc}^{1}\left({\Omega }^{\prime }\right)$  . We designate the general class of $k$  -convex functions in $\Omega$  by ${\phi }^{k}\left(\Omega \right)$  . The estimates of the preceeding sections then extend as regularity properties and imbeddings of ${\phi }^{k}\left(\Omega \right)$  . In particular if the system $X$  satisfies conditions  (i) and  (ii) , we have from Theorem  3.1 that the distributional derivatives $Xu\in {L}_{loc}^{q}\left(\Omega \right)$  for
 $\begin{array}{c}1\le q<\frac{Qk\left(m-1\right)}{\left(Q-1\right)\left(m-k\right)}\end{array}$ (5.1)
and that ${\phi }^{k}\left(\Omega \right)$  imbeds continuously in the Sobolev space,
 $\begin{array}{c}{S}_{loc}^{1,q}\left(\Omega \right)=\left\{u\in {L}_{loc}^{1}\left(\Omega \right)|Xu\in {L}_{loc}^{q}\left(\Omega \right)\right\},\end{array}$ (5.2)
as defined, for example, in [2, [17. Moreover for
 $\begin{array}{c}k>\frac{\left(Q-1\right)m}{\left(Q+m-2\right)},\end{array}$ (5.3)
we see from ( 3.18 ) that ${\phi }^{k}\left(\Omega \right)$  imbeds continuously in the Hölder space ${C}^{0,\alpha }\left(\Omega \right)$  where
 $\begin{array}{c}\alpha =\frac{k\left(Q+m-2\right)-m\left(Q-1\right)}{k\left(m-1\right)},\end{array}$ (5.4)
if $k  ( $\alpha <1$  if $k=m$  ). Note that for smaller $\alpha$  , this follows from ( 3.13 ) but we need the full strength of the capacity estimates in [3to get the sharp exponent ( 5.4 ). It follows from [10that we can also take $\alpha =1$  when $k=m$  and $X$  generates the Lie algebra of a Carnot group.
For $k\ge 2$  , the symmetric Hessian ${X}_{s}^{2}u$  consists of signed Radon measures.
This follows exactly as in the Euclidean case [4, [15but it may also be observed directly from the degenerate ellipticity of ${F}_{2}$  ,
 $\begin{array}{c}{\Delta }_{X}uI-{X}_{s}^{2}u\ge 0.\end{array}$ (5.5)
By taking limits, we see that the above quantities are measures if $u\in {\phi }^{k}\left(\Omega \right)$  .
If we also assume $X$  satisfies  (iii) , we obtain from Lemma  4.1 , that the commutators $\left[{X}_{i},{X}_{j}\right]u\in {L}_{loc}^{2}\left(\Omega \right)$  for $u\in {\phi }^{k}\left(\Omega \right)$  , $k\ge 2$  . Hence the full Hessian ${X}^{2}u$  comprises Radon measures. From Theorem  3.1 , we also infer that ${\phi }^{k}\left(\Omega \right)$  , $k\ge 2$  , imbeds continuously into the Euclidean Sobolev space ${W}_{loc}^{1,2}\left(\Omega \right)$  if condition  (iii) is strengthened to all second order commutators vanishing. Also if $X$  generates the Lie algebra of a Carnot group as in [1, then from the Hölder estimate ( 3.18 ) and the weak differentiability result in [1we conclude that functions in ${\phi }^{k}\left(\Omega \right)$  will be twice differentiable almost everywhere in $\Omega$  , with respect to the system $X$  , if $k$  satisfies ( 5.3 ). This extends the corresponding result for the convex case $k=m$  in [5, [9, [10and the Euclidean case $k>n/2$  in [4.
Finally we may define the Hessian measure ${\mu }_{2}\left[u\right]$  , with respect to $X$  , for any function $u\in {\phi }^{2}\left(\Omega \right)$  by
 $\begin{array}{c}{\int }_{\Omega }\eta d{\mu }_{2}\left[u\right]={lim}_{m\to \infty }{\int }_{\Omega }\eta d{\mu }_{2}\left[{u}_{m}\right]\end{array}$ (5.6)
where $\eta \in {C}_{0}^{0}\left(\Omega \right)$  , $\text{supp}\eta \subset {\Omega }^{\prime }\subset \subset \Omega$  and $\left\{{u}_{m}\right\}\subset {C}^{2}\left({\Omega }^{\prime }\right)$  is a sequence of $k$  -convex functions converging to $u$  in ${L}_{loc}^{1}\left({\Omega }^{\prime }\right)$  . By Theorem  1.1  ${\mu }_{2}$  is well defined and weakly continuous with respect to convergence in ${L}_{loc}^{1}\left(\Omega \right)$  , that is if $\left\{{u}_{m}\right\}\subset {\phi }^{2}\left(\Omega \right)$  , converges to $u\in {\phi }^{2}\left(\Omega \right)$  in ${L}_{loc}^{1}\left(\Omega \right)$  , the corresponding sequence of measures $\left\{{\mu }_{2}\left[{u}_{m}\right]\right\}$  converges weakly to ${\mu }_{2}\left[u\right]$  . Note that the case $\alpha =\frac{1}{2}$  in Theorem  1.1 also shows that the sum of the principal $2×2$  minors of the full Hessian ${X}^{2}u$  also extends as a weakly continuous measure on ${\phi }^{2}\left(\Omega \right)$  .
From Lemma  2.1 , we also conclude a more general monotonicity property, namely that if $u,v\in {\phi }^{2}\left(\Omega \right)$  satisfy $u\le v$  in $\Omega$  , $u=v$  continuously on $\partial \Omega$  , then
 $\begin{array}{c}{\mu }_{2}\left[v\right]\left(\Omega \right)+\frac{3}{4}{\mathfrak{E}}_{2}\left[v\right]\le {\mu }_{2}\left[u\right]\left(\Omega \right)+\frac{3}{4}{\mathfrak{E}}_{2}\left[u\right].\end{array}$ (5.7)
General subharmonic functions.
More generally we may define subharmonic functions along the lines of [15, [16, [17. In particular we define an upper-semicontinuous function $u:\Omega ↦\left[-\infty ,\infty \right)$  to be subharmonic with respect to the operator ${F}_{k}$  if $u$  satisfies ${F}_{k}\left[u\right]\ge 0$  in the viscosity sense, that is for any quadratic polynomial $q$  for which the difference $u-q$  has a finite local maximum at a point $y\in \Omega$  , we have ${F}_{k}\left[q\right]\left(y\right)\ge 0$  . For smooth vector fields and $k=1$  , this is also equivalent to our definition in [17, which corresponds to the traditional definition of subharmonicity. A $k$  -convex function, as defined above by approximation, will be equivalent to a subharmonic function and moreover the estimates of Section 3 extend to the class of proper subharmonic functions. For Carnot groups of step 2, it follows from [10that proper subharmonic functions will also be $k$  -convex but we would expect this characterization to hold more generally.
References

1. L. Ambrosio and V. Magnani. Weak differentiability of BV functions on stratified groups. Math. Zeit., (245):123 – 153, 2003.
2. L. Capogna, D. Danielli, and N. Garofalo. An embedding theorem and the Harnack inequality for nonlinear subelliptic equations. Comm. Part. Diff. Eqns, (18):1765 – 1794, 1993.
3. L. Capogna, D. Danielli, and N. Garofalo. Capacitary estimates and the local behaviour of solutions of nonlinear subelliptic equations. Amer. J. Math, (118):1153 – 1196, 1997.
4. N. Chaudhuri and N.S. Trudinger. An Alexsandrov type theorem for k-convex functions. Bull. Aust. Math. Soc., 71:305–314, 2005.
5. D. Danielli, N. Garofalo, D.M. Nhieu, and F. Tournier. The theorem of Busemann-Feller-Alexandrov in Carnot groups. Comm. Anal. Geom, 12(4):853–886, 2004.
6. G.B. Folland and E.M. Stein. Hardy spaces on homogeneous groups. Princeton Univ. Press, 1982.
7. N. Garofalo and F. Tournier. New properties of convex functions in the Heisenberg group. Trans. Amer. Math. Soc., to appear.
8. C.E. Gutiérrez and A. Montanari. Maximum and comparison principles for convex functions on the Heisenberg group. Comm. Part. Diff. Eqns., 29:1305–1334, 2004.
9. C.E. Gutiérrez and A. Montanari. On the second order derivatives of convex functions on the Heisenberg group. Ann. Scuola Norm. Sup. Pisa Cl. Sci (5), (3):349–366, 2004.
10. V. Magnani. Lipschitz continuity, Aleksandrov theorem and characterizations for H-convex functions. preprint.
11. A. Nagel, E.M. Stein, and S. Wainger. Balls and metrics defined by vector fields. Acta. Math., (155):103–147, 1985.
12. A.V. Pogorelov. Monge-Ampére equations of elliptic type. Noordhoff, Groningen, 1964.
13. E.M. Stein. Harmonic Analysis. Princeton Univ Press, 1993.
14. N.S. Trudinger and X.J. Wang. Hessian Measures I. Topol. Methods Nonlinear Anal., (10):225–239, 1997.
15. N.S. Trudinger and X.J. Wang. Hessian Measures II. Ann. Math, (150):579–604, 1999.
16. N.S. Trudinger and X.J. Wang. Hessian Measures III. J. Funct. Anal., (193):1–23, 2002.
17. N.S. Trudinger and X.J. Wang. On the weak continuity of elliptic operators and applications to potential theory. Amer. J. Math, (124):369 – 410, 2002.

Centre for Mathematics and its Applications, Australian National University, Canberra, ACT 0200, Australia. Email: neil.trudinger@anu.edu.au