Graphene bilayer with a twist: electronic structure

Scanning probe analysis of twisted graphene grown on a graphene/silicon carbide template

Overlayer growth of graphene on an epitaxial graphene/silicon carbide (SiC) as a solid template by ethanol chemical vapor deposition is performed over a wide growth temperature range from 900 °C to 1450 °C. Structural analysis using atomic force and scanning tunneling microscopies reveal that graphene islands grown at 1300 °C form hexagonal twisted bilayer graphene as a single crystal. When the growth temperature exceeds 1400 °C, the grown graphene islands show a circular shape. Moreover, moiré patterns with different periods are observed in a single graphene island. This means that the graphene islands grown at high temperature are composed of several graphene domains with different twist angles. From these results, we conclude that graphene overlayer growth on the epitaxial graphene/SiC solid at 1300 °C effectively synthesizes the twisted few-layer graphene with a high crystallinity.

Dirac Magic and Lifshitz Transitions in AA-Stacked Twisted Multilayer Graphene

We uncover a new type of magic-angle phenomena when an AA-stacked graphene bilayer is twisted relative to another graphene system with band touching. In the simplest case this constitutes a trilayer system formed by an AA-stacked bilayer twisted relative to a single layer of graphene. We find multiple anisotropic Dirac cones coexisting in such twisted multilayer structures at certain angles, which we call "Dirac magic." We trace the origin of Dirac magic angles to the geometric structure of the twisted AA-bilayer Dirac cones relative to the other band-touching spectrum in the moiré reciprocal lattice. The anisotropy of the Dirac cones and a concomitant cascade of saddle points induce a series of topological Lifshitz transitions that can be tuned by the twist angle and perpendicular electric field. We discuss the possibility of direct observation of Dirac magic as well as its consequences for the correlated states of electrons in this moiré system.

First-principles study on the heterostructure of twisted graphene/hexagonal boron nitride/graphene sandwich structure

In recent years, the discovery of 'magic angle' graphene has given new inspiration to the formation of heterojunctions. Similarly, the use of hexagonal boron nitride, known as white graphene, as a substrate for graphene devices has more aroused great interest in the graphene/hexagonal boron nitride heterostructure system. Based on the first principles method of density functional theory, the band structure, density of states, Mulliken population, and differential charge density of a tightly packed model of twisted graphene/hexagonal boron nitride/graphene sandwich structure have been studied. Through the establishment of heterostructure models twisted bilayer-graphene inserting hBN with different twisted angles, it was found that the band gap, Mulliken population, and charge density, exhibited specific evolution regulars with the rotation angle of the upper graphene, showing novel electronic properties and realizing metal-insulator phase transition. We find that the particular value of the twist angle at which the metal-insulator phase transition occurs and propose a rotational regulation mechanism with angular periodicity. Our results have guiding significance for the practical application of heterojunction electronic devices.

Correlations in twisted double-layer graphene with virtual photons in a microcavity

We analyze the entanglement generation of a system composed of two decoupled rotated graphene layers inside a planar microcavity. By considering the electromagnetic field of the cavity in the vacuum state and using time-dependent perturbation theory it is possible to obtain the range of geometric parameters at which the quantum states of electrons in different layers are entangled. By employing the negativity measure, correlations between layers are obtained for time scales smaller than the light-crossing time of the layers. It is shown that the negativity measure is modulated by the rotation angle between layers, allowing manipulation ofXstates. Finally, an experimental protocol is analyzed in order to detect non-causal effects between layers, by allowing back-voltage switching functions in the two layers with supports that do not overlap in time. By turning off the second-back voltage at a time smaller than the light-crossing time, it is possible to obtain correlations between layers through the independent interaction with virtual photons. The exchange of virtual photons implies that the propagator can be nonzero outside the light cone and this non-causal propagation can create entangled quantum states.

Interlayer Interactions in 1D Van der Waals Moiré Superlattices

Studying two-dimensional (2D) van der Waals (vdW) moiré superlattices and their interlayer interactions have received surging attention after recent discoveries of many new phases of matter that are highly tunable. Different atomistic registry between layers forming the inner and outer nanotubes can also form one-dimensional (1D) vdW moiré superlattices. In this review, experimental observations and theoretical perspectives related to interlayer interactions in 1D vdW moiré superlattices are summarized. The discussion focuses on double-walled carbon nanotubes (DWNTs), a model 1D vdW moiré system, and the authors highlight the new optical features emerging from the non-trivial strong interlayer coupling effect and the unique physics in 1D DWNTs. Future directions and questions in probing the intriguing physical phenomena in 1D vdW moiré superlattices such as, correlated physics in different 1D moiré systems beyond DWNTs are proposed and discussed.

Science of 2.5 dimensional materials: paradigm shift of materials science toward future social innovation

The past decades of materials science discoveries are the basis of our present society - from the foundation of semiconductor devices to the recent development of internet of things (IoT) technologies. These materials science developments have depended mainly on control of rigid chemical bonds, such as covalent and ionic bonds, in organic molecules and polymers, inorganic crystals and thin films. The recent discovery of graphene and other two-dimensional (2D) materials offers a novel approach to synthesizing materials by controlling their weak out-of-plane van der Waals (vdW) interactions. Artificial stacks of different types of 2D materials are a novel concept in materials synthesis, with the stacks not limited by rigid chemical bonds nor by lattice constants. This offers plenty of opportunities to explore new physics, chemistry, and engineering. An often-overlooked characteristic of vdW stacks is the well-defined 2D nanospace between the layers, which provides unique physical phenomena and a rich field for synthesis of novel materials. Applying the science of intercalation compounds to 2D materials provides new insights and expectations about the use of the vdW nanospace. We call this nascent field of science '2.5 dimensional (2.5D) materials,' to acknowledge the important extra degree of freedom beyond 2D materials. 2.5D materials not only offer a new field of scientific research, but also contribute to the development of practical applications, and will lead to future social innovation. In this paper, we introduce the new scientific concept of this science of '2.5D materials' and review recent research developments based on this new scientific concept.

