# Geometric interpretations for resonances of plasmonic nanoparticles

## Abstract

The field of plasmonics can be roughly categorized into two branches: surface plasmon polaritons (SPPs) propagating in waveguides and localized surface plasmons (LSPs) supported by scattering particles. Investigations along these two directions usually employ different approaches, resulting in more or less a dogma that the two branches progress almost independently of each other, with few interactions. Here in this work we interpret LSPs from a Bohr model based geometric perspective relying on SPPs, thus establishing a connection between these two sub-fields. Besides the clear explanations of conventional scattering features of plasmonic nanoparticles, based on this geometric model we further demonstrate other anomalous scattering features (higher order modes supported at lower frequencies, and blueshift of the resonance with increasing particle sizes) and multiple electric resonances of the same order supported at different frequencies, which have been revealed to originate from backward SPP modes and multiple dispersion bands supported in the corresponding plasmonic waveguides, respectively. Inspired by this geometric model, it is also shown that, through solely geometric tuning, the absorption of each LSP resonance can be maximized to reach the single channel absorption limit, provided that the scattering and absorption rates are tuned to be equal.

## Introduction

Fueled by the observation of extraordinary transmission through thin metal films [1] and the rapid expansion of the field of metamaterials [2], the old subject of plasmonics has gained a strong renewed impetus and has been experiencing an unprecedented explosive growth [3, 4, 5, 6]. Great progress has been made not only in the branch of propagating SPPs, such as plasmonic circuitry [6], plasmonic nanolaser [7], plasmonic beam shaping [8, 9, 10] and so on, but also in the other branch of LSPs, e.g. superscattering [11, 12], efficient scattering pattern shaping [13, 14, 15] and even LSPs based biological and medical applications [4, 16], to name but a few. However, due to the contrastingly different characteristics (propagating or localized) of the excited states of SPPs and LSPs, and that investigations conducted in these two sub-fields usually employ quite different approaches and techniques, usually SPPs and LSPs are discussed separately [3, 5], which consequently leads to more or less a dogma that the two branches develop almost independently, with few interactions between each other.

Recently there have been some demonstrations relying on the interactions of these two sub-fields, such as the achievements in superscattering [11, 12] and transformation optics [17], nevertheless general and comprehensive studies on the links between them are still not available. Such investigations have been further necessitated by the need of more intuitive understanding of the LSP resonances supported by various plasmonic particles. It is well known that the topic of particle scattering plays a fundamental role in countless related applications [18, 3, 5, 16], however an intuitive understanding of LSPs is not so direct, even for the simplest structures of homogeneous spheres and cylinders. Although Mie theory can give all the information required for the descriptions of LSP resonances of basic spherical or cylindrical plasmonic nanoparticles, the basic physical mechanism is somehow shadowed by the kind of complicated formulas and thus an intuitive and clear picture is still not accessible. Moreover, the topic of light scattering by plasmonic particles is merging rapidly with the current booming fields of graphene [19] and topological states of light and matter [20, 21], making the searching for an intuitive geometric physical model even more urgent.

Here in this paper we provide an intuitive geometric picture for LSPs based on the Bohr model, which requires an integral number of phase accumulation along an enclosed propagating loop to support a well defined localized resonance [22, 23]. The propagating SPPs play a fundamental role in this geometric model, which can explain simply and clearly many conventional scattering features of plasmonic nanoparticles, including well defined LSP resonance in the subwavelength regime, higher (lower) order modes supported at higher (lower) frequencies, and the redshift (blueshift) of the resonances with increasing (decreasing) particle sizes. Based on this geometric model, we further demonstrate many other anomalous scattering features of plasmonic nanoparticles: higher (lower) order modes supported at lower (higher) frequencies and the blueshift (redshift) of the resonances with increasing (decreasing) particle sizes. Those anomalous scattering features have been revealed to originate from the backward SPP modes supported by the corresponding plasmonic waveguides. At the same time, we also show the existence of multiple electric resonances of the same order supported at different frequencies, which have been proved to originate from the multiple dispersion bands of the corresponding waveguides. At the end to further exemplify the effectiveness of the geometric model, we demonstrate that at a fixed resonance frequency, the absorption of a LSP resonance can be maximized to reach the single channel absorption limit through pure geometric tuning, when the scattering and absorption rates are tuned to be equal.

## Results

### Geometric model for Localized surface plasmons

#### The application of Bohr condition to localized surface plasmons

