## Flat-band FFLO State from Quantum Geometric Discrepancy

Zi-Ting Sun, 1, ∗ Ruo-Peng Yu, 1, ∗ Shuai A. Chen, 1 Jin-Xin Hu, 1, † and K. T. Law 1, ‡

1 Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong, China (Dated: January 3, 2025)

We propose a new scheme for the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) Cooper pairing states within ultraflat bands. Central to our approach is the concept of 'quantum geometric discrepancy' (QGD), which characterizes the discrepancy in the quantum geometry of paired electrons, giving rise to the flat-band FFLO instability. Remarkably, we find that this instability is directly related to a new quantum geometric quantity we term 'anomalous quantum distance', which formally captures QGD. To model both QGD and the anomalous quantum distance, we examine a flat-band electronic Hamiltonian with tunable spin-dependent quantum metrics. Utilizing the band-projection method, we analyze the QGD-induced FFLO instability from pair susceptibility near the superconducting critical temperature. Furthermore, we perform self-consistent mean-field calculations to obtain the phase diagram of the BCS-FFLO transition driven by QGD, which aligns well with our analytical results. We emphasize that QGD serves as a distinctive protocol for stabilizing the flat-band FFLO phase.

Introduction. -Understanding unconventional superconductivity is a central theme in modern condensed matter physics, with important implications for both theoretical and experimental aspects [1-4]. Among the various mechanisms that extend beyond the conventional Bardeen-Cooper-Schrieffer (BCS) theory of pairing [5], the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase is notable for its unique characterization of Cooper pairs with finite center-of-mass momentum Q [6-9]. Various mechanisms have been explored to realize these exotic FFLO states in systems characterized by dispersive bands [10-19].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-3af06a44-fe2b-4f16-a649-615e21e44b5e.jpg)

In the original formulation of the FFLO phase, a spinsinglet pairing state with finite center-of-mass momentum Q ( | Q | ∝ B/v F ) is favored in the presence of a large Zeeman field B [see Fig. 1(a)]. However, the applicability of this microscopic mechanism becomes problematic in flat-band systems, where the Fermi velocity v F is nearly zero. This results in an unphysical divergence of Q as v F → 0, in the presence of a finite magnetic field. Consequently, two critical questions arise: First, what mechanism can stabilize the FFLO state as v F → 0? Second, what determines the center-of-mass momentum Q of Cooper pairs in the flat-band FFLO state? In dispersive bands, this is typically governed by Fermi surface nesting; however, in flat bands, the poorly defined Fermi surface weakens this relationship.

Lying in the context of the recent studies of the pronounced quantum geometric effects in flat-band systems [20-43], in this work, we address the aforementioned questions by linking the FFLO instability to the quantum geometry of Bloch wave functions. For concreteness, we propose a new scheme for FFLO pairing states within an isolated flat band, driven by 'quantum geometric discrepancy (QGD)'. This concept refers to the discrepancy in the quantum geometry of the wave functions of the two electrons forming a Cooper pair, which explicitly breaks time reversal ( T ) symmetry. Based

̸

FIG. 1. Schematic illustrations of (a) conventional FFLO from Zeeman effect, and (b) flat-band FFLO from QGD. In (a), an external B field induces spin-split Fermi surfaces, leading to the formation of Cooper pairs with finite momentum Q , where | Q | ∝ B/v F . In (b), the finite momentum Q is stabilized when the quantum metrics of the paired electrons are different ( G ↑ = G ↓ ) within flat bands.

on the band-projection formalism, we analyse the static pair susceptibility and establish that the critical temperature for finite-momentum pairing, which leads to FFLO instability, is directly associated with the band average of a quantum geometric quantity [as described in Eq. (6)]-'anomalous quantum distance (AQD)' ¯ d k ′ , k , formally defined as

$$\bar { d } _ { k ^ { \prime }, k } ^ { 2 } = 1 - \left | \langle u _ { k ^ { \prime } \uparrow } | u _ { - k \downarrow } ^ { * } \rangle \right | ^ { 2 }. \text{ \quad \ \ } ( 1 )$$

̸

Here, u k σ represents the Bloch wave function at momentum k and spin σ . AQD quantifies the overlap between the wave functions of paired electrons. In the presence of T symmetry, where u k ↑ = u ∗ - ↓ k , the quantity ¯ d k k , ′ becomes the well-known Hilbert-Schmidt quantum distance [44-47], which measures the quantum mechanical distance between two Bloch states. In this scenario, the BCS (zero-momentum) pairing state is favored. In the absence of T symmetry, u k ↑ = u ∗ - ↓ k generally makes this relation unestablished, and Eq. (1) will exhibit anomalous behaviors due to the QGD, enabling a

finite-momentum pairing state.

̸

As we will demonstrate in the following sections, by examining a flat-band model with tunable spin-dependent quantum metrics, we find that even in the flat-band limit, the finite momentum Q can be stabilized when pairing is frustrated by QGD, specifically when G ↑ = G ↓ [see Fig. 1(b)], where G ↑ and G ↓ denote the quantum metrics of the Bloch functions of the paired electrons. Our findings provide a viable route toward realizing flat-band FFLO states.

Effective action from band projection. -To investigate the FFLO instability on the flat band, we first derive the static superconducting Gaussian fluctuations using the band-projection formalism [38, 48, 49], and this analysis is rather general when considering conventional pairing. Without loss of generality, we consider the local Hubbard interaction H int = -U ∑ iα ˆ c † iα ↑ ˆ c † iα ↓ ˆ c iα ↓ ˆ c iα ↑ , where U &gt; 0, α and i denote the orbital and site index.

In order to analyze the low-energy behavior of superconducting pairings on the relevant isolated band using field theory, after the Fourier transformation (ˆ c iασ = 1 / √ N c ∑ k e i k r · i ˆ c α, k σ , where N c is the number of unit cells) and the uniform pairing approximation [30, 31], we project the interaction onto the target flat band by ˆ c α, k σ → u ∗ k σ ( α a )ˆ k σ , where ˆ c α, k σ and ˆ a k σ denotes the orbital and band-basis electron annihilation operators. We can obtain the effective attractive interaction

$$H _ { \text{int} } = - g \sum _ { q } \hat { \theta } _ { q } ^ { \dagger } \hat { \theta } _ { q }, \quad \begin{array} { c } \Gamma \omega \\ \beta \sum \\ \hat { \theta } _ { q } = \sum _ { k } \Lambda ^ { * } ( k, q ) \hat { a } _ { - k \downarrow } \hat { a } _ { k + q \uparrow }, \end{array} \quad ( 2 ) \quad \begin{array} { c } \Gamma \omega \\ \beta \sum \\ \text{filii} \end{array}$$

in which Λ( k q , ) = ∑ α u k + q ↑ ( α u ) - ↓ k ( α ) is the form factor, which encodes the information of quantum geometry [38, 39]. Here the effective coupling constant g ≈ U/N orb N c , where N orb is the number of orbitals in each unit cell. For convenience, we naturally introduce the ideal condition such that the attractive interaction strength g is large compared to the bandwidth W of the isolated flat band but small compared to the band gap ∆ gap ( W ≪ g ≪ ∆ gap ). From Eq. (2), we can arrive at an effective one-band description in the functional integral formalism [50], and the effective action can be expressed as

