[guided]The goal is to compute the variance of position, so we need the second moment and then subtract the square of the first moment. Before expanding $X^2$, we check that this expression is defined on the coherent state. Since $D$ is a linear subspace and $a,a^*:D\to D$, the vector $a\psi+a^*\psi$ lies in $D$ for every $\psi\in D$, and therefore $X\psi=c_X(a\psi+a^*\psi)$ also lies in $D$. In particular $X^2\psi_\alpha$ is well-defined. Expanding $X^2$ on $D$ gives
\begin{align*}
X^2 = c_X^2(a+a^*)^2 = c_X^2(a^2+aa^*+a^*a+(a^*)^2).
\end{align*}
The coherent-state eigenvalue equation gives the pure $a$ term:
\begin{align*}
a^2\psi_\alpha = a(\alpha\psi_\alpha)=\alpha^2\psi_\alpha.
\end{align*}
Because $\|\psi_\alpha\|_H=1$, this implies
\begin{align*}
(a^2\psi_\alpha,\psi_\alpha)_H = \alpha^2.
\end{align*}
For the pure $a^*$ term, we must use the adjoint relation carefully because $a$ and $a^*$ are unbounded operators. The theorem supplies the needed domain information: $\psi_\alpha \in D$, and $D$ is invariant under both $a$ and $a^*$. Therefore $a\psi_\alpha$, $a^*\psi_\alpha$, $a^2\psi_\alpha$, and $(a^*)^2\psi_\alpha$ are all defined. Applying the adjoint relation first to $a\psi_\alpha$ and $\psi_\alpha$, and then to $\psi_\alpha$ and $a^*\psi_\alpha$, gives
\begin{align*}
(a^2\psi_\alpha,\psi_\alpha)_H=(a\psi_\alpha,a^*\psi_\alpha)_H=(\psi_\alpha,(a^*)^2\psi_\alpha)_H.
\end{align*}
Taking complex conjugates then gives
\begin{align*}
((a^*)^2\psi_\alpha,\psi_\alpha)_H = \overline{(\psi_\alpha,(a^*)^2\psi_\alpha)_H} = \overline{(a^2\psi_\alpha,\psi_\alpha)_H} = \bar{\alpha}^2.
\end{align*}
For the normally ordered mixed term, the adjoint relation first gives
\begin{align*}
(\psi_\alpha,a^*a\psi_\alpha)_H = (a\psi_\alpha,a\psi_\alpha)_H = (\alpha\psi_\alpha,\alpha\psi_\alpha)_H = |\alpha|^2.
\end{align*}
Because the inner product is linear in the first argument, this is the adjoint identity in the correct direction. To obtain the expectation with $a^*a\psi_\alpha$ in the first argument, we use conjugate symmetry and the fact that $|\alpha|^2$ is real:
\begin{align*}
(a^*a\psi_\alpha,\psi_\alpha)_H = \overline{(\psi_\alpha,a^*a\psi_\alpha)_H} = |\alpha|^2.
\end{align*}
The remaining mixed term is not normally ordered. Using the commutation relation $aa^*=a^*a+I$ and $\|\psi_\alpha\|_H=1$, we get
\begin{align*}
(aa^*\psi_\alpha,\psi_\alpha)_H = (a^*a\psi_\alpha,\psi_\alpha)_H+(\psi_\alpha,\psi_\alpha)_H = |\alpha|^2+1.
\end{align*}
Substituting these four scalar expectations into the expansion of $X^2$ gives
\begin{align*}
\mathbb E_{\psi_\alpha}[X^2] = c_X^2\bigl(\alpha^2+\bar{\alpha}^2+2|\alpha|^2+1\bigr).
\end{align*}
The square of the first moment is
\begin{align*}
\mathbb E_{\psi_\alpha}[X]^2 = c_X^2(\alpha+\bar{\alpha})^2 = c_X^2\bigl(\alpha^2+\bar{\alpha}^2+2|\alpha|^2\bigr).
\end{align*}
Subtracting this from the second moment cancels every term depending on $\alpha$:
\begin{align*}
(\Delta_{\psi_\alpha}X)^2 = \mathbb E_{\psi_\alpha}[X^2]-\mathbb E_{\psi_\alpha}[X]^2 = c_X^2 = \frac{\hbar}{2m\omega}.
\end{align*}[/guided]