I fitted the dominant peak in each spectrum over 1.25–1.42 eV with a Lorentzian plus a linear background. Intensity is the fitted Lorentzian height. The allometric fits used log–log linear regression; the reported uncertainties are standard errors on the exponent.

| Power (µW) | Peak intensity | Peak position \(x_c\) (eV) |
|---:|---:|---:|
| 5 | 125.1 | 1.316842 |
| 10 | 230.0 | 1.321914 |
| 15 | 325.5 | 1.324114 |
| 20 | 413.9 | 1.325922 |
| 25 | 505.7 | 1.327342 |
| 30 | 605.1 | 1.328307 |
| 35 | 680.4 | 1.328912 |
| 40 | 805.9 | 1.330269 |
| 45 | 867.3 | 1.330492 |
| 50 | 965.6 | 1.331583 |
| 60 | 1158.1 | 1.332092 |
| 70 | 1354.8 | 1.333414 |
| 80 | 1542.8 | 1.333605 |
| 90 | 1772.0 | 1.334302 |
| 100 | 2046.8 | 1.334944 |

- [Peak intensity vs power plot](peak_intensity_vs_power.png): \(b=0.9233\pm0.0121\), \(R^2_{\log}=0.9978\).
- [Peak position vs power plot](peak_position_vs_power.png): \(b=0.004466\pm0.000082\), \(R^2_{\log}=0.9956\).
- [Fitted values and parameter uncertainties](lorentzian_peak_fits.csv)

RESULT: b = 0.9233 ± 0.0121