, ,

In this paper, we study the homogeneous wave equation

\begin{matrix} {\begin{matrix} \partial_{t}^{2} u_{0} - Δ u_{0} & = 0 in R^{n} \times (0, \infty) \\ u_{0} (x, 0) & = φ (x) \\ \partial_{t} u_{0} (x, 0) & = ψ (x) \end{matrix} \end{matrix}

(1.1)

with radially symmetric data in

n

space dimensions. Our goal is to obtain a priori decay estimates for the solution when generic assumptions are imposed on the initial data. Such estimates have been used by several authors in the existence theory for the semilinear wave equation

\partial_{t}^{2} u - Δ u = | u |^{p}

, where

p > 1

. In a later paper, we are going to utilize our results to address the existence of global solutions to the more general nonlinear wave equation with potential

\partial_{t}^{2} u - Δ u = | u |^{p} - V (x) \cdot u

, where

p > 1

and

V (x)

is radially symmetric.

First, consider the solution

u_{0}

of 1.1 when the initial data are such that

\sum_{| α | \leq 3} | \partial_{x}^{α} φ (x) | + \sum_{| α | \leq 2} | \partial_{x}^{α} ψ (x) | \leq ɛ (1 + | x |)^{- k - 1}

for some

ɛ > 0

and

k \geq 0

. In low space dimensions

n = 2, 3

, sharp decay estimates for

u_{0}

were obtained by Kubota [6] and independently by Tsutaya [11, 12] ; see also the earlier work of Asakura [1] . In these papers, the assumption of radial symmetry was not needed to control the

L^{\infty}

-norm of

u_{0}

. To obtain similar decay estimates in higher dimensions

n \geq 4

, however, one has to consider initial data that are either more regular or else radially symmetric. In what follows, we focus on the latter case and allow a generic regularity assumption. Thus, our goal is to study the homogeneous problem

\begin{matrix} {\begin{matrix} \partial_{t}^{2} u_{0} - \partial_{r}^{2} u_{0} - \frac{n - 1}{r} \cdot \partial_{r} u_{0} & = 0 in R_{+}^{2} = (0, \infty)^{2} \\ u_{0} (r, 0) & = φ (r) \\ \partial_{t} u_{0} (r, 0) & = ψ (r) \end{matrix} \end{matrix}

(1.2)

for a fixed integer

n \geq 4

Before we state our main result, however, let us first introduce some notation. Given an integer

n \geq 4

, we define the parameters

a, m

according to the formula

\begin{matrix} (a, m) = {\begin{matrix} (1, \frac{n - 3}{2}) & if n is odd \\ (\frac{1}{2}, \frac{n - 2}{2}) & if n is even. \end{matrix} \end{matrix}

(1.3)

Note that

m \geq 1

whenever

n \geq 4

and that the sum

a + m = (n - 1) / 2

is independent of the parity of

n

. Also, we shall frequently use the bracket notation

〈 λ 〉 = 1 + | λ |

for each

λ \in R

. Now, consider the solution to 1.2 when the initial data are such that

\begin{matrix} \sum_{s = 0}^{l + 1} λ^{s} | φ^{(s)} (λ) | + \sum_{s = 0}^{l} λ^{s + 1} | ψ^{(s)} (λ) | \leq ɛ λ^{l - m} {〈 λ 〉}^{m - l - k} \end{matrix}

(1.4)

for some

ɛ > 0

k \geq 0

and some integer

1 \leq l \leq m

. Note that this condition allows the data to be singular at the origin and automatically holds in the case that

\sum_{s = 0}^{l + 1} {〈 λ 〉}^{s} | φ^{(s)} (λ) | + \sum_{s = 0}^{l} {〈 λ 〉}^{s + 1} | ψ^{(s)} (λ) | \leq ɛ {〈 λ 〉}^{- k} .

Our main result in this paper is the following

Theorem 1.1. Let

n \geq 4

be an integer and define

a, m

by 1.3 . Fix an integer

1 \leq l \leq m

and consider functions

φ \in C^{l + 1} (R_{+})

and

ψ \in C^{l} (R_{+})

which aresubject to 1.4 . Then the homogeneous equation 1.2 admits a unique solution

u_{0} \in C^{l} (R_{+}^{2})

that satisfies the following estimates when

D = (\partial_{r}, \partial_{t})

and

β

is any multi-index with

| β | \leq l

(a) When $0 \leq k < m + a = (n - 1) / 2$ , $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t - r 〉}^{- | β |} {〈 t + r 〉}^{| β | - l + m - k} .$
(b) When $k = m + a = (n - 1) / 2$ , $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t - r 〉}^{- | β |} {〈 t + r 〉}^{| β | - l - a} (1 + ln \frac{〈 t + r 〉}{〈 t - r 〉}) .$
(c) When $m + a < k < 2 (m + a) = n - 1$ , $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t - r 〉}^{m + a - k - | β |} {〈 t + r 〉}^{| β | - l - a} .$
(d) When $k = 2 (m + a) = n - 1$ , $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t - r 〉}^{- m - a - | β |} {〈 t + r 〉}^{| β | - l - a} \cdot (1 + ln 〈 t - r 〉) .$
(e) When $k > 2 (m + a) = n - 1$ , $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t - r 〉}^{- m - a - | β |} {〈 t + r 〉}^{| β | - l - a} .$

Besides, the constant

C_{0}

that appears above depends solely on

k

and

n

When it comes to the decay rates

k < (n - 1) / 2

of part (a), our method leads to slightly sharper conclusions which we establish separately, in Corollary 4.3 . Also, when

n

is odd, our conclusions may be further improved due to the strong form of Huygens' principle. Although we shall not bother to prove this explicitly, one has

u_{0} \in C^{l + 1} (R_{+}^{2})

for odd values of

n

, while the estimate of part (c) holds for any decay rate

k > (n - 1) / 2

whatsoever.

Let us briefly describe the role of Theorem 1.1 in the existence theory for the nonlinear wave equation

\partial_{t}^{2} u - Δ u = | u |^{p}

, where

p > 1

. Since the nonlinearity

| u |^{p}

is a

C^{1}

function, we require the same regularity on the initial data. Then the special case

l = 1

of our theorem asserts that the solution

u_{0}

of the homogeneous equation is

O (r^{1 - m})

r \to 0

. To control this singularity at the origin, one is forced to impose an upper bound on

p

, as in [2, 4, 5] . Our theorem actually refines the estimates of [2, 4, 5] in any space dimension

n \geq 4

, however, one may avoid such precise estimates by utilizing the method-driven upper bound on

p

. Our main contribution here is that, unlike [5] , we are placing no restrictions on the decay rate

k

of the initial data in even space dimensions.

Similarly, let us consider the nonlinear wave equation

\partial_{t}^{2} u - Δ u = | u |^{p}

when

p \geq m

. If we require the same regularity on the initial data, the special case

l = m

of our theorem asserts that the solution

u_{0}

of the homogeneous equation is not singular at the origin. In this case, the upper bound on

p

is redundant, but the estimates of Corollary 4.3 are almost necessary for an iteration argument to go through. Although the slightly weaker estimates of [3] do suffice, those were only obtained for odd space dimensions. The interesting feature in our method is that our approach is the same regardless of the parity of

n

. The plausibility of such an approach is not evident from previous considerations which depended on various representations of the Riemann operator for the wave equation.

Here, some new representations will be established to facilitate this kind of an approach.

Definition 1.2. Given an integer

m \geq 1

, define the

m

th Legendre polynomial by

\begin{matrix} P_{m} (x) = \frac{1}{2^{m} m!} \cdot \frac{d^{m}}{d x^{m}} (x^{2} - 1)^{m} \end{matrix}

(1.5)

and the

m

th Tchebyshev polynomial by

\begin{matrix} T_{m} (x) = \frac{(- 1)^{m}}{(2 m - 1)!!} \cdot \sqrt{1 - x^{2}} \cdot \frac{d^{m}}{d x^{m}} (1 - x^{2})^{m - 1 / 2} . \end{matrix}

(1.6)

Lemma 1.3 (The Riemann operator). Letting

z (λ, r, t)

be the rational function

\begin{matrix} z (λ, r, t) = \frac{λ^{2} + r^{2} - t^{2}}{2 r λ}, \end{matrix}

(1.7)

we define the Riemann operator

L

as follows. When

n

is odd, we set

\begin{matrix} [L f] (r, t) = \frac{1}{2 r^{(n - 1) / 2}} \int_{| t - r |}^{t + r} λ^{(n - 1) / 2} f (λ) \cdot P_{m} (z (λ, r, t)) d λ \end{matrix}

(1.8)

with

m = (n - 3) / 2

. When

n

is even, on the other hand, we set

\begin{matrix} [L f] (r, t) & = \frac{\sqrt{2}}{2 π r^{(n - 1) / 2}} \int_{| t - r |}^{t + r} \int_{z (λ, r, t)}^{1} \frac{λ^{(n - 1) / 2} f (λ)}{\sqrt{σ - z (λ, r, t)}} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ d λ \end{matrix}

\begin{matrix} + \frac{\sqrt{2}}{2 π r^{(n - 1) / 2}} \int_{0}^{max (t - r, 0)} \int_{- 1}^{1} \frac{λ^{(n - 1) / 2} f (λ)}{\sqrt{σ - z (λ, r, t)}} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ d λ \end{matrix}

(1.9)

\begin{matrix}  \end{matrix}

with

m = (n - 2) / 2

. A solution to the Cauchy problem 1.2 is then provided by the formula

\begin{matrix} u_{0} (r, t) = [L ψ] (r, t) + \partial_{t} [L φ] (r, t) . \end{matrix}

(1.10)

When

ψ \in C^{l} (R_{+})

and

φ \in C^{l + 1} (R_{+})

for some integer

l \geq 1

, this solution belongs to

C^{l} (R_{+}^{2})

Proof. Our assertion that 1.10 defines a

C^{l}

function will be established through the proof of Theorem 1.1 , where the derivatives of

u_{0}

will be estimated. Thus, we need only establish the integral representations 1.8 and 1.9 for the Riemann operator. One may derive those using the equivalent representations of Lamb [7] or Rammaha [9] ; see also [10] . For the sake of completeness, however, we shall include their derivation.

The odd-dimensional representation 1.8 is precisely the one that appears in [9] . To prove the even-dimensional representation 1.9 , we manipulate the formula

[L f] (r, t) = \frac{1}{π r^{n / 2}} \int_{0}^{t} \frac{ρ}{\sqrt{t^{2} - ρ^{2}}} \int_{| ρ - r |}^{ρ + r} λ^{(n - 2) / 2} f (λ) \cdot \frac{T_{m} (z (λ, r, ρ))}{\sqrt{1 - z (λ, r, ρ)^{2}}} d λ d ρ

which appears as equation (6b) in [9] . Here,

m = (n - 2) / 2

is an integer and

T_{m}

is the

m

th Tchebyshev polynomial 1.6 . Switching the order of integration, one arrives at

\begin{matrix} L f & = \frac{1}{π r^{n / 2}} \int_{| t - r |}^{t + r} \int_{| λ - r |}^{t} \frac{ρ}{\sqrt{t^{2} - ρ^{2}}} \cdot λ^{(n - 2) / 2} f (λ) \cdot \frac{T_{m} (z (λ, r, ρ))}{\sqrt{1 - z (λ, r, ρ)^{2}}} d ρ d λ \end{matrix}

\begin{matrix} + \frac{1}{π r^{n / 2}} \int_{0}^{max (t - r, 0)} \int_{| λ - r |}^{λ + r} \frac{ρ}{\sqrt{t^{2} - ρ^{2}}} \cdot λ^{(n - 2) / 2} f (λ) \cdot \frac{T_{m} (z (λ, r, ρ))}{\sqrt{1 - z (λ, r, ρ)^{2}}} d ρ d λ . \end{matrix}

\begin{matrix}  \end{matrix}

Note that

z (λ, r, | λ \pm r |) = \mp 1

\partial_{ρ} z (λ, r, ρ) = - ρ / (r λ)

and also

t^{2} - ρ^{2} = 2 r λ \cdot (z (λ, r, ρ) - z (λ, r, t)) .

Once we now use the substitution

σ = z (λ, r, ρ)

in the integrals above, we obtain 1.9 .

In the remaining of this paper, we shall proceed as follows. In section 2 , we establish some basic facts about the various functions that appear in the previous lemma. In section 3 , we combine these facts to estimate the Riemann operator and its derivatives. Finally, section 4 is devoted to the proof of Theorem 1.1 and its improved version, Corollary 4.3 .

2 Basic Estimates and Facts

Our main goal in this section is to collect a few basic facts about the Riemann operator of Lemma 1.3 . Its integral representation involves the rational function 1.7 regardless of the parity of

n

, so we intend to focus on this function first. In our next lemma, we estimate its derivatives using a rather painstaking approach. However, we do need such an approach in order to gain certain cancellations in our subsequent treatise of the Riemann operator.

Lemma 2.1. With

(r, t) \in R_{+}^{2}

arbitrary and

z = z (λ, r, t)

as in 1.7 , the sharp estimate

\begin{matrix} | \partial_{λ}^{i} z | \leq C (i) \cdot \frac{λ^{2} + | t^{2} - r^{2} |}{r λ^{i + 1}}, λ \geq 0 \end{matrix}

(2.1)

holds for each integer

i \geq 0

, and the general estimate

