Say I have $X \sim \mathcal N(a, b)$ and $Y\sim \mathcal N(c, d)$.  Is $XY$ also normally distributed?

Is the answer any different if we know that $X$ and $Y$ are independent?

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.
