The GW approximation

The $GW$ approximation is one of the oldest and most widely used approximations for the self-energy of an interacting many-body system. Its name refers to its defining structure: the self-energy is given by the product of a propagator $G$ and a screened interaction $W$. In the language of the single boson exchange (SBE) decomposition, it arises very naturally as the approximation in which the Hedin vertices are replaced by their lowest-order contribution, $\gamma^r \simeq \mathbf{1}^r$ and $\overline{\gamma}^r \simeq \mathbf{1}^r$.

The purpose of this page is to spell out what this means concretely, in the conventions used throughout this documentation. This is worth doing carefully, because “$GW$” is used for a whole family of related but inequivalent approximations in the literature. Two choices in particular have to be made explicit, and neither of them is visible at second order in the bare interaction:

1. Which two-particle channel is resummed. The Schwinger-Dyson equation (SDE) can be written in three ways, one per channel, which are equivalent for the exact vertex $F$ but not for a truncated one. As shown below, the two particle-hole forms turn out to yield the same $GW$ self-energy, whereas the particle-particle form yields a genuinely different approximation.

2. How the spin structure is handled. For SU(2)-symmetric systems, the density and magnetic ladders in the $ph$ channel differ already at third order in $F_0$, and the self-energy picks them up with different weights.

We work in the Matsubara formalism throughout, and we specialize to a system with SU(2) spin symmetry and a local and instantaneous bare interaction, i.e. the Hubbard interaction of the Hubbard model example. This is the case in which the SBE machinery simplifies most drastically, and it is the case relevant for most applications.

Spin structure of the bare interaction¶

We first record the spin components of the antisymmetrized bare vertex of the Hubbard model, which were derived in the Hubbard model example of the section on the starting point,

where we suppressed the (trivial) frequency and momentum dependence. Reading off the three independent spin components defined in the section on spin parametrizations gives

and hence, in the density/magnetic and singlet/triplet bases,

We will also need the spin components of the $ph$ channel unit. Two related objects appear in what follows: the channel identity operator $\mathbb{1}^{ph}_{1234} = \delta_{1,2}\delta_{3,4}$, which is the unit with respect to the full connector ∘, and the unit vertex $\mathbf{1}^{ph}$, which is the unit with respect to $\bullet$, i.e. with respect to all degrees of freedom except frequency and momentum. Both carry the same spin part, $\delta_{\sigma_1, \sigma_2}\delta_{\sigma_3, \sigma_4}$, so that

From here on we consistently use $\mathbf{1}^{ph}$, since that is the object appearing in the SBE equations.

Spin projection of the loop products¶

The self-energy is obtained from a four-point object by closing two of its legs with a propagator. As discussed in the section on the Schwinger-Dyson equation, there are two such loop products,

Because the propagator is diagonal in spin, $G_{\tilde{2}\tilde{1}} \sim \delta_{\tilde{\sigma}_2, \tilde{\sigma}_1}$, and because the self-energy is diagonal in spin as well, each loop product projects the four-point object $X$ onto one specific spin combination. These two combinations are different:

The two projections look inequivalent, but for a crossing-symmetric $X$ they are not: the identity $[X\cdot G] = \zeta [G \cdot X]$ quoted in the section on the Schwinger-Dyson equation follows from $X_{1 \tilde{2} \tilde{1} 2} = \zeta X_{1 2 \tilde{1} \tilde{2}}$. For the bare vertex, which is crossing symmetric and frequency independent, we can check this on the spin components directly:

Polarization and screened interaction in the $ph$ channel¶

We now specialize to the $ph$ channel, which is the channel in which the density and magnetic fluctuations live, as is evident from its interpretation in terms of the density-density correlator discussed in the section on two-particle channels.

The $GW$ approximation consists of replacing the Hedin vertices in the SBE equations by the unit vertex,

Everything else follows. The polarization then reduces to the bare bubble,

Its spin structure is particularly simple: inserting the unit vertices into the $ph$ contraction gives $P^{ph}_{1234} = \delta_{\sigma_1, \sigma_2}\delta_{\sigma_3, \sigma_4} P^{ph}$, i.e. the polarization is proportional to the $ph$ unit vertex in spin space, with the spin-independent scalar

With the polarization at hand, the screened interaction follows from the bosonic Dyson equation,

Since $P^{ph}$ is proportional to the $ph$ unit vertex in spin space, the $\bullet$ contraction carries the same spin algebra as the ∘ contraction. That algebra is worked out explicitly in the section on spin parametrizations — the spin sums there never touch the frequency arguments, which is precisely why the result is already stated with $\bullet$ between the two factors — and it decouples in the density/magnetic basis. The corresponding spin table for all three channels is collected in the section on second-order perturbation theory. The bosonic Dyson equation therefore reduces to two independent scalar equations,

with the closed-form solutions

It is often convenient to split off the bare interaction and work with the decaying part of the screened interaction, which we denote by a tilde,

where we identified the susceptibilities in the two channels,

The $GW$ self-energy¶

We start from the $ph$ form of the Schwinger-Dyson equation,

