File inspected: single workbook **Book1** with 16 columns — col A is the energy axis (eV, 1.196–1.606 eV) and columns 1–15 hold the 15 power spectra (5 → 100 µW). Data is already consolidated in one workbook (no separate Si/reference folder present), so no cross-folder merge is needed. Suggested fit window: **x_min = 1.23, x_max = 1.45 eV**.

Despiking is mandatory before any fit — running it now (σ=4, write-back).



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

Pipeline complete on **Book1** (15 power spectra, 5→100 µW). After despiking (130 spikes removed), each spectrum was Lorentzian-fitted in the 1.23–1.45 eV window; all 15 converged to **HIGH tier** (R² = 0.887–0.984, median 0.981), so every point is used in both panels.

**Allometric fits (I = a·P^b):**
- **Intensity vs power:** a = 16.28 ± 1.81, **b = 1.005 ± 0.026**, R² = 0.995 → essentially linear (excitonic, no saturation/biexciton quenching over this range).
- **Peak position vs power:** monotonic **blueshift** from 1.3206 eV (5 µW) to 1.3356 eV (100 µW); ≈ +5.0 meV/decade. FWHM narrows from ~62 → ~54 meV with increasing power.

**Per-power extracted values (Lorentzian fit):**

| Power (µW) | Intensity (counts) | xc (eV) |
|---:|---:|---:|
| 5 | 98.8 | 1.3206 |
| 10 | 189.9 | 1.3244 |
| 15 | 267.9 | 1.3264 |
| 20 | 353.8 | 1.3278 |
| 25 | 425.2 | 1.3290 |
| 30 | 516.0 | 1.3300 |
| 35 | 564.8 | 1.3304 |
| 40 | 694.6 | 1.3315 |
| 45 | 738.5 | 1.3318 |
| 50 | 818.3 | 1.3328 |
| 60 | 935.9 | 1.3332 |
| 70 | 1140.4 | 1.3343 |
| 80 | 1264.7 | 1.3345 |
| 90 | 1538.0 | 1.3351 |
| 100 | 1714.6 | 1.3356 |

📥 Full data: [power_fit_results.csv]([file link removed])

RESULT: b = 1.0048 ± 0.0261