Pair density waves (PDWs) are superconducting states formed by Cooper pairs of electrons containing a nonzero center-of-mass momentum. They are characterized by a spatially modulated order parameter and may occur in a variety of emerging quantum materials such as cuprates, transition-metal dichalcogenides (TMDs), and Kagome metals. Despite extensive theoretical and numerical studies seeking PDWs in a variety of lattices and interacting settings, there is currently no exact mechanism that spontaneously favors a modulated solution of the superconducting order parameter. Here, we study the problem of two electrons subject to an anisotropic attractive potential. We solve the two-body Schr\"odinger wave equation exactly to determine the pair binding energy as a function of the center-of-mass momentum. We find that a modulated (finite momentum) pair is favored over a homogeneous (zero momentum) solution above a critical, intermediate interaction strength. Hence our exact result justifies previous mean-field approximations that obtain modulated ground states at finite but large interactions. Using this insight from the exact two-body solution, we construct a variational many-body wave function and show that the conclusions of the two-body problem are robust in the many-body limit. Our results thus lay the theoretical and microscopic foundation for the existence of PDWs.
Pair density waves (PDWs) are superconducting states formed by Cooper pairs of electrons containing a nonzero center-of-mass momentum. They are characterized by a spatially modulated order parameter and may occur in a variety of emerging quantum materials such as cuprates, transition-metal dichalcogenides (TMDs), and Kagome metals. Despite extensive theoretical and numerical studies seeking PDWs in a variety of lattices and interacting settings, there is currently no exact mechanism that spontaneously favors a modulated solution of the superconducting order parameter. Here, we study the problem of two electrons subject to an anisotropic attractive potential. We solve the two-body Schr"odinger wave equation exactly to determine the pair binding energy as a function of the center-of-mass momentum. We find that a modulated (finite momentum) pair is favored over a homogeneous (zero momentum) solution above a critical, intermediate interaction strength. Hence our exact result justifies previous mean-field approximations that obtain modulated ground states at finite but large interactions. Using this insight from the exact two-body solution, we construct a variational many-body wave function and show that the conclusions of the two-body problem are robust in the many-body limit. Our results thus lay the theoretical and microscopic foundation for the existence of PDWs.
A robust theory of the mechanism of pair density wave (PDW) superconductivity (i.e., where Cooper pairs have nonzero center of mass momentum) remains elusive. Here we explore the triangular lattice $t\text{\ensuremath{-}}J\text{\ensuremath{-}}V$ model, a low-energy effective theory derived from the strong-coupling limit of the Holstein-Hubbard model, by large-scale variational Monte Carlo simulations. When the electron density is sufficiently low, the favored ground state is an s-wave PDW, consistent with results obtained from previous studies in this limit. Additionally, a PDW ground state with nematic d-wave pairing emerges in the intermediate range of electron densities and phonon frequencies. For these s-wave and d-wave PDWs arising in states with spontaneous breaking of time-reversal and inversion symmetries, PDW formation derives from valley-polarization and intrapocket pairing.
Abstract Flat-band systems provide an ideal platform for studying strong electron correlations, in which multiple quasi-degenerate charge orders frequently coexist. Motivated by this context, we investigate the pair-density-wave (PDW) superconducting state in a quasi-one-dimensional Creutz lattice featuring two strictly flat bands. Using a renormalized mean-field approach for the strongly correlated t - J model, we obtain self-consistent solutions that demonstrate a PDW phase—characterized by finite-center-of-mass-momentum Cooper pairing—coexisting with charge density wave (CDW) and bond density wave (BDW) orders. Physically, this coexistence phenomenon originates from the perfect band nesting inherent to the flat-band system, which allows for changes in the particle momentum distributions without energy cost. Here, a pure PDW solution without CDW and BDW orders emerges in the underdoped regime when magnetic effects are neglected, while strong Hubbard correlations favor a pure CDW phase. We conclude that a robust PDW phase, accompanied by CDW and BDW orders, occupies the majority of the parameter space. Furthermore, both the PDW order parameter and its critical temperature scale linearly with the superexchange coupling strength—a hallmark of flat-band superconductivity. Overall, our results reveal that a PDW superconducting state can be realized on the Creutz lattice with strong correlations.
A pair density wave (PDW) is a superconductor whose order parameter is a periodic function of space, without an accompanying spatially uniform component. Since PDWs are not the outcome of a weak-coupling instability of a Fermi liquid, a generic pairing mechanism for PDW order has remained elusive. We describe and solve models having robust PDW phases. To access the intermediate coupling limit, we invoke large-N limits of Fermi liquids with repulsive BCS interactions that admit saddle point solutions. We show that the requirements for long-range PDW order are that the repulsive BCS couplings must be nonmonotonic in space and that their strength must exceed a threshold value. We obtain a phase diagram with both finite temperature transitions to PDW order and a T=0 quantum critical point, where non-Fermi liquid behavior occurs.
Pair density waves (PDW) are novel forms of superconducting states that exhibit periodically modulated pairing. A remaining challenge is to elucidate how intrinsic PDW order can emerge robustly in strongly correlated electrons. Here we propose that PDW is prone to form in strongly coupled multiband superconductors simply with interband Cooper pairing between electrons from oppositely dispersing bands. This scenario is heuristically motivated by the observation that uniform interband pairing in such systems would exhibit negative superfluid weight -- a signature of an instability towards pairing modulation, implying that PDW emerges naturally in the true ground state. Using large-scale density-matrix-renormalization-group calculations with finite-size scaling analysis, we demonstrate this PDW mechanism in a minimal model with strong interband attractions. Our simulations reveal power-law superconducting correlations characterized by incommensurate modulations. The exponent Ksc of the power-law PDW correlation decreases systematically with increasing ladder width, confirming a genuine long-range PDW order in the 2D limit. Our study therefore demonstrates a promising route to robust PDW states in multiband systems.
An indispensable ingredient for pair density wave (PDW) superconductivity is the presence of an attractive pairing interaction at finite momentum. Here, we show how this condition can be met with straightforward electron-density interactions in multiband systems. The electron-density interaction, when projected to the band basis, acquires form factors with nontrivial momentum dependence and thereby exhibits a potential tendency to a finite-momentum pairing instability. By applying a mean field analysis to two simple multiband models, the checkerboard lattice and three-band Hubbard model, we find that PDW order can indeed become the leading instability if a strong nearest-neighbor attraction is present. Moreover, the condition for the transition from a uniform superconductor to a PDW superconductor is shown via a simple quantum geometric argument.
The following thesis is divided in two main parts. Chapters 2, 3 and 4 are devoted to the study of the so called pair-density-wave (PDW) superconducting state and some of its connections to electronic liquid crystal (ELC) phases, its topological aspects in a one dimensional model and its appearance in a quasi-one dimensional system. On the other hand, chapter 5 is focused on the investigation of the classical statistical mechanics properties of dimers, in particular, the dimer model on the Aztec diamond graph and its relation with the octahedron equation. In chapter 2 we present a theory of superconducting states where the Cooper pairs have a nonzero center-of-mass momentum, inhomogeneous superconducting states known as a pair-density-waves (PDWs) states. We show that in a system of spin-1/2 fermions in two dimensions in an electronic nematic spin-triplet phase where rotational symmetry is broken in both real and spin space PDW phases arise naturally in a theory that can be analysed using controlled approximations. We show that several superfluid phases that may arise in this phase can be treated within a controlled BCS mean field theory, with the strength of the spin-triplet nematic order parameter playing the role of the small parameter of this theory. We find that in a spin-triplet nematic phase, in addition to a triplet p-wave and spin-singlet d-wave (or s depending on the nematic phase) uniform superconducting states, it is also possible to have a d-wave (or s) PDW superconductor. The PDW phases found here can be either unidirectional, bidirectional, or tridirectional depending on the spin-triplet nematic phase and which superconducting channel is dominant. In addition, a triple-helix state is found in a particular channel. We show that these PDW phases are present in the weak-coupling limit, in contrast to the usual Fulde-Ferrell-Larkin-Ovchinnikov phases, which require strong coupling physics in addition to a large magnetic field (and often both). In chapter 3 we show that the pair-density-wave (PDW) superconducting state emergent in extended Heisenberg-Hubbard models in two-leg ladders is topological in the presence of an Ising spin symmetry and supports a Majorana zero mode (MZM) at an open boundary and at a junction with a uniform d-wave one-dimensional superconductor. Similarly to a conventional finite-momentum paired state, the order parameter of the PDW state is a charge-2e field with finite momentum. However, the order parameter here is a quartic electron operator and conventional mean-field theory cannot be applied to study this state. We use bosonization to show that the 1D PDW state has a MZM at a boundary. This superconducting state is an exotic topological phase supporting Majorana fermions with finite-momentum pairing fields and charge-4e superconductivity. In chapter 4 we provide a quasi-one-dimensional model which can support a PDW state. The model consists of an array of strongly-interacting one-dimensional systems, where the one-dimensional systems are coupled to each other by local interactions.Within the interchain mean-field theory (MFT), we find several SC states from the model, including a conventional uniform SC state, PDW SC state, and a coexisting phase of the uniform SC and PDW states. In this quasi-1D regime we can treat the strong correlation physics essentially exactly using bosonization methods and the crossover to the 2D system by means of interchain MFT. The resulting critical temperatures of the SC phases generically exhibit a power-law scaling with the coupling constants of the array, instead of the essential singularity found in weak-coupling BCS-type theories. Electronic excitations with an open Fermi surface, which emerge from the electronic Luttinger liquid systems below their crossover temperature to the Fermi liquid, are then coupled to the SC order parameters via the proximity effect. From the Fermi surface thus coupled to the SC order parameters, we calculate the quasiparticle spectrum in detail. We show that the quasiparticle spectrum can be fully gapped or nodal in the uniform SC phase and also in the coexisting phase of the uniform SC and PDW parameters. In the pure PDW state, the excitation spectrum has a reconstructed Fermi surface in the form of Fermi pockets of Bogoliubov quasiparticles. In chapter 5 we study the octahedron relation (also known as the A∞ T-system), obeyed in particular by the partition function for dimer coverings of the Aztec Diamond graph. For a suitable class of doubly periodic initial conditions, we find exact solutions with a particularly simple factorized form. For these, we show that the density function that measures the average dimer occupation of a face of the Aztec graph, obeys a system of linear recursion relations with periodic coefficients. This allows us to explore the thermodynamic limit of the corresponding dimer models and to derive exact “arctic” curves separating the various phases of the system.