Graphene’s ‘magic angle’ reveals quantum stable states at key electron counts

Researchers at the University of California, Berkeley and Lawrence Berkeley National Laboratory have developed a quantum embedding workflow to model the behavior of magic-angle twisted bilayer graphene (MATBG). The work reveals robust insulating states at specific electron counts and a fragile semimetal state at others, with energy splittings at the meV scale making correlated phases highly sensitive to material parameters.

Solving models using Hartree-Fock and coupled cluster methods, the team recovered Kramers intervalley coherent states at charge neutrality and electron doping, and observed a weak Kekulé modulation at hole doping. Their approach connects microscopic electronic structure to scanning tunneling microscopy signatures in MATBG.

Ab Initio Quantum Embedding for Magic-Angle Graphene

The constrained random phase approximation (cRPA) is central to a new workflow designed to map the electronic properties of magic-angle twisted bilayer graphene (MATBG), allowing researchers to derive interacting flat-band Hamiltonians from density functional theory calculations. This approach addresses the sensitivity of MATBG’s correlated phases to both material parameters and modeling choices, a consequence of the material’s narrow moiré bands exhibiting meV-scale energy splittings. The workflow combines cRPA screening with a controlled double-counting subtraction procedure, refining the accuracy of the resulting Hamiltonian.

A key element of this method is an automated gauge-fixing procedure based on the selected columns of the density matrix (SCDM), ensuring compatibility with symmetry-resolved many-body calculations. This automated process is critical because subtle variations in the initial conditions can significantly impact the resulting electronic structure, particularly at the low energies characteristic of MATBG.

The researchers focused on even-integer fillings, specifically ν = 0, +2, and -2, to identify and characterize candidate quantum phases through mean-field energetics, band structures and potential signatures in scanning tunneling microscopy. Results from Hartree-Fock (HF) and coupled cluster singles and doubles (CCSD) calculations reveal robust insulating Kramers intervalley coherent (KIVC) states at charge neutrality (ν = 0) and at electron doping (ν = +2).

However, the most notable finding emerges at hole doping (ν = -2), where the model predicts a fragile semimetal characterized by a weak $\sqrt{3}$× $\sqrt{3}$ Kekulé modulation. This semimetallic state also exhibits enhanced intervalley-scattering peaks in the Fourier-transformed local density of states, suggesting a unique scattering mechanism within the material. The study acknowledges that establishing the robustness of this hole-side semimetal requires further investigation, including enlarged active-space calculations with refined gauge fixing and subtraction techniques.

Researchers suggest that extrinsic effects, such as heterostrain, may play a decisive role in stabilizing gapped states in experimental devices, and that the choice of subtraction parameters can qualitatively affect the hole-side spectrum at meV energy scales. “Ab Initio Quantum Embedding of Twisted Bilayer Graphene” details the first-principles simulations and cRPA screening used in the analysis, as well as implementation details of the adaptive gauge fixing procedure.

Supplementary materials also include structured initializations and symmetry diagnostics, providing a comprehensive overview of the methodology. Future work, the authors note, will explore strained geometries, a more detailed comparison with continuum models, and ab initio treatment of electron-phonon interactions. The team’s approach offers a pathway toward more accurate and reliable modeling of MATBG’s complex electronic behavior, potentially bridging the gap between theoretical predictions and experimental observations.

Constrained RPA Screening and Double-Counting Subtraction

Constructing accurate models of magic-angle twisted bilayer graphene (MATBG) requires careful attention to the interplay between material properties and computational methods, as the narrow moiré bands exhibit energy splittings on the scale of just a few milli-electronvolts. This approach combines constrained random phase approximation (cRPA) screening with a controlled subtraction of double-counted interactions, an important step in accurately representing the system’s electronic behavior. The researchers found that the choice of reference density for the double-counting subtraction significantly impacts the resulting band structure, revealing a particle-hole asymmetry that influences the stability of different quantum states.

Specifically, using a reference density derived from the charge-neutral state proved quantitatively and qualitatively important for accurately modeling the system’s behavior. This careful handling of double-counting is essential, as it removes interaction effects already implicitly included in the initial density functional theory band energies through the Hartree and exchange-correlation potentials.

The model’s construction begins with defining an interacting Hamiltonian, projecting onto the active manifold and embedding interactions within an ab initio environment. This projection-and-screening construction functions as a quantum embedding, explicitly treating the flat-band subspace while incorporating the remaining electronic degrees of freedom through screened interactions and subtraction terms. “We find that this semimetallic behavior is largely controlled by the double-counting subtraction procedure in the effective Hamiltonian,” the paper notes, highlighting the critical role of this step in determining the system’s properties.

The team’s approach to double-counting subtraction uses a mean-field Hamiltonian of the form $H_{sub}(k)$ = J[$P_{ref}$](k) + K[$P_{ref}$](k), where J and K represent the Coulomb and exchange contributions constructed from a reference one-body density $P_{ref}$. This careful consideration of interactions and their removal allows for a more accurate representation of the correlated electronic behavior within the twisted bilayer graphene system.

SCDM Gauge-Fixing Aligns with MATBG Symmetry

