The file is a single workbook `Book1`/`Sheet1` where column A is the energy axis (eV) and columns 1–15 are already-merged spectra at powers 5→100 µW. No Si/reference folder exists, so I skip merging. Per protocol I despike first (mandatory before any fit), then run the Lorentzian power-dependence fit.



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

All 15 merged spectra (5–100 µW) were despiked (130 cosmic-ray spikes removed) and Lorentzian-fitted in the 1.23–1.45 eV window. Every column converged to a **HIGH**-tier fit (R² = 0.887–0.984, median 0.981), so all 15 points enter both panels.

**Extracted per-power values (Lorentzian 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 fit — intensity (left panel), I = a·Pᵇ**
- a = 16.28 ± 1.81
- **b = 1.0048 ± 0.0261**
- R² = 0.9951 → essentially linear (b ≈ 1), i.e. excitonic regime with no measurable saturation over 5–100 µW.

**Peak-position trend (right panel), xc vs P**
- Monotonic **blueshift** from 1.3206 eV (5 µW) to 1.3356 eV (100 µW): total shift ≈ 15 meV, ≈ 5 meV per decade of power. FWHM narrows correspondingly (62 → 54 meV), consistent with the weak power-induced blueshift.

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

RESULT: b = 1.0048 ± 0.0261