# Change of structure of the Cherenkov emission at modulated source in dispersive metamaterials.

1. INTRODUCTIONNanophotonics is the study of the behavior of light on the nanometer scale with involving the interaction of light with nano-structures. Novel optical properties of materials results from their extremely small size that have a variety of applications in nanophotonics and plasmonics. The investigations of optical negative-index [1] metamaterials (NIM) using the nanostructured metal-dielectric composites already have led to both fundamental and applied achievements that have been realized in various structures [2-26]. The main applications of negative index metamaterials (or left-handed materials (LHM)) are connected with a remarkable property: the direction of the energy flow and the direction of the phase velocity are opposite in NIM that results unusual properties of electromagnetic waves propagating in these mediums.

Cherenkov radiation by a charged source that moves in a left-handed material and has not the own frequency has been studied in number of works [17-23]. Both experimental and theoretical frameworks are investigated, see review [18] and references therein. However, in some important cases a moving particle has the internal own frequency [27] [[omega].sub.0], for instance the ion oscillating at the transition frequency [[omega].sub.0]. The Cherenkov emission in metamaterials for such situation is insufficiently investigated yet, although it is a logical extension of other works in this area.

In this paper, the Cherenkov optical radiation in 3D metamaterials by a nonrelativistic modulated source having the own frequency [[omega].sub.0] with an emphasis on the dispersive properties of the medium is numerically studied. We performed the FDTD simulations with the use of the material parameters at various modulating frequencies [[omega].sub.0], however, without references to the operational frequency range.

2. BASIC EQUATIONS

In metamaterials, it is necessary to treat electromagnetic wave interactions with a metal ingredient using a dispersive formulation that allows correct description of the internal electron dynamics. It is well-known that realistic LHM (where for a range of frequency the negative refraction index n is expected) can be characterized by either Lorentz or Drude dispersion models, see, e.g., Ref. [25], Chapters 5, 6 and references therein. In case of a small collision frequency the study of nonlinear dynamical excitations is allowed [26]. In our paper, the opposite case is considered when the low frequency edge of magnetic forbidden band [25] [[omega].sub.m0] is small enough with respect to the magnetic collision frequency [[gamma].sub.m] [much greater than] [[omega].sub.m0]. The latter allows reducing the Lorenz model to the Drude model that is widely used for modeling in complex materials. The Maxwell equations read

[nabla] x E = -[[mu].sub.0][[mu].sub.h] [[partial derivative]H/[partial derivative]t] ~ [J.sub.m] - [[sigma].sub.m]H, (1)

[nabla] x H = [[epsilon].sub.0][[epsilon].sub.h] [[partial derivative]E/[partial derivative]t] + q[v.sub.0]f (r,t) cos ([[omega].sub.0]t) + [J.sub.e] + [[sigma].sub.e]E, (2)

where a radiating particle (bunch) has modulating frequency [[omega].sub.0], [J.sub.e] is the electrical current and [J.sub.m] is the magnetic current which obey the following material equations

[[??].sub.e] + [[gamma].sub.e][J.sub.e] = [b.sub.e]E, [[??].sub.m] + [[gamma].sub.m][J.sub.m] = [b.sub.m]H, (3)

