Coupled cluster

Coupled cluster method is a technique used for description of the many-body systems. The method was initially developed by Fritz Coester and Hermann Kümmel in 1950's for studying nuclear physics phenomena but it became more frequently used after Jiři Čížek and Josef Paldus reformulated the method for studying electronic correlation in atoms and molecules in 1960's. It is now one of the most prevalent methods in quantum chemistry that include electronic correlation. The method is based on exponential Ansatz:

$\vert{\Psi}\rangle = e^\hat{T} \vert{\Phi_0}\rangle$ ,

where $\vert{\Psi}\rangle$ is the wave function, $\vert{\Phi_0}\rangle$ is the reference function (e.g. Hartree-Fock function), and $\hat{T}$ is the cluster operator:

$\hat{T}=\hat{T}_1 + \hat{T}_2 + \cdots$ ,

where in the formalism of second quantization:

$\hat{T}_1=\sum_{i}\sum_{a} t_{i}^{a} \hat{a}_{i}\hat{a}^{\dagger}_{a},$

$\hat{T}_2=\frac{1}{4}\sum_{i,j}\sum_{a,b} t_{ij}^{ab} \hat{a}_{i}\hat{a}_j\hat{a}^{\dagger}_{a}\hat{a}^{\dagger}_{b}.$

In the above formulae $\hat{a}$ and $\hat{a}^{\dagger}$ denote the creation and annihilation operators and i,j stand for occupied and a,b for unoccupied orbitals. The coupled cluster equations are usually derived using diagrammatic technique and result in nonlinear equations which can be solved in an iterative way.

In the simplest version one considers only $\hat{T}_2$ operator (double excitations). This method is called coupled cluster with doubles (CCD in short).

The method gives exact non-relativistic solution of the Schrödinger equation of the n-body problem if one includes up to the $\hat{T}_n$ cluster operator. However, the computational effort of solving the equations grows steeply with the order of the cluster operator and in practical applications the method is limited to the first few orders.

Categories: Quantum mechanics