Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2608.19044

Emergent BEC mechanism in flat-band superconductors

V. Berger, I. S. Tupitsyn, B. Currie, Boris V. Svistunov, Evgeny Kozik, Nikolay V. Prokof'ev

Published Aug 19, 2026• 0 citations
arXivHTMLPDF

Abstract

The formation of bound bosonic pairs of fermions, followed by their Bose-Einstein (quasi) condensation (BEC), is a foundational mechanism of superconductivity. At low filling, flat-band superconductivity is well captured by this mechanism provided the flat band is separated from the occupied lower band by an energy gap. However, particularly high $T_c$ values are anticipated when the non-interacting flat and lower bands touch---as in the prototypical attractive Lieb-lattice model studied here---invalidating the conventional picture: while interactions might protect the bound state by opening a gap, no small parameter guarantees the separation of the bound-state energy from this gap or occupied-band excitations, leaving the pair's fate uncertain. Based on a controlled-precision numerical protocol---which we demonstrate to be essential in this fundamentally non-perturbative problem---we show that the BEC mechanism, underpinned by an interaction-induced gap, is generically robust and remarkably efficient: fermions doped into the flat band form bound pairs within this gap with an anomalously light effective mass, enabling an exceptionally high $T_c$.

DOI
10.48550/arxiv.2608.19044

Related papers

Numerically Exact Study of Flat-Band Superconductivity

I. S. Tupitsyn, B. Currie, Boris V. Svistunov +2 more

Current theories of high-temperature superconductivity in flat-band systems predict a linear dependence of the transition temperature on the attractive interaction, $T_c(U) = c|U|$. However, the value of $c$ and the full nonlinear $T_c(U)$ curve---with a maximum at large $|U|$---are known beyond mean-field and quantum geometry estimates only for systems with isolated flat bands. Using a controlled diagrammatic Monte Carlo technique, we trace the onset of superfluid response in the Lieb lattice with attractive Hubbard interaction. Focusing on the half-filled flat-band case, where the ordering mechanism differs fundamentally from both BCS and BEC (preformed Cooper pair) scenarios, we find that the pairing response diverges linearly with decreasing temperature over a broad range of $U$, leading to a sharp crossover to long-range correlations at a characteristic temperature $T_*$, which provides a controlled upper bound on $T_c$. The highest $T_*$ occurs when all three bands touch at a single momentum point, potentially corresponding to high $T_c$ values.

arXiv:2604.05997•Apr 7, 2026•0 citations

Superconductivity versus bound-state formation in a two-band superconductor with small Fermi energy: Applications to Fe pnictides/chalcogenides and doped SrTiO3

Andrey V. Chubukov, Ilya Eremin, Dmitri V. Efremov

We analyze the interplay between superconductivity and the formation of bound pairs of fermions (BCS-BEC crossover) in a 2D model of interacting fermions with small Fermi energy ${E}_{F}$ and weak attractive interaction, which extends to energies well above ${E}_{F}$. The 2D case is special because a two-particle bound state forms at arbitrary weak interaction, and already at weak coupling, one has to distinguish between the bound-state formation and superconductivity. We briefly review the situation in the one-band model and then consider two different two-band models: one with one hole band and one electron band and another with two hole or two electron bands. In each case, we obtain the bound-state energy $2{E}_{0}$ for two fermions in a vacuum and solve the set of coupled equations for the pairing gaps and the chemical potentials to obtain the onset temperature of the pairing ${T}_{\mathrm{ins}}$ and the quasiparticle dispersion at $T=0$. We then compute the superfluid stiffness ${\ensuremath{\rho}}_{s}(T=0)$ and obtain the actual ${T}_{c}$. For definiteness, we set ${E}_{F}$ in one band to be near zero and consider different ratios of ${E}_{0}$ and ${E}_{F}$ in the other band. We show that at ${E}_{F}\ensuremath{\gg}{E}_{0}$, the behavior of both two-band models is BCS-like in the sense that ${T}_{c}\ensuremath{\approx}{T}_{\mathrm{ins}}\ensuremath{\ll}{E}_{F}$ and $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{T}_{c}$. At ${E}_{F}\ensuremath{\ll}{E}_{0}$, the two models behave differently: in the model with two hole/two electron bands, ${T}_{\mathrm{ins}}\ensuremath{\sim}{E}_{0}/ln\frac{{E}_{0}}{{E}_{F}}, \mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{({E}_{0}{E}_{F})}^{1/2}$, and ${T}_{c}\ensuremath{\sim}{E}_{F}$, like in the one-band model. In between ${T}_{\mathrm{ins}}$ and ${T}_{c}$, the system displays a preformed pair behavior. In the model with one hole and one electron bands, ${T}_{c}$ remains of order ${T}_{\mathrm{ins}}$, and both remain finite at ${E}_{F}=0$ and of the order of ${E}_{0}$. The preformed pair behavior still does exist in this model because ${T}_{c}$ is numerically smaller than ${T}_{\mathrm{ins}}$. For both models, we reexpress ${T}_{\mathrm{ins}}$ in terms of the fully renormalized two-particle scattering amplitude by extending to the two-band case (the method pioneered by Gorkov and Melik-Barkhudarov back in 1961). We apply our results for the model with a hole and an electron band to Fe pnictides and Fe chalcogenides in which a superconducting gap has been detected on the bands that do not cross the Fermi level, and to FeSe, in which the superconducting gap is comparable to the Fermi energy. We apply the results for the model with two electron bands to Nb-doped ${\mathrm{SrTiO}}_{3}$ and argue that our theory explains the rapid increase of ${T}_{c}$ when both bands start crossing the Fermi level.

