The file is a single workbook (`Book1`) with one column per power level (5–100 µW, 15 powers), energy axis in eV, and `suggested_x_range = [1.23, 1.45]`. Data is already merged per power — no separate Si folder. Now I'll despike (mandatory before fitting), then run the Lorentzian power-dependence fit.



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

All 15 merged power spectra (5–100 µW) were despiked (130 cosmic-ray spikes removed), Lorentzian-fitted over 1.23–1.45 eV, and every column landed in the **HIGH** tier (R² = 0.887–0.984, median 0.981). Peak intensity and position per power:

| Power (µW) | Intensity (counts) | x_c (eV) | FWHM (meV) |
|---|---|---|---|
| 5 | 98.8 | 1.3206 | 62.4 |
| 10 | 189.9 | 1.3244 | 64.4 |
| 20 | 353.8 | 1.3278 | 62.7 |
| 30 | 516.0 | 1.3300 | 61.3 |
| 40 | 694.6 | 1.3315 | 59.3 |
| 50 | 818.3 | 1.3328 | 59.0 |
| 70 | 1140.4 | 1.3343 | 57.1 |
| 90 | 1538.0 | 1.3351 | 55.3 |
| 100 | 1714.6 | 1.3356 | 54.4 |

**Allometric fits:**
- **Intensity:** $I = a\,P^{b}$ → $a = 16.28 \pm 1.81$, **$b = 1.0048 \pm 0.0261$**, R² = 0.995 — essentially linear (excitonic regime, no saturation or Auger quenching over this range).
- **Peak position:** x_c rises monotonically from 1.3206 → 1.3356 eV with power (~5 meV/decade), i.e. a mild blue-shift; the dashed curve tracks it smoothly.

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

RESULT: b = 1.0048 ± 0.0261