I fitted the dominant peak in each spectrum with a Lorentzian and linear baseline. The table reports the fitted Lorentzian intensity above baseline and peak position \(x_c\).

| Power (µW) | Intensity (a.u.) | \(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 |

The allometric model was \(y=aP^b\), fit by linear regression in log space. The intensity fit gives \(b=0.9233\pm0.0121\) (\(R^2_{\log}=0.9978\)); the peak-position fit gives \(b=0.004466\pm0.000082\) (\(R^2_{\log}=0.9956\)). Uncertainties are standard errors of the fitted exponents.

- [Peak intensity vs power plot](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-14/peak_intensity_vs_power.png>)
- [Peak position vs power plot](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-14/peak_position_vs_power.png>)
- [Fitted peak values CSV](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-14/peak_fits.csv>)

RESULT: b = 0.9233 ± 0.0121