Let
- $d\in\mathbb N$
- $\Lambda\subseteq\mathbb R^d$ be bounded and open
- $\mathcal D(A):=\left\{u\in H_0^1(\Lambda):u\text{ admits a weak Laplacian }\Delta u\in L^2(\Lambda)\right\}$ and $$A:=-\Delta u\;\;\;\text{for }u\in\mathcal D(A)$$
$A$ is a densely-defined, positive definite and symmetric operator on $L^2(\Lambda)$.There is an ortonormal basis $(e_n)_{n\in\mathbb N}\subseteq\mathcal D(A)$ of $L^2(\Lambda)$ (which is a complete orthogonal system in $H_0^1(\Lambda)$) with $$Ae_n=\lambda_ne_n\tag1$$ for some $(\lambda_n)_{n\in\mathbb N}\subseteq(0,\infty)$ with $$\lambda_n\le\lambda_{n+1}\;\;\;\text{for all }n\in\mathbb N\tag2\;.$$
Let $$\mathcal D(A^\alpha):=\left\{u\in L^2(\Lambda):\sum_{n\in\mathbb N}\lambda_n^{2\alpha}\langle u,e_n\rangle_{L^2(\Lambda)}^2<\infty\right\}$$ for $\alpha\in\mathbb R$.
Question: How can we show that $\mathcal D(A^{1/2})=H_0^1(\Lambda)$?
Note that $$\left|u\right|:=\left\|\nabla u\right\|_{L^2(\Lambda)}\;\;\;\text{for }u\in H_0^1(\Lambda)$$ is a norm equivalent to $\left\|\;\cdot\;\right\|_{H^1(\Lambda)}$ on $H_0^1(\Lambda)$. Moreover, $$\sum_{n\in\mathbb N}\lambda_n\langle u,e_n\rangle_{L^2(\Lambda)}^2=\sum_{n\in\mathbb N}\langle u,Ae_n\rangle_{L^2(\Lambda)}\langle u,e_n\rangle_{L^2(\Lambda)}\stackrel{(\ast)}=\sum_{n\in\mathbb N}\langle Au,e_n\rangle_{L^2(\Lambda)}\langle u,e_n\rangle_{L^2(\Lambda)}=\langle Au,u\rangle_{L^2(\Lambda)}=\left|u\right|^2\tag 3$$ for all $u\in\mathcal D(A)$.
The problematic thing is that $(\ast)$ doesn't hold, unless $u\in\mathcal D(A)\subseteq\mathcal D(A^{1/2})$. Moreover, I don't see why $u\in\mathcal D(A^{1/2})\setminus\mathcal D(A)$ should admit a weak gradient; by definition we only know that $u\in L^2(\Lambda)$ with $\sum_{n\in\mathbb N}\lambda_n\langle u,e_n\rangle_{L^2(\Lambda)}^2<\infty$. So, what am I missing?