no, it is not positive everywhere; notice that the integral only depends on $|x|$, so we may orient the vector $x$ along the $x_1$ axis; going to hyper spherical coordinates we have
$$\Phi_n(x)=\frac{2\pi^{(n-1)/2}}{\Gamma[\tfrac{1}{2}(n-1)]}\int_0^\infty dr \int_0^\pi d\phi\; e^{-r^n}\cos\bigl(|x|r\cos\phi\bigr)r^{n-1}(\sin\phi)^{n-2}$$
I checked that it becomes negative for $n=3$.