Finite-difference time-domain simulations of exciton-polariton resonances in quantum-dot arrays

The optical properties of nanosize quantum-dot (QD) arrays are found to vary significantly around the exciton resonance frequency of the QDs. In order to simulate the interactions between electromagnetic waves and QD arrays, a general auxiliary-differential-equation, finite-difference time-domain approach is introduced and utilized in this article. Using this numerical method, the exciton-polariton resonances of single-layer and double-layer GaAs QD arrays are studied. The optical properties of a single-layer QD array are found to be characterized by the Mie resonance of its constituent QDs, while a double-layer QD array is characterized by the quasi-dipole formed by two QDs positioned in each of the two layers. © 2008 Optical Society of America OCIS codes: (160.4670) Optical materials; (260.5740) Physical Optics, resonance. References and links 1. A. V. Zayats, I. I. Smolyaninov, A. A. Maradudin, “Nano-optics of surface plasmon polaritons,” Phys. Rep. 408, 131 (2005). 2. E. Ozbay, “Plasmonics: Merging photonics and electronics at nanoscale dimensions,” Science 311, 189-93 (2006). 3. Y. Fu, M. Willander, E. L. Ivchenko, “Photonic dispersions of semiconductor-quantum-dot-array-based photonic crystals in primitive and face-centered cubic lattices,” Superlatt. Microstruct. 27, 255 (2000). 4. Y. Fu, E. Berglind, L. Thylén, H. Ågren, “Optical transmission and waveguiding by excitonic quantum dot lattices,” J. Opt. Soc. Am. B 23, 2441 (2006). 5. Y. Zeng, Y. Fu, X. Chen, W. Lu, H. Ågren, “Complete band gaps in three-dimensional quantum dot photonic crystals,” Phys. Rev. B 74, 115325 (2006). 6. Y. Zeng, X. Chen, W. Lu, Y. Fu, “Exciton polaritons of nano-spherical-particle photonic crystals in compound lattices,” Eur. Phys. J. B 49, 313 (2006). 7. H. Ajiki, T. Tsuji, K. Kawano, K. Cho, “Optical spectra and exciton-light coupled modes of a spherical semiconductor nanocrystal,” Phys. Rev. B 66, 245322 (2002). 8. H. Ajiki, T. Kaneno, H. Ishihara, “Vacuum-field Rabi splitting in semiconducting core-shell microsphere,” Phys. Rev. B 73, 155322 (2006). 9. H. Mertens, J. S. Biteen, H. A. Atwater, A. Polman, “Polarization-Selective Plasmon-Enhanced Silicon QuantumDot Luminescence,” Nano Lett. 6, 2622 (2006). 10. C. Wang, L. Huang, B. A. Parviz, L. Y. Lin, “Subdiffraction Photon Guidance by Quantum-Dot Cascades,” Nano Lett. 6, 2549 (2006). 11. G. Ya. Slepyan, S. A. Maksimenko, V. P. Kalosha, A. Hoffmann, D. Bimberg, “Effective boundary conditions for planar quantum dot structures,” Phys. Rev. B 64, 125326 (2001). #90132 $15.00 USD Received 27 Nov 2007; revised 23 Jan 2008; accepted 23 Jan 2008; published 18 Mar 2008 (C) 2008 OSA 31 March 2008 / Vol. 16, No. 7 / OPTICS EXPRESS 4507 12. O. Voskoboynikov, C. M. J. Wijers, J. L. Liu, C. P. Lee, “Magneto-optical response of layers of semiconductor quantum dots and nanorings,” Phys. Rev. B 71, 245332 (2005). 13. V. Bondarenko, M. Zaluzny, Y. Zhao, “Interlevel electromagnetic response of systems of spherical quantum dots,” Phys. Rev. B 71, 115304 (2005). 14. L. Belleguie, S. Mukamel, “Nonlocal electrodynamics of arrays of quantum dots,” Phys. Rev. B 52, 1936 (1995). 15. X. Zhang, P. Sharma, “Size dependency of strain in arbitrary shaped anisotropic embedded quantum dots due to nonlocal dispersive effects,” Phys. Rev. B 72, 195345 (2005). 16. F. Thiele, C. Fuchs, R. Baltz, “Optical absorption in semiconductor quantum dots: Nonlocal effects,” Phys. Rev. B 64, 205309 (2001). 17. S. Hughes, “Coupled-cavity QED using planar photonic crystals,” Phys. Rev. Lett. 98, 083603 (2007). 18. L. C. Andreani, D. Gerace, M. Agio, “Exciton-polaritons and nanoscale cavities in photonic crystal slab,” Phys. Status Solidi B. 242, 2197 (2005). 19. C. Wang, L. Y. Lin, B. A. Parviz, “Modeling and simulation for a nano-photonic quantum dot waveguide fabricated by DNA-directed self-assembly,” J. Sel. Top. Quantum Electron. 11, 500 (2005). 20. L. Huang, C. Wang, L. Y. Lin, “A comparison of crosstalk effects between colloidal quantum dot waveguides and conventional waveguides,” Opt. Lett. 32, 235 (2007). 21. S. Nojima, “Optical response of excitonic polaritons in photonic crystals,” Phys. Rev. B 59, 5662 (1999). 22. H. Haug, S. W. Koch, Quantum Theory of the Optical and Electronic Properties of Semiconductors (World Scientific, Singapore, 1986). 23. R. Zeyher, J. L. Birman, and W. Brenig, “Spatial Dispersion Effects in Resonant Polariton Scattering. I. Additional Boundary Conditions for Polarization Fields,” Phys. Rev. B 6, 4613 (1972) 24. M. Born, E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light, Pergamon press 1959. 25. A. Taflove, S. C. Hagness, Computational Electrodynamics: the finite-difference time-domain method, Second Edition, Artech House Boston 2000. 26. I. Vurgaftman, J. R. Meyer, L. R. Ram-Mohan, “Band parameters for III-V compound semiconductors and their alloys,” J. Appl. Phys. 89, 5815 (2001). 27. T. Iida, H. Ishihara, “Force control between quantum dots by light in polaritonic molecule states,” Phys. Rev. Lett. 97, 117402 (2006). 28. G. Ya. Slepyan, S. A. Maksimenko, V. P. Kalosha, J. Herrmann, N. N. Ledentsov, I. L. Krestnikov, Zh. I. Alferov, D. Bimberg, “Polarization splitting of the gain band in quantum wire and quantum dot arrays,” Phys. Rev. B 59, 12275 (1999). 29. K. Kempa, R. Ruppin, J. B. Pendry, “Electromagntic response of a point-dipole crystal,” Phys. Rev. B 72, 205103 (2005). 30. J. A. Klugkist, M. Mostovoy, J. Knoester, “Mode Softening, Ferroelectric Transition, and Tunable Photonic Band Structures in a Point-Dipole Crystal,” Phys. Rev. Lett. 96, 163903 (2006). 31. F. J. Taylor, Principles of signals and systems, McGraw-Hill, New York, 1994. 32. J. A. Pereda, L. A. Vielva, A. Vegas, A. Prieto, “Computation of resonant frequencies and quality factors of open dielectric resonators by a combination of the finite-difference time-domain (FDTD) and Prony’s methods,” IEEE Microwave Guid. Wave Lett. 2, 431 (1992). 33. Y. Hua, T. K. Sarkar, “Generalized pencil-of-function method for extracting poles of an EM system from its transient response,” IEEE Trans. Antennas Propag. 37, 229 (1989). 34. S. Dey, R. Mittra, “Efficient computation of resonant frequencies and quality factors of cavities via a combination of the finite-difference time-domain technique and the Pade approximation,” IEEE Microwave Guid. Wave Lett. 8, 415 (1998). 35. W. Guo, W. Li, Y. Huang, “Computation of resonant frequencies and quality factors of cavities by FDTD technique and Pade approximation,” IEEE Microwave and Wireless Components Lett. 11, 223 (2001). 36. G. A. Baker, Essentials of Padé approximants, Academic Press, New York, 1975. 37. V. S. C. Manga Rao, S. Hughes, “Single quantum-dot Purcell factor and β factor in a photonic crystal waveguide,” Phys. Rev. B 75, 205437 (2007).


