---
title: Superfluid Weight Tensor
url: https://www.emergentmind.com/topics/superfluid-weight-tensor
type: topic
---

# Superfluid Weight Tensor

The superfluid weight tensor is a fundamental observable in superconductivity, characterizing the phase stiffness of a superconductor and determining its response to electromagnetic fields and phase gradients. It plays a central role in phenomena such as the Meissner effect, London penetration depth, and the Berezinskii-Kosterlitz-Thouless (BKT) transition in two-dimensional systems. The superfluid weight incorporates both conventional contributions associated with band dispersion and geometric contributions arising from the quantum geometry of electronic or bosonic Bloch states, including the quantum metric and, in nontrivial cases, non-Abelian generalizations. Recent work firmly establishes its decomposition, symmetry constraints, and unique physical consequences in both conventional and unconventional superconductors, multiorbital and flat-band systems, as well as quasicrystals.

## 1. Definition and Physical Interpretation

The superfluid weight tensor, commonly denoted $D_{ij}$, quantifies the second derivative of the equilibrium free energy $F$ of a superconductor with respect to a uniform vector potential $\mathbf{A}$ (or equivalently a phase twist or flux insertion) at zero field:

\[
D_{ij} = \frac{1}{V} \frac{\partial^2 F}{\partial A_i \partial A_j}\biggr|_{\mathbf{A}=0}
\]
where $V$ is the system volume. Equivalently, in linear response,
\[
\langle J_i \rangle = D_{ij} A_j
\]
links the induced supercurrent to the applied vector potential. In the London gauge, $D_{ij}$ directly determines the tensorial inverse squared penetration depth: $(\lambda^{-2})_{ij} = \mu_0 D_{ij}$, controlling the Meissner screening of magnetic fields [2204.12163, 2603.10955].

Alternative, but equivalent, formulations use the Kubo formula: $D_{ij}$ is the static, long-wavelength limit of the current-current response kernel,
\[
D_{ij} = K_{ij}(\mathbf{q} \to 0, \omega=0)
\]
with $K_{ij}$ comprising diamagnetic and paramagnetic terms, distinguished by their microscopic origins [1610.01803, 2505.17349].

In bosonic systems such as flat-band Bose-Einstein condensates (BECs), analogous definitions arise from imposing a long-wavelength phase twist and expanding the free energy, including necessary corrections for condensate-density and chemical-potential changes induced by the twist [2210.11906].

## 2. Decomposition: Conventional and Geometric Contributions

A general multiband formalism shows that at the mean-field level, $D_{ij}$ decomposes as:
\[
D_{ij} = D_{ij}^\text{conv} + D_{ij}^\text{geom}
\]
where:

- $D_{ij}^\text{conv}$ (conventional term) is controlled by the band curvature (Fermi velocity and its derivatives) and vanishes if the bands are perfectly flat,
- $D_{ij}^\text{geom}$ (geometric term) arises from quantum geometric properties of the Bloch states—specifically the quantum metric tensor—and can remain finite even for exactly flat bands provided the quantum metric is nonzero [1610.01803, 2110.14663].

The microscopic origin is elucidated by expanding the current operator matrix elements in the Bloch basis. The diagonal (intraband) part corresponds to $D_{ij}^\text{conv}$, while the off-diagonal (interband) matrix elements yield $D_{ij}^\text{geom}$. Closed-form expressions include
\[
D_{ij}^\text{conv} = \frac{1}{V}\sum_{m,\mathbf{k}} \frac{\Delta^2}{\left(\epsilon_{m\mathbf{k}}^2 + \Delta^2\right)^{3/2}} \left(\partial_i \epsilon_{m\mathbf{k}}\right)\left(\partial_j \epsilon_{m\mathbf{k}}\right)
\]
and
\[
D_{ij}^\text{geom} = \frac{2}{V} \sum_{m \neq n, \mathbf{k}} \frac{\Delta^2}{\sqrt{\epsilon_{m\mathbf{k}}^2+\Delta^2}\;\left(\sqrt{\epsilon_{n\mathbf{k}}^2+\Delta^2}-\sqrt{\epsilon_{m\mathbf{k}}^2+\Delta^2}\right)} \frac{g_{mn}^{ij}(\mathbf{k})}{\epsilon_{n\mathbf{k}}+\epsilon_{m\mathbf{k}}}
\]
where $g_{mn}^{ij}(\mathbf{k})$ are components of the (non-Abelian) quantum metric tensor [2603.10955, 2501.16965].

