Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2609.38247

Accidental accuracy and formal consistency in $GW$+BSE: Exact benchmarks and regime-dependent error cancellation

Michael O. Atambo

Published Sep 29, 2026• 0 citations
arXivHTMLPDF

Abstract

Standard many-body perturbation theory for optical excitations combines the $GW$ approximation for quasiparticle energies with the Bethe-Salpeter equation (BSE) for the electron-hole response. Formally, the BSE kernel must equal the functional derivative of the self-energy, $K=δΣ/δG$. In routine practice, however, the kernel is built from a statically screened direct interaction and a bare exchange interaction, which breaks this consistency. Using exact diagonalization of the extended Hubbard dimer, we construct a controlled benchmark that isolates the price of this inconsistency. We show that the internally inconsistent practical $GW$+BSE construction frequently yields accurate optical gaps through accidental cancellation between an underestimated $GW$ gap and an oversized bare-exchange kernel, and we map the regions of parameter space where this accidental accuracy occurs. Upon switching on the nearest-neighbor interaction and binding the exciton, the practical construction remains the accurate one throughout the weak-binding regime, with cancellation optima at parts in $10^{3}$, while the frozen-$W$ construction fails at every binding strength. Beyond a binding threshold both static constructions deteriorate, showing that deeply bound excitons require physics beyond any static kernel. Accidental accuracy is therefore a property of the weak-binding regime, and the systematic failure of the frozen-$W$ derivative kernel within the explored parameter space traces to its overscreening of the exchange channel. The accuracy of the practical kernel is an accidental cancellation of two independent errors, an underestimated quasiparticle gap and an oversized exchange term, not a hidden consistency. These results provide a controlled taxonomy of approximation failure in $GW$+BSE and a precise diagnostic for when standard workflows can be trusted.

DOI
10.48550/arxiv.2609.38247

Related papers

Accidental accuracy and vertex corrections in $GW$: Exact benchmarks for the extended Hubbard model

Michael O. Atambo

The $GW$ approximation is the standard tool for quasiparticle predictions in materials, yet its regime of validity in correlated systems remains poorly quantified, because \textit{ab initio} vertex corrections are computationally prohibitive. Using exact diagonalization of the half-filled extended Hubbard model on finite rings as a numerically exact reference, we construct the corresponding model-space $GW$ theory on the identical Hilbert space and quantify its error as a function of local ($U$) and non-local ($V$) interaction strength. We find that the required vertex correction changes character across the phase diagram: in the weak-coupling regime the effective vertex $Γ_{\rm eff} < 1$, reflecting the suppression of RPA charge fluctuations by exact short-range correlations, whereas in the Mott regime $Γ_{\rm eff}$ grows monotonically (to $\sim 3$ at $U=8t$ for $N=6$), reflecting the local dynamical self-energy structure required to open the Hubbard gap. Vertex corrections in the electron-hole (polarizability) channel are shown to \emph{worsen} the gap error, indicating that reproducing the Mott gap requires dynamical self-energy structure rather than improved screening. For $V=0$, static COHSEX is accidentally exact at a single crossover $U^* \approx 3.5\,t$; finite $V$, through non-local Fock exchange, splits this point into a double-crossover window that collapses toward weak coupling. Even at the crossover, however, the exact spectral function retains Hubbard-band structure that no static functional reproduces, so gap agreement does not imply functional accuracy. These results yield quantitative diagnostics for the reliability of $GW$ in correlated materials.

arXiv:2608.24838•Aug 25, 2026•0 citations

Scrutinizing GW-Based Methods Using the Hubbard Dimer

S. Di Sabatino, Pierre‐François Loos, Pina Romaniello

Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.

arXiv:2108.00311•Frontiers in Chemistry•Oct 29, 2021•26 citations

Correlation beyond the random phase approximation: A consistent many-body perturbation theory approach

F. Bechstedt

The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function $L$. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential $W$ in the BSE for $L$ leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.

Physical review. B./Physical review. B•Jun 15, 2018•5 citations

Total Energy beyond GW: Exact Results and Guidelines for Approximations

Abdallah El-Sahili, Francesco Sottile, Lucia Reining

The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.

arXiv:2312.04154•Journal of Chemical Theory and Computation•Feb 7, 2024•5 citations

Exact benchmarks for the plasmon-pole approximation: Multi-scale screening and the breakdown of the quasiparticle picture

Michael O. Atambo

