Excitonically renormalised coupled-cluster theory

Excitonic Configuration Interaction: Going Beyond the Frenkel Exciton Model

We present the excitonic configuration interaction (ECI) method─a fragment-based analogue of the CI method for electronic structure calculations of multichromophoric systems. It can also be viewed as a generalization of the exciton approach, with the following properties: (i) It constructs the effective Hamiltonian exclusively from monomer calculations. (ii) It employs the strong orthogonality assumption and is exact within McWeeny's group function theory, thus requiring only one-electron density matrices of the monomer states. (iii) It is agnostic of the monomer electronic structure method, allowing us to use/combine different methods. (iv) It includes an embedding point charge scheme (called excitonic Hartree-Fock, EHF) to improve the accuracy of the monomer states, but such that the effective full-system Hamiltonian is not explicitly dependent on the embedding. (v) It is systematically improvable, by expanding the set of monomer states and by including configurations where two or more monomers are excited (defining the ECIS, ECISD, etc., methods). The performance of ECI is assessed by computing the absorption spectrum of two exemplary multichromophoric systems, using CIS as the monomer electronic structure method. The accuracy of ECI significantly depends on the chosen embedding charges and the ECI expansion. The most accurate assessed combinations─ECIS or ECISD with EHF embedding─yield spectra that agree qualitatively and quantitatively with full-system direct calculations, with deviations of the excitation energies below 0.1 eV. We also show that ECISD based on CIS monomer calculations can predict states where two monomers are excited simultaneously (e.g., triplet-triplet double-local excitations) that are inaccessible in a full-system CIS calculation.

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.

Cluster many-body expansion: A many-body expansion of the electron correlation energy about a cluster mean field reference

The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. In the past, applications of MBE to correlation energy have been unfeasible for large systems, but recent improvements to computing resources have sparked renewed interest in capturing the correlation energy using the generalized nth order Bethe-Goldstone equation. In this work, we extend this approach, originally proposed for a Slater determinant, to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean field reference TPS configuration and includes the correlation contribution of the excited TPSs using the MBE. This method, named cluster MBE (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from a RHF reference. We present numerical results for strongly correlated systems, such as the one- and two-dimensional Hubbard models and the chromium dimer. The performance of the cMBE method is also tested by partitioning the extended π space of several large π-conjugated systems, including a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 10⁶⁶ determinants for the exact FCI solution.