Multipolar scattering analysis of hybrid metal-dielectric nanostructures
Abstract
We perform a systematic study showing the evolution of the multipoles along with the spectra for a hybrid metal-dielectric nanoantenna, a Si cylinder and an Ag disk stacked one on top of another, as its dimensions are varied one by one. We broaden our analysis to demonstrate the "magnetic light" at energies above 1 eV by varying the height of the Ag on the Si cylinder and below 1 eV by introducing insulating spacing between them. We also explore the appearance of the anapole state along with some exceptionally narrow spectral features by varying the radius of the Ag disk.
1 Introduction
Metallic nanostructures can produce a strong field enhancement by the resonant oscillation of free electrons and have thus been a popular research topic for several decades [1]. There are however two main limitations to those resonances: the intrinsic losses associated with metals, which limit the quality factor that can be achieved [2]; and their essentially electric dipolar character, which restricts the degrees of freedom for engineering arbitrary electromagnetic responses [3]. These issues can be solved by replacing metals with low loss dielectrics, which have recently emerged as a vivid research topic [4, 5, 6, 7, 8]. Dielectric nanostructures exhibit both electric and magnetic resonances, which enable several interesting phenomena like unidirectional scattering [9, 10] Fano resonances [11], or even anapole states that do not radiate in the far field and are formed by the destructive interference between electric and toroidal dipoles [12]. Unfortunately, despite all their advantages, the use of dielectrics reduces the field enhancement in the vicinity of the structure, which is detrimental for applications like fluorescence [13] and sensing [14]. In order to circumvent the drawbacks of both metals and dielectrics, hybrid metal-dielectric nanostructures have recently come up as a promising option. These nanostructures are expected to provide the best of both worlds by combining low losses from dielectric and high field enhancement from metal.
The significant progress made in various facets of hybrid nanophotonics has been efficiently collated in a review by Lepeshov et al. [15], which illustrates experimentally achievable hybrid geometries based on different working principles like the whispering gallery modes in Si-Al2O3-Ag dimers leading to directional radiation [16] and the coupling of a hybrid mode at the Si-Ag interface of a Si nanoparticle on an Ag film to the propagating surface plasmon-polariton [17]. Unconventional and newer hybrid core-shell geometries like an Ag rod covered by a semiconductor layer for broadband absorption enhancement [18] and a gold rod covered by a Cu2O shell for generation of third harmonic [19], to name a few, has also been mentioned. The hybrid nanoantenna, due to the presence of dielectrics, boasts of some unique features like the pure magnetic dipole resonance or the anapole states, which will also be seen in the later sections of this work. Predominantly Magnetic Fano resonances have also been observed in hybrid oligomers formed from a combination of Si cones and Au particles, which can be tuned by laser melting of the Au particles [20, 21]. Additionally there are several other recent works on hybrid nanoantenna that also need to be acknowledged such as the ones that discuss the efficient generation of the second [22, 23], and third harmonics [24]; nonlinear optical properties [25], subtractive color filtering [26, 27], sensing [28], unidirectional radiation[29] and asymmetric absorption [30].
Upon careful considerations of these articles it can be noted that core-shell or stacked cylindrical geometries are the most popular. While the core-shell geometry has analytical solutions, the stacked cylindrical geometry still lacks a consistent theoretical tool for the analysis. One possible way to fulfill this gap is to perform a comprehensive systematic parametric study on how the geometrical parameters of the system influence the scattering characteristics of these structures, which represents the focus of this paper.
To the author’s knowledge, only a few publications have unveiled some of these dependencies. For example, the scattering characteristics of dielectric core coated with metal shell have been implemented experimentally and studied analytically for spherical and numerically for co-axial cylindrical configurations [8, 31, 32, 33]. A hybrid geometry with metal nanorods on top of Si cylinders was studied, noting the influence of a spacer between a metal rod and a dielectric cylinder to enhance directional radiation [34, 35]. The effects of thickness and permittivity of the dielectric spacer and the metal nanodisk in a metal-dielectric-metal sandwiched geometry on the directional radiation and manipulation of multipoles were also investigated [36, 37].
In principle, the coupled electric and magnetic dipole theory that works for spherical objects [38] could be applied to predict scattering from cylindrical geometries. Unfortunately, from our experience, this approach does not give sufficiently accurate results and we have to resort to numerical methods to explain how various features of these hybrid geometries appear and disappear as we perform a systematic study of one specified hybrid geometry by varying its dimensions. For our work we chose a metal-dielectric nanostructure made of Si cylinder and Ag disk placed on top of one another. We start by varying the geometrical parameters like height and radius of a Si cylinder, followed by Ag disk and finally that of a hybrid structure. Since the parameters considered here are experimentally meaningful, this systematic analysis showing how each geometrical parameter affects the spectral response is also very beneficial from an engineering perspective for designing hybrid geometry. This analysis and the underlying physics can also be extended to several other similar hybrid geometries like core-shell particle [39, 40, 41] and cylinder-disk [26, 27] and cylinder-nanorod geometries [34, 35].
We perform the numerical simulations with the surface integral method, which has been proven a very accurate frequency domain technique [42, 43]. To avoid the singularities of the Green’s tensor close to the scatterer’s surface, we use the singularity subtraction technique [44, 45]. Refs. [43] and [44] discuss in detail the performance of this technique and the accuracy of the algorithm in computing the response from various nanostructures. Numerical simulations are challenging and choosing a method that is best suited to the needs of the analysis with minimal trade offs is always an art. In this work we have presented all the results of the parametric study using a well-established SIE method. The use of a single method throughout the paper helps in building a coherent understanding as the geometry progresses from a single Si cylinder to a more complicated Ag disk-Si cylinder hybrid nanoantenna. The Cartesian multipoles decomposition is performed with the vector spherical harmonic basis [46], which is available online [47, 48, 49]. The resultant multipole for any structure strongly depends on the choice of origin of the coordinate system [50]. We make the choice of origin at the geometrical center of the mesh to provide a coherent analysis. Aiming at clarifying how different geometrical parameters affect the overall spectrum, we vary the size of each component and monitor the evolution of the multipolar response traces. The system is illuminated with a planewave propagating along the symmetry axis z, leading to a polarization insensitive scattering response. Such a geometry can be fabricated with a high level of control [30, 28]. We have chosen Si as the dielectric material since it has high permittivity and modest losses in the visible and has been widely used experimentally [5, 51, 52] and Ag as the metal for its low losses [53, 54, 55]. For our simulations we use available data for Si from Ref.[56]. To put an effort to be consistent with our recent experimental research on hybrid antenna for sensing, we use our experimentally measured data for Ag [57] that, in comparison to the data from Johnson and Christy [58], show almost the same real part and slightly lower imaginary part of the dielectric permittivity of silver in the region 1100-1700 nm.
2 Results and discussion
The results are divided into several sections to explore in details the influence of the different geometrical parameters on the spectral features supported by the system. We first explore the scattering cross-section (SCS) for a bare Si cylinder, then add a metal disk of varying thickness and diameter. Finally, we investigate the coupling between the dielectric and metallic parts as a function of their separation.
2.1 Response of the dielectric cylinder
To understand the advantages of a hybrid nanoantenna and the additional control it offers, we begin by analyzing the SCS of an isolated Si cylinder illuminated with a planewave. Although this has been reported by others [59, 60, 12, 61, 62], it will help us form a coherent understanding as we proceed to the subsequent sections, where we modify the geometry gradually by adding Ag to the Si cylinder.
We first vary the radius of the Si cylinder from 235 nm until 47 nm in steps of 47 nm, as shown in Fig. 1(a); for the time being, we retain the height of the Si cylinder as 220 nm. Already in this simple system, we observe a rather rich spectrum with many different features, especially for the larger cylinders. The entire SCS can be decomposed into an electric dipole (ED), a magnetic dipole (MD) and an electric quadrupole (EQ) Cartesian multipoles. It can be observed that with the decrease in radius, both the ED and MD blueshift. Interestingly, the ED is shifting more than the MD, such that the ED leaves the MD far away behind, on the left side of the spectrum, as the radius decreases, as seen for nm in Fig. 1(a). This faster shift of the ED is actually not surprising and has been observed previously for the case of a dielectric sphere (see e.g. Fig. 1(c) in Ref. [15]), and is explained by analyzing the electric and magnetic Mie coefficients and in the framework of the Padé expansion[7]. The magnitudes of both ED peaks decrease as the cylinder shrinks, while the MD resonance retains its magnitude.

Let us now vary the height of the Si cylinder from 220 nm to 44 nm in steps of 44 nm, Fig. 1(b). Again, the ED and MD resonance positions blueshift with decreasing sizes. Opposite to what was observed for a change of radius, now the ED magnitudes remain rather constant, while the MD magnitude drastically decreases as we diminish the height. Also, we notice that the position of the MD peak blueshifts much faster than the ED ones, leading to significant "overtaking" of the MD resonance over the ED, as seen for nm in Fig. 1(b). Actually, this observation is important for engineering directional scattering [63, 64, 9, 10], or producing a far-field response constituting of only one particular mode, e.g. the MD for "magnetic light" [65, 31]. The ED and MD resonances of dielectric structures can have strong spectral overlap between them; a difference in their shift rate enables their efficient spectral separation. For example, by varying the Si height, we can control the ED and MD responses and achieve a good spectral separation for nm, Fig. 1(b).
The different spectral behaviour for ED and MD stems from the physical nature of these multipoles. The MD response is known to be produced by a current loop induced in this particular case in the xz-plane. With the decrease in height, the effective space volume pertaining to the current loop decreases, reducing the resonance energy of the MD mode. On the contrary, the ED energy is related to the bulk polarization distributed mostly in the yx-plane inside the structure. Consequently, because of the asymmetry of the structure with respect to the z-axis, these two multipoles are affected differently as we vary radius and height [61].
In order to get a complete picture, we also present the spectra attributed to the simultaneous variation of height and radius in Fig. 1(c). In this case, the blueshift rate is almost the same for the ED and MD. Also, we observe a simultaneous magnitude reduction for both resonances as we decrease the cylinder volume.
Figure SI 1 shows the position of the peaks obtained from Fig. 1 for the ED and MD resonances as a function of the geometrical parameters. Here, the circles stand for the resonance peaks at lower energies while the diamonds stand for the resonance peaks at higher energies. It should also be noted that in this plot we have no data pertaining to nm and nm since there is no resonance peak for the ED or the MD in the considered energy range, Fig. 1(c). At 1.07 eV we observe a dip in the dipolar trace for nm nm, which is associated with the excitation of the anapole state [66]. This rather unusual mode stems from the destructive interference between the ED and the toroidal dipole (TD)[12, 67, 68]. This is confirmed by plotting the SCS of the ED and TD, which indeed intersect at this wavelength with a phase difference of (data not shown). Interestingly, the anapole state also blueshifts more with decreasing radii, as compared to the decrease in height, since the anapole state appears due to the destructive interference of ED and TD; since the ED shifts more with change in radius, so does the anapole state.
2.2 Influence of the metal disk

