Modeling of electronic dynamics in twisted bilayer graphene

Tianyu Kong,Diyi Liu,Mitchell Luskin,Alexander B. Watson

2024-02-26

Abstract:We consider the problem of numerically computing the quantum dynamics of an electron in twisted bilayer graphene. The challenge is that atomic-scale models of the dynamics are aperiodic for generic twist angles because of the incommensurability of the layers. The Bistritzer-MacDonald PDE model, which is periodic with respect to the bilayer's moiré pattern, has recently been shown to rigorously describe these dynamics in a parameter regime. In this work, we first prove that the dynamics of the tight-binding model of incommensurate twisted bilayer graphene can be approximated by computations on finite domains. The main ingredient of this proof is a speed of propagation estimate proved using Combes-Thomas estimates. We then provide extensive numerical computations which clarify the range of validity of the Bistritzer-MacDonald model.

Mathematical Physics,Mesoscale and Nanoscale Physics,Analysis of PDEs,Numerical Analysis,Quantum Physics

What problem does this paper attempt to address?

This paper mainly discusses the modeling problem of electron dynamics in twisted bilayer graphene. The challenge in the study lies in the non-periodicity of the atomic-scale model under general twist angles due to the incompatibility between the two layers. The paper introduces the Bistritzer-MacDonald partial differential equation model, which can rigorously describe the dynamic behavior within a specific parameter range.

The author first proves that the dynamics of the incompatible twisted bilayer graphene tight-binding model can be approximated by calculations on a finite domain. This proof relies on the velocity propagation estimate of Combes-Thomas estimation. Then, they use numerical calculations to elucidate the validity range of the Bistritzer-MacDonald model.

The main findings of the study include:
1. Within a specific parameter range (e.g., when the ratio of twist angle to wave packet spectral width is smaller than a certain threshold), the Bistritzer-MacDonald model indeed captures the key dynamic characteristics of the tight-binding model, such as predicting the group velocity of wave packets.
2. Changes in parameters generally follow the expected error scaling relation, but in some cases (e.g., with increasing twist angle), the error is smaller than expected.
3. Even with relatively large twist angles, as long as the wave packet spectral width and certain parameters remain small, the Bistritzer-MacDonald model can maintain high accuracy.

These findings contribute to understanding the observed quantum many-body phase transitions in experiments, especially superconductivity in twisted bilayer graphene, and provide a deeper understanding of the applicable range of effective continuous models.