Emergence and Tuning of Multiple Flat Bands in Twisted Bilayer γ-Graphyne

Using twisted bilayer γ-graphyne (TBGY) as an example, we show that it is possible to create multiple flat bands in carbon allotropes without the requirement of a specified magic angle. The origin of the flat bands can be understood by a simple two-level coupling model. The narrow bandwidth and strong localization of the flat band states might lead to strong correlation effects, which make TBGY a good platform for studying correlation physics. On the basis of the two-level coupling model, we further propose that the width and extent of localization of flat bands can be tuned by an energy mismatch ΔE between the two layers of TBGY, which can be realized by either applying a perpendicular electric field or introducing a heterostrain. This allows continuous modulation of TBGY from the strong-correlation regime to the medium- or weak-correlation regime, which could be utilized to study the quantum phase transition.

Oscillations of the Spacing between van Hove Singularities Induced by sub-Ångstrom Fluctuations of Interlayer Spacing in Graphene Superlattices

Physical properties of two-dimensional van der Waals (vdWs) structures depend sensitively on both stacking orders and interlayer interactions. Yet, in most cases studied to date, the interlayer interaction is considered to be a "static" property of the vdWs structures. Here we demonstrate that applying a scanning tunneling microscopy (STM) tip pulse on twisted bilayer graphene (TBG) can induce sub-Ångstrom fluctuations of the interlayer separation in the TBG, which are equivalent to dynamic vertical external pressure of about 10 GPa on the TBG. The sub-Ångstrom fluctuations of the interlayer separation result in large oscillations of the energy separations between two van Hove singularities (VHSs) in the TBG. The period of the oscillations of the VHSs spacing is extremely long, about 500-1000 sec, attributing to tip-induced local stress in the atomic-thick TBG. Our result provides an efficient method to tune and measure the physical properties of the vdWs structures dynamically.

Emerging properties of carbon based 2D material beyond graphene

Graphene turns out to be the pioneering material for setting up boulevard to a new zoo of recently proposed carbon based novel two dimensional (2D) analogues. It is evident that their electronic, optical and other related properties are utterly different from that of graphene because of the distinct intriguing morphology. For instance, the revolutionary emergence of Dirac cones in graphene is particularly hard to find in most of the other 2D materials. As a consequence the crystal symmetries indeed act as a major role for predicting electronic band structure. Since tight binding calculations have become an indispensable tool in electronic band structure calculation, we indicate the implication of such method in graphene's allotropes beyond hexagonal symmetry. It is to be noted that some of these graphene allotropes successfully overcome the inherent drawback of the zero band gap nature of graphene. As a result, these 2D nanomaterials exhibit great potential in a broad spectrum of applications, viz nanoelectronics, nanooptics, gas sensors, gas storages, catalysis, and other specific applications. The miniaturization of high performance graphene allotrope based gas sensors to microscopic or even nanosized range has also been critically discussed. In addition, various optical properties like the dielectric functions, optical conductivity, electron energy loss spectra reveal that these systems can be used in opto-electronic devices. Nonetheless, the honeycomb lattice of graphene is not superconducting. However, it is proposed that the tetragonal form of graphene can be intruded to form new hybrid 2D materials to achieve novel superconducting device at attainable conditions. These dynamic experimental prospects demand further functionalization of these systems to enhance the efficiency and the field of multifunctionality. This topical review aims to highlight the latest advances in carbon based 2D materials beyond graphene from the basic theoretical as well as future application perspectives.

Dichotomy of Electron-Phonon Coupling in Graphene Moiré Flat Bands

Graphene moiré superlattices are outstanding platforms to study correlated electron physics and superconductivity with exceptional tunability. However, robust superconductivity has been measured only in magic-angle twisted bilayer graphene (MA-TBG) and magic-angle twisted trilayer graphene (MA-TTG). The absence of a superconducting phase in certain moiré flat bands raises a question on the superconducting mechanism. In this work, we investigate electronic structure and electron-phonon coupling in graphene moiré superlattices based on atomistic calculations. We show that electron-phonon coupling strength λ is dramatically different among graphene moiré flat bands. The total strength λ is very large (λ>1) for MA-TBG and MA-TTG, both of which display robust superconductivity in experiments. However, λ is an order of magnitude smaller in twisted double bilayer graphene (TDBG) and twisted monolayer-bilayer graphene (TMBG) where superconductivity is reportedly rather weak or absent. We find that the Bernal-stacked layers in TDBG and TMBG induce sublattice polarization in the flat-band states, suppressing intersublattice electron-phonon matrix elements. We also obtain the nonadiabatic superconducting transition temperature T_{c} that matches well with the experimental results. Our results clearly show a correlation between strong electron-phonon coupling and experimental observations of robust superconductivity.

Exact Diagonalization for Magic-Angle Twisted Bilayer Graphene

We report on finite-size exact-diagonalization calculations in a Hilbert space defined by the continuum-model flat moiré bands of magic angle twisted bilayer graphene. For moiré band filling 3>|ν|>2, where superconductivity is strongest, we obtain evidence that the ground state is a spin ferromagnet. Near |ν|=3, we find Chern insulator ground states that have spontaneous spin, valley, and sublattice polarization, and demonstrate that the anisotropy energy in this order-parameter space is strongly band-filling-factor dependent. We emphasize that inclusion of the remote band self-energy is necessary for a reliable description of magic angle twisted bilayer graphene flat band correlations.

Magic-angle lasers in nanostructured moiré superlattice

Conventional laser cavities require discontinuity of material property or disorder to localize a light field for feedback. Recently, an emerging class of materials, twisted van der Waals materials, have been explored for applications in electronics and photonics. Here we propose and develop magic-angle lasers, where the localization is realized in periodic twisted photonic graphene superlattices. We reveal that the confinement mechanism of magic-angle lasers does not rely on a full bandgap but on the mode coupling between two twisted layers of photonic graphene lattice. Without any fine-tuning in structure parameters, a simple twist can result in nanocavities with strong field confinement and a high quality factor. Furthermore, the emissions of magic-angle lasers allow direct imaging of the wavefunctions of magic-angle states. Our work provides a robust platform to construct high-quality nanocavities for nanolasers, nano light-emitting diodes, nonlinear optics and cavity quantum electrodynamics at the nanoscale.

