Boron nitride nanoresonators for phonon-enhanced molecular vibrational spectroscopy at the strong coupling limit

Enhanced light-matter interactions are the basis of surface-enhanced infrared absorption (SEIRA) spectroscopy, and conventionally rely on plasmonic materials and their capability to focus light to nanoscale spot sizes. Phonon polariton nanoresonators made of polar crystals could represent an interesting alternative, since they exhibit large quality factors, which go far beyond those of their plasmonic counterparts. The recent emergence of van der Waals crystals enables the fabrication of high-quality nanophotonic resonators based on phonon polaritons, as reported for the prototypical infrared-phononic material hexagonal boron nitride (h-BN). In this work we use, for the first time, phonon-polariton-resonant h-BN ribbons for SEIRA spectroscopy of small amounts of organic molecules in Fourier transform infrared spectroscopy. Strikingly, the interaction between phonon polaritons and molecular vibrations reaches experimentally the onset of the strong coupling regime, while numerical simulations predict that vibrational strong coupling can be fully achieved. Phonon polariton nanoresonators thus could become a viable platform for sensing, local control of chemical reactivity and infrared quantum cavity optics experiments.


INTRODUCTION
Infrared spectroscopy is a powerful tool for label-free and nondestructive characterization of materials via their specific vibrational fingerprints 1 . However, the extremely small infrared absorption crosssections are challenging the detection and characterization of thin layers and small amounts of molecules. In order to overcome this limitation, surface-enhanced infrared absorption (SEIRA) spectroscopy has been implemented [2][3][4][5] . SEIRA relies on enhancing the interaction of infrared light with molecules via the strongly confined near fields at the surface of plasmon-resonant metallic structures, such as gratings 6 , nanoparticles or antennas 4,[7][8][9] . The same approach has been realized also with metal oxides 10 , highly doped semiconductors 11,12 and graphene [13][14][15] .
SEIRA experiments might also allow for exploring strong lightmatter coupling. In this regime, the light and matter states hybridize to form two polaritonic states that coherently exchange energy faster than the decay rate of the original states 16 . Recent studies showed that strong coupling of light and molecular vibrations can be applied to change and control the molecules' chemical reactivity 17 . Yet, molecular vibrational strong coupling (VSC) in the mid-IR spectral range is challenging because of the weak infrared dipole moment of molecular vibrations. So far, it has been achieved essentially with infrared microcavities 18,19 , and only very recently with mid-IR surface plasmons coupled to a micrometer-thick molecular layer 20 . However, molecular VSC at the nanoscale using plasmonic nanostructures has not been realized yet, owing to the high losses of plasmonic materials in the mid-IR range.
Strong infrared field enhancement and field confinement can be also achieved with surface phonon polaritons (SPhPs) in polar crystals, such as SiC [21][22][23][24][25][26] and quartz 27 . However, SEIRA based on SPhPs has been barely investigated so far 28,29 . SPhPs originate from the coupling of electromagnetic radiation to crystal lattice vibrations (phonons) in the so-called reststrahlen bands, where the real part of the dielectric permittivity is negative. They exhibit field enhancements and quality factors (Q) far beyond their mid-IR plasmonic counterparts 21,30 . Particularly, the recent emergence of van der Waals (vdW) materials opens new possibilities for phonon-polariton-based nanophotonics 31 , owing to the easy exfoliation of high-quality thin layers and subsequent nanofabrication via standard lithography techniques, as recently demonstrated with the prototypical vdW material hexagonal boron nitride (h-BN) 32,33 . With h-BN nanocones 33 and nanorods 34 , extremely narrow resonances (Q up to 283) 33 have been already demonstrated experimentally. In addition, vdW materials exhibit uniaxial anisotropy due to their layered crystal structure. This leads to a natural hyperbolic optical response in their reststrahlen bands, characterized by hyperbolic phonon polaritons (HPhPs) confined in the volume rather than at the surface of the material, which propagate with very high momenta 35,36 . The behavior of HPhPs in h-BN nanostructures and the creation of electric field hotspots on their surface have been recently studied 33,34 . However, their application to sensing and strong coupling studies has not been reported so far.
Here, we demonstrate SEIRA using, for the first time, well-defined infrared-phononic resonators. In particular, we use h-BN ribbons (HPhP Fabry-Pérot resonators) to detect small amounts of organic molecules. Moreover, we show that strong coupling between the HPhPs and the molecular vibrations can be reached with these simple resonating structures.

