The file contains one workbook (`Book1`) with the x-axis in eV (range 1.196–1.606 eV) and 15 power columns from 5 µW to 100 µW. Suggested fit window is **1.23–1.45 eV**. The spectra are already organized one-per-power, so I'll despike (mandatory before fitting) and then run the Lorentzian power-dependence analysis.



Despiking removed 130 cosmic-ray spikes across all 15 columns (written back). The main emission sits near **1.32 eV**, well inside the 1.23–1.45 eV window. Now running the Lorentzian power-dependence analysis with both panels (intensity and peak position vs power, each with an allometric fit).



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

**Pipeline:** 15 power-resolved spectra (5→100 µW) were already one column per power in `Book1`; I despiked them (130 cosmic-ray spikes removed) and ran a **Lorentzian** fit on every spectrum within the 1.23–1.45 eV window. All 15 fits converged to the **HIGH** tier (R² = 0.887–0.984, median 0.981).

**Per-power extracted values (from the Lorentzian fits):**

| Power (µW) | Intensity (counts) | xc (eV) |
|---|---|---|
| 5 | 98.8 | 1.3206 |
| 10 | 189.9 | 1.3244 |
| 15 | 267.9 | 1.3264 |
| 20 | 353.8 | 1.3278 |
| 25 | 425.2 | 1.3290 |
| 30 | 516.0 | 1.3300 |
| 35 | 564.8 | 1.3304 |
| 40 | 694.6 | 1.3315 |
| 45 | 738.5 | 1.3318 |
| 50 | 818.3 | 1.3328 |
| 60 | 935.9 | 1.3332 |
| 70 | 1140.4 | 1.3343 |
| 80 | 1264.7 | 1.3345 |
| 90 | 1538.0 | 1.3351 |
| 100 | 1714.6 | 1.3356 |

**Allometric fits (y = a·xᵇ):**
- **Peak intensity vs power:** I = a·Pᵇ with **b = 1.0048 ± 0.0261**, a = 16.28 ± 1.81, R² = 0.995 → essentially **linear (excitonic)** scaling, no saturation across 5–100 µW.
- **Peak position vs power:** a smooth, monotonic **blueshift** from 1.3206 eV (5 µW) to 1.3356 eV (100 µW), ≈ **5.0 meV per decade** of power (FWHM narrows 62→54 meV). This is a gentle sub-linear shift rather than a strong power law, so its effective exponent is much smaller than the intensity one.

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

The headline allometric exponent (intensity scaling) is:

RESULT: b = 1.005 ± 0.026