$$S = \sum _ { k, \sigma } \bar { a } _ { k \sigma } \left ( - i \omega _ { n } + \xi _ { k \sigma } \right ) a _ { k \sigma } - g \sum _ { q } \bar { \theta } _ { q } \theta _ { q }. \quad ( 3 ) \quad \text{In} _ { \quad \text{bein} _ { \downarrow } }$$

Here we have used the notations k = ( k , ω n ) and q = ( q , Ω m ) with ω n (Ω m ) being the fermionic (bosonic) Matsubara frequencies. Then we adopt the standard Hubbard-Stratonovich decoupling of Eq. (3) to derive the effective action S [∆ ] (for details see Supplemental Mateq rial Note I [51]) of the superconducting order parameter field ∆ . q By retaining the second order terms of ∆ q in

S [∆ ], the Gaussian action around the normal state sadq dle point ∆ q = 0 reads

$$S _ { G } \left [ \Delta _ { q } \right ] = \sum _ { q } \Gamma _ { q } ^ { - 1 } | \Delta _ { q } | ^ { 2 },$$

where the coefficient Γ -1 q = g -1 -χ c q is the pairing propagator, and the static pair susceptibility χ c q ≡ χ c q =( q , 0) reads

$$\underset { \underset { \underset { \underset { \text{i.e.} } { \iota } } { \iota } } { \iota } \quad \chi _ { q } ^ { c } = \frac { 1 } { N _ { c } } \sum _ { k } \frac { 1 - n _ { \text{F} } \left ( \xi _ { k + q } \right ) - n _ { \text{F} } \left ( \xi _ { k } \right ) } { \xi _ { k + q } + \xi _ { k } } | \Lambda ( k, q ) | ^ { 2 }. \quad ( 5 )$$

Here we use the notation as ξ k = ξ k ↑ = ξ - ↓ k and n F is the Fermi-Dirac distribution. The number equation near the critical temperature is N = -β -1 ∂S G [∆ q = 0] /∂µ , where N is the total particle number. This yields 2 A ν = ∑ k [1 -tanh ( βξ k / 2)], where A is the area of the first Brillouin zone, β -1 = k B T , and ν is the filling factor of the target band (0 ≤ ν ≤ 1).

̸

Critical temperature of FFLO instability. -By solving the linearized gap equation gχ c q = 1 together with the number equation, we can derive the critical temperature T c q for the instability of a static, spatially non-uniform pairing state characterized by a finite q index. Restricting the analysis to singleq states, ∆ q ∝ δ q Q , , the value of Q is found at the highest T c q . If Q = 0 ( Q = 0), the pairing instability is associated with the FFLO (BCS) phase.

Importantly, in the flat-band limit, we find χ c q = β ∑ k | Λ( k q , ) | 2 (1 / 2 -ν / N ) [ c ln( ν -1 -1)]. Rewriting Eq. (1) as ¯ d 2 k + q k , = 1 - | Λ( k q , ) | 2 , T c q with fixed band filling ν and coupling constant g can be derived as

$$T _ { c q } / T _ { 0 } = 1 - \bar { d } _ { \mathcal { A } } ^ { 2 } ( q ),$$

in which the average of the AQD over the Brillouin zone within the band is defined as ¯ d 2 A ( q ) = ∑ k ¯ d 2 k + q k , /N c , and T 0 is a constant which is independent of q . The highest critical temperature should be determined by T c Q = T 0 { 1 -min[ ¯ d 2 A ( q )] } over different q in the first Brillouin zone, and the corresponding optimal momentum is q = Q . We have clearly seen that the flat-band FFLO instability associated with finite-momentum Q is solely determined by the quantum geometric quantity [ ¯ d 2 A ( q )] from Eq. (6), and is independent of the interaction strength and band filling. This is the main result of our work.

̸

In a T -symmetric system, u k ↑ = u ∗ - ↓ k renders ¯ d k + q k , being the ordinary quantum distance. Since ¯ d 2 k k , = 0, the quantity ¯ d 2 A ( q ) always has a global minimum at q = 0. However, as we will discuss, the presence of QGD explicitly breaks the T symmetry and shifts the large pair susceptibility away from q = 0, allowing the global minimum of ¯ d 2 A ( q ) to occur at q = 0. Therefore, with Eq. (6) in hand, it becomes clear that the concept of QGD is both natural and useful for stabilizing the FFLO state in an isolated flat band.

Model Hamiltonian. -To study the Cooper pairings with QGD, now we examine a concrete model [32, 52-54] that allows for continuous tuning of the quantum geometry. As shown in Fig. 2(a), the ζ -lattice model features two orbitals (A and B) per site on a square lattice, described by the Hamiltonian H = ∑ k σ ˆ c † k σ ( H k σ -µ σ )ˆ c k σ , where

$$\mathcal { H } _ { k \sigma } ( \zeta _ { \sigma } ) = - t \left [ \lambda _ { x } \sin ( \zeta _ { \sigma } p _ { k } ) + m _ { \sigma } \lambda _ { y } \cos ( \zeta _ { \sigma } p _ { k } ) \right ], \ \ ( 7 ) \quad \ \ _ { - }$$

with p k = cos k x + cos k y . For brevity we set the lattice constant a = 1 and µ ↑ = µ ↓ = µ . With the spin index σ = ↑ ↓ , , the fermion operators are defined as ˆ c k σ = (ˆ c A, k σ , c ˆ B, k σ ) T . The Pauli matrices λ i act on orbital space, and time-reversal flavors are related by m σ = ↑ / ↓ = ± 1. The periodic function ζ σ p k contains spindependent long range hoppings [see Fig. 2(a)]. As plotted in Fig. 2(b), this model Hamiltonian possesses a pair of ideal isolated flat bands with band gap ∆ gap = 2 . t For each spin flavor, the band dispersion ξ k σ as well as Bloch function u k σ can be solved as ξ k σ = ± -t µ and

$$u _ { k \sigma } = \frac { 1 } { \sqrt { 2 } } \begin{pmatrix} \pm 1 \\ e ^ { i m _ { \sigma } \zeta _ { \sigma } p _ { k } } \end{pmatrix}. \quad \quad \ ( 8 ) \quad \begin{matrix} \ F I C \\ \text{pin} \\ \text{Th} \end{matrix}.$$

̸

̸

̸

Given the Bloch function u k σ , the spin-dependent quantum metrics can be expressed as G ab ( k , σ ) = ζ 2 σ sin( k a ) sin( k b ) / 4. When considering electrons with opposite momentum k that yield spin-singlet pairing, the time-reversal partners ( k , ↑ ) and ( -k , ↓ ) form a Cooper pair. Crucially, ( k , ↑ ) and ( -k , ↓ ) are degenerate in energy as ξ k ↑ = ξ - ↓ k , but are polarized in quantum metrics when ζ ↑ = ζ ↓ , meaning G ( k , ↑ ) = G -( k , ↓ ). In this case, T symmetry of the system is explicitly broken, which is evidenced by the u k ↑ = u ∗ - ↓ k , rather than by an energy difference. For the ζ -lattice model, we define the dimensionless parameter η = | ζ ↑ -ζ ↓ | / ζ ( ↑ + ζ ↓ ) as a global measure of QGD. In the following, we will illustrate the role of QGD in the formation of flat-band FFLO states.

̸

QGD-stabilized flat-band FFLO state. -In terms of Eq. (6), in Fig. 2(c) we plot T c Q ( η ) (dashed blue) as well as the corresponding Q ( η ) (solid green). As a comparison, T c 0 ( η ) (solid red) is also depicted. In the BCS phase where η &lt; η c = 0 22, . T c Q ( η ) is identical to T c 0 ( η ), and Q ( η ) = 0. However, when η &gt; η c , T c Q ( η ) = T c 0 ( η ), and a sharp transition occurs in Q ( η ) [from Q η ( c -) = 0 to Q η ( c +) ≈ 0 31 . π ], implying a transition from BCS to FFLO phase. The tri-critical point (orange dot) of BCS, FFLO, and the normal phase is the so-called Lifshitz point, which is located at ( η c = 0 22 . , T = 0 57 . T 0 ). We also find that Q is four-fold degenerate, which are located at ± Q x (ˆ + ˆ) and y ± Q ( -ˆ + ˆ) [51], x y so in our calculations we take Q = Q x (ˆ+ˆ) throughout this work. y

To further understand the QGD-driven FFLO transition and the role of AQD, we plot the T c q ( η ) evaluated from Eq. (6) in Fig. 2(d), for η = 0 (pink curve) and η = 0 25 (orange curve), respectively. . The highest critical temperature is determined by the minimum of ¯ d 2 A ( q ),

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-08b5006b-efe6-4862-be9b-8b7b56d11bd0.jpg)

