The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

Cumulant Green's function methods for molecules

The cumulant expansion of the Green's function is a computationally efficient beyond-GW approach renowned for its significant enhancement of satellite features in materials. In contrast to the ubiquitous GW approximation of many-body perturbation theory, ab initio cumulant expansions performed on top of GW (GW + C) have demonstrated the capability to handle multi-particle processes by incorporating higher-order correlation effects or vertex corrections, yielding better agreements between experiment and theory for satellite structures. While widely employed in condensed matter physics, very few applications of GW + C have been published on molecular systems. Here, we assess the performance of this scheme on a series of 10-electron molecular systems (Ne, HF, H2O, NH3, and CH4) where full configuration interaction estimates of the outer-valence quasiparticle and satellite energies are available.

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.

Proton Conducting Neuromorphic Materials and Devices

Neuromorphic computing and artificial intelligence hardware generally aims to emulate features found in biological neural circuit components and to enable the development of energy-efficient machines. In the biological brain, ionic currents and temporal concentration gradients control information flow and storage. It is therefore of interest to examine materials and devices for neuromorphic computing wherein ionic and electronic currents can propagate. Protons being mobile under an external electric field offers a compelling avenue for facilitating biological functionalities in artificial synapses and neurons. In this review, we first highlight the interesting biological analog of protons as neurotransmitters in various animals. We then discuss the experimental approaches and mechanisms of proton doping in various classes of inorganic and organic proton-conducting materials for the advancement of neuromorphic architectures. Since hydrogen is among the lightest of elements, characterization in a solid matrix requires advanced techniques. We review powerful synchrotron-based spectroscopic techniques for characterizing hydrogen doping in various materials as well as complementary scattering techniques to detect hydrogen. First-principles calculations are then discussed as they help provide an understanding of proton migration and electronic structure modification. Outstanding scientific challenges to further our understanding of proton doping and its use in emerging neuromorphic electronics are pointed out.

Quantum-Quantum and Quantum-Quantum-Classical Schemes for Near-Gap Excitations with Projection-Based-Embedded GW-Bethe-Salpeter Equation

We present quantum-quantum and quantum-quantum-classical schemes based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-BSE) employing projection-based-embedding (PbE). Such approaches allow defining active and inactive subsystems of larger, complex molecular systems, with only the smaller active subsystem being explicitly treated by GW-BSE offering significant computational advantages. However, as PbE can modify the single-particle states in the Kohn-Sham (KS) ground state calculation and screening effects from the inactive region are not automatically included in GW-BSE, results from such PbE-GW-BSE calculations can deviate from a full-system reference. Here, we scrutinize in detail, e.g., the individual and combined effects of different choices of active regions, the influence of omitting the screening from the inactive region, and strategies for basis set truncation on frontier orbital and near-gap electron-hole excitation energies. As prototypical systems, we consider a diketopyrrolopyrrole bicyclic ring including side-chains, a polarity-sensitive dye (prodan) in aqueous environment, and a π-stacked dimer of benzene and tetracyanoethylene in water, respectively, covering a variety of excitation characters in molecular systems with complex chemical environments and photoinduced processes. Our results suggest that to obtain agreement of approximately 0.1 eV between near-gap excitation energies from embedded and full calculations, the active region should be chosen based on the Mulliken population of the full highest-occupied molecular orbital and that careful benchmarking should be done on the KS level before the actual GW-BSE steps when basis set truncation is used. We find that PbE-GW-BSE offers significant reductions in computation times and, more importantly, memory requirements, making calculations for considerably larger systems tractable.

Chemical Potentials and the One-Electron Hamiltonian of the Second-Order Perturbation Theory from the Functional Derivative Approach

