## Geometric contribution to the Goldstone mode in spin-orbit coupled Fermi superfluids

## M. Iskin

Department of Physics, Ko¸ c University, Rumelifeneri Yolu, 34450 Sarıyer, Istanbul, Turkey (Dated: November 11, 2019)

The so-called quantum metric tensor is a band-structure invariant whose measure corresponds to the quantum distance between nearby states in the Hilbert space, characterizing the geometry of the underlying quantum states. In the context of spin-orbit coupled Fermi gases, we recently proposed that the quantum metric has a partial control over all those superfluid properties that depend explicitly on the mass of the superfluid carriers, i.e., the effective-mass tensor of the corresponding (two- or many-body) bound state. Here we scrutinize this finding by analyzing the collective phase and amplitude excitations at zero temperature. In particular to the Goldstone mode, we present extensive numerical calculations for the Weyl and Rashba spin-orbit couplings, revealing that, despite being small, the geometric contribution is solely responsible for the nonmonotonic evolution of the sound velocity in the BCS-BEC crossover.

## I. INTRODUCTION

The quantum geometric tensor [1-3] lies at the heart of modern solid-state and condensed-matter physics: while it describes the local geometry of the underlying Bloch states, it can be related to the global properties of some systems in most peculiar ways. For instance, its imaginary part corresponds to the Berry curvature measuring the emergent gauge field in momentum space, and it can be associated with the topological invariants that are used to classify quantum Hall states or topological insulators [4, 5]. In contrast, its real part corresponds to the quantum metric tensor measuring the quantum distance/fidelity between nearby Bloch states, and it can be associated with, e.g., the density and the effectivemass tensor of the superfluid (SF) carriers in spin-orbit coupled Fermi gases [6, 7], as well as with many other systems [8-10].

ity of the Goldstone mode has two contributions as well: in addition to the conventional intraband one that is controlled by the helicity spectrum, the geometric interband contribution is controlled by the quantum metric. We present extensive numerical calculations for the Weyl and Rashba SOCs revealing that, despite being much smaller than the conventional one, the geometric contribution is solely responsible for the nonmonotonic evolution of the sound velocity in the BCS-BEC crossover.

In the context of multi-band superconductivity, the new surge of theoretical interest in quantum-metric effects dates back only to 2015, soon after its direct association with the geometric origin of superfluidity in topologically nontrivial flat bands [11]. It was found that the SF weight (also known as the SF stiffness tensor) has two distinct contributions based on the physical mechanisms involved. In addition to the conventional contribution that is controlled by the electronic spectra of the Bloch bands, there is also the so-called geometric contribution coming from the interband processes [11, 12]. It later turns out that the quantum geometry has a partial control not only over the SF weight but also over all those SF properties that depend explicitly on the effective-mass tensor of the corresponding (two- or many-body) bound state [7]. There have been many studies on the subject exploring a variety of multi-band Hamiltonians, including most recently that of the twisted bilayer graphene [1315].

Motivated by the ongoing experimental progress [1620] and given our recent findings [6, 7, 21], here we consider a two-component Fermi gas with an arbitrary spinorbit coupling (SOC), and analyze its collective excitations in the BCS-BEC crossover. We show that the veloc-

The rest of this paper is organized as follows. In Sec. II, after introducing the many-body Hamiltonian, we first describe the effective-Gaussian action for the phase- and amplitude-only degrees of freedoms, and then derive analytical expressions for the collective modes with a special emphasis on the distinction between the conventional and geometric contributions. Much of the theoretical details on the effective-action approach are deferred to the Appendix A. The resultant Goldstone (phase) mode is analyzed Sec. III, where we present numerical results for a 3D Fermi gas with a Weyl or a Rashba SOC, and for a 2D Fermi gas with a Rashba SOC for completeness. We end the paper with a brief summary of our conclusions in Sec. IV.

## II. COLLECTIVE EXCITATIONS

Having a two-component ( ↑ and ↓ ) Fermi gas with an arbitrary SOC that is driven across an s -wave Feshbach resonance in mind, we consider the Hamiltonian

$$\underset { \text{less} } { \text{use} } \quad H & = \sum _ { k } ( a _ { \uparrow k } ^ { \dagger } a _ { \downarrow k } ^ { \dagger } ) \left [ \xi _ { k } \sigma _ { 0 } + \mathbf d _ { k } \cdot \sigma \right ] \left ( \begin{array} { c } a _ { \uparrow k } \\ a _ { \downarrow k } \end{array} \right ) \\ \text{ject} \quad \ \ & - U \sum _ { k k ^ { \prime } \mathbf q } a _ { \uparrow, k + \mathbf q / 2 } ^ { \dagger } a _ { \downarrow, - k + \mathbf q / 2 } ^ { \dagger } a _ { \downarrow, - k ^ { \prime } + \mathbf q / 2 } a _ { \uparrow, k ^ { \prime } + \mathbf q / 2 }, \quad ( 1 )$$

where a † σ k ( a σ k ) creates (annihilates) a pseudospinσ fermion with momentum k , i.e., in units of glyph[planckover2pi1] → 1 the Planck constant. In the first line, ξ k = glyph[epsilon1] k -µ is the usual free-particle dispersion glyph[epsilon1] k = k / 2 (2 m ) shifted by the chemical potential µ , d k = ∑ i d i k ̂ i is the SOC field

with ̂ i denoting the unit vector along the i = ( x, y, z ) direction, σ 0 is a 2 × 2 identity matrix, and σ = ∑ i σ i ̂ i is a vector of Pauli spin matrices. Our notation is such that d i k = α k i i corresponds to the Weyl SOC when α i = α for all i , and to the Rashba SOC when α z = 0. We choose α ≥ 0 without the loss of generality. Thus, the one-body problem is described by the Hamiltonian density h 0 k = glyph[epsilon1] k σ 0 + d k · σ , determining the helicity spectrum glyph[epsilon1] s k = glyph[epsilon1] k + sd k and the helicity eigenstates | s k 〉 via the wave equation h 0 k | s k 〉 = glyph[epsilon1] s k | s k 〉 . Here, s = ± labels the helicity bands, d k = | d k | , and the unit vector ̂ d k = d k /d k plays an important role throughout the paper.

