Preprint
Article

Numerical Investigation of Localized Surface Plasmons in Gold Nanoridges Dimer-on-Mirror

Altmetrics

Downloads

88

Views

30

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

13 August 2024

Posted:

13 August 2024

You are already at the latest version

Alerts
Abstract
The study of localized surface plasmons (LSPs) in nanoscale structures is an essential step towards identifying optimal plasmonic modes that can facilitate robust optomechanical coupling and deepen our understanding of light-matter interactions at the nanoscale. This paper investigates numerically, using the finite element method, LSPs modes in a design comprising two coupled nano-ridges deposited on a gold layer, with an interposing polymer spacer layer. Such a structure, usually referred to as particles-on-mirror, shows exquisite optical properties at the nanoscale. We first examine the LSPs modes of a single nano-ridge through the analysis of its scattering cross-section in the visible and infrared range. To enhance the plasmonic response, a thin polymer layer is placed at the middle of the ridge, which introduces additional LSPs modes confined within the former. Then, we extend the analysis to the dimer configuration, which exhibits more complex and enhanced plasmonic behavior compared to a single nano-ridge. In particular, the dimer configuration yields LSPs resonances with a quality factor enhancement of approximately threefold relative to single nano-ridge. Furthermore, the presence of the polymer layer within the ridges significantly improves plasmon field localization and quality factor. These findings underscore the potential of nano-ridge-based structures in advancing optomechanical coupling and plasmon-phonon interaction studies, offering valuable insights for the development of high-performance acousto-plasmonic devices. In particular, the proposed device could help significantly improve the design of nano-acousto-optic modulators operating in the visible or in the near infrared, that require enhanced light-photon coupling rate.
Keywords: 
Subject: Physical Sciences  -   Optics and Photonics

1. Introduction

Plasmonics and its applications is based on the optical processes that occur in metallic nanostructures involving the coupling of light and collective oscillations of conduction band electrons (known as plasmons) at the metal–dielectric interface. The size, shape, and composition of the nanostructures can be engineered to manipulate the properties of the plasmons, such as their resonance frequency and localization, which, in turn, can enable a wide range of applications in various fields [1,2,3,4,5].
In addition, the application of plasmonics represents a good opportunity for nano-optomechanics to explore the interaction between light and mechanical motion at the nanoscale, enabling ultra-sensitive mass detection, precision measurement of forces, and advancements in quantum information processing [6,7,8,9]. These systems typically utilize nanoscale resonators where light controls and measures mechanical vibrations. Achieving strong optomechanical coupling relies on the ability to localize and enhance electromagnetic fields at the nanoscale. Metallic nanostructures possess the ability to confine light below the diffraction limit via the excitation of localized surface plasmons (LSPs), the oscillations of free electrons at the surface of metallic nanostructures, offering a powerful means of confining light to sub-wavelength dimensions [10,11]. LSPs can interact with mechanical vibrations, making them ideal for studying nanoscale mechanical oscillations and enhancing optomechanical interactions [12,13,14]. Investigating systems with high LSP confinement is crucial for advancing nano-optomechanical applications.
Among the nanostructures supporting highly confined LSPs, one can cite nanoparticles-on-mirror (NPsOM) (or film-coupled nanoparticles) [15,16], where a metal nanoparticle is separated from a metal film by a thin dielectric layer acting as a spacer; these structures have attracted significant attention due to their ease of fabrication [17] and the ability to trap plasmon modes in ultrasmall volumes [18,19,20]. The response of NPsOM can be significantly influenced by the thickness of the metal film [21,22], as thinner metal films tend to result in stronger coupling between the nanoparticle and the film, leading to more pronounced plasmonic effects, while thicker films may reduce this interaction and broaden the resonance linewidths. Similarly, the interspacing between the gold film and the nanoparticle [23] and the shape of the nanoparticles plays a crucial role in determining the plasmonic properties [24]. Also, NPsOM with dimers or trimers atop the metal film, as opposed to monomers, are expected to exhibit more interesting and/or improved features such as narrower plasmon linewidths and more pronounced scattered intensity [13,25].
This study focuses on investigating numerically based on the finite element method (FEM), LSPs in a structure design in the category of NPsOM devices, consisting of nano-ridges deposited on a gold layer, with a polymer layer acting as a spacer. It is worth mentioning that such a design could be readily fabricated experimentally given the state-of-the-art technology [26]. We first investigate the LSPs modes of a single ridge by examining its scattering cross section in the visible and near-infrared range. To introduce additional LSPs modes, we insert a thin polymer layer at the middle of the ridge, which results in new LSPs confined within such a layer. Then, we explore the dimer configuration, which has been proven to show more interesting characteristics [13,25]. Our findings indicate that dimerization, achieved by positioning two ridges in close proximity, enhances the interaction between the modes, leading to LSPs resonances with significantly improved quality factors—approximately a threefold enhancement compared to a single monomer. Additionally, the polymer layer at the middle of the ridges plays an important role in enhancing the localization of the plasmon field and the quality factor of the different modes as compared to previous reports [13]. In a recent paper [27], we have studied localized phonons in a similar monomer and dimer structures. Here, we show prospectively based on the dimer LSPs modes effective volumes, how such modes would couple with the localized phonons in the dimer, thereby enabling their classification in terms of optomechanical coupling intensity.
This paper is organized as follows: In Section II, we present a detailed investigation of the LSPs of a single gold nano-ridge positioned atop a multilayer structure. In Sec. III, we extend this analysis to a coupled gold nanoridges, examining how the dimer configuration influences the localization and enhancement of the plasmonic modes. In Sec. IV, we provide an analysis of the mode volume for each LSP mode, enabling the classification of the different LSP modes and discussing their potential applications in acousto-plasmonic interactions.

2. Localized Surface Plasmons of a Single Gold Nano-Ridge

