All-electron Quasiparticle Self-consistent GW for Molecules and Periodic Systems within the Numerical Atomic Orbital Framework
Pith reviewed 2026-05-22 10:45 UTC · model grok-4.3
The pith
NAO-based quasiparticle self-consistent GW produces ionization potentials and band gaps consistent with reference implementations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We report an all-electron implementation of the quasiparticle self-consistent GW (QSGW) method for molecular and periodic systems within the framework of numerical atomic orbitals (NAOs). The implementation rests on the space-time formalism combined with the localized resolution-of-identity approximation. Analytical continuation of the self-energy matrix together with the Mode B QSGW scheme yields stable self-consistent quasiparticle energy spectra. Benchmark calculations on molecules and crystalline solids demonstrate that the resulting molecular ionization potentials and quasiparticle band gaps are consistent with reference values from established implementations.
What carries the argument
The NAO-based QSGW scheme that employs space-time formalism and localized resolution-of-identity to treat two-electron quantities while performing self-consistent quasiparticle updates.
If this is right
- Molecular ionization potentials become accessible at the QSGW level with all-electron NAO accuracy.
- Quasiparticle band gaps of semiconductors and wide-gap insulators match reference values from plane-wave or other established codes.
- The low-scaling NAO infrastructure already used for G0W0 can now be applied to full QSGW calculations.
- Large systems that were previously out of reach for self-consistent GW become computationally tractable.
Where Pith is reading between the lines
- The same NAO basis and resolution-of-identity machinery could be reused to embed QSGW corrections inside larger-scale hybrid or embedding schemes.
- Extension to systems with hundreds of atoms would test whether the observed consistency with reference methods persists when system size increases.
- Because the method is all-electron, it offers a direct route to core-level quasiparticle energies that are often approximated in pseudopotential codes.
Load-bearing premise
Analytical continuation of the self-energy matrix combined with the Mode B QSGW scheme produces stable self-consistent quasiparticle energy spectra.
What would settle it
A direct comparison on a small molecule or semiconductor where the self-consistent quasiparticle energies obtained with analytical continuation and Mode B diverge from or show larger errors than results from contour-deformation or other stable continuation methods.
Figures
read the original abstract
We report an all-electron implementation of the quasiparticle self-consistent GW (QSGW) method for molecular and periodic systems within the framework of numerical atomic orbitals (NAOs), as implemented in the LibRPA software package. Our implementation is based on the space-time formalism, combined with the localized resolution-of-identity approximation to treat two-electron quantities. We found that analytical continuation of the self-energy matrix, in combination with the ``Mode B" QSGW scheme, can yield stable self-consistent quasiparticle energy spectra. Systematic benchmark calculations on molecules and crystalline solids (including typical semiconductors and wide-gap insulators) demonstrate that our NAO-based QSGW scheme yields molecular ionization potentials and quasiparticle band gaps for periodic solids that are consistent with reference results from established implementations. Our work opens the way for large-scale QSGW calculations, taking advantage of the NAO-based low-scaling algorithm previously developed for the G0W0 method.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports an all-electron implementation of quasiparticle self-consistent GW (QSGW) for molecules and periodic systems using numerical atomic orbitals (NAOs) within the LibRPA package. It employs the space-time formalism together with a localized resolution-of-identity approximation for two-electron integrals. The authors state that analytical continuation of the self-energy matrix combined with the Mode B QSGW scheme produces stable self-consistent quasiparticle spectra, and they present systematic benchmarks on molecular ionization potentials and quasiparticle band gaps of semiconductors and wide-gap insulators that are consistent with results from established reference codes. The work emphasizes the potential for large-scale calculations via the low-scaling NAO algorithm previously developed for G0W0.
Significance. If the implementation and benchmarks are robust, the approach would enable efficient all-electron QSGW calculations on larger molecular and periodic systems, extending the reach of parameter-free quasiparticle methods in materials science. The combination of NAO basis sets with space-time techniques and localized RI is a concrete technical contribution that could be adopted by other codes.
major comments (1)
- §4 (Benchmark calculations) and the discussion of self-consistency: the assertion that analytical continuation of the self-energy matrix with Mode B yields stable spectra is central to the reliability of the reported quasiparticle energies, yet the manuscript provides no quantitative metrics such as number of iterations to convergence, sensitivity of final eigenvalues to the fitting range or number of poles in the continuation, or direct comparison against contour-deformation or real-frequency evaluations for the same systems. For wide-gap insulators, where the self-energy frequency structure is pronounced, small continuation artifacts could propagate through the self-consistent loop and affect the claimed consistency with reference results.
minor comments (2)
- The abstract and introduction would benefit from explicit numerical error bars or mean absolute deviations for the benchmarked ionization potentials and band gaps rather than the qualitative statement of 'consistency'.
- Notation for the localized RI approximation and the precise definition of Mode B should be cross-referenced to the earlier G0W0 paper to avoid ambiguity for readers.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comment on the validation of our self-consistent procedure. We address the point in detail below and are happy to revise the manuscript accordingly.
read point-by-point responses
-
Referee: §4 (Benchmark calculations) and the discussion of self-consistency: the assertion that analytical continuation of the self-energy matrix with Mode B yields stable spectra is central to the reliability of the reported quasiparticle energies, yet the manuscript provides no quantitative metrics such as number of iterations to convergence, sensitivity of final eigenvalues to the fitting range or number of poles in the continuation, or direct comparison against contour-deformation or real-frequency evaluations for the same systems. For wide-gap insulators, where the self-energy frequency structure is pronounced, small continuation artifacts could propagate through the self-consistent loop and affect the claimed consistency with reference results.
Authors: We agree that additional quantitative metrics would strengthen the presentation. In the revised manuscript we will add to §4 explicit data on the number of iterations to convergence for all benchmark systems, together with tests of the sensitivity of the final quasiparticle energies to the fitting range and the number of poles retained in the analytical continuation. These checks confirm rapid convergence (typically fewer than ten iterations) and only weak dependence on the continuation parameters for the systems studied, including wide-gap insulators. Direct side-by-side comparisons with contour-deformation or real-frequency evaluations are not currently feasible within our space-time implementation; however, the close agreement of our results with multiple independent reference codes that employ different frequency treatments provides supporting evidence that continuation artifacts remain small. We will include a brief discussion of this point in the revision. revision: yes
Circularity Check
No significant circularity: implementation benchmarked against independent references
full rationale
The paper reports an all-electron NAO-based implementation of the established QSGW method using the space-time formalism and localized RI approximation. The central results consist of benchmark comparisons showing that molecular IPs and solid quasiparticle gaps are consistent with outputs from independent, established QSGW codes. No equations or procedures reduce a claimed prediction to a fitted parameter defined from the same data, nor does any load-bearing step rest on a self-citation chain whose validity is presupposed by the present work. The prior NAO-G0W0 algorithm is cited only as a computational building block; the QSGW self-consistency and stability claims are validated externally rather than by construction. The derivation chain therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The space-time formalism combined with localized resolution-of-identity accurately treats two-electron quantities in the NAO basis.
- domain assumption Analytical continuation of the self-energy matrix with Mode B QSGW yields stable quasiparticle spectra.
Reference graph
Works this paper leans on
-
[1]
Start with an initial guess for the effective potential V xc, usually taken from a preceding DFT calcula- tion
-
[2]
Construct the effective HamiltonianH QSGW using the currentV xc and solve the eigenvalue problem to obtain the QSGWorbitalsψ p(r) and eigenvalues ϵp. 4
-
[3]
Build the non-interacting Green’s functionG 0 from the QSGWorbitals and eigenvalues
-
[4]
Calculate the self-energy Σ(r,r ′, ω) using the GWapproximation based onG 0 as described in Eqs. (2)–(6)
-
[5]
Update the effective potentialV xc using the chosen QSGWmode (A or B) to obtain a new static and Hermitian approximation of the self-energy
-
[6]
Check for convergence. If the change in the quasi- particle energies is below a predefined threshold, the calculation is considered converged. Otherwise, return to step 2. The above-outlined QSGWscheme has been imple- mented in the LibRPA package [61, 63]. LibRPA was originally developed for the RPA correlation energy [63] and one-shotG 0W 0 calculations ...
-
[7]
NAO basis set representation First, the single-particle orbitals (KS orbitals in the initial iteration, or the QSGWorbitals in subsequent it- erations) are expanded in a basis of numerical atomic orbitals (NAOs): ψpk(r) = X i ci pkφik(r) (11) Here,c i pk are the expansion coefficients for thep-th state at Bloch wavevectorkin terms of the NAO basis. The Bl...
-
[8]
Localized resolution of identity technique For QSGW, we employ an auxiliary basis set to rep- resent the two-electron quantities and LRI to decompose the four-index tensors into three- and two-index ones, similar to the canonicalG 0W 0 implementation for peri- odic systems in FHI-aims [59] and our low-scaling RPA andG 0W 0 implementation in LibRPA [60, 61...
-
[9]
Space-time formalism of QSGWwith NAO and LRI The non-interacting Green’s function in real space and imaginary time is expressed in the NAO basis as: G0(r,r ′,iτ) = X ij X R,R′ φiR(r)G0 ij(R′ −R,iτ)φ jR′(r′), (15) with the matrix elements given by G0 ij(R,iτ) = 1 Nk X pk e−ik·Rci pkcj∗ pkΞpk(τ),(16) 5 where Ξpk(τ) =e −(ϵpk−ϵF)τ [θ(τ)(1−f pk)−θ(−τ)f pk].(17...
-
[10]
Basis sets For the benchmark calculations presented in this work, we adopt different basis sets for periodic solids and molecules. For periodic systems, we employ the intermediate gwsetting, which offers a practical com- promise between accuracy and computational cost for the present test set. For molecular systems, we use the really tightsetting, i.e., t...
work page 2020
-
[11]
and to other applications such as core-level binding energies [89], relativistic two-componentGW[90], mag- netic systems [91], and data-driven modeling [92]. Going beyond theGWapproximation to include higher-order correlations [93] remains an active research direction. Last but not least, our present QSGWis built upon the existing low-scalingG 0W 0 implem...
- [12]
- [13]
- [14]
-
[15]
G. Strinati, H. J. Mattausch, and W. Hanke, Physical Review B25, 2867 (1982)
work page 1982
-
[16]
M. S. Hybertsen and S. G. Louie, Physical Review B35, 5585 (1987)
work page 1987
-
[17]
R. W. Godby, M. Schl¨ uter, and L. J. Sham, Physical Review B37, 10159 (1988)
work page 1988
-
[18]
F. Aryasetiawan and O. Gunnarsson, Physical Review Letters74(1995), 10.1103/physrevlett.74.3221
-
[19]
W. G. Aulbur, L. J¨ onsson, and J. W. Wilkins, inSolid State Physics, Vol. 54 (Elsevier, 2000) pp. 1–218
work page 2000
- [20]
- [21]
- [22]
-
[23]
D. Golze, M. Dvorak, and P. Rinke, Frontiers in Chem- istry7(2019), 10.3389/fchem.2019.00377
- [24]
-
[25]
M. F. Lange and T. C. Berkelbach, J. Chem. Theory Comput.14, 4224 (2018)
work page 2018
-
[26]
J. P. Perdew, K. Burke, and M. Ernzerhof, Physical Review Letters77, 3865 (1996)
work page 1996
-
[27]
F. Bruneval and M. A. L. Marques, J. Chem. Theory Comput.9, 324 (2013)
work page 2013
-
[28]
T. K¨ orzd¨ orfer and N. Marom, Physical Review B86, 041110 (2012)
work page 2012
- [29]
-
[30]
S. E. Gant, J. B. Haber, M. R. Filip, F. Sagredo, D. Wing, G. Ohad, L. Kronik, and J. B. Neaton, Phys. Rev. Mater.6, 053802 (2022)
work page 2022
- [31]
-
[32]
W. Ku and A. G. Eguiluz, Physical Review Letters89 (2002), 10.1103/physrevlett.89.126401
-
[33]
A. Stan, N. E. Dahlen, and R. van Leeuwen, Europhysics Letters (EPL)76, 298 (2006)
work page 2006
-
[34]
C. Rostgaard, K. W. Jacobsen, and K. S. Thygesen, Phys. Rev. B81(2010), 10.1103/physrevb.81.085103
- [35]
- [36]
-
[37]
A. L. Kutepov, Phys. Rev. B94, 155101 (2016)
work page 2016
-
[38]
A. L. Kutepov, Phys. Rev. B95, 195120 (2017)
work page 2017
- [39]
- [40]
- [41]
- [42]
-
[43]
M. Wen, V. Abraham, G. Harsha, A. Shee, K. B. Whaley, and D. Zgid, Journal of Chemical Theory and Computa- tion20, 3109 (2024)
work page 2024
-
[44]
E. Maggio and G. Kresse, Journal of chemical theory and computation13, 4765 (2017)
work page 2017
-
[45]
A. L. Kutepov, S. Y. Savrasov, and G. Kotliar, Phys. Rev. B80, 041103 (2009)
work page 2009
- [46]
- [47]
-
[48]
M. Shishkin, M. Marsman, and G. Kresse, Physical re- view letters99, 246403 (2007)
work page 2007
-
[49]
S. V. Faleev, M. van Schilfgaarde, and T. Kotani, Phys. Rev. Lett.93, 126406 (2004)
work page 2004
-
[50]
M. van Schilfgaarde, T. Kotani, and S. Faleev, Physical Review Letters96, 226402 (2006)
work page 2006
- [51]
-
[52]
A. Marie and P.-F. c. Loos, Journal of Chemical Theory and Computation19, 3943 (2023)
work page 2023
-
[53]
B. Cunningham, M. Gr¨ uning, D. Pashov, and M. van Schilfgaarde, Phys. Rev. B108, 165104 (2023)
work page 2023
-
[54]
A. L. Kutepov, Phys. Rev. B105, 045124 (2022)
work page 2022
- [55]
-
[56]
H. Sakakibara, T. Kotani, M. Obata, and T. Oda, Phys- ical Review B101, 205120 (2020)
work page 2020
- [57]
-
[58]
W. R. Lambrecht, C. Bhandari, and M. Van Schilf- gaarde, Physical Review Materials1, 043802 (2017)
work page 2017
-
[59]
S. Das, J. E. Coulter, and E. Manousakis, Physical Re- view B91, 115105 (2015)
work page 2015
-
[60]
N. Salas-Illanes, D. Nabok, and C. Draxl, Physical Re- view B106, 045143 (2022)
work page 2022
-
[61]
D. Deguchi, K. Sato, H. Kino, and T. Kotani, Japanese Journal of Applied Physics55, 051201 (2016)
work page 2016
- [62]
- [63]
-
[64]
M. Feyereisen, G. Fitzgerald, and A. Komornicki, Chem- ical physics letters208, 359 (1993)
work page 1993
-
[65]
O. Vahtras, J. Alml¨ of, and M. Feyereisen, Chemical Physics Letters213, 514 (1993)
work page 1993
-
[66]
F. Weigend, M. H¨ aser, H. Patzelt, and R. Ahlrichs, Chemical physics letters294, 143 (1998)
work page 1998
-
[67]
A. F¨ orster and L. Visscher, J. Chem. Theory Comput. 16, 7381 (2020)
work page 2020
-
[68]
M. J. van Setten, F. Caruso, S. Sharifzadeh, X. Ren, M. Scheffler, F. Liu, J. Lischner, L. Lin, J. R. Deslippe, and S. G. Louie, J. Chem. Theory Comput.11, 5665 (2015)
work page 2015
-
[69]
X. Ren, P. Rinke, V. Blum, J. Wieferink, A. Tkatchenko, A. Sanfilippo, K. Reuter, and M. Scheffler, New Journal of Physics14, 053020 (2012)
work page 2012
- [70]
-
[71]
R. Shi, P. Lin, M.-Y. Zhang, L. He, and X. Ren, Physical Review B109, 035103 (2024)
work page 2024
-
[72]
Low-scalinggwcalculation of quasi-particle energies within numerical atomic orbital framework,
M.-Y. Zhang, P. Lin, R. Shi, M. Scheffler, and X. Ren, “Low-scalinggwcalculation of quasi-particle energies within numerical atomic orbital framework,” (2026), arXiv preprint, in preparation
work page 2026
-
[73]
Ismail-Beigi, Journal of Physics: Condensed Matter 29, 385501 (2017)
S. Ismail-Beigi, Journal of Physics: Condensed Matter 29, 385501 (2017)
work page 2017
- [74]
-
[75]
P. Lin, X. Ren, and L. He, The Journal of Physical Chemistry Letters11, 3082 (2020)
work page 2020
- [76]
- [77]
-
[78]
P. Lin, X. Ren, and L. He, J. Chem. Theory Comput. 17, 222 (2021)
work page 2021
-
[79]
P. Lin, Y. Ji, L. He, and X. Ren, Journal of Chemical Theory and Computation21, 3394 (2025)
work page 2025
-
[80]
H. Vidberg and J. Serene, Journal of Low Temperature Physics29, 179 (1977)
work page 1977
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.