Particle-in-cell simulation of surface plasmon polaritons excited by external introduction of electron density

Abstract

Plasmonic devices work at visible and near-infrared frequencies, where a large number of bound electrons (polarization) in metals are excited. However, the existing particle-in-cell (PIC) software cannot take into account these bound electrons. In this paper, a new PIC simulation method is developed to study plasmonic devices, which has the ability to model both free electrons and bound electrons in metals. In this study, surface plasmon polaritons (SPPs) are excited by depositing excess electrons at one end of the metal film of thickness 140 nm (surrounded by air), thus initiating nonequilibrium in the electron density to start the oscillations. The spectrum, mode size, and propagation length of SPPs are calculated to confirm the electrical excitation of SPPs and characterize their properties. The excited SPPs are of a broad range of frequencies. Unlike SPPs excited by inelastic electron tunneling, the spectrum of SPPs excited in this way is not determined by eVbias but only by the properties of metals and dielectrics, and therefore the excited SPPs are intrinsic. The mode sizes of the intrinsic SPPs in the Ag–air geometry are 30 nm in Ag and 120 nm in air. The propagation length is ∼730 nm. The excited SPPs are converted into radiation modes by a sub-wavelength slit etched on the back side of the film to enable the detection of the radiation modes.