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

Pipeline: despike (σ=4, 130 spikes removed) → Lorentzian fit on all 15 spectra (5–100 µW) → allometric fits. All 15 powers landed in the **HIGH** tier (R² 0.887–0.984, median 0.981); none rejected.

**Per-power peak parameters (Lorentzian fit)** — representative values:

| P (µW) | Intensity (counts) | xc (eV) | FWHM (meV) | R² |
|---|---|---|---|---|
| 5 | 98.8 | 1.3206 | 62.4 | 0.887 |
| 25 | 425.2 | 1.3290 | 61.7 | 0.973 |
| 50 | 818.3 | 1.3328 | 59.0 | 0.981 |
| 75-ish | — | — | — | — |
| 100 | 1714.6 | 1.3356 | 54.4 | 0.984 |

(Full table for all 15 powers is in the CSV.)

**Plot 1 — Intensity vs power (allometric I = a·Pᵇ):**
- a = 16.28 ± 1.81, **b = 1.0048 ± 0.0261**, R² = 0.995
- Interpretation: ~linear (b ≈ 1) → single excitonic regime, no saturation or biexciton quenching over 5–100 µW.

**Plot 2 — Peak position vs power:**
- xc rises monotonically from 1.3206 eV (5 µW) to 1.3356 eV (100 µW) — a **blueshift** of ~15 meV across the range (~5.0 meV per decade of power). This is consistent with bandgap renormalization / carrier-density filling rather than thermal broadening.

📥 Full per-power data: [power_fit_results.csv]([file link removed])

RESULT: b = 1.0048 ± 0.0261