7.11 Correlated Excited State Methods: The ADC(n) Family

7.11.2 IP- and EA-ADC

(May 16, 2021)

Similar to ADC scheme of the polarization propagator (ADC), the procedure can also be applied to the (N-1)- and (N+1)-parts of the electron propagator, which is given in its spectral representation as

Gpq(ω)=nΨ0N|cp|ΨnN+1ΨnN+1|cq|Ψ0Nω+E0N-EnN+1(N+1)-electron (EA) part+nΨ0N|cq|ΨnN-1ΨnN-1|cp|Ψ0Nω+EnN-1-E0N(N-1)-electron (IP) part. (7.89)

Doing so, the (non-Dyson) IP- and EA-ADC methods up to third order of perturbation theory have been derived. 967 Schirmer J., Trofimov A. B., Stelter G.
J. Chem. Phys.
(1998), 109, pp. 4734–4744.
Link
, 1102 Trofimov A. B., Schirmer J.
J. Chem. Phys.
(2005), 123, pp. 144115.
Link
As in the case of the ADC scheme of the polarization propagator, the derivation of the same working equations is possible via the ISR formalism, the only difference to the procedure shown in the previous section 7.11.1 being the excitation operators

{C^JN-1} ={ci;cacicj,i<j;} IP-ADC (7.90)
{C^JN+1} ={ca;cacbci,a<b;} EA-ADC (7.91)

replacing the electron number-conserving one C^J in Eq. (7.82).

Diagonalization of the IP- and EA-ADC secular matrices 𝐌 yields electron-detachment energies (or ionization potentials, IPs) and electron-attachment energies (or negative electron affinities, EAs), respectively. In addition, relative spectral intensities of (N-1)- and (N+1)-transitions are accessible as pole strengths Pn, which are computed according to

Pn=p|xpn|2, (7.92)

where the xpn are the spectroscopic factors computed by means of the IP- and EA-ADC eigenvectors 𝐗 and the matrix of effective transition amplitudes 𝐟 using the relations

xpnN-1 =JXJnfJpN-1=JXJnΨ~nN-1|cp|Ψ0N IP-ADC (7.93)
xpnN+1 =JXJnfJpN+1=JXJnΨ~JN+1|cp|Ψ0N EA-ADC. (7.94)

When requesting electron-detached or electron-attached states, the pole strengths are automatically computed. For IP- and EA-ADC(2) calculations, second-order pole strengths are used [i.e., the IP- and EA-ADC(2) 𝐟 matrix is employed in their computation]. As suggested in Ref. 967, for computational reasons the same second-order pole strengths are computed in case of strict IP- and EA-ADC(3) calculations, i.e., using the second-order ground state density throughout the 𝐌 matrix equations, ADC_DENSITY_ORDER = 2. When requesting a higher-order ground state density to be used, e.g., by setting ADC_DENSITY_ORDER = 3 (corresponding to the Σ(4) scheme 1102 Trofimov A. B., Schirmer J.
J. Chem. Phys.
(2005), 123, pp. 144115.
Link
) or ADC_DENSITY_ORDER = 4 (corresponding to the Σ(4+) scheme, 1102 Trofimov A. B., Schirmer J.
J. Chem. Phys.
(2005), 123, pp. 144115.
Link
also denoted as standard IP- and EA-ADC(3) schemes), the corresponding pole strengths are used, i.e., third-order pole strenghts in case of Σ(4) and improved third-order pole strengths in case of Σ(4+).

The spectroscopic factors also allow for computing the Dyson orbitals |ϕn connected to electron detachment and attachment processes

|ϕn=pxpn|φp, (7.95)

where in the latter relation the φp denote HF orbitals. Dyson orbital output is triggered by ADC_DO_DYSON = TRUE together with STATE_ANALYSIS = TRUE. Also see Section 10.2.6 for further details.