We develop a functional derivative approach to calculate the chemical potentials of second-order perturbation theory (MP2). In the functional derivative approach, the correlation part of the MP2 chemical potential, which is the derivative of the MP2 correlation energy with respect to the occupation number of frontier orbitals, is obtained from the chain rule via the noninteracting Green's function. First, the MP2 correlation energy is expressed in terms of the noninteracting Green's function, and its functional derivative to the noninteracting Green's function is the second-order self-energy. Then, the derivative of the noninteracting Green's function to the occupation number is obtained by including the orbital relaxation effect. We show that the MP2 chemical potentials obtained from the functional derivative approach agree with that obtained from the finite difference approach. The one-electron Hamiltonian, defined as the derivative of the MP2 energy with respect to the one particle density matrix, is also derived using the functional derivative approach, which can be used in the self-consistent calculations of MP2 and double-hybrid density functionals. The developed functional derivative approach is promising for calculating the chemical potentials and the one-electron Hamiltonian of approximate functionals and many-body perturbation approaches dependent explicitly on the noninteracting Green's function.

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.

Electronic Couplings and Conversion Dynamics between Localized and Charge Transfer Excitations from Many-Body Green's Functions Theory

We investigate the determination of electronic coupling between localized excitations (LEs) and charge-transfer (CT) excitations based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-BSE). Using a small molecule dimer system, we first study the influence of different diabatization methods, as well as different model choices within GW-BSE, such as the self-energy models or different levels of self-consistency, and find that these choices affect the LE-CT couplings only minimally. We then consider a large-scale low-donor morphology formed from rubrene and fullerene and evaluate the LE-CT couplings based on coupled GW-BSE-molecular mechanics calculations. For these disordered systems of bulky molecules, we observe differences in the couplings based on the Edmiston-Ruedenberg diabatization compared to the more approximate Generalize Mulliken-Hush and fragment charge difference diabatization formalisms. In a kinetic model for the conversion between LE and CT states, these differences affect the details of state populations in an intermediate time scale but not the final populations.

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.

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.

The NOMAD mini-apps: A suite of kernels from ab initio electronic structure codes enabling co-design in high-performance computing

This article introduces a suite of mini-applications (mini-apps) designed to optimise computational kernels in ab initio electronic structure codes. The suite is developed from flagship applications participating in the NOMAD Center of Excellence, such as the ELPA eigensolver library and the GW implementations of the exciting, Abinit, and FHI-aims codes. The mini-apps were identified by targeting functions that significantly contribute to the total execution time in the parent applications. This strategic selection allows for concentrated optimisation efforts. The suite is designed for easy deployment on various High-Performance Computing (HPC) systems, supported by an integrated CMake build system for straightforward compilation and execution. The aim is to harness the capabilities of emerging (post)exascale systems, which necessitate concurrent hardware and software development - a concept known as co-design. The mini-app suite serves as a tool for profiling and benchmarking, providing insights that can guide both software optimisation and hardware design. Ultimately, these developments will enable more accurate and efficient simulations of novel materials, leveraging the full potential of exascale computing in material science research.

Ab Initio Modeling of Mixed-Dimensional Heterostructures: A Path Forward

Understanding the electronic structure of mixed-dimensional heterostructures is essential for maximizing their application potential. However, accurately modeling such interfaces is challenging due to the complex interplay between the subsystems. We employ a computational framework integrating first-principles methods, including GW, density functional theory (DFT), and the polarizable continuum model, to elucidate the electronic structure of mixed-dimensional heterojunctions formed by free-base phthalocyanines and monolayer molybdenum disulfide. We assess the impact of dielectric screening across various scenarios, from isolated molecules to organic films on a substrate-supported monolayer. Our findings show that while polarization effects cause significant renormalization of molecular energy levels, band energies and alignments in the most relevant setup can be accurately predicted through DFT simulations of the individual subsystems. Additionally, we analyze orbital hybridization, revealing potential pathways for interfacial charge transfer. This study offers new insights into hybrid inorganic/organic interfaces and provides a practical computational protocol suitable for scaled-up studies.

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.

AB-G0W0: A practical G0W0 method without frequency integration based on an auxiliary boson expansion