̸

FIG. 2. (a) The ζ -lattice with spin-dependent long-range hoppings. ζ σ denotes the quantum metric of spinσ sector. (b) The energy dispersions in flat-band limit ( W = 0) and finite bandwidth ( W = 0). (c) T c Q ( η ) (dashed blue), T c 0 ( η ) (solid red), and the corresponding Q η ( ) (solid green) from Eq. (6) as a funtion of η . The orange dot denotes the Lifshitz point. The Q ( η ) calculated from the self-consistent mean-field (MF) theory is plotted in the dashed green line as a comparison. (d) Calculated T c q ( η /T ) 0 from Eq. (6) on q = (ˆ+ˆ) q x y , q ∈ [0 , π ], for η = 0 and η = 0 25, respectively. . In our calculations, we set ζ ↓ = 3 and ζ ↑ = ζ ↓ (1 + η / ) (1 -η ).

̸

̸

which corresponds to the maximum of T c q . As we mentioned before, T symmetry ( η = 0) ensures the favorable BCS pairing state ( Q = 0) because of the global maximum of T c q always occurs at q = 0 (dashed green circle on the pink curve). Note that a local maximum of T c q (0) can be seen at nonzero q in Fig. 2(d). As QGD ( η = 0) decreases T c 0 ( η ) significantly (the green dashed circle on the orange curve), it is possible that a nonzero Q can be stabilized. Upon increasing η , the local maximum finally transitions to global maxima (the cyan dashed circle on the orange curve) at Q η ( ) = 0, where a FFLO phase can naturally manifest. Here our calculations are based on the square lattice model, and in Supplemental Material Note III [51] we also show a triangular lattice model, where our theory still applies.

Self-consistent mean field calculations -The pair susceptibility analysis enables us to find FFLO instabilities in the single-band description. Here we numerically show, based on self-consistent mean field study [40], there is a first-order phase transition from BCS to FFLO phase driven by QGD within the previously introduced ζ -lattice model [Eq. (S23)]. Near the critical temperature, the numerical results perfectly match Eq. (6). The

FIG. 3. The BCS-FFLO phase diagram, without (a) and with (b) band dispersion. The dashed black line labels T c 0 ( η ), while the solid red line labels T c Q ( η ). The BCS and FFLO phases are colored in light blue and yellow, respectively. Parameters used here: U = 0 4 , . t ν = 0 5, . ζ ↓ = 3, W/U = 0 in (a) and W/U = 0 12 in (b). . (c) (d) (e): The free energy density landscapes of three highlighted dots in (a) at T = 0 1 . T 0 and different η ( η = 0 10, 0 23, and 0 29, respectively). . . .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-00ac22ac-65cf-4f16-abd3-e5488d2538ec.jpg)

mean-field decoupling of the attractive Hubbard interaction yields orbital-dependent order parameters given by ∆ iα = -U c ⟨ ˆ iα ↓ ˆ c iα ↑ ⟩ = ∆ ( α q ) e i q r · i with q index. In the canonical ensemble, the Helmholtz free energy density F is written as

$$F ( q ) = & - \frac { 1 } { \beta N _ { c } } \sum _ { k } \text{Tr} \left [ \ln \left ( 1 + e ^ { - \beta H _ { B d G } ( k, q ) } \right ) \right ] \quad \begin{pmatrix} \text{To $sl$} \\ \text{trans} \\ \text{Block} \\ \end{pmatrix} \\ & + \frac { 1 } { U } \sum _ { \alpha } | \Delta _ { \alpha } ( q ) | ^ { 2 }, \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \\ \text{hopp} \\ \text{Flag.\}$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $sl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \ \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end</#} \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gls} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatarray} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{</#} \begin{pmatrix} 0 \end{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{\end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{\end{} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{\end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{\end{} \text{To $gl$} \quad \begin{pmatrix} 0 \end{\end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} \text{To $gl$} \quad \begin{pmatrix} 0 \end{pmatrix} 0 \end$$

where the Bogoliubov-de Gennes (BdG) Hamiltonian H BdG ( k q , ) reads

$$H _ { \text{BDG} } ( k, q ) = \left ( \begin{array} { c c } \mathcal { H } _ { k + q \uparrow } - \mu _ { q } & \hat { \Delta } _ { q } \\ \hat { \Delta } _ { q } ^ { \dagger } & - \mathcal { H } _ { - k \downarrow } ^ { * } + \mu _ { q } \end{array} \right ), \quad ( 1 0 ) \quad \text{band} \\ \text{ deeper}$$

with the order parameters ˆ ∆ q = diag[∆ A ( q ) , ∆ ( )]. B q ˆ ∆ and q µ q can be determined self-consistently by solving the gap and number equations, respectively. Substituting ∆ ˆ q and µ q into Eq. (9), we can obtain F ( q ) for each q (see more details in Supplemental Material Note II [51]).

The results from mean-field calculations are summarized in Fig. 3. As shown in Fig. 3(a), the BCS as well as FFLO phases are colored in light blue and yellow, respectively. T c 0 ( η ) and T c Q ( η ) are plotted in the dashed black and the solid red curves. The most favored pairing momentum Q ( η ) is independent of the temperature, resulting in a straight line for the boundary between the BCS and FFLO phases, which implies that the phase transition can only be driven by QGD, rather than the thermal

̸

fluctuations. The calculated T c 0 ( η ), T c Q ( η ) and Q ( η ) [ Q ( η ) is also plotted in Fig. 2(a) to compare with the one evaluated from the band projection method] exhibits quantitatively good agreements with the corresponding analytical results in Fig. 2(c). To provide a clear picture of the phase transition, in Fig. 3(c) to (e), we show the landscapes of free energy F ( q ) for the three points (blue, green, and pink dots) in Fig. 3(a). Initially, in the BCS region ( η &lt; η c ) [Fig. 3(c)], the global minimum of the free energy F ( q ) is at q = 0, which is stable under weak fluctuations. Some local minima exist at q = 0 but are metastable. When η &gt; η c [Fig. 3(d)], q = 0 turns into a local minimum but remains a metastable state, while some q = Q points become global minima (the model respects C 4 rotational symmetry, so the minima have four-fold degeneracy in the whole Brillouin zone, see Supplemental Material Note III [51]), the FFLO phase can manifest. This indicates that the QGD energetically penalizes BCS pairing more than FFLO pairing. A firstorder phase transition from BCS to FFLO occurs near η c , which is consistent with the results obtained from the band-projection method. As η continues to increase [Fig. 3(e)], the local minimum at q = 0 eventually nearly disappears.

Robustness at finite bandwidth. -The QGD-induced BCS-FFLO phase transition has been examined in the ideal flat-band case, both analytically and numerically. However, a natural concern arises about whether the theory will break down at finite band dispersion. To show the robustness of QGD-induced BCS-FFLO transition with finite band dispersion, we modify the Bloch Hamiltonian in Eq. (S23) as H ′ k σ = H k σ -W (cos k x +cos k y ) λ / 0 4, where an additional nearest hopping term brings a finite bandwidth W as shown in Fig. 2(b). In this case, the Bloch wave functions remain unchanged. We keep W ≪ U to ensure a narrow bandwidth compared to the interaction. In Fig. 3(b), the calculated phase diagram at W/U = 0 12 . ( U = 0 4 ) is shown. . t Obviously, compared to the exactly flat band case, the transition boundary line is temperaturedependent and becomes curved due to the nonzero W , requiring a larger QGD to reach the FFLO region compared with Fig. 3(a), especially at lower temperatures. This fact indicates that increasing the temperature may also drive the BCS-FFLO phase transition (e.g., along the dashed pink line at η = 0 25 in . Fig. 3(b)). The reason is, in Eq. (5), the form factor gets Fermi-surface averaged for a dispersive band, which weakens the effects of QGD. And the thermal excitation energy k B T influences the Fermi-surface average in Eq. (5), reshaping the transition boundary line. On the other hand, for conventional superconductors with a hkighly dispersive band, we have W ≫ U ≫ k B T in the weak coupling regime, where the QGD-driven mechanism for the FFLO state may not be effective (more discussions, see Supplemental Material Note II [51]), which is consistent with our

intuition.

Conclusion and Discussion. -In this work, we have identified the QGD as the quantum-geometric origin of the FFLO pairing phase for flat-band electrons. This mechanism differs from other pair-breaking mechanisms, such as the Zeeman effects in helical superconductors [55] or the orbital effects in layered superconductors [14]. Although we only consider the spin-singlet pairing channel, there is no obstacle in generalizing it to other pseudospin degrees of freedom (such as valley and layer) and unconventional pairing channels.

Besides, throughout this work we expect the LO-type pairing state to be more stable than the FF type due to the ± Q degeneracy (see Fig. 3) for the inversionsymmetric model Hamiltonian [Eq. (S23)]. By turning on an extra inversion-breaking term, the ± Q degeneracy is lifted, resulting in the FF-type pairing state and emergent superconducting diode effect (for more details, see Supplemental Material IV[51]), driven by QGD. This is an analog to the helical superconductors with Rashba spin-orbit interaction and in-plane magnetic fields [8, 56].

Several possible routes can manifest QGD in solid-state systems. For example, optical cavities have been proposed to engineer the topology and quantum geometry of the electronic states [57, 58]. Another option is that in ultracold atom systems, the hopping parameters can be controlled artificially [8, 59]. These schemes remain the topics for future study.

Acknowledgements -We thank Patrick A. Lee, XiaoGang Wen, Xilin Feng, Ying-Ming Xie, Akito Daido, and Qing-Dong Jiang for inspiring discussions. K. T. L. acknowledges the support of the Ministry of Science and Technology, China, and Hong Kong Research Grant Council through Grants No. 2020YFA0309600, No. RFS2021-6S03, No. C6025-19G, No. AoE/P701/20, No. 16310520, No. 16307622, and No. 16309223.

- ∗ These authors contributed equally to this work.
- † jhuphy@ust.hk
- ‡ phlaw@ust.hk
- [1] M. Sigrist and K. Ueda, Reviews of Modern physics 63 , 239 (1991).
- [2] V. P. Mineev, K. Samokhin, and L. Landau, Introduction to unconventional superconductivity (CRC Press, 1999).
- [3] E. Bauer and M. Sigrist, Non-centrosymmetric superconductors: introduction and overview , Vol. 847 (Springer Science &amp; Business Media, 2012).
- [4] D. J. Scalapino, Reviews of Modern Physics 84 , 1383 (2012).
- [5] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Physical review 108 , 1175 (1957).
- [6] P. Fulde and R. A. Ferrell, Physical Review 135 , A550 (1964).
- [7] A. Larkin and Y. N. Ovchinnikov, Soviet Physics-JETP 20 , 762 (1965).
- [8] J. J. Kinnunen, J. E. Baarsma, J.-P. Martikainen, and P. T¨ orm¨, a Reports on Progress in Physics 81 , 046401 (2018).
- [9] D. F. Agterberg, J. S. Davis, S. D. Edkins, E. Fradkin, D. J. Van Harlingen, S. A. Kivelson, P. A. Lee, L. Radzihovsky, J. M. Tranquada, and Y. Wang, Annual Review of Condensed Matter Physics 11 , 231 (2020).
- [10] J. Lu, O. Zheliuk, I. Leermakers, N. F. Yuan, U. Zeitler, K. T. Law, and J. Ye, Science 350 , 1353 (2015).
- [11] X. Xi, Z. Wang, W. Zhao, J.-H. Park, K. T. Law, H. Berger, L. Forr´, o J. Shan, and K. F. Mak, Nature Physics 12 , 139 (2016).
- [12] B. T. Zhou, N. F. Yuan, H.-L. Jiang, and K. T. Law, Physical Review B 93 , 180501 (2016).
- [13] Y. Saito, Y. Nakamura, M. S. Bahramy, Y. Kohama, J. Ye, Y. Kasahara, Y. Nakagawa, M. Onga, M. Tokunaga, T. Nojima, et al. , Nature Physics 12 , 144 (2016).
- [14] C.-X. Liu, Physical Review Letters 118 , 087001 (2017).
- [15] K. W. Song and A. E. Koshelev, Physical Review X 9 , 021025 (2019).
- [16] Y.-M. Xie and K. T. Law, Physical Review Letters 131 , 016001 (2023).
- [17] N. F. Q. Yuan, Phys. Rev. Res. 5 , 043122 (2023).
- [18] K. Nakamura, A. Daido, and Y. Yanase, Physical Review B 109 , 094501 (2024).
- [19] P. Wan, O. Zheliuk, N. F. Yuan, X. Peng, L. Zhang, M. Liang, U. Zeitler, S. Wiedmann, N. E. Hussey, T. T. Palstra, et al. , Nature , 1 (2023).
- [20] S. Peotta and P. T¨rm¨, Nature communications o a 6 , 8944 (2015).
- [21] A. Julku, S. Peotta, T. I. Vanhala, D.-H. Kim, and P. T¨rm¨ a, Physical Review Letters o 117 , 045303 (2016).
- [22] L. Liang, T. I. Vanhala, S. Peotta, T. Siro, A. Harju, and P. T¨rm¨ a, Physical Review B o 95 , 024515 (2017).
- [23] X. Hu, T. Hyart, D. I. Pikulin, and E. Rossi, Physical Review Letters 123 , 237002 (2019).
- [24] A. Julku, T. J. Peltonen, L. Liang, T. T. Heikkil¨, a and P. T¨rm¨ a, Physical Review B o 101 , 060505 (2020).
- [25] F. Xie, Z. Song, B. Lian, and B. A. Bernevig, Physical review letters 124 , 167002 (2020).
- [26] V. Peri, Z.-D. Song, B. A. Bernevig, and S. D. Huber, Physical review letters 126 , 027002 (2021).
- [27] P. He, H.-T. Ding, and S.-L. Zhu, Physical Review A 103 , 043329 (2021).
- [28] J. Herzog-Arbeitman, V. Peri, F. Schindler, S. D. Huber, and B. A. Bernevig, Physical review letters 128 , 087002 (2022).
- [29] P. T¨ orm¨, a S. Peotta, and B. A. Bernevig, Nature Reviews Physics 4 , 528 (2022).
- [30] P. T¨rm¨ a, L. Liang, o and S. Peotta, Physical Review B 98 , 220511 (2018).
- [31] K.-E. Huhtinen, J. Herzog-Arbeitman, A. Chew, B. A. Bernevig, and P. T¨rm¨, Physical Review B o a 106 , 014518 (2022).
- [32] J. S. Hofmann, E. Berg, and D. Chowdhury, Physical Review Letters 130 , 226001 (2023).
- [33] G. Jiang and Y. Barlas, Physical Review Letters 131 , 016002 (2023).
- [34] W. Chen and W. Huang, Science China Physics, Mechanics &amp; Astronomy 66 , 287212 (2023).
- [35] Z. Han, J. Herzog-Arbeitman, B. A. Bernevig, and S. A. Kivelson, arXiv preprint arXiv:2401.04163 (2024).
- [36] N. S. Ticea, S. Raghu, and Y.-M. Wu, arXiv preprint arXiv:2403.00156 (2024).

- [37] H.-X. Wang and W. Huang, arXiv preprint arXiv:2406.17187 (2024).
- [38] S. A. Chen and K. T. Law, Phys. Rev. Lett. 132 , 026002 (2024).
- [39] J.-X. Hu, S. A. Chen, and K. T. Law, arXiv preprint arXiv:2308.05686 (2023).
- [40] T. Kitamura, A. Daido, and Y. Yanase, Physical Review B 106 , 184507 (2022).
- [41] P. T¨rm¨ a, Physical Review Letters o 131 , 240001 (2023).
- [42] H. Tian, X. Gao, Y. Zhang, S. Che, T. Xu, P. Cheung, K. Watanabe, T. Taniguchi, M. Randeria, F. Zhang, et al. , Nature 614 , 440 (2023).
- [43] J. Yu, C. J. Ciccarino, R. Bianco, I. Errea, P. Narang, and B. A. Bernevig, Nature Physics , 1 (2024).
- [44] A. Shapere and F. Wilczek, Geometric phases in physics , Vol. 5 (World scientific, 1989).
- [45] V. Buˇ zek and M. Hillery, Physical Review A 54 , 1844 (1996).
- [46] V. Dodonov, O. Man'Ko, V. Man'Ko, and A. W¨ unsche, Journal of Modern Optics 47 , 633 (2000).
- [47] J.-W. Rhim, K. Kim, and B.-J. Yang, Nature 584 , 59 (2020).
- [48] C. Setty, L. Fanfarillo, and P. Hirschfeld, Nature Communications 14 , 3181 (2023).
- [49] M. Iskin, Physical Review B 107 , 224505 (2023).
- [50] A. Altland and B. Simons, Condensed Matter Field Theory , 3rd ed. (Cambridge University Press, 2023).
- [51] The Supplemental material of 'Flat-band FFLO State from Quantum Geometric Discrepancy'. Including five parts, I. Derivation of the Gaussian fluctuation in the band projection formalism; II. Details of the selfconsistent mean-field theory; III. Another flat-band model on the triangular lattice; IV. Some discussions on the flat-band helical pairing state.
- [52] J. S. Hofmann, D. Chowdhury, S. A. Kivelson, and E. Berg, npj Quantum Materials 7 , 83 (2022).
- [53] D. Mao and D. Chowdhury, Proceedings of the National Academy of Sciences 120 , e2217816120 (2023).
- [54] M. Thumin and G. Bouzerar, arXiv e-prints , arXiv (2024).
- [55] A. Daido, Y. Ikeda, and Y. Yanase, Physical Review Letters 128 , 037001 (2022).
- [56] M. Smidman, M. Salamon, H. Yuan, and D. Agterberg, Reports on Progress in Physics 80 , 036501 (2017).
- [57] H. H¨bener, U. De Giovannini, C. Sch¨fer, J. Andberger, u a M. Ruggenthaler, J. Faist, and A. Rubio, Nature materials 20 , 438 (2021).
- [58] F. Schlawin, D. M. Kennes, and M. A. Sentef, Applied Physics Reviews 9 (2022).
- [59] G. Jotzu, M. Messer, F. G¨rg, D. Greif, R. Desbuquois, o and T. Esslinger, Physical review letters 115 , 073002 (2015).
- [60] S. Peotta, K.-E. Huhtinen, and P. T¨rm¨, arXiv preprint o a arXiv:2308.08248 (2023).
- [61] M. Nadeem, M. S. Fuhrer, and X. Wang, Nature Reviews Physics 5 , 558 (2023).

## Supplemental Material for 'Flat-band FFLO State from Quantum Geometric Discrepancy'

Zi-Ting Sun, 1 Ruo-Peng Yu, 1 Shuai A. Chen, 1 Jin-Xin Hu, 1 K. T. Law 1

1 Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong, China

## I. EFFECTIVE ACTION IN THE BAND-PROJECTION FORMALISM

To study the physics in a multi-band system with pairing instability, it is convenient to project the attractive interaction onto the target flat band. In this section, we will derive the effective action for the single-band description and we will show the full numerical results in the next section. To begin with, we consider a general local attractive Hubbard interaction:

$$\hat { H } _ { \text{int} } = - U \sum _ { i \alpha } \hat { c } _ { i \alpha \uparrow } ^ { \dagger } \hat { c } _ { i \alpha \downarrow } ^ { \dagger } \hat { c } _ { i \alpha \downarrow } \hat { c } _ { i \alpha \uparrow }.$$

Here U denotes the attractive interaction strength which can give rise to Cooper pairing states. i and α are the site and orbital index, respectively. The fermion operator c † α, k σ in the momentum space can be expressed by the Fourier transform:

$$\hat { c } _ { \alpha, k \sigma } ^ { \dagger } = \frac { 1 } { \sqrt { N _ { c } } } \sum _ { i } e ^ { - i k \cdot r _ { i \alpha } } \hat { c } _ { i \alpha \sigma } ^ { \dagger }$$

$$\hat { c } _ { i \alpha \sigma } ^ { \dagger } = \frac { 1 } { \sqrt { N _ { c } } } \sum _ { k } ^ {. } e ^ { i k \cdot r _ { i \alpha } } \hat { c } _ { \alpha, k \sigma } ^ { \dagger }.$$

In this way, we can rewrite the interacting Hamiltonian in the momentum space as

$$\hat { H } _ { \text{int} } = - \frac { U } { N _ { c } } \sum _ { \alpha k k ^ { \prime } q } \hat { c } ^ { \dagger } _ { \alpha, k + q \uparrow } \hat { c } ^ { \dagger } _ { \alpha, - k \downarrow } \hat { c } _ { \alpha, - k ^ { \prime } \downarrow } \hat { c } _ { \alpha, k ^ { \prime } + q \uparrow }.$$

