# Honeycomb-like aluminum antennas for surface-enhanced infrared absorption sensing

Melissa Najem, Franck Carcenac, Luka Coutaud, Mohamed Mouhibi, Thierry Taliercio, Fernando Gonzalez-Posada

<figure><img src="image_1.png" /></figure>

## To cite this version:

Melissa Najem, Franck Carcenac, Luka Coutaud, Mohamed Mouhibi, Thierry Taliercio, et al. Honeycomb-like aluminum antennas for surface-enhanced infrared absorption sensing. Nanophotonics, 2023, 12 (12), pp.2199-2212. 10.1515/nanoph-2023-0131. hal-05089054

HAL Id: hal-05089054

https://hal.science/hal-05089054v1

Submitted on 28 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Research Article

Melissa Najem*, Franck Carcenac, Luka Coutaud, Mohamed Mouhibi, Thierry Taliercio and Fernando Gonzalez-Posada

# Honeycomb-like aluminum antennas for surface-enhanced infrared absorption sensing

https://doi.org/10.1515/nanoph-2023-0131

Received February 27, 2023; accepted March 30, 2023;

published online April 10, 2023

**Abstract:** Surface-enhanced infrared absorption (SEIRA) spectroscopy is a competent method to detect trace quantity of molecules and even protein conformational flexibility by enhancing their vibrational modes. To improve the spectroscopy features, we propose a surface with honeycomb-like (HC) arrangement of aluminum equilateral triangles within a metal-insulator-metal configuration. With adjustable geometric parameters, the HC nanoantennas allow a tunable and wide spectral coverage in the IR. The reflectance measurements correlate extremely well with the numerical simulations. Being compact and insensitive to the light polarization, the HC are appealing for boosting the signal-to-noise ratio and increasing the number of hotspots as required for sensing applications. These nanoantennas are thus suitable for accurate and broadband SEIRA sensing via a spectral overlap between the large plasmonic resonances and the narrow IR vibrational modes of our analyte (vanillin). In line with our previously studied bowties nanoantennas, we demonstrate, using HC, SEIRA enhancement factors greater than 10⁷ achieved at a tuning ratio below 1 stating the best spectral overlap. Around 10⁴ molecules are sensed per HC tip. The investigation results are matching the best-reported SEIRA studies. These findings pave the way toward sensitive, adaptable, and miniaturized IR spectroscopy devices for vital applications like biosensing and environmental monitoring.

**Keywords:** aluminum; honeycomb nanoantennas; MIM; plasmonics; SEIRA spectroscopy; wide tunability.

## 1 Introduction

Customized plasmonic nanostructures are building blocks for various optical applications, among those sensing. Surface-enhanced infrared absorption (SEIRA) is an appealing technique for different domains, such as electrochemistry [1], bioelectrochemistry [2], and sensing [3]. The term “SEIRA” has been first announced by chemists in the early 90s [4]. Since then, it has been massively introduced as a molecule detection and identification method [5]. It capitalizes on the excitation of the localized surface plasmons as collective oscillations of the free-electron cloud of a metallic structure [6]. The molecular fingerprints are identified by resonant detection of their vibrational modes through the coupling to the localized surface plasmons resonance (LSPR) of the nanostructures [7]. Once the LSPR ($\omega_{res}$) is situated at the same frequency as a molecular vibration ($\omega_{vib}$), the nanostructure-molecule system can couple [8], resulting in spectral features with Fano-like lineshapes [5, 9], characteristics of a coupling between a broad and a narrow energy states [10]. The Fano-like profiles depend on a tuning ratio between the plasmonic and molecular mode frequencies. Following such coupling, their interaction cross-section increases [5]. Hence, SEIRA overcomes the limitation of conventional IR absorption techniques that cannot afford good sensitivity due to small IR absorption cross-sections such as for gases, or very thin samples as monolayers.

The enhanced IR absorption can only originate from the molecules located within the hotspots of plasmonic nanostructures [5, 11]. Sharp and closely-placed nanoantennas are thus desired to enhance the near electric field (E-field) in-plane of the antennas through both the lightning-rod and the gap effects [12]. The characteristics of metallic nanoantennas are widely studied in the literature, like cross antennas [13] and fan-shaped antennas by Brown et al. [14],

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<footnote>*Corresponding author: Melissa Najem, University of Montpellier, Institut d'Electronique et des Systèmes, Montpellier, Occitanie, France, E-mail: melissa.najem@umontpellier.fr. https://orcid.org/0000-0002-7125-517X</footnote>

<footnote>Franck Carcenac, CNRS Laboratory for Systems Analysis and Architecture, Toulouse, Occitanie, France</footnote>

<footnote>Luka Coutaud, Mohamed Mouhibi, Thierry Taliercio, and Fernando Gonzalez-Posada, University of Montpellier, Institut d'Electronique et des Systèmes, Montpellier, Occitanie, France</footnote>
nanogapped BT antennas by Dong et al. [15], and dendritic
resonators by Wallace et al. [16], to name a few. Further-
more, the size of the nanoantennas is a promising degree of
freedom to tune their LSPR frequency which goes along with
the near E-field localization and enhancement [17]. Based on
this, a barcode-like wide spectral coverage can be foreseen
across targeted frequency ranges of the electromagnetic
spectrum [18, 19]. In addition, a compact and symmetric
arrangement could be a smart engineering to ameliorate
the plasmonic nanoantennas behavior. Inspired by the most
efficient manner of filling space as studied in crystallogra-
phy [20], equilateral triangles will be arranged in a hexag-
onal packing that mimics the natural honeycombs. This
special organization will allow a symmetry invariance that
might be convenient to ease the incident light implementa-
tion as well as to improve the signal-to-noise ratio during the
measurements.

In this article, we investigate the plasmonic behav-
ior of equilateral triangles organized as honeycomb-like
(HC) arrays. These plasmonic nanoantennas are fabricated
with aluminum (Al). Lately, Al has demonstrated great
promises as alternative material to noble metals, being
affordable, abundant, and compatible with silicon tech-
nologies [21, 22]. Herein, the HC antennas are incorporated
within a metal-insulator-metal (MIM) configuration. The
MIM allows a strong lateral confinement of the E-field in
the insulator (dielectric) [23, 24]. The combination of both
the lateral and the in-plane E-field confinements leads to
maximize the E-field intensity in the gap region between
the neighboring triangles as required for SEIRA application.
LSPR positions of the HC nanoantennas are extended over
a broad IR spectral range by controlling the side length of
the triangles. This platform is thus capable of enhancing
multiple IR vibrational features offering a broadband SEIRA
sensing.

Recently, Aslan et al. have employed honeycomb-
shaped plasmonic nanoantennas as photonic metamaterials
within a MIM structure to demonstrate a potential SEIRA
sensor through the detection of a polymer nanolayer [25].
A large electromagnetic near-field enhancement was pro-
vided at the activated hotspots that are located at the sharp
extremities of a completely-filled HC building-block. How-
ever, in our present article, an advanced highly-compact and
highly-activated structure is promoted: more hotspots are
triggered by filling the HC unit-cell with six sharp equilateral
triangles. HC arrays with different sizes are studied exper-
imentally by single detector Fourier transform infrared
(FTIR) measurements and by finite-difference time-domain
(FDTD) simulations. Their corresponding LSPR are defined
in a wide IR range. The invariance of the HC structure to

the light polarization is proved. The advantage of the com-
pact HC arrangement is also underlined and supported by
a brief comparison to our formerly published bowties (BT)
results [18]. SEIRA results are as well compared for both
nanoantennas. As a case study, we demonstrate a broad-
band SEIRA using vanillin as a trial analyte. The vanillin
molecule exhibits several sharp modes characteristics of
the vibrational analysis of its aromatic ring. The analysis
of the resulting Fano-like vibrational signals of a µ-droplet
of vanillin/IPA solution shows a maximum enhancement
factor higher than 10⁷ reached at a frequency tuning ratio
below 1. The number of sensed molecules was analytically
found around 10⁴ molecules per HC tip. These results are
in agreement with the stated results in the literature as
well as in our earlier publication about Al-BT [18]. These
outcomes are encouraging for the conception of a com-
pact, cost-effective, and sensitive biosensor based on silicon
technologies.

**2 Methods**

2.1 Simulation

The electromagnetic response of the HC antenna integrated in
a MIM configuration has been simulated by three-dimensional
FDTD calculations using Ansys Lumerical 2021 R2.3 Finite difference
solver.

The MIM structure has been modeled as a 50 nm-thick aluminum
(Al) HC with variable side length *L* and gap *g*, on top of a 100 nm thick
Al layer (mirror), with a 20 nm thick dielectric (SiO₂) layer stacked in
between. For the Al, we referred to the dielectric constants reported
in CRC Handbook of Chemistry and Physics by E. Shiles [26]. The SiO₂
material data are based on the Handbook of Optical Constants of Solids
I – III by E. Palik [27], but are slightly modified to match the real quality
of the deposited material. An axiomatic refractive index equal to 1
(representing air) is maintained in the structure surroundings to create
an asymmetric environment.

