“…The induced electron density ∆n DP (r) = n(r, t * ) − n(r, t = 0) is obtained from the time-dependent electron density n(r, t) by exposing the system to the linearly z-polarized electromagnetic wave with frequency Importantly, when Q electrons are added to the system, the change of the absorption spectra predicted by the classical theory, is not confirmed by the TDDFT calculations even qualitatively. Indeed, for the 2D-systems it is commonly assumed that the extra charge is homogeneously distributed over the nanoobject [26][27][28][47][48][49][50][51] , so that the plasmon frequency changes as ∆w p /w p = 0.5Q/N e . Thus, within this assumption, in the classical calculations the plasmon modes experience a blue-shift with increasing Q as shown in Fig.…”