An efficient self-consistent field method for large systems of weakly interacting components.

An efficient method for removing the self-consistent field (SCF) diagonalization bottleneck is proposed for systems of weakly interacting components. The method is based on the equations of the locally projected SCF for molecular interactions (SCF MI) which utilize absolutely localized nonorthogonal molecular orbitals expanded in local subsets of the atomic basis set. A generalization of direct inversion in the iterative subspace for nonorthogonal molecular orbitals is formulated to increase the rate of convergence of the SCF MI equations. Single Roothaan step perturbative corrections are developed to improve the accuracy of the SCF MI energies. The resulting energies closely reproduce the conventional SCF energy. Extensive test calculations are performed on water clusters up to several hundred molecules. Compared to conventional SCF, speedups of the order of (N/O)2 have been achieved for the diagonalization step, where N is the size of the atomic orbital basis, and O is the number of occupied molecular orbitals.

[1]  Steve Scheiner,et al.  Hydrogen Bonding: A Theoretical Perspective , 1997 .

[2]  Martin Head-Gordon,et al.  A tensor formulation of many-electron theory in a nonorthogonal single-particle basis , 1998 .

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

[4]  T. Nagata,et al.  Perturbation expansion theory corrected from basis set superposition error. I. Locally projected excited orbitals and single excitations. , 2004, The Journal of chemical physics.

[5]  M. Tuckerman.  Ab initio molecular dynamics: basic concepts, current trends and novel applications , 2002 .

[6]  Jon Baker,et al.  Q‐Chem 2.0: a high‐performance ab initio electronic structure program package , 2000, J. Comput. Chem..

[7]  Jing Kong,et al.  Fast and accurate Coulomb calculation with Gaussian functions. , 2005, The Journal of chemical physics.

[8]  Hermann Stoll,et al.  On the use of local basis sets for localized molecular orbitals , 1980 .

[9]  Martin Head-Gordon,et al.  Quantum chemistry and molecular processes , 1996 .

[10]  G. Halász,et al.  BSSE-free SCF theories: a comment , 2000 .

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

[12]  J. Kinoshita.  Simple Life Satisfies This Grad Student , 1996, Science.

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

[14]  M. Raimondi,et al.  Extension of the SCF-MI Method to the Case of K Fragments one of which is an Open-Shell System. , 1998 .

[15]  Martin Head-Gordon,et al.  Approaching the Basis Set Limit in Density Functional Theory Calculations Using Dual Basis Sets without Diagonalization , 2004 .

[16]  Ralf Ludwig,et al.  Water: From Clusters to the Bulk. , 2001, Angewandte Chemie.

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

[18]  Mario Raimondi,et al.  Modification of the Roothaan equations to exclude BSSE from molecular interaction calculations , 1996 .

[19]  Christian Ochsenfeld,et al.  Locality and Sparsity of Ab Initio One-Particle Density Matrices and Localized Orbitals , 1998 .

[20]  G. A. Jeffrey,et al.  An Introduction to Hydrogen Bonding , 1997 .

[21]  J. Cullen,et al.  An examination of the effects of basis set and charge transfer in hydrogen-bonded dimers with a constrained Hartree–Fock method , 1991 .

[22]  M. Head‐Gordon,et al.  Curvy steps for density matrix based energy minimization: tensor formulation and toy applications , 2003 .

[23]  M. Raimondi,et al.  Application of the DIIS technique to improve the convergence properties of the SCF-MI algorithm , 2001 .

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

[25]  Yihan Shao,et al.  Efficient evaluation of the Coulomb force in density-functional theory calculations , 2001 .

[26]  A. Szabó,et al.  Modern quantum chemistry : introduction to advanced electronic structure theory , 1982 .

[27]  F. Keutsch,et al.  The water trimer. , 2003, Chemical reviews.

[28]  Osamu Takahashi,et al.  Basis set superposition error free self-consistent field method for molecular interaction in multi-component systems : Projection operator formalism , 2001 .

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

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

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

[32]  Martin Head-Gordon,et al.  Absolute and Relative Energies From Polarized Atomic Orbital Self-consistent Field Calculations and a Second Order Correction.: Convergence with Size and Composition of the Secondary Basis , 2000, Comput. Chem..

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

[34]  Martin Head-Gordon,et al.  A J matrix engine for density functional theory calculations , 1996 .

[35]  Ross D. Adamson,et al.  Empirical density functionals , 1998 .

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

[37]  A. Castleman,et al.  van der Waais Molecules II: Introduction , 1994 .

[38]  Mark E. Tuckerman,et al.  Ab Initio Molecular Dynamics Simulations , 1996 .

[39]  David R. Bowler,et al.  Recent progress in linear scaling ab initio electronic structure techniques , 2002 .

[40]  T. Dunning,et al.  AB initio characterization of water and anion-water clusters , 1998 .

[41]  R. Parr.  Density-functional theory of atoms and molecules , 1989 .

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