arXiv:1601.01678•Physical review. B./Physical review. B•May 25, 2016•90 citations

Superconducting transitions in flat-band systems

Vladimir Iglovikov, F. Hébert, Benoît Grémaud +2 more

The physics of strongly correlated quantum particles within a flat band was originally explored as a route to itinerant ferromagnetism and, indeed, a celebrated theorem by Lieb rigorously establishes that the ground state of the repulsive Hubbard model on a bipartite lattice with an unequal number of sites in each sublattice must have nonzero spin $S$ at half filling. Recently, there has been interest in Lieb geometries due to the possibility of topological insulator, nematic, and Bose-Einstein condensed (BEC) phases. In this paper, we extend the understanding of the attractive Hubbard model on the Lieb lattice by using determinant quantum Monte Carlo to study real space charge and pair correlation functions not addressed by the Lieb theorems. Specifically, our results show unusual charge and charge transfer signatures within the flat band, and a reduction in pairing order at $\ensuremath{\rho}=2/3$ and $\ensuremath{\rho}=4/3$, the points at which the flat band is first occupied and then completely filled. We compare our results to the case of flat bands in the Kagome lattice and demonstrate that the behavior observed in the two cases is rather different.

arXiv:1404.5482•Physical Review B•Sep 11, 2014•111 citations

Cooper pair insulators and theory of correlated superconductors

Predrag Nikolić, Zlatko Tešanović

The pseudogap state of cuprate high-temperature superconductors has been often viewed as either a yet unknown competing order or a precursor state to superconductivity. While awaiting the resolution of the pseudogap problem in cuprates, we demonstrate that local pairing fluctuations, vortex liquid dynamics, and other precursor phenomena can emerge quite generally whenever fermionic excitations remain gapped across the superconducting transition, regardless of the gap origin. Our choice of a tractable model is a lattice band insulator with short-range attractive interactions between fermions in the $s$-wave channel. An effective crossover between Bardeen-Cooper-Schrieffer (BCS) and Bose-Einstein condensate (BEC) regimes can be identified in any band insulator above two dimensions, while in two dimensions only the BEC regime exists. The superconducting transition is ``unconventional'' (non-pair-breaking) in the BEC regime, identified by either the bosonic mean field or $\mathit{XY}$ universality class. The insulator adjacent to the superconductor in the BEC regime is a bosonic Mott insulator of Cooper pairs, which may be susceptible to charge density wave ordering. We construct a function of the many-body excitation spectrum whose nonanalytic changes define a sharp distinction between band and Mott insulators. The corresponding ``second-order transition'' can be observed out of equilibrium by driving a Cooper pair laser in the Mott insulator. We explicitly show that the gap for charged bosonic excitations lies below the threshold for Cooper pair breakup in any BEC regime, despite quantum fluctuations. Our discussion ends with a view of possible consequences for cuprates, where antinodal pair dynamics has certain features in common with our simple $s$-wave picture.

arXiv:1008.4369•Physical Review B•Feb 2, 2011•16 citations

Designer flat bands: Topology and enhancement of superconductivity

