Valleytronics Transition Metal Dichalcogenides Tmdc Berry Curvature Valley Hall Effect

# Valleytronics in 2D Transition Metal Dichalcogenides: Valley Pseudospin, Berry Curvature, and Optical Selection Rules

## 1. Introduction to Valleytronics

Valleytronics exploits the valley degree of freedom—distinct, inequivalent minima in the electronic band structure—as a platform for information processing and quantum control, analogous to spintronics but with valley pseudospin replacing electron spin. In 2D transition metal dichalcogenides (TMDCs) such as MoS₂, WSe₂, and WS₂, the valence band maximum and conduction band minimum are located at the K and K' points at opposite corners of the hexagonal Brillouin zone. These two valleys represent nonequivalent but energetically degenerate quantum states, and electrons can be selectively excited into one valley via circularly polarized light (optical circular dichroism). The combination of valley degeneracy, Berry curvature effects, and strong spin-orbit coupling creates a rich platform for manipulating quantum information and engineering novel topological properties.

This article develops the theory of valleytronics from first principles, covering the low-energy Dirac Hamiltonian at K and K' valleys, Berry curvature and anomalous velocity, optical selection rules for circularly polarized light, the valley Hall effect, and valley coherence and relaxation. A companion Python solver computes Berry curvature maps, anomalous velocity fields, and valley Hall conductivity tensors.

## 2. Crystal Structure and Electronic Dispersion in TMDCs

Monolayer TMDCs possess a hexagonal lattice with D₃h point-group symmetry. Each unit cell contains one transition metal atom (Mo or W) sandwiched between two chalcogen atoms (S or Se), forming a trigonal prismatic coordination. The band structure near the Fermi level consists of d-orbitals of the metal and p-orbitals of the chalcogen, hybridized by strong spin-orbit coupling.

The conduction band minimum (CBM) and valence band maximum (VBM) are located at the K and K' points (inequivalent corners of the hexagonal Brillouin zone), separated by a reciprocal lattice vector Q = K - K'. Time-reversal symmetry relates band structures at K and K': E(mathbfk, K) = E(-mathbfk, K') for states of opposite spin-valley quantum numbers.

## 3. Low-Energy Dirac Hamiltonian at K and K' Valleys

Near the band edges at the K and K' points, the electronic structure can be described by an effective Dirac Hamiltonian. At the K valley, the low-energy Hamiltonian (within 2×2 dimension for spin-up conduction and valence states) is:

