Independent amplitude approximations in coupled cluster valence bond theory: Incorporation of 3-electron-pair correlation and application to spin frustration in the low-lying excited states of a ferredoxin-type tetrametallic iron-sulfur cluster.

Coupled cluster valence bond (CCVB) is a simple electronic structure method based on a perfect pairing (PP) reference with 2-pair recouplings for strong electron correlation problems. CCVB is spin-pure, size-consistent, and can exactly (in its active space) separate any molecule into atoms for which unrestricted Hartree-Fock (UHF) at dissociation is the sum of the ground state UHF energies of the atoms. However CCVB is far from a complete description of strong correlations. Its first failure to exactly describe spin-recouplings arises at the level of 3 electron pairs, such as the recoupling of 3 triplet oxygen atoms in the dissociation of singlet ozone. Such situations are often associated with spin frustration. To address this limitation, an extension of CCVB, termed CCVB+i3, is reported here that includes an independent (i) amplitude approximation to the 3-pair recouplings. CCVB+i3 thereby has the same basic computational requirements as those of CCVB, which has previously been shown to be an efficient method. CCVB+i3 correctly separates molecules that CCVB cannot. As a by-product, an independent 2-pair amplitude approximation to CCVB, called PP+i2, is also defined. Remarkably, PP+i2 can also correctly separate systems that CCVB cannot. CCVB+i3 is validated on the symmetric dissociation of D _{3 h} ozone. CCVB+i3 is then used to explore the role of 3-pair recouplings in an [Fe ₄ S ₄ (SCH ₃ ) ₄ ] ^2- cluster that has been used to model the iron-sulfur core of [Fe ₄ S ₄ ] ferredoxins. Using localized PP orbitals, such recouplings are demonstrated to be large in some low-lying singlet excited states of the cluster. Significant 3 pair recoupling amplitudes include the usual triangular motif associated with spin frustration and other geometric arrangements of the 3 entangled pairs across the 4 iron centers.