Complex bound and leaky modes in chains of plasmonic nanospheres

doi:10.1364/OE.19.018345

« Show journal navigation

Complex bound and leaky modes in chains of plasmonic nanospheres

Salvatore Campione, Sergiy Steshenko, and Filippo Capolino »View Author Affiliations

Optics Express, Vol. 19, Issue 19, pp. 18345-18363 (2011)
http://dx.doi.org/10.1364/OE.19.018345

View Full Text Article

Acrobat PDF (1284 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (55)
Cited By
Figures (9)
Metrics

Abstract

Bound and leaky modes with complex wavenumber in chains (linear arrays) of plasmonic nanospheres are characterized for both longitudinal and transverse polarization states (with respect to the array axis). The proposed method allows for the description of each mode evolution when varying frequency. As a consequence, full characterization of the guided modes with complex wavenumber is provided in terms of propagation direction, guidance or radiance, proper or improper, and physical or nonphysical conditions. Each nanosphere is modeled according to the single dipole approximation, and the metal permittivity is described by the Drude model. Modal wavenumbers are obtained by computing the complex zeroes of the homogeneous equation characterizing the field in the one dimensional periodic array. The required periodic Green’s function is analytically continued into the complex wavenumber space by using the Ewald method. Furthermore, a parametric analysis of the mode wavenumbers is performed with respect to the geometrical parameters of the array.

1. Introduction

Periodic arrays of plasmonic nanospheres have been studied because of their interesting performance in several practical applications. For example, to achieve near-field enhancement and super-resolution at optical wavelengths [1

P. Alitalo, C. Simovski, A. Viitanen, and S. Tretyakov, “Near-field enhancement and subwavelength imaging in the optical region using a pair of two-dimensional arrays of metal nanospheres,” Phys. Rev. B 74(23), 235425 (2006). [CrossRef]

–5

S. Steshenko, F. Capolino, P. Alitalo, and S. A. Tretyakov, “Effective model and investigation of the near-field enhancement and subwavelength imaging properties of multilayer arrays of plasmonic nanospheres,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 84(1), 016607 (2011). [CrossRef] [PubMed]

]; to achieve directive radiation [6

A. F. Koenderink, “Plasmon nanoparticle array waveguides for single photon and single plasmon sources,” Nano Lett. 9(12), 4228–4233 (2009). [CrossRef] [PubMed]

X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B 82(14), 144305 (2010). [CrossRef]

]; to perform microscopy and spectroscopy using SERS [8

J. Beermann, S. M. Novikov, K. Leosson, and S. I. Bozhevolnyi, “Surface enhanced Raman imaging: periodic arrays and individual metal nanoparticles,” Opt. Express 17(15), 12698–12705 (2009). [CrossRef] [PubMed]

F. Liu, Z. Cao, C. Tang, L. Chen, and Z. Wang, “Ultrathin diamond-like carbon film coated silver nanoparticles-based substrates for surface-enhanced Raman spectroscopy,” ACS Nano 4(5), 2643–2648 (2010). [CrossRef] [PubMed]

]; to design biosensors [10

I. Firkowska, S. Giannona, J. A. Rojas-Chapana, K. Luecke, O. Brustle, and M. Giersig, “Biocompatible Nanomaterials and Nanodevices Promising for Biomedical Applications ” in Nanomaterials for Application in Medicine and Biology , M. Giersig, and G. B. Khomutov, eds. (Springer, Berlin, 2008), p. I.1.

,11

A. J. Haes and R. P. Van Duyne, “A nanoscale optical biosensor: sensitivity and selectivity of an approach based on the localized surface plasmon resonance spectroscopy of triangular silver nanoparticles,” J. Am. Chem. Soc. 124(35), 10596–10604 (2002). [CrossRef] [PubMed]

], and DNA detectors [12

I. E. Sendroiu and R. M. Corn, “Nanoparticle diffraction gratings for DNA detection on photopatterned glass substrates,” Biointerphases 3(3), FD23–FD29 (2008). [CrossRef] [PubMed]

]. However, despite these innovative applications, it is still not clear how to associate determinate characteristics with a particular set of excitation wavelengths. Several difficulties are still encountered in the analysis and interpretation of mode propagation in one dimensional (1D-periodic, linear chains) and two dimensional (2D-periodic) arrays of metallic nanospheres when varying frequency, for example. Based on what it has been published so far it is somehow hard to discern if a mode is physical (excitable) or not. The scope of this paper is the characterization of the bound (non radiating) and leaky (radiating) modes with complex wavenumber in linear chains of plasmonic nanospheres (Fig. 1 ), accounting for metal losses, as we recently did in [13

A. L. Fructos, S. Campione, F. Capolino, and F. Mesa, “Characterization of complex plasmonic modes in two-dimensional periodic arrays of metal nanospheres,” J. Opt. Soc. Am. B 28(6), 1446–1458 (2011). [CrossRef]

] for 2D-periodic arrays.

Fig. 1 Linear chain of plasmonic nanospheres embedded in a homogeneous medium with permittivity

ε_{h}

. The radius of each nanosphere is r, and the array spatial period is d.

Modes of linear chains of metallic nanospheres, object of this paper, have been studied in [7

X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B 82(14), 144305 (2010). [CrossRef]

,14

S. A. Maier, P. G. Kik, and H. A. Atwater, “Optical pulse propagation in metal nanoparticle chain waveguides,” Phys. Rev. B 67(20), 205402 (2003). [CrossRef]

–26

M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]

], assuming either lossless or lossy metal and using the single dipole approximation (SDA), described in [27

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).

,28

S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 8.1.

]. In particular, in [14

S. A. Maier, P. G. Kik, and H. A. Atwater, “Optical pulse propagation in metal nanoparticle chain waveguides,” Phys. Rev. B 67(20), 205402 (2003). [CrossRef]

] it has been shown that the dispersion relation k-β (where k is the host medium wavenumber, and β is the wavenumber of the guided wave) for a linear chain of gold nanospheres, obtained from finite-difference time-domain (FDTD) calculations, closely agreed with that predicted by the SDA in case of real wavenumbers. In [15

S. Y. Park and D. Stroud, “Surface-plasmon dispersion relations in chains of metallic nanoparticles: an exact quasistatic calculation,” Phys. Rev. B 69(12), 125418 (2004). [CrossRef]

], the authors have calculated the dispersion relation for a linear chain of gold spherical particles in a uniform host in the case of particle dimensions much smaller than the wavelength, including all the multipolar terms within the quasistatic approximation, finding that the dipole approximation was reasonably accurate for a periodicity larger than three times the nanosphere radius, becoming increasingly inaccurate for smaller periodicity. In [16

R. A. Shore and A. D. Yaghjian, “Travelling electromagnetic waves on linear periodic arrays of lossless spheres,” Electron. Lett. 41(10), 578–580 (2005). [CrossRef]

], the SDA with the expression for the nanosphere polarizability taken from Mie theory has been applied to compute the dispersion diagram (the real part of the wavenumber) of lossless chains of gold nanospheres for longitudinal and transverse polarizations, comparing the results with the ones reported in [14

S. A. Maier, P. G. Kik, and H. A. Atwater, “Optical pulse propagation in metal nanoparticle chain waveguides,” Phys. Rev. B 67(20), 205402 (2003). [CrossRef]

,15

S. Y. Park and D. Stroud, “Surface-plasmon dispersion relations in chains of metallic nanoparticles: an exact quasistatic calculation,” Phys. Rev. B 69(12), 125418 (2004). [CrossRef]

]. In [7

X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B 82(14), 144305 (2010). [CrossRef]

,18

D. S. Citrin, “Plasmon-polariton transport in metal-nanoparticle chains embedded in a gain medium,” Opt. Lett. 31(1), 98–100 (2006). [CrossRef] [PubMed]

–21

R. A. Shore and A. D. Yaghjian, “Complex Waves on 1D, 2D, and 3D Periodic Arrays of Lossy and Lossless Magnetodielectric Spheres,” AFRL-RY-HS-TR-2010–0019, (Air Force Research Laboratory, Hanscom, MA 2010).

,26

M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]

], the authors used polylogarithmic functions to analytically continue the validity of the dispersion relation into the complex wavenumber plane, and they also accounted for metal losses. However, conditions to classify proper and improper modes (proper/improper modes are defined as those that vanish/grow getting away from the 1D array axis, and this concept will be better clarified in Sec. 2.2), as well as to distinguish the subset of physical waves allowed in the linear chain from the complete mathematical set of modes, have not been developed yet when using polylogarithmic functions. In particular in [26

M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]

], the authors derived the 1D dispersion relation by considering all multipole field contributions, showing that in the case of a linear chain of plasmonic nanospheres the dispersion relation is well approximated by the one obtained through the SDA with the expression for the nanosphere polarizability taken from Mie theory, as we do in this paper, assuming small nanospheres with respect to the wavelength. Experimentally, modes have been analyzed in linear chains of gold nanoparticles comparing the dispersion diagram results from measurements (real part of the wavenumber) with corresponding ones obtained from theoretical fully retarded solutions and numerical FDTD simulations [22

K. B. Crozier, E. Togan, E. Simsek, and T. Yang, “Experimental measurement of the dispersion relations of the surface plasmon modes of metal nanoparticle chains,” Opt. Express 15(26), 17482–17493 (2007). [CrossRef] [PubMed]

,23

T. Yang and K. B. Crozier, “Dispersion and extinction of surface plasmons in an array of gold nanoparticle chains: influence of the air/glass interface,” Opt. Express 16(12), 8570–8580 (2008). [CrossRef] [PubMed]

]. In [24

A. Alù and N. Engheta, “Effect of small random disorders and imperfections on the performance of arrays of plasmonic nanoparticles,” N. J. Phys. 12(1), 013015 (2010). [CrossRef]

], the effect of small random disorder (e.g., due to fabrication limitations) in linear chains of metallic nanoparticles has been analyzed theoretically. Specifically, guidance properties of these aperiodic arrays (in fact, the periodicity has been broken by a small disorder) have been analyzed, resulting into additional radiation losses for the guided mode. A review of the state of the art of studies performed for both 1D- and 2D-periodic arrays of plasmonic nanospheres, including mode analysis, fabrication possibilities, and possible applications, has been outlined in [25

S. Campione and F. Capolino, “Linear and Planar Periodic Arrays of Metallic Nanospheres: Fabrication, Optical Properties and Applications,” in Selected Topics in Metamaterials and Photonic Crystals , A. Andreone, A. Cusano, A. Cutolo, and V. Galdi, eds. (World Scientific Publishing, 2011), pp. 141–194.

]. To conclude, we want to emphasize that all the aforementioned excellent referenced work did not provide yet a complete characterization and classification of the modes in a linear chain of plasmonic nanospheres, in terms of leaky/bound, proper/improper, forward/backward, including the conditions for their physical existence, which is the goal of the present paper. We want to stress that for a physical mode we mean a mode that can be excited (launched) into the array by a localized source in proximity of the array, or by a defect or by an array truncation. Other modes that are not classified as physical here may, however, be excited by much more complicated source distributions, though this study is not within the scope of this paper. The analysis in this paper is a necessary step towards a thorough understanding of modal description in periodic arrays. Several results and observations provided in this paper have not been provided elsewhere.

The novelty of the approach presented here is that it allows for the tracking and especially for the characterization of the evolution of all the modes varying frequency. Each nanosphere is modeled as a single dipole, by using the SDA [25

,27

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).

,28

