Is the product of two Gaussian random variables also a Gaussian?
Solution 1:
The product of two Gaussian random variables is distributed, in general, as a linear combination of two Chi-square random variables:
$$ XY \,=\, \frac{1}{4} (X+Y)^2 - \frac{1}{4}(X-Y)^2$$
Now, $X+Y$ and $X-Y$ are Gaussian random variables, so that $(X+Y)^2$ and $(X-Y)^2$ are Chi-square distributed with 1 degree of freedom.
If $X$ and $Y$ are both zero-mean, then
$$ XY \sim c_1 Q - c_2 R$$
where $c_1=\frac{Var(X+Y)}{4}$, $c_2 = \frac{Var(X-Y)}{4}$ and $Q, R \sim \chi^2_1$ are central.
The variables $Q$ and $R$ are independent if and only if $Var(X) = Var(Y)$.
In general, $Q$ and $R$ are noncentral and dependent.
Solution 2:
As @Yemon Choi showed in the first question, without any hypothesis the answer is negative since $P(X^2<0)=0$ whereas $P(U<0)\neq 0$ if $U$ is Gaussian.
For the second question the answer is also no. Take $X$ and $Y$ two Gaussian random variables with mean $0$ and variance $1$. Since they have the same variance, $X-Y$ and $X+Y$ are independent Gaussian random variables. Put $Z:=\frac{X^2-Y^2}2=\frac{X-Y}{\sqrt 2}\frac{X+Y}{\sqrt 2}$. Then $Z$ is the product of two independent Gaussian, but the characteristic function of $Z$ is $\varphi_Z(t)=\frac 1{\sqrt{1+t^2}}$, which is not the characteristic function of a Gaussian.
Solution 3:
You can use moments to see that the product $XY$ of independent normals cannot be normal except in trivial cases. By trivial, I mean $\mathbb{V}(X)\mathbb{V}(Y)=0.$
Suppose that $X,Y$ are independent normals so that $XY$ is normal.
Case 1: Suppose that $\mathbb{E}(X)=0$. By independence, $\mathbb{E}(XY)=\mathbb{E}(X)\mathbb{E}(Y)=0$, so $XY$ is mean zero normal, and hence
$$\mathbb{E}((XY)^4)=3\mathbb{E}((XY)^2)^2.$$ By independence we get
$$\mathbb{E}(X^4)\mathbb{E}(Y^4)=3\mathbb{E}(X^2)^2\mathbb{E}(Y^2)^2.$$
Either $\mathbb{V}(X)=0$, or dividing by $\mathbb{E}(X^4)$ gives $ \mathbb{E}(Y^4)= \mathbb{E}(Y^2)^2.$
This shows that $Y^2$, and hence $Y$, has zero variance.
Case 2: Suppose that $\mathbb{E}(X^2)>0$ and $\mathbb{E}(Y^2)>0$. Then, without loss of generality, $\mathbb{E}(X^2)=1$ and $\mathbb{E}(Y^2)=1$. In this case, we also have $\mathbb{E}((XY)^2)=1$, so $$\begin{eqnarray} \mathbb{E}(X^3)&=&\mathbb{E}(X)(3-2\mathbb{E}(X)^2)\\ \mathbb{E}(Y^3)&=&\mathbb{E}(Y)(3-2\mathbb{E}(Y)^2)\\ \mathbb{E}((XY)^3)&=&\mathbb{E}(XY)(3-2\mathbb{E}(XY)^2) \end{eqnarray} $$ Subtracting the product of the first two lines from the third line gives $$0=6\mathbb{E}(X)\,\mathbb{E}(Y)\,\mathbb{V}(X)\,\mathbb{V}(Y).$$ Either we are in a trivial case, or back to Case 1.
Thus, the product cannot be normal except in trivial cases.