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

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.

Relativistic Fully Self-Consistent GW for Molecules: Total Energies and Ionization Potentials

The fully self-consistent GW (scGW) method with an iterative solution of the Dyson equation provides a consistent approach for describing the ground and excited states without any dependence on the mean-field reference. In this work, we present a relativistic version of scGW for molecules containing heavy elements using the exact two-component (X2C) Coulomb approximation. We benchmark SOC-81 data set containing closed shell heavy elements for the first ionization potential using the fully self-consistent GW as well as one-shot GW. The self-consistent GW provides superior results compared to G0W0 with PBE reference and comparable results to G0W0 with PBE0 while also removing the starting point dependence. The photoelectron spectra obtained at the X2C level demonstrate very good agreement with the experimental spectra. We also observe that scGW provides very good estimation of ionization potential for the inner d-shell orbitals. Additionally, using the well-conserved total energy, we investigate the equilibrium bond length and harmonic frequencies of a few halogen dimers using scGW. Overall, our findings demonstrate the applicability of the fully self-consistent GW method for accurate ionization potential, photoelectron spectra, and total energies in finite systems with heavy elements with a reasonable computational scaling.

Efficient Full-Frequency GW Calculations Using a Lanczos Method

The GW approximation is widely used for reliable and accurate modeling of single-particle excitations. It also serves as a starting point for many theoretical methods, such as its use in the Bethe-Salpeter equation (BSE) and dynamical mean-field theory. However, full-frequency GW calculations for large systems with hundreds of atoms remain computationally challenging, even after years of efforts to reduce the prefactor and improve scaling. We propose a method that reformulates the correlation part of the GW self-energy as a resolvent of a Hermitian matrix, which can be efficiently and accurately computed using the standard Lanczos method. This method enables full-frequency GW calculations of material systems with a few hundred atoms on a single computing workstation. We further demonstrate the efficiency of the method by calculating the defect-state energies of silicon quantum dots with diameters up to 4 nm and nearly 2,000 silicon atoms using only 20 computational nodes.

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.

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.

Spatial Decay and Limits of Quantum Solute-Solvent Interactions

Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.

Slater transition methods for core-level electron binding energies

Methods for computing core-level ionization energies using self-consistent field (SCF) calculations are evaluated and benchmarked. These include a "full core hole" (or "ΔSCF") approach that fully accounts for orbital relaxation upon ionization, but also methods based on Slater's transition concept in which the binding energy is estimated from an orbital energy level that is obtained from a fractional-occupancy SCF calculation. A generalization that uses two different fractional-occupancy SCF calculations is also considered. The best of the Slater-type methods afford mean errors of 0.3-0.4 eV with respect to experiment for a dataset of K-shell ionization energies, a level of accuracy that is competitive with more expensive many-body techniques. An empirical shifting procedure with one adjustable parameter reduces the average error below 0.2 eV. This shifted Slater transition method is a simple and practical way to compute core-level binding energies using only initial-state Kohn-Sham eigenvalues. It requires no more computational effort than ΔSCF and may be especially useful for simulating transient x-ray experiments where core-level spectroscopy is used to probe an excited electronic state, for which the ΔSCF approach requires a tedious state-by-state calculation of the spectrum. As an example, we use Slater-type methods to model x-ray emission spectroscopy.

Efficient All-Electron Time-Dependent Density Functional Theory Calculations Using an Enriched Finite Element Basis

We present an efficient and systematically convergent approach to all-electron real-time time-dependent density functional theory (TDDFT) calculations using a mixed basis, termed as enriched finite element (EFE) basis. The EFE basis augments the classical finite element basis (CFE) with a compactly supported numerical atom-centered basis, obtained from atomic ground-state DFT calculations. Particularly, we orthogonalize the enrichment functions with respect to the classical finite element basis to ensure good conditioning of the resultant basis. We employ the second-order Magnus propagator in conjunction with an adaptive Krylov subspace method for efficient time evolution of the Kohn-Sham orbitals. We rely on a priori error estimates to guide our choice of an adaptive finite element mesh as well as the time step to be used in the TDDFT calculations. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretizations. Notably, we attain a 50-100 times speedup for the EFE basis over the CFE basis. We also demonstrate the efficacy of the EFE basis for both linear and nonlinear responses by studying the absorption spectra in sodium clusters, the linear to nonlinear response transition in the green fluorescence protein chromophore, and the higher harmonic generation in the magnesium dimer. Lastly, we attain good parallel scalability of our numerical implementation of the EFE basis for up to ∼1000 processors, using a benchmark system of a 50-atom sodium nanocluster.

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. We present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