Common G0W0 implementations rely on numerical or analytical frequency integration to determine the G0W0 self-energy, which results in a variety of practical complications. Recently, we have demonstrated an exact connection between the G0W0 approximation and equation-of-motion quantum chemistry approaches [J. Tölle and G. Kin-Lic Chan, J. Chem. Phys. 158, 124123 (2023)]. Based on this connection, we propose a new method to determine G0W0 quasiparticle energies, which completely avoids frequency integration and its associated problems. To achieve this, we make use of an auxiliary boson (AB) expansion. We name the new approach AB-G0W0 and demonstrate its practical applicability in a range of molecular problems.

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.

The effect of particle size on the optical and electronic properties of hydrogenated silicon nanoparticles

We use a combination of many-body perturbation theory and time-dependent density functional theory to study the optical and electronic properties of hydrogen terminated silicon nanoparticles. We predict that the lowest excited states of these silicon nanoparticles are excitonic in character and that the corresponding excitons are completely delocalised over the volume of the particle. The size of the excitons is predicted to increase proportionally with the particle size. Conversely, we predict that the fundamental gap, the optical gap, and the exciton binding energy increase with decreasing particle size. The exciton binding energy is predicted to counter-act the variation in the fundamental gap and hence to reduce the variation of the optical gap with particle size. The variation in the exciton binding energy itself is probably caused by a reduction in the dielectric screening with decreasing particle size. The intensity of the excited state corresponding to the optical gap and other low energy excitations are predicted to increase with decreasing particle size. We explain this increase in terms of the 'band structure' becoming smeared out in reciprocal space with decreasing particle size, increasing the 'overlap' between the occupied and unoccupied quasiparticle states and thus, the oscillator strength. Fourier transforms of the lowest excitons show that they inherit the periodicity of the frontier quasiparticle states. This, combined with the delocalisation of the exciton and the large exciton binding energy, means that the excitons in silicon nanoparticles combine aspects of Wannier-Mott, delocalisation and effect of periodicity of the underlying structure, and Frenkel, large exciton binding energy, excitons.

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.

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.

Accelerating Analytic-Continuation GW Calculations with a Laplace Transform and Natural Auxiliary Functions

We present a simple and accurate GW implementation based on a combination of a Laplace transform (LT) and other acceleration techniques used in post-self-consistent field quantum chemistry, namely, natural auxiliary functions and the frozen-core approximation. The LT-GW approach combines three major benefits: (a) a small prefactor for computational scaling, (b) easy integration into existing molecular GW implementations, and (c) significant performance improvements for a wide range of possible applications. Illustrating these advantages for systems consisting of up to 352 atoms and 7412 basis functions, we further demonstrate the benefits of this approach combined with an efficient implementation of the Bethe-Salpeter equation.

Low-Scaling GW Algorithm Applied to Twisted Transition-Metal Dichalcogenide Heterobilayers

The GW method is widely used for calculating the electronic band structure of materials. The high computational cost of GW algorithms prohibits their application to many systems of interest. We present a periodic, low-scaling, and highly efficient GW algorithm that benefits from the locality of the Gaussian basis and the polarizability. The algorithm enables G0W0 calculations on a MoSe2/WS2 bilayer with 984 atoms per unit cell, in 42 h using 1536 cores. This is 4 orders of magnitude faster than a plane-wave G0W0 algorithm, allowing for unprecedented computational studies of electronic excitations at the nanoscale.

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Electronic and optical properties of core-shell InAlN nanorods: a comparative study via LDA, LDA-1/2, mBJ, HSE06, G0W0 and BSE methods

