International Journal of PhotoenergyVolume 2012 (2012), Article ID 702637, 7 pagesdoi:10.1155/2012/702637

Research Article

Negative Propagation Effects in Two-Dimensional Silicon Photonic Crystals

Ping Jiang,1,2 Kang Xie,² Huajun Yang,¹ and Zhenhai Wu²

¹School of Physical Electronics, University of Electronic Science and Technology of China, Chengdu 610054, China²School of Optoelectronic Information, University of Electronic Science and Technology of China, Chengdu 610054, China

Received 31 May 2011; Revised 29 August 2011; Accepted 12 September 2011

Academic Editor: Leonid Khriachtchev

Copyright © 2012 Ping Jiang et al. Ping Jiang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We demonstrated negative refraction effects of light propagating in two-dimensional square and hexagonal-lattice silicon photonic crystals (PhCs). The plane wave expansion method was used to solve the complex eigenvalue problems, as well as to find dispersion curves and equal-frequency contour (EFC). The finite-difference time-domain (FDTD) method was used to simulate and visualize electromagnetic wave propagation and scattering in the PhCs. Theoretical analyses and numerical simulations are presented. Two different kinds of negative refractions, namely, all-angle negative refraction (AANR) without a negative index and negative refraction with effective negative index, have been verified and compared.

1. Introduction

At the beginning of the 21st century, a new type of artificial material that displays anomalous properties in their interactions with electromagnetic waves has attracted a significant degree of scientific interests. These media are commonly known as “negative index materials” (NIMs) (or left-handed materials, LHMs), since its permittivities and permeabilities are simultaneously negative. NIMs can be exploited for use in a wide variety of guiding and radiating structures, such as absorber, filter, coupler, antenna, superlens, and superprism. The direction of phase velocity of light wave propagating inside such a material is opposite to the direction of its energy flow. It means that the Poynting vector and wave vector are antiparallel; therefore negative refraction occurs at the interface [1]. Veselago first pointed out the negative refraction phenomena in 1968 [2]. Years later, many theoretical and experimental groups have investigated left-handed materials [3–5].

Since PhCs are more feasible to fabricate in optical region than negative index materials [6], effective-negative-index PhCs attracted tremendous interest of electrical engineers and physicists. One of the most prominent properties of negative index PhC is its capability of amplification of evanescent waves. Light propagating in such PhCs may exhibit negative refraction near the photonic bandgaps, such as mirror-like effect, image-transfer effect, and superprisming effect [7–12]. The effective refraction index of such PhC, which is analogous to the effective mass in solid state physics [13], is controllable by the band structure of the material. Notomi studied negative refraction in strongly modulated PhCs near the photonic bandgaps [14]. Moussa et al. investigated superlensing behavior in 2D PhC [15]. Cubukcu et al. demonstrated experimentally single beam negative refraction and superlensing in the valence band of 2D PhCs operating in the microwave regime [16]. Berrier et al. presented an experimental verification of negative refraction for infrared wavelength in an effective-negative-index PhC [17]. Zhang proposed a 2D metal-core PhC, which can achieve effective negative indexes for both TE and TM modes [18]. Lu et al. demonstrated an imaging experiment by using a slab of an effective-negative-index PhC in microwave region [19].

It was demonstrated that negative refraction phenomena in PhCs are possible in regimes of negative group velocity and negative effective index above the first band near the Brillouin zone center Γ [20]. Negative refraction can also be resulted from some special EFC in the first band of PhC [21]. The frequencies in the lowest band may be more desirable in high-resolution superlensing.

In this paper we study two different kinds of negative refraction at the first band and the second band, respectively. Firstly we study AANR without a negative index, which is resulted from a special EFC near the Brillouin zone corner M in the first band of a square-lattice silicon PhC. Secondly, the negative refraction with effective negative index is investigated, which happens in regime of negative group velocity and negative effective index near the Brillouinzone center Γ in the second band of a hexagonal-lattice silicon PhC. The plane wave expansion method is used to solve complex eigenvalue problems. The photonic band structures and equifrequency surfaces are calculated, and two different negative refraction frequency ranges are found in two different types of PhC lattices. The finite-difference time-domain (FDTD) method is employed to simulate and visualize the subwavelength imaging effects.

