Convergent recursive O(N) calculations for ab initio tight-binding
A theory is presented for a novel recursion method for O(N) ab initio tight-binding calculations. A long-standing problem of generalizing the recursion method to a non-orthogonal basis, which is a crucial step to make the recursion method applicable to ab initio calculations, is solved by the introduction of the hybrid representation and the two-sided block Lanczos algorithm. As a test of efficiency of the new method, convergence properties in energy and force of insulators, semiconductors, metals, and molecules are studied for not only simple model systems but also some real materials within the density functional theory. The present O(N) method possesses good convergence properties for metals as well as insulators.