Electrostatic dipole polarizability and plasmon resonances of multilayer nanoshells

We propose a generalized formula for calculating the dipole polarizability of spherical multilayer nanoshells (MNSs) within the long-wavelength approximation (LWA). Given a MNS with a finite number of concentric layers, radii, and dielectric properties, embedded in a dielectric medium, in the presence of a uniform electric field, we show that its frequency-dependent and complex dipole polarizability can be expressed in terms of the dipole polarizability of the preceding MNS. This approach is different from previous more involved methods where the LWA polarizability of a MNS is usually derived from scattering coefficients. Using both finite-element method- and Mie theory-based simulations, we show that our proposed formula reproduces the usual LWA results, when it is used to predict absorption spectra, by comparing the results to simulated spectra obtained from MNSs with $n$ number of layers up to $n = 6$ layers. An iterative algorithm for calculating the dipole polarizability of a MNS based on the generalized formula is presented. A Fr\"{o}hlich function whose zeroes correspond to the dipolar localized surface plasmon resonances (LSPRs) supported by the MNS is proposed. We identify a pairing behaviour by some LSPRs in the Fr\"{o}hlich function that might also be useful for mode characterization.


Introduction
The optical properties of nanoshells have been studied since the works of Kerker [1], Wu and Yang [2], and others [3][4][5], using Mie theory [6][7][8].The approach presented in this work is a surprisingly simple alternative for obtaining the dipole polarizability of nanoshells using the long-wavelength approximation (LWA) of Maxwell's equations, also known as the quasi-static limit [6].Hence, the formulae proposed herein have the usual assumptions associated with the LWA.In the LWA, we will show that if the electrostatic polarizability is considered as the physical quantity of interest instead of the ubiquitous approach that is based on scattering coefficients [9][10][11][12][13][14][15], a generalized analytical form of the dipole polarizability can be derived for the nth nanoshell given a multilayer nanoshell (MNS) with n number of concentric layers.In addition, a function that predicts the localized surface plasmon resonances (LSPRs) supported by the nth nanoshell, which we will refer to as the Fröhlich function, can be derived using this general form of the dipole polarizability of the nth nanoshell when the Fröhlich condition [7,16,17] is applied to the polarizability.
Several authors have shown that the dipole polarizability of a MNS can be calculated using the LWA -an approximation that leads to the de-coupling of the electric and magnetic field intensities and allows retardation effects such as radiation damping and dynamic depolarization to be ignored, as long as the sizes of the nanoparticles being considered are within the Rayleigh regime [6,10,12,14,15,[18][19][20].This regime considers particle sizes that are small compared to the excitation wavelength, usually within one-tenth of the excitation wavelength or less [6,7,16,17,21].However, at such sub-100 nm particle sizes, usually around 10 nm or less, metal nanoparticles (MNPs) have been shown to display non-local effects [16,17,22].These effects, which are mostly due to increased scattering of the quasi-free electrons near the metal surface, lead to a dependence of the dielectric function of the metal on the longitudinal wavevector of the incident field [16,17].In most studies, this dependence is ignored via the introduction of the local response approximation (LRA) [11,21,23].While the non-local response leads to spectral broadening and suppression of the scattered intensity due to the increased plasmon damping as well as size-induced LSPR shift, it predicts the same number of LSPRs supported by a given MNP when compared to the LRA [16,17,22].
Previously, Daneshfar and Bazyari [9] had proposed a generalized electrostatic dipole polarizability for an n−layer nanoshell that depends on the scattering coefficients of the electric fields in the core, shell, and medium regions of the nanoshell, given the material properties, and sizes of the core and shells, which they used to calculate the spectral properties of a MNS with n = 3 shells.However, Wu and Yang [2] were the first to propose and implement a recursive algorithm for calculating the multipole optical response of an n-layered sphere in terms of the scattering coefficients of the electric and magnetic field intensities in the core, shell, and host regions of the nanoshell, based on Mie theory [8], following the works of Kerker [1].Pal et al. [3][4][5] improved Yang's work and have used it to predict the dipolar extinction efficiencies of a MNS with n = 3 layers [5].The approach by Wu and Yang [2], Pal et al. [3][4][5], and more recent work [33][34][35][36][37] should always be the standard, especially because they account for retardation effects and multipolar response.However, in the LWA, these approaches [9,10,12,14,19]-where the dipole polarizability of the nth nanoshell is obtained from scattering coefficients-are not necessary if the polarizability of the preceding (n − 1)th nanoshell is known, as we will show by reproducing the usual LWA results using the proposed formula.Our model is less complicated to work with, and also very insightful, especially when it is used to predict LSPRs via the proposed Fröhlich function.
Given the electrostatic dipole polarizability of a MNP, its imaginary part is proportional to the absorption efficiency of the MNP [18,25].We will show that the dipole polarizability of a MNS with n shells can be calculated from the dipole polarizability of the preceding (n − 1)th shell and the material properties and size of the nth shell.Here, we will employ both the LWA and the LRA to predict the dipolar LSPRs and the absorption spectra of concentric MNS with n layers or shells up to n = 6 shells.We will show that the formula is more straightforward to work it, and that it reduces to the dipole polarizability of a spherical MNP when there are no shell(s), i.e., when n = 0. We will also show that the implemented iterative algorithm for predicting the dipole polarizability of the MNS reproduces the usual LWA results by comparing our model to Mie theory-based simulations using the open source python software, scattnlay [3][4][5], as well as finite-element method-based, classical electrodynamics simulations in 3D using Comsol Multiphysics ® [38,39].

Electrostatic dipole polarizability
An x−polarized electric field, E = E 0 e −ikz êx , with amplitude, E 0 , is incident on the MNS as shown in Fig. 1.The magnitude of the wavevector is k = 2π √ ϵ n+1 /λ, where ε n+1 is the permittivity of the medium surrounding the nth shell, and λ is the excitation wavelength.Since z = r cos θ, given that r n is the radius of the nth shell, and θ is the polar angle, the electric field in each region of a MNS depends on the size parameter, kr n , as shown in Refs.[3,5].However, in the LWA, this is not the case, since r n << λ, and kr n −→ 0, leading to the following electrostatic dipole polarizability of a sphere (i.e., no shells, n = 0) [6,7]: where Eq. ( 1) is the dipole polarizability normalized by 4πϵ 1 , and ϵ 0 and ϵ 1 are the permittivities of the sphere and the medium surrounding the sphere, respectively, and r 0 is the sphere radius.When there is only one concentric shell (n = 1), the dipole polarizability of the nanoshell in the LWA, derived in Ref. [6], and normalized by 4πϵ 2 , is obtained as where Equation ( 2) can be re-written as: where r 1 is the nanoshell radius, ϵ 1 is the permittivity of the shell, ϵ 2 is the permittivity of the medium surrounding the shell, and When there are two concentric shells (n = 2), the dipole polarizability of the nanoshell in the LWA, as derived in Refs.[10,12,40], and normalized by 4πϵ 3 , is and can be re-written as: where r 2 is the nanoshell radius, ϵ 2 is the permittivity of the shell, ϵ 3 is the permittivity of the medium surrounding the shell, and By inspection, there is a common pattern which α 1 and α 2 follow, and so do f 0 and f 1 .We can utilize these patterns to obtain α 3 and f 2 by induction, for a MNS with three concentric shells (n = 3), as follows: where α 3 is the dipole polarizability of the nanoshell, ϵ 3 is the permittivity of the shell, ϵ 4 is the permittivity of the medium surrounding the shell, and r 3 is the radius of the nanoshell.Hence, by induction, we can re-write Eqs.(9a) and (9b) for the nth shell to obtain: where r n is the radius of the nth nanoshell, ϵ n is the permittivity of the nth shell, ϵ n+1 is the permittivity of the medium surrounding the nth shell, and α n−1 is the normalized (normalized by 4πϵ n ) dipole polarizability of the (n − 1)th shell, and α n is the normalized (normalized by 4πϵ n+1 ) dipole polarizability of the nth shell.An iterative algorithm proposed for calculating α n can be found in Algorithm 1 of the Appendix section.Although substituting n = 0, 1, 2, ..., in Eqs.(10a) and (10b) will reproduce α 0 , α 1 , α 2 , ..., we now need to use these equations to calculate the spectral properties of a given MNS, and compare our results to simulations, in order to ascertain their validity.One such property is absorption.In the LWA, the absorption efficiency of a MNS in the presence of an incident electric field can be calculated using the equation [6,10,12,16]:

The Fröhlich function
Given a MNP with a certain polarizability, the Fröhlich condition [6,7,19] states that the LSPRs supported by the MNP are the poles of the polarizability.In other words, the LSPRs correspond to frequencies where the real part of the denominator of the polarizability vanishes.Let us call this real part of the denominator of the polarizability the Fröhlich function.Starting with Eq. ( 2), the dipole polarizability of a spherical nanoshell with n = 1, the denominator of α 1 simplifies to: We can re-write Eq. ( 12) as: with and so that the Fröhlich function for the n = 1 nanoshell is: For the MNS with n = 2, the denominator of α 2 given in Eq. ( 6) simplifies to: with We can re-write Eqs.(18a) and (18b) in terms of N 0 and D 0 to obtain: so that with , we obtain the Fröhlich function for the n = 2 nanoshell: ) Hence, by induction, the Fröhlich function for the nth shell is: In Eq. ( 21), where, for n ≥ 0, The zeroes of the Fröhlich function for the nth shell, i.e., frequencies (or excitation wavelengths) where F n = 0, correspond to the LSPRs supported by the MNS.An iterative algorithm proposed for calculating F n can be found in Algorithm 2 of the Appendix section.

