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.

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

AB-G$_0$W$_0$: A practical G$_0$W$_0$ method without frequency integration based on an auxiliary boson expansion

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Precision benchmarks for solids: G0W0 calculations with different basis sets

Low Rank Approximation in $G_0W_0$ Approximation

Low rank approximation in G 0 W 0 calculations

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

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

Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

Solving multi-pole challenges in the GW100 benchmark enables precise low-scaling GW calculations

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

Fully dynamic G3W2 self-energy for finite systems: Formulas and benchmark

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

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

Efficient Full-Frequency GW Calculations Using a Lanczos Method

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

Generalized Discrete Truncated Wigner Approximation for Nonadiabtic Quantum-Classical Dynamics

Accelerating Nonequilibrium Green functions simulations: the G1-G2 scheme and beyond

The $GW$ Approximation: A Quantum Chemistry Perspective

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW