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:2209.00007 · Communications Physics

Many-body superconductivity in topological flat bands

Jonah Herzog-Arbeitman, Aaron Chew, Kukka-Emilia Huhtinen, Päivi Törmä, Bogdan Andrei Bernevig

Published Jul 4, 2026• 9 citations
arXivHTMLPDF

Abstract

In a flat band superconductor, bosonic excitations disperse while unpaired electrons are immobile. To study this strongly interacting system, we construct a family of multi-band Hubbard models with exact eta-pairing ground states in all space groups. We analytically compute their many-body excitations and find that the Cooper pair bound states obey an effective single-boson hopping Hamiltonian written in terms of the non-interacting wavefunctions. These bound states possess a unique zero-energy excitation whose quadratic dispersion is determined by the minimal quantum metric. The rest of the bound state spectrum is classified by topological quantum chemistry, which we use to identify Cooper pairs with Weyl nodes, higher angular momentum pairing, and fragile topology. Broadening the band dispersion, we find the strongest pairing occurs at half filling and not at the highest density of states. Our analytical results prove that the elementary fermionic and bosonic excitations, corresponding respectively to pairing and coherence, are governed by distinct properties of the flat band, rigorously establishing its unique phenomenology. Superconductivity in flat bands is a central feature of quantum geometric effects, supported by mean-field theory and numerics. Here, the authors exhibit a large family of flat band Hubbard models whose ground states, Bogoliubov-de-Gennes particles, and Cooper pairs are exactly solvable with analytically calculable dispersions.

DOI
10.1038/s42005-026-02732-2

Related papers

Obstructed Cooper pairs in flat band systems -- weakly-coherent superfluids and exact spin liquids

Tamaghna Hazra, Verma, Nishchhal, Schmalian, Jörg

Superconductivity in a partially filled flat band presents a vexing conceptual hurdle because the absence of a Fermi surface precludes a weak-coupling regime where one can extend insights from the Bardeen-Cooper-Schrieffer picture of a Fermi surface instability. We approach the strongly correlated problem of flat band superconductivity from the strong coupling limit of local attractive interactions on line-graph lattices, whose non-interacting bandstructures host exactly flat bands due to frustrated hopping. In this limit, the pair kinetic energy which sets the superfluid stiffness is expected to scale inversely with the pair binding interaction. Here we demonstrate a striking counterexample. We show that when doped charges propagate on the line-graph of a lattice with strong pairing interaction and broken time-reversal symmetry, they bind into obstructed Cooper pairs whose motion is frustrated by destructive interference. As a result, the leading-order pair kinetic energy vanishes identically in the strong-coupling expansion, producing a flat bosonic band of compact localised pair states, zero superfluid stiffness at leading order, and an extensively degenerate many-body ground state manifold. At quarter filling, the frustrated pair dynamics maps onto a quantum dimer model which has a $d$-wave resonating-valence-bond ground state when time-reversal is broken. The pairing Hamiltonian in this limit thus has a topologically ordered spin liquid ground state which becomes exact at the analytically solvable Rokhsar-Kivelson point with long-range entanglement and deconfined holon excitations. Interestingly, we find exact compact localised eigenstates and extensive degeneracies in the many-body eigenstates of this emergent dimer model. Our results establish a disorder-free mechanism for interaction-driven localisation, in which strong pairing collapses the kinetic energy of Cooper pairs.

arXiv:2411.17815•Nov 26, 2024•0 citations

Superconductivity in geometrically and topologically nontrivial lattice models

Aleksi Julku

