4
$\begingroup$

Question:

I want to evaluate the decay estimate of the integral

$I^d(t; v) = \int_0^{\sqrt{d}\pi} dr \, r^{d-2} \int_0^\pi \sin(tr) e^{i\sqrt{d}vtr\cos\theta} \sin^{d-2}\theta \, d\theta $

for sufficiently large $ t $, where $ d \in \mathbb{Z}_{\ge 2} $ and $ v \in \mathbb{R}_{\ge 0}$.

It is hypothesized that the inequality

$ |I^d(t; v)| \le C t^{-\frac{d-1}{2}} $

may hold with constant $C$.

How can I evaluate this estimate?

update

The integral is a solution for discrete linear wave equation. I am interested in parallel structures between continuous and discrete system.

$\endgroup$
6
  • $\begingroup$ Could you give some context ? $\endgroup$ Commented Jul 31 at 10:01
  • $\begingroup$ @DavidLehavi I solve an initial value problem of discrete linear wave equation. The integral is given by the solution in polar coordinate. I want the decay estimate of it. $\endgroup$
    – Ko Hey
    Commented Jul 31 at 10:06
  • $\begingroup$ The claimed inequality cannot hold. Consider the special case when $v = 0$. The $\theta$ integral yields a positive constant (the integrand is positive), so you are reduced to considering the integral $\int_0^{\sqrt{d}\pi} r^{d-2} \sin(tr) ~\mathrm{d}r$. This can be explicitly evaluated, and you find that asymptotic decay is no better than $1/t$. But I am wondering whether you actually want the expression you wrote down, or whether you want to impose a radial cut-off. $\endgroup$ Commented Jul 31 at 16:52
  • $\begingroup$ @WillieWong Thank you. The condition $v > 0$ may need. I think the integrand is similar to a solution for the continuous linear wave equation. The difference is only integral range. The claimed inequality holds in the continuous case. $\endgroup$
    – Ko Hey
    Commented Jul 31 at 17:04
  • $\begingroup$ Your integral $I^d(t;v)$ is a continuous function in $v$ for any fixed $t$ and parameter $d$. The explicit values of $v = 0$ means that you cannot have a uniform in $v$ estimate of the form you claimed. The best you can do is to have the constant $C$ degenerating as $v \searrow 0$. On the other hand, if you impose a cut-off (replace $r^{d-2}$ by $r^{d-2} \phi(r)$ where $\phi:\mathbb{R}\to\mathbb{R}$ is smooth and compactly supported), then I think I can prove a uniform bound. $\endgroup$ Commented Jul 31 at 17:46

1 Answer 1

4
$\begingroup$

TL;DR

As stated the desired bound cannot hold. However, the bound holds when $v > 0$

  • for $d \leq 3$;
  • for $d > 3$, for a constant $C$ depending on $v$ that blows up as $v \to 0$;
  • for any $d$, for a uniform constant $C$ independent of $v$ if we assume additional cut-offs in $r$.

For convenience I will prove the claims for $d = 2k+1$ odd, with $k \geq 1$; the even cases requires a little bit more work but follow the same principles.

Re-normalization

First I wish to do some change of variables to make typing a bit easier. If we replace $\rho = \sqrt{d}\pi r$, and $\tau = t / \sqrt{d}{\pi}$, the integral can be written as

$$ \frac{1}{(\sqrt{d}\pi)^{d-1}}\int_0^1 \rho^{d-2} \mathrm{d}\rho \int_0^\pi \sin(\tau\rho) \exp(i \nu \tau \rho \cos\theta) (\sin\theta)^{d-2} ~\mathrm{d}\theta $$

(here we also set $\nu = \sqrt{d}v$). Our goal reduces to studying the integral

$$ J(\tau;\nu,k) := \int_0^1 \rho^{2k-1} \mathrm{d}\rho \int_0^\pi \sin(\tau\rho) \exp(i \nu \tau \rho \cos\theta) (\sin\theta)^{2k-1} ~\mathrm{d}\theta $$

It is convenient to take another change of variables $z = \cos\theta$ to write $$ J(\tau;\nu,k) = \int_0^1 \rho^{2k-1} \mathrm{d}\rho \int_{-1}^1 \sin(\tau\rho) \exp(i \nu \tau\rho z)(1 - z^2)^{k-1} ~\mathrm{d}z $$

$\nu = 0$

