[Paper Review] Plasmonic Hot Spots in Triangular Tapered Graphene Microcrystals

A. J. Rodin, Z Fei+5 more|arXiv (Cornell University)|Sep 7, 2013

Plasmonic and Surface Plasmon Research3 citations

TL;DR

This paper develops a simplified analytical model for plasmonic hot spots in triangular graphene microcrystals by replacing long-range Coulomb interactions with a short-range approximation, reducing the problem to a solvable Helmholtz equation. The model successfully reproduces experimentally observed unevenly spaced bright spots along the edges due to mode interference and reflection, demonstrating that geometric confinement alone can explain the key features of plasmonic hot spots in these systems.

ABSTRACT

Recently, plasmons in graphene have been observed experimentally using scattering scanning near-field optical microscopy. In this paper, we develop a simplified analytical approach to describe the behavior in triangular samples. Replacing Coulomb interaction by a short-range one reduces the problem to a Helmholtz equation, amenable to analytical treatment. We demonstrate that even with our simplifications, the system still exhibits the key features seen in the experiment.

Motivation & Objective

• To understand the origin of experimentally observed plasmonic hot spots in triangular graphene microcrystals, which appear as unevenly spaced bright spots along the edges.
• To address the complexity of electron-electron interactions and damping in graphene plasmonics by simplifying the Coulomb interaction to a short-range form.
• To demonstrate that geometric confinement in triangular shapes, rather than material inhomogeneities, is sufficient to explain the formation of plasmonic hot spots.
• To validate the simplified model against finite element simulations and experimental observations, showing qualitative agreement in signal amplitude and phase patterns.

Proposed method

• Replace long-range Coulomb interaction with a short-range interaction kernel to transform the problem into a Helmholtz-type equation amenable to analytical treatment.
• Use the hydrodynamic approximation and Ohm’s law to relate the induced charge density to the total potential via the conductivity σ and plasmon wave vector qp.
• Model the SNOM tip as a delta-function or point-dipole perturbation to compute the total potential response using Green’s functions and Fourier-Bessel expansions.
• Introduce a cutoff parameter γ to regularize the divergent response near the apex, simulating the finite angular size of the tip’s potential.
• Include a small imaginary part in qp (Im[qp]/Re[qp] ≈ 0.1) to account for experimental damping, enabling realistic amplitude decay away from edges.
• Solve the resulting eigenmode expansion numerically up to n=8 and compare results with finite element method simulations for validation.

Experimental results

Research questions

• RQ1 What causes the formation of unevenly spaced bright spots along the edges of triangular graphene microcrystals in s-SNOM experiments?
• RQ2 Can a simplified analytical model based on short-range interactions reproduce the key features of plasmonic hot spots observed in experiments?
• RQ3 To what extent is the plasmonic response in triangular graphene microcrystals governed by geometric confinement rather than material inhomogeneities or disorder?
• RQ4 How does the inclusion of damping and finite tip size affect the predicted signal distribution and hot spot pattern?
• RQ5 Can the Helmholtz equation with position-dependent qp and edge reflection conditions explain the interference patterns seen in the experimental signal maps?

Key findings

• The simplified model, based on a short-range interaction approximation, successfully reproduces the experimentally observed unevenly spaced bright spots along the edges of triangular graphene microcrystals.
• The signal amplitude map shows alternating bright and dark lines parallel to the edges, consistent with interference patterns from multiple reflected plasmon modes.
• The inclusion of a small imaginary part in the plasmon wave vector (Im[qp]/Re[qp] ≈ 0.1) leads to realistic amplitude decay away from the edges, matching experimental observations.
• The use of a cutoff parameter γ to model the finite size of the SNOM tip prevents unphysical divergence at the apex and improves agreement with experimental signal strength.
• Both the delta-function and point-dipole perturbation models yield qualitatively similar signal patterns, confirming that the core physics of mode interference and geometric confinement is robust.
• The finite element method simulation confirms the analytical results, validating that the simplified model captures the essential physics of plasmonic hot spot formation in triangular graphene microcrystals.