🔬 mesoscale physics

Relaxation effects on Hartree-Fock ground states in twisted bilayer graphene at even integer fillings

This paper presents a novel systematic derivation of a relaxed single-particle continuum model for magic angle twisted bilayer graphene, demonstrating that incorporating structural relaxation drives the Hartree-Fock ground state at filling $\nu = -2$ into a semi-metallic phase due to particle-hole asymmetry and enhanced Hartree potentials, thereby corroborating recent *ab initio* findings while offering explanations for the discrepancy with experimental observations.

Original authors: Tianyu Kong, Alexander B. Watson, Lin Lin, Mitchell Luskin, Kevin D. Stubbs

Published 2026-09-03

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Tianyu Kong, Alexander B. Watson, Lin Lin, Mitchell Luskin, Kevin D. Stubbs

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

Deep within the world of materials science, there is a class of substances known as twisted bilayer graphene. Imagine two sheets of carbon atoms, arranged in a honeycomb pattern like a chicken wire fence, stacked directly on top of one another. When researchers twist the top sheet by a very specific, tiny angle—about 1.1 degrees—the atoms in the two layers no longer line up perfectly. Instead, they create a new, larger pattern called a moiré pattern, which acts like a giant, repeating landscape of hills and valleys for electrons moving through the material. At this specific "magic" angle, the electrons slow down dramatically, becoming heavy and sluggish. This slowdown allows the electrons to interact with one another in powerful ways, leading to strange and exciting behaviors, such as superconductivity, where electricity flows with zero resistance. Understanding exactly how these electrons arrange themselves and move is crucial for unlocking the potential of these materials for future technology.

For years, scientists have tried to predict how these electrons behave using computer models. A standard approach involves simplifying the complex physics of the atoms into a manageable set of rules, focusing only on the electrons that are moving the slowest. However, these models often assume that the atoms stay perfectly rigid in their positions, ignoring the fact that atoms are actually free to wiggle and shift to find the most comfortable, low-energy arrangement. This shifting, known as structural relaxation, changes the landscape the electrons travel through. In a new study, researchers have developed a more precise way to include this atomic movement in their models. They created a detailed simulation that accounts for how the atoms relax and shift, and then used this improved map to calculate how the electrons settle into their lowest energy states at specific filling levels, where the number of electrons matches the number of available spots in the material's energy bands.

The team focused on two specific scenarios: when the material has two extra electrons per unit of the pattern and when it is missing two electrons. In the older, rigid models, scientists predicted that the material would act as an insulator in both cases, meaning the electrons would get stuck in place and electricity could not flow. However, when the researchers applied their new model that included the atomic shifts, the results changed dramatically for the case with two missing electrons. Instead of becoming an insulator, the material behaved as a semi-metal, a state where electrons can move freely, albeit with some resistance. This difference was not a small fluctuation; it was a fundamental change in the material's character. The researchers found that this shift occurred because the atomic relaxation caused the electron clouds to squeeze into tighter, more concentrated spots within the pattern. This concentration strengthened the electrical forces between the electrons, effectively closing the gap that usually stops them from moving.

Interestingly, this effect was not symmetrical. When the researchers looked at the case with two extra electrons, the material remained an insulator, just as the older models had predicted. The reason for this difference lies in the uneven nature of the energy landscape created by the relaxed atoms. The energy levels for electrons moving in one direction were different from those moving in the opposite direction, a property known as particle-hole asymmetry. This asymmetry meant that the squeezing of the electron clouds helped the material stay insulating when it had extra electrons, but pushed it toward a conducting state when it was missing electrons. The researchers confirmed that their findings held true regardless of the specific mathematical adjustments used to handle the complex interactions between electrons, suggesting the result is robust.

These findings align with recent, more complex computer simulations that used different methods to study the same material, both of which also found that the material becomes semi-metallic when missing two electrons. This stands in contrast to many previous studies that relied on the simpler, rigid models and predicted an insulating state. The new work provides a clearer explanation for why these different models yield different results, pointing directly to the physical movement of the atoms and the resulting concentration of electron clouds. However, the researchers note that real-world experiments have mostly observed insulating behavior at these filling levels, not the semi-metallic state their simulations predict. This suggests that other physical factors, such as uneven stretching of the material or the influence of other electron bands that were not included in this specific study, might be at play in actual experiments. The study does not claim to have solved the mystery of why experiments differ from simulations, but it offers a vital piece of the puzzle by showing how the physical relaxation of the atomic structure fundamentally alters the electronic properties of twisted bilayer graphene.

Technical Summary: Relaxation Effects on Hartree-Fock Ground States in Twisted Bilayer Graphene at Even Integer Fillings