Simulations
The first set of simulations were based on the finite-element method (FEM).These FEM-based simulations were performed in the Wave Optics module of COMSOL Multiphysics ® software [38] using a spherically-symmetric perfectly-matched layer (PML) and scattering boundary conditions applied in the internal PML surface.An x−polarized incident electric field, as described earlier, was applied to the MNS.The absorption efficiency of the nanoshell is calculated through the following expression [10,39]: The integral in Eq. ( 24) is a volume integral of the total power dissipation density of the nanoshell, P diss , where P 0 is the power density of the incident electric field, and A is the area of the nanoshell obtained from a surface integral over the nanoshell surface [39].
In the second set of simulations, we used scattnlay [5]-an open source code for investigating the scattering of electromagnetic radiation by a multilayered spheredeveloped by Pal et al. [3][4][5] based on Mie theory (MT) [1,6,8].In these MT-based simulations, the absorption efficiency of a given MNS is calculated from the extinction efficiency, Q ext , and the scattering efficiency, Q sca , using the following equation [3]: where [3,4] In Eqs.(26a) and (26b), l = 1, 2, 3, ... is an angular momentum number, denoting the number of multipoles over which the summation is done, a l and b l are scattering coefficients, and x n is the size parameter of the nth layer [3].
As an example, we investigated nanoshells containing only one type of metal as well as one type of dielectric, as illustrated in Fig. 1.To account for both intra-and inter-band electron damping in the metal, we used the Drude-Lorentz model for the permittivity of metals, given in Ref. [41] as : where ε(ω) in the case of Fig. 1, is the permittivity of the metallic core (or shell), ω is the frequency of the incident field, ε ∞ is the high-frequency permittivity of the positive ion core in the metal, f 0 , ω 0 , and γ 0 are the oscillator strength, plasma frequency, and damping rate of the intra-band electrons in the metal, respectively, s is the number of Lorentz oscillators with oscillator strength f s , damping rate γ s , and frequency ω s , associated with inter-band electrons in the metal.The values of these parameters for the most common metals used in plasmonics are given in Ref. [41].

