Operator ordering as an emergent gauge field in twisted bilayer graphene: singular spectral signatures at the magic angle

Technical Summary: Operator Ordering as an Emergent Gauge Field in Twisted Bilayer Graphene

Problem Statement
Scanning-tunneling spectroscopy (STS) of magic-angle twisted bilayer graphene (TBG) has established Van Hove singularities (VHS) at AA stacking sites consistent with Bistritzer–MacDonald (BM) predictions. However, the local density of states (LDOS) at the AB and BA stacking points—where the interlayer coupling $|T(r)|$ vanishes and coincides with the flat-band Dirac points of the chiral-limit spectrum—remains unanalyzed for fine structure. The standard BM spectrum predicts a single degenerate state at these points. This paper addresses a potential inconsistency in the continuum Dirac description of TBG: the standard Hamiltonian, when projected onto a spatially varying mass $m(r)$, fails to satisfy functional-analytic Hermiticity (unitarity of time evolution) despite being algebraically Hermitian. This failure manifests as a non-vanishing residual in the continuity equation, suggesting that the standard BM Hamiltonian is incomplete in the real-space continuum limit.

Methodology
The authors approach the problem by treating the interlayer coupling $T(r)$ as a position-dependent mass in a low-energy Dirac Hamiltonian.

1. Operator Ordering Derivation: In Section II, the authors derive a necessary Hermitian ordering correction, $H_{ord}$, for a Dirac fermion with spatially varying mass. By imposing three natural conditions—linearity in mass derivatives, membership in the Clifford algebra, and cancellation of the continuity equation residual—they uniquely determine the correction term:
   $$H_{ord} = -\frac{i\hbar v_F}{2} \boldsymbol{\alpha} \cdot \nabla \ln m(r)$$
   This term is shown to be distinct from gauge degrees of freedom and is required to restore functional-analytic Hermiticity in the $L^2(\mathbb{R}^2)$ Hilbert space.
2. Local Diagonalization: In Section III, the BM Hamiltonian in the chiral limit is locally diagonalized via a spatially dependent unitary transformation. This maps the interlayer coupling to a scalar effective mass $m_{eff}(r) = w|f(r)|$. The ordering term transforms into $H_{ord} = -\frac{i}{2} \boldsymbol{\sigma} \cdot \nabla \ln |f(r)|$.
3. Singularity Analysis: The authors analyze the behavior of $H_{ord}$ near the zeros of $f(r)$ (the AB/BA points). They demonstrate that $\nabla \ln |f(r)|$ develops a $1/r$ singularity, which dominates over regular Berry connection terms.
4. Topological and Spectral Analysis: Section IV establishes that this singularity acts as an emergent Aharonov–Bohm flux of $h/(2e)$, shifting the effective angular momentum $j \to j + 1/2$ and lifting the $j \leftrightarrow -j$ degeneracy. The zero-mode wavefunction is identified as a parabolic cylinder function $D_{-1/2}$.
5. Numerical Verification: Section V employs adiabatic interpolation and numerical diagonalization of the radial system to calculate the spectral splitting $\Delta_{split}$, confirming convergence and stability against geometric corrections.

Key Contributions and Results

• Emergent Flux and Splitting: The primary result is that the ordering correction generates an emergent Aharonov–Bohm flux of $h/(2e)$ at each AB/BA zero. This lifts the degeneracy of the lowest angular momentum states ($j = \pm 1/2$), creating a spectral splitting of $\Delta_{split} \approx 171$ meV at the magic angle.
• Scaling Law: The splitting scales as $\Delta_{split} \propto \sqrt{\theta}$ (where $\theta$ is the twist angle). This distinguishes the effect from correlation-driven gaps (scaling as $\theta$) and certain topological contributions (scaling as $1/\theta$).
• Spatial Localization: The lower energy peak ($j = -1/2$) is spatially confined within a radius $r_c \approx 2.1$ nm of the AB/BA points, a scale accessible to current STS instruments.
• Quantitative Consistency: The paper calculates that the cell-averaged broadening resulting from this splitting is $\approx 14$ meV. This value is quantitatively consistent with the $\approx 16$ meV discrepancy observed between existing STM data and the standard BM tight-binding prediction, offering a single-particle geometric explanation for this residual.
• Topological Protection: The emergent flux is protected by a $Z_2$ holonomy invariant ($W = -1$), ensuring it is a physical observable rather than a gauge artifact.

Significance and Claims
The paper claims that operator ordering is not merely a formal consistency requirement but a physically active mechanism capable of generating observable local spectral structures in moiré systems.

• Experimental Falsifiability: The authors propose that the predicted two-peak structure at AB/BA points, characterized by $\sqrt{\theta}$ scaling and gate-voltage independence, constitutes an immediate experimental test. This distinguishes the effect from correlation-driven phenomena.
• Resolution of Discrepancy: The work suggests that the ordering mechanism provides a natural explanation for the broadening of Van Hove singularities in global STS measurements, which was previously attributed solely to correlation effects.
• Theoretical Necessity: The paper argues that without the ordering term, the candidate zero mode at the AB/BA points fails to be square-integrable in the continuum formulation. Thus, the ordering term is necessary for a well-defined zero-mode subspace, complementing the topological analysis of Tarnopolsky et al.
• Broader Applicability: The framework is presented as applicable to any system with a spatially varying Dirac mass possessing isolated zeros, such as twisted transition-metal dichalcogenides and moiré-modulated topological insulators.

The authors emphasize that the effect was likely overlooked because standard BM derivations operate in momentum space where algebraic Hermiticity is automatic, and the functional-analytic inconsistency only arises in the real-space continuum limit with spatially varying mass.