“…While, within TDLDA, the Kohn-Sham wave functions are directly determined by the (self-consistent) confining potential, the local fields are, instead, defined solely by the induced charge density via the (screened) Coulomb potential and, therefore, depend on the confining potential indirectly. Hence, the deviation of E n from its classical behavior across the interface is determined by the electron density spillover over the classical (hard wall) boundary [73], and, therefore, is largely independent of the system overall shape. Furthermore, recent TDLDA studies of relatively large (up to 10 nm) nanoparticles revealed [55,56] that the main impact on plasmon linewidth comes precisely from the electron density tail and dielectric environment, implying that it is the plasmon local field near the interface, rather than electron wave functions, that chiefly determines the plasmon decay rate magnitude in real structures.…”