MATERIALS AND METHODS
We studied both experimentally and theoretically the infrared response of h-BN resonators without and with a thin molecular layer of the organic semiconductor 4,4′-bis(N-carbazolyl)-1,1′-biphenyl (CBP) 37 . The details on simulations, fabrication, experimental techniques and modeling are described below.

Electromagnetic simulations
Full-wave numerical simulations using the finite-elements method in frequency domain (COMSOL) were performed to study the spectral response of h-BN and Au resonators on a CaF 2 substrate, bare and with CBP molecules on top. The dielectric permittivity of Au was taken from Olmon et al. 38 , CaF 2 was described by a constant ε CaF2 ¼ 1:882, while h-BN and CBP dielectric permittivities were modeled as described below.

Dielectric function of h-BN
The dielectric permittivity tensor of h-BN is modeled according to where a = ||, ⊥ indicates the component parallel or perpendicular to the anisotropy axis. We use the parameters ε ||,∞ = ε ⊥,∞ = 4.52, o TO Refs 32,39).

Dielectric function of CBP
The transmission spectrum of a 100-nm-thick CBP layer was used to extract the dielectric function of CBP. To that end, we used a fitting procedure that uses the standard formula for the transmission of thin films on a substrate 40 . As shown in the Supplementary Information, the CBP transmission exhibits three dips in the range 1400-1550 cm − 1 , each one associated with a molecular vibration. To fit the transmission data, we consequently modeled the CBP permittivity ε CBP by three Lorentzians, according to The complete set of parameters obtained from the fit is listed in the Supplementary Information. For the dielectric non-dispersive background ε ∞ , we assumed the value ε ∞ = 2.8. This value provided the best fitting of the data, and is consistent with ellipsometry measurements reported in the literature 37 . For the central frequency and the width of the C-H deformation bond, which is the CBP vibration coupled to the HPhP resonance in this work, we obtained ω CBP = 1450 cm − 1 and γ CBP = 8.3 cm − 1 , respectively.

Fabrication of h-BN ribbon arrays
Large and homogeneous h-BN flakes were isolated and deposited on a CaF 2 substrate. To that end, we first performed mechanical exfoliation of commercially available h-BN crystals (HQ graphene Co, N2A1) using blue Nitto tape (Nitto Europe NV, Genk, Belgium). We then performed a second exfoliation of the h-BN flakes from the tape onto a transparent polydimethylsiloxane stamp. Using optical inspection and atomic force microscope (AFM) characterization of the h-BN flakes on the stamp, we identified high-quality flakes with large areas and appropriate thickness. These flakes were transferred onto a CaF 2 substrate using the deterministic dry transfer technique 41 . h-BN nanoribbon arrays of different widths were fabricated from the h-BN flakes by high-resolution electron beam lithography using poly(methyl methacrylate) (PMMA) as a positive resist and subsequent chemical dry etching. First, PMMA resist was spin-coated over the substrate. Second, arrays of ribbons of different widths were patterned on top of the flakes (i.e., the inverse of the final h-BN ribbon array pattern was exposed to the electron beam). Third, the sample was developed in MIBK:IPA (3:1), leaving PMMA ribbons of the desired width and length on top of the flakes. Fourth, using these ribbons as a mask, the uncovered h-BN areas were chemically etched in a RIE Oxford Plasmalab 80 Plus reactive ion etcher (Oxford Instruments Plasma Technology, Bristol, UK) in an SF6/Ar 1:1 plasma mixture at 20 sccm flow, 100 mTorr pressure and 100 W power for 20 s. Finally, the PMMA mask was removed by immersing the sample overnight in acetone, rinsing it in IPA and drying it using a N 2 gun.

Fourier transform infrared microspectroscopy measurements
Transmission spectra of the bare and molecule-coated h-BN arrays were recorded with a Bruker Hyperion 2000 infrared microscope (Bruker Optics GmbH, Ettlingen, Germany) coupled to a Bruker Vertex 70 FTIR spectrometer (Bruker Optics GmbH, Ettlingen, Germany). The normal-incidence infrared radiation from a thermal source (globar) was linearly polarized via wire grid polarizer. The spectral resolution was 2 cm − 1 .

