Kinetic obstruction to pairing in the doped Kitaev-Heisenberg ladder

Using density-matrix renormalization group simulations on a two-leg ladder, this study reveals that hole doping in the Kitaev-Heisenberg model exhibits strong kinetic obstruction to pairing at high hopping strengths while establishing a doping-dependent phase diagram where superconducting, charge-density-wave, and spin-density-wave correlations dominate in distinct regimes defined by the interplay of kinetic energy and magnetic interactions.

The Gist

Imagine a microscopic dance floor made of a honeycomb pattern. On this floor, tiny dancers (electrons) are holding hands in a very specific, rigid way. This is the Kitaev model, a special state of matter known as a "Quantum Spin Liquid." In this liquid state, the dancers are so entangled that they never settle into a fixed pattern, even when it's freezing cold. They are constantly jiggling and swapping partners in a chaotic but harmonious way.

Now, imagine we invite some new dancers onto the floor—holes (missing electrons). These new dancers want to move around. The big question scientists have been asking is: When these new dancers arrive, do they pair up to dance together (superconductivity), or do they just run around and push each other away?

This paper investigates that question using a "two-lane highway" version of the dance floor (a ladder geometry) to see what happens when we add these extra dancers. Here is the story they found, broken down simply:

1. The "Speed Limit" of the Dancers

The most important discovery is that speed matters.

• The Slow Dancers (Low Energy): When the new dancers move slowly, they are like shy people at a party. They stick close together. In this "slow" regime, they find each other and form pairs. This is the "good news" for superconductivity.
• The Fast Dancers (High Energy): When the dancers move too fast (high kinetic energy), they become like rowdy kids running down a hallway. They crash into the existing dance partners, breaking the delicate "Quantum Spin Liquid" connections. Instead of pairing up, they start pushing each other apart or forming rigid lines (magnetic order).

The Analogy: Think of the Quantum Spin Liquid as a delicate glass sculpture.

• If you gently place a marble on it (slow hole), the marble might roll into a groove and pair up with another marble.
• If you throw the marble at it (fast hole), the glass shatters. The structure breaks, and the marbles just scatter.

2. The "Kinetic Obstruction"

The paper calls this phenomenon "Kinetic Obstruction." It means that the motion of the holes actually stops them from pairing up.

In the specific "Antiferromagnetic" version of this dance floor (where neighbors like to face opposite directions), pairing only happens if the dancers move slower than a specific speed limit (about 65% of the maximum speed). If they go faster, the pairing breaks, and the dancers start organizing into a rigid, marching band formation (Spin Density Wave) instead of dancing in pairs.

In the "Ferromagnetic" version (where neighbors like to face the same direction), the speed limit is even stricter. The dancers must move very slowly to pair up. If they speed up even a little, they immediately turn into a "Nagaoka Ferromagnet"—a state where they all align in a rigid, magnetic line, killing any chance of superconductivity.

3. The "Flux" Detective Work

How did they know the dancers were breaking the floor? They looked at a special "floor sensor" called the Plaquette Operator.

• The Metaphor: Imagine the dance floor is made of hexagonal tiles. In the perfect liquid state, every tile is "happy" (value of 1).
• The Observation: When a slow dancer moves, they disturb one or two tiles, but the floor mostly stays happy.
• The Breakdown: When the dancers move too fast, they smash up the tiles. The researchers saw that the "happiness" of the floor dropped significantly.
• The Clue: They also noticed that the pattern of the broken tiles looked exactly like the pattern of the dancers themselves. If the dancers were paired up, the broken tiles formed one big puddle. If the dancers were running apart, the broken tiles split into two separate puddles. This proved that the dancers' speed directly controlled whether they stayed together or ran away.

4. The Map of the Dance Floor

The researchers drew a map (Phase Diagram) showing what happens in different zones:

• The "Rung-Singlet" Zone: This is the sweet spot. Here, the dancers naturally want to pair up across the rungs of the ladder. When you add slow dancers here, you get Superconductivity (they flow without resistance).
• The "Stripy" Zone: Here, the dancers form rigid lines (Charge Density Waves). They aren't pairing up; they are just lining up.
• The "Kitaev" Zones (Pure Spin Liquid): This is where the magic happens, but it's fragile. If the dancers are too fast, the liquid turns into a solid magnetic block.

Why Does This Matter?

