Calculation of densities of states and spectral functions by Chebyshev recursion and maximum entropy
Abstract
We present an efficient algorithm for calculating spectral properties of large sparse Hamiltonian matrices such as densities of states and spectral functions. The combination of Chebyshev recursion and maximum entropy achieves highenergy resolution without significant roundoff error, machine precision, or numerical instability limitations. If controlled statistical or systematic errors are acceptable, CPU and memory requirements scale linearly in the number of states. The inference of spectral properties from moments is much better conditioned for Chebyshev moments than for power moments. We adapt concepts from the kernel polynomial method, a linear Chebyshev approximation with optimized Gibbs damping, to control the accuracy of Fourier integrals of positive nonanalytic functions. We compare the performance of kernel polynomial and maximum entropy algorithms for an electronic structure example.
 Publication:

Physical Review E
 Pub Date:
 October 1997
 DOI:
 10.1103/PhysRevE.56.4822
 arXiv:
 arXiv:condmat/9703229
 Bibcode:
 1997PhRvE..56.4822S
 Keywords:

 02.70.c;
 71.15.m;
 71.20.b;
 Computational techniques;
 simulations;
 Methods of electronic structure calculations;
 Electron density of states and band structure of crystalline solids;
 Condensed Matter
 EPrint:
 8 pages RevTex, 3 postscript figures