Introduction
Miniaturization and high-density integration constitute important factors in contemporary fabrication of photonic components.The main obstacle for further progress of these factors is governed by the diffraction limit of light, something that has led to a wealth of approaches proposed to circumvent this problem.One of these have concerned the utilization of surface plasmon polaritons, that is coupled polaritons formed by photons and free electrons in a metal [1,2].Recently, quantum-dot (QD) arrays have been found to strongly manipulate electromagnetic (EM) waves whose wavelengths are one or two orders of magnitude the radius of the QDs [3,4,5,6].For instance, a single-layer array of GaAs/Al x Ga 1−x As QDs with a radius of 20 nm is found to significantly reflect an EM wave with a wavelength of about 800 nm.The QD-based optical devices could therefore provide a possible approach to decrease the size of devices and further increase the integral density of optical circuits [4,5].This possibility has now sparked an abundance of research, both experimentally [7,8,9,10] and theoretically [11,12,13,14,15,16,17,18,19,20].
The underlying physical mechanism of QD-based devices is provided by the excitonpolariton resonances of the QDs [21].Under an incident EM wave, a confined exciton can be optically excited in the QD.When a photon and an exciton interact in the dispersion-crossover region, a combined quasi-particle, normally known as an exciton polariton, is formed [22].Because excitons in a macroscopic system spread through the whole structure with an appropriate dispersion, they must be treated with a nonlocal theory [21,23].Linear response theory constitutes such a theory, giving the dielectric polarization as P(r, ω) = T (ω)Φ(r) Φ(r )E(r , ω)dr , (1) where Φ(r) := Φ(r, r) is the ground-state wavefunction of the exciton excited in QD (whose center is denoted by a), In addition, the coefficient T (ω) is given by Here ω LT and a B are the exciton longitudinal-transverse splitting and Bohr radius in the corresponding bulk semiconductor, respectively.ε b is the dielectric index of the well material, j 2 = −1, δ is a phenomenological parameter describing the decay of single-QD exciton, and R is the radius of the quantum dots.ω 0 is the ground-state exciton resonance frequency of the QDs, which is given by when the QD is assumed as a spherical square quantum well.Here E g is the band gap in the corresponding bulk semiconductor, e is the elementary charge, m e and m h are the electronic effective mass and hole effective mass, respectively.
To numerically simulate the interaction between EM waves and QD-based devices, planewave-expansion methodologies have frequently been employed.Such an approach was proposed in Ref. [3], and further extended in Ref. [5,6].However, it has turned out very difficult to simulate complicated structures and calculate the resonance modes.In this article, we introduce a general auxiliary-differential-equation, finite-difference time-domain (FDTD) approach, which is more powerful and applicable than the original plane-wave-expansion method.This method is further employed for simulations of exciton-polariton resonances in quantum-dot arrays.It should be emphasized that, although the auxiliary-differential-equation FDTD method is very popular in modeling the pulse dynamics of dispersive and nonlinear media, it is the first time, to the best of our knowledge, that this method utilized to simulate quantum dots with nonlocal dielectric polarizations.
Our paper is organized as follows: In Section II we present the numerical calculation method.Numerical results and analysis are presented in Section III.We further discuss our calculation approach and results in Section IV.Section V contains our conclusions.The quantum dot has a radius of R, the square array has a period of L, and the dielectric film has a thickness of 2R.(b) A cross section of the computational domain consisting of a single unit cell of the quantum-dot array.Periodic boundary conditions are imposed on the four surfaces perpendicular to the dielectric film, while perfect matched layers are imposed at the top and bottom surfaces.The input light wave is polarized along the z direction and propagates to the top dielectric surface along the x direction.The transmitted electric field is collected at the detector point.

