Let $X_i$ be independent normally distributed random variables with zero mean and variance $\sigma^2 \neq 1$. What is the probability density function of the random variable formed by the sum of their squares?
Here is my attempt: Let $Y = \sum_{i=1}^{k}X_i^2$. Then $$\begin{align}Y &= \sigma^2\sum_{i=1}^{k}(\frac{X_i}{\sigma^2})^2\\ &\sim \sigma^2\chi^2_{k}\\ &\sim \sigma^2\Gamma_{k/2}(\theta = 2)\\ &\sim \Gamma_{k/2}(\theta = 2\sigma^2) \end{align} $$
where $\chi^2_{k}$ is the central chi-squared distribution, $\Gamma_{k}(\theta)$ is the Gamma distribution with Scaling parameter $\theta$.
Therefore, the mean of $Y$ is $\frac{k}{2}2\sigma^2 = k \sigma^2$ and variance $\frac{k}{2}(2\sigma^2)^2 = 2k\sigma^4$.
Also, is it possible to derive the mean and variance without going back to the Gamma distribution?