Having studied the effects caused by increasing the dimensions of the Si cylinder, we now study the effects caused by adding a metallic component in the system, i.e. by appending a few nm thick metallic disk on top of the Si cylinder. Indeed, quite significant level of control over the position of the resonances can be achieved by introducing a metal in the system [31, 69]. We choose a Si cylinder with nm and nm since it has prominent additional multipoles, like MD and EQ. The advantage of involving these multipoles in the analysis can be appreciated in the following where "magnetic light" is achieved. We start with the Ag thickness nm, as this is experimentally achievable [70, 71, 72], and increase the thickness of Ag up to 60 nm [28], as shown in Fig. 2.
An addition of metal, albeit extremely thin, drastically modifies the spectrum. In contrast to the dielectric cylinder, where we were able to only manipulate the position and the magnitude of the resonances, we are now able to modify the profile of the resonance itself. Indeed, by inspecting the dipolar trace of the spectrum achieved by adding only 10 nm of Ag we notice sudden dips at 0.91, 1.18 and 1.34 eV, compare and in Fig. 2(b). The dip at 0.91 eV is fairly sharp with a full width at half-maximum of 0.01 eV. The emergence of these dips can be explained by the interference between two -shifted Cartesian electric dipoles excited in the metal and dielectric parts, giving rise to the effective excitation of the EQ and MD resonances in the far-field [28]. Quite interestingly, this effect of efficient translation of two ED Cartesian multipoles into EQ and MD is seen only at particular frequencies. Moreover, it becomes less noticeable as the metal thickness increases; possibly because the intrinsic losses in the metal take over.
From the findings for a dielectric cylinder, it is expected to observe a redshift of the MD resonance as the height of the structure is increased. Surprisingly, here the overall profile of the MD spectrum retains its position and magnitude, once the Ag thickness is increased above 15 nm. In the meanwhile the position and the profile of the ED resonance is subject to more noticeable changes. Thanks to these features, almost pure MD resonance can be achieved in this geometry at 1.11 eV, Fig. 2 for nm. Similar results have been reported for the geometries of two co-axial metal-dielectric cylinders [8] and two metallic cylinders placed one on top of another and separated by a dielectric [36]. The decoupling of the ED and MD in the range of 1-1.2 eV is also very beneficial for the generation of an excellent anapole state. For example, for Ag thicknesses between nm and nm in Fig. 2 we observe around 1.1 eV negligible ED and MD components, but still a sufficiently strong MQ contribution. It has been shown that, by using radially polarized excitation instead of linearly polarized one, the magnetic contributions can be suppressed [31, 73]. This will result in an excellent anapole state making the structure non-radiating at this wavelength. The broad dip in the ED at 1.1 eV for Ag disk of height nm (Fig. 2) is attributed to the anapole state, which leads to the minimum in the ED response and so only a contribution from the MQ is seen [12]. This anapole state suffers negligible blueshifting with an increase in the thickness of Ag. Also, we noticed that increasing the height of silver above 60 nm barely affects the spectral features of the system. Thus, we take the Si cylinder with 60 nm thick Ag disk as a final structure on which we carry out further modifications to understand the interaction between the Si cylinder and the Ag disk.