Nanoscale Fourier transform infrared spectroscopy measurements
We used nanoscale Fourier transform infrared (nano-FTIR) spectroscopy to characterize the near-field response of the h-BN ribbons and to identify the nature of the observed HPhP resonance. We used a commercial scattering-type scanning near-field optical microscope setup equipped with a nano-FTIR module (Neaspec GmbH, Martinsried, Germany), in which the oscillating (at a frequency Ω ≅ 270 KHz) metal-coated (Pt/Ir) AFM tip (Arrow-NCPt-50, Nanoworld, Nano-World AG, Neuchâtel, Switzerland) is illuminated by p-polarized mid-IR broadband radiation generated by a supercontinuum laser (Femtofiber pro IR and SCIR; Toptica, Gräfelfing, Germany; average power of about 1 mW; frequency range 1200-1700 cm − 1 ). The spectral resolution was set to 4 cm − 1 . The spatial step size of the spectral linescan was set to 20 nm per pixel. To suppress background scattering from the tip shaft and sample, the detector signal was demodulated at a frequency 3Ω.

Classical model for coupled harmonic oscillators
In order to analyze the transmission spectra of the CBP-coated h-BN ribbon array, we phenomenologically described the coupling of the molecular vibrations and the phonon polaritons via a classical model of coupled harmonic oscillators. The equations of motion for the two coupled harmonic oscillators are given by 42 x :: where x HPhP , ω HPhP and γ HPhP represent the displacement, frequency and damping of the HPhP oscillator, respectively. F HPhP represents the effective force that drives its motion and is proportional to the external electromagnetic field. The corresponding notation is valid for the CBP vibration. g indicates the coupling strength and o ¼ o HPhP þoCBP 2 . Each damping term, γ HPhP and γ CBP , includes both radiative and dissipative losses. Note that radiative losses may even be negligible, because of the deep subwavelength-scale size of both the HPhP resonators and the molecules (i.e., the scattering cross-sections are much smaller than the absorption cross-sections). The extinction E ð Þ of such system can be calculated according to Ep/F HPhP t ð Þx : 43). We fit this model to the transmission spectra In the fitting procedure, γ CBP = 8.3 cm − 1 was fixed according to the CBP dielectric function, while ω CBP was limited within a few wavenumbers from its initial value (ω CBP = 1450.0 cm − 1 ), to allow for an eventual Lamb shift of the molecular vibration 44,45 . Both ω HPhP and γ HPhP were considered as free parameters, since the HPhP resonance undergoes a significant frequency shift and a potential linewidth modification compared to the bare ribbon array, due to the change of the dielectric environment once the molecules are placed on top of the ribbons.
Note that Equation (1) in Results and discussion section can be obtained by solving for the eigenvalues of the described model for F HPhP = F CBP = 0, with the approximation ω − ω i ≪ ω i and therefore It is the standard equation obtained in cavity quantum electrodynamics for strong coupling studies 46 .

