G B Folland Lectures on Partial Differential Equations

background image

Lectures on

Partial Di

fferential Equations

By

G.B. Folland

Tata Institute of Fundamental Research

Bombay

1983

background image

Lectures on

Partial Di

fferential Equations

By

G.B. Folland

Lectures delivered at the

Indian Institute of Science, Bangalore

under the

T.I.F.R. – I.I.Sc. Programme in Applications of

Mathematics

Notes by

K.T. Joseph and S. Thangavelu

Published for the

Tata Institute of Fundamental Research Bombay

Springer-Verlag

Berlin Heidelberg New York

1983

background image

Author

G.B. Folland

University of Washington

Seattle, Washington 98175

U.S.A.

©Tata Institute of Fundamental Research, 1983

ISBN 3-540-12280-X Springer-Verlag, Berlin, Heidelberg. New York
ISBN 0-387-12280-X Springer-Verlag, New York. Heidelberg. Berlin

No part of this book may be reproduced in any
form by print, microfilm or any other means with-
out written permission from the Tata Institute of
Fundamental Research, Colaba, Bombay 400 005

Printed by N.S. Ray at The Book Center Limited,

Sion East, Bombay 400 022 and published by H. Goetze,

Springer-Verlag, Heidelberg, West Germany

Printed in India

background image
background image

Preface

This book consists of the notes for a course I gave at the T.I.F.R. Center
in Bangalore from September 20 to November 20, 1981. The purpose
of the course was to introduce the students in the Programme in Appli-
cation of Mathematics to the applications of Fourier analysis-by which I
mean the study of convolution operators as well as the Fourier transform
itself-to partial differential equations. Faced with the problem of cover-
ing a reasonably broad spectrum of material in such a short time, I had to
be selective in the choice of topics. I could not develop any one subject
in a really thorough manner; rather, my aim was to present the essential
features of some techniques that are well worth knowing and to derive
a few interesting results which are illustrative of these techniques. This
does not mean that I have dealt only with general machinery; indeed,
the emphasis in Chapter 2 is on very concrete calculation with distribu-
tions and Fourier transforms-because the methods of performing such
calculations are also well worth knowing.

If these notes suffer from the defect of incompleteness, they posses

the corresponding virtue of brevity. They may therefore be of value to
the reader who wishes to be introduced to some useful ideas without
having to plough through a systematic treatise. More detailed accounts
of the subjects discussed here can be found in the books of Folland [1],
Stein [2], Taylor [3], and Treves [4].

No specific knowledge of partial differential equations or Fourier

Analysis is presupposed in these notes, although some prior acquittance
with the former is desirable. The main prerequisite is a familiarity with
the subjects usually gathered under the rubic “real analysis”: measure

v

background image

vi

and integration, and the elements of point set topology and functional
analysis. In addition, the reader is expected to be acquainted with the
basic facts about distributions as presented, for example, in Rudin [7].

I wish to express my gratitude to professor K.G. Ramanathan for

inviting me to Bangalore, and to professor S.Raghavan and the staff of
the T.I.F.R. Center for making my visit there a most enjoyable one. I also
wish to thank Mr S. Thangavelu and Mr K.T. Joseph for their painstak-
ing job of writing up the notes.

background image

Contents

Preface

v

1

Preliminaries

1

1

General Theorems About Convolutions

. . . . . . . . .

1

2

The Fourier Transform . . . . . . . . . . . . . . . . . .

5

3

Some Results From the Theory of Distributions . . . . . 12

2

Partial Di

fferential Operators with Constant Coefficients

15

1

Local Solvability and Fundamental Solution . . . . . . . 15

2

Regularity Properties of Differential Operators . . . . . . 20

3

Basic Operators in Mathematical Physics

. . . . . . . . 22

4

Laplace Operator . . . . . . . . . . . . . . . . . . . . . 25

5

The Heat Operator

. . . . . . . . . . . . . . . . . . . . 41

6

The Wave Operator . . . . . . . . . . . . . . . . . . . . 44

3

L

2

Sobolev Spaces

53

1

General Theory of L

2

Sobolev Spaces . . . . . . . . . . 53

2

Hypoelliptic Operators With Constant Coefficients . . . 64

4

Basic Theory of Pseudo Di

fferential Operators

69

1

Representation of Pseudo differential Operators . . . . . 69

2

Distribution Kernels and the Pseudo Local Property . . . 73

3

Asymptotic Expansions of Symbols . . . . . . . . . . . 77

4

Properly Supported Operators

. . . . . . . . . . . . . . 79

5

ψdo

s Defined by Multiple Symbols . . . . . . . . . . . 81

vii

background image

viii

Contents

6

Products and Adjoint of ψDO

S

. . . . . . . . . . . . . 87

7

A Continuity Theorem for ψ Do on Sobolev Spaces . . . 90

8

Elliptic Pseudo Differential Operators . . . . . . . . . . 92

9

Wavefront Sets . . . . . . . . . . . . . . . . . . . . . . 95

10

Some Further Applications... . . . . . . . . . . . . . . . 99

5

L

P

and Lipschitz Estimates

107

background image

Chapter 1

Preliminaries

IN THIS CHAPTER, we will study some basic results about convolu-

1

tion and the Fourier transform.

1 General Theorems About Convolutions

We will begin with a theorem about integral operators.

Theorem 1.1. Let K be a measurable function on R

n

× R

n

such that, for

some c > 0,

R |K(x, y)|dy c, R |K(x, y)|dx c, for every x, y in R

n

.

If 1

p ≤ ∞ and f L

p

(R

n

), then the function T f , defined by

T f (x) =

R

K(x, y) f (y)dy for almost every x in R

n

, belongs to L

p

(R

n

)

and further,

||T f ||

p

c|| f ||

p

.

Proof. If p =

∞, the hypothesis

R |K(x, y)|dx c is superfluous and

the conclusion of the theorem is obvious. If p <

∞, let q denote the

conjugate exponent. Then, by H ¨older ’s inequality,

|T f (x)| ≤

(Z

|K(x, y)|dy

)

1/q

(Z

|K(x, y)| | f (y)|

p

dy

)

1/p

1

background image

2

1. Preliminaries

c

1/q

(Z

|K(x, y)| | f (y)|

p

dy

)

1/p

.



From this we have,

Z

|T f (x)|

p

dx

c

p/q

x

|K(x, y)|| f (y)|

p

dydx

c

1+p/q

Z

| f (y)|

p

dx = c

1+p/q

|| f ||

p
p

.

Therefore

||T f ||

p

c|| f ||

p

.

2

Next, we define the convolution of two locally integrable functions.

Definition 1.2. Let f and g be two locally integrable functions. The
convolution of f and g, denoted by f

g, is define by

( f

g)(x) =

Z

f (x

y)g(y)dy =

Z

f (y)g(x

y)dy = (g f )(x),

provided that the integrals in question exist.

The basic theorem on convolution is the following theorem, called

Young’s inequality.

Theorem 1.3 (Young’s Inequality). Let f

L

1

(R

n

) and g

L

p

(R

n

), for

1

p ≤ ∞. Then f g L

p

(R

n

) and

|| f g||

p

≤ ||g||

p

|| f ||

1

.

Proof. Take K(x, y) = f (x

y) in Theorem 1.1. Then K(x, y) satisfies all

the conditions of Theorem 1.1. and the conclusion follows immediately.



The next theorem underlies one of the most important uses of con-

volution. Before coming to the theorem, let us prove the following

Lemma 1.4. For a function f defined on R

n

and x in R

n

, we define a

function f

x

by f

x

(y) = f (y

x). If f L

p

, 1

p < ∞, then lim

x

O

|| f

x

f

||

p

= 0.

background image

1. General Theorems About Convolutions

3

Proof. If g is a compactly supported continuous function, then g is uni-

3

formly continuous, and so g

x

converges to g uniformly as x tends to

0.



Further, for

|x| ≤ 1, g

x

and g are supported in a common compact

set. Therefore, lim

x

→0

||g

x

g||

p

= 0. Given f

L

p

, we can find a function

g which is continuous and compactly supported such that

|| f g||

p

<

∈ /3

for

∈> 0. But then ||g

x

f

x

||

p

<

∈ /3 also holds. Therefore

|| f

x

f ||

p

≤ || f

x

g

x

||

p

+

||g

x

g||

p

+

||g f ||

p

≤ 2 ∈ /3 + ||g

x

g||

p

.

Since lim x

→ 0||g

x

g||

p

= 0, we can choose x close to 0 so that

||g

x

g||

p

<

∈ /3. Then || f

x

f || <∈ and this proves the lemme since ∈ is

arbitrary.

Remark 1.5. The above lemma is false when p =

∞. Indeed, “ f

x

f

in L

” means precisely that f is uniformly continuous.

Let us now make two important observations about convolutions

which we shall use without comment later on.

i) Supp( f

g) ⊂ Supp f + Supp g, where

A + B =

{x + 1 : x A, y B}.

ii) If f is of class C

k

and ∂

α

(

|α| ≤ k) and g satisfy appropriate condi-

tions so that differentiation under the integral sign is justified, then

f

g is of class C

k

and ∂

α

( f

g) = (∂

α

f )

g.

Theorem 1.6. Let g

L

1

(R

n

) and

R

g(x)dx = a. Let g

(x) =

n

g(x/

∈)

4

for

∈> 0. Then, we have the following:

i) If f

L

p

(R

n

), p <

∞, f g

converges to a f in L

p

as

tends to 0.

ii) If f is bounded and continuous, then f

g

converge to a f uniformly

on compact sets as

tends to 0.

background image

4

1. Preliminaries

Proof. By the change of variable x

→∈ x, we see that

R

g

(x)dx = a for

all

∈> 0. Now

( f

g

)(x)

a f (x) =

Z

f (x

y)g

(y)dy

Z

f (x)g

(y)dy

=

Z

[ f (x

y) − f (x)] g

(y)dy

=

Z

[ f (x

− ∈ y) − f (x)]g(y)dy

=

Z

[ f

y

(x)

f (x)]g(y)dy.



If f

L

p

and p <

∞, we apply Minkoswski’s inequality for integrals

to obtain

|| f g

a f ||

p

Z

|| f

y

f ||

p

|g(y)|dy.

The function y

→ || f

y

f ||

p

is bounded by 2

|| f ||

p

and tends to 0 as

tends to 0 for each y, by lemma 1.4. Therefore, we can apply Lebesgue
Dominated Convergence theorem to get the desired result.

On the other hand, suppose f is bounded and continuous. Let K be

any compact subset of R

n

. Given δ > 0, choose a compact set G R

n

such

that

Z

R

n

G

|g(y)|dy < δ.

5

Then

Sup

x

k

(

| f g

ǫ

)(x)

a f (x)| ≤ 2δ|| f ||

+ Sup

(x,yK

×G

| f (x − ǫy) − f (x)|

Z

G

|g|dy.

Since f is uniformly continuous on the compact set K the second

term tends to 0 as

∈ to 0. Since δ is arbitrary, we see that

sup x

K|( f g

ǫ

)(x)

a f (x)| → 0 as ǫ → 0.

Hence the theorem is proved.

background image

2. The Fourier Transform

5

Corollary 1.7. The space C

o

(R

n

) is dense in L

p

(R

) for 1

p < ∞.

Proof. Let

φ(x) = e

−1/(1−|x|

2

)

for

|x| < 1

= 0

for

|x| ≤ 1

Then φǫC

o

(R

n

) and

R

φ(x)dx = 1/a > 0. If f ǫL

p

and has compact

support, then a( f

∗ φ

ǫ

C

o

(R

n

) and by theorem 1.6, a( f

∗ φ

ǫ

) converges

to f in L

p

as ǫ tends to 0. Since L

p

functions with compact support are

dense in L

p

, this completes the proof.



Proposition 1.8. Suppose K

⊂ R

n

is compact and

K be an open

subset of R

n

. Then there exists a C

o

function φ such that φ(x) = 1 for

xǫK and Supp φ

⊂ Ω.

Proof. Let V =

{x ∈ Ω : d(x, K) ≤

1

2

δ

} where δ = d(K, R

n

\Ω). Choose

6

a φ

o

ǫC

o

such that Supp φ

o

B

0,

1

2

δ

!

and

R

φ

o

(x)dx = 1. Define

φ(x) =

Z

V

φ

o

(x

y)dy = (φ

0

X

V

)(x).

Then φ(x) is a function with the required properties.



2 The Fourier Transform

In this section, we will give a rapid introduction to the theory of the
Fourier transform.

For a function f ǫL

1

(R

n

), the Fourier transform of the function f ,

denoted by ˆ

f , is defined by

ˆ

f (ξ) =

Z

e

−2πix

f (x)dx, ξǫR

n

.

Remark 1.9. Our definition of ˆ

f differs from some other in the place-

ment of the factor of 2π.

background image

6

1. Preliminaries

BASIC PROPERTIES OF THE FOURIER TRANSFORM

For

(1.10)

f ǫL

1

,

|| ˆf||

≤ || f ||

1

.

The proof of this is trivial.
For

(1.11)

f, g

L

1

, ( f

g)

ˆ

(ξ) = ˆf(ξ)ˆg(ξ).

Indeed,

( f

g)

ˆ

(ξ) =

x

e

−2πix

f (y)g(x

y)dydx

=

x

e

−2πi(xy)·ξ

g(x

y)e

−2πiY

f (y)dydx

=

Z

e

−2πi(xy)·ξ

g(x

y)dx

Z

f (y)e

−2πiy

dy

= ˆf(ξ)ˆg(ξ).

Let us now consider the Fourier transform in the Schwartz class S =

7

S (R

n

).

Proposition 1.12. For f

S , we have the following:

i) ˆ

f ǫC

(R

n

) and

β

ˆ

f = ˆg where g(x) = (

−2πix)

β

f (x).

ii) (∂

β

f )

ˆ

(ξ) = (2πiξ)

β

ˆ

f (ξ).

Proof.

i) Differentiation under the integral sign proves this.

ii) For this, we use integration by parts.

(∂

β

f )

ˆ

(ξ) =

Z

e

−2πix

(∂

β

f )(x)dx

= (

−1)

|β|

Z

β

[e

−2πix

] f (x)dx

= (

−1)

|β|

(

−2πiξ)

β

Z

e

−2πix·ξ

f (x)dx

= (2πiξ)

β

ˆ

f (ξ).



background image

2. The Fourier Transform

7

Corollary 1.13. If f ǫS , then ˆ

f

S also.

Proof. For multi-indices α and β, using proposition 1.12, we have

ξ

α

(∂

β

ˆ

f )(ξ) = ξ

α

((

−2πix)

β

f (x))

ˆ

(ξ)

= (2πi)

−|α|

[∂

α

((

−2πix)

β

f (x))

ˆ

(ξ)]

= (

−1)

|β|

(2πi)

|β|−|α|

(∂

α

(x

β

f (x)))

ˆ

(ξ)

Since f ǫS , ∂

α

(x

β

f (x))ǫL

1

and hence (∂

α

(x

β

f (x)))

ˆ

ǫL

. Thus ξ

α

(∂

β

ˆ

f ) is bounded. Since α and β are arbitrary, this proves that ˆ

f ǫS .



Corollary 1.14 (RIEMANN-LEBESGUE LEMMA). If f ǫL

1

, then ˆ

f is

continuous and vanishes at

.

Proof. By Corollary 1.13, this is true for f ǫS . Since S is dense in L

1

8

and

|| ˆf||

≤ || f ||

1

, the same in true for all f ǫL

1

.



Let us now compute the Fourier transform of the Gaussian.

Theorem 1.15. Let f (x) = e

−πa|x|

2

, Re a > 0. Then, ˆ

f (ξ) = a

n/2

e

a

−1

π

|ξ|

2

.

Proof.

ˆ

f (ξ) =

Z

e

−2πix

e

aπ|x|

2

dx

i.e.,

ˆ

f (ξ) =

n

Y

j=1

Z

−∞

e

−2πix

j

ξ

j

e

aπx

2

j

dx

j

.

Thus it suffices to consider the case n = 1. Further, we can take

a = 1 by making the change of variable x

a

−1/2

x.



Thus we are assuming f (x) = e

−πx

2

, xǫR. Observe that f

(x) +

x f (x) = 0. Taking the Fourier transform, we obtain

iξ ˆ

f (ξ) + i ˆf

′′

(ξ) = 0.

Hence

ˆf

′′

(ξ)/ f (ξ) =

−2πξ

background image

8

1. Preliminaries

which, on integration, gives f (ξ) = ce

−πξ

2

, c being a constant.

The constant c is given by

c = ˆf(0) =

Z

−∞

e

−πx

2

dx = 1.

Therefore ˆ

f (ξ) = e

−πξ

2

, which completes the proof.

We now derive the Fourier inversion formula for the Schwartz class

S .

Let us define f

V

(ξ) =

R

e

ix

·ξ

f (x)dx = ˆf(

−ξ).

9

Theorem 1.16 (Fouries Inversion Theorem). For f ǫ( ˆ

f )

= f .

Proof. First, observe that for f, gǫL

1

,

R

f ˆg =

R

ˆ

f g. In fact,

Z

ˆ

f (x)g(x)dx =

x

e

−2πiy·x

f (y)g(x)dydx

=

Z

"Z

e

−2πiy·x

g(x)dx

#

f (y)dy

=

Z

ˆg(y) f (y)dy.



Given ǫ > 0 and x in R

n

, take the function φ defined by

φ(ξ) = e

−2πix·ξ−πǫ

2

|ξ|

2

.

Now

ˆ

φ(y) =

Z

e

−2πiy·ξ

e

−2πix·ξ−πǫ

2

|ξ|

2

dξ

=

Z

e

−2πi(yx)·ξ

e

−πǫ

2

|ξ|

2

dξ

=

n

e

−π∈

−2

|xy|

2

.

If we take g(x) = e

−π|x|

2

and define g

ǫ

(x) = ǫ

n

g(x/ǫ), then

ˆ

φ(y) = g

ǫ

(x

y).

background image

2. The Fourier Transform

9

Therefore

Z

e

ix

·ξ

ˆ

f (ξ)e

−πǫ

2

|ξ|

2

dξ =

Z

ˆ

f (ǫ)φ(ξ)dξ

=

Z

f (y) ˆ

φ(y)dy

=

Z

f (y)g

ǫ

(x

y)dy

= ( f

g

ǫ

)(x)

But as ǫ tends to 0, ( f

g

ǫ

) converges to f , by Theorem 1.6 and

10

clearly

Z

e

i

·xξ

ˆ

f (ξ)e

−πǫ

2

|ξ|

2

dξ

Z

ˆ

f (ξ)e

i

·x

dξ

Therefore ( ˆ

f )

v

= f .

Corollary 1.17. The Fourier transform is an isomorphism of S onto S .

Next, we prove the Plancherel Theorem.

Theorem 1.18. The Fourier transform uniquely extends to a unitary
map of L

2

(R

n

) onto itself.

Proof. For f ǫS , define ˜

f (x) = f (

x). Then it is easily checked that

ˆ˜f = ¯ˆf, so that

|| f ||

2
2

=

Z

| f (x)|

2

dx

=

Z

f (x) ˜

f (

x)dx

= ( f

∗ ˜f)(0)

=

Z

( f

˜

f )

ˆ

(ξ)dξ

=

Z

ˆ

f (ξ) ˆ˜

f (ξ)dξ

=

Z

ˆ

f (ξ) ¯ˆ

f (ξ)dξ =

|| ˆf||

2
2

.

background image

10

1. Preliminaries

Therefore, the Fourier transform extends continuously to an isom-

etry of L

2

. It is a unitary transformation, since its image S is dense in

L

2

.



Let us observe how the Fourier transform interacts with translations,

rotations and dilations.
(1.19)

The Fourier transform and translation: If f

x

(y) f (y

x) then

11

ˆ

f

x

(ξ) =

Z

e

iy

·ξ

f (y

x)dy

=

Z

e

i(z+x)

·ξ

f (z)dz ( by putting y

x = z)

=

Z

e

ix

·ξ

ˆ

f (ξ).

(1.20)

The Fourier transform and rotations (orthogonal transfor-

mations):

Let T : ’

n

→ ’

n

be an orthogonal transformation. Then

( f

T )

ˆ

(ξ) =

Z

e

−2πix·ξ

( f oT )(x)dx

=

Z

e

−2πiT

−1

y

·ξ

f (y)dy ( by putting y = T x)

=

Z

e

−2πiy·T ξ

f (y)dy

= ˆf(T ξ) = ( ˆfoT )(ξ).

Thus, ( f oT )

ˆ

= ˆfoT i.e. ˆ commutes with rotations.

(1.21)

The Fourier transform and dilation: Let f

r

(x) = r

n

f (x/r).

Then

ˆ

f

r

(ξ) =

Z

e

−2πix·ξ

r

n

f (x/r)dx

=

Z

e

−2πiy·ξ

f (y)dy = ˆf(rξ).

The last equation suggests, roughly: the more spread out f is, the

more ˆ

f will be concentrated at the origin and vice versa. This notion

can be put in a precise form as follows.

background image

2. The Fourier Transform

11

(1.22)

HEISENBERG INEQUALITY (n = 1: For f ǫS (R), we have

||x f (x)||

2

||ξ ˆf(ξ)||

2

≥ (1/4π)|| f ||

2
2

.

Proof. Observe that

12

d

dx

(x f (x)) = x

d f

dx

(x) + f (x).

Thus

|| f ||

2
2

=

Z

f (x) f (x)dx

=

Z

f (x)

" d

dx

(x f (x))

x

d f

dx

(x)

#

dx

=

Z

x f (x)

d ¯

f

dx

(x)dx

Z

x

d f

dx

(x) ¯

f (x)dx

=

−2Re

Z

x f (x)

d ¯

f

dx

(x)dx.

≤ 2||x f (x)||

2

||

d ¯

f

dx

||

2

(by Cauchy-Schwarz)

i.e.,

|| f ||

2
2

≤ 2||x f (x)||

2

||

d f

dx

||

2

.

But

(

d ¯

f

dx

)

ˆ

(ξ) = 2πiξ ˆf(ξ).

Therefore,

|| f ||

2
2

≤ 2.2π||x f (x)||

2

||ξ ˆf(ξ)||

2

or

||x f (x)||

2

||ξ f (ξ)||

2

≥ (1/4π)|| f ||

2
2

.



background image

12

1. Preliminaries

A GENERALISATION TO n VARIABLES

We replace x by x

j

,

d

dx

by

x

j

. Also for any a

j

, b

j

ǫ’, we can re-

place x

j

and

x

j

by x

j

a

j

and

x

j

b

j

respectively. The same proof

then yields:

(1.23)

||(x

j

a

j

) f (x)

||

2

||(ξ

j

b

j

) ˆ

f (ξ)

||

2

≥ (1/4π)|| f ||

2
2

.

13

Let us now take

|| f ||

2

= 1 and A = (a

1

, a

2

, . . . , a

n

)ǫ’

n

. Let f be

small outside a small nighbourhood of A.

In this case,

||(x

j

a

j

) f (x)

||

2

will be small. Consequently, the other

factor on the left in (1.23) has to be large. That is, if the mass of f is
concentrated near one point, the mass of ˆ

f cannot be concentrated near

any point.

Remark 1.24. If we take

a

j

=

Z

x

i

| f (x)|

2

dx, b

j

=

Z

ξ

j

| ˆf(ξ)|

2

dξ,

then inequality (1.23) is the mathematical formulation of the position-
momentum uncertainty relation
in Quantum Mechanics.

3 Some Results From the Theory of Distributions

In this section, let us recall briefly some results from the theory of dis-
tributions. (For a more detailed treatment, see, for example [7] or [8]).

In the sequel,

D

(Ω) will denote the space of distributions on the

open set Ω

⊂ R

n

which is the dual space of C

o

(Ω). When Ω = ’

n

, we

will simply write

D

instead of

D

n

). In the same way, S

= S

(R

n

)

will denote the space of tempered distributions with and E

= E

n

)

will stand for the space of distributions with compact support.

The value of a distribution uǫ

D

at a function φǫC

o

will be denoted

by < u, φ >. If u is a locally integrable function, then u defines a distri-
bution by < u, φ >=

R

u(x)φ(x)dx. It will sometimes be convenient to

14

background image

3. Some Results From the Theory of Distributions

13

write

R

u(x)φ(x)dx for < u, φ >, when u is an arbitrary distribution.

The convergence in

D

is the weak convergence defined by the fol-

lowing:

u

n

, uǫ

D

, u

n

→ in means < u

n

, φ >

→< u, φ > for every φ in C

o

.

Let us now recall briefly certain operations on distributions.

(1.25)

We can multiply a distribution uǫ

D

by a C

function φ to

get another distribution φu which is defined by

< φu, ψ >=< u, φψ >

A C

function ψ is said to be tempered if, for every multiindex

α, ∂

α

φ grows at most polynomially at

∞. We can multiply an u, ǫS

by

a tempered function to get another tempered distribution. The definition
is same as in the previous case.
(1.26)

If uǫ

D

and f ǫC

o

, we define the convolution u

f by

(u

f )(x) =< u, f

x

> where f

x

(x) = f (x

y).

The function u

f is C

and when uǫE, u

f is in C

o

. The convo-

lution

∗ : D

× C

o

C

can be extended to a map from

D

× E

to

D

.

Namely, if uǫ

D

, vǫE

and φǫC

o

, < u

v, φ >=< u, ˜v ∗ φ > where ˜v is

defined by < ˜v, ψ >=

R

v(x)ψ(

x)dx. The associative law

u

∗ (v w) − (u w) ∗ w holds for u, v, wǫD

provided that at most one of them does not have compact support. For

15

uǫS

and f ǫS , u

f can also defined in the same way and u f is a

tempered C

function.

(1.27)

Since the Fourier transform is an isomorphism of S onto

S . and

R

f ˆg =

R

ˆ

f g, the Fourier transform extends by duality to an

isomorphism of S

onto S

.

For uǫE

S

, we have ˆu(ξ) =< u, e

−2πi(·)·ξ

> which is an entire

analytic function.

background image
background image

Chapter 2

Partial Di

fferential

Operators with Constant
Coe

fficients

1 Local Solvability and Fundamental Solution

For the sake of convenience in taking the Fourier transforms, from now

16

onwards, we will use the differential monomials

D

α

= (2πi)

|α|

α

. Thus (D

α

)

ˆ

(ξ) = ξ

α

ˆ

f (ξ).

By a partial differential operator with constant coefficients,we mean

a differential operator L of the form

L =

X

|α|≤k

a

α

D

α

, a

α

ǫƒ

Further, we assume that

P

|α|≤k

|a

α

| , 0. In this case, we say that the

operator L is of order k. If we write P(ξ) =

P

|α|≤k

a

α

ξ

α

, then we have

L = P(D).

For uǫS , taking the Fourier transform, we see that

(P(D)u)

ˆ

(ξ) =

X

|α|≤k

a

α

ξ

α

ˆu(ξ) = P(ξ)ˆu(ξ).

15

background image

16

2. Partial Differential Operators with Constant Coefficients

Let us now consider the following problem:
Given f in C

, we want to find a distribution u such that P(D)u = f .

We say that the differential operator L is locally solvable at x

ǫ’

n

,

if there is a solution of the above problem in some neighbourhood of the
point x

for any f in C

.

Remark 2.1. We can assume that f has compact support. To see this,

17

we can take any φ in C

o

such that φ = 1 in some nighbourhood of the

point x

o

. If solve the problem P(D)u = f φ near x

, then u is a solution

of our original problem since we have f φ = f in a neighbourhood of x

o

.

In the following theorem we give an affirmative answer to the ques-

tion of local solvability of L = P(D). The simple proof exhibited here is
due to L. Nirenberg.

Theorem 2.2. Let L =

P

|α|≤k

a

α

D

α

be a differential operator with con-

stant coefficients. If uǫC

o

, there is a C

function u satisfying Lu = f

on ’

n

.

Proof. Taking the Fourier transform of the equation Lu = P(D)u = f ,
we see that P(ξ)ˆu(ξ) = ˆf(ξ). It is natural to try to define u by the formula

u(x) =

Z

e

ix

·ξ

ˆ

f (ξ)/P(ξ)ξ.

In general, P will have many zeros: so there will be a problem in ap-
plying the inverse Fourier transform to ˆ

f /P. But things are not so bad.

Since f ǫC

, ˆ

f is an entire function of ξǫƒ

n

and P is obviously entire.

Hence we can deform the contour of integration to avoid the zeros of
P(ξ).



To make this precise, let us choose a unit vector η so that

P

|α|≤k

a

α

η

α

,

0. By a rotation of coordinates, we can assume that η = (0, 0, . . . , 0, 1).
Multiplying by a constant, we can also assume that a

α0

= 1 where

18

α

o

= (0, 0, . . . , 0, k). Then we have P(ξ) = ξ

k

n

+(lower order terms in

ξ

n

). Denote ξ = (ξ

, ξ

n

) with ξ

= (ξ

1

, . . . , ξ

n

−1

) in ’

n

−1

.

background image

1. Local Solvability and Fundamental Solution

17

Consider P(ξ) = P

, ξ

n

) as a polynomial in the last variable ξ

n

in

C

for ξ

in R

n

−1

. Let λ

1

), . . . , λ

k

) be the roots of P

, ξ

n

) = 0

arranged so that if i

j.

ℑλ

i

)

≤ ℑλ

j

) and Re λ

i

)

≤ Re λ

i

) when

ℑλ

i

) =

ℑλ

j

).

Since the roots of a polynomial depends continuously on the coef-

ficients we see that λ

j

) are continuous in ξ

. To proceed further, we

need the help of two Lemmas.

Lemma 2.3. There is a measurable function φ : ’

n

−1

→ [−k − 1, k + 1]

such that for all ξ

in ’

n

−1

min

1

jk

{|φ(ξ

)

Imλ

j

)

|} ≥ 1.

Proof. Left as an exercise to the reader : (cf. G.B. Folland [1]).

(The idea is that at least one of the k+1 intervals [

k−1, k+1], [−k

1, k+1] . . . , [

k−1, k+1] must contain none of the k points Im λ

j

), j =

1, 2, . . . , k).



Lemma 2.4. Let P(ξ) = ξ

k

n

+( lower order terms) and let N(P) be the set

{ζǫƒ

n

: P(ζ) = 0

}. Let d(ξ, N(P)). Then, we have

|P(ξ)| ≥ (d(ξ)/2)

k

.

Proof. Take ξ in R

n

such that P(ξ) , 0. Let η = (0, 0, . . . 0, 1) and

define g(z) = P(ξ + zη) for z in C. This g is a polynomial in one complex
variable z. Let λ

1

, λ

2

, . . . , λ

k

be the zeros of the polynomial g. Then

19

g(z) = c(z

− λ

1

)

· · · (z − λ

k

)

so that

|

g(z)

g(0)

| =

k

Y

j=1

|1 −

z

λ

j

|.

Since ξ + λ

j

η ǫ N(P)

j

| ≥ d(ξ) so that when |z| ≤ d(ξ), |

g(z)

g(0)

| ≤ 2

k

. Also

|g

(k)

(0)

| = |

k!

i

Z

|ζ|=d(ξ)

q(ζ)ζ

k−1

dζ

| ≤

k!

2

|g(0)|