S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 8.1.

] and the metal permittivity is described by the Drude model. The computation of the complex zeroes of the dispersion relation of the 1D-periodic array is performed by using a periodic Green’s function (GF), analytically continued into the complex wavenumber space by using the Ewald representation. The computation of the complex roots of complex valued functions is a difficult problem in its own rights. We have (i) implemented the Muller’s method [29

D. E. Muller, “A Method for Solving Algebraic Equations Using an Automatic Computer,” Math. Tables Other Aids Comput. 10(56), 208–215 (1956). [CrossRef]

,30

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipes: the art of scientific computing (Cambridge University Press, 2007).

] (which can be used to find complex zeroes of analytic functions) and (ii) used the built-in IMSL function ZANLY for Fortran [31

“IMSL Fortran Numerical Library,” (Visual Numerics Corporate Headquarters, 2500 Wilcrest Drive, Suite 200, Houston, TX), www.vni.com.

] (which uses Muller’s method) for verification purposes. The two obtained solutions are in excellent agreement with each other. The use of the periodic GF in the computation of the complex zeroes will be explained in Sec. 2.1. As a consequence, full characterization of the modes with complex wavenumber in the linear chain is provided, in terms of guidance/radiance, proper/improper, physical/nonphysical, and propagation direction conditions. In this paper, a new numerical procedure that uses the Ewald representation for the periodic dyadic GF for linear arrays has been introduced, based on previous scalar developments [32

F. Capolino, D. R. Wilton, and W. A. Johnson, “Efficient computation of the 3D Green's function for the Helmholtz operator for a linear array of point sources using the Ewald method,” J. Comput. Phys. 223(1), 250–261 (2007). [CrossRef]

–34

F. T. Celepcikay, D. R. Wilton, D. R. Jackson, and F. Capolino, “Choosing splitting parameters and summation limits in the numerical evaluation of 1-D and 2-D periodic Green's functions using the Ewald method,” Radio Sci. 43(6), RS6S01 (2008). [CrossRef]

]. Our method is alternative to the one used in [18

D. S. Citrin, “Plasmon-polariton transport in metal-nanoparticle chains embedded in a gain medium,” Opt. Lett. 31(1), 98–100 (2006). [CrossRef] [PubMed]

,20

A. Alù and N. Engheta, “Theory of linear chains of metamaterial/plasmonic particles as subdiffraction optical nanotransmission lines,” Phys. Rev. B 74(20), 205436 (2006). [CrossRef]

,21

,26

M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]

] based on polylogarithms, and future studies will be devoted to relate these two techniques.

The structure of the paper is as follows. The mathematical description for the simulation model adopted to perform the analyses described in this paper is briefly summarized in Sec. 2, including a discussion regarding conditions of guidance, radiance and physical existence of modes and their constituent Floquet waves. Then, the characterization of the complex bound and leaky modes in linear chains of plasmonic nanospheres is reported in Sec. 3, considering nanospheres made of silver, and accounting for metal losses. Though our examples relate to an array with a precise set of representative parameters, the analysis, wave classification and results are general. The analysis of guidance and radiance is reported in Sec. 4. Parameterization of modes with respect to the geometrical parameters of the array is shown in Sec. 5. Criteria for physical existence (excitability) of modes are presented in Appendix A. The dyadic form of the Ewald representation of the periodic GF evaluation used in this paper is provided in Appendix B.

2. Theoretical background

The monochromatic time harmonic convention

\exp (- i ω t)

is assumed. Moreover, in the following equations, bold letters refer to vector quantities, a caret on top of a bold letter refers to unit vector quantities, and a bar under a bold letter refers to dyadic quantities.

2.1 Simulation model

As stated in Sec. 1, according to SDA, each plasmonic nanosphere of the linear chain, at optical frequencies, acts as a single electric dipole and can be described by the induced dipole moment

p_{n} = α_{ee} E_{n}^{loc},

(1)

where

α_{ee}

is the electric polarizability of the nanosphere, and

E_{n}^{loc}

is the local field produced by all the nanospheres of the array except the considered n-th nanosphere plus the external incident field to the array [25

,28

S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 8.1.

]. SDA is a good approximation when the metallic nanospheres are used close to their fundamental plasmonic frequency and particle dimensions are much smaller than the wavelength, when

d \geq 3 r

; however, for small period comparable with the spheres’ radius (i.e.,

2 r < d < 3 r

), more accurate results would involve multipole field contributions [27

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).

,35

J. D. Jackson, Classical Electrodynamics (Wiley, 1998).

–37

D. W. Mackowski, “Calculation of total cross sections of multiple-sphere clusters,” J. Opt. Soc. Am. A 11(11), 2851–2861 (1994). [CrossRef]

]. The polarizability of spherical nanoparticles can be expressed, according to Mie theory [27

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).

], by the following equation [25

,27

C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).

,28

S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 8.1.

,38

R. Ruppin, “Evaluation of extended Maxwell-Garnett theories,” Opt. Commun. 182(4–6), 273–279 (2000). [CrossRef]

,39

W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B Condens. Matter 39(14), 9852–9858 (1989). [CrossRef] [PubMed]

]

α_{ee} = \frac{6 π ε_{0} ε_{h} i}{k^{3}} \frac{m ψ_{1} (m k r) ψ_{1}^{'} (k r) - ψ_{1} (k r) ψ_{1}^{'} (m k r)}{m ψ_{1} (m k r) ξ_{1}^{'} (k r) - ξ_{1} (k r) ψ_{1}^{'} (m k r)},

(2)

where r is the nanosphere radius,

ψ_{1} (ρ) = ρ j_{1} (ρ) = \sin ρ / ρ - \cos ρ

ξ_{1} (ρ) = ρ h_{1}^{(1)} (ρ) = (- i / ρ - 1) e^{i ρ}

are the Riccati-Bessel functions [40

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).

ε_{h}

is the relative permittivity of the host medium,

k = ω \sqrt{ε_{h}} / c = k_{0} \sqrt{ε_{h}}

is the host wavenumber, where

k_{0}

denotes the free space wavenumber,

ε_{0}

is the free space absolute permittivity,

m = \sqrt{ε_{m} / ε_{h}}

is the relative refractive index, and

ε_{m}

is the relative permittivity of the metal, which is described by the Drude model

ε_{m} = ε_{\infty} - ω_{p}^{2} / [ω (ω + i γ)]

, where

ε_{\infty}

is a high-frequency fitting parameter,

ω_{p}

is the plasma frequency of the metal (expressed in rad/s) andγis the damping factor (expressed in 1/s).

Consider now a linear chain of nanospheres as in Fig. 1, for which each nanosphere is placed at positions

z_{n} = n d

, with

n = 0, \pm 1, \pm 2, \dots

, and d is the spatial period of the chain. Suppose that the array is immersed in a homogeneous background, with relative permittivity

ε_{h}

, and that the array is then excited by a plane wave or by a quasi-periodic excitation with wavevector

k_{B} = k_{z} \hat{z}

[28

S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 8.1.

]. Consequently, each spherical nanoparticle will have a dipole moment equal to

p_{n} = p_{0} exp (i n k_{z} d)

, with

p_{0}

defined as in Eq. (1) by using n = 0. Then, the local field acting on a nanosphere at position

z_{0} = 0

is given by

E_{0}^{loc} (k_{z}) = E^{inc} (k_{z}) + {\underline{\overset{⌣}{G}}}^{\infty} (k_{z}) \cdot p_{0},

(3)

where

{\underline{\overset{⌣}{G}}}^{\infty} (k_{z})

represents the regularized electric field dyadic GF for the phased periodic array of nanospheres, defined as

{\underline{\overset{⌣}{G}}}^{\infty} (k_{z}) = \sum_{\begin{array}{l} n = - \infty \\ (n \neq 0) \end{array}}^{\infty} {\underline{G}}_{n} e^{i n k_{z} d},

(4)

and

{\underline{G}}_{n}

is the dyadic GF of the background medium. Since the array is placed in a homogeneous medium, the GF (field produced at

z = 0

by a unit dipole at

z = n d

) simply is

{\underline{G}}_{n} = \frac{e^{i k | n | d}}{4 π ε_{0} ε_{h}} [(\frac{k^{2}}{| n | d} + \frac{i k}{{(| n | d)}^{2}} - \frac{1}{{(| n | d)}^{3}}) \underline{I} - (\frac{k^{2}}{| n | d} + \frac{3 i k}{{(| n | d)}^{2}} - \frac{3}{{(| n | d)}^{3}}) \hat{z} \hat{z}],

(5)

where

\underline{I}

is the identity dyad. Substituting the expression for the local field given in Eq. (3) into Eq. (1) for the nanosphere at position

z_{0} = 0

, it then follows that

p_{0} = α_{ee} E^{inc} (k_{z}) + α_{ee} {\underline{\overset{⌣}{G}}}^{\infty} (k_{z}) \cdot p_{0},

which leads to the linear system

\underline{A} (k_{z}) \cdot p_{0} = α_{ee} E^{inc} (k_{z}), \underline{A} (k_{z}) = \underline{I} - α_{ee} {\underline{\overset{⌣}{G}}}^{\infty} (k_{z}) .

(6)

Mode analysis in the linear chain is then performed by computing the complex zeroes of the homogeneous linear system in Eq. (6), after having imposed no excitation source (i.e.,

E^{inc} (k_{z}) = 0

). In other words, the complex

k_{z}

-zeroes of the determinant of

\underline{A} (k_{z})

have to be determined (i.e.,

k_{z}

such that

\det \underline{A} (k_{z}) = 0

) using the methods described in Sec. 1. The series inside the GF definition in Eq. (4) is slowly convergent, and it is even divergent for complex values of the wavenumber

k_{z}

, and some technique must be adopted to allow analytic continuation of the GF into the complex

k_{z}

plane. The GF representation adopted in this work makes use of Ewald’s method for linear arrays, which, besides providing analytic continuation, it also provides series with Gaussian convergence and only a handful of terms are needed [32

–34

]. The dyadic form of the Ewald representation of the periodic GF evaluation used in this paper is provided in Appendix B.

2.2 Floquet waves representation

The electric field in the linear chain can be described as

E (r + n d \hat{z}, k_{z}) = E (r, k_{z}) exp (i n k_{z} d)

, where d is the spatial period of the periodic structure,

n = 0, \pm 1, \pm 2, \dots

k_{z} = β_{z} + i α_{z}

is the wavenumber along the array axis, and

r = x \hat{x} + y \hat{y} + z \hat{z}

. The field relative to a mode in the linear chain is expressed in term of a Fourier series expansion, and thus represented as the superposition of Floquet wave (FW) spatial harmonics as

E^{mode} (r, k_{z}) = \sum_{p = - \infty}^{\infty} e^{i k_{z, p} z} e_{p}^{mode} (x, y, k_{z}),

(7)

where

k_{z, p} = k_{z} + 2 π p / d

, with p the order of the Floquet spatial harmonic, and

e_{p}^{mode}

is the transverse eigenvector of the p-th harmonic. Each longitudinal Floquet wavenumber

k_{z, p} = β_{z, p} + i α_{z}

has the same attenuation constant

α_{z}

, whereas the real part

β_{z, p}

is periodic in the complex wavenumber domain, with period

2 π p / d

. Due to this periodicity, a mode will be described by the wavenumber of its FW in the fundamental Brillouin zone (BZ), defined as

- π / d < β_{z} < π / d

. Therefore here, by definition, the wavenumber

k_{z}

lies in the fundamental BZ. Furthermore, the radial wavenumber is

k_{ρ, p} = \sqrt{k^{2} - k_{z, p}^{2}} = β_{ρ, p} + i α_{ρ, p}

. Proper FWs are defined as those with

α_{ρ, p} > 0

, i.e., those that vanish for

ρ = \sqrt{x^{2} + y^{2}} \to \infty

. Improper FWs are those with

α_{ρ, p} < 0

, i.e., those that grow for

ρ \to \infty

. A proper mode has all its spatial harmonics proper, whereas an improper mode has at least one spatial harmonic improper. Note that certain leaky waves are improper, and physical, but they never grow indefinitely because they exist only in certain regions of space [41

F. Capolino, D. R. Jackson, and D. R. Wilton, “Field Representations in Periodic Artificial Materials Excited by a Source,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 12.1.

,42

L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).

] (in other words, their domain of existence is limited). This means that the mode itself is not present in certain regions, though the radiation caused by a leaky mode is present everywhere as explained in [42

L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).

]).

2.3 Physical excitation conditions by a source, a defect or a truncation

Among all the mathematical wavenumber solutions of

\det \underline{A} (k_{z}) = 0

only a subset represents physical waves, i.e., waves that can be excited by a localized source, or defect, or array truncation, as described in Appendix A. The physical modes that can be excited in linear chains of plasmonic nanospheres, according to Appendix A and what follows, are summarized in Table 1 , similar to what was previously reported for the 2D-periodic case in [41

,43

P. Baccarelli, S. Paulotto, and C. Di Nallo, “Full-wave analysis of bound and leaky modes propagating along 2D periodic printed structures with arbitrary metallisation in the unit cell,” IET Proc. Microwaves, Antennas Propag. 1(1), 217–225 (2007). [CrossRef]

], based on the complex wavenumber of the fundamental harmonic (with p = 0). Table 1 applies to the case of chain period shorter than half wavelength, since it implies that only wavenumbers in the fundamental BZ may represent fast waves.

Table 1 Classification of Physical Modes with Complex Wavenumbers, Based on the Wavenumber of the Mode-Harmonic in the Fundamental BZ

	Forward Wave $β_{z, 0} α_{z} > 0$	Backward Wave $β_{z, 0} α_{z} < 0$
