Linear-scaling implementation of molecular electronic self-consistent field theory.

A linear-scaling implementation of Hartree-Fock and Kohn-Sham self-consistent field (SCF) theories is presented and illustrated with applications to molecules consisting of more than 1000 atoms. The diagonalization bottleneck of traditional SCF methods is avoided by carrying out a minimization of the Roothaan-Hall (RH) energy function and solving the Newton equations using the preconditioned conjugate-gradient (PCG) method. For rapid PCG convergence, the Lowdin orthogonal atomic orbital basis is used. The resulting linear-scaling trust-region Roothaan-Hall (LS-TRRH) method works by the introduction of a level-shift parameter in the RH Newton equations. A great advantage of the LS-TRRH method is that the optimal level shift can be determined at no extra cost, ensuring fast and robust convergence of both the SCF iterations and the level-shifted Newton equations. For density averaging, the authors use the trust-region density-subspace minimization (TRDSM) method, which, unlike the traditional direct inversion in the iterative subspace (DIIS) scheme, is firmly based on the principle of energy minimization. When combined with a linear-scaling evaluation of the Fock/Kohn-Sham matrix (including a boxed fitting of the electron density), LS-TRRH and TRDSM methods constitute the linear-scaling trust-region SCF (LS-TRSCF) method. The LS-TRSCF method compares favorably with the traditional SCF/DIIS scheme, converging smoothly and reliably in cases where the latter method fails. In one case where the LS-TRSCF method converges smoothly to a minimum, the SCF/DIIS method converges to a saddle point.

[1]  Weitao Yang,et al.  A density‐matrix divide‐and‐conquer approach for electronic structure calculations of large molecules , 1995 .

[2]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[3]  Byron H. Lengsfield,et al.  General second order MCSCF theory: A density matrix directed algorithm , 1980 .

[4]  Eric Schwegler,et al.  Linear scaling computation of the Fock matrix , 1997 .

[5]  R. Baer,et al.  Electronic structure of large systems: Coping with small gaps using the energy renormalization group method , 1998 .

[6]  Weitao Yang,et al.  An algorithm for 3D numerical integration that scales linearly with the size of the molecule , 1995 .

[7]  Paweł Sałek,et al.  The trust-region self-consistent field method: towards a black-box optimization in Hartree-Fock and Kohn-Sham theories. , 2004, The Journal of chemical physics.

[8]  Matt Challacombe,et al.  A simplified density matrix minimization for linear scaling self-consistent field theory , 1999 .

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

[10]  M. Ratner.  Molecular electronic-structure theory , 2000 .

[11]  A. St.-Amant,et al.  Linear scaling for the charge density fitting procedure of the linear combination of Gaussian-type orbitals density functional method , 1996 .

[12]  B. C. Carlson,et al.  Orthogonalization Procedures and the Localization of Wannier Functions , 1957 .

[13]  Stefan Goedecker,et al.  Linear scaling electronic structure methods in chemistry and physics , 2003, Comput. Sci. Eng..

[14]  Martin Head-Gordon,et al.  Energy renormalization-group method for electronic structure of large systems , 1998 .

[15]  Florian Weigend,et al.  Auxiliary basis sets for main row atoms and transition metals and their use to approximate Coulomb potentials , 1997 .

[16]  Trygve Helgaker,et al.  Direct optimization of the AO density matrix in Hartree-Fock and Kohn-Sham theories , 2000 .

[17]  Poul Jo,et al.  A direct approach to second‐order MCSCF calculations using a norm extended optimization scheme , 1984 .

[18]  John R. Sabin,et al.  On some approximations in applications of Xα theory , 1979 .

[19]  Trygve Helgaker,et al.  Direct optimization of the atomic-orbital density matrix using the conjugate-gradient method with a multilevel preconditioner , 2001 .

[20]  Yihan Shao,et al.  An improved J matrix engine for density functional theory calculations , 2000 .

[21]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[22]  E. Hylleraas,et al.  Numerische Berechnung der 2S-Terme von Ortho- und Par-Helium , 1930 .

[23]  J. L. Whitten,et al.  Coulombic potential energy integrals and approximations , 1973 .

[24]  Martin Head-Gordon,et al.  Advances in methodologies for linear-scaling density functional calculations , 1996 .

[25]  Michael J. Frisch,et al.  Achieving linear scaling in exchange-correlation density functional quadratures , 1996 .

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

[27]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[28]  Vanderbilt,et al.  Generalization of the density-matrix method to a nonorthogonal basis. , 1994, Physical review. B, Condensed matter.

[29]  Yihan Shao,et al.  Curvy steps for density matrix-based energy minimization: Application to large-scale self-consistent-field calculations , 2003 .

[30]  Evert Jan Baerends,et al.  Towards an order , 1998 .

[31]  G. Scuseria,et al.  A black-box self-consistent field convergence algorithm: One step closer , 2002 .

[32]  J. MacDonald,et al.  Successive Approximations by the Rayleigh-Ritz Variation Method , 1933 .

[33]  Benny G. Johnson,et al.  THE CONTINUOUS FAST MULTIPOLE METHOD , 1994 .

[34]  Eric Schwegler,et al.  Linear scaling computation of the Hartree–Fock exchange matrix , 1996 .

[35]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[36]  Joseph E. Subotnik,et al.  Linear scaling density fitting. , 2006, The Journal of chemical physics.

[37]  Christian Ochsenfeld,et al.  Linear and sublinear scaling formation of Hartree-Fock-type exchange matrices , 1998 .

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

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

[40]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[41]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[42]  Eric Schwegler,et al.  Linear scaling computation of the Fock matrix. II. Rigorous bounds on exchange integrals and incremental Fock build , 1997 .

[43]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials (Chem. Phys. Letters 240 (1995) 283-290) , 1995 .

[44]  Isaiah Shavitt,et al.  Comparison of the convergence characteristics of some iterative wave function optimization methods , 1982 .

[45]  Eric Schwegler,et al.  Linear scaling computation of the Fock matrix. IV. Multipole accelerated formation of the exchange matrix , 1999 .

[46]  J. Almlöf,et al.  Integral approximations for LCAO-SCF calculations , 1993 .

[47]  E. Zaremba.  Some recent developments in density functional theory , 1983 .

[48]  Michael J. Frisch,et al.  Achieving Linear Scaling for the Electronic Quantum Coulomb Problem , 1996, Science.

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

[50]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .

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

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

[53]  R. Mcweeny,et al.  Methods Of Molecular Quantum Mechanics , 1969 .

[54]  Paweł Sałek,et al.  The trust-region self-consistent field method in Kohn-Sham density-functional theory. , 2005, The Journal of chemical physics.

[55]  Michael J. Frisch,et al.  A linear scaling method for Hartree–Fock exchange calculations of large molecules , 1996 .

[56]  Martin Head-Gordon,et al.  Auxiliary basis expansions for large-scale electronic structure calculations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.