d(ξ)

k+1

2

k

d(ξ)

background image

18

2. Partial Differential Operators with Constant Coefficients

i.e.

|g

(k)

(0)

| ≤ k!|g(0)|2

k

d(ξ)

k

. But g(0) = P(ξ) and

|g

(k)

(0)

| =

k

∂ξ

k

n

P(ξ) = k!

Therefore,

k!

k!|g(0)|2

k

d(ξ)

k

or

|P(ξ)| ≥ (d(ξ)/2)

k

Hence the lemma is proved.



Returning to the proof of the theorem, consider the function

u(x) =

Z

R

n

−1

Z

IMξ

n

=φ(ξ

)

e

ix

·ξ

( ˆ

f (ξ)/P(ξ))dξ

n

dξ.

By Lemmas 2.3 and 2.4 we have

|P(ξ)| ≥ (d(ξ)/2)

k

≥ 2

k

along Imξ

n

= φ(ξ

).

Since f ǫC

, ˆ

f (ξ) is rapidly decreasing as

|Reξ| tends to ∞ as along

20

as

|Imξ| stays bounded, so the integral converges absolutely and uni-

formly together with all derivatives defining a C

function u.

Finally, by Cauchy’s theorem,

P(D)u(x) =

Z

R

n

−1

Z

IMξ

n

=φ(ξ

)

e

ix

·ξ

( ˆ

f (ξ)/P(ξ))dξ

n

dξ

=

Z

R

n

e

ix

·ξ

( ˆ

f (ξ)dξ = f (x).

This completes the proof of Theorem 2.2.
Let us now consider the local solvability of L = P(D) in the case

when f is a distribution.

As before, we remark that it suffices to take f ǫE

. Further, it is

enough to consider the case where f = δ. Indeed, if K satisfies P(D)K =
δ, then we have for any f ǫE

, P(D)(K

f ) = P(D)K f = δ ∗ f = f .

background image

1. Local Solvability and Fundamental Solution

19

Definition 2.5. A distribution K satisfying P(D)K = δ is called a fun-
damental solution or elementary solution of the di
fferential operator
L
= P(D).

A remarkable theorem due to MALGRANGE AND EHRENPREIS

states that every differential operator P(D) with constant coefficients has
a fundamental solution. In fact, we can prove this result by a simple
extension of the preceding argument.

Theorem 2.6. Every partial differential operator P(D) with constant co-

21

efficients has a fundamental solution.

Proof. Proceeding as in the previous theorem, we try to define

K(x) =

Z

’

n

−1

Z

Imξ

n

=φ(ξ

)

e

ix

·ξ

(p(ξ))

−1

dξ

n

dξ

.

Here, however, the integral may diverge at infinity. So, we consider

the polynomial

P

N

(ξ) = P(ξ)(1 + 4Π

2

n

X

j=1

ξ

2

j

)N

where N is a large positive integer. Let

K

N

(x) =

Z

’

n

−1

Z

Imξ

n

=φ(ξ

)

e

ix

·ξ

(P

N

(ξ))

−1

dξ

n

dξ

where φ is chosen appropriately for the polynomial P

N

.Then on the re-

gion of integration, we have

|P

N

(ξ)

| ≥ c(1 + |ξ|

2

)

N

;

so the integral will converge when N > n/2. We claim that P

N

(D)K

N

=

δ.



To see this, take φǫC

o

and observe that the transpose of P

N

(D) is

P

N

(

D). Thus

< P

N

(D)K

N

, φ > =< K

N

, P

N

(

D)φ >

background image

20

2. Partial Differential Operators with Constant Coefficients

=

Z

’

n

Z

’

n

−1

Z

Imξ

n

=φ(ξ

)

e

ix

·ξ

(P

N

(

D)φ)(x)

P

N

(ξ)

dξ

n

dξ

dx

=

Z

’

n

−1

Z

Imξ

n

=φ(ξ

)

(P

N

(ξ))

−1

dξ

n

dξ

Z

’

n

e

ix

·ξ

P

N

(

D)φ(x)dx

=

Z

’

n

−1

Z

Imξ

n

=φ(ξ

)

ˆ

φ(

−ξ)dξ

n

dξ

=

Z

’

n

ˆ

φ(

−ξ)dξ = φ(0) =< δ, φ > .

22

Thus,

δ = P

N

(D)K

N

= P(D)(1

− ∆)

N

K

N

so if we put K = (1

− ·∆)

N

K

N

we have P(D)K = δ. Thus K is a

fundamental solution of P(D).

2 Regularity Properties of Di

fferential Operators

Definition 2.7. The singular support of a distribution f ǫD

is defined to

be the complement of the largest open set on which f is a C

function.

The singular support of f will be denoted by sing supp f .

Definition 2.8. Let L =

P

|α|≤K

a

α

(x)D

α

where a

α

ǫC

be a differential

operator. L is said to be hypoelliptic, if and only if, for any uǫ

D

, sing

supp u

sing supp Lu. In other words, L is hypoelliptic if and only

if for any open set

⊂ ’

n

and any uǫ

D

(Ω),

LuǫC

(Ω) implies

uǫC

(Ω)”.

Remark 2.9. The operator L is said to be elliptic at xǫ’

n

if

X

|α|=k

a

α

(x

α

,

0 for every ξ in ǫ’

n

\{0}.

L is said to be elliptic if L is elliptic at every xǫ’

n

. Elliptic operators

are hypoelliptic, as we shall prove later. This accounts for the name

23

hypoelliptic.

background image

2. Regularity Properties of Differential Operators

21

We know that an ordinary differential operator with C

coefficients

is hypoelliptic as long as the top order coefficient is non-zero. But this is
not the case with partial differential operators as seen from the following

Example 2.10. Take the operator L =

2

xy

in ’

2

. The general solution

of the equation Lu = 0 is given by u(x, y) = f (x) + g(y) where f and g
are arbitrary C

1

functions. This shows that L is not hypoelliptic.

We observe that if L is hypoelliptic, then every fundamental solution

of L is a C

function in ’

n

\{0}.

In the case of partial differential operators with constant coefficients,

this is also sufficient for hypoellipticity. Indeed, we have the following

Theorem 2.11. Let L be a partial differential operator with constant
coefficients. Then the following are equivalent:

a) L is hypoelliptic.

b) Every fundamental solution of L is C

in’

2

\{0}.

c) At least one fundamental solution of L is C

in ’

2

\{0}.

Proof. That (a) simple (b) follows from the above observation and that
(b) implies (c) is completely trivial. The only nontrivial part we need to
prove is that (c) implies (a). To prove this implication, we need



Lemma 2.12. Suppose f ǫ

D

is such that f is C

in ’

n

\{0} and gǫE

.

24

Then sing supp( f

g) ⊂ supp g.

Proof. Suppose x < supp g. We will show that f

g is C

in a neigh-

bourhood of x.



Since x < supp g, there exists an ǫ > 0 such that B(x, ǫ)

∩supp g = φ.

Choose φǫC

0

(B(0, ǫ/2)) such that φ = 1 on B(0, ǫ/4). Now, f

g =

f )

g + (1 − φ) f g. Since (1 − φ) f is a C

function, (1

− φ) f g is

C

. Also

Supp(φ f

g) ⊂ Supp φ f + Supp g

background image

22

2. Partial Differential Operators with Constant Coefficients

⊂ {y : d(y, supp g) ≤ ǫ/2}

which does not intersect B(x, ǫ/2). Therefore, in B(x, ǫ/2), f

g = (1 −

φ) f

G which is C

.

Proof of theorem 2.11. Let K be a fundamental solution of L such that
K is C

in ’

n

\{0}. Suppose uǫD

and LuǫC

(Ω) where Ω

⊂ ’

n

is

open.

For xǫΩ, pick ǫ > 0 small enough so that ¯B(x, ǫ)

⊂ Ω. Choose

φǫC

o

(B(x, ǫ)) so that φ = 1 on B(x, ǫ/2). Then Lu) = φLu + v where

v = 0 on B(x, ǫ/2) and also outside B(x, ǫ). We write

K

Lu) = K ∗ φLu + K v.

φLu is a C

o

function so that K

Lu is a C

function. Also K

v is

a C

function on the ball B(x, ǫ/2) by Lemma 2.12.

Therefore K

Lu) is a C

function on B(x, ǫ/2). But K

Lu) =

25

LK

∗ φu = δ ∗ φu = φu. Thus, φu is a C

function on B(x, ǫ/2). Since

φ = 1 on B(x, ǫ/2), u is a C

function on B(x, ǫ/2). Since x is arbitrary,

this completes the proof.

3 Basic Operators in Mathematical Physics

In this section, we introduce the three basic operators in Mathematical
Physics. In the following sections, we shall compute fundamental so-
lutions for these operators and show how they can be applied to solve
boundary value problems and to yield other information.

(i) LAPLACE OPERATOR: ∆ =

n

P

j=1

2

x

2

j

in ’

n

. If u(x) represents

electromagnetic potential (or gravitational potential ) and ρ de-
notes the charge (resp. mass) density, then they are connected by
the equation ∆u =

−4πρ. If the region contains no charge, i.e., if

ρ = 0, ∆u = 0. This equation is called the homogeneous Laplace
equation
, and its solutions are called harmonic functions.

background image

3. Basic Operators in Mathematical Physics

23

(ii) HEAT OPERATOR: L =

t

− ∆ in ’

n+1

. If u(x, t) represents the

temperature of a homogeneous body at the position x and time t,

then u satisfies the heat equation

u

t

− ∆u = 0 in ’

n+1

.

(iii) WAVE OPERATOR: @ =

2

t

2

− ∆ in ’

n+1

. If u(x, t) represents the

amplitude of an electromagnetic wave in vacuum at position x and
time t, then u satisfies the wave equation @u = 0 in ’

n+1

.

26

The equation u = 0 can also be used to describe other types of wave
phenomena although in most cases, it is only an approximation valid for
small amplitudes. More generally, the equation u = f describes waves
subject to a driving force f .

The Laplace operator ∆ is an ingredient in all the above examples.

The reason for this is that the basic laws of Physics are invariant un-
der translation and rotation of coordinates which severely restricts the
differential operators which can occur in them. Indeed, we have

Theorem 2.13. Let L be a differential operator which is invariant under
rotations and translations. Then there exists a polynomial Q in one
variable with constant coe
fficients such that L = Q(∆).

Proof. Let L = P(D) =

P

|α|≤K

a

α

D

α

. Since P(D) is invariant under trans-

lations, a

α

are constants. Let T be any rotation. We have

P(T ξ) ˆ

φ(T ξ) = (P(D)φ)ˆ(T ξ)

= (P(D

T )ˆ(ξ)

= (P(D

T ))ˆ(ξ)

= P(ξ) ˆφ(T ξ), for every φ.

Thus P(T ξ) = P(ξ) so that P is rotation - invariant.



Write P(ξ) =

k

P

j=0

P

j

(ξ) where P

j

(ξ) is the part of P which is ho-

27

mogeneous of degree j. We claim that each P

j

is rotation - invariant.

Indeed, if t is any real number,

background image

24

2. Partial Differential Operators with Constant Coefficients

P(tξ) =

k

P

j=0

t

j

P

j

(T ξ). For a rotation T , since P(ξ) = P(T ξ), we have

P(tξ) =

k

X

j=0

t

j

P

j

(T ξ)

so that

k

X

j=0

t

j

(P

j

(T ξ)

P

j

(ξ)) = 0.

This being true for all t, it follows that P

j

T = P

j

. But the only

rotation-invariant functions which are homogeneous of degree j are of
the form P

j

(ξ) = c

j

|ξ|

j

where c

j

are constants.

Indeed, since P

j

is rotation-invariant, P

j

(ξ) depends only on

|ξ| and

thus, for

|ξ| , 0,

P

j

(ξ) = P

j

|ξ|

ξ

|ξ|

!

=

|ξ|

j

P

j

ξ

|ξ|

!

= c

j

|ξ|

j

.

Since

|ξ|

j

is not a polynomial when j is odd, c

j

= 0 in that case.

Therefore,

P(ξ) =

X

c

2

j

|ξ|

2 j

.

Taking Q(x) =

P(c

2 j

/(

−4π

2

)

j

)x

j

, we get P(D) = Q(∆).

Remark 2.14. This theorem applies to scalar differential operators. If
one considers operators on vector or tensor functions, there are first or-
der operators which are translation - and rotation - invariant, namely, the

28

familiar operators grad, curl, div of 3-dimensional vector analysis and
their n-dimensional generalisations.

Definition 2.15. A function F(x) is said to be radial , if there is a
function f of one variable such that F
(x) = f (

|x|) = f (r), r = |x|.

When F is radial and FǫL

1

n

), we have

Z

’

n

F(x)dx =

Z

0

Z

|x|=1

f (r)r

n

−1

dσ(x)dr

background image

4. Laplace Operator

25

where dσ(x) is the surface measure on S

n

−1

, the unit sphere in ’

n

.

Thus

Z

’

n

F(x)dx = ω

n

Z

0

f (r)r

n

−1

dr

where ω

n

is the area of S

n

−1

.

Let us now calculate ω

n

. We have

Z

e

−π|x|

2

dx = 1.

Now

Z

e

−π|x|

2

dx = ω

n

Z

0

e

−πr

2

r

n

−1

dr

=

ω

n

Z

0

e

s

(s/π)

n/2

−1

ds, s = πr

2

=

ω

n

n/2

Z

0

e

s

s

n/2

−1

ds

= ω

n

Γ(n/2)/(2π

n/2

).

Therefore

ω

n

= 2π

n/2

/Γ(n/2).

4 Laplace Operator

First, let us find a fundamental solution of ∆, i.e., we want to find a

29

distribution K such that ∆K = δ.

Since ∆ commute with rotations and δ has the same property, we

observe that if u is a fundamental solution of ∆ and T is a rotation,
then u

T is also a fundamental solution. We therefore expect to find a

fundamental solution K which is radial.

background image

26

2. Partial Differential Operators with Constant Coefficients

Let us tarry a little to compute the Laplacian of a radial function F.

LetF(x) = f (r), r =

|x|. Then,

F(x) =

n

X

j=1

x

j

[ f

(r)x

j

/r]

=

n

X

j=1

f

′′

(r)

x

2

j

r

2

+

f

(r)

r

f

(r)

r

3

x

2

j

= f

′′

(r) + ((n

− 1)/r) f

(r).

Now set K(x) = f (r). If K is to be a fundamental solution of ∆, we

must have

f

′′

(r) + ((n

− 1)/r) f

(r) = 0 on (0,

∞).

From this equation

f

′′

(r)/ f

′′

(r) =

−(n − 1)/r.

Integrating, we get

f

(r) = c

r

1

n

with a constant c

. One more integration yields

f (r) = c

1

r

2

n

+ c

2

when n , 2

= c

1

log r + c

2

when n = 2.

30

Since constants are solutions of the homogeneous Laplace equation,

we may assume that c

2

= 0. Thus if we set F(x) =

|x|

2

n

(n , 2) or

F(x) = log

|x|(n = 2), we expect to find that ∆F = cδ for some c , 0,

and then K = c

−1

F will be our fundamental solution.

In fact, we have

Theorem 2.16. If F(x) =

|x|

2

n

on ’

n

(n , 2), then F = (2

n

n

δ,

where ω

n

is the area of the unit sphere in ’

n

,

background image

4. Laplace Operator

27

Proof. For ǫ > 0, we define F

ǫ

(x) = (ǫ

2

+

|x|

2

)

(2

n)/2

. Then F

ǫ

is a

C

function on ’

n

, and an application of Lebesgue’s Dominated Con-

vergence Theorem reveals that F

ǫ

converges to F as ǫ tends to 0, in

the sense of distributions. Therefore ∆F

ǫ

converges to ∆F, in the same

sense. Let us now compute ∆F

ǫ

.

F

ǫ

(x) =

n

X

j=1

x

j

n

(2

n)(|x|

2

+ ǫ

2

)

n/2

x

j

o

=

n

X

j=1

n

(2

n)(|x|

2

+ ǫ

2

)

n/2

+ (2

n)(−n)(|x|

2

+ ǫ

2

)

(

n/2)−1

x

2

j

o

= (2

n)(|x|

2

+ ǫ

2

)

(

n/2)−1

nǫ

2

.

Thus we see that ∆F

ǫ

ǫL

1

n

). A simple computation shows that

F

ǫ

(x) = ǫ

n

F

1

(x/ǫ); so, by Theorem 1.6, F

ǫ

tends to (

R

F

1

)δ as

ǫ tends to 0. Therefore ∆F = (

R

F

1

)δ, and we need only to compute

R

F

1

.

Z

F

1

= n(2

n)

Z

(1 +

|x|

2

)

(

n/2)−1

dx

