Combined sub-sampling and analytical integration for efficient large-scale $GW$ calculations for 2D systems

Abstract

Accurate and efficient predictions of the quasiparticle properties of complex materials remain a major challenge due to the convergence issue and the unfavorable scaling of the computational cost with respect to the system size. Quasiparticle $GW$ calculations for two dimensional (2D) materials are especially difficult. The unusual analytical behaviors of the dielectric screening and the electron self-energy of 2D materials make the conventional Brillouin zone (BZ) integration approach rather inefficient and require an extremely dense $k$-grid to properly converge the calculated quasiparticle energies. In this work, we present a combined non-uniform sub-sampling and analytical integration method that can drastically improve the efficiency of the BZ integration in 2D $GW$ calculations. Our work is distinguished from previous work in that, instead of focusing on the intricate dielectric matrix or the screened Coulomb interaction matrix, we exploit the analytical behavior of various terms of the convolved self-energy $\Sigma (\mathbf{q})$ in the small $\mathbf{q}$ limit. This method, when combined with another accelerated $GW$ method that we developed recently, can drastically speed-up (by over three orders of magnitude) $GW$ calculations for 2D materials. Our method allows fully converged $GW$ calculations for complex 2D systems at a fraction of computational cost, facilitating future high throughput screening of the quasiparticle properties of 2D semiconductors for various applications. To demonstrate the capability and performance of our new method, we have carried out fully converged $GW$ calculations for monolayer C${}_2$N, a recently discovered 2D material with a large unit cell, and investigate its quasiparticle band structure in detail.

References in corpus (16)

• Optical Spectrum of MoS${}_2$: Many-body Effects and Diversity of Exciton States
• BerkeleyGW: A Massively Parallel Computer Package for the Calculation of the Quasiparticle and Optical Properties of Materials and Nanostructures
• Quasiparticle band structures and optical properties of strained monolayer MoS2 and WS2
• Dielectric screening in two-dimensional insulators: Implications for excitonic and impurity states in graphane
• An exact Coulomb cutoff technique for supercell calculations
• Effect of spin-orbit interaction on the excitonic effects in single-layer, double-layer, and bulk MoS2
• Screening and Many-Body Effects in Two-Dimensional Crystals: Monolayer MoS${}_2$
• Truncation of Periodic Image Interactions for Confined Systems
• X-ray Absorption Near-Edge Structure calculations with pseudopotentials. Application to K-edge in diamond and alpha-quartz
• How dielectric screening in two-dimensional crystals affects the convergence of excited-state calculations: Monolayer MoS${}_2$
• Non-uniform sampling schemes of the Brillouin zone for many-electron perturbation-theory calculations in reduced dimensionality
• Efficient many-body calculations of 2D materials using exact limits for the screened potential: Band gaps of MoS${}_2$, hBN, and phosphorene
• Towards Fully Converged GW Calculations for Large Systems
• Screening in 2D: GW calculations for surfaces and thin films using the repeated-slab approach
• A many-body GW+BSE investigation of electronic and optical properties of C2N
• Tunable Half-metallic Magnetism in Atom-thin Holey Two-dimensional C${}_2$N Monolayer