Microscopic Quantum Theory of SPPs

• The microscopic quantum theory of surface plasmon polaritons is a rigorous framework integrating quantum electrodynamics and many-body theory to describe plasmonic modes at metal-dielectric interfaces.
• It utilizes the Power–Zienau–Woolley Hamiltonian and derived quantum dispersion relations to capture geometry-dependent renormalizations and ultrastrong coupling effects.
• The theory underpins advances in quantum plasmonic scattering, molecule–plasmon interactions in nanocavities, and the design of nonclassical plasmonic devices for nanophotonics.

Surface plasmon polaritons (SPPs) are bosonic quasiparticles arising from the nonperturbative coupling between electronic collective oscillations (plasmons) in a metal and the electromagnetic field at metal–dielectric interfaces. Their fully quantum-mechanical description is essential for nanophotonics, quantum plasmonics, and the exploration of ultrastrong light–matter regimes. The microscopic quantum theory of SPPs integrates canonical quantum electrodynamics, many-body theory of electron gases, and boundary-conditioned electromagnetic quantization, providing rigorous predictive frameworks for the structure, dynamics, and fluctuations of SPP modes in arbitrary geometries.

1. Power–Zienau–Woolley Hamiltonian Formulation

The Power–Zienau–Woolley (PZW) transformation yields a decomposition of the total Hamiltonian for the plasmonic system, separating electronic and photonic degrees of freedom. For a metal of volume $V$ with free-electron density $\rho$ and bulk plasma frequency $\omega_p$,

$H = H_\text{el} + H_\text{ph} + H_\text{int}$

where

• Matter Hamiltonian (Bulk Plasmon):

$$H_\text{el} = \sum_{\mu} \hbar\omega_p \left(B^\dagger_\mu B_\mu + \tfrac12 \right)$$

with $B_\mu$ quantizing the collective polarization $\mathbf P(\mathbf r)$ of the electron gas.

• Photon Hamiltonian:

$$H_\text{ph} = \sum_{\mathbf{k},\sigma} \hbar c k \left(a^\dagger_{\mathbf{k},\sigma} a_{\mathbf{k},\sigma} + \tfrac12\right)$$

for free photon modes $a_{\mathbf{k},\sigma}$.

• Interaction Hamiltonian:

$$H_\text{int} = -i \sum_{\mu,\mathbf{k},\sigma} \hbar \left[C_{\mu,\mathbf{k},\sigma} B_\mu a_{\mathbf{k},\sigma} - \tilde{C}_{\mu,\mathbf{k},\sigma} B_\mu a^\dagger_{\mathbf{k},\sigma}\right] + \text{h.c.}$$

with geometry- and boundary-dependent coupling strengths $C_{\mu,\mathbf{k},\sigma}$ (Maurer et al., 16 Jan 2026).

The bulk plasmon is the unique matter oscillator. All interface and cavity effects enter through the electrostatic mode functions $g_\mu$ and the photonic coupling coefficients, which encode boundary conditions and geometry. The Hamiltonian remains valid for spheres, cylinders, planar films, and more complex stacks.

2. Geometry-Dependent Plasmon Renormalization

Plasmon–photon interaction induces a geometry-specific renormalization of the plasma frequency, revealing the microscopic origin of confined (localized and propagating) plasmonic resonances. For each mode $\mu$ determined by the geometry $\mathbf{g}$,

$\omega_\mu = \omega_p \sqrt{1 + \lambda_\mu}$

where $\lambda_\mu$ arises directly from electrostatic boundary conditions. Alternatively,

$$\tilde{\omega}_p(\mathbf{g}) = \omega_p \sqrt{1 + f(\mathbf{g},\epsilon_d)}$$

with $f(\mathbf{g},\epsilon_d)$ reflecting geometry and dielectric constant (Maurer et al., 16 Jan 2026). For a planar interface, one recovers the familiar surface plasmon resonance,

$$\omega_{\rm sp} = \frac{\omega_p}{\sqrt{\epsilon_\infty + \epsilon_d}}$$

These expressions underpin the quantum derivation of surface-plasmon dispersions in generic structures.

3. Quantum Dispersion Relations and Mode Structure

The exact quantum dispersion relation for propagating SPPs at flat metal–dielectric interfaces is extracted from the Green's-function pole condition or diagonalization of the coupled Hamiltonian:

$$k_\parallel(\Omega) = \frac{\Omega}{c} \sqrt{\frac{\epsilon_m(\Omega)\,\epsilon_d}{\epsilon_m(\Omega) + \epsilon_d}}$$

with $\epsilon_m(\Omega)$ the (possibly dispersive) metallic permittivity. The quantization of field and matter, and the explicit solution of Laplace's and Helmholtz equations via mode functions, yield the exact polaritonic eigenfrequencies including radiative and nonlocal effects (Maurer et al., 16 Jan 2026, Ballester et al., 2010). For planar and multilayer geometries, microscopic jellium models further resolve spatial profiles and eigenmode parities, allowing for explicit calculation of symmetric/antisymmetric branches and their hybridization (Cloots et al., 2024).

Table: Plasmonic Eigenmodes in Representative Geometries

Geometry	Quantum Mode Label $\mu$	Renormalized $\omega_\mu$
Planar interface	$\mu(k_\parallel)$	$\omega_p / \sqrt{\epsilon_\infty + \epsilon_d}$
Sphere ($l$-pole)	$\mu = l$	$$\omega_p \sqrt{ \frac{l}{ l(\epsilon_\infty + \epsilon_d ) + \epsilon_d} }$$
Multilayer slabs	$(S/A)$, hybrid	Mode splitting by $2\Lambda(q)$