In this case the inner exponential drops, and we get no oscillation in $\theta$. Let $\gamma_k = \int_0^\pi (\sin\theta)^{2k-1} ~\mathrm{d}\theta > 0$, we have after integrating by parts $$ J(\tau;0,k) = \gamma_k \int_0^1 \rho^{2k-1} \sin(\tau\rho) ~\mathrm{d}\rho = \frac{\gamma_k}{\tau} \left( - \cos(\tau) + \int_0^1 (2k-1) \rho^{2k-2} \cos(\tau\rho) ~\mathrm{d}\rho \right) $$ Integrating by parts again we find $$ J(\tau;0,k) = -\frac{\gamma_k}{\tau}\cos(\tau) + O_{\tau\to\infty}(\tau^{-2}) $$

As $J(\tau;\nu,k)$ is continuous in $\nu$ for fixed $k$ and $\tau$, the asymptotics above shows that it is impossible to have any bound of the form $$ |J(\tau;\nu,k)| \leq \frac{C_k}{\tau^{k}} $$ with $C_k$ dependent on $k$ but independent of $\nu$, when $k > 1$.

Angular integration

We may write $$ \exp(i \nu \tau \rho z) = \frac{1}{i \nu \tau \rho} \frac{\mathrm{d}}{\mathrm{d}z}\exp(i \nu \tau \rho z) $$ Using that $(1-z^2)^{k-1}$ vanishes to order $k-1$ at $z = \pm 1$, we can integrate by parts $k$ times to get $$ J(\tau;\nu,k) = \int_0^1 \rho^{k-1} \left( \frac{-1}{i \nu \tau }\right)^{k} \sin(\tau\rho) \mathrm{d}\rho \left[ \int_{-1}^1 \exp(i \nu\tau \rho z) \frac{\mathrm{d}^k}{\mathrm{d}z^k} (1-z^2)^{k-1} ~\mathrm{d}z - \left. \exp(i \nu \tau\rho z) \frac{\mathrm{d}^{k-1}}{\mathrm{d}z^{k-1}} (1-z^2)^{k-1} \right|_{z = -1}^{z = 1} \right]$$

This shows that there exists a constant $C_k$ such that $$ | J(\tau;\nu,k) | \leq \frac{C_k}{\nu^k \tau^k} $$ which is the desired decay rate except for the fact that the rate depends on $\nu$.

Radial integration for small $\nu$

It turns out that the $\tau^{-1}$ decay at $\nu = 0$ can be extended to all $\nu$ small in a uniform way. We may write $$ \sin(\tau\rho) = \frac{1}{2i} (e^{i\tau\rho} - e^{-i\tau\rho}) $$ and so $$J(\tau;\nu,k) = J_+(\tau;\nu,k) - J_-(\tau;\nu,k)$$ where $$J_\pm(\tau;\nu,k) = \int_0^1 \rho^{2k-1} \mathrm{d}\rho \int_{-1}^1 \exp(i \tau\rho (\nu z \pm 1)) (1-z^2)^{k-1} ~\mathrm{d} z $$ Integrating by parts in $\rho$ we find $$J_\pm(\tau;\nu,k) = \int_{-1}^1 \frac{1}{i \tau (\nu z \pm 1)} (1 - z^2)^{k-1} \left[ \exp(i \tau(\nu z \pm 1)) - \int_0^1 (2k-1)\rho^{2k-2} \exp(i \tau\rho(\nu z \pm 1)) ~\mathrm{d}\rho \right] ~\mathrm{d} z $$ Using that $\nu \in [0,\frac12] \implies |\nu z \pm 1| \geq \frac12$ we conclude that there exists a constant $\tilde{C}_k$ dependent on $k$ but not on $\nu$ such that $$\nu \in [0,\frac12] \implies |J_\pm(\tau;\nu,k)| \leq \frac{\tilde{C}_k }{\tau} $$

Result for $d = 3$ ($k = 1$)

Apply the result from the angular integration to $\nu \geq \frac12$, and the result from the radial integration to $\nu \leq \frac12$, we obtain that

$$ |J(\tau;\nu,1)| \leq \frac{2\max( \tilde{C}_1,C_1)}{\tau} $$

as hypothesized.

Additional radial cut-off

We remark again that the analysis at $\nu = 0$ already rules out the hypothesized bound when $k > 1$. The issue is firmly with the radial integral having an amplitude $\rho^{d-2}$ that does not vanish at the upper limit $\rho = 1$.

If $\rho^{d-2}$ is replaced by $\phi(\rho)$ which vanishes to order at least $\kappa$ (with $\kappa\in \mathbb{Z}_{\geq 1}$) at both $\rho = 0$ and $\rho = 1$, then we would be able to repeat the radial integration by parts $\kappa + 1$ times (exactly as the angular integration), and obtain that for $\nu \in [0,\frac12]$ the upper bound $\tilde{\tilde{C}}_{\kappa,k} / \tau^{\kappa + 1}$

$\endgroup$

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .