Mixing of moiré-surface and bulk states in graphite.

Van der Waals assembly enables the design of electronic states in two-dimensional (2D) materials, often by superimposing a long-wavelength periodic potential on a crystal lattice using moiré superlattices1-9. This twistronics approach has resulted in numerous previously undescribed physics, including strong correlations and superconductivity in twisted bilayer graphene10-12, resonant excitons, charge ordering and Wigner crystallization in transition-metal chalcogenide moiré structures13-18 and Hofstadter's butterfly spectra and Brown-Zak quantum oscillations in graphene superlattices19-22. Moreover, twistronics has been used to modify near-surface states at the interface between van der Waals crystals23,24. Here we show that electronic states in three-dimensional (3D) crystals such as graphite can be tuned by a superlattice potential occurring at the interface with another crystal-namely, crystallographically aligned hexagonal boron nitride. This alignment results in several Lifshitz transitions and Brown-Zak oscillations arising from near-surface states, whereas, in high magnetic fields, fractal states of Hofstadter's butterfly draw deep into the bulk of graphite. Our work shows a way in which 3D spectra can be controlled using the approach of 2D twistronics.

Electronic phase separation in multilayer rhombohedral graphite.

Of the two stable forms of graphite, hexagonal and rhombohedral, the former is more common and has been studied extensively. The latter is less stable, which has so far precluded its detailed investigation, despite many theoretical predictions about the abundance of exotic interaction-induced physics1-6. Advances in van der Waals heterostructure technology7 have now allowed us to make high-quality rhombohedral graphite films up to 50 graphene layers thick and study their transport properties. Here we show that the bulk electronic states in such rhombohedral graphite are gapped8 and, at low temperatures, electron transport is dominated by surface states. Because of their proposed topological nature, the surface states are of sufficiently high quality to observe the quantum Hall effect, whereby rhombohedral graphite exhibits phase transitions between a gapless semimetallic phase and a gapped quantum spin Hall phase with giant Berry curvature. We find that an energy gap can also be opened in the surface states by breaking their inversion symmetry by applying a perpendicular electric field. Moreover, in rhombohedral graphite thinner than four nanometres, a gap is present even without an external electric field. This spontaneous gap opening shows pronounced hysteresis and other signatures characteristic of electronic phase separation, which we attribute to emergence of strongly correlated electronic surface states.

Resonantly hybridized excitons in moiré superlattices in van der Waals heterostructures.

Atomically thin layers of two-dimensional materials can be assembled in vertical stacks that are held together by relatively weak van der Waals forces, enabling coupling between monolayer crystals with incommensurate lattices and arbitrary mutual rotation1,2. Consequently, an overarching periodicity emerges in the local atomic registry of the constituent crystal structures, which is known as a moiré superlattice3. In graphene/hexagonal boron nitride structures4, the presence of a moiré superlattice can lead to the observation of electronic minibands5-7, whereas in twisted graphene bilayers its effects are enhanced by interlayer resonant conditions, resulting in a superconductor-insulator transition at magic twist angles8. Here, using semiconducting heterostructures assembled from incommensurate molybdenum diselenide (MoSe2) and tungsten disulfide (WS2) monolayers, we demonstrate that excitonic bands can hybridize, resulting in a resonant enhancement of moiré superlattice effects. MoSe2 and WS2 were chosen for the near-degeneracy of their conduction-band edges, in order to promote the hybridization of intra- and interlayer excitons. Hybridization manifests through a pronounced exciton energy shift as a periodic function of the interlayer rotation angle, which occurs as hybridized excitons are formed by holes that reside in MoSe2 binding to a twist-dependent superposition of electron states in the adjacent monolayers. For heterostructures in which the monolayer pairs are nearly aligned, resonant mixing of the electron states leads to pronounced effects of the geometrical moiré pattern of the heterostructure on the dispersion and optical spectra of the hybridized excitons. Our findings underpin strategies for band-structure engineering in semiconductor devices based on van der Waals heterostructures9.

Author Correction: Resonantly hybridized excitons in moiré superlattices in van der Waals heterostructures.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Competition of Moiré Network Sites to Form Electronic Quantum Dots in Reconstructed MoX2/WX2 Heterostructures.

