Finite-Momentum Cooper Pairing in Superconductors

• Finite-momentum Cooper pairing is a superconducting state where Cooper pairs carry nonzero momentum, resulting in a spatially modulated order parameter.
• It emerges when broken symmetries from Zeeman fields, spin-orbit coupling, or anisotropic interactions force a shift in the conventional zero-momentum pairing paradigm.
• Experimental techniques such as Josephson junction studies, STM/QPI imaging, and 2RDM analysis provide key evidence for exotic FFLO, pair-density wave, and helical superconducting phases.

Finite-momentum Cooper pairing is the unconventional superconducting phenomenon where the center-of-mass momentum of the Cooper pair is nonzero, resulting in spatially modulated superconducting states. This concept generalizes the Bardeen-Cooper-Schrieffer (BCS) paradigm and underlies a broad spectrum of exotic phases, including Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states, pair-density waves (PDWs), and emergent superconducting orders in spin-orbit-coupled, noncentrosymmetric, and topologically nontrivial systems. The recent advances in theory and experiment have enabled detailed studies of its microscopic origins, manifestations in strongly correlated and engineered heterostructures, and the connection to odd-frequency pairing and topological superconductivity.

1. Theoretical Foundations of Finite-Momentum Cooper Pairing

In conventional BCS theory, the superconducting condensate wavefunction is spatially uniform, with Cooper pairs formed by electrons of opposite momentum $(\mathbf{k}, -\mathbf{k})$ and zero net center-of-mass momentum. When either an external perturbation or intrinsic many-body effects break this symmetry, a pairing state with finite momentum $\mathbf{Q}$ emerges: the condensate attains a spatial modulation $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$, leading to an order parameter characterized by a nonzero wavevector.

The minimal microscopic description is provided by the Bogoliubov–de Gennes (BdG) Hamiltonian with shifted momenta: $$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$ where $\Psi_{\mathbf{k}}$ is the Nambu spinor, and $\xi_{\mathbf{k}}$ is the single-particle dispersion. Diagonalization yields Bogoliubov quasiparticles with energy spectra shifted by $\mathbf{Q}$, and the gap equation must be solved self-consistently for the spatially modulated pairing amplitude $\Delta_{\mathbf{k}}(\mathbf{Q})$ (Zhu et al., 2020, Setty et al., 2022).

The canonical realization of finite-momentum pairing is the FFLO state, predicted to appear in spin-imbalanced superconductors under Zeeman fields, where the Fermi surfaces of opposite spins are mismatched and favor pairing at finite total momentum (Chakraborty et al., 2022, Qu et al., 2013). In the presence of crystalline anisotropy or strong correlations, more intricate spatial patterns can arise, including stripe-like LO states and PDWs.

2. Microscopic Mechanisms and Phase Instability

The transition from a uniform ($\mathbf{Q}=0$) to a finite-momentum pairing state arises from various microscopic mechanisms:

• Zeeman-induced Fermi surface mismatch: In magnetic fields, spin-split Fermi surfaces prevent $k$-space nesting at $\mathbf{Q}$0 and $\mathbf{Q}$1, favoring FFLO states with $\mathbf{Q}$2 (Chakraborty et al., 2022).
• Spin-orbit coupling and inversion symmetry breaking: Rashba or "Ising"-type SOCs induce momentum-dependent spin splitting; in noncentrosymmetric materials with strong SOC, functional forms of the order parameter can directly favor states with nonzero $\mathbf{Q}$3 even in the absence of applied fields (Asaba et al., 2024, Yang et al., 2024, Steinbok et al., 2016).
• Magnetoelectric or Lifshitz invariants: In systems with local inversion symmetry breaking, free energy gradients such as $\mathbf{Q}$4 produce an energetic minimum at $\mathbf{Q}$5, with $\mathbf{Q}$6, giving rise to orbital-driven finite-momentum pairing (Yang et al., 2024).
• Anisotropic and multi-band interactions: For anisotropic dispersions and $\mathbf{Q}$7-wave or interband interactions, the two-particle scattering vertex $\mathbf{Q}$8 can favor pair formation at a finite $\mathbf{Q}$9 once interaction strength exceeds a critical value, as shown for $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$0-wave separable interactions (Setty et al., 2022).

The linearized stability analysis involves evaluating the mean-field pairing susceptibility (pair bubble) $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$1. In gapped band insulators or multiband systems, nontrivial behavior in the vertex function and band structure gives rise to a linear-in- $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$2 term in $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$3, energetically favoring a PDW at $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$4 (Nikolic et al., 2010). Renormalization group analyses reveal that in 2D, fluctuations can be enhanced, leading to transitions into bosonic Cooper-pair insulating phases separating uniform and finite-momentum superconductivity.

Odd-frequency pairing components further destabilize the uniform phase when the Meissner kernel becomes negative, guaranteeing that the first instability at $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$5 occurs at nonzero $$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$6 whenever the odd-frequency contribution dominates (Sato et al., 3 Mar 2025).

3. Experimental Realizations and Diagnostic Signatures

Recent advances have enabled the experimental identification and characterization of finite-momentum pairing in proximity-coupled heterostructures, intrinsic quantum materials, and engineered devices:

• Josephson Junctions and Proximitized Systems: In systems such as Bi$$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$7Se$$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$8 and HgTe quantum wells, the application of in-plane magnetic field or tuning of spin-orbit coupling yields finite-momentum pairing, evidenced by spatially modulated Josephson interference, trident-like Fraunhofer patterns, and supercurrent nodes at discrete field values. The total pair momentum is quantitatively extracted from features such as node positions and side-branch slopes, in agreement with predictions based on Zeeman and orbital contributions (Chen et al., 2018, Hart et al., 2015, Pal1 et al., 2021).
• Quasiparticle Interference (QPI) Imaging: STM-based QPI, as demonstrated in Bi$$\Delta(\mathbf{r}) = \Delta_0\,e^{i\mathbf{Q}\cdot\mathbf{r}}$$9Te$$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$0/NbSe$$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$1 films, can directly resolve the Fermi surface segmentation imposed by Cooper pair momentum. The transition from a gapped to “segmented” (arcs) Fermi surface, with gapless Bogoliubov quasiparticles restricted to angular intervals set by Doppler shift compensation, provides a real-space and momentum-resolved diagnostic for finite-momentum pairing (Zhu et al., 2020).
• Non-reciprocal Transport and the Superconducting "Diode" Effect: In tricolor superlattices of CeCoIn$$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$2, the supercurrent diode effect provides a sensitive probe of finite-momentum (helical) pairing. A second-harmonic voltage response ($$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$3) tracks the Cooper-pair momentum and reveals phase transitions into helical regimes, as does direct measurement of diode efficiency in topological semimetal junctions (Asaba et al., 2024, Pal1 et al., 2021).
• Rotational and Symmetry Breaking in Bulk Heterostructures: In 3D Ising superconductors such as 4Hb-TaS$$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$4, transport anisotropy and rotational symmetry breaking (from $$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$5 to $$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$6) under field reflect the emergence and orientation locking of finite-$$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$7 pairing induced by orbital-magnetoelectric coupling (Yang et al., 2024).

These experimental platforms provide both direct and indirect signatures of the underlying spatial modulation, with critical field, symmetry, and device-dependent observables used to extract the key momentum-dependent parameters.

4. Many-Body and Correlated Electron Perspectives

In the context of interacting electronic systems, particularly strongly correlated materials and lattice models, finite-momentum pairing is most rigorously diagnosed and characterized by the two-particle reduced density matrix (2RDM) formalism under the Penrose-Onsager criterion. The eigenvalues and eigenvectors of the 2RDM projected onto center-of-mass momentum sectors yield the condensate fraction and internal wavefunction structure, enabling unbiased detection and visualization of finite-momentum (FFLO-type) order. Large-scale AFQMC and DMRG studies in the 2D attractive Hubbard model under Zeeman field reveal finite-size scaling to a dominant nonzero-$$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$8 condensate fraction, spatially oscillating pair order, and momentum-space peaks at $$H_{\text{BdG}} = \sum_{\mathbf{k}} \Psi^{\dagger}_{\mathbf{k}} \begin{bmatrix} \xi_{\mathbf{k} + \mathbf{Q}/2} & \Delta_{\mathbf{k}} \ \Delta^*_{\mathbf{k}} & -\xi_{-\mathbf{k} + \mathbf{Q}/2} \ \end{bmatrix} \Psi_{\mathbf{k}},$$9 (Karlsson et al., 26 Jan 2026).

In bosonic or pseudogap systems, the transition into PDW or Cooper-pair insulating phases can be mapped via fluctuations in the ground state, particle-number scaling, and the collapse of the gap across the normal-to-PDW transition. Quantum fluctuation effects in 2D promote non-BCS, fluctuation-dominated, and even bosonic superfluid-insulator character (Nikolic et al., 2010).

5. Topological, Orbital, and Exotic Finite-Momentum States

The interplay between finite-momentum Cooper pairing and topology opens direct pathways to topological superconductors and the realization of Majorana fermion edge modes. In FFLO states with strong spin-orbit coupling and Zeeman fields, as in engineered cold-atom and solid-state heterostructures, the self-consistent solution yields regions of topological FFLO superfluids characterized by nontrivial invariants and chiral Majorana bound states (zero modes) at the edges or defects (Qu et al., 2013). Precise criteria for topologically nontrivial gapped FFLO require that the generalized gap-closing condition is satisfied while maintaining finite momentum and the nonvanishing order parameter.

In materials where the magnetization vanishes but prominent anisotropic spin splitting exists (altermagnets), spatially modulated superconductivity is possible via orbital effects alone. Here, finite-momentum pairing exhibits strong orientation dependence, leading to field-controllable $\Psi_{\mathbf{k}}$0-$\Psi_{\mathbf{k}}$1 Josephson transitions and new classes of interference behavior without net magnetization (Zhang et al., 2023).

Fractional Chern insulators and quantum Hall systems provide a fundamentally different realization: the magnetic translation symmetry of the lattice and flux attachment physics enforce the fragmentation of the composite-fermion Fermi surface, such that pairing of these composite fermions naturally occurs at a set of nonzero commensurate finite-momenta. The resulting phases intertwine PDW/stripe order and topological order, supporting non-Abelian statistics in the presence of Majorana zero modes (Sohal et al., 2020).

6. Spectroscopic Probes and Diagnostic Tools

Probing the momentum structure of the superconducting condensate and distinguishing finite-momentum states demands experimental techniques with both momentum and spin resolution:

• Two-electron ARPES: By measuring the joint emission spectrum, two-electron ARPES can directly probe the Cooper pair momentum by mapping the sum momentum and spin of simultaneously emitted electron pairs. The appearance of a zero-energy peak at $\Psi_{\mathbf{k}}$2 differentiates between uniform, PDW, or FFLO order, and by scanning this sum-momentum, the full structure of $\Psi_{\mathbf{k}}$3 can be reconstructed (Mahmood et al., 2021).
• STM/QPI and ARPES: Momentum-resolved STM techniques, including QPI imaging, have demonstrated the ability to observe Doppler-shifted Fermi arcs (“segmented” Fermi surface), directly revealing finite-momentum pairing from the evolution of quasiparticle interference with field or current (Zhu et al., 2020).

Table: Key Theoretical and Diagnostic Aspects

Mechanism / Observable	Theoretical Origin / Formula	Diagnostic Tool / Signature
Zeeman-driven FFLO	$\Psi_{\mathbf{k}}$4	Josephson interference, 2$\Psi_{\mathbf{k}}$5-ARPES
SOC-induced finite $\Psi_{\mathbf{k}}$6	Rashba/SIA splitting, Lifshitz invariants	Nonreciprocal $\Psi_{\mathbf{k}}$7, symmetry breaking
QPI "segmented" Fermi surface	$\Psi_{\mathbf{k}}$8	STM QPI, ARPES, density of states
2RDM Penrose-Onsager	Leading 2RDM eigenvalue at $\Psi_{\mathbf{k}}$9	DMRG, QMC, real/momentum-pair functions
Odd-frequency-driven instability	$\xi_{\mathbf{k}}$0: paramagnetic Meissner kernel	$\xi_{\mathbf{k}}$1-selected $\xi_{\mathbf{k}}$2, magnetotransport
Topological FFLO and Majorana	Gap closing at $\xi_{\mathbf{k}}$3: topological invariant $\xi_{\mathbf{k}}$4	Edge spectroscopy, transport

7. Broader Implications and Outlook

Finite-momentum Cooper pairing exposes the rich interplay between broken symmetries, spin-orbit coupling, band structure, interactions, topology, and dimensionality. It underpins nontrivial ground states, as in pair-density-waves, helical superconductors, and topological FFLO phases supporting Majorana modes. These phenomena are now being visualized directly and controlled in systems ranging from proximitized quantum wells and topological materials to tricolor superlattice heterostructures and engineered Dirac semimetals.

Their realization has prompted proposals for momentum-tunable superconducting diodes, nonreciprocal devices, and quantum computation architectures (Pal1 et al., 2021, Asaba et al., 2024). The cross-fertilization between momentum-resolved spectroscopies, advanced simulation (2RDM-based frameworks), and materials engineering positions finite-momentum pairing at the convergence of modern condensed matter theory and experiment, with ongoing efforts to elucidate its direct relevance to correlated unconventional superconductors, intertwined orders in Chern insulators, and emergent phases protected by crystalline or topological symmetries.

References: (Zhu et al., 2020, Karlsson et al., 26 Jan 2026, Setty et al., 2022, Yang et al., 2024, Asaba et al., 2024, Pal1 et al., 2021, Zhang et al., 2023, Hart et al., 2015, Chen et al., 2018, Chakraborty et al., 2022, Mahmood et al., 2021, Nikolic et al., 2010, Sato et al., 3 Mar 2025, Qu et al., 2013, Steinbok et al., 2016, Sohal et al., 2020).