Convergence can be quite a lot faster if we have better starting guesses than the unscreened hydrogenic orbitals we have now. Either we can implement some heuristic method that counts the number of electrons "below" the orbital to be initialized, or we use the Roothan-Hartree-Fock wavefunctions as tabulated by