I fit the dominant peak in each power spectrum with a Lorentzian plus a linear baseline over 1.25–1.40 eV. “Intensity” is the fitted Lorentzian peak height above the baseline. The allometric fits use \(y=aP^b\), fit by linear regression in log space.

| Power (µW) | Peak intensity | Peak position \(x_c\) (eV) |
|---:|---:|---:|
| 5 | 129.3 | 1.317343 |
| 10 | 235.3 | 1.322246 |
| 15 | 334.9 | 1.324550 |
| 20 | 426.0 | 1.326352 |
| 25 | 517.8 | 1.327701 |
| 30 | 621.3 | 1.328711 |
| 35 | 697.6 | 1.329293 |
| 40 | 825.2 | 1.330629 |
| 45 | 888.9 | 1.330868 |
| 50 | 991.0 | 1.331973 |
| 60 | 1190.7 | 1.332512 |
| 70 | 1399.4 | 1.333897 |
| 80 | 1588.9 | 1.334047 |
| 90 | 1821.6 | 1.334719 |
| 100 | 2100.4 | 1.335333 |

Plots: [peak intensity vs power](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-20/peak_intensity_vs_power.png>) · [peak position vs power](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-20/peak_position_vs_power.png>). The fitted data are in [lorentzian_peak_fits.csv](<<BENCH>/raw/commercial/luna-high/bench/t-pl/run-20/lorentzian_peak_fits.csv>).

The peak-position allometric exponent is \(0.00446 \pm 0.00008\) (1σ). For intensity, the exponent is \(0.923 \pm 0.012\) (1σ).

RESULT: b = 0.923 ± 0.012