Currently, self-induced InAlN core-shell nanorods enjoy an advanced stage of accumulation of experimental data from their growth and characterization as well as a comprehensive understanding of their formation mechanism by the ab initio modeling based on Synthetic Growth Concept. However, their electronic and optical properties, on which most of their foreseen applications are expected to depend, have not been investigated comprehensively. GW and the Bethe-Salpeter equation (BSE) are regarded as the state-of-the-art ab initio methodologies to study these properties. However, one of the major drawbacks of these methods is the computational cost, much higher than density-functional theory (DFT). Therefore, in many applications, it is highly desirable to answer the question of how well approaches based on DFT, such as e.g. the local density approximation (LDA), LDA-1/2, the modified Becke-Johnson (mBJ) and the Heyd-Scuseria-Ernzerhof (HSE06) functionals, can be employed to calculate electronic and optical properties with reasonable accuracy. In the present paper, we address this question, investigating how effective the DFT-based methodologies LDA, LDA-1/2, mBJ and HSE06 can be used as approximate tools in studies of the electronic and optical properties of scaled down models of core-shell InAlN nanorods, thus, avoiding GW and BSE calculations.

Excitons in metal-halide perovskites from first-principles many-body perturbation theory

Metal-halide perovskites are a structurally, chemically, and electronically diverse class of semiconductors with applications ranging from photovoltaics to radiation detectors and sensors. Understanding neutral electron-hole excitations (excitons) is key for predicting and improving the efficiency of energy-conversion processes in these materials. First-principles calculations have played an important role in this context, allowing for a detailed insight into the formation of excitons in many different types of perovskites. Such calculations have demonstrated that excitons in some perovskites significantly deviate from canonical models due to the chemical and structural heterogeneity of these materials. In this Perspective, I provide an overview of calculations of excitons in metal-halide perovskites using Green's function-based many-body perturbation theory in the GW + Bethe-Salpeter equation approach, the prevalent method for calculating excitons in extended solids. This approach readily considers anisotropic electronic structures and dielectric screening present in many perovskites and important effects, such as spin-orbit coupling. I will show that despite this progress, the complex and diverse electronic structure of these materials and its intricate coupling to pronounced and anharmonic structural dynamics pose challenges that are currently not fully addressed within the GW + Bethe-Salpeter equation approach. I hope that this Perspective serves as an inspiration for further exploring the rich landscape of excitons in metal-halide perovskites and other complex semiconductors and for method development addressing unresolved challenges in the field.

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.

Efficient Quasiparticle Determination beyond the Diagonal Approximation via Random Compression

Calculations of excited states in the Green's function formalism often invokes the diagonal approximation, in which the quasiparticle states are taken from a mean-field calculation. In this paper, we extend the stochastic approaches applied in the many-body perturbation theory and overcome this limitation for large systems in which we are interested in a small subset of states. We separate the problem into a core subspace whose coupling to the remainder of the system environment is stochastically sampled. This method is exemplified on computing hole injection energies into CO2 on an extended gold surface with nearly 3000 electrons. We find that in the extended system the size of the problem can be compressed up to 95% using stochastic sampling. This result provides a way forward for self-consistent stochastic methods and determination of Dyson orbitals in large systems.

Benchmarking the accuracy of the separable resolution of the identity approach for correlated methods in the numeric atom-centered orbitals framework

Four-center two-electron Coulomb integrals routinely appear in electronic structure algorithms. The resolution-of-the-identity (RI) is a popular technique to reduce the computational cost for the numerical evaluation of these integrals in localized basis-sets codes. Recently, Duchemin and Blase proposed a separable RI scheme [J. Chem. Phys. 150, 174120 (2019)], which preserves the accuracy of the standard global RI method with the Coulomb metric and permits the formulation of cubic-scaling random phase approximation (RPA) and GW approaches. Here, we present the implementation of a separable RI scheme within an all-electron numeric atom-centered orbital framework. We present comprehensive benchmark results using the Thiel and the GW100 test set. Our benchmarks include atomization energies from Hartree-Fock, second-order Møller-Plesset (MP2), coupled-cluster singles and doubles, RPA, and renormalized second-order perturbation theory, as well as quasiparticle energies from GW. We found that the separable RI approach reproduces RI-free HF calculations within 9 meV and MP2 calculations within 1 meV. We have confirmed that the separable RI error is independent of the system size by including disordered carbon clusters up to 116 atoms in our benchmarks.

