Graphene's 'magic angle' reveals quantum stable states at specific 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). This study reveals robust insulating states at specific electron counts and a fragile semimetal state at others, with energy splittings at the millielectronvolt scale making correlated phases highly sensitive to material parameters.

The team employed Hartree-Fock and coupled cluster methods to recover Kramers intervalley coherent states at charge neutrality and electron doping, observing a weak Kekulé modulation at hole doping. Their approach connects microscopic electronic structure to scanning tunneling microscopy signatures in MATBG.

Central to this workflow is the constrained random phase approximation (cRPA), which maps the electronic properties of 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 material parameters and modeling choices, due to the material’s narrow moiré bands exhibiting meV-scale energy splittings.

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

The study 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).

At hole doping (ν = -2), the model predicts a fragile semimetal characterized by a weak 3 × 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. Further investigation is needed to establish the robustness of this hole-side semimetal.

Researchers suggest that extrinsic effects, such as heterostrain, may play a decisive role in stabilizing gapped states in experimental devices. 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.

Constructing accurate models of 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 cRPA screening with a controlled subtraction of double-counted interactions, an important step in accurately representing the system’s electronic behavior.

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.

The automated gauge-fixing procedure aligns the computational framework with the inherent symmetries of 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.

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 cRPA screening and controlled double-counting subtraction. The resulting framework allows for symmetry-resolved many-body calculations, 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 ψ (σ , τ), 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 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.

The behavior of electrons in 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.