where:
- mathbfk = (k_x, k_y) is the deviation from the K point
- sigma_x,y,z are Pauli matrices for band/orbital index (conduction vs valence)
- tau_z is the valley pseudospin operator (±1 for K/K')
- s_z is the true electron spin (up/down)
- v_F is the Fermi velocity (~10⁶ m/s for MoS₂)
- Delta is the band gap (1.6–2.8 eV for monolayer TMDCs)
- lambda is the spin-orbit coupling strength (10–100 meV)

The key feature is the term lambda tau_z s_z: spin-orbit coupling couples valley pseudospin to true electron spin. This spin-valley locking means spin-up electrons preferentially occupy one valley (say K), while spin-down electrons occupy the other (K'). This coupling is a consequence of the broken inversion symmetry in monolayer TMDCs and is absent in graphene.

At the K' valley, the Hamiltonian transforms as:

reflecting the combined effect of time-reversal symmetry and the inversion-breaking crystal structure.

## 4. Berry Curvature and Topological Properties

The electronic band structure in TMDCs exhibits nontrivial Berry curvature—a geometric property of the quantum mechanical wavefunction that acts like a fictitious magnetic field in reciprocal space. For a band n with Bloch wavefunction |psi_n(mathbfk)rangle, the Berry connection is:

and the Berry curvature is:

In 2D systems, the z-component Omega_z(mathbfk) is most relevant. For a simple 2D Dirac Hamiltonian (as in TMDCs near the K/K' points), the Berry curvature exhibits a striking monopole-like structure: large positive Omega_z near the conduction band minimum at K, and large negative Omega_z at K'. The integrated Berry curvature (Berry phase) around a closed loop enclosing a band touching point is π, signifying a topological winding number of ±1/2 for each valley-band pair.

## 5. Anomalous Velocity and Anomalous Hall Effect

The Berry curvature modifies the electron dynamics in external electric fields. The equation of motion for an electron in band n under electric field mathbfE becomes:

where the second term in the first equation is the anomalous velocity, proportional to the Berry curvature crossed with the applied force. This term arises from the geometric phase of the wavefunction and has no counterpart in classical mechanics.

For a finite density of carriers in a given band, the anomalous velocity produces a transverse current perpendicular to the applied field—the anomalous Hall effect (AHE):

where f(E) is the Fermi-Dirac distribution. In TMDCs with filled or partially filled valleys, this can produce sizable anomalous Hall conductivity.

## 6. Optical Selection Rules and Circular Dichroism

The optical response of TMDCs to circularly polarized light is governed by dipole selection rules that emerge from band symmetry. In the convention used here, right-handed circular polarization selectively drives the direct valence-to-conduction transition within K, with dipole matrix element:

where M is the dipole moment amplitude. Photon momentum is negligible compared with the valley separation, so momentum conservation keeps the transition within the same valley:

The orbital angular momentum of the photon combines with the valley-dependent Bloch orbital moment to select a specific same-valley interband transition.

For left-handed circular polarization, the selection is reversed and the direct transition occurs within K′. This optical valley polarization is a direct probe of valley physics: one helicity preferentially excites K excitons and the opposite helicity excites K′ excitons, subject to the stated convention.

The optical circular dichroism (the difference in absorption for sigma^+ versus sigma^-) quantifies this valley selectivity:

where the valley Hall effect arises from the spin-orbit coupling–induced splitting of the valence band into spin-polarized subbands at each valley.

## 7. Valley Hall Effect and Valley Hall Conductivity

The valley Hall effect is a transverse transport phenomenon where a charge current driven along one direction produces a valley current (net flow of electrons in one valley minus the other) perpendicular to the charge current. The valley Hall conductivity tensor is:

summed over the conduction band contributions from K and K' valleys. Since Omega_z has opposite signs at K and K' (monopole structure), the valley Hall conductivity can be substantial. Typical values for monolayer TMDC systems at room temperature are sigma_xy^textvalley sim 10^2 to 10^3 Ω⁻¹.

The valley Hall effect manifests experimentally as valley-dependent transverse voltages and can be exploited to separate and manipulate electrons in different valleys. This forms the basis for valleytronic devices: valley filters (structures that preferentially transmit one valley), valley modulators (electrically or optically switching between valleys), and valley-based logic gates.

## 8. Valley Coherence and Valley Decoherence

While valleys are robust quantum numbers protected by inversion and time-reversal symmetry, real materials break these symmetries through disorder, strain, and lattice defects. Valley-impurity scattering—scattering that mixes the K and K' valleys—causes valley coherence to decay. The typical valley decoherence time T_v at room temperature is in the picosecond range (a few picoseconds to tens of picoseconds), limiting the timescale on which valley information can be maintained.

Sources of valley decoherence include:

1. Disorder-induced scattering: Impurities and vacancies couple valleys via nonmagnetic scattering.
2. Lattice vibrations (phonons): Acoustic and optical phonons can scatter electrons between valleys, with coupling strength depending on electron-phonon matrix elements.
3. Magnetic moments: Paramagnetic impurities or substrate-induced magnetism breaks time-reversal symmetry.
4. Substrate effects: Dielectric environment of the substrate modulates the effective band structure and can suppress or enhance valley splitting.

Strategies to extend valley coherence include:
- Encapsulation in hBN (hexagonal boron nitride), which reduces disorder and provides a cleaner electrostatic environment
- Resonant excitation tuned to valley-active exciton lines
- Application of perpendicular magnetic fields to lift valley-orbit coupling

## 9. Valley Excitons and Exciton Valley Polarization

When an electron is excited to the conduction band via photon absorption, it leaves behind a hole in the valence band. The electron-hole pair forms a Coulomb-bound state—an exciton. In TMDCs, valley selection rules apply to exciton formation: sigma^+ photons create excitons with electrons in the K valley and holes in the K valley (valley-matched excitons), while sigma^- photons create K'-valley excitons.

These valley excitons are the primary optical excitations in TMDCs and exhibit strong valley polarization under circularly polarized excitation. The exciton binding energy is exceptionally large (~300–500 meV) due to the reduced dielectric screening in 2D (the dielectric constant ε_r ~ 4–6 is much lower than in bulk semiconductors ε_r ~ 10). This makes excitons highly stable even at room temperature.

Valley-exciton dynamics—their formation, diffusion, decay, and depolarization—are rich research topics. Exciton valley polarization can be probed via time-resolved circular dichroism or valley-sensitive Kerr rotation measurements.

## 10. Spin-Valley Locking and Coupled Spin-Valley Physics

The spin-orbit coupling term lambda tau_z s_z in the TMDC Hamiltonian creates spin-valley locking: the spin angular momentum of an electron is correlated with its valley index. Specifically, at the K valley, spin-up electrons have lower energy than spin-down electrons (spin-up state in K is the top of the valence band). At K', the situation reverses due to time-reversal symmetry.

This spin-valley coupling has profound implications:

1. Selection rules: Optical transitions must conserve both spin and valley, leading to strong circular dichroism.
2. Spin Hall effect and valley Hall effect: The spin Hall effect and valley Hall effect are intimately related; the Rashba-like spin-orbit term couples them.
3. Spin relaxation: Spin-valley locking suppresses certain spin-relaxation channels that would be active in conventional semiconductors. This can lead to remarkably long spin relaxation times in some TMDC systems.

## 11. Numerical Solver: Berry Curvature and Valley Hall Conductivity

The following Python code computes Berry curvature distributions and valley Hall conductivity:

import numpy as np
import matplotlib.pyplot as plt
from scipy.linalg import eigh
from matplotlib.colors import SymLogNorm

def dirac_hamiltonian_k_valley(kx, ky, vF, Delta, lambda_so):
    """
    Low-energy Dirac Hamiltonian at K valley.
    2x2 matrix: [conduction, valence] x [up, down] spin
    Simplified: ignoring explicit spin index, considering spin-up at K valley.
    """
    H = np.array([
        [Delta/2, vF * (kx - 1j*ky)],
        [vF * (kx + 1j*ky), -Delta/2]
    ], dtype=complex)
    return H

def compute_berry_curvature(kx, ky, vF, Delta, lambda_so, dk=0.01):
    """
    Numerically compute Berry curvature at (kx, ky).
    Using finite differences of Berry connection.
    """
    # Compute eigenstates at kx, ky
    H_center = dirac_hamiltonian_k_valley(kx, ky, vF, Delta, lambda_so)
    E_c, V_c = eigh(H_center)
    # V_c[:, 0] is the conduction band eigenstate (lower energy)
    psi_c = V_c[:, 0]
    
    # Compute Berry connection components numerically
    # A_x = <psi | i d/dx psi>
    H_dx = dirac_hamiltonian_k_valley(kx + dk, ky, vF, Delta, lambda_so)
    E_dx, V_dx = eigh(H_dx)
    psi_dx = V_dx[:, 0]
    
    # Fix phase to minimize |psi_dx - e^{i*alpha} psi_c|
    phase_x = np.angle(np.dot(np.conj(psi_c), psi_dx))
    psi_dx_aligned = psi_dx * np.exp(-1j * phase_x)
    
    A_x = np.imag(np.dot(np.conj(psi_c), (psi_dx_aligned - psi_c) / dk))
    
    # Similarly for y-direction
    H_dy = dirac_hamiltonian_k_valley(kx, ky + dk, vF, Delta, lambda_so)
    E_dy, V_dy = eigh(H_dy)
    psi_dy = V_dy[:, 0]
    
    phase_y = np.angle(np.dot(np.conj(psi_c), psi_dy))
    psi_dy_aligned = psi_dy * np.exp(-1j * phase_y)
    
    A_y = np.imag(np.dot(np.conj(psi_c), (psi_dy_aligned - psi_c) / dk))
    
    # Omega_z = dA_y/dx - dA_x/dy
    H_dxy = dirac_hamiltonian_k_valley(kx + dk, ky + dk, vF, Delta, lambda_so)
    E_dxy, V_dxy = eigh(H_dxy)
    psi_dxy = V_dxy[:, 0]
    
    dA_y_dx = (A_y - np.imag(np.dot(np.conj(psi_c), (psi_dx_aligned - psi_c) / dk))) / dk
    dA_x_dy = (A_x - np.imag(np.dot(np.conj(psi_c), (psi_dy_aligned - psi_c) / dk))) / dk
    
    Omega_z = dA_y_dx - dA_x_dy
    
    return Omega_z, A_x, A_y

def compute_valley_hall_conductivity(kx_range, ky_range, vF, Delta, lambda_so, T=300):
    """
    Compute valley Hall conductivity from Berry curvature.
    """
    kx_vals = np.linspace(-kx_range, kx_range, 50)
    ky_vals = np.linspace(-ky_range, ky_range, 50)
    
    sigma_valley_xy = 0
    
    for kx in kx_vals:
        for ky in ky_vals:
            try:
                Omega_z, _, _ = compute_berry_curvature(kx, ky, vF, Delta, lambda_so, dk=0.005)
                # Fermi-Dirac factor (simplified)
                E_k = np.sqrt((vF * np.sqrt(kx**2 + ky**2))**2 + (Delta/2)**2)
                f_k = 1.0 / (1.0 + np.exp(E_k / (0.025 * 300)))  # k_B T ~ 0.025 meV at 300K
                
                sigma_valley_xy += -f_k * Omega_z
            except:
                pass
    
    return sigma_valley_xy

def plot_valley_physics(vF=1.0, Delta=1.6, lambda_so=0.05):
    """
    Plot Berry curvature, anomalous velocity, and related quantities.
    """
    # Grid for contour plots
    kx_vals = np.linspace(-0.3, 0.3, 100)
    ky_vals = np.linspace(-0.3, 0.3, 100)
    KX, KY = np.meshgrid(kx_vals, ky_vals)
    
    Omega_z_map = np.zeros_like(KX)
    
    for i, kx in enumerate(kx_vals):
        for j, ky in enumerate(ky_vals):
            try:
                Omega_z, _, _ = compute_berry_curvature(kx, ky, vF, Delta, lambda_so, dk=0.01)
                Omega_z_map[j, i] = Omega_z
            except:
                Omega_z_map[j, i] = 0
    
    fig, axes = plt.subplots(2, 2, figsize=(14, 12))
    
    # Panel 1: Berry curvature at K valley
    ax = axes[0, 0]
    norm = SymLogNorm(linthresh=0.01, vmin=-0.5, vmax=0.5)
    contour = ax.contourf(KX, KY, Omega_z_map, levels=30, cmap='RdBu_r', norm=norm)
    ax.contour(KX, KY, Omega_z_map, levels=[0], colors='black', linewidths=1)
    cbar = plt.colorbar(contour, ax=ax, label='Berry curvature Ω_z')
    ax.set_xlabel('k_x (1/Å)', fontsize=11)
    ax.set_ylabel('k_y (1/Å)', fontsize=11)
    ax.set_title('Berry Curvature Monopole Structure at K Valley', fontsize=12)
    
    # Panel 2: Band dispersion
    ax = axes[0, 1]
    k_vals = np.linspace(0, 0.4, 100)
    E_c = np.sqrt((vF * k_vals)**2 + (Delta/2)**2)
    E_v = -np.sqrt((vF * k_vals)**2 + (Delta/2)**2)
    
    ax.plot(k_vals, E_c, 'b-', linewidth=2.5, label='Conduction band')
    ax.plot(k_vals, E_v, 'r-', linewidth=2.5, label='Valence band')
    ax.fill_between(k_vals, E_v, 0, alpha=0.2, color='red', label='Filled states (T=0)')
    ax.set_xlabel('|k| (1/Å)', fontsize=11)
    ax.set_ylabel('Energy (eV)', fontsize=11)
    ax.set_title('TMDC Band Dispersion Near K Point', fontsize=12)
    ax.legend(fontsize=10)
    ax.grid(True, alpha=0.3)
    ax.axhline(0, color='k', linewidth=0.8)
    
    # Panel 3: Optical selection rules
    ax = axes[1, 0]
    photon_types = ['σ⁺ (RCP)', 'σ⁻ (LCP)']
    valley_populations = [0.9, 0.1]  # K and K' valley occupations
    
    bars = ax.bar(photon_types, valley_populations, color=['blue', 'red'], alpha=0.7, edgecolor='black', linewidth=2)
    ax.set_ylabel('Valley polarization', fontsize=11)
    ax.set_title('Optical Circular Dichroism: Valley-Selective Excitation', fontsize=12)
    ax.set_ylim([0, 1])
    for i, (photon, pop) in enumerate(zip(photon_types, valley_populations)):
        ax.text(i, pop + 0.05, f'{pop:.0%}', ha='center', fontsize=11, fontweight='bold')
    ax.grid(True, alpha=0.3, axis='y')
    
    # Panel 4: Summary of valley physics
    ax = axes[1, 1]
    ax.axis('off')
    summary_text = f"""
    TMDC Valleytronics Summary

    Crystal structure: Hexagonal monolayer
    Transition metal: Mo or W
    Chalcogen: S or Se

    Valley locations:
      K point: (4π/3a, 0)
      K' point: (-4π/3a, 0)

    Band gap: Δ = {Delta:.2f} eV

    Low-energy physics:
    • Dirac-like dispersion at K, K'
    • Berry curvature monopoles
    • Spin-valley locking
    • Time-reversal symmetry

    Optical selection rules:
    • σ⁺ → K valley (circular dichroism)
    • σ⁻ → K' valley

    Valley Hall Effect:
    • Transverse valley current
    • σ_xy^valley ~ 100 Ω⁻¹

    Applications:
    ✓ Valley-based qubits
    ✓ Valleytronic switches
    ✓ Chiral optical modulators
    ✓ Topological electronics
    
    Valley decoherence time:
      T_v ~ 1-10 ps (room T)
    """
    ax.text(0.05, 0.95, summary_text, transform=ax.transAxes, fontsize=10,
            verticalalignment='top', family='monospace',
            bbox=dict(boxstyle='round', facecolor='lightyellow', alpha=0.8))
    
    plt.tight_layout()
    plt.show()

# Parameters for MoS2
vF_MoS2 = 3.6e5  # Fermi velocity in m/s -> normalized to 1.0
Delta_MoS2 = 1.66  # Band gap in eV
lambda_so_MoS2 = 0.052  # Spin-orbit coupling in eV

print("TMDC Valleytronics Analysis")
print(f"Material: MoS₂")
print(f"Band gap: {Delta_MoS2:.2f} eV")
print(f"Fermi velocity: ~{vF_MoS2/1e5:.1f} × 10⁵ m/s")
print(f"Spin-orbit coupling: {lambda_so_MoS2*1000:.1f} meV")

plot_valley_physics(vF=1.0, Delta=Delta_MoS2, lambda_so=lambda_so_MoS2)

## 12. Exciton Valley Physics and Brightdark Exciton Splittings

In direct bandgap TMDCs (MoS₂, WSe₂), the lowest-energy exciton states are dark excitons—the optically forbidden 1s triplet state (parallel spins on electron and hole)—and bright excitons—the optically allowed 1s singlet state (antiparallel spins). The exchange interaction between the electron and hole creates an energy splitting of 1–10 meV between bright and dark excitons, depending on material and dielectric environment. Dark excitons have exceptionally long lifetimes (nanoseconds) compared to bright excitons (picoseconds), making them candidates for long-lived quantum memories in valleytronic systems.

## 13. Valleytronics Beyond Monolayers

While monolayer TMDCs provide the most direct platform for valleytronics, valley physics also emerges in:

1. Bilayer TMDCs: Moiré patterns in twisted bilayers create flat bands and tunable valley polarization.
2. Heterostructures: TMDCs on graphene or other 2D materials inherit or modify valley properties.
3. Topological insulators: Valley-momentum locking in topological edge states offers protection against certain types of disorder.

## 14. Experimental Probes and Detection

Valley physics is probed experimentally via:

• Time-resolved Kerr rotation: Monitors valley polarization decay
• Circular dichroism spectroscopy: Measures valley-selective optical absorption
• Transverse voltage measurements: Detects valley Hall effect
• Photoluminescence under circularly polarized excitation: Directly images valley-exciton population

## 15. Future Directions and Applications

### 15.1 Valleytronic Devices
Valley-based transistors, modulators, and memory devices operating at lower power than conventional electronics.

### 15.2 Quantum Information
Valley qubits encoded in superpositions of K and K' states, with long coherence times potentially exceeding conventional spin qubits.

### 15.3 Topological Protection
Valley topology offers protection for certain transport phenomena against disorder, enabling robust quantum devices.

## Conclusion

Valleytronics in 2D TMDCs exploits the valley degree of freedom—nonequivalent band extrema in momentum space—as a platform for quantum information and novel transport phenomena. Berry curvature monopoles at the K and K' points endow conduction electrons with anomalous velocity, producing transverse Hall currents (valley Hall effect). Optical circular dichroism enables valley-selective excitation: right-handed photons preferentially excite electrons into the K valley, while left-handed photons excite K'. Spin-valley locking couples the true spin to valley pseudospin, creating rich opportunities for manipulating both degrees of freedom simultaneously. The presented numerical solver demonstrates Berry curvature mapping and valley Hall conductivity calculation. TMDCs continue to unveil new valleytronic phenomena and remain at the frontier of 2D materials physics.

## 16. Unified TMDC Hamiltonian and Valley Geometry

With lattice scale a, valley label τ_z = +1 at K and −1 at K′, orbital Pauli matrices σ, and spin s_z, the requested low-energy model is

The last term is algebraically λτ_z(1−σ_z)s_z/2 and primarily splits the valence band in this minimal basis. Time reversal changes k, τ_z, and spin together, preserving equal total energy at the paired time-reversed states.

For the corresponding massive Dirac sector, the stated valley-contrasting curvature is

Opposite signs at K and K′ cause charge Hall currents to cancel when the valleys are equally populated, while the valley-weighted transverse current adds.

Wave-packet dynamics in band n obeys

The second term is the anomalous velocity. An in-plane electric field therefore separates opposite valley populations toward opposite sample edges.

## 17. Optical Circular Dichroism Without Intervalley Momentum Transfer

Because photon momentum is negligible on the Brillouin-zone scale, a direct optical transition is intravalley. In the usual convention, right-circular σ+ light addresses one valley, commonly K, while left-circular σ− light addresses K′; reversing convention for helicity or layer orientation reverses the labels but not the one-to-one selection.

The transition does not move an electron from K to K′. Instead, the orbital angular momentum of the Bloch states and photon helicity determine which same-valley interband matrix element is nonzero. Intervalley scattering requires large momentum supplied by a defect, edge, or suitable phonon.

Excitons dominate optical spectra in monolayer TMDCs. Measured circular polarization therefore reflects excitation selectivity together with intervalley exchange, phonon scattering, recombination, and dark-state conversion. Optical dichroism alone is not a direct measurement of dc valley Hall conductivity.

## 18. Python Berry-Curvature and Valley-Hall Solver

The following model evaluates the requested curvature in both valleys, plots the momentum-space landscapes, and integrates a valley-weighted Hall response up to a tunable Fermi wave vector.

import numpy as np
import matplotlib.pyplot as plt
from scipy.constants import hbar, elementary_charge as e

a = 0.319e-9
vF = 5.3e5
Delta = 1.75*e
kmax = 6.0e9
grid = np.linspace(-kmax, kmax, 501)
kx, ky = np.meshgrid(grid, grid, indexing="xy")
k2 = kx**2+ky**2

def berry_curvature(k_squared, tau):
    numerator = tau*2*a**2*Delta*hbar**2*vF**2
    denominator = (4*hbar**2*vF**2*k_squared+Delta**2)**1.5
    return numerator/denominator

Omega_K = berry_curvature(k2, +1)
Omega_Kp = berry_curvature(k2, -1)

# Radial integration of the valley-weighted intrinsic conductivity.
kr = np.linspace(0.0, kmax, 5000)
Omega_radial = berry_curvature(kr**2, +1)
kF_values = np.linspace(0.0, kmax, 300)
sigma_v = []
for kF in kF_values:
    use = kr <= kF
    integral_K = np.trapz(2*np.pi*kr[use]*Omega_radial[use], kr[use])/(2*np.pi)**2
    # Valley weighting tau*Omega makes K and K' contributions add.
    integral_valley = 2*integral_K
    sigma_v.append(-(e**2/hbar)*integral_valley)
sigma_v = np.asarray(sigma_v)
conductance_quantum = e**2/(2*np.pi*hbar)

# Anomalous velocity field for an electric field along +x.
Ex = 1.0e4
v_anom_K_y = e*Ex*Omega_K/hbar

fig, axes = plt.subplots(2, 2, figsize=(12, 10))
extent = [grid[0]/1e9, grid[-1]/1e9, grid[0]/1e9, grid[-1]/1e9]
im0 = axes[0, 0].imshow(Omega_K*1e18, origin="lower", extent=extent)
axes[0, 0].set_title("Berry curvature at K (nm²)")
fig.colorbar(im0, ax=axes[0, 0])
im1 = axes[0, 1].imshow(Omega_Kp*1e18, origin="lower", extent=extent)
axes[0, 1].set_title("Berry curvature at K′ (nm²)")
fig.colorbar(im1, ax=axes[0, 1])

axes[1, 0].plot(kF_values/1e9, sigma_v/conductance_quantum)
axes[1, 0].set(xlabel="Fermi wave vector (nm⁻¹)",
               ylabel="Valley Hall conductivity / (e²/h)",
               title="Intrinsic valley Hall response")

mid = grid.size//2
axes[1, 1].plot(grid/1e9, v_anom_K_y[mid, :], label="K")
axes[1, 1].plot(grid/1e9, -v_anom_K_y[mid, :], label="K′")
axes[1, 1].set(xlabel="kx (nm⁻¹)", ylabel="Transverse anomalous velocity (m/s)",
               title="Opposite valley deflection")
axes[1, 1].legend()
for ax in axes.flat:
    ax.set_xlabel(ax.get_xlabel() or "kx (nm⁻¹)")
    ax.set_ylabel(ax.get_ylabel() or "ky (nm⁻¹)")
    ax.grid(alpha=0.25)
plt.tight_layout()
plt.show()

This continuum integral is valid near the valleys and should not be extended across the entire Brillouin zone. Quantitative comparison requires the correct band sign, spin branch, Fermi occupation, excitonic effects for optical experiments, disorder and side-jump terms where relevant, and a clear convention for valley current.

## 19. Transport and Optical Verification

A strong valley Hall experiment distinguishes edge accumulation from ordinary charge Hall signals, contact asymmetry, strain gradients, and thermoelectric effects. Useful controls include helicity reversal, electric-field reversal, spatially resolved edge emission, time-resolved depolarization, and comparison of monolayer with inversion-symmetric multilayer controls.

Report the helicity convention, crystal orientation, excitation energy relative to A and B excitons, carrier density, temperature, intervalley lifetime, device geometry, contact resistance, and whether conductivity is charge-, spin-, or valley-weighted. These definitions prevent sign and factor-of-two ambiguities.