The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy

Extrapolating Unconverged GW Energies up to the Complete Basis Set Limit with Linear Regression

The GW approximation to the electronic self-energy is now a well-recognized approach to obtain the electron quasiparticle energies of molecules and, in particular, their ionization potential and electron affinity. Though much faster than the corresponding wavefunction methods, the GW energies are still affected by slow convergence with respect to the basis completeness. This limitation hinders a wider application of the GW approach. Here, we show that we can reach the complete basis set limit for the cumbersome GW calculations solely based on fast preliminary calculations with an unconverged basis set. We introduce a linear model that correlates the molecular orbital characteristics and the basis convergence error for a large database of approximately 600 states in 104 organic molecules that contain H, C, O, N, F, P, S, and Cl. The model employs molecular-orbital-based non-linear descriptors that encode efficiently the chemical space offering outstanding transferability. Using a low number of descriptors (17) the performance of this extrapolation procedure is superior to that of the earlier more physically motivated approaches. The predictive power of the method is finally demonstrated for a selection of large acceptor molecules.

Computation of Molecular Ionization Energies Using an Ensemble Density Functional Theory Method

Computation of the ionization energies and of the respective Dyson orbitals based on the use of the extended Koopmans theorem (EKT) is implemented in connection with an ensemble density functional theory (eDFT) method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method. The new methodology enables fast computation of the ionization energies and evaluation of the respective Dyson orbitals, the square norms of which are related with the ionization probabilities, in the ground and excited electronic states of molecules. As the application of EKT recycles the intermediate quantities from the SSR analytical energy gradient, evaluation of the ionization energies and probabilities can be carried out on-the-fly during the nonadiabatic molecular dynamics simulations. This opens up a perspective for fast theoretical simulation of the time-resolved photoelectron spectroscopy observations. In the present work, the new methodology is tested in the computation of the ionization energies and Dyson orbitals of several molecules in the ground and excited electronic states, including strongly correlated species, such as the ozone molecule, dissociating chemical bonds, and conical intersections.

CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

Excited-state electronic structure of molecules using many-body Green's functions: Quasiparticles and electron-hole excitations with VOTCA-XTP

We present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green's function theory in the GW approximation with the Bethe-Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA-XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW-BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron-hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW-BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene-fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron-hole excited states of various characters on an equal footing.

Accurate Absolute and Relative Core-Level Binding Energies from GW

We present an accurate approach to compute X-ray photoelectron spectra based on the GW Green's function method that overcomes the shortcomings of common density functional theory approaches. GW has become a popular tool to compute valence excitations for a wide range of materials. However, core-level spectroscopy is thus far almost uncharted in GW. We show that single-shot perturbation calculations in the G0W0 approximation, which are routinely used for valence states, cannot be applied for core levels and suffer from an extreme, erroneous transfer of spectral weight to the satellite spectrum. The correct behavior can be restored by partial self-consistent GW schemes or by using hybrid functionals with almost 50% of exact exchange as a starting point for G0W0. We also include relativistic corrections and present a benchmark study for 65 molecular 1s excitations. Our absolute and relative GW core-level binding energies agree within 0.3 and 0.2 eV with experiment, respectively.

Atomic structures and orbital energies of 61,489 crystal-forming organic molecules

Data science and machine learning in materials science require large datasets of technologically relevant molecules or materials. Currently, publicly available molecular datasets with realistic molecular geometries and spectral properties are rare. We here supply a diverse benchmark spectroscopy dataset of 61,489 molecules extracted from organic crystals in the Cambridge Structural Database (CSD), denoted OE62. Molecular equilibrium geometries are reported at the Perdew-Burke-Ernzerhof (PBE) level of density functional theory (DFT) including van der Waals corrections for all 62 k molecules. For these geometries, OE62 supplies total energies and orbital eigenvalues at the PBE and the PBE hybrid (PBE0) functional level of DFT for all 62 k molecules in vacuum as well as at the PBE0 level for a subset of 30,876 molecules in (implicit) water. For 5,239 molecules in vacuum, the dataset provides quasiparticle energies computed with many-body perturbation theory in the G0W0 approximation with a PBE0 starting point (denoted GW5000 in analogy to the GW100 benchmark set (M. van Setten et al. J. Chem. Theory Comput. 12, 5076 (2016))).

From flat to tilted: gradual interfaces in organic thin film growth

We investigate domain formation and local morphology of thin films of α-sexithiophene (α-6T) on Au(100) beyond monolayer coverage by combining high resolution scanning tunneling microscopy (STM) experiments with electronic structure theory calculations and computational structure search. We report a layerwise growth of highly-ordered enantiopure domains. For the second and third layer, we show that the molecular orbitals of individual α-6T molecules can be well resolved by STM, providing access to detailed information on the molecular orientation. We find that already in the second layer the molecules abandon the flat adsorption structure of the monolayer and adopt a tilted conformation. Although the observed tilted arrangement resembles the orientation of α-6T in the bulk, the observed morphology does not yet correspond to a well-defined surface of the α-6T bulk structure. A similar behavior is found for the third layer indicating a growth mechanism where the bulk structure is gradually adopted over several layers.

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.