HC-like pattern is obtained by creating three groups of periodic
lines rotated from the horizontal x-axis with an angle of 0°, 60°, and
120° degrees within the 50 nm thick Al layer. The lines are defined
as etched material. The lines width *w* and periodicity *p* are affecting
the sculpted Al triangles size (triangle's height equal to *L* × sin 60°)
and their spacing distance (that is half of the eventual tip-to-tip gap
*g*). Symmetric and anti-symmetric boundaries are employed along the
*x*- and *y*-direction, and perfect matching layers (PMLs) are along the
*z*-direction. A mesh grid of 5 nm is used around the intersection gap.
An incident plane-wave is co-or cross-polarized (E-filed component
along or across the tip-to-tip facing pair of triangles, respectively)
and launched backward along the *z*-axis for illumination. Thereafter,
optical responses of Al-HC are collected using a monitor placed above
the excitation source, so-called reflectance monitor. A field monitor is
introduced in the *xy*-plane to represent the electric field cartography.
Throughout both monitors, the influence of geometrical parameters on
the plasmonic behavior is studied.
## 2.2 Fabrication

The main substrate fabrication is detailed in our previous work [18]. In addition, onto the 20 nm-thick SiO₂ layer, a 50 nm-thick Al layer is thermally evaporated at 1 Å/s.

On top of this layer, the HC patterns are engineered by electron beam lithography (EBL) using dilute positive tone resist (CSAR) within the extent of 100 µm × 100 µm. A Schottky-emitter field emission gun (FEG) EBL system (RAITH150 equipment) is used with a 10 µm aperture, operating at 20 kV (beam current of 40 pA), a working distance of 7 mm. An area stepsize of 4 nm with an area nominal dose of 50 µC/cm²; and a line stepsize of 2 nm with a line nominal dose of 100 pC/cm² are selected. During the exposure, the electron beam is intentionally scanned over the sample to define a matrix of 5 × 11 zones, in each of which L and g are fixed. After electron exposure, the exposed resist is developed in an amyl acetate solution for 60 s, and then the sample is immediately rinsed in IPA (stopper solution) for 30 s.

To transfer the patterns into the metal layer, the 50 nm-thick Al layer is etched via the opened lines using argon (Ar) sputtering into an ICP-Corial machine (150 W as RF and 200 W as LF, with a voltage bias of 170 V) for 180 s. At the end, the unexposed resist was cleaned by an inductive plasma O₂ (100 % O₂) for 5 min.

## 2.3 Characterization

The etching depth is surveyed using atomic force microscopy (AFM) analysis to calibrate the etching rate. The morphology of the tailored HC antennas is examined by a Hitachi S4800 scanning electron microscope (SEM) at 10 keV. The optical response of the fabricated resonators is characterized through a Hyperion 3000 microscope coupled to a Fourier transform IR (FTIR) spectrometer (Bruker Vertex 70). A broadband IR light source, namely Globar (silicon carbide SiC-based-MIR source) with potassium bromide (KBr) beam splitter is implemented for the measurements. The light is focused on the sample through a 36× Cassegrain objective (numerical aperture N.A. = 0.5) and collected onto a liquid N₂-cooled mercury cadmium telluride detector (MCT or Hg:Cd:Te) with a single element detector. The accessible spectral band is ranging between 400 and 8500 cm⁻¹. The acquired spectra are background-corrected to the spectrum of a gold mirror serving as background. Background and spectral measurements are systematically executed using single element configuration with a resolution of 4 cm⁻¹, and interferograms from 1000 scans are typically averaged to obtain one background-corrected spectrum, achieving a high signal-to-noise ratio. Co- and cross-polarized illuminations are subsequently introduced as 0° and 90° angles are respectively read on the polarizer goniometer. Unpolarized light is also implemented by removing the polarizer from the optical pathway.

# 3 Results and discussion

Figure 1 illustrates the MIM structure and the HC array design consisting of equilateral triangles. All the geometrical parameters are denoted in the caption of Figure 1.

A 100 nm mirror layer is sufficiently thick to suppress the light transmission in the IR spectral range. The SiO₂ insulator thickness was determined by FDTD calculations to maximize the E-field enhancement in agreement with

<figure><img src="image_6.png" /><figcaption>**Figure 1:** Schematic of the HC nanoantenna arrays integrated in a MIM structure (Al-HC [50 nm]/SiO₂-spacer [20 nm]/Al-mirror [100 nm]), stacked on a MOS-compatible substrate (SiO₂/Si). L is the equilateral triangle side length; g is the tip-to-tip gap. Each HC is considered as an FDTD unit cell and the periodicity is defined by the boundary conditions. Red, black, and grey arrows illustrate, respectively, the directions of the wavevector of the incident-light, the co-polarized electric field, and the magnetic field.</figcaption></figure>

previous studies of Al and Au BT within a MIM configuration, investigated by Wang et al. [28] and Lin et al. [23], respectively.

## 3.1 Nanoengineering

The HC antennas are fabricated by combing an EBL and a dry Ar etching procedure. By this strategy, HC nanoantennas with controlled L and g are created consisting of hexagonally-arranged equilateral triangles. Overall, L values are varied from 2.0 to 1.0 µm with a step of 0.1 µm; and g from 0.1 to 0.02 µm with a step of 0.02 µm. The technological process flow is summarized in Figure 2 as the following: (A) positive tone (CSAR) resist spin-coating on top of the 50 nm-thick Al layer, (B) EBL linear exposure, (C) exposed-resist development, (D) metallic dry etching by argon (Ar) sputtering for 180 s, (E) unexposed-resist cleaning in a 100 % plasma O₂ environment, and (F) the final structure of 50 nm-thick Al-HC arrangement.

The etching rate is proportional to the opening width that is used to transfer the patterns towards the metallic layer. Hence, the etching step is challenging to optimize as all the *g* values (line width) are engaged on the same sample surface. Nevertheless, to prove that the desired depth of 50 nm is achieved with Ar etching, AFM measurements are done as detailed in Figure 3 for *L* = 2.0 µm and opening of 0.1 and 0.02 µm. Figure 3A underlines the verticality of
<figure><img src="image_7.png" /><figcaption>**Figure 2:** Technological process for HC starts from (A) CSAR resist spin-coating, followed by (B) an EBL exposure and (C) a standard development, then (D) an argon etching and ending with (E) a resist pure plasma O₂ cleaning process, to retrieve (F) the final Al-HC.</figcaption></figure>

the triangle's profiles following the etching step by Ar sputtering. Based on Figure 3C and D, the RMS values of the triangle's surface are alike for both cases and respectively equal 2.5 ± 0.3 nm and 2.6 ± 0.7 nm. To evaluate the etching depth, two cross-sections are sketched in Figure 3C and D along (in red) and across (in black) the tip-to-tip facing triangles. In Figure 3C, a global depth value is found around 55 ± 2 nm through the in-between lines, and a slightly deeper value is found at the intersection point. However, for a smaller opening (Figure 3D), the expected depth is only found at the intersection point. In principle, at the intersection of the exposed lines, the width of the opening is wider than that of the lines between the edges of the triangles. This is due to the triply exposed resist at the intersection sites. Except for these spots, the depth is found around 27 ± 2 nm. Hence, it would have been better to separate the large and the small opening on two different samples and to etch the small openings for a longer time.

Complementary to the AFM images, those collected by SEM highlight the effect of the oxide grains on the final shape of the transferred triangles. As revealed in the SEM images (Figure 4A and B) for $L = 1.0$ µm and $g = 0.1$ µm, a damaging etch occurs due to the oxide grain's resistance against the Ar bombardment. Their resistance is redirecting the Ar ions towards the side walls of the Al triangles. Thus, it is recommended to improve the smoothness of the multi-stacked layers.

Although the triangle's sharpness is not impeccably reached, these HC antennas are robust and fairly adequate to retrieve a good optical behavior compared to the simulation and to successfully cover a wide IR range as proved in the upcoming sections.

## 3.2 Optical characterization

### 3.2.1 Wide spectral coverage

Al-HC is illuminated under a co-polarized (co-pol) plane-wave light source, and the E-field is parallel to the vertical axis of the tip-to-tip facing triangles (black arrow in Figure 1). The reflectance spectra of various Al-HC are simultaneously collected through an xy-plane-monitor placed above the incident source level. Similar to their Al-BT counterparts, the Al-HC are able to cover a wide IR spectral range thanks to their tunable LSPR. Hence, a wide spectral coverage from NIR to MIR is insured using Al-HC from 1000 to 10,000 $cm^{-1}$ as plotted in Figure 5A. In particular, the LSPR of the HC cover the range between 1.0 and 6.6 µm (i.e. 10,000 to 1540 $cm^{-1}$ with a step of ~70 ± 30 $cm^{-1}$), by varying $L$ from 0.3 to 2.0 µm with a step of 0.1 µm. For all HC, the first-order of LSPR is always established at the lowest wavenumber (i.e. lowest energy), and then higher orders of the resonance appear at higher wavenumbers (Figure 5A). This behavior corresponds to the dipolar mode which demonstrates a wavelength shift as $L$ varies. This well-known
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

**Figure 3:** AFM images at the end of the process for *L* = 2.0 µm and *g* equal to (A, C) 0.1 and (B, D) 0.02 µm. Ar etching time is set to 180 s. Color scale intensity is varied between 50 and 70 in (A) and (B) for better visualization. The surface of the Al triangles is normalized likewise in both cases, and the lines definition in (A) is stating a very clean and directional etch. The depth in (C) is standardized around 55 ± 2 nm, but in (D) a deeper etch occurs at the intersection spots, and the lines are not perfectly open.

<figure><img src="image_12.png" /></figure>

**Figure 4:** SEM images of the HC patterns. *L* = 1.0 µm and *g* = 0.1 µm in both (A, B), but the magnification is changed from ×20.0 k (A) to ×40.0 k (B). The triangle finishing is damaged during the Ar etching due to the presence of oxide grains on the sample surface.
<figure><img src="image_13.png" /><figcaption>**Figure 5:** Influence of L value on the LSPR spectral position. (A) Numerical (FDTD) reflectance results for HC illuminated under co-polarized light, where g is kept constant at 100 nm and L is varied from 0.3 to 2 µm with a step of 0.1 µm. Note the spectra are stacked along the reflectance axis for ease of viewing. (B) LSPR wavenumber values are plotted in function of 1/L and fitted by a linear fit (black line), and their corresponding width (FWHM) is showing smaller values as the L increases. The black arrow refers to the LSPR wavenumbers scale (to the left). The red arrow refers to the FWHM values scale (to the right).</figcaption></figure>

tendency for metallic objects is generally defined in the literature [17, 29, 30]. Some higher-order radiative modes are visible only for L > 0.7 µm with a similar tendency. The spectral position of the diffraction peaks Λ_y are resulting from their periodic arrangement which is triggered herein along the y-axis and Λ_y has been analytically deduced in Equation (1) as it follows:

$$ \Lambda_y = L \times \sin 60^{\circ} + \sqrt{3} \times g. \quad (1) $$

Considering the behavior of the LSPR wavelength (λ_res) to red-shift as soon as L increases (Figure 5A), the LSPR wavenumbers ω_res (as proportional to 1/λ_res) are plotted as a function of 1/L in Figure 5B. In analogy to the behavior of a vibrating string, the anticipated linear dependency is confirmed by the automatic linear fit with the coefficient of determination R² = 0.99 and a slope coefficient of 0.2470. Though, as a multiresonant MIM system is hosted, one should discuss the properties of such a plasmonic waveguide. To get more physical insight, the MIM resonator is analyzed as a Fabry-Perot cavity. The thickness of the dielectric (insulator) layer, denoted as t_d, is a key factor in MIM cavities. For a t_d >> λ, the waveguide is formed by uncoupled surface plasmons (SPs) propagating at the two semi-infinite metal/dielectric interfaces. However, as soon as t_d decreases (t_d << λ), the SPs propagating at both interfaces are coupled, i.e. these modes interact more strongly, thus, the effective index n-eff in this narrow cavity rapidly increases [31]. The effective index for the fundamental waveguide mode is given by Equation (2) as:

$$ n_{\text{eff}}^{t_d \ll \lambda} = \epsilon_d^{1/2} \left( 1 - \frac{i\lambda}{\pi t_d \epsilon_m^{1/2}} \sqrt{1 - \frac{\epsilon_d}{\epsilon_m}} \right)^2 \quad (2) $$

For such a subwavelength cavity the effective index of reflection will be written as n eff = ñ eff + i k̃ eff, and with good approximation, ñ eff = ε 1/2 d = 1.3 (SiO₂ as dielectric).

As one of the metallic layers has a finite length (L) and based on the Fabry-Perot phase matching condition, ñ eff is expressed at the resonance wavelength (λ_res) as in Equation (3):

$$ \tilde{n}_{\text{eff}} = \frac{\lambda_{\text{res}}}{2} \frac{1}{L} \left( m - \frac{\phi_r}{\pi} \right) \quad (3) $$

in which, m is an integer equal to 1, 2, ... representing the mode number, and φ_r is the phase of the modal reflection coefficient at the air/MIM interface by having r = √Re iφ_r. The reflection phase φ_r varies depending on the cavity dimensions. As detailed in Equation (3), the coefficient of proportionality between 1/λ_res and 1/L depends on ñ_eff, m and φ_r as 1/λ res = (m - φr/π) 1/L. For the fundamental resonance mode (m = 1), and by recalling the value of the slope coefficient calculated in Figure 5B, the reflection phase φ_r of our multiresonant MIM system is around π/3.

Furthermore, the simulated FWHM (full width at half maximum of the resonance) are revealing a decrease while going deeper in the IR, i.e. for bigger HC antennas (with larger L). The LSPR is evaluated as well through a Q-factor expressed in Equation (4) as:

$$ \text{Q-factor} = \lambda_{\text{res}}/\text{FWHM} \quad (4) $$

where, λ_res is the central wavelength of the LSPR.
For large HC, $1.4 \ \mu\text{m} < L < 2.0 \ \mu\text{m}$, the Q-factor mean value is $4.2 \pm 0.7$. It is higher than the one found for smaller NA for $0.3 \ \mu\text{m} < L < 1.0 \ \mu\text{m}$ (FWHM mean value is $2.5 \pm 0.2$). This reduction of the Q-factor is most likely affected by the Al internal losses (interband transition mode) situated between $\lambda = 0.8$ and $0.9 \ \mu\text{m}$ (fairly close to $10,000 \ \text{cm}^{-1}$). Al losses are due not only to interband transitions but also to Drude scattering. In both cases, the corresponding Q-factors are considered in the same order of magnitude.

The influence of $L$ on the LSPR is also observed using the FTIR under a co-polarized (co-pol) light. Here $L$ is changed from $1.0$ and $2.0 \ \mu\text{m}$ with a step of $0.1 \ \mu\text{m}$. The corresponding LSPR reflectance spectra are plotted in Figure 6. Figure 6A and B show respectively the measured (FTIR) and the simulated reflectance spectra (FDTD). Experimentally, as $L$ increases, the energy of the main LSPR peaks is progressively reduced from $2800$ to $1500 \ \text{cm}^{-1}$ (red-shifting from $3.57$ to $6.6 \ \mu\text{m}$). In particular, while changing the $L$ with a step of $0.1 \ \mu\text{m}$, the LSPR is red-shifted by $\sim 100 \pm 50 \ \text{cm}^{-1}$. A good agreement is deduced by comparing the measured and the simulated LSPR. Though, a slight spectral shift of $\sim 60 \pm 20 \ \text{cm}^{-1}$ is separating the FTIR-resonance positions and the ones retrieved by FDTD. This shift is mainly due to the deposited material quality.

Despite their similitude, the spiky diffraction peaks that we can observe in the simulation spectra, $\Lambda_y$, are suppressed under the FTIR detector having a wide N.A. objective in the microscope. Moreover, in Figure 6A the tangible belongings of the nearby environment, e.g. absorption lines of water and $CO_2$, are witnessed since the optical measurements are acquired in an $N_2$ medium and not in vacuum. The signature of weak oscillations spanning between $1090$ and $1800 \ \text{cm}^{-1}$, $3500$ and $3900 \ \text{cm}^{-1}$, and above $5000 \ \text{cm}^{-1}$ corresponds to water lines. The ro-vibrational modes of $CO_2$ are determined around $2350 \ \text{cm}^{-1}$, and barely envisaged unless it is perfectly intersecting with the resonance of $L = 1.3 \ \mu\text{m}$ situated at $2283 \ \text{cm}^{-1}$.

### 3.2.2 Polarization effect on the LSPR

The influence of the polarization on the HC LSPR is studied as shown in Figure 7. HC with $L = 2.0 \ \mu\text{m}$ and $g = 0.1 \ \mu\text{m}$ are illuminated under co-pol, cross-pol, and un-pol light sources. Both plasmonic modes, i.e. the dipole and multipole, are found insensitive to light polarization. In the simulation, the LSPR positions (black and red circles) of the dipole (at $1630 \ \text{cm}^{-1}$) and the multipole mode (at $2970 \ \text{cm}^{-1}$) under co- and cross-pol light are situated at the same frequency. Similarly, the measured LSPR positions (black and

<figure><img src="image_14.png" /><figcaption>**Figure 6:** Plasmonic activation of the optical transducer over a broad spectral band in the IR according to the side length of the hexagonally placed equilateral triangles in each HC array zone. (A) Experimental (FTIR) versus (B) numerical (FDTD) reflectance results for HC arrays illuminated under co-polarized light (E-field parallel to the axis of the tip-to-tip facing triangles), where $g$ is kept constant at 100 nm and $L$ is varied from 1.0 to $2.0 \ \mu\text{m}$ with a step of $0.1 \ \mu\text{m}$. Note the spectra are stacked along the reflectance axis for ease of viewing.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 7:** Polarization influence on the LSPR of the HC NA illuminated under co-pol, cross-pol, and un-pol light sources. The inset shows the direction of the E-field in each case. FTIR measurements are plotted in lines and the FDTD simulation results are presented in circles.</figcaption></figure>

red lines) of the dipole (at 1535 cm⁻¹) and the multipole mode (at 2755 cm⁻¹) under co- and cross-pol light as well as under an un-pol light (green line) are located at the same frequency. Note that the LSPR spectrum under un-pol light slightly exceeds the spectra under a co-pol and a cross-pol. This is due to the fabrication imperfection as shown in the SEM images of Figure 4. Plus, the measured LSPR are red-shifted compared to the simulated ones as previously discussed. The LSPR amplitude is decreased by 30 %. Both spectral results are the consequence of technological imperfection and the inaccuracy of the material parameters between the simulation and the real measurements.

The results of Figure 7 highlight the rotation invariance of the HC antennas so that the incident light polarization will not affect the position of their corresponding LSPR. This feature will be peculiarly useful to improve the signal-to-noise ratio, especially for the subsequently performed sensing experiments, since the totality of the incident light intensity reaches the zone of interest onto the sample without being attenuated by a factor of 1/2 upon navigating through a polarizer.

### 3.2.3 Local near-field enhancement

Besides the far-field interpretation, the LSPR are accompanied by strong near-field localization and enhancement at the triangle's pointed extremities forming a hotspot [11]. Indeed, the MIM cavity absorbs the incoming light at the

resonance, yielding a large field enhancement because the field gets “squeezed” into the cavity at the edges of the finite antenna, resulting in a lightning-rod effect. By promoting the strong interaction between the oscillating charges on the closely-placed sharp apexes, the near E-field intensity, $|E|^2/|E_0|^2$, is accordingly improved inside the activated gap. Such a coupling is forced to occur along the polarization direction. Consequently, under a co-polarized light, the near E-field enhancement and intensity are boosted along the tip-to-tip facing triangles and are affected by this gap value.

In the following, the gap influence on the plasmonic response of the Al-HC is investigated under a co-pol light. The structure's side length is $L = 0.8$ µm and the gaps $g$ are selected to be equal to 0.1, 0.06, and 0.02 µm. In principle, based on our previous work regarding the BT [18], the LSPR red-shifts, and the $|E|^2/|E_0|^2$ increases as soon as the value of $g$ decreases. To evoke the E-field dependency on the $g$ value, a finer mesh override of 2 nm is adapted instead of 5 nm for more accurate E-field mapping. The E-field intensity is collected in the xy-plane and represented in Figure 8. Under a co-pol light, the E-field is pinched at the extremities of the tip-to-tip facing triangles.

For $g = 0.1$ µm (Figure 8A), the max $|E|^2/|E_0|^2$ equals 250 and appears at $3030.3$ cm⁻¹ ($\lambda = 3.3$ µm). It is slightly red-shifted with respect to the far-field LSPR response, i.e. reflectance minimum found at $3336.7$ cm⁻¹ ($\lambda_{res} = 2.97$ µm). The maximum near E-field is found always at lower energy than the far-field LSPR energy because our plasmonic resonators are conceived as damped harmonic oscillators with natural material losses.

Furthermore, the max value of $|E|^2/|E_0|^2$ should increase as far as $g$ decreases. However, for a $g = 0.06$ µm, LSPR energetic position appears at $3278.7$ cm⁻¹ ($\lambda_{res} = 3.05$ µm), and $|E|^2/|E_0|^2$ equal to 233 is found at lower energy, $2967.35$ cm⁻¹, (i.e. $\lambda = 3.37$ µm) that is slightly lower than the one found for $g = 0.1$ µm. This result seems illogical until looking at the corresponding E-field mapping in Figure 8B that shows a considerable distribution along the triangle's neighboring edges. This behavior is emphasized when having $g = 0.02$ µm (Figure 8C). Herein, a max $|E|^2/|E_0|^2$ of 216 is found at $2857.15$ cm⁻¹ ($\lambda = 3.5$ µm) (red-shifted in respect to the far-field LSPR found at $2906.97$ cm⁻¹, i.e. $\lambda_{res} = 3.44$ µm) and is lower than the one found for a larger gap. In fact, for such a small gap, the near E-field is not only enhanced at the sharp tips but also at the corners of the triangles. This is lowering the possibility of the near E-field being largely enhanced by only being pinched at the sharp facing tips. As shown in Figure 8C, this near E-field distribution creates a kind of long-way cavity where the near E-field could be probably delocalized.
<figure><img src="image_16.png" /><figcaption>**Figure 8:** Gap value influence on the local field enhancement and intensity for *L* = 0.8 µm with *g* equal to 0.1, 0.06, and 0.02 µm. Under a co-pol light, the enhanced electric field mapping is collected in the *xy*-plane (A) for *g* = 0.1 µm calculated at *λ* = 3.3 µm, (B) for *g* = 0.06 µm calculated at *λ* = 3.37 µm, (C) for *g* = 0.02 µm calculated at *λ* = 3.5 µm. The color scale ranging from blue to red indicates the minimum to maximum field intensity |*E*|²/|*E*₀|², max (|*E*|²/|*E*₀|²) equals 250, 233, 216 while decreasing the *g* value. For better contrast, the maximum color scale is fixed at 50. All scale bars are 0.1 µm long.</figcaption></figure>

In comparison to the Al-BT structure, for $g = 0.02$ µm, $|E|^2/|E_0|^2$ of 1250 is found instead of 216, so, in terms of effective E-field enhancement, the BT structure is superior to the HC structure. However, in terms of near E-field enhanced surface, which is vital for molecule-antenna coupling, such a long-way cavity would boost the SEIRA outcome, *e.g.* the enhanced surface in the enhancement factor calculation. In addition, the enhanced near E-field using HC is driven by the light polarization and thus localized at more than 2 tips (surpassing the BT) by turning the polarization direction. Herewith, the enhanced selected-sites will be as well beneficial for sensitive SEIRA sensing.

### 3.2.4 Arrangement influence on the LSPR

HC structure might be more interesting than a BT antenna by affording a higher signal-to-noise ratio (SNR). At this stage, a brief comparison between both BT and HC nanoantennas is worth mentioning, as presented in Figure 9. Here, *L* is fixed at 2.0 µm and *g* at 0.1 µm and a co-pol light is applied. For a HC compact arrangement of the triangles (pink curve), one can remark that the amplitude of the LSPR is improved by 20 % compared to BT as simple tip-to-tip triangle arrays (brown curve). This underlines the advantage of the compact nanoantennas organization. Moreover, the compact arrangement is translated into a slight red-shift of the LSPR as the interaction between the adjacent triangles becomes more effective, and thus, arises at lower energy. The LSPR of the HC has shown a slight red-shift ($\Delta\omega_{res} \sim 20 \pm 4 \text{ cm}^{-1}$) compared to the BT resonance for the same *L* and *g* values. The interpretation of the measured spectra agrees fairly well with the simulation results taking into account the fabrication restrictions.

Plus, the LSPR spectra are large enough, thus, are beneficial for the upcoming sensing application that anticipates

<figure><img src="image_17.png" /><figcaption>**Figure 9:** The arrangement influence on the reflectance of both BT (brown curve) and HC (pink curve) nanoantennas (insets), illuminated under a co-pol light source. In both cases, *L* and *g* are respectively fixed at 2.0 µm and 0.1 µm. For a compact arrangement, the LSPR is more pronounced in terms of amplitude and it is slightly red-shifted compared to the one of the BT, as more resonators are interacting.</figcaption></figure>

a hybrid coupling between the wide $\omega_{res}$ and narrow $\omega_{vib}$ frequencies.

## 4 Sensing demonstration

The SEIRA investigation is first presented from a far-field perspective in Figure 10. As a starting point, the vanillin absorption lines from 1100 to 2100 cm⁻¹ are obtained from a transmittance measurement of a highly-concentrated vanillin solution (50 mg/mL) in IPA, deposited on an IR-transparent KBr crystal, as revealed in Figure 10A.

To evaluate the features of the SEIRA spectroscopy using HC nanoantennas, a 2 mg/mL solution of vanillin/IPA is prepared. A 1 µL droplet is then released onto the HC sample surface. Around 10 nmol are hence dropped on
<figure><img src="image_18.png" /><figcaption>**Figure 10:** Resonant SEIRA sensing demonstration of vanillin using Al-HC under a co-pol light. The gap value is fixed at 0.1 µm. (A) Transmittance spectrum of vanillin (50 mg/mL) normalized to the one of cleaned KBr platelet. Gray vertical lines eye-guide the IR absorption lines of vanillin. Al-HC of different *L* varied from 2.0 to 1.4 µm with a step of 0.1 µm are covering a wide spectral range from 1400 to 2100 cm⁻¹. The corresponding reflectance spectra are acquired before (B) and after (C) the deposition of vanillin. The dashed curves are the corresponding baseline (BL) fitted using Eiler's smoothing algorithm, excluding the vibrational features (Fano-profiles). (D) Vibrational signals are normalized by the BL for each HC array (side length *L* given in microns). The plasmonic resonance positions *ω*<sub>res</sub> are marked with thick-colored vertical lines. The signal strengths of the enhanced molecular vibrational modes are obtained as a peak-to-peak value defined by the horizontal ticks.</figcaption></figure>

the whole transducer surface (having the vanillin with a
molecular mass of 152.1 g/mol). The distribution of vanillin
molecules is assumed to be homogenous on the 1 cm² sam-
ple surface. Indeed, the adsorption of vanillin on Al follows
a Langmuir isotherm model where the adsorbed molecules
do not interact and mono-layer coverage could be assumed
once the vanillin is dropped onto the surface [32]. So, 0.2 ng
of vanillin are located within the zone of 100 × 100 µm².

The SEIRA behavior using the HC nanoantennas is val-
ued once far-field measurements are done using the single-
element detector under a 36 × IR objective (NA = 0.5).
The reflectance area of the optical transducer is selected
by a manual knife-edge aperture delimiting the 100 ×
100 µm² tailored zone full of HC patterns. Figure 10 presents
the achieved reflectance measurements under a co-pol

light, on the same sample position without (Figure 10B)
and in presence of (Figure 10C) the vanillin/IPA (2 mg/mL)
droplet. The coupling between the wide ωres and narrow
vibrational modes ωvib frequencies is studied through the
obtained vibrational signals (Figure 10D), known as asym-
metric Fano-like profiles. Five narrow IR absorption lines
of the vanillin for wavenumbers higher than 1400 cm⁻¹
are interfering within broad ωres corresponding to Al-HC
arrangement with L varying from 2.0 to 1.4 µm with a step
of 0.1 µm. The five absorption features of vanillin are sit-
uated at 1429.8, 1464.1, 1508.9, 1587.5, and 1664.6 cm⁻¹ (fol-
lowing the vertical gray lines along Figure 10). These lines
correspond closely to C=O, CH, CC, CO, and CH₃ stretching,
in-plane bending, and asymmetric deformation vibrational
lines [33]. Note, the interaction of the molecules with the
substrate and their medium (solid and liquid) is directly affecting their binding energy. In line with this, the experimentally found vibrational lines are slightly shifted compared with the analytical ones.

In Figure 10D, the vibrational signals are calculated thanks to the Eiler's least square smoothing algorithm [34], by normalizing the reflectance curve for each HC array to their corresponding baseline (BL, dashed lines in Figure 10C). The efficiency of the occurring coupling is interpreted through the strength (peak-to-peak value) and the shape of the Fano-like vibrational signals and is modulated by a tuning ratio $\omega_{vib}/\omega_{res}$. The vibrational signal strength tends to be maximized for the different L values once the tuning ratio $\omega_{vib}/\omega_{res}$ is around 1 (typically around the colored vertical lines) as shown in Figure 10D. The maximum amplitude of the vibrational signal states the best spectral overlap between $\omega_{res}$ and $\omega_{vib}$. Plus, the tuning ratio also affects the shape of the vibrational signals. For example, selecting $\omega_{vib} = 1664 \text{ cm}^{-1}$, for coupling tunings above 1 (Figure 10D red curve), the line shape keeps on its asymmetric silhouette but reverses its minimum and maximum compared to the tunings below 1 (Figure 10D blue curve) [9]. It is crucial to mention that similar results are found under un-pol light.

To evaluate the SEIRA sensitivity of the Al-HC arrays, the enhancement factor, EF, for the five IR lines of the vanillin were examined according to the spectral tuning ratio. Equation (5) defines EF as the following:

$$EF = (S_{SEIRA}/S_0) \times (A_0/A_{SEIRA}) \quad (5)$$

where $A_{SEIRA}$ and $A_0$ denote, respectively, the areas covered with molecules in SEIRA and the reference measurements. $S_{SEIRA}$ and $S_0$ are, respectively, the enhanced and unenhanced signal strengths. From the measurements on a metallic surface, the unenhanced signal strengths, $S_0$ are found in the order of $2 - 5 \times 10^{-3}$ for the absorption lines. The signal strengths, $S_{SEIRA}$, are in the order of $2.5 \times 10^{-2}$, as measured in Figure 10. This value is slightly higher than the retrieved values for Al-BT under a co-pol light ($S_{SEIRA}$ in the order of $2 \times 10^{-2}$), and it is related to a signal improvement thanks to compacting the nanoantennas into the characterized zone. Under un-pol light, $S_{SEIRA}$ are in the order of $4 \times 10^{-2}$, two times greater than the $S_{SEIRA}$ values measured under a co-pol light, and it is due to the improvement of the SNR once the totality of the incident illumination intensity reaches the zone of interest onto the sample without being attenuated.

$A_0$ is defined as the repetitive surface entity (unit cell) and is analytically calculated as developed in Equation (6):

$$A_0 = 3 \times \sqrt{3}/2 \times a^2; \quad a = L + \sqrt{3} \times g \quad (6)$$

The total antenna area, $A_0$, yields 0.2 pg per HC antenna, the half of the mass that is covering the $A_0$ for BT, as the characterized surface of HC is half of the BT's one ($S_{objective,HC} = 100 \times 100 \text{ µm}^2 < S_{objective,BT} = 230 \times 90 \text{ µm}^2$). For an entirely covered antenna, $7 \times 10^8$ molecules per antenna (lower than $15 \times 10^8$ molecules per BT antenna) are estimated. However, the probed molecules are outlined in the antenna intense near-field hotspots triggered by the light polarization (see the inset of Figure 11). Hence, an effective area, $A_{eff}$, that is accessible to the analyte is defined as the 5/6 of an 8 nm-radius circular surface centered around the tip, where the local E-field is intensified, excluding thus the internal area of the equilateral triangle tip. Under a co-pol light, $A_{SEIRA}$ is counting 6 times the $A_{eff}$. $A_{SEIRA}$ represents 0.12% of the unit cell area $A_0$. Therefore, approximately $1.3 \times 10^4$ molecules are triggered per tip ($A_{eff}$). Under un-pol light, the same amount of molecules is retrieved per tip, yet, $A_{SEIRA}$ is counting 18 times the $A_{eff}$. It is three orders of magnitude lower than the value reported by Barho et al. for vanillin detection (~$5 \times 10^7$ molecules per antenna) with highly doped semiconductor periodic rectangular nanoantenna arrays [30]. $A_{SEIRA}$ represents 0.37% of the unit cell area $A_0$. Despite

<figure><img src="image_19.png" /><figcaption>Figure 11: Overview of the normalized enhancement factor (EF) of the five absorption features of vanillin versus the tuning ratio $\omega_{vib}/\omega_{res}$. Under a co-pol light, the EF are studied for HC with a $g = 0.1 \text{ µm}$. The inset shows both $A_0$ and $A_{SEIRA}$ under a co-pol light. Note, under un-pol light all the triangles tips will be activated. The vibrational signal enhancement is generally maximized for a tuning ratio below 1 (black line) as found with a Lorentz fit (LF), declaring the perfect overlap between vibrational modes and plasmonic resonances. The thick error bar represents the systematic error related to $A_{SEIRA}$ estimation, to be introduced to the maximum EF value. Additional EF sequences are presented for HC with $g = 0.1$ and 0.08 µm illuminated under un-pol light. The maximum EF is reduced by a factor of 3 under un-pol light which is inversely proportional to the $A_{SEIRA}$ prompted by the light polarization.</figcaption></figure>
of the difference in $A_0$ definition and $A_{SEIRA}$ quantification, the retrieved numbers of triggered molecules using Al-HC antennas are in the same order of magnitude as for the Al-BT antennas [18]. It is successfully matching the state-of-art as reported for self-assembled monolayers detection by Halas et al. [22], Neubrech et al. [5], and Brown et al. [14].

Moreover, the normalized enhancement factor *EF* depends on the tuning ratio $\omega_{vib}/\omega_{res}$ as illustrated in Figure 11 for all the selected *L* values. *EF* is higher than $\sim 10^7$ under a co-pol light as shown in the y-axis of Figure 11. It is identical to the order of magnitude achieved with the Al-BT antennas and it is agreeing with the greatest stated value reported by Dong et al. using 3-nm gapped BT antennas onto an optical cavity [15]. The *EF* maximum value of $1.2 \times 10^7$ is found for a tuning ratio $\omega_{vib}/\omega_{res} = 0.98$. It is achieved for the absorption line $\omega_{vib}$ situated at $1587 \text{ cm}^{-1}$ with an Al-HC array of $L = 1.9 \text{ µm}$ resonating at a $\omega_{vib} = 1533 \text{ cm}^{-1}$. It is four orders of magnitude higher than the reported value by Vogt et al. for IR sensing (*EF* $\sim 10^3$ using Au nanorods) [9].

Under un-pol light, the maximum *EF* is established around $3 \times 10^6$ for $g = 0.1 \text{ µm}$ (black square in Figure 11) and around $3.6 \times 10^6$ for $g = 0.08 \text{ µm}$ (red circle in Figure 11). It is reduced by a factor of 3 compared to the previously mentioned *EF*. From a near-field point of view, this result seems to correlate with Equation (5) as *EF* is inversely proportional to $A_{SEIRA}$. Hence, the maximum *EF* is reduced under un-pol light (triggering 18 tips) versus the one retrieved under a co-pol light (triggering 6 tips). Moreover, despite of the SNR improvement by a factor of 2 when removing the polarizer, the anticipated results under a cross-pol light should be taken into consideration as the un-pol light is embracing both polarization cases. Here, the hypothesis reclaims a drop-off in the *EF* values under cross-pol light that is translated into a total defeat of those under un-pol light. Based on a simulation under cross-pol, the near E-field intensity, where the SEIRA mainly arises, is about half of the intensity value that is triggered under a co-pol light. This near-field result supports the previous hypothesis and justifies the *EF* slight dropping under an un-pol light.

# 5 Conclusion and perspective

In this work, a multiresonant optical transducer of Al-HC is integrated within a MIM configuration. The nanonan- tennas are tailored by combining a novel patterning tac- tic using electron-beam lithography and a challenging dry etching. The dependency of the nanoantenna's resonance ($\omega_{res}$) on their geometric parameters, e.g. the size and the gap spacing, is studied in detail. A linear dependency is found between $\omega_{res}$ and $1/L$. By considering the Fabry-Perot

phase matching condition, a reflection phase $\phi_r$ around $\frac{\pi}{3}$ is found for the fundamental mode of our MIM system. By varying the size of the HC nanoantennas, a wide IR spectral range is successfully covered. The far-field optical response of Al-HC has been numerically modeled and exper- imentally verified. Good agreement was found between both approaches. Thanks to the symmetry invariance of this hexagonal arrangement, its corresponding resonance is proved to be insensitive to the light polarization, lead- ing thus to an improvement of the SNR by a factor of 2 while removing the light polarizer. Nevertheless, the near- field localization will be driven by the light polarization. Accordingly, we can conceive a plasmon-mediated and site- selective immobilization of biomolecules (DNA, proteins, amino acids) [35] and nanoparticles (quantum dots) [36]. Besides this appealing feature, the HC arrangement, due to its compactness, surpasses the BT arrays by exhibiting an extra increase of the SNR ratio. To conclude, HC nanoan- tennas are well-adapted for an accurate vibrational SEIRA sensing demonstration. Indeed, the far-field and near-field results are matching the SEIRA state-of-art. In addition, the arrangement of the antennas boosts their reflectance ampli- tude, and in consequence, the amplitude of the vibrational signals $S_{SEIRA}$. Under a co-pol light, an enhancement factor, *EF*, higher than $10^7$ is achieved at a tuning ratio below 1 and is perfectly matching the state-of-art of SEIRA using Al-BT. Under un-pol light, despite of the SNR improvement, the *EF* is decreased by a factor of 3 due to the increase of $A_{SEIRA}$ (18 triggered hotspots per HC antenna). This valuable activation of $A_{SEIRA}$ stimulates an encouraging analysis of a minimal detectable concentration of our analyte. Finally, to improve the SEIRA readout, smaller gaps can be finely designed to offer a greater near-field enhancement inside the gap as shown theoretically. In principle, a stronger E-field intensity and a highly “hotspot-populated” HC unit will impeccably lead to higher SEIRA enhancements. Hence, our propitious Al-HC nanoantennas within a MIM configuration antici- pates a mass production of site-selective, sensitive, compact and cost-effective biosensing devices for delicate and broad- band IR detection based on smart Silicon technologies.

**Author contributions:** All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.

**Research funding:** This work was achieved in collab- oration with the micro and nanotechnologies platform at the LAAS-CNRS, a member of the Renatech French national network. M.N. acknowledges a doctoral grant from the National French Ministry for University Education, Research, and Innovation. Additional funding is acknowl- edged from the following projects: SEA (Occitanie French
Region – ESR-PREMAT-238), EXTRA (ANR 11-EQPX-0016), NanoElastir (ASTRID 2020–2023), and ENVIRODISORDERS (MUSE UM 2021–2023). The authors acknowledge technical support from Delphine Jourdain, Ariane Meguekam, Jean-Marc Aniel, and Guilhem Boissier. The authors acknowledge the technical work of the cleanroom staff, Frederic Pichot, Jean-Marie Peiris, and Renaud Felix at Montpellier Technological Center (CTM). The authors thank Michel Ramonda for the AFM measurements at the CTM. The authors would like to extend the thanks to Dr. Mikhael Bechelany from the IEM of Montpellier for the SEM images. The authors thank Dr. Frank Neubrech for sharing his Eiler's smoothing program for the evaluation of SEIRA signals and Prof. Annemarie Pucci for the fruitful discussion. The authors thank Dr. Aude Mezy and Dr. Celine Schmitt for collaborating on the IR spectral lines interpretation.

**Conflict of interest statement:** The authors declare no conflicts of interest regarding this article.

## References

[1] M. Osawa, "In-situ surface-enhanced infrared spectroscopy of the electrode/solution interface," in *Advances in Electrochemical Science and Engineering*, Weinheim, Wiley VCH, 2006, pp. 269–314.

[2] K. Ataka, T. Kottke, J. T. Heberle, and F. Smaller, "IR techniques to probe the functionality of biological and biomimetic systems," *Angew. Chem. Int. Ed.*, vol. 49, no. 32, pp. 5416–5424, 2010.

[3] H. L. Wang, E. M. You, R. Panneerselvam, S. Y. Ding, and Z. Q. Tian, "Advances of surface-enhanced Raman and IR spectroscopies: from nano/microstructures to macro-optical design," *Light Sci. Appl.*, vol. 10, no. 1, pp. 1–19, 2021.

[4] M. Osawa and M. Ikeda, "Surface-enhanced infrared absorption of p-nitrobenzoic acid deposited on silver island films: contributions of electromagnetic and chemical mechanisms," *J. Phys. Chem.*, vol. 95, no. 24, pp. 9914–9919, 1991.

[5] F. Neubrech, C. Huck, K. Weber, A. Pucci, and H. Giessen, "Surface-enhanced infrared spectroscopy using resonant nanoantennas," *Chem. Rev.*, vol. 117, no. 7, pp. 5110–5145, 2017.

[6] S. A. Maier, *Plasmonics: Fundamentals and Applications*, New York, Springer, 2007.

[7] A. Tittl, A. John-Herpin, A. Leitis, E. R. Arvelo, and H. Altug, "Metasurface-based molecular biosensing aided by artificial intelligence," *Angew. Chem. Int. Ed.*, vol. 58, no. 42, pp. 14810–14822, 2019.

[8] R. Adato, A. Artar, S. Erramilli, and H. Altug, "Engineered absorption enhancement and induced transparency in coupled molecular and plasmonic resonator systems," *Nano Lett.*, vol. 13, no. 6, pp. 2584–2591, 2013.

[9] J. Vogt, C. Huck, F. Neubrech, A. Toma, D. Gerbert, and A. Pucci, "Impact of the plasmonic near- and far-field resonance-energy shift on the enhancement of infrared vibrational signals," *Phys. Chem. Chem. Phys.*, vol. 17, no. 33, pp. 21169–21175, 2014.

[10] U. Fano, "The theory of anomalous diffraction gratings and of quasi-stationary waves on metallic surfaces (sommerfeld's waves)," *J. Opt. Soc. Am.*, vol. 31, no. 3, pp. 213–222, 1941.

[11] S. A. Maier and H. A. Atwater, "Plasmonics: localization and guiding of electromagnetic energy in metal/dielectric structures," *Appl. Phys. Rev. – Focus. Rev.*, vol. 98, no. 011101, pp. 1–10, 2005.

[12] M. Urbieta, M. Barbry, Y. Zhang, et al., "Atomic-scale lightning rod effect in plasmonic picocavities: a classical view to a quantum effect," *ACS Nano*, vol. 12, no. 1, pp. 585–595, 2018.

[13] L. V. Brown, K. Zhao, N. King, H. Sobhani, P. Nordlander, and N. J. Halas, "Surface-enhanced infrared absorption using individual cross antennas tailored to chemical moieties," *J. Am. Chem. Soc.*, vol. 135, no. 9, pp. 3688–3695, 2013.

[14] L. V. Brown, X. Yang, K. Zhao, B. Y. Zheng, P. Nordlander, and N. J. Halas, "Fan-shaped gold nanoantennas above reflective substrates for surface-enhanced infrared absorption (SEIRA)," *Nano Lett.*, vol. 15, no. 2, pp. 1272–1280, 2015.

[15] L. Dong, X. Yang, C. Zhang, et al., "Nanogapped Au antennas for ultrasensitive surface-enhanced infrared absorption spectroscopy," *Nano Lett.*, vol. 17, no. 9, pp. 5768–5774, 2017.

[16] G. Q. Wallace, H. C. Foy, S. M. Rosendahl, and F. Lagugné-Labarthet, "Dendritic plasmonics for mid-infrared spectroscopy," *J. Phys. Chem. C*, vol. 121, no. 17, pp. 9497–9507, 2017.

[17] P. Biagioni, J. S. Huang, and B. Hecht, "Nanoantennas for visible and infrared radiation," *Rep. Prog. Phys.*, vol. 75, no. 2, pp. 1–40, 2012.

[18] M. Najem, F. Carcenac, T. Taliercio, and F. Gonzalez-posada, "Aluminum bowties for plasmonic-enhanced infrared sensing," *Adv. Opt. Mater.*, vol. 10, no. 13, pp. 1–12, 2022.

[19] A. E. Cetin, D. Etezadi, and H. Altug, "Accessible nearfields by nanoantennas on nanopedestals for ultrasensitive vibrational spectroscopy," *Adv. Opt. Mater.*, vol. 2, no. 9, pp. 866–872, 2014.

[20] C. Kittel, *Introduction to Solid State Physics*, 8th ed., New York, Wiley, 2005.

[21] J. Martin and J. Plain, "Fabrication of aluminium nanostructures for plasmonics," *J. Phys. D. Appl. Phys.*, vol. 48, no. 18, p. 184002, 2015.

[22] B. Cerjan, X. Yang, P. Nordlander, and N. J. Halas, "Asymmetric aluminum antennas for self-calibrating surface-enhanced infrared absorption spectroscopy," *ACS Photonics*, vol. 3, no. 3, pp. 354–360, 2016.

[23] L. Lin and Y. Zheng, "Optimizing plasmonic nanoantennas via coordinated multiple coupling," *Sci. Rep.*, vol. 5, no. 14788, pp. 1–11, 2015.

[24] R. Nicolas, G. Lévêque, J. Marae-Djouda, et al., "Plasmonic mode interferences and Fano resonances in metal-insulator- metal nanostructured interface," *Sci. Rep.*, vol. 5, no. 144119, pp. 1–11, 2015.

[25] E. Aslan, E. Aslan, O. G. Saracoglu, and M. Turkmen, "An effective triple-band enhanced-infrared-absorption detection by honeycomb-shaped metamaterial-plasmonic absorber," *Sensors Actuators, A Phys.*, vol. 288, pp. 149–155, 2019.

[26] E. Shiles, T. Sasaki, M. Inokuti, and D. Y. Smith, "Self-consistency and sum-rule tests in the Kramers–Kronig analysis of optical data: applications to aluminum," *Phys. Rev. B*, vol. 22, no. 4, pp. 1612–1628, 1980.

[27] E. D. Palik, *Handbook of Optical Constants of Solids*, vol. 2, San Diego, USA, Academic Press, 1991.
[28] B. Wang, S. C. Singh, H. Lu, and C. Guo, "Design of aluminum bowtie nanoantenna array with geometrical control to tune LSPR from UV to near-IR for optical sensing," *Plasmonics*, vol. 15, no. 3, pp. 609–621, 2020.

[29] M. J. Milla, F. Barho, F. González-Posada, et al., "Localized surface plasmon resonance frequency tuning in highly doped InAsSb/GaSb one-dimensional nanostructures," *Nanotechnology*, vol. 27, no. 42, pp. 1–8, 2016.

[30] F. B. Barho, F. Gonzalez-Posada, M. J. Milla, et al., "Highly doped semiconductor plasmonic nanoantenna arrays for polarization selective broadband surface-enhanced infrared absorption spectroscopy of vanillin," *Nanophotonics*, vol. 7, no. 2, pp. 507–516, 2017.

[31] J. Yang, C. Sauvan, A. Jouanin, S. Collin, J.-L. Pelouard, and P. Lalanne, "Ultrasmall metal-insulator-metal nanoresonators: impact of slow-wave effects on the quality factor," *Opt. Express*, vol. 20, no. 15, p. 16880, 2012.

[32] A. Y. El-Etre, "Inhibition of acid corrosion of aluminum using vanillin," *Corros. Sci.*, vol. 43, no. 6, pp. 1031–1039, 2001.

[33] V. Balachandran and K. Parimala, "Vanillin and isovanillin: comparative vibrational spectroscopic studies, conformational stability and NLO properties by density functional theory calculations," *Spectrochim. Acta A Mol. Biomol. Spectrosc.*, vol. 95, nos. 1386–1425, pp. 354–368, 2012.

[34] P. H. C. Eilers, "A perfect smoother," *Anal. Chem.*, vol. 75, no. 14, pp. 3631–3636, 2003.

[35] L. N. Liu and T. Liedl, "DNA-assembled advanced plasmonic architectures," *Chem. Rev.*, vol. 118, no. 6, pp. 3032–3053, 2018.

[36] I. Kherbouche, Y. Luo, N. Félidj, and C. Mangeney, "Plasmon-mediated surface functionalization: new horizons for the control of surface chemistry on the nanoscale," *Chem. Mater.*, vol. 32, no. 13, pp. 5442–5454, 2020.