This research helps us understand why some materials might become superconductors (conduct electricity with zero loss) and others won't.

• Real-world implication: Many real materials (like certain iridates or ruthenates) are close to this "Kitaev" state. The paper suggests that if we try to dope them (add holes) in a way that makes the holes move too fast, we will never get superconductivity. The kinetic energy will destroy the pairing.
• The Solution: To get superconductivity in these exotic materials, we might need to slow the holes down. This could be done by creating special "sandwich" structures (heterostructures) or using quantum simulators to control the speed of the particles.

Summary

In short: Motion can kill magic.
In these exotic quantum materials, the ability of electrons to pair up and superconduct depends entirely on how fast they are moving. If they move too fast, they break the delicate quantum connections and refuse to pair up. To get the "super" in superconductivity, you have to keep the dancers slow and steady.

Technical Summary

1. Problem Statement

The paper addresses the fundamental question of how quantum spin liquids (QSLs), specifically the Kitaev spin liquid, respond to charge doping. While the undoped Kitaev model on a honeycomb lattice is exactly solvable and hosts a gapless QSL ground state, the introduction of mobile holes (doping) creates a complex interplay between kinetic energy, spin, charge, and $Z_2$ flux degrees of freedom.

Key open questions include:

• Does doping a Kitaev QSL lead to exotic superconductivity, as hypothesized in some theoretical works?
• How does the kinetic energy of the doped holes (hopping strength $t$) influence the stability of the QSL background and the formation of hole pairs?
• What are the resulting phases in the doped Kitaev-Heisenberg ($t-J-K$) model, and how do they differ between the antiferromagnetic (AFK) and ferromagnetic (FK) Kitaev limits?

Previous studies on wider cylinders (3- and 4-leg) suggested that superconducting correlations decay exponentially in the "fast-hole" regime (large $t$), but a comprehensive phase diagram and the microscopic mechanism of "kinetic obstruction" to pairing remained to be fully elucidated.

2. Methodology

The authors employ the Density Matrix Renormalization Group (DMRG) method to study the doped Kitaev-Heisenberg model on a two-leg ladder geometry. This geometry is chosen for its numerical tractability while maintaining strong correspondence to the 2D honeycomb limit.

• Model Hamiltonian: The study utilizes the $t-J-K$ Hamiltonian:
  $$H = K \sum_{\langle i,j \rangle} S_i^\gamma S_j^\gamma + J \sum_{\langle i,j \rangle} \mathbf{S}_i \cdot \mathbf{S}_j - t \sum_{\langle i,j \rangle, \sigma} P (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.}) P$$
  where $K = \sin \phi$ (Kitaev interaction), $J = \cos \phi$ (Heisenberg exchange), and $t$ is the hole hopping amplitude. The projector $P$ enforces no double occupancy.
• System Parameters:
  • System sizes up to $N = 126$ sites ($2 \times 63$).
  • Open boundary conditions (OBC) along the leg direction and periodic boundary conditions (PBC) along the rung direction.
  • Doping levels ($n_h$) ranging from 1 to 16 holes.
• Observables: The authors calculate:
  • Binding Energy ($\Delta E$): To determine if holes form bound pairs ($\Delta E < 0$).
  • Correlation Functions: Spin-spin ($S(r)$), charge-charge ($C(r)$), and singlet/triplet pair-pair correlations ($P_S(r), P_T(r)$).
  • Plaquette Operator ($W_p$): A conserved quantity in the pure Kitaev limit ($J=t=0$) used to track the integrity of the $Z_2$ flux background.
  • Static Structure Factors: To identify magnetic ordering (SDW, FM, etc.).

3. Key Contributions

1. Identification of Kinetic Obstruction: The paper establishes that the ability of holes to pair in a Kitaev QSL is strictly controlled by their kinetic energy. Pairing is only possible in the "slow-hole" regime ($t \lesssim 0.65|K|$ for AFK and $t \lesssim 0.1|K|$ for FK).
2. Microscopic Mechanism via Plaquette Operator: The authors demonstrate a direct link between the spatial profile of the plaquette operator (flux) and charge density. They show that the "obstruction" to pairing is marked by a distinct change in the slope of the average plaquette value and a transition from a single-minimum to a two-minimum spatial profile, indicating the breakdown of the QSL background.
3. Comprehensive Doped Phase Diagram: A detailed phase diagram is constructed for the doped system at intermediate hopping ($t=1$), identifying distinct phases including superconducting (SC), charge-density-wave (CDW), and various spin-density-wave (SDW) orders.
4. Distinction between AFK and FK Regimes: The study reveals that the FK limit is significantly more fragile to doping than the AFK limit, with pairing vanishing at much lower hopping strengths due to the rapid emergence of Nagaoka-type ferromagnetism.

