Interplay between many-body correlations, strain and lattice relaxation in twisted bilayer graphene

This paper presents a unified theoretical framework demonstrating that the temperature-dependent electronic spectra and thermodynamic properties of magic-angle twisted bilayer graphene arise from the cooperative interplay between electron correlations and external symmetry-breaking effects induced by strain and lattice relaxation.

The Gist

Imagine you have a piece of graphene, which is essentially a single layer of carbon atoms arranged in a honeycomb pattern (like chicken wire). Now, take two of these sheets and stack them on top of each other, but twist one slightly relative to the other. This creates a new, giant pattern called a "moiré pattern," similar to the rippling effect you see when you hold two window screens slightly offset from each other.

When you twist them at a very specific "magic" angle, something magical happens: the electrons inside stop zooming around freely and get stuck, moving incredibly slowly. This "slow-motion" state allows the electrons to start talking to each other, creating exotic new states of matter, like superconductors (materials that conduct electricity with zero resistance) or insulators (materials that block electricity).

For years, scientists have been trying to build a perfect map of how these electrons behave. They had a theoretical map (the "Bistritzer-MacDonald model"), but when they compared it to real experiments, the map didn't match the terrain. The experiments showed strange features that the map couldn't explain.

The Problem: The Map Was Missing the "Real World" Mess
The authors of this paper realized that the theoretical map was too perfect. It assumed the graphene sheets were perfectly flat and perfectly aligned. But in the real world, things are messy:

1. Strain: When you make these samples, they get stretched or squished slightly, like a rubber band that's been pulled unevenly.
2. Relaxation: The atoms don't just sit still; they wiggle and shift to find the most comfortable position, like people in a crowded room shifting to make more space.

These "messy" real-world effects break the perfect symmetry of the system. The authors asked: What if we stop pretending the system is perfect and actually include these messy effects in our calculations?

The Solution: A Unified Theory
The team built a new, more realistic computer model that included:

• Electron Correlations: How the electrons push and pull on each other.
• Strain: The stretching and squishing of the lattice.
• Relaxation: The atoms shifting to their comfortable spots.

When they ran this new model, it suddenly matched the experimental data perfectly. Here is how they explained the three biggest mysteries using simple analogies:

1. The "Ghost" Signal (The Persistent Feature)

The Mystery: In experiments, scientists saw a strange "ghost" signal in the energy spectrum that appeared at the same energy level (about 10 meV) no matter how many electrons they added to the system. It was like a lighthouse that stayed on even when the tide changed. Previous models said this shouldn't exist.

The Explanation: Imagine the electrons live in a two-story apartment building (the flat bands). In the perfect model, the building is symmetrical. But when you add strain (stretching the building), the floors tilt.

• The tilt splits the building into two distinct sections: a "Bonding" section and an "Anti-bonding" section.
• Depending on whether you add electrons (electron-doped) or remove them (hole-doped), the electrons fill up one section completely and leave the other empty.
• The empty section acts like a "ghost." Even though no electrons are living there, the potential for them to jump there creates a signal. Because the tilt (strain) is fixed, this signal always appears at the same energy level, regardless of how many people (electrons) are in the building. This explains the "filling-independent" ghost signal.

2. The "Fidgety" Electrons (Entropy and Degeneracy)

The Mystery: Scientists measured the "entropy" (a measure of disorder or how many ways the electrons can arrange themselves). They found that at high temperatures, the electrons acted like they had 8 different "personas" (degeneracy), but as the temperature dropped, they suddenly settled into only 4 personas.

The Explanation: Think of the electrons as dancers.

• High Temperature: The music is loud and fast. The dancers are spinning wildly, and it's hard to tell who is who. They have 8 possible moves (high entropy).
• Low Temperature: The music slows down. Because of the strain (the tilted floor), the dance floor is now split into two separate zones. The dancers are forced to pick a zone.
• Once the temperature drops below a certain point, the dancers in one zone get "frozen" in place (they become inactive). They stop dancing. The remaining active dancers only have 4 moves left.
• The paper shows that the "freezing" of one group of electrons explains why the entropy drops so sharply. The strain acts like a gatekeeper, closing off half the dance floor.

