GW Vertex Corrected Calculations for Molecular Systems

Reference Energies for Valence Ionizations and Satellite Transitions

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.

Going Beyond the GW Approximation Using the Time-Dependent Hartree-Fock Vertex

The time-dependent Hartree-Fock (TDHF) vertex of many-body perturbation theory (MBPT) makes it possible to extend TDHF theory to charged excitations. Here we assess its performance by applying it to spherical atoms in their neutral electronic configuration. On a theoretical level, we recast the TDHF vertex as a reducible vertex, highlighting the emergence of a self-energy expansion purely in orders of the bare Coulomb interaction; then, on a numerical level, we present results for polarizabilities, ionization energies (IEs), and photoemission satellites. We confirm the superiority of THDF over simpler methods such as the random phase approximation for the prediction of atomic polarizabilities. We then find that the TDHF vertex reliably provides better IEs than GW and low-order self-energies do in the light-atom, few-electron regime; its performance degrades in heavier, many-electron atoms instead, where an expansion in orders of an unscreened Coulomb interaction becomes less justified. New relevant features are introduced in the satellite spectrum by the TDHF vertex, but the experimental spectra are not fully reproduced due to a missing account of nonlinear effects connected to hole relaxation. We also explore various truncations of the self-energy given by the TDHF vertex, but do not find them to be more convenient than low-order approximations such as GW and second Born (2B), suggesting that vertex corrections should be carried out consistently both in the self-energy and in the polarizability.

Fully Dynamic G3W2 Self-Energy for Finite Systems: Formulas and Benchmark

Over the years, Hedin's GW self-energy has been proven to be a rather accurate and simple approximation to evaluate electronic quasiparticle energies in solids and in molecules. Attempts to improve over the simple GW approximation, the so-called vertex corrections, have been constantly proposed in the literature. Here, we derive, analyze, and benchmark the complete second-order term in the screened Coulomb interaction W for finite systems. This self-energy named G3W2 contains all the possible time orderings that combine 3 Green's functions G and 2 dynamic W. We present the analytic formula and its imaginary frequency counterpart, with the latter allowing us to treat larger molecules. The accuracy of the G3W2 self-energy is evaluated on well-established benchmarks (GW100, Acceptor 24, and Core 65) for valence and core quasiparticle energies. Its link with the simpler static approximation, named SOSEX for static screened second-order exchange, is analyzed, which leads us to propose a more consistent approximation named 2SOSEX. In the end, we find that neither the G3W2 self-energy nor any of the investigated approximations to it improve over one-shot G0W0 with a good starting point. Only quasi-particle self-consistent GW HOMO energies are slightly improved by addition of the G3W2 self-energy correction. We show that this is due to the self-consistent update of the screened Coulomb interaction, leading to an overall sign change of the vertex correction to the frontier quasiparticle energies.

Comparing Self-Consistent GW and Vertex-Corrected G0W0 (G0W0Γ) Accuracy for Molecular Ionization Potentials

We test the performance of self-consistent GW and several representative implementations of vertex-corrected G0W0 (G0W0Γ). These approaches are tested on benchmark data sets covering full valence spectra (first ionization potentials and some inner valence shell excitations). For small molecules, when comparing against state-of-the-art wave function techniques, our results show that full self-consistency in the GW scheme either systematically outperforms vertex-corrected G0W0 or gives results of at least comparative quality. Moreover, G0W0Γ results in additional computational cost when compared to G0W0 or self-consistent GW. The dependency of G0W0Γ on the starting mean-field solution is frequently more dominant than the magnitude of the vertex correction itself. Consequently, for molecular systems, self-consistent GW performed on the imaginary axis (and then followed by modern analytical continuation techniques) offers a more reliable approach to make predictions of ionization potentials.

Total Energy beyond GW: Exact Results and Guidelines for Approximations

The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.

The three channels of many-body perturbation theory: GW, particle-particle, and electron-hole T-matrix self-energies

