OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 19, Iss. 22 — Oct. 24, 2011
  • pp: 21432–21444
« Show journal navigation

Geometrical Mie theory for resonances in nanoparticles of any shape

F. Papoff and B. Hourahine  »View Author Affiliations


Optics Express, Vol. 19, Issue 22, pp. 21432-21444 (2011)
http://dx.doi.org/10.1364/OE.19.021432


View Full Text Article

Acrobat PDF (7099 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We give a geometrical theory of resonances in Maxwell’s equations that generalizes the Mie formulae for spheres to all scattering channels of any dielectric or metallic particle without sharp edges. We show that the electromagnetic response of a particle is given by a set of modes of internal and scattered fields that are coupled pairwise on the surface of the particle and reveal that resonances in nanoparticles and excess noise in macroscopic cavities have the same origin. We give examples of two types of optical resonances: those in which a single pair of internal and scattered modes become strongly aligned in the sense defined in this paper, and those resulting from constructive interference of many pairs of weakly aligned modes, an effect relevant for sensing. This approach calculates resonances for every significant mode of particles, demonstrating that modes can be either bright or dark depending on the incident field. Using this extra mode information we then outline how excitation can be optimized. Finally, we apply this theory to gold particles with shapes often used in experiments, demonstrating effects including a Fano-like resonance.

© 2011 OSA

1. Introduction

The interaction of light with wavelength sized particles has been intensely investigated for more than a century, continuing to provide interesting and surprising results. This coupling is essential in single molecule spectroscopy and single photon processes [1

1. D. Graham and R. Goodacre, “Chemical and bioanalytical applications of surface enhanced Raman scattering spectroscopy,” Chem. Soc. Rev. 37, 883–884 (2008). [CrossRef] [PubMed]

, 2

2. T. Aoki, B. Dayan, E. Wilcut, W. P. Bowen, A. S. Parkins, T. J. Kippenberg, K. J. Vahala, and H. J. Kimble, “Observation of strong coupling between one atom and a monolithic microresonator,” Nature 443, 671–674 (2006). [CrossRef] [PubMed]

], while interference between different scattering channels of particles produces classical analogues of quantum processes [3

3. B. Lukyanchuk, N. Zheludev, S. Maier, N. Halas, P. Nordlander, H. Giessen, and C. Tow Chong, “The Fano resonance in plasmonic nanostructures and metamaterials,” Nat. Mater. 9, 707–715 (2010). [CrossRef]

, 4

4. N. Liu, L. Langguth, T. Weiss, J. Kästel, M. Fleischhauer, T. Pfau, and H. Giessen, “Plasmonic analogue of electromagnetic induced transparency at the Drude damping limit,” Nat. Mater. 8, 758–762 (2009). [CrossRef] [PubMed]

], and carefully designed particles underpin the physical realization of metamaterials [5

5. J. Schuller, E. Barnard, W. Cai, Y. C. Jun, J. White, and M. Brongersma, “Plasmonics for extreme light concentration and manipulation,” Nat. Mater. 9, 193–204 (2010). [CrossRef] [PubMed]

7

7. J. Pendry, D. Schuring, and D. Smith, “Controlling electromagnetic fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

]. The basis all of these effects is that the particle-light interaction, which depends on the composition and shape of the particle and on the properties the incident light, can become very strong around resonances. For particles much larger than the wavelength of light, resonances are described by closed orbits of light rays [8

8. G. Roll and G. Schweiger, “Geometrical optics model of Mie resonances,” J. Opt. Soc. Am. A 17, 1301–1311 (2000). [CrossRef]

] inside the particle. This geometric approach becomes less and less effective as the size of the particle decreases, eventually requiring the solution of Maxwell’s equations. Mie-type solutions [9

9. G. Mie, “Beiträge zur optik trüber medien, speziell kolloidaler metallösungen,” Ann. Phys. 330, 377–445 (1908). [CrossRef]

, 10

10. Y. Han and Z. Wu, “Scattering of a spheroidal particle illuminated by a gaussian beam,” Appl. Opt. 40, 2501–2509 (2001). [CrossRef]

], based on symmetry and coordinate separability, provide analytical description of resonances for a few specific shapes of particle. For spheres, internal and scattered fields are expanded by electric and magnetic multipoles and each multipole in one field is coupled exclusively to the corresponding multipole in the other field. Resonances are independent of the incident field and occur when the coefficients of one pair of multipoles reaches a maximum for the particular values of the particle radius, permittivity and susceptibility.

For particles of arbitrary shape, resonances are always (to the best of our knowledge) defined implicitly by maxima in properties such as the far field extinction or scattering efficiency spectra. This is because, unlike the sphere, there is no decomposition of the internal and scattered fields into partial waves that are pairwise coupled. Hence monitoring of calculated specific properties of the scattered field are instead used empirically to define resonances. While there are several methods that can find spectra and their maxima [11

11. M. I. Mishchenko, J. H. Hovernier, and L. D. Travis, eds., Light scattering by nonspherical particles: Theory, Measurements and Applications (Academic Press, 2000).

], this approach to resonances is unsatisfying because it depends on the incident field; further it fails to recognize the consequences of resonances associated with fields which are strong only in one region of space, and finally requires an a priori choice of the property which being monitored to determine the resonance. For example a resonance associated with a strong surface field, that itself is not efficient at transporting energy to infinity, would not appear as an obvious feature in any far field efficiency spectra used to define resonance; nevertheless such resonances can be extremely important in near field applications or through interference with other channels which themselves are able to transport energy into the far field.

Here we introduce a hermitian operator to define field expansions and resonances for any particle, where Mie’s treatment of the sphere is a special case of this more general theory. We use the mathematical framework of the angles between subspaces of functions [12

12. C. Jordan, “Essai sur la géométrie à n dimension,” Bul. Soc. Math. France 3, 103–174 (1875).

, 13

13. A. Knyazev, A. Jujusnashvili, and M. Argentati, “Angles between Infinite Dimensional Subspaces with Applications to the Rayleigh-Ritz and Alternating Projectors Methods,” J. Func. Anal. 259, 1323–1345 (2010). [CrossRef]

] to reveal the geometrical nature of resonances in Maxwell’s equations. This allows us to predict resonances of particles independently of the particular incident field, obtain the near and far features of particular modes and choose incident fields to optimize their excitation.

2. Theory

The theory we develop applies to general metallic and dielectric particles without sharp edges in cases where the interaction between light and matter inside the particle is described by a local macroscopic permittivity and susceptibility. For metallic particles, this means that the free propagation length of carriers is smaller than the skin depth and all the characteristic lengths of the particle. The tangential components of electric and magnetic fields E,H are continuous on passing through the particle boundaries, and the energy scattered by a particle flows towards infinity. The interaction of the particle with an incident field is determined by finding appropriate solutions of the Maxwell’s equations in the internal and the external media that satisfy these boundary conditions. We use [14

14. K. Holms, B. Hourahine, and F. Papoff, “Calculation of internal and scattered fields of axisymmetric nanoparticles at any point in space,” J. Opt. A, Pure Appl. Opt. 11, 054009 (2009). [CrossRef]

] six component vectors F = [E,H]T for electromagnetic fields. Their projections, f, onto the boundary of the particle are surface fields each with four components, two electric and two magnetic, that form a space ℋ where scalar products are defined in terms of overlap integrals on the surface of the particle, fg=Sfj*gjds, where the index j labels the components, f* is the complex conjugate of f and we sum over repeated indexes. In this formalism the boundary conditions become
f0=fifs,
(1)
which has a simple geometrical meaning in ℋ : the projection, f0, of the incident field, F0(x), onto the surface is equal to the difference between the projections of the internal and scattered fields, fi and fs. This suggests that an incident field with small tangent components can excite large internal and scattered surface fields provided that these two fields closely match. This happens when the “angle” between these two fields, and therefore their difference, is small. The angles in question can be rigorously defined as the angles between standing and outgoing waves that are solutions of the Maxwell’s equations for the internal and external media respectively, and that form two subspaces of ℋ. For each particle, these angles and the associated waves characterize completely the particle’s electromagnetic response, which can be determined with arbitrary precision from any complete set of solutions of the Maxwell equations for the internal and external media.

There exist several sets of exact solutions of the Maxwell equations that are linearly independent and complete [15

15. A. Aydin and A. Hizal, “On the completeness of the spherical vector wave functions,” J. Math. Anal. Appl. 117, 428–440 (1986). [CrossRef]

, 16

16. A. Doicu, T. Wriedt, and Y. Eremin, Light Scattering by Systems of Particles (Springer, 2006). [CrossRef]

] on surfaces without sharp edges [17

17. Complete sets of functions exist on surfaces (Lyapunov surfaces) that are mathematically characterized by three conditions: the normal is well defined at every point; the angle between the normals at any two points on the surface is bounded from above by a function of the distance between these points; all the lines parallel to a normal at an arbitrary point on the surface intercept only once the patches of surface contained in balls centered at the point and smaller than a critical value [18].

]. We choose two sets of electric and magnetic multipoles, {i˜n}n=1 for internal fields and {s˜n}n=1 for scattered fields, centered at different positions within the particle [19

19. A. Doicu and T. Wriedt, “Calculation of the T matrix in the null-field method with discrete sources,” J. Opt. Soc. Am. A 16, 2539–2544 (1999). [CrossRef]

, 20

20. A. Doicu and T. Wriedt, “Extended boundary condition method with multipole sources located in the complex plane,” Opt. Commun. 139, 85–91 (1997). [CrossRef]

]. Any function in ℋ can be approximated to arbitrary precision by a sufficiently large, but finite, number of multipoles [21

21. T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res. 38, 47–95 (2002). [CrossRef]

]; that is, {i˜n}n=1{s˜n}n=1 is complete and no function in this set is the closure of the linear combinations of all the remaining functions. We remark that both the internal and scattered fields exist in real metallic and dielectric particles and fulfill the boundary conditions in Eq. (1) for the electric and magnetic tangential components. For this reason the interaction of light with these particles is determined by surface fields f with four components, and completeness in the space ℋ of the surface fields f is provided by the union of internal and scattered fields and not by either the scattered or the internal field separately. This point is illustrated by the spherical particles considered in Mie theory, where both internal and scattered modes are necessary to form a complete basis.

One can show [21

21. T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res. 38, 47–95 (2002). [CrossRef]

] that the coefficients of the internal and scattered fields, {a˜ni,a˜ns}, that minimize the discrepancy between an incident field and the expansion of internal and scattered fields, |f0+n=1Na˜nss˜na˜nii˜n|, are the solutions of
[ϒ˜ϒ˜ϒ˜Σ˜Σ˜ϒ˜Σ˜Σ˜][a˜ia˜s]=[ϒ˜f0Σ˜f0],
(2)
with a˜i/s=[a˜1i/s,,a˜Ni/s]T and ϒ̃, Σ̃ being matrices whose columns are the functions {i˜n}n=1N and {s˜n}n=1N. The linear independence and completeness of the functions used guarantees that the Gram matrix in Eq. (2), i.e. the matrix of all possible inner products, can be inverted. However, this is numerical challenging [20

20. A. Doicu and T. Wriedt, “Extended boundary condition method with multipole sources located in the complex plane,” Opt. Commun. 139, 85–91 (1997). [CrossRef]

] and for this reason this approach has received little attention. Here we instead take advantage of the block structure of the matrix in Eq. (2) to provide a complete characterization of the geometry of the internal and scattered fields and analytical expression for the coefficient of the expansions of the fields.

The first step is to find orthogonal modes for the scattering and internal fields: for any number of multipoles, N, we achieve this through the matrix decomposition
ϒ˜=UiQi,
(3)
Σ˜=UsQs,
(4)
where Qi,Qs are invertible matrices that can be found through SVD or QR decomposition [22

22. A. Knyazev and M. Argentati, “Principal angles between subspaces in an A-based scalar product: algorithms and perturbation estimates,” SIAM J. Sci. Comput. 23, 2008–2040 (2002). [CrossRef]

] and Ui,Us are unitary matrices whose columns are the orthogonal internal and scattering modes respectively. Scalar products between internal and scattering modes form a matrix with decomposition
UiUs=ViCVs,
(5)
where C a diagonal matrix with positive elements, and Vi,Vs are unitary matrices acting on the internal and scattered fields, respectively. These identities enable us to simplify the Gram matrix through the transformation
[ViQi100VsQs1][ϒ˜ϒ˜ϒ˜Σ˜Σ˜ϒ˜Σ˜Σ˜][Qi1Vi00Qs1Vs]=[1CC1],
(6)
which leads to the matrix equation
[1CC1][aias]=[ϒf0Σf0],
(7)
On the left hand side of Eq. (7), 1 is the identity matrix and ϒ = UiVi and Σ = UsVs are matrices whose columns are formed by the so called principal internal and scattering modes {in} and {sn}, which are one of the main tools in this theory. ai, as are the coefficients of the principal modes in the field’s expansion. The most important part of our theory is that, because matrix C is diagonal, principal modes are coupled pairwise, i.e., each mode is orthogonal to all but at most one function in the other space. This is the essential feature of the multipoles used in Mie’s theory for spheres. The positive diagonal elements of C define the principal angles, ξn, between sn and in as follows
insn=cos(ξn).
(8)
The terms on the right-hand side of Eq. (8) are the principal cosines [12

12. C. Jordan, “Essai sur la géométrie à n dimension,” Bul. Soc. Math. France 3, 103–174 (1875).

]: cos(ξn) and sin (ξn) are the statistical correlation [23

23. E. Hannan, “The general theory of canonical correlation and its relation to functional analysis,” J. Aust. Math. Soc. 2, 229–242 (1961/1962). [CrossRef]

] and the orthogonal distance between sn and in.

From the general theory of angles between subspaces [13

13. A. Knyazev, A. Jujusnashvili, and M. Argentati, “Angles between Infinite Dimensional Subspaces with Applications to the Rayleigh-Ritz and Alternating Projectors Methods,” J. Func. Anal. 259, 1323–1345 (2010). [CrossRef]

] and the definition above, the angles ξ are invariant under unitary transformation of the multipoles and they completely characterize the geometry of the subspaces of the internal and scattered solutions in ℋ. This geometry is induced by the particular scattering particle through the surface integrals of the scalar products; its relevance to scattering and resonances has not been previously realized. The importance of the principal cosines is twofold: Theoretically they provide analytic equations for the coefficients of the internal and scattered principal modes, generalizing the Mie formulae and clarifying the nature of all scattering channels of a particle. Numerically, they allow us to reduce large matrices to their sub-blocks and eliminate the need for numerical inversion to determination of the mode coefficients. For spherical particles, each pair of modes corresponds to a pair of electric or magnetic multipoles of Mie theory. For non-spherical particles, principal modes are instead combinations of different multipoles (although in some cases there can be dominant contributions from a specific multipole).

The interaction of particles with light can now be interpreted in terms of eigenvalues and orthogonal eigenvectors, wn±=(in±sn)/2, of the hermitian operator in Eq. (7); providing interesting analogies between the electromagnetic response of a classical particles with the response of atoms or molecules. However, away from the surface one measures either internal or scattered fields, so we transform the eigenfunctions, { wn±}, to find the coefficients of the principal modes:
ani=incos(ξn)snsin2(ξn)f0=inininf0,
(9)
ans=sncos(ξn)insin2(ξn)f0=snsnsnf0.
(10)
Here in = in – cos(ξn)sn, sn = sn – cos(ξn)in are bi-orthogonal to in, sn (in · sn = sn · in = 0) with in · in = sn · sn = sin2 (ξn). Both the principal or the bi-orthogonal modes fully specify the response of the particle at any point outside and inside the particle. This is shown by recasting the expansions of internal and scattered field as
Fs/i(x)=GS(x,s)f0(s)=(𝒯i(x)In(x)in(s)inin𝒯s(x)Sn(x)sn(s)snsn)f0(s),
(11)
where GS(x, s) is the surface Green’s function [14

14. K. Holms, B. Hourahine, and F. Papoff, “Calculation of internal and scattered fields of axisymmetric nanoparticles at any point in space,” J. Opt. A, Pure Appl. Opt. 11, 054009 (2009). [CrossRef]

,21

21. T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res. 38, 47–95 (2002). [CrossRef]

] of the particle and 𝒯i(x) (𝒯s(x)) is 1 inside (outside) the particle and null elsewhere. In practice, because the principal modes are combinations of known solutions of the Maxwell’s equations, propagation of the fields away from the surface (I(x) and S(x)) is performed for Eq. (11) by evaluation of Bessel or Hankel functions and vector spherical harmonics (all at a very low computational cost). Eq. (11) shows that the convergence of principal modes and principal angles as N → ∞ is a consequence of the convergence of the surface Green’s function [21

21. T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res. 38, 47–95 (2002). [CrossRef]

] for any complete set of solutions of the Maxwell equations. This convergence can be monitored by the surface residual |f0+n=1Nanssnaniin|, which provides an upper bound for the maximum error of scattered and internal fields that decreases with the distance from the surface [16

16. A. Doicu, T. Wriedt, and Y. Eremin, Light Scattering by Systems of Particles (Springer, 2006). [CrossRef]

, 24

24. B. Hourahine, K. Holms, and F. Papoff, “Accurate near and far field determination for non spherical particles from Mie-type theory,” submitted (2011).

]. Furthermore, the form of Eq. (9) remains unchanged as N → ∞, even if θn, in, sn change [25

25. The angles relevant to this work are the point angles 0 < ξ < π/2 of the infinite dimensional theory [13], together with the corresponding subspaces (principal modes) and their orthogonal complements (bi-orthogonal modes).

].

To fully characterize the interaction of particles with electromagnetic fields, we need also to evaluate the ability of principal modes to transport energy. This is determined by the integral of the Poynting vector of each mode over the surface of the particle,
Φns/i=sRe(n^s/iEns/i×Hns/i*)ds,
(12)
where s/i is the outward (inward) pointing normal to the surface for scattered (internal) modes and Ens/i, Hns/i are the electric and magnetic components of the principal mode with |ani/s|=1. Large values of this “intrinsic” mode flux, Φns, corresponds to radiative modes that are very effective at transporting energy from the particle surface to infinity, while small values correspond to modes that are mainly confined near the surface and can only transport energy effectively close to the surface, but not into the far field region. Analogously, large (small) values of Φni correspond to strongly (weakly) absorbing modes.

We now show that this approach gives us an unambiguous way to generalize the definition of Mie resonances to any particle without sharp edges, while also being computationally efficient. Principal modes and scalar products are functions of the frequency dependent permittivity and susceptibility of the particle, and therefore the principal angles ξn change with the frequency of incident light. Internal and scattered coefficients diverge when the denominators of Eqs. (9) and (10) vanish. This happens when a pair of normalized internal and scattering modes are parallel. For a sphere the angular dependence of internal and scattered modes can be factored out and the condition in = sn can be recast in terms of the amplitude of the electric and magnetic components giving the usual Mie resonance condition, which can be interpreted geometrically in terms of alignment between internal and scattered modes. For spheres, the condition in = sn occurs at complex wavelengths; for real wavelengths, resonances correspond to minima of the principal angles. This is also generally true for any smooth particle because the linear independence and completeness of the principal modes makes perfect alignment impossible. So, as with spherical particles [29

29. M. I. Tribelsky and B. S. Lukyanchuk, “Anomalous light scattering by small particles,” Phys Rev. Lett. 97, 263902 (2006). [CrossRef]

], actual resonances correspond to minimum angles (ξn ≠ 0) of pairs in ℋ, which are also minima of the eigenvalues of the hermitian operator in Eq. (7).

We also need to consider the energy transported by resonances. A resonant pair can have very different Φns and Φni, which means that it can be strong or weak depending on whether it is observed in the near or far field, or if scattering or absorption are measured. Furthermore, one important difference between spherical and non-spherical structures is that the total flux of energy scattered or absorbed (integrals of the Pointing vectors over all directions) is given by the sum of principal mode contributions plus interference terms between modes, which are absent for spherical particles. This is due to the orthogonality of the multipoles used in Mie theory being equivalent to orthogonality of the modes’ fluxes. Hence for non-spherical particles, efficiencies can have strong peaks caused by constructive interference within a group of modes, as well as sharp asymmetric features resulting from Fano-like interference between broad and narrow resonances. In the following, we give examples of both of these phenomena. For all particles presented here we use a fitted dielectric function [30

30. P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125, 164705 (2006). [CrossRef] [PubMed]

, 31

31. P. G. Etchegoin, E. C. Le Ru, and M. Meyer, Erratum: “An analytic model for the optical properties of gold”. J. Chem. Phys. 127, 189901 (2007). [CrossRef]

] for gold.

3. Numerical validation

From the point of view of numerical calculations, this approach is a surface method that relies on the the fact that n=1NaniIn, n=1NansSn converge to the exact fields at any point inside and outside the particle [14

14. K. Holms, B. Hourahine, and F. Papoff, “Calculation of internal and scattered fields of axisymmetric nanoparticles at any point in space,” J. Opt. A, Pure Appl. Opt. 11, 054009 (2009). [CrossRef]

, 21

21. T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res. 38, 47–95 (2002). [CrossRef]

] as N→ ∞. We calculate |f0+n=1Nanssnaniin| to determine the numerical error in the evaluation of the surface fields; also the scattered power both at infinity and on the surface of the particle is evaluated to check that the error in the propagation, i.e. in the evaluation of the special functions, is negligible. Comparison with Mie theory [32

32. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1998). [CrossRef]

] for metallic spheres is shown in Fig. 1 and demonstrates that this method is numerically very accurate for particles with radii between [10−2,4] times the wavelength of light. The timings for these calculations are shown in Table 1.

