🔬 mesoscale physics

Commensuration torques in double-moiré twisted trilayer hexagonal boron nitride and graphene heterostructures

This study utilizes large-scale atomistic relaxations to demonstrate that double-moiré commensuration in twisted trilayer hBN and graphene/hBN heterostructures induces local energy minima and torque sign reversals, establishing a system-dependent mechanism for twist-angle stabilization driven by enhanced stacking domain overlap and Coulomb interactions.

Original authors: Youngju Park, Nicolas Leconte, Prathap Kumar Jharapla, Md Shaifullah, E. H. Hwang, Jeil Jung

Published 2026-05-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Youngju Park, Nicolas Leconte, Prathap Kumar Jharapla, Md Shaifullah, E. H. Hwang, Jeil Jung

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have three thin, flat sheets of material stacked on top of each other, like a very delicate sandwich. In this study, scientists looked at two types of sandwiches: one made entirely of hexagonal boron nitride (hBN) and another made of alternating layers of graphene and hBN.

These sheets aren't perfectly aligned; they are slightly twisted relative to one another. When you twist two flat sheets, they create a giant, repeating pattern called a "moiré" pattern (think of the wavy lines you see when you hold two window screens slightly askew).

The researchers were investigating what happens when you have two of these twist interfaces in a three-layer stack. They wanted to know: Do these layers slide around freely, or do they get "stuck" in specific positions?

Here is the breakdown of their findings using simple analogies:

1. The "Double-Moiré" Effect

In a standard two-layer twist, the layers might slide easily or get stuck at one specific angle. But in this three-layer "double-moiré" system, the scientists found a special rule: The layers prefer to lock together when the twist angle of the top layer matches the twist angle of the bottom layer.

Think of it like a dance. If the bottom dancer spins clockwise at a certain speed, and the top dancer spins clockwise at the exact same speed, they find a "sweet spot" where they feel most comfortable and stable. If they spin at different speeds, they feel unstable and want to adjust until they match up again.

2. The "Rubber Band" Torque

The paper uses the concept of "torque" to explain this locking. Imagine the layers are connected by invisible rubber bands.

When the angles match: The rubber bands are relaxed. This is the "energy minimum" (the most comfortable spot).
When the angles don't match: The rubber bands stretch. This creates a force (torque) that pulls the layers back toward the matching angle.
The "Sign Reversal": If you twist the top layer slightly too far to the left, the rubber band pulls it right. If you twist it slightly too far to the right, the rubber band pulls it left. This "pulling back" is what the scientists call angle locking.

3. The Two Types of Sandwiches

The researchers tested two different "recipes" for their three-layer sandwiches, and they behaved slightly differently:

The "All-BN" Sandwich (Homolayer):
In a stack made entirely of boron nitride, the layers are naturally very similar. Here, the "matching angle" (where the top and bottom twists are equal) creates a local energy minimum.
- Analogy: Imagine a valley in a mountain range. The layers love to sit in this valley because it's comfortable. However, if you pushed them hard enough, they could roll all the way down to the very bottom of the mountain (perfect alignment, zero twist). The "matching angle" is just a very comfortable resting spot, but not the absolute lowest point.
The "Mixed" Sandwich (Heterolayer):
In stacks mixing graphene and boron nitride, the atoms don't line up perfectly because the materials are slightly different sizes (lattice mismatch).
- Analogy: Here, the "matching angle" valley is so deep that it becomes the bottom of the mountain. In some cases, the layers actually prefer to stay twisted at this specific angle (around 0.6 degrees) rather than aligning perfectly straight. It's as if the "sweet spot" became the only place the layers wanted to live.

4. Why Do They Lock? (The Puzzle Piece Analogy)

Why does this locking happen? The scientists looked at the atomic level.

In the locked (commensurate) state: The "low-energy" spots (where atoms fit together nicely, like puzzle pieces) on the bottom interface line up perfectly with the "low-energy" spots on the top interface. This creates a large, continuous zone of comfort.
In the unlocked (incommensurate) state: The puzzle pieces on the top and bottom don't line up. The comfortable spots are scattered and mixed with uncomfortable spots. The system has to "average out" the discomfort, making it less stable overall.

5. The Role of Electricity

Since boron nitride is a polar material (it has a slight electrical charge), the researchers checked if electricity changed the game. They found that while electrical forces make the "locking" even stronger (deeper valleys), the basic mechanism remains the same. The layers still want to match their twist angles to find stability.

Summary

The paper concludes that in these three-layer twisted systems, there is a strong, natural tendency for the layers to "lock" their twist angles together.

If the materials are the same, this locking creates a stable resting spot, though perfect alignment is still the ultimate goal.
If the materials are different, this locking can become the most stable state of all, preventing the layers from ever aligning perfectly.

This discovery helps scientists understand how to control these materials, suggesting that by twisting them to specific angles, we can create stable structures that stay put, rather than sliding around randomly.

Technical Summary: Commensuration Torques in Double-Moiré Twisted Trilayer Hexagonal Boron Nitride and Graphene Heterostructures