Slow Wave	(A) $\begin{array}{l} \| β_{z, 0} \| > k \\ α_{ρ, 0} > 0 \end{array}$ (proper, bound)	(B) $\begin{array}{l} \| β_{z, 0} \| > k \\ α_{ρ, 0} > 0 \end{array}$ (proper, bound)
Fast Wave	(C) $\begin{array}{l} \| β_{z, 0} \| < k \\ α_{ρ, 0} < 0 \end{array}$ (improper, leaky)	(D) $\begin{array}{l} \| β_{z, 0} \| < k \\ α_{ρ, 0} > 0 \end{array}$ (proper, leaky)

Assume that a point source is near the origin (Fig. 2 ) and an observer is located along the positive z axis of the array. The field at the observer is represented by a spectral path integral on the real axis as shown in the figure in Appendix A. Deformation onto the steepest descent paths encounters one or more periodic sets of poles that represent the wavenumbers of the physical FWs [41

,44

F. Capolino, D. R. Jackson, D. R. Wilton, and L. B. Felsen, “Comparison of methods for calculating the field excited by a dipole near a 2-D periodic material,” IEEE Trans. Antenn. Propag. 55(6), 1644–1655 (2007). [CrossRef]

,45

F. Capolino, D. R. Jackson, and D. R. Wilton, “Fundamental properties of the field at the interface between air and a periodic artificial material excited by a line source,” IEEE Trans. Antenn. Propag. 53(1), 91–99 (2005). [CrossRef]

]. A mode is represented by an infinite set of wavenumbers, according to Eq. (7). Periodic sets of poles that are not captured in the path deformation do not contribute to the final field value and therefore represent modes that are not excited. In order to be captured in the path deformation in Appendix A, a pole representing a wave with complex wavenumber

k_{z} = β_{z} + i α_{z}

should have

α_{z} > 0

. Therefore, considering the modal harmonic in the fundamental BZ (often the one with p = 0), a physical forward wave belongs to the I quadrant of the complex

k_{z}

plane (i.e.,

β_{z} α_{z} > 0

), whereas a physical backward one to the II quadrant (i.e.,

β_{z} α_{z} < 0

). Moreover, to be physical, not only the wave needs to have a wavenumber with

α_{z} > 0

, but also its wavenumber has to respect the conditions stated in Table 1, which are also valid in the case of an observer along the negative z axis for which physical waves have

α_{z} < 0

. Further discussion about integration path and field representation is provided in Appendix A and references therein.

Fig. 2 Illustration of a physical complex wave for an observer along the positive z axis, according to its classification provided in Table 1. (A) and (B) represent forward and backward bound waves. (C) and (D) represent forward and backward leaky waves.

Figure 2 shows how each physical wave outlined in Table 1 evolves in the analyzed linear chain, for the case in which the observer is along the positive z-direction. If the wave is proper, as defined in the previous section, it vanishes for

ρ \to \infty

(as can be observed in Figs. 2(A), 2(B) and 2(D), where the wave decays along the radial direction). If the wave is improper, instead, it grows for

ρ \to \infty

, as shown in Fig. 2(C). Note that in this case, the leaky wave existence domain is limited in space as described at the end of Sec. 2.2, and therefore it does not violate any physical principle [41

,42

L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).

]. If the wave is forward, the phase travels towards the positive z-direction, as highlighted by the red arrows in Figs. 2(A) and 2(C), and attenuates while it travels. Whereas, if the wave is backward, the phase travels towards negative z-direction, as highlighted by the red arrows in Figs. 2(B) and 2(D), but still attenuates along positive z (for conservation of energy, a wave cannot grow along positive z). If the wave is bound, then no energy is radiated away (as in Figs. 2(A) and 2(B)); radiation phenomena are present when the wave is leaky, as shown in Figs. 2(C) and 2(D) by the black arrows.

2.4 Conditions for guidance and radiance

An important part of this work lies in the fact that we show which modes are physical (i.e., that are excitable by a localized source), and consequently employable in practical applications. For this reason, the conditions to discriminate physical bound and radiating (leaky) modes are here stated based on Table 1 and the discussion in Appendix A.

A bound mode is a mode that has all its Floquet wavenumbers outside the “radiating spectral region”, i.e., the spectral region for which

| β_{z, p} | < k

. Therefore, a bound mode is a mode with

| β_{z, p} | > k

for any p. Taking into account the physical condition, considering the wave harmonic with

β_{z}

in the fundamental BZ, an observer along the positive z-direction would be reached by a bound wave (excited by a localized source near the origin) either when (i) it is proper and

k < β_{z}

and

α_{z}

is positive (forward wave), or when (ii) it is proper and

β_{z} < - k

and

α_{z}

is positive (backward wave), i.e., the wavenumber is either in the first or in the second quadrant of the complex

k_{z}

plane. Improper bound modes cannot be launched by a source in Fig. 2.

A leaky mode is a mode that has at least one Floquet wavenumber within the radiating spectral region, or

| β_{z, p} | < k

for some p. Taking into account the physical condition, and considering the harmonic

β_{z}

in the fundamental BZ, an observer along the positive z-direction would be reached by a radiating wave (excited by a localized source near the origin) either when (i) it is improper and

0 < β_{z} < k

and

α_{z}

is positive (forward wave), or when (ii) it is proper and

- k < β_{z} < 0

and

α_{z}

is positive (backward wave).

A guided mode (leaky or bound) travels a long distance when the imaginary part of its wavenumber is small with respect to the host wavenumber, or

α_{z} < < k

(i.e., low decay).

For completeness, we want to emphasize that the higher order Floquet harmonics outside the radiating spectral region (with

| β_{z, p} | > k

p \neq 0

) can be either proper or improper non radiating waves, and only the ones represented by poles captured from the path deformation described in Appendix A can be excited by a source in Fig. 2, and thus be physical. Also, as mentioned in Sec. 2.3, in the following sections we restrict our analysis to the case of chain period shorter than half wavelength, where only wavenumbers in the fundamental BZ may represent fast (leaky) waves.

3. Mode analysis

In this section, complex modes in a linear chain of silver nanospheres embedded in free space (i.e.,

k = k_{0}

) are shown for both transverse and longitudinal polarization states (T- and L-pol, respectively), accounting also for metal losses. The radius of the nanospheres is r = 25 nm and the spatial period is d = 75 nm. The adopted Drude model parameters for silver are

ε_{\infty} = 5

, plasma frequency

ω_{p} = 1.37 \times 10^{16} rad/s

and damping factor

γ = 27.3 \times 10^{12} s^{- 1}

[46

A. Alù, A. Salandrino, and N. Engheta, “Negative effective permeability and left-handed materials at optical frequencies,” Opt. Express 14(4), 1557–1567 (2006). [CrossRef] [PubMed]

,47

I. El-Kady, M. M. Sigalas, R. Biswas, K. M. Ho, and C. M. Soukoulis, “Metallic photonic crystals at optical wavelengths,” Phys. Rev. B 62(23), 15299–15302 (2000). [CrossRef]

]. In the L-pol case, each nanosphere has an induced dipole moment along the z-direction (i.e.,

p = p_{z} \hat{z}

). In the T-pol case, each nanosphere has an induced dipole moment orthogonal to z.

For both polarizations, if

k_{z}

is a complex wavenumber solution of

\det \underline{A} (k_{z}) = 0

, then also

- k_{z}

is a solution. In absence of losses (case not reported here for the sake of brevity) a four-quadrant symmetry holds, thus the correspondent complex conjugate solutions,

\pm k_{z}^{*}

, would be solutions as well. When losses are present, the complex conjugate solutions are modified from

\pm k_{z}^{*}

correspondingly to the amount of losses introduced in the model. Because of space limitation, and without deterring information, we show here only the fundamental p = 0 FW harmonic of complex modes. The periodic condition

k_{z, p} = k_{z} + 2 π p / d

would determine the behavior of FW harmonics with wavenumbers in other BZs. Therefore, the focus is on the wavenumber solutions belonging to the I and II quadrants of the complex

k_{z}

plane in the dispersion diagrams in the fundamental BZ, i.e., those solutions with

α_{z} \geq 0

, which are the only ones that can be captured in the path deformation described in Appendix A when observing along the positive z axis. Note that this subset has the same information as the complete set of modes (i.e., a complete basis for the field) in the linear chain because of the symmetry property stating that

k_{z}

and

- k_{z}

are both solutions.

3.1 Transverse polarization (T-pol)

Figure 3 shows the dispersion diagrams for the real and imaginary parts of the wavenumber of the FW harmonic in the fundamental BZ of the modes with T-pol and relatively small imaginary part. Figure 4 shows the wavenumber evolutions in the complex

k_{z}

plane varying frequency, using two different normalizations. We want to clarify that the reported curves in Figs. 3 and 4 are relative to both physical and nonphysical solutions (because together they form a complete basis for the field); physical ones are then distinguished by using Table 1.

Fig. 3 Dispersion diagram for T-pol. (a) Real part and (b) imaginary part of the wavenumber of the fundamental FW

k_{z} = β_{z} + i α_{z}

Fig. 4 Modal wavenumbers in the complex relative

k_{z}

plane for T-pol, normalized with respect to (a) the spatial period d, and (b) the host wavenumber

k = k_{0} \sqrt{ε_{h}}

(here

ε_{h} = 1

). Note that in (b), crossing “1” in the abscissa means crossing the light line. The shown curves are parametric curves of frequency, and the arrows indicate direction of increasing frequency.

Modes ‘Proper 1’, ‘Proper 2’, ‘Improper 1’ and ‘Improper 2’ in Figs. 3 and 4 are computed using the method described in [33

V. R. Komanduri, F. Capolino, D. R. Jackson, and D. R. Wilton, “Computation of the one-dimensional free-space periodic green's function for leaky waves using the ewald method,” in Proc. URSI Gen. Ass.(Chicago, IL, 2008).

], where only m = 0 and −1 are considered, and adapted to the present case as in Appendix B. When the wavenumber dispersion curves change complex quadrant, they can be continued with dispersion curves computed with different m-values (see appendix B). We also show the solution ‘Improper 3′, computed with m = +1, as a continuation of the m = 0 red curve into the II and IV quadrants (Fig. 4), which also represents the wavenumber close to the complex conjugate of the ‘Improper 1’, according to the symmetry properties of the modal solutions described in Sec. 3. Though these waves computed with m values different from 0 and −1 are correct mathematical solutions, their physical validity requires further investigation, as explained in Appendix B. For an observer along the positive z-direction, according to Table 1 and Sec. 2.4, and looking at the complex

k_{z}

plane I and II quadrants (the upper half) in Fig. 4 (i.e., the only part in which physical solutions could be found), the mode ‘Proper 1’ is backward and physical in the entire II quadrant (specifically, it is a physical leaky mode when

- k < β_{z} < 0

, and physical bound mode when

β_{z} < - k

); similarly, the mode ‘Improper 1’ is a physical forward leaky mode in the I quadrant when

β_{z} < k

, whereas it is a nonphysical bound mode when

β_{z} > k

. The mode ‘Proper 2’ is a physical forward bound mode only in the I quadrant when

β_{z} > k

, whereas the mode ‘Improper 2’ is always nonphysical, being in the II quadrant (remember that improper modes in the II quadrant cannot be captured, according to Appendix A).

From the point of view of modal wavenumber evolution when increasing frequency, the wavenumber of the blue mode ‘Proper 1’ goes from the II to the I quadrant of the complex

k_{z}

plane (or from the IV to the III due to the symmetry of the modal wavenumber solutions with respect to the origin stated in Sec. 3, i.e.,

k_{z}

and

- k_{z}

are both solutions, even in lossy cases), thus it crosses the imaginary axis (exhibiting its lowest

α_{z}

before this transition) and so a branch cut positioned there (as shown in Appendix A). Therefore, this modal wavenumber passes from the top to the bottom Riemann sheet, i.e., it becomes the one of the blue mode ‘Improper 1’. This transition is well observed in Fig. 4. It is noteworthy that due to this change in the Riemann sheet of the modal wavenumber, the physical backward leaky mode ‘Proper 1’ (II quadrant), when frequency increases, transitions into the physical forward leaky mode ‘Improper 1’ (I quadrant). The wavenumber of the red mode ‘Proper 2’ is almost over the light line

β_{z} = k

at low frequencies. Increasing frequency, it departs from the light line in the bound region until its trajectory is inverted (Fig. 4), and then it goes from the I to the II quadrant (or from the III to the IV due to symmetry), and thus, similarly to what described for the mode ‘Proper 1’, crosses a branch cut, transitioning into the wavenumber of the mode ‘Improper 3′. With smaller losses, the wavenumber of the mode ‘Proper 2’ would be able to depart farther away from the light line in the I quadrant, reaching large values of

β_{z}

