Tuning Graphene Oligomer Electromagnetic Response via Local Chemical Potential Control

Abstract

This study demonstrates that the electromagnetic characteristics of a graphene oligomer can be dramatically reshaped by locally tuning the chemical potentials of its constituent nanodisks. Variations in the chemical potential at different positions produce distinct impacts on extinction spectra and field distributions, enabling flexible localization of electromagnetic hotspots. The engineered nanostructures open new pathways for graphene‑based plasmonic devices in nanosensing, light trapping, and photodetection.

Introduction

Subwavelength metamaterials (MMs) have garnered significant attention for their ability to manipulate electromagnetic (EM) behavior, yielding phenomena such as negative refractive index, extraordinary optical transmission, and electromagnetically induced transparency. Plasmonic MMs, which exploit surface plasmons (SPs) in metal‑dielectric interfaces, enable strong field confinement and are essential for applications in sensing, spectroscopy, and nonlinear optics. However, conventional noble‑metal plasmonic structures suffer from high ohmic losses and limited tunability once fabricated.

Graphene, a two‑dimensional carbon lattice, supports surface plasmon polaritons (SPPs) across terahertz to mid‑infrared frequencies with lower losses and greater confinement. When arranged in oligomeric configurations—plasmonic molecules (PMs)—graphene nanodisks interact similarly to atoms in a molecule, producing collective resonances. Unlike prior work that uniformly adjusts the chemical potential across a structure, we investigate the effects of selectively varying the chemical potential of individual nanodisks within a graphene oligomer, providing unprecedented control over EM responses.

Simulation Methods and Model

Graphene is modeled as an atom‑thick film (Δ = 0.334 nm) with complex permittivity

ε = 1 + iσ_gη_0/(k_0Δ),

where σ_g is the surface conductivity obtained from Kubo’s formulation: σ_g = σ_intra + σ_inter. The intraband and interband conductivities are given by

σ_intra = (2e^2k_BT/πℏ^2)·i/(ω + iτ^{-1})·ln[2cosh(μ_c/k_BT)],

σ_inter = (e^2/4ℏ)[sinh(ℏω/2k_BT)/(cosh(μ_c/k_BT)+cosh(ℏω/2k_BT)) - (i/2π)ln((ℏω+2μ_c)^2/((ℏω-2μ_c)^2+4k_BT^2))].

Parameters: T = 300 K, τ = 0.5 ps, μ_c variable. The effective refractive index is

n_eff = 2iε_effε_0c/σ_g,

which depends strongly on μ_c.

The graphene oligomer consists of 13 identical nanodisks arranged with D_12h symmetry: a central disk surrounded by 12 disks on a dodecagon. Disk radius R_1 = 50 nm, outer radius R_0 = 240 nm. The oligomer rests on a silica substrate (n_2 = 1.5) and is surrounded by air (n_1 = 1). Light is incident normal to the plane, polarized along y. Finite‑element simulations are performed using COMSOL Multiphysics (RF Module), with perfect matched layers (PMLs) to absorb outgoing waves. The extinction cross‑section σ_ext = σ_sc + σ_abs is computed from surface and volume integrals of the Poynting vector and power‑loss density.

Results and Discussion

Local Chemical‑Potential Tuning of Plasmonic Modes

Figure 1 shows extinction spectra for μ_c ranging from 0.4 to 0.6 eV. Two resonances (A_0 and B_0) emerge, corresponding to Y‑mode (vertical) and X‑mode (horizontal) field distributions. Increasing μ_c of the entire oligomer enhances both peaks and shifts them to higher frequencies, owing to increased carrier density and reduced Landau damping.

To explore spatial selectivity, we partition peripheral disks into three groups (μ_c1, μ_c2, μ_c3) as illustrated in Figure 2. Keeping μ_c1 (intersection part) and the central disk at 0.5 eV, we vary μ_c2 or μ_c3 to 0.6 eV. This selective tuning modifies the relative intensities and positions of A_0 and B_0: raising μ_c2 weakens Y‑mode while strengthening X‑mode, whereas raising μ_c3 has the opposite effect. These observations confirm that local chemical potential changes reconfigure the EM landscape according to the symmetry of each mode.

Enhancement via Intersection‑Disk Tuning

By increasing μ_c of the four intersection disks (μ_c2) while holding others fixed, both Y‑ and X‑mode resonances are amplified (Figure 3). A new peak emerges adjacent to the Y‑mode as μ_c2 exceeds 0.6 eV, indicating mode splitting and degeneration. The two emergent peaks exhibit orthogonal Ey components, confirming they arise from the same parent mode. This demonstrates that adjusting the intersection disks can simultaneously boost multiple resonances.

Impact of the Central Disk

The central disk, distant from the peripheral ones, does not participate in the native modes. When its μ_c is increased to 0.8 eV, a dominant new resonance appears, suppressing the original peaks (Figure 4). This central mode dominates the extinction spectrum, illustrating that local tuning of non‑coupled disks can introduce entirely new resonances.

Fabrication Outlook

High‑quality monolayer graphene is typically grown via chemical vapor deposition (CH_4 source) and verified by Raman spectroscopy. Patterning employs electron‑beam lithography with PMMA resist, followed by O_2 plasma etching and acetone lift‑off. Chemical doping (e.g., HNO_3 vapor) or electrostatic gating (top‑gate voltage) can selectively alter μ_c of targeted disks, enabling the dynamic control demonstrated in simulations.

Conclusions

Selective tuning of local chemical potentials in a graphene oligomer affords unprecedented flexibility in shaping EM responses. By targeting specific disks—Y‑mode, X‑mode, intersection, or central—we can enhance, suppress, or create resonances without altering geometry. This capability opens new avenues for designing graphene‑based plasmonic devices in nanosensing, light trapping, and photodetection.