Si Min Chan, Benoît Grémaud, G. G. Batrouni

We construct quasi-one-dimensional topological and nontopological three-band lattices with a tunable band gap and winding number of the flat band. Using full multiband mean field (MF) and exact density matrix renormalization group (DMRG) calculations, we show explicitly how the band gap affects pairing and superconductivity (SC) in a Hubbard model with attractive interactions. We obtain excellent agreement between MF and DMRG calculations. The SC weight ${D}_{s}$ on the gapped topological, $W\ensuremath{\ne}0$, flat band increases linearly with interaction strength $U$ for low values and with a slope that depends on the details of the compact localized state at $U=0$. As $U\ensuremath{\rightarrow}0$ for the gapped nontopological flat band ($W=0$), ${D}_{s}$ decays with a power law faster than quadratically but slower than exponentially. In the gapless case (flat band touching the band above it), we find at low $U$ (for both $W=0$ and $W\ensuremath{\ne}0$) that ${D}_{s}\ensuremath{\propto}{U}^{\ensuremath{\varphi}}$, with $\ensuremath{\varphi}<1$. In other words, ${D}_{s}$ increases faster than linearly for low $U$, thus favoring SC at weak interaction more than the gapped case. Our results reestablish that the BCS mean field and quantum metric alone are insufficient to characterize SC at weak coupling.

arXiv:2206.10651•Physical review. B./Physical review. B•Sep 22, 2022•25 citations

Band-edge BCS–BEC crossover in a two-band superconductor: physical properties and detection parameters

Andrea Guidini, Andrea Perali

Superconductivity in iron-based, magnesium diborides, and other novel superconducting materials has a strong multi-band and multi-gap character. Recent experiments support the possibillity for a BCS-BEC crossover induced by strong-coupling and proximity of the chemical potential to the band edge of one of the bands. Here we study the simplest theoretical model which accounts for the BCS-BEC crossover in a two-band superconductor, considering tunable interactions and tunable energy separations between the bands. Mean-field results for condensate fraction, correlation length, and superconducting gap are reported in typical crossover diagrams to locate the boundaries of the BCS, crossover, and BEC regimes. When the superconducting gap is of the order of the local chemical potential, superconductivity is in the crossover regime of the BCS-BEC crossover and the Fermi surface of the small band is smeared by the gap opening. In this situation, small and large Cooper pairs coexist in the total condensate, which is the optimal condition for high-Tc superconductivity. The ratio between the gap and the Fermi energy in a given band results to be the best detection parameter for experiments to locate the system in the BCS-BEC crossover. Using available experimental data, our analysis shows that iron-based superconductors have the partial condensate of the small Fermi surface in the crossover regime of the BCS-BEC crossover, supporting the recent ARPES findings.