The proposed design is shown in Figure 1a. It consists of a gold nano-ridge with a thin polymer layer lying in the middle, deposited above a polymer thin spacer layer made of the same material, that covers a gold thick film beneath. Experimentally speaking, the whole design could be grown on a silicon dioxide substrate (SiO2), but in terms of simulations the thick gold film absorbs and reflects all incident light before reaching SiO2. Hence, the latter medium is excluded from simulations. Due to the system’s symmetry along the y-axis, all calculations were performed in the x-z plane.
The modeled cross-section of the system is depicted in Figure 1a. The parameters of the studied system are as follows: W is the nano-ridge width while L is its height, dpol is the thickness of the polymer layer inserted in the middle of the ridge, Lpol is the thickness of the polymer spacer layer and LAu is the thickness of the gold layer beneath. The simulation results presented in this work were obtained using electromagnetic frequency domain analysis based on the FEM, using COMSOL Multiphysics software. In order to excite the system, a plane light-wave with TM polarization is launched from the incidence medium at the top, at oblique incidence onto the latter, as shown in Figure 1a. To eliminate undesirable reflections at the edges of the system, perfectly matched layers (PMLs) were placed on all boundaries. The materials used include a polymer with the refractive index, n p o l = 1.5 , while for dispersive materials namely, the gold metal, a frequency-dependent complex-valued relative permittivity based on the Lorentz-Drude model was adopted [28]:
ε r = ε + m = 1 M f m   ω p 2 ω m 2 ω 2 + i Γ m ω
where the sum over m corresponds to different oscillators, ε is the high-frequency contribution to the relative permittivity, ω p is the plasma frequency of free electrons, f m is the oscillator strength, ω m is the resonance frequency, and Γ m is the damping coefficient. We use M = 6 damped harmonic oscillators in order to take the inter-band transitions of gold. The values of the different parameters were taken from [29] and are listed in Table 1. Throughout the paper, unless stated otherwise, the geometrical parameters are taken as follows: W = 320   n m , L = 320   n m , d p o l = 300   n m , L p o l = 50   n m , L A u = 100   n m .
Figure 1. (a) Design of the 3D studied system composed of a gold nano-ridge with a thin polymer layer of thickness d p o l = 50 nm in the middle, placed above a gold layer with a polymer layer acting as a spacer in between. A zoom of the structure studied in 2D space is shown at the right side. (b) Simulated scattering cross section of the system in (a), as a function of wavelength λ and the angle of incidence θ . (c) Same as in (b), but for the fixed incidence angle, θ = 40 ° .. (d)-(f) Field maps of the scattered electric field within and around the nano-ridge corresponding to modes M1, M2, and M3 shown in (c) at wavelengths ~975 nm, ~1493 nm, and ~1689 nm, respectively.
Figure 1. (a) Design of the 3D studied system composed of a gold nano-ridge with a thin polymer layer of thickness d p o l = 50 nm in the middle, placed above a gold layer with a polymer layer acting as a spacer in between. A zoom of the structure studied in 2D space is shown at the right side. (b) Simulated scattering cross section of the system in (a), as a function of wavelength λ and the angle of incidence θ . (c) Same as in (b), but for the fixed incidence angle, θ = 40 ° .. (d)-(f) Field maps of the scattered electric field within and around the nano-ridge corresponding to modes M1, M2, and M3 shown in (c) at wavelengths ~975 nm, ~1493 nm, and ~1689 nm, respectively.
Preprints 115070 g001aPreprints 115070 g001b
To obtain the dispersion relation of the system, we have simulated the scattering cross-section σ s c as a function of both the angle of incidence in the range [ 0 ° 85 ° ], and the wavelength in the range of interest [ 600   n m 2000   n m ], that covers the visible and near-infrared domains. The corresponding results are shown in Figure 1b in color scale. Let us indicate that red color indicates maximum scattering intensity while blue indicates minimum (scattering) intensity. Three distinctive modes are observed, highlighted with horizontal white lines, as a guide for the eye. We refer to these modes as M1 for the lower wavelength (or higher frequency) mode, M2 for the middle wavelength mode and M3 for the higher wavelength mode. It is worth noting that modes M2 and M3 are excited for nearly all values of θ 0 ° ,   ~ 70 ° , while mode M1 is only excited in the range of θ between approximately ~ 20 ° and ~ 60 ° . To give an idea about the behavior of these modes in the scattering spectra, Figure 1c shows σ s c as a function of wavelength at = 40 ° , indicating three peaks corresponding to LSPs resonances at wavelengths 975 nm (mode M1), 1493 nm (mode M2), and 1689 nm (mode M3). Figure 1d–f display the maps of the scattered electric field at these wavelengths, respectively. The field localization shows the presence of two types of modes: one confined underneath the ridges (modes M1 and M3) and one at the middle of the ridge (mode M2). A common characteristic between the modes is the fact that they are trapped in between two metal materials within a dielectric, which suggests that they correspond to metal-insulator-metal (MIM) type modes. To explore in more details this aspect of the modes and their general behavior versus the key geometrical parameters, we analyzed the effect of the nano-ridge width in Figure 2a, where σ s c is plotted as a function of both the wavelength and the nano-ridge width, W . The remaining parameters are kept as in Figure 1, while the incidence angle is set to, θ = 40 ° , ensuring the excitation of the three modes. One notes that the resonance wavelengths of the three modes (corresponding to the scattering maximum/peaks) increase linearly with the increase of W. This trend is typical in an effective plasmonic cavity acting as a Fabry- Pérot cavity, where MIM-type resonances correspond to confined plasmonic standing waves [13,30]. Accordingly, the thin dielectric layers either underneath or at the middle of the nano-ridge, where the modes are localized, acts as an effective plasmonic Fabry-Pérot cavity (PFPC) of nominal length W .
In particular, the plasmonic-MIM standing waves build-up due to reflections off the PFPC edges on each side. The resonance wavelength of the (plasmonic) modes can be expressed as [30]:
λ r = 2 n e f f W m
where, n e f f = k s p p / k 0 is the effective refractive index of the of the propagating MIM-like plasmon in the cavity, k s p p is its wave vector, k 0 is the wavevector in the vacuum, and m = 1 , 2, 3, … is an integer representing the mode number. This equation is derived from the resonance condition corresponding to a round-trip of the wave in the cavity, namely, 2   k s p p W + 2 ϕ r e f = π m , where ϕ r e f is the phase-shift due to the reflection off the edges of the nano-ridge, whose value can be approximated as ϕ r e f k s p p δ W [30], where δ W is the effective length extension of the mode off the nano-ridge edges. Therefore, W is actually replaced by the length, W e f f = W + δ W . By fitting the LSP resonance wavelengths (shown with open white circles in Figure 2a), we estimated δ W to be approximately δ W 142   n m for mode M1, δ W 152   n m for mode M2, and δ W 48   n m for mode M3. As for the values of m   for each mode, these can be determined from the plot of the localized field inside the PFPC along the x axis. From Figure 1d–f, we can deduce that m = 2 for mode M1, while for modes M2 and M3 m = 1. The solid white lines in Figure 2a represents the fits of the numerical data (white open circles) using Eq. (3). We observe using the PFPC model based on the latter equation, that we have successfully fitted the numerical results (Figure 2a), which confirms the mechanism behind the LSPs modes, M1, M2 and M3 formation.
In order to better understand the physics behind these modes, we have examined in Figure 2b the effect of the polymer layer thickness d p o l in the middle of the nano-ridge, for θ = 40 ° . The rest of the parameters are set similar to Figure 1. Let us recall that the field is essentially confined inside the polymer layer for mode M2, whereas for modes M1 and M3 the field is localized underneath the ridge. This, means that only mode M2 should be sensitive to d p o l , which is evident in Figure 2b as only mode M2 is affected by changes in d p o l . In particular, one remarks that the resonance wavelength of mode M2, λ r 2 , drops nearly hyperbolically as d p o l increases. This trend reminds us that of the effective refractive index ( n e f f ) of a metal-insulator-metal type waveguide versus its width, whereby n e f f changes as the inverse of the guide width d i.e., n e f f 1 / d [31]. This relationship can be exploited to explain the behavior of λ r 2 versus d p o l considering what we mentioned above namely, mode M2 (or M1/M3) arises out of the trapping of a MIM-type mode in an effective plasmonic cavity (in the form of a MIM-standing wave). Hence, the effective refractive index of such MIM-plasmonic modes should changes as n e f f   1 / d p o l , which enables to estimate the resonance wavelength as:
λ r = 2 W m × C d p o l
where, C is a fitting parameter. Applying equation (3), we have fitted the numerically calculated (FEM) resonance wavelength of mode M2 in Figure 2b. Note that modes M1 or M3 are not fitted given their insensitivity to d p o l . The curve based on equation (3) is depicted in a white solid line while the FEM results are shown in white open circles; we can notice the very good agreement achieved. The fitting provided in Figure 2 using equations (2) and (3), enables us to definitely validate the physical mechanism behind the formation of modes M1, M2 and M3, mentioned above.