We derive the explicit expression of the three self-energies that one encounters in many-body perturbation theory: the well-known GW self-energy, as well as the particle-particle and electron-hole T-matrix self-energies. Each of these can be easily computed via the eigenvalues and eigenvectors of a different random-phase approximation linear eigenvalue problem that completely defines their corresponding response function. For illustrative and comparative purposes, we report the principal ionization potentials of a set of small molecules computed at each level of theory. The performance of these schemes on strongly correlated systems (B2 and C2) is also discussed.

A Similarity Renormalization Group Approach to Green's Function Methods

The family of Green's function methods based on the GW approximation has gained popularity in the electronic structure theory thanks to its accuracy in weakly correlated systems combined with its cost-effectiveness. Despite this, self-consistent versions still pose challenges in terms of convergence. A recent study [Monino and Loos J. Chem. Phys. 2022, 156, 231101.] has linked these convergence issues to the intruder-state problem. In this work, a perturbative analysis of the similarity renormalization group (SRG) approach is performed on Green's function methods. The SRG formalism enables us to derive, from first-principles, the expression of a naturally static and Hermitian form of the self-energy that can be employed in quasiparticle self-consistent GW (qsGW) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qsGW calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code.

Embedding vertex corrections in GW self-energy: Theory, implementation, and outlook

The vertex function (Γ) within the Green's function formalism encapsulates information about all higher-order electron-electron interaction beyond those mediated by density fluctuations. Herein, we present an efficient approach that embeds vertex corrections in the one-shot GW correlation self-energy for isolated and periodic systems. The vertex-corrected self-energy is constructed through the proposed separation-propagation-recombination procedure: the electronic Hilbert space is separated into an active space and its orthogonal complement denoted as the "rest;" the active component is propagated by a space-specific effective Hamiltonian different from the rest. The vertex corrections are introduced by a rescaled time-dependent nonlocal exchange interaction. The direct Γ correction to the self-energy is further updated by adjusting the rescaling factor in a self-consistent post-processing cycle. Our embedding method is tested mainly on donor-acceptor charge-transfer systems. The embedded vertex effects consistently and significantly correct the quasiparticle energies of the gap-edge states. The fundamental gap is generally improved by 1-3 eV upon the one-shot GW approximation. Furthermore, we provide an outlook for applications of (embedded) vertex corrections in calculations of extended solids.

Natural virtual orbitals for the GW method in the random-phase approximation and beyond

The increasingly popular GW method is becoming a convenient tool to determine vertical ionization energies in molecular systems. However, depending on the formalism used and the range of orbitals investigated, it may be hampered by a steep computational scaling. To alleviate this issue, correlated natural virtual orbitals (NVOs) based on second-order Møller-Plesset (MP2) and direct MP2 correlation energies are implemented, and the resulting correlated NVOs are tested on GW quasiparticle energies. Test cases include the popular GW variants G₀W₀ and evGW₀ as well as more elaborate vertex corrections. We find that for increasingly larger molecular systems and basis sets, NVOs considerably improve efficiency. Furthermore, we test the performance of the truncated (frozen) NVO ansatz on the GW100 test set. For the latter, it is demonstrated that, using a carefully chosen truncation threshold, NVOs lead to a negligible loss in accuracy while providing speedups of one order of magnitude. Furthermore, we compare the resulting quasiparticle energies to very accurate vertical ionization energies obtained from coupled-cluster theory with singles, doubles, and noniterative triples [CCSD(T)], confirming that the loss in accuracy introduced by truncating the NVOs is negligible compared to the methodical errors in the GW approximation. It is also demonstrated that the choice of basis set impacts results far more than using a suitably truncated NVO space. Therefore, at the same computational expense, more accurate results can be obtained using NVOs. Finally, we provide improved reference CCSD(T) values for the GW100 test set, which have been obtained using the def2-QZVPP basis set.

Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies

