When deriving the Hartee-Fock method, we minimize the electronic energy with respect to all molecular orbitals with the contraint of orthonormality of the molecular orbitals by using the method of Lagrance multipliers. Is there a fundamental reason why the molecular orbitals need to be orthogonal? Does it ensure a lower energy compared to any non-orthogonal set of molecular orbitals?
Thank you very much for your help