Problem Statement
Magic angle twisted bilayer graphene (MATBG) exhibits strongly correlated electronic phases, including superconductivity and insulating states. A standard approach to modeling these phenomena involves projecting Coulomb interactions onto effective single-particle models involving flat bands and remote bands. However, the predictions of these effective models are sensitive to their derivation, particularly regarding mechanical structural relaxation. In relaxed MATBG, atoms displace from their monolayer equilibrium positions to minimize total energy, altering the stacking configuration (enlarging "AB" regions relative to "AA"). While previous studies have attempted to incorporate relaxation by phenomenologically tuning hopping parameters (e.g., the ratio of "AA" to "AB" hopping in the Bistritzer-MacDonald model), there is a need for a systematic derivation of a continuum model that accounts for relaxation while preserving the Lagrangian frame. Furthermore, recent ab initio density functional theory (DFT) studies and specific scanning tunneling microscopy (STM) observations suggest semi-metallic behavior at filling $\nu = -2$ , a result that contradicts many effective model predictions which typically find insulating phases at even integer fillings ( $\nu = \pm 2$ ).

Methodology
The authors develop a novel, systematic derivation of an effective single-particle continuum model for MATBG that incorporates structural relaxation using multiple-scale analysis.

Derivation in the Lagrangian Frame: Unlike previous approaches that may switch frames or rely on ad-hoc parameter tuning, this work remains entirely in the Lagrangian frame. The authors model structural relaxation by minimizing a continuum energy functional coupling linear elasticity to a generalized stacking fault energy (GSFE) parametrized by DFT.
Multiple-Scale Expansion: By introducing a small parameter $\epsilon$ $ϵ$ representing the ratio of the atomic lattice scale to the moiré scale, the authors expand the tight-binding Hamiltonian. They derive a relaxed continuum Hamiltonian ( $H_{relax}$ $H_{r e l a x}$ ) that includes:
- Intralayer terms with second-order corrections and linearized rotation corrections.
- Interlayer potentials that depend on relaxation-induced displacements.
- Non-local terms (momentum-dependent) arising from the relaxation, which are identified as the primary source of particle-hole asymmetry in the single-particle dispersion.
Many-Body Modeling: The Coulomb interaction is projected onto the flat bands of this relaxed model. The authors compute Hartree-Fock (HF) many-body ground states at even integer fillings ( $\nu = 0, \pm 2$ ).
Comparison: The results are compared against the standard Bistritzer-MacDonald (BM) model where relaxation is simulated by tuning the hopping ratio $\kappa = w_0/w_1$ . The study employs both "charge neutral" and "average" double-counting subtraction schemes to ensure robustness.

Key Contributions

Systematic Derivation: The paper provides a rigorous derivation of a relaxed single-particle continuum model for moiré materials that stays in the Lagrangian frame, simplifying the treatment of structural relaxation compared to prior methods.
Identification of Mechanisms: The authors isolate two distinct mechanisms by which relaxation alters the many-body ground state:
1. Single-Particle Dispersion: The relaxed model exhibits significant particle-hole asymmetry in the single-particle band structure due to non-local terms.
2. Wavefunction Concentration: The flat band wavefunctions in the relaxed model are more spatially concentrated than in the BM model. This concentration enhances the Hartree potential.
Resolution of Discrepancies: The work clarifies why recent ab initio DFT studies found semi-metallic phases at $\nu = -2$ , whereas tuning the BM model often predicts insulating phases.

Results

Particle-Hole Asymmetry: The relaxed model displays strong particle-hole asymmetry in its single-particle dispersion, unlike the symmetric BM model.
Ground State at $\nu = -2$ : At filling $\nu = -2$ , the relaxed model yields a semi-metallic Hartree-Fock ground state. This is driven by the enhanced Hartree potential resulting from wavefunction concentration. Specifically, the down-spin Hartree term is negative and larger in magnitude for the relaxed model (away from the $\Gamma_m$ point), shifting the bands downward and causing them to overlap near $\Gamma_m$ .
Ground State at $\nu = +2$ : At filling $\nu = +2$ , the system remains insulating. The same wavefunction concentration effect occurs, but the Hartree potential shift is positive (due to the different occupation), and this upward shift is canceled by the particle-hole asymmetric single-particle dispersion, preserving the gap.
Consistency: These findings hold across different double-counting subtraction schemes ("charge neutral" and "average").
Comparison with Experiments: The paper notes that while their results align with recent ab initio DFT studies and DFT-parametrized continuum models, they contradict experimental observations which generally report insulating behavior at $\nu = \pm 2$ .

Significance and Claims
The authors claim that their work provides a clearer understanding of the origin of semi-metallic phases at $\nu = -2$ in MATBG, attributing it fundamentally to wavefunction concentration and particle-hole asymmetry induced by structural relaxation. They argue that previous effective models based solely on tuning hopping parameters in the BM model fail to capture these specific relaxation-induced effects.

The paper modestly acknowledges that the discrepancy between their theoretical prediction (semi-metallic at $\nu = -2$ ) and experimental observation (insulating) remains unresolved. They suggest that effects omitted from their current model—such as out-of-plane layer corrugation, the inclusion of remote bands, heterostrain (layer-specific strain), and twist angle disorder—may be necessary to recover the insulating behavior seen in experiments. The authors propose that their derived framework is capable of incorporating these additional physical effects in future work.

More like this