SIAM Digital Library
 
 
 

You are not logged in Logged Out Log In

SIAM J. Sci. Comput. 34, pp. B1-B23 (23 pages)

Controlling Errors in Recursive Fermi–Dirac Operator Expansions with Applications in Electronic Structure Theory

Emanuel H. Rubensson

Full Text: Download PDF | Buy PDF (US$25) | View Cart
Recursive Fermi–Dirac operator expansion is an efficient way to compute one-electron density matrices in electronic structure theory. The convergence is rapid and depends only weakly on the conditioning of the problem and, for many systems, the computational cost increases only linearly with system size. In this article, errors introduced when evaluating the recursive expansion are analyzed and schemes to control the forward error are proposed. The error has previously been analyzed for explicit schemes working at zero electronic temperature [J. Chem. Phys., 128 (2008), 074106]. Here, implicit schemes [Phys. Rev. B, 68 (2003), 233104] working at zero or finite temperature are treated. The proposed schemes for error control are demonstrated by tight-binding as well as density functional theory electronic structure calculations on several test systems. Condition numbers for the problem of computing the density matrix are derived, giving quantitative insight into under what circumstances a temperature dependent formulation results in better conditioning. It is shown that for the considered recursive expansions, the number of matrix-matrix multiplications needed to compute the density matrix increases only with the squared logarithm of the condition number of the problem.

© 2012 Society for Industrial and Applied Mathematics

RELATED DATABASES

To view database links for this article, you need to log in.

PUBLICATION DATA

ISSN

1064-8275 (print)  
1095-7197 (online)

ARTICLE DATA

History
Received May 10, 2011
Accepted September 30, 2011
Published online January 31, 2012

For access to fully linked references, you need to log in.

For access to citing articles, you need to log in.


Close

close