Fig. 1 (Color online) Comparison between this theory (points) and exact Mie results (lines) for gold spheres of radius 10, 25, 50 and 400 nm. Here the refractive index of the external medium is n = 1.3

Table 1. Total Calculation Time Using this Theory for the Spectra Shown in Fig. 1, Evaluated at 81 Wavelengths Using an Intel Core 2 Duo 2.13 GHz System

table-icon
View This Table

4. Resonances in gold nanorods and nanodiscs

We now consider light interaction with experimentally relevant gold particles in vacuum, without considering the possible presence of supporting substrates.

Fig. 2 (Color online) Calculated modes of the gold nanodisc. a) Scattering, absorption and extinction efficiencies, showing the strong resonance at 613 nm. b) Differential scattering cross section (DSCS) as a function of wavelength, illustrating the dipole nature of this particular resonance. c) and d) The induced fields of the bright resonance and the most strongly aligned dark mode, both at 613 nm.

Figure 3(a) shows the landscape of a group of principal mode angles as a function of the order of the mode’s principal cosines and wavelength of incident light. The modes discussed in this section do not change significantly when the number of sources is increased. The “height” of this landscape is sin−1(ξ), i.e. the largest value of Eqs. (9) and (10) for |f0| = 1; while the shading of the traces overlaid on top show, for each wavelength, the values of the intrinsic mode fluxes Φns normalized into the range [0, 1]. Figure 3(d) shows the same information as Fig. 3(a), but instead for Φni. These two figures do not depend on the incident field and provide a visualization of the properties of the surface Green’s function and principal modes of this particle. We can see that there is only one mode pair with comparable intrinsic fluxes, Φi and Φs, which has a resonance at 613 nm. Most of the other mode pairs are strongly absorbing and not able to effectively scatter energy away from the particle. However, one poorly aligned mode pair is capable of strongly radiating if excited, but very weakly absorbs. We also clearly identify a subset of absorbing pairs that become resonant at short wavelengths around 525 nm.

