Plasmonic nanocavity-enabled universal detection of layer-breathing vibrations in two-dimensional materials

Conventional Raman spectroscopy faces inherent limitations in detecting interlayer layer breathing (LB) vibrations with inherently weak electron-phonon coupling or Raman inactivity in two-dimensional materials, hindering insights into interfacial coupling and stacking dynamics. Here we demonstrate a universal plasmon-enhanced Raman spectroscopy strategy using gold or silver nanocavities to strongly enhance and detect LB modes in multilayer graphene, hBN, and their van der Waals heterostructures. Plasmonic nanocavities even modify the linear and circular polarization selection rules of the LB vibrations. By developing an electric-field-modulated interlayer bond polarizability model, we quantitatively explain the observed intensity profiles and reveal the synergistic roles of localized plasmonic field enhancement and interfacial polarizability modulation. This model successfully describes the behavior across different material systems and nanocavity geometries. This work not only overcomes traditional detection barriers but also provides a quantitative framework for probing interlayer interactions, offering a versatile platform for investigating hidden interfacial phonons and advancing the characterization of layered quantum materials.

💡 Research Summary

The authors address a long‑standing challenge in Raman spectroscopy: the detection of interlayer layer‑breathing (LB) vibrations in two‑dimensional (2D) materials that are either Raman inactive or have intrinsically weak electron‑phonon coupling (EPC). By coupling multilayer graphene, hexagonal boron nitride (hBN), and van der Waals heterostructures with gold (Au) or silver (Ag) nanocavities, they create a universal plasmon‑enhanced Raman spectroscopy (PERS) platform that boosts the Raman signal of LB modes by orders of magnitude.

Fabrication involves depositing an 8 nm Au (or Ag) film onto the 2D flakes, which self‑assembles into irregular nano‑islands separated by nanogaps. Dark‑field scattering shows a plasmon resonance around 633 nm. When the excitation laser matches this resonance, low‑frequency Raman spectra reveal strong LB peaks that are absent in pristine samples. The enhancement persists across a range of layer numbers (N = 1–30) but diminishes for very thick flakes, reflecting reduced interfacial coupling.

A conventional linear chain model (LCM) predicts N − 1 LB modes for an N‑layer crystal. However, the authors observe N LB modes for N = 1–3, indicating additional “interface modes” (IM) arising from the graphene/substrate and graphene/nanocavity interfaces. By extending the LCM to include interfacial spring constants (k_Gr/Sub and k_Au/Gr) alongside the nearest‑ and next‑nearest‑neighbor interlayer force constants, they fit the measured LB frequencies and extract k_Gr/Sub ≈ 0.2 k_1st Gr and k_Au/Gr ≈ 0.3 k_1st Gr. This demonstrates that the two interfaces significantly modify the vibrational spectrum.

Polarization‑dependent measurements further reveal that, unlike conventional LB modes which appear only in co‑polarized (VV) or same‑handed circular (σ⁺σ⁺) configurations, the plasmon‑enhanced LB modes in AuNCs/6LG show comparable intensity in all four configurations (VV, HV, σ⁺σ⁺, σ⁺σ⁻). This indicates that the nanocavities not only amplify the local electromagnetic field but also alter the Raman tensor by modulating the interlayer bond polarizability.

To quantitatively describe this effect, the authors develop an electric‑field‑modulated interlayer bond polarizability model (E‑IBPM). Finite‑difference time‑domain (FDTD) simulations provide the spatial distribution of the local field |E_loc| across each graphene layer; the field enhancement varies from 1.5 to 3× depending on the layer’s distance from the nanocavity. The model introduces a field‑dependent polarizability term α′_i,xx that multiplies the interlayer spring constant, yielding an expression for the Raman intensity of each LB or IM mode. Using Lorentzian line shapes with layer‑dependent linewidths, the E‑IBPM reproduces the experimental spectra for N = 5, 10, 15, and 20 with excellent agreement in both peak positions and relative intensities.

Thermal annealing experiments (600 °C for AuNCs, 400 °C for AgNCs) demonstrate that reshaping the nanocavities reduces the field enhancement and consequently weakens the LB peaks, confirming the critical role of nanocavity geometry in the model parameters.

Overall, the work establishes four key advances: (1) a universal PERS approach that makes previously invisible LB vibrations observable across diverse 2D systems; (2) a modified LCM that incorporates interfacial coupling to accurately predict LB frequencies; (3) the E‑IBPM that explains the breakdown of conventional Raman selection rules and quantifies the synergistic effect of field enhancement and bond‑polarizability modulation; and (4) experimental validation that nanocavity material (Au vs. Ag) and morphology directly control the Raman response. This integrated framework opens a pathway for probing hidden interfacial phonons, weakly coupled excitations, and other low‑energy quasiparticles in layered quantum materials, thereby expanding the analytical toolbox for nanomaterials research.