I fitted the dominant peak in each spectrum with a Lorentzian plus a linear baseline, using the 1.25–1.42 eV region. Peak intensity is the fitted Lorentzian amplitude; peak position is its center, \(x_c\). The allometric fits use \(y=aP^b\), with uncertainty in \(b\) from the log-linear fit.

- [Peak intensity vs power](peak_intensity_vs_power.png)
- [Peak position vs power](peak_position_vs_power.png)
- [Fitted peak parameters for each power](lorentzian_peak_fits.csv)
- [Analysis script](analyze_spectra.py)

RESULT: b = 0.9233 ± 0.0121