Robust fitting techniques in the chain of spheres approximation to the Fock exchange: The role of the complementary space.

The chain of spheres seminumerical integration algorithm is reexamined from the point of view of various numerical error cancellation schemes. The two sources of ideas are the robust fitting technique of Dunlap and the dealiasing concept of the pseudospectral methods of Friesner. The most efficient scheme is identified and is found to be a generalization of the overlap fitting procedure discussed earlier. It will be shown that due to the conceptual similarity between complementary auxiliary basis functions of explicitly correlated approaches and the dealiasing functions of pseudospectral methods, the former can be used in a similar manner as the latter to cancel numerical errors. In contrast to the dealiasing functions, the complementary auxiliary basis functions are not specific to the integration grid. The scheme obtained can dramatically improve reaction energies even using small grid setups at a moderate additional cost.

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

[2]  Frank Neese,et al.  Efficient Structure Optimization with Second-Order Many-Body Perturbation Theory: The RIJCOSX-MP2 Method. , 2010, Journal of chemical theory and computation.

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

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

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

[6]  Frank Neese,et al.  An overlap fitted chain of spheres exchange method. , 2011, The Journal of chemical physics.

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

[8]  Automated selection of optimal Gaussian fits to arbitrary functions in electronic structure theory , 1998 .

[9]  Edward F. Valeev.  Improving on the resolution of the identity in linear R12 ab initio theories , 2004 .

[10]  Stefan Grimme,et al.  A General Database for Main Group Thermochemistry, Kinetics, and Noncovalent Interactions - Assessment of Common and Reparameterized (meta-)GGA Density Functionals. , 2010, Journal of chemical theory and computation.

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

[12]  William L. Jorgensen,et al.  PDDG/PM3 and PDDG/MNDO: Improved semiempirical methods , 2002, J. Comput. Chem..

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

[14]  Wim Klopper,et al.  Explicitly correlated second-order Møller–Plesset methods with auxiliary basis sets , 2002 .

[15]  Robert M Parrish,et al.  Communication: Tensor hypercontraction. III. Least-squares tensor hypercontraction for the determination of correlated wavefunctions. , 2012, The Journal of chemical physics.

[16]  Robert M Parrish,et al.  Quartic scaling second-order approximate coupled cluster singles and doubles via tensor hypercontraction: THC-CC2. , 2013, The Journal of chemical physics.

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

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

[19]  Robert M Parrish,et al.  Tensor hypercontraction. II. Least-squares renormalization. , 2012, The Journal of chemical physics.

[20]  Matthias Krack,et al.  AN ADAPTIVE NUMERICAL INTEGRATOR FOR MOLECULAR INTEGRALS , 1998 .

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

[22]  Benny G. Johnson,et al.  A standard grid for density functional calculations , 1993 .

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

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

[25]  Frank Neese,et al.  An improvement of the resolution of the identity approximation for the formation of the Coulomb matrix , 2003, J. Comput. Chem..

[26]  Hans-Joachim Werner,et al.  Systematically convergent basis sets for explicitly correlated wavefunctions: the atoms H, He, B-Ne, and Al-Ar. , 2008, The Journal of chemical physics.

[27]  F. Neese,et al.  Efficient time-dependent density functional theory approximations for hybrid density functionals: analytical gradients and parallelization. , 2011, The Journal of chemical physics.

[28]  Frank Neese,et al.  The ORCA program system , 2012 .

[29]  B. Dunlap.  ACCURATE DENSITY FUNCTIONAL CALCULATIONS ON LARGE SYSTEMS , 1996 .

[30]  F. Neese,et al.  Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange , 2009 .

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

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

[33]  Robert M Parrish,et al.  Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory. , 2012, The Journal of chemical physics.

[34]  Emily A. Carter,et al.  PSEUDOSPECTRAL METHODS APPLIED TO THE ELECTRON CORRELATION PROBLEM , 1995 .

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

[36]  Wim Klopper,et al.  Wave functions with terms linear in the interelectronic coordinates to take care of the correlation cusp. I. General theory , 1991 .

[37]  Robert M Parrish,et al.  Discrete variable representation in electronic structure theory: quadrature grids for least-squares tensor hypercontraction. , 2013, The Journal of chemical physics.

[38]  F. Neese,et al.  Speeding up spin-component-scaled third-order pertubation theory with the chain of spheres approximation: the COSX-SCS-MP3 method , 2013 .

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

[40]  W. Klopper,et al.  A composite ‘density fitting + numerical integration’ approximation for electron-repulsion integrals , 2013 .

[41]  J. Whitten,et al.  Electron repulsion integral approximations and error bounds: Molecular applications , 1974 .

[42]  F. Weigend,et al.  Gaussian basis sets of quadruple zeta valence quality for atoms H–Kr , 2003 .

[43]  F. Neese,et al.  The resolution of identity and chain of spheres approximations for the LPNO-CCSD singles Fock term , 2012 .