The $GW$ approximation is the gold standard for calculating quasiparticle band structures, yet its computational cost frequently necessitates the use of the plasmon-pole approximation (PPA). While PPA is known to be highly accurate for simple metals and weakly correlated semiconductors, its regime of validity in materials with competing energy scales remains poorly quantified. Here, we construct an exact, numerical benchmark of the PPA using the Lehmann representation of the density response on finite one-dimensional Hubbard clusters. By analytically convolving the exact non-interacting Green's function $G_0$ with the exact pole representation of the screened interaction $W$, we compute the $GW$ self-energy without numerical frequency integration. We demonstrate that in single-scale Mott insulators, the moment-conserving PPA is essentially exact. However, in multi-band semiconductors where interband transitions introduce a low-energy screening channel that competes with high-energy Mott fluctuations, the PPA systematically misjudges the quasiparticle gap by several electron-volts. Furthermore, we show that strong multi-pole screening can drive the exact quasiparticle weight $Z \to 0$, destroying the quasiparticle picture-an effect entirely missed by the PPA, which artificially stabilizes sharp quasiparticles. Finally, we propose a computationally inexpensive diagnostic based on the polydispersity of the loss function's spectral weight, which accurately predicts PPA failure \textit{ab initio}. Our results provide a rigorous framework for assessing the validity of dynamical screening approximations in strongly correlated materials.

arXiv:2609.29635•Aug 29, 2026•0 citations

A formally exact one-frequency-only Bethe-Salpeter-like equation. Similarities and differences between GW +BSE and self-consistent RPA

Valério Olevano, Julien Toulouse, Peter Schuck

A formally exact Bethe-Salpeter-like equation for the linear-response function is introduced with a kernel which depends only on the one frequency of the applied field. This is in contrast with the standard Bethe-Salpeter equation (BSE) which involves multiple-frequency integrals over the kernel and response functions. From the one-frequency kernel, known approximations are straightforwardly recovered. However, the present formalism lends itself to more powerful approximations. This is demonstrated with the exact analytical solution of the Hubbard molecule. Similarities and differences of the GW + BSE approach with the self-consistent random-phase approximation are also discussed.

