No Arabic abstract
The structural and electronic properties of twisted bilayer graphene are investigated from first principles and tight binding approach as a function of the twist angle (ranging from the first magic angle $theta=1.08^circ$ to $theta=3.89^circ$, with the former corresponding to the largest unit cell, comprising 11164 carbon atoms). By properly taking into account the long-range van der Waals interaction, we provide the patterns for the atomic displacements (with respect to the ideal twisted bilayer). The out-of-plane relaxation shows an oscillating (buckling) behavior, very evident for the smallest angles, with the atoms around the AA stacking regions interested by the largest displacements. The out-of-plane displacements are accompanied by a significant in-plane relaxation, showing a vortex-like pattern, where the vorticity (intended as curl of the displacement field) is reverted when moving from the top to the bottom plane and viceversa. Overall, the atomic relaxation results in the shrinking of the AA stacking regions in favor of the more energetically favorable AB/BA stacking domains. The measured flat bands emerging at the first magic angle can be accurately described only if the atomic relaxations are taken into account. Quite importantly, the experimental gaps separating the flat band manifold from the higher and lower energy bands cannot be reproduced if only in-plane or only out-of-plane relaxations are considered. The stability of the relaxed bilayer at the first magic angle is estimated to be of the order of 0.5-0.9 meV per atom (or 7-10 K). Our calculations shed light on the importance of an accurate description of the vdW interaction and of the resulting atomic relaxation to envisage the electronic structure of this really peculiar kind of vdW bilayers.
We present a simple model that we believe captures the key aspects of the competition between superconducting and insulating states in twisted bilayer graphene. Within this model, the superconducting phase is primary, and arises at generic fillings, but is interrupted by the insulator at commensurate fillings. Importantly, the insulator forms because of electron-electron interactions, but the model is agnostic as to the superconducting pairing mechanism, which need not originate with electron-electron interactions. The model is composed of a collection of crossed one-dimensional quantum wires whose intersections form a superlattice. At each superlattice point, we place a locally superconducting puddle which can exchange Cooper pairs with the quantum wires. We analyze this model assuming weak wire-puddle and wire-wire couplings. We show that for a range of repulsive intrawire interactions, the system is superconducting at `generic incommensurate fillings, with the superconductivity being `interrupted by an insulating phase at commensurate fillings. We further show that the gapped insulating states at commensurate fillings give way to gapless states upon application of external Zeeman fields. These features are consistent with experimental observations in magic-angle twisted bilayer graphenes despite the distinct microscopic details. We further study the full phase diagram of this model and discover that it contains several distinct correlated insulating states, which we characterize herein.
We introduce and analyze a model that sheds light on the interplay between correlated insulating states, superconductivity, and flavor-symmetry breaking in magic angle twisted bilayer graphene. Using a variational mean-field theory, we determine the normal-state phase diagram of our model as a function of the band filling. The model features robust insulators at even integer fillings, occasional weaker insulators at odd integer fillings, and a pattern of flavor-symmetry breaking at non-integer fillings. Adding a phonon-mediated inter-valley retarded attractive interaction, we obtain strong-coupling superconducting domes, whose structure is in qualitative agreement with experiments. Our model elucidates how the intricate form of the interactions and the particle-hole asymmetry of the electronic structure determine the phase diagram. It also explains how subtle differences between devices may lead to the different behaviors observed experimentally. A similar model can be applied with minor modifications to other moir{e} systems, such as twisted trilayer graphene.
We investigate the topological properties of Floquet-engineered twisted bilayer graphene above the magic angle driven by circularly polarized laser pulses. Employing a full Moire-unit-cell tight-binding Hamiltonian based on first-principles electronic structure we show that the band topology in the bilayer, at twisting angles above 1.05$^circ$, essentially corresponds to the one of single-layer graphene. However, the ability to open topologically trivial gaps in this system by a bias voltage between the layers enables the full topological phase diagram to be explored, which is not possible in single-layer graphene. Circularly polarized light induces a transition to a topologically nontrivial Floquet band structure with the Berry curvature of a Chern insulator. Importantly, the twisting allows for tuning electronic energy scales, which implies that the electronic bandwidth can be tailored to match realistic driving frequencies in the ultraviolet or mid-infrared photon-energy regimes. This implies that Moire superlattices are an ideal playground for combining twistronics, Floquet engineering, and strongly interacting regimes out of thermal equilibrium.
Twisted double bilayer graphene has recently emerged as an interesting moire material that exhibits strong correlation phenomena that are tunable by an applied electric field. Here we study the atomic and electronic properties of three different graphene double bilayers: double bilayers composed of two AB stacked bilayers (AB/AB), double bilayers composed of two AA stacked bilayers (AA/AA) as well as heterosystems composed of one AB and one AA bilayer (AB/AA). The atomic structure is determined using classical force fields. We find that the inner layers of the double bilayer exhibit significant in-plane and out-of-plane relaxations, similar to twisted bilayer graphene. The relaxations of the outer layers depend on the stacking: atoms in AB bilayers follow the relaxations of the inner layers, while atoms in AA bilayers attempt to avoid higher-energy AA stacking. For the relaxed structures, we calculate the electronic band structures using the tight-binding method. All double bilayers exhibit flat bands at small twist angles, but the shape of the bands depends sensitively on the stacking of the outer layers. To gain further insight, we study the evolution of the band structure as the outer layers are rigidly moved away from the inner layers, while preserving their atomic relaxations. This reveals that the hybridization with the outer layers results in an additional flattening of the inner-layer flat band manifold. Our results establish AA/AA and AB/AA twisted double bilayers as interesting moire materials with different flat band physics compared to the widely studied AB/AB system.
We calculate the interactions between the Wannier functions of the 8-orbital model for twisted bilayer graphene (TBG). In this model, two orbitals per valley centered at the AA regions, the AA-p orbitals, account for the most part of the spectral weight of the flats bands. Exchange and assisted-hopping terms between these orbitals are found to be small. Therefore, the low energy properties of TBG will be determined by the density-density interactions. These interactions decay with the distance much faster than in the two orbital model, following a 1/r law in the absence of gates. The magnitude of the largest interaction in the model, the onsite term between the flat band orbitals, is controlled by the size of the AA regions and is estimated to be ~ 40 meV. To screen this interaction, the metallic gates have to be placed at a distance smaller than 5 nm. For larger distances only the long-range part of the interaction is substantially screened. The model reproduces the band deformation induced by doping found in other approaches within the Hartree approximation. Such deformation reveals the presence of other orbitals in the flat bands and is sensitive to the inclusion of the interactions involving them.