Machine Learning Many-Body Green's Functions for Molecular Excitation Spectra

We present a machine learning (ML) framework for predicting Green's functions of molecular systems, from which photoemission spectra and quasiparticle energies at quantum many-body level can be obtained. Kernel ridge regression is adopted to predict self-energy matrix elements on compact imaginary frequency grids from static and dynamical mean-field electronic features, which gives direct access to real-frequency many-body Green's functions through analytic continuation and Dyson's equation. Feature and self-energy matrices are represented in a symmetry-adapted intrinsic atomic orbital plus projected atomic orbital basis to enforce rotational invariance. We demonstrate good transferability and high data efficiency of the proposed ML method across molecular sizes and chemical species by showing accurate predictions of density of states (DOS) and quasiparticle energies at the level of many-body perturbation theory (GW) or full configuration interaction. For the ML model trained on 48 out of 1995 molecules randomly sampled from the QM7 and QM9 data sets, we report the mean absolute errors of ML-predicted highest occupied and lowest unoccupied molecular orbital energies to be 0.13 and 0.10 eV, respectively, compared to GW@PBE0. We further showcase the capability of this method by applying the same ML model to predict DOS for significantly larger organic molecules with up to 44 heavy atoms.

G0W0 Ionization Potentials of First-Row Transition Metal Aqua Ions

We report computations of the vertical ionization potentials within the GW approximation of the near-complete series of first-row transition metal (V-Cu) aqua ions in their most common oxidation states, i.e., V3+, Cr3+, Cr2+, Mn2+, Fe3+, Fe2+, Co2+, Ni2+, and Cu2+. The d-orbital occupancy of these systems spans a broad range from d2 to d9. All of the structures were first optimized at the density functional theory level using a large cluster of explicit water molecules that are embedded in a continuum solvation model. Vertical ionization potentials were computed with the one-shot G0W0 approach on a range of transition metal ion clusters (6, 18, 40, and 60 explicit water molecules), wherein the convergence with respect to the basis set size was evaluated using the systems with 40 water molecules. We assess the results using three different density functional approximations as starting points for the vertical ionization potential calculations, namely, G0W0@PBE, G0W0@PBE0, and G0W0@r2SCAN. While the predicted ground-state structures are similar to all three exchange-correlation functionals, the vertical ionization potentials were in closer agreement with experiment when using the G0W0@PBE0 and G0W0@r2SCAN approaches, with the r2SCAN-based calculations being significantly less expensive. Computed bond distances and vertical ionization potentials for all structures are in good agreement with available experimental data.

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

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

Excess and excited-state dipole moments of real-life dyes: a comparison between wave-function, BSE/GW, and TD-DFT values

In this work, we assess the accuracy of the Bethe-Salpeter equation (BSE) many-body Green's function formalism, adopting the eigenvalue-self-consistent evGW exchange-correlation kernel, for the calculation of the excited-state (μES) and excess dipole moments (Δμ), the latter ones being the changes of dipole amplitude between the ground and excited states (ES), in organic dyes. We compare the results obtained with wave-function methods [ADC(2), CC2, and CCSD], time-dependent density functional theory (TD-DFT), and BSE/evGW levels of theory. First, we compute the evolution of the dipole moments of the two lowest singlet excited states of 4-(dimethylamino)benzonitrile (DMABN) upon twisting of the amino group. Next, we use a set of 25 dyes having ES characters ranging from locally excited to charge transfer to determine both μES and Δμ. For DMABN our results show that BSE/evGW provides Δμ values closer to the CCSD reference and more consistent trends than TD-DFT. Moreover, a statistical analysis of both Δμ and μES for the set of 25 dyes shows that the BSE/evGW accuracy is comparable or sometimes slightly better than that of TD-M06-2X and TD-CAM-B3LYP, BSE/evGW outperforming TD-DFT in challenging cases (zwitterionic and cyanine transitions). Finally, the starting point dependency of BSE/evGW seems to be larger for Δμ, ES dipoles, and oscillator strengths than for transition energies.

