No Arabic abstract
Van der Waals (vdW) heterostructures ---formed by stacking or growing two-dimensional (2D) crystals on top of each other--- have emerged as a new promising route to tailor and engineer the properties of 2D materials. Twisted bilayer graphene (tBLG), a simple vdW structure where the interference between two misaligned graphene lattices leads to the formation of a moire pattern, is a test bed to study the effects of the interaction and misalignment between layers, key players for determining the electronic properties of these stackings. In this chapter, we present in a pedagogical way the general theory used to describe lattice mismatched and misaligned vdW structures. We apply it to the study of tBLG in the limit of small rotations and see how the coupling between the two layers leads both to an angle dependent renormalization of graphenes Fermi velocity and appearance of low-energy van Hove singularities. The optical response of this system is then addressed by computing the optical conductivity and the dispersion relation of tBLG surface plasmon-polaritons.
The generalized tight-binding model is developed to investigate the magneto-electronic properties in twisted bilayer graphene system. All the interlayer and intralayer atomic interactions are included in the Moire superlattice. The twisted bilayer graphene system is a zero-gap semiconductor with double-degenerate Dirac-cone structures, and saddle-point energy dispersions appearing at low energies for cases of small twisting angles. There exist rich and unique magnetic quantization phenomena, in which many Landau-level subgroups are induced due to specific Moire zone folding through modulating the various stacking angles. The Landau-level spectrum shows hybridized characteristics associated with the those in monolayer, and AA $&$ AB stackings. The complex relations among the different sublattices on the same and different graphene layers are explored in detail.
Experiments on bilayer graphene unveiled a fascinating realization of stacking disorder where triangular domains with well-defined Bernal stacking are delimited by a hexagonal network of strain solitons. Here we show by means of numerical simulations that this is a consequence of a structural transformation of the moir{e} pattern inherent of twisted bilayer graphene taking place at twist angles $theta$ below a crossover angle $theta^{star}=1.2^{circ}$. The transformation is governed by the interplay between the interlayer van der Waals interaction and the in-plane strain field, and is revealed by a change in the functional form of the twist energy density. This transformation unveils an electronic regime characteristic of vanishing twist angles in which the charge density converges, though not uniformly, to that of ideal bilayer graphene with Bernal stacking. On the other hand, the stacking domain boundaries form a distinct charge density pattern that provides the STM signature of the hexagonal solitonic network.
The electronic and optical response of Bernal stacked bilayer graphene with geometry modulation and gate voltage are studied. The broken symmetry in sublattices, one dimensional periodicity perpendicular to the domain wall and out-of-plane axis introduces substantial changes of wavefunctions, such as gapless topological protected states, standing waves with bonding and anti-bonding characteristics, rich structures in density of states and optical spectra. The wavefunctions present well-behaved standing waves in pure system and complicated node structures in geometry-modulated system. The optical absorption spectra show forbidden optical excitation channels, prominent asymmetric absorption peaks, and dramatic variations in absorption structures. These results provide that the geometry-modulated structure with tunable gate voltage could be used for electronic and optical manipulation in future graphene-based devices.
We study the electronic properties of twisted bilayers graphene in the tight-binding approximation. The interlayer hopping amplitude is modeled by a function, which depends not only on the distance between two carbon atoms, but also on the positions of neighboring atoms as well. Using the Lanczos algorithm for the numerical evaluation of eigenvalues of large sparse matrices, we calculate the bilayer single-electron spectrum for commensurate twist angles in the range $1^{circ}lesssimthetalesssim30^{circ}$. We show that at certain angles $theta$ greater than $theta_{c}approx1.89^{circ}$ the electronic spectrum acquires a finite gap, whose value could be as large as $80$ meV. However, in an infinitely large and perfectly clean sample the gap as a function of $theta$ behaves non-monotonously, demonstrating exponentially-large jumps for very small variations of $theta$. This sensitivity to the angle makes it impossible to predict the gap value for a given sample, since in experiment $theta$ is always known with certain error. To establish the connection with experiments, we demonstrate that for a system of finite size $tilde L$ the gap becomes a smooth function of the twist angle. If the sample is infinite, but disorder is present, we expect that the electron mean-free path plays the same role as $tilde L$. In the regime of small angles $theta<theta_c$, the system is a metal with a well-defined Fermi surface which is reduced to Fermi points for some values of $theta$. The density of states in the metallic phase varies smoothly with $theta$.
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.