Problem Statement
The study addresses the energetic stability and angle-locking mechanisms in "supermoiré" systems, specifically double-moiré trilayer heterostructures composed of graphene (G) and hexagonal boron nitride (hBN). While single-moiré systems (twisted bilayers) are well-characterized, the energetics of stacking three layers with two distinct moiré interfaces remain less explored. Previous work on twisted trilayer graphene (t3G) suggested that double-moiré commensuration (where the twist angles between adjacent layers satisfy specific relations, e.g., $\theta_{12} = -\theta_{23}$ ) leads to local energy minima and torque sign reversals, indicating a restoring tendency toward commensurate configurations. However, it was unclear how these mechanisms translate to systems involving hBN, which possesses a polar lattice structure (B and N atoms) and different lattice constants compared to graphene. The central problem is to determine whether double-moiré commensuration acts as a universal mechanism for twist-angle stabilization across various G/hBN combinations and to quantify the resulting binding energies and torque constants.

Methodology
The authors employ large-scale atomistic relaxations using the LAMMPS software package with the FIRE minimization scheme. The total energy ( $E_{tot}$ ) is minimized by balancing elastic (intralayer) and potential (interlayer) energy components.

Potentials: Intralayer interactions for hBN and graphene are modeled using the extended Tersoff potential (BN-ExTeP) and the REBO2 potential (CH.rebo), respectively. Interlayer interactions are described using DRIP parametrizations informed by EXX-RPA calculations.
System Construction: The study constructs commensurate supercells and large rational approximants for incommensurate twist angles using six integers to enforce commensurability conditions.
Control Parameters: To isolate double-moiré effects, the bottom-layer rotation angle ( $\theta_{12}$ ) is fixed, while the top-layer rotation angle ( $\theta_{32}$ ) is varied as a continuous control parameter.
Torque and Binding Energy: The torque ( $k$ ) is defined as the angular derivative of the total energy with respect to $\theta_{32}$ . Sign reversals in torque across a commensurate angle indicate a local energy minimum. Binding energies ( $E_b$ ) are calculated relative to a smooth reference curve derived from cubic-spline interpolation of data points away from the commensuration dip.
Electrostatics: Benchmark calculations include pairwise Coulomb interactions (using coul/shield and coul/long pair styles) to assess the impact of hBN's polarity on the energetics.

Key Contributions and Results

Twisted Trilayer hBN (t3BN) Homolayers:
- In t3BN systems (where lattice constants are matched), double-moiré commensuration occurs when $\theta_{12} = \theta_{32}$ .
- The system exhibits local energy minima at these commensurate angles, accompanied by torque sign reversals (negative on the left, positive on the right of the commensurate angle).
- Binding energies are found to be $\sim 0.2–0.3$ meV/atom. These minima arise from the enhanced spatial overlap of low-energy stacking domains (e.g., AB and BA regions) between the two moiré interfaces.
- In incommensurate phases, the overlap of low-energy domains is reduced, leading to spatial averaging of interfacial energies and enhanced superlubricity.
- The global energy minimum for t3BN remains at zero twist (perfect alignment), with commensurate angles providing only local minima.
Graphene/hBN Heterolayers (Lattice Mismatched):
- Symmetric Systems (G/BN/G, BN/G/NB, BN/G/BN): For these systems, commensuration also occurs at $\theta_{12} = \theta_{32}$ . However, a distinct behavior emerges at small twist angles ( $\sim 0.6^\circ$ ). Here, the double-moiré commensurate configuration becomes the global energy minimum, surpassing the perfectly aligned state. Binding energies in this regime are $\sim 0.1–0.2$ meV/atom.
- Asymmetric Systems (G/BN/BN, G/BN/NB): Due to structural asymmetry, commensuration occurs at $\theta_{12} \neq \theta_{32}$ (specifically $\theta_{32} \approx 0.58^\circ$ for $\theta_{12} = \pm 1.20^\circ$ ). These systems also exhibit local energy minima with binding energies of $\sim 0.1–0.15$ meV/atom.
- Torque Characteristics: In lattice-mismatched systems, the torque sign reversal persists over a broader angular range compared to lattice-matched systems, indicating enhanced angular stability around the commensurate configuration.
Role of Electrostatics:
- Benchmark calculations including short-range (coul/shield) and long-range (coul/long) Coulomb interactions confirm that the local energy minima at double-moiré commensuration persist.
- Long-range interactions slightly deepen the energy dips and alter the corrugation profiles (shifting stability between "breathing" and "bending" modes), but the fundamental mechanism of angle locking remains unchanged.

Significance
The paper establishes double-moiré commensuration as a general, system-dependent mechanism for twist-angle stabilization in trilayer heterostructures. The key findings are:

Universal Mechanism: The torque sign reversal and local energy minima driven by the overlap of low-energy stacking domains are robust features observed across both homolayer (t3BN) and heterolayer (G/hBN) systems.
System-Dependent Global Stability: Unlike t3G and t3BN, where commensuration yields only local minima, certain G/hBN heterostructures (specifically at small twist angles near $0.6^\circ$ ) exhibit the double-moiré commensurate state as the global energy minimum. This suggests that stable, angle-locked bilayer-like states can be realized in trilayer systems without perfect alignment.
Experimental Relevance: The calculated binding energies ( $\sim 0.1–0.3$ meV/atom) are sufficiently large to suggest that these commensurate phases should be experimentally accessible. The results provide a theoretical basis for interpreting metastable states observed in graphene flakes on hBN substrates near $0.6^\circ$ , attributing them to double-moiré effects rather than edge-induced coincidences or rigid-substrate mechanisms.
Guidance for Design: The study posits that calculating torque constants and binding energies is a necessary step to determine the robustness and experimental realizability of specific double-moiré configurations.