2. All-Angle Negative Refraction without a Negative Index

In this part we consider an AANR effect that does not employ a negative effective index of refraction in the PhCs. In the first band of a square-lattice PhC, the EFC is not necessarily a circle around Γ point, but a convexity around M point for certain frequency (as shown in Figure 2). So no effective refractive index can be defined in such PhCs. This kind of PhCs is referred to as AANR PhCs.

We begin with TE modes (in-plane electric field vector) and consider a square lattice of air holes in silicon host material with . The hole radius , where a is the lattice constant. Firstly, we employ the plane wave expansion method to calculate the AANR frequency range in our model system. The band curves for TE polarization is shown in Figure 1. It is identified that the AANR region occurs in the vicinity of . The upper frequency limit is determined directly by the intersection of the first band curve with the light line. The lower frequency limit is determined by the frequency at which radius of the curvature of the contours along ΓM becomes divergent (as shown in Figure 2(b)).

Figure 1: Band curves of square lattice for TE polarization. The 2D PhC is made of air holes of radius (where is the lattice constant) drilled in a silicon host substrate with . The light line shifted to M point is shown in black. The red line indicates the frequency range of negative refraction.

(a)

(b)

Figure 2: (a) Wave vector diagram of negative refraction in the PhC. (b) EFCs of the first band of the model PhC in AANR. Value of frequency is in unit of .

The contours of dispersion of the first band, drawn in the first Brillouin zone, are shown in Figure 2(a), where the thick arrows represent directions of group velocity, and thin arrows stand for directions of phase velocity. The EFCs of the actual PhC are shown in Figure 2(b). AANR occurs at frequencies that correspond to all-convex contours with a negative photonic effective mass. If the air-PhC interface normal is along the Γ-M direction, and the contour is everywhere convex, an incoming plane wave from air will couple to a single mode that propagates in this PhC on the negative side of the boundary normal.

Figure 3 shows the negative refraction of a light beam across a square-lattice 2D PhC slab. The incoming wave with normalized frequency of is incident from up right side with an angle of 45° with respect to the normal (as shown by the arrows). The refracted light inside the slab is on the right side of the normal. On its arrival at the lower surface of the slab, negative refraction happens again. The outgoing beam has the same propagation direction as the incoming beam but with a negative transverse shift.

Figure 3: The negative refraction phenomena of PhC.

AANR is essential for superlensing applications. The existence of AANR effect implies that there is essentially no physical limit on aperture of the imaging system. Such a superlens possesses several key advantages over conventional lenses. To realize AANR for superlensing, the required conditions for our model system are that the EFCs should be both convex and larger than the EFCs for air (circles with radius ). Incident beams at any incident angle will then experience negative refraction when entering the PhC. Thickness and finishing of the surfaces of the PhC slab should be chosen carefully to minimize reflections. Moreover, when air holes (or dielectric rods) are cut in half at the terminating surfaces, surface mode can come to exist on the slab. The transmission of evanescent waves can be amplified through the excitation of surface modes in a slab of PhC, which plays a very important role in subwavelength imaging.

Shown in Figure 4 are the FDTD simulations of a continuous-wave point source in front of a PhC slab of square lattice. The radiation is TE polarized at frequency. The slabs in Figures 4(a)–4(c) are, respectively, 8 layers, 12 layers, and 20 layers thick. The distances of point sources to the left-hand side surfaces of the slabs are respectively, , , and . It is seen that an image of the point source is formed in each case at, respectively, a distance of , , and away from the right-hand side surface of the superlens. Although small aberrations are visible in the field patterns, the simulations clearly demonstrate that superlensing effect occurs in these PhC slab systems. An infinitely long slab can focus all the Fourier components of a subject within this frequency range (including both propagation waves and evanescent waves), so that to produce a perfect image [22]. The flat lens property of AANR PhCs has many important potential applications, such as high-capacity optical storage systems [23].

