The dominant peak in each spectrum was fit with a Lorentzian plus a linear baseline over 1.24–1.42 eV. The fitted intensity is the Lorentzian amplitude above baseline. The high-energy shoulder blends into the rising edge, so its center could not be isolated reliably across all powers.

| Power (µW) | Fitted intensity (a.u.) | xc (eV) |
|---:|---:|---:|
| 5 | 123.2 | 1.31708 |
| 10 | 225.6 | 1.32224 |
| 15 | 319.8 | 1.32441 |
| 20 | 409.0 | 1.32612 |
| 25 | 500.6 | 1.32752 |
| 30 | 597.6 | 1.32853 |
| 35 | 673.2 | 1.32910 |
| 40 | 798.7 | 1.33043 |
| 45 | 858.8 | 1.33067 |
| 50 | 956.3 | 1.33176 |
| 60 | 1147.2 | 1.33226 |
| 70 | 1343.4 | 1.33357 |
| 80 | 1529.9 | 1.33376 |
| 90 | 1757.6 | 1.33445 |
| 100 | 2031.4 | 1.33508 |

The allometric fits use \(y=aP^b\), with log–log linear regression. The intensity exponent is \(b=0.9268\pm0.0120\) (1-SE, \(R^2=0.9978\)); the peak-position exponent is \(b=0.004426\pm0.000085\) (\(R^2=0.9953\)).

- [Peak intensity vs power plot](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-17/peak_intensity_vs_power.png>)
- [Peak position vs power plot](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-17/peak_position_vs_power.png>)
- [Lorentzian fit values CSV](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-17/lorentzian_peak_fits.csv>)

RESULT: b = 0.9268 ± 0.0120