Connections between many-body perturbation and coupled-cluster theories

Here, we build on the works of Scuseria et al. [J. Chem. Phys. 129, 231101 (2008)] and Berkelbach [J. Chem. Phys. 149, 041103 (2018)] to show connections between the Bethe-Salpeter equation (BSE) formalism combined with the GW approximation from many-body perturbation theory and coupled-cluster (CC) theory at the ground- and excited-state levels. In particular, we show how to recast the GW and Bethe-Salpeter equations as non-linear CC-like equations. Similitudes between BSE@GW and the similarity-transformed equation-of-motion CC method are also put forward. The present work allows us to easily transfer key developments and the general knowledge gathered in CC theory to many-body perturbation theory. In particular, it may provide a path for the computation of ground- and excited-state properties (such as nuclear gradients) within the GW and BSE frameworks.

Scrutinizing GW-Based Methods Using the Hubbard Dimer

Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.

Modeling core-level excitations with variationally optimized reduced-density matrices and the extended random phase approximation

The information contained within ground-state one- and two-electron reduced-density matrices (RDMs) can be used to compute wave functions and energies for electronically excited states through the extended random phase approximation (ERPA). The ERPA is an appealing framework for describing excitations out of states obtained via the variational optimization of the two-electron RDM (2-RDM), as the variational 2-RDM (v2RDM) approach itself can only be used to describe the lowest-energy state of a given spin symmetry. The utility of the ERPA for predicting near-edge features relevant to x-ray absorption spectroscopy is assessed for the case that the 2-RDM is obtained from a ground-state v2RDM-driven complete active space self-consistent field (CASSCF) computation. A class of killer conditions for the CASSCF-specific ERPA excitation operator is derived, and it is demonstrated that a reliable description of core-level excitations requires an excitation operator that fulfills these conditions; the core-valence separation (CVS) scheme yields such an operator. Absolute excitation energies evaluated within the CASSCF/CVS-ERPA framework are slightly more accurate than those obtained from the usual random phase approximation (RPA), but the CVS-ERPA is not more accurate than RPA for predicting the relative positions of near-edge features. Nonetheless, CVS-ERPA is established as a reasonable starting point for the treatment of core-level excitations using variationally optimized 2-RDMs.

Restricted second random phase approximations and Tamm-Dancoff approximations for electronic excitation energy calculations

In this article, we develop systematically second random phase approximations (RPA) and Tamm-Dancoff approximations (TDA) of particle-hole and particle-particle channels for calculating molecular excitation energies. The second particle-hole RPA/TDA can capture double excitations missed by the particle-hole RPA/TDA and time-dependent density-functional theory (TDDFT), while the second particle-particle RPA/TDA recovers non-highest-occupied-molecular-orbital excitations missed by the particle-particle RPA/TDA. With proper orbital restrictions, these restricted second RPAs and TDAs have a formal scaling of only O(N(4)). The restricted versions of second RPAs and TDAs are tested with various small molecules to show some positive results. Data suggest that the restricted second particle-hole TDA (r2ph-TDA) has the best overall performance with a correlation coefficient similar to TDDFT, but with a larger negative bias. The negative bias of the r2ph-TDA may be induced by the unaccounted ground state correlation energy to be investigated further. Overall, the r2ph-TDA is recommended to study systems with both single and some low-lying double excitations with a moderate accuracy. Some expressions on excited state property evaluations, such as ⟨Ŝ(2)⟩ are also developed and tested.

Equivalence of particle-particle random phase approximation correlation energy and ladder-coupled-cluster doubles

The recent proposal to determine the (exact) correlation energy based on pairing matrix fluctuations by van Aggelen et al. ["Exchange-correlation energy from pairing matrix fluctuation and the particle-particle random phase approximation," preprint arXiv:1306.4957 (2013)] revived the interest in the simplest approximation along this path: the particle-particle random phase approximation (pp-RPA). In this paper, we present an analytical connection and numerical demonstrations of the equivalence of the correlation energy from pp-RPA and ladder-coupled-cluster doubles. These two theories reduce to identical algebraic matrix equations and correlation energy expressions. The numerical examples illustrate that the correlation energy missed by pp-RPA in comparison with coupled-cluster singles and doubles is largely canceled out when considering reaction energies. This theoretical connection will be beneficial to design density functionals with strong ties to coupled-cluster theories and to study molecular properties at the pp-RPA level relying on well established coupled cluster techniques.