Fast spectral projection algorithms for density-matrix computations

被引:39
作者
Beylkin, G [1 ]
Coult, N [1 ]
Mohlenkamp, MJ [1 ]
机构
[1] Univ Colorado, Dept Math Appl, Boulder, CO 80309 USA
关键词
spectral projectors; density matrix; fast algorithms; wavelets; partitioned SVD;
D O I
10.1006/jcph.1999.6215
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present a fast algorithm for the construction of a spectral projector. This algorithm allows us to compute the density matrix, as used in, e.g., the Kohn-Sham iteration. and so obtain the electron density. We compute the spectral projector by constructing the matrix sign function through a simple polynomial recursion. We present several matrix representations for fast computation within this recursion, using bases with controlled space-spatial-frequency localization. In particular we consider wavelet and local cosine bases. Since spectral projectors appear in many contexts, we expect many additional applications of our approach. (C) 1999 Academic Press.
引用
收藏
页码:32 / 54
页数:23
相关论文
共 32 条