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

## Power-dependence analysis complete

**Pipeline:** inspect → despike (130 cosmic-ray spikes removed, σ=4) → Lorentzian fit on all 15 spectra → allometric fits on both panels.

**Lorentzian fits (per power):** All 15 columns (5–100 µW) converged at **HIGH tier**, R² = 0.887–0.984 (median 0.981). Extracted per spectrum:

| Quantity | Range across power |
|---|---|
| Peak intensity | 99 → 1715 counts |
| Peak position xc | 1.3206 → 1.3356 eV |
| FWHM | 62 → 54 meV |

**Allometric fits (I = a·Pᵇ):**

- **Intensity panel (left):** $I = a\cdot P^{b}$, with $a = 16.28 \pm 1.81$, $b = 1.0048 \pm 0.0261$, R² = 0.995. The exponent b ≈ 1 means the PL intensity scales **linearly** with excitation power — an excitonic regime with no saturation or bimolecular recombination over 5–100 µW.
- **Peak-position panel (right):** xc **blueshifts** monotonically from 1.3206 eV (5 µW) to 1.3356 eV (100 µW), a total shift of ~15 meV (~5 meV/decade), consistent with thermal broadening / band-filling as carrier density rises.

📥 **[power_fit_results.csv]([file link removed])** — full per-power table (intensity, xc, FWHM, R², SNR, tier).

RESULT: b = 1.0048 ± 0.0261