In the second line of Eq. (1), U ≥ 0 corresponds to the strength of the zero-ranged density-density interactions between ↑ and ↓ particles, and we tackle this term by making use of the imaginary-time functional pathintegral formalism for the paired fermions. This formalism offers a systematic approach to go beyond the meanfield theory as the fluctuation effects can be incorporated into the effective action order by order as a power-series expansion around the stationary saddle-point [22-25]. As shown in the Appendix, the zeroth-order (i.e., the saddlepoint) action coincides with the mean-field one, and can be written as

$$S _ { 0 } = \frac { \Delta _ { 0 } ^ { 2 } } { T U } + \sum _ { s \mathbf k } \left \{ \frac { \xi _ { s \mathbf k } - E _ { s \mathbf k } } { 2 T } + \ln [ f ( - E _ { s \mathbf k } ) ] \right \}, \quad ( 2 ) \quad M _ { a. C } ^ { 1 1 }$$

where ∆ 0 = U a 〈 ↑ k a ↓ -, k 〉 is the saddle-point order parameter with 〈 . . . 〉 denoting a thermal average, T is the temperature in units of k B → 1 the Boltzmann constant, E s k = √ ξ 2 s k +∆ is the quasiparticle energy 2 0 spectrum associated with ξ s k = glyph[epsilon1] s k -µ, and f ( x ) = 1 / e ( x/T + 1) is the Fermi-Dirac distribution. Here, we take ∆ 0 to be a real parameter without the loss of generality, and determine its value by imposing the saddlepoint condition ∂S /∂ 0 ∆ 0 = 0 on the order parameter, leading to 1 = U ∑ s k [1 -2 f ( E s k )] / (4 E s k ) . Furthermore, µ is determined by the thermodynamic relation N 0 = -T∂S /∂µ 0 for the number of particles, leading to N 0 = ∑ s k { 1 -( ξ s k /E s k )[1 -2 f ( E s k )] } / 2 . These saddlepoint equations constitute the self-consistent BCS-BEC crossover theory for pairing, and they coincide with those of the mean-field approximation. The coupled solutions for ∆ 0 and µ have served as a reliable guide for the coldatom community, providing a qualitative description of the SF state at sufficiently low temperatures across an s -wave Feshbach resonance [22-25].

Going beyond S 0 in the Appendix, we first introduce the bosonic field ∆ q = ∆ 0 + Λ q , where Λ q corresponds to the fluctuations of the order parameter around its stationary value ∆ 0 with the collective index q = ( q , i ν n ) denoting both the pair momentum q and the bosonic Matsubara frequency ν n = 2 πnT , and then derive the effective-Gaussian action S Gauss = S 0 + S 2 . Note that the first order term S 1 vanishes thanks to the saddlepoint condition on the order parameter. Furthermore, it is physically insightful to separate the phase and amplitude degrees of freedoms through the unitary transformation Λ q = ( λ q + i θ q ) / √ 2 , leading to a phase-only action S θ 2 and to a amplitude-only one S λ 2 that are of the form

$$\underset { \underset { \underset { \underset { \underset { \underset { \underset { \underset { \text{ecl} } { \epsilon } } } } } { \text{the} } } } { \text{the} } \\ \underset { \underset { \text{ecl} } } { \text{the} } } \underset { \underset { \text{ecl} } } { S } ^ { \theta } _ { 2 } = \frac { 1 } { 2 T } \sum _ { q } \theta _ { q } ^ { * } \left [ M _ { q, E } ^ { 1 1 } - M _ { q } ^ { 1 2 } - \frac { ( M _ { q, O } ^ { 1 1 } ) ^ { 2 } } { M _ { q, E } ^ { 1 1 } + M _ { q } ^ { 1 2 } } \right ] \theta _ { q }, \ \ ( 3 )$$

