Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis.

A production level implementation of the closed-shell local quadratic configuration interaction and coupled cluster methods with single and double excitations (QCISD and CCSD) based on the concept of pair natural orbitals [local pair natural orbital LPNO-QCISD and LPNO-CCSD) is reported, evaluated, and discussed. This work is an extension of the earlier developed LPNO coupled-electron pair approximation (LNPO-CEPA) method [F. Neese et al., Chem. Phys. 130, 114108 (2009)] and makes extended use of the resolution of the identity (RI) or density fitting (DF) approximation. Two variants of each method are compared. The less accurate approximations (LPNO2-QCISD/LPNO2-CCSD) still recover 98.7%-99.3% of the correlation energy in the given basis and have modest disk space requirements. The more accurate variants (LPNO1-QCISD/LPNO1-CCSD) typically recover 99.75%-99.95% of the correlation energy in the given basis but require the Coulomb and exchange operators with up to two-external indices to be stored on disk. Both variants have comparable computational efficiency. The convergence of the results with respect to the natural orbital truncation parameter (T(CutPNO)) has been studied. Extended numerical tests have been performed on absolute and relative correlation energies as function of basis set size and T(CutPNO) as well as on reaction energies, isomerization energies, and weak intermolecular interactions. The results indicate that the errors of the LPNO methods compared to the canonical QCISD and CCSD methods are below 1 kcal/mol with our default thresholds. Finally, some calculations on larger molecules are reported (ranging from 40-86 atoms) and it is shown that for medium sized molecules the total wall clock time required to complete the LPNO-CCSD calculations is only two to four times that of the preceding self-consistent field (SCF). Thus these methods are highly suitable for large-scale computational chemistry applications. Since there are only three thresholds involved that have been given conservative default values, the methods can be confidentially used in a "black-box" fashion in the same way as their canonical counterparts.

[1]  Guntram Rauhut,et al.  Integral transformation with low‐order scaling for large local second‐order Møller–Plesset calculations , 1998 .

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

[3]  Edward F. Valeev,et al.  Simple coupled-cluster singles and doubles method with perturbative inclusion of triples and explicitly correlated geminals: The CCSD(T)R12 model. , 2008, The Journal of chemical physics.

[4]  Wilfried Meyer,et al.  Configuration Expansion by Means of Pseudonatural Orbitals , 1977 .

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

[6]  T Daniel Crawford,et al.  Potential energy surface discontinuities in local correlation methods. , 2004, The Journal of chemical physics.

[7]  Guntram Rauhut,et al.  Analytical energy gradients for local second-order Mo/ller–Plesset perturbation theory , 1998 .

[8]  Tatiana Korona,et al.  Local CC2 electronic excitation energies for large molecules with density fitting. , 2006, The Journal of chemical physics.

[9]  Peter Pulay,et al.  Parallel Calculation of Coupled Cluster Singles and Doubles Wave Functions Using Array Files. , 2007, Journal of chemical theory and computation.

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

[11]  A. Savin,et al.  Coupled-cluster calculations using local potentials , 2002 .

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

[13]  M. Frisch,et al.  The Theoretical Calculation of Vibrational Circular Dichroism Spectra. , 1995 .

[14]  R. Ahlrichs,et al.  The Method of Self Consistent Electron Pairs. A Matrix Oriented Direct CI , 1984 .

[15]  Curtis L. Janssen,et al.  An efficient reformulation of the closed‐shell coupled cluster single and double excitation (CCSD) equations , 1988 .

[16]  Juana Vázquez,et al.  HEAT: High accuracy extrapolated ab initio thermochemistry. , 2004, The Journal of chemical physics.

[17]  Guntram Rauhut,et al.  Accurate calculation of anharmonic vibrational frequencies of medium sized molecules using local coupled cluster methods. , 2007, The Journal of chemical physics.

[18]  P. Taylor.  A rapidly convergent CI expansion based on several reference configurations, using optimized correlating orbitals , 1981 .

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

[20]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

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

[22]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[23]  K. P. Lawley,et al.  Ab initio methods in quantum chemistry , 1987 .

[24]  Hans-Joachim Werner,et al.  Calculation of smooth potential energy surfaces using local electron correlation methods. , 2006, The Journal of chemical physics.

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

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

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

