Increasing spin-flips and decreasing cost: Perturbative corrections for external singles to the complete active space spin flip model for low-lying excited states and strong correlation

An approximation to the spin-flip extended configuration interaction singles method is developed using a second-order perturbation theory approach. In addition to providing significant efficiency advantages, the new framework is general for an arbitrary number of spin-flips, with the current implementation being applicable for up to around 4 spin-flips. Two new methods are introduced: one which is developed using non-degenerate perturbation theory, spin-flip complete active-space (SF-CAS(S)), and a second quasidegenerate perturbation theory method, SF-CAS(S)1. These two approaches take the SF-CAS wavefunction as the reference, and then perturbatively includes the effect of single excitations. For the quasidegenerate perturbation theory method, SF-CAS(S)1, the subscripted "1" in the acronym indicates that a truncated denominator expansion is used to obtain an energy-independent down-folded Hamiltonian. We also show how this can alternatively be formulated in terms of an extended Lagrangian, by introducing an orthonormality constraint on the first-order wavefunction. Several numerical examples are provided, which demonstrate the ability of SF-CAS(S) and SF-CAS(S)1 to describe bond dissociations, singlet-triplet gaps of organic molecules, and exchange coupling parameters for binuclear transition metal complexes.