RESULTS AND DISCUSSION
We first demonstrate the potential of HPhP nanoresonators for nanophotonics and sensing applications by comparing the spectral response of linear Au and h-BN antennas. To that end, we performed full-wave numerical simulations of the extinction cross-section for a canonical antenna structure ( Figure 1a): a rod, placed on top of a CaF 2 substrate, is illuminated by a plane wave with electric field polarized parallel to its main axis. The height and width were fixed to 100 nm, while the length was set to 245 nm for the h-BN rod and to 2300 nm for the Au rod, in order to obtain the same longitudinal dipolar resonance frequency of about 1450 cm − 1 (energy 180 meV, wavelength 6.9 μm). Note that the h-BN antenna is much shorter because of the extremely small wavelengths of the HPhP modes in the h-BN rod compared to the surface plasmon polariton modes in the Au rod 32 . In Figure 1c and 1e, we show the extinction cross-sections σ ext for both the h-BN (red line) and the Au (blue line) antenna normalized to their respective geometrical cross-section σ geo . Comparison of the spectra clearly shows the dramatically larger quality factor of the h-BN antenna. It exceeds that of the Au antenna by almost two orders of magnitude (Q Au ≈ 4 and Q h-BN ≈ 230), unveiling the great potential for low-loss metamaterials 47 , enhanced light-matter interaction 21,48 , sensing 4 , thermal emission 23,25,49,50 and cavity quantum optics 51 .
To explore the sensing capability of the h-BN antenna, we calculated the extinction spectra of the antennas covered by a 5-nmthick layer of the organic semiconductor 4,4′-bis(N-carbazolyl)-1,1′biphenyl 37 . We choose CBP molecules because of their well-defined vibrational modes located within the h-BN reststrahlen band. Further, CBP can be homogeneously deposited onto the antennas by simple thermal evaporation, and with sub-nanometer control of the layer thickness. The dielectric function of CBP was obtained by infrared spectroscopy (see Materials and methods section) and can be described by ε CBP (ω) = ε ∞ +ε vib (ω), where ε ∞ = 2.8 is a non-dispersive background and ε vib (ω) is the contribution owing to molecular vibrations. The calculated extinction spectra of the CBP-covered antennas are shown in Figure 1d and 1f. For the Au antenna (blue solid curves), we observe a small dip on top of the antenna resonance at 1450 cm − 1 , which can be associated with the absorption of the C-H vibration of CBP molecules. This kind of signature has been described in previous SEIRA experiments with metal antennas in terms of a Fano-type interference between the molecular vibration and the resonant surface plasmon polariton oscillations in the antenna (weak light-matter interaction) 52,53 . In stark contrast, the spectrum of the CBP-covered h-BN antenna (red solid curves) is spectrally shifted by nearly its full-width at half-maximum (due to the change in the dielectric environment; see discussion below) and its lineshape is dramatically modified. To highlight the effect of the molecular vibration, we calculated the spectrum of the h-BN rod covered by a 5-nm-thick layer with dielectric function ε = ε ∞ , that is, neglecting the contribution of the molecular vibrations (red dashed curve). The comparison of the results with and without the molecular vibrations (indicated by gray shaded areas in Figure 1f) clearly shows that the coupling between the molecule vibrations and the HPhPs results in a strong suppression and broadening of the HPhP resonance peak, which is accompanied by a clear dip. This large coupling makes h-BN antennas promising candidates for enhanced molecular vibrational spectroscopy.
In the following, we experimentally demonstrate SEIRA with an array of h-BN ribbons (we fabricated ribbons rather than rods for simplicity of the array design and sample processing), which were designed to exhibit a transverse dipolar HPhP resonance in the frequency range of the molecular vibrations of CBP. The FTIR transmission spectrum of the array without molecules (sketch in Figure 2a) clearly shows the HPhP resonance at 1465 cm − 1 (Figure 2b). From a Lorentzian fit we obtain a quality factor of Q exp ≅70, which is remarkably higher than the ones achieved for plasmon-resonant ribbons of metals or CVD graphene 13 .
To experimentally verify that the dip at 1465 cm − 1 in the transmission spectrum of Figure 2b is associated to the fundamental transverse HPhP mode, we performed a nano-FTIR spectroscopic linescan 54 across one ribbon. In nano-FTIR spectroscopy (Figure 2c), a metalized AFM tip is illuminated with p-polarized broadband infrared laser radiation. The tip acts as an antenna and generates strongly concentrated near fields at its apex, which locally excite HPhPs 34 in the h-BN ribbon. By recording the tip-scattered light with a Fourier transform spectrometer, we obtain local infrared amplitude and phase spectra of the h-BN ribbon with nanoscale spatial resolution 55 . In Figure 2d (lower panel) we show the nano-FTIR amplitude spectra s 3 (y,ω) recorded across one ribbon (along the dashed blue line in Figure 2c). We observe a resonance around 1460 cm − 1 (closely matching the dip in the far-field spectrum in Figure 2b) 56,57 , which can be only excited when the tip is located close to the ribbon edges (manifested by the two vertically aligned bright spots). We can thus identify this resonance as a transverse first-order (fundamental) dipolar HPhP resonance of the h-BN ribbon (illustrated in the upper panel of Figure 2d).
To demonstrate SEIRA, we measured infrared transmission spectra (blue curve in Figure 3a) of a ribbon array covered with a 20-nm-thick CBP layer (thermally evaporated, see Materials and methods section). Compared to the spectrum of the bare ribbons (black curve in Figure 3a), we observe a shift of the HPhP resonance dip by almost its full-width at half-maximum, as well as a dramatic change of its shape. We clearly recognize a significant modulation at the spectral position of the C-H bond, which is combined with a resonance broadening (analogous to Figure 1f). We attribute this remarkable redistribution of spectral weight to the coupling of the CBP vibration to the HPhP resonance, due to the strongly localized fields in the proximity of the h-BN ribbons. For better visibility of the spectral modification, we also show (as guide to the eye) the spectrum of the bare ribbon array (dark grey curve), which was artificially frequencyshifted in order to mimic the effect of ε ∞ . We stress that the C-H vibration can be barely observed in the spectrum of the 20-nm-thick CBP layer on the bare CaF 2 substrate (lower brown curve), that is, without the presence of any HPhP field enhancement.
We also study the effect of the CBP layer thickness. To that end, we removed the 20-nm-thick CBP layer and evaporated thinner layers of 1, 3 and 10 nm on top the ribbons. As the CBP thickness increases, we observe an increasing red-shift of the HPhP resonance (Figure 3a). This effect, corroborated by full-wave electromagnetic simulations (Figure 3b), is caused by the change of the dielectric environment of the h-BN ribbons (from ε air to ε ∞ ). We stress that the narrow linewidth of the HPhP resonance allows for resolving even small shifts, which could be applied to determine the thickness of a thin molecular layer placed on top of the ribbons. In addition, we observe that the vibrational feature starts to affect the lineshape of the HPhP resonance dip already for the 3-nm-thick CBP layer (yellow curve), and becomes more clearly visible for a 10-nm-thick layer (green curve). We stress that with our initial phonon-enhanced SEIRA experiments we achieved already femtomolar sensitivity, although we performed the experiments Boron nitride nanoresonators for SEIRA spectroscopy M Autore et al with a thermal source, a rather small ribbon area (20 × 20 μm 2 ) and at ambient pressure. In the future, the sensitivity could be further improved by matching the HPhP resonance to the molecular vibration for each layer thickness, and by optimizing the field enhancement and confinement by further engineering the HPhP resonators' geometry. The striking modification of the HPhP resonance caused by the 20nm-thick CBP layer indicates a rather strong interaction between the HPhPs and the molecular vibrations, which could go beyond a Fano interference. In order to explore the possibility to reach strong coupling between the CBP vibrations and the HPhPs of the h-BN ribbons, we studied the interaction between the two modes as a function of their spectral detuning in order to extract the coupling strength 16,58 . To that end, we tuned the HPhP resonance across the molecular vibration by changing the h-BN ribbon width w from 85 to 162 nm, while fixing the ribbon period D = 400 nm. In Figure 4a we plot the relative transmission spectra of the bare ribbon arrays of different ribbon widths w (black lines). We see that the resonance (dip) in the transmission spectrum shifts to higher frequencies as w decreases, as expected for a Fabry-Pérot resonator. To verify the resonances, we performed numerical simulations of the normalized transmission spectra T/T 0 , which are plotted in Figure 4b as a function of the inverse ribbon width, w − 1 (proportional to the mode momentum). We observe a well-defined dip in the spectra (indicated by brown color in Figure 4b) shifting to higher frequencies as w − 1 increases, in good agreement with the experimental dip positions (green dots, extracted by Lorentzian fits of the experimental spectra).
From the calculated electric field profile in the inset of Figure 4b, we estimate that 85% of the electric field outside the ribbon is confined within the first 30 nm from the surface. We thus evaporate 30 nm of CBP onto the ribbon arrays, in order to achieve a large overlap between the HPhP and the molecular vibrational mode, which is key to maximize their coupling strength and eventually reach the strong coupling regime 59 . Relative FTIR transmission spectra of the CBPcovered h-BN ribbons (black curves in Figure 4c) show two dips. We extract the central frequency of each dip by a Lorentzian fit (see Supplementary Information) and plot it as a function of w − 1 in Figure 4d (green dots). A good agreement with simulations of the relative transmission spectra (T/T 0 , contour plot in Figure 4d) is found. Interestingly, we observe that the transmission minima (dips) do not cross, indicating that the HPhP and the CBP vibrational modes are coupled 42 . In order to determine whether we observe weak or strong coupling, a quantitative analysis of the spectral separation of the eigenmodes of the coupled system and comparison to the damping of the uncoupled modes is required 16 . To that end, we modeled the HPhP resonance and the CBP vibrational mode as two classical harmonic oscillators interacting with a coupling strength g, each characterized by a central frequency and a damping parameter (ω HPhP and γ HPhP for the HPhP resonance, and ω CBP and γ CBP for the CBP vibration, respectively; see Materials and methods section for details on the model). The fits to the transmission curves are shown as red dotted lines in Figure 4c. The extracted values of the coupling strength for each spectra are plotted in the inset of Figure 4e, yielding an average value of g = 7.0 cm − 1 . With the parameters extracted from the fits we are able to calculate the frequencies ω ± of the new eigenmodes of the coupled system, according to 60 where δ = ω HPhP − ω CBP is the detuning between the HPhP resonance and the CBP vibration. In Figure 4e we plot ω ± as a function of ω HPhP (green dots). We clearly observe anti-crossing, which is a necessary condition for strong coupling. It is equivalent to the mathematical condition C 1 ¼ def jgj=jg HPhP 2g CBP j40: 25 (Refs 16,58). With our experimental data we indeed find C 1 ≅0.37. However, the more restrictive criterion C 2 ¼ def jgj=jg HPhP þ g CBP j]0:25 is required to see emerging strong coupling effects 16 . Measuring C 2 ≅0.19 thus indicates that we reach the onset of strong coupling between the HPhP mode and the molecular vibrations. We note that in many strongly coupled systems the mode splitting (ω + − ω − ) can be increased with the number N of molecules inside the mode volume V, that is, by increasing the density of molecules 16 . In our experiment, however, layers with a thickness beyond 30 nm did not increase the splitting. This can be explained by additional molecules being deposited in a region where the electromagnetic field associated with the HPhPs vanishes (due to their evanescent nature), that is, the additional molecules do not increase the molecule density within the mode volume.
In order to corroborate that strong coupling between HPhPs and molecular vibrations can be achieved, we performed numerical simulations of the transmission spectra of the bare ( Figure 5, red curve) and CBP-covered (blue curve) h-BN ribbon array, using literature data for the h-BN permittivity (as described in Materials and methods section). From a Lorentzian fit to the calculated transmission of the bare ribbons we obtain Q calc ≅120. By analyzing the calculated spectra for several ribbon arrays with different values of w and a 30-nm-thick CBP layer on top (using the coupled oscillators , for a CBP-covered h-BN ribbon array (same color notation of a) and for a bare CBP layer (brown curves). The calculated spectrum for the bare ribbons is repeatedly shown in the background of the other spectra (gray curves, shifted along the frequency axis). The gray shaded areas visualize the difference between the gray and colored spectra.
model described above), we obtain C 1 ≅1.4 and C 2 ≅0.31 (see Supplementary Information for details). According to the definitions provided above, both C 1 and C 2 values indicate strong coupling between the CBP molecular vibrations and the HPhPs. We support this finding by calculating the absorption of the CBP molecules (green curve in Figure 5, obtained by integration over the whole layer volume). We clearly see two absorption maxima, which further supports the mixed-state character of a strongly coupled system. 59,61 Our numerical calculations predict that strongly coupling in principle could be achieved in our experiment. We explain the discrepancy to the experiment (i.e., that the values C 1 and C 2 are smaller in the experiment than in the numerical simulations) by the reduced experimental quality factor. From the Lorentzian fits of the experimental data in Figure 4a we obtain an average quality factor Q exp ≅50, which is more than a factor of 2 worse than in the numerical simulations. We explain this finding by fabrication-induced roughness and width variations along the h-BN ribbons, defects, and residual impurities on the ribbons, which lead to increased HPhP scattering and inhomogeneous resonance broadening. By improving the quality of the ribbons we envision the full achievement of VSC in the future.

