If we have an infinite $L^2$ basis that is determined by the eigenfunctions of the eigenvalue problem $$ -\Delta \psi=\lambda \psi$$
How can we orthonormalize this basis with respect to the $H^1_0$ inner product $(\nabla \cdot, \nabla \cdot)$? Or why is this orthonormalization possible ?