Let $f(x,y)$ be a non-negative function with $x,y \in \mathbb R^3$ that satisfies $$ I_1(f) := \iint_{\mathbb R^3 \times \mathbb R^3 } f(x,y) \, dx \ dy < \infty $$ and $$ I_2(f) := \iint_{\mathbb R^3 \times \mathbb R^3 } |y|^2 f(x,y) \, dx \ dy < \infty. $$ Also, define the sequence by $$ f_n(x,y) := \left(\frac n{4\pi} \right)^3 \iint_{|x'|<n,|y'|<n} e^{-\frac n4 (|x-x'|^2+|y-y'|^2)} f(x',y') \, dx' \, dy'. $$ Show that $$ I_2(f_n) \le 4I_2(f) + 16 \Gamma(5/2) \frac 1n I_1(f), $$ where $\Gamma$ is the Gamma function.
Honestly, all I can do is rewrite the expression of $I_2$ as follows: $$ I_2(f_n) = \left(\frac n{4\pi} \right)^3 \iint_{\mathbb R^3 \times \mathbb R^3} |y|^2 \left(\iint_{|x'|<n,|y'|<n} e^{-\frac n4 (|x-x'|^2+|y-y'|^2)} f(x',y') \, dx' \, dy'\right) \, dx \, dy. $$ There's a product-to-sum inequality being used here, I think.
Also, I'm okay if any answer to this question does not provide the explicit constants $4$ and $16 \Gamma(5/2)$; they are just positive constants.
Remark: This is the third item on the list of assertions of a lemma on page 142 of this paper.