Basically LSPs are localized resonances and are characterized by a set of discrete resonant frequencies [18, 3, 5]. Similar to most localized resonances, the existence of the resonance requires the satisfaction of the Bohr condition [22, 23]. According to the Bohr model, to support a well defined resonance, the length of the enclosed orbit should contain an integral number of the de Broglie wavelengths. It was exactly through Bohr condition that the quantum states of Hydrogen atom were geometrized. The Bohr model can be applied in optics and the Bohr condition can be expressed as:

(1) |

where is the angular wave-number in free space, is position dependent effective refractive index, and m is an integer, which corresponds to the order of the modes supported ( corresponds to dipole, quadrupole, hexapole…). This condition basically means that a resonance requires that the phase accumulation along an enclosed optical path should be an integral number of . One challenge for the application of this model is to decide the loop along which the phase is accumulated. For example, related analysis has been performed for the whispering-gallery modes of dielectric spheres [24]. However, the process to decide the phase accumulation loop is so complicated, especially for the lower order modes, rendering the calculations as complicated as those of rigourous Mie theory and thus losing the simplicity of this model.

In contrast, SPP modes in plasmonic waveguides are propagating coupled states of electrons and photons, with a well defined wave-vector along the propagation direction and with most of the energy confined at the metal-dielectric boundary [3, 5, 6]. Such features of SPP modes make it easier to decide both the effective wave-number and the phase accumulation path [25, 26], and as a result Eq. (1) can then be simplified as:

(2) |

where is the angular wave-number for the SPPs in a corresponding plasmonic waveguide. Similar analysis, such as the Wentzel-Kramers-Brillouin approximation [27], have been performed for plasmonic waveguides in the transverse direction, and explains intuitively the origins of propagating plasmonic modes of different orders excited on the surface of plasmonic nanowires [28, 29].

#### Normal dispersions and conventional scattering features

To clearly explain the links between SPPs and LSPs, and to verify the validity of Eq. (2) for LSPs, we start with the simplest case in Fig.1: single-layered metal nanowire [inset of Fig.1(b)] and its corresponding plasmonic waveguide of semi-infinite metal-dielectric structure [inset of Fig.1(a)]. To perform proof-of-concept demonstrations of the model and at the same time not necessarily lose the generality, throughout this work, we characterize the metal with the Drude model: , where is the angular frequency of light, is the plasma frequency (with corresponding plasmon wavelength of ), is the collision frequency and the loss factor is defined as . The dispersion relation for the SPP mode of the semi-infinite metal-dielectric structure [inset of Fig.1(a)] is: , where and are the permittivities of the metal and dielectric respectively, and is the effective wavenumber of the SPP mode supported. In Fig.1(a) we show the dispersion curve (), where the horizontal axis of angular wave-number is replaced by the phase accumulation along the circumference of the corresponding nanowire in Fig.1(b): . The resonant positions of the first three electric modes predicted by the Bohr model are indicated by points A, B and C, which corresponds to respectively in Eq. (2). For comparison, Fig.1(b) shows the extinction efficiency spectra of the corresponding nanowire () with incident plane wave polarized in plane, to guarantee the excitation of LSP resonances. We note that for simplicity, throughout this work for the dispersion calculation we adopt the lossless Drude model (), and for the scattering calculation we adopt the lossy Drude model (). The positions of extinction peaks for the first three electric modes in Fig.1(b) (marked by by points A’, B’ and C’ respectively) agree well with the points pinpointed in Fig.1(a) (points A, B, C, respectively), verifying the validity of the Bohr model.

For dielectric particles, a well defined resonance can be supported only when the size of the particle is comparable with the effective wavelength of incident waves, in order to provide long enough optical path for sufficient phase accumulation [18], as is indicated by Eq. (1). In sharp contrast to dielectric particles, even deep subwavelength plasmonic nanoparticles are able to support well defined resonances [18, 3, 5]. According to the geometric picture described above, this is due to the fact that the SPP mode can have much larger effective wavenumber and thus sufficient phase can be accumulated over much shorter optical paths. Other well known conventional scattering features of plasmonic nanoparticles are quite similar to those of dielectric particles, including that: (i) higher (lower) order modes are supported at higher (lower) frequencies [as shown in Fig.1(b)]; (ii) the resonances are redshifted (blueshifted) with increasing (decreasing) particle sizes. We exemplify the latter point in Fig.1(c), where the extinction efficiency spectra of the electric dipole modes () are shown for three nanowires of different radii . According to Eq. (2) it is easy to figure out that all those conventional scattering features originate from the normal dispersion of the SPP modes () in the plasmonic waveguides as shown in Fig.1(a).