, and thus causing short plasmonic wavelengths.

3.2 Longitudinal polarization (L-pol)

Figure 5 shows the dispersion diagrams for the real and imaginary parts of the wavenumber of the FW harmonic in the fundamental BZ of the modes with L-pol and relatively small imaginary part. Figure 6 shows the wavenumber evolutions in the complex

k_{z}

plane varying frequency, using two different normalizations. We want to clarify again that the reported curves in Figs. 5 and 6 are relative to both physical and nonphysical solutions; physical ones are then distinguished by using Table 1.

Fig. 5 As in Fig. 3, for L-pol.

Fig. 6 As in Fig. 4, for L-pol. Inset in (b) shows the modal wavenumbers zero crossing.

Modes ‘Proper 1’ and ‘Improper 1’ are computed using the method described in [33

], and adapted to the present case as in Appendix B, where only m = 0 and −1 are considered. When the wavenumber dispersion curves change complex quadrant, they can be continued with different m-values (see appendix B). We also show the curve related to the ‘Improper 2’ mode obtained with m = +1, which represents the continuation of the ‘Proper 1’ solution into the complex

k_{z}

plane quadrants II and IV (Fig. 6). This wavenumber also represents the wavenumber close to the complex conjugate solution symmetric to the ‘Improper 1’ solution, according to the symmetry properties described in Sec. 3. The physical validity of the waves found with m-values different than 0 and −1, requires further investigation, as explained in Appendix B. For an observer along the positive z-direction, according to Table 1 and Sec. 2.4, and looking at quadrants I and II in the upper half part of the complex

k_{z}

plane in Fig. 6, the mode ‘Proper 1’ is a physical forward bound mode in the I quadrant when

β_{z} > k

; similarly, the mode ‘Improper 1’ is a physical forward leaky mode in the I quadrant when

0 < β_{z} < k

. The mode ‘Improper 1’ in the II quadrant is always nonphysical.

From the point of view of modal wavenumber evolution, at high frequencies, the wavenumber of the blue mode ‘Proper 1’ is in the bound region in the I quadrant, as in Fig. 6(b). Decreasing frequency, the real part of its wavenumber increases while its imaginary part decreases up to about

k d / π \approx 0.45

, then both the real and imaginary parts of this modal wavenumber decrease. Around

k d / π \approx 0.39

, this wavenumber enters in the radiating spectral region (i.e.,

| β_{z} | < k

). There, its imaginary part transitions from positive to negative values (see the inset in Fig. 6). This transition translates into the wavenumber passing from the I to the IV quadrant of the complex

k_{z}

plane (or from the III to the II due to symmetry), thus crossing the

Re (k_{z})

axis and so a branch cut located parallel to the

Re (k_{z})

axis (in Appendix A, this branch cut is slightly above the axis when small ambient losses are present; it would be exactly on the axis in absence of ambient losses), and it changes Riemann sheet, i.e., it becomes the wavenumber of the blue mode ‘Improper 2’. This transition is well observed in Fig. 6. The wavenumber of the forward green mode ‘Improper 1’ starts in the I quadrant of the complex plane in Fig. 6. Increasing frequency, its imaginary part goes to zero while its real part remains under the light line (i.e.,

| β_{z} | > k

) as in the inset in Fig. 6, thus it does not encounter any branch cut, as can be seen in Appendix A, and remains improper with negative growing imaginary part (this modal wavenumber transitions from the I to the IV quadrant, or from the III to the II due to symmetry).

Notice for example that in [21

] modes ‘Improper 1’ and ‘Proper 1’ have been described as the same mode, discontinuous for a small frequency range; whereas in [26

M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]

] they have been described as one continuous longitudinal mode.

4. Analysis of modal guidance and radiance

4.1 Bound modes for transverse and longitudinal polarizations

Accounting for the presence of losses in the silver nanospheres, the only bound mode for T-pol that travels without large decay is the mode ‘Proper 2’ in Figs. 3(a) and 3(b). This mode is bound, forward and physical (as described in Sec. 3.1); it presents growing imaginary part for increasing frequency, with very small positive attenuation constant near the light line for

k d / π \leq 0.4

(this mode lies on top of the light line for

k d / π < 0.37

). In the entire domain of its physical existence, this mode can be guided very well, in the sense that it can travel long distances from the excitation source. The other physical bound mode for T-pol is the ‘Proper 1’ that has the imaginary part of its wavenumber which is large in comparison to the host wavenumber k. This mode is bound, backward and physical; it presents first decreasing imaginary part while frequency increases, a minimum at

k d / π \approx 0.4

at which

α_{z} d / π \approx 0.17

, and then increasing imaginary part.

