Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods

Molecular Ionization Energies from GW and Hartree-Fock Theory: Polarizability, Screening, and Self-Energy Vertex Corrections

Accurate prediction of electron removal and addition energies is essential for reproducing neutral excitation spectra in molecules using Bethe-Salpeter equation methods. A Hartree-Fock starting point for GW/BSE calculations, combined with a random phase approximation (RPA) polarizability in the screened interaction, W, is well-known to overestimate neutral excitation energies. Using a Hartree-Fock starting point, we apply several different approximations for W to molecules in the Quest-3 database [Loos et al. J. Chem. Theory Comput. 2020, 16, 1711]. W is calculated using polarizabilities in RPA and time-dependent HF approximations. Inclusion of screened electron-hole attraction in the polarizability yields valence ionization energies in better agreement with experimental values and ADC(3) calculations than the more commonly applied RPA polarizability. Quasiparticle weights are also in better agreement with ADC(3) values when electron-hole attraction is included in W. Shake-up excitations for the 1π levels in benzene and azines are indicated only when electron-hole attraction is included. Ionization energies derived from HF eigenvalues via Koopmans theorem for molecules with nitrogen or oxygen lone pairs have the largest differences from experimental values in the molecules considered, leading to incorrect ordering of nonbonding and π bonding levels. Inclusion of electron-hole attraction in the polarizability results in correct ordering of ionization energies and marked improvement in agreement with experimental data. Vertex corrections to the self-energy further improve agreement with experimental ionization energies for these localized states.

Reference CC3 Excitation Energies for Organic Chromophores: Benchmarking TD-DFT, BSE/GW, and Wave Function Methods

To expand the QUEST database of highly accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly accurate excitation energies. For several singlet transitions, we have been able to determine CCSDT transition energies with a compact basis set, finding minimal deviations from the CC3 values for most states. Subsequently, we employ these reference values to benchmark a series of lower-order wave function approaches, including the popular ADC(2) and CC2 schemes, as well as time-dependent density-functional theory (TD-DFT), both with and without applying the Tamm-Dancoff approximation (TDA). At the TD-DFT level, we evaluate a large panel of global, range-separated, local, and double hybrid functionals. Additionally, we assess the performance of the Bethe-Salpeter equation (BSE) formalism relying on both G0W0 and evGW quasiparticle energies evaluated from various starting points. It turns out that CC2 and ADC(2.5) are the most accurate models among those with respective O ( N 5 ) and O ( N 6 ) scalings with system size. In contrast, CCSD does not outperform CC2. The best performing exchange-correlation functionals include BMK, M06-2X, M06-SX, CAM-B3LYP, ωB97X-D, and LH20t, with average deviations of approximately 0.20 eV or slightly below. Errors on vertical excitation energies can be further reduced by considering double hybrids. Both SOS-ωB88PP86 and SOS-ωPBEPP86 exhibit particularly attractive performances with overall quality on par with CC2, whereas PBE0-DH and PBE-QIDH are only slightly less efficient. BSE/evGW calculations based on Kohn-Sham starting points have been found to be particularly effective for singlet transitions, but much less for their triplet counterparts.

Polarizable Continuum Models and Green's Function GW Formalism: On the Dynamics of the Solvent Electrons

The many-body GW formalism, for the calculation of ionization potentials or electronic affinities, relies on the frequency-dependent dielectric function built from the electronic degrees of freedom. Considering the case of water as a solvent treated within the polarizable continuum model, we explore the impact of restricting the full frequency-dependence of the solvent electronic dielectric response to a frequency-independent (ϵ∞) optical dielectric constant. For solutes presenting small to large highest-occupied to lowest-unoccupied molecular orbital energy gaps, we show that such a restriction induces errors no larger than a few percent on the energy level shifts from the gas to the solvated phase. We further introduce a remarkably accurate single-pole model for mimicking the effect of the full frequency dependence of the water dielectric function in the visible-UV range. This allows a fully dynamical embedded GW calculation with the only knowledge of the cavity reaction field calculated for the ϵ∞ optical dielectric constant.

BSE@GW Prediction of Charge Transfer Exciton in Molecular Complexes: Assessment of Self-Energy and Exchange-Correlation Dependence

The Bethe-Salpeter equation using the GW approximation to the self-energy (BSE@GW) is a computationally attractive method for studying electronic excitation from first principles within the many-body Green's function theory framework. We examine its dependence on the underlying exchange-correlation (XC) approximation as well as on the GW approximation for predicting the charge transfer exciton formation at representative type-II interfaces between molecular systems of tetrachloro-1,2-benzoquinone (TCBQ) and acene derivatives. For the XC approximation, we consider several widely used generalized gradient approximation (GGA) and hybrid GGA functionals. For the GW self-energy approximation, we examine the recently proposed renormalized singles approach by Yang and coauthors [J. Phys. Chem. Lett. 2019, 10 (3), 447-452; J. Chem. Theory Comput. 2022, 18, 7570-7585] in addition to other commonly employed approximated GW schemes. We demonstrate a reliable prediction of the charge transfer exciton within the BSE@GW level of theory.

Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods

Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.