Calculational method
The interaction between light and nonlocal QDs can be described by the time-dependent Maxwell's equations that are coupled to an equation for the light-induced excitonic polarization current in the QD [24], Here E is the electric field vector, H is the magnetic field vector, J is a current density corresponding to the nonlocal polarization P of the QD.To solve the above curl equations, Yee's discretization scheme is here employed [25].All field variables are defined on a cubic grid.Electric and magnetic fields are temporally separated by one-half time-step and spatially interlaced by half a grid cell.Based on this scheme, center differences in both space and time are applied to approximate Maxwell's equations.
In order to obtain the relation between the current J and the polarization P, Eq.( 1) is rewritten as where A = πε 0 ε b ω LT ω 0 , and the new variable E new (ω) is defined as with sinc(x) = sin(x)/x.It should be stressed that the dielectric permittivity of the QD is then very similar to that of a Lorentz-type medium.
Next the polarization current density J(ω) is introduced as with j 2 = −1.Fourier transforming the above equation, its time-domain analog can be written as With the discrete time step Δt, and notation J n ≡ J(nΔt), a time-difference expression is then found as where On the other hand, following the Ampere's law, the finite-difference expression can be written as Equation 10 with Eq. 13 therefore can be utilized to simulate the nonlocal polarization of the QDs in two steps: (1) From J n , H n and E n−1/2 obtain E n+1/2 ; (2) From E n−1/2 and E n+1/2 obtain E n−1/2 new and E n+1/2 new and therefore J n+1 , meanwhile, from E n+1/2 and H n obtain H n+1 .Note that the electric field E is synchronous with J in Ref. [25] while it here is separated by Δt/2.Through numerical validation, we found that our equations are stable and effective.Moreover, they can further be in favor of simplifying the following finite-difference expression of the electric field E.
In order to obtain the frequency spectrum as well as the resonance mode, a time-dependent signal is generally transformed into the frequency domain by using discrete Fourier transformation.A detailed description of this transform in FDTD is presented in Appendix A. However, due to that the exciton-polariton resonances of the QD arrays are very sharp, the frequency resolution of the spectrum obtained by the discrete Fourier transform is too poor to distinguish the resonances.We then employ the Padé approximation instead of the discrete Fourier transform to improve the accuracy of the frequency response (Appendix B).Moreover, to alleviate the calculational burden of the Padé approximation, the original time-dependent FDTD output is decimated by a factor of 1/50, to reduce the number of time samples.