#### Backward modes and anomalous scattering features

It is well known that the dispersions of SPP modes in plasmonic waveguides can be easily engineered through geometric tuning [3, 5, 30, 31, 32, 33]. Besides the normal dispersion, both flat band () and anomalous dispersion (backward modes, ) can be obtained. According to Eq. (2), unconventional scattering features can be obtained through dispersion engineering of the modes within the corresponding plasmonic waveguides. For example, for a flat band of SPP mode, it is quite natural to expect that modes of different orders can be overlapped at the same frequency and the resonant frequency is fixed and independent of the particle size. Such unusual scattering features have already been observed, with a noticeable example of the recent demonstration of superscattering [11, 12].

Beyond the flat-band dispersion of SPP modes, even backward modes can be obtained in plasmonic structures when the direction of overall energy flow is contradirectional to that of the wave-vector [30, 31, 32]. The inset of Fig.2(a) shows the three-layered plasmonic waveguide we investigate. The green layer denotes dielectric layer, here for which we assign refractive index of (e.g. GaAs, Si or Ge). The bottom metal layer is supposed to be semi-infinite and the background is filled with air of . The corresponding scattering nanowire of this three-layered waveguides is shown in the inset of Fig.2(b) [scattering configuration is the same as that in Fig.1(b)] and the size parameters satisfy: and . Here we define normalized size parameters for both the waveguide and the nanowire: , . For the waveguide, as there are two metal-dielectric boundaries, we approximate the effective length of the phase accumulation path as . The specific size parameter for the waveguide is and for the nanowire is . The dispersion of backward SPP mode (TM modes) supported in the waveguide is shown in Fig.2(a), with two points A and B pinpointed, where the Bohr condition is satisfied with respectively. The extinction efficiency spectra of the corresponding nanowire is shown in Fig.2(b), with contributions from the first two electric modes shown. As expected, the positions of the spectra peaks agree well with the positions of pinpointed points in Fig.2(a). In sharp contrast to the results in Fig.1(b), here the higher order mode (, quadruple mode) is supported at a lower frequency, and lower order mode (, dipole mode) is supported at a higher frequency. It is obvious that such anomalous scattering comes from the anomalous dispersion shown in Fig.2(a), where the backward mode plays the fundamental role. We further investigate the size dependence of the resonances originating from backward modes, and show the results in Fig.2(c). The extinction efficiency spectra of the dipole modes () supported by three multi-layered nanowires with different size parameter are shown with: , (r=[0.10, 0.17, 0.23]) and (r=[0.20, 0.27, 0.33]). It is clear that the resonances are blueshifted (redshifted) with increasing (decreasing) particle sizes, which is opposite to the size-scaling trend shown in Fig.1(c). This anomalous feature also has its origin in the existence of the backward mode. We note here that although similar unusual scattering features mentioned above have been partly observed before [34, 35], but unfortunately the physical mechanism behind had not been clearly revealed or systematically investigated.

#### Multiple dispersion bands induced multiple modes of the same order

Up to now, we have discussed only the LSP resonances based on SPP modes of a single dispersion band. According to Eq. (2), when several dispersion bands are simultaneously involved, different sets of modes of the same order can be supported. The number of the modes of the same order is decided by the number of the dispersion bands. Fig.3(a) shows the dispersion curves of the three-layered metal-dielectric waveguide of the size parameter of , where two dispersion bands are supported. Points A and B predict the spectral position of the localized modes of the corresponding three-layered nanowire of the the size parameter of . At both points the Bohr condition is satisfied with the same integer of , indicating that two electric dipoles should be supported at different frequencies. To verify the prediction shown in Fig.3(a), in Fig.3(b) we show the extinction efficiency spectra of the dipole modes supported by the corresponding nanowire. As predicted, it is clear that two dipoles are supported at different frequencies: dipoles I and II correspond to the points A and B in Fig.3(a) respectively. To further clarify the different origins of the two dipole modes, we show the near-field distributions [out-of-plane magnetic field intensity, insets of Fig.3(b)] of the two modes at the resonant frequencies of and respectively. It is clear that though the two modes are of the same order () and are both electric dipole modes, they show quite different near-field distributions. This is due to the fact that the two dipole modes originate from different SPP modes of different dispersion bands [see Fig.3(a)], which themselves correspond to different near-field distributions within the three-layered waveguides. Multiple electric modes of the same order play a fundamental role in the recent interpretation of scattering dark states, but unfortunately the physical origin had not been revealed [36]. We emphasize that some other models that deal with the LSP modes of plasmonic nanoparticles, such as the the hybridization model put forward in Ref.[37], can explain the existence of multiple modes of the same order, but in comparison the Bohr model presented in this paper is more direct and simpler. We also note that in Fig.3 we discuss only the dipole modes and this principle can be easily applied to other higher order modes.