Ab Initio Electron Propagators with an Hermitian, Intermediately Normalized Superoperator Metric Applied to Vertical Electron Affinities

New-generation ab initio electron propagator methods for calculating electron detachment energies of closed-shell molecules and anions have surpassed their predecessors' accuracy and computational efficiency. Derived from an Hermitian, intermediately normalized superoperator metric, these methods contain no adjustable parameters. To assess their versatility, a standard set (NIST-50-EA) of 50 vertical electron affinities of small closed-shell molecules based on NIST reference data has been created. Errors with respect to reference data on 23 large, conjugated organic photovoltaic (OPV23) molecules have also been analyzed. All final states are valence anions that correspond to electron affinities between 0.2 and 4.2 eV. For a given scaling of the arithmetic bottleneck, the new-generation methods obtain the lowest mean absolute errors (MAEs). The best methods with fifth-power arithmetic scaling realize MAEs below 0.1 eV. Composite models comprising cubically and quintically scaling calculations executed with large and small basis sets, respectively, produce OPV23 MAEs near 0.05 eV. The accuracy of quintically scaling methods executed with large basis sets is thereby procured with reduced computational effort. New-generation results obtained with and without the diagonal self-energy approximation in the canonical Hartree-Fock basis have been compared. These results indicate that Dyson orbitals closely resemble canonical Hartree-Fock orbitals multiplied by the square root of a probability factor above 0.85.

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.

The C1s core levels of polycyclic aromatic hydrocarbons and styrenic polymers: A first-principles study

Understanding core level shifts in aromatic compounds is crucial for the correct interpretation of x-ray photoelectron spectroscopy (XPS) of polycyclic aromatic hydrocarbons (PAHs), including acenes, as well as of styrenic polymers, which are increasingly relevant for the microelectronic industry, among other applications. The effect of delocalization through π aromatic systems on the stabilization of valence molecular orbitals has been widely investigated in the past. However, little has been reported on the impact on the deeper C1s core energy levels. In this work, we use first-principles calculations at the level of many body perturbation theory to compute the C1s binding energies of several aromatic systems. We report a C1s red shift in PAHs and acenes of increasing size, both in the gas phase and in the molecular crystal. C1s red shifts are also calculated for stacked benzene and naphthalene pairs at decreasing intermolecular distances. A C1s red shift is in addition found between oligomers of poly(p-hydroxystyrene) and polystyrene of increasing length, which we attribute to ring-ring interactions between the side-chains. The predicted shifts are larger than common instrumental errors and could, therefore, be detected in XPS experiments.

Benchmarking Ionization Potentials from pCCD Tailored Coupled Cluster Models

The ionization potential (IP) is an important parameter providing essential insights into the reactivity of chemical systems. IPs are also crucial for designing, optimizing, and understanding the functionality of modern technological devices. We recently showed that limiting the CC ansatz to the seniority-zero sector proves insufficient in predicting reliable and accurate ionization potentials within an IP equation-of-motion coupled-cluster formalism. Specifically, the absence of dynamical correlation in the seniority-zero pair coupled cluster doubles (pCCD) model led to unacceptably significant errors of approximately 1.5 eV. In this work, we aim to explore the impact of dynamical correlation and the choice of the molecular orbital basis (canonical vs localized) in CC-type methods targeting 230 ionized states in 70 molecules, comprising small organic molecules, medium-sized organic acceptors, and nucleobases. We focus on pCCD-based approaches as well as the conventional IP-EOM-CCD and IP-EOM-CCSD. Their performance is compared to the CCSD(T) or CCSDT equivalent and experimental reference data. Our statistical analysis reveals that all investigated frozen-pair coupled cluster methods exhibit similar performance, with differences in errors typically within chemical accuracy (1 kcal/mol or 0.05 eV). Notably, the effect of the molecular orbital basis, such as canonical Hartree-Fock or natural pCCD-optimized orbitals, on the IPs is marginal if dynamical correlation is accounted for. Our study suggests that triple excitations are crucial in achieving chemical accuracy in IPs when modeling electron detachment processes with pCCD-based methods.

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe-Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

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.

Static versus dynamically polarizable environments within the many-body GW formalism

