Massive Dirac Fermion Magnetoexciton in TMDCs

• Massive Dirac fermion magnetoexciton is a quasiparticle formed by an electron-hole pair in Dirac Landau levels of TMDCs, showcasing unique relativistic dispersion and spin‐orbit coupling.
• The phenomenon is characterized by quantized Landau levels, optical selection rules, and experimental signatures such as binding energy shifts and linewidth broadening under varying magnetic fields.
• Many-body interactions, including self-energy and vertex corrections, drive observable redshifts and polarization effects, providing insights into valley/spin ferromagnetism and excitonic behavior in quantum Hall states.

A massive Dirac fermion magnetoexciton is an emergent quasi-particle formed by the Coulomb-bound state of an electron and a hole, each residing in massive Dirac Landau levels, in two-dimensional semiconductors such as monolayer transition metal dichalcogenides (TMDCs) subjected to strong perpendicular magnetic fields. Unlike conventional magnetoexcitons in parabolic bands, these entities manifest unique spectroscopic signatures due to the relativistic (Dirac-like) dispersion, strong spin‐orbit coupling, multivalley structure, and quantization into Landau levels (LLs). Theoretical and experimental investigation of these magnetoexcitons directly probes many-body interactions, valley and spin polarization, and collective excitations in TMDC quantum Hall states (Sadecka et al., 17 Jan 2026, Kim et al., 2016).

1. Massive Dirac Fermion Model in TMDC Monolayers

The optoelectronic properties of monolayer TMDCs are dominated by carriers governed by a two-band massive Dirac fermion (mDF) Hamiltonian near the $K$ and $K'$ valleys. Including spin–orbit coupling (SOC), the low-energy Hamiltonian in the basis $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$ takes the form: $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$ where $v_F$ is the Fermi velocity, $\Delta$ the direct gap, $\lambda$ the SOC parameter, $\tau = \pm1$ the valley index, and $s = \pm1$ the spin. This Hamiltonian captures the direct band gap, valley and spin splitting, and the chiral nature of the Bloch states (Kim et al., 2016).

Experimental extraction of parameters from ARPES data, corroborated by tight-binding theory, yields (for MoS$_2$) values: $K'$0 eV, $K'$1 eV, $K'$2 m/s, effective mass $K'$3, lattice constant $K'$4 Å. Analogous parameter sets are obtained for MoSe$K'$5, WS$K'$6, and WSe$K'$7.

2. Landau Levels and Magnetic Quantization

In a perpendicular magnetic field $K'$8, minimal coupling leads to quantization of the massive Dirac spectrum into LLs. For $K'$9, the Hamiltonian becomes: $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$0 with $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$1 and magnetic length $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$2. Diagonalization yields LL energies for $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$3: $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$4 where $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$5 are effective $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$6-factors for electron and hole, and Zeeman effects are explicitly included. The $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$7 LL is a nondegenerate, highly spin and valley selective state dictated by the Dirac mass and Zeeman energies, a prominent feature in mDF systems (Sadecka et al., 17 Jan 2026, Kim et al., 2016).

3. Magnetoexciton Formation and Emission Selection Rules

A neutral massive Dirac magnetoexciton is a correlated excitation comprising an electron in CB-LL $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$8 and a hole in VB-LL $$\{\mathrm{CB}\downarrow,\mathrm{VB}\downarrow,\mathrm{CB}\uparrow,\mathrm{VB}\uparrow\}$$9, created on a reference ground state $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$0 as: $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$1 The many-body Hamiltonian for such configurations includes single-particle LL energies, electron–electron Coulomb interactions $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$2, and, if needed, a positive background subtraction. Solving the resultant Bethe–Salpeter–like eigenproblem in the restricted exciton subspace determines the correlated spectrum (Sadecka et al., 17 Jan 2026).

Optical transitions are governed by polarization- and valley-dependent selection rules: $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$3 for transitions between VB-LL~$$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$4 and CB-LL~$$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$5. Only such channels contribute to the magneto-photoluminescence (PL) spectra.

