10. Coupled cluster, size extensivity and diagnostics
Position in the course: Lesson 10 of 10. Complete the preceding derivation and use the explained exercises to check understanding.
Prerequisites: Second quantization and electron correlation
1. Exponentiate connected excitation operators
Coupled cluster uses \(|\Psi\rangle=e^T|\Phi_0\rangle\), where \(T=T_1+T_2+\cdots\) contains excitation operators relative to a reference determinant. CCSD retains connected singles and doubles in \(T\). Expanding the exponential gives
Consequently CCSD contains higher determinant excitation ranks through disconnected products. It is not identical to CISD and not equivalent to full triple or quadruple amplitudes. Connected and disconnected refer to the algebraic organization of contributions, not to literal fragments floating apart.
2. Derive the projected equations
Insert the ansatz into \(He^T|\Phi_0\rangle=Ee^T|\Phi_0\rangle\) and left-multiply by \(e^{-T}\). Define \(\bar H=e^{-T}He^T\). Projection gives
The Baker-Campbell-Hausdorff expansion is \(\bar H=H+[H,T]+\tfrac12[[H,T],T]+\cdots\). For the standard finite-basis two-body Hamiltonian and pure excitation operators, the connected expansion terminates at finite commutator order. The resulting nonlinear amplitude equations require iteration; a small energy change alone does not establish small amplitude residuals.
3. Understand size extensivity by factorization
For noninteracting fragments with factorizing reference and \(T=T_A+T_B\), the operators commute. Thus
The disconnected cross-products needed for simultaneous fragment correlation appear automatically. The energy is additive under the corresponding assumptions. This property matters for comparing molecules of different sizes and for dissociation limits. It does not ensure that a chosen single reference is appropriate for each fragment or for the entire bond-breaking path.
4. Worked two-fragment comparison
Let each fragment's approximate state be \((1+tD_A)|\Phi_A\rangle\) or \((1+tD_B)|\Phi_B\rangle\), where \(D\) makes a double substitution. Their product is \(1+tD_A+tD_B+t^2D_AD_B\). A global CISD truncation omits the last term. A CC doubles exponential generates it from \(\tfrac12(T_A+T_B)^2\), because the two commuting cross terms cancel the one-half. This explicit algebra explains the structural advantage without claiming any numerical molecular benchmark was run.
5. Accuracy, triples and response
CCSD(T) adds a particular perturbative triples correction to CCSD; it is not full iterative CCSDT. Conventional formal costs increase from roughly sixth power for CCSD to seventh for the triples correction. Strong near-degeneracy, poorly behaved amplitudes or an unstable reference can invalidate the perturbative triples assumption. Inspect reference character and diagnostics, but use no universal diagnostic cutoff as proof of accuracy.
The similarity-transformed Hamiltonian is generally non-Hermitian. Standard projected CC energies are not variational upper bounds. Property derivatives require the appropriate left-state or Lagrangian response, rather than simply treating the reference bra as the conjugate of the correlated ket. Excited-state equation-of-motion variants require their own truncation and state-character checks.
6. Exercises and explained answers
Does \(T_2^2/2\) describe a connected quadruple amplitude? No. It represents a product of connected doubles. Full connected quadruple flexibility requires \(T_4\).
Can a CC energy below a benchmark be accepted because lower is better? No. Projected CC lacks the variational ordering. Compare systematic truncations, basis refinements and reference suitability.
Why report residual norms? Nonlinear iterations can show a nearly unchanged energy while amplitudes remain nonstationary. Residuals test the actual projected equations.
Analytical teaching schematic, not simulation data.
Further conceptual check
Correlation amplitudes depend on orbital choices, whereas a properly converged observable must be interpreted with the whole approximation. Rotating occupied and virtual subspaces can alter the numerical organization without adding the missing connected excitation ranks. Local-correlation thresholds require convergence tests in addition to the canonical truncation choice.
7. References and study connections
- DePrince group academic HF tutorial
- DePrince group academic tutorial collection
- Psi4 official coupled-cluster documentation
- Helgaker, Jørgensen and Olsen: Molecular Electronic-Structure Theory
The derivations and toy arithmetic are original teaching constructions. No molecular simulation is reported here.
8. Related theory and practice
Quantum mechanics · Density functional theory · Quantum Monte Carlo · Gaussian