$$\stackrel { | s \mathbf k \rangle } { = } \pm \sum _ { \substack { S _ { 2 } ^ { \lambda } = \frac { 1 } { 2 T } \sum _ { q } \lambda _ { q } ^ { * } \left [ M _ { q, E } ^ { 1 1 } + M _ { q } ^ { 1 2 } - \frac { ( M _ { q, O } ^ { 1 1 } ) ^ { 2 } } { M _ { q, E } ^ { 1 1 } - M _ { q } ^ { 1 2 } } \right ] } \lambda _ { q }. \ ( 4 )$$

Here, since M 11 q,O couples the phase and amplitude modes, the phase-only action is derived by integrating out the amplitude fields, and vice versa [23]. While the thermal expressions for the matrix elements M ij q of the inversefluctuation propagator M q are given in the Appendix, here we restrict ourselves to the collective excitations on top of the thermal ground state. Thus, setting T = 0, and denoting ξ s, k + q / 2 by ξ ; ξ s , ′ -k + q / 2 by ξ ′ ; E s, k + q / 2 by E E ; s , ′ -k + q / 2 by E ′ ; and ̂ d ± k + q / 2 by ̂ d ± , a compact way to write the matrix elements is [23]

$$\underset { ( 2 ) } { \underset { M _ { q }, E } { \text{$\Lambda$} } } = \frac { 1 } { U } + \sum _ { s s ^ { \prime } \mathfrak { k } } \frac { ( \xi \xi ^ { \prime } + E E ^ { \prime } ) ( E + E ^ { \prime } ) ( 1 - s s ^ { \prime } \widehat { \mathbf d } _ { + } \cdot \widehat { \mathbf d } _ { - } ) } { 8 E E ^ { \prime } [ ( \mathfrak { i } \nu _ { n } ) ^ { 2 } - ( E + E ^ { \prime } ) ^ { 2 } ] }, \\ \wedge \quad \wedge$$

$$\Gamma _ { \Gamma } ^ { ( 2 ) } \quad M _ { q, O } ^ { 1 1 } = \sum _ { s s ^ { \prime } \mathbf k } \frac { ( \xi E ^ { \prime } + E \xi ^ { \prime } ) \mathbf i \nu _ { n } ( 1 - s s ^ { \prime } \widehat { \mathbf d } _ { + } \cdot \widehat { \mathbf d } _ { - } ) } { 8 E E ^ { \prime } [ ( \mathbf i \nu _ { n } ) ^ { 2 } - ( E + E ^ { \prime } ) ^ { 2 } ] }, \quad \ \quad ( 6 )$$

$$\frac { \ p a _ { i } } { T \text{ is} } \sum _ { \text{s} n } M _ { q } ^ { 1 2 } = - \sum _ { s s ^ { \prime } \mathbf k } \frac { \Delta _ { 0 } ^ { 2 } ( E + E ^ { \prime } ) ( 1 - s s ^ { \prime } \widehat { \mathbf d } _ { + } \cdot \widehat { \mathbf d } _ { - } ) } { 8 E E ^ { \prime } [ ( \mathbf i \nu _ { n } ) ^ { 2 } - ( E + E ^ { \prime } ) ^ { 2 } ] }. \quad \quad ( 7 )$$

Here, M 11 q,E and M 12 q are even functions of i ν n , but M 11 q,O is an odd one.

glyph[negationslash]

glyph[negationslash]

The dispersions ω q for the collective modes are determined by the poles of the propagator matrix M -1 q for the pair fluctuation fields, i.e., ω θ q for the phase mode and ω λ q for the amplitude one are determined, respectively, by setting [ . . . ] = 0 in Eqs. (3) and (4) after an analytic continuation i ν n → ω + 0 i + to the real axis [22-25]. For this purpose, it is sufficient to retain terms that are up to quadratic orders in their small q and ω expansions [26], i.e., M 11 q,E + M 12 q = A + ∑ ij C q q ij i j -Dω 2 + · · · ; M 11 q,E -M 12 q = ∑ ij Q q q ij i j -Rω 2 + · · · ; and M 11 q,O = -Bω + · · · . This calculation leads to a gapless phase (Goldstone) mode ( ω θ q ) 2 = ∑ ij x ij q q i j with x ij = Q / R ij ( + B /A 2 ) for the phase-only action, and to a gapped amplitude (Higgs) mode ( ω λ q ) 2 = ω 2 0 + ∑ ij y ij q q i j with ω 0 = √ ( A + B /R /D 2 ) and y ij = C /D ij + Q / R ij ( + R A/B 2 2 ) for the amplitude-only action. Notice that the velocities of the modes have a tensor structure in general, and can be anisotropic depending on the k -space structure of the SOC. We are interested only in the low-energy Goldstone mode at T = 0 since this mode is not damped for all U = 0 provided that the quasiparticle excitations E s k are gapped for all ∆ 0 = 0.

It turns out that the general expressions for the nonkinetic-expansion coefficients can simply be written as a sum over their well-known counterparts for the usual one-band problem, i.e., A = ∑ s k ∆ 2 0 / (4 E 3 s k ) , B = ∑ s k ξ s k / (8 E 3 s k ) , D = ∑ s k ξ 2 s k / (16 E 5 s k ) , and R = ∑ s k 1 / (16 E 3 s k ) are due entirely to intraband mechanisms [22-25]. Motivated by the recent findings [6, 7, 21], and to gain more physical insight, we divide the kineticexpansion coefficients into two distinct contributions Q ij = Q intra ij + Q inter ij and C ij = C intra ij + C inter ij , based on their physical origins, i.e., depending on whether intraband or interband processes are involved. For instance, a compact way to write their intraband (i.e., conventional) contributions is

$$Q _ { i j } ^ { \text{intra} } = \sum _ { \text{sk} } \frac { 1 } { 1 6 E _ { \text{sk} } ^ { 3 } } \frac { \partial \xi _ { \text{sk} } } { \partial k _ { i } } \frac { \partial \xi _ { \text{sk} } } { \partial k _ { j } }, \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \collections$$

$$C _ { i j } ^ { \text{intra} } = \sum _ { s k } ^ { \dots } \frac { 1 } { 1 6 E _ { s k } ^ { 3 } } \left ( 1 - \frac { 5 \Delta _ { 0 } ^ { 2 } \xi _ { s k } ^ { 2 } } { E _ { s k } ^ { 4 } } \right ) \frac { \partial \xi _ { s k } } { \partial k _ { i } } \frac { \partial \xi _ { s k } } { \partial k _ { j } }. \quad ( 9 )$$

We note that these terms can be transformed into their more familiar but more involved forms [22, 27] through first writing ∂ [( ∂ξ s k /∂k j ) /E 3 s k ] /∂k i = [ ∂ ξ 2 s k / ∂k ∂k ( i j )] /E 3 s k -3( ∂ξ s k /∂k i )( ∂ξ s k /∂k j ) ξ s k /E 5 s k , and then identifying ∑ k ξ s k ∂ [( ∂ξ s k /∂k j ) /E 3 s k ] /∂k i ≡ -∑ k ( ∂ξ s k /∂k i )( ∂ξ s k /∂k j ) /E 3 s k after an integration by parts, and hence the name conventional. On the other hand, a compact way to write their interband (i.e., geometric) contributions is

$$Q _ { i j } ^ { \text{inter} } = - \sum _ { s \mathbf k } \frac { d _ { \mathbf k } } { 8 s \xi _ { \mathbf k } E _ { s \mathbf k } } g _ { \mathbf k } ^ { i j }, \quad \quad ( 1 0 ) \quad \frac { \mu \cdots } { a _ { s } \, \text{ar} } \\ \sinh _ { i \cdots \cdots } ^ { a _ { s } }$$

$$C _ { i j } ^ { \text{inter} } = - \sum _ { s \mathbf k } ^ { \text{---} } \frac { d _ { \mathbf k } } { 8 s \xi _ { \mathbf k } E _ { s \mathbf k } } \left ( 1 + \frac { 2 \Delta _ { 0 } ^ { 2 } } { d _ { \mathbf k } ^ { 2 } } \right ) g _ { \mathbf k } ^ { i j }, \quad ( 1 1 ) ^ { \quad \stackrel { \infty } { x _ { 0 } ^ { \text{inter} } } } } \\ \sum _ { \mathcal { Q } _ { \mathcal { Q } } }$$

