Abstract:The $GW$ approximation is a widely used method for computing electron addition and removal energies of molecules and solids. The computational effort of conventional $GW$ algorithms increases as $O(N^4)$ with the system size $N$, hindering the application of $GW$ to large and complex systems. Low-scaling $GW$ algorithms are currently very actively developed. Benchmark studies at the single-shot $G_0W_0$ level indicate excellent numerical precision for frontier quasiparticle energies, with mean absolute deviations $<10$ meV between low-scaling and standard implementations for the widely used $GW100$ test set. A notable challenge for low-scaling $GW$ algorithms remains in achieving high precision for five molecules within the $GW100$ test set, namely O$_3$, BeO, MgO, BN, and CuCN, for which the deviations are in the range of several hundred meV at the $G_0W_0$ level. This is due to a spurious transfer of spectral weight from the quasiparticle to the satellite spectrum in $G_0W_0$ calculations, resulting in multi-pole features in the self-energy and spectral function, which low-scaling algorithms fail to describe. We show in this work that including eigenvalue self-consistency in the Green's function ($\text{ev}GW_0$) achieves a proper separation between satellite and quasiparticle peak, leading to a single solution of the quasiparticle equation with spectral weight close to one. $\text{ev}GW_0$ quasiparticles energies from low-scaling $GW$ closely align with reference calculations; the mean absolute error is only 12 meV for the five molecules. We thus demonstrate that low-scaling $GW$ with self-consistency in $G$ is well-suited for computing frontier quasiparticle energies.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to improve the accuracy of the low - scaling GW algorithm when processing specific molecules, especially for the multipole challenges of five molecules (O₃, BeO, MgO, BN, and CuCN) in the GW100 benchmark set. Specifically: 1. **The computational complexity problem of the traditional GW algorithm**: The computational complexity of the traditional GW algorithm grows as O(N⁴) with the system size N, which limits its application in large and complex systems. Therefore, the development of low - scaling GW algorithms has become a research hotspot. 2. **The accuracy problem caused by multipole characteristics**: Although the low - scaling GW algorithm performs well in most cases, when processing the above five molecules, the calculation results at the G₀W₀ level deviate significantly (several hundred meV) from the standard implementation. This is because in the G₀W₀ calculation, an incorrect transfer of spectral weights occurs between quasiparticle peaks and satellite spectra, resulting in multipole characteristics in the self - energy and spectral functions, and the low - scaling algorithm cannot correctly describe these characteristics. 3. **Methods to solve multipole characteristics**: To overcome this problem, the author proposes to achieve the correct separation between quasiparticle peaks and satellite peaks by introducing eigenvalue self - consistency (evGW₀), that is, partial eigenvalue self - consistency, in the Green's function. This method can obtain a single quasiparticle equation solution and make the spectral weight close to 1. 4. **Verify the effectiveness of the improved method**: The author shows that the low - scaling GW algorithm using evGW₀ can significantly improve the calculation accuracy of these five molecules, reducing the mean absolute error to only 12 meV, thus proving that the low - scaling GW algorithm is effective in processing frontier quasiparticle energies. ### Summary of key formulas - **GW quasiparticle equation**: \[ \varepsilon^{\text{G0W0}}_n=\varepsilon^{\text{DFT}}_n + \operatorname{Re} \Sigma^c_n(\varepsilon^{\text{G0W0}}_n)+\Sigma^x_n - v_{xc,n} \] - **Multipole characteristics of self - energy**: \[ \Sigma^c_n(\omega)=\sum_m\sum_s\frac{\langle \psi_n \psi_m | P_s | \psi_m \psi_n \rangle}{\omega - \varepsilon^{\text{DFT}}_m+(\Omega^{\text{RPA}}_s - i\eta)\operatorname{sgn}(\varepsilon_F - \varepsilon^{\text{DFT}}_m)} \] - **Pole positions after eigenvalue self - consistency correction**: \[ \omega_{\text{pole}}=\varepsilon^{\text{DFT}}_i+\Delta^{\text{QP}}_i - \Omega^{\text{RPA}}_s\quad\text{(occupied orbit)} \] \[ \omega_{\text{pole}}=\varepsilon^{\text{DFT}}_a+\Delta^{\text{QP}}_a+\Omega^{\text{RPA}}_s\quad\text{(unoccupied orbit)} \] Through these improvements, the paper successfully solves the multipole characteristic problem of the low - scaling GW algorithm when processing specific molecules and improves the calculation accuracy.

Solving multi-pole challenges in the GW100 benchmark enables precise low-scaling GW calculations

Benchmark of $\boldsymbol{GW}$ Methods for Core-Level Binding Energies

Two-component $GW$ calculations: Cubic scaling implementation and comparison of vertex corrected and partially self-consistent $GW$ variants

GW100: Benchmarking G0W0 for Molecular Systems

Scalable Molecular GW Calculations: Valence and Core Spectra

Precision benchmarks for solids: G0W0 calculations with different basis sets

Quasi-Particle Self-Consistent $GW$ for Molecules

Efficient Full-frequency GW Calculations using a Lanczos Method

Exploiting a derivative discontinuity estimate for accurate $G_0W_0$ ionization potentials and electron affinities

Joint Approximate Diagonalization approach to Quasiparticle Self-Consistent $GW$ calculations

Benchmark of GW Methods for Core-Level Binding Energies

Accelerating core-level $GW$ calculations by combining the contour deformation approach with the analytic continuation of $W$

A Density-Based Basis-Set Incompleteness Correction for GW Methods

Why does the $GW$ approximation give accurate quasiparticle energies? The cancellation of vertex corrections quantified

Quasiparticle and fully self-consistent GW methods: an unbiased analysis using Gaussian orbitals

$GW$ density matrix to estimate self-consistent $GW$ total energy in solids

Fully dynamic G3W2 self-energy for finite systems: Formulas and benchmark

A 'moment-conserving' reformulation of GW theory

Automated Workflow for Accurate High-Throughput GW Calculations

FEAST nonlinear eigenvalue algorithm for $GW$ quasiparticle equations