= n(2

n

n

Z

0

(1 + r

2

)

(

n/2)−1

r

n

−1

dr.

Putting r

2

+ 1 = s, we see that

31

Z

F

1

=

1

2

n(2

n

n

Z

1

s

(

n/2)−1

(s

− 1)

(n

−1)/2

ds

=

1

2

n(2

n

n

Z

1

s

−2

(1

1

s

)

(n/2)

−1

ds

=

1

2

n(2

n

n

Z

0

(1

− σ)

(n/2)

−1

dσ, σ =

1

s

=

1

2

n(2

n

n

(2/n) = (2

n

n

background image

28

2. Partial Differential Operators with Constant Coefficients

which completes the proof.



Exercise. Show that if F(x) = log

|x| on ’

2

, then F = 2πδ.

Corollary 2.17. Let K(x) =


|x|

2

n

(2

n

n

(n , 2)

1

log

|x|(n = 2)

Then K is a fundamental solution of the Laplacian.

Corollary 2.18. is hypoelliptic.

Proof. Follows from Theorem 2.11.



It is also instructive to compute the fundamental solution of ∆ by the

Fourier transform method.

If K is a fundamental solution of ∆, we have ∆(k

f ) = ∆K f = f .

Since (∆g)ˆ(ξ) =

−4π

2

|ξ|

2

ˆg(ξ), we get

−4π

2

|ξ|

2

ˆ

K(ξ) ˆ

f (ξ) = ˆf(ξ),

so that, at least formally,

32

ˆ

K(ξ) =

−1/(4π

2

|ξ|

2

).

We observe that when n > 2, the function

−1/(4π

2

|ξ|

2

) is locally

integrable and so defines a tempered distribution. We want to show that
its inverse Fourier transform is our fundamental solution K. For that
purpose, we will prove a more general theorem.

Theorem 2.19. For 0 < α < n, let F

α

be the locally integrable function

F

α

(ξ) =

|ξ|

−α

. Then

F

α

(x) =

Γ

1

2

(n

− α)

!

1

2

α

!!

π

α

n/2

|x|

α

n

.

Before proving this theorem, let us pause a minute to observe that it

implies

F

2

(x) = Γ((n/2)

− 1)π

2

n/2

|x|

2

n

= (Γ(n/2)/(n/2

− 1))π

2

n/2

|x|

2

n

= (

−4π

2

/((2

n

n

))

|x|

2

n

, so that (

F

2

/4π

2

))

= K

as desired.

background image

4. Laplace Operator

29

Proof of theorem 2.19. The idea of the proof is to express F

α

as a

weighted average of Gaussian functions, whose Fourier transforms we
can compute.

To begin with,

Z

0

e

rt

t

α

−1

dt =

Z

0

e

s

s

α

−1

r

−α

ds = r

−α

Γ(α), (r > 0, α > 0).

In other words, for any r > 0 and α > 0,

33

r

−α

=

1

Γ(α)

Z

0

e

rt

t

α

−1

dt.

Taking r = π

| ξ |

2

and replacing α by α/2, then

| ξ |

−α

=

π

α/2

Γ(α/2)

Z

0

e

−π|ξ|

2t

t

(α/2)

−1

dt

which is the promised formula for F

α

as a weighted average of Gaus-

sians. Formally, we can write

Z

R

n

e

ix

| ξ |

−α

dξ =

π

α/2

Γ(α/2)

Z

R

n

Z

0

e

ix

e

−π|ξ|

2t

t

α/2

−1

dtdξ

=

π

α/2

Γ(α/2)

Z

0

(e

π

|ξ|

2t

)

v

(x)t

(α/2)

−1

dt

=

π

α/2

Γ(α/2)

Z

0

e

−π(|x|

2

)/t

t

n/2

t

α/2

−1

dt

=

π

α/2

Γ(α/2)

Z

0

e

−π|x|

2

S

S

(n

−α)/2−1

ds

=

π

α/2

Γ(α/2)

Γ((n

− α)/2)

n

(n

−α)/2

| x |

α

n

.

background image

30

2. Partial Differential Operators with Constant Coefficients

In this computation, the change of order of the integration is unfor-

tunately not justified, because the double integral is not absolutely con-
vergent. This is not surprising, since F

α

is not an L

1

function; instead,

we must use the definition of the Fourier transform for distributions.

For every φǫS ,

< e

−π|ξ|

2t

, ˆ

φ > =< e

−π|ξ|

2t

)ˆ,φ >,

i.e.,

Z

R

n

e

−π|ξ|

2t

ˆ

φ(ξ)dξ =

Z

R

n

t

n/2

e

−(π/t)|x|

2

φ(x)dx.

34

Now multiply both sides by t

(α/2)

−1

dt and integrate from 0 to

∞.

Z

0

Z

R

n

e

−π|ξ|

2t

t

(α/2)

−1

ˆ

φ(ξ)dξ dt =

Z

0

Z

R

n

e

−π/t|x|

2

t

n)/2−1

φ(x)dx dt.

The change of order in the integral is permitted now and integrating,

we obtain,

Γ(α/2)

π

α/2

Z

R

n

ˆ

φ(ξ)

| ξ |

−α

dξ =

Γ((n

− α)/2)

π

(n

−α)/2

Z

R

n

| x |

α

n

φ(x) dx

i.e.,

<

| ξ |

−α

, ˆ

φ > =<

Γ((n

− α)/2)

(Γ(α/2)

π

α

n/2

| x |

α

n

, φ > .

This shows that

F

α

(x) =

Γ((n

− α)/2)

Γ(α/2)

π

α

n/2

| x |

α

n

as desired.

This analysis does not suffice to explain

−(4π

2

| ξ |

2

)

−1

as (in some

sense) the Fourier transform of (2π)

−1

log

| x | in the case n = 2. One

way to proceed is to define G

α

(ξ) = (2π

| ξ |)

−α

and to study the be-

haviour of G

α

and G

α

as α tends to 2. This analysis can be carried out

just as easily in n dimensions where the problem is to study G

α

and G

v

α

as α tends to n.

background image

4. Laplace Operator

31

Let R

α

= G

α

for 0 < α < n. (R

α

called the Riesz potential of order

α). By the preceding theorem, we have

R

α

(x) =

(Γ((n

− α)/2)

2

α

π

n/2

Γ(α/2)

| x |

α

n

.

35

Moreover, if f ǫS ,

(

−∆)

α/2

( f

R

α

) = f,

in the sense that (2π

| ξ |)

α

( f

R

α

)

ˆ

(ξ) = ˆf(ξ).

As α tends to n, the Γ-function in R

α

blows up. However, since

(2π

| ξ |)

α

δ

= 0, we see that (

−∆)

α/2

1 = 0; so we can replace R

α

by

R

α

c, c being a constant, still having

(

−∆)

α/2

( f

∗ (R

α

c)) = f.

If we choose c = c

α

appropriately, we can arrange that R

α

c

α

will

have a limit as α tends to n. In fact, let us take

c

α

=

Γ((n

− α)/2)

2

α

π

n/2

Γ(α/2)

and define R

α

= R

α

c

α

.

Then

R

α

(x) =

Γ((n

− α)/2)

2

α

n

n/2

Γ(a/2)

(

| x |

α

n

−1)

=

2Γ((n

− α)/2) + 1

2

α

π

n/2

Γ(α/2)

(

| x |)

α

n

− 1

n

− α

.

Letting α tend to n, we get in the limiting case

R

n

(x) =

−2

1

n

π

n/2

Γ(n/2)

log

| x | .

when n = 2,

R

2

(x)

1

log

| x | and ∆( f ∗ (−R

2

)) = f , so we recover our

fundamental solution for the case n = 2.

It remains to relate the function G

n

(ξ) = (2π

| ξ |)

n

to the Fourier

36

background image

32

2. Partial Differential Operators with Constant Coefficients

transform of the tempered distribution F

n

. One way to make G

n

into a

distribution is the following.

Define a functional F on S by

< F, φ >=

Z

|ξ|≤1

φ(ξ)

− φ(0)

(2π

| ξ |)

n

dξ +

Z

|ξ|>1

φ(ξ)

(2π

| ξ |)

n

dξ.

Note that the first integral converges, since, in view of the mean

value theorem,

| φ(ξ) − φ(0) |≤ c | ξ |. It is easy to see that F is

indeed a tempered distribution. Further, we observe that when φ(0) = 0,
< F, φ >=

R

φ(ξ)G

n

(ξ)dξ, i.e., F agrees with G

n

on R

n

\{0}.

Just as R

n

was obtained from R

n

by subtracting off an infinite con-

stant, F is obtained from G

n

by subtracting off an infinite multiple of the

δ function. This suggests that F is essentially the Fourier transform of
R

n

. In fact, one has the following

Exercise. Show that F

− (R

n

is a multiple of δ.

We now examine a few of the basic properties of harmonic functions,

i.e., functions satisfying u = 0. We shall need the following results from
advanced calculus.

Theorem 2.20 (DIVERGENCE THEOREM). Let Ω be a bounded open
set in R

n

with smooth boundary ∂Ω. Let v be the unit outward normal

vector on ∂Ω. Let F : R

n

→ R

n

be a C

1

function on ¯

, the closure of

. Then we have

Z

∂Ω

F. v dσ =

Z

div F dx =

Z

n

X

j=1

F

j

x

j

dx.

37

CONSEQUENCES OF THEOREM 2.20

(2.21)

When we take F = grad u, F.v = grad u.v =

u

v

, the normal

derivative of u, and div F = div grad u = ∆u. Therefore

Z

∂Ω

u

v

dσ =

Z

u dx.

background image

4. Laplace Operator

33

(2.22)

When we take F = u grad v

v grad u, div F = uv vu.

Therefore

Z

∂Ω

u

v
v

v

u

v

!

dσ =

Z

(uv

vu)dx.

This formula is known as Green’s formula.
For harmonic functions, we have the following mean value theorem.

Theorem 2.23. Let u be harmonic in B(x

0

, r). Then,for any ρ < r, we

have

u(x

0

) =

1

ω

n

ρ

n

−1

Z

|xx

0

|=ρ

u(x)dσ(x).

That is, u(x

0

) is the mean value of u on any sphere centred at x

0

Proof. Without loss of generality, assume that x

0

= 0. Since ∆u =

0 and ∆ is hypoelliptic, u is a C

function. Now formally,

u(0) =< δ, u >

=

Z

|x|<ρ

δ(x)u(x)dx

=

Z

|x|<ρ

K(x)u(x)dx

=

Z

|x|<ρ

(uK

Ku)dx

=

Z

|x|=ρ

u

K

v

K

u

v

!

dσ, by (2.22).



38

Of course, we are cheating here by applying Green’s formula to the

non-smooth function K. Nonetheless the result is correct, and we leave it
as an exercise for the reader to justify it rigorously (either approximate
K by C

functions as in the proof of Theorem 2.16 or apply Green’s

formula to the region ǫ <

| x |< ρ and let ǫ tend to 0.)

background image

34

2. Partial Differential Operators with Constant Coefficients

On the circle

| x |= ρ, K is a constant and so

R

|x|=ρ

K

u

v

dσ =

Const

R

|x|=ρ

u

v

dσ = Const .

R

|x|≤ρ

u dσ = 0.

Further on

| x |= ρ.

v

=

r

, the radial derivative. Therefore

K

v

= 1/(ω

n

ρ

n

−1

) on

| x |= ρ.

Thus we have

u(0) =

1

ω

n

ρ

n

−1

Z

|x|=ρ

u(x)dσ

as desired.

As a corollary to the mean value theorem, we can derive the maxi-

mum principle for harmonic functions.

Theorem 2.24. Suppose is a connected open set in R

n

. Let u be a real-

39

valued function harmonic in . If A = sup

u(x), then either u(x) < A for

all xǫΩ or u(x)

A on .

Proof. Suppose that u(x

0

) = A for some x

0

ǫΩ. From the mean value

theorem

u(x

0

) =

1

ω

n

ρ

n

−1

Z

|xx

0

|=ρ

u(x)dσ(x)

where ρ is small enough so that

{x :| x x

o

|≤ ρ} ⊂ Ω. By our as-

sumption, the integral does not exceed A. If u(x

1

) < A for some point

of B(x

0

, ρ) then u(x) < A in some neighbourhood of x

1

by continuity.

Taking r =

| x

1

x

0

|< ρ, we have

u(x

0

) =

1

ω

n

r

n

−1

Z

|xx

o

|=r

u(x)dσ(x) < A,

a contradiction. Therefore, if we set Ω

1

=

{xǫΩ : u(x) = A}, then Ω

1

is

an open subset of Ω and Ω

1

,

φ. Further Ω

2

=

{xǫΩ : u(x) < A} is also

on open subset of Ω and Ω

1

∪ Ω

2

= Ω. The connectedness of Ω and the

non-emptiness of Ω

1

force Ω

2

to be empty. Thus u(x)

A on Ω.



background image

4. Laplace Operator

35

Corollary 2.25. If is a bounded open set in R

n

, uǫC( ¯

Ω) and u =

0 in , then

sup

| u(x) |= sup

∂Ω

| u(x) | .

Proof. Since the function u is continuous on the compact set ¯

Ω,

| u(x) |

attains its supremum on ¯

Ω at some point x

0

. By multiplying u by a

constant, we may assume that u(x

0

) =

| u(x

0

)

|.



If x

0

ǫ∂Ω, the corollary is proved. Otherwise x

0

ǫΩ and by the previ-

40

ous theorem, on the connected component of containing x

0

, Re u(x) =

u(x

0

) and hence Im u(x) = 0. Since u is continuous on ¯

Ω, u(x) = u(x

0

)

for all xǫ∂Ω. Thus sup

| u(x) |= sup

∂Ω

| u(x) |.

Corollary 2.26. If u and v are in C( ¯

Ω), Ω is bounded u = ∆v = 0 on

and u = v on ∂Ω, then u = v everywhere.

Proof. Apply the previous corollary to u

v.



The following boundary value problem for Laplace’s equation,

known as the Dirichlet problem, is of fundamental importance: given
a function f on ∂Ω, find a function u such that ∆u = 0 on Ω and u =

f on ∂Ω.

When Ω is bounded and f ǫC(∂Ω) the uniqueness of the solution,

if it exists at all, is assured by corollary 2.26. The problem of proving
the existence of solution of the Dirichlet problem is highly nontrivial.
We shall solve a special case, namely, when Ω is a half-space and then
indicate how similar ideas can be applied for other regions.

First of all a word about notation: we will replace R

n

by R

n+1

with

coordinates (x

1

, x

2

, . . . , x

n

, t). Our Laplacian in R

n+1

will be denoted by

2

t

+ ∆ where ∂

t

=

t

and ∆ =

n

P

j=1

2

x

j

. Now our Dirichlet Problem is

the following :

given a function f on R

n

, find a function u such that

(∂

2
t

+ ∆) u(x, t) = 0, xǫR

n

, t > 0

(2.27)

u(x, o) = f (x), xǫR

n

.

background image

36

2. Partial Differential Operators with Constant Coefficients

41

Since the half space is unbounded, the uniqueness argument given

above does not apply. Indeed, without some assumption on the be-
haviour of u at

∞, uniqueness does not hold. For example, if u(x, t)

is a solution, so is u(x, t) + t. However, we have

Theorem 2.28. Let u be a continuous function on R

n

× [0, ∞) satisfying

i) (∆ + ∂

2

t

)u = 0 on R

n

× (0, ∞),

ii) u(x, o) = 0 for xǫR

n

, and

iii) u vanishes at

.

Then U

≡ 0 on R

n

× [0, ∞).

Proof. Applying the maximum principle for u on B(0, R)

× (0, T ), we

see that maximum of

| u | on B(0, R) × (0, T ) tends to zero as R, T tend

∞. Hence u ≡ 0.



Remark 2.29. Hypothesis (iii) in Theorem 2.28 can be replaced by (iii):
u is bounded on R

n

× [0, ∞). But this requires a deeper argument (See

Folland [1]).

To solve the Dirichelt problem, we apply the Fourier transform in

the variable x. We denote by ˜u(ξ, t) the Fourier transform

˜u(ξ, t) =

Z

e

−2πix

u(x, t)dx.

If we take Fourier transform of (2.27), we obtain

42

(∂

2
t

− 4π

2

| ξ |

2

u(ξ, t) = o on R

n

× (0, ∞)

˜u(ξ, o) = ˆf(ξ) on R

n

.

The general solution of the ordinary differential equation (ODE)

(∂

2

t

− 4π

2

| ξ |

2

) ˜u(ξ, t) = 0 is given by ˜u(ξ, t) = a(ξ)e

−2π|ξ|t

+ b(ξ)e

|ξ|t

and the initial condition is a(ξ) + b(ξ) = ˆf(ξ).

This formula for ˜u will define a tempered distribution in ξ, provided

that as

| ξ | tends to ∞, | a(ξ) | grows at most polynomially and | b(ξ) |

background image

4. Laplace Operator

37

decreases faster than exponentially. We remove the non- uniqueness by
requiring that u should satisfy good estimates in terms of f which are
uniform in t-for example, that u should be bounded if f is bounded. This
clearly forces b(ξ) = 0, since e

|ξ|t

blows up as t tends to

∞. Thus we

take ˜u(ξ, t) = ˆf(ξ)e

−2π|ξ|

t. Taking the inverse Fourier transform in the

variable ξ, u(x, t) = ( f

P

t

)(x) where P

t

(x) = [e

t(

|.|)

]

(x) is the Poisson

Kernel for R

n

× [0, ∞).

We now compute P

t

explicitly. When n = 1, this is easy:

P

t

(x) =

0

Z

−∞

e

ixξ

e

2πξt

dξ +

Z

0

e

ixξ

e

−2πξt

dξ

=

1

[(t + ix)

−1

+ (t

ix)

−1

] =

t

π

(x

2

+ t

2

)

−1

.

To compute P

t

(x) for the case n > 1, as in the proof of Theorem

2.19, the idea is to express e

−2π|ξ|t

as a weighted average of Gaussians.

Here’s how!

43

Lemma 2.30. When β > 0, e

−β

=

R

0

e

s

e

−β

2

/4s

s)

ds.

Proof. First, observe that

e

−β

=

1

π

Z

−∞

e

iβτ

τ

2

+ 1

dτ.



This can be proved by using the residue theorem or by applying the

Fourier inversion theorem to P

1

(x) =

1

π

(1 + x

2

)

−1

on R

1

. We also have

1

τ

2

+ 1

=

Z

0

e

−(τ

2

+1)s

ds.

background image

38

2. Partial Differential Operators with Constant Coefficients

Therefore,

e

−β

=

1

π

Z

−∞

Z

0

e

i2πτ

e

−(τ

2

+1)s

ds dτ.

Putting τ = 2πσ and changing the order of integration,

e

=

1

π

Z

0

Z

−∞

e

i2πσβ

e

s

e

−4π

2

σ

2

s

dσ ds

i.e.,

e

−β

= 2

Z

0

e

−2

(4πs)

−1/2

e

β

2

/4s

ds

=

Z

0

e

s

e

β

2

/4s

s)

ds

which is the required expression.

Returning to the computation of P

t

(x), we have

P

t

(x) =

Z

e

ix

e

−2π|ξ|t

dξ.

Taking β = 2π

| ξ | t in the lemma, we get

44

P

t

(x) =

Z

Z

0

e

ix

e

s

πs

e

(4π

2

|ξ|

2

t

2

)/4s

dsdξ

=

Z

0

e

s

πs

Z

e

ixξ

e

−(π

2

|ξ|

2

t

2

)/S

dξ.ds

=

Z

0

e

s

πs

πt

2

s

!

n/2e

−(s|x|

2

)/t

2

ds

= π

−(n+1)/2

t

n

Z

0

e

s(1+(|x|

2

/t

2

))

s

(n

−1)/2

ds

background image

4. Laplace Operator

39

=

π

−(n+1)/2

t

n

Γ(n + 1)/2



1 +

|x|

2

t

2



(n + 1/2)

.

Thus we have

(2.31)

P

t

(x) = Γ((n + 1)/2)π

(n+1)/2

t/(t

2

+

| x |

2

)

(n+1)/2

In particular, we see that P

t

ǫL

1

L

, so that f

P

t

is well defined if,

for example, f ǫL

p

, 1

p ≤ ∞. That (∆ + ∂

2

t

)P

t

= 0 follows by taking

the Fourier transform, and hence

(∆ + ∂

2
t

)( f

P

t

) = f

∗ (∆ + ∂

2
t

)P

t

= 0 for any f.

Further note that

P

t

(x) = t

n

P

1

(x/t) and

Z

P

1

(x)dx = ˆ

P

1

(0) = 1.

Therefore, by Theorem 1.6, when f ǫL

p

, 1

p < ∞, f P

t

tends to

f in the L

p

norm and when f is bounded and continuous f

P

t

tends to

f uniformly on compact sets as t tends to 0.

If we take f continuous and bounded, then for u(x, t) = P

t

f , we

have lim

t

→0

u(x, t) = u(x, 0) = f (x). Thus the function u(x, t) = ( f

P

t

)(x)

45

is a solution of the Dirichlet Problem for the half space.

Remark 2.32. The Poisson Kernel P

t

is closely related to the fundamen-

tal solution K

n+1

of the Laplacian in R

n+1

. Indeed

K

n+1

(x, t) =

1

(1

n

n+1

(

|x|

2

+ t

2

)

−((n−1)/2)

and hence

P

t

(x) = 2∂

t

K

n+1

(x, t).

Exercise. Check the above equation using the Fourier transform. (Start
with (

2

t

+ ∆)K

n+1

= δ(x)δ(t) and take the Fourier transform in the vari-

able x to obtain

(∂

2
t

− 4π

2

|ξ|

2

) ˜

K

n+1

= δ(t).

background image

40

2. Partial Differential Operators with Constant Coefficients

Solve this equation to obtain

˜

K

n+1

(ξ, t) =

e

|ξ||t|

|ξ|

.

Then

t

˜

K

n+1

=

1

2

e

−2π|ξ|t

, t > 0, so that

t

˜

K

n+1

=

1

2

ˆ

P

t

).

Our formula (2.31) for P

t

makes sense even when t < 0 and we have

P

t

(x) =

P

t

(x), so that lim

t

o±

f

P

t

=

± f . We further observe that

f

(x)

P

t

= f (x)δ(t)

(x,t)

2

t

K

n+1

(x, t)

where

(

x) and

(x,t)

mean convolution on R

n

and R

n+1

respectively.

Therefore,

46

f

(x)

P

t

= f (x

(t)

(x,t)

2K

n+1

(x, t)

i.e.,

u(x, t) = 2 f (x

(t)

(

x, t)K

n+1

(x, t).

Form this, we

2
t

+ ∆)u(x, t) = 2 f (x

(t)

∗ δ(x)δ(t) = 2 f (x

(t).

Exercise. Show directly that if

i) (∆ + ∂

2

t

)u = 0 on R

n+1

\R

n

× {0}

ii) u(x,

t) = −u(x, t)

iii) lim

t

→0±

u(x, t) =

± f (x), then u is a distribution solution of (∂

2

t

+ ∆)u =

2 f (x

(t).

[To avoid technicalities, assume f is sufficiently smooth so that

lim

t

→0±

u

t

exists; e.g. f ǫC

2

is sufficient ].

We now indicate, without giving any proof, how these ideas can be

used to solve the Dirichlet Problem in a bounded open set Ω in R

2

. We

assume that the boundary ∂Ω of Ω is of class C

2

.

background image

5. The Heat Operator

41

We know that the solution of the Dirichlet Problem in the case when

Ω = R

n

−1

× (0, ∞) is given by u = f ∗ 2

K

x

n

(the convolution is in R

n

−1

).

we note that

K

x

n

is the inward normal derivative of K. For bounded Ω

with C

2

boundary let us consider

u(x) = 2

Z

∂Ω

f (y)

v

y

K(x

y)dσ(y).

Here σ(y) is the surface measure on the boundary, ν is the unit out-

47

ward normal on ∂Ω, and

∂φ

∂ν

y

(x, y) = Σν

j

∂φ

y

j

. The minus sign in K(x

y)

compensates the switch from inward to outward normal.

Since ∆K(x

y) = δ(x y) we see that ∆u = 0 in Ω. What about the

behaviour of u on ∂Ω? It turns out that if u is defined as above, then

uǫC ¯

(Ω) and u

|

∂Ω

= f + T f

where T is a compact integral operator on L

2

(∂Ω) or C(∂Ω). Hence if

we can find φ in C(∂Ω) such that φ + T φ = f , and we set

u(x) = 2

Z

∂Ω

φ(y)

∂ν

y

K(x

y)dσ(y)

then u satisfies ∆u = 0 in Ω and further u

|

∂Ω

= φ + T φ = f . Hence

the Dirichlet problem is reduced to solving the equation φ + T φ = f ,
and for this purpose, the classical Fredholm - Riesz theory is available.
The upshot is that a solution to the Dirichlet problem always exists (pro-
vided, as we have assumed, that ∂Ω is of class C

2

). See Folland [1] for

a detailed treatment.

5 The Heat Operator

The Heat operator is given by ∂

t

− ∆. We want to find a distribution

48

K such that (∂

t

− ∆)K = δ(x)δ(t). Taking the Fourier transform in both

variables we have

(2.33)

ˆ

K(ξ, τ) = (2πi + 4π

2

|ξ|

2

)

−1

.

background image

42

2. Partial Differential Operators with Constant Coefficients

Exercise . Show that ˆ

K is locally integrable near the origin, and hence

defines a tempered distribution.

ˆ

K is not globally integrable, however; so computing its inverse Fou-

rier transform directly is a bit tricky. Instead, if we take the Fourier
transform in the variable x only,

(∂

t

+ 4π

2

|ξ|

2

) ˜

K(ξ, t) = δ(t).

we can solve this ODE get

˜

K(ξ, t) =


a(ξ)e

−4π

2

|ξ|

2

t

,

t > 0

b(ξ)e

−4π

2

|ξ|

2

t

,

t < 0

with a(ξ)

b(ξ) = 1.

As in the previous section, there is some latitude in the choice of a

and b, but the most natural choice is the one which makes ˜

K tempered

in t as well as ξ, namely a = 1, b = 0. So,

˜

K(ξ, t) =


e

−4π

2

|ξ|

2

t

, t > 0

0 otherwise .

Taking the inverse Fourier transform,

49

K(x, t) =


(4πt)

n/2

e

(

|x|

2

/4t)

,

t > 0

0,

t > 0.

This is a fundamental solution of the heat operator.

Remark 2.34. If we take the Fourier transform of ˜

K(ξ, t) in the t vari-

able, we obtain

˜

K(ξ, τ) =

Z

−∞

˜

K(ξ, t)e

2π i t τ

dt

=

Z

0

e

t(4π

2

|ξ|

2

i τ)

dt

= (4π

2

|ξ|

2

+ 2πiτ)

−1

,

thus recovering formula (2.33).

background image

5. The Heat Operator

43

Exercise . Show that this really works, i.e., the iterated Fourier trans-
form of K first in x, then in t, is the distribution Fourier transform of K
in all variables.

OBSERVATIONS: From the formula for K, we get

(1) K(x, t) vanishes to infinite order as t tends to 0 when x , 0 and

hence is C

on R

n+1

\{(0, 0)}. Therefore by Theorem 2.11, ∂

t

− ∆ is

hypoelliptic.

(2) K(x, t) = t

n/2

K



x

t

, 1



. Thus if we set K(x, 1) = φ(x), then K(x, ǫ

2

)

= ǫ

n

φ(x/ǫ) = φ

ǫ

(x); so, as t tends to 0, K(x, t) tends to δ, Theorem

1.6. From this, we infer

50

(3) If f ǫL

p

and if we set u(x, t) = f

(x)

K(x, t), then

(∂

− ∆)u = 0 for t > 0

u(x, 0) = f (x).

i.e., as t tends to 0,

||u(., t) − f ||

P

converges to 0. Thus we have

solved the initial value problem for the homogeneous heat equation,
when f ǫL

p

. Actually, since K(x, t) decreases so rapidly as

|x| tends

to

∞, this works for much wider classes of f

s. The convolution

f

(x)

K(., t) make sense, for example, if

| f (x)| < Ce

|x|

2

−ǫ

. If f is

also continuous, it is not hard to see that f

(x)

K(., t) converges to f

uniformly on compact sets as t tends to 0.

It is now a simple matter to solve the inhomogeneous initial value

problem:

(∂

t

− ∆)u = F on R

n

× (0, ∞),

u(x, 0) = f (x) on R

n

.

If we take u

1

= F

(x,t)

K and

u

2

= ( f

u

1

(., 0))

(x)

K,

background image

44

2. Partial Differential Operators with Constant Coefficients

then we see that

(∂

t

− ∆)u

1

= F on R

n

× (0, ∞), (∂

t

− ∆)u

2

= 0

on R

n

× (0, ∞), and (u

1

+ u

2

)(x, 0) = f (x). Thus u = u

1

+ u

2

solves the

51

problem.

As another application of the fundamental solution K, we can derive

the Weierstrass approximation theorem.

Theorem 2.35. (WEIERSTRASS) If f is continuous with compact
support on R

n

, then, for any compact Ω

⊂ R

n

, there exists a sequence

(P

m

) of polynomials such that P

m

converges to f uniformly on Ω

Proof. Let u(x, t) = ( f

(x)

K(., t))(x). Then u(x, t) converges to f (x)

uniformly as t tends to 0. Moreover, for any t > 0,

u(x, t) =

Z

f (y)(4πt)

n/2

e

n

P

j=1

(x

j

y

j

)

2

/4t

dy

is an entire holomorphic function of x

∈ ƒ

n

. So u(., t) can be uniformly

approximated on any compact set by partial sums of its Taylor series.



6 The Wave Operator

If we take the Fourier transform of the equation

(∂

2
t

− ∆)K = δ(x)δ(t)

in both variables x and t, we have, formally

(2.36)

ˆ

K(ξ, τ) = (4π

2

|ξ|

2

− 4π

2

τ

2

)

−1

.

This function ˆ

K is not locally integrable, so it is not clear how to

52

interpret it as a distribution. Again, it is better to take the Fourier trans-
form in the x variable obtaining

(∂

2
t

+ 4π

2

|ξ|

2

) ˜

K(ξ, t) = δ(t).

background image

6. The Wave Operator

45

Solving this ODE,

˜

K(ξ, t) = a

±

(ξ)e

i

|ξ|t

+ b

±

(ξ)e

−2πi|ξ|t

for t/

|t| = ±1

and the coefficients a

±

, b

±

must be determined so that

˜

K(ξ, 0+) = ˜

K(ξ, 0

−), ∂

t

˜

K(ξ, 0+)

− ∂

t

K(ξ, 0 =) = 1.

This gives two equations in four unknowns. In contrast to the situa-

tion with the Dirichlet problem and the heat operator, there is no way to
narrow down the choices further by imposing growth restrictions on K
as

|t| tends to ∞. Rather, it is a characteristic feature of the wave opera-

tor that one can adapt the choice of fundamental solution to the problem
at hand. The two which we shall use, called K

+

and K

,

are the ones sup-

ported in the half - space t

≥ 0 and t ≤ 0. K

+

and K

are thus determined

by the requirements a

= b

= 0 and a

+

= b

+

= 0 respectively, from

which one easily sees that

˜

K

+

(ξ, t) = H(t)

sin 2π

|ξ|t

|ξ|

˜

K

(ξ, t) =

H(−t)

sin 2π

|ξ|t

|ξ|

= ˜

K

+

(ξ,

t)

where H is the Heaviside function, i.e., the characteristic function of

53

[0,

∞). Let us compute the Fourier transforms of ˜

K

+

and ˜

K

in t to see

how to make sense out of (2.36). ˜

K

+

and ˜

K

are not integrable on t, but

it is easy to approximate them in the distribution topology by integrable
functions whose Fourier transforms we can calculate.Indeed, if we set

˜

K

ǫ

+

(ξ, t) = e

−2πǫt

H(t)

sin 2π

|ξ|t

|ξ|

, ǫ > 0,

then ˜

K

ǫ

+

is an integrable function and ˜

K

ǫ

+

converges to ˜

K

+

is S

as ǫ tends

to 0. Therefore ˜

K

+

= lim

ǫ

→0

˜

K

ǫ

+

where

ˆ

K

+

(ξ, τ) =

Z

0

e

2πǫt

−2π it τ

sin 2π

|ξ|t

|ξ|

dt

background image

46

2. Partial Differential Operators with Constant Coefficients

=

Z

0

e

2πǫt

−2π it τ

(e

πi

|ξ|t

e

πi

|ξ|t

)

i

|ξ|

dt

= (4π

2

)

−1

(

|ξ|

2

− (τ − iǫ)

2

)

−1

.

Exercise. Prove that ˆ

K

= lim

ǫ

→0

ˆ

K

ǫ

in S

where

ˆ

K

ǫ

(ξ, τ) = (4π

2

)

−1

(

|ξ|

2

− (τ + iǫ)

2

)

−1

.

Thus we have two distinct ways of making the function [4π

2

(

|ξ|

2

τ

2

)]

−1

into a tempered distribution. The difference ˜

K

+

− ˜

K

is of course

a distribution supported on the cone

|ξ| = |τ|.

We now propose to use the fundamental solutions K

+

and K

to

54

solve the Initial Value Problem or Cauchy Problem, for the operator:

(2.37)

(∂

2
t

− ∆)u = f on R

n+1

u(x, o) = u

o

(x),

t

u(x, 0) = u

1

(x),

where u

0

, u

1

, f are given functions.

For the moment, we assume that u

0

, u

1

ǫS and f ǫC

(R

n

x

)) ( That is,

t

f (., t) is a C

function with values in S (R

n

)).

Taking the Fourier transform in the variable x in (2.37),

(∂

2
t

+ 4π

2

|ξ|

2

) ˜u(ξ, t) = ˜f(ξ, t)

˜u(ξ, 0) = ˜u

0

(ξ)

t

˜u(ξ, 0) = ˜u

0

(ξ).

When f = 0, the general solution of the ODE is

˜u(ξ, t) = A(ξ) sin 2π

|ξ|t + B(ξ) cos 2π|ξ|t,

A(ξ) =

ˆu

1

(ξ)

|ξ|

and B(ξ) = ˜u

0

(ξ).

when u

0

= u

1

= 0, the solution is

˜u(ξ, t) =

t

Z

0

˜

f (ξ, s)

sin(2π

|ξ|(t s))

|ξ|

ds

background image

6. The Wave Operator

47

(This may be derived by variation of parameters; in any case it is

55

easy to check that this u is in fact the solution). Therefore, the solution
for the general case is given by

˜u(ξ, t) = ˆu

0

(ξ) cos 2π

|ξ|t +

ˆu

1

(ξ)

|ξ|

sin 2π

|ξ|t+

+

t

Z

0

˜

f (ξ, s)

sin(2π

|ξ|(t s))

|ξ|

ds.

For t > 0, we can rewrite this as

˜u(ξ, t) = ˆu

1

(ξ) ˜

K

+

(ξ, t) + ˆu

0

(ξ)∂

t

˜

K

+

(ξ, t) +

t

Z

0

˜

f (ξ, s) ˜

f (ξ, s)K

+

(ξ, t

s)ds.

Since K

+

(ξ, t

s) = 0 for t < s,

˜u(ξ, t) = ˆu

1

(ξ) ˜

K

+

(ξ, t) + ˆu

0

(ξ)∂

t

˜

K

+

(ξ, t) +

t

Z

0

˜

f (ξ, s) ˜

f (ξ, s)K

+

(ξ, t

s)ds.

Take the inverse Fourier transform:

u(x, t) = (u

1

(x)

K

+

) + u

0

(x)

t

K

+

) + (H(t) f

(x,t)

K

+

).

Likewise, for t < 0

u(x, t) =

−(u

1

(x)

K

)

u

0

(x)

t

K

) + ((H(

t) f ) ∗

(x,t)

K

).

So, for arbitrary t, our solution u can be expressed as

u(x, t) = (u

1

(x)

(K

+

K

)

+ u

0

(x)

t

(K

+

K

)

+ (H(t) f

(x,t)

K

+

) + (H(

t) f

(x,t)

K

)

).

(2.38)

So far, we have avoided the question of computing K

+

and K

ex-

56

plicitly. Indeed, since ˜

K

+

and ˜

K

are not L

1

functions, it is not an easy

matter to calculate their inverse Fourier transform,

background image

48

2. Partial Differential Operators with Constant Coefficients

However, for the case n = 1, we can find K

+

and K

by solving the

wave wave equation directly. We have

2
t

− ∆ = ∂

2
t

− ∂

2

x

.

If we make the change of variables ξ = x + t, n = x

t, then the wave

operator becomes ∂

2

t

− ∂

2

x

=

−4

2

∂ξ∂η

. The general solution of

2

u

∂ξ∂η

= 0

is given by

u(ξ, n) = f (ξ) + g(η)

where f and g are arbitrary functions. Therefore

u(x, t) = f (x + t) + g(x

t).

To solve (∂

2

t

− ∂

2

x

)u = 0, u(x, 0) = u

0

(x), ∂

t

u(x, 0) = u

1

(x) we must

have

u

0

(x) = f (x) + g(x)

u

1

(x) = f

(x)

g

(x).

From these equations, we have

f

(x) =

1

2

(u

0

(x) + u

1

(x))

g

(x) =

1

2

(u

0

(x)

u

1

(x)).

Thus u(x, t) is given by

57

u(x, t) =

1

2

(u

0

(x + t) + u

0

(x

t)) +

1

2

x+t

Z

x

t

u

1

(s)ds.

Comparing with the previous formula (2.38), we find that

K

+

(x, t) =

1

2

H(t

− |x|), K

(x, t) =

1

2

H(

t − |x|).

Exercise. Compute ˜

K

±

directly from these formulas.

background image

6. The Wave Operator

49

It turns out that, for n = 2,

K

±

(x, t) =

1

t

2

x

2

H(

±t − |x|),

and, for n = 3,

K

±

(x, t) =

±1

t

δ(

±t − |x|).

For n = 1, 2, K

+

, K

are functions; for n = 3, K

+

, K

are not func-

tions but they are measures. For n > 3, K

+

, K

are neither functions nor

measures; they are more singular distributions. The exact formula for
K

±

is rather messy and we shall not write it out; it may be read off from

Theorem 5.13 and 5.14 of Folland [1] in view of our formula (2.38). The
most important qualitative feature of K

±

, however, is that it is always

supported in the light cone

{(x, t) : ±t > |x|} and we shall now prove this

as a consequence of the following result.

Theorem 2.39. Suppose u is a C

2

function on

{(x, t) : t o} such that

58

(∂

2

t

−△)u = 0 for t > 0 and u = ∂

t

u = 0 on the set B

0

=

{(x, 0) : |xx

0

| ≤

t

0

}. Them u vanishes on Ω = {(x, t) : 0 ≤ t t

0

,

|x x

0

| ≤ t

0

t}.

Proof. Assume u is real valued; otherwise, we can consider the real and
imaginary parts separately. Let

B

t

=

{x : |x x

0

| ≤ t

0

t} and E(t) =

1

2

Z

B

t

| grad

x,t

u

|

2

dx

i.e.,

E(t) =

1

2

Z

B

t

(∂

t

u)

2

+

n

X

1

u

x

j

!

2

dx.

2

Then

dE

dt

=

1

2

Z

B

t

2

u

t

.

2

u

2

t

+

n

X

1

u

x

j

2

u

x

j

t

dx

1

2

Z

B

t

u

t

!

2

+

n

X

1

u

x

j

!

2

dσ(x).

background image

50

2. Partial Differential Operators with Constant Coefficients

(The second term comes from the change in the region B

t

. If this is

not clear, write the derivative as a limit of difference quotients and work
it our).

Since

X

u

x

j

2

u

x

j

t

+

X

2

u

x

2

j

u

t

= div

u

t

grad

x

u

!

,

applying the divergence theorem and using (∂

2

t

− △)u = 0, we obtain

dE

dt

=

Z

B

t

" ∂u

t

u

∂ν

1

2

| grad

x,t

u

|

2

#

dσ

where ν is the unit normal to B

t

in R

n

. Now

59

|

u

t

u

∂ν

| <

1

2

"

|

u

t

|

2

+

|

u

∂ν

|

2

#

<

1

2

"

|

u

t

|

2

+

| grad

x

u

|

2

=

1

2

| grad

x,t

u

|

2

#

.

Thus we see that the integrand is non-positive, and hence

dE

dt

≤ 0.

Also E(0) = 0 since u = ∂

t

u = 0 on B

0

, so E(t)

≤ 0. But E(t) ≥ 0 by

definition, so E(t) = 0. This implies that grad

x,t

u = 0 on Ω =

S

t

t

0

B

t

and

since u = 0 on B

0

, we conclude that u = 0 on Ω.

Corollary 2.40. Suppose uǫC

2

on R

n

× [0, ∞), (∂

2

t

− △)u = 0 for t >

0, u(x, 0) = u

0

(x), ∂

t

u(x, 0) = u

1

(x). If

0

= (supp u

0

)

∪ (supp u

1

), then

supp u

⊂ Ω = {(x, t) : d(x, Ω

0

)

t}.

(The set Ω is the union of the forward light cones with vertices in

0

).

Proof. Suppose (x, t

0

) < Ω. Then for some ǫ > 0, the set

B

0

=

{x : d(x, x

0

)

t

0

+ ǫ

} is disjoint from Ω

0

.

background image

6. The Wave Operator

51

Therefore, by Theorem 2.39, u = 0 on the cone

{(x, t) : |x x

0

| ≤ t

0

+ ǫ

t, 0 ≤ t t

0

+ ǫ

}.

In particular, u = 0 on a neighbourhood of (x

0

, t

0

), i.e., (x

0

, t

0

) is not in

the support of u.



Corollary 2.41.

SuppK

+

⊂ {(x, t) : t ≥ |x|}

SuppK

+

⊂ {(x, t) : −t ≥ |x|}.

60

Proof. Pick aφǫC

α
0

(B(0, 1)) such that

R

φ = 1. Put

φ

ǫ

(x) = ǫ

n

φ(ǫ

−1

x). Let u

ǫ

(x, t) = φ

ǫ

(x)

K

+

, t > 0.



Then (∂

2

t

− △)u

ǫ

= 0, u(x, 0) = 0, ∂

t

u(x, 0) = φ

ǫ

(x) and u

ǫ

is C

α

.

By the previous corollary

supp u

ǫ

⊂ {(x, t) : |x| ≤ t + ǫ}.

Now u

ǫ

converges to K

+

in S

, as ǫ tends to 0. Therefore,

supp K

+

⊂ {(x, t) : |x| ≤ t}.

The result for K

follows then since K

(x, t) = K

+

(x,

t).

Remarks 2.42.

(i) One could also deduce the above result from our

formulas for ˜

K

±

by using the Paley-Wiener theorem.

(ii) Actually for n = 3, 5, 7, . . ., it turns out that supp K

±

=

{(x, t) :

|x| = ±t}. This is known as the Huygens principle. See Folland
[1].

background image

52

2. Partial Differential Operators with Constant Coefficients

(iii) The distributions K

±

are smooth functions of t (except at t = 0)

with values in E

(R

n

). Therefore, we now see that our formula

(2.38) for the solution of the Cauchy problem makes sense even
when u

0

, u

1

, ǫD

(R

n

) and f ǫC(R

t

, D

(R

n

x

)), and it is easily checked

61

that u thus defined still solves the Cauchy problem in the sense
of distributions. Corollary 2.40 remains valid in this more general
setting, as can seen by an approximation argument as in the proof
of Corollary 2.41.

EXERCISE If (∂

2

t

−△)u = f on R

n+1

, u(x, 0) = u

0

(x), and ∂

t

u(x, 0)

= u

1

(x), figure out how supp u is related to supp f, supp u

0

and supp u

1

.

background image

Chapter 3

L

2

Sobolev Spaces

THERE ARE MANY ways of measuring smoothness properties of func-

62

tions in terms of various norms. Often it is convenient to use L

2

norms,

since L

2

interacts nicely with the Fourier transforms. In this chapter, we

set up a precise theory of L

2

differentiability and use it to prove Horman-

der’s theorem on the hypoellipticity of constant coefficient differential
operators.

1 General Theory of L

2

Sobolev Spaces

Definition 3.1. For a non-negative integer k, the Sobolev space H

k

is

defined to be the space of all tempered distributions all of whose deriva-
tives of order less than or equal to k are in L

2

.

Thus

H

k

=

{ f ǫS

: D

α

f ǫL

2

(R

n

) for 0

≤ |α| ≤ k}.

From the definition, we note that f ǫH

k

if and only if ξ

α

ˆ

f (ξ)ǫL

2

for

0

≤ |α| ≤ k.

Proposition 3.2. f ǫH

k

if and only if (1 +

|ξ|

2

)

k/2

ˆ

f ǫL

2

.

Proof. First assume that (1 +

|ξ|

2

)

k/2

ˆ

f ǫL

2

. Since, for

|ξ| ≥ 1,

α

| ≤ |ξ|

k

≤ (1 + |ξ|

2

)

k/2

for all

|α| ≤ k

53

background image

54

3. L

2

Sobolev Spaces

and for

|ξ| < 1,

63

α

| ≤ 1 ≤ (1 + |ξ|

2

)

k/2

for all

|α| ≤ k,

we have

α

ˆ

f (ξ)

| ≤ (1 + |ξ|

2

)

k/2

ˆ

f (ξ)

and hence

ξ

α

ˆ

f ǫL

2

for

|α| ≤ k which implies f ǫH

k

.

2



Conversely, assume that f ǫH

k

. Since

|ξ|

k

and

n

P

j=1

j

|

k

are homoge-

neous of degree k and nonvanishing for ξ , 0, we have

(1 +

|ξ|

2

)

k/2

c

0

(1 +

|ξ|

k

)

c

0

1 + c

n

X

j=1

j

|

k

so that

||(1 + |ξ|

2

)

k/2

ˆ

f

||

2

c

0

|| f ||

2

+ c

0

c

n

X

j=1

||∂

k

j

f

||

2

which shows that (1 +

|ξ|

2

)

k/2

ˆ

f ǫL

2

.

The characterisation of H

k

given in the above proposition immedi-

ately suggests a generalisation to non-integral values of k which turns
out to be very useful.

Definition 3.3. For sǫR, we define the operator Λ

s

: S

S by

s

f )ˆ(ξ) = (1 +

|ξ|

2

)

s/2

ˆ

f (ξ).

In other words, Λ

s

=



I

2



s/2

. Clearly Λ

s

maps S continuously

64

onto itself. We can therefore extend Λ

s

continuously from S

onto itself.

The Sobolev space of order s is defined by

H

s

=

{ f ǫS

: Λ

s

f ǫL

2

}.

We equip H

s

with the norm

|| f ||

(s)

=

||Λ

s

f

||

2

. If s is a positive integer,

the proof of Proposition 3.2 shows that

|| ||

(s)

is equivalent to the norm

|| f || = Σ

0

≤|α|≤s

||D

α

f

||

2

.

PROPERTIES OF H

s

background image

1. General Theory of L

2

Sobolev Spaces

55

(i) H

s

is a Hilbert space with the scalar product defined by (u, v)

(s)

=

s

u, Λ

s

v). Here the scalar product on the right is that in L

2

. The

Fourier transform is a unitary isomorphism between H

s

and the

space of functions which are square integrable with respect to the
measure (1 +

|ξ|

2

)

s

dξ.

(ii) For every sǫR, S is dense in H

s

.

(iii) If s > t, H

s

H

t

with continuous imbedding. In fact, for uǫH

s

,

||u||

(t)

≤ ||u||

(s)

. In particular, H

s

L

2

for s > 0.

(iv) D

α

is a bounded operator from H

s

into H

s

−|α|

sǫR.

65

(v) If f ǫH

s

, then f as a linear functional on S extends continuously

to H

s

and

|| f ||

(

s)

is the norm of f in (H

s

)

. So we can identify

(H

s

)

with H

(

s)

. For f ǫH

s

, gǫH

s

the pairing is given by

< f, g >=< Λ

s

f, Λ

s

g >=

Z

ˆ

f ˆg.

(If s = 0, this identification of H

0

= L

2

with its dual is the complex

conjugate of the usual one).

(vi) The norm

||.||

(s)

is translation invariant. Indeed, if g(x) = f (x

x

0

),

then ˆg(ξ) = e

ix

0

ˆ

f (ξ) and hence

|| f ||

(s)

=

||g||

(s)

.

Proposition 3.4. For s > n/2, we have δǫH

s

.

Proof. (1 +

|ξ|

2

)

s/2

ˆδ = (1 + |ξ|

2

)

s/2

ǫL

2

, whenever s > n/2. This is a

consequence of the following observation :



R

(1 +

|ξ|

2

)

s

dξ

∼ 1 +

R

1

r

−2s

r

n

−1

dr <

∞ if and only if s > n/2.

As an immediate consequence of this proposition, we have

Corollary 3.5. For s > n/2 +

|α|, D

α

δǫH

s

.

background image

56

3. L

2

Sobolev Spaces

Theorem 3.6 (SOBOLEV IMBEDDING THEOREM). For s >

n

2

+

66

k, H

s

C

k

. Further, we have

(3.7)

X

|α|≤k

sup

R

N

|D

α

f

| ≤ C

sk

|| f ||

(s)

.

Proof. Since S is dense in H

s

, it suffices to prove (3.7) for f ǫS . Let δ

x

denote the Dirac measure at x. For

|α| ≤ k, since s >

n

2

+

|α|, D

α

δ

x

ǫH

s

.



Since

< D

α

δ

x

, f >= (

−1)

|α|

< δ

x

, D

α

f >= (

−1)

|α|

(D

α

f )(x)

and

||D

α

δ

x

||

(

s)

is independent of x,

X

|α|≤k

sup

R

n

|D

α

f (x)

| =

X

|α|≤k

sup

R

n

|D

α

δ

x

, f >

|

X

|α|≤k

sup

R

n

|D

α

δ

x

||

(

s)

|| f ||

(s)

= C

sk

|| f ||

(s)

Now, given uǫH

s

, choose a sequence (u

j

) in S such that

||u u

j

||

(s)

converges to 0, as j tends to

∞. The above inequality with f = u

i

u

j

shows that (D

α

u

j

) is a Cauchy sequence in the uniform norm for

|α| ≤ k;

so its limit D

α

u is continuous.

Corollary 3.8. If uǫH

s

for all sǫR, then uǫC

i.e.,

T

s

H

s

C

.

This argument can be extended to show that if s >

n

2

+ k, then

elements of H

s

and their derivatives of order less than or equal to k are

67

not merely continuous but actually H ¨older continuous.

Proposition 3.9. If 0 < α < 1 and s =

n

2

+ α, then

||δ

x

− δ

y

||

(

s)

C

α

|x y|

α

Proof. We have

||δ

x

− δ

y

||

2
(

s)

=

Z

|e

−2πix

e

−2πiy

|

2

(1 +

|ξ|

2

)

s

dξ.

background image

1. General Theory of L

2

Sobolev Spaces

57

Let R be a positive number, to be fixed later. When

|ξ| ≤ R we use

the estimate

|e

−2πix

e

−2πiy

| ≤ 2π|ξ||xy| (by the mean value theorem)

and when

|ξ| > R, we use |e

−2πix

e−2πiy.ξ| ≤ 2. Then we have

||δ

x

− δ

y

||

2
(

s)

≤ 4π

2

|x y|

2

Z

|ξ|≤R

|ξ|

2

(1 +

|ξ|

2

)

s

dξ + 4

Z

|ξ|>R

(1 +

|ξ|

2

)

s

dξ

c

|x y|

2

R

Z

0

(1 + r

2

)

s

r

n+1

dr +

R

Z

0

r

−2s+n−1

dr

c

h

|x y|

2

R

−2s+n+2

+ R

−2s+n

i

as

n

2

< s <

n

2

+ 1

= c

h

|x y|

2

R

2

−2α

+ R

−2α

i

When we take R =

|x y|

−1

we get our result.



Exercise . Show that the above argument does not work when α = 1.
Instead, we get

||δ

x

− δ

y

||

(

s)

c|x y|| log |x y||

1/2

when x is near y.

What happens when α > 1?

Corollary 3.10. Let 0 < α < 1 and Λ

α

= bounded functions g : sup

x,y

68

|g(x) − g(y)|

|x y|

α

<

. If s =

n

2

+ α + l and f ǫH

s

, then D

β

f ǫΛ

α

for

|β| ≤ k.

Remark 3.11. We shall obtain an analogue of this result for L

p

norms

in Chapter 5.

The following lemma will be used in several arguments hereafter.

Lemma 3.12. For all ξ, nǫR

n

and sǫR, we have

" 1 + |ξ|

2

1 +

|η|

2

#

s

< 2

|s|

(1 +

|ξ − η|

2

)

|s|

.

Proof.

|ξ| ≤ |η| + |ξ − η| gives

|ξ|

2

< 2(

|η|

2

+

|ξ − η|

2

) so that

(1 +

|ξ|

2

) < 2(1 +

|η|

2

)(1 +

|ξ − η|

2

).



background image

58

3. L

2

Sobolev Spaces

If s > 0, then raise both sides to the s

th

power. If s < 0, interchange

ξ and η and raise to the

s

th

power.

Proposition 3.13. If φǫS , then the operator f

→ φ f is bounded for all

s.

Proof. The operator f

→ φ f is bounded on H

s

if and only if the opera-

tor of g

→ Λ

s

φΛ

s

g is bounded on L

2

, as one sees by



setting g =

s

f . But

69

(

s

φ

s

g)

(ξ) = (1 +

|ξ|

2

)

s/2

s

g)ˆ(ξ)

