Low-scaling GW algorithm applied to twisted transition-metal dichalcogenide heterobilayers

Abstract

The $GW$ method is widely used for calculating the electronic band structure of materials. The high computational cost of $GW$ algorithms prohibits their application to many systems of interest. We present a periodic, low-scaling and highly efficient $GW$ algorithm that benefits from the locality of the Gaussian basis and the polarizability. The algorithm enables $G_0{W}_0$ calculations on a MoSe${}_2$/WS${}_2$ bilayer with 984 atoms per unit cell, in 42 hours using 1536 cores. This is four orders of magnitude faster than a plane-wave $G_0{W}_0$ algorithm, allowing for unprecedented computational studies of electronic excitations at the nanoscale.

References in corpus (10)

• Resonantly hybridised excitons in moiré superlattices in van der Waals heterostructures
• Many-body perturbation theory calculations using the yambo code
• The GW compendium: A practical guide to theoretical photoemission spectroscopy
• The Bethe-Salpeter Equation Formalism: From Physics to Chemistry
• Cubic scaling $GW$: towards fast quasiparticle calculations
• Towards GW Calculations on Thousands of Atoms
• 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
• Fast evaluation of solid harmonic Gaussian integrals for local resolution-of-the-identity methods and range-separated hybrid functionals
• Periodic Coupled-Cluster Green's Function for Photoemission Spectra of Realistic Solids

Cited by in corpus (4)

• Linear and Nonlinear Optical Properties of Molecules from Real-Time Propagation Based on the Bethe-Salpeter Equation
• Parquet theory for molecular systems: Formalism and static kernel parquet approximation
• Low-scaling \textit{GW} calculations of quasi-particle energies for extended systems within the numerical atomic orbital framework
• Efficient $GW$ band structure calculations using Gaussian basis functions and application to atomically thin transition-metal dichalcogenides