where g ij k = (1 / 2) lim q → 0 ∂ 2 ( ̂ d + · ̂ d -) / ∂q ∂q ( i j ) or equivalently g ij k = ( ∂ ̂ d k /∂k i ) · ( ∂ ̂ d k /∂k j ) / 2 is the total quantum metric tensor of the helicity bands [6, 7]. Alternatively, the quantum metric can be written as g ij k = [ ∑ glyph[lscript] ( ∂d /∂k glyph[lscript] k i )( ∂d /∂k glyph[lscript] k j ) -( ∂d /∂k k i )( ∂d /∂k k j )] / (2 d 2 k ) , reducing to g ij k = α α i j ( d 2 k δ ij -d i k d j k ) / (2 d 4 k ) for the specific case when d i k = α k i i . Since the interband contributions are directly controlled by the geometry of the underlying quantum states in the parameter ( k ) space, we refer to them as the geometric ones.

All of the coefficients A B , , C ij , etc., are derived at T = 0 for arbitrary U and d k , and can be used both in two and three dimensions. Next we present their specific application to three systems that are of immediate experimental and/or theoretical interest: (i) a 3D Fermi gas with a Weyl SOC, (ii) a 3D Fermi gas with a Rashba SOC, and (iii) a 2D Fermi gas with a Rashba SOC.

## III. NUMERICAL RESULTS

In this paper, our primary interest is in the cooperation and/or competition between the conventional intraband contribution x intra ij = Q intra ij / R ( + B /A 2 ) and the geometric interband one x inter ij = Q inter ij / R ( + B /A 2 ) on the velocity of the Goldstone mode. For this purpose, one first needs to obtain a self-consistent solution for ∆ 0 and µ by coupling the saddle-point condition and the saddle-point number equation, and then plug these solutions into the coefficients.

FIG. 1: (color online) A 3D Fermi gas with Weyl SOC is mapped at T = 0 in the plane of two-body scattering length a s and SOC strength α . Square of the sound velocity x /v 0 2 F is shown (upper panel) along with its geometric contribution x inter 0 /v 2 F (lower panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00818v2-figure-913f82c3-11e8-4a28-abfc-a9241155e7ab.jpg)

For instance, in Fig. 1, we show x ij = x δ 0 ij for a 3D Fermi gas with a Weyl SOC, where d k = α , k and therefore, x 0 = x intra 0 + x inter 0 is isotropic in the entire q space with δ ij a Kronecker-delta. In accordance with the coldatom literature, the two-body s -wave scattering length a s in vacuum is substituted in this figure for U via the usual relation V/U = -mV/ (4 πa s ) + ∑ k 1 / (2 glyph[epsilon1] k ) , where V is the volume of the system. Furthermore, we choose N = k 3 F V/ (3 π 2 ) for the number of particles, and use the Fermi momentum k F and the Fermi velocity v F = k F /m as the relevant length and velocity scales in our numerical calculations. This is such that increasing U from 0 changes the dimensionless parameter 1 / k ( F a s ) continuously from -∞ to 0 to + ∞ in the BCS-BEC crossover.

In the weak-SOC limit when α/v F glyph[lessmuch] 1, we recover the usual problem where x /v 0 2 F is known to be a monotonically decreasing function of 1 / k ( F a s ) with the wellknown limits, i.e., 1 / 3 in the BCS regime and k F a / s (3 π ) in the BEC regime [22, 23]. Note that since the geometric contribution x inter 0 is originating from the interband pro-

cesses, it plays a negligible role when α/v F → 0. On the other hand, for intermediate α/v F values, we find that x /v 0 2 F is a nonmotonic function of 1 / k ( F a s ), exhibiting a noticeable peak in the crossover region around the unitarity [24, 25]. This result is quite a surprise given that all of the ground-state properties are found to be monotonic functions of 1 / k ( F a s ) for the usual BCS-BEC crossover problem with s -wave interactions [22]. In comparison to the dominant contribution from x intra 0 , our study reveals that, despite being small, the contribution from x inter 0 is sizeable enough to cause such a significant difference in x 0 . Thus, even though a direct observation of x inter 0 may not be possible, realization of a nonmonotonic x 0 may be considered as the ultimate evidence for its subleading role.

FIG. 2: (color online) A 3D Fermi gas with Rashba SOC is mapped at T = 0 in the plane of two-body scattering length a s and SOC strength α . Square of the sound velocities x ⊥ /v 2 F and x zz /v 2 F are shown (upper and lower panels) along with the geometric contribution x inter ⊥ /v 2 F (middle panel) of the former.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00818v2-figure-88d84715-449c-4c06-bd5f-874d1caf17cd.jpg)