The automated gauge-fixing procedure aligns the computational framework with the inherent symmetries of magic-angle twisted bilayer graphene (MATBG). Researchers found that an uncontrolled Bloch-orbital gauge can obstruct both the initialization and interpretation of these solutions, hindering accurate predictions of material properties. By exploiting the translational symmetry of the moiré superlattice, simulations use a sampling of $N_k$ k-points within the moiré Brillouin zone, enabling a detailed examination of the system’s electronic structure.

This method builds upon a downfolded interacting Hamiltonian derived from Kohn-Sham density functional theory (KS-DFT) of a relaxed, unstrained structure, combining it with constrained random phase approximation (cRPA) screening and controlled double-counting subtraction. The resulting framework allows for symmetry-resolved many-body calculations, which are important for understanding the subtle interplay of quantum effects in MATBG.

The SCDM procedure functions by transforming Bloch orbitals into sublattice- and valley-polarized orbitals, denoted ${\psi }_{\text{(σ},\text{τ)},k}$ where σ represents the sublattice (An or B) and τ represents the valley (K or K’). These symmetries dictate that the Hamiltonian satisfies certain conditions, which are reflected in the transformed basis. The work also details the construction of the downfolded interacting Hamiltonian and the associated double-counting subtraction, outlining the implementation details in supplementary material. The process begins by defining creation and annihilation operators for each orbital, then transforming the Bloch orbitals to the sublattice- and valley-polarized basis.

KIVC Insulating States at ν = 0 and ν = +2

The stability of these KIVC states is particularly notable given the narrow moiré bands within MATBG, exhibiting energy splittings measured in the millielectronvolt (meV) range, which makes correlated phases unusually susceptible to subtle shifts in material properties or computational choices. Further analysis of ν = + 2 mirrors the behavior observed at ν = 0, with KIVC remaining the lowest-energy state and CCSD calculations yielding only minor corrections to the HF results.

However, the electronic landscape shifts considerably at ν = -2, where the calculations point to a fragile semimetallic state. Unlike predictions from some continuum model studies, which typically anticipate a robust insulator at this doping level, the research indicates bands crossing the Fermi level, resulting in gapless behavior. Our work provides a first-principles route for connecting microscopic electronic structure, screened interactions, subtraction choices, and scanning tunneling microscopy signatures in MATBG.

While previous studies have predicted the KIVC insulating state as the ground state at even-integer filling in the absence of strain or electron-phonon interactions, these results suggest a more detailed relationship between doping and the emergence of distinct quantum phases within MATBG. The observed contrast between the behavior at ν = -2 and the robust insulating states at ν = 0 and ν = + 2 highlights the importance of precise modeling and the sensitivity of these correlated electron systems to subtle changes in their environment.

Particle-Hole Asymmetry at ν = ±2 Fillings

The behavior of electrons in magic-angle twisted bilayer graphene (MATBG) diverges significantly depending on whether the material is doped with electrons or holes, with calculations revealing a fragile semimetallic state at a hole-doping level of ν = -2. While the underlying band structure appears nearly symmetrical, the effective interacting Hamiltonian exhibits a pronounced particle-hole asymmetry at both ν = +2 and ν = -2, stemming from momentum-dependent single-particle renormalizations introduced during the computational process.

This asymmetry is particularly evident when comparing the behavior at ν = +2, where robust insulating states are predicted, to the semimetallic solution found at ν = -2. The research demonstrates that the double-counting subtraction procedure within the effective Hamiltonian largely controls this semimetallic behavior, a finding that contrasts with some previous continuum and downfolded descriptions of unstrained MATBG.

The team traced the discrepancy to the construction of the subtraction Hamiltonian from an active-space reference density consistent with the KS-DFT filling, which introduces momentum-dependent Hartree renormalization. This renormalization shifts the valence manifold upward near the moiré γ point, reducing the indirect gap on the hole side and driving the semimetallic phase. The choice of reference density is therefore an important factor influencing the predicted low-energy spectrum at the meV energy scales characteristic of MATBG.

Previous studies often employ a particle-hole symmetric “average” reference density, simplifying the mean-field problem, but this approach deviates significantly from the projected KS density within the active subspace when the conduction flat bands are empty in the KS reference state. The findings emphasize that even at these small energy scales, the choices made during the computational process, specifically, the subtraction procedure and the selection of reference densities, can qualitatively alter the predicted electronic properties of MATBG.

MATBG’s Kekulé Modulation at ν = -2 Doping

The observed particle-hole asymmetry at ν = ± 2, despite an initially nearly particle-hole symmetric Kohn-Sham density functional theory (ks-dft) band structure, stems from momentum-dependent single-particle renormalizations. These renormalizations arise from subtraction terms used in constructing the effective interacting Hamiltonian, terms derived from reference densities consistent with the ks-dft filling.

The resulting ν = -2 semimetal exhibits a weak $\sqrt{3}$× $\sqrt{3}$ Kekulé modulation and enhanced intervalley-scattering peaks in the Fourier-transformed local density of states. While KIVC remains the lowest-energy state among the structured ansatzes, the CCSD correction remains small, indicating a relatively stable, though delicate, configuration.

More like this

PhysicsCaSb₂ superconductor’s quantum vortex clusters defy theory

Quantum Research NewsQuantum asymmetry can grow faster in broken systems

Quantum Research NewsEntanglement analysis matches order found in Rydberg atom chains

Quantum Research NewsNoisy quantum cooling