arXiv:1812.08421•The Journal of Chemical Physics•Feb 27, 2019•22 citations
Accidental accuracy and formal consistency in GWGWGW+BSE: Exact benchmarks and regime-dependent error cancellation
Standard many-body perturbation theory for optical excitations combines the GWGWGW approximation for quasiparticle energies with the Bethe-Salpeter equation (BSE) for the electron-hole response. Formally, the BSE kernel must equal the functional derivative of the self-energy, K=δΣ/δGK=δΣ/δGK=δΣ/δG. In routine practice, however, the kernel is built from a statically screened direct interaction and a bare exchange interaction, which breaks this consistency. Using exact diagonalization of the extended Hubbard dimer, we construct a controlled benchmark that isolates the price of this inconsistency. We show that the internally inconsistent practical GWGWGW+BSE construction frequently yields accurate optical gaps through accidental cancellation between an underestimated GWGWGW gap and an oversized bare-exchange kernel, and we map the regions of parameter space where this accidental accuracy occurs. Upon switching on the nearest-neighbor interaction and binding the exciton, the practical construction remains the accurate one throughout the weak-binding regime, with cancellation optima at parts in 10310^{3}103, while the frozen-WWW construction fails at every binding strength. Beyond a binding threshold both static constructions deteriorate, showing that deeply bound excitons require physics beyond any static kernel. Accidental accuracy is therefore a property of the weak-binding regime, and the systematic failure of the frozen-WWW derivative kernel within the explored parameter space traces to its overscreening of the exchange channel. The accuracy of the practical kernel is an accidental cancellation of two independent errors, an underestimated quasiparticle gap and an oversized exchange term, not a hidden consistency. These results provide a controlled taxonomy of approximation failure in GWGWGW+BSE and a precise diagnostic for when standard workflows can be trusted.
Accidental accuracy and vertex corrections in GWGWGW: Exact benchmarks for the extended Hubbard model
The GWGWGW approximation is the standard tool for quasiparticle predictions in materials, yet its regime of validity in correlated systems remains poorly quantified, because \textit{ab initio} vertex corrections are computationally prohibitive. Using exact diagonalization of the half-filled extended Hubbard model on finite rings as a numerically exact reference, we construct the corresponding model-space GWGWGW theory on the identical Hilbert space and quantify its error as a function of local (UUU) and non-local (VVV) interaction strength. We find that the required vertex correction changes character across the phase diagram: in the weak-coupling regime the effective vertex Γeff<1Γ_{\rm eff} < 1Γeff​<1, reflecting the suppression of RPA charge fluctuations by exact short-range correlations, whereas in the Mott regime ΓeffΓ_{\rm eff}Γeff​ grows monotonically (to ∼3\sim 3∼3 at U=8tU=8tU=8t for N=6N=6N=6), reflecting the local dynamical self-energy structure required to open the Hubbard gap. Vertex corrections in the electron-hole (polarizability) channel are shown to \emph{worsen} the gap error, indicating that reproducing the Mott gap requires dynamical self-energy structure rather than improved screening. For V=0V=0V=0, static COHSEX is accidentally exact at a single crossover U∗≈3.5 tU^* \approx 3.5\,tU∗≈3.5t; finite VVV, through non-local Fock exchange, splits this point into a double-crossover window that collapses toward weak coupling. Even at the crossover, however, the exact spectral function retains Hubbard-band structure that no static functional reproduces, so gap agreement does not imply functional accuracy. These results yield quantitative diagnostics for the reliability of GWGWGW in correlated materials.
Using the simple (symmetric) Hubbard dimer, we analyze some important features of the GW approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot GW method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard GW approximation and find, in particular, that 1) some neutral excitation energies become complex when the electron-electron interaction U increases, which can be traced back to the approximate nature of the GW quasiparticle energies; 2) the BSE formalism yields accurate correlation energies over a wide range of U when the trace (or plasmon) formula is employed; 3) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); 4) the trace formula has the correct behavior for weak (i.e., small U) interaction, unlike the ACFDT expression.
The correlation energy of an inhomogeneous electron gas is ruled by the density correlation function LLL. It obeys a Bethe-Salpeter equation (BSE), whose kernel is the functional derivative of the particle self-energy. It is solved on equal footing for two-particle and single-particle Green's functions within Hedin's GW approximation. According to the Hartree and exchange-correlation contributions to the kernel, the correlation expression decays into random-phase-approximation (RPA) -like and screened-exchange-like parts. The actual treatment of the screened potential WWW in the BSE for LLL leads to different second-order screened exchange corrections (SOSEX). Including correlation in the kernel even yields approximations beyond SOSEX. The progress beyond RPA is illustrated for spinless and spin-polarized homogeneous electron gases, indicating that the screening treatment is substantial.
The total energy and electron addition and removal spectra can, in principle, be obtained exactly from the one-body Green's function (GF). In practice, the GF is obtained from an approximate self-energy. In the framework of many-body perturbation theory, we derive different expressions that are based on an approximate self-energy, but that yield nevertheless, in principle, the exact exchange-correlation contribution to the total energy for any interaction strength. Response functions play a crucial role, which explains why, for example, ingredients from time-dependent density functional theory can be used to build these approximate self-energies. We show that the key requirement for obtaining exact results is the consistent combination of ingredients. Also when further approximations are made, as it is necessary in practice, this consistency remains the key to obtain good results. All findings are illustrated using the exactly solvable symmetric Hubbard dimer.
The GWGWGW approximation is the gold standard for calculating quasiparticle band structures, yet its computational cost frequently necessitates the use of the plasmon-pole approximation (PPA). While PPA is known to be highly accurate for simple metals and weakly correlated semiconductors, its regime of validity in materials with competing energy scales remains poorly quantified. Here, we construct an exact, numerical benchmark of the PPA using the Lehmann representation of the density response on finite one-dimensional Hubbard clusters. By analytically convolving the exact non-interacting Green's function G0G_0G0​ with the exact pole representation of the screened interaction WWW, we compute the GWGWGW self-energy without numerical frequency integration. We demonstrate that in single-scale Mott insulators, the moment-conserving PPA is essentially exact. However, in multi-band semiconductors where interband transitions introduce a low-energy screening channel that competes with high-energy Mott fluctuations, the PPA systematically misjudges the quasiparticle gap by several electron-volts. Furthermore, we show that strong multi-pole screening can drive the exact quasiparticle weight Z→0Z \to 0Z→0, destroying the quasiparticle picture-an effect entirely missed by the PPA, which artificially stabilizes sharp quasiparticles. Finally, we propose a computationally inexpensive diagnostic based on the polydispersity of the loss function's spectral weight, which accurately predicts PPA failure \textit{ab initio}. Our results provide a rigorous framework for assessing the validity of dynamical screening approximations in strongly correlated materials.
A formally exact Bethe-Salpeter-like equation for the linear-response function is introduced with a kernel which depends only on the one frequency of the applied field. This is in contrast with the standard Bethe-Salpeter equation (BSE) which involves multiple-frequency integrals over the kernel and response functions. From the one-frequency kernel, known approximations are straightforwardly recovered. However, the present formalism lends itself to more powerful approximations. This is demonstrated with the exact analytical solution of the Hubbard molecule. Similarities and differences of the GW + BSE approach with the self-consistent random-phase approximation are also discussed.