File inspected. Single merged workbook **Book1** with 15 power levels (5–100 µW), energy axis in eV, suggested fit range 1.23–1.45 eV. No separate Si folder — data is already merged.

Despiking first (mandatory before any fit):



Done. All 15 merged spectra (5–100 µW) were despiked, then each was fit with a **Lorentzian profile** over 1.23–1.45 eV. Every column converged at **HIGH** quality (R² = 0.887–0.984, median 0.981). The extracted peak intensities and centers were plotted against power and given allometric fits.

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

**Left — Peak intensity vs power** (allometric $I = a\,P^{b}$):
- $a = 16.28 \pm 1.81$ counts
- $b = 1.0048 \pm 0.0261$ → essentially **linear (excitonic)** scaling, R² = 0.995
- Intensity spans 99 (5 µW) → 1715 counts (100 µW)

**Right — Peak position (xc) vs power:**
- Blueshift from **1.3206 eV** (5 µW) to **1.3356 eV** (100 µW), i.e. ~15 meV total
- Trend slope ≈ 5.0 meV per decade of power (allometric fit)

Per-power centers run monotonically 1.3206 → 1.3356 eV; FWHM narrows slightly from ~62 → 54 meV as power rises.

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

RESULT: b = 1.005 ± 0.026