When two-dimensional van der Waals materials are stacked to build heterostructures, moir\'e patterns emerge from twisted interfaces or from mismatch in lattice constant of individual layers. Relaxation of the atomic positions is a direct, generic consequence of the moir\'e pattern, with many implications for the physical properties. Moir\'e driven atomic relaxation may be naively thought to be restricted to the interfacial layers and thus irrelevant for multi-layered heterostructures. However, we provide experimental evidence for the importance of the three dimensional nature of the relaxation in two types of van der Waals heterostructures: First, in multi-layer graphene twisted on graphite at small twist angles ($\theta\approx0.14^\circ$) we observe propagation of relaxation domains even beyond 18 graphene layers. Second, we show how for multi-layer PdTe$_2$ on Bi$_2$Se$_3$ the moir\'e lattice constant depends on the number of PdTe$_2$ layers. Motivated by the experimental findings, we developed a continuum approach to model multi-layered relaxation processes based on the generalized stacking fault energy functional given by ab-initio simulations. Leveraging the continuum property of the approach enables us to access large scale regimes and achieve agreement with our experimental data for both systems. Furthermore it is well known that the electronic structure of graphene sensitively depends on local lattice deformations. Therefore we study the impact of multi-layered relaxation on the local density of states of the twisted graphitic system. We identify measurable implications for the system, experimentally accessible by scanning tunneling microscopy. Our multi-layered relaxation approach is not restricted to the discussed systems, and can be used to uncover the impact of an interfacial defect on various layered systems of interest.

Twisted van der Waals multilayers are widely regarded as a rich platform to access novel electronic phases thanks to the multiple degrees of freedom available for controlling their electronic and chemical properties. Here, we propose that the stacking domains that form naturally due to the relative twist between successive layers act as an additional "knob" for controlling the behavior of these systems and report the emergence and engineering of stacking domain-dependent surface chemistry in twisted few-layer graphene. Using mid-infrared near-field optical microscopy and atomic force microscopy, we observe a selective adhesion of metallic nanoparticles and liquid water at the domains with rhombohedral stacking configurations of minimally twisted double bi- and trilayer graphene. Furthermore, we demonstrate that the manipulation of nanoparticles located at certain stacking domains can locally reconfigure the moiré superlattice in their vicinity at the micrometer scale. Our findings establish a new approach to controlling moiré-assisted chemistry and nanoengineering.

Magnetic fields can have profound effects on the motion of electrons in quantum materials. Two-dimensional electron systems subject to strong magnetic fields are expected to exhibit quantized Hall conductivity, chiral edge currents and distinctive collective modes referred to as magnetoplasmons and magnetoexcitons. Generating these propagating collective modes in charge-neutral samples and imaging them at their native nanometre length scales have thus far been experimentally elusive. Here we visualize propagating magnetoexciton polaritons at their native length scales and report their magnetic-field-tunable dispersion in near-charge-neutral graphene. Imaging these collective modes and their associated nano-electro-optical responses allows us to identify polariton-modulated optical and photo-thermal electric effects at the sample edges, which are the most pronounced near charge neutrality. Our work is enabled by innovations in cryogenic near-field optical microscopy techniques that allow for the nano-imaging of the near-field responses of two-dimensional materials under magnetic fields up to 7 T. This nano-magneto-optics approach allows us to explore and manipulate magnetopolaritons in specimens with low carrier doping via harnessing high magnetic fields. Dirac magnetoexcitons with non-trivial nanoscale electrodynamics are formed from the excitation of Landau levels in charge-neutral graphene. Here, the Dirac magnetoexciton dispersion is directly imaged up to 7 T via a magneto cryogenic near-field microscope.

