I fit the dominant peak in each spectrum over 1.245–1.425 eV with a Lorentzian plus a linear background. The CSV contains the fitted peak intensity, position \(x_c\), their fit uncertainties, and FWHM for all 15 powers. The plots show the allometric fits on log–log axes.

- [Peak intensity vs power plot](peak_intensity_vs_power.png)
- [Peak position vs power plot](peak_position_vs_power.png)
- [Fitted peak values](lorentzian_peak_fits.csv)

The intensity fit is \(I = 26.56\,P^{0.9231}\); the position fit is \(x_c = 1.3082\,P^{0.004446}\), with \(P\) in µW. The uncertainty on each exponent is the standard error from the log–log regression.

RESULT: b = 0.9231 ± 0.0121