[Paper Review] Minimal model for low-energy electronic states of twisted bilayer graphene

Stephen Carr, Shiang Fang+2 more|arXiv (Cornell University)|Jan 10, 2019

Graphene research and applicationsMaterials Science3 citations

TL;DR

This paper proposes a minimal ten-band effective model for twisted bilayer graphene (tBLG) that incorporates lattice relaxation via $k \cdot p$ perturbation theory, combining DFT accuracy with continuum model efficiency. The inclusion of relaxation eliminates a second magic angle and accurately captures flat bands near the Fermi level at the first magic angle, fundamentally altering the low-energy electronic structure.

ABSTRACT

We introduce a physically motivated minimal model for the electronic structure of twisted bilayer graphene (tBLG), which incorporates the crucial role of lattice relaxation. Our model, based on $k \cdot p$ perturbation theory, combines the accuracy of DFT calculations through effective tight-binding Hamiltonians with the computational efficiency and complete control of the twist angle offered by continuum models. The inclusion of relaxation significantly changes the bandstructure at the first magic-angle twist corresponding to flat bands near the Fermi level (the states), and eliminates the appearance of a second magic-angle twist. We argue that the minimal model for the low-energy states of tBLG consists of ten bands, necessary to capture the changes in electronic states as a function of twist angle. We also provide information on the nature of these bands through their wavefunctions, which is closely tied to the features of the atomic relaxation.

Motivation & Objective

• To develop a minimal, physically motivated model for low-energy electronic states in twisted bilayer graphene (tBLG) that accurately captures the effects of lattice relaxation.
• To reconcile the accuracy of DFT calculations with the computational efficiency and tunability of continuum models in tBLG.
• To resolve discrepancies in prior models by including atomic relaxation, particularly around the first magic angle.
• To identify the minimal set of bands required to describe the evolution of electronic states across twist angles.
• To characterize the nature of the bands through their wavefunctions and their connection to atomic relaxation patterns.

Proposed method

• Formulating an effective $k \cdot p$ perturbation theory approach to derive low-energy Hamiltonians from first principles.
• Constructing effective tight-binding models calibrated to DFT results to ensure accuracy in band structure.
• Incorporating lattice relaxation effects into the continuum model through geometric and electronic structure corrections.
• Using twist-angle-dependent parameters to maintain full control over the system's symmetry and band dispersion.
• Deriving a minimal ten-band model that captures the essential physics across the full range of twist angles.
• Analyzing the wavefunctions of the bands to link their character to atomic relaxation features.

Experimental results

Research questions

• RQ1 How does lattice relaxation alter the band structure of twisted bilayer graphene near the first magic angle?
• RQ2 What is the minimal number of bands required to accurately describe the low-energy electronic states in tBLG across varying twist angles?
• RQ3 Does the inclusion of relaxation eliminate the spurious second magic-angle peak observed in prior models?
• RQ4 How do the wavefunctions of the low-energy bands relate to the atomic relaxation patterns in the system?
• RQ5 Can a minimal model combining DFT accuracy and continuum model efficiency be constructed for tBLG?

Key findings

• The inclusion of lattice relaxation eliminates the appearance of a second magic-angle twist, resolving a long-standing discrepancy in prior models.
• The minimal model for low-energy states in tBLG consists of exactly ten bands, necessary to capture the full evolution of electronic structure with twist angle.
• The bandstructure near the first magic angle exhibits flat bands near the Fermi level, consistent with DFT results and stabilized by relaxation effects.
• The wavefunctions of the bands are directly linked to the features of atomic relaxation, revealing a strong correlation between geometry and electronic character.
• The model achieves a balance between computational efficiency and accuracy, enabling precise control over the twist angle while retaining DFT-level fidelity.
• The ten-band model successfully describes the low-energy physics across the full range of twist angles, including the critical region near the first magic angle.