No Arabic abstract
In this letter, we examine the role of Coulomb interactions in the emergence of macroscopically ordered states in graphene supported on hexagonal boron nitride substrates. Due to incommensuration effects with the substrate, graphene can develop gapped low energy modes that spatially conform into a triangular superlattice of quantum rings. In the presence of these modes, we show that Coulomb interactions lead to spontaneous formation of chiral loop currents in bulk and to macroscopic spin-valley order at zero temperature. We show that this exotic state breaks time reversal symmetry and can be detected with interferometry and polar Kerr measurements.
The effect of an hexagonal boron nitride (hBN) layer close aligned with twisted bilayer graphene (TBG) is studied. At sufficiently low angles between twisted bilayer graphene and hBN, $theta_{hBN} lesssim 2^circ$, the graphene electronic structure is strongly disturbed. The width of the low energy peak in the density of states changes from $W sim 5 - 10$ meV for a decoupled system to $sim 20 - 30$ meV. Spikes in the density of states due to van Hove singularities are smoothed out. We find that for a realistic combination of the twist angle in the TBG and the twist angle between the hBN and the graphene layer the system can be described using a single moire unit cell.
The Coulomb interaction is widely known to enhance the effective mass of interacting particles and therefore tends to favor a localized state at commensurate filling. Here, we will show that, in contrast to this consensus, in a van der Waals heterostructure consisting of graphene and hexagon boron nitride (h-BN), the onsite Coulomb repulsion will at first destroy the localized state. This is due to the fact that the onsite Coulomb repulsion tends to suppress the asymmetry between neighboring carbons induced by h-BN substrate. We corroborate this surprising phenomenon by solving a tight-binding model with onsite Coulomb repulsion treated within coherent potential approximation, where hopping parameters are derived from density functional theory calculations based on the graphene/h-BN heterostructure. Our results indicate that both gapless and gapped states observed experimentally in graphene/h-BN heterostructures can be understood after a realistic value of the onsite Coulomb repulsion as well as different interlayer distances are taken into account. Finally, we propose ways to enhance the gapped state which is essential for potential application of graphene to next-generation electronics. Furthermore, we argue that band gap suppressed by many-body effect should happen in other van der Waals heterostructures.
Graphene has demonstrated great promise for future electronics technology as well as fundamental physics applications because of its linear energy-momentum dispersion relations which cross at the Dirac point. However, accessing the physics of the low density region at the Dirac point has been difficult because of the presence of disorder which leaves the graphene with local microscopic electron and hole puddles, resulting in a finite density of carriers even at the charge neutrality point. Efforts have been made to reduce the disorder by suspending graphene, leading to fabrication challenges and delicate devices which make local spectroscopic measurements difficult. Recently, it has been shown that placing graphene on hexagonal boron nitride (hBN) yields improved device performance. In this letter, we use scanning tunneling microscopy to show that graphene conforms to hBN, as evidenced by the presence of Moire patterns in the topographic images. However, contrary to recent predictions, this conformation does not lead to a sizable band gap due to the misalignment of the lattices. Moreover, local spectroscopy measurements demonstrate that the electron-hole charge fluctuations are reduced by two orders of magnitude as compared to those on silicon oxide. This leads to charge fluctuations which are as small as in suspended graphene, opening up Dirac point physics to more diverse experiments than are possible on freestanding devices.
The recent observation [R. V. Gorbachev et al., Science {bf 346}, 448 (2014)] of nonlocal resistance $R_mathrm{NL}$ near the Dirac point (DP) of multiterminal graphene on aligned hexagonal boron nitride (G/hBN) has been interpreted as the consequence of topological valley Hall currents carried by the Fermi sea states just beneath the bulk gap $E_g$ induced by the inversion symmetry breaking. However, the valley Hall conductivity $sigma^v_{xy}$, quantized inside $E_g$, is not directly measurable. Conversely, the Landauer-B{u}ttiker formula, as numerically exact approach to observable nonlocal transport quantities, yields $R_mathrm{NL} equiv 0$ for the same simplistic Hamiltonian of gapped graphene that generates $sigma^v_{xy} eq 0$. We combine ab initio with quantum transport calculations to demonstrate that G/hBN wires with zigzag edges host dispersive edge states near the DP that are absent in theories based on the simplistic Hamiltonian. Although such edge states exist also in isolated zigzag graphene wires, aligned hBN is required to modify their energy-momentum dispersion and generate $R_mathrm{NL} eq 0$ near the DP persisting in the presence of edge disorder. Concurrently, the edge states resolve the long-standing puzzle of why the highly insulating state of G/hBN is rarely observed. We conclude that the observed $R_mathrm{NL}$ is unrelated to Fermi sea topological valley currents conjectured for gapped Dirac spectra.
When hexagonal boron nitride (hBN) and graphene are aligned at zero or small twist angle, a moire structure is formed due to the small lattice constant mismatch between the two structures. In this work, we analyze magnetic ordering tendencies, driven by onsite Coulomb interactions, of encapsulated bilayer graphene (BG) forming a moire structure with one (hBN-BG) or both hBN layers (hBN-BG-hBN), using the random phase approximation. The calculations are performed in a fully atomistic Hubbard model that takes into account all $pi$-electrons of the carbon atoms in one moire unit cell. We analyze the charge neutral case and find that the dominant magnetic ordering instability is uniformly antiferromagnetic. Furthermore, at low temperatures, the critical Hubbard interaction $U_c$ required to induce magnetic order is slightly larger in those systems where the moire structure has caused a band gap opening in the non-interacting picture, although the difference is less than 6%. Mean-field calculations are employed to estimate how such an interaction-induced magnetic order may change the observable single-particle gap sizes.