Results and analysis
The geometry of the system studied computationally is shown in Fig. 1.A freestanding dielectric slab embedded with QD arrays is placed in the middle of the space with its top and bottom surfaces positioned perpendicular to the x direction.Plane waves propagating along the x axis are generated by exciting a plane of identical dipoles in phase.Perfect matched absorbing boundary conditions are applied at the top and bottom of the computational space whereas periodic boundary conditions apply on other boundaries [25].By placing one unit cell of the periodic QD array in the computational space, we can simulate the temporal transmission of the plane waves normally incident on the QD array which extends infinitely in the y and z directions.
The size of the spatial grid cell Δ is 2 nm in the numerical calculations.The total number of time steps involved in the numerical analysis is 100000 with each time step Δt = Δ/2c ≈ 3.34 attosecond, where c is the speed of light in vacuum.A time-discrete unit pulse is given as the initial excitation.The dielectric constant of the dielectric slab ε b is fixed as 12.40, and the square QD array has a period of L in the yz plane.The radius of the QDs is assumed to be uniformly 20 nm.GaAs QDs are considered here.The effective masses of electrons and holes are m e = 0.067 and m h = 0.51, respectively, in the unit of electron rest mass, and the band gap is E g = 1.51914 eV [26].The corresponding ground-state exciton resonance frequency ω 0 of QDs is then given as 1.535 eV [5].The exciton longitudinal-transverse splitting ω LT is assumed to be 0.01 eV, and the δ is assumed to be zero.
We first consider the effect of the period L on the resonance frequency of single-layer QD arrays.Three different structures are studied and the amplitude spectra (at the center of the QD) are plotted in Fig. 2. In each structure, a strong exciton-polariton resonance is observed around the ground-state exciton resonance of QD.Three important conclusions can be deduced: (1) The quality factor Q of the exciton-polariton resonance, defined by the ratio between the resonance frequency and its bandwidth ω/Δω, is very big, indicating a strong localization of the EM field inside the QD array.For instance, Q is about as high as 760 for the structure with L = 60 nm.It is for this huge Q that the Padé approximation instead of the discrete Fourier transform is employed in this paper.
(3) The bandwidth of the resonance decreases with the increase of L. This has also been observed from the temporal behavior of the polarization current density J(t) (the results are not presented here).The amplitude of the current density remains almost unchanged for the structure with L = 150 nm, while it decreases fast for L = 60 nm.The transfer of energy stored inside the QD is, therefore, much slower for the larger L.
The above observations can be explained mainly by two facts, (1) the EM field is strongly localized inside the QD but not within the dielectric slab, and (2) the coupling of two neighboring QDs, which can be approximated as dipole-dipole coupling, exponentially decreases with the increase of their interval (the period L).The optical properties of a single-layer QD array is, therefore, mainly determined by its constituent QD.Since the increase of the period L weakens the interaction of the QDs, the resonance frequency of the QD array tends to the Mie resonance frequency of an individual QD, meanwhile, the bandwidth of array decreases to the bandwidth of the Mie resonance of a single QD.Notice that Mie resonance generally means the resonance of a spherical particle under EM radiation.Here we use this term in broader contexts.
To further investigate the underlying physics of exciton-polariton resonances of single-layer QD arrays, the corresponding resonance modes are calculated.The distribution of |E 2 z | at both horizontal (x-y plane) and vertical (y-z plane) cross sections are presented in Fig. 3. Two interesting properties can be observed: (1) the |E 2 z | inside QD increases with the increase of L, due to a diminished coupling between the neighboring QDs; and (2) the EM field is observed to be strongly localized inside the QD.The optical properties of a single QD then significantly affects that of the QD array.Furthermore, the resonance modes for three structures with different L are very similar, indicating an almost identical physical mechanism.This mechanism should be the Mie resonance of an individual QD, as can be observed from the variation of |E 2 z | along the z axis (the polarization direction of the incident light), as plotted in Fig. 4. Two nodes appear at z = ±14 nm, indicating the formation of a standing wave as well as a resonance along this direction.It is important to point out that this strong localization of visible light inside QD has been already numerically observed and further employed to manipulate collective dynamics of CuCl QD (with a radius of 20 nm) [27].
The experimentally observable transmission spectrum of a single-layer quantum-dot array is also calculated and plotted in Fig. 5.The transmission coefficients are almost 100% for frequency far away from the exciton-polariton resonance, while they rapidly drop to almost zero around the resonant frequency.Similar phenomena have been numerically observed in Ref. [4], and an effective classical refractive index of quantum dot has further been employed to explain the dramatic variation of transmission spectrum.More specifically, due to the strong localization of EM field inside quantum dot around the resonance (shown in Fig. 3), the quantum dot phenomenally has a refractive index which is much bigger than that of a non-resonant quantum dot.The light is then almost totally reflected when its frequency is close to the exciton-polariton resonance frequency.
The resonance in a double-layer QD array is here studied using a period of L = 60 nm in the y − z plane, and a distance between two identical layers along the x direction (the propagation direction of the incident light) of 40 nm.Its spectrum is plotted in Fig. (6a).As a comparison, the spectrum of the corresponding single-layer array is also presented.Their resonance frequencies are very different: ω 0 + 0.726ω LT for the double-layer array, while only ω 0 + 0.404ω LT for the single-layer array.Furthermore, the resonance bandwidth of the double-layer structure is much narrower than that of the single-layer array.It should be pointed out that the QDs in the double-and single-layer structures have different local environments, namely QDs in the double-layer are completely embedded in the dielectric while those in the single-layer are not (see Fig. (6b) and Fig. (6c)).We believe that this is the main cause for the difference in the resonant frequencies and bandwidths of the QD arrays.
The resonant mode of such a double-layer array is calculated and the results are plotted in Fig. 7 and Fig. 8.The distribution of |E 2 z | is very similar to that of a dipole in character.The  electric field is very strong inside the QD of the first layer, while it is extremely weak inside the QD of the second layer.Furthermore, no nodes have been found in the profile of |E 2 z | along the z axis.The underlying mechanism of the resonances of double-layer arrays are, therefore, very different from the Mie resonance of single-layer arrays.It be largely determined by the quasi-dipole formed by two QDs, as shown by the resonance mode.
Based on the results presented above, we can conclude that, due to the excitation of exciton polariton, the optical properties of a QD array vary strongly around the ground-state exciton resonance frequency ω 0 of the constituent QDs.For example, if we consider a single array of GaAs/Al x Ga 1−x As QDs with a radius of 20 nm, the corresponding resonant wavelength is about 800 nm.An electromagnetic wave with a wavelength of about 800 nm will then be significantly manipulated by such a QD array whose thickness is as small as 40 nm.These facts have a bearing on the use of exciton polaritons for beating the diffraction limit of light, and QD arrays of the kind analyzed here may form active constituents in nanoscale photonic circuits and subwavelength components.