Mode-specific renormalizations and splittings (e.g., for two different slabs: $\Delta\omega \sim 2\Lambda(q)$) are direct quantum effects (Cloots et al., 2024).

4. Ultrastrong Coupling and Ground-State Fluctuations

At metal–dielectric boundaries, the plasmon–photon coupling parameter $g_\mu$ is intrinsically nonperturbative:

$$g_\mu^2 = \sum_{\mathbf{k},\sigma} |C_{\mu,\mathbf{k},\sigma}|^2 / \omega_p$$

with $g_\mu / \omega_p \lesssim O(1)$, placing the system in the ultrastrong coupling (USC) regime. Counter-rotating terms are not negligible; the sum-rule renormalization from $\omega_p$ to $\omega_\mu$ is substantial.

Ground-state quantum fluctuations are a prominent manifestation. The mean bulk-plasmon population in the ground state, for mode $\mu$,

$n_\text{pl}^{(G)} = |x|^2$

with, for the dipolar mode of a sphere ($\epsilon_d = \epsilon_\infty = 1$),

$$n_\text{pl}^{(G)} = \frac{1}{\sqrt{3}} - \frac12 \approx 0.077$$

and for planar SPPs in the quasistatic limit,

$$n_\text{pl}^{(G)}(k_\parallel \gg \omega_p/c) = \frac{3}{4\sqrt{2}} - \frac12 \approx 0.03$$

These ground-state occupations are universal (typically $<0.1$) but can be tuned via $\epsilon_d$ or $\epsilon_\infty$ (Maurer et al., 16 Jan 2026).

Ground-state polarization and electric field fluctuations scale with $(|x|^2 + \frac12) \hbar\omega_p \epsilon_0 L_\mu$, dependent on dielectric refractive index, enabling active control over vacuum fluctuations at the interface.

5. Quantum Corrections in Jellium and Hydrodynamic Models

Microscopic jellium and hydrodynamic models extend the quantum theory to spatially resolved systems, capturing nonlocal and quantum statistical pressure effects. The hydrodynamic (Bohm–Staver) equation integrates with Maxwell’s equations, yielding quantum-modified SPP dispersions and field profiles (Bekshaev et al., 2021):

• Nonlocal corrections yield a blue shift in the SPP frequency at high $k_s$,
• Field confinement at boundaries tightens with increasing $k_s$,
• The energy and momentum densities split cleanly into “field” and “material” (electron motion) contributions, with orbital and spin decompositions.

Quantum hydrodynamic effects induce negative group velocity in certain modes and support rigorous analysis of backward wave propagation in thin-film SPP systems.

6. Plasmonic Nanostructures, Supercrystals, and Polariton Band Theory

In nanoparticle assemblies and supercrystals, the full quadratic minimal-coupling Hamiltonian includes intra-particle dipole and quadrupole operators, inter-particle Coulomb couplings, and photonic modes (Barros et al., 2021). Plasmon–photon interaction matrices are block-diagonalized via Bloch and Bogoliubov–Valatin transformations. The resulting plasmon–polariton bands exhibit:

• Rabi splittings and anticrossings,
• Dipole–quadrupole–photon mixing essential for accurate bandstructure,
• Reduced coupling strengths $\eta = \Omega_R / \omega_0$ entering the deep strong coupling (DSC) regime for metal fill fractions $f > 0.8$.

In DSC, the eigenmodes “freeze out” with distinct dipole, quadrupole, or photon character, indicating the regime of light–matter decoupling and breakdown of Purcell enhancement.

7. Quantum Plasmonic Scattering and Applications

Quantum theory for SPP scattering at interfaces employs a second-quantized description and transfer-matrix formalism (Ballester et al., 2010):

• SPP and radiation modes are coupled via interface coefficients,
• Transfer matrix relates input and output operators, enabling quantum beamsplitter designs,
• At balance ($\rho=\tau=1/2$), one realizes a 50:50 beamsplitter for SPPs,
• First-order quantum interference, analogous to the Hong–Ou–Mandel effect, directly reveals bosonic statistics for SPPs.

Microscopic quantum methods underpin applications in nonclassical state generation, quantum information processing, and surface-enhanced transport phenomena.

8. Quantum Strong-Coupling with Molecules and Prospective Directions

Plasmon–molecule coupling in nanocavities is addressed via quantum-electrodynamical coupled-cluster theory (QED-CCSD-1) (Fregoni et al., 2021):

• Localized surface plasmon modes are quantized, and coupled nonperturbatively to molecular states via transition charges,
• Full correlated polaritonic eigenstates recover, extend, and benchmark the Jaynes–Cummings predictions,
• Rabi splittings, ground-state correlation energies, and spatial charge redistribution are rigorously quantifiable,
• Nanocavity-induced changes in molecular charge density and oscillator strength can be manipulated via modal engineering.

This approach sets a standard for predictive, atomistic-level modeling of quantum polaritonics and supports extensions to multimode, ensemble, and dynamical analyses.

---

In sum, the microscopic quantum theory of surface plasmon polaritons provides a canonical and system-independent formalism for treating coupled electron–photon excitations in arbitrary metal–dielectric geometries. The PZW Hamiltonian isolates the bulk plasmon oscillator as the universal matter degree of freedom and nonperturbatively couples it to photonic continua, giving rise to geometry-dependent spectral renormalizations, exact dispersions, quantum fluctuations, and ultrastrong coupling phenomena (Maurer et al., 16 Jan 2026, Cloots et al., 2024, Bekshaev et al., 2021, Ballester et al., 2010, Fregoni et al., 2021, Barros et al., 2021). This framework enables the prediction and engineering of quantum plasmonic states and lays a robust foundation for nanophotonics, quantum optics, strong-coupling chemistry, and condensed-matter quantum technologies.