Inspect shows a single workbook `Book1` with 15 power-level spectra (5–100 µW, columns B–P) on an energy axis (eV). No subfolders and no Si/background reference exist, so the per-power spectra are already consolidated — no merge step is needed. Suggested fit window is 1.23–1.45 eV.

Proceeding to the mandatory despike before fitting.






Despike complete — 130 cosmic-ray spikes removed across all 15 columns and written back. The spectra show a dominant exciton peak near ~1.32–1.33 eV. Now running the Lorentzian fit on every spectrum with cascade extraction of xc and intensity, plus the allometric fits on both panels.





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

**Pipeline:** 15 per-power spectra (5–100 µW) already consolidated in `Book1` → despiked (130 cosmic-ray spikes removed) → Lorentzian fit on every spectrum (window 1.23–1.45 eV) → allometric fits on both panels. All 15 columns passed at **HIGH** tier (R² = 0.887–0.984, median 0.981).

**Per-power Lorentzian results** (intensity from wing-baseline peak height; xc from fit):

| Power (µW) | Intensity (counts) | xc (eV) | FWHM (meV) | R² |
|---:|---:|---:|---:|---:|
| 5   | 98.8    | 1.32062 | 62.4 | 0.887 |
| 10  | 189.9   | 1.32437 | 64.4 | 0.950 |
| 15  | 267.9   | 1.32636 | 64.1 | 0.958 |
| 20  | 353.8   | 1.32783 | 62.7 | 0.968 |
| 25  | 425.2   | 1.32905 | 61.7 | 0.973 |
| 30  | 516.0   | 1.32998 | 61.3 | 0.975 |
| 35  | 564.8   | 1.33044 | 60.8 | 0.978 |
| 40  | 694.6   | 1.33152 | 59.3 | 0.981 |
| 45  | 738.5   | 1.33179 | 59.1 | 0.981 |
| 50  | 818.3   | 1.33279 | 59.0 | 0.981 |
| 60  | 935.9   | 1.33317 | 58.0 | 0.982 |
| 70  | 1140.4  | 1.33430 | 57.1 | 0.983 |
| 80  | 1264.7  | 1.33446 | 56.3 | 0.983 |
| 90  | 1538.0  | 1.33509 | 55.3 | 0.984 |
| 100 | 1714.6  | 1.33562 | 54.4 | 0.984 |

**Allometric fits:**
- **Intensity vs power** — I = a·P^b: a = 16.28 ± 1.81, **b = 1.0048 ± 0.0261**, R² = 0.995. b ≈ 1 → linear (excitonic) PL response, no saturation over 5–100 µW.
- **Peak position vs power** — monotonic **blueshift**: xc rises 1.3206 → 1.3356 eV (+15 meV), slope ≈ +5.0 meV/decade (band filling / bandgap renormalization signature); fitted curve shown in the right panel.

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

RESULT: b = 1.0048 ± 0.0261