Inter-Ribbon Plasmon Hybridization

• Inter-ribbon plasmon hybridization is the Coulombic coupling and mode mixing of charge oscillations in neighboring nanoribbons, producing distinct bonding and antibonding modes.
• Analytical and numerical models, including TDDFT and RPA, reveal that geometric and compositional parameters control exponential spectral splitting and anisotropic dispersion.
• Tunable coupling in superlattice arrays enables engineered mid-IR resonances, advancing applications in infrared photonics, metasurfaces, and quantum plasmonics.

Inter-ribbon plasmon hybridization refers to the Coulombic coupling and mode mixing of collective oscillations of charge carriers (plasmons) localized in or around neighboring nanoribbon structures. This phenomenon is essential in extended arrangements such as superlattice arrays, where the spectral characteristics and field profiles of plasmons are strongly impacted by geometric, material, and compositional parameters. It differs from simple particle–particle hybridization by the richness of its spectral features and the interplay of multiple coupling geometries, especially in systems comprising graphene, transition metal dichalcogenide, and anisotropic 2D ribbons. Plasmon hybridization in ribbons is central to optical field enhancement and nanoscale spectral engineering, with direct relevance to infrared photonics, metasurface design, and quantum plasmonics.

1. Physical Mechanisms of Inter-Ribbon Plasmon Coupling

Hybridization arises fundamentally from Coulomb interactions between the oscillating charge densities of adjacent ribbons. Two primary configurations are established:

Transverse (side-by-side) coupling: Parallel ribbons spaced at center-to-center distance $d$ exhibit near-field overlap of dipole oscillations ($p_1$, $p_2$) perpendicular to the ribbon axes. This produces an energy red-shift for the in-phase (“bonding”) mode and a blue-shift for the out-of-phase (“antibonding”) mode. The coupling strength $g$ displays exponential or Bessel function decay, $g \sim e^{-q d}$ or $g \sim K_0(q d)$, where $q \approx \pi/W$ is the characteristic plasmon wavevector for ribbon width $W$ (Rodrigo et al., 2015).

Longitudinal (end-to-end) coupling: Ribbons aligned tip-to-tip produce collective modes from dipoles oriented along the ribbon axes. Field leakage along the chain leads to more substantial splitting, with longitudinal coupling roughly $2\times$ stronger than transverse for equal spacing.

Rectangular nanoparticles (“ribbons”) present more complex hybridization. Each rectangle supports two orthogonal dipolar plasmon branches (horizontal/vertical), with distinct Lorentzian polarizabilities ($\alpha_H$, $\alpha_V$). The coupled-mode splitting is mediated by dyadic Green’s function elements, leading to anisotropic plasmon dispersion and, in rectangular symmetry, separation-dependent split modes that can diverge with increasing $d$ (Kuntman et al., 2020).

2. Analytical Models and Dispersion Relations

The coupled oscillators (bonding/antibonding) paradigm provides a general analytic framework. For two identical ribbons,

$$H = \sum_{i=1}^2 \hbar \omega_i^{(0)} a_i^\dagger a_i + \hbar g (a_1 a_2^\dagger + a_1^\dagger a_2)$$

yields hybrid plasmon frequencies: $$\omega_\pm = \frac{\omega_1^{(0)} + \omega_2^{(0)}}{2} \pm \sqrt{ \left( \frac{\omega_1^{(0)} - \omega_2^{(0)}}{2} \right)^2 + g^2 }$$ with $g$ derived from geometric and material parameters, $g \propto \omega_p K_0(q d)$.

For edge plasmons in 2D isotropic ribbons, the general bonding/antibonding dispersions are

$$\omega_{AB/B}(q;w) = \sqrt{ \frac{2}{\pi} \int_0^\infty du\, e^{-u}[K_0(u) \pm K_0(u + q w)] }\; \Omega_{bulk}(q)$$

(where $w$ is the ribbon width), with “+” for the antibonding (even) and “–” for bonding (odd) modes (Gonçalves et al., 2017).

Table: Splitting of Hybridized Plasmon Modes vs. Center-to-Center Separation

$d$ (nm)	$f_{bonding}$ (THz)	$f_{antibonding}$ (THz)
20	27 (–10%)	33 (+10%)
50	28.5 (–5%)	31.5 (+5%)
100	29.7 (–1%)	30.3 (+1%)

For graphene ribbons with $W = 100\,\text{nm}$, $E_F = 0.35\,\text{eV}$, this highlights exponential decay of mode splitting with $d$ (Rodrigo et al., 2015).

3. Superlattice Arrays and Extended Hybridization

Superlattice nanoribbon arrays, including systems of alternating doping or layer-number, exhibit plasmonic modes that reflect both local segment properties and multicellular band structure effects. In such arrays, the hybridization Hamiltonian generalizes to: $$H = \sum_j \hbar \omega_j^{(0)} a_j^\dagger a_j + \sum_{\langle i, j \rangle} \hbar g_{ij} (a_i a_j^\dagger + a_i^\dagger a_j)$$ with site energies set by width, doping, or layer-dependent $\omega_j^{(0)}$ and inter-segment couplings $g_{ij}$.