The only bound mode for L-pol that travels without large decay is the mode ‘Proper 1’ in Figs. 5(a) and 5(b). This mode is bound, forward and physical (as described in Sec. 3.2); it presents growing imaginary part for increasing frequency, and can be excited in all its domain of physical existence, although it would be slowly decaying while traveling only when it has very small attenuation constant (e.g., near the light line at about

k d / π \approx 0.39

4.2 Leaky modes for transverse and longitudinal polarizations

The radiating leaky modes for T-pol are described next. The physical modes found in the array of silver nanospheres here analyzed radiate, but they cannot provide very directive radiation since the imaginary part of their wavenumber is large with respect to the host wavenumber k [48

A. A. Oliner and D. R. Jackson, “Leaky-wave antennas,” in Antenna Engineering Handbook , J. Volakis, ed. (McGraw Hill, 2007), p. 11.1.

]. In Figs. 3 and 4, the mode ‘Proper 1’ is a leaky backward mode, whereas the mode ‘Improper 1’ is a leaky forward mode. We will show in Sec. 5 the possibility of modifying the imaginary part of the wavenumber of these two modes by varying geometrical parameters of the array.

The only radiating leaky mode for L-pol is the mode ‘Improper 1’ in Figs. 5 and 6. This mode is leaky and forward. However, notice that this mode lies very close to the light line; since the dipoles approximating the nanospheres are longitudinally polarized (i.e., along the z-direction), their far-field radiation pattern would have the classical

\sin θ

shape, with a far-field null corresponding to the direction along which the dipoles are aligned. Indeed, in the case of the analyzed linear chain, when

β_{z} = k

the radiation should take place along the z-direction and this is not possible due to the aforementioned null, whereas when

β_{z} < k

for

β_{z} \approx k

, a very small radiation (still due to the dipole

\sin θ

radiation pattern) could take place.

A detailed analysis of how to adopt these leaky modes to obtain directive radiation at optical frequencies has been recently shown in [7

X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B 82(14), 144305 (2010). [CrossRef]

5. Parametric analysis of bound and leaky modes

In this section, we analyze the evolution of bound and leaky modes by varying the radius r of each nanosphere and the spatial period d of the array.

For the T-pol case, the dependence with respect to the array geometrical parameters of the backward ‘Proper 1’ (always physical, bound for

β_{z} < - k

, leaky for

- k < β_{z} < 0

) and the forward ‘Improper 1’ (physical leaky for

0 < β_{z} < k

, nonphysical bound for

β_{z} > k

) is illustrated in Fig. 7 . The linear chain under analysis here is the same as the one in Sec. 3 where we vary one parameter at the time. In Fig. 7(a), d = 75 nm and r varies from 15 nm to 35 nm; whereas in Fig. 7(b), r = 25 nm and d varies from 55 nm to 100 nm. It is observed in Fig. 7(a) that increasing the radius of each nanosphere to 35 nm (blue lines), both the ‘Proper 1’ and ‘Improper 1’ modes present lower losses for a wider frequency range than the array in Sec. 3 (green lines). Also, notice that by increasing the nanospheres’ radius, the mode ‘Proper 1’ is “shifted” towards lower frequencies, whereas the mode ‘Improper 1’ is “shifted” towards higher frequencies in Fig. 7(a). It is noteworthy that for r = 35 nm the backward proper mode has a very small imaginary part in correspondence of

k d / π \approx 0.38

(where this mode is bound); if excited at this particular frequency, this mode would travel long distances from the excitation source. Moreover,

β_{z}

is quite large, thus causing short plasmonic wavelengths. In Fig. 7(b), we show the mode dependence on the array period: it can be noticed that the imaginary part

α_{z}

increases with increasing periodicity. Notice that this is in agreement with what it has been previously shown in [20

A. Alù and N. Engheta, “Theory of linear chains of metamaterial/plasmonic particles as subdiffraction optical nanotransmission lines,” Phys. Rev. B 74(20), 205436 (2006). [CrossRef]

Fig. 7 Dependence of the ‘Proper 1’ bound mode (solid line) and ‘Improper 1’ leaky mode (dashed line) for T-pol in Fig. 3 with respect to the geometrical parameters of the linear array. The analyzed array is the same as in Sec. 3. In (a) d = 75 nm and r varies from 15 nm to 35 nm. In (b) r = 25 nm and d varies from 55 nm to 100 nm.

In the L-pol case, the dependence with respect to the array geometrical parameters of the forward ‘Improper 1’ (physical leaky only for

0 < β_{z} < k

, nonphysical bound for

β_{z} > k

) is illustrated in Fig. 8 . The performed parameterization is the same as in Fig. 7. As for the T-pol case previously analyzed, it can be observed in Fig. 8(a) that increasing the radius of each nanosphere to 35 nm (blue lines), the ‘Improper 1’ mode presents lower losses for a wider frequency range than the array in Sec. 3 (green lines). Furthermore, notice that for large nanospheres’ radius the frequency range for which the mode ’Improper 1’ is physical is wider (this range is about 0.01 for the red curve and about 0.21 for the blue curve, using the scale

k d / π

in Fig. 8(a)) and “shifted” towards lower frequencies. In Fig. 8(b), we modify the periodicity of the array with respect to the case in Fig. 5, and it is noticed that the imaginary part

α_{z}

decreases for decreasing d. Again, this result is in agreement with what it has been previously shown in [20

A. Alù and N. Engheta, “Theory of linear chains of metamaterial/plasmonic particles as subdiffraction optical nanotransmission lines,” Phys. Rev. B 74(20), 205436 (2006). [CrossRef]

]. Also, notice that, similar to the case in Fig. 8(a), for small periodicity d, the frequency range for which the mode ’Improper 2’ is physical is wider than the other cases with larger periodicity.

Fig. 8 Parameter analysis as in Fig. 7, but for the ‘Improper 1’ leaky mode for L-pol in Fig. 5.

6. Conclusion

This work provides an extensive and comprehensive description of modal analysis in a representative case of a 1D-periodic array of metallic nanospheres in terms of bound/leaky, proper/improper, physical/nonphysical wave conditions. The subsets of physical and nonphysical waves have been identified and the general principle of their identification has been provided. A physical wave is one that can be excited by a realistic localized source, and it is classified in terms of path deformation in the complex wavenumber domain. Though the theory is involved, a simple Table summarizing the classification in terms of longitudinal and transverse wavenumbers, based on spectral theory, has been given. This is the first time such a complete physical characterization in a chain of plasmonic nanospheres has been presented, and to the authors‘ knowledge, this analysis for complex modes is the most comprehensive carried out so far. The formulation and classification are general. The high subwavelength regime and resonance of the array nanospheres are responsible for the interesting modal dispersion described in this paper; for example, it is very interesting to note that the T-pol physical backward leaky wave “Proper 1” in Fig. 4, increasing frequency, transitions into a physical forward leaky wave ‘Improper 1’. Knowledge of complex mode propagation is useful for the development of innovative sensors, field confinement and field enhancement structures, and highly directive antennas. Mode diagrams can be engineered by modifying the array parameters, thus allowing the design of waveguiding structures by using physical bound modes (such as ‘Proper 2’ in T-pol, and ‘Proper 1’ in L-pol) and radiating structure by using physical leaky modes (such as ‘Improper 1’ in L-pol).

Appendices

Appendix A: on the physical existence of modes excited by a localized source

Assume that the linear array of nanospheres is excited by a single dipolar source,

P^{i} (r') = {\hat{p}}^{i} δ (r' - r_{0})

r_{0} = x_{0} \hat{x} + y_{0} \hat{y} + z_{0} \hat{z}

, where

{\hat{p}}^{i}

, having units [Cm], is the impressed unit-amplitude electric dipole source. The total field produced by

P^{i}

and the chain, at r, can be represented as

E (r, r_{0}) = \frac{d}{2 π} \int_{- \infty}^{\infty} e^{i k_{z} z} e_{0}^{\infty} (x, y, r_{0}, k_{z}) d k_{z},

(8)

where

e_{0}^{\infty}

is the 0-th Floquet expansion term of the total field

E^{\infty} (r, r_{0}, k_{z}) = \sum_{p} e^{i k_{z, p} z} e_{p}^{\infty} (x, y, r_{0}, k_{z})

, generated by the chain excited by a phased array of point sources

p^{i, \infty} (r^{'}, k_{z}) = {\hat{p}}^{i} \sum_{n = - \infty}^{\infty} δ (r^{'} - r_{0} - n d \hat{z}) e^{i k_{z} n d}

(phased replicas of the real source

P^{i}

r_{0}

) according to the array scanning method (ASM), outlined in [41

,44

,45

], where a reader is referred to also for a direct expression of

e_{0}^{\infty}

. The integrand term

e_{0}^{\infty}

possesses spectral singularities: poles and branch points. Branch points are in correspondence of

k_{z b, p} = \pm k - 2 π p / d

such that

k_{ρ, p} = \sqrt{k^{2} - k_{z b, p}^{2}} = 0

. The term with p = 0 is within the fundamental BZ.

Figure 9 shows the path deformation onto the steepest descent paths (SDPs) and the capturable regions for poles representing physical modes in the complex

k_{z}

plane when an observer is along positive z. For an observer along negative z, the path on the real axis should be deformed onto the lower half of the complex

k_{z}

plane. These SDPs, in correspondence of branch points at

k_{z b, p} = k - 2 π p / d

, represent the “space wave” or continuous spectrum terms, and a precise discussion of these terms is not within the scope of this paper (more details in [41

,42

L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).

,45

]). The thick dark blue lines in Fig. 9 represent branch cuts, and when the integration path crosses a cut, a change in the Riemann sheet is performed (i.e., from proper to improper sheet and vice versa). The residues of the captured poles represent the physical waves (those that can be excited by a source

P^{i}

r_{0}

). Note that each mode is represented by an infinite set of poles (Floquet spatial harmonics) periodic along the

Re (k_{z})

axis. As can be noticed in Fig. 9, looking only at the spatial harmonics in the fundamental BZ, for an observer along the positive z-direction, only poles relative to backward proper (i.e., belonging to the top Riemann sheet) modes (leaky or bound) can be captured in the II quadrant (so called region R1) of the complex

k_{z}

plane. In the I quadrant, instead, only poles relative to forward proper bound (region R3, bound region, below the light line in the

k β plane

) and forward improper (i.e., belonging to the bottom Riemann sheet) leaky modes (region R2, radiating region, above the light line in the

k β plane

) can be captured (see for example Appendix A in [41

] or [42

L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).

] for more details regarding spectral singularities and spectral path deformation).

Fig. 9 Complex

k_{z}

wavenumber plane (assuming

k < π / d

). The field produced by a dipole source is represented by the integration path along the real axis. For an observer along the positive z axis, the path can be deformed in the upper half space. The poles encountered represent waves (each mode is represented by an infinite set of residues). The SDPs around the branch points at

k_{z} = k_{z b, p} = k - 2 π p / d

represent the “space wave” or continuous spectrum terms [41

,44

,45

,49

F. Capolino, D. R. Jackson, and D. R. Wilton, “Mode excitation from sources in two-dimensional EBG waveguides using the array scanning method,” IEEE Microw. Wirel. Compon. Lett. 15(2), 49–51 (2005). [CrossRef]

]. Physical modes (only the harmonic in the fundamental BZ, i.e., p = 0, is indicated) are those relative to poles captured in the three indicated regions (R1, R2 and R3). In each region, it is indicated the type of physical pole that can be captured (physical wave harmonic). PB = proper bound, PL = proper leaky, and IL = improper leaky wave.

The procedure outlined above leads to the following field representation for the total field at r excited by a source at

r_{0}

, as a sum of two wave species:

E (r, r_{0}) = \sum_{n} E_{n}^{mode} (r, r_{0}, k_{z}^{n}) + E_{t o t}^{sp} (r, r_{0}) .

(9)

Here,

E_{t o t}^{sp} (r, r_{0})

denotes a “spatial wave” collecting the contributions from all the SDPs, and

E_{n}^{mode} (r, r_{0}, k_{z}^{n})

corresponds to the n-th mode with wavenumber

k_{z}^{n}

[41

,44

,45

], which is expressed as in Eq. (7) and its characterization is the main purpose of this paper. Note that the physical property of the modes depends on the location of the corresponding pole and its capability of being captured. In summary physical waves are classified according to their wavenumbers as in Table 1.

Appendix B: Dyadic Ewald representation of the periodic GF for linear chains

The electric-field dyadic Green’s function

{\underline{\overset{⌣}{G}}}^{\infty} (k_{z})

for the phased periodic array of nanospheres in Eq. (4) is represented according to the Ewald method [32

–34

] as

{\underline{\overset{⌣}{G}}}^{\infty} (k_{z}) = {\underline{G}}_{spectral} (k_{z}) + {\underline{\overset{⌣}{G}}}_{spatial} (k_{z}),

(10)

where the spectral and the regularized spatial parts are

{\underline{G}}_{spectral} (k_{z}) = \frac{1}{ε_{h} ε_{0}} [k^{2} G_{spectral} (k_{z}) \underline{I} + {\underline{H}}_{spectral} (k_{z})],

(11)

{\underline{\overset{⌣}{G}}}_{spatial} (k_{z}) = \frac{1}{ε_{h} ε_{0}} [k^{2} {\overset{⌣}{G}}_{spatial} (k_{z}) \underline{I} + {\underline{\overset{⌣}{H}}}_{spatial} (k_{z})] .

(12)

The Ewald representation for the periodic scalar GF [32

,33

] (specialized for ρ = 0 and z = 0, and with regularized spatial term) is

G_{^{spectral}} (k_{z}) = \frac{1}{4 π d} \sum_{p = - \infty}^{\infty} E_{1}^{G} (- \frac{k_{ρ, p}^{2}}{4 E^{2}}),

(13)

{\overset{⌣}{G}}_{^{spatial}} (k_{z}) = \frac{1}{8 π} \sum_{\begin{array}{l} n = - \infty \\ (n \neq 0) \end{array}}^{\infty} \frac{e^{i k_{z} n d}}{r_{n}} f (r_{n}) + \frac{f^{'} (0) - 2 i k}{8 π},

(14)

with

f (r_{n}) = e^{- i k r_{n}} erfc (β_{-}) + e^{i k r_{n}} erfc (β_{+}), β_{\pm} = r_{n} E \pm i \frac{k}{2 E},

(15)

where

f^{'}

is the prime derivative of f with respect to its argument and E is the Ewald parameter as in [32

,33

], which is in general chosen as

E = \sqrt{π} / d

(other choices are discussed in [32

]). Furthermore,

E_{1}^{G} (ν) = E_{1} (ν) - i 2 π m

is the 1st-order generalized exponential integral, where

E_{1} (ν) = \int_{ν}^{\infty} t^{- 1} e^{- t} d t

- π \leq \arg ν \leq π

is the standard exponential integral [33

] and

m = 0 / \pm 1

is an integer, associated to the proper/improper character of the cylindrical p-th harmonic. In this particular implementation of the GF,

r_{n} = | n | d

is the distance between the point where we evaluate the GF (at

r = 0

) and the n-th dipolar source at

r_{n}^{'} = n d \hat{z}

, and

erfc (ν) = 2 π^{- 1 / 2} \int_{ν}^{\infty} e^{- t^{2}} d t

is the complementary error function [32

,33

]. Moreover

{\underline{H}}_{spectral} (k_{z}) = \frac{1}{4 π d} \sum_{p = - \infty}^{\infty} {\underline{F}}_{_{spectral, p}} (k_{z}),

(16)

{\underline{\overset{⌣}{H}}}_{spatial} (k_{z}) = \frac{1}{8 π} \sum_{\begin{array}{l} n = - \infty \\ (n \neq 0) \end{array}}^{\infty} e^{i k_{z} n d} {\underline{F}}_{_{spatial, n}} (k_{z}) + \frac{1}{8 π} \frac{f^{‴} (0) + 2 i k^{3}}{3} \underline{I},

(17)

where

{\underline{F}}_{_{spectral, p}} (k_{z}) = [- 2 E^{2} E_{2}^{G} (\frac{- k_{ρ, p}^{2}}{4 E^{2}}) (\underline{I} - \hat{z} \hat{z}) - E_{1}^{G} (\frac{- k_{ρ, p}^{2}}{4 E^{2}}) k_{z, p}^{2} \hat{z} \hat{z}],

(18)

{\underline{F}}_{_{spatial, n}} (k_{z}) = (\frac{f^{'} (r_{n})}{r_{n}^{2}} - \frac{f (r_{n})}{r_{n}^{3}}) \underline{I} + (\frac{f^{″} (r_{n})}{r_{n}} - \frac{3 f^{'} (r_{n})}{r_{n}^{2}} + \frac{3 f (r_{n})}{r_{n}^{3}}) \hat{z} \hat{z},

(19)

and

f^{″}

and

f^{‴}

are the second and third derivatives of f in Eq. (15) with respect to its argument, respectively, and

E_{2}^{G} (ν) = e^{- ν} - ν E_{1}^{G} (ν)

The discussion that follows is related to the choice of the m value in the generalized exponential integral

E_{1}^{G} (ν)

. We first recall that a purely spectral representation of the periodic Green’s function is in terms of outgoing cylindrical waves [32

] described by the Hankel function of the first kind

H_{0}^{(1)}

; therefore, to discern between proper and improper Riemann sheets of the

k_{z}

plane, we need to analyze the properties of the Hankel function with argument

k_{ρ, p} ρ

, describing the p-th spectral harmonic, whose asymptotic expansion for large argument is given by

H_{0}^{(1)} (k_{ρ, p} ρ) \sim \sqrt{\frac{2}{π k_{ρ, p} ρ}} e^{i (k_{ρ, p} ρ - \frac{π}{4})} \sim \frac{e^{i k_{ρ, p} ρ}}{\sqrt{k_{ρ, p} ρ}} (- π < arg k_{ρ, p} < 2 π),

(20)

as described in [40

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).

], page 364, where also the argument range

(- π, 2 π)

is provided.

From Eq. (20), considering

k_{ρ, p} = β_{ρ, p} + i α_{ρ, p}

and assuming that the general argument of the Hankel function, including argument phases in other Riemann sheets, is dictated by

k_{ρ, p} = ξ e^{m π i}

with

0 < arg ξ < π

, we identify the following mapping (the first two conditions have been previously given in [33

])

{\begin{matrix} - π < \arg (ξ e^{m π i}) < 0, m = - 1 Improper \\ 0 < \arg (ξ e^{m π i}) < π, m = 0 Proper \\ π < \arg (ξ e^{m π i}) < 2 π, m = + 1 Improper \end{matrix} .

(21)

The indexes m = 0 and −1 refer to the principal Riemann sheet of the Hankel function [40

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).

] with argument

k_{ρ, p} ρ = ξ e^{m π i} ρ

decaying (proper) and growing (improper) radially. The choice with m = +1 in Eq. (21) is relative to cylindrical waves growing radially (improper). Notice that when analyzing the FW harmonic with modal longitudinal wavenumber

k_{z}

in the fundamental BZ, i.e., with p = 0, improper modes are those with

- π < \arg k_{ρ, 0} < 0

π < \arg k_{ρ, 0} < 2 π

, and are thus computed by using m = −1 or m = +1 only for p = 0, as it has been described in [41

,45

] for planar 1D-periodic structures. Accordingly, when the wavenumber

k_{z}

is in the fundamental BZ, all other radial wavenumbers

k_{ρ, p}

with p ≠ 0 in Eq. (13) should be always such that

0 < \arg k_{ρ, p} < π

(proper choice), and thus evaluated with m = 0. Furthermore, for m values other than 0 and −1, according to the formula

π m < \arg k_{ρ, p} < π (m + 1)

, the

k_{z}

solutions are always found to be improper according to the analytic continuation of the Hankel function and its asymptotic behavior for large argument [40

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).

]

\begin{array}{l} H_{0}^{(1)} (k_{ρ, p} ρ) = H_{0}^{(1)} (ξ ρ) - 2 m J_{0} (ξ ρ) \\ \sim \frac{e^{i Re (ξ) ρ} e^{- Im (ξ) ρ}}{\sqrt{ξ ρ}} (1 - m) - m \frac{e^{- i Re (ξ) ρ} e^{Im (ξ) ρ}}{\sqrt{ξ ρ}}, \end{array}

(22)

where each index

m \neq 0, - 1

refers to non-principal Riemann sheets of the Hankel function.

The physical validity of these last improper solutions related to non-principal Riemann sheets of the Hankel function involved in cylindrical problems need further analysis, therefore they are not shown in this paper and will be the focus of future work. Nevertheless, here we present a detailed physical description of the modes previously reported (for example in [21

,26

M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]

]), distinguishing between proper and improper ones, forward and backward, bound and leaky, in a well defined framework. We would like to point out that the study of longitudinal complex modal propagation in open or non-uniform cylindrical domains, like a dielectric rod or a waveguide with an inner dielectric rod (simpler cases than the present one because they are not periodic), is an old subject [50

P. J. B. Clarricoats and K. R. Slinn, “Complex modes of propagation in dielectric-loaded circular waveguide,” Electron. Lett. 1(5), 145–146 (1965). [CrossRef]

–53

T. Rozzi, L. Pierantoni, and M. Farina, “General constraints on the propagation of complex waves in closed lossless isotropic waveguides,” IEEE Trans. Microw. Theory Tech. 46(5), 512–516 (1998). [CrossRef]

]. However, not much has been done in the study of modal solutions in cylindrical problems involving other Riemann sheets of the Hankel function. Several papers have been discussing the physical validity (individual excitability) of complex modes, even reaching contradicting conclusions [50

P. J. B. Clarricoats and K. R. Slinn, “Complex modes of propagation in dielectric-loaded circular waveguide,” Electron. Lett. 1(5), 145–146 (1965). [CrossRef]

–53

], as discussed in [54

R. Islam and G. V. Eleftheriades, “On the Independence of the Excitation of Complex Modes in Isotropic Structures,” IEEE Trans. Antenn. Propag. 58(5), 1567–1578 (2010). [CrossRef]

]. Recently, the authors of [54

R. Islam and G. V. Eleftheriades, “On the Independence of the Excitation of Complex Modes in Isotropic Structures,” IEEE Trans. Antenn. Propag. 58(5), 1567–1578 (2010). [CrossRef]

] claimed they sorted out old misconceptions regarding the physical validity of certain complex modes. Here we have provided a complete description regarding the principal Riemann sheet of the Hankel function but as already said further studies are required to analyze the solutions obtained from non-principal Riemann sheets. These are necessary to establish the wavenumber symmetry in the

k_{z}

complex plane, and indeed, according to [55

J. R. James, “Leaky waves on a dielectric rod,” Electron. Lett. 5(11), 252–254 (1969). [CrossRef]

], the complex conjugate symmetry of improper solutions can be obtained only considering solutions in adjacent Riemann sheets of the Hankel function (i.e., with different m values). What described in [55

J. R. James, “Leaky waves on a dielectric rod,” Electron. Lett. 5(11), 252–254 (1969). [CrossRef]

] has been confirmed here for the mode ‘Improper 1’ in T-pol in Figs. 3 and 4 (where we also found the wavenumber close to the complex conjugate solution correspondent to the mode ‘Improper 2’ by using m = +1, not shown), and for the mode ‘Improper 1’ in L-pol in the I and III quadrants of the complex

k_{z}

plane in Figs. 5 and 6 (where we also found the wavenumber close to the complex conjugate solution correspondent to the mode ‘Improper 1’ in the II and IV quadrants of the complex

k_{z}

plane by using m = +1, not shown). In this paper, in accordance with [55

J. R. James, “Leaky waves on a dielectric rod,” Electron. Lett. 5(11), 252–254 (1969). [CrossRef]

], symmetry of the proper solutions has been found in the proper plane (i.e., with m = 0).

Acknowledgments

Support by NSF Award #CMMI-1101074 is acknowledged. The authors would also like to thank Prof. D.R. Jackson, University of Houston, TX, Dr. R. Shore and Dr. A. Yaghjian, HANSCOM-AFRL, MA, for useful discussions.

References and links

1.	P. Alitalo, C. Simovski, A. Viitanen, and S. Tretyakov, “Near-field enhancement and subwavelength imaging in the optical region using a pair of two-dimensional arrays of metal nanospheres,” Phys. Rev. B 74(23), 235425 (2006). [CrossRef]
2.	S. Steshenko, F. Capolino, S. A. Tretyakov, and C. R. Simovski, “Super-Resolution and Near-Field Enhancement with Layers of Resonant Arrays of Nanoparticles,” in Applications of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 4.1.
3.	C. R. Simovski, A. J. Viitanen, and S. A. Tretyakov, “Sub-wavelength resolution in linear arrays of plasmonic particles,” J. Appl. Phys. 101(12), 123102 (2007). [CrossRef]
4.	C. Simovski, S. Tretyakov, and A. Viitanen, “Subwavelength imaging in a superlens of plasmon nanospheres,” Tech. Phys. Lett. 33(3), 264–266 (2007). [CrossRef]
5.	S. Steshenko, F. Capolino, P. Alitalo, and S. A. Tretyakov, “Effective model and investigation of the near-field enhancement and subwavelength imaging properties of multilayer arrays of plasmonic nanospheres,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 84(1), 016607 (2011). [CrossRef] [PubMed]
6.	A. F. Koenderink, “Plasmon nanoparticle array waveguides for single photon and single plasmon sources,” Nano Lett. 9(12), 4228–4233 (2009). [CrossRef] [PubMed]
7.	X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B 82(14), 144305 (2010). [CrossRef]
8.	J. Beermann, S. M. Novikov, K. Leosson, and S. I. Bozhevolnyi, “Surface enhanced Raman imaging: periodic arrays and individual metal nanoparticles,” Opt. Express 17(15), 12698–12705 (2009). [CrossRef] [PubMed]
9.	F. Liu, Z. Cao, C. Tang, L. Chen, and Z. Wang, “Ultrathin diamond-like carbon film coated silver nanoparticles-based substrates for surface-enhanced Raman spectroscopy,” ACS Nano 4(5), 2643–2648 (2010). [CrossRef] [PubMed]
10.	I. Firkowska, S. Giannona, J. A. Rojas-Chapana, K. Luecke, O. Brustle, and M. Giersig, “Biocompatible Nanomaterials and Nanodevices Promising for Biomedical Applications ” in Nanomaterials for Application in Medicine and Biology , M. Giersig, and G. B. Khomutov, eds. (Springer, Berlin, 2008), p. I.1.
11.	A. J. Haes and R. P. Van Duyne, “A nanoscale optical biosensor: sensitivity and selectivity of an approach based on the localized surface plasmon resonance spectroscopy of triangular silver nanoparticles,” J. Am. Chem. Soc. 124(35), 10596–10604 (2002). [CrossRef] [PubMed]
12.	I. E. Sendroiu and R. M. Corn, “Nanoparticle diffraction gratings for DNA detection on photopatterned glass substrates,” Biointerphases 3(3), FD23–FD29 (2008). [CrossRef] [PubMed]
13.	A. L. Fructos, S. Campione, F. Capolino, and F. Mesa, “Characterization of complex plasmonic modes in two-dimensional periodic arrays of metal nanospheres,” J. Opt. Soc. Am. B 28(6), 1446–1458 (2011). [CrossRef]
14.	S. A. Maier, P. G. Kik, and H. A. Atwater, “Optical pulse propagation in metal nanoparticle chain waveguides,” Phys. Rev. B 67(20), 205402 (2003). [CrossRef]
15.	S. Y. Park and D. Stroud, “Surface-plasmon dispersion relations in chains of metallic nanoparticles: an exact quasistatic calculation,” Phys. Rev. B 69(12), 125418 (2004). [CrossRef]
16.	R. A. Shore and A. D. Yaghjian, “Travelling electromagnetic waves on linear periodic arrays of lossless spheres,” Electron. Lett. 41(10), 578–580 (2005). [CrossRef]
17.	R. A. Shore and A. D. Yaghjian, ““Traveling electromagnetic waves on linear periodic arrays of lossless penetrable spheres,” IEICE Trans. Commun. E88B(6), 2346–2352 (2005). [CrossRef]
18.	D. S. Citrin, “Plasmon-polariton transport in metal-nanoparticle chains embedded in a gain medium,” Opt. Lett. 31(1), 98–100 (2006). [CrossRef] [PubMed]
19.	A. F. Koenderink and A. Polman, “Complex response and polariton-like dispersion splitting in periodic metal nanoparticle chains,” Phys. Rev. B 74(3), 033402 (2006). [CrossRef]
20.	A. Alù and N. Engheta, “Theory of linear chains of metamaterial/plasmonic particles as subdiffraction optical nanotransmission lines,” Phys. Rev. B 74(20), 205436 (2006). [CrossRef]
21.	R. A. Shore and A. D. Yaghjian, “Complex Waves on 1D, 2D, and 3D Periodic Arrays of Lossy and Lossless Magnetodielectric Spheres,” AFRL-RY-HS-TR-2010–0019, (Air Force Research Laboratory, Hanscom, MA 2010).
22.	K. B. Crozier, E. Togan, E. Simsek, and T. Yang, “Experimental measurement of the dispersion relations of the surface plasmon modes of metal nanoparticle chains,” Opt. Express 15(26), 17482–17493 (2007). [CrossRef] [PubMed]
23.	T. Yang and K. B. Crozier, “Dispersion and extinction of surface plasmons in an array of gold nanoparticle chains: influence of the air/glass interface,” Opt. Express 16(12), 8570–8580 (2008). [CrossRef] [PubMed]
24.	A. Alù and N. Engheta, “Effect of small random disorders and imperfections on the performance of arrays of plasmonic nanoparticles,” N. J. Phys. 12(1), 013015 (2010). [CrossRef]
25.	S. Campione and F. Capolino, “Linear and Planar Periodic Arrays of Metallic Nanospheres: Fabrication, Optical Properties and Applications,” in Selected Topics in Metamaterials and Photonic Crystals , A. Andreone, A. Cusano, A. Cutolo, and V. Galdi, eds. (World Scientific Publishing, 2011), pp. 141–194.
26.	M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B 27(8), 1576–1582 (2010). [CrossRef]
27.	C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).
28.	S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 8.1.
29.	D. E. Muller, “A Method for Solving Algebraic Equations Using an Automatic Computer,” Math. Tables Other Aids Comput. 10(56), 208–215 (1956). [CrossRef]
30.	W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipes: the art of scientific computing (Cambridge University Press, 2007).
31.	“IMSL Fortran Numerical Library,” (Visual Numerics Corporate Headquarters, 2500 Wilcrest Drive, Suite 200, Houston, TX), www.vni.com.
32.	F. Capolino, D. R. Wilton, and W. A. Johnson, “Efficient computation of the 3D Green's function for the Helmholtz operator for a linear array of point sources using the Ewald method,” J. Comput. Phys. 223(1), 250–261 (2007). [CrossRef]
33.	V. R. Komanduri, F. Capolino, D. R. Jackson, and D. R. Wilton, “Computation of the one-dimensional free-space periodic green's function for leaky waves using the ewald method,” in Proc. URSI Gen. Ass.(Chicago, IL, 2008).
34.	F. T. Celepcikay, D. R. Wilton, D. R. Jackson, and F. Capolino, “Choosing splitting parameters and summation limits in the numerical evaluation of 1-D and 2-D periodic Green's functions using the Ewald method,” Radio Sci. 43(6), RS6S01 (2008). [CrossRef]
35.	J. D. Jackson, Classical Electrodynamics (Wiley, 1998).
36.	V. A. Markel, V. N. Pustovit, S. V. Karpov, A. V. Obuschenko, V. S. Gerasimov, and I. L. Isaev, “Electromagnetic density of states and absorption of radiation by aggregates of nanospheres with multipole interactions,” Phys. Rev. B 70(5), 054202 (2004). [CrossRef]
37.	D. W. Mackowski, “Calculation of total cross sections of multiple-sphere clusters,” J. Opt. Soc. Am. A 11(11), 2851–2861 (1994). [CrossRef]
38.	R. Ruppin, “Evaluation of extended Maxwell-Garnett theories,” Opt. Commun. 182(4–6), 273–279 (2000). [CrossRef]
39.	W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B Condens. Matter 39(14), 9852–9858 (1989). [CrossRef] [PubMed]
40.	M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).
41.	F. Capolino, D. R. Jackson, and D. R. Wilton, “Field Representations in Periodic Artificial Materials Excited by a Source,” in Theory and Phenomena of Metamaterials , F. Capolino, ed. (CRC Press, 2009), p. 12.1.
42.	L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).
43.	P. Baccarelli, S. Paulotto, and C. Di Nallo, “Full-wave analysis of bound and leaky modes propagating along 2D periodic printed structures with arbitrary metallisation in the unit cell,” IET Proc. Microwaves, Antennas Propag. 1(1), 217–225 (2007). [CrossRef]
44.	F. Capolino, D. R. Jackson, D. R. Wilton, and L. B. Felsen, “Comparison of methods for calculating the field excited by a dipole near a 2-D periodic material,” IEEE Trans. Antenn. Propag. 55(6), 1644–1655 (2007). [CrossRef]
45.	F. Capolino, D. R. Jackson, and D. R. Wilton, “Fundamental properties of the field at the interface between air and a periodic artificial material excited by a line source,” IEEE Trans. Antenn. Propag. 53(1), 91–99 (2005). [CrossRef]
46.	A. Alù, A. Salandrino, and N. Engheta, “Negative effective permeability and left-handed materials at optical frequencies,” Opt. Express 14(4), 1557–1567 (2006). [CrossRef] [PubMed]
47.	I. El-Kady, M. M. Sigalas, R. Biswas, K. M. Ho, and C. M. Soukoulis, “Metallic photonic crystals at optical wavelengths,” Phys. Rev. B 62(23), 15299–15302 (2000). [CrossRef]
48.	A. A. Oliner and D. R. Jackson, “Leaky-wave antennas,” in Antenna Engineering Handbook , J. Volakis, ed. (McGraw Hill, 2007), p. 11.1.
49.	F. Capolino, D. R. Jackson, and D. R. Wilton, “Mode excitation from sources in two-dimensional EBG waveguides using the array scanning method,” IEEE Microw. Wirel. Compon. Lett. 15(2), 49–51 (2005). [CrossRef]
50.	P. J. B. Clarricoats and K. R. Slinn, “Complex modes of propagation in dielectric-loaded circular waveguide,” Electron. Lett. 1(5), 145–146 (1965). [CrossRef]
51.	J. D. Rhodes, “General constraints on propagation characteristics of electromagnetic waves in uniform inhomogeneous waveguides,” in Proc. Inst. Electr. Eng . 118, 849–856 (1971).
52.	T. F. Jablonski, “Complex modes in open lossless dielectric waveguides,” J. Opt. Soc. Am. A 11(4), 1272–1282 (1994). [CrossRef]
53.	T. Rozzi, L. Pierantoni, and M. Farina, “General constraints on the propagation of complex waves in closed lossless isotropic waveguides,” IEEE Trans. Microw. Theory Tech. 46(5), 512–516 (1998). [CrossRef]
54.	R. Islam and G. V. Eleftheriades, “On the Independence of the Excitation of Complex Modes in Isotropic Structures,” IEEE Trans. Antenn. Propag. 58(5), 1567–1578 (2010). [CrossRef]
55.	J. R. James, “Leaky waves on a dielectric rod,” Electron. Lett. 5(11), 252–254 (1969). [CrossRef]

