Cholesky decomposition of the two-electron integral matrix in electronic structure calculations.

A standard Cholesky decomposition of the two-electron integral matrix leads to integral tables which have a huge number of very small elements. By neglecting these small elements, it is demonstrated that the recursive part of the Cholesky algorithm is no longer a bottleneck in the procedure. It is shown that a very efficient algorithm can be constructed when family type basis sets are adopted. For subsequent calculations, it is argued that two-electron integrals represented by Cholesky integral tables have the same potential for simplifications as density fitting. Compared to density fitting, a Cholesky decomposition of the two-electron matrix is not subjected to the problem of defining an auxiliary basis for obtaining a fixed accuracy in a calculation since the accuracy simply derives from the choice of a threshold for the decomposition procedure. A particularly robust algorithm for solving the restricted Hartree-Fock (RHF) equations can be speeded up if one has access to an ordered set of integral tables. In a test calculation on a linear chain of beryllium atoms, the advocated RHF algorithm nicely converged, but where the standard direct inversion in iterative space method converged very slowly to an excited state.

[1]  Emily A. Carter,et al.  Pseudospectral full configuration interaction , 1992 .

[2]  Martin Schütz,et al.  Low-order scaling local electron correlation methods. III. Linear scaling local perturbative triples correction (T) , 2000 .

[3]  Richard A. Friesner,et al.  Pseudospectral localized Mo/ller–Plesset methods: Theory and calculation of conformational energies , 1995 .

[4]  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.

[5]  Michael W. Schmidt,et al.  Effective convergence to complete orbital bases and to the atomic Hartree–Fock limit through systematic sequences of Gaussian primitives , 1979 .

[6]  Martin Schütz,et al.  Low-order scaling local electron correlation methods. V. Connected triples beyond (T): Linear scaling local CCSDT-1b , 2002 .

[7]  Richard A. Friesner,et al.  Pseudospectral localized generalized Mo/ller–Plesset methods with a generalized valence bond reference wave function: Theory and calculation of conformational energies , 1997 .

[8]  E. Cancès,et al.  On the convergence of SCF algorithms for the Hartree-Fock equations , 2000 .

[9]  Steven A. Orszag,et al.  Comparison of Pseudospectral and Spectral Approximation , 1972 .

[10]  Claude Le Bris,et al.  Mathematical models and methods for ab initio quantum chemistry , 2000 .

[11]  Peter Pulay,et al.  Localizability of dynamic electron correlation , 1983 .

[12]  Alistair P. Rendell,et al.  COUPLED-CLUSTER THEORY EMPLOYING APPROXIMATE INTEGRALS : AN APPROACH TO AVOID THE INPUT/OUTPUT AND STORAGE BOTTLENECKS , 1994 .

[13]  Stephen Wilson,et al.  Universal basis sets and Cholesky decomposition of the two-electron integral matrix , 1990 .

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

[15]  N. H. Beebe,et al.  Simplifications in the generation and transformation of two‐electron integrals in molecular calculations , 1977 .

[16]  G. G. Hall.  The molecular orbital theory of chemical valency VIII. A method of calculating ionization potentials , 1951, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[17]  Hans-Joachim Werner,et al.  Local treatment of electron correlation in coupled cluster theory , 1996 .

[18]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[19]  Richard A. Friesner,et al.  Pseudospectral Hartree–Fock theory: Applications and algorithmic improvements , 1990 .

[20]  R. K. Nesbet,et al.  Electronic Correlation in Atoms and Molecules , 2007 .

[21]  E. Davidson,et al.  One- and two-electron integrals over cartesian gaussian functions , 1978 .

[22]  Hans-Joachim Werner,et al.  Local perturbative triples correction (T) with linear cost scaling , 2000 .

[23]  Hans-Joachim Werner,et al.  Low-order scaling local electron correlation methods. IV. Linear scaling local coupled-cluster (LCCSD) , 2001 .

[24]  Philippe Y. Ayala,et al.  Linear scaling second-order Moller–Plesset theory in the atomic orbital basis for large molecular systems , 1999 .

[25]  I. Røeggen.  Analytic functions for the three-body potential of the helium trimer. , 2007, The Journal of chemical physics.

[26]  I. R. eggen.  Derivation of an extended geminal model , 1983 .

[27]  Frederick R. Manby,et al.  Fast linear scaling second-order Møller-Plesset perturbation theory (MP2) using local and density fitting approximations , 2003 .

[28]  Clemens C. J. Roothaan,et al.  New Developments in Molecular Orbital Theory , 1951 .

[29]  Richard A. Friesner,et al.  Solution of the Hartree–Fock equations for polyatomic molecules by a pseudospectral method , 1987 .

[30]  Richard A. Friesner,et al.  Solution of self-consistent field electronic structure equations by a pseudospectral method , 1985 .

[31]  V. R. Saunders,et al.  A “Level–Shifting” method for converging closed shell Hartree–Fock wave functions , 1973 .

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

[33]  T. Martínez,et al.  Pseudospectral multireference single and double excitation configuration interaction , 1995 .

[34]  Martin Head-Gordon,et al.  Closely approximating second-order Mo/ller–Plesset perturbation theory with a local triatomics in molecules model , 2000 .

[35]  Martin W. Feyereisen,et al.  Use of approximate integrals in ab initio theory. An application in MP2 energy calculations , 1993 .

[36]  Peter Pulay,et al.  The local correlation treatment. II. Implementation and tests , 1988 .

[37]  R. Lindh,et al.  Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals. , 2007, The Journal of chemical physics.

[38]  Peter Pulay,et al.  Fourth‐order Mo/ller–Plessett perturbation theory in the local correlation treatment. I. Method , 1987 .

[39]  I. Røeggen.  An ab initio study of the fcc and hcp structures of helium. , 2006, The Journal of chemical physics.

[40]  Yixiang Cao,et al.  Correlated ab Initio Electronic Structure Calculations for Large Molecules , 1999 .

[41]  I. Røeggen,et al.  On the Beebe-Linderberg two-electron integral approximation , 1986 .

[42]  Peter Pulay,et al.  Local configuration interaction: An efficient approach for larger molecules , 1985 .

[43]  Thomas Bondo Pedersen,et al.  Reduced scaling in electronic structure calculations using Cholesky decompositions , 2003 .

[44]  O. Sǐnanoğlu,et al.  Many‐Electron Theory of Atoms, Molecules and Their Interactions , 2007 .

[45]  Richard A. Friesner,et al.  Solution of the Hartree–Fock equations by a pseudospectral method: Application to diatomic molecules , 1986 .

[46]  E. Carter,et al.  Local weak-pairs pseudospectral multireference configuration interaction , 2002 .

[47]  Eric Cancès,et al.  SCF algorithms for HF electronic calculations , 2000 .

[48]  J. Simons,et al.  Application of cholesky-like matrix decomposition methods to the evaluation of atomic orbital integrals and integral derivatives , 1989 .

[49]  I. Røeggen.  An Extended Group Function Model for Intermolecular Interactions , 2006 .