\begin{matrix} | \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} z | \leq C (i, j, k) \cdot λ^{- 1 - i} r^{- 1 - j} \cdot (t + r)^{2 - k}, 0 \leq λ \leq t + r \end{matrix}

(2.2)

holds for all integers

i, j, k \geq 0

. In the case that

t \geq r

, one also has the auxiliary estimate

\begin{matrix} | D^{β} z | \leq C (β) \cdot \frac{max (r^{2} + | t^{2} - λ^{2} |, r t)}{λ r^{| β | + 1}}, 0 \leq λ \leq t + r \end{matrix}

(2.3)

with

D = (\partial_{r}, \partial_{t})

and

β

an arbitrary multi-index.

Proof. Our first assertion 2.1 follows trivially once we explicitly compute

\begin{matrix} z = \frac{λ^{2} - t^{2} + r^{2}}{2 r λ}; \partial_{λ} z = \frac{λ^{2} + t^{2} - r^{2}}{2 r λ^{2}}; \partial_{λ}^{i} z = C (i) \cdot \frac{t^{2} - r^{2}}{r λ^{i + 1}} if i \geq 2 . \end{matrix}

(2.4)

Let us now turn to our second assertion 2.2 . According to Leibniz' rule, we have

| \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} z | \leq C \sum_{i_{1} = 0}^{i} \sum_{j_{1} = 0}^{j} λ^{- 1 - i_{1}} r^{- 1 - j_{1}} \cdot | \partial_{λ}^{i - i_{1}} \partial_{r}^{j - j_{1}} \partial_{t}^{k} (λ^{2} - t^{2} + r^{2}) | .

Since

λ, r, t \leq t + r

by assumption, we may thus deduce the desired estimate

\begin{matrix} | \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} z | & \leq C \sum_{i_{1} = 0}^{i} \sum_{j_{1} = 0}^{j} λ^{- 1 - i_{1}} r^{- 1 - j_{1}} \cdot (t + r)^{2 - (i - i_{1}) - (j - j_{1}) - k} \end{matrix}

\begin{matrix} \leq C λ^{- 1 - i} r^{- 1 - j} \cdot (t + r)^{2 - k} . \end{matrix}

\begin{matrix}  \end{matrix}

Finally, we prove our last assertion 2.3 . In the case that

β = (j, 0)

, the inequality

| D^{β} z | = | \partial_{r}^{j} z | \leq C (j) \cdot \frac{r^{2} + | t^{2} - λ^{2} |}{λ r^{j + 1}}, λ \geq 0

follows by 2.1 because

z

remains unchanged when the roles of

λ

and

r

are interchanged.

To settle the remaining case

β = (j, k)

with

k \geq 1

, we resort to our general estimate 2.2 .

Since

t \geq r

by assumption,

t + r

is equivalent to

t

, so we get

| D^{β} z | = | \partial_{r}^{j} \partial_{t}^{k} z | \leq C λ^{- 1} r^{- 1 - j} \cdot t^{2 - k} \leq C λ^{- 1} r^{- j - k} \cdot t

because

t \geq r

and

k \geq 1

. This also completes the proof of our last assertion 2.3 .

Corollary 2.2. With

D_{*} = (\partial_{λ}, \partial_{r}, \partial_{t})

and

α

any multi-index, one has

| D_{*}^{α} z (λ, r, t) | \leq C (α) \cdot λ^{- | α |}

whenever

0 < t \leq 2 r

and

| t - r | \leq λ \leq t + r

Proof. In the case that $α = (i, j, k)$ with $j + k \geq 1$ , our general estimate 2.2 gives $| D_{*}^{α} z (λ, r, t) | \leq C λ^{- 1 - i} r^{- 1 - j} \cdot (t + r)^{2 - k} \leq C λ^{- 1 - i} r^{1 - j - k} \leq C λ^{- i - j - k}$ because $r \leq t + r \leq 3 r$ and $λ \leq t + r \leq 3 r$ by assumption.
In the case that $α = (i, 0, 0)$ , on the other hand, our sharp estimate 2.1 gives $| D_{*}^{α} z (λ, r, t) | = | \partial_{λ}^{i} z (λ, r, t) | \leq C (i) \cdot \frac{λ^{2} + | t^{2} - r^{2} |}{r λ^{i + 1}} \leq C (i) \cdot \frac{λ + t + r}{r λ^{i}} \leq C λ^{- i}$ since $| t - r | \leq λ$ and $λ + t + r \leq 2 (t + r) \leq 6 r$ . In either case then, the result follows.

In our next corollary, we concern ourselves with the derivatives of

1 / (\partial_{λ} z)

. This function will arise as soon as we integrate by parts the integrals of Lemma 1.3 , namely the ones that appear in the explicit representation of the Riemann operator. To estimate its derivatives in a rather precise manner, we first recall 2.1 and 2.4 , according to which

(\partial_{λ} z)^{- 1} = \frac{2 r λ^{2}}{λ^{2} + t^{2} - r^{2}}, | \partial_{λ}^{i} z | \leq C (i) \cdot \frac{λ^{2} + | t^{2} - r^{2} |}{r λ^{i + 1}} .

Keeping such expressions intact, one does not have to estimate them separately while dealing with their product. This is also the main idea in the following proof, where special emphasis is laid on the derivatives with respect to

λ

Corollary 2.3. With

z = z (λ, r, t)

as in 1.7 , one has

\begin{matrix} | \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} (\partial_{λ} z)^{- 1} | \leq C (i, j, k) \cdot \frac{λ^{2 - i} r^{1 - j} t^{j - 1}}{(t - r)^{1 + j + k}} \end{matrix}

(2.5)

whenever

0 < r \leq t

and

0 \leq λ \leq t + r

Proof. Let us set

D = (\partial_{r}, \partial_{t})

and

β = (j, k)

for ease of notation. By repeated applications of the chain rule, one obtains an identity of the form

\begin{matrix} D^{β} (\partial_{λ} z)^{- 1} = \sum_{l = 1}^{| β |} (\partial_{λ} z)^{- 1 - l} \sum_{β_{1} + \dots + β_{l} = β} C (β_{1}, \dots, β_{l}, l) \cdot^{l} \prod_{q = 1} (D^{β_{q}} \partial_{λ} z), \end{matrix}

(2.6)

where each

β_{q}

is a multi-index. Further differentiating with respect to

λ

, one then obtains

\begin{matrix} | \partial_{λ}^{i} D^{β} (\partial_{λ} z)^{- 1} | \leq C \sum_{l = 1}^{| β |} \sum_{s = 0}^{i} | \partial_{λ}^{i - s} (\partial_{λ} z)^{- 1 - l} | \sum_{β_{1} + \dots + β_{l} = β} | {\partial_{λ}^{s}}^{l} \prod_{q = 1} (D^{β_{q}} \partial_{λ} z) | \end{matrix}

(2.7)

by Leibniz' rule. Once we write

β_{q} = (j_{q}, k_{q})

for each

q

, the innermost sum becomes

\sum_{β_{1} + \dots + β_{l} = β} | {\partial_{λ}^{s}}^{l} \prod_{q = 1} (D^{β_{q}} \partial_{λ} z) | \leq \sum_{j_{1} + \dots + j_{l} = j} \sum_{k_{1} + \dots + k_{l} = k} \sum_{s_{1} + \dots + s_{l} = s} C^{l} \prod_{q = 1} | \partial_{λ}^{s_{q} + 1} \partial_{r}^{j_{q}} \partial_{t}^{k_{q}} z | .

Meanwhile, our general estimate 2.2 ensures that

| \partial_{λ}^{s_{q} + 1} \partial_{r}^{j_{q}} \partial_{t}^{k_{q}} z | \leq C λ^{- 2 - s_{q}} r^{- 1 - j_{q}} \cdot (t + r)^{2 - k_{q}}, 0 \leq λ \leq t + r

so we may combine the last two equations to arrive at

\sum_{β_{1} + \dots + β_{l} = β} | {\partial_{λ}^{s}}^{l} \prod_{q = 1} (D^{β_{q}} \partial_{λ} z) | \leq C λ^{- 2 l - s} r^{- l - j} \cdot (t + r)^{2 l - k} .

Next, we insert this fact in 2.7 . Since

D = (\partial_{r}, \partial_{t})

and

β = (j, k)

by above, we find

\begin{matrix} | \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} (\partial_{λ} z)^{- 1} | \leq C \sum_{l = 1}^{j + k} \sum_{s = 0}^{i} | \partial_{λ}^{i - s} (\partial_{λ} z)^{- 1 - l} | \cdot λ^{- 2 l - s} r^{- l - j} \cdot (t + r)^{2 l - k} . \end{matrix}

(2.8)

To handle the latter derivatives, we will employ our sharp estimate 2.1 instead of 2.2 . First, we use repeated applications of the chain rule to establish the inequality

| \partial_{λ}^{i - s} (\partial_{λ} z)^{- 1 - l} | \leq C \sum_{m = 1}^{i - s} | \partial_{λ} z |^{- 1 - l - m} {\sum_{a_{1} + \dots + a_{m} = i - s}}^{m} \prod_{q = 1} | \partial_{λ}^{a_{q} + 1} z |

in analogy with 2.6 . When it comes to the product, 2.1 and 2.4 combine to give

^{m} \prod_{q = 1} | \partial_{λ}^{a_{q} + 1} z | \leq {C_{1}}^{m} \prod_{q = 1} (\frac{λ^{2} + t^{2} - r^{2}}{r λ^{a_{q} + 2}}) = 2^{m} C_{1} | \partial_{λ} z |^{m} \cdot^{m} \prod_{q = 1} λ^{- a_{q}},

so we get

| \partial_{λ}^{i - s} (\partial_{λ} z)^{- 1 - l} | \leq C | \partial_{λ} z |^{- 1 - l} \cdot λ^{s - i} .

Inserting this fact in 2.8 , we now arrive at

| \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} (\partial_{λ} z)^{- 1} | \leq C \sum_{l = 1}^{j + k} | \partial_{λ} z |^{- 1 - l} \cdot λ^{- 2 l - i} r^{- l - j} \cdot (t + r)^{2 l - k} .

According to our computation 2.4 , we also have

| \partial_{λ} z |^{- 1 - l} = {(\frac{2 λ^{2} r}{λ^{2} + t^{2} - r^{2}})}^{l + 1} \leq C λ^{2 l + 2} r^{l + 1} \cdot (t + r)^{- l - 1} (t - r)^{- l - 1},

whence

| \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} (\partial_{λ} z)^{- 1} | \leq C \sum_{l = 1}^{j + k} λ^{2 - i} r^{1 - j} \cdot (t + r)^{l - j - k + j - 1} (t - r)^{- l - 1} .

Since

l \leq j + k

within the last sum, this trivially gives

| \partial_{λ}^{i} \partial_{r}^{j} \partial_{t}^{k} (\partial_{λ} z)^{- 1} | \leq C λ^{2 - i} r^{1 - j} \cdot (t + r)^{j - 1} (t - r)^{- 1 - j - k} .

Besides,

t + r

is equivalent to

t

whenever

t \geq r

, so the desired estimate 2.5 follows.

The last fact we need to treat the Riemann operator in odd dimensions is also the most crucial one and appears in our next lemma. Here, we use the integral representation 1.8 to derive some new representations 2.9 . The latter are only valid in the interior of the light cone, however, they are less singular at

r = 0

Lemma 2.4. Let

n \geq 5

be an odd integer and

a, m

be as in 1.3 . Suppose that

f_{0} \in C^{l} (R_{+})

and

f_{1} \in C^{l + 1} (R_{+})

for some integer

1 \leq l \leq m

. When

t \geq r

, the Riemann operator 1.8 is then subject to an identity of the form

\begin{matrix} \partial_{t}^{i} [L f_{i}] (r, t) = \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} [H_{i j} f_{i}] (λ, r, t) \cdot r^{- m - a} P_{j m} (z (λ, r, t)) d λ \end{matrix}

(2.9)

for

i = 0, 1

and all integers

i \leq j \leq l

. Here,

P_{j m}

denotes a polynomial 2.11 which vanishes with order

j

at each of

\pm 1

, while

H_{i j}

denotes a linear operator which acts on functions of

λ

and is defined by either 2.12 or 2.14 . If we let

D = (\partial_{r}, \partial_{t})

and

β

be any multi-index, then we also have the estimate

\begin{matrix} | D^{β} [H_{i j} f] | \leq C (β, j) \cdot \frac{λ^{j} r^{j - | β |} t^{| β | - j}}{(t - r)^{i + j + | β |}} \cdot \sum_{s = 0}^{i + j} λ^{m + a + s} | f^{(s)} (λ) |, 0 \leq λ \leq t + r \end{matrix}

(2.10)

for

i = 0, 1

and any integer

j \geq 0

, provided that

t \geq r

as above.

Proof. First, we shall use induction on

j

to prove 2.9 for the case

i = 0

and

0 \leq j \leq l

Explicitly, the polynomials

P_{j m}

of interest are given by the formula

\begin{matrix} P_{j m} (x) \equiv \frac{1}{2^{m} m!} \cdot \frac{d^{m - j}}{d x^{m - j}} (x^{2} - 1)^{m}, 0 \leq j \leq m . \end{matrix}

(2.11)

Note that

P_{0 m}

is merely the Legendre polynomial 1.5 . Besides,

P_{j m}

vanishes with order

j

at each of

\pm 1

, and we also have