Carrier transport theory for twisted bilayer graphene in the metallic regime

Understanding the normal-metal state transport in twisted bilayer graphene near magic angle is of fundamental importance as it provides insights into the mechanisms responsible for the observed strongly correlated insulating and superconducting phases. Here we provide a rigorous theory for phonon-dominated transport in twisted bilayer graphene describing its unusual signatures in the resistivity (including the variation with electron density, temperature, and twist angle) showing good quantitative agreement with recent experiments. We contrast this with the alternative Planckian dissipation mechanism that we show is incompatible with available experimental data. An accurate treatment of the electron-phonon scattering requires us to go well beyond the usual treatment, including both intraband and interband processes, considering the finite-temperature dynamical screening of the electron-phonon matrix element, and going beyond the linear Dirac dispersion. In addition to explaining the observations in currently available experimental data, we make concrete predictions that can be tested in ongoing experiments.

Purcell-like Enhancement of Electron-Phonon Interactions in Long-Period Superlattices: Linear-Temperature Resistivity and Cooling Power

In the Purcell effect, the efficiency of optical emitters is enhanced by reducing the optical mode volume. Here we predict an analogous enhancement for electron-phonon (el-ph) scattering, achieved by compressing the electronic Wannier orbitals. Reshaped Wannier orbitals are a prominent attribute of graphene moiré superlattices, where the orbital size is tunable by the twist angle. A reduction in the orbital size leads to an enhancement in the el-ph interaction strength, yielding the values considerably larger than those in pristine monolayer graphene. The enhanced coupling boosts the el-ph scattering rates, pushing them above the values expected for the flat-band-enhanced density of electronic states. The enhanced phonon emission and scattering rates are manifested through the observables such as the electron-lattice cooling and the linear-temperature (T) resistivity, both of which are directly tunable by the moiré twist angle.

Heterostrain Determines Flat Bands in Magic-Angle Twisted Graphene Layers

The moiré of twisted graphene bilayers can generate flat bands in which charge carriers do not possess enough kinetic energy to escape Coulomb interactions with each other, leading to the formation of novel strongly correlated electronic states. This exceptionally rich physics relies on the precise arrangement between the layers. Here, we survey published scanning tunneling microscope measurements to prove that near the magic-angle, native heterostrain, the relative deformations between the layers, dominates twist in determining the flat bands as opposed to the common belief. This is demonstrated at full filling where electronic correlations have a weak effect and where we also show that tip-induced strain can have a strong influence. In the opposite situation of zero doping, we find that electronic correlation further renormalizes the flat bands in a way that strongly depends on experimental details.

Moiré, Euler and self-similarity - the lattice parameters of twisted hexagonal crystals

A real-space approach for the calculation of the moiré lattice parameters for superstructures formed by a set of rotated hexagonal 2D crystals such as graphene or transition-metal dichalcogenides is presented. Apparent moiré lattices continuously form for all rotation angles, and their lattice parameter to a good approximation follows a hyperbolical angle dependence. Moiré crystals, i.e. moiré lattices decorated with a basis, require more crucial assessment of the commensurabilities and lead to discrete solutions and a non-continuous angle dependence of the moiré-crystal lattice parameter. In particular, this lattice parameter critically depends on the rotation angle, and continuous variation of the angle can lead to apparently erratic changes of the lattice parameter. The solutions form a highly complex pattern, which reflects number-theoretical relations between formation parameters of the moiré crystal. The analysis also provides insight into the special case of a 30° rotation of the constituting lattices, for which a dodecagonal quasicrystalline structure forms.

Resonance Raman enhancement by the intralayer and interlayer electron-phonon processes in twisted bilayer graphene

Twisted bilayer graphene is a fascinating system due to the possibility of tuning the electronic and optical properties by controlling the twisting angle \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ between the layers. The coupling between the Dirac cones of the two graphene layers gives rise to van Hove singularities (vHs) in the density of electronic states, whose energies vary with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ. Raman spectroscopy is a fundamental tool to study twisted bilayer graphene (TBG) systems since the Raman response is hugely enhanced when the photons are in resonance with transition between vHs and new peaks appear in the Raman spectra due to phonons within the interior of the Brillouin zone of graphene that are activated by the Moiré superlattice. It was recently shown that these new peaks can be activated by the intralayer and the interlayer electron–phonon processes. In this work we study how each one of these processes enhances the intensities of the peaks coming from the acoustic and optical phonon branches of graphene. Resonance Raman measurements, performed in many different TBG samples with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\theta$$\end{document}θ between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4^{\circ }$$\end{document}4∘ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$16^{\circ }$$\end{document}16∘ and using several different laser excitation energies in the near-infrared (NIR) and visible ranges (1.39–2.71 eV), reveal the distinct enhancement of the different phonons of graphene by the intralayer and interlayer processes. Experimental results are nicely explained by theoretical calculations of the double-resonance Raman intensity in graphene by imposing the momentum conservation rules for the intralayer and the interlayer electron–phonon resonant conditions in TBGs. Our results show that the resonant enhancement of the Raman response in all cases is affected by the quantum interference effect and the symmetry requirements of the double resonance Raman process in graphene.

Twist dependent magneto-optical response in twisted bilayer graphene

By employing a linearised Boltzmann equation, we calculate the magneto-optical properties of twisted bilayer graphene using non-magnetic wave functions. Both transverse and longitudinal responses are calculated up to the second order in applied magnetic field with their twist angle and Fermi level dependence examined. We find that increasing the twist angle increases the transverse metallic response so long as the Fermi level remains below the upper conduction band. Interlayer transitions provide an appreciable enhancement when the Fermi level traverses the gap between the two conduction bands. Interlayer transitions are also responsible for a nonzero anomalous Hall conductivity in this model. As the Fermi level moves towards zero, the longitudinal response begins to dominate and a highly anisotropic negative magneto-resistance is observed.

