---
abstract: |
    Consider the following stochastic heat equation, $$\begin{aligned}
    \frac{\partial u_t(x)}{\partial t}=-\nu(-\Delta)^{\alpha/2} u_t(x)+\sigma(u_t(x))\dot{F}(t,\,x), \quad t>0, \; x \in {\mathbf{R}}^d.\end{aligned}$$ Here $-\nu(-\Delta)^{\alpha/2}$ is the fractional Laplacian with $\nu>0$ and $\alpha \in (0,2]$, $\sigma: {\mathbf{R}}\rightarrow {\mathbf{R}}$ is a globally Lipschitz function, and $\dot{F}(t,\,x)$ is a Gaussian noise which is white in time and colored in space. Under some suitable additional conditions, we prove a strong comparison theorem and explore the effect of the initial data on the spatial asymptotic properties of the solution. This constitutes an important extension over a series of works most notably [@CJK], [@CJKS2], [@CK] and [@CDKh].\
    [*Keywords:*]{} Stochastic PDEs, comparison theorems, colored noise.\
    [*AMS 2010 subject classification:*]{} Primary 60H15; Secondary: 35K57.
author:
- |
    Mohammud Foondun\
    University of Strathclyde
- |
    Eulalia Nualart\
    Pompeu Fabra University\
bibliography:
- 'Foon.bib'
title: 'Spatial asymptotics and strong comparison principle for some fractional stochastic heat equations.'
---

Introduction and main results
=============================

Consider the following stochastic heat equation, $$\begin{aligned}
\label{main-eq}
\frac{\partial u_t(x)}{\partial t}=-\nu(-\Delta)^{\alpha/2} u_t(x)+\sigma(u_t(x))\dot{F}(t,\,x),  \quad t>0, \; x \in {\mathbf{R}}^d,\end{aligned}$$ where $-\nu(-\Delta)^{\alpha/2}$ is the fractional Laplacian, that is, the infinitessimal generator of a symmetric $\alpha$-stable process with density $p_t(x)$, where $\alpha \in (0,2]$, and $\nu>0$ is a viscosity constant. The noise $\dot{F}(t,\,x)$ is white in time and colored in space satisfying $$\operatorname{Cov}(\dot{F}(t,\,x), \dot{F}(s,\,y))=\delta_0(t-s)f(x-y),$$ where $f$ is the spatial correlation function which we take to be the Riesz kernel $$\begin{aligned}
f(x):=\frac{1}{|x|^\beta},\qquad 0<\beta<d. \end{aligned}$$ The function $\sigma:{\mathbf{R}}\rightarrow {\mathbf{R}}$ is a globally Lipschitz continuous function with $\sigma(0)=0$, that is, there exists a constant $L_{\sigma}>0$ such that $$\vert \sigma(x)\vert {\leqslant}L_{\sigma} \vert x\vert,\quad \text{  for all } x \in {\mathbf{R}}^d.$$ The initial condition $u_0$ is always going to be a nonnegative function in ${\mathbf{R}}^d$ such that $$\bar{u}_0:=\sup_{x\in {\mathbf{R}}^d}u_0(x)<\infty.$$ Following Walsh [@Walsh], if one further assume that $$\begin{aligned}
\beta<\min(\alpha,\,d),\end{aligned}$$ then has a unique mild solution $\{u_t(x), t {\geqslant}0, x \in {\mathbf{R}}^d\}$ which is adapted, jointly measurable and satisfies $$\label{mild}
u_t(x)=(p_t\ast u_0)(x)+\int_0^t\int_{{\mathbf{R}}^d}p_{t-s}(x-y)\sigma(u_s(y))\,F({{\rm d}}s\,{{\rm d}}y),$$ where $$(p_t\ast u_0)(x)=\int_{{\mathbf{R}}^d} p_t(x-y) u_0(y) dy,$$ and $$\sup_{x\in {\mathbf{R}}^d, t\in [0,\,T]}\E|u_t(x)|^k<\infty\quad\text{for all}\quad k{\geqslant}2\quad\text{and}\quad T<\infty.$$ For more information about existence-uniqueness considerations, consult [@Walsh], [@minicourse] and [@Davar-CBMS]. This paper is motivated by two important results proved recently in [@FLJ]. The first one is the following weak comparison principle.

\[weakc\][[@FLJ]]{.nodecor} Suppose that $u$ and $v$ are two solutions to with initial conditions $u_0$ and $v_0$ respectively such that $u_0{\leqslant}v_0$. Then $${\mathrm{P}}(u_t(x){\leqslant}v_t(x)\quad \text{for all}\quad x\in {\mathbf{R}}^d,\quad t{\geqslant}0)=1.$$

Theorem \[weakc\] ensures nonnegativity of the solution, since the initial condition is assumed to be nonnegative. For the special case $\sigma(x)=x$ (known as the Parabolic Anderson model), this fact can be deduced from the Feynman-Kac representation of the solution. However, for the general non-linear case, this property for the solution to (\[mild\]) was unknown until the work of [@FLJ].

The first aim of this paper is to use Theorem \[weakc\] in order to show the following strong comparison principle.

\[strongcomparison\] Suppose that $u$ and $v$ are two solutions to with initial conditions $u_0$ and $v_0$ respectively such that $u_0< v_0$. Assume $\alpha {\geqslant}1$. Then $${\mathrm{P}}(u_t(x)< v_t(x)\quad \text{for all}\quad x\in {\mathbf{R}}^d,\quad t{\geqslant}0)=1.$$

The (strong) comparison principle for equation (\[main-eq\]) with space-time white noise and $\alpha=2$ is the well-known Mueller’s comprison principle (see [@Mueller]). Recently, several extensions have been developed. In [@CK] the authors extend Mueller’s result when the initial data is more general and there is a more general fractional differential operator than the fractional Laplacian. In [@ChenHuang] the authors consider the non-linear heat equation in ${\mathbf{R}}^d$ with a general spatial covariance and measured-valued initial data. The proof of our strong comparison principle uses the same strategy as in the papers mentioned above. But the presence of the fractional Laplacian and the colored noise makes it that we have to work a bit harder to prove our result. For the sake of conciseness, we only consider the Riesz kernel spatial covariance, but we believe that our method could be extended to general spatial covariances as in [@ChenHuang].

As another consequence of the weak comparison principle (Theorem \[weakc\]), we show the next quantitative result on the strict positivity of the solution, which is an extension of [@CJKCorr Theorem 5.1] (space-time white noise and $\alpha=2$). See also [@CK Theorem 1.4] and [@ChenHuang Theorem 1.6]. Note that $\alpha$ is not required to be bigger than $1$.

\[positivity\] Let $T>0$ and $K\subset {\mathbf{R}}^d$ be a compact set contained in the support of the initial condition $u_0$. Then, there exist constants $c_1$ and $c_2$ depending on $T$ and $K$ such that for all $\epsilon>0$, we have

$${\mathrm{P}}\left( \inf_{t\in [0,\,T]}\inf_{x\in K}u_t(x) < \epsilon \right) {\leqslant}c_2 \exp \left( -c_1 \{\vert \log \epsilon \vert  \log \vert \log \epsilon \vert\}^{\frac{2\alpha-\beta}{\alpha}}\right).$$

Let us now state the second motivation of this paper, which is the following moment comparison theorem.

\[mct\][[@FLJ]]{.nodecor} Let $u$ and $v$ two solutions to , the first one with $\sigma$, the other with another globally Lipschitz continuous function $\bar{\sigma}$ such that $\bar{\sigma}(0)=\sigma(0)=0$ and $\sigma(x) {\geqslant}\bar{\sigma}(x){\geqslant}0$ for all $x\in {\mathbf{R}}_+$. Then for any $k \in \mathbb{N}$, $x \in {\mathbf{R}}^d$, and $t {\geqslant}0$, $$\E[u_t(x)^k] {\geqslant}\E[v_t(x)^k].$$

An important consequence of this result are the following sharp estimates on the moments of the solution to (\[mild\]), when the initial condition is bounded below and under the additional assumption that there exists a constant $l_\sigma>0$ such that $$\begin{aligned}
 \label{sig}
\quad  \sigma(x){\geqslant}l_\sigma|x|,\quad \text{for all}\quad x\in {\mathbf{R}}^d.\end{aligned}$$ This was unknown until the work of [@FLJ].

\[momentbounds\] Let $u$ be the solution to [(\[mild\])]{.nodecor}. Assume [(\[sig\])]{.nodecor} and $$\label{infsup}
0<\underline{u}_0:=\inf_{x\in {\mathbf{R}}^d}u_0(x).$$ Then there exists a positive constant $A$ such that for all $x\in {\mathbf{R}}^d$, $t>0$, and $k {\geqslant}2$, $$\frac{\underline{u}_0^k}{A^{k}} \exp\left(\frac{1}{A}k^{\frac{2\alpha-\beta}{\alpha-\beta}}t \nu^{-\frac{\beta}{\alpha-\beta}}\right){\leqslant}\E|u_t(x)|^k{\leqslant}A^k\bar{u}_0^k\exp\left(Ak^{\frac{2\alpha-\beta}{\alpha-\beta}}t\nu^{-\frac{\beta}{\alpha-\beta}}\right).$$

For the Parabolic Anderson model, the above is given by [@KK Lemma 4.1]. The scaling property of the heat kernel gives the dependence of the bounds on the parameter $\nu$. An immediate consequence of Theorem \[momentbounds\] is that the solution to (\[main-eq\]) is fully intermittent meaning that for all $k {\geqslant}2$, the function $$k  \rightarrow \frac{1}{k} \gamma(k):=\frac{1}{k}\limsup_{t \rightarrow \infty} \log \E|u_t(x)|^k
\quad \text{ is strictly increasing}.$$ Intuitively, this means that the solution develops many high peaks distributed over small $x$-intervals when $t$ is large (see [@FK] and the references therein). The fact that the solution to (\[main-eq\]) is intermittent was already known (see e.g.[@FLO] and the references therein) but was shown by showing $$\gamma(2)>0 \quad \text{ and } \quad \gamma(k)<\infty, \text{ for all } k {\geqslant}2.$$ The previous results concern the moments of the solution to (\[main-eq\]), but much less is known about the almost sure asymptotic behaviour of the solution, which is crucial to understand better its chaotic behaviour. This brings us to the second purpose of this paper is to explore how the almost surely spatial asymptotic behaviour of the solution to (\[main-eq\]) depends on the initial function $u_0$. We start with the case that $u_0$ is bounded below as in Theorem \[momentbounds\]. A first observation is that, since $u_0$ is also bounded above, then we can easily see that $$\E u_t(x){\leqslant}c,$$ where $c$ is the upper bound of $u_0$. Since $u_0$ is bounded below, it is not trivial to say more about $\liminf_{|x|\rightarrow \infty}u_t(x)$ other than it is almost surely bounded above. This is in sharp contrast with the $\limsup$ behaviour of the solution described by the next theorem.

\[asym\] Let $u$ be the unique solution to [(\[mild\])]{.nodecor}, and assume that [(\[sig\])]{.nodecor} and [(\[infsup\])]{.nodecor} hold. Then there exist positive constants $c_1,c_2$ such that for every $t > 0$, $$\begin{aligned}
c_1\frac{t^{(\alpha-\beta)/(2\alpha-\beta)}}{\nu^{\beta/(2\alpha-\beta)}} {\leqslant}\lim_{R \rightarrow \infty}\frac{\log \sup_{x\in B(0,\,R)}u_t(x)}{\left(\log R\right)^{\alpha/(2\alpha-\beta)}}{\leqslant}c_2 \frac{t^{(\alpha-\beta)/(2\alpha-\beta)}}{\nu^{\beta/(2\alpha-\beta)}}\quad\text{a.s.}\end{aligned}$$