OCIS Codes
(260.2110) Physical optics : Electromagnetic optics
(160.3918) Materials : Metamaterials
(250.5403) Optoelectronics : Plasmonics

ToC Category:
Metamaterials

History
Original Manuscript: June 3, 2011
Revised Manuscript: August 13, 2011
Manuscript Accepted: August 17, 2011
Published: September 6, 2011

Citation
Salvatore Campione, Sergiy Steshenko, and Filippo Capolino, "Complex bound and leaky modes in chains of plasmonic nanospheres," Opt. Express 19, 18345-18363 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-19-18345

Sort: Author | Year | Journal | Reset

References

P. Alitalo, C. Simovski, A. Viitanen, and S. Tretyakov, “Near-field enhancement and subwavelength imaging in the optical region using a pair of two-dimensional arrays of metal nanospheres,” Phys. Rev. B74(23), 235425 (2006). [CrossRef]
S. Steshenko, F. Capolino, S. A. Tretyakov, and C. R. Simovski, “Super-Resolution and Near-Field Enhancement with Layers of Resonant Arrays of Nanoparticles,” in Applications of Metamaterials, F. Capolino, ed. (CRC Press, 2009), p. 4.1.
C. R. Simovski, A. J. Viitanen, and S. A. Tretyakov, “Sub-wavelength resolution in linear arrays of plasmonic particles,” J. Appl. Phys.101(12), 123102 (2007). [CrossRef]
C. Simovski, S. Tretyakov, and A. Viitanen, “Subwavelength imaging in a superlens of plasmon nanospheres,” Tech. Phys. Lett.33(3), 264–266 (2007). [CrossRef]
S. Steshenko, F. Capolino, P. Alitalo, and S. A. Tretyakov, “Effective model and investigation of the near-field enhancement and subwavelength imaging properties of multilayer arrays of plasmonic nanospheres,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys.84(1), 016607 (2011). [CrossRef] [PubMed]
A. F. Koenderink, “Plasmon nanoparticle array waveguides for single photon and single plasmon sources,” Nano Lett.9(12), 4228–4233 (2009). [CrossRef] [PubMed]
X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B82(14), 144305 (2010). [CrossRef]
J. Beermann, S. M. Novikov, K. Leosson, and S. I. Bozhevolnyi, “Surface enhanced Raman imaging: periodic arrays and individual metal nanoparticles,” Opt. Express17(15), 12698–12705 (2009). [CrossRef] [PubMed]
F. Liu, Z. Cao, C. Tang, L. Chen, and Z. Wang, “Ultrathin diamond-like carbon film coated silver nanoparticles-based substrates for surface-enhanced Raman spectroscopy,” ACS Nano4(5), 2643–2648 (2010). [CrossRef] [PubMed]
I. Firkowska, S. Giannona, J. A. Rojas-Chapana, K. Luecke, O. Brustle, and M. Giersig, “Biocompatible Nanomaterials and Nanodevices Promising for Biomedical Applications ” in Nanomaterials for Application in Medicine and Biology, M. Giersig, and G. B. Khomutov, eds. (Springer, Berlin, 2008), p. I.1.
A. J. Haes and R. P. Van Duyne, “A nanoscale optical biosensor: sensitivity and selectivity of an approach based on the localized surface plasmon resonance spectroscopy of triangular silver nanoparticles,” J. Am. Chem. Soc.124(35), 10596–10604 (2002). [CrossRef] [PubMed]
I. E. Sendroiu and R. M. Corn, “Nanoparticle diffraction gratings for DNA detection on photopatterned glass substrates,” Biointerphases3(3), FD23–FD29 (2008). [CrossRef] [PubMed]
A. L. Fructos, S. Campione, F. Capolino, and F. Mesa, “Characterization of complex plasmonic modes in two-dimensional periodic arrays of metal nanospheres,” J. Opt. Soc. Am. B28(6), 1446–1458 (2011). [CrossRef]
S. A. Maier, P. G. Kik, and H. A. Atwater, “Optical pulse propagation in metal nanoparticle chain waveguides,” Phys. Rev. B67(20), 205402 (2003). [CrossRef]
S. Y. Park and D. Stroud, “Surface-plasmon dispersion relations in chains of metallic nanoparticles: an exact quasistatic calculation,” Phys. Rev. B69(12), 125418 (2004). [CrossRef]
R. A. Shore and A. D. Yaghjian, “Travelling electromagnetic waves on linear periodic arrays of lossless spheres,” Electron. Lett.41(10), 578–580 (2005). [CrossRef]
R. A. Shore and A. D. Yaghjian, ““Traveling electromagnetic waves on linear periodic arrays of lossless penetrable spheres,” IEICE Trans. Commun.E88B(6), 2346–2352 (2005). [CrossRef]
D. S. Citrin, “Plasmon-polariton transport in metal-nanoparticle chains embedded in a gain medium,” Opt. Lett.31(1), 98–100 (2006). [CrossRef] [PubMed]
A. F. Koenderink and A. Polman, “Complex response and polariton-like dispersion splitting in periodic metal nanoparticle chains,” Phys. Rev. B74(3), 033402 (2006). [CrossRef]
A. Alù and N. Engheta, “Theory of linear chains of metamaterial/plasmonic particles as subdiffraction optical nanotransmission lines,” Phys. Rev. B74(20), 205436 (2006). [CrossRef]
R. A. Shore and A. D. Yaghjian, “Complex Waves on 1D, 2D, and 3D Periodic Arrays of Lossy and Lossless Magnetodielectric Spheres,” AFRL-RY-HS-TR-2010–0019, (Air Force Research Laboratory, Hanscom, MA 2010).
K. B. Crozier, E. Togan, E. Simsek, and T. Yang, “Experimental measurement of the dispersion relations of the surface plasmon modes of metal nanoparticle chains,” Opt. Express15(26), 17482–17493 (2007). [CrossRef] [PubMed]
T. Yang and K. B. Crozier, “Dispersion and extinction of surface plasmons in an array of gold nanoparticle chains: influence of the air/glass interface,” Opt. Express16(12), 8570–8580 (2008). [CrossRef] [PubMed]
A. Alù and N. Engheta, “Effect of small random disorders and imperfections on the performance of arrays of plasmonic nanoparticles,” N. J. Phys.12(1), 013015 (2010). [CrossRef]
S. Campione and F. Capolino, “Linear and Planar Periodic Arrays of Metallic Nanospheres: Fabrication, Optical Properties and Applications,” in Selected Topics in Metamaterials and Photonic Crystals, A. Andreone, A. Cusano, A. Cutolo, and V. Galdi, eds. (World Scientific Publishing, 2011), pp. 141–194.
M. Conforti and M. Guasoni, “Dispersive properties of linear chains of lossy metal nanoparticles,” J. Opt. Soc. Am. B27(8), 1576–1582 (2010). [CrossRef]
C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).
S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials, F. Capolino, ed. (CRC Press, 2009), p. 8.1.
D. E. Muller, “A Method for Solving Algebraic Equations Using an Automatic Computer,” Math. Tables Other Aids Comput.10(56), 208–215 (1956). [CrossRef]
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipes: the art of scientific computing (Cambridge University Press, 2007).
“IMSL Fortran Numerical Library,” (Visual Numerics Corporate Headquarters, 2500 Wilcrest Drive, Suite 200, Houston, TX), www.vni.com .
F. Capolino, D. R. Wilton, and W. A. Johnson, “Efficient computation of the 3D Green's function for the Helmholtz operator for a linear array of point sources using the Ewald method,” J. Comput. Phys.223(1), 250–261 (2007). [CrossRef]
V. R. Komanduri, F. Capolino, D. R. Jackson, and D. R. Wilton, “Computation of the one-dimensional free-space periodic green's function for leaky waves using the ewald method,” in Proc. URSI Gen. Ass.(Chicago, IL, 2008).
F. T. Celepcikay, D. R. Wilton, D. R. Jackson, and F. Capolino, “Choosing splitting parameters and summation limits in the numerical evaluation of 1-D and 2-D periodic Green's functions using the Ewald method,” Radio Sci.43(6), RS6S01 (2008). [CrossRef]
J. D. Jackson, Classical Electrodynamics (Wiley, 1998).
V. A. Markel, V. N. Pustovit, S. V. Karpov, A. V. Obuschenko, V. S. Gerasimov, and I. L. Isaev, “Electromagnetic density of states and absorption of radiation by aggregates of nanospheres with multipole interactions,” Phys. Rev. B70(5), 054202 (2004). [CrossRef]
D. W. Mackowski, “Calculation of total cross sections of multiple-sphere clusters,” J. Opt. Soc. Am. A11(11), 2851–2861 (1994). [CrossRef]
R. Ruppin, “Evaluation of extended Maxwell-Garnett theories,” Opt. Commun.182(4–6), 273–279 (2000). [CrossRef]
W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B Condens. Matter39(14), 9852–9858 (1989). [CrossRef] [PubMed]
M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).
F. Capolino, D. R. Jackson, and D. R. Wilton, “Field Representations in Periodic Artificial Materials Excited by a Source,” in Theory and Phenomena of Metamaterials, F. Capolino, ed. (CRC Press, 2009), p. 12.1.
L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).
P. Baccarelli, S. Paulotto, and C. Di Nallo, “Full-wave analysis of bound and leaky modes propagating along 2D periodic printed structures with arbitrary metallisation in the unit cell,” IET Proc. Microwaves, Antennas Propag.1(1), 217–225 (2007). [CrossRef]
F. Capolino, D. R. Jackson, D. R. Wilton, and L. B. Felsen, “Comparison of methods for calculating the field excited by a dipole near a 2-D periodic material,” IEEE Trans. Antenn. Propag.55(6), 1644–1655 (2007). [CrossRef]
F. Capolino, D. R. Jackson, and D. R. Wilton, “Fundamental properties of the field at the interface between air and a periodic artificial material excited by a line source,” IEEE Trans. Antenn. Propag.53(1), 91–99 (2005). [CrossRef]
A. Alù, A. Salandrino, and N. Engheta, “Negative effective permeability and left-handed materials at optical frequencies,” Opt. Express14(4), 1557–1567 (2006). [CrossRef] [PubMed]
I. El-Kady, M. M. Sigalas, R. Biswas, K. M. Ho, and C. M. Soukoulis, “Metallic photonic crystals at optical wavelengths,” Phys. Rev. B62(23), 15299–15302 (2000). [CrossRef]
A. A. Oliner and D. R. Jackson, “Leaky-wave antennas,” in Antenna Engineering Handbook, J. Volakis, ed. (McGraw Hill, 2007), p. 11.1.
F. Capolino, D. R. Jackson, and D. R. Wilton, “Mode excitation from sources in two-dimensional EBG waveguides using the array scanning method,” IEEE Microw. Wirel. Compon. Lett.15(2), 49–51 (2005). [CrossRef]
P. J. B. Clarricoats and K. R. Slinn, “Complex modes of propagation in dielectric-loaded circular waveguide,” Electron. Lett.1(5), 145–146 (1965). [CrossRef]
J. D. Rhodes, “General constraints on propagation characteristics of electromagnetic waves in uniform inhomogeneous waveguides,” in Proc. Inst. Electr. Eng. 118, 849–856 (1971).
T. F. Jablonski, “Complex modes in open lossless dielectric waveguides,” J. Opt. Soc. Am. A11(4), 1272–1282 (1994). [CrossRef]
T. Rozzi, L. Pierantoni, and M. Farina, “General constraints on the propagation of complex waves in closed lossless isotropic waveguides,” IEEE Trans. Microw. Theory Tech.46(5), 512–516 (1998). [CrossRef]
R. Islam and G. V. Eleftheriades, “On the Independence of the Excitation of Complex Modes in Isotropic Structures,” IEEE Trans. Antenn. Propag.58(5), 1567–1578 (2010). [CrossRef]
J. R. James, “Leaky waves on a dielectric rod,” Electron. Lett.5(11), 252–254 (1969). [CrossRef]

