Twist-configured moire-moire reconstruction governs diverse commensurate double-moire phases in twisted bilayer graphene on h-BN

Problem and Motivation
Van der Waals heterostructures featuring multiple moiré superlattices present a fundamental challenge: how do distinct moiré patterns interact and reconstruct when coupled through shared atomic layers? While structural reconstruction is well-established in single-interface systems like twisted bilayer graphene (tBG) or graphene on hexagonal boron nitride (Gr/h-BN), the physics of "multi-moiré" systems remains less understood. Specifically, in tBG aligned with h-BN, two qualitatively different moiré patterns coexist—a triangular lattice from the graphene–graphene interface and a hexagonal lattice from the graphene–h-BN interface. Previous experiments have hinted at local commensuration in limited parameter spaces, but it remains unclear whether this is a special case or a general principle. Key questions regarding the roles of global twist configuration (helical vs. alternate), twist angle, and strain in determining local stacking registry and domain formation have largely remained unexplored.

Methodology
The authors investigate tBG/h-BN hetero-trilayers using a combination of experimental imaging and theoretical modeling:

• Sample Fabrication: They utilize a dry-transfer method with gold-assisted flipping to create tBG/h-BN stacks in two distinct global twist configurations: "helical" (successive twists in the same direction) and "alternate" (successive twists in opposite directions). Twist angles are controlled via AFM-based cutting of graphene flakes.
• Experimental Imaging: Conductive Atomic Force Microscopy (C-AFM) is employed in contact mode. By applying a bias voltage, the technique maps stacking-dependent local electronic properties, allowing for the visualization of moiré superlattices and the identification of specific stacking registries (e.g., AA, AB, BA, AA', AB', BA').
• Theoretical Modeling: Continuum-model simulations are used to calculate binding-energy landscapes, displacement fields, and electronic band structures. These models incorporate lattice relaxation, hetero-strain, and homo-strain to predict stacking registries and commensurate phases.

Key Contributions and Results

1. Discovery of Twist-Configuration-Dependent Registry:
   The study reveals that the global twist configuration deterministically selects the local spatial registry between the two moiré lattices.

   • In the helical configuration, the vertices of the triangular tBG moiré lattice (AA stacking) align with the centers of the hexagonal Gr/h-BN moiré cells (AB' stacking).
   • In the alternate configuration, the tBG AA vertices align with the vertices of the Gr/h-BN moiré (AA'/BA' stacking).

2. Mechanism of Local Rotation Matching:
   The authors identify "local rotation matching" as the governing mechanism. Through continuum simulations, they demonstrate that the shared graphene layer undergoes rotational relaxation to minimize energy. The system selects a registry where the local rotational directions induced by the tBG interface and the Gr/h-BN interface reinforce each other.

   • In the helical case, both interfaces drive a counterclockwise rotation near specific sites, favoring the AA-on-AB' registry.
   • In the alternate case, the reversed twist of the h-BN layer inverts the rotational field, favoring an AA-on-AA'/BA' registry to maintain cooperative relaxation.

3. Diverse Commensurate Double-Moiré Phases:
   The research establishes that moiré–moiré reconstruction stabilizes a rich family of commensurate domains across multiple length scales, governed by twist angles and strain:

   • C3z-Symmetric Phases: Quantized integer period ratios (e.g., 1:1, 2:1, 3:1 in helical; $\sqrt{7}$:1 in alternate) form large, ordered domains.
   • Strained Symmetry-Broken Phases: Even when strain breaks the $C_{3z}$ symmetry of the lattice, commensurate domains are preserved. The local rotation matching mechanism remains robust, allowing the system to form symmetry-modified structures that act as counterparts to the symmetric phases.
   • Phase Diagrams: The authors construct phase diagrams mapping the formation of these commensurate domains as a function of twist angles ($\theta_{tBG}$, $\theta_{Gr/BN}$) and strain, showing that reconstruction extends the window of commensurability beyond discrete, exact matches.

4. Mesoscopic Order and Dynamics:
   Large-area C-AFM images reveal well-ordered sub-micrometer domains separated by straight boundaries. The Gr/h-BN moiré acts as a periodic template organizing the tBG moiré. Furthermore, the domain boundaries exhibit collective sliding dynamics, where chains of tBG AA domains shift coherently parallel to the boundary without disrupting the surrounding structure.

5. Electronic Consequences:
   Theoretical calculations of the electronic band structure show that the twist-configured moiré–moiré stacking significantly modifies the electronic properties. The Gr/h-BN potential acts as a perturbation that opens gaps at charge neutrality and stabilizes topological flat bands with distinct Chern numbers at twist angles well below the magic angle. The global twist configuration (helical vs. alternate) serves as a tunable parameter for controlling the topology of these minibands.

Significance
The paper establishes "moiré–moiré reconstruction" based on the local rotation matching mechanism as a general framework for engineering structural and electronic order in multilayer moiré materials. It demonstrates that the interaction between coexisting moiré lattices is not merely a superposition but a coupled phenomenon that dictates local registry, stabilizes diverse commensurate phases under strain, and enables the design of topological flat bands. This work provides a predictive foundation for understanding and manipulating complex multilayer van der Waals heterostructures, moving beyond single-interface moiré physics to a regime where distinct moiré lattices can be deliberately coupled to create new quantum states.