1 Introduction
Two-dimensional (2D) vdW layered materials [
1,
2], including graphene [
3,
4], silicene [
5], hexagonal boron nitride (hBN) [
6,
7], phosphorene [
8,
9], and transition metal dichalcogenides (TMDs) [
10-
12], have aroused wide attention in recent years due to their exotic properties [
11,
13] and promising applications in high-performance devices [
9,
14-
20]. The relatively weak interlayer vdW force enables vdW materials to be exfoliated into atomically thin layers with outstanding physical properties as well as highly adjustable 2D structures [
12]. A variety of methods [
21,
22] have been shown to efficiently regulate material properties, such as external field [
23-
26], doping [
26-
28], strain [
29-
33], integration with optical microcavity [
34], and changing dielectric environment [
35-
37]. In semiconducting 2D materials, the significantly enhanced Coulomb interactions due to reduced dielectric screening and increased quantum confinement give rise to rich excitonic species, which are formed by tightly bound electron and hole pairs and could persist up to room temperature, showing strong light−matter couplings and many-body effects [
38]. Specifically, TMDs possess fascinating optical and electronic properties, embracing widely tunable direct bandgap [
39], large oscillator strength [
40], substantial exciton binding energy [
41], and strong spin−orbit coupling that endows TMDs with spin-valley contrasting physics [
42-
44]. More importantly, different layered materials can be stacked together through vdW interactions to form various vertical heterostructures or homostructures, leading to enriched physical properties with high flexibility and tunability beyond the simple combination of the constituent monolayers [
45,
46]. Consequently, 2D vdW heterostructures offer an excellent platform for exploring modern solid-state physics, such as Bose−Einstein condensation (BEC) [
47] and fractional quantum Hall effect [
48]. The high adjustability of vdW heterostructures allows customized properties for the new generation of electronic and optoelectronic applications including quantum emitters [
49] and excitonic devices [
26,
50,
51].
The atomic alignment between the constituent layers of vdW heterostructure changes with interlayer lattice constant mismatch and relative orientation angle. Thus, the twist angle allows continuous modulation of strain distribution, interlayer coupling, electronic band structure and excitonic dynamics due to the incommensurate stacking [
38,
52,
53], which prompts a new topic “twistronics” [
54] that emphasizes the interplay between lattice rotational mismatch and interlayer electronic coupling. When the incommensurate heterostructures are stacked with a small twist angle or lattice mismatch, the atomic registries will change spatially with large-scale periodicity, forming so-called moiré patterns [
49,
55]. Such moiré patterns can be directly observed via nanoscale techniques [
56-
58] and attract abundant research attention due to the fascinating moiré effects in theoretical predictions. The superlattice structure of moiré pattern can generate strongly correlated phases, such as insulating states [
59] and unconventional superconductivities [
60] in magic-angle graphene bilayer, and Hofstadter’s butterfly pattern in graphene/hBN twisted heterostructure [
48,
61]. Recently, TMD vdW heterostructures with moiré superlattice show novel excitonic and optical physics including excitons localized by periodic moiré potentials [
62,
63], moiré exciton minibands [
64,
65] and spatially modulated optical selection rules [
49,
66,
67]. Furthermore, as a highly tunable system, moiré superlattice in vdW heterostructures can be feasibly combined with external factors, which will promote the exploration of frontier physics, including new spin optoelectronics, many-body interactions and potential applications in quantum information processing [
49].
Here, we focus on the interlayer-twist-modified electronic and optical physics in semiconducting vdW heterostructures and the burgeoning area of moiré superlattice. This review is organized as follows. In Section 2, we first summarize the optoelectronic properties affected by interlayer twist in vdW heterostructures, including the exciton properties, the electronic band structures and the exciton dynamics. In Section 3, moiré pattern and the emergent superlattice effects in vdW heterostructures are systematically explained. Specifically, we discuss the theoretical analysis and experimental observations of diverse moiré excitons, as well as how the properties of these excitons (localization and spin-valley configuration) are modulated by moiré superlattice. On this basis, the influence of external conditions such as electric field, doping, and strain on the moiré system is further discussed. In Section 4, we introduce the recently emerging field of correlated electronic states reported in TMDs-based moiré superlattices that possess several advantages compared with graphene-based twisted bilayers. The conclusions and future outlooks are summarized in Section 5.
2 Twist-angle-tunable optoelectronic properties in van der Waals heterostructure
Early experimental observations on vdW heterostructures have revealed that changing relative orientation angle or lattice mismatch could extend the optical and electronic properties of vdW bilayers. Most of these works focused on bilayers stacked in high-symmetry patterns (Fig.1), including A−A (point group D
3h), A−A′ (point group D
3d), A′−B (point group D
3d), A−B (point group C
3v) and A−B′ (point group D
3d) [
68]. Then the interlayer twist of bilayer structures is employed to realize continuous tunning, of which excitonic physics is the most concerned [
22,
41,
65]. It is theoretically predicted that relative twist between constituent layers can affect the excitonic phenomena and exciton dynamics through interlayer couplings [
52] and electronic band structures [
46,
68,
69], which has been confirmed and further explored in subsequent experiments [
45,
51,
70-
86].
2.1 Exciton properties controlled by twist-angle-dependent interlayer interactions
Interlayer coupling in multilayer system depends on the relative positions of atoms between layers. In 2H and 3R patterns [see Fig.1(b)], Mo and S atoms with opposite charges from different layers are stacked together, resulting in the strongest interlayer coupling. On the contrary, if atoms with the same charge in different layers are stacked together, the Coulomb repulsion between them increases the interlayer distance and reduces the coupling strength, among which the 3R-like pattern exhibits the weakest coupling due to most of the S atoms overlapped [
77]. Various theoretical calculations have reported the effects of stacking patterns on interlayer interactions. Early studies mainly focused on the effect of the interlayer distance on coupling strength under different atomic arrangements [
68,
69], which can induce tunable optical responses such as Raman [
75,
78,
87] and photoluminescence (PL) [
51,
77,
88] spectroscopy. Liu
et al. [
69] evaluated the geometries and stabilities of bilayer MoS
2 with five typical stacking patterns using density functional theory (DFT) calculations based on the local density approximation (LDA) function. They found that A−B and A−A′ in Fig.1(a) had shorter interlayer distances, lower relative energies and higher binding energies compared with the other patterns, as these two patterns exhibited the weakest interlayer repulsive interaction. Using random phase approximation (RPA), He
et al. [
68] studied the binding energy and interlayer spacing of TMD bilayers under different stacking patterns. They found that high symmetry A−A′ (bulklike 2H symmetry) was the most stable type with short interlayer distance, while the A−B stacking (3R type) followed closely behind, ~5 meV/formula unit less stable. The stacking order also affects the valence band splitting and the relative energy of the K point to Г point, where the exciton binding energy and corresponding optical response are influenced. In the case of incommensurate bilayer systems (with lattice constant mismatch and/or interlayer twist), numerical calculation for the estimation of interlayer coupling is nontrivial due to the lack of periodic features. For this, Wang
et al. [
52] adopted the analytical effective perturbative treatment to study the interlayer coupling represented by interlayer hopping integral in different valleys in both commensurate and incommensurate TMDCs bilayers. They show that the coupling strengths in ±K valleys are sensitive to the interlayer distance determined by the stacking pattern.
The stacking-dependent interlayer coupling is usually probed by Raman and PL spectra, which indicate the lattice vibration modes and optical transitions of the bilayer structures, respectively. Homojunctions of bilayer TMD with varying twist angles show remarkable modulations of optical transition energies and Raman modes [
75]. In stacked MoS
2 bilayer with CVD-grown monolayers, van der Zande and collaborators revealed the twist-angle-controlled strong tunning of the indirect optical transition energy and second-harmonic generation (SHG), as well as the weak tuning of direct optical transition energy and Raman frequency [
87], as shown in Fig.2(a). The shift of indirect optical transition energies was attributed to the angle-dependent interlayer separation due to repulsion between S atoms [see Fig.2(b)]. The results suggested that the effective interlayer interaction was maximized in the aligned (
θ = 0°) or anti-aligned (
θ = 60°) conformation, and minimized at
θ = 30°. Zheng
et al. [
74] reported the modulations of bandgap structures and phonon vibrations in a CVD-as-grown WS
2 bilayer. These bilayer WS
2 with random twist angles formed at high growth temperatures would exhibit a much stronger PL intensity of peak A and an additional peak A
1 compared with aligned stacking cases (
θ = 0° or 60°), accompanied with the absence of indirect transition peak I (interlayer exciton) [see Fig.2(e)]. Such modulations were attributed to two mechanisms. On one hand, a large interlayer distance weakened the coupling strength, which was proved by the broadening and redshift of the out-of-plane A
1g mode in Raman spectra. On the other hand, symmetry breaking in the random-twisted bilayers led to a higher possibility of interlayer hoping, which would be suppressed by spin-layer locking in aligned case. Zhang
et al. [
76] also observed an extra peak in PL spectra at 77 K in as-grown MoS
2/WS
2 heterostructure with type II band alignment, which was identified as the recombination of interlayer carriers [see Fig.2(c) and (d)]. Twisted heterostructures also show angle-tunable Raman modes. Wang
et al. [
78] reported a redshift of A
1g by ~2 cm
−1 in WSe
2/WS
2 bilayer compared with the case of monolayer WSe
2, along with the emergence of a new peak at 309.4 cm
−1 that was regarded as a general feature for strong coupling in bilayer system. The peak separations between A
1g(Г) and E
2g(Г) as well as the linewidths of A
1g(Г) in WS
2/MoS
2 heterostructures were periodically modulated by twist angles, caused by the interlayer coupling and charge transfer effects [
80]. Apart from TMD vdW heterostructures, Raman spectra of bilayer graphene are also affected by interlayer twist [
89,
90]. These twist-angle-sensitive Raman peaks demonstrate that the interlayer interactions and lattice dynamics can be regulated by interlayer rotation. Therefore, it is possible to control the electron-phonon coupling strength and thermal transport properties of 2D vdW heterostructures via twist angle [
91].
Since the optical responses of 2D semiconducting materials are mostly determined by their excitonic physics, lots of researches on modulations by stacking and twist focus on the excitonic properties [
88] and dynamics [
78,
79]. Huang
et al. [
88] concentrated on the trions and excitons under varying twist angles between 0° and 120° in the stacked bilayer of CVD-grown MoS
2. In PL spectra, both the ratio of the trion (
) to exciton (
) and the trion binding energy showed periodic and oscillatory features as a function of twist angle
θ, as shown in Fig.2(f). These quantities reached their maxima in aligned cases (
θ = 0° or 60°) and minima at
θ = 30° or 90°, consistent with the angle-dependent interlayer coupling strength. Using micro-absorption spectroscopy, Wang
et al. [
78] showed that in WSe
2/WS
2 bilayer after annealing, the intralayer WSe
2 A exciton presented significant PL quenching and spectral broadening due to angle-insensitive interlayer charge separation under strong coupling condition, while the interlayer excitons formed by ultrafast interlayer charge transfer exhibit sensitively twist-angle-dependent nature. The PL intensity of the interlayer exciton in MoSe
2/WSe
2 heterostructure was enhanced at coherently stacking (
θ = 0° and 60°) while vanished at incoherent intermediate angles. Utilizing vdW-corrected DFT calculations, Nayak
et al. [
70] found that the states at valence band edges were affected by the twist angle, which exhibited a long tail stretching over another layer in coherently stacked MoSe
2/WSe
2, but localized in an individual layer in incoherent case, as shown in Fig.2(g). The lattice mismatch controlled by twist angle would also significantly affect the properties of devices based on vdW heterostructures. Choi
et al. [
71] fabricated p-WSe
2/n-MoSe
2 junctions with
θ = 0°, 15° and 30°, and measured their
I−
V curves under illumination with different wavelengths. In the aligned case (
θ = 0°) where strong coupling is acquired, the sample showed a diodelike behavior under dark condition, functioning as a well-established p−n junction and displaying higher sensitivity under IR illumination. Whereas in the cases of lattice mismatch (
θ = 15°, 30°), interlayer potential barrier was built due to the vdW gap under weak coupling, and such barrier could be overcome under ultraviolet (UV) illumination, allowing photocurrent to flow in both forward and reverse directions. Apart from rotation mismatch, lattice relaxation and strain distribution of vdW twisted bilayers can also influence the formation of excitonic species and consequently the optical responses. Recently, Dai et al. [
92] have reported a four-fold enhancement of phonon-assisted photon up-conversion (UPC) in 5.5° twisted WSe
2 bilayer, along with PL intensity reduced by half, which was not observed in bilayers under other rotation angles such as 1.1° and 13.8°. Such UPC enhancement was attributed to the significant lattice relaxation at a mediate twist angle (5.5°) between the relaxed (1.1°) and rigid (13.8°) regimes. This mediate angle would cause increased strain and lattice deformation, leading to considerable interlayer separation. As a result, the formation of indirect interlayer exciton through non-radiative relaxation of direct neutral exciton was suppressed, which increased the UPC efficiency. On the other hand, the reduction of required assisted phonons due to the spectral redshift at the mediate angle also improved the UPC efficiency.
Light with circular polarization can be utilized to manipulate the valley degree of freedom in TMDCs, where the strong spin−orbit coupling (SOC) and the broken inversion symmetry lead to spin-valley-locked excitons at ±K valleys, resulting in the valley contrasting optical selection rule [
42-
44,
93-
96]. The spin-valley-locked excitons can be modulated by external electric [
23,
25,
97,
98] and magnetic fields [
99], allowing desirable potential applications [
16]. In bilayer systems, the valley-contrasting physics can also be regulated by the interlayer coupling and structural symmetry, depending on the stacking orders. Jiang
et al. [
77] fabricated MoS
2 bilayers with different stacking patterns by folding exfoliated monolayers to break the inversion symmetry and suppress interlayer hopping, resulting in strong valley and spin polarizations. Circularly polarized PL spectra showed the different degree of polarization (DOP) for different stacking bilayers, where the DOP in folded bilayers was comparable to that in monolayer, and much higher than that in natural 2H bilayer. The remarkable DOP of folded bilayers confirmed that the suppression of interlayer hopping and broken inversion symmetry generated spin-valley configuration locked to the layer index. On the contrary, in the centrosymmetric 2H-stacking bilayer, the interlayer hopping amplitude of holes was estimated to be comparable with the SOC strength [
100], giving rise to the mixing of layer pseudospins and suppressed valley- and spin-selective circular dichroism, as illustrated in Fig.2(h) and (i).
2.2 Electronic band structures and relative twist in momentum space
The variety of stacking patterns and relative twist angles between crystal orientations in 2D multilayer system provides a vast realm to modify the electronic band structures through angle- and mismatch-dependent vdW interactions at the atomic level, flourishing novel optoelectronic properties such as spatial- and/or momentum- indirect excitons with twist-angle-tunable nature [
72]. Unlike the monolayer TMD possessing a direct bandgap at the K point, bilayer TMD is a bulk-like indirect gap semiconductor. Liu
et al. [
69] used the LDA method to examine the band structures of bilayer MoS
2 with different stacking conformations. The indirect bandgap was found to be determined by the energy levels of K point in the conduction band and Г point in the valence band. The states at the K point in the conduction band were constitutive of the strongly localized d orbitals at the Mo atom sites (Mo-d), whose double degeneracy remained intact under weak interaction. Whereas, the states at Г point in the valence band originated from the linear combination of outer p orbitals of the S atom sites (S-p) and internal Mo-d orbitals, so the double level degeneracy at the Г point was lifted by the stronger interaction between the S-p from different layers, leading to a higher energy level at the Г point in valence band than that at the K point. Consequently, the band structure transformed from direct bandgap to indirect bandgap. Stronger interlayer coupling would enlarge the splitting at the Г point in the valence band and lead to a smaller indirect bandgap. Therefore, the direct A-exciton transition was insensitive to the interlayer interaction, while the indirect transition was sensitively dependent on stacking order [
77]. To gain a more accurate understanding of the stacking-induced modulation on electronic structures, excitonic effects and optical properties of bilayer TMDs, He
et al. [
68] utilized a many-body framework adopting the RPA and the GW approximation combined with the Bethe-Salpeter equation (BSE). The G
0W
0 quasiparticle band structure was represented by three valleys at T, K and Σ of the conduction band and two valleys at K and Г of the valence band, as shown in Fig.3(a). The conduction band minimum (CBM) was always located at the T point, while the valence band maximum (VBM) was related to the stacking pattern and vdW interlayer distance, which is mainly determined by the valence band splitting at K and the energy difference between K and Г. For A excitons, although their binding energy were sensitive to stacking patterns (about 60−100 meV higher for A−A′ than A−B pattern) the optical transition energies were not significantly dependent on stacking, in accordance with the experimental results [see Fig.2(a)] [
87]. Apart from that, the theoretically predicted stacking-induced strong tunning on Г valley in the valence band was also demonstrated by the observation of twist-angle-dependent indirect exciton transition.
Electronic structure and optical properties for twist heterostructures have also been well-explored. Jin and co-workers [
45] used cathode lens microscopy and microprobe angle-resolved photoemission spectroscopy (μ-ARPES) to realize a direct measurement of the electronic structure in twisted graphene/MoS
2 vdW heterostructures. The observed band structures indicated that the electronic structure of graphene (Gr) remained intact, while the MoS
2-derived band exhibited twist angle dependence, as shown in Fig.3(b) and (c). In particular, the relative energy between Г and K decreased from ~0.2 to ~0 eV as the twist angle varied from 5° to 28°, suggesting an indirect bandgap character except for twist angle near 30°. These results showed a similar trend to the DFT calculation prediction, corroborating the view that the electron-transfer-induced strain between graphene and MoS
2 would affect the Mo−S bond and lead to the twist-angle-tunable band structure [
101]. In TMD vdW heterostructures, the band structures exhibit direct- or indirect-gap determined by twist angles through tuning the relative energy levels between K and Г but have robust type II band alignments [see Fig.4(a)]. For twisted MoSe
2/WSe
2, first-principle DFT calculations reveal a tunable band structure shown in Fig.3(f) [
70], where the direct K−K transitions and indirect K−I transitions are almost degenerate, and the conduction band at Г point varying with twist angle presents a trend similar to that in homostructure TMDs [
68,
69]. However, the observed VBM is located at K instead of Г due to the large SOC. The VBM and CBM at the K-point mainly come from WSe
2 and MoSe
2, respectively. Since direct transition possesses a larger amplitude than that of indirect transition which involves with second-order process, the interlayer transitions observed experimentally are expected to occur between band edge states from two different layers at K point with twist-angle-dependence [Fig.2(g)]. For example, in twisted WSe
2/WS
2 heterostructure, VBM mostly originates from orbitals of the WSe
2 layer, while CBM is mainly contributed by the WS
2 layer. Such strong type-II band alignments are independent of the twist angle and bandgap nature [
78]. The spatially indirect nature of interlayer exciton makes its transition dipole smaller than that of the intralayer exciton by one to two magnitudes and consequently weaker oscillator strength [
54], which has been theoretically predicted [
46,
65] and experimentally demonstrated [
18,
78].
Such type II band alignment can easily promote the formation of interlayer exciton in which electrons and holes reside at K valleys from different layers. The electron and hole valleys of interlayer excitons will shift away from each other under a small twist and/or lattice mismatch, leading to a momentum mismatch of the Brillouin zone, shown in Fig.3(h). The exciton’s kinematic momentum is defined as
, where
and
represent kinematic momentum of carriers from different layers, as shown in Fig.3(i). The transition dipole of interlayer excitons is sensitive to
, forming sixfold degenerate light cones at finite center-of-mass (COM) velocities on the parabolic energy dispersion [
46]. The so-called “light cone” means the region for allowed interconversion with photons in the exciton phase-space. When the electrons and holes relax to the bottom of the ±
K valleys in momentum space with
Q = 0, their transitions must involve phonon or impurity scattering to enter the finite-velocity light cones, resulting in long lifetimes of the interlayer excitons due to this inefficiency scattering process [
102]. However, when the small mismatch can be compensated by the exciton’s kinematic momentum
, direct transition for radiative recombination is allowed. Here,
represents the displacement from the Brillouin zone corner
of one constituent layer to the nearest Brillouin zone corner
of the other layer, where
denotes the electron-hole valley configuration. Such interlayer excitons are within the light cones, whose transition dipole strength would decay rapidly with an increase in the magnitude of momentum. Only in commensurate stacking at
θ = 0° or 60° with negligible lattice mismatch, all the main and Umklapp light cones can merge into a single cone at
Q = 0 and get close enough to the bottom of energy dispersion, leading to strong interlayer exciton PL [
46,
103].
Unlike the interlayer excitons with spatially separated electron and hole both located at K valleys, Kunstmann and co-workers [
72] have identified a new type of interlayer excitons in MoS
2/WSe
2 heterostructure, where the hole and the electron reside at the hybridized Г valley and localized K valley, respectively, as shown in Fig.3(d) and (e). This momentum-space indirect interlayer exciton was revealed by combining temperature-dependent PL spectroscopy and first-principle calculations. The emission energy of this kind of interlayer excitons could be tuned via relative twist by controlling the hybridization [
72]. Likewise, Sung
et al. [
82] have also observed and investigated the Г−K interlayer excitons at temperature
T = 4 K. These momentum-indirect interlayer excitons possessed opposite out-of-plane electric dipole moments supported by broken mirror and inversion symmetry due to the spatially alternating stacking domains, which originated from atomic reconstruction. The electric dipole moments of Г−K interlayer excitons were identified by the linear Stark effect in PL spectra, measured in a hexagonal boron nitride (hBN) encapsulated near-0°-twisted MoSe
2 bilayer under an out-of-plane electric field, as shown in Fig.3(g). Additionally, K−K intralayer excitons with higher energy were also observed, which exhibited field-asymmetric hybridization with interlayer K−K excitons under large fields, manifesting as field-dependent avoided crossings between the two corresponding peaks in reflectance spectra [Fig.3(g)]. The electric field dependence of these K−K interlayer excitons also varied spatially under alternating crystal symmetry. This result indicates the profound effects of crystal symmetry on excitonic properties and potentials for realizing alternating exciton arrays with high tunability.
In TMD heterostructures, interlayer excitons with spin-valley polarization inherited from the monolayers can directly interconvert with a photon with elliptic polarization [
46,
54]. The interlayer exciton spin-valley polarization originates from valley-polarized intralayer excitons through interlayer charge transfer. Such valley-contrasting physics can be modified by interlayer twist, which can expand fascinating and complicated properties through changing mismatch between the Brillouin zone of constituent layers in momentum space. In exciton momentum space, light cones of interlayer exciton possess helicity determined by valley configurations of electron and hole, which allow tunable resonant optical injection of valley excitons and valley currents under linearly polarized light excitation, as illustrated in Fig.5(a) [
46]. For the commensurate case (i.e.,
θ = 0° or 60° with neglectable lattice mismatch ~0.1%) where all light cones merge into a single one at
Q = 0 and interlayer exciton PL gets enhanced, the brightness and the polarization selection rules depend on the interlayer translation. This interlayer exciton system allows observation of valley Hall and inverse valley Hall effects, which is regarded as an ideal platform to investigate Berry phase effects in the Bloch band [
46]. Experimentally, Rivera and collaborators [
51] observed interlayer excitons with valley-contrasting physics in WSe
2/MoSe
2 vertical heterostructures at near-0° stacking with slight twist and lattice mismatch at room temperature. In this case, the existence of light cones at small kinematic momenta allows interlayer excitons to interconvert with photons and radiatively recombine through scattering such as exciton‒phonon or exciton‒exciton interactions. Interlayer hopping between valleys with the same electron spin from two constituent layers becomes the dominant relaxation channel under small momentum mismatch, leading to the co-polarized optical selection rules [see Fig.5(b) and (c)]. Whereas, the interlayer hopping between spin-contrasting valleys is strongly suppressed under large momentum difference. Interlayer exciton exhibits a long valley lifetime of about 40 nanoseconds, a few orders of magnitude longer than that of intralayer excitons in monolayer, enabling the visualization of the valley-polarized interlayer exciton cloud expanding, including lateral drift and diffusion over several micrometers [
51]. Therefore, twisted van der Waals heterostructures act as a promising playground for exploring valley exciton physics. In recent years, several works have reported on the chiral optical properties and valley polarization physics in vdW bilayers via twist modulation. Scuri
et al. [
73] controlled the spin-valley properties of interlayer excitons in twisted WSe
2 bilayers encapsulated in hBN by tunning the momentum separation between valleys of different layers [Fig.5(d)]. The interlayer excitons presented long valley lifetimes (> 40 ns) due to the indirect nature in both real and momentum space [
65,
103-
106]. Michl
et al. [
83] reported a high degree of intrinsic circular polarization of exciton PL driven by non-resonant linearly polarized laser in nonmagnetic twisted hBN/MoSe
2/WSe
2 heterostructure (
θ = 5.3 ± 1.4°) at 1.8 K [Fig.5(e)], which was attributed to the kinetic spin polarization of excitons arising from new pyromagnetic symmetry of chiral stack structures. Furthermore, it has been predicted theoretically that valley splitting of WSe
2 would be significantly enhanced in twisted WSe
2/CrI
3 heterostructure in contrast to non-twisted heterostructure, which was equivalent to the case with a perpendicular external magnetic field of about 20 T [
107].
2.3 Twist-angle-dependent exciton dynamics
Lattice mismatch and relative twist between layers in vdW heterostructures with the ability to modify the interlayer interaction and electronic structure efficiently, are supposed to play a key role in exciton dynamics, including ultrafast charge transfer, diffusion, and relaxation processes. With the help of transient spectroscopy such as pump-probe technique, time-resolved SHG and ultrafast electron diffraction (UED), a large range of works have focused on how exciton dynamics in heterostructures are affected by twist angle over these years [
52,
78,
79,
84,
85]. Physical properties of semiconductor heterostructures are highly influenced by band alignment, which can be classified into three types, including type-I (straddling gap), type-II (staggered gap) and type-III (broken gap) [
108], as depicted in Fig.4(a). Theoretical calculations have found that the type-II interface band alignment is the most typical and common case for TMD vdW heterostructures [
13,
109]. Meanwhile, the type-II band alignment with CBM and VBM of the heterostructures located in different layers supports the electron and hole to transfer to the spatial-separated band edges [see Fig.4(b)].
Early experiments have demonstrated the rapid interlayer charge transfer (within 100 fs) by various methods such as micro-absorption spectroscopy, PL mapping and femtosecond pump-probe spectroscopy [
78,
110]. Such efficient interfacial charge transfer (< 50 fs) dynamics in TMD vdW heterostructures do not exhibit twist-angle dependence [
79,
111-
113]. Ji
et al. [
79] found the charge transfer in MoS
2/WS
2 heterostructures was robust in a time scale of about 90 fs under varying twist angles with different interlayer coupling strengths, measured by optical two-color ultrafast pump-probe spectroscopy, as shown in Fig.4(c) and (d). The twist-angle-independent charge transfer with high efficiency violates the common perception that the interlayer charge transfer dynamics would sensitively depend on the relatively weak interlayer coupling affected by twist angle [
54]. Moreover, the momentum mismatch between valleys from different constituent monolayers is expected to reduce the efficiency of interlayer charge transfer. Utilizing the STEM and time-dependent DFT simulations, the atomic-level local structural inhomogeneity caused by interlayer stretching and shifting is considered to be the origin of the observed robustness of ultrafast interlayer charge transfer. Such interlayer sliding provides extra parallel multi-channels for the efficient transfer process [
79]. From the perspective of momentum mismatch in twisted heterostructures, phonon scattering coupled with interlayer hybridization mediates the twist-insensitive ultrafast interlayer charge transfer [
52,
54]. Unlike the layer-localized ±
K valleys whose coupling strengths exhibit sensitive dependence on the interlayer twist and translations, the
Гv (
Г point of valence band as the VBM) and
Qc (six CBMs near the middle of the
Г−
path) show strong interlayer coupling comparable to the corresponding band offsets in TMD heterostructures, as shown in Fig.4(f). Such strong interlayer mixing is attributed to the nonignorable effects of
orbital of chalcogen atoms, as well as the relatively larger hopping integral with smaller momentums compared with those in ±
K valleys [
52]. These highly layer-hybridized valleys allow the twist-independent ultrafast interlayer charge transfer of type-II TMD heterobilayers via a two-step mechanism involving the
Гv and
Qc valleys as intermediary states. Specifically, the photoexcited electron with high energy in
valley can be scattered to the
Qc valleys through coupling with phonons, defects, or other carriers, and subsequently relax to the
, as depicted in Fig.4(g). The high energy hole in
can similarly transfer to the
mediated by the strongly layer mixed
Гv valley. The efficient scatterings are supported by the strong couplings between carriers in ±
K valleys and the related phonons. Wang and co-workers [
52] predicted this two-step interlayer electron (hole) transfer process assisted by phonons, and the calculated intervalley scattering rates accorded with experimental results. Twist in bilayers does not influence interlayer coupling and band offset of
Гv valley, while
Qc valleys are always within the strongly interlayer hybridized ring region under arbitrary twist angle. Therefore, the ultrafast (< 50 fs) interlayer charge transfer in TMD heterostructures is unimpeded by interlayer rotational mismatch [
52]. This two-step mechanism has been demonstrated by the observation of ultrafast electron intervalley scattering in MoS
2/WS
2 heterostructure with a small twist angle of 5° using the time- and angle-resolved photoemission spectroscopy (TR-ARPES), where electrons in
K valleys are scattered to intermediary valleys (M, M/2 and Q) within 70 fs and relaxed to
and
valleys on timescales of ~0.4 ps [
114].
However, such robustness against interlayer twist in TMD heterostructures does not persist in all circumstances. Recent experimental researches on how the twist angle affects rapid charge transfer show diverse results. For a longer timescale of hundreds of femtoseconds in twisted WS
2/WSe
2 heterostructures, the interlayer tunneling rate decreases significantly as the rotational mismatch rises, affecting the formation and relaxation dynamics of interlayer excitons [
115]. On the other hand, a more complicated stacking angle modulation mechanism of ultrafast dynamics has been revealed by more precise measurements (e.g., much higher time resolutions). Zimmermann
et al. [
85] investigated the ultrafast charge transfer dynamics affected dramatically by stacking angle in MoS
2/WSe
2 heterostructures at an ultrafast time resolution of about 6 fs. They conducted the measurement by the time- and polarization-resolved second-harmonic (SH) imaging microscopy, which could directly identify signals from different regions (two kinds of monolayers and their overlap) and control the SH intensity of different regions. For 1.70 eV pump which only resonantly excited the intralayer A-excitons of WSe
2, the extracted transfer time Δ
tCT from WSe
2 to MoS
2 layer of probed electron sensitively depended on stacking angle: which are comparable in quasi-3R (
θ = 9°, Δ
tCT = 65 ± 10 fs) and misaligned (
θ = 16°, Δ
tCT = 84 ± 10 fs) cases, while become much faster in quasi-2H (
θ = 52°, Δ
tCT = 12 ± 6 fs) case [see Fig.4(e)]. The twist-angle dependence was attributed to the different wave function overlap at
K valleys between constituent layers at different stackings. Such distinct wave function overlap also influenced the recombination lifetimes, about hundreds of picoseconds, where stronger spatial coincidence was found in quasi-2H case [
85]. When both materials were excited, the transient SH response of quasi-2H heterostructure presented a rapid recovery on femtosecond timescale slightly slower than those in the corresponding monolayer, which was not observed in the quasi-3R and misaligned samples. This discrepancy was explained by the competition between the coherent (intralayer) radiative recombination and additional relaxation tunnels, where holes from the MoS
2 layer could be scattered to the WSe
2 layer in
K valleys or strongly hybridized states in
Г valley. As shown in Fig.4(e), hole transfer from
in MoS
2 to
in WSe
2 was spin forbidden in the quasi-2H case, and the misaligned case had a larger momentum mismatch than that in quasi-3R stacking sample. Therefore, the efficiency of the hole transfer decreased from quasi-3R, misaligned to quasi-2H, which corresponded to the increasing coherent recombination [
85]. This result indicates the possibility of modifying the layer mixed states through stacking configuration, thus affecting the exciton dynamics. Apart from the vdW heterostructures composed of the TMD family, the charge transfer and relaxation dynamics of TMD-Gr heterostructures are also strongly modulated through the twist-angle-tunable overlap between interlayer electronic states. Recently, Luo and co-workers [
86] have realized the simultaneous visualization of charge transfer and electron−phonon coupling using the UED method and demonstrated the sensitive dependence of interlayer charge transfer on the twist angle [
116].
In vdW heterostructures, intralayer and interlayer Coulomb interaction leads to the formation of diverse quasi-particle species such as intralayer excitons, interlayer excitons and trions. Twist angles in vdW heterostructures also influence the optical behaviors of these “many-body” species. Recently, Shi
et al. [
84] observed increased luminescence ratio of positive trions compared with excitons in steady-state PL spectra of WS
2/WSe
2 heterostructure under twist angle of 30° and 60° but not 0°. They used ultrafast transient absorption spectroscopy to reveal the faster formation time for trions than intralayer excitons at
θ = 60°, and the higher formation probability of trions due to accumulation of residual free carriers at
θ = 30°, as illustrated in Fig.4(h).
There are several theoretical discussions on moiré physics and a considerable part of experimental observations of the moiré quasi-periodic structures using high-resolution methods such as annular dark-field scanning transmission electron microscopy (ADF-STEM) [
56]. Nevertheless, comprehensive investigations focusing on moiré pattern remained scarce in early experimental researches. Since the moiré pattern is easily disordered in a heterostructure without hexagonal boron nitride encapsulation [
93], limitations of the quality of the sample mainly account for this lack. Besides, the effects arising from moiré potential would be submerged at room temperature.
3 Moiré excitonic physics
Twistronics investigates the physics including electronic, optical and magnetic properties [
117] of twisted 2D van der Waals heterostructures, and has promoted various studies as mentioned above. In the wake of experimental technology progress, an increasing number of researchers have set their sights on the moiré superlattices (MSL), which are formed by heterostructures stacked with a finite relative twist angle or lattice mismatch. MSL can trigger various exotic properties, such as flat bands [
49,
50,
64,
118-
120], excitons modulated by moiré periodic potential [
57,
121] and many-body effects [
122] in such artificial superlattice. The magic angles in twisted bilayer graphene generate fascinating properties including correlated insulator states and unconventional superconductivity originating from the flat bands [
59-
61,
123]. Since then correlated electronic physics has been investigated in various graphene-based structures such as ABC trilayer graphene/BN [
61,
124-
126] and twisted double bilayer graphene (TDBG) [
127]. In the last few years, moiré effects and correlated phenomena have also been found in other 2D materials-based systems such as TMDs [
128-
131]. TMDs possess unique advantages compared with graphene due to their desirable properties, such as sizeable bandgap and large spin-orbit coupling. Therefore, MSL in TMDs would perform thermal stability and feasibility for topological band engineering and spatially alternating spin-valley configuration [
93]. There are several theoretical predictions on the moiré modulated excitons and corresponding optical response of TMDs-based heterostructures, like localized excitons and nanodot quantum emission arrays [
132]. In addition, TMD MSLs have flat bands that exist over a wide range of stacking angles instead of certain magic angles [
133], making MSL easier to realize and providing an extra rotation degree of freedom for modulating the moiré effects.
When vertically stacked vdW layers present a complex periodic pattern consisting of varying interlayer atomic alignment, the in-plane moiré superlattice will take shape, as shown in Fig.6(a) and (b). Moiré superlattice will be formed in two ways mainly determined by the magnitude of relative twist, including the rigid-lattice and atomic reconstruction [
132]. The rigid-lattice moiré pattern would naturally occur under a relative rotation angle or lattice mismatch between constituent layers, and such misalignment serves as a crucial node for controlling the electronic structure and optical properties [
134-
137]. The periodic MSL can be represented by a moiré wavelength
λm [
138]:
where lattice mismatch is presented as
, and
,
denote the lattice constants of the two constituent layers, with
θ as the relative twist angle. The corresponding moiré period (also known as moiré lattice constant) under small lattice constant mismatch δ can be approximately represented as
[
65], with length scales from a few nanometers to tens of nanometers [
54,
93]. As mentioned above, there are different types of commensurate stacking patterns, which can be subsumed under two categories: near-rhombohedral stacking (R-stacking) and near-hexagonal stacking (H-stacking). The high-symmetry local atomic registries existing in one moiré supercell preserve three-fold rotational symmetry
and generate charming properties such as spatially alternating optical selection rules [
49]. As depicted in Fig.6(a) and (b), these high symmetry points can be labeled by the notation
(A = h, X, M) for R-(H-)stacking, indicating the “A” in the electron layer aligned to the “h” in the hole layer. The superscript “A” here refers to the lattice center (h), chalcogen atom (X) or metal atom (M). Varying interlayer atomic registries in the moiré supercell result in the spatial modulation of interlayer distance, exciton potential energy and spin-valley characteristics [
49,
93,
132]. For stacking order at near-commensurate twist angle (
θ ≈ 0° or 60°), atomic-level reconstruction would be induced through interlayer coupling and lattice relaxation. Instead of a smoothly varying rigid-lattice moiré pattern, discrete commensurate regions with narrow boundaries will be formed [
139-
142]. These discrete domains with constant conductivity could be observed using conductive atomic force microscopy (CAFM), as shown in Fig.6(c). Furthermore, the existence of strain variation across the moiré supercell leads to substantial strain accumulations at region boundaries and domains with relatively unstable stacking orders, which introduces lattice deformation and extra modulation on band structures [
143]. DFT calculated band structures revealed the spatial modulation of the relative electronic band edges induced by atomic reconstruction [
140]. Consequently, such spatial distribution of considerable strain affects the formation of excitonic species and optical properties, which are also highly tunable by relative twist [
92].
As the fabrication, detection and observation of moiré structures become increasingly achievable, multiple experimental studies on moiré excitons have been performed through electronic and optical measurements such as STEM [
56], PL [
63,
144,
145] and Reflectance spectroscopy [
145,
146]. The mechanisms of how moiré structure modulates intralayer and interlayer excitons have attracted much attention but still lack a comprehensive explanation. For the moiré superlattice with a periodicity of about 10 nm in a twisted semiconducting vdW heterostructure, excitons have the Bohr radius between the moiré wavelength and the lattice constant of monolayer [
132]. Therefore, exciton in superlattice can be regarded as a composite quasi-particle modulated by exciton potential with spatially periodical variation, manifesting in the locally changing oscillator strength and optical transition dipole [
49]. The center-of-mass motion could be described by [
147]:
where
is a constant representing the energy of the lowest exciton state without moiré potential energy,
is the center-of-mass kinetic energy, and
denotes the exciton moiré potential energy with
as the exciton effective mass.
varies smoothly and can be approximated by the lowest-order harmonic expansion [
64,
65]:
where
defines the reciprocal lattice vectors of the moiré superlattice under a finite
θ [
65]. Near
points,
has the form of a harmonic oscillator with quantized energy levels modulated by twist angle through changing moiré period [
147]. This harmonic-potential model is widely accepted in most TMD heterostructures but cannot dominate all systems [
132]. In the case of band alignment under small band-edge offset, the near-resonant hybridization between interlayer and intralayer excitons would occur through interlayer tunneling in moiré minibands, which generates hybridized excitons with intriguing features [
146,
148-
150].
3.1 Theoretical predictions and experimental observations of new excitons in moiré structure
As mentioned in the previous section, early experimental methods and sample quality cannot provide sufficient conditions for the observation and further investigation of moiré modulated exciton properties. Whereas, several theoretical prediction works [
49,
52,
64-
66,
118] on the new physics of excitons in moiré superlattice have been reported. In the superlattice with a large-scale moiré period, a mini Brillouin zone (BZ) would be formed, whose corners are located at
−
and its π/3 rotations [
52], as depicted in Fig.7(a). The original monolayer bands are folded into a series of closely spaced mini bands inside the mini BZ, where interlayer hopping can happen between states with the same wave vector
k from different layers [Fig.7(b)]. The local atomic registry is spatially varying in a moiré supercell much larger than the monolayer unit cell, resulting in a periodic modulation of the local band structure [
52]. Yu
et al. [
49] described the moiré modulated exciton as an exciton wave packet moving adiabatically in the periodic potential, which was constructed by the basis of kinematic momentum eigenstates under misaligned heterobilayer. In the near R-type stacking case, exciton wave packets at the three high symmetry points [
−
,
−
and
−
in Fig.7(c)] have distinct
transformations due to the different
eigenvalues of the electron Bloch function:
Here, χ is the exciton wave packets and
denotes the spin-valley index for exciton at
K (−
K) valley with up (down) spin. As shown in Fig.7(d) and (e), the convertible photons must possess the same rotational symmetry, which leads to the spatially alternating optical selection rules. The interlayer exciton emissions are circularly polarized at A and B points with opposite helicity and forbidden at C, while the emissions are elliptically polarized at other locals in the moiré supercell. The distinct rotational symmetries in the three points also make them the extrema of exciton potential [Fig.7(e)]. The strong deep traps of exciton potential (~100 meV) around the global minima at A points can realize perfect arrays of nanodot confinement, which can act as uniform quantum emitters tunable by perpendicular electric field and strain [
49]. The moiré pattern with periodic varying interlayer separation was observed experimentally by scanning tunneling microscope (STM) measurement in MoS
2/WSe
2 heterobilayer, as shown in Fig.7(f) [
56]. The STM and scanning tunneling spectra (STS) also revealed the formation of an electronic superlattice and detailed information on the local electronic structures affected by interlayer coupling. The experimentally observed and DFT-calculated site-dependent electronic structures showed consistent trends as illustrated in Fig.7(g). Wu and co-workers [
65] adopted the theory of moiré potential modulated interlayer exciton and predicted that moiré potential affected interlayer excitons in two aspects: First, the periodic potential mixed and split the three independent exciton COM eigenstates that were otherwise degenerate at the moiré Brillouin zone (MBZ) corners, resulting in new exciton eigenstates with rotational symmetries of the moiré pattern. Thus, the circular optical selection rule remained, which was sensitive to the spatially varying potential energy and local optical absorption strength, but not locked to the valley. The moiré potential energy mixed COM eigenstates, and consequently split the spectral weight associated with interlayer excitons across a variety of different states. Such energy splitting can be tuned by twist angles through changing moiré periodicity and the MBZ [Fig.7(h)]. Second, the exciton states could be confined in potential minima [Fig.7(i)], allowing the realization of nanoscale 2D arrays of identical quantum dots. The calculated exciton band structure in the moiré Brillouin zone had flat-band features for the first few exciton bands, showing the mixing of excitons with different momenta and strong confinement for these Wannier excitons states. The exciton moiré band structure would exhibit Dirac dispersion within the light cone when the moiré potential was tuned to have an additional two-fold rotational symmetry besides the
symmetry [
65]. Wu and co-workers also focused on the intralayer excitons under periodic moiré potentials and found that moiré superlattice in the TMD bilayer generated the twist-angle-dependent satellite excitonic peaks in the optical absorption spectrum [
64]. These theoretical works also point out the topological features of the moiré exciton bands that would further relate to the unique Hubbard models with twist-angle-dependent hopping and interaction strength [
118]. Twisted bilayers with moiré structure provide an accessible platform to simulate the Hubbard models for exploring many-body phenomena such as spin-liquid states and quantum anomalous Hall insulators [
49,
64].
With the development of experimental methods, multiple observations of moiré effects in vdW heterostructures have been reported and showed diverse characteristics and fascinating physics. In most cases, the periodic moiré potential was regarded to play a crucial role in the emergence and modulation of excitons in a twisted heterostructure, supporting the harmonic-potential model based on perturbation theory. Pan
et al. [
57] observed the moiré pattern in MoS
2/WSe
2 on epitaxial graphene using STM and STS at temperatures of 5 and 80 K. STS showed significant energy shifts of about 0.2 eV between the maxima and minima of the moiré corrugation [Fig.8(a)]. Sharp peaks (linewidth <10 meV) occurred near band-edge energies of the spectra for both the valence band (VB) and the conduction band (CB) at 5 K, which were attributed to the quantum-confined states associated with the Г
W−VB and K
Mo−CB under moiré potential. These confined states showed obvious temperature dependence and the corresponding narrow peaks vanished at 80 K. The spatially periodic in-plane potential leads to splittings of both intralayer exciton and trion energies, which was observed in the emission and absorption spectra of MoSe
2/MoS
2 heterostructure encapsulated in hBN at 5 K, as shown in Fig.8(b)‒(d) [
121]. The energy differences between the split peaks were consistent with the theoretical predictions mentioned above. The intensities of the split peaks of excitons and trions also exhibited temperature dependence, suggesting that excitons can be detrapped from moiré potential minima through thermally activated phonon scattering. Moiré exciton states were also observed from the multiple emergent peaks around the original WSe
2 A exciton resonance in the absorption spectra of closely aligned WSe
2/WS
2 heterostructures at 10 K [
151]. These peaks were considered to arise from the strongly coupled heterobilayer rather than from several separated domains, since the photoluminescence excitation (PLE) spectrum showed the strong enhancement of the interlayer exciton emission under excitation at the energies of each of the three emergent peaks (labeled as resonances I, II and III), as shown in Fig.8(e). Different moiré exciton peaks in WSe
2 exhibited distinct dependences on electron doping in WS
2 [Fig.8(g)]. The discrepancy indicates that peaks I, II and III arose from exciton states with different real-space distributions in the moiré superlattice. These phenomena were well explained by the theoretical model in which the periodic moiré potential was much stronger than the exciton kinetic energy and generated multiple flat exciton minibands [
151]. Besides, excitons modulated by moiré potential have also been observed in TMD vdW homostructures [
63,
141,
145,
152]. Wu
et al. [
63] observed moiré intralayer excitons in WSe
2/WSe
2 bilayers with small twist angles in the PL spectrum at 8 K. The extracted moiré potential depth (inversely proportional to the moiré superlattice period) increased with the twist angle, manifesting as more splitting and broader energy peaks [Fig.8(h)]. As the excitation power increased, moiré excitons turned into a main peak dominated by intralayer excitons with a larger peak width in the PL spectrum. The power dependence of moiré exciton peaks was attributed to the saturation of flat minibands filled by trapped excitons [
63]. Moiré excitons were also identified in the twisted WS
2/WS
2 homostructure with a twist angle of 1.46° ± 0.07° determined by SHG [
145]. The multiple emergent peaks in the PL spectrum at 8 K exhibited a nearly constant peak spacing of 19 ± 3 meV and showed similar pump power dependence, which also resulted from the flat electronic minibands generated by the reconstructed moiré potential.
Although the effective harmonic-potential model has been utilized to describe moiré effects successfully in TMD heterostructures with considerable interlayer band edges offsets like WS
2/WSe
2 and MoSe
2/MoS
2, this theoretical model cannot explain moiré excitonic behaviors in all systems. In the case of TMD heterostructures with a slight band-edges offset featuring nearly degenerate carrier bands, such as MoSe
2/WS
2 and some twisted homobilayers, resonant interlayer hybridization must be explicitly considered. The interlayer hybridization would enhance the moiré superlattice effects on the electronic structures and lead to the formation of hybridized excitons (hXs) through interlayer carrier tunneling [
146,
150]. The hybridized excitons arise out of the band hybridization effects [Fig.8(j)], consisting of strongly mixed electron or hole states involved in the formation of intralayer excitons (Xs) and interlayer excitons (IXs). The hXs inherit the brightness of intralayer excitons and the polarity of interlayer excitons, leading to their identical optical selection rules to those of intralayer excitons in constituent monolayer and sensitivity to vertical electrical bias, respectively [
132,
150]. The theoretical absorption spectra of resonant heterobilayers MoSe
2/WS
2 varied sharply with twist angle, as shown in Fig.8(l), in which a bright line at higher energies appeared when closely aligned due to moiré umklapp electron-photon interactions. The umklapp processes fold exciton states with finite-momentum onto zero-momentum, regarded as the evidence for moiré superlattice minibands for the excitons. Such features were absent in non-resonant cases due to the weak mixing between IXs and X [
150]. Alexeev
et al. [
146] investigated the interplay between moiré superlattice effects and twist-angle dependent hybridization of X and IX in MoSe
2/WS
2 heterostructures with near-degenerate conduction band edges. The exciton energy and linewidth in PL spectra change periodically with the interlayer twist angle at a period of
θ ≈ 60°, suggesting a connection to the valley rotational mismatch between the Brillouin zones of two constituent monolayers
(
). The conduction-band states had plane-wave projections onto all Bragg vectors of their corresponding reciprocal lattices
and
. As a result, the interlayer tunneling between states satisfied momentum conservation
. Here,
and
were the corresponding first-Brillouin-zone wave vectors from the K valley. The
denoted the first Bragg vectors of moiré superlattice, which defined the twist-angle-dependent mini BZ, as depicted in Fig.8(i). In closely aligned (or anti-aligned) MoSe
2/WS
2 heterobilayers, the minimized momentum mismatch for interlayer conduction-band tunneling enabled hybridization between the interlayer excitons and the MoSe
2 intralayer excitons (both A and B excitons), leading to shared hole band of the hybridized excitons (hXs) [Fig.8(i)] [
146]. The resulting hXs (labeled as hX
1) and higher-energy peaks arising from the moiré mini-band states of hX (labeled as hX
2 and hX
3) appeared at closely aligned case in reflectance contrast spectra at low temperatures (
T = 10 − 105 K). Apart from the bright excitons, momentum-dark excitons in moiré superlattice have also been reported to be affected by interlayer hybridization. Berm
et al. [
149] investigated the twisted WSe
2 bilayers and showed the twist-angle-dependent electron hybridization at the Λ point of momentum-indirect K−Λ exciton in phonon-assisted PL (
T = 4 K). The K−Λ excitons exhibited a prominent hybridization effect due to their strong delocalization over both layers, in contrast to the K-K excitons, as shown in Fig.8(f) and (k).
Recently, several experimental methods apart from optical measurements for observation of moiré superlattice and moiré exciton have been reported. Andersen and co-workers [
141] realized the direct stacking domains imaging in vdW heterostructure devices. They developed a secondary electron microscope (SEM) technique to investigate reconstructed moiré patterns in twisted WSe
2/WSe
2 bilayers. They directly demonstrated the considerable variation of local moiré periodicity in a small twist angle bilayer and correlated the local structure with distinct optical properties. Karni
et al. [
153] utilized the time-resolved and angle-resolved photoemission spectroscopy (TR-μ-ARPES) to capture images of the electron and hole distributions of interlayer exciton in WSe
2/MoS
2 heterostructure. The interlayer exciton was found to have a diameter comparable with the moiré period, whereas its center of mass was highly confined in a region much smaller than the size of the exciton itself.
3.2 Localization and modulation of excitons by moiré superlattice
Excitons in TMD moiré superlattices under the spatially varying potentials show in-plane localization [
153] and possess the flat-band feature. The flat bands are related to the in-plane localization due to the vanishing of group velocity [
154]. The moiré potential traps are analogous to the localized defect levels in semiconductors, where the moiré periodicity corresponds to the defect−defect distance. The “defect” bands would become less dispersive as the periodicity increased and their interactions reduced [
50]. Since the amplitude and period of moiré potential are sensitive to the interlayer stacking and twisting in vdW heterostructures, the spatial distributions, resonant energies and oscillation strengths of moiré excitons can also be continuously tunned by twist angle. Guo
et al. [
50] theoretically predicted the distributions of moiré excitons in twisted MoS
2/WSe
2 heterostructures under different angles. Based on time-dependent density functional theory (TDDFT), they predicted that exciton with lower energy could be more localized and would experience deeper moiré potential under a small twist angle relative to
θ = 0° or 60°, as shown in Fig.9(a). The holes showed stronger in-plane localization than electrons, consistent with the narrower bandwidth of VBM than that of CBM [
50]. The calculated charge density showed that the hole with stronger in-plane localization extended to both layers (modulated by interlayer moiré potential), while the delocalized electron is confined at the MoSe
2 layer (modulated by intralayer moiré potential). This result suggested a more localized interlayer exciton than intralayer exciton because the bandgap modulation of interlayer exciton was larger than that of the intralayer exciton [
50]. Numbers of studies have observed the intralayer excitons in moiré superlattices. Under moiré potential with a relatively large period, the moiré intralayer excitons and trions in MoS
2/MoSe
2 heterostructure exhibited unchanged decay times compared with those in monolayer at 5 K due to the relatively weak localization effect [
121]. In the reconstructed moiré patterns of WSe
2/WSe
2 homobilayer with spatially alternating AB- and BA-stacked domains, a higher-energy intralayer exciton peak emerged in the p-doped regime of PL spectra (
T = 4 K) at regions with larger
λm, in contrast to those with
λm < 10 nm. As shown in Fig.9(b), the two peaks showed different doping dependence and distinct valley coherence represented by the degree of linear polarization (DOLP) [
141]. These results could be explained by the formation of two exciton species (type I and II excitons). The two intralayer exciton species were distributed in different intralayer locations with different bandgaps in the reconstructed superlattice, as depicted in Fig.9(c). Type II excitons could only be observed in perfectly reconstructed moiré patterns with aligned domains size larger than the exciton diffusion length (estimated to be 10−100 nm), because the interlayer tunneling at the K point was symmetry forbidden in AB/BA-stacked domains, which would prevent the excitons from relaxing to the lower-energy state (type I) before recombining. Type I excitons could bind strongly to extra holes under p-doped to form charged excitons due to the higher VBM of exciton I than that of type II excitons. Whereas, the less-layer-localized electrons were expected to interact with both exciton species under n-doping, since the global CBM was located at the Q point [Fig.9(c) and (d)] [
141]. The doped hole in type I location and exciton I could form an entangled state
under the Pauli exclusion principle. The orthogonal
and
made linearly polarized emission impossible because of the broken symmetry caused by correlated resident holes. On the contrary, the resident hole of type II excitons was in the other layer and carried no information about the exciton valley, allowing a non-zero DOLP. DOLP was non-zero for both exciton species in the n-doped regime, because the electron would be located in a valley (Q/Q′) different from that of the exciton (K/K′) [
141]. The results indicate the possibility of realizing exciton arrays with tunable quantum coherence properties based on reconstructed moiré superlattices.
Compared with intralayer excitons, interlayer excitons with spatial separation and permanent dipole are expected to be more deeply affected by moiré potential, which could serve as interacting bosonic lattices suitable for experimental realization of Bose−Hubbard model and study on bosonic supersolid and superfluid states [
155-
157]. Given the relatively long lifetime and diffusion length of interlayer excitons, modulation of the interlayer excitons diffusion dynamics by moiré superlattice is of great interest [
50,
54]. The exciton diffusion in the moiré system can be described as a range of incoherent hops between sites with local minimum of moiré potential. The confined moiré excitons are expected to have shorter diffusion lengths compared with the delocalized excitons. The observed diffusion length of the moiré exciton was one order of magnitude smaller than that in the untwisted WSe
2/MoSe
2 heterostructure [
26,
50]. Using transient absorption (TA) microscopy, Yuan
et al. [
134] investigated interlayer exciton dynamics and transport in WS
2/WSe
2 heterobilayers under twist angles of 0° and 60° with a 4% lattice mismatch, which had a moiré period of about 7.6 nm. The diffusions of interlayer excitons with non-linear temporal dependence of mean square distance
deviated from the normal diffusions of intralayer excitons in the constituent monolayers [Fig.9(e)]. The deviations suggested the interplay between moiré potentials and strong repulsive exciton−exciton interactions of interlayer excitons. Diffusions of interlayer exciton exhibited twist-angle dependence, corroborating the twist-angle tunable moiré potentials. Specifically, the DFT calculation indicated deeper potential and stronger localization for interlayer excitons at
θ ≈ 0° than those at
θ ≈ 60°. A larger energy difference between the K and Q valleys of the conduction band
in the 0°-twisted bilayer suppressed the efficient electron scattering from Q to K, compared with the 60° twisted case. Consequently, WS
2/WSe
2 heterobilayers presented a shorter K−K exciton decay time and higher temperature for thermally excited phonons to assist the Q−K scattering at near-0° than at near-60°, as shown in Fig.9(f) [
134]. Choi and co-workers [
158] also reported that interlayer exciton diffusion in WSe
2/MoSe
2 heterobilayers could be highly controlled by twist angle, using spatially and time-resolved PL measurements at
T = 4 K, as shown in Fig.9(g). The presence of moiré superlattice with a large supercell and deep potential at small twisting could completely localize interlayer exciton and lead to a reduced lifetime and velocity. Consequently, the diffusion length was decreased in comparison with that in commensurate heterostructure. At a larger twist angle leading to a reduced moiré period, excitons may tunnel between supercells and have a longer interlayer exciton lifetime [
46,
158]. Apart from the low-temperature case, the temperature-activated diffusivity is also dependent on the twist angle due to the deep moiré potentials [
136]. Moiré potential effects can be effectively screened at a sufficiently high exciton density or temperature [
134]. The moiré interlayer excitons of MoSe
2/WSe
2 heterobilayers showed a transition temperature, above which the confined IXs turned to be delocalized, manifesting as significantly decreased PL intensity. The temperature dependence of IXs was considered to result from IX−IX exchange interactions exceeding the moiré potential traps at a sufficiently high temperature [
159]. Brem and Malic [
155] have newly reported the theoretical prediction that the strong delocalization of moiré-confined interlayer excitons would occur at intermediate densities (10
11‒10
12 cm
−2) due to the combination of strong repulsive IX−IX interactions and their bosonic nature. The delocalization causes a blueshift and the collapse of the split series of optical moiré resonances, which would merge into a single free exciton peak. The interlayer exciton lifetime has also been demonstrated to be influenced by moiré potential [
160]. Based on a low-energy continuum mode, the twist-angle-dependent IX lifetimes in twisted TMD bilayers should consider two main mechanisms: One is the shift to indirect transitions due to the rotation momentum mismatch under an increasing twist angle. The other mechanism is that the presence of moiré potential introduces an interaction among the IXs with different COM momenta and opens additional radiative recombination channels, resulting in the relaxed momentum conservation and consequently a softened twist angle dependence [
160].
In addition to interlayer twist and temperature, various factors in twisted vdW heterostructures have shown the ability to modulate the localization of moiré excitons in recent researches. Cai and co-workers [
161] focused on the formation mechanism of numerous PL peaks from localized interlayer excitons in moiré superlattices. They constructed an interlayer donor-acceptor pair exciton (DAP IX) dynamic model to explain the origins of these narrow excitonic peaks, where the localized DAP exciton has an electron trapped in the donor site and a hole trapped in the acceptor site. They also fabricated an hBN-encapsulated MoSe
2/WSe
2 bilayer (
θ ≈ 60° ± 1°) on a nanopillar array, and experimentally confirmed the DAP IX model through observations of the emission energy-lifetime correlation and unique nonmonotonic power dependence of the energy-resolved lifetime. The nanopillars here produced strain to provide PL emission enhancement through strain gradient-induced funnelling as well as the increased coupling between the band-edge and DAP IX, indicating the modulation of moiré excitons from out-of-plain strain. Coincidentally, Lin
et al. [
143] have demonstrated the influence of atomic-level in-plane strain on moiré potential. They measured the moiré potential for excitons in a large amount of MoSe
2/MoS
2 heterobilayers with various twist angles. Surprisingly, the moiré potential depth for intralayer excitons was found to be up to ~130 meV, which was larger than the predicted value of most DFT calculations and comparable to the potential depth for interlayer excitons. The deepened moiré potential depth for intralayer excitons was attributed to the remarkable in-plane strain distribution caused by atomic reconstruction, which would induce band gap variation up to ~100 meV in MoSe
2. Besides vdW heterobilayers, multilayer-based systems can provide new degree of freedom for moiré superlattice regulations. Zheng and collaborators reported that interlayer excitons in both twisted WSe
2/WSe
2/WSe
2 and WSe
2/WS
2/WSe
2 trilayers were confined by deeper moiré potentials compared with bilayer systems, manifesting as multiple dense and sharp emission peaks in PL spectra [
162,
163].
For the circumstance suitable for the resonant hybridization model, twist angle modulate the properties and dynamics of interlayer and intralayer hybrid excitons (hXs) by changing moiré reciprocal lattices [
135]. The properties of the hXs such as oscillator strengths, resonance energies, and mixing of intralayer and interlayer components, are sensitively dependent on the moiré period (twist angle) [
135,
146], as shown in Fig.9(i). It is notable that though large detuning between interlayer and intralayer states would suppress hybridization over a range of angles deviating from aligned or anti-aligned cases, pronounced hybridized moiré excitons with strong interlayer tunneling occur under twist angles near 21.8° and 38.2°, These high-angle moiré excitons were observed in WS
2/MoSe
2 bilayers by reflectance contrast spectra at 4 K. The features and valley configurations of moiré hXs at
θ = 21.8° and
θ = 38.2° resembled those at
θ = 0° and 60°, respectively [
135]. Moiré reciprocal lattice at these angles are commensurate with the monolayer lattices in Brillouin zones and therefore the angularly mismatched valleys of constituent monolayers are connected by primitive moiré reciprocal lattices vectors. The commensurate moiré lattices allow nearly resonant tunneling between the intra- and interlayer states, which contributes to the observable hX features, as depicted in Fig.9(h) [
135]. Furthermore, the moiré reciprocal lattices at
θ = 21.8° and 38.2° have the largest reciprocal lattice constant equivalent to
of the monolayer reciprocal lattice constant, corresponding to the smallest moiré lattice with Σ = 7. Σ is the coincidence site lattice (CSL), defined as the smallest unit cell volume shared by the commensurately twisted bilayers. The unit cell maintains the in-plane rotational and translational gliding symmetry [
137]. Zhao and co-workers [
137] systematically investigated the optical properties of TMD homobilayers with various twist angles. For the MoS
2 homobilayers with Σ = 7 (
θ = 21.8°), an intense peak of indirect moiré excitons (1.6 eV) appeared in the PL spectra at 10 K. The indirect moiré excitons showed PL intensity stronger than that of direct intravalley excitons (1.9 eV), which was attributed to a smaller effective mass at the multiple valley points in the conduction bands. Similar commensurate structures at Σ = 13 (
θ = 27.8°) also showed a remarkable moiré exciton peak, but its relative intensity with respect to the intravalley exciton decreased due to a flattened conduction valley. These high-angle commensurate twisted bilayers exhibit optoelectronic properties different from those in the small twist angle cases and are commonly found in as-prepared samples without deliberate alignment, which would motivate further studies into dynamics and valley properties of these high-angle moiré excitons [
137].
Moiré system based on twisted TMD heterostructures can form large-scale superlattice with spatially modulated bandgap confining excitons in localized orbitals [
52], providing a strongly interacting regime [
122]. The moiré superlattice offers a compelling platform for exotic many-body interactions phenomena, such as exciton BEC, especially for the indirect interlayer excitons with interlayer electric dipole and relatively longer lifetimes [
122]. Recently, Remez
et al. [
122] have shown that the strong interactions in twisted TMD bilayers can drive “leaky emission” originate from “dark condensate” where momentum-indirect excitons form a Bose−Einstein condensate. These “leaky emissions” originate from strong interactions and dominate the photoluminescence of TMD moiré excitons at low temperatures with distinctive features.
3.3 Spin-valley configuration of moiré excitons
Excitons in moiré superlattice experience spatially varying periodic potential due to the local bandgaps, determined by the atomic registry. As shown in Fig.7(c) and (e), intralayer excitons and interlayer excitons confined in the local energy minima have different landscapes, accompanied with distinct spin-valley optical selection rules [
49]. In near-0° twisted heterostructures and homostructures, valley optical selection rules for moiré intralayer excitons were maintained, presenting a degree of circular polarization (DOP) very similar to that in the monolayer [Fig.10(a)]. In conclusion, the moiré pattern has no influence on the intralayer optical transitions due to the unchanged rotational symmetry in the plane of the monolayer [
63,
121,
145].
On the contrary, optical selection rules for interlayer excitons exhibit extra complexity. The intralayer excitons excited by circularly polarized light will then relax through the interlayer transfer process to form interlayer excitons. The spin and valley indices of carriers would be conserved in this process [
65]. Unlike intralayer excitons, interlayer excitons with certain spin-valley configurations may generate emissions with circular polarizations opposite to those of the excitation light. Take R-type stacking MoX
2/WX
2 heterobilayer (X denotes chalcogen) as an example. As mentioned in Section 3.1, moiré interlayer exciton wave packets at the high-symmetry locals have distinct
transformations described by Eq. (3.4). Hence the corresponding convertible photons must have the same rotational symmetry [
65]. Therefore, the calculated optical selection rules would vary spatially, as depicted in Fig.7(e): The spin-up exciton wave packet at local A (B) only couples to light with
circular polarization, while light coupling at C is forbidden. At other locals without
symmetry, optical selection rules manifest as elliptical polarization [
49]. Wu
et al. [
63] investigated the valley polarization of interlayer excitons in WSe
2/WSe
2 twisted homobilayer by measuring the polarized-resolved PL spectra at 8 K. Here, interlayer excitons showed two emission peaks with opposite polarizations, as shown in Fig.10(b) and (c). The distinct behaviours of the two IX peaks were attributed to the spin−orbit coupling effect. The conduction band split into spin-up and spin-down energy levels, where the two kinds of interlayer excitons with different transition energies exhibit distinct spin configurations, depicted in Fig.10(c) [
63]. In contrast to the monolayer TMD where the spin-triplet excitons are considered dark, the optical transition dipole of the spin-triplet interlayer exciton is expected to be comparable to that of the spin-singlet exciton in TMD heterobilayer due to the breaking of the out-of-plane symmetry [
66]. Yu
et al. [
66] reported that as the local interlayer atomic registry changes, the optical transition dipole of interlayer exciton wave packet varies between two in-plane transition dipoles of opposite circular polarizations and the out-of-plane one, which enables the spatial modulation of IX light-coupling including the strength of optical transition dipole and the optical selection rule in a large scale moiré pattern. Additionally, the spin-triplet and spin-singlet excitons have distinct valley polarization selection rules at a certain atomic registry, as shown in Fig.10(f). As a result, the selective optical addressing of both the valley and spin configurations of IXs would be achievable. Zhang and collaborators [
164] identified the spin-singlet and spin-triplet exciton states with very high valley polarizations (over 80%) and opposite helicities in high-quality hBN encapsulated WSe
2/MoSe
2 heterobilayers with near 60° (58.7° ± 0.7°) twist angle at 4 K [Fig.10(g)]. The two exciton states with high valley polarization and short lifetimes originated from direct band-gap transition, which were regarded to be localized at the
atomic registry (local point A as the potential minimum). In addition, the long interlayer exciton valley depolarization time suggested strongly suppressed intervalley exchange interactions, on account of both spatial separation of interlayer exciton as well as the discretized spin singlet and triplet states [
164]. Jin
et al. [
67] also observed that the coexisting multiple interlayer exciton states localized at different sites in moiré supercell and exhibited opposite optical selection rules in a near-0° twisted WSe
2/WS
2 heterostructure at 10 K. In this work, PLE and resonant pump-probe spectroscopy were utilized to detect the relative oscillator strengths and spin-valley configurations of the moiré interlayer excitons. The conception of exciton quasi-angular momentum (QAM) was adopted here, which describes the phase change of the exciton wave function after the
transformation [Eq. (3.4)]. The total QAM of an interlayer exciton state directly determines its optical selection rule, which is related to the spin, valley, and moiré pattern. The spatially varying moiré QAM associated with different interlayer lattice registrations explained the opposite selection rules of multiple interlayer excitons in the WSe
2/WS
2 heterostructure [
67].
The presence of moiré confinement for interlayer excitons can realize perfect arrays of quantum emitters, which are highly tunable by several external methods [
49]. Seyler
et al. [
62] studied the valley-polarized nature of the trapped interlayer excitons of MoSe
2/WSe
2 heterobilayers at
T = 1.6 K. The moiré trapped IXs exhibited strong valley polarizations with selection rules depending on the twist angle, as shown in Fig.10(d), corroborating the retention of three-fold rotational symmetry under a sufficient smooth moiré potential. They performed magneto-photoluminescence spectroscopy to determine the Landé
g-factor of the trapped interlayer excitons through the valley Zeeman effect [Fig.10(e)]. The extracted
g-factors were homogeneous across the same sample and took two values, −15.89 ± 0.03 and 6.72 ± 0.02 for near 60° (57°) and 0° (2°), respectively. The
g-factor of an interlayer exciton represents its valley configuration and valley magnetic moment. Thus the two different values indicated the distinct valley-pairing configurations at the two twist angles. At
θ ≈ 60°, the spin-flip optical transitions would lead to emissions with much larger Zeeman shifts, whereas the spin-singlet excitonic Zeeman splitting at
θ ≈ 0° was similar to that of the monolayer. In particular, the trapped interlayer excitons at
θ = 20° (near the commensurate stacking angle 21.8° providing the shortest superlattice periodicity Σ = 7) had a
g-factor close to that at
θ = 57° because of the same valley configurations, while the emission intensities at
θ = 20° were reduced by two orders of magnitude compared with those at
θ = 57° due to the weak optical dipole of the umklapp recombination [
62].
Beyond the moiré system based on heterobilayer, introducing the layer degree of freedom would open up more possibility of tuning the spin-valley configuration through spin-layer locking. Brotons-Gisbert and co-workers [
165] fabricated a trilayer (TL) heterostructure consisting of a monolayer (ML) WSe
2 and bilayer (BL) 2H-MoSe
2 (anti-aligned) and observed the spin-layer locking of moiré-trapped excitons at cryogenic temperature. Two types of interlayer excitons IX
H and IX
R were identified, with electrons localized in the lower and upper MoSe
2 layer, respectively [Fig.10(j)]. In magneto-optical spectroscopy measurements, the two IX species exhibited distinct spin-valley-locked magnetic moments and exhibited different optical selection rules determined by the local atomic registry of the moiré trapping sites, as shown in Fig.10(h) and (i). Such spin-layer locking induces the layer index as a potential knob to engineer few-level quantum systems in vdW heterostructures, providing new opportunities for quantum information applications [
165].
3.4 Moiré system under external tunning
As the novel electronic and optical physics of moiré superlattice based on van der Waals heterostructures has been widely investigated theoretically and experimentally, several researches on control of moiré systems by external tunning methods have been reported in recent years, such as electric fields [
148], carriers doping [
154,
166], magnetic field [
62,
165], temperature [
136,
159], strain [
144,
167,
168] and so on. Thanks to the remarkable sensitivity to external factors, modulations on moiré systems based on 2D vdW heterostructures have shown desirable performance and novel physics.
An applied external electric field can effectively tune the properties of interlayer excitons due to their static electric dipole perpendicular to the plane (
, where
is the charge quantity and
is the vector denoting charge separation from negative charge to positive charge) [
50,
54]. The coupling between the perpendicular electric field and the electric dipole would shift energy levels due to the Stark effect, which can realize modulations of the resonant energy and polarity of interlayer exciton. Consequently, the in-plane distributions of moiré interlayer excitons are controllable by external electric fields. Guo
et al. [
50] used first-principle calculations to predict the charge density distributions tuned by the electric field in the twisted MoS
2/WS
2 heterostructure, as shown in Fig.11(a). Furthermore, the shifts of CBM and VBM from different constituent monolayers would control the interlayer resonant tunneling of the electron or hole and therefore the formation of hybridized excitons. In addition, an electric field can change both the electric dipole orientations and dipole moments of hXs [
50]. Tang
et al. [
148] employed a dual-gate device [Fig.11(b)] to independently apply a vertical electric field and doping to the WSe
2/WS
2 (type II band alignment) and MoSe
2/WS
2 (type I band alignment) heterobilayers at 5 K. They identified layer-hybridized moiré excitons through an emergent feature iX in the doping-dependent reflectance contrast spectrum, which was only observed in anti-aligned (
θ = 60°) WSe
2/WS
2 heterostructures. The configuration of anti-aligned WSe
2/WS
2 bilayers [Fig.11(b)] enabled extra Bragg reflection from the moiré potential to help conserve quasi-momentum and promote strong coupling. Therefore, the resonant interlayer tunneling was regarded as spin-conserving. Under an out-of-plane electric field, strong exciton coupling was observed in WSe
2/WS
2 heterobilayers with coupling constants comparable to the intralayer moiré potential (~25 meV), manifesting as energy-level anticrossing and electric-field-dependent oscillator redistribution between the interlayer and intralayer resonances, as shown in Fig.11(c). The IX peak exhibited linear electric-field response and weak oscillator strength when it was decoupled from other excitons under a sufficiently large positive perpendicular electric field, indicating the interlayer exciton nature of IX. The angle-aligned MoSe
2/WS
2 heterobilayers had much weaker coupling with the extracted coupling constants smaller than the exciton linewidths (about 5−10 meV), but still showed discernible multiple layer-hybridized moiré excitons with two distinct electric-field dispersions. In addition to layer hybridization, electric fields are also expected to control moiré exciton transport. Based on the Förster theory, the relative orientation of the dipole moments plays an important role in the modulation of moiré exciton diffusion [
50]. Lately, several works have reported the electric-field-tunable spatial distribution and transport of excitonic species in 2D vdW heterostructures. Tagerlli
et al. [
170] have achieved electrostatic control over layer-hybridized exciton states represented by effective out-of-plane dipole lengths and revealed the dipole-dependent micrometre-scale transport of hybrid exciton. They modulated the layer hybridization and many-body interactions of excitonic species in 2H-stacked WSe
2 homobilayers by applying a vertical electric field at 4 K. The moiré-less 2H-WSe
2 structures eliminated the potential traps and type-II band alignment, which allowed tunable exciton−exciton interactions and long-range propagating. In the presence of moiré superlattice, the electric-field-dependent hybridization would allow continuous tuning of interlayer exciton hopping between different moiré sites, and generate a mixture of these interlayer excitons as a superposition moiré states, which has been demonstrated in a trilayer WSe
2/monolayer WS
2 heterostructure reported by Lian and co-workers [
171]. Besides, in the R-stacked TMD vdW heterostructures accompanied with breaking out-of-plane mirror symmetry, the electric field can act as a switch to flip the out-of-plane polarization through in-plane sliding motion, which would bring about interfacial ferroelectricity [
81,
172,
173]. Utilizing piezoelectric force microscopy (PFM), Wang
et al. [
81] have achieved the visualization of moiré ferroelectric domains and observed the electric-field-induced domain wall motion in several nearly parallel-stacked R-type bilayers.
Doping introduces charged moiré excitons that obey fermionic statistics and exhibit distinct properties such as the long-range Coulomb interaction in contrast to bosonic neutral excitons, providing an accessible platform to study fermionic many-body effects [
166]. Recently, Dandu
et al. [
154] have realized the control of carrier concentration by gate-induced doping in H-type stacked twisted MoS
2 bilayer, and investigated the coexisting localized and delocalized intralayer moiré excitonic species at 4.5 K. As
increased, the energies of intralayer exciton peaks stayed relatively unchanged, whereas their intensities (except the localized exciton X
I) quenched rapidly accompanied with emergent peaks with lower energy, as shown in Fig.11(d). The results suggested the unchanged spatial modulation of intralayer moiré potential and the formation of trions from the delocalized higher energy excitons. Through a combination of gating (modulation of carrier density and potential depth) and varying excitation power (exciton density), multiple intralayer moiré trion species would be formed with tunable localization phases, depicted in Fig.11(e). They also demonstrated the gate tunable valley coherence and exciton−phonon coupling strength for the moiré excitons. The charging effects would influence the spin-valley physics of interlayer moiré exciton as well. Wang and collaborators [
166] observed moiré trions as the moiré traps were filled with either electrons or holes in H-stacked WSe
2/MoSe
2 heterobilayer at 1.6 K, manifesting as new sets of narrow peaks in PL spectra, as shown in Fig.11(f). These emergent peaks were about 7 meV below the energies of the neutral interlayer moiré excitons. The similar g-factors of these excitons and trions (about ‒16) implied their same spin-triplet configuration and the smooth moiré potential traps. When positively charged, the valley polarization of moiré trion was inverted to cross-circular polarization under the σ
+ excitation at WSe
2 monolayer exciton resonance, in contrast to the co-circularly PL in the neutral and negative cases. Such inversion was explained by the competition between spin-flip valley-conserved and spin-conserved valley-flip relaxation channels of electrons during the formation of interlayer trions. In both negative and neutrally charged cases (
and
) where the optical selection rules were determined by the spin-valley index of the photo-excited hole, the excited K-valley-polarized hole could form a spin-triplet recombining pair with −K′ valley electron. The resulting high co-circular PL polarization suggested that the valley polarization of the excited hole was well-protected by spin-valley locking and that these multiple election−hole pairs were all trapped at
locals. Under hole-doping (
), the emission polarization was determined by the electron valley configuration since the doped holes would distribute at both opposite valleys. The relaxation channels of photo-excited electrons involved either a spin flip or a valley flip, where the latter was found to dominate the process, resulting in the smaller cross-circular polarization of
[
166]. Moreover, the properties of layer-hybridized excitons including energies, coupling constants and effective dipole moment can also be significantly modulated by doping through electron-exciton interactions under considerable localization of moiré potential [
148]. The rich and complex phenomena of moiré excitonic species under electrical modulation indicate the promising potential of engineering both bosonic and fermionic many-body physics. Furthermore, their dynamic adjustable features suggest the possibilities of photonic and valleytronic devices based on moiré systems.
The excellent inherent mechanical flexibility and remarkable strength of two-dimensional materials make strain engineering an ideal tuning method for manipulating their electronic and optical properties. For example, Both graphene and TMD possess maximum elastic strains larger than 10% [
29]. Strain can realize continuous band structure modulation [
30,
32,
174,
175], and inhomogeneous strain can creat local variation of bandgap [
31,
176-
178]. The stacking structures of twisted 2D vdW heterostructures enable sensitive and dynamic strain engineering of moiré potential landscape [
167] and topological properties [
33]. Moreover, the moiré system only requires small strains to achieve considerable modulation. By applying proper differential (inequivalent strains on the two constituent layers in the order of the lattice mismatch) and uniaxial tensile strain to twisted WSe
2/MoSe
2 heterobilayers, Bai
et al. [
167] transformed the zero-dimensional (0D) moiré pattern into a one-dimensional (1D) structure with parallel stripes and presented real-space imaging of these two patterns using PFM, as shown in Fig.11(g). In the 1D moiré pattern, the photoluminescence emission at 4 K from interlayer exciton showed a single broad peak and significantly higher intensity in contrast to the 0D case. The enhanced broad peak in PL spectra was attributed to the higher density of states with densely spaced excitonic levels and radiative recombination from multiple levels. The emissions of 1D structures were linear polarized and their directions were aligned with the direction of one of the primary structures obtained from 2D-FFT (Fast Fourier Transforms) image analysis [Fig.11(i)]. The linear polarization here was qualitatively attributed to the hybridization of the two valley configurations through electron-hole exchange interaction, because the electron-hole wave functions were highly anisotropic in the 1D moiré potential from the breaking
rotational symmetry. The strain-field-induced 1D moiré potential, analogous to the quantum wires, is attractive for studying strongly correlated and anisotropic charge transport [
167]. Apart from the differential strain, in-situ deformation is another approach to realize dynamic control of the moiré potential and correlated states by directly altering the interlayer separation of the vdW heterostructures, which can help gain a clear understanding of the exciton physical behaviour in moiré system [
144,
168]. Li and co-workers [
144] have performed the
in-situ dynamic modulation of moiré potential for twisted WS
2/WSe
2 heterostructures by employing diamond anvil cell (DAC) hydrostatic pressure engineering at room temperature [Fig.11(h)]. When the pressure was lower than 1.1 GPa, the exciton response of heterostructures (HS) had distinct pressure dependence under different twist angles. For HS-7.5°, PL peaks exhibited blue shift trends due to the changed band structure, in which the energy of direct transition from the K valley would increase under larger pressure. Whereas in HS-1.15°, the predominant moiré effect led to different exciton responses. The PL peaks of HS-1.15° showed continuous redshift with approximately linear dependence for intralayer excitons, which could be explained by the increased moiré potential under smaller interlayer spacing. The exciton dispersion band of the mini-Brillouin zone shifted to lower energy levels under the deepened moiré potential. At pressure higher than 1.1 GPa, effective enhancement of interfacial coupling in all twisted heterostructures dominated the exciton response in PL spectra [
144]. However, conventional experiment methods to realize in-situ deformation like DAC are challenging to apply to moiré heterostructures at low temperatures. Zhao
et al. [
168] utilized an optical fiber tip in a cryogenic scanning near-field optical microscope (SNOM) to locally deform the sample and measure its background-free near-field optical response simultaneously, as illustrated in Fig.11(h). The moiré potential of the near-zero-twisted WSe
2/WS
2 heterostructure increased under compressive deformation at 25 K, resulting in redshift of both intralayer and interlayer excitons. Under the maximum stress of 1.4 GPa, the sample obtained a modulation of 7 meV on the moiré potential coefficient
in Eq. (3.3). The maximum stress was limited by the damage threshold of the optical fiber tip. Therefore, this method requires a tip made from stronger materials to achieve a higher limit [
168]. Cai
et al. [
161] utilized a nanopillar array placed under the hBN-encapsulated MoSe
2/WSe
2 heterostructures to introduce strain-enhanced coupling and strain gradient-induced funneling, which helped them to corroborate and develop the DAP IX dynamic model for understanding the primary origin of moiré localized interlayer excitons.
Moiré system combined with other elements would facilitate more fertile physics and deserve further investigation. Anđelković
et al. [
169] discussed the supermoiré (SM) resulting from two slightly misaligned interfering moiré patterns in the trilayer system. They concentrated on graphene encapsulated by hBN and quantitatively described the angle-dependent SM effects that led to electronic spectrum reconstruction, velocity suppression and spatial electron redistribution. Fitzgerald
et al. [
34] studied twist-angle-tunable moiré exciton polaritons resulting from hybridization between moiré intralayer excitons and photons in the strong coupling regime. They focused on twisted AA-stacked MoSe
2/WSe
2 heterostructures integrated within a Fabry‒Pérot cavity. They theoretically predicted the delocalized nature of the new hybrid light-exciton states and the optical responses of moiré exciton polaritons including exciton-light coupling, polariton energy and the number of polariton branches, which could be modulated by twist angle and cavity [
34]. Moiré patterns have also prompted new progress in diverse fields. For example, moiré phonon physics derived from the atomic stress redistribution and folding effects of phonon dispersion in twisted 2D materials can drastically regulate the lattice vibration modes and thermal transport properties [
91]. Besides, a variety of researches involved with moiré structures have been reported like gap solitons in parity−time symmetric moiré optical lattice [
179], and moiré bilayer hyperbolic metasurface with highly controllable electromagnetic responses [
180]. These research efforts are expected to pave the way to realize flat bands and localized states at arbitrary frequencies for multifunctional devices.
4 Correlated electronic states in moiré superlattice
Moiré superlattice in vdW heterostructures provides a general and promising solid-state platform for exploring and engineering correlated electronic states, supported by the strong periodic moiré potential with the flat electronic band. The core of the flat band lies in the comparability between the electron−electron interaction
and the electronic bandwidth
under sufficient narrow electronic bands, enabling the effects of Coulomb interactions to become pronounced in the system [
181]. We can consider a simplified model of a two-dimensional electron gas in a periodic potential with a period
, accompanied with the formation of a set of electronic minibands. For the lowest electronic miniband, it can be estimated that
and
, where
and
are the effective dielectric constant and electron effective mass, respectively. Strong correlation would become achievable when the ratio of electron-electron interaction to electronic bandwidth satisfies
, which needs sufficient large
[
181]. The theoretical analysis of such systems must consider the correlated electronic effects instead of one electron approximation, which can trigger abundant quantum many-body states such as quantum anomalous Hall insulators and high-temperature superconductivity [
118,
182]. A comprehensive understanding of such exotic physics based on strong electronic correlations requires both in-depth theory and experiment researches that complement each other. In a moiré superlattice, the electronic correlated effects would be enhanced under a large
with reduced bandwidth
, which can be achieved by a small twist angle θ. Since the correlated insulator [
59] and unconventional superconductivity [
60] were observed in magic-angle graphene superlattice by Cao and co-workers, a wealth of moiré-system-based correlated electronic states in twisted vdW heterostructures has been reported [
61,
129-
131,
181,
183-
188].
A real moiré system based on van der Waals heterostructures possesses more complexity beyond the simplified model shown above. Flat bands in the graphene-based heterostructures appear only at a discrete set of magic angles [
123,
127], limiting the tunability of the strongly correlated systems. A sufficiently strong moiré effects require effective interlayer coupling between the constituent materials. Besides, a small enough twist angle provides a large
but simultaneously brings about atomic reconstruction [
132,
140], which would decrease the homogeneity of the moiré superlattice. Given these factors, TMD moiré heterostructures are supposed to offer an ideal system for the exploration of quantum many-body phenomena thanks to their strong light−matter interactions, considerable moiré potential and the reduced degeneracy arising from the coupled spin-valley index [
128,
131,
181]. The electronic band structures of TMDs with relatively large effective mass (
, with
as the free electron mass) facilitate flat band features and many-body effects [
181]. Wu
et al. [
118] have presented that TMD moiré flat bands could exist over a range of twist angles (smaller than around 3.5°), enabling continuous tuning of correlation effects. They used the continuum model to describe the low-energy electrons in the moiré superlattice, in which the Hamiltonians were spatially periodic under effective moiré potential [
64,
118]. The moiré band Hamiltonian for twisted TMDs heterobilayers had only two low-energy valence bands when the chemical potential was within the topmost valence bands, and could be mapped to triangular-lattice single-band Hubbard models with the Wannier wave function localized in one of the moiré potential extremum positions. The estimated hopping parameters
(related to the bandwidth
) and on-site Coulomb repulsion energy
(sensitive to the effective dielectric constant
) both varied with the moiré period, resulting in a correlation strength that could be modulated over a wide range by twist angle. The relatively strong on-site repulsion would suppress double occupation of one moiré site when satisfying
, and therefore the Mott insulator ground state was predicted at hole half-filling (filling factor
considering the Hubbard model spin degeneracy, signifying one electron per moiré unit cell) [
118]. Hubbard model is a simple theoretical model of interacting quantum particles in a lattice, which is considered to be paradigmatic in the physics of complex quantum many-body effects [
129,
189,
190]. The TMDs moiré superlattices are regarded as an analogue quantum simulator of the Hubbard model with widely tunable parameters like bandwidth, interaction strength and band filling, which could play an important role in revealing the physics of various strongly interacting quantum particles [
118,
128-
131,
183,
191]. Wang
et al. [
128] have reported the observations of the correlated insulator and zero-resistance regions in bilayer WSe
2 under twist angles ranging from 4 to 5.1°, indicating the continuously tunable collective phases in the TMD-based moiré system. Utilizing transport measurements, the effective electronic interactions
here were demonstrated to be tunable by external electric field and doping, driving a continuous metal-insulator transition (MIT) at near half-filling of the first moiré subband (one hole/electron per superlattice site) [
183,
191]. The realization of the highly controllable triangular-lattice Hubbard model offers new avenues in exploring the correlation insulating phase transitions [
191] that would potentially relate to exotic states such as quantum spin liquids [
192,
193] and unconventional superconductivity [
194]. Furthermore, TMDs moiré superlattice are expected to provide a promising platform for engineering novel correlated states, such as emerging correlated valley and spin physics thanks to the layer-valley-spin locking in TMDs [
165,
185,
195].
The strong light−matter interactions of TMD heterostructures enable optical spectroscopy measurements to detect the strongly correlated phases, complementary to the traditional transport measurements [
93,
181]. Regan
et al. [
181] have realized the optical detection of strongly correlated phases in near-zero twisted WSe
2/WS
2 heterostructures. They developed a novel measuring technique called optically detected resistance and capacitance (ODRC) [Fig.12(a)], which could avoid the large contact resistances that impeded direct electrical transport measurements. As shown in Fig.12(b), they observed the Mott insulating state at hole doping factor
(
;
was the moiré supercell density) and additional insulating states at
(
) and
(
) in the heterostructures with moiré period λ ≈ 8 nm at 3 K. The additional insulating states were attributed to generalized Wigner crystallization of holes in the moiré superlattice, where holes tried to avoid simultaneous occupation of adjacent sites [Fig.12(c)]. The emergence of generalized Wigner crystallization arises from the relatively large inter-site interaction compared with the moiré miniband bandwidth, indicating the necessary of an extended Hubbard model that considers inter-site (long-range) interactions in addition to on-site (short-range) interactions in TMDs moiré superlattice [
181]. Tang et al. [
129] have also reported the observation of Mott insulating state in aligned WSe
2/WS
2 heterobilayers using the close-cycle cryostat magneto-optical measurements. The Mott insulating state occurred at half-filling of the first hole moiré band in the limit of
with antiferromagnetic Curie−Weiss behaviour [Fig.12(d) and (e)], suggesting the antiferromagnetic interaction between the moiré band holes. Xu and co-workers [
130] have observed nearly two dozen correlated insulating states at several fractional fillings of WSe
2/WS
2 moiré superlattice. They utilized a new optical sensing method based on dielectric-sensitive excited-state excitons in WS
2 monolayer as a probe of correlated states. As depicted in Fig.12(f), the cascade insulating states with energy ordering nearly symmetric about filling factor
(half a particle per moiré supercell) were attributed to a series of charge-ordered states ranging from generalized Wigner crystals to charge density waves, which further corroborated the remarkable long-range interactions in TMD-based moiré systems.
Shimasaki
et al. [
131] have reported the incompressible electronic states probed by the exciton-polariton resonances in a long-period moiré superlattice. They fabricated a MoSe
2/hBN/MoSe
2 homobilayer with
λ = 24 nm and a large
parameter (
, the ratio of interaction energy to kinetic energy). The inserted hBN layer weakened the moiré potential and increased the on-site interaction strength, making the homobilayer moiré system a unique regime beyond the standard Fermi-Hubbard model. They controlled the hybridization as well as the doping of the top and bottom layers independently. The hybridization between intralayer and interlayer excitons with opposite dipole moments was regulated through electric-field-dependent coherent carrier tunneling. The charging configuration of the two layers could be simultaneously detected by the attractive polaron (AP) and repulsive polaron (RP) branches in the differential reflectance spectrum due to the dynamical screening of excitons by itinerant carriers. Negative compressibility was observed as the top gate voltage increased while the back gate fixed, manifesting as the depletion of electrons from the bottom layer due to the intralayer exchange interaction [Fig.12(g)], indicating the prominent electron−electron interactions at high electron densities [
131]. In the low-electron-density regime at half-filling of the lowest moiré subband (
), abrupt collective interlayer transfer of electrons was observed with resilience towards electric-field-induced charge redistribution, as shown in Fig.12(i). This result suggested the existence of Mott-like incompressible electronic states, which could be accommodated in either the top or the bottom MoSe
2 layer. Similar interaction-induced incompressible electronic states with abrupt phase transitions also appeared for other doping levels equal to an integer multiple of the moiré supercell density. This kind of moiré system was expected to act as a new platform for studying the interaction between dipolar excitons and electrons confined to flat bands, like Bose−Fermi mixtures [
131].
As mentioned above, TMDs moiré superlattice have emerged as an ideal system for studying correlated electronic physics with high tunability, where the onsite Coulomb interaction
and the nearest-neighbour hopping parameter
are controllable by varying twist angle [
128] and electric field [
183,
191]. Apart from these methods, Wang
et al. [
187] and Chen
et al. [
185] have introduced other factors that could modulate the electronic correlation in TMDs moiré superlattice. The former studied the impact of optical excitation on the correlated electronic states in WS
2/WSe
2 heterobilayers. They showed that optical excitation enhanced magnetic interaction between moiré-trapped carriers and consequently led to ferromagnetic order [
187]. The hysteresis loop at low filling factor
measured by the reflective magnetic circular dichroism (RMCD) at
T = 1.6 K, signified the emergence of ferromagnetic order due to long-range spin-spin interaction in generalized Wigner crystal state. The power dependence of RMCD responses [Fig.12(h)] indicated the remarkable optical tuning of spin−spin interaction strength between moiré-trapped holes mediated by itinerant photoexcited excitons, which could realize the modulation of magnetic phases [
187]. Chen and co-workers [
185] introduced the layer degree of freedom in the angle-aligned WSe
2/WS
2 heterostructure with WSe
2 varying from monolayer to trilayer, in order to systematically tune both electronic flat bands and moiré exciton bands. The systematic changes including emergent exciton resonances and increased resonance energy separations in the differential reflectance spectra suggested that moiré coupling was strongly confined at the WSe
2/WS
2 interface with limited out-of-plane extension. In the multilayer WSe
2/1L WS
2, the added WSe
2 layer(s) could affect the hybridization between moiré excitons and the interlayer-like hybrid exciton. Microwave impedance microscopy (MIM) measurements demonstrated the tunability of the electron correlation by layer degree of freedom, as illustrated in Fig.12(n). As the layer number of WSe
2 increased, the correlated insulating states became fewer and less pronounced in the MIM spectra accompanied with reduced Mott transition temperature. The reduction of transition temperature indicated the weaker inter-site and on-site electron interaction strengths, which could be explained by reduced interfacial Coulomb interaction under enhanced dielectric screening and the increased kinetic energy originated from the larger bandwidth of the flat band in multilayer WSe
2/1L WS
2. It is worth noting that even the weakened electron correlation (with the Mott transition temperature about 60 K for 3L WSe
2/1L WS
2) reported here was still stronger than that in the graphene moiré system (with transition temperature of about 4 K).
Taking advantage of the rich excitonic physics in TMD vdW heterostructures, several new correlated insulating phases of excitons have been reported in recent years. Zhang
et al. [
188] integrated the moiré superlattice into a “double-layer system” composed of a near-aligned WS
2/WSe
2 heterobilayer (moiré superlattice) and a WSe
2 monolayer separating by an ultrathin hBN layer, aiming at realizing flat electronic bands and strong interlayer electron−hole interactions simultaneously. They controlled the total hole concentration
and the average electric field independently. The electric field determined the hole distribution between the WSe
2 monolayer
and the WSe
2 from moiré bilayer
, following a relation of
. As finite hole density was distributed in both WSe
2 layers, a higher total
should be required to achieve the moiré bilayer Mott insulator state (
, with
denoting one hole per moiré lattice site), while actually a stable correlated insulator state appeared at
[Fig.12(m)]. The unexpectedly low
originated from the strong interlayer Coulomb interaction preventing the holes in the WSe
2 monolayer from occupying the same sites with holes from the WS
2/WSe
2 moiré bilayer. This correlated insulating phase was interpreted as an interlayer exciton insulator through particle-hole transformation, where holes in the WSe
2 monolayer and the doped electrons in the WS
2/WSe
2 moiré bilayer could form tightly bound interlayer excitons, as depicted in Fig.12(l). This interlayer exciton insulator would evolve to other correlated phases such as metal-like optical response under sufficient high
, and could potentially form exciton condensate and counterflow superfluidity [
188]. Recently, Xiong
et al. [
196] have also reported an exciton incompressible state as a bosonic correlated insulator in WSe
2/WS
2 moiré superlattices persisting above 30 K. Using a unique pump-probe spectroscopy method to independently tune the charge and exciton density through electrostatic gating and pump light respectively, they achieved continuous transitions between the bosonic correlated insulating phase and the electron correlated insulator where excitons avoid electron-occupied sites under on-site electron‒exciton repulsion. These excitonic correlated states observed in semiconducting moiré superlattices open up new ways for manipulating emergent bosonic phases and many-body interactions involving both fermions and bosons [
196].
The moiré excitons behaviours under optical measurements with limited spatial resolutions are commonly described by a continuum model [
118,
197], which considers the center-of-mass motion of the pristine exciton modulated by the effective moiré potential, assuming that the internal spatial electron-hole correlations from the pristine single layers are retained [
186]. The prevailing effective continuum model is limited to large moiré periods and cannot capture the microscopic nature of exciton states in real moiré superlattice, where electrons and holes could experience different moiré effects and exhibit distinct spatially modulated moiré flat bands [
133,
139,
198]. Combining large-scale first-principles GW-BSE predictions with micro-reflection spectroscopy, Naik and co-workers [
186] have identified various exciton states with distinct electron−hole correlations in rotationally aligned WSe
2/WS
2 reconstructed moiré superlattice at
T = 1.6 K, including a Wannier-type exciton with lowest energy and a new-found intralayer charge-transfer exciton with higher resonance energy. Different from the Wannier-type moiré exciton with tightly correlated electron and hole, the intralayer charge transfer exciton had electron and hole densities distributed with about 5 nm separation in the WSe
2 layer. The distinct spatial electron−hole distributions resulted in different carrier densities and magnetic-field dependence. As shown in Fig.12(k), the moiré peak I (Wannier exciton) was almost unaffected by doping except for the formation of Mott state at one doped hole per moiré unit cell, while the intralayer charge exciton resonance (the moiré peak III) was strongly modulated by carrier doping due to the Coulomb interactions and large dipole of charge-transfer exciton [Fig.12(j)]. The attractive electron−hole interaction between doped electrons at local A led to the emergent peak III
A close to peak I. As electron doping increased, the peak III
A showed blueshifts and consequently reduced the resonant hybridization between peak I and peak III
A, resulting in the weakened oscillator strength. The brief enhanced oscillator strength of the peak III
A at around 1/3 and 2/3 filling was ascribed to the reduced screening of correlated states. The peak II exhibited a mixed character. The non-trivial intralayer charge-transfer exciton state with specific spatial characters can be more easily modified and dissociated, suggesting a new approach to engineering many-body physics and desirable optoelectronic nanodevices based on moiré system [
186]. Such charge-transfer excitons with long lifetimes have lately been reported in monolayer WSe
2/twisted bilayer graphene (TBG) by Hu
et al. [
199], due to strongly asymmetric interlayer Coulomb interactions under a relatively large moiré wavelength compared to the exciton size. At a quite small twist angle (
θ = 0.6°), they observed Rydberg moiré excitons where the accumulated charges in TBG strongly attracted the loosely bound 2s excitons in WSe
2, resembling the Rydberg atom trapped by arrays of optical tweezers. The moiré system brings about highly tunable properties and mutual interactions of the spatially confined Rydberg excitons, demonstrating great potential for the applications of quantum information processing and quantum simulations based on the excitonic Rydberg states in solid-state systems [
199].
5 Conclusions and outlook
This review systematically summarizes the recent progress in twistronic and moiré exciton physics in van der Waals heterostructures that consist of 2D semiconducting materials. The interlayer twist angle can act as an effective knob to modulate the electronic structures and optical properties through interlayer interactions and momentum-space mismatches. The most concerning aspect is exciton physics, where the exciton properties such as binding energy, static electric dipole, emission intensities of different excitonic species [
87,
88] and degree of spin-valley polarization [
77], are all controlled by the interlayer twist angle. Moreover, the interlayer twist in momentum space also affect the excitonic dynamics including formation, transport and recombination processes with complex mechanisms that involve localized and hybridized valleys [
52,
79,
85]. The twist-angle-tunable physics is significantly enriched when the moiré effects are considered, which still lacks a comprehensive model to understand moiré excitons in various systems. Early theoretical analyses have predicted the Moiré Brillouin zone of superlattice structures as well as the moiré patterns-induced spatial modulations of exciton energy and wave function symmetry, which will lead to some novel phenomena such as splitting resonance energies, spatially modulated optical selection rules and localization of moiré excitons, etc. [
49,
52]. The harmonic moiré potential model [
65] describes the moiré excitons as particles moving under approximated potentials with spatial periodicity. This model is suitable for heterostructures with relatively large band edge offsets, but is not sufficient to capture the excitonic physics in systems with near-resonant band edges, where interlayer and intralayer hybridization through carrier tunneling should be considered [
146]. Apart from the moiré superlattice under small twist angles, the enhanced moiré exciton features also appear in high-angle twisted heterostructures due to the umklapp processes in commensurate moiré reciprocal lattice [
135,
137]. Moiré systems based on vdW heterostructures exhibit expanded properties and dynamic tunability in conjunction with external factors, such as electric fields, carrier doping and strain, which would promote an in-depth understanding of moiré physics and pave the way to promising applications. Recently, the enhanced electron−electron interactions and flat band features in moiré vdW heterostructures have facilitated research on strongly correlated electronic states at a few certain filling factors, such as Mott insulating states, generalized Wigner crystal and other correlated phases. These studies are expected to experimentally realize novel quantum phenomena like fractional quantum anomalous Hall effects and spin-liquid states [
118].
In this review, we mainly focus on the excitonic behaviours in twisted vdW heterostructures revealed by optical spectroscopy, which only covers a part of new physics in twistronics and moiré superlattice. There are still a large variety of controversies and challenges in current research that need to be tackled, even within this certain topic. A comprehensive model to understand the excitonic physics modulated by moiré pattern has yet to be established. Excitons presented distinct behaviours in different twisted heterostructure systems. On one hand, such variations are possibly caused by the randomness of twisted heterostructures, which require preparation and detection technology with high precision. On the other hand, the impacts of the moiré pattern on electronic band structures and optical properties in vdW heterostructures are still not clear, such as strain distributions and atomic reconstructions under lattice mismatches [
92,
143]. Besides, the twist-angle-dependent electron-phonon coupling strength in combination with moiré phonon physics in twisted 2D vdW materials might deserve further study [
91]. The carrier dynamics and transport influenced by interlayer twist show much more complicated mechanisms other than interlayer momentum mismatch. For example, the ultrafast interlayer charge transfer processes are relatively robust due to the mixed valley serving as an intermediate in some TMD heterostructures, whereas they show obvious twist angle dependence in other cases [
85,
86]. Interlayer stretching and shift caused by atomic reconstructions are also expected to affect the carrier dynamics through generating extra channels [
79]. In addition, moiré excitons exhibit splitting resonance energies and localization introduced by moiré periodic potential. The cascade transitions between moiré interlayer excitons with different energies have been reported recently, suggesting the spatial overlap between the multiple IX minibands in TMDs-based moiré superlattice [
200]. The links between excitons distributed in different moiré sites and different energies still need in-depth exploration. Moreover, the study on strong moiré excitonic features observed in high-angle commensurately twisted bilayers at several twist angles is still preliminary but interesting, since these cases are widely found in randomly stacked TMD bilayers without the need for deliberate alignment [
137]. Apart from the moiré exciton properties and dynamics, correlated electronic states in TMDs moiré superlattice have spurred great interest and generated many open questions in recent years. The TMDs-based moiré systems working as designed quantum simulation platforms provide many opportunities for new quantum states. It is inspiring that the strong light−matter interactions of TMDs allow optical probing methods such as ODRC [
181] and differential reflectance spectrum [
131,
188], which utilize excitons as sensors. However, the microscopic mechanisms of the interplay between correlated electrons and excitons remain obscure. The nature of various correlated states identified at certain filling factors needs further investigation.
The optical measurement techniques used for probing moiré excitons so far still have limited spatial resolutions in the far-field limit, which are unable to detect the excitonic physics within an individual moiré supercell [
93]. It has been theoretically predicted that electrons, holes, and excitons with different energies have distinct spatial distributions and modulated band structures in moiré superlattice, which are challenging for conventional optical probing methods, especially considering the spatial inhomogeneity and defects. Therefore, non-invasive experimental techniques with sufficiently high spatial- and momentum-resolutions are indispensable for mapping the exact charge distributions and revealing their microscopic nature and mutual correlations in moiré systems. Techniques like scanning near-field optical microscopy [
201], PFM [
58], STM spectroscopy [
202], and nano-ARPES [
153,
203] are promising for future explorations that might eliminate ambiguity in early researches. Besides, the heterogeneity and intricate many-body effects of moiré excitons can possibly be captured by the optical 2D coherent spectroscopy (2DCS), leaving a wealth of opportunities for exploring the interaction between moiré superlattice and excitonic quantum dynamics [
204].
Moiré system based on vdW heterostructure inherits the high tunability and ability to be customized in multiple ways. Moiré vdW heterostructures combined with various external factors such as cavity and strain are interesting and intriguing for further investigations. For instance, strain engineering of twisted vdW heterostructures can cause efficient tunning of optoelectronic properties due to the nano-scale structures of the moiré pattern [
167]. Moiré superlattice can also couple to multiple 2D vdW materials with different stacking sequences and orders to construct vdW heterostructure possessing new electronic and photonic physics, including controllable ferroelectricity [
20,
81,
173,
205] and magnetism [
206]. Beyond the most studied TMDs and graphene reviewed here, the material choice to fabricate a vdW heterostructure is wider regardless of lattice constants matching. The research on twisted vdW heterostructures will prompt engineering artificial materials with moiré superlattice, providing a fascinating platform for exploring emerging quantum phenomena and valuable applications.