Attractive electron–electron interactions within robust local fitting approximations

An analysis of Dunlap's robust fitting approach reveals that the resulting two‐electron integral matrix is not manifestly positive semidefinite when local fitting domains or non‐Coulomb fitting metrics are used. We present a highly local approximate method for evaluating four‐center two‐electron integrals based on the resolution‐of‐the‐identity (RI) approximation and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair‐atomic resolution‐of‐the‐identity (PARI) approach, atomic‐orbital (AO) products are expanded in auxiliary functions centered on the two atoms associated with each product. Numerical tests indicate that in 1% or less of all Hartree–Fock and Kohn–Sham calculations, the indefinite integral matrix causes nonconvergence in the self‐consistent‐field iterations. In these cases, the two‐electron contribution to the total energy becomes negative, meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the standard Coulomb‐metric RI method. The speed‐up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple‐zeta quality basis set. A positive semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky‐decomposition techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calculations. © 2013 Wiley Periodicals, Inc.

[1]  FRANCESCO AQUILANTE,et al.  MOLCAS 7: The Next Generation , 2010, J. Comput. Chem..

[2]  Frederick R. Manby,et al.  Fast Hartree–Fock theory using local density fitting approximations , 2004 .

[3]  James A Platts,et al.  Auxiliary basis sets for density fitting-MP2 calculations: Nonrelativistic triple-zeta all-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2008, The Journal of chemical physics.

[4]  G. Friesecke,et al.  Decay behavior of least-squares coefficients in auxiliary basis expansions. , 2005, The Journal of chemical physics.

[5]  F. Weigend.  Accurate Coulomb-fitting basis sets for H to Rn. , 2006, Physical chemistry chemical physics : PCCP.

[6]  Martin Head-Gordon,et al.  Scaled opposite-spin second order Møller-Plesset correlation energy: an economical electronic structure method. , 2004, The Journal of chemical physics.

[7]  Florian Weigend,et al.  A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency , 2002 .

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

[9]  Florian Weigend,et al.  Hartree–Fock exchange fitting basis sets for H to Rn † , 2008, J. Comput. Chem..

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

[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]  F. Manby,et al.  Explicitly correlated second-order perturbation theory using density fitting and local approximations. , 2006, The Journal of chemical physics.

[13]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials , 1995 .

[14]  R. Rein,et al.  Integral approximations for molecular orbital theory , 1966 .

[15]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[16]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[17]  Poul Jørgensen,et al.  The second-order approximate coupled cluster singles and doubles model CC2 , 1995 .

[18]  J. Connolly,et al.  On first‐row diatomic molecules and local density models , 1979 .

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

[20]  Roland Lindh,et al.  Unbiased auxiliary basis sets for accurate two-electron integral approximations. , 2007, The Journal of chemical physics.

[21]  Frederick R. Manby,et al.  Linear scaling local coupled cluster theory with density fitting. Part I: 4-external integrals , 2003 .

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

[23]  Martin Head-Gordon,et al.  Hartree-Fock exchange computed using the atomic resolution of the identity approximation. , 2008, The Journal of chemical physics.

[24]  Notker Rösch,et al.  Variational fitting methods for electronic structure calculations , 2010 .

[25]  E. V. Lenthe,et al.  Hybrid density functional calculations of nuclear magnetic shieldings using Slater-type orbitals and the zeroth-order regular approximation , 2009 .

[26]  Julian Yarkony,et al.  Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration. , 2010, The Journal of chemical physics.

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

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

[29]  Trygve Helgaker,et al.  Excitation energies in density functional theory: an evaluation and a diagnostic test. , 2008, The Journal of chemical physics.

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

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

[32]  Frank P. Billingsley,et al.  Limited Expansion of Diatomic Overlap (LEDO): A Near‐Accurate Approximate Ab Initio LCAO MO Method. I. Theory and Preliminary Investigations , 1971 .

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

[34]  Brett I. Dunlap,et al.  Robust and variational fitting , 2000 .

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

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

[37]  D. Foerster.  Elimination, in electronic structure calculations, of redundant orbital products. , 2008, The Journal of chemical physics.

[38]  Arnim Hellweg,et al.  Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to Rn , 2007 .

[39]  Evert Jan Baerends,et al.  Self-consistent molecular Hartree—Fock—Slater calculations I. The computational procedure , 1973 .

[40]  B. Dunlap,et al.  Robust and variational fitting: Removing the four-center integrals from center stage in quantum chemistry , 2000 .

[41]  D. York,et al.  Density-functional expansion methods: generalization of the auxiliary basis. , 2011, The Journal of chemical physics.

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

[43]  J. V. Ortiz,et al.  Integral approximations in ab initio, electron propagator calculations. , 2009, The Journal of chemical physics.

[44]  Roland Lindh,et al.  Density fitting with auxiliary basis sets from Cholesky decompositions , 2009 .

[45]  Christof Hättig,et al.  CC2 excitation energy calculations on large molecules using the resolution of the identity approximation , 2000 .

[46]  Roland Lindh,et al.  Analytic derivatives for the Cholesky representation of the two-electron integrals. , 2008, The Journal of chemical physics.

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

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

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

[50]  Roland Lindh,et al.  Atomic Cholesky decompositions: a route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency. , 2009, The Journal of chemical physics.

[51]  L. Cederbaum,et al.  On the Cholesky decomposition for electron propagator methods: General aspects and application on C(60). , 2009, The Journal of chemical physics.

[52]  Paweł Sałek,et al.  Linear-scaling implementation of molecular electronic self-consistent field theory. , 2007, The Journal of chemical physics.

[53]  F. Weigend,et al.  Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations , 2002 .

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

[55]  Ranbir Singh,et al.  J. Mol. Struct. (Theochem) , 1996 .

[56]  Marco Häser,et al.  CALCULATION OF EXCITATION ENERGIES WITHIN TIME-DEPENDENT DENSITY FUNCTIONAL THEORY USING AUXILIARY BASIS SET EXPANSIONS , 1997 .

[57]  B. I. Dunlap,et al.  Robust variational fitting: Gáspár's variational exchange can accurately be treated analytically , 2000 .

[58]  Ronald H. Felton,et al.  A new computational approach to Slater’s SCF–Xα equation , 1975 .

[59]  Patrizia Calaminici,et al.  Robust and efficient density fitting. , 2009, The Journal of chemical physics.

[60]  Christof Hättig,et al.  Optimization of auxiliary basis sets for RI-MP2 and RI-CC2 calculations: Core–valence and quintuple-ζ basis sets for H to Ar and QZVPP basis sets for Li to Kr , 2005 .

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

[62]  Holger Patzelt,et al.  RI-MP2: optimized auxiliary basis sets and demonstration of efficiency , 1998 .

[63]  T. Helgaker,et al.  Variational and robust density fitting of four-center two-electron integrals in local metrics. , 2008, The Journal of chemical physics.