Basis Set Selection for Molecular Core-Level GW Calculations

The GW approximation has been recently gaining popularity among the methods for simulating molecular core-level X-ray photoemission spectra. Traditionally, Gaussian-type orbital GW core-level binding energies have been computed using either the cc-pVnZ or def2-nZVP basis set families, extrapolating the obtained results to the complete basis set limit, followed by an element-specific relativistic correction. Despite achieving rather good accuracy, it has been previously stated that these binding energies are chronically underestimated. In the present work, we show that those previous studies obtained results that were not well-converged with respect to the basis set size and that, once basis set convergence is achieved, there seems to be no such underestimation. Standard techniques known to offer a good cost-accuracy ratio in other theories demonstrate that the cc-pVnZ and def2-nZVP families exhibit contraction errors and might lead to unreliable complete basis set extrapolations for absolute binding energies, often deviating about 200-500 meV from the putative basis set limit found in this work. On the other hand, uncontracted versions of these basis sets offer vastly improved convergence. Even faster convergence can be obtained using core-rich property-optimized basis set families like pcSseg-n, pcJ-n, and ccX-nZ. Finally, we also show that the improvement observed for core properties using these specialized basis sets does not degrade their description of valence excitations: vertical ionization potentials and electron affinities computed with these basis sets converge as fast as the ones obtained with the aug-cc-pVnZ family, thus offering a balanced description of both core and valence regions.

Quantum embedding theories to simulate condensed systems on quantum computers

Quantum computers hold promise to improve the efficiency of quantum simulations of materials and to enable the investigation of systems and properties that are more complex than tractable at present on classical architectures. Here, we discuss computational frameworks to carry out electronic structure calculations of solids on noisy intermediate-scale quantum computers using embedding theories, and we give examples for a specific class of materials, that is, solid materials hosting spin defects. These are promising systems to build future quantum technologies, such as quantum computers, quantum sensors and quantum communication devices. Although quantum simulations on quantum architectures are in their infancy, promising results for realistic systems appear to be within reach.

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).

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.

Combined First-Principles Calculations of Electron-Electron and Electron-Phonon Self-Energies in Condensed Systems

We present a method to efficiently combine the computation of electron-electron and electron-phonon self-energies, which enables the evaluation of electron-phonon coupling at the G₀W₀ level of theory for systems with hundreds of atoms. In addition, our approach, which is a generalization of a method recently proposed for molecules [J. Chem. Theory Comput. 2018, 14, 6269-6275], enables the inclusion of nonadiabatic and temperature effects at no additional computational cost. We present results for diamond and defects in diamond and discuss the importance of numerically accurate G₀W₀ band structures to obtain robust predictions of zero point renormalization (ZPR) of band gaps, and of the inclusion of nonadiabatic effects to accurately compute the ZPR of defect states in the band gap.

Integrating Computation and Experiment to Investigate Photoelectrodes for Solar Water Splitting at the Microscopic Scale