Based on the results shown in Fig.3, we could comment further on our approximation of the effective radius of the phase accumulation path: , which actually assumes that the center of field distribution is in the middle of the two metal-dielectric interfaces. However, according to the near-field distributions of the two electric dipoles shown in Fig.3(b), at points A and B the approximation underestimates and overestimates the effective length of the phase accumulation path respectively. This is because for dipole I (II), most of the fields are confined closer to the outer (inner) metal-dielectric interface. Consequently points A () and B () over-predicts and under-predicts the actual resonant frequencies of the two dipoles respectively (which are actually at and ).

It is worth mentioning that up to now we have restricted ourselves to two-dimensional cylindrical structures. The Bohr model we have employed is very fundamental and the extension of this model to three-dimensional spherical structures is direct [12]. Specifically, Eq. (1)-Eq. (2) can be directly applied: is the same for cylindrical and spherical structures as they share the same corresponding waveguide; also we can employ the same approximation for the effective phase accumulation path of cylindrical and spherical structures. As a simple demonstration that the Bohr model is also applicable to three-dimensional particles, in Fig.3(c) we show the extinction efficiency spectra of the dipole modes supported by the metal-dielectric-metal nanosphere []. The scattering configuration and the definition of the geometric parameters are the same as that of the nanowire (see Fig.1 and Fig.2). As expected, two dipole modes are observed at different spectral positions, and the physical origin of the two modes is the same as that of the nanowire as we have discussed above.

### Absorption maximization for localized surface plasmons

From the discussions presented above, we know that the Bohr model offers us a new perspective for the understanding of LSP resonances based on SPP modes. As a specific application of this model, we investigate the absorption maximization of LSP resonances based on those new understandings obtained. The ideal case of loss maximization is the recent demonstration of coherent perfect absorption [38, 39, 40]. However, to achieve coherent perfect absorption for scattering particles, there are harsh restrictions on both the shape of the incident waves and the particles themselves [38, 39, 40]. It is hard to meet all those requirements simultaneously and thus more practical to discuss how to maximize the absorption of particles with plane wave incidence. This is certainly one of the most general problems, which might find applications in various fields and applications when the incident waves can be approximated as pane waves.

#### Single-channel absorption limit

We still restrict our discussions to the simple structure of multi-layered nanowires. For the two-dimensional cylindrical structures with plane wave incidence, the absorption cross section of a single resonance channel can be expressed as [41, 11]:

(3) |

where and are the wavelength and frequency of the incident wave respectively; is the resonant frequency; and are the intrinsic loss rate (due to Ohmic loss) and scattering loss rate (due to free-space scattering) respectively. At the resonant frequency, the expression can be simplified as;

(4) |

where is the loss ratio. Obviously the absorption cross sections can be maximized as when , which requires that the intrinsic loss and scattering loss rates are equal [42, 43, 44, 45]. is basically the single channel absorption limit. However, when the wavelength of the incident wave is fixed, to tune the loss ratio to in order to reach the absorption limit without shifting the resonant wavelength is usually very challenging, which could involve complicated material engineering such as doping to change the intrinsic loss rate.

#### Absorption maximization through pure geometric tuning

According to the geometric picture of the LSP resonances we have discussed, a fixed resonant frequency only requires that an integral number of phase has been accumulated at this frequency, which can be satisfied for different phase accumulation paths with different dispersions of the SPP modes. Consequently the loss ratio could still be tuned through pure geometric tuning, without even shifting the resonant frequency. We note here that in this paper our discussions have been restricted to loss maximization of an individual resonance channel, and the losses can be further enhanced through overlapping resonances of different orders [46, 45].

