The Field Confinement , Narrow Transmission Resonances , and Green Function of a MultilayeredMicrosphere with Metamaterial Defects

We numerically investigate the optical transmission through a compound spherical stack with conventional and metamaterial (MM) layers and also embedded MM defect. A formation of extremely narrow resonant peak with nearly complete transmittance in area of a band gap is found. We demonstrate that photon fields of certain frequencies can be strongly confined by a left-handed (LH) defect. The influence of a random deviation in the width of compound spherical layers as well the transit to the whispering gallery mode (WGM) is also discussed.


Introduction
Controlling the electromagnetic properties of materials, going beyond the limit that is attainable with naturally existing substances, has become a reality with the advent of metamaterials [1,2].The range of various structured artificial structures has promised a vast variety of otherwise unexpected physical phenomena [3][4][5] among which the experimental realization of a negative refractive index has been one of the main advantage.
A typical metamaterial combines continuous metal films (about 50 nm) with nanostrip magnetic resonators that finally yield a negative-index material.The film can be formed with a mixture of dielectric (e.g., silica) and metal, such as silver or gold, forming a semicontinuous metal film.Such structures are fabricated using techniques as the evaporation of a metal onto a dielectric substrate (see more details in review [5]).Preparation of metamaterial structure in the optical range still is an advanced task.In recent experiment [6] it has been demonstrated that the incorporation of the gain material (rhodamine 800) in the metamaterial makes it possible to fabricate an extremely low-loss and an active optical negative-index metamaterial (NIM) that is not limited by the inherent loss in its metal constituent.In this experiment the optical (NIM) structure (silver layers + rhodamine 800) is the fishnet with periodicity about 300 nm.Epoxy doped with rhodamine 800 (Rh800) and dye is used as a gain medium.The loss compensation mechanism in the sample is straightforward to understand.When Rh800 is excited by a pump pulse with sufficiently high power, a population inversion is formed inside the dye molecules.This provides amplification for a properly delayed probe pulse whose wavelength is coincident with the stimulated emission wavelength of the dye molecules (see details in [6]).
Natural extension of such directions is the analysis of a compacted spherical multilayered system with embedded metamaterial layers.It is well known that except for the whispering gallery mode (WGM) regime [7][8][9][10], a bare dielectric sphere has a complex spectrum of the electromagnetic low-quality (Q factor) eigenoscillations because of the energy leakage into the outer space [11].The case of the compound structure, when the dielectric sphere is coated by an alternative stack, is much richer.The Q-factor of optical oscillations has a large value in the frequency regions of weak transmittance and, beyond these regions Q, remains small [12][13][14].This gives rise to a large variety of optical properties of microspheres coated with a multilayer stack.Such a system can serve as a spherical symmetric photonic band gap structure, which possesses strong selective transmittance properties [15,16], with the incorporated nanometer-sized photon emitters.These possibilities essentially allow the expansion of the operational properties of microspheres with attractive artificial light sources for advanced optical technologies.The energy confinement into small volumes has applications in many fields such as photonics and quantum electrodynamics.In this sense, the microstructures that have the greatest ability to store energy for long periods of time are the dielectric spherical resonators.Recently, it has been feasible to construct such a microsphere accurately, and the parameters may be precisely controlled and measured [17].New materials are used to extend the optical properties of layered microspheres [18].
In studying such an effect in multilayered microspheres, the following question emerges: how does such an effect will be changed at a deviation from strict periodicity of a stack?Generally such a case can be interpreted as a translation symmetry breaking due to placing of defect layers in the periodic stack.The question rises whether it is possible to obtain narrow transmittance peaks in an alternating spherical stack with inserting a defect layer?
A translation symmetry breaking can be achieved in a quasiperiodic system.In [19][20][21] optical properties of layered microsphere with a dielectric stack, in which optical layers are constructed following the Fibonacci sequence, were investigated.Such a quasiperiodic structure serves as an intermediate case between a strictly periodic and cleanly casual sequence of alternative layers in a stack.It was found that when the number of layers (Fibonacci order) increases, the structure of the spectrum acquires a fractal form.If the metamaterial layers (with negative refraction index) are included in such a stack then the width of the resonant peaks in the frequency spectrum becomes extremely narrow [22].It is found the formation and coupling of defect modes in two-dimensional photonic crystal (PC) band gaps associated with degenerate edges [23].In very recent article [24] it has been shown that by varying the radii and permittivity of the defects, total transmission or reflection of the impinging electromagnetic wave can be achieved.
In this paper, we study the photon transmissions through the compound spherical stack containing conventional and metamaterial (MM) layers and also an embedded MM defect that are transparent and can bend light in opposite directions.A narrow well-separated resonant peak with nearly complete transmittance in a bandgap arises when such a MM defect layer is embedded to the center of the spherical stack.We demonstrate that photon fields of certain frequencies can be strongly confined by a left-handed (LH) defect.The influence of a random deviation in the width of compound spherical layers as well the transit to whisperinggallery mode (WGM) is also discussed.