In systems with time-reversal symmetry, the geometric term can be expressed as an integral over the trace of the quantum metric:
\[
D_{ij}^\text{geom} = -\Delta^2 \int_{\text{BZ}}\! \frac{d^d k}{(2\pi)^d} \operatorname{Tr} R_{ij}(\mathbf{k}) \frac{n_F(\epsilon_{\mathbf{k}}) - n_F(-\epsilon_{\mathbf{k}})}{\epsilon_{\mathbf{k}}}
\]
with $R_{ij}$ the (non-Abelian) real quantum metric, and $n_F$ the Fermi distribution [2501.16965].

This decomposition persists for bosonic BECs, where $D_{ij}^{\text{geom}} \propto U n_0 g_{ij}(\mathbf{k}_c)$ (with $g_{ij}$ the Fubini-Study metric at the condensate momentum) [2210.11906], as well as in superconductors with unconventional (momentum-dependent) pairing, although additional functional terms appear in the latter [2505.09249].

## 3. Quantum Geometry, Topological Bounds, and Flat Bands

The geometric contribution is fundamentally tied to the quantum metric of Bloch bands. In paradigmatic flat-band systems with uniform pairing and strong isolation from other bands, the geometric term dominates and can be written as:
\[
D_{ij} = \mathcal{C} \Delta^2 \int_{\text{BZ}}\! d^d k\; g_{ij}(\mathbf{k})
\]
where $\mathcal{C}$ depends on microscopic details [2110.14663, 1906.02213].

Several rigorous lower bounds have been established:
- In 2D, the integrated quantum metric is bounded below by the modulus of the Berry curvature, and, for topologically nontrivial bands, by the Wilson-loop winding number or Chern number:
\[
\int_{\text{BZ}} \operatorname{Tr} g(k)\,dk \ge 2\pi |C|
\]
leading to a nonzero lower bound on $D_{ij}$, and thus on the BKT transition temperature $T_{BKT}$ [1906.02213, 2110.14663, 2405.11260]. Bands with obstructed Wannier centers—i.e., Wannier functions displaced from atomic sites—can be strictly identified using real space invariants (RSIs) and are guaranteed a positive superfluid weight even when Berry curvature vanishes identically [2110.14663].

For composite bands (multiple nearly flat and isolated bands), the geometric contribution includes both intra- and inter-band quantum metric terms. Provided the total Chern number is nonzero, a topological lower bound on $D_{ij}$ persists. Detailed perturbative expressions for these "composite quantum metrics" and their dependence on band representatives are available [2405.11260].

Notably, in systems where one or more of isolation, flatness, or symmetry conditions are relaxed (e.g., multiorbital, spin-polarized, or three-dimensional systems), these bounds can be lost, and $D_{ij}$ can vanish or become negative within certain parameter regimes [2308.13400].

## 4. Anisotropy, Symmetry, and Extensions

The superfluid weight tensor is intrinsically anisotropic in systems lacking full rotational or crystalline symmetry, with $D_{ij}$ elements reflecting the underlying symmetry group. For instance, orthorhombic or lower-symmetry crystals possess distinct diagonal and, potentially, off-diagonal components. The tensor structure is inherited from both the microscopic band structure and the symmetry of the superconducting order parameter [2603.10955, 2210.11906].

Pairing symmetry profoundly influences the tensorial structure. For $s$-wave (on-site) pairing, especially in the presence of $C_3$ symmetry, $D_{ij}$ is isotropic, while in nematic (nearest-neighbor) or unconventional (e.g., chiral $d$-wave) pairing, anisotropy and nonvanishing off-diagonal elements emerge—an effect directly observable in the kinetic inductance or penetration depth anisotropy [1906.06313, 2505.09249].

For systems with degenerate bands, the appropriate quantum geometric tensor is non-Abelian, with all physically relevant quantities depending only on the trace of the real (symmetric) part. A crucial result is that a vanishing net Chern number does not guarantee a vanishing geometric contribution; in realistic systems, the metric can be responsible for a significant fraction of $D_{ij}$, as seen in DFT calculations for MoS$_2$ and TiSe$_2$ [2501.16965].

## 5. Inhomogeneity, Quasicrystals, and Real-Space Formulation

