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

**Pipeline:** inspect → despike (σ=4, 130 cosmic-ray spikes removed) → Lorentzian power-dependence fit (window 1.23–1.45 eV, screen-reader peak metric). All 15 powers (5–100 µW) converged to **HIGH tier** (Lorentzian R² = 0.887–0.984, median 0.981).

**Per-power results (Lorentzian fit):**

| Power (µW) | Intensity (counts) | xc (eV) | FWHM (meV) |
|---|---|---|---|
| 5 | 98.8 | 1.3206 | 62.4 |
| 10 | 189.9 | 1.3244 | 64.4 |
| 20 | 353.8 | 1.3278 | 62.7 |
| 30 | 516.0 | 1.3300 | 61.3 |
| 50 | 818.3 | 1.3328 | 59.0 |
| 70 | 1140.4 | 1.3343 | 57.1 |
| 100 | 1714.6 | 1.3356 | 54.4 |

(Full 15-point table: [power_fit_results.csv]([file link removed]))

**Allometric fits:**
- **Intensity vs power** — $I = a\,P^{b}$: $a = 16.28$, **$b = 1.0048 \pm 0.0261$**, R² = 0.9951 → essentially **linear (excitonic)** emission, no saturation or biexciton quenching across 5–100 µW.
- **Peak position vs power** — $x_c$ **blueshifts** monotonically from 1.3206 eV (5 µW) to 1.3356 eV (100 µW), a total shift of **~15.0 meV** (slope ≈ 5.0 meV/decade of power), consistent with bandgap renormalization / many-body effects at higher excitation.

RESULT: b = 1.0048 ± 0.0261