ConspectusPhotoelectrochemical water-splitting is a promising and sustainable way to store the energy of the sun in chemical bonds and use it to produce hydrogen gas, a clean fuel. The key components in photoelectrochemical cells (PECs) are photoelectrodes, including a photocathode that reduces water to hydrogen gas and a photoanode that oxidizes water to oxygen gas. Materials used in photoelectrodes for PECs must effectively absorb sunlight, yield photogenerated carriers, and exhibit electronic properties that enable the efficient shuttling of carriers to the surface to participate in relevant water-splitting reactions. Discovering and understanding the key characteristics of optimal photoelectrode materials is paramount to the realization of PEC technologies.Oxide-based photoelectrodes can satisfy many of these materials requirements, including stability in aqueous environments, band edges with reasonable alignment with the redox potentials for water splitting, and ease of synthesis. However, oxide photoelectrodes generally suffer from poor charge transport properties and considerable bulk electron-hole separation, and they have relatively large band gaps. Numerous strategies have been proposed to improve these aspects and understand how these improvements are reflected in the photoelectrochemical performance. Unfortunately, the structural and compositional complexity of multinary oxides accompanied by the inherent complexity of photoelectrochemical processes makes it challenging to understand the individual effects of composition, structure, and defects in the bulk and on the surface on a material's photoelectrochemical properties. The integration of experiment and theory has great potential to increase our atomic-level understanding of structure-composition-property relationships in oxide photoelectrodes.In this Account, we describe how integrating experiment and theory is beneficial for achieving scientific insights at the microscopic scale. We highlight studies focused on understanding the role of (i) bulk composition via solid-state solutions, intercalation, and comparison with isoelectronic compounds, (ii) dopants for both the anion and cation and their interactions with oxygen vacancies, and (iii) surface/interface structure in the photocurrent generation and photoelectrochemical performance in oxide photoelectrodes. In each instance, we outline strategies and considerations for integrating experiment and theory and describe how this integration led to valuable insights and new directions in uncovering structure-composition-property relationships. Our aim is to demonstrate the unique value of combining experiment and theory in studying photoelectrodes and to encourage the continued effort to bring experiment and theory in closer step with each other.

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.

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.

