---
title: Generalized Eliashberg Approach
url: https://www.emergentmind.com/topics/generalized-eliashberg-approach
type: topic
---

# Generalized Eliashberg Approach

Generalized Eliashberg approach designates a family of extensions of conventional Migdal–Eliashberg theory in which the self-consistent treatment of normal and anomalous self-energies is preserved while one or more standard assumptions—phonon-only pairing, static Coulomb repulsion, adiabaticity, equilibrium, isotropy, Fermi-liquid normal states, or mean-field-only treatment—are relaxed. In the cited literature, the term covers ab initio Coulomb-inclusive formalisms, nonequilibrium Keldysh–Usadel theories, multiband and full-bandwidth spin-fluctuation frameworks, quantum-critical large-\(N\) constructions, external-pair-potential models, and several numerical or functional reformulations [2106.05751][1911.04959][2004.01539][2506.11952].

## 1. Formal basis

The common formal starting point is the Nambu–Gorkov Green’s function and Dyson equation,
\[
{\bm G}={\bm G}_0+{\bm G}_0{\bm \Sigma}{\bm G},
\]
together with a decomposition of the superconducting self-energy into diagonal and anomalous components. A representative form is
\[
{\bm \Sigma}(i\omega_j)=i\omega_j(1-Z(i\omega_j))\sigma_0+\phi(i\omega_j)\sigma_1+(\chi(i\omega_j)-V_{\rm eff})\sigma_3,
\]
where \(Z\) is the renormalization function, \(\phi\) the anomalous pairing self-energy, and \(\chi\) the diagonal energy-shift or band-renormalization term [2106.05751]. In time-reversal-symmetric cases, \(Z\) is identified with the odd-in-frequency part of the normal self-energy and \(\chi\) with the even-in-frequency part:
\[
i\omega_j(1-Z(i\omega_j))=\frac{1}{2}\left[\Sigma(i\omega_j)-\Sigma(-i\omega_j)\right],
\]
\[
\chi(i\omega_j)-V_{\rm eff}=\frac{1}{2}\left[\Sigma(i\omega_j)+\Sigma(-i\omega_j)\right].
\]
This separation is not merely notational. Several generalized formulations turn on the fact that \(Z\) multiplies frequency while \(\chi\) shifts the dispersion, so they renormalize observables in inequivalent ways [2106.05751].