Fig. 3 (Color online) Principal angle landscapes: the order of the mode pairs according to their principal cosines is on the x-axis, the wavelength of the incident field is on the y-axis and sin−1(ξ), the largest possible values of |ani/s| for |f0| = 1, is on the z-axis. Landscapes are overlaid with traces shaded according to: a) the intrinsic mode fluxes Φns, normalized for each wavelength to the range [0,1], b) the amplitudes |ans| (dependent on incident field), c) the mode fluxes |ans|2Φns again normalized to [0, 1]. d)–f) as a)–c) but for internal modes. Discontinuities in the traces are due to crossing between different modes which occur when the values of their principal cosines become the same, this is visible since the ordering of the mode indices then change. a) and d) show that most pairs either do not transport energy or only show strong absorbance, except the one pair that is resonant at 613 nm. c) and f) show that this one resonant pair accounts for most of the energy absorbed and transported into the far region for this incident field. A resonance in the absorbing pairs is visible at short wavelengths (at ∼ 525 nm) in d), but not in e) and f) showing that those modes are not able to coupled to the particular illuminating axial incident field, but would be strongly absorbing if excited by other fields.

Figures 3(b) and 3(e) show the same landscape, but with the amplitudes, |ani/s|, of the internal and scattered principal modes due to the specific axial incident field overlaid on top. This illustrates the excitation paths of the modes as the wavelength changes. We can see that the short-wavelength absorbing resonance around 525 nm is not effectively excited by the specific incident field used; similarly the weakly aligned mode pair which does not pass through a resonance in this range and is capable of strong scattering but weak absorption shows excitation, but only of its internal mode. In contrast the the internal and scattered amplitudes of the pair which reaches resonance at 613 nm is strongly excited by this field.