Interplay between structural deformations and flat band phenomenology in twisted bilayer antimonene

In this work we apply first principles calculations to investigate the flat band phenomenology in twisted antimonene bilayer. We show that the relatively strong interlayer interactions which characterize this compound have profound effects in the emergence and properties of the flat bands. Specifically, when the moiré length becomes large enough to create well defined stacking patterns along the structure, out-of-plane displacements take place and are stabilized in the regions dominated by the AB stacking, leading to the emergence of flat bands. The interplay between structural and electronic properties allows for detection of flat bands in higher twist angles comparable to other two-dimensional materials. We also show that their energy position may be modulated by noncovalent functionalization with electron acceptor molecules.

Out-of-Plane Dielectric Susceptibility of Graphene in Twistronic and Bernal Bilayers

We describe how the out-of-plane dielectric polarizability of monolayer graphene influences the electrostatics of bilayer graphene-both Bernal (BLG) and twisted (tBLG). We compare the polarizability value computed using density functional theory with the output from previously published experimental data on the electrostatically controlled interlayer asymmetry potential in BLG and data on the on-layer density distribution in tBLG. We show that monolayers in tBLG are described well by polarizability α_exp = 10.8 ų and effective out-of-plane dielectric susceptibility ϵ_z = 2.5, including their on-layer electron density distribution at zero magnetic field and the interlayer Landau level pinning at quantizing magnetic fields.

Coulomb interaction, phonons, and superconductivity in twisted bilayer graphene

The polarizability of twisted bilayer graphene, due to the combined effect of electron-hole pairs, plasmons, and acoustic phonons, is analyzed. The screened Coulomb interaction allows for the formation of Cooper pairs and superconductivity in a significant range of twist angles and fillings. The tendency toward superconductivity is enhanced by the coupling between longitudinal phonons and electron-hole pairs. Scattering processes involving large momentum transfers, Umklapp processes, play a crucial role in the formation of Cooper pairs. The magnitude of the superconducting gap changes among the different pockets of the Fermi surface.

Strong Interminivalley Scattering in Twisted Bilayer Graphene Revealed by High-Temperature Magneto-Oscillations

Twisted bilayer graphene (TBG) provides an example of a system in which the interplay of interlayer interactions and superlattice structure impacts electron transport in a variety of nontrivial ways and gives rise to a plethora of interesting effects. Understanding the mechanisms of electron scattering in TBG has, however, proven challenging, raising many questions about the origins of resistivity in this system. Here we show that TBG exhibits high-temperature magneto-oscillations originating from the scattering of charge carriers between TBG minivalleys. The amplitude of these oscillations reveals that interminivalley scattering is strong, and its characteristic timescale is comparable to that of its intraminivalley counterpart. Furthermore, by exploring the temperature dependence of these oscillations, we estimate the electron-electron collision rate in TBG and find that it exceeds that of monolayer graphene. Our study demonstrates the consequences of the relatively small size of the superlattice Brillouin zone and Fermi velocity reduction on lateral transport in TBG.

Emulating Heavy Fermions in Twisted Trilayer Graphene

Twisted van der Waals materials have been shown to host a variety of tunable electronic structures. Here we put forward twisted trilayer graphene (TTG) as a platform to emulate heavy fermion physics. We demonstrate that TTG hosts extended and localized modes with an electronic structure that can be controlled by interlayer bias. In the presence of interactions, the existence of localized modes leads to the development of local moments, which are Kondo coupled to coexisting extended states. By electrically controlling the effective exchange between local moments, the system can be driven from a magnetic into a heavy fermion regime, passing through a quantum critical point, allowing one to electrically explore a generalized Doniach phase diagram. Our results put forward twisted graphene multilayers as a platform for the realization of strongly correlated heavy fermion physics in a purely carbon-based platform.

Flat band carrier confinement in magic-angle twisted bilayer graphene

Magic-angle twisted bilayer graphene has emerged as a powerful platform for studying strongly correlated electron physics, owing to its almost dispersionless low-energy bands and the ability to tune the band filling by electrostatic gating. Techniques to control the twist angle between graphene layers have led to rapid experimental progress but improving sample quality is essential for separating the delicate correlated electron physics from disorder effects. Owing to the 2D nature of the system and the relatively low carrier density, the samples are highly susceptible to small doping inhomogeneity which can drastically modify the local potential landscape. This potential disorder is distinct from the twist angle variation which has been studied elsewhere. Here, by using low temperature scanning tunneling spectroscopy and planar tunneling junction measurements, we demonstrate that flat bands in twisted bilayer graphene can amplify small doping inhomogeneity that surprisingly leads to carrier confinement, which in graphene could previously only be realized in the presence of a strong magnetic field.

The role of the 2D-to-3D transition in x-ray diffraction analysis of crystallite size

The diffraction behaviour of stacked layers of graphene and hexagonal boron nitride are studied computationally by direct calculation of the diffraction pattern using the Debye scattering equation. Analysis of the position and profile of the diffraction peaks show that while single-layer graphene is unambiguously a 2D material, ordered stacks of three or more layers diffract as bulk material. Following the Scherrer equation, we correlate the known crystallite size with the diffraction peak parameters and observe strong affine relationships which exist separately for single-layer, bi-layer and three or more layer (bulk) structures. We determine a series of expressions to calculate the crystallite size which do not suffer the well-known size-dependence or rely on assumptions about the shape. We present a detailed workflow showing how these expressions can be applied to experimental data.

Flat Bands in Magic-Angle Bilayer Photonic Crystals at Small Twists

