Coupled-cluster linear response theory

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,... 

Abstract The modified equation-of-motion coupled cluster approach of Sekino and Bartlett is extended to computations of the mixed electric-dipole/magnetic-dipole polarizability tensor associated with optical rotation in chiral systems. The approach - referred to here as a linearized equation-of-motion coupled cluster (EOM-CCl) method - is a compromise between the standard EOM method and its linear response counterpart, which avoids the evaluation of computationally expensive terms that are quadratic in the field-perturbed wave functions, but still yields properties that are size-extensive/intensive. Benchmark computations on five representative chiral molecules, including (P)-hydrogen peroxide, (5)-methyloxirane, (5 )-2-chloropropioniuile, (/ )-epichlorohydrin, and (75,45)-norbornenone, demonstrate typically small deviations between the EOM-CCl results and those from coupled cluster linear response theory, and no variation in the signs of the predicted rotations. In addition, the EOM-CCl approach is found to reduce the overall computing time for multi-wavelength-specific rotation computations by up to 34%. 

The effects of including the triple excitations in coupled cluster linear response theory for evaluating the dynamic polarizabilities have been assessed for a set of closed-shell (Ne, HF, N2, CO) and open-shell (CN, CO, O2) systems, in view of exploring a new accuracy regime for molecular properties. The main conclusions include that i) for systems with little or no static correlation, CC3 is nearly identical to CCSDT, ii) CC3 and PS(T) [pole shifted technique where the CCSD-LR poles are corrected by adding a noniterative correction due to the triples] methods perform better than CCSD but their relative accuracy is not determined yet, iii) differences between CCSD and CC3 results as well as the errors with respect to CCSDT drop when the basis set is increased, and iv) ROHF-based CC-LR approaches should be favored over their UHF counterparts while the dilfer-ences between the ROHF and UHF appear as an appropriate criterion for determining whether higher-order UHF-based CC calculations can be used. 

Coriani S, Fransson T, Christiansen O, Norman P. Asymmetric-Lanczos-Chain-Driven Implementation of Electronic Resonance Convergent Coupled-Cluster Linear Response Theory. J Chem Theory Comput. 2012 8 1616-1628. 

Another category of approaches that avoids the symmetry breaking problem of the orbitals for the target state is based on using a related, well-behaved HF reference instead. Examples of such techniques include quasi-restricted Hartree-Fock coupled-cluster (QRHF CC) (11,19), symmetry adapted cluster configuration interaction (SAC-CI) (22,23), coupled-cluster linear response theory (CCLRT) (24-26) or equivalently equation-of-motion coupled-cluster (EOM-CC) (27-32), Fock space multi-reference coupled-cluster (FSMRCC) (33-37), and similarity transformed equation-of-motion coupled-cluster (STEOM-CC) (38-40). In the case of NO3 and N03, the restricted Hartree-Fock (RHF) orbitals of the closed-shell N03 anion system can be used as the reference. The anion orbitals are stable with respect to symmetry perturbations, and the system does not suffer from the artifactual symmetry breaking of the neutral and cation. 

A wide range of first-order properties of the electronic excited-states of molecules in solution can be computed in terms of the gradients of the excitation energies coK with respect to external or internal perturbations. In this section, we discuss the analytical theory for the gradient of the excitation energies (4.11) computed from the coupled-cluster linear response. 

The last method used in this study is CCSD linear response theory [37]. The frequency-dependent polarizabilities are again identified from the time evolution of the corresponding moments. However, in CCSD response theory the moments are calculated as transition expectation values between the coupled cluster state l cc(O) and a dual state... 

However, until today no systematic comparison of methods based on MpUer-Plesset perturbation (MP) and Coupled Cluster theory, the SOPPA or multiconfigurational linear response theory has been presented. The present study is a first attempt to remedy this situation. Calculations of the rotational g factor of HF, H2O, NH3 and CH4 were carried out at the level of Hartree-Fock (SCF) and multiconfigurational Hartree-Fock (MCSCF) linear response theory, the SOPPA and SOPPA(CCSD) [40], MpUer-Plesset perturbation theory to second (MP2), third (MP3) and fourth order without the triples contributions (MP4SDQ) and finally coupled cluster singles and doubles theory. The same basis sets and geometries were employed in all calculations for a given molecule. The results obtained with the different methods are therefore for the first time direct comparable and consistent conclusions about the performance of the different methods can be made. 

