Let $\epsilon_n$ denote a real-valued discrete-time stochastic process of residuals, the ARCH($p$) specification is given by
\begin{equation} \label{1.1} \epsilon_n=Z_n\sqrt{\sigma_n} \end{equation} \begin{equation}\label{1.2} \sigma_n=\alpha_0+\sum \limits_{i=1}^p\alpha_i\epsilon_{n-i}^2\,, \end{equation}
where $\alpha_0, \alpha_1, \dots,\alpha_p$ are scalar parameters to be estimated, $\mu_n$ is the fitted model. $Z_n$, are a sequence of independent, identically distributed random variables with mean zero and variance one.
Now, can we say that $E(Z_n^2\sigma_n)=E(\sigma_n)$, since $E(Z_n^2)=1$?