For the convenience of the analytical treatment, we approximate the interaction term [30, 31, 60] as

$$\hat { H } _ { \text{int} } \approx - \frac { U } { N _ { o r b } N _ { c } } \sum _ { \alpha \beta k k ^ { \prime } q } \hat { c } ^ { \dagger } _ { \alpha, k + q \uparrow } \hat { c } ^ { \dagger } _ { \alpha, - k \downarrow } \hat { c } _ { \beta, - k ^ { \prime } \downarrow } \hat { c } _ { \beta, k ^ { \prime } + q \uparrow }.$$

We use g = U/ N ( orb N c ) to characterize the effective interaction. We project the interaction onto the target flat band by ˆ c α, k σ → u ∗ k σ ( α a )ˆ k σ , where ˆ a k σ denotes the band-basis electron annihilation operator. Then from the effective ˆ H int , on the target flat band, the effective action after the band projection can be expressed as

$$S = \sum _ { k, \sigma } \bar { a } _ { k \sigma } \left ( - i \omega _ { n } + \xi _ { k \sigma } \right ) a _ { k \sigma } - g \sum _ { q } \bar { \theta } _ { q } \theta _ { q }.$$

where

$$\theta _ { q } = \sum _ { k } \Lambda ^ { * } ( k, q ) a _ { - k \downarrow } a _ { k + q \uparrow },$$

in which Λ( k q , ) = ∑ α u k + q ↑ ( α u ) - ↓ k ( α ) = ⟨ u k + q ↑ | u ∗ - ↓ k ⟩ is introduced as the form factor. Here we have used the notations k = ( k , ω n ) and q = ( q , Ω m ) with ω n (Ω m ) being the fermionic (bosonic) Matsubara frequencies.

After a standard Hubbard-Stratonovich decoupling, one can find the action in the Nambu basis ψ k,q = ( a k + q ↑ , a ¯ - ↓ k ) T takes the form as

$$S = \sum _ { q } \frac { | \Delta _ { q } | ^ { 2 } } { g } + \sum _ { k, q } \bar { \psi } _ { k, q } \hat { \mathbb { G } } _ { k, q } ^ { - 1 } \psi _ { k, q },$$

where

$$\hat { \mathbb { G } } _ { k, q } ^ { - 1 } = \left ( \begin{array} { c c } - i \omega _ { n } - i \Omega _ { m } + \xi _ { k + q \uparrow } & - \Lambda ( k, q ) \Delta _ { q } \\ - \Lambda ^ { * } ( k, q ) \Delta _ { q } ^ { * } & - i \omega _ { n } - \xi _ { - k \downarrow } \end{array} \right ).$$

After integrating out the fermionic field, we can get

$$S [ \Delta _ { q } ] = \sum _ { q } \frac { | \Delta _ { q } | ^ { 2 } } { g } - \text{Tr} \ln \hat { G } _ { k, q } ^ { - 1 },$$

and expanding it up to the second order in ∆ , the Gaussian action around the trivial saddle point ∆ q q = 0 reads

$$S _ { G } \left [ \Delta _ { q } \right ] = \sum _ { q } \Gamma _ { q } ^ { - 1 } | \Delta _ { q } | ^ { 2 },$$

in which the coefficient

$$\Gamma _ { q } ^ { - 1 } = \frac { \delta ^ { 2 } S [ \Delta _ { q } ] } { \delta \Delta _ { q } ^ { * } \delta \Delta _ { q } } \Big | _ { \Delta _ { q } = \Delta _ { q } ^ { * } = 0 } = g ^ { - 1 } - \chi _ { q } ^ { c }$$

is the pairing propagator, and the pair susceptibility χ c q ≡ χ c q =( q , 0) at the static limit (Ω m = 0) reads

$$\chi _ { q } ^ { c } = \frac { 1 } { N _ { c } } \sum _ { k } \frac { 1 - n _ { F } \left ( \xi _ { k + q \uparrow } \right ) - n _ { F } \left ( \xi _ { - k \downarrow } \right ) } { \xi _ { k + q \uparrow } + \xi _ { - k \downarrow } } | \Lambda ( k, q ) | ^ { 2 }.$$

Near the critical temperature, we also have the number equations (at the weak coupling limit), which read

$$N _ { \sigma } = - \beta ^ { - 1 } \frac { \partial S [ \Delta _ { q } = 0 ] } { \partial \mu _ { \sigma } } = \frac { 1 } { 2 } \sum _ { k } \left [ 1 - \tanh \left ( \frac { \beta \xi _ { k \sigma } } { 2 } \right ) \right ].$$

In the ideal flat-band limit, we have ξ k σ = -µ . One can find that χ c q = β ∑ k | Λ( k q , ) | 2 (1 / 2 -ν / N ) [ c ln( ν -1 -1)], where ν is the filling factor. At half filling ν = 1 2, / χ c q → β ∑ k | Λ( k q , ) | 2 / (4 N c ).

FIG. S1. (a)-(c): The ∆ -T relation for the η = 0 1 (in BCS region) states paired with the finite momenta of . q = (ˆ+ˆ), and q x y q = 0, 0 19 . π , 0 24 . π respectively. The favoured finite-momentum Q = 0. (d)-(f): The ∆ -T relation for the η = 0 23 (in FFLO . region) states paired with the same q in (a)-(c), and the favoured finite-momentum Q = 0 24 . π . In (a)-(f), ∆ A is donated by the red solid curve and ∆ B is donated by the black dashed curve.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-21f54918-41dc-4033-9052-bee39ccb976f.jpg)

## II. SELF-CONSISTENT MEAN-FIELD THEORY

In this section, we provide the details of the mean-field numerical calculations in the main text. To be specific, we do the decoupling of the on-site attractive Hubbard interaction within the mean-field approximation

$$\hat { H } _ { \text{int} } & = U \sum _ { i \alpha } \hat { c } ^ { \dagger } _ { i \alpha \uparrow } \hat { c } ^ { \dagger } _ { i \alpha \downarrow } \hat { c } _ { i \alpha \uparrow } \\ & \approx \sum _ { \alpha k } \left ( \Delta _ { \alpha } \hat { c } ^ { \dagger } _ { \alpha, k + q \uparrow } \hat { c } ^ { \dagger } _ { \alpha, - k \downarrow } + h. c. \right ) + \frac { N _ { c } } { U } \sum _ { \alpha } | \Delta _ { \alpha } | ^ { 2 },$$

where we have used ∆ α = -U c ⟨ ˆ α, - ↓ k ˆ c α, k + q ↑ ⟩ and ∆ α = ∆ ( α q ) can be determined by the self-consistent gap equation for a given q . The full mean-field Hamiltonian can be rewritten in the matrix form in Nambu basis

$$\hat { \mathcal { H } } _ { \text{MF} } = \sum _ { k, q } \psi ^ { \dagger } ( k, q ) H _ { \text{BdG} } ( k, q ) \psi ( k, q ) + \sum _ { k, q } \text{Tr} [ H _ { - k \downarrow } - \mu ( q ) \mathbb { I } _ { 2 \times 2 } ] + \frac { N _ { c } } { U } \sum _ { \alpha } | \Delta _ { \alpha } | ^ { 2 },$$

in which