here [[gamma].sub.e] and [[gamma].sub.m] are the electrical and magnetic collision frequencies respectively, be = [[epsilon].sub.0][[omega].sup.2.sub.pe]; [b.sub.m] = [[mu].sub.0][[omega].sup.2.sub.pm], [[omega].sub.pe] and [[omega].sub.pm] are frequencies of electric and magnetic plasmons respectively; [[sigma].sub.e] and [[sigma].sub.m] are conductivities; [[epsilon].sub.h], [[mu].sub.h] are dielectric and magnetic functions of the host medium respectively [28,29]. For metals such as silver, gold, copper and aluminum the density of the free electrons is on the order of [10.sup.23] [cm.sup.-3]. The typical value [[omega].sub.pe] [approximately equal to] 2 x [10.sup.16][s.sup.-1] ([30], p.44). In a metamaterial with fishnet structure [14] we consider the charge particle (charge q) moving with a uniform velocity parallel to x direction: [v.sub.0][parallel][[??].sub.x] and the density of the particle is defined by the Gaussian as f (r,t) = [W.sup.-3]exp{-[[(x - [[upsilon].sub.0]t).sup.2] + [y.sup.2] + [z.sup.2]]/[W.sup.2]}, where W is the width; at W [right arrow] 0 such a distribution is simplified to the isotropic point-source distribution f(r,t) [right arrow] [([pi]).sup.3/2][delta](x - [[upsilon].sub.0]t)[delta](y)[delta](z). In the following for simulations we use dimensionless variables, where for renormalization are used: the vacuum light velocity c = [([[epsilon].sub.0][[mu].sub.0]).sup.-0.5] and the typical spatial scale for nanooptics objects [l.sub.0] = 75 nm. With such a normalization above indicated the metal plasma frequency becomes [[omega].sub.pe] = 5. The electrical and magnetic fields are renormalized with the electrical scale [E.sub.0] = [ql.sub.0][[epsilon].sub.0] and magnetic scale [H.sub.0] = [([[epsilon].sub.0]/[[mu].sub.0]).sup.0.5][E.sub.0] respectively. Some metamaterials exhibit anisotropic properties with tensor permittivity and permeability. To seek for simplicity, in this paper, we concentrated in the isotropic geometry. Modeling anisotropic medium is a straightforward extension of this model, see details in [18, 23].

In our FDTD simulations (in time domain) we use the following idea. In optical experiments normally one refers only to the parameters of material ([[gamma].sub.e], [[sigma].sub.e], [[omega].sub.pe], and [[gamma].sub.m], [[sigma].sub.m], [[omega].sub.pm]). Therefore, we consider the Cherenkov emission produced by a particle with modulating frequency [[omega].sub.0], with the use only the material parameters and without of references to the operational frequency band. In this approach, the frequency spectrum of internal excitations [omega] is left as a free parameter that is defined from simulations by a self-consistent way. For 3D dispersion material such a problem becomes too difficult for analytical consideration. In this paper, the well-known numerical algorithms in the time domain [31] were used, for details see Ref. [11].

General 3D case in Cartesian coordinates is considered since such a geometry normally is used on the optical investigations [14]. We examine a spatially averaged metamaterial composition: nanostructured metal-dielectric composites (fishnet), similarly that was used in the experiment [14]. In this case, the spatial average scale is less then the infrared and visual wavelengths, so we can deem that the material (dielectric and magnetic) dispersion is allowed by the Drude model and the role of the active dielectric ingredient is reduced to a compensation of losses due to the metal ingredient. In our simulations was used the numerical grid [L.sup.3], L = 100, 120, 150, for more details see Ref. [11].

The following steps have been used in our approach: (i) First, we calculated the time-spatial field dynamics that is raised by the crossing radiating particle for different modulating frequencies [[omega].sub.0], Eqs. (1)-(3). (ii) In second step we applied the Fourier analysis for the time and spatial dependencies of the field calculated in the first step in order to reveal the dynamics and spectrum of internal excitations. The following dimensionless parameters were used in our simulations: [[omega].sub.pe] = 5, [[omega].sub.pm] = 7, [[epsilon].sub.h] = 1.44, [[mu].sub.h] = 1, [[gamma].sub.e] = [[gamma].sub.m] = [10.sup.-4], [[sigma].sub.e] = [[sigma].sub.m] = [10.sup.-7], W = 4, [q.sub.2] = 2. We also varied the particle velocity [[upsilon].sub.0] for different [[omega].sub.0] to study regimes of the Cherenkov emission. Our results are shown in Figs. 1-5.

Figure 1 shows the spatial structure of [E.sub.x](r, t) field for the particle with modulating frequency [[omega].sub.0] that moves with uniform velocity [[upsilon].sub.0] = 0.52. Fig. 1(a) depicts the case [[omega].sub.0] = 0; in this situation only the spectrum of the periodic plasmonic excitations with the frequency [[omega].sub.C] and negative phase velocity is observed (see details in Ref. [11]). At the increase of [[omega].sub.0] up to 6 (see Figs. 1(b), (c) and (d) respectively) the spatial field structure changes considerably. For [[omega].sub.0] = 2 and [[omega].sub.0] = 6 (see Figs. 1(b) and (d) respectively) the field has oscillating shape, but for [[omega].sub.0] = 4 in Fig. 1(c) the field acquires a near monotonic form.

