No Arabic abstract
Multilayer graphene with rhombohedral and Bernal stacking are supposed to be metallic, as predicted by density functional theory calculations using semi-local functionals. However recent angular resolved photoemission and transport data have questioned this point of view. In particular, rhombohedral flakes are suggested to be magnetic insulators. Bernal flakes composed of an even number of layers are insulating, while those composed of an odd number of layers are pseudogapped. Here, by systematically benchmarking with plane waves codes, we develop very accurate all-electron Gaussian basis sets for graphene multilayers. We find that, in agreement with our previous calculations, rhombohedral stacked multilayer graphene are gapped for and magnetic. However, the valence band curvature and the details of the electronic structure depend crucially on the basis set. Only substantially extended basis sets are able to correctly reproduce the effective mass of the valence band top at the K point, while the popular POB-TZVP basis set leads to a severe overestimation. In the case of Bernal stacking, we show that exact exchange gaps the flakes composed by four layers and opens pseudogaps for N = 3, 6, 7, 8. However, the gap or pseudogap size and its behaviour as a function of thickness are not compatible with experimental data. Moreover, hybrid functionals lead to a metallic solution for 5 layers and a magnetic ground state for 5, 6 and 8 layers. Magnetism is very weak with practically no effect on the electronic structure and the magnetic moments are mostly concentrated in the central layers. Our hybrid functional calculations on trilayer Bernal graphene multilayers are in excellent agreement with non-magnetic GW calculations. For thicker multilayers, our calculations are a benchmark for manybody theoretical modeling of the low energy electronic structure.
We study the electronic structure of graphene with a single substitutional vacancy using a combination of the density-functional, tight-binding, and impurity Greens function approaches. Density functional studies are performed with the all-electron spin-polarized linear augmented plane wave (LAPW) method. The three $sp^2 sigma$ dangling bonds adjacent to the vacancy introduce localized states (V$sigma$) in the mid-gap region, which split due to the crystal field and a Jahn-Teller distortion, while the $p_z pi$ states introduce a sharp resonance state (V$pi$) in the band structure. For a planar structure, symmetry strictly forbids hybridization between the $sigma$ and the $pi$ states, so that these bands are clearly identifiable in the calculated band structure. As for the magnetic moment of the vacancy, the Hunds-rule coupling aligns the spins of the four localized V$sigma_1 uparrow downarrow$, V$sigma_2 uparrow $, and the V$pi uparrow$ electrons resulting in a S=1 state, with a magnetic moment of $2 mu_B$, which is reduced by about $0.3 mu_B$ due to the anti-ferromagnetic spin-polarization of the $pi$ band itinerant states in the vicinity of the vacancy. This results in the net magnetic moment of $1.7 mu_B$. Using the Lippmann-Schwinger equation, we reproduce the well-known $sim 1/r$ decay of the localized V$pi$ wave function with distance and in addition find an interference term coming from the two Dirac points, previously unnoticed in the literature. The long-range nature of the V$pi$ wave function is a unique feature of the graphene vacancy and we suggest that this may be one of the reasons for the widely varying relaxed structures and magnetic moments reported from the supercell band calculations in the literature.
In van der Waals bonded or rotationally disordered multilayer stacks of two-dimensional (2D) materials, the electronic states remain tightly confined within individual 2D layers. As a result, electron-phonon interactions occur primarily within layers and interlayer electrical conductivities are low. In addition, strong covalent in-plane intralayer bonding combined with weak van der Waals interlayer bonding results in weak phonon-mediated thermal coupling between the layers. We demonstrate here, however, that Coulomb interactions between electrons in different layers of multilayer epitaxial graphene provide an important mechanism for interlayer thermal transport even though all electronic states are strongly confined within individual 2D layers. This effect is manifested in the relaxation dynamics of hot carriers in ultrafast time-resolved terahertz spectroscopy. We develop a theory of interlayer Coulomb coupling containing no free parameters that accounts for the experimentally observed trends in hot-carrier dynamics as temperature and the number of layers is varied.
We have used scanning tunneling microscopy and spectroscopy to resolve the spatial variation of the density of states of twisted graphene layers on top of a highly oriented pyrolytic graphite substrate. Owing to the twist a moire pattern develops with a periodicity that is substantially larger than the periodicity of a single layer graphene. The twisted graphene layer has electronic properties that are distinctly different from that of a single layer graphene due to the nonzero interlayer coupling. For small twist angles (about 1-3.5 degree) the integrated differential conductivity spectrum exhibits two well-defined Van Hove singularities. Spatial maps of the differential conductivity that are recorded at energies near the Fermi level exhibit a honeycomb structure that is comprised of two inequivalent hexagonal sub-lattices. For energies |E-E_F|>0.3 eV the hexagonal structure in the differential conductivity maps vanishes. We have performed tight-binding calculations of the twisted graphene system using the propagation method, in which a third graphene layer is added to mimic the substrate. This third layer lowers the symmetry and explains the development of the two hexagonal sub-lattices in the moire pattern. Our experimental results are in excellent agreement with the tight-binding calculations.
Since its inception in 2001, the science and technology of epitaxial graphene on hexagonal silicon carbide has matured into a major international effort and is poised to become the first carbon electronics platform. A historical perspective is presented and the unique electronic properties of single and multilayered epitaxial graphenes on electronics grade silicon carbide are reviewed. Early results on transport and the field effect in Si-face grown graphene monolayers provided proof-of-principle demonstrations. Besides monolayer epitaxial graphene, attention is given to C-face grown multilayer graphene, which consists of electronically decoupled graphene sheets. Production, structure, and electronic structure are reviewed. The electronic properties, interrogated using a wide variety of surface, electrical and optical probes, are discussed. An overview is given of recent developments of several device prototypes including resistance standards based on epitaxial graphene quantum Hall devices and new ultrahigh frequency analog epitaxial graphene amplifiers.
We apply the semi-classical quantum Boltzmann formalism for the computation of transport properties to multilayer graphene. We compute the electrical conductivity as well as the thermal conductivity and thermopower for Bernal-stacked multilayers with an even number of layers. We show that the window for hydrodynamic transport in multilayer graphene is similar to the case of bilayer graphene. We introduce a simple hydrodynamic model which we dub the multi-fluid model and which can be used to reproduce the results for the electrical conductivity and thermopower from the quantum Boltzmann equation.