Figures 3(c) and 3(f) show the effective flux carried by the modes for the incident field used ( Φni/s|ai/s|2) again normalized to [0,1]. By comparing all of these figures, we see that the observed resonance is a Mie-like single principal cosine pair, which reaches maximum alignment (its smallest value of ξ) at the resonance. The internal field also contains a second more weakly aligned excited mode, where its counterpart in the scattering field is not excited. This mode does not absorb much energy on its own, as shown in Fig. 3(f), but affects the surface current, which can be found using these fields and the Ohm equation. The internal and scattered near field of the resonant pair, an electric dipole, is shown in Fig. 2(c). Other strongly aligned mode pairs are not excited, i.e. are dark: this is because they rapidly vary at the surface as in Fig. 2(d), so do not couple to the smoothly varying incident field. The appreciable asymmetry of the absorption and scattering efficiencies are explained by the asymmetry in the principal cosine of the resonant mode as a function of wavelength. For non axial incident light, the main peak in Fig. 2(a) becomes smaller, due to a weaker coupling with this resonant mode (with a corresponding decrease in the amplitudes of Figs. 3(b) and 3(e).

Figure 4(a) and 4(b) shows the calculated optical efficiencies and DSCS of a 480 nm long rod with diameter 40 nm, illuminated axially by plane-wave light. There is a strong resonance at 205 nm and a weaker absorption peak at ∼515 nm. The DSCS demonstrates that this particle strongly scatters this incident light forward, particularly at short wavelengths. Figures 5(a) and 5(d) show that most mode pairs are either strongly radiative or absorbing, except for one weakly aligned pair that is both absorbing and radiating. The peak at 205 nm is not a single mode, but instead the excitation amplitudes of the principal modes, Figs. 5(b) and 5(c), show it to be due to constructive interference between a group of several weakly aligned principal mode pairs.

Fig. 4 Rounded nanorod when illuminated axially. a) (Multimedia online) Scattering, absorption and extinction efficiencies showing the 205 nm mode and weaker absorption peak at 486 nm. Excitation of other resonances for different angles of incidence is shown online. b) (Color online) DSCS showing the strong forward scattering of light by the particle.
Fig. 5 (Color online) Principal angle landscapes as in Fig. 3 for the the rounded nanorod illuminated axially shown in Fig. 4. b) and c) show that the peak at 205 nm is a multimode resonance due to the excitation amplitudes of several principal modes that contribute similarly to the radiative flow, although two are also absorbing modes. The weak absorption peak at 515 nm is also due to a group of internal modes, which become more strongly aligned at around this resonance. In these pairs the coupling of internal and scattering modes is different and the far field is dominated by the radiative mode with the weakest alignment, as shown in c).

