Gap structure and phase diagram of twisted bilayer cuprates from a microscopic perspective

This paper employs a tight-binding lattice model to investigate the gap structure and phase diagram of twisted bilayer cuprates, revealing that the time-reversal symmetry breaking $d+id'$ state is correlated with the Van Hove singularity's position—which varies with doping and interlayer tunneling—and discussing how these findings reconcile conflicting experimental results.

The Gist

Imagine you have two sheets of superconducting material (like a special kind of metal that conducts electricity with zero resistance). Now, imagine stacking one sheet directly on top of the other, but you twist the top sheet slightly, like turning a page in a book.

This simple act of twisting creates a complex, repeating pattern called a Moiré pattern (think of the wavy lines you see when you hold two window screens slightly out of alignment). Scientists have been trying to figure out what happens to the electricity flowing through these twisted sheets, specifically at a "magic" twist angle of 45 degrees.

Here is the story of what this paper discovered, explained simply:

1. The Great Mystery: The Conflicting Experiments

For a while, scientists had a big disagreement.

• Team A twisted two big crystals and found that at 45 degrees, the electricity flow (called "critical current") almost stopped. They thought this meant a weird, exotic state of matter had formed that broke the rules of time symmetry (like a clock that runs backward).
• Team B took tiny flakes of the same material, twisted them, and found that the electricity kept flowing just fine, even at 45 degrees. They saw no sign of that exotic state.

Why the difference? Was one team wrong? Or was the material behaving differently?

2. The Simulation: A Digital Playground

The authors of this paper didn't build a new lab experiment. Instead, they built a massive digital simulation on a computer. They created a virtual model of these twisted layers, down to the individual atoms.

Think of their model as a giant, complex dance floor.

• The Dancers: Electrons.
• The Music: The twist angle, how many electrons are on the floor (doping), and how strongly the two layers talk to each other (tunneling).
• The Goal: To see what kind of "dance routine" (superconducting state) the electrons naturally fall into.

3. The Discovery: It's All About the "Tunnel"

The paper found that the answer isn't just about the angle (45 degrees or not). It's about how strongly the two layers are connected.

Imagine the two layers are like two rooms separated by a door.

• Weak Connection (Small Door): If the door is barely open (weak tunneling), the electrons in the top room and bottom room don't mix much. In this case, the "Team A" result happens: at 45 degrees, the dance stops, and the current vanishes. This is the exotic, time-breaking state.
• Strong Connection (Big Open Door): If the door is wide open (strong tunneling), the electrons mix freely. In this case, the "Team B" result happens: the dance continues, and the current keeps flowing, even at 45 degrees.

The Analogy:
Think of the electrons as a crowd trying to cross a bridge.

• If the bridge is narrow and shaky (weak tunneling), the crowd gets stuck at a specific angle, and traffic stops.
• If the bridge is wide and sturdy (strong tunneling), the crowd flows smoothly, regardless of the angle.

4. The "Van Hove" Sweet Spot

The paper also discovered a hidden factor: the Van Hove Singularity.
Imagine the electrons are like water in a landscape. Sometimes, the landscape has a "dip" or a "valley" where the water naturally pools. The authors found that the exotic, time-breaking state only appears when the "water level" (doping) and the "shape of the valley" (tunneling) line up perfectly so that the electrons are sitting right in that sweet spot. If you miss that spot, the exotic state disappears.

5. Why the Experiments Disagreed

The paper suggests that the two experimental teams likely had different "door sizes" (tunneling strengths) without realizing it.

• Team A (Crystals): Their surfaces were rough, or the pressure they applied created a very specific, weak connection. This kept them in the "weak tunneling" zone, where the current stops.
• Team B (Flakes): Their thin flakes might have had rougher surfaces or different impurities that accidentally created a stronger connection (a bigger door). This pushed them into the "strong tunneling" zone, where the current keeps flowing.

The Bottom Line

This paper acts like a detective solving a mystery. It tells us that the "exotic" state isn't a one-size-fits-all phenomenon. It is a delicate balance.

• Twist angle matters.
• How many electrons are there matters.
• How strongly the layers touch matters most.

If you tweak any of these, you can switch the material from a state where electricity stops, to a state where it flows freely. This explains why different labs saw different things: they were likely looking at slightly different versions of the same material, where the "tunnel" between the layers was just a little bit wider or narrower.

In short: The universe of twisted superconductors is a chameleon. It changes its colors (behavior) based on how tightly you squeeze the layers together and how you twist them. This paper gives us the map to predict exactly what color it will be.

Technical Summary

Here is a detailed technical summary of the paper "Gap structure and phase diagram of twisted bilayer cuprates from a microscopic perspective" by Panda et al.

1. Problem Statement

Since the theoretical prediction by Can et al. (2021) that twisted bilayer cuprate superconductors could host a time-reversal symmetry breaking (TRSB) $d + id'$ topological state, experimental results have been conflicting.

• Zhao et al. (2023) observed a strong suppression of the Josephson critical current ($J_c$) at a $45^\circ$ twist angle, consistent with the predicted TRSB state.
• Zhu et al. (2021) observed a critical current that was largely independent of the twist angle, with no significant suppression at $45^\circ$, suggesting the absence of the topological state.

The central problem is to identify the physical parameters (twist angle, doping, interlayer tunneling strength) that control the stability of TRSB states and explain the discrepancies between these experimental setups. Previous continuum models lacked microscopic detail, while earlier lattice models often relied on simplified ansatzes for the order parameter.

2. Methodology

The authors employ a microscopic tight-binding lattice model to simulate two twisted square lattices (representing CuO$_2$ layers) with commensurate twist angles.