Continuum- or discrete-polarizable models for the study of optoelectronic processes in embedded subsystems rely mostly on the restriction of the surrounding electronic dielectric response to its low frequency limit. Such a description hinges on the assumption that the electrons in the surrounding medium react instantaneously to any excitation in the central subsystem, thus treating the environment in the adiabatic limit. Exploiting a recently developed embedded GW formalism with an environment described at the fully ab initio level, we assess the merits of the adiabatic limit with respect to an environment where the full dynamics of the dielectric response are considered. Furthermore, we show how to properly take the static limit of the environment's susceptibility by introducing the so-called Coulomb-hole and screened-exchange contributions to the reaction field. As a first application, we consider a C60 molecule at the surface of a C60 crystal, namely, a case where the dynamics of the embedded and embedding subsystems are similar. The common adiabatic assumption, when properly treated, generates errors below 10% on the polarization energy associated with frontier energy levels and associated energy gaps. Finally, we consider a water molecule inside a metallic nanotube, the worst case for the environment's adiabatic limit. The error on the gap polarization energy remains below 10%, even though the error on the frontier orbital polarization energies can reach a few tenths of an electronvolt.

Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels

We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.

New-Generation Electron-Propagator Methods for Molecular Electron-Binding Energies

A new generation of electron-propagator methods for the calculation of electron binding energies has surpassed its antecedents with respect to accuracy, efficiency, and interpretability. No adjustable parameters are introduced in these fully ab initio procedures. Numerical tests versus several databases of valence, vertical electron binding energies of closed-shell molecules and atoms have been performed. Easily interpreted self-energy approximations with cubic arithmetic scaling produce mean absolute errors (MAEs) of 0.2 and 0.3 eV for electron detachments and attachments, respectively. The most accurate explicitly renormalized methods with fifth-power arithmetic scaling yield MAEs below 0.1 eV for detachments and attachments. Approximate renormalization leads to more efficient fifth-power alternatives for electron detachments that achieve similar accuracy with fewer bottleneck operations. Composite protocols generate excellent predictions versus highly accurate basis-extrapolated standards and experiments. The validity of the diagonal self-energy approximation and the accuracy of the approximate renormalizations are confirmed. The success of these perturbative methods based on canonical Hartree-Fock orbitals rests on a Hermitized, intermediately normalized superoperator metric. The results of all of the new-generation calculations may be analyzed in terms of final-state orbital relaxation and differential correlation effects.

Reaching High Accuracy for Energetic Properties at Second-Order Perturbation Cost by Merging Self-Consistency and Spin-Opposite Scaling

Quantum chemical methods dealing with challenging systems while retaining low computational costs have attracted attention. In particular, many efforts have been devoted to developing new methods based on second-order perturbation that may be the simplest correlated method beyond Hartree-Fock. We have recently developed a self-consistent perturbation theory named one-body Møller-Plesset second-order perturbation theory (OBMP2) and shown that it can resolve issues caused by the noniterative nature of standard perturbation theory. In this work, we extend the method by introducing spin-opposite scaling to the double-excitation amplitudes, resulting in the O2BMP2 method. We assess the O2BMP2 performance on the triple-bond N2 dissociation, singlet-triplet gaps, and ionization potentials. O2BMP2 performs much better than standard MP2 and reaches the accuracy of coupled-cluster methods in all cases considered in this work.

Peak Broadening in Photoelectron Spectroscopy of Amorphous Polymers: The Leading Role of the Electrostatic Landscape

The broadening in photoelectron spectra of polymers can be attributed to several factors, such as light source spread, spectrometer resolution, the finite lifetime of the hole state, and solid-state effects. Here, for the first time, we set up a computational protocol to assess the peak broadening induced for both core and valence levels by solid-state effects in four amorphous polymers by using a combination of density functional theory, many-body perturbation theory, and classical polarizable embedding. We show that intrinsic local inhomogeneities in the electrostatic environment induce a Gaussian broadening of 0.2-0.7 eV in the binding energies of both core and semivalence electrons, corresponding to a full width at half-maximum (FWHM) of 0.5-1.7 eV for the investigated systems. The induced broadening is larger in acrylate-based than in styrene-based polymers, revealing the crucial role of polar groups in controlling the roughness of the electrostatic landscape in the solid matrix.

New-Generation Electron-Propagator Methods for Calculations of Electron Affinities and Ionization Energies: Tests on Organic Photovoltaic Molecules

A new generation of ab initio electron-propagator self-energies recently superseded its antecedents' accuracy and computational efficiency in calculating vertical ionization energies (VIEs) of closed-shell molecules. (See J. Chem. Phys. 2021, 155, 204107, J. Chem. Theory Comput. 2022, 18, 4927, J. Chem. Phys. 2023, 159, 124109.) No adjustable parameters were introduced in the generation of reference orbitals or in the construction of self-energies. The same approach has been extended in this work to vertical electron affinities (VEAs). Calculations were performed on 24 conjugated, organic photovoltaic molecules with diverse functional groups. These molecules are considerably larger than those studied in previous tests on VIEs. Several new-generation self-energies produce mean absolute errors (MAEs) below 0.1 eV versus ΔCCSD(T) (i.e., total energy differences from the coupled-cluster singles, doubles, and perturbative triples method) VIEs and VEAs obtained with identical basis sets. A composite model employs cubically and quintically scaling algorithms and power-law basis-set extrapolations based on augmented double-triple or triple-quadruple ζ data. Its MAEs are near 0.05 eV versus benchmark values, with 0.03 eV error bars for the lowest VIE and the highest VEA of each molecule. A more efficient and equally accurate composite model for calculating VIEs avoids full transformations of electron repulsion integrals to the molecular orbital basis. High probability factors support the diagonal self-energy approximation, wherein Dyson orbitals are proportional to canonical, Hartree-Fock orbitals.