= (1 +

|ξ|

2

)

s/2

h ˆφ(ξ) ∗ (∧

s

g)

ˆ

(ξ)

i

= (1 +

|ξ|

2

)

s/2

Z

ˆ

φ(ξ

− η)(1 + |η|

2

g(η)dη

=

Z

ˆg(η)K(ξ, η)dη

where

K(ξ, η) = (1 + η

|

2

)

s/2

(1 +

|ξ|

2

)

s/2

ˆ

φ(ξ

− η).

By lemma 3.12,

|K(ξ, η)| < 2

|s|/2

(1 +

|ξ − η|

2

)

|s|/2

| ˆφ(ξ − η)|.

Therefore, since ˆ

φ is rapidly decreasing at

∞,

Z

|K(ξ, η)|dξ ≤ c for every ηǫR

n

,

Z

|K(ξ, η)|dη ≤ c for every ηǫR

n

.

Thus from Theorem 1.1, the operator with kernel K is bounded on

L

2

. Hence our proposition is proved.

The spaces H

s

are defined on R

n

globally by means of the Fourier

transform. Frequently, it is more appropriate to consider the following
versions of these spaces.

background image

1. General Theory of L

2

Sobolev Spaces

59

Definition 3.14. If

⊂ R

n

is open and sǫR, we define H

loc

s

Ω =

{ f ǫD

(Ω) :

∀Ω

⋐ Ω,

g

, ǫH

s

such that

g

= f on

}.

Proposition 3.15. f ǫH

loc

s

(Ω) if and only if φ f ǫH

s

for φǫC

0

(Ω).

70

Proof. If f ǫH

loc

s

(Ω) and φǫC

0

(Ω), then there exists gǫH

s

such that f =

g on supp φ. Therefore φ f = φgǫH

s

by proposition 3.13.



Conversely, if φ f ǫH

s

for all φǫC

0

(Ω

), and Ω

⋐ Ω, choose

φǫC

0

(

∞) with φ ≡ 1 on Ω

. Then φ f ǫH

s

and f = φ f on Ω

.

Corollary 3.16. If L =

P

|∝|≤k

a

(x)D

with a

ǫC

(Ω), then L maps

H

loc

s

(Ω) into H

loc
s

k

(Ω) for all sǫR.

It is a consequence of the Arzela-Ascoli theorem that if (u

j

) is a

sequence of C

k

functions such that

|u

j

| and |∂

α

u

j

|(|α| ≤ k) are bounded

on compact set uniformly in j, there exists a subsequence (v

j

) of (u

j

)

such that (∂

α

v

j

) converges uniformly on compact set for

|α| ≤ k − 1. In

particular, if the u

j

s are supported in a common compact set, then (∂

α

v

j

)

converges uniformly.

There is an analogue of this result for H

s

spaces.

Lemma 3.17. Suppose (u

k

) is a sequence of C

functions supported in

a fixed compact set such that sup

k

||u

k

||

(s)

<

. Then there exists a

subsequence which converges in the H

t

norm for all t < s.

Proof. Pick a φǫC

0

such that φ = 1 on Ω so that u

k

= φu

k

and hence

71

ˆu

k

= ˆφ

ˆu

k

. Then

(1 +

|ξ|

2

)

s/2

u

k

(ξ)

| = (1 + |ξ|

2

)

s/2

|

Z

ˆ

φ(ξ

− η)ˆu

k

(η)dη

|

Z

| ˆφ(ξ − η)||ˆu

k

(η)

|2

|s|/2

(1 +

|η|

2

)

s/2

(1 +

|ξ − η|

2

)

|s|/2

dη

≤ 2

|s|/2

||φ||

(

|s|)

||u

k

||

(s)

c

1

independent of k.

Likewise, we have

(1 +

|ξ|

2

)

s/2

|∂

j

ˆu

k

(ξ)

| ≤ 2

|s|/2

||2πx

j

φ(x)

||

|s|

||u

k

||

(s)

c

2

background image

60

3. L

2

Sobolev Spaces

independently of k. Therefore, by the Arzela-Ascoli theorem there ex-
ists a subsequence (ˆv

k

) of (ˆu

k

) which converges uniformly on compact

sets. For t

s,

||v

j

v

k

||

2
(t)

=

Z

(1 +

|ξ|

2

)

t

v

j

− ˆv

k

|

2

dξ

=

Z

|ξ|≤R

(1 +

|ξ|

2

)

t

v

j

− ˆv

k

|

2

dξ +

Z

|ξ|>R

(1 +

|ξ|

2

)

t

v

j

− ˆv

k

|

2

dξ

< (1 + R

2

)

max(t,0)

sup

|ξ|≤R

v

j

(ξ)

− ˆv

k

(ξ)

|

2

Z

|ξ|≤R

dξ+

+ (1 + R

2

)

t

s

Z

|ξ|>R

[(1 +

|ξ|

2

)

s

v

j

(ξ)

− ˆv

k

(ξ)

|

2

dξ]

c(l + R

2

)

n+

|t|

sup

|ξ|≤R

v

j

(ξ)

− ˆv

j

(ξ)

|

2

+ (1 + R

2

)

t

s

||v

j

v

k

||

2
(s)

.



Given ǫ > 0, choose R large enough so that the second term is less

than ǫ/2 for all j and k. This is possible since

||v

j

v

k

||

(s)

c and

72

t

s < 0. Then for j and k large enough the first term is less than ǫ/2,

since (ˆv

k

) converges uniformly on compact sets. Thus we see that (v

k

)

is a Cauchy sequence in H

t

and since H

t

is complete we are done.

Remark 3.12. Lemma 3.17 is false, if we do not assume that all the u

k

s

have support in a fixed compact set. For example, for uǫC

0

and x

k

ǫR

n

with

|x

k

| tending to ∞, define u

k

(x) = u(x

x

k

). Then the invariance of H

s

norms shows that

||u

k

||

(s)

=

||u||

(s)

. But no subsequence of (u

k

) converges,

in any H

t

. For, if a subsequence (v

k

) of (u

k

) converges, it must converge

to 0, since u

k

converges to 0 in S

. But then lim

||v

k

||

(t)

= 0 which is not

the case.

Theorem 3.20 (RELLICH THEOREM). Let H

0

s

(Ω) be the closure of

C

0

(Ω) in H

s

. If is bounded and t < s, the inclusion H

0

s

(Ω) ֒

H

t

is

compact, i.e., bounded sets in H

0

s

(Ω) are relatively compact in H

t

.

Proof. Let (u

k

) be a sequence in H

0

s

(Ω). To each k, find a v

k

ǫC

0

(Ω)

background image

1. General Theory of L

2

Sobolev Spaces

61

such that

||u

k

v

k

||

(s)

1

k

. Then we have

||v

k

||

(s)

≤ ||u

k

||

(s)

+

| ≤ c

(independent of k.) Therefore, by lemma 3.17, a subsequence (w

k

) of

(v

k

) exists such that (w

k

) converges in H

t

. If (u

k

) is the subsequence of

(u

k

) corresponding to the sequence (w

k

), we have

73

||u

i

u

j

||

(t)

<

||u

i

w

i

||

(t)

+

||w

i

w

j

||

(t)

+

||w

j

u

j

||

(t)

<

1

i

+

1

j

+

||w

i

w

j

||

(t)

→ 0 as ||i, j → ∞.

Hence (u

k

) converges in H

t

.



I the proof of the next theorem, we will use the technique of complex

interpolation, which is based on the following result from elementary
complex analysis called ‘Three lines lemma’.

Lemma 3.20. Suppose F(z) is analytic in o < ReZ < 1, continuous and
bounded on
0

ReZ ≤ 1. If |F(1 + it)| ≤ c

0

and

|F(l + it)| ≤ c

1

, then

F(s + it)

c

1

s

0

c

s

1

, for 0 < s < 1.

Proof. If ǫ > 0, the function

g

ǫ

(z) = c

z

−1

0

c

z

−1

1

e

ǫ(z

2

z)

f (z)

satisfies the hypotheses with c

0

and c

1

replaced by 1, and also

|g

ǫ

(z)

|

converges to 0 as

|Imz| → ∞ for 0 ≤ ReZ ≤ 1. From the maximum

modulus principle, it follows that

|g(z)| ≤ 1 for 0 ≤ Rez ≤ 1 and letting

ǫ tend to 0, we obtain the desired result.



Theorem 3.21. Suppose that

−∞ < s

0

< s

1

<

and T is a bounded

linear operator H

s

0

such that T

|H

s

1

is bounded on H

s

1

. Then T

|H

s

is

bounded on H

s

for all s with s

0

s s

1

.

Proof. Our hypothesis means that

74

s

0

T

s

0

and

s

1

T

s

1

are bounded operators on L

2

. For 0

Rez ≤ 1 we define

s

z

= (l

z)s

0

+ zs

1

and T

z

=

s

z

T

s

z

.



background image

62

3. L

2

Sobolev Spaces

Then what we wish to prove is that T

z

is bounded on L

2

for 0

z

1. Observe that when w = x + iy,

w

,

x

iy

and

(

iy

f )ˆ(ξ) = (1 +

|ξ|

2

)

iy/2

ˆ

f

2

(ξ) so that

|(∧

iy

f )ˆ(ξ) =

| ˆf(ξ)|.

Thus

iy

is unitary on H

s

for all s.

For φ, ψǫS , we define

F(z) =

Z

(T

z

φ)ψ =<

s

z

T

s

z

φ, ψ > .

Then

|F(z)| = | < ∧

s

z

T

s

z

φ, ψ >

|

=

| < T

s

z

φ,

s

z

ψ >

|

≤ ||T

s

z

φ

||

(

s

0

)

|| ∧

s

z

ψ

||

(

s

0

)

c|| ∧

s

z

φ

||

s

0

|| ∧

s

z

||

(

s

0

)

c||φ||

(s

0

s

1

)Re z

||ψ||

(s

1

s

0

)Rez

c||φ||

(0)

||ψ||

s

1

s

0

F(z) is clearly an analytic function of z for 0 < Re z < 1. Further, by our

75

hypothesis on T , when Rez = 0, we have

|F(z)| ≤ c

0

||φ||

(0)

||ψ||

(0)

and when Rez = 1, we have

|F(z)| ≤ c

1

||φ||

(0)

||ψ||

(0)

.

Therefore, by the Three lines lemma

|F(z)| ≤ c

1

z

0

c

z
0

||φ||

(0)

||ψ||

(0)

for 0 < z < 1.

Finally, by the self duality of H

0

= L

2

, this gives

||T

z

φ

|| ≤ c

1

z

0

c

z
0

||φ||

(0)

which completes the proof.

background image

1. General Theory of L

2

Sobolev Spaces

63

Remark 3.22. The same proof also yields the following more general
result:

Suppose

∞ < s

0

< s

1

<

∞ < t

0

< t

1

<

∞. If T is a bounded

linear operator from H

s

0

to H

t

0

whose restriction to H

s

1

to H

t

1

, then the

restriction of T to H

t

θ

is bounded from H

t

θ

to H

t

θ

for 0 < θ < 1 where

s

θ

= (1

− θ)s

0

+ θs

1

and t

θ

= (1

− θ)t

0

+ θt

1

.

As a consequence of this result, we obtain an easy proof that H

loc

s

is

invariant under smooth coordinate changes.

Theorem 3.23. Suppose and

are open subsets of R

n

and φ : Ω

is a C

diffeomorphism. Then the mapping f

f oφ maps H

loc

s

(Ω

).

76

continuously onto H

loc

s

(Ω).

Proof. The statement of the theorem is equivalent to the assertion that
for any φǫC

0

(Ω

), the map T f = (φ f )

◦ φ is bounded on H

s

for sǫR. If

s = 0, 1, 2, . . . , this follows from the chain rule and the fact that H

s

=

{ f : D

f ǫL

2

for

| ∝ | ≤ s}. By Theorem 3.21, it is true for all s ≥ 0. But

the adjoint of T is another map of the same form :



T

g = (ψg)

ψ where ψ = θ

−1

and ψ = φ

|J| ◦ Φ, J being the Jacobian

determinant of Φ

−1

. Hence T

is bounded on H

s

for all s

≥ 0 and by

duality of H

s

and H

s

, this yields the boundedness of T on H

s

for s < 0.

Finally, we ask to what extend the H

s

spaces include all distribu-

tions. Globally they do not, since, if f ǫH

s

, then f is tempered and ˆ

f is

a function. But locally they do, as we see from the following result.

Proposition 3.24. Every distribution with compact support lies in some
H

s

: i.e., E

S

s

∈R

H

s

.

Proof. If f ǫE

, then it is a continuous linear functional on C

. There-

fore, there exists a constant c > 0, a compact set K, and a nonnegative
integer k such that

| < f, φ > | ≤ c

X

|∝|≤k

sup

k

|D

φ

| for all φǫC

,

i.e.,

| < f, φ > | ≤ c P

|α|≤|k

sup

R

n

|D

α

φ

| for all φǫC

.

77

background image

64

3. L

2

Sobolev Spaces

By the Sobolev imbedding theorem,

X

|∝|≤k

sup

R

n

|D

α

φ

| ≤ c

||φ||

(k+

n
2

+ǫ)

for ǫ > 0.

Therefore,

| < f, φ > | ≤ c

′′

||(

k+

n
2

) for all φǫS since S is dense in

H

(k+

n
2

+ǫ)

, this shows that f is a continuous linear functional on H(

k+

n
2

).

Hence f ǫH

n
2

k−ǫ

.



Corollary 3.25. If f ǫD

(Ω) and

has compact closure in then there

exists s in R such that f ǫH

loc

s

(Ω

).

2 Hypoelliptic Operators With Constant

Coe

fficients

We now apply the machinery of Sobolev spaces to derive a criterion for
the hypoellipticity of constant coefficient differential operators. First,
we have a few preliminaries.

Definition 3.26. Let P be a polynomial in n variables. For a multi-index

∝, P

(

∝)

will ne defined by P

(

∝)

(ξ) = (

∂ξ

)

P(ξ).

Proposition 3.27. LEIBNIZ RULE When f ǫC

, gǫD

and P(D) is a

constant-coefficient partial differential operator of order k, we have

P(D)( f g) =

X

|∝|≤k

1

∝!

(P

(

∝)

(D)g)D

f ).

The proof of this proposition is left as an exercise to the reader.

Definition 3.28. We say that a polynomial P satisfies condition (H) if

78

there exists a δ > 0 such that

|P

(

∝)

(ξ)

|

|P(ξ)|

= 0 (

|ξ|

−δ|α|

) as

|δ| → ∞, ∀ ∝ .

Theorem 3.29 (H ¨

ORMANDER). If P satisfies condition (H), then P(D)

is hypoelliptic. More precisely, if f is in D

Ω)and P(D) f H

loc
s+kδ

(Ω),

where δ is as in condition (H) and k is the degree of P.

background image

2. Hypoelliptic Operators With Constant Coefficients

65

Proof. We first observe that the second assertion implies the first since
C

(Ω) =

T

sǫR

H

loc

s

(Ω) by Corollary 3.8.



Suppose therefore that P(D) f ǫH

loc

s

(Ω), f ǫD

(Ω). Given φǫC

(Ω)

we have to prove that φ f ǫH

s+kδ

. Let Ω

, be an open set such that

⋐ Ω and supp φ

⊂ Ω

. By Corollary 3.25, therefore exists t in R

such that f ǫH

loc

t

(Ω

). By decreasing t, we can some that t = s + k

1

mδ for assume positive integer m. Set φ

m

= φ and then choose

φ

m

−1

, φ

m

−2

, . . . , φ

0

, φ

−1 in C

0

(Ω

) such that φ

j

= 1 on supp φ

j+1

.

Then φ

j

P(D) f ǫH

S

H

t

k+l+ jδ

for 0

j m and φ

−1

f ǫH

t

. Now

P(D)(φ

0

f ) = φ

0

P(D) f +

X

∝,0

1

∝!

P

(D)(φ

−1

f )D

φ0

since θ

−1

= 1 on the support of φ

0

