Preconditioned iterative minimization for linear-scaling electronic structure calculations

Linear-scaling electronic structure methods are essential for calculations on large systems. Some of these approaches use a systematic basis set, the completeness of which may be tuned with an adjustable parameter similar to the energy cut-off of plane-wave techniques. The search for the electronic ground state in such methods suffers from an ill-conditioning which is related to the kinetic contribution to the total energy and which results in unacceptably slow convergence. We present a general preconditioning scheme to overcome this ill-conditioning and implement it within our own first-principles linear-scaling density functional theory method. The scheme may be applied in either real space or reciprocal space with equal success. The rate of convergence is improved by an order of magnitude and is found to be almost independent of the size of the basis.

[1]  R. Mcweeny.  Some Recent Advances in Density Matrix Theory , 1960 .

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

[3]  Peter D. Haynes,et al.  Corrected penalty-functional method for linear-scaling calculations within density-functional theory , 1999 .

[4]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[5]  C. Y. Fong,et al.  Real-space local polynomial basis for solid-state electronic-structure calculations: A finite-element approach , 1999, cond-mat/9903313.

[6]  David R. Bowler,et al.  Length-scale ill conditioning in linear-scaling DFT , 1998 .

[7]  B. Alder,et al.  THE GROUND STATE OF THE ELECTRON GAS BY A STOCHASTIC METHOD , 2010 .

[8]  J. D. Cloizeaux,et al.  Energy Bands and Projection Operators in a Crystal: Analytic and Asymptotic Properties , 1964 .

[9]  Hernández,et al.  Linear-scaling density-functional-theory technique: The density-matrix approach. , 1996, Physical review. B, Condensed matter.

[10]  J. Fattebert,et al.  Towards grid-based OÑNÖ density-functional theory methods: Optimized nonorthogonal orbitals and multigrid acceleration , 2000 .

[11]  Peter D. Haynes,et al.  Preconditioned conjugate gradient method for the sparse generalized eigenvalue problem in electronic structure calculations , 2001 .

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

[13]  Arash A. Mostofi,et al.  Nonorthogonal generalized Wannier function pseudopotential plane-wave method , 2002 .

[14]  Benny G. Johnson,et al.  Linear scaling density functional calculations via the continuous fast multipole method , 1996 .

[15]  A. Zunger,et al.  Self-interaction correction to density-functional approximations for many-electron systems , 1981 .

[16]  Leonard Kleinman,et al.  Efficacious Form for Model Pseudopotentials , 1982 .

[17]  Gustavo E. Scuseria,et al.  Linear scaling conjugate gradient density matrix search as an alternative to diagonalization for first principles electronic structure calculations , 1997 .

[18]  J. Pask,et al.  Finite-element methods in electronic-structure theory , 2001 .

[19]  D. Vanderbilt,et al.  Exponential decay properties of Wannier functions and related quantities. , 2001, Physical review letters.

[20]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[21]  M. Gillan.  Calculation of the vacancy formation energy in aluminium , 1989 .

[22]  Y. Saad,et al.  Finite-difference-pseudopotential method: Electronic structure calculations without a basis. , 1994, Physical review letters.

[23]  C. M. Goringe,et al.  Basis functions for linear-scaling first-principles calculations , 1997 .

[24]  Walter Kohn,et al.  Analytic Properties of Bloch Waves and Wannier Functions , 1959 .

[25]  Sullivan,et al.  Large-scale electronic-structure calculations with multigrid acceleration. , 1995, Physical review. B, Condensed matter.

[26]  Artacho,et al.  Nonorthogonal basis sets in quantum mechanics: Representations and second quantization. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[27]  Martin Head-Gordon,et al.  The tensor properties of energy gradients within a non-orthogonal basis , 1997 .

[28]  Gustavo E. Scuseria,et al.  Linear Scaling Density Functional Calculations with Gaussian Orbitals , 1999 .

[29]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .

[30]  O. Sankey,et al.  Ab initio multicenter tight-binding model for molecular-dynamics simulations and other applications in covalent systems. , 1989, Physical review. B, Condensed matter.

[31]  Sohrab Ismail-Beigi,et al.  LOCALITY OF THE DENSITY MATRIX IN METALS, SEMICONDUCTORS, AND INSULATORS , 1999 .

[32]  Matt Probert,et al.  First-principles simulation: ideas, illustrations and the CASTEP code , 2002 .

[33]  Nicola Marzari,et al.  Ensemble density-functional theory for ab initio molecular dynamics of metals and finite-temperature insulators , 1997 .

[34]  Arash A. Mostofi,et al.  Total-energy calculations on a real space grid with localized functions and a plane-wave basis , 2002 .

[35]  Giulia Galli,et al.  Linear scaling methods for electronic structure calculations and quantum molecular dynamics simulations , 1996 .

[36]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[37]  Gene H. Golub,et al.  Matrix computations , 1983 .

[38]  Peter D. Haynes,et al.  Localised spherical-wave basis set for O(N) total-energy pseudopotential calculations , 1997 .

[39]  R. O. Jones,et al.  The density functional formalism, its applications and prospects , 1989 .