Consider the heat equation $$ (\partial_t + \Delta + V)u = 0$$ on a complete (open) Riemannian manifold with bounded geometry, where $V$ is a smooth and bounded potential.

Consider the semigroup generated by $A:=\Delta + V$ on the spaces $X = C_0(M)$ or $X = L^p(M)$ for $1\leq p < \infty$. The domain of $A$ is in each case $$ \mathrm{dom}(A) =\{ u \in X \mid Au \in X \},$$ where $Au$ denotes the distributional derivative.

My question is the following: It is clear by parabolic regularity that $e^{-tA}u$ is smooth for any $u \in X$. But is it also true that solutions $e^{-tA}u$ satisfy $$ |\nabla^ku| \in X ~~~~~~~\forall ~k$$ for $u \in \mathrm{dom}(A)$? The local statement is clear (i.e. solutions are in $W^{k, p}_{\mathrm{loc}}(M)$ and in $C^k(M)$ for every $k$), but is it also true that we have this decay at infinity?

What about the special case of vanishing $V$?