I'm interested in diffusion, a.k.a. the heat kernel driven by the Laplace-Beltrami operator, on the $n$-dimensional sphere. There are lots of bounds showing that, for small times, it behaves in a way close to the heat kernel in $\mathbb{R}^n$: that is, the probability $p\_t(\theta)$ that we have moved an angle $\theta$ away from the starting point, at time $t$, is bounded by a Gaussian of variance $t$. But I need some control on how the errors or leading constants in these bounds depend on $n$.


The most precise estimates I've found (where the sphere has radius $1$) is due to Molchanov, 
$$
p\_t(\theta) \sim \frac{\mathrm{e}^{-\theta^2/2t}}{(2\pi t)^{n/2}} \left( \frac{\theta}{\sin \theta} \right)^{(n-1)/2}
$$
This is the first term of an asymptotic series, and if I understand his paper correctly, for fixed $n$ the next term would give a multiplicative error of $1+O(t)$. But does the constant hidden in $O(t)$ grow rapidly with $n$?


Another family of bounds gives
$$
p\_t(\theta) \le C \,t^{-n/2} \,\mathrm{e}^{-\theta^2/(4+\delta)t}
$$
Any constant value of $\delta$ would be fine with me. But we need $C$ to decay roughly as $(4\pi)^{-n/2}$ to match the normalization of the flat-space Gaussian.


All we really need for our application is the following. If we think of the heat kernel as a stochastic process, let $\theta$ be the angular distance from the initial point. We know that for small enough $t$ we have 
$$
\mathbb{E}(\theta) \le C' \sqrt{nt}
$$
for some constant $C'$, as it would be in flat space; but we need to know this holds for all $t$ up to $1/n$ or so. In other words, we need to know that the error term doesn't do something horrible like
$$
\mathbb{E}(\theta) = C' \sqrt{nt} + O(2^n t) \, . 
$$
This would follow, for instance, from a bound on the error term in Molchanov's estimate above, although this might be overkill.


Thanks!


Cris



heat kernel on n-sphere
Yes, it is true that $\theta$ on $S\_n$ is dominated by $\theta$ on $\mathbb{R}^n$.


Let $(B\_t)$ be a Brownian motion on the sphere. The radial process $\theta\_t=d(x,B\_t)$ is a Jacobi process, that is a Markov process with generator


$
L=\frac{n-1}{2} \text{cotan} (r) \frac{d}{dr} +\frac{1}{2} \frac{d^2}{dr^2 }
$


Since $ \text{cotan} (r) \le \frac{1}{r}$, we deduce from the comparison theorem for stochastic differential equations that


$
\theta\_t \le \beta\_t
$


where $\beta\_t$ is a Bessel process, that is a Markov process with generator


$
L=\frac{n-1}{2r} \frac{d}{dr} +\frac{1}{2} \frac{d^2}{dr^2 }
$


The Bessel process is the radial part of the Euclidean Brownian motion in $\mathbb{R}^n$, so you get the desired stochastic domination.


