Let $H$ be a Hilbert space over $\mathbb{C}$ and $A$ be a positive linear operator on $H$ (not necessarily bounded). Suppose that $A$'s eigenvectors $\{e_j\}$ form an ONB for $H$ and $\{c_j\}$s are the corresponding eigenvalues of $A$. We know from spectral theory that $A$'s square root $\sqrt{A}$ exists and is positive since $A$ is positive. But is it now true that $\{\sqrt{c_j}\}$s are eigenvalues of $\sqrt{A}$?
If I recall correctly, in the cast that $A$ is bounded we have that if $\sqrt{A}$ exists, it is self-adjoint. Therefore we know that for every $e_j$, $c_j = \left<Ae_j, e_j\right> = ||\sqrt{A}e_j||^2$ whence $\sqrt{c_j}\leq ||\sqrt{A}||$. But I don't really know how to turn this approach to a proof that $\sqrt{A}e_j = \sqrt{c_j}e_j$.
How about reasoning along the following lines. Let $v$ be an eigenvector of $A$ with eigenvalue $\lambda$, and let $$ w = \sqrt{A}v.$$ By looking at the matrix-representation of $\sqrt{A}|_{\operatorname{span}(v,w)}$ we easily find the eigenvectors $\sqrt{\lambda}v+w$ and $\sqrt{\lambda}-w$. It is easily check that $$\begin{align*} \sqrt{A}(\sqrt{\lambda}v+w) &= \sqrt{\lambda}\sqrt{A}v + Av \\ &= \sqrt{\lambda}w+\lambda v \\ &= \sqrt{\lambda}(\sqrt{\lambda}v + w) \end{align*}$$