A functional-integral formulation makes the same structure explicit from another direction. After integrating out phonons and decoupling the retarded electron-electron interaction, both the anomalous pairing function and the normal self-energy appear as Hubbard–Stratonovich fields. The mean-field Eliashberg equations then arise as saddle-point conditions of a bosonic effective action, with schematic self-consistency relations
\[
\phi(k)=\sum_{k'}\lambda(k-k')F(k'), \qquad
\Sigma(k)=\sum_{k'}\lambda(k-k')G(k'),
\]
and Gaussian fluctuations organized around Cooper and density/self-energy channels [2012.12287]. In that sense, generalized Eliashberg theory is often less a single model than a self-consistent architecture for frequency-dependent normal and anomalous sectors.

## 2. Major directions of generalization

The phrase “generalized Eliashberg approach” is used in several distinct but related ways.

| Direction | Characteristic modification | Representative work |
|---|---|---|
| External pair potential | \(\widetilde W=W+\frac{W}{|W|}\frac{\upsilon}{2}\) and \(\widetilde\Delta=\Delta+\frac{\upsilon}{2}\frac{\Delta}{|\Delta|}\frac{1}{Z}\) | [1705.08871] |
| Microwave-driven dirty superconductors | Generalized self-consistency \({\cal F}_{\rm eq}+\alpha{\cal F}_{\rm neq}=0\) at arbitrary \(T,\omega,\Gamma,\gamma_{\rm in}\) for weak power | [1911.04959] |
| Full-bandwidth spin fluctuations | Momentum-, band-, and Matsubara-dependent \(Z,\Gamma,\phi\) over the full Brillouin zone | [2004.01539] |
| Quantum-critical and NFL pairing | Coupled \(\Sigma\), \(\Phi\), and \(\Pi\) equations with critical bosons and large-\(N\) control | [2506.11952] |
| Nonadiabatic finite-bandwidth extensions | Separation of pairing \(\lambda\) from mass-renormalization \(\tilde\lambda\) | [1908.00718] |
| Dynamical Coulomb/plasmonic kernels | Frequency-dependent \(W(i\nu)\) and hybrid Eliashberg–SCDFT treatments | [2007.12886] |

In the external-pair-potential construction, the parameter \(\upsilon\) is not itself a pairing interaction in the usual sense. Rather, it shifts the anomalous self-energy or gap sector directly. The resulting equations remain Eliashberg-like, but the anomalous combination entering the propagators is \(\widetilde W\) or \(\widetilde\Delta\), not the standard \(W\) or \(\Delta\). Within the model analyzed there, inclusion of the Coulomb pseudopotential restores a finite \(T_c\), and the gap vanishes linearly as \(T\to T_c^-\), \(|\Delta|\propto T_c-T\), rather than with the standard mean-field square-root law [1705.08871].

In multiorbital unconventional superconductors, the generalization is often toward full momentum and orbital resolution. For FeSe, the interaction is constructed from multiorbital Hubbard–Hund terms, RPA spin and charge susceptibilities, and then transformed from orbital to band space before entering coupled equations for \(Z_{\mathbf{k}n}\), \(\Gamma_{\mathbf{k}n}\), and \(\phi_{\mathbf{k}n}\) [2004.01539]. In dirty superconductors under microwave irradiation, the same phrase denotes a nonequilibrium weak-power extension in which spectral depairing and quasiparticle redistribution are treated on equal footing [1911.04959]. In quantum-critical settings, it denotes self-consistent pairing out of non-Fermi liquids rather than quasiparticle metals [2506.11952].

## 3. Diagonal self-energies, Coulomb renormalization, and normalization

A major contemporary use of the generalized Eliashberg approach concerns the explicit treatment of the screened Coulomb interaction. A central claim of Akashi’s homogeneous-electron-gas analysis is that a properly “normed” ab initio Eliashberg theory must reproduce the homogeneous electron gas limit, and in that limit the Coulomb self-energy cannot be treated consistently if one keeps only the usual \(Z\)-type renormalization while dropping the diagonal \(\chi\) term [2106.05751]. The generalized Eliashberg equations already show the asymmetry:
\[
Z_i(i\omega_j)\ \text{renormalizes } \omega_j,\qquad
\chi_i(i\omega_j)\ \text{shifts } \xi_i.
\]
In the homogeneous electron gas, both contributions are individually large and largely cancel in the quasiparticle velocity, but they do not cancel in the pairing kernel. The paper therefore concludes that neglect of \(\chi\) breaks the normalization of the theory already in the reference electron gas [2106.05751].

That distinction between diagonal and anomalous sectors reappears in a different setting in the pseudogap study based on combined electron-phonon and electron-electron-phonon interactions. There, the density of states is argued to be gapped primarily by the diagonal self-energy components \(Z(\omega)\) and \(\chi(\omega)\), while the direct contribution of the anomalous self-energy \(\varphi(\omega)\) is negligible. In this construction the pseudogap persists above \(T_c\) because \(Z\) and \(\chi\) remain strongly frequency dependent in the anomalous normal state [1503.06932]. This suggests a broader lesson: generalized Eliashberg theories often cease to identify “the gap” with the anomalous sector alone.

Several recent ab initio extensions reach a similar conclusion from the Coulomb side. A plasmonic extension of Eliashberg theory replaces the static Coulomb kernel by a dynamically screened interaction \(W(i\nu)\), represented through an electronic spectral density \(\mathcal S\). This introduces Coulomb contributions both to the anomalous kernel and to the mass renormalization \(Z^c\). Numerically, however, fully plasmonic Eliashberg theory substantially overestimates plasmon-mediated pairing and \(T_c\), whereas a hybrid scheme that keeps Eliashberg for the electron-phonon sector and uses an SCDFT-like treatment for the dynamical Coulomb sector gives results close to SCDFT and in excellent agreement with experiment [2007.12886]. A related beyond-RPA extension replaces the usual test-charge screened Coulomb interaction by the Kukkonen–Overhauser effective interaction. In a set of conventional superconductors, the resulting vertex corrections systematically decrease \(T_c\), from a few percent in bulk lead to more than \(40\%\) in compressed lithium [2304.14682]. Taken together, these works constrain generalized Eliashberg practice: making Coulomb dynamics explicit is necessary, but doing so within a naive \(GW\)-like anomalous self-energy can be quantitatively unreliable.

## 4. Nonadiabatic, antiadiabatic, and weak-coupling limits

Another major axis of generalization concerns the breakdown of the adiabatic assumption. In the nonadiabatic or antiadiabatic regime, where characteristic phonon frequencies are comparable to or larger than the Fermi energy or bandwidth, the standard identification of the Eliashberg coupling with both pairing strength and mass renormalization ceases to hold. The generalized mass-renormalization coupling becomes
\[
\tilde\lambda = 2\int_0^\infty \frac{d\omega}{\omega}\,\alpha^2(\omega)F(\omega)\frac{D}{\omega+D},
\]
while the pairing strength remains controlled by \(\lambda\) [1908.00718]. In the antiadiabatic limit, \(\tilde\lambda\sim \lambda D/\Omega_0\) is parametrically small, and the relevant perturbative parameter becomes
\[
\lambda \frac{D}{\Omega_0}\sim \lambda \frac{E_F}{\Omega_0}\ll 1.
\]
The transition temperature correspondingly crosses over from a phonon-scale prefactor to an electronic one, with
\[
T_c\sim D\,\exp\!\left(-\frac{1+\tilde\lambda}{\lambda}\right)
\]
in the antiadiabatic regime [1908.00718]. The general implication is that pairing strength, quasiparticle mass renormalization, and vertex-correction control parameter need not coincide once the phonon energy ceases to be the smallest scale.

The weak-coupling limit of conventional Eliashberg theory also differs from the BCS limit in a precise way. For an Einstein phonon with \(\mu^*=0\), the renormalized weak-coupling transition temperature is
\[
T_c=\frac{1.13}{\sqrt e}\,\omega_E \exp\!\left(-\frac{1+\lambda}{\lambda}\right),
\]
not the usual BCS prefactor. Moreover, the Matsubara-axis gap tends toward
\[
\Delta(i\omega_n)\to \frac{1}{1+(\omega_n/\omega_E)^2},
\]
so retardation survives at arbitrarily small coupling and the weak-coupling limit remains non-BCS in functional form [1807.04907]. This is a useful benchmark for generalized theories: weak coupling does not by itself justify replacing a retarded kernel by an instantaneous one.

## 5. Nonequilibrium, quantum criticality, and pairing without quasiparticles

Generalized Eliashberg theory has also been extended beyond equilibrium and beyond the Fermi-liquid normal state. In dirty superconductors subject to microwaves, the microscopic weak-power theory is formulated in Keldysh–Usadel language and organized around
\[
{\cal F}_{\rm eq}(\Delta,\Gamma,T,\gamma_{\rm in})+\alpha\,{\cal F}_{\rm neq}(\Delta,\Gamma,T,\omega,\gamma_{\rm in})=0,
\]
where \({\cal F}_{\rm neq}\) splits into a spectral part and a kinetic part [1911.04959]. The spectral correction encodes direct depairing by the ac field; the kinetic correction encodes nonequilibrium quasiparticle redistribution. The resulting phase diagram contains only a finite region of superconductivity enhancement near \(T_c\); at sufficiently high frequencies and low temperatures, direct depairing prevails and superconductivity is suppressed [1911.04959].

In quantum-critical metals and SYK-like large-\(N\) models, the generalized Eliashberg approach addresses pairing out of non-Fermi liquids. The common structure is again a coupled set of self-consistent equations, now for fermionic self-energy, anomalous self-energy, and bosonic polarization. In the Yukawa-SYK formulation, the large-\(N\) saddle-point equations take the Eliashberg form
\[
\Sigma\sim DG,\qquad \Phi\sim DF,\qquad \Pi\sim GG-FF,
\]
with the boson itself becoming critical or overdamped [2506.11952]. In related analyses of critical Fermi surfaces and SYK-like lattice models, Eliashberg-type equations are asymptotically exact because the dominant diagrams are self-consistent rainbow or melonic contributions while vertex corrections are suppressed by large \(N\), small \(\epsilon\), or locality [1912.07646]. This is the sense in which the formalism extends to “pairing without quasiparticles.”

The comparative study of phonons and soft nematic fluctuations sharpens the issue of validity. For the electron-phonon problem, the Migdal parameter
\[
\lambda_E=\lambda\,\frac{\omega_D}{E_F}
\]
controls the neglect of vertex corrections. For the nematic quantum critical point, there is no analogous microscopic fast-boson/slow-electron separation; nevertheless, the two-loop vertex corrections are argued to be not small parametrically but small numerically, while perturbative evaluation of the one-loop self-energy remains rigorously justified when the fermion-boson coupling is small compared to \(E_F\) [2404.11820]. The paper’s conclusion is unusually strong: Eliashberg theory for the nematic case can hold all the way to a quantum critical point, whereas in the electron-phonon case it breaks down at some distance from where the dressed Debye frequency would vanish [2404.11820].

## 6. Multiband, fluctuation, and computational formulations

A further meaning of “generalized Eliashberg approach” is methodological. In the functional-integral reformulation, the standard Eliashberg saddle point is embedded in a bosonic effective action whose Gaussian fluctuations generate both Cooper-channel and density-channel propagators. The inverse fluctuation kernels take the form
\[
\Gamma_\eta^{-1}=\lambda^{-1}-G_nG_n,\qquad
\Gamma_\theta^{-1}=\lambda^{-1}-2G_nG_n,
\]
making explicit that density/self-energy fluctuations are as systematic as pairing fluctuations in the generalized field theory [2012.12287]. This extension was used to analyze strong-coupling fluctuation diamagnetism near \(T_c\) [2012.12287].

A different reformulation maps nonlinear Migdal–Eliashberg theory onto a classical spin chain in Matsubara-frequency space,
\[
\mathcal H_{\rm spin}=-2\pi\sum_n \omega_n S_n^z-\pi^2 T\sum_{n,m}V_{n,m}\,\vec S_n\cdot\vec S_m.
\]
The mapping yields both a heat-bath Monte Carlo solver for the full nonlinear equations and a high-frequency decimation procedure that reproduces the Morel–Anderson pseudopotential reduction,
\[
\mu^*=\frac{\mu}{1+\mu\ln(\Lambda/\Lambda^*)},
\]
for the Bogoliuov–Tomachov–Morel–Anderson kernel [2311.12094]. Here the generalization lies in turning Eliashberg theory into a classical statistical problem whose renormalization can be carried out directly in Matsubara space.

Material-specific implementations show the same breadth. In \(\mathrm{Ba}_{0.68}\mathrm{K}_{0.32}\mathrm{Fe}_2\mathrm{As}_2\), a four-band strong-coupling Eliashberg model was reduced to an effective two-band description to explain superconductivity-induced infrared anomalies. The theory accounts for two gap scales, \(2\Delta_A\approx 6\,k_BT_c\) and \(2\Delta_B\approx 2.2\,k_BT_c\), and attributes optical suppression up to \(14\,k_BT_c\) to spin-fluctuation-assisted processes in the clean limit [1103.0938]. In FeSe, a full-bandwidth, multiband, anisotropic Eliashberg theory based on RPA spin and charge fluctuations yields \(T_c\approx 6\) K and an \(s_\pm+s\) gap for bulk FeSe, but only a \(d\)-wave solution with \(T_c\le 11\) K for monolayer FeSe/SrTiO\(_3\), thereby favoring interfacial electron-phonon coupling as the dominant pairing mediator in the monolayer system [2004.01539].

Across these formulations, the invariant core is the self-consistent treatment of frequency-dependent normal and anomalous sectors. The main unresolved issues are equally consistent across the literature: the treatment of the diagonal \(\chi\) channel, Coulomb vertex corrections beyond RPA/\(GW\), high-energy incoherent sectors, and the loss of a universal Migdal control parameter outside ordinary phonon problems [2106.05751][2304.14682][2004.01539]. The generalized Eliashberg approach is therefore best understood not as a single extension, but as a broad research program for preserving Eliashberg’s dynamical self-consistency while reworking its assumptions to fit strongly renormalized, multiscale, or nonequilibrium superconducting systems.

Source: https://www.emergentmind.com/topics/generalized-eliashberg-approach