The spectra reported in Fig. 2 have multiple prominent resonant features, which distinguish these hybrid nanostructures from their pure metal or dielectric counterparts. However, sometimes there is a need for a spectrum with only one or two sharp features in the SCS, e.g. for sensing when the change in the refractive index of the background analyte is measured by tracking one or two spectral features [74]. In this case, the sharper the resonance, the higher is their figure of merit. This can be achieved in the nanostructure at hand by varying the radius of the Ag disk as shown in Fig. 3 for Ag disks with radii between nm and 235 nm (for clarity, the baseline for the Ag disk SCS is shifted by 1m2). We notice that the MD position can be adjusted within the range of 1.08-1.35 eV by changing the Ag disk radius. Specifically, the magnitude of both MD resonances within this region is much stronger than in the remaining part of the spectrum. For nm, we find the MDs at 1.15 and 1.24 eV, Fig. 3. For larger Ag radii, the lower energy MD remains around 1.15 eV, while the higher energy MD slightly redshifts for nm, and then begins to fuse with the low energy MD, producing a very strong spectral feature for nm. This feature becomes even more prominent for nm, where it dominates the entire SCS, Fig. 3. This peak can also be explained by the interference between two -shifted Cartesian electric dipoles excited in the metal and dielectric parts giving rise to effective excitation of EQ and MD resonances in the far field. The MD resonance arising from this phenomenon interferes with the MD from the Si cylinder and reduces the net MD of the system [28]. Such a very sharp feature could be very useful for sensing applications as it will give high figure of merit [74, 75]. Beyond nm, this sharp feature at 1.15 eV starts again to split into two peaks as shown in Fig. 3 for and 235 nm.
2.3 Coupling between the dielectric cylinder and the metal disk
Finally, we investigate an additional control over the spectrum offered by this hybrid nanostructure by introducing a spacing of varying thickness between the Si cylinder and the Ag disk. For the ease of understanding, we consider an air spacer; a dielectric spacer will merely increase the coupling as is the case for a dielectric loaded plasmonic structure [76]. This study complements the one done previously for Si cylinder sandwiched between two Ag disks [36]. In Fig. 4 we slowly increase the spacing from nm, which is achievable with current standards of nanofabrication [77], to nm.