We report the epitaxial growth of PdTe 2 ultrathin films on a topological insulator Bi 2 Se 3 . A prominent moiré pattern was observed in scanning tunneling microscope measurements. The moiré periodicity increases as film thickness decreases, indicating a lattice expansion of epitaxial PdTe 2 thin films at lower thicknesses. In addition, our simulations based on a multilayer relaxation technique reveal uniaxial lattice strains at the edge of PdTe 2 domains, and anisotropic strain distributions throughout the moiré supercell with a net change in lattice strain up to ∼2.9%. Our density functional theory calculations show that this strain effect leads to a narrowing of the band gap at Γ point near the Fermi level. Under a strain of ∼2.8%, the band gap at Γ closes completely. Further increasing the lattice strain makes the band gap reopen and the order of conduction band and valence bands inverted in energy. The experimental and theoretical results shed light on a method for constructing quantum grids of topological band structure under the modulation of moiré potentials.

Ferroelectricity, a spontaneous and reversible electric polarization, is found in certain classes of van der Waals (vdW) materials. The discovery of ferroelectricity in twisted vdW layers provides new opportunities to engineer spatially dependent electric and optical properties associated with the configuration of moiré superlattice domains and the network of domain walls. Here, we employ near-field infrared nano-imaging and nano-photocurrent measurements to study ferroelectricity in minimally twisted WSe 2 . The ferroelectric domains are visualized through the imaging of the plasmonic response in a graphene monolayer adjacent to the moiré WSe 2 bilayers. Specifically, we find that the ferroelectric polarization in moiré domains is imprinted on the plasmonic response of the graphene. Complementary nano-photocurrent measurements demonstrate that the optoelectronic properties of graphene are also modulated by the proximal ferroelectric domains. Our approach represents an alternative strategy for studying moiré ferroelectricity at native length scales and opens promising prospects for (opto)electronic devices.

We report the epitaxial growth of PdTe 2 ultrathin films on topological insulator Bi 2 Se 3 . A prominent moiré pattern was observed in STM measurements. The moiré periodicity increases as film thickness decreases, indicating a lattice expansion of epitaxial PdTe 2 thin films with lower thicknesses. In addition, our simulations based on a multilayer relaxation technique reveal uniaxial lattice strains at the edge of PdTe 2 domains, and anisotropic strain distributions throughout the moiré supercell with a net change in lattice strain up to ~2.9%. Our DFT calculations show that this strain effect leads to a narrowing of the band gap at Γ point near the Fermi level. Under a strain of ~2.8%, the band gap at Γ closes completely. Further increasing the lattice strain makes the band gap reopen and the order of conduction band and valence bands inverted in energy. The results offer a proof of concept for constructing quantum grids of topological materials under the modulation of moiré potentials.

When two atomic layers are brought into contact at a relative twist angle, a large-scale pattern, called a moiré superlattice, emerges due to the (angular or lattice) mismatch between the layers. This has profound consequences in terms of the Hamiltonian of the system but was also considered in several publications as a means to extract the local strain tensor. While extracting the twist angle based on knowledge of the periodicity of the moiré is trivial in the case of a regular moiré pattern, in many examples in the literature, that is not the case. In particular, extracting the strain tensor and twist angle maps from a spatially varying moiré pattern is not straightforward. This article aims to provide a practical tool to extract the strain tensor and twist angle from an experimentally observable pattern. It further addresses the limitation of any such approach in the absence of additional experimental information beyond the moiré superlattice pattern.

The electronic and structural properties of atomically thin materials can be controllably tuned by assembling them with an interlayer twist. During this process, constituent layers spontaneously rearrange themselves in search of a lowest energy configuration. Such relaxation phenomena can lead to unexpected and novel material properties. Here, we study twisted double trilayer graphene (TDTG) using nano-optical and tunneling spectroscopy tools. We reveal a surprising optical and electronic contrast, as well as a stacking energy imbalance emerging between the moiré domains. We attribute this contrast to an unconventional form of lattice relaxation in which an entire graphene layer spontaneously shifts position during assembly, resulting in domains of ABABAB and BCBACA stacking. We analyze the energetics of this transition and demonstrate that it is the result of a non-local relaxation process, in which an energy gain in one domain of the moiré lattice is paid for by a relaxation that occurs in the other.