In this thesis, superconductivity of fermions in two-dimensional Hubbard lattice models is theoretically studied. Each of the studied lattice structures features intriguing properties that provide valuable insight to understand superconductivity in general. While real condensed matter systems are too complicated to tackle with theoretical tools, the Hubbard model is one of the simplest quantum many-body models to possess a wide spectrum of physically relevant phenomena, from superconductivity to magnetism and topological order. Hubbard models are also interesting as they can be experimentally realized and accurately controlled by utilizing ultracold gases. The first part of the thesis explains the required theoretical background of superconductivity in the level of mean-field and linear response theories. Specifically, we primarily focus on the superfluid weight which implies the existence of the Meissner effect and dissipationless current, defining properties of superconductors. We further discuss how the Berezinskii-Kosterlitz-Thouless (BKT) transition, related to the loss of superconductivity due to the thermal fluctuations of the superconducting order parameter in high enough temperatures, can be accessed via the superfluid weight. In the second part, the main results are discussed. We show how the superconductivity of lattice models featuring dispersionless Bloch bands, so-called flat bands, can be large and originating from the quantum geometric properties of the Bloch states of the flat bands. This is in contrast to conventional dispersive bands where superconductivity is determined by the effective mass of the particles, i.e. by the dispersion relation of the Bloch bands alone and not by the Bloch states. This is well demonstrated in the Lieb lattice geometry, studied in Publication I, which features both a strictly flat and two dispersive Bloch bands, allowing direct comparison of superconductivity in dispersive and flat bands. We show how the superconductivity of the flat band is more robust than that of the dispersive bands and how the superfluid weight behaves non-monotonically as a function of the interaction strength. This non-monotonic behavior should be experimentally accessible. Flat band superconductivity is also the topic of Publication III, where the experimentally realized twisted bilayer graphene (TBG) is studied with local and non-local attractive interactions. We predict qualitative differences between local and non-local interaction schemes which could be distinguished experimentally and reveal the quantum geometric origin of superconductivity in TBG systems in case of both the assumed interaction schemes. Exotic Fulde-Ferrell (FF) superfluid phases, characterized by the non-zero momentum of Cooper pairs, were explored in Publication II. We show how the FF phases can be stabilized against thermal fluctuations in the presence of spin-orbit coupling (SOC) and that topologically non-trivial FF states, created with SOC and Zeeman fields, can be found at finite temperatures.; Tassa…

•Jan 1, 2019•0 citations

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

Pairing and non-Fermi liquid behavior in partially flat-band systems: Beyond nesting physics

Sharareh Sayyad, Edwin W. Huang, Motoharu Kitatani +4 more

While many-body effects in flat-band systems are receiving renewed hot interest in condensed-matter physics for superconducting and topological properties as well as for magnetism, studies have primarily been restricted to multiband systems (with coexisting flat and dispersive bands). Here we focus on one-band systems where a band is ``partially flat,'' comprising flat and dispersive portions in $k$ space to reveal whether intriguing correlation effects can already arise on the simplest possible one-band level. For that, the two-dimensional repulsive Hubbard model is studied for two models having different flat areas in an intermediate-coupling regime with the dynamical mean-field theory combined with the fluctuation exchange approximation. We have a crossover from ferromagnetic to antiferromagnetic spin fluctuations as the band filling is varied and this triggers, for the model with a wider flat portion, a triplet-pair superconductivity favored over an unusually wide filling region, which is taken over by a sharply growing singlet pairing. For the model with a narrower flat portion, ${T}_{C}$ against filling exhibits an unusual double-peaked ${T}_{C}$ dome, associated with different numbers of nodes in the gap function having remarkably extended pairs in real space. We identify these as a manifestation of the physics outside the conventional nesting physics where only the pair scattering across the Fermi surface in designated (hot) spots is relevant. Another correlation effect arising from the flattened band is found in a non-Fermi-liquid behavior as detected in the momentum distribution function, frequency dependence of the self-energy, and spectral function. These indicate that unusual correlation physics can indeed occur in flat-band systems.

arXiv:1903.09888•Physical review. B./Physical review. B•Jan 6, 2020•35 citations

Superconductivity in repulsively interacting fermions on a diamond chain: Flat-band-induced pairing

Keita Kobayashi, Masahiko Okumura, Susumu Yamada +2 more

To explore whether a flat-band system can accommodate superconductivity, we consider repulsively interacting fermions on the diamond chain, a simplest possible quasi-one-dimensional system that contains a flat band. Exact diagonalization and the density-matrix renormalization group are used to show that we have a significant binding energy of a Cooper pair with a long-tailed pair-pair correlation in real space when the total band filling is slightly below 1/3, where a filled dispersive band interacts with the flat band that is empty but close to ${E}_{F}$. Pairs selectively formed across the outer sites of the diamond chain are responsible for the pairing correlation. At exactly 1/3-filling an insulating phase emerges, where the entanglement spectrum indicates the particles on the outer sites are highly entangled and topological. These come from a peculiarity of the flat band in which ``Wannier orbits'' are not orthogonalizable.

arXiv:1608.00125•Physical review. B./Physical review. B•Dec 1, 2016•77 citations

Pair size and quantum geometry in a multiband Hubbard model

M. Iskin