This theorem is a major improvement of [@CJK Theorem 1.3] (space-time white noise case) and [@CJKS2 Theorem 2.6] (Riesz kernel spatial covariance). See also [@XC] for exact spatial asymptotics when then noise is fractional in time and correlated in space. All these papers deal with the Parabolic Anderson model and the usual Laplacian ($\alpha=2$). Moreover, in [@CJK; @CJKS2] the dependence in time of the bounds is not explicit. The case $\sigma(x)=x$, fractional Laplacian and Riesz kernel spatial covariance is considered in the preprint [@KK Theorem 1.2], without the dependence on $\nu$ and constant initial data. Obtaining the exact dependence on the viscosity constant $\nu$ is important to understand in which universality class the equation can be associated to. See [@CJK Remark 1.5].

A key ingredient of the proof of Theorem \[asym\] are the moment bounds of Theorem \[momentbounds\], that will allow to obtain some tail estimates for the solution. Let us now consider the next observation where $u_0$ is not bounded below.

If $u_0(x):=1_{B(0,1)}(x)$, then one can show that for $x\in B(0,\,R)^c$ and $R$ large enough, we have $$\E u_t(x)=(p_t\ast u_0)(x){\leqslant}\frac{ct}{R^{d+\alpha}}.$$ Then, a Borel-Cantelli argument shows that almost surely $$\lim_{|R|\rightarrow \infty}\inf_{B(0,\,R^c)}u_t(x)=0.$$ This motivated our next result.

The above remark can be seen as a motivation for us to drop the assumption that the initial function is bounded below. We have the following trichotomy result, that studies the amount of decay that the initial conditions needs to ensure that the solution is a bounded function a.s.

\[trichotomy\] Let $u$ be the unique solution to [(\[mild\])]{.nodecor}. Assume [(\[sig\])]{.nodecor} and that $u_0(x)$ is a radial function satisfying $$\lim_{x\rightarrow \infty}u_0(x)=0\quad \text{and}\quad u_0(x){\leqslant}u_0(y)\quad\text{whenever}\quad x{\geqslant}y.$$ Set $$\Lambda:=\lim_{|x|\rightarrow \infty}\frac{|\log u_0(x)|}{(\log |x|)^{\alpha/(2\alpha-\beta)}}.$$ Then, if $0<\Lambda<\infty$, there exists a random variable $T$ such that $$\begin{aligned}
{\mathrm{P}}\left(\sup_{x\in {\mathbf{R}}^d} u_t(x)<\infty, \quad \forall t<T\quad \text{and}\quad \sup_{x\in {\mathbf{R}}^d} u_t(x)=\infty, \quad \forall t>T\right)=1.\end{aligned}$$ Moreover, if $\Lambda=\infty$, then $$\begin{aligned}
{\mathrm{P}}\left(\sup_{x\in {\mathbf{R}}^d} u_t(x)<\infty, \quad \forall t>0\right)=1.\end{aligned}$$ Finally, if $\Lambda=0$, then $$\begin{aligned}
{\mathrm{P}}\left(\sup_{x\in {\mathbf{R}}^d} u_t(x)=\infty, \quad \forall t>0\right)=1.\end{aligned}$$

This result is a major extension of [@CDKh Theorem 1.1], where the case $\alpha=2$ and space-time white noise is considered. The proof of their result is based on the technical Lemma [@CDKh Lemma 2.3] which follows the ideas of [@Chen-Dalang]. Here, we improve the method of the proof, and our insensitivity theorem is based on a Gronwall’s type result (see Proposition \[renew\] below). The latter result is one of the technical innovations of this paper. Theorem \[trichotomy\] shows precisely the effect of the fractional Laplacian and a smoother noise. In fact we can make the following observation.

We assume $\beta=1$. Observe that for $\alpha<2$, $\frac{2}{3}<\frac{\alpha}{2\alpha-1}$. Choose $\epsilon$ such that $\frac{2}{3}<\epsilon<\frac{\alpha}{2\alpha-1}$, and $u_0$ such that $|\log u_0(x)| \sim (\log |x|)^{\epsilon}$ as $|x|\rightarrow \infty$. Then $\Lambda=\infty$ if $\alpha=2$ but $\Lambda=0$ if $\alpha<2$. Thus, for the same initial data and noise, the solution is bounded for the usual Laplacian but unbounded for the fractional Laplacian.

Observe that when $u_0$ has compact support corresponds to the case where $\Lambda=\infty$, and Theorem \[trichotomy\] shows that the solution is bounded for all times a.s.

We now give a plan of the article. In Section 2 we give some preliminary results needed throughout the paper. Section 3 is devoted to an approximation result needed for the proof of Theorems [\[asym\]]{.nodecor} and [\[trichotomy\]]{.nodecor}. These theorems are proved in Section 4. Finally, Section 5 gives the proof of Theorems [\[strongcomparison\]]{.nodecor} and [\[positivity\]]{.nodecor}.

Preliminary results
===================

Let $X_t$ be the symmetric $\alpha$-stable process associated with the fractional Laplacian $-\nu (-\Delta)^{\alpha/2}$ and let $p_t(x)$ denote its heat kernel. We will frequently use the following properties.

-   [*Scaling property*]{}: For any positive constant $a$, we have $$\begin{aligned}
    p_t(x)=a^dp_{a^\alpha t}(ax), \quad \text{ for all } x \in {\mathbf{R}}^d, t >0. \end{aligned}$$ This property follows from $$p_t(x)=(2\pi)^{-d}\int_{{\mathbf{R}}^d}e^{-ix\cdot z}e^{-t\nu|z|^\alpha}\,{{\rm d}}z.$$

-   [*Heat kernel estimates*]{} (see [@kolo] and references therein): For $0<\alpha<2$, there exist positive constants $c_1$ and $c_2$ such that for all $x \in {\mathbf{R}}^d$ and $t>0$, $$\begin{aligned}
    c_1\left(\frac{1}{t^{d/\alpha}}\wedge\frac{t}{|x|^{d+\alpha}}\right){\leqslant}p_t(x){\leqslant}c_2\left(\frac{1}{t^{d/\alpha}}\wedge\frac{t}{|x|^{d+\alpha}}\right).\end{aligned}$$

The proofs of Theorems [\[asym\]]{.nodecor} and [\[trichotomy\]]{.nodecor} will only use the upper bound $$\label{ee2}
p_t(x) {\leqslant}c_2\frac{t}{\vert x\vert^{d+\alpha}},\quad  \text{ for sufficiently large } \vert x \vert,$$ while the proof of Theorems [\[strongcomparison\]]{.nodecor} and [\[positivity\]]{.nodecor} will only use the lower bound $$\label{ee1}
p_t(x) {\geqslant}c_1 \frac{1}{t^{d/\alpha}},\quad  \text{ for sufficiently small } \vert x \vert.$$ Both are also valid for $\alpha=2$.

The next result provides some estimates that involve the above heat kernel and the correlation function $f$. These estimates will be useful for the proof of our ‘insensitivity’ result; see Theorem \[sensitivity\].

\[heat-corr\] There exist positive constants $c_1,c_2$ and $c_3$ such that for all $t>0$, $x\in {\mathbf{R}}^d$, and $R>0$, we have $$\begin{aligned}
 \label{a1}
\int_{B(x,\,R)^c\times B(x,\,R)^c}p_{t}(x-y)p_{t}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w &{\leqslant}c_1 \frac{t^2}{R^{2\alpha+\beta}},\\ \label{a2}
\int_{B(x,\,R)^c\times B(x,\,R)}p_{t}(x-y)p_{t}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w&{\leqslant}c_2 \frac{t^{1-\beta/\alpha}}{R^\alpha},\\ \label{a3}
\int_{{\mathbf{R}}^d\times {\mathbf{R}}^d}p_{t}(x-y)p_{t}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w&{\leqslant}c_3 t^{-\beta/\alpha}.\end{aligned}$$

