An efficient and accurate decomposition of the Fermi operator.

We present a method to compute the Fermi function of the Hamiltonian for a system of independent fermions based on an exact decomposition of the grand-canonical potential. This scheme does not rely on the localization of the orbitals and is insensitive to ill-conditioned Hamiltonians. It lends itself naturally to linear scaling as soon as the sparsity of the system's density matrix is exploited. By using a combination of polynomial expansion and Newton-like iterative techniques, an arbitrarily large number of terms can be employed in the expansion, overcoming some of the difficulties encountered in previous papers. Moreover, this hybrid approach allows us to obtain a very favorable scaling of the computational cost with increasing inverse temperature, which makes the method competitive with other Fermi operator expansion techniques. After performing an in-depth theoretical analysis of computational cost and accuracy, we test our approach on the density functional theory Hamiltonian for the metallic phase of the LiAl alloy.

[1]  R. Baer,et al.  Fast methods for resumming matrix polynomials and Chebyshev matrix polynomials , 2004 .

[2]  Martin Head-Gordon,et al.  Chebyshev expansion methods for electronic structure calculations on large molecular systems , 1997 .

[3]  Alavi,et al.  Ab initio molecular dynamics with excited electrons. , 1994, Physical review letters.

[4]  Galli,et al.  Large scale electronic structure calculations. , 1992, Physical review letters.

[5]  Matthias Krack,et al.  Pseudopotentials for H to Kr optimized for gradient-corrected exchange-correlation functionals , 2005 .

[6]  Yihan Shao,et al.  Improved Fermi operator expansion methods for fast electronic structure calculations , 2003 .

[7]  Michele Parrinello,et al.  Linear scaling electronic structure calculations and accurate statistical mechanics sampling with noisy forces , 2006 .

[8]  P. Löwdin.  On the Non‐Orthogonality Problem Connected with the Use of Atomic Wave Functions in the Theory of Molecules and Crystals , 1950 .

[9]  S. Goedecker,et al.  Relativistic separable dual-space Gaussian pseudopotentials from H to Rn , 1998, cond-mat/9803286.

[10]  W. Kohn,et al.  Nearsightedness of electronic matter. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[11]  M. Teter,et al.  Tight-binding electronic-structure calculations and tight-binding molecular dynamics with localized orbitals. , 1994, Physical review. B, Condensed matter.

[12]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

[13]  P. Giannozzi,et al.  Towards Very Large-Scale Electronic-Structure Calculations , 1992 .

[14]  Matthias Krack,et al.  Efficient and accurate Car-Parrinello-like approach to Born-Oppenheimer molecular dynamics. , 2007, Physical review letters.

[15]  David E. Manolopoulos,et al.  Canonical purification of the density matrix in electronic-structure theory , 1998 .

[16]  A. Tustin.  Automatic Control , 1951, Nature.

[17]  D. Frenkel,et al.  Grand-canonical simulations of solvated ideal fermions. Evidence for phase separation , 1992 .

[18]  M. Parrinello,et al.  Linear scaling for quasi-one-dimensional systems , 2006 .

[19]  Victor Y. Pan,et al.  Efficient parallel solution of linear systems , 1985, STOC '85.

[20]  Li,et al.  Density-matrix electronic-structure method with linear system-size scaling. , 1993, Physical review. B, Condensed matter.

[21]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[22]  Colombo,et al.  Efficient linear scaling algorithm for tight-binding molecular dynamics. , 1994, Physical review letters.

[23]  Teter,et al.  Separable dual-space Gaussian pseudopotentials. , 1996, Physical review. B, Condensed matter.

[24]  Michele Parrinello,et al.  Stochastic linear scaling for metals and nonmetals , 2005 .

[25]  Michele Parrinello,et al.  Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach , 2005, Comput. Phys. Commun..

[26]  M. Parrinello,et al.  Linear scaling electronic structure Monte Carlo method for metals , 2007 .

[27]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[28]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[29]  C. Loan.  A note on the evaluation of matrix polynomials , 1979 .