With increasing interelectronic distance, the screening of the electron–electron interaction by the presence of other electrons becomes the dominant source of electron correlation. This effect is described by the random phase approximation (RPA) which is therefore a promising method for the calculation of weak interactions. The success of the RPA relies on the cancellation of errors, which can be traced back to the violation of the crossing symmetry of the 4-point vertex, leading to strongly overestimated total correlation energies. By the addition of second-order screened exchange (SOSEX) to the correlation energy, this issue is substantially reduced. In the adiabatic connection (AC) SOSEX formalism, one of the two electron–electron interaction lines in the second-order exchange term is dynamically screened (SOSEX(W, v_c)). A related SOSEX expression in which both electron–electron interaction lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, v_c), the evaluation of this correlation energy expression does not require an expensive numerical frequency integration and is therefore advantageous from a computational perspective. We compare the accuracy of the statically screened variant to RPA and RPA+SOSEX(W, v_c) for a wide range of chemical reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, v_c) for charged excitations and noncovalent interactions where they lead to major improvements over RPA.

Numerically Precise Benchmark of Many-Body Self-Energies on Spherical Atoms

We investigate the performance of beyond-GW approaches in many-body perturbation theory by addressing atoms described within the spherical approximation via a dedicated numerical treatment based on B-splines and spherical harmonics. We consider the GW, second Born (2B), and GW + second order screened exchange (GW+SOSEX) self-energies and use them to obtain ionization potentials from the quasi-particle equation (QPE) solved perturbatively on top of independent-particle calculations. We also solve the linearized Sham–Schlüter equation (LSSE) and compare the resulting xc potentials against exact data. We find that the LSSE provides consistent starting points for the QPE but does not present any practical advantage in the present context. Still, the features of the xc potentials obtained with it shed light on possible strategies for the inclusion of beyond-GW diagrams in the many-body self-energy. Our findings show that solving the QPE with the GW+SOSEX self-energy on top of a PBE or PBE0 solution is a viable scheme to go beyond GW in finite systems, even in the atomic limit. However, GW shows a comparable performance if one agrees to use a hybrid starting point. We also obtain promising results with the 2B self-energy on top of Hartree–Fock, suggesting that the full time-dependent Hartree–Fock vertex may be another viable beyond-GW scheme for finite systems.

Simplified GW/BSE Approach for Charged and Neutral Excitation Energies of Large Molecules and Nanomaterials

Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [Grimme, S. J. Chem. Phys. 2013, 138, 244104], we describe a simplified approach to excited-state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from ab initio GW calculations by only 0.2 eV. The performance of sBSE (based on the sGW input) is tested on the excitation energies of molecules in Thiel's set, for which it differs from ab initio GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms.

Assessing the G0W0Γ0(1) Approach: Beyond G0W0 with Hedin's Full Second-Order Self-Energy Contribution

We present and benchmark a self-energy approach for quasiparticle energy calculations that goes beyond Hedin's GW approximation by adding the full second-order self-energy (FSOS-W) contribution. The FSOS-W diagram involves two screened Coulomb interaction (W) lines, and adding the FSOS-W to the GW self-energy can be interpreted as first-order vertex correction to GW (GWΓ⁽¹⁾). Our FSOS-W implementation is based on the resolution-of-identity technique and exhibits better than O(N⁵) scaling with system size for small- to medium-sized molecules. We then present one-shot GWΓ⁽¹⁾ (G₀W₀Γ₀⁽¹⁾) benchmarks for the GW100 test set and a set of 24 acceptor molecules. For semilocal or hybrid density functional theory starting points, G₀W₀Γ₀⁽¹⁾ systematically outperforms G₀W₀ for the first vertical ionization potentials and electron affinities of both test sets. Finally, we demonstrate that a static FSOS-W self-energy significantly underestimates the quasiparticle energies.

Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation

The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have, therefore, been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham bandgaps and dielectric constants, and comment on the issue of self-consistency.

Cubic-Scaling All-Electron GW Calculations with a Separable Density-Fitting Space-Time Approach

