Electromagnetic scattering from a graphene disk: Helmholtz-Galerkin technique and surface plasmon resonances

The surface plasmon resonances of a monolayer graphene disk, excited by an impinging plane wave, are studied by means of an analytical-numerical technique based on the Helmholtz decomposition and the Galerkin method. An integral equation is obtained by imposing the impedance boundary condition on the disk surface, assuming the graphene surface conductivity provided by the Kubo formalism. The problem is equivalently formulated as a set of one-dimensional integral equations for the harmonics of the surface current density. The Helmholtz decomposition of each harmonic allows for scalar unknowns in the vector Hankel transform domain. A fast-converging Fredholm second-kind matrix operator equation is achieved by selecting the eigenfunctions of the most singular part of the integral operator, reconstructing the physical behavior of the unknowns, as expansion functions in a Galerkin scheme. The surface plasmon resonance frequencies are simply individuated by the peaks of the total scattering cross-section and the absorption cross-section, which are expressed in closed form. It is shown that the surface plasmon resonance frequencies can be tuned by operating on the chemical potential of the graphene and that, for orthogonal incidence, the corresponding near field behavior resembles a cylindrical standing wave with one variation along the disk azimuth.