The new physics of magic-angle twisted bilayer graphene (TBG) motivated extensive studies of flat bands hosted by moiré superlattices in van der Waals structures, inspiring the investigations into their photonic counterparts with potential applications including Bose-Einstein condensation. However, correlation between photonic flat bands and bilayer photonic moiré systems remains unexplored, impeding further development of moiré photonics. In this work, we formulate a coupled-mode theory for low-angle twisted bilayer honeycomb photonic crystals as a close analogy of TBG, discovering magic-angle photonic flat bands with a non-Anderson-type localization. Moreover, the interlayer separation constitutes a convenient degree of freedom in tuning photonic moiré bands without high pressure. A phase diagram is constructed to correlate the twist angle and separation dependencies to the photonic magic angles. Our findings reveal a salient correspondence between fermionic and bosonic moiré systems and pave the avenue toward novel applications through advanced photonic band or state engineering.

Trigonal quasicrystalline states in [Formula: see text] rotated double moiré superlattices

We study the lattice configuration and electronic structure of a double moiré superlattice, which is composed of a graphene layer encapsulated by two other layers in a way such that the two hexagonal moiré patterns are arranged in a dodecagonal quasicrystalline configuration. We show that there are between 0 and 4 such configurations depending on the lattice mismatch between graphene and the encapsulating layer. We then reveal the resonant interaction, which is distinct from the conventional 2-, 3-, 4-wave mixing of moiré superlattices, that brings together and hybridizes twelve degenerate Bloch states of monolayer graphene. These states do not fully satisfy the dodecagonal quasicrystalline rotational symmetry due to the symmetry of the wave vectors involved. Instead, their wave functions exhibit trigonal quasicrystalline order, which lacks inversion symmetry, at the energies much closer to the charge neutrality point of graphene.

Emerging flat bands in large-angle twisted bi-layer graphene under pressure

Recent experiments on magic-angle twisted bi-layer graphene have attracted intensive attention due to exotic properties such as unconventional superconductivity and correlated insulation. These phenomena were often found at a magic angle less than 1.1°. However, the preparation of precisely controlled bi-layer graphene with a small magic angle is challenging. In this work, electronic properties of large-angle twisted bi-layer graphene (TBG) under pressure are investigated with density functional theory. We demonstrate that large-angle TBG can display flat bands nearby the Fermi level under pressure, which may also induce interesting properties such as superconductivity which have only been found in small-angle TBG at ambient pressure. The Fermi velocity is found to decrease monotonously with pressure for large twisted angles, e.g., 21.8°. Our work indicates that applying pressure provides opportunities for flat-band engineering in larger angle TBG and supports further exploration in related investigations.

Electronic Structures of Twisted Bilayer InSe/InSe and Heterobilayer Graphene/InSe

Building vertical van der Waals heterojunctions between two-dimensional layered materials has become a promising strategy for modulating the properties of two-dimensional materials. Herein, we investigate the electronic structures of non-twisted/twisted bilayer InSe/InSe and heterobilayer graphene/InSe (Gr/InSe) by employing density functional theory calculations. For twisted bilayer InSe/InSes, their interlayer distances and band gaps are almost identical but a bit larger than those of the AB-stacking one due to the spontaneous polarization. Differently, the band gaps of twisted Gr/InSe are found to vary with the rotation angles. Our results provide an effective way to tune the electronic properties of two-dimensional materials.

Quantifying the Charge Carrier Interaction in Metallic Twisted Bilayer Graphene Superlattices

The mechanism of charge carrier interaction in twisted bilayer graphene (TBG) remains an unresolved problem, where some researchers proposed the dominance of the electron-phonon interaction, while the others showed evidence for electron-electron or electron-magnon interactions. Here we propose to resolve this problem by generalizing the Bloch-Grüneisen equation and using it for the analysis of the temperature dependent resistivity in TBG. It is a well-established theoretical result that the Bloch-Grüneisen equation power-law exponent, p, exhibits exact integer values for certain mechanisms. For instance, p = 5 implies the electron-phonon interaction, p = 3 is associated with the electron-magnon interaction and p = 2 applies to the electron-electron interaction. Here we interpret the linear temperature-dependent resistance, widely observed in TBG, as p→1, which implies the quasielastic charge interaction with acoustic phonons. Thus, we fitted TBG resistance curves to the Bloch-Grüneisen equation, where we propose that p is a free-fitting parameter. We found that TBGs have a smoothly varied p-value (ranging from 1.4 to 4.4) depending on the Moiré superlattice constant, λ, or the charge carrier concentration, n. This implies that different mechanisms of the charge carrier interaction in TBG superlattices smoothly transition from one mechanism to another depending on, at least, λ and n. The proposed generalized Bloch-Grüneisen equation is applicable to a wide range of disciplines, including superconductivity and geology.

Recent Advances in 2D Superconductors

The emergence of superconductivity in 2D materials has attracted much attention and there has been rapid development in recent years because of their fruitful physical properties, such as high transition temperature (T_c ), continuous phase transition, and enhanced parallel critical magnetic field (B_c ). Tremendous efforts have been devoted to exploring different physical parameters to figure out the mechanisms behind the unexpected superconductivity phenomena, including adjusting the thickness of samples, fabricating various heterostructures, tuning the carrier density by electric field and chemical doping, and so on. Here, different types of 2D superconductivity with their unique characteristics are introduced, including the conventional Bardeen-Cooper-Schrieffer superconductivity in ultrathin films, high-T_c superconductivity in Fe-based and Cu-based 2D superconductors, unconventional superconductivity in newly discovered twist-angle bilayer graphene, superconductivity with enhanced B_c , and topological superconductivity. A perspective toward this field is then proposed based on academic knowledge from the recently reported literature. The aim is to provide researchers with a clear and comprehensive understanding about the newly developed 2D superconductivity and promote the development of this field much further.

Theoretical Prediction of Two-Dimensional Materials, Behavior, and Properties

Predictive modeling of two-dimensional (2D) materials is at the crossroad of two current rapidly growing interests: 2D materials per se, massively sought after and explored in experimental laboratories, and materials theoretical-computational models in general, flourishing on a fertile mix of condensed-matter physics and chemistry with advancing computational technology. Here the general methods and specific techniques of modeling are briefly overviewed, along with a somewhat philosophical assessment of what "prediction" is, followed by selected practical examples for 2D materials, from structures and properties, to device functionalities and synthetic routes for their making. We conclude with a brief sketch-outlook of future developments.

