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

Polarity-variable birefringence on hyperlens structure

Open Access Open Access

Abstract

A variable birefringence effect has been observed with 1D PMMA surface gratings on a gold film substrate. By changing the operation wavelength on the Au film, the birefringence value Δneff changes from positive to negative. The result verified that this uniaxial crystal-like plasmonic surface gratings showed good superlensing effect at 515nm when PMMA width:Air width=1:1 where the absolute value of Δneff is relatively large.

©2010 Optical Society of America

1. Introduction

The optical superlensing effect has been caused attentions since J. Pendry reports the negative lens concept [1]. Various kinds of superlens have been discussed [25]. It is interesting to note that the two dimensional hyperlens structure proposed by Narimanov et.al. in 2006 has hyperbolic dispersion characteristic [2]. This hyperlens can carry very large k-vector and achieve super-resolution effect by tolerating the infinite propagation wave-vectors. Except the imaging function, the proposed hyperlens structure acts like a uniaxial plasmonic crystal. A more detailed explanation can be referred to [610]. The Dyakonov states on uniaxial crystal film [68] were adopted to explain the super-resolution effect demonstrated by Smolyaninov [4] more specifically. The 1D grating structures inherit birefringence property and surface plasmon polariton (SPP) experiences different in-plane dielectric constants. Here we measured the birefringence Δneff of the 1D straitened hyperlens structure as shown in Fig. 1 .(a). A similar analysis on anomalous phase from periodic birefringent material without SPP has been discussed in [11,12].

 figure: Fig. 1

Fig. 1 (a) Structure of the sample with Au thickness 50nm and PMMA thickness 150nm. (b) Grating pitch is 500nm. Polarization of the incident light is rotated along x-y plane where z is the direction of light incidence.

Download Full Size | PDF

The interesting point of our birefringence observation is that Δneff is usually either positive or negative once crystal orientation is fixed. Here we report Δneff can change from positive to negative simply by changing the operation wavelength without changing other parameters like periodicity and grating orientation. If this crystal is made from Poly(methyl methacrylate) (PMMA) material, there should be no any chromatic dispersion which causes Δneff flip sign unless SPP waves at the interface between the Au film and the dielectric gratings play a role inside. An elliptical index sphere is used to describe uniaxial crystal material and it is naturally to believe that it should have fixed long and short axes once the grating orientation is determined. The expression of the dielectric constants on a finely dielectric layered structure can be found in [13]. It shows that the roles of long axis and short axis never swap when the dielectric constants of both layers are positive.

2. Experiment

PMMA stripe gratings with 150nm thickness have been patterned on a 50nm thick gold film on a cover glass substrate. The Au film is intact under the PMMA gratings which periodicity is 500nm as shown in Fig. 1(a). 515nm Ar laser light is normally incident from the substrate side with polarization angle θ varied from 0° to 180°. TM and TE modes are defined as E-field along x and y axes respectively as shown in Fig. 1(b).

Some different PMMA stripe widths with same periodicity have been fabricated and shown in Fig. 2 . The optical measurement setup is shown in Fig. 3 . The open aperture (A) is around 2mm*2mm placed behind the objective lens. The magnified image has been projected onto the photo-detector and no other scattering light from non-grating area was collected. In order to get the four Stokes parameters, we use a quarter wave plate, half wave plate and polarizers to get the following operations:

S0=I0+I90
S1=I0-I90
S2=I45-I135
S3=IRHC-ILHC
(S1, S2, S3) is the (x,y,z) position of the transmitted signal on the Poincare sphere. RHC and LHC present right-hand circular and left-hand circular polarizations, respectively. Ix is the optical intensity of polarization at x-degree. The polarization state of a light can be determined after the three orthogonal indexes (S1, S2, S3) being determined.

 figure: Fig. 2

Fig. 2 PMMA width: Air width (a) 1:1 (b) 1:2 (c) 1:3.

Download Full Size | PDF

 figure: Fig. 3

Fig. 3 Optical setup for Poincare sphere measurement. λ/2: half wave-plate, P1: polarizer, Obj: 20X objective, λ/4: quarter wave-plate, P2: analyzer, D: detector.

