---
title: Gap Renormalization Model
url: https://www.emergentmind.com/topics/gap-renormalization-model
type: topic
---

# Gap Renormalization Model

Searching arXiv for recent and relevant papers on gap renormalization models across the usages present in the data block.
{"query":"all:\"gap renormalization\" OR ti:\"gap renormalization\" OR abs:\"band gap renormalization\" OR abs:\"superconducting gap renormalization\"","max_results":10,"sort_by":"submittedDate","sort_order":"descending"}
{"query":"id:1912.10453 OR id:2503.09715 OR id:1511.05685 OR id:1004.4834","max_results":10,"sort_by":"relevance","sort_order":"descending"}
Gap Renormalization Model denotes a set of constructions in which a gap is recomputed after incorporating renormalizing effects from self-energies, screening, quasiparticle weights, running couplings, or block projections. In the usage anchored by iron-based superconductors, the model supplements the usual linearized gap equation by a Fermi-surface–dependent quasiparticle weight \(Z_k\) extracted from a one-loop spin-fluctuation self-energy, so that quasiparticle renormalization enters directly into the pairing kernel at \(T_c\) [1912.10453]. In the literature summarized here, the same expression is also attached to band-gap narrowing in photo-excited semiconductors [2503.09715], zero-point and temperature-dependent band-gap renormalization [2007.09603; 2011.12765; 1604.02394], doping- and environment-driven quasiparticle-gap shifts in low-dimensional materials [1408.5815; 1712.04854; 2306.04759; 1308.6610], renormalization-group treatments of mass gaps and dense spectral gaps [1004.4834; 1511.05685], and rigorous renormalization procedures for gap mappings and gapped spin chains [1907.07630; 2111.09358].

## 1. Superconducting gap renormalization via momentum-dependent quasiparticle weight

In the iron-based-superconductor formulation, one works in the normal state at \(T\approx T_c\) and computes the lowest-order spin-fluctuation self-energy
\[
\Sigma(k,i\omega_n)  =  T \sum_{q,\Omega_m}  V_{\rm sf}(q,i\Omega_m)\, G_0(k-q,i\omega_n-i\Omega_m)\,.
\]
After analytic continuation and expansion for small \(\omega\),
\[
\Sigma(k,\omega+i0^+) \simeq \Sigma(k,0) + (1-Z_k^{-1})\omega + \cdots,
\qquad
Z_k = \Bigl[1-\partial_\omega \Sigma(k,\omega)\big|_{\omega\to 0}\Bigr]^{-1}.
\]
The static part \(\Sigma(k,0)\) is absorbed into a re-fitted Fermi surface, while the residual \(Z_k<1\) varies around each Fermi pocket and suppresses spectral weight most strongly where spin-fluctuation scattering is peaked [1912.10453].

The corresponding dimensionless pairing-strength functional is written as
\[
\Lambda[\Delta]
=
-
\sum_{\mu k}\sum_{\nu k'}
\Delta_\mu^*(k)\,V_{\mu\nu}(k,k')\,Z_\nu(k')\,\Delta_\nu(k')
\times
\frac{1}{2|\xi_\nu(k')|},
\]
subject to
\[
\sum_{\mu k} Z_\mu(k)\,\frac{|\Delta_\mu(k)|^2}{2|\xi_\mu(k)|}=1.
\]
Its stationarity condition yields the generalized gap equation
\[
\lambda\,\Delta_\mu(k)
=
-
\sum_{\nu,k'}
V_{\mu\nu}(k,k')\,
\frac{Z_\nu(k')\,\Delta_\nu(k')}{2|\xi_\nu(k')|},
\]
with \(\lambda=1\) at \(T=T_c\). In weak-coupling language this is often rewritten in the more familiar form
\[
\lambda\,\Delta_\mu(k)
=
-
\sum_{\nu,k'}
V_{\mu\nu}(k,k')\,
\frac{Z_\nu(k')\,\Delta_\nu(k')}{2E_\nu(k')}
\tanh\!\bigl[E_\nu(k')/(2T_c)\bigr].
\]
The central new feature is the extra factor \(Z_\nu(k')\) inside the sum.

The derivation uses a Fermi-surface restriction in momentum integrals, a static pairing interaction, one-loop self-energy only, and intra-band pairing only. Its physical interpretation is direct: since the scattering rate is largest for Fermi-surface hot spots connected by \(Q=(\pi,0)/(0,\pi)\), one obtains
\[
Z_k = [1+\Gamma(k)/\omega_0]^{-1},
\]
so regions with strong nesting acquire the smallest \(Z_k\) and contribute less to the pairing kernel. Applied to LiFeAs, the framework reproduces the observed anisotropic gap structure: ARPES and quasiparticle-interference imaging find a gap varying by up to \(\sim 30\%\) around each pocket, and inclusion of \(Z_k\) yields an angular dependence
\[
\Delta(\phi)\simeq \Delta_0[1+A\cos 4\phi]
\]
with the sign and magnitude of the observed four-fold anisotropy on the \(\Gamma\)-centered pockets.

## 2. Dyson-equation and electron-phonon formulations of band-gap renormalization

For photo-excited semiconductors, the model begins from a two-band periodic Hamiltonian with long-range Coulomb interactions and treats dynamical correlations in the non-equilibrium Green’s-function formalism using the real-time Dyson expansion. The dressed quasiparticle energies satisfy
\[
E_{c/v}^*(k)=E_{c/v}^0(k)+\Re\bigl[\Sigma_{e/h}(k,\omega=E_{c/v}^*(k))\bigr],
\]
and the renormalized gap is
\[
E_g^*=E_g^0+\Re\,\Sigma_e(E_c^*)-\Re\,\Sigma_h(E_v^*).
\]
The self-energies depend on photodoping through the non-equilibrium occupations, while screening is controlled by
\[
\Pi(q,0)=\sum_k \frac{f_{k+q}-f_k}{E_{k+q}-E_k},
\qquad
\varepsilon(q)=1+V(q)\Pi(q,0).
\]
A heavier effective mass yields weaker polarization, reduced screening, and stronger self-energy shifts. The resulting \(E_g^*(n_{ph})\) is non-monotonic, with a turnover point defined by \(dE_g^*/dn_{ph}=0\); numerically, the real-time Dyson-expansion calculations locate the critical density \(n_c\) in the range \(3\)–\(5\%\) photodoping, and increasing effective mass leads to as much as a \(6\%\) enhancement in band-gap renormalization. The optical gap follows the same trend but is altered by the Moss-Burstein effect, with up to a \(\sim 9\%\) total variation in \(\Delta_{\rm opt}\) between heavy and light mass limits [2503.09715].

In crystalline naphthalene, the relevant self-energy is the electron-phonon self-energy in lowest-order many-body perturbation theory,
\[
\Sigma^{ep}_{n\mathbf k}(\omega,T)=\Sigma^{\rm Fan}_{n\mathbf k}(\omega,T)+\Sigma^{\rm DW}_{n\mathbf k}(T).
\]
Within the on-the-mass-shell approximation,
\[
\varepsilon_{n\mathbf k}(T)=\varepsilon^0_{n\mathbf k}+\Re\,[\Sigma^{ep}_{n\mathbf k}(\varepsilon^0_{n\mathbf k},T)],
\]
and the zero-point renormalization of the fundamental gap is
\[
\Delta E_g^{\rm ZPR}
=
\Re\,\Sigma^{ep}_{c,\mathbf k_c}(\varepsilon^0_{c,\mathbf k_c},0)
-
\Re\,\Sigma^{ep}_{v,\mathbf k_v}(\varepsilon^0_{v,\mathbf k_v},0).
\]
An eigenvalue-self-consistent cycle sums high-order non-crossing diagrams implicitly. For the gap edges, the evSC shifts differ only moderately from the one-shot OMS values (\(\Delta<20\) meV), and the calculated phonon renormalization of the \(GW\)-corrected quasiparticle band structure predicts a fundamental band gap of \(5\) eV for naphthalene at room temperature [2007.09603].

Large-scale first-principles calculations beyond the adiabatic approximation show that zero-point band-gap renormalization is often larger than \(0.3\) eV, and up to \(0.7\) eV, in materials with light elements. For infrared-active materials, global agreement with available experimental data is obtained only when non-adiabatic effects are taken into account, and a generalized Fröhlich model with multiple phonon branches, anisotropic extrema, and degenerate extrema reproduces a large fraction of the full non-adiabatic shift in highly ionic and polar materials [2011.12765].

A complementary deterministic route is the one-shot Williams-Lax method, in which the stochastic sampling of nuclear wavefunctions is replaced by a single optimal supercell configuration. In the large-supercell limit, one calculation captures the temperature-dependent band gap renormalization including quantum nuclear effects in direct and indirect-gap semiconductors, while also encompassing phonon-assisted optical absorption and the degenerate-band-extrema correction identified for small-cell calculations [1604.02394].

## 3. Doping, dielectric environment, and image-charge gap renormalization in low-dimensional materials

Within a first-principles-based effective-mass model in the \(GW\) approximation for doped two-dimensional materials, the self-energy is decomposed into intrinsic, occupation/exchange, screening-induced correlation, and mixed terms,
\[
\Sigma
=
\Sigma_{\rm int}
+\Sigma_1
+\Sigma_2
+\Sigma_3
=
i\,G_{\rm int}W_{\rm int}
+i\,\delta G\,W_{\rm int}
+i\,G_{\rm int}\delta W
+i\,\delta G\,\delta W.
\]
The band-gap renormalization is written as
\[
\Delta E_g(n)=\Delta\Sigma_{\rm CH}(n)+\Delta\Sigma_{\rm SEX}(n).
\]
The Coulomb-hole contribution dominates at low doping densities, while the screened-exchange contribution dominates at high doping densities. The model was applied to h-BN, MoS\(_2\), and black phosphorus: h-BN reaches \(\Delta E_g\simeq 1\) eV at \(n\sim 10^{13}\,\mathrm{cm}^{-2}\), MoS\(_2\) gives \(\Delta E_g\sim 400\) meV at \(n\sim 10^{13}\,\mathrm{cm}^{-2}\), and anisotropic black phosphorus exhibits a particularly large renormalization because \(\Sigma_{\rm SEX}\propto 1/m_{\rm DOS}\) [1712.04854].

A related plasmon-pole treatment for monolayer MoS\(_2\) isolates the doping-induced change in the head of the inverse dielectric function,
\[
\delta \varepsilon^{-1}_{00}(q,\omega)=\frac{\Omega_d^2(q)}{\omega^2-\omega_d^2(q)},
\]
with an acoustic-like plasmon branch
\[
\omega_d(q)\simeq \sqrt{2\pi e^2 n\,q/m^*}.
\]
This approach predicts an enhanced band-gap renormalization around \(400\) meV and an unusual nonlinear evolution of the gap with doping. In the dilute limit,
\[
\Delta E_g(n)\simeq -230\,\mathrm{meV}\times (n/10^{13}\,\mathrm{cm}^{-2})^{1/2}
\qquad (n\lesssim 5\times 10^{12}\,\mathrm{cm}^{-2}),
\]
while the renormalization saturates to \(\sim -400\) meV for \(n\gtrsim 10^{13}\,\mathrm{cm}^{-2}\) [1408.5815].

When a planar 2D material is placed near an interface between two dielectrics, pseudo quantum electrodynamics and a large-\(N\) RPA treatment yield the renormalized gap as a function of carrier density,
\[
\Delta_R(n;\epsilon_1,\epsilon_2,z_0)
=
m_0
\exp\!\left[
\int_{n_0}^{n}\frac{dn'}{n'}\,C_{\bar\lambda}(\sqrt{\pi n'})
\right].
\]
The limiting cases are explicit. For \(\epsilon_1=\epsilon_2\), one recovers the single-medium exponent; for \(d\to 0\), the two-dielectric system reduces to a single effective medium with \(\epsilon_{\rm eff}=(\epsilon_1+\epsilon_2)/2\); and for \(d\to\infty\), one recovers the single-medium result with \(\epsilon=\epsilon_1\) [2306.04759].

For semiconducting carbon nanotubes near a metallic surface, the quasiparticle band-gap renormalization follows the electrostatic scaling law
\[
\Delta E_{\rm gap}^{QP}(h_a)\approx -\frac{1}{2h_a},
\]
with \(h_a\) the apparent nanotube height. The binding energy of excitons is reduced dramatically, by as much as \(75\%\), near the surface. Compensation between quasiparticle and excitonic effects leads to small changes in the optical gap,
\[
E_{\rm opt}(h_a)=E_{\rm gap}^{QP}(h_a)-E_B(h_a),
\]
typically of order \(50\)–\(100\) meV [1308.6610].

## 4. Renormalization-group and variational mass-gap constructions

In the interacting Aubry-André model, rigorous multiscale renormalization-group analysis shows that the infinitely many gaps of the single-particle spectrum persist in the presence of weak many-body interactions. For each Diophantine label \(n\), the physical gap is identified with the effective mass scale and obeys
\[
\Delta_n(U,u)\simeq \bigl(\Delta_n^{(0)}+F_n(u,U)\bigr)^{\alpha_n(U)},
\qquad
\alpha_n(U)=\frac{1}{1-\eta_\mu(U)}.
\]
In the weak-potential limit,
\[
\Delta_n(U)=Z_n(U)\,[\Delta_n^{(0)}]^{\alpha_n(U)}.
\]
The exact scaling relations include
\[
\alpha_n=1/(1-\eta_\mu),\qquad X_n\equiv \alpha_n^{-1}=1-\eta_\mu,\qquad X_n=2-K_n.
\]
Repulsive interactions give \(\alpha_n<1\) and attractive interactions give \(\alpha_n>1\); no insulator-metal transition occurs [1511.05685].

For the two-dimensional \(O(2N)\) Gross-Neveu model, RG-improved optimized perturbation theory introduces an interpolated Lagrangian and fixes the variational mass and coupling simultaneously through the principle of minimal sensitivity and the reduced RG equation. At two-loop order, the ordinary perturbative pole mass
\[
M_{\rm pert}^{(2)}(m,g)
=
m\Bigl[
1+g(c_1-\gamma_0L)
+g^2\bigl(c_2+(\gamma_0^2-\gamma_1-c_1(\gamma_0+b_0))L+\tfrac{\gamma_0}{2}(\gamma_0+b_0)L^2\bigr)
\Bigr]
\]
is converted into optimized gap equations for \((\tilde m_0,\tilde g)\). At large \(N\), the exact result is reproduced already at the very first order of the modified perturbation. For arbitrary \(N\), using perturbative information known at two-loop order, the method yields a discrepancy of order \(1\!\sim\!2\%\) relative to the exact Bethe–Ansatz mass gap [1004.4834].

In the gapped Kondo model, Wegner’s flow-equation renormalization group is applied to a conduction density of states with a hard gap,
\[
\rho_0(\omega)=\rho_0\,\Theta(|\omega|-\Delta).
\]
The projected one-loop flow satisfies
\[
\frac{dJ(l)}{dl}=-C\,f(\Delta,l)\,J(l)^2,
\]
and the flow is cut off by the gap. For \(\Delta<\Delta_c\) the coupling still runs to strong coupling, while for \(\Delta>\Delta_c\) it saturates and the system flows to the local-moment fixed point. As the gap is increased, the spin susceptibility crosses over to Curie behavior [1711.10997].

A distinct rigorous use of renormalization appears in proofs that frustration-free local spin chains are gapped. The basic map sends
\[
H=H^S+H^L
\quad\longmapsto\quad
R[H]=h=P^S\,H^L\,P^S,
\]
and the central theorem states
\[
H\text{ is gapped}
\quad\Longleftrightarrow\quad
h\text{ is gapped.}
\]
Under decay of correlations and lower bounds on the local block gaps, the method produces a uniform lower bound on the many-body spectral gap. The construction was applied to the teleportation chain and the swap chain [2111.09358].

## 5. Localized and disordered gap renormalization

For an \(s\)-wave superconductor with two magnetic impurities, the full Hamiltonian combines kinetic, pairing, and exchange terms, and the gap is determined self-consistently by
\[
\Delta(\mathbf r)=-g\,\langle c_\downarrow(\mathbf r)c_\uparrow(\mathbf r)\rangle_H.
\]
In the weak-coupling Shiba regime, the local gap renormalization at an impurity site reduces to
\[
\delta\Delta(r_1)\approx -3\,\alpha_1^2\,\Delta_0,
\qquad
\alpha_i=\pi\nu_F J_i S_{iz}.
\]
For two impurities, the orientation-dependent contribution behaves as
\[
\delta\Delta^{(12)}_{\rm orient.}\propto
\alpha_1\alpha_2\,e^{-2r_{12}/\xi}\,(k_Fr_{12})^{-2}\cos\theta.
\]
At strong exchange coupling, a self-consistent lattice Bogoliubov–de Gennes calculation shows that the local gap under an impurity can change sign; when this occurs, a sub-gap Shiba state turns into a supra-gap Andreev bound state [1501.07901].

In the random Kitaev spin ladder, the strong-disorder renormalization group distinguishes a spin gap
\[
\Delta_s=E_1-E_0
\]
from a local flux gap
\[
\Delta_f(j)=E(W_j=+1)-E(W_j=-1).
\]
In the Ising limit, the minimum spin gap obeys a Fréchet law,
\[
P(\ln u)=\frac{1}{z}\,u^{1/z}\exp[-u^{1/z}],
\qquad
u\equiv u_0\,\Delta_s\,L^z,
\]
and the infinite-disorder critical point has
\[
-\ln\Delta_s\sim L^\psi,\qquad \psi=1/2.
\]
The raw flux-gap distribution is non-universal, but after dividing out the trivial \(J_y\) prefactor the scaled distribution collapses to the universal infinite-disorder form. In the XX limit, the flux sector exhibits
\[
-\ln \Delta_f\sim L^\psi,\qquad \psi=1,
\]
because the \(x\)- and \(y\)-couplings are renormalized simultaneously while the \(z\)-couplings are not renormalized drastically [2205.06829].

## 6. Gap mappings, stochastic gap regularization, and recurring structure

In one-dimensional dynamics, a dissipative gap map is a \(C^3\) Lorenz map with two strictly increasing branches, a nonempty gap in the range, and derivative bounded by \(0<f'(x)\le \nu<1\). Renormalization is defined by taking a suitable first-return map to a central interval and rescaling,
\[
\mathcal R f(x)=\frac{1}{|I'|}\,R(|I'|x),
\qquad x\in [-1,1]\setminus\{0\}.
\]
For infinitely renormalizable maps, the tangent space splits as
\[
T_f\mathcal D = E_f^u\oplus E_f^s,
\]
with a one-dimensional unstable direction and a codimension-one stable direction. The resulting stable and unstable manifolds are \(C^1\) embedded submanifolds, and topological-conjugacy classes of infinitely renormalizable gap mappings are \(C^1\) manifolds [1907.07630].

In planar Brownian self-intersection local time, “gap renormalization” refers to a time-gap cutoff that removes the logarithmic divergence at \(t_2\downarrow t_1\). The regularized local time is
\[
L(\Lambda)
=
\int_{\substack{0<t_1<t_2<T\\ t_2-t_1>\Lambda}}
d^2 t\;
\delta\!\bigl(\mathbf B(t_2)-\mathbf B(t_1)\bigr),
\]
with divergent expectation
\[
\mathbb E[L(\Lambda)]
=
(2\pi)^{-1}\int_\Lambda^T d\tau\,\frac{T-\tau}{\tau}
=
-\frac{T}{2\pi}\ln\Lambda + O(1).
\]
Subtracting this counterterm gives
\[
L^R(\Lambda)=L(\Lambda)-\mathbb E[L(\Lambda)],
\]
and \(L^R(\Lambda)\) converges in \(L^2\) to a finite limit \(L_c\). In the planar Edwards model, this refines Varadhan’s renormalization argument and extends the integrability bound from \(g<\pi/T\) to \(g<2\pi/T\) [1505.07099].

Taken together, these constructions suggest a recurrent architecture. A bare gap equation or gap observable is first supplemented by a renormalizing quantity—\(Z_k\), \(\Sigma\), \(W\), running couplings, a projector, or an explicit counterterm. The renormalized gap is then extracted from a stationarity condition, a Dyson equation, an RG flow, a first-return operator, or a convergent regularization. This suggests that “gap renormalization model” is best understood not as a single formalism, but as a technically varied class of procedures unified by the decision to treat the gap itself as the object on which renormalization acts.

Source: https://www.emergentmind.com/topics/gap-renormalization-model