I'm having trouble with a few concepts dealing with slater determinants and many electron systems. In particular, when dealing with 3 electron systems, I know that
$$ \psi_0 = \frac{1}{\sqrt{6}} \det \begin{bmatrix} 1s(1)\alpha(1) & 1s(2)\alpha(2) & 1s(3)\alpha(3) \\ 1s(1)\beta(1) & 1s(2)\beta(2) & 1s(3)\beta(3) \\ 2s(1)\alpha(1) & 2s(2)\alpha(2) & 2s(3)\alpha(3) \end{bmatrix}, $$
where the $1s$ and $2s$ are spatial functions and $\alpha$, $\beta$ are spin functions.
Performing the determinant expansion yields 6 terms, and it turns out that $\psi_0$ is an eigenfunction for the unperturbed Hamiltonian (in atomic units) $$ \hat{H}_0 = -\frac{1}{2} \nabla^{2}_{1} -\frac{3}{r_1} -\frac{1}{2}\nabla^{2}_{2} -\frac{3}{r_2} -\frac{1}{2}\nabla^{2}_{3} -\frac{3}{r_3} $$ and that the resulting eigenenergy is $$E^{0}=E_{1s}+E_{1s}+E_{2s}.$$ How does one show this?
My confusion mainly stems from the difficulty in acting the Hamiltonian on the determinant expression above. For example, how does one compute (in bra-ket notation) $$ \frac{1}{\sqrt{6}} \hat{H}_0 | 1s(1)\alpha(1) 1s(2)\beta(2) 2s(3)\alpha(3)\rangle? $$ This is only the first term in the expansion, and once I have some intuition on how to do this properly, then the general pattern should emerge quite clearly. Any help and guidance would be greatly appreciated!