4. Many-Body Effects: Binding Energies, Redshift, and Broadening

The interplay between single-particle quantization and electron–hole interactions leads to distinct many-body phenomena:

• Self-energy corrections ($$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$6) cause a blue-shift in the pair continuum, scaling as $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$7.
• Vertex corrections (direct and exchange) provide the magnetoexciton binding (redshift), $$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$8. For MoS$$\hat H_{s,\tau}(\mathbf{q}) = v_F\left(\tau q_x \sigma_x + q_y \sigma_y\right) + \frac{\Delta}{2}\sigma_z + \frac{\lambda s \tau}{2}(I_2 - \sigma_z)$$9 at $v_F$0 T, $v_F$1 meV.
• State manifold broadening arises from configuration mixing, giving finite linewidths that increase with $v_F$2 and temperature. For a negatively charged trion ($v_F$3), the linewidth $v_F$4 meV at $v_F$5 K, decreasing to $v_F$6 meV at $v_F$7 K.
• Polarization: For a $v_F$8 quantum Hall ferromagnet in valley $v_F$9, emission is fully $\Delta$0 polarized, reflecting the underlying spin and valley polarization.

5. Quantitative Estimates and Experimental Comparison

Numerical diagonalization for monolayer MoS$\Delta$1 provides:

$\Delta$2 (T)	$\Delta$3 (meV)	$\Delta$4 (meV)	$\Delta$5 (meV)	$\Delta$6 (meV)
5	33	–11	–6	3
10	47	–16	–9	5
20	65	–23	–13	8

Here, $\Delta$7 is the neutral magnetoexciton binding energy, $\Delta$8 is the redshift for the exciton interacting with the $\Delta$9 state, $\lambda$0 the trion redshift, and $\lambda$1 the trion linewidth. All energies increase with $\lambda$2 as $\lambda$3 (Sadecka et al., 17 Jan 2026).

These quantitative features match experimental observations: $\lambda$4-polarized emission, redshifts of $\lambda$5–$\lambda$6 meV at $\lambda$7 T, and linewidth broadening of a few meV with rising $\lambda$8 or $\lambda$9 as in the reports by Finley (2021), Oreszczuk (2023), and Stier (2016).

6. Material Dependence and Model Parameters

Parameter sets for the mDF model across various TMDCs, established via ARPES and tight-binding calculations, are:

Material	$\tau = \pm1$0 [eV]	$\tau = \pm1$1 [eV]	$\tau = \pm1$2 [$\tau = \pm1$3 m/s]	$\tau = \pm1$4
MoS$\tau = \pm1$5	1.90	0.16	4.20	0.95
MoSe$\tau = \pm1$6	1.67	0.20	3.93	0.95
WS$\tau = \pm1$7	1.86	0.44	5.18	0.61
WSe$\tau = \pm1$8	2.04	0.48	4.88	0.75

These fundamental quantities determine LL quantization, magnetoexciton binding, and selection rules. At $\tau = \pm1$9 T, the single-particle 0→1 LL transition shift $s = \pm1$0 is $s = \pm1$1 meV (MoS$s = \pm1$2), $s = \pm1$3 meV (MoSe$s = \pm1$4), $s = \pm1$5 meV (WS$s = \pm1$6), and $s = \pm1$7 meV (WSe$s = \pm1$8), prior to inclusion of Coulomb binding (Kim et al., 2016).

7. Relevance and Outlook

The many-body treatment of massive Dirac fermion magnetoexcitons in monolayer TMDCs provides a quantitative framework for interpreting magneto-optical measurements at finite carrier density and high magnetic field. The ab initio–parametrized mDF model combined with exact diagonalization captures Landau level splitting, exciton and trion binding energies, field- and carrier-dependent redshifts, linewidths, and valley polarization (Sadecka et al., 17 Jan 2026). These findings establish magneto-spectroscopy as a direct probe of electronic correlations, valley/spin ferromagnetism, and nontrivial topology in atomically thin, strongly interacting, relativistic electron systems.