A local framework for calculating coupled cluster singles and doubles excitation energies (LoFEx-CCSD)

ABSTRACT The recently developed Local Framework for calculating Excitation energies (LoFEx) is extended to the coupled cluster singles and doubles (CCSD) model. In the new scheme, a standard CCSD excitation energy calculation is carried out within a reduced excitation orbital space (XOS), which is composed of localised molecular orbitals and natural transition orbitals determined from time-dependent Hartree–Fock theory. The presented algorithm uses a series of reduced second-order approximate coupled cluster singles and doubles (CC2) calculations to optimise the XOS in a black-box manner. This ensures that the requested CCSD excitation energies have been determined to a predefined accuracy compared to a conventional CCSD calculation. We present numerical LoFEx-CCSD results for a set of medium-sized organic molecules, which illustrate the black-box nature of the approach and the computational savings obtained for transitions that are local compared to the size of the molecule. In fact, for such local transitions, the LoFEx-CCSD scheme can be applied to molecular systems where a conventional CCSD implementation is intractable.

[1]  J. Olsen,et al.  An efficient algorithm for solving nonlinear equations with a minimal number of trial vectors: applications to atomic-orbital based coupled-cluster theory. , 2008, The Journal of chemical physics.

[2]  M. Ratner.  Molecular electronic-structure theory , 2000 .

[3]  T. Daniel Crawford,et al.  Locally correlated equation-of-motion coupled cluster theory for the excited states of large molecules , 2002 .

[4]  P. Jørgensen,et al.  Large-scale calculations of excitation energies in coupled cluster theory: The singlet excited states of benzene , 1996 .

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

[6]  H. Koch,et al.  The extended CC2 model ECC2 , 2013 .

[7]  Hideo Sekino,et al.  A linear response, coupled‐cluster theory for excitation energy , 1984 .

[8]  Anna I Krylov,et al.  Equation-of-motion coupled-cluster methods for open-shell and electronically excited species: the Hitchhiker's guide to Fock space. , 2008, Annual review of physical chemistry.

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

[10]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[11]  John D. Watts,et al.  Economical triple excitation equation-of-motion coupled-cluster methods for excitation energies , 1995 .

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

[13]  Henrik Koch,et al.  The multilevel CC3 coupled cluster model. , 2016, The Journal of chemical physics.

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

[15]  Weitao Yang,et al.  Insights into Current Limitations of Density Functional Theory , 2008, Science.

[16]  P. Jørgensen,et al.  Characterization and Generation of Local Occupied and Virtual Hartree-Fock Orbitals. , 2016, Chemical reviews.

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

[18]  Jochen Schirmer,et al.  Beyond the random-phase approximation: A new approximation scheme for the polarization propagator , 1982 .

[19]  Luca Frediani,et al.  The Dalton quantum chemistry program system , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[20]  Henrik Koch,et al.  Coupled cluster response functions , 1990 .

[21]  R. Mata,et al.  An incremental correlation approach to excited state energies based on natural transition/localized orbitals. , 2011, The Journal of chemical physics.

[22]  T. Etienne.  Transition matrices and orbitals from reduced density matrix theory. , 2015, The Journal of chemical physics.

[23]  M. Head‐Gordon,et al.  Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes. , 2004, Journal of the American Chemical Society.

[24]  Trygve Helgaker,et al.  Recent advances in wave function-based methods of molecular-property calculations. , 2012, Chemical reviews.

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

[26]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[27]  John D. Watts,et al.  Iterative and non-iterative triple excitation corrections in coupled-cluster methods for excited electronic states: the EOM-CCSDT-3 and EOM-CCSD(T̃) methods , 1996 .

[28]  Kasper Kristensen,et al.  LoFEx - A local framework for calculating excitation energies: Illustrations using RI-CC2 linear response theory. , 2016, The Journal of chemical physics.

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

[30]  K. Emrich,et al.  An extension of the coupled cluster formalism to excited states: (II). Approximations and tests , 1981 .

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

[32]  John F. Stanton,et al.  Perturbative treatment of the similarity transformed Hamiltonian in equation‐of‐motion coupled‐cluster approximations , 1995 .

[33]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

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

[35]  Karol Kowalski,et al.  New coupled-cluster methods with singles, doubles, and noniterative triples for high accuracy calculations of excited electronic states. , 2004, The Journal of chemical physics.

[36]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

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

[38]  H. Monkhorst,et al.  Some aspects of the time-dependent coupled-cluster approach to dynamic response functions , 1983 .

[39]  Richard L. Martin.  NATURAL TRANSITION ORBITALS , 2003 .

[40]  K. Hirao,et al.  A generalization of the Davidson's method to large nonsymmetric eigenvalue problems , 1982 .

[41]  Trygve Helgaker,et al.  Excitation energies from the coupled cluster singles and doubles linear response function (CCSDLR). Applications to Be, CH+, CO, and H2O , 1990 .

[42]  Christof Hättig,et al.  Local pair natural orbitals for excited states. , 2011, The Journal of chemical physics.

[43]  Poul Jørgensen,et al.  Discarding Information from Previous Iterations in an Optimal Way To Solve the Coupled Cluster Amplitude Equations. , 2015, Journal of chemical theory and computation.

[44]  Christof Hättig,et al.  CC2 excitation energy calculations on large molecules using the resolution of the identity approximation , 2000 .

[45]  Martin Schütz,et al.  A multistate local coupled cluster CC2 response method based on the Laplace transform. , 2009, The Journal of chemical physics.

[46]  J. Cizek.  On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in Ursell-Type Expansion Using Quantum-Field Theoretical Methods , 1966 .

[47]  John F. Stanton,et al.  The equation of motion coupled‐cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties , 1993 .

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

[49]  Ove Christiansen,et al.  Response functions in the CC3 iterative triple excitation model , 1995 .

[50]  H. Monkhorst,et al.  Calculation of properties with the coupled-cluster method , 2009 .

[51]  Christof Hättig,et al.  Transition moments and excited-state first-order properties in the coupled-cluster model CC2 using the resolution-of-the-identity approximation , 2002 .

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

[53]  C. Hättig,et al.  A pair natural orbital implementation of the coupled cluster model CC2 for excitation energies. , 2013, The Journal of chemical physics.

[54]  H. Koch,et al.  Multi-level coupled cluster theory. , 2014, The Journal of chemical physics.

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

[56]  Frank Neese,et al.  Towards a pair natural orbital coupled cluster method for excited states. , 2016, The Journal of chemical physics.

[57]  R. Bartlett,et al.  Simplified methods for equation-of-motion coupled-cluster excited state calculations , 1996 .

[58]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[59]  Poul Jørgensen,et al.  Perturbative triple excitation corrections to coupled cluster singles and doubles excitation energies , 1996 .