I fit the resolved band near 1.33 eV in each power spectrum with a Lorentzian plus linear background. The fitted peak intensity is the Lorentzian amplitude. The allometric fit to intensity gives \(b = 0.92295 \pm 0.01237\) (slope standard error, log–log fit).

Plots: [peak intensity vs power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-30/intensity_vs_power.png>) · [peak position vs power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-30/position_vs_power.png>)

Fit results for every power: [lorentzian_peak_fits.csv](<<BENCH>/raw/commercial/luna/bench/t-pl/run-30/lorentzian_peak_fits.csv>)

RESULT: b = 0.92295 ± 0.01237