A localized basis that allows fast and accurate second-order Moller-Plesset calculations

We present a method for computing a basis of localized orthonormal orbitals (both occupied and virtual), in whose representation the Fock matrix is extremely diagonal dominant. The existence of these orbitals is shown empirically to be sufficient for achieving highly accurate second-order Moller-Plesset (MP2) energies, calculated according to Kapuy's method. This method (which we abbreviate KMP2) involves a different partitioning of the n-electron Hamiltonian and scales at most quadratically, with potential for linearity, in the number of electrons. As such, we believe the KMP2 algorithm presented here could be the basis of a viable approach to local-correlation calculations.