Many-body GW calculations with very large scale polarizable environments made affordable: A fully ab initio QM/QM approach

We present a many-body GW formalism for quantum subsystems embedded in discrete polarizable environments containing up to several hundred thousand atoms described at a fully ab initio random phase approximation level. Our approach is based on a fragment approximation in the construction of the Green's function and independent-electron susceptibilities. Further, the environing fragments susceptibility matrices are reduced to a minimal but accurate representation preserving low order polarizability tensors through a constrained minimization scheme. This approach dramatically reduces the cost associated with inverting the Dyson equation for the screened Coulomb potential W, while preserving the description of short to long-range screening effects. The efficiency and accuracy of the present scheme is exemplified in the paradigmatic cases of fullerene bulk, surface, subsurface, and slabs with varying number of layers.

NWChem: Recent and Ongoing Developments

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

koopmans: An Open-Source Package for Accurately and Efficiently Predicting Spectral Properties with Koopmans Functionals

Over the past decade we have developed Koopmans functionals, a computationally efficient approach for predicting spectral properties with an orbital-density-dependent functional framework. These functionals impose a generalized piecewise linearity condition to the entire electronic manifold, ensuring that orbital energies match the corresponding electron removal/addition energy differences (in contrast to semilocal DFT, where a mismatch between the two lies at the heart of the band gap problem and, more generally, the unreliability of Kohn-Sham orbital energies). This strategy has proven to be very powerful, yielding molecular orbital energies and solid-state band structures with comparable accuracy to many-body perturbation theory but at greatly reduced computational cost while preserving a functional formulation. This paper reviews the theory of Koopmans functionals, discusses the algorithms necessary for their implementation, and introduces koopmans, an open-source package that contains all of the code and workflows needed to perform Koopmans functional calculations and obtain reliable spectral properties of molecules and materials.

Linear Scaling Calculations of Excitation Energies with Active-Space Particle-Particle Random-Phase Approximation

We developed an efficient active-space particle-particle random-phase approximation (ppRPA) approach to calculate accurate charge-neutral excitation energies of molecular systems. The active-space ppRPA approach constrains both indexes in particle and hole pairs in the ppRPA matrix, which only selects frontier orbitals with dominant contributions to low-lying excitation energies. It employs the truncation in both orbital indexes in the particle-particle and the hole-hole spaces. The resulting matrix, whose eigenvalues are excitation energies, has a dimension that is independent of the size of the systems. The computational effort for the excitation energy calculation, therefore, scales linearly with system size and is negligible compared with the ground-state calculation of the (N - 2)-electron system, where N is the electron number of the molecule. With the active space consisting of 30 occupied and 30 virtual orbitals, the active-space ppRPA approach predicts the excitation energies of valence, charge-transfer, Rydberg, double, and diradical excitations with the mean absolute errors (MAEs) smaller than 0.03 eV compared with the full-space ppRPA results. As a side product, we also applied the active-space ppRPA approach in the renormalized singles (RS) T-matrix approach. Combining the non-interacting pair approximation that approximates the contribution to the self-energy outside the active space, the active-space GRSTRS@PBE approach predicts accurate absolute and relative core-level binding energies with the MAEs around 1.58 and 0.3 eV, respectively. The developed linear scaling calculation of excitation energies is promising for applications to large and complex systems.

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.

Complex Oxides for Brain-Inspired Computing: A Review

