Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Incident-angle-insensitive and polarization independent polarization rotator

Open Access Open Access

Abstract

This paper proposes a method to design an incident-angle-insensitive polarization-independent polarization rotator. This polarization rotator is composed of layers of impedance-matched anisotropic metamaterial (IMAM) with each layer’s optical axes gradually rotating an angle. Numerical simulation based on the generalized 4 × 4 transfer matrix method is applied, and the results reveal that the IMAM rotator is not only polarization-independent but also insensitive to the angle of incidence. A 90° polarization rotation with tiny ellipticity variation is still available at a wide range of incident angles from 0 to 40°, which is further confirmed with a microwave bi-split-ring resonator (bi-SRR) rotator. This may be valuable for the design of optoelectronic and microwave devices.

©2010 Optical Society of America

1. Introduction

As an important property of transverse EM waves, polarization has been widely applied in engineering and scientific researches. Many approaches have been employed to manipulate the polarization of light. Traditional anisotropic crystals and chiral liquid crystals [1] are commonly used as wave retarders and polarization rotators. However, the polarization rotator made of anisotropic crystal is generally polarization-dependent and incident-angle-sensitive, and the polarization rotator composed of liquid crystals is also incident-angle-sensitive. These shortcomings limit their applications in microwave and optoelectronic devices. To our knowledge, how to design and fabricate an incident-angle-insensitive and polarization-independent polarization rotator still remains a challenge.

Thanks to the advent of metamaterials, this gives a chance to overcome the challenge. Metamaterials are artificial structures, which usually have periodic arrangements and exhibit exotic electromagnetic properties [28]. These manmade structures provide completely new mechanisms and novel methods to control light. Early and ongoing researches on metamaterials have shown that it is possible to obtain strong anisotropy or chirality via deliberate design and fabrication. This can be used in the design of polarization devices. One of them is polarization beam splitter achieved by anomalous reflection and transmission [912]. Another application is the polarization rotator based on modes coupling, extraordinary optical transmission (EOT), chirality of the structure, or phase mutation at resonance frequency [5,1317]. Up till now, wide-angle polarizer [18], splitter [10,11] and absorber [19] have been realized with these artificial structures. In this paper, we proposed a method to design an incident-angle-insensitive and polarization-independent polarization rotator. The polarization rotator is composed of layered impendence-matched anisotropic metamaterials (IMAM) with the crystal axes rotated, as shown in Fig. 1 . The cross polarization conversion becomes polarization-independent and insensitive to the incident direction of the light beam. It can be proved that only two layers of IMAMs are enough to construct a cross-polarization rotator, which greatly alleviates the difficulty in fabrication. It might be valuable for the design of optoelectronic and microwave devices.

 figure: Fig. 1

Fig. 1 Schematic structure of a multilayered metamaterial rotator. The plane of incidence is in x-z plane. ϕ is the angle between the axes of the adjacent layers. φ and θ are the polarization angle and incident angle of the incident wave, respectively.

Download Full Size | PDF

2. General formalism of IMAM polarization rotator

To study the transmission and reflection of plane waves in an anisotropic slab, a generalized 4 × 4 transfer-matrix method is applied [20]. By transforming the permittivity and permeability of IMAM with a rotation matrix: ε2=ε11, μ2=μ11, where

=(cosϕsinϕ0sinϕcosϕ0001),
one can easily obtain the transfer matrix of an anisotropic slab with its crystal axes rotated along z axis by an angle ϕ. When many layers are aligned together, the total transfer matrix is given by
P¯¯total(H)=P¯¯m(hm)P¯¯m1(hm1)P¯¯1(h1),
where H=i=1mhi,which is the total thickness of the stratified slabs. The 2 × 2 Jones transmission coefficient matrix can be acquired after some tedious algebraic operation that is not presented here [21]. To simplify the discussion, we restrict the incidence in x-z plane, i.e. ky=0,and one of the crystal axes is chosen as z axis.

When the polarization of the transmitted wave is perpendicular to a linearly-polarized incident wave (EiM,EiE)=E0(cosφ,sinφ),one can obtain

Et(sinφcosφ)=E0(TMMTMETEMTEE)(cosφsinφ),
where φ is the polarization angle of the incident wave. TMMand TEEare the co-polarization transmission coefficients of TM and TE waves, respectively; while the cross terms are the cross-polarization transmission coefficients. By multiplying (cosφ,sinφ) on both sides of the equation, we have

0=E0[TMMcos2φ+TEEsin2φ+(TME+TEM)sinφcosφ].

For a polarization-independent rotator, the following condition must be satisfied: TMM=TEE0, TEM+TME0.Generally, it cannot be fulfilled with a single anisotropic medium but with a layered structure. To give an intuitionistic discussion, we choose an IMAM slab with the following permittivity and permeability tensors:

εr=(εxx000εyy000εzz),
μr=(μxx000μyy000μzz),
where
εyy/μxx=εxx/μyy=εsur/μsur=1/Z.
Zis the wave impedance, andεxxεyy,μxxμyy, εsur and μsurare permittivity and permeability of the surrounding material. When the light strikes normally at the interface between the lossless surroundings and an IMAM slab with thicknessh,the transmission coefficients of TE and TM waves can be simplified as
TE=exp(iεyy'k0hZ)exp(εyy''k0hZ),
TM=exp(iμyy'k0h/Z)exp(μyy''k0h/Z),
wherek0is the wave vector in vacuum, andε' (μ') and ε'' (μ'') are the real and imaginary parts of permittivity (permeability). When
εyy''=εxx''=μxx''/Z2=μyy''/Z2,
TE and TM waves experience the same attenuation, i.e. |TE|=|TM|1.

Since the impendence-matched condition is satisfied, reflections at the interface either between slab and surroundings or between slabs disappear. Then we can give an explicit expression on the transmission coefficient of the stratified slabs. When the plane wave (EiM,EiE) is normally incident onto x-y plane, the projections of the transmitted wave polarization on the crystal axes (denoted by suffix 1, 2) of the (n + 1)th slab can be written as

(E1<n+1>E2<n+1>)=(TEcosϕTEsinϕTMsinϕTMcosϕ)n×{(TEcosδTEsinδTMsinδTMcosδ)(EiMEiE)}.

As shown in Fig. 2 , δ is the angle between the x axis and the crystal axis 1 of the first layer; ϕ is the angle of the crystal axes between adjacent layers. We once again project (E1<n+1>,E2<n+1>) onto the coordinate axes and the Jones matrix can be obtained:

 figure: Fig. 2

Fig. 2 Schematic view of the polarization vector projections on x-y plane when the light is normally incident. The crystal axes of the anisotropic slabs are rotated along z direction by an angleϕin sequence.

Download Full Size | PDF

(TMMTMETEMTEE)=(cos(nϕ+δ)sin(nϕ+δ)sin(nϕ+δ)cos(nϕ+δ))×{(TEcosϕTEsinϕTMsinϕTMcosϕ)n(TEcosδTEsinδTMsinδTMcosδ)}.

By carefully tuning the thickness of IMAM slab that makes it as a half-wave retarder, one can obtainTE=TM=T. When n is odd, the Jones matrix can be simplified as follows:

(TMMTMETEMTEE)=Tn+1(cos[(n+1)ϕ]sin[(n+1)ϕ]sin[(n+1)ϕ]cos[(n+1)ϕ]).

It is clear that, if T is not considered, the Jones matrix is a coordinate transformation matrix by rotating an angle of (n+1)ϕ anticlockwise. Hence the polarization rotation angle(n+1)ϕ can be tuned dynamically and is independent of φ and δ. When (n+1)ϕ=±π/2,

(TMMTMETEMTEE)=Tn+1(0±110).

This is the condition for the cross polarization conversion. It is clear that only two (n = 1) IMAM slabs are enough to construct a polarization-independent rotator.

Equation (14) indicates that the eigen polarization states of the rotator are left and right circular polarizations. Although it is the same for planar chiral structures with four-fold symmetry, it is very difficult to achieve linear polarization rotation in planar metallic chiral structures due to dichroism; while in dielectric chiral structures, one can achieve linear polarization rotation but this property strongly depends on incident direction, as has been pointed out [22].

3. Incident-angle-insensitivity and polarization independency

It is more interesting to find that the IMAM rotator is not only polarization-independent, but also incident-angle-insensitive. This is distinct from conventional anisotropic or chiral rotators. Without loss of generality, we choose air as surroundings, and a lossless IMAM bilayered rotator with parameters asεxx=μyy=2, εyy=μxx=1, εzz=μzz=1.The thickness h of each layer is chosen ash/λ0=1.5 so that TE=TM=exp[iϑ(kx)] is satisfied at normal incidence and each layer of slab acts as a half-wave retarder. To construct a polarization-independent rotator, two half-wave plates are aligned together in the above mentioned way, i.e. 2ϕ=90°. To study the effects of incident direction on the performance of the polarization rotator, incident angle θ (denoted bykx/k0), angle δ and incident polarization angle φ should be considered. Figures 3(a) to 3(d) show the polarization changes of the transmitted light. The ellipticity Δ and the polarization rotation angleψtof the transmitted light are defined as

Δ=Ea/Eb,
ψt=φtφ,
where the electric field amplitudesEa and Eb are the minor and major polarization components of the transmitted light, respectively; φt is the transmitted polarization angle.

 figure: Fig. 3

Fig. 3 The projections of (a) the ellipticity Δ and (b) the polarization rotation angle ψt with respect to kx/k0 and δ when φ = 45°. The projection views of (c) Δ and (d) ψt with respect to δ and φ at 40° incident angle, in which φ ranges from −90° to 90°. The insets demonstrate the three-dimensional relations. The lossless IMAM bilayered polarization rotator is characterized withεxx=μyy=2, εyy=μxx=1, εzz=μzz=1,and a total thickness 2h/λ0=3.

Download Full Size | PDF

Since we have restricted the incident beam in x-z plane, the variation of δ is actually equivalent to the change of incidence in polar direction. Now we look about how the performance of the polarization rotator depends on δ when φ = 45°. Figure 3(a) is the projection of the three-dimensional (3D) curvature of ellipticity Δ as a function of δ and kx/k0 onto the Δ~kx/k0 plane; Fig. 3(b) is the projection of the 3D curvature of polarization rotation angle ψt as a function of δ and kx/k0 onto the ψt~kx/k0 plane. It is clear that, at normal incidence (kx/k0=0), Δ0 andψt90, the bilayered structure acts as a perfect cross-polarization rotator. As the incident angle increases, the variation ranges of Δ and ψt become larger but still acceptable within a broad range of incident angles. It can be seen in Fig. 3(a) and Fig. 3(b) that when δ ranges from 0° to 360°, Δ(0.077,0.077) (indicating a value of 0.006 as an intensity contrast of two polarization components) and ψt (90.69°,89.05°) even at a 40° incident angle (kx/k00.64). As the incident angle decreases to 30°, variations can drop toΔ(0.024,0.024) and ψt (90.4°,89.6°).

Figures 3(c) and 3(d) depict the effect of incident polarization angle φ on the polarization rotation at the oblique incidence of 40°. The insets plot the variations of Δ and ψt with respect to δ and φ at a 40° incident angle. The projections of these curvatures onto Δ~δ and ψt~δ planes are shown in these figures. When φ changes from −90° to 90°, the maximum variation range of Δ is (0.077,0.077),while the maximum variation of ψt is (90.69°,89.05°).After studying the effects of incident direction under variouskx/k0, φ and δ, we confirm that the IMAM rotator is more insensitive to incident direction when the incident angle is smaller.

Moreover, we studied the effect of material loss. As is expected, the loss has a small impact on the performance of the IMAM rotator that satisfies Eq. (10). We set εyy''=εxx''=μxx''=μyy''=0.2 while other parameters maintain the same. The results show that for a lossy rotator, the corresponding curvatures of Figs. 3(a), 3(b) and 3(c) are similar to the lossless ones; while for the one correspond to Fig. 3(d), the variation range of ψt is enlarged to (−91.66°,-88.37°). Even so, a less than 2° deviation in polarization rotation at 40° incidence is still quite acceptable.

4. Comparison and discussion

To make a comparison, we choose three lossless anisotropic slabs (surrounded by air) with parameters as εxx=2b, εyy=εzz=1, μyy=2/b, μxx=μzz=1, and the thickness h of each layer is chosen as h/λ0=1.5, where b=1,2,0.5 respectively. Only the case b = 1 is impedance-matched among the three cases. Since (εxxμyy)1/2k0h= 6π for TM waves and (εyyμxx)1/2k0h = 3π for TE waves, TE=TM=exp[iϑ(kx)] is satisfied for all three half-wave plates at normal incidence; when two identical plates are stacked together in the aforementioned way, i.e. 2ϕ=90°,they all act as polarization-independent rotators at normal incidence. However, the impedance-matched conditions are not satisfied at oblique incidence, and this will bring in negative influences on the performance of the rotators. Nevertheless, as is shown above, the structure composed of IMAM slabs (b = 1) is more insensitive to incident angle.

Figures 4(a) to 4(d) demonstrate the cross-polarization transmission coefficients and polarization rotations with respect to kx/k0 whenδ=45°, φ = 45°. It is clear that for the impendence-matched polarization rotator, the incident angle region with |TEM||TME|1 (|TEE|and |TMM| tend to zero, which are not shown here) and TEMTME=±π is much wider than the impendence-mismatched ones. Thus, the polarization rotation angle |ψt|90° and the polarization ellipticity |Δ|0 can be realized within a much wider incident angle as well. For the IMAM polarization rotator, when the incident angle varies from 0° to 40° (kx/k00.64), Δ(0.077,0.077) and ψt (90.69°,89.05°). While for the impendence-mismatched rotators, the polarization rotation and ellipticity become very instable as the incident angle increases. To take the nonmagnetic rotator with b = 2 as an example, the polarization rotation becomes −105° and the ellipticity is around-0.6 at a 40° incidence.

 figure: Fig. 4

Fig. 4 (a) The amplitudes of transmission coefficients |TEM|and |TME|, (b) the phase difference TEMTME, (c) the polarization rotation angle ψt and (d) the ellipticity Δ with respect to the incident transversal vector kx/k0 for three different anisotropic bilayered structures characterized as εxx=2b, εyy=εzz=1, μyy=2/b, μxx=μzz=1, and the total thickness 2h/λ0=3, when δ=45°, φ=45°. (e) Longitudinal wave vector difference (kMZkEZ) between TM and TE components and (f) p + 1/p as a function of kx/k0 for different b values.

Download Full Size | PDF

The basic mechanism of the incident-angle-insensitivity of the IMAM polarization rotator can be understood by a simple argument. A lossless incident-angle-insensitive polarization rotator actually requires that TE(kx)TM(kx)exp[iϑ(kx)] can be satisfied for each single slab even at oblique incidence. Whenkx0, the transmission coefficient of TM wave through a homogeneous slab surrounded by air (the transmission coefficient of TE wave can be obtained by duality) can be written as

TM=2×[2cos(kMZd)i(pM+1/pM)sin(kMZd)]1,
where

pM=[(μyy/εxx)k02kx2/(εxxεzz)]1/2/(k02kx2)1/2.

The key issue is that for the impendence-mismatched slab, (μyy/εxx)1 and thuspM+pM1>2.Calculations show that (not shown here), the reflection increases quickly and become fluctuating whenkx>0;meanwhile, TE(kx)and TM(kx)change periodically with respect to kx, indicating the existence of high-order Fabry-Perot interference. While for the IMAM slab, (μyy/εxx)=1, pM1 (i.e. pM+1/pM2) can be maintained at a much larger kx, as depicted in Fig. 4(f), and thus TE(kx) and TM(kx) vary smoothly as a function of kx.Moreover, it can be seen from the dispersion curves in Fig. 4(e) that kMZkEZk0 is maintained better in the IMAM polarization rotator as kx/k0increases, which indicates a more stable phase difference. Hence, TE(kx)TM(kx) exp[iϑ(kx)] can be obtained, and Eq. (14) can be fulfilled even at a large-angle incidence. This is the reason that a bilayered polarization rotator composed of such impedance-matched half-wave retarders will thus exhibit incident-angle-insensitivity.

It should be noted that we mainly consider about the effects of transversal parameters (εxx,εyy,μxxandμyy) in the above discussion. Our calculations disclosed that the ratio between εzzand μzzwill also affect the property of angle-insensitivity. Numerical simulations show that the optimized ratio generally locates around one when the surroundings is air, i.e. εzzμzz.SincekEZ=(εyyμxxk02μxxkx2/μzz)1/2andkMZ=(μyyεxxk02εxxkx2/εzz)1/2, as the longitudinal parameters εzzand μzzincrease, the dispersion curve kx~kzof the metamaterial will become flatter. Thus, the IMAM polarization rotator will become more insensitive to incident angle and even exhibits a self-collimating property when εxx/εzzand μxx/μzz are sufficiently small [5]. Then, the ratio deviation betweenεzzand μzzcan be larger within the same tolerable variation range of Δ and ψt. To take the lossless IMAM polarization rotator discussed in Fig. 3 as an example, within the same tolerable variation range of Δ and ψt (i.e. |Δ| < 0.077 and ψt (90.69°,89.05°) at 40° incidence), when εzzμzz=5,the maximum ratio deviation can be up to 20% (i.e. εzz=6.0,μzz=5.0).

5. Construct a microwave IMAM rotator with bi-SRR structure

To demonstrate an incident-angle-insensitive and polarization-independent polarization rotator, we need to construct an IMAM half-wave retarder first, and then stack two retarders together with their optical axes rotated 45° to form such a polarization rotator. Various metamaterial structures can be employed to fabricate an IMAM retarder; however, metamaterials with both electric and magnetic resonances are preferred, since the effective permeability of nonmagnetic structures generally lies around one, which indicates that high refractive index and impedance-matching condition can hardly be achieved simultaneously. Moreover, spatial dispersion and anisotropy are inevitable because metamaterials are artificial mesostructures [23]. In order to diminish such effects, electrically small non-bianisotropic microwave bi-split-ring resonator (bi-SRR) is selected to construct an IMAM half-wave retarder [24].

Figures 5(a) and 5(b) are the schematic layouts of the IMAM retarder. The metallic bi-SRR patterns are fabricated on one side of the FR-4 (lossless) substrate characterized as ε=4.9ε0, μ=μ0,the substrate thickness t = 0.5 mm. The other dimensions of a unit cell are as follows: lattice constant az = ay = 5 mm; the length of metal slices (perfect electric conductor) in Z and Y directions z = y = 4 mm, the separation distance of the metal slices p = 0.12 mm, the gap g = 0.2 mm, the separation distance between adjacent unit cells s = 1.0 mm, and the width w and thickness of metal are 0.2 mm and 0.08mm respectively. The design principle of the IMAM half-wave retarder is that only TM component can excite magnetic resonances while TE wave propagates “quietly” through the structure. Thus nearly full transmission and low effective refractive index can be obtained for TE wave, while high effective index is available for TM component at the impedance-matched frequency near the resonance. By tuning the dimensions or the substrate material, one can change the impedance-matched frequency and the phase difference between the two orthogonal components of the transmitted waves.

 figure: Fig. 5

Fig. 5 (a) and (b) are schematic layouts of the bi-SRR half-wave retarder. az = ay = 5 mm, z = y = 4 mm, w = 0.2 mm, p = 0.12 mm, g = 0.2 mm, the substrate thickness t = 0.5 mm, the separation distance between adjacent layers s = 1.0 mm, and the thickness of metal (perfect electric conductor) is 0.08mm. (c) The structure of the polarization rotator, where yellow arrow indicates the incident waves, pink arrows indicate the direction of electric field.

Download Full Size | PDF

The simulation is performed using the software package CST Microwave Studio, in which periodic boundary conditions are applied. Figure 6(a) depicts the amplitudes of co-polarization terms of transmission coefficientT¯¯1(kx)and reflection coefficientR¯¯1(kx)for the first retarder, where the cross-polarization terms are negligible. There are two dips in the TM wave reflection coefficient curve corresponding to two possible impedance-matched frequencies at 3.12 GHz and 3.25 GHz. Here the working frequency is chosen to be 3.25GHz. Because the ratio of unit cell size to wavelength is less than 1/9, the effective medium description is valid, and the retrieval effective parameters at normal incidence are: nTE=1.156, nTM=3.457, ZTE=0.854 andZTM=1.001. Though the structure is not perfectly impedance-matched for TE wave at the working frequency, the transmission is still very high (>99%) and it will have less influence on the transmitted polarization state.

 figure: Fig. 6

Fig. 6 (a) The amplitudes of the transmission and reflection coefficients, (b) the phase of the transmission coefficient of TE and TM waves at normal incidence for polarization rotator with lossless substrate. (c) The transmission coefficients and (d) the phase difference between the two orthogonal polarized waves with respect to kx/k0 at 3.25 GHz.(e) and (g) Polarization rotation angle, (f) and (h) ellipticity variations with respect to incident wave polarization under different angles of incidence, where (e) and (f) are for lossless substrate withεr=4.9, (g) and (h) are for lossy substrate withεr=4.9+0.01.

Download Full Size | PDF

Note that in Fig. 6(b), the phase difference between two polarization states at 3.25GHz is 180°, which indicates that the structure acts as a half-wave retarder. Figure 6(c) reveals the amplitudes of transmission coefficients for TE and TM waves with respect to kx/k0 at 3.25GHz, and Fig. 6(d) shows the transmitted phase difference between the two orthogonal components, in which a 174° phase difference is still available at a 40° incidence (kx/k00.64). It is clear that high transmission and a near 180° phase difference can be maintained within a large incident angle.

Then, we construct a polarization rotator by aligning two such half-wave retarders in the aforementioned way, i.e. 2ϕ=90°, as depicted in Fig. 5(c). In the simulation, the polar angle of incident plane is rotated by 45° to perform an equivalent rotation of the IMAM retarder. Thus, we can get the transmission and reflection Jones matrices of the two IMAM retardersT¯¯1(kx),R¯¯1(kx)andT¯¯2(kx),R¯¯2(kx)directly from the simulation rather than from the effective parameters. The cross-polarization terms of the transmission coefficient for the second retarder T¯¯2(kx)satisfy TEM(kx)TME(kx) exp[iϑ(kx)] at 3.25GHz, as expected, and are not shown here. The total transmission Jones matrix of the rotator can then be written asT¯¯R(kx)=T¯¯2(kx)T¯¯1(kx)+O(T)(kx),whereO(T)=T¯¯2(kx)R¯¯1(kx)R¯¯2(kx)T¯¯1(kx)+, which is the small quantity due to multiple reflections between the slabs. The first term of the expression of O(T) represents the first order approximation. The ellipticity and polarization rotation angle with respect to incident wave polarization under different incident angles can then be worked out. For clarity but without loss of generality, we only give out the results for incident polar angle δ=0oin Figs. 6(e) and 6(f). The results under zeroth order and first order approximations are compared, it is clear that the variation ranges of ellipticity and polarization rotation can be well described under zeroth order approximation. The variations of ellipticity and polarization rotation are still quite acceptable even at the angle of incidenceθ=40°, showing that the structure can act as an incident-angle-insensitive and polarization-independent polarization rotator.

When a lossy substrate with a permittivity εr=4.9+εr'' is considered in the polarization rotator, it is found that the imaginary term may lead to a stronger absorption for TM wave near resonant frequencies. However, we can still find the two dips in reflection curve of the half-wave retarder. Although they are not as sharp as that in the lossless case, but they are small enough to suppress the influence of multiple reflections between the two retarders in a polarization rotator. For a commercial available low loss substrate withεr''0.01,our simulations show that the phase difference has a tiny change (less than 2°), and there is a small decrease in the transmission amplitude of TM wave due to absorption. The influences of the substrate loss are illustrated in Figs. 6(g) and 6(h) when εr''=0.01and the polarization rotator still works at 3.25GHz. Compared with the lossless rotator, the variations of ellipticity and polarization rotation angle are slightly increased. But we have to emphasize that the depicted results for the lossy rotator above are just for comparison, without any structural adjustment according to the loss, so the performance can be further improved after a proper structural optimization.

6. Conclusion

In conclusion, we proposed a design method to realize broad–angle and polarization-independent polarization rotator based on impedance-matched anisotropic metamaterial (IMAM). The mechanism and the influencing factors of the polarization rotator are discussed analytically, in conjunction with the generalized 4 × 4 transfer matrix method. Also, we illustrate how to construct a microwave IMAM rotator by using a bi-SRR structure. Compared with traditional polarization rotators, the IMAM polarization rotator is more insensitive to incident direction. This may offer the possibility to unyoke the limitation of the narrow working angle and enhance the compactness and stability of microwave and optoelectronic systems.

Acknowledgement

The authors acknowledge the financial support from the National Natural Science Foundation of China (NSFC) under the contracts 10774195, U0834001 and 10974263. The work is also partially supported by Program for New Century Excellent Talents in University and the Chinese National Key Basic Research Special Fund (2010CB923200).

References and links

1. D. Y. Yu, and H. Y. Tan, Engineering Optics(in chinese) (China Machine Press, Beijing, 2006).

2. J. Yao, Z. Liu, Y. Liu, Y. Wang, C. Sun, G. Bartal, A. M. Stacy, and X. Zhang, “Optical negative refraction in bulk metamaterials of nanowires,” Science 321(5891), 930 (2008). [CrossRef]   [PubMed]  

3. S. A. Ramakrishna, “Physics of negative refractive index materials,” Rep. Prog. Phys. 68(2), 449–521 (2005). [CrossRef]  

4. A. Salandrino and N. Engheta, “Far-field subdiffraction optical microscopy using metamaterial crystals: Theory and simulations,” Phys. Rev. B 74(7), 075103 (2006). [CrossRef]  

5. W. Zhang, J. Liu, W. P. Huang, and W. Zhao, “Self-collimating photonic-crystal wave plates,” Opt. Lett. 34(17), 2676–2678 (2009). [CrossRef]   [PubMed]  

6. J. Zhao, Y. Feng, B. Zhu, and T. Jiang, “Sub-wavelength image manipulating through compensated anisotropic metamaterial prisms,” Opt. Express 16(22), 18057–18066 (2008). [CrossRef]   [PubMed]  

7. J. Valentine, S. Zhang, T. Zentgraf, E. Ulin-Avila, D. A. Genov, G. Bartal, and X. Zhang, “Three-dimensional optical metamaterial with a negative refractive index,” Nature 455(7211), 376–379 (2008). [CrossRef]   [PubMed]  

8. J.-M. Lourtioz, “Photonic crystals and metamaterials,” C. R. Phys. 9(1), 4–15 (2008). [CrossRef]  

9. M. Beruete, M. Navarro-Cía, M. Sorolla, and I. Campillo, “Polarization selection with stacked hole array metamaterial,” J. Appl. Phys. 103(5), 053102 (2008). [CrossRef]  

10. J. Zhao, Y. Chen, and Y. Feng, “Polarization beam splitting through an anisotropic metamaterial slab realized by a layered metal-dielectric structure,” Appl. Phys. Lett. 92(7), 071114 (2008). [CrossRef]  

11. H. Luo, Z. Ren, W. Shu, and F. Li, “Construct a polarizing beam splitter by an anisotropic metamaterial slab,” Appl. Phys. B 87(2), 283–287 (2007). [CrossRef]  

12. V. Zabelin, L. A. Dunbar, N. Le Thomas, R. Houdré, M. V. Kotlyar, L. O’Faolain, and T. F. Krauss, “Self-collimating photonic crystal polarization beam splitter,” Opt. Lett. 32(5), 530–532 (2007). [CrossRef]   [PubMed]  

13. J. Hao, Y. Yuan, L. Ran, T. Jiang, J. A. Kong, C. T. Chan, and L. Zhou, “Manipulating electromagnetic wave polarizations by anisotropic metamaterials,” Phys. Rev. Lett. 99(6), 063908 (2007). [CrossRef]   [PubMed]  

14. K. Bayat, S. K. Chaudhuri, and S. Safavi-Naeini, “Ultra-compact photonic crystal based polarization rotator,” Opt. Lett. 17, 7145–7158 (2009).

15. J. Y. Chin, J. N. Gollub, J. J. Mock, R. Liu, C. Harrison, D. R. Smith, and T. J. Cui, “An efficient broadband metamaterial wave retarder,” Opt. Express 17(9), 7640–7647 (2009). [CrossRef]   [PubMed]  

16. T. Li, H. Liu, S. M. Wang, X. G. Yin, F. M. Wang, S. N. Zhu, and X. Zhang, “Manipulating optical rotation in extraordinary transmission by hybrid plasmonic excitations,” Appl. Phys. Lett. 93(2), 021110 (2008). [CrossRef]  

17. M. Kuwata-Gonokami, N. Saito, Y. Ino, M. Kauranen, K. Jefimovs, T. Vallius, J. Turunen, and Y. Svirko, “Giant optical activity in quasi-two-dimensional planar nanostructures,” Phys. Rev. Lett. 95(22), 227401 (2005). [CrossRef]   [PubMed]  

18. S. K. Awasthi and S. P. Ojha, “Wide-angle, broadband plate polarizer with 1D photonic crystal,” Prog. Electromag. Res. PIER 88, 321–335 (2008). [CrossRef]  

19. Y. Avitzour, Y. A. Urzhumov, and G. Shvets, “Wide-angle infrared absorber based on a negative-index plasmonic metamaterial,” Phys. Rev. B 79(4), 045131 (2009). [CrossRef]  

20. J. L. Tsalamengas, “Interaction of electromagnetic waves with general bianisotropicslabs,” IEEE Trans. Microwave Theory Tech. 40(10), 1870–1878 (1992). [CrossRef]  

21. R. M. A. Azzam, and N. M. Bashara, Ellipsometry and Polarized Light (North-Holland Pub. Co., New York, 1977).

22. B. Bai, Y. Svirko, J. Turunen, and T. Vallius, “Optical activity in planar chiral metamaterials: Theoretical study,” Phys. Rev. A 76(2), 023811 (2007). [CrossRef]  

23. C. Menzel, C. Rockstuhl, T. Paul, F. Lederer, and T. Pertsch, “Retrieving effective parameters for metamaterials at oblique incidence,” Phys. Rev. B 77(19), 195328 (2008). [CrossRef]  

24. J. D. Baena, L. Jelinek, R. Marques, and J. Zehentner, “Electrically small isotropic three-dimensional magnetic resonators for metamaterial design,” Appl. Phys. Lett. 88(13), 134108 (2006). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 Schematic structure of a multilayered metamaterial rotator. The plane of incidence is in x-z plane. ϕ is the angle between the axes of the adjacent layers. φ and θ are the polarization angle and incident angle of the incident wave, respectively.
Fig. 2
Fig. 2 Schematic view of the polarization vector projections on x-y plane when the light is normally incident. The crystal axes of the anisotropic slabs are rotated along z direction by an angleϕin sequence.
Fig. 3
Fig. 3 The projections of (a) the ellipticity Δ and (b) the polarization rotation angle ψ t with respect to k x / k 0 and δ when φ = 45°. The projection views of (c) Δ and (d) ψ t with respect to δ and φ at 40° incident angle, in which φ ranges from −90° to 90°. The insets demonstrate the three-dimensional relations. The lossless IMAM bilayered polarization rotator is characterized with ε x x = μ y y = 2 , ε y y = μ x x = 1 , ε z z = μ z z = 1 , and a total thickness 2 h / λ 0 = 3.
Fig. 4
Fig. 4 (a) The amplitudes of transmission coefficients | T EM | and | T M E | , (b) the phase difference T E M T M E , (c) the polarization rotation angle ψ t and (d) the ellipticity Δ with respect to the incident transversal vector k x / k 0 for three different anisotropic bilayered structures characterized as ε x x = 2 b , ε y y = ε z z = 1 , μ y y = 2 / b , μ x x = μ z z = 1 , and the total thickness 2 h / λ 0 = 3 , when δ = 45 ° , φ = 45 ° . (e) Longitudinal wave vector difference ( k MZ k EZ ) between TM and TE components and (f) p + 1/p as a function of k x / k 0 for different b values.
Fig. 5
Fig. 5 (a) and (b) are schematic layouts of the bi-SRR half-wave retarder. az = ay = 5 mm, z = y = 4 mm, w = 0.2 mm, p = 0.12 mm, g = 0.2 mm, the substrate thickness t = 0.5 mm, the separation distance between adjacent layers s = 1.0 mm, and the thickness of metal (perfect electric conductor) is 0.08mm. (c) The structure of the polarization rotator, where yellow arrow indicates the incident waves, pink arrows indicate the direction of electric field.
Fig. 6
Fig. 6 (a) The amplitudes of the transmission and reflection coefficients, (b) the phase of the transmission coefficient of TE and TM waves at normal incidence for polarization rotator with lossless substrate. (c) The transmission coefficients and (d) the phase difference between the two orthogonal polarized waves with respect to k x / k 0 at 3.25 GHz.(e) and (g) Polarization rotation angle, (f) and (h) ellipticity variations with respect to incident wave polarization under different angles of incidence, where (e) and (f) are for lossless substrate with ε r = 4.9 , (g) and (h) are for lossy substrate with ε r = 4.9 + 0.01.

Equations (18)

Equations on this page are rendered with MathJax. Learn more.

= ( cos ϕ sin ϕ 0 sin ϕ cos ϕ 0 0 0 1 ) ,
P ¯ ¯ t o t a l ( H ) = P ¯ ¯ m ( h m ) P ¯ ¯ m 1 ( h m 1 ) P ¯ ¯ 1 ( h 1 ) ,
E t ( sin φ cos φ ) = E 0 ( T M M T M E T E M T E E ) ( cos φ sin φ ) ,
0 = E 0 [ T M M cos 2 φ + T E E sin 2 φ + ( T M E + T E M ) sin φ cos φ ] .
ε r = ( ε x x 0 0 0 ε y y 0 0 0 ε z z ) ,
μ r = ( μ x x 0 0 0 μ y y 0 0 0 μ z z ) ,
ε y y / μ x x = ε x x / μ y y = ε s u r / μ s u r = 1 / Z .
T E = exp ( i ε y y ' k 0 h Z ) exp ( ε y y ' ' k 0 h Z ) ,
T M = exp ( i μ y y ' k 0 h / Z ) exp ( μ y y ' ' k 0 h / Z ) ,
ε y y ' ' = ε x x ' ' = μ x x ' ' / Z 2 = μ y y ' ' / Z 2 ,
( E 1 < n + 1 > E 2 < n + 1 > ) = ( T E cos ϕ T E sin ϕ T M sin ϕ T M cos ϕ ) n × { ( T E cos δ T E sin δ T M sin δ T M cos δ ) ( E i M E i E ) } .
( T M M T M E T E M T E E ) = ( cos ( n ϕ + δ ) sin(n ϕ + δ ) sin ( n ϕ + δ ) cos ( n ϕ + δ ) ) × { ( T E cos ϕ T E sin ϕ T M sin ϕ T M cos ϕ ) n ( T E cos δ T E sin δ T M sin δ T M cos δ ) } .
( T M M T M E T E M T E E ) = T n + 1 ( cos [ ( n + 1 ) ϕ ] sin [ ( n + 1 ) ϕ ] sin [ ( n + 1 ) ϕ ] cos [ ( n + 1 ) ϕ ] ) .
( T M M T M E T E M T E E ) = T n + 1 ( 0 ± 1 1 0 ) .
Δ = E a / E b ,
ψ t = φ t φ ,
T M = 2 × [ 2 cos ( k M Z d ) i ( p M + 1 / p M ) sin ( k M Z d ) ] 1 ,
p M = [ ( μ yy / ε x x ) k 0 2 k x 2 / ( ε x x ε zz ) ] 1 / 2 / ( k 0 2 k x 2 ) 1 / 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.