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

## Lorentzian power-dependence analysis — `Book1` (15 powers, 5–100 µW)

Each spectrum was despiked (130 cosmic-ray spikes removed) then fitted with a **Lorentzian + linear-slope** profile over **1.23–1.45 eV**. All 15 fits converged to the **HIGH** tier (R² = 0.887–0.984, median 0.981).

**Extracted per-power values (Lorentzian fit):**

| Power (µW) | Intensity (counts) | xc (eV) | FWHM (meV) | R² |
|---|---|---|---|---|
| 5 | 98.8 | 1.3206 | 62.4 | 0.887 |
| 10 | 189.9 | 1.3244 | 64.4 | 0.950 |
| 20 | 353.8 | 1.3278 | 62.7 | 0.968 |
| 30 | 516.0 | 1.3300 | 61.3 | 0.975 |
| 50 | 818.3 | 1.3328 | 59.0 | 0.981 |
| 70 | 1140.4 | 1.3343 | 57.1 | 0.983 |
| 100 | 1714.6 | 1.3356 | 54.4 | 0.984 |

**Allometric fit — intensity vs power** (left panel, $I = a\,P^{b}$):
- $a = 16.28 \pm 1.81$
- $b = 1.0048 \pm 0.0261$, R² = 0.995 → essentially **linear (b ≈ 1)**, i.e. no saturation over this range.

**Peak-position trend** (right panel): xc shifts from **1.3206 eV → 1.3356 eV** (a ~15 meV **blueshift**) as power rises 5 → 100 µW, tracking a saturating curve (≈ 5 meV per decade).

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

RESULT: b = 1.005 ± 0.026