Degenerate coupled-cluster theory

A size-extensive, converging, black-box, ab initio coupled-cluster ($Δ$CC) ansatz is introduced that computes the energies and wave functions of stationary states from any degenerate or nondegenerate Slater-determinant references with any numbers of $α$- and $β$-spin electrons, any patterns of orbital occupancy, any spin multiplicities, and any spatial symmetries. For a nondegenerate reference, it reduces to the single-reference coupled-cluster ansatz. For a degenerate multireference, it is a natural coupled-cluster extension of degenerate Rayleigh-Schrödinger perturbation ($Δ$MP) theory. For ionized and electron-attached references, it can be viewed as a coupled-cluster Green's function, although the present theory is convergent toward the full-configuration-interaction (FCI) limits, while Feynman-Dyson many-body Green's function (MBGF) theory generally is not. Additionally, a new state-universal multireference coupled-cluster theory for general model spaces is developed by slightly modifying the $Δ$CC ansatz. This quasidegenerate coupled-cluster (QCC) theory is size-extensive, converging, but not black-box, which is expected to be well suited for strong correlation. Determinant-based, general-order algorithms of $Δ$CC and QCC theories are implemented, which are compared with configuration-interaction (CI) and equation-of-motion coupled-cluster (EOM-CC) theories through octuple excitations and with $Δ$MP and MBGF theories up to the nineteenth order. For transition energies, the order of performance is: QCC $\approx$ $Δ$CC $>$ EOM-CC $>$ CI at the same excitation order or QCC $\approx$ $Δ$CC $>$ $Δ$MP $>$ MBGF at the same cost scaling.