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

Done. All 15 power spectra (5–100 µW) were Lorentzian-fitted after despiking; every column landed in the **HIGH** tier (R² = 0.887–0.984, median 0.981; FWHM ≈ 54–64 meV), so all 15 points went into both panels.

**Peak intensity vs power — allometric fit I = a·Pᵇ**
- b = 1.0048 ± 0.0261, a = 16.28 ± 1.81, R² = 0.995
- Interpretation: essentially **linear (excitonic)** — no saturation or biexciton recombination over this range.

**Peak position (xc) vs power — allometric/trend fit**
- Blueshift from **1.3206 eV (5 µW) → 1.3356 eV (100 µW)**, ≈ 15 meV total, ~5 meV/decade.

Per-power extraction (intensity / xc): 5 µW → 98.8 / 1.3206 eV · 30 µW → 516 / 1.3300 eV · 60 µW → 936 / 1.3332 eV · 100 µW → 1715 / 1.3356 eV (full table in the CSV).

📥 [power_fit_results.csv]([file link removed])

RESULT: b = 1.005 ± 0.026