(a)

(b)

(c)

Figure 4: TE polarization of a point source and its image across a PhC slab of square lattice. The normalized frequency is . (a) 8-layer PhC slab, (b) 12-layer PhC slab, and (c) 20-layer PhC slab.

3. Negative Refraction with Effective Negative Index

Distinguished from AANR PhCs, there is another kind of PhCs that can also give rise to negative refraction. It is referred to as effective-negative-index PhCs. In this part we consider the negative refraction property of PhC of 2D hexagonal lattices that are made of air holes. The background material is silicon (), and the radius of the air holes is 0.4a, where a is the lattice constant. The band curves for TM polarization (in-plane magnetic field vector) are shown in Figure 5. A negative refraction region is identified in the vicinity of where the second band curve and the light line intersect.

Figure 5: Band curves of 2D PhC of hexagonal lattice for TM polarization. The lattices are made of air holes of radius in a silicon host material with . The light line shifted to Γ point is shown in black. The red line indicates the frequency range of negative refraction.

Shown in Figure 6 are EFCs of TM mode for the normal frequency range between and . The shapes of EFCs for the second band are almost circular in the frequency range between and , where higher frequency contours concentrate near the Γ point (the center of circles). The blue dotted circle is the line of intersection of the second band surface and the light cone. The radius of the dotted circle is equal to the radius of EFC of air cone at .

(a)

(b)

Figure 6: Band structure and EFCs of the second band of the 2D PhC of hexagonal lattice for TM polarization. Values of frequency are in units of . (a) Surfaces of the second band and light cone. (b) EFCs of the 2nd band. The blue dotted circle is the line of intersection of the second band surface and the light cone.

The EFC of frequency is plotted in Figure 7. Light of this frequency is incident from air () into the model PhC that has an effective negative index (as shown in Figure 7). In the PhC, the energy velocity vector coincides with the group velocity vector of the Bloch mode . A Bloch mode denoted by A is excited in the PhC along the direction perpendicular to the EFC and point to Γ. Given the conservation of tangential wave vector, negative refraction occurs at the interface between air and PhC.

Figure 7: The EFC of frequency . The air-PhC interface is normal to the Γ-M direction. is the incident angle of the plane wave. A Bloch wave denoted by A is excited in the PhC. is the group velocity of the Bloch wave, which has a refraction angle of against the normal of the interface.

Using Snell’s law, light propagation satisfies where is the incident angle in air and is the refraction angle in the PhC. Because of the negative refraction, which leads to . It can be inferred that the absolute value of approximates to the ratio of the radius of the PhC contour to the radius of the air contour. For the frequency , the radius of the PhC contour equals the radius of the air contour (represented by the blue dotted line in Figure 6(b)), and the effective negative index is therefore . It can be derived that when frequency is higher than , . For frequency lower than , . The effective negative index of the PhC varies with the working frequency.

The homogeneity of refractive index of a PhC with effective negative index can be judged by its property of imaging as a flat lens. An isotropic homogeneous medium with refractive index of creates an image of a point source on the other side of the slab at a distance which depends on the distance of the source to the slab. The relation between source distance and image distance is illustrated directly from the simple ray diagram in Figure 8. There is no reflection at the interface between a negative index material of and air (). The coupling coefficient for any incident angle is always 100%. On the other hand, for a PhC with effective negative index of , the situation is slightly different. The FDTD simulations of imaging capability of a PhC slab with effective negative index of are shown in Figure 9. The imaging property is similar to that of the isotropic homogeneous medium with refractive index of . It is found that the image distance varies with the source distance in a similar way, with the sum of source distance, and image distance equals thickness of the PhC slab. On the other hand the transmission rate at the air-PhC interfaces in this case is highly angular dependent [24]. Incident wave cannot couple completely to the PhC. This is evident in Figure 9, where the field on the right-hand side of the slab is weaker than that on the left-hand side. From these simulations it can be concluded that the proposed PhC is not an isotropic homogeneous medium.