Results and Discussion
Table 1 shows the material compositions of the MNS we investigated using the proposed dipole polarizability formula.Here, we will consider glass, with a permittivity of 2.25 as the dielectric as well as the medium surrounding the outermost nanoshell, and gold, modelled with the Drude-Lorentz parameters given in Ref. [41], as the metal.We Table 1 List of the nanoshell configurations we investigated.The letter "D" stands for dielectric and "M" stands for metal.The starting configurations differ but we have kept the motif uniform.For the nanoshells in Fig. 1(a), the number of layers, n, are odd, i.e., n = 1, 3, 5, while for the nanoshells in Fig. 1 chose to use only two materials for simplicity and ease of analysis of the results but the formula is also valid when different material compositions are used.The radius of the core of the nanoshell in each of the nanoshell configurations in Table 1 (or in Fig. 1) is kept constant at 15 nm.The shells were modelled as equidistant shells, each of thickness, 5 nm.

Absorption spectra
Using the MT-based results as the benchmark-since it is based on exact multipole expansion of the field intensities in the core, shell, and host mediums of the nanoshellthe usual limitations of the LWA, i.e., under-prediction of the absorption efficiencies (and other efficiency factors [16]) and under-prediction of the LSPR shift [16,17] (less redshift), can be seen in Fig. 2.However, in the MNS configurations in Fig. 2, the LWA does not under-predict all the absorption peaks-it appears that one of the absorption peaks is always in good agreement with the MT-based simulations.For instance, as can be seen in Fig. 2, their absorption peaks are comparable at: 693 nm (for n = 1), 644 nm (for n = 3), and 628 nm (for n = 5), for the MNS with a dielectric core, and at: 958 nm (for n = 2), 835 nm (for n = 4), and 784 nm (for n = 6), for the the MNS with a metallic core-under-going a blueshift as the number of shells increases.The FEMbased simulations also follow a similar trend as the LWA except that it over-predicts some of the absorption peaks.Though the most redshifted peaks (Figs.2(e) and 2(f)) are more suppressed in the MT-based results, the number of peaks are always the same when compared to the LWA results, unlike the FEM-based simulations, showing that the LWA model is in good agreement with the MT-based simulations.
On the other hand, when compared to the MT-based simulations, the FEM-based simulations are very accurate in predicting the LSPR shift, as shown by the agreement in the wavelengths corresponding to the absorption peaks in Fig. 2.This is because the FEM-based simulations account for phase retardation effects between oscillating modes of the scattered fields in the MNP [39].We therefore attribute the underprediction of the LSPRs by the LWA to mostly retardation effects.The difference in the absorption peaks of the FEM-based and MT-based simulations most likely stems from the approximation of the spherical Bessel functions in the former while the latter uses the exact functions.It is likely that these approximations accumulate as the number of shells increase causing a noticeable discrepancy between the two simulations.The accuracy of the FEM-based simulations depends largely on the properties of the mesh elements used as well as on the PML features, so that in principle, the simulation can be improved.However, this work is the first to compare comsol multiphysics and scattnlay simulations, using MNSs as the target nanostructure, to the best of our knowledge.
As we will show in Section 3.2 using near-field intensity plots, the LSPRs in Fig. 2 (green curves) all have a dipole character, which means that their origin is most likely to have been dominated by dipolar modes, due to the small particle size of the nanoshells we studied.In addition, Eq. ( 10a), which we used to calculate the absorption spectra of the MNS in Fig. 2 (red curves) is a dipole polarizability.Thus, the LSPRs in the LWA spectra are all due to dipole hybridization between the core and shell plasmons of the MNS.Hence, in the LWA, these LSPRs all have a dipole character.
A noticeable difference between the absorption spectra of the nanoshells with a dielectric core (Figs.2(a), (c), and (e) with n = 1, 3 and 5, respectively) and those with a metallic core (Figs.2(b), (d), and (f), with n = 2, 4 and 6, respectively), is the enhancement versus suppression of the leading absorption peak.This can be attributed to the presence of cavity plasmons in the core (in the case of Figs.2(a), (c), and (e)) and solid plasmons in the core (in the case of Figs.2(b), (d), and (f)).In the MNS in Figs.2(a), (c), and (e), hybridization of the cavity plasmons of the core with the solid plasmon of the innermost shell results in the formation of the leading LSPR, while in the MNS in Figs.2(b), (d), and (f)), the solid plasmons of the core hybridize with the cavity plasmon of the innermost shell to form the leading LSPR.Previous studies have reported similar trends for n = 2 [4,10] and n = 3 [4,9].As we will discuss in the next section, any one of the LSPRs supported by the MNS is either due to a bonding mode, an anti-bonding mode, or a non-bonding mode [42][43][44][45].