In what follows we use the complementary Cherenkov emission angle [[theta].sub.1] = [pi]/2 - [theta], where [28] cos([theta]) = c/n[[upsilon].sub.0]. For such angle we have [11] cos([theta]) = sin([[theta].sub.1]) = c/n[[upsilon].sub.0], so for conventional material (with Re(n) > 0) [[theta].sub.1] is positive [[theta].sub.1] = [[theta].sub.1+] > 0, while for negative refraction index metamaterial NIM (with Re(n) < 0) [[theta].sub.1] is negative [[theta].sub.1] = [[theta].sub.1-] < 0. From Fig. 1 we observe that for frequencies [[omega].sub.0] [less than or equal to] 3 the angle [[theta].sub.1] is negative what corresponds to the reverse Cherenkov radiation. But for larger [[omega].sub.0] [greater than or equal to] 6 the Cherenkov radiation already acquires conventional structure with positive [[theta].sub.1] > 0, see Fig. 1(d). More interesting is found the intermediate case with [[omega].sub.0] [equivalent] 4 where the field acquires nearly monotonic shape, see Fig. 1(c). Such a behavior can be interpreted as a transition in the spectra of plasmonic-polariton excitations. Further, we study such field behavior with details.

The interesting question emerges, why the structure of field in Fig. 1 changes so considerably with the change of the modulating frequency [[omega].sub.0]? To study that we explore the variation of the field spatial structure in the wavenumber domain (k) at the change [[omega].sub.0]. It is instructively to compare the field spatial structure with [[omega].sub.0] = 0 in NIM with other modulated particle case when [[omega].sub.0] [not equal to] 0 but in a dispersiveless material (with [[omega].sub.pe] = 0, [[omega].sub.pm] = 0). (Such an evaluation in principle could be done similarly to the Lienard-Wiechert potentials approach [28], but analytical calculations of the field structure are difficult, so in this paper, we evaluate the field dependencies numerically in the framework of our model.) Such comparison is presented in Fig. 2 for the time t = [t.sub.0] ~ 190 when the particle reaches the output side of the system. From Fig. 2 we observe that for dispersiveless and dispersive cases (left and right panels in Fig. 2 respectively) the fields [E.sub.x](x,t) have well-defined resonances in the wavenumber domain [E.sub.x](k,t) at [k.sub.s] [equivalent] [k.sub.s0] [equivalent] 4.33 (resonance point). One can conclude that in this point a resonant interplay arises between the plasmonic-polariton spectrum and the modulating particle field. Corresponding changes are resonantly accumulated within the system length that finally leads to considerable modification of the field shape along the particle path.

Now we turn to other cases of the velocity [[upsilon].sub.0] and modulation frequency [[omega].sub.0]. To understand the effect deeper it is important to investigate the Doppler shift of the frequency [[omega].sub.0] with respect of the rest metamaterial. For dispersiveless case (n = const) the shifted frequency is calculated explicitly as [[omega].sub.s0] = [[omega].sub.0]/(1 + [[upsilon].sub.0]n/c). However, in dispersive case the shifted frequency [[omega].sub.s0] acquires dependence on the dispersive refraction index of medium n([omega]). For this situation, the shift already has to be evaluated from the complete equation [[omega].sub.s0] = [[omega].sub.0]/(1 + [[upsilon].sub.0]n([[omega].sub.s0])/c). The latter already requires the formulation of the dispersive material model for n([omega]). But in used FDTD approach (the time domain) such the shift should be formed by a self-consistent way in NIM model, where in some frequency ranges the negativity of Re(n([omega])) can be reached. (In general the Doppler frequency shift has a nonlinear dependence on the modulating frequency, see details in Ref. [12].)