Similarly, in Fig. 2, we show x ij = x ii δ ij for a 3D Fermi gas with a Rashba SOC, where d k = α k ⊥ with k ⊥ = k x ̂ x + k y ̂ y , and therefore, x ⊥ = x intra ⊥ + x inter ⊥ is still isotropic in q x q y plane but it differs from that x zz = x intra zz of the q z direction with x inter zz = 0. Note that since g iz k = 0 for a Rashba SOC, all of the matrix elements x inter iz trivially vanish, leading to the aforementioned anisotropy. For intermediate α/v F values, we again find that x ⊥ /v 2 F is a nonmonotonic function of 1 / k ( F a s ), exhibiting a noticeable peak in the crossover region around the unitarity. However, x zz /v 2 F turns out to be a monotonic one no matter what α/v F is. Thus, the dramatic difference between x ⊥ and x zz may be used as an ultimate evidence for the presence of an interband contribution to the former, reflecting the geometry of the underlying quantum states.

FIG. 3: (color online) A 2D Fermi gas with Rashba SOC is mapped at T = 0 in the plane of two-body binding energy glyph[epsilon1] b and SOC strength α . Square of the sound velocity x /v 0 2 F is shown (upper panel) along with its geometric contribution x inter 0 /v 2 F (lower panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00818v2-figure-cf29b2d6-ba53-4410-bbb3-bba3ea48fbf7.jpg)

For completeness, in Fig. 3, we show x ij = x δ 0 ij for a 2D Fermi gas with a Rashba SOC, where d k = α , k and therefore, x 0 = x intra 0 + x inter 0 is isotropic in the entire q space. In accordance with the cold-atom literature, the two-body binding energy ε b ≥ 0 in vacuum is substituted in this figure for U via the usual relation A/U = ∑ k 1 / (2 glyph[epsilon1] k + ε b ) , where A is the area of the system. Furthermore, we choose N = k 2 F A/ (2 π , ) and use the Fermi energy glyph[epsilon1] F = k 2 F / (2 m ) as the relevant energy scale in our numerical calculations. This is such that increasing U from 0 increases ε /glyph[epsilon1] b F continuously from 0 to ∞ in the BCS-BEC crossover. In the weak-SOC limit when α/v F glyph[lessmuch] 1, we again recover the usual problem where x /v 0 2 F is known to be 1 / 2 for all ε /glyph[epsilon1] b F , i.e., it remains a constant in the entire BCS-BEC crossover [23]. In contrast to the 3D results shown in Figs. 1 and 2 where the geometric contribution x inter 0 plays important roles only in the crossover region with a rapid decay in the BCS and BEC regimes, here x inter 0 retains its strength for most ε /glyph[epsilon1] b F values except for the extreme BCS limit when ε /glyph[epsilon1] b F glyph[lessmuch] 1. Even then we find that the fraction x inter 0 /x 0 does not grow much, and saturates around a quarter.

## IV. CONCLUSIONS

In summary, here we considered a two-component Fermi gas with an arbitrary SOC, and analyzed its collective excitations in the BCS-BEC crossover. Motivated by the recent findings [6, 7], we divided the velocity of the Goldstone mode into two distinct contributions based

on the physical mechanisms involved. While the conventional (intraband) contribution is controlled by the k -space derivatives of the helicity (one-body) spectrum, the geometric (interband) one is controlled by the k -space derivatives of the helicity eigenstates or the socalled quantum metric tensor characterizing the geometry of the underlying quantum states. We presented extensive numerical calculations for the Weyl and Rashba SOCs revealing that, despite being much smaller than the conventional one, the geometric contribution is solely responsible for the nonmonotonic evolution of the mode velocity in the BCS-BEC crossover.

These findings for the Goldstone mode in spin-orbit coupled Fermi gases are in complete agreement with our recent works that the quantum metric tensor has a subleading control over all those SF properties that depend explicitly on the effective-mass tensor of the corresponding (twoor many-body) bound state [6, 7]. We note that whether the quantum-metric contribution can be isolated from the total is far from being clear with the current experimental techniques involving spin-orbit coupling. However, in the same spirit, given that the starting Hamiltonian Eq. (1) is quite generic, our analytical expressions for the Goldstone mode may also find direct applications in other contexts that exhibit a two-band band structure. In particular, it is desirable to study the Goldstone mode in a narrow- or a flat-band system [28, 29], for which the geometric contribution is expected to take the lead. As a final remark, our work suggests that an analogous contribution to the collective excitations must be present in many other multi-band systems including the twisted bilayer graphene.

## Acknowledgments

This work is supported by the funding from T ¨ B ˙ ITAK U Grant No. 1001-118F359.

## Appendix A: Effective-Action Approach

The imaginary-time functional path-integral formalism allows a systematic derivation of the effective action in the form of a power-series expansion of the fluctuations of the order parameter around its saddle-point value [2225]. In this approach, to decouple the interaction term that is quartic in the fermionic degrees of freedom, one first uses a Hubbard-Stratanovich transformation at the expense of introducing a bosonic field ∆ , and then inq tegrates out the remaining terms that are quadratic in the fermionic degrees of freedom. Finally, by decomposing ∆ q = ∆ +Λ 0 q in terms of a q -independent stationary field ∆ 0 and q -dependent fluctuations around it, one may in principle obtain an effective action at the desired order in Λ . q Here we include the first nontrivial term and obtain the effective-Gaussian action S Gauss = S 0 + S 2 , as the first-order term S 1 trivially vanishes thanks to the saddle-point condition given in the main text.

For our interests, the zeroth-order action in Λ q can be written as S 0 = ∆ 2 0 / TU ( ) + ∑ s k ξ s k / (2 T ) -(1 / 2) ∑ k ln[det( G -1 k /T )] , where the collective index k = ( k , i ω glyph[lscript] ) denotes both the particle momentum k and the fermionic Matsubara frequency ω glyph[lscript] = (2 glyph[lscript] +1) πT [22-25]. Here, G -1 k = i ω glyph[lscript] 1 -H 0 k is the inverse Green's function for the saddle-point Hamiltonian H 0 k , i.e.,

