Characteristic charge density pattern of Bennett multipole surface plasmon mode revealed by TDDFT simulation

The multipole branch of surface plasmon modes was theoretically predicted [1] and experimentally observed [2] long time ago, yet its charge and field patterns do not manifest easily in classical electromagnetism simulation [3]. We utilized the ab initio time dependent density functional theory (TDDFT) [4] and simulated the electromagnetic response of a jellium sphere. The multipole surface plasmon mode was successfully identified, and its dispersion relation and characteristic charge density pattern near the surface were clearly revealed.

References: [1] A. J. Bennett, Phys. Rev. B 1, 203 (1970); [2] K.-D. Tsuei, E. W. Plummer, A. Liebsch, K. Kempa, and P. Bakshi, Phys. Rev. Lett. 64, 44 (1989); [3] N. A. Mortensen, Nanophotonics 10(10): 2563-2616 (2021); [4] M.A.L. Marques, Alberto Castro, George F. Bertsch, and Angel Rubio, Comput. Phys. Commun. 151 60-78 (2003)