P_{j m} = P_{j + 1, m}^{'}

for each

j

. As for the linear operators

H_{0 j}

we are going to need, those are given by the formula

\begin{matrix} [H_{0 j} f] (λ, r, t) = {[\frac{\partial}{\partial λ} \frac{1}{\partial_{λ} z (λ, r, t)}]}^{j} (λ^{m + a} f (λ)), \end{matrix}

(2.12)

where

z (λ, r, t)

is the rational function 1.7 . We may introduce them for any integer

j \leq l

and any function

f \in C^{l} (R_{+})

. Let us now proceed to the induction argument. Since

t \geq r

by assumption, the Riemann operator 1.8 takes the form

[L f_{0}] (r, t) = \frac{1}{2} \int_{t - r}^{t + r} λ^{(n - 1) / 2} f_{0} (λ) \cdot r^{- (n - 1) / 2} P_{0 m} (z (λ, r, t)) d λ .

Moreover,

(n - 1) / 2 = m + a

by our definition 1.3 , so the desired identity

\begin{matrix} [L f_{0}] (r, t) = \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} [H_{0 j} f_{0}] (λ, r, t) \cdot r^{- m - a} P_{j m} (z (λ, r, t)) d λ \end{matrix}

(2.13)

does hold when

j = 0

. Suppose it holds for some

j \geq 0

. As we have already remarked, the polynomials 2.11 are such that

P_{j m} (z (λ, r, t)) = P_{j + 1, m}^{'} (z (λ, r, t)) = \frac{1}{\partial_{λ} z (λ, r, t)} \cdot \partial_{λ} P_{j + 1, m} (z (λ, r, t)) .

Integrating 2.13 by parts, we thus get to replace the leftmost factor in the integrand by

- \frac{\partial}{\partial λ} (\frac{1}{\partial_{λ} z (λ, r, t)} \cdot [H_{0 j} f_{0}] (λ, r, t)) = - [H_{0, j + 1} f_{0}] (λ, r, t)

and the polynomial

P_{j m}

P_{j + 1, m}

. This allows us to finish the inductive proof of 2.13 as long as no boundary terms arise in the process. In fact,

P_{j + 1, m}

vanishes at each of

\pm 1

and our definition 1.7 gives

z (t \pm r, r, t) = \pm 1

, hence no boundary terms arise, indeed.

Next, we establish 2.9 for the case

i = 1

and

1 \leq j \leq l

. The linear operators

H_{1 j}

we are going to need are given by the formula

\begin{matrix} [H_{1 j} f] (λ, r, t) = \partial_{t} [H_{0 j} f] (λ, r, t) - \partial_{λ} ((\partial_{λ} z)^{- 1} \cdot \partial_{t} z \cdot [H_{0 j} f] (λ, r, t)) . \end{matrix}

(2.14)

We may introduce them for any integer

j \leq l

and any function

f \in C^{l + 1} (R_{+})

. Let us now employ the identity 2.13 we just proved to write

[L f_{1}] (r, t) = \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} [H_{0 j} f_{1}] (λ, r, t) \cdot r^{- m - a} P_{j m} (z (λ, r, t)) d λ .

Since

j \geq 1

, we have

P_{j m} (z (t \pm r, r, t)) = P_{j m} (\pm 1) = 0

by above, so we find that

\begin{matrix} \partial_{t} [L f_{1}] (r, t) & = \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} \partial_{t} [H_{0 j} f_{1}] (λ, r, t) \cdot r^{- m - a} P_{j m} (z (λ, r, t)) d λ \end{matrix}

\begin{matrix} + \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} \partial_{t} z \cdot [H_{0 j} f_{1}] (λ, r, t) \cdot r^{- m - a} P_{j m}^{'} (z (λ, r, t)) d λ . \end{matrix}

\begin{matrix}  \end{matrix}

Recalling our definition 2.14 , we may then integrate the latter integral by parts to get

\partial_{t} [L f_{1}] (r, t) = \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} [H_{1 j} f_{1}] (λ, r, t) \cdot r^{- m - a} P_{j m} (z (λ, r, t)) d λ .

This is precisely the desired identity 2.9 for the case

i = 1

and

1 \leq j \leq l

Finally, we turn our attention to the estimate 2.10 , which asserts that

\begin{matrix} | \partial_{r}^{j_{0}} \partial_{t}^{k_{0}} [H_{i j} f] | \leq C (j_{0}, k_{0}, j) \cdot \frac{λ^{j} r^{j - j_{0} - k_{0}} t^{j_{0} + k_{0} - j}}{(t - r)^{i + j + j_{0} + k_{0}}} \cdot \sum_{s = 0}^{i + j} λ^{m + a + s} | f^{(s)} (λ) | \end{matrix}

(2.15)

for

i = 0, 1

and each

0 \leq λ \leq t + r

. Allowing the remaining indices to be all arbitrary, our plan is to establish the more general estimate

\begin{matrix} | \partial_{λ}^{i_{0}} \partial_{r}^{j_{0}} \partial_{t}^{k_{0}} [H_{i j} f] | \leq C (i_{0}, j_{0}, k_{0}, j) \cdot \frac{λ^{j - i_{0}} r^{j - j_{0}} t^{j_{0} - j}}{(t - r)^{i + j + j_{0} + k_{0}}} \cdot \sum_{s = 0}^{i + i_{0} + j} λ^{m + a + s} | f^{(s)} (λ) | \end{matrix}

(2.16)

for

i = 0, 1

and each

0 \leq λ \leq t + r

. When

i_{0} = 0

, this actually improves 2.15 by an extra factor

r^{k_{0}} t^{- k_{0}}

, which is less than

1

since

r \leq t

by assumption. In what follows, we may thus focus on the derivation of 2.16 , instead.

To prove 2.16 for the case

i = 0

, we need to check that

\begin{matrix} | \partial_{λ}^{i_{0}} \partial_{r}^{j_{0}} \partial_{t}^{k_{0}} [H_{0 j} f] | \leq C (i_{0}, j_{0}, k_{0}, j) \cdot \frac{λ^{j - i_{0}} r^{j - j_{0}} t^{j_{0} - j}}{(t - r)^{j + j_{0} + k_{0}}} \cdot \sum_{s = 0}^{i_{0} + j} λ^{m + a + s} | f^{(s)} (λ) | \end{matrix}

(2.17)

for each

0 \leq λ \leq t + r

. Let us proceed using induction on

j

. When

j = 0

and

j_{0} + k_{0} \geq 1

, our task is trivial since

H_{00} f = λ^{m + a} f (λ)

depends only on

λ

. When

j = 0

and

j_{0} = k_{0} = 0

, on the other hand, an application of Leibniz' rule gives

| \partial_{λ}^{i_{0}} [H_{00} f] | \leq C (i_{0}) \cdot \sum_{s = 0}^{i_{0}} λ^{m + a + s - i_{0}} | f^{(s)} (λ) | .

This proves 2.17 when

j = 0

, so suppose the same estimate holds for some

j \geq 0

. In view of our definition 2.12 , we may then write

[H_{0, j + 1} f] (λ, r, t) = \partial_{λ} ((\partial_{λ} z)^{- 1} \cdot [H_{0 j} f] (λ, r, t)) .

The derivatives of

(\partial_{λ} z)^{- 1}

were treated in Corollary 2.3 , while those of

H_{0 j} f

are subject to our induction hypothesis 2.17 . Resorting to Leibniz' rule, one may easily use these facts to finish the inductive proof of 2.17 , so we shall omit the details.

To establish 2.16 for the case

i = 1

, we need to check that

\begin{matrix} | \partial_{λ}^{i_{0}} \partial_{r}^{j_{0}} \partial_{t}^{k_{0}} [H_{1 j} f] | \leq C (i_{0}, j_{0}, k_{0}, j) \cdot \frac{λ^{j - i_{0}} r^{j - j_{0}} t^{j_{0} - j}}{(t - r)^{1 + j + j_{0} + k_{0}}} \cdot \sum_{s = 0}^{i_{0} + 1 + j} λ^{m + a + s} | f^{(s)} (λ) | \end{matrix}

(2.18)

for each

0 \leq λ \leq t + r

. Here, we may proceed directly starting with the definition

[H_{1 j} f] (λ, r, t) = \partial_{t} [H_{0 j} f] (λ, r, t) - \partial_{λ} ((\partial_{λ} z)^{- 1} \cdot \partial_{t} z \cdot [H_{0 j} f] (λ, r, t))

we introduced in 2.14 . The derivatives of

H_{0 j} f

are subject to 2.17 , those of

(\partial_{λ} z)^{- 1}

were treated in Corollary 2.3 , while our general estimate 2.2 applies for the derivatives of

z

. In view of these facts, our last assertion 2.18 is now easy to deduce using Leibniz' rule.

Our next step is to establish an even-dimensional analogue of the previous lemma. Here, our task is similar but harder, as the representation 1.9 for the Riemann operator is more subtle. Consider, for instance, the function

