The file is a single merged workbook `Book1` with 15 power spectra (5–100 µW) plus the energy axis in column A. Suggested fit range is 1.23–1.45 eV. No Si reference folder exists, so the data is already merged. I'll despike (mandatory before any fit), then run the Lorentzian power-dependence fit.



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

## Lorentzian power-dependence analysis (Book1, 5–100 µW, window 1.23–1.45 eV)

All 15 spectra were despiked first (130 cosmic-ray spikes removed), then each was fitted with a **Lorentzian** profile. Every column reached the **HIGH** quality tier (R² = 0.887–0.984, median 0.981; SNR 7.6–24.0), so all 15 points enter both panels.

**Peak intensity vs power — allometric fit $I = a\,P^{b}$**
- $a = 16.28 \pm 1.81$ counts
- $b = 1.0048 \pm 0.0261$
- $R^2 = 0.995$ (15 points)
- Interpretation: **~linear (b ≈ 1)** → excitonic regime, no saturation or Auger quenching over 5–100 µW.

**Peak position (xc) vs power**
- Blueshifts monotonically from **1.3206 eV** (5 µW) to **1.3356 eV** (100 µW) — a total blue shift of ~15 meV (~+5.0 meV per decade of power).
- FWHM narrows from 62.4 meV → 54.4 meV across the same range.

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

RESULT: b = 1.0048 ± 0.0261