We study the size of two-body bound states and Cooper pairs within a multiband Hubbard model that features time-reversal symmetry and uniform pairing on a generic lattice. Our analysis involves (i) an exact calculation of the localization tensor to determine the size of lowest-lying two-body bound state in vacuum, and (ii) an evaluation of the analogous tensor to estimate the average size of Cooper pairs within the mean-field Bardeen-Cooper-Schrieffer to Bose-Einstein condensate (BCS-BEC) crossover theory at zero temperature. Beyond the conventional intraband contribution that depends on Bloch bands, we show that pair size also has a geometric contribution governed by the quantum-metric tensor of the Bloch states and their band-resolved quantum-metric tensors. As a concrete example, we investigate the pyrochlore-Hubbard model numerically and demonstrate that, while the pair size diverges in the weakly interacting BCS regime of dispersive bands, it remains finite and relatively small in the flat-band regime, even for infinitesimal interaction, perfectly matching the exact two-body result in the dilute limit.

arXiv:2409.14921•Physical review. B./Physical review. B•Jan 3, 2025•6 citations
In a flat band superconductor, bosonic excitations disperse while unpaired electrons are immobile. To study this strongly interacting system, we construct a family of multi-band Hubbard models with exact eta-pairing ground states in all space groups. We analytically compute their many-body excitations and find that the Cooper pair bound states obey an effective single-boson hopping Hamiltonian written in terms of the non-interacting wavefunctions. These bound states possess a unique zero-energy excitation whose quadratic dispersion is determined by the minimal quantum metric. The rest of the bound state spectrum is classified by topological quantum chemistry, which we use to identify Cooper pairs with Weyl nodes, higher angular momentum pairing, and fragile topology. Broadening the band dispersion, we find the strongest pairing occurs at half filling and not at the highest density of states. Our analytical results prove that the elementary fermionic and bosonic excitations, corresponding respectively to pairing and coherence, are governed by distinct properties of the flat band, rigorously establishing its unique phenomenology. Superconductivity in flat bands is a central feature of quantum geometric effects, supported by mean-field theory and numerics. Here, the authors exhibit a large family of flat band Hubbard models whose ground states, Bogoliubov-de-Gennes particles, and Cooper pairs are exactly solvable with analytically calculable dispersions.
Superconductivity in a partially filled flat band presents a vexing conceptual hurdle because the absence of a Fermi surface precludes a weak-coupling regime where one can extend insights from the Bardeen-Cooper-Schrieffer picture of a Fermi surface instability. We approach the strongly correlated problem of flat band superconductivity from the strong coupling limit of local attractive interactions on line-graph lattices, whose non-interacting bandstructures host exactly flat bands due to frustrated hopping. In this limit, the pair kinetic energy which sets the superfluid stiffness is expected to scale inversely with the pair binding interaction. Here we demonstrate a striking counterexample. We show that when doped charges propagate on the line-graph of a lattice with strong pairing interaction and broken time-reversal symmetry, they bind into obstructed Cooper pairs whose motion is frustrated by destructive interference. As a result, the leading-order pair kinetic energy vanishes identically in the strong-coupling expansion, producing a flat bosonic band of compact localised pair states, zero superfluid stiffness at leading order, and an extensively degenerate many-body ground state manifold. At quarter filling, the frustrated pair dynamics maps onto a quantum dimer model which has a ddd-wave resonating-valence-bond ground state when time-reversal is broken. The pairing Hamiltonian in this limit thus has a topologically ordered spin liquid ground state which becomes exact at the analytically solvable Rokhsar-Kivelson point with long-range entanglement and deconfined holon excitations. Interestingly, we find exact compact localised eigenstates and extensive degeneracies in the many-body eigenstates of this emergent dimer model. Our results establish a disorder-free mechanism for interaction-driven localisation, in which strong pairing collapses the kinetic energy of Cooper pairs.
In this thesis, superconductivity of fermions in two-dimensional Hubbard lattice models is theoretically studied. Each of the studied lattice structures features intriguing properties that provide valuable insight to understand superconductivity in general. While real condensed matter systems are too complicated to tackle with theoretical tools, the Hubbard model is one of the simplest quantum many-body models to possess a wide spectrum of physically relevant phenomena, from superconductivity to magnetism and topological order. Hubbard models are also interesting as they can be experimentally realized and accurately controlled by utilizing ultracold gases. The first part of the thesis explains the required theoretical background of superconductivity in the level of mean-field and linear response theories. Specifically, we primarily focus on the superfluid weight which implies the existence of the Meissner effect and dissipationless current, defining properties of superconductors. We further discuss how the Berezinskii-Kosterlitz-Thouless (BKT) transition, related to the loss of superconductivity due to the thermal fluctuations of the superconducting order parameter in high enough temperatures, can be accessed via the superfluid weight. In the second part, the main results are discussed. We show how the superconductivity of lattice models featuring dispersionless Bloch bands, so-called flat bands, can be large and originating from the quantum geometric properties of the Bloch states of the flat bands. This is in contrast to conventional dispersive bands where superconductivity is determined by the effective mass of the particles, i.e. by the dispersion relation of the Bloch bands alone and not by the Bloch states. This is well demonstrated in the Lieb lattice geometry, studied in Publication I, which features both a strictly flat and two dispersive Bloch bands, allowing direct comparison of superconductivity in dispersive and flat bands. We show how the superconductivity of the flat band is more robust than that of the dispersive bands and how the superfluid weight behaves non-monotonically as a function of the interaction strength. This non-monotonic behavior should be experimentally accessible. Flat band superconductivity is also the topic of Publication III, where the experimentally realized twisted bilayer graphene (TBG) is studied with local and non-local attractive interactions. We predict qualitative differences between local and non-local interaction schemes which could be distinguished experimentally and reveal the quantum geometric origin of superconductivity in TBG systems in case of both the assumed interaction schemes. Exotic Fulde-Ferrell (FF) superfluid phases, characterized by the non-zero momentum of Cooper pairs, were explored in Publication II. We show how the FF phases can be stabilized against thermal fluctuations in the presence of spin-orbit coupling (SOC) and that topologically non-trivial FF states, created with SOC and Zeeman fields, can be found at finite temperatures.; Tassa…
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.
While many-body effects in flat-band systems are receiving renewed hot interest in condensed-matter physics for superconducting and topological properties as well as for magnetism, studies have primarily been restricted to multiband systems (with coexisting flat and dispersive bands). Here we focus on one-band systems where a band is ``partially flat,'' comprising flat and dispersive portions in kkk space to reveal whether intriguing correlation effects can already arise on the simplest possible one-band level. For that, the two-dimensional repulsive Hubbard model is studied for two models having different flat areas in an intermediate-coupling regime with the dynamical mean-field theory combined with the fluctuation exchange approximation. We have a crossover from ferromagnetic to antiferromagnetic spin fluctuations as the band filling is varied and this triggers, for the model with a wider flat portion, a triplet-pair superconductivity favored over an unusually wide filling region, which is taken over by a sharply growing singlet pairing. For the model with a narrower flat portion, TC{T}_{C}TC​ against filling exhibits an unusual double-peaked TC{T}_{C}TC​ dome, associated with different numbers of nodes in the gap function having remarkably extended pairs in real space. We identify these as a manifestation of the physics outside the conventional nesting physics where only the pair scattering across the Fermi surface in designated (hot) spots is relevant. Another correlation effect arising from the flattened band is found in a non-Fermi-liquid behavior as detected in the momentum distribution function, frequency dependence of the self-energy, and spectral function. These indicate that unusual correlation physics can indeed occur in flat-band systems.
To explore whether a flat-band system can accommodate superconductivity, we consider repulsively interacting fermions on the diamond chain, a simplest possible quasi-one-dimensional system that contains a flat band. Exact diagonalization and the density-matrix renormalization group are used to show that we have a significant binding energy of a Cooper pair with a long-tailed pair-pair correlation in real space when the total band filling is slightly below 1/3, where a filled dispersive band interacts with the flat band that is empty but close to EF{E}_{F}EF​. Pairs selectively formed across the outer sites of the diamond chain are responsible for the pairing correlation. At exactly 1/3-filling an insulating phase emerges, where the entanglement spectrum indicates the particles on the outer sites are highly entangled and topological. These come from a peculiarity of the flat band in which ``Wannier orbits'' are not orthogonalizable.
We study the size of two-body bound states and Cooper pairs within a multiband Hubbard model that features time-reversal symmetry and uniform pairing on a generic lattice. Our analysis involves (i) an exact calculation of the localization tensor to determine the size of lowest-lying two-body bound state in vacuum, and (ii) an evaluation of the analogous tensor to estimate the average size of Cooper pairs within the mean-field Bardeen-Cooper-Schrieffer to Bose-Einstein condensate (BCS-BEC) crossover theory at zero temperature. Beyond the conventional intraband contribution that depends on Bloch bands, we show that pair size also has a geometric contribution governed by the quantum-metric tensor of the Bloch states and their band-resolved quantum-metric tensors. As a concrete example, we investigate the pyrochlore-Hubbard model numerically and demonstrate that, while the pair size diverges in the weakly interacting BCS regime of dispersive bands, it remains finite and relatively small in the flat-band regime, even for infinitesimal interaction, perfectly matching the exact two-body result in the dilute limit.