I want to know how the analytical solutions of Hartree-Fock can come handy during computation of large systems although the formal scaling order of HF is N^4 compared to N^3 scaling of DFT (where N is an approximate measure of system size). I also want to know how it addresses the convergence problems that might arise with use of DFT for such large systems which can lead to an incorrect description of HOMO-LUMO band gaps and hence convergence issues.