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.

39 pages, 9 figures

References in corpus (12)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• NWChem: Past, Present, and Future
• 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
• Quasi-Particle Self-Consistent $GW$ for Molecules
• Towards GW Calculations on Thousands of Atoms
• 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

Cited by in corpus (1)

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