In-plane magnetic response and Maki parameter of alternating-twist multilayers

This paper analytically demonstrates that alternating-twist multilayer graphene systems exhibit a unique hierarchy of magic angles and distinct in-plane orbital magnetic responses—negligible for odd-layer systems but highly angle-dependent and significantly enhanced for even-layer systems—leading to large in-plane Maki parameters that suggest the potential for diverse superconducting phases.

The Gist

Imagine you have a stack of thin, magical sheets of carbon (graphene). If you twist two of these sheets slightly relative to each other, something amazing happens: the electrons inside slow down so much they form a "flat" energy landscape. This is called Magic-Angle Twisted Bilayer Graphene (TBG), and it's famous for becoming a superconductor (conducting electricity with zero resistance) at very low temperatures.

Now, imagine stacking four or five of these sheets, but twisting them in an alternating pattern (like a zig-zag: twist left, twist right, twist left...). This is what the authors of this paper studied. They wanted to know: How do these multi-layer stacks react when you push on them with a magnetic field from the side?

Here is the breakdown of their findings using simple analogies:

1. The "Magic" Decoupling Trick

The authors used a mathematical trick (a "unitary transformation") to simplify the problem. Think of a complex 4-layer sandwich. Instead of analyzing the whole messy sandwich, they realized they could mathematically "unzip" it.

• For 4 layers (Tetralayer): It acts like two separate, independent 2-layer sandwiches sitting next to each other, but each with a different twist angle.
• For 5 layers (Pentalayer): It acts like two separate 2-layer sandwiches plus one lonely, single sheet floating in the middle.

2. The Magnetic "Push" (In-Plane Response)

Usually, when you push a magnet sideways against a superconductor, the electrons try to fight back by creating their own magnetic field (orbital response). In standard magic-angle graphene, this fight is huge and loud.

The authors asked: Does this happen in our multi-layer stacks?

• The 5-Layer Case (The Odd One Out):
  Because there is a "lonely" single sheet in the middle, the magnetic pushes from the top and bottom layers cancel each other out perfectly. It's like two people pushing a car from opposite sides with equal force; the car doesn't move.

  • Result: The magnetic response is negligible (almost zero). This holds true for any stack with an odd number of layers.

• The 4-Layer Case (The Split Personality):
  This is where it gets interesting. The 4-layer stack has two different "magic angles" where superconductivity happens. The authors found that the stack behaves completely differently depending on which angle you tune it to:

  • At the "Large" Magic Angle: The two internal 2-layer systems are out of sync. One pushes hard, the other pushes weakly, and they cancel each other out. The result is a tiny magnetic response (about 1% of the usual strength).
  • At the "Small" Magic Angle: The two internal systems sync up perfectly. They both push in the same direction, amplifying the effect. The result is a huge magnetic response (about 3.6 times stronger than a standard 2-layer sheet).