$$H _ { \text{BdG} } ( k, q ) = \left ( \begin{array} { c c } H _ { k + q \uparrow } - \mu ( q ) \mathbb { I } _ { 2 \times 2 } & \text{diag} ( \Delta _ { \alpha } ) \\ \text{diag} ( \Delta _ { \alpha } ^ { \dagger } ) & - H _ { - k \downarrow } ^ { T } + \mu ( q ) \mathbb { I } _ { 2 \times 2 } \end{array} \right ),$$

where diag(∆ α ) = diag(∆ A , ∆ ), the momentum B q = q x ˆ + x q y ˆ, and the Nambu spinor reads y

$$\psi ( k, q ) = \left ( \hat { c } _ { A, k + q \uparrow }, \hat { c } _ { B, k + q \uparrow }, \hat { c } _ { A, - k \downarrow } ^ { \dagger }, \hat { c } _ { B, - k \downarrow } ^ { \dagger } \right ) ^ { T }.$$

The grand potential is defined as

$$\nolimits _ { \mu } & \nolimits [ \mu ( q ), \Delta _ { \alpha } ( q ), q ] & = - \frac { 1 } { \beta } \log Z = - \frac { 1 } { \beta } \ln \text{Tr} [ e ^ { - \beta \hat { H } _ { M P } } ] \\ & = \sum _ { k } \text{Tr} ( H _ { - k _ { l } } - \mu ( q ) \mathbb { I } _ { 2 \times 2 } ) + \frac { N _ { c } } { U } \sum _ { \alpha } | \Delta _ { \alpha } | ^ { 2 } - \frac { 1 } { \beta } \sum _ { k } \text{Tr} \left [ \ln \left ( 1 + e ^ { - \beta H _ { B a c } ( k, q ) } \right ) \right ] \\ & = \sum _ { k } \text{Tr} ( H _ { - k _ { l } } - \mu ( q ) \mathbb { I } _ { 2 \times 2 } ) + \frac { N _ { c } } { U } \sum _ { \alpha } | \Delta _ { \alpha } | ^ { 2 } - \frac { 1 } { \beta } \sum _ { n k } \ln \left ( 1 + e ^ { - \beta E _ { n } ( k, q ) } \right ),$$

1

FIG. S2. The superconducting phase diagram of the ζ -lattice model. The red solid line labels T c 0 ( η ), the black dashed line labels T c Q ( η ), where Q = Q x (ˆ + ˆ). y The BCS regime and the FFLO regime are labeled in light blue and yellow, respectively. Parameters used here: U = 0 4 , . t ζ ↓ = 3, W/U = 0 1 in (a) and . W/U = 0 2 in (b). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-b04a7ee0-2b6b-4a5d-bf02-5de9e45124dd.jpg)

where E n ( k q , ) is the n-th eigenvalue of the BdG Hamiltonian H BdG ( k q , ). By minimizing the grand potential with respect to the pairing potential ∆ α we can get the self-consistent equation:

$$& \frac { \partial \Omega } { \partial \Delta _ { \alpha } ^ { * } } = \frac { N _ { c } } { U } \Delta _ { \alpha } ( \mathbf q ) + \sum _ { k } \langle \hat { c } _ { \alpha, - k \downarrow } \hat { c } _ { \alpha, k + q \uparrow } \rangle = 0, \\ & \Delta _ { \alpha } ( \mathbf q ) = - \frac { U } { N _ { c } } \sum _ { k } \langle \hat { c } _ { \alpha, - k \downarrow } \hat { c } _ { \alpha, k + q \uparrow } \rangle = - \frac { U } { N _ { c } } \sum _ { k } \text{Tr} \left [ \mathcal { U } _ { k, q } ^ { \dagger } \mathcal { M } _ { \alpha } \mathcal { U } _ { k, q } \hat { \rho } \right ].$$

And the filling factor ν can be obtained from the thermodynamic relation:

$$N & = 2 N _ { c } \nu = - \frac { \partial \Omega } { \partial \mu } = \sum _ { k } \langle \hat { c } _ { a, k \uparrow } ^ { \dagger } \hat { c } _ { a, k \uparrow } + \hat { c } _ { a, k \downarrow } ^ { \dagger } \hat { c } _ { a, k \downarrow } + \hat { c } _ { b, k \uparrow } ^ { \dagger } \hat { c } _ { b, k \uparrow } + \hat { c } _ { b, k \downarrow } ^ { \dagger } \hat { c } _ { b, k \downarrow } \rangle, \\ 2 \nu & = 1 + \sum _ { k } \text{Tr} \left [ \mathcal { U } _ { k, q } ^ { \dagger } \tau _ { z } \otimes \mathbb { I } _ { 2 \times 2 } \mathcal { U } _ { k, q } \hat { \rho } \right ],$$

where U k q , is the unitary operator that diagonalizes the BdG Hamiltonian as U † k q , H BdG ( k q , ) U k q , = diag( E n ), M A =

$$\left ( \begin{array} { c } 0 & 0 & 1 & 0 \\ 0 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{array} \right ), \, \mathcal { M } _ { B } = \left ( \begin{array} { c } 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{array} \right ), \, \text{and the density matrix } \hat { \rho } = \frac { 1 } { 2 } \, \text{diag} \left ( \frac { 1 } { 1 + \epsilon ^ { B E } ( x, \epsilon ) } \right ). \, \text{ In turn, if the particle number}$$