Basic Equation
The geometry of the multilayered microsphere with the embedded defect is depicted in Figure 1.
Despite that the peculiarities of the propagation of optical waves in 1D plane systems with defects are investigated in a number of articles the results of such studies cannot be applied directly to a spherical stack for the following reasons.(i) In a spherical case, the transfer matrices M j have the determinant depending on number j of a layer as det(M j ) = (r j+1 /r j ) 2 [14] (is not unimodular), where r i+1 and r j are radii of the external and internal boundaries of the layer, respectively.Such a structure of the transfer matrix is physically caused by preservation of the energy flux in a radial direction in a solid spherical angle.(ii) In a spherical stack, the transfer matrix M j is quite involved since it is written through the complex Hankel functions.Physically, it is due to the preferential role of the pole (center) of a microsphere that breaks the translational symmetry in such a system.Furthermore, in a spherical geometry, the transmittance coefficient T depends on the spherical quantum number m (angular momentum), for great m 1 that leads to a whispering gallery mode (WGM) regime with practically zero transmittance in a frequency range.
In order to study the optical properties of a spherical stack with embedded defect layer, let us first formulate the transfer matrix method exploited in spherical multilayered geometry.Various approaches were proposed [14,[25][26][27][28] (see [29,30] and references therein).Here we follow [25].
In a multilayered microsphere, the set of the Maxwell equation is where E and B are electric and magnetic fields and ε(ω) is a dielectric permittivity of a layer.We use the complex exponential multiplier in the form exp(iωt).Equation (1) in the spherical coordinate frame (r, θ, ϕ) usually reduces to the Helmholtz equation for a scalar function called the Debye potential Π(r, θ, ϕ) [31].Since our calculations did not register significant difference between TM and TE wave cases further we will concentrate mainly in the TM case.The equation for the radial part of the Debye-potential Π = Π(r) in a layer can be readily obtained from (1) and is given by where k 0 = ω/c.Equation ( 2) is easily solved in terms of spherical Hankel functions.In each layer of the stack, we use the next matrix presentation for the fields where a and b are arbitrary constants, and matrix D = D(y) is given by inP (2)  m y e iy inP (1)  m y e −iy G (2)  m y e iy G (1)  m y e Here, n = n(ω) = ε(ω) is the refractive index of a particular layer (that can be positive or negative for metamaterial layers); P (1,2)  m (y) is the rational part of Hankel spherical functions h (1,2)  m (y) = P (1,2)  m (y)e ±iy ; G (1,2)  m (y) is the rational part of derivative of Hankel spherical functions (∂/∂y)h (1,2)  m (y) = G (1,2)  m (y)e ±iy ; m is the number of a spherical harmonic, y = ωn(ω)r/c.The recursive relations for calculations of P (1,2)  m (y) and G (1,2)  m (y) are given in [26].Below, we provide the matrix D, as well as vectors q and u by indices according to the number of layers in the stack.As the vector q in each layer is constant, for any two points r 1 and r 2 of a k-layer we obtain from (3) the following: To start the calculation, we assume that the points r 1 and r 2 belong to the boundaries of layers; r 2 = r 1 + d 1 , d 1 is the thickness of the layer.At the boundary between layers, k and k + 1, the continuity of fields gives u k (r 2 ) = u k+1 (r 2 ).With the latter (5) can be rewritten as This relation can be extended as follows.Let us start from the bottom layer of the stack with the number k = 1.Then from (6) we have where is the transfer matrix between inner and outer layers in the spherical stack.
In the simplest case of the spherical mode m = 1, the elements of matrix M k in (8) are where is the thickness of the kth layer.From ( 9), one can see that in the spherical case, the transition matrix M k depends explicitly on the distance to the center of the microsphere.
Note that in the case where r k d k or y 1 ∼ y 2 1 in ( 8), one gets the well-known plane case expression for M k [32]: Only if such an approximation is fulfilled, the matrix elements in (10) are independent of the absolute values of the coordinates and dependent only on the thicknesses of layers d k through value y 2 − y 1 = Φ k .This means that in general, the local properties of electromagnetic oscillations depend on the place of a layer with respect to the center of the spherical cavity.Equation ( 7) connects the values of the fields on both external and internal boundaries of the spherical stack.Therefore, solving the Maxwell equations is completely equivalent to calculating products of the transfer matrices M. We have used the Sommerfeld radiation conditions, where there is only an outgoing wave at the external boundary that gives q N in (3) b N = 0.As a result, the amplitudes of waves in the internal substrate become where the matrix After solving (11), one can easily obtain R and T in the next form: In (12), two equations interrelate three variables: R, T, and frequency ω.Defining ω, one can calculate the frequency dependence of R(ω), T(ω) for the stack.
We can use ( 3)-( 12) to calculate the reflectance, transmittance, eigenfrequencies, and eigenfields in an arbitrary spherical stack.Since no specific properties of a spherical stack were used, one can conclude that (3)-( 11) are valid for any structure of the spherical stack.Now, we exploit such a technique for studying a combined alternating spherical stack with conventional and metamaterial layers.
Further we consider metamaterial layers characterized by a (relative) permittivity ε(r, ω) and a (relative) permeability μ(r, ω), both of which are spatially varying, complex functions of frequency satisfying the relations [33,34] In the following, we refer to the material of a layer as being left handed (LH or metamaterial) if the real part of its refractive index is negative.In order to allow a dependence on the frequency of the refractive index, let us restrict our attention to a single-resonance permittivity and a single-resonance permeability where ω Pe , ω Pm are the coupling strengths, ω Te , ω Tm are the transverse resonance frequencies, and γ e , γ m are the absorption parameters.In our case we used the following typical values: ω Tm /2π = f Tm = 159.2THz, γ m /2π = f γ = 0.001592 THz, f Te = 163.9THz, f Pe = 119.4THz, f Pm = 68.44THz.
In the optical range the performance of all NIM applications is significantly limited by the inherent and strong energy dissipation in metals, especially in the near-infrared and visible wavelength ranges [5].In our study we have explored mainly the case where the loss compensation can be attained.It is worth noting that in general case both dependencies ε(ω) and μ(ω) may be more complicated than used in ( 14) and (15).However, since the metamaterials are artificial structures there are no strict expressions for ε(ω) and μ(ω) derived from the first principles up to now.Some analytical approximations were obtained numerically in order to fit the experiment data smoothly.In this paper we use a single (or double) resonance structure for ε(ω) and μ(ω) similarly that was detected in [6, Figure 4(b)].On the other hand we use the relations ( 14)-( 15) that already were used in a number of articles [33,35,36] and are in agreement with the principle of causality as expressed in the Kramers-Kronig relations (see recent discussion in [37,38]).We believe that such an approach cannot lead to significant errors.This allows the use in ( 14)-( 15) of the dissipation parameter γ with small value γ/2π = 0.0016 THz.
Figure 2(b) shows the LH refractive index n(r, ω) = Re n(r, ω) + i Im n(r, ω) (ω = 2π f , with the permittivity ε(ω) and the permeability μ(ω) being, resp., given by ( 14) and ( 15)).In the inset, the details of n(r, ω) are shown in the frequency interval from 164 THz where Re n(r, ω) < 0. It is worth noting that the negative real part of the refractive index is typically observed together with strong dispersion, so that absorption cannot be disregarded in general.However, in a very recent work [6], it was experimentally demonstrated that the incorporation of gain material in a metamaterial makes it possible to fabricate extremely low-loss and active optical devices.Thus, the original loss-limited negative refractive index can be drastically improved with the loss compensation in the visible wavelength range.
It is worth noting that the features of the use of tensorial or scalar ε(ω) and μ(ω) expressions for metamaterials in spherical geometry (as well for other types of curved boundaries) up to now are an open problem.Synthesized technologies are mainly elaborated for preparation of hybrid spherical layered structures (see [18] and references therein).To the best of our knowledge the separation of component ε(ω) and μ(ω) for spherical interfaces is not studied in details.However, as far as the ratio of typical sizes of the metamaterial structures (period of fishnet) and microsphere is small (about 0.1 an less) and to seek for simplicity we use the average values ε(ω) and μ(ω) that are scalar in the first approximation.

Numerical Results
Our numerical results are shown in Figures ( 2)- (10).The following parameters have been used in our calculations: the geometry of system is A, {B, C}, . . ., G, . . .{B, C}, D, where each n i a small imaginary part that corresponds to a material dissipation.
Figure 2 compares the structure of the frequency spectrum of the transmission coefficient for the spherical stack containing conventional materials, Figures 2(a) and 2(c); the case when Si is replaced by LH metamaterial is shown in Figures 2(b) and 2(d).In the latter figure it is shown the frequency range where the real part of the refraction index n(ω) of a metamaterial is negative; see an insert Figure 2(b).We observe that the transmittance spectra in Figures 2(c) and 2(d) have very different structure.Such a distinction leads to an essential change also in a frequency spectrum of spherical system with a defect layer.Further we mainly study the details of a transmittance spectrum for a spherical system with the LH defect embedded.
Figure 3(a) shows the case of the λ/4 spherical stack without a defect, whereas in Figure 3(b) the situation when a LH defect is embedded in the center of the spherical stack is depicted.We observe from Figures 3(a) and 3(b) that the bandgaps become less for frequencies that are closer to the dispersive area with larger values of the refraction index |n|.(It is worth noting that in very recent experiment [40] demonstrated that a broadband, extremely high index of refraction can be realized from large-area, free-standing, flexible terahertz metamaterials composed of strongly coupled unit cells.)Besides, from Figure 3(b) we observe that the narrow well-separated resonant peaks with nearly complete transmittance arise in the center of the forbidden frequency zone.
Figure 3 represents only main spectral features of multilayered spherical system with an embedded LH defect.Further we study the details of the transmittance peak with respect to the change of width, position, and the field spherical quantum number.Another important feature considered here is the dependence of the peak details on the random deviation in the widths of the stack layers.Figure 4 shows the details of transmittance peak in the frequency band gap when a defect is embedded in different positions in the stack.The system consists of 15 alternating λ/4 layers LH + SiO 2 and the embedded LH defect layer.We observe from Figure 4 that the transmittance resonance reaches the maximum value for a case when the defect is placed in center of the stack.
The behavior of transmittance peak in the frequency band gap of the alternating 15-layered λ/4 spherical stack, at variations of the defect layer width (the defect is embedded in the center of the stack, layer 8) is shown in Figure 5.Our calculations have shown that the position of the peak strongly depends on the defect layer width d and does not depend significantly on the material.From Figure 5 we observe that for defect thickness d = λ/4 the peak is situated in the center of the band gap, and it migrates to the lower frequencies range when d increases.
In order to investigate the features of such an effect closely to the WGM regime we studied the dependence of the peak width on the spherical quantum number m; see Figure 6.The alternation spherical stack is composed by 15 alternating layers LH, SiO 2 , and the LH defect layer is embedded in the center of the stack.Such a behavior can be easily explained as follows.At a large m, there is a transition to the WGM (whispering gallery mode) which occurs when the optical field is exponentially localized near the surface of the microsphere.In this regime, the field does not "feel" the details of the spherical stack structure.Thus, with the increase m the amplitude of the peak will decrease rapidly since at a large m, the spatial structure of the layers does not essentially affect the spectrum of the optical oscillations.
Figure 7(a) studies the vicinity of some "defect" resonances with narrow spectral width that have transmittance T close to 1. On the other hand it is interesting to explore The localization of the defect photon mode can be explained as following [42].Once a photon enters the defect region, it encounters two λ/4 Bragg reflectors (the periodic parts of the stack) before and behind it.This leads to the fact that the photon will be strongly reflected back to the defect region and thus remains long time in the defect area.Such a long dwell time results in very high energy field density around the defect.
In above we have paid more attention to the purely quarter-wave stack case; however, we now study what happens when the thicknesses of the layers in the stack have some random deviations (technological roughness) from the ordered λ/4 width.In this case, the far-order periodicity in the system is broken, and each particular layer has an influence on the total properties of the stack.
Figure 8 shows that random deviations in the width layer of the λ/4 stack results in the change of the resonance position.One can see that despite the value of the random deviation level (a) 0.5%, (b) 2%, (c) 3%, and (d) 9% the peak band remains very narrow.We can explain that as follows.For a range of frequencies in the band gap, the field is confined already in the vicinity of the defect layer and farther is considerably reduced.As a result, the field in the other layers of the stack practically does not affect the structure of the confined field in the defect layer.Such behavior appreciably differs from the case of a stack with strong alternating layers (without defects).In the later case the formation of the spectral peaks is a collective effect due to the wave re-reflections from the great number of individual layers.In such structures small random deviations in a layers' width can be statistically recompensed and result in the destruction of the frequency resonances.
From the above, we can conclude that the transmittance spectrum of spherical structures with an embedded LH defect is rather stable with respect to random deviations.The considered confinement of the field by the defect layer is a typical internal property of such nonuniform metamaterial structures.

Green Function and Spectrum of Nanoemitters Radiation
In the above we have studied the transmittance properties of the compound spherical stack with a defect, particularly at separated m.Other question of fundamental interest is the details of optical spectrum of nanoemitter radiating from such a coated microsphere.The latter is formed by a superposition of re-reflecting optical waves with different m.Such a study already requires more advanced formalism, namely, the Green function approach.In this section we apply the Green function technique to study the frequency spectrum of a nanoemitter radiation.The spatial scale of the nanoemitter objects (∼1-100 nm) is at least one order of magnitude smaller than the spatial scale of microspheres (∼10 3 -10 4 nm).Therefore in the coated microsphere (Figure 1), we can represent the nanoemitter structure as a point source placed at r and having a dipole moment d 0 .It is well known that the solution of the wave equation for the radiated electromagnetic field E due to a general source J(r ) is [33,43,44] where G(r, r , ω) is the dyadic Green function (DGF), which depends on the type of boundary conditions imposed on E(r) and contains all the physical information necessary for describing the multilayered structure (the time dependence is assumed to be e iωt ).Equation ( 16) is complemented by the standard boundary conditions: limitation of the fields in the center of the microsphere and continuity of the tangential components of the fields at the interfaces of layers.We also use Sommerfeld's radiation conditions, where there is only an outgoing wave in the external boundary of the microsphere.In this case, the electromagnetic field E in the coated structure consists of the sum of the waves radiating in the surrounding medium and the multiple wave reflections due to the interfaces between layers.Substituting the nanoemitter source in the form J = iωd, d = d 0 δ(r − r ) in ( 16), we obtain where p 0 = (μd 0 /ε 0 )(ω 2 /c 2 ).In such a situation, the nanoemitter frequency spectrum is identical to the dyadic Green's function (DGF) spectrum.Thus, the equation of the field generated by a nanoemitter assumes the form of the DGF G(r, r , ω) equation and is given by [43,44] ( where κ(r, ω) = μ −1 (r, ω); r is the point where the field is observed, while r is the nanoemitter (point source) location.
The dyadic Green's function G(r, r , ω) also satisfies the boundary conditions at the interfaces of spherical layers [43,45] Let us consider the multilayered spherical structure: a concentric system of spherical layers contacting with the sphere (concentric stack) deposited onto the surface of the microsphere with nanoemitters placed in such a structure (see Figure 1).The layers are localized at the distance R k from the center, where d k = R k − R k+1 is the width of a kth layer.
Let us first specify some details of the Green's function technique for multilayered microspheres and introduce our notations.Following the approach [45], we write down DGF of such a system as follows: where G V 0e represents the contribution of the direct waves from the radiation sources in the unbounded medium, whereas G ( f s) es describes the contribution of the multiple reflection and transmission waves due to the layer interfaces.The dyadic Green's tensor G V 0e (r, r , ω) in ( 20) is given by where tensor G V q,nm (r, r , ω) in ( 21) has a form q,nm(r, k s )M q,nm (r , k s ) + N (1)  q,nm(r, k s )N q,nm (r , k s ), r > r M q,nm (r, k s )M (1)  q,nm(r , k s ) + N q,nm (r, k 1 )N (1)  q,nm(r , k s ), r < r In (22), vectors M and N represent TE-and TM-waves, where the prime denotes the nanoemitter coordinates (r , θ , ϕ ), n and m are spherical and azimuthal quantum numbers, respectively, while k s is the wave number of the medium where the radiated nanoemitters are located.It is worth noting that due to the dyad r r, the δ-function in (21) contributes to the radial (longitudinal) part [43].Due to the equality r • (θ θ + ϕ ϕ) = 0 (where r, θ, and ϕ are the unitary vectors of the spherical coordinates) such a singularity does not contribute to the field (17) for case of the tangential dipole.Since nanoemitters (nanorods) are highly polarized objects, we pay more attention to the case when the dipole orientation of nanoemitter is d = d ϕ, so only the tangential components of the Green's tensor G ϕϕ contribute.Description of the Green function approach and complete expression for the scattering Green tensor G ( f s) es (r, r , ω) is rather bulky and can be found in [22,28,45].
It is well known that due to the fluctuation-dissipation theorem the correlation function of the photon states in the absorbing environment for temperature T can be written down with the help of the macroscopic Green-function as follows [46]: In general a number of spherical harmonics contribute to the spectrum of the emitter's radiation in a microsphere.We investigate the frequency distribution of the field energy W = Im( G ϕϕ (r, r, f )) at partial spherical numbers m in order to see the transit to the WGM regime.In Figure 10 the structure of W peaks at different spherical numbers m, in a vicinity of the metamaterial resonances, is shown.The multilayered spherical stack is composed by 12 alternating layers LH, SiO 2 , and the LH defect layer is embedded in the center of the stack.We observe from Figure 10 the existence of a strong frequency peak in area f ∼ 170.5 THz.However, the number of such peaks decreases for larger m 1 that corresponds to the transit to the WGM regime.Such a behavior can be explained as follows.At large m (that is an index of the spherical Bessel function j m (z) in DGF) the following asymptotic expression takes place [47] j m (z) (2πm) −0.5 (ez) m /2m.Thus, at some z = r/r ext < 1 (r ext is the external radius) the optical amplitudes will decrease rapidly for deeper r r ext and large m.This results in concentration of the photonic field closely to the external boundary of microsphere; therefore field does not "feel" the details of multilayered spherical structure.This leads to a diminution of the boundary re-reflections and phase interferences, which normally are observed as frequency resonant peaks.

Conclusion
We have investigated optical transmission through a compound spherical stack with conventional and metamaterial (MM) layers and also embedded MM defect.We have found that incorporation of such a defect allows a formation of extremely narrow resonant peaks with nearly complete transmittance in area of a band gap.It is found a strong dependence of the spectral position of the resonance on the width of the defect layer.The latter in principle opens new possibilities creating optical filters with extremely narrow passbands in the desired frequency range.We demonstrate that the optical fields of certain frequencies can be arrested by a left-handed (LH) defect.The latter can be used to obtain resonant field structures by engineering defects in microspheres with a compound conventional and metamaterial dielectric stack.

Figure 2 :
Figure 2: (Color online) Comparison of the spectral properties of the transmission coefficients T for 15-layered alternating stacks containing conventional (Si + SiO 2 ) and (LH + SiO 2 ) materials correspondingly (spherical number m = 1): (a) the refractive index n(ω) ( f = ω/2π) for the Si (dotted line) and SiO 2 (solid line); (b) the real (dotted line) and imaginary part (solid line) of the refraction coefficient n(ω) for LH metamaterial.The spectrum T(ω) is shown in (c) for conventional stack, while (d) shows T(ω) for a combined stack when Si is replaced with LH layers.See details in text.

Figure 3 : 2 Figure 4 :Figure 5 :
Figure 3: (Color online) The change of the spectrum shape of transmission coefficient T for alternating 15-layered λ/4 spherical stack containing LH and SiO 2 layers (spherical number is m = 1) without and with a defect: (a) strictly periodic case; (b) the stack with a defect LH placed in the center of the stack.The frequency range in both panels corresponds to the frequency range where the real part of LH refraction index is negative; see the inset of Figure 2(b).

Figure 6 :
Figure 6: (Color online) The dependence of the transmission coefficients T on the spherical quantum numbers m for the periodic multilayers microsphere system composed of 15 alternating layers of LH and SiO 2 and a defect of LH placed in the center of the system: (a) m = 2; (b) m = 3; (c) m = 4; and (d) m = 6.

Figure 7 :
Figure 7: (Color online) The field structure for various resonances (numerated as f 1,2,3 ) produced by a defect in the forbidden frequency gaps.(a) shows the spectrum of transmission coefficient T for the periodic multilayers microsphere system composed of 15 alternating layers of LH and SiO 2 and a LH defect placed in the center of the stack.In (b), (c), and (d) show the radial dependence of the field (arbitrary units) on r/r 1 corresponding to the resonant peaks marked in (a) as f 1 , f 2 , f 3 , respectively.The dashed lines in (b), (c), and (d) show the structures of the spherical stack (refractive indices).The spherical number is m = 1, r 1 is radius of the internal bottom microsphere.

Figure 8 :
Figure 8: (Color online.)The transmittance peak position for a microsphere coated by 15 alternating layers and the LH defect embedded in the center for cases: (a) 0.5%, (b) 2%, (c) 5%, and (d) 9% random deviations in the width of all the layers of the spherical stack from the exact λ/4 case.

Figure 9 :
Figure 9: (Color online) (a) The frequency spectrum of the imaginary part of the Green function W = Im(G ϕϕ (r, r, f )), (20) (arbitrary units); (b) real and imaginary parts of the LH material refractive index n in the frequency area closely to the metamaterial resonances.

Figure (
Figure (Color online) The frequency spectrum of the imaginary part of the Green function W = Im(G ϕϕ (r, r, f )) (arbitrary units) for partial spherical quantum numbers m in the area closely to the metamaterial resonances.
r)E(r ) ω = Dω 2 c 2 coth Dω 2T Im G(r, r , ω) .(23) In particular, from (23) one can see that the case r = r yields the inequality Im( G(r, r)) > 0 that corresponds to the energy of a fluctuating electromagnetic field E(r) at small dissipation [46]: E(r) 2 ∼ Im( G(r, r)), where G(r, r) is the Green function, that is, G ϕϕ (r, r) in our case.The typical structure of the Green function spectrum is shown in Figure 9(a).We observe various peak frequency resonances in the area where real part of the LH material is negative Re(n) < 0 (see Figure 9(b)).However no resonances one observes in the frequency zone f ≥ 173 THz, where Re(n) ≥ 0.