I fitted a Lorentzian to the dominant peak in each spectrum over 1.26–1.42 eV, with a sloping baseline. Peak intensity is the fitted Lorentzian amplitude. Both plots use log axes and show allometric fits of the form \(y=aP^b\).

- [Peak intensity vs power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-23/peak_intensity_vs_power.png>)
- [Peak position vs power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-23/peak_position_vs_power.png>)
- [Lorentzian fit results for each power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-23/lorentzian_peak_fits.csv>)

The intensity allometric exponent is \(b=0.9188\), with standard error \(0.0124\) (\(R^2=0.9976\)). For peak position, the exponent is \(0.004512 \pm 0.000085\) (\(R^2=0.9954\)).

RESULT: b = 0.9188 ± 0.0124