It is readily to see that [k.sub.s] peak is calculated from the Doppler shift as [k.sub.s] = [[omega].sub.s0]n/c = [[omega].sub.0]n/c(1 + [[upsilon].sub.0]n/c). To verify if the observed [k.sub.s] is related to the Doppler shifted frequency [[omega].sub.0] in Fig. 3 the dependence [k.sub.s] = [k.sub.s]([[upsilon].sub.0]) for various frequencies [[omega].sub.0] in the dispersiveless case is shown. From Fig. 3 we observe that (in used normalization) the wavenumber peak is about [k.sub.s] [equivalent] 4.33 for [[omega].sub.0] [equivalent] 4, [[upsilon].sub.0] = 0.52 that coincides with peak [k.sub.s0] = 4.33 for no-modulated source at [[upsilon].sub.0] = 0.52 for NIM in Fig. 1(a). Thus, we conclude that for [[omega].sub.0] [equivalent] 4 and [[upsilon].sub.0] [equivalent] 0.52 a spatial resonance arises due to Doppler effect around of [absolute value of [k.sub.s] - [k.sub.s0]].

Finally, we study the frequency spectra of the Cherenkov emission for different values of [[omega].sub.0] and [[upsilon].sub.0]. In the frequency domain the dispersive permittivity and permeability have the following form ([31] Chapter 9) [epsilon]([omega]) = [[epsilon].sub.h] - [[omega].sup.2.sub.pe]/([[omega].sup.2] + i[[gamma].sub.e][omega]) and [mu]([omega]) = [[mu].sub.h] - [[omega].sup.2.sub.pm]/([[omega].sup.2] + i[[gamma].sub.m][omega]). For such [epsilon]([omega]) and [mu]([omega]) the complex refraction index for the NIM metamaterial can be written as [33]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

In this case, a peak frequency [omega] = [[omega].sub.C] has to be substituted into [epsilon]([omega]), [mu]([omega]) and then in n([omega]) Eq. (4). Fig. 4 (left and right panels respectively) show the structure of field [E.sub.x](r,t) and the complex refractive indices n([omega]) in NIM for the cases different particle modulating frequencies [[omega].sub.0] and velocity [[upsilon].sub.0] = 0.52. We observe from Fig. 4(a) that the field structures corresponds to the inverse Cherenkov emission. Fig. 4(c) shows that for [[omega].sub.0] = 2 the plasmon-polaritons excitation (PPE) are generated at the peak frequency [[omega].sub.C] = 3.77 where the complex refractive index is n([[omega].sub.C]) = -0.87 - i8.10 x [10.sup.-5]. For larger frequency [[omega].sub.0] = 3.9 the field [E.sub.x] has monotonic spatial structure, see Fig. 4(b) that allows exciting PPE at the peak frequency [[omega].sub.C] = 3.83 where n([[omega].sub.C]) = -0.77 - i8.09 x [10.sup.-5]. In both cases, Re(n) < 0, thus, the optical waves have negative phase velocity that corresponds to reverse Cherenkov emission in NIM.

The situation for smaller v0 is shown in Fig. 5. Fig. 5 (left and right panels respectively) show the structure of field [E.sub.x](r,t) and the complex refractive indices n([omega]) in NIM for the cases of particle with different modulating frequencies [[omega].sub.0] and particle velocity [[upsilon].sub.0] = 0.35. We observe from Fig. 5(c) that for [[omega].sub.0] = 2 the plasmon-polaritons excitation (PPE) are generated at the peak frequency [[omega].sub.C] = 3.76 where the complex refractive index is n = -0.89 - i8.0 x [10.sup.-5].

Figures 4 and 5 show that the considered spatial resonance between the modulating frequency [[omega].sub.0] and the plasmonic-polariton excitations: (i) happens in the frequency range where Re(n) is negative and the inverse Cherenkov emission occurs, and (ii) the value of the metamaterial refraction index (and thus the polariton phase velocity) depends on the particle velocity [[upsilon].sub.0]. The latter in principle allows to control the property of PPE in metamaterial.

3. CONCLUSION

We numerically studied the Cherenkov optical emission by a nonrelativistic modulated source crossing 3D dispersive metamaterial. It is found that the resonant interaction of the field produced by the modulated source with the spectrum of the periodic plasmonic-polariton excitations leads to considerable change of spatial structure for the Cherenkov emission. The field acquires monotonic shape in the frequency range where the dispersive refractive index of metamaterial is negative and the reversed Cherenkov radiation is generated. This effect opens new interesting possibilities in various applications metamaterials in nanophotonics with the potential for creating and control light confining structures to considerably enhance light-matter interactions.

