next up previous contents index
Next: Running excitation energy calculations: Up: Second-Order Approximate Coupled-Cluster (CC2) Previous: Diagnostics:   Contents   Index

Calculation of Excitation Energies

With the ricc2 program excitation energies can at present be calculated with the RI variants of the methods CCS/CIS, CIS(D), CIS(D), ADC(2) and CC2. The CC2 excitation energies are obtained by standard coupled-cluster linear response theory as eigenvalues of the Jacobian, defined as derivative of the vector function with respect to the cluster amplitudes.

$\displaystyle \bf A^{{\rm CC2}}_{{\mu\nu}}$ = $\displaystyle {\frac{{d \Omega_{\mu}}}{{d t_{\nu}}}}$ = $\displaystyle \left(\vphantom{ \begin{array}{cc} \langle \mu_{1} \vert [[\hat{H...
...le \mu_{2} \vert [F,\tau_{\nu_2}] \vert\mathrm{HF} \rangle \end{array} }\right.$$\displaystyle \begin{array}{cc} \langle \mu_{1} \vert [[\hat{H}+[\hat{H},T_2],\...
...e & \langle \mu_{2} \vert [F,\tau_{\nu_2}] \vert\mathrm{HF} \rangle \end{array}$$\displaystyle \left.\vphantom{ \begin{array}{cc} \langle \mu_{1} \vert [[\hat{H...
...le \mu_{2} \vert [F,\tau_{\nu_2}] \vert\mathrm{HF} \rangle \end{array} }\right)$ (9.8)

Since the CC2 Jacobian is a non-symmetric matrix, left and right eigenvectors are different and the right (left) eigenvectors Eiν ( $ \bar{{E}}^{i}_{{\mu}}$) are not orthogonal among themselves, but form a biorthonormal basis (if properly normalized):

$\displaystyle \bar{{E}}^{{i}}_{}$Ej = $\displaystyle \bar{{E}}^{i}_{{\mu_1}}$Ejν1 + $\displaystyle \bar{{E}}^{i}_{{\mu_2}}$Ejν2 = δij  . (9.9)

To obtain excitation energies only the right or the left eigenvalue problem needs to be solved, but for the calculation of transition strengths and first-order properties both, left and right, eigenvectors are needed (see below). A second complication that arises from the non-symmetric eigenvalue problem is that in the case of close degeneracies within the same irreducible representation (symmetry) it can happen that instead of two close lying real roots a degenerate complex conjugated pair of excitation energies and eigenvectors is obtained. CC2 (and also other standard coupled-cluster response methods) are thus not suited for the description of conical intersections etc. For the general theory behind coupled cluster response calculations see e.g. ref. [104,105] or other reviews.

The ricc2 program exploits that the doubles/doubles block of the CC2 Jacobian is diagonal and the (linear) eigenvalue problem in the singles and doubles space can be reformulated as a (non-linear) eigenvalue problem in single-substitution space only:

$\displaystyle \bf A^{{eff}}_{{\mu_1\nu_1}}$(t, ω) = $\displaystyle \bf A^{{\rm CC2}}_{{\mu_1\nu_1}}$(t) - $\displaystyle \bf A^{{\rm CC2}}_{{\mu_1\gamma_2}}$(t)($\displaystyle \bf A_{{\gamma_2\gamma_2}}^{}$ - ω)$\displaystyle \bf A^{{\rm CC2}}_{{\gamma2\nu_1}}$(t)

$\displaystyle \bf A^{{eff}}_{{\mu_1\nu_1}}$(tCC2, ωCC2)Eν1 = ωCC2Eν1

This allows to avoid the storage of the double-substitution part of the eigen- or excitation vectors Eν2, $ \bar{{E}}_{{\nu_2}}^{}$. The algorithms are described in refs. [10,11], about the RI error see ref. [103].

The solution of the CC2 eigenvalue problem can be started from the solutions of the CCS eigenvalue problem (see below) or the trial vectors or solutions of a previous CC2 excitation energy calculation. The operation count per transformed trial vector for one iteration for the CC2 eigenvalue problem is about 1.3-1.7 times the operation count for one iteration for the cluster equations in the ground-state calculation--depending on the number of vectors transformed simultaneously. The disk space requirements are about O(V + N)Nx double precision words per vector in addition to the disk space required for the ground state calculation.

CCS excitation energies are obtained by the same approach, but here double-substitutions are excluded from the expansion of the excitation or eigenvectors and the ground-state amplitudes are zero. Therefore the CCS Jacobian,

$\displaystyle \bf A^{{\rm CCS}}_{{\mu\nu}}$ = $\displaystyle {\frac{{d \Omega_{\mu}}}{{d t_{\nu}}}}$ = 〈μ1|[H, τν1]| HF〉  , (9.10)

is a symmetric matrix and left and right eigenvectors are identical and form an orthonormal basis. The configuration interaction singles (CIS) excitation energies are identical to the CCS excitation energies. The operation count for a RI-CIS calculation is $ \cal {O}$(ON2Nx) per iteration and transformed trial vector.

The second-order perturbative correction CIS(D) to the CIS excitation energies is calculated from the expression

ωCIS(D) = ωCIS + ω(D) = $\displaystyle \bf E^{{\rm CIS}}_{}$$\displaystyle \bf A^{{eff}}_{}$(tMP1, ωCIS)$\displaystyle \bf E^{{\rm CIS}}_{}$ (9.11)

(Note that tMP1 are the first-order double-substitution amplitudes from which also the MP2 ground-state energy is calculated; the first-order single-substitution amplitudes vanish for a Hartree-Fock reference due to the Brillouin theorem.) The operation count for a RI-CIS(D) calculation is similar to that of a single iteration for the CC2 eigenvalue problem. Also disk space requirements are similar.

next up previous contents index
Next: Running excitation energy calculations: Up: Second-Order Approximate Coupled-Cluster (CC2) Previous: Diagnostics:   Contents   Index