“…In this section, the Monte Carlo method for simulating plasmon excitations is summarized. The plasmon scattering depth ( s ), polar ( θ ), and azimuthal ( ϕ ) scattering angles are estimated using the following formulas (Mendis, 2019; Barthel et al, 2020):where R 1 , R 2 , and R 3 are computer-generated linear random variables within the range [0,1], λ p is the plasmon mean free path, and θ E and θ c are the characteristic and critical plasmon scattering angles, respectively (Egerton, 1996). θ E can be calculated from E p /(2 E o ), where E p is the plasmon energy [17 eV for silicon (Mendis, 2019)] and E o is the primary beam energy.…”