Each spectrum was fit over 1.25–1.42 eV with a Lorentzian plus a linear background. Intensity is the fitted Lorentzian height. The allometric fits used log–log regression; uncertainties are 1 standard error.

| Power (µW) | Peak intensity (a.u.) | Peak position, xc (eV) |
|---:|---:|---:|
| 5 | 125.08 | 1.316842 |
| 10 | 229.97 | 1.321914 |
| 15 | 325.46 | 1.324114 |
| 20 | 413.89 | 1.325922 |
| 25 | 505.72 | 1.327342 |
| 30 | 605.11 | 1.328307 |
| 35 | 680.38 | 1.328912 |
| 40 | 805.90 | 1.330269 |
| 45 | 867.31 | 1.330492 |
| 50 | 965.56 | 1.331583 |
| 60 | 1158.08 | 1.332092 |
| 70 | 1354.78 | 1.333414 |
| 80 | 1542.83 | 1.333605 |
| 90 | 1771.96 | 1.334302 |
| 100 | 2046.81 | 1.334944 |

The fitted exponents are **b = 0.9233 ± 0.0121** for intensity and **b = 0.004466 ± 0.000082** for peak position.

- [Intensity vs power plot](intensity_vs_power.png)
- [Peak position vs power plot](peak_position_vs_power.png)
- [Fitted values CSV](spectral_fits.csv)

RESULT: b = 0.9233 ± 0.0121