Linear-scaling tight binding from a truncated-moment approach
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (United States)
We present an approximation to the total-energy tight-binding (TB) method based on use of the kernel polynomial method within a truncated subspace. Chebyshev polynomial moments of the Hamiltonian matrix are generated in a stable and efficient manner through recursive matrix-vector multiples. To compute the energy, either the electronic density of states (DOS) or the zero-temperature Fermi function is smeared by convolution with the kernel polynomial, with Jackson damping to minimize Gibbs oscillations while maintaining the positivity of the DOS. These are shown to give approximate lower and upper bounds, respectively, on the exact TB energy, and are averaged to obtain an improved energy estimate. The scaling of the computational work is made linear in the number of atoms by truncating the moment computation at a certain range about each atom. Energy derivatives necessary for molecular dynamics are obtained via a matrix-polynomial derivative relation. The method converges to exact TB as the number of moments and the truncation range are increased. We demonstrate the convergence properties and viability of the method for materials simulations in an examination of defects in silicon. We also discuss the relative importance of truncation range versus number of moments. {copyright} {ital 1996 The American Physical Society.}
- OSTI ID:
- 284354
- Journal Information:
- Physical Review, B: Condensed Matter, Vol. 53, Issue 19; Other Information: PBD: May 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
Calculation of densities of states and spectral functions by Chebyshev recursion and maximum entropy
A sparse matrix–vector multiplication based algorithm for accurate density matrix computations on systems of millions of atoms