ACKNOWLEDGMENT

The work is partially supported by CONACyT grant 169496.

REFERENCES

[1.] Veselago, V. G., "The electrodynamics of substances with simultaneously negative values of [epsilon] and [mu]," Sov. Phys. Usp., Vol. 10, 509, 1968; Usp. Fiz. Nauk, Vol. 92, 517-526, 1967.

[2.] Shalaev, V. M., "Optical negative-index metamaterials," Nature Photonics, Vol. 1, 41-48, 2007.

[3.] Soukoulis, C. M. and M. Wegener, "Past achievements and future challenges in the development of three-dimensional photonic metamaterials," Nature Photonics, Vol. 5, 523-530, 2011.

[4.] Hess, O., J. B. Pendry, S. A. Maier, R. F. Oulton, J. M. Hamm, and K.L. Tsakmakidis, "Active nanoplasmonic metamaterials," Nature Materials, Vol. 11, 573-584, 2012.

[5.] Chen, H., C. T. Chan, and P. Sheng, "Transformation optics and metamaterials," Nature Materials, Vol. 9, 387-396, 2010.

[6.] Gordon, J. A. and R. W. Ziolkowski, "CNP optical metamaterials," Opt. Express, Vol. 16, 6692-6716, 2008.

[7.] Milton, G. W., "Realizability of metamaterials with prescribed electric permittivity and magnetic permeability tensors," New Journal of Physics, Vol. 12, 033035, 2010.

[8.] Podolskiy, V., A. Sarychev, and V. Shalaev, "Plasmon modes and negative refraction in metal nanowire composites," Opt. Express, Vol. 11, 735-745, 2003.

[9.] Shalaev, V. M., W. Cai, U. K. Chettiar, H.-K. Yuan, A. K. Sarychev, V. P. Drachev, and A. V. Kildishev, "Negative index of refraction in optical metamaterials," Opt. Lett., Vol. 30, No. 24, 3356-3358, 2005.

[10.] Burlak, G., A. D-de-Anda, R. S. Salgado, and J. P. Ortega, "Narrow transmittance peaks in a multilayered microsphere with a quasiperiodic left-handed stack," Optics Commun., Vol. 283, No. 19, 3569-3577, 2010.

[11.] Burlak, G., "Spectrum of Cherenkov radiation in dispersive metamaterials with negative refraction index," Progress In Electromagnetics Research, Vol. 132, 149-158, 2012.

[12.] Burlak, G. and V. Rabinovich, "Time-frequency integrals and the stationary phase method in problems of waves propagation from moving sources," Symmetry, Integrability and Geometry: Methods and Applications, Vol. 8, 096, 21, 2012.

[13.] Burlak, G. and A. D-de-Anda, "The field confinement, narrow transmission resonances and Green function of a multilayered microsphere with metamaterial defects," Journal of Atomic, Molecular, and Optical Physics, Article ID 217020, 1-13, 2011.

[14.] Xiao, S., V. P. Drachev, A. V. Kildishev, X. Ni, U. K. Chettiar, H.-K. Yuan, and V. M. Shalaev, "Loss-free and active optical negative-index metamaterials," Nature, Vol. 466, 735-738, 2010.

[15.] Deb, S. and S. D. Gupta, "Absorption and dispersion in metamaterials: Feasibility of device applications," J. Phys., Vol. 75, No. 5, 837-854, 2010.

[16.] Cherenkov, P. A., "Visible emission of clean liquids by action of [gamma]-radiation," Dokl. Akad. Nauk., Vol. 2, 451-454, 1934.

[17.] Averkov, Y. O. and V. M. Yakovenko, "Cherenkov radiation by an electron particle that moves in a vacuum above a left-handed material," Phys. Rev. B, Vol. 79, 193402-193412, 2005.

[18.] Duan, Z. Y., B. I. Wu, S. Xi, H. S. Chen, and M. Chen, "Research progress in reversed Cherenkov radiation in double-negative metamaterials," Progress In Electromagnetics Research, Vol. 90, 75-87, 2009.

