Separable resolution-of-the-identity with all-electron Gaussian bases: Application to cubic-scaling RPA.

We explore a separable resolution-of-the-identity (RI) formalism built on quadratures over limited sets of real-space points designed for all-electron calculations. Our implementation preserves, in particular, the use of common atomic orbitals and their related auxiliary basis sets. The setup of the present density fitting scheme, i.e., the calculation of the system specific quadrature weights, scales cubically with respect to the system size. Extensive accuracy tests are presented for the Fock exchange and MP2 correlation energies. We finally demonstrate random phase approximation (RPA) correlation energy calculations with a scaling that is cubic in terms of operations, quadratic in memory, with a small crossover with respect to our standard RI-RPA implementation.

[1]  David Pines,et al.  Correlation Energy of a Free Electron Gas , 1958 .

[2]  Georg Kresse,et al.  Low Scaling Algorithms for the Random Phase Approximation: Imaginary Time and Laplace Transformations. , 2014, Journal of chemical theory and computation.

[3]  D. Pines.  A Collective Description of-Electron Interactions: III. Coulomb Interactions in a Degenerate Electron Gas , 1953 .

[4]  Wim Klopper,et al.  Implementation of the Bethe−Salpeter equation in the TURBOMOLE program , 2017, J. Comput. Chem..

[5]  G. Cuniberti,et al.  Decacene: On-Surface Generation. , 2017, Angewandte Chemie.

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

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

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

[9]  Samia M. Hamed,et al.  Evaluating the GW Approximation with CCSD(T) for Charged Excitations Across the Oligoacenes. , 2016, Journal of chemical theory and computation.

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

[11]  Denis Jacquemin,et al.  Benchmarking the Bethe–Salpeter Formalism on a Standard Organic Molecular Set , 2015, Journal of chemical theory and computation.

[12]  Lexing Ying,et al.  Fast algorithm for periodic density fitting for Bloch waves , 2015, 1512.00432.

[13]  Carlo Adamo,et al.  Extensive TD-DFT Benchmark: Singlet-Excited States of Organic Molecules. , 2009, Journal of chemical theory and computation.

[14]  Fabien Bruneval,et al.  A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules. , 2015, The Journal of chemical physics.

[15]  R A Friesner,et al.  New methods for electronic structure calculations on large molecules. , 1991, Annual review of physical chemistry.

[16]  T. Van Voorhis,et al.  Fluctuation-dissipation theorem density-functional theory. , 2005, The Journal of chemical physics.

[17]  Denis Jacquemin,et al.  The Bethe-Salpeter equation in chemistry: relations with TD-DFT, applications and challenges. , 2018, Chemical Society reviews.

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

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

[20]  J. Mintmire,et al.  Fitting the Coulomb potential variationally in linear-combination-of-atomic-orbitals density-functional calculations , 1982 .

[21]  Christian Ochsenfeld,et al.  Communication: An effective linear-scaling atomic-orbital reformulation of the random-phase approximation using a contracted double-Laplace transformation. , 2016, The Journal of chemical physics.

[22]  Eran Rabani,et al.  Expeditious Stochastic Calculation of Random-Phase Approximation Energies for Thousands of Electrons in Three Dimensions. , 2013, The journal of physical chemistry letters.

[23]  Jianfeng Lu,et al.  Cubic scaling algorithms for RPA correlation using interpolative separable density fitting , 2017, J. Comput. Phys..

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

[25]  K. Burke,et al.  Describing static correlation in bond dissociation by Kohn-Sham density functional theory. , 2004, The Journal of chemical physics.

[26]  Christian Ochsenfeld,et al.  Vanishing-Overhead Linear-Scaling Random Phase Approximation by Cholesky Decomposition and an Attenuated Coulomb-Metric. , 2017, Journal of chemical theory and computation.

[27]  J. Moussa,et al.  Cubic-scaling algorithm and self-consistent field for the random-phase approximation with second-order screened exchange. , 2013, The Journal of chemical physics.

[28]  Tjerk P. Straatsma,et al.  NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations , 2010, Comput. Phys. Commun..

[29]  X. Blase,et al.  Combining the Many-Body GW Formalism with Classical Polarizable Models: Insights on the Electronic Structure of Molecular Solids. , 2016, The journal of physical chemistry letters.

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

[31]  C. Van Alsenoy,et al.  Ab initio calculations on large molecules: The multiplicative integral approximation , 1988 .

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

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

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

[35]  Lexing Ying,et al.  Compression of the electron repulsion integral tensor in tensor hypercontraction format with cubic scaling cost , 2014, J. Comput. Phys..

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

[37]  Chao Yang,et al.  molgw 1: Many-body perturbation theory software for atoms, molecules, and clusters , 2016, Comput. Phys. Commun..

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

[39]  R. Needs,et al.  Space-time method for ab initio calculations of self-energies and dielectric response functions of solids. , 1995, Physical review letters.

[40]  A. Tikhonov,et al.  Numerical Methods for the Solution of Ill-Posed Problems , 1995 .

[41]  W. Thiel,et al.  Basis set effects on coupled cluster benchmarks of electronically excited states: CC3, CCSDR(3) and CC2 , 2010 .

[42]  R. Ahlrichs,et al.  Efficient molecular numerical integration schemes , 1995 .

[43]  A. Tkatchenko,et al.  Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions , 2012, 1201.0655.

[44]  X. Blase,et al.  Hybrid and Constrained Resolution-of-Identity Techniques for Coulomb Integrals. , 2017, Journal of chemical theory and computation.

[45]  Mihály Kállay,et al.  Linear-scaling implementation of the direct random-phase approximation. , 2015, The Journal of chemical physics.

[46]  Decay properties of the one-particle Green function in real space and imaginary time , 2000, cond-mat/0008399.

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

[48]  V. Lebedev.  Values of the nodes and weights of ninth to seventeenth order gauss-markov quadrature formulae invariant under the octahedron group with inversion☆ , 1975 .

[49]  Jan Almlöf,et al.  Elimination of energy denominators in Møller—Plesset perturbation theory by a Laplace transform approach , 1991 .

[50]  B. Lundqvist,et al.  Exchange and correlation in atoms, molecules, and solids by the spin-density-functional formalism , 1976 .

[51]  John P. Perdew,et al.  Exchange-correlation energy of a metallic surface: Wave-vector analysis , 1977 .

[52]  Matthias Scheffler,et al.  Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory , 2015 .

[53]  F. Furche.  Developing the random phase approximation into a practical post-Kohn-Sham correlation model. , 2008, The Journal of chemical physics.

[54]  D. Pines.  A Collective Description of Electron Interactions: II. Collective vs Individual Particle Aspects of the Interactions , 1952 .

[55]  Marco Häser,et al.  Møller-Plesset (MP2) perturbation theory for large molecules , 1993 .

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

[57]  Walter Thiel,et al.  Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3. , 2008, The Journal of chemical physics.

[58]  Patrick Seewald,et al.  Large-Scale Cubic-Scaling Random Phase Approximation Correlation Energy Calculations Using a Gaussian Basis. , 2016, Journal of chemical theory and computation.