Scalable Molecular GW Calculations: Valence and Core Spectra

Abstract

We present a scalable implementation of the $GW$ approximation using Gaussian atomic orbitals to study the valence and core ionization spectroscopies of molecules. The implementation of the standard spectral decomposition approach to the screened Coulomb interaction, as well as a contour deformation method are described. We have implemented both of these approaches using the robust variational fitting approximation to the four-center electron repulsion integrals. We have utilized the MINRES solver with the contour deformation approach to reduce the computational scaling by one order of magnitude. A complex heuristic in the quasiparticle equation solver further allows a speed-up of the computation of core and semi-core ionization energies. Benchmark tests using the GW100 and CORE65 datasets and the carbon 1{\it s} binding energy of the well-studied ethyl trifluoroacetate, or ESCA molecule, were performed to validate the accuracy of our implementation. We also demonstrate and discuss the parallel performance and computational scaling of our implementation using a range of water clusters of increasing size.

arXiv comment39 pages, 9 figures

References in corpus (31)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• A real-space grid implementation of the Projector Augmented Wave method
• The GW method
• Recent developments in the PySCF program package
• Yambo: an \textit{ab initio} tool for excited state calculations
• BerkeleyGW: A Massively Parallel Computer Package for the Calculation of the Quasiparticle and Optical Properties of Materials and Nanostructures
• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• NWChem: Past, Present, and Future
• All-electron self-consistent GW approximation: Application to Si, MnO, and NiO
• Many-body perturbation theory calculations using the yambo code
• The GW compendium: A practical guide to theoretical photoemission spectroscopy
• Efficient implementation of the GW approximation within the all-electron FLAPW method
• Quasiparticle GW calculations for solids, molecules and 2D materials
• Unified description of ground and excited states of finite systems: the self-consistent GW approach
• GW quasi-particle spectra from occupied states only
• Levels of self-consistency in the GW approximation
• Self-consistent GW: All-electron implementation with localized basis functions
• Quasi-Particle Self-Consistent $GW$ for Molecules
• Towards GW Calculations on Thousands of Atoms
• Questaal: a package of electronic structure methods based on the linear muffin-tin orbital technique
• Applicability of the strongly constrained and appropriately normed density functional to transition metal magnetism
• Accurate absolute and relative core-level binding energies from $GW$
• Low-order Scaling $G_0{W}_0$ by Pair Atomic Density Fitting
• Meta-GGA Performance in Solids at Almost GGA Cost
• Density Functional Methods for the Magnetism of Transition Metals: SCAN in Relation to Other Functionals
• Low-scaling $GW$ with benchmark accuracy and application to phosphorene nanosheets
• Analysis of over-magnetization of elemental transition metal solids from the SCAN Density Functional
• Optimizations of the Eigensolvers in the ELPA Library
• A mathematical analysis of the GW0 method for computing electronic excited energies of molecules

Cited by in corpus (1)

• On the basis set selection for molecular core-level $GW$ calculations