We present an implementation of the GW space-time approach that allows cubic-scaling all-electron calculations with standard Gaussian basis sets without exploiting any localization or sparsity considerations. The independent-electron susceptibility is constructed in a time representation over a nonuniform distribution of real-space locations {r_k} optimized within a separable resolution-of-the-identity framework to reproduce standard Coulomb-fitting calculations with meV accuracy. The compactness of the obtained {r_k} distribution leads to a crossover with the standard Coulomb-fitting scheme for system sizes below a few hundred electrons. The needed analytic continuation follows a recent approach that requires the continuation of the screened Coulomb potential rather than the much more structured self-energy. The present scheme is benchmarked over large molecular sets, and scaling properties are demonstrated on a family of defected hexagonal boron-nitride flakes containing up to 6000 electrons.

Electronic and optical properties of fluorinated graphene within many-body Green's function framework

In this article, a systematic examination of the electronic and optical properties of partially fluorinated graphene is presented. In order to capture a large variety of fluorination degrees and configurations, different sizes of the supercell combining with various degrees of fluorination are considered. On top of periodic density functional theory, the G₀W₀ method and the G₀W₀Γ method within many-body Green's function framework are employed. Including the description of electron-hole interactions, the optical spectra based on the Bethe-Salpeter equation are calculated. Two-sided fluorination with compact fluorination arrangements is energetically most favorable. The fluorination degree has a determined impact on the bandgap value in the system, while the fluorination pattern strongly influences the characteristics of the bands in the electronic structures. Depending on the polarization of the applied electromagnetic field, the optical absorption spectra of the same structure could vary significantly. These interesting results suggest the potential applications of partially fluorinated graphene as optoelectronic materials.

Advanced First-Principle Modeling of Relativistic Ruddlesden—Popper Strontium Iridates

In this review, we provide a survey of the application of advanced first-principle methods on the theoretical modeling and understanding of novel electronic, optical, and magnetic properties of the spin-orbit coupled Ruddlesden–Popper series of iridates Srn+1IrnO3n+1 (n = 1, 2, and ∞). After a brief description of the basic aspects of the adopted methods (noncollinear local spin density approximation plus an on-site Coulomb interaction (LSDA+U), constrained random phase approximation (cRPA), GW, and Bethe–Salpeter equation (BSE)), we present and discuss select results. We show that a detailed phase diagrams of the metal–insulator transition and magnetic phase transition can be constructed by inspecting the evolution of electronic and magnetic properties as a function of Hubbard U, spin–orbit coupling (SOC) strength, and dimensionality n, which provide clear evidence for the crucial role played by SOC and U in establishing a relativistic (Dirac) Mott–Hubbard insulating state in Sr2IrO4 and Sr3Ir2O7. To characterize the ground-state phases, we quantify the most relevant energy scales fully ab initio—crystal field energy, Hubbard U, and SOC constant of three compounds—and discuss the quasiparticle band structures in detail by comparing GW and LSDA+U data. We examine the different magnetic ground states of structurally similar n = 1 and n = 2 compounds and clarify that the origin of the in-plane canted antiferromagnetic (AFM) state of Sr2IrO4 arises from competition between isotropic exchange and Dzyaloshinskii–Moriya (DM) interactions whereas the collinear AFM state of Sr3Ir2O7 is due to strong interlayer magnetic coupling. Finally, we report the dimensionality controlled metal–insulator transition across the series by computing their optical transitions and conductivity spectra at the GW+BSE level from the the quasi two-dimensional insulating n = 1 and 2 phases to the three-dimensional metallic n=∞ phase.

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

We derive a low-scaling G0W0 algorithm for molecules using pair atomic density fitting (PADF) and an imaginary time representation of the Green's function and describe its implementation in the Slater type orbital (STO)-based Amsterdam density functional (ADF) electronic structure code. We demonstrate the scalability of our algorithm on a series of water clusters with up to 432 atoms and 7776 basis functions and observe asymptotic quadratic scaling with realistic threshold qualities controlling distance effects and basis sets of triple-ζ (TZ) plus double polarization quality. Also owing to a very small prefactor, a G0W0 calculation for the largest of these clusters takes only 240 CPU hours with these settings. We assess the accuracy of our algorithm for HOMO and LUMO energies in the GW100 database. With errors of 0.24 eV for HOMO energies on the quadruple-ζ level, our implementation is less accurate than canonical all-electron implementations using the larger def2-QZVP GTO-type basis set. Apart from basis set errors, this is related to the well-known shortcomings of the GW space-time method using analytical continuation techniques as well as to numerical issues of the PADF approach of accurately representing diffuse atomic orbital (AO) products. We speculate that these difficulties might be overcome by using optimized auxiliary fit sets with more diffuse functions of higher angular momenta. Despite these shortcomings, for subsets of medium and large molecules from the GW5000 database, the error of our approach using basis sets of TZ and augmented double-ζ (DZ) quality is decreasing with system size. On the augmented DZ level, we reproduce canonical, complete basis set limit extrapolated reference values with an accuracy of 80 meV on average for a set of 20 large organic molecules. We anticipate our algorithm, in its current form, to be very useful in the study of single-particle properties of large organic systems such as chromophores and acceptor molecules.

The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as measured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW, a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physicists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

GW100: Comparison of Methods and Accuracy of Results Obtained with the WEST Code

The reproducibility of calculations carried out within many-body perturbation theory at the G₀ W₀ level is assessed for 100 closed shell molecules and compared to that of density functional theory. We consider vertical ionization potentials (VIP) and electron affinities (VEA) obtained with five different codes: BerkeleyGW, FHI-aims, TURBOMOLE, VASP, and WEST. We review the approximations and parameters that control the accuracy of G₀ W₀ results in each code, and we discuss in detail the effect of extrapolation techniques for the parameters entering the WEST code. Differences between the VIP and VEA computed with the various codes are within ∼60 and ∼120 meV, respectively, which is up to four times larger than in the case of the best results obtained with DFT codes. Vertical ionization potentials are validated against experiment and CCSD(T) quantum chemistry results showing a mean absolute relative error of ∼4% for data obtained with WEST. Our analysis of the differences between localized orbitals and plane-wave implementations points out molecules containing Cu, I, Ga, and Xe as major sources of discrepancies, which call for a re-evaluation of the pseudopotentials used for these systems in G₀ W₀ calculations.

Assessing GW Approaches for Predicting Core Level Binding Energies

Here we present a systematic study on the performance of different GW approaches: G₀W₀, G₀W₀ with linearized quasiparticle equation (lin-G₀W₀), and quasiparticle self-consistent GW (qsGW), in predicting core level binding energies (CLBEs) on a series of representative molecules comparing to Kohn-Sham (KS) orbital energy-based results. KS orbital energies obtained using the PBE functional are 20-30 eV lower in energy than experimental values obtained from X-ray photoemission spectroscopy (XPS), showing that any Koopmans-like interpretation of KS core level orbitals fails dramatically. Results from qsGW lead to CLBEs that are closer to experimental values from XPS, yet too large. For the qsGW method, the mean absolute error is about 2 eV, an order of magnitude better than plain KS PBE orbital energies and quite close to predictions from ΔSCF calculations with the same functional, which are accurate within ∼1 eV. Smaller errors of ∼0.6 eV are found for qsGW CLBE shifts, again similar to those obtained using ΔSCF PBE. The computationally more affordable G₀W₀ approximation leads to results less accurate than qsGW, with an error of ∼9 eV for CLBEs and ∼0.9 eV for their shifts. Interestingly, starting G₀W₀ from PBE0 reduces this error to ∼4 eV with a slight improvement on the shifts as well (∼0.4 eV). The validity of the G₀W₀ results is however questionable since only linearized quasiparticle equation results can be obtained. The present results pave the way to estimate CLBEs in periodic systems where ΔSCF calculations are not straightforward although further improvement is clearly needed.

The Bethe–Salpeter equation in chemistry: relations with TD-DFT, applications and challenges

We review the Bethe–Salpeter formalism and analyze its performances for the calculation of the excited state properties of molecular systems.