U_{0 m} (λ, r, t) \equiv \int_{z (λ, r, t)}^{1} [σ - z (λ, r, t)]^{- 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ, | t - r | \leq λ \leq t + r

that appears in 1.9 as an inner integral. At the endpoints

λ = r \pm t

, this function behaves rather nicely and actually attains the value

π / \sqrt{2}

. At the endpoint

λ = t - r

, however, the integral above happens to diverge. Namely, our definition 1.7 gives

z (t - r, r, t) = - 1

and this creates a singularity that is not integrable near

σ = - 1

. It is worth noting that we only have to deal with the troublesome endpoint when

t \geq r

. To handle this case, we shall need to introduce and study the functions

U_{j m} (λ, r, t) \equiv \int_{z (λ, r, t)}^{1} [σ - z (λ, r, t)]^{j - 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ, | t - r | \leq λ \leq t + r

for all integers

0 \leq j \leq m

. These will play a role similar to that of the polynomials

P_{j m}

we introduced in odd dimensions 2.11 .

Lemma 2.5. Fix an integer

m \geq 1

and let

T_{m}

be the polynomial 1.6 . With

z = z (λ, r, t)

given by 1.7 whenever

(r, t) \in R_{+}^{2}

, set

\begin{matrix} U_{j m} (λ, r, t) \equiv \int_{z}^{1} (σ - z)^{j - 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ \end{matrix}

(2.19)

for each

| t - r | \leq λ \leq t + r

and each integer

0 \leq j \leq m

(a) At the endpoint $λ = t + r$ , the function $U_{j m}$ has a zero of order $j$ .

(b) Letting

D_{*} = (\partial_{λ}, \partial_{r}, \partial_{t})

and

α

be any multi-index, one has

\begin{matrix} | D_{*}^{α} U_{0 m} (λ, r, t) | \leq C (α, m) \cdot {(\frac{λ}{r - t + λ})}^{1 / 2 + | α |} \cdot λ^{- | α |} \end{matrix}

(2.20)

whenever

0 < t \leq 2 r

and

| t - r | \leq λ \leq t + r

t \geq r

. Letting

D = (\partial_{r}, \partial_{t})

, one then has

\begin{matrix} | D^{β} U_{j m} (λ, r, t) | \leq C (m) \cdot {(\frac{λ r}{t + r})}^{1 / 2 - | β |} \cdot \frac{1}{\sqrt{r - t + λ}} \end{matrix}

(2.21)

for each

| t - r | \leq λ \leq t + r

and each multi-index

β

with

| β | \leq j \leq m

Proof. For the first two parts, we use the substitution

ν = (σ - z) / (1 - z)

to write

U_{j m} (λ, r, t) = (1 - z)^{j} \int_{0}^{1} \frac{ν^{j - 1 / 2}}{\sqrt{1 - ν}} \cdot \frac{T_{m} (σ)}{\sqrt{1 + σ}} d ν, σ = ν + (1 - ν) z .

In view of 1.7 , we have

z = 1

when

λ = t + r

, so part (a) is clear. To prove the estimate of part (b), we take

j = 0

in the last equation and differentiate to find that

\begin{matrix} D_{*}^{α} U_{0 m} (λ, r, t) = \int_{0}^{1} \frac{1}{\sqrt{ν (1 - ν)}} \cdot D_{*}^{α} (\frac{T_{m} (σ)}{\sqrt{1 + σ}}) d ν, σ = ν + (1 - ν) z . \end{matrix}

(2.22)

When it comes to the rightmost factor in the integrand, we have

D_{*}^{α} (\frac{T_{m} (σ)}{\sqrt{1 + σ}}) = \sum_{s = 1}^{| α |} \sum_{α_{1} + \dots + α_{s} = α} C (α_{1}, \dots, α_{s}, s) \cdot \frac{d^{s}}{d σ^{s}} (\frac{T_{m} (σ)}{\sqrt{1 + σ}}) \cdot^{s} \prod_{q = 1} (D_{*}^{α_{q}} σ) .

Since

λ

r

and

t

are as in Corollary 2.2 , an estimate of the form

| D_{*}^{α_{q}} z | \leq C (α_{q}) \cdot λ^{- | α_{q} |}

holds for

z

, hence also for

σ = ν + (1 - ν) z

because

0 \leq ν \leq 1

. In particular,

σ

is bounded and

| D_{*}^{α} (\frac{T_{m} (σ)}{\sqrt{1 + σ}}) | \leq C (α, m) \cdot (1 + σ)^{- 1 / 2 - | α |} \cdot λ^{- | α |}

by above. Inserting this fact in 2.22 , we then arrive at

| D_{*}^{α} U_{0 m} (λ, r, t) | \leq C λ^{- | α |} \cdot \int_{0}^{1} \frac{1}{\sqrt{ν (1 - ν)}} \cdot (1 + σ)^{- 1 / 2 - | α |} d ν, σ = ν + (1 - ν) z .

It is easy to check that

| z | \leq 1

for each

| t - r | \leq λ \leq t + r

, whence

1 + σ = 1 + ν + (1 - ν) z = 1 + z + (1 - z) ν \geq 1 + z

within the region of integration. Meanwhile, the definition 1.7 of

z

ensures that

1 + z = \frac{(r + t + λ) \cdot (r - t + λ)}{2 r λ} \geq \frac{r - t + λ}{λ}

for each

λ \geq | r - t |

. Once we now combine the last three equations, we obtain 2.20 .

Next, we turn to part (c). Here, we wish to differentiate

U_{j m}

a total of

| β | \leq j

times. As the integrand in 2.19 vanishes with order

j - 1 / 2

at the lower limit of integration, we get

\begin{matrix} D^{β} U_{j m} (λ, r, t) = \int_{z}^{1} D^{β} (σ - z)^{j - 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ . \end{matrix}

(2.23)

For the leftmost factor in the integrand, repeated applications of the chain rule give

| D^{β} (σ - z)^{j - 1 / 2} | \leq C (β) \cdot \sum_{s = 1}^{| β |} \sum_{β_{1} + \dots + β_{s} = β} (σ - z)^{j - s - 1 / 2} \cdot^{s} \prod_{q = 1} | D^{β_{q}} z | .

Moreover, one has

| z | \leq 1

whenever

| t - r | \leq λ \leq t + r

, so one easily finds

(σ - z)^{j - s - 1 / 2} \leq (1 - z)^{j - s} \cdot (σ - z)^{- 1 / 2} \leq 2^{j - s} \cdot (σ - z)^{- 1 / 2}

whenever

z \leq σ \leq 1

and

s \leq | β | \leq j

. Combining the last two equations, we then obtain

\begin{matrix} | D^{β} (σ - z)^{j - 1 / 2} | \leq C \sum_{s = 1}^{| β |} \sum_{β_{1} + \dots + β_{s} = β} (σ - z)^{- 1 / 2} \cdot^{s} \prod_{q = 1} | D^{β_{q}} z | . \end{matrix}

(2.24)

Since

t \geq r

by assumption, our auxiliary estimate 2.3 is applicable here. In fact, we have

r^{2} + | t^{2} - λ^{2} | \leq r (r + t + λ) \leq 2 r (r + t), | t - r | \leq λ \leq t + r

because

| t - λ | \leq r

for such

λ

, hence 2.3 ensures that

\begin{matrix} | D^{β_{q}} z | \leq C (β_{q}) \cdot \frac{r (t + r)}{λ r^{| β_{q} | + 1}} = C (β_{q}) \cdot \frac{t + r}{λ r^{| β_{q} |}}, | t - r | \leq λ \leq t + r . \end{matrix}

(2.25)

Once we now combine the last equation with 2.24 , we get

| D^{β} (σ - z)^{j - 1 / 2} | \leq C r^{- | β |} \sum_{s = 1}^{| β |} (σ - z)^{- 1 / 2} {(\frac{t + r}{λ})}^{s} \leq C (σ - z)^{- 1 / 2} {(\frac{t + r}{λ r})}^{| β |}

because

λ \leq t + r

. Inserting this fact in 2.23 , we thus arrive at

\begin{matrix} | D^{β} U_{j m} (λ, r, t) | \leq C {(\frac{t + r}{λ r})}^{| β |} \cdot \int_{z}^{1} (σ - z)^{- 1 / 2} \cdot (1 - σ^{2})^{- 1 / 2} \cdot | T_{m} (σ) | d σ . \end{matrix}

(2.26)

Note that the polynomial factor is bounded, as

| z | \leq 1

by above. If we also neglect a factor of

(1 + σ)^{- 1 / 2}

from the integrand, the resulting integral becomes

\int_{z}^{1} (σ - z)^{- 1 / 2} \cdot (1 - σ)^{- 1 / 2} d σ = \int_{0}^{1} ν^{- 1 / 2} (1 - ν)^{- 1 / 2} d ν = π

by means of the substitution

ν = (σ - z) / (1 - z)

. In particular, 2.26 leads to

| D^{β} U_{j m} (λ, r, t) | \leq C {(\frac{t + r}{λ r})}^{| β |} \cdot (1 + z)^{- 1 / 2} .

Besides, the definition 1.7 of

z

is such that

1 + z = \frac{(r + t + λ) \cdot (r - t + λ)}{2 λ r} \geq (\frac{t + r}{2 λ r}) \cdot (r - t + λ)

for each

λ \geq | t - r |

, so we may combine the last two equations to finally deduce 2.21 .

The functions of the previous lemma are related to the first component of the Riemann operator in even dimensions 1.9 . To treat the second component, we shall need to study some similar functions that we introduce below 2.27 . Although their definition resembles the one we had before 2.19 , we are now interested in the values

0 \leq λ \leq t - r

. For such values, the rational function 1.7 is no longer bounded and a much more delicate approach is needed. In fact, it is only here that the Tchebyshev polynomial 1.6 becomes important.

Lemma 2.6. Fix an integer

m \geq 1

and let

T_{m}

be the polynomial 1.6 . With

z = z (λ, r, t)

given by 1.7 whenever

(r, t) \in R_{+}^{2}

, set

\begin{matrix} W_{i m} (λ, r, t) \equiv \int_{- 1}^{1} (σ - z)^{i - 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ \end{matrix}

(2.27)

for each

0 \leq λ \leq t - r

and each integer

0 \leq i \leq m

(a) At the endpoint $λ = t - r$ , the function $W_{i m}$ agrees with the function $U_{i m}$ of 2.19 up to any derivative of order $i$ .

(b) Assuming that

0 \leq λ \leq t - r

, the estimate

\begin{matrix} | W_{i m} (λ, r, t) | \leq C (m) \cdot {(\frac{r}{t + r} \cdot \frac{λ}{t - r - λ})}^{i_{0} - i + 1 / 2} \end{matrix}

(2.28)

holds for all integers

0 \leq i \leq i_{0} \leq m

Proof. The only difference between the definitions of $U_{i m}$ and $W_{i m}$ is that the lower limit of integration is $z (λ, r, t)$ for the former and $- 1$ for the latter. By 1.7 , these two quantities agree when $λ = t - r$ , so part (a) follows easily. To prove the estimate of part (b), we use the definition 1.6 of the Tchebyshev polynomial $T_{m}$ to first write $W_{i m} (λ, r, t) = C (m) \int_{- 1}^{1} (σ - z)^{i - 1 / 2} \cdot \frac{d^{i_{0}}}{d σ^{i_{0}}} [\frac{d^{m - i_{0}}}{d σ^{m - i_{0}}} (1 - σ^{2})^{m - 1 / 2}] d σ .$ In the case that $i_{0} \geq 1$ , the expression in square brackets vanishes at $\pm 1$ . This means that we may integrate by parts $i_{0}$ times to obtain $W_{i m} (λ, r, t) = C (i_{0}, i, m) \int_{- 1}^{1} (σ - z)^{i - i_{0} - 1 / 2} \cdot [\frac{d^{m - i_{0}}}{d σ^{m - i_{0}}} (1 - σ^{2})^{m - 1 / 2}] d σ .$ Here, the expression in square brackets is integrable on $[- 1, 1]$ whenever $i_{0} \geq 0$ , so we get $| W_{i m} (λ, r, t) | \leq C (- 1 - z)^{i - i_{0} - 1 / 2}$ for all integers $0 \leq i \leq i_{0}$ . Since the definition 1.7 of $z$ provides the inequality $- 1 - z = \frac{r + t + λ}{2 r} \cdot \frac{t - r - λ}{λ} \geq \frac{t + r}{2 r} \cdot \frac{t - r - λ}{λ}$ for each $0 \leq λ \leq t - r$ , we may then combine the last two equations to deduce 2.28 .

Corollary 2.7. Assume that

m \geq 1

and

t \geq r

. Letting

D = (\partial_{r}, \partial_{t})

, one then has

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | \leq C (m) \cdot {(\frac{r}{t + r})}^{1 / 2 - | β |} \cdot \frac{λ^{m - j + 1 / 2}}{(t - r)^{m - j + | β |}} \cdot \frac{1}{\sqrt{t - r - λ}} \end{matrix}

(2.29)

for each

0 \leq λ \leq t - r

and each multi-index

β

with

| β | \leq j \leq m

Proof. Let us first differentiate 2.27 to get

D^{β} W_{j m} (λ, r, t) = \int_{- 1}^{1} D^{β} (σ - z)^{j - 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ .

Using the chain rule repeatedly, one may write the leftmost factor in the integrand as

D^{β} (σ - z)^{j - 1 / 2} = \sum_{s = 1}^{| β |} \sum_{β_{1} + \dots + β_{s} = β} C (β_{1}, \dots, β_{s}, s, j) \cdot (σ - z)^{j - s - 1 / 2} \cdot^{s} \prod_{q = 1} (D^{β_{q}} z) .

Recalling our definition 2.27 , one may then combine the last two equations to write

D^{β} W_{j m} (λ, r, t) = \sum_{s = 1}^{| β |} \sum_{β_{1} + \dots + β_{s} = β} C (β_{1}, \dots, β_{s}, s, j) \cdot W_{j - s, m} (λ, r, t) \cdot^{s} \prod_{q = 1} (D^{β_{q}} z) .

This trivially gives

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | \leq C \sum_{s = 1}^{| β |} \sum_{β_{1} + \dots + β_{s} = β} | W_{j - s, m} (λ, r, t) | \cdot^{s} \prod_{q = 1} | D^{β_{q}} z | \end{matrix}

(2.30)

and we shall now resort to our auxiliary estimate 2.3 . It is easy to check that

r^{2} + | t^{2} - λ^{2} | \geq r^{2} + r (t + λ) \geq r t, 0 \leq λ \leq t - r

because

t - λ \geq r

for such

λ

. In particular, our auxiliary estimate 2.3 gives

| D^{β_{q}} z | \leq C (β_{q}) \cdot \frac{r^{2} + t^{2} - λ^{2}}{λ r^{| β_{q} | + 1}} \leq \frac{C (β_{q})}{r^{| β_{q} |}} \cdot (\frac{t + r + λ}{λ} \cdot \frac{t + r - λ}{r}) .

Employing this fact in 2.30 , we thus arrive at

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} | W_{j - s, m} (λ, r, t) | \cdot {(\frac{t + r + λ}{λ} \cdot \frac{t + r - λ}{r})}^{s} . \end{matrix}

(2.31)

Here, let us note that

0 \leq | β | - s \leq j - s \leq j \leq m

since

| β | \leq j \leq m

by assumption.

Case 1:When

0 \leq λ \leq (t - r) / 2

, it suffices to show that

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | \leq C (m) \cdot {(\frac{r}{t + r})}^{1 / 2 - | β |} \cdot \frac{λ^{m - j + 1 / 2}}{(t - r)^{m - j + | β | + 1 / 2}} \end{matrix}

(2.32)

because

t - r - λ

is equivalent to

t - r

. For the exact same reason, we have

| W_{j - s, m} (λ, r, t) | \leq C (m) \cdot {(\frac{r}{t + r} \cdot \frac{λ}{t - r})}^{m - j + s + 1 / 2}

by Lemma 2.6 with

i = j - s

and

i_{0} = m

. Inserting this fact in 2.31 , we then get

| D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} {(\frac{r}{t + r} \cdot \frac{λ}{t - r})}^{m - j + 1 / 2} {(\frac{t + r + λ}{t + r} \cdot \frac{t + r - λ}{t - r})}^{s} .

Since

t + r \pm λ \leq 2 (t + r)

whenever

λ \leq t - r

, this also implies

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | & \leq C r^{- | β |} \cdot {(\frac{r}{t + r} \cdot \frac{λ}{t - r})}^{m - j + 1 / 2} {(\frac{t + r}{t - r})}^{| β |} \end{matrix}

\begin{matrix} = C {(\frac{r}{t + r})}^{m - j + 1 / 2 - | β |} \cdot \frac{λ^{m - j + 1 / 2}}{(t - r)^{m - j + | β | + 1 / 2}} . \end{matrix}

\begin{matrix}  \end{matrix}

In view of our assumption that

j \leq m

, we may now deduce the desired estimate 2.32 .

Case 2:When

(t - r) / 2 \leq λ \leq t - r

and

r \leq t \leq 2 r

, we need only show that

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | \leq \frac{C λ^{1 / 2 - | β |}}{\sqrt{t - r - λ}} . \end{matrix}

(2.33)

Since

t + r \pm λ \leq 2 (t + r) \leq 6 r

under the present assumptions, we easily get

| D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} | W_{j - s, m} (λ, r, t) | \cdot {(\frac{r}{λ})}^{s}

by means of 2.31 . Moreover,

λ \leq t - r \leq r

for this case, so we find that

| D^{β} W_{j m} (λ, r, t) | \leq C λ^{- | β |} \cdot \sum_{s = 1}^{| β |} | W_{j - s, m} (λ, r, t) | .

Applying Lemma 2.6 with

i = j - s = i_{0}

, we also have

| W_{j - s, m} (λ, r, t) | \leq C (m) \cdot {(\frac{r}{t + r} \cdot \frac{λ}{t - r - λ})}^{1 / 2} \leq \frac{C λ^{1 / 2}}{\sqrt{t - r - λ}}

so we may simply combine the last two equations to deduce the desired estimate 2.33 .

Case 3:When

(t - r) / 2 \leq λ \leq t - r

and

t \geq 2 r

, it suffices to show that

\begin{matrix} | D^{β} W_{j m} (λ, r, t) | \leq \frac{C r^{1 / 2 - | β |}}{\sqrt{t - r - λ}} \end{matrix}

(2.34)

because

λ

is equivalent to

t \pm r

. For the exact same reason, 2.31 reduces to

| D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} | W_{j - s, m} (λ, r, t) | \cdot {(\frac{t + r - λ}{r})}^{s},

and then Lemma 2.6 applies to give

| D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} {(\frac{r}{t - r - λ})}^{i_{0} - j + s + 1 / 2} {(\frac{t + r - λ}{r})}^{s}

for any integer

i_{0}

with

j - s \leq i_{0} \leq m

Subcase 3a:If it happens that

t - 3 r / 2 \leq λ

, the admissible choice

i_{0} = j - s

yields

| D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} {(\frac{r}{t - r - λ})}^{1 / 2} {(\frac{t + r - λ}{r})}^{s} \leq \frac{C r^{1 / 2 - | β |}}{\sqrt{t - r - λ}}

