Orbital-optimized linearized coupled-cluster doubles with density-fitting and Cholesky decomposition approximations: an efficient implementation.

An efficient implementation of the orbital-optimized linearized coupled-cluster double method with the density-fitting (DF-OLCCD) and Cholesky decomposition (CD-OLCCD) approximations is presented. The DF-OLCCD and CD-OLCCD methods are applied to a set of alkanes to compare the computational cost with the conventional orbital-optimized linearized coupled-cluster doubles (OLCCD) [U. Bozkaya and C. D. Sherrill, J. Chem. Phys., 2013, 139, 054104]. Our results demonstrate that the DF-OLCCD method provides substantially lower computational costs than OLCCD, and there are more than 9-fold reductions in the computational time for the largest member of the alkane set (C8H18). For barrier heights of hydrogen transfer reaction energies, the DF-OLCCD method again exhibits a substantially better performance than DF-LCCD, providing a mean absolute error of 0.9 kcal mol(-1), which is 7 times lower than that of DF-LCCD (6.2 kcal mol(-1)), and compared to MP2 (9.6 kcal mol(-1)) there is a more than 10-fold reduction in errors. Furthermore, the MAE value of DF-OLCCD is also lower than that of CCSD (1.2 kcal mol(-1)). For open-shell noncovalent interactions, the performance of DF-OLCCD is significantly better than that of MP2, DF-LCCD, and CCSD. Overall, the present application results indicate that the DF-OLCCD and CD-OLCCD methods are very promising for challenging open-shell systems as well as closed-shell molecular systems.

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

[2]  Uğur Bozkaya,et al.  Orbital-optimized coupled-electron pair theory and its analytic gradients: accurate equilibrium geometries, harmonic vibrational frequencies, and hydrogen transfer reactions. , 2013, The Journal of chemical physics.

[3]  M. Duncan,et al.  IR photodissociation spectroscopy of O4+, O6+ and O8+ cluster ions , 2009 .

[4]  U. Bozkaya,et al.  Assessment of Orbital-Optimized MP2.5 for Thermochemistry and Kinetics: Dramatic Failures of Standard Perturbation Theory Approaches for Aromatic Bond Dissociation Energies and Barrier Heights of Radical Reactions. , 2015, Journal of chemical theory and computation.

[5]  H. Schaefer,et al.  Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method: improving upon CCSD(T) and CCSD(T)(Λ): preliminary application. , 2012, The Journal of chemical physics.

[6]  Thom H. Dunning,et al.  Gaussian basis sets for use in correlated molecular calculations. V. Core-valence basis sets for boron through neon , 1995 .

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

[8]  M. Jacox,et al.  Evidence for the stabilization of rectangular O+4 in solid neon , 1994 .

[9]  John D. Watts,et al.  Non-iterative fifth-order triple and quadruple excitation energy corrections in correlated methods , 1990 .

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

[11]  Uğur Bozkaya,et al.  Assessment of the extended Koopmans' theorem for the chemical reactivity: Accurate computations of chemical potentials, chemical hardnesses, and electrophilicity indices , 2016, J. Comput. Chem..

[12]  Uğur Bozkaya,et al.  Quadratically convergent algorithm for orbital optimization in the orbital-optimized coupled-cluster doubles method and in orbital-optimized second-order Møller-Plesset perturbation theory. , 2011, The Journal of chemical physics.

[13]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[14]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[15]  R. Lindh,et al.  Symmetry breaking in O+4: an application of the Brueckner coupled-cluster method , 1994 .

[16]  Rohini C. Lochan,et al.  Exploring the competition between localization and delocalization of the neutral soliton defect in polyenyl chains with the orbital optimized second order opposite spin method. , 2012, The Journal of chemical physics.

[17]  M. Jacox,et al.  The vibrational spectra of molecular ions isolated in solid neon. II. O+4 and O−4 , 1989 .

[18]  T. Daniel Crawford,et al.  Spin-restricted Brueckner orbitals for coupled-cluster wavefunctions , 1997 .

[19]  Thomas Bondo Pedersen,et al.  Reduced scaling in electronic structure calculations using Cholesky decompositions , 2003 .

[20]  Ugur Bozkaya,et al.  Assessment of the orbital‐optimized coupled‐electron pair theory for thermochemistry and kinetics: Improving on CCSD and CEPA(1) , 2014, J. Comput. Chem..