Spatial inhomogeneity (e.g., vortex lattices, pinning, strong disorder, quasicrystal structures) requires more sophisticated formalisms, as translational symmetry is lost and momentum-space geometry is not directly accessible. In these situations, the superfluid weight retains a definition in terms of boundary-twist or real-space Kubo formulas:
\[
D_{ij} = \frac{1}{V} \frac{\partial^2 F}{\partial \Phi_i \partial \Phi_j}\biggr|_{\Phi=0}
\]
where $\Phi_i$ are fluxes inserted through the holes of a torus or ring [2507.20540, 2606.09989]. A key development is the introduction of a "flux-space quantum metric", defined via derivatives with respect to fluxes:
\[
g_{ij} = \Re \langle \partial_i \Psi(0) | [1 - |\Psi(0)\rangle\langle\Psi(0)|] | \partial_j \Psi(0) \rangle
\]
with $|\Psi(\Phi)\rangle$ the many-body ground state at flux $\Phi$.

In the flat-band, weak-coupling limit, the conventional component vanishes and the geometric term, proportional to the averaged local Wannier spread or its fluctuations, dictates the entire superfluid weight [2507.20540].

For quasicrystals hosting critical states, which are neither extended nor exponentially localized, the geometric contribution to $D_{ij}$ dominates at all relevant couplings. Real-space Kubo and momentum-space approaches reveal that $D_{ij}$ strictly traces the critical nature of the wavefunctions. This is an essential demonstration that geometric superfluidity extends far beyond periodic crystals [2606.09989].

## 6. Experimental Manifestations and Numerical Evaluation

The superfluid weight determines macroscopic observables such as the London penetration depth $\lambda_L$ and the BKT transition temperature in two dimensions. Explicitly,
\[
\lambda_L^{-2} = \mu_0 D
\]
(for isotropic $D$), and for a 2D superconductor,
\[
k_B T_{BKT} = \frac{\pi}{8} \sqrt{\det D_{ij}(T_{BKT})}
\]
[2603.10955, 1906.02213, 1906.06313]. Experimentally, in conventional superconductors with wide bands such as Al, Pb, Nb, the geometric contribution is negligible (3–4 orders of magnitude smaller than the band-curvature term), and numerical evaluation using density functional theory band structures achieves agreement within 10–20% of measured London depths [2603.10955]. In topological or flat-band candidates (e.g., twisted bilayer graphene, copper-doped apatite, $\alpha$–𝒯$_3$ lattice), the geometric term dominates or contributes significantly, with its magnitude tunable via band-structure engineering, electron filling, and interaction strength [2512.09901, 2308.13400].

For detailed numerical evaluation, dense $k$-grid sampling, analytic derivatives (e.g., via Wannier interpolation), and careful separation of band curvature and geometric contributions are standard. Real-space “superfluid-weight markers” can now be constructed to capture local diamagnetic responses even in the presence of strong disorder or spatial inhomogeneity, revealing direct physical signatures such as turbulence in the local supercurrent and modifications to the penetration depth [2505.17349].

## 7. Generalizations, Theoretical Developments, and Open Questions

The functional approach provides a general framework for arbitrary pairing symmetries. In unconventional superconductors, a nontrivial functional term appears that encodes the response of the gap function to gauge fields. This is associated with Wilczek-Zee connections and two-point fidelity tensors, only vanishing for $s$-wave pairing. In realistic topological materials the geometric superfluid weight can thus be governed by non-Abelian quantum geometry and topological invariants. In composite or multi-orbital systems, the lattice geometric (orbital-position dependent) term cancels in the physical observable, ensuring gauge invariance and independence from microscopic model details [2505.09249, 2405.11260].

The true minimal bounds, universality, and the fate of the geometric superfluid weight under strong correlation, disorder, or topological transitions remain active topics of investigation. The generality of the real-space quantum metric approach enables studies in amorphous, fractal, and other non-crystalline media, offering a unified description of geometric effects in superfluidity.

| Contribution    | Dependency (Summary)                  | Dominance Regime                                           |
|-----------------|--------------------------------------|------------------------------------------------------------|
| Conventional    | Band curvature ($\partial_{k}^2 \epsilon$) | Dispersive/wide-band, conventional superconductors          |
| Geometric       | Quantum metric of Bloch states         | Flat/topological bands, quasicrystals, critical states      |
| Functional      | Gap function and Wilczek-Zee connection | Unconventional ($k$-dependent) pairing, topological bands   |

In summary, the superfluid weight tensor synthesizes a comprehensive, symmetry-respecting, and quantum-geometric account of superconducting phase stiffness, providing a bridge between condensed-matter topology, strong correlations, and observable electromagnetic responses. Its refined understanding, particularly the interplay of geometric, topological, and conventional effects, underpins ongoing developments in materials discovery and theoretical superconductivity [1610.01803, 2603.10955, 2110.14663, 2507.20540, 2405.11260, 2501.16965, 2505.17349].

Source: https://www.emergentmind.com/topics/superfluid-weight-tensor