The fields of brain-inspired computing, robotics, and, more broadly, artificial intelligence (AI) seek to implement knowledge gleaned from the natural world into human-designed electronics and machines. In this review, the opportunities presented by complex oxides, a class of electronic ceramic materials whose properties can be elegantly tuned by doping, electron interactions, and a variety of external stimuli near room temperature, are discussed. The review begins with a discussion of natural intelligence at the elementary level in the nervous system, followed by collective intelligence and learning at the animal colony level mediated by social interactions. An important aspect highlighted is the vast spatial and temporal scales involved in learning and memory. The focus then turns to collective phenomena, such as metal-to-insulator transitions (MITs), ferroelectricity, and related examples, to highlight recent demonstrations of artificial neurons, synapses, and circuits and their learning. First-principles theoretical treatments of the electronic structure, and in situ synchrotron spectroscopy of operating devices are then discussed. The implementation of the experimental characteristics into neural networks and algorithm design is then revewed. Finally, outstanding materials challenges that require a microscopic understanding of the physical mechanisms, which will be essential for advancing the frontiers of neuromorphic computing, are highlighted.

Accelerating Core-Level GW Calculations by Combining the Contour Deformation Approach with the Analytic Continuation of W

In recent years, the GW method has emerged as a reliable tool for computing core-level binding energies. The contour deformation (CD) technique has been established as an efficient, scalable, and numerically stable approach to compute the GW self-energy for deep core excitations. However, core-level GW calculations with CD face the challenge of higher scaling with respect to system size N compared to the conventional quartic scaling in valence-state algorithms. In this work, we present the CD-WAC method [CD with W analytic continuation (AC)], which reduces the scaling of CD applied to the inner shells from O(N5) to O(N4) by employing an AC of the screened Coulomb interaction W. Our proposed method retains the numerical accuracy of CD for the computationally challenging deep core case, yielding mean absolute errors <5 meV for well-established benchmark sets, such as CORE65, for single-shot GW calculations. More extensive testing for different GW flavors proves the reliability of the method. We have confirmed the theoretical scaling by performing scaling experiments on large acene chains and amorphous carbon clusters, achieving speedups of up to 10× for structures of only 116 atoms. This improvement in computational efficiency paves the way for more accurate and efficient core-level GW calculations on larger and more complex systems.

Stochastic Real-Time Second-Order Green's Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green's function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(Ne3), where Ne is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

Robust relativistic many-body Green's function based approaches for assessing core ionized and excited states

A two-component contour deformation (CD) based GW method that employs frequency sampling to drastically reduce the computational effort when assessing quasiparticle states far away from the Fermi level is outlined. Compared to the canonical CD-GW method, computational scaling is reduced by an order of magnitude without sacrificing accuracy. This allows for an efficient calculation of core ionization energies. The improved computational efficiency is used to provide benchmarks for core ionized states, comparing the performance of 15 density functional approximations as Kohn-Sham starting points for GW calculations on a set of 65 core ionization energies of 32 small molecules. Contrary to valence states, GW calculations on core states prefer functionals with only a moderate amount of Hartree-Fock exchange. Moreover, modern ab initio local hybrid functionals are also shown to provide excellent generalized Kohn-Sham references for core GW calculations. Furthermore, the core-valence separated Bethe-Salpeter equation (CVS-BSE) is outlined. CVS-BSE is a convenient tool to probe core excited states. The latter is tested on a set of 40 core excitations of eight small inorganic molecules. Results from the CVS-BSE method for excitation energies and the corresponding absorption cross sections are found to be in excellent agreement with those of reference damped response BSE calculations.

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.

Exploring Bethe-Salpeter Excited-State Dipoles: The Challenging Case of Increasingly Long Push-Pull Oligomers

The change of molecular dipole moment induced by photon absorption is key to interpret the measured optical spectra. Except for compact molecules, time-dependent density functional theory (TD-DFT) remains the only theory allowing to quickly predict excited-state dipoles (μ^ES), albeit with a strong dependency on the selected exchange-correlation functional. This Letter presents the first assessment of the performances of the many-body Green's function Bethe-Salpeter equation (BSE) formalism for the evaluation of the μ^ES. We explore increasingly long push-pull oligomers as they present an excited-state nature evolving with system size. This work shows that BSE's μ^ES do present the same evolution with oligomeric length as their CC2 and CCSD counterparts, with a dependency on the starting exchange-correlation functional that is strongly decreased as compared to TD-DFT. This Letter demonstrates that BSE is a valuable alternative to TD-DFT for properties related to the excited-state density and not only for transition energies and oscillator strengths.

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.

