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

Pipeline: inspected `Book1` (15 merged spectra, 5→100 µW, one column each) → despiked (130 cosmic-ray spikes removed, σ=4, written back) → Lorentzian fit on every spectrum over $x \in [1.23,\,1.45]$ eV → allometric fit $I = a\,P^{b}$.

**Per-power extraction (all 15 columns, HIGH tier, $R^2$ 0.887–0.984):**

| Power (µW) | Intensity (counts) | xc (eV) | FWHM (meV) |
|---|---|---|---|
| 5 | 98.8 | 1.3206 | 62.4 |
| 10 | 189.9 | 1.3244 | 64.4 |
| 20 | 353.8 | 1.3278 | 62.7 |
| 40 | 694.6 | 1.3315 | 59.3 |
| 60 | 935.9 | 1.3332 | 58.0 |
| 80 | 1264.7 | 1.3345 | 56.3 |
| 100 | 1714.6 | 1.3356 | 54.4 |

**Allometric fit — Intensity vs Power (left panel):**
$$I = a\,P^{b},\quad a = 16.28,\quad b = 1.0048 \pm 0.0261,\quad R^2 = 0.9951$$
$b \approx 1$: intensity scales linearly with excitation power — an **excitonic (non-saturating)** regime.

**Peak position vs Power (right panel):**
xc rises monotonically from 1.3206 eV (5 µW) to 1.3356 eV (100 µW) — a **blueshift** of ~15 meV across the range (~+5.0 meV/decade). FWHM narrows slightly (62→54 meV) as power increases.

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

RESULT: b = 1.0048 ± 0.0261