LINEAR RESPONSE THEORY IN CONNECTION TO DENSITY FUNCTIONAL THEORY/MOLECULAR DYNAMICS AND COUPLED CLUSTER/MOLECULAR DYNAMICS METHODS... 

Trimerized organic conductors are of special interest, because two electrons per three sites constitute the simplest situation, where both electronic transitions resulting in single- and double-site occupation take place [21]. As one considers larger n-mers, two complications arise. First, the number of equations that should be solved sharply increases. The second complication is the increase in the number of n-meric normal modes, which are coupled to an external electromagnetic field. Recently, Yartsev et al. [22] have proposed using the linear response theory for several variables to describe the optical properties of trimers with arbitrary equilibrium charge density distribution. This approach can be extended to any cluster—the size is limited only by computer facilities. 

Debashis Mukherjee is a Professor of Physical Chemistry and the Director of the Indian Association for the Cultivation of Science, Calcutta, India. He has been one of the earliest developers of a class of multi-reference coupled cluster theories and also of the coupled cluster based linear response theory. Other contributions by him are in the resolution of the size-extensivity problem for multi-reference theories using an incomplete model space and in the size-extensive intermediate Hamiltonian formalism. His research interests focus on the development and applications of non-relativistic and relativistic theories of many-body molecular electronic structure and theoretical spectroscopy, quantum many-body dynamics and statistical held theory of many-body systems. He is a member of the International Academy of the Quantum Molecular Science, a Fellow of the Third World Academy of Science, the Indian National Science Academy and the Indian Academy of Sciences. He is the recipient of the Shantiswarup Bhatnagar Prize of the Council of Scientihc and Industrial Research of the Government of India. 

Keywords Coupled cluster theory linear response theory equation-of-motion coupled cluster theory optical rotation chirality... 

Pedersen, T. B., Sanchez de Meras, A. M.and Koch, H. (2004). Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions, /. Chem. Phys. 120, pp. 8887-8897, doiilO.1063/1. 1705575. 

At the wavefunction level methods based on coupled cluster (CC) theory are among the most reliable ones. For ground-state energetics the CCSD(T) approach is the gold standard of chemistry, whereas for excited states one can use the equation-of-motion (EOM) CC (EOM-CC) method or CC linear response theory (CC-LRT) [4] approaches. Note that the CC-LRT is size-extensive for both energies and properties such as intensities, however for EOM-CC this is true only for energies (unless one uses the closely related similarity transformed (ST)EOM-CC method [18]). As the computational cost of fully iterative (e.g., CCSD, CCSDT, etc.) methods can quickly become prohibitive, perturbative methods [4]... 

In the CC based linear response theory, the excited functions are written in terms of the ground state function via the action of a suitable excitation operator. In the spirit of the single-reference coupled-cluster based Unear response theory (SR-CCLRT) [38, 44, 45] (or SAC-CI [46]) we posit on our excited state, [V fc), the foUow-ing Ansatz ... 

Static charge-density susceptibilities have been computed ab initio by Li et al (38). The frequency-dependent susceptibility x(r, r cd) can be calculated within density functional theory, using methods developed by Ando (39 Zang-will and Soven (40 Gross and Kohn (4I and van Gisbergen, Snijders, and Baerends (42). In ab initio work, x(r, r co) can be determined by use of time-dependent perturbation techniques, pseudo-state methods (43-49), quantum Monte Carlo calculations (50-52), or by explicit construction of the linear response function in coupled cluster theory (53). Then the imaginary-frequency susceptibility can be obtained by analytic continuation from the susceptibility at real frequencies, or by a direct replacement co ico, where possible (for example, in pseudo-state expressions). 

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,... 

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... 

Our present focus is on density functional theory and coupled cluster methods for describing molecular systems interacting with a structured environment, and we focus on the derivation of linear response properties and compare the expressions that we obtain for the two different electronic structure methods. Based on linear response... 

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

---