Twisted bilayers of two-dimensional semiconductors offer a versatile platform for engineering quantum states for charge carriers using moiré superlattice effects. Among the systems of recent interest are twistronic MoX2/WX2 heterostructures (X = Se or S), which undergo reconstruction into preferential stacking domains and highly strained domain wall networks, determining the electron/hole localization across moiré superlattices. Here, we present a catalogue of options for the formation of self-organized quantum dots and wires in lattice-reconstructed marginally twisted MoX2/WX2 bilayers with a relative lattice mismatch δ ≪ 1 for twist angles ranging from perfect alignment to θ ∼ 1°. On the basis of multiscale modeling taking into account twirling of domain wall networks, we analyze bilayers with both parallel and antiparallel orientations of their unit cells and describe crossovers between different positioning of band edges for electrons and holes across moiré superlattices when θ < δ and θ > δ.

Twirling and Spontaneous Symmetry Breaking of Domain Wall Networks in Lattice-Reconstructed Heterostructures of Two-Dimensional Materials.

Lattice relaxation in twistronic bilayers with close lattice parameters and almost perfect crystallographic alignment of the layers results in the transformation of the moiré pattern into a sequence of preferential stacking domains and domain wall networks. Here, we show that reconstructed moiré superlattices of the perfectly aligned heterobilayers of same chalcogen transition metal dichalcogenides have broken-symmetry structures featuring twisted nodes ("twirls") of domain wall networks. The analysis of twist-angle dependence of strain characteristics for the broken-symmetry structures shows that the formation of twirl reduces the amount of hydrostatic strain around the nodes, potentially weakening their influence on the band edge energies of electrons and holes.

Mixed-Stacking Few-Layer Graphene as an Elemental Weak Ferroelectric Material.

Ferroelectricity (Valasek, J. Phys. Rev. 1921, 17, 475), a spontaneous formation of electric polarization, is a solid state phenomenon, usually, associated with ionic compounds or complex materials. Here we show that, atypically for elemental solids, few-layer graphenes can host an equilibrium out-of-plane electric polarization, switchable by sliding the constituent graphene sheets. The systems hosting such effect include mixed-stacking tetralayers and thicker (5-9 layers) rhombohedral graphitic films with a twin boundary in the middle of a flake. The predicted electric polarization would also appear in marginally (small-angle) twisted few-layer flakes, where lattice reconstruction would give rise to networks of mesoscale domains with alternating value and sign of out-of-plane polarization.

Mid-Infrared Mapping of Four-Layer Graphene Polytypes Using Near-Field Microscopy.

The mid-infrared (MIR) spectral region attracts attention for accurate chemical analysis using photonic devices. Few-layer graphene (FLG) polytypes are promising platforms, due to their broad absorption in this range and gate-tunable optical properties. Among these polytypes, the noncentrosymmetric ABCB/ACAB structure is particularly interesting, due to its intrinsic bandgap (8.8 meV) and internal polarization. In this study, we utilize scattering-scanning near-field microscopy to measure the optical response of all three tetralayer graphene polytypes in the 8.5-11.5 µm range. We employ a finite dipole model to compare these results to the calculated optical conductivity for each polytype obtained from a tight-binding model. Our findings reveal a significant discrepancy in the MIR optical conductivity response of graphene between the different polytypes than what the tight-binding model suggests. This observation implies an increased potential for utilizing the distinct tetralayer polytypes in photonic devices operating within the MIR range for chemical sensing and infrared imaging.

ARPES Signatures of Few-Layer Twistronic Graphenes.

Diverse emergent correlated electron phenomena have been observed in twisted-graphene layers. Many electronic structure predictions have been reported exploring this new field, but with few momentum-resolved electronic structure measurements to test them. We use angle-resolved photoemission spectroscopy to study the twist-dependent (1° < θ < 8°) band structure of twisted-bilayer, monolayer-on-bilayer, and double-bilayer graphene (tDBG). Direct comparison is made between experiment and theory, using a hybrid k·p model for interlayer coupling. Quantitative agreement is found across twist angles, stacking geometries, and back-gate voltages, validating the models and revealing field-induced gaps in twisted graphenes. However, for tDBG at θ = 1.5 ± 0.2°, close to the magic angle θ = 1.3°, a flat band is found near the Fermi level with measured bandwidth Ew = 31 ± 5 meV. An analysis of the gap between the flat band and the next valence band shows deviations between experiment (Δh = 46 ± 5 meV) and theory (Δh = 5 meV), indicative of lattice relaxation in this regime.

