Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals.

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in quantum chemistry and beyond.

[1]  Emanuel H. Rubensson,et al.  Bringing about matrix sparsity in linear‐scaling electronic structure calculations , 2011, J. Comput. Chem..

[2]  Martin Head-Gordon,et al.  A near linear-scaling smooth local coupled cluster algorithm for electronic structure. , 2006, The Journal of chemical physics.

[3]  Wilfried Meyer,et al.  PNO–CI Studies of electron correlation effects. I. Configuration expansion by means of nonorthogonal orbitals, and application to the ground state and ionized states of methane , 1973 .

[4]  G. Scuseria,et al.  Atom pair partitioning of the correlation energy , 2000 .

[5]  Georg Hetzer,et al.  Multipole approximation of distant pair energies in local MP2 calculations , 1998 .

[6]  Jörg Kussmann,et al.  Linear‐scaling self‐consistent field methods for large molecules , 2013 .

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

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

[9]  A. Becke A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[10]  Frederick R Manby,et al.  Optimization of orbital-specific virtuals in local Møller-Plesset perturbation theory. , 2012, The Journal of chemical physics.

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

[12]  Gunnar Schmitz,et al.  A scaling PNO–MP2 method using a hybrid OSV–PNO approach with an iterative direct generation of OSVs† , 2013 .

[13]  Dimitrios G Liakos,et al.  Weak Molecular Interactions Studied with Parallel Implementations of the Local Pair Natural Orbital Coupled Pair and Coupled Cluster Methods. , 2011, Journal of chemical theory and computation.

[14]  Hans-Joachim Werner,et al.  Comparison of explicitly correlated local coupled-cluster methods with various choices of virtual orbitals. , 2012, Physical chemistry chemical physics : PCCP.

[15]  Peter Pulay,et al.  Efficient calculation of canonical MP2 energies , 2001 .

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

[17]  Peter Pulay,et al.  Comparison of the boys and Pipek–Mezey localizations in the local correlation approach and automatic virtual basis selection , 1993, J. Comput. Chem..

[18]  Benny G. Johnson,et al.  A simple yet powerful upper bound for Coulomb integrals , 1994 .

[19]  Hans-Joachim Werner,et al.  Scalable electron correlation methods I.: PNO-LMP2 with linear scaling in the molecular size and near-inverse-linear scaling in the number of processors. , 2015, Journal of chemical theory and computation.

[20]  Martin Head-Gordon,et al.  A Resolution-Of-The-Identity Implementation of the Local Triatomics-In-Molecules Model for Second-Order Møller-Plesset Perturbation Theory with Application to Alanine Tetrapeptide Conformational Energies. , 2005, Journal of chemical theory and computation.

[21]  J. D. Cloizeaux,et al.  Energy Bands and Projection Operators in a Crystal: Analytic and Asymptotic Properties , 1964 .

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

[23]  Frank Neese,et al.  Natural triple excitations in local coupled cluster calculations with pair natural orbitals. , 2013, The Journal of chemical physics.

[24]  Georg Hetzer,et al.  Low-order scaling local electron correlation methods. I. Linear scaling local MP2 , 1999 .

[25]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[26]  Jan M. L. Martin,et al.  Conformational equilibria in butane-1,4-diol: a benchmark of a prototypical system with strong intramolecular H-bonds. , 2014, The journal of physical chemistry. A.

[27]  Wilfried Meyer,et al.  Ionization energies of water from PNO‐CI calculations , 2009 .

[28]  C. Hättig,et al.  A pair natural orbital based implementation of ADC(2)-x: Perspectives and challenges for response methods for singly and doubly excited states in large molecules , 2014 .

[29]  D. Tew,et al.  Explicitly correlated PNO-MP2 and PNO-CCSD and their application to the S66 set and large molecular systems. , 2014, Physical chemistry chemical physics : PCCP.

[30]  Frank Neese,et al.  Chemical applications carried out by local pair natural orbital based coupled-cluster methods. , 2014, Chemical Society reviews.

[31]  Philippe Y. Ayala,et al.  Linear scaling coupled cluster and perturbation theories in the atomic orbital basis , 1999 .

[32]  B. Paulus The method of increments—a wavefunction-based ab initio correlation method for solids , 2006 .

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

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

[35]  Matt Challacombe,et al.  An N-Body Solution to the Problem of Fock Exchange , 2014, ArXiv.

[36]  H. Lischka,et al.  PNO–CI (pair natural orbital configuration interaction) and CEPA–PNO (coupled electron pair approximation with pair natural orbitals) calculations of molecular systems. II. The molecules BeH2, BH, BH3, CH4, CH−3, NH3 (planar and pyramidal), H2O, OH+3, HF and the Ne atom , 1975 .

[37]  Frederick R Manby,et al.  Explicitly correlated local second-order perturbation theory with a frozen geminal correlation factor. , 2006, The Journal of chemical physics.