[28]  Frederick R Manby,et al.  Analytical energy gradients for local second-order Møller-Plesset perturbation theory using density fitting approximations. , 2004, The Journal of chemical physics.

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

[30]  Georg Hetzer,et al.  Low-order scaling local correlation methods II: Splitting the Coulomb operator in linear scaling local second-order Møller–Plesset perturbation theory , 2000 .

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

[32]  Donald G. Truhlar,et al.  Basis-set extrapolation , 1998 .

[33]  M. Krauss,et al.  Pseudonatural Orbitals as a Basis for the Superposition of Configurations. I. He2 , 1966 .

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

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

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

[37]  H. Werner,et al.  The effect of local approximations in coupled-cluster wave functions on dipole moments and static dipole polarisabilitiesDedicated to Prof. W. Meyer on the occasion of his 65th birthday. , 2004 .

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

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

[40]  S. F. Boys.  Construction of Some Molecular Orbitals to Be Approximately Invariant for Changes from One Molecule to Another , 1960 .

[41]  Peter Pulay,et al.  Efficient elimination of basis set superposition errors by the local correlation method: Accurate ab initio studies of the water dimer , 1993 .

[42]  M. Nooijen,et al.  Dynamically screened local correlation method using enveloping localized orbitals. , 2006, The Journal of chemical physics.

[43]  Rodney J. Bartlett,et al.  Correlation energy estimates in periodic extended systems using the localized natural bond orbital coupled cluster approach , 2003 .

[44]  Peter Pulay,et al.  An efficient reformulation of the closed‐shell self‐consistent electron pair theory , 1984 .

[45]  Frederick R Manby,et al.  Local explicitly correlated second-order perturbation theory for the accurate treatment of large molecules. , 2009, The Journal of chemical physics.

[46]  Martin Head-Gordon,et al.  Quadratic configuration interaction. A general technique for determining electron correlation energies , 1987 .

[47]  T. Crawford,et al.  An Introduction to Coupled Cluster Theory for Computational Chemists , 2007 .

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

[49]  R. Bartlett,et al.  Localized correlation treatment using natural bond orbitals , 2003 .

[50]  R. Mcweeny,et al.  Methods Of Molecular Quantum Mechanics , 1969 .

[51]  J. Malrieu,et al.  Eigenvalue problem formulation of coupled-cluster expansions through intermediate Hamiltonians , 1995 .

[52]  An Overview of Molecular Quantum Mechanics , 1992 .

[53]  J. Gauss,et al.  Analytical Energy Gradients in Møller-Plesset Perturbation and Quadratic Configuration Interaction Methods: Theory and Application , 1992 .

[54]  Martin Schütz,et al.  A new, fast, semi-direct implementation of linear scaling local coupled cluster theory , 2002 .

[55]  R. Bartlett,et al.  Localized orbitals in the coupled cluster singles and doubles model , 2009 .

[56]  H. Stoll.  On the correlation energy of graphite , 1992 .

[57]  Michael Dolg,et al.  Fully automated implementation of the incremental scheme: application to CCSD energies for hydrocarbons and transition metal compounds. , 2007, The Journal of chemical physics.

[58]  Guntram Rauhut,et al.  Analytical energy gradients for local coupled-cluster methods , 2001 .

[59]  P. Kollman,et al.  Encyclopedia of computational chemistry , 1998 .

[60]  Trygve Helgaker,et al.  Quantitative quantum chemistry , 2008 .

[61]  Paul G. Mezey,et al.  A fast intrinsic localization procedure applicable for ab initio and semiempirical linear combination of atomic orbital wave functions , 1989 .

[62]  R. Ahlrichs,et al.  Many body perturbation calculations and coupled electron pair models , 1979 .

[63]  R. Bartlett,et al.  Can simple localized bond orbitals and coupled cluster methods predict reliable molecular energies , 1985 .

[64]  R. Bartlett,et al.  Coupled-cluster theory in quantum chemistry , 2007 .

[65]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[66]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[67]  W. Kutzelnigg,et al.  Ab initio calculations on small hydrides including electron correlation , 1968 .

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

[69]  T. Daniel Crawford,et al.  Local correlation in coupled cluster calculations of molecular response properties , 2004 .

[70]  Martin Head-Gordon,et al.  Fast localized orthonormal virtual orbitals which depend smoothly on nuclear coordinates. , 2005, The Journal of chemical physics.