3. The "One-Sided" Behavior (Particle-Hole Asymmetry)

The Mystery: The system behaves differently when you add electrons (positive charge) versus when you remove them (negative charge). For example, the material is more stable as a superconductor on one side, and the "squishiness" (compressibility) of the electron cloud changes more drastically on the other side.

The Explanation: Imagine a seesaw. In a perfect world, the seesaw is balanced. But lattice relaxation (the atoms shifting) puts a heavy weight on one side of the seesaw.

• This weight shifts the energy levels so that the "hole" side (removing electrons) and the "electron" side (adding electrons) are no longer mirror images.
• On the side where the atoms shifted closer together, the electrons can mix more easily with the surrounding material (hybridization). This makes that side "softer" and more metallic.
• On the other side, the electrons are more isolated and "harder."
• This explains why experiments show different behaviors for adding vs. removing electrons. The "relaxation" of the atoms broke the perfect symmetry, making the system lopsided.

The Big Picture

This paper is a triumph of "getting real." It tells us that to understand complex quantum materials, we can't just look at the ideal, mathematical version. We have to account for the fact that materials are physical objects that get stretched, squished, and relaxed.

By combining the messy reality of strain and relaxation with the complex dance of electron correlations, the authors finally cracked the code. They showed that these three factors working together create a unified story that explains why twisted bilayer graphene behaves the way it does, solving a puzzle that had stumped physicists for years. It's like finally realizing that the reason a car won't start isn't just a dead battery, but a combination of a flat tire, a cold engine, and a loose belt—and fixing all three makes it run perfectly.

Technical Summary

1. Problem Statement

Twisted bilayer graphene (TBLG) at the magic angle exhibits rich correlated electronic phases, including superconductivity and Mott-like insulating states. However, a unified theoretical framework that quantitatively explains diverse experimental observations has remained elusive. Specifically, existing models struggle to simultaneously account for:

• Persistent Spectral Features: The observation of a filling-independent spectral peak at $\sim \pm 10$ meV in Scanning Tunneling Microscopy (STM) and Quantum Twisting Microscopy (QTM) data.
• Local Moment Degeneracy: A temperature-dependent transition from 8-fold to 4-fold degenerate local moments, evidenced by entropy measurements.
• Particle-Hole Asymmetry: Significant differences in charge compressibility and superconducting stability between electron-doped and hole-doped sides, which standard symmetric models fail to reproduce.

The authors posit that these discrepancies arise because previous theoretical efforts often neglected the simultaneous effects of heterostrain (asymmetric lattice stretching) and lattice relaxation (atomic reconstruction), which break key symmetries of the system.

2. Methodology

The authors employ a comprehensive theoretical framework combining a microscopic model with advanced many-body techniques:

• Topological Heavy Fermion (THF) Model: They utilize the THF model, which maps the Bistritzer-MacDonald (BM) continuum model of TBLG onto a generalized periodic Anderson model. This model uses maximally localized Wannier orbitals:
  • $f$-orbitals: Localized in AA-stacked regions, representing the flat bands ($\sim 96\%$ of spectral weight).
  • $c$-orbitals: Dispersive orbitals representing remote bands, ensuring correct topological band representation via hybridization at the $\Gamma_M$ point.
• Incorporation of Strain and Relaxation:
  • Strain ($\epsilon$): Modeled as a uniaxial heterostrain (typically $\epsilon \approx 0.15\%$) that breaks $C_{3z}$ and $C_{2x}$ symmetries. This is introduced via a perturbation term $\delta H_\epsilon$ that splits the flat-band manifold.
  • Lattice Relaxation ($\Lambda$): Modeled via non-local tunneling terms and orbital-dependent chemical potential shifts ($\delta H_\Lambda$). This breaks particle-hole ($P$) symmetry, causing the upper flat bands to become more dispersive than the lower ones.