arXiv:1407.3109•Superconductor Science and Technology•Oct 31, 2014•66 citations
The formation of bound bosonic pairs of fermions, followed by their Bose-Einstein (quasi) condensation (BEC), is a foundational mechanism of superconductivity. At low filling, flat-band superconductivity is well captured by this mechanism provided the flat band is separated from the occupied lower band by an energy gap. However, particularly high TcT_cTc​ values are anticipated when the non-interacting flat and lower bands touch---as in the prototypical attractive Lieb-lattice model studied here---invalidating the conventional picture: while interactions might protect the bound state by opening a gap, no small parameter guarantees the separation of the bound-state energy from this gap or occupied-band excitations, leaving the pair's fate uncertain. Based on a controlled-precision numerical protocol---which we demonstrate to be essential in this fundamentally non-perturbative problem---we show that the BEC mechanism, underpinned by an interaction-induced gap, is generically robust and remarkably efficient: fermions doped into the flat band form bound pairs within this gap with an anomalously light effective mass, enabling an exceptionally high TcT_cTc​.
Current theories of high-temperature superconductivity in flat-band systems predict a linear dependence of the transition temperature on the attractive interaction, Tc(U)=c∣U∣T_c(U) = c|U|Tc​(U)=c∣U∣. However, the value of ccc and the full nonlinear Tc(U)T_c(U)Tc​(U) curve---with a maximum at large ∣U∣|U|∣U∣---are known beyond mean-field and quantum geometry estimates only for systems with isolated flat bands. Using a controlled diagrammatic Monte Carlo technique, we trace the onset of superfluid response in the Lieb lattice with attractive Hubbard interaction. Focusing on the half-filled flat-band case, where the ordering mechanism differs fundamentally from both BCS and BEC (preformed Cooper pair) scenarios, we find that the pairing response diverges linearly with decreasing temperature over a broad range of UUU, leading to a sharp crossover to long-range correlations at a characteristic temperature T∗T_*T∗​, which provides a controlled upper bound on TcT_cTc​. The highest T∗T_*T∗​ occurs when all three bands touch at a single momentum point, potentially corresponding to high TcT_cTc​ values.
We analyze the interplay between superconductivity and the formation of bound pairs of fermions (BCS-BEC crossover) in a 2D model of interacting fermions with small Fermi energy EF{E}_{F}EF​ and weak attractive interaction, which extends to energies well above EF{E}_{F}EF​. The 2D case is special because a two-particle bound state forms at arbitrary weak interaction, and already at weak coupling, one has to distinguish between the bound-state formation and superconductivity. We briefly review the situation in the one-band model and then consider two different two-band models: one with one hole band and one electron band and another with two hole or two electron bands. In each case, we obtain the bound-state energy 2E02{E}_{0}2E0​ for two fermions in a vacuum and solve the set of coupled equations for the pairing gaps and the chemical potentials to obtain the onset temperature of the pairing Tins{T}_{\mathrm{ins}}Tins​ and the quasiparticle dispersion at T=0T=0T=0. We then compute the superfluid stiffness ${\ensuremath{\rho}}_{s}(T=0)$ and obtain the actual Tc{T}_{c}Tc​. For definiteness, we set EF{E}_{F}EF​ in one band to be near zero and consider different ratios of E0{E}_{0}E0​ and EF{E}_{F}EF​ in the other band. We show that at ${E}_{F}\ensuremath{\gg}{E}_{0}$, the behavior of both two-band models is BCS-like in the sense that ${T}_{c}\ensuremath{\approx}{T}_{\mathrm{ins}}\ensuremath{\ll}{E}_{F}$ and $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{T}_{c}$. At ${E}_{F}\ensuremath{\ll}{E}_{0}$, the two models behave differently: in the model with two hole/two electron bands, ${T}_{\mathrm{ins}}\ensuremath{\sim}{E}_{0}/ln\frac{{E}_{0}}{{E}_{F}}, \mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{({E}_{0}{E}_{F})}^{1/2}$, and ${T}_{c}\ensuremath{\sim}{E}_{F}$, like in the one-band model. In between Tins{T}_{\mathrm{ins}}Tins​ and Tc{T}_{c}Tc​, the system displays a preformed pair behavior. In the model with one hole and one electron bands, Tc{T}_{c}Tc​ remains of order Tins{T}_{\mathrm{ins}}Tins​, and both remain finite at EF=0{E}_{F}=0EF​=0 and of the order of E0{E}_{0}E0​. The preformed pair behavior still does exist in this model because Tc{T}_{c}Tc​ is numerically smaller than Tins{T}_{\mathrm{ins}}Tins​. For both models, we reexpress Tins{T}_{\mathrm{ins}}Tins​ in terms of the fully renormalized two-particle scattering amplitude by extending to the two-band case (the method pioneered by Gorkov and Melik-Barkhudarov back in 1961). We apply our results for the model with a hole and an electron band to Fe pnictides and Fe chalcogenides in which a superconducting gap has been detected on the bands that do not cross the Fermi level, and to FeSe, in which the superconducting gap is comparable to the Fermi energy. We apply the results for the model with two electron bands to Nb-doped SrTiO3{\mathrm{SrTiO}}_{3}SrTiO3​ and argue that our theory explains the rapid increase of Tc{T}_{c}Tc​ when both bands start crossing the Fermi level.
The physics of strongly correlated quantum particles within a flat band was originally explored as a route to itinerant ferromagnetism and, indeed, a celebrated theorem by Lieb rigorously establishes that the ground state of the repulsive Hubbard model on a bipartite lattice with an unequal number of sites in each sublattice must have nonzero spin SSS at half filling. Recently, there has been interest in Lieb geometries due to the possibility of topological insulator, nematic, and Bose-Einstein condensed (BEC) phases. In this paper, we extend the understanding of the attractive Hubbard model on the Lieb lattice by using determinant quantum Monte Carlo to study real space charge and pair correlation functions not addressed by the Lieb theorems. Specifically, our results show unusual charge and charge transfer signatures within the flat band, and a reduction in pairing order at $\ensuremath{\rho}=2/3$ and $\ensuremath{\rho}=4/3$, the points at which the flat band is first occupied and then completely filled. We compare our results to the case of flat bands in the Kagome lattice and demonstrate that the behavior observed in the two cases is rather different.
The pseudogap state of cuprate high-temperature superconductors has been often viewed as either a yet unknown competing order or a precursor state to superconductivity. While awaiting the resolution of the pseudogap problem in cuprates, we demonstrate that local pairing fluctuations, vortex liquid dynamics, and other precursor phenomena can emerge quite generally whenever fermionic excitations remain gapped across the superconducting transition, regardless of the gap origin. Our choice of a tractable model is a lattice band insulator with short-range attractive interactions between fermions in the sss-wave channel. An effective crossover between Bardeen-Cooper-Schrieffer (BCS) and Bose-Einstein condensate (BEC) regimes can be identified in any band insulator above two dimensions, while in two dimensions only the BEC regime exists. The superconducting transition is ``unconventional'' (non-pair-breaking) in the BEC regime, identified by either the bosonic mean field or XY\mathit{XY}XY universality class. The insulator adjacent to the superconductor in the BEC regime is a bosonic Mott insulator of Cooper pairs, which may be susceptible to charge density wave ordering. We construct a function of the many-body excitation spectrum whose nonanalytic changes define a sharp distinction between band and Mott insulators. The corresponding ``second-order transition'' can be observed out of equilibrium by driving a Cooper pair laser in the Mott insulator. We explicitly show that the gap for charged bosonic excitations lies below the threshold for Cooper pair breakup in any BEC regime, despite quantum fluctuations. Our discussion ends with a view of possible consequences for cuprates, where antinodal pair dynamics has certain features in common with our simple sss-wave picture.
We construct quasi-one-dimensional topological and nontopological three-band lattices with a tunable band gap and winding number of the flat band. Using full multiband mean field (MF) and exact density matrix renormalization group (DMRG) calculations, we show explicitly how the band gap affects pairing and superconductivity (SC) in a Hubbard model with attractive interactions. We obtain excellent agreement between MF and DMRG calculations. The SC weight Ds{D}_{s}Ds​ on the gapped topological, $W\ensuremath{\ne}0$, flat band increases linearly with interaction strength UUU for low values and with a slope that depends on the details of the compact localized state at U=0U=0U=0. As $U\ensuremath{\rightarrow}0$ for the gapped nontopological flat band (W=0W=0W=0), Ds{D}_{s}Ds​ decays with a power law faster than quadratically but slower than exponentially. In the gapless case (flat band touching the band above it), we find at low UUU (for both W=0W=0W=0 and $W\ensuremath{\ne}0$) that ${D}_{s}\ensuremath{\propto}{U}^{\ensuremath{\varphi}}$, with $\ensuremath{\varphi}<1$. In other words, Ds{D}_{s}Ds​ increases faster than linearly for low UUU, thus favoring SC at weak interaction more than the gapped case. Our results reestablish that the BCS mean field and quantum metric alone are insufficient to characterize SC at weak coupling.
Superconductivity in iron-based, magnesium diborides, and other novel superconducting materials has a strong multi-band and multi-gap character. Recent experiments support the possibillity for a BCS-BEC crossover induced by strong-coupling and proximity of the chemical potential to the band edge of one of the bands. Here we study the simplest theoretical model which accounts for the BCS-BEC crossover in a two-band superconductor, considering tunable interactions and tunable energy separations between the bands. Mean-field results for condensate fraction, correlation length, and superconducting gap are reported in typical crossover diagrams to locate the boundaries of the BCS, crossover, and BEC regimes. When the superconducting gap is of the order of the local chemical potential, superconductivity is in the crossover regime of the BCS-BEC crossover and the Fermi surface of the small band is smeared by the gap opening. In this situation, small and large Cooper pairs coexist in the total condensate, which is the optimal condition for high-Tc superconductivity. The ratio between the gap and the Fermi energy in a given band results to be the best detection parameter for experiments to locate the system in the BCS-BEC crossover. Using available experimental data, our analysis shows that iron-based superconductors have the partial condensate of the small Fermi surface in the crossover regime of the BCS-BEC crossover, supporting the recent ARPES findings.