A Scalable Network Model for Electrically Tunable Ferroelectric Domain Structure in Twistronic Bilayers of Two-Dimensional Semiconductors.

Moiré structures in small-angle-twisted bilayers of two-dimensional (2D) semiconductors with a broken-symmetry interface form arrays of ferroelectric (FE) domains with periodically alternating out-of-plane polarization. Here, we propose a network theory for the tunability of such FE domain structure by applying an electric field perpendicular to the 2D crystal. Using multiscale analysis, we derive a fully parametrized string-theory-like description of the domain wall network (DWN) and show that it undergoes a qualitative change, after the arcs of partial dislocation (PD) like domain walls merge (near the network nodes) into streaks of perfect screw dislocations (PSD), which happens at a threshold displacement field dependent on the DWN period.

Band Gap Opening in Bilayer Graphene-CrCl3/CrBr3/CrI3 van der Waals Interfaces.

We report experimental investigations of transport through bilayer graphene (BLG)/chromium trihalide (CrX3; X = Cl, Br, I) van der Waals interfaces. In all cases, a large charge transfer from BLG to CrX3 takes place (reaching densities in excess of 1013 cm-2), and generates an electric field perpendicular to the interface that opens a band gap in BLG. We determine the gap from the activation energy of the conductivity and find excellent agreement with the latest theory accounting for the contribution of the σ bands to the BLG dielectric susceptibility. We further show that for BLG/CrCl3 and BLG/CrBr3 the band gap can be extracted from the gate voltage dependence of the low-temperature conductivity, and use this finding to refine the gap dependence on the magnetic field. Our results allow a quantitative comparison of the electronic properties of BLG with theoretical predictions and indicate that electrons occupying the CrX3 conduction band are correlated.

Angstrom-Scale Electrochemistry at Electrodes with Dimensions Commensurable and Smaller than Individual Reacting Species.

In nature and technologies, many chemical reactions occur at interfaces with dimensions approaching that of a single reacting species in nano- and angstrom-scale. Mechanisms governing reactions at this ultimately small spatial regime remain poorly explored because of challenges to controllably fabricate required devices and assess their performance in experiment. Here we report how efficiency of electrochemical reactions evolves for electrodes that range from just one atom in thickness to sizes comparable with and exceeding hydration diameters of reactant species. The electrodes are made by encapsulating graphene and its multilayers within insulating crystals so that only graphene edges remain exposed and partake in reactions. We find that limiting current densities characterizing electrochemical reactions exhibit a pronounced size effect if reactant's hydration diameter becomes commensurable with electrodes' thickness. An unexpected blockade effect is further revealed from electrodes smaller than reactants, where incoming reactants are blocked by those adsorbed temporarily at the atomically narrow interfaces. The demonstrated angstrom-scale electrochemistry offers a venue for studies of interfacial behaviors at the true molecular scale.

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 Å3 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.

Design of van der Waals interfaces for broad-spectrum optoelectronics.

Van der Waals (vdW) interfaces based on 2D materials are promising for optoelectronics, as interlayer transitions between different compounds allow tailoring of the spectral response over a broad range. However, issues such as lattice mismatch or a small misalignment of the constituent layers can drastically suppress electron-photon coupling for these interlayer transitions. Here, we engineered type-II interfaces by assembling atomically thin crystals that have the bottom of the conduction band and the top of the valence band at the Γ point, and thus avoid any momentum mismatch. We found that these van der Waals interfaces exhibit radiative optical transitions irrespective of the lattice constant, the rotational and/or translational alignment of the two layers or whether the constituent materials are direct or indirect gap semiconductors. Being robust and of general validity, our results broaden the scope of future optoelectronics device applications based on two-dimensional materials.

Control of Giant Topological Magnetic Moment and Valley Splitting in Trilayer Graphene.