because

t + r - λ \leq 5 r / 2

whenever

t - 3 r / 2 \leq λ

Subcase 3b:If it happens that

λ \leq t - 3 r / 2

, the admissible choice

i_{0} = j

yields

| D^{β} W_{j m} (λ, r, t) | \leq C r^{- | β |} \cdot \sum_{s = 1}^{| β |} {(\frac{r}{t - r - λ})}^{s + 1 / 2} {(\frac{t + r - λ}{r})}^{s} \leq \frac{C r^{1 / 2 - | β |}}{\sqrt{t - r - λ}}

because

t + r - λ \leq 5 (t - r - λ)

whenever

λ \leq t - 3 r / 2

This establishes the desired 2.34 and also concludes the proof of the corollary.

We are finally in a position to prove an even-dimensional analogue of Lemma 2.4 .

Lemma 2.8. Let

n \geq 4

be an even integer and

a, m

be as in 1.3 . Suppose

f_{0} \in C^{l} (R_{+})

and

f_{1} \in C^{l + 1} (R_{+})

for some integer

1 \leq l \leq m

. Also, assume the singularity condition

\begin{matrix} \sum_{s = 0}^{i + l - 1} λ^{s} | f_{i}^{(s)} (λ) | = O (λ^{- 2 m - 2 + δ}) as λ \to 0 \end{matrix}

(2.35)

for

i = 0, 1

and some fixed

δ > 0

. When

t \geq r

, one then has the analogue

\begin{matrix} \partial_{t}^{i} [L f_{i}] (r, t) & = C_{j} \int_{t - r}^{t + r} [H_{i j} f_{i}] (λ, r, t) \cdot r^{- m - a} U_{j m} (λ, r, t) d λ \end{matrix}

\begin{matrix} + C_{j} \int_{0}^{t - r} [H_{i j} f_{i}] (λ, r, t) \cdot r^{- m - a} W_{j m} (λ, r, t) d λ \end{matrix}

(2.36)

\begin{matrix}  \end{matrix}

of 2.9 for

i = 0, 1

and all integers

i \leq j \leq l

. Here, the functions

U_{j m}

and

W_{j m}

are given by 2.19 and 2.27 , respectively, while the linear operators

H_{i j}

are the linear operators that we introduced in Lemma 2.4 .

Proof. For the case

i = 0

and

0 \leq j \leq l

, we shall use induction on

j

to verify the desired

\begin{matrix} [L f_{0}] (r, t) & = C_{j} \int_{t - r}^{t + r} [H_{0 j} f_{0}] (λ, r, t) \cdot r^{- m - a} U_{j m} (λ, r, t) d λ \end{matrix}

\begin{matrix} + C_{j} \int_{0}^{t - r} [H_{0 j} f_{0}] (λ, r, t) \cdot r^{- m - a} W_{j m} (λ, r, t) d λ . \end{matrix}

(2.37)

\begin{matrix}  \end{matrix}

When

j = 0

, this follows by our various definitions 1.3 , 1.9 , 2.12 , 2.19 and 2.27 . Let us now assume the last equation holds for some

j \geq 0

. In view of the definition

U_{j m} (λ, r, t) = \int_{z (λ, r, t)}^{1} [σ - z (λ, r, t)]^{j - 1 / 2} \cdot \frac{T_{m} (σ)}{\sqrt{1 - σ^{2}}} d σ

we introduced in 2.19 , it is easy to see that

U_{j m}

and

U_{j + 1, m}

are related by the formula

U_{j m} (λ, r, t) = - \frac{1}{j + 1 / 2} \cdot \frac{1}{\partial_{λ} z (λ, r, t)} \cdot \partial_{λ} U_{j + 1, m} (λ, r, t) .

Moreover, the same formula relates

W_{j m}

and

W_{j + 1, m}

because of our definition 2.27 . If we use this fact to integrate by parts the integrals of 2.37 , we then get to replace the leftmost factor in each integrand by

\frac{1}{j + 1 / 2} \cdot \frac{\partial}{\partial λ} (\frac{1}{\partial_{λ} z (λ, r, t)} \cdot [H_{0 j} f_{0}] (λ, r, t)) = \frac{1}{j + 1 / 2} \cdot [H_{0, j + 1} f_{0}] (λ, r, t)

in view of our definition 2.12 . Since the integration by parts also replaces

U_{j m}

U_{j + 1, m}

and

W_{j m}

W_{j + 1, m}

, we may thus complete the inductive proof of 2.37 by merely showing that no boundary terms arise in the process. Now, the boundary term at

λ = t + r

is zero because

U_{j + 1, m}

vanishes there by part (a) of Lemma 2.5 . At

λ = t - r

, we have two boundary terms with opposite signs. In fact, the functions

U_{j + 1, m}

and

W_{j + 1, m}

agree at that point by part (a) of Lemma 2.6 , so these two boundary terms plainly cancel one another. As for the boundary term at

λ = 0

, that one is given by

{[- \frac{C_{j}}{j + 1 / 2} \cdot \frac{[H_{0 j} f_{0}] (λ, r, t)}{\partial_{λ} z (λ, r, t)} \cdot r^{- m - a} W_{j + 1, m} (λ, r, t)]}_{λ = 0} .

Using 2.5 , 2.10 and 2.29 together with our singularity condition 2.35 on

f_{0}

, one can easily check that this expression vanishes as well. Thus, no boundary terms arise, indeed. Finally, to establish 2.36 for the case

i = 1

and

1 \leq j \leq l

, one proceeds as in the proof of Lemma 2.4 . Since only minor modifications are needed here, we shall omit the details.

3 The Riemann Operator in High Dimensions

Using the results of the previous section, we shall now study the Riemann operator, which is defined by either 1.8 or 1.9 , according to the parity of

n

. To estimate its derivatives, we divide our analysis into two cases that we treat separately in our next two propositions.

First, we focus on the interior region

t \geq 2 r

and provide an estimate that does not involve any boundary terms. For odd values of

n

, a similar estimate appears in [2] . The main idea we use to handle this region coincides with the one in [8] , where the wave equation with the inverse-square potential is studied. However, the approach in [8] is quite different from ours.

Proposition 3.1. Let

n \geq 4

be an integer and define

a, m

by 1.3 . Suppose

f_{0} \in C^{l} (R_{+})

and

f_{1} \in C^{l + 1} (R_{+})

for some integer

1 \leq l \leq m

. Also, assume the singularity condition 2.35 for

i = 0, 1

and some fixed

δ > 0

. When

D = (\partial_{r}, \partial_{t})

and

t \geq 2 r

, the Riemann operator is then such that

\begin{matrix} | D^{β} \partial_{t}^{i} [L f_{i}] (r, t) | & \leq C_{1} (n) \cdot r^{j - | β | - m - a} \int_{t - r}^{t + r} \frac{λ^{m - j}}{(r - t + λ)^{1 - a}} \cdot \sum_{s = 0}^{i + j} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | d λ \end{matrix}

\begin{matrix} + \frac{C_{2} (n) \cdot r^{j - | β | - m}}{(t - r)^{j + m + a}} \int_{0}^{t - r} \frac{λ^{2 m}}{(t - r - λ)^{1 - a}} \cdot \sum_{s = 0}^{i + j} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | d λ \end{matrix}

(3.1)

\begin{matrix}  \end{matrix}

for

i = 0, 1

and each

max (i, | β |) \leq j \leq l

. Moreover, one has

C_{2} (n) = 0

when

n

is odd.

Proof. We divide our analysis into two cases, according to the parity of

n

Case 1:For odd values of

n

, an application of Lemma 2.4 allows us to write

\partial_{t}^{i} [L f_{i}] (r, t) = \frac{(- 1)^{j}}{2} \int_{t - r}^{t + r} [H_{i j} f_{i}] (λ, r, t) \cdot r^{- m - a} P_{j m} (z (λ, r, t)) d λ, i \leq j \leq l .

Recall that

P_{j m}

is a polynomial with a zero of order

j

at each of

\pm 1

, while

z (t \pm r, r, t) = \pm 1

by 1.7 . In particular, we may differentiate the last equation

| β | \leq j

times without creating any boundary terms to obtain

\begin{matrix} | D^{β} \partial_{t}^{i} [L f_{i}] | \leq C \int_{t - r}^{t + r} \sum_{β_{1} + β_{2} + β_{3} = β} | D^{β_{1}} [H_{i j} f_{i}] | \cdot r^{- m - a - | β_{2} |} \cdot | D^{β_{3}} P_{j m} (z (λ, r, t)) | d λ . \end{matrix}

(3.2)

When it comes to the operators

H_{i j}

, the general estimate

\begin{matrix} | D^{β_{1}} [H_{i j} f_{i}] | \leq \frac{C λ^{j} r^{j - | β_{1} |} t^{| β_{1} | - j}}{(t - r)^{i + j + | β_{1} |}} \cdot \sum_{s = 0}^{i + j} λ^{m + a + s} | f_{i}^{(s)} (λ) |, 0 \leq λ \leq t + r \end{matrix}

(3.3)

is provided by Lemma 2.4 . Under our assumption that

t \geq 2 r

, it actually implies

\begin{matrix} | D^{β_{1}} [H_{i j} f_{i}] | \leq C λ^{- i - j} r^{j - | β_{1} |} \cdot \sum_{s = 0}^{i + j} λ^{m + a + s} | f_{i}^{(s)} (λ) |, t - r \leq λ \leq t + r \end{matrix}

(3.4)

because

λ

t \pm r

and

t

are all equivalent here. For the exact same reason, 2.25 gives

| D^{γ} z (λ, r, t) | \leq C (γ) \cdot \frac{t + r}{λ r^{| γ |}} \leq \frac{C}{r^{| γ |}}, t - r \leq λ \leq t + r

for any multi-index

γ

. Moreover, one has

| z (λ, r, t) | \leq 1

for such

λ

, so this easily leads to

| D^{β_{3}} P_{j m} (z (λ, r, t)) | \leq C r^{- | β_{3} |} .

Using the last equation and 3.4 , we may then estimate 3.2 as

| D^{β} \partial_{t}^{i} [L f_{i}] (r, t) | \leq C r^{j - | β | - m - a} \int_{t - r}^{t + r} λ^{m - j} \cdot \sum_{s = 0}^{i + j} λ^{s + a - i} | f_{i}^{(s)} (λ) | d λ .

Since

a = 1

when

n

is odd, this is precisely the desired estimate 3.1 with

C_{2} (n) = 0

Case 2:For even values of

n

, an application of Lemma 2.8 allows us to write

\begin{matrix} \partial_{t}^{i} [L f_{i}] (r, t) & = C_{j} \int_{t - r}^{t + r} [H_{i j} f_{i}] (λ, r, t) \cdot r^{- m - a} U_{j m} (λ, r, t) d λ \end{matrix}

\begin{matrix} + C_{j} \int_{0}^{t - r} [H_{i j} f_{i}] (λ, r, t) \cdot r^{- m - a} W_{j m} (λ, r, t) d λ, i \leq j \leq l . \end{matrix}

\begin{matrix}  \end{matrix}

By Lemma 2.5 , the function

U_{j m}

vanishes with order

j

at the endpoint

λ = t + r

. Also, its derivatives up to order

j

agree with those of

W_{j m}

at the endpoint

λ = t - r

by Lemma 2.6 .

Thus, no boundary terms arise upon

| β | \leq j

differentiations of the last equation, so we get

\begin{matrix} | D^{β} \partial_{t}^{i} [L f_{i}] | & \leq C \int_{t - r}^{t + r} \sum_{β_{1} + β_{2} + β_{3} = β} | D^{β_{1}} [H_{i j} f_{i}] | \cdot r^{- m - a - | β_{2} |} \cdot | D^{β_{3}} U_{j m} (λ, r, t) | d λ \end{matrix}

\begin{matrix} + C \int_{0}^{t - r} \sum_{β_{1} + β_{2} + β_{3} = β} | D^{β_{1}} [H_{i j} f_{i}] | \cdot r^{- m - a - | β_{2} |} \cdot | D^{β_{3}} W_{j m} (λ, r, t) | d λ \end{matrix}

\begin{matrix} \equiv A_{1} + A_{2} . \end{matrix}

(3.5)

