Convergent recursive O(N) method for ab initio tight-binding calculations

A theory is presented for a recursion method for $\mathrm{O}(N)$ ab initio tight-binding calculations based on the density-functional theory. 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 hybrid representation and a two-sided block Lanczos algorithm. As a test of efficiency of the proposed 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 $\mathrm{O}(N)$ method possesses good convergence properties for metals as well as insulators.