We employ again the three-layered nanowire shown in the inset of Fig.2(b). This time we substitute the dielectric layer with the dielectric layer. To decrease the refractive index of the dielectric layer will make smaller the effective wave-vector mismatch inside and outside the three-layered nanowire, and thus can increase the scattering loss rate [35], rendering it comparable to . To make the two loss rates equal to each other requires extra geometric tuning of the structure. Instead of calculating separately and to decide the loss ratio of , we calculate firstly the scattering and absorption cross sections of the nanowire through Mie theory and then the loss ratio can be obtained directly as the ratio of the absorption cross section to the scattering cross section [41, 11]:

(5) |

This is because that similar to Eq. (3), the scattering cross section can be expressed as: [41, 11]. Fig.4(a) shows both the scattering efficiency (red curves) and absorption efficiency (blue curves) of the dipole modes () supported by two three-layered nanowires of different size parameters: (dashed curves) and (solid curves). When the size parameters are tuned, it is clear from Fig.4(a) that: (1) the resonance position of the absorption can still be fixed () and (ii) as the ratio of the absorption cross section to the scattering cross section, the loss ratio can be effectively tuned. Fig.4(b) shows the change of absorption cross section (normalized by the single channel absorption limit at ) of the electric dipole modes () at the resonant frequency of with different loss ratios . The solid curve corresponds to Eq. (4) and the crosses correspond to the calculated results of nanowires (through Mie theory) with different size parameters. Three points A-C are pinpointed with the corresponding size parameters specified. Further details about absorption cross section spectra at those three points are shown in Fig.4(c). It is clear that the absorption can be effectively maximized through pure geometric tuning, which reaches the single channel absorption limit at the size parameter of . We note here that there are some discrepancies between the results from Eq. (4) and those from Mie theory [Fig.4(a)]. This is due to that the central resonant positions of the absorption and scattering spectra do not fully overlap [Fig.4(a)], and consequently the relation in Eq. (5) is not exactly rigorous. Nevertheless this approximation does not at all obscure the physical mechanism and the evolution trend of absorption with changing loss ratios has been clearly revealed, as shown in Fig.4(b).

## Discussion

In this work, based on the Bohr model, we give a geometric picture for the LSP resonances supported by plasmonic nanoparticles relying on the SPP modes propagating in the corresponding plasmonic waveguides. This geometric model actually establishes a connection between the two sub-fields of plasmonics, and can be firstly applied to explain directly the well known scattering features of plasmonic nanoparticles, such as well defined LSP resonance in the subwavelength regime, higher (lower) order modes supported at higher (lower) frequencies and redshift (blueshift) of the resonances with increasing (decreasing) particle sizes. As a next step, through this geometric model, it is further demonstrated that other anomalous scattering features can be obtained within geometrically engineered multi-layered plasmonic nanoparticles, including higher (lower) order modes supported at lower (higher) frequencies and blueshift (redshift) of the resonances with increasing (decreasing) particle sizes. We reveal that those anomalous scattering features originate from the existence of backward SPP modes in the corresponding waveguides. At the same time, it is shown that the multi-layered nanowires can be engineered to support multiple electric resonances of the same order at different frequencies with different near-field distributions. This unusual feature is induced by the multiple dispersion bands of the corresponding plasmonic waveguides. The number of the modes of the same order is decided by the number of the dispersion bands. At the end we show that at a fixed resonant frequency, the absorption of a plasmonic resonance can be effectively maximized through pure geometric tuning to reach the single channel absorption limit, when the scattering and absorption rates are made equal.

The intuitive geometric picture should not be confined to plasmonic particles and can be actually applied to many other localized resonances in other structures. Despite the simplicity and generality of this model, however there do exist some challenges with it: (i) As is shown in the discussions associated with Fig.3, it is not so simple or direct to rigorously decide the effective length of the phase accumulation path. This challenge becomes even tougher for structures with sharp corners. There have been some related attempts [24, 47], which however have neither the simplify nor clarity. Consequently they show no superiority to the Mie theory or direct simulations, thus offering very limited physical insights. (ii) A lot of work has to be done to make this model applicable and compatible with structures which incorporate materials with nonlocal, nonlinear and/or quantum effects [48, 49, 50, 51]. (iii) Our geometric model applies only to electric plasmonic resonances, and all the resonances discussed in this paper are electric resonances. The geometric model for the optically-induced magnetic resonances [52, 13] is still not available. (iv) The Bohr model interpretation would lead to the conclusion that for all plasmonic particles with sizes getting smaller and smaller into the quasi-static spectral regime, the resonance frequency of all electric resonances would converge to the surface plasmon frequency (). This is exactly the case for nanowire scattering configuration we have investigated in our paper, but not the case for other three-dimensional structures (e.g., for spherical plasmonic nanoparticles, the electric dipole resonant frequency would converge to and for electric resonances of other orders it would converge to higher frequencies). But it is worth noticing that the Bohr model works well for spherical structures with particle sizes not far smaller than the wavelength [see Fig.3(c)]. At the end we would like to note here that in this paper we characterize resonances through far-field scattering spectra and we should keep in mind that resonances can be also described by near-field properties [53, 54].