We start with (\[a1\]). $$\begin{aligned}
\int_{B(x,\,R)^c\times B(x,\,R)^c}&p_{t}(x-y)p_{t}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w\\
&{\leqslant}\int_{B(0,\,R)^c\times B(0,\,R)^c}p_t(y)p_t(w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w.\end{aligned}$$ From (\[ee2\]), the above quantity is bounded by a constant times $$\begin{aligned}
\frac{t^2}{R^{2\alpha+\beta}}\int_{B(0,\,1)^c\times B(0,\,1)^c}\frac{1}{\vert y\vert^{d+\alpha}\vert w\vert^{d+\alpha}}\frac{1}{\vert
y-w\vert^\beta}\,{{\rm d}}w\,{{\rm d}}y.\end{aligned}$$ The above integral is finite so the proof of (\[a1\]) is complete. For (\[a2\]), we write $$\begin{aligned}
\int_{B(x,\,R)^c\times B(x,\,R)}&p_{t}(x-y)p_{t}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w\\
&{\leqslant}\int_{B(0,\,R)^c\times {\mathbf{R}}^d}p_t(y)p_t(w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w.\end{aligned}$$ By the scaling property, $$\begin{aligned}
\int_{{\mathbf{R}}^d}p_t(w)f(y-w)\,{{\rm d}}w=\E^y|X_t|^{-\beta}{\leqslant}t^{-\beta/\alpha}\E^0|X_1|^{-\beta}.\end{aligned}$$ Finally, proceeding as before, we get $$\int_{B(0,\,R)^c}p_t(y) dy {\leqslant}c \frac{t}{R^\alpha}.$$ Combining the above estimates, we obtain (\[a2\]). For (\[a3\]), it suffices to use the semigroup property $$\begin{aligned}
\int_{{\mathbf{R}}^d\times {\mathbf{R}}^d}p_{t}(x-y)p_{t}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w=\int_{{\mathbf{R}}^d}p_{2t}(w)f(w)\,{{\rm d}}w,\end{aligned}$$ and using the scaling property as before we obtain the desired bound.

We now return to $u_t(x)$ the solution to [(\[mild\])]{.nodecor}. Our next property can be read from [@BEM2010a]. For any $k {\geqslant}2$, there exists a positive constant $c:=c(k)$ such that for all $s,t>0$, and $x,\,y\in{\mathbf{R}}^d $ $$\begin{aligned}
\E|u_s(x)-u_t(y)|^k{\leqslant}c\left(\vert x-y\vert^{\eta k}+|s-t|^{\tilde{\eta}k}\right),\end{aligned}$$ where $\eta=\frac{\alpha-\beta}{2}\wedge \frac{1}{2}$ and $\tilde{\eta}=\frac{\alpha-\beta}{2\alpha}$. The above together with the upper moment bound of Theorem \[momentbounds\] has the following consequence.

\[cty\] $u_t(x)$ has a continuous version, that is, for any $k {\geqslant}2$, there exist positive constants $c_1,c_2:=c_2(k)$ such that $$\begin{aligned}
\E\left[\sup_{x\not=y, s\not= t}\frac{|u_s(x)-u_t(y)|^k}{\vert x-y\vert^{\eta k}+|s-t|^{\tilde{\eta}k}}\right]{\leqslant}c_2 e^{c_1k^{(2\alpha-\beta)/(\alpha-\beta)} \nu^{-\beta/(\alpha-\beta)}t}.\end{aligned}$$

The proof is very similar to Theorem 4.3 of [@minicourse] and is therefore omitted.

We also have the following property.

\[markov\] Fix $x\in {\mathbf{R}}^d$, then the solution $u_t(x)$ satisfies the strong Markov Property.

We omit the proof since it is very similar to [@Mueller Lemma 3.3].

As mentioned earlier, a key ingredient of the proof of Theorem \[trichotomy\] is an [*insensitivity*]{} theorem; see Theorem \[sensitivity\]. Its proof hinges on the following proposition which is one of the main technical innovations of this paper. The proof follows that of [@Henry Lemma 7.1.1].

\[renew\] Suppose that for $R>0$, the function $f_R(\cdot)$ is a non-negative non-decreasing locally integrable function on $[0,\,T]$ satisfying the following $$\begin{aligned}
f_R(t){\leqslant}A_R(t)+B\int_0^t \frac{f_{2R}(s)}{(t-s)^{\gamma}}\,\,ds,\end{aligned}$$ where $A_R(\cdot)$ is also a locally integrable non-decreasing function $[0,\,T]$, $B$ is a positive constant and $\gamma<1.$ If $$\begin{aligned}
\sup_{n{\geqslant}1}f_{2^nR}(t)<\infty \quad \text{and}\quad A_{(n+1)R}(t){\leqslant}A_{nR}(t)\quad \text{for}\quad n{\geqslant}1,\end{aligned}$$ then there exist positive constants $c_1,c_2$ such that for all $t \in [0,T]$, $$\begin{aligned}
f_R(t){\leqslant}c_2 A_R(t)e^{c_1t}.\end{aligned}$$

We begin by defining $$\begin{aligned}
L\phi(t):=B\int_0^t\frac{\phi(s)}{(t-s)^\gamma}\,{{\rm d}}s\end{aligned}$$ and $L^{(n)}\phi(t):=L(L^{(n-1)}\phi(t))$ for $n{\geqslant}1$. The inequality stated in the statement of the proposition can therefore be written as $$\begin{aligned}
f_R(t){\leqslant}A_R(t)+Lf_{2R}(t).\end{aligned}$$ Iterating the above and using the facts that $A_R(s)$ is non-decreasing in $s$ and satisfies $A_{nR}(s){\leqslant}A_{(n-1)R}(s)$, we have $$\begin{aligned}
f_R(t){\leqslant}A_{R}(t)\sum_{k=0}^{n-1} L^{(k)}1(t)+R_nf_R(t),\end{aligned}$$ where $1(t):=1$. $R_nf_R(t)$ is the remainder term given by $$R_nf_R(t):= \frac{1}{\Gamma(n(1-\gamma))}\int_0^t [B\Gamma((1-\gamma))]^n(t-s)^{n(1-\gamma)-1} f_{2^nR}(s)\,{{\rm d}}s.$$ Since $\sup_{n{\geqslant}1}f_{2^nR}(t)<\infty$, we have that as $n\rightarrow \infty$, $R_nf_R(t)\rightarrow 0$. Some computations show that

$$\begin{aligned}
 \sum_{k=0}^{\infty} L^{(k)}1(t)&=\int_0^t\sum_{k=0}^\infty\frac{(B\Gamma(1-\gamma))^k(t-s)^{k(1-\gamma)-1}}{\Gamma(k(1-\gamma))}\,{{\rm d}}s\\
&{\leqslant}c_2 e^{c_1t},\end{aligned}$$

where $c_1:=(B\Gamma(1-\gamma))^{1/(1-\gamma)}.$ Putting these estimates together, we obtain the desired result.

An approximation result
=======================

Theorems [\[asym\]]{.nodecor} and [\[trichotomy\]]{.nodecor} are almost sure limit theorems and rely on some Borel-Cantelli type arguments. To be able to carry out the proof, we will need to find an appropriate independent sequence of random variables and it is apriori not clear how to find such a sequence. We follow [@CJK] and [@CJKS2] where this issue was successfully resolved.

Let $n{\geqslant}1$ and consider the following approximation $F^{(n)}$ of the measure $F$ appearing in (\[mild\]). Recall that the covariance of $\dot{F}$ is given by $f(x)=\frac{1}{|x|^\beta}$, $x \in {\mathbf{R}}^d$. This can be written as $f=h\ast\tilde h$, where $h(x)=\frac{1}{|x|^{\frac{d+\beta}{2}}}$ and $\tilde{h}(x):=h(-x)$. Define $h_n(x):=h(x)Q_n(x)$ and $f_n(x)=(h-h_n)\ast(\tilde{h}-\tilde{h}_n)$, where $$Q_n(x)=\prod_{j=1}^d \left(1-\frac{\vert x_j \vert}{n}\right)_+.$$ We take $\dot F^{(n)}$ to be the noise satisfying $$\operatorname{Cov}(\dot{F}^{(n)}(t,\,x), \dot{F}^{(n)}(s,\,y))=\delta_0(t-s)g_n(x-y),$$ where $g_n=h_n \ast \tilde{h}_n$.

By an argument similar to that of the proof of [@CJKS2 Lemma 9.3], we have that for any $\gamma\in(0,\,\beta \wedge 1)$ there exists a positive constant $c$ such that for all $s>0$ and $n {\geqslant}1$, $$(p_s\ast f_n)(0){\leqslant}c \frac{1}{n^{\gamma}}\frac{1}{s^{(\beta-\gamma)/\alpha}}.$$ As a consequence, using the semigroup property, we obtain that $$\label{rty} \begin{split}
 &\int_0^t\int_{{\mathbf{R}}^d\times {\mathbf{R}}^d}p_{t-s}(x-y)p_{t-s}(x-z)
f_n(y-z)\, {{\rm d}}y \, {{\rm d}}z \, {{\rm d}}s \\ 
&=\int_0^t\int_{{\mathbf{R}}^d} p_{2(t-s)}(z)
f_n(z)\,  {{\rm d}}z \,  {{\rm d}}s =\int_0^t (p_{2r} \ast f_n)(0) {{\rm d}}r \\
&{\leqslant}c t^{1-\frac{\beta-\gamma}{\alpha}} \frac{1}{n^{\gamma}}.
\end{split}$$

Next, consider the following integral equation, $$\label{utn}
U_t^{(n)}(x)=(p_t\ast u_0)(x)+\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)\sigma(U_s^{(n)}(y))F^{(n)}({{\rm d}}s\,{{\rm d}}y).$$ The unique solution to this integral equation can be found via a standard fixed point argument. Fix $n{\geqslant}1$. Set $U_t^{(n, 0)}:=u_0$ and for each $j{\geqslant}1$, the $j$th Picard iteration is given by $$\begin{aligned}
U_t^{(n, j)}(x)=\int_{{\mathbf{R}}^d}p_t(x-y)u_0(y)\,{{\rm d}}y+\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)\sigma(U_s^{(n, j-1)}(y))F^{(n)}({{\rm d}}s\,{{\rm d}}y).\end{aligned}$$ Moreover, one can show that under our current standing conditions, the unique solution satisfies for all $t>0$ and $k {\geqslant}2$, $$\sup_{n {\geqslant}1}\sup_{s\in [0,\,t]} \sup_{x \in {\mathbf{R}}^d}\E|U_s^{(n)}(x)|^k{\leqslant}c_2 e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t},$$ for some positive constants $c_1,c_2(k)$. As a consequence, for all $t>0$, $k {\geqslant}2$, and sufficiently large $n$, $$\label{usup}
\sup_{s\in [0,\,t]} \sup_{x \in {\mathbf{R}}^d}\E|U_s^{(n,n-1)}(x)|^k{\leqslant}c_2 e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t},$$ for some positive constants $c_1,c_2(k)$. We also have the following result which gives us the independent quantities we need.

\[ind\] Let $t>0$ and $n{\geqslant}1$. Suppose that $\{x_i \}_{i=1}^\infty\subset {\mathbf{R}}^d$ with $\vert x_i-x_j\vert{\geqslant}2n^{1+1/\alpha}t^{1/\alpha}$ for all $i\not=j$. Then $\{U^{(n,n)}_t(x_i)\}_{i=1}^\infty$ are independent random variables.

The proof is similar to that of [@CJKS2 Lemma 5.4] and is omitted.

We will also need the fact that the random variables defined above approximate the solution to . We provide a proof of this next. Recall that $\gamma<\beta\wedge 1.$

\[approx\] For all $T >0$ and $k {\geqslant}2$, there exist positive constants $c_1$ and $c_2$ such that for large enough $n$, $$\begin{aligned}
\sup_{t \in [0,T]} \sup_{x \in {\mathbf{R}}^d}\E|u_t(x)-U_t^{(n, n)}(x)|^k{\leqslant}c_2  \frac{1}{n^{\gamma k/2}}e^{c_1k^{(2\alpha-\beta)/(\alpha-\beta)}t}.\end{aligned}$$

Consider the following integral equation, $$V_t^{(n)}(x)=(p_t\ast u_0)(x)+\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)\sigma(V_s^{(n)}(y))F({{\rm d}}s\,{{\rm d}}y).$$ We first look at $V^{(n)}_t(x)-U_t^{(n, n)}(x)$ and its moments. $$\begin{aligned}
V^{(n)}_t(x)-U_t^{(n,n)}(x)&=\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)\sigma(V_s^{(n)}(y))F({{\rm d}}s\,{{\rm d}}y)\\
&-\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)\sigma(U_s^{(n,n-1)}(y))F^{(n)}({{\rm d}}s\,{{\rm d}}y).\end{aligned}$$ We rewrite the above as $$\begin{aligned}
V^{(n)}_t(x)-U_t^{(n,n)}(x)&=\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)[\sigma(V_s^{(n)}(y))-\sigma(U_s^{(n,n-1)}(y))]F({{\rm d}}s\,{{\rm d}}y)\\
&-\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)\sigma(U_s^{(n,n-1)}(y))[F^{(n)}({{\rm d}}s\,{{\rm d}}y)-F({{\rm d}}s\,{{\rm d}}y)]\\
&:=I_1+I_2.\end{aligned}$$ We start by bounding $I_2$. Using Burkholder-Davis-Gundy inequality and Minkowski’s inequalities together with (\[usup\]), we get that $$\begin{split}
\E \vert I_2 \vert^k{\leqslant}c_2 e^{c_1 k^{(2\alpha-\beta)/\alpha-\beta}t}\left( \int_0^t\int_{{\mathbf{R}}^d\times {\mathbf{R}}^d}p_{t-s}(x-y)p_{t-s}(x-z)
f_n(y-z)\, {{\rm d}}y \, {{\rm d}}z \, {{\rm d}}s\right)^{k/2}.
\end{split}$$ Appealing to (\[rty\]), we conclude that $$\sup_{t \in [0,T]} \sup_{x \in {\mathbf{R}}^d}\E \vert I_2 \vert^k{\leqslant}c_2(T)\frac{e^{c_1 k^{(2\alpha-\beta)/\alpha-\beta}t}}{n^{\gamma k/2}}.$$

We next treat $I_1$. We look at $U^{(n,n)}_t(x)-U_t^{(n, n-1)}(x)$. Using Burkholder-Davis-Gundy and Minkowski’s inequalities together with Lemma \[heat-corr\](c), we obtain $$\begin{split}
\E|U_t^{(n,n)}(x)-U_t^{(n,n-1)}(x)|^k{\leqslant}c(k) \sup_{s \in [0,t]}\sup_{x \in {\mathbf{R}}^d} \E|U_s^{(n,n-1)}(x)-U_s^{(n,n-2)}(x)|^k \left( \int_0^t  s^{-\beta/\alpha}{{\rm d}}s\right)^{k/2}.
\end{split}$$ Iterating $n$ times this procedure and choosing $T {\leqslant}1/2$, we get $$\begin{split}
\sup_{t \in [0,T]} \E|U_t^{(n,n)}(x)-U_t^{(n,n-1)}(x)|^k&{\leqslant}c(k) T^{nk/2} {\leqslant}c(k) \left(\frac{1}{2}\right)^{nk/2} \\
&{\leqslant}c(k) \frac{1}{n^{\gamma k/2}}.
\end{split}$$ Splitting the interval $[0,T]$ into subintervals of length $\frac12$, we deduce that for all $T>0$, $$\begin{split}
\sup_{t \in [0,T]} \E|U_t^{(n,n)}(x)-U_t^{(n,n-1)}(x)|^k{\leqslant}c(k,T) \frac{1}{n^{\gamma k/2}}.
\end{split}$$ We next set $$\mathcal{D}^n_t:=\sup_{ x\in {\mathbf{R}}^d}\E|V^{(n)}_t(x)-U_t^{(n,n)}(x)|^k.$$ Using Burkholder-Davis-Gundy inequality and Minkowski’s inequalities, together with Lemma \[heat-corr\](c), and adding and substracting the term $U_s^{(n,n)}(y)$, we obtain $$\E \vert I_1 \vert^k{\leqslant}c(k,T) \int_0^t\frac{\mathcal{D}^n_s+n^{-\gamma k/2}}{(t-s)^{\beta/\alpha}}\,{{\rm d}}s$$ Using Proposition \[renew\], $$\E \vert I_1 \vert^k{\leqslant}c_2 \left(\int_0^t\frac{\mathcal{D}^n_s}{(t-s)^{\beta/\alpha}}\,{{\rm d}}s+
\frac{e^{c_1 k^{(2\alpha-\beta)/\alpha-\beta}t}}{n^{\gamma k/2}}\right).$$