A new generation of non-diagonal, renormalized self-energies for calculation of electron removal energies

A new generation of diagonal self-energies for the calculation of electron removal energies of molecules and molecular ions that has superseded its predecessors with respect to accuracy, efficiency, and interpretability is extended to include non-diagonal self-energies that permit Dyson orbitals to be expressed as linear combinations of canonical Hartree-Fock orbitals. In addition, an improved algorithm for renormalized methods eliminates the convergence difficulties encountered in the first studies of the new, diagonal self-energies. A dataset of outer-valence, vertical ionization energies with almost full-configuration-interaction quality serves as a standard of comparison in numerical tests. The new non-diagonal, renormalized methods are slightly more accurate than their diagonal counterparts, with mean absolute errors between 0.10 and 0.06 eV for outer-valence final states. This advantage is procured at the cost of an increase in the scaling of arithmetic bottlenecks that accompany the inclusion of non-diagonal self-energy terms. The new, non-diagonal, renormalized self-energies are also more accurate and efficient than their non-diagonal predecessors.

Two-Component GW Calculations: Cubic Scaling Implementation and Comparison of Vertex-Corrected and Partially Self-Consistent GW Variants

We report an all-electron, atomic orbital (AO)-based, two-component (2C) implementation of the GW approximation (GWA) for closed-shell molecules. Our algorithm is based on the space-time formulation of the GWA and uses analytical continuation (AC) of the self-energy, and pair-atomic density fitting (PADF) to switch between AO and auxiliary basis. By calculating the dynamical contribution to the GW self-energy at a quasi-one-component level, our 2C-GW algorithm is only about a factor of 2-3 slower than in the scalar relativistic case. Additionally, we present a 2C implementation of the simplest vertex correction to the self-energy, the statically screened G3W2 correction. Comparison of first ionization potentials (IPs) of a set of 67 molecules with heavy elements (a subset of the SOC81 set) calculated with our implementation against results from the WEST code reveals mean absolute deviations (MAD) of around 70 meV for G0W0@PBE and G0W0@PBE0. We check the accuracy of our AC treatment by comparison to full-frequency GW calculations, which shows that in the absence of multisolution cases, the errors due to AC are only minor. This implies that the main sources of the observed deviations between both implementations are the different single-particle bases and the pseudopotential approximation in the WEST code. Finally, we assess the performance of some (partially self-consistent) variants of the GWA for the calculation of first IPs by comparison to vertical experimental reference values. G0W0@PBE0 (25% exact exchange) and G0W0@BHLYP (50% exact exchange) perform best with mean absolute deviations (MAD) of about 200 meV. Explicit treatment of spin-orbit effects at the 2C level is crucial for systematic agreement with experiment. On the other hand, eigenvalue-only self-consistent GW (evGW) and quasi-particle self-consistent GW (qsGW) significantly overestimate the IPs. Perturbative G3W2 corrections increase the IPs and therefore improve the agreement with experiment in cases where G0W0 alone underestimates the IPs. With a MAD of only 140 meV, 2C-G0W0@PBE0 + G3W2 is in best agreement with the experimental reference values.

Accurate GW frontier orbital energies of 134 kilo molecules

HOMO and LUMO energies are critical molecular properties that typically require high accuracy computations for practical applicability. Until now, a comprehensive dataset containing sufficiently accurate HOMO and LUMO energies has been unavailable. In this study, we introduce a new dataset of HOMO/LUMO energies for QM9 compounds, calculated using the GW method. The GW method offers adequate HOMO/LUMO prediction accuracy for diverse applications, exhibiting mean unsigned errors of 100 meV in the GW100 benchmark dataset. This database may serve as a benchmark of HOMO/LUMO prediction, delta-learning, and transfer learning, particularly for larger molecules where GW is the most accurate but still numerically feasible method. We anticipate that this dataset will enable the development of more accurate machine learning models for predicting molecular properties.

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.

Assessing the Role of the Kohn-Sham Density in the Calculation of the Low-Lying Bethe-Salpeter Excitation Energies