Broken symmetries induce strong even-order nonlinear optical responses in materials and at interfaces. Unlike conventional covalently bonded nonlinear crystals, van der Waals (vdW) heterostructures feature layers that can be stacked at arbitrary angles, giving complete control over the presence or lack of inversion symmetry at a crystal interface. Here, we report highly tunable second harmonic generation (SHG) from nanomechanically rotatable stacks of bulk hexagonal boron nitride (BN) crystals and introduce the term twistoptics to describe studies of optical properties in twistable vdW systems. By suppressing residual bulk effects, we observe SHG intensity modulated by a factor of more than 50, and polarization patterns determined by moiré interface symmetry. Last, we demonstrate greatly enhanced conversion efficiency in vdW vertical superlattice structures with multiple symmetry-broken interfaces. Our study paves the way for compact twistoptics architectures aimed at efficient tunable frequency conversion and demonstrates SHG as a robust probe of buried vdW interfaces.

Significance Micrometer-sized uniform four-layer (ABCA) rhombohedral graphene is realized by introducing a small twist angle between two bilayers of Bernal graphene. By means of scanning tunneling spectroscopy we observe an extremely sharp van Hove singularity of 3–5-meV half-width and a correlated many-body gap of 9.5 meV at neutrality, thus making small twisted double-bilayer graphene a unique platform to realize electronic correlations in the absence of a moiré potential. Furthermore, ABCA graphene domain walls display tunable topological edge states, of great interest in Floquet engineering.

The emerging field of twistronics, which harnesses the twist angle between two-dimensional materials, represents a promising route for the design of quantum materials, as the twist-angle-induced superlattices offer means to control topology and strong correlations. At the small twist limit, and particularly under strain, as atomic relaxation prevails, the emergent moiré superlattice encodes elusive insights into the local interlayer interaction. Here we introduce moiré metrology as a combined experiment-theory framework to probe the stacking energy landscape of bilayer structures at the 0.1 meV/atom scale, outperforming the gold-standard of quantum chemistry. Through studying the shapes of moiré domains with numerous nano-imaging techniques, and correlating with multi-scale modelling, we assess and refine first-principle models for the interlayer interaction. We document the prowess of moiré metrology for three representative twisted systems: bilayer graphene, double bilayer graphene and H-stacked MoSe_2/WSe_2. Moiré metrology establishes sought after experimental benchmarks for interlayer interaction, thus enabling accurate modelling of twisted multilayers.

In twisted bilayers of semiconducting transition metal dichalcogenides, a combination of structural rippling and electronic coupling gives rise to periodic moire potentials that can confine charged and neutral excitations(1-5). Here we show that the moire potential in these bilayers at small angles is unexpectedly large, reaching values above 300 meV for the valence band and 150 meV for the conduction band-an order of magnitude larger than theoretical estimates based on interlayer coupling alone. We further demonstrate that the moire potential is a non-monotonic function of moire wavelength, reaching a maximum at a moire period of similar to 13 nm . This non-monotonicity coincides with a change in the structure of the moire pattern from a continuous variation of stacking order at small moire wavelengths to a one-dimensional soliton-dominated structure at large moire wavelengths. We show that the in-plane structure of the moire pattern is captured by a continuous mechanical relaxation model, and find that the moire structure and internal strain, rather than the interlayer coupling, are the dominant factors in determining the moire potential. Our results demonstrate the potential of using precision moire structures to create deeply trapped carriers or excitations for quantum electronics and opto-electronics.