Resonance modes depend on the relationship between segment length $L$ and SPP wavelength $\lambda_{SPP}$:

Section $L$	Mode(s)	Frequency $f$ (THz)
$L \ll \lambda_{SPP}$	$m_0$ only	29 (intermediate)
$L \approx \lambda_{SPP}$	$m_0$, $m_1$	25, 35 (FP-type peaks)
$L \gg \lambda_{SPP}$	$m_0$, $m_1$	27 ($\mu_1$-like), 33 ($\mu_2$-like)

Short-period arrays yield a single intermediate-frequency mode not described by simple two-site hybridization. A plausible implication is that rapid modulation creates an effective chemical potential, supporting new hybrid modes between those of the isolated components (Rodrigo et al., 2015).

4. Numerical and Ab Initio Results

Time-dependent density functional theory (TDDFT) and random-phase approximation (RPA) computations provide spectral prediction and physical insight:

• Zigzag GNRs (semimetallic) exhibit 2D-like intraband plasmons with $\omega_p(q) \sim \sqrt{q n^*}$, tunable via room-T carrier populations or extrinsic doping.
• Armchair GNRs (semiconducting) display both intraband and interband plasmons. Hybridization, especially at intermediate doping ($\Delta E_F \sim 0.4\,\text{eV}$), produces a merged peak with splitting up to 0.2 eV at $q \sim 0.05\,\text{\AA}^{-1}$.
• The bonding/antibonding splitting and collective loss maps directly visualize hybridization, with mode mixing modeled by a dielectric matrix: $$\boldsymbol{\epsilon}(\omega,q) = \begin{pmatrix} \epsilon_{intra} & C(q) \ C(q) & \epsilon_{inter} \end{pmatrix}$$ yielding branches separated by $\Delta \omega \approx 2 C(q)$, where $C(q) \propto \exp(-q a)$. (Gomes et al., 2016).

5. Field Profiles, Spatial Character, and Anisotropy Effects

Hybrid modes in ribbons possess distinct charge distributions and field localizations, which are central to their optical response:

• For superlattices, field schematics reveal dipolar arrangements in high-doping segments and multipolar (node-rich) structures in low-doping sections, with parallel/antiparallel node parity influencing red and blue shifts and mini-bandgap formation (Rodrigo et al., 2015).
• Semi-analytic solutions in 2D ribbons show eigenmodes confined near the physical edges, with both bonding and antibonding modes possessing exponentially localized fields over $\sim 1/q$ scales. The splitting vanishes exponentially as $q w \gg 1$ (Gonçalves et al., 2017).
• Anisotropic materials (e.g., black phosphorus) modify hybridization via tensor conductivity, entering dispersion via $w_{eff} = w \sqrt{\sigma_{xx}/\sigma_{yy}}$ and yielding direction-dependent modal splitting.

6. Design Principles and Tunability

Control over inter-ribbon hybridization enables extensive spectral engineering in mid-IR and beyond:

• Single intermediate-frequency resonance: Achievable via short segments ($L \ll \lambda_{SPP}$) with alternated chemical potential or layers.
• Multiple Fabry–Perot resonances: Realized by matching segment length to SPP wavelength ($L \sim \lambda_{SPP}$).
• Isolated component resonances: Occur in the large-$L$ limit.
• Spectral tuning: Adjust via electrostatic gating ($\omega \propto \sqrt{E_F}$), ribbon width ($\omega \propto 1/\sqrt{W}$), inter-ribbon spacing (modifies splitting), and total layer count (augments chemical potential as $\sigma_{total} = \sum \sigma_i$).
• Band structure engineering: Introduce supercells with more than two segments or modulate doping/layer arrangements to synthesize minibands and gaps, following tight-binding analogies (Rodrigo et al., 2015).

In arrays, inter-ribbon coupling strength and modal mixing can be tuned smoothly, providing a flexible toolbox for designing metasurfaces with targeted plasmon resonances spanning the THz–visible domain (Gomes et al., 2016).

7. Unique Features and Practical Implications

Unlike simple particle–particle hybridization, inter-ribbon plasmon hybridization in extended systems introduces:

• Hybrid modes whose energies can lie strictly between those of isolated segments, contrary to predictions of the orbital hybridization model (Rodrigo et al., 2015).
• Unusual separation-dependent splitting in anisotropic or rectangular ribbons, where hybridized mode branches can diverge with increasing distance (Kuntman et al., 2020).
• Exponential sensitivity of coupling strength and spectral splitting to geometrical parameters such as ribbon width, inter-ribbon separation, and segment configuration.
• The possibility to create engineered, tunable mid-IR spectra and highly localized electromagnetic fields for photonic and optoelectronic devices, with spectral signatures readable via loss-function imaging and field-distribution mapping.

A plausible implication is that inter-ribbon plasmon hybridization will play a central role in future atomically thin photonic devices, offering both broad and sharply-defined spectral features, strong field localization, and tunable modal properties driven by nanoscale design.