\section{Introduction and Main Results} \label{sec:intro}

The main topic of this paper is the study of the Besov regularity of L\'evy white noises.
We are especially interested in identifying the critical local smoothness and the critical asymptotic growth rate of those random processes for any integrability parameter.
In a nutshell, our contributions are as follows: 

\begin{itemize}

	\item \emph{Wavelet methods for L\'evy white noises.} First appearing in the eighties, especially in the works of Y. Meyer~\cite{MeyerWaO} and I. Daubechies~\cite{Daubechies92}, wavelet techniques became established are primary tools in functional analysis.
	As such, they can naturally be used for the study of random processes, as it was done, for instance, for fractional Brownian motion~\cite{Meyer1999wavelets}, S$\alpha$S processes~\cite{Pad2015optimality}, and for the study of stochastic partial differential equations~\cite{Hairer2014theory,Hairer2017reconstruction}.
	In this paper, we demonstrate that such wavelet techniques are also adapted to the analysis of L\'evy white noises.
	In particular, all our results are derived using the wavelet characterization of weighted Besov spaces.
	
	\item \emph{New positive results.} We determine in which weighted Besov spaces is a given L\'evy white noise $w$ is. We are able to improve known results; in particular, for the growth properties of the noise. This requires the identification of a new index associated to $w$, characterized by moments properties.
	
	\item \emph{Negative regularity results.} We also determine in which Besov spaces a L\'evy white noise is \emph{not}. This requires a more evolved analysis as compared to positive results and had, to the best of our knowledge, only received limited attention of researchers in the past. 
	
	\item \emph{Critical local smoothness and asymptotic rate.} The combination of positive and negative results allows   determining the critical Besov exponents of the L\'evy white noises, both for the local smoothness and the asymptotic behavior. The results are summarized in Theorem \ref{maintheo:regdecaylimit}. Two consequences are the characterization of the Sobolev and H\"older-Zigmund regularities of L\'evy white noises in Corollary \ref{coro:SHZ}.
	
\end{itemize}

	\subsection{Local Smoothness and Asymptotic Rate of Tempered Generalized Functions}

We construct random processes as random elements in the space ${\mathcal{S}}'( \mathbb{R}^d)$ of tempered generalized functions from $ \mathbb{R}^d$ to $ \mathbb{R}$ (see Section \ref{subsec:noises}). 
We will therefore describe their local and asymptotic properties as we would do for a (deterministic) tempered generalized function.
To do so, we rely on the family of weighted Besov spaces, which are embedded in ${\mathcal{S}}'( \mathbb{R}^d)$ and allows to jointly study the local smoothness and the asymptotic behavior of a generalized function. 
 
Besov spaces are denoted by $B_{p,q}^{\tau}( \mathbb{R}^d)$, where $\tau \in  \mathbb{R}$ is the \emph{smoothness}, $p \in (0,\infty]$ the \emph{integrability}, and $q \in (0,\infty]$ a secondary parameter.
In this paper, we focus on the case $p=q$. 
We use the simplified notation $B_{p,p}^{\tau}( \mathbb{R}^d) = B_p^{\tau}( \mathbb{R}^d)$ for this spaces which are sometimes refer to as Sobolev-Sobodeckij spaces. 
We say that $f$ is in the weighted Besov space $B_p^{\tau}( \mathbb{R}^d;\rho)$ with \emph{weight} $\rho \in  \mathbb{R}$ if $\langle \cdot \rangle^{\rho} \times f$ is in the classical Besov space $B_p^{\tau}( \mathbb{R}^d)$, with the notation $\langle \bm{x} \rangle = (1 + \lVert \bm{x} \rVert^2)^{1/2}$. We precisely define weighted Besov spaces in Section \ref{subsec:besovspaces} in terms of wavelet expansions.  For the time being, it is sufficient to remember that the space of tempered generalized functions satisfies \cite[Proposition 1]{Kabanava2008tempered}
\begin{equation} \label{eq:decomposeSprime}
	{\mathcal{S}}'( \mathbb{R}^d) = \bigcup_{\tau, \rho \in  \mathbb{R}} B_p^{\tau}( \mathbb{R}^d;\rho).
\end{equation} 
Ideally,  we aim at identifying in which weighted Besov space  a given $f \in {\mathcal{S}}'( \mathbb{R}^d)$ is. The relation \eqref{eq:decomposeSprime} implies that for any $p$, there exists some $\tau,\rho$ for which this is true.
For $p$ fixed, Besov spaces are continuously embedded in the sense that, for $\tau, \tau_1,  \tau_2$ and $\rho, \rho_1 , \rho_2$  such that $\tau_1 \leq \tau_2$ and $\rho_1 \leq \rho_2$, we have
\begin{align*}
	B_{p}^{\tau_2}( \mathbb{R}^d;\rho)  \subseteq B_p^{\tau_1}( \mathbb{R}^d;\rho)  \quad \text{and} \quad
	B_{p}^{\tau}( \mathbb{R}^d;\rho_2)  \subseteq B_p^{\tau}( \mathbb{R}^d;\rho_1).
\end{align*}
To characterize the properties of $f \in {\mathcal{S}}'( \mathbb{R}^d)$, the key is then to determine the two critical exponents $\tau_p(f) \in (-\infty, \infty]$ and $\rho_p(f) \in (-\infty, \infty]$ such that:
\begin{itemize}
	\item if $\tau < \tau_p(f)$ and $\rho < \rho_p(f)$, then $f \in B_p^\tau( \mathbb{R}^d;\rho)$; while
	\item  if $\tau > \tau_p(f)$ or $\rho > \rho_p(f)$, then $f \notin B_p^\tau( \mathbb{R}^d;\rho)$.
\end{itemize}
The case $\tau_p(f) = \infty$ corresponds to smooth functions, and \mbox{$\rho_p(f) = \infty$} means that $f$ is rapidly decaying. 
The quantity $\tau_p(f)$ measures the \emph{local smoothness} and $\rho_p(f)$ the \emph{asymptotic rate} of $f$ for the integrability   $p$. When $\rho_p(f)<0$ (which will be the case for L\'evy white noise), we talk about the \emph{asymptotic growth rate} of $f$.

	\subsection{Local Smoothness and Asymptotic Rate of L\'evy White Noises} \label{subsec:taurhointro}

The complete family of L\'evy white noises, defined as random elements in the space of generalized functions, was introduced by I.M. Gel'fand and N.Y. Vilenkin \cite{GelVil4}. We briefly  recap the concepts required to state our main results. A more complete exposition is given in Section \ref{sec:maths}.

Our main contributions concern the localization of a  L\'evy white noise $w$ in weighted Besov spaces. It includes positive   ($w$ is almost surely \emph{in} a given Besov space) and negative  ($w$ is almost surely \emph{not in} a given Besov space) results. 
In order to characterize the  local smoothness $\tau_p(w)$ and the asymptotic rate $\rho_p(w)$, we now introduce some notations.  
Let $X = \langle w , \One_{[0,1]^d} \rangle$ be the random variable corresponding to the integration of the L\'evy white noise $w$ over the domain $[0,1]^d$. The \emph{L\'evy exponent} $\Psi$ of $w$ is the logarithm of the characteristic function of $x$; that is, for every $\xi \in  \mathbb{R}$, 
\begin{equation}
	\Psi(\xi) = \log \mathbb{E} \left[ \mathrm{e}^{\mathrm{i} \xi  \langle w , \One_{[0,1]^d} \rangle } \right]
 =  \log \mathbb{E} \left[ \mathrm{e}^{\mathrm{i} \xi X } \right] .
\end{equation}

	We associate to a L\'evy white noises its \emph{Blumenthal-Getoor indices}, defined as 
\begin{align}
	\beta_\infty &= \inf \setb{p>0}{\underset{\abs{\xi} \rightarrow \infty}{\lim} \frac{\abs{\Psi(\xi)}}{\abs{\xi}^p}= 0},  \label{eq:indass} \\
	\ubar{\beta}_\infty &= \inf \setb{p>0}{\underset{\abs{\xi} \rightarrow \infty}{\lim\inf} \frac{\abs{\Psi(\xi)}}{\abs{\xi}^p}= 0}. \label{eq:indasi}	
\end{align}
The distinction is that $\beta_\infty$ considers the limit, while $\ubar{\beta}_\infty$ deals with the inferior limit. In general, one has that $0 \leq \ubar{\beta}_\infty \leq \beta_\infty \leq 2$. The Blumenthal-Getoor indices are linked to the local behavior of L\'evy processes and L\'evy white noises (see Section \ref{subsec:indices} for more details and references to the literature).
	In addition, we introduce the \emph{moment index} of $w$
\begin{align}
	p_{\max} &= \sup \setb{p>0}{\mathbb{E}[\lvert\langle w ,\One_{[0,1]^d}\rangle \rvert^p] < \infty},	 \label{eq:indmom}
\end{align}
white noise which is closely related---but not identical--- to the Pruitt index in general (see Section \ref{subsec:indices}).
As we shall see, $p_{\max}$ fully characterizes the asymptotic rate of the L\'evy.


The class of L\'evy white noise is very general. Specific examples are the Gaussian and  compound Poisson white noises, which will receive a special treatment thereafter. We summarize the results of this paper in Theorem \ref{maintheo:regdecaylimit}. 

\begin{theorem} \label{maintheo:regdecaylimit}
Consider a L\'evy white noise $w$ with Blumenthal-Getoor indices $0\leq \ubar{\beta}_\infty \leq \beta_\infty \leq 2$ and moment index $0 < p_{\max} \leq \infty$. We fix $0<p\leq \infty$.
\begin{itemize}
	\item If $w$ is Gaussian, then, almost surely,
	\begin{equation} \label{eq:GNtaurho}
	\tau_p(w) = - \frac{d}{2} \text{ and } \rho_p(w) = - \frac{d}{p}.
	\end{equation}
	\item If $w$ is compound Poisson, then, almost surely,
	\begin{equation} \label{eq:PNtaurho}
	\tau_p(w) =  \frac{d}{p} -d   \text{ and } \rho_p(w) = - \frac{d}{\min(p,p_{\max})}.
	\end{equation}		
	\item If $w$ is non-Gaussian, then, almost surely, if $p\in (0,2)$, $p$ is an even integer, or $p = \infty$,
	\begin{equation} \label{eq:LNtaurho}
  \frac{d}{\max(p,\beta_\infty)} -d    \leq \tau_p(w) \leq  \frac{d}{\max(p,\ubar{\beta}_\infty)} -d   \text{ and } \rho_p(w) = - \frac{d}{\min(p,p_{\max})}.
	\end{equation}					
\end{itemize}
\end{theorem}

Theorem \ref{maintheo:regdecaylimit} provides the full answer to the non-Gaussian scenario with $0<p\leq 2$ when $\ubar{\beta}_\infty=\beta_\infty$. This covers the case for most of the L\'evy white noises used in practice, including symmetric-$\alpha$-stable or symmetric-gamma.
For $p>2$, we have full  results for even integers. Note that the smoothness is  $\tau_p(w) = d / p - d$ in this case (because $\beta_\infty \leq 2$). 
In Section \ref{sec:general}, we refine Theorem  \ref{maintheo:regdecaylimit} by showing that the relation $ \tau_p(w) \leq  \frac{d}{\max(p,\ubar{\beta}_\infty)} -d$ and $\rho_p(w) = - \frac{d}{\min(p,p_{\max})}$ are valid for any $p>0$. We conjecture that \eqref{eq:LNtaurho} is actually true without restriction on $p$.
		Two direct consequences are the identification of the Sobolev ($p=2$) and H\"older-Zigmund ($p= \infty$) regularity of L\'evy white noises.
	
	\begin{corollary} \label{coro:SHZ}
		Let $w$ be a L\'evy white noise in ${\mathcal{S}}'( \mathbb{R}^d)$ with indices $p_{\max}, \beta_\infty, \ubar{\beta}_\infty$. Then, the Sobolev local smoothness and asymptotic growth rate ($p=2$) are 
		\begin{equation}
			\tau_2(w) = - \frac{d}{2} \quad \text{and} \quad \rho_2(w) = - \frac{d}{\min ( 2, p_{\max} )}.
		\end{equation}
	Moreover, the H\"older-Zigmund local smoothness and asymptotic growth rate ($p=\infty$) are 
	\begin{align}
			\tau_\infty(w) = - \frac{d}{2} \quad &\text{and} \quad \rho_\infty(w) = 0  \quad \text{if } w \text{ is Gaussian, and} \\
			\tau_\infty(w) = -  d \quad &\text{and} \quad \rho_\infty(w) = - \frac{d}{ p_{\max} } \quad \text{otherwise}.
	\end{align}
	\end{corollary}
	
	\begin{proof}
	The case $p=2$ is directly deduced from Theorem \ref{maintheo:regdecaylimit} and the relation $\beta_\infty \leq 2$. For $p= \infty$, we use again Theorem \ref{maintheo:regdecaylimit} with $p=2m$, together with the embedding relations (Proposition \ref{prop:besovembed}) from which we know that $\tau_w(\infty) = \lim_{m\rightarrow \infty} \tau_w(2m)$ and $\rho_w(\infty) = \lim_{m\rightarrow \infty} \rho_w(2m)$. 
	\end{proof}	
	
		Theorem \ref{maintheo:regdecaylimit} also allows deducing the following results on the local smoothness of L\'evy processes.

	\begin{corollary} \label{coro:levyprocess}
		Let $X$ be a L\'evy process with Blumenthal-Getoor indices $0 \leq \ubar{\beta}_\infty \leq \beta_\infty \leq 2$. Then, we have almost surely that, for any $0<p\leq \infty$,
	\begin{align}
			\tau_p(X) & = \frac{1}{2} \quad \text{if } X \text{ is Gaussian, and } \\
		 	\tau_p(X) &=\frac{1}{p}  \quad  \text{if } X \text{ is compound Poisson}.
	\end{align}
	In the general case, we have almost surely that, for $p \in (0,2)$, $p$ is an even integer, or $p=\infty$, 
	\begin{equation}
	\frac{1}{\max(p,\beta_\infty)}\leq 	\tau_p(X)  \leq \frac{1}{\max(p,\ubar{\beta}_\infty)} .
	\end{equation}
	\end{corollary}

\begin{proof}
	A L\'evy white noise $w$ is the weak derivative of the corresponding L\'evy process $X$ with identical L\'evy exponent. 
	This well-known fact has been rigorously shown in the sense of generalized random processes in \cite{Dalang2015Levy}.
	A direct consequence is that $\tau_p(X) = \tau_p(w) + 1$, where $w = X'$. 
	Then, Corollary \ref{coro:levyprocess} is a reformulation of the local smoothness results of Theorem \ref{maintheo:regdecaylimit} with $d=1$. 
\end{proof}
	
	

	\subsection{Related Works on L\'evy Processes and L\'evy White Noises}

In this section, for comparison purposes, we reinterpret all the results in terms of the critical smoothness  and asymptotic rate of the considered random processes.

	\textit{L\'evy processes.}
	Most of the attention has been  so far  devoted to classical L\'evy processes $s :  \mathbb{R} \rightarrow  \mathbb{R}$. 
	The Brownian motion was studied in \cite{Sjogren1982riemann,Roynette1993,Ciesielski1993quelques}. The work of \cite{Ciesielski1993quelques} also contains results on the Besov regularity of fractional Brownian motions and S$\alpha$S processes. By exploiting the self-similarity of the stable processes, Ciesielski \emph{et al.} obtained the following results  the Gaussian \cite[Theorem IV.3]{Ciesielski1993quelques} and  stable non-Gaussian \cite[Theorem VI.1]{Ciesielski1993quelques} scenarios:
	\begin{align}
		\frac{1}{2} \leq \tau_p(X_{\mathrm{Gauss}})  \leq \frac{1}{\min(2,p)} , \label{eq:BMtau} \\
		\frac{1}{\max(p,\alpha)} \leq \tau_p(X_{\alpha} ) \leq  \frac{1}{\min(p,\alpha)}, \label{eq:SaStau}
	\end{align}	
for any $1 \leq p < \alpha$, where $X_{\mathrm{Gauss}}$ is the Brownian motion and $X_{\alpha}$ is the S$\alpha$S process with parameter $1 < \alpha < 2$. 

The complete family of L\'evy processes---and more generally of L\'evy-type processes---has been considered by R. Schilling in a series of papers~\cite{Schilling1997Feller,Schilling1998growth,Schilling2000function} synthesized in \cite[Chapter V]{Bottcher2014levy} and by V. Herren~\cite{Herren1997levy}.
To summarize, Schilling has shown that, for a L\'evy process $X$ with indices $0 \leq \beta_\infty \leq 2$ and $0<p_{\max}$,
	\begin{align}
		\frac{1}{\max(p,\beta_\infty)} &\leq \tau_{p}(X) \leq \frac{1}{p}, \label{eq:boundtauXSchilling} \\
		 -\frac{1}{p} - \frac{1}{\min(p_{\max},2)} &\leq \rho_p(X). \label{eq:boundrhoXSchilling}
	\end{align}	

	 We see that  \eqref{eq:BMtau}, \eqref{eq:SaStau}, and \eqref{eq:boundtauXSchilling}  are consistent with Corollary \ref{coro:levyprocess}. 
	Moreover, our results provide an improvement by showing that the lower bounds of  \eqref{eq:BMtau} and \eqref{eq:SaStau} are actually sharp.
		Finally, we significantly improve the upper bound of \eqref{eq:boundtauXSchilling} for general L\'evy processes. 
	

 In contrast to the smoothness, the growth rate \eqref{eq:boundrhoXSchilling} of the L\'evy process $X$  does not seem to be related to the one of the L\'evy white noise $w = X'$   by a simple relation. In particular, the rate of $X$ is expressed in terms of the Pruitt index  $ \beta_0 = \min ( p_{\max} , 2)$, conversely to $p_{\max}$ for $w$ (see Section \ref{subsec:indices}). 
		This needs to be confirmed by a precise estimation of $\tau_p(X)$ for which only a lower bound is known. Our conjecture is that \eqref{eq:boundrhoXSchilling} is sharp; that is, $\rho_p(X) = -\frac{1}{p} - \frac{1}{\min(p_{\max},2)}$.
		
		\textit{L\'evy white noises.}
Veraar extensively studied  the local Besov regularity of the $d$-dimensional Gaussian white noise. As a corollary of  \cite[Theorem 3.4]{Veraar2010regularity}, one then deduces that $\tau_{p}(w_{\mathrm{Gauss}}) = - d /2$. 
This work is based on the Fourier series expansion of the process, which is specific to the Gaussian case.

In our own works, we have investigated the question for general L\'evy white noises in dimension $d$ in the periodic \cite{Fageot2017besov} and global settings \cite{Fageot2017multidimensional}. We had obtained the lower bounds
	\begin{align}
		\frac{d}{\max(p,\beta_\infty)} - d \leq &\tau_{p}(w)  \text{ and }\label{eq:boundtauwus} \\
		- \frac{d}{\min(p,p_{\max},2)} \leq &\rho_p(w). \label{eq:boundrhowus}
	\end{align}	
	These estimates are improved by	Theorem \ref{maintheo:regdecaylimit}, which now provides an upper bound for $\tau_{p}(w)$ and shows that \eqref{eq:boundtauwus} is sharp when $\beta_\infty = \ubar{\beta}_\infty$. 
	It is also worth noticing that the lower bound of \eqref{eq:boundrhowus} is sharp if and only if $p_{\max} \leq 2$; in particular, when the L\'evy white noise has an infinite variance ($p_{\max} < 2$).

	\subsection{Sketch of Proof and the Role of Wavelet Methods}
	
	Our techniques are based on the wavelet characterization of Besov spaces as presented by H. Triebel in~\cite{Triebel2008function}. We shall wee that the wavelets are especially relevant to the analysis of L\'evy white noises.
	
	We  briefly present the strategy of the proof of Theorem \ref{maintheo:regdecaylimit} in the simplified scenario of $d=1$.  	
	The more general case, $d\geq 1$, is analogous (up to some normalization factors) and will be comprehensively introduced in Section \ref{sec:maths}.
	Let $(\psi_M,\psi_F)$ be the mother and father Daubechies wavelet of a fixed order (the choice of the order has no influence on the results), respectively. 
	For $j\geq 0$ and $k\in  \mathbb{Z}$, we define the rescaled and shifted functions $\psi_{F,k} = \psi_F (\cdot - k)$ and $\psi_{j,M,k} = 2^{j/2} \psi_M (2^j \cdot - k)$.
	Then, the family $(\psi_{F,k})_{k \in  \mathbb{Z}} \cup (\psi_{j , M, k})_{j\geq 0 , k\in  \mathbb{Z}}$ forms an orthonormal basis of $L_2( \mathbb{R}^d)$.  
	For a given one-dimensional L\'evy white noise $w$, one considers the family of random variables 
	\begin{equation} \label{eq:waveletrandomvariable}
	(\langle w , \psi_{F,k}\rangle )_{k\in  \mathbb{Z}} \cup (\langle w , \psi_{j,M,k} \rangle)_{j \geq 0 , k\in  \mathbb{Z}}.
	\end{equation}
	We then have that $w = \sum_{k\in  \mathbb{Z}} \langle w ,\psi_{F,k}\rangle \psi_{F,k}+\sum_{j \geq 0}\sum_{k\in  \mathbb{Z}} \langle w ,\psi_{j,M,k}\rangle \psi_{j,M,k}$, for which the convergence is almost sure in ${\mathcal{S}}'( \mathbb{R})$.
	
	Then, for $0<p\leq \infty$ and $\tau, \rho\in \mathbb{R}$, the random variable 
\begin{equation} \label{eq:methodsbesov}
	\lVert w \rVert_{B_p^\tau( \mathbb{R} ;\rho)} =  \left(
		\sum_{k\in  \mathbb{Z}}  \langle k \rangle^{\rho}  \lvert \langle w , \psi_{F,k} \rangle \rvert^p
		+
		\sum_{j\geq 0} 2^{j(\tau p -1 + \frac{p}{2})} 
		\sum_{k\in  \mathbb{Z}} \langle 2^{-j} k \rangle^{\rho}    \lvert \langle w ,\psi_{j,M,k}\rangle \rvert^p \right)^{1/p}
\end{equation}
is well-defined, and takes values in $[0,\infty]$. We adapt \eqref{eq:methodsbesov} in the usual manner when $p=\infty$. We recognize the Besov (quasi-)norm that characterizes the Besov localization of the L\'evy white noise. This means that $w$ is a.s. (almost surely) in $B_p^\tau( \mathbb{R} ;\rho)$ if and only if $\lVert w \rVert_{B_p^\tau( \mathbb{R} ;\rho)} < \infty$  a.s., and  a.s.  not in  $B_p^\tau( \mathbb{R} ;\rho)$ if and only if $\lVert w \rVert_{B_p^\tau( \mathbb{R} ;\rho)}  = \infty$ a.s.

We then fix $0<p\leq \infty$. We assume that we have guessed the values $\tau_p(w)$ and $\rho_p(w)$  introduced in Section \ref{sec:intro}. Here are the main steps leading to the proof that these values are the correct ones.

\begin{itemize}
	
	\item For $ \tau < \tau_p(w) $ and $\rho < \rho_p(w)$, we show that $\lVert w \rVert_{B_p^\tau( \mathbb{R};\rho)} < \infty$ a.s. For $p < p_{\max}$ (see \eqref{eq:indmom}), we establish the strongest result  $\mathbb{E} [\lVert w \rVert_{B_p^\tau( \mathbb{R};\rho)}^p] < \infty$. This requires moment estimates for the wavelet coefficients of L\'evy white noise; that is, a precise estimation of the behavior of $\mathbb{E}[  \lvert \langle w , \psi_{j,G, k} \rangle \rvert^p ]$ as $j$ goes to infinity. When $p > p_{\max}$,  the random variables $\langle w , \psi_{j,G, k} \rangle$ have an infinite $p$th moment and the present method is not applicable. In that case, we actually deduce the result using   embeddings relations between Besov spaces. It appears that this approach is sufficient to obtain sharp   results. 
	
	\item For $\tau > \tau_p(w)$, we show that $\lVert w \rVert_{B_p^\tau( \mathbb{R};\rho)} = \infty$ a.s.
		To do so, we only consider the mother wavelet and truncate the sum over $k$ to yield the lower bound
		\begin{equation} \label{eq:principlesmoothness}
		\lVert w \rVert_{B_p^\tau( \mathbb{R};\rho)}^p \geq 
		 C \sum_{j\geq 0} 2^{j(\tau p -d + \frac{dp}{2})} 
		\sum_{0\leq k_1, \ldots , k_d < 2^{j}}  \lvert \langle w , \psi_{j,M, k} \rangle \rvert^p
		\end{equation}
		for some constant $C$ such that   $\langle 2^{-j} k  \rangle^{\rho p} \geq C$ for every $j \geq 0$ and $0\leq k < 2^{j}$.
 		We then need to show  that the wavelet coefficients $\langle w , \psi_{j,G, k} \rangle$ cannot be too small altogether using Borel-Cantalli-type arguments. Typically, this requires to control the evolution  of quantities such as $\mathbb{P}(\abs{\langle w ,\psi_{j,G,k} \rangle } > x)$ with respect to $j$ and is again based on moment estimates.  
		
	\item For $\rho > \rho_p(w)$, we show again that $\lVert w \rVert_{B_p^\tau( \mathbb{R};\rho)} = \infty$ a.s.
		This time, we only consider the father wavelet  in \eqref{eq:methodsbesov}  and use the lower bound 
\begin{equation}  \label{eq:principledecay}
	\lVert w \rVert_{B_p^\tau( \mathbb{R};\rho)}^p  \geq	\sum_{k\in  \mathbb{Z}} \langle  k  \rangle^{\rho p} \lvert \langle w , \psi_{0, F,k} \rangle \rvert^p.
\end{equation}			
	A Borel-Cantelli-type argument is again used to show that the $ \lvert \langle w , \psi_{0, F,k} \rangle \rvert $ cannot be too small altogether, and that the Besov norm is a.s.  infinite.
\end{itemize}

	The rest of the paper is dedicated to the proof of Theorem \ref{maintheo:regdecaylimit}. The required  mathematical concepts---L\'evy white noises as generalized random processes and weighted Besov spaces---are laid out in Section \ref{sec:maths}.
	In Sections \ref{sec:Gaussian}, \ref{sec:Poisson}, and \ref{sec:finitemoments}, we consider  the case of Gaussian noises, compound Poisson noises, and finite moments L\'evy white noises, respectively. 
	The general case for any L\'evy white noise is deduced in Section \ref{sec:general}. Finally, we discuss our results and give important examples in Section \ref{sec:blabla}. 
	
\section{Preliminaries: L\'evy White Noises and Weighted Besov Spaces} \label{sec:maths}

	\subsection{L\'evy White Noises as Generalized Random Processes} \label{subsec:noises}
	
Let	${\mathcal{S}}( \mathbb{R}^d)$ be the space of rapidly decaying smooth functions from $ \mathbb{R}^d$ to $ \mathbb{R}$. It is endowed with its natural Fr\'echet nuclear topology \cite{Treves1967}. Its topological dual is the space of tempered generalized functions ${\mathcal{S}}'( \mathbb{R}^d)$.
	We shall define random processes as random elements of the space ${\mathcal{S}}'( \mathbb{R}^d)$. This allows for a proper definition of L\'evy white noises even if they do not have a pointwise interpretation (they can only be described by their effect into test functions).
	
	The space ${\mathcal{S}}'( \mathbb{R}^d)$ is endowed with the strong topology and $\mathcal{B}({\mathcal{S}}'( \mathbb{R}^d))$ denotes the Borelian $\sigma$-field for this topology. Throughout the paper, we fix a probability space $(\Omega, \mathcal{F}, \mathscr{P})$.
	
	\begin{definition} \label{def:GRP}
		A \emph{generalized random process} $s$ is a measurable function from $(\Omega, \mathcal{F})$ to $({\mathcal{S}}'( \mathbb{R}^d), \mathcal{B}({\mathcal{S}}'( \mathbb{R}^d))$. Its \emph{probability law} is the probability measure on ${\mathcal{S}}'( \mathbb{R}^d)$ defined for $B \in \mathcal{B}({\mathcal{S}}'( \mathbb{R}^d))$ by 
		\begin{equation}
			\mathscr{P}_s ( B ) = \mathscr{P} \{ \omega \in \Omega , \ s(\omega) \in B \}.
		\end{equation}
		The \emph{characteristic functional} of $s$ is the functional ${\widehat{\mathscr{P}}}_s : {\mathcal{S}}( \mathbb{R}^d) \rightarrow  \mathbb{C}$ such that 
		\begin{equation}
			{\widehat{\mathscr{P}}}_s(\varphi) = \int_{{\mathcal{S}}'( \mathbb{R}^d)} \mathrm{e}^{\mathrm{i} \langle u , \varphi \rangle } \mathrm{d} \mathscr{P}_s (u).
		\end{equation}
	\end{definition}
	
The characteristic functional is continuous, positive-definite over ${\mathcal{S}}( \mathbb{R}^d)$, and normalized such that ${\widehat{\mathscr{P}}}_s(0) = 1$. The converse of this result is also true: if ${\widehat{\mathscr{P}}}$ is a continuous and positive-definite functional over ${\mathcal{S}}( \mathbb{R}^d)$ such that ${\widehat{\mathscr{P}}}(0) = 1$, it is the characteristic functional of a generalized random process in ${\mathcal{S}}'( \mathbb{R}^d)$. This  is known as the Minlos-Bochner theorem \cite{GelVil4,Minlos1959generalized}. 
It means in particular that one can define generalized random processes via the specification of their characteristic functional. Following Gelfand and Vilenkin, we use this principle to introduce L\'evy white noises.

We consider functionals of the form ${\widehat{\mathscr{P}}}(\varphi) = \exp \left( \int_{ \mathbb{R}^d} \Psi ( \varphi ( \bm{x} ) ) \mathrm{d} \bm{x} \right)$.
It is known \cite{GelVil4} that ${\widehat{\mathscr{P}}}$ is a characteristic functional over the space ${\mathcal{D}}( \mathbb{R}^d)$  of compactly supported smooth functions, if and only if the function $\Psi :  \mathbb{R} \rightarrow  \mathbb{C}$ is continuous, conditionally positive-definite, with $\Psi(0) = 0$ \cite{GelVil4}. A function $\Psi$ satisfying these conditions is called a \emph{L\'evy exponent} and can be decomposed according to the L\'evy-Khintchine theorem \cite{Sato1994levy} as 
		\begin{equation}
		 	\Psi(\xi) = \mathrm{i} \mu \xi - \frac{\sigma^2 \xi^2}{2} + \int_{ \mathbb{R}} (\mathrm{e}^{\mathrm{i} \xi t} - 1 - \mathrm{i} \xi t \One_{\abs{t} \leq 1} ) \mathrm{d} \nu (t)
		\end{equation}
		where $\mu \in  \mathbb{R}$, $\sigma^2 \geq 0$, and $\nu$ is a \emph{L\'evy measure}; that is, a positive measure on $ \mathbb{R}$ such that $\nu\{0\} = 0$ and $\int_{ \mathbb{R}}  \inf( 1 , t^2 ) \mathrm{d} \nu ( t) < \infty$. The triplet $(\mu,\sigma^2,\nu)$ is unique and called the \emph{L\'evy triplet} of $\Psi$. 

In our case, we are only interested by defining random processes, especially L\'evy white noises, over ${\mathcal{S}}'( \mathbb{R}^d)$. This requires an adaptation of the theory of Gelfand and Vilenkin. 
We say that the L\'evy exponent $\Psi$ satisfies the \emph{$\epsilon$-condition} if there exists some $\epsilon>0$ such that $\int_{ \mathbb{R}}  \inf( \abs{t}^\epsilon , t^2 ) \mathrm{d} \nu ( t) < \infty$, with $\nu$ the L\'evy measure of $\Psi$. This condition is   fulfilled by all the L\'evy measures encountered in practice. Then, the functional ${\widehat{\mathscr{P}}}(\varphi) = \exp \left( \int_{ \mathbb{R}^d} \Psi ( \varphi ( \bm{x} ) ) \mathrm{d} \bm{x} \right)$ is a  characteristic functional over ${\mathcal{S}}( \mathbb{R}^d)$ if and only if $\Psi$ is a L\'evy exponent satisfying the $\epsilon$-condition. The sufficiency is proved in \cite{Fageot2014} and the necessity in \cite{Dalang2015Levy}. 

\begin{definition} \label{def:levynoise}
	A \emph{L\'evy white noise in ${\mathcal{S}}'( \mathbb{R}^d)$} (or simply L\'evy white noise) is a generalized random process $w$ with characteristic functional of the form
	\begin{equation}
		{\widehat{\mathscr{P}}}_w(\varphi) = \exp\left(  \int_{ \mathbb{R}^d} \Psi( \varphi(\bm{x} ) ) \mathrm{d} \bm{x}\right) 
	\end{equation}
	for every $\varphi \in {\mathcal{S}}( \mathbb{R}^d)$, where $\Psi$ is a L\'evy exponent satisfying the $\epsilon$-condition.
	
	The L\'evy triplet of $w$ is denoted by $(\mu,\sigma^2, \nu)$. Then, we say that $w$ is a \emph{Gaussian white noise} if $\nu = 0$, \emph{compound Poisson white noise} if $\sigma^2 = 0$ and $\nu = \lambda \mathscr{P}$, with $\lambda > 0$ and $\mathscr{P}$ a probability measure such that $\mathscr{P} \{0\}= 0$, and a  \emph{finite moment white noise} if $\mathbb{E}[\abs{\langle w, \varphi \rangle}^p] < \infty$ for any $\varphi \in {\mathcal{S}}( \mathbb{R}^d)$ and $p>0$.
\end{definition}

L\'evy white noises are stationary and independent at every point, meaning that $\langle w , \varphi_1 \rangle$ and $\langle w , \varphi_2 \rangle$  are independent as soon as $\varphi_1$ and $\varphi_2 \in {\mathcal{S}}( \mathbb{R}^d)$ have disjoint supports. 

	One can extend the space of test functions that can be applied to the L\'evy white noise. This is done by approximating a test function $\varphi$ with functions in ${\mathcal{S}}( \mathbb{R}^d)$ and showing that the underlying sequence of random variables converges in probability to a random variable that we denote by $\langle w , \varphi \rangle$. This principle is developed with more generality in \cite{Fageot2017unified} by connecting the theory of generalized random process to independently scattered random measures in the sense of B.S. Rajput and J. Rosinski \cite{Rajput1989spectral}. 
	In particular, as soon as $\varphi \in L_2( \mathbb{R}^d)$ is compactly supported, the random variable $\langle w, \varphi \rangle$ is well-defined. Daubechies wavelets or indicator functions of measurable sets with finite Lebesgue measures are in this case. This was implicitely used in Section \ref{subsec:taurhointro} when considering the random variable $\langle w,  \One_{[0,1]^d}\rangle$.

	\subsection{Indices of L\'evy White Noises} \label{subsec:indices}
	
	First of all, we will only consider L\'evy white noises whose L\'evy exponent satisfies the \emph{sector condition}; that is, there exists $M >0$ such that
	\begin{equation} \label{eq:sectorcondition}
		\abs{\Im\{ \Psi\} (\xi)} \leq M \abs{\Re \{\Psi\}(\xi)}.
	\end{equation}
	This condition ensures that no drift is dominating the L\'evy white noise (a drift appears as purely imaginary in the L\'evy exponent) and is needed for linking the indices with the L\'evy measure~\cite{Schilling1998growth,Deng2015shift}.  \\	
	
	In Theorem \ref{maintheo:regdecaylimit}, the smoothness and growth rate of L\'evy white noises is characterized in terms of the indices \eqref{eq:indass}, \eqref{eq:indasi}, and \eqref{eq:indmom}. We give here some additional insight about these quantities. The two former are classical, while the latter has never been considered to characterize the behavior of L\'evy processes or noises for the best of our knowledge. For the notations of the different indices, with the exception of $p_{\max}$, we follow \cite{Bottcher2014levy}.

The index $\beta_\infty$ was introduced by R. Blumenthal and R. Getoor \cite{Blumenthal1961sample} in order to characterize the behavior of L\'evy process at the origin. This quantity appears to characterize many local properties of random processes driven by L\'evy white noises, including the Hausdorff dimension of the image set \cite{Bottcher2014levy}, the spectrum of singularities \cite{Jaffard1999multifractal,Durand2012multifractal}, the Besov regularity \cite{Schilling1997Feller,Schilling2000function,Fageot2017multidimensional,Bottcher2014levy}, the local self-similarity \cite{Fageot2016scaling}, or the local compressiblity \cite{Fageot2017nterm}. Finally, the   index  $\ubar{\beta}_\infty$   plays a crucial role in the specification of negative results; that is, to identify the Besov spaces in which the L\'evy white noises are not. It satisfies moreover the relation $\ubar{\beta}_\infty \leq \beta_\infty$.

In \cite{Pruitt1981growth}, W. Pruitt proposed the index
\begin{equation}
\beta_0 = \sup \setb{p>0}{\underset{\abs{\xi} \rightarrow 0}{\lim} \frac{\abs{\Psi(\xi)}}{\abs{\xi}^p}= 0}   \label{eq:indpruitsup}
\end{equation}
 as the asymptotic counterpart of $\beta_\infty$. This quantity appears in the asymptotic growth rate of the supremum of L\'evy(-type) processes \cite{Schilling1998growth} and the asymptotic self-similarity of random processes driven by L\'evy white noises \cite{Fageot2016scaling}. The Pruitt index differs from the index $p_{\max}$ appearing in Theorem \ref{maintheo:regdecaylimit}. Actually, the two quantities are linked by the relation $\beta_0 = \inf(p_{\max}, 2)$. This is shown by linking $\beta_0$ to the L\'evy measure~\cite{Bottcher2014levy} and knowing that $\beta_0 \leq 2$ (see the appendix of \cite{Deng2015shift} for a short and elegant proof). This means that $\beta_0 < p_{\max}$ when the L\'evy white noise has finite moments of order bigger than $2$, and one cannot recover $p_{\max}$ from $\beta_0$ in this case. It is therefore required to introduce the index $p_{\max}$  in addition to the Pruitt index in our analysis. \\

The first part of the next result is the famous L\'evy-It\^o decomposition, usually given for L\'evy processes and reformulated here for L\'evy white noises. We add a proof in the context of generalized random processes for the sake of completeness.
The second part presents how does the indices of the corresponding noise behave.

\begin{proposition} \label{prop:decomposition}
	A L\'evy white noise $w$ can be decomposed as
	\begin{equation*}
		w = w_{1} +  w_2 + w_3
	\end{equation*}	
	with  $w_{1}$  a Gaussian noise, $w_2$ a compound Poisson noise, and $w_3$ a L\'evy white noise with finite moments, the three noises   being independent.
	Moreover, we have the following relations between the indices:
	\begin{itemize}
		\item If $w$ has no Gaussian part ($w_1 = 0$), then
		\begin{align*}		
			\ubar{\beta}_\infty  (w) &= \ubar{\beta}_\infty (w_3) \geq \ubar{\beta}_\infty (w_2) 
			 = 0; \\
		\beta_\infty (w) &= \beta_\infty  (w_3) \geq \beta_\infty (w_2) 
		  = 0; \\
			p_{\max} (w) &= p_{\max} (w_2) \leq p_{\max} (w_3) = \infty.
		\end{align*}
			\item If $w$ has a Gaussian part, then
		\begin{align*}		
			\ubar{\beta}_\infty(w) &= \beta_\infty(w) = 2  ; \\
			p_{\max}(w) &= p_{\max}(w_2) \leq p_{\max}(w_1) = p_{\max}(w_3) = \infty.
		\end{align*}
	\end{itemize}
\end{proposition}

\begin{proof}
	Two generalized random processes $s_1$ and $s_2$ are independent if and only if one has, for every $\varphi \in {\mathcal{S}}( \mathbb{R}^d)$,
	\begin{equation}
		{\widehat{\mathscr{P}}}_{s_1 + s_2}(\varphi) = {\widehat{\mathscr{P}}}_{s_1}(\varphi) {\widehat{\mathscr{P}}}_{s_2}(\varphi).
	\end{equation}
	We split the L\'evy exponent $\Psi$ of $w$ with L\'evy triplet $(\mu,\sigma^2,\nu)$ as
	\begin{align}
		\Psi(\xi) 	&= \left( \mathrm{i} \mu \xi + \frac{\sigma^2\xi^2}{2} \right)  + \left(  \int_{ \mathbb{R}} (\mathrm{e}^{\mathrm{i} \xi t} - 1 )  \One_{\abs{t} > 1} \mathrm{d} \nu (t)  \right)+ \left(  \int_{ \mathbb{R}} (\mathrm{e}^{\mathrm{i} \xi t} - 1 - \mathrm{i} \xi t  ) \One_{\abs{t} \leq 1} \mathrm{d} \nu (t)
 \right) \\
				&:= \Psi_1 (\xi) + \Psi_2(\xi) + \Psi_3(\xi).
	\end{align}
	Then, $\Psi_{1}$, $\Psi_{2}$, and $\Psi_3(\xi)$ are three L\'evy exponents with respective triplets $(\mu, \sigma^2, 0)$, $(0,0, \One_{\abs{\cdot} \leq 1} \nu)$, and $(0,0, \One_{\abs{\cdot} > 1} \nu)$. We denote by $w_1 $, $w_2$, and $w_3$ the respective underlying L\'evy white noises. We easily see that
	\begin{equation}
		{\widehat{\mathscr{P}}}_w(\varphi) = {\widehat{\mathscr{P}}}_{w_1}(\varphi) {\widehat{\mathscr{P}}}_{w_2}(\varphi) {\widehat{\mathscr{P}}}_{w_3}(\varphi)
	\end{equation}
	showing the independence of the noises. It is clear that $w_1$ is Gaussian. Moreover, $\int_{ \mathbb{R}} \mathrm{d} \nu_2(t) = \int_{\abs{t} > 1} \mathrm{d} \nu (t) < \infty$, hence $w_2$ is a compound Poisson noise in the sense of Definition \ref{def:levynoise}. Finally, the moments of $w_3$ are finite since $\int_{\abs{t}>1} \abs{t}^p \mathrm{d} \nu_3 (t) = 0 < \infty$ for every $p>0$ and using \cite[Theorem 25.3]{Sato1994levy}. 
	
	The relations between the Blumenthal-Getoor indices are easily deduced from the fact that $\Psi(\xi) = \Psi_1(\xi) + \Psi_2(\xi) + \Psi_3(\xi)$ and $\Psi_1(\xi) = - \sigma_0^2 \xi^2 / 2$, with $\sigma_0^2$ the variance of $w_1$. For the moment indice, we recall that for any $p>0$, there exists $c_p>0$ such that $\abs{a+b}^p \leq c_p (\abs{a}^p + \abs{b}^p)$. Therefore, if $X$ and $Y$ are two random variables such that $\mathbb{E}[\abs{Y}^p] < \infty$, we have
	\begin{equation} \label{eq:XplusY}
		c_p^{-1} (\mathbb{E}[\abs{X}^p] - \mathbb{E}[\abs{Y}^p])
		\leq \mathbb{E}[\abs{X + Y}^p]
		\leq
		c_p(\mathbb{E}[\abs{X}^p]  + \mathbb{E}[\abs{Y}^p]). 
	\end{equation}
	Applying this to $X  = \langle w_2 , \One_{[0,1]^d} \rangle$ and $Y = \langle w_1 + w_3 , \One_{[0,1]^d} \rangle$, which has finite $p$th moments for any $p>0$, we deduce from \eqref{eq:XplusY} that
	\begin{equation}
	\mathbb{E} \left[ \abs{\langle w ,  \One_{[0,1]^d} \rangle}^p\right] =\mathbb{E}[\abs{X+Y}^p]   < \infty \ \Longleftrightarrow \ \mathbb{E} \left[ \abs{\langle w_2 ,  \One_{[0,1]^d} \rangle}^p\right] = \mathbb{E}[\abs{X}^p]  < \infty.
	\end{equation}
	Hence, $w$ and $w_2$ have the same moment index. 
\end{proof}

	\subsection{Weighted Besov Spaces} \label{subsec:besovspaces}
	
	We define the family of weighted Besov spaces based on wavelet methods, as exposed in \cite{Triebel2008function}. Essentially, Besov spaces are subspaces of ${\mathcal{S}}'( \mathbb{R}^d)$ that are characterized by weighted sequence norms of the wavelet coefficients.
Following Triebel, we use Daubechies wavelets, which we introduce first.

The scale and shift parameters of the wavelets are respectively denoted by $j\geq 0$ and $\bm{k}\in  \mathbb{Z}^d$. 
The letters $M$ and $F$ refer to the \emph{gender} of the wavelet ($F$ for the father wavelets and $G$ for the mother wavelet). 
Consider two  functions $\psi_{M}$ and $\psi_{F} \in L_2( \mathbb{R})$.
We set $\mathrm{G}^0 = \{ M, F \}^{d}$ and, for $j\geq 1$, $\mathrm{G}^j = \mathrm{G}^0 \backslash \{ F^d \}$. For $G = (G_1,\ldots, G_d) \in \mathrm{G}^0$, called a gender, we set, for every $\bm{x} = (x_1, \ldots , x_d) \in  \mathbb{R}^d$, $\psi_G(\bm{x} ) = \prod_{i=1}^d \psi_{G_i}(x_i)$. 
For $j\geq 0$, $G \in \mathrm{G}^j$, and $\bm{k} \in   \mathbb{Z}^d$, we define 
\begin{equation}
	\psi_{j,G,\bm{k}} := 2^{jd/2} \psi_G(2^j \cdot - \bm{k}).
\end{equation}

It is known that, for any regularity parameter $r_0 \geq 1$, there exists two functions $\psi_M, \psi_F \in L_2( \mathbb{R})$ that are compactly supported, with at least $r_0$ continuous derivatives such that the family 
	\begin{equation}
	\{ \psi_{j,G,\bm{k}} \}_{(j,G,\bm{k})\in  \mathbb{N} \times \mathrm{G}^j \times  \mathbb{Z}^d}
	\end{equation}
is an orthonormal basis of $L_2( \mathbb{R}^d)$ \cite{Triebel2008function}. Concretely, one can consider the family of {Daubechies wavelets} \cite{Daubechies1988orthonormal}.

We now introduce the family of weighted Besov spaces $B_{p}^{\tau}( \mathbb{R}^d;\rho)$. Traditionally, Besov spaces also depends on the additional parameter $q \in (0,\infty]$ (see for instance \cite[Definition X]{Triebel2008function}).  We should only consider the case $q=p$ in this paper, so that we do not refer to this parameter.

The following definition of weighted Besov spaces is based on wavelets. It is equivalent with the more usual Fourier-based definitions. This equivalence is proved in \cite{Triebel2008function}.  

\begin{definition}\label{def:besovspaces}
	Let $\tau, \rho \in  \mathbb{R}$ and $0< p \leq \infty$. Fix $r_0 > \max ( \tau , d(1 / p - 1)_+ - \tau )$ and consider a family of compactly supported wavelets $\{ \psi_{j,G,\bm{k}} \}_{(j,G,\bm{k})\in  \mathbb{N} \times \mathrm{G}^j \times  \mathbb{Z}^d}$ with at least $r_0$ continuous derivatives.

	The \emph{weighted Besov space} $B_p^{\tau}( \mathbb{R}^d;\rho)$ is the collection of tempered generalized functions $ f \in {\mathcal{S}}'( \mathbb{R}^d)$ that can be written as
	\begin{equation} \label{eq:fconvuncond}
		f = \sum_{(j,G,\bm{k}) \in  \mathbb{N}\times \mathrm{G}^j \times  \mathbb{Z}^d} c_{j,G,\bm{k}} \psi_{j,G,\bm{k}}.
	\end{equation}
	where the convergence holds unconditionally on ${\mathcal{S}}'( \mathbb{R}^d)$. 
\end{definition}

The parameter $r_0$ in Definition \ref{def:besovspaces} is chosen such that the wavelet is regular enough to be applied to a function of $B_p^{\tau}( \mathbb{R}^d;\rho)$.
When the convergence \eqref{eq:fconvuncond} occurs, the duality product $\langle f, \psi_{j,G,\bm{k}} \rangle$ is well defined and we have $c_{j,G,\bm{k}} = \langle f, \psi_{j,G,\bm{k}} \rangle$. Moreover, the quantity
\begin{equation}
	\lVert f \rVert_{B_p^\tau( \mathbb{R}^d;\rho)} := \left(  \sum_{j\geq 0} 2^{j(\tau p -d + \frac{dp}{2})} 
	 	\sum_{G \in \mathrm{G}^j}
		\sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^{-j} \bm{k}  \rangle^{\rho p}  \lvert \langle f , \psi_{j,G, \bm{k}} \rangle \rvert^p \right)^{1/p}	
\end{equation}
is finite for $f \in B_p^\tau( \mathbb{R}^d;\rho)$ and specifies a norm (a quasi-norm, respectively) on $ B_p^\tau( \mathbb{R}^d;\rho)$ for $p\geq 1$ ($p<1$, respectively). The space $ B_p^\tau( \mathbb{R}^d;\rho)$ is a Banach (a quasi-Banach, respectively) for this norm (quasi-norm, respectively). 

\begin{proposition}[Embeddings between weighted Besov spaces] \label{prop:besovembed}
	Let $0 < p_0 \leq p_1 \leq \infty$ and $\tau_0 , \tau_1 , \rho_0 , \rho_1 \in  \mathbb{R}$.
	\begin{itemize}
		\item We have the embedding $B_{p_0}^{\tau_0}( \mathbb{R}^d;\rho_0) \subseteq B_{p_1}^{\tau_1}( \mathbb{R}^d;\rho_1)$ as soon as
		\begin{equation} \label{eq:embedd1}
			\tau_0 - \tau_1 > \frac{d}{p_0} - \frac{d}{p_1}  \ \text{ and } \ \rho_0 > \rho_1.
		\end{equation}
		\item We have the embedding $B_{p_1}^{\tau_1}( \mathbb{R}^d;\rho_1) \subseteq B_{p_0}^{\tau_0}( \mathbb{R}^d;\rho_0)$ as soon as
		\begin{equation} \label{eq:embedd2}
			\rho_0 - \rho_1 > \frac{d}{p_0} - \frac{d}{p_1}  \ \text{ and } \ \tau_0 > \tau_1.
		\end{equation}
	\end{itemize}	
\end{proposition}

	The embedding for the relation \eqref{eq:embedd1} is proved in \cite[Section 4.2.3]{Edmunds2008function} and the one for \eqref{eq:embedd2} in \cite[Section 2.2.2]{Fageot2017multidimensional}. 
	Proposition \ref{prop:besovembed} is summarized in the two Triebel diagrams\footnote{The representation of the smoothness properties in diagrams with axis $(1/p,\tau)$ is inherited from the work of H. Triebel. It is very convenient, since the smoothness has often a simple formulation in terms of $1/p$. This is also valid for the asymptotic rate $\rho$.} of Figure \ref{fig:embeddings}. 
	
\begin{figure}[t!] 
\centering
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.34]


\fill[green, opacity= 0.3] (0,-1) -- (2,1) -- (3,1) -- (3,-3/2) -- (0,-3/2) -- cycle; 
\draw[green, thick, dashed,->](2,1) --(3,1);
\draw[green, thick, dashed,- ](0,-1) -- (2,1);
\draw[green, thick,->](0,-1) -- (0,-3/2);

\fill[red, opacity= 0.3] (0,1) -- (2,1) -- (3,2) -- (0,2) -- cycle; 
\draw[red, thick, dashed,->](2,1) --(3,2);
\draw[red, thick, dashed,- ](2,1) -- (0,1);
\draw[red, thick,->](0,1) -- (0,2);

\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-3/2)--(0,2) node[circle,above] {$\tau$} ;

\draw[ thick,color=black]
(-0.05,1) -- (0.05,1)  node[black,left] { $\tau_0$};

\draw[ thick,color=black]
(2,-0.05) -- (2,0.05)  node[black,below] { $\frac{1}{p_0}$};



\end{tikzpicture}
\caption{$(1/p,\tau)$-diagram for fixed $\rho_0$}
\end{subfigure}
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.34]

\fill[green, opacity= 0.3] (0,1) -- (1,1) -- (3,-1) -- (3,-3/2) -- (0,-3/2) -- cycle; 
\draw[green, thick, dashed,->](1,1) --(3,-1);
\draw[green, thick, dashed,- ](0,1) -- (1,1);
\draw[green, thick,->](0,1) -- (0,-3/2);

\fill[red, opacity= 0.3] (0,2) -- (1,1) -- (3,1) -- (3,2) -- cycle; 
\draw[red, thick, dashed,->](1,1) --(3,1);
\draw[red, thick, dashed,- ](0,2) -- (1,1);

\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-3/2)--(0,2) node[circle,above] {$\rho$} ;

\draw[ thick,color=black]
(-0.05,1) -- (0.05,1)  node[black,left] { $\rho_0$};

\draw[ thick,color=black]
(1,-0.05) -- (1,0.05)  node[black,below] { $\frac{1}{p_0}$};

\end{tikzpicture}
\caption{$(1/p,\rho)$-diagram for fixed $\tau_0$}
\end{subfigure}

\caption{Representation of the embeddings between Besov spaces: If $f \in B_{p_0}^{\tau_0}( \mathbb{R}^d;\rho_0)$, then $f$ is in any Besov space that is in the shaded green regions. Conversely, if $f \notin B_{p_0}^{\tau_0}( \mathbb{R}^d;\rho_0)$, then $f$ is not in any of the Besov spaces of the shaded red regions.}
\label{fig:embeddings}
\end{figure}
	

	As a simple example, we obtain the Besov localization of the Dirac distribution. Of course, this result is known, and an alternative proof can be found for instance in \cite{schmeisser87}. We  believe that it is interesting to give our own proof here. First, it illustrates how to use the wavelet-based characterization of Besov spaces, and second, it will be used to obtain sharp results for compound Poisson processes, justifying to include the proof for the sake of completeness.
	
	\begin{proposition} \label{prop:DiracBesov}
 	The Dirac impulse $\delta$ is in $B_p^\tau( \mathbb{R}^d;\rho)$ if and only if $\tau < \frac{d}{p} - d$.
	\end{proposition}
	
	\begin{proof}
		The definition of the Besov (quasi-)norm easily gives
		\begin{equation}
			\lVert \delta \rVert^p_{B_p^\tau( \mathbb{R}^d;\rho)} 
			= 
			\sum_{j\geq 0} 2^{j(\tau p - d + dp)} \sum_{G \in \mathrm{G}^j} \sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^{-j} \bm{k} \rangle^{\rho p} \lvert \psi_{G}(k) \rvert^p.
		\end{equation}
		The common support $K_{\psi}$ of the $\psi_G$ is compact. Therefore, only finitely many $\psi_G(\bm{k})$ are non zero, and for such $\bm{k}$ and every $j$ we have
		\begin{equation}
			0 < \min_{\bm{x} \in K_\psi} \langle \bm{x} \rangle^{\rho p} =   \langle 2^{-j} \bm{k} \rangle^{\rho p}  \leq  \max_{\bm{x} \in K_\psi} \langle \bm{x} \rangle^{\rho p}  < \infty.
		\end{equation}	
		It is then easy to find $0<A\leq B < \infty$ such that
		\begin{equation}
			A \sum_{j\geq 0} 2^{j (\tau - d + dp)} \leq \lVert \delta \rVert^p_{B_p^\tau( \mathbb{R}^d;\rho)} \leq B  \sum_{j\geq 0} 2^{j (\tau - d + dp)}.
		\end{equation}
		The sum converges for $\tau - d + dp < 0$ and diverges otherwise, implying the result.
	\end{proof}
	
\section{Gaussian White Noise} \label{sec:Gaussian}

Our goal in this section is to prove the Gaussian part of Theorem \ref{maintheo:regdecaylimit}.
Without loss of generality, we will always assume that the noise has variance $1$. 

The Gaussian case is much simpler than the general one since the wavelet coefficients of the Gaussian noise are independent and identically distributed.
We present it separately for three reasons: (i) it can be considered as an instructive toy problem, containing already some of the technicalities that will appear for the general case, (ii) it cannot be deduced from the other sections, where the results are based on a careful study of the L\'evy measure of the noise (the L\'evy measure of the Gaussian noise is zero), and (iii) we are not aware that the localization of the Gaussian noise in \emph{weighted} Besov spaces has been addressed in the literature (for the local Besov regularity, an in-depth  answer is given in \cite{Veraar2010regularity}). 

We first state a simple lemma that will be useful throughout the paper to show negative results ($w$ \emph{not in} a given weighted Besov space) for the asymptotic rate $\rho$.

\begin{lemma} \label{lemma:sumiid}
	Assume that $(X_{\bm{k}})_{\bm{k} \in  \mathbb{Z}^d}$, is a sequence of i.i.d. nonzero random variables. Then,
	\begin{equation}
		\sum_{\bm{k} \in  \mathbb{Z}^d} \frac{\lvert X_{\bm{k}} \rvert}{\langle \bm{k}\rangle^d} = +\infty \text{ a.s.}
	\end{equation}
\end{lemma}

\begin{proof}
	First of all, the result for any dimension $d$ is easily deduced from the one-dimensional case. Moreover, $\lvert k \rvert$ and $\langle k \rangle$ are equivalent asymptotically, so that it is equivalent to show that $\sum_{k \geq 1} \frac{\lvert X_k \rvert}{k} = \infty$ for $X_k$ i.i.d. For $k \geq 1$, we set $Z_k = \frac{1}{2^{k}} \sum_{l =2^{k-1}}^{2^{k}-1} \lvert X_l \rvert$, so that
	\begin{equation}
	\sum_{k \geq 1} \frac{\lvert X_k \rvert}{k}  = \sum_{k\geq 1} \sum_{l =2^{k-1}}^{2^{k}-1} \frac{\lvert X_l \rvert}{l}  \geq \sum_{k\geq 1} \frac{1}{2^{k}} \sum_{l =2^{k-1}}^{2^{k}-1} \lvert X_l \rvert = \sum_{k\geq 1} Z_k
	\end{equation} 
The $Z_k$ are independent because the $X_k$ are. Moreover, for all $k$, we have $\mathbb{E}[Z_k] = \mathbb{E}[\lvert X_1\rvert]$. The weak law of large numbers ensures that $\mathscr{P}(Z_k >  \mathbb{E}[\lvert X_1\rvert] / 2)$ goes to $1$, therefore $\sum_{k\geq 1 } \mathscr{P}(Z_k >  \mathbb{E}[\lvert X_1\rvert] / 2) = \infty$. Since the events $\{ Z_k >  \mathbb{E}[\lvert X_1\rvert] / 2 \}$ are independent, we can apply the Borel-Cantelli lemma to deduce that infinitely many $Z_k$ are bigger than $\mathbb{E}[\lvert X_1\rvert] / 2$ almost surely. Finally, this implies that $\sum_{k\geq 1}Z_k = \infty$ almost surely and the result is proved. 
\end{proof}


\begin{proposition} \label{prop:GaussBesov}
	The Gaussian white noise $w$ is 
	\begin{itemize}
	\item almost surely in $B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau < - d/2$ and $\rho < - d /p$, and 
	\item almost surely not in $B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau \geq  - d/2$ or $\rho \geq - d /p$.
	\end{itemize}
\end{proposition}


\begin{proof
	\textbf{If $\tau < - d/2$ and $\rho < - d /p$.}
	For $p > 0$, we denote by $C_p$ the $p$-th moment of a Gaussian random variable with $0$ mean and variance $1$.
	For the Gaussian noise, $\langle w, \varphi_1 \rangle$ and $\langle w, \varphi_2\rangle$ are independent if and only if $\langle \varphi_1 ,\varphi_2 \rangle = 0$, and $\langle w, \varphi \rangle$ is a Gaussian random variable with variance $\lVert \varphi \rVert_2^2$ \cite{Unser2014sparse}. The family of functions $( \psi_{j,G,\bm{k}} )_{j,G,\bm{k}}$ being orthonormal, the random variables $\langle w , \psi_{j,G,\bm{k}} \rangle$  are therefore i.i.d. with law $\mathcal{N}(0,1)$. 
	We then have
	\begin{align}
	\mathbb{E}[\|w  \|_{B_p^\tau(\mathbb{R}^d;\rho)}^p] 
	&= 
	 \sum_{j\geq 0} 2^{j(\tau p -d + \frac{dp}{2})} 
	 	\sum_{G \in \mathrm{G}^j}
		\sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^{-j} \bm{k}  \rangle^{\rho p}  \mathbb{E} [\lvert \langle w , \psi_{j,G, \bm{k}} \rangle \rvert^p ]	 
 	\nonumber   \\ 	
	& = 
	C_p \sum_{j\geq 0} 2^{j(\tau p -d + \frac{dp}{2})}  \mathrm{Card}(\mathrm{G}^j) \sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^{-j} \bm{k}  \rangle^{\rho p} \nonumber \\
	& \leq 2^d C_p \sum_{j\geq 0} 2^{j(\tau p -d + \frac{dp}{2})}  \sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^{-j} \bm{k}  \rangle^{\rho p}.
	\end{align}
	The last inequality is due to $\mathrm{Card} ( \mathrm{G}^j ) \leq 2^d$. 
	Since $\rho p < -d$ and 
	$\langle 2^{-j} \bm{k} \rangle \underset{\lVert \bm{k}\rVert \rightarrow \infty}{\sim} 2^{-j} \lVert \bm{k} \rVert $, 
	we have that 
	$\sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^j \bm{k} \rangle^{\rho p} < \infty$. 
	Moreover, we recognize a Riemann sum and have the convergence
	\begin{equation}
		2^{-jd} \sum_{\bm{k}\in  \mathbb{Z}^d} \langle 2^{-j} \bm{k} \rangle^{\rho p}  \underset{j \rightarrow \infty}{\longrightarrow} \int_{ \mathbb{R}^d} \langle \bm{x} \rangle^{\rho p} \mathrm{d} \bm{x} < \infty.
	\end{equation}
	In particular, the series $\sum_j 2^{j(\tau p  + \frac{dp}{2})}  \left( 2^{-jd} \sum_{\bm{k}} \langle 2^{-j} \bm{k}  \rangle^{\rho p}\right)$ converges if and only if the series 
	$\sum_j 2^{j(\tau p   + \frac{dp}{2})}$ does; that is, if and only if $\tau < d / 2$. Finally, if $\tau < d / 2$ and $\rho < - d / p$, we have  shown that $\mathbb{E}[\|w  \|_{B_p^\tau(\mathbb{R}^d;\rho)}^p]  < \infty$ and therefore $w \in {B_p^\tau(\mathbb{R}^d;\rho)}$ almost surely. \\
	
	\textbf{If $\tau \geq - d/2 $.}
	Then, we have $2^{j (\tau - d + dp/2)} \geq 2^{-jd}$. 
	We aim at lower bound the Besov norm of $w$ and we restrict to the wavelet with gender $G = M^d \in \mathrm{G}^j$ for any $j \geq 0$. 
	Hence, we have, since $\langle 2^{-j} \bm{k} \rangle \geq 1$, 
	\begin{align}
		\lVert w \rVert^p_{B_p^{\tau}( \mathbb{R}^d;\rho)} \geq   \sum_{j \geq 0} 2^{-jd} \sum_{0\leq k_1, \ldots , k_d < 2^{j}} \lvert \langle w , \psi_{j, M^d ,\bm{k}} \rangle \rvert^p :=   \sum_{j\geq 0} Z_j.
	\end{align}
	The random variables $Z_j = 2^{-jd} \sum_{0\leq k_1, \ldots , k_d < 2^{j}} \lvert \langle w , \psi_{j, M^d ,\bm{k}} \rangle \rvert^p$ are independent, non-negative, and have the same average $\mathbb{E}[Z_j]=C_p$ equals to the $p$th-moment of a Gaussian random variable with variance $1$. 
	By the weak law of large numbers, we know that $\mathscr{P} ( Z_j \geq C_p/2)$ converges to $1$, 
	and then $\sum_{j\geq0} \mathscr{P} (Z_j \geq C_p/2) = \infty$. 
	The events $\{Z_j \geq C_p /2\}$ being independent, the Borel-Cantelli lemma ensures that $Z_j \geq C_p /2$ for infinitely many $j$ almost surely. 
	This ensures that $\lVert w \rVert^p_{B_p^{\tau}( \mathbb{R}^d;\rho)} \geq    \sum_{j\geq 0} Z_j = \infty$ almost surely. \\
	
	\textbf{If $\rho \geq - d/p $.}
	Keeping only the father wavelet $\phi = \psi_{F^d}$ and the scale $j = 0$ and exploiting the relation $\rho p \geq -d$, we have the lower bound 
	\begin{equation}
	\lVert w \rVert^p_{B_p^{\tau}( \mathbb{R}^d;\rho)} 
	\geq 
	\sum_{\bm{k}\in  \mathbb{Z}^d} \langle \bm{k} \rangle^{\rho p} \lvert \langle w, \phi_{\bm{k}} \rangle \rvert^p 
	\geq
	 \sum_{\bm{k} \in  \mathbb{Z}^d} \frac{\lvert\langle w, \phi_{\bm{k}} \rangle \rvert^p}{\langle \bm{k} \rangle^d }
	\end{equation}
	Finally, the random variables $\langle w, \phi_{\bm{k}} \rangle$ being i.i.d., Lemma \ref{lemma:sumiid} applies and $\lVert w \rVert^p_{B_p^{\tau}( \mathbb{R}^d;\rho)}  = \infty$ almost surely. 	
\end{proof}


\section{Compound Poisson Noises} \label{sec:Poisson}
	
	\subsection{Moment estimates for compound Poisson noises}
	
	Our positive results for compound Poisson noises are based on a careful estimation of moments for compound Poisson noise.
	
	\begin{lemma} \label{lemma:momentPoisson}
		Let $w$ be a compound Poisson process with moment index $p_{\max}$ and $p < p_{\max}$.
		Then, there exists a constant $C$ such that
		\begin{equation}
			\mathbb{E} [ \abs{\langle w, \psi_{j,G, \bm{k}} \rangle}^p ] \leq C 2^{jpd/2 - jd}
		\end{equation}
		for every $j\geq 0$, $G\in \mathrm{G}^j$, and $\bm{k}\in  \mathbb{Z}^d$.
	\end{lemma}

	\begin{proof}
	We denote by $\lambda >0$ the sparsity parameter of $w$. Being a compound Poisson noise, $w$ is equal in law to
	\begin{equation}
		w \overset{(\mathcal{L})}{=} \sum_{k\in  \mathbb{N}} a_k \delta (\cdot - \bm{x}_k )
	\end{equation}
	according to \cite{Unser2011stochastic}. For $M$ big enough, the support of the $\Psi_{G}$ is included in $[-M/2,M/2]^d$. Hence, the support of $\Psi_{j,G,\bm{k}}$ is included in 
	\begin{equation} \label{eq:Ijbmk}
	I_{j,\bm{k}} := \prod_{i=1}^d [2^{-j}(k_i - M/2) , 2^{-j}(k_i + M/2) ].
	\end{equation}
	We set
	\begin{equation} 
		A_{j,\bm{k}} = \mathrm{Card} \{ n \in  \mathbb{N}, \ \bm{x}_n \in I_{j,\bm{k}}\}.
	\end{equation}
	It is a Poisson random variable with parameter $\lambda L(I_{j,\bm{k}}) = \frac{\lambda M^d}{2^{jd}}$. 
	Then, we have the equality in law
	\begin{equation}
		\langle w , \psi_{j,\bm{k}} \rangle \overset{(\mathcal{L})}{=} \sum_{n=1}^{A_{j,\bm{k}}} a'_n \psi_{j,\bm{k}}  (\bm{x}'_n )
	\end{equation}
	where the $a_k'$ are i.i.d. with the same law than the $a_k$. The law of the $x'_k$ can be specified explicitly but will play no role in the sequel. 	
	By conditioning on  $A_{j,\bm{k}}$ and using the relation $ \abs{\sum_{n=1}^N x_n}^p \leq N^{\max( 0, p-1)} \sum_{n=1}^N \abs{x_n}^p$, 
	we deduce that
	\begin{align} \label{eq:tructructruc}
		\mathbb{E} [ \abs{\langle w, \psi_{j,G, \bm{k}} \rangle}^p ]  & = 
			\sum_{N=1}^{\infty} \mathscr{P} ( A_{j,\bm{k}} = N ) \mathbb{E} \left[ \abs{\langle w, \psi_{j,G, \bm{k}} \rangle}^p  | A_{j,\bm{k}} = N \right] 	 \nonumber \\
			&=   
			\sum_{N=1}^{\infty} \mathscr{P} ( A_{j,\bm{k}} = N )
			 \mathbb{E} \left[ \abs{ \sum_{n=1}^{N} a'_n \psi_{j, G, \bm{k}}  (\bm{x}'_n ) }^p \right] \nonumber \\
			&\leq   
			\sum_{N=1}^{\infty} \mathscr{P} ( A_{j,\bm{k}} = N ) N^{\max(0,p-1)} 
			 \mathbb{E}  \left[	\sum_{n=1}^{N} \abs{ a'_n \psi_{j, G, \bm{k}}  (\bm{x}'_nn ) }^p\right] \nonumber \\
			&\leq   
			\lVert \psi_{j, G, \bm{k}} \rVert_\infty^p \sum_{N=1}^{\infty} \mathscr{P} ( A_{j,\bm{k}} = N ) N^{\max(1,p)} 			  \mathbb{E}  \left[ \abs{ a_1 }^p\right] \nonumber \\
			&=  2^{jdp/2} \lVert \psi_{G} \rVert_\infty^p   \mathbb{E} \left[ \abs{ a_1 }^p\right]   \sum_{N=1}^{\infty} N^{\max(1,p)}  \mathscr{P} ( A_{j,\bm{k}} = N ) ,		 
	\end{align}
	where we used at the end the obvious relation  $\lVert \psi_{j, G, \bm{k}} \rVert_\infty^p = 2^{jdp/2} \lVert \psi_{G} \rVert_\infty^p$. 
	Knowing the law of $A_{j,\bm{k}}$, we then have
	\begin{equation}
		  \sum_{N=1}^{\infty}  N^{\max(1,p)}  \mathscr{P} ( A_{j,\bm{k}} = N ) 
		  	 =   \sum_{N=1}^{\infty} N^{\max(1,p)}  \frac{1}{N!} (M^d \lambda)^N 2^{-jd N} \mathrm{e}^{- \lambda M^d 2^{-jd}}.
	\end{equation}
	Then, $2^{-jdN} \leq 2^{-jd}$ for every $N\geq 1$ and $\mathrm{e}^{- \lambda M^d 2^{-jd}} \leq 1$, hence,
	\begin{equation}  \label{eq:tructruc}
	\sum_{N=1}^{\infty}  N^{\max(1,p)}  \mathscr{P} ( A_{j,\bm{k}} = N ) \leq \tilde{C} 2^{-jd} 
	\end{equation}
	where $\tilde{C} =  \sum_{N=1}^{\infty} N^{\max(1,p)}  \frac{1}{N!} (M^d \lambda)^N < \infty$. Finally, including \eqref{eq:tructruc} into  \eqref{eq:tructructruc}, we deduce the result with $C= \tilde{C} \mathbb{E} [\abs{a_1}^p] \lVert \psi_{G} \rVert_\infty^p$.
	\end{proof}
	
	\subsection{The Besov regularity of compound Poisson noises}
	
	We shall prove the following result.
	
	\begin{proposition} \label{prop:PoissonBesov}
	Let $w$ be a compound Poisson noise with index $p_{\max} \in (0,\infty]$.
		Then, $w$ is 
		\begin{itemize}
		\item almost surely in $B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau < d/p - d$ and $\rho < - d /\min(p,p_{\max})$, and 
		\item almost surely not in $B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau \geq d/p - d$ or $\rho > - d /\min(p,p_{\max})$.	
		\end{itemize}
	\end{proposition}

\begin{proof}
	\textbf{If $p < p_{\max}$, $\tau < d/p - d$, and $\rho < - d / p$.}
	Under this assumption, we apply Lemma \ref{lemma:momentPoisson} to deduce that
	\begin{align} \label{eq:intermposCPreg}
		\mathbb{E} [ \lVert w \rVert_{B_p^\tau( \mathbb{R}^d;\rho)} ] &= \sum_{j\geq 0} 2^{j (\tau p - d + dp/2)} \sum_{G,\bm{k}} \langle 2^{-j} \bm{k} \rangle^{\rho p} \mathbb{E} [ \abs{\langle w, \psi_{j,G, \bm{k}} \rangle}^p ] \nonumber \\
		& \leq 
		C 2^d \sum_{j\geq 0} 2^{j (\tau p - d + dp)} \frac{1}{2^{jd}} \sum_{\bm{k} \in  \mathbb{Z}^d}  \langle 2^{-j} \bm{k} \rangle^{\rho p}.
	\end{align}
	The sum over the gender was removed using that $\mathrm{Card}(\mathrm{G}^j) \leq 2^d$. Then, $\frac{1}{2^{jd}} \sum_{\bm{k} \in  \mathbb{Z}^d}  \langle 2^{-j} \bm{k} \rangle^{\rho p} \underset{j \rightarrow \infty}{\longrightarrow} \int_{ \mathbb{R}^d} \langle \bm{x} \rangle^{\rho p}\mathrm{d} \bm{x} < \infty$ as soon as $\rho < - d / p$. Assuming this condition on $\rho$, the sum in \label{eq:intermposCPreg} is finite if and only if $\sum_{j} 2^{j (\tau p - d + dp)} < \infty$, that is, if and only if $\tau p - d + d/p<0$, as expected. \\
	
\textbf{If $p \geq p_{\max}$, $\tau < d/p - d$, and $\rho < - d / p_{\max} $.}
	From the conditions on $\tau$ and $\rho$, one can find $p_0$, $\rho_0$, and $\tau_0$ such that
	\begin{align}
		\rho <   \rho_0  < - \frac{d}{p_0} < - \frac{d}{p_{\max}}, \\
		\tau + \frac{d}{p_0} - \frac{d}{p} <   \tau_0 < \frac{d}{p_0} - d.
	\end{align}
	Then, in particular, $p_0 < p$, $\tau_0 - \tau > d / p_0 - d / p$, and $\rho_0 > \rho$,
	so that $B_{p_0}^{\tau_0}( \mathbb{R}^d;\rho_0) \subset B_p^\tau( \mathbb{R}^d;\rho)$ (according to \eqref{eq:embedd1}).
	Moreover, $p_0 < p_{\max}$, $\tau_0 < d / p_0 -d$, and $\rho_0 < - d / p_0$. We are therefore back to the first case, for which have shown that $w \in B_{p_0}^{\tau_0}( \mathbb{R}^d;\rho_0)$ almost surely.
	In conclusion, $w \in  B_p^\tau( \mathbb{R}^d;\rho) $ almost surely. \\
	
	Combining this first two cases, we obtain that $w \in  B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau < d/p - d$ and $\rho < - d /\min(p,p_{\max})$. \\


	\textbf{If $\tau \geq d/p - d$.}
The compound Poisson noise $w$ is equal in law to
	 \begin{equation} \label{eq:formPoisson}
	 	\sum_{n \in  \mathbb{Z}} a_n \delta (\cdot - \bm{x}_n)
	 \end{equation}
	 with $(a_n)$ i.i.d. with law of jumps $\mathscr{P}_{\mathrm{J}}$ and $(\bm{x}_n)$ the jumps location, independent of $(a_n)$ and such that $\mathrm{Card} \{ n , \bm{x}_n \in B \}$ is a Poisson random variable of parameter $\lambda \mathrm{Leb}(B)$, with $\mathrm{Leb}(B)$ the Lebesgue measure of any Borelian set $B \subset  \mathbb{R}^d$. 
	Assume that $w$ is in   $B_p^\tau( \mathbb{R}^d;\rho)$ for some $\rho \in  \mathbb{R}$. 
	Then, the product of $w$ by any compactly supported smooth test function $\varphi$ is also in $B_p^\tau( \mathbb{R}^d;\rho)$. Choosing $\varphi$ such that $\varphi(\bm{x}_0) = 1$ and $\varphi(\bm{x}_n) = 0$ for $n \neq 0$, we get 
	\begin{equation}
		\varphi \cdot w = a_0 \delta(\cdot - \bm{x}_0) \in B_p^\tau( \mathbb{R}^d;\rho),
	\end{equation}
	what is absurd due to Proposition  \ref{prop:DiracBesov}. This proves that $w \notin  B_p^\tau( \mathbb{R}^d;\rho)$ for any $\rho \in  \mathbb{R}$. \\

	\textbf{If  $\rho \geq - d /p$.}
	By restricting to the scale $j \geq 0$, with only the father wavelet $\phi$, and selecting $k_0$ such that the $\phi_{\bm{k}}$ have disjoint supports two by two for $\bm{k} \in k_0  \mathbb{Z}^d$, we have the lower bound
	\begin{equation} \label{eq:lowerbesovwfather}
		\lVert w \rVert^p_{B_p^\tau( \mathbb{R}^d;\rho)} \geq \sum_{\bm{k} \in k_0  \mathbb{Z}^d} \langle \bm{k} \rangle^{\rho p} \lvert \langle w , \phi_{\bm{k}} \rangle\rvert^p \geq \sum_{\bm{k} \in k_0  \mathbb{Z}^d} \frac{\lvert \langle w , \phi_{\bm{k}} \rangle\rvert^p}{\langle \bm{k}\rangle^d}.
	\end{equation}
	The $\phi_{\bm{k}}$ having disjoint supports, the random variables $\langle w, \phi_{\bm{k}} \rangle$ are i.i.d. when $\bm{k} \in k_0  \mathbb{Z}^d$. Lemma \ref{lemma:sumiid} hence applies and $\lVert w \rVert^p_{B_p^\tau( \mathbb{R}^d;\rho)} = \infty$, as expected.  
	
	\textbf{If  $p > p_{\max}$ and $\rho > - d / p_{\max}$.}
	This means in particular that $p_{\max} < \infty$. 
	We treat the case $\rho < 0$, the extension for $\rho \geq 0$ coming easily by embedding.
	We set $q = - d / \rho$.
	Proceeding as for \eqref{eq:lowerbesovwfather}, we have that
\begin{equation} \label{eq:lowerbesovwfather2}
		\lVert w \rVert^p_{B_p^\tau( \mathbb{R}^d;\rho)}  
		\geq 
		\sum_{\bm{k} \in k_0  \mathbb{Z}^d} \frac{\lvert \langle w , \phi_{\bm{k}} \rangle\rvert^p}{\langle \bm{k}\rangle^{dp/q}}.
	\end{equation}	
	Consider the events $A_{\bm{k}} = \{ \abs{\langle w, \phi_{\bm{k}} \rangle} \geq \langle \bm{k} \rangle^{d/q} \}$	for $\bm{k} \in k_0  \mathbb{Z}^d$. The $A_{\bm{k}}$ are independent because the $X_{\bm{k}} = \langle w, \phi_{\bm{k}} \rangle$ are. Moreover, the $X_{\bm{k}}$ have the same law, since $w$ is stationary. Set $Y = \abs{X_{\bm{0}}}^{q}$. Then, 
	\begin{equation} \label{eq:sumAk}
		\sum_{\bm{k} \in k_0 \mathbb{Z}^d} \mathscr{P}( A_{\bm{k}} ) = \sum_{\bm{k} \in k_0 \mathbb{Z}^d} \mathscr{P}( Y \geq \langle \bm{k} \rangle^d ) \geq \sum_{m \geq 1}\mathscr{P} (Y \geq m k_0).
	\end{equation}
	Moreover, we have, exploiting that $\mathscr{P} ( Y \geq x )$ is decreasing in $x$, 
	\begin{equation} \label{eq:momentY}
	\mathbb{E} [Y]  = \int_0^\infty \mathscr{P}(Y \geq x) \mathrm{d} x  
	 = \sum_{m \geq 1} \int_{m k_0}^{(m+1) k_0} \mathscr{P}( Y \geq x) \mathrm{d} x  	 \leq \sum_{m \geq 1}  \mathscr{P}( Y \geq mk_0 )   		\end{equation}
	The relation $q = \frac{d}{-\rho} > p_{\max}$ implies that $\mathbb{E} [Y] = \mathbb{E} [\abs{X_{\bm{0}}}^{q}] =  \infty$.
	Hence, from \eqref{eq:sumAk} and \eqref{eq:momentY}, we deduce that  $	\sum_{\bm{k} \in k_0 \mathbb{Z}^d} \mathscr{P}( A_{\bm{k}} )  = \infty$. The Borel-Cantelli lemma implies that $ \abs{X_{\bm{k}} }^p \geq \langle \bm{k} \rangle^{ p d / q}$ for infinitely many $\bm{k}$ almost surely. Due to \eqref{eq:lowerbesovwfather2}, this implies that $\lVert w \rVert^p_{B_p^\tau( \mathbb{R}^d;\rho)}  = \infty$ almost surely and the result is proved.
\end{proof}


\section{L\'evy White Noises with Finite Moments} \label{sec:finitemoments}

	\subsection{Moment estimates for L\'evy white noise}
	
As a preliminary result, we estimate the evolution of the even moments of the wavelet coefficients of $w$ with the scale $j$. This result is adapted to the case where the noise has finite moments of order bigger than $2$. Many moment estimates in the litterature usually deal with $p$th-moments where $p\leq 2$~\cite{Luschgy2008moment,Deng2015shift,Kuhn2017existence,Fageot2017multidimensional}. 

	\begin{lemma} \label{lemma:boundcoefficientfinitemoments}
		Let $w$ be a L\'evy white noise with finite moments. Then, for every integer $m\geq 1$, there exists a constant $C>0$ such that, for every $j\geq 0$, $G \in \mathrm{G}^j$, and $\bm{k}\in  \mathbb{Z}^d$,
		\begin{equation} \label{eq:boundcoefficientsfinitemoment}
			\mathbb{E} [\langle w, \psi_{j,G,\bm{k}} \rangle^{2m}] \leq C 2^{j d (m-1)}.
		\end{equation}
	\end{lemma}
	
	\begin{proof}
		Consider a test function $\varphi \in {\mathcal{S}}( \mathbb{R}^d)$ and set $X = \langle w , \varphi \rangle$.
		The characteristic function of $X$ is 
		\begin{equation} \label{eq:CFX}
			{\widehat{\mathscr{P}}}_X(\xi) = \exp\left( \int_{ \mathbb{R}^d} \Psi(\xi  \varphi (\bm{x}) )\mathrm{d} \bm{x} \right) := \exp( \Psi_{\varphi}(\xi)).
		\end{equation}
		The functions ${\widehat{\mathscr{P}}}_X$ and $\Psi_{\varphi}$ are infinitely differentiable because all the moments of $X$ are finite. 
		Their Taylor expansions give respectively the moments and the cumulants of $X$. In particular, we have $\mathbb{E}[X^{2m}] = (-1)^m {\widehat{\mathscr{P}}}_X^{(2m)}(0)$. 
		The $(2m)$th derivative of ${\widehat{\mathscr{P}}}_X$ is deduced for instance from the Fa\`a di Bruno formula, and is
		\begin{equation}
			{\widehat{\mathscr{P}}}_X^{(2m)}(\xi) = \left( \sum_{n_1,\ldots n_{2m} : \sum_u u n_u = 2m} \frac{(2m)!}{n_1 ! \ldots n_{2m}!} \prod_{v = 1}^{2m} \left( \frac{\Psi_{\varphi}^{(v)}(\xi)}{v!}\right)^{n_v}\right) {\widehat{\mathscr{P}}}_X(\xi).
		\end{equation}
		Exploiting that $\Psi^{(v)}_\varphi (0) = \left( \int_{ \mathbb{R}^d} (\varphi(\bm{x}))^v \mathrm{d} \bm{x}\right) \Psi^{(v)}(0)$ we obtain the bound, for $\xi= 0$,
		\begin{equation} \label{eq:boundCFXvarphi}
			\abs{{\widehat{\mathscr{P}}}_X^{(2m)}(0)} \leq C'  \sum_{n_1,\ldots n_{2m} : \sum_u u n_u = 2m} \prod_{v=1}^{2m} \left(\int_{ \mathbb{R}^d} \abs{\varphi(\bm{x})}^v \mathrm{d} \bm{x}\right)^{n_v}		
		\end{equation}
with $C'>0$ a constant. 

We now apply \eqref{eq:boundCFXvarphi} to $\varphi = \psi_{j,G}$. 
Since we have
	\begin{equation}
		\int_{ \mathbb{R}^d}  \abs{\psi_{j,G}(\bm{x})}^v \mathrm{d} \bm{x} = 2^{j d v / 2} \int_{ \mathbb{R}^d} \abs{\psi_G(2^j \bm{x} - \bm{k})}^v \mathrm{d} \bm{x} = 2^{jd ( v/2 - 1)} \int_{ \mathbb{R}^d} \abs{ \psi_G(\bm{x})}^v \mathrm{d} \bm{x},
	\end{equation}
	we deduce from \eqref{eq:boundCFXvarphi} the new bound
	\begin{align}
	\mathbb{E} [\langle w, \psi_{j,G,\bm{k}} \rangle^{2m}] &= \abs{{\widehat{\mathscr{P}}}_{\langle w, \psi_{j,G,\bm{k}}\rangle}^{(2m)} (0) }  \nonumber \\
				&\leq  C''  \sum_{n_1,\ldots n_{2m} : \sum_u u n_u = 2m} \prod_{v=1}^{2m} 
						2^{jd( v/2 - 1)n_v} \nonumber \\
				& = C'' \sum_{n_1,\ldots n_{2m} : \sum_u u n_u = 2m} 2^{j d  \sum_v \left(n_v ( v/2 - 1)\right)}.
	\end{align}	
	Finally, since $\sum_v v n_v = 2m$ and $\sum_v n_v \geq 1$, we have $ \sum_v \left(n_v ( v/2 - 1)\right) \leq m -1$, and therefore
	\begin{equation}
		\mathbb{E} [\langle w, \psi_{j,G,\bm{k}} \rangle^{2m}]  \leq C 2^{j d (m-1)}
	\end{equation}	
	for an adequate $C>0$, as expected.
	\end{proof}

	Until now, we gave moment estimates by proposing upper bounds for the quantity $\mathbb{E}[ \abs{\langle w, \varphi\rangle}^p]$ (see Lemmas \ref{lemma:momentPoisson} and \ref{lemma:boundcoefficientfinitemoments}, but also Theorem 2 in \cite{Fageot2017multidimensional}). This allows in particular to identify in which Besov space is $w$. We now address the following problem: Can we lower bound $\mathbb{E}[ \abs{\langle w, \varphi\rangle}^p]$ with the moments of $\varphi$? Lemma \ref{lemma:lowerbound} answers positively to this question.
	
	\begin{lemma} \label{lemma:lowerbound}
		Let $w$ be a L\'evy white noise whose  indices satisfy $0 < \ubar{\beta}_\infty$ and $\beta_\infty  < p_{\max}$.
		We also fix a nontrivial bounded measurable function $\psi$ (typically a Daubechies wavelet).  
		Consider an integrability   $0<p< \beta_\infty$ satisfying the conditions
		Then, for $\epsilon>0$ small enough, there exists $A, B > 0$ such that
		\begin{equation} \label{eq:boundsonmoments}
			A 2^{- j\epsilon} 2^{j d p (1 / 2 - 1 / \ubar{\beta}_\infty)} \leq \mathbb{E} [ \abs{\langle w , \psi_{j,\bm{k}} \rangle}^p] \leq B 2^{j\epsilon}  2^{j d p (1 / 2 - 1 / \beta_\infty)} 
		\end{equation}		
		for any $j\geq 0$ and $\bm{k} \in  \mathbb{Z}^d$.
	\end{lemma}

	\begin{proof}
		As a preliminary remark, the shift parameter $\bm{k}$ in \eqref{eq:boundsonmoments} can be omitted since $w$ is stationary. We also remark that the upper bound of \eqref{eq:boundsonmoments} has already been proven~\cite[Corollary 1]{Fageot2017multidimensional}, where the conditions $p<\beta_\infty < p_{\max}$ are required. Actually, \cite{Fageot2017multidimensional} does not consider the index $p_{\max}$ and distinguishes between the conditions $\beta_\infty < \beta_0<2$ and $\beta_\infty \leq \beta_0 = 2$ with finite variance. These two scenarios are covered by the condition $\beta_\infty < p_{\max}$. 		
		Hence we focus on the lower bound.
		
		Because $p<2$, one can use the following representation of the $p$th-moment of $\langle w, \varphi\rangle$, which can be found for instance in \cite{Laue1980remarks,Deng2015shift}:
		\begin{equation}\label{eq:pthmomentRV}
		\mathbb{E}[ \abs{\langle w, \varphi\rangle}^p] = c_p \int_{ \mathbb{R}} \frac{1 - \Re {\widehat{\mathscr{P}}}_{\langle w, \varphi\rangle} (\xi) }{\abs{\xi}^{p+1}} \mathrm{d}\xi
		\end{equation}
		for some explicit constant $c_p>0$.
	For a short proof, see \cite[Proposition 6]{Fageot2017multidimensional}.
	We then remark that
	\begin{equation}
		 \Re {\widehat{\mathscr{P}}}_{\langle w, \varphi\rangle} (\xi)\leq \abs{{\widehat{\mathscr{P}}}_{\langle w, \varphi\rangle} (\xi)} = \abs{{\widehat{\mathscr{P}}}_w(\xi \varphi)} = \abs{\exp\left(\int_{ \mathbb{R}^d} \Psi(\xi\varphi(\bm{x})) \mathrm{d}\bm{x} \right) } = \exp\left( \int_{ \mathbb{R}^d} \Re\Psi (\xi \varphi(\bm{x})) \mathrm{d}\bm{x} \right).
	\end{equation}
	The sector condition \eqref{eq:sectorcondition} implies that one can found $c_{\Psi}>0$ such that $- \Re \Psi(\xi) =  \abs{\Re \Psi} (\xi)  \geq c_{\Psi} \abs{\Psi(\xi)}$. Thus, one has that $\abs{{\widehat{\mathscr{P}}}_{\langle w, \varphi\rangle} (\xi)} \leq \exp\left( -c_{\Psi} \int_{ \mathbb{R}^d} \abs{\Psi(\xi \varphi(\bm{x})) }\mathrm{d}\bm{x} \right)$, and then	
	\begin{equation} \label{eq:lowerbound1}
		\mathbb{E}[ \abs{\langle w, \varphi\rangle}^p]  \geq   c_p \int_{ \mathbb{R}} \frac{1 - \mathrm{e}^{ -c_{\Psi} \int_{ \mathbb{R}^d} \abs{\Psi(\xi \varphi(\bm{x})) }\mathrm{d}\bm{x}} }{\abs{\xi}^{p+1}} \mathrm{d}\xi.
	\end{equation}
	Now, by definition of the index $\ubar{\beta}_\infty$, the function $\abs{\Psi}$ is dominating 
	$\abs{\xi}^{\ubar{\beta}_\infty - \delta}$ at infinity 
	for $0< \delta < \ubar{\beta}_\infty$ arbitrarily small. 
	This domination together with the continuity of the functions $\abs{\Psi}$ and $\abs{\cdot}^{\ubar{\beta}_\infty - \delta} \One_{\abs{\cdot}>1}$ over $[1,\infty)$,  imply  the existence of $C>0$ such that, for any $\xi \in  \mathbb{R}$, $\abs{\Psi(\xi)} \geq C \abs{\xi}^{\ubar{\beta}_\infty - \delta} \One_{\abs{\xi}>1}$. In particular, we have that
	\begin{equation}
		\int \abs{\Psi(\xi(\varphi(\bm{x}))} \mathrm{d}\bm{x} \geq C \abs{\xi}^{\ubar{\beta}_\infty - \delta} \int_{\abs{\xi \varphi(\bm{x})}>1} \abs{\varphi(\bm{x})}^{\ubar{\beta}_\infty - \delta}.
	\end{equation}
	For simplicity, we write $\beta = \ubar{\beta}_\infty - \delta$ in this proof. 
	We now set $\varphi = 2^{jd/2} \psi(2^{j} \cdot)$, then we have by chance of variable $\bm{x} = 2^j \bm{x}$ that 
	\begin{equation} \label{eq:expressmythismoment}
	\int_{\abs{\xi \varphi(\bm{x})}>1} \abs{\varphi(\bm{x})}^{\beta} = 2^{j d \beta \left( \frac{1}{2}- \frac{1}{\beta}\right)} \int_{\abs{\psi(\bm{y})} > 2^{-jd/2} \abs{\xi}^{-1}} \abs{\psi(\bm{y})}^{\beta} \mathrm{d}\bm{y}.
	\end{equation}  
	
	Let $\gamma \geq 0$ and set $F(\gamma) = \int_{\abs{\psi(\bm{y})} > \gamma} \abs{\psi(\bm{y})}^{\beta} \mathrm{d}\bm{y}$. Then, $F$ is clearly decreasing, hence for every $0\leq \gamma \leq \lVert \psi \rVert_\infty / 2$ (this last quantity is finite by assumption on $\psi$), we have that $F(\gamma) \geq F(\lVert \psi\rVert_\infty / 2)$. Applying this with $\gamma = \xi 2^{jd/2}$ to \eqref{eq:expressmythismoment}, one deduces that 
	\begin{equation} \label{eq:crazystuff}
	  \int_{\abs{\psi(\bm{y})} > 2^{-jd/2} \abs{\xi}^{-1}} \abs{\psi(\bm{y})}^{\beta} \mathrm{d}\bm{y} \geq   
	   \rVert \psi \One_{\abs{\psi} > \lVert \psi \rVert_\infty / 2} \rVert_{\beta}^{\beta} \One_{\abs{\xi} \geq 2^{-jd/2} 2 / \lVert \psi \rVert_\infty}.
	\end{equation}
	
	Restarting from \eqref{eq:lowerbound1}, we then have
	\begin{align}
	\mathbb{E}[ \abs{\langle w, 2^{jd/2} \psi(2^j \cdot) \rangle}^p] &\geq 
	c_p \int_{ \mathbb{R}} \frac{1 - \exp\left( - c_{\Psi} C \abs{\xi}^\beta 2^{j d \beta  \left( \frac{1}{2}- \frac{1}{\beta}\right)} \int_{\abs{\psi(\bm{y})} > 2^{-jd/2} \abs{\xi}^{-1}} \abs{\psi(\bm{y})}^{\beta} \mathrm{d}\bm{y} \right)}{\abs{\xi}^{p+1}} \mathrm{d}\xi \nonumber \\
	&\overset{(i)}{\geq}
	c_p \int_{ \mathbb{R}} \frac{1 - \exp\left( - c_{\Psi} C \abs{\xi}^\beta  2^{j d \beta  \left( \frac{1}{2}- \frac{1}{\beta}\right)} 
   \rVert \psi \One_{\abs{\psi} > \lVert \psi \rVert_\infty / 2} \rVert_{\beta}^{\beta} \One_{\abs{\xi} \geq 2^{-jd/2} 2 / \lVert \psi \rVert_\infty}	
	 \right)}{\abs{\xi}^{p+1}} \mathrm{d}\xi \nonumber \\
	 &\overset{(ii)}{=}
	c_p c_{\Psi}^{p / \beta} C^{p/\beta} \rVert \psi \One_{\abs{\psi} > \lVert \psi \rVert_\infty / 2} \rVert^p_{\beta} 2^{j d \left( \frac{1}{2}- \frac{1}{\beta}\right)} 
	\int_{ \mathbb{R}} \frac{1 - \mathrm{e}^{ - \abs{u}^\beta} \One_{\abs{u} > 2^{- j d / \beta} 2 / \lVert \psi\rVert_\infty}}{\abs{u}^{p+1}} \mathrm{d}u \nonumber \\
	& \overset{(iii)}{\geq}
	c_p c_{\Psi}^{p / \beta} C^{p/\beta} \rVert \psi \One_{\abs{\psi} > \lVert \psi \rVert_\infty / 2} \rVert^p_{\beta} 2^{j d \left( \frac{1}{2}- \frac{1}{\beta}\right)} 
	\int_{ \mathbb{R}} \frac{1 - \mathrm{e}^{ - \abs{u}^\beta} \One_{\abs{u} >  2 / \lVert \psi\rVert_\infty}}{\abs{u}^{p+1}} \mathrm{d}u 
	\end{align}
	where we used \eqref{eq:crazystuff} in $(i)$, the change of variable $u = \xi c_{\Psi}^{1 / \beta} \rVert \psi \One_{\abs{\psi} > \lVert \psi \rVert_\infty / 2} \rVert_{\beta} 2^{j d  \left( \frac{1}{2}- \frac{1}{\beta}\right)}$ in $(ii)$, and the inequality $ \One_{\abs{u} > 2^{- j d / \beta} 2 / \lVert \psi\rVert_\infty} \geq  \One_{\abs{u} >  2 / \lVert \psi\rVert_\infty}$ in $(iii)$, valid for any $j\geq 0$. 
		Finally, we have shown that $\mathbb{E}[ \abs{\langle w, 2^{jd/2} \psi(2^j \cdot) \rangle}^p] \geq A 2^{j d  \left( \frac{1}{2}- \frac{1}{\beta}\right)}$ with $A>0$ a constant independent on $j$.
		
		To conclude, we remark that, for $\epsilon >0$ fixed, one can find $\delta>0$ small enough such that $2^{-j\epsilon} 2^{j d  \left( \frac{1}{2}- \frac{1}{\ubar{\beta}_\infty}\right)} \leq2^{j d  \left( \frac{1}{2}- \frac{1}{\ubar{\beta}_\infty - \delta}\right)}$ for any $j \geq 0$, giving the lower bound in \eqref{eq:boundsonmoments}. 
	\end{proof}

	
	\subsection{The Besov regularity of finite moment white noises}	
	
	We have all the tools to deduce the following results on the Besov regularity of finite moment white noises. 
	
		\begin{proposition} \label{prop:finiteMomentsBesov}
	Let $w$ be a finite moments L\'evy white noise with Blumenthal-Getoor indices $0 \leq \ubar{\beta}_\infty \leq \beta_\infty \leq 2$.
		Then, $w$ is 
\begin{itemize}
	\item almost surely in $B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau < d/\max(p,\beta_\infty) - d$ and $\rho < - d /p$, for $0<p\leq 2$, $p$ an even integer, or $p=\infty$; and
	\item  almost surely not in $B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau >  d/\max(p,\ubar{\beta}_\infty) - d$ or $\rho \geq  - d /p$ for every $0<p \leq \infty$
\end{itemize}
	\end{proposition}
	
	\begin{proof}
	We only treat the case $p< \infty$. For $p=\infty$, the result is obtained using embeddings (with $p=2m$, $m\rightarrow \infty$ for positive results and $p\rightarrow \infty$ for negative results). \\
	
	\textbf{ If $\tau < d/\max(p,\beta_\infty) - d$ and $\rho < - d /p$.}
	We first remark that for $p \leq 2$, we have $\rho < - \frac{d}{\min(p,2, p_{\max})} =   - \frac{d}{\min(p, p_{\max})}$, and the result is a consequence of our previous work \cite[Theorem 3]{Fageot2017multidimensional}. We can therefore assume that $p = 2m$  with $m \geq 1$ an integer. Then, $p=2m \geq 2 \geq \beta_\infty$, and the conditions on $\tau$ and $\rho$ become  $\tau < \frac{d}{2m} - d$ and $\rho < - \frac{d}{2m}$.
	Thanks to Lemma \ref{lemma:boundcoefficientfinitemoments}, we have
	\begin{align} \label{eq:seriestobefinite}
		\mathbb{E} \left[\lVert w \rVert_{B_{2m}^\tau( \mathbb{R}^d;\rho)}^{2m}  \right]
			&= \sum_{j\geq 0} 2^{j (2m\tau - d + dm)} \sum_{G,\bm{k}} \langle 2^{-j}\bm{k} \rangle^{2m\rho} \mathbb{E} [\langle w, \psi_{j,G,\bm{k}} \rangle^{2m}] \nonumber \\
			& \leq C 2^d \sum_{j \geq 0} 2^{j (2m\tau - d + 2dm)} \frac{1}{2^jd} \sum_{\bm{k} \in  \mathbb{Z}^d}  \langle 2^{-j}\bm{k} \rangle^{2m\rho}. 
	\end{align}
	Then, 
	$ \frac{1}{2^jd} \sum_{\bm{k} \in  \mathbb{Z}^d}  \langle 2^{-j}\bm{k} \rangle^{2m\rho} \underset{j \rightarrow \infty}{\longrightarrow} \int_{ \mathbb{R}^d} \langle \bm{x} \rangle^{2 m \rho} \mathrm{d} \bm{x} < \infty$ 
	since $2 m \rho < - d$, 
	and the series in \eqref{eq:seriestobefinite} is finite if and only if 
	$2 m \tau - d + 2dm < 0$, as we assumed. 
	Finally, we have shown that $\mathbb{E} \left[\lVert w \rVert_{B_{2m}^\tau( \mathbb{R}^d;\rho)}^{2m} \right] < \infty$, so that $w \in B_{2m}^\tau( \mathbb{R}^d;\rho)$ almost surely.  \\
	
	\textbf{If $\tau > d / p - d$.}
	This part of the proof is actually valid for any L\'evy white noise. It uses the decomposition $w= w_1 + w_2$ with $w_1$ a non-trivial compound Poisson noise and $w_2$ a white noise with finite moments (see Proposition \ref{prop:decomposition}). The main idea is that the jumps of the compound Poisson part are enough to make explode the Besov norm. 

	One writes $w_1 = \sum_{n\in \mathbb{Z}} a_n \delta(\cdot - \bm{x}_n)$ as in \eqref{eq:formPoisson}. Then, almost surely, all the $a_n$ are nonzero. Let $m>0$ be such that $\mathscr{P} (\abs{a_0} \geq m) > 0$. One can assume without loss of generality that $m=1$ (otherwise, consider the white noise $w/m$ that has the same Besov regularity than $w$). Then, there exists almost surely $n \in  \mathbb{Z}$ such that $\abs{a_n} \geq 1$ (Borel-Cantelli argument using that the $a_n$ are independent). For simplicity, one reorders the jumps such that this is achieved for $n=0$; that is, $\abs{a_0} \geq 1$. 
	
	Fix $I^d = [a,a+1)^d$ a set such that the Daubechies wavelet with gender $M^d$ satisfies $\abs{\psi_{M^d} (\bm{x}) } \geq C >0$ for some constant $C>0$ and every $\bm{x} \in I^d$. Due to the size of $I^d$, there exists for any $j \geq 0$ a unique $\bm{k}_j \in  \mathbb{Z}^d$ such that $2^j \bm{x}_0 - \bm{k}_j \in I^d$. Moreover, there exists almost surely a (random) integer $j_0 \geq 0$ such that for any $j \geq j_0$, $\psi_{M^d}(2^j \bm{x}_n - \bm{k}_j) = 0$ for any $n \neq 0$. This is due to the fact that there are (almost surely) finitely many $\bm{x}_n$ in $I^d$. From these preparatory considerations, one has that, for $j \geq j_0$,
	\begin{equation} \label{eq:lowerboundpoissonpart}
		\abs{ \langle w_1 , \psi_{M^d}(2^j \cdot - \bm{k}_j ) \rangle}
			= \abs{\sum_{n\in  \mathbb{Z}} a_n \psi_{M^d}(2^j \bm{x}_n - \bm{k}_j ) }
			= \abs{ a_0 \psi_{M^d}(2^j \bm{x}_0 - \bm{k}_j ) }
			\geq C
	\end{equation}
	almost surely.
	
	Since $w_2$ has a finite variance, we have, using the Markov inequality
	\begin{equation}
		\mathscr{P} 
		\left( \abs{\langle w_2 , \psi_{M^d}(2^j \cdot - \bm{k}_j) \rangle } \geq \frac{C}{2} \right) 
		\leq \frac{4 \mathbb{E} \left[\langle w_2 , \psi_{M^d}(2^j \cdot - \bm{k}_j) \rangle^2 \right]}{C^2} = \frac{4 \sigma_0^2 \lVert \psi_{M^d} \rVert_2^2}{C^2} 2^{-jd},
	\end{equation}
	where $\sigma_0^2$ is the variance of $w_2$ such that $\mathbb{E}[\langle w_2 , \varphi \rangle^2 ] = \sigma_0^2 \lVert \varphi \rVert_2^2$. In particular, $$\sum_{j\geq 0} \mathscr{P} \left( \abs{\langle w_2 , \psi_{M^d}(2^j \cdot - \bm{k}_j) \rangle } \geq {C} / {2} \right) < \infty.$$ From a new Borel-Cantelli argument, we readily deduce the existence of integers $j$ arbitrarily large such that $ \abs{\langle w_2 , \psi_{M^d}(2^{j} \cdot -  \bm{k}_j) \rangle } < {C} / {2}$. In particular, for $j \geq j_1$, we have, due to \eqref{eq:lowerboundpoissonpart}, $\abs{\langle w , \psi_{M^d} (2^{j} \cdot -  \bm{k}_j )\rangle } \geq C/ 2$.
	
	Putting the pieces together, we can now lower bound the Besov norm of $w$ by keeping only the mother wavelet $\psi_{M^d}$, the scale $j \geq j_1$ and the corresponding shift parameter $ \bm{k}_j$. Then, we obtain the almost sure lower bound
	\begin{align}
		\lVert w \rVert_{B_p^\tau( \mathbb{R}^d;\rho)}^p &\geq 2^{j (\tau p - d + dp)} \abs{\langle w , \psi_{M^d}(2^{j} \cdot -  \bm{k}_j) \rangle }^p   \geq   2^{j (\tau p - d + dp)} \left(C/2 \right)^p.
	\end{align} 
	Finally, this is valid for any $j \geq 	j_1$ and because $\tau p - d + dp > 0$, one concludes that $\lVert w \rVert_{B_p^\tau( \mathbb{R}^d;\rho)}^p = \infty$ almost surely. \\
	
	
	\textbf{If $0 < p < \ubar{\beta}_\infty$ and $d / \ubar{\beta}_\infty - d < \tau < d / p - d$.} Together with the  case $\tau > d / p - d$ above and by embeddings, we then deduce the result for $\tau < d / \max(\ubar{\beta}_\infty , p )  - d$.
	Note that the case $\ubar{\beta}_\infty=0$ is already covered by the previous case.
	Moreover, as soon as $f \notin B_p^{\tau}( \mathbb{R}^d;\rho)$ for some $p>0$, we also have that $f \notin B_q^{\tau + \epsilon}( \mathbb{R}^d;\rho)$ for any $q>p$ and $\epsilon > 0$ (see Figure \ref{fig:embeddings}). A crucial consequence for us is that it suffices to work with arbitrarily small $p$ in order to obtain the negative result we expect. We assume here that
	\begin{equation} \label{eq:troisconditionsurp}
	 p < \ubar{\beta}_\infty / 2, \quad p < \frac{\ubar{\beta}_\infty \beta_\infty}{2(\beta_\infty - \ubar{\beta}_\infty)}
	\end{equation}
	Note that the latter condition is simply $p< \infty$ (that is, no restriction) when $\ubar{\beta}_\infty = \beta_\infty$. 

	Let $k_0 \geq 1$ be such that the families of random variables $(\langle w , \Psi_{j, G, \bm{k}}\rangle)_{\bm{k} \in k_0  \mathbb{Z}^d}$ are independent when $j \geq 0$ and $G \in \mathrm{G}^j$ are fixed. This is possible because the wavelets are compactly supported and because $\langle w , \varphi_1 \rangle$ and $\langle w , \varphi_2 \rangle$ are independent as soon as $\varphi_1$ and $\varphi_2$ have disjoint support.
	By restricting the range of $\bm{k}$ and the gender to $G = M^d$, we have that
	\begin{equation} \label{eq:lowernormbesov}
		\lVert w \rVert_{B_p^\tau( \mathbb{R}^d;\rho)}^p \geq  C \sum_{j \geq 0} 2^{j ( \tau p - d + dp /2)} \sum_{\bm{k}\in k_0 \mathbb{Z}^d, 0 \leq k_i < k_0 2^j} \abs{\langle w,  \psi_{j,M^d,\bm{k}} \rangle}^p,
	\end{equation}
	with $C = \inf_{\lVert \bm{x}\rVert_\infty \leq k_0} \langle \bm{x} \rangle^{\rho} > 0$ a constant such that   $\langle 2^{-j} \bm{k} \rangle \geq C$.
	We set $X_{j,\bm{k}} = 2^{jd \left( \frac{1}{\ubar{\beta}_\infty} - \frac{1}{2}\right)}    {\langle w, \psi_{j,M^d,\bm{k}}\rangle}$ and
	$$M_{j,p} := 2^{-jd} \sum_{\bm{k}\in k_0 \mathbb{Z}^d, 0 \leq k_i < k_0 2^j} \abs{X_{j,\bm{k}}}^p,$$
	which is an average among $2^{jd}$ random variables. 
	
	Recall that $p < p_{\max} / 2, \ubar{\beta}_\infty / 2$. Moreover, since all the moments are finite, $p_{\max} = \infty > \beta_\infty$. Hence, one can apply Lemma \ref{lemma:lowerbound} with integrability   $q = p$ or $q = 2p$. There exists $\epsilon>0$ that can be choosen arbitrarily small and constants $m_q, M_q$ such that
	\begin{equation}
	m_q 2^{-j \epsilon} 2^{jqd \left(\frac{1}{2} - \frac{1}{\ubar{\beta}_\infty}\right)}
	\leq 
	\mathbb{E} \left[\abs{\langle w , \psi_{j,M^d,\bm{k}}\rangle}^q\right]
	\leq 
	M_q 2^{j \epsilon} 2^{jqd \left(\frac{1}{2} - \frac{1}{\beta_\infty}\right)}
	\end{equation}
	for any $j\geq 0, \bm{k}\in  \mathbb{Z}^d$.
	In particular, with our notations, we have that
	\begin{equation} \label{eq:boundtoprove}
		m_q 2^{-j \epsilon} \leq \mathbb{E} [M_{j,q}] \leq 
			M_q 2^{j \epsilon} 2^{jqd \left(\frac{1}{\ubar{\beta}_\infty} - \frac{1}{\beta_\infty}\right)}
	\end{equation}
	for any $j, \bm{k}$ and for $q=p$ or $2p$.
	
	Then, we control the variance of $M_{j,q}$ as follows (here, the sum over $\bm{k}$ are as in \eqref{eq:lowernormbesov}; that is, over $2^{jd}$ terms):
	\begin{align}\label{eq:boundtouse}
		\mathrm{Var}(M_{j,p}) 
		&= 
		\mathbb{E} [(M_{j,p} - \mathbb{E}[M_{j,q}])^2 
		\overset{(i)}{=}
2^{-jd} \mathbb{E} \left[2^{-jd}  \left( \sum_{\bm{k}} (\abs{X_{j,\bm{k}}}^p - \mathbb{E}[\abs{X_{j,\bm{k}}}^p]) \right)^2 \right] \nonumber  \\
		&\overset{(ii)}{=}
		 2^{-jd} \mathbb{E} \left[2^{-jd}   \sum_{\bm{k}} \left( (\abs{X_{j,\bm{k}}}^p - \mathbb{E}[\abs{X_{j,\bm{k}}}^p]) \right)^2 \right] \nonumber 
		 \overset{(iii)}{\leq}		 
		 2^{-jd} \mathbb{E} \left[2^{-jd}   \sum_{\bm{k}}  (\abs{X_{j,\bm{k}}}^{2p} \right] \nonumber \\
		 &= 2^{-jd} \mathbb{E}[M_{j,2p}],
	\end{align}
	where we used that $	 \mathbb{E}[M_{j,p}] =   \mathbb{E}[ \abs{X_{j,\bm{k}}}^p ]$ for every $\bm{k}$ in $(i)$, the independence of the $X_{j,\bm{k}}$ in $(ii)$, and the relation $\mathrm{Var}(X) \leq \mathbb{E} [X^2]$ in $(iii)$. 
	
	We now apply the Chebyshev's inequality $\mathscr{P}( \abs{X- \mathbb{E}[X]} \geq x )\leq \mathrm{Var} X / x^2$ to $X=M_{j,p}$ and $x = 2^{-j\epsilon} m_p / 2$ to get
	\begin{align} \label{eq:trucmachin}
	\mathscr{P}( \abs{M_{j,p}- \mathbb{E}[M_{j,p}]} \geq 2^{-j\epsilon} m_p / 2 )
	&\leq \frac{4 \mathrm{Var}(M_{j,p})}{m_p^2} 2^{2j \epsilon} \nonumber \\
	&\leq \frac{4 }{m_p^2} 2^{-jd} \mathbb{E}[M_{j,2p}]2^{2j \epsilon} \nonumber \\
	&	\leq 
	\frac{4 M_{2p}}{m_p^2} 2^{j \left( d ( 2p ( 1 / \ubar{\beta}_\infty - 1 / \beta_\infty) - 1 ) +  3\epsilon \right)}
	\end{align}
	where we used \eqref{eq:boundtouse} and \eqref{eq:boundtoprove} (applied to $q=2p$). 
	Due to the second relation in \eqref{eq:troisconditionsurp}, the exponent in \eqref{eq:trucmachin} is strictly negative for $\epsilon$ small enough (what we assume from now). Hence, $\sum_j \mathscr{P}( \abs{M_{j,p}- \mathbb{E}[M_{j,p}]}\geq 2^{-j\epsilon} m_p / 2 ) < \infty$. From the lemma of Borel-Cantelli, only a finite number of the $j$ can therefore satisfy the relation $\abs{M_{j,p}- \mathbb{E}[M_{j,p}]}\geq 2^{-j\epsilon} m_p / 2$. A direct consequence is then that there exists almost surely $J \geq 0$ such that for every $j \geq J$, 
		\begin{equation}
			M_{j,p} \geq \mathbb{E} [M_{j,p}] - \abs{M_{j,p} - \mathbb{E} [M_{j,p}] } \geq m_p 2^{-j\epsilon} - \frac{m_p}{2} 2^{-j\epsilon} = \frac{m_p}{2} 2^{-j \epsilon}.
		\end{equation} 
	Using the latter inequality, we deduce that
\begin{equation*}
	\lVert w \rVert_{B_p^\tau( \mathbb{R}^d;\rho)}^p 
	\geq 
	C \sum_{j\geq J} 2^{j p (\tau - d + d / \ubar{\beta}_\infty)} M_{j,p} 
	\geq 
	\frac{C m_p}{2} \sum_{j\geq J} 2^{j p (\tau - d + d / \ubar{\beta}_\infty- \epsilon)}.
	\end{equation*}		
	For $\epsilon$ small enough, we have that $\tau + d - d /\ubar{\beta}_\infty - \epsilon > 0$ and therefore $\lVert w \rVert_{B_p^\tau( \mathbb{R}^d;\rho)}^p  = \infty$ almost surely. 	\\
	
	\textbf{ If $\rho \geq  - d /p$.}
	The proof is identical to the one of compound Poisson noises where we also have treat the case $\rho \geq - d / p$. 
	\end{proof}
	

\section{Non-Gaussian L\'evy White Noise} \label{sec:general}

This section is the occasion to synthesized the results and to deduce the general case from the previous ones. We consider non-Gaussian L\'evy white noises, the Gaussian white noise having already been treated in Section \ref{sec:Gaussian}.
We conclude with the proof of Theorem \ref{maintheo:regdecaylimit}.

\begin{proposition} \label{prop:sparsenoise}
	Consider a non-Gaussian L\'evy white noise $w$  with indices $0\leq \ubar{\beta}_\infty \leq \beta_\infty \leq 2$ and $0 < p_{\max} \leq \infty$.
		Then, $w$ is
		\begin{itemize}
		\item  almost surely in $B_p^\tau( \mathbb{R}^d;\rho)$ 
		if $\tau < d/\max(p,\beta_\infty) - d$ 
		and $\rho < - d /\min(p,p_{\max})$, for $0<p\leq 2$, $p$ an even integer, or $p=\infty$; and 
		\item almost surely not in $B_p^\tau( \mathbb{R}^d;\rho)$
		if $\tau > d/\max(p,\ubar{\beta}_\infty) - d$ 
		or $\rho > - d /\min(p,p_{\max})$ for every $0 < p \leq \infty$.
		\end{itemize}
\end{proposition}

\begin{proof}	
	According to Proposition \ref{prop:decomposition}, $w = w_1 + w_2$ with $w_1$ a compound Poisson noise, $w_2$ a L\'evy white noise with finite moments (that can include a Gaussian part). 
	Moreover, we have
		\begin{align}
		\beta_\infty 	&= \beta_\infty(w_2) \geq \beta_\infty(w_1) = 0, \\
		p_{\max} 	&= p_{\max} (w_1) \leq p_{\max}(w_2) =   \infty.
	\end{align}
	
	\textbf{If  $\tau < \left( \frac{d}{\max(p,\beta_\infty)} - d \right) $ and $\rho < - \frac{d}{\min(p,p_{\max})}$. }
	Then, $\tau < \left( \frac{d}{p} - d \right) $ and $\rho < - \frac{d}{\min(p,p_{\max}(w_1))}$ so that $w_1 \in B_{p}^{\tau}( \mathbb{R}^d;\rho)$ due to Proposition \ref{prop:PoissonBesov}. Similarly, $w_2 \in B_{p}^{\tau}( \mathbb{R}^d;\rho)$ thanks to Proposition \ref{prop:finiteMomentsBesov}. Finally, $w = w_1 + w_2 \in B_{p}^{\tau}( \mathbb{R}^d;\rho)$. \\
	
\textbf{If  $\tau > \left( \frac{d}{\max(p,\beta_\infty)} - d \right) $.} The arguments of Proposition \ref{prop:finiteMomentsBesov} for this case are still valid for $w$.	\\	
	
	\textbf{If  $\rho >  - \frac{d}{\min(p,p_{\max})}$. }
	If $\rho \geq - d / p$, then the same proof as for the compound Poisson case implies that $w \notin B_p^\tau( \mathbb{R}^d;\rho)$.
	We already know that $w \notin B_p^\tau( \mathbb{R}^d;\rho)$ if $\tau \geq d / \max(p,\beta_\infty) - d$ and if $\rho \geq - d /p$. The only remaining case is when $p < p_{\max}$, $ \tau < d / \max(p,\beta_\infty) - d$, and $- d / p_{\max} < \rho < - d / p$. In this case, $w_2 \in B_p^\tau( \mathbb{R}^d;\rho)$ from Proposition \ref{prop:finiteMomentsBesov}, while $w_1 \notin B_p^\tau( \mathbb{R}^d;\rho)$ with Proposition \ref{prop:PoissonBesov} due to the condition $\rho > - d / p_{\max}$. Finally, $w \notin B_p^\tau( \mathbb{R}^d;\rho)$ as the sum of a element of $B_p^\tau( \mathbb{R}^d;\rho)$ and a function not in $B_p^\tau( \mathbb{R}^d;\rho)$.
	\end{proof}

Finally, we can translate our results in terms of the local smoothness $\tau_p(w)$ and the asymptotic growth rate $\rho_p(w)$ of L\'evy white noises.

\begin{proof}[Proof of Theorem \ref{maintheo:regdecaylimit}]
	The values of $\tau_p(w)$ and $\rho_p(w)$ are directly deduced from Propositions \ref{prop:GaussBesov}, \ref{prop:PoissonBesov}, and  \ref{prop:sparsenoise}.  
	Positive results ($w \in B_p^{\tau}( \mathbb{R}^d;\rho)$) directly give lower bounds for $\tau_p(w)$ and $\rho_p(w)$, while negative results  ($w \notin B_p^{\tau}( \mathbb{R}^d;\rho)$) provide upper bounds. 
	For the Gaussian and compound Poisson cases, studied separatly, there are no restrictions on $0<p \leq \infty$. The case of a general non-Gaussian L\'evy white noise is deduced from Proposition \ref{prop:sparsenoise}
\end{proof}

\section{Discussion and Examples} \label{sec:blabla}

	\subsection{Application to some Subfamilies of L\'evy White Noises} \label{subsec:examples}
	
		We apply Theorem \ref{maintheo:regdecaylimit} to deduce the local smoothness and asymptotic growth rate of notorious L\'evy white noise. This includes Gaussian, symmetric-$\alpha$-stable, symmetric Gamma (including Laplace), compound Poisson, inverse Gaussian and layered stable white L\'evy nosies. All the underlying laws are known to be infinitely divisible \cite{Sato1994levy,Houdre2007layered}.
		In Table \ref{table:noises}, we define the different families in terms of the characteristic functional of $X= \langle w , \One_{[0,1]^d} \rangle$ and give adequate references.
		Most of these families together with the convention we are following in this paper are introduced with more details  in \cite[Section 5.1]{Fageot2016scaling}. 
		This is especially relevant for the layered stable white noises, which have the particularity of describing the complete spectrum of possible couples  $(\alpha_1,\alpha_2) = (\beta_\infty,p_{\max}) \in (0,2)^2$, and for which the L\'evy exponent is
		\begin{equation}\label{eq:layeredpsi}
		\Psi_{\alpha_1,\alpha_2} (\xi) = \int_{ \mathbb{R}} (\cos(t \xi) - 1 ) \left( \One_{\abs{t}\leq1} \abs{t}^{-(\alpha_1+1)} + \One_{\abs{t}>1} \abs{t}^{-(\alpha_2+1)} \right) \mathrm{d} t.
		\end{equation}	
	
We also provide a visualization of our results in terms of Triebel diagrams. In Figures \ref{fig:taurhogauss} to \ref{fig:taurhogenerallevy}, we plot the local smoothness $\frac{1}{p} \mapsto \tau_p(w)$ and asymptotic growth rate $\frac{1}{p} \mapsto \rho_p(w)$ for different L\'evy white noises (with the exception of $\tau_p(w)$ that is not fully determined for the general case in Figure \ref{fig:taurhogenerallevy}; here, we represent the lower and upper bounds of \eqref{eq:LNtaurho}). A given noise is almost surely in a Besov space $B_p^\tau( \mathbb{R}^d;\rho)$ is the points $(1/p,\tau)$ and $(1/p,\rho)$ are in the green regions. A contrario, the noise is almost surely not in $B_p^\tau( \mathbb{R}^d;\rho)$ if $(1/p,\tau)$ or $(1/p,\rho)$ are in the red region. In Figure~\ref{fig:taurhogenerallevy}, the white region corresponds to the case where we do not now if the L\'evy white noise is or not in the corresponding Besov spaces.

\begin{table*}[t!] 
\centering
\caption{L\'evy White Noises and their Indices}

\begin{tabular}{cccccc} 
\hline
\hline\\[-2ex] 
White noise & Parameters &  ${\widehat{\mathscr{P}}}_X(\xi)$  & $\beta_\infty = \ubar{\beta}_\infty$ & $p_{\max}$   \\
\\
[-2ex] 
\hline\\[-1.5ex]
Gaussian    & $\sigma^2>0$ & $\mathrm{e}^{- \sigma^2 \omega^2 /2}$ & $2$ & $\infty$  \\[+1ex]
Cauchy   \cite{Taqqu1994stable}  & $\gamma > 0$ & $\mathrm{e}^{- \gamma \lvert \xi \rvert}$    & $1$ & $1$   \\[+1ex]
S$\alpha$S  \cite{Taqqu1994stable}  & $0<\alpha<2$ & $\mathrm{e}^{-\lvert \xi \rvert^\alpha}$    & $\alpha$ & $\alpha$   \\[+1ex]
sum of S$\alpha$S	& $0 < \alpha_1 , \alpha_2 \leq 2$ & $\mathrm{e}^{-\lvert \xi \rvert^{\alpha_1} - \lvert \xi \rvert^{\alpha_2}}$ & $\max(\alpha_1, \alpha_2)$ & $\min(\alpha_1, \alpha_2)$   \\[+1ex]
Laplace \cite{Koltz2001laplace}  & $\sigma^2 > 0$ & $(1 + \sigma^2\xi/2)^{-1}$ & $0$ & $\infty$ \\[+1ex]
symmetric Gamma  \cite{Koltz2001laplace} & $\sigma^2 , \lambda > 0$ & $(1 + \sigma^2\xi/2)^{-\lambda}$ & $0$ & $\infty$ \\[+1ex]
finite-moments & $\lambda >0$ & $\mathrm{e}^{\lambda (\widehat{\mathscr{P}} (\xi) -1 )}$ & $0$ & $\infty$ \\
compound Poisson  & $\mathscr{P}$ & &  \\[+1ex]
layered stable \cite{Houdre2007layered} & $0 < \alpha_1 , \alpha_2 < 2$ &    see \eqref{eq:layeredpsi} & $\alpha_1$ & $\alpha_2$ \\[+1ex]
inverse Gaussian \cite{Barndorff1997processes} &  -  & $\mathrm{e}^{1 - (1 - 2 \mathrm{i} \xi)^{1/2}}$ & $1/2$ & $\infty$  \\[+1ex]
finite-moments & - & - & $2$ & $\infty$ \\
with Gaussian part  &   & &  \\
\hline
\hline
\end{tabular} \label{table:noises}
\end{table*}

\begin{figure}[h!] 
\centering
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,-0.75) -- (3,-0.75) -- (3,-2)  -- (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-0.75) -- (3,-0.75) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\tau$} ;
\draw[ thick,color=black] (-0.05,0) -- (0.05,0)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-1.5) -- (0.1,-1.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,-0.05) -- (1.5,0.1)  node[black,above] { $1$};
\draw[black, thick, ->](0,-0.75) --(3,-0.75);
\draw[ thick,color=black] (-0.05,-0.75) -- (0.05,-0.75)  node[black,left] { $-\frac{d}{2}$};

\end{tikzpicture}
\end{subfigure}
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,1) -- (3,-2)  -- (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,1) -- (3,-2) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,1)--(3,1) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\rho$} ;
\draw[ thick,color=black] (-0.05,1) -- (0.05,1)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-0.5) -- (0.1,-0.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,0.9) -- (1.5,1.1)  node[black,above] { $1$};
\draw[black, thick, ->](0,1) --(3,-2);

\end{tikzpicture}
\end{subfigure}
\caption{Gaussian noise} \label{fig:taurhogauss}

\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,-1.5) -- (2.25,0.75) -- (3,0.75)  -- (3,-2) --  (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-1.5) -- (2.25,0.75) -- (3,0.75) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\tau$} ;
\draw[ thick,color=black] (-0.05,0) -- (0.05,0)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-1.5) -- (0.1,-1.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,-0.05) -- (1.5,0.1)  node[black,above] { $1$};
\draw[ thick,color=black](2.25,-0.05) -- (2.25,0.1)  node[black,below] { $\frac{1}{\alpha}$};
\draw[black, thick, -](0,-1.5) --(2.25,0.75);
\draw[black, thick, ->](2.25,0.75) --(3,0.75);


\end{tikzpicture}
\end{subfigure}
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,-1.25) -- (2.25,-1.25) -- (3,-2)  -- (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-1.25) -- (2.25,-1.25) -- (3,-2) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,1)--(3,1) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\rho$} ;
\draw[ thick,color=black] (-0.05,1) -- (0.05,1)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.05,-1.25) -- (0.05,-1.25)  node[black,left] { $-\frac{d}{\alpha}$};
\draw[ thick,color=black](2.25,0.9) -- (2.25,1.1)  node[black,above] { $\frac{1}{\alpha}$};
\draw[ thick,color=black] (-0.1,-0.5) -- (0.1,-0.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,0.9) -- (1.5,1.1)  node[black,above] { $1$};
\draw[black, thick, -](0,-1.25) --(2.25,-1.25);
\draw[black, thick, ->](2.25,-1.25) --(3,-2);

\end{tikzpicture}
\end{subfigure}
\caption{S$\alpha$S noise with $\alpha = 2/3$}

\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,-1.5) -- (3,1.5)  -- (3,-2) --  (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-1.5)-- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\tau$} ;
\draw[ thick,color=black] (-0.05,0) -- (0.05,0)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-1.5) -- (0.1,-1.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,-0.05) -- (1.5,0.1)  node[black,above] { $1$};
\draw[black, thick, ->](0,-1.5)--(3,1.5);

\end{tikzpicture}
\end{subfigure}
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,0.30) -- (0.75,0.30)-- (3,-2)  -- (0,-2) -- cycle; 
\fill[red, opacity= 0.3](0,0.30) -- (0.75,0.30)-- (3,-2) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,1)--(3,1) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\rho$} ;
\draw[ thick,color=black] (-0.05,1) -- (0.05,1)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.05,0.30) -- (0.05,0.30)  node[black,left] { $-\frac{d}{p_{\max}}$};
\draw[ thick,color=black](0.75,0.9) -- (0.75,1.1)  node[black,above] { $\frac{1}{p_{\max}}$};
\draw[ thick,color=black] (-0.1,-0.5) -- (0.1,-0.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,0.9) -- (1.5,1.1)  node[black,above] { $1$};
\draw[black, thick, -](0,0.30) --(0.75,0.30);
\draw[black, thick, ->](0.75,0.30) --(3,-2);

\end{tikzpicture}
\end{subfigure}
\caption{Compound Poisson noise with $p_{\max} = 2$}
\end{figure}

\begin{figure}[h!] 
\centering
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]


\fill[green, opacity= 0.3] (0,-1.5) -- (3,1.5)  -- (3,-2) --  (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-1.5)-- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\tau$} ;
\draw[ thick,color=black] (-0.05,0) -- (0.05,0)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-1.5) -- (0.1,-1.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,-0.05) -- (1.5,0.1)  node[black,above] { $1$};
\draw[black, thick, ->](0,-1.5)--(3,1.5);

\end{tikzpicture}
\end{subfigure}
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,1) -- (3,-2)  -- (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,1) -- (3,-2) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,1)--(3,1) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\rho$} ;
\draw[ thick,color=black] (-0.05,1) -- (0.05,1)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-0.5) -- (0.1,-0.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,0.9) -- (1.5,1.1)  node[black,above] { $1$};
\draw[black, thick, ->](0,1) --(3,-2);

\end{tikzpicture}
\end{subfigure}
\caption{Symmetric-Gamma noise}

\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]


\fill[green, opacity= 0.3] (0,-1.5) -- (0.75,-0.75) -- (3,-0.75)  -- (3,-2) --  (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-1.5) -- (2.25,0.75) -- (3,0.75) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,0)--(3,0) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\tau$} ;
\draw[ thick,color=black] (-0.05,0) -- (0.05,0)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.1,-1.5) -- (0.1,-1.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,-0.05) -- (1.5,0.1)  node[black,above] { $1$};
\draw[ thick,color=black](2.25,-0.05) -- (2.25,0.1)  node[black,below] { $\frac{1}{\underline{\beta}_\infty}$};
\draw[ thick,color=black](0.75,-0.05) -- (0.75,0.1)  node[black,above] { $\frac{1}{\beta_\infty}$};
\draw[black, thick, -](0,-1.5) --(2.25,0.75);
\draw[black, thick, ->](2.25,0.75) --(3,0.75);
\draw[black,  thick, ->](0.75,-0.75) --(3,-0.75);
\end{tikzpicture}
\end{subfigure}
\begin{subfigure}[b]{0.30\textwidth}
\begin{tikzpicture}[x=2cm,y=2cm,scale=0.42]

\fill[green, opacity= 0.3] (0,-1.25) -- (2.25,-1.25) -- (3,-2)  -- (0,-2) -- cycle; 
\fill[red, opacity= 0.3] (0,-1.25) -- (2.25,-1.25) -- (3,-2) -- (3,1.5) -- (0,1.5) -- cycle; 
\draw[thick, ->] (0,1)--(3,1) node[circle,right] {$\frac{1}{p}$} ;
\draw[thick, ->] (0,-2)--(0,1.5) node[circle,above] {$\rho$} ;
\draw[ thick,color=black] (-0.05,1) -- (0.05,1)  node[black,left] { $0$};
\draw[ thick,color=black] (-0.05,-1.25) -- (0.05,-1.25)  node[black,left] { $-\frac{d}{p_{\max}}$};
\draw[ thick,color=black](2.25,0.9) -- (2.25,1.1)  node[black,above] { $\frac{1}{p_{\max}}$};
\draw[ thick,color=black] (-0.1,-0.5) -- (0.1,-0.5)  node[black,left] { $-d$};
\draw[ thick,color=black](1.5,0.9) -- (1.5,1.1)  node[black,above] { $1$};
\draw[black, thick, -](0,-1.25) --(2.25,-1.25);
\draw[black, thick, ->](2.25,-1.25) --(3,-2);

\end{tikzpicture}
\end{subfigure}
\label{fig:taurhogenerallevy}
\caption{L\'evy white noise with $\underline{\beta}_\infty = 2/3 < \beta_\infty = 2$ and $p_{\max}  = 2/3$}
\end{figure}


	\subsection{L\'evy White Noises with Distinct Blumenthal-Getoor Indices}
	\label{subsec:handcrafted} 
	
	In Theorem \ref{maintheo:regdecaylimit}, we obtained lower and upper bounds for the local regularity of a L\'evy white noise that are equal for every $p$ if and only if $\ubar{\beta}_\infty = \beta_\infty$. This equality is valid for all the examples presented in Section \ref{subsec:examples}.  It is however possible to find counterexamples, and it is even possible to define L\'evy exponents with indices taking any values $(\ubar{\beta}_\infty ,\beta_\infty) \in [0,2]^2$ with the obvious constraint that $\ubar{\beta}_\infty \leq \beta_\infty$. This is done in \cite[Examples 1.1.14, 1.1.15]{Farkas2001function}, where the authors introduce
	\begin{align}
		\Psi_{\beta_2,M}(\xi) &=    \sum_{k \geq 1} 2^{\beta_2 M^k - k} (\cos( 2^{-M^k}\xi) - 1), \\
		\Psi_{\beta_1,\beta_2,M}(\xi) &= \int_{\abs{t}\leq1} \frac{\cos(t \xi) - 1}{\abs{t}^{\beta_1+1}} \mathrm{d}t + \Psi_{\beta_2,M}(\xi).
	\end{align}
	with $0<\beta_1 \leq \beta_2 <2$ and $M > 2 / (2-\beta_2)$, and show that $\Psi_{\beta_M}$ ($\Psi_{\beta_1,\beta_2,M}$, respectively) is a L\'evy exponent with Blumenthal-Getoor indices $\ubar{\beta}_\infty = 0$ and $\beta_\infty = \beta_2$ ($\ubar{\beta}_\infty = \beta_1$ and $\beta_\infty = \beta_2$, respectively).

	\subsection{Conclusive Remarks and Open Questions}

	We have obtained new results on the localization of L\'evy white noises in weighted Besov spaces. This includes the identification of the local smoothness in many cases, and lower and upper bounds for the general case, together with the identification of the asymptotic growth rate. 
	
	Some open questions remain for a complete answer on the Besov regularity of L\'evy white noise. 
	\begin{itemize}
		\item Our general results are valid for integrability   $0<p<2$,  $p = 2m$ with $m \geq 1$ integer, or $p=\infty$. We conjecture that the derived formulas remind true for any $0<p \leq \infty$. 
		
		\item More importantly, we have distinct lower and upper bound for the local smoothness when $\ubar{\beta}_\infty \neq \beta_\infty$ in Theorem \ref{maintheo:regdecaylimit} for $\ubar{\beta}_\infty \leq p \leq \beta_\infty$. 
		For worst-case examples  of Section \ref{subsec:handcrafted}, this results in the least sharp bound $d / \max(p,2) - d \leq \tau_p(w) \leq d / p - d$. 
		The identification of $\tau_p(w)$ for these cases is unknown at this stage, and it is not sure that it can be expressed in terms of the indices considered so far.
		
		\item Finally, we did not investigate the localization of a general L\'evy white noise for the critical value $\tau = \tau_p(w)$ or $\rho = \rho_p(w)$. However, partial answers have been given for compound Poisson white noises (Proposition \ref{prop:PoissonBesov}) and finite moments white noises (Proposition \ref{prop:finiteMomentsBesov}), and a complete characterization was given in the Gaussian case (Proposition \ref{prop:GaussBesov}). For the general case, we conjecture that $w \notin B_{p}^{\tau}( \mathbb{R}^d;\rho)$ as soon as $\tau = \tau_p(w)$ or $\rho = \rho_p(w)$, in accordance with the known results.
	\end{itemize}

\section*{Acknowledgment}
This work has benefited from exchanges with John Paul Ward, Alireza Fallah and Virginie Uhlmann. We are also grateful to Robert Dalang,  St\'ephane Jaffard, and Ren\'e Schilling for enlightening discussions.
The research leading to these results has received funding from the European Research Council under Grant H2020-ERC (ERC grant agreement $\text{n}^\circ$ 692726 - GloblBioIm). 

\begin{small}
\bibliographystyle{plain}