Photoelectron spectra of early 3d-transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3d-transition metal dioxide anions, ScO₂ ^-, TiO₂ ^-, VO₂ ^-, CrO₂ ^-, and MnO₂ ^-, are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point (G₀W₀@PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV-0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew-Burke-Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G₀W₀@PBE0, the G-only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3d-transition metal dioxide anions. Adding eigenvalue self-consistency also in W, on the other hand, worsens the agreement with experiment overall. Our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe-Salpeter equation (BSE) for electron-hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

Quantum Information and Algorithms for Correlated Quantum Matter

Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.

Density Functional Prediction of Quasiparticle, Excitation, and Resonance Energies of Molecules With a Global Scaling Correction Approach

Molecular quasiparticle and excitation energies determine essentially the spectral characteristics measured in various spectroscopic experiments. Accurate prediction of these energies has been rather challenging for ground-state density functional methods, because the commonly adopted density function approximations suffer from delocalization error. In this work, by presuming a quantitative correspondence between the quasiparticle energies and the generalized Kohn–Sham orbital energies, and employing a previously developed global scaling correction approach, we achieve substantially improved prediction of molecular quasiparticle and excitation energies. In addition, we also extend our previous study on temporary anions in resonant states, which are associated with negative molecular electron affinities. The proposed approach does not require any explicit self-consistent field calculation on the excited-state species, and is thus highly efficient and convenient for practical purposes.

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.

Spin-Triplet Excitonic Insulator: The Case of Semihydrogenated Graphene

While various excitonic insulators have been studied in the literature, due to the perceived too-small spin splitting, spin-triplet excitonic insulator is rare. In two-dimensional systems such as a semihydrogenated graphene (known as graphone), however, it is possible, as revealed by first-principles calculations coupled with Bethe-Salpeter equation. The critical temperature, given by an effective Hamiltonian, is 11.5 K. While detecting excitonic insulators is still a daunting challenge, the condensation of triplet excitons will result in spin superfluidity, which can be directly measured by a transport experiment. Nonlocal dielectric screening also leads to an unexpected phenomenon, namely, an indirect-to-direct transition crossover between single-particle band and exciton dispersion in the semihydrogenated graphene, which offers yet another test by experiment.

Approximate equation-of-motion coupled-cluster methods for electron affinities of closed-shell molecules

We investigate performance of the equation-of-motion coupled-cluster method at the single and doubles level (EOM-CCSD) and a series of approximate methods based on EOM-CCSD on electron affinities (EA) of closed-shell cations and neutral molecules with positive and negative EAs in this work. Our results confirm that P-EOM-MBPT2 can provide reasonable EAs when molecules with significant multireference character are not considered and its mean absolute error on EAs of these molecules is around or less than 0.2 eV. Its accuracy is comparable to that of the more expensive EOM-CCSD(2) method. Results of EOM-CCSD(2), P-EOM-MBPT2, and CIS(D_∞) indicate that the [[H, a_c ⁺], T₂] term in the 1h2p-1h block is more important on EAs than the term neglected in the 1h2p-1h2p block in P-EOM-MBPT2. We proposed an economical method where EAs from CIS(D_∞) are corrected by treating this [[H, a_c ⁺], T₂] term in the 1h2p-1h block perturbatively [corr-CIS(D_∞)]. EAs with corr-CIS(D_∞) agree very well with those of P-EOM-MBPT2 with a difference of less than 0.02 eV. Computational scaling of this method is N⁴ for the iterative part and N⁵ for some non-iterative steps. Its storage requirement is only of OV³. Corr-CIS(D_∞) is an economical and reliable method on EAs, and it can be applied to EAs of large molecules.

Ab Initio Many-Body Perturbation Theory Calculations of the Electronic and Optical Properties of Cyclometalated Ir(III) Complexes

Cyclometalated Ir(III) compounds are the preferred choice as organic emitters in organic light-emitting diodes. In practice, the presence of the transition metal surrounded by carefully designed ligands allows fine-tuning of the emission frequency as well as good efficiency of the device. To support the development of new compounds, experimental measurements are generally compared with absorption and emission spectra obtained from ab initio calculations. The standard approach for these calculations is time-dependent density functional theory (TDDFT) with a hybrid exchange-correlation functional like B3LYP. Because of the size of these compounds, the application of more complex quantum chemistry approaches can be challenging. In this work, we used many-body perturbation theory approaches, in particular the GW approximation with the Bethe-Salpeter equation (BSE) implemented in Gaussian basis sets, to calculate the quasiparticle properties and the absorption spectra of six cyclometalated Ir(III) complexes, going beyond TDDFT. In the presented results, we compared standard TDDFT simulations with BSE calculations performed on top of perturbative G₀W₀ and accounting for eigenvalue self-consistency. Moreover, in order to investigate in detail the effect of the DFT starting point, we concentrated on Ir(ppy)₃ and performed GW-BSE simulations starting from different DFT exchange-correlation potentials.

Improving the efficiency of G0W0 calculations with approximate spectral decompositions of dielectric matrices

Recently, it was shown that the calculation of quasiparticle energies using the G₀W₀ approximation can be performed without computing explicitly any virtual electronic states, by expanding the Green function and screened Coulomb interaction in terms of the eigenstates of the static dielectric matrix. Avoiding the evaluation of virtual electronic states leads to improved efficiency and ease of convergence of G₀W₀ calculations. Here, we propose a further improvement of the efficiency of these calculations, based on an approximation of density-density response functions of molecules and solids. The approximation relies on the calculation of a subset of eigenvectors of the dielectric matrix using the kinetic operator instead of the full Hamiltonian, and it does not lead to any substantial loss of accuracy for the quasiparticle energies. The computational savings introduced by this approximation depend on the system, and they become more substantial as the number of electrons increases.

Practical GW scheme for electronic structure of 3d-transition-metal monoxide anions: ScO-, TiO-, CuO-, and ZnO. J Chem Phys 2019;151:134305. [PMID: 31594362 DOI: 10.1063/1.5118671]

The GW approximation to many-body perturbation theory is a reliable tool for describing charged electronic excitations, and it has been successfully applied to a wide range of extended systems for several decades using a plane-wave basis. However, the GW approximation has been used to test limited spectral properties of a limited set of finite systems (e.g., frontier orbital energies of closed-shell sp molecules) only for about a decade using a local-orbital basis. Here, we calculate the quasiparticle spectra of closed- and open-shell molecular anions with partially and completely filled 3d shells (shallow and deep 3d states, respectively), ScO^-, TiO^-, CuO^-, and ZnO^-, using various levels of GW theory, and compare them to experiments to evaluate the performance of the GW approximation on the electronic structure of small molecules containing 3d transition metals. We find that the G-only eigenvalue self-consistent GW scheme with W fixed to the PBE level (G_nW₀@PBE), which gives the best compromise between accuracy and efficiency for solids, also gives good results for both localized (d) and delocalized (sp) states of 3d-transition-metal oxide molecules. The success of G_nW₀@PBE in predicting electronic excitations in these systems reasonably well is likely due to the fortuitous cancellation effect between the overscreening of the Coulomb interaction by PBE and the underscreening by the neglect of vertex corrections. Together with the absence of the self-consistent field convergence error (e.g., spin contamination in open-shell systems) and the GW multisolution issue, the G_nW₀@PBE scheme gives the possibility to predict the electronic structure of complex real systems (e.g., molecule-solid and sp-d hybrid systems) accurately and efficiently.

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.

Finite-Field Approach to Solving the Bethe-Salpeter Equation

We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation and the calculation of the screened Coulomb interaction in a finite field. The method does not require either the explicit evaluation of dielectric matrices or of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition, it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wave functions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.

Koopmans-Compliant Functionals and Potentials and Their Application to the GW100 Test Set

Koopmans-compliant (KC) functionals have been shown to provide accurate spectral properties through a generalized condition of piecewise linearity of the total energy as a function of the fractional addition/removal of an electron to/from any orbital. We analyze the performance of different KC functionals on a large and standardized set of 100 molecules, the GW100 test set, comparing vertical ionization potentials (taken as opposite of the orbital energies) to those obtained from accurate quantum chemistry methods, and to experimental results. We find excellent agreement, with a mean absolute error of 0.20 eV for the KIPZ functional on the first ionization potential, which is state-of-the-art for both density functional theory (DFT)-based calculations and many-body perturbation theory. We highlight similarities and differences between KC functionals and other electronic-structure approaches, such as dielectric-dependent hybrid functionals and Green's function methods, both from a theoretical and from a practical point of view, arguing that KC potentials can be considered as local and orbital-dependent approximations to the electronic self-energy, already including approximate vertex corrections.

Coupling First-Principles Calculations of Electron-Electron and Electron-Phonon Scattering, and Applications to Carbon-Based Nanostructures

We report first-principles calculations of electronic gaps, lifetimes, and photoelectron spectra of a series of molecules, performed by efficiently combining the computation of electron-electron and electron-phonon self-energies. The dielectric matrix is represented in terms of dielectric eigenpotentials, utilized for both the calculation of G₀ W₀ quasi-particle energies and the diagonalization of the dynamical matrix; virtual electronic states are never explicitly computed and all self-energies are evaluated over the full frequency spectrum. Our formulation enables electronic structure calculations at the many-body perturbation theory level, inclusive of electron-phonon coupling, for systems with hundreds of electrons.

Green's Function Coupled-Cluster Approach: Simulating Photoelectron Spectra for Realistic Molecular Systems

In this paper, we present an efficient implementation for the analytical energy-dependent Green's function coupled-cluster with singles and doubles (GFCCSD) approach with our first practice being computing spectral functions of realistic molecular systems. Because of its algebraic structure, the presented method is highly scalable and is capable of computing spectral function for a given molecular system in any energy region. Several typical examples have been given to demonstrate its capability of computing spectral functions not only in the valence band but also in the core-level energy region. Satellite peaks have been observed in the inner valence band and core-level energy region where a many-body effect becomes significant and the single particle picture of ionization often breaks down. The accuracy test has been carried out by extensively comparing the computed spectral functions by our GFCCSD method with experimental photoelectron spectra as well as the theoretical ionization potentials obtained from other methods. It turns out the GFCCSD method is able to provide a qualitative or semiquantitative level of description of ionization processes in both the core and valence regimes. To significantly improve the GFCCSD results for the main ionic states, a larger basis set can usually be employed, whereas the improvement of the GFCCSD results for the satellite states needs higher-order many-body terms to be included in the GFCC implementation.