\begin{matrix}  \end{matrix}

Within the first integral,

λ

is equivalent to

t \pm r

as before, so an estimate of the form

| D^{β_{3}} U_{j m} (λ, r, t) | \leq \frac{C r^{1 / 2 - | β_{3} |}}{\sqrt{r - t + λ}}, t - r \leq λ \leq t + r

holds by part (c) of Lemma 2.5 . Using this fact along with 3.4 , we then find that

A_{1} \leq C r^{j - | β | - m - a + 1 / 2} \int_{t - r}^{t + r} \frac{λ^{m - j}}{\sqrt{r - t + λ}} \cdot \sum_{s = 0}^{i + j} λ^{s + a - i} | f_{i}^{(s)} (λ) | d λ .

Moreover, we have

r \leq t - r \leq λ

within the region of integration, hence also

A_{1} \leq C r^{j - | β | - m - a} \int_{t - r}^{t + r} \frac{λ^{m - j}}{\sqrt{r - t + λ}} \cdot \sum_{s = 0}^{i + j} λ^{s + a + 1 / 2 - i} | f_{i}^{(s)} (λ) | d λ .

Since

a = 1 / 2

when

n

is even, the desired estimate 3.1 is thus satisfied by

A_{1}

To treat the second integral

A_{2}

, we recall that

t

is equivalent to

t \pm r

because

t \geq 2 r

by assumption. An immediate consequence of 3.3 is then

| D^{β_{1}} [H_{i j} f_{i}] | \leq \frac{C λ^{j} r^{j - | β_{1} |}}{(t - r)^{i + 2 j}} \cdot \sum_{s = 0}^{i + j} λ^{m + a + s} | f_{i}^{(s)} (λ) |, 0 \leq λ \leq t - r

and an immediate consequence of Corollary 2.7 is

| D^{β_{3}} W_{j m} (λ, r, t) | \leq \frac{C r^{1 / 2 - | β_{3} |} \cdot λ^{m - j + 1 / 2}}{(t - r)^{m - j + 1 / 2}} \cdot \frac{1}{\sqrt{t - r - λ}}, 0 \leq λ \leq t - r .

Employing these two facts in 3.5 , we now get

A_{2} \leq \frac{C r^{j - | β | - m - a + 1 / 2}}{(t - r)^{i + j + m + 1 / 2}} \cdot \int_{0}^{t - r} \frac{λ^{2 m}}{\sqrt{t - r - λ}} \cdot \sum_{s = 0}^{i + j} λ^{s + a + 1 / 2} | f_{i}^{(s)} (λ) | d λ .

Moreover,

λ \leq t - r

within the region of integration, so this also implies

A_{2} \leq \frac{C r^{j - | β | - m - a + 1 / 2}}{(t - r)^{j + m + 1 / 2}} \cdot \int_{0}^{t - r} \frac{λ^{2 m}}{\sqrt{t - r - λ}} \cdot \sum_{s = 0}^{i + j} λ^{s + a + 1 / 2 - i} | f_{i}^{(s)} (λ) | d λ .

Since

a = 1 / 2

when

n

is even, the desired estimate 3.1 is thus satisfied by

A_{2}

as well.

In our next proposition, we provide a similar analysis for the exterior region

t \leq 2 r

. The most notable difference from our previous estimate is the presence of some boundary terms.

A more subtle difference is that the regularity of the initial data is not so crucial anymore.

Namely, our sharpest conclusions occur for the smallest possible value of

j

in the following

Proposition 3.2. When

t \leq 2 r

, the assumptions of the previous proposition imply