Many of these mode pairs are important mainly in the near field, but collectively they very efficiently extract energy from the incident field. Such multimode resonances are potentially very useful for sensing applications because interference between different scattering channels leads to enhanced sensitivity to perturbations near to the scattering surface. The absorption feature at ∼515 nm is due to a small group of modes, two of which are weakly aligned, both absorbing comparable amounts of energy and becoming more strongly aligned at around this resonance. The scattering mode of the pair with the weakest alignment is instead dominant, as shown in Fig. 5(c), where the resonance at ∼515 nm is enhanced, but in absolute terms this resonance is barely observable in the far field scattering efficiency of Fig. 4(a). Figure 6 shows the same particle illuminated equatorially with an incident light polarization of 45° with respect to its long axis. Analyzing other particles with the same diameter but varying length, we find that the broad feature at around 200–450 nm, containing structures similar to the composite modes of Fig. 4, is insensitive to the particle length. It also shows no clear hot or cold spots, see Fig. 6(c). The sharp resonance at 676 nm shifts with rod length, and its surface field as shown in Fig. 6(d), has the strong nodal local structure of a “waveguide” mode on the long axis, remarkably similar to the experimental results of Ref. [35

35. H. Okamoto and K. Imura, “Near field optical imaging of enhanced electric fields and plasmon waves in metal nanostructures,” Prog. Surf. Sci. 84, 199–229 (2009). [CrossRef]

]. Figure 7 shows only mode pairs that cannot be excited by symmetry for axial incidence; the broad features at short wavelength are multimode resonances similar to the one discussed for axial incidence, and most of the corresponding modes are not shown as a result. Similar to the disc, in Fig. 7 most mode pairs are either absorbing or radiating and the three resonances at around 550–700 nm are all absorbing. The excitation paths in Figs. 6(b) and 6(e) show that only the best aligned of the three resonant pairs is excited for this particular field, this resonance is not visible for axial incident light, but becomes excited as the angle of incidence is rotated towards the equator at 90° (Fig. 6, Media 1); the more weakly aligned mode pair at ∼600 nm becomes excited instead for incident angles of ∼ 50°. Figure 6(d) shows that near the sharp resonance at 676 nm, the energy is transported into the far field by the resonant mode and by a weakly aligned, non-resonant, scattering mode. This leads to a Fano-like asymmetric feature in the total scattering cross-section that is sharper and more asymmetric than the single pair resonance, this is due to the interference between these two modes. On the contrary, absorption being due only to the resonant mode, produces a feature in the absorption cross-section which is symmetric with its peak coinciding with the maximum alignment of the mode pair.