[19.] Xi, S., H. Chen, T. Jiang, L. Ran, J. Huangfu, B.-I. Wu, J. A. Kong, and M. Chen, "Experimental verification of reversed Cherenkov radiation in left-handed metamaterial," Phys. Rev. Lett., Vol. 103, 194801, 2009.

[20.] Averkov, Y. O., A. V. Kats, and V. M. Yakovenko, "Electron beam excitation of left-handed surface electromagnetic waves at artificial interfaces," Phys. Rev. B, Vol. 72, 205110-205114, 2005.

[21.] Zhou, J., Z. Duan, Y. Zhang, M. Hu, W. Liu, P. Zhang, and S. Liu, "Numerical investigation of Cherenkov radiations emitted by an electron beam particle in isotropic double-negative metamaterials," Nuclear Instruments and Methods in Physics Research Section A, Vol. 654, No. 1, 475-480, 2011.

[22.] Duan, Z. Y., Y. S. Wang, X. T. Mao, W. X. Wang, and M. Chen, "Experimental demonstration of double-negative metamaterials partially filled in a circular waveguide," Progress In Electromagnetics Research, Vol. 121, 215-224, 2011.

[23.] Duan, Z., C. Guo, and M. Chen, "Enhanced reversed Cherenkov radiation in a waveguide with double-negative metamaterials," Opt. Express, Vol. 19, 13825-13830, 2011.

[24.] Zhu, L., F.-Y. Meng, F. Zhang, J. Fu, Q. Wu, X. M. Ding, and J. L.-W. Li, "An ultra-low loss split ring resonator by suppressing the electric dipole moment approach," Progress In Electromagnetics Research, Vol. 137, 239-254, 2013.

[25.] Hao, Y. and R. Mittra, FDTD Modeling of Metamaterials: Theory and Applications, Artech House, 2009.

[26.] Gabitov, I. R., R. A. Indik, N. A. Litchinitser, et al., "Double-resonant optical materials with embedded metal nanoparticles," J. Opt. Soc. Am. B, Vol. 23, 535-542, 2006.

[27.] Ginzburg, V. L., "Radiation by uniformly moving sources (Vavilov-Cherenkov effect, transition radiation, and other phenomena)," Phys. Usp., Vol. 39, 973-982, 1996.

[28.] Jackson, J. D., Classical Electrodynamics, John Willey and Sons, 1975.

[29.] Oughstun, K. E., Electromagnetic and Optical Pulse Propagation 2: Temporal Pulse Dynamics in Dispersive, Attenuative Media (Springer Series in Optical Sciences), Springer, 2009.

[30.] Yeh, P., Optical Waves in Layered Media, John Wiley and Sons, New York, 1988.

[31.] Taflove, A. and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method, Artech House, Boston, 2005.

[32.] Afanasiev, G. N., Cherenkov Radiation in a Dispersive Medium, Vavilov-Cherenkov and Synchrotron Radiation, Fundamental Theories of Physics, Kluwer Academic Publishers, 2004.

[33.] Ziolkowski, R. W., "Superluminal transmission of information through an electromagnetic metamaterial," Phys. Rev. E, Vol. 63, 046604, 2001.

Gennadiy Burlak * and Erika Martinez-Sanchez

Centro de Investigacion en Ingenieria y Ciencias Aplicadas, Universidad Autonoma del Estado de Morelos, Cuernavaca, Mor., Mexico

Received 20 March 2013, Accepted 20 April 2013, Scheduled 28 April 2013

* Corresponding author: Gennadiy Burlak (gburlak@uaem.mx).

Printer friendly Cite/link Email Feedback | |

Author: | Burlak, Gennadiy; Martinez-Sanchez, Erika |
---|---|

Publication: | Progress In Electromagnetics Research |

Article Type: | Report |

Geographic Code: | 1MEX |

Date: | Jul 1, 2013 |

Words: | 3706 |

Previous Article: | A compact inductively loaded monopole antenna for future UWB applications. |

Next Article: | A novel power divider integrated with SIW and DGS technology. |

Topics: |