Abstract: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.

$G_0W_0$ Ionization Potentials of First-Row Transition Metal Aqua Ions

Ab Initio Investigation of the Lower Energy Candidate Structures for (h2o)(5)+ Water Cluster

Ab Initio Investigation of Possible Candidate Structures and Properties of Water Cluster (H2O)7+ Via Particle Swarm Optimization Method

Ab Initio Investigation Of Structure, Stability, Thermal Behavior, Bonding, And Infrared Spectra Of Ionized Water Cluster (H2o)(6)(+) (Vol 145, 154307, 2016)

Exploiting a derivative discontinuity estimate for accurate $G_0W_0$ ionization potentials and electron affinities

Photoemission spectra and density functional theory calculations of 3d transition metal-aqua complexes (Ti-Cu) in aqueous solution

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

Two-component $GW$ calculations: Cubic scaling implementation and comparison of vertex corrected and partially self-consistent $GW$ variants

Comparing self-consistent GW and vertex corrected G0W0 (G0W0Γ) accuracy for molecular ionization potentials

Precision benchmarks for solids: G0W0 calculations with different basis sets

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

Relativistic fully self-consistent $GW$ for molecules: Total energies and ionization potentials

Vertical Ionization Energies, Generalized Kohn-Sham Orbital Energies, and the Curious Case of the Copper Oxide Anions

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics.

All-electron first-principles GWΓ simulations for accurately predicting core-electron binding energies considering first-order three-point vertex corrections

Fully selfconsistent GW calculations for molecules

Electron affinities of water clusters from density-functional and many-body-perturbation theory

q-AQUA: a many-body CCSD(T) water potential, including 4-body interactions, demonstrates the quantum nature of water from clusters to the liquid phase

Low-order Scaling $G_0W_0$ by Pair Atomic Density Fitting

Ab initio studies of ionization potentials of hydrated hydroxide and hydronium