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

## Results

**Lorentzian fits:** All 15 spectra (5–100 µW) fitted successfully in the window 1.23–1.45 eV. Every column reached **HIGH tier** quality (R² = 0.887–0.984, median 0.981).

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

**Allometric fit (I = a·P^b):**
- $a = 16.28 \pm 1.81$
- $b = 1.005 \pm 0.026$
- $R^2 = 0.995$

The exponent $b \approx 1.0$ indicates **linear (excitonic)** scaling of PL intensity with excitation power — no bi-exciton or Auger saturation regime is observed over this range.

**Peak position trend:** xc shifts from 1.3206 eV (5 µW) to 1.3356 eV (100 µW), a total blueshift of ~15 meV (~5 meV per decade of power), consistent with bandgap renormalization / exciton–exciton interactions at higher densities.

📥 [Download full results CSV]([file link removed])

RESULT: b = 1.005 ± 0.026