$$\text{th our} _ { \text{a sub-} } \quad G _ { k } ^ { - 1 } = \left [ \begin{array} { c c } ( i \omega _ { \ell } - \xi _ { k } ) \sigma _ { 0 } - \mathrm d _ { k } \cdot \sigma & i \Delta _ { 0 } \sigma _ { y } \\ - i \Delta _ { 0 } \sigma _ { y } & ( i \omega _ { \ell } + \xi _ { k } ) \sigma _ { 0 } - \mathrm d _ { k } \cdot \sigma ^ { * } \end{array} \right ],$$

determining the quasiparticle and quasihole energies through det G -1 k = [(i ω glyph[lscript] ) 2 -E 2 + , k ][(i ω glyph[lscript] ) 2 -E 2 -, k ] = 0 , leading to E s k = √ ξ 2 s k +∆ for the former. Performing 2 0 the summation over the fermionic Matsubara frequency leads to the saddle-point action given in Eq. (2) of the main text.

The second-order term in Λ q can be written as S 2 = ∑ q | Λ q | 2 / TU ( )+(1 4)Tr / ∑ kq G Σ G k q k + q Σ -q , where Tr denotes a trace over the particle/hole and spin sectors, and the collective index q = ( q , i ν n ) denotes both the pair momentum q and the bosonic Matsubara frequency ν n = 2 πnT [23]. Here, the matrix elements of the Green's function G k can be written as

$$G _ { k } ^ { 1 1 } = \frac { 1 } { 2 } \sum _ { s } \frac { \mathrm i \omega _ { \ell } + \xi _ { s k } } { ( \mathrm i \omega _ { \ell } ) ^ { 2 } - E _ { s k } ^ { 2 } } \left ( \sigma _ { 0 } + s \widehat { \mathrm d } _ { k } \cdot \sigma \right ), \quad ( \mathrm A 1 )$$

$$G _ { k } ^ { 2 2 } = \frac { 1 } { 2 } \sum _ { s } \frac { \text{i} \omega _ { \ell } - \xi _ { s k } } { ( \text{i} \omega _ { \ell } ) ^ { 2 } - E _ { s k } ^ { 2 } } \left ( \sigma _ { 0 } - s \widehat { \mathbf d } _ { \mathbf k } \cdot \sigma ^ { * } \right ), \quad ( \text{A2} )$$

$$G _ { k } ^ { 1 2 } = \frac { 1 } { 2 } \sum _ { s } \frac { - \mathrm i \Delta _ { 0 } \sigma _ { y } } { ( \mathrm i \omega _ { \ell } ) ^ { 2 } - E _ { s \mathrm k } ^ { 2 } } \left ( \sigma _ { 0 } - s \widehat { \mathrm d } _ { \mathrm k } \cdot \sigma ^ { * } \right ), \quad ( \mathrm A 3 )$$

$$G _ { k } ^ { 2 1 } = \frac { 1 } { 2 } \sum _ { s } \frac { \mathrm i \Delta _ { 0 } \sigma _ { y } } { ( \mathrm i \omega _ { \ell } ) ^ { 2 } - E _ { s k } ^ { 2 } } \left ( \sigma _ { 0 } + s \widehat { \mathrm d } _ { \mathbf k } \cdot \sigma \right ), \quad ( \Lambda 4 )$$

and those of the fluctuations Σ q are Σ 11 q = Σ 22 q = 0 , Σ 12 q = iΛ q σ y and Σ 21 q = -iΛ ∗ -q σ y [30, 31]. After a lengthy but a straightforward summation over the fermionic Matsubara frequency, we find S 2 = [1 / (2 T )] ∑ q ¯ Λ M ¯ † q q Λ q , where ¯ Λ † q = (Λ ∗ q Λ -q ) and M q is the inverse fluctuation propagator [22-25]. Simplifying the notation by denoting ξ s, k + q / 2 by ξ ; ξ s , ′ -k + q / 2 by ξ ′ ; E s, k + q / 2 by E ; E s , ′ -k + q / 2 by E ′ ; and ̂ d ± k + q / 2 by ̂ d ± , and defining the functions u 2 = (1 + ξ/E / ) 2; u ′ 2 = (1+ ξ /E ′ ′ ) / 2; v 2 = (1 -ξ/E / ) 2; v ′ 2 = (1 -ξ /E ′ ′ ) / 2; f = 1 ( / e E/T +1); and f ′ = 1 ( / e E /T ′ +1), a compact way

to write the matrix elements of M q is [22, 23]

$$\begin{array} { l l } \text{to write the matrix elements of } M _ { q } \text{ is [22, 23]} & \text{ones i} \\ M _ { q } ^ { 1 1 } = \frac { 1 } { U } + \frac { 1 } { 4 } \sum _ { s s ^ { \prime } \mathbf k } \left ( 1 - s s ^ { \prime } \widehat { d } _ { + } \cdot \widehat { d } _ { - } \right ) & \text{as a unit and } \theta \\ \lambda _ { q } ^ { \prime } \text{ an} \\ \times \left [ ( 1 - f - f ^ { \prime } ) \left ( \frac { u ^ { 2 } u ^ { \prime 2 } } { i \nu _ { n } - E - E ^ { \prime } } - \frac { v ^ { 2 } v ^ { \prime 2 } } { i \nu _ { n } - E + E ^ { \prime } } \right ) \\ + ( f - f ^ { \prime } ) \left ( \frac { v ^ { 2 } u ^ { \prime 2 } } { i \nu _ { n } + E - E ^ { \prime } } - \frac { u ^ { 2 } v ^ { \prime 2 } } { i \nu _ { n } - E + E ^ { \prime } } \right ) \right ], & ( A 5 ) \\ \\ M _ { q } ^ { 1 2 } = \frac { 1 } { 4 } \sum _ { s s ^ { \prime } \mathbf k } \left ( 1 - s s ^ { \prime } \widehat { d } _ { + } \cdot \widehat { d } _ { - } \right ) \\ \times \left [ ( 1 - f - f ^ { \prime } ) \left ( \frac { \ u v u ^ { \prime } v ^ { \prime } } { i \nu _ { n } + E + E ^ { \prime } } - \frac { \ u v u ^ { \prime } v ^ { \prime } } { i \nu _ { n } - E - E ^ { \prime } } \right ) & \text{where} \\ + ( f - f ^ { \prime } ) \left ( \frac { \ u v u ^ { \prime } v ^ { \prime } } { i \nu _ { n } + E - E ^ { \prime } } - \frac { \ u v u ^ { \prime } v ^ { \prime } } { i \nu _ { n } - E + E ^ { \prime } } \right ) \right ]. & ( A 6 ) & \text{and } u \text{given} \\ \text{The remaining elements are related via } M _ { q } ^ { 2 2 } = M _ { q } ^ { 1 1 } \text{ and } \text{ due } \text{since}. \end{array}$$

