Ground State Properties, Gradients and Geometries

For CC2, one distinguishes between orbital-relaxed and unrelaxed properties. Both are calculated as first derivatives of the respective energy with respect to an external field corresponding to the calculated property. They differ in the treatment of the SCF orbitals. In the orbital-relaxed case the external field is (formally) already included at the SCF stage and the orbitals are allowed to relax in the external field; in the orbital-unrelaxed case the external field is first applied after the SCF calculation and the orbitals do not respond to the external field. Orbital-unrelaxed CC2 properties are calculated as first derivatives of the real part of the unrelaxed Lagrangian [82]

$$L^{\rm ur~CC2}(t,\bar{t},\beta) = \langle\mathrm{HF}\vert H\vert\mathrm{CC}\rangle + \sum_{\mu_1} \bar{t}_{\mu_1} \langle\mu_1\vert\hat{H} +[\hat{H},T_2]\vert\mathrm{HF}\rangle$$ (9.12)

$$+ \sum_{\mu_2} \bar{t}_{\mu_2} \langle\mu_2\vert\hat{H} +[F_0+\beta\hat{V},T_2]\vert\mathrm{HF}\rangle$$

with $H = H_0 + \beta V$ mathend000#--where $V$ mathend000# is the (one-electron) operator describing the external field, $\beta$ mathend000# the field strength, and $H_0$ mathend000# and $F_0$ mathend000# are the Hamiltonian and Fock operators of the unperturbed system--by the expression:

$$\langle V \rangle^{\rm ur~CC2} = \Re \left(\frac{ \partial L^{\rm ur~CC2}(t,\bar{t},\beta) }{\partial \beta}\right)_0 ~=~ \sum_{pq} D^{\rm ur}_{pq} ~V_{pq} ~~,$$ (9.13)

$$= \Re\bigg( \langle\mathrm{HF}\vert\hat{V}\vert\mathrm{HF}\rangle ... ...\mu_1} \bar{t}_{\mu_1} \langle\mu_1\vert\hat{V} +[V,T_2]\vert\mathrm{HF}\rangle$$ (9.14)

$$+ \sum_{\mu_2} \bar{t}_{\mu_2} \langle\mu_2\vert[\hat{V},T_2]\vert\mathrm{HF}\rangle \bigg)~~,$$

where $\Re$ mathend000# indicates that the real part is taken. Relaxed CC2 properties (and gradients) are calculated from the the full variational density including the contributions from the orbital response to the external perturbation, which are derived from the Lagrangian [88,13]

$$L^{\rm rel~CC2}(t,\bar{t}) = \langle\mathrm{HF}\vert H\vert\mathrm{CC}\rangle + \sum_{\mu_1} \bar{t}_{\mu_1} \langle\mu_1\vert\hat{H} +[\hat{H},T_2]\vert\mathrm{HF}\rangle$$ (9.15)

$$+ \sum_{\mu_2} \bar{t}_{\mu_2} \langle\mu_2\vert\hat{H} +[F,T_2]\vert\mathrm{HF}\rangle + \sum_{\mu_0} \bar{\kappa}_{\mu_0} F_{\mu_0} ~~,$$

where $F$ mathend000# is the Fock operator corresponding to the Hamiltonian of the perturbed system $H = H_0 + \beta V$ mathend000#. One-electron properties are then obtained as:

$$\langle V \rangle^{\rm rel~CC2} = \Re\bigg( \langle\mathrm{HF}\vert\hat{V}\vert\mathrm{HF}\rangle ... ...\mu_1} \bar{t}_{\mu_1} \langle\mu_1\vert\hat{V} +[V,T_2]\vert\mathrm{HF}\rangle$$ (9.16)

$$+ \sum_{\mu_2} \bar{t}_{\mu_2} \langle\mu_2\vert[V,T_2]\vert\mathrm{HF}\rangle + \sum_{\mu_0} \bar{\kappa}_{\mu_0} V_{\mu_0} \bigg)~~,$$

$$= \sum_{pq} D^{\rm rel}_{pq}~ V_{pq} ~~.$$ (9.17)

The calculation of one-electron first-order properties requires that in addition to the cluster equations also the linear equations for the Lagrangian multipliers $\bar{t}_\mu$ mathend000# are solved, which requires similar resources (CPU, disk space, and memory) as the calculation of a single excitation energy. For orbital-relaxed properties also a CPHF-like linear equation for the Lagrangian multipliers $\bar{\kappa}_{\mu_0}$ mathend000# needs to be solved and the two-electron density has to be build, since it is needed to set up the inhomogeneity (right-hand side). The calculation of relaxed properties is therefore somewhat more expensive--the operation count for solving the so-called Z-vector equations is similar to what is needed for an SCF calculation--and requires also more disk space to keep intermediates for the two-electron density--about $O(2V+2N)N_x + N_x^2$ mathend000# in addition to what is needed for the solution of the cluster equations. For ground states, orbital-relaxed first-order properties are standard in the literature.

The calculation of the gradient implies the calculation of the same variational densities as needed for relaxed one-electron properties and the solution of the same equations. The construction of the gradient contributions from the densities and derivative integrals takes about the same CPU time as 3-4 SCF iterations and only minor extra disk space. For details of the implementation of CC2 relaxed first-order properties and gradients and a discussion of applicability and trends of CC2 ground-state equilibrium geometries see ref. [13]. The following is in example input for a MP2 and CC2 single point calculation of first-order properties and gradients:

$ricc2
  mp2
  cc2
$response
  static relaxed operators=diplen,qudlen
  gradient

A different input is required for geometry optimizations: in this case the model for which the geometry should be optimized must be specified in the data group $ricc2 by the keyword geoopt:

$ricc2
  mp2
  cc2
  geoopt model=cc2

For CC2 calculations, the single-substitution part of the Lagrangian multipliers $\bar{t}_\mu$ mathend000# are saved in the file CCL0--1--1---0 and can be kept for a restart (for MP2 and CCS, the single-substitution part $\bar{t}_\mu$ mathend000# vanishes).

For MP2 only relaxed first-order properties and gradients are implemented (unrelaxed MP2 properties are defined differently than in CC response theory and are not implemented). For MP2, only the CPHF-like Z-vector equations for $\bar{\kappa}_{\mu_0}$ mathend000# need to be solved, no equations have to be solved for the Lagrangian multipliers $\bar{t}_\mu$ mathend000#. CPU time and disk space requirements are thus somewhat smaller than for CC2 properties or gradients.

For SCF/CIS/CCS it is recommended to use the modules GRAD and RDGRAD for the calculation of, ground state gradients and first-order properties.