Exact relationships between the GW approximation and equation-of-motion coupled-cluster theories through the quasi-boson formalism

We describe the relationship between the GW approximation and various equation-of-motion (EOM) coupled-cluster (CC) theories. We demonstrate the exact equivalence of the G₀W₀ approximation and the propagator theory for an electron-boson problem in a particular excitation basis. From there, we establish equivalence within the quasi-boson picture to the IP+EA-EOM unitary CC propagator. We analyze the incomplete description of screening provided by the standard similarity-transformed IP+EA-EOM-CC and the recently introduced G₀W₀ Tamm-Dancoff approximation. We further consider the approximate decoupling of IP and EA sectors in EOM-CC treatments and devise the analogous particle-hole decoupling approach for the G₀W₀ approximation. Finally, we numerically demonstrate the exact relationships and magnitude of the approximations in the calculations of a set of molecular ionization potentials and electron affinities.

Excited state potential energy surfaces of N-phenylpyrrole upon twisting: reference values and comparison between BSE/GW and TD-DFT

The puzzling case of the mixing between the charge transfer (CT) and local excited (LE) characters upon twisting of the geometry of N-phenylpyrrole (N-PP) is investigated considering the six low-lying singlet excited states (ES). The theoretical calculations of the potential energy surfaces (PES) have been performed for these states using a Coupled Cluster method accounting for the impact of the contributions from the triples, many-body Green's function GW and Bethe-Salpeter equation (BSE) formalisms, as well as Time-Dependent Density Functional Theory (TD-DFT) using various exchange-correlation functionals. Our findings confirm that the BSE formalism is more reliable than TD-DFT for close-lying ES with mixed CT/LE nature. More specifically, BSE/GW yields a more accurate evolution of the excited state PES than TD-DFT when compared to the reference coupled cluster values. BSE/GW PES curves also show negligible exchange-correlation functional starting point dependency in sharp contrast with their TD-DFT counterparts.

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.

Inverse Design of Tetracene Polymorphs with Enhanced Singlet Fission Performance by Property-Based Genetic Algorithm Optimization

The efficiency of solar cells may be improved by using singlet fission (SF), in which one singlet exciton splits into two triplet excitons. SF occurs in molecular crystals. A molecule may crystallize in more than one form, a phenomenon known as polymorphism. Crystal structure may affect SF performance. In the common form of tetracene, SF is experimentally known to be slightly endoergic. A second, metastable polymorph of tetracene has been found to exhibit better SF performance. Here, we conduct inverse design of the crystal packing of tetracene using a genetic algorithm (GA) with a fitness function tailored to simultaneously optimize the SF rate and the lattice energy. The property-based GA successfully generates more structures predicted to have higher SF rates and provides insight into packing motifs associated with improved SF performance. We find a putative polymorph predicted to have superior SF performance to the two forms of tetracene, whose structures have been determined experimentally. The putative structure has a lattice energy within 1.5 kJ/mol of the most stable common form of tetracene.

Interpolative Separable Density Fitting for Accelerating Two-Electron Integrals: A Theoretical Perspective

Low-rank approximations have long been considered an efficient way to accelerate electronic structure calculations associated with the evaluation of electron repulsion integrals (ERIs). As an accurate and efficient algorithm for compressing the ERI tensor, the interpolative separable density fitting (ISDF) decomposition has recently attracted great attention in this context. In this perspective, we introduce the ISDF decomposition from the theoretical aspects and technique details. The ISDF decomposition can construct a fully separable low-rank approximation (tensor hypercontraction factorization) of ERIs in real space with a cubic cost, offering great flexibility for accelerating high-scaling electronic structure calculations. We review the typical applications of ISDF in hybrid functionals, time-dependent density functional theory, and GW approximation. Finally, we discuss the promising directions for future development of ISDF.