Bloch states of electrons in honeycomb two-dimensional crystals with multivalley band structure and broken inversion symmetry have orbital magnetic moments of a topological nature. In crystals with two degenerate valleys, a perpendicular magnetic field lifts the valley degeneracy via a Zeeman effect due to these magnetic moments, leading to magnetoelectric effects which can be leveraged for creating valleytronic devices. In this work, we demonstrate that trilayer graphene with Bernal stacking (ABA TLG), hosts topological magnetic moments with a large and widely tunable valley g factor (g_{ν}), reaching a value g_{ν}∼1050 at the extreme of the studied parametric range. The reported experiment consists in sublattice-resolved scanning tunneling spectroscopy under perpendicular electric and magnetic fields that control the TLG bands. The tunneling spectra agree very well with the results of theoretical modeling that includes the full details of the TLG tight-binding model and accounts for a quantum-dot-like potential profile formed electrostatically under the scanning tunneling microscope tip.

Raman spectroscopy of GaSe and InSe post-transition metal chalcogenides layers.

III-VI post-transition metal chalcogenides (InSe and GaSe) are a new class of layered semiconductors, which feature a strong variation of size and type of their band gaps as a function of number of layers (N). Here, we investigate exfoliated layers of InSe and GaSe ranging from bulk crystals down to monolayer, encapsulated in hexagonal boron nitride, using Raman spectroscopy. We present the N-dependence of both intralayer vibrations within each atomic layer, as well as of the interlayer shear and layer breathing modes. A linear chain model can be used to describe the evolution of the peak positions as a function of N, consistent with first principles calculations.

Spectroscopic Signatures of Electronic Excitations in Raman Scattering in Thin Films of Rhombohedral Graphite.

Rhombohedral graphite features peculiar electronic properties, including persistence of low-energy surface bands of a topological nature. Here, we study the contribution of electron-hole excitations toward inelastic light scattering in thin films of rhombohedral graphite. We show that, in contrast to the featureless electron-hole contribution toward Raman spectrum of graphitic films with Bernal stacking, the inelastic light scattering accompanied by electron-hole excitations in crystals with rhombohedral stacking produces distinct features in the Raman signal which can be used both to identify the stacking and to determine the number of layers in the film.

Topologically Nontrivial Valley States in Bilayer Graphene Quantum Point Contacts.

We present measurements of quantized conductance in electrostatically induced quantum point contacts in bilayer graphene. The application of a perpendicular magnetic field leads to an intricate pattern of lifted and restored degeneracies with increasing field: at zero magnetic field the degeneracy of quantized one-dimensional subbands is four, because of a twofold spin and a twofold valley degeneracy. By switching on the magnetic field, the valley degeneracy is lifted. Because of the Berry curvature, states from different valleys split linearly in magnetic field. In the quantum Hall regime fourfold degenerate conductance plateaus reemerge. During the adiabatic transition to the quantum Hall regime, levels from one valley shift by two in quantum number with respect to the other valley, forming an interweaving pattern that can be reproduced by numerical calculations.

Oscillating Magnetoresistance in Graphene p-n Junctions at Intermediate Magnetic Fields.

We report on the observation of magnetoresistance oscillations in graphene p-n junctions. The oscillations have been observed for six samples, consisting of single-layer and bilayer graphene, and persist up to temperatures of 30 K, where standard Shubnikov-de Haas oscillations are no longer discernible. The oscillatory magnetoresistance can be reproduced by tight-binding simulations. We attribute this phenomenon to the modulated densities of states in the n- and p-regions.

Realization of Valley and Spin Pumps by Scattering at Nonmagnetic Disorders.

The recent success in optical pumping of valley polarization in two-dimensional transition metal dichalcogenides (TMDs) has greatly promoted the concept of valley-based informatics and electronics. However, between the demonstrated valley polarization of transient electron-hole pair excitations and practical valleytronic operations, there exist obvious gaps to fill, among which is the valley pump of long-lived charge carriers. Here we discover that the quested valley pump of electrons or holes can be realized simply by scattering at the ubiquitous nonmagnetic disorders, not relying on any specific material property. The mechanism is rooted in the nature of the valley as a momentum space index: the intervalley backscattering in general has a valley contrasted rate due to the distinct momentum transfers, causing a net transfer of population from one valley to another. As examples, we numerically demonstrate the sizable valley pump effects driven by charge current in nanoribbons of monolayer TMDs, where the spin-orbit scattering by nonmagnetic disorders also realizes a spin pump for the spin-valley locked holes. Our finding points to a new opportunity towards valley spintronics, turning disorders from a deleterious factor to a resource of valley and spin polarization.