Coupled cluster response equation

In this section we will introduce some wavefunction-based methods to calculate photoabsorption spectra. The Hartree-Fock method itself is a wavefunction-based approach to solve the static Schrodinger equation. For excited states one has to account for time-dependent phenomena as in the density-based approaches. Therefore, we will start with a short review of time-dependent Hartree-Fock. Several more advanced methods are available as well, e.g. configuration interaction (Cl), multireference configuration interaction (MRCI), multireference Moller-Plesset (MRMP), or complete active space self-consistent field (CASSCF), to name only a few. Also flavours of the coupled-cluster approach (equations-of-motion CC and linear-response CQ are used to calculate excited states. However, all these methods are applicable only to fairly small molecules due to their high computational costs. These approaches are therefore discussed only in a more phenomenological way here, and many post-Hartree-Fock methods are explicitly not included. 

In recent years, these methods have been greatly expanded and have reached a degree of reliability where they now offer some of the most accurate tools for studying excited and ionized states. In particular, the use of time-dependent variational principles have allowed the much more rigorous development of equations for energy differences and nonlinear response properties [81]. In addition, the extension of the EOM theory to include coupled-cluster reference fiuictioiis [ ] now allows one to compute excitation and ionization energies using some of the most accurate ab initio tools. 

In the next section we derive the Taylor expansion of the coupled cluster cubic response function in its frequency arguments and the equations for the required expansions of the cluster amplitude and Lagrangian multiplier responses. For the experimentally important isotropic averages 7, 7i and yx we give explicit expressions for the A and higher-order coefficients in terms of the coefficients of the Taylor series. In Sec. 4 we present an application of the developed approach to the second hyperpolarizability of the methane molecule. We test the convergence of the hyperpolarizabilities with respect to the order of the expansion and investigate the sensitivity of the coefficients to basis sets and correlation treatment. The results are compared with dispersion coefficients derived by least square fits to experimental hyperpolarizability data or to pointwise calculated hyperpolarizabilities of other ab inito studies. 

The field- and time-dependent cluster operator is defined as T t, ) = nd HF) is the SCF wavefunction of the unperturbed molecule. By keeping the Hartree-Fock reference fixed in the presence of the external perturbation, a two step approach, which would introduce into the coupled cluster wavefunction an artificial pole structure form the response of the Hartree Fock orbitals, is circumvented. The quasienergy W and the time-dependent coupled cluster equations are determined by projecting the time-dependent Schrodinger equation onto the Hartree-Fock reference and onto the bra states (HF f[[exp(—T) ... 

In two recent publications we have tried to characterize the excited state properties of 1 and 3 in order to facilitate their detection by LIF-spectroscopy. Our main tool in this effort has been equation of motion coupled cluster theory (EOM-CC). The EOM-CCSD method, which is equivalent to linear response CCSD, has been shown to provide an accurate description of both valence and excited states even in systems where electron correlation effects play an important role [39]. Computed transition energies for excitations that are of mainly single substitution character are generally accurate to within 0.1 eV. We have found the EOM-CCSD method to perform particularly well in combination with the doubly-augmented cc-pVDZ (d-aug-cc-pVDZ) basis set. This basis seems to provide equally balanced descriptions of ground and excited states,... 

M. Nooijen, J.G. Snijders, Int. J. Quantum Chem. 48 (1993) 15. For ionization potentials, equation-of-motion coupled-cluster, coupled-cluster linear response theory,... 

Also in response theory the summation over excited states is effectively replaced by solving a system of linear equations. Spin-orbit matrix elements are obtained from linear response functions, whereas quadratic response functions can most elegantly be utilized to compute spin-forbidden radiative transition probabilities. We refrain from going into details here, because an excellent review on this subject has been published by Agren et al.118 While these authors focus on response theory and its application in the framework of Cl and multiconfiguration self-consistent field (MCSCF) procedures, an analogous scheme using coupled-cluster electronic structure methods was presented lately by Christiansen et al.124... 

More demanding calculations conducted by Daniel et al. using equations-of-motion coupled cluster methods (EOM-CC) provided a better description of the contribution of MLCT andMC transitions to the electronic spectrum of Cr(CO)6 [13]. These calculations confirm that the lowest-energy transition, which is responsible for the shoulder on the low-energy side of the first absorption maximum, corresponds to... 

AN INTRODUCTION TO EQUATION-OF-MOTION AND LINEAR-RESPONSE COUPLED-CLUSTER METHODS FOR ELECTRONICALLY EXCITED STATES OF MOLECULES... 

Keywords Equation-of-Motion Coupled-Cluster Theory, Linear-Response Coupled-Cluster Theory,... 

The approach outlined above combines the calculation of response functions (i.e. of frequency-dependent properties) with the theory of analytic derivatives developed for static higher-order properties. In the limit of a static perturbation all equations above reduce to the usual equations for (unrelaxed) coupled cluster energy derivatives. This is an invaluable advantage for the implementation of frequency-dependent properties in quantum chemistry programs. 

---