Kernel Polynomial Method for non-orthogonal electronic structure calculations

摘要

The Kernel Polynomial Method (KPM) has been successfully applied to tight-binding electronic structure calculations as an O(N) method. Here we extend this method to nonorthogonal basis sets with a sparse overlap matrix S and a sparse Hamiltonian H. Since the KPM method utilizes matrix vector multiplications it is necessary to apply S(sup -1) H onto a vector. The multiplication of S(sup -1) is performed using a preconditioned conjugate gradient method and does not involve the explicit inversion of S. Hence the method scales the same way as the original KPM method, i.e. O(N), although there is an overhead due to the additional conjugate gradient part. We show an application of this method to defects in a titanate/platinum interface and to a large scale electronic structure calculation of amorphous diamond.