Magnetic response of twisted bilayer graphene

Simon Becker,Jihoi Kim,Xiaowen Zhu

2024-01-04

Abstract:In this article, we analyse the Bistritzer--MacDonald (BM) model (also known as the continuum model) of twisted bilayer graphene (TBG) with an additional external magnetic field. We provide an explicit semiclassical asymptotic expansion of the density of states (DOS) in the limit of strong magnetic fields. The explicit expansion of the DOS enables us to study magnetic response properties such as magnetic oscillations which includes Shubnikov-de Haas and de Haas-van Alphen oscillations as well as the integer quantum Hall effect. In particular, we elucidate the role played by different types of interlayer tunnelings ($AA^{\prime}$/$BB^{\prime}$ vs. $AB^{\prime}$/$BA^{\prime}$) in the study of the DOS, and magnetic properties.

Mathematical Physics,Mesoscale and Nanoscale Physics,Strongly Correlated Electrons,Spectral Theory,Quantum Physics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to study the magnetic response characteristics of twisted bilayer graphene (TBG) in an external magnetic field. Specifically, the authors analyzed the semi - classical asymptotic expansion of the density of states (DOS) in the Bistritzer - MacDonald (BM) model under a strong magnetic field. Through this expansion, they were able to study magnetic oscillation phenomena (such as Shubnikov - de Haas and de Haas - van Alphen oscillations) and the integer quantum Hall effect. In addition, they also explored the influence of different types of inter - layer tunneling (AA′/BB′ vs. AB′/BA′) on DOS and magnetic properties.

### Key point overview

1. **Research background**:
- Twisted bilayer graphene exhibits a superconducting phase at certain "magic angles", which is one of the most exciting discoveries in condensed matter physics in recent years.
- The BM model is an effective 4×4 matrix - valued Hamiltonian used to describe two twisted graphene layers and predicts the first magic angle.

2. **BM model**:
- The BM model consists of two twisted Dirac operators \( H_\theta^D \) and \( H_{-\theta}^D \), representing two isolated graphene layers respectively.
- The model also includes a tunneling potential term \( T_\theta(x) \), representing two different types of inter - layer tunneling:
- **Chiral tunneling** \( U(x/\lambda_\theta) \) and \( U(-x/\lambda_\theta) \), located near the vertices of each moiré hexagon, with a tunneling strength of \( \alpha_1 \), similar to AB′ and BA′ stackings.
- **Antichiral tunneling** \( V(x/\lambda_\theta) \), located near the center of each moiré hexagon, with a tunneling strength of \( \alpha_0 \), similar to AA′/BB′ stackings.

3. **Asymptotic expansion of DOS under a strong magnetic field**:
- The authors derived the explicit asymptotic expansion of DOS for the chiral and antichiral models under a strong magnetic field.
- For the chiral model, the peaks of DOS at Landau levels are enhanced, which leads to the enhancement of magnetic oscillations and the quantum Hall effect.
- For the antichiral model, the peaks of DOS at Landau levels are weakened.

4. **Physical applications**:
- These results can be used to understand the Landau level dependence at small twist angles and the influence of strong pseudo - magnetic fields (generated by physical strain).
- The results also provide a strict theoretical basis for DOS peaks at small twist angles.

### Mathematical formulas

- **Hamiltonian of the BM model**:
\[
H_\theta=\begin{pmatrix}
H_\theta^D & 0 \\
0 & H_{-\theta}^D
\end{pmatrix}+\begin{pmatrix}
0 & T_\theta \\
(T_\theta)^* & 0
\end{pmatrix}
\]
where \( T_\theta(x)=\begin{pmatrix}
\alpha_0 V(x/\lambda_\theta) & \alpha_1 U(-x/\lambda_\theta) \\
\alpha_1 U(x/\lambda_\theta) & \alpha_0 V(x/\lambda_\theta)
\end{pmatrix} \).

- **Asymptotic expansion of DOS**:
- **Chiral model**:
\[
\tilde{\text{Tr}}(f(H_c))=\left[\frac{B}{\pi} f(\lambda_n, B)+\frac{|n|}{2\pi} \text{...}
\]