We adopt the GW many-body perturbation theory in conjunction with the Bethe-Salpeter equation (BSE) to compute 57 excitation energies of a set of 37 molecules. By using the PBEh global hybrid functional and a self-consistent scheme on the eigenvalues in GW, we show a strong dependence of the BSE energy on the starting Kohn-Sham (KS) density functional. This arises from both the quasiparticle energies and the spatial localization of the frozen KS orbitals employed to compute the BSE. In order to address the arbitrariness in the mean field choice, we adopt an orbital-tuning scheme where the amount of Fock exchange, α, is tuned to impose that the KS HOMO matches the GW quasiparticle eigenvalue, thus fulfilling the ionization potential theorem in DFT. The performance of the proposed scheme yields excellent results and it is similar to M06-2X and PBEh with α = 75%, consistent with tuned values of α ranging between 60% and 80%.

True molecular conformation and structure determination by three-dimensional electron diffraction of PAH by-products potentially useful for electronic applications

The true molecular conformation and the crystal structure of benzo[e]dinaphtho[2,3-a;1',2',3',4'-ghi]fluoranthene, 7,14-diphenylnaphtho[1,2,3,4-cde]bisanthene and 7,16-diphenylnaphtho[1,2,3,4-cde]helianthrene were determined ab initio by 3D electron diffraction. All three molecules are remarkable polycyclic aromatic hydrocarbons. The molecular conformation of two of these compounds could not be determined via classical spectroscopic methods due to the large size of the molecule and the occurrence of multiple and reciprocally connected aromatic rings. The molecular structure of the third molecule was previously considered provisional. These compounds were isolated as by-products in the synthesis of similar products and were at the same time nanocrystalline and available only in very limited amounts. 3D electron diffraction data, taken from submicrometric single crystals, allowed for direct ab initio structure solution and the unbiased determination of the internal molecular conformation. Detailed synthetic routes and spectroscopic analyses are also discussed. Based on many-body perturbation theory simulations, benzo[e]dinaphtho[2,3-a;1',2',3',4'-ghi]fluoranthene may be a promising candidate for triplet-triplet annihilation and 7,14-diphenylnaphtho[1,2,3,4-cde]bisanthene may be a promising candidate for intermolecular singlet fission in the solid state.

Accurate Prediction of Vertical Ionization Potentials and Electron Affinities from Spin-Component Scaled CC2 and ADC(2) Models

The CC2 and ADC(2) wave function models and their spin-component scaled modifications are adopted for predicting vertical ionization potentials (VIPs) and electron affinities (VEAs). The ionic solutions are obtained as electronic excitations in the continuum orbital formalism, making possible the use of existing, widespread quantum chemistry codes with minimal modifications, in full consistency with the treatment of charge transfer excitations. The performance of different variants is evaluated via benchmark calculations on various sets from previous works, containing small- and medium-sized systems, including the nucleobases. It is shown that with the spin-scaled approximate methods, in particular the scaled opposite-spin variant of the ADC(2) method, the accuracy of EOM-CCSD is achievable at a fraction of the computational cost, also outperforming many common electron propagator approaches.

Renormalized Singles with Correlation in GW Green's Function Theory for Accurate Quasiparticle Energies

We apply the renormalized singles with the correlation (RSc) Green function in the GW approximation for accurate quasiparticle (QP) energies and orbitals. The RSc Green function includes singles contributions from the associated density functional approximation (DFA) and considers correlation contributions perturbatively. G_RScW_RSc uses the RSc Green function as the new starting point and in the formulation of the screened interaction. G_RScW₀ fixes the screened interaction at the DFA level. For the calculations of ionization potentials, G_RScW_RSc and G_RScW₀ significantly reduce the starting point dependence and provide accurate results with errors around 0.2  eV. For the calculations of core-level binding energies, G_RScW_RSc slightly overestimates the results because of underscreening, but G_RScW₀ with GGA functionals provides the optimal accuracy with errors of 0.40  eV. We also show that G_RScW_RSc predicts accurate dipole moments. G_RScW_RSc and G_RScW₀, are computationally favorable compared with any self-consistent GW methods. The RSc approach is promising for making GW and other Green function methods efficient and robust.

Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems

The GW-Bethe–Salpeter equation (BSE) method is promising for calculating the low-lying excitonic states of molecular systems. However, so far it has only been applied to rather small molecules and in the commonly implemented diagonal approximations to the electronic self-energy, it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency can lead to an unfaithful description of some excitonic states for chlorophyll dimers while the qsGW-BSE vertical excitation energies (VEEs) are in excellent agreement with spectroscopic experiments for chlorophyll monomers and dimers measured in the gas phase. Furthermore, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use the new qsGW-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11,000 (6000) basis functions, the calculation could be completed in less than 5 (2) days on a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters in the low-energy spectrum of the hexameric complex. Earlier works with RSH kernels have demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environmental electrostatics.

Investigation of Ionization Potential in Quantum Dots Using the Stratified Stochastic Enumeration of Molecular Orbitals Method

