He aquí una solución elemental:
Dejemos que $v\in D(\sqrt{A})$ . Obsérvese que la función $f: t \mapsto \frac{\lVert \sqrt{A}(u_0 + tv) \rVert^2 }{\lVert u_0 + tv \rVert^2}$ se minimiza en $t=0$ por supuesto, por lo que $f'(0) = 0$ . Ahora ampliando todo se obtiene $$f'(t) = \frac{(2 \operatorname{Re } \langle \sqrt{A}u_0, \sqrt{A}v \rangle + 2t\lVert \sqrt{A} v \rVert^2)\lVert u_0 +tv \rVert^2 - (2 \operatorname{Re } \langle u_0, v\rangle + 2t \Vert v \rVert^2)\lVert \sqrt{A} (u_0 + tv) \rVert^2}{\lVert u_0 + tv \rVert^4}$$ y evaluando en $t = 0$ da $$0 = 2 \operatorname{Re } \langle \sqrt{A} u_0 , \sqrt{A} v \rangle - 2\operatorname{Re } \langle u_0, v \rangle \lVert \sqrt{A} u_0 \rVert ^2,$$ o de forma equivalente $$\operatorname{Re} \langle \sqrt{A}u_0, \sqrt{A}v\rangle = E \operatorname{Re }\langle u_0, v \rangle.$$ Haciendo lo mismo con $iv$ en lugar de $v$ obtenemos el mismo resultado para la parte imaginaria y juntos $$\langle \sqrt{A} u_0, \sqrt{A} v \rangle = \langle Eu_0, v \rangle.$$ Desde $A$ es autoadjunto, $D(A) = D(A^*)$ por lo que basta con demostrar $u_0 \in D(A^*).$ De hecho, tenemos $$\sup_{\substack{v\in D(A) \\ \lVert v \rVert = 1}} \lvert \langle u_0, Av \rangle \rvert \leq \sup_{\substack{v\in D(\sqrt{A}) \\ \lVert v \rVert = 1}} \lvert \langle u_0, Av\rangle \rvert = \sup_{\substack{v\in D(\sqrt{A}) \\ \lVert v \rVert = 1}} \lvert \langle \sqrt{A}u_0, \sqrt{A}v\rangle \rvert = \sup_{\substack{v\in D(\sqrt{A}) \\ \lVert v \rVert = 1}} \lvert \langle Eu_0, v \rangle \rvert = E < \infty,$$ que implica $u_0 \in D(A^*) = D(A)$ . Entonces tenemos para cada $v\in D(A)$ $$\langle Eu_0 , v \rangle = \langle \sqrt{A} u_0 , \sqrt{A} v \rangle = \langle A u_0 , v \rangle$$ y como $D(A)$ es denso $Eu_0 = Au_0$ , según se desee.