Lecture Notes On Elliptic PDE Luigi Ambrosio
Lecture Notes On Elliptic PDE Luigi Ambrosio
⇤
Luigi Ambrosio
Contents
1 Some basic facts concerning Sobolev spaces 3
4 Regularity Theory 33
4.1 Nirenberg method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
10 Schauder theory 61
11 Regularity in Lp spaces 65
⇤
PhD course given in 2009-2010 and then in 2012-2013, 2014-2015, lectures typed by A.Carlotto and
A.Massaccesi
1
12 Some classical interpolation theorems 68
14 Calderón-Zygmund decomposition 72
Preface
Prerequisites: basic knowledge of Functional Analysis and Measure Theory, preferably
also a basic knowledge of Sobolev spaces of functions of one independent variable.
2
Br (x) Ball with center x and radius r (also Br = Br (0), B = B1 )
A⇢B Inclusion in the weak sense
AbB A ⇢ B (typically used for pairs of open sets)
Ln Lebesgue measure in Rn
C k (⌦) Functions continuously k-di↵erentiable in ⌦
Lp (⌦) Lebesgue Lp space
@u
@i u, @xi u, ri u, @xi
i-th partial derivative (weak or classical)
ru
R Gradient of u R
⌦
f dµ Mean integral value, namely ⌦ f dµ/µ(⌦)
We define the space H 1,p (⌦) to be the completion with respect to the W 1,p norm of C 1 (⌦).
For unbounded domains, including the whole space Rn , the definition is similar and
based on the completion of
3
Whenever such g1 , . . . , gn exist, we say that u is di↵erentiable in weak sense and we write
gi = @i u and
ru = @1 u, . . . , @n u).
For 1 p 1 we define the space W 1,p (⌦) as the subset of Lp (⌦) whose elements u are
weakly di↵erentiable with corresponding derivatives @i u also belonging to Lp (⌦).
It is clear that if gi exists, it must be uniquely determined up to Lebesgue negligible
sets, since h 2 L1loc (⌦) and
Z
h' dx = 0 8' 2 Cc1 (⌦)
⌦
implies h = 0. This implication can be easily proved by approximation, showing that the
property above is stable under convolution, namely h" = h ⇤ ⇢" satisfies
Z Z
h" ' dx = h' ⇤ ⇢" = 0 8' 2 Cc1 (⌦" ),
⌦" ⌦
⇢" (x) = " n ⇢(x/") with ⇢ smooth, even and compactly supported in the unit ball and we
used the simmetry property (a consequence of Fubini’s theorem).
Z Z
(a ⇤ ⇢" )b dx = a(b ⇤ ⇢" ) dx. (1.4)
Obviously, classical derivatives are weak derivatives and the notation @i u (or, equivalently,
@u
@xi u, ri u or even @x i
) is justified.
Another classical way to relate weak and strong derivatives is via convolution: namely
if u has weak i-th derivative equal to g, then
Knowing the identity (1.5) for smooth functions, its validity can be easily extended con-
sidering both sides as weak derivatives and using (1.4):
Z Z Z Z Z
(u⇤⇢" )@i ' dx = u(@i ')⇤⇢" dx = u@i ('⇤⇢" ) dx = g'⇤⇢" dx = g ⇤⇢" ' dx
⌦ ⌦ ⌦ ⌦ ⌦
for all ' 2 Cc1 (⌦" ). Now, the smoothness of u ⇤ ⇢" tells us that the derivative in the left
hand side of (1.5) is (equivalent to) a classical one.
Another consequence of (1.5) is:
4
Theorem 1.3 (Constancy theorem). If u 2 L1loc (⌦) satisfies ru = 0 in the weak sense,
then for any ball B ⇢ ⌦ there exists a constant c 2 R such that u = c L n -a.e. in B. In
particular, if ⌦ is connected, u = c L n -a.e. in ⌦ for some c 2 R.
Proof. Again we argue by approximation, using the fact that (1.5) ensures that the func-
tion u ⇤ ⇢" are locally constant in ⌦" and taking the L1loc limit as " ! 0.
Notice also that Definition 1.2 covers the case p = 1, while it is not immediately clear
how to adapt Definition 1.1 to cover this case: usually H Sobolev spaces are defined for
p < 1 only.
In the next proposition we consider the relation of W 1,1 with Lipschitz functions. We
omit, for brevity, the simple proof, based once more on convolutions.
Since H 1,p (⌦) is defined by means of approximation by regular functions, for which
(1.2) is just the elementary “integration by parts formula”, it is clear that H 1,p (⌦) ⇢
W 1,p (⌦); in addition, the same argument shows that the weak derivative of u 2 H 1,p (⌦),
in the sense of W Sobolev spaces, is precisely the strong Lp (⌦, Rn ) limit of ruh , where
uh 2 C 1 (⌦) are strongly convergent to u. This allows to show by approximation some
basic calculus rules in H Sobolev spaces for weak derivatives, as the chain rule
0
r( u) = (u)ru 2 C 1 (R) Lipschitz with (0) = 0, u 2 H 1,p (⌦) , (1.7)
and, with a little more e↵ort (because one has first to show using the chain rule that
bounded H 1,p functions can be strongly approximated in H 1,p by equibounded C 1 (⌦)
functions) the Leibniz rule
On the other hand, we don’t have to prove the same formulas for the W Sobolev spaces:
indeed, using convolutions and a suitable extension operator described below (in the case
⌦ = Rn the proof is a direct application of (1.5), since in this case ⌦" = Rn ), one can
prove the following result:
5
With the word regular we mean that ⌦ is the epigrah of a Lipschitz function of (n 1)-
variables, written in a suitable system of coordinates, near to any boundary point.
However the equality H = W is not true in general, as the following example shows.
Example 1.6. In the Euclidean plane R2 , consider the open unit ball x2 +y 2 < 1 deprived
of one of its radii, say for instance the segment ⌃ given by ( 1, 0] ⇥ {0} . We can define on
this domain ⌦ a function v having values in ( ⇡, ⇡) and representing the angle in polar
coordinates. Fix an exponent 1 p < 2. It is immediate to see that v 2 C 1 (⌦) and that
its gradient is p-integrable, hence v 2 W 1,p (⌦). On the other hand, v 2 / H 1,p (⌦) because
the definition we have given would require the existence of regular approximations for
v up to the boundary: more precisely, one can easily show, using Fubini’s theorem and
polar coordinates, that any u 2 H 1,p (⌦) satisfies
1,p
! 7! u(rei! ) 2 Wloc (R) for L 1 -a.e. r 2 (0, 1). (1.10)
0 1 1,p
Indeed, if un 2 C
P(⌦) \ C (⌦) converge to u strongly in H (⌦) and (possibly extracting
a subsequence) n krun+1 run kp < 1, for all 2 (0, 1) the inequality |@✓ v| |rv|/r
gives
Z 1 X ✓Z ⇡ ◆1/p X
@un+1 @un p
| | d✓ dr 1 1/p krun+1 run kp < 1.
n ⇡ @✓ @✓ n
Since > 0 is arbitrary, it follows that for L 1 -a.e. r 2 (0, 1) the 2⇡-periodic continuous
functions ✓ 7! un (rei✓ ) have derivatives strongly convergent in Lploc (R), and therefore
(by the fundamental theorem of calculus) are equicontinuous. Any limit point of these
functions in Lploc (R) must then be 2⇡-periodic, continuous and W 1,p . If, by contradiction,
we take u = v, a similar Fubini argument shows that un (rei✓ ) converge in Lp ( ⇡, ⇡) to
the function v for L 1 -a.e. r 2 (0, 1). But, the function v(r, ✓) = ✓ 2 ( ⇡, ⇡) has no
continuous 2⇡-periodic extension. Therefore we get a contradiction and v can’t be in
H 1,p (⌦).
Remark 1.7. Taking into account the example above, we mention the Meyers-Serrin
theorem [24], ensuring that, for any open set ⌦ ⇢ Rn and 1 p < 1, the identity
W 1,p
C 1 (⌦) \ W 1,p (⌦) = W 1,p (⌦) (1.11)
6
be true unconditionally. In the case p = 1, the construction in the Meyers-Serrin theorem
provides for all u 2 W 1,1 (⌦) a sequence (un ) ⇢ C 1 (⌦) converging to u uniformly in ⌦,
with sup⌦ |run | convergent to kruk1 . Again, this might lead to a definition of H 1,1 for
which H 1,1 = W 1,1 unconditionally.
As it will be clear soon, we also need to define an appropriate subspace of H 1,p (⌦) in
order to work with functions vanishing at the boundary.
Definition 1.8. Given ⌦ ⇢ Rn open, we define H01,p (⌦) to be the completion of Cc1 (⌦)
with respect to the W 1,p norm.
It is clear that H01,p (⌦), being complete, is a closed subspace of H 1,p (⌦). Notice also
that H 1,p (Rn ) coincides with H01,p (Rn ). To see this, suffices to show that any function
u 2 C 1 (Rn ) with both |u| and |ru| in Lp (Rn ) belongs to H01,p (Rn ). We can indeed
approximate any such function u, strongly in H 1,p norm, by the compactly supported
functions R u, where R : Rn ! [0, 1] are smooth, 2-Lipschitz, identically equal to 1 on
B R and identically equal to 0 on Rn \ B R+1 .
We now turn to some classical inequalities.
Theorem 1.9 (Poincaré inequality, first version). Let ⌦ ⇢ Rn be an open bounded set
and p 2 [1, 1). Then there exists a constant C(⌦, p), depending only on ⌦ and p, such
that
kukLp C(⌦, p) krukLp 8u 2 H01,p (⌦). (1.12)
In addition C(⌦) C(n, p)diam(⌦).
The proof of this result can be simplified by means of these properties:
• H01,p (⌦) ⇢ H01,p (⌦0 ) if ⌦ ⇢ ⌦0 (monotonicity);
7
Hölder’s inequality gives
Z 1 p
p 0 @u
|u| (x1 , x ) 2 p 1
(t, x0 )dt (1.15)
1 @x1
and hence we just need to integrate w.r.t. x1 to get
Z 1 Z 1 p
p 0 @u
|u| (x1 , x ) dx1 2 p
(t, x0 )dt. (1.16)
1 1 @x1
Now, integrating w.r.t. x0 , repeating the previous argument for all the variables xj , j =
1, . . . , n and summing all such inequalities we obtain the thesis with C(Q, p) 2/n1/p .
⇤
Theorem 1.10 (Rellich). Let ⌦ be an open bounded subset with regular boundary and let
p 2 [1, 1). Then the immersion W 1,p (⌦) ,! Lp (⌦) is compact.
We do not give a complete proof of this result. Instead, we observe that it can be
obtained using an appropriate linear and continuous extension operator
T : W 1,p (⌦) ! W 1,p (Rn ) (1.17)
such that 8
< Tu = u in ⌦;
:
supp(T u) ⇢ ⌦0 ,
being ⌦0 a fixed bounded domain in Rn containing ⌦. When ⌦ is an halfspace the operator
can be achieved simply by a reflection argument; in the general case the construction relies
on the fact that the boundary of @⌦ is regular and so can be locally straightened by means
of Lipschitz maps (we will use these ideas later on, treating the boundary regularity of
solutions to elliptic PDE’s). The global construction is then obtained thanks to a partition
of unity.
The operator T allows basically a reduction to the case ⌦ = Rn , considered in the
next theorem.
Theorem 1.11. The immersion W 1,p (Rn ) ,! Lploc (Rn ) is compact, namely if (uk ) ⇢
W 1,p (Rn ) is bounded, then (un ) has limit points in the Lploc (Rn ) topology, and any limit
point belongs to W 1,p (Rn ).
Remark 1.12. It should be noted that the immersion W 1,p (Rn ) ,! Lp (Rn ) is obviously
continuous, but certainly not compact: just take a fixed element in W 1,p (Rn ) and sup-
ported in the unit square and consider the sequence of its translates along vectors ⌧h
with |⌧h | ! 1. Of course this is a bounded sequence in W 1,p (Rn ) but no subsequence
converges in Lp (Rn ) (indeed, all functions have the same Lp norm, while it is easily seen
that their Lploc limit is 0).
8
Let us now briefly sketch the main points of the proof of this theorem, since some of
the ideas we use here will be often considered in the sequel.
Proof. Basically, it is enough to prove that a bounded family F ⇢ W 1,p (Rn ) is totally
bounded in Lploc (Rn ). To obtain this, observe firstly that given any Borel domain A ⇢ Rn
and any ' 2 W 1,p (A|h| ) we have
where A|h| is the |h| neighbourhood of the set A and ⌧h '(x) = '(x + h). By approxima-
tion, we can assume with no loss of generality that ' 2 C 1 (A|h| ). The inequality (1.18)
follows by the elementary representation
Z 1
(⌧h ' ')(x) = hr'(x + sh), hi ds (1.19)
0
since
Z Z 1
k⌧h ' 'kpLp (A) |hr'(x + sh), hi|p ds dx (1.20)
A 0
Z 1Z
p
|h| |r'(y)|p dy ds = |h|p kr'kpLp (A|h| ) (1.21)
0 A|h|
where we used the Cauchy-Schwarz inequality and Fubini’s theorem. Hence, denoting by
(⇢" )">0 any rescaled family of smooth mollifiers such that supp(⇢" ) ⇢ B(0, "), we have
that for any R > 0
sup k' ⇤ ⇢" 'kLp (BR ) ! 0 (1.22)
'2F
for " ! 0. In fact, since ' ⇤ ⇢" is a mean, weighted by ⇢" , of translates of '
Z
' ⇤ ⇢" = ⌧ y '⇢" (y) dy ,
To conclude we just need to observe that the regularized family {' ⇤ ⇢" , ' 2 F} is rela-
tively compact in Lploc (Rn ) for any fixed " > 0. But this is easy since the Young inequality
implies
sup |' ⇤ ⇢" | k'kL1 (BR+" ) k⇢" k1 (1.24)
BR
9
and similarly
sup |r(' ⇤ ⇢" )| k'kL1 (BR+" ) kr⇢" k1 (1.25)
BR
so the claim follows by means of the Ascoli-Arzelá theorem. Notice that we used the
gradient bounds on elements of F only in (1.23). ⇤
We also need to mention another inequality due to Poincaré.
Theorem 1.13 (Poncaré inequality, second version). Let us consider a bounded, regular
and connected domain ⌦ ⇢ Rn and an exponent 1 p < 1, so that by Rellich’s theorem
we have the compact immersion W 1,p (⌦) ,! Lp (⌦). Then, there exists a constant C(⌦, p)
such that Z Z
p
|u u⌦ | dx C |ru|p dx 8u 2 W 1,p (⌦) (1.26)
⌦ ⌦
R
where u⌦ = ⌦
u dx.
Proof. By contradiction, if the desired inequality were not true, exploiting its homogene-
ity and translation invariance we could find a sequence (un ) ⇢ W 1,p (⌦) such that
Remark 1.14. It should be observed that the previous proof, even though very simple, is
far from giving the sharp constant for the Poincaré inequality (1.26). The determination
of the sharp constant is a difficult problem, solved only in very special cases (for instance
on intervals of the real line and p = 2, by Fourier analysis). Many more results are instead
available for the sharp constant in the Poincaré inequality (1.12).
10
2 Variational formulation of some PDEs
After the introductory section, whose main purpose was to fix the notation and recall
some basic tools of the theory of Sobolev spaces, we are now ready to discuss some basic
elliptic PDEs.
Let us consider the generalised Poisson equation
8 P
< u=f ↵ @↵ f↵ in ⌦;
(2.1)
:
u 2 H01,2 (⌦)
with data f, f↵ 2 L2 (⌦) for some fixed bounded and regular domain ⌦. This equation
has to be intended in a weak sense, that is, we look for u 2 H01,2 (⌦) satisfying
Z Z X
hru, r'i dx = (f ' + f↵ @↵ ') dx 8' 2 Cc1 (⌦). (2.2)
⌦ ⌦ ↵
Equivalently, by continuity of the bilinear form and density of Cc1 (⌦), the previous con-
dition could be requested for any ' 2 H01,2 (⌦).
In order to obtain existence
R weP just need to apply Riesz’s theorem to the associated
linear functional F (v) = ⌦ (f v + ↵ f↵ @↵ v) dx on the Hilbert space H01,2 (⌦) endowed
with the scalar product Z
(u, v) = hru, rvi dx (2.3)
⌦
which is equivalent to the usual one thanks to the Poincaré inequality (first version) proved
in Theorem 1.9.
We can consider many variants of the previous problem, basically by introduction of
one or more of the following elements:
• more general di↵erential operators instead of ;
• inhomogeneous or mixed boundary conditions;
• systems instead of single equations.
Our purpose now is to briefly discuss each of these situations.
:
u 2 H01,2 (⌦)
11
where, as before f, f↵ 2 L2 (⌦), and A 2 Rn⇥n is a constant matrix satisfying the following
requirements:
(i) A↵ 2 Rn⇥n is symmetric, that is A↵ = A ↵ ;
(ii) A has only positive eigenvalues, equivalently, A cI for some c > 0, in the sense of
quadratic forms.
Here and in the sequel we use the capital letter I to denote the identity n ⇥ n matrix. It
is not difficult to show that a change of independent variables, precisely u(x) = v(A 1 x),
transforms this problem into one of the form (2.1). For this reason it is convenient to deal
immediately with the case of a non-constant matrix A(x) 2 Rn⇥n satisfying:
(i) A is a Borel and L1 function defined on ⌦;
As indicated above, the previous problem has to be intended in weak sense and precisely
Z Z X
hAru, r'i dx = (f ' + f↵ @↵ ') dx 8' 2 Cc1 (⌦). (2.5)
⌦ ⌦ ↵
By continuity and density, also in this case it is equivalent to require the validity of the
identity above for all ' 2 H01,2 (⌦). In order to obtain existence we could easily modify
the previous argument when |A| 2 L1 (⌦), using the equivalent scalar product
Z X
hu, vi := A↵ @↵ u@ v dx .
⌦ ↵,
However, in order to include also unbounded A’s, thus dropping assumption (i), we prefer
here to proceed di↵erently and introduce some ideas that belong to the so-called direct
method of the Calculus of Variations. Let us consider the functional F : H01,2 (⌦) ! R
Z Z XZ
1
F (v) = hArv, rvi dx f v dx f↵ @↵ v dx. (2.6)
⌦ 2 ⌦ ↵ ⌦
First we note that, thanks to the assumption (2.4) on A, for all " > 0 it holds
Z Z X Z
c 2 1 2 2 ✏
F (v) |rv| dx (|f | + |f↵ | ) dx v 2 + |rv|2 dx.
2 ⌦ 2" ⌦ ↵
2 ⌦
12
Choosing " < c/2 we get
Z Z X Z
c 1 ✏
F (v) |rv|2 dx 2
(|f | + 2
|f↵ | ) dx v 2 dx
4 ⌦ 2" ⌦ ↵
2 ⌦
and now, thanks to the Poincaré inequality, we can choose possibly " even smaller to get
Z Z X
c 2 1
F (v) |rv| dx (|f |2 + |f↵ |2 ) dx.
8 ⌦ 2" ⌦ ↵
and consequently, in order to look for its minimum points it is enough to consider a closed
ball of H01,2 (⌦). Now, take any minimizing sequence (un ) of F : since H01,2 (⌦) is (being
Hilbert) a reflexive space we can assume, possibly extracting a subsequence, that un * u
for some u 2 H01,2 (⌦). Using Fatou’s lemma and the fact that uh ! u in H 1,2 implies
ruh(k) ! ru a.e. in ⌦ for a suitable subsequence h(k), it is not difficult to prove that
F is lower semicontinuous (we shall also prove this in Theorem 3.2, in a more general
framework). In addition, F is convex, being the sum of a linear and a convex functional.
It follows that F is also weakly lower semicontinuous, hence
and the desired result follows. Even in the case when |A| 2 L1loc we can still di↵erentiate
the functional, but a priori only along directions in ' 2 Cc1 (⌦), and recover the weak
formulation of our PDE.
13
2.2 Inhomogeneous boundary conditions
We now turn to study the boundary value problem for u 2 H 1,2 (⌦)
8 P
< u=f ↵ @↵ f↵ in ⌦;
:
u=g on @⌦
with f, f↵ 2 L2 (⌦) and a suitable class of functions g 2 L2 (@⌦). Since the immersion
H 1,2 (⌦) ,! C(⌦) does not hold if n 2, the boundary condition has to be considered in
the weak sense described below.
Here and in the sequel, unless otherwise stated, we indicate with ⌦ an open, bounded
and regular subset of Rn .
Theorem 2.1. For any p 2 [1, 1) the restriction operator
satisfies kT ukLp (@⌦) C(p, ⌦)kukW 1,p (⌦) . Therefore it can be uniquely extended to a linear
and continuous operator from W 1,p (⌦) to Lp (@⌦).
Proof. We prove the result only in the case when ⌦ is the subgraph of a Lipschitz function
f inside the rectangle ⌦0 ⇥ (a, b), with ⌦ ⇢ Rn 1 open, with a0 = inf f > a, proving the
estimate on the portion
:= {(x0 , f (x0 )) : x0 2 ⌦0 }
of its boundary (here we use the notation x = (x0 , t) with x0 2 ⌦0 and t 2 (a, b)). The
general case can be easily achieved by a partition of unity argument.
By the fundamental theorem of calculus, for all t 2 (0, a0 a) we have
Z f (x0 ) p Z f (x0 )
0 0 0 0 0
p
|u(x , f (x ) t) u(x , f (x ))| @xn u(x , r) dr (b a) p 1
|@xn u(x0 , r)|p dr .
f (x0 ) t a
14
Now we average this estimate with respect to t 2 (0, a0 a), together with the fact
that the determinant of the gradient of the map (x0 , t) 7! (x0 , f (x0 ) t) is identically equal
to 1, to get
Z Z Z
1 p 2p 1 p p 1 p 1
p |u| d 0 |u| dx + 2 (b a) |@xn u|p dx .
1 + L2 a a ⌦ ⌦
Because of the previous result, for u 2 W 1,p (⌦) we will interpret the boundary condi-
tion u|@⌦ = g as
T u = g. (2.11)
It can also be easily proved that T u is characterized by the identity
Z Z Z
@' @u
u dx = ' dx + ' T u ⌫i d 8' 2 C 1 (⌦) (2.12)
⌦ @xi ⌦ @xi @⌦
where ⌫ = (⌫1 , . . . , ⌫n ) is the unit normal vector, pointing out of ⌦. Indeed, using the
equality H 1,p (⌦) = W 1,p (⌦) of Theorem 1.5 one can start from the classical divergence
theorem with u 2 C 1 (⌦) and then argue by approximation.
Remark 2.2. It is possible to show that the previously defined restriction operator T is
not surjective if p > 1 and that its image can be described in terms of fractional Sobolev
spaces W s,p , characterized by the finiteness of the integral
Z Z
|u(x) u(y)|p
dxdy ,
|x y|n+sp
see [1], with s = 1 1/p. The borderline case p = 1 is special, and in this case Gagliardo
proved in [13] the surjectivity of T .
We can now mimic the argument described in the previous section in order to achieve
an existence result, provided the function g belongs to the image of T, that is there exists
e 2 W 1,2 (⌦) such that T u
a function u e = g. Indeed, if this is the case, our problem is
reduced to show existence of a solution for the equation
8 P
< v = fe e
↵ @↵ f↵ in ⌦;
:
v 2 H01,2 (⌦) .
where fe = f and fe↵ = f↵ @↵ u e. This is precisely the first problem we have discussed
above and so, denoted by v its unique solution, the function u = v + ue will satisfy both
our equation and the required boundary condition.
15
Finally, let us discuss the so-called Neumann boundary conditions, involving the be-
haviour of the normal derivative of u on the boundary. We consider a problem of the
form 8 P P
↵
< ↵, @↵ (A @ u) + u = f ↵ @↵ f↵ in ⌦;
:
A↵ @ u⌫↵ = g on @⌦
with A↵ a real matrix and > 0 a fixed constant. For the sake of brevity, we just discuss
the case A↵ = ↵ so that the problem above becomes
8
>
< u + u = f in ⌦;
: @u = g
>
on @⌦.
@⌫
In order to give it a clear meaning, note that if u, v 2 C 1 (⌦) then
Z Z Z
@u
hru, rvi dx = v u dx + v d (2.13)
⌦ ⌦ @⌦ @⌫
and so in this case it is natural to ask that for any v 2 C 1 (⌦) the desired solution u
satisfies Z Z Z
[hru, rvi + uv] dx = vf dx + vg d . (2.14)
⌦ ⌦ @⌦
In order to obtain existence (and uniqueness) for this problem when g 2 L2 (@⌦), it is
enough to apply Riesz’s theorem to the bilinear form on H 1,2 (⌦)
Z
a(u, v) = [hru, rvi + uv] dx (2.15)
⌦
16
the form A cI for some c > 0) to the vector-valued case. The first idea is to define the
Legendre condition X ↵
Aij ⇠↵i ⇠ j c |⇠|2 8⇠ 2 Rm⇥n (2.16)
↵, ,i,j
where R m⇥n
indicates, as above, the space of m ⇥ n real matrices. Let us apply it in order
to obtain existence and uniqueness for the system
8 P ↵ j
P ↵
< ↵, ,j @↵ (Ai,j @ u ) = fi ↵ @ ↵ fi i = 1, . . . , m
:
u 2 H01 (⌦; Rm )
with data fi , fi↵ 2 L2 (⌦).1 The weak formulation of the problem is obviously
Z X Z "X X
#
A↵ij @ uj @↵ 'i dx = f i 'i + fi↵ @↵ 'i dx (2.17)
⌦ i,j,↵, ⌦ i i,↵
m
for every ' 2 [Cc1 (⌦)] and i = 1, . . . , m. Now, if the matrix A↵ij is symmetric with respect
to the transformation (↵, i) ! ( , j) (which is implied, for instance, by the symmetries
in (↵, ) and (i, j)), then it defines a scalar product on H01 (⌦; Rm ) by the formula
Z X
(', ) = A↵ij @↵ 'i @ j dx. (2.18)
⌦ i,j,↵,
If, moreover, A satisfies the Legendre condition (2.16) for some c > 0, it is immediate to
see that this scalar product is equivalent to the standard one (with A↵i,j = ↵ ij ) and so
we are led to apply again Riesz’s theorem to conclude the proof.
From now on, we will often adopt Einstein’s summation convention on repeated in-
dices, using it without explicit mention.
It should be noted that in the proof of some existence result (and, in particular, in
the scalar case) the symmetry hypothesis w.r.t. the transformation (↵, i) ! ( , j) is not
necessary, since we can exploit the following:
Theorem 2.3 (Lax-Milgram). Let H be a (real) Hilbert space and let a : H ⇥ H ! R a
bilinear, continuous and coercive form so that
a(u, u) |u|2 8u 2 H ,
0
for some > 0. Then for any F 2 H there exists uF 2 H such that a(uF , v) = F (v) for
all v 2 H.
1
Note that we sometimes omit the Sobolev exponent when this is equal to two: for instance H01 (⌦)
stands for H01,2 (⌦).
17
Proof. By means of the standard Riesz’s theorem it is possible to define a linear operator
T : H ! H such that
a(u, v) = hT u, vi 8u, v 2 H
and such T is continuous since
where C is a constant of continuity for a(·, ·) and hence kT k C. Now we introduce the
auxiliary bilinear form
a(u, v) = hT T ⇤ u, vi = hT ⇤ u, T ⇤ vi ,
e
which is obviously symmetric and continuous. Moreover, thanks to the coercivity of a we
have that e
a is coercive too, because
p
kuk2 a(u, u) = hT u, ui = hu, T ⇤ ui kuk kT ⇤ uk = kuk ea(u, u)
and so e
a(u, u) 2
kuk2 . Since e
a determines an equivalent scalar product on H we can
apply again Riesz Theorem to obtain a vector u0F 2 H such that
a(u0F , v) = F (v) 8 v 2 H .
e
⇤
As indicated above, we now want to formulate a di↵erent notion of ellipticity for the
vector case. To this aim, it is useful to analyse more in detail the scalar case. We have
the two following conditions:
It is obvious by integration that (E) ) (C) and we may wonder about the converse
implication. As we will see below, this holds in the scalar case (m = 1) and fails in the
vectorial case (m > 1).
Proposition 2.4. Let (C) and (E) as above. Then, (C) is equivalent to (E).
18
Proof. Let is prove that (C) implies (E). The computations become more transparent if
we work with functions having complex values, and so let us define for any u, v 2 H01 (⌦; C)
Z Z n
X
⌦ ↵
aA (u, v) = Aru, rv dx = A↵ @↵ u@ u dx .
⌦ ⌦ ↵, =1
A simple computation shows that (here ru 2 Cn stands for r<u + ir=u, where <u and
=u are respectively the real and imaginary part of u)
Now consider a function ' 2 Cc1 (⌦) and define u⌧ (x) = '(x)ei⌧ x·⇠ . We have that
Z Z
1 2 ↵ ↵
<aA (u⌧ , u⌧ ) = ' A ⇠↵ ⇠ dx + o⌧ = A ⇠↵ ⇠ '2 dx + o⌧
⌧2 ⌦ ⌦
which immediately implies the thesis (it is enough to choose ' not identically zero).
Actually, every single part of our discussion is still true in the case when A↵ = A↵ (x)
is Borel and L1 function in ⌦ and we can conclude that (E) holds for a.e. x 2 ⌦: we just
need to choose, in the very last step, an appropriate sequence of rescaled and normalized
mollifiers concentrating around x0 , for any Lebesgue point x0 of A. The conclusion comes,
in this situation, by Lebesgue di↵erentiation theorem.
For the reader’s convenience we recall here some basic facts concerning Lebesgue points
(see also Section 13). Given f 2 L1loc (Rn ) and x0 2 Rn , we say that x0 is a Lebesgue point
for f if there exists 2 R such that
Z
lim |f (y) | dy = 0. (2.21)
r#0 Br (x0 )
19
In this case is unique and it is sometimes written
= fe(x0 ) = g
lim f (x). (2.22)
x!x0
Notice that both the notion of Lebesgue point and fe are invariant in the Lebesgue
equivalence class of f . The Lebesgue di↵erentiation theorem says that for L n -a.e. x0 2 Rn
the following two properties hold: x0 is a Lebesgue point of f and fe(x0 ) = f (x0 ). Notice
however that the validity of the second property at a given x0 does depend on the choice
of a representative of f in the Lebesgue equivalence class.
Going back to the previous discussion, it is very interesting to note that the argument
above does not give a complete equivalence when m > 1: in fact, the coercivity condition
Z
aA (u, u) |ru|2 dx u 2 H 1 (⌦; Rm ) (2.23)
⌦
can be applied to test functions having the form u⌧ (x) = '(x)bei⌧ x·a with a 2 Rn and
b 2 Rm and implies the Legendre-Hadamard condition
that is the Legendre condition restricted to rank one matrices ⇠↵i = a↵ bi . Explicit exam-
ples show that the Legendre-Hadamard condition is in general strictly weaker than the
Legendre condition.
with " 0. Since t 7! det(M + tN ) is linear for any rank one matrix N , the Legendre-
Hadamard condition with = " is fulfilled. On the other hand our quadratic form,
restricted to diagonal matrices with eigenvalues t and t, equals
t2 + 2t2 " .
It follows that the Legendre condition with = 0 fails when 2" < 1.
Theorem 2.6 (Gårding). Assume that A↵ij satisfies the Legendre-Hadamard condition
for some positive constant and the symmetry condition A↵ij = A ↵ ji. Then aA (u, u)
R
|ru|2 dx for all u 2 H 1 (Rn ; Rm ).
20
In the proof of Gårding’s theorem, we denote by S(Rn ) the Schwartz space of smooth
C-valued functions that decay at 1 faster than any polynomial, and by ' b and ' e the
Fourier transform of ' and its inverse, respectively
Z
n/2
b = (2⇡)
'(⇠) '(x)e ix·⇠ dx (2.26)
and Z
n/2
e
'(x) = (2⇡) '(⇠)eix·⇠ d⇠ . (2.27)
Proof. By density it is enough to prove the result when u 2 [Cc1 (Rn )]m . In this case we
use the representation Z
n/2
u(⇠) = (2⇡) '(x)e ix·⇠ dx ,
Rn
m
b
that is u(⇠) = '(⇠) for some ' 2 [S(Rn )] . Consequently,
@↵ uj (⇠) = ix[ j
↵' ,
hence
Z Z Z
@uj @ul
aA (u, u) = A↵jl d⇠ = i 2
A↵jl x[ j dl
↵ ' x ' d⇠ = A↵jl (x↵ 'j )(x 'l ) dx ,
Rn @⇠↵ @⇠ Rn Rn
the last passage being due to the Plancherel identity (2.28). Now we can apply our
hypothesis to get
A↵jl a↵ bj a bl |a|2 |b|2
with a = x 2 Rn and b = ' 2 Cn , so that
Z
aA (u, u) |x|2 |'(x)|2 dx . (2.29)
Rn
If we perform the same steps with ↵ jl in place of A↵jl we see at once that
Z Z
2
|ru| (⇠) d⇠ = |x|2 |'(x)|2 dx . (2.30)
Rn Rn
21
Remark 2.7. Gårding’s theorem marks in some sense the di↵erence between pointwise
and integral inequalities. It is worth mentioning some related inequalities that are typi-
cally nonlocal, namely, they do not arise from the integration of a pointwise inequality.
An important example is Korn’s inequality:
Z Z p
p ru + (ru)t
|ru| dx c(n, p) dx for all u 2 Cc1 (Rn , Rn ) , (2.31)
R n R n 2
for p 2 (1, 1). A variant of this example is the Korn-Poincaré inequality: if ⌦ is an open,
bounded and regular set in Rn and p 2 (1, 1), then
Z Z p
p ru + (ru)t
inf |u(x) Bx c| dx C(⌦, p) dx . (2.32)
c2Rm , B t = B ⌦ ⌦ 2
We will make the following standard assumptions on the Lagrangian: we assume that
L : ⌦ ⇥ Rm ⇥ Rm⇥n ! R is Borel and, denoting the variables as (x, s, p), we assume that
L is of class C 1 in (s, p) with
for any domain K = ⌦0 ⇥ {(s, p)| |s| + |p| R} with R > 0 and ⌦0 b ⌦. In this case it is
possible to show that the map
Z
t 7! L(x, u + t'ru + tr') dx
⌦0
1,1
is of class C 1 for all u, ' 2 Wloc (⌦; Rm ) and ⌦0 b ⌦, and its derivative equals
Z
Ls (x, u + t', ru + tr') · ' + Lp (x, u + t', ru + tr') · r' dx
⌦0
22
(the assumption (2.34) is needed to di↵erentiate under the integral sign). As a conse-
quence, if a locally Lipschitz function u is a local minimizer and ' 2 Cc1 (⌦0 ; Rm ), since
F (u, ⌦0 ) F (u + t', ⌦0 ) we can di↵erentiate at t = 0 to obtain
Z "X X i
#
@'
Lsi (x, u, ru)'i + Lp↵i (x, u, ru) dx = 0 . (2.35)
⌦ 0
i ↵,i
@x ↵
Hence, exploiting the arbitrariness of ', we obtain the Euler-Lagrange equations in the
weak sense:
@
Lp↵ (x, u, ru) = Lsi (x, u, ru) i = 1, 2, . . . , m .
@x↵ i
Exploiting this idea, we can associate to many classes of PDEs appropriate energy func-
tionals, so that the considered problem is nothing but the Euler-Lagrange equation for
the corresponding functional. For instance, neglecting the boundary conditions (that can
actually be taken into account by an appropriate choice of the ambient functional space),
equations having the form
u = g(x, u) (2.36)
derive from the functional
Z s
1
L(x, s, p) = |p|2 g(x, r) dr . (2.37)
2 0
Adding stronger hypotheses on the Lagrangian L, in analogy with what has been done
above, i.e. requiring that
for any domain K = ⌦0 ⇥{(s, p)| |s| + |p| R} with ⌦0 b ⌦, we can find another necessary
minimality condition corresponding to
d2
F (u + t') 0,
dt2 t=0
namely
Z
0 (', ') = [Ar'r' + Br' · ' + C' · '] dx 8' 2 Cc1 (⌦; Rm ) , (2.38)
⌦
where the dependence on x and all indices are omitted for brevity and
8
>
<A(x) = Lpp (x, u(x), ru(x)) ;
B(x) = Lps (x, u(x), ru(x)) ; (2.39)
>
:
C(x) = Lss (x, u(x), ru(x)) .
23
We can finally obtain pointwise conditions on the local minimizer u by means of the
following theorem, whose proof can be obtained arguing as in the proof that coercivity
implies ellipticity (Proposition 2.4).
where A = A↵ij (x), B = Bij↵ (x) and C = Cij (x) are Borel and L1 functions. If ⇥(u, u)
0 for all u 2 H01 (⌦; Rm ) then A(x) satisfies the Legendre-Hadamard condition with = 0
for a.e. x 2 ⌦.
Hence, in our case, we find that Lpp (x, u(x), ru(x)) satisfies the Legendre-Hadamard
condition with = 0 for a.e. x 2 ⌦.
Theorem 3.1 (Dunford-Pettis). Let (X, A, µ) be a finite measure space and let F ⇢
L1 (X, A, µ). Then the following facts are equivalent:
(i) the family F is sequentially relatively compact w.r.t. the weak-L1 topology;
(ii) there exists : [0, 1) ! [0, 1], with (t)/t ! +1 as t ! 1, such that
Z
(|f |) dµ 1 8f 2 F ;
X
24
Theorem 3.2 (Tonelli). Let L(x, s, p) : ⌦ ⇥ Rm ⇥ Rm⇥n be a Lagrangian satisfying the
following properties:
(1) L is non-negative;
(2) L is lower semicontinuous w.r.t. s and the partial derivatives Lp↵i exist and are
continuous w.r.t. s;
(3) p 7! L(x, s, ·) is convex2 .
Then any sequence (uh ) ⇢ W 1,1 (⌦; Rm ) converging to u in L1 (⌦; Rm ) with |ruh | uniformly
integrable satisfies the lower semicontinuity inequality
lim inf F (uh ) F (u) . (3.2)
h!1
Proof. We start by noticing that there is a subsequence uh(k) such that
lim inf F (uh ) = lim F (uh(k) )
h!1 k!1
25
Hence, the weak convergence ruh(k) * ru ensures that
Z
lim inf F (uh ) L(x, u(x), ru(x)) dx
h!1 K"
This proves that the dependence from ⌦ of this notion is only seeming. Another way to
see this relies on the observation that whenever (3.3) is valid for ⌦, then:
26
• it is valid for every ⌦0 ⇢ ⌦, thanks to the previous observation;
• it is valid for x0 + ⌦, for x0 2 Rn and > 0, considering the transformation
'(x) 7! '(x0 + x)/ .
Finally, a simple approximation argument gives
Z
F (A + r') dx F (A) 8 ' 2 Cc1 (⌦; Rm ) . (3.4)
⌦
Notice that the inequality above makes sense: either F ⌘ +1 or it is finite at least
one point. In the second case the negative part of F has at most linear growth and the
integral in the left hand side makes sense, finite or infinite.
Now, let f 2 L1 (⌦, Rm⇥n ) and consider the law µ of the map f with respect to the
rescaled Lebesgue measure L n /L n (⌦). If F : Rm⇥n ! R[{+1} is lower semicontinuous
and convex, thanks to Jensen’s inequality one has
Z Z ✓Z ◆ ✓Z ◆
F (f (x)) dx = F (y) dµ(y) F y dµ(y) = F f dx . (3.6)
⌦ Rm⇥n Rm⇥n ⌦
27
Remark 3.8. The following chain of implications holds:
Hence, Example 2.5 provides a quasiconvex function that is not convex when n =
m = 2, and considering the determinant of a 2 ⇥ 2 minor the example fits also the
case min{m, n} 2;
Definition 3.9 (w⇤ -convergence in W 1,1 ). Let us consider an open set ⌦ ⇢ Rn and
fk 2 W 1,1 (⌦). We write fk ! f in w⇤ -W 1,1 (⌦) if
• fk ! f uniformly in ⌦;
This is a direct consequence of the fact that (rfk ) is sequentially compact in the w⇤ -
topology of L1 , and any weak⇤ limit provides a weak derivative of f (hence f 2 W 1,1 ,
the limit is unique and the whole sequence of derivatives w⇤ -converges). Obviously an
analogous statement holds for Rm -valued maps.
Before stating Morrrey’s lower semicontinuity theorem we give a quick proof of Rademacher’s
di↵erentiability theorem.
28
R
Proof. Fix a Lebesgue point x0 of the weak gradient ru, namely Br (x0 ) |ru L| dx ! 0
as r # 0 for some linear map L : Rn ! Rm . We shall prove that f is di↵erentiable at x0
and that the (classical) gradient rf at x0 is equal to L.
First of all, it is easy to see that this property can be equivalently stated as follows:
where fr (y) = (f (x0 + ry) f (x0 ))/r are the rescaled maps. Notice, that fr are equi-
Lipschitz in B 1 and equi-bounded (because fr (0) = 0), hence fr is relatively compact in
C 0 (B 1 ) as r # 0. Hence, suffices to show that any limit point f0 (y) = limi fri (y) coincides
with L(y). A simple change of variables shows that (understanding here gradients as weak
gradients!) Z Z
|rfr L| dy = |rf L| dx .
B1 Br (x)
29
• each Fr is lower semicontinuous at 0 with respect to the w⇤ -W 1,1 (Q; Rm ) conver-
gence, indeed
Z
1
Fr (v) = n L (x, u(x) + rv ((x x0 )/r) , ru(x) + rv ((x x0 )/r)) dx
r Qr (x0 )
✓ Z ◆
1
= n F (u + rv ((x x0 )/r)) L(x, u(x), ru(x)) dx ;
r ⌦\Qr (x0 )
• being x0 a Lebesgue point for ru, for any bounded sequence (vh ) ⇢ W01,1 (Q; Rm )
it is easily checked that the continuity of L gives
Given a test function ' 2 Cc1 (Q, Rm ), we work with the 1-periodic function such that
|Q = ' and the sequence of highly oscillating (because h 1 -periodic) functions
1
vh (x) := (hx) ,
h
⇤
which obviously converge uniformly to 0. Since rvh (x) = r (hx) we have also vh * 0
in W 1,1 (Q; Rm ), so that thanks to the lower semicontinuity of F0 at 0 one has
Z Z
n
H(0) = F0 (0) lim inf H(rvh (x)) dx = lim inf h H(r (y)) dy
h!1 Q h!1 Qh
Z Z
= H(r (y)) dy = H(r'(y)) dy ,
Q Q
30
Now, considering the decomposition
⇥ ⇤
L(x, uh (x), ruh (x)) = L(x, uh (x), ruh (x)) L(x, u(x), ruh (x)) + L(x, u(x), ruh (x))
we obtain
Z XZ p
n`
lim inf L1 (x, ruh (x)) dx L1 (xi , ru(x)) dx !( )
h!1 Q i Qi 2k
Z p
n`
L1 (x, ru(x)) dx 2!( ).
Q 2k
31
: ⌦ ! [0, 1] with compact support. We can apply the quasiconvexity inequality to
' = (uh u) and the local Lipschitz property with R(p) = L2 (p) L2 (0) to get
Z
R(A) R((1 )A + ruh + (uh u) ⌦ r ) dx
⌦ Z Z Z
C(|A| + kruh k1 ) (1 ) dx + Ckr k1 |uh u| dx + R(ruh ) ,
⌦ ⌦ ⌦
so that passing to the limit first as h ! 1 and then as " 1 gives the result.
Finally, we consider the generalR case, using Rademacher’s theorem and a blow-up
argument. Assume that the lim inf ⌦ L2 (ruh ) dx is a limit, that we call L, and consider
the family of measures µh := L2 (ruh )L n . Being this family bounded, we can assume
with no loss of generality that µh weakly converge, in the duality with Cc (⌦), to some
measure µ. Recall that the evaluations on compact sets K and open sets A are respectively
upper and lower semicontinuous w.r.t. weak convergence, i.e.
µ(K) lim sup µh (K), µ(A) lim inf µh (A) . (3.7)
h!1 h!1
32
The previous result shows that quasiconvexity of the Lagrangian is equivalent to se-
quential lower semicontinuity of the integral functional in the weak⇤ -W 1,1 convergence.
However, in many problems of Calculus of Variations only L↵ bounds, with ↵ < 1, are
available on the gradient. A remarkable improvement of Morrey’s result is the following:
Theorem 3.13 (Acerbi-Fusco). Suppose that a Borel Lagrangian L(x, s, p) is continuous
in (s, p) and satisfies
for some ↵ > 1 and some constant C. Suppose also that the map p 7! L(x, s, p) is
quasiconvex for all (x, s). Then F is sequentially lower semicontinuous w.r.t. the weak
W 1,↵ (⌦; Rm )-topology.
4 Regularity Theory
We begin studying the local behaviour of (weak) solutions of the system of equations
⇢
@↵ A↵ij @ uj = fi @↵ Fi↵ i = 1, . . . , m
(4.1)
u 2 Hloc (⌦; R )
1 m
with A↵ij 2 L1 (⌦), fi 2 L2loc (⌦) and Fi↵ 2 L2loc (⌦). From now on we shall use | · | for the
Hilbert-Schmidt norm of matrices and tensors, even though some estimates would still be
valid with the (smaller) operator norm.
Theorem 4.1 (Caccioppoli-Leray inequality). If the Borel coefficients A↵ij satisfy the
Legendre condition (L) with > 0 and
then there exists a positive constant c = c( , ⇤) such that for any ball BR (x0 ) b ⌦ and
any k 2 Rm it holds
Z Z Z Z
2 2 2 2 2
c |ru| dx R |u(x) k| dx+R |f (x)| dx+ |F (x)|2 dx .
BR/2 (x0 ) BR (x0 ) BR (x0 ) BR (x0 )
(4.2)
Before proceeding to the proof, some remarks are in order.
Remark 4.2. (1) The validity of (4.2) for all k 2 Rm depends on the translation in-
variance of the PDE. Moreover, the inequality (and the PDE as well) has a natural
scaling invariance: if we think of u as an adimensional quantity, then all sides have
dimension lengthn 2 , because f ⇠ length 2 and F ⇠ length 1 .
33
(2) The Caccioppoli-Leray inequality is meaningful because for a general function u the
gradient ru can not be controlled by the variance of u! Precisely because of this
fact we can expect that several useful (regularity) informations can be drawn from
it. We will see indeed that CL inequalities are very “natural” and useful in the
context of regularity theory.
Remark 4.3 (Absorbtion scheme). In the regularity theory it often happens that one
can estimate, for some ↵ < 1,
A BA↵ + C .
The absorption scheme allows to bound A in terms of B, C and ↵ only and works as
follows: by the Young inequality
b "p ap bq 1 1
ab = "a + q (with + = 1)
" p "q p q
for p = 1/↵ one obtains
"p A B q
A BA↵ + C + q +C.
p "q
"p 1
Now, if we choose " = "(p) sufficiently small, so that , we get
p 2
Bq
A2 + 2C .
"q q
Let us prove Theorem 4.1.
Proof. Without loss of generality, we can consider x0 = 0 and k = 0. As typical in
regularity theory, we choose test functions depending on the solution u itself, namely
:= u⌘ 2
where ⌘ 2 Cc1 (BR ), ⌘ ⌘ 1 in BR/2 , 0 ⌘ 1 and |r⌘| 4/R.
Since u solves (4.1), we have that
Z Z Z
Arur f F ·r =0 (4.3)
34
• By the Legendre condition
Z Z
⌘ 2
A↵ij i
@↵ u @ u j
⌘ 2 |ru|2 .
BR BR
• We have
Z Z Z
8⇤
2 ⌘Aru (u ⌦ r⌘) 2 ⌘|A||ru||u||r⌘| (⌘|ru|) |u|
R
Z Z
4⇤" 2 2 4⇤
⌘ |ru| + |u|2 ,
R R"
where the first estimate is due to Schwarz inequality, the second one relies on the
boundedness of coefficients A↵ij and the estimate on |r⌘|, and the third one is based
on the Young inequality.
• Again by the same arguments (Schwarz inequality, estimate on |r⌘| and Young
inequality)
Z Z Z Z
↵ i 8 2 4
2 ⌘|Fi u |@↵ ⌘| |F ||u| 4 |F | + 2 |u|2 .
BR R BR BR R BR
35
By choosing " sufficiently small, in such a way that 4⇤"/R = /4, one can absorb line
(4.6), and the thesis follows noticing that
Z Z
2 2
⌘ |ru| |ru|2 .
BR BR/2
Remark 4.4 (Widman’s hole-filling technique). There exists a sharper version of the
Caccioppoli-Leray inequality, let us illustrate it in the simpler case f = 0, F = 0. Indeed,
since
4
|r⌘| B \B ,
R R R/2
following the proof of Theorem 4.1 one obtains
Z Z
c
|ru(x)|2 dx 2 |u(x) k|2 dx . (4.7)
BR/2 R BR \BR/2
R
Setting k := BR/2 u, the Poincaré inequality in the domain B1 \ B 1/2 and a scaling
argument give Z Z
|ru(x)|2 dx c |ru(x)|2 dx . (4.8)
BR/2 BR \BR/2
R
Adding to (4.8) the term c BR/2
|ru(x)|2 dx, we get
Z Z
2
(c + 1) |ru(x)| dx c |ru(x)|2 dx .
BR/2 BR
k 1
Iterating (4.7) and interpolating (i.e. considering the integer k such that 2 R<r
2 k R), it is not difficult to get
Z ⇣ r ⌘↵ Z
2 ↵
|ru(x)| dx 2 |ru(x)|2 dx 0<rR (4.9)
Br R BR
with (1/2)↵ = ✓, i.e. ↵ = ln2 (1/✓). When n = 2, this implies that u 2 C 0,↵/2 , as we will
see.
The following is another example of “unnatural” inequality, which provides additional
informations on functions that satisfy it.
36
Definition 4.5 (Reverse Hölder’s inequality). Let ↵ 2 (1, 1). A non-negative function
f 2 L↵loc (⌦) satisfies a reverse Hölder’s inequality with exponent ↵ if there exists a constant
c > 0 such that Z ✓Z ◆↵
↵
f c f 8BR (x) b ⌦ .
BR/2 (x) BR (x)
For the sake of completeness, we now recall the Sobolev inequalities. Detailed proofs
will be provided later on: concerning the cases p = n and p > n, we will see them in the
more general context of Morrey’s theory. We will treat the case p < n while dealing with
De Giorgi’s solution of Hilbert’s XIX problem, since slightly more general versions of the
Sobolev inequality are needed there.
Theorem 4.6 (Sobolev inequalities). Let ⌦ be either the whole space Rn or a bounded
regular domain.
• If p < n, denoting with p⇤ := nnpp > p the Sobolev conjugate exponent (characterized
also by p1⇤ = p1 n1 ), we have the continuous immersion
⇤
W 1,p (⌦) ,! Lp (⌦) .
• If p > n,
W 1,p (⌦) ,! C 0,1 n/p
(⌦) .
⇤
coming from the continuity of the embedding W 1,p ,! Lp , we get
✓Z ◆1/p⇤ ✓Z ◆1/p
p⇤ p
|v v| c B1 |rv| ,
B1
3
The result is basically sharp, as the example of ( ln |x|)↵ 2 W 1,n (B1 ) for n > 1 and ↵ 2 (0, 1 1/n)
shows.
37
If u solves (4.1) with f = F = 0, combining (4.10) with the CL inequality when p⇤ = 2
(that is, p = 2n/(n + 2) < 2), we write
Z !1/2 ✓ Z ◆1/2 ✓Z ◆1/p
2 2 p
cCL R |ru| |u u| cR |ru| .
BR/2 BR BR
This way we proved that |ru|p satisfies a reverse Hölder’s inequality with exponent ↵ =
2/p > 1 and C = c/cCL , that is
Z !1/2 ✓Z ◆1/p
2 p
|ru| C |ru| .
BR/2 BR
Remark 4.8 (Embedding for higher order Sobolev spaces). Recall first that higher order
Sobolev spaces W k,p (⌦) are recursively defined (k 1 integer, 1 p 1)
W k,p (⌦) := u 2 W 1,p (⌦) : ru 2 W k 1,p
(⌦; Rn ) .
Together with the Sobolev embedding in Theorem 4.6, with p > n, another way to gain
continuity is using the Sobolev spaces W k,p , with k sufficiently large. In fact, we can
arbitrarily expand the chain
⇤ ⇤ )⇤
W 2,p ,! W 1,p ,! L(p .
Iterating the ⇤ operation k-times we get
1 1 k
= ,
p⇤···⇤ p n
therefore if k > [ np ] (where [·] denotes the integer part) we obtain W k,p ⇢ C 0,↵ with any
positive ↵ with ↵ < 1 n/p + [n/p].
38
in a weak sense. For any ball BR (x0 ) b ⌦, by the Caccioppoli-Leray inequality we get,
Z ✓ ◆ 2 Z 2 Z
@u c @u
r 2 + |f |2 . (4.11)
BR/2 (x0 ) @xi R BR (x0 ) @xi BR (x0 )
We have chosen the Poisson equation because constant coefficients di↵erential operators
commute with convolution, so in this case the a priori regularity assumption can be a
posteriori removed. Indeed, estimate (4.11) applies to u ⇤ ⇢" with f ⇤ ⇢" in place of f ,
since u ⇤ ⇢" satisfies
(u ⇤ ⇢" ) = f ⇤ ⇢" .
2
Passing to the limit as " ! 0 we obtain that u 2 Hloc (⌦) and that the same inequality
1
holds for u, starting from the assumption u 2 Hloc (⌦).
The situation is much more complex when the coefficients A↵ij are not constant and
therefore di↵erentiation provides a worse right hand side in the PDE. Nirenberg’s idea is
to introduce partial discrete derivatives
In the next lemma we show that membership to W 1,p with p > 1 can be characterized in
terms of uniform Lp bounds on h,i u; notice that one implication was already established
in (1.18).
Lemma 4.10. Consider u 2 Lploc (⌦), with 1 < p 1 and fix i 2 {1, . . . , n}. The partial
@u
derivative @x i
belongs to Lploc (⌦) if and only if
Z
8⌦ b ⌦
0
9 c(⌦ )0
s.t. ( h,i u) ' c(⌦0 )k'kLp0 (⌦0 ) 8 ' 2 Cc1 (⌦0 ) ,
⌦0
39
Proof. The first implication has been proved in (1.18), because we know that h,i u is
bounded in Lploc (⌦) when h ! 0, so we can conclude with Hölder’s inequality.
Now fix ⌦0 b ⌦,
Z Z Z
@'
u dx = lim u h,i ' dx = lim ( h,i u) ' dx c(⌦0 )k'kLp0 (⌦0 ) ;
⌦0 @xi h!0 ⌦0 h!0 ⌦0
0
because of the duality relation between Lp (⌦0 ) and Lp (⌦0 ), the weak derivative @xi u exists
and belongs to Lploc (⌦). ⇤
Let us see how Lemma 4.10 contributes to regularity theory, still in the simplified case
1
of the Poisson equation. Suppose f 2 Hloc (⌦) in the Poisson equation, then translation
invariance and linearity allow us to write
Thanks to Lemma 4.10, h,i f is bounded in L2loc (⌦), then by the Caccioppoli-Leray
inequality |r h,i u| is bounded in L2loc (⌦).
As h,i (ru) = r h,i u is bounded in L2loc (⌦; Rn ), thanks to Lemma 4.10 again (applied
componentwise) we get
@
(ru) 2 L2loc (⌦; Rn ) .
@xi
After these preliminaries about Nirenberg’s method, we are now ready to prove the
main result concerning H 2 regularity.
1
for all u 2 Hloc (⌦; Rm ) weak solution of the equation
div(Aru) = f div(F )
In order to simplify the notation, in the following proof let s denote the unit vector
corresponding to a given fixed direction and consequently ⌧h := ⌧h,s and h := h,s .
40
Remark 4.12. Although the thesis concerns a generic domain ⌦0 b ⌦, it is enough to
prove it for balls inside ⌦. More precisely, if 2R < dist(⌦0 , @⌦), we just need to prove the
inequality
Z ⇢Z Z
2 2 2
⇥ 2 ⇤
|r u| dx c |u| dx + |f | + |rF |2 dx
BR/2 (x0 ) B2R (x0 ) B2R (x0 )
for any x0 2 ⌦0 . The general result can be easily obtained by a compactness and covering
argument.
Notice also that the statement as given is redundant, since the term div(F ) can always
be absorbed into f . We will see however that the optimal estimate is obtained precisely
doing the opposite, i.e. considering heuristically f as a divergence.
Proof. We assume x0 = 0 and, by the previous remark, F = 0 (possibly changing f ). In
addition, we prove the result under the stronger assumption that the Legendre condition
with constant holds uniformly in ⌦.
First note that the given equation is equivalent, by definition, to the identity
Z Z
Arur' dx = f ' dx (4.12)
⌦ ⌦
for all ' 2 Cc1 (⌦; Rm ). If we apply it to the test function ⌧ h ' with |h| ⌧ 1 and we do
a change of variable, we find
Z Z
⌧h (Aru)r' dx = ⌧h f ' dx . (4.13)
⌦ ⌦
Subtracting (4.12) to equation (4.13) and dividing by h, we get (thanks to the discrete
Leibniz property)
Z Z Z
(⌧h A)r( h u)r' dx = ( h f )' dx ( h A)rur' dx ,
⌦ ⌦ ⌦
41
To this aim, take a cut-o↵ function ⌘ compactly supported in BR , with 0 ⌘ 1,
identically equal to 1 on BR/2 and such that |r⌘| 4/R, and insert in (4.14) the test
function := ⌘ 2 h u = ⌘ 2 v with |h| < R/2.
Using Young inequality as in Theorem 4.1 (see (4.6)), we get
Z Z
3 2 2 4⇤"
⌘ |rv| ⌘ 2 |rv|2
4 BR R BR
✓ ◆Z Z ✓ ◆Z
4⇤ 4 2 2 1
+ + |v| + ⌘ v hf + +4 |G|2 ,
R" R2 BR BR BR
with ⇤ depending only on A. As in the proof of Theorem 4.1, we absorb the term with
k⌘rvk2L2 (BR ) in the left side of the inequality, so that, up to some constant c > 0 depending
on ( , ⇤, R), we get
Z Z Z Z
2 2 2
c ⌘ |rv| dx |v| dx + |G| dx + ⌘ 2 v h f dx .
2
(4.15)
B2R BR BR
by means of (1.18). The right hand side can in turn be estimated using the classical
Caccioppoli-Leray inequality for u between the balls B3R/2 and B2R : it gives an upper
bound of the desired form.
R
Concerning the term ⌘ 2 v h f dx, by means of discrete integration by parts and Young
inequality, we can write
Z Z Z
2 2 2 1
⌘ v h f dx "˜ | h (⌘ v)| dx + |f |2 dx . (4.16)
BR BR "˜ BR
The first term in the right hand side of (4.16) can be estimated with (since |r⌘|2 64/R2 )
Z Z Z
2 2 4 2 128
|r(⌘ v)| dx 2 ⌘ |rv| dx + 2 |v|2 dx ,
BR+h BR+h R BR+h
4 2
R ⌘ ⌘2 we can absorb the
so that choosing " sufficiently small and using the inequality
first term and use once more the CL inequality to estimate BR+h |v| dx.
The term involving the integral kGk2L2 (BR ) can be estimated in the very same way, using
this time also the local Lipschitz assumption on A to bound h A, so that finally we put
together all the corresponding estimates to obtain the thesis (the conclusion comes from
Lemma 4.10 and then letting h ! 0 in the estimate involving v = h u). ⇤
42
Remark 4.13. It should be clear from the proof that the previous result only concerns
interior regularity and cannot be used in order to get information about the behaviour
of the function u near the boundary @⌦. In other terms, we can not guarantee that the
constant c remains bounded as ⌦0 invades ⌦ (so that R ! 0), even if global regularity
assumptions on A, u, f and F are made. The issue of boundary regularity requires
di↵erent techniques that will be described later on.
43
for some constant ck > 0. Consequently, thanks to our choice of the integer k, we can
find a constant such that
Z
sup |u|2 dx |u|2 dx .
B2 k (x0 ) B1 (x0 )
where !n denotes the Lebesgue measure of the unit ball in Rn . Hence, for this
R case we have
the thesis, provided c( , ⇤) !n . If r 2 (2 k , 1), then it is clear that Br (x0 ) |u|2 dx
R
B1 (x0 )
|u|2 dx and so, since we have a lower bound for r, we just need to choose c( , ⇤)
such that c( , ⇤) 2kn .
We can now prove the second inequality, that concerns the notion of variance of the
function u on a ball.
Proof of (5.2). Again, it is useful to study two cases separately. If r R/2, then by
the Poincaré inequality there exists a constant c(n) such that
Z Z
2
|u ux0 ,r | dx c(n)r 2
|ru|2 dx
Br (x0 ) Br (x0 )
and so
Z ◆n Z ✓
2 r 2
|u ux0 ,r | dx c(n)r |ru|2 dx
Br (x0 ) R/2 BR/2 (x0 )
✓ ◆n+2 Z
r
c(n, , ⇤) |u uR,x0 |2 dx
R/2 BR (x0 )
respectively by the previous result applied to the gradient ru and finally by the Cacciop-
poli-Leray inequality. For the case R/2 < r R we need to use the following fact, that
will be discussed below: the mean value ux0 ,r is a minimizer for the function
Z
m7 ! |u m|2 dx . (5.3)
Br (x0 )
44
Let us go back to the study of
Z
inf |u m|p dx
m2R ⌦
Since the problem is clearly translation invariant, it is sufficient to prove inequality (5.4)
for m = 0. But in this case
Z Z Z Z
p p 1 p p 1 p p
|u u⌦ | dx 2 |u| dx + 2 |u⌦ | dx 2 |u|p dx ,
⌦ ⌦ ⌦ ⌦
thanks to the elementary inequality |a + b|p 2p 1 |a|p + |b|p and to the fact that
Z Z
p
|u⌦ | dx |u|p dx
⌦ ⌦
45
rectangle R0 = ( a/2, a/2)n 1 ⇥ (0, a/2) is not relatively compact in R, nevertheless via
Nirenberg’s method we may find estimates of the form
Z Z Z
2 c 2
|@xs ru| dx 2 |ru| dx + c |f |2
R 0 a R R
• f 2 L2 (⌦; Rm );
• F 2 H 1 (⌦; Rm⇥n );
2 ⇥n2
• A 2 C 0,1 (⌦; Rm );
• ⌦ has a C 2 boundary, in the sense that it is, up to a rigid motion, locally the
subgraph of a C 2 function.
Theorem 6.1. Under the previous assumptions, the function u belongs to H 2 (⌦; Rm ) and
⇥ ⇤
kukH 2 c(⌦, A, n) kf k2 + kF kH 1
46
Proof. Since we already have the interior regularity result at our disposal, it suffices to
show that for any x0 2 @⌦ there exists a neighbourhood U of x0 in ⌦ such that u 2 H 2 (U ).
Without loss of generality we assume x0 = 0. We consider first the case of a flat boundary.
Step 1. (Flat boundary) By applying Nirenberg’s method as described above for the
case of the constant coefficient operator we get @x↵ ui 2 H 1 (R0 ) for ↵ = 1, 2, . . . , n 1
and i = 1, 2, . . . , m, and
✓ ✓ ◆◆ ✓ ◆ ✓ ◆
@u @f @F @A
div Ar = + div + div ru . (6.2)
@x↵ @x↵ @x↵ @x↵
Anyway, we cannot include in the previous conclusion the second derivatives @x2n xn ui and
here we really need to refine a bit the strategy seen above for the Poisson equation. Ac-
tually, this is not complicated because the equation readily implies that @xn (Ann j
ij @xn u ) 2
L2 (R0 ) for any i 2 {1, 2, . . . , m}. Formally this implies, by the Leibniz rule, that Ann 2
ij @xn xn u
j
2 0
belong to L (R ); this is formal because one of the factors is only a distribution (not yet
a function). To make this rigorous, we use the di↵erence quotients in the xn direction
and the discrete Lebniz rule: since by Lemma 4.10 the di↵erence quotients h (Ann j
ij @xn u )
have uniformly bounded L2 norm in Rh0 = {x 2 R : dist(x, @R0 ) > h}, we obtain that the
same is true for Ann ij
j nn
h @xn u . Since the matrix Aij is invertible with detAij
nn m
(as a
consequence of the Legendre-Hadamard condition) we get
Z
lim sup | h @xn uj |2 dx < 1
h!0+ 0
Rh
47
where of course the boundary condition has to be interpreted in the weak sense and
⇥ ⇤
fe = f G, Fe = (F · DH) G, A e = DH · A · (DH)t G
(here contractions are understood with respect to the greek indices, the only ones involved
in the change of variables, see (6.4) below). These formulas can be easily derived by an
elementary computation, starting from the weak formulation of the problem and apply-
ing a change of variables in order to express the di↵erent integrals in terms of the new
coordinates. For instance
Z Z
i
fi (x)' (x) dx = fi G(y)'i G(y) det(rG(y)) dy
⌦0 R
just letting x = G(y), but then det(rG) = 1 and we can set ' = H so that equivalently
= ' G and Z Z
i
fi (x)' (x) dx = fei (y) i (y) dy .
⌦0 R
The computation for Fe or A e is less trivial, but there is no conceptual difficulty. We just
see the first one:
Z Z
↵ @'i @'i
Fi (x) (x) dx = Fi↵ (G(y)) (G(y)) det(rG(y)) dy
⌦0 @x↵ R @x↵
Z
@ i @H
= Fi↵ (G(y)) (y) (G(y)) dy
R @y @x↵
which leads to the conclusion. Note that here and above the arbitrary test function '
has been replaced by the arbitrary test function . However, we should ask whether the
conditions on A (for instance, the Legendre-Hadamard condition) still hold true for A.e
e above. In
This is the case and we can verify it directly by means of the expression of A
fact, ✓ 0 0◆
e↵0 0 @H ↵ ↵ @H
Aij = A G (6.4)
@x↵ ij @x
a 2 Rn and b 2 Rm
and so, for any e
✓ 0 ◆✓ 0 ◆
e↵0 0 (y)e @H ↵ @H
Aij a 0 bi bj = A↵ij (G(y))
a↵ 0 e (G(y))e
a↵ 0 (G(y))e
a 0 bi bj
@x↵ @x
2
a|2 |b|2
|rH(G(y))e (rH(G(y))) 1
a|2 |b|2
|e
since clearly
1 2
a|2 (rH(G(y)))
|e a| 2 .
|rH(G(y))e
48
Hence, A e satisfies the Legendre-Hadamard condition for an appropriate constant 0 > 0
depending on and H, and of course A e 2 C 0,1 (R).
Through this transformation of the domain, we can finally apply Step 1 and find that
v 2 H 2 (R0 ). Coming back to the original variables we obtain the H 2 regularity of u.
⇤
If both the boundary and the data are sufficiently regular, this method can be iterated
to get the following theorem.
Theorem 6.2. Assume, in addition to the hypotheses above, that f 2 H k (⌦; Rm ) and
2 2
also F 2 H k+1 (⌦; Rm⇥n ), A 2 C k,1 (⌦, Rm ⇥n ) with ⌦ such that @⌦ 2 C k+2 . Then u 2
H k+2 (⌦; Rm ).
We are not going to present the detailed proof of the previous result, but the basic
idea consists in di↵erentiating the starting equation with respect to each fixed direction
to get an equation having the form of (6.3), as in (6.2), provided we set Fe = @x
@F
↵
@A
+ @x ↵
ru.
49
Recall that u is a local minimizer for I if, for any
Z Z
u 2 Hloc (⌦; R ), spt(u u ) ⇢ ⌦ b ⌦ =)
1 m 0 0
F (ru0 ) dx F (ru) dx .
⌦0 ⌦0
If this is the case, we have already seen how the Euler-Lagrange equation can be ob-
tained: considering perturbations of the form u0 = u + tr' with ' 2 Cc1 (⌦, Rm ) one can
prove (using the fact that the regularity assumptions on F allow di↵erentiation under the
integral sign) that
Z Z
d @'i
0= F (ru + tr') dx = Bi↵ (ru) dx .
dt ⌦ t=0 ⌦ @x↵
Now, suppose s is a fixed coordinate direction (and let es be the corresponding unit vector)
and h > 0 a small positive increment: if we apply the previous argument to a test function
having the form ⌧ h ', we get
Z
@'i
⌧h (Bi↵ (ru)) dx = 0
⌦ @x↵
and consequently, subtracting this to the previous one
Z
↵ @'i
h,s (B i (ru)) dx = 0 .
⌦ @x↵
However, as a consequence of the regularity of F , we can write
Z 1
↵ ↵ d ↵
Bi (ru(x + hes )) Bi (ru(x)) = Bi (tru(x + hes ) + (1 t)ru(x)) dt
0 dt
Z 1 j
↵ @u @uj
= Aij (tru(x + hes ) + (1 t)ru(x)) dt (x + hes ) (x)
0 @x @x
and setting Z 1
e↵ (x)
A := A↵ij (tru(x + hes ) + (1 t)ru(x)) dt
ij,h
0
we rewrite the previous condition as
Z j
e↵ (x) @ h,s u @'i
A ij,h (x) (x) dx = 0 .
⌦ @x @x↵
eh rw) = 0 .
div(A (7.1)
50
It is obvious by the definition that A e↵ (x) satisfies both the Legendre condition for
ij,h
the given constant > 0 and a uniform upper bound on the L1 norm. Therefore we can
apply the Caccioppoli-Leray inequality to the problem (7.1) to obtain constants C1 and
C2 , not depending on h, such that
Z Z
2 C1
|r( h,s u)| dx 2 | h,s u|2 dx C2
BR (x0 ) R B2R (x0 )
for any BR (x0 ) ⇢ B2R (x0 ) b ⌦. Consequently, by Lemma 4.10 we deduce that
2
u 2 Hloc (⌦; Rm ). (7.2)
51
Definition 8.1 (Hölder spaces). Given A ⇢ Rn , u : A ! Rm and ↵ 2 (0, 1] we define
the ↵-Hölder semi-norm on A as
|u(x) u(y)|
kuk↵,A := sup .
x6=y2A |x y|↵
Remark 8.2. The spaces C k,↵ (⌦; Rm ) are Banach when endowed with the norm
X
kukC k,↵ = ri u C 0,↵ .
|i|k
Remark 8.4. We mention here some of the basic properties of the Morrey spaces Lp, :
(i) Lp, (⌦; R) are Banach spaces, for any 1 p < 1 and 0;
52
Note that the condition (n )/p (n µ)/q can also be expressed by asking c
with the critical value c defined by the equation (n c )/p = (n µ)/q. The proof of the
first result is standard, the second statement is trivial, while the third and fourth ones are
immediate applications of Lebesgue Di↵erentiation Theorem. Finally the last one relies
on Hölder inequality:
✓Z ◆ ✓Z ◆p/q
p q
|f | dx |f | dx (!n rn )(1 p/q)
⌦(x,r) ⌦(x,r)
The mean fx0 ,r defined in (8.2) might not be optimal in the calculation of the sort of
p-variance in (8.1), anyway it gives equivalent results, thanks to (5.4).
Remark 8.6. As in Remark 8.4, we briefly highlight the main properties of Campanato
spaces.
(i) As defined in (8.1), k · kLp, is merely a seminorm because constants have null Lp,
norm. If ⌦ is connected, then Lp, modulo constants is a Banach space.
We will see that a converse statement holds (namely functions in these Campanato
spaces have Hölder continuous representatives in their Lebesgue equivalence class),
and this is very useful: we can replace the pointwise definition of Hölder spaces with
an integral one.
53
Actually, Campanato spaces are interesting only when n, exactly because of
their relationship with Hölder spaces. On the contrary, if < n, Morrey spaces and
Campanato spaces are basically equivalent. In the proof of this and other results we need
a mild regularity assumption on ⌦, namely the existence of c⇤ > 0 satisfying
For instance this assumption includes domains which are locally subgraphs of Lipschitz
functions, while it rules out domains with outer cusps.
Theorem 8.7. Let ⌦ ⇢ Rn be an open bounded set satisfying (8.3) and let 0 < n.
Then the spaces Lp, and Lp, are equivalent, i.e.
2 p 1
kf kpLp, r +⇢ 2 p
kf kpLp, ⇢ ,
54
thus we obtained that
n
⇣ ⇢ ⌘ np n
|fx0 ,r fx0 ,⇢ | ckf kLp, r p ⇢ p = ckf kLp, ⇢ p . (8.4)
r
(k+1)
Now fix a radius R > 0: if r = 2 R and ⇢ = 2 k R, inequality (8.4) means that
✓ ◆ n
R p
|fx0 ,R/2k+1 fx0 ,R/2k | ckf kLp, , (8.5)
2k
and, adding up when k = 0, . . . , N 1, it means that
n ✓ ◆ n
n 2N p 1 R p
|fx0 ,R/2N fx0 ,R | ckf kLp, R p
n ckf kLp, . (8.6)
2 p 1 2N
Let us go back to our purpose of estimating |fx0 ,r |p : we choose R 2 (d⌦ /2, d⌦ ) and
N 2 N such that r = R/2N . By triangular inequality
|fx0 ,r |p 2p 1
(|fx0 ,r fx0 ,R |p + |fx0 ,R |p ) ;
since
|fx0 ,R | c(d⌦ )kf kLp ,
the only thing left to conclude is to apply inequality (8.6) in this case:
|fx0 ,r fx0 ,R |p ckf kpLp, r n
,
that is all we needed to conclude that
Z
r |f |p c kf kpLp, + dn⌦ kf kpLp .
⌦(x0 ,r)
Remark 8.8. When the dimension of the domain space is n, the Campanato space L1,n is
very important in harmonic analysis and elliptic regularity theory: after John-Nirenberg
seminal paper, this space is called BM O (bounded mean oscillation). It consists of the
space of all functions f : ⌦ ! R such that there exists a constant C satisfying the
inequality Z
|f (x) fx0 ,r | dx Crn 8 r 2 (0, d⌦ ), 8 x0 2 ⌦ .
⌦(x0 ,r)
Notice that L (⌦) ( BM O(⌦): for example, consider ⌦ = (0, 1) and f (x) = ln x. For
1
55
R a+r
hence ln x 2 BM O(⌦). For simplicity, we replaced the mean a ln s ds with ln(a + r),
but, up to a multiplicative factor 2, this does not make a di↵erence. On the contrary
/ L1 (⌦).
ln x 2
Theorem 8.9 (Campanato). With the previous notation, when n < n+p Campanato
spaces Lp, are equivalent to Hölder spaces C 0,↵ with ↵ = ( n)/p. Moreover, if ⌦ is
connected and > n + p, then Lp, is equivalent to the set of constants.
Proof. As in the proof of Theorem 8.7, the letter c denotes a generic constant depending
on the exponents, the space dimension n and the constant c⇤ in (8.3).
Let = n + ↵p. We already observed in Remark 8.6 that C 0,↵ ⇢ Lp, , so we need to
prove the converse inclusion: given a function f 2 Lp, , we are looking for a representative
in the Lebesgue equivalence class of f which belongs to C 0,↵ .
Recalling inequality (8.5) with fixed radius R > 0 and x 2 ⌦, we obtain that the
sequence (fx,R/2k ) has the Cauchy property. Hence we define
Z
˜
f (x) := lim f (y) dy .
k!1 ⌦(x,R/2k )
Clearly
Z Z
|f (y) p
fx,R/2k | dy ! 0 =) |f (y) f˜(x)|p dy ! 0 , (8.7)
⌦(x,R/2k ) ⌦(x,R/2k )
In particular, f˜ does not depend on the chosen initial radius R. Let us prove that
f˜ 2 C 0,↵ (⌦) .
We employ again an inequality from the proof of Theorem 8.7: letting N ! 1 in (8.6),
we get that
|f˜(x) fx,R | ckf kLp, R↵
with ↵ = ( n)/p; consequently, given x, y 2 ⌦ and choosing R = 2|x y|,
|f˜(x) f˜(y)| |f˜(x) fx,R | + |fx,R fy,R | + |fy,R f˜(y)| c|x y|↵ + |fx,R fy,R | .
56
The theorem will be proved if we can estimate |fx,R fy,R |. To this aim, we use the
inclusion ⌦(y, R/2) ⇢ ⌦(x, R) to get
Z
n n p
c⇤ 2 R |fx,R fy,R | |fx,R fy,R |p ds
⌦(y,R/2)
✓Z Z ◆
p 1 p p
2 |f (s) fx,R | ds + |f (s) fy,R |
⌦(x,R) ⌦(y,R)
p p
2 kf kLp, R ,
and finally
n
|fx,R fy,R | ckf kLp, R p c|x y|↵ .
⇤
The following inclusions follow by the Hölder and the Poincaré inequalities, respec-
tively.
Proposition 8.10 (Inclusions between Lebesgue and Morrey spaces, Morrey and Cam-
1,n/p0
panato spaces). For all p 2 (1, 1), Lploc (⌦) ⇢ Lloc (⌦). In addition
|ru| 2 Lp,
loc (⌦) =) u 2 Lp,
loc
+p
(⌦) . (8.8)
0,↵
Corollary 8.11 (Sobolev embedding for p > n). If p > n, then W 1,p (⌦) ⇢ Cloc (⌦), with
1,p 0,↵
↵ = 1 n/p. If ⌦ is bounded and regular, then W (⌦) ⇢ C (⌦).
Proof. By the previous proposition we get
1,p 1,n/p0
u 2 Wloc =) |ru| 2 Lloc (⌦) = L1,n n/p
(⌦) = L1,n 1+↵
(⌦) . (8.9)
57
We have seen that u satisfies the Euler-Lagrange equations, for (9.1) they are
@
Fp↵i (ru) = 0 i = 1, . . . , m . (9.2)
@x↵
We have also seen in Section 7 how, di↵erentiating (9.2) along the direction xs , one can
obtain ✓ ◆
@ @ 2 uj
Fp↵ p (ru) =0 i = 1, . . . , m . (9.3)
@x↵ i j @x @xs
In the spirit of Hilbert’s XIX problem, we are interested in the regularity properties
of u. Fix s 2 {1, . . . , n}, let us call
@u
w(x) := (x) 2 L2 (⌦, Rm ) ,
@xs
A(x) := r2 F (ru(x)) ,
thus (9.3) can be written as
✓ ◆
@ @ 2 uj
div (Arw) = Fp↵ p (ru) =0. (9.4)
@x↵ i j @x @xs
1
Since w 2 Hloc (⌦; Rm ) by (7.2), we can use the Caccioppoli-Leray inequality for w, in
the sharp version of Remark 4.4. Combining it with the Poincaré inequality (choosing k
equal to the mean value of w on the ball BR (x0 ) \ BR/2 (x0 )), we obtain
Z Z Z
2 2 2
|rw| dx cR |w k| dx c |rw|2 dx ,
BR/2 (x0 ) BR (x0 )\BR/2 (x0 ) BR (x0 )\BR/2 (x0 )
R
thus, adding c |rw|2 dx to both sides, we get
BR/2 (x0 )
Z Z
2 c
|rw| dx |rw|2 dx .
BR/2 (x0 ) c + 1 BR (x0 )
Now, if ✓ := c/c + 1 < 1 and ↵ = log2 ✓, we can write the previous inequality as
Z ✓ ◆↵ Z
2 1
|rw| dx |rw|2 dx . (9.5)
BR/2 (x0 ) 2 BR (x0 )
In order to get a power decay inequality from (9.5), we state this basic iteration lemma.
Lemma 9.1. Consider a non-decreasing function f : (0, R0 ] ! R satisfying
⇣ ⇢ ⌘ ✓ 1 ◆↵
f f (⇢) 8 ⇢ R0 .
2 2
Then ⇣ r ⌘↵
↵
f (r) 2 f (R) 8 0 < r R R0 .
R
58
Proof. Fix r < R R0 and choose a number N 2 N such that
R R
<r .
2N +1 2N
It is clear from the iteration of the hypothesis that
✓ ◆ ✓ ◆↵N
R 1
f f (R) ,
2N 2
thus, by monotonicity,
N ↵N
f (r) f 2 R 2 f (R) = 2↵ 2 ↵(N +1)
f (R) < 2↵ (r/R)↵ f (R) .
⇤
Thanks to Lemma 9.1, we are ready to transform (9.5) in
Z ⇣ ⇢ ⌘↵ Z
2
|rw| dx c |rw|2 dx 80 < ⇢ R ,
B⇢ (x0 ) R BR (x0 )
As long as F is sufficiently regular, the iteration of this argument solves XIX Hilbert’s
regularity problem in the C 1 category.
We close this section with a more technical but useful iteration lemma in the same
spirit of Lemma 9.1.
59
If
✓ ◆ ↵↵
1
" (9.7)
2A
for some 2 ( , ↵), then
h⇣ ⇢ ⌘ i
f (⇢) c(↵, , , A) 8 0 < ⇢ R R0 .
f (R) + B⇢ (9.8)
R
Proof. Without loss of generality, we assume A > 1/2. We choose ⌧ 2 (0, 1) such that
2A⌧ ↵ = ⌧ , (9.9)
thus (9.7) gives the inequality
" ⌧↵ . (9.10)
The following basic estimate uses the hypothesis (9.6) jointly with (9.9) and (9.10):
f (⌧ R) A(⌧ ↵ + ")f (R) + BR
2A⌧ ↵ f (R) + BR = ⌧ f (R) + BR . (9.11)
The iteration of (9.11) easily gives
f (⌧ 2 R) ⌧ f (⌧ R) + B⌧ R ⌧ 2 f (R) + ⌧ BR + B⌧ R
= ⌧ 2 f (R) + BR ⌧ (1 + ⌧ ).
It now can be easily proven by induction that
N
X1
N N (N 1) 1 ⌧ N( )
f (⌧ R) ⌧ f (R) + BR ⌧ ⌧ k( )
= ⌧ N f (R) + BR ⌧ (N 1)
.
k=0
1 ⌧( )
60
Remark 9.3. The fundamental gain in Lemma 9.2 is the passage from R to ⇢ and the
removal of ", provided that " is small enough. These improvements can be obtained at
the price of passing from the power ↵ to the worse power < ↵.
10 Schauder theory
We are treating Schauder theory in a local form in ⌦ ⇢ Rn , just because it would be
too long and technical to deal also with boundary regularity (some ideas are analogous
to those used in Section 6). We shall describe first a model result for constant coefficient
operators, and then we will consider the case of Hölder continuous coefficients.
We recall the usual PDE we are studying, in a divergence form:
8
< div (Aru) = divF in ⌦ ;
(10.1)
: 1
u 2 Hloc (⌦; Rm ) .
Theorem 10.1. If A↵ij are constant and satisfy the Legendre-Hadamard condition for
some > 0, then for all µ < n + 2 it holds
F 2 L2,µ
loc (⌦) =) ru 2 L2,µ
loc (⌦) .
Proof. In this proof, c = c(n, , |A|) and its value can change from line to line. Since
the estimates we make are local, we assume with no loss of generality that F 2 L2,µ (⌦).
Let us fix a ball BR b ⌦ with center x0 2 ⌦ and compare with u the solution v of the
homogeneous problem 8
< div(Arv) = 0 in BR ;
(10.2)
:
v=u in @BR .
Since rv belongs to H 1 for previous results concerning H 2 regularity and its components
@v
@x↵
solve the same problem (because we supposed to have constant coefficients), we can
use the decay estimates (5.1) and (5.2).
So, if 0 < ⇢ < R, (5.2) provides us with the following inequality:
Z ⇣ ⇢ ⌘n+2 Z
2
|rv(x) (rv)⇢ | dx c |rv(x) (rv)R |2 dx . (10.3)
B⇢ R BR
Now we try to employ (10.3) to get some estimate for u, the original “non-homogeneous”,
solution of (10.1). Obviously, we can write
u=w+v ,
61
where w 2 H01 (BR ; Rm ). Thus (first using ru = rv + rw, then the minimality of the
mean and (10.3), eventually rv = ru rw and (rw)R = 0)
Z
|ru(x) (ru)⇢ |2 dx
B⇢
Z Z !
2 |rw(x) (rw)⇢ |2 dx + |rv(x) (rv)⇢ |2 dx
B⇢ B⇢
Z ⇣ ⇢ ⌘n+2Z
2 |rw(x) (rw)R |2 dx + c |rv(x) (rv)R |2 dx
B R BR
Z ⇢ ⇣ ⇢ ⌘n+2Z
c |rw(x)|2 dx + c |ru(x) (ru)R |2 dx .
BR R BR
is non decreasingR because of the minimality property of the mean (ru)⇢ , when one min-
imizes m 7! B⇢ |ru(x) m|2 dx. In order to get that f satisfies the hypothesis of
R
Lemma 9.2, we have to estimate BR |rw|2 dx. We can consider w as a function in
H 1 (Rn ) (null out of ⌦) so, by Gårding inequality (choosing the test function ' = w),
Z Z
2
|rw(x)| dx c Arw(x)rw(x) dx
BR BR
Z Z
= c F (x)rw(x) dx = c (F (x) FR )rw(x) dx (10.4)
BR
because div(Arw)
R = divF by linearity. Applying Young inequality to (10.4) and then
absorbing BR |rw|2 dx in the left side of (10.4), we get
Z Z
2
|rw(x)| dx c |F (x) FR |2 dx ckF k2L2,µ Rµ ,
BR BR
because F 2 L2,µ .
Therefore we obtained the decay inequality of Lemma 9.2 for f with ↵ = n + 2, =µ
and " = 0, then ⇣ ⇢ ⌘µ
f (⇢) c f (R) + c⇢µ ,
R
that is ru 2 L2,µ . ⇤
Corollary 10.2. With the previous notation, when µ = n + 2↵, Theorem 10.1 and Cam-
panato Theorem 8.9 yield that
F 2 C 0,↵ =) ru 2 C 0,↵ .
62
In the next theorem we consider the case of variable, but continuous, coefficients,
proving in this case a Lp,µ regularity of |ru| with µ < n; as we have seen, the Poincaré
inequality then provides Hölder regularity at least of u if µ + p > n.
Theorem 10.3. Considering again (10.1), suppose that A↵ij 2 C(⌦) and A satisfies a
(locally) uniform Legendre-Hadamard condition for some > 0. If F 2 L2,µ
loc with µ < n,
then |ru| 2 L2,µ
loc .
Naturally, since µ < n, Campanato spaces and Morrey spaces coincide, so that we
used Morrey spaces for simplicity.
Proof. Here is an example of Korn’s technique of freezing of coefficients. We use the
same convention on c of the previous proof, namely c = c(n, , sup |A|).
Fix a point x0 2 ⌦ and define
F̃ (x) := F (x) + (A(x0 ) A(x)) ru(x) ,
so that the solution u of (10.1) solves
div(A(x0 )ru(x)) = divF̃ (x) with F̃ (x) := F (x) + (A(x0 ) A(x))ru(x) .
Write u = v + w, where v solves the homogeneous PDE (10.2) with frozen coefficients
A(x0 ). Using (5.1) for v we obtain
Z ⇣ ⇢ ⌘n Z Z
2 2
|ru(x)| dx c |rv(x)| dx + c |rw(x)|2 dx
B⇢ R BR B
⇣ ⇢ ⌘n Z Z R
c |rv(x)|2 dx + c |F̃ (x) FR |2 dx .
R BR BR
Consequently, as F 2 L2,µ
loc ,
Z Z
2 µ 2
|F̃ (x) FR | dx c̃R + 2! (R) |ru(x)|2 dx
BR BR
with c̃ depending only on kF kL2,µ . We are ready to use Lemma 9.2 with f (⇢) :=
R 2
loc
B⇢
|ru(x)| dx, ↵ = n, = µ < n and " = ! 2 (R): it tells us that if R is under a
threshold depending only on c, ↵, , ! and kF kL2,µ we have
loc
⇣ ⇢ ⌘µ
f (⇢) c f (R) + c⇢µ ,
R
so that |ru| 2 L2,µ
loc . ⇤
63
We can now prove Schauder theorem for elliptic PDE’s in divergence form. In the
non-divergence form the result is (in the scalar case)
X @ 2u
A↵ 2 C 0,↵ =) u 2 C 2,↵ , (10.6)
↵,
@x ↵ @x
if A is of class C 0,↵ . The proof follows similar lines, i.e. starting for second derivative decay
estimates for constant coefficient operators, and then freezing the coefficients. Notice also
that both (10.6) and Theorem 10.4 below are easily seen to be optimal, considering 1-
dimensional ODE’s au00 = f or (au0 )0 = f 0 .
Theorem 10.4 (Schauder). Suppose that the coefficients A↵ij (x) of the PDE (10.1) belong
to C 0,↵ (⌦) and A satisfies a (locally) uniform Legendre-Hadamard in ⌦ for some > 0.
Then the following implication holds
0,↵ 0,↵
F 2 Cloc =) ru 2 Cloc ,
that is to say
F 2 L2,n+2↵
loc =) ru 2 L2,n+2↵
loc .
Proof. With the same idea of freezing coefficients (and the same notation, too), we
estimate by (5.1)
Z ⇣ ⇢ ⌘n+2 Z Z
2 2
|ru(x) (ru)⇢ | dx c |ru(x) (ru)R | dx + c |F̃ (x) FR |2 dx .
B⇢ R BR BR
(10.7)
Additionally, the Hölder property of A makes us rewrite (10.5) as
Z Z Z
2 2 2↵
|F̃ (x) FR | dx 2 |F (x) FR | dx + cR |ru(x)|2 dx . (10.8)
BR BR BR
0,↵
Since F 2 Cloc , we obtain
Z Z
2 n+2↵ 2↵
|F̃ (x) FR | dx cR + cR |ru(x)|2 dx .
BR BR
Adding (10.9) to (10.7) and applying Lemma 9.2 with exponents n + 2 and n + ↵, we
get ru 2 L2,n+↵ , so that ru 2 C 0,↵/2 , in particular |ru| is locally bounded. Using this
information we can improve (10.9) as follows:
Z
|F̃ (x) FR |2 dx cRn+2↵ .
BR
Now we reach the conclusion, again by Lemma 9.2 with exponents n + 2 and n + 2↵.
⇤
64
11 Regularity in Lp spaces
In this section we deal with elliptic regularity in the category of Lp spaces, obviously a
natural class of spaces besides Morrey, Hölder and Campanato spaces.
65
The Markov inequality inspires the definition of a space which is weaker than Lp , but
still keeps (11.3).
Definition 11.3 (Marcinkiewicz space). Given a measure space (⌦, F, µ) and an exponent
1 p < 1, the Marcinkiewicz space Lpw (⌦, µ) is defined by
q<p =) Lp ⇢ Lpw ⇢ Lq .
The first inclusion is due to Markov inequality (11.2), on the other hand, if f 2 Lpw , then
Z Z 1 ✓Z 1 Z 1 ◆
q q 1 q 1 q 1
|f | dµ(x) = q t F (t) dt q t F (t) dt + t F (t) dt
⌦ 0 0 1
Z 1
q
qµ(⌦) + q tq 1 kf kpLpw t p dt = qµ(⌦) + kf kpLpw .
1 p q
Definition 11.5 (Maximal operator). When f 2 L1loc (Rn ) we define the maximal function
Mf by Z
Mf (x) := sup |f (y)| dy , (11.5)
Qr (x) Qr (x)
where Qr (x) is the n-dimensional cube with center x and side length r.
1
Mf (x) = when |x| 1,
2|x|
66
However, if f 2 L1 , the maximal operator Mf belongs to the weaker Marcinkiewicz
space L1w , as we are going to see in Theorem 11.8. We first recall the Vitali covering
theorem, in a version valid in any metric space.
Lemma 11.7 (Vitali). Let F be a finite family of balls in a metric space (X, d). Then,
there exists G ⇢ F, made of disjoint balls, satisfying
[ [
B⇢ B̂ .
B2F B2G
Here, for B ball, B̂ denotes the ball with the same center and triple radius.
Proof. The initial remark is that if B1 and B2 are intersecting balls then B1 ⇢ B c2 ,
provided the radius of B2 is larger than the radius of B1 . Assume that the family of balls
is ordered in such a way that their radii are non-increasing. Pick the first ball B1 , then
pick the first ball among those that do not intersect B1 and continue in this way, until
either there is no ball left or all the balls left intersect one of the chosen balls. The family
G of chosen balls is, by construction, disjoint. If B 2 F \ G, then B has not been chosen
because it intersects one of the balls in G; the first of these balls Bf has radius larger
than the radius of B (otherwise B would have been chosen before Bf ), hence B ⇢ B cf .
⇤
67
then Vitali’s lemma stated for the distance induced by the sup norm in Rn allows us to
find J ⇢ I such that the cubes Qr(xj ) (xj ), j 2 J, are pairwise disjoint and
[ [
Q3r(xj ) (xj ) Qri (xi ) K .
j2J i2I
We conclude that
X Z
n n n 3n X 3n
L (K) 3 (r(xj )) f (y) dy kf kL1 .
j2J
t i2I Qr(x ) (xi ) t
i
68
independent of u.
If s < 1, we say that T is of weak type (s, s) if
kukss
µ ({x : |T u(x)| > ↵}) C s 8↵ > 0, u 2 D
↵
for some constant C independent of u and ↵. Finally, by convention, T is called of weak
type (1, 1) if it is of strong type (1, 1).
We can derive an appropriate interpolation theorem even in the case of weak continuity.
Theorem 12.3 (Marcinkiewicz Interpolation Theorem). Assume that p, q 2 [1, 1] with
p < q, D ⇢ Lp (X, µ) \ Lq (X, µ) is a linear space stable under truncations and T : D !
Lp (X, µ) \ Lq (X, µ) is Q-subadditive, of weak type (p, p) and of weak type (q, q).
Then T is of strong type (r, r) for all r 2 (p, q).
Remark 12.4. The most important application of the previous result is perhaps the
study of the boundedness of maximal operators (see the next Remark). In that case, one
typically works with p = 1 and q = 1 and we limit ourselves to prove the theorem under
this additional hypothesis.
Proof. We can truncate f 2 D as follows:
f = g + h, g(x) = f (x) {|f | s} (x), h(x) = f (x) {|f |> s} (x) ,
where A1 is the constant appearing in the weak (1, 1) estimate. By integration of the
previous inequality, we get
Z 1 Z 1Z
p 1
p s µ ({|T (f )| > s}) ds 2A1 Qp sp 2 |f | dµ ds
0 0 {|f | s}
69
Remark 12.5 (The limit case p = 1). In the limit case p = 1 we can argue similarly to
find
Z 1
µ ({|T (f )| > s}) ds
1
Z Z |f (x)|/ ! Z
1
2A1 Q s ds |f (x)| dµ(x) = 2A1 Q |f | log |f | dµ.
{|f | } 1 {f }
Proof. Let ⇢ Z
⇤t := x 2 X| lim sup |f (y) f (x)|p dµ(y) > t .
r#0 Br (x)
The thesis can be achieved showing that for any t > 0 we have µ(⇤t ) = 0, since the stated
property holds out of [n ⇤1/n . Now, we can exploit the metric structure of X in order to
approximate f in L1 (µ) norm by means of continuous and bounded functions: for any
" > 0 we can write f = g + h with g 2 Cb (X) and khkpLp t". Hence, it is enough to
prove that for any t > 0 we have µ(At ) = 0 where
⇢ Z
At := x 2 X| lim sup |h(y) h(x)|p dµ(y) > t .
r#0 Br (x)
70
This is easy, because by definition
⇢ ⇢
t t
At ⇢ |h|p > [ M(|h|p ) >
2p+1 2p+1
and, if we consider the corresponding measures, we have (taking Remark 12.6 into account)
2p+1 2p+1
µ(At ) khkpLp + M khkpLp 2p+1 (1 + M )"
t t
where M is the constant in the weak (1, 1) bound. Since " > 0 is arbitrary we get the
thesis. ⇤
Remark 13.2. All the previous results have been derived for the maximal operator
defined in terms of centered balls, that is
Z
M f (x) = sup f (y) dy
r>0 Br (x)
and the Lebesgue di↵erentiation theorem has been stated according to this setting. How-
ever, it is clear that we can generalize everything to any metric space (X, d, µ) with a
finite doubling measure and a suitable family of sets F := [x2X Fx with
Z
MF f (x) = sup f (y) dy
A2Fx A
for a.e. x 2 Rn , Notice that requiring |Q| ! 0 (i.e. diam(Q) ! 0) is essential to “factor”
continuous functions as in the proof of Theorem 13.1.
71
14 Calderón-Zygmund decomposition
We need to introduce another powerful tool, that will be applied to the study of the BM O
spaces. Here and below Q will indicate an open cube in Rn and similarly Q0 or Q00 .
Theorem
R 14.1. Let f 2 L1 (Q), f 0 and consider a positive real number ↵ such that
Q
f dx ↵. Then, there exists a finite or countable family of open cubes {Qi }i2I with
Qi ⇢ Q and sides parallel to the ones of Q, such that
(i) Qi \ Qj = ; if i 6= j;
R
(ii) ↵ < Qi f dx 2n ↵ 8i;
(iii) f ↵ a.e. on Q \ [i Qi .
Remark 14.2. The remarkable (and useful) aspect of this decomposition is that the
“bad” set {f > ↵} is almost all packed inside a family of cubes, carefully chosen in such
a way that still the mean values inside the cubes is of order ↵. As a consequence of the
existence of this decomposition, we have
X XZ
n
↵ L (Qi ) < f dx kf k1 .
i i Qi
At each step we collect the cubes that verify the first condition and put together all
such cubes, thus forming a countable family. The first two properties are obvious
R by
construction: indeed, if Qi is a chosen cube then its parent cube Q̃i satisfies Q̃i f ↵,
R
which gives easily Qi f 2n ↵. For the third one, note that if x 2 Q \ [i Qi , then there
R
exists a sequence of subcubes (Q ej ) with x 2 \j Q
ej and L n Q ej ! 0, e f dx ↵.
Qj
Thanks to the Lebesgue di↵erentiation theorem we get f (x) ↵ for a.e. x 2 Q \ [i Qi .
⇤
72
Remark 14.3 (Again in the limit case p = 1). Using the Calderon-Zygmund decompo-
sition, for ↵ > kf k1 we can reverse somehow the weak (1, 1) estimate:
Z XZ
|f | dx |f | dx 2n ↵L n (Qi ) 2n ↵L n {M |f | > ↵/2n } ,
{|f |>↵} i Qi
because the cubes Qi are contained in {M |f | > ↵/2n }. Using this inequality we can also
reverse theRimplication of Remark 12.5, namely assuming with no loss of generality that
f 0 and f dx = 1:
Z Z 1Z Z 1 Z
1
f log f dx = f dxdt = f dx
{f >1} 0 {log f >t} 1 s {f >s}
Z 1 Z
n s 1 +
2 L {M f > } ds = 2 (M f ) dx .
1 2 2
where uQ0 denotes the mean value of u on Q0 . We also define the seminorm kukBM O
as the supremum in the right hand side. An elementary argument replacing balls with
concentric cubes shows that BM O(Q) ⇠ L1,n , that is the two spaces consist of the same
elements and the corresponding semi-norms are equivalent. Here we recall the inclusion
already discussed in Remark 8.8.
73
However, it should be clear that the previous inclusion is far from being an equality
as elementary examples show, see Remark 8.8. We shall extend now to n-dimensional
spaces the example in Remark 8.8, stating first a simple sufficient (and necessary, as we
will see) condition for BMO.
Proposition 15.2. Let u : Q ! R be a measurable function such that, for some b > 0,
B 0, the following property holds:
Then u 2 BM O(Q).
Example 15.3. Thanks to Proposition 15.2 we can check that ln |x| 2 BM O (0, 1)n .
Indeed, ln |x| satisfies (15.1) (the parameters b and B will be made precise later). To see
this, fix a cube C, with h the length of the side of C. We define, respectively,
so that ✓ ◆
⇠
aC u = ln 0.
|x|
We estimate the Lebesgue measure of C \ {⇠ |x|e }: naturally we can assume that
⇠ ⌘e , otherwise there is nothing to prove, so
p
⇠e ⌘ ⇠ diam(C) ⇠ nh ,
then p
nh
⇠ .
1 e
Finally p
1 n 1 ( n)n !n
L (C \ {|u aC | }) n L n B⇠e e n
,
hn h (1 e )n
so that distinguishing
p n the cases 1 and > 1 we see that (15.1) holds with b = n and
n 1 n
B = max{e , ( n) !n (1 e ) }.
The following theorem by John and Nirenberg was first proved in [21].
74
Theorem 15.4 (John-Nirenberg, first version). There exist constants c1 , c2 depending
only on the dimension n such that
Remark 15.5. In the proof we present here, we will find explicitly c1 = e and c2 =
1/(2n e). However, these constants are not sharp.
Before presenting the proof, we discuss here two very important consequences of this
result.
Corollary 15.6 (Exponential integrability of BM O functions). For any c < c2 there
exists K(c, c1 , c2 ) such that
Z
ec|u uQ |/kukBM O dx K(c, c1 , c2 ) 8u 2 BM O(Q) \ {0} .
Q
where we assumed kukBM O(Q) = 1, L n (Q) = 1 and we used the John-Nirenberg inequal-
ity. ⇤
Remark 15.7 (Better integrability of W 1,n functions). The previous theorem tells that
the class BM O (and hence also W 1,n ) has exponential integrability properties. This result
can be in part refined by the celebrated Moser-Trudinger inequality, that we quote here
without proof.
1/(n 1)
For any n > 1 set ↵n := n!n 1 . and consider a bounded domain ⌦ in Rn , with
n > 1. Then
⇢Z Z
n/(n 1) 1,n
C(⌦) := sup exp ↵n |u| dx : u 2 W0 (⌦), |ru|n dx 1 < 1 .
⌦ ⌦
75
The proof of Theorem 15.8 relies on a simple and standard computation, similar
to the one presented before in order to get exponential integrability. Indeed, assuming
kukBM O = 1, (15.2) gives
Z Z 1 Z 1
p n p 1
|u uQ | dx = p L {|u uQ | > s} s ds c1 p e c2 s sp 1 ds .
Q 0 0
We can now conclude this section, by proving the John-Nirenberg inequality (15.2).
Proof. By homogeneity, we can assume without loss of generality that kukBM O = 1. Let
↵ > 1 be a parameter, to be specified later. We claim that it is possible to define, for any
k 1 a countable family of subcubes Qki i2I contained in Q such that
k
The combination of linear growth in (i) and geometric decay in (ii) leads to the exponential
decay of the repartition function: indeed, choose k such that 2n ↵k t < 2n ↵(k + 1), then
By the first relation c2 = log ↵/(2n ↵) and we maximize with respect to ↵ > 1 to find
1
↵ = e, c1 = e, c2 = .
2n e
Now we just need to prove the claim. If k = 1 we simply apply the Calderón-Zygmund
decomposition to f = |u uQ | for the level ↵ and get a collection {Q1i }i2I1 . We have to
verify that the required conditions are verified. Condition (ii) follows by Remark 14.2,
while (i) is obvious since |u(x) uQ | ↵ a.e. out of the union of Q1i by construction.
But, since kukBM O = 1, we also know that
Z
8i 2 I1 |u uQ1i | dx 1 < ↵ ,
Q1i
76
hence we can iterate the construction, by applying the Calderón-Zygmund decomposition
to each of the functions |u uQ1i | with respect to the corresponding cubes Q1i . In this way,
we find a family of cubes Q2i,l , each contained in one of the previous ones. Moreover
Remark 14.2 and the induction assumption give
X X1Z X1 1
n 2
L (Qi,l ) |u uQ1i | dx L n (Q1i ) 2 L n (Q) ,
i,l i
↵ Q1i i
↵ ↵
so for the first set in the inclusion the thesis is obvious by the case k = 1. For the second
one, we first observe that
Z
|uQ uQ1i | |uQ u| dx 2n ↵
Q1i
With minor changes, we can deal with the general case k > 1 and this is what we need
to conclude the argument and the proof. ⇤
The John-Nirenberg theorem stated in Theorem 15.4 can be extended considering the
p
L norms, so that the case of BM O maps corresponds to the limit as p ! 1.
Theorem 15.9 (John-Nirenberg, second version). For any p 2 [1, 1) and u 2 Lp (Q)
define
( ✓Z ◆p )
X
p n
Kp (u) := sup L (Qi ) |u(x) uQi | dx | {Qi } partition of Q .
i Qi
The proof of Theorem 15.9 is basically the same as Theorem 15.4, the goal being to
prove the polynomial decay
c(p, n)
|{|u uQ | > t}| Kp (u) t>0
tp
77
instead of an exponential decay.
The following important result improves the classical interpolation theorems in Lp
spaces, replacing L1 with BM O. This is crucial for the application to elliptic PDE’s, as
we will see.
is continuous for every r 2 [p, 1], and its continuity constant c can be bounded indepen-
dently of the chosen partition {Qi }.
78
In order to get information from Theorem 15.9, for r 2 [p, 1), we estimate
X ✓Z ◆r
r n
Kr (T u) = sup L (Qi ) |T u(y) (T u)Qi | dy = sup kT̃{Qi } ukrLr ckukLr ,
{Qi } i Qi {Qi }
where we used the continuity property of T̃ : Lr (Q) ! Lr (Q0 ) stated in (15.4). Therefore,
by Theorem 15.9, we get
kT u (T u)Q kLrw c(r, n, T )kukLr 8u 2 D .
Since u 7! (T u)Q obviously satisfies a similar Lrw estimate, we conclude that kT ukLrw
c(r, n, T )kukLr for all u 2 D. Again, thanks to Marcinkiewicz theorem, with exponents
0 0
p and r, we have that the continuity Lr ! Lr holds for every r0 2 [p, r). Since r is
arbitrary, we got our conclusion. ⇤
We are now ready to employ these harmonic analysis tools to the study of regularity
in Lp spaces for elliptic PDEs, considering first the case of constant coefficients. Suppose
that ⌦ ⇢ Rn is an open, bounded set with Lipschitz boundary @⌦, suppose that the
coefficients A↵ij satisfy the Legendre-Hadamard condition with > 0 and consider the
divergence form of the PDE
⇢
div(Aru) = divF
(15.5)
u 2 H01 (⌦; Rm ) .
In the spirit of Theorem 15.10, we define
T F := ru.
Thanks to Campanato regularity theory, we already got the continuity of T : L2, ! L2,
when 0 < n + 2, thus choosing = n and using the isomorphism (15.3) we see that
T is continuous as an operator
T : L1 (⌦; Rm⇥n ) ! BM O(⌦; Rm⇥n ) . (15.6)
Remark 15.11. Let us remark the importance of weakening the norm in the target
space in (15.6): we passed from L1 (for which, as we will see, no estimate is possible) to
BM O. For BM O the regularity result for PDEs is true and Theorem 15.10 allows us to
interpolate between 2 and 1.
We are going to apply Theorem 15.10 with D = L1 (⌦; Rs ) and s = m ⇥ n. By the
global Caccioppoli-Leray inequality (see Theorem 6.1) we obtain the second hypothesis
of Theorem 15.10: T : L2 (⌦; Rm⇥n ) ! L2 (⌦; Rm⇥n ) is continuous. Therefore
T : D ! Lp (⌦; Rm⇥n ) (15.7)
is (Lp , Lp )-continuous if p 2 [2, 1). Since the (unique) extension of T to the whole of Lp
still maps F into ru, with u solution to (15.5), we have proved the following result:
79
Theorem 15.12. For all p 2 [2, 1) the operator F 7! ru in (15.5) maps Lp (⌦; Rm⇥n )
into Lp (⌦; Rm⇥n ) continuously.
Our intention is now to extend the previous result for p 2 (1, 2), by a duality argument.
G = Br + G̃ , (15.8)
div(G̃) = 0 in ⌦
div(Br ) = divG .
80
the right hand side is Lq only. Notice that the solutions
R obtained in this way have no
variational character anymore, since their energy Aruru dx is infinite (for this reason
they are sometimes called very weak solutions). Since the variational characterization
is lacking, the uniqueness of these solutions needs a new argument, based on Helmholtz
decomposition.
Theorem 15.14. For all q 2 (1, 2) there exists a continuous operator T : Lq (⌦; Rm⇥n ) !
H01,q (⌦; Rm ) mapping F to the unique weak solution u to (15.5).
Proof. We already illustrated the construction of a solution u, by a density argument
and uniform Lq bounds. To show uniqueness, it suffices to show that u 2 H01,q and that
div(Aru) = 0 implies u = 0. To this aim, we define G = |ru|q 2 ru 2 Lp and we apply
Helmholtz decomposition G = A⇤ r + G̃ with 2 H01,p and G̃ 2 Lp divergence-free. By
aR density argumentRw.r.t. u and w.r.t. (notice that the exponents are dual) we have
G̃ru dx = 0 and Arur dx = 0, hence
Z Z Z Z
q ⇤
|ru| dx = Gru dx = A r ru dx = Arur dx = 0.
⌦ ⌦ ⌦ ⌦
81
but its gradient is patently not bounded.
So, it remains to prove that T is necessarily discontinuous, which we will do restricting
our discussion to the scalar case for the sake of simplicity. By the same duality argument
used before, if T were continuous we would get an estimate of the form
and
|Du|(⌦) c|µ|(⌦), (15.11)
where |µ| (resp. |Du|) denote the total variation of the measure µ (resp. Du). On the
other hand, we claim that the inequality (15.11) can’t be true. In fact, when n = 2 and
m = 1, consider the identity matrix A↵ := ↵ and the corresponding Laplace equation
v= 0 , (15.12)
µ1 = 0 + v ⌘ + 2hrv, r⌘i, µ2 = 0,
82
we have that the function w = @x1 (⌘v) 2 L1 (R2 ) is a distributional solution in R2 to the
equation X
w= @ x↵ µ ↵ .
↵
div(Aru) = f + divF
2 2
with coefficients A 2 C(⌦; Rn m ) which satisfy a uniform Legendre-Hadamard condition
for some > 0. Moreover, if p 2 (1, 1), let us suppose that F 2 Lploc and f 2 Lqloc , where
the Sobolev conjugate exponent q ⇤ = qn/(n q) coincides with p. Then |ru| 2 Lploc (⌦).
Proof. We give the proof for p 2 (the other cases come again by duality). Let us fix
s 2 and let us show that
⇤
|ru| 2 Ls^p
loc (⌦) =) |ru| 2 Lsloc^p (⌦) . (15.13)
Proving (15.13) ends the proof because |ru| 2 L2loc (⌦) (case s = 2) and in finitely many
steps s⇤ becomes larger than p.
Fix a point x0 2 ⌦ and a radius R > 0 such that BR (x0 ) b ⌦: we choose a cut-o↵
function ⌘ 2 Cc1 (BR (x0 )), with 0 ⌘ 1 and ⌘ ⌘ 1 in BR/2 (x0 ).
⇤
We claim that ⌘u belongs
⇤
to H01,s ^p (BR (x0 )) if R ⌧ 1, as it is the unique fixed point of a
contraction in H01,s ^p (BR (x0 )), that we are going to define and study in some steps. This
⇤
implies, in particular, that |ru| 2 Ls ^p (BR/2 (x0 )).
(1) We start localizing the equation. Replacing ' with ⌘' in the PDE, by algebraic
83
computations we obtain
Z
A(x)r(⌘u)(x)r'(x) dx
BR (x0 )
Z
= A(x) (⌘(x)ru(x) + u(x) ⌦ r⌘(x)) r'(x) dx
BR (x0 )
Z
= A(x) (ru(x)r(⌘')(x) + u(x) ⌦ r⌘(x)r'(x) ru(x) (r⌘(x)'(x))) dx
BR (x0 )
Z
= f (x)⌘(x)'(x)+F (x)r(⌘')(x)+A(x) (u(x)⌦r⌘(x)r'(x) ru(x)r⌘(x)'(x)) dx
BR (x0 )
Z
= f˜(x)'(x) + F̃ (x)r'(x) dx ,
BR (x0 )
defining
f˜(x) := f (x)⌘(x) + F (x)r⌘(x) A(x)ru(x)r⌘(x)
and
F̃ (x) := F (x)⌘(x) + A(x)u(x) ⌦ r⌘(x) .
Thus ⌘u satisfies
Thanks to the previous Lp regularity result for constant coefficients PDEs, since f˜ 2 Ls^q
loc
(because we assumed that |ru| 2 Ls^p 2 s^q
loc ), we have |r w| 2 Lloc (see also Remark 15.16).
(s^q)⇤
By Sobolev immersion we get |rw| 2 Lloc , hence
⇤ ⇤ ⇤
|rw| 2 Lsloc^q = Lsloc^p .
Now we define ⇤
F ⇤ (x) := F̃ (x) + rw(x) 2 Lsloc^p .
⇤
(3) Let E = H01,s ^p (BR (x0 ); Rm ) and let us define the operator ⇥ : E ! E which
associates to each V 2 E the function v 2 E that solves
84
⇤
The operator ⇥ is well-defined because |F ⇤ | 2 Ls ^p (BR (x0 )) (we saw this in step (2))
and we can take advantage of regularity theory for constant coefficients operators. The
operator ⇥ is a contraction, in fact
1
kr(v1 v2 )kE ck (A(x0 ) A) r(V1 V2 )kE kr(V1 V2 )kLs⇤ ^p (BR (x0 ))
2
if R is sufficiently small, according to the continuity of A. Here we use the fact that the
constant c in the first inequality is scale invariant, so it can be “beaten” by the oscillation
of A in BR (x0 ), if R is small enough.
Let us call v⇤ 2 E the unique fixed point of (15.15). According to (15.14), ⌘u already
1,s^p 1,s⇤ ^p
solves (15.15), but in the larger space H0 . Thus ⌘u 2 H0 if we are able to show
that v⇤ = ⌘u, and to see this it suffices to show that uniqueness holds in the larger space
as well.
Consider the di↵erence v 0 := v⇤ ⌘u 2 H01,s^p (BR (x0 ); Rm ) ⇢ H01 (BR (x0 ); Rm ): v 0 is a
weak solution of
div (A(x)rv 0 ) = 0 ,
hence v 0 ⌘ 0 (we can indeed use the variational characterization of the solution). This
concludes the proof. ⇤
where F 2 C 2, (Rn ) (at least, for some > 0) satisfies the following ellipticity property:
there exist two positive constants ⇤ such that ⇤I r2 F (p) I for all p 2 Rn (this
implies in particular that |r2 F | is uniformly bounded). We have already seen that under
these assumptions it is possible to derive the Euler-Lagrange equations divFp (rv) = 0.
By di↵erentiation, for any direction s 2 {1, . . . , n}, the equation for u := @v/@xs is
✓ ◆
@ @
Fp↵ p (rv) u =0.
@x↵ @x
85
Recall also the fact that, in order to obtain this equation, we needed to work with the
approximation h,s v and with the interpolating operator
Z 1
Aeh (x) := Fpp (trv(x + hes ) + (1 t)rv(x)) dt
0
and to exploit the Caccioppoli-Leray inequality.
One of the striking ideas of De Giorgi was basically to split the problem, that is to
deal with u and v separately, as rv is only involved in the coefficients of the equation for
u. The key point of the regularization procedure is then to show that under no regularity
assumption on rv (i.e. not more than measurability), if u is a solution of this equation,
0,↵
then u 2 Cloc (⌦), with ↵ depending only on n and on the ellipticity constants , ⇤. If
this is true, we can proceed as follows:
u 2 C 0,↵ ) v 2 C 1,↵ ) Fpp (rv) 2 C 0,↵ ) u 2 C 1,↵ ,
where the implications rely upon the fact that Fpp is Hölder continuous and on the
Schauder estimates of Theorem 10.4. Since u is any partial derivative of v, we eventually
get v 2 C 2,↵ . If F is more regular, by continuing this iteration (now using Shauder
regularity for PDE’s whose coefficients are C 1, , C 2, and so on) we obtain
F 2 C1 ) v 2 C1
and also, by the tools developed in [20], that F 2 C ! ) v 2 C ! , which is the complete
solution of the problem raised by Hilbert.
Actually, we have solved this problem in the special case n = 2, since, by means of
Widman’s technique, we could prove that |ru| 2 L2,↵ and hence u 2 L2,↵+2 for some
↵ > 0. This is enough, if n = 2, to conclude that u 2 C 0,↵/2 .
First of all, let us fix our setting. Let ⌦ be an open domain in Rn , 0 < ⇤ < 1
and let A↵ be a Borel symmetric matrix satisfying a.e. the condition I A(x) ⇤I.
1
We want to show that if u 2 Hloc solves the problem
div (A(x)ru(x)) = 0
0,↵
then u 2 Cloc . Some notation is needed: for B⇢ (x) ⇢ ⌦ we define
A(k, ⇢) := {u > k} \ B⇢ (x) ,
where the dependence on the center x can be omitted. This should not create confusion,
since we will often work with a fixed center. In this section, we will derive many func-
tional inequalities, but typically we are not interested in finding the sharpest constants,
but only on the functional dependence of these quantities. Therefore, in order to avoid
complications of the notation we will use the same symbol (generally c) to indicate di↵er-
ent constants, possibly varying from one passage to the next one. However we will try to
indicate the functional dependence explicitly whenever this is appropriate and so we will
use expressions like c(n) or c(n, , ⇤) many times.
86
Theorem 16.1 (Caccioppoli inequality on level sets). For any k 2 R and B⇢ (x) ⇢
BR (x) b ⌦ we have
Z Z
2 c
|ru| dy (u k)2 dy (16.1)
A(k,⇢) (R ⇢)2 A(k,R)
with c = 16⇤2 / 2 .
Remark 16.2. It should be noted that the previous theorem generalizes the Caccioppoli-
Leray inequality, since we do not ask ⇢ = R/2 and we introduce the sublevels.
1,1
Theorem 16.3 (Chain rule). If u 2 Wloc (⌦), then for any k 2 R the function (u k)+
1,1
belongs to Wloc (⌦). Moreover we have that r(u k)+ = ru a.e. on {u > k} , while
r(u k)+ = 0 a.e. on {u k} .
Proof. Since this theorem is rather classical, we just sketch the proof. By the arbitrari-
ness of u, the problem is clearly translation-invariant and we can assume
p without loss of
generality k = 0. Consider the family of functions defined by '" (t) := t + "2 " for t 0
2
and identically zero elsewhere, whose derivatives are uniformly bounded and converge to
1 1,1
{t>0} . Moreover, let (un ) be a sequence of Cloc functions approximating u in Wloc . We
have that for any n 2 N and " > 0 the classical chain rule gives r ['" (un )] = '0" (un )run .
Passing to the limit as n ! 1 gives r ['" (u)] = '0" (u)ru. Now, we can pass to the limit
as " # 0 and use the dominated convergence theorem to conclude that ru+ = {u>0} ru.
⇤
We can come to the proof of the Caccioppoli inequality on level sets.
Proof. Let ⌘ be a cut-o↵ function supported in BR (x), with ⌘ ⌘ 1 on B ⇢ (x) and |r⌘|
2/(R ⇢). If we apply the weak form of our equation to the test function ' := ⌘ 2 (u k)+
we get
Z Z
2
⌘ Aruru dy = 2 ⌘Arur⌘(u k)+ dy
A(k,R) B (x)
Z R Z
⇤ 2 2 4"⇤
⌘ |Du| dy + (u k)2 dy
" A(k,R) (R ⇢)2 A(k,R)
for any " > 0, by our upper bound and by Young inequality. Here we set " = 2⇤/ so
that, thanks to the uniform ellipticity assumption, we obtain
Z Z Z
8⇤2
2
⌘ Aruru dy 2 2
⌘ |ru| dy + 2
(u k)2 dy .
A(k,R) 2 A(k,R) (R ⇢) A(k,R)
Since on the smaller ball ⌘ is identically equal to 1, we eventually get
Z Z
16⇤2
2
|ru| dy 2 2
(u k)2 dy ,
A(k,⇢) (R ⇢) A(k,R)
87
The second great idea of De Giorgi was that (one-sided) regularity could be achieved
for all functions satisfying the previous functional inequality, regardless of the fact that
these were solutions to an elliptic equation. For this reason he introduced a special class
of objects.
Definition 16.4 (De Giorgi’s class). We define the De Giorgi class DG+ (⌦) as follows:
DG+ (⌦) := {u | 9 c 2 R s.t. 8 k 2 R, Br (x) b BR (x) b ⌦, u satisfies (16.1) } .
In this case, we also define c+
DG (u) to be the minimal constant larger than 1 for which
the condition (16.1) is verified.
Remark 16.5. From the previous proof, it should be clear that we do not really require
u to be a solution, but just a sub-solution of our problem. In fact, we have proved that
16⇤2
div (Aru) 0 in D0 (⌦) =) u 2 DG(⌦), c+
DG (u) 2
.
In a similar way, the class DG (⌦) (corresponding to supersolutions) and cDG (u) could
be defined by
Z Z
2 c
|ru| dy 2
(u k)2 dy
{u<k}\B⇢ (x) (R ⇢) {u<k}\BR (x)
The main part of the program by De Giorgi can be divided into two steps:
(i) If u 2 DG+ (⌦), then it satisfies a strong maximum principle in a quantitative form
(more precisely the L2 to L1 estimate in Theorem 16.8);
0,↵
(ii) If both u and u belong to DG+ (⌦), then u 2 Cloc (⌦).
Let us start by discussing the first point. We define these two crucial quantities:
Z
U (h, ⇢) := (u h)2 dy, V (h, ⇢) := L n (A(h, ⇢)) .
A(h,⇢)
88
Proof. The first statement the first inequality in the second statement are trivial, since
Z Z
2 2
(h k) V (h, ⇢) = (h k) dy (u k)2 dy
A(h,⇢) A(h,⇢)
Z
(u k)2 dy = U (k, ⇢) .
A(k,⇢)
For the second inequality, let us introduce a Lipschitz cut-o↵ function ⌘ supported in
B(R+⇢)/2 (x) with ⌘ ⌘ 1 on B ⇢ (x) and |r⌘| 4/(R ⇢). We need to note that
Z Z
4c+
DG (u)
⌘ 2 |r(u k)+ |2 dy (u k)2 dy
B(R+⇢)/2 (R ⇢)2 A(k,R)
and Z Z
+ 2 2 16
((u k) ) |r⌘| dy (u k)2 dy .
B(R+⇢)/2 (R ⇢)2 A(k,R)
and by the Sobolev embedding inequality with the function ⌘(u k)+ this implies
✓Z ◆2/2⇤ Z
2⇤ c(n) · c+
DG (u)
(u k) dy (u k)2 dy
A(k,⇢) (R ⇢)2 A(k,R)
for some constant c(n) depending on the dimension n. In order to conclude, we just need
to apply Hölder’s inequality, in fact
Z ✓Z ◆2/2⇤
2 2⇤
U (k, ⇢) = (u k) dy (u k) dy V (k, ⇢)2/n
A(k,⇢) A(k,⇢)
89
We can view these inequalities as joint decay properties of U and V ; in order to get
the decay of a single quantity, it is convenient to define ' := U ⇠ V ⌘ for some choice of the
(positive) real parameters ⇠, ⌘ to be determined. We obtain:
C⇠ 1
U ⇠ (h, ⇢)V ⌘ (h, ⇢) U ⇠+⌘ (k, R)V 2⇠/n (k, R).
(h k) (R ⇢)2⇠
2⌘
where C := c(n) · c+DG (u), a convention that will be systematically adopted in the sequel.
Since we are looking for some decay inequality for ', we look for solutions (✓, ⇠, ⌘) to the
system
2⇠
⇠ + ⌘ = ✓⇠, = ✓⌘ .
n
Setting ⌘ = 1 (by homogeneity this choice is not restrictive), we get ⇠ = n✓/2 and we can
use the first equation to get r
1 1 2
✓= + + . (16.2)
2 4 n
Note that ✓ > 1 : this fact will play a crucial role in the following proof. In any case, we
get the decay relation
C⇠ 1
'(h, ⇢) '✓ (k, R) .
(h k) (R ⇢)2⇠
2⌘
Theorem 16.7. Let u 2 DG+ (⌦), BR0 (x) b ⌦. For any h0 2 R there exists d =
d(h0 , R0 , c+
DG (u)) such that '(h0 + d, R0 /2) = 0. Moreover, we can take
'(h0 , R0 )✓ 1
d2 = c0 (n)[c+
DG (u)]
n✓/2
,
R0n✓
Corollary 16.8 (L2 to L1 estimate). If u 2 DG+ (⌦), then for any BR0 (x) ⇢ ⌦ and for
any h0 2 R
✓ Z ◆1/2 ✓ ◆(✓ 1)/2
00 1 V (h0 , R0 )
ess sup u h0 + c (n)[c+
DG (u)]
n✓/4
(u 2
h0 ) dy .
BR0 /2 (x) !n R0n A(h0 ,R0 ) R0n
Proof. This corollary comes immediately from Theorem 16.7, once we express ' in terms
of U and V and recall that ⇠ + 1 = ✓⇠ (that is ⇠(✓ 1) = 1), by means of simple algebraic
computations. ⇤
90
Remark 16.9. From Corollary 16.8 with h0 = 0, we can get the maximum principle for
u, as anticipated above. In fact
Z
+ 2 + n✓/2
ess sup (u ) q(n)[cDG (u)] u2 dy
BR0 /2 (x) BR0 (x)
This is true for any µ 2 R but we fix it so that (2⇠ +2) = µ(✓ 1), leading to a cancellation
of two factors in the previous inequality. Having chosen µ, if we choose d as follows
2⇠
2µ C ⇠ 24⇠+2 ✓ 1
0 R0 d
2
=1
then 1 0 . Hence, 2µ C ⇠ 24⇠+2 1✓ 1 R0 2⇠ d 2 1 and the decay inequality yields 2 1 .
By induction, it follows that p 0 , 8p 2 N. In that case, '(kp , Rp ) 2 µp '(h0 , R0 ) !
0 and, since by monotonicity
'(h0 + d, R0 /2) '(kp , R0 /2) '(kp , Rp ) ,
we get the thesis. But the previous condition on d is satisfied if
d2 c0 (n)[c+
DG (u)]
n✓/2
R0 2⇠ ✓ 1
0
When no confusion arises, we will omit the explicit dependence on the center of the ball,
thus identifying !(r) = !(Br (x)).
91
It is an immediate consequence of the previous results that if u 2 DG+ (⌦) \ DG (⌦),
then ✓Z ◆ ✓Z ◆
1/2 1/2
2 2
ess sup u ⇣ u dy , ess inf u ⇣ u dy
Br/2 (x) Br (x) Br/2 (x) Br (x)
for a constant ⇣, which is a function of the dimension n and of cDG (u). Here and in the
sequel we shall denote by cDG (u) the maximum of c+ DG (u) and cDG (u) and by DG(u) the
intersection of the spaces DG+ (⌦) and DG (⌦).
Consequently, under the same assumptions,
✓Z ◆1/2
2
!(Br/2 (x))(u) 2⇣ u dy .
Br (x)
Let us see the relation between the decay of the oscillation of u and the Hölder regu-
larity of u. We prove this result passing through the theory of Campanato spaces (a more
elementary proof is based on the observation that the Lebesgue representative defined at
approximate continuity points is Hölder continuous).
Theorem 16.11. Let ⌦ ⇢ Rn be open, c 0, ↵ 2 (0, 1] and let u : ⌦ ! R be a measurable
0,↵
function such that for any Br (x) ⇢ ⌦ we have !(Br (x)) cr↵ . Then u 2 Cloc (⌦), that
0,↵
is, there exists in the Lebesgue equivalence class of u a Cloc representative.
Proof. By definition of essential extrema, for L n -a.e. y 2 Br (x) we have that
ess inf Br (x) u u(y) ess supBr (x) u. These inequalities imply ess inf Br (x) u uBr (x)
ess supBr (x) and hence that L n -a.e. in Br (x) the inequality |u uBr (x) | cr↵ holds. We
0,↵
have proved that u 2 L2,n+2↵ (⌦), but this gives u 2 Cloc (⌦) (regularity is local since no
assumption is made on ⌦), which is the thesis. ⇤
This theorem motivates our interest in the study of oscillation of u, that will be carried
on by means of some tools we have not introduced so far.
92
Theorem 16.12 (Isoperimetric inequality). Let E ⇢ Rn be a regular set such that
n 1 (@E) < 1. Then
1⇤
min {L n (E), L n (Rn \ E)} c(n) [ n 1 (@E)]
93
Hence, by the Gauss-Green theorem (with ⌫t outer normal to {u > t}) we obtain
Z Z 1 Z ! Z 1
|ru| dx sup h', ⌫t i d n 1 dt n 1 (⌦ \ {u = t}) dt ,
⌦ 0 '2Cc1 , |'|1 ⌦\{u=t} 0
again exploiting the fact that for a.e. t the set {u = t} is the (regular) boundary of
{u > t}.
Let us consider the converse inequality, namely
Z Z 1
|ru| dx n 1 (⌦ \ {u = t}) dt .
⌦ 0
Lemma 16.16. Let G : [0, 1) ! [0, 1) a non-increasing measurable function. Then for
any ↵ 1 we have Z 1 ✓Z 1 ◆↵
↵ 1 1/↵
↵ t G(t) dt G (t) dt .
0 0
Proof. It is sufficient to prove that for any T > 0 we have the finite time inequality
Z T ✓Z T ◆↵
↵ 1 1/↵
↵ t G(t) dt G (t) dt . (16.3)
0 0
which is equivalent to
✓Z t ◆↵ 1
↵ 1 1/↵
t G(t) G (s) ds G1/↵ (t) .
0
94
We are now ready to derive the Sobolev inequalities stated in Theorem 4.6.
Theorem 16.17 (Sobolev embedding, p = 1). For any u 2 W 1,1 (Rn ) we have that
✓Z ◆1/1⇤ Z
1⇤
|u| dx c(n) |ru| dx .
Rn Rn
thanks to Lemma 16.16. Consequently, the isoperimetric inequality and the coarea for-
mula give
Z ✓Z 1 ◆1 ⇤ ✓Z ◆1 ⇤
1⇤
u dx c(n) n 1 {u = t} dt = c(n) |ru| dx .
Rn 0 Rn
The continuous embedding in (2) follows by the global one in (1) applied to an extension
of u (recall that regularity of ⌦ yields the existence of a continuous extension operator
from W 1,1 (⌦) to W 1,1 (Rn )). ⇤
Theorem 16.18 (Sobolev embeddings, 1 < p < n). For any u 2 W 1,p (Rn ) we have that
✓Z ◆1/p⇤ ✓Z ◆1/p
p⇤ p
|u| dx c(n, p) |ru| dx .
Rn Rn
95
and, by Hölder’s inequality, the right hand side can be estimated from above with
Z 1/p0 Z 1/p
(↵ 1)p0
c(n)↵ u dx |ru|p dx .
Rn Rn
but 1/1⇤ 1/p0 = 1/p⇤ , ↵1⇤ = p⇤ and the claim follows. The second part of the statement
can be obtained as in Theorem 16.17. ⇤
We will also make use of the following refinement of the Poincaré inequality in W01,1 :
even though no assumption is made on the behaviour of u at the boundary of the domain,
⇤
it is still possible to control the L1 norm with the gradient.
Theorem 16.19. Let u 2 W 1,1 (BR ) with u 0 and suppose that L n ({u = 0})
n
L (BR )/2. Then
✓Z ◆1/1⇤ Z
1⇤
u dx c(n) |ru| dx .
BR BR
⇤
Proof. This result is the local version of the embedding W 1,1 ,! L1 . Hence, in order
to give the proof, it is just needed to mimic the previous argument substituting the
isoperimetric inequality with the relative isoperimetric inequality, that is, here
⇤
L n (BR \ {u > t}) c(n) n 1 [L n (BR \ {u = t})]1 .
96
Proof. Take two levels h, k such that M h k k0 and define u := u ^ h u ^ k =
(u ^ h k)+ . By construction u 0 and since u 2 W 1,1 (⌦) we also have u 2 W 1,1 (⌦). It
is also clear that ru 6= 0 only on A(k, r) \ A(h, r). Notice that
1 n
L n ({u = 0} \ Br ) L n ({u k} \ Br ) L n ({u k0 } \ Br ) L (Br )
2
and so we can apply the relative version of the critical Sobolev embedding and Hölder’s
inequality to get
Z ✓Z ◆1⇤
1⇤ n 1⇤
(h k) L (A(h, r)) = u dy c(n) |ru| dy
A(h,r) Br
Z
= c(n) |ru| dy
A(k,r)\A(h,r)
✓Z ◆1⇤ /2
2 ⇤ /2
c(n) |ru| dy L n (A(k, r) \ A(h, r))1 .
A(k,r)
in order to obtain
⇤
(h k)2 L n (A(h, r))2/1 c(n)c+
DG (u)(M k)2 rn 2 (V (k, r) V (h, r)) . (16.5)
Here we can conclude the proof by applying (16.5) for h = ki+1 and k = ki , so that
⌫
X
2/1⇤ ⇤
⌫V (k⌫ , r) V (ki , r)2/1
i=1
⌫
X
4c(n)c+
DG (u)r
n 2
[V (ki , r) V (ki+1 , r)]
i=1
4c(n)c+
DG (u)!n r
2n 2
.
Theorem 16.21 (C 0,↵ regularity). Let ⌦ ⇢ Rn be open and let u 2 DG(⌦). Then
0,↵
u 2 Cloc (⌦), with 2↵ = log2 1 2 (⌫+2) ,
97
Proof. Pick an R > 0 such that B2R (x) b ⌦ and consider for any r R the functions
m(r) := ess inf Br (x) (u) and M (r) := ess supBr (x) (u). Moreover, set !(r) = M (r) m(r)
and µ(r) := (m(r) + M (r)) /2. We apply the previous lemma to the sequence k⌫ :=
M (2r) !(2r)
2⌫+1
, but to do this we should check the hypothesis (16.4), which means
1
L n ({u > µ(2r)} \ Br (x)) L n (Br (x)).
2
Anyway, either L n ({u > µ(2r)} \ Br (x)) 12 L n (Br (x)) or L n ({u < µ(2r)} \ Br (x))
1
2
L n (Br (x)). The second case is analogous, provided we work with u instead of u, and
it is precisely here that we need the assumption that both u and u belong to DG+ (⌦).
Using Lemma 16.20 it is easily seen that the choice of ⌫ as in (16.6), with c(n) large
enough, provides
✓ ◆
00
⇥ + ⇤n✓/4 V (k⌫ , r) (✓ 1)/2 1
c (n) cDG (u) ,
rn 2
where c00 (n) is the dimensional constant in Theorem 16.8. Moreover, this choice of ⌫ has
been made independently of of r and R (this is crucial for the validity of the scheme
below).
Now apply the maximum principle in Theorem 16.8 to u with radii r/2 and r and
h0 = M (2r) !(2r)2⌫+1
= k⌫ (for the previous choice of ⌫) to obtain
⇣r⌘ ✓ ◆(✓ 1)/2
00
⇥ ⇤n✓/4 V (h0 , r)
M h0 + c (n) c+
DG (u) (M (2r) h0 )
2 rn
which is the desired decay estimate. By the standard iteration argument5 , we find
✓ ◆↵
↵ r
!(r) 4 !(R) 0<rR
R
5
We refer to Lemma 9.1, with the obvious changes.
98
17 Regularity for systems
17.1 De Giorgi’s counterexample to regularity for systems
In the previous section we saw De Giorgi’s regularity result for solutions u 2 H 1 (⌦) of
the elliptic equation
div (A(x)ru(x)) = 0
0,↵
with bounded Borel coefficients A satisfying I A ⇤I. It turned out that u 2 Cloc (⌦),
with ↵ = ↵(n, , ⇤).
It is natural to investigate about similar regularity properties for systems, still under
no regularity assumption on A (otherwise, Schauder theory is applicable). In 1968, in [8],
Ennio De Giorgi provided a counterexample showing that the scalar case is special. De
Giorgi’s example not only solves an elliptic PDE, but it is also the minimum of a convex
variational problem.
When m = n, consider
u(x) := x|x|↵ . (17.1)
We will show in (17.7), (17.8) and (17.9) that, choosing
!
n 1
↵= 1 p , (17.2)
2 (2n 2)2 + 1
If n / L1 (B1 ), because
3 then |u| 2
! ✓ ◆
n 1 3 1
↵= 1 p 1 p >1
2 (2n 2)2 + 1 2 17
and this provides a counterexample not only to Hölder regularity, but also to local boun-
dedness of solutions. In the case n = 2 we already know from Widman’s technique (see
Remark 4.4) that u is locally Hölder continuous, nevertheless De Giorgi’s example will
show that this regularity cannot be improved Rto local Lipschitz.
Calling A(x) the matrix such that L(u) = B1 hA(x)ru, rui dx, we remark that A has
a discontinuity at the origin (determined by the term x ⌦ x/|x|2 ).
99
The Euler-Lagrange equation associated to (17.3) is the following (in the weak distri-
butional sense): for every h = 1, . . . , n it must be
n n
!
@ X @ut X xs xt @ut
0 = (n 2) (n 2) +n (17.4)
@xh t=1
@xt s,t=1
|x|2 @xs
n
" n n
!#
X @ xh xk X @ut X xs xt @ut
+ n (n 2) +n (17.5)
k=1
@xk |x|2 t=1
@xt s,t=1
|x|2 @xs
n
X @ 2 uh
+ . (17.6)
k=1
@x2k
We are going to prove in a few steps that u is the unique minimizer of L, with boundary
data given by u itself, and that u solves the Euler-Lagrange equations. The steps are:
(i) u, as defined in (17.1), belongs to C 1 (B1 \ {0}; Rn ) and solves in B1 \ {0} the
Euler-Lagrange equations;
(ii) u 2 H 1 (B1 ; Rn ) and is also a weak solution in B1 of the system.
Let us perform step (i). Fix h 2 {1, . . . , n}, and use extensively the identity
@
|x|↵ = ↵xh |x|↵ 2
.
@xh
Then |x|↵ = (n↵ + ↵2 2↵)|x|↵ 2
and
@
(xh |x|↵ ) = xh |x|↵ + |x|↵ = (↵n + ↵2 )xh |x|↵ 2
(17.7)
@xh
and this is what we need to put in (17.6) when u is given by (17.1). For both (17.4) and
(17.5) we have to calculate
Xn
@
(xt |x|↵ ) = (n + ↵)|x|↵ ,
t=1
@x t
and n n
X xs xt @ut X xs xt
2
= 2
↵xs xt |x|↵ 2
+ st |x|
↵
= (↵ + 1)|x|↵ .
s,t=1
|x| @x s s,t=1
|x|
Therefore (17.4) is given by
n n
!
@ X @ut X xs xt @ut
(n 2) (n 2) +n = ↵(n 2)[(n 2)(n+↵)+n(↵+1)]xh |x|↵ 2
.
@xh t=1
@xt s,t=1
|x|2 @xs
(17.8)
100
In order to compute the term (17.5) we first get
Xn
@
xh xk |x|↵ 2
= (n + ↵ 1)xh |x|↵ 2
,
k=1
@x k
Putting together (17.7), (17.8) and (17.9), u(x) = x|x|↵ solves the Euler-Lagrange equa-
tion if and only if ⇣ n ⌘2
(2n 2)2 ↵ + + ↵n + ↵2 = 0 ,
2
which leads to the choice (17.2) of ↵.
Let us now perform step (ii), checking first that u 2 H 1 . As |ru(x)| ⇠ |x|↵ and
2↵ > n, it is easy to show that |ru| 2 L2 (B1 ). Moreover, for every ' 2 Cc1 (B1 \ {0})
we have classical integration by parts formula
Z Z
ru(x)'(x) dx = u(x)r'(x) dx . (17.10)
Thanks to Lemma 17.1 below, we are allowed to approximate in H 1 (B1 ; Rn ) norm every
' 2 Cc1 (B1 ) with a sequence ('k ) ⇢ Cc1 (B1 \ {0}). Then we can pass to the limit in
(17.10) because |ru| 2 L2 (B1 ) to obtain u 2 H 1 (B1 ; Rm ). Now, using the fact that the
Euler-Lagrange PDE holds in the weak sense in B1 \ {0} (because it holds in the classical
sense), we have Z
A(x)ru(x)r'(x) dx = 0 (17.11)
B1
for every ' 2 Cc1 (B1 \ {0}; Rn ). Using Lemma 17.1 again, we can extend (17.11) to every
' 2 Cc1 (B1 ; Rn ), thus obtaining the validity of the Euler-Lagrange PDE in the weak sense
in the whole ball.
Finally, since the functional L in (17.3) is convex, the Euler-Lagrange equation is
satisfied by u if and only if u is a minimizer of L(u) with boundary condition
u(x) = x in @B1 .
This means that De Giorgi’s counterexample holds not only for solution of the Euler-
Lagrange equation, but also for minimizers.
101
Lemma 17.1. Assume that n > 2. For every ' 2 Cc1 (B1 ) there exists 'k 2 Cc1 (B1 \{0})
such that 'k tends to ' strongly in W 1,2 (B1 ).
Proof. Consider 2 Cc1 (Rn ) with ⌘ 1 on B 1 , then rescale setting k (x) := (kx).
2
Set 'k := '(1 k ); in L topology we have ' 'k = ' k ! 0 and (r') k ! 0. Since
r(' 'k ) = (r') k + 'r k ,
the thesis is equivalent to verify that
Z
'(x)2 |r k (x)|
2
dx ! 0 ,
B1
but
Z Z
2 2 2 2
'(x) |r k (x)| dx (sup ' )k |r (kx)|2 dx
B1 B1
Z
2 2 n
(sup ' )k |r (x)|2 dx ! 0 ,
Rn
102
18 Partial regularity for systems
As we have seen with De Giorgi’s counterexample, it is impossible to expect an “every-
where” regularity result for elliptic systems: the main idea is to pursue a di↵erent goal, a
“partial” regularity result, away from a singular set. This strategy goes back to De Giorgi
himself, and it was implemented for the first time in the study of minimal surfaces.
Definition 18.1 (Regular and singular sets). For a generic function u : ⌦ ! R we call
regular set of u the set
The set ⌦reg (u) is obviously the largest open subset A of ⌦ such that u 2 C 1 (A).
Briefly, let us recall here the main results we have already obtained for elliptic systems.
then
@
Fp↵i (ru) = 0 8 i = 1, . . . , m .
@x↵
(b) If F satisfies a uniform Legendre condition for some > 0, by Nirenberg method we
have ru 2 Hloc 1
(⌦; Rm⇥n ) and (by di↵erentiation of the (EL) equations with respect
to x )
✓ ◆
@ @ 2 uj
Fp↵ p (ru) = 0 8 i = 1, . . . , m, = 1, . . . , n . (18.1)
@x↵ i j @x @x
In this section we shall provide a fairly complete proof of the following result, following
with minor variants the original proof in [10].
103
Theorem 18.3 (Evans). If F 2 C 2 (Rm⇥n ) is -uniformly quasiconvex with > 0 and
satisfies
|r2 F (p)| ⇤ 8p 2 Rm⇥n , (18.2)
for some ⇤ > 0. Then any local minimizer u belongs to C 1, (⌦reg ) for some =
(n, m, , ⇤) and
L n (⌦ \ ⌦reg ) = 0 .
The following list summarizes some results in the spirit of Theorem 18.3. At this stage
we should point out that the growth condition (18.2) is a bit restrictive if we want to allow
the standard examples of quasiconvex functions, i.e. convex functions of determinants of
minors of ru; it includes for instance functions of the form
s X
2
F (ru) := |ru| + 1 + (M ru)2
M
(ii) If r2 F I for some > 0 and it is globally uniformly continuous, then we have
even H n 2
(⌃(u)) = 0.
(iii) If u is locally Lipschitz, then Kristensen and Mingione proved in [23] that there
exists > 0 such that
H n (⌃(u)) = 0 .
(iv) On the contrary, when n = 2 and m = 3, there exists a Lipschitz solution u for
the system @x@↵ Fp↵i (ru) (with F smooth and satisfying the Legendre-Hadamard
condition), provided in [25], such that
⌦reg (u) = ; .
This last result clarifies once for all that partial regularity can be expected for (local)
minimizers only. We will see how local minimality (and not only the validity of the
Euler-Lagrange equations) plays a role in the proof of Evans’ result.
104
We will start with a decay lemma relative to constant coefficients operators.
Lemma 18.4. There exists a constant C⇤ = C⇤ (n, m, , ⇤) 2 (0, 1) such that, for every
constant matrix A satisfying the Legendre-Hadamard condition with and the inequality
|A| ⇤, any solution u 2 H 1 (Br ; Rm ) of
div(Aru) = 0 in Br
satisfies
Z Z
2 2
|ru(x) (ru)B↵r | dx C⇤ ↵ |ru(x) (ru)Br |2 dx 8↵ 2 (0, 1) .
B↵r Br
Proof. As a consequence of what we proved in the section about decay estimates for
systems with constant coefficients, considering (5.2) with ⇢ = ↵r and ↵ < 1, we have that
Z ⇣ ↵r ⌘n+2 Z
2
|ru(x) (ru)B↵r | dx c(n, m, , ⇤) |ru(x) (ru)Br |2 dx . (18.4)
B↵r r Br
1
Definition 18.5 (Excess). For any function u 2 Hloc (⌦; Rm ) and any ball B⇢ (x) b ⌦ the
excess of u in B⇢ (x) is defined by
Z !1/2
Exc (u, B⇢ (x)) := |ru(y) (ru)B⇢ (x) |2 dy .
B⇢ (x)
When we consider functions F satisfying the more general growth condition (18.3),
then we should modify the definition of excess as follows, see [10]:
Z
2
Exc (u, B⇢ (x)) = 1 + |ru(y) (ru)B⇢ (x) |q 2 |ru(y) (ru)B⇢ (x) |2 dy.
B⇢ (x)
Remark 18.6 (Properties of the excess). We list here the basic properties of the excess,
they are trivial to check.
105
(i) Any additive perturbation by an affine function p(x) does not change the excess,
that is
Exc (u + p, B⇢ (x)) = Exc (u, B⇢ (x)) .
Remark 18.7. The name “excess” is inspired by De Giorgi’s theory of minimal surfaces,
presented in [6] and [7], see also [15] for a modern presentation. The excess of a set E at
a point is defined (for regular sets) by
Z
Exc (E, B⇢ (x)) := |⌫E (y) ⌫E (x)|2 d n 1 (y) ,
B⇢ (x)\@E
where ⌫E is the inner normal of the set E. The correspondence between Exc (u, B⇢ (x))
and Exc (E, B⇢ (x)) can be made more evident seeing near x the set @E as the graph
associated to a function
p u, in a coordinate system where ru(x) = 0. Indeed, the identity
⌫E = ( ru, 1)/ 1 + |ru|2 and the area formula for graphs give
Z Z p Z
2
|⌫E (y) ⌫E (x)| d n 1 (y) = 2 2
1 + |ru(z)| 1 dz ⇠ |ru(z)|2 dz ,
B⇢ (x)\@E ⇡(B⇢ (x)\@E) B⇢ (z)
where ⇡(B⇢ (x) \ @E) denotes the projection of the B⇢ (x) \ @E on the hyperplane.
The main ingredient in the proof of Evans’ theorem will be the decay property of the
excess: there exists a critical threshold such that, if the decay in the ball is below the
threshold, then decay occurs in the smaller balls.
Theorem 18.8 (Excess decay). Let F be as in Theorem 18.3. For every M 0 and all
↵ 2 (0, 1/4) there exists "0 = "0 (n, m, , ⇤, M, ↵) > 0 satisfying the following implication:
if
R
(a) u 2 H 1 (Br (x); Rm ) is a local minimizer in Br (x) of v 7! F (rv) dx,
106
then
Exc (u, B↵r (x)) Ce ↵ Exc (u, Br (x))
with Ce depending only on (n, m, , ⇤). When r2 F is uniformly continuous, condition
(b) is not needed for the validity of the implication and "0 is independent of M .
Proof. We choose Ce in such a way that 16C⇤2 CP C ⇤ < Ce2 , where C⇤ is the constant
of Lemma 18.4, CP is the constant in the Poincaré inequality and C ⇤ is the constant of
Proposition 18.9 below.
The proof is by contradiction, assuming that the statement fails for some ↵ and M
(for simplicity we keep F fixed in the contradiction argument, but a slightly more complex
proof would give the stronger result): in step (ii) we will normalize the excesses, obtaining
functions wk with Exc (wk , B↵ (0)) Ce ↵ while Exc (wk , B1 (0)) = 1. Each wk is a solution
of
@
Fp↵i (rwk ) = 0 .
@x↵
We will see in step (iii) that, passing through the limit as k ! 1, any limit point w1
w.r.t. the weak H 1 topology solves
⇣ ⌘
div Fp↵ p (p1 )rw1 = 0.
i j
Using Lemma 18.4 in combination with Proposition 18.9 we will reach the contradiction.
(i) By contradiction, we have M 0, ↵ 2 (0, 1/4) and local minimizers uk : ⌦ ! Rm in
Brk (xk ) with
"k := Exc (uk , Brk (xk )) ! 0
satisfying
(ruk )Brk (xk ) M (18.5)
but
Exc (uk , B↵rk (xk )) > Ce ↵ Exc (uk , Brk (xk )) 8k 2 N .
(ii) Suitably rescaling and translating the functions uk , we can assume that xk = 0,
rk = 1 and (uk )B1 = 0 for all k. Setting pk := (ruk )B1 , the hypothesis (18.5) gives, up to
subsequences,
pk ! p1 2 Rm⇥n . (18.6)
We start here a parallel and simpler path through this proof, in the case when r2 F is
uniformly continuous: in this case no uniform bound on pk is needed and we can replace
(18.6) with
2 2
r2 F (pk ) ! A1 2 Rm ⇥n . (18.7)
Notice that (18.7) holds under (18.6), simply with A1 = r2 F (p1 ). Notice also that,
in any case, A1 satisfies a (LH) condition with constant (this can be achieved using
107
oscillating test functions, as we did to show that quasi-convexity implies the Legendre-
Hadamard condition) and |A1 | ⇤
We do a second translation in order to annihilate the mean of the gradients of mini-
mizers: let us define
vk (x) := uk (x) pk (x) ,
so that (vk )B1 = 0 and (rvk )B1 = 0. According to property (i) of Remark 18.6 the excess
does not change, so still
Exc (vk , B1 ) = "k ! 0
and
Exc (vk , B↵ ) > Ce ↵ "k .
During these operations, we need not lose sight of the variational problem we are solving,
for example every function vk minimizes the integral functional associated to
p 7! F (p + pk ) F (pk ) rF (pk )p .
In order to get some contradiction, our aim is to find a “limit problem” with some decaying
property. Let us define
vk
wk := k 2 N.
"k
It is trivial to check that (wk )B1 = (rwk )B1 = 0, moreover
1
Fk (p) := [F ("k p + pk ) F (pk ) rF (pk )"k p] .
"2k
Here we used the fact local minimality w.r.t. to an integrand F is preserved if we multiply
F by a positive constant or add to F an affine function.
(iii) We now study both the limit of Fk and the limit of wk , as k ! 1. Since Fk 2
C 2 (Rm⇥n ), by Taylor expansion we are able to identify a limit Lagrangian, given by
1
F1 (p) = hA1 p, pi ,
2
to which Fk (p) converge uniformly on compact subsets of Rm⇥n . Indeed, this is clear
with A1 = r2 F (p1 ) in the case when pk ! p1 ; it is still true with A1 given by
(18.7) when r2 F is uniformly continuous, writing Fk (p) = 12 hr2 F (pk + ✓"k p)p, pi with
✓ = ✓(k, p) 2 (0, 1).
108
Once we have the limit problem defined by F1 , we drive our attention to wk : it is
a bounded sequence in H 1,2 (B1 ; Rm ) because the excesses are constant, so by Rellich
theorem we have that (possibly extracting one more subsequence)
wk ! w1 in L2 (B1 ; Rm )
and, as a consequence,
rwk * rw1 in L2 (B1 ; Rm ) . (18.9)
The analysis of the limit problem now requires the verification that w1 solves the Euler
equation associated to F1 . We need just to pass to the limit in the (EL) equation satisfied
by wk , namely
X Z 1 ✓ @F @F
◆ i
@
↵
(pk + "k rwk (x)) ↵
(pk ) (x) dx = 0 8 ' 2 Cc1 (B1 ; Rm ) .
" @pi
↵,i B1 k
@pi @x↵
Writing the di↵erence quotient of rF with the mean value theorem and using r2 F (pk ) !
A1 we obtain
Z
hA1 rw1 (x), r'(x)i dx = 0 8 ' 2 Cc1 (B1 ; Rm ) , (18.10)
B1
This can be obtained splitting the integral into the regions {|rwk | L} and {|rwk | > L},
with L fixed. The first contribution goes to zero, thanks to the convergence of pk to p1
or, when pk is possibly unbounded, thanks to the uniform continuity of r2 F . The second
contribution tends to 0 as L " 1 uniformly in k, since |r2 F | ⇤ and krwk k2 1.
(iv) Equality (18.10) means that
in a weak sense: since the equation has constant coefficients we can apply Lemma 18.4 to
get
Z Z
2 2 2
|rw1 (x) (rw1 )B2↵ | dx 4C⇤ ↵ |rw1 (x)|2 dx 4C⇤2 ↵2 . (18.11)
B2↵ B1
109
hence passing to the limit as k ! 1 gives
Z
Ce2 4
↵ |w1 (w1 )2↵ (rw1 )2↵ (x)|2 dx .
C⇤ B2↵
⇤
The following proposition can be considered as a nonlinear Caccioppoli inequality. It
can be derived without using the Euler-Lagrange equation (which would not help) and
using the minimality instead.
Proposition 18.9 (Caccioppoli inequality for minimizers). There exists C ⇤ = C(n, m, , ⇤)
such that if F is -quasiconvex with |r2 F | ⇤ and if u is a local minimizer in ⌦, then
Z Z
2 C⇤
|ru A| dx 2 |u a A(x x0 )|2 dx
Br/2 (x0 ) r Br (x0 )
110
Combining (18.12) with (18.13) we get
Z Z
2
|r | dx |rF (A) rF (ru)||r | + C|r |2 dx ,
Bs Bs
with C = C( , ⇤).
Now, since |r | |ru A| + 2|u Ax|/(s t), we get
Z Z Z
2 2 D
|ru A| dx D |ru A| dx + |u Ax|2 dx
Bt Bs \Bt (s t)2 Br
for some new constant D = D( , ⇤). Now we apply the hole-filling technique to get
Z Z Z
2 2 D
|ru A| dx ✓ |ru A| dx + |u Ax|2 dx .
Bt Bs (s t)2 Br
with ✓ = D/(D +1) < 1. At this point, since the inequality is true for all r/2 t s r,
a standard iteration scheme gives the result. Indeed, let ⌧ 2 (0, 1) with ✓ < ⌧ 2 and define
ti = (1 ⌧ i /2)r, so that t0 = r/2, ti " r and ti+1 ti = r(1 ⌧ )⌧ i /2. By iteration of the
inequality Z Z
2 4D
|ru A| dx ✓ |ru A|2 dx + 2 ⌧ 2i
Bt i Bti+1 r (1 ⌧ )2
we get
Z Z N
X1
2 N 2 4D
|ru A| dx ✓ |ru A| dx + 2 (✓/⌧ 2 )i
Bt 0 BtN r (1 ⌧ )2 i=0
Z
N 4D⌧ 2
✓ |ru A|2 dx +
Br r2 (1 ⌧ )2 (⌧ 2 ✓)
for any integer N 1. As N ! 1 we get the result.
111
• L n ⌃(u) = 0;
• H n 2+" ⌃(u) = 0 for all " > 0 in the uniformly convex case and H n 2
(⌃(u)) = 0
if r2 F is also uniformly continuous.
In order to exploit Theorem 18.8 and prove that L n ⌃(u) = 0, we fix once for all
the constant ↵ 2 (0, 1/4) in such a way that Ce ↵ < 1/2 (recall that Ce depends only on
the dimensions and on the ellipticity constants). Then, we fix M 0, so that there is an
associated "0 = "0 (n, m, , ⇤, M ) for which the decay property of the excess applies with
halving of the excess from the scale r to the scale ↵r.
Definition 18.10. We will call
⌦M (u) := x 2 ⌦ 9Br (x) b ⌦ with (ru)Br (x) < M1 and Exc (u, Br (x)) < "1
where
M1 := M/2 (18.14)
and "1 verifies
2n/2 "1 "0 (18.15)
and for ↵ 2 (0, 1/4) fixed, chosen in such a way that Ce ↵ < 1/2,
(2n+1 + ↵ n 1+n/2
2 )"1 M . (18.16)
Remark 18.11. The set ⌦M (u) ⇢ ⌦ of Definition 18.10 is open, since the inequalities
are strict. Moreover, by Lebesgue approximate
R continuity theorem (that is, if f 2 Lp (⌦),
then for L -almost every x one has Br (x) |f (y) f (x)|p dy ! 0 as r # 0), it is easy to
n
see that
L n ({|ru| < M1 } \ ⌦M (u)) = 0 . (18.17)
Finally, using (18.17), we realize that
! !
[ [
Ln ⌦\ ⌦M (u) = L n ⌦\ {|ru| < M1 } =0. (18.18)
M 2N M 2N
we obtain L n (⌃(u)) = 0. So, the rest of this section will be devoted to the proof of the
inclusion above, with M fixed.
Fix x 2 ⌦M (u), according to Definition 18.10 there exists r > 0 such that Br (x) b ⌦,
|(ru)Br (x) | < M1 and Exc (u, Br (x)) < "1 . We will prove that
112
so let us fix y 2 Br/2 (x).
(1) Thanks to our choice of "1 (see property (18.15) of Definition 18.10) we have
Z !1/2
Exc u, Br/2 (y) = |ru(z) (ru)Br/2 (y) |2 dz
Br/2 (y)
Z !1/2
|ru(z) (ru)Br (x) |2 dz
Br/2 (y)
✓Z ◆1/2
n/2 2
2 |ru(z) (ru)Br (x) | dz = 2n/2 Exc (u, Br (x)) < "0
Br (x)
so, momentarily ignoring the hypothesis that |(ru)Br/2 (y) | should be bounded by M (we
are postponing this to point (2) of this proof), Theorem 18.8 gives tout court
1 1
Exc u, B↵r/2 (y) Exc u, Br/2 (y) < "0 ,
2 2
thus, just iterating Theorem 18.8, we get
As we have often seen through these notes, we can apply an interpolation argument to a
sequence of radii with ratio ↵ to obtain
✓ ◆µ ✓ ◆µ
µ ⇢ µ ⇢
Exc (u, B⇢ (y)) ↵ Exc u, Br/2 (y) ↵ "0 8⇢ 2 (0, r/2], y 2 Br/2 (x)
r/2 r/2
with µ = (log2 (1/↵)) 1 . We conclude that the components of ru belong to the Cam-
panato space L2,n+2µ (Br/2 (x)) and then u belongs to C 1,µ (Br/2 (x)).
(2) Now that we have explained how the proof runs through the iterative application of
Theorem 18.8, we deal with the initially neglected hypothesis, that is |(ru)Br/2 (y) | < M
and, at each subsequent step, |(ru)B↵k r/2 (y) | < M . Remember that in part (1) of this
proof we never used (18.14) and (18.16).
Since x 2 ⌦M (u) and r fulfills Definition 18.10, for the first step it is sufficient to use the
113
triangular inequality in (18.21) and Hölder’s inequality in (18.22): in fact we can estimate
Z
(ru)Br/2 (y) = ru(z) (ru)Br (x) dz + (ru)Br (x)
Br/2 (y)
Z
ru(z) (ru)Br (x) dz + (ru)Br (x) (18.21)
Br/2 (y)
✓ Z ◆
2n
ru(z) (ru)Br (x) dz + (ru)Br (x)
!n rn Br (x)
✓Z ◆1/2
n 2
2 ru(z) (ru)Br (x) dz + (ru)Br (x) (18.22)
Br (x)
for every k 1.
The first step (k = 1) follows from (18.23), because, estimating as in (18.21) and (18.22),
we immediately get
Z
(ru)B↵r/2 (y) ru(z) (ru)Br/2 (y) dz + (ru)Br/2 (y)
B↵r/2 (y)
n
↵ Exc u, Br/2 (y) + (ru)Br/2 (y)
n n/2
↵ 2 "1 + 2 n "1 + M 1 .
Being the first step already proved, we fix our attention on the (k + 1)th step. With the
same procedure, we estimate again
Z
(ru)B↵k+1 r/2 (y) ru(z) (ru)B↵k r/2 (y) dz + (ru)B↵k r/2 (y)
B↵k+1 r/2 (y)
n
↵ Exc u, B↵k r/2 (y) + (ru)B↵k r/2 (y)
k 1
X
n n/2 k n n n/2 j
↵ 2 "1 + M 1 + 2 " 1 + ↵ "1 2 2 (18.25)
j=0
114
where (18.25) is obtained joining the estimate on the excess (18.20) with the inductive
hypothesis (18.24).
In order to carry out our second goal, namely to prove that
Hn 2+"
⌃(u) = 0 8" > 0 ,
while H k is defined by
H k (B) := lim H k (B) , (18.26)
!0
the limit in (18.26) being well defined because 7! H k (B) is non-increasing. The constant
ck 2 (0, 1) will be conveniently fixed in Remark 18.14.
It is easy to check that H k is the counting measure when k = 0 (provided c0 = 1)
and H k is identically 0 when k > n.
The spherical Hausdor↵ measure S k has a definition analogous to Definition 18.12,
but only covers made with balls are allowed, so that
H k S k 2k H k , H k S k 2k H k . (18.27)
H k (B + h) = H k (B) 8 B ⇢ Rn , 8 h 2 Rn ,
H k ( B) = k
H k (B) 8 B ⇢ Rn , 8 > 0 .
(ii) The Hausdor↵ measures are countably subadditive, which means that whenever we
have a countable cover of a subset B, namely B ⇢ [i2I Bi , then
X
H k (B) H k (Bi ) .
i2I
115
(iii) For every set A ⇢ Rn the map B 7! H k (A \ B) is -additive on Borel sets, which
means that whenever we have a countable pairwise disjoint cover of a Borel set B
by Borel sets Bi , we have
X
H k (A \ B) = H k (A \ Bi ) .
i2I
116
Proposition 18.15. Consider a locally finite measure µ 0 on the family of Borel sets
B(Rn ) and, fixing t > 0, set
⇢
µ(B r (x))
B := x | lim sup >t , (18.31)
r!0 !k rk
Theorem 18.16 (Besicovitch). There exists an integer ⇠ = ⇠(n) with the following prop-
erty: if A ⇢ Rn is bounded and ⇢ : A ! (0, 1), there exist sets A1 , . . . , A⇠(n) ⇢ A such
that
(a) for all j = 1, . . . , ⇠, the balls in {B⇢(x) (x)}x2Aj are pairwise disjoint;
Definition 18.17 (Fine cover). A family F of closed balls in a metric space (X, d) is a
fine cover of a set A ⇢ X if
Theorem 18.18. Fix k 0, consider a fine cover F of A ⇢ X, with (X, d) metric space.
Then there exists a countable and pairwise disjoint subfamily F 0 = {B i }i 1 ⇢ F such that
at least one of the following conditions holds:
P
1
(i) [r(Bi )]k = 1,
i=1
✓ ◆
k
S
1
(ii) H A\ Bi = 0.
i=1
117
Proof. The subfamily F 0 is chosen inductively, beginning with F0 := F. Surely, there
exists a closed ball, let us call it B 1 , such that
1
r B1 > sup r(B)| B 2 F0 .
2
Now put
F1 := {B 2 F0 | B \ B 1 = ;} ,
and choose among them a ball B 2 2 F1 such that
1
r(B 2 ) > sup r(B)| B 2 F1 .
2
If we try to go on analogously, the only chance by which the construction has to stop
is that for some l 2 N the family Fl = ;, so we are getting (because the cover is fine)
that the union of the chosen balls covers the whole of A and therefore option (ii) in the
statement.
Otherwise, assuming that the construction does not stop, we get a family F 0 = {B i }i 1 =
{B ri (yi )}i 1 . We prove that if (i) does not hold, and in particular diam(B i ) ! 0, then
we have to find (ii) again. S
Fix an index i0 2 N: for every x 2 A \ i10 B i there exists a ball B r(x) (x) 2 F such
that
i0
[
B r(x) (x) \ Bi = ; ,
i=1
because F is a fine cover of A and the complement of [i10 B i is open in X. On the other
hand, we claim that there exists an integer i(x) > i0 such that
In fact if
B r(x) (x) \ B i = ; 8 i > i0 , (18.33)
then
r(x)
ri 88 i > i0 (18.34)
2
but ri ! 0, so (18.34) leads to a contradiction. Without loss of generality, we can
think that i(x) is the first index larger than i0 for which (18.32) holds, too. Since, by
construction, ri(x) > 12 sup{r(B)|B 2 Fi(x) 1 } (and B r(x) (x) 2 Fi(x) 1 by the minimality
of i(x)), then r(x) 2ri(x) .
Since the balls intersect, the inequality d(x, yi(x) ) r(x) + ri(x) 3ri(x) gives
118
and therefore
i0
[ 1
[
A\ Bi ⇢ B 5ri (yi ) . (18.35)
i=1 i=i0 +1
Choosing i0 such that 10ri < for every i > i0 , (18.35) says that
1
! i0
! 1
[ [ X
Hk A\ Bi H k
A\ Bi !k (10ri )k .
i=1 i=1 i=i0 +1
⇤
Now we are able to prove Proposition 18.15.
Proof. Intersecting B with balls, one easily reduces to the case of a bounded set B.
Hence, we can assume B bounded and µ finite measure. Fix > 0, an open set A B
and consider the family
Since (ii) holds and we can compare H k with S k via (18.27), to get
1
! 1 1
[ X 1X µ(A)
S k (B) S k Bi !k rik < µ(B i ) , (18.37)
i=1 i=1
t i=1
t
As # 0 we get tS k (B) µ(A) and the outer regularity of µ gives tS k (B) µ(B).
Finally, the last statement of the proposition can be achieved noticing that the in-
equality (18.37) gives that S k (B) is finite; if we assume that µ vanishes on sets with
finite k-dimensional measure we obtain that µ(B) = 0; applying once more the inequality
we get S k (B) = 0. ⇤
119
18.3 Partial regularity for systems: H n 2+"
(⌃(u)) = 0
Aware of the usefulness of Proposition 18.15 for our purposes, we are now ready to obtain
that if F 2 C 2 (Rm⇥n ) satisfies the Legendre condition for some > 0 and satisfies also
then we have a stronger upper bound on the size of the singular set, namely
Hn 2+"
(⌃(u)) = 0 8" > 0 , (18.38)
and ⇢
⌃2 (u) := x 2 ⌦ lim sup (ru)Br (x) = +1 .
r!0
• there exists M1 < 1 such that (ru)Br (x) < M1 for arbitrarily small radii r > 0;
120
• thanks to Poincaré inequality
Z
2 2 n
Exc (u, Br (x)) C(n)r |r2 u(y)|2 dy ! 0 ;
Br (x)
thus for some M = M (M1 , n, m, , ⇤) > 0 we have that x 2 ⌦M (u), where ⌦M (u) has
been specified in Definition 18.10, and ⌦M (u) ⇢ ⌦reg due to (18.19).
The second part of the statement can be achieved noticing that, in the case when r2 F
is uniformly continuous, no bound on |(ru)Br (x) | is needed in the decay theorem and in
the characterization of the regular set. ⇤
2,2
Proposition 18.20. If u 2 Wloc (⌦), we have that
Hn 2
(⌃1 (u)) = 0 .
Proof. Let us employ Proposition 18.15 with the absolutely continuous measure µ :=
|r2 u|2 L n . Obviously we choose k = (n 2) and we have that µ vanishes on sets with
finite H n 2 -measure. The thesis follows when we observe that
[1 ⇢
µ(B r (x)) 1
⌃1 (u) = x 2 ⌦ lim sup n 2
> .
⌫=1
r!0 !n 2 r ⌫
⇤
By the second part of the statement of Proposition 18.19 we get:
Corollary 18.21. If we add the uniform continuity of D2 F to the hypotheses of Propo-
sition 18.20, we can conclude that
Hn 2
(⌃(u)) = 0 . (18.39)
The estimate on the Hausdor↵ measure of ⌃2 (u) is a bit more complex and passes
through the estimate of the Hausdor↵ measure of the so-called approximate discontinuity
set Sv of a function v.
Definition 18.22. Given a function v 2 L1loc (⌦), we put
⇢ Z
⌦ \ Sv := x 2 ⌦ 9 z 2 R s.t. lim |v(y) z| dy = 0 .
r#0 Br (x)
When such a z exists, it is unique and we will call it approximate limit of v at the point
x.
Theorem 18.23. If v 2 W 1,p (⌦), 1 p n, then
Hn p+"
(Sv ) = 0 8" > 0 .
121
Notice that the statement is trivial in the case p > n, by the Sobolev Embedding
Theorem (i.e. Sv = ;): as p increases the Hausdor↵ dimension of the approximate
discontinuity set moves from n 1 to 0.
Applying this theorem to v = ru 2 H 1,2 (⌦; Rm⇥n ), p = 2, we get that H n 2+" (⌃2 (u)) =
0.
Proof. (1) Fix 0 < ⌘ < ⇢, we claim that
Z ⇢ Z Z
n!n (v)B⌘ (x) (v)B⇢ (x) (n 1) t n |rv(y)| dy dt + ⇢ (n 1)
|rv(y)| dy ;
0 Bt (x) B⇢ (x)
(18.40)
we will show this in the part (3) of this Rproof.
Suppose that x is a point for which Bt (x) |rv(y)| dy = o(tn 1+" ) for some " > 0, then
R
it is also true that ⇢ (n 1) B⇢ (x) |rv(y)| dy ! 0 and the sequence (v)Br (x) admits a limit
z as r ! 0 because it is a Cauchy sequence. Thanks to the Poincaré inequality
Z Z
(n 1) r!0
|v(y) (v)Br (x) | dy C(n)r |rv(y)| dy ! 0 ,
Br (x) Br (x)
therefore Z
r!0
|v(y) z| dy ! 0 ,
Br (x)
that is to say, x 2
/ Sv . This chain of implications means that, for all " > 0,
⇢ Z
⌦ \ Sv x2⌦ |rv(y)| dy = o(tn 1+" ) . (18.41)
Bt (x)
122
(3) This third part is devoted to the proof of (18.40); for the sake of simplicity we put
x = 0. Let us consider the characteristic function B1 ; since we would like to di↵erentiate
the map Z ✓ ◆
n y
⇢ 7! ⇢ v(y) dy ,
⇢
a possible proof of (18.40) is based on a regularization of , di↵erentiation and passage
to the limit.
We produce instead a direct proof based on a ad hoc calibration: we need a vector
field with supp ⇢ B ⇢ whose divergence almost coincides with the operator acting on
v in left member of (18.40), that is
n n
div = n ⌘ B⌘ ⇢ B⇢ . (18.43)
Therefore,
n n n n
(x) := x ⌘ ^ |x| ⇢ ^ |x|
verifies (18.43) and, with the notation µ = |rv| B⇢ L n , there holds
Z Z Z
n n
v(y) dy v(y) dy = v(y)div (y) dy (18.44)
⌘ n B⌘ ⇢n B⇢
Z Z Z
= (y) · rv(y) dy | (y)||rv(y)| dy |y| (n 1) dµ(y) (18.45)
B⇢ R n
Z 1 Z 1
(n 1)
= µ |y| > t dt = (n 1) s n µ(Bs ) ds (18.46)
0 0
Z ⇢ Z Z 1 Z
n n
= (n 1) s |rv(y)| dy ds + (n 1) s |rv(y)| dy ds
0 Bs ⇢ B⇢
Z ⇢ Z Z
n (n 1)
= (n 1) s |rv(y)| dy ds + ⇢ |rv(y)| dy ,
0 Bs B⇢
where we pass from (18.44) to (18.45) by the divergence theorem, from (18.45) to (18.46)
by Cavalieri’s principle and then it is all change of variables and Fubini’s theorem. ⇤
Remark 18.24. In the case p = 1 it is even possible to prove that Sv is -finite with
respect to H n 1 , so the measurement of the discontinuity set with the scale of Hausdor↵
measures is sharp. On the contrary, in the case p > 1 the right scale for the measurement
of the approximate discontinuity set are the so-called capacities.
123
19 Some tools from convex and nonsmooth analysis
19.1 Subdi↵erential of a convex function
In this section we briefly recall some classical notions and results from convex and nons-
mooth analysis, which will be useful in dealing with uniqueness and regularity results for
viscosity solutions to partial di↵erential equations.
In the sequel we consider a convex open subset ⌦ of Rn and a convex function u : ⌦ !
R. Recall that u is convex if
If u 2 C 2 (⌦) this is equivalent to say that r2 u(x) 0, in the sense of symmetric operators,
for all x 2 ⌦.
Definition 19.1 (Subdi↵erential). For each x 2 ⌦, the subdi↵erential @u(x) is the set
u(x + t0 v) u(x) (1 t0 /t) u(x) + (t0 /t)u(x + tv) u(x) u(x + tv) u(x)
= ,
t0 t0 t
(19.2)
for any 0 < t0 < t. Hence, for every y 2 ⌦, we have (choosing t = 1, v = y x)
124
(i) The graph of the subdi↵erential, i.e. {(x, p)|p 2 @u(x)} ⇢ ⌦ ⇥ Rn , is closed, in
fact convex functions are continuous (suffices, by (ii) below, to show that they are
locally bounded to obtain even local Lipschitz continuity).
(ii) Convex functions are locally Lipschitz in ⌦; to see this, fix a point x0 2 ⌦ and
x, y 2 Br (x0 ) b BR (x0 ) b ⌦. Thanks to the monotonicity of di↵erence quotients
seen in (19.2), we can estimate
where yR 2 @BR (x0 ) is on the halfline starting from x and containing y. Reversing
the roles of x and y we get
osc(u, B R (x0 ))
Lip(u, Br (x0 )) .
R r
This proves the local Lipschitz continuity and we can use this information to replace
B r (x0 ) by BR (x0 ), or even Br (x0 ) by B r (x0 ) in the inequality above. Equivalently
osc(u, BR (x0 ))
ess sup |ru| ,
Br (x0 ) R r
because of (1.6).
p 2 @f (x) .
125
Theorem 19.4 (Nonsmooth mean value theorem). Consider a convex function f : ⌦ ! R
and a couple of points x, y 2 ⌦. There exist z in the closed segment between x and y and
p 2 @f (z) such that
f (x) f (y) = hp, x yi .
Proof. Choose a positive convolution kernel ⇢ with support contained in B 1 and define
the sequence of functions f" := f ⇤ ⇢" , which are easily seen to be convex in the set ⌦" in
(1.3), because
Z
f" ((1 t)x + ty) = f ((1 t)x + ty "⇠)⇢(⇠) d⇠
Z ⌦
moreover f" ! f locally uniformly. Thanks to the classical mean value theorem for
regular functions, for every " > 0 there exists z" = (1 ✓" )x + ✓" y, with ✓" 2 (0, 1), such
that
f" (x) f" (y) = hp" , x yi .
with p" = rf" (z" ) 2 @f" (z" ). Since (z" , p" ) are uniformly bounded as " ! 0, we can find
"k ! 0 with ✓"k ! ✓ 2 [0, 1] and p"k ! p. Remark 19.3(iv) allows us to conclude that
p 2 @f ((1 ✓)x + ✓y) and
f (x) f (y) = hp, x yi .
⇤
As an application of the nonsmooth mean value theorem, we can derive a pointwise
version of Remark 19.3(iii). Notice that we will follow a similar idea to achieve second
order di↵erentiability.
126
Remark 19.6. Recall that a continuous function f : ⌦ ! R is convex if and only if its
Hessian r2 f is non-negative, i.e. for every non-negative ' 2 Cc1 (⌦) and every ⇠ 2 Rn
there holds Z
@ 2'
f (x) 2 (x) dx 0 .
⌦ @⇠
This result is easily obtained by approximation by convolution, because, still in the weak
sense,
r2 (f ⇤ ⇢" ) = r2 f ⇤ ⇢" .
Although we shall not need this fact in the sequel, except in Remark 19.17, let us
mention, for completeness, that the positivity condition on the weak derivative r2 f im-
plies that this derivative is representable by a symmetric matrix-valued measure. To see
this, it suffices to apply the following result to the second derivatives r2⇠⇠ f :
Lemma 19.7. Consider a positive distribution T 2 D 0 (⌦), i.e.
Proof. Fix an open set ⌦0 b ⌦, define K := ⌦0 and choose a non-negative cut-o↵ function
' 2 Cc1 (⌦) with '|K ⌘ 1. For every test function 2 Cc1 (⌦0 ), since (k kL1 ' ) 0
and T is a positive distribution, we have
where C(⌦0 ) := hT, 'i. Replacing by , the same estimate holds with |hT, i| in the
left hand side. By Riesz representation theorem we obtain the existence of µ. ⇤
for every non-negative ' 2 Cc1 (⌦) and for every ⇠ 2 Rn (in short r2 f I). We say
also that
• f is uniformly convex if > 0;
• f is semiconvex if 0.
127
Notice that, with the notation of Definition 19.8, a function f is -convex if and only
if f (x) |x|2 /2 is convex.
Analogous concepts can be given in the concave case, namely -concavity (i.e. r2 f
I), uniform concavity, semiconcavity. An important class of semiconcave functions is
given by squared distance functions:
Example 19.9. Given a closed set E ⇢ Rn , the square of the distance from E is 2-
concave. In fact,
dist2 (x, E) |x|2 = inf (x y)2 |x|2 = inf |y|2 2hx, yi ; (19.3)
y2E y2E
since the functions x 7! |y|2 2hx, yi are affine, their infimum over y 2 E, that is (19.3),
is concave.
Particularly in the duality theory of convex functions, it is useful to extend the concept
and convexity to functions f : Rn ! R[{+1}. The concept of subdi↵erential at points x
where f (x) < 1, extends immediately and, in the interior of the convex set {f < 1}, we
recover all the properties stated before (mean value theorem, local Lipschitz continuity).
Conversely, given f : ⌦ ! R convex with ⌦ convex, a canonical extension f˜ of f to the
whole of Rn is n o
f˜(x) := inf lim inf f (xh ) : xh 2 ⌦, xh ! x .
h!1
hp q, x yi 0 8p 2 @f (x), 8q 2 @f (y).
Proof. It is sufficient to sum the inequalities satisfied, respectively, by p and q, i.e.
hp q, x yi |x y|2 . (19.4)
128
(ii) If > 0, for every p 2 Rn no more than one x 2 {f < 1} can satisfy p 2 @f (x),
because, through (19.4), we get
In particular, setting [
L := @f (x) ,
f (x)<1
there exists a single-valued and onto map (@f ) 1 : L ! {x : @f (x) 6= ;} such that
p 2 @f ((@f ) 1 (p)). In addition, L = Rn : given p, to find x such that p 2 @f (x) it
suffices to minimize y 7! f (y) hp, yi and to take x as the (unique) minimum point.
1 1
(iii) Moreover, (@f ) is a Lipschitz map: rewriting (19.4) for (@f ) we get
thus Lip((@f ) 1 ) 1/ .
The conjugate of a function f : Rn ! R [ {+1}, not identically equal to +1, is
defined as
f ⇤ (x⇤ ) := sup hx⇤ , xi f (x) ;
x2Rn
we immediately point out that f ⇤ is convex and lower semicontinuous, because it is the
supremum of a family of affine functions. The assumption that f (x) < 1 for at least one
x ensures that f ⇤ : Rn ! R [ {+1}. Equivalently, f ⇤ is the smallest function satisfying
129
order-preserving and coincides, as it is easily seen, with the identity on affine functions
`(x) = hp, xi+c (notice that `⇤ is finite only at x⇤ = p and `⇤ (p) = c). Since convex lower
semicontinuous functions are supremum of affine functions (again as an application of the
Hahn-Banach theorem), these two facts yield (f ⇤ )⇤ f on convex lower semicontinuous
functions, completing the proof. ⇤
A byproduct of the previous proof is that (f ⇤ )⇤ = f in the class of convex and lower
semicontinuous functions f : Rn ! R [ {+1}, not identically equal to +1. This way
(19.5) becomes completely symmetric and it is easily seen that (19.6) gives
In particular, in the case when f is -convex for some > 0, from the quadratic
growth of f we obtain that f ⇤ is finite and that @f ⇤ = (@f ) 1 is single-valued and
Lipschitz, therefore f ⇤ 2 C 1,1 (Rn ).
Defining
1
fr (y) :=(f (x0 + ry) f (x0 ))
r
and noticing that rfr (y) = rf (x0 + ry) (still in the distributional sense), we are able to
rewrite (19.8) as Z
r!0
|rfr (y) rf (x0 )| dy ! 0,
B1 (0)
where (fr ) is a sequence of functions with equibounded Lipschitz constant and fr (0) = 0
for every r > 0. Thanks to the Ascoli-Arzelà theorem, as r # 0, this family of functions
130
has limit points in the uniform topology. Any limit point g obviously satisfies g(0) = 0,
and since rg is a limit point of rfr in the weak⇤ topology, the strong convergence of rfr
to rf (x0 ) gives rg ⌘ rf (x0 ), still in the weak sense. We conclude that g(x) = rf (x0 )x,
so that g is uniquely determined and
1 r!0
fr (y) = (f (x0 + ry) f (x0 )) ! rf (x0 )y
r
uniformly in B 1 (0). This convergence property is immediately seen to be equivalent to
the classical di↵erentiability of f at x0 , with gradient equal to rf (x0 ). ⇤
The proof of the following classical result can be found, for instance, in [11] and [12].
is L n -measurable6 and
Z Z
| det rf (x)| dx = N (y, A) dy L n (f (A)) .
A Rn
131
We briefly recall, from Remark 19.11, that @f associates to each x 2 Rn the subdif-
ferential set, on the contrary is a single-valued map which associates to each p 2 Rn
the point x such that p 2 @f (x). Let us define the set of “bad” points
⌃ := {p | @r (p) or 9r (p) and det r (p) = 0} .
Since is Lipschitz, Rademacher Theorem 19.13 and the area formula 19.14 give
Z
n
L ( (⌃)) | det r | dp = 0 .
⌃
We shall prove that the stated di↵erentiability property holds at all points x 2 / (⌃).
Let us write x = (p) with p 2 / ⌃, so that rf (x) = p, there exists the derivative r (p)
and, since it is invertible, we can name
1
S(x) := (r (p)) .
If y = (q), we get
1
S(x) (q p S(x)(y x)) = (y x r (p)(q p))
= ( (q) (p) r (p)(q p))
= o(|p q|) = o(|x y|) .
Therefore
|q rf (x) S(x)(y x)|
lim =0. (19.9)
y!x
q2@f (y)
|y x|
The result got in (19.9), together with the nonsmooth mean value Theorem 19.4, give
us the second order expansion. In fact, let
1
f˜(y) := f (y) f (x) hrf (x), (y x)i hS(x)(y x), (y x)i .
2
Since
@ f˜(y) = @f (y) rf (x) S(x)(y x)
we can read (19.9) as lim |q|/|y x| = 0. Now, choose ✓ 2 [0, 1] and a vector
q2@ f˜(x), y!x
q 2 @ f˜((1 ✓)y + ✓x) such that f˜(y) = hq, y xi (since f˜(x) = 0) to find
f˜(y) = hq, y xi = o(|y x|2 ) .
By the very definition of f˜, the statement follows. ⇤
Remark 19.17 (Characterization of S). A blow-up analysis, analogous to the one per-
formed in the proof of Rademacher’s theorem, shows that the matrix S(x) in Alexandrov’s
theorem is the density of the measure r2 f with respect to L n , see [2] for details.
132
20 Viscosity solutions
20.1 Basic definitions
In this section we want to give the notion of viscosity solution for general equations having
the form
E(x, u(x), ru(x), r2 u(x)) = 0 (20.1)
where u is defined on some locally compact domain A ⇢ Rn . This topological assumptions
is actually very useful, because we can deal at the same time with open and closed domains,
and also domain of the form Rn 1 ⇥ [0, 1), which typically occur in parabolic problems.
We first need to recall two classical ways to regularize a function.
Definition 20.1 (u.s.c. and l.s.c. regularizations). Let A0 ⇢ A be a dense subset and
u : A0 ! R. We define its upper regularization u⇤ on A by one of the following equivalent
formulas:
⇢
u (x) := sup lim sup u(xh ) | (xh ) ⇢ A0 , xh ! x
⇤
h
= inf sup u
r>0 B (x)\A0
r
133
(ii) for any x 2 A, if ' is C 1 in a neighbourhood of x and u⇤ ' has a local maximum
at x, then
E⇤ (x, u⇤ (x), r'(x), r2 '(x)) 0 . (20.2)
It is obvious from the definition that the property of being a subsolution is invariant
under u.s.c. regularization, i.e. u is a subsolution if and only if u⇤ is a subsolution.
The geometric idea in this definition is to use a local comparison principle, since
assuming that u⇤ ' has a maximum at x implies, if u is smooth, that ru⇤ (x) = r'(x)
and r2 u⇤ (x) r2 '(x). So, while in the classical theory of PDEs an integration by parts
formula allows to transfer derivatives from u to the test function ', here the comparison
principle allows to transfer (to some extent, since only an inequality holds for second order
derivatives) the derivatives from u to the test function '.
Similarly, we give the following:
Definition 20.4 (Supersolution). A function u : A ! R is a supersolution for the
equation (20.1) (and we write E 0) if the two following conditions hold:
(i) u⇤ is a real-valued function;
(ii) for any x 2 A, if ' is C 1 in a neighbourhood of x and u⇤ ' has a local minimum
at x, then
E ⇤ (x, u⇤ (x), r'(x), r2 '(x)) 0. (20.3)
This example shows that some continuity assumption is needed, in order to hope for
reasonable existence and uniqueness results.
Remark 20.7. Rather surprisingly, a solution of E = 0 in the viscosity sense does not
necessarily solve E = 0 in the viscosity sense. To show this, consider the equations
|f 0 | 1 = 0 and 1 |f 0 | = 0 and the function f (t) = min {1 t, 1 + t} . In this case, it
is immediate to see that f is a subsolution of the first problem (and actually a solution,
as we will see), but it is not a subsolution of the second problem, since we can choose
134
identically ' = 1 to find that the condition 1 |'0 (0)| 0, corresponding to (20.2), is
violated.
We have instead the following parity properties:
(a) Let E be odd in (u, p, S). If u verifies E 0, then u verifies E 0.
(b) Let E be even in (u, p, S). If u verifies E 0, then u verifies E 0.
We now spend some words on the ways of simplifying the conditions that have to be
checked in order prove the subsolution or supersolution property. We just examine the
case of subsolutions, the case of supersolutions being the same (with obvious variants).
We have already seen in Remark 20.5 that we can assume without loss of generality
that u⇤ ' has a strict local maximum, equal to 0, at x. We can also work equivalently with
the larger class of C 2 functions ', in a neighbourhood of x. One implication is trivial, let us
see the converse one. Let ' 2 C 2 and assume u⇤ (y) '(y) 0 for y 2 B r (x), with equality
only when y = x. By appropriate mollifiers, we can build a sequence ('k ) ⇢ C 1 (B r (x))
with 'k ! ' in C 2 (B r (x)). Let then xk be a maximum in B r (x) of the function u⇤ 'k .
Since 'k ! ' uniformly, it is easy to see that any limit point of (xk ) has to be a maximum
for u⇤ ', hence it must be x; in addition the convergence of the maximal values yields
u⇤ (xk ) ! u⇤ (x). The subsolution property, applied with 'k at xk , gives
E⇤ (xk , u⇤ (xk ), r'k (xk ), r2 'k (xk )) 0
and we can now let k ! 1 and use the lower semicontinuity of E⇤ to get the thesis.
Actually, it is rather easy now to see that the subsolution property is even equivalent
to
E⇤ (x, u⇤ (x), p, S) 0 8 (p, S) 2 J2+ u⇤ (x)
where J2+ u⇤ is the second-order super jet of u, namely
J2+ u⇤ (x) := (p, S) u⇤ (y) u⇤ (x) + hp, y xi + 12 hS(y x), y xi + o(|y x|2 ) .
Indeed, let P (y) := u⇤ (x)+hp, y xi+ 12 hS(y x), y xi, so that u⇤ (y) P (y)+o(|y x|2 ),
with equality when y = x. Hence, for any " > 0 we have u⇤ (y) P (y) + "|y x|2 on a
sufficiently small neighbourhood of x with equality at y = x and we can apply (20.2) to
this smooth function to get
E⇤ (x, u⇤ (x), p, S + 2"I) = E⇤ (x, u⇤ (x), rP (x), r2 P (x) + 2"I) 0
and by lower semicontinuity we can let " ! 0 and prove the claim. Of course, if we are
dealing with first order equations, only the first order super jet is needed.
Remark 20.8. After these preliminary facts, it should be clear that this theory, despite
its elegance, has two main restrictions: on the one hand it is only suited to first or second
order equations (since no information on third derivatives comes from local comparison),
on the other hand it cannot be generalized to vector-valued functions.
135
20.2 Viscosity versus classical solutions
We first observe that a classical solution is not always a viscosity solution. To see this,
consider on R the problem u00 2 = 0. The function f (t) = t2 is clearly a classical solution,
but it is not a viscosity solution, because it is not a viscosity supersolution (take ' ⌘ 0
and study the situation at the origin).
Since we can always take u = ' if u is at least C 2 , the following theorem is trivial:
136
Proof. Assume as usual that u⇤ ' has a strict local maximum at x, equal to 0, and
denote by K the compact set B r (x) \ A for some r to be chosen sufficiently small, so that
x is the unique maximum of u ⇤ ' on K.
By a diagonal argument can find a sequence (xh ) inside K, convergent to x, and a
sequence of functions (vh ) ⇢ F such that u⇤ (x) = limh u(xh ) = limh vh (xh ). Hence, if we
call yh the maximum of vh⇤ ' on K, then
u⇤ (yh ) '(yh ) vh⇤ (yh ) '(yh ) vh⇤ (xh ) '(xh ) vh (xh ) '(xh ).
Since by our construction we have vh (xh ) '(xh ) ! 0 for h ! 1, we get that every limit
point y of (yh ) satisfies
u⇤ (y) '(y) 0.
Hence y is a maximum in K of u⇤ ', u⇤ (y) '(y) = 0 and y must coincide with x.
Consequently yh ! x, lim suph (u⇤ (yh ) '(yh ) u⇤ (x) '(x) and, by comparison, the
same is true for the intermediate terms, so that vh⇤ (yh ) ! u⇤ (x). In order to conclude, we
just need to consider the viscosity condition at the points yh , which reads
E⇤ (yh , vh⇤ (yh ), r'(yh ), r2 '(yh )) 0 ,
and let h ! 1 to get
E⇤ (x, u⇤ (x), r'(x), r2 '(x)) 0.
⇤
We can now state a first existence result.
Theorem 20.12 (Perron). Let f and g be respectively a subsolution and a supersolution
of (20.1), such that f⇤ > 1 and g ⇤ < +1 on A. If f g on A and the functions
E⇤ (x, u, p, ·) and E ⇤ (x, u, p, ·) are non-increasing, then there exists a solution u of (20.1)
satisfying f u g.
Proof. Call
F := {v | v is a subsolution of (20.1) and v g } .
We know that f 2 F, so that this set is not empty. Hence, we can define u :=
sup {v| v 2 F} . By our definition of F, we have that u g and therefore u⇤ g ⇤ < +1.
Since u⇤ u⇤ f⇤ > 1, in A, by Theorem 20.11 u is a subsolution on A. Consequently,
we just need to prove that it is also a supersolution on the same domain.
Pick a test function ' such that u⇤ ' has a relative minimum, equal to 0, at x0 .
Without loss of generality, we can assume that
u⇤ (x) '(x) |x x0 | 4 on A \ B r (x) (20.4)
for some sufficiently small r > 0. Assume by contradiction that
E ⇤ (x0 , u⇤ (x0 ), r'(x0 ), r2 '(x0 )) < 0 (20.5)
4
and define a function w := max{' + , u} for some parameter > 0. We claim that:
137
(a) w is a subsolution of (20.1);
(b) {w > u} =
6 ;;
E ⇤ (x, '(x) + 4
, r'(x), r2 '(x)) 0 on B 2 (x0 ) \ A .
This means that '+ 4 is a classical subsolution of (20.1) on this domain and hence, by our
monotonicity hypothesis, it has to be also a viscosity subsolution. Consequently, by a very
special case of the previous theorem, we get that the function w is a viscosity subsolution
of (20.1) on B 2 (x0 )\A. Moreover, by (20.4), we know that w = u on (A \ Br (x))\B (x0 ).
Since the notions of viscosity subsolution and supersolution are clearly local, w is a global
subsolution on A.7
To prove that {w > u} = 6 ; we just need to observe that, for any > 0, u⇤ (x0 ) =
4
'(x0 ) < '(x0 ) + , and on any sequence (xh ) such that u(xh ) ! u⇤ (x0 ), we must have
for h sufficiently large the inequality u(xh ) < '(xh ) + 4 .
Finally, we have to show that w g: this completes the proof of the claim and gives the
desired contradiction. To this aim, it is enough to prove that there exists > 0 such that
' + 4 g on A \ B (x0 ). But this readily follows, by an elementary argument, showing
that '(x0 ) = u⇤ (x0 ) < g⇤ (x0 ). Again, assume by contradiction that u⇤ (x0 ) = g⇤ (x0 ) : if
this were the case, the function g⇤ ' would have a local minimum at x0 and so, since
g⇤ is a viscosity supersolution, we would get
138
that, for any y 2 C, the function x 7! |x y| is a classical supersolution in A (because
y2/ A) and hence a viscosity supersolution of our problem.
The fact that u is also a subsolution follows by the general implication:
Indeed, let x be a local maximum for f ', so that f (y) '(y) f (x) '(x)
for any y 2 Br (x) (and r small enough). This is equivalent, on the same domain, to
'(y) '(x) f (y) f (x) |y x| and, by the Taylor expansion, we finally get
for f continuous (or at least upper semicontinuous), which is proved by means of the
regularizations f " (x) := supy f (y) |x y|2 /" that we will study more in detail later
on. We just sketch here the structure of the argument:
(2) |rf " |2 1 0 pointwise L n -a.e., because f " is semiconcave, hence locally Lipschitz,
and therefore the inequality holds at any di↵erentiaiblity point by the super-jet
characterization of viscosity subsolutions;
139
It remains to show that w u. Assume first that A is bounded: we will show later on
that this is not restrictive. By contradiction, assume that w(x0 ) > u(x0 ) for some x0 ; in
this case there exist 0 > 0 and 0 > 0 such that
⇢
1
sup w(x) (1 + )u(y) |x y|2 0
x,y 2"
for all " > 0 and 2 (0, 0 ). Indeed, it suffices to bound from below the supremum with
w(x0 ) (1 + )u(x0 ), which is larger than 0 := (w(x0 ) u(x0 ))/2 for > 0 small enough.
Moreover, for " > 0 and 2 (0, 0 ), the supremum is actually a maximum because it
is clear that we can localize x in A (otherwise the whole sum is non-positive) and y in a
bounded set of Rn (because w is bounded on A, and again for |y x| large the whole sum
is non-positive). So, call (x, y) a maximizing couple, omitting for notational simplicity the
dependence on the parameters ", . The function x 7! w(x) 2" 1
|x y|2 has a maximum
at x = x and so we can exploit the fact that w(·) is a viscosity solution of our equation
(with respect to the test function '(x) = |x y|2 /(2")) to derive |r'|2 (x) 1, that is
|x y|
1.
"
We also claim that necessarily y 2 A, if " is sufficiently small, precisely " < 0 . Indeed,
assume by contradiction that y 2/ A, so that w(y) = 0, then by the triangle inequality
1 1
0 w(x) |x y|2 |x y| |x y|2 |x y| .
2" 2"
As a consequence, we get 0 |x y| ", which gives a contradiction.
Now, choosing " > 0 so that y 2 A, the function y 7! (1 + )u(y) + 1
2"
|x y|2 has a
minimum at y = y and arguing as above we obtain
x y
(1 + ) ,
"
which is not compatible with |x y| ". Hence, at least when A is bounded, we have
proved that w = u.
In the general case, fix a constant R > 0 and define uR (x) := u(x) ^ dist(x, Rn \ B R ) :
this is a supersolution of our problem on A \ BR , since u(x) is a supersolution on A
and dist(x, Rn \ B R ) is a supersolution on BR (by the infimum property). Moreover,
Lip(uR ) 1 implies that uR is a global subsolution and we can apply the previous result
(special case) to the function uR to get
uR (x) = d(x, Rn \ (A \ BR )).
Letting R ! 1 we first exclude C = ; since in that case uR " 1 which is not admissible
since uR u and then (by C 6= ;) we obtain u(x) = dist(x, C). ⇤
140
Remark 20.14. We can also give a di↵erent interpretation of the result above. In the
spirit of the classical Liouville’s theorems we can say that “the equation |ru|2 1 = 0 does
not have entire viscosity solutions on Rn that are bounded from below”. Nevertheless,
there exist trivial examples of functions that solve this equation in the viscosity sense and
are unbounded from below (e.g. take u(x) = xi for some i 2 {1, . . . , n}).
On the other hand, f (resp. g) is also a regular subsolution (resp. supersolution) so that
Hence, if we compare (20.6) with (20.7), we find a contradiction as soon as one of the two
inequalities in (20.7) is strict.
In order to hope for a comparison principle, this argument shows the necessity to
approximate subsolutions (or supersolutions) with strict subsolutions, and this is always
linked to some form of strict monotonicity of the equation, variable from case to case (of
course in the trivial case F ⌘ 0 no comparison principle is possible). To clarify this point,
let us consider the following example. Consider the space-time coordinates x = (y, t) and
a parabolic problem
F (ry,t u, r2y,t u) = @t u G(r2y u)
with G non-decreasing, in the appropriate sense. In this case, we can reduce ourselves to
strict inequalities by performing the transformation u e t u.
In order to get a general uniqueness result for viscosity solution, we cannot just argue
as in the case of the distance function and we need to follow a strategy introduced by
Jensen. The first step is to obtain a refined versions of the maximum principle. We start
with an elementary observation.
141
Remark 20.15. If (p, S) 2 J2 u(x) and u has a relative maximum at x, then necessarily
p = 0 and S 0. To see this, it is enough to apply the definitions: by our two hypotheses
1
0 u(y) u(x) hp, y xi + hS(y x), y xi + o(|y x|2 )
2
and hence
y x
hp, i o(|y x|) ) p = 0 ,
|y x|
hS(y x), y xi
o(1) ) S 0 .
|y x|2
We are now ready to state and prove Jensen’s maximum principle for semiconvex
functions.
Theorem 20.16 (Jensen’s maximum principle). Let u : ⌦ ! R be semiconvex and let
x0 2 ⌦ a local maximum for u. Then, there exist a sequence (xk ) converging to x0 and
"k # 0 such that u is pointwise second order di↵erentiable at xk and
The proof is based on the following lemma. In the sequel we shall denote by sc(u, ⌦)
the least nonnegative constant C such that u is ( C)-convex, i.e. u + C|x|2 /2 is convex
(recall Definition 19.8).
Theorem 20.17. Let B ⇢ Rn be a ball of radius R centered at the origin and u 2 C(B)
semiconvex, with8
max u > max u .
B @B
Then, if we let
it must be
!n n
L n (G ) (20.8)
[sc(u, B)]n
for 0 < < (maxB u minB u) /(2R).
Proof. We assume first that u is also in C 1 (B). Pick a > 0, so small that 2R <
maxB u max@B u, and consider a perturbation u(y) + hp, yi with |p| . We claim that
such function necessarily attains its maximum in B. Indeed, this immediately comes from
the two inequalities
max (u + hp, yi) max u + R
@B @B
8
Notice that this implies sc(u, B) > 0, since maxB = max@B for convex functions.
142
and
max(u + hp, yi) max u R.
B B
Consequently, there exists x 2 B such that ru(x) = p. This shows that ru(G ) = B .
To go further, we need the area formula. In this case, it gives
Z Z
2
| det r u| dx = card ({x | ru(x) = p }) dp !n n
G B
because the points in G are maxima for the function u(y)+hp, yi : this implies r2 u(x) 0
for any x 2 G and, by semiconvexity, r2 u(x) sc(u, B)I. If we combine these two
inequalities, we get (20.8).
In the general case we argue by approximation, finding radii rh " R and smooth functions
uh in B rh such that uh ! u locally uniformly in B and lim suph sc(uh , Brh ) sc(u, B); to
conclude, it suffices to notice that any limit of points in G (uh ) \ Brh belongs to G (u),
hence L n (G (u)) lim suph L n (G (uh ) \ Brh ). ⇤
We can now prove Jensen’s maximum principle. As a preliminary remark, observe
that, in Definition 19.8 one has (for our u) = 0 then the claim is trivial, so that we can
without loss of generality assume that < 0 and Theorem 20.17 applies.
Proof. Let x0 be a local maximum of u. We can choose R > 0 sufficiently small so
that u u(x0 ) in B R (x0 ) and, without loss of generality, we can assume u(x0 ) = 0.
This becomes a strict local maximum for the function u e(x) = u(x) |x x0 |4 . It is also
e is semiconvex in B R (x0 ). We now apply Theorem 20.17 to u
easy to verify that u e: for any
= 1/k with k large enough we obtain that L n (G1/k ) > 0 and (thanks to the Alexandrov
theorem) this means that there exists a sequence of points (xk ) such that u e is pointwise
second order di↵erentiable at xk and, for appropriate vectors pk with |pk | 1/k, the
function u e(y) hpk , yi has a local maximum at xk . Since |pk | ! 0, any limit point of
(xk ) for k ! 1 has to be a local maximum for u e, but in B R (x0 ) this necessarily implies
2
xk ! x0 . Moreover pk = re u(xk ) ! 0 and r u e(xk ) 0. As a consequence
ru(xk ) = re
u(xk ) + 4|xk x0 |2 (xk x0 ) ! 0
and the identity
r2 |z|4 = 4|z|2 I + 8z ⌦ z (20.9)
gives
r2 u(xk ) = r2 u
e(xk ) + 8(xk x0 ) ⌦ (xk x0 ) + 4|xk x0 | 2 I
2
ru e(xk ) + 12|xk x0 |2 I .
Setting "k = 12|xk x0 |2 we get the thesis. ⇤
143
We now introduce another important tool in the theory of viscosity solutions.
Definition 20.18 (Inf and sup-convolutions). Given u : A ! R and a parameter " > 0,
we can build the regularized functions
⇢
" 1
u (x) := sup u(y) |x y|2 (20.10)
y2A "
(ii) u" u and u" # u pointwise in A. If u is continuous, then u" # u locally uniformly;
and the functions in the right hand side are affine with respect to x. It follows that the
left hand side is convex, which means sc(u" , Rn ) 2/".
(ii) The inequality u" u and the monotonicity in " are trivial. In addition, we can take
quasi-maxima (y" ) satisfying
2 2 2
" " "
u" (x) u(y" ) + " K(1 + |y" |) + " K(1 + |x| + | " |) +" .
" " "
144
with " = |y" x|. Via these two inequalities, one first sees that y" ! x so that, exploiting
the upper semicontinuity of u and neglecting the quadratic term in the first inequality we
get
u(x) lim sup u(y" ) lim sup u" (x) .
"!0 "!0
If u is continuous, the claim comes from Dini’s monotone convergence theorem and the
local compactness of A.
(iii) Let x0 2 A" and let y0 2 A be a corresponding maximum, so that u" (x0 ) = u(y0 )
|x0 y0 |2 /". Let then ' be a smooth function such that u" ' has a local maximum in
x0 and, without loss of generality, we can take u" (x0 ) = '(x0 ). Let us call r the radius
such that u" ' on Br (x0 ).
Define (x) := '(x y0 + x0 ) : we claim that u has a local maximum at y0 with value
|x0 y0 |2 /". If we prove this claim, then it must be
F (r (y0 ), r2 (y0 )) 0
F (r'(x0 ), r2 '(x0 )) 0.
Remark 20.20. We will also need an x-dependent version of the previous result, that
reads as follows: if F (x, ru, r2 u) 0 in the sense of viscosity solutions on A, then for
all > 0 there holds F (x, ru" , r2 u" ) 0 on A" , where
145
20.5 Existence and uniqueness results
In this section we will collect some existence and uniqueness results for second order
equations. The main tool is the comparison principle, stated below. Throughout the
section we shall always assume that A is a bounded open set in Rn .
Proposition 20.21 (Comparison principle). Let F : A ⇥ Symn⇥n ! R be continuous and
satisfying, for some > 0, the strict monotonicity condition
F (x, S + tI) F (x, S) + t 8t 0
and the uniform continuity assumption
F (·, S), S 2 Symn⇥n , are equi-continuous in A.
Let u, u : A ! R be respectively a bounded u.s.c. subsolution and a bounded l.s.c. super-
solution to F (x, r2 u) = 0 in A, with (u)⇤ (u)⇤ on @A. Then u u on A.
Notice that the uniform continuity assumption, though restrictive, covers equations of
the form G(r2 u) + f (x) with f continuous in A.
A direct consequence of the comparison principle (take u = u = u) is the following
uniqueness result:
Theorem 20.22 (Uniqueness of continuous solutions). Let F be as in Proposition 20.21
and h 2 C(@A). Then the problem
8
< F (x, r2 u(x)) = 0 in A;
(20.13)
:
u=h on @A
146
Lemma 20.24. Let F, u and u be as in Proposition 20.21 and set
v , := u + |x|2 .
2
Hence:
(iii) if the comparison principle holds for v , for any > ( , A), that is
v, u on A, 8 > ( , A) , (20.14)
then u u on A.
Proof. Statements (i) follows by the translation invariance w.r.t. u of the equation, and
by r2 v = r2 u + I. Statement (ii) follows by the fact that u < u on @A.
If (20.14) holds, then
u v, u on A ,
and the comparison principle for u follows letting # 0, which allows to choose arbitrarily
small in view of (ii). ⇤
Proof. (of Proposition 20.21) Thanks to Lemma 20.24, without loss of generality we can
assume that u satisfies the stronger property
F (x, r2 u) 0
147
of u; since u" u and u" u we have
✓ ◆
" 1 4
max u (x) u" (y) |x y| u" (x0 ) u" (x0 ) u(x0 ) u(x0 ) = d0
A⇥A 4"
and we shall denote by (x" , y" ) 2 A ⇥ A a maximizing pair, so that
1
d0 + |x" y" |4 u" (x" ) u" (y" ) sup u inf u . (20.17)
4"
Also, we denote by x0" 2 A and y"0 2 A maximizers and minimizers respectively in
(20.15) and (20.16).
Now we claim that:
(a) lim inf dist(x" , @A) > 0 and lim inf dist(y" , @A) > 0;
"#0 "#0
(b) setting M = max{osc(i), osc(u)}, for " small enough, the supremum in (20.15) with
any x 2 A satisfying |x x" | < " is attained at a point x0 2 A with |x0 x|2 M "
and the infimum in (20.16) with any y 2 A satisfying |y y" | < " is attained at a
point y 0 2 A with |y 0 y|2 M ".
To prove (a), notice that, if (x̄, ȳ) is any limit point of (x" , y" ) as " # 0, then (20.17)
gives x̄ = ȳ and
✓ ◆
⇤ 0 0 |x" x0" |2 + |y" y"0 |2
d0 lim sup (u) (x" ) (u)⇤ (y" ) .
"#0 "
Since the supremum of (u)⇤ (u)⇤ is finite, this implies that |x" x0" | ! 0, |y" y"0 | ! 0,
hence (x0" , y"0 ) ! (x̄, x̄) as well and semicontinuity gives d0 (u)⇤ (x̄) (u)⇤ (x̄). By
assumption (u)⇤ (u)⇤ on @A, therefore x̄ 2 A and this proves (a).
To prove (b), it suffices to choose, thanks to (a), "0 > 0 and 0 > 0 small enough, so
that dist(x" , @A) 0 for " 2 (0, "0 ). In general, for x 2 A we have
1 0
u(x0 )|x x|2 u(x) u" (x)
"
which implies that the supremum in the definition
p of u" (x) is unchanged if we maximize
in the ball B x centered at x with radius M ". If |x x" | < ✏ and " < "0 , since
dist(x" , @A) 0 , this implies that the ball B x is contained in A for " small enough,
hence the supremum is attained. The argument for y" is similar.
Let us fix " small enough so that (b) holds and both x0" and y"0 belong to A, and let
us apply Jensen’s maximum principle to the (locally) semiconvex9 function
1
w(x, y) := u" (x) u" (y) |x y|4
4"
9
The local semiconvexity of w follows from Proposition 20.19.
148
to find zn := (x",n , y",n ) ! (x" , y" ) and n # 0 such that w is pointwise second order
di↵erentiable at zn , rw(zn ) ! 0 and r2 w(zn ) n I. By statement (b) and Remark 20.20,
for n large enough we have
On the other hand, the upper bound on r2 w(zn ) together with (20.9) give
⇢
r2 u" (x",n ) 2" (x",n y",n ) ⌦ (x",n y",n ) 1" |x",n y",n |2 I n I ,
(20.19)
r2 u" (y",n ) 2" (x",n y",n ) ⌦ (x",n y",n ) 1" |x",n y",n |2 I n I .
By (20.19) we obtain that r2 u" (x",n ) are uniformly bounded above, and they are also uni-
formly bounded below, since u" is semiconvex. Since similar remarks apply to r2 u" (y",n ),
we can assume with no loss of generality that r2 u" (x",n ) ! X" and r2 u" (y",n ) ! Y" . If
we now di↵erentiate w along a direction (⇠, ⇠) with ⇠ 2 Rn , we may use the fact that
along these directions the fourth order term is constant to get
Taking limits, this proves that X" Y" . On the other hand, from (20.18) we get
Hence
sup F (x, Y" ) inf F (y, Y" ) .
x2B pM " (x" ) y2B pM " (y" )
Since and are fixed positive constants independent of ", and since |x" y" | ! 0, this
contradicts the uniform continuity of F (·, S) for " sufficiently small. ⇤
149
for some c 2 R, p 2 Rn and S 2 Symn⇥n . We say that P is a paraboloid with opening
M 2 R if S = M I, namely
M 2
P (x) = c + hp, xi +
|x| .
2
It will be occasionally convenient to center a paraboloid P with opening M at some point
x0 , writing P (x) = P (x0 ) + hrP (x0 ), x x0 i + M2 |x x0 |2 .
Definition 20.25 (Tangent paraboloids). Given a function u : ⌦ ! R and a subset
A ⇢ ⌦ ⇢ Rn , we denote
✓(x0 , A, u) := inf {M |there exists P with opening M , u(x0 ) = P (x0 ) and u P on A } .
Moreover, we set
✓(x0 , A, u) := sup {M |there exists P with opening M , u(x0 ) = P (x0 ) and u P on A } ,
so that ✓(x0 , A, u) = ✓(x0 , A, u). Finally, denoting by ± the positive and negative
parts, we set n o
+
✓(x0 , A, u) := max ✓ (x0 , A, u), ✓ (x0 , A, u) 0.
Given a function u : ⌦ ! R and h > 0, let us consider the symmetric di↵erence
quotient in the direction ⇠ 2 Rn
2 u(x0 + h⇠) + u(x0 h⇠) 2u(x0 ) @ 2u
h,⇠ u(x0 ) := h,⇠ ( h,⇠ u)(x0 ) = ⇠ (x0 ) ,
h2 @⇠ 2
well defined if h|⇠| < dist(x0 , @⌦) and identically equal to M on paraboloids with open-
ing M . Notice that the symmetric di↵erence quotient satisfies, by applying twice the
integration by parts formula for h,⇠ ,
Z Z
2 2
u h,⇠ dx = h,⇠ u dx (20.20)
⌦ ⌦
whenever u 2 L1loc (⌦), 2 L1 (⌦) has compact support, |⇠| = 1 and the h-neighbourhood
of supp is contained in ⌦.
Remark 20.26 (Maximum principle for 2⇠ ). If a paraboloid P with opening M “touches”
u from above (i.e. P (x0 ) = u(x0 ) and P (x) u(x) in some ball Br (x0 )), then
2 2
h,⇠ u(x0 ) h,⇠ P (x0 ) =M whenever |⇠| = 1 and |h| r ,
and a similar property holds for paraboloids touching from below. Thus, passing to the
infimum from above and the supremum from below, we deduce the inequalities
2
✓(x0 , Br (x0 ), u) h,⇠ u(x0 ) ✓(x0 , Br (x0 ), u) whenever |⇠| = 1 and |h| r , (20.21)
and
2
| h,⇠ u(x0 )| ✓(x0 , Br (x0 ), u) whenever |⇠| = 1 and |h| r . (20.22)
150
Proposition 20.27. If u : ⌦ ! R satisfies
for some " > 0 and 1 < p 1, then u belongs to W 2,p (⌦) and, more precisely,
where we pass from the first to the second line with (20.20) and the inequality follows
from (20.22). Thanks to Riesz representation theorem, we know that the map ' 7!
R 2
⌦
u(x) @@⇠'2 (x) dx admits a representation with an element of Lp (⌦), which represents the
derivative r2⇠⇠ u in the sense of distributions and which satisfies (20.23). ⇤
In the space of n ⇥ n matrices we will consider the operator norm | · |L and, in the
subspace of symmetric matrices, the norm k · k provided by the largest modulus of the
eigenvalues in the spectrum (M ). Obviously these two norms coincide on Symn⇥n . From
(20.21) we get
kr2 u(x0 )k ✓(x0 , B" (x0 ), u) for all " > 0 (20.24)
at any point x0 where u has a second order Taylor expansion.
Corollary 20.29. If ⌦ ⇢ Rn is convex and ✓" 2 L1 (⌦) for some " > 0, then
151
because, when v is continuously di↵erentiable, there holds
Z 1 Z 1
v(x) v(y) = Dv((1 t)x + ty)(x y) dt |x y| |rv|L ((1 t)x + ty) dt .
0 0
152
Analogously, since bii ⇤ one has
n
X n
X
F (M + N ) F (M ) = tr(BN ) = bii ⇢i ⇤ ⇢i n⇤⇢max .
i=1 i=1
After this introductory part about definitions and notation, we enter in the core of
the matter of the Hölder regularity for viscosity solutions: as in De Giorgi’s work on the
XIX Hilbert problem, the regularity will be deduced only from inequalities derived from
ellipticity, without a specific attention to the original equation.
We will omit the dependence on and ⇤, when clear from the context.
As a matter of fact, denoting with (bij ) the coefficients of the matrix B 2 A ,⇤ in the
system of coordinates where M is diagonal, with M = diag(⇢1 , . . . , ⇢n ) we get
n
X X X
tr(BM ) = bii ⇢i ⇢i + ⇤ ⇢i (20.31)
i=1 ⇢i >0 ⇢i <0
153
(b) for every M, N it is simple to obtain from (20.29) and (20.30) that
M+ (M ) + M (N ) M+ (M + N ) M+ (M ) + M+ (N )
and, similarly,
M (M ) + M (N ) M (M + N ) M (M ) + M+ (N ) ;
(c) M± are elliptic (i.e., they satisfy (20.27)) with constants , n⇤, because of Exam-
ple 20.32 and (20.29), (20.30) which represent M± as an envelope of a family of
functionals with ellipticity constants , n⇤.
We also set
Sol ,⇤ (f ) := Sub /n,⇤ ( |f |) \ Sup ,n⇤ (|f |) . (20.33)
Remark 20.37. Roughly speaking, the classes defined above correspond to De Giorgi’s
classes DG± (⌦), since u being a solution to (20.26) with F having ellipticity constants
and ⇤ implies u 2 Sol ,⇤ (f ); thus, if we are able to infer regularity of functions in
Sol ,⇤ (f ) then we can “forget” thanks to Remark 20.35(d) the specific equation.
where M± are Pucci’s extremal operators, and we shall not emphasize from now on the
dependence on the ellipticity coefficients and ⇤. Notice that, since M+ M , the
intersection of the two sets can be nonempty.
154
The estimate we want to prove is named after Alexandrov, Bakelman and Pucci and is
therefore called ABP weak maximum principle. It plays the role in this regularity theory
played by the Caccioppoli inequality in the standard linear elliptic theory.
In the sequel we call “universal” a constant which depends only on the space dimension
n and on the ellipticity constants , ⇤.
Since f + measures, in some sense, how far u is from being concave, the estimate above
can be seen as a quantitative formulation of the fact that a concave function in a ball
attains its minimum on the boundary of the ball.
In order to prove the ABP estimate we set M := maxB r u and assume with no loss
of generality that M > 0.
The following facts are either trivial consequences of the definitions or easy applications
of the tools introduced in the convex analysis part: firstly M u 0, as a consequence
1,1
u 2 Wloc (B2r ) and finally since u is di↵erentiable a.e. by Rademacher’s theorem and
the graph of the subdi↵erential is closed, we get @ u (x) 6= ; for all x 2 B2r . We will use
this last property to provide a supporting hyperplane to u at any point in B r .
We need some preliminary results, here is the first one.
with c = c(n).
Proof. Let x1 2 Br be such that u (x1 ) = M . Fix ⇠ with |⇠| < M/(3r) and denote by L↵
the affine function L↵ (x) = ↵+hx, ⇠i. It is obvious that if ↵ 1, then the corresponding
hyperplane lies below the graph of u and there is a minimum value of ↵ such that this
happens, that is u L↵ on B 2r . The graph of u will then meet the corresponding
155
hyperplane at some point, say x0 2 B 2r . If it were |x0 | > r, then L↵ (x0 ) = 0, but on the
other hand |L↵ (x1 )| M and, since |x0 x1 | 3r, L↵ would have slope |⇠| M/3r,
which is a contradiction. Hence any contact point x1 must lie inside the ball Br ; from
u u L↵ we get r u (x1 ) = ⇠ and therefore BM/(3r) ⇢ r u (Br ). If we measure
the corresponding volumes and use the area formula, we get
✓ ◆n Z
M
!n det r2 u dx
3r Br
or, equivalently,
✓Z ◆1/n
1/n 2
M 3!n r det r u dx .
Br
1/n
This proves the claim with c = 3!n . ⇤
Remark 21.4. The previous theorem implies the ABP estimate, provided we show that
• u 2 C 1,1 (Br ), as a consequence of u 2 Sup(f );
Since u is convex, the corresponding quantity ✓ is null. Recall also that we have already
proved that ✓, ✓ 2 L1 implies u 2 C 1,1 in Corollary 20.29.
Theorem 21.6 (Regularity at contact points). Consider v 2 Sup(f ) in B , ' convex in
B with 0 ' v and v(0) = '(0) = 0. Then '(x) C supB f + |x|2 in B⌫ , where ⌫
and C are universal constants.
156
We can get a naive interpretation of this lemma (or, better, of its infinitesimal version
as # 0) by this formal argument: v ' having a local minimum at 0 implies, by the
assumption v 2 Sup(f ) M (r2 '(0)) f (0). Formally, M (r2 '(0)) M (r2 v(0))
f (0).
Now it is possible to see how these tools allow to prove the ABP estimate.
Proof. [of Theorem 21.1] Pick a point x0 2 Br \ {u = u } and let L be a supporting
hyperplane for u at x0 , so that u L and u (x0 ) = L (x0 ) . Recalling Theorem 21.6,
define ' := u L, v := u L (and notice that v is a supersolution because v 2
Sup(f Br )). Now, '(x0 ) = v(x0 ) implies, by means of Theorem 21.6,
By Theorem 21.5 we get u 2 C 1,1 and det r2 u = 0 a.e. in the non-contact region.
Finally, in order to get the desired estimate, we have to show that a.e. in the contact
region one has det r2 u c(f + )n . But this comes at once by passing to the limit as ! 0
in (21.2) at any di↵erentiability point x0 of u . In fact, all the eigenvalues of r2 u (x0 )
do not exceed Cf + (x0 ) and the conclusion follows. ⇤
Now we prove Theorem 21.6.
Proof. Let r 2 (0, /4) and call c := supB r ' /r2 . Let then x̄ 2 @Br be a maximum
point of ' on B r (by convexity the maximum is attained at the boundary). By means of
a rotation, we can write x = (x0 , xn ), x0 2 Rn 1 , xn 2 R, and assume x̄ = (0, r). Consider
the intersection A of the closed strip defined by the hyperplanes xn = r and xn = r
with the ball B /2 . We clearly have that @A = A1 [ A2 [ A3 , where A1 = B /2 \ {xn = r},
A2 = B /2 \ {xn = r} and A3 = @B /2 \ {|xn | < r}.
We claim that ' '(x̄) on A1 . To this aim, we first prove that '(y) '(x̄)+o(|y x̄|)
for y ! x̄, y 2 H := {xn = r}. In fact, this comes from '(ry/|y|) '(x̄) and observing
that '(y) '(ry/|y|) = o(|y x̄|), because ' is Lipschitz continuous. On the other hand,
we have that ⇠ 2 @'|H (x̄) implies '(y) '(x̄) + h⇠, y x̄i for all y 2 H. Hence, by
comparison, it must be ⇠ = 0 and so '(y) '(x̄) on A1 (this can be seen as a nonsmooth
version of the Lagrange multipliers theorem).
As a second step, set
c c
p(x) := (xn + r)2 4 2
r2 |x0 |2
8
and notice that the following properties hold:
157
(a) on A1 , p(x) c/(2r2 ) = '(x̄)/2 '(x)/2;
(b) on A2 , p(x) 0 '(x) (and in particular p(x) v(x));
(c) on A3 , 2 /4 = |x0 |2 + x2n |x0 |2 + r2 |x0 |2 + 2 /16, which implies |x0 |2 (3/16) 2 .
By means of the last estimate we get p(x) (c/2)r2 (3/4)cr2 0 '.
Combining (a), (b), (c) above we get p v on @A. Since p(0) = cr2 /8 > 0 = '(0)
we can rigidly move down this paraboloid until we get a limit paraboloid p0 = p ↵ (for
some translation parameter ↵ > 0) lying below the graph of v and touching it at some
point, say y. Since p v on @A, the point y is internal to A.
By the supersolution property M (r2 p) f (y) supB f we get (since we have an
explicit expression for p)
c r2
8(n 1)⇤c 2 sup f.
4 B
2 2
pwe can fix r such that 8(n 1)⇤cr / 8 c/8 (it is done by taking r so that
But now
8r /((n 1)⇤)): we have therefore c supB f . The statement then follows
p
with C = 8/ and ⌫ := 18 /((n 1)⇤)). ⇤
It remains to prove Theorem 21.5.
Proof. Recall first that we are assuming the existence of a uniform estimate
✓(x, B" (x), u) M 8x 2 B r \ {u = u }.
Thanks to Proposition 20.27, we are able to obtain C 1,1 regularity of u as soon we are
able to propagate this estimate also to non-contact points.
Consider now any point x0 2 B r \ {u > u } and call L a supporting hyperplane for
u at x0 . Notice that x0 2 { u = L} ⇢ {u = u }. We claim that:
(a) There exist n + 1 points x1 , . . . , xn+1 such that x0 2 S := co(x1 , . . . , xn+1 ) (here
and in the sequel co stands for convex hull) and, moreover, all such points belong
to B r \ { u = L} with at most one exception lying on @B2r . In addition u ⌘ L
on S;
P
(b) x0 = n+1 i=1 ti xi with at least one index i verifying both xi 2 B r \ { u = L} and
ti 1/(3n).
To show the utility of this claim, just consider how these two facts imply the thesis: on
the one hand, if r u is di↵erentiable at x0 , we get det r2 u (x0 ) = 0 because u = L on
S and dim(S) 1. On the other hand we may assume, without loss of generality that
x1 2 {u = u } \ Br and t1 (1/3n) so that, since
✓ ◆
h
x0 + h = t 1 x1 + + t2 x2 + · · · + tn+1 xn+1 ,
t1
158
one has
and this estimate is clearly uniform since we only require |h/t1 | ", which is implied by
|h| "/(3n).
Hence, the problem is reduced to prove the two claims above. This is primarily based
on a standard result in convex analysis (first proved by Carathéodory for closed sets),
which is recalled here for completeness.
Hence, we would have (L + L0 )(x0 ) > L(x0 ) and, at the same time,
L + L0 u on B 2r ,
159
Let us now prove that u (x) = L(x) on S := co(x1 , . . . , xn+1 ). The implication is
trivial, the converse one is clear for each x = xi , since L u u , and it is obtained
by means of the convexity of u at all points in S.
Now we prove part (b) of the claim. If all points xj verify |xj | r, then max ti
1 1
n+1
> 3n . Otherwise, if one point, say xn+1 , satisfies |xn+1 | = 2r, then ti < 1/(3n) for all
i = 1, . . . , n implies tn+1 > 2/3 and therefore
n
X 4 n
r |x0 | 2tn+1 r ti |xi | > r r=r.
i=1
3 3n
Let us show how (21.5) leads to the Hölder regularity result for viscosity solutions of
the fully nonlinear elliptic PDE
F (r2 u(x)) + f (x) = 0 . (21.6)
Step 1. As usual, we need to control the oscillation (now on cubes), defined by
!r := Mr mr
10
Notice that Sup(f ) ⇢ Sup(|f |) and Sub(f ) ⇢ Sub( |f |).
160
with Mr := supQr u and mr := inf Qr u.
With the same notation of Theorem 21.8, there exists a universal constant µ 2 (0, 1) such
that
!1/2 µ !1 + 2kf kLn (Q1 ) . (21.7)
Indeed, we apply the Harnack inequality (21.5)
• to the function u m1 , so that
M1/2 m1 CH m1/2 m1 + kf kLn (Q1 ) ; (21.8)
161
and with C depending only on µ and kf kLn (Q1 ) , thus we have Hölder regularity.
In order to prove the Harnack inequality, we will pass through the following reformu-
lation of Theorem 21.8.
Theorem 21.9. There exist universal positive constants "0 , C such that if u : Q4pn !
[0, 1) belongs to Sol(f ) \ C(Q4pn ) on Q4pn , then
provided
kf kLn (Q4pn ) "0 .
Remark 21.10. Theorem 21.8 and Theorem 21.9 are easily seen to be equivalent: since
we will prove the second one, it is more important for us to check that Theorem 21.8
follows from Theorem 21.9.
For some positive > 0 (needed to avoid a potential division by 0) consider the function
u
v := .
+ inf Q1/4 u + kf kLn (Q4pn ) /"0
Denoting by fv the source term associated with v, the homogenity of Pucci’s operators
gives kfv kLn (Q4pn ) "0 . Since inf Q1/4 v 1 we have supQ1/4 v C, hence
✓ ◆
sup u C inf u + + kf kLn (Q4pn ) /"0 .
Q1/4 Q1/4
We let ! 0 and we obtain Harnack inequality with the cubes Q1/4 , Q4pn ; by the same
scaling argument we already used, this means
✓ ◆
sup u C inf u + rkf kLn (Q16rpn (x0 ) ) . (21.11)
Qr (x0 ) Qr (x0 )
Now, we pass to the cubes Q1/2 , Q1 with a simple covering argument: there exists an
integer N = N (n) such that for all x 2 Q1/2 , y 2 Q1 we can find points xi , 1 i N ,
with xi = x, xN = y and xi+1 2 Qr (xi ) for 1 i < N , with r = r(n) so small that all
cubes Q16rpn (xi ) are contained in Q1 . Applying repeatedly (21.11) we get (21.5) with
CH ⇠ C N .
We describe the strategy of the proof of Theorem 21.9, even if the full proof will be
completed at the end of this section.
We will study the map
t 7! L n ({u > t} \ Q1 )
in order to prove:
162
• a decay estimate of the form L n ({u > t} \ Q1 ) dt " , thanks to the fact that
u 2 Sup(|f |) (see Lemma 21.13),
• the full thesis of Theorem 21.9 using the fact that u 2 Sol(f ) ⇢ Sub( |f |), too.
The first goal will be achieved using the Alexandrov-Bakelman-Pucci inequality of the
previous section. The structure of the proof remembers that of De Giorgi’s regularity
theorem, as we said, and we will complete it through the following lemmas and remarks.
The first lemma is a particular case of Calderón-Zygmund decomposition.
Lemma 21.11 (Dyadic Lemma). Consider Borel sets A ⇢ B ⇢ Q1 with L n (A) < 1.
If the implication
L n (A \ Q) > L n (Q) =) Q̃ ⇢ B , (21.12)
holds for any dyadic cube Q ⇢ Q1 , with Q̃ being the predecessor of Q, then
L n (A) L n (B) .
Proof. We apply the construction of Calderón-Zygmund (seen in the proof of The-
orem 14.1) to f = A with ↵ = : there exists a countable family of cubes {Qi }i2I ,
pairwise disjoint, such that
[
A L n -a.e. on Q1 \ Qi (21.13)
i2I
and L n (A \ Qi ) > L n (QSi ) for all i 2 I. Since < 1 and A is a characteristic function,
(21.13) means that A ⇢ i2I Qi up to Lebesgue negligible sets. Moreover, if Q̃i are the
predecessors of Qi , from (21.12) we get Q̃i ⇢ B for all i and
This is due to the fact that a cube Q, in the Calderón-Zygmund construction, is divided
in subcubes as long as L n (A \ Q) L n (Q). Thus (note that we sum on Q̃i rather than
on i, because di↵erent cubes might have the same predecessor)
X X
L n (A) L n (A \ Q̃i ) L n (Q̃i ) L n (B) .
Q̃i Q̃i
⇤
It is bothering, but necessary to go on with the proof, to deal at the same time with
balls and cubes: balls emerge from the radial construction in the next lemma and cubes
are needed in Calderón-Zygmund Theorem.
Lemma 21.12 (Truncation Lemma). There exists a universal function ' 2 C 1 (Rn ) such
that
163
(i) ' 0 on Rn \ B2pn ;
Since ' is an increasing function of |x|, we can find M1 = M1 (↵) > 0 and M2 =
M2 (↵) > 0 such that
After choosing a smooth extension for ' on B1/2 , still less than 2, we conclude check-
ing that there exists an exponent ↵ that is suitable to verify the third property of the
statement, that needs to be checked only on . We compute
↵ x⌦x
r2 |x| ↵
= ↵+2
I + ↵(↵ + 2) ↵+4 ,
|x| |x|
thus the eigenvalues of r2 ' when |x| 1/2 are M2 ↵|x| (↵+2) with multiplicity n 1
and M2 ↵(↵ + 1)|x| (↵+2) with multiplicity 1 (this is the eigenvalue due to the radial
direction). Hence, when |x| 1/4 we have
M2
M+ (r2 ') = (⇤(n 1)↵ ↵(↵ + 1))
|x|↵+2
Lemma 21.13 (Decay Lemma). There exist universal constants "0 > 0, M > 1 and
µ 2 (0, 1) such that if u 2 Sup(|f |), u 0 on Q4pn , inf Q3 u 1 and kf kLn (Q4pn ) "0 ,
then for every integer k 1
164
Proof. We prove the first step, that is
L n ({u > M } \ Q1 ) (1 µ) , (21.16)
with M := max ' , ' given by Lemma 21.12, and µ and "0 are respectively given by
n 1
µ := (2CABP C' ) , "0 = , (21.17)
2CABP
where CABP is the universal constant of the Alexandrov-Bakelman-Pucci estimate of The-
orem 21.1. Since u is nonnegative, in order to obtain a meaningful result from the ABP
estimate, we apply the estimate in the ball B2pn for the function w, defined as the function
u additively perturbed with the truncation function '. If w := u + ', then
(i)
w 0 on @B2pn (21.18)
because u 0 on Q4pn B2pn and ' 0 on Rn \ B2pn ;
(ii)
inf w inf w 1 (21.19)
B2pn Q3
because Q3 ⇢ B2pn and ' 2 on B2pn , and, at the same time, we are assuming
that inf Q3 u 1;
(iii) directly from the definition of Sup(|f |) we get M (r2 u) + |f | 0, moreover
M+ (r2 ') C' Q1 . Since in general M (A + B) M (A) + M+ (B) (see Re-
mark 20.35), then
M (r2 w) + (|f | + C' Q1 ) ( M (r2 u) + |f |)+( M+ (r2 ') + C' Q1 ) 0.
(21.20)
The inequality (21.20) means that w 2 Sup(|f | + C' Q1 ).
Thanks to the ABP estimate (which we can apply to w thanks to (21.18) and (21.20))
we get
✓Z ⇣ ⌘n ◆1/n
max w (x) CABP |f (y)| + C' Q1 (y) dy . (21.21)
x2B 2pn {w= w}
165
where we pass from line (21.22) to line (21.23) by Minkowski inequality and from line
(21.23) to line (21.24) because, if w(x) 0, then u(x) '(x) and then u(x) M .
Using our choice of "0 we obtain from (21.24) the lower bound
1
L n (Q1 \ {u M })1/n . (21.25)
2CABP C'
Concerning (21.12), suppose by contradiction that for some dyadic cube Q ⇢ Q1 we have
that
L n (A \ Q) > L n (Q) (21.26)
but Q̃ 6⇢ B, Q̃ being the predecessor of Q, as usual: there exists z 2 Q̃ such that
u(z) M k 1 . Let us rescale and translate the problem, putting ũ(y) := u(x)M (k 1)
with x = x0 + 2 i y if Q has edge length 2 i and centre x0 (so that, in this transformation
Q becomes the unit cube and Q̃ is contained in Q3 ). Because of the rescaling technique,
we need to adapt f , that is define a new datum
f (x)
f˜(y) := .
22i M k 1
The intention of this definition of f˜ is to ensure that ũ 2 Sup(|f˜|), in fact
1
M (r2 ũ) + |f˜| = M (r2 u) + |f | 0.
22i M k 1
Since the point corresponding to z belongs to Q3 , we ge
u(z)
inf ũ(y) 1.
y2Q3 Mk 1
If kf˜kLn (Q4pn ) "0 , then, applying what we already saw in (21.25) to ũ instead of u,
µ L n ({ũ M } \ Q1 ) = 2ni L n {u M k } \ Q ,
166
which contradicts (21.26).
In order to complete our proof, we show that e↵ectively kf˜kLn (Q4pn ) "0 . In general,
let us remark that the rescaling technique does not cause any problem at the level of the
source term f . Indeed
1
kf˜kLn (Q4pn ) = kf kLn (Q4pn/2i (x0 )) kf kLn (Q4pn ) "0 .
M k 1 2i
⇤
Corollary 21.14. There exist universal constants " > 0 and d 0 such that if u 2
Sup(|f |), u 0 on Q4pn , inf Q3 u 1 and kf kLn (Q4pn ) "0 , then
Lemma 21.15. Suppose that u 2 Sub( |f |) is nonnegative on Q4pn and kf kLn (Q4pn )
"0 , with "0 given by the decay Lemma 21.13. Assume that (21.27) holds. Then there exist
universal constants M0 > 1 and > 0 such that if
then
9 x1 2 Qlj (x0 ) such that u(x1 ) M0 ⌫ j ,
"/n "j/n
where ⌫ := 2M0 /(2M0 1) > 1 and lj := M0 ⌫ .
Proof. First of all, we fix a large universal constant > 0 such that
1 n
p > d2" (21.28)
2 4 n
167
and then we choose another universal constant M0 so large that
1
dM0 " < (21.29)
2
and p
"/n
M0 <2 n. (21.30)
We first estimate the superlevels
where we used condition (21.28) on and the definition of lj , as given in the statement
of the lemma.
By contradiction, assume that for some j 1 we have
Under this assumption, we claim that the superlevel can be estimated as follows:
1
L n {u < ⌫ j M0 /2} \ Qlj /(4pn) (x0 ) < L n Qlj /(4pn) . (21.33)
2
Obviously the validity of (21.31) and (21.33) is the contradiction that will conclude the
proof, so we need only to show (21.33).
Define the auxiliary function
⌫M0 u(x)⌫ (j 1)
u(x)
v(y) := = 2 M0 ,
(⌫ 1)M0 ⌫j
l
where x = x0 + 4pj n y and the second equality is a consequence of the relation M0 =
⌫/[2(⌫ 1)]. Since y 2 Q4pn () x 2 Qlj (x0 ), by (21.32) the function v is defined
and positive on Q4pn . In addition, using the first equality in the definition of v, we
immediately see that u(x0 ) M0 ⌫ j 1 implies inf Q4pn v 1.
Using the second equality we see that (modulo the change of variables)
168
Moreover,p if we compute the datum fv which corresponds to v, since the rescaling radius
is lj /(4 n), we get
2lj2
fv (y) = j f (x)
⌫
so that
2lj
kfv kLn (Q4pn ) = p j kf kLn (Qlj (x0 )) "0 (21.34)
4 n⌫
because
"/n
2l M0
pj j = p ⌫ "j/n j < 1
4 n⌫ 2 n
thanks to (21.30). The estimate in (21.34) allows us to use Corollary 21.14 for v, that is
L n ({v > M0 } \ Q1 ) dM0 " ,
and we can use this, together with (21.29), to obtain that (21.33) holds:
1
L n {u < ⌫ j M0 /2} \ Qlj /(4pn) (x0 ) dM0 " L n Qlj /(4pn) < L n Qlj /(4pn) .
2
⇤
We can now complete the proof of Theorem 21.9, using Lemma 21.15. Notice that
in Theorem 21.9 we made all assumptions needed to apply Lemma 21.15, taking also
Corollary 21.14 into account, which ensures the validity of (21.27).
Roughly speaking, if we assume, by (a sort of) contradiction, that u is not bounded from
above by M ⌫ k0 on Q1/4 for k0 sufficiently large, then, thanks to Lemma 21.15, we should
be able to find recursively a sequence (xj ) with the property that
u(xj ) M0 ⌫ j and xj+1 2 Qlj (xj );
P
since j lj < 1, the sequence (xj ) admits a converging subsequence, and in the limit
point we find a contradiction. However, in order to iterate Lemma 21.15 we have to
confine the sequence in the cube Q1/2 (for this purpose it is convenient to use the distance
induced by the L1 norm in Rn , whose balls are cubes). P
To achieve this, we fix a universal positive integer j0 such that j j0 lj < 1/4 and we
assume, by contradiction, that there exists a point x0 2 Q1/4 with u(x0 ) M0 ⌫ j0 1 . This
time, the sequence (xk ) we generate iterating Lemma 21.15 is contained in Q1/2 and
u(xk ) M0 ⌫ j0 +k 1
. (21.35)
When k ! 1 in (21.35) we obtain the contradiction. This way, we obtained also an
“explicit” expression of the universal constant in (21.10), in fact we proved that
sup u(x) M0 ⌫ j0 1
.
x2Q1/4
References
[1] R.Adams: Sobolev spaces. Academic Press, 1975.
[3] L.Ambrosio, N.Fusco, D.Pallara: Functions of bounded variation and free dis-
continuity problems. Oxford University Press, 2000.
[6] E.De Giorgi: Complementi alla teoria della misura (n 1)-dimensionale in uno
spazio n-dimensionale. Seminario di Matematica della Scuola Normale Superiore di
Pisa, (1960-61), Editrice Tecnico Scientifica, Pisa.
[7] E.De Giorgi: Frontiere orientate di misura minima. Seminario di Matematica della
Scuola Normale Superiore di Pisa, (1960-61), Editrice Tecnico Scientifica, Pisa.
[9] E.De Giorgi: Sulla di↵erenziabilità e l’analicità degli estremali degli integrali mul-
tipli regolari. Mem. Acc. Sc. Torino, 3 (1957), 25–43.
[10] L.C.Evans: Quasiconvexity and partial regularity in the calculus of variations. Arch.
Rational Mech. Anal. 95, 3 (1986), 227–252.
[13] E.Gagliardo: Caratterizzazione delle tracce sulla frontiera relative ad alcune class
di funzioni in piú variabili. Rend. Sem. Mat. Univ. Padova, 27 (1957), 284–305.
170
[16] M.Giaquinta, E.Giusti: On the regularity of the minima of variational integrals.
Acta Math. 148, (1982), 31–46.
[17] M.Giaquinta, E.Giusti: Quasiminima. Ann. Inst. H. Poincaré Anal. Non Linéaire
1, 2 (1984), 79–107.
[18] M.Giaquinta, E.Giusti: The singular set of the minima of certain quadratic func-
tionals. Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4) 11, 1 (1984), 45–55.
[19] M.Giaquinta: Multiple integrals in the Calculus of Variations and Nonlinear elliptic
systems. Princeton University Press, 1983.
[20] E.Hopf: Über den funktionalen, insbesondere den analytischen Charakter der
Lösungen elliptischer Di↵erentialgleichungen zweiter Ordnung. Math. Zeitschrift,
Band 34 (1932), 194–233.
[23] The singular set of lipschitzian minima of multiple integrals. Arch. Ration. Mech.
Anal. 184 (2007), 341–369.
[24] N.G.Meyers, J.Serrin: H = W . Proc. Nat. Acad. Sci. U.S.A., 51 (1964), 1055–
1056.
[25] S.Müller, V.Sverak: Convex integration for Lipschitz mappings and counterex-
amples to regularity. Ann. of Math., 157 (2003), 715–742.
[27] N.Trudinger: On embedding into Orlicz spaces and some applications. J. Math.
Mech., 17 (1967), 473–483.
[28] K.Yosida: Functional Analysis. Mathematical surveys and monographs, 62, Amer-
ican Mathematical Society, 1998.
171