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

**Pipeline:** despike (σ=4, 130 spikes removed) → Lorentzian fit on all 15 power spectra (5–100 µW) → extract peak intensity & center per power → allometric fit on both panels.

**Lorentzian fits:** all 15 columns converged at HIGH quality (R² = 0.887–0.984, median 0.981; SNR 7.6–24.0).

**Intensity vs power** (left panel), allometric fit $I = a\,P^{b}$:
- $a = 16.28 \pm 1.81$
- $b = 1.0048 \pm 0.0261$, R² = 0.995 → **~linear (excitonic)** regime, no saturation over 5–100 µW.

**Peak position vs power** (right panel): monotonic **blueshift** from 1.3206 eV (5 µW) to 1.3356 eV (100 µW); total shift ≈ 15 meV (~5 meV/decade), consistent with exciton–exciton / band-filling blue-shift.

📥 Per-power data (intensity, xc, FWHM, R²): [power_fit_results.csv]([file link removed])

RESULT: b = 1.0048 ± 0.0261