Incomplete selected inversion for linear-scaling electronic structure calculations
Pole Expansion and Selected Inversion (PEXSI) is an efficient scheme for evaluating selected entries of functions of large sparse matrices as required e.g. in electronic structure algorithms. We show that the triangular factorizations computed by the PEXSI scheme exhibit a localization property similar to that of matrix functions, and we present a modified PEXSI algorithm which exploits this observation to achieve linear scaling. To the best of our knowledge, the resulting incomplete PEXSI (iPEXSI) algorithm is the first linear-scaling algorithm which scales provably better than cubically even in the absence of localization, and we hope that this will help to further lower the critical system size where linear-scaling algorithms begin to outperform the diagonalization algorithm.
READ FULL TEXT