The influence of orbital rotation on the energy of closed-shell wavefunctions

The orbital dependence of closed-shell wavefunction energies is investigated by performing doubly-occupied configuration interaction (DOCI) calculations, representing the most general class of these wavefunctions. Different local minima are examined for planar hydrogen clusters containing two, four, and six electrons applying (spin) symmetry-broken restricted, unrestricted, and generalised orbitals with real and complex coefficients. Contrary to Hartree–Fock (HF), restricted DOCI is found to properly break bonds and thus unrestricted orbitals, while providing a quantitative improvement of the energy, are not needed to enforce a qualitatively correct bond dissociation. For the beryllium atom and the BH diatomic, the lowest possible HF energy requests symmetry-broken generalised orbitals, whereas accurate results for DOCI can be obtained within a restricted formalism. Complex orbital coefficients are shown to increase the accuracy of HF and DOCI results in certain cases. The computationally inexpensive AP1roG geminal wavefunction is proven to agree very well with all DOCI results of this study.