The remaining elements are related via M 22 q = M 11 -q and M 21 q = M 12 q . We note that while M 12 q is even both under q → -q and i ν n → -i ν n , M 11 q is even only under q → -q .

While calculating the collective excitations, it is useful to represent Λ q = α e q i γ q in terms of real functions α q and γ q , and associate λ q = √ 2 α q cos( γ q ) with the amplitude degrees of freedom and θ q = √ 2 α q sin( γ q ) with the phase

- [1] J. P. Provost and G. Vallee, Riemannian structure on manifolds of quantum states, Commun. Math. Phys. 76 , 289 (1980).
- [2] M. V. Berry, The quantum phase, five years after in Geometric Phases in Physics, edited by A. Shapere and F. Wilczek (World Scientific, Singapore, 1989).
- [3] D. J. Thouless, Topological Quantum Numbers in Nonrelativistic Physics, (World Scientific, Singapore,1998)
- [4] D. Xiao, M.-C. Chang, and Q. Niu, Berry phase effects on electronic properties, Rev. Mod. Phys. 82 ,1959 (2010).
- [5] X.-L. Qi and S.-C. Zhang, Topological insulators and superconductors, Rev. Mod. Phys. 83 , 1057 (2011).
- [6] M. Iskin, Exposing the quantum geometry of spin-orbit coupled Fermi superfluids, Phys. Rev. A 97 , 063625 (2018).
- [7] M. Iskin, Quantum metric contribution to the pair mass in spin-orbit coupled Fermi superfluids, Phys. Rev. A 97 , 033625 (2018).
- [8] R. Resta, The insulating state of matter: A geometrical theory, Eur. Phys. J. B 79 , 121 (2011).
- [9] X. Tan, D.-W. Zhang, Z. Yang, J. Chu, Y.-Q. Zhu, D. Li, X. Yang, S. Song, Z. Han, Z. Li, Y. Dong, H.-F. Yu, H. Yan, S.-L. Zhu, and Y. Yu, Experimental Measurement of the Quantum Metric Tensor and Related Topological Phase Transition with a Superconducting Qubit, Phys. Rev. Lett. 122 , 210401 (2019).
- [10] M. Yu, P. Yang, M. Gong, Q. Cao, Q. Lu, H. Liu, M. B. Plenio, F. Jelezko, T. Ozawa, N. Goldman, S. Zhang, and J. Cai, Experimental measurement of the complete quantum geometry of a solid-state spin system, arXiv:1811.12840 (2018).

ones in the small γ q limit. This procedure corresponds to a unitary transformation Λ q = ( λ q +i θ q ) / √ 2 , where λ q and θ q are real functions. Furthermore, assuming λ -q = λ ∗ q and θ -q = θ ∗ q , we find

$$\sum _ { 0 5 } ^ { \dagger } S _ { 2 } = \frac { 1 } { 2 T } \sum _ { q } \Theta _ { q } ^ { \dagger } \left ( \begin{array} { c c } M _ { q, E } ^ { 1 1 } + M _ { q } ^ { 1 2 } & \mathrm i M _ { q, O } ^ { 1 1 } \\ - \mathrm i M _ { q, O } ^ { 1 1 } & M _ { q, E } ^ { 1 1 } - M _ { q } ^ { 1 2 } \end{array} \right ) \Theta _ { q },$$

where we denote Θ † q = ( λ ∗ q θ ∗ q ), and separate M 11 q = M 11 q,E + M 11 q,O in terms of an even function M 11 q,E = ( M 11 q + M 22 q ) / 2 in i ν n and an odd one M 11 q,O = ( M 11 q -M 22 q ) / 2 , and use M 12 q = M 21 q [22-25]. The phase-only action given in Eq. (3) is derived by integrating out the amplitude fields, and vice versa for the Eq. (4). We note that since the quasiparticle-quasihole terms with the prefactor ( f -f ′ ) have the usual Landau singularity for q → ( 0 , 0) causing the collective modes to decay into the two-quasiparticle continuum, a small q expansion is possible only in two cases: (i) just below the critical SF transition temperature provided ∆ 0 → 0 glyph[lessmuch] i ν n → ω [7], and (ii) at T = 0. In this work, we are interested in the latter case [21].

