Layer-multiple-scattering method for photonic structures of general scatterers based on a discrete-dipole approximation/t-matrix point-matching method
|
|
- Iris Page
- 6 years ago
- Views:
Transcription
1 V. Yannopapas Vol. 31, No. 3 / March 2014 / J. Opt. Soc. Am. B 631 Layer-multiple-scattering method for photonic structures of general scatterers based on a discrete-dipole approximation/t-matrix point-matching method Vassilios Yannopapas Department of Physics, National Technical University of Athens, GR Athens, Greece (vyannop@mail.ntua.gr) Received October 23, 2013; revised December 5, 2013; accepted December 30, 2013; posted January 23, 2014 (Doc. ID ); published February 27, 2014 We present a hybrid discrete-dipole approximation (DDA)/layer-multiple-scattering (LMS) method for treating photonic structures consisting of general scatterers. The present method is a major extension of the existing LMS formalism and code that combine the merits of both the DDA and LMS technique. Namely, the new hybrid DDA/LMS technique treats, in principle, scatterers of general (nonspherical) shape that might be anisotropic and/ or inhomogeneous, thanks to the DDA component while, at the same time, it incorporates theoretical tools provided by the LMS method, such as the doubling-layer process and the complex frequency band structure that are not met in contemporary electromagnetic solvers. The merging of both techniques is accomplished via a pointmatching module that provides the scattering T-matrix of an arbitrary scatterer. To demonstrate the applicability of the new method, we study the optical properties of 2D and 3D plasmonic lattices of gold nanocubes Optical Society of America OCIS codes: ( ) Multiple scattering; ( ) Metamaterials; ( ) Photonic crystals; ( ) Surface plasmons; ( ) Nanomaterials INTRODUCTION Man-made periodic structures, depending on their functionality and scope, fall within three main categories: photonic crystals, plasmonic materials, and metamaterials. Photonic crystals possess an absolute bandgap, i.e., a frequency region where no electromagnetic (EM) Bloch modes are allowed as a result of destructive interference, allowing for the passive control of light emission and flow within their volume [1]. Plasmonic materials and metamaterials operate both in the subwavelength regime, i.e., around the center of the Brillouin zone. Plasmonic materials are mainly used for boosting the electric field in small volumes and tailoring light absorption [2]. Metamaterials have more unconventional features, such as negative refractive index and/or magnetic permeability with application in optical microscopy, imaging, and cloaking [3]. The fabrication of photonic structures is made either by laser- and electron-beam lithography or by self-assembly. By lithography, one usually realizes planar structures with repeatable units characterized by sharp edges, such as rods, pillars, disks, etc., whereas by self-assembly the repeated units are of spherical or nearly spherical shape. The modeling tools for simulating photonic structures depend on the fabrication technique, i.e., lithography or self-assembly. For lithographic designs purely numerical EM solvers are employed, such as the transfer-matrix method [4], the finite element method (FEM) [5], the finite-difference time-domain (FDTD) method [6], the pseudo-spectral time-domain method (PSTD) [7,8] the Fourier modal method [9], and the discrete-dipole approximation (DDA) [10]. In principle, such methods can be also used for self-assembled structures, such as arrays of spheres requiring, however, huge amounts of memory and computing time. The state-of-the-art technique in simulating two- and three-dimensional arrays of spheres is the layermultiple scattering (LMS) method [11 13]. There, the EM field within a given 2D lattice (plane) of spheres is obtained via the multiple-scattering technique. In between a pair of two consecutive planes of spheres, the EM field is written in a plane wave basis where the multiple-scattering between the planes is fully taken into account. A slab of macroscopic thickness, i.e., containing several hundreds of planes of spheres, can be easily built up via a doubling-layer process [11 13]. By imposing boundary conditions along the growth direction of the crystal slab, one can also calculate the complex frequency band structure [11 13]. In this work, we aim at merging the DDA with the LMS method to a common theoretical framework and computing code (a hybrid DDA/LMS method). By doing so, we combine the merits of both techniques so that the new method is capable of treating photonic structures of general scatterers. Namely, the new hybrid DDA/LMS technique treats, in principle, scatterers of general (nonspherical) shape that might be anisotropic and/or inhomogeneous, thanks to the DDA component while, at the same time, it incorporates theoretical tools provided by the LMS method, such as the doubling-layer process and the complex frequency band structure, which are not met in contemporary EM solvers. The proposed hybrid DDA/LMS method modularizes the computational process into rapidly converging programming units which addresses the shortcomings of the existing simulators. Currently, the LMS method has been extended so as to treat nonspherical scatterers with axial symmetry [14,15], anisotropic spherical scatterers [16,17] as well as clusters of spherical scatterers [18]. The extension of the LMS to axisymmetric scatterers [14,15] has been achieved by the boundary element method /14/ $15.00/ Optical Society of America
2 632 J. Opt. Soc. Am. B / Vol. 31, No. 3 / March 2014 V. Yannopapas (BEM), which is faster than the DDA in treating homogeneous, isotropic, nonspherical scatterers [19]. However, it lacks the generality of the DDA, which treats anisotropic and inhomogeneous scatterers on equal footing. Originally, the DDA had been mainly applied to individual low-index particles [10,20 22] and 2D arrays of such [23]. However, the rapid increase in computing power, allowed for the study of individual high-index scatterers such as silicon [24,25] and metallic nanoparticles [26] as well as their corresponding 2D arrays [27 29]. Extension of the DDA to 3D arrays is still lacking due to the expected computational overhead. From this point of view, the hybrid DDA/LMS proposed here can be seen as an extension of the DDA to 3D crystals although it possesses much more functionalities than a mere latticetype DDA. 2. THEORY A. Outline of the Hybrid DDA/LMS Method The DDA is among the well-established methods for solving EM scattering from an individual or multiple scatterers (targets) of general shape and moderate size. The central idea of the DDA is the replacement of the actual scatterer by a set of point dipoles interacting with each other and an incident EM field, giving rise to a linear system of equations that provides the polarizations at each point dipole. Due to its conceptual simplicity and transparency, the DDA is very efficient in treating homogeneous or inhomogeneous scatterers with elaborate shape and/or anisotropic EM response. The LMS method is a rigorous theory for solving the problem of scattering of EM waves by a slab containing a finite number of 2D planes of spherical particles and/or homogeneous slabs. The LMS method simulates an actual experimental setup since it provides the transmission, reflection and/or absorption of light incident (with arbitrary angle and polarization) on a finite slab made from several planes of spheres and/or homogeneous slabs. A central quantity in the LMS method is the so-called scattering T-matrix. The latter connects the Fourier coefficients in the spherical-wave expansion of an EM wave scattered off of a finite object with those of the incident wave. The aim of the present paper is to generalize the LMS method along the above lines with hope of superseding the current state-of-start in electrodynamic simulations of metamaterials and photonic crystals. The generalization can be realized by the combination of the DDA and the LMS method within the framework of the point-matching method. According to this (see Fig. 1), the first step is the application of the DDA, wherein the actual scatterer [Fig. 1(a)] is substituted by a finite grid of interacting point dipoles [Fig. 1(b)], which is illuminated by an externally incident EM field. Solution of the corresponding DDA linear system of equations provides the dipole moments at each grid point, as well as the scattered electric field at a given point in space. The scattered electric field calculated by DDA is matched to that written in terms of vector spherical harmonics. The matching of the fields obtained by both methods takes place at points around the spherical scatterer, e.g., on a spherical surface surrounding the scatterer [Fig. 1(c)]. Having ensured the convergence in the calculation of the T-matrix, the latter is embedded into the LMS programming code [Fig. 1(d)], which provides the EM response of the metamaterial containing the scatterers in question [Fig. 1(e)]. B. Discrete-Dipole Approximation Next, we study the optical response of a generally anisotropic scatterer of arbitrary shape using the DDA [10,20 22]. The scattering object is considered as an array of point dipoles (i 1; ;N), each of which is located at the position r i and corresponds to a dipole moment P i and a (positiondependent) polarizability tensor ~α i. The above quantities are connected by where E i is the electric field at ith dipole, P i ~α i E i ; (1) E i E 0 i X j i A ij P j ; (2) which is the sum of the directly incident field E 0 i as well as the field scattered by all the other dipoles j i and it is incident on the ith dipole [second term of Eq. (2)]. The interaction matrix A ij is given from A ij expikr ij r ij k 2 ˆr ij ˆr ij 1 3 ikr ij 1 r 2 ij 3ˆr ij ˆr ij 1 3 ; i j; where 1 3 is the 3 3 unit matrix, r ij r i r j, ˆr ij r ij jr ij j. By combining Eqs. (1) (3), we obtain a linear system of equations, i.e., X N j1 (3) A ij P j E 0 i ; (4) where the diagonal elements of the interaction matrix are essentially the inverse of the polarizability tensor of each dipole, i.e., A ii ~α i 1 : (5) For an anisotropic sphere characterized by a dielectric tensor ~ϵ s and is immersed within an isotropic host of dielectric constant ϵ h, the polarizability tensor of the sphere is given by the Clausius Mossoti formula for anisotropic spheres; that is, Fig. 1. Modeling steps in the hybrid DDA/LMS method. ~α i V s 3ϵ h 4π ~ϵ s ϵ h 1 3 ~ϵ s 2ϵ h : (6) Equation (4) is preferentially solved by conjugate-gradienttype solvers for fast convergence [30]. Having determined the dipole moment P i at each point dipole, one can calculate
3 V. Yannopapas Vol. 31, No. 3 / March 2014 / J. Opt. Soc. Am. B 633 quantities such as the scattering, extinction, and absorption cross sections [10,20 22]. The scattered field E r p at a given point r p in space is the sum of the (secondary) field emitted by each dipole into which the actual scatterer is discretized, i.e., E r p X i A pi P i ; (7) where A pi is provided by Eq. (3) and the point in space does not coincide with one of the dipoles (r p r i ). The polarization vectors P i at each point dipole are provided by Eq. (4). C. Multipole Expansion of the EM Field Let us consider a harmonic EM wave, of angular frequency ω, which is described by its electric-field component: ~Er;tReEr exp iωt: (8) In a homogeneous medium characterized by a dielectric function ϵωϵ 0 and a magnetic permeability μωμ 0, where ϵ 0, μ 0 are the electric permittivity and magnetic permeability of vacuum, Maxwell equations imply that Er satisfies a vector Helmholtz equation, subject to the condition E 0, with p a wave number q ω c, where c 1 p μϵμ 0 ϵ 0 c 0 μϵ is the velocity of light in the medium. The spherical-wave expansion of Er is given by [31] Er X X l l1 m l a H lm f lqrx lm ˆra E i lm q f lqrx lm ˆr ; where a P lm (P E, H) are coefficients to be determined. X lm ˆr are the so-called vector spherical harmonics [31] and f l may be any linear combination of the spherical Bessel function, j l, and the spherical Hankel function, h l. The corresponding magnetic induction, Br, can be readily obtained from Er;t using Maxwell s equations, p ϵμ Br c 0 X X l l1 m l (9) a E lm f lqrx lm ˆr a H i lm q f lqrx lm ˆr ; and we shall now write it down explicitly in what follows. (10) D. Scattering by a Single Scatterer In this subsection, we present a brief summary of the solution to the problem of EM scattering from a single scatterer. We will make use of the compact notation of Ref. [18] for the eigenfunctions and the angular-momentum indices, which allow for easier computer coding. We consider a sphere of radius S, with its center at the origin of coordinates, and assume that its electric permittivity ϵ s and/or magnetic permeability μ s are different from those, ϵ h. μ h of the surrounding homogeneous medium. An EM plane wave incident on this scatterer is described, respectively, by Eq. (9) with f l j l (since the plane wave is finite everywhere), and appropriate coefficients a 0 L, where L denotes collectively the indices Plm. That is, where E 0 r X a 0 L J Lr; (11) L J Elm r i j q l q h rx lm ˆr; J Hlm r j l q h rx lm ˆr; h (12) p and q h ϵ h μ h ω c 0. The coefficients a 0 L depend on the amplitude, polarization and propagation direction of the incident EM plane wave [11 13,31]. Similarly, the wave that is scattered from the sphere is described by Eq. (9) with f l h l, which has the asymptotic form appropriate to an outgoing spherical wave: h l i l expiq h r iq h r as r, and appropriate expansion coefficients a L. Namely, where E r X a L H Lr; (13) L H Elm r i h l q q h rx lm ˆr; h H Hlm r h l q h rx lm ˆr: (14) By applying the requirement that the tangential components of E and H be continuous at the surface of the scatterer, we obtain a relation between the expansion coefficients of the incident and the scattered field, as follows: a L X T LL 0a 0 L0; (15) L 0 where T LL 0 are the elements of the so-called scattering transition T-matrix [32,33]. Eq. (15) is valid for any shape of scatterer; for spherically symmetric scatterers each spherical wave scatters independently of all others, which leads to a transition T-matrix that does not depend on m and is diagonal in l, i.e., T LL 0 T L δ LL 0 [33]. E. Calculation of the T-matrix via Point-Matching The point-matching (PM) method for calculating the T-matrix of general scatterers was first proposed in Ref. [34]. Here, we restate the DDA-based PM method in the spirit of the LMS formalism [18]. Note in passing that a PM method for the calculation of the T-matrix has also been presented in conjunction with the FDTD method [35]. Moreover, the PM method has been employed for obtaining the T-matrix and the multipole moments for a collection of spherical scatterers [36]. We assume that a single spherical EM wave is incident on an arbitrary scatterer, i.e., E 0 r J L0 r. This means that, in Eq. (15), we set a 0 L δ LL 0 which becomes a L T LL 0 : (16) Eq. (13) then becomes X T LL0 H L r p E r p : (17) L
4 634 J. Opt. Soc. Am. B / Vol. 31, No. 3 / March 2014 V. Yannopapas The matrix elements of the L 0 -column of the T-matrix can be calculated from the above equation, provided that we know the scattered field E at a sufficient number of points r p in space. Equation (17) must be solved for different incident spherical waves J L0 to obtain all the T-matrix columns. Namely, we calculate the scattered field via DDA [Eq. (7)] at several points on a spherical surface surrounding the scatterer, so that Eq. (17) becomes a linear system of equations for the T-matrix elements. However, in practice, we calculate the scattered field E via DDA at a large number of points, so that the unknowns T LL0 are fewer than the system equations. In this case, Eq. (17) is solved by seeking a least-squares solution. The scattering/extinction and absorption cross sections of light scattered off a single scatterer can be calculated in a spherical wave expansion by use of the obtained T-matrix [see Eq. (15) of Ref. 14]. F. Hybrid DDA/LMS Method Having calculated the T-matrix via Eq. (17), we insert it within the existing LMS formalism and computer code (see Fig. 2) which allows us to perform light-transmission or band structure calculations for 2D and 3D periodic arrays (photonic and plasmonic crystals, metamaterials, etc.) of scatterers described by the given T-matrix. In case where the 2D unit cell consists of more than one scatterer, e.g., a cluster of scatterers, one can first calculate the total scattering T-matrix of all the scatterers within the unit cell via real-space multiplescattering and then insert the resulting T-matrix within the LMS code, in the manner of Ref. [18]. The latter possibility renders the hybrid DDA/LMS method a multiscale EM solver, as it allows for the study of hierarchically organized photonic structures, where more than one-length scales are involved, e.g., the period of the photonic structure and the characteristic inter-particle distance within the cluster of scatterers. 3. TEST EXAMPLE We consider a 2D square array with lattice constant a 150 nm consisting of 100 nm gold nanocubes. The dielectric function of gold is taken from experiment [37]. Figure 3 shows the cross sections for linearly polarized light incident normally on a nanocube face. The curves are obtained by discretizing the cubes into dipoles while the cutoff in the angular momentum expansion is taken l max 7. We note that the cross sections already converge for a lower cutoff (l max 4); however, for the periodic photonic structures of Fig. 2. Flowchart of the hybrid DDA/LMS method. Fig. 3. Extinction (solid), scattering (dotted), and absorption (dashed line) cross sections (in arbitrary units) for light incident normally on one of the faces of gold nanocube with 100 nm edge. The curves are calculated by assuming 512 dipoles in the DDA and l max 7 in the angular momentum expansion. The inset shows the relative error in the scattering spectrum when assuming angular momentum cutoffs l max 6 and l max 7 (black line), as well as the DDA convergence, when assuming and point dipoles within the nanocube (gray line). gold nanocubes studied right below, a higher angular momentum cutoff (l max 7) is needed to achieve convergence [11 13]. The computer time needed for the calculation of the cross sections depicted in Fig. 3 was 9.25 h for 101 frequency points on a single core i7-3930k 3.2 GHz while it required about 620 MB of RAM. We note here that a computer time per frequency is meaningless in our case since the iterations needed for solving the DDA Eq. (4) via a conjugategradient-type solver [30] depends critically on frequency [26]. The evident peak in all three curves is attributed to the excitation of the surface plasmon resonance at 540 nm (when defined from the absorption cross section). In the inset of Fig. 3 we show the convergence (spectrum of relative error) of the scattering cross section in terms of the angular momentum expansion (black line) and of the number of dipoles (gray line) into which the nanocube is discretized in the DDA. In Fig. 4 we show the spectra of light transmittance/ reflectance/absorbance from a 2D square array of 100 nm gold nanocubes (those of Fig. 3) with lattice constant a 150 nm. In both Figs. 3 and 4, convergence is achieved for l max 7 in the EM spherical wave expansion and for 37 reciprocal-lattice vectors in the EM plane wave expansion [11 13] for Fig. 4. For these parameters, the relative error in the transmittance averaged over the studied spectrum ( nm) is 1.3%, in the reflectance curve 0.3%, and in the absorbance curve 0.8%. The computer time needed for the calculation of the transmittance/reflectance/absorbance spectra depicted in Fig. 4 was about 1.3 s per frequency on a single core i7-3930k 3.2 GHz while it required around 400 MB of RAM. As the only quantity that can be compared with the case of a single nanocube (Fig. 3) is absorbance, we observe a strong redshift of the surface plasmon peak, now located at 576 nm, which is a result of the interaction of the surface plasmons among the nanocubes of the square lattice [38]. Finally, Fig. 5 studies the case of a 3D plasmonic lattice; namely, a simple cubic lattice of 100 nm gold nanocubes with
5 V. Yannopapas Vol. 31, No. 3 / March 2014 / J. Opt. Soc. Am. B 635 Fig. 4. Transmittance, reflectance, and absorbance for light incident normally on an infinitely periodic 2D square array of 100 nm gold nanocubes with lattice constant a 150 nm. lattice constant of 150 nm, viewed as a succession of layers parallel to the (001) crystallographic plane. Figure 5(a) shows the dispersion curves (band structure) for an infinite crystal and Figs. 5(b) and 5(c) show the reflectance and absorbance spectra from finite slabs of the crystal, respectively. Figure 5(a) depicts both the real Rk z and imaginary Ik z parts of the wavevector k z, i.e., normal to the (001) crystallographic plane; from Rk z and Ik z, the effective-medium properties of the structure can be obtained. In Figs. 5(b) and 5(c), we only depict plasmonic-crystal slabs up to 4 layers as both the reflectance and absorbance spectra saturate very fast with increasing slab thickness. This means that adding more than 4 layers of the crystal does not practically change the depicted spectra. From applications point of view, even a single layer of a lattice of nanocubes may serve as an efficient broadband light absorber thanks to the wide absorbance plateau from 300 to about 500 nm. In the case studies presented here, the nanocubes are separated from each other by 50 nm. As the inter-cube distance, Fig. 5. (a) Dispersion curves along the [001] crystallographic direction of simple cubic crystal of 100 nm gold nanocubes with lattice constant a 150 nm. The solid (broken) line depicts the dispersion of the real (imaginary) part of k z. (b) Reflectance and (c) absorbance of light incident normally on finite slabs of various thicknesses (number of layers) of the simple cubic crystal of (a). i.e., the lattice constant of the square lattice, becomes smaller and the nanocubes approach each other, the number of terms needed in the spherical-wave expansion of the EM field increases rapidly. This is caused by the unnatural choice of a local dielectric function for metallic scatterers when the latter come very close to each other. When a nonlocal dielectric function is taken into account this effect is mitigated and an angular-momentum cutoff l max 15 is sufficient for converged results in the cases under study. However, the DDA/ LMS programming code scales linearly with l max as the time-consuming part is the solution of DDA system, i.e., Eq. (4). As stated above, Eq. (4) must be solved for different spherical waves whose number scales as (2l max 1). So, small inter-particle distances within a single plane of particles does not prohibitively increase the computer time as the corresponding programming code of the hybrid DDA/LMS method scales linearly with l max. A last point on the convergence of the presented method. Since the DDA method is ideally suited for scattering targets of small to medium size compared with the operating wavelength, so does the hybrid DDA/LMS method. For scatterers of large size, time-domain methods (FDTD and PSTD) are superior [39,40] in terms of computer time and memory suggesting that a hybrid FDTD/LMS [41] or PSTD/LMS code is more favorable in this case. 4. CONCLUSION We have presented a novel electromagnetic solver, namely, a hybrid discrete-dipole approximation/layermultiple-scattering method that can treat photonic structures consisting scatterers, which, in general, may be nonspherical, inhomogeneous, and/or anisotropic. The present method is a major extension of the existing layer-multiple-scattering formalism and code that currently treats isotropic [12,13], anisotropic [16,17] spheres, and isotropic axisymmetric [14] as scatterers. The developed method is superior to existing freeware and commercial EM packages (FEM, FDTD) since, due to its modular formalism and structure, it is very fast and accurate in treating 3D metamaterials and photonic/ plasmonic crystals of arbitrary scatterers while it provides the complex band structure of an infinite crystal. At the same time, the presented method possesses multiscale functionalities as it solves Maxwell s equations in different multiplescattering stages (the DDA method is generally considered as such) allowing for the treatment of complex photonic structures, such as 3D hierarchically organized metamaterials characterized by more than one length scale. REFERENCES 1. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals (Princeton, 1995). 2. F. G. J. de Abajo, Colloquium: light scattering by particle and hole arrays, Rev. Mod. Phys. 79, (2007). 3. L. Solymar and E. Shamonina, Waves in Metamaterials (Oxford, 2009). 4. P. M. Bell, J. B. Pendry, L. Martín-Moreno, and A. J. Ward, A program for calculating photonic band structures and transmission coefficients of complex structures, Comput. Phys. Commun. 85, (1995). 5. P. Monk, Finite Element Methods for Maxwell s Equations (Clarendon, 2003). 6. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method (Artech House, 2005).
6 636 J. Opt. Soc. Am. B / Vol. 31, No. 3 / March 2014 V. Yannopapas 7. Q. H. Liu, The PSTD algorithm: a time-domain method requiring only two cells per wavelength, Microw. Opt. Technol. Lett. 15, (1997). 8. C. Liu, R. L. Panetta, and P. Yang, Application of the pseudospectral time-domain method to compute particle singlescattering properties for size parameters up to 200, J. Quant. Spectrosc. Radiat. Transfer 113, (2012). 9. S. G. Tikhodeev, A. L. Yablonksii, E. A. Muljarov, N. A. Gippius, and T. Ishihara, Quasiguided modes and optical properties of photonic crystal slabs, Phys. Rev. 66, (2002). 10. E. M. Purcell and C. R. Pennypacker, Scattering and absorption of light by nonspherical dielectric grains, Astrophys. J. 186, (1973). 11. N. Stefanou, V. Karathanos, and A. Modinos, Scattering of electromagnetic waves by periodic structures, J. Phys. Condens. Mat. 4, 7389 (1992). 12. N. Stefanou, V. Yannopapas, and A. Modinos, Heterostructures of photonic crystals: frequency bands and transmission coefficients, Comput. Phys. Commun. 113, 49 (1998). 13. N. Stefanou, V. Yannopapas, and A. Modinos, MULTEM 2: a new version of the program for transmission and band-structure calculations of photonic crystals, Comput. Phys. Commun. 132, 189 (2000). 14. G. Gantzounis and N. Stefanou, Layer-multiple-scattering method for photonic crystals of nonspherical particles, Phys. Rev. B 73, (2006). 15. G. Gantzounis, N. Stefanou, and N. Paanikolaou, Optical properties of periodic structures of metallic nanodisks, Phys. Rev. B 77, (2008). 16. A. Christofi and N. Stefanou, Nonreciprocal optical response of helical periodic structures of plasma spheres in a static magnetic field, Phys. Rev. B 87, (2013). 17. A. Christofi and N. Stefanou, Nonreciprocal photonic surface states in periodic structures of magnetized plasma nanospheres, Phys. Rev. B 88, (2013). 18. V. Yannopapas and A. G. Vanakaras, Layer-multiple scattering theory for metamaterials made from clusters of nanoparticles, Phys. Rev. B 84, (2011). 19. V. Myroshnychenko, J. Rodríguez-Fernádez, I. Pastoriza-Santos, A. M. Funston, C. Novo, P. Mulvaney, L. M. Liz-Marzán, and F. J. G. de Abajo, Modeling the optical response of gold nanoparticles, Chem. Soc. Rev. 37, (2008). 20. B. T. Draine, The discrete-dipole approximation and its application to stellar graphite grains, Astrophys. J. 333, (1988). 21. B. T. Draine and P. J. Flatau, Discrete-dipole approximation for scattering calculations, J. Opt. Soc. Am. A 11, (1994). 22. M. A. Yurkin and A. G. Hoekstra, The discrete dipole approximation: an overview and recent developments, J. Quant. Spectrosc. Radiat. Transfer 106, (2007). 23. B. T. Draine and P. J. Flatau, Discrete-dipole approximation for periodic targets: theory and tests, J. Opt. Soc. Am. A 25, (2008). 24. A. B. Evlyukhin, C. Reinhardt, and B. N. Chichkov, Multipole light scattering by nonspherical nanoparticles in the discrete dipole approximation, Phys. Rev. B 84, (2011). 25. A. B. Evlyukhin, C. Reinhardt, E. Evyukhin, and B. N. Chichkov, Multipole analysis of light scattering by arbitrary-shaped nanoparticles on a plane surface, J. Opt. Soc. Am. B 30, (2013). 26. M. A. Yurkin, D. de Kanter, and A. G. Hoekstra, Accuracy of the discrete dipole approximation for simulation of optical properties of gold nanoparticles, J. Nanophoton. 4, (2010). 27. A. B. Evlyukhin, C. Reinhardt, A. Seidel, B. S. Luk yanchuk, and B. N. Chichkov, Optical response features of Si-nanoparticle arrays, Phys. Rev. B 82, (2010). 28. A. B. Evlyukhin, C. Reinhardt, U. Zywietz, and B. N. Chichkov, Collective resonances in metal nanoparticle arrays with dipole-quadrupole interactions, Phys. Rev. B 85, (2012). 29. S. Droulias and V. Yannopapas, Broad-band giant circular dichroism in metamaterials of twisted chains of metallic nanoparticles, J. Phys. Chem. C 117, (2013). 30. P. C. Chaumet and A. Rahmani, Efficient discrete dipole approximation for magnetoelectric scatterers, Opt. Lett. 34, (2009). 31. J. D. Jackson, Classical Electrodynamics (Wiley, 1999). 32. P. C. Waterman, Symmetry, unitarity, and geometry in electromagnetic scattering, Phys. Rev. D 3, (1971). 33. C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983). 34. V. L. Y. Loke, T. A. Nieminen, N. R. Heckenberg, and H. Rubinsztein-Dunlop, T-matrix calculation via discrete dipole approximation, point matching and exploiting symmetry, J. Quant. Spectrosc. Radiat. Transfer 110, (2009). 35. V. L. Y. Loke, T. A. Nieminen, S. J. Parkin, N. R. Heckenberg, and H. Rubinsztein-Dunlop, FDTD/T-matrix hybrid method, J. Quant. Spectrosc. Radiat. Transfer 106, (2007). 36. C. Menzel, S. Mühlig, C. Rockstuhl, and F. Lederer, Multipole analysis of meta-atoms, Metamaterials 5, (2011). 37. R. B. Johnson and R. W. Christy, Optical constants of the noble metals, Phys. Rev. B 6, (1972). 38. V. Yannopapas, A. Modinos, and N. Stefanou, Optical properties of metallodielectric photonic crystals, Phys. Rev. B 60, (1999). 39. M. A. Yurkin, A. G. Hoekstra, R. S. Brock, and J. Q. Lu, Systematic comparison of the discrete dipole approximation and the finite difference time-domain method for large dielectric scatterers, Opt. Express 15, (2007). 40. C. Liu, L. Bi, R. L. Panetta, P. Yang, and M. A. Yurkin, Comparison between the pseudo-spectral time-domain method and the discrete dipole approximation for light scattering simulations, Opt. Express 20, (2012). 41. A. Deinega, S. Belousov, and I. Valuev, Transfer-matrix approach for finite-difference time-domain simulation of periodic structures, Phys. Rev. E 88, (2013).
Scattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
More informationarxiv:astro-ph/ v1 13 Nov 2003
Discrete-dipole approximation with polarizabilities that account for both finite wavelength and target geometry Matthew J. Collinge and B. T. Draine Princeton University Observatory, Princeton, New Jersey
More informationPhotonic crystals: a novel class of functional materials
Materials Science-Poland, Vol. 23, No. 4, 2005 Photonic crystals: a novel class of functional materials A. MODINOS 1, N. STEFANOU 2* 1 Department of Physics, National Technical University of Athens, Zografou
More informationBackscattering enhancement of light by nanoparticles positioned in localized optical intensity peaks
Backscattering enhancement of light by nanoparticles positioned in localized optical intensity peaks Zhigang Chen, Xu Li, Allen Taflove, and Vadim Backman We report what we believe to be a novel backscattering
More informationSupplementary Figure 1. Schematics of light transmission and reflection from a slab confined between
Supplementary Figures: Supplementary Figure. Schematics of light transmission and reflection from a slab confined between two infinite media. Supplementary Figure. Reflectivity of a magneto-electric slab
More informationPhotonic/Plasmonic Structures from Metallic Nanoparticles in a Glass Matrix
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Photonic/Plasmonic Structures from Metallic Nanoparticles in a Glass Matrix O.Kiriyenko,1, W.Hergert 1, S.Wackerow 1, M.Beleites 1 and
More informationSupplementary Figure S1 SEM and optical images of Si 0.6 H 0.4 colloids. a, SEM image of Si 0.6 H 0.4 colloids. b, The size distribution of Si 0.
Supplementary Figure S1 SEM and optical images of Si 0.6 H 0.4 colloids. a, SEM image of Si 0.6 H 0.4 colloids. b, The size distribution of Si 0.6 H 0.4 colloids. The standard derivation is 4.4 %. Supplementary
More informationSCATTERING OF ELECTROMAGNETIC WAVES ON METAL NANOPARTICLES. Tomáš Váry, Juraj Chlpík, Peter Markoš
SCATTERING OF ELECTROMAGNETIC WAVES ON METAL NANOPARTICLES Tomáš Váry, Juraj Chlpík, Peter Markoš ÚJFI, FEI STU, Bratislava E-mail: tomas.vary@stuba.sk Received xx April 2012; accepted xx May 2012. 1.
More informationLight transmission through a single cylindrical hole in a metallic film
Light transmission through a single cylindrical hole in a metallic film F. J. García de Abajo Centro Mixto CSIC-UPV/EHU and Donostia International Physics Center (DIPC), Aptdo. 1072, 20080 San Sebastián,
More informationLight Localization in Left-Handed Media
Vol. 112 (2007) ACTA PHYSICA POLONICA A No. 4 Proceedings of the 3rd Workshop on Quantum Chaos and Localisation Phenomena Warsaw, Poland, May 25 27, 2007 Light Localization in Left-Handed Media M. Rusek,
More informationPhotonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique
Photonic nanojet enhancement of backscattering of light by nanoparticles: a potential novel visible-light ultramicroscopy technique Zhigang Chen and Allen Taflove Department of Electrical and Computer
More informationME equations. Cylindrical symmetry. Bessel functions 1 kind Bessel functions 2 kind Modifies Bessel functions 1 kind Modifies Bessel functions 2 kind
Δϕ=0 ME equations ( 2 ) Δ + k E = 0 Quasi static approximation Dynamic approximation Cylindrical symmetry Metallic nano wires Nano holes in metals Bessel functions 1 kind Bessel functions 2 kind Modifies
More informationAngular and polarization properties of a photonic crystal slab mirror
Angular and polarization properties of a photonic crystal slab mirror Virginie Lousse 1,2, Wonjoo Suh 1, Onur Kilic 1, Sora Kim 1, Olav Solgaard 1, and Shanhui Fan 1 1 Department of Electrical Engineering,
More informationTime Domain Modeling of All-Optical Switch based on PT-Symmetric Bragg Grating
Time Domain Modeling of All-Optical Switch based on PT-Symmetric Bragg Grating Sendy Phang 1, Ana Vukovic 1, Hadi Susanto 2, Trevor M. Benson 1, and Phillip Sewell 1 1 School of Electrical and Electronic
More informationON THE IMPORTANCE OF LOCAL-FIELD CORRECTIONS FOR POLARIZABLE PARTICLES ON A FINITE LATTICE: APPLICATION TO THE DISCRETE DIPOLE APPROXIMATION
The Astrophysical Journal, 607:873 878, 2004 June 1 # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. ON THE IMPORTANCE OF LOCAL-FIELD CORRECTIONS FOR POLARIZABLE PARTICLES
More informationBehavior of light at photonic crystal interfaces
Behavior of light at photonic crystal interfaces Emanuel Istrate, Alexander A. Green, and Edward H. Sargent Department of Electrical and Computer Engineering, University of Toronto, 10 King s College Road,
More informationEnhancing and suppressing radiation with some permeability-near-zero structures
Enhancing and suppressing radiation with some permeability-near-zero structures Yi Jin 1,2 and Sailing He 1,2,3,* 1 Centre for Optical and Electromagnetic Research, State Key Laboratory of Modern Optical
More informationMie theory for light scattering by a spherical particle in an absorbing medium
Mie theory for light scattering by a spherical particle in an absorbing medium Qiang Fu and Wenbo Sun Analytic equations are developed for the single-scattering properties of a spherical particle embedded
More informationComparison between the pseudo-spectral time domain method and the discrete dipole approximation for light scattering simulations
Comparison between the pseudo-spectral time domain method and the discrete dipole approximation for light scattering simulations Chao Liu, 1 Lei Bi, 1 R. Lee Panetta, 1 Ping Yang, 1,* and Maxim A. Yurkin
More informationAbsorption suppression in photonic crystals
PHYSICAL REVIEW B 77, 442 28 Absorption suppression in photonic crystals A. Figotin and I. Vitebskiy Department of Mathematics, University of California at Irvine, Irvine, California 92697, USA Received
More informationScattering of light from quasi-homogeneous sources by quasi-homogeneous media
Visser et al. Vol. 23, No. 7/July 2006/J. Opt. Soc. Am. A 1631 Scattering of light from quasi-homogeneous sources by quasi-homogeneous media Taco D. Visser* Department of Physics and Astronomy, University
More informationDesign of a Multi-Mode Interference Crossing Structure for Three Periodic Dielectric Waveguides
Progress In Electromagnetics Research Letters, Vol. 75, 47 52, 2018 Design of a Multi-Mode Interference Crossing Structure for Three Periodic Dielectric Waveguides Haibin Chen 1, Zhongjiao He 2,andWeiWang
More informationScattering of ECRF waves by edge density fluctuations and blobs
PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA
More informationDefect-based Photonic Crystal Cavity for Silicon Laser
Defect-based Photonic Crystal Cavity for Silicon Laser Final Term Paper for Nonlinear Optics PHYC/ECE 568 Arezou Khoshakhlagh Instructor: Prof. M. Sheikh-Bahae University of New Mexico karezou@unm.edu
More informationPlasmonic metamaterial cloaking at optical frequencies
Plasmonic metamaterial cloaking at optical frequencies F. Bilotti *, S. Tricarico, and L. Vegni Department of Applied Electronics, University Roma Tre Via della Vasca Navale 84, Rome 146, ITALY * Corresponding
More informationSupplementary information for. plasmonic nanorods interacting with J-aggregates.
Supplementary information for Approaching the strong coupling limit in single plasmonic nanorods interacting with J-aggregates. by Gülis Zengin, Göran Johansson, Peter Johansson, Tomasz J. Antosiewicz,
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationOn the signs of the imaginary parts of the effective permittivity and permeability in metamaterials
1016 J. Opt. Soc. Am. B/ Vol. 27, No. 5/ May 2010 J. Woodley and M. Mojahedi On the signs of the imaginary parts of the effective permittivity and permeability in metamaterials J. Woodley 1, * and M. Mojahedi
More informationLocalized surface plasmons (Particle plasmons)
Localized surface plasmons (Particle plasmons) ( Plasmons in metal nanostructures, Dissertation, University of Munich by Carsten Sonnichsen, 2001) Lycurgus cup, 4th century (now at the British Museum,
More informationSURFACE PLASMONS AND THEIR APPLICATIONS IN ELECTRO-OPTICAL DEVICES
SURFACE PLASMONS AND THEIR APPLICATIONS IN ELECTRO-OPTICAL DEVICES Igor Zozouleno Solid State Electronics Department of Science and Technology Linöping University Sweden igozo@itn.liu.se http://www.itn.liu.se/meso-phot
More informationElectromagnetic Metamaterials
Photonic Bandgap and Electromagnetic Metamaterials Andrew Kirk andrew.kirk@mcgill.ca ca Department of Electrical and Computer Engineering McGill Institute for Advanced Materials A Kirk 11/24/2008 Photonic
More informationNegative Index of Refraction in Optical Metamaterials
1 Negative Index of Refraction in Optical Metamaterials V. M. Shalaev, W. Cai, U. Chettiar, H.-K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev School of Electrical and Computer Engineering,
More informationOverview. 1. What range of ε eff, µ eff parameter space is accessible to simple metamaterial geometries? ``
MURI-Transformational Electromagnetics Innovative use of Metamaterials in Confining, Controlling, and Radiating Intense Microwave Pulses University of New Mexico August 21, 2012 Engineering Dispersive
More informationLarge omnidirectional band gaps in metallodielectric photonic crystals
PHYSICAL REVIEW B VOLUME, NUMBER 16 15 OCTOBER 1996-II Large omnidirectional band gaps in metallodielectric photonic crystals Shanhui Fan, Pierre R. Villeneuve, and J. D. Joannopoulos Department of Physics,
More informationInside-out electromagnetic cloaking
Inside-out electromagnetic cloaking Nina A. Zharova 1,2, Ilya V. Shadrivov 1, and Yuri S. Kivshar 1 1 Nonlinear Physics Center, Research School of Physical Sciences and Engineering, Australian National
More informationNanoscale shift of the intensity distribution of dipole radiation
Shu et al. Vol. 26, No. 2/ February 2009/ J. Opt. Soc. Am. A 395 Nanoscale shift of the intensity distribution of dipole radiation Jie Shu, Xin Li, and Henk F. Arnoldus* Department of Physics and Astronomy,
More informationExtinction properties of a sphere with negative permittivity and permeability
PERGAMON Solid State Communications 116 (2000) 411 415 www.elsevier.com/locate/ssc Extinction properties of a sphere with negative permittivity and permeability R. Ruppin* Department of Physics and Applied
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationSimulated Study of Plasmonic Coupling in Noble Bimetallic Alloy Nanosphere Arrays
CHAPTER 4 Simulated Study of Plasmonic Coupling in Noble Bimetallic Alloy Nanosphere Arrays 4.1 Introduction In Chapter 3, the noble bimetallic alloy nanosphere (BANS) of Ag 1-x Cu x at a particular composition
More informationCircularly polarized thermal emission from chiral metasurface in the absence of magnetic field
Journal of Physics: Conference Series PAPER OPEN ACCESS Circularly polarized thermal emission from chiral metasurface in the absence of magnetic field To cite this article: S.A. Dyakov et al 2018 J. Phys.:
More informationFrequency Dependence Effective Refractive Index of Meta Materials by Effective Medium Theory
Advance in Electronic and Electric Engineering. ISSN 31-197, Volume 3, Number (13), pp. 179-184 Research India Publications http://www.ripublication.com/aeee.htm Frequency Dependence Effective Refractive
More informationElectromagnetic Interactions within Small Particles The Discrete Dipole Approximation Tristan Sharp
Electromagnetic Interactions within Small Particles The Discrete Dipole Approximation Tristan Sharp Small particles can interact with light in non-intuitive ways, scattering much less light than expected
More informationRouting of Deep-Subwavelength Optical Beams and Images without Reflection and Diffraction Using Infinitely Anisotropic Metamaterials
Peter B. Catrysse * and Shanhui Fan Routing of Deep-Subwavelength Optical Beams and Images without Reflection and Diffraction Using Infinitely Anisotropic Metamaterials Media that are described by extreme
More informationSuper-reflection and Cloaking Based on Zero Index Metamaterial
Super-reflection and Cloaking Based on Zero Index Metamaterial Jiaming Hao, Wei Yan, and Min Qiu Photonics and Microwave ngineering, Royal Institute of Technology (KTH), lectrum 9, 164 4, Kista, Sweden
More informationThe Broadband Fixed-Angle Source Technique (BFAST) LUMERICAL SOLUTIONS INC
The Broadband Fixed-Angle Source Technique (BFAST) LUMERICAL SOLUTIONS INC. 1 Outline Introduction Lumerical s simulation products Simulation of periodic structures The new Broadband Fixed-Angle Source
More informationlimitations J. Zhou, E. N. Economou and C. M. Soukoulis
Mesoscopic Physics in Complex Media, 01011 (010) DOI:10.1051/iesc/010mpcm01011 Owned by the authors, published by EDP Sciences, 010 Optical metamaterials: Possibilities and limitations M. Kafesaki, R.
More informationLight scattering by irregular interplanetary dust particles
Earth Planets Space, 50, 577 585, 1998 Light scattering by irregular interplanetary dust particles Hajime Okamoto 1 and Yu-lin Xu 2 1 Kashima Space Research Center, Communications Research Laboratory,
More informationAuthor(s) Tamayama, Y; Nakanishi, T; Sugiyama. Citation PHYSICAL REVIEW B (2006), 73(19)
Observation of Brewster's effect fo Titleelectromagnetic waves in metamateri theory Author(s) Tamayama, Y; Nakanishi, T; Sugiyama Citation PHYSICAL REVIEW B (2006), 73(19) Issue Date 2006-05 URL http://hdl.handle.net/2433/39884
More informationNanomaterials and their Optical Applications
Nanomaterials and their Optical Applications Winter Semester 2013 Lecture 02 rachel.grange@uni-jena.de http://www.iap.uni-jena.de/multiphoton Lecture 2: outline 2 Introduction to Nanophotonics Theoretical
More informationSupplementary Material for. Resonant Transparency and Non-Trivial Excitations in Toroidal Metamaterials
Supplementary Material for Resonant Transparency and Non-Trivial Excitations in Toroidal Metamaterials V. A. Fedotov 1, A. V. Rogacheva 1, V. Savinov 1, D. P. Tsai 2,3, N. I. Zheludev 1, 4 1 Optoelectronics
More informationarxiv: v2 [physics.optics] 5 Feb 2009
All-Dielectric Rod-Type Metamaterials at Optical Frequencies K. Vynck, D. Felbacq, E. Centeno, A. I. Căbuz, D. Cassagne, and B. Guizal Groupe d Etude des Semiconducteurs, UMR 5650 CNRS-UM2, CC074, Place
More informationEnergy transport in metal nanoparticle plasmon waveguides
Energy transport in metal nanoparticle plasmon waveguides Stefan A. Maier, Pieter G. Kik, and Harry A. Atwater California Institute of Technology Thomas J. Watson Laboratory of Applied Physics, Pasadena,
More informationSupplementary Figure 1 SEM images and corresponding Fourier Transformation of nanoparticle arrays before pattern transfer (left), after pattern
Supplementary Figure 1 SEM images and corresponding Fourier Transformation of nanoparticle arrays before pattern transfer (left), after pattern transfer but before pattern shrinkage (middle), and after
More informationtransmission reflection absorption
Optical Cages V. Kumar*, J. P. Walker* and H. Grebel The Electronic Imaging Center and the ECE department at NJIT, Newark, NJ 0702. grebel@njit.edu * Contributed equally Faraday Cage [], a hollow structure
More informationSub-wavelength electromagnetic structures
Sub-wavelength electromagnetic structures Shanhui Fan, Z. Ruan, L. Verselegers, P. Catrysse, Z. Yu, J. Shin, J. T. Shen, G. Veronis Ginzton Laboratory, Stanford University http://www.stanford.edu/group/fan
More informationDemonstration of Near-Infrared Negative-Index Materials
Demonstration of Near-Infrared Negative-Index Materials Shuang Zhang 1, Wenjun Fan 1, N. C. Panoiu 2, K. J. Malloy 1, R. M. Osgood 2 and S. R. J. Brueck 2 1. Center for High Technology Materials and Department
More informationPlasmonics: elementary excitation of a plasma (gas of free charges) nano-scale optics done with plasmons at metal interfaces
Plasmonics Plasmon: Plasmonics: elementary excitation of a plasma (gas of free charges) nano-scale optics done with plasmons at metal interfaces Femius Koenderink Center for Nanophotonics AMOLF, Amsterdam
More information7. Localized surface plasmons (Particle plasmons)
7. Localized surface plasmons (Particle plasmons) ( Plasmons in metal nanostructures, Dissertation, University of Munich by Carsten Sonnichsen, 2001) Lycurgus cup, 4th century (now at the British Museum,
More informationChapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice
Chapter 5 Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice In chapter 3 and 4, we have demonstrated that the deformed rods, rotational rods and perturbation
More informationNanoscale optical circuits: controlling light using localized surface plasmon resonances
Nanoscale optical circuits: controlling light using localized surface plasmon resonances T. J. Davis, D. E. Gómez and K. C. Vernon CSIRO Materials Science and Engineering Localized surface plasmon (LSP)
More informationNew Aspects of Old Equations: Metamaterials and Beyond (Part 2) 신종화 KAIST 물리학과
New Aspects of Old Equations: Metamaterials and Beyond (Part 2) 신종화 KAIST 물리학과 Metamaterial Near field Configuration in Periodic Structures New Material Material and Metamaterial Material Metamaterial
More informationSpontaneous emission rate of an electric dipole in a general microcavity
PHYSICAL REVIEW B VOLUME 60, NUMBER 7 15 AUGUST 1999-I Spontaneous emission rate of an electric dipole in a general microcavity Jeong-Ki Hwang, Han-Youl Ryu, and Yong-Hee Lee Department of Physics, Korea
More informationNegative Refraction and Subwavelength Lensing in a Polaritonic Crystal
Negative Refraction and Subwavelength Lensing in a Polaritonic Crystal X. Wang and K. Kempa Department of Physics, Boston College Chestnut Hill, MA 02467 We show that a two-dimensional polaritonic crystal,
More informationA Novel Design of Photonic Crystal Lens Based on Negative Refractive Index
PIERS ONLINE, VOL. 4, NO. 2, 2008 296 A Novel Design of Photonic Crystal Lens Based on Negative Refractive Index S. Haxha 1 and F. AbdelMalek 2 1 Photonics Group, Department of Electronics, University
More informationModeling Focused Beam Propagation in scattering media. Janaka Ranasinghesagara, Ph.D.
Modeling Focused Beam Propagation in scattering media Janaka Ranasinghesagara, Ph.D. Teaching Objectives The need for computational models of focused beam propagation in scattering media Introduction to
More informationResearch on the Wide-angle and Broadband 2D Photonic Crystal Polarization Splitter
Progress In Electromagnetics Research Symposium 2005, Hangzhou, China, August 22-26 551 Research on the Wide-angle and Broadband 2D Photonic Crystal Polarization Splitter Y. Y. Li, P. F. Gu, M. Y. Li,
More informationUSAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION
October 4-6, 2007 - Chiinu, Rep.Moldova USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION Ionu- P. NICA, Valeriu Gh. DAVID, /tefan URSACHE Gh. Asachi Technical University Iai, Faculty
More informationAn efficient way to reduce losses of left-handed metamaterials
An efficient way to reduce losses of left-handed metamaterials Jiangfeng Zhou 1,2,, Thomas Koschny 1,3 and Costas M. Soukoulis 1,3 1 Ames Laboratory and Department of Physics and Astronomy,Iowa State University,
More informationDensity of modes maps for design of photonic crystal devices
RESEARCH Revista Mexicana de Física 62 (2016) 193 198 MAY-JUNE 2016 Density of modes maps for design of photonic crystal devices I. Guryev, I. Sukhoivanov, R.I. Mata Chavez, N. Gurieva, J.A. Andrade Lucio,
More informationTaking cascaded plasmonic field enhancement to the ultimate limit in silver nanoparticle dimers S. Toroghi* a, P. G. Kik a,b
Taking cascaded plasmonic field enhancement to the ultimate limit in silver nanoparticle dimers S. Toroghi* a, P. G. Kik a,b a CREOL, The College of Optics and Photonics, University of Central Florida,
More informationSupporting Information
Supporting Information Light emission near a gradient metasurface Leonard C. Kogos and Roberto Paiella Department of Electrical and Computer Engineering and Photonics Center, Boston University, Boston,
More informationSCATTERING CROSS SECTION OF A META-SPHERE
Progress In Electromagnetics Research Letters, Vol. 9, 85 91, 009 SCATTERING CROSS SECTION OF A META-SPHERE A. Alexopoulos Electronic Warfare and Radar Division Defence Science and Technology Organisation
More informationInvestigation of one-dimensional photonic bandgap structures containing lossy double-negative metamaterials through the Bloch impedance
Shi et al. Vol. 3, No. 6 / June 23 / J. Opt. Soc. Am. B 473 Investigation of one-dimensional photonic bandgap structures containing lossy double-negative metamaterials through the Bloch impedance Fenghua
More informationNon-Rayleigh limit of the Lorenz-Mie solution and suppression of scattering by spheres of negative refractive index
PHYSICAL REVIEW A 80, 013808 2009 Non-Rayleigh limit of the Lorenz-Mie solution and suppression of scattering by spheres of negative refractive index Andrey E. Miroshnichenko* Nonlinear Physics Centre
More informationAAPP Atti della Accademia Peloritana dei Pericolanti Classe di Scienze Fisiche, Matematiche e Naturali ISSN
DOI: 1.1478/C1V89S1P33 AAPP Atti della Accademia Peloritana dei Pericolanti Classe di Scienze Fisiche, Matematiche e Naturali ISSN 1825-1242 Vol. 89, Suppl. No. 1, C1V89S1P33 (211) ON RAYLEIGH APPROXIMATION
More informationEXTENSIONS OF THE COMPLEX JACOBI ITERATION TO SIMULATE PHOTONIC WAVELENGTH SCALE COMPONENTS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux Eds c TU Delft, The Netherlands, 2006 EXTENSIONS OF THE COMPLEX JACOBI ITERATION TO SIMULATE PHOTONIC
More informationFundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation
Fundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation Klaus Jockers November 11, 2014 Max-Planck-Institut für Sonnensystemforschung
More informationSpring 2009 EE 710: Nanoscience and Engineering
Spring 009 EE 710: Nanoscience and Engineering Part 10: Surface Plasmons in Metals Images and figures supplied from Hornyak, Dutta, Tibbals, and Rao, Introduction to Nanoscience, CRC Press Boca Raton,
More informationPHYSICAL REVIEW B 71,
Coupling of electromagnetic waves and superlattice vibrations in a piezomagnetic superlattice: Creation of a polariton through the piezomagnetic effect H. Liu, S. N. Zhu, Z. G. Dong, Y. Y. Zhu, Y. F. Chen,
More informationPhotonic devices for quantum information processing:
Outline Photonic devices for quantum information processing: coupling to dots, structure design and fabrication Optoelectronics Group, Cavendish Lab Outline Vuckovic s group Noda s group Outline Outline
More informationAs a partial differential equation, the Helmholtz equation does not lend itself easily to analytical
Aaron Rury Research Prospectus 21.6.2009 Introduction: The Helmhlotz equation, ( 2 +k 2 )u(r)=0 1, serves as the basis for much of optical physics. As a partial differential equation, the Helmholtz equation
More informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 26 Aug 2003
Circular dichroism simulated spectra of chiral gold nanoclusters: A dipole approximation arxiv:cond-mat/0308552v1 [cond-mat.mtrl-sci] 26 Aug 2003 Carlos E. Román-Velázquez, 1 Cecilia Noguez, 2, and Ignacio
More informationCLASSICAL ELECTRICITY
CLASSICAL ELECTRICITY AND MAGNETISM by WOLFGANG K. H. PANOFSKY Stanford University and MELBA PHILLIPS Washington University SECOND EDITION ADDISON-WESLEY PUBLISHING COMPANY Reading, Massachusetts Menlo
More informationSuperconductivity Induced Transparency
Superconductivity Induced Transparency Coskun Kocabas In this paper I will discuss the effect of the superconducting phase transition on the optical properties of the superconductors. Firstly I will give
More informationNanomaterials and their Optical Applications
Nanomaterials and their Optical Applications Winter Semester 2012 Lecture 08 rachel.grange@uni-jena.de http://www.iap.uni-jena.de/multiphoton Outline: Photonic crystals 2 1. Photonic crystals vs electronic
More information1. Reminder: E-Dynamics in homogenous media and at interfaces
0. Introduction 1. Reminder: E-Dynamics in homogenous media and at interfaces 2. Photonic Crystals 2.1 Introduction 2.2 1D Photonic Crystals 2.3 2D and 3D Photonic Crystals 2.4 Numerical Methods 2.5 Fabrication
More informationPhotonic band gap engineering in 2D photonic crystals
PRAMANA c Indian Academy of Sciences Vol. 67, No. 6 journal of December 2006 physics pp. 1155 1164 Photonic band gap engineering in 2D photonic crystals YOGITA KALRA and R K SINHA TIFAC-Center of Relevance
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationElectromagnetic cloaking by layered structure of homogeneous isotropic materials
Electromagnetic cloaking by layered structure of homogeneous isotropic materials Ying Huang, Yijun Feng, Tian Jiang Department of Electronic Science and Engineering, Nanjing University, Nanjing, 210093,
More informationPlasmonics: elementary excitation of a plasma (gas of free charges) nano-scale optics done with plasmons at metal interfaces
Plasmonics Plasmon: Plasmonics: elementary excitation of a plasma (gas of free charges) nano-scale optics done with plasmons at metal interfaces Femius Koenderink Center for Nanophotonics AMOLF, Amsterdam
More informationSupporting information. Unidirectional Doubly Enhanced MoS 2 Emission via
Supporting information Unidirectional Doubly Enhanced MoS 2 Emission via Photonic Fano Resonances Xingwang Zhang, Shinhyuk Choi, Dake Wang, Carl H. Naylor, A. T. Charlie Johnson, and Ertugrul Cubukcu,,*
More informationGuided and defect modes in periodic dielectric waveguides
Fan et al. Vol. 12, No. 7/July 1995/J. Opt. Soc. Am. B 1267 Guided and defect modes in periodic dielectric waveguides Shanhui Fan, Joshua N. Winn, Adrian Devenyi, J. C. Chen, Robert D. Meade, and J. D.
More informationPhotonic band structure in periodic dielectric structures
Photonic band structure in periodic dielectric structures Mustafa Muhammad Department of Physics University of Cincinnati Cincinnati, Ohio 45221 December 4, 2001 Abstract Recent experiments have found
More informationModeling of Kerr non-linear photonic components with mode expansion
Modeling of Kerr non-linear photonic components with mode expansion Björn Maes (bjorn.maes@intec.ugent.be), Peter Bienstman and Roel Baets Department of Information Technology, Ghent University IMEC, St.-Pietersnieuwstraat
More informationOptical cavity modes in gold shell particles
9 Optical cavity modes in gold shell particles Gold (Au) shell particles with dimensions comparable to the wavelength of light exhibit a special resonance, with a tenfold field enhancement over almost
More informationProgress In Electromagnetics Research M, Vol. 20, 81 94, 2011
Progress In Electromagnetics Research M, Vol. 2, 8 94, 2 PHOTONIC BAND STRUCTURES AND ENHANCE- MENT OF OMNIDIRECTIONAL REFLECTION BANDS BY USING A TERNARY D PHOTONIC CRYSTAL IN- CLUDING LEFT-HANDED MATERIALS
More informationare the eigenvalues of the permittivity tensor in Cartesian coordinates, written as, 3 " xx + i" xy ( , which gives
Supplemental Material for Faraday rotation enhancement of gold coated Fe2O3 nanoparticles: Comparison of experiment and theory Raj Kumar Dani, Hongwang Wang, Stefan H. Bossmann, Gary Wysin and Viktor Chikan,
More informationSUPPLEMENTARY INFORMATION
On-chip zero-index metamaterials Yang Li 1, Shota Kita 1, Philip Muñoz 1, Orad Reshef 1, Daryl I. Vulis 1, Mei Yin 1,, Marko Lončar 1 *, and Eric Mazur 1,3 * Supplementary Information: Materials and Methods
More informationSupplementary Figure 1 Schematics of an optical pulse in a nonlinear medium. A Gaussian optical pulse propagates along z-axis in a nonlinear medium
Supplementary Figure 1 Schematics of an optical pulse in a nonlinear medium. A Gaussian optical pulse propagates along z-axis in a nonlinear medium with thickness L. Supplementary Figure Measurement of
More information