Download Full Size | PDF

Referring to Fig. 1(b), suppose the intrinsic polarization of the laser output is horizontal (θ = 0°), which is expressed as θ(Εlaser) = 0°. If we want to check what the response of the sample is when the polarization of the incident light Ein at the sample surface is rotated to θ = 15° (which is expressed as θ(Εin) = 15°), the optical elements before the sample should be arranged as θ(Εlaser) = 0°, θ(λ/2−fast axis) = 7.5°, and θ(P1 transmission axis) = 15°. This gives us θ(Εin) = 15°. The λ/2 wave plate is like a mirror which flips the polarization of the laser output from 0° to 15° without losing optical power. P1 is used to guarantee that the polarization of the light is the same as what we expected. The rotational angles of the analysis elements after the sample are arranged as shown in Table 1 . IRHC is measured when the analyzer P2 is at θ(P2) = θ(Ein) + 135° while ILHC is acquired when P2 is at θ(P2) = θ(Ein) + 45°. When the relative angle between the fast axis of λ/4 wave plate and P2 is fixed at 135°, this optical set λ/4 -P2 should be viewed as a part which can completely block out LHC polarization light. In other words, when λ/4 -P2 has 45° difference, the optical set is orthogonal to RHC polarization. The detailed polarization measurement can be found in [14] or any textbooks related to Jones vectors.

Tables Icon

Table 1. The orientations of the analysis components λ/4 and P2 when θ(Εin) = 15°.

In the measurement setup and above descriptions, we actually fix the sample and rotate the coordinates of the input and output optical components simultaneously in counterclockwise direction. It is same as the sample being rotated 15° clockwise. In fact, the data points are taken with every 15° rotation of the optical elements from 0° to 180°.

3. Results and discussions

The degree of polarization (DOP) is (S1+S2+S3)/S0. In our measurement, DOP ranges from 0.98 to 1.0 which shows that the sample is a linear crystal. If a crystal is isotropic, its Poincare path will just stick to the same point and does not move due to no birefringence. In our case, if the sample is air, the Poincare path will just stay at (1,0,0) after we took the data from 0° to 180°. To understand further what the Poincare sphere represents, we first look at the Stokes parameters of an elliptical polarized light. The elliptical light has E field components on x and y axis:

Ex=Ex0cos(ωtkzδx)
Ey=Ey0cos(ωtkzδy)
The definitions of the four Stokes parameters are:

S0=(Ex0)2+(Ey0)2
S1=(Ex0)2-(Ey0)2
S2=2(Ex0)(Ey0)cos(δy-δx)
S3=2(Ex0)(Ey0)sin(δy-δx)

The handness of the polarization is determined by the sign of phase difference δyx between Ey and Ex. The ellipticity is determined by the ratio of Ex to Ey. By calculating the four Stokes parameters, the polarization state of the light can be determined and marked on the Poincare sphere surface as shown in Fig. 4(a) . Each position represents different handness, phase difference, and ellipticity. To illustrate further, Fig. 4(b) shows the Poincare paths for λ/8, λ/4, and λ/2 wave plates by rotating them every 15° from 0° to 90°. If the wave plate provides enough phase difference, the locus of the light can circle around the equator; otherwise, it will just go around a loop. If it is λ/4 wave plate, the light can reach RHC polarization point at θ =45°.

 figure: Fig. 4

Fig. 4 (a) The corresponding polarization states on the Poincare sphere (b) The Poincare paths for λ/8, λ/4, and λ/2 wave plates by rotating them every 15° from 0° to 90°.

Download Full Size | PDF

The measurement on polarization state of the transmission light through the PMMA surface gratings has been depicted in Fig. 5 . It is interesting to note that δyx(=Δneff*ko*d) flips signs when λ changes from 515nm to 633nm on Au sample. In our case, δyx > 0 means right-hand elliptical polarization. Four conditions in Fig. 5 all show birefringence responses. It is interesting to note that the birefringence of the PMMA gratings on Au substrate is more eminent than that on ITO substrates. There is still some birefringence effect (with same polarity) on ITO substrate as shown in Fig. 5 (c) and (d). Light at shorter wavelength shows stronger Δneff which can be explained by the band behavior of 1D periodic layered structure [15].

 figure: Fig. 5