(a)

(b)

(c)

(d)

Figure 8: Ray tracing illustration of the flat lens made of an isotropic homogeneous medium with index of refraction . is the distance between the source and the first interface, and is the distance between the image point and the second interface, d is the thickness of the slab. (a) , (b) , (c) , (d) source is placed inside of the slab.

(a)

(b)

(c)

Figure 9: Imaging of a point source by a slab of 2D PhC of hexagonal lattice. The interface is normal to the Γ-M direction, and the normalized frequency is . The distance between the source and the first interface is , the distance between the image and the second interface is , and the thickness of the slab is . (a) , (b) , (c) .

To examine the frequency dependence of effective negative index , we employ the FDTD method to simulate imaging of a point source of different frequency by a finite-width PhC slab. The results shown in Figure 10 indicate that source of different frequency is focused to a different position. The point source locates at a fixed position in air, with a distance to the slab being half of the slab thickness. For the dispersive case of which corresponds to the aforementioned , the image distance outside the slab is equal to the source distance as shown in Figure 10(a). For a lower frequency which corresponds to , the image distance is slightly shorter than the source distance as shown in Figure 10(b). On the other hand, for a higher frequency which corresponds to , the image distance is slightly longer than the source distance as shown in Figure 10(c).

(a)

(b)

(c)

Figure 10: The FDTD simulations of point sources of three different frequencies placed in front of a 2D PhC slab of hexagonal lattices. The thickness of the slab is , and the distance between source and the first slab interface is . The air-PhC interface is perpendicular to the Γ-M direction, and the normalized frequencies are, respectively, (a) , (b) , and (c) .

4. Conclusions

In this paper we demonstrated two types of negative refraction mechanisms of PhC, namely the all-angle negative refraction without a negative index and the negative refraction with effective negative index. Light focusing effect occurs in both PhC slab systems. The AANR PhC slab can focus all Fourier components of an image within the frequency range of negative refraction. In the PhC with effective negative index, negative refraction has high angular dependence. The AANR PhC is more suitable for high-resolution superlensing application.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant no. 60588502. One of the authors (P. Jiang) would like to thank her parents and husband for understanding and spiritual support.

References

S. Foteinopoulou, E. N. Economou, and C. M. Soukoulis, “Refraction in media with a negative refractive index,” Physical Review Letters, vol. 90, no. 10, Article ID 107402, 4 pages, 2003.

V. G. Veselago, “Photomagnetism,” Soviet Physics—Uspekhi, vol. 10, p. 509, 1968.

G. Gbur, H. F. Schouten, and T. D. Visser, “Achieving superresolution in near-field optical data readout systems using surface plasmons,” Applied Physics Letters, vol. 87, no. 19, Article ID 191109, pp. 1–3, 2005.

P. Markoš and C. M. Soukoulis, “Transmission studies of left-handed materials,” Physical Review B, vol. 65, no. 3, Article ID 033401, 4 pages, 2001.

J. D. Wilson and Z. D. Schwartz, “Multifocal flat lens with left-handed metamaterial,” Applied Physics Letters, vol. 86, no. 2, Article ID 021113, 3 pages, 2005.

M. Rajput, R. K. Sinha, and S. Rawal, “Negative refraction in visible region using nano-structured metallo-dielectric photonic crystal,” in Photonic Fiber and Crystal Devices: Advances in Materials and Innovations in Device Applications III, vol. 7420 of Proceedings of SPIE, p. 74, San Diego, Calif, USA, August 2009.

G. Sun, A. S. Jugessur, and A. G. Kirk, “Saturated super-resolution imaging of photonic crystal with negative refraction,” in Optical Data Storage, vol. 6282 of Proceedings of SPIE, Bellingham, Wash, USA, April 2006.

X. Zhang, “Image resolution depending on slab thickness and object distance in a two-dimensional photonic-crystal-based superlens,” Physical Review B, vol. 70, no. 19, Article ID 195110, 9 pages, 2004.

