I found one resolved peak in each spectrum and fitted it over 1.24–1.42 eV with a Lorentzian plus linear background. The fitted intensity is the Lorentzian peak amplitude; \(x_c\) is the fitted center.

| Power (µW) | Intensity (a.u.) | \(x_c\) (eV) |
|---:|---:|---:|
| 5 | 123.2 | 1.317081 |
| 10 | 225.6 | 1.322238 |
| 15 | 319.8 | 1.324412 |
| 20 | 409.0 | 1.326125 |
| 25 | 500.6 | 1.327519 |
| 30 | 597.6 | 1.328525 |
| 35 | 673.2 | 1.329097 |
| 40 | 798.7 | 1.330428 |
| 45 | 858.8 | 1.330666 |
| 50 | 956.3 | 1.331756 |
| 60 | 1147.2 | 1.332260 |
| 70 | 1343.4 | 1.333570 |
| 80 | 1529.9 | 1.333758 |
| 90 | 1757.6 | 1.334451 |
| 100 | 2031.4 | 1.335083 |

The allometric fits use \(y=aP^b\), with uncertainty reported as the standard error of the fitted exponent.

- [Peak intensity vs power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-18/intensity_vs_power.png>): \(b=0.9268 \pm 0.0120\)
- [Peak position vs power](<<BENCH>/raw/commercial/luna/bench/t-pl/run-18/position_vs_power.png>): \(b=0.004426 \pm 0.000085\)

RESULT: b = 0.9268 ± 0.0120