. So, P(D)(φ

0

f H

r

k+1

. This means

that

Z

(1 +

|ξ|

2

)

t

k+1

|P(ξ)(φ

0

f )ˆ(ξ)

|

2

dξ <

∞.

By condition (H)

79

Z

(1 +

|ξ|

2

)

t

k+l+δ|δ|

|P

∝|

(ξ)(φ

0

f )

(

ξ)

|

2

dξ <

This implies that P

(α)

(D)(φ

f H

t

k + l + δ|δ|.

Next,

P(D)(φ

1

f ) = φ

1

P(D) f +

X

∝,0

1

∝!

p

(

∝)

(D)(φ

o

f )D

1

f )

since φ

0

= 1 on the support of φ

1

, so P(D)(φ

1

f H

t

k+1+δ

. By the same

argument same argument as above,

P

(D)

1

f H

t+k+l+δ(1+

|∝|)

.

Continuing inductively, we obtain P(D)(φ

j

f H

t

k+l+ jδ

, which im-

plies that

P

(

∝)

(D)(φ

j

f H

t

k+l+δ( j+|∝|)

.

background image

66

3. L

2

Sobolev Spaces

For j = m, we above

P

(

∝)

(D)(φ

m

f H

t

k+l+δ(m+(δ))

= H

s

|∝|

.

If

P(ξ) =

X

|∝|≤k

a

ξ

,

choose

∝ with | ∝ | = k such that a

,

0. Then P

(ξ) =

∝!a

,

0.

whence φ f = φ

m

f ǫH

s+kδ

and we are done.

Remark 3.30. The condition (H) is equivalent to the following appar-

80

ently weaker condition

(H

) :

|p

(

∝)(ξ)

|

|P(ξ)|

→ 0 as |ξ| → ∞ for ∝, 0.

Condition (H

) is in turn equivalent to

(H

′′

) :

|Imζ| → ∞, in the set {ζǫƒ

n

: P(ζ) = 0

}.

The converse of H ¨ormander’s theorem is also true, i.e., hypoelliptic-

ity implies condition (H).

The proofs of these assertions can be found in H ¨ormander [6]. The

logical order of the proofs is

(H)

hypoellipticity ⇒ (H

′′

)

⇒ (H

)

⇒ (H).

The implication (H

)

⇒ (H) requires the use of some results from

(semi) algebraic geometry.

Definition 3.31. P(ξ) =

P

|∝|≤k

a

ξ

is called elliptic if

P

|∝|=k

a

ξ

,

0 for

every ξ , 0.

EXERCISES

1. Prove that P is elliptic if and only if

|P(ξ)| ≥ c

|ξ|

k

for large

|ξ|.

2. Prove that P is elliptic if and only if P satisfies condition (H) with

δ = 1.

background image

2. Hypoelliptic Operators With Constant Coefficients

67

3. Prove that no P satisfies condition (H) with δ > 1. (Hint: If

| ∝ | = k, P

(

∝)

is a constant).

4. Let P be elliptic and real valued. Define Q on R

n+1

by Q(ξ, τ) =

81

iτ + P(ξ), so that Q(D

x

, D

t

) = ∂

t

+ P(D

x

). Show that Q satisfies

condition (H) with δ = 1/k where k is the degree of P and that
1/k is the best possible value of δ.

(Hint : Consider the regions

|ξ|

k

≤ |τ| and |τ| ≤ |ξ|

k

separately).

background image
background image

Chapter 4

Basic Theory of Pseudo
Di

fferential Operators

1 Representation of Pseudo di

fferential Operators

Let L =

P

|α|≤k

a

α

(x)D

α

be a partial differential operator with C

coeffi-

82

cients on Ω. Using the Fourier transform, we can write

(Lu)(x) =

X

|α|≤k

a

α

(x)

Z

e

ix

ˆu(ξ)ξ

α

dξ

=

Z

e

ix

p(x, ξ)ˆu(ξ)dξ

where p(x, ξ) =

P

|α|≤k

a

α

(x

α

. This representation suggests that p(x, ξ)

can be replaced by more general functions. So, we make the following
definition.

Definition 4.1. For an open set

⊂ R

n

and a real number m we define

S

m

(Ω), the class of symbols of order m on , by

S

m

(Ω) =

{pǫC

(Ω

× R

n

) :

∀α, β, V

⊂ Ω, c = c

αβΩ

such that sup

xǫΩ

|D

β

x

D

α
ξ

p(x, ξ)

| ≤ c(1 + |ξ|)

m

−|α|

}.

69

background image

70

4. Basic Theory of Pseudo Differential Operators

We note that S

m

(Ω)

S

m

whenever m < m

, and we set S

−∞

(Ω) =

T

mǫR

S

m

(Ω).

Examples

(i) Let p(x, ξ) =

P

|α|≤k

a

α

(x

α

with a

α

ǫC

(Ω). Then pǫS

k

(Ω).

(ii) Let p(x, ξ) =

N

P

j=1

a

j

(x) f

j

(ξ) where a

f

ǫC

(Ω) and f

j

ǫC

(R

n

) is

83

homogeneous of degree m

j

large ξ : that is,

f

i

(rξ) = r

m

j

f

j

(ξ) for

|ξ| ≥ c, r ≥ 1.

In this case, pǫS

m

(Ω), where m = max

1

jN

{m

j

}.

(iii) Let p(x, ξ) = (1 +

|ξ|

2

)

s/2

. This p belongs to S

s

(R

n

).

(iv) Let p(x, ξ) = φ(ξ) sin log

|ξ| with φǫC

, φ = 0 near the origin and

φ = 1 when

|ξ| ≥ 1. Then pǫS

o

(R

n

).

Remark 4.2. We observe that when pǫS

m

(Ω), D

β

x

D

α
ξ

pǫS

m

−|α|

(Ω).

Further, if pǫS

m

1

(Ω) and qǫS

m

2

(Ω),then p + qǫS

m

(Ω), where m =

max

{m

1

, m

2

} and pqǫS

m

1

+m

2

(Ω).

Remark 4.3. Our symbol classes S

m

(Ω) are special cases of Horman-

der’s classes S

m
ρ,δ

(Ω). Namely, for 0

≤ δ ≤ ρ ≤ 1, and mǫR,

S

m
ρ,δ

(Ω) =

{pǫC

(Ω

× R

n

) :

∀α, β, V

⋐ Ω, c = c

αβΩ

such that

sup

xǫΩ

|D

β

x

D

α
ξ

p(x, ξ)

| ≤ c(1 + |ξ|)

m

−ρ|α|+δ|β|

}.

In this terminology, S

m

(Ω) = S

m
1,0

(Ω).

Definition 4.4. For pǫS

m

(Ω), we define the operator p(x.D) on the do-

main C

o

(Ω) by

p(x, D) u(x) =

Z

e

ix

p(x, ξ) ˆu(ξ)dξ, uǫC

o

(Ω).

background image

1. Representation of Pseudo differential Operators

71

(Sometimes, we shall denote p(x, D) by p). Operators of the form

p(x, D) with pǫS

m

(Ω) are called pseudo differential operators of order

84

m on Ω.

The set of all pseudo differential operators of order m on Ω will de-

noted by Ψ

m

(Ω). For brevity, we will sometimes write “ΨDO

′′

instead

of “pseudo differential operators”.

The next theorem states that p(x, D) is a continuous linear map of

C

o

(Ω) into C

(Ω which extends to E

(Ω). For the proof, we need a

result which depends on

Lemma 4.5. Let pǫS

m

(Ω) and φǫC

o

(Ω). Then, for each positive inte-

ger N, there exists c

N

> 0 such that for all ξ, η in R

n

,

|

Z

e

ix

p(x, ξ)φ(x)dx

| ≤ c

N

(1 +

|ξ|)

m

(1 +

|η|)

N

.

Proof. For any ξ and η in R

n

,

α

R

e

ix

p(x, ξ)φ(x)dx

|

=

|

Z

D

α

x

e

ix

p(x, ξ)φ(x)dx

|.



Integrating by parts, we have

α

Z

e

ix

p(x, ξ)φ(x)dx

| = |

Z

e

ix

D

α

x

(p(x, ξ)φ(x))dx

|

C

α

(1 +

|ξ|)

m

for all α.

Therefore

X

|α|≤N

α

Z

e

ix

p(x, ξ)φ(x)dx

| ≤ c

N

(1 +

|ξ|)

m

for all N.

Since (1 +

|η)

N

c P

|α|≤N

α

|, the required result follows.

85

background image

72

4. Basic Theory of Pseudo Differential Operators

Corollary 4.6. If pǫS

m

(Ω) and φǫC

o

(Ω), then the function

g(ξ) =

Z

e

ix

p(x, ξ)φ(x)dx

is rapidly decreasing as ξ tends to

.

Proof. Set ξ = η in the lemma.



Theorem 4.7. If pǫS

m

(Ω), then p(x, D) is a continuous linear map of

C

o

(Ω) into C

(Ω) which can be extended as a linear map from E

(Ω)

into D

(Ω).

Proof. For uǫC

o

(Ω),

p(x, D)u(x) =

Z

e

ix

p(x, ξ)ˆu(ξ)dξ.

The integral converges absolutely and uniformly on compact sets, as

do the integrals

Z

D

α

x

(e

ix

p(x, ξ) ˆu(ξ)dξ for all α,

since pǫS

m

and ˆuǫ

S.



This proves that p(x, D)uǫC

(Ω), and continuity of p(x, D) from

C

o

(Ω) to C

(Ω) is an easy exercise.

To prove the rest of the theorem, we will make use of Corollary 4.6.
For uǫE

(Ω), we define p(x, D)u as a functional on C

o

(Ω) as fol-

lows:

< p(x, D)u, φ > =

x

p(x, ξ)ˆu(ξ)e

ix

φ(x)dxdξ

=

Z

g(ξ)ˆu(ξ)dξ

where

86

g(ξ) =

Z

e

ix

p(x, ξ)φ(x)dx, for φǫC

o

(Ω).

background image

2. Distribution Kernels and the Pseudo Local Property

73

By Corollary 4.6, g(ξ) is rapidly decreasing, while ˆu is of polyno-

mial growth; so the last integral is absolutely convergent and the func-
tional on C

o

thus defined is easily seen to be continuous. Moreover, if

uǫC

o

, the double integral is absolutely convergent, and by interchang-

ing the order of integration, we see that this definition of p(x, D)u coin-
cides with the original one. Hence we have extended p(x, D) to a map
from E

(Ω) to

D(Ω).

Remark 4.8. It follows easily from the above argument that p(x, D) is
sequentially from E

(Ω) to

D(Ω), i.e., if u

k

converges to u in E

(Ω), then

p(x, D)u

k

converges to p(x, D)u in

D

(Ω). Actually, p(x, D) is contin-

uous from E

(Ω) to

D

(Ω) : this follows from the fact (which we shall

prove later) that the transpose of a pseudo differential operator is again
a pseudo differential operator. Thus p(x, D)

t

: C

o

(Ω)

C

(Ω) is con-

tinuous, and so, by duality, p(x, D) = (p(x, D)

t

)

t

: E

(Ω)

→ D

(Ω) is

continuous.

2 Distribution Kernels and the Pseudo Local Prop-

erty

Definition 4.9. Let T be an operator from C

o

(Ω) to C

(Ω). If there

exists a distribution K on

× Ω such that

< T u, v >=< k, v

u > for u, vǫC

o

(Ω),

we say that K is the distribution kernel of the operator T .

In this definition, v

u is defined by (vu)(x, y) = v(x)u(y). Formally,

this definition says that.

T u(x) =

Z

K(x.y)u(y)dy.

K is uniquely determined since linear combinations of functions of the

87

form v

u are dense in C

o

(Ω

× Ω).

If pǫS

m

(Ω), it is easy to compute the distribution kernel of p(x, D).

In fact,

< p(x, D)u, v > =

x

p(x, ξ)ˆu(ξ)e

ix

v(x)dξdx

background image

74

4. Basic Theory of Pseudo Differential Operators

=

x

p(x, ξ)e

ix

(v

u)

ˆ

2

(x, ξ)dξdx

where (v

u)

ˆ

2

(x, ξ) means the Fourier transform in the second variable.

It follows immediately from the definition of K that

< K, w > =

x

e

ix

p(x, ξ)w

ˆ

2

(x, ξ)dξdx,

wǫC

o

(Ω

× Ω).

or < K, w > =

y

e

i(x

y).ξ

p(x, ξ)w(x, y)dy dξ dx.

From this it is easy to see that K(x, y) = p

v
2

(x, x

y) where p

v
2

(x, .) is

the inverse Fourier transform of the tempered distribution p(x, .). In par-
ticular, this shows that p is uniquely determined by the operator p(x, D).
If PǫΨ

m

(Ω), we shall sometimes denote the unique pǫS

m

(Ω) such that

p = p(x, D) by σ

p

.

The following theorem gives precise results on the kernel K of

p(x, D).

Theorem 4.10. The distribution kernel K of p(x, D) with pǫS

m

(Ω) is in

C

on (Ω

× Ω)/∆ where ∆{(x, x) : xǫΩ} is the diagonal. More precisely,

if

|α| > M + n + j for a positive integer j, then (x y)

α

K(x, yC

j

(Ω

× Ω).

Proof. For wǫC

o

(Ω

× Ω), let us compute < (x y)

α

K, w >.

88

< (x

y)

α

K, w > =< K, (x

y)

α

w >

=

x

e

ix

p(x, ξ)(x + D

ξ

)

α

w

ˆ

2

(x, ξ)dξ dx

=

x

w

ˆ

2

(x, ξ)(x

D

ξ

)

α

{e

ix

p(x, ξ)

}dξ dx.



Using the Leibniz formula, it is easily seen that

(x

D

ξ

)

α

{e

ix

p(x, ξ)

} = (−D

ξ

)

α

p(x, ξ).e

ix

Therefore

< (x

y)

α

K, w > =

x

w

ˆ

2

(x, ξ)e

ix

(

D

ξ

)

α

p(x, ξ)dξ dx

=

y

w(x, y)e

i(x

y).ξ

(

D

ξ

)

α

p(x, ξ)dy dξ dx.

background image

2. Distribution Kernels and the Pseudo Local Property

75

From the above expression, we infer that

(x

y)

α

K(x, y) =

Z

e

i(x

y).ξ

(

D

ξ

)

α

p(x, ξ)dξ.

Since

|(−D

ξ

)

α

p(x, ξ)

| ≤ c(1 + |ξ|)

m

−|α|

, the integral converges abso-

lutely and uniformly on compact sets whenever m

− |α| < −n. Also

we can differentiate with respect to x and y

j times under the integral

sign provided m

− |α| + j < −n or |α| > m + n + j. Thus, we see that

(x

y)

α

KǫC

j

(Ω

× Ω).

Corollary 4.11 (PSEUDO LOCAL PROPERTY OF ΨDO). If PǫΨ

m

(Ω), then, for all uǫE

(Ω), sing supp Pu is contained in sing supp u.

Proof. Let uǫE

and let V be an arbitrary neighbourhood of sing supp u.

Take a φǫC

o

(V) such that φ = 1 on sing supp u. Then u = φu+(1

−φ)u =

u

1

+ u

2

; u

2

is a C

o

function and supp u

1

ǫV.



Therefore, Pu = Pu

1

+ Pu

2

and Pu

2

is a C

function. Moreover,

when x

o

<

V,

89

Pu

1

(x) =

Z

V

K(x, y)u

1

(y)dy

is also a C

function in a neighbourhood of x

o

since u

1

(y) = 0 for y near

x

o

. This implies that sing supp Pu

V. Since V is any neighbourhood

of sing supp u, the corollary is proved.

Corollary 4.12. If pǫS

−∞

(Ω), then the distribution kernel K of p(x, D)

is in C

(Ω

× Ω).

Proof. Follows from the theorem if we take

|α| = 0.



Corollary 4.13. If pǫS

−∞

(Ω), then p(x, D), maps E

(Ω) continuously

into C

(Ω).

Proof. If uǫE

(Ω), in view of Corollary 4.12, it is easily seen that

p(x, D)u is a smooth function defined by

P(x, D)u(x) =

Z

K(x, y)u(y)dy =< u, K(x, .) >,

whence the result follows.



background image

76

4. Basic Theory of Pseudo Differential Operators

Definition 4.14. A smoothing operator is a linear operator T which
maps E

(Ω) continuously into C

(Ω).

If KǫC

(Ω

× Ω), then the operator T defined by

(T f )(x) =< K(x, .), f >=

Z

K(x, y) f (y)dy

is a smoothing operator. Conversely, every smoothing operator T can
be given in the above form with K
(x, y) = (T δ

y

)(x).

As we have already remarked if pǫS

−∞

, then the corresponding

ΨDO is smoothing. However, not every smoothing operator is a ΨDO.
For example, we have

Proposition 4.15. Suppose p(x, .) is a C

function of x with values in

90

E

. Then p(x, D) (defined in the same way as in the case of pǫS

m

(Ω)) is

a smoothing operator.

Proof. The distribution kernel K of p(x, D) is given by

K(x, y) =

Z

e

i(x

y).ξ

p(x, ξ)dξ

=< p(x, .), e

i(x

y).(.)

>

and hence K(x, yC

. Thus K defines a smoothing operator.



Remark 4.16. Sometimes, it is convenient to enlarge the class of
pseudo-differential operators of order m by including operators of the
form P + S where PǫΨ

m

(Ω) and S is smoothing. However, the general

philosophy is the following:

1. On the level of operators, smoothing operators are negligible.

2. On the level of symbols, what counts is the asymptotic behaviour

at

∞, so that symbols of order −∞ are negligible.

background image

3. Asymptotic Expansions of Symbols

77

3 Asymptotic Expansions of Symbols

Definition 4.17. Suppose m

0

> m

1

> m

2

>

· · · m

j

ǫR, m

j

tends to

−∞,

and p

j

ǫS

m

j

(Ω), pǫS

m

o

(Ω). We say that p

P

j=0

p

j

if p

P

j<k

p

j

ǫS

m

k

(Ω)

for all k.

Proposition 4.18. Suppose m

j

tends to

−∞ and p

j

ǫS

m j

(Ω). Then there

exists a p in S

m

o

(Ω) such that p

P j = 0p

j

. This p is unique modulo

S

−∞

(Ω).

Proof. Let (Ω

n

) be an increasing sequence of compact subsets of Ω

whose union is Ω. Fix φǫC

with φ = 1 for

|ξ| ≥ 1 > and φ = 0

for

|ξ| ≤

1

2

.



Claim

There exists a sequence (t

j

),t

j

≥ 0 and t

j

tending to

∞ so

91

rapidly that we have
(4.19)
|D

β

x

D

α
ξ

(φ(ξ/t

j

)p

j

(x, ξ))

| ≤ 2

j

(1+

|ξ|)

m

j

−1

−|α|

for xǫΩ

i

, and

|α|+|β|+i j.

Granted this, we define

p(x, ξ) =

X

j=0

φ(ξ/t

j

)p

j

(x, ξ).

Note that for each x and ξ, the sum is finite. Using (4.19), it is

straightforward to show that pǫS

m

o

(Ω) and p

P

j=0

p

j

. Moreover, sup-

pose that qǫS

m

o

(Ω) and q

P

j=0

p

j

. Then

p

q = (p

X

j<k

p

j

)

− (q

X

j<k

p

j

S

m

k

(Ω) for all k.

Hence p

qǫS

−∞

.

we now briefly indicate the steps involved in proving the claim :

i) First observe that

|D

ν
ξ

φ(ξ/t

j

)

| ≤ c|ξ|

−|ν|

uniformly in j.

background image

78

4. Basic Theory of Pseudo Differential Operators

ii) Hence we have

|D

β

x

D

α
ξ

(φ(ξ/t

j

)p

j

(x, ξ))

| ≤ c

j

(1 +

|ξ|)

m

j

−|α|

for xǫΩ

i

,

|α| + |β| + i j.

iii) Finally, pick t

j

so large that “

|ξ| ≥ t

j

/2

c

j

(1 +

|ξ|)

m

j

m

j

−1

≤ 2

j

Details are left to the reader.

The following theorem provides a useful criterion for the asymptotic

relation p

P

j=0

p

j

to hold.

Theorem 4.20. Suppose p

j

ǫS

m

j

, m

j

tends to

−∞, and pǫC

(Ω

× R

n

)

92

satisfies following conditions:

i) for all α and β and all

⋐ Ω, there exists c > 0, µǫRsuch that

|D

β

x

D

α
ξ

p(x, ξ)

| ≤ c(1 + |ξ|)

µ

, xǫΩ, and

ii) there exists a sequence

k

), µ

k

tending to

−∞ so that |p(x, ξ) −

P

j<k

p

j

(x, ξ)

| ≤ C

, (1 +

|ξ|)

µ

k

for xǫΩ

. Then pǫS

m

o

(Ω) and p

P

j=0

p

j

.

To prove this theorem, we need the following

Lemma 4.21. Let

1

and

2

be two compact subsets of R

n

such that

1

is contained in the interior of

2

. Then there exists constants c

1

> 0 and

c

2

> 0 such that for all f ǫC

2

(Ω

2

),

Sup

1

|∂

j

f

|

2

c

1

(Sup

2

| f |

2

) + c

2

(Sup

2

| f |)(Sup

2

|∂

2

j

f

|).

Proof. It suffices to assume that n = 1 and f is real valued. With this
reduction, the proof becomes an exercise in elementary calculus. This
idea is roughly as follows: We wish to show that if

k f k

and

k f

′′

k

are both small, then so is

k f

k

is small and

| f

(x

o

)

| is large, then

| f

| will be large in some sizable interval [a, b] containing x

o

. But then

| f (b) − f (a)| and hence k f k

is large. We leave it to the reader to work

out the quantitative details of this argument.



background image

4. Properly Supported Operators

79

Proof of the theorem . By Proposition 4.18, there exists q in S

m

o

such

93

that q

P

j=0

p

j

; so, it will suffice to show that p

qǫS

−∞

. First,

|p(x, ξ) − q(x, ξ)| = |(p(x, ξ) −

X

j<k

p

j

(x, ξ))

− (q(x, ξ) −

X

j<k

p

j

(x, ξ))

|.

For

⊂⊂ Ω, |p(x, ξ) −

X

j<k

p

j

(x, ξ)

| ≤ c

, (1 +

|ξ|)

µ

k

, xǫΩ

.

Also q

P

j=0

p

j

. These show that, for any N,

|p(x, ξ)| − q(x, ξ)| ≤ C

N

(1 +

|ξ|)

N

, xǫΩ

.

We want to prove such an estimate for D

β

x

D

α
ξ

(p

q) also. To this

end, we will apply Lemma 4.21 to the function (x, η)

→ (p q)(x, ξ +

η) considering ξ as a parameter, and taking Ω

1

= Ω

× 0 and Ω

2

a

small neighbourhood of Ω

1

. If

|α| + |β| = 1, we use the estimate just

established for p

q, together with the hypothesis (i) on the second

order derivatives of p; the lemma implies that D

β

x

D

α
ξ

(p

q) is rapidly

decreasing for

|α| + |β| = 1. Combining this with the hypothesis (i) on

the third order derivatives of p, we see that D

β

x

D

α
ξ

is rapidly decreasing

for

|α| + |β| = 2. Proceeding by induction on |α| + |β| we get the required

result.

4 Properly Supported Operators

Since pseudo differential operators map C

o

to C

rather then C

o

, it

is generally not possible to compose two of them. The problem may
be remedied by considering a more restricted class of operators, the so
called “properly supported” ones.

Definition 4.22. A subset K of

× Ω is said to be proper if, for any

94

compact set

⊂ Ω, both π

−1

1

(Ω

)

K and π

−1

2

(Ω

)

K are compact.

background image

80

4. Basic Theory of Pseudo Differential Operators

Here π

1

and π

2

are projections of Ω

× Ω onto the first and second

factors.

For example the diagonal ∆ =

{(x, x) : xǫΩ} is proper. Most of the

proper subsets we will be considering are neighbourhoods of subsets of
the diagonal.

Figure 4.1: A proper set

Definition 4.23. An operator T : C

o

(Ω)

C

(Ω) is said to be prop-

erly supported, if its distribution kernel K has proper support.

Exercise. Let T =

P a

α

(x)D

α

be a differential operator on . Computer

the distribution kernel K of T and show that supp K is a subset of the
diagonal. Thus T is properly supported.

If T is properly supported then it maps C

o

into itself, since supp

T u

⊂ π

1

−1

2

(supp u)

∩ supp K), as is easily seen from the formula

T u(x) =

R

K(x, y)u(y)dy. More generally, for any Ω

⊂⊂ Ω, there exists

′′

⊂⊂ Ω such that the values of T u on Ω

depend only on the values of

u on Ω

′′

namely, Ω

′′

= π

2

−1

(Ω

)

∩ supp K). From this it follows that

T can be extended to a map from C

(Ω) to itself. In fact, if uǫC

(Ω)

and Ω

, Ω

′′

are as above, we define Tu on Ω

by

T u

|

= T u)

|

,

where φǫC

(Ω) and φ = 1 on Ω

′′

. This definition is independent of the

95

choice of φ and gives the same result on the intersection of two Ω

s; so

T u is well defined on all of Ω.

If T is a properly supported pseudo differential operator, so that T

extends to a map from E

(Ω) to

D

(Ω), the same arguments show that T

maps E

(Ω) into itself and extends further to a map of

D

(Ω) to itself.

Suppose T and S are properly supported operators on C

(Ω) with

distribution kernels K and L which are C

off the diagonal. Then T S

background image

5. ψdo

s

Defined by Multiple Symbols

81

is an operator on C

(Ω) with distribution kernel M formally given by

M(x, y) =

R

K(x, z)L(z, y)dz. In fact, if x , y, since K(x, .) and L(., Y)

are smooth except at x and y, the product K(x, .)L(., Y) is a well defined
element of E

, and the formula M(x, y) =< K(x, .)L(., y), 1 > displays M

as a C

function off the diagonal.

Proposition 4.24. Supp M is proper. Thus, T S is properly supported.

Proof. Clearly,

Supp M

⊂ {(x, y) : π

2

−1

1

(x)

∩ supp K) ∩ π

1

−1

2

(y)

∩ supp L) , φ}

Suppose A

⊂ Ω is compact and set B = π

2

−1

1

(A)

∩ supp K), then B is

compact and

π

−1

1

(A)

∩ supp M A ×

n

y : B

∩ π

1

−1

2

(y)

∩ supp L) , φ

o

= A

×

n

y : π

−1

1

(B)

∩ π

−1

2

(y)

∩ supp L , φ

o

= A

× π

2

−1

1

(B)

∩ supp L)

which is compact. Likewise π

−1

2

(A)

∩ supp M is also compact.



Exercise . Suppose A

⊂ (Ω × Ω) is proper. Show that there exists a

96

properly supported φǫC

(Ω

×Ω) such that φ = 1 on A. (Hint: Let

j

} ⊂

C

o

(Ω

× Ω) be a partition of unity on Ω × Ω and let φ =

P

A

∩supp φ

j

= φ

φ

j

).

5 ψdo

s Defined by Multiple Symbols

We have arranged our definition of pseudo differential operators to agree
with the usual convention for differential operators, according to which
differentiations are performed first, followed by multiplication by the
coefficients. However, in some situations (for example, computing ad-
joints), it is convenient to have a more flexible setup which allows mul-
tiplication operators both before and after differentiations. We there-
fore, introduce the following apparently more general class of operators
(which, however, turns out to coincide with the class of ψDO, modulo
smoothing operators).

background image

82

4. Basic Theory of Pseudo Differential Operators

Definition 4.25. For open in R

n

and m real, we define the class of

multiple symbols of order m on Ω,

S

m

(Ω

× Ω) = {aǫC

(Ω

× R

n

× Ω) : ∀α, β, ∀Ω

⊂⊂ Ω, ∃c = c

αβΩ

such that sup

x,yǫΩ

|D

β

x,y

D

α
ξ

a(x, ξ, y)

| ≤ c(1 + |ξ|)

m

−|α|

}

When aǫS

m

(Ω

× Ω) we define the operator A = a(x, D, y) by

Au(x) =

x

e

i(x

y).ξ

a(x, ξ, y)u(y)dydξ

Here the integral must be interpreted as an iterated integral with in-

tegration performed first in y then in ξ as it is not absolutely convergent
as a double integral.

We observe that if a(x, ξ, y) = a(x, ξ) is independent of y, then

97

A = a(x, D); thus this class of operators include the ψDO

s. We also

observe that different a

s may give rise to same operator. For example,

if a(x, ξ, y) = φ(x)ψ(y) with ψ, φǫC

o

(Ω) and supp ψ

∩ supp φ = φ, then

aǫ s

o

(Ω

× Ω) and a(x, D, y) = 0.

Definition 4.26. Given aǫS

m

(Ω

× Ω) we define

X

a

=

{(x, y)ǫΩ × Ω : (x, ξ, y)ǫ supp a for some ξǫR

n

}.

Proposition 4.27. Let aǫS

m

(Ω

× Ω) and let K be the distribution kernel

of A = a(x, D, y). Then

i) Supp K

⊂ P

a

, and

ii) if the support of K is proper, then there exists a

ǫS

m

(Ω

×Ω) such that

a

(x, ξ, y) = a(x, ξ, y) when (x, y) is near the diagonal in

× Ω, P

a,

is proper and a(x, D, y) = a

(x, D, y).

Proof. The kernel K is given by

< K, w >=

y

e

i(x

y).ξ

a(x, ξ, y)w(x, y)dydξdx.

From this, we see that < K, w > . = 0 whenever supp w

∩ P

a

= φ.

Therefore, supp K

⊂ P

a

. This proves (i).



background image

5. ψdo

s

Defined by Multiple Symbols

83

To prove (ii), choose a properly supported φǫC

(Ω

× Ω) with φ = 1

on ∆

∪ supp K, ∆ being the diagonal. Set a(x, ξ, y) = φ(x, y)a(x, ξ, y).

Then

P

a

,

⊂ supp φ and hence P

a

, is proper. Also a

(x, ξ, y) = a(x, ξ, y)

when (x, y) is near the diagonal. Now,

< a(x, D, y)u, v > =< K, v

u >

=< φK, v

u >

=< K, φ(v

u) >

=

y

e

i(x

y)·ξ

a(x, ξ, y)φ(x, y)

× v(x)u(y)dydξdx

=

y

e

i(x

y).ξ

a

(x, ξ, y)u(y)dydξdx

=< a

(x, D, y)u, v > .

98

Since v is arbitrary, we have a(x, D, y) = a

(x, D, y).

Theorem 4.28. Suppose aǫS

m

(Ω

× Ω) and A = a(x, D, y) is prop-

erly supported. Let p(x, ξ) = e

−2πix

A(e

i(.).ξ

)(x). Then pǫS

m

(Ω) and

p(x, D) = A. Further,

p(x, ξ)

X

a

1

α!

α
ξ

D

α
y

a(x, ξ, y)

|

y=x

.

Proof. For u in C

o

(Ω), we have

u(x) =

Z

e

ix

ˆu(ξ)dξ.

Therefore, by linearity and continuity of A,

Au(x) =

Z

A(e

i(.).ξ

)(xu(ξ)dξ

=

Z

e

ix

p(x, ξ)ˆu(ξ)dξ = p(x, D)u(x).



background image

84

4. Basic Theory of Pseudo Differential Operators

By the proposition above, we can modify a so that

P

a

is proper with-

out affecting A or the behaviour of a along the diagonal x = y (which
is what enters into the asymptotic expansion of p). Set b(x, ξ, y) =
a(x, ξ, x + y). Then as a function of y, b has compact support. Indeed,

99

for fixed x and ξ, if yǫ

{y

: b(x, ξ, y

) , 0

} then x + yǫπ

2

−1

1

(x)

∩ P

a

)

which is compact. Now,

p(x, η) = e

−2πix

x

e

i(x

y).ξ

b(x, ξ, y

x)e

iy

dydξ

= e

−2πix

x

e

iz

b(x, ξz)e

i(z+x).η

dzdξ

=

Z

ˆb

3

(x, ξ, ξ

− η)dξ =

Z

ˆb

3

(x, ξ + η, ξ)dξ,

where ˆb

3

is the Fourier transform of b in the third variable. Since

b(x, ξ, .) is in C

o

, this calculation is justified and ˆb

3

(x, η, ξ) is rapidly

decreasing in the variable ξ.

More precisely, since aǫS

m

(Ω

× Ω), we have

(4.29)

|D

β

x

D

α
η

ˆb

3

(x, η, ξ)

| ≤ c

αβN

(1 +

|η|)

m

−|α|

(1 +

|ξ|)

N

.

Since

1 + |η|

1 +

|ξ|

!

s

≤ (1 + |ξ − η|)

|s|

, taking η + ξ instead of η

we have (1 +

|ξ + η|)

s

≤ (1 + |ξ|)

s

(1 +

|η|)

|s|

. If we use this in

|D

β

x

D

α
η

ˆb

3

(x, ξ + η, ξ)

| ≤ c

αβN

(1 +

|ξ + η|)

m

−|α|

(1 +

|ξ|)

N

we get

|D

β

x

D

α
η

ˆb

3

(x, ξ + η, ξ)

| ≤ c

αβN

(1 +

|ξ|)

N+m−|α|

(1 +

|η|)

|m|+|α|

.

If we take N so that

N + m − |α| < −n we have

|D

β

x

D

α
η

p(x, η)

| = |

Z

D

β

x

D

α
η

ˆb

3

(x, ξ + η, ξ)dξ

|

(4.30)

c

αβ

(1 +

|η|)

|m|+|α|

background image

5. ψdo

s

Defined by Multiple Symbols

85

On the other hand, taking the Taylor expansion of ˆb

3

(x, ξ + η, ξ) in

100

the middle argument about the point η, (4.29) gives

b

3

(x, ξ + η, ξ)

X

|α|<k

α
η

ˆb

3

(x, η, ξ)

ξ

α

α!

|

C|ξ|

k

sup

|α|=k

0

t≤1

|∂

α
η

ˆb

3

(x, η + tξ, ξ)

|

C

N

|ξ|

k

sup

0

t≤1

(1 +

|η + tξ|)

m

k

(1 +

|ξ|)

N

.

When

|ξ| <

1
2

|η|, taking N = k we get

b

3

(x, ξ + η, ξ)

X

|α|<k

α
η

ˆb

3

(x, η, ξ)

ξ

α

α!

| ≤ c

k

(1 +

|η|)

m

k

.

When

|ξ| ≥

1

2

|η|, we have

b

3

(x, ξ + η, ξ)

X

|α|<k

α
η

ˆb

3

(x, η, ξ)

ξ

α

α!

| ≤ C

N

(1 +

|ξ|)

m+k

N

.

(Actually, the exponent of (1 +

|ξ|)) can be taken as m N when

m

k ≥ 0 and when m k < 0). Also we have

Z

α
η

ˆb

3

(x, η, ξ)ξ

α

dξ = D

α
Y

Z

e

iY

·ξ

α
η

ˆb

3

(x, η, ξ)dξ

|

Y=0

= D

α
Y

α
η

b(x, η, y)

|

Y=0

= D

α
Y

α
η

a(x, η, y)

|

y=x

.

So finally

|p(x, η) −

X

|α|<k

D

α
y

α
η

a(x, η, y)

1

α!

|

y=x

|

=

|

Z

ˆb

3

(x, ξ + η, ξ)dξ

X

|α|<k

1

α!

Z

α
η

ˆb

3

(x, η, ξ)ξ

α

dξ

|

background image

86

4. Basic Theory of Pseudo Differential Operators

C

k

Z

|ξ|<

|η|

2

(1 +

|η|)

m

k

dξ + C

N

Z

|ξ|≥

|η|

2

(1 +

|ξ|)

m+k

N

dξ.

Now

101

Z

|ξ|<

|η|

2

C

k

(1 +

|η|)

m

k

dξ = C

k

(1 +

|η|)

m

k

|η|

n

C

k

(1 +

|η|)

m

k+n

In the second integral, we take N > n + m + k so that

Z

|ξ|

≥ |η|

2

C

N

(1 +

|ξ|)

m

N

dξ

C

N

(1 +

|η|)

m

N+n

,

by the usual integration in polar coordinates. Therefore, since N

k we

obtain

|p(x, η) −

X

|α|<k

1

α!

D

α
y

σ

α
η

a(x, η, y)

|

x=y

| ≤ C

k

(1 +

|η|)|

µ

k

with µ

k

= m

k + n.

Combining this with (4.30), we see that all the conditions of Theo-

rem 4.20 are satisfied, so we are done.

Corollary 4.31. If Pǫψ

m

, there exists a Qǫψ

m

such that Q is properly

supported and P

Q is smoothing.

Proof. If P = p(x, D) choose φǫC

(Ω

× Ω) which is properly sup-

ported and φ = 1 near the diagonal. Set a(x, ξ, y) = φ(x, y)p(x, ξ). Then
aǫS

m

(Ω

× Ω) and by construction

P

a

is proper. So a(x, D, y) = Qǫψ

m

.

If K

p

and K

Q

denote the distribution kernels of P and Q then it is eas-

ily seen that K

Q

= φK

p

. This shows that K

Q

K

p

vanishes near the

diagonal: so, by Theorem 4.10, K

Q

K

p

is C

. Therefore, Q

p is

smoothing.



Remark 4.32. The idea of multiple symbols can clearly be generalised.
For example, if aǫC

(Ω

× R

n

× Ω × R

n

) satisfies estimates of the form

|D

v

x,y

D

α
ξ

D

β
η

a(x, ξ, y, η)

| ≤ C(1 + η|)

m

2

−|β|

(1 +

|ξ|)

m

1

−|α|

, one can define an

102

operator A = a(x, D, y, D) by

background image

6. Products and Adjoint of ψDO

S

87

Au(x) =

y

e

i(x

y.ξ+y.η)

a(x, ξ, y, η)ˆu(η)dηdydξ

which agree with our previous definitions, if a is independent of the last
one or two variables. As one would expect, under reasonable restriction
on supp a, it turns out that Aǫψ

m

1

+m

2

(Ω). The reader may wish to amuse

himself by working this out and computing the symbol of A.

6 Products and Adjoint of ψDO

S

We are now in a position to compute products and adjoints of ψDO

S .

First we clarify our terminology. Let T S be linear operators from
C

o

(Ω) into C

(Ω). We say that S = T

t

if < T u, v >=< u, S u >

for u, vǫC