Discussion
The fact that QD arrays can manipulate strongly EM waves whose vacuum wavelengths are as large as almost two orders of magnitude the QD's radius may suggest the application of  effective-medium theory.For instance, the Maxwell-Garnett model with the Clausius-Mossotti correction has been utilized to study the polarization splitting of the gain band in QD arrays [28].In the framework of effective-medium theory, a complex composite material is replaced by a homogeneous medium with effective constitutive parameters.These effective parameters depend upon the generic and the geometrical parameters of the composite material.In other words, the space-dependent dielectric polarization of the original composite material is substituted by the space-independent polarization of the effective homogeneous medium.The EM field accordingly averages over material inhomogeneities.However, as shown in our results, the nonlocal polarization of the QD is very strong around resonance, and so as the EM field inside the QD.Therefore, depolarizing fields, induced by the significant difference between the dielectric constants of the QD and the host semiconductor, must be taken into account in the effective-medium theory.In contrast, the depolarizing field is inherently included in our approach, owing to the application of the exact dielectric polarization P (Eq.( 1)) of the QD.
We now move on and discuss the relation between our method and the point dipole model [29,30].When QD has a dimension much smaller than the wavelength of the EM radiation, and the inter-dot separation was large enough to avoid the wave function overlap, the QD can be approximated as a point dipole centered at its site.In other words, we neglect the spatial dependence of the QD wavefunction, and replace it by a delta function.It should be emphasized that our approach can be modified to simulate point-dipole-model QDs.Consider a QD centered at r 0 , replacing its ground-state wavefunction Φ(r) by delta function δ (r − r 0 ), its dielectric polarization is accordingly simplified as P(r, ω) = T (ω)E(r, ω)δ (r − r 0 ).Because this simplified polarization has the same characteristics of a point dipole, we can then safely expect that all the physical effects deduced from the point dipole model can be repeated by our method.
Next we discuss the effects of disorder of the QD structures.In our numerical simulations, only ideal systems are considered, that is, the radii of all QDs are assumed to be exactly identical and equal to R, and so the period L. However, certain degrees of structural disorder inevitably exist in actual experiments.The statistic distribution of the sizes of QDs can be described by a Gaussian function centered around R, and so as the resonant frequencies ω 0 .Consequently, the structural disorders lead to nonhomogeneous broadenings of the exciton-polariton resonances.The bandwidths of these resonances shown in Fig. 2 and Fig. 5 should, therefore, be much narrower than that of realistic systems.These resonances may be even completely suppressed when the structural disorders are strong enough.
To close this section we note that the work presented here, a general FDTD method, is just the first part of our whole QD-based-nanophotonics project.Many useful complex structures, such as the coupled QD-photonic-crystal-cavity [17] (and waveguide [37]) systems, will be investigated in the near future by using the numerical method introduced here.Because the optical responses of QDs in our method can be simulated more correctly than by other semiclassical models (the point-dipole model as an example), we believe that the obtained results will be quite important from both fundamental physics and applied perspectives.

