I have the function: $$ G_k(x)_=\frac{1}{(4\pi)^{k/2}\Gamma(k/2)}\int_0^\infty e^{-\pi|x|^2/\delta}e^{-\delta/4\pi}\delta^{-(n-k)/2}\,\frac{d\delta}{\delta}, $$ for all $x\in\mathbb{R}^n$ and $k>0$. I know that $G_k\in L^1(\mathbb{R}^n)$ and: $$ \mathcal{F}(G_k)(\xi)=(1+4\pi^2|\xi|^2)^{-k/2},\quad\forall\xi\in\mathbb{R}^n,k>0,$$ i want to prove that: $$ \mathcal{F}(G_{2k})(\xi/2\pi)=(2\pi)^k\mathcal{F}(G_{2k})(2\pi\xi),\quad\forall\xi\in\mathbb{R}^n,k>0.$$ I don't know how to go on. Any help would be appreciated.
$\mathcal{F}$ is the Fourier transform: for all $f\in L^1(\mathbb{R}^n)$, we set: $$\mathcal{F}(f)(\xi)=\int_{\mathbb{R}^n} f(x)e^{-2\pi i x\cdot \xi}\,dx\quad\forall\xi\in\mathbb{R}^n.$$