o

(Ω) and we say that S = T

if < T u, ¯v >=< u, ¯

S v > for

u, vǫC

o

(Ω).

Remark 4.33. If T is properly supported, then T

t

and T

are also prop-

erly supported. Indeed, if K is the distribution kernel of T , then the dis-
tribution kernel of T

t

is K

t

and that of T

is K

where K

t

(x, y) = K(y, x)

and K

(x, y) = K(y, x).

We recall that if Pǫψ

m

(Ω), we denote the corresponding symbol in

S

m

(Ω) by σ

p

.

Theorem 4.34. If Pǫψ

m

(Ω) is properly supported, then P

t

, P

ǫψ

m

(Ω)

and σ

p

t(x, ξ)

P

α

(

−1)

|α|

α!

α
ξ

D

α

x

σ

p

(x,

−ξ),

σ

p

(x, ξ)

X

a

1

α!

α
ξ

D

α

x

¯

σ

p

(x, ξ).

Proof. For u, vǫC

o

(Ω), we have

103

< Pu, v > =

x

e

ix

p(x, ξ)ˆu(ξ)v(x)dξdx

=

Z

Z

e

ix

p(x, ξ)v(x)dx

!

ˆu(ξ)dξ

background image

88

4. Basic Theory of Pseudo Differential Operators

=

Z

g(ξ)ˆu(ξ)dξ

where

g(ξ) =

Z

e

ix

p(x, ξ)v(x)dx.



By Corollary 4.6, g us a rapidly decreasing function. Thus

< Pu, v >=

Z

gˆu =

Z

uˆg =< u, P

t
v

>

so that

P

t

v(y) = ˆg(y) =

x

e

i(x

y).ξ

p(x, ξ)v(x)dxdξ

(4.35)

=

x

e

i(y

x).ξ

p(x,

−ξ)v(x)dxdξ

= a(x, D, y)v(y)

where a(x, ξy) = p(y.

− ξ). Therefore, by Theorem 4.28,

P

t

ǫψ

m

(Ω) and σ

P

t

(x, ξ)

X

a

(

−1)

|α|

α!

α
ξ

D

α

x

σ

P

(x,

−ξ).

The assertions about P

follows along similar lines

Theorem 4.36. If Pǫψ

m

1

(Ω) and Qǫψ

m

2

(Ω) are properly supported,

then QPǫψ

m

1

+m

2

(Ω) and

σ

QP

(x, ξ)

X

α

1

α!

α
ξ

σ

Q

(x, ξ)D

α

x

σ

p

(x, ξ).

Proof. Since P = (P

t

)

t

, we have

104

(Pu)(x) =

x

e

i(x

y)·ξ

σ

t
P

(y,

−ξ)u(y)dydξ

by (4.35).



background image

6. Products and Adjoint of ψDO

S

89

In other words

(Pu)ˆ(ξ) =

Z

e

−2πiy·ξ

σ

P

t

(y,

xi)u(y)dy.

Therefore,

(QP)u(x) =

x

e

i(x

y)·ξ

σ

Q

(x, ξ)σ

P

t

(y,

xi)u(y)dydξ

= a(x, D, y)u(x)

where a(x, ξ, y) = σ

Q

(x, ξ)σ

P

t

(y,

−ξ). Clearly aǫS

m

1

+m

2

(Ω

× Ω) which

implies QPǫψ

m

1

+m

2

(Ω). Moreover

σ

QP

(x, ξ)

X

α

1

α!

α
ξ

D

α
y

Q

(x, ξ)σ

p

t

(y,

−ξ))|

y=x

X

α

1

α!

X

β+ν=α

α

β!ν!

β
ξ

σ

Q

(x, ξ)∂

ν
ξ

D

α
Y

σ

P

t

(x,

−ξ)

X

α,ν

X

δ

(

−1)

|δ|

ν!β!δ!

β
ξ

σ

Q

(x, ξ)D

β+ν+δ

ν+δ

x

ν+δ

ξ

σ

P

(x,

−ξ)

X

β,λ

X

ν+δ=λ

(

−1)

|δ|

ν!β!δ!

1

β!

β
ξ

σ

Q

(x, ξ)D

β+λ

x

λ
ξ

σ

p

(x, ξ)

But

P

v+δ=λ

(

−1)

|δ|

v!δ!

= (x

0

x

0

)

λ

with x

0

= (1, 1, . . . , 1)

=


1 if λ = 0
0 if λ , 0

Therefore σ

QP

(x, ξ)

P

β

1

β!

β
ξ

σ

Q

(x, ξ)D

β

x

σ

p

(x, ξ).

Corollary 4.37 (PRODUCT RULES FOR ψDO

S ).

If Q = q(x, D)ǫψ

m

(λ) and f ǫC

(Ω), then for any positive integer N,

105

q(x, D)( f u) =

P

|α|<N

1

α!

D

α

f [(∂

α
ξ

q)(x, D)u] + T

N

u where T

N

ǫψ

m

N

(Ω).

background image

90

4. Basic Theory of Pseudo Differential Operators

Corollary 4.38. The correspondence p

p(x, D) and P → σ

p

are *

homomorphisms modulo lower order term i.e., if pǫS

m

1

(Ω) and qǫS

m

2

(Ω), then (pq)(x, D)

p(x, D)q(x, D)ǫψ

m

1

+m

2

−1

(Ω) and p(x, D)

− ¯p(x, D)

ǫψ

m

1

−1

and also, if Pǫψ

m

1

(Ω) and Qǫψ

m

2

(Ω), then σ

PQ

− σ

P

σ

Q

ǫS

m

1

+m

2

−1

(Ω), σ

P

− ¯σ + pǫS

m

1

−1

(Ω).

Corollary 4.39. If Pǫψ

m

1

(Ω) and Qǫψ

m

2

(Ω), then [P, Q] = PQ

QP

ǫ ψ

m

1

m

2

−1

(Ω) and σ

[P,Q]

1

i

P

, σ

Q

S

m

1

m

2

−1

(Ω) where

{ f, q} stands

for the Poisson bracket defined by

{ f, q} =

X

f

∂ξ

j

g

x

j

f

x

j

g

∂ξ

j

!

.

Following the philosophy that smoothing operators are negligible

it is a trivial matter to extend these results to non-properly supported
operators.

Suppose Pǫψ

m

(Ω) is not properly supported. By corollary 4.31, we

can write P = P

1

+ s where P

1

is properly supported and S is smoothing.

Then P

t

= P

t
1

+ S where P

t
1

is given by Theorem 4.34 and S

t

is again

smoothing (c.f. Remark 4.33). Likewise for P

. If Q is a properly

supported ψDO, the products PQ and QP are well defined as operators
from C

o

(Ω) to C

(Ω). Again, we have PQ = P

1

Q + S Q and QP =

QP

1

+ Qs; P

1

Q and QP

1

are described by Theorem 4.36, while S Q and

QS are smoothing.

7 A Continuity Theorem for ψ Do on Sobolev

Spaces

We now state and prove a continuity theorem for pseudo defferential

106

operators acting on Sobolev spaces. We are indebted to Dr P.N. Srikanth
for simplifying our original argument.

Theorem 4.40. Suppose P = p(x, D)ǫψ

m

(Ω). Then

i) P : H

s

H

loc

s

m

(Ω) continuous for all sǫR.

background image

7. A Continuity Theorem for ψ Do on Sobolev Spaces

91

ii) If P is properly supported, P : H

loc

s

(Ω)

H

loc

s

m

(Ω) continuously

for all sǫR

Proof. To established (i), we must show that the map u

→ φPu is

bounded from H

s

H

s

m

for any φǫC

o

(Ω). Replacing p(x, ξ) by

φ(x)p(x, ξ) we must show:



(4.41)

if pǫ s

m

(R) and p(x, ξ) = 0 for x outside a compact set, then

p = p(x, D) is bounded from H

s

and H

s

m

for every s in R.

Suppose then that uǫH

s

. Then puǫE

, and from the definition of P

on distributions, we see that

(Pu)

ˆ

(η) =< Pu, e

−2πiη.(.)

>

=

x

e

i

−η).x

p(x, ξ)ˆu(ξ)dxdξ

=

Z

ˆp

1

− ξ, ξ)ˆu(ξ)dξ

and hence, if VǫS ,

< Pu, ¯v > =

Z

(Pu)

ˆ

ˆv =

x

K(η, ξ) f (ξ)¯g(η)dξdη

where

f (ξ) = (1 +

|ξ|

2

)

s/2

ˆu(ξ)

g(η) = (1 +

|η|

2

)

(m

s)/2

ˆv(η)

and K(η, ξ) = ˆp

1

− ξ, ξ)(1 + |ξ|

2

)

s/2

(1 +

|η|

2

)

(s

m)/2

.

107

We wish to estimate K(η, ξ).

For any multi-index α, since

p(., ξ)ǫC

o

, we have

α

ˆp

1

(ζ, ξ)

| = |

Z

(D

α

x

e

−2πix.ζ)

p(x, ξ)dx

|

=

|

Z

e

−2πix.ζ)

(D

α

x

p(x, ξ)dx

|

C

α

(1 +

|ξ|)

m

,

so that for any positive integer N,

| ˆp

1

(ζ, ξ)

| ≤ C

N

(1 +

|ξ|)

m

(1 +

|ζ|)

N

background image

92

4. Basic Theory of Pseudo Differential Operators

and hence

|K(η, ξ)| ≤ C

N

(1 +

|ξ|)

m

s

(1 +

|η|)

s

m

(1 +

|ξ − η|)

N

C

N

(1 +

|ξ|)

|ms|−N

If we take N > n +

|m s|, we see that

Z

|K(η, ξ)|dη ≤ C,

Z

|K(η, ξ)|dξ ≤ C.

Therefore by Theorem 1.1 and the Schwarz inequality,

| < Pu, ¯v > | ≤ C|| f ||

2

||g||

2

= C

||u||

(s)

||v||

(m

s)

.

From this, it follows that

||Pu||

(s

m)

C||u||

(s)

for uǫH

s

which estab-

lishes (4.41) and hence (i).

Now suppose P is properly supported and uǫH

loc

s

(Ω). If φǫC

o

(Ω)

there exists Ω

⋐ Ω such that the values of Pu on supp φ depend only

the values of u on Ω. Thus if we pick φ

ǫC

o

(Ω) with φ

= 1 on Ω

, we

have φPu = φP

u). But φ

uǫH

s

and (4.41) φP is bounded from H

s

to

108

H

s

m

. This establishes (ii) and completes the proof.

8 Elliptic Pseudo Di

fferential Operators

Definition 4.42. Pǫψ

m

(Ω) is said to be elliptic of order m if

P

(x, ξ)

| ≥

C

,

|ξ|

m

for large

|ξ|, for all xǫΩ

, Ω

⋐ Ω.

Definition 4.43. If Pǫψ

m

(Ω), a left (resp. right ) parametrix of P is a

ψDOQ such that QP

I(resp. PQ -I) is smoothing.

Theorem 4.44. If Pǫψ

m

is elliptic of order m and properly supported,

then there exists a properly supported Qǫψ

m

which is a two-sided para-

metrix for P.

Proof. Let P = p(x, D). We will obtain Q = q(x, D) with q

P

j=0

q

j

where the q

j

s are defined recursively. Let ζ(x, ξ) be a C

function such

that ζ(x, ξ) = 1 for large ξ and ζ(x, ξ) = 0 in a neighbourhood of the

background image

8. Elliptic Pseudo Differential Operators

93

zeros of p. Define q

o

(x, ξ) =

ζ(x, ξ)

p(x, ξ)

. Then q

o

ǫ s

m

. Let Q

o

= q

o

(x, D).

We have

σ

Q

p

O

= q

p

0

mod S

−1

= 1 mod S

−1

= 1 + r

1

with r

1

ǫS

−1

.



Let q

1

=

r

1

ζ

p

, Q

1

= q

1

(x, D). Then

σ

(Q

0

+Q

1

)P

= σ

Q

p

O

+ σ

Q

p

1

= 1 + r

1

r

1

ζ

mod S

−2

= 1 + r

2

with r

2

ǫS

−2

.

Let q

1

=

r

2

ζ

p

etc. Having determined q

1

, q

2

, . . . q

j

so that

σ

(Q

0

+Q

1

+

···Q

j

)P

= 1 + r

j+1

, r

j+1

ǫS

−( j+1)

set

q

j+1

=

r

j+1

ζ

p

.

109

Let Q be a properly supported operator with symbol q

P

j=0

q

j

.

Then σ

QP

= 1 mod S

−∞

, i.e., QP

I is smoothing. Thus Q is a left

parametrix. In the same way, we can construct a right parametrix Q

.

Then if S = QP

I and S

= PQ

I, we have

QPQ

= S Q

+ Q

= QS

+ Q.

So Q

Q

= S Q

QS ’ is smoothing. Hence Q is a two sided

parametrix.

Exercise. What happens if P is not properly supported? (cf. the remarks
at the end of §6)
.

The left parametrices are used to prove regularity theorems and the

right parametrices are used to prove existence theorems.

Indeed, we have:

background image

94

4. Basic Theory of Pseudo Differential Operators

Corollary 4.45 (ELLIPTIC REGUIARITY THEOREM). If P is elliptic
of order m
, uǫ

D

(Ω), PuǫH

loc

s

(Ω) implies uǫH

loc

s+m

(Ω). In particular, P is

hypoelliptic.

Proof. Let Q be a left parametrix and set S = QP

I. Then u = QPu

S u. SincePuǫH

loc

s

(Ω) and Q is properly supported, Qǫψ

m

we have

QPuǫH

loc

s+m

(Ω) by Theorem 4.40. Also S uǫC

since S is smoothing.

Hence uǫH

loc

s+m

(Ω).



Theorem 4.46. Every elliptic differential operator is locally solvable.
In fact, if P is elliptic on
Ω, f ǫ

D

(Ω) and x

o

ǫΩ, there exists uǫ

D

(Ω)

110

such that Pu = f in a neighbourhood of x

o

. (Of course, if f ǫC

, then

uǫC

near x

o

, by the previous corollary).

Proof. BY cutting f off away from x

o

, we can assume that f ǫE

and

hence f ǫH

s

for some s. Let Q be a properly supported parametrix for P

and let S = PQ

I. Then S is also properly supported. If Ω

1

⊂⊂ Ω is

a neighbourhood of x

o

then there exists Ω

2

⊂⊂ Ω such that the values

of Su on Ω

1

depend only on the values of u on Ω

2

. Pick φ

1

, φ

2

ǫC

o

(Ω)

such that φ

j

≡ 1 on Ω

j

, j = 1, 2 and set T u = φ

1

S

2

u).



Now observe the following :

i) T u = su on Ω

1

.

ii) T : H

s

C

o

(suup φ

1

) continuously.

From (ii) and the Arzela -Ascoli theorem, it follows that if (u

k

)

is a bounded sequence in H

s

(T u

k

) has a convergent subsequence in

C

o

(suup φ

1

) and hence in H

s

. Therefore, T is compact on H

s

. So

the equation (T + I)u = f can be solved if f is orthogonal (with respect
to the pairing of H

s

and H

s

) to the space N =

{g : (T

+ I)g = 0

}. This

space N is a finite dimensional space of C

functions so we can always

make this happen by modifying f outside a small neighbourhood of x

o

.

Indeed, pick a basis g

1

, g

2

, . . . , g

ν

for N and pick a neighbourhood U

of x

o

so small that g

1

, . . . , g

ν

are linearly independent as functionals on

C

o

(Ω

\U). (Such a U exists; otherwise, by a limiting argument using the

local compactness of N, we could find a nontrivial linear combination of

background image

9. Wavefront Sets

95

g

1

, . . . , g

ν

supported at

{x

o

}, which is absurd). Then we can make f N

by adding to f a function in C

o

(Ω

\U). But then

PQu = (I + S )u = (I + T )u on Ω

1

and (I + T )u = f

in a neighbourhood of x

o

. So Qu solves the problem.

111

9 Wavefront Sets

We now introduce the notion of wavefront sets, which provides a precise
way of describing the singularities of distributions: it specifies not only
the points at which a distribution is not smooth but the directions in
which it is not smooth.

All pseudo differential operators encountered in this section will be

presumed to be properly supported, and by “ψDO” we shall always
mean “properly supported ψDO”.

Definition 4.47.

(i) Let be an open set in R

n

. Then we define

T

o

Ω = Ω

× (R

n

\{0}). (In coordinate - invariant terms, T

o

is

the cotangent bundle of with the zero section removed).

(ii) A set S

T

o

is called conic if (x, ξ)ǫS

⇒ (x, rξ)ǫS, ∀r > 0

(iii) Suppose P = p(x, D)ǫψ

m

(Ω). Then (x

o

, ξ

o

T

o

is said to be non-

characteristic for P if

|p(x, ξ)| ≥ c|ξ|

m

for

|ξ| large and (x, ξ) in

some conic neighbourhood of (x

o

, ξ

o

).

(iv) The characteristic variety of P, denoted by char P, is defined by

char P =

{(x, ξ)ǫT

o

Ω : (x, ξ) is characteristic for p

}.

(v) Let uǫ

D

(Ω). The wavefront set of u, denoted by WF(u) is defined

by

W F(U) =

∩{ charP : Pǫψ(Ω), PuǫC

(Ω)

}.

The restriction Pǫψ

0

(Ω) in the definition of WF(u) is merely a con-

112

venient normalisation. We could allow ψDO of arbitrary order with-
out changing anything, for if Pǫψ

m

(Ω) and Qǫψ

m

(Ω) is elliptic, then

QPǫψ

0

(Ω), char (QP) = char (P)(by Theorem 4.36), and QPuǫC

if

and only if PuǫC

(by Corollaries 4.11 and 4.45).

background image

96

4. Basic Theory of Pseudo Differential Operators

Exercise . When p =

P

|α|=≤m

a

α

(x

α

show that the characteristic variety

of the differential operator p = P(x, D) is

char P =

{(x, ξ) :

X

|α|=m

a

α

(x

α

= 0

}.

The motivation for W F(u) is as follows. If uǫE

, to say that u is not

smooth in the direction ξ

o

should mean that ˆu is not rapidly decreasing

on the through ξ

o

. On the other hand, if PuǫC

o

then (Pu)

ˆ

should be

rapidly decreasing everywhere.

If these conditions both hold, then ξ

o

must be characteristic for P.

Localising these ideas, we arrive at W F(u).

Thus (x

o

ξ

o

W F(u) means roughly that u fails to be C

at x

o

in the

direction ξ

o

. For another interpretation of this statement, see Theorem

4.56 below. For the present, we show that W F(u) is related to sing
supp u as it should be. We denote the projection T

o

Ω onto Ω by π.

Theorem 4.48. For uǫ

D

(Ω), sing supp u = π(WF(u)).

Proof. Suppose x

0

<

sin g supp u. Then there exists φǫC

o

with φ = 1

near x

o

and φuǫC

o

. But multiplication by φ is a ψDO of order 0. Call

113

it P

φ

. Then char P

φ

= π

−1

−1

(0)). Since φ = 1 near x

o

, π

−1

−1

(0))

is disjoint from π

−1

(x

o

). Therefore, π

−1

(x

o

)

WF(u) = φ, i.e., x

o

<

π(W F(u)).



Conversely, suppose x

o

<

π(W F(u)). Then for each ξ with

|ξ| = 1,

there exists Pǫψ

0

(Ω) with PuǫC

(Ω) and (x

o

, ξ) < char P.

Each char P is a closed conic set; so, by the compactness of the unit

sphere, there exists a finite number of ψDO

s say, P

1

, . . . P

N

ǫψ

0

(Ω) with

P

j

uǫC

(Ω) and

N

\

j=1

char P

j

∩ π

−1

(x

o

) = φ. Set P =

N

X

j=1

P

j

P

j

.

Then P is elliptic near x

o

and PuǫC

(Ω). Therefore by Corollary

4.45, u is C

near x

o

, i.e., x

o

<

sin g supp u.

background image

9. Wavefront Sets

97

Definition 4.49. Let P = p(x, D)ǫψ

m

(Ω) and U be an open conic set in

T

o

. We say that P has order

−∞ on U, if, for all closed conic sets

K

U with π(K) compact for every positive integer N and multi-indices

α and β, there exist constants C

αβKN

such that

|D

β

x

D

α
ξ

p(x, ξ)

| ≤ C

αβKN

(1 +

|ξ|)

N

, (x, ξ)ǫK.

The essential support of P is defined to be the smallest closed conic

set outside of which P has order

−∞.

Exercise . If P is a differential operator whose coefficients do not all
vanish at any point of
, show that the essential support of P is T .

Proposition 4.50. Ess. supp PQ

Ess. supp PEss. supp Q.

114

Proof. This follows immediately from the expansion

2

σ

PQ

X

1

α!

α
ξ

σ

P

D

α

x

σ

C

.

Lemma 4.51. Let (x

o

, ξ

o

T

o

and U be any conic neighbourhood of

(x

o

, ξ

o

). Then there exists a Pǫψ

0

(Ω) such that (x

o

, ξ

o

) < char P but Ess.

supp P

U.

Proof. Choose p

o

(ξ)ǫC

with the following properties:

i) p

o

(ξ) = 1 for ξ near ξ

o

, p

o

(ξ) = 0 outside

{ξ : (x

o

ξ)ǫU

} and

ii) p

o

is homogeneous of degree 0 for large

|ξ|.



Then take φǫC

(π(U)) with φ = 1 near x

o

and put p(x, ξ) = p

o

(ξ)

φ(x). This will do the job.

The following theorem and its corollary are refinements of Corol-

lary 4.11 (pseudo local property of ψDO) and Corollary 4.45 (elliptic
regularity theorem).

Theorem 4.52. If Pǫψ

m

(Ω) and uǫD

(Ω), then WF(Pu)

WF(u)∩ Ess.

supp P.

background image

98

4. Basic Theory of Pseudo Differential Operators

Proof. If (x

o

, ξ

o

) < Ess. supp P, by Lemma 4.51, we can find a Qǫψ

0

such that (x

o

, ξ

o

) < char Q and Ess. supp P

∩ Ess. supp Q = φ Then

QPǫψ

−∞

by Proposition 4.50, so that QPǫC

. But this implies that

(x

o

, ξ

o

) < W F(Pu).



Suppose now that (x

o

, ξ

o

) < W F(u). Then there exists Aǫψ

0

with

AuǫC

and (x

o

, ξ

o

) < char A.

Claim . There exist operators B, Cǫψ

0

with (x

o

, ξ

o

) < char B and BP =

115

CA mod ψ

−∞

.

Granted this, BPu = CAu + (C

function) ǫC

which implies that

(x

o

, ξ

o

) < W F(Pu).

To prove the claim, let A

o

be elliptic with σ

A

= σ

A

on a conic

neighbourhood U of (x

o

, ξ

o

). Then Ess. supp(A

A

o

)

U = φ. By

Lemma 4.51, choose B with (x

o

, ξ

o

) < char B and Ess. supp B

U.

Let E

o

be a parametrix for A

o

and set C = BPE

o

.Then CA = BPE

o

A =

BPE

o

(A

A

o

) + BPE

o

A

o

. Since E

o

A

o

= I mod ψ

−∞

and BPE

o

(A

A

o

)

also belongs to ψ

−∞

(by Proposition 4.50 again), CA = BP mod ψ

−∞

.

This completes the proof.

Corollary 4.53. If P is elliptic, then W F(u) = WF(Pu)).

Proof. By the theorem, W F(Pu)

WF(u). If E is a parametrix for

P, then W F(u) = WF(EPu + C

function)

= WF(EPu)

WF(Pu), by the theorem again

Hence

W F(u) = WF(Pu).

2

Theorem 4.56. Let uǫ

D

(D). Then (x

o

, ξ

o

) < W F(u) if and only if there

exists φǫC

o

(R

n

) such that φ(x

o

) = 1 and u)

ˆ

is rapidly decreasing on

a conic neighbourhood of ξ

o

.

Proof. (Sufficiency) If such a φ exists, choose pǫS

0

, p(x, ξ) = p(ξ)

with p(ξ) = 1 near ξ

o

and p(ξ) = 0 outside the region where (φu)

ˆ

is

rapidly decreasing. Then pu)

ˆ

ǫS and hence p(D)(φuS where p(D) =

p(x, D). The operator Pu = φp(D)(φu) is a pseudo differential operator

116

of order 0 with symbol φ(x)

2

p(ξ) modulo S

−1

, so (x

o

, ξ

o

) < char P. This

implies that (x

o

, ξ

o

) < W F(u).



background image

10. Some Further Applications...

99

(Necessity) Suppose (x

o

, ξ

o

) < W F(u). Then there exists a neigh-

bourhood U of x

o

in Ω such that (x

o

, ξ

o

) < W F(u) for all xǫU. Choose

φǫC

o

(U) such that φ(x

o

) = 1. Then (x

o

, ξ

o

) < W Fu) for all xǫR

n

. Let

X

=

{ξ : (x, ξ)ǫWFu) for somex}.

This

P does not contain ξ

o

and is a closed conic set. There exists

p(ξ)ǫS

0

, p = 1 on a conic neighbourhood of ξ

o

and p = 0 on a neigh-

bourhood of

P, say P

o

. Since p(ξ) = 0 on

P

o

, Ess. supp p(D)

∩ (R

n

×

P

o

)

c

and hence Ess. supp p(D)

WFu) = φ. But this gives Ess.