[21]  H. Koch,et al.  Gauge invariant coupled cluster response theory , 1999 .

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

[23]  Thomas Bondo Pedersen,et al.  Gauge invariant coupled cluster response theory using optimized nonorthogonal orbitals , 2001 .

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

[25]  Alexander Yu Sokolov,et al.  Orbital-optimized density cumulant functional theory. , 2013, The Journal of chemical physics.

[26]  Uğur Bozkaya,et al.  Orbital-Optimized Second-Order Perturbation Theory with Density-Fitting and Cholesky Decomposition Approximations: An Efficient Implementation. , 2014, Journal of chemical theory and computation.

[27]  Alistair P. Rendell,et al.  COUPLED-CLUSTER THEORY EMPLOYING APPROXIMATE INTEGRALS : AN APPROACH TO AVOID THE INPUT/OUTPUT AND STORAGE BOTTLENECKS , 1994 .

[28]  R. Bartlett,et al.  On the choice of orbitals for symmetry breaking problems with application to NO3 , 1992 .

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

[30]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

[31]  R. Bartlett,et al.  Property evaluation and orbital relaxation in coupled cluster methods , 1987 .

[32]  Uğur Bozkaya,et al.  Accurate Electron Affinities from the Extended Koopmans' Theorem Based on Orbital-Optimized Methods. , 2014, Journal of chemical theory and computation.

[33]  U. Bozkaya,et al.  Assessment of Orbital-Optimized Third-Order Møller-Plesset Perturbation Theory and Its Spin-Component and Spin-Opposite Scaled Variants for Thermochemistry and Kinetics. , 2013, Journal of chemical theory and computation.

[34]  R. Lindh,et al.  The fraternal twins of quartet O+4 , 1994 .

[35]  D R Yarkony,et al.  Modern electronic structure theory , 1995 .

[36]  Frank Neese,et al.  Correlated ab initio spin densities for larger molecules: orbital-optimized spin-component-scaled MP2 method. , 2010, The journal of physical chemistry. A.

[37]  Matthew L. Leininger,et al.  Psi4: an open‐source ab initio electronic structure program , 2012 .

[38]  R. Shepard.  The Multiconfiguration Self‐Consistent Field Method , 2007 .

[39]  Trygve Helgaker,et al.  Mo/ller–Plesset energy derivatives , 1988 .

[40]  Poul Jo,et al.  Optimization of orbitals for multiconfigurational reference states , 1978 .

[41]  Anna I. Krylov,et al.  Energies and analytic gradients for a coupled-cluster doubles model using variational Brueckner orbitals: Application to symmetry breaking in O4+ , 1998 .

[42]  Uğur Bozkaya,et al.  Derivation of general analytic gradient expressions for density-fitted post-Hartree-Fock methods: an efficient implementation for the density-fitted second-order Møller-Plesset perturbation theory. , 2014, The Journal of chemical physics.

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

[44]  Anna I. Krylov,et al.  Second order perturbation corrections to singles and doubles coupled-cluster methods: General theory and application to the valence optimized doubles model , 2000 .

[45]  Anna I. Krylov,et al.  Size-consistent wave functions for nondynamical correlation energy: The valence active space optimized orbital coupled-cluster doubles model , 1998 .

[46]  Peter J Knowles,et al.  Rigorously extensive orbital-invariant renormalized perturbative triples corrections from quasi-variational coupled cluster theory. , 2013, The Journal of chemical physics.

[47]  P. Knowles,et al.  Benchmark Quasi-Variational Coupled Cluster Calculations of Multiple Bond Breaking. , 2012, Journal of chemical theory and computation.

[48]  A Eugene DePrince,et al.  Accuracy and Efficiency of Coupled-Cluster Theory Using Density Fitting/Cholesky Decomposition, Frozen Natural Orbitals, and a t1-Transformed Hamiltonian. , 2013, Journal of chemical theory and computation.

[49]  Rohini C. Lochan,et al.  Orbital-optimized opposite-spin scaled second-order correlation: an economical method to improve the description of open-shell molecules. , 2007, The Journal of chemical physics.

