Strain-induced excitonic instability in twisted bilayer graphene

Hector Ochoa

DOI: https://doi.org/10.1103/PhysRevB.102.201107

2020-10-23

Abstract:The low-energy bands of twisted bilayer graphene form Dirac cones with approximate electron-hole symmetry at small rotation angles. These crossings are protected by the emergent symmetries of moiré patterns, conferring a topological character to the bands. Strain accumulated between layers (heterostrain) shifts the Dirac points both in energy and momentum. The overlap of conduction and valence bands favors an excitonic instability of the Fermi surface close to the neutrality point. The spontaneous condensation of electron-hole pairs breaks time reversal symmetry and the separate conservation of charge within each valley sector. The order parameter describes interlayer circulating currents in a Kekulé-like orbital magnetization density wave. Vortices in this order parameter carry fermion numbers owing to the underlying topology of the bands. This mechanism may explain the occurrence of insulating states at neutrality in the most homogenous samples, where uniform strain fields contribute both to stabilizing the relative orientation between layers and to the formation of an excitonic gap.

Mesoscale and Nanoscale Physics

What problem does this paper attempt to address?

This paper aims to explore the strain - induced excitonic instability in twisted bilayer graphene (TBG). Specifically, the research focuses on the phenomena near the charge neutral point. The following are the main problems that the paper attempts to solve:

1. **Explain the formation mechanism of the insulating state**:
- When the twist angle is close to the "magic angle", it has been experimentally observed that TBG exhibits an insulating state near the charge neutral point. The origin of these insulating states is not fully understood. This paper proposes a possible mechanism: a uniform strain field causes the Fermi surface to be unstable, which in turn forms an excitonic gap, thus explaining these insulating states.

2. **The influence of strain on the electronic structure**:
- Strain will introduce significant changes in the band structure, especially in the low - energy region. The article explores how strain affects the position and energy of the Dirac cone, and points out that such changes may cause the Fermi surface to be unstable, making electron - hole pairs more likely to pair and form an excitonic condensate.

3. **Excitonic instability and symmetry breaking**:
- The article analyzes in detail how excitonic instability breaks time - reversal symmetry and inter - valley charge conservation. Specifically, the spontaneous electron - hole pair condensation breaks these symmetries, forming a Kekulé - type orbital magnetization density wave and, in some cases, carrying a fermion number.

4. **Topological effects and vortex phenomena**:
- Due to the topological nature of the band gap, vortices in the excitonic order parameter can carry a fermion number. This phenomenon is similar to the vortices in the Kekulé bond order, but has different quantum number characteristics.

5. **The possibility of experimental observation**:
- The proposed theoretical model predicts some phenomena that can be directly observed in scanning tunneling microscopy (STM) and other experimental techniques, such as orbital magnetization density wave and neutral valley supercurrents. In addition, this mechanism should also exhibit specific symmetry - breaking features in magnetotransport measurements.

### Summary of mathematical formulas

- **Energy shift of the Dirac cone**:
\[
\hat{H}_{\text{str}} = g_A^2 \omega_{ii} \hat{\Lambda}_z \hat{\Gamma}_z + g_E^{(1)} \left[ (\omega_{xx} - \omega_{yy}) \hat{\Sigma}_x \hat{\Gamma}_z - 2\omega_{xy} \hat{\Sigma}_y \hat{\Gamma}_z \right] + g_E^{(2)} \left[ 2\omega_{xy} \hat{\Sigma}_x \hat{\Lambda}_z + (\omega_{xx} - \omega_{yy}) \hat{\Sigma}_y \hat{\Lambda}_z \right]
\]

- **Expression for the excitonic gap**:
\[
\Delta_0(T = 0)\approx\frac{2g_A^2\bar{u}}{\sinh\left(\frac{2}{\nu_F V}\right)}
\]
where \(\nu_F\) is the density of states per valley and spin, and \(V\) is the screening potential of the Coulomb interaction.

- **Critical temperature**:
\[
T_c\approx\frac{2e^\gamma E g_A^2\bar{u}}{\pi k_B}e^{-\frac{2}{\nu_F V}}
\]
where \(\gamma_E\) is the Euler constant.

Through these formulas and theoretical analysis, the paper proposes a reasonable mechanism to explain twisted bilayer graphene.