We consider the Gaussian heat kernel $p_t$ on $\mathbb R^d$, i.e., $$ p_t (x) := (4\pi t)^{-\frac{d}{2}} e^{-\frac{|x|^2}{4t}}, \quad t>0, x \in \mathbb R^d, $$ and define the operator $P_t$ by $$ P_t f (x) := \int_{\mathbb R^d} p_t ( x- y) f(y) \, \mathrm d y. $$
We fix $t>0,\delta \in (0, 1), y \in \mathbb R^d$ and $i,j \in \{1, 2, \ldots, d\}$. Let $$ I := \int_{\mathbb R^d} | \partial_i \partial_j (1- \Delta)^{- \frac{\delta}{2}} p_t ( \cdot - y) (x) | \, \mathrm d x. $$
At page $584$ of this paper, the author obtains an upper bound $$ I \lesssim t^{\frac{\delta}{2}-1} \tag{$*$} \label{*} $$ by using an auxiliary result (in the same paper)
Lemma 5.2(2) For any $\alpha, \beta, k \geq 0$, there exists a constant $c>0$ such that $$ \left\|(1-\Delta)^{-k} P_t f\right\|_{C_b^{\alpha+\beta}} \leq c t^{-\left(\frac{\alpha}{2}-k\right)^{+}}\|f\|_{C_b^\beta}, \quad t>0 . $$
The related function space is defined as
For any $n \in \mathbb{Z}^{+}$and $\alpha \in(0,1), C_b^{n+\alpha}\left(\mathbb{R}^d\right)$ is the space of functions $f \in C_b^n\left(\mathbb{R}^d\right)$ such that $$ \|f\|_{C_b^{n+\alpha}}:=\|f\|_{C_b^n}+\sup _{x \neq y} \frac{\left|\nabla^n f(x)-\nabla^n f(y)\right|}{|x-y|^\alpha}<\infty . $$
Unfortunately, I could not see how to get (\ref{*}). We have $$ \begin{align*} (1- \Delta)^{-k} P_t f &:= (1- \Delta)^{-k} (P_t f) \\ & = (1- \Delta)^{-k} \left ( \int_{\mathbb R^d} p_t ( \cdot - y) f(y) \, \mathrm d y \right ). \end{align*} $$
So the operator $\partial_i \partial_j (1- \Delta)^{- \frac{\delta}{2}}$ is inside the integral in $I$ whereas $(1- \Delta)^{-k}$ is outside the integral in $(1- \Delta)^{-k} P_t f$. I don't know how to put $I$ in a form that Lemma 5.2(2) is applicable.
Could you explain how to get (\ref{*})?
Thank you so much for your help!