The overarching goal of this work is to investigate the size-dependent characteristics of the ionization potential of PbS and CdS quantum dots. The ionization potentials of quantum dots provide critical information about the energies of occupied states, which can then be used to quantify the electron-removal characteristics of quantum dots. The energy of the highest-occupied molecular orbital is used to understand electron-transfer processes when invesigating the energy-level alignment between quantum dots and electron-accepting ligands. Ionization potential is also important for investigating and interpreting electron-detachment processes induced by light (photoelectron spectra), external voltage (chemiresistance), and collision with other electrons (impact ionization). Accurate first-principles calculations of ionization potential continue to be challenging because of the computational cost associated with the construction of the frequency-dependent self-energy operator and the numerical solution of the associated Dyson equation. The computational cost becomes prohibitive as the system size increases because of the large number of 2particle-1hole (2p1h) and 1particle-2hole (1p2h) terms needed for the calculation. In this work, we present the Stratified Stochastic Enumeration of Molecular Orbitals (SSE-MO) method for efficient construction of the self-energy operator. The SSE-MO method is a real-space method and the central strategy of this method is to use stochastically enumerated sampling of molecular orbitals and molecular-orbital indices for the construction of the 2p1h and 1p2h terms. This is achieved by first constructing a composite MO-index Cartesian coordinate space followed by transformation of the frequency-dependent self-energy operator to this composite space. The evaluation of both the real and imaginary components of the self-energy operator is performed using a stratified Monte Carlo technique. The SSE-MO method was used to calculate the ionization potentials and the frequency-dependent spectral functions for a series of PbS and CdS quantum dots by solving the Dyson equation using both single-shot and iterative procedures. The ionization potentials for both PbS and CdS quantum dots were found to decrease with increasing dot size. Analysis of the frequency-dependent spectral functions revealed that for PbS quantum dots the intermediate dot size exhibited a longer relative lifetime whereas in CdS the smallest dot size had the longest relative lifetime. The results from these calculations demonstrate the efficacy of the SSE-MO method for calculating accurate ionization potentials and spectral functions of chemical systems.

Modeling X-ray Photoelectron Spectroscopy of Macromolecules Using GW

We propose a simple additive approach to simulate X-ray photoelectron spectra (XPS) of macromolecules based on the GW method. Single-shot GW (G₀W₀) is a promising technique to compute accurate core-electron binding energies (BEs). However, its application to large molecules is still unfeasible. To circumvent the computational cost of G₀W₀, we break the macromolecule into tractable building blocks, such as isolated monomers, and sum up the theoretical spectra of each component, weighted by their molar ratio. In this work, we provide a first proof of concept by applying the method to four test polymers and one copolymer and show that it leads to an excellent agreement with experiments. The method could be used to retrieve the composition of unknown materials and study chemical reactions, by comparing the simulated spectra with experimental ones.

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

Unphysical discontinuities, intruder states and regularization in GW methods

By recasting the non-linear frequency-dependent GW quasiparticle equation into a linear eigenvalue problem, we explain the appearance of multiple solutions and unphysical discontinuities in various physical quantities computed within the GW approximation. Considering the GW self-energy as an effective Hamiltonian, it is shown that these issues are key signatures of strong correlation in the (N ± 1)-electron states and can be directly related to the intruder state problem. A simple and efficient regularization procedure inspired by the similarity renormalization group is proposed to avoid such issues and speed up the convergence of partially self-consistent GW calculations.

Recommendation of Orbitals for G0W0 Calculations on Molecules and Crystals

The many-body GW approximation, especially the G₀W₀ method, has been widely used for condensed matter and molecules to calculate quasiparticle energies for ionization, electron attachment, and band gaps. Because G₀W₀ calculations are well-known to have a strong dependence on the orbitals, the goal of the present work is to provide guidance on the choice of density functional used to generate orbitals and to recommend a choice that gives the most broadly accurate results. We have systematically investigated the dependence of G₀W₀ calculations on the orbitals for 100 molecules and 8 crystals by considering orbitals obtained with a diverse set of Kohn-Sham (KS) and generalized KS (GKS) functionals (63 functionals plus Hartree-Fock). The percentage of Hartree-Fock exchange employed in density functionals has been found to have strong influence on the predicted molecular ionization energy and crystal fundamental band gaps (with optimum values between 40 and 56%), but to have less effect on predicting molecular electron affinities. The low cost of the Gaussian implementation, even with hybrid functionals in periodic calculations, the better performance of global hybrids as compared to range-separated hybrids of either than screened exchange or long-range-corrected type, and the relatively low cost of global-hybrid-functional periodic calculations using Gaussians means that one can employ global-hybrid functionals at a very reasonable cost and obtain more accurate band gaps of semiconductors than are obtained by the methods currently widely employed, namely local gradient approximations. We single out three global-hybrid functionals that give especially good results for both molecules (100 in the test set) and crystals (8 in the test set, for all of which our benchmark data are the proper band gap rather than an optical band gap uncorrected for exciton effects).

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.

The GW Miracle in Many-Body Perturbation Theory for the Ionization Potential of Molecules