Z. Li, B. Liang, H. Guo, J. Chen, and S. Zhuang, “Negative refraction and birefringence in a two-dimensional flat perfect photonic crystal,” in Optoelectronic Materials and Devices II, vol. 6782 of Proceedings of SPIE, Wuhan, China, November 2007.

S. Feng, Z. Y. Li, Z. F. Feng, B. Y. Cheng, and D. Z. Zhang, “Engineering the imaging properties of a metallic photonic-crystal slab lens,” Applied Physics Letters, vol. 88, no. 3, Article ID 031104, 3 pages, 2006.

S. Feng, Z. Y. Li, Z. F. Feng, K. Ren, B. Y. Cheng, and D. Z. Zhang, “Focusing properties of a rectangular-rod photonic-crystal slab,” Journal of Applied Physics, vol. 98, no. 6, Article ID 063102, 6 pages, 2005.

C. W. Wong, R. Chatterjee, K. Liu, C. J. Chen, and C. A. Husko, “Negative refraction and nonlinearities in photonic band gap nanostructures,” in The International Society for Optical Engineering, vol. 6327 of Proceedings of SPIE, San Diego, Calif, USA, August 2006.

S. Anantha Ramakrishna, “Physics of negative refractive index materials,” Reports on Progress in Physics, vol. 68, no. 2, pp. 449–521, 2005.

M. Notomi, “Negative refraction in photonic crystals,” Optical and Quantum Electronics, vol. 34, no. 1–3, pp. 133–143, 2002.

R. Moussa, S. Foteinopoulou, L. Zhang, G. Tuttle, K. Guven, E. Ozbay, and C. M. Soukoulis, “Negative refraction and superlens behavior in a two-dimensional photonic crystal,” Physical Review B, vol. 71, no. 8, pp. 1–85106, 2005.

E. Cubukcu, K. Aydin, E. Ozbay, S. Foteinopolou, and C. M. Soukoulis, “Subwavelength resolution in a two-dimensional photonic-crystal-based superlens,” Physical Review Letters, vol. 91, no. 20, Article ID 207401, 4 pages, 2003.

A. Berrier, M. Mulot, M. Swillo, M. Qiu, L. Thylén, A. Talneau, and S. Anand, “Negative refraction at infrared wavelengths in a two-dimensional photonic crystal,” Physical Review Letters, vol. 93, no. 7, Article ID 10.1103/PhysRevLett.93.073902, 4 pages, 2004.

X. Zhang, “Tunable non-near-field focus and imaging of an unpolarized electromagnetic wave,” Physical Review B, vol. 71, no. 23, Article ID 235103, 5 pages, 2005.

Z. Lu, C. Chen, C. A. Schuetz, S. Shi, J. A. Murakowski, G. J. Schneider, and D. W. Prather, “Subwavelength imaging by a flat cylindrical lens using optimized negative refraction,” Applied Physics Letters, vol. 87, no. 9, Article ID 091907, 3 pages, 2005.

C. Luo, S. G. Johnson, M. Soljačić, J. D. Joannopouls, and J. B. Pendry, “Novel optical phenomena with photonic crystals,” in The International Society for Optical Engineering, vol. 5166 of Proceedings of SPIE, pp. 207–209, 2004.

C. Luo, S. G. Johnson, J. D. Joannopoulos, and J. B. Pendry, “All-angle negative refraction without negative effective index,” Physical Review B, vol. 65, no. 20, Article ID 201104, 4 pages, 2002.

J. B. Pendry, “Negative refraction makes a perfect lens,” Physical Review Letters, vol. 85, no. 18, pp. 3966–3969, 2000.

L. Liu and S. He, “Near-field optical storage system using a solid immersion lens with a left-handed material slab,” Optics Express, vol. 12, no. 20, pp. 4835–4840, 2004.

Z. Ruan, Dispersion engineering: negative refraction and designed surface plasmons in periodic structures, Ph.D. thesis, Microelectronics and Applied Physics, Stockholm, Sweden, 2007.