Combining the bound for $I_2$ and $I_1$, we obtain $$\begin{aligned}
\mathcal{D}^n_t{\leqslant}c_2 \left( \frac{e^{c_1 k^{(2\alpha-\beta)/\alpha-\beta}t}}{n^{\gamma k/2}}+\int_0^t\frac{\mathcal{D}^n_s}{(t-s)^{\beta/\alpha}}\,{{\rm d}}s\right).\end{aligned}$$ By an appropriate use of Proposition \[renew\], we conclude that $$\label{d1}
\mathcal{D}^n_t{\leqslant}c_2 \frac{e^{c_1k^{(2\alpha-\beta)/\alpha-\beta}t}}{n^{\gamma k/2}}.$$ We now look at $u_t(x)-V_t^{(n)}(x)$ to obtain $$\begin{aligned}
u_t(x)-V_t^{(n)}(x)&=\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)[\sigma(u_s(y)-V_s^{(n)}(y))]F({{\rm d}}s\,{{\rm d}}y)\\
&+\int_0^t\int_{B(x,\,(nt)^{1/\alpha})^c}p_{t-s}(x-y)\sigma(u_s(y))F({{\rm d}}s\,{{\rm d}}y),\end{aligned}$$ which gives us $$\begin{aligned}
\E|u_t(x)-V_s^{(n)}(x)|^k&{\leqslant}c\bigg( \E\left|\int_0^t\int_{B(x,\,(nt)^{1/\alpha})}p_{t-s}(x-y)[\sigma(u_s(y)-V_s^{(n)}(y))]F({{\rm d}}s\,{{\rm d}}y)\right|^k\\
&+\E\left|\int_0^t\int_{B(x,\,(nt)^{1/\alpha})^c}p_{t-s}(x-y)\sigma(u_s(y))F({{\rm d}}s\,{{\rm d}}y)\right|^k \bigg)\\
&:=I_1+I_2.\end{aligned}$$ We bound the second term first. Using the bound on the moments of the solution together with Lemma \[heat-corr\](a), we obtain $$\begin{aligned}
I_2&{\leqslant}c_2 e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}\left[ \int_0^t\int_{B(x,\,(nt)^{1/\alpha})^{2,c}}p_{t-s}(x-y)p_{t-s}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w\, {{\rm d}}s\right]^{k/2}\\
&{\leqslant}c_2 \frac{1}{n^{(2+\beta/\alpha)k/2}}e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}.\end{aligned}$$ We now consider the first term. $$\begin{aligned}
I_1&{\leqslant}c \int_0^t\sup_{y\in {\mathbf{R}}^d}\E|u_s(y)-V_s^{(n)}(y)|^k\int_{{\mathbf{R}}^d\times{\mathbf{R}}^d}p_{t-s}(x-y)p_{t-s}(x-w)f(y-w)\,{{\rm d}}y\,{{\rm d}}w\, {{\rm d}}s\\
&{\leqslant}c \int_0^t \sup_{y\in {\mathbf{R}}^d}\E|u_s(y)-V_s^{(n)}(y)|^k\frac{1}{(t-s)^{\beta/\alpha}}\,{{\rm d}}s.\end{aligned}$$ Putting these two bounds together and using Proposition \[heat-corr\], we obtain $$\label{d2}
\sup_{s \in [0,T]} \sup_{x \in {\mathbf{R}}^d}\E|u_s(x)-V_s^{(n)}(x)|^k{\leqslant}c_2 \frac{1}{n^{(2+\beta/\alpha)k/2}}e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}.$$ Combining the estimates (\[d1\]) and (\[d2\]) and using the fact that $\gamma<2+\beta/\alpha$ we obtain the required result.

Proof of the spatial asymptotic results
=======================================

In this section we give the proofs of Theorems [\[asym\]]{.nodecor} and [\[trichotomy\]]{.nodecor}. We start with several preliminary results.

Tail estimates I
----------------

This subsection is devoted to the proof of two tail estimates which are a consequence of the sharp moment estimates in Theorem \[momentbounds\].

\[uppertail\] There exists a constant $c_{A,\alpha,\beta}>0$ such that for all $\lambda>0$ and $t>0$, $$\begin{aligned}
\sup _{x \in {\mathbf{R}}^d} {\mathrm{P}}(u_t(x)>\lambda){\leqslant}\exp\left(-\frac{c_{A,\alpha,\beta} \nu^{\beta/\alpha}}{t^{(\alpha-\beta)/\alpha}}\left|\log \frac{\lambda}{A\bar{u}_0}\right|^{(2\alpha-\beta)/\alpha}\right),\end{aligned}$$ where $A$ and $\bar{u}_0$ are defined in Theorem [\[momentbounds\]]{.nodecor}.

We start by using Chebyshev’s inequality to obtain, $$\begin{aligned}
{\mathrm{P}}(u_t(x)>\lambda)&{\leqslant}\frac{1}{\lambda^k}\E|u_t(x)|^k\\
&{\leqslant}A^k\bar{u}_0^k\lambda^{-k}e^{Ak^{(2\alpha-\beta)/(\alpha-\beta)} \nu^{-\beta/(\alpha-\beta)}t}\\
&{\leqslant}\exp\left(Ak^{\frac{2\alpha-\beta}{\alpha-\beta}} \nu^{-\beta/(\alpha-\beta)}t-k\log \frac{\lambda}{A\bar{u}_0}\right).\end{aligned}$$ The function $F(k):=Ak^{\frac{2\alpha-\beta}{\alpha-\beta}} \nu^{-\beta/(\alpha-\beta)}t-k\log \frac{\lambda}{A\bar{u}_0}$ is optimised at the point $$k^\ast=\left[\frac{ \nu^{\beta/(\alpha-\beta)}}{At}\left(\frac{\alpha-\beta}{2\alpha-\beta}\right)\left(\log \frac{\lambda}{A\bar{u}_0}\right) \right]^{(\alpha-\beta)/\alpha}.$$ Some computations then give $${\mathrm{P}}(u_t(x)>\lambda){\leqslant}\exp\left(-\frac{c_{A,\alpha,\beta} \nu^{\beta/\alpha}}{t^{(\alpha-\beta)/\alpha}}\left|\log \frac{\lambda}{A\bar{u}_0}\right|^{(2\alpha-\beta)/\alpha}\right),$$ where $c_{A,\alpha,\beta}=\frac{\alpha}{2\alpha-\beta}\left[\frac{1}{A}\left(\frac{\alpha-\beta}{2\alpha-\beta}\right) \right]^{(\alpha-\beta)/\alpha}.$

\[lowertail\] Fix $t>0$. Set $\lambda:=\frac{\underline{u}_0}{2A} e^{tk^{\alpha/(\alpha-\beta)} \nu^{-\beta/(\alpha-\beta)}/A}$ for $k{\geqslant}2$. Then there exists a constant $\tilde{c}_{A,\alpha,\beta}>0$, $$\begin{aligned}
\inf _{x \in {\mathbf{R}}^d} {\mathrm{P}}(u_t(x)>\lambda){\geqslant}\frac14 \exp\left(-\frac{\tilde{c}_{A,\alpha,\beta} \nu^{\beta/\alpha}}{t^{(\alpha-\beta)/\alpha}} \left(\log \frac{2\lambda A}{\underline{u}_0}\right)^{(2\alpha-\beta)/\alpha}\left(1+\frac{1}{\log \frac{2\lambda A}{\underline{u}_0}} \right)\right),\end{aligned}$$ where the quantities $A$ and $\underline{u}_0$ are defined in Theorem [\[momentbounds\]]{.nodecor}.

By Paley-Zygmund inequality, we have for all $k {\geqslant}2$, $$\begin{aligned}
{\mathrm{P}}(u_t(x){\geqslant}\frac{1}{2}\|u_t(x)\|_{L^{2k}(\Omega)}){\geqslant}\frac{(\E|u_t(x)|^{2k})^2}{4\E|u_t(x)|^{4k}}.\end{aligned}$$ Set $\lambda:=\frac{\underline{u}_0}{2A}e^{tk^{\alpha/(\alpha-\beta)} \nu^{-\beta/(\alpha-\beta)}/A}$. Taking into account the bounds on the moments, we obtain $$\begin{aligned}
 \label{aux1}
{\mathrm{P}}(u_t(x){\geqslant}\lambda){\geqslant}\frac{1}{4}\exp\left(-c_{A, \alpha,\beta}k^{(2\alpha-\beta)/(\alpha-\beta)} \nu^{-\beta/(\alpha-\beta)}t+k\log \left(\frac{\tilde{u}_0^4}{A^8}\right)\right),\end{aligned}$$ where $c_{A, \alpha,\beta}:=2^{(2\alpha-\beta)/(\alpha-\beta)}[A2^{(2\alpha-\beta)/(\alpha-\beta)}-\frac{2}{A}]$ and $\tilde{u}_0=\frac{\underline{u}_0}{\bar{u}_0}.$ Finally, some computations we get the desired bound.

Insensitivity analysis
----------------------

The next theorem is crucial in the proof the of spatial asymptotic result when the initial condition is not bounded below. Intuitively, we study how the solution is sensible to changes to the initial data, and we conclude that when $R$ is large, the values of the solution in a given ball of radius $R$ are insensitive to the changes of the initial value outside the ball.

\[sensitivity\] Let $a\in {\mathbf{R}}^d$ and $R>1$. Let $u$ and $v$ be the solution to with respective initial conditions $u_0$ and $v_0$. Suppose that on $B(a,\,2R)$, $u_0(x)=v_0(x)$. Then for all $k{\geqslant}2$ there exist positive constants $c_1,c_2$ such that for all $t>0$, $$\sup_{x\in B(a,\,R)}\E|u_t(x)-v_t(x)|^k{\leqslant}c_2 \|u_0-v_0\|_{L^\infty({\mathbf{R}}^d)}\frac{1}{R^{\alpha k/2}}e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}.$$

