A minimalistic approach to static and dynamic electron correlations: Amending generalized valence bond method with extended random phase approximation correlation correction

Geminal-Based Strategies for Modeling Large Building Blocks of Organic Electronic Materials

We elaborate on unconventional electronic structure methods based on geminals and their potential to advance the rapidly developing field of organic photovoltaics (OPVs). Specifically, we focus on the computational advantages of geminal-based methods over standard approaches and identify the critical aspects of OPV development. Examples are reliable and efficient computations of orbital energies, electronic spectra, and van der Waals interactions. Geminal-based models can also be combined with quantum embedding techniques and a quantum information analysis of orbital interactions to gain a fundamental understanding of the electronic structures and properties of realistic OPV building blocks. Furthermore, other organic components present in, for instance, dye-sensitized solar cells (DSSCs) represent another promising scope of application. Finally, we provide numerical examples predicting the properties of a small building block of OPV components and two carbazole-based dyes proposed as possible DSSC sensitizers.

Benchmarking ionization potentials using the simple pCCD model

The electron-detachment energy is measured by the ionization potential (IP). As a result, it is a fundamental, observable and important molecular electronic signature in photoelectron spectroscopy. A precise theoretical prediction of electron-detachment energies or ionization potentials is essential for organic optoelectronic systems like transistors, solar cells, or light-emitting diodes. In this work, we benchmark the performance of the recently presented IP variant of the equation-of-motion pair coupled cluster doubles (IP-EOM-pCCD) model to determine IPs. Specifically, the predicted ionization energies are compared to experimental results and higher-order coupled cluster theories based on statistically assessing 201 electron-detached states of 41 organic molecules for three different molecular orbital basis sets and two sets of particle-hole operators. While IP-EOM-pCCD features a reasonable spread and skewness of ionization energies, its mean error and standard deviation differ by up to 1.5 eV from reference data. Our study, thus, highlights the importance of dynamical correlation to reliably predict IPs from a pCCD reference function in small organic molecules.

Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry

In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.

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.

Systematically Improvable Generalization of Self-Interaction Corrected Density Functional Theory

