Interplay between Quantum Metric and Hybridized Collective Modes in Flat-Band Superfluids

Abstract: We investigate collective excitations in flat-band superfluids by incorporating the coupled dynamics of pairing (phase and amplitude) and density fluctuations. We demonstrate that for any time-reversal symmetric superfluid system with an isolated flat band, only a single low-energy collective mode emerges in the long-wavelength limit. In contrast to the linearly dispersive Goldstone mode in conventional superfluids, this hybridized mode is gapless at zero momentum but exhibits a quadratic dispersion ($ω\propto q<sup>2$) at small momenta. Analytically, we reveal that the dispersion coefficient of this collective mode is governed by the normal-state quantum metric of the flat band. These analytical predictions are in excellent agreement with numerical calculations. Our results are universally applicable to any generic $s$-wave flat-band superfluid, provided the flat band is energetically well separated from other dispersive bands.

• The paper demonstrates that retaining phase, amplitude, and density fluctuations produces one gapless hybrid mode with quadratic dispersion, ω ∝ q², in any doped, time-reversal-symmetric isolated flat band.
• A generalized Hubbard–Stratonovich treatment and Schur-complement analysis show that the mode curvature is fixed by the projected quantum metric, κ = (2E/Nₖ)Σₖgᵩᵩ(k), rather than by an effective mass.
• Lieb-lattice calculations reproduce an exponent of approximately 2.0 and confirm that tuning staggered hopping changes the mode curvature through the average quantum metric, while broken time-reversal symmetry remains an open case.

Motivation and central claim

Flat-band superfluids occupy a special position in the theory of superconductivity: because the single-particle kinetic energy is quenched and the density of states diverges, the conventional BCS relation between superfluid weight and effective mass breaks down, and the quantum geometry of the Bloch wavefunctions—specifically the normal-state quantum metric—sets a lower bound on the superfluid weight (2606.01235). The paper under discussion extends this geometric picture from static transport to dynamics: it computes the low-energy collective excitations of a generic two-dimensional $s$-wave flat-band superfluid while treating the pairing (phase and amplitude) and particle-density channels on an equal footing. Its central claims are twofold. First, for any time-reversal-symmetric (TRS) isolated flat band at finite doping ($\mu \neq 0$), only one gapless collective mode survives in the long-wavelength limit, and it disperses quadratically, $\omega \propto q^2$, rather than linearly as in conventional BCS superfluids. Second, the curvature coefficient of this hybridized mode is analytically determined by the normal-state quantum metric, $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$, with $E = \sqrt{\mu^2 + \Delta_0^2}$.

The authors argue that prior treatments are incomplete because they typically decouple only the pairing channel and integrate out density fluctuations. In dispersive superconductors this is justified by the kinetic-energy penalty for charge fluctuations and approximate particle-hole symmetry near the Fermi surface; in a flat band neither protection exists, so the density mode becomes anomalously soft and must be retained. This is consistent with exact numerical evidence for macroscopic phase separation tendencies in flat-band superfluids cited in the paper. The work thereby unifies earlier model-specific observations—a quadratic Goldstone dispersion in lowest-Landau-level functional integrals and in pseudo-magnetic flat bands—into a general statement valid for any TRS $s$-wave flat-band superfluid whose flat band is energetically separated from dispersive bands.

Model and mean-field projection

The microscopic setting is a generic multi-sublattice lattice Hamiltonian $H = H_0 + H_{\text{int}}$ with an attractive on-site interaction $-U n_{\gamma\uparrow}n_{\gamma\downarrow}$ and sublattice-dependent real $s$-wave order parameters $\Delta_\gamma$. When $\mu \neq 0$0 lies within an isolated flat band of energy $\mu \neq 0$1 and $\mu \neq 0$2 is smaller than the band gap, the BdG Hamiltonian projects onto the flat-band subspace as

$\mu \neq 0$3

where $\mu \neq 0$4 is the flat-band projector and $\mu \neq 0$5. The quasiparticle spectrum is strictly flat at $\mu \neq 0$6, so all momentum dependence of collective dynamics enters exclusively through the wavefunction fidelity $\mu \neq 0$7. This factorization is the structural reason quantum geometry ultimately controls the dispersion.

Generalized Hubbard–Stratonovich formalism

A key methodological contribution is the use of the time-sliced Generalized Hubbard–Stratonovich (GHS) transformation, following Kerman et al., which decouples the interaction simultaneously in the density (particle-hole) and pairing (particle-particle) channels via auxiliary fields $\mu \neq 0$8 and $\mu \neq 0$9, with gauge parameters $\omega \propto q^2$0 distributing the interaction strength between channels. The residual cross-channel term scales as $\omega \propto q^2$1 and vanishes in the continuous-time limit, avoiding the overcounting of microscopic interactions inherent in standard continuous HS decouplings applied to a single channel. After rescaling fields into the symmetric basis $\omega \propto q^2$2, the Gaussian fluctuation action takes the form $\omega \propto q^2$3, where the $\omega \propto q^2$4 response matrix is $\omega \propto q^2$5, with bubble elements given by Matsubara sums over Nambu transition matrix elements weighted by the fidelity $\omega \propto q^2$6. Time-reversal symmetry guarantees $\omega \propto q^2$7, i.e., equality of electron- and hole-sector fidelities, which underpins several cancellations below.

