Plasmons in N-layer systems

Abstract

In multilayer structures, the coupling between layers gives rise to unique plasmon modes, but analytic solutions are typically available only for bilayers due to the increasing complexity as the number of layers increases. We investigate plasmons in multilayer structures, including the effects of interlayer tunneling. By introducing the Coulomb eigenvector basis for multilayer systems, which can be solved exactly using Kac-Murdock-Szeg\H{o} Toeplitz matrices, we analytically derive the long-wavelength plasmon dispersions both with and without interlayer tunneling. In the $N$-layer systems, we find that, in the absence of interlayer tunneling, the out-of-phase acoustic or charge neutral plasmon modes with linear dispersions ($\omega_\alpha\propto q/\sqrt{{1-\cos{\left(\frac{\alpha-1}{N}\pi\right)}}}$ for $\alpha = 2, 3, \cdots, N$) exist, while the in-phase classical plasmon mode exhibits its conventional dispersion ($\omega_1\propto \sqrt{q}$). When interlayer tunneling is present, the out-of-phase modes develop plasmon gaps that are governed by specific interband transitions, whereas the classical mode remains unaffected. These findings have broad applicability to general coupled-layer structures.