CONCLUSIONS
In summary, we used h-BN ribbon arrays for SEIRA spectroscopy of molecular vibrations, exploiting the field enhancement due to HPhP resonances. We proved at least femtomolar sensitivity, showing that phonon-enhanced molecular vibrational spectroscopy could be an interesting alternative to conventional plasmon-based SEIRA for the identification and monitoring of small amounts of molecules. Further increase of the sensitivity could be achieved by placing smaller quantities of molecules directly on top of the hotspots of the infraredphononic nanoresonators, or in the gap of more sophisticated antenna and resonating structures. In addition, we found that in our experiment the interaction between the HPhPs and the molecular vibrations already reached the onset of strong coupling, while numerical simulations predict the full realization of VSC. We envision full VSC in future experiments by improving the Q-factor of the h-BN nanoresonators through optimization of the fabrication process. An even deeper strong coupling regime could be achieved with more sophisticated HPhP resonator geometries (such as gap antennas or split ring resonators) or by using isotopically enriched h-BN supporting ultralow-loss HPhPs 62 .
Our findings could open novel exciting perspectives towards local control of chemical reactivity 17,63 , site-selective catalysis and cavity quantum optics 51,64 applications.

CONFLICT OF INTEREST
R Hillenbrand is co-founder of Neaspec GmbH, a company producing scattering-type scanning near-field optical microscope systems, such as the one used in this study. The other authors declare no conflict of interests.