3. Localized Surface Plasmons of Coupled Gold Nano-ridges (Dimer)

In this section, we investigate the effect of two identical coupled nano-ridges that form a dimer. This configuration is expected to enhance the localization of the plasmonic field and the quality factors of the modes, leading to stronger light-matter interaction [13]. The two nano-ridges are set to be separated by a distance d m , as shown in the 2D (x-z) plane of the simulated system in Figure 3a. Note that d m is a key parameter as it enables to control the coupling magnitude in-between the nano-ridges (forming the dimer). Figure 3b, shows the simulated scattering cross section as a function of wavelength and incidence angle, with d m = 50   n m , while the other parameters are set as earlier (Figure 1). The plot reveals that for θ < 30 ° , the modes are indistinguishable and not excited, but rather correspond to large scattering bands. However, for θ > 30 ° , four modes labelled D1, D2, D3 and D4, show up. To examine these modes more closely, Figure 3c shows the scattering cross section at θ = 55° as a function of wavelength (blue line), where four peaks are clearly identified at roughly 816 nm (D1), 975 nm (D2), 1483 nm (D3), and 1690 nm (D4). A significant observation is that the quality factor of the modes in the dimer case is improved by approximately a factor of 3 compared to the single monomer case (red dashed line). For instance, if we focus on modes D3 and D4 which are the most strongly confined within the dimer, their quality factors are estimated as Q D 3 14.5 and Q D 4 18.2 , respectively. In comparison, mode M3 (of the monomer), which is the closest in terms of wavelength, has a quality factor of Q M 3 6.25 . This improvement is attributed to the coupling between the ridges, which enhances the localization of the electric field as it is evident from the field maps of the scattered field shown in Figure 4d–g. This enhanced localization is particularly evident for modes D2 and D4, where the field is strongly confined underneath and between the ridges, as well as inside the polymer layers in the middle of the ridges. For mode D2, the field is primarily localized in the polymer layers at the middle, but also in-between the ridges. In contrast, mode D1 exhibits weaker localization, having antenna-like components associated with multipolar localizations at the edges, that causes the field to radiate. Let us mention that the enhanced quality factor of the dimer modes comparing with the monomer case has been reported elsewhere [13,25]. As a matter of comparison, we have simulated the scattering cross section as a function of wavelength for θ=55°, for the dimer-structure without the polymer layer inserted in the middle of the ridges (blue dashed line). The curve is overlaid with that of a dimer (blue solid line) and monomer (red dashed line) structures with the polymer layer set in the middle as shown in Figure 3c. We can observe that the insertion of the thin polymer layers within the dimer ridges results in improved quality factor of the modes. In addition, the spectrum is quite richer in the latter case as additional localized surface plasmons (with four modes) show up as compared to the case of no layers in the ridges (three modes). Furthermore, the localization of the field is significantly better with the polymer layer as seen in Figure 3d–g, compared to the field maps in the case where no layer is inserted (Figure 3h–j).
To analyze the impact of geometrical parameters on the dimer modes and the physical origin of the latters, we examine how variations in the width of the dimer L d i m e r , and the thickness of the polymer layer in the middle of the ridges, d p o l , affect the modes. The dimer width is defined as, L d i m e r = 2   W + d m ,   where W represents the single ridge width and d m is the separation distance between the ridges. For the sake of simplicity, we alter L d i m e r by changing only W and keeping d m , that controls coupling intensity within the dimer, fixed. Also, all remaining parameters of the system are set equivalent to Figure 1a (as mentioned above). Figure 4a, presents the scattering cross section as a function of wavelength and L d i m e r . We note a near linear redshift in the resonance wavelengths as L d i m e r is increased, which is similar to what we observed above in case of the monomer (Figure 2a) width increasing effect. This hints at the fact that modes, D1-D4, of the dimer could be considered as plasmonic MIM-type standing waves confined in an effective plasmonic cavity. The latter is either made of the thin layers in the middle of the ridges plus the small space in between (like for mode D3) or of the gap layer underneath the dimer, where mode D4 for instance is primarily localized. Let us call the latter cavity as the lower cavity, and the former as the upper cavity. Note that some modes (D1, D2 and D4) are confined in both the lower and upper cavities (with a different rate) as well as in-between the dimer. To confirm this, we consider that the MIM-type plasmonic standing waves corresponding to the modes are localized in either the lower or upper cavities (or both) of nominal length, L d i m e r , and hence, the resonance wavelength is given by:
λ r = 2 n e f f   L d i m e r m
where all the parameters are defined in an equivalent manner as in Eq. (3) for the monomer. It should be stressed that the nominal length is replaced by, L e f f = L d i m e r + δ L d i m e r , where the latter parameter is a fitting parameter. Such a parameter is meant to account for the mode extension off the dimer outer edges as already explained.
As for the modes order, m, it is determined based on the localized field in the effective plasmonic cavity (lower or upper). From Figure 3d–g, we estimate m = 4 for the mode D2, m = 2 for both modes D3 and D4. We should mention that because mode D1 has complex features like antenna-like components, with no neat/simple localization features either in the lower or upper cavity, we do not consider fitting its resonance wavelength. Figure 4b shows the fitting provided by Eq. (4) (solid lines) to the position of the modes D2, D3 and D4, resonance wavelengths (indicated by symbols) as a function of L d i m e r . The fitting parameter δ L d i m e r has been determined for each mode as: 58.35 nm for mode D2, 11 nm for mode D3, and 286.76 nm for mode D4. We note that the fitting agrees perfectly with the FEM results for all the three modes which clearly underscores the MIM-character of such modes, and that they correspond to MIM-type plasmon standing waves trapped in the effective lower/upper cavities.
We also explored the effect of varying the thickness of the polymer layers inserted in the middle of the ridges d p o l . Figure 4c shows the scattering cross section as a function of wavelength and d p o l . This plot indicates that all modes are affected by changes in d p o l , but the impact is significantly more pronounced for mode D3, which is understandable since mode D3 is primarily localized within such layers. Following mode D3, we see that mode D2 is the most affected as compared to mode D4, which is consistent with their field localization as the latter is more confined beneath the ridges than in their middle while for the former it is the opposite. Note that for mode D1 due to its complex features, it is difficult to analyze its behaviour versus changes in d p o l , however, it can be remarked that for d p o l 50   n m , the mode is nearly cut-off as its scattering section becomes quite weak. For d p o l 50   n m , nevertheless, the mode shows a sensitivity with the increase of d p o l as its wavelength drops with a non-negligeable magnitude, reflecting its localization in the (middle) layers. Using a similar equation to Eq. (4) and the same reasoning as in the case of the monomer, one can capture the modes, D3, D2 and D4 behaviour as d p o l is changed. Equation (4) can be written in case of the dimer as, λ r = 2 L e f f m × C d p o l , where the ridge width W is simply replaced by the dimer effective width, L e f f , and C is a fitting parameter. Based on such equation, we provide the fitting of the three modes resonance wavelengths as a function of d p o l in Figure 4d, where we can observe the excellent agreement with the FEM results. Hence, we can be quite confident that modes D3, D2 and D4 arise out of the trapping of MIM-type surface plasmon guided mode in the lower/upper effective cavity, giving rise to resonating MIM plasmonic standing waves.

