Low-scaling \textit{GW} calculations of quasi-particle energies for extended systems within the numerical atomic orbital framework

Abstract

The many-body perturbation theory within the $GW$ approximation is a widely used method for describing the electronic band structures in real materials. Its application to large-scale systems is, however, impeded by its high computational cost. The rate-limiting steps in a typical $GW$ implementation are the evaluation of the polarization function under the random phase approximation (RPA) and the evaluation of the $GW$ self-energy, both of which have a canonical $O(N^4)$ scaling with $N$ being the system size. The conventional space-time algorithm within the plane-wave basis sets reduces the scaling from $O(N^4)$ to $O(N^3)$, albeit with a large prefactor and increased memory cost. Here, we present a space-time algorithm within the numerical atomic orbital (NAO) basis-set framework, for which the evaluation of the polarization function and self-energy is formally reduced to $O(N^2)$ or better with respect to system size. This is achieved by computing these quantities in real space, where low-scaling algorithms can be formulated by leveraging the localized resolution of identity (LRI) technique. The resulting NAO-based, LRI-enhanced space-time $GW$ algorithm has been implemented in the LibRPA library interfaced with the FHI-aims code package. Benchmark calculations for crystalline solids show that the low-scaling implementation yields quasi-particle energies in close agreement with the conventional $O(N^4)$ k-space formalism previously implemented in FHI-aims. For the systems studied here, the observed overall scaling is substantially reduced relative to the canonical approach, and the low-scaling implementation becomes advantageous already for systems containing fewer than 100 atoms.

59 pages, 11 figures

References in corpus (19)

• The GW compendium: A practical guide to theoretical photoemission spectroscopy
• Predictive GW calculations using plane waves and pseudopotentials
• Cubic scaling $GW$: towards fast quasiparticle calculations
• Towards GW Calculations on Thousands of Atoms
• Hybrid-functional calculations with plane-wave basis sets: The effect of the singularity correction on total energies, energy eigenvalues, and defect energy levels
• Cubic-scaling all-electron GW calculations with a separable density-fitting space-time approach
• Low-order Scaling $G_0{W}_0$ by Pair Atomic Density Fitting
• All-electron Gaussian-based $G_0{W}_0$ for Valence and Core Excitation Energies of Periodic Systems
• Low-scaling $GW$ with benchmark accuracy and application to phosphorene nanosheets
• All-electron periodic $G_0{W}_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks
• GW100: A Slater Type Orbital Perspective
• Efficient All-electron Hybrid Density Functionals for Atomistic Simulations Beyond 10,000 Atoms
• Towards Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy
• Spin-orbit torque in single-molecule junctions from ab initio
• Low-scaling GW algorithm applied to twisted transition-metal dichalcogenide heterobilayers
• Efficient Evaluation of Exact Exchange for Periodic Systems via Concentric Atomic Density Fitting
• Roadmap on Advancements of the FHI-aims Software Package
• Finite-size Effects in periodic EOM-CCSD for Ionization Energies and Electron Affinities: Convergence Rate and Extrapolation to the Thermodynamic Limit
• Benchmarking the accuracy of the separable resolution of the identity approach for correlated methods in the numeric atom-centered orbitals framework

Cited by in corpus (1)

• All-electron Quasiparticle Self-consistent GW for Molecules and Periodic Systems within the Numerical Atomic Orbital Framework