supp(p(Du) = φ. So p(DuǫC

. We now claim that p(DuǫS .

Accepting this, we see that pu)

ˆ

ǫS and in particular, (φu)

ˆ

is rapidly

decreasing near ξ

o

. To prove the claim, observe that

D

β

α

p(ξ)) = 0(1 +

|ξ|)

|α|−|β|

ǫL

2

if

|α| − |β| < −

n

2

.

So if we put K(x) = ˜p(x), x

β

D

α

K(xL

2

if

|α| − |β| < −

n
2

. An application

of Leibniz’s rule and the the Sobolev imbedding theorem shows that

x

β

D

α

K(xL

if

|β| − |α| > −

3n

2

+ 2. Hence K and all its derivatives are

rapidly decreasing at

∞ and since φuǫE

, the same is true of p(Du =

φu

K.

10 Some Further Applications of Pseudo De

fferen-

tial Operators

We conclude our discussion of pseudo differential operators by giving

117

brief and informal descriptions of some further applications. We recall
that an m

th

order ordinary differential equation u

(m)

= F(x, u, u

(1)

, . . .

u

(m

−1)

) can be reduced to the first order system

d

dx

u

1

u

2

..

.

u

m

=

u

2

u

3

..

.

F(x, u

1

, . . . , u

m

)

simply by introducing the derivatives of u of order < m as new vari-
ables, and this reduction is frequently a useful technical device. To do

background image

100

4. Basic Theory of Pseudo Differential Operators

something similar for m

th

order partial differential equations, however,

is more problematical. If one simply introduces all partial derivatives of
u of order < m as new variables and writes down the first order differen-
tial relations they satisfy, one usually obtains more equations than there
are unknowns, because of the equality of the mixed partials. Moreover,
such a reduction usually does not preserve the character of the original
equation. For example take the Laplace equation in two variables :

2

u

x

2

+

2

u

y

2

= f.

putting u

1

= u, u

2

=

u

x

, u

3

=

y

we have a 3

× 3 system given by

u

1

x

u

2

= 0

u

1

y

u

3

= 0

u

2

x

+

u

3

y

= f.

118

The original equation is elliptic. But consider the matrix of the top

order symbols of the 3

× 3 system obtained above. The matrix is given

by

i

ξ

0

0

η

0

0

0

ξ

η

which is not invertible. This means that the first order system is not
elliptic.

There is, however, a method, due to A.P. CALDERON, of reducing

an m

th

order linear partial differential equation to an m

× m system of

first order pseudo differential equations which preserves the character-
istic variety of the equation in a sense which we shall make precise be-
low. In this method, one of the variables is singled out to play a special
role, so we shall suppose that we are working on R

n+1

with coordinates

(x

1

, x

2

, . . . , x

n

, t).

background image

10. Some Further Applications...

101

Let L be a partial differential operator of order m on R

n+1

such that

the coefficient of ∂

m

t

is nowhere vanishing. Dividing throughout by this

coefficient, we can assume that L is of the following form:

L = ∂

m
t

m

−1

X

j=0

A

m

j

(x, t, D

x

)∂

j

t

.

119

Here A

m

j

is a differential operator of order

m j in the x variable

with coefficients depending on t.

We want to reduce the equation Lv = f to a first order system. To

this end we proceed as follows:

Using the familiar operators Λ

s

with symbol (1 +

|ξ|)

s/2

(acting in

the x variables ) we put

u

1

= Λ

m

−1

v

u

2

= Λ

m

−2

t

v

u

3

= Λ

m

−3

2
t

v

..

.

u

m

= ∂

m

−1

t

v.

For j < m, we observe that ∂

t

u

j

= Λu

j+1

. Therefore,

t

u

m

= f +

m

−1

X

j=0

A

m

j

(x, t, D

x

)∂

j

t

v

= f +

m

−1

X

j=0

A

m

j

(x, t, D

x

j

m+1

u

j+1

= f +

m

X

j=1

A

m

j

(x, t, D

x

j

m

u

j

Set B

j

(x, t, D

x

) = A

m

j+1

(x, t, D

x

j

m

. This B

j

is a ΨDO of order

1 in the variable x. Then the equation Lv = f is equivalent to

t

u = Ku + ˜f

background image

102

4. Basic Theory of Pseudo Differential Operators

where u = (u

1

, . . . , u

m

)

t

, ˜

f = (0, 0, . . . , f )

t

. (Here (

· · · )

t

denotes the

120

transpose of the row vector (

· · · )) and K is the matrix given by

K =

0

Λ

0

· · · 0

0

0

Λ

· · · 0

.

.

· · ·

.

.

· · ·

0

0

0

· · · Λ

B

1

B

2

B

3

· · · B

m

This K is a matrix of first order pseudo differential operators in x,

with coefficients depending on t.

Let us now make a little digression to fix on some notations which

will be used in the further development. For a pǫS

m

, a function p

m

(x, ξ)

homogeneous of degree m in ξ is said to be the Principal symbol of

p(x, D) if p

p

m

agrees with an element of S

m

−1

for large

|ξ|. We

remark that not all ΨDO

s have a principal symbol; but most of them

that arise in practice do. Moreover, the principal symbol, if it exists, is
clearly unique.

Examples

(i) If p is a polynomial,

p(x, ξ) =

X

|α|≤m

a

α

(x

α

,

then

p

m

(x, ξ) =

X

|α|=m

a

α

(x

α

.

(ii) If p(x, ξ) = (1 +

|ξ|

2

)

s/2

, p

s

(x, ξ) =

|ξ|

s

.

121

Returning to our discussion, let a

k

(x, t, ξ) be the principal symbol

of A

k

(x, t, d

x

) (Here we are considering A

k

as an operator of order k, so

if it happens to be of lower order, its principal symbol is zero). Then

background image

10. Some Further Applications...

103

the principal symbol of B

j

is b

j

(x, t ξ) = a

m

j+1

(x, t, ξ)

|ξ|

j

m

and the

principal symbol of k is the matrix K

1

K

1

(x, t, ξ) =

0

|ξ|

0

0

· · ·

0

0

0

|ξ|

0

· · ·

0

0

0

0

|ξ| · · ·

0

.

.

.

· · ·

.

.

.

· · ·

.

.

.

· · ·

0

0

0

0

· · · |ξ|

b

1

b

2

b

3

b

4

· · · b

m

The principal symbol of L, on the other hand, is

L

m

(x, t, ξ, τ) = (2πiτ)

m

m

−1

X

j=0

a

m

j

(x, t, ξ)(2πiτ)

j

The characteristic variety of L, i.e., the set of zeros of L

m

, can easily

be read off from the matrix K as follows:

Proposition 4.57. The eigenvalues of K

1

(x, t, ξ) are precisely iτ

1

, . . .,

iτ

m

where τ

1

, τ

2

, . . . , τ

m

are the roots of the polynomial L

m

(x, t, ξ, .).

Proof. The characteristic polynomial p(λ) of K

1

is the determinant of

λ

−|ξ|

0

· · ·

0

0

λ

−|ξ| · · ·

0

..

.

0

0

0

· · ·

−|ξ|

b

1

b

2

b

3

˙˙

λ

b

m

122

Expanding in minors along the last row, we get

p(λ) =

m

−1

X

j=1

(

−1)

m+ j

(

b

j

j

−1

(

−|ξ|)

m

j

+ λ

m

−1

b

m

)

= λ

m

− λ

m

−1

b

m

m

−1

X

j=1

b

j

|ξ|

m

−1

λ

j

−1

background image

104

4. Basic Theory of Pseudo Differential Operators

= λ

m

m

X

j=1

a

m

j+1

λ

j

−1

= λ

m

m

−1

X

j=0

a

m

j

λ

j

= L

m

(x, t, ξ, (2πi)

−1

λλ.

This proves the proposition.



It is also easy to incorporate boundary conditions into this scheme.

Suppose we want to solve the equation Lv = f with the boundary con-
ditions B

j

v

|

t=0

= g

j

, 1, 2, . . . , ν where

B

j

=

m

j

X

k=0

b

k

j

(x, t, D

x

)∂

k
t

, b

k

j

is of order m

j

k and m

j

m − 1.

If we set

B

k

j

= Λ

m

m

j

−1

b

k

−1

j

(x, 0, D

x

k

m

,

φ

j

= Λ

m

m

j

−1

g

j

then with u

j

= Λ

m

j

j

−1

t

v as above and u

0

j

(x) = u

j

(x, 0), the boundary

conditions will become

m

j

+1

X

k=1

B

k

j

u

0
k

= φ

j

, j = 1, 2, . . . ν.

123

This is a system of zeroth order pseudo differential equations.
The above method of reduction is useful, for example, in the follow-

ing two important problems.

(1) Cauchy Problem for Hyperbolic Equations

Lv = f, ∂

j

t

v

|

t=0

= g

j

, j = 0, 1, . . . , m

− 1.

Here the equation is said to be hyperbolic when the eigenvalues of the
matrix K

1

(i.e., the principal symbol of K occurring in the linear system

corresponding to Lv = f ) are purely imaginary.Equivalently, the roots
of L

m

(x, t, ξ, .) are real.

background image

10. Some Further Applications...

105

(2) Elliptic Boundary Value Prob Lems

Lu = f on Ω, B

j

u = g

j

on ∂Ω

where L is elliptic of order 2m and j = 1, 2, . . . m. Here one works
locally near a point x

0

ǫ∂Ω and makes a change of coordinates so that ∂Ω

becomes a hyperplane near x

o

. One then can apply Calderon’s reduction

technique, taking t as the variable normal to ∂Ω. (See Michael Taylor
[3]).

We shall now sketch an example of a somewhat different technique

for applying ψDO

S to elliptic boundary value problems.

Let Ω be a bounded open set of R

n

with C

boundary ∂Ω Consider

the problem

u = 0 in Ω, ∂

X

u + au = g on ∂Ω, aǫC

(∂Ω).

Here X is a real a vector field on the boundary which is nowhere

124

tangent to the boundary and ∂

X

u = grad u.χ.

If ν is the unit outward normal to ∂Ω, by normalising we can assume

that χ = ν + τ where τ is tangent to ∂Ω. We want to use ΨDO to
reduce this to the Dirichlet problem. Pretend for the moment that ∂Ω =

R

n

−1

x

{0} and Ω = {x : x

n

< 0

} so that ∂

ν

=

x

n

.

If ∆u = 0 then ∂

2

ν

u =

n

−1

P

j=1

2

u

x

2

j

=

−∆

b

u, ∆

b

being the Laplacian on

the boundary. So formally ∂

ν

u =

±

(∆

b

u). Indeed, if we compare this

with our discussion of the Poisson kernel in Chapter 2(taking account
of the fact that Ω is now the lower rather than the upper half space), we
see that the equation ∂

ν

u =

(

−∆

b

)u is correct if interpreted in terms of

the Fourier transform in the variables x

= (x

1

, . . . x

n

−1

);

ν

˜u

, x

n

) = 2π

u

, x

n

).

It turns out that something similar works for our original domain Ω.

Namely, if u is smooth on ¯

Ω and ∆u = 0 on Ω then ∂

ν

u

|

∂Ω

= P(u

|

∂Ω

)

where P is a pseudo differential operator of order 1 on ∂Ω which equals

(

−∆

b

) modulo terms of order

≤ 0 where ∆

b

is the Laplace -Beltrami

operator on ∂Ω. In particular, P is elliptic and has real principal symbol.

background image

106

4. Basic Theory of Pseudo Differential Operators

The boundary conditions become Pv +

|∂

t

v + av = g, v = u

|

∂Ω

. This

is a first order pseudo differential equation for ν. It is elliptic (because

τ

has imaginary symbol). So it can be solved modulo smoothing oper-

ators. Since ∂Ω is compact. smoothing operators on ∂Ω are compact, so

125

our boundary equation can be solved provided g is orthogonal to some
finite dimensional space of smooth functions.

Having done this, we have reduced our original problem to the fa-

miliar Dirichlet problem ∆u = 0 in Ω and u

|

∂Ω

= ν which has a unique

solution.

background image

Chapter 5

L

P

and Lipschitz Estimates

OUR AIM IN this chapter is to study how to measure the smoothing

126

properties of pseudo differential operators of non positive order in terms
of various important function spaces. Most of the interesting results
are obtained by considering operators of order

−λ with 0 ≤ λ ≤ n.

Indeed, if PǫΨ

−λ

with λ > n and K(x, y) is the distribution kernel of

P then, by Theorem 4.10 we know that KǫC

j

× Ω) when j < λ −

n. So these operators can be studied by elementary methods. What
is more, when PǫΨ

−λ

, D

α

PǫΨ

−λ+|α|

. So, by a proper choice of α, we

can make 0

≤ λ − |α| ≤ n and then study D

α

P rather than P itself.

Actually, we shall restrict attention to operators of order

−λ where 0 ≤

λ < n. The transitional case λ = n requires a separate treatment to obtain
sharp results; however, for many purpose, it suffices to make the trite
observation that an operator of order

n can be regarded as an operator

of order

n + ǫ. Further, we restrict our attention to ΨDO

s whose

symbols have asymptotic expansions

p(x, ξ)

X

j=0

p

j

(x, ξ)

where P

j

is homogeneous of degree

−λ− j. (These p

j

s no longer belong

to our symbol classes, being singular at ξ = 0, but we can still consider
the operators p

j

(x, D).) By the preceding remarks, it will suffice to con-

sider the operators corresponding to the individual terms in the series

127

107

background image

108

5. L

P

and Lipschitz Estimates

P p

j

whose degrees of homogeneity are between 0 and

n. Thus we are

looking at p(x, D) where p(x, ξ) is C

on R

n

×R

n

\{0} and homogeneous

of degree

−λ in ξ where 0 ≤ λ < n.

Since λ < n, for each x, p(x, .) is locally integrable at the origin

and hence defines a tempered distribution. Denoting the inverse Fourier
transform of this distribution p

2

(x, .) we then have

p(x, D)u(x) =

Z

e

ix

p(x, ξ)ˆu(ξ)dξ

= (p

2

(x, .)

u)(x)

=

Z

p

2

(x, x

y)u(y)dy.

Thus we see that the distribution kernel of p(x, D) is given by K(x, y)

= p

v
2

(x, x

y).

Let us digress a little to make some remarks on homogeneous dis-

tributions.

Definition 5.1. A distribution f ǫS

is said to be homogeneous of degree

µ, i f < f, φ

r

>= r

µ

< f, φ > for all φǫS where φ

r

is the function defined

by φ

r

(x) = r

n

φ(x/r).

Exercise.

1. Show that the above definition agrees with the usual definition of

homogeneity when f is a locally integrable function.

2. Show that if f ǫS

is homogeneous of degree µ, then D

α

f is ho-

mogeneous of degree µ

− |α|.

3. Show that D

α

δ is homogeneous of degree

n − |α| (where δ is the

128

Dirac measure at 0).

Let us now prove a proposition concerning the Fourier transform of

homogeneous distributions. It turns out that the Fourier transform of
such a distribution is also a homogeneous one. Precisely, we have the
following

background image

109

Proposition 5.2. If f is a tempered distribution, homogeneous of degree
µ, then ˆ

f is homogeneous of degree

−µ − n. If f is also C

away from

the origin, then the same is true of ˆ

f .

Proof If f is homogeneous of degree µ, then for φǫS

< ˆ

f , φ

r

> =< f, (φ

r

)ˆ >=< f, r

n

ˆ

φ

1/r

>= r

n−µ

< f, ˆ

φ >

= r

−µ−n

< ˆ

f , φ > .

2

This proves the first assertion. To prove the second assertion, choose

φǫC

0

with φ = 1 in a neighbourhood of the origin and write f = φ f +

(1

− φ) f . Since φ f ǫE

, (φ f )ˆ is C

everywhere. On the other hand,

(1

− φ) f is C

and homogeneous of degree µ for large ξ, and hence

lies in S

µ

(R

n

). Let p(x, ξ) = (1

− φ)(ξ) f (ξ). The distribution kernel of

p(x, D) is given by

K(x, y) =

Z

e

−2πi(yx),ξ

(1

− φ)(ξ) f (ξ)dξ

= ((1

− φ) f )ˆ(y x).

Since K is C

away from the diagonal, we see that ((1

− φ) f )ˆis C

away from the origin. Hence ˆ

f is C

on R

n

\{0}.

Returning to our operators p(x, D) with p(x, ξ) homogeneous of de-

129

gree

−λ in ξ, 0 ≤ λ < n, by the proposition above, we have

p(x, D)u =

Z

K(x, x

y)u(y)dy

where K = P

2

is C

on R

n

× R

n

\{0} and homogeneous of degree λ − n

in the second variable.

Finally, we restrict attention to constant coefficient case, i.e., K(x, x

y) = K(x

y). The essential ideas are already present in this case and

the results we shall obtain can be generalised to the variable coefficient
case in a rather routine fashion. Let us give a name to the objects we are
finally going to study.

Definition 5.3. A tempered distribution K which is homogeneous of de-
gree
λ

n and C

away from the origin is called a kernel of type λ. If

K is a kernel of type λ, the operator T f = K

f is called an operator of

type λ.

background image

110

5. L

P

and Lipschitz Estimates

We now classify the kernels of type λ

≥ 0.

Proposition 5.4. Suppose λ > 0 and f ǫC

(R

n

\{0}) is homogeneous of

degree λ

n in the sense of functions. Then f is locally integrable and

defines a distribution F which is a kernel of type λ.

Conversely, if F is a kernel of type λ and f is the function with which

F agrees on R

n

\{0},then < F, φ >=

R

f φ for every φ in S .

Proof. Since λ > 0, f is locally integrable and of (at most ) polynomial
growth at

∞, so f defines an F in S

. It is easy to check that F is

homogeneous in the sense defined above, i.e.< F, φ

r

>= r

λ

n

< F, φ >,

130

and so F is a kernel of type λ.



For the converse, define G by < G, φ >=< F, φ >

R

f φ, φ

S .

Then G is a distribution supported at 0 and hence we have G =

P c

α

D

α

δ.

Therefore,

< G, φ

r

>=

X

c

α

r

n−|α|

(D

α

φ)(0) = O(r

n

), as r

→ ∞.

On the other hand,

< G, φ

r

>= r

n−λ

< G, φ >

and this is not O(r

n

) as r tends to

∞ unless < G, φ >= 0. Hence G = 0

and this completes the proof.

Suppose that f ǫC

(R

n

\{0}) is homogeneous of degree −n. Then f

is not locally integrable near 0 and so does not define a distribution in a
trivial way. However, let us define

µ

f

=

Z

|x|=1

f (x)dσ(x).

If µ

f

= 0 there is a canonical distribution associated with f which is

called principal value of f, PV( f ), defined by

< PV( f ), φ >= lim

ǫ

→0

Z

|x|>ǫ

f (X)dx.

background image

111

To see that this limit exists, we observe that

Z

ǫ<

|x|<1

f (x)dx = µ

f

1

Z

ǫ

r

−1

dr = 0.

Hence

< PV( f ), φ >= lim

ǫ

→0

Z

ǫ<

|x|<1

f (x)(φ(x)

− φ(0))dx +

Z

|x|≥1

f (x)φ(x)dx,

where the last integrals are absolutely convergent, since

131

|φ(x) − φ(0)| ≤ c|x|,

so that

Z

|x|<1

| f (x)||φ(x) − φ(0)|dx c

Z

|x|<1

|x|

n+1

dx <

∞.

Further, the estimate on φ(x)

− φ(0) depends only on the first deriva-

tives of φ via the mean value theorem, so it is easily verified that the
functional PV( f ) is continuous on S . Finally, we observe that

< PV( f ), φ

r

> = lim

ǫ

→0

Z

|x|>ǫ

f (x)r

n

φ(x/r)dx

= lim

ǫ

→0

Z

|y|>(ǫ/r)

f (y)r

n

φ(y)dy

= r

n

< PV( f ), φ > .

Thus PV( f ) is homogeneous of degree -n and so is a kernel of type

0.

The following theorem gives a sort of converse to the above results.

Theorem 5.5. Suppose F is a kernel of type 0 and f is the function with
which F agrees on
R

n

\{0}. Then µ

f

= 0 and F = PV( f ) + cδ, for some

constant c.

background image

112

5. L

P

and Lipschitz Estimates

Proof. Define the functional G on S by

2

< G, φ >=

Z

|x|≤1

f (x)(φ(x)

− (φ(0))dx +

Z

|x|>1

f (x)φ(x)dx.

By the same argument as above, G is a tempered distribution and

G = F on R

n

/

{0}. Therefore, G F =

P c

α

D

α

δ. Now, < F, φ

r

>= r

n

<

P, φ > so that

< G, φ

r

>

r

n

< G, φ >=< G

Fφ

r

>

r

n

< G

F, φ >= 0(r

n

)

as r tends to

∞. On the other hand, for r > 1,

132

< G, φ

r

>

r

n

< G, φ >=

Z

|x|≤1/r

r

n

f (x)(φ(x)

− φ(0))dx

+

Z

|x|>1/r

r

n

f (x)φ(x)dx

r

n

Z

|x|≤1

f (x)(φ(x)

− φ(0))dx r

n

Z

|x|>1

f (x)φ(x)dx

= r

n

φ(0)

Z

1r

≤|x|≤1

f (x)dx = r

n

φ(0) log rµ

f

for every φ.

This is not 0(r

n

) as r tends to

∞ unless µ

f

= 0. Finally, F

PV( f )

is a kernel of type 0 which is supported at the origin and hence is a
multiple of δ.

To study the boundedness of operators of type λ on L

p

spaces, we

need some concepts from measure theory.

Definition 5.6. Let F be a measurable function defined on R

n

. Then the

distribution function of f is the function δ

f

: (0,

∞) → [0, ∞] given by

δ

f

(t) =

|E

t

|, where E

t

=

{x : | f (x)| > t}.

From the distribution function of f , we can get a large amount of

information regarding f . For example, we have

Z

| f (x)|

p

dx =

Z

0

t

p

dδ

f

(t)

background image

113

(To see this, observe that the Riemann sums for the Stieltjes integral on
the right are approximating sums for the Lebesque integral on the left).

If t

p

δ

f

(t) converges to 0 as t tends to 0 and t tends to

∞, we can

integrate by parts to obtain

R | f (x)|

p

dx = p

R

0

t

p

−1

δ

f

(t)dt.

Using the concept of distribution functions, we will now define weak

133

L

p

spaces.

Definition 5.7. For 1

p < ∞ we define weak L

p

as

{ f : δ

f

(t)

≤ (c/t)

p

for some constant c

}. For f ǫ weak L

p

, the smallest such constant c will

be denoted by [ f ]

p

. Thus, if f ǫweak L

p

, we have

δ

f

(t)

≤ ([ f ]

p

/t)

p

.

For p =

∞, we set weak L

= L

.

Proposition 5.8. CHEBYSHEV ’S INEQUALITY) L

p

weak L

p

and

[ f ]

p

≤ || f ||

p

.

Proof. For f ǫL

p

, if E

t

=

{x : | f (x)| > t},

|| f ||

p
p

=

Z

R

n

| f (x)|

p

dx

Z

E

t

| f (x)|

p

dx

t

p

|E

t

|

i.e.,

|E

t

| ≤ (|| f ||

p

/t)p.

From this, it follows that f ǫ weak L

p

and [ f ]

p

≤ || f ||

p

.



Remark 5.9. It is not true that L

p

= weak L

p

for 1

p < ∞. For

example, f (x) =

|x|

n/p

belongs to weak L

p

but not L

p

.

Remark 5.10. The function f

→ [ f ]

p

satisfies [c f ]

p

=

|c|[ f ]

p

.

But it fails to satisfy triangle inequality and hence is not a norm.

However, since

{| f + g| > t} ⊂ {| f |t/2} ∪ {|g| > t/2},

we have

δ

f +g

(t)

≤ δ

f

(t/2) + δ

g

(t/2)

background image

114

5. L

P

and Lipschitz Estimates

which gives, when f and g are in weak L

p

,

134

δ

f +g

(t)

≤ (2

p

[ f ]

p
p

+ 2

p

[g]

p
p

)/t

p

so that

[ f + g]

p

≤ 2([ f ]

p
p

+ [g]

p
p

)

1/p

≤ 2([ f ]

p

+ [g]

p

).

The functional [.]

p

thus defines a topology on weak L

p

which will

turn it into a (non-locally convex) topological vector space.

Definition 5.11. A linear operator T defined on a space of functions is
said to be of
weak type (p, q), 1

p, q ≤ ∞, if T is a bounded linear

operator from L

p

into weak L

q

, i.e. for every f ǫL

p

, [T f ]

q

c|| f ||

p

for

some constant c independent of f .

We will now state the Marcinkiewicz interpolation theorem and use

it to prove generalisations of Young’s inequality (Theorem 1.3).

Theorem 5.12 (J. Marcinkiewicz). Suppose T is of weak types (p

o

, q

o

)

and (p

1

, q

1

) with 1

p

i

q

i

≤ ∞, p

0

< p

1

, q

0

,

q

1

, i.e., [

Γ f ]

q

i

c

i

|| f ||

p

i

for i = 0, 1. Then, if

1/p

θ

= (1

− θ)/p

0

+ (θ/p

1

), 1/q

θ

= (1

− θ)/q

0

+ (θ/q

1

), 0 < θ < 1,

T is bounded from L

p

θ

to L

q

θ

i.e.,

||T f ||

q

θ

c

θ

|| f ||

p

θ

where the constant

c

θ

depends only on p

0

, q

0

, p

1

, q

1

, c

0

, c

1

and θ.

For the proof of this theorem, see A. Zygmund [5] or E.M.Stein [2].

Theorem 5.13 (GENERAL FROM OF YOUNG’S INEQUALITY). If
(1/p) + (1/q)

− 1 = 1/r, 1 ≤ p, q, r < ∞, Then, for f ǫL

p

, gǫL

q

f

gǫL

r

and we have

|| f g||

r

c

pq

|| f ||

p

||g||

q

.

(In fact, c

pq

≤ 1 for all p, q, although our proof does not yield this

135

estimate).

Proof. Fixing f ǫL

p

, consider the convolution operator g

T

q

= f

g.

We know that for gǫL

1

, T gǫL

p

and

||T g||

p

≤ || f ||

p

||g||

1

. Further if p

is the conjugate of p, then, by H ¨older’s inequality

||T g||

≤ ||g||

p

,

|| f ||

p

for gǫL

p

. Thus we see that the operator T is of weak types (1, p) and

background image

115

(p

,

∞). Therefore, by the Marcinkiewicz interpolation theorem, T maps

L

p

θ

boundedly into L

q

θ

for every 0 < θ < 1 with

(1/p

θ

) = 1

− θ + (θ/p

) = 1

− (θ/p) and 1/q

θ

= (1/p)

− θ/p.



Given r with (1/p) + (1/q)

− 1 = (1/r), set q

θ=r

. Then (1/p

θ

) =

1 + (1/r)

− (1/p) = (1/q). Hence we get the required result.

Theorem 5.14 (WEAK TYPE YOUNG’S INEQUALITY). Let 1

p <

q <

, (1/p) + (1/q) > 1 and (1/p) + (1/q) − 1 = (1/r). Suppose f ǫL

p

and gǫ weak L

q

. Then, we have

a) f

g exists a.e. and is in weak L

r

, and

[ f

g]

r

c

pq

|| f ||

p

[g]

q

b) If p > 1, then f

gǫL

r

and

|| f g||

r

c

pq

|| f ||

p

[g]

q

.

Proof. a) It suffices to assume that

|| f ||

p

= [g]

p

= 1 and to show that

[ f

g]

r

c

pq

. Given α > 0, let

M = (α/2)

p

/(p

q)

(p

/(p

q)

−1/(p

q)

where p

is as usual the conjugate of p. Define

g

1

(x) = g(x), if

|g(x)| > M

= 0 otherwise

and

g

2

(x) = g(x)

g

1

(x).

Then

136

δ

f

g

(α)

≤ δ

f

g

1

(α/2) + δ

f

g

2

(α/2)

and we shall estimate the quantities on the right separately.

By

Holder’s inequality,

|| f g

2

||

≤ || f ||

p

||g

2

||

p

=

||g

2

||

p

.

background image

116

5. L

P

and Lipschitz Estimates

Since (1/q)

− (1/p

) = 1/r > 0, we see that p

q > 0, and hence

||g

2

||

p

p

= p

Z

0

t

p

−1

δ

g

2

(t)dt

= p

M

Z

0

t

p

−1

g

(t)

− δ

g

(M))dt

p

M

Z

0

t

p

−1

t

q

dt

= (p

/(p

q))M

p

q

= (α/2)

p

.

Thus ( f

g

2

)(x) exists at all points and

|| f g

2

||

≤ α/2.

Consequently δ

f

g

2

(α/2) = 0. Next consider

||g

1

||

1

=

Z

0

δ

g

1

(t)dt =

M

Z

0

δ

g

(M)dt +

Z

M

δ

g

(t)dt

M

Z

0

M

q

dt +

Z

M

t

q

dt.

The integral

R

M

t

q

dt converges since q > 1, and we obtain

||g

1

||

1

M

1

q

+ M

1

q

/(q

− 1) = (q/(q − 1))M

1

q

.

Also, by Chebyshev inequality,

δ

f

g

1



α

1



|| f g

1

||

α/2

p

!

p

≤ (2

p

p

)(q/(q

− 1))

p

(α/2)

p

p

q

p(1

q)

(p

/(p

q))

p(1

q)

p

q

= c

pq

α

r

137

Hence (a) is proved.

background image

117

b) The operator f

f g is of weak type (1, q) by (a). Also, if

we choose ¯p > p with (1/ ¯p) + (1/q) > 1 and put (1¯r) = (1/ ¯p) +
(1/q)

− 1 then by (a), f f g is of weak type ( ¯p, ¯r). Therefore, by

Marcinkiewicz, f

f g is bounded from L

p

into L

q

θ

with

1/p

θ

= 1

− θ + (θ/ ¯p), (1/q

θ

) = ((1

− θ)/q) + (θ/¯r).

Putting q

θ

= r we get p

θ

= p. Hence T maps L

p

continuously to L

r

and

|| f g||

r

c

pq

|| f ||

p

[g]

q

.



Corollary 5.15. If T is an operator of type λ, 0 < λ < n, then T is
bounded from L

p

into L

q

, whenever 1 < p < q <

and 1/q = 1/p−λ/n.

Also, T is of weak type (1, n/(n

− λ)).

Proof. If K is a kernel of type λ, we have

|K(x)| = |x|

λ

n

K(x/

|x|)| ≤ c|x|

λ

n

,

which implies that Kǫ weak L

n/(n

−λ)

. Therefore, by Theorem 5.14(b), if

f ǫL

p

and T f = K

f then T f ǫL

q

where

(1/q) = (1/p) + ((n

− λ(/n) − 1 = (1/p) − (λ/n) and ||T f ||

q

c[ f ]

p

c|| f ||

p

.



Also we see that T is of weak type (1, n/(n

− λ)) by Theorem 5.14

(a).

The limiting case of this result with λ = 0 is also true, but this is a

138

much deeper theorem:

Theorem 5.16. (Calderon - Zygmund) Operators of type 0 are bounded
on L

p

, 1 < p <

.

Proof. Let T be an operator of type 0 with kernel K i.e., T f = K

f . To

begin with, we can regard T as a map from c

o

to c

and we shall show

that

||T f ||

p

c

p

|| f ||

p

for f ǫc

o

, 1 < p <

∞, so that T extends uniquely

to a bounded operator on L

p

.



We have K = PV(k) + cδ so that T f = PV(k)

f + c f . Since the

identity operator is continuous, we shall assume that c = 0 and also we
shall identify K with k. The proof now proceed in several steps.

background image

118

5. L

P

and Lipschitz Estimates

Step 1. T is bounded on L

2

. Indeed, we have (T f t = ˆk ˆfk is smooth

away from 0 and homogeneous of degree 0, hence is bounded on R

n

Therefore,

||T f ||

2

=

||(T f )ˆ||

2

≤ ||ˆk||

|| ˆf||

2

=

||ˆk||

|| f ||

2

.

Step 2. Fix a radial function φǫc

o

with φ(x) = 1 for

|x| ≤

1

2

and φ(x) = 0

for

|x| ≥ 1. For ǫ > 0, we define k

ǫ

(x) = k(x)(1

−φ(x/ǫ)) and T

ǫ

f = k

ǫ

f .

Then we claim that T

ǫ

is bounded on L

2

uniformly in ǫ. To see this, we

observe that

(T

ǫ

f )ˆ = ˆfˆk

ǫ

= ˆf(k

kφ(x/ǫ))ˆ

= ˆfˆk

− ˆfk ∗ (φ(x/ǫ))ˆ)

which gives

||T

ǫ

f

||

2

=

||(T

ǫ

f

||

2

≤ || ˆf||

2

||ˆk||

{1 + ǫ

n

Z

| ˆφ(ǫξ)|dξ}

=

|| f ||

2

||ˆk||

{1 + || ˆφ||

1

}

Step 3. T

ǫ

is of weak type (1,1) uniformly in ǫ. The proof of this is more

139

involved and will be given later.

Step 4. By steps 2 and 3, using Marcinkiewicz, we get that T

ǫ

is bounded

on L

p

for 1 < p < 2 uniformly in ǫ.

Step 5. T

ǫ

is bounded on L

p

, for 2 < p <

∞, uniformly in ǫ. Indeed, for

f, gǫC

o

,

Z

(T

ǫ

f ) (x) g(x) dx =

Z

f (x) ( ˜

T

ǫ

g) (x) dx

with ˜

T

ǫ

g = ˜k

ǫ

g, ˜k

ǫ

= k

ǫ

(

x) = k(−x)(1 − φ(x/ǫ)). Since ˜k

ǫ

satisfies

the same conditions as k

ǫ

, we see that ˜

T

ǫ

is bounded on L

q

for 1 < q < 2

uniformly in ǫ. So if (1/p) + (1/q) = 1, 1 < q < 2,

||T

ǫ

f

||

p

=

sup

g ǫ C

o

|

R

(T

ǫ

f )g

|

||g||

q

≤ sup

gǫC

o

|| ˜T

ǫ

g

||

q

||g||

q

|| f ||

p

= c

|| f ||

p

.

background image

119

Step 6. If f ǫC

o

, T

ǫ

f converges to T f in the L

p

norm, 1

p ≤ ∞ as ǫ

tends to 0. Since φ is radial and µ

k

= 0,

Z

|x|=r

φ(x/ǫ)k(x)dσ(x) = c

Z

|x|=r

k(x)dσ(x) = 0;

so, for ǫ

≤ 1, we have

(T

ǫ

f )(x) =

Z

|y|≤1

( f (x

y) − f (x)) k(y)(1 − φ(y/ǫ))dy +

Z

|y|>1

f (x

y) k(y)dy

and hence

(T

ǫ

f )(x)

− (T f )(x) =

Z

|y|≤ǫ

f (x

y) − f (x)k(y)φ(y/ǫ)dy.

140

Now supp(T

ǫ

f

T f ) ⊂ {x : d(x, supp f ) ≤ ǫ} ⊂ A, a fixed compact

set. Since, on compact sets, the uniform norm dominates all L

p

norms,

it suffices to show that (T

ǫ

f

T f ) converges to 0 uniformly on A. But

||T

ǫ

f

f ||

≤ || grad f ||

Z

|y|≤ǫ

c

|y||y|

n

dy

c

|| grad f ||

ǫ

Z

0

r

1

n

r

n

−1

dr = c

′′

ǫ

→ 0, as ǫ → 0.

Step 7. If f ǫL

p

and η > 0, choose gǫC

o

with

||g f ||

p

< η.

Then

||T

ǫ

f

T

δ

f

||

p

≤ ||T

ǫ

( f

g)||

p

+

||T

ǫ

g

T

δ

g

||

p

+

||T

δ

(g

f )||

p

2cη +

||T

ǫ

g

T

δ

g

||

p

.

Since η is arbitrary and

||T

ǫ

g

T

δ

g

||

p

converges to 0 as ǫ, δ tend to 0

by step 6, we see that (T

ǫ

f ) is Cauchy in the L

p

norm. Setting

T f = lim

ǫ

→0

T

ǫ

f,

||T f ||

p

c|| f ||

p

.

Thus the theorem is proved modulo Step 3.
Let us now proceed to the proof of Step 3. First, we need a lemma

which will be used in the proof.

background image

120

5. L

P

and Lipschitz Estimates

Lemma 5.17. Suppose FǫL

1

, F

≥ 0 and α > 0. Then there exists a

sequence (Q

k

) of closed cubes with sides parallel to the coordinate axes

141

and disjoint interiors such that

a) α <

1

|Q

k

|

R

Q

k

F

≤ 2

n

α for all k,

b) If Ω =

S

1

Q

k

, then

|Ω| ≤ (1/α)||F||

1

,

c) F(x)

≤ α for a.e. x < Ω.

Proof. Let r = (

||F||

1

/α)

1/n

and for j = 1, 2, . . . let Q

j

be the collection

of closed cubes of side length r/2

j

and vertices in (r/2

j

)Z

n

. Our se-

quence will be constructed in the following way. Put those cubes QǫQ

1

in the sequence which satisfy α <

1

|Q|

R

Q

F. Then

1

|Q|

Z

Q

F

1

|Q|

||F||

1

= (2

n

α/

||F||

1

)

||F||

1

= 2

n

α

so that the first condition is satisfied.



Put those cubes QǫQ

2

into the sequence which are not contained in

one of the previously accepted cubes and satisfy

α <

1

|Q|

Z

Q

F.

Inductively, put those QǫQ

j

which are not contained in one of the

previously accepted cubes and satisfy α <

1

|Q|

R

Q

F. If QǫQ

j

is in the

sequence and Q

is the cube in Q

j

−1

containing Q then

1

|Q|

Z

Q

F

1

|Q|

Z

Q

, F =

2

n

|Q

|

Z

Q

F.

Since Q

Q

, Q

cannot be in the sequence and hence

142

background image

121

1

|Q

|

Z

Q

F

≤ α which then yields

1

|Q|

Z

Q

F

≤ 2

n

α.

Thus the condition (a) is satisfied. Also,

|Ω| =

X

1

|Q

k

| ≤

1

α

Z

Q

k

F

≤ (1/α)||F||

1

; so (b) follows .

Finally, by the Lebesgue differentiation theorem,

lim

xǫ QǫQ

j

→∞

1

|Q|

Z

Q

F = F(x) for a.e. x.

So, if x < Ω,

1

|Q|

R

Q

F

≤ α for all those Q

s and hence F(x)

≤ α a.e.

on R

n

\Ω.

Coming back to the proof of Step 3, given f ǫL

1

and α > 0, let (Q

k

)

be the sequence of cubes as in the lemma with F =

| f |. We write

f = g +

X

k=1

b

k

with

b

k

(x) =


f (x)

1

|Q

k

|

R

Q

k

f (y)dy, for xǫQ

k

0, otherwise

and

g(x) =


1

|Q

k

|

R

Q

k

f (y)dy, for xǫQ

k

f (x), for x < Ω.

Now me make the following observations :

a) Supp b

k

Q

k

,

R

b

k

= 0 and

(5.18)

X

1

||b

k

||

1

X

1

2

Z

Q

k

| f | ≤ 2

X

1

2

n

α

|Q

k

| ≤ 2

n+1

|| f ||

1

143

background image

122

5. L

P

and Lipschitz Estimates

b)

|g(x)| ≤ 2

n

α for xǫΩ and

|g(x)| ≤ α for a.e. x < Ω.

Therefore,

|g(x)| ≤ 2

n

α a.e. and

||g||

2
2

=

Z

|g|

2

+

Z

R

n

−Ω

|g|

2

(5.19)

≤ (2

n

α)

2

|Ω| + α

Z

R

n

−Ω

|| f ||

≤ (2

2n

+ 1)α

|| f ||

1

.

We put

P

1

b

k

= b so that T

ǫ

f = T

ǫ

g + T

ǫ

b and

|{|T

ǫ

f

| > α}| ≤ |{|T

ǫ

g

| > α/2}| + |{|T

ǫ

b

| > α/2}|.

We shall show that both terms on the right are dominated by

|| f ||

1/α

uniformly in ǫ.

To estimate the first term on the right, we use Chebyshev inequality,

Step 2 and the estimate (5.19) obtaining

|{|T

ǫ

g

| > α/2}| ≤ (2/α)||T

ǫ

g

||

2

)

2

c(||g||

2

/α)

2

c

1

(

|| f ||

1

α).

To estimate the second term, let y

k

be the center of the cube Q

k

and

˜

Q

k

be the cube centred at y

k

but with length side 2

n times that of Q

k

.

We put

S

1

˜

Q

k

= ˜

Ω. Then

| ˜Ω| ≤

X

1

| ˜

Q

k

| =



2

n



n

X

1

|Q

k

| ≤



2

n



n

|| f ||

1

/α = c(

|| f ||

1

/α).

So

|{|T

ǫ

b

| > α/2} ≤ | ˜Ω| + |{|T

ǫ

b

| > α/2}\ ˜Ω|

c(|| f ||

1

/α) +

|{|T

ǫ

b

| > α/2}\ ˜Ω|,

and it suffices to estimate

|{|T

ǫ

| > α/2}/ ˜Ω. Since

R

Q

k

b = 0,

144

T

ǫ

b(x) =

Z

k

ǫ

(x

y)b(y)dy

background image

123

=

X

1

Z

Q

k

k

ǫ

(x

y)b(y)dy

=

X

1

Z

Q

k

(k

ǫ

(x

y) − k

ǫ

(x

y

k

))b(y)dy

Therefore, we have

|{|T

ǫ

b

| > α/2}\ ˜Ω|

≤ (2/α)

Z

R

n

\ ˜Ω

|T

ǫ

b

|dx

≤ (2/α)

X

1

Z

R

n

\ ˜Ω

Z

Q

k

|(k

ǫ

(x

y) − k

ǫ

(x

y

k

))

||b(y)|dy dx

≤ (2/α)

X

1

Z

Q

k

Z

R

n

\ ˜

Q

k

|(k

ǫ

(x

y) − k

ǫ

(x

y

k

))

||b(y)|dx dy

We now claim that

R

R

n

\ ˜

Q

k

|(k

ǫ

(x

y) − k

ǫ

(x

y

k

))

|dx c independent of k and ǫ for yǫQ

k

.

Accepting this claim, by the estimate (5.18), we have

|{|T

ǫ

| > α/2}\ ˜Ω| ≤ (2/α)c

X

1

Z

Q

k

|b(y)|dy ≤ 2

n+2

c(

|| f ||

1

/α).

Thus

|{|T

ǫ

b

| > α}| ≤ c

o

(

|| f ||

1

/α) which completes the proof of Step

3.

Returning to the claim, we observe that if xǫR

n

\ ˜

Q

k

and yǫQ

k

, then

|x y

k

| ≥ 2|y y

k

|.

So, if we set z = x

y

k

, w = y

y

k

we must show that

145

Z

|z|>2|w|

|(k

ǫ

(z

w) − k

ǫ

(z))

|dz c

independent of w and ǫ.

background image

124

5. L

P

and Lipschitz Estimates

First consider the case ǫ = 1. Now k

1

(z) is a C

function which

is homogeneous of degree

n for large z. Therefore, |gradk

1

(z)

| ≤

c

′′

|z|

n−1

. By the mean value theorem,

|k

1

(z

w) − k

1

(z)

| ≤ c

|w| sup

0<t<1

|z tw|

n−1

c

′′

|w||z|

n−1

for

|z| > 2|w|.

So

Z

|z|>2|w|

|(k

1

(z

w) − k

1

(z))

|dz

c

′′

Z

|z|>2|w|

|w||z|

n−1

dz

c

′′′

|w|

Z

2

|w|

r

−2

dr = c.

Now for general ǫ,

k

ǫ

(z) = (1

− φ(z/ǫ))k(z/ǫ)ǫ

n

by the homogeneity of k.

Therefore, if we set z

= ǫ

−1

z and w

= ǫ

−1

w, we see that

Z

|z|>2|w|

|(k

ǫ

(z

w) − k

ǫ

(z))

|dz =

Z

|z

|>2|w

|

|(k

1

(z

w

)

k

1

(z

))

|dz

which is bounded by a constant, by the result for ǫ = 1. Hence the claim
above is proved.

To complete the picture, we should observe that operators of type

146

background image

125

0 are not bounded on L

1

(and hence, by duality, not bounded on L

).

Indeed, if T is an operator of type 0 with kernel k, (T f )ˆ = ˆk ˆf. Since ˆf
is homogeneous of degree 0, it has a discontinuity at 0 (unless k = cδ).
Thus if f ǫL

1

, (T f )ˆ is not continuous at 0 whenever ˆ

f (0) , 0 and this

implies that T f is not in L

1

.

This reflects the fact that if

R

f = ˆf(0) , 0, then T f will not be

integrable near

∞, because k is itself not integrable at ∞. However,

there are also problems with the local integrability of T f caused by the
singularity of k at the origin. In fact, let φǫC

o

be a radial function such

that φ = 1 near 0, and set S f = f

∗ (φk). Then the argument used to

prove Theorem 5.16 shows that S is bounded on L

p

for 1 < p <

∞. In

this case (φk)ˆ = ˆφ

∗ ˆk is in C

, but still S is not bounded on L

1

.

This follows from the following general fact.

Proposition 5.20. If kǫS

and the operator f

k f is bounded on L

1

,

then k is necessarily a finite Borel measure.

Proof. Choose φǫC

o

with

R

φ = 1 and put φ

ǫ

(x) = ǫ

n

φ(x/ǫ). Then

||φ

ǫ

||

1

is independent of ǫ, so

||φ

ǫ

k||

1

c. Therefore there exists a

sequence ǫ

k

tending to 0 such that φ

ǫ

k

k converges to a finite Borel

measure µ in the weak* topology of measures and hence φ

ǫ

k

k converges

to µ in S

also. On the other hand, since (φ

ǫ

k

) is an approximate identity,

φ

ǫ

k

k converges to k in S

.



Hence µ = k and the proposition is proved.

147

It is easy to see that kernels of type 0 are not measures even when

truncated away from the origin, as their total variation in any neighbour-
hood of 0 is infinite. One can also see directly that they do not define
bounded functionals on C

o

.

Exercise . Let k(x) = 1/x on R and let f be a continuous compactly

supported function such that f (x) = (log x)

−1

for 0 < x <

1

2

and f (x) =

0 for x

≤ 0. Show that

lim

x

→0−

( f

PV(k))(x) = ∞.

background image

126

5. L

P

and Lipschitz Estimates

Generalise this to kernels of type 0 on R

n

, n > 1.

This example also shows that operators of type 0 do not map con-

tinuous functions into continuous functions. However, they do preserve
Lipschitz or H ¨older continuity, as we shall now see.

Definition 5.21. For 0 < α < 1, we define

| f |

α

= sup

x,y

| f (x + y) − f (x)|

|y|

α

Λ

α

=

{ f : || f ||

Λ

α

=

de f

|| f ||

+

| f |

α

<

∞}.

Λ

α

is called Lipschitz class of order α.

Remark 5.22. The definition makes perfectly good sense for α = 1
(When α > 1 it is an easy exercise to show that if

| f |

α

<

∞ then f is

constant ). However, we shall not use this definition for α = 1, because
the theorems we wish to prove are false in this case.

We are going to prove, essentially, that operators of type 0 are boun-

ded on Λ

α

. However, if k is a kernel of type 0 and f ǫΛ

α

, the integral

148

defining k

f will usually diverge f need not decay at ∞. Consequently,

we shall work instead with Λ

α

L

p

(1 < p <

∞), concerning which we

have the following useful result.

Proposition 5.23. If f ǫL

p

, 1

p < ∞ and | f |

α

<

, then f ǫΛ

α

and

|| f ||

c(|| f ||

p

+

| f |

α

). Consequently, L

p

∩ Λ

α

is a Banach space with

norm

|| f ||

p

+

| f |

α

.

Proof. Let

A

x

= (

| f (x)|

2

| f |

α

)

1/α

for xǫR

n

.

Then

| f (y)| ≥ | f (x)|/2 for all y such that |x y| ≤ A

x

. So

Z

| f (y)|

p

dy

Z

|xy|≤A

x

| f (y)|

p

dy

≥ | f (x)|

p

2

p

c

A

n

x

= c

′′

| f (x)|

p+(n/α)

| f |

n

α

or

| f (x)|

1+(np)

c

′′′

|| f ||

p

| f |

n/alphap
α

.

2

background image

127

Since this is true for all x, setting θ = n/pα we have

|| f ||

c

′′′

|| f ||

1/(1+θ)

p

| f |

θ/(1+θ)

α

c(|| f ||

p

+

| f |

α

).

Theorem 5.24. Operators of type 0 are bounded on Λ

α

L

p

(0 < α <

1, < 1 < p <

∞).

Proof. Let T : f

K f be an operator of type 0. Since we know that

||T f ||

p

c

p

|| f ||

p

, by proposition 5.23, it will suffice to show

|T f |

α

c

α

| f |

α

for 0 < α < 1, f ǫL

p

∩ Λ

alpha

. As in the proof of Theorem 5.16,

149

we may assume that K = PV(k) and identify K with k.



Given yǫR

n

\{0} and f ǫL

p

∩ Λ

α

define

g(x) =

Z

|z|≤3|y|

k(z) f (x

z)dz

h(x) =

Z

|z|≤3|y|

k(z) f (x

z)dz

so that T f = g + h. Since µ

k

= 0, we have

|g(x)| = |

Z

|z|≤3|y|

k(z)( f (x

z) − f (x))dz|

Z

|z|≤3|y|

c

|z|

n

| f |

α

|z|

α

dz

c

1

| f |

α

|y|

α

.

Since this is true for all x,

|g(x + y) − g(x)| ≤ 2c

1

| f |

α

|y|

α

.

Next

h(x + y) = lim

η

→∞

Z

3

|y|<|z|<η

k(z)( f (x + y

z) − f (x))dz

background image

128

5. L

P

and Lipschitz Estimates

= lim

η

→∞

Z

3

|y|<|z+y|<η

k(z + y)( f (x

z) − f (x))dz.

Therefore,

h(x + y)

h(x) = lim

η

→∞

Z

3

|y|<|z|<η

(k(z + y)

k(z))( f (x z) − f (x))dz + ǫ

1

+ ǫ

2

where ǫ

1

and ǫ

2

are errors coming from difference between the regions

of integration.

ǫ

1

is the error coming from the difference between the regions

|z| < η

and

|z + y| < η.

The symmetric difference between these two regions is contained in

150

the annulus η

− |y| < |z| < η + |y|.

If z is in this region and η

≫ |y| then |z| ≈ |z + y| ≈ η so that

1

| ≤ c

Z

η

−|y|<|z|<η+|y|

n

|| f ||

dz

c

′′

η

n

|| f ||

((η +

|y|)

n

− (η − |y|)

n

)

= 0(η

−1

)

→ 0 as η → ∞.

The term ǫ

2

comes from the symmetric difference of the regions

|z| >

3

|y| and |z + y| > 3|y| which is contained in the annulus 2|y| < |z| < 4|y|.

In this region

|z + y| ≈ |z| ≈ |y|.

Therefore

2

| ≤ c

′′′

Z

2

|y|<|z|<4|y|

|y|

n

| f |

α

dz =

= c

2

| f |

α

|y|

n

((4

|y|)

n

− (2|y|)

n

)

= c

2

| f |

α

|y|

α

.

background image

129

Finally coming to the main term, we have

|k(z + y) − k(z)| ≤ |y| sup

0<t<1

| grad k(z + ty)|

≤ |y| sup

0<t<1

|z + ty|

n−1

c

o

|y||z|

n−1

for

|z| ≥ 3|y|.

Hence, since α < 1 so that

n − 1 + α < −n,

151

|

Z

3

|y|<|z|<η

(k(z + y)

k(z))( f (x z) − f (x))dz|

c

o

Z

3

|y|<|z|<η

|y||z|

n−1

| f |

α

|z|

α

dz

c

o

Z

3

|y|<|z|

|y||z|

n−1+α

| f |

α

dz

c

3

|y|| f |

α

|y|

α

−1

= c

3

|y|

α

| f |

α

.

Therefore,

|T f (x + y) − T f (x)|

|y|

α

c| f |

α

and consequently

|T f |

α

c| f |

α

.

For kernels of positive type, we have the following result.

Theorem 5.25. Suppose 0 < λ < n, 1 < p < n/λ < q < Q where
Q
=

if λ ≤ 1, Q = n/(λ − 1), if λ < 1. Let 1/r = (1/p) − (λ/n) and

α = λ

− (n/q).(Thus r < ∞ and 0 < α < 1). Thus operators of type λ

are bounded from L

p

L

q

into L

r

∩ Λ

α

.

Proof. Let T f = k

f be an operator of type λ. By Corollary 5.15, T is

bounded from L

p

to L

r

, so by proposition 5.23, it is enough to show that

|T f |

α

c|| f ||

q

.

We have (T f )(x) =

Z

f (x

z)k(z)dz

background image

130

5. L

P

and Lipschitz Estimates

(T f )(x + y) =

Z

f (x + y

z)k(z)dz

=

Z

f (x

z)k(z + y)dz,

so that

(T f )(x + y)

T f (x) =

Z

f (x

z)(k(z + y) − k(z))dz

=

Z

|z|≤2|y|

f (x

z)(k(z + y) − k(z))dz

+

Z

|z|>2|y|

f (x

z)(k(z + y) − k(z))dz.



152

If q

is the conjugate exponent of q,

|

Z

|z|≤2|y|

f (x

z)(k(z + y) − k(z))dz|

≤ || f ||

q


Z

|w|≤3|y|

|k(w)

q

dw

1/q

+

Z

|z|≤2|y|

|k(z)

q

dz

1/q


≤ 2|| f ||

q

Z

|z|≤3|y|

|k(z)|

q

dz

1/q

c

1

|| f ||

q

Z

|z|≤3|y|

|z|

n)q