• Many-Body Solver:
  • Charge Self-Consistent Dynamical Mean-Field Theory (DMFT): Used to treat strong electronic correlations non-perturbatively.
  • CT-QMC Solver: The continuous-time quantum Monte Carlo solver w2dynamics is used to solve the impurity problem.
  • Interaction Parameters: Coulomb interactions are modeled with a double-gate screened potential. The dominant on-site $f$-$f$ interaction is $U \approx 58$ meV, significantly larger than the flat band width ($\sim 10$ meV), placing the system in the strong correlation regime.
  • Approximations: Spontaneous symmetry breaking is suppressed by working at $T \ge 11.6$ K (above the ordering temperature). Off-diagonal hybridization terms induced by strain are handled via a basis rotation to mitigate the fermionic sign problem in QMC, validated to be a negligible perturbation.

3. Key Contributions and Results

A. Origin of the Persistent Spectral Feature ($\sim 10$ meV)

• Mechanism: Strain induces a crystal-field-like splitting of the 8-fold degenerate flat-band manifold into two 4-fold degenerate subsets: bonding ($f_-$) and anti-bonding ($f_+$) bands. The splitting magnitude is $\sim 7\text{--}10$ meV for typical strain levels.
• Result: Depending on the doping (electron vs. hole), one of these split sectors becomes "inactive" (fully occupied or empty) and pinned away from the Fermi level. The other remains "active" and crosses the Fermi level.
• Observation: The "persistent feature" observed in STM/QTM at $\pm 10$ meV corresponds to excitations within this inactive sector. Because the inactive sector is pinned by the strain-induced gap, its spectral weight does not shift with the chemical potential, creating a filling-independent peak. This matches experimental data where the peak position depends on strain magnitude but not filling.

B. Entropy and Local Moment Degeneracy

• Mechanism: The strain-induced splitting creates an energy gap between the bonding and anti-bonding flat bands.
• Result:
  • At low temperatures ($T < \text{strain gap}$), the system effectively freezes out one sector, reducing the local moment degeneracy from 8-fold to 4-fold. This leads to entropy vanishing at the Charge Neutrality Point (CNP).
  • At higher temperatures ($T > \text{strain gap}$), thermal excitation restores the 8-fold degeneracy, resulting in finite entropy at CNP.
• Agreement: The calculated entropy curves as a function of filling and temperature show excellent quantitative agreement with recent experimental measurements, resolving the controversy over the degeneracy of local moments.

C. Particle-Hole Asymmetry in Compressibility

• Mechanism: Lattice relaxation breaks particle-hole symmetry by shifting the $f$-orbital energies relative to the $c$-Dirac bands. This results in a higher density of states (DOS) for $c$-electrons on the hole-doped side compared to the electron-doped side.
• Result:
  • Hybridization: The higher $c$-DOS on the hole-doped side enhances $f$-$c$ hybridization, which suppresses local moments and "softens" the charge compressibility variations.
  • Electron-doped side: Lower $c$-DOS leads to weaker hybridization, stronger localization, and sharper, more pronounced peaks in inverse compressibility.
• Agreement: The model reproduces the experimentally observed strong particle-hole asymmetry, where charge compressibility varies more strongly on the electron-doped side, and explains why superconductivity is often more stable on the hole-doped side (due to softer compressibility and reduced correlation effects).

4. Significance

This work provides a unified microscopic framework that resolves several long-standing puzzles in TBLG physics:

1. Symmetry Breaking as a Driver: It establishes that external symmetry breaking (strain) and internal structural relaxation are not mere perturbations but essential ingredients that dictate the emergent electronic behavior.
2. Inactive Sector Freezing: The concept of an "inactive" charge sector that freezes out at low temperatures offers a new perspective on the low-energy effective theories of TBLG, suggesting that collective instabilities (like superconductivity) must be understood within a constrained low-energy landscape.
3. Quantitative Agreement: By quantitatively matching spectroscopy, compressibility, and entropy data without adjustable parameters beyond standard material values, the study validates the THF model extended with strain and relaxation as the correct effective theory for magic-angle TBLG.
4. General Implications: The findings highlight the extreme sensitivity of correlated quantum materials to lattice deformations, suggesting that controlling strain and relaxation could be a viable route for tuning electronic phases in moiré materials.

In conclusion, the paper demonstrates that the interplay between many-body correlations, strain, and lattice relaxation is the key to unlocking the complex phenomenology of twisted bilayer graphene, reconciling previously conflicting experimental signatures into a single coherent picture.