Agency for Science, Technology, and Research (A*STAR), Institute of High Performance Computing, 1 Fusionopolis Way, 16-16 Connexis, Singapore, 138632, Singapore
Lim Jeremy Zhen Jie
SUTD-MIT International Design Center and Science and Math Cluster, Singapore University of Technology and Design (SUTD), 8 Somapah Road, Singapore, 487372, Singapore
Do Hue Thi Bich
NUS Graduate School for Integrative Sciences and Engineering, National University of Singapore, 21 Lower Kent Ridge, Singapore, 119077, Singapore
Xiong Xiao
Agency for Science, Technology, and Research (A*STAR), Institute of High Performance Computing, 1 Fusionopolis Way, 16-16 Connexis, Singapore, 138632, Singapore
Mahfoud Zackaria
Institute of Materials Research and Engineering, Agency for Science, Technology and Research (A*STAR), 2 Fusionopolis Way, Singapore, 138634, Singapore
Png Ching Eng
Agency for Science, Technology, and Research (A*STAR), Institute of High Performance Computing, 1 Fusionopolis Way, 16-16 Connexis, Singapore, 138632, Singapore
Bosman Michel
Department of Materials Science and Engineering, National University of Singapore, 9 Engineering Drive 1, Singapore, 117575, Singapore
Ang Lay Kee
SUTD-MIT International Design Center and Science and Math Cluster, Singapore University of Technology and Design (SUTD), 8 Somapah Road, Singapore, 487372, Singapore
Wu Lin
Agency for Science, Technology, and Research (A*STAR), Institute of High Performance Computing, 1 Fusionopolis Way, 16-16 Connexis, Singapore, 138632, Singapore
Particle simulation has been widely used in studying plasmas. The technique follows the motion of a large assembly of charged particles in their self-consistent electric and magnetic fields. Plasmons, collective oscillations of the free electrons in conducting media such as metals, are connected to plasmas by very similar physics, in particular, the notion of collective charge oscillations. In many cases of interest, plasmons are theoretically characterized by solving the classical Maxwell’s equations, where the electromagnetic responses can be described by bulk permittivity. That approach pays more attention to fields rather than motion of electrons. In this work, however, we apply the particle simulation method to model the kinetics of plasmons, by updating both particle position and momentum (Newton–Lorentz equation) and electromagnetic fields (Ampere and Faraday laws) that are connected by current. Particle simulation of plasmons can offer insights and information that supplement those gained by traditional experimental and theoretical approaches. Specifically, we present two case studies to show its capabilities of modeling single-electron excitation of plasmons, tracing instantaneous movements of electrons to elucidate the physical dynamics of plasmons, and revealing electron spill-out effects of ultrasmall nanoparticles approaching the quantum limit. These preliminary demonstrations open the door to realistic particle simulations of plasmons.