From the mild solution, we have $$\begin{aligned}
u_t(x)-v_t(x)=(\cG u_0)_t(x)-(\cG v_0)_t(x) +\int_0^t\int_{{\mathbf{R}}^d}p_{t-s}(x-y)[\sigma(u_s(y))-\sigma(v_s(y))]\,F({{\rm d}}y\,{{\rm d}}s),\end{aligned}$$ where $(\cG u_0)_t(x):=(p_t\ast u_0)(x)$. We obtain $$\begin{aligned}
\|u_t(x)-v_t(x)\|^{2}_{L^k(\Omega)}
&{\leqslant}c \big(\|(\cG u_0)_t(x)-(\cG v_0)_t(x)\|^{2}_{L^k(\Omega)}\\
&+\|\int_0^t\int_{{\mathbf{R}}^d}p_{t-s}(x-y)[\sigma(u_s(y))-\sigma(v_s(y))]\,F({{\rm d}}y\,{{\rm d}}s)\|^{2}_{L^k(\Omega)}\big)\\
&:=I_1+I_2.\end{aligned}$$ We bound $I_1$ first. Noting that $x\in B(a,\,R)$ and $y\in B(a,\,2R)^c$, we have $$\begin{aligned}
I_1&{\leqslant}\|u_0-v_0\|^2_{L^\infty({\mathbf{R}}^d)}\left[\int_{B(a,\,2R)^c}p_t(x-y)\,{{\rm d}}y\right]^2\\
&{\leqslant}c_1 \|u_0-v_0\|^2_{L^\infty({\mathbf{R}}^d)}\frac{t^2}{R^{2\alpha}}\\
&{\leqslant}c_1 \|u_0-v_0\|^2_{L^\infty({\mathbf{R}}^d)}\frac{e^{c_2 k^{(2\alpha-\beta)/(\alpha-\beta)}t}}{R^{2\alpha}}.\end{aligned}$$ We use Burkholder-Davis-Gundy inequality and Minkowski’s inequalities to bound the second term as follows $$\begin{aligned}
I_2&{\leqslant}c k  \left\{\E\left[\int_0^t\int_{{\mathbf{R}}^d\times {\mathbf{R}}^d}p_{t-s}(x-y_1)p_{t-s}(x-y_1)[u_s(y_1)-v_s(y_1)][u_s(y_2)-v_s(y_2)]f(y_1,\,y_2){{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\right]^{k/2}\right\}^{2/k}\\
&{\leqslant}k \int_0^t\int_{{\mathbf{R}}^d\times {\mathbf{R}}^d}p_{t-s}(x-y_1)p_{t-s}(x-y_1)\|[u_s(y_1)-v_s(y_1)][u_s(y_2)-v_s(y_2)]\|_{L^{k/2}(\Omega)}f(y_1,\,y_2){{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s.\end{aligned}$$ We split the integral on the right hand side as follows $$\begin{aligned}
 & \E\left[\int_0^t\int_{B(x,\,R)\times B(x,\,R)}\left[\cdots\right]{{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\right]+ \E\left[\int_0^t\int_{B(x,\,R)^c\times B(x,\,R)}\left[\cdots\right]{{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\right]\\
&+\E\left[\int_0^t\int_{B(x,\,R)\times B(x,\,R)^c}\left[\cdots\right]{{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\right]+\E\left[\int_0^t\int_{B(x,\,R)^c\times B(x,\,R)^c}\left[\cdots\right]{{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\right]\\
&=I_3+I_4+I_5+I_6.\end{aligned}$$ Since $x\in B(a,\,R)$, we can bound $I_3$ as follows, $$\begin{aligned}
I_3&{\leqslant}\int_0^t\sup_{y\in B(x,\,R)}\|u_s(y)-v_s(y)\|_{L^k(\Omega)}^{2}\int_{B(x,\,R)\times B(x,\,R)}p_{t-s}(x-y_1)p_{t-s}(x-y_1)f(y_1,\,y_2){{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\\
&{\leqslant}c \int_0^t\sup_{y\in B(a,\,2R)}\|u_s(y)-v_s(y)\|_{L^k(\Omega)}^{2}\frac{1}{(t-s)^{\beta/\alpha}}\,{{\rm d}}s.\end{aligned}$$ We now use Lemma \[heat-corr\](b) to obtain $$\begin{aligned}
I_4&{\leqslant}\int_0^t\sup_{y\in {\mathbf{R}}^d}\|u_s(y)-v_s(y)\|_{L^k(\Omega)}^{2}\int_{B(x,\,R)^c\times B(x,\,R)}p_{t-s}(x-y_1)p_{t-s}(x-y_1)f(y_1,\,y_2){{\rm d}}y_1\, {{\rm d}}y_2\,{{\rm d}}s\\
&{\leqslant}c_1 \frac{e^{c_2k^{(2\alpha-\beta)/(\alpha-\beta)}t}}{R^\alpha}. \end{aligned}$$ Similarly, we can use Lemma \[heat-corr\](b) again to $$I_5{\leqslant}c_1 \frac{e^{c_2k^{(2\alpha-\beta)/(\alpha-\beta)}t}}{R^\alpha}\quad\text{and}\quad I_6{\leqslant}c_1\frac{e^{c_2k^{(2\alpha-\beta)/(\alpha-\beta)}t}}{R^\alpha}.$$ Combining those bounds, we obtain $$\begin{aligned}
&\sup_{y\in B(a,\,R)}\|u_t(y)-v_t(y)\|_{L^k(\Omega)}^{2}\\
&{\leqslant}c_1 \|u_0-v_0\|^2_{L^\infty({\mathbf{R}}^d)}\frac{e^{c_2k^{(2\alpha-\beta)/(\alpha-\beta)}t}}{R^\alpha}+  k c_3\int_0^t\sup_{y\in B(a,\,2R)}\|u_s(y)-v_s(y)\|_{L^k(\Omega)}^{2}\frac{1}{(t-s)^{\beta/\alpha}}\,{{\rm d}}s.\end{aligned}$$ We set $$f_R(t):=\sup_{y\in B(a,\,R)}\|u_t(y)-v_t(y)\|_{L^k(\Omega)}^{2}$$ and $$A_R(t):=\|u_0-v_0\|^2_{L^\infty({\mathbf{R}}^d)}\frac{e^{c_2k^{(2\alpha-\beta)/(\alpha-\beta)}t}}{R^\alpha}.$$

We now use Proposition \[renew\] to arrive at the result.

Tail Estimates II
-----------------

In this subsection we are going to prove tail estimates when the initial condition is not bounded below. The next result is an extension of [@CDKh Theorem 2.4]

\[tail:asymp\] Suppose that $u$ and $u_0$ are as in Theorem [\[trichotomy\]]{.nodecor}. Then, there exist positive constants $K_1,K_2$ such that for all $\lambda>0$, $$\begin{aligned}
-K_1\frac{\Lambda^{(2\alpha-\beta)/\alpha}}{t^{(\alpha-\beta)/\alpha}}{\leqslant}\liminf_{|x|\rightarrow \infty}\frac{\log {\mathrm{P}}(u_t(x)>\lambda)}{\log |x|}{\leqslant}\limsup_{|x|\rightarrow \infty}\frac{\log {\mathrm{P}}(u_t(x)>\lambda)}{\log |x|}{\leqslant}-K_2\frac{\Lambda^{(2\alpha-\beta)/\alpha}}{t^{(\alpha-\beta)/\alpha}},\end{aligned}$$ uniformly for all $t$ in every fixed compact subset of $(0, \infty)$.

We prove the lower bound first. Fix $a \in {\mathbf{R}}^d$. Let $w_t$ be the solution to when the initial condition is given by the following $$\begin{aligned}
w_0(x):=u_0(|x|\vee(3|a|))\quad \text{for all}\quad x\in {\mathbf{R}}^d.\end{aligned}$$ Since $w_0{\leqslant}u_0$, the weak comparison principle Theorem \[weakc\] tells us that for all $t>0$ and $x\in {\mathbf{R}}^d$, $$\begin{aligned}
w_t(x){\leqslant}u_t(x). \end{aligned}$$ This means that finding a lower bound on the tail distribution of $u_t(x)$ amounts to finding a lower bound for the corresponding distribution of $w_t(x)$.

Now, let $u_t^a(x)$ be the solution to when initial condition $u_0(3 \vert a \vert)$. Fix $\lambda>0$. Then, by Theorem \[sensitivity\], whenever $R=|a|>1$, $$\begin{aligned}
\sup_{x\in B(a,\,|a|)}{\mathrm{P}}(|u_t^a(x)-w_t(x)|>\lambda)&{\leqslant}\sup_{x\in B(a,\,|a|)}\frac{\E|u_t^a(x)-w_t(x)|^k}{\lambda^k}\\
&{\leqslant}c_2 \frac{1}{\lambda^k \vert a \vert^{\alpha k/2}}e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}.\end{aligned}$$

Recall that $u_0(3\vert a\vert)$ is decreasing in $a$ and $$\lim_{a\rightarrow \infty}u_0(3\vert a \vert)=0.$$ Therefore, we can then take $\vert a \vert$ large enough so that $$k:=\left(\frac{A\nu^{\frac{\beta}{\alpha-\beta}}}{t}\log \left(\frac{4A \lambda}{u_0(3\vert a \vert)}\right)\right)^{\frac{\alpha-\beta}{\alpha}} {\geqslant}2,$$ and $$\frac{1}{2}{\leqslant}\left|\log \frac{4\lambda A}{u_0(3\vert a \vert)}\right|.$$ We now use Lemma \[lowertail\] to obtain $$\begin{aligned}
\inf_{x\in B(a,\,|a|)}{\mathrm{P}}(u^a_t(x){\geqslant}2\lambda){\geqslant}\frac14 \exp\left(-\frac{\tilde{c}_{A,\alpha,\beta} \nu^{\beta/\alpha}}{t^{(\alpha-\beta)/\alpha}} \left(\log \frac{2\lambda A}{\underline{u}_0}\right)^{(2\alpha-\beta)/\alpha}\right).\end{aligned}$$ Upon taking $\vert a \vert$ larger if required so that $$\left|\log \frac{4\lambda A}{u_0(3\vert a \vert)}\right|{\leqslant}2|\log u_0(3\vert a \vert)|,$$ we can use the above together with the definition of $\Lambda$ to write $$\begin{split}
\inf_{x\in B(a,\,|a|)}{\mathrm{P}}(u_t(x){\geqslant}\lambda)&{\geqslant}\inf_{x\in B(a,\,|a|)}{\mathrm{P}}(u^a_t(x){\geqslant}2\lambda)-\sup_{x\in B(a,\,|a|)}{\mathrm{P}}(|u_t^a(x)-w_t(x)|>\lambda)\\
&{\geqslant}\frac14  \exp\left(-\frac{c_{A,\alpha,\beta} \nu^{\beta/\alpha}\Lambda^{(2\alpha-\beta)/\alpha}}{t^{(\alpha-\beta)/\alpha}}\log |a|\right) - c_2 \frac{1}{\lambda^k \vert a \vert^{\alpha k/2}}e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}.
\end{split}$$ For any fixed $t>0$, we choose $k$ large enough so that for $a$ large, we obtain $$\begin{aligned}
\inf_{x\in B(a,\,|a|)}{\mathrm{P}}(u_t(x){\geqslant}\lambda)&{\geqslant}\frac{1}{8|a|^{K_1\frac{\Lambda^{(2\alpha-\beta)/\alpha}}{t^{(\alpha-\beta)/\alpha}}}}.\end{aligned}$$ The above immediately gives the lower bound needed. We now turn our attention to the upper bound. The proof uses a similar strategy as the above. We look at $w_t$, the solution to but this time, the initial condition is defined by $$\begin{aligned}
w_0(x):=u_0(|x|\wedge 2|a|),\end{aligned}$$ so that now we have $w_0(x){\geqslant}u_0(x)$ which gives us $w_t(x){\geqslant}u_t(x)$ by the weak comparison principle. Now consider $u^a_t$ a solution with constant initial condition given by $$\begin{aligned}
z_0(x):=u_0(2|a|).\end{aligned}$$ We choose $a$ large enough such that $$\left|\log \frac{\lambda}{2Au_0(2|a|)}\right|{\geqslant}\frac{|\log u_0(2|a|)|}{2}.$$ Then, by Theorem \[sensitivity\] and Lemma \[uppertail\], for $\vert a \vert$ large enough $$\label{tail2}\begin{split}
\sup_{x\in B(a,\,|a|)}{\mathrm{P}}(u_t(x)>\lambda)&{\leqslant}\sup_{x\in B(a,\,|a|)} {\mathrm{P}}(w_t(x)>\lambda)\\
 &{\leqslant}\sup_{x\in B(a,\,|a|)}{\mathrm{P}}(u_t^a(x)>\lambda/2)+\sup_{x\in B(a,\,|a|)}{\mathrm{P}}(|u_t^a(x)-w_t(x)|>\lambda/2)\nonumber\\
&{\leqslant}\exp\left(-\frac{c_{A,\alpha,\beta}\Lambda^{(2\alpha-\beta)/\alpha}}{
t^{(\alpha-\beta)/\alpha}}\log |a|\right)+c_2 \frac{1}{\lambda^k \vert a \vert^{\alpha k/2}}e^{c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t}.
\end{split}$$ By choosing $k$ large, we obtain upon taking $a$ large enough, $$\begin{aligned}
\sup_{x\in B(a,\,|a|)}{\mathrm{P}}(u_t(x)>\lambda){\leqslant}\frac{1}{|a|^{K_2\frac{\Lambda^{(2\alpha-\beta)/\alpha}}{t^{(\alpha-\beta)/\alpha}}}}.\end{aligned}$$ This finishes the proof.

Proof of Theorem \[asym\]
-------------------------

Let $t>0$ and set $$\begin{aligned}
L:= \frac{\underline{u}_0}{6A}\exp\left[\delta_1t^{(\alpha-\beta)/(2\alpha-\beta)}\left|\log R\right|^{\alpha/(2\alpha-\beta)} \nu^{-\beta/(2\alpha-\beta)}\right],\end{aligned}$$ where $\delta_1$ be a positive constant. Then, we choose $$k=(A \delta_1)^{(\alpha-\beta)/\alpha}\left(\frac{\vert \log(R) \vert}{t}\right)^{(\alpha-\beta)/(2\alpha-\beta)}\nu^{\beta/(2\alpha-\beta)},$$ so that $L$ becomes $$\begin{aligned}
L:= \frac{\underline{u}_0}{6A}e^{tk^{\alpha/(\alpha-\beta)} \nu^{-\beta/(\alpha-\beta)}/A}.\end{aligned}$$ We now apply inequality (\[aux1\]) to obtain for sufficiently large $R$, $$\begin{aligned}
&{\mathrm{P}}(u_t(x)>3L)\\
&{\geqslant}\frac{1}{4}\exp\left(-c_{A, \alpha,\beta}(A\delta_1)^{(2\alpha-\beta)/\alpha} \vert \log R \vert-\log \left(A^8\right)(A\delta_1)^{(\alpha-\beta)/\alpha} \left(\vert \log R \vert/t\right)^{(\alpha-\beta)/(2\alpha-\beta}\nu^{\beta/(2\alpha-\beta)}\right) \\
&{\geqslant}\frac{1}{4 R^{\delta_2}},  \end{aligned}$$ where $\delta_1$ is chosen such that $\delta_2<2$.

Let $N>0$ and choose $x_1,x_2,\ldots, x_N \in {\mathbf{R}}^d$ such that $\vert x_i-x_j\vert {\geqslant}2n^{1+1/\alpha}t^{1/\alpha}$ for $i\not=j$. Lemma \[ind\] then implies that the $U^{(n,n)}_t(x_i)$’s are independent for large enough $n$. We have $$\begin{aligned}
{\mathrm{P}}(\max_{1{\leqslant}i{\leqslant}N}u_t(x_i)<L)&{\leqslant}{\mathrm{P}}(\max_{1{\leqslant}i{\leqslant}N}|U^{(n,n)}_t(x_i)|<2L)\\
&\qquad +{\mathrm{P}}(|u_t(x_i)-U_t^{(n,n)}(x_i)|>L\quad\text{for some}\quad1{\leqslant}i{\leqslant}N)\\
&:=I_1+I_2.\end{aligned}$$ We will look at the second term first. By Lemma \[approx\], for all $k {\geqslant}2$ and large $n$, $$\begin{aligned}
I_2{\leqslant}\frac{N\E|u_t(x_i)-U_t^{(n,n)}(x_i)|^k}{L^k}{\leqslant}c_2  \frac{N}{n^{\gamma k/2}}e^{c_1k^{(2\alpha-\beta)/(\alpha-\beta)}t},\end{aligned}$$ where we have chosen $R$ large enough such that $L {\geqslant}1$.

We now choose $n{\geqslant}N^{10/(3\gamma)}$ so that we have $$\begin{aligned}
I_2{\leqslant}c_2\frac{1}{N^{2/3}}e^{c_1k^{(2\alpha-\beta)/(\alpha-\beta)}t}.\end{aligned}$$ Upon choosing $N$ to be an integer greater than $R^3$, we obtain $$\begin{aligned}
I_2 {\leqslant}c(T,k) \frac{1}{R^2}.\end{aligned}$$ To bound $I_1$, we have for large enough $R$, $$\begin{aligned}
{\mathrm{P}}(U^{(n,n)}_t(x_i){\geqslant}2L)&{\geqslant}{\mathrm{P}}(|u_t(x_i)|{\geqslant}3L)-{\mathrm{P}}(|u_t(x_i)-U_t^{(n,n)}(x_i)|{\geqslant}L)\\
&{\geqslant}c\left(\frac{1}{R^{\delta_2}} -\frac{1}{R^2}\right){\geqslant}\frac{c}{R^2}.\end{aligned}$$ By independence, we have $$\begin{aligned}
I_1 {\leqslant}\left(1-{\mathrm{P}}(U^{(n,n)}_t(x_i){\geqslant}2L)\right)^N.\end{aligned}$$ Combining the above and bearing in mind that $N$ is larger than $R^3$, we obtain $$\begin{aligned}
{\mathrm{P}}(\max_{1{\leqslant}i{\leqslant}N}u_t(x_i)<L)&{\leqslant}\left(1-{\mathrm{P}}(U^{(n,n)}_t(x_i){\geqslant}2L)\right)^N+\frac{c}{R^2}\\
&{\leqslant}\frac{c}{R^2},\end{aligned}$$ for $R$ large enough. And hence by a standard monotonicity argument, we have $$\begin{aligned}
{\mathrm{P}}\left(\sup_{x\in B(0,\,R)}u_t(x){\leqslant}\frac{\underline{u}_0}{6A}\exp\left[\delta_1t^{(\alpha-\beta)/(2\alpha-\beta)}\left(\log R\right)^{\frac{\alpha}{(2\alpha-\beta)}}\nu^{-\beta/(2\alpha-\beta)}\right]\right){\leqslant}\frac{c}{R^2}.\end{aligned}$$ We now use Borel Cantelli lemma to obtain that almost surely, for $R\rightarrow \infty$, we have $$\begin{aligned}
\sup_{x\in B(0,\,R)}u_t(x){\geqslant}\frac{\underline{u}_0}{6A}\exp\left[\delta_1t^{(\alpha-\beta)/(2\alpha-\beta)}\left(\log R\right)^{\frac{\alpha}{2\alpha-\beta}}\right],\end{aligned}$$ which concludes the proof of the lower bound. We now prove the upper bound. Set $$\begin{aligned}
U:= A\bar{u}_0\exp\left[\delta_3t^{(\alpha-\beta)/(2\alpha-\beta)}\left(\log R\right)^{\frac{\alpha}{2\alpha-\beta}}\nu^{-\beta/(2\alpha-\beta)}\right],\end{aligned}$$ for some positive constant $\delta_3$. For $x\in \bZ^d$, denote the cube of side length $1$ by $Q_x$. Let $R$ be a positive integer and decompose $[-R,\,R]^d$ into cubes of the form $Q_x$ so that $[-R,\,R]^d=\cup_{x\in S}Q_x$ where $S$ is some finite set. By Proposition \[cty\], for any $x\in {\mathbf{R}}^d$ and $k {\geqslant}2$, we have $$\begin{aligned}
\label{spatial-cty}
\E\left[\sup_{w,\,y\in Q_x}|u_t(w)-u_t(y)|^{k}\right]{\leqslant}c_2 \exp\left[c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t\right].\end{aligned}$$ We can now write $$\begin{aligned}
{\mathrm{P}}\left(\sup_{x\in [-R,\,R]^d}u_t(x){\geqslant}2 U\right)&{\leqslant}{\mathrm{P}}\left(\max_{x\in S}u_t(x){\geqslant}U\right)+{\mathrm{P}}\left(\max_{x\in S}\sup_{y\in Q_x}|u_t(y)-u_t(x)|{\geqslant}U \right)\\
&:=I_1+I_2.\end{aligned}$$ To bound $I_1$, we use Lemma \[uppertail\] to obtain $$\begin{aligned}
I_1{\leqslant}|S|{\mathrm{P}}(u_t(x){\geqslant}U){\leqslant}\frac{c}{R^{\delta_4}},\end{aligned}$$ where the constant $\delta_3$ is chosen so that $\delta_4>1$. We now bound $I_2$ by making use of , $$\begin{aligned}
I_2{\leqslant}|S|{\mathrm{P}}\left(\sup_{y\in Q_x}|u_t(y)-u_t(x)|{\geqslant}U \right){\leqslant}\frac{c_2 |S|\exp(c_1 k^{(2\alpha-\beta)/(\alpha-\beta)}t)}{\exp(k\delta_3 t^{(\alpha-\beta)/(2\alpha-\beta)}(\log R)^{\alpha/(2\alpha-\beta)}\nu^{-\beta/(2\alpha-\beta)})}.\end{aligned}$$ We now set $k=\delta_5\left(\frac{\log R}{t}\right)^{(\alpha-\beta)/(2\alpha-\beta)}\nu^{\beta/(2\alpha-\beta)}$ to obtain $$\begin{aligned}
I_2{\leqslant}\frac{c}{R^{\delta_6}},\end{aligned}$$ where we choose $\delta_5$ so that $\delta_6>1$. We can conclude that $$\begin{aligned}
\sum_{R=1}^\infty{\mathrm{P}}\left(\sup_{x\in [-R,\,R]^d}u_t(x)>2U\right)<\infty.\end{aligned}$$ We can now use Borel-Cantelli and the fact that $B(0,\,R)\subset [-R,\,R]^d$ to finish the proof.

Proof of Theorem \[trichotomy\].
--------------------------------

We split the proof into two parts. In the first part we assume that $\Lambda>0$. Consider the sequence $\{x_n\}_{n {\geqslant}1}\subset {\mathbf{R}}^d$ such that $\vert x_n\vert=n^{1/2}$ and all $x_n$ lie on a straight line through the origin. We next choose $$\lambda\in (0,\,\Lambda),$$ and consider $$\quad t(j,\,n):=\frac{jT}{n}, \quad \text{ for } \quad j\in [\frac{n\tau}{T},\,n]\cap {\mathbf{Z}}.$$ We look at the following parameters $\tau$ and $T$ such that $$0<\tau<T:=\frac{K_2^{\alpha/(\alpha-\beta)} \lambda^{(2\alpha-\beta)/(\alpha-\beta)}}{2^{\alpha/(\alpha-\beta)} },$$ where $K_2$ is the constant in the statement of Theorem \[tail:asymp\]. Then, by Theorem \[tail:asymp\], for all $\theta>0$, $t\in (\tau,\,T)$ and large enough $n$, $$\begin{aligned}
{\mathrm{P}}\left(\max_{j\in [\frac{n\tau}{T},\,n]}u_{t(j,\,n)}(x_n)>\theta\right)&{\leqslant}\sum_{j\in [\frac{n\tau}{T},\,n] \cap {\mathbf{Z}}}{\mathrm{P}}\left(u_{t(j,\,n)}(x_n)>\theta\right)\\
&{\leqslant}c \sum_{j\in [\frac{n\tau}{T},\,n] \cap {\mathbf{Z}}} \exp \left( -\frac{K_2 \lambda^{(2\alpha-\beta)/\alpha} }{t(j,n)^{(\alpha-\beta)/\alpha}}\log \vert x_n \vert\right)\\
&{\leqslant}c n \exp \left( -\frac{K_2 \lambda^{(2\alpha-\beta)/\alpha} }{T^{(\alpha-\beta)/\alpha}}\log n^{1/2}\right)
\\
&{\leqslant}\frac{c}{n^3}.\end{aligned}$$ An application of Borel-Cantelli lemma gives us $$\begin{aligned}
\lim_{n\rightarrow \infty} \max_{j\in [\frac{n\tau}{T},\,n]\cap {\mathbf{Z}}}u_{t(j,\,n)}(x_n)=0\quad \text{a.s}.\end{aligned}$$ We now use Proposition \[cty\] to obtain for all $\theta>0$, $$\begin{aligned}
{\mathrm{P}}\Big\{\sup_{t\in (\tau,\,T)}&\min_{j\in [\frac{n\tau}{T},\,n]\cap {\mathbf{Z}}} |u_{t(j,\,n)}(x_n)-u_t(x_n)|>\theta\Big\}\\
&{\leqslant}{\mathrm{P}}\Big\{\sup_{t\in (\tau,\,T): |s-t|{\leqslant}T/n}|u_s(x_n)-u_t(x_n)|>\theta\Big\}\\
&{\leqslant}\frac{c_{T,k}}{n^{\tilde{\eta} k}}.\end{aligned}$$ By choosing $k$ large enough, we can apply Borel-Cantelli and use the above to see that $$\begin{aligned}
 \label{aux5}
\lim_{n\rightarrow \infty}\sup_{t\in (\tau,\,T)}u_t(x_n)=0\quad \text{a.s}.\end{aligned}$$ We next use Proposition \[cty\], to get for all $\theta>0$, $$\begin{aligned}
{\mathrm{P}}\Big\{\sup_{t\in (\tau,\,T)}\sup_{x\in [x_n,\,x_{n+1}]} &|u_t(x_n)-u_t(x)|{\geqslant}\theta\Big\}\\
&{\leqslant}c n^{4/\alpha} {\mathrm{P}}\Big\{\sup_{t\in (\tau,\,T)}\sup_{\vert x-y\vert {\leqslant}\frac{1}{n}} |u_t(x)-u_t(y)|{\geqslant}\theta\Big\}\\
&{\leqslant}c\frac{n^{1/2}}{n^{\eta k}}.\end{aligned}$$ We then take $k$ large enough, use Borel-Cantelli again and (\[aux5\]) to conclude that $$\begin{aligned}
\lim_{\vert x\vert \rightarrow \infty}\sup_{t\in (\tau,\,T)}u_t(x)=0\quad \text{a.s}, \end{aligned}$$ where in the above, $x$ tends to infinity along a fixed straight line. Since the line is arbitrary and $u$ is almost surely continuous (Proposition \[cty\]), it follows that $${\mathrm{P}}\left( \sup_{x \in {\mathbf{R}}^d} u_t(x)<\infty, \text{ for all } t \in (\tau, T)\right)=1.$$ The above is valid when $\Lambda=\infty$ in which case we can take $T$ as large as we want. Taking $\tau$ close to zero finishes the second part of the proof. We next assume that $\Lambda<\infty$. Fix $\theta>0$ and set $$\begin{aligned}
E_t(x):=\{\omega\in \Omega: u_t(x){\leqslant}\theta \}\quad \text{for every}\quad t>0,  x \in {\mathbf{R}}^d.\end{aligned}$$ We will show that solution is almost surely unbounded for large enough times. Let $$\tau> (\Lambda^{(2\alpha-\beta)/\alpha} 2K_1/\alpha)^{\alpha/(\alpha-\beta)}\quad \text{and} \quad T>\tau.$$ According to Theorem \[tail:asymp\], for every $\lambda \in (\Lambda, (\alpha\tau^{(\alpha-\beta)/\alpha}/(2K_1) )^{\alpha/(2\alpha-\beta)}]$, we can find a real number $n(\lambda, \theta)>1$ such that $$\begin{aligned}
 \label{k1}
{\mathrm{P}}(E_t(x)){\leqslant}\left( 1-|x|^{-K_1 \lambda^{(2\alpha-\beta)/\alpha}/t^{(\alpha-\beta)/\alpha}}\right){\leqslant}\left( 1-\frac{1}{|x|^{\alpha/2}}\right),\end{aligned}$$ uniformly for all $\vert x \vert {\geqslant}n(\lambda, \theta)$ and $t\in (\tau,\,T)$. Consider the events $$\begin{aligned}
E^{(n)}_t(x):=\{\omega\in \Omega: U^{(n,\,n)}_t(x){\leqslant}2\theta \}\quad \text{for every}\quad x \in {\mathbf{R}}^d, n {\geqslant}1.\end{aligned}$$ By Lemma \[approx\], we get $$\label{k2} \begin{split}
\sup_{t\in(\tau,\,T)}{\mathrm{P}}(E_t(x)\backslash E_t^{n}(x))&{\leqslant}\sup_{t\in(\tau,\,T)}{\mathrm{P}}\left(\left|u_t(x)-U^{(n,\,n)}_t(x) \right|{\geqslant}\theta\right)\\
&{\leqslant}\frac{c_{T,k}}{n^{\gamma k/2}}.
\end{split}$$ Therefore, $$\begin{aligned}
{\mathrm{P}}\left(\bigcap_{x \in [n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d} E_t(x)\right)&{\leqslant}{\mathrm{P}}\left(\bigcap_{\ell\in [n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d\cap {\mathbf{Z}}^d}E_t(\ell)\right)\\
&{\leqslant}{\mathrm{P}}\left(\bigcap_{\ell\in [n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d\cap {\mathbf{Z}}^d}E^{(n)}_t(\ell)\right)+\frac{c}{n^{\gamma k/2}},\end{aligned}$$ uniformly for all $n {\geqslant}1$ and $t \in (\tau, T)$, where $a>1$ is constant. We will now look at the first term of the above display. Set $x_1:=(n^{a+\frac{1}{\alpha}},\ldots, n^{a+\frac{1}{\alpha}})\in {\mathbf{R}}^d$ and define iteratively for $j{\geqslant}1$, $$\begin{aligned}
x_{j+1}:=x_j+(2n^{1+\frac{1}{\alpha}}t^{1/\alpha},\ldots,2n^{1+\frac{1}{\alpha}}t^{1/\alpha}).\end{aligned}$$ Let $$\begin{aligned}
\gamma_n:=\max\left\{j{\geqslant}1: x_{j,i}{\leqslant}2n^{a+\frac{1}{\alpha}}, \text{ for all } i=1,\ldots, d\right\},\end{aligned}$$ where $x_j=(x_{j,1}, \ldots, x_{j,d})$. Observe that $$\begin{aligned}
\gamma_n{\geqslant}\frac{n^{a-1}}{2T^{1/\alpha}}.\end{aligned}$$ By independence (Lemma \[ind\]), (\[k2\]) and (\[k1\]), we get $$\begin{aligned}
{\mathrm{P}}\left(\bigcap_{\ell\in [n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d\cap {\mathbf{Z}}^d}E^{(n)}_t(\ell)\right)&{\leqslant}{\mathrm{P}}\left(\bigcap_{j=1}^{\gamma_n}E^{(n)}_t(x_j)\right)=\prod_{j=1}^{\gamma_n}{\mathrm{P}}\left(E^{(n)}_t(x_j)\right)\\
&{\leqslant}\prod_{j=1}^{\gamma_n}[{\mathrm{P}}\left(E_t(x_j)\right)+\frac{c}{n^{\gamma k/2}}]\\
&{\leqslant}\left[1-\frac{1}{n^{\frac{a\alpha}{2}+\frac{1}{2}}}+\frac{c}{n^{\gamma k/2}}\right]^{\gamma_n}.\end{aligned}$$ We now take $k$ larger if necessary to obtain $$\begin{aligned}
{\mathrm{P}}\left(\bigcap_{\ell\in [n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d\cap {\mathbf{Z}}^d}E^{(n)}_t(\ell)\right)&{\leqslant}\exp(-c_1n^{-\frac{\alpha a}{2}-\frac{1}{2}+a-1}).\end{aligned}$$ We choose $a>\frac{5}{(2-\alpha)}$, so that we have $$\begin{aligned}
{\mathrm{P}}\left(\bigcap_{\ell\in [n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d\cap {\mathbf{Z}}^d}E^{(n)}_t(\ell)\right)&{\leqslant}\exp(-c_1n).\end{aligned}$$ Combining the above estimates, we have for large enough $k$, $$\begin{aligned}
 \label{k3}
\sup_{t\in (\tau,\,T)}{\mathrm{P}}\left(\sup_{x \in [n^{a+\frac{1}{\alpha}}, 2n^{a+\frac{1}{\alpha}}]^d}u_t(x){\leqslant}\theta\right){\leqslant}\frac{c}{n^{\gamma k/2}}.\end{aligned}$$ We next write $$\begin{aligned}
{\mathrm{P}}&\left(\inf_{t\in (\tau,\,T)}\sup_{x\in[n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d}u_t(x){\leqslant}\theta \right){\leqslant}{\mathrm{P}}\left(\inf_{1{\leqslant}i{\leqslant}n}\sup_{x\in[n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d}u_{t_i}(x){\leqslant}2\theta \right)\\
&\qquad \qquad \qquad \qquad +{\mathrm{P}}\left(\inf_{|t-s|<1/n}\sup_{x\in[n^{a+\frac{1}{\alpha}},\,2n^{a+\frac{1}{\alpha}}]^d}\left|u_s(x)-u_t(x)\right|{\geqslant}\theta \right):=I_1+I_2.\end{aligned}$$ From (\[k3\]), we can bound the first term as follows $$\begin{aligned}
I_1{\leqslant}\frac{c n}{n^{\gamma k/2}}.\end{aligned}$$ We now look at the second term. Using Proposition \[cty\] we obtain that $$\begin{aligned}
I_2 {\leqslant}\sum_{k=1}^{1+n^{a+1/\alpha}}{\mathrm{P}}\left(\inf_{|t-s|<1/n}\sup_{x\in (k,k+1)^d}\left|u_s(x)-u_t(x)\right|{\geqslant}\theta \right){\leqslant}\frac{c}{n^\kappa},\end{aligned}$$ where $\kappa$ can be made as large as possible. Combining the above estimates, we conclude that $$\begin{aligned}
{\mathrm{P}}\left(\inf_{t\in (\tau,\,T)}\sup_{x\in {\mathbf{R}}^d} u_t(x)<\theta \right)=0.\end{aligned}$$ If $\Lambda=0$, then we can choose $\tau$ as close to zero as we want and hence the final part of the theorem is proved. For each $N{\geqslant}1$, set $$\begin{aligned}
\mathcal{T}_N:=\inf\{t>0: \sup_{x\in {\mathbf{R}}^d} u_t(x){\geqslant}N \} \end{aligned}$$ and let $\mathcal{T}:=\lim_{N\rightarrow \infty} \mathcal{T}_N$. The first part of the theorem is then proved by using exactly the same argument as in the main theorem of [@CDKh]. We leave it to the reader to fill in the details.

Proof of the comparison principle and strict positivity
=======================================================

In order to prove the strong comparison principle (Theorem \[strongcomparison\]), we need the next two preliminary results which are extensions of [@ChenHuang Lemmas 7.1 and 7.2] (see also [@CK Lemmas 4.1 and 4.3]). In particular, the proof of the next proposition is new compared with that of [@ChenHuang Lemma 7.1] or [@CK Lemma 4.1].

\[initial\] Let $M>0$. For all $R>0$ and $t>0$, there exist constants $0<c_R<1$ and $1<m_0(t,R)<\infty$ such that for all $m {\geqslant}m_0$, $$\begin{aligned}
\int_{B(0,\,R)}p_{s}(x-y)\,{{\rm d}}y{\geqslant}c_R \quad \text{for all }\quad (s,x) \in A_{m,t,R},\end{aligned}$$ where $$A_{m,t,R}:=\{(s,x): x\in B(0,\,R+M(t/m)^{1/\alpha}) \quad\text{and} \quad  \frac{t}{2m}{\leqslant}s{\leqslant}\frac{t}{m}\}.$$

We take $m$ large enough so that $(\frac{t}{m})^{1/\alpha} {\leqslant}R$. Then, using the the lower bound (\[ee2\]), we obtain $$\begin{aligned}
\int_{B(0,\,R)}p_s(x-y)\,{{\rm d}}y
&{\geqslant}\int_{B(0,\,R)\cap B(x,\,2M(t/m)^{1/\alpha})}p_s(x-y)\,{{\rm d}}y\\
&{\geqslant}c \left( \frac{t}{m}\right)^{d/\alpha} s^{-d/\alpha} \\
&{\geqslant}c \end{aligned}$$ where the constant $c$ might depend on $R$ and $M$ but can be chosen to be strictly less than $1$.

\[stoch\] Fix $R>0$, $t>0$ and $M>0$ and assume that $$u_0(x) {\geqslant}1_{B(0,\,R)}(x).$$ Then there exist positive constants $c_1(R)$, $c_2(R)$, and $m_0(t,R)$ such that for all $m {\geqslant}m_0$, $$\begin{aligned}
{\mathrm{P}}(u_{s}(x){\geqslant}c_1 1_{B(0,\,R+M(t/m)^{1/\alpha})}(x)&\quad\text{for all}\quad \frac{t}{2m} {\leqslant}s{\leqslant}\frac{t}{m}\quad\text{and} \quad x\in {\mathbf{R}}^d )\\
&{\geqslant}1-c_m,\end{aligned}$$ where $$c_m:=\exp\left(-c_2m^{(\alpha-\beta)/\alpha}[\log m]^{(2\beta-\alpha)/\alpha}\right).$$

From the mild formulation of the solution of the equation, we have $$\begin{aligned}
u_{s}(x)=\int_{{\mathbf{R}}^d}p_{s}(x-y)u_0(y)\,{{\rm d}}y+\int_0^{s}\int_{{\mathbf{R}}^d}p_{s-l}(x-y)\sigma(u_l(y)) F({{\rm d}}y\,{{\rm d}}l).\end{aligned}$$ By Proposition \[initial\], there exists a $0<c_1<1$ such that for large enough $m$, $$\begin{aligned}
\int_{{\mathbf{R}}^d}p_{s}(x-y)u_0(y)\,{{\rm d}}y{\geqslant}2c_1 1_{B(0,\,R+M(t/m)^{1/\alpha})}(x)\quad \text{for all }\quad x\in {\mathbf{R}}^d \quad\text{and} \quad \frac{t}{2m} {\leqslant}s{\leqslant}\frac{t}{m}.\end{aligned}$$ By using the mild formulation and the above, we obtain $$\begin{aligned}
&{\mathrm{P}}(u_{s}(x){\leqslant}c_1 1_{B(0,\,R+M(t/m)^{1/\alpha})}(x) \quad\text{for some}\quad \frac{t}{2m}{\leqslant}s{\leqslant}\frac{t}{m})\\
&{\leqslant}{\mathrm{P}}\Big(\int_0^{s}\int_{{\mathbf{R}}^d}p_{s-l}(x-y)\sigma(u_l(y)) F({{\rm d}}y\,{{\rm d}}l)<- c_1  \quad\text{for some}\quad (s,x) \in A_{m,t,R}\Big)\\
&{\leqslant}{\mathrm{P}}\Big(\left|\int_0^{s}\int_{{\mathbf{R}}^d}p_{s-l}(x-y)\sigma(u_l(y)) F({{\rm d}}y\,{{\rm d}}l)\right|>c_1\quad\text{for some}\quad  (s,x) \in A_{m,t,R}\Big).\end{aligned}$$ The term in the above display can now be bounded by $$\begin{aligned}
c_1^{-k}\E\sup_{(s,x) \in A_{m,t,R} }\left|\int_0^{s}\int_{{\mathbf{R}}^d}p_{s-l}(x-y)\sigma(u_l(y)) F({{\rm d}}y\,{{\rm d}}l) \right|^k.\end{aligned}$$ The above in turn can be bounded using Proposition \[cty\] to obtain $$\begin{aligned}
&\E\sup_{(s,x) \in A_{m,t,R}}\left|\int_0^{s}\int_{{\mathbf{R}}^d}p_{s-l}(x-y)\sigma(u_l(y)) F({{\rm d}}y\,{{\rm d}}l) \right|^k\\
&{\leqslant}c \rho^{\tilde{\eta}k}\exp(Ak^{(2\alpha-\beta)/(\alpha-\beta)}\rho),\end{aligned}$$ where $\rho:=t/m$ and $\tilde{\eta}:=\frac{\alpha-\beta}{2\alpha}$. We now optimise the above quantity with respect to $k$ and combine all our estimates to end up with the result. See [@ChenHuang] and [@CK] for details.

Proof of Theorem \[strongcomparison\]
-------------------------------------

We next prove the strong comparison principle. We leave it to the reader ro consult [@Mueller] for the original idea and to [@ChenHuang] and [@CK] for further details.

It suffices to show that if $u_0$ has finite support then $u_t(x)>0$ for all $t>0$ and $x \in {\mathbf{R}}^d$ a.s. The general case will follow as in [@ChenHuang] and [@CK]. Assume that $u_0(x)=1_{B(0,\, R)}(x)$, for some $R>0$. Choose $M>0$, $t>0$ and $m>0$. Define for $k{\geqslant}0$, $$\begin{aligned}
A_k:=\left\{u_{s}(x){\geqslant}c_1^{k+1}1_{B^m_k}(x)\quad\text{for all} \quad s\in \Big[\frac{(2k+1)t}{2m},\,\frac{(k+1)t}{m}\Big] \quad\text{and}\quad x\in {\mathbf{R}}^d\right\}.\end{aligned}$$ For $k{\geqslant}1$, set $$\begin{aligned}
\tilde{A}_k:=\left\{u_{s}(x){\geqslant}c_1^{k+1}1_{B^m_k}(x)\quad\text{for all} \quad s\in \Big[\frac{kt}{m},\,\frac{(2k+1)t}{2m}\Big] \quad\text{and}\quad x\in {\mathbf{R}}^d\right\}.\end{aligned}$$ Finally, we define $$\begin{aligned}
\tilde{A}_0:=\left\{u_{\frac{t}{2m}}(x){\geqslant}c_11_{B^m_0}(x)\quad\text{for all} \quad x\in {\mathbf{R}}^d\right\},\end{aligned}$$ where $B^m_k=B(0, \, R+kM(t/m)^{1/\alpha})$ and $c_1$ is as in Proposition \[stoch\]. It is clear that if $\alpha>1$, then as $m$ gets large, the sets $B^m_k$ cover the whole space. For $\alpha=1$, the sets $B^m_k$ cover $B(0, \, R+ Mt^{1/\alpha})$. Since $M$ is arbitrary, we have for any $M$ large enough, $$\begin{split}\label{h1}
&{\mathrm{P}}(u_{s}(x)>0\quad\text{for all}\quad t/2{\leqslant}s{\leqslant}t\quad\text{and}\quad x\in B(0,\, M/2))\\
&\qquad \qquad {\geqslant}\lim_{m\rightarrow \infty}{\mathrm{P}}(\cap_{0{\leqslant}k{\leqslant}m-1}A_k\cap \tilde{A}_k)\\
\end{split}$$ Proposition \[stoch\] can be used to obtain $$\label{h2}
{\mathrm{P}}(A_0){\geqslant}1-c_m,$$ whenever $m$ is large enough. On the other hand, on the event $A_{k-1}$, $k {\geqslant}1$, $$u_{\frac{kt}{m}}(x){\geqslant}c_1^{k}1_{B^m_{k-1}}(x), \quad \text{ for all } \; x\in {\mathbf{R}}^d.$$ By the Markov property, $\{u_{s+\frac{kt}{m} }(x), s {\geqslant}0, x \in {\mathbf{R}}^d\}$ solves (\[main-eq\]) with the time-shifted noise $\dot{F}_{k}(s,\,x):=\dot{F}(s+\frac{kt}{m},\,x)$ starting from $u_{\frac{kt}{m}}(x)$. Let $\{v^{(k)}_{s}(x), s {\geqslant}0, x \in {\mathbf{R}}^d\}$ be the solution to (\[main-eq\]) with the time-shifted noise $\dot{F}_{k}(s,\,x)$, $\sigma$ replaced by $\sigma_k(x)=c_1^{-k} \sigma(c_1^kx)$, and initial condition $1_{B^m_{k-1}}(x)$. On one hand, by Proposition \[stoch\] we get that $$\begin{aligned}
{\mathrm{P}}\big(v_{s}^{(k)}(x){\geqslant}c_1 1_{B^m_k}(x)\quad\text{for all}\quad s\in \big[\frac{t}{2m},\frac{t}{m}\big]\quad \text{and}\quad x\in {\mathbf{R}}^d\big){\geqslant}1-c_m,\end{aligned}$$ whenever $m$ is large enough. On the other hand, by Markov property and the weak comparison principle (Theorem \[weakc\]) we see that on $A_{k-1}$, $u_{s+kt/m}(x){\geqslant}c_1^kv_{s}^{(k)}(x)$ for all $x\in {\mathbf{R}}^d$ and $s{\geqslant}0$. We therefore have for all $k=1,\ldots, m-1$, $$\begin{aligned}
{\mathrm{P}}(A_k|\mathcal{F}_{kt/m}){\geqslant}1-c_m\quad \text{on}\quad A_{k-1}.\end{aligned}$$ And hence $$\label{h3}
{\mathrm{P}}(A_k|A_{k-1}\cap\cdots \cap A_0){\geqslant}1-c_m.$$ Similarly, we have $$\label{h4}
{\mathrm{P}}(\tilde{A}_k|\tilde{A}_{k-1}\cap\cdots \cap \tilde{A}_0){\geqslant}1-c_m.$$ Combining the above estimates as in [@ChenHuang] or [@CK], we have $$\begin{aligned}
{\mathrm{P}}(u_{s}(x)>0\quad\text{for all}\quad t/2{\leqslant}s{\leqslant}t\quad\text{and}\quad x\in B(0,\, M/2)){\geqslant}2(1-c_m)^{m}-1 \rightarrow 1,\end{aligned}$$ as $m \rightarrow \infty$. This finishes the proof since $t>0$ is arbitrary and $M$ can be taken as large as possible.

Proof of Theorem \[positivity\]
-------------------------------

The proof is very similar to those in [@ChenHuang], [@CK] and [@CJKCorr], using the strong Markov property (Lemma \[markov\]) and the weak comparison principle (Theorem \[weakc\]). So we omit it.

****

As mentioned in the introduction, the above comparison theorem and strict positivity results are shown under the assumption that the initial conditions are bounded functions. But this can be relaxed to a wider class of initial conditions as studied in [[@ChenHuang]]{.nodecor} and [[@CK]]{.nodecor}. We leave it to the reader to check the details.