Conclusion
In conclusion, an auxiliary-differential-equation, finite-difference, time-domain approach is proposed for studies of exciton-polariton resonances in quantum-dot (QD) arrays.The approach is here used to study the effect of the period of the array.Due to the excitation of exciton polaritons, QD arrays are shown to significantly manipulate light with a wavelength around the ground-state exciton resonance of the constituent QDs.The optical properties of a single-layer QD array is found to be largely affected by the Mie resonance of the constituent QDs.On the other hand, the optical properties of double-layer QD arrays is characterized by the quasi-dipole formed by two QDs positioned in each of the two layers.Because the radius of a QD is very smaller than the wavelength of its ground-state exciton resonance, exciton polaritons may offer a solution to the diffraction limit of light, and serve as a basis for constructing nanoscale photonic circuits as well as for the design and fabrication of subwavelength components.

Fig. 1 .
Fig.1.(a) Schematic drawing of a freestanding dielectric film embedded with a square array of quantum dots.The quantum dot has a radius of R, the square array has a period of L, and the dielectric film has a thickness of 2R.(b) A cross section of the computational domain consisting of a single unit cell of the quantum-dot array.Periodic boundary conditions are imposed on the four surfaces perpendicular to the dielectric film, while perfect matched layers are imposed at the top and bottom surfaces.The input light wave is polarized along the z direction and propagates to the top dielectric surface along the x direction.The transmitted electric field is collected at the detector point.

