Linear scaling electronic structure methods

Abstract

Methods exhibiting linear scaling with respect to the size of the system, the so-called O(N) methods, are an essential tool for the calculation of the electronic structure of large systems containing many atoms. They are based on algorithms that take advantage of the decay properties of the density matrix. In this article the physical decay properties of the density matrix will first be studied for both metals and insulators. Several strategies for constructing O(N) algorithms will then be presented and critically examined. Some issues that are relevant only for self-consistent O(N) methods, such as the calculation of the Hartree potential and mixing issues, will also be discussed. Finally some typical applications of O(N) methods are briefly described.