Moiré Patterns in 2D Materials: A Review

Quantum materials have attracted much attention in recent years due to their exotic and incredible properties. Among them, van der Waals materials stand out due to their weak interlayer coupling, providing easy access to manipulating electrical and optical properties. Many fascinating electrical, optical, and magnetic properties have been reported in the moiré superlattices, such as unconventional superconductivity, photonic dispersion engineering, and ferromagnetism. In this review, we summarize the methods to prepare moiré superlattices in the van der Waals materials and focus on the current discoveries of moiré pattern-modified electrical properties, recent findings of atomic reconstruction, as well as some possible future directions in this field.

Hetero-site nucleation for growing twisted bilayer graphene with a wide range of twist angles

Twisted bilayer graphene (tBLG) has recently attracted growing interest due to its unique twist-angle-dependent electronic properties. The preparation of high-quality large-area bilayer graphene with rich rotation angles would be important for the investigation of angle-dependent physics and applications, which, however, is still challenging. Here, we demonstrate a chemical vapor deposition (CVD) approach for growing high-quality tBLG using a hetero-site nucleation strategy, which enables the nucleation of the second layer at a different site from that of the first layer. The fraction of tBLGs in bilayer graphene domains with twist angles ranging from 0° to 30° was found to be improved to 88%, which is significantly higher than those reported previously. The hetero-site nucleation behavior was carefully investigated using an isotope-labeling technique. Furthermore, the clear Moiré patterns and ultrahigh room-temperature carrier mobility of 68,000 cm² V^-1 s^-1 confirmed the high crystalline quality of our tBLG. Our study opens an avenue for the controllable growth of tBLGs for both fundamental research and practical applications.

Tunable Optical Properties of Thin Films Controlled by the Interface Twist Angle

Control of materials properties has been the driving force of modern technologies. So far, materials properties have been modulated by their composition, structure, and size. Here, by using cathodoluminescence in a scanning transmission electron microscope, we show that the optical properties of stacked, >100 nm thick hexagonal boron nitride (hBN) films can be continuously tuned by their relative twist angles. Due to the formation of a moiré superlattice between the two interface layers of the twisted films, a new moiré sub-band gap is formed with continuously decreasing magnitude as a function of the twist angle, resulting in tunable luminescence wavelength and intensity increase of >40×. Our results demonstrate that moiré phenomena extend beyond monolayer-based systems and can be preserved in a technologically relevant, bulklike material at room temperature, dominating optical properties of hBN films for applications in medicine, environmental, or information technologies.

Fabrication Strategies of Twisted Bilayer Graphenes and Their Unique Properties

Twisted bilayer graphene (tBLG) exhibits a host of innovative physical phenomena owing to the formation of moiré superlattice. Especially, the discovery of superconducting behavior has generated new interest in graphene. The growing studies of tBLG mainly focus on its physical properties, while the fabrication of high-quality tBLG is a prerequisite for achieving the desired properties due to the great dependence on the twist angle and the interfacial contact. Here, the cutting-edge preparation strategies and challenges of tBLG fabrication are reviewed. The advantages and disadvantages of chemical vapor deposition, epitaxial growth on silicon carbide, stacking monolayer graphene, and folding monolayer graphene methods for the fabrication of tBLG are analyzed in detail, providing a reference for further development of preparation methods. Moreover, the characterization methods of twist angle for the tBLG are presented. Then, the unique physicochemical properties and corresponding applications of tBLG, containing correlated insulating and superconducting states, ferromagnetic state, soliton, enhanced optical absorption, tunable bandgap, and lithium intercalation and diffusion, are described. Finally, the opportunities and challenges for fabricating high-quality and large-area tBLG are discussed, unique physical properties are displayed, and new applications inferred from its angle-dependent features are explored, thereby impelling the commercialization of tBLG from laboratory to market.

Chern insulators, van Hove singularities and topological flat bands in magic-angle twisted bilayer graphene

Magic-angle twisted bilayer graphene exhibits intriguing quantum phase transitions triggered by enhanced electron-electron interactions when its flat bands are partially filled. However, the phases themselves and their connection to the putative non-trivial topology of the flat bands are largely unexplored. Here we report transport measurements revealing a succession of doping-induced Lifshitz transitions that are accompanied by van Hove singularities, which facilitate the emergence of correlation-induced gaps and topologically non-trivial subbands. In the presence of a magnetic field, well-quantized Hall plateaus at a filling of 1,2,3 carriers per moiré cell reveal the subband topology and signal the emergence of Chern insulators with Chern numbers, C = 3,2,1, respectively. Surprisingly, for magnetic fields exceeding 5 T we observe a van Hove singularity at a filling of 3.5, suggesting the possibility of a fractional Chern insulator. This van Hove singularity is accompanied by a crossover from low-temperature metallic, to high-temperature insulating behaviour, characteristic of entropically driven Pomeranchuk-like transitions.

Fermi Velocity Reduction of Dirac Fermions around the Brillouin Zone Center in In2 Se3 -Bilayer Graphene Heterostructures

Emergent phenomena such as unconventional superconductivity, Mott-like insulators, and the peculiar quantum Hall effect in graphene-based heterostructures are proposed to stem from the superlattice-induced renormalization of (moiré) Dirac fermions at the graphene Brillouin zone corners. Understanding the corresponding band structure commonly demands photoemission spectroscopy with both sub-meV resolution and large-momentum coverage, beyond the capability of the current state-of-the-art. Here the realization of moiré Dirac cones around the Brillouin zone center in monolayer In₂ Se₃ /bilayer graphene heterostructure is reported. The renormalization is evidenced by reduced Fermi velocity (≈23%) of the moiré Dirac cones and the reshaped Dirac point at the Γ point where they intersect. While there have been many theoretical predictions and much indirect experimental evidence, the findings here are the first direct observation of Fermi velocity reduction of the moiré Dirac cones. These features suggest strong In₂ Se₃ /graphene interlayer coupling, which is comparable with that in twisted bilayer graphene. The strategy expands the choice of materials in the heterostructure design and stimulates subsequent broad investigations of emergent physics at the sub-meV energy scale.

