skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, April 16 until 2:00 AM ET on Friday, April 17 due to maintenance. We apologize for the inconvenience.


Title: Electronic Observables for Relaxed Bilayer Two-Dimensional Heterostructures in Momentum Space
Momentum space transformations for incommensurate two-dimensional electronic structure calculations are fundamental for reducing computational cost and for representing the data in a more physically motivating format, as exemplified in the Bistritzer--MacDonald model [Proc. Natl. Acad. Sci. USA, 108 (2011), pp. 12233--12237]. However, these transformations can be difficult to implement in more complex systems such as when mechanical relaxation patterns are present. In this work, we aim for two objectives. First, we strive to simplify the understanding and implementation of this transformation by rigorously writing the transformations between the four relevant spaces, which we denote real space, configuration space, momentum space, and reciprocal space. This provides a straightforward algorithm for writing the complex momentum space model from the original real space model. Second, we implement this for twisted bilayer graphene with mechanical relaxation effects included. We also analyze the convergence rates of the approximations and show the tight-binding coupling range increases for smaller relative twists between layers, demonstrating that the 3-nearest neighbor coupling of the Bistritzer--MacDonald model is insufficient when mechanical relaxation is included for very small angles. We quantify this and verify with numerical simulation.  more » « less
Award ID(s):
1906129
PAR ID:
10472345
Author(s) / Creator(s):
; ;
Publisher / Repository:
SIAM
Date Published:
Journal Name:
Multiscale Modeling & Simulation
Volume:
21
Issue:
4
ISSN:
1540-3459
Page Range / eLocation ID:
1344 to 1378
Subject(s) / Keyword(s):
momentum space, real space, 2D, electronic structure, density of states, conductivity, heterostructure, mechanical relaxation, moir\'e patterns
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Twisted bilayer graphene (TBG) has drawn significant interest due to recent experiments which show that TBG can exhibit strongly correlated behavior such as the superconducting and correlated insulator phases. Much of the theoretical work on TBG has been based on analysis of the Bistritzer–MacDonald model which includes a phenomenological parameter to account for lattice relaxation. In this work, we use a newly developed continuum model which systematically accounts for the effects of structural relaxation. In particular, we model structural relaxation by coupling linear elasticity to a stacking energy that penalizes disregistry. We compare the impact of the two relaxation models on the corresponding many-body model by defining an interacting model projected to the flat bands. We perform tests at charge neutrality at both the Hartree–Fock and coupled cluster singles and doubles level of theory and find the systematic relaxation model gives quantitative differences from the simplified relaxation model. 
    more » « less
  2. Abstract Magic angles in the chiral model of twisted bilayer graphene are parameters for which the chiral version of the Bistritzer–MacDonald Hamiltonian exhibits a flat band at energy zero. We compute the sums over powers of (complex) magic angles and use that to show that the set of magic angles is infinite. We also provide a new proof of the existence of the first real magic angle, showing also that the corresponding flat band has minimal multiplicity for the simplest possible choice of potentials satisfying all symmetries. These results indicate (though do not prove) a hidden integrability of the chiral model. 
    more » « less
  3. 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\'e 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. 
    more » « less
  4. Coupling between exciton states across the Brillouin zone in monolayer transition metal dichalcogenides can lead to ultrafast valley depolarization. Using time- and angle-resolved photoemission, we present momentum- and energy-resolved measurements of exciton coupling in monolayer WS2. By comparing full 4D (kx,ky,E,t) data sets after both linearly and circularly polarized excitation, we are able to disentangle intervalley and intravalley exciton coupling dynamics. Recording in the exciton binding energy basis instead of excitation energy, we observe strong mixing between the B1s exciton and An>1 states. The photoelectron energy and momentum distributions observed from excitons populated via intervalley coupling (e.g. K− → K+) indicate that the dominant valley depolarization mechanism conserves the exciton binding energy and center-of-mass momentum, consistent with intervalley Coulomb exchange. On longer timescales, exciton relaxation is accompanied by contraction of the momentum space distribution. 
    more » « less
  5. We calculate the weak-coupling resistivity of the two-dimensional Bose Hubbard model, comparing with the more familiar fermionic case. At high temperature the resistivity is linear in T , while in the low-temperature normal state it is exponentially suppressed. We explore the density dependence and calculate the momentum relaxation rate. 
    more » « less