Table Of ContentA Liouville Theorem for the Fractional
4
Laplacian
1
0
2
n Ran Zhuo Wenxiong Chen Xuewei Cui Zixia Yuan
a
J
9
2
Abstract
]
P Inthispaper,weconsiderthefollowingfractionalLaplaceequation
A
. (−∆)α/2u(x) = 0, in Rn,
h (1)
t ( u(x) ≥0, in Rn,
a
m
where n ≥ 2 and α is any real number between 0 and 2. We prove
[
that the only solution for (1) is constant. Or equivalently,
1
Every α-harmonic function bounded either above or below in all of
v
2 Rn must be constant.
0 Thisextends theclassical Liouville TheoremfromLaplacian tothe
4
fractional Laplacian.
7
. As an immediate application, we use it to obtain an equivalence
1
between a semi-linear pseudo-differential
0
4
1 (−∆)α/2u= up(x), x ∈ Rn (2)
:
v
i and the corresponding integral equation
X
r 1
a u(x) = dy, x ∈ Rn.
ZRn |x−y|n−α
Combining this with the existing results on the integral equation, one
can obtained much more general results on the qualitative properties
of the solutions for (2).
1 Introduction
The well-known Liouville’s Theorem states that
1
Any harmonic function bounded either above or below in all of Rn is
constant.
One of its important applications is the proof of the Fundamental Theo-
rem of Algebra. It is also a key ingredient in deriving a priori estimates for
solutions in PDE analysis.
The main purpose of this article is to extend this classical theorem to the
fractional Laplacian.
Essentially different from the Laplacian, the fractional Laplacian in Rn is
a nonlocal operator, taking the form
u(x)−u(z)
(−∆)α/2u(x) = C PV dz (3)
n,α
ZRn |x−z|n+α
where α is any real number between 0 and 2 and PV stands for the Cauchy
principal value. This operator is well defined in S, the Schwartz space of
rapidly decreasing C∞ functions in Rn. In this space, it can also be defined
equivalently in terms of the Fourier transform
(−∆)α/2u(ξ) = |ξ|αuˆ(ξ)
d
where uˆ is the Fourier transform of u. One can extend this operator to a
wider space of distributions as the following.
Let
|u(x)|
L = {u : Rn → R | < ∞}.
α ZRn (1+|x|n+α)
For u ∈ L , we define (−∆)α/2u as a distribution:
α
< (−∆)α/2u(x),φ >=< u,(−∆)α/2φ >, ∀φ ∈ C∞(Rn).
0
This defines the operator in a weak sense. In this paper, we consider the
class of functions where the fractional Laplacian is defined in a little bit
more stronger sense, that is
u ∈ L , such that the right hand side of (3) is well defined for every
α
x ∈ Rn.
One can verify that, all the above definitions coincides when u is in S.
We say that u is an α-harmonic function if u ∈ L , such that the right
α
hand side of (3) is well defined for every x ∈ Rn and equals zero. In this
sense, we have
2
Theorem 1 Every α-harmonic function bounded either above or below in all
of Rn for n ≥ 2 must be constant.
This is the main result of the paper. To prove it, we study
(−∆)α/2u(x) = 0, in Rn,
(4)
u(x) ≥ 0, in Rn.
(
We say that u ≥ 0 is a strong solution of (4), if u ∈ L , such that the right
α
hand side of (3) is well defined for every x ∈ Rn and equals zero.
Apparently, Theorem 1 is equivalent to the following
Theorem 2 Assume that n ≥ 2. Let u be a strong solution of (4), then
u ≡ C.
As an immediate application of the Liouville theorem for α harmonic
functions, we prove an equivalence between a psedodifferential equation and
an integral equation.
Theorem 3 Assume that n ≥ 2 and u ∈ L is a nonnegative strong solution
α
of
(−∆)α/2u(x) = up(x), x ∈ Rn, (5)
then u also satisfies
c
u(x) = n up(y)dy,
ZRn |x−y|n−α
and vice versa.
Remark 1 i) Actually, the right hand side of equation (5) can be a much
more general function f(x,u), such that, for any constant c > 0,
1
f(y,c)dy = ∞. (6)
ZRn |x−y|n−α
ii) The idea of proof can be extended to establish the equivalence between
a general system of m equations in Rn
(−∆)α/2u (x) = f (x,u (x),···u (x)), i = 1,···,m,
i i 1 m
u ≥ 0, i = 1,···,m,
( i
and the corresponding integral system
u (x) = cn f (y,u (y),···,u (y)), i = 1,···,m,
i Rn |x−y|n−α i 1 m
( ui(x) ≥ R0, i = 1,···,m.
3
Combining Theorem 3 with the qualitative properties established for the
integral equations in [CLO] and [CLO1], one obtain immediately that
Theorem 4 Assume that n ≥ 2 and u is a nonnegative strong solution of
(5) for 0 < α < 2. Then
i) In the critical case when p = n+α, it must assume the form
n−α
t
u(x) = c( )(n−α)/2
t2 +|x−x |2
o
for some t > 0, x ∈ Rn.
o
ii) In the subcritical case when 1 < p < n+α, we must have u ≡ 0.
n−α
Remark 2 i) In [CLO] and [CLO1], in order the results in Theorem 4 to
hold, one requires u to be in Hα/2(Rn). Here we only requires u ∈ L , a
α
much weaker restriction.
ii) In [BCPS], by using the extension method to obtain the same results
as in Theorem 4, the authors require that 1 ≤ α < 2 and u be bounded.
Obviously, our condition here is much weaker.
In Section 2, we prove the Liouville Theorem 2 and hence Theorem 1. In
Section 3, we establish the equivalence and hence prove Theorem 3 and 4.
2 The proof of the Liouville Theorem
In this section, we prove Theorem 2.
Proof. First, we define
u(x), |x| ≤ k,
u (x) = (7)
k P (y,x)u(y)dy, |x| > k,
( Bk k
R
where P (y,x) is a Poisson kernel in the exterior of the ball B with radius
k k
k and centered at the origin:
n n πα (|x|2 −k2)α/2 1
Pk(y,x) = Γ( )π−2−1sin( ) , |y| < k, |x| > k.
2 2 (k2 −|y|2)α/2 |x−y|n
(8)
4
Obviously, for each x,
lim u (x) = u(x).
k
k→∞
One can also verify (see [L]) that
(−∆)α/2u (x) = 0, for |x| > k, and u (x) ≤ u(x), ∀x ∈ Rn. (9)
k k
Moreover, by Taylor expansion, it is easy to derive that
c 1
1
u (x) = +O( ). (10)
k
|x|n−α |x|n−α+1
In order to prove that u is constant, it suffice to show that for any ψ ∈
C∞(Rn), satisfying the condition
0
ψ(x)dx = 0, (11)
ZRn
we have
u(x)ψ(x)dx = 0.
ZRn
Actually, we only need to prove
lim u (x)ψ(x)dx = 0.
k
k→∞ZRn
We divided the proof into two steps.
Step 1. Let
ψ(y)
ϕ(x) = .
ZRn |x−y|n−α
Combining Taylor expansion with (11), we deduce that
1
ϕ(x) = O( ), as|x| → ∞. (12)
|x|n−α+1
It follows that ϕ(x) ∈ L2(Rn) for n ≥ 2 , and
(−∆)α/2ϕ(x) = ψ(x), x ∈ Rn. (13)
In this step, we will show that
u (x)(−∆)α/2ϕ(x)dx = (−∆)α/2u (x)ϕ(x)dx. (14)
k k
ZRn ZRn
5
Let
c
1
v (x) = u (x)− ,
k k
|x|n−α
then v (x) ∈ L2(Rn) due to (10).
k
Applying the Parseval’s formula to one part of the left hand side of (14),
we derive that
c
u (x)(−∆)α/2ϕ(x)dx = (v (x)+ 1 )(−∆)α/2ϕ(x)dx
ZRn k ZRn k |x|n−α
c
= v (x)(−∆)α/2ϕ(x)dx+ 1 (−∆)α/2ϕ(x)dx
ZRn k ZRn |x|n−α
= v (ξ)|ξ|αϕ(ξ)+cϕ(0). (15)
k
ZRn
Herewehaveusedaresultin[L]thact,inthesbenseofdistributions, theFourier
transform of c1 is a constant multiple of |ξ|−α, and by the definition of
|x|n−α
Fourier transform on distributions (see [L]), we have
c
1 (−∆)α/2ϕ(x)dx = |ξ|−α|ξ|αϕ(ξ)dξ = cϕ(0).
ZRn |x|n−α ZRn
b
Also note that we are not able to apply the Parseval’s formula directly to
u (x)(−∆)α/2ϕ(x)dx because u may not be in L2(Rn).
Rn k k
For the right hand side of (14), we have
R
c
(−∆)α/2u (x) ϕ(x)dx = (−∆)α/2(v (x)+ 1 )ϕ(x)dx
ZRn k ZRn k |x|n−α
c
= (−∆)α/2v (x) ϕ(x)dx+ (−∆)α/2( 1 ) ϕ(x)dx
k
ZRn ZRn |x|n−α
= |ξ|αv (ξ)ϕ(ξ)dξ +cϕ(0). (16)
k
ZRn
Here we have used a well-known fact thcat b1 is a constant multiple of the
|x|n−α
fundamental solution of (−∆)α/2.
Now from (15) and (16), we arrive at (14).
Step 2. We prove
(−∆)α/2u (x) ϕ(x)dx → 0, ask → ∞. (17)
k
ZRn
6
By elementary calculation, we separate the integral in(17) into two parts,
u (x)−u (y)
(−∆)α/2u (x)ϕ(x)dx = c k k dyϕ(x)dx
ZRn k ZRnZRn |x−y|n+α
u (x)−u (y)
k k
= c dyϕ(x)dx
ZBr(0)ZRn |x−y|n+α
u (x)−u (y)
k k
+ c dyϕ(x)dx
ZRn\Br(0)ZRn |x−y|n+α
= I +I ,
1 2
where r < k.
First, we consider I .
1
From the first equation of (4), we have
u(x)−u(y)
0 = c dy
ZRn |x−y|n+α
u(x)−u(y) u(x)−u(y)
= c dy +c dy.
ZBk(0) |x−y|n+α ZRn\Bk(0) |x−y|n+α
It follows that
u (x)−u (y)
k k
I = c dyϕ(x)dx
1 ZBr(0)ZRn |x−y|n+α
u(x)−u (y)
k
= c dyϕ(x)dx
|x−y|n+α
ZBr(0)ZBk(0)
u(x)−u (y)
k
+ c dyϕ(x)dx
ZBr(0)ZRn\Bk(0) |x−y|n+α
u(y)−u (y)
k
= c dyϕ(x)dx. (18)
ZBr(0)ZRn\Bk(0) |x−y|n+α
Applying (18), we obtain, for each fixed r,
|u(y)|
|I | ≤ c dy|ϕ(x)|dx
1 ZBr(0)ZRn\Bk(0) |x−y|n+α
|u(y)|
≤ c dy|ϕ(x)|dx → 0, ask → ∞. (19)
ZBr(0)ZRn\Bk(0) (1+|y|)n+α
Here we have used the facts that u ∈ L and 0 ≤ u ≤ u.
α k
7
Next we will show that
I = c (−∆)α/2u ϕ(x)dx → 0, asr → ∞, (20)
2 k
ZRn\Br(0)
uniformly in k.
Let
f (x) = (−∆)α/2u (x).
k k
We first show that
f (y)
k
dy = u (x). (21)
ZRn |x−y|n−α k
To this end, denote
f (y)
k
g (x) = dy.
k
ZRn |x−y|n−α
For any ψ(x) ∈ C∞(Rn), we show that
0
g (x)ψ(x)dx = u (x)ψ(x)dx. (22)
k k
ZRn ZRn
Actually,
(−∆)α/2u (y)
k
g (x)ψ(x)dx = dyψ(x)dx
ZRn k ZRn ZRn |x−y|n−α
(−∆)α/2(v (y)+c /|y|n−α)
k 1
= dyψ(x)dx
ZRn ZRn |x−y|n−α
(−∆)α/2v (y) c
k 1
= dyψ(x)dx+ ψ(x)dx
ZRn ZRn |x−y|n−α ZRn |x|n−α
= I +I .
1 2
Since both (−∆)α/2v (x) and ψ(x) are compactly supported, one can
k
exchange the order of integration to derive
ψ(x)
I = (−∆)α/2v (y) dxdy
1 ZRn k ZRn |x−y|n−α
= |ξ|αv (ξ)ψ(ξ)|ξ|−αdξ
k
ZRn
= v (ξc)ψ(ξ)bdξ
k
ZRn
= vc(x)ψb(x)dx.
k
ZRn
8
Herewehaveusedaresult in[LL]thattheFouriertransformof ψ(x) dx
Rn |x−y|n−α
is a constant multiple of ψ(ξ)|ξ|−α. It follows that R
b
g (x)ψ(x)dx = I +I
k 1 2
ZRn
c
1
= v (x)+ ψ(x)dx
ZRn k |x|n−α!
= u (x)ψ(x)dx.
k
ZRn
This proves (22) and hence (21), and from which, we arrive immediately that
f (x)
k
dx ≤ u (0) ≤ u(0). (23)
ZRn\Br(0) |x|n−α k
Combining (12) with (23), we derive that
u (x)−u (y)
k k
I = c dy ϕ(x)dx
2 ZRn\Br(0)ZRn |x−y|n+α
= c f (x)ϕ(x)dx
k
ZRn\Br(0)
f (x)
k
≤ c dx
ZRn\Br(0) |x|n−α+1
c f (x)
k
≤ dx
r ZRn\Br(0) |x|n−α
c
≤ u(0) → 0, asr → ∞, uniformly in k. (24)
r
(19) and (24) imply that (17) holds. Hence we have
lim u (x)ψ(x)dx = lim u (x)(−∆)α/2ϕ(x)dx
k k
k→∞ZRn k→∞ZRn
= lim (−∆)α/2u (x)ϕ(x)dx = 0. (25)
k
k→∞ZRn
That is
u(x)ψ(x)dx = 0.
ZRn
Therefore we come to the conclusion that u ≡ C.
This complete the proof of Theorem 2.
9
The Proof of Theorem 1.
To see that Theorem 2 implies Theorem 1, let v be any α-harmonic func-
tion that is bounded from above by a constant M in Rn. Take u(x) =
M −v(x), then
(−∆)α/2u(x) = 0, x ∈ Rn,
u(x) ≥ 0, x ∈ Rn.
(
By Theorem 2, u must be constant, hence so does v. Similarly, if v is any α-
harmonic function that is bounded from below by a constant M in Rn, then
we let u(x) = v(x)−M to derive that v must be constant. This completes
the proof of Theorem 1.
3 Applications
The Proof of Theorem 3.
Assume u ∈ L is a nonnegative locally bounded strong solution of
α
(−∆)α/2u(x) = up(x), x ∈ Rn, (26)
Let
v (x) = G (x,y)up(y)dy, (27)
R R
ZBR
where G (x,y) is Green’s function on the ball B (0):
R R
(−∆)α/2G (x,y) = δ(x−y), x,y ∈ B (0),
R R
G (x,y) = 0, x or y ∈ Bc(0).
( R R
Thanks to Kulczycki [Ku], one can write
s (n−α)
A B (s−tb)
n,α n,α t 2
G (x,y) = 1− db , x,y ∈ B (0),
R s(n−2α) (s+t)(n−2α) Z0 bα2(1+b) R
wheres = |x−y|2,t = (1−|x|2)(1−|y|2). A andB areconstantsdepending
R2 R2 R2 n,α n,α
on n and α.
It is easy to verify that
(−∆)α/2v (x) = up(x), in B (0),
R R (28)
v = 0, in Bc(0).
( R R
10