Spin-Twisted Optical Lattices: Tunable Flat Bands and Larkin-Ovchinnikov Superfluids

Moiré superlattices in twisted bilayer graphene and transition-metal dichalcogenides have emerged as a powerful tool for engineering novel band structures and quantum phases of two-dimensional quantum materials. Here we investigate Moiré physics emerging from twisting two independent hexagonal optical lattices of atomic (pseudo-)spin states (instead of bilayers) that exhibit remarkably different physics from twisted bilayer graphene. We employ a momentum-space tight-binding calculation that includes all range real-space tunnelings and show that all twist angles θ≲6° can become magic and support gapped flat bands. Because of the greatly enhanced density of states near the flat bands, the system can be driven to superfluidity by weak attractive interaction. Strikingly, the superfluid phase corresponds to a Larkin-Ovchinnikov state with finite momentum pairing that results from the interplay between flat bands and interspin interactions in the unique single-layer spin-twisted lattice. Our work may pave the way for exploring novel quantum phases and twistronics in cold atomic systems.

Moiré-Potential-Induced Band Structure Engineering in Graphene and Silicene

A moiré pattern results from the projection of one periodic pattern to another with relative lattice constant or misalignment and provides great periodic potential to modify the electronic properties of pristine materials. In this Review, recent research on the effect of the moiré superlattice on the electronic structures of graphene and silicene, both of which possess a honeycomb lattice, is focused on. The moiré periodic potential is introduced by the interlayer interaction to realize abundant phenomena, including new generation of Dirac cones, emergence of Van Hove singularities (vHs) at the cross point of two sets of Dirac cones, Mott-like insulating behavior at half-filling state, unconventional superconductivity, and electronic Kagome lattice and flat band with nontrivial edge state. The role of interlayer coupling strength, which is determined by twist angle and buckling degree, in these exotic properties is discussed in terms of both the theoretical prediction and experimental measurement, and finally, the challenges and outlook for this field are discussed.

Strain-controlled electronic and magnetic properties of tVS2/hVS2 van der Waals heterostructures

The structural, electronic and magnetic properties of the T-phase and H-phase of the VS₂ monolayer and their heterobilayers are studied by means of first-principles calculations. We find that the two phases of the VS₂ monolayer are both ferromagnetic (FM) semiconductors and that these two monolayers form a typical van der Waals (vdW) heterostructure with a weak interlayer interaction. By comparing the energy of different magnetic configurations, the FM state of the tVS₂/hVS₂ heterostructure is found to be in the ground state under normal conditions or biaxial strains. Under compressive strains, the anti-FM (AFM) and FM states degenerate. Based on the band structure obtained and the work function, it is found that the T-phase and H-phase are capable of forming an efficient p-n heterostructure. Due to spontaneous charge transfer at the interface, a gapless semiconductor is formed in our HSE06 calculations. We also find that the twist angle between the monolayers has a negligible impact on the band structure of the heterostructure in its spin-down channel. Moreover, the tVS₂/hVS₂ heterostructure is found to switch from a gapless semiconductor to a metal or a half-metal under some given biaxial or uniaxial strains. Therefore, the heterostructure could be a half-metallic property with strains, realizing 100% polarization at the Fermi level. Our study provides the possibility of realizing 100% spin-polarization at the Fermi level in these FM vdW heterostructures, which is significant for further spin transport exploration.

Chirality-Induced Giant Unidirectional Magnetoresistance in Twisted Bilayer Graphene

Twisted bilayer graphene (TBG) exhibits fascinating correlation-driven phenomena like the superconductivity and Mott insulating state, with flat bands and a chiral lattice structure. We find by quantum-transport calculations that the chirality leads to a giant unidirectional magnetoresistance (UMR) in TBG, where the unidirectionality refers to the resistance change under the reversal of the direction of current or magnetic field. We point out that flat bands significantly enhance this effect. The UMR increases quickly upon reducing the twist angle, and reaches about 20% for an angle of 1.5° in a 10 T in-plane magnetic field. We propose the band structure topology (asymmetry), which leads to a direction-sensitive mean free path, as a useful way to anticipate the UMR effect. The UMR provides a probe for chirality and band flatness in the twisted bilayers.

Micro-photoluminescence of Carbon Dots Deposited on Twisted Double-Layer Graphene Grown by Chemical Vapor Deposition

Carbon-based nanomaterials, such as carbon dots (CDs) and graphene (Gr), feature outstanding optical and electronic properties. Hence, their integration in optoelectronic and photonic devices is easier thanks to their low dimensionality and offers the possibility to reach high-quality performances. In this context, the combination of CDs and Gr into new nanocomposite materials CDs/Gr can further improve their optoelectronic properties and eventually create new ones, paving the way for the development of advanced carbon nanotechnology. In this work, we have thoroughly investigated the structural and emission properties of CDs deposited on single-layer and bilayer graphene lying on a SiO₂/Si substrate. A systematic Raman analysis points out that bilayer (BL) graphene grown by chemical vapor deposition does not always respect the Bernal (AB) stacking, but it is rather a mixture of twisted bilayer (t-BL) featuring domains with different twist angles. Moreover, in-depth micro-photoluminescence measurements, combined with atomic force microscopy (AFM) morphological analysis, show that CD emission efficiency is strongly depleted by the presence of graphene and in particular is dependent on the number of layers as well as on the twist angle of BL graphene. Finally, we propose a model which explains these results on the basis of photoinduced charge-transfer processes, taking into account the energy levels of the hybrid nanosystem formed by coupling CDs with t-BL/SiO₂.

Higher-Order Band Topology in Twisted Moiré Superlattice