Fig. 5 Poincare sphere plot of PMMA surface gratings on different substrates with different width ratio (a) Au film substrate, 515nm, (b) Au substrate, 633nm(c) ITO substrate, 515nm, (d) ITO substrate, 633nm.

Download Full Size | PDF

To analyze the system in an easier way, we fix the coordinates of the optical elements but rotate the sample gratings in clockwise direction. Assume the matrix of our PMMA surface gratings is ABCD matrix. The matrix expressions for our system setup are as follows:

InputP1=[10]
SamplematrixW(θ)=[cos(θ)sin(θ)-sin(θ)cos(θ)][ABCD][cos(θ)-sin(θ)sin(θ)cos(θ)]
AnalyzerP2(ψ)=[cos(ψ)sin(ψ)-sin(ψ)cos(ψ)][1000][cos(ψ)-sin(ψ)sin(ψ)cos(ψ)]=[cos2(ψ)-sin(ψ)cos(ψ)-sin(ψ)cos(ψ)sin2(ψ)]
QWP=[ejπ/4o0e-jπ/4]fastaxisonx^axiswithej(ωtkz)convention
θ and ψ are the rotational angles of the sample and the analyzer P2 with respect to the polarizer P1. For example, if the sample is rotated θ degrees clockwise, the intensity of the right-hand circular polarization light can be measured as:

Analyzer P2(135°)*QWP*Sample (-θ) *Input P1(0°), i.e.

IRHC=[cos2(135)-sin(135)cos(135)-sin(135)cos(135)sin2(135)][ejπ/4o0e-jπ/4]*W(-θ)*[10]

