Abstract
We investigate the phenomena of negative refraction and backward wave in pseudochiral mediums, with illustrations of Gaussian beams. Due to symmetry breaking intrinsic in pseudochiral mediums, there exist two elliptically polarized eigenwaves with different wave vectors. As the chirality parameter increases from zero, the two waves begin to split from each other. For a wave incident from vacuum onto a pseudochiral medium, negative refraction may occur for the right-handed wave, whereas backward wave may appear for the left-handed wave. These features are illustrated with Gaussian beams based on Fourier integral formulations for the incident, reflected, and transmitted waves. Negative refraction and backward wave are manifest, respectively, on the energy flow in space and wavefront movement in time.
© 2013 Optical Society of America
1. Introduction
Pseudochiral mediums are a special type of bianisotropic materials [1], whose properties are characterized by the magnetoelectric tensors [2]. The chirality parameters are purely imaginary as in chiral mediums. Their diagonal elements, however, are zero as in ordinary dielectrics. Therefore, the magnetoelectric couplings in pseudochiral mediums occur at mutually perpendicular directions rather than the parallel direction [3]. This feature is related to some optical properties that are different from either isotropic chiral or ordinary dielectric materials. In practice, the underlying magnetoelectric couplings in pseudochiral mediums can be achieved by surface currents induced on Ω-shaped metal wires [4, 5].
Due to symmetry breaking intrinsic in bianisotropic mediums, there exit two eigenwaves with elliptically polarizations. These waves are considered the most general form of eigenwaves in complex mediums [6]. From the perspective of constitutive relations, the chirality parameters in pseudochiral mediums play dual roles: magnetoelectric coupling and anisotropy. The former is the cross coupling effect that the polarization (magnetization) in the medium can be induced by the magnetic (electric) field. The latter is a directionally dependent property that causes the electric field no longer perpendicular to the wave vector, or the wave vector no longer parallel to the Poynting vector. The wave propagation in pseudochiral mediums thus behave in part as isotropic chiral mediums and in part as anisotropic dielectrics. In particular, negative refraction and backward wave that occur in anisotropic [7–9], biisotropic (or isotropic chiral) [10–12], and bianisotropic (or anisotropic chiral) materials [13, 14] may appear as well in pseudochiral mediums.
In this study, we investigate the phenomena of negative refraction and backward wave in pseudochiral mediums. In order to characterize these features in a more concise manner, we study the propagation of Gaussian beams based on Fourier integral formulations [15]. The eigenwaves and Fresnel equations for the reflection and transmission coefficients are employed to formulate the beams in terms of their spectrsal components. In practice, the resulting infinite integrals are performed by the Fourier transforms and accelerated by the FFT algorithm. Due to the magnetoelectric couplings in pseudochiral mediums, the transmitted waves are split into two waves with different handednesses. For a beam incident from vacuum onto a pseudochiral medium, negative refraction may occur for the right-handed beam, whereas backward wave may appear for the left-handed beam. These features are manifest on the power flow of beam centers in space and the movement of beam wavefronts in time.
2. Basic equations
2.1. Eigenwaves
Bianisotropic materials are an important class of complex mediums [6], which are characterized by the following constitutive relations:
where ε, μ, ξ, and ζ are in general tensors of complex quantities. In particular, the magneto-electric tensors for pseudochiral mediums have the following forms: with γ being a real quantity. Consider the case where the anisotropy of the permittivity and permeability be small compared with the effect of chirality, and the material parameters be approximated by ε = εI and μ = μI, where I is the identity tensor. This is the case when the magnetoelectric coupling is much more significant than the dielectric-magnetic property of the medium. The dispersion relation for pseudochiral mediums is given as (see Appendix A for details) where and . Assume that ε > 0, μ > 0, and γ2 < ε μ. Therefore, 0 < γ̃ < 1 and 0 < ρ < 1. For a single frequency ω, there exist two wave numbers k± = |k±| in the pseudochiral medium.Let the xy plane be an interface between vacuum and a pseudochiral medium. Consider a wave incident from vacuum with the wave vector lying on the xz plane. For a given tangential wave vector component kx, there are two solutions for the normal wave vector component [cf. Eq. (4)]:
where . The eigenwaves in the pseudochiral medium are given as (see Appendix A for details) where and are the electric field magnitudes for right-handed elliptically polarized (REP) and left-handed elliptically polarized (LEP) waves, respectively, , , and The time-averaged Poynting vectors, , where * denotes the complex conjugate, are given as The angles of wave vectors k± and Poynting vectors 〈S±〉 with respect to the interface normal (the z axis) are given, respectively, as For a wave incident from vacuum onto a pseudochiral medium at an angle of incidence θ, negative refraction will occur for the REP wave (ϕ+ < 0°) when θ < θNR, where is the threshold angle for negative refraction. Note that when |kx| = γ̃h0, the time-averaged Poynting vector 〈S+〉 is directed to the interface normal (σ+ = 0) [cf. Eq. (8)]. On the other hand, backward wave will occur for the LEP wave (θ− > 90°) when θ > θBW, where is the threshold angle for backward wave. Note that when |kx| = ρh0, the wave vector k− is oriented along the interface ( = 0) [cf. Eq. (5)].2.2. Gaussian beams
Suppose that the xy plane is an interface between vacuum on the left (z < 0) and a pseudochiral medium on the right (z > 0), characterized by the permittivity ε, permeability μ, and the chirality parameter γ. Consider a Gaussian beam incident from vacuum, with the beam center making an angle θ with respect to the interface normal, as schematically shown in Fig. 1. Let the xz plane be the plane of incidence and kx be the wave vector component along the interface. The Gaussian beam with the center located at x = x0 and z = −h is well approximated by the Fourier integral on kx as [15]
where is the wave number component normal to the interface, and w0 is the waist size of the Gaussian beam. Based on this formulation, the incident beams are formulated as where are the electric field magnitudes for the p- or s-polarized incident waves, , , and , are the basis vectors of the incident electric and magnetic fields, respectively.The reflected beams are formulated as
where , and , are the basis vectors of the reflected electric and magnetic fields, respectively, and Rp,s are the reflection coefficients, given as and . Detailed formulas for rpp, rps, rsp, and rss are listed in Appendix B.In the pseudochiral medium, there are two wave numbers: k± [cf. Eq. (4)], corresponding to the REP and LEP waves. The transmitted beams are formulated as
where , and are the basis vectors of the transmitted electric and magnetic fields, respectively [cf. Eq. (6)], and T± are the transmission coefficients, given as and . Detailed formulas for t+p, t+s, t−p, and t−s are listed in Appendix B.If the incidence beams are circularly polarized, we have
where and are the electric field magnitudes for the right-handed circularly polarized (RCP) and left-handed circularly polarized (LCP) incident waves, respectively, and and are the basis vectors of the incident electric and magnetic fields, respectively. The corresponding reflected beams are formulated as where and are the basis vectors of the reflected electric and magnetic fields, respectively, and R± are the reflection coefficients, given as and . Detailed formulas for r++, r+−,r−+, and r−− are listed in Appendix B. The corresponding transmitted beams have the same forms as Eqs. (18) and (19), except that the transmission coefficients are changed to and . Detailed formulas for t++, t+−, t−+, and t−− are listed in Appendix B.2.3. Fourier transforms
The Gaussian beams formulated as Fourier integrals in Eqs. (14)–(23) can be carried out by performing the Fourier transform (FT) and inverse Fourier transform (IFT) as
where is the field distribution of the Gaussian beam at the interface (z = 0), ϕ (kx) is the function of kx that contains the components of basis vectors ( , , , ) and the reflection and transmission coefficients (Rp,s, R±, T±), qz = kz (−kz) for the incident (reflected) beam, and for the transmitted beam. FT and IFT are defined, respectively, asIn Eq. (24), FT[f0 (x)] is to obtain the spectral components of the fields at the interface, and IFT is to restore these components to the physical quantity f (x,z), taking into account the wave propagation in the z direction and the influence of the interface (reflection and transmission). In practice, the calculations of FT and IFT can be significantly accelerated by employing the FFT algorithm.
After obtaining the fields of the incident, reflected, and transmitted beams, the power intensity is given as I = |〈S〉|, where is the time-averaged Poynting vector for the total electric field E and total magnetic field H in vacuum or the pseudochiral medium.
3. Results and discussion
Figure 2(a) shows the normalized power intensity for a p-polarized (TM) Gaussian beam incident from vacuum onto a pseudochiral medium with ε/ε0 = εr = 2, μ/μ0 = μr = 1, and γ = 0.2. The incident beam center makes an angle θ = 20° with respect to the interface normal (the z axis). There are two transmitted beams in the pseudochiral medium, one is the REP wave with a smaller angle of the beam center, and the other is the LEP wave with a larger angle. The beam centers coincide with the time-averaged Poynting vectors 〈S±〉 [cf. Eq. (8)], with ϕ+ ≈ 6° and ϕ− ≈ 21.2°. Note that in this case, θ > θNR ≈ 11.5° [cf. Eq. (10)] and ordinary refraction is expected to occur. As the chirality parameter γ increases, the REP beam tends to move toward the interface, whereas the LEP beam tends to move away from it. In this configuration, the reflection coefficients rpp and rsp are relatively small. Compared to the incident beam, the intensity of the reflected beam is insignificant.
In Fig. 2(b), the chirality parameter is increased to γ = 0.8, with the other material parameters unchanged. The REP (lower) beam is located on the same side of the incident beam, with a negative angle of refraction: ϕ+ ≈ −19.8°. Note that in this case, θ < θNR ≈ 53.1° [cf. Eq. (10)] and negative refraction is expected to occur. The LEP (upper) beam has an angle of refraction: ϕ− ≈ 37.6°, which is larger than that in Fig. 2(a). For a s-polarized incident Gaussian beam, the features stated above are similar.
Figure 3(a) shows the incidence of a RCP Gaussian beam from vacuum at θ = 20° onto a pseudochiral medium with the same material parameters as in Fig. 2(b). In this case, the transmission is dominated by the REP beam, which is negatively refracted with the same ϕ+ as in Fig. 2(b). Compared to the REP beam, the power intensity of the LEP beam is not significant (0.037%). In Fig. 3(b), the incident beam is changed to a LCP Gaussian beam. The transmission, on the other hand, is dominated by the LEP beam, which is positively refracted with the same ϕ− as in Fig. 2(b). Compared to the LEP beam, the power intensity of the REP beam is not significant (0.037%).
A special feature of the Gaussian beam incident from vacuum at θ = 0° onto a pseudochiral medium is shown in Fig. 4(a). The transmitted wave is split into two beams even at normal incidence. This feature does not appear in either isotropic chiral or anisotropic dielectric materials, and can be manifest on the time-averaged Poynting vectors 〈S±〉 at kx = 0. In this situation, σ± = ∓γ̃ [cf. Eq. (7)] and ϕ± = arctan(∓γ̃) [cf. Eq. (9)]. In the presence of chirality parameter γ̃, ϕ+ ≠ ϕ− at θ = 0°. The effect of γ̃ on ϕ± is plotted in Fig. 4(b). The splitting angle between the two transmitted beams increases with γ̃. This feature is further analyzed by examining the basis wave functions that compose the Gaussian beams:
Using Eq. (5) for , the above equation can be integrated to give where and erfi (x) = erf (ix)/i is the imaginary error function. In Eq. (29), f± (x,z) attain their local maximal values along x = ±γ̃z, which are the locations of beam centers at θ = 0°.Finally, an instance of backward wave is shown in Fig. 5 ( Media 1) for a s-polarized Gaussian beam incident from vacuum at θ = 25° onto the pseudochiral medium with εr = 1, μr = 1, and γ̃ = 0.95. In this case, the angles of wave vectors with respect to the interface normal for the REP and LEP waves are θ+ ≈ 31.7° and θ− ≈ 105.7°, respectively. The latter represents a backward wave. The wavefronts of the LEP (upper) beam move toward the interface rather than away from it, which is indicated by the orientation of k−. The angles of time-averaged Poynting vectors with respect to the interface normal for the REP and LEP waves are ϕ+ ≈ −38.8° and ϕ− ≈ 47.6°, respectively. The former corresponds to negative refraction. Note that the threshold angles are θNR ≈ 71.8° and θBW ≈ 18.2° [cf. Eqs. (10) and (11)]. Accordingly, both negative refraction and backward wave occur in the same configuration.
4. Concluding remarks
In conclusion, we have investigated the phenomena of negative refraction and backward wave in pseudochiral mediums, with illustrations of Gaussian beams. The physical origin of these features comes from the symmetry breaking intrinsic in pseudochiral mediums. There exist two elliptically polarized eigenwaves, with the wave vectors no longer parallel to the Poynting vectors. The right- and left-handed waves behave in a very different manner, which is more evident as the chirality parameter increases. These features are well illustrated with Gaussian beams based on Fourier integral formulations for the incident, reflected, and transmitted waves.
Appendix
A. Dispersion relation and eigenwaves
Using the constitutive relations (1) and (2) in Maxwell’s equations: ∇ ×E = iωB, ∇ × H = −iωD and eliminate D and B, we may obtain the following separate equations for E and H fields:
which are regarded as the wave equations for bianisotropic mediums. Assume that E and H are of the form eik·r, the wave equations are rewritten as ME = 0 and NH = 0, where with I being the identity tensor. The existence of nontrivial solutions for E and H fields requires that |M| = 0 and |N| = 0. Using k = (kx, ky, kz) and the chirality parameters (3), the zero determinant gives rise to the characteristic equation: where and . The above equation is then solved to give the dispersion relation: The nullspace of M or N gives the eigenwaves as where , , , andB. Reflection and transmission coefficients
Consider a planar interface between vacuum and a pseudochiral medium. The continuity of the tangential electric and magnetic fields at the interface gives rise to the following relations between the incident, reflected, and transmitted waves as
where with , , and . If the incident wave is circularly polarized, the relations between the incident, reflected, and transmitted waves are given by whereAcknowledgments
This work was supported in part by National Science Council of the Republic of China under Contract No. NSC 99-2221-E-002-121-MY3.
References and links
1. J. A. Kong, “Theorems of bianisotropic media,” Proc. IEEE 60, 1036–1046 (1972). [CrossRef]
2. A. Serdyukov, I. Semchenko, S. Tretyakov, and A. Sihvola, Electromagnetics of Bi-anisotropic Materials: Theory and Applications (Gordon and Breach, 2001).
3. S. A. Tretyakov, A. H. Sihvola, A. A. Sochava, and C. R. Simovski, “Magnetoelectric interactions in bi-anisotropic media,” J. Electromagn. Waves Appl. 12, 481–497 (1998). [CrossRef]
4. M. M. I. Saadoun and N. Engheta, “A reciprocal phase shifter using novel pseudochiral or s medium,” Microwave and Optical Technology Letters 5, 184–188 (1992). [CrossRef]
5. S. A. Matos, C. R. Paiva, and A. M. Barbosa, “Surface and proper leaky-modes in a lossless grounded pseudochiral omega slab,” Microw. Opt. Technol. Lett. 50, 814–818 (2008). [CrossRef]
6. W. S. Weiglhofer and A. Lakhtakia, Introduction to Complex Mediums for Optics and Electromagnetics (SPIE, 2003). [CrossRef]
7. P. A. Belov, “Backward waves and negative refraction in uniaxial dielectrics with negative dielectric permittivity along the anisotropy axis,” Microw. Opt. Technol. Lett. 37, 259–263 (2003). [CrossRef]
8. L. Yonghua, W. Pei, Y. Peijun, X. Jianping, and M. Hai, “Negative refraction at the interface of uniaxial anisotropic media,” Opt. Comm. 246, 429–435 (2005). [CrossRef]
9. H. Chen, S. Xu, and J. Li, “Negative reflection of waves at planar interfaces associated with a uniaxial medium,” Opt. Lett. 34, 3283–3285 (2009). [CrossRef] [PubMed]
10. J. B. Pendry, “A chiral route to negative refraction,” Science 306, 1353–1355 (2004). [CrossRef] [PubMed]
11. C. Monzon and D. W. Forester, “Negative refraction and focusing of circularly polarized waves in optically active media,” Phys. Rev. Lett. 95, 123904 (2005). [CrossRef] [PubMed]
12. S. Tretyakov, A. Sihvola, and L. Jylha, “Backward-wave regime and negative refraction in chiral composites,” Photonics Nanostruct. 3, 107–115 (2005). [CrossRef]
13. Q. Cheng and T. J. Cui, “Negative refractions in uniaxially anisotropic chiral media,” Phys. Rev. B 73, 113104 (2006). [CrossRef]
14. S. A. Tretyakov, C. R. Simovski, and M. Hudlicka, “Bianisotropic route to the realization and matching of backward-wave metamaterial slabs,” Phys. Rev. B 75, 153104 (2007). [CrossRef]
15. B. R. Horowitz and T. Tamir, “Lateral Displacement of a light beam at a dielectric interface,” J. Opt. Soc. Am. 61, 586–594 (1971). [CrossRef]