The two-dimensional (2D) twisted bilayer materials with van der Waals coupling have ignited great research interests, paving a new way to explore the emergent quantum phenomena by twist degree of freedom. Generally, with the decreasing of twist angle, the enhanced interlayer coupling will gradually flatten the low-energy bands and isolate them by two high-energy gaps at zero and full filling, respectively. Although the correlation and topological physics in the low-energy flat bands have been intensively studied, little information is available for these two emerging gaps. In this Letter, we predict a 2D second-order topological insulator (SOTI) for twisted bilayer graphene and twisted bilayer boron nitride in both zero and full filling gaps. Employing a tight-binding Hamiltonian based on first-principles calculations, three unique fingerprints of 2D SOTI are identified, that is, nonzero bulk topological index, gapped topological edge state, and in-gap topological corner state. Most remarkably, the 2D SOTI exists in a wide range of commensurate twist angles, which is robust to microscopic structure disorder and twist center, greatly facilitating the possible experimental measurement. Our results not only extend the higher-order band topology to massless and massive twisted moiré superlattice, but also demonstrate the importance of high-energy bands for fully understanding the nontrivial electronics.

Spontaneous Valley Spirals in Magnetically Encapsulated Twisted Bilayer Graphene

Van der Waals heterostructures provide a rich platform for emergent physics due to their tunable hybridization of layers, orbitals, and spin. Here, we find that twisted bilayer graphene stacked between antialigned ferromagnetic insulators can feature flat electronic bands due to the interplay between twist, exchange proximity, and spin-orbit coupling. These flat bands are nearly degenerate in valley only and are effectively described by a triangular superlattice model. At half filling, we find that interactions induce spontaneous valley correlations that favor spiral order and derive a low-energy valley-Heisenberg model with symmetric and antisymmetric exchange couplings. We also show how electric interlayer bias broadens the bands and tunes these couplings. Our results put forward magnetic van der Waals heterostructures as a platform to explore valley-correlated states.

Raman Spectroscopy of Twisted Bilayer Graphene

When two periodic two-dimensional structures are superposed, any mismatch rotation angle between the layers generates a Moiré pattern superlattice, whose size depends on the twisting angle θ. If the layers are composed by different materials, this effect is also dependent on the lattice parameters of each layer. Moiré superlattices are commonly observed in bilayer graphene, where the mismatch angle between layers can be produced by growing twisted bilayer graphene (TBG) samples by CVD or folding the monolayer back upon itself. In TBG, it was shown that the coupling between the Dirac cones of the two layers gives rise to van Hove singularities (vHs) in the density of electronic states, whose energies vary with θ. The understanding of the behavior of electrons and their interactions with phonons in atomically thin heterostructures is crucial for the engineering of novel 2D devices. Raman spectroscopy has been often used to characterize twisted bilayer graphene and graphene heterostructures. Here, we review the main important effects in the Raman spectra of TBG discussing firstly the appearance of new peaks in the spectra associated with phonons with wavevectors within the interior of the Brillouin zone of graphene corresponding to the reciprocal unit vectors of the Moiré superlattice, and that are folded to the center of the reduced Brillouin Zone (BZ) becoming Raman active. Another important effect is the giant enhancement of G band intensity of TBG that occurs only in a narrow range of laser excitation energies and for a given twisting angle. Results show that the vHs in the density of states is not only related to the folding of the commensurate BZ, but mainly associated with the Moiré pattern that does not necessarily have a translational symmetry. Finally, we show that there are two different resonance mechanisms that activate the appearance of the extra peaks: the intralayer and interlayer electron–phonon processes, involving electrons of the same layer or from different layers, respectively. Both effects are observed for twisted bilayer graphene, but Raman spectroscopy can also be used to probe the intralayer process in any kind of graphene-based heterostructure, like in the graphene/h-BN junctions.

Correlation-driven topological phases in magic-angle twisted bilayer graphene

Magic-angle twisted bilayer graphene (MATBG) exhibits a range of correlated phenomena that originate from strong electron-electron interactions. These interactions make the Fermi surface highly susceptible to reconstruction when ±1, ±2 and ±3 electrons occupy each moiré unit cell, and lead to the formation of various correlated phases^1-4. Although some phases have been shown to have a non-zero Chern number^5,6, the local microscopic properties and topological character of many other phases have not yet been determined. Here we introduce a set of techniques that use scanning tunnelling microscopy to map the topological phases that emerge in MATBG in a finite magnetic field. By following the evolution of the local density of states at the Fermi level with electrostatic doping and magnetic field, we create a local Landau fan diagram that enables us to assign Chern numbers directly to all observed phases. We uncover the existence of six topological phases that arise from integer fillings in finite fields and that originate from a cascade of symmetry-breaking transitions driven by correlations^7,8. These topological phases can form only for a small range of twist angles around the magic angle, which further differentiates them from the Landau levels observed near charge neutrality. Moreover, we observe that even the charge-neutrality Landau spectrum taken at low fields is considerably modified by interactions, exhibits prominent electron-hole asymmetry, and features an unexpectedly large splitting between zero Landau levels (about 3 to 5 millielectronvolts). Our results show how strong electronic interactions affect the MATBG band structure and lead to correlation-enabled topological phases.

Topological flat bands in twisted trilayer graphene

Twisted trilayer graphene (TLG) may be the simplest realistic system so far, which has flat bands with nontrivial topology. Here, we give a comprehensive calculation about its band structures and the band topology, i.e., valley Chern number of the nearly flat bands, with the continuum model. With realistic parameters, the magic angle of twisted TLG is about 1.12°, at which two nearly flat bands appears. Unlike the twisted bilayer graphene, a small twist angle can induce a tiny gap at all the Dirac points, which can be enlarged further by a perpendicular electric field. The valley Chern numbers of the two nearly flat bands in the twisted TLG depends on the twist angle θ and the perpendicular electric field E_⊥. Considering its topological flat bands, the twisted TLG should be an ideal experimental platform to study the strongly correlated physics in topologically nontrivial flat band systems. And, due to its reduced symmetry, the correlated states in twisted TLG should be quite different from that in twisted bilayer graphene and twisted double bilayer graphene.