• Hamiltonian: The model includes:
  • In-layer nearest-neighbor ($t$) and next-nearest-neighbor ($t'$) hopping.
  • On-site chemical potential ($\mu$) tuned to fix the electron filling ($\nu$).
  • Attractive nearest-neighbor density-density interaction ($V$) in the singlet Cooper channel.
  • Interlayer tunneling ($g_{ij}$): Modeled phenomenologically as an exponential decay with distance ($g_{ij} = g_0 e^{-(r_{ij}-c)/\rho}$), allowing control over the tunneling strength ($g_0$).
• Symmetry Analysis: The twisted bilayer system possesses $D_4$ point group symmetry (lower than the $D_{4h}$ of isolated layers). The authors project the self-consistently solved order parameter onto the irreducible representations (irreps) of the $D_4$ group ($A_1, A_2, B_1, B_2, E$) to rigorously classify the superconducting states.
• Calculation: They solve the mean-field BCS gap equation self-consistently. To ensure unbiased convergence, they start with random complex initial guesses and apply gauge transformations to verify that any remaining imaginary components in real space are intrinsic (indicating TRSB) rather than artifacts.
• Observables:
  • Order parameter symmetry and magnitude.
  • Density of States (DOS) to identify gap opening and Van Hove singularities.
  • Sublattice spectral functions and Fermi surface topology.
  • Josephson critical current ($J_c$) calculated via the phase difference between layers.

3. Key Contributions

• Microscopic Phase Diagram: The authors map out the phase diagram as a function of twist angle ($\theta$), filling ($\nu$), and interlayer tunneling strength ($g_0$).
• Symmetry Classification: They move beyond simple $d+id'$ or $d+is$ labels, identifying specific $D_4$ symmetry combinations (e.g., $B_1 + iA_1$, $B_1 + iB_2$) and demonstrating how the symmetry of the order parameter changes with parameters.
• Mechanism for Experimental Discrepancy: They propose that the conflicting experimental results arise from differences in effective interlayer tunneling strength.
  • Bulk crystal junctions (Zhao et al.) likely have weaker effective tunneling, stabilizing a TRSB state with suppressed $J_c$.
  • Exfoliated flake junctions (Zhu et al.) likely have stronger effective tunneling (due to surface roughness or strain), driving the system into a TRSB state with finite $J_c$ or even a trivial Time-Reversal Symmetric (TRS) $s$-wave ($A_1$) state.
• Role of Van Hove Singularities: They establish a direct correlation between the position of the Van Hove singularity (VHS) relative to the chemical potential and the stability of TRSB states.

4. Key Results

A. Order Parameter and Symmetry

• 0° and 90° Twists: At weak tunneling, the system exhibits a pure $B_1$ ($d_{x^2-y^2}$) state. As tunneling ($g_0$) increases or filling changes, it transitions to a TRSB $B_1 + iA_1$ state (analogous to $d+is$) or a pure $A_1$ ($s$-wave) state.
• Near 45° Twists (e.g., 36.87°, 43.60°):
  • At low $g_0$, the system stabilizes a $B_1 + iB_2$ state (analogous to $d+id'$).
  • As $g_0$ increases, the system transitions to a $B_1 + iA_1$ state.
  • At very high $g_0$, the system can enter a fully gapped, trivial $A_1$ (s-wave) state.
• TRSB Identification: A true TRSB state is identified by a non-zero imaginary part of the order parameter in real space that cannot be removed by gauge transformation, accompanied by a fully gapped DOS (lifting of nodes).

B. Phase Diagram Evolution

• Tunneling Dependence: Increasing $g_0$ generally drives the system from TRSB states ($B_1 + iB_2$ or $B_1 + iA_1$) toward TRS states ($B_1$ or $A_1$).
• Filling Dependence: The stability of TRSB states is highly sensitive to doping. The VHS position shifts with doping; TRSB states are most stable when the VHS is close to the chemical potential.
• Temperature Dependence: TRSB states (multi-component) generally evolve into single-component TRS states as temperature increases, as different irreps have different critical temperatures ($T_c$).

C. Josephson Critical Current ($J_c$)

• Low Tunneling ($g_0 \approx 0.03$): $J_c$ decreases monotonically and significantly as the twist angle approaches $45^\circ$, reproducing the results of Zhao et al. This corresponds to the$B_1 + iB_2$ state.
• High Tunneling ($g_0 \approx 0.16$): $J_c$ remains significant and shows weak angle dependence, reproducing the results of Zhu et al. In this regime, the system stabilizes a $B_1 + iA_1$ state (which has a finite $J_c$) or a trivial $A_1$ state.
• Mechanism: The suppression of $J_c$ at $45^\circ$is linked to the specific symmetry mismatch in the$B_1 + iB_2$state, whereas the$B_1 + iA_1$and$A_1$ states allow for finite tunneling currents.

5. Significance

This work provides a unified microscopic explanation for the contradictory experimental observations in twisted bilayer cuprates. It demonstrates that the "topological" nature of the superconducting state is not fixed but is a tunable property dependent on the interlayer coupling and doping.

• Experimental Guidance: It suggests that the lack of a $45^\circ$suppression in flake-based experiments (Zhu et al.) is likely due to enhanced interlayer coupling (possibly from surface roughness or strain) pushing the system into a different phase ($B_1 + iA_1$or$A_1$) rather than the absence of TRSB physics entirely.
• Theoretical Advancement: By using a full microscopic lattice model without restrictive ansatzes, the authors reveal a richer phase diagram with multiple competing TRSB and TRS states, highlighting the critical role of the Van Hove singularity and Fermi surface topology in determining the ground state symmetry.
• Future Directions: The paper emphasizes the need for precise calculations of tunneling matrix elements (e.g., via Wannier functions) to fully validate the hypothesis that structural variations in different experimental setups drive the system into distinct regions of the phase diagram.