Exact two-electron spectrum and interaction effects in the Lieb lattice
Phys. Rev. B 114, 105103 – Published 3 August, 2026
DOI: https://doi.org/10.1103/1fdj-pbwz
Abstract
We present an exact solution of the two-electron problem on the Lieb lattice and analyze how the combination of a flat band, sublattice structure, and local interactions shapes the two-particle excitation spectrum. The separable form of the on-site interaction reduces the problem to a secular equation whose roots yield a triplet of discrete modes. At the point these modes form bonding, nonbonding, and antibonding combinations of sublattice amplitudes, while at the point they become strictly site selective due to the vanishing of intersublattice correlations. Along high-symmetry lines, destructive interference suppresses the bandwidth of one branch to order , providing a clear signature of flat-band locality. The interaction-induced redistribution of spectral weight is reflected in the two-electron density of states, where resonances evolve into discrete bound states that constitute the low-density precursors of the upper Hubbard band. The flat-band peak remains pinned by symmetry, but virtual transitions into the dispersive bands deplete its intensity and generate an asymmetric infrared continuum with a logarithmic divergence. Together, these results provide a complete and analytically controlled description of two-electron physics on the Lieb lattice and illustrate how flat-band singularities and local interactions conspire to produce discrete excitations, anomalously weak dispersion, and marginal infrared behavior.