dz

1/q

c

1

|| f ||

q

|y|

n+λ+(n/q

)

= c

1

|| f ||

q

|y|

α

,

by the definition of α

For the second integral, we estimate k(z +y)

k(z) by the mean value

theorem:

|k(z + y) − k(z)| ≤ |y| sup

0<t<1

|gradk(z + ty)|

background image

131

c|y||z|

λ

n−1

for

|z| ≥ 2|y|.

Thus

|

Z

|z|>2|y|

f (x

z)(k(z + y) − k(z))dz|

c|| f ||

q


Z

|z|>2|y|



|y||z|

λ

n−1



q

dz


1/q

c

2

|| f ||

q

|y||y|

λ

n−1+(n/q

)

since (λ

n − 1)q

<

n

= c

2

|| f ||

q

|y|

α

.

Hence

153

|T f (x + y) − T f (x)|

|y|

α

c|| f ||

q

which gives

|T f |

α

c|| f ||

q

.

Remark 5.26. As in the preceding theorem, the reason for taking the
domain of T to be L

p

L

q

instead of just L

q

is that the integral defining

T f will usually diverge when f is merely in L

q

.

Nonetheless, the point of these theorems is that operators of type 0

are in essence, bounded from L

q

to Λ

α

(for appropriate q and α). To

make this precise without losing simplicity, one can observe that opera-
tors of type 0 map Λ

α

E

into Λ

α

while operators of type λ > 0 map

L

q

E

into Λ

α

.

We now introduce spaces of functions whose derivatives upto a cer-

tain order are in L

p

or Λ

α

.

Definition 5.27. Suppose 1

p ≤ ∞ and k is a positive integer. We

define

L

p

k

=

{ f : D

β

f ǫL

p

for0

≤ |β| ≤ k}.

We equip L

p

k

with the norm

|| f ||

k,p

=

P

|β|≤k

||D

β

f

||

p

. (Thus L

2
k

= H

k

in

the notation of Chapter 3).

background image

132

5. L

P

and Lipschitz Estimates

Definition 5.28. Suppose k is a positive integer and k < α < k + 1. We
define

Λ

α

=

{ f : D

β

f ǫΛ

α

k

for0

≤ |β| ≤ k}.

We equip Λ

α

with the norm

|| f ||

Λ

α

=

P

|β|≤k

||D

β

f

||

Λ

α

k

.

Remarks 5.29.
(i) f ǫΛ

α

if and only if D

β

f is bounded and continuous for 0

≤ |β| ≤ k

154

and D

β

f ǫΛ

α

k

for

|β| = k. Indeed, if |β| < k, for |y| ≤ 1,

|D

β

f (x + y)

D

β

f (x)

|

|y|

α

k

|D

β

f (x + y)

D

β

f (x)

|

|y|

c

X

|v|=|β|+1

||D

v

f

||

(by the mean value theorem) and for

|y| > 1,

|D

β

f (x + y)

D

β

f (x)

|

|y|

α

k

≤ 2||D

β

f

||

.

This shows that D

β

f ǫΛ

α

k

for 0

≤ |β| ≤ k i.e., f ǫΛ

α

.

(ii) If k < α < k + 1 then L

p

k

∩ Λ

α

is a Banach space with the norm

X

|β|≤k

(

||D

β

f

||

p

+

|D

β

f

|

α

k

).

This follows from the corresponding fact that L

p

∩ Λ

α

is a Banach space

(Proposition 5.23).

Theorem 5.30. Suppose 0

≤ λ < n, 1 < p < n/λ, 1/r = (1/p) − (λ/n)

and k = 0, 1, 2, . . .. Then we have:

a) Operators of type λ are bounded from L

p

k

into L

r
k

.

b) If λ is an integer, λ = 0, 1, 2, . . . and k < α < k + 1, then operators of

type λ are bounded from L

p

k

∩ Λ

α

to L

r
k

∩ Λ

α+λ

.

Proof. a) This is an easy consequence of Corollary 5.15 and Theorem
5.16, since convolution commutes with differentiation. In the same way
(b) follows Theorem 5.24 when λ = 0.



background image

133

We now proceed by induction on λ. Therefore, assume that λ

≥ 1.

Suppose f ǫL

p

k

∩ Λ

α

. Then ∂

j

( f

k) = f ∗ ∂

j

k and ∂

j

k is a kernel of

155

type λ

− 1. So, if 1/s = (1/p) − (λ − 1)/n, by our induction hypothesis

j

( f

kL

s

k

−1

∩ Λ

α+λ

−1

. Since s < r <

∞, L

s

L

L

r

, so ∂

j

( f

kL

r
k

−1

∩ Λ

α+λ

−1

. Also f

kǫL

r

C

1

, hence in Λ

α+λ

provided it is

bounded. But ∂

j

( f

k) ∈ Λ

α+λ

−1

implies that ∂

j

( f

k) is bounded. By

the mean value theorem, we then have

|( f k)(x + y) − ( f k)(x)|

|y|

c.

This together with f

kǫL

r

implies that f is bounded (by definition

5.21). Hence f

kǫL

r
k

∩ Λ

α+λ

.

The above theorem can be generalised. For example, if 0

≤ λ < n

one can show that operators of type λ map Λ

α

E

in to Λ

α+λ

.

Also generalisations of the L

p

Sobolev spaces L

p

k

can be given for

non-integral values of k. In fact, a theorem due to CALDERON says
that for 1 < p <

∞, f ǫL

p

k

if and only if Λ

k

f ǫL

p

. (Here Λ = (1

− ∆)

1/2

).

Therefore, we can define L

p
s

for any real s by

L

p
s

=

{ f : Λ

s

f ǫL

p

} with the norm || f ||

s,p

=

||Λ

s

f

||

p

.

Then part (b) of the above theorem is still true for 0

≤ λ < n, λ not

necessarily an integer in this case.

Refer to E.M. Stein [2].
We will now prove the Sobolev imbedding theorem for L

p

Sobolev

spaces L

p

k

with positive integral k. This theorem can also be generalised

156

to L

p
s

for s in R.

Theorem 5.31 (SOBOLEV IMBEDDING THEOREM). Suppose 1 <

p <

and k a positive integer. If k < n/p, then L

p

k

L

r

for 1/r =

(1/p)

− (k/n) (Hence also L

p

k+ j

L

r

j

for any j). If k > n/p, and α =

k

n/p is not integer, then L

p

k

⊂ Λ

α

.

Proof. Let N be the fundamental solution of ∆ given by

N(x =)


(2

n)

−1

ω

−1

n

|x|

2

n

for n , 2

(2π)

−1

log

|x| for n = 2.

Then K

j

(x) = ∂

j

N(x) = ω

n

x

j

|x|

n

( true for all n) is a kernel of type

1.



background image

134

5. L

P

and Lipschitz Estimates

Now, if f

L

p

k

E

f = f

∗ δ = f ∗ ∆N = f

n

P

j=1

j

k

j

=

n

P

j=1

(∂

j

f

k

j

).

Suppose k = 1. If 1 < n/p, ∂

j

f ǫL

p

⇒ ∂

j

f

k

j

ǫL

r

where 1/r =

(1/p)

− (1/n), by Theorem 5.25 and hence f ǫL

r

. If 1 > n/p, ∂

j

f

k

j

ǫΛ

1

−(n/p)

by Theorem 5.25 which implies that f ǫΛ

1(n/p)

. Thus the

theorem is true for k = 1.

For k > 1, we proceed by induction. Now

f ǫL

p

k

E

f and ∂

j

f are in L

p

k

−1

E

.

Therefore, if p < n/(k

− 1) and 1/q = (1/p) − (k − 1)/n, we have f ,

j

f ǫL

q

, i.e., f

L

q
1

, while if p > n/(k

− 1), we have f , ∂

j

f ǫΛ

k

−1−(n/p)

,

i.e., f ǫΛ

k

−(n/p)

. In the second case, we are done, and in the first case, we

apply the result for k = 1 to see that f is in the required space. Finally,

157

it is easy to check that if we keep track of the norm inequalities that are
implicit in the above arguments, we obtain

|| f ||

r

c|| f ||

k,p

or

|| f ||

Λ

α

c|| f ||

k,p

,

as appropriate, for f ǫL

p

k

E

. Since L

p

k

E

is clearly dense in L

p

k

, the

desired result follows immediately.

We can summarise these theorems in an elegant way using the fol-

lowing picture.

For

n < α < 0 we define x

α

= L

n/

|α|

and when α > 0, α not an

integer, we define x

α

= Λ

α

.

0

1

2

3

(The small circle represent missing spaces)
In this terminology, we have

Theorem 5.32. Operators of type λ map x

α

E

into x

α+λ

o

≤ λ < n.

Theorem 5.33 (SOBOLEV IMBEDDING THEOREM). If D

β

f ǫ x

α

for

0

≤ |β| ≤ k, then f x

α+k

.

background image

135

We now indicate how to fill the gaps in this picture at α = 0, 1, 2, . . .
For α = 1, we define
Λ

1

=

{ f : f is continuous, bounded and

sup

x,y

| f (x + y) + f (x y) − 2 f (x)|

|y|

<

∞}.

For k = 2, 3, 4, . . . we define Λ

k

=

{ f : D

β

f ǫΛ

1

for

|β| ≤ k − 1}. The

158

sudden jump from first differences in the definition of Λ

α

for α < 1 to

the second differences at α = 1 is less mysterious than it seems at first,
because it can be shown that if 0 < α < 2 then f ǫΛ

α

if and only if f is

bounded, continuous and satisfies

sup

x,y

| f (x + y) + f (x y) − 2 f (x)|

|y|

α

<

∞.

To fill the gap at α = 0 we use the space BMO (“ bounded mean

oscillation”) first introduced by F. JOHN and L. NIRENBERG in 1961,
which is defined as follows.

For f ǫL

1
loc

(R

n

), we denote by m

E

f mean value of f over a measur-

able set E

⊂ R

n

, that is,

m

E

f =

1

|E|

Z

E

f.

Let

Q denote the collection of all cubes in R

n

with sides parallel to

the axes.

Definition 5.34. BMO =

{ f ǫL

1
loc

(R

n

) : sup

QǫQ

m

Q

(

| f m

Q

f

|) < ∞}.

Clearly L

BMO, for, f ǫL

⇒ |m

Q

f

| ≤ || f ||

for every QǫQ

and consequently

m

Q

(

| f m

Q

f

|) ≤ 2|| f ||

.

It can be shown that BMO

L

q
loc

for every q <

∞.

If we define x

α

= Λ

α

for α = 1, 2, . . . and X

0

= BMO then Theorem

5.32 remains valid for all α

ǫ

(

n, ∞), except that the Sobolev imbedding

theorem for α = 0 must be slightly modified as follows:

If D

β

f is in the closure of BMO

E

in BMO for

|β| ≤ k then f ǫΛ

k

.

159

background image

136

5. L

P

and Lipschitz Estimates

WARNING. BMO is not an interpolation space between L

p

and Λ

α

i.e.,

it is not true that if T is a linear operator which is bounded on L

p

for

some p <

and on Λ

α

some α > 0, then T is bounded on BMO.

For proofs of the foregoing assertions, see E.M. Stein [2] and also

the following papers :

(i) E.M Stein and A. Zygmund: Boundedness of translation invari-

ant operators on Holder spaces and L

p

spaces, Ann. Math 85

(11967)337-349, and

(ii) C. Fefferman and E.M. Stein : H

p

spaces of several variables Acta

Math. 129, (1972), 137-193.

As we indicated at the beginning of this chapter, the arguments we

have developed can be extended in a fairly straightforward way to give
estimates for ψDO with variable coefficients. We conclude by sum-
marising the result in

Theorem 5.35. [Let] P = p(x, D) ba a properly supported ψDO of

order

−λ on , where λ ≥ 0 and p

P

j=0

p

j

with p

j

(x, ξ) homogeneous

of degree λ

j for large |ξ|. Then P maps L

p

k

(Ω, loc) into L

p

k

(Ω, loc) for

1 < p <

, and in the terminology of Theorem 5.32 P maps x

α

(Ω, loc)

into x

α+λ

(Ω, loc) for

n < α < ∞.

background image

Bibliography

[1] G. B. FOLLAND: Introduction to partial differential equations

160

Princeton University press, Princeton, N. J., 1975.

[2] E.M. STEIN : Singular integrals and differentiability properties of

functions, Priceton University press, Princeton, N.J., 1970.

[3] M. TAYLOR: (i) Pseudo differential operators,

Lecture Notes in Math # 416, Springer-Verlag, New York, 1974.
(ii) Pseudo differential operators, Princeton University press,
Princeton, N,J., 1981

[4] F. TREVES: Basic linear partial differential equations, Academic

press, New York, 1975.

[5] A. ZYGMUND : Trigonometric series, Cambridge University

press Cambridge, U.K. 1959.

[6] L. HORMANDER : Linear partial differential operators, Springer-

Verlag, New York, 1963.

[7] W. RUDIN : Functional Analysis, McGraw-Hill, New York, 1973.

[8] F. TREVES : Topological Vectors spaces, Distributions and Ker-

nels, Academic press, New York, 1967.

137


Document Outline


Wyszukiwarka

Podobne podstrony:
Whittaker E T On the Partial Differential Equations of Mathematical Physics
Grigoryan V Partial differential equations (draft, UCSB, 2010)(O)(96s) MCde
Partial Differential Equations (excerpt from larger) (2001) WW
Pinchover Y , Rubinstein J An introduction to partial differential equations Extended solutions for
differential equations lecture
Robert Terrell Notes on Differential Equations
Kushner A G Three lectures on contact geometry of Monge Ampere equations (Univ de Los Andes, Bogota,
CALC1 L 11 12 Differenial Equations
Evans L C Introduction To Stochastic Differential Equations
Feynman Lectures on Physics Volume 1 Chapter 04
Complex Numbers and Ordinary Differential Equations 36 pp
Crowley A Lecture on the Philosophy of Magick
Mathematics HL paper 3 series and differential equations 001
Feynman Lectures on Physics Volume 1 Chapter 13
Mathematics HL paper 3 series and differential equations
Feynman Lectures on Physics Volume 1 Chapter 05
Lectures on Language
Feynman Lectures on Physics Volume 1 Chapter 02

więcej podobnych podstron