Fig. 2 .
Fig. 2. The normalized amplitude spectra (at the center of the QD) of single-layer quantumdot arrays with different period L. The dielectric film has a thickness of 40 nm.

2 Fig. 3 .
Fig. 3.The distribution of |E 2z | of single-layer quantum-dot arrays at resonant frequencies, at vertical (x-y plane) cross section and cross section (y-z plane), respectively.(a) L = 60 nm and the resonant frequency is ω 0 + 0.404ω LT .(b) L = 100 nm and the resonant frequency is ω 0 + 0.363ω LT .(c) L = 150 nm and the resonant frequency is ω 0 + 0.345ω LT .The position of quantum dot is marked by dotted lines.

Fig. 4 .
Fig. 4. Variation of |E 2 z | along (a) z axis and (b) x axis including the center of quantum dot, for single-layer quantum-dot arrays with different period L. The position of quantum dot is marked by dash dotted lines.

Fig. 5 .
Fig. 5.The transmission spectrum of a single-layer quantum-dot array.The square array has a period of L = 60 nm.

Fig. 6 .
Fig. 6.(a) The normalized spectrum of a double-layer quantum-dot array as well as that of the corresponding single-layer array.The square array has a period of L = 60 nm, and the separation between the two layers is s = 40 nm.Schematic drawing of a double-layer geometry (b) and a single-layer geometry (c).The quantum dot has a radius of R.

Fig. 7 .
Fig. 7.The distribution of |E 2 z | at a resonant frequency of ω 0 + 0.726ω LT , at (a,b) cross section (y-z plane) and (c) vertical (x-y plane) cross section, respectively.The double-layer quantum-dot array is same as that in Fig.(6b).The position of quantum dot is marked by dotted lines.

Fig. 8 .
Fig. 8. Variation of |E 2 z | along (a) x axis and (b) z axis including the center of quantum dot.The double-layer quantum-dot array is same as that in Fig.(6b).The position of quantum dot is marked by dash dotted lines.