Moiré-driven equilibrium of perturbations in moiré systems

Technical Summary: Moiré-driven equilibrium of perturbations in moiré systems

Problem Statement
Twisted bilayer graphene (TBG) and other moiré materials exhibit remarkable electronic properties, such as unconventional superconductivity and strongly correlated phases, which are intrinsically linked to the formation of flat bands near the "magic angle." While the properties of these flat bands (topology, origin, and recurrence of magic angles) have been extensively studied, the behavior of perturbations at the magic angle remains largely unexplored. In experimental settings, moiré heterostructures are ubiquitous with perturbations arising from substrates (e.g., hexagonal boron nitride), strain, or relaxation fields. These perturbations can significantly alter electronic properties, such as gapping Dirac cones or shifting them in momentum and energy. A key open question is how these layer-specific perturbations behave when the layers are strongly coupled by the moiré potential, particularly near the magic angle where interlayer hybridization is maximal. Recent observations suggest that electronic properties near the magic angle become insensitive to interlayer differences, but the underlying mechanism was not fully understood.

Methodology
The authors employ a combination of analytical first-order perturbation theory and extensive numerical calculations based on the continuum model of TBG.

1. Continuum Model: The system is modeled using a Hamiltonian that includes Dirac Hamiltonians for the top ($H_t$) and bottom ($H_b$) layers coupled by a moiré potential $U$. The moiré potential is expanded to the leading order, characterized by effective hopping energies $u$ (AA stacking) and $u'$ (AB/BA stacking).
2. Perturbation Analysis: Generic perturbations $P_\ell$ (scalar potential $V_\ell$, gauge potential $A_\ell$, and mass term $m_\ell$) are introduced into individual layers. The study focuses on the regime where perturbations are small compared to the moiré energy scale ($|P_\ell| \ll \hbar v k_\theta$).
3. First-Order Perturbation Theory: The Hamiltonian is truncated at the first shell of the moiré Brillouin zone. The authors derive renormalized expressions for the Fermi velocity and the perturbations themselves. This involves mixing the perturbations from both layers via the moiré coupling terms.
4. Numerical Validation: The analytical predictions are tested against numerical simulations of the full continuum model, varying the moiré coupling strength $\alpha$ (related to the twist angle) and the ratio of hopping parameters $\lambda = u/u'$. The study also extends to strain-extended continuum models to account for perturbations that simultaneously affect the moiré potential itself.

Key Contributions and Results
The central finding is the identification of a "moiré-driven equilibrium" where perturbations initially localized on one layer are redistributed between the two layers, eventually reaching an equilibrium state near the magic angle.

• Renormalization of Perturbations: The moiré coupling renormalizes layer-specific perturbations. For a perturbation $P_\ell$ in layer $\ell$, the effective perturbation $P^*_\ell$ becomes a linear combination of the initial perturbations in both layers ($P_t$ and $P_b$).
• Equilibrium Conditions:
  • Mass Perturbations: For mass terms (which open a gap), the renormalized gaps at the two Dirac points become equal when $3\alpha^2(1-\lambda^2) = 1$. This condition is met near the first magic angle ($\alpha \approx 0.58$). At this point, the equilibrium gap is a weighted average of the initial gaps, approaching a perfect average $(\Delta_t + \Delta_b)/2$ in the chiral limit ($\lambda=0$).
  • Scalar Perturbations: For scalar potentials (energy shifts), the equilibrium condition is $3\alpha^2(1+\lambda^2) = 1$. The equilibrium shift is exactly the average of the initial shifts ($V^*_{eq} = (V_t + V_b)/2$), independent of $\lambda$.
  • Gauge Perturbations: For gauge potentials (momentum shifts), the Dirac points can reach an equilibrium where their momentum shifts are equal in magnitude but opposite in sign ($\delta k_t = -\delta k_b$), leading to a collapse of the Dirac points within the moiré Brillouin zone. This collapse occurs specifically at the first magic angle.
• Numerical Observations: Numerical results confirm that as the moiré coupling $\alpha$ increases towards the magic angle, the system naturally evolves toward these equilibrium states.
  • The equilibrium tendency is robust even when the perturbation magnitude varies.
  • Beyond the first magic angle, the redistribution of perturbations does not vanish but follows a non-trivial oscillatory behavior, with bandgaps and energy shifts exhibiting complex closing and reopening dynamics.
  • The layer polarization of electronic states vanishes at the equilibrium point, indicating that the layers are strongly hybridized and the initial layer-specific nature of the perturbation is effectively masked.
• Robustness under Moiré Perturbation: The study demonstrates that this equilibrium behavior persists even when the moiré potential itself is perturbed (e.g., by strain or hBN substrates). In these cases, the perturbations still flow toward an equilibrium state, explaining why layer-resolved properties are difficult to disentangle in experiments.

Significance
The paper claims that the "magic angle" regime is part of a broader phenomenon termed moiré-driven equilibrium. This mechanism provides a natural explanation for the experimental observation that electronic properties near the magic angle appear insensitive to interlayer differences or the specific layer on which a perturbation is applied. By showing that strong moiré coupling inherently mixes and equalizes perturbations, the work extends the concept of the magic angle to a general regime governed by this equilibrium. The findings suggest that the masking of layer-resolved properties is a generic feature of moiré materials arising from strong interlayer hybridization, rather than a specific artifact of pristine TBG. This framework helps interpret recent experiments on strained TBG and TBG on hBN, where the electronic response is identical regardless of whether strain or potential differences are applied to one or both layers.