\begin{matrix} | D^{β} \partial_{t}^{i} [L f_{i}] (r, t) | & \leq C_{1}^{'} (n) \cdot r^{- m - a} \int_{| t - r |}^{t + r} \frac{λ^{m - | β |}}{(r - t + λ)^{1 - a}} \cdot \sum_{s = 0}^{i + j} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | d λ \end{matrix}

\begin{matrix} + \frac{C_{2}^{'} (n) \cdot r^{- m - a}}{(t - r)^{m + | β |}} \int_{0}^{max (t - r, 0)} \frac{λ^{2 m}}{(t - r - λ)^{1 - a}} \cdot \sum_{s = 0}^{i + j} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | d λ \end{matrix}

\begin{matrix} + C_{1}^{'} (n) \cdot r^{- m - a} \sum_{s = 0}^{i + j - 1} {[λ^{m + a - | β | + s + 1 - i} | f_{i}^{(s)} (λ) |]}_{λ = | t \pm r |} \end{matrix}

(3.6)

\begin{matrix}  \end{matrix}

for

i = 0, 1

and each

max (i, | β |) \leq j \leq l

. Moreover, one has

C_{2}^{'} (n) = 0

when

n

is odd.

Proof. We divide our analysis into three cases.

Case 1:For odd values of

n

, the Riemann operator 1.8 is of the form

\begin{matrix} [L f_{i}] (r, t) = \frac{1}{2} \int_{| t - r |}^{t + r} λ^{m + a} f_{i} (λ) \cdot r^{- m - a} P_{m} (z (λ, r, t)) d λ . \end{matrix}

(3.7)

To obtain the desired estimate in this case, we shall not have to distinguish between radial and time derivatives, so it is convenient to introduce a multi-index

γ

of order

| γ | = | β | + i

Once we differentiate the last equation, we then get an identity of the form

\begin{matrix} D^{γ} [L f_{i}] (r, t) & = \frac{1}{2} \int_{| t - r |}^{t + r} λ^{m + a} f_{i} (λ) \cdot D^{γ} (r^{- m - a} P_{m} (z (λ, r, t))) d λ \end{matrix}

\begin{matrix} + \sum_{| γ_{1} | = | γ | - 1} C_{\pm} (γ_{1}) \cdot D^{γ_{1}} {[λ^{m + a} f_{i} (λ) \cdot r^{- m - a} P_{m} (z (λ, r, t))]}_{λ = | t \pm r |} \end{matrix}

\begin{matrix} \equiv B_{1} + B_{2} . \end{matrix}

(3.8)

\begin{matrix}  \end{matrix}

Since

| t - r | \leq λ \leq t + r

within these terms, our assumption

t \leq 2 r

makes Corollary 2.2 applicable. Letting

D_{*} = (\partial_{λ}, \partial_{r}, \partial_{t})

, we thus have the estimate

| D_{*}^{α} z (λ, r, t) | \leq C (α) \cdot λ^{- | α |}, | t - r | \leq λ \leq t + r

for any multi-index

α

. Besides,

| z (λ, r, t) | \leq 1

for such

λ

, so we easily get

\begin{matrix} | D_{*}^{α} P_{m} (z (λ, r, t)) | \leq C (α, m) \cdot λ^{- | α |}, | t - r | \leq λ \leq t + r \end{matrix}

(3.9)

by repeated applications of the chain rule. This is actually the only fact we need to control the right hand side of 3.8 . Due to our assumption that

t \leq 2 r

, it also implies

\begin{matrix} | D_{*}^{α} (r^{- m - a} P_{m} (z (λ, r, t))) | \leq C r^{- m - a} \cdot λ^{- | α |}, | t - r | \leq λ \leq t + r \end{matrix}

(3.10)

for any multi-index

α

, as

λ \leq 3 r

here. Applying the last inequality to the integral term

B_{1}

, let us now recall that

| γ | = | β | + i

to find

\begin{matrix} | B_{1} | \leq C r^{- m - a} \int_{| t - r |}^{t + r} λ^{m + a - | γ |} \cdot | f_{i} (λ) | d λ = C r^{- m - a} \int_{| t - r |}^{t + r} λ^{m - | β |} \cdot λ^{a - i} | f_{i} (λ) | d λ . \end{matrix}

\begin{matrix}  \end{matrix}

Since

a = 1

when

n

is odd, the desired estimate 3.6 thus follows for

B_{1}

. Applying 3.10 to the boundary terms

B_{2}

, we similarly get

\begin{matrix} | B_{2} | & \leq \sum_{| γ_{1} | = | β | + i - 1} \sum_{s_{1} + s_{2} = | γ_{1} |} C r^{- m - a} \cdot {[λ^{m + a - s_{1}} | f_{i}^{(s_{2})} (λ) |]}_{λ = | t \pm r |} \end{matrix}

\begin{matrix} = C r^{- m - a} \sum_{s_{2} = 0}^{i + | β | - 1} {[λ^{m + a - | β | + s_{2} + 1 - i} | f_{i}^{(s_{2})} (λ) |]}_{λ = | t \pm r |} . \end{matrix}

\begin{matrix}  \end{matrix}

Moreover,

| β | \leq j

by assumption, so the desired estimate 3.6 follows for

B_{2}

as well.

Case 2:When

n

is even and

t \leq r

, the Riemann operator 1.9 takes the form

\begin{matrix} [L f_{i}] (r, t) = C_{0}^{'} \int_{| t - r |}^{t + r} λ^{m + a} f_{i} (λ) \cdot r^{- m - a} U_{0 m} (λ, r, t) d λ \end{matrix}

(3.11)

with

U_{0 m}

given by 2.19 . This closely resembles the equation 3.7 we used in the previous case, although the factor

P_{m} (z)

has been replaced by

U_{0 m}

and the value of

a

has changed due to the change in the parity of

n

. Differentiating directly as before, we now get

D^{γ} [L f_{i}] (r, t) = B_{1}^{'} + B_{2}^{'},

where each

B_{k}^{'}

denotes the corresponding

B_{k}

of 3.8 with

U_{0 m}

instead of

P_{m} (z)

. Since

t \leq 2 r

by assumption, part (b) of Lemma 2.5 applies to give

\begin{matrix} | D_{*}^{α} U_{0 m} (λ, r, t) | \leq C λ^{- | α |} \cdot {(\frac{λ}{r - t + λ})}^{1 / 2 + | α |}, | t - r | \leq λ \leq t + r \end{matrix}

(3.12)

for any multi-index

α

. This provides the first fact that we are going to need; the second one is the inequality

λ \leq r - t + λ

, which trivially holds for the case

t \leq r

under consideration.

Combining these two facts, we see that 3.9 remains valid when

P_{m} (z)

is replaced by

U_{0 m}

In particular, the argument of Case 1 applies verbatim, except for the part where the exact value of

a

was invoked. Said differently, we need only worry about the integral term

\begin{matrix} B_{1}^{'} = C_{0}^{'} \int_{| t - r |}^{t + r} λ^{m + a} f_{i} (λ) \cdot D^{γ} (r^{- m - a} U_{0 m} (λ, r, t)) d λ . \end{matrix}

(3.13)

If we combine 3.12 with our inequality

λ \leq r - t + λ

once again, we get

| D_{*}^{α} U_{0 m} (λ, r, t) | \leq \frac{C λ^{1 / 2 - | α |}}{\sqrt{r - t + λ}}, | t - r | \leq λ \leq t + r

for any multi-index

α

. As long as

λ \leq t + r \leq 3 r

, we then get

| D_{*}^{α} (r^{- m - a} U_{0 m} (λ, r, t)) | \leq \frac{C r^{- m - a} \cdot λ^{1 / 2 - | α |}}{\sqrt{r - t + λ}}, | t - r | \leq λ \leq t + r

for any multi-index

α

. Applying this fact to the integral term 3.13 , we now find

\begin{matrix} | B_{1}^{'} | \leq C r^{- m - a} \int_{| t - r |}^{t + r} \frac{λ^{1 / 2 - | γ |}}{\sqrt{r - t + λ}} \cdot λ^{m + a} | f_{i} (λ) | d λ . \end{matrix}

\begin{matrix}  \end{matrix}

Since

a = 1 / 2

when

n

is even and

| γ | = | β | + i

by above, we thus obtain

\begin{matrix} | B_{1}^{'} | \leq C r^{- m - a} \int_{| t - r |}^{t + r} \frac{λ^{m - | β |}}{(r - t + λ)^{1 - a}} \cdot λ^{1 - i} | f_{i} (λ) | d λ . \end{matrix}

\begin{matrix}  \end{matrix}

In particular, we obtain the desired estimate 3.6 for this case as well.

Case 3:Suppose now that

n

is even and

r < t \leq 2 r

. Although our estimate 3.12 behaves nicely at the boundary terms

λ = r \pm t

we had before, this is not the case for the boundary term

λ = t - r

which may now emerge. One may overcome this difficulty using our approach in the previous proposition to avoid the boundary terms altogether. However, that approach does not lead to the desired estimate, so we shall need to resort to a combination of the two.

Proceeding to the details, let us fix some

ζ \in C^{\infty} (R)

such that

\begin{matrix} ζ (x) = {\begin{matrix} 1 & if x \leq 2 \\ 0 & if x \geq 5 / 2 \end{matrix}} \end{matrix}

(3.14)

and introduce the functions

\begin{matrix} ζ_{1} (λ) = 1 - ζ (\frac{λ}{t - r}), ζ_{2} (λ) = ζ (\frac{λ}{t - r}) . \end{matrix}

(3.15)

Then the linearity of the Riemann operator allows us to write

[L f_{i}] (r, t) = [L (ζ_{1} f_{i})] (r, t) + [L (ζ_{2} f_{i})] (r, t)

as a sum of two terms to be treated separately.

In order to treat

[L (ζ_{1} f_{i})]

, we proceed as in Case 2. Since

ζ_{1}

vanishes when

λ \leq 2 (t - r)

by definition, the Riemann operator 1.9 takes the form

[L (ζ_{1} f_{i})] (r, t) = C_{0}^{'} \int_{2 (t - r)}^{t + r} λ^{m + a} ζ_{1} (λ) f_{i} (λ) \cdot r^{- m - a} U_{0 m} (λ, r, t) d λ .

This closely resembles the equation 3.11 we had in Case 2, although the extra factor

ζ_{1} (λ)

is now present in the integrand. When it comes to the first fact 3.12 we used in Case 2, it is provided by Lemma 2.5 whenever

t \leq 2 r

, so it is still valid. As for the inequality we used, its analogue

λ \leq 2 (r - t + λ)

now holds since

λ \geq 2 (t - r)

within the region of integration. The main question is then whether the extra factor

ζ_{1} (λ)

will bring any changes to our previous approach. Although this function is bounded by our definition 3.15 , its derivatives involve powers of

t - r

we did not have before. Nevertheless, the only place where such derivatives may occur is at the endpoint

λ = t + r

, and one has

t + r \geq 3 (t - r)

whenever

t \leq 2 r

. Once we now recall our definition 3.14 - 3.15 , we see that all derivatives of

ζ_{1}

vanish at

t + r

. In particular, one does not have to deal with such derivatives at all.

Next, we focus on

[L (ζ_{2} f_{i})]

. For

λ

sufficiently close to zero,

ζ_{2} (λ) \equiv 1

by 3.14 - 3.15 , so the functions

ζ_{2} f_{i}

and

f_{i}

satisfy the same singularity condition as

λ \to 0

. This means that we may employ Lemma 2.8 and argue as in the previous proposition to find

\begin{matrix} | D^{β} \partial_{t}^{i} [L (ζ_{2} f_{i})] | & \leq C \int_{t - r}^{t + r} \sum_{β_{1} + β_{2} + β_{3} = β} | D^{β_{1}} [H_{i j} (ζ_{2} f_{i})] | \cdot r^{- m - a - | β_{2} |} \cdot | D^{β_{3}} U_{j m} | d λ \end{matrix}

\begin{matrix} + C \int_{0}^{t - r} \sum_{β_{1} + β_{2} + β_{3} = β} | D^{β_{1}} [H_{i j} (ζ_{2} f_{i})] | \cdot r^{- m - a - | β_{2} |} \cdot | D^{β_{3}} W_{j m} | d λ \end{matrix}

\begin{matrix} \equiv B_{1}^{''} + B_{2}^{''} \end{matrix}

(3.16)

\begin{matrix}  \end{matrix}

for

i = 0, 1

and each

max (i, | β |) \leq j \leq l

. Namely, the condition

i \leq j \leq l

allows us to invoke the identity of Lemma 2.8 and the condition

| β | \leq j

to differentiate that identity

| β |

times without introducing any boundary terms.

Let us note that

3 (t - r) \leq t + r

whenever

t \leq 2 r

and temporarily assume the estimate

\begin{matrix} | D^{β_{1}} [H_{i j} (ζ_{2} f_{i})] | \leq \frac{C λ^{j}}{(t - r)^{i + j + | β_{1} |}} \cdot \sum_{s_{2} = 0}^{i + j} λ^{m + a + s_{2}} | f_{i}^{(s_{2})} (λ) |, 0 \leq λ \leq 3 (t - r) \end{matrix}

(3.17)

whose proof is given below. For the remaining values

λ \geq 3 (t - r)

, the left hand side is zero because

ζ_{2}

vanishes for such

λ

by 3.14 - 3.15 . Employing 3.17 , one thus arrives at

B_{1}^{''} \leq C r^{- m - a} \int_{t - r}^{3 (t - r)} \sum_{β_{1} + β_{2} + β_{3} = β} r^{- | β_{2} |} \cdot | D^{β_{3}} U_{j m} | \cdot λ^{- i - | β_{1} |} \sum_{s_{2} = 0}^{i + j} λ^{m + a + s_{2}} | f_{i}^{(s_{2})} (λ) | d λ .

Since

r

is equivalent to

t + r

whenever

t \leq 2 r

, part (c) of Lemma 2.5 ensures that

| D^{β_{3}} U_{j m} (λ, r, t) | \leq \frac{C λ^{1 / 2 - | β_{3} |}}{\sqrt{r - t + λ}}, t - r \leq λ \leq 3 (t - r)

because

3 (t - r) \leq t + r

by above. Once we now combine the last two equations, we get

B_{1}^{''} \leq C r^{- m - a} \int_{t - r}^{3 (t - r)} \sum_{β_{1} + β_{2} + β_{3} = β} \frac{r^{- | β_{2} |} λ^{m - | β_{1} | - | β_{3} |}}{\sqrt{r - t + λ}} \cdot \sum_{s_{2} = 0}^{i + j} λ^{s_{2} + a + 1 / 2 - i} | f_{i}^{(s_{2})} (λ) | d λ .

Moreover,

r \geq t - r \geq λ / 3

within the region of integration, so this also gives

B_{1}^{''} \leq C r^{- m - a} \int_{t - r}^{t + r} \frac{λ^{m - | β |}}{\sqrt{r - t + λ}} \cdot \sum_{s_{2} = 0}^{i + j} λ^{s_{2} + a + 1 / 2 - i} | f_{i}^{(s_{2})} (λ) | d λ .

Since

a = 1 / 2

when

n

is even, the desired estimate 3.6 is thus satisfied by

B_{1}^{''}

Next, we focus on

B_{2}^{''}

. Noting that

r \leq t + r \leq 3 r

by assumption, we get

| D^{β_{3}} W_{j m} (λ, r, t) | \leq \frac{C λ^{m - j + 1 / 2}}{(t - r)^{m - j + | β_{3} |}} \cdot \frac{1}{\sqrt{t - r - λ}}, 0 \leq λ \leq t - r

by means of Corollary 2.7 . Together with 3.17 , this allows us to estimate

B_{2}^{''}

B_{2}^{''} \leq \frac{C r^{- m - a}}{(t - r)^{i + m}} \int_{0}^{t - r} \sum_{β_{1} + β_{2} + β_{3} = β} \frac{r^{- | β_{2} |}}{(t - r)^{| β_{1} | + | β_{3} |}} \cdot \frac{λ^{m + 1 / 2}}{\sqrt{t - r - λ}} \cdot \sum_{s_{2} = 0}^{i + j} λ^{m + a + s_{2}} | f_{i}^{(s_{2})} (λ) | d λ .

Moreover,

r^{- | β_{2} |} (t - r)^{- i} \leq (t - r)^{- | β_{2} |} λ^{- i}

whenever

λ \leq t - r \leq r

, so we find

B_{2}^{''} \leq \frac{C r^{- m - a}}{(t - r)^{m + | β |}} \int_{0}^{t - r} \frac{λ^{2 m}}{\sqrt{t - r - λ}} \cdot \sum_{s_{2} = 0}^{i + j} λ^{s_{2} + a + 1 / 2 - i} | f_{i}^{(s_{2})} (λ) | d λ .

Since

a = 1 / 2

when

n

is even, the desired estimate 3.6 is thus satisfied by

B_{2}^{''}

as well.

To finish the proof, it remains to establish 3.17 . According to Lemma 2.4 , we do have

| D^{β_{1}} [H_{i j} (ζ_{2} f_{i})] | \leq \frac{C λ^{j} r^{j - | β_{1} |} t^{| β_{1} | - j}}{(t - r)^{i + j + | β_{1} |}} \cdot \sum_{s = 0}^{i + j} λ^{m + a + s} \cdot | \frac{d^{s}}{d λ^{s}} (ζ_{2} (λ) f_{i} (λ)) |

for each

0 \leq λ \leq 3 (t - r)

because

3 (t - r) \leq t + r

by above. For the case

r < t \leq 2 r

we are presently considering, our remaining assertion 3.17 will then follow once we show that

\begin{matrix} \sum_{s = 0}^{i + j} λ^{s} \cdot | \frac{d^{s}}{d λ^{s}} (ζ_{2} (λ) f_{i} (λ)) | \leq C \sum_{s_{2} = 0}^{i + j} λ^{s_{2}} \cdot | f_{i}^{(s_{2})} (λ) |, 0 \leq λ \leq 3 (t - r) . \end{matrix}

(3.18)

Recalling our definition 3.14 - 3.15 , we easily get this when

2 (t - r) \leq λ \leq 3 (t - r)

since

\sum_{s = 0}^{i + j} λ^{s} \cdot | \frac{d^{s}}{d λ^{s}} (ζ_{2} (λ) f_{i} (λ)) | \leq C \sum_{s = 0}^{i + j} \sum_{s_{2} = 0}^{s} λ^{s} (t - r)^{s_{2} - s} \cdot | f_{i}^{(s_{2})} (λ) | \leq C \sum_{s_{2} = 0}^{i + j} λ^{s_{2}} \cdot | f_{i}^{(s_{2})} (λ) |

by Leibniz' rule. When

0 \leq λ < 2 (t - r)

, on the other hand, our definition 3.14 - 3.15 is such that all derivatives of

ζ_{2}

vanish at

λ

, whence

\sum_{s = 0}^{i + j} λ^{s} \cdot | \frac{d^{s}}{d λ^{s}} (ζ_{2} (λ) f_{i} (λ)) | = \sum_{s = 0}^{i + j} λ^{s} \cdot | ζ_{2} (λ) | \cdot | f_{i}^{(s)} (λ) | \leq C \sum_{s = 0}^{i + j} λ^{s} \cdot | f_{i}^{(s)} (λ) | .

In any case whatsoever, the desired estimate 3.18 follows and the proof is complete.

4 Estimates for the Free Solution

Using the results of the previous section, we now study the solution

u_{0}

of the homogeneous equation 1.2 . To prove the estimates that Theorem 1.1 asserts for its derivatives, we shall also need two elementary facts that we list separately in our next two lemmas. Their proofs follow those of similar results in [1, 11] , for instance, so we are going to omit them. As for the parameter

a > 0

we introduce below, we are merely interested in the special cases

a = 1

and

a = 1 / 2

that we had before 1.3 .

Lemma 4.1. Let

(r, t) \in R_{+}^{2}

be arbitrary. Assuming that

a > 0

, one has

\int_{| t - r |}^{t + r} \frac{{〈 λ 〉}^{b} d λ}{(r - t + λ)^{1 - a}} \leq {\begin{matrix} C r^{a} {〈 t + r 〉}^{b} & if & b > - a \\ C r^{a} {〈 t + r 〉}^{- a} (1 + ln \frac{〈 t + r 〉}{〈 t - r 〉}) & if & b = - a \\ C r^{a} {〈 t + r 〉}^{- a} {〈 t - r 〉}^{a + b} & if & b < - a \end{matrix}}

for some constant

C

depending solely on

a

and

b

Lemma 4.2. Let

a > 0

be arbitrary. Assuming that

b \geq 0

and

t \geq r > 0

, the integral

ℬ (b, c) \equiv \int_{0}^{t - r} \frac{λ^{b} {〈 λ 〉}^{c} d λ}{(t - r - λ)^{1 - a}}

satisfies an estimate of the form

ℬ (b, c) \leq {\begin{matrix} C (t - r)^{a + b} \cdot {〈 t - r 〉}^{c} & if & b + c > - 1 \\ C (t - r)^{a + b} \cdot {〈 t - r 〉}^{- b - 1} \cdot (1 + ln 〈 t - r 〉) & if & b + c = - 1 \\ C (t - r)^{a + b} \cdot {〈 t - r 〉}^{- b - 1} & if & b + c < - 1 \end{matrix}},

where the constant

C

depends solely on

a, b

and

c

Proof of Theorem 1.1 . Let

u_{0} (r, t)

be the solution of the homogeneous equation 1.2 provided by Lemma 1.3 . Setting

f_{0} = ψ

and

f_{1} = φ

, we may then write

D^{β} u_{0} (r, t) = D^{β} [L ψ] (r, t) + D^{β} \partial_{t} [L φ] (r, t) = \sum_{i = 0}^{1} D^{β} \partial_{t}^{i} [L f_{i}] (r, t) .

This is also the expression we wish to estimate when

| β | \leq l

. Using the notation above, let us now write our assumption 1.4 on the initial data simply as

\begin{matrix} \sum_{i = 0}^{1} \sum_{s = 0}^{i + l} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | \leq ɛ λ^{l - m} {〈 λ 〉}^{m - l - k} . \end{matrix}

(4.1)

Here, one may readily check that the singularity condition 2.35 holds for

i = 0, 1

Step 1:In the interior region

t \geq 2 r

, we use Proposition 3.1 with

j = l

to obtain

\begin{matrix} | D^{β} u_{0} (r, t) | & \leq C_{1} r^{l - | β | - m - a} \int_{t - r}^{t + r} \frac{λ^{m - l}}{(r - t + λ)^{1 - a}} \cdot \sum_{i = 0}^{1} \sum_{s = 0}^{i + l} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | d λ \end{matrix}

\begin{matrix} + \frac{C_{2} r^{l - | β | - m}}{(t - r)^{l + m + a}} \int_{0}^{t - r} \frac{λ^{2 m}}{(t - r - λ)^{1 - a}} \cdot \sum_{i = 0}^{1} \sum_{s = 0}^{i + l} λ^{s + 1 - i} | f_{i}^{(s)} (λ) | d λ, \end{matrix}

\begin{matrix}  \end{matrix}

where

C_{2} = 0

for odd values of

n

. Once we now combine this with 4.1 , we get

\begin{matrix} | D^{β} u_{0} | & \leq C_{1} ɛ r^{l - | β | - m - a} \int_{t - r}^{t + r} \frac{{〈 λ 〉}^{m - l - k} d λ}{(r - t + λ)^{1 - a}} + \frac{C_{2} ɛ r^{l - | β | - m}}{(t - r)^{l + m + a}} \int_{0}^{t - r} \frac{λ^{l + m} {〈 λ 〉}^{m - l - k} d λ}{(t - r - λ)^{1 - a}} \end{matrix}

\begin{matrix} \equiv A_{1} + A_{2}, \end{matrix}

\begin{matrix}  \end{matrix}

where

A_{2} = 0

for odd values of

n

In what follows, we only deal with the decay rates

0 \leq k < 2 (m + a)

that correspond to parts (a) through (c), omitting the similar approach that one needs for parts (d) and (e).

Since

t + r

and

t - r

are equivalent whenever

t \geq 2 r

, an inequality of the form

\begin{matrix} | D^{β} u_{0} | \leq C ɛ r^{l - | β | - m} \cdot {〈 t + r 〉}^{m - l - k} \end{matrix}

(4.2)

implies each of the three desired estimates. One can easily obtain this upper bound for

A_{1}

since

λ

is equivalent to

t \pm r

within the region of integration. As for the integral term

A_{2}

, we may assume that

n

is even and thus focus on decay rates

0 \leq k < 2 m + 1

. Write

A_{2} = C_{2} ɛ r^{l - | β | - m} \cdot (t - r)^{- l - m - a} \cdot ℬ (l + m, m - l - k)

using the notation of Lemma 4.2 . Since

2 m - k > - 1

by above, the lemma ensures that

\begin{matrix} ℬ (l + m, m - l - k) \leq C (t - r)^{l + m + a} \cdot {〈 t - r 〉}^{m - l - k} . \end{matrix}

(4.3)

Once we now combine the last two equations, we may deduce the desired 4.2 .

Step 2:In the exterior region

t \leq 2 r

, Proposition 3.2 with

j = max (i, | β |)

similarly gives

\begin{matrix} | D^{β} u_{0} (r, t) | & \leq C_{1}^{'} ɛ r^{- m - a} \int_{| t - r |}^{t + r} \frac{λ^{l - | β |} {〈 λ 〉}^{m - l - k} d λ}{(r - t + λ)^{1 - a}} + C_{1}^{'} ɛ r^{- m - a} | t \pm r |^{a + l - | β |} {〈 t \pm r 〉}^{m - l - k} \end{matrix}

\begin{matrix} + \frac{C_{2}^{'} ɛ r^{- m - a}}{(t - r)^{m + | β |}} \int_{0}^{max (t - r, 0)} \frac{λ^{l + m} {〈 λ 〉}^{m - l - k} d λ}{(t - r - λ)^{1 - a}} \end{matrix}

\begin{matrix} \equiv B_{1} + B_{2} + B_{3}, \end{matrix}

\begin{matrix}  \end{matrix}

where

B_{3} = 0

for odd values of

n

. Let us only consider the decay rates

0 \leq k < 2 (m + a)

, as before. In view of 4.3 , the integral term

B_{3}

is then such that

B_{3} \leq C ɛ r^{- m - a} \cdot | t - r |^{a + l - | β |} \cdot {〈 t - r 〉}^{m - l - k} .

Since the right hand side appears in the boundary terms

B_{2}

, it thus suffices to treat

\begin{matrix} B_{1} + B_{2} = C_{1}^{'} ɛ r^{- m - a} [\int_{| t - r |}^{t + r} \frac{λ^{l - | β |} {〈 λ 〉}^{m - l - k} d λ}{(r - t + λ)^{1 - a}} + | t \pm r |^{a + l - | β |} {〈 t \pm r 〉}^{m - l - k}] . \end{matrix}

(4.4)

Case 1:If

t \leq 2 r

and

r \leq 1

, then

t + r

is bounded and we need only show that

| D^{β} u_{0} (r, t) | \leq C ɛ r^{l - | β | - m} .

In fact,

| t \pm r | \leq 3 r

for this case, so we easily get the desired estimate

B_{1} + B_{2} \leq C ɛ r^{l - | β | - m - a} [\int_{| t - r |}^{t + r} (r - t + λ)^{a - 1} d λ + r^{a}] \leq C ɛ r^{l - | β | - m}

because

| β | \leq l

and

a > 0

by assumption.Case 2:Suppose now that

t \leq 2 r

and

r \geq 1

. Since

| β | \leq l

by assumption, 4.4 leads to

\begin{matrix} B_{1} + B_{2} & \leq C ɛ r^{- m - a} [\int_{| t - r |}^{t + r} \frac{{〈 λ 〉}^{m - k - | β |} d λ}{(r - t + λ)^{1 - a}} + {〈 t \pm r 〉}^{m + a - k - | β |}] \end{matrix}

\begin{matrix} \leq C ɛ r^{- m - a} {〈 t - r 〉}^{- | β |} [\int_{| t - r |}^{t + r} \frac{{〈 λ 〉}^{m - k} d λ}{(r - t + λ)^{1 - a}} + {〈 t \pm r 〉}^{m + a - k}] \end{matrix}

\begin{matrix}  \end{matrix}

and we may use Lemma 4.1 to handle the integral.

Subcase 2a:When

0 \leq k < m + a

, we find that

B_{1} + B_{2} \leq C ɛ r^{- m - a} {〈 t - r 〉}^{- | β |} [r^{a} {〈 t + r 〉}^{m - k} + {〈 t \pm r 〉}^{m + a - k}] .

Moreover,

m + a - k

is positive, so we get

B_{1} + B_{2} \leq C ɛ r^{- m - a} {〈 t - r 〉}^{- | β |} {〈 t + r 〉}^{m + a - k} .

This does imply the desired estimate, as

r

is equivalent to

〈 t + r 〉

when

r \geq max (t / 2, 1)

Subcase 2b:When

k = m + a

, our previous approach applies verbatim, although an extra logarithmic factor is now included in the estimate that Lemma 4.1 provides.

Subcase 2c:When

m + a < k < 2 (m + a)

, Lemma 4.1 yields

\begin{matrix} B_{1} + B_{2} & \leq C ɛ r^{- m - a} {〈 t - r 〉}^{- | β |} [r^{a} {〈 t + r 〉}^{- a} {〈 t - r 〉}^{m + a - k} + {〈 t \pm r 〉}^{m + a - k}] \end{matrix}

\begin{matrix} \leq C ɛ r^{- m - a} {〈 t - r 〉}^{m + a - k - | β |} . \end{matrix}

\begin{matrix}  \end{matrix}

Invoking the equivalence of

r

with

〈 t + r 〉

, we may thus deduce the desired estimate.

Finally, the uniqueness assertion of our theorem can be established using the same energy argument as in [4] . Namely, the condition that [4] imposed was of the form

\begin{matrix} | \partial_{r} u_{0} (r, t) | + | \partial_{t} u_{0} (r, t) | = O (r^{- (n - 1) / 2 + δ}) as r \to 0 \end{matrix}

(4.5)

for some

δ > 0

. In our case, the estimates we just proved are such that

| \partial_{r} u_{0} (r, t) | + | \partial_{t} u_{0} (r, t) | = O (r^{l - 1 - m}) as r \to 0 .

However,

l \geq 1

and

m

is either

(n - 3) / 2

(n - 2) / 2

, so the condition 4.5 does hold.

Corollary 4.3. Using the assumptions and notation of Theorem 1.1 , one may improve the estimates provided for the decay rates

0 \leq k < m + a = (n - 1) / 2

of part (a) as follows.

(i) When $0 \leq k < m + a - l$ , we actually have $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t + r 〉}^{m - l - k} .$
(ii) When $k_{0} - 1 \leq k - m - a + l < k_{0}$ for some integer $1 \leq k_{0} \leq l$ , we have $| D^{β} u_{0} (r, t) | \leq C_{0} ɛ r^{l - | β | - m} \cdot {〈 t - r 〉}^{min (l - | β | - k_{0}, 0)} {〈 t + r 〉}^{m - l - k - min (l - | β | - k_{0}, 0)} .$

Proof. If either

t \geq 2 r

r \leq 1

, then

〈 t + r 〉

is equivalent to

〈 t - r 〉

and each of the desired estimates follows from part (a) of Theorem 1.1 . Let us now assume that

t \leq 2 r

and

r \geq 1

As in Case 2 of the previous proof, our task reduces to the estimation of

\begin{matrix} B_{1} + B_{2} = C_{1}^{'} ɛ r^{- m - a} [\int_{| t - r |}^{t + r} \frac{λ^{l - | β |} {〈 λ 〉}^{m - l - k} d λ}{(r - t + λ)^{1 - a}} + | t \pm r |^{a + l - | β |} {〈 t \pm r 〉}^{m - l - k}] . \end{matrix}

(4.6)

In what follows, we only concern ourselves with part (ii) because part (i) is easier to settle.

Given any

| t - r | \leq λ \leq t + r

and any

0 \leq k_{0} \leq l

, let us write

λ^{l - | β |} = λ^{min (l - | β |, k_{0})} \cdot λ^{max (l - | β | - k_{0}, 0)} .

Since

| β | \leq l

by assumption, the exponents on the right hand side are non-negative, so

\begin{matrix} λ^{l - | β |} & \leq {〈 λ 〉}^{min (l - | β | - k_{0}, 0) + k_{0}} \cdot (t + r)^{max (l - | β | - k_{0}, 0)} \end{matrix}

\begin{matrix} \leq {〈 λ 〉}^{k_{0}} \cdot {〈 t - r 〉}^{min (l - | β | - k_{0}, 0)} \cdot (3 r)^{max (l - | β | - k_{0}, 0)} \end{matrix}

\begin{matrix}  \end{matrix}

whenever

| t - r | \leq λ \leq t + r

and

t \leq 2 r

. Inserting this fact in 4.6 , we thus arrive at

\begin{matrix} B_{1} + B_{2} & \leq C ɛ r^{- m - a + max (l - | β | - k_{0}, 0)} \cdot {〈 t - r 〉}^{min (l - | β | - k_{0}, 0)} \times \end{matrix}

\begin{matrix} [\int_{| t - r |}^{t + r} \frac{{〈 λ 〉}^{m - l - k + k_{0}} d λ}{(r - t + λ)^{1 - a}} + {〈 t \pm r 〉}^{m - l - k + k_{0} + a}] . \end{matrix}

\begin{matrix}  \end{matrix}

For the decay rates of part (ii), we have

m - l - k + k_{0} > - a

, so Lemma 4.1 ensures that

\int_{| t - r |}^{t + r} \frac{{〈 λ 〉}^{m - l - k + k_{0}} d λ}{(r - t + λ)^{1 - a}} + {〈 t \pm r 〉}^{m - l - k + k_{0} + a} \leq C {〈 t + r 〉}^{m - l - k + k_{0} + a} .

Since

r

is equivalent to

〈 t + r 〉

when

r \geq max (t / 2, 1)

, the desired estimate follows easily.

Acknowledgement

The author would like to express his gratitude to Walter A. Strauss, under whose direction this research was conducted as part of the author's doctoral dissertation at Brown University. References

F. Asakura, Existence of a global solution to a semilinear wave equation with slowly decreasing initial data in three space dimensions, Comm. Partial Differential Equations, 11 (1986), pp. 1459–1487.
H. Kubo, On the critical decay and power for semilinear wave equations in odd space dimensions, Discrete Contin. Dynam. Systems, 2 (1996), pp. 173–190.
, Slowly decaying solutions for semilinear wave equations in odd space dimensions, Nonlinear Anal., 28 (1997), pp. 327–357.
H. Kubo and K. Kubota, Asymptotic behaviors of radially symmetric solutions of $□ u = | u |^{p}$ for super critical values $p$ in odd space dimensions, Hokkaido Math. J., 24 (1995), pp. 287–336.
, Asymptotic behaviors of radially symmetric solutions of $□ u = | u |^{p}$ for super critical values $p$ in even space dimensions, Japan. J. Math. (N.S.), 24 (1998), pp. 191–256.
K. Kubota, Existence of a global solution to a semi-linear wave equation with initial data of noncompact support in low space dimensions, Hokkaido Math. J., 22 (1993), pp. 123–180.
H. Lamb, On the theory of waves propagated vertically in the atmosphere, Proc. London Math. Soc. (2), 7 (1909), pp. 122–141.
F. Planchon, J. Stalker, and A. S. Tahvildar-Zadeh, $L^{p}$ estimates for the wave equation with the inverse-square potential, Discrete Contin. Dynam. Systems, 9 (2003), pp. 427–442.
M. Rammaha, Finite-time blow-up for nonlinear wave equations in high dimensions, Comm. Partial Differential Equations, 12 (1987), pp. 677–700.
H. Takamura, Blow-up for semilinear wave equations with slowly decaying data in high dimensions, Differential Integral Equations, 8 (1995), pp. 647–661.
K. Tsutaya, A global existence theorem for semilinear wave equations with data of noncompact support in two space dimensions, Comm. Partial Differential Equations, 17 (1992), pp. 1925–1954.
, Global existence and the life span of solutions of semilinear wave equations with data of noncompact support in three space dimensions, Funkcial. Ekvac., 37 (1994), pp. 1–18.

Brown University, Box 1917, Providence, RI 02912 E-mail address : petekara@math.brown.edu

Paschalis Karageorgis