ACS Nano

F. Liu, Z. Cao, C. Tang, L. Chen, and Z. Wang, “Ultrathin diamond-like carbon film coated silver nanoparticles-based substrates for surface-enhanced Raman spectroscopy,” ACS Nano4(5), 2643–2648 (2010). [CrossRef] [PubMed]

Biointerphases

I. E. Sendroiu and R. M. Corn, “Nanoparticle diffraction gratings for DNA detection on photopatterned glass substrates,” Biointerphases3(3), FD23–FD29 (2008). [CrossRef] [PubMed]

Electron. Lett.

R. A. Shore and A. D. Yaghjian, “Travelling electromagnetic waves on linear periodic arrays of lossless spheres,” Electron. Lett.41(10), 578–580 (2005). [CrossRef]
P. J. B. Clarricoats and K. R. Slinn, “Complex modes of propagation in dielectric-loaded circular waveguide,” Electron. Lett.1(5), 145–146 (1965). [CrossRef]
J. R. James, “Leaky waves on a dielectric rod,” Electron. Lett.5(11), 252–254 (1969). [CrossRef]

IEEE Microw. Wirel. Compon. Lett.

F. Capolino, D. R. Jackson, and D. R. Wilton, “Mode excitation from sources in two-dimensional EBG waveguides using the array scanning method,” IEEE Microw. Wirel. Compon. Lett.15(2), 49–51 (2005). [CrossRef]