- [11] S. Peotta and P. T¨rm¨ a, o Superfluidity in topologically nontrivial flat bands, Nat. Commun. 6 , 8944 (2015).
- [12] L. Liang, T. I. Vanhala, S. Peotta, T. Siro, A. Harju, and P. T¨rm¨ a, Band geometry, Berry curvature, and supero fluid weight, Phys. Rev. B 95 , 024515 (2017).
- [13] X. Hu, T. Hyart, D. I. Pikulin, and E. Rossi, Geometric and conventional contribution to superfluid weight in twisted bilayer graphene, arXiv:1906.07548 (2019).
- [14] A. Julku, T. J. Peltonen, L. Liang, T. T. Heikkil¨, a and P. T¨ orm¨, a Superfluid weight and BerezinskiiKosterlitz-Thouless transition temperature of twisted bilayer graphene, arXiv:1906.06313 (2019).
- [15] F. Xie, Z. Song, B. Lian, and B. A. Bernevig, Topology-Bounded Superfluid Weight In Twisted Bilayer Graphene, arXiv:1906.02213 (2019).
- [16] P. Wang, Z. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. Zhai, and J. Zhang, Spin-orbit coupled degenerate Fermi gases, Phys. Rev. Lett. 109 , 095301 (2012).
- [17] L. W. Cheuk, A. T. Sommer, Z. Hadzibabic, T. Yefsah, W. S. Bakr, and M. W. Zwierlein, Spin-Injection Spectroscopy of a Spin-Orbit Coupled Fermi Gas, Phys. Rev. Lett. 109 , 095302 (2012).
- [18] R. A. Williams, M. C. Beeler, L. J. LeBlanc, K. Jim´neze Garc´ ıa, and I. B. Spielman, Raman-induced interactions in a single-component Fermi gas near an s-wave Feshbach resonance, Phys. Rev. Lett. 111 , 095301 (2013).
- [19] L. Huang, Z. Meng, P. Wang, P. Peng, S.-L. Zhang, L. Chen, D. Li, Q. Zhou, and J. Zhang, Experimental realization of a two-dimensional synthetic spin-orbit coupling in ultracold Fermi gases, Nat. Phys. 12 , 540 (2016).
- [20] Z. Meng, L. Huang, P. Peng, D. Li, L. Chen, Y. Xu,

- C. Zhang, P. Wang, and J. Zhang, Experimental Observation of a Topological Band Gap Opening in Ultracold Fermi Gases with Two-Dimensional Spin-Orbit Coupling, Phys. Rev. Lett. 117 , 235304 (2016).
- [21] This work is complementary to our recent works. In [6], we used a linear response theory and identified the SFdensity tensor from the response of the thermodynamic potential to an infinitesimal superfluid flow. In [7], we derived the time-dependent Ginzburg-Landau theory near the critical temperatures, and identified the effectivemass tensor of the Cooper pairs by making an analogy with the Gross-Pitaevskii theory of the corresponding molecular Bose gas.
- [22] J. R. Engelbrecht, M. Randeria, and C. A. R. S´ de Melo, a BCS to Bose crossover: Broken-symmetry state, Phys. Rev. B 55 , 15153 (1997)
- [23] M. Iskin and C. A. R. S´ de Melo, BCS-BEC crossover of a collective excitations in two-band superfluids, Phys. Rev. B 72 , 024512 (2005).
- [24] L. He and X.-G. Huang, BCS-BEC crossover in threedimensional Fermi gases with spherical spin-orbit coupling, Phys. Rev. B 86 , 014511 (2012).
- [25] J. P. Vyasanakere and V. B. Shenoy, Collective excitations, emergent Galilean invariance, and boson-boson interactions across the BCS-BEC crossover induced by a synthetic Rashba spin-orbit coupling, Phys. Rev. A 86 ,

053617 (2012).

- [26] To be more precise, one needs to derive the characteristic equation up to fourth order in the expansion, e.g., β ω 5 4 + ∑ ijkl β ijkl 4 q q i j q k q l + ∑ ij β ij 3 ω q q 2 i j + β ω 2 2 + ∑ ij β ij 1 q q i j = 0 , such that x ij = -β ij 1 /β , 2 ω 2 0 = -β /β 2 5 and y ij = -β ij 3 /β 5 + β ij 1 /β . 2 While our quadratic expansion fully determines x ij , we neglect the additional corrections to ω 2 0 and y ij coming from the higher order terms.
- [27] The familiar expressions for the one-band coefficients can be written as [22] Q ij = Q intra ij = ∑ k [ ξ k ∂ ξ 2 k / ∂k ∂k ( i j ) -( ∂ξ k /∂k i )( ∂ξ k /∂k j )(1 -3∆ 2 0 /E 2 k )] / (8 E 3 k ) and C ij = C intra ij = ∑ k [ ξ k (1 -3∆ 2 0 /E 2 k ) ∂ ξ 2 k / ∂k ∂k ( i j ) -( ∂ξ k /∂k i )( ∂ξ k /∂k j )(1 -10∆ 2 0 ξ 2 k /E 4 k )] / (8 E 3 k ) , and our intraband coefficients recover them in the absence of a SOC when d k → 0.
- [28] M. Iskin, Superfluid stiffness for the attractive Hubbard model on a honeycomb optical lattice, Phys. Rev. A 99 , 023608 (2019).
- [29] M. Iskin, Origin of flat-band superfluidity on the Mielke checkerboard lattice, Phys. Rev. A 99 , 053608 (2019).
- [30] V. M. Edelstein, Magnetoelectric Effect in Polar Superconductors, Phys. Rev. Lett. 75 , 2004 (1995).
- [31] L. P. Gorkov and E. I. Rashba, Superconducting 2D System with Lifted Spin Degeneracy: Mixed Singlet-Triplet State, Phys. Rev. Lett. 87 , 037004 (2001).
M. Iskin
Geometric contribution to the Goldstone mode in spin-orbit coupled Fermi superfluids
The so-called quantum metric tensor is a band-structure invariant whose
measure corresponds to the quantum distance between nearby states in the
Hilbert space, characterizing the geometry of the underlying quantum states. In
the context of spin-orbit coupled Fermi gases, we recently proposed that the
quantum metric has a partial control over all those superfluid properties that
depend explicitly on the mass of the superfluid carriers, i.e., the
effective-mass tensor of the corresponding (two- or many-body) bound state.
Here we scrutinize this finding by analyzing the collective phase and amplitude
excitations at zero temperature. In particular to the Goldstone mode, we
present extensive numerical calculations for the Weyl and Rashba spin-orbit
couplings, revealing that, despite being small, the geometric contribution is
solely responsible for the nonmonotonic evolution of the sound velocity in the
BCS-BEC crossover.