and use the exact SBE identity $F \circ \chi_0^{r} \circ F_0 = \overline{\gamma}^{r} \bullet W^{r} - F_0$ derived in the section on the SDE in SBE form. Setting $\overline{\gamma}^{ph} \simeq \mathbf{1}^{ph}$ gives $F \circ \chi_0^{ph} \circ F_0 \simeq W^{ph} - F_0$, and therefore

which is precisely the third line of the SBE form of the SDE with the Hedin vertex set to unity. Projecting onto spin components with $[G \cdot X] \rightarrow \tfrac{1}{2}(X_d + 3X_m)$, the factor $\tfrac{1}{2}$ combines with the $\tfrac{1}{2}$ of the SDE into the weights $\tfrac{1}{4}$ and $\tfrac{3}{4}$:

where the spin sum has been carried out, so that the remaining $\cdot$ denotes only the frequency and momentum contraction of the loop.

Frequency parametrization¶

The loop closes the two legs of the screened interaction that carry the bosonic transfer frequency of the $ph$ channel. Writing out the frequency structure gives

Finally, splitting off the static contribution and using $\widetilde{W}_{d/m} = U^2 \chi_{d/m}$ yields the form in which the $GW$ self-energy is most easily recognized,

with the density per spin $n = \int_\nu G(\nu) e^{i\nu 0^+}$.

Limits and checks¶

Reduction to second-order perturbation theory¶

Expanding the screened interactions to leading order, $\widetilde{W}_{d/m} = U^2 P^{ph} + \mathcal{O}(U^3)$, both spin channels contribute identically, and the weights add up to one. For the dynamic part of the self-energy we thus obtain

Written out, and using $\zeta = -1$,

which is the familiar second-order (dynamic) self-energy of the Hubbard model, with the static Hartree term $Un$ accounted for separately.

Stoner criterion¶

The magnetic screened interaction $W_m = U / (1 - U P^{ph})$ diverges when $U P^{ph}(\omega) = 1$, while the density one, $W_d = -U/(1 + U P^{ph})$, stays finite for $P^{ph} > 0$. At $\omega = 0$ and for fermions,

which, if evaluated with the bare propagator, equals $\partial n / \partial \mu$ and is positive. The divergence condition then becomes $U \, \partial n / \partial\mu = 1$, which is the Stoner criterion for a magnetic instability. That the instability appears in the magnetic and not in the density channel is a direct consequence of the relative sign $F_{0,m} = -F_{0,d}$, and is a useful sanity check on the signs above.

Comparison with the literature¶

The result agrees with the $GW$ equations for the Hubbard atom given in Sec. 4.2 of Kiese et al., SciPost Phys. Codebases 24 (2024), who write

Their polarization does not carry the closed-loop sign, so $P = -P^{ph}$, and consequently their screened interactions are related to ours by an overall sign,

With $\tfrac{1}{4}F_{0,d} + \tfrac{3}{4}F_{0,m} = U/2$ accounting for their explicit $Un/2$, the two expressions are identical term by term.

Which form of the Schwinger-Dyson equation?¶

The SDE can be written in three equivalent ways, and the SBE identities turn each of them into a $GW$-like expression once the Hedin vertices are set to unity,

For the exact vertex the three forms of the SDE are identical. The three expressions above, however, are not related by an exact identity: each of them involves the replacement $\gamma^r \simeq \mathbf{1}^r$ made in a specific channel, and discarding $T^r \circ \chi_0^r \circ \mathbf{1}^r$ is a different approximation for different $r$. So the question of which form to start from is a real one.

It has a partly reassuring answer:

• the $ph$ form resums density and magnetic fluctuations, and is the one derived above;

• the $\overline{ph}$ form resums the crossed particle-hole ladder. It decouples not in the density/magnetic basis but in the combinations $X_\pm = X_{\uparrow\uparrow} \pm X_{\overline{\uparrow\downarrow}}$, for which the bare interaction takes the values $F_{0,+} = +U$ and $F_{0,-} = -U$. Nevertheless it gives exactly the same self-energy as the $ph$ form, to all orders;

• the $pp$ form resums the particle-particle ladder, i.e. it is a $T$-matrix approximation. In the singlet/triplet basis the bare interaction is $F_{0,s} = -2U$ and $F_{0,t}=0$, so only the singlet channel is screened. It agrees with the other two at $\mathcal{O}(U^2)$, as it must, but differs from $\mathcal{O}(U^3)$ onwards.

Variants: $G_0 W_0$ and self-consistent $GW$¶

The equations above do not yet specify which propagator enters the bubble $\chi_0^{ph}$ and the self-energy loop. Two common choices are:

• $G_0 W_0$ (one-shot): the bare propagator $G_0$ is used everywhere, so $\chi_0^{ph}$ is built from $[\chi_0^0]^{ph}$ and the self-energy is evaluated once. The result is not a conserving approximation, but it is cheap and avoids the self-consistency loop.

• self-consistent $GW$: the full propagator $G$, obtained from the Dyson equation with the current self-energy, is used in both the bubble and the loop, and the equations are iterated to convergence. This is a $\Phi$-derivable, conserving approximation in the sense of Baym and Kadanoff.

Intermediate variants (self-consistency in the self-energy loop but not in the polarization, or vice versa) exist as well and are sometimes labelled $GW_0$ or $G_0W$.