Excited states coupled-clusters

Because of its size-extensivity and faster convergence with respect to excitation level Coupled cluster theory has replaced Cl theory as the dominant approach in ab initio correlation calculations. Like MBPT the theory is still mainly applied in cases where the exact wave function is dominated by a single determinant, but multireference methods have been formulated and begin to enter mainstream quantum chemistry. Generalization of the algorithms to the relativistic no-pair level can again be achieved via the spinorbital formulation of the methods. I will first discuss the single reference method and then consider the Fock space method [40] that uses multi-reference wavefiinctions for ionized or excited states. 

The last term of AG ° in O Eq. 28.52 introduces the constraint for the ground state coupled-cluster wavefunction, and contains the de-excitation operator Z given by... 

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. 

There is a variation on the coupled cluster method known as the symmetry adapted cluster (SAC) method. This is also a size consistent method. For excited states, a Cl out of this space, called a SAC-CI, is done. This improves the accuracy of electronic excited-state energies. 

Since the singly excited determinants effectively relax the orbitals in a CCSD calculation, non-canonical HF orbitals can also be used in coupled cluster methods. This allows for example the use of open-shell singlet states (which require two Slater determinants) as reference for a coupled cluster calculation. 

The use of Cl methods has been declining in recent years, to the profit of MP and especially CC methods. It is now recognized that size extensivity is important for obtaining accurate results. Excited states, however, are somewhat difficult to treat by perturbation or coupled cluster methods, and Cl or MCSCF based methods have been the prefen ed methods here. More recently propagator or equation of motion (Section 10.9) methods have been developed for coupled cluster wave functions, which allows calculation of exited state properties. 

Regarding the emission properties, AM I/Cl calculations, performed on a cluster containing three stilbene molecules separated by 4 A, show that the main lattice deformations take place on the central unit in the lowest excited state. It is therefore reasonable to assume that the wavefunction of the relaxed electron-hole pair extends at most over three interacting chains. The results further demonstrate that the weak coupling calculated between the ground state and the lowest excited state evolves in a way veiy similar to that reported for cofacial dimers. 

Figure 4-13. Evolution with the size of the sexithienyl cluster of the excitation energies from the ground stale to the lowest excited state (open circles), to the high-lying excited slate strongly coupled to the ground state (open squares), and to the lowest charge transfer-excited stale (open triangles). In all cases, only inlralayer interactions have been considered.

The two iron ions of the Rieske cluster are antiferromagnetically coupled therefore, the ground state has a spin S = while excited states of the spin ladder S = I, i, I, and, are at energies -3J, 8J,... 

The CC2 method [74] is an approximation to coupled cluster with singles and doubles (CCSD), and the excited state energies calculated have MP2 quality. An implementation that employs the resolution of identity (RI) approximation for two-electron integrals to reduce the CPU time is also available, RI-CC2 [75], which is suitable for large scale integral-direct calculations. This method has been implemented in TURBOMOLE [76],... 

Solvatochromic shifts for cytosine have also been calculated with a variety of methods (see Table 11-7). Shukla and Lesczynski [215] studied clusters of cytosine and three water molecules with CIS and TDDFT methods to obtain solvatochromic shifts. More sophisticated calculations have appeared recently. Valiev and Kowalski used a coupled cluster and classical molecular dynamics approach to calculate the solvatochromic shifts of the excited states of cytosine in the native DNA environment. Blancafort and coworkers [216] used a CASPT2 approach combined with the conductor version of the polarizable continuous (CPCM) model. All of these methods predict that the first three excited states are blue-shifted. S i, which is a nn state, is blue-shifted by 0.1-0.2 eV in water and 0.25 eV in native DNA. S2 and S3 are both rnt states and, as expected, the shift is bigger, 0.4-0.6eV for S2 and 0.3-0.8 eV for S3. S2 is predicted to be blue-shifted by 0.54 eV in native DNA. 

Stanton JF, Bartlett RJ (1993) The equation of motion coupled-cluster method - a systematic biorthogonal approach to molecular-excitation energies, transition-probabilities, and excited-state properties. J Chem Phys 98 7029... 

Krylov AI (2006) Spin-flip equation-of-motion coupled-cluster electronic structure method for a description of excited states, bond breaking, diradicals, and triradicals. Acc Chem Res 39 83-91... 

---