4. Effective Volume of the Dimer Localized Plasmon Modes and Prospective Optomechanical Coupling

Here, we focus on the analysis of the effective mode volume of the localized surface plasmon modes of the dimer, aiming to classify these modes for an optimal enhanced optomechanical coupling with localized phonons supported by the dimer. As a matter of fact, the present study which is about a thorough analysis of the plasmonic properties of the dimer-based device, comes after a previous study about phononic characteristics of the system [27]. The general realm being to achieve enhanced plasmon-phonon coupling using surface acoustic waves (Figure 5a). In that framework, we would like to provide in this last paragraph a prospective qualitative analysis of enhanced optomechanical coupling harnessing the dimer plasmonic modes effective volume. Note that the optical mode volume is one of the key parameters that govern light-matter interactions in general, as it is related to field enhancement and spatial localization [32,33]. Having an estimation of the effective mode volume can help quantify the spatial extent of the electromagnetic energy confinement, directly impacting the interaction strength between light and matter. To get an estimate of the effective mode volumes V m , we use the quasi-static approximation to calculate the mode-volume integral for a leaky plasmonic cavity given [34,35,36]:
V m = W   d r 3 max W
where, the integral is performed over of the volume (surface times unit length) of both nano-ridges, W is the electromagnetic energy density and max W is its maximum value in the simulation domain. For a dielectric medium, the energy density takes the simple form, W = ε E 2 , while in the case of a dispersive and lossy medium like a metal, an appropriate relation is as follows [31]:
W = ε 1 + 2 ω ε 2 / Γ E 2
where, ε ω = ε 1 ω + i ε 2 ω is the complex dielectric function of the metal. In our case, the latter is given by the Lorentz-Drude model (Eq. 1). Γ on the other hand, is the mode damping taken as the linewidth of the respective LSP mode. As already mentioned, the simulation is performed in 2D due to symmetry along y-axis, so to calculate the mode volume one should perform surface integrals in the 2D-cross section system, and then multiply the result by a unit length along the y-direction. Figure 5b, shows the normalized effective mode volume V m / V 0 calculated for each mode of the dimer (D1, D2, D3 and D4), where V 0 = λ / n p o l 2 ( × 1   m ), is the volume occupied by the associated resonance wavelength of the mode in the polymer layer. As can be seen, the effective volumes of modes are well below λ / 2 n p o l , reflecting the strong confinement of the dimer plasmonic modes far below diffraction limit of light. From the values of the mode volumes, we observe that mode D1 has the largest value of 0.0117   λ / n p o l 2 , indicating relatively lower field localization rate compared to the other modes. This relates to such mode localization features like its tendency to radiate energy, due to its antenna-like components associated with multipolar localizations at the edges. In contrast, mode D4 with the lowest value of 0.0016   λ / n p o l 2 , is relatively the most strongly confined mode. Mode D2 and D3 fall in-between with values 0.00306   λ / n p o l 2 and 0.00704   λ / n p o l 2 respectively. This result shows that for enhanced light-matter coupling purposes, mode D4 should offer the highest interaction rate as the coherent coupling rate between light and matter scales as g 1 / V m [36].
Regarding the optomechanical coupling specifically, we can affirm that the plasmon-phonon interaction is qualitatively related to the overlap degree between the plasmon and the phonon in a given cavity. This is estimated quantitively via the so-called optomechanical coefficient whose expression is determined based on the photoelastic volume and the surface moving-boundary effects [13,37]. Given that the dimer plasmonic modes with the lowest volume magnitude V m / V 0 are the most highly trapped plasmons, they should manifest the higher overlap rate with a given phonon localized in the dimer. As a result, such modes should lead to an optimal optomechanical coupling intensity. In Figure 5c–h, we recall the localized phonons of the dimer that we deemed to be optomechanically active based on symmetry arguments in a previous work [27], and the design enabling to use surface acoustic waves to excite such phonons. Considering for instance the in-phase compression mode A2 (Figure 5d), that deforms the polymer layers both under and in the middle of the ridges (see Figure 5g), where modes D2, D3 and D4 are localized, we can roughly classify such modes (based on their volume) in terms of optomechanical coupling rate g as:
g D 4 g D 2 g D 3 .
The same reasoning approximately applies if we consider another phonon of the dimer, to classify the LSPs modes in terms of acousto-plasmonic coupling magnitude. Obviously, a thorough study based on the rigorous simulations of the optomechanical coefficient g is needed to determine the best localized phonons and plasmons that yield the highest coupling magnitude. However, this qualitative argument based on the modes volume gives already a good hint regarding the most promising plasmonic mode for the enhanced optomechanical purpose namely, mode D4 that shows the lowest volume or equivalently the highest confinement rate.

