Coupled Electron Pair-Type Approximations for Tensor Product State Wave Functions

Size extensivity, defined as the correct scaling of energy with system size, is a desirable property for any many-body method. Traditional configuration interaction (CI) methods are not size extensive, hence the error increases as the system gets larger. Coupled electron pair approximation (CEPA) methods can be constructed as simple extensions of a truncated CI that ensures size extensivity. One of the major issues with the CEPA and its variants is that singularities arise in the amplitude equations when the system starts to be strongly correlated. In this work, we extend the traditional Slater determinant based coupled electron pair approaches like CEPA-0, averaged coupled-pair functional, and average quadratic coupled-cluster to a new formulation based on tensor product states (TPS). We show that a TPS basis can often be chosen such that it removes the singularities that commonly destroy the accuracy of CEPA based methods. A suitable TPS representation can be formed by partitioning the system into separate disjoint clusters and forming the final wave function as the tensor product of the many body states of these clusters. We demonstrate the application of these methods on simple bond breaking systems such as CH₄ and F₂ where determinant based CEPA methods fail. We further apply the TPS-CEPA approach to stillbene isomerization and few planar π-conjugated systems. Overall, the results show that the TPS-CEPA method can remove the singularities and provide improved numerical results compared to common electronic structure methods.