The most prominent feature that can be noticed as we increase the spacing is the appearance of a broader resonance in the SCS at 0.74 eV due to the emergence in the spectral window of a very strong MD resonance for a separation larger than nm. Such a feature in the near infrared can be useful for applications in surface enhanced infrared absorption [78]. It can also be used in applications requiring only the magnetic response from the structure. However, overall, the spectral features observed for large spacings are relatively broad and will not be well suited for practical applications. Another feature in the spectrum attracted our attention and is worth mentioning. The two peaks at 1.09 eV and 1.16 eV in the SCS shown in Fig. 4 for , are relatively sharp. By increasing the spacing between dielectric and metal, the peak at higher energy becomes less pronounced while the magnitude of that at lower energy significantly increases. This is in contrast to the effect shown in Fig. 2, where adding Ag on top of the Si cylinder introduced sharp small dips in the ED spectra in the range of 1-1.5 eV which affected the SCS in this range likewise.
Nanostructure | Parameter | ED | MD | ||
Energy | Amplitude | Energy | Amplitude | ||
Si cylinder | |||||
& | |||||
Hybrid structure | |||||
1 Sharp optical features.
2 Both MD fuse to give rise to a sharp peak at nm and then split again.
3 Strong interaction between ED and MD.
3 Conclusion
To conclude, a hybrid dielectric–metal nanostructure provides many degrees of freedom to control different spectral features, including pure magnetic modes and anapoles, as well as extremely narrow resonances. The results of the geometrical variations are summarised in Table. 1. We see that, depending on the geometry, the ED and MD resonances can be spectrally overlapped or separated and the knowledge gathered in this study could be combined with advanced design approaches to develop systems that support a given set of multipolar characteristics[79, 37]. Overall, hybrid metal-dielectric nanostructures provide both a very high degree of freedom and a good control for engineering resonant systems with tailored spectral characteristics.
Funding Funding from the Swiss National Science Foundation (project 200021-162453) and the European Research Council (ERC-2015-AdG-695206 Nanofactory) is gratefully acknowledged.
Acknowledgments It is a pleasure to acknowledge stimulating discussions with G. D. Bernasconi during the initial phase of this project.
Disclosures The authors declare no conflicts of interest.
Data availability Datasets underlying the results presented in this paper are available in Refs. [47, 57].
Supplemental document See Supplement 1 for supporting content.
References
- [1] S. Maier, Plasmonics: Fundamentals and Applications (Springer US, 2007).
- [2] J. P. Kottmann, O. J. F. Martin, D. R. Smith, and S. Shultz, “Field polarization and polarization charge distributions in plasmon resonant nanoparticles,” \JournalTitleNew Journal of Physics. 2 (2000).
- [3] Y. Ekinci, A. Christ, M. Agio, O. J. F. Martin, H. H. Solak, and J. F. Löffler, “Electric and magnetic resonances in arrays of coupled gold nanoparticle in-tandem pairs,” \JournalTitleOpt. Express 16, 13287–13295 (2008).
- [4] Y. Yang, I. I. Kravchenko, D. P. Briggs, and J. Valentine, “All-dielectric metasurface analogue of electromagnetically induced transparency,” \JournalTitleNature Communications 5, 5753 (2014).
- [5] A. I. Kuznetsov, A. E. Miroshnichenko, M. L. Brongersma, Y. S. Kivshar, and B. Luk’yanchuk, “Optically resonant dielectric nanostructures,” \JournalTitleScience 354 (2016).
- [6] O. Yavas, M. Svedendahl, P. Dobosz, V. Sanz, and R. Quidant, “On-a-chip biosensing based on all-dielectric nanoresonators,” \JournalTitleNano Letters 17, 4421–4426 (2017).
- [7] D. Tzarouchis and A. Sihvola, “Light scattering by a dielectric sphere: Perspectives on the mie resonances,” \JournalTitleApplied Sciences 8 (2018).
- [8] T. Feng, Y. Xu, W. Zhang, and A. E. Miroshnichenko, “Ideal magnetic dipole scattering,” \JournalTitlePhys. Rev. Lett. 118, 173901 (2017).
- [9] A. Kiselev, K. Achouri, and O. J. F. Martin, “Multipole interplay controls optical forces and ultra-directional scattering,” \JournalTitleOpt. Express 28, 27547–27560 (2020).
- [10] K. Achouri, A. Kiselev, and O. J. F. Martin, “Multipolar origin of electromagnetic transverse force resulting from two-wave interference,” \JournalTitlePhys. Rev. B 102, 085107 (2020).
- [11] M. Limonov, M. Rybin, A. Poddubny, and Y. Kivshar, “Fano resonances in photonics,” \JournalTitleNature Photonics 11, 543–554 (2017).
- [12] A. E. Miroshnichenko, A. B. Evlyukhin, Y. F. Yu, R. M. Bakker, A. Chipouline, A. I. Kuznetsov, B. Luk’yanchuk, B. N. Chichkov, and Y. S. Kivshar, “Nonradiating anapole modes in dielectric nanoparticles,” \JournalTitleNature Communications 6 (2015).
- [13] D. Darvill, A. Centeno, and F. Xie, “Plasmonic fluorescence enhancement by metal nanostructures: shaping the future of bionanotechnology,” \JournalTitlePhys. Chem. Chem. Phys. 15, 15709–15726 (2013).
- [14] X. Han, K. Liu, and C. Sun, “Plasmonics for biosensing,” \JournalTitleMaterials 12, 1411 (2019).
- [15] S. I. Lepeshov, A. E. Krasnok, P. A. Belov, and A. E. Miroshnichenko, “Hybrid nanophotonics,” \JournalTitlePhysics-Uspekhi 61, 1035–1050 (2018).
- [16] Q. Wang, H. Zhao, X. Du, W. Zhang, M. Qiu, and Q. Li, “Hybrid photonic-plasmonic molecule based on metal/si disks,” \JournalTitleOpt. Express 21, 11037–11047 (2013).
- [17] Y. Yang, O. D. Miller, T. Christensen, J. D. Joannopoulos, and M. Soljačić, “Low-loss plasmonic dielectric nanoresonators,” \JournalTitleNano Letters 17, 3238–3245 (2017). PMID: 28441499.
- [18] S. A. Mann and E. C. Garnett, “Extreme light absorption in thin semiconductor films wrapped around metal nanowires,” \JournalTitleNano Letters 13, 3173–3178 (2013). PMID: 23758555.
- [19] O. Bar-Elli, E. Grinvald, N. Meir, L. Neeman, and D. Oron, “Enhanced third-harmonic generation from a metal/semiconductor core/shell hybrid nanostructure,” \JournalTitleACS Nano 9, 8064–8069 (2015). PMID: 26197215.
- [20] S. I. Lepeshov, D. A. Zuev, S. V. Makarov, V. A. Milichko, I. S. Mukhin, A. E. Krasnok, and P. A. Belov, “Manipulating fano resonance via fs-laser melting of hybrid oligomers at nanoscale,” \JournalTitleJournal of Physics: Conference Series 741, 012140 (2016).
- [21] S. Lepeshov, A. Krasnok, I. Mukhin, D. Zuev, A. Gudovskikh, V. Milichko, P. Belov, and A. Miroshnichenko, “Fine-tuning of the magnetic fano resonance in hybrid oligomers via fs-laser-induced reshaping,” \JournalTitleACS Photonics 4, 536–543 (2017).
- [22] V. F. Gili, L. Ghirardini, D. Rocco, G. Marino, I. Favero, I. Roland, G. Pellegrini, L. Duò, M. Finazzi, L. Carletti, A. Locatelli, A. Lemaître, D. Neshev, C. De Angelis, G. Leo, and M. Celebrano, “Metal–dielectric hybrid nanoantennas for efficient frequency conversion at the anapole mode,” \JournalTitleBeilstein Journal of Nanotechnology 9, 2306–2314 (2018).
- [23] H.-G. Hu, M. Yang, P. Yue, Y.-T. Bai, W.-J. Wang, and S.-D. Liu, “Second-harmonic generation with metal/dielectric/metal hybridized nanoantennas: enhanced efficiency, reduced mode volume and ideal magnetic/electric dipole scattering,” \JournalTitleJournal of Physics D: Applied Physics 53, 215101 (2020).
- [24] T. Shibanuma, G. Grinblat, P. Albella, and S. A. Maier, “Efficient third harmonic generation from metal–dielectric hybrid nanoantennas,” \JournalTitleNano Letters 17, 2647–2651 (2017).
- [25] M. Hentschel, B. Metzger, B. Knabe, K. Buse, and H. Giessen, “Linear and nonlinear optical properties of hybrid metallic–dielectric plasmonic nanoantennas,” \JournalTitleBeilstein Journal of Nanotechnology 7, 111–120 (2016).
- [26] W. Yue, G. Song, L. Sang-Shin, K. Eun-Soo, and C. Duk-Yong, “Subtractive color filters based on a silicon-aluminum hybrid-nanodisk metasurface enabling enhanced color purity,” \JournalTitleScientific Reports 6, 29756 (2016).
- [27] Y. Nagasaki, I. Hotta, M. Suzuki, and J. Takahara, “Metal-masked mie-resonant full-color printing for achieving free-space resolution limit,” \JournalTitleACS Photonics 5, 3849–3855 (2018).
- [28] D. Ray, T. V. Raziman, C. Santschi, D. Etezadi, H. Altug, and O. J. F. Martin, “Hybrid metal-dielectric metasurfaces for refractive index sensing,” \JournalTitleNano Letters 20, 8752–8759 (2020).
- [29] J. Ho, Y. H. Fu, Z. Dong, R. Paniagua-Dominguez, E. H. H. Koay, Y. F. Yu, V. Valuckas, A. I. Kuznetsov, and J. K. W. Yang, “Highly directive hybrid metal–dielectric yagi-uda nanoantennas,” \JournalTitleACS Nano 12, 8616–8624 (2018).
- [30] J.-H. Yang and K.-P. Chen, “Hybridization of plasmonic and dielectric metasurfaces with asymmetric absorption enhancement,” \JournalTitleJournal of Applied Physics 128, 133101 (2020).
- [31] P. Li, K. Du, F. Lu, K. Gao, F. Xiao, W. Zhang, and T. Mei, “Excellent anapole by decoupling electric multipoles of ag/si core–shell nanoparticles,” \JournalTitleThe Journal of Physical Chemistry C 124, 19252–19258 (2020).
- [32] C. Xu, K. Cheng, Q. Li, X. Shang, C. Wu, Z. Wei, X. Zhang, and H. Li, “The dual-frequency zero-backward scattering realized in a hybrid metallo-dielectric nanoantenna,” \JournalTitleAIP Advances 9, 075121 (2019).
- [33] W. Liu, A. E. Miroshnichenko, D. N. Neshev, and Y. S. Kivshar, “Broadband unidirectional scattering by magneto-electric core–shell nanoparticles,” \JournalTitleACS Nano 6, 5489–5497 (2012).
- [34] E. Rusak, I. Staude, M. Decker, J. Sautter, A. E. Miroshnichenko, D. A. Powell, D. N. Neshev, and Y. S. Kivshar, “Hybrid nanoantennas for directional emission enhancement,” \JournalTitleApplied Physics Letters 105, 221109 (2014).
- [35] R. Guo, E. Rusak, I. Staude, J. Dominguez, M. Decker, C. Rockstuhl, I. Brener, D. N. Neshev, and Y. S. Kivshar, “Multipolar Coupling in Hybrid Metal Dielectric Metasurfaces,” \JournalTitleACS Photonics 3, 349–353 (2016).
- [36] Y. Zhang, P. Yue, J.-Y. Liu, W. Geng, Y.-T. Bai, and S.-D. Liu, “Ideal magnetic dipole resonances with metal-dielectric-metal hybridized nanodisks,” \JournalTitleOpt. Express 27, 16143–16155 (2019).
- [37] F. Qin, D. Zhang, Z. Liu, Q. Zhang, and J. Xiao, “Designing metal-dielectric nanoantenna for unidirectional scattering via bayesian optimization,” \JournalTitleOpt. Express 27, 31075–31086 (2019).
- [38] S. Sun, D. Wang, Z. Feng, and W. Tan, “Highly efficient unidirectional forward scattering induced by resonant interference in a metal–dielectric heterodimer,” \JournalTitleNanoscale 12, 22289–22297 (2020).
- [39] R. S. Savelev, O. N. Sergaeva, D. G. Baranov, A. E. Krasnok, and A. Alù, “Dynamically reconfigurable metal-semiconductor yagi-uda nanoantenna,” \JournalTitlePhys. Rev. B 95, 235409 (2017).
- [40] W. Liu, A. E. Miroshnichenko, D. N. Neshev, and Y. S. Kivshar, “Broadband unidirectional scattering by magneto-electric core–shell nanoparticles,” \JournalTitleACS Nano 6, 5489–5497 (2012). PMID: 22545872.
- [41] H. Wang, P. Liu, Y. Ke, Y. Su, L. Zhang, N. Xu, S. Deng, and H. Chen, “Janus magneto–electric nanosphere dimers exhibiting unidirectional visible light scattering and strong electromagnetic field enhancement,” \JournalTitleACS Nano 9, 436–448 (2015). PMID: 25554917.
- [42] B. Gallinet, J. Butet, and O. J. F. Martin, “Numerical methods for nanophotonics: standard problems and future challenges,” \JournalTitleLaser & Photonics Reviews 9, 577–603 (2015).
- [43] A. M. Kern and O. J. F. Martin, “Surface integral formulation for 3d simulations of plasmonic and high permittivity nanostructures,” \JournalTitleJ. Opt. Soc. Am. A 26, 732–740 (2009).
- [44] T. V. Raziman, W. R. C. Somerville, O. J. F. Martin, and E. C. Le Ru, “Accuracy of surface integral equation matrix elements in plasmonic calculations,” \JournalTitleJournal of the Optical Society of America B 32, 485–492 (2015).
- [45] A. M. Kern and O. J. F. Martin, “Pitfalls in the determination of optical cross sections from surface integral equation simulations,” \JournalTitleIEEE Transactions on Antennas and Propagation 58, 2158–2161 (2010).
- [46] S. Mühlig, C. Menzel, C. Rockstuhl, and F. Lederer, “Multipole analysis of meta-atoms,” \JournalTitleMetamaterials 5, 64 – 73 (2011).
- [47] A. Kiselev, G. Bernasconi, and O. J. F. Martin, “Vector spherical decomposition code and nonlinear scattering data,” (2019).
- [48] A. Kiselev, G. D. Bernasconi, and O. J. F. Martin, “Modes interplay and dynamics in the second harmonic generation of plasmonic nanostructures,” \JournalTitleOptics Express 27, 38708–38720 (2019).
- [49] C. Yan, X. Wang, T. V. Raziman, and O. J. F. Martin, “Twisting fluorescence through extrinsic chiral antennas,” \JournalTitleNano Letters 17, 2265–2272 (2017).
- [50] L. M. Burko, “Multipole structure and coordinate systems,” \JournalTitleEuropean Journal of Physics 28, 593–602 (2007).
- [51] M. Decker, I. Staude, M. Falkner, J. Dominguez, D. N. Neshev, I. Brener, T. Pertsch, and Y. S. Kivshar, “High-efficiency dielectric huygens’ surfaces,” \JournalTitleAdvanced Optical Materials 3, 813–820 (2015).
- [52] I. Staude and J. Schilling, “Metamaterial-inspired silicon nanophotonics,” \JournalTitleNature Photonics 11, 274–284 (2017).
- [53] P. West, S. Ishii, G. Naik, N. Emani, V. Shalaev, and A. Boltasseva, “Searching for better plasmonic materials,” \JournalTitleLaser & Photonics Reviews 4, 795–808 (2010).
- [54] B. Dastmalchi, P. Tassin, T. Koschny, and C. M. Soukoulis, “A new perspective on plasmonics: Confinement and propagation length of surface plasmons for different materials and geometries,” \JournalTitleAdvanced Optical Materials 4, 177–184 (2016).
- [55] J. Khurgin, “Replacing noble metals with alternative materials in plasmonics and metamaterials: how good an idea?” \JournalTitlePhilosophical Transactions of The Royal Society A Mathematical Physical and Engineering Sciences 375 (2016).
- [56] D. T. Pierce and W. E. Spicer, “Electronic structure of amorphous si from photoemission and optical studies,” \JournalTitlePhys. Rev. B 5, 3017–3029 (1972).
- [57] D. Ray, A. Kiselev, and O. J. F. Martin, “Refractive index data for 150 nm Ag measured with ellipsometry,” (2021).
- [58] P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” \JournalTitlePhys. Rev. B 6, 4370–4379 (1972).
- [59] I. Staude, A. E. Miroshnichenko, M. Decker, N. T. Fofang, S. Liu, E. Gonzales, J. Dominguez, T. S. Luk, D. N. Neshev, I. Brener, and Y. Kivshar, “Tailoring directional scattering through magnetic and electric resonances in subwavelength silicon nanodisks,” \JournalTitleACS Nano 7, 7824–7832 (2013).
- [60] A. B. Evlyukhin, C. Reinhardt, and B. N. Chichkov, “Multipole light scattering by nonspherical nanoparticles in the discrete dipole approximation,” \JournalTitlePhys. Rev. B 84, 235429 (2011).
- [61] J. van de Groep and A. Polman, “Designing dielectric resonators on substrates: Combining magnetic and electric resonances,” \JournalTitleOpt. Express 21, 26285–26302 (2013).
- [62] M. A. van de Haar, J. van de Groep, B. J. Brenny, and A. Polman, “Controlling magnetic and electric dipole modes in hollow silicon nanocylinders,” \JournalTitleOpt. Express 24, 2047–2064 (2016).
- [63] R. Alaee, R. Filter, D. Lehr, F. Lederer, and C. Rockstuhl, “A generalized kerker condition for highly directive nanoantennas,” \JournalTitleOpt. Lett. 40, 2645–2648 (2015).
- [64] W. Liu and Y. S. Kivshar, “Generalized kerker effects in nanophotonics and meta-optics [invited],” \JournalTitleOpt. Express 26, 13085–13105 (2018).
- [65] A. Kuznetsov, A. Miroshnichenko, Y. H. Fu, J. Zhang, and B. Luk’yanchuk, “Magnetic light,” \JournalTitleScientific reports 2, 492 (2012).
- [66] F. Monticone, D. Sounas, A. Krasnok, and A. Alù, “Can a nonradiating mode be externally excited? nonscattering states versus embedded eigenstates,” \JournalTitleACS Photonics 6, 3108–3114 (2019).
- [67] V. Savinov, N. Papasimakis, D. P. Tsai, and N. Zheludev, “Optical anapoles,” \JournalTitleCommunications Physics 2, 69 (2019).
- [68] A. Krasnok, D. Baranov, H. Li, M.-A. Miri, F. Monticone, and A. Alú, “Anomalies in light scattering,” \JournalTitleAdv. Opt. Photon. 11, 892–951 (2019).
- [69] A. B. Evlyukhin and S. I. Bozhevolnyi, “Resonant unidirectional and elastic scattering of surface plasmon polaritons by high refractive index dielectric nanoparticles,” \JournalTitlePhys. Rev. B 92, 245419 (2015).
- [70] K. Thyagarajan, C. Santschi, P. Langlet, and O. J. F. Martin, “Highly improved fabrication of ag and al nanostructures for uv and nonlinear plasmonics,” \JournalTitleAdvanced Optical Materials 4, 871–876 (2016).
- [71] S. Wack, P. Lunca Popa, N. Adjeroud, J. Guillot, B. R. Pistillo, and R. Leturcq, “Large-scale deposition and growth mechanism of silver nanoparticles by plasma-enhanced atomic layer deposition,” \JournalTitleThe Journal of Physical Chemistry C 123, 27196–27206 (2019).
- [72] M. Mäkelä, T. Hatanpää, K. Mizohata, K. Meinander, J. Niinistö, J. Räisänen, M. Ritala, and M. Leskelä, “Studies on thermal atomic layer deposition of silver thin films,” \JournalTitleChemistry of Materials 29, 2040–2045 (2017).
- [73] T. Das, P. P. Iyer, R. A. DeCrescent, and J. A. Schuller, “Beam engineering for selective and enhanced coupling to multipolar resonances,” \JournalTitlePhys. Rev. B 92, 241110 (2015).
- [74] A. Unger and M. Kreiter, “Analyzing the performance of plasmonic resonators for dielectric sensing,” \JournalTitleThe Journal of Physical Chemistry C 113, 12243–12251 (2009).
- [75] B. Špačková, P. Wrobel, M. Bocková, and J. Homola, “Optical biosensors based on plasmonic nanostructures: A review,” \JournalTitleProceedings of the IEEE 104, 2380–2408 (2016).
- [76] B. Abasahl, C. Santschi, and O. J. F. Martin, “Quantitative extraction of equivalent lumped circuit elements for complex plasmonic nanostructures,” \JournalTitleACS Photonics 1, 403–407 (2014).
- [77] H.-P. Ma, J.-H. Yang, J.-G. Yang, L.-Y. Zhu, W. Huang, G.-J. Yuan, J.-J. Feng, T.-C. Jen, and L. Hong, “Systematic study of the siox film with different stoichiometry by plasma-enhanced atomic layer deposition and its application in siox/sio2 super-lattice,” \JournalTitleNanomaterials 9, 55 (2019).
- [78] R. Adato and H. Altug, “In-situ ultra-sensitive infrared absorption spectroscopy of biomolecule interactions in real time with plasmonic nanoantennas,” \JournalTitleNature communications 4, 2154 (2013).
- [79] A.-P. Blanchard-Dionne and O. J. F. Martin, “Teaching optics to a machine learning network,” \JournalTitleOpt. Lett. 45, 2922–2925 (2020).
4 Supplementary Information: Variation of peaks of ED and MD resonances from Fig. 1