[38]  Marco Häser,et al.  Improvements on the direct SCF method , 1989 .

[39]  F. Neese,et al.  A comparative study of single reference correlation methods of the coupled-pair type , 2008 .

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

[41]  Christian Ochsenfeld,et al.  Rigorous integral screening for electron correlation methods. , 2005, The Journal of chemical physics.

[42]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

[43]  Branislav Jansík,et al.  Linear scaling coupled cluster method with correlation energy based error control. , 2010, The Journal of chemical physics.

[44]  Christian Ochsenfeld,et al.  A Linear-Scaling MP2 Method for Large Molecules by Rigorous Integral-Screening Criteria , 2010 .

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

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

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

[48]  Joachim Friedrich,et al.  Incremental CCSD(T)(F12*)|MP2: A Black Box Method To Obtain Highly Accurate Reaction Energies. , 2013, Journal of chemical theory and computation.

[49]  P. Jørgensen,et al.  A perspective on the localizability of Hartree–Fock orbitals , 2013, Theoretical Chemistry Accounts.

[50]  Dimitrios G Liakos,et al.  Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis. , 2009, The Journal of chemical physics.

[51]  Michele Benzi,et al.  Decay Properties of Spectral Projectors with Applications to Electronic Structure , 2012, SIAM Rev..

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

[53]  Jörg Kussmann,et al.  Linear-scaling atomic orbital-based second-order Møller-Plesset perturbation theory by rigorous integral screening criteria. , 2009, The Journal of chemical physics.

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

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

[56]  Edoardo Aprà,et al.  Parallel Implementation of Multireference Coupled-Cluster Theories Based on the Reference-Level Parallelism. , 2012, Journal of chemical theory and computation.

[57]  Manoj K. Kesharwani,et al.  Exploring the Accuracy Limits of Local Pair Natural Orbital Coupled-Cluster Theory. , 2015, Journal of chemical theory and computation.

[58]  John F. Stanton,et al.  A massively parallel tensor contraction framework for coupled-cluster computations , 2014, J. Parallel Distributed Comput..

[59]  R. Ahlrichs,et al.  Direct determination of pair natural orbitals , 1975 .

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

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

[62]  Martin Head-Gordon,et al.  Non-iterative local second order Møller–Plesset theory , 1998 .

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

[64]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[65]  M. Challacombe A general parallel sparse-blocked matrix multiply for linear scaling SCF theory , 2000 .

[66]  Frank Neese,et al.  An efficient and near linear scaling pair natural orbital based local coupled cluster method. , 2013, The Journal of chemical physics.

[67]  Christian Ochsenfeld,et al.  Distance-dependent Schwarz-based integral estimates for two-electron integrals: reliable tightness vs. rigorous upper bounds. , 2012, The Journal of chemical physics.

[68]  Kasper Kristensen,et al.  The divide–expand–consolidate MP2 scheme goes massively parallel , 2013 .

[69]  Jan M. L. Martin,et al.  The melatonin conformer space: benchmark and assessment of wave function and DFT methods for a paradigmatic biological and pharmacological molecule. , 2013, The journal of physical chemistry. A.

[70]  Thomas Kjærgaard,et al.  MP2 energy and density for large molecular systems with internal error control using the Divide-Expand-Consolidate scheme. , 2012, Physical chemistry chemical physics : PCCP.

[71]  Frank Neese,et al.  Geminal-spanning orbitals make explicitly correlated reduced-scaling coupled-cluster methods robust, yet simple. , 2014, The Journal of chemical physics.

[72]  F. Neese,et al.  Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. , 2009, The Journal of chemical physics.

[73]  Ajith Perera,et al.  An ab initio study of the (H2O)20H+ and (H2O)21H+ water clusters , 2009 .

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

[75]  S. F. Boys,et al.  Canonical Configurational Interaction Procedure , 1960 .

[76]  Michael Klemm,et al.  Efficient Implementation of Many-Body Quantum Chemical Methods on the Intel® Xeon Phi Coprocessor , 2014, SC14: International Conference for High Performance Computing, Networking, Storage and Analysis.

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

[78]  N. Hush,et al.  The coupled-pair approximation in a basis of independent-pair natural orbitals , 1976 .

[79]  Sohrab Ismail-Beigi,et al.  LOCALITY OF THE DENSITY MATRIX IN METALS, SEMICONDUCTORS, AND INSULATORS , 1999 .

[80]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .

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

[82]  Peter Pulay,et al.  Orbital-invariant formulation and second-order gradient evaluation in Møller-Plesset perturbation theory , 1986 .

[83]  Shuhua Li,et al.  An efficient linear scaling procedure for constructing localized orbitals of large molecules based on the one-particle density matrix. , 2011, The Journal of chemical physics.

[84]  Stinne Høst,et al.  Local orbitals by minimizing powers of the orbital variance. , 2011, The Journal of chemical physics.