Perdew-Zunger self-interaction correction (PZSIC) reintroduces an exact constraint to approximate density functional theory (DFT). However, PZSIC can paradoxically degrade performance, and standard DFT approximations (with or without PZSIC) are not systematically improvable. We use the adiabatic projection formalism (Janesko, B. G. J. Chem. Phys. 2022, 156, 014111, https://doi.org/10.1063/5.0076144) to derive PZSIC in terms of a reference system experiencing only electron self-interaction. Generalization to a "self-and-some-others" interaction introduces correlation into the reference system, systematically bridging from PZSIC to exact wave function theory without the double counting of correlation. Minimal active spaces accurately treat nearly one electron, near-equilibrium, and strongly correlated model systems at modest computational expense.

An efficient implementation of time-dependent linear-response theory for strongly orthogonal geminal wave function models

We present an implementation of time-dependent linear-response equations for strongly orthogonal geminal wave function models: the generalized valence bond perfect-pairing theory (TD-GVB) and the antisymmetrized product of strongly orthogonal geminals (TD-APSG). The geminal wave functions are optimized using a restricted-step second-order algorithm suitable for handling many geminals, and the linear-response equations are solved in an efficient way using a direct iterative approach. The wave function optimization algorithm features an original scheme to create initial orbitals for the geminal functions in a black-box fashion. The implementation is employed to examine the accuracy of the geminal linear response for singlet excitation energies of small and medium-sized molecules. In systems dominated by dynamic correlation, geminal models constitute only a minor improvement with respect to time-dependent Hartree-Fock. Compared to linear-response complete active space self-consistent field (LR-CASSCF), TD-GVB either misses or gives large errors for states dominated by double excitations.

Dispersion Interactions between Molecules in and out of Equilibrium Geometry: Visualization and Analysis

The London dispersion interactions between systems undergoing bond breaking, twisting, or compression are not well studied due to the scarcity and the high computational cost of methods being able to describe both the dynamic correlation and the multireference character of the system. Recently developed methods based on the Generalized Valence Bond wave function, such as EERPA-GVB and SAPT(GVB) (SAPT = symmetry-adapted perturbation theory), allow one to accurately compute and analyze noncovalent interactions between multireference systems. Here, we augment this analysis by introducing a local indicator for dispersion interactions inspired by Mata and Wuttke's Dispersion Interaction Density [ J. Comput. Chem. 2017, 38, 15-23] applied on top of an EERPA-GVB computation. Using a few model systems, we show what insights into the nature and evolution of the dispersion interaction during bond breaking and twisting such an approach is able to offer. The new indicator can be used at a minimal cost additional to an EERPA-GVB computation and can be complemented by an energy decomposition employing the SAPT(GVB) method. We explain the physics behind the initial increase, followed by a decrease in the interaction of linear molecules upon bond stretching. Namely, the elongation of covalent bonds leads to the enhancement of attractive dispersion interactions. For even larger bond lengths, this effect is canceled by the increase of the repulsive exchange forces resulting in a suppression of the interaction and finally leading to repulsion between monomers.

Ring coupled cluster doubles at the multireference level

A ring approximation within an internally contracted multireference (MR) Coupled Cluster (CC) framework is worked out and tested. Derivation of equations utilizes MR based, generalized normal ordering and the corresponding generalized Wick-theorem (MR-GWT). Contractions among cluster operators are avoided by adopting a normal ordered exponential ansatz. The original version of the MR ring CC doubles (MR-rCCD) equations [Á. Szabados and Á. Margócsy, Mol. Phys. 115, 2731 (2017)] is rectified in two aspects. On the one hand, over-completeness of double excitations is treated by relying on the concept of frames. On the other hand, restriction on the maximal cumulant rank is lifted from two to four. This is found essential for obtaining reliable correlation corrections to the energy. The MR function underlying the approach is provided by the Generalized Valence Bond (GVB) model. The pair structure of the reference ensures a fragment structure of GVB cumulants. This represents a benefit when evaluating cumulant contractions appearing as a consequence of MR-GWT. In particular, cumulant involving terms remain less expensive than their traditional, pair-contracted counterpart, facilitating an O(N⁶) eventual scaling of the proposed MR-rCCD method. Pilot applications are presented for covalent bond breaking, deprotonation energies, and torsional potentials.

Second-Order Dispersion Energy Based on Multireference Description of Monomers

We propose a method for calculating a second-order dispersion energy for weakly interacting multireference systems in arbitrary electronic states. It is based on response properties obtained from extended random phase approximation equations. The introduced formalism is general and requires only one- and two-particle reduced density matrices of monomers. We combine the new method with either generalized valence bond perfect pairing (GVB) or complete active space (CAS) self-consistent field description of the interacting systems. In addition to a general scheme, three approximations, leading to significant reduction of the computational cost, are developed by exploiting Dyall partitioning of the monomer Hamiltonians. For model multireference systems (H₂···H₂ and Be···Be) the method is accurate, unlike its single-reference-based counterpart. Neither GVB nor CAS description of single-reference monomers improves the dispersion energy with respect to the Hartree-Fock-based results.

Automatic Construction of the Initial Orbitals for Efficient Generalized Valence Bond Calculations of Large Systems

We propose an efficient general strategy for generating initial orbitals for generalized valence bond (GVB) calculations which makes routine black-box GVB calculations on large systems feasible. Two schemes are proposed, depending on whether the restricted Hartree-Fock (RHF) wave function is stable (scheme I) or not (scheme II). In both schemes, the first step is the construction of active occupied orbitals and active virtual orbitals. In scheme I, active occupied orbitals are composed of the valence orbitals (the inner core orbitals are excluded), and the active virtual orbitals are obtained from the original virtual space by requiring its maximum overlap with the virtual orbital space of the same system at a minimal basis set. In scheme II, active occupied orbitals and active virtual orbitals are obtained from the set of unrestricted natural orbitals (UNOs), which are transformed from two sets of unrestricted HF spatial orbitals. In the next step, the active occupied orbitals and active virtual ones are separately transformed to localized orbitals. Localized occupied and virtual orbital pairs are formed using the Kuhn-Munkres (KM) algorithm and are used as the initial guess for the GVB orbitals. The optimized GVB wave function is obtained using the second-order self-consistent-field algorithm in the GAMESS program. With this procedure, GVB energies have been obtained for the lowest singlet and triplet states of polyacenes (up to decacene with 96 pairs) and the singlet ground state of two di-copper-oxygen-ammonia complexes. We have also calculated the singlet-triplet gaps for some polyacenes and the relative energy between two di-copper-oxygen-ammonia complexes with the block-correlated second-order perturbation theory based on the GVB reference.

Electron Correlation from the Adiabatic Connection for Multireference Wave Functions

An adiabatic connection (AC) formula for the electron correlation energy is derived for a broad class of multireference wave functions. The AC expression recovers dynamic correlation energy and assures a balanced treatment of the correlation energy. Coupling the AC formalism with the extended random phase approximation allows one to find the correlation energy only from reference one- and two-electron reduced density matrices. If the generalized valence bond perfect pairing model is employed a simple closed-form expression for the approximate AC formula is obtained. This results in the overall M^{5} scaling of the computation cost making the method one of the most efficient multireference approaches accounting for dynamic electron correlation also for the strongly correlated systems.