5. Conclusion

In this paper, thorough finite element-based numerical calculations have been carried-out to investigate localized surface plasmons resonances in a film-coupled nanoparticles system. The proposed design is made out of two gold nano-ridges forming a dimer atop a metal plate in-between which a thin polymer film acting as spacer is inserted. We have first studied optical scattering properties of a monomer ridge, and found three MIM-type localized surface plasmons in the frequency range of interest namely, the visible and near infrared domains. Our results indicate that the presence of the polymer layer in the middle of the ridge gives rise to additional LSPs modes confined within. In the dimer structure composed of two identical nanoridges, we observed that the plasmonic modes exhibit narrower linewidths and enhanced field localization compared to the monomer configuration. The enhanced localization is attributed to the strong coupling between the ridges leading to a higher rate of modes trapping. We have also shown that the inserted polymer layers within the ridges, in their middle, further enhance the modes quality factor as compared to the ridges with no layers. The former case, adds means to confine even more the dimer LSPs modes (in the middle layers), which results in narrower linewidths. Also, in that case a richer plasmonic spectrum based on the dimer is found as additional modes trapped in the ridges middle layers show up. We have successfully demonstrated that the dimer/monomer modes arise out of the trapping of MIM-type guided modes in the polymer layers beneath or in the middle of the ridges acting as plasmonic effective cavities.
Furthermore, based on the simulation of the LSPs modes effective volume of the dimer, we have been able to roughly class the modes in terms of their optomechanical coupling magnitude, with a given phonon of the dimer. Specifically, mode D4 with its smallest normalized mode volume, stands out as the most promising LSP mode for enhanced plasmon-phonon interaction. The showcased dimer-device in this work, could help design new acousto-optic components for high frequency nanomechanical sensing, optical switching or to perform infrared/visible light modulation operations.

Supplementary Materials

No supplementary materials are available for this article.

Author Contributions

Conceptualization, A.N., E.H.B. and M.E.G.; methodology, M.E.G, M.A. and M.M.; validation, A.N., E.H.B and M.E.G; investigation, M.E.G., M.A. and M.M.; data curation, M.E.G. and A.N.; writing—original draft preparation, M.E.G. and A.N.; writing—review and editing, A.N., E.H.B and M.E.G.; funding acquisition, E.H.B. and A.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

