![Power Fit]([file link removed])

**Pipeline executed:** merged spectra (15 powers, 5–100 µW) → despiked (130 cosmic-ray spikes removed) → Lorentzian fit on every spectrum (R² = 0.887–0.984, all 15 HIGH-tier) → intensity & xc extracted → allometric fit I = a·Pᵇ.

**Allometric fit (intensity vs power):**
- **b = 1.0048 ± 0.0261**, a = 16.28 ± 1.81, R² = 0.995
- Interpretation: essentially **linear (excitonic)** — no saturation up to 100 µW.

**Peak position (xc) vs power:**
- Blueshifts from **1.3206 eV (5 µW) → 1.3356 eV (100 µW)**, i.e. ≈ +15 meV across the range (~5.0 meV/decade).
- FWHM narrows from ~62 meV (5 µW) to ~54 meV (100 µW).

Full per-power table (power, intensity, xc, FWHM, R²): [power_fit_results.csv]([file link removed])

RESULT: b = 1.0048 ± 0.0261