The intuitive Bohr model based geometric picture we discussed in this work offers new insights into the understanding of LSP resonances and other localized resonances. It might shed new light to investigations related to the particle scattering problem and moreover can possibly accelerate the merging of this field with other rapidly developing fields, such as graphene [19] and topological states of light and matter [20, 21]. This could bring out a lot of extra opportunities, for not only fundamental studies, but also for a lot of particle scattering related applications, including biosensing, nanoantennas, photovoltaic devices, switching and even medical treatments.

## Methods

### Dispersions of the three-layered metal-dielectric plasmonic waveguides

The dispersion relation of the three-layered metal-dielectric plasmonic waveguide [shown in the inset of Fig.2(a)] can be obtained through solving the scalar Helmholtz equation combined with the continuity equations for the fields at the boundaries. As we are interested only in the SPP modes and thus only the TM modes have been discussed. The dispersion relation for TM modes within this waveguide, which links the free-space frequency and the SPP wave-vector , can be expressed as [55]:

(6) |

where ; ; is the hyperbolic tangent function; ; ; denotes the relative permittivity; is the wave-vector in vacuum; subscripts , , denote the metal layer, dielectric layer, and background media respectively. The dispersion curve can be obtained through solving the above equation.

### Scattering of the plasmonic nanowires and nanospheres

For the three-layered plasmonic nanowire (and nanosphere) shown in the inset of Fig.2(b), we obtain its scattering properties through applying the generalized Mie theory [18]:

(1) For nanowires, the incident plane wave is polarized in plane (perpendicular to the nanowire axis). The expressions for the cross sections of extinction, scattering and absorption are [18]:

(7) | |||

(8) | |||

(9) |

where means the real part; is the wavelength in the background material (vacuum in this study); and are the scattering coefficients, whose expressions are long and tedious, which will not be shown here but can be found in Refs. [18]. More specifically and correspond to the optically-induced magnetic dipole and electric modes of order , respectively [15].

(2) For nanospheres, the cross sections are independent of the direction of the polarization:

(10) | |||

(11) | |||

(12) |

where and are the scattering coefficients, whose expressions can be found in Refs. [18]. More specifically and correspond to the electric and optically-induced magnetic modes of order , respectively [13, 14, 56].

It is obvious that the contributions from modes of specific orders can be calculated directly and separately. Our discussion in this paper are restricted to electric modes only. The corresponding efficiencies are defined as: the ratios of the cross sections to the diameter of the outmost layer of the nanowire (for nanowires) and the ratios of the cross sections to the area of the outmost layer of the nanosphere (for nanospheres).

## Acknowledgements

We are indebted to A. E. Miroshnichenko, D. N. Neshev, and J. Hou for many useful discussions and suggestions. This work has been supported by the National Natural Science Foundation of China (Grant No.: ), the Australian Research Council through the Center of Excellence CUDOS, the Leverhulme Trust of UK, and the Basic Research Scheme of College of Optoelectronic Science and Engineering, National University of Defense Technology.

## Author contributions statement

W. L. conceived the idea and performed the calculations; W. L., R. F. O. and Y. S. K. discussed the results and analysed the data; W. L. wrote the paper with a lot of input from R. F. O. and Y. S. K..

## Additional information

Competing financial interests: The authors declare no competing financial interests.

### References

