In this paper, Lemma 6, Pinsky proves that $$H(x) =(4\pi t)^{-d/2} \int_{\mathbb{R}^d} e^{\frac{-|x-y|^2}{4t}}(1+|y|)^m \, dy$$ attains its maximum in $x=0$ for $m<0$. This can also be proven using the [Riesz rearrangement inequality][2]. My question at first was whether $$G(x) =(4\pi t)^{-d/2} \int_{\mathbb{R}^d} e^{\frac{-|x-y|^2}{4t}}|y|^k \, dy,$$ $k>0$, there is $x_0\in \mathbb{R}^d$ at which the integral reaches its maximum. This time it is not possible to use Riesz rearrangement inequality, as $|y|^k$, $k>0$ is non-decreasing.
In a comment of Christian Remling below, it was noted that $$(4\pi t)^{-d/2} \int_{\mathbb{R}^d} e^{\frac{-|x-y|^2}{4t}}|y|^k \, dy \rightarrow \infty \text{ as } |x|\rightarrow \infty.$$ So a question arose, whether there is any upper bound for $$\int_{\mathbb{R}^d} e^{\frac{-|x-y|^2}{4t}}|y|^k \, dy,$$ $k>0$, of polynomial type, so that I can control the growth of $G(x)$, possibly multiplying by $|x|^\beta$, for some $\beta <0$. Is there any result in this direction, or any suggestions?