WS₂/WSe₂ Hetero-Bilayers & Moiré Physics

• WS₂/WSe₂ hetero-bilayers are atomically thin TMD systems where vertically stacked monolayers generate moiré superlattices that tune electronic and excitonic properties.
• They exhibit type-II band alignment with robust interlayer exciton formation and twist-angle–dependent modulations in binding energies and diffusion behaviors.
• Advanced simulation techniques, including DFT, GW, and hybrid-DFT, reveal Bose–Fermi–Hubbard physics, opening pathways for quantum emitter and device applications.

WS$_2$/WSe$_2$ hetero-bilayers are vertically stacked, atomically thin transition metal dichalcogenide (TMD) systems, wherein a monolayer of tungsten disulfide (WS$_2$) and a monolayer of tungsten diselenide (WSe$_2$) form a van der Waals superlattice with distinctive electronic, optical, and many-body correlated properties. The 4% lattice constant mismatch and variable twist angle between layers generate moiré patterns that fundamentally alter the landscape of the electronic band structure, exciton physics, and collective electronic phases. These systems exhibit type-II band alignment, strong interlayer exciton formation, twist-angle-controllable moiré minima, and host lattice-encoded Bose–Fermi–Hubbard many-body Hamiltonians.

1. Structural and Electronic Ground State Properties

In WS$_2$/WSe$_2$ hetero-bilayers, both edge termination and stacking registry govern thermodynamic stability and electronic band structure. Chalcogen-terminated (c-term) zigzag edges, in both AA and AB stacking, yield robust bilayer formation, whereas metal-terminated (m-term) edges result in significant bond distortions and unstable, frequently metallic, configurations. Density Functional Theory (DFT) indicates that for c-term flakes, AA stacking produces the most negative binding energy (–2.645 eV for triangular W$_{20}$S$_{30}$Se$_{30}$), followed by AB (–1.906 eV), while m-term stackings are thermodynamically unfavorable—AA(m-term) is strictly endothermic, precluding its equilibrium formation (Ghatak et al., 2019).

Interlayer separations vary with edge type: c-term AA and AB are tightly bound (6.25–6.46 Å), whereas m-term are expanded (6.67–6.84 Å). The partial density of states for the c-term, AA-stacked bilayer shows a semiconducting gap with chalcogen (S 2p, Se 4p) valence and W 5d conduction band edge character. GW and hybrid-DFT corrections increase the magnitude of this gap beyond the PBE value of ~1 eV, eliminating low-energy metallic states and supporting robust semiconducting behavior.

Crucially, type-II band alignment arises from the offset between the conduction band minimum (CBM) of WS$_2$ and the valence band maximum (VBM) of WSe$_2$0. DFT and photoluminescence confirm spatially indirect exciton formation, with electrons preferentially localized in WS$_2$1 and holes in WSe$_2$2, yielding interlayer exciton emission at 1.45 eV at cryogenic temperature (Ghatak et al., 2019, Yuan et al., 2019).

2. Moiré Potentials, Twist Angle, and Lattice Reconstruction

A lattice constant mismatch (~4%) and arbitrary twist angle θ between layers establish a moiré superlattice characterized by a periodicity λ$_2$3 ≃ 8.5 nm near zero twist. Full atomistic relaxation using Stillinger–Weber and Kolmogorov–Crespi potentials yields commensurate moiré supercells for 0°, 60°, and other high-symmetry θ (Vitale et al., 2021, wang et al., 10 Feb 2026).

Under rotation, the moiré potential $_2$4 arises from three dominant mechanisms (wang et al., 10 Feb 2026):

1. Lattice Reconstruction: Nonuniform local strain $_2$5 shifts the conduction and valence band edges by Δ$_2$6 ($_2$7 eV for conduction, $_2$8 eV for valence), yielding modulations $_2$9200 meV for electrons and $_2$020 meV for holes.
2. Piezopotentials: Strain fields induce piezoelectric charges, leading to a screened potential $_2$1 of up to 90 meV for conduction and 40 meV for valence carriers.
3. Interlayer Charge Transfer: Stacking-dependent charge redistribution introduces a local vacuum-level difference $_2$2, modulating the band edges by $_2$3. This effect modulates conduction bands by 80 meV (R-type) or 40 meV (H-type).

The combined moiré potential localizes wavefunctions at specific high-symmetry stackings. R-type moiré traps both electrons and holes at the same dot (enhancing correlation), while H-type segregates them to different stackings (wang et al., 10 Feb 2026).

The table below summarizes energetics at key stackings calculated for high-symmetry registries:

Registry	Binding Energy (eV / cell)	Interlayer Spacing (Å)
0° c-term AA	–2.645	6.25
60° c-term AB	–1.906	6.46
0° m-term AA	+0.831	6.84

3. Band Structure, Spin–Orbit Effects, and Twist Engineering

First-principles calculations (HSE06+D3) for commensurate 0° and 60° stackings reveal (Lin et al., 2024):

• 0° stacking (AA′-like): Indirect gap $_2$4 eV (VBM at Γ, CBM at K), spin–orbit split Δ$_2$5 ≃ 0.45 eV at K.
• 60° stacking (AB′-like): Indirect gap $_2$6 eV (VBM at K, CBM at Λ), Δ$_2$7 ≃ 0.44 eV.
• Both registries retain type-II alignment; electrons and holes reside on WS$_2$8 and WSe$_2$9, respectively.