Fig. 6 (Color online) The particle from Fig. 4, but illuminated from the side with an incident light polarization of 45° with respect to its long axis. a) Far field scattering efficiencies showing the presence of both a broad feature similar to Fig. 4 at around 200–450 nm and also a sharp resonance at 676 nm, the effect of changing the incident angle on these efficiencies is shown in Media 1. Note that the scattering efficiency shows a strongly asymmetric Fano-like resonance that is sharper than the resonance in the extinction spectrum. b) The DSCS for equatorial illumination. c) and d) Near field for the 480 nm rod, shown at the broad feature (207 nm) and at the “waveguide” Mie-like mode at 676 nm.
Fig. 7 (Color online) Principal angle landscapes as in Fig. 3 for the the rounded nanorod illuminated from the side as in Fig. 6. The sharp mode is again a single well aligned principal pair, the best aligned of three (and the only one excited by this particular field). c) Shows that the energy is transported into the far field by the resonant mode and by a weakly aligned, non-resonant, scattering mode. These modes also interfere to give the the total scattering cross section and produce the characteristic asymmetric and sharpened Fano-like resonance shown in the scattering efficiencies of Fig. 6(a).

5. Conclusion

In summary we provide a theory that generalizes Mie’s formulae to all particles without sharp edges, and reveals the common mechanism behind resonances in nanoparticles and excess noise in macroscopic cavities. This enables a very detailed analysis of the near and far field properties of the electromagnetic response of particles by finding expansions of the internal and scattered fields in terms of modes that are coupled pairwise. A particle can absorb more energy than it scatters or vice-versa, depending on the ability of the scattering and internal modes to transport energy and on their coupling with the incident field. We further show that there are sharp resonances caused by strong alignment between one pair of internal and one scattered mode, and broad resonances due to several pairs of modes having weak, but similar, alignment. We also provide an example of Fano-type resonance in the total scattering cross section of a particle due to interference effects between a resonant and a non-resonant mode, this effect is not possible for spherical particles. We also find that even simple particles that are far smaller than the wavelength of light posses many modes that are dark with respect to incident fields not matching their surface structure. These modes however play an important role for other incident fields, especially when applying near field excitation. This approach can be generalized to assemblies of particles and to complex particles, and is applicable to any process where the interaction between the system and the environment is described by internal and external functions complete at the boundary, including scattering of acoustic or electron waves and coupling to optical cavities.

References and links

1.

D. Graham and R. Goodacre, “Chemical and bioanalytical applications of surface enhanced Raman scattering spectroscopy,” Chem. Soc. Rev. 37, 883–884 (2008). [CrossRef] [PubMed]

2.

T. Aoki, B. Dayan, E. Wilcut, W. P. Bowen, A. S. Parkins, T. J. Kippenberg, K. J. Vahala, and H. J. Kimble, “Observation of strong coupling between one atom and a monolithic microresonator,” Nature 443, 671–674 (2006). [CrossRef] [PubMed]

3.

B. Lukyanchuk, N. Zheludev, S. Maier, N. Halas, P. Nordlander, H. Giessen, and C. Tow Chong, “The Fano resonance in plasmonic nanostructures and metamaterials,” Nat. Mater. 9, 707–715 (2010). [CrossRef]

4.

N. Liu, L. Langguth, T. Weiss, J. Kästel, M. Fleischhauer, T. Pfau, and H. Giessen, “Plasmonic analogue of electromagnetic induced transparency at the Drude damping limit,” Nat. Mater. 8, 758–762 (2009). [CrossRef] [PubMed]

5.

J. Schuller, E. Barnard, W. Cai, Y. C. Jun, J. White, and M. Brongersma, “Plasmonics for extreme light concentration and manipulation,” Nat. Mater. 9, 193–204 (2010). [CrossRef] [PubMed]

6.

Q. Zhao, J. Zhou, F. Zhang, and D. Lippens, “Mie resonance-based dielectric metamaterials,” Mater. Today 12, 60–69 (2009). [CrossRef]

7.

J. Pendry, D. Schuring, and D. Smith, “Controlling electromagnetic fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

8.

G. Roll and G. Schweiger, “Geometrical optics model of Mie resonances,” J. Opt. Soc. Am. A 17, 1301–1311 (2000). [CrossRef]

9.

G. Mie, “Beiträge zur optik trüber medien, speziell kolloidaler metallösungen,” Ann. Phys. 330, 377–445 (1908). [CrossRef]

10.

Y. Han and Z. Wu, “Scattering of a spheroidal particle illuminated by a gaussian beam,” Appl. Opt. 40, 2501–2509 (2001). [CrossRef]

11.