[71]  Clifford E. Dykstra,et al.  Advanced theories and computational approaches to the electronic structure of molecules , 1984 .

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

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

[74]  Kazuo Kitaura,et al.  Coupled-cluster theory based upon the fragment molecular-orbital method. , 2005, The Journal of chemical physics.

[75]  Electron correlation in large molecular systems using the atomic orbital formalism. The case of intermolecular interactions in crystalline urea as an example , 2000 .

[76]  Joseph G. Hoffman,et al.  Quantum Theory of Atoms, Molecules and the Solid State: A Tribute to John C. Slater , 1967 .

[77]  Werner Kutzelnigg,et al.  Direct Calculation of Approximate Natural Orbitals and Natural Expansion Coefficients of Atomic and Molecular Electronic Wavefunctions. II. Decoupling of the Pair Equations and Calculation of the Pair Correlation Energies for the Be and LiH Ground States , 1968 .

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

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

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

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

[82]  Hans-Joachim Werner,et al.  Correlation regions within a localized molecular orbital approach. , 2008, The Journal of chemical physics.

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

[84]  Rodney J Bartlett,et al.  A natural linear scaling coupled-cluster method. , 2004, The Journal of chemical physics.

[85]  M. Krauss,et al.  Configuration‐Interaction Calculation of H3 and H2 , 1965 .

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

[87]  P. Taylor,et al.  Application of Gaussian-type geminals in local second-order Møller-Plesset perturbation theory. , 2006, The Journal of chemical physics.

[88]  Reinhold F. Fink,et al.  A multi-configuration reference CEPA method based on pair natural orbitals , 1993 .

[89]  P Pulay,et al.  Local Treatment of Electron Correlation , 1993 .

[90]  Hans-Joachim Werner,et al.  Eliminating the domain error in local explicitly correlated second-order Møller-Plesset perturbation theory. , 2008, The Journal of chemical physics.

[91]  H. Schaefer.  Methods of Electronic Structure Theory , 1977 .

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

[93]  Jing Ma,et al.  Linear scaling local correlation approach for solving the coupled cluster equations of large systems , 2002, J. Comput. Chem..

[94]  CEPA calculations on open-shell molecules. I. Outline of the method , 1981 .

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

[96]  Tatiana Korona,et al.  Transition strengths and first-order properties of excited states from local coupled cluster CC2 response theory with density fitting. , 2007, The Journal of chemical physics.

[97]  Trygve Helgaker,et al.  Basis-set convergence in correlated calculations on Ne, N2, and H2O , 1998 .

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

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

[100]  F. Manby,et al.  Explicitly correlated second-order perturbation theory using density fitting and local approximations. , 2006, The Journal of chemical physics.

[101]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[102]  Hans-Joachim Werner,et al.  Simplified CCSD(T)-F12 methods: theory and benchmarks. , 2009, The Journal of chemical physics.

[103]  H. Werner,et al.  Chapter 4 On the Selection of Domains and Orbital Pairs in Local Correlation Treatments , 2006 .

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

[105]  Frank Neese,et al.  Accurate theoretical chemistry with coupled pair models. , 2009, Accounts of chemical research.

[106]  Peter R. Taylor,et al.  Unlinked cluster effects in molecular electronic structure. I. The HCN and HNC molecules , 1978 .

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

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

[109]  Guntram Rauhut,et al.  Impact of local and density fitting approximations on harmonic vibrational frequencies. , 2006, The journal of physical chemistry. A.

[110]  A. Dutoi,et al.  Accurate local approximations to the triples correlation energy: formulation, implementation and tests of 5th-order scaling models , 2005 .

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

[112]  Joseph E. Subotnik,et al.  The limits of local correlation theory: electronic delocalization and chemically smooth potential energy surfaces. , 2008, The Journal of chemical physics.

[113]  Hans-Joachim Werner,et al.  A comparison of the efficiency and accuracy of the quadratic configuration interaction (QCISD), coupled cluster (CCSD), and Brueckner coupled cluster (BCCD) methods , 1992 .

[114]  Wilfried Meyer,et al.  Theory of self‐consistent electron pairs. An iterative method for correlated many‐electron wavefunctions , 1976 .

[115]  H. Werner,et al.  Local treatment of electron excitations in the EOM-CCSD method , 2003 .