Critical assessment of $G_0W_0$ calculations for 2D materials: the example of monolayer MoS$_2$

Ronaldo Rodrigues Pela,Cecilia Vona,Sven Lubeck,Ben Alex,Ignacio Gonzalez Oliva,Claudia Draxl
DOI: https://doi.org/10.1038/s41524-024-01253-2
2024-04-23
Abstract:Two-dimensional (2D) materials combine many fascinating properties that make them more interesting than their three-dimensional counterparts for a variety of applications. For example, 2D materials exhibit stronger electron-phonon and electron-hole interactions, and their energy gaps and effective carrier masses can be easily tuned. Surprisingly, published band gaps of several 2D materials obtained with the $GW$ approach, the state-of-the-art in electronic-structure calculations, are quite scattered. The details of these calculations, such as the underlying geometry, the starting point, the inclusion of spin-orbit coupling, and the treatment of the Coulomb potential can critically determine how accurate the results are. Taking monolayer MoS$_2$ as a representative material, we employ the linearized augmented planewave + local orbital method to systematically investigate how all these aspects affect the quality of $G_0W_0$ calculations, and also provide a summary of literature data. We conclude that the best overall agreement with experiments and coupled-cluster calculations is found for $G_0W_0$ results with HSE06 as a starting point including spin-orbit coupling, a truncated Coulomb potential, and an analytical treatment of the singularity at $q=0$.
Materials Science,Other Condensed Matter,Chemical Physics
What problem does this paper attempt to address?
This paper attempts to address the inconsistency issues in calculating band gaps of two-dimensional (2D) materials using the GW method. Specifically: 1. **Research Background**: 2D materials have garnered significant attention due to their unique electronic, thermal, and optical properties. For instance, their band structures can be easily tuned by external electric fields or strain. However, the band gap values of different 2D materials calculated using the GW method vary greatly, posing challenges for theoretical research. 2. **Core Issue**: When using the GW method (especially the G0W0 approximation) to calculate the band gap of 2D materials such as molybdenum disulfide (MoS2), factors like different geometric structures, initial point selection, treatment of spin-orbit coupling, and the way Coulomb potential is truncated can lead to significant discrepancies in the results. 3. **Research Objective**: The paper systematically investigates how the aforementioned factors affect the quality of G0W0 calculations using monolayer MoS2 as an example and provides a summary of literature data. The study finds that using HSE06 as the initial point, including spin-orbit coupling, truncating the Coulomb potential, and analytically treating the singularity at q=0 can yield results that best match experimental and coupled-cluster calculations. Through these analyses, the paper aims to provide more reliable methods and benchmark datasets for the electronic structure calculations of 2D materials.