M. I. Mishchenko, J. H. Hovernier, and L. D. Travis, eds., Light scattering by nonspherical particles: Theory, Measurements and Applications (Academic Press, 2000).

12.

C. Jordan, “Essai sur la géométrie à n dimension,” Bul. Soc. Math. France 3, 103–174 (1875).

13.

A. Knyazev, A. Jujusnashvili, and M. Argentati, “Angles between Infinite Dimensional Subspaces with Applications to the Rayleigh-Ritz and Alternating Projectors Methods,” J. Func. Anal. 259, 1323–1345 (2010). [CrossRef]

14.

K. Holms, B. Hourahine, and F. Papoff, “Calculation of internal and scattered fields of axisymmetric nanoparticles at any point in space,” J. Opt. A, Pure Appl. Opt. 11, 054009 (2009). [CrossRef]

15.

A. Aydin and A. Hizal, “On the completeness of the spherical vector wave functions,” J. Math. Anal. Appl. 117, 428–440 (1986). [CrossRef]

16.

A. Doicu, T. Wriedt, and Y. Eremin, Light Scattering by Systems of Particles (Springer, 2006). [CrossRef]

17.

Complete sets of functions exist on surfaces (Lyapunov surfaces) that are mathematically characterized by three conditions: the normal is well defined at every point; the angle between the normals at any two points on the surface is bounded from above by a function of the distance between these points; all the lines parallel to a normal at an arbitrary point on the surface intercept only once the patches of surface contained in balls centered at the point and smaller than a critical value [18].

18.

V. S. Vladimirov, Equations of mathematical physics (MIR, Moscow, 1984).

19.

A. Doicu and T. Wriedt, “Calculation of the T matrix in the null-field method with discrete sources,” J. Opt. Soc. Am. A 16, 2539–2544 (1999). [CrossRef]

20.

A. Doicu and T. Wriedt, “Extended boundary condition method with multipole sources located in the complex plane,” Opt. Commun. 139, 85–91 (1997). [CrossRef]

21.

T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res. 38, 47–95 (2002). [CrossRef]

22.

A. Knyazev and M. Argentati, “Principal angles between subspaces in an A-based scalar product: algorithms and perturbation estimates,” SIAM J. Sci. Comput. 23, 2008–2040 (2002). [CrossRef]

23.

E. Hannan, “The general theory of canonical correlation and its relation to functional analysis,” J. Aust. Math. Soc. 2, 229–242 (1961/1962). [CrossRef]

24.

B. Hourahine, K. Holms, and F. Papoff, “Accurate near and far field determination for non spherical particles from Mie-type theory,” submitted (2011).

25.

The angles relevant to this work are the point angles 0 < ξ < π/2 of the infinite dimensional theory [13], together with the corresponding subspaces (principal modes) and their orthogonal complements (bi-orthogonal modes).

26.

G. New, “The origin of excess noise,” J. Mod. Opt. 42, 799–810 (1995). [CrossRef]

27.

W. J. Firth and A. Yao, “Giant excess noise and transient gain in misaligned laser cavities,” Phys. Rev. Lett. 95, 073903 (2005). [CrossRef] [PubMed]

28.

F. Papoff, G. D’Alessandro, and G.-L. Oppo, “State dependent pseudoresonances and excess noise,” Phys. Rev. Lett. 100, 123905 (2008). [CrossRef] [PubMed]

29.

M. I. Tribelsky and B. S. Lukyanchuk, “Anomalous light scattering by small particles,” Phys Rev. Lett. 97, 263902 (2006). [CrossRef]

30.

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys. 125, 164705 (2006). [CrossRef] [PubMed]

31.

P. G. Etchegoin, E. C. Le Ru, and M. Meyer, Erratum: “An analytic model for the optical properties of gold”. J. Chem. Phys. 127, 189901 (2007). [CrossRef]

32.

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1998). [CrossRef]

33.

Evaluation of the data for the disc at 81 wavelengths required 418 seconds using the same machine as described in the caption of Table 1.

34.

J. Aizpurua, P. Hanarp, D. Sutherland, M. Kall, G. Bryant, and F. J. G. de Abajo, “Optical properties of gold nanorings,” Phys. Rev. Lett. 90, 057401 (2003). [CrossRef] [PubMed]

35.

H. Okamoto and K. Imura, “Near field optical imaging of enhanced electric fields and plasmon waves in metal nanostructures,” Prog. Surf. Sci. 84, 199–229 (2009). [CrossRef]

OCIS Codes
(160.3900) Materials : Metals
(290.0290) Scattering : Scattering
(160.4236) Materials : Nanomaterials
(290.5825) Scattering : Scattering theory

ToC Category:
Scattering

History
Original Manuscript: July 20, 2011
Revised Manuscript: September 23, 2011
Manuscript Accepted: September 28, 2011
Published: October 17, 2011