4. Key Results

A. Kinetic Obstruction to Pairing

• AFK Limit ($K>0, J=0$): Hole pairing (negative binding energy) persists only for $t \lesssim 0.65|K|$. Beyond this critical hopping ($t_c^{AFK}$), the binding energy becomes positive, and holes repel.
• FK Limit ($K<0, J=0$): Pairing is extremely fragile, surviving only for $t \lesssim 0.1|K|$. At slightly higher $t$, the system transitions into a Nagaoka-type ferromagnetic state (SDW1).
• Mechanism: In the slow-hole limit, holes act as quasi-static vacancies that minimally disturb the flux background, allowing them to bind to minimize broken magnetic bonds. As $t$ increases, holes become mobile, disrupting the $Z_2$ flux structure (measured by the plaquette operator $W_p$). This "kinetic obstruction" destroys the QSL background necessary for pairing.
• Evidence: The spatial profile of the plaquette operator $\langle W_p(l) \rangle$ mimics the charge density. For bound pairs ($t < t_c$), a single broad minimum is observed. For unbound holes ($t > t_c$), two distinct minima appear, signaling independent hole motion.

B. Phase Diagram at Intermediate Hopping ($t=1$)

At $t=1$ (where pairing is generally obstructed in the pure Kitaev limits), the phase diagram as a function of interaction angle $\phi$ and doping $n_h$ reveals:

• Rung-Singlet (RS) Phase: Upon doping, this phase develops dominant singlet superconducting correlations (algebraic decay of pair-pair correlations). This is the only region where superconductivity is robust.
• Stripy Phase (ST): Evolves into an incommensurate spin-density-wave order (SDW2).
• AFK Limit: At weak doping, a disordered phase (DS) with power-law decay of both spin and charge correlations is found. At higher doping, it transitions to an incommensurate SDW3 phase.
• FK Limit: Evolves immediately into a commensurate SDW1 (Nagaoka ferromagnetism) even at low doping.
• Transition Region (ST-RS): Near the transition between the stripy and rung-singlet phases, a narrow region dominated by charge-density-wave (CDW) correlations is identified at weak doping.

C. Correlation Analysis

• Superconductivity: Strongest in the rung-singlet phase. Pair-pair correlations decay algebraically, while spin-spin correlations decay exponentially.
• Magnetism: In the AFK and FK limits, and the stripy phase, spin-spin correlations dominate or become quasi-long-ranged, suppressing superconductivity.
• Charge Order: CDW correlations dominate near the ST-RS boundary and in the AFK regime at intermediate doping.

5. Significance and Implications

• Resolution of the "Fast vs. Slow Hole" Debate: The paper provides robust numerical evidence that the physics of doped Kitaev systems is fundamentally different for slow and fast holes. Fast holes ($t \gtrsim |K|$) destroy the QSL background, preventing pairing, whereas slow holes allow for bound states.
• Material Relevance:
  • Realistic Kitaev candidate materials (e.g., $\alpha$-RuCl$_3$, Na$_2$IrO$_3$) typically have significant Heisenberg exchange and potentially large effective hopping, placing them in the "fast hole" regime where superconductivity is unlikely.
  • The authors suggest that superconductivity in Kitaev systems might only be accessible in the slow-hole regime, which could be engineered via heterostructures or quantum simulation platforms.
  • Materials closer to the Heisenberg limit (e.g., YbCl$_3$) are more likely to exhibit superconductivity upon doping.
• Methodological Insight: The study highlights the utility of the plaquette operator as a diagnostic tool for the stability of QSLs under perturbation, showing that its spatial profile directly reflects the localization/delocalization of doped carriers.

In conclusion, the paper demonstrates that kinetic energy acts as a primary obstruction to pairing in doped Kitaev spin liquids. Superconductivity is not a generic outcome of doping these systems but is highly sensitive to the ratio of hopping to interaction strength, surviving only in specific parameter regimes (Rung-Singlet) and low hopping limits.