The excitation energies of the randomly selected single configurations of the PNP and NR systems are shown in Figure and 15, respectively. It is immediately apparent that model B outperforms model A. The maximum errors from model A are 0.03, 0.20, and -0.16 eV for PNP1, PNP2, and PNP3, respectively, and 0.05 and 0.11 eV for NR-in-BLG and NR-in-POPC, respectively. With model B, these are lowered to 0.01, 0.05, and -0.05 eV for PNP1, PNP2, and PNP3, respectively, and 0.03 eV for both NR-in-BLG and NRin-POPC. This observation agrees with the analysis of the electrostatic potentials, where we found that the permanent potentials from model A differ substantially from model B (and C) and converge rather slowly. Also, in line with the previous analysis, we see that the cutoff distances at which the excitation energies from model B converge are roughly the