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

**Pipeline:** inspect → despike (130 cosmic‑ray spikes removed, σ=4) → Lorentzian fit on all 15 powers (x‑window 1.23–1.45 eV) → allometric fit I = a·Pᵇ.

**Lorentzian fit per power** (all HIGH tier, R² = 0.88–0.98):

| P (µW) | xc (eV) | FWHM (meV) | Intensity (cts) | R² |
|---|---|---|---|---|
| 5 | 1.3206 | 61.8 | 99 | 0.881 |
| 10 | 1.3244 | 64.4 | 190 | 0.950 |
| 20 | 1.3278 | 62.7 | 354 | 0.968 |
| 30 | 1.3300 | 61.3 | 516 | 0.975 |
| 40 | 1.3315 | 59.4 | 695 | 0.981 |
| 50 | 1.3328 | 59.0 | 818 | 0.981 |
| 70 | 1.3343 | 57.1 | 1140 | 0.983 |
| 100 | 1.3356 | 54.4 | 1707 | 0.984 |

**Allometric fit — intensity vs power** (left panel, I = a·Pᵇ):
- **b = 0.9973 ± 0.0254**, a = 16.82 ± 1.82, R² = 0.995 → essentially **linear (excitonic)** emission, no saturation over 5–100 µW.

**Peak position vs power** (right panel): monotonic **blue shift** from 1.3206 eV (5 µW) to 1.3356 eV (100 µW), total **+15 meV** (~5 meV/decade), following a saturating (log‑like) trend — the signature of band‑filling / Burstein–Moss. FWHM narrows slightly (62→54 meV) with power.

📥 Data: [power_fit_results.csv]([file link removed]) · [fit_results.csv]([file link removed])

RESULT: b = 0.9973 ± 0.0254