We use the GW100 benchmark set to systematically judge the quality of several perturbation theories against high-level quantum chemistry methods. First of all, we revisit the reference CCSD(T) ionization potentials for this popular benchmark set and establish a revised set of CCSD(T) results. Then, for all of these 100 molecules, we calculate the HOMO energy within second and third-order perturbation theory (PT2 and PT3), and, GW as post-Hartree-Fock methods. We found GW to be the most accurate of these three approximations for the ionization potential, by far. Going beyond GW by adding more diagrams is a tedious and dangerous activity: We tried to complement GW with second-order exchange (SOX), with second-order screened exchange (SOSEX), with interacting electron-hole pairs (W_TDHF), and with a GW density-matrix (γ^GW). Only the γ^GW result has a positive impact. Finally using an improved hybrid functional for the non-interacting Green’s function, considering it as a cheap way to approximate self-consistency, the accuracy of the simplest GW approximation improves even more. We conclude that GW is a miracle: Its subtle balance makes GW both accurate and fast.

Automated assessment of redox potentials for dyes in dye-sensitized photoelectrochemical cells

Sustainable solutions for hydrogen production, such as dye-sensitized photoelectrochemical cells (DS-PEC), rely on the fundamental properties of its components whose modularity allows for their separate investigation. In this work, we design and execute a high-throughput scheme to tune the ground state oxidation potential (GSOP) of perylene-type dyes by functionalizing them with different ligands. This allows us to identify promising candidates which can then be used to improve the cell's efficiency. First, we investigate the accuracy of different theoretical approaches by benchmarking them against experimentally determined GSOPs. We test different methods to calculate the vertical oxidation potential, including GW with different levels of self-consistency, Kohn-Sham (KS) orbital energies and total energy differences. We find that there is little difference in the performance of these methods. However, we show that it is crucial to take into account solvent effects as well as the structural relaxation of the dye after oxidation. Other thermodynamic contributions are negligible. Based on this benchmark, we decide on an optimal strategy, balancing computational cost and accuracy, to screen more than 1000 dyes and identify promising candidates which could be used to construct more robust DS-PECs.

A new generation of diagonal self-energies for the calculation of electron removal energies

A new generation of diagonal self-energy approximations in ab initio electron propagator theory for the calculation of electron removal energies of molecules and molecular ions has been derived from an intermediately normalized, Hermitized super-operator metric. These methods and widely used antecedents such as the outer valence Green's function and the approximately renormalized partial third order method are tested with respect to a dataset of vertical ionization energies generated with a valence, triple-ζ, correlation-consistent basis set and a converged series of many-body calculations whose accuracy approaches that of full configuration interaction. Several modifications of the diagonal second-order self-energy, a version of G₀W₀ theory based on Tamm-Dancoff excitations and several non-diagonal self-energies are also included in the tests. All new methods employ canonical Hartree-Fock orbitals. No adjustable or empirical parameters appear. A hierarchy of methods with optimal accuracy for a given level of computational efficiency is established. Several widely used diagonal self-energy methods are rendered obsolete by the new hierarchy whose members, in order of increasing accuracy, are (1) the opposite-spin non-Dyson diagonal second-order or os-nD-D2, (2) the approximately renormalized third-order quasiparticle or Q3+, (3) the renormalized third-order quasiparticle or RQ3, (4) the approximately renormalized linear third-order or L3+, and (5) the renormalized linear third-order or RL3 self-energies.

High-throughput virtual screening for organic electronics: a comparative study of alternative strategies

We present a review of the field of high-throughput virtual screening for organic electronics materials focusing on the sequence of methodological choices that determine each virtual screening protocol. These choices are present in all high-throughput virtual screenings and addressing them systematically will lead to optimised workflows and improve their applicability. We consider the range of properties that can be computed and illustrate how their accuracy can be determined depending on the quality and size of the experimental datasets. The approaches to generate candidates for virtual screening are also extremely varied and their relative strengths and weaknesses are discussed. The analysis of high-throughput virtual screening is almost never limited to the identification of top candidates and often new patterns and structure-property relations are the most interesting findings of such searches. The review reveals a very dynamic field constantly adapting to match an evolving landscape of applications, methodologies and datasets.

Exploring vibronic coupling in the benzene radical cation and anion with different levels of the GW approximation

The linear vibronic coupling constants of the benzene radical cation and anion have been obtained with different levels of the GW approximation, including G₀W₀, eigenvalue self-consistent GW, and quasiparticle self-consistent GW, as well as DFT with the following exchange-correlation functionals: BLYP, B3LYP, CAM-B3LYP, tuned CAM-B3LYP, and an IP-tuned CAM-B3LYP functional. The vibronic coupling constants were calculated numerically using the gradients of the eigenvalues of the degenerate HOMOs and LUMOs of the neutral benzene molecule for DFT, while the numerical gradients of the quasiparticle energies were used in the case of GW. The results were evaluated against those of high level wave function methods in the literature, and the approximate self-consistent GW methods and G₀W₀ with long-range corrected functionals were found to yield the best results on the whole.