$$\begin{smallmatrix} \text{where} \, \varkappa _ { k, q } \, is \, \text{me with any or per-avon } \, \text{max} \,, \\ \left ( \begin{smallmatrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{smallmatrix} \right ), \, \mathcal { M } _ { B } = \left ( \begin{smallmatrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{smallmatrix} \right ), \, \text{all}$$

(filling factor ν ) is fixed, the chemical potential µ ( q ) can be determined from this equation. Finally, the free energy density (up to a constant given by the summation 1 Nc ∑ k Tr( H - ↓ k )) is:

$$F ( q ) = \frac { 1 } { N _ { c } } \left ( \Omega [ \mu ( q ), \Delta _ { \alpha } ( q ), q ] + \mu ( q ) N \right ) = \frac { 1 } { U } \sum _ { \alpha } | \Delta _ { \alpha } | ^ { 2 } - \frac { 1 } { \beta N _ { c } } \sum _ { n k } \ln \left ( 1 + e ^ { - \beta E _ { n } ( k, q ) } \right ).$$

In our numerical calculations, we consider the trivial irreducible representation as the intra-orbital pairing channel, i.e., ∆ αβ ∝ δ αβ is picked up as the only relevant pairing channel, with the other pairing channels neglected. We note that in the previous section, the band-projected interaction Hamiltonian is approximated under the uniform pairing approximation. To show this approximation is accurate near the critical temperature in this model, we plot the ∆ -T relation at different QGD η and finite-momentum q in Fig. S1. For different η and q , at low temperatures, ∆ A can differ from ∆ B . But near the critical temperature, we find ∆ A = ∆ B , which satisfies the uniform pairing condition. This fact explains why our critical temperature analysis based on the uniform pairing condition matches the numerical calculations well.

We also provide more numerical results for the case of finite band dispersion. We keep W ≪ U ≪ 2 t to ensure a narrow bandwidth compared to the band gap and the interaction. In Fig. S2 We show the phase diagram at W/U = 0 1 and . W/U = 0 2. . Compared to Fig.3(b) in the main text, it is clear that the area of the FFLO phase is shrinking with the increase of W/U . Especially in Fig. S2(b), the FFLO states can only be found in a very narrow region. It is straightforward to infer that if we keep increasing W/U , the FFLO states may vanish, which supports our claim in the main text.

FIG. S3. Averaged anomalous quantum distance ¯ ( d 2 η q ) for the square and triangular ζ -lattice model, in the absence (presence) of QGD. (a) Square lattice for η = 0. (b) Square lattice for η = 0 6. . (c) Triangular lattice for η = 0. (d) Triangular lattice for η = 0 6. . Other parameters used here: ζ ↑ = 1, t = 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-75cb585b-a46b-49fc-b673-1b33b130079b.jpg)

## III. A FLAT-BAND MODEL ON TRIANGULAR LATTICE

In the main text, we study a ζ -lattice model on a square lattice, which respects C 4 symmetry. To further show the generality of our approach, in this section, we show similar results for the triangular lattice that has orbitals ( A and B ) per site. Similarly, we can write down the Bloch Hamiltonian as

$$& \mathcal { H } _ { k \sigma } = - t \left ( \lambda _ { x } \sin \zeta _ { \sigma } p _ { k } + m _ { \sigma } \lambda _ { y } \cos \zeta _ { \sigma } p _ { k } \right ) - \mu _ { \sigma } \lambda _ { 0 }, \\ & p _ { k } = \sum _ { j = 1 } ^ { 3 } \cos ( k \cdot R _ { j } ).$$

where R 1 = a (1 , 0) , R 1 = a ( -1 / 2 , √ 3 / 2) , R 3 = a ( -1 / 2 , -√ 3 / 2). To demonstrate the general validity of our theory in the main text, we calculate the averaged anomalous quantum distance ¯ ( d 2 η q ) for this model in the absence (presence) of QGD in Fig. S3(c) and (d). As a comparison, we also make the plots for the square lattice model presented in the main text in Fig. S3(a) and (b). In all of these figures, the dark blue region denotes the global minimum of ¯ ( d 2 η q ). We can see that in both (a) and (c), when η = 0, the global minimum of ¯ ( d 2 η q ) is at q = 0. But for (b) and (d), the global minima are not at q = 0, which implies the FFLO instability, as we explained in the main text. It is worth noting that in Fig. S3(b) the FFLO states are four-fold degenerate, while in Fig. S3(d) the FFLO states are six-fold degenerate, due to the different crystalline symmetries of the Bloch function. Our theory shows universality for various model Hamiltonians.

FIG. S4. (a) The superconducting phase diagram of ζ -lattice model with an inversion breaking term γ . The Q &gt; 0 (weak helical superconducting state) and Q &lt; 0 (strong helical superconducting state) regimes are colored in blue and yellow, respectively. Parameters used here: U = 0 4 , . t ζ ↑ = 3, γ = 3. (b) The red curve shows the finite-momentum Q = Q (1 1) versus QGD ( , η ). Meanwhile, the blue curve is the corresponding diode efficiency Θ. When increasing η to across the phase boundary η = 0 26, . the diode efficiency Θ is enhanced dramatically, and the sign of Q is reversed. (c) and (d) The free energy (red curve) and supercurrent (blue curve) of 2 highlighted points in (a) with the temperature T = 0 1 . T 0 but different η [ η = 0 24 in (c) and . 0 29 in (d)]. . Both of the free energy and supercurrent are normalized to [-1, 1]. The green point is in the Q &gt; 0 phase, and the red point is in the Q &lt; 0 phase.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00548v2-figure-123fd1a1-ff84-4bee-8f21-7fd46eed4122.jpg)

## IV. FLAT-BAND HELICAL PAIRING STATE AND SUPERCONDUCTING DIODE EFFECT

The helical superconductivity can be realized in noncentrosymmetric superconductors with the aid of in-plane magnetic field [55]. To achieve an analog of the helical pairing state in a flat-band system with QGD, we manually turn on an extra inversion breaking term in the ζ -lattice model by modifying (recall the periodic function is α k ,σ = ζ σ p k )

$$\alpha _ { k, \sigma } = \zeta _ { \sigma } p _ { k } + m _ { \sigma } \gamma ( \sin k _ { x } + \sin k _ { y } ).$$

̸

Here the term with parameter γ explicitly breaks the inversion symmetry in the Bloch wave function such that α -k ,σ = α k ,σ , while still preserving the time-reversal symmetry as α k , ↑ = α -k , ↓ when η = 0. However, when η = 0, this term entirely lifts ± Q degeneracy in the free energy [see Fig. S4(d)], resulting in FF-type pairing state. Here the Cooper pair momentum is along the Γ -M line as q = (ˆ ˆ). q x, y

̸

The mean-field results are summarized in Fig. S4. In this case, the finite momentum pairing state is favored for any nonzero η . Interestingly, Q ( η ) shows a sharp change near a critical η c , manifested as the sign change [Fig. S4 (b)], which suggests a first-order phase transition. In Fig. S4 (a), we color the Q &gt; 0 and Q &lt; 0 regimes in blue and yellow, respectively. This is an analog to the weak and strong helical superconducting states [55].

Recently, the superconducting diode effect (SDE), in which the magnitudes of the critical supercurrents differ in opposite directions, has been studied in helical superconductors [55, 61]. Inspired by this, here we also calculate the diode efficiency Θ = ( J c + -| J c -| ) / J ( c + + | J c -| ) in Fig. S4 (b), where J c ± is the upper and lower critical supercurrent in the (1 , 1)-direction. It shows that a dramatic enhancement of Θ appears near η c , which is similar to what has been proposed in Ref. [55].

We provide some details in evaluating the supercurrent. The supercurrent can be evaluated by J q ( ) = 2[ ∂F q ( )] /∂q , and the minimized free energy can be calculated from Eq. (S22). To give a clearer picture of the phase transition, in Fig. S4(c) and (d), we show the free energy F q ( ) (red curves) as well as the supercurrent J q ( ) (blue curves) for the two points (green and red dots) in Fig. S4(a). In the blue region [Fig. S4 (a)], the global minimum of the free energy F q ( ) is at Q &gt; 0 [Fig. S4 (c)], but | Q | is small, associated with the small diode efficiency. When η &gt; η c [Fig. S4 (d)], the pairing state with Q &lt; 0 becomes the new global minimum. In this case, an enhanced diode efficiency can be observed in Fig. S4(b).
Zi-Ting Sun Ruo-Peng Yu Shuai A. Chen Jin-Xin Hu K. T. Law
Flat-band FFLO State from Quantum Geometric Discrepancy
We propose a new scheme for the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)
Cooper pairing states within ultraflat bands. Central to our approach is the
concept of "quantum geometric discrepancy" (QGD), which characterizes the
discrepancy in the quantum geometry of paired electrons, giving rise to the
flat-band FFLO instability. Remarkably, we find that this instability is
directly related to a new quantum geometric quantity we term "anomalous quantum
distance", which formally captures QGD. To model both QGD and the anomalous
quantum distance, we examine a flat-band electronic Hamiltonian with tunable
spin-dependent quantum metrics. Utilizing the band-projection method, we
analyze the QGD-induced FFLO instability from pair susceptibility near the
superconducting critical temperature. Furthermore, we perform self-consistent
mean-field calculations to obtain the phase diagram of the BCS-FFLO transition
driven by QGD, which aligns well with our analytical results. We emphasize that
QGD serves as a distinctive protocol for stabilizing the flat-band FFLO phase.