Plasmon resonances
The LSPRs supported by the nth nanoshell can be further investigated using graphical solutions of the Fröhlich function given in Eq. ( 21).Electric field enhancement plots produced from simulations in scattnlay [5] reveal the hybridized nature of these LSPRs, in agreement with plasmon hybridization theory (PHT) [42].We will also characterize the LSPRs as either due to a bonding, an anti-bonding, or a non-bonding dipole mode, using the behaviour of the electric field distributions in the core and shell regions of the MNS.
Fig. 3(a) shows the plasmon spectrum of the MNS produced by plotting the Fröhlich function.There are more LSPRs located at shorter wavelengths (between 400 nm and 1200 nm) than at longer wavelengths (between 1200 nm and 2400 nm), which can be attributed to particle size effects.The number of LSPRs increase with increase in the number of metallic shells in the MNS.This is due to an increase in the number of wavelength-dependent permittivity terms in the polarizability.Also, the spacing between the LSPRs increases as the number of shells increases.For n = 1, 2, and n = 4, the number of LSPRs predicted by the Fröhlich function agrees with the absorption spectra (Figs.2(a), (b) and (d), respectively).However, for n = 3, 5, and 6, respectively, the Fröhlich function predicts an extra LSPR in each case compared to the absorption spectra (Figs.2(c), (e), and (f), respectively).This is most likely due to the dipole moment of the extra LSPR being too weak to contribute to the absorption spectra in Fig. 2, making it to appear invisible in the spectra.According to Figs. 3(a) and (b), the extra LSPR is located at a shorter wavelength than the leading LSPR in the absorption spectra in Figs.2(c), (e), and (f), respectively.
The n = 1 nanoshell supports two LSPRs-a bonding dipole (BD) LSPR and an anti-bonding dipole (ABD) LSPR, according to PHT [42,43].The BD LSPR (ω − ) is a longer-wavelength mode formed due to the symmetric coupling between solid sphere plasmons of the shell and cavity sphere plasmons of the core while the ABD LSPR (ω + ) is a shorter-wavelength mode formed as a result of antisymmetric coupling between cavity sphere plasmons of the core and the solid sphere plasmons of the shell.In the BD mode, solid plasmons contribute more than cavity plasmons while the opposite is the case in the ABD mode [43].However, the dipole moment of the ABD mode is very weak, and therefore, it can only be visible in spectral calculations when interband damping is ignored [43,44], which is not the case in this work.Fig. 4(a) shows that the LSPR at 685 nm is the BD mode since the incident electric field experiences the most enhancement outside the nanoshell in contrary to the ABD mode, where the electric field has been shown to be more enhanced inside the nanoshell (see Ref. [44], Fig. 1).Hence, the LSPR in Fig. 3 (red line/curve), for n = 1, is a BD mode.
The n = 2 nanoshell supports three LSPRs -a BD LSPR, an ABD LSPR, and a non-bonding dipole (NBD) LSPR [27,45].The NBD LSPR (ω − + ) is a short-wavelength mode formed by symmetric coupling between the solid sphere plasmon of the core and the ABD mode (ω + ) of the nanoshell.Due to its weak dipole moment, it is only visible in the spectra of silver nanoshells [29], where plasmon damping is significantly lower, or in the absence of inter-band effects.Hence, in our case, the two LSPRs supported by the n = 2 MNS, as shown in Fig. 3, are an ABD mode (ω + − ) and a BD mode (ω − − ).The longer-wavelength, BD (ω − − ) mode, is due to anti-symmetric coupling between the BD (ω − ) nanoshell plasmon and the solid sphere plasmon of the core, while the shorter-wavelength, ABD (ω + − ) mode, is due to symmetric coupling between the BD (ω − ) nanoshell plasmon and the solid sphere plasmon of the core [27,29,45].As shown in Fig. 4(b), for n = 2, the enhanced electric field is situated entirely in the dielectric shell at 952 nm (the BD mode) compared to the electric field distribution at 548 nm (the ABD mode), in agreement with previous work [45].
In Fig. 3(b), a plot of the natural logarithm of the Fröhlich function is shown.Here, the positions of the LSPRs are much more visible, and the plot also reveals an intricate detail about the plasmonic behaviour of the LSPRs-the pairing of certain LSPRs beyond n = 1 due to the convex curves in the Fröhlich function.We propose that this pairing behaviour can be used in determining whether a given LSPR supported by the MNS has a BD or an ABD character.For instance, when n = 2, we have the blueviolet curve in Fig. 3(b) which contains a pair of LSPRs corresponding to an ABD mode at 548 nm and a BD mode at 952 nm, in agreement with Fig. 2(b).However, we will use the rest of this section to show that any given pair of LSPRs in the Fröhlich function does not always consist of an ABD and a BD mode.On the other hand, subsequent analysis of the unpaired LSPRs, shows that a shorter-wavelength LSPR (the leftmost LSPR in the Fröhlich function) that does not participate in the pairing behaviour is always an ABD mode, while a longer-wavelength LSPR (the rightmost LSPR in the Fröhlich function) that does not form a pair has a BD character.These results are summarized in Table 2.
To support the above claim, consider the rest of the electric field intensity plots shown in Figs.4((c)-(f)).The n = 3 MNS supports up to four LSPRs in the absence of inter-band damping in gold, according to PHT [43].These LSPRs are formed as a result of hybridization between the inner-and outer-nanoshell BD and ABD modes.They include: ω − − , a BD mode due to symmetric coupling between the inner-and outer-nanoshell BD modes, ω + − , an ABD mode due to anti-symmetric coupling between the inner-and outer-nanoshell BD modes, ω − + , a NBD mode due to anti-symmetric coupling between the inner-and outer-nanoshell ABD modes, and ω + + , a NBD mode due to anti-symmetric coupling between the inner-and outer-nanoshell ABD modes.The dipole moment of the ω + + mode is very weak, and thus, it can only be visualized with a Drude model for gold.Likewise, the ω − + mode is also a dark mode, since its contribution is not visible in the absorption spectra in Fig. 2(c), but it does appear in the Fröhlich function at 482 nm (green line/curve in Fig. 3).The weak-dipole moment explanation for the absence of ω − + (NBD) mode in Fig. 2(c) is also supported by Fig. 4(c), where the electric field is weakly distributed in the core and inside the second shell at 482 nm.On the other hand, the ω + − and ω − − modes are visible in Fig. 2(c) due to their large dipole moments.This gave rise to the electric field distributions shown in Fig. 4(c), where at 637 nm, the electric field is distributed in both the dielectric core and metallic shells, compared to the field distribution at 1227 nm, where the field is only enhanced inside one of the dielectric shells.Therefore, compared to Fig. 3 (green curve), we have a LSPR pair corresponding to an ABD mode at 482 nm and another ABD mode at 637 nm, as well as a rightmost LSPR corresponding to the BD mode at 1227 nm.
To characterize the rest of the LSPRs in Figs. 2 and 3, we utilized the electric field enhancement plots in Fig. 4((d) -(f)).Using the same analogy as above, for n = 4 (Fig. 4((d)), we assigned the LSPRs at both 823 nm and 1530 nm a BD character, since the enhanced electric field is mostly distributed in the dielectric shells, i.e., the electric field enhancements in the metallic regions are negligible compared to the LSPR at 560 nm, which we assigned an ABD character.Thus, for n = 5, the LSPRs at both 1050 nm and 1830 nm were assigned a BD character while the LSPRs at both 506 nm and 623 nm, where some significant field enhancement exist in the metallic regions, were assigned an ABD character.Similarly for n = 6, the LSPRs at 767 nm, 1295 nm, and 2204 nm were assigned a BD character, while those at both 483 nm and 571 nm were assigned an ABD character.As shown in Table 2, a LSPR pair located in the shortwavelength region of the plasmon spectrum is most likely an ABD-ABD pair, while a pair located in the long-wavelength region consists of a BD-BD pair.However, a LSPR pair in-between the short-and long-wavelength regions is most likely an ABD-BD pair.

Conclusion
The dipole polarizability formula we have proposed for spherical MNSs is straightforward to implement.It is an alternative for obtaining the dipole polarizability of MNSs in the LWA without going through the conventional approach based on scattering coefficients.The formula works well especially when used to predict the LSPRs supported by a MNS via the proposed Fröhlich function.We have shown that the LSPRs correspond to the zeroes of the Fröhlich function for the nth shell, and that a pairing behaviour by some of the LSPRs identified in the Fröhlich function might be useful for mode characterization.However, we maintain that it is an approximate formula strictly valid in the LWA, and should not be used otherwise.The limit of validity of the formula was revealed by comparison with numerical simulations, showing that it reproduces the usual LWA results.Therefore, for the analytical modelling of the optical response of MNSs, the proposed formula is neither meant to replace the standard non-LWA approaches nor to eliminate the limitations of the LWA, but rather to highlight its simplicity and applicability when compared to conventional LWA.

Fig. 1
Fig.12D schematic of the model geometries of a MNS with: (a) 5 layers and (b) 6 layers, in the presence of an x−polarized incident electric field, E, propagating along the z−direction with a wavevector k.The permittivity of the core is ϵ 0 , that of the medium surrounding the: (a) 5th shell is ϵ 6 , and (b) 6th shell is ϵ 7 , and ϵn, n = 1, 2, 3, ..., 6, are the shell permittivities.The radius of the core is r 0 , and rn, n = 1, 2, 3, ..., 6, are the nanoshell radii.As we will explain in detail in Section 3, we only consider two kinds of MNSs-in (a), the core is a dielectric and in (b), the core is a metal.These are the two most common MNS configurations based on material composition but the model is also valid when semiconductors are used.

Fig. 3
Fig. 3 (a) A plot of the Fröhlich function showing the plasmon spectrum of the MNSs.The LSPRs correspond to wavelengths at which Fn(λ) = 0.The Fröhlich function is either an n degree polynomial or an (n − 1) degree polynomial in λ.The plot has been zoomed-in to show only the parts of the polynomials that pass through Fn(λ) = 0. (b) Natural logarithm of the Fröhlich function, ln[Fn(λ)], showing a more detailed version of the plasmon spectrum of the MNSs.The paired and/or unpaired LSPRs are respectively due to the convex and/or concave curves in the Fröhlich function of each MNS.

Fig. 4
Fig. 4 Contour plots of the electric field enhancement, |E|/|E 0 |, inside and outside the MNS, in the xy-plane, using scattnlay [3, 5], where |E| is the magnitude of the incident and scattered fields and |E 0 | is the magnitude of the incident field.The contour plots were obtained at the LSPRs supported by each MNS in the LWA via the Fröhlich function: (a) n = 1, (b) n = 2, (c) n = 3, (d) n = 4, (e) n = 5, and (f) n = 6.These LSPRs are indicated (in nm units) in the core of each nanoshell.

Table 2
List of the nanoshell configurations we investigated and the LSPRs they support, according to the Fröhlich function.Here, " -" stands for a LSPR pair in the plasmon spectrum.