Analytical results: gapless pole, vanishing velocity, and metric-controlled curvature

Three analytic steps produce the main result. First, evaluating the static determinant using the gap equation yields $\omega \propto q^2$8, so $\omega \propto q^2$9 is the unique root: a gapless Goldstone boson exists despite full phase–amplitude–density hybridization. The null space of $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$0 is two-dimensional, spanned by a pure phase vector $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$1 and an amplitude–density hybrid vector $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$2; the latter shows that particle-hole asymmetry locks amplitude and density fluctuations at a fixed ratio $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$3, forming a zero-energy composite. The remaining massive amplitude–density sector acts as a virtual polarization background that is integrated out via a Schur complement, downfolding its geometric content into the soft-mode dynamics.

Second, expanding the secular equation to first order in $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$4 gives a skew-Hermitian block $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$5, while the momentum gradient of $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$6 vanishes identically because the derivative of the normalization identity $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$7 forces the linear term of the fidelity expansion to zero. The resulting algebraic condition forces $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$8: the linear acoustic term is absent, independent of spatial symmetries. This is a strong, parameter-free claim—it rules out any linear Goldstone velocity in any TRS isolated flat-band superfluid.

Third, at second order the projected secular equation $$\kappa = \frac{2E}{N_k}\sum_{\bm{k}} g_{qq}(\bm{k})$$9 with $E = \sqrt{\mu^2 + \Delta_0^2}$0 yields

$E = \sqrt{\mu^2 + \Delta_0^2}$1

where $E = \sqrt{\mu^2 + \Delta_0^2}$2 is the quantum metric projected along the propagation direction. A notable corollary is the cancellation of explicit doping dependence: although $E = \sqrt{\mu^2 + \Delta_0^2}$3 drives the hybridization, the combination $E = \sqrt{\mu^2 + \Delta_0^2}$4 absorbs it, so $E = \sqrt{\mu^2 + \Delta_0^2}$5 depends on filling only through the constraint that $E = \sqrt{\mu^2 + \Delta_0^2}$6 remains inside the flat band. The derivation also makes clear where generality ends: the factorization of electron and hole fidelities relies on TRS ($E = \sqrt{\mu^2 + \Delta_0^2}$7), so both $E = \sqrt{\mu^2 + \Delta_0^2}$8 and the metric formula are restricted to TRS-preserved systems; in Chern bands Berry curvature may revive a linear term or modify $E = \sqrt{\mu^2 + \Delta_0^2}$9.

Numerical verification on the Lieb lattice

The predictions are tested on the Lieb lattice with staggered hoppings, using $s$0, $s$1, $s$2, and $s$3. Solving the full three-channel secular equation $s$4 produces a gapless parabolic dispersion surface, and log-log scaling fits along $s$5 yield power-law exponents $s$6 across all values of $s$7, confirming the absence of the linear term. Fitted curvatures track the analytical prediction $s$8 as the staggered hopping tunes the average quantum metric, providing quantitative support for the geometric origin of the dispersion coefficient.

Limitations and open questions

Several restrictions should be stated plainly. The analysis assumes an energetically isolated flat band with $s$9 smaller than the band gap, so interband screening and hybridization with dispersive bands are excluded; the applicability to realistic materials such as magic-angle twisted bilayer graphene, where bands are nearly but not exactly flat and remote bands lie nearby, is not established here. Long-range Coulomb interactions, which would gap or shift the density channel and could alter the hybridization structure, are not included—the interaction is strictly local and attractive. The results are derived at $H = H_0 + H_{\text{int}}$0 within Gaussian fluctuation theory, and their extension to finite temperature holds only for $H = H_0 + H_{\text{int}}$1; non-Gaussian or strong-coupling corrections are unaddressed. Finally, the fate of the collective mode when TRS is broken—for instance whether Berry curvature generates a finite linear velocity or non-reciprocal dispersion—is explicitly left open by the authors.

Conclusion

This paper establishes that in any time-reversal-symmetric, isolated flat-band superfluid at finite doping, the Higgs, Goldstone, and density channels hybridize into a single gapless collective mode with quadratic dispersion whose curvature is fixed by the normal-state quantum metric. The multi-channel GHS framework, the Schur-complement downfolding, and the Lieb-lattice numerical verification together provide a controlled and quantitatively confirmed derivation, extending the known role of quantum geometry from the static superfluid weight to the dynamical excitation spectrum of flat-band superfluids.