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

Pipeline executed end-to-end on `powerdep_manuscript.opj` (single merged sheet, 15 power levels 5–100 µW, energy axis in eV):

1. **Despike** (σ=4, cosmic-ray removal) — 130 spikes cleaned across all 15 columns before fitting.
2. **Lorentzian fit** on every spectrum (integration window 1.23–1.45 eV) — all 15 peaks converged to the HIGH tier (individual R² = 0.887–0.984, median 0.981).
3. Extracted **intensity** and **peak position xc** per power, then built both power-dependence panels and ran the allometric (power-law) fits.

**Left panel — Peak intensity vs power, I = a·Pᵇ**
- a = 16.28, **b = 1.0048 ± 0.0261**, R² = 0.995
- Interpretation: ~linear (b ≈ 1) → excitonic regime, no saturation over 5–100 µW.

**Right panel — Peak position xc vs power**
- Mild **blueshift**: xc rises from 1.3206 eV (5 µW) to 1.3356 eV (100 µW), a net ≈ 15 meV shift (~5 meV/decade).
- Allometric fit to the positions gives a near-zero exponent (b_xc ≈ 0.004, a ≈ 1.313) — i.e. the peak position is essentially power-independent within this range.

📥 Full per-power table (power, intensity, xc, FWHM, R²): [power_fit_results.csv]([file link removed])

RESULT: b = 1.0048 ± 0.0261