Chiral Weyl fermions with linear energy-momentum dispersion in the bulk accompanied by Fermi-arc states on the surfaces prompt a host of enticing optical effects. While new Weyl semimetal materials keep emerging, the available optical probes are limited. In particular, isolating bulk and surface electrodynamics in Weyl conductors remains a challenge. We devised an approach to the problem based on near-field photocurrent imaging at the nanoscale and applied this technique to a prototypical Weyl semimetal TaIrTe4 As a first step, we visualized nano-photocurrent patterns in real space and demonstrated their connection to bulk nonlinear conductivity tensors through extensive modeling augmented with density functional theory calculations. Notably, our nanoscale probe gives access to not only the in-plane but also the out-of-plane electric fields so that it is feasible to interrogate all allowed nonlinear tensors including those that remained dormant in conventional far-field optics. Surface- and bulk-related nonlinear contributions are distinguished through their "symmetry fingerprints" in the photocurrent maps. Robust photocurrents also appear at mirror-symmetry breaking edges of TaIrTe4 single crystals that we assign to nonlinear conductivity tensors forbidden in the bulk. Nano-photocurrent spectroscopy at the boundary reveals a strong resonance structure absent in the interior of the sample, providing evidence for elusive surface states.

Graphene-based heterostructures display a variety of phenomena that are strongly tunable by electrostatic local gates. Monolayer graphene (MLG) exhibits tunable surface plasmon polaritons, as revealed by scanning nano-infrared experiments. In bilayer graphene (BLG), an electronic gap is induced by a perpendicular displacement field. Gapped BLG is predicted to display unusual effects such as plasmon amplification and domain wall plasmons with significantly larger lifetime than MLG. Furthermore, a variety of correlated electronic phases highly sensitive to displacement fields have been observed in twisted graphene structures. However, applying perpendicular displacement fields in nano-infrared experiments has only recently become possible [Li, H.; Nano Lett. 2020, 20, 3106-3112]. In this work, we fully characterize two approaches to realizing nano-optics compatible top gates: bilayer MoS2 and MLG. We perform nano-infrared imaging on both types of structures and evaluate their strengths and weaknesses. Our work paves the way for comprehensive near-field experiments of correlated phenomena and plasmonic effects in graphene-based heterostructures.

Twisted two-dimensional van der Waals (vdW) heterostructures have unlocked a new means for manipulating the properties of quantum materials. The resulting mesoscopic moiré superlattices are accessible to a wide variety of scanning probes. To date, spatially-resolved techniques have prioritized electronic structure visualization, with lattice response experiments only in their infancy. Here, we therefore investigate lattice dynamics in twisted layers of hexagonal boron nitride (hBN), formed by a minute twist angle between two hBN monolayers assembled on a graphite substrate. Nano-infrared (nano-IR) spectroscopy reveals systematic variations of the in-plane optical phonon frequencies amongst the triangular domains and domain walls in the hBN moiré superlattices. Our first-principles calculations unveil a local and stacking-dependent interaction with the underlying graphite, prompting symmetry-breaking between the otherwise identical neighboring moiré domains of twisted hBN.

Topological spin textures are field arrangements that cannot be continuously deformed to a fully polarized state. In particular, merons are topological textures characterized by half-integer topological charge ±1/2 and vortex-like swirling patterns at large distances. Merons have been studied previously in the context of cosmology, fluid dynamics, condensed matter physics and plasmonics. Here, we visualized optical spin angular momentum of phonon polaritons that resembles nanoscale meron spin textures. Phonon polaritons, hybrids of infrared photons and phonons in hexagonal boron nitride, were excited by circularly polarized light incident on a ring-shaped antenna and imaged using infrared near-field techniques. The polariton field reveals a half-integer topological charge determined by the handedness of the incident beam. Our phonon polaritonic platform opens up new pathways to create, control, and visualize topological textures.

Quasi-periodic moiré patterns and their effect on electronic properties of twisted bilayer graphene have been intensely studied. At small twist angle θ , due to atomic reconstruction, the moiré superlattice morphs into a network of narrow domain walls separating micron-scale AB and BA stacking regions. We use scanning probe photocurrent imaging to resolve nanoscale variations of the Seebeck coefficient occurring at these domain walls. The observed features become enhanced in a range of mid-infrared frequencies where the hexagonal boron nitride substrate is optically hyperbolic. Our results illustrate the capabilities of the nano-photocurrent technique for probing nanoscale electronic inhomogeneities in two-dimensional materials.