I fitted each spectrum with a Lorentzian plus a sloped local background over 1.27–1.40 eV. The peak intensity is the fitted Lorentzian amplitude. I then fit each set of peak values to \(y=aP^b\) by linear regression in log space; uncertainties are 1σ standard errors on \(b\).

The allometric exponents are \(b=0.9126\pm0.0140\) for peak intensity and \(b=0.00457\pm0.00010\) for peak position.

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

RESULT: b = 0.9126 ± 0.0140