We would like to acknowledge the financial support under the framework of the Research Excellence Grants Program from the Moroccan National Center for Scientific and Technical Research (CNRST).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Stockman, M.I.; Kneipp, K.; Bozhevolnyi, S.I.; Saha, S.; Dutta, A.; Ndukaife, J.; Kinsey, N.; Reddy, H.; Guler, U.; Shalaev, V.M. Roadmap on plasmonics. J. Opt. 2018, 20, 043001. [Google Scholar] [CrossRef]
  2. Baumberg, J.J.; Aizpurua, J.; Mikkelsen, M.H.; Smith, D.R. Extreme nanophotonics from ultrathin metallic gaps. Nat. Mater. 2019, 18, 668–678. [Google Scholar] [CrossRef] [PubMed]
  3. Amrani, M.; Khattou, S.; Rezzouk, Y.; Mouadili, A.; Noual, A.; El Boudouti, E.H.; Djafari-Rouhani, B. Analytical and numerical study of T-shaped plasmonic demultiplexer based on Fano and induced transparency resonances. J. Phys. D Appl. Phys. 2021, 55, 075106. [Google Scholar] [CrossRef]
  4. Duan, Q.; Liu, Y.; Chang, S.; Chen, H.; Chen, J.H. Surface plasmonic sensors: Sensing mechanism and recent applications. Sensors 2021, 21, 5262. [Google Scholar] [CrossRef] [PubMed]
  5. Noual, A.; Akiki, R.; Pennec, Y.; El Boudouti, E.H.; Djafari-Rouhani, B. Surface acoustic waves-localized plasmon interaction in pillared phononic crystals. Phys. Rev. Appl. 2020, 13, 024077. [Google Scholar] [CrossRef]
  6. Koya, A.N.; Cunha, J.; Guerrero--Becerra, K.A.; Garoli, D.; Wang, T.; Juodkazis, S.; Proietti Zaccaria, R. Plasmomechanical systems: Principles and applications. Adv. Funct. Mater. 2021, 31, 2103706. [Google Scholar] [CrossRef]
  7. Barzanjeh, S.; Xuereb, A.; Gröblacher, S.; Paternostro, M.; Regal, C.A.; Weig, E.M. Optomechanics for quantum technologies. Nat. Phys. 2022, 18, 15–24. [Google Scholar] [CrossRef]
  8. Liu, X.; Liu, W.; Ren, Z.; Ma, Y.; Dong, B.; Zhou, G.; Lee, C. Progress of optomechanical micro/nano sensors: A review. Int. J. Optomechatronics 2021, 15, 120–159. [Google Scholar] [CrossRef]
  9. Sansa, M.; Defoort, M.; Brenac, A.; Hermouet, M.; Banniard, L.; Fafin, A.; Gely, M.; Masselon, C.; Favero, I.; Jourdan, G. Optomechanical mass spectrometry. Nat. Commun. 2020, 11, 3781. [Google Scholar] [CrossRef]
  10. Sansa, M.; Defoort, M.; Brenac, A.; Hermouet, M.; Banniard, L.; Fafin, A.; Gely, M.; Masselon, C.; Favero, I.; Jourdan, G. Optomechanical mass spectrometry. Nat. Commun. 2020, 11, 3781. [Google Scholar] [CrossRef]
  11. Ramos, D.; Malvar, O.; Davis, Z.J.; Tamayo, J.; Calleja, M. Nanomechanical plasmon spectroscopy of single gold nanoparticles. Nano Lett. 2018, 18, 7165–7170. [Google Scholar] [CrossRef] [PubMed]
  12. Simone, G. Strong plasmon-mechanical coupling through standing acoustic waves and antenna enhancement, toward mass transduction. J. Appl. Phys. 2024, 135. [Google Scholar] [CrossRef]
  13. Noual, A.; Akiki, R.; Lévêque, G.; Pennec, Y.; Djafari-Rouhani, B. Enhanced phonon-plasmon interaction in film-coupled dimer nanoridges mediated by surface acoustic waves. Opt. Express 2021, 29, 43104–43123. [Google Scholar] [CrossRef]
  14. Ahmidayi, N.; d’Orsonnens, W.; Maurer, T.; Lévêque, G. Mechanical Enhancement of the Strain--Sensor Response in Dimers of Strongly Coupled Plasmonic Nanoparticles. Ann. Der Phys. 2023, 535, 2300319. [Google Scholar] [CrossRef]
  15. Ahmidayi, N.; d’Orsonnens, W.; Maurer, T.; Lévêque, G. Mechanical Enhancement of the Strain--Sensor Response in Dimers of Strongly Coupled Plasmonic Nanoparticles. Ann. Der Phys. 2023, 535, 2300319. [Google Scholar] [CrossRef]
  16. Maurer, T.; Adam, P.-M.; Lévêque, G. Coupling between plasmonic films and nanostructures: From basics to applications. Nanophotonics 2015, 4, 363–382. [Google Scholar] [CrossRef]
  17. Goncalves, M.R. Plasmonic nanoparticles: Fabrication, simulation and experiments. J. Phys. D Appl. Phys. 2014, 47, 213001. [Google Scholar] [CrossRef]
  18. Charchi, N.; Li, Y.; Huber, M.; Kwizera, E.A.; Huang, X.; Argyropoulos, C.; Hoang, T. Small mode volume plasmonic film-coupled nanostar resonators. Nanoscale Adv. 2020, 2, 2397–2403. [Google Scholar] [CrossRef] [PubMed]
  19. Huang, S.; Ming, T.; Lin, Y.; Ling, X.; Ruan, Q.; Palacios, T.; Wang, J.; Dresselhaus, M.; Kong, J. Ultrasmall mode volumes in plasmonic cavities of nanoparticle--on--mirror structures. Small 2016, 12, 5190–5199. [Google Scholar] [CrossRef]
  20. Deacon, W.M.; Lombardi, A.; Benz, F.; del Valle-Inclan Redondo, Y.; Chikkaraddy, R.; de Nijs, B.; Kleemann, M.-E.; Mertens, J.; Baumberg, J.J. Interrogating nanojunctions using ultraconfined acoustoplasmonic coupling. Phys. Rev. Lett. 2017, 119, 023901. [Google Scholar] [CrossRef]
  21. Yeshchenko, O.A.; Kozachenko, V.V.; Naumenko, A.P.; Berezovska, N.I.; Kutsevol, N.V.; Chumachenko, V.A.; Haftel, M.; Pinchuk, A.O. Gold nanoparticle plasmon resonance in near-field coupled Au NPs layer/Al film nanostructure: Dependence on metal film thickness. Photonics Nanostructures-Fundam. Appl. 2018, 29, 1–7. [Google Scholar] [CrossRef]
  22. Armstrong, R.E.; Van Liempt, J.C.; Zijlstra, P. Effect of film thickness on the far-and near-field optical response of nanoparticle-on-film systems. J. Phys. Chem. C 2019, 123, 25801–25808. [Google Scholar] [CrossRef]
  23. Mock, J.J.; Hill, R.T.; Degiron, A.; Zauscher, S.; Chilkoti, A.; Smith, D.R. Distance-dependent plasmon resonant coupling between a gold nanoparticle and gold film. Nano Lett. 2008, 8, 2245–2252. [Google Scholar] [CrossRef]
  24. Lee, H.-Y.; Seong, H.K.; Yoo, J. Optimal shape design of the film-coupled nanoparticle using the phase field design method. Int. J. Precis. Eng. Manuf. 2016, 17, 453–460. [Google Scholar] [CrossRef]
  25. Li, G.-C.; Zhang, Y.-L.; Lei, D.Y. Hybrid plasmonic gap modes in metal film-coupled dimers and their physical origins revealed by polarization resolved dark field spectroscopy. Nanoscale 2016, 8, 7119–7126. [Google Scholar] [CrossRef]
  26. Kunert, B.; Guo, W.; Mols, Y.; Tian, B.; Wang, Z.; Shi, Y.; Van Thourhout, D.; Pantouvaki, M.; Van Campenhout, J.; Langer, R. III/V nano ridge structures for optical applications on patterned 300 mm silicon substrate. Appl. Phys. Lett. 2016, 109. [Google Scholar] [CrossRef]
  27. El Ghafiani, M.; Noual, A.; Rezzouk, Y.; Moutaouekkil, M.; El Boudouti, E.H. Surface acoustic waves interaction with confined acoustic phonons in a coupled nanoridges dimer set atop of a multilayer design. Arab J. Basic Appl. Sci. 2023, 30, 593–600. [Google Scholar] [CrossRef]
  28. Raman, A.; Fan, S. Photonic Band Structure of Dispersive Metamaterials Formulated as a Hermitian Eigenvalue Problem. Phys. Rev. Lett. 2010, 104, 087401. [Google Scholar] [CrossRef]
  29. Palik, E.D. Handbook of optical constants of solids; Academic press, 1998; Volume 3. [Google Scholar]
  30. Dorfmüller, J.; Vogelgesang, R.; Weitz, R.T.; Rockstuhl, C.; Etrich, C.; Pertsch, T.; Lederer, F.; Kern, K. Fabry-Pérot resonances in one-dimensional plasmonic nanostructures. Nano Lett. 2009, 9, 2372–2377. [Google Scholar] [CrossRef] [PubMed]
  31. Maier, S.A. Plasmonics: Fundamentals and applications; Springer, 2007; Volume 1. [Google Scholar]
  32. Maier, S.A. Plasmonic field enhancement and SERS in the effective mode volume picture. Opt. Express 2006, 14, 1957–1964. [Google Scholar] [CrossRef]
  33. Baumberg, J.J.; Aizpurua, J.; Mikkelsen, M.H.; Smith, D.R. Extreme nanophotonics from ultrathin metallic gaps. Nat. Mater. 2019, 18, 668–678. [Google Scholar] [CrossRef] [PubMed]
  34. Koenderink, A.F. On the use of Purcell factors for plasmon antennas. Opt. Lett. 2010, 35, 4208–4210. [Google Scholar] [CrossRef] [PubMed]
  35. Derom, S.; Vincent, R.; Bouhelier, A.; des Francs, G.C. Resonance quality, radiative/ohmic losses and modal volume of Mie plasmons. Europhys. Lett. 2012, 98, 47008. [Google Scholar] [CrossRef]
  36. Srinivasan, K.; Borselli, M.; Painter, O.; Stintz, A.; Krishna, S. Cavity Q, mode volume, and lasing threshold in small diameter AlGaAs microdisks with embedded quantum dots. Opt. Express 2006, 14, 1094–1105. [Google Scholar] [CrossRef]
  37. Ahmed, A.; Pelton, M.; Guest, J.R. Understanding how acoustic vibrations modulate the optical response of plasmonic metal nanoparticles. ACS Nano 2017, 11, 9360–9369. [Google Scholar] [CrossRef]