IEEE Trans. Antenn. Propag.

R. Islam and G. V. Eleftheriades, “On the Independence of the Excitation of Complex Modes in Isotropic Structures,” IEEE Trans. Antenn. Propag.58(5), 1567–1578 (2010). [CrossRef]
F. Capolino, D. R. Jackson, D. R. Wilton, and L. B. Felsen, “Comparison of methods for calculating the field excited by a dipole near a 2-D periodic material,” IEEE Trans. Antenn. Propag.55(6), 1644–1655 (2007). [CrossRef]
F. Capolino, D. R. Jackson, and D. R. Wilton, “Fundamental properties of the field at the interface between air and a periodic artificial material excited by a line source,” IEEE Trans. Antenn. Propag.53(1), 91–99 (2005). [CrossRef]

IEEE Trans. Microw. Theory Tech.

T. Rozzi, L. Pierantoni, and M. Farina, “General constraints on the propagation of complex waves in closed lossless isotropic waveguides,” IEEE Trans. Microw. Theory Tech.46(5), 512–516 (1998). [CrossRef]

IEICE Trans. Commun.

R. A. Shore and A. D. Yaghjian, ““Traveling electromagnetic waves on linear periodic arrays of lossless penetrable spheres,” IEICE Trans. Commun.E88B(6), 2346–2352 (2005). [CrossRef]

IET Proc. Microwaves, Antennas Propag.

P. Baccarelli, S. Paulotto, and C. Di Nallo, “Full-wave analysis of bound and leaky modes propagating along 2D periodic printed structures with arbitrary metallisation in the unit cell,” IET Proc. Microwaves, Antennas Propag.1(1), 217–225 (2007). [CrossRef]

J. Am. Chem. Soc.

A. J. Haes and R. P. Van Duyne, “A nanoscale optical biosensor: sensitivity and selectivity of an approach based on the localized surface plasmon resonance spectroscopy of triangular silver nanoparticles,” J. Am. Chem. Soc.124(35), 10596–10604 (2002). [CrossRef] [PubMed]

J. Appl. Phys.

C. R. Simovski, A. J. Viitanen, and S. A. Tretyakov, “Sub-wavelength resolution in linear arrays of plasmonic particles,” J. Appl. Phys.101(12), 123102 (2007). [CrossRef]

J. Comput. Phys.

F. Capolino, D. R. Wilton, and W. A. Johnson, “Efficient computation of the 3D Green's function for the Helmholtz operator for a linear array of point sources using the Ewald method,” J. Comput. Phys.223(1), 250–261 (2007). [CrossRef]

J. Opt. Soc. Am. A

J. Opt. Soc. Am. B

Math. Tables Other Aids Comput.

D. E. Muller, “A Method for Solving Algebraic Equations Using an Automatic Computer,” Math. Tables Other Aids Comput.10(56), 208–215 (1956). [CrossRef]

N. J. Phys.

A. Alù and N. Engheta, “Effect of small random disorders and imperfections on the performance of arrays of plasmonic nanoparticles,” N. J. Phys.12(1), 013015 (2010). [CrossRef]

Nano Lett.

A. F. Koenderink, “Plasmon nanoparticle array waveguides for single photon and single plasmon sources,” Nano Lett.9(12), 4228–4233 (2009). [CrossRef] [PubMed]

Opt. Commun.

R. Ruppin, “Evaluation of extended Maxwell-Garnett theories,” Opt. Commun.182(4–6), 273–279 (2000). [CrossRef]

Opt. Express

Opt. Lett.

D. S. Citrin, “Plasmon-polariton transport in metal-nanoparticle chains embedded in a gain medium,” Opt. Lett.31(1), 98–100 (2006). [CrossRef] [PubMed]

Phys. Rev. B

A. F. Koenderink and A. Polman, “Complex response and polariton-like dispersion splitting in periodic metal nanoparticle chains,” Phys. Rev. B74(3), 033402 (2006). [CrossRef]
A. Alù and N. Engheta, “Theory of linear chains of metamaterial/plasmonic particles as subdiffraction optical nanotransmission lines,” Phys. Rev. B74(20), 205436 (2006). [CrossRef]
S. A. Maier, P. G. Kik, and H. A. Atwater, “Optical pulse propagation in metal nanoparticle chain waveguides,” Phys. Rev. B67(20), 205402 (2003). [CrossRef]
S. Y. Park and D. Stroud, “Surface-plasmon dispersion relations in chains of metallic nanoparticles: an exact quasistatic calculation,” Phys. Rev. B69(12), 125418 (2004). [CrossRef]
X. Liu and A. Alu, “Subwavelength leaky-wave optical nanoantennas: directive radiation from linear arrays of plasmonic nanoparticles,” Phys. Rev. B82(14), 144305 (2010). [CrossRef]
P. Alitalo, C. Simovski, A. Viitanen, and S. Tretyakov, “Near-field enhancement and subwavelength imaging in the optical region using a pair of two-dimensional arrays of metal nanospheres,” Phys. Rev. B74(23), 235425 (2006). [CrossRef]
V. A. Markel, V. N. Pustovit, S. V. Karpov, A. V. Obuschenko, V. S. Gerasimov, and I. L. Isaev, “Electromagnetic density of states and absorption of radiation by aggregates of nanospheres with multipole interactions,” Phys. Rev. B70(5), 054202 (2004). [CrossRef]
I. El-Kady, M. M. Sigalas, R. Biswas, K. M. Ho, and C. M. Soukoulis, “Metallic photonic crystals at optical wavelengths,” Phys. Rev. B62(23), 15299–15302 (2000). [CrossRef]

Phys. Rev. B Condens. Matter

W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B Condens. Matter39(14), 9852–9858 (1989). [CrossRef] [PubMed]

Phys. Rev. E Stat. Nonlin. Soft Matter Phys.

S. Steshenko, F. Capolino, P. Alitalo, and S. A. Tretyakov, “Effective model and investigation of the near-field enhancement and subwavelength imaging properties of multilayer arrays of plasmonic nanospheres,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys.84(1), 016607 (2011). [CrossRef] [PubMed]

Radio Sci.

F. T. Celepcikay, D. R. Wilton, D. R. Jackson, and F. Capolino, “Choosing splitting parameters and summation limits in the numerical evaluation of 1-D and 2-D periodic Green's functions using the Ewald method,” Radio Sci.43(6), RS6S01 (2008). [CrossRef]

Tech. Phys. Lett.

C. Simovski, S. Tretyakov, and A. Viitanen, “Subwavelength imaging in a superlens of plasmon nanospheres,” Tech. Phys. Lett.33(3), 264–266 (2007). [CrossRef]

Other

S. Steshenko, F. Capolino, S. A. Tretyakov, and C. R. Simovski, “Super-Resolution and Near-Field Enhancement with Layers of Resonant Arrays of Nanoparticles,” in Applications of Metamaterials, F. Capolino, ed. (CRC Press, 2009), p. 4.1.
R. A. Shore and A. D. Yaghjian, “Complex Waves on 1D, 2D, and 3D Periodic Arrays of Lossy and Lossless Magnetodielectric Spheres,” AFRL-RY-HS-TR-2010–0019, (Air Force Research Laboratory, Hanscom, MA 2010).
J. D. Jackson, Classical Electrodynamics (Wiley, 1998).
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical recipes: the art of scientific computing (Cambridge University Press, 2007).
“IMSL Fortran Numerical Library,” (Visual Numerics Corporate Headquarters, 2500 Wilcrest Drive, Suite 200, Houston, TX), www.vni.com .
V. R. Komanduri, F. Capolino, D. R. Jackson, and D. R. Wilton, “Computation of the one-dimensional free-space periodic green's function for leaky waves using the ewald method,” in Proc. URSI Gen. Ass.(Chicago, IL, 2008).
I. Firkowska, S. Giannona, J. A. Rojas-Chapana, K. Luecke, O. Brustle, and M. Giersig, “Biocompatible Nanomaterials and Nanodevices Promising for Biomedical Applications ” in Nanomaterials for Application in Medicine and Biology, M. Giersig, and G. B. Khomutov, eds. (Springer, Berlin, 2008), p. I.1.
S. Campione and F. Capolino, “Linear and Planar Periodic Arrays of Metallic Nanospheres: Fabrication, Optical Properties and Applications,” in Selected Topics in Metamaterials and Photonic Crystals, A. Andreone, A. Cusano, A. Cutolo, and V. Galdi, eds. (World Scientific Publishing, 2011), pp. 141–194.
C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).
S. Steshenko and F. Capolino, “Single Dipole Approximation for Modeling Collections of Nanoscatterers,” in Theory and Phenomena of Metamaterials, F. Capolino, ed. (CRC Press, 2009), p. 8.1.
M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965).
F. Capolino, D. R. Jackson, and D. R. Wilton, “Field Representations in Periodic Artificial Materials Excited by a Source,” in Theory and Phenomena of Metamaterials, F. Capolino, ed. (CRC Press, 2009), p. 12.1.
L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves (IEEE Press, 1994).
A. A. Oliner and D. R. Jackson, “Leaky-wave antennas,” in Antenna Engineering Handbook, J. Volakis, ed. (McGraw Hill, 2007), p. 11.1.
J. D. Rhodes, “General constraints on propagation characteristics of electromagnetic waves in uniform inhomogeneous waveguides,” in Proc. Inst. Electr. Eng. 118, 849–856 (1971).

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.

Click here to see a list of articles that cite this paper