Abstract
We present a comprehensive study using finite element numerical analysis of the acoustic localized phonons supported by a gold nanoridges dimer-based multilayer design. The latter consists in a SiO2-substrate over which a gold film covered with a thin polymer is deposited. We investigate first the mechanical eigen-modes analysis of a single monomer ridge, where we find flexural and compressional type modes in the sub-GHz frequency range. This is realized by either setting the ridge in a periodic structure, which enables to get the dispersion curve of the modes, or by considering an isolated system bounded by perfect matched layers, where we use the equivalent of the local density of states to track the modes. A good agreement is obtained between the two methods. Similarly, we find in case of the coupled dimer ridges hybridized modes, namely the flexural and compressional modes of a monomer split-up into in- and -out-of-phase type modes. We demonstrate efficient coupling between the monomer/dimer localized phonons with surface acoustic waves (SAWs) as the simulated transmission spectra show dips at the frequencies of the monomer/dimer eigenmodes. For symmetry reasons, some of the dimer modes are expected to be optomechanically active. The proposed SAW-based device is meant to help design acousto-optic modulators or ultrasensitive sensors.
1. Introduction
Control and modulation of light properties in nanoscale photonic structures based on the design of acousto-optic nanodevices has concentrated a lot of attention lately (Laude et al., Citation2017; Gao et al., Citation2023). Such components can be useful in quantum information processing as well as to significantly enhance light-matter interaction at the nanoscale (Sansa et al., Citation2020; Barzanjeh et al., Citation2022). Physically speaking, the acousto-optical coupling (AOC) is described by two mechanisms namely, i) the photoelastic and ii) the moving boundary effects. In the latter, the boundaries deformations between two adjacent materials induce a local variation of the dielectric constant all along the interface (Royer & Dieulesaint, Citation2000) while in the former effect, the dielectric is modulated locally within the volume by the acoustic strain field (Balram, Davanço, Lim, Song, & Srinivasan, Citation2014). In order to reach a meaningful AOC magnitude, an appreciable overlap between the phonon field and the optical mode is essential, and this is often realized by colocalizing the latter in the same volume such as in an optomechanical cavity (Arregui et al., Citation2023). In particular, a very promising new kind of OM cavities has been proposed based on the so-called phoxonic crystals that exhibit band gaps for both photons and phonons, giving rise to highly localized light and sound waves, and subsequently to enhanced AOC rates (Aram & Khorasani, Citation2018).
Lately, nanoscale multilayered structured designs have attracted a great deal of interest due to their exquisite properties in terms of highly localized acoustic and optical modes they exhibit (Dobrzynski et al., Citation2017; Noual, Akiki, Pennec, El Boudouti, & Djafari-Rouhani, Citation2020). A special class of these structures is often referred to as film-coupled nanoparticles or particles on mirror due to the coupling of metallic nanoparticles to their image within a metal-film underneath (Li, Zhang, Maier, & Lei, Citation2018). One key feature of such designs is the strongly localized surface plasmons resonances modes they support, which have been harnessed in various applications (Saada et al, Citation2020; Oo, Silva et al, Citation2016; Cheng & Xu, Citation2021). In particular, localized surface plasmons have been shown useful in order to enhance inelastic Brillouin scattering (Noual, Kang, et al., Citation2021; Vasileiadis et al., Citation2022). Besides, the use of surface acoustic waves (SAWs) with high amplitudes (), helps to increase vibration magnitudes of localized phonons in metallic nanoparticles, which enhances the AOC rate meaningfully (Lin, Lin, & Hsu, Citation2015; Noual, Akiki, Lévêque, Pennec, & Djafari-Rouhani, Citation2021).
To design a SAW-device for enhanced AOC, we thoroughly investigate in this work localized acoustic phonons in a gold dimer ridges-based component using SAW. Specifically, numerical simulations using Comsol-Multiphysics are performed, so to find localized mechanical eigen-modes supported by the nanostructure, with a particular focus on the phonons considered potentially useful for AOC. The proposed device consists in SiO2 substrate over which a gold film covered with a thin polymer is deposited. Atop the structure surface, a gold nano-ridge dimer is placed, where each ridge contains a thin polymer layer lying in the middle. This feature can be used for the excitation of localized plasmons within the polymer layer that could be strongly modulated by localized phonons in the ridge. Prior studies have only considered localized plasmons mostly under the particles, meaning that the current design offers a supplemental degree of freedom. In order to understand the physical origin of the dimer-based system eigenmodes, we analyze the modes of a single monomer first. We show efficient excitation of compressional-like and flexural type eigen-modes of the dimer using SAWs in the sub-GHz range, where the coupled ridges hybridize their modes, yielding in- and -out-of-phase type modes. The symmetry of the modes indicate that compressional modes and anti-phase flexural modes should be optomechanically active (Noual, Akiki, et al., Citation2021). The reported findings in this paper could help design SAW-based acousto-optic modulators or ultrasensitive nanosensors, due to the easily accessible fabrication technology of SAW-devices nowadays (Mei, Zhang, & Friend, Citation2020; Kelly, Northfield, Rashid, Bao, & Berini, Citation2022).
2. Single monomer ridge eigenmodes
In this section, we analyse the mechanical eigenmodes of a single monomer ridge. We shall do this using two complementary approaches, namely, by studying first a periodic structure and then by performing a thorough eigen-analysis of the isolated ridge.
2.1. Periodic structure
We consider a periodic structure in which the single ridge atop the multilayer is repeated periodically along the direction as depicted in . More precisely, the structure design is made of a gold ridge of height h and width w, with a thin layer of polyamide of height hpol placed at the middle. The ridge monomer is deposited over a thin polyamide layer of thickness epol which coats a gold film measuring eAu. We chose a SiO2 semi-infinite substrate to support the multilayer (). It is worth indicating that the method followed here, enables us to compute the dispersion curve of the system which gives a global insight regarding the modes supported by the structure, including the localized type modes and their properties in terms of dispersion and symmetry. For that matter, we perform an eigenfrequency study of the system as the wave vector spans the first Brillouin zone, where is the period. Numerically, the periodicity is set using Bloch boundary conditions, and the eigenvalue equation is solved using structural mechanics module of Comsol 5.2a. At the bottom of the unit-cell of the lattice we use perfectly matched layer (PML) which absorb incident waves without reflections, in order to mimic a semi-infinite (SiO2) substrate. The elastic properties of the materials are defined based on young’s modulus E, Poisson ratio and the density as follows: Poison ratio and the density for the polyamide, and for the glass substrate SiO2 (Royer & Dieulesaint, Citation2000). As for the gold layer which is an anisotropic material with face centered cubic crystal, its elastic constants were taken from the same reference. shows the simulated dispersion curve for and while was taken large enough to mimic a semi-infinite substrate. As mentioned above, we focus our analysis on the modes at the edge of the first Brillouin zone, i.e. situated under the SiO2-sound line labelled as A1, A2, A3, and A4, which do not propagate within the substrate. The field maps of the displacement field norm associated with such modes near the monomer are given in . According to the field maps, the four branches modes (A1, A2, A3, and A4) are localized within the ridge with different rates, meaning that these are rather eigenmodes of the latter. As a matter of fact, the first branch A1 is a first order flexural mode at the mode A2 at is a compressional mode, whilst modes (A3, A4) at and correspond to the second order flexural mode which splits into two branches. This, happens due to its coupling with the gold/polymer films vibrations possibly induced by the periodicity effect, namely, the interaction with neighboring cells.
Indeed, the ridge-films coupling is apparent in the field maps () of modes A3 and A4. Note in particular the strong Sezawa type mode in the films interacting with the flexural mode A4 of the ridge, where the two modes vibrate in-phase. Conversely, for mode A3 the films and the ridge seem to vibrate out-of-phase, meaning that they tend to cancel-out. This might explain why the branch A3 is less dispersive comparing to A4, where the film and the ridge add-up their vibrations. In other words, there is more coupling with neighboring cells in this latter mode, hence the dispersion. The branch A1 is pretty much flat, due to a strong localization of such mode in the ridge as compared to the rest of the modes. In such mode, the ridge is clearly decoupled from its surrounding i.e. the films and substrate.
2.2. Isolated single ridge
Here, we study the eigen-vibrations of a single monomer ridge. In other words, the particle is completely isolated which enables to retrieve the modes of the ridge independently of any exterior perturbation by opposition to the periodic structure. This is performed by surrounding the ridge atop the multilayer with PMLs all along except the top surface which is set to vibrate freely; this enables to absorb incoming waves towards the boundaries (PMLs) with no reflections. Numerically speaking, we solve the same eigen-value problem where the Bloch periodic conditions are replaced by PMLs as boundary conditions. To track the modes, we compute the integral of the square modulus of the displacement field in the monomer normalized to its value in the whole system, namely: (1) (1) where, represents the displacement field with components along the spatial axes and respectively (which is adopted throughout the paper). Noting that, is the equivalent of the local density of states (LDOS) (Economou, Citation2006), the computed quantity is proportional to the normalized integral of LDOS within the ridge, which yields the localized eigenmodes of the latter. In , we give the simulations results of versus the eigen-frequency. The plot shows the existence of three peaks at the frequencies Such peaks match respectively the modes (A1, A2, A3/A4) studied previously (). In particular, it is interesting to note that modes A3 and A4 now appear as one mode (at ) namely the second order flexural mode of the ridge, this is due to the non-coupling with neighboring cells as the ridge is now isolated.
One observes that the frequencies of the modes are slightly different comparing with the periodic structure except for the first flexural mode (). As mentioned earlier, there is some coupling in modes A2 and A3/A4 with neighboring cells while for mode A1 the coupling is weak. As a result, for the isolated ridge, the frequency of the latter is almost the same whereas for the other modes there is a slight shift. Besides, the relative height of the peaks shows that the first-order flexural mode (P1) is the most localized mode (i.e. the elastic energy if localized solely in the ridge) followed by the higher-order flexural mode and then the compressional mode where part of the elastic energy is diffused in the surface layers and the substrate under the ridge.
Now, we use the SAW-platform in order to excite the ridge mechanical eigen-modes. The structure design is depicted in . The geometrical parameters are the same as in . It should be pointed-out that we perform a frequency response analysis of the structure, where the excitation of SAW is realized using a harmonic vertical line force at the system inlet on top (see ); where is the frequency and is the time. shows an example of SAW excitation at around in which we have depicted map, where represents the maximum value of the displacement norm in the whole system. We see the effective excitation of the SAWs and their complete absorption inside the PMLs regions preventing any undesired reflections. Let us indicate that the normalized transmission coefficient of the structure is computed by taking a line average of the registered displacement field magnitude near the output. The normalization is performed with respect to the system in absence of the ridges. shows the simulated transmission spectrum computed in the frequency range The curve shows the excitation of three modes at the frequencies (M1), (M2) and (M3). The frequencies match very well the position of the peaks shown in . These modes correspond to the first-order flexural mode (M1), the compressional mode (M2), and the higher-order flexural mode (M3). These results show clearly, the efficient excitation of the localized phonons of the single ridge using SAWs.
3. Dimer ridges eigenmodes
3.1. Periodic structure
We analyse in this paragraph the eigenmodes of a coupled dimer made of two identical ridges having the same dimensions as in , and separated with the distance For that, we follow the same approaches as in the case of the monomer ( and ) namely, in a first step, the ridges are placed atop the system in a unit-cell repeated periodically along axis, as shown in . shows the simulated dispersion curves (of the system in ). We notice the existence of six branches, labelled as B1 up to B6, underneath the SiO2-line with the corresponding frequencies at the edge of the first Brillouin zone (), respectively, and One observes a split-up of the modes () associated with a single monomer, yielding the modes (B1, B2), (B3, B4) and (B5, B6); these are the in- and out-of-phase modes of the first-order flexural mode (A1), the compressional mode (A2) and the second-order flexural mode (A3), respectively. This splitting occurs due to the strong coupling within the dimer, which results in the hybridization of the individual modes associated with each monomer. This can be seen in the displacement field norm maps near the dimer at the frequencies of the latter modes (B1-B6), as depicted in . By opposition with the monomer, the periodicity effect (coupling with neighbouring cells) at the second order flexural modes (B5, B6) of the dimer is different as no splitting is observed in the frequency range of interest. Also, the out-of-phase compressional mode B4 is influenced by periodicity as it is clearly entrained in a flexion-like movement induced by the films vibrations, in a Sezawa type mode.
3.2. Isolated dimer ridges eigenmodes
Similar to the previous case of a monomer, the eigen modes of the dimer are studied using the geometry model sketched in , except that the Bloch periodic conditions are replaced by PMLs, and the top surface is set free. The determination of the eigenmodes is done as earlier, that is, by computing the normalized integral of the equivalent of LDOS namely: (2) (2)
Note that the integration is obviously naturally performed within the coupled dimer, that is, the two ridges at the same time. In , we present the simulated results, where we notice the peaks labelled as P1 up to P6 at the frequencies, respectively. In the inset of the figure (near the top), we plot the spectrum of the average value of the square modulus of the displacement field in the dimer in the frequency range, This slightly different way of finding localized modes, enabled us to highlight the in-phase compressional mode P3, which is otherwise almost missing in the plot of The reason for this is that such a mode strongly radiates and couples with the surrounding films and substrate, making it a less localized mode to the point that its peak height is nearly not visible in . Overall, one notes a good match of modes P1-P6 with the modes obtained with the periodic lattice, B1-B6. There is slight shift in frequency depending on the mode which can be attributed as before to the periodicity effect or to the coupling with neighboring cells.
At this stage, we would like to show the possibility of exciting the localized acoustic phonons of the dimer using SAWs. The SAWs-based platform, where a dimer is atop the multilayer is shown in . The frequency response analysis of the system is performed as previously using vertical line force on top near the inlet (see ), to excite SAWs. The normalized transmission spectrum of the dimer-based SAWs design is computed is using the same procedure as mentioned earlier. For a matter of comparison, we overlay in the transmission of the monomer (blue curve) and that of the dimer (red curve).
The dips observed in case of the dimer (red curve) are labelled D1 up to D6; they fall at the frequencies, and respectively. According to the displacement field norm maps at these frequencies, the dips correspond, respectively, to the in-and out-of-phase first flexural mode, compressional mode and second order flexural mode. The modes’ frequencies match well the frequencies positions of the peaks in . The comparison with the transmission of the monomer shows that each dip associated with an eigenmode of the monomer is split-up into two dips. Indeed, the elastic coupling in the dimer induces a split of the modes M1, M2 and M3 of the monomer, respectively, into the modes, (D1, D2), (D3, D4), and (D5, D6), of the dimer. These results show clearly the efficient excitation of the dimer mechanical eigen-modes using SAWs.
4. Phonon modes selection for acousto-plasmonic coupling purposes
In this last section, we want to provide insights regarding the dimer-structure localized acoustic phonons modes that may be beneficial for acousto-plasmonic coupling. It is important to note that dimers-based designs have been shown to have superior optical (Li, Zhang, & Lei, Citation2016), elastic and optomechanical properties (Matheny, Citation2018) comparing with monomers, making them the best choice for enhanced phonon-plasmon interaction. As a matter of fact, we have demonstrated in this work that dimers display a richer spectrum of localized phonons in contrast to monomers, due to modes splitting of the latters. Besides, it has been found that the out-of-phase flexural modes of the dimer turned out to be optomechanically active while they become inactive in case of monomers, and the in-phase flexural and out-of-phase compressional modes are proven inactive (Noual, Akiki, et al., Citation2021). In the light of this, we shall analyse the acousto-plasmonic properties of the following localized phonons of the proposed dimer-system namely: the out-of-phase flexural modes (B2 and B6) and the in-phase compressional mode (mode B3).
In that regard, we would like to mention that the polymer material underneath the dimer has been shown to act as an effective optical cavity (see the encircled region in ) where localized surface plasmons could be excited (Mrabti et al., Citation2016). Thus, we examine the deformation induced by the dimer phonons modes on this specific region. The modes that deform this area volume shall cause important local strain, which by virtue of the photoelastic and/or moving boundary effects should cause a shift in the plasmon frequency. In other words, such phonons can be considered as optomechanically active. In , we plot and components of the displacement field normalized to the max value of the latter in the system, and associated with modes, B2, B3 and B6, respectively. In both first order () and second order out-of-phase flexural () modes (B2 and B6), the cavity undergoes a symmetrical elongation with respect to the symmetry plane at the middle, for both and resulting into an important deformation of its overall volume; hence such modes should be optomechanically active. In the in-phase compressional B3 (), the cavity is compressed according to and symmetrically elongated based on the -field, which indicates that this mode should also be active, in-line with earlier works (Noual et al., Citation2020). We have omitted to discuss the dimer phonons effect on probable localized surface plasmons within the polymer layer in the middle of each ridge. The reason is that except for the in-phase compressional mode which we believe should deform such a polymer layer where the plasmon would be localized, the other modes require further investigations, including the optical characterization of such plasmons. We intend to do this in a forthcoming work.
5. Conclusion
In conclusion, we have performed numerical simulations to investigate localized acoustic phonons supported by a gold dimer nanoridges-based multilayer design using finite element method. The structure design consists of a coupled gold ridges dimer set atop a thin polymer layer coating a thick gold film that is deposited over a semi-infinite SiO2-substrate. We have performed eigen-frequency analysis of the system where the ridges are either isolated in a finite system bounded by perfect matched layers or set to form a periodic lattice which enables to get the dispersion curve. In both cases, we have analysed the dimer eigenmodes by comparison with a single monomer. In the latter case, we found flexural and compressional type modes in the sub-GHz region. The elastic coupling in the dimer results in the split of the monomer modes, yielding in- and out-of-phase flexural and compressional modes. We have demonstrated their efficient excitation using a SAW-based device as important dips associated with the dimer modes show in the transmission coefficient of the system. Performing a symmetry analysis of the modes based on their and components, we can deduce that out-of-phase flexural and in-phase compressional modes should be optomechanically active. The proposed dimer-based SAW-device can be used for designing acousto-optic modulators or ultrasensitive sensors.
Disclosure statement
No potential conflict of interest was reported by the authors.
Correction Statement
This article has been corrected with minor changes. These changes do not impact the academic content of the article.
References
- Aram, M. H., & Khorasani, S. (2018). Optomechanical coupling strength in various triangular phoxonic crystal slab cavities. Journal of the Optical Society of America B, 35(6), 1390–1396. doi:10.1364/JOSAB.35.001390
- Arregui, G., Ng, R. C., Albrechtsen, M., Stobbe, S., Sotomayor-Torres, C. M., & García, P. D. (2023). Cavity Optomechanics with Anderson-Localized Optical Modes. Physical Review Letters, 130(4), 043802. doi:10.1103/PhysRevLett.130.043802
- Balram, K. C., Davanço, M., Lim, J. Y., Song, J. D., & Srinivasan, K. (2014). Moving boundary and photoelastic coupling in GaAs optomechanical resonators. Optica, 1(6), 414–420. doi:10.1364/OPTICA.1.000414
- Barzanjeh, S., Xuereb, A., Gröblacher, S., Paternostro, M., Regal, C. A., & Weig, E. M. (2022). Optomechanics for quantum technologies. Nature Physics, 18(1), 15–24. doi:10.1038/s41567-021-01402-0
- Cheng, J., Cen, Z., & Xu, S. (2021). Learning Nonlinear Waves in Plasmon-induced Transparency. ArXiv, Abs/2108.01508. 10.48550/arXiv.2108.01508.
- Dobrzynski, L., El Boudouti, E. H., Akjouj, A., Pennec, Y., Al-Wahsh, H., Leveque, G., & Djafari-Rouhani, B. (2017). Phononics: interface transmission tutorial book series. Amsterdam: Elsevier.
- Economou, E. N. (2006). Green’s functions in quantum physics (Vol. 7). 3rd Edition, Springer-Verlag, Berlin. doi:10.1007/3-540-28841-4
- Gao, R., He, Y., Zhang, D., Sun, G., He, J. X., Li, J. F., … Yang, Z. (2023). Gigahertz optoacoustic vibration in Sub-5 nm tip-supported nano-optomechanical metasurface. Nature Communications, 14(1), 485. doi:10.1038/s41467-023-36127-6
- Kelly, L., Northfield, H., Rashid, S., Bao, X., & Berini, P. (2022). Fabrication of high frequency SAW devices using tri-layer lift-off photolithography. Microelectronic Engineering. 253, 111671. doi:10.1016/j.mee.2021.111671
- Laude, V., Belkhir, A., Alabiad, A. F., Addouche, M., Benchabane, S., Khelif, A., & Baida, F. I. (2017). Extraordinary nonlinear transmission modulation in a doubly resonant acousto-optical structure. Optica, 4(10), 1245–1250. doi:10.1364/OPTICA.4.001245
- Li, G. C., Zhang, Q., Maier, S. A., & Lei, D. (2018). Plasmonic particle-on-film nanocavities: A versatile platform for plasmon-enhanced spectroscopy and photochemistry. Nanophotonics, 7(12), 1865–1889. doi:10.1515/nanoph-2018-0162
- Li, G.-C., Zhang, Y.-L., & Lei, D. Y. (2016). Hybrid plasmonic gap modes in metal film-coupled dimers and their physical origins revealed by polarization resolved dark field spectroscopy. Nanoscale, 8(13), 7119–7126. doi:10.1039/C5NR09260D
- Lin, T. R., Lin, C. H., & Hsu, J. C. (2015). Strong optomechanical interaction in hybrid plasmonic-photonic crystal nanocavities with surface acoustic waves. Scientific Reports, 5(1), 13782. doi:10.1038/srep13782
- Matheny, M. H. (2018). Enhanced photon-phonon coupling via dimerization in one-dimensional optomechanical, crystals. Appl. Phys. Lett 18, 112(25), 253104. doi:10.1063/1.5030659
- Mei, J., Zhang, N., & Friend, J. (2020). Fabrication of surface acoustic wave devices on lithium niobate. Journal of Visualized Experiments, 160(160), e61013. doi:10.3791/61013
- Mrabti, A., Lévêque, G., Akjouj, A., Pennec, Y., Djafari-Rouhani, B., Nicolas, R., … Adam, P. M. (2016). Elastoplasmonic interaction in metal-insulator-metal localized surface plasmon systems. Physical Review B, 94(7), 075405. doi:10.1103/PhysRevB.94.075405
- Noual, A., Akiki, R., Lévêque, G., Pennec, Y., & Djafari-Rouhani, B. (2021). Enhanced phonon-plasmon interaction in film-coupled dimer nanoridges mediated by surface acoustic waves. Optics Express, 29(26), 43104–43123. doi:10.1364/OE.444430
- Noual, A., Akiki, R., Pennec, Y., El Boudouti, E. H., & Djafari-Rouhani, B. (2020). Surface acoustic waves-localized plasmon interaction in pillared phononic crystals. Physical Review Applied, 13(2), 024077. doi:10.1103/PhysRevApplied.13.024077
- Noual, A., Kang, E., Maji, T., Gkikas, M., Djafari-Rouhani, B., & Fytas, G. (2021). Optomechanic coupling in ag polymer nanocomposite films. The Journal of Physical Chemistry. C, Nanomaterials and Interfaces, 125(27), 14854–14864. doi:10.1021/acs.jpcc.1c04549
- Oo, S. Z., Silva, G., Carpignano, F., Noual, A., Pechstedt, K., Mateos, L., … Melvin, T. (2016). A nanoporous gold membrane for sensing applications. Sensing and Bio-sensing Research 7, 133–140. doi:10.1016/j.sbsr.2016.01.001
- Royer, D., & Dieulesaint, E. (2000). Elastic waves in solids. New York: Springer. doi:10.1002/9781119902942
- Saada, T. N., Marques da Silva, A. G., Subramanian, P., Pang, L., Adnane, N., Djafari-Rouhani, B., … Szunerits, S. (2020). Plasmon-enhanced electrocatalytic oxygen reduction in alkaline media on gold nanohole electrodes. Journal of Materials Chemistry A, 8(20), 10395–10401. doi:10.1039/C9TA14174J
- Sansa, M., Defoort, M., Brenac, A., Sansa, M., Defoort, M., Brenac, A., … Hentz, S. (2020). Optomechanical mass spectrometry. Nature Communications, 11(1), 3781. doi:10.1038/s41467-020-17592-9
- Vasileiadis, T.,Noual, A.,Wang, Y.,Graczykowski, B.,Djafari-Rouhani, B.,Yang, S., &Fytas, G. (2022). Optomechanical hot-spots in metallic nanorod-polymer nanocomposites. ACS Nano, 16(12), 20419–20429. 10.1021/acsnano.2c06673. 36475620