- Ebbesen, T. W., Lezec, H. J., Ghaemi, H. F., Thio, T. & Wolff, P. A. Extraordinary optical transmission through sub-wavelength hole arrays. Nature 391, 667 (1998).
- Zheludev, N. I. & Kivshar, Y. S. From metamaterials to metadevices. Nat. Mater. 11, 917 (2012).
- Zayats, A. V., Smolyaninov, I. I. & Maradudin, A. A. Nano-optics of surface plasmon polaritons. Phys. Rep. 408, 131 (2005).
- Atwater, H. A. The promise of plasmonics. Sci. Am. 296, 56 (2007).
- Maier, S. A. Plasmonics : fundamentals and applications (Springer, New York, 2007).
- Gramotnev, D. K. & Bozhevolnyi, S. I. Plasmonics beyond the diffraction limit. Nat. Photon. 4, 83 (2010).
- Oulton, R. F. et al. Plasmon lasers at deep subwavelength scale. Nature 461, 629–632 (2009).
- Liu, Y. M., Zentgraf, T., Bartal, G. & Zhang, X. Transformational plasmon optics. Nano Lett. 10, 1991 (2010).
- Huidobro, P. A., Nesterov, M. L., Martin-Moreno, L. & Garcia-Vidal, F. J. Transformation optics for plasmonics. Nano Lett. 10, 1985 (2010).
- Liu, W., Neshev, D. N., Miroshnichenko, A. E., Shadrivov, I. V. & Kivshar, Y. S. Polychromatic nanofocusing of surface plasmon polaritons. Phys. Rev. B 83, 073404 (2011).
- Ruan, Z. C. & Fan, S. H. Superscattering of light from subwavelength nanostructures. Phys. Rev. Lett. 105, 013901 (2010).
- Ruan, Z. C. & Fan, S. H. Design of subwavelength superscattering nanospheres. Appl. Phys. Lett. 98, 043101 (2011).
- Liu, W., Miroshnichenko, A. E. & Kivshar, Y. S. Control of light scattering by nanoparticles with optically-induced magnetic responses. Chin. Phys. B 23, 047806 (2014).
- Liu, W., Miroshnichenko, A. E., Neshev, D. N. & Kivshar, Y. S. Broadband unidirectional scattering by magneto-electric core-shell nanoparticles. ACS Nano 6, 5489 (2012).
- Liu, W. et al. Scattering of core-shell nanowires with the interference of electric and magnetic resonances. Opt. Lett. 38, 2621 (2013).
- Huschka, R. et al. Light-induced release of dna from gold nanoparticles: Nanoshells and nanorods. J. Am. Chem. Soc. 133, 12247 (2011).
- Pendry, J. B., Aubry, A., Smith, D. R. & Maier, S. A. Transformation optics and subwavelength control of light. Science 337, 549 (2012).
- Kerker, M. The scattering of light, and other electromagnetic radiation (Academic Press, New York, 1969).
- Geim, A. K. & Novoselov, K. S. The rise of graphene. Nat. Mater. 6, 183 (2007).
- Yin, X. B., Ye, Z. L., Rho, J., Wang, Y. & Zhang, X. Photonic spin hall effect at metasurfaces. Science 339, 1405 (2013).
- Khanikaev, A. B. et al. Photonic topological insulators. Nat. Mater. 12, 233 (2013).
- Landau, L. D. & Lifshits, E. M. Quantum mechanics : non-relativistic theory (Pergamon Press, Oxford; New York, 1977), 3d edn.
- Yang, C. N. Quantum numbers, chern classes, and a bodhisattva. Phys. Today 65, 33 (2012).
- Roll, G. & Schweiger, G. Geometrical optics model of mie resonances. J. Opt. Soc. Am. A 17, 1301 (2000).
- Beck, F. J., Verhagen, E., Mokkapati, S., Polman, A. & Catchpole, K. R. Resonant spp modes supported by discrete metal nanoparticles on high-index substrates. Opt. Express 19, A146 (2011).
- Hasan, S. B. et al. Relating localized nanoparticle resonances to an associated antenna problem. Phys. Rev. B 84, 195405 (2011).
- Marcuse, D. Light transmission optics (Van Nostrand Reinhold, New York, 1972).
- Schmidt, M. A. & Russell, P. S. J. Long-range spiralling surface plasmon modes on metallic nanowires. Opt. Express 16, 13617 (2008).
- Catrysse, P. B. & Fan, S. H. Understanding the dispersion of coaxial plasmonic structures through a connection with the planar metal-insulator-metal geometry. Appl. Phys. Lett. 94, 231111 (2009).
- Shin, H. & Fan, S. All-angle negative refraction for surface plasmon waves using a metal-dielectric-metal structure. Phys. Rev. Lett. 96, 073907 (2006).
- Lezec, H. J., Dionne, J. A. & Atwater, H. A. Negative refraction at visible frequencies. Science 316, 430 (2007).
- Dionne, J. A., Verhagen, E., Polman, A. & Atwater, H. A. Are negative index materials achievable with surface plasmon waveguides? a case study of three plasmonic geometries. Opt. Express 16, 19001 (2008).
- Archambault, A., Besbes, M. & Greffet, J.-J. Superlens in the time domain. Phys. Rev. Lett. 109, 097405 (2012).
- Peng, S., McMahon, J. M., Schatz, G. C., Gray, S. K. & Sun, Y. Reversing the size-dependence of surface plasmon resonances. Proc. Natl. Acad. Sci. U.S.A. 107, 14530 (2010).
- Yang, X., Yao, J., Rho, J., Yin, X. & Zhang, X. Experimental realization of three-dimensional indefinite cavities at the nanoscale with anomalous scaling laws. Nat. Photon. 6, 450 (2012).
- Hsu, C. W., DeLacy, B. G., Johnson, S. G., Joannopoulos, J. D. & Soljacic, M. Theoretical criteria for scattering dark states in nanostructured particles. Nano Lett. 14, 2783 (2014).
- Prodan, E., Radloff, C., Halas, N. J. & Nordlander, P. A hybridization model for the plasmon response of complex nanostructures. Science 302, 419–422 (2003).
- Chong, Y. D., Ge, L., Cao, H. & Stone, A. D. Coherent perfect absorbers: Time-reversed lasers. Phys. Rev. Lett. 105, 053901 (2010).
- Wan, W. et al. Time-reversed lasing and interferometric control of absorption. Science 331, 889 (2011).
- Noh, H., Chong, Y., Stone, A. D. & Cao, H. Perfect coupling of light to surface plasmons by coherent absorption. Phys. Rev. Lett. 108, 186805 (2012).
- Hamam, R. E., Karalis, A., Joannopoulos, J. D. & Soljacic, M. Coupled-mode theory for general free-space resonant scattering of waves. Phys. Rev. A 75, 053801 (2007).
- Seok, T. J. et al. Radiation engineering of optical antennas for maximum field enhancement. Nano Lett. 11, 2606 (2011).
- Fleury, R., Soric, J. & Alu, A. Physical bounds on absorption and scattering for cloaked sensors. Phys. Rev. B 89, 045122 (2014).
- Tretyakov, S. Maximizing absorption and scattering by dipole particles. Plasmonics 9, 935 (2014).
- Mohammadi Estakhri, N. & Alu, A. Minimum-scattering superabsorbers. Phys. Rev. B 89, 121416 (2014).
- Hasegawa, K., Rohde, C. & Deutsch, M. Enhanced surface-plasmon resonance absorption in metal-dielectric-metal layered microspheres. Opt. Lett. 31, 1136 (2006).
- Papoff, F. & Hourahine, B. Geometrical mie theory for resonances in nanoparticles of any shape. Opt. Express 19, 21432 (2011).
- Boardman, A. D. Electromagnetic surface modes (Wiley, Chichester; New York, 1982).
- Scholl, J. A., Koh, A. L. & Dionne, J. A. Quantum plasmon resonances of individual metallic nanoparticles. Nature 483, 421 (2012).
- Kauranen, M. & Zayats, A. V. Nonlinear plasmonics. Nat. Photon. 6, 737 (2012).
- Esteban, R., Borisov, A. G., Nordlander, P. & Aizpurua, J. Bridging quantum and classical plasmonics with a quantum-corrected model. Nat. Commun. 3, 825 (2012).
- Kuznetsov, A. I., Miroshnichenko, A. E., Fu, Y. H., Zhang, J. B. & Lukyanchuk, B. S. Magnetic light. Sci. Rep. 2, 492 (2012).
- Zuloaga, J. & Nordlander, P. On the energy shift between near-field and far-field peak intensities in localized plasmon systems. Nano Lett. 11, 1280 (2011).
- Alonso-Gonzalez, P. et al. Experimental verification of the spectral shift between near- and far-field peak intensities of plasmonic infrared nanoantennas. Phys. Rev. Lett. 110, 203902 (2013).
- Kudo, M. & Mushiake, Y. Mode selecting characteristics of four-layer dielectric slab waveguide. Radio Sci. 17, 125 (1982).
- Liu, W., Miroshnichenko, A. E., Neshev, D. N. & Kivshar, Y. S. Polarization-independent fano resonances in arrays of core-shell nanoparticles. Phys. Rev. B 86, 081407(R) (2012).