Low-Order Scaling Quasiparticle Self-Consistent GW for Molecules

Low-order scaling GW implementations for molecules are usually restricted to approximations with diagonal self-energy. Here, we present an all-electron implementation of quasiparticle self-consistent GW for molecular systems. We use an efficient algorithm for the evaluation of the self-energy in imaginary time, from which a static non-local exchange-correlation potential is calculated via analytical continuation. By using a direct inversion of iterative subspace method, fast and stable convergence is achieved for almost all molecules in the GW100 database. Exceptions are systems which are associated with a breakdown of the single quasiparticle picture in the valence region. The implementation is proven to be starting point independent and good agreement of QP energies with other codes is observed. We demonstrate the computational efficiency of the new implementation by calculating the quasiparticle spectrum of a DNA oligomer with 1,220 electrons using a basis of 6,300 atomic orbitals in less than 4 days on a single compute node with 16 cores. We use then our implementation to study the dependence of quasiparticle energies of DNA oligomers consisting of adenine-thymine pairs on the oligomer size. The first ionization potential in vacuum decreases by nearly 1 electron volt and the electron affinity increases by 0.4 eV going from the smallest to the largest considered oligomer. This shows that the DNA environment stabilizes the hole/electron resulting from photoexcitation/photoattachment. Upon inclusion of the aqueous environment via a polarizable continuum model, the differences between the ionization potentials reduce to 130 meV, demonstrating that the solvent effectively compensates for the stabilizing effect of the DNA environment. The electron affinities of the different oligomers are almost identical in the aqueous environment.

GW100: A Slater-Type Orbital Perspective

We calculate complete basis set (CBS) limit-extrapolated ionization potentials (IPs) and electron affinities (EA) with Slater-type basis sets for the molecules in the GW100 database. To this end, we present two new Slater-type orbital (STO) basis sets of triple-(TZ) and quadruple-ζ (QZ) quality, whose polarization is adequate for correlated-electron methods and which contain extra diffuse functions to be able to correctly calculate EAs of molecules with a positive lowest unoccupied molecular orbital (LUMO). We demonstrate that going from TZ to QZ quality consistently reduces the basis set error of our computed IPs and EAs, and we conclude that a good estimate of these quantities at the CBS limit can be obtained by extrapolation. With mean absolute deviations (MAD) from 70 to 85 meV, our CBS limit-extrapolated IP are in good agreement with results from FHI-AIMS, TURBOMOLE, VASP, and WEST, while they differ by more than 130 meV on average from nanoGW. With a MAD of 160 meV, our EA are also in good agreement with the WEST code. Especially for systems with positive LUMOs, the agreement is excellent. With respect to other codes, the STO-type basis sets generally underestimate EAs of small molecules with strongly bound LUMOs. With 62 meV for IPs and 93 meV for EAs, we find much better agreement with CBS limit-extrapolated results from FHI-AIMS for a set of 250 medium to large organic molecules.

Exact Second-Order Corrections and Accurate Quasiparticle Energy Calculations in Density Functional Theory

We develop a second-order correction to commonly used density functional approximations (DFAs) to eliminate the systematic delocalization error. The method, based on the previously developed global scaling correction (GSC), is an exact quadratic correction to the DFA for the fractional charge behavior and uses the analytical second derivatives of the total energy with respect to fractional occupation numbers of the canonical molecular orbitals. For small and medium-size molecules, this correction leads to ground-state orbital energies that are a highly accurate approximation to the corresponding quasiparticle energies. It provides excellent predictions of ionization potentials, electron affinities, photoemission spectrum, and photoexcitation energies beyond previous approximate second-order approaches, thus showing potential for broad applications in computational spectroscopy.

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.

Intermediate state representation approach to physical properties of molecular electron-attached states: Theory, implementation, and benchmarking

Computational schemes for comprehensive studies of molecular electron-attached states and the calculation of electron affinities (EAs) are formulated and implemented employing the intermediate state representation (ISR) formalism and the algebraic-diagrammatic construction approximation for the electron propagator (EA-ADC). These EA-ADC(n)/ISR(m) schemes allow for a consistent treatment of not only electron affinities and pole strengths up to third-order of perturbation theory (n = 3) but also one-electron properties of electron-attached states up to second order (m = 2). The EA-ADC/ISR equations were implemented in the Q-Chem program for Ŝ_z-adapted intermediate states, allowing also open-shell systems to be studied using unrestricted Hartree-Fock references. For benchmarking of the EA-(U)ADC/ISR schemes, EAs and dipole moments of various electron-attached states of small closed- and open-shell molecules were computed and compared to full configuration interaction data. As an illustrative example, EA-ADC(3)/ISR(2) has been applied to the thymine-thymine (6-4) DNA photolesion.