Citation
F. Papoff and B. Hourahine, "Geometrical Mie theory for resonances in nanoparticles of any shape," Opt. Express 19, 21432-21444 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-22-21432


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. D. Graham and R. Goodacre, “Chemical and bioanalytical applications of surface enhanced Raman scattering spectroscopy,” Chem. Soc. Rev.37, 883–884 (2008). [CrossRef] [PubMed]
  2. T. Aoki, B. Dayan, E. Wilcut, W. P. Bowen, A. S. Parkins, T. J. Kippenberg, K. J. Vahala, and H. J. Kimble, “Observation of strong coupling between one atom and a monolithic microresonator,” Nature443, 671–674 (2006). [CrossRef] [PubMed]
  3. B. Lukyanchuk, N. Zheludev, S. Maier, N. Halas, P. Nordlander, H. Giessen, and C. Tow Chong, “The Fano resonance in plasmonic nanostructures and metamaterials,” Nat. Mater.9, 707–715 (2010). [CrossRef]
  4. N. Liu, L. Langguth, T. Weiss, J. Kästel, M. Fleischhauer, T. Pfau, and H. Giessen, “Plasmonic analogue of electromagnetic induced transparency at the Drude damping limit,” Nat. Mater.8, 758–762 (2009). [CrossRef] [PubMed]
  5. J. Schuller, E. Barnard, W. Cai, Y. C. Jun, J. White, and M. Brongersma, “Plasmonics for extreme light concentration and manipulation,” Nat. Mater.9, 193–204 (2010). [CrossRef] [PubMed]
  6. Q. Zhao, J. Zhou, F. Zhang, and D. Lippens, “Mie resonance-based dielectric metamaterials,” Mater. Today12, 60–69 (2009). [CrossRef]
  7. J. Pendry, D. Schuring, and D. Smith, “Controlling electromagnetic fields,” Science312, 1780–1782 (2006). [CrossRef] [PubMed]
  8. G. Roll and G. Schweiger, “Geometrical optics model of Mie resonances,” J. Opt. Soc. Am. A17, 1301–1311 (2000). [CrossRef]
  9. G. Mie, “Beiträge zur optik trüber medien, speziell kolloidaler metallösungen,” Ann. Phys.330, 377–445 (1908). [CrossRef]
  10. Y. Han and Z. Wu, “Scattering of a spheroidal particle illuminated by a gaussian beam,” Appl. Opt.40, 2501–2509 (2001). [CrossRef]
  11. M. I. Mishchenko, J. H. Hovernier, and L. D. Travis, eds., Light scattering by nonspherical particles: Theory, Measurements and Applications (Academic Press, 2000).
  12. C. Jordan, “Essai sur la géométrie à n dimension,” Bul. Soc. Math. France3, 103–174 (1875).
  13. A. Knyazev, A. Jujusnashvili, and M. Argentati, “Angles between Infinite Dimensional Subspaces with Applications to the Rayleigh-Ritz and Alternating Projectors Methods,” J. Func. Anal.259, 1323–1345 (2010). [CrossRef]
  14. K. Holms, B. Hourahine, and F. Papoff, “Calculation of internal and scattered fields of axisymmetric nanoparticles at any point in space,” J. Opt. A, Pure Appl. Opt.11, 054009 (2009). [CrossRef]
  15. A. Aydin and A. Hizal, “On the completeness of the spherical vector wave functions,” J. Math. Anal. Appl.117, 428–440 (1986). [CrossRef]
  16. A. Doicu, T. Wriedt, and Y. Eremin, Light Scattering by Systems of Particles (Springer, 2006). [CrossRef]
  17. Complete sets of functions exist on surfaces (Lyapunov surfaces) that are mathematically characterized by three conditions: the normal is well defined at every point; the angle between the normals at any two points on the surface is bounded from above by a function of the distance between these points; all the lines parallel to a normal at an arbitrary point on the surface intercept only once the patches of surface contained in balls centered at the point and smaller than a critical value [18].
  18. V. S. Vladimirov, Equations of mathematical physics (MIR, Moscow, 1984).
  19. A. Doicu and T. Wriedt, “Calculation of the T matrix in the null-field method with discrete sources,” J. Opt. Soc. Am. A16, 2539–2544 (1999). [CrossRef]
  20. A. Doicu and T. Wriedt, “Extended boundary condition method with multipole sources located in the complex plane,” Opt. Commun.139, 85–91 (1997). [CrossRef]
  21. T. Rother, M. Kahnert, A. Doicu, and J. Wauer, “Surface Green’s Function of the Helmholtz Equation in Spherical Coordinates,” Prog. Electromagn. Res.38, 47–95 (2002). [CrossRef]
  22. A. Knyazev and M. Argentati, “Principal angles between subspaces in an A-based scalar product: algorithms and perturbation estimates,” SIAM J. Sci. Comput.23, 2008–2040 (2002). [CrossRef]
  23. E. Hannan, “The general theory of canonical correlation and its relation to functional analysis,” J. Aust. Math. Soc.2, 229–242 (1961/1962). [CrossRef]
  24. B. Hourahine, K. Holms, and F. Papoff, “Accurate near and far field determination for non spherical particles from Mie-type theory,” submitted (2011).
  25. The angles relevant to this work are the point angles 0 < ξ < π/2 of the infinite dimensional theory [13], together with the corresponding subspaces (principal modes) and their orthogonal complements (bi-orthogonal modes).
  26. G. New, “The origin of excess noise,” J. Mod. Opt.42, 799–810 (1995). [CrossRef]
  27. W. J. Firth and A. Yao, “Giant excess noise and transient gain in misaligned laser cavities,” Phys. Rev. Lett.95, 073903 (2005). [CrossRef] [PubMed]
  28. F. Papoff, G. D’Alessandro, and G.-L. Oppo, “State dependent pseudoresonances and excess noise,” Phys. Rev. Lett.100, 123905 (2008). [CrossRef] [PubMed]
  29. M. I. Tribelsky and B. S. Lukyanchuk, “Anomalous light scattering by small particles,” Phys Rev. Lett.97, 263902 (2006). [CrossRef]
  30. P. G. Etchegoin, E. C. Le Ru, and M. Meyer, “An analytic model for the optical properties of gold,” J. Chem. Phys.125, 164705 (2006). [CrossRef] [PubMed]
  31. P. G. Etchegoin, E. C. Le Ru, and M. Meyer, Erratum: “An analytic model for the optical properties of gold”. J. Chem. Phys.127, 189901 (2007). [CrossRef]
  32. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1998). [CrossRef]
  33. Evaluation of the data for the disc at 81 wavelengths required 418 seconds using the same machine as described in the caption of Table 1.
  34. J. Aizpurua, P. Hanarp, D. Sutherland, M. Kall, G. Bryant, and F. J. G. de Abajo, “Optical properties of gold nanorings,” Phys. Rev. Lett.90, 057401 (2003). [CrossRef] [PubMed]
  35. H. Okamoto and K. Imura, “Near field optical imaging of enhanced electric fields and plasmon waves in metal nanostructures,” Prog. Surf. Sci.84, 199–229 (2009). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Supplementary Material


» Media 1: AVI (284 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited