The errors in the excitation energies calculated using model C (and models A and B with long cutoff distances) along the MD trajectory are shown in Figures and. The errors range between -0.03-0.01 eV for PNP1, -0.04-0.04 eV for PNP2, -0.08-0.06 eV for PNP3, -0.04-0.02 eV for NR-in-BLG, and -0.02-0.02 eV for NR-in-POPC. As for the induced potentials, the errors in the excitation energies fluctuate around zero and therefore cancel out for the average excitation energies; the average error for PNP1, PNP2, and NR-in-BLG is zero, and roughly -0.01 and 0.01 eV for PNP3 and NR-in-BLG, respectively.