I fitted the dominant peak in each spectrum over 1.27–1.40 eV with a Lorentzian and sloped baseline. Its intensity follows an allometric fit with exponent \(b=0.9985\pm0.0182\); the peak position fit gives \(b=0.00457\pm0.00010\).

- [Peak intensity vs power](intensity_vs_power.png)
- [Peak position vs power](peak_position_vs_power.png)

RESULT: b = 0.9985 ± 0.0182