3. The "Maki Parameter" (The Superconductor's Shield)

Superconductors have a limit: if the magnetic field gets too strong, the superconductivity breaks. This limit is usually set by how the electron spins react (the "Pauli Limit").

However, if the electrons also have a strong orbital (circular motion) reaction, they can actually help the superconductor survive stronger magnetic fields. The authors introduced a new "Maki Parameter" to measure this.

• Think of it like a shield: A high Maki parameter means the orbital motion is building a strong shield against the magnetic field.
• Their Finding: In the 4-layer system at the "Small" magic angle, this shield is incredibly strong (Maki parameter $\approx$ 7). This means the superconductivity here is very robust and might survive magnetic fields that would destroy other superconductors.

The Big Picture Takeaway

This paper reveals that alternating-twist multilayers are not just bigger versions of the original 2-layer graphene; they are entirely new creatures.

• Odd layers (3, 5, etc.): They are "quiet" regarding magnetic fields. This is great for experiments because it means scientists can easily see the subtle effects of electron spins without the "noise" of orbital motion masking them.
• Even layers (4, 6, etc.): They are "loud" and tunable. You can switch them between a "quiet" mode and a "super-shield" mode just by changing the twist angle slightly.

Why does this matter?
It suggests that we can engineer superconductors with specific, custom-made properties. We might be able to create a material that is superconducting and incredibly resistant to magnetic fields, or one that allows us to study the fundamental nature of electron pairing without interference. It's like having a volume knob for the magnetic behavior of a superconductor.

Technical Summary

1. Problem Statement

The discovery of superconductivity in magic-angle twisted bilayer graphene (MATBG) has established flat-band engineering as a powerful tool for inducing unconventional phase transitions. A critical question in these systems is the nature of the superconducting pairing mechanism. In typical experiments, an in-plane magnetic field is applied to test the Pauli limit (the critical field $B_P \approx 1.86 T_c$ required to break singlet Cooper pairs via spin alignment).

However, in two-dimensional systems, the in-plane magnetic field also induces an orbital magnetic response. If this orbital susceptibility ($\chi_{orb}$) is large, it can mask the spin susceptibility ($\chi_{spin}$), making it difficult to determine if the violation of the Pauli limit is due to unconventional pairing (e.g., triplet pairing) or simply strong orbital effects. While MATBG exhibits a large orbital response, it is unclear how this behaves in alternating-twist multilayer graphene (e.g., trilayer, tetralayer, pentalayer), where the electronic structure can be mapped onto decoupled twisted bilayer systems. Specifically, the authors investigate whether the orbital response in these multilayers is suppressed (allowing direct observation of spin effects) or enhanced, and how this varies across the multiple "magic angles" present in these systems.

2. Methodology

The authors employ an analytical approach based on the unitary transformation introduced by Khalaf et al. (2019), which maps alternating-twist $N$-layer graphene systems onto a hierarchy of decoupled twisted bilayer graphene (TBG) systems and, for odd $N$, a single decoupled layer.

• Theoretical Framework:

  • Layer-Resolved Conductivity: The electromagnetic response is decomposed into sheet currents for individual layers using the Kubo formula and Ohm's law in the frequency domain.
  • Mapping: The $N$-layer Hamiltonian is block-diagonalized into effective TBG systems with distinct renormalized twist angles ($\theta_{N,m}$).
  • Magnetic Field Definition: An in-plane magnetic field $B$ is defined via the layer-discriminated electric fields ($E_\ell$) using the discrete Maxwell-Faraday law. This creates a linear gradient of electric fields across the layers.
  • Maki Parameter: The authors introduce an in-plane Maki parameter ($\alpha_M$), defined as the ratio of the difference in orbital susceptibility between the normal and superconducting states to the Pauli spin susceptibility ($\alpha_M = \Delta\chi_{orb} / \chi_P$). This quantifies the deviation from the standard Pauli limit.

• Specific Systems Analyzed:

  • Odd $N$ (3, 5): Mapped to decoupled TBGs + 1 Single Layer Graphene (SLG).
  • Even $N$ (4): Mapped to two decoupled TBGs with different effective twist angles.

3. Key Contributions and Results

A. Orbital Magnetic Response in Odd-Layer Systems ($N=3, 5$)

• Result: The in-plane orbital magnetic response is negligibly small for systems with an odd number of layers.
• Mechanism: The response arises solely from "cross terms" (coupling between the decoupled TBG subsystem and the detached SLG). Due to the large mismatch in Fermi velocities between the flat bands of the TBG and the linear bands of the SLG, the phase space for these cross-correlations is restricted, rendering the contribution negligible.
• Implication: In odd-layer systems, the measured critical magnetic field directly reflects the spin susceptibility of Cooper pairs, as the orbital contribution does not mask it.

B. Orbital Magnetic Response in Even-Layer Systems ($N=4$)

The tetralayer system exhibits a strongly angle-dependent response, behaving differently at its two distinct magic angles:

1. First Magic Angle ($\theta_{4,m}^{(1)} \approx 1.70^\circ$):
   • Corresponds to the larger effective twist angle ($\phi \theta_m$).
   • The orbital susceptibility is suppressed to $\approx 0.01$ times the value of standard TBG.
   • The magnetization profile changes sign across the central layers, leading to regions of vanishing magnetization.
   • Significance: Similar to odd layers, the spin susceptibility can be isolated here.
2. Second Magic Angle ($\theta_{4,m}^{(2)} \approx 0.65^\circ$):
   • Corresponds to the smaller effective twist angle ($\phi^{-1} \theta_m$).
   • The orbital susceptibility is enhanced, reaching $\approx 3.6$ times the value of magic-angle TBG.
   • The magnetization is uniform and large.
   • Significance: The large orbital response will dominate the magnetic measurement, potentially masking the spin susceptibility and complicating the interpretation of the pairing symmetry.

C. The In-Plane Maki Parameter ($\alpha_M$)

The authors calculate $\alpha_M$ to quantify the modification of the Pauli limit:

• TBG: $\alpha_M$ reaches values up to 2 near the magic angle, indicating a significant orbital contribution.
• Tetralayer ($N=4$):
  • At the first magic angle: $\alpha_M \approx 0.02$ (Pauli limit holds; spin effects dominate).
  • At the second magic angle: $\alpha_M \approx 7$ (Pauli limit is heavily violated; orbital effects dominate).
• Conclusion: The same alternating-twist multilayer system can host distinct magnetic behaviors and potentially distinct superconducting phases depending on which magic angle is tuned.

4. Significance

• Experimental Guidance: The paper provides a clear theoretical prediction for experimentalists: in alternating-twist multilayers, the ability to probe the intrinsic spin pairing mechanism via the Pauli limit depends critically on the layer number ($N$) and the specific twist angle tuned.
• Novel Physics in Even Layers: It reveals that even-layer systems are not monolithic; they can exhibit diamagnetic-like suppression of orbital response at one magic angle and giant enhancement at another.
• Superconducting Phases: The drastic difference in orbital response suggests that the two magic angles in a tetralayer system might support contrasting superconducting phases (e.g., one dominated by spin fluctuations, the other by orbital mechanisms).
• Generalization: The finding that odd-layer systems have negligible in-plane orbital response offers a robust platform for studying pure spin-susceptibility effects in moiré superconductors.

In summary, this work analytically demonstrates that alternating-twist multilayers offer a tunable platform where the orbital magnetic response can be switched from negligible to dominant, fundamentally altering the interpretation of magnetic field experiments and the nature of superconductivity in these materials.