Figure 2. (a) Scattering cross section as a function of wavelength and the width of the nano-ridge (W). White open circles indicate the position of the LSP resonances for various values of W, while white solid lines represent a fit of these points using MIM-like waveguide model (Eq. (2)). (b) Similar to Figure 2a, though this time, the scattering cross section is plotted as a function of wavelength and the thickness of the polymer layer at the middle of nano-ridge showing the effect of varying the thickness of the dielectric layer at the middle of the ridge (Eq. (3)).
Figure 2. (a) Scattering cross section as a function of wavelength and the width of the nano-ridge (W). White open circles indicate the position of the LSP resonances for various values of W, while white solid lines represent a fit of these points using MIM-like waveguide model (Eq. (2)). (b) Similar to Figure 2a, though this time, the scattering cross section is plotted as a function of wavelength and the thickness of the polymer layer at the middle of nano-ridge showing the effect of varying the thickness of the dielectric layer at the middle of the ridge (Eq. (3)).
Preprints 115070 g002
Figure 3. (a) 2D schematic representation of the simulated system consisting of two nano-ridges sitting atop a multilayer structure, separated by a distance dm. (b) Scattering cross section of the system in (a), as a function of wavelength λ and the incidence angle ϴ. (c) Same as in (b), but at θ = 55 ° , for the monomer structure (red dashed lines), the dimer structure with the polymer layer set at the middle of the ridges (blue solid line) and the dimer structure without the polymer layer at the middle (blue dashed line). (d)-(g) Scattered field map of the modes D1, D2, D3, and D4 in the case of the nano-ridges with a thin polymer layer placed at the middle at wavelengths, λ D 1 ~   816   n m , λ D 2 ~   975   n m , λ D 3 ~   1483   n m , and λ D 4 ~   1690   n m , respectively. (h)-(j) Same as in (d), but for the modes A1, A2 and A3 of the nano-ridges without the polymer layer at wavelengths, λ A 1 ~   735   n m , λ A 2 ~   1005   n m , and λ A 3 ~   1740   n m , respectively.
Figure 3. (a) 2D schematic representation of the simulated system consisting of two nano-ridges sitting atop a multilayer structure, separated by a distance dm. (b) Scattering cross section of the system in (a), as a function of wavelength λ and the incidence angle ϴ. (c) Same as in (b), but at θ = 55 ° , for the monomer structure (red dashed lines), the dimer structure with the polymer layer set at the middle of the ridges (blue solid line) and the dimer structure without the polymer layer at the middle (blue dashed line). (d)-(g) Scattered field map of the modes D1, D2, D3, and D4 in the case of the nano-ridges with a thin polymer layer placed at the middle at wavelengths, λ D 1 ~   816   n m , λ D 2 ~   975   n m , λ D 3 ~   1483   n m , and λ D 4 ~   1690   n m , respectively. (h)-(j) Same as in (d), but for the modes A1, A2 and A3 of the nano-ridges without the polymer layer at wavelengths, λ A 1 ~   735   n m , λ A 2 ~   1005   n m , and λ A 3 ~   1740   n m , respectively.
Preprints 115070 g003aPreprints 115070 g003b
Figure 4. (a) Scattering cross section of the dimer-based system (Figure 3a) as a function of the wavelength and L d i m e r . (b) Plot of the resonance wavelengths of modes D2, D3 and D4 as a function of L d i m e r (open symbols). The plot for each mode is fitted (solid lines) with Eq. (5) giving the resonance wavelength of the effective (lower or upper) plasmonic cavity. Note that mode D3 is confined in the upper cavity while D4 is mostly trapped in the lower cavity, and D2 is localized in the upper cavity and in the lower to a lesser extent. (c) Similar to (a), but as a function of the polymer layers thickness, d p o l . (d) Same as in (b), but as a function of d p o l , and using Eq. (6) as for the provided fitting to the modes.
Figure 4. (a) Scattering cross section of the dimer-based system (Figure 3a) as a function of the wavelength and L d i m e r . (b) Plot of the resonance wavelengths of modes D2, D3 and D4 as a function of L d i m e r (open symbols). The plot for each mode is fitted (solid lines) with Eq. (5) giving the resonance wavelength of the effective (lower or upper) plasmonic cavity. Note that mode D3 is confined in the upper cavity while D4 is mostly trapped in the lower cavity, and D2 is localized in the upper cavity and in the lower to a lesser extent. (c) Similar to (a), but as a function of the polymer layers thickness, d p o l . (d) Same as in (b), but as a function of d p o l , and using Eq. (6) as for the provided fitting to the modes.
Preprints 115070 g004
Figure 5. (a) Effective mode volume of the dimer LSPs modes normalized to V 0 = λ / n p o l 2 × 1   m , the volume occupied by the associated resonance wavelengths in the polymer layer. (b) Design of the system used to excite confined acoustic phonons by means of surface acoustic waves (SAWs). The device consists in the dimer structure of Figure 3a deposited above an SiO2 substrate. (c)-(e) Zoom in on the mechanical displacement field in and around the dimer at its localized phonons resonance frequencies namely at, f A 1 = 0.238   G H z for the symmetric flexural mode, f A 2 = 0.456   G H z for the in-phase compression mode and f A 3 = 0.619   G H z for the high order symmetric flexural mode. (f)-(h) Plot of the displacement field in the polymer layers, both underneath and in the middle of the ridges, associated with the phonon modes in (c)-(e), respectively.
Figure 5. (a) Effective mode volume of the dimer LSPs modes normalized to V 0 = λ / n p o l 2 × 1   m , the volume occupied by the associated resonance wavelengths in the polymer layer. (b) Design of the system used to excite confined acoustic phonons by means of surface acoustic waves (SAWs). The device consists in the dimer structure of Figure 3a deposited above an SiO2 substrate. (c)-(e) Zoom in on the mechanical displacement field in and around the dimer at its localized phonons resonance frequencies namely at, f A 1 = 0.238   G H z for the symmetric flexural mode, f A 2 = 0.456   G H z for the in-phase compression mode and f A 3 = 0.619   G H z for the high order symmetric flexural mode. (f)-(h) Plot of the displacement field in the polymer layers, both underneath and in the middle of the ridges, associated with the phonon modes in (c)-(e), respectively.
Preprints 115070 g005
Table 1. Lorentz-Drude parameters values adopted for modeling the relative permittivity of gold.
Table 1. Lorentz-Drude parameters values adopted for modeling the relative permittivity of gold.
Oscillator m f m ω p   r a d / s ω m r a d / s Γ m r a d / s
1 0.76 1.37188 × 10 16 0 8.05202 × 10 13
2 0.024 1.37188 × 10 16 6.30488 × 10 14 3.66139 × 10 14
3 0.01 1.37188 × 10 16 1.26098 × 10 15 5.24141 × 10 14
4 0.71 1.37188 × 10 16 4.51065 × 10 15 1.32175 × 10 15
5 0.601 1.37188 × 10 16 6.53885 × 10 15 3.78901 × 10 15
6 4.384 1.37188 × 10 16 2.02364 × 10 16 3.36362 × 10 15
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated