Incommensurate pair-density-wave correlations in two-leg ladder $t$--$J$--$J_\perp$ model

By combining density-matrix renormalization group simulations with bosonization analysis, this study identifies a robust spin-gapped phase in the two-leg $t$-$J$-$J_\perp$ ladder model featuring incommensurate pair-density-wave correlations driven by distinct interlayer and intralayer pairing mechanisms, with potential relevance to bilayer nickelates and optical lattice experiments.

The Gist

Imagine a microscopic world made of two parallel train tracks (called "legs") running side-by-side. On these tracks, tiny particles called electrons are trying to move around. In this specific experiment, the scientists created a special setup where electrons can hop back and forth along their own track, and they can also "talk" to each other across the gap between the tracks, but they cannot jump directly from one track to the other.

The researchers discovered a strange and beautiful new state of matter that happens when the two tracks are not equally crowded.

The Setup: An Imbalanced Crowd

Think of the two tracks as two lanes of a highway.

• Lane 1 is somewhat empty.
• Lane 2 is more crowded.

This imbalance is called "polarization." In the past, scientists mostly studied these systems when both lanes had the exact same number of cars. But here, the authors asked: "What happens if one lane is busier than the other?"

The Discovery: A Wavy Dance

When the lanes are imbalanced, the electrons don't just pair up and move smoothly like a normal superconductor (where electricity flows with zero resistance). Instead, they start doing a complex, wavy dance called a Pair-Density Wave (PDW).

The paper identifies two specific types of this dance happening at the same time:

1. The "Mismatched" Dance (Interlayer PDW):
   Imagine a dancer on the left track trying to hold hands with a dancer on the right track. Because the tracks have different densities of people, the "steps" (momentum) of the dancers don't match up perfectly.

   • The Result: They form pairs, but these pairs are constantly moving forward in a wave pattern. It's like a wave of hand-holding that travels down the tracks. The scientists call this an "incommensurate" wave because the rhythm doesn't fit neatly into the background grid of the tracks. It's driven by the fact that the two lanes are different sizes.

2. The "Echo" Dance (Intralayer PDW):
   Now, look at the dancers on just one track. Even though they are on the same track, they are influenced by the dancers on the other track.

   • The Result: The dancers on the crowded track start to pair up in a rhythm that is actually a "mirror image" or an echo of the rhythm on the empty track. It's as if the empty track is whispering a beat, and the crowded track is dancing to that beat, creating a wave pattern that is distinct from the first type.

Why This Matters (According to the Paper)

The authors found that this "wavy dance" state is very stable and robust. It exists across a wide range of conditions, as long as the two tracks remain imbalanced.

• The "Goldilocks" Zone: If the tracks are perfectly balanced (no imbalance), the dance is smooth and uniform. If one track is completely empty, the dance changes again. But in the middle, where there is a partial imbalance, this special "incommensurate" wave state appears.
• The "Spin Gap": In this state, the "spins" (a quantum property like a tiny internal magnet) of the electrons get locked in place and stop fluctuating wildly. This is a key feature that makes this state unique.

The Catch: A Tiny Leak

The paper also tested what happens if you allow the electrons to jump directly between the tracks (a "leak" or tunneling).

• The Result: Even a tiny bit of jumping between tracks starts to destabilize this special wavy dance. Eventually, if the jumping is strong enough, the dance changes into a different, simpler pattern (called "charge-4e correlations"). However, the paper notes that for very small amounts of jumping, the special wavy dance is surprisingly tough and can survive for a long time before changing.

Real-World Connection

The authors suggest that this model isn't just a math game. It could be built in the real world using optical lattices (traps made of laser light) where scientists can control the number of atoms in each "lane" with lasers.

They also mention a connection to a real material called La3Ni2O7 (a type of nickelate), which is a high-temperature superconductor. The behavior of electrons in this material might be similar to the "wavy dance" described in this paper, especially under high pressure.

Summary

In short, the paper describes a new, stable state of matter where electrons on two parallel tracks form a complex, wavy pairing pattern because the tracks are unevenly crowded. It's a delicate balance between two different types of rhythmic dancing, driven by the difference in crowd size, which creates a unique state that is hard to destroy but fragile if the tracks start mixing too much.

Technical Summary

Technical Summary: Incommensurate Pair-Density-Wave Correlations in Two-Leg Ladder t–J–J⊥ Model

Problem Statement
Pair-density-wave (PDW) superconductivity, characterized by Cooper pairing at finite center-of-mass momentum, is a central theme in strongly correlated systems. While the Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase (driven by external Zeeman fields) and commensurate PDW phases (e.g., in one-dimensional Kondo–Heisenberg models) are well-studied, incommensurate PDW (iC-PDW) correlations remain theoretically less explored and have not been unambiguously demonstrated in numerical studies. Furthermore, the effects of symmetry breaking (specifically layer polarization) in models relevant to bilayer nickelates like La$_3$Ni$_2$O$_7$ have been relatively under-explored compared to the mirror-symmetric regime. This work addresses the search for iC-PDW correlations in a two-leg ladder system with population imbalance.

Methodology
The authors investigate a two-leg ladder $t-J-J_\perp$ model, consisting of two $t-J$ chains coupled solely by an interlayer spin-exchange interaction $J_\perp$, with no interlayer charge hopping ($t_\perp = 0$). The system is partially polarized along the rung direction, creating a population imbalance between the legs.

• Numerical Approach: The study employs Density-Matrix Renormalization Group (DMRG) simulations using the TeNPy package. Calculations were performed on systems with lengths $L=64$ to $128$ and bond dimensions $\chi$ up to 5500. Both finite-DMRG (to analyze real-space correlations and spin gaps) and infinite-DMRG (iDMRG, to extract central charges and correlation lengths via transfer-matrix techniques) were utilized.
• Analytical Approach: Abelian bosonization was used to analytically clarify the nature of the PDW phase. The model was treated by decoupling the legs and then coupling them via $J_\perp$, transforming fields into symmetric (even) and antisymmetric (odd) sectors to identify gapped and gapless modes.

Key Contributions and Results
The paper reports the discovery of a generalized Luther-Emery liquid phase characterized by two distinct types of iC-PDW correlations driven by layer polarization $P$:

1. Interlayer iC-PDW (Dominant):

   • Mechanism: Arises from pairing between layers with mismatched Fermi momenta ($k_{F,1} \neq k_{F,2}$) due to charge conservation per leg.
   • Characteristics: Exhibits FFLO-like oscillations with a pairing momentum $q = \delta k_F = |k_{F,1} - k_{F,2}|$.
   • Evidence: DMRG results show algebraic decay of interlayer pair correlations with the smallest power-law exponent ($\alpha_\perp$), indicating quasi-long-range order. The Fourier transform of these correlations shows a sharp peak at $q = \delta k_F$.

2. Intralayer iC-PDW (Subdominant):

   • Mechanism: Arises from the coupling between charges on one layer and spin fluctuations on the opposite layer. It is understood as a descendant of the commensurate PDW found in the generalized Kondo model (where one leg is half-filled).
   • Characteristics: The pairing on the less populated leg (Leg 2) oscillates with momentum $q = 2k_{F,1}$ (twice the Fermi momentum of the opposite leg).
   • Evidence: These correlations decay faster than the interlayer ones ($\alpha_{\parallel,2} > \alpha_\perp$) but are still algebraic. The peak at $2k_{F,1}$ is sharp at high polarization but broadens as $P$ decreases.

3. Phase Diagram and Robustness:

   • Spin Gap: A finite spin gap exists for $J > J_c \approx 0.78$ across the entire polarization range $0 \le P \le 1$.
   • Central Charge: iDMRG analysis reveals a central charge $c=2$ for the intermediate polarization regime ($0 < P < 1$), indicating two gapless charge modes and gapped spin modes. This contrasts with $c=1$ at $P=0$ (Luther-Emery liquid) and $P=1$ (commensurate PDW).
   • Effect of Interlayer Hopping ($t_\perp$): The introduction of a small $t_\perp$ breaks individual leg particle number conservation. Bosonization suggests this drives an instability toward a $c=1$ phase with charge-4e correlations. However, numerical simulations on accessible system sizes show that the charge-2e PDW order remains robust with power-law decay, and the induced gap is too small to be detected in finite-size simulations.

Significance and Claims
The authors conclude that the two-leg $t-J-J_\perp$ ladder serves as a unified platform for realizing distinct pairing regimes controlled by polarization:

• $P=0$: Uniform Luther-Emery liquid.
• $0 < P < 1$: Interlayer iC-PDW with a finite spin gap.
• $P=1$: Commensurate PDW.

The paper claims that this model provides a clear mechanism for iC-PDW correlations without requiring external magnetic fields, driven instead by intrinsic population imbalance. The authors emphasize the relevance of this model to bilayer nickelates (La$_3$Ni$_2$O$_7$), where strong Hund's coupling generates effective interlayer spin exchange, and suggest that the model is realizable in optical lattice platforms where polarization is easily tunable. The work motivates future experiments targeting the finite-polarization regime to probe polarization-driven PDW transitions. The authors remain modest regarding the stability of this phase in higher dimensions or with finite $t_\perp$ in the thermodynamic limit, noting that while the phase is robust in their simulations, theoretical arguments suggest a flow to a charge-4e phase at infinite correlation lengths.