[50]  Uğur Bozkaya,et al.  Orbital-Optimized MP3 and MP2.5 with Density-Fitting and Cholesky Decomposition Approximations. , 2016, Journal of chemical theory and computation.

[51]  Uğur Bozkaya,et al.  The extended Koopmans' theorem for orbital-optimized methods: accurate computation of ionization potentials. , 2013, The Journal of chemical physics.

[52]  F. Neese,et al.  An orbital-invariant and strictly size extensive post-Hartree-Fock correlation functional. , 2011, The Journal of chemical physics.

[53]  Hans-Joachim Werner,et al.  An efficient local coupled cluster method for accurate thermochemistry of large systems. , 2011, The Journal of chemical physics.

[54]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[55]  Uğur Bozkaya,et al.  Orbital-optimized third-order Møller-Plesset perturbation theory and its spin-component and spin-opposite scaled variants: application to symmetry breaking problems. , 2011, The Journal of chemical physics.

[56]  Rodney J. Bartlett,et al.  The expectation value coupled-cluster method and analytical energy derivatives☆ , 1988 .

[57]  M. Jacox.  Vibrational and Electronic Energy Levels of Polyatomic Transient Molecules. Supplement B , 1990 .

[58]  M. Head‐Gordon,et al.  Violations of N-representability from spin-unrestricted orbitals in Møller–Plesset perturbation theory and related double-hybrid density functional theory , 2009 .

[59]  Uğur Bozkaya,et al.  Accurate Open-Shell Noncovalent Interaction Energies from the Orbital-Optimized Møller-Plesset Perturbation Theory: Achieving CCSD Quality at the MP2 Level by Orbital Optimization. , 2013, Journal of chemical theory and computation.

[60]  R. J. Boyd,et al.  A comparative study of the hyperfine structures of neutral nitrogen oxides: DFT vs CISD results , 1994 .

[61]  Trygve Helgaker,et al.  Analytical Calculation of Geometrical Derivatives in Molecular Electronic Structure Theory , 1988 .

[62]  Anna I. Krylov,et al.  Excited states theory for optimized orbitals and valence optimized orbitals coupled-cluster doubles models , 2000 .

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

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

[65]  K. Kawaguchi,et al.  Infrared spectroscopy of the NO3 radical , 1998 .

[66]  M. Head‐Gordon,et al.  Approaching closed-shell accuracy for radicals using coupled cluster theory with perturbative triple substitutions , 2003 .

[67]  Frank Neese,et al.  Assessment of Orbital-Optimized, Spin-Component Scaled Second-Order Many-Body Perturbation Theory for Thermochemistry and Kinetics. , 2009, Journal of chemical theory and computation.

[68]  Uğur Bozkaya,et al.  Analytic Energy Gradients and Spin Multiplicities for Orbital-Optimized Second-Order Perturbation Theory with Density-Fitting Approximation: An Efficient Implementation. , 2014, Journal of chemical theory and computation.

[69]  Uğur Bozkaya,et al.  Orbital-optimized MP2.5 and its analytic gradients: approaching CCSD(T) quality for noncovalent interactions. , 2014, The Journal of chemical physics.

[70]  U. Bozkaya.  Analytic energy gradients for the orbital-optimized third-order Møller-Plesset perturbation theory. , 2013, The Journal of chemical physics.

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

[72]  Martin Head-Gordon,et al.  The Theoretical Prediction of Molecular Radical Species: a Systematic Study of Equilibrium Geometries and Harmonic Vibrational Frequencies , 2001 .

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

[74]  C David Sherrill,et al.  High-level ab initio studies of hydrogen abstraction from prototype hydrocarbon systems. , 2006, The journal of physical chemistry. A.

[75]  P. Knowles,et al.  Approximate variational coupled cluster theory. , 2011, The Journal of chemical physics.

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

[77]  A. Hesselmann,et al.  The role of orbital transformations in coupled-pair functionals , 2010 .

[78]  J. Olsen,et al.  Orbital-optimized coupled-cluster theory does not reproduce the full configuration-interaction limit. , 2005, The Journal of chemical physics.

[79]  G. Scuseria,et al.  The optimization of molecular orbitals for coupled cluster wavefunctions , 1987 .

[80]  Peter J Knowles,et al.  Quasi-variational coupled cluster theory. , 2012, The Journal of chemical physics.