Unlike homobilayer WS$_2$0/WS$_2$1, where “magic” twist angles (17.9°, 42.1°) permit indirect-direct gap transitions and pronounced miniband flattening, WS$_2$2/WSe$_2$3 heterobilayers do not host direct band gaps in accessible stackings, and the gap variation with θ is modest (Lin et al., 2024). The valence band top remains K-derived for all θ ≥ 4.5°, and Γ-flat bands do not emerge (Vitale et al., 2021), precluding strongly correlated flat-band Mott physics accessible in certain homobilayer settings.

4. Moiré Exciton Physics and Many-Body Effects

WS$_2$4/WSe$_2$5 moiré heterobilayers support multiple classes of spatially modulated excitons, elucidated by combined GW-BSE theory and micro-reflection spectroscopy (Naik et al., 2022):

1. Modulated Wannier Excitons (Peak I, 1.68 eV): Highly localized (spatial width ∼2 nm), with strong oscillator strength ($_2$6) and binding energy ≈350 meV.
2. Intralayer Charge-Transfer Excitons (Peak III, 1.75 eV): Electron and hole are localized at separated moiré sites (Δr ≈ 5 nm), with reduced oscillator strength and binding ≈200 meV.
3. Mixed-Character Excitons (Peak II, 1.71 eV): Intermediate localization and oscillator strength.

Unique doping and magnetic field signatures further discriminate these states: Peak I is robust under hole doping, while Peak III red-shifts and broadens with electron filling, consistent with wavefunction spatial separation. These observations are beyond the reach of simplified continuum models and highlight the importance of atomically resolved moiré potentials in determining moiré exciton physics (Naik et al., 2022).

5. Exciton Transport, Diffusion Anomalies, and Quantum Emitter Prospects

Moiré potentials of depth $_2$7 ≈ 0.15 eV (0° twist) and 0.11 eV (60° twist), as quantified via DFT and transient absorption microscopy, localize interlayer excitons and generate twist-angle-dependent exciton superlattice arrays (Yuan et al., 2019). Exciton dynamics obey a generalized equation incorporating both Fickian diffusion and dipole-dipole repulsion:

$_2$8

with $_2$9. At low density and temperature, excitons exhibit sub-diffusive, spatially localized motion ($_2$0 nm at 1 ns for $_2$1 K); at high density, screening by dipolar interactions delocalizes them ($_2$2 up to 250 nm) (Yuan et al., 2019). Thus, twist angle and carrier density act as key control knobs for engineering excitonic quantum emitter lattices, single-photon sources, and long-range excitonic circuits.

Table: Diffusion Length and Moiré Trap Depth

Twist Angle (θ)	$_2$3 (meV)	$_2$4 (nm, 1 ns, 295 K)
0°	150	150
60°	110	200

K–Q interlayer excitons are energetically favored ($_2$5 = 60–90 meV over K–K) and constitute the true ground state (Yuan et al., 2019), compelling a revision of device concepts for interlayer excitonic devices.

6. Moiré Hubbard Lattices and Mott Physics

WS$_2$6/WSe$_2$7 heterobilayers realize moiré-trapped Bose–Fermi mixtures, enabling direct simulation of the Bose–Fermi–Hubbard model with unprecedented interaction strengths (Gao et al., 2023). Device architectures—60° H-stacked, dual-gated and hBN encapsulated—allow independent control of electron (WS$_2$8 CB miniband) and exciton filling.

Key results include:

• Observation of discrete photoluminescence lines (X₁, X₂) corresponding to single-exciton occupancy and double occupancy (electron + exciton or two excitons) per moiré site.
• Extraction of repulsion strengths: $_2$9 meV, $_2$0–29 meV.
• Exciton diffusion suppression under high filling (smoking-gun incompressibility), consistent with a Mott insulating state.

Effective low-energy physics is captured by a two-band Bose–Fermi–Hubbard Hamiltonian,

$_2$1

with strong interactions $_2$2. The phase diagram spans mixed-coherent, fermionic, and bosonic Mott insulating regions (Gao et al., 2023).

This platform extends quantum simulation of Hubbard models from ultracold atoms to solid state, facilitating studies of non-equilibrium, spinful, and driven many-body phases, including perspectives on Kondo lattice and quantum light emitters.

7. Future Prospects and Engineering Pathways

Twist angle and layer composition present powerful control means for engineering the electronic and excitonic landscape in WS$_2$3/WSe$_2$4 hetero-bilayers. Notably:

• Magic angles in homobilayer WS$_2$5 can realize flat-band and direct-gap transitions; in heterobilayers, analogous phenomena are implied but experimentally unresolved at present twist ranges (Lin et al., 2024).
• Tunability of moiré potential depth ($_2$6) and miniband widths by twist and gating/strain enables programmable exciton superlattices and prospects for Bose–Hubbard realization, Wigner crystallization, and designer quantum emitter arrays (Yuan et al., 2019, wang et al., 10 Feb 2026, Gao et al., 2023).
• Optical selection rules, valley/spin control, and stacking-dependent many-body regimes can be accessed by integrating gate-controlled charge doping and exciton injection, with device architectures suited for quantum communication, valleytronic circuits, and nonequilibrium condensate physics.

A plausible implication is that further reduction of twist angle or control of dielectric screening may yield more pronounced flat-band correlation, opening new routes to quantum simulation and device functionalities in the moiré landscape of TMD heterobilayers.