For simplicity, Fig. 5 just shows the data with the sample rotated from 0° to 90°. Table 2 shows the fitting parameters of the ABCD sample matrix. The mean deviation of the 7 data points is expressed as: i=1i=7[(S1S1')2+(S2S2')2+(S3S3')2]1/2/7. For better understanding, the off-diagonal coefficients are set zero. Except the 1:1 case on Au substrate, others show little birefringence. At 515nm, the amplitudes of A and D coefficients on 1:1 Au substrate are different. It means that except the phase difference along x and y axis, the scattering strength is also different on both axes.

Tables Icon

Table 2. Fitting data of ABCD matrix of the sample element.

From [6], when the operation wavelength is at 633nm, the dispersion curve of the Dyakonov SPP waves is a closed ellipse, just like a normal dispersion curve on a layer of uniaxial crystal film. But the mode fields are Dyakonov SPP waves with wave vectors defined as q|| andq. In Fig. 1 (b), q||means the wave vector is along x-axis. The ellipse has specific long and short axes. The polarization states of the transmission light are like the ensemble of the Dyakonov SPP waves excited on the grating structure. These SPP waves get scattered out in z-direction and the output phase we measured in z-direction is the combination of in-plane waves Ex and Ey with SPP wave vectors(q||,q). On average, (q||q)is almost zero if the dispersion curve is a closed circle. Figure 6(a) shows that if the phase difference is small, the closed dispersion curve in center is almost a circle. We assume the incident light excites (q||,q) in the first quadrant in Fig. 6 due to no perfect illumination on the sample. With the aid of Au film, the anisotropy has been enhanced as shown in Fig. 5 (a) and (b). According to [6], the anisotropy at 633nm shown in Fig. 5(b) should come from the existence of both hyperbolic dispersion and circular dispersion curves as is shown in Fig. 6(a). In Fig. 5(a) the anisotropy comes from the hyperbolic dispersion curve only and it is shown in Fig. 6(b). Our measurement in Fig. 5(a) and (b) show that the roles of the long and short axes of the ellipse swap when wavelength is tuned from 515nm to 633nm. That means Δneff*ko*d changed sign or the handness of the sample changed due to wavelength change. It also means that the sign of the total integral of(q||q)at 515nm is opposite to that at 633nm.

 figure: Fig. 6

Fig. 6 The dispersion curves of the 1D PMMA gratings on Au film according to our experimental data. (a) Δneff*ko*d ~-2π/16 at 633nm. (b) Δneff*ko*d ~2π/7 at 515nm.

Download Full Size | PDF

It can be deducted that the anisotropy of the hyperlens system demonstrated in [4] is much more effective at 515nm as well as the cone of the energy flow is much well-defined compared to the device at 633nm. In [4], the superresolution is achieved when the energy propagates along x-axis. With different width ratio between air gap and PMMA stripe, Δneff changes and large Δneff can be achieved at 1:1 in our demonstration. In fact, maximum Δneff can be determined by judging from the band diagrams of this system. If the operation wavelength is tuned, slight change in Δneff will result in different ellipse in Fig. 6(b). To get maximum ellipticity, the optimization on the width ratio might be needed. According to our measurement, 1:1 width ratio is good for imaging function in [4]. With larger index contrast between stripes and air, the practical superresolution limit can be pushed further.

4. Conclusion

A variable optical birefringence effect on surface gratings has been observed. The refractive index difference Δneff of the uniaxial crystal-like grating covers positive and negative regions by changing the operation wavelength. The superresolution effect in [4] and the Dyakonov states in [6] have been discussed. The measured birefringence suggests that when PMMA width: Air width =1:1, we can get large birefringence from our sample. The resolution capability can be improved.

Acknowledgments

We thank Professor Yi-Jen Chiu for the assistance on Au film deposition, Professor Tsung-Hsien Lin for polarimetry measurement discussion. Both are from Dept. of Photonics, National Sun Yat-sen University. We thank Professor Christopher C. Davis from University of Maryland College Park for suggestions on polarimetry measurement. This work is supported by National Science Council of the Republic of China (Taiwan), (Contract No. NSC 99-2112-M-100-001-MY3)

References and links

1. J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. 85(18), 3966–3969 (2000). [CrossRef]   [PubMed]  

2. Z. Jacob, L. V. Alekseyev, and E. E. Narimanov, “Optical Hyperlens: far-field imaging beyond the diffraction limit,” Opt. Express 14(18), 8247 (2006). [CrossRef]   [PubMed]  

3. Z. Liu, H. Lee, Y. Xiong, C. Sun, and X. Zhang, “Far-field optical hyperlens magnifying sub-diffraction-limited objects,” Science 315(5819), 1686 (2007). [CrossRef]   [PubMed]  

4. I. I. Smolyaninov, Y. J. Hung, and C. C. Davis, “Magnifying superlens in the visible frequency range,” Science 315(5819), 1699–1701 (2007). [CrossRef]   [PubMed]  

5. W. Cai, D. A. Genov, and V. M. Shalaev, “Superlens based on metal-dielectric composites,” Phys. Rev. B 72(19), 193101 (2005). [CrossRef]  

6. Z. Jacob and E. E. Narimanov, “Optical hyperspace for plasmons: Dyakonov states in metamaterials,” Appl. Phys. Lett. 93(22), 221109 (2008). [CrossRef]  

7. R. Li, C. Cheng, F. F. Ren, J. Chen, Y. X. Fan, J. Ding, and H. T. Wang, “Hybridized surface plasmon polaritons at an interface between a metal and a uniaxial crystal,” Appl. Phys. Lett. 92(14), 141115 (2008). [CrossRef]  

8. M. I. D’yakonov, “New type of electromagnetic wave propagating at an interface,” Sov. Phys. JETP 67, 714 (1988).

9. A. Drezet, C. Genet, and T. W. Ebbesen, “Miniature plasmonic wave plates,” Phys. Rev. Lett. 101(4), 043902 (2008). [CrossRef]   [PubMed]  

10. L. Feng, Z. Liu, V. Lomakin, and Y. Fainman, “Form birefringence metal and its plasmonic anisotropy,” Appl. Phys. Lett. 96(4), 041112 (2010). [CrossRef]  

11. A. Mandatori, C. Sibilia, M. Bertolotti, S. Zhukovsky, J. W. Haus, and M. Scalora, “Anomalous phase in one-dimensional, multilayer, periodic structures with birefringent materials,” Phys. Rev. B 70(16), 165107 (2004). [CrossRef]  

12. S. Foteinopoulou and C. M. Soukoulis, “Electromagnetic wave propagation in two-dimensional photonic crystals:A study of anomalous refractive effects,” Phys. Rev. B 72(16), 165112 (2005). [CrossRef]  

13. S. M. Tytov, “Electromagnetic Properties of a Finely Stratified Medium,” Sov. Phys. JETP 2, 466 (1956).

14. P. Yeh, and C. Gu, Optics of Liquid Crystal Displays (Wiley, 1999), pp. 22–41.

15. P. Yeh, Optical waves in layered media (Wiley, 1991), pp. 118–138.

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 (a) Structure of the sample with Au thickness 50nm and PMMA thickness 150nm. (b) Grating pitch is 500nm. Polarization of the incident light is rotated along x-y plane where z is the direction of light incidence.
Fig. 2
Fig. 2 PMMA width: Air width (a) 1:1 (b) 1:2 (c) 1:3.
Fig. 3
Fig. 3 Optical setup for Poincare sphere measurement. λ/2: half wave-plate, P1: polarizer, Obj: 20X objective, λ/4: quarter wave-plate, P2: analyzer, D: detector.
Fig. 4
Fig. 4 (a) The corresponding polarization states on the Poincare sphere (b) The Poincare paths for λ/8, λ/4, and λ/2 wave plates by rotating them every 15° from 0° to 90°.
Fig. 5
Fig. 5 Poincare sphere plot of PMMA surface gratings on different substrates with different width ratio (a) Au film substrate, 515nm, (b) Au substrate, 633nm(c) ITO substrate, 515nm, (d) ITO substrate, 633nm.
Fig. 6
Fig. 6 The dispersion curves of the 1D PMMA gratings on Au film according to our experimental data. (a) Δneff*ko*d ~-2π/16 at 633nm. (b) Δneff*ko*d ~2π/7 at 515nm.

Tables (2)

Tables Icon

Table 1 The orientations of the analysis components λ/4 and P2 when θ(Εin) = 15°.

Tables Icon

Table 2 Fitting data of ABCD matrix of the sample element.

Equations (15)

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

S 0 = I 0 + I 90
S 1 = I 0 - I 90
S 2 = I 45 - I 135
S 3 = I R H C - I L H C
E x = E x 0 cos ( ω t k z δ x )
E y = E y 0 cos ( ω t k z δ y )
S 0 = ( E x 0 ) 2 + ( E y 0 ) 2
S 1 = ( E x 0 ) 2 - ( E y 0 ) 2
S 2 = 2 ( E x 0 ) ( E y 0 ) cos ( δ y - δ x )
S 3 = 2 ( E x 0 ) ( E y 0 ) sin ( δ y - δ x )
Input P 1 = [ 1 0 ]
Sample matrix W(θ) = [ cos (θ) sin (θ) - sin (θ) cos (θ) ] [ A B C D ] [ cos (θ) - sin (θ) sin (θ) cos (θ) ]
Analyzer P 2 (ψ) = [ cos (ψ) sin (ψ) - sin (ψ) cos (ψ) ] [ 1 0 0 0 ] [ cos (ψ) - sin (ψ) sin (ψ) cos (ψ) ] = [ cos 2 (ψ) - sin (ψ) cos (ψ) - sin (ψ) cos (ψ) sin 2 (ψ) ]
QWP = [ e jπ/4 o 0 e -jπ/4 ] fast axis on x ^ axis with e j ( ω t k z ) convention
I R H C = [ cos 2 ( 135 ) - sin ( 135 ) cos ( 135 ) - sin ( 135 ) cos ( 135 ) sin 2 ( 135 ) ] [ e j π / 4 o 0 e - j π / 4 ] * W ( - θ ) * [ 1 0 ]
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.