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

Extremely small wavevector regime in a one-dimensional photonic crystal heterostructure for angular transmission filtering

Open Access Open Access

Abstract

One-dimensional photonic crystals (PCs), when operating near the band edge in the dispersion diagram, inherently possess nearly polarization-independent angular selectivity—an angular transmission window around the normal direction with reflection for other angles. However, the incident light is mostly reflected at the PC–air interface due to large impedance mismatch. We show that the reflection may be sufficiently suppressed by utilizing a specially designed antireflection structure consisting of a PC having a different pitch from that of the host PC. The underlying mechanism is that the interfaces of the antireflection PC with the host PC and the air structure are selected such that the transverse impedance has a real value, which is positioned at the center of the thickness of a material film. Moreover, our structure provides a high-throughput wide angular transmission window, including the normal direction in both s and p polarizations. We develop an analytical model that captures the angular selectivity observed in numerical results.

© 2016 Optical Society of America

One-dimensional (1D) photonic crystals (PCs) that consist of an alternative stack of different material films provide an opportunity to manipulate the dispersion relation [1] and have found numerous optical and thermal applications [27]. Particularly, angular selectivity is among the important functionalities of PCs and has been extensively investigated [8,9]. In Ref. [10], an angular filter for p polarization has been experimentally realized in the full visible spectrum by utilizing Brewster modes, providing an angular transmission window of 8° in an oblique direction of 55°. In Refs. [11,12], a transmission window for p polarization has been set to the normal direction through the use of anisotropic layers. In Ref. [13], an ultrabroadband plasmonic Brewster transmission for p polarization has been presented by utilizing a metallic grating. In general, a polarization-independent angular transmission window in the normal direction is desirable for many application scenarios, such as optical sensor systems [14] and communication systems [15]. However, to the best of our knowledge, there has been no work on a polarization-independent angular filter having a high-throughput transmission window in the normal direction.

PCs exhibit polarization-independent omnidirectional reflection when operating in the bandgap of the dispersion diagram [16]. Using such bandgap operation, a defect layer in PCs allows light transmission at a specific frequency as a Fabry–Perot resonator [1,1721], and the angular transmission window is very narrow due to the intrinsic nature of resonance. PCs inherently possess a nearly polarization-independent angular transmission window that includes the normal direction when operating near the band edge in the dispersion diagram. However, the incident light is mostly reflected back at the PC–air interface due to large impedance mismatch. The wave impedance of PCs is dependent on local observation points, and, thus, the antireflection in the extremely small wavevector regime is challenging.

In this Letter, we show analytically and numerically that the reflection of the incident light is sufficiently suppressed in the extremely small wavevector regime by utilizing an antireflection PC having a different pitch from that of the host PC, resulting in a high-throughput angular transmission filter. The underlying mechanism is that the antireflection PC needs to have the interfaces with the host PC and the air structure such that the transverse impedance at each interface between two semi-infinite media has a real value, which allows us to design the antireflection PC by using the general methodology of impedance matching. In Ref. [22], we presented an antireflection structure for a two-dimensional PC in the extremely small wavevector regime. However, it showed only a sharp angular and frequency response.

We start by considering the dispersion diagram of a PC, as shown in Figs. 1(a) and 1(b), obtained by using the plane wave expansion (PWE) method [23]. The PWE method is widely used for analyses of PCs in which Maxwell’s equations are exactly solved by expanding the electromagnetic fields with a plane-wave basis set, and, thus, the wavevectors and electromagnetic fields are accurately evaluated. Two material films are alternatively and infinitely stacked in the PC [inset of Fig. 1(a)]. We assume TiO2 and SiO2 for two materials with refractive indices of nA=2.4 and nB=1.45, respectively. The thicknesses of the two materials are selected at tA=214nm and tB,h=90nm such that the PC has an angular transmission window of 2θd=12° [oblique black dashed lines closer to the vertical axis in Fig. 1(a)] for s and p polarizations at the design frequency of fd=500THz (horizontal pink solid line). We note that extended modes that are allowed to propagate exist within angles of less than θd=6° (shaded region) for both polarizations at fd. In other angles at fd, no modes are allowed to propagate (white region) except for large angles θ>70° in p polarization. Therefore, we expect that the PC may work as an angular filter having a nearly polarization-independent transmission window of 2θd around the normal direction, and a reflection in other angles. The effect of the leakage of electromagnetic waves in such large angles (θ>70°) in p polarization may be reduced by optimizing the refractive indices and the thicknesses of materials.

 figure: Fig. 1.

Fig. 1. (a) Dispersion diagram of a 1D PC in s and p polarizations. Extended modes exist in the shaded regions. The oblique black dashed lines represent the design angle of θd=6° and light lines, respectively. The horizontal pink solid line represents the design frequency of fd=500THz. The diagram is enlarged around fd in (b). (c) Configuration of a PC having angular selectivity. (d) Two-dimensional map of transmission coefficients obtained numerically by using Eqs. (1a)–(1d) and the transfer matrix method. Angular dependence of transmission coefficients at (e) 500 THz [horizontal pink solid line in (d)] and at (f) 515 THz [horizontal blue dashed line in (d)]. [blue solid lines, numerical results, Eqs. (1a)–(1d); pink dots, analytical results, Eqs. (3a)–(3d)] (parameters: nA=2.4, nB=1.45, tA=214nm, tB,h=90nm, ah=tA+tB,h=304nm, tB,ar=93nm, NAB=97, m=15).

Download Full Size | PDF

Figure 1(c) shows the configuration of a PC having the angular selectivity function. Two material films with refractive indices of nA and nB are alternatively stacked. The films are numbered from j=1 to NAB, where NAB is the total number of layers. From the viewpoint of the underlying mechanism, the heterostructure consists of the host PC (j=m to NABm+1) and the antireflection PCs at both sides (j=1 to m and j=NABm+1 to NAB), where each of two layers j=m and NABm+1 has the interface between the antireflection PC and the host PC positioned at the center of the thickness of the layer [indicated by red dotted lines F in Fig. 1(c)] and is counted for the layer numbers of both the host and antireflection PCs. From the geometrical point of view, material A has the same thickness tA=214nm in both the host PC and the antireflection PC except for the first layer j=1 and the last layer j=NAB. The two layers at both ends have the half-thickness tA/2 such that the antireflection PC has the same impedance at the interface with the air structure and at the interface with the host PC [indicated by red dotted lines F in Fig. 1(c)]. Material B, however, has two different thicknesses in the entire structure: tB,h=90nm for the host PC and tB,ar=93nm for the antireflection PC.

The monochromatic electromagnetic waves with s polarization or p polarization are impinged on the entire PC from the air, where the electric field or the magnetic field is perpendicular to the incident plane, respectively. The Fresnel reflection and transmission coefficients of the electric field for s polarization and the magnetic field for p polarization at the interface between layers j and j+1 are given by [24]

rs,j,j+1=kx,jkx,j+1kx,j+kx,j+1,
ts,j,j+1=2kx,jkx,j+kx,j+1,
rp,j,j+1=kx,jnj2kx,j+1nj+12kx,jnj2+kx,j+1nj+12,
tp,j,j+1=2kx,jnj2kx,jnj2+kx,j+1nj+12,
where kx,j=nj2k02ky2 and ky=k0sin(θ) are the x (normal to the surface of the layers) and the y (parallel to the surface of the layers) components of the wavevector, k0 is the free-space wavenumber, and nj is the refractive index of the jth layer. The air regions are labeled by j=0 and NAB+1 in Fig. 1(c). The subscripts s and p represent the polarization. The reflection and transmission coefficients R, T of the entire structure of PC in air can be calculated by using Eqs. (1a)–(1d) and the transfer matrix method [25], which is widely used for analyses of multilayer structures. We assume the lossless model, |R|2+(kx,NAB+1/kx,0)(τ0/τNAB+1)|T|2=1, where τj=1 for s polarization and τj=nj2 for p polarization, and plot the power transmission coefficient (kx,NAB+1/kx,0)(τ0/τNAB+1)|T|2 for numerical results throughout the Letter.

Figure 1(d) shows the two-dimensional map of the transmission coefficients of the entire structure of Fig. 1(c) with variations of the incident angle and the operating frequency. Due to the properly designed antireflection structure, high transmission (bright-color region) is obtained for both polarizations in the propagation-allowed region that is characterized by the dispersion diagram (shaded regions) in Figs. 1(a) and 1(b). It should be emphasized that the reflection is sufficiently suppressed in a wide range of frequency in the propagation-allowed region; our structure provides not only a transmission window of 2θd at fd=500THz but also wide angular transmission window at high frequency. When θ>70° for p polarization around fd=500THz in Fig. 1(d), the incident light is leaked out, as observed from the dispersion diagram.

We plot the transmission coefficient as a function of the incident angle at fd=500THz in Fig. 1(e). We clearly observe that our structure has efficient angular selectivity; light passes through around the normal direction and is reflected over other angles. The half-power angular width 2θ is 11° for s polarization and 13° for p polarization, respectively. In addition, at a frequency of 515 THz in Fig. 1(f), a wide transmission window appears: 57° for s- polarization and 69° for p polarization. Interestingly, sufficient antireflection is maintained when the frequency goes up. This wideband characteristic of the sufficient antireflection comes from the fact that the antireflection structure consists of a PC; the antireflection structure has a trend in the frequency response similar to that of the host medium.

We next elucidate how sufficient antireflection, as discussed above, is obtained. The underlying mechanism can be understood in the impedance behavior. Consider the impedance behavior at fd=500THz in Fig. 2. The transverse impedance for s or p polarization, respectively, is defined as

Zs(x)=Ez(x)Hy(x),
Zp(x)=Ey(x)Hz(x).

 figure: Fig. 2.

Fig. 2. Transverse impedance of the host PC at fd=500THz. The transverse impedance is normalized with respect to the free-space intrinsic impedance Z0. (a) Locally dependent distribution of the real part (blue solid line) and the imaginary part (pink dashed line) along the x axis at normal incidence. (b) Angular dependence of the real part of the transverse impedance Zh of the host PC (blue solid line) at the location F. Zar2/Zair, which comes from the antireflection condition, is plotted as the green dashed line. The imaginary parts of the impedances Zh, Zar, and Zair are zero (not presented).

Download Full Size | PDF

We assume infinite PCs for impedance analysis throughout this Letter and thus obtain each transverse impedance in Eqs. (2a) and (2b) by the PWE method. Figure 2(a) shows the locally dependent transverse impedance of the host PC at normal incidence. The transverse impedance is normalized with respect to the free-space intrinsic impedance Z0. We point out that the imaginary part of the impedance is zero at the center of the thickness of each film [indicated by red dotted lines F and H in Fig. 2(a)]. This allows us to obtain sufficient antireflection by using the general methodology of impedance matching. However, at interface G, which is a regular interface between layers, the impedance has a non-zero value in the imaginary part, and, thus, antireflection is in general difficult. In our case, we have selected interface F since the real part of the impedance at interface F is smaller than that at interface H.

Figure 2(b) shows the angular dependence of the transverse impedance Zh (blue solid line) of the host PC at interface F. The real part of the transverse impedance is presented since the imaginary part is zero even for oblique incidence. The impedance has a slowly varying characteristic within ±5° for both polarizations. In addition, we calculate the transverse impedance Zar of the antireflection PC that has a different thickness of tB,ar=93nm for SiO2 films with the same thickness of tA=214nm for TiO2 films and plot Zar2/Zair in the same figure, where Zair=Z0/cos(θ) for s polarization and Zair=Z0cos(θ) for p polarization. We verify that the general antireflection condition Zar2ZhZair [25] is satisfied within ±5° in Fig. 2(b). The length of the antireflection PC has been selected at a quarter-guided-wavelength (number of layers: m=15) at fd according to the general antireflection condition [25]. Therefore, the reflection at the interface of the entire structure has been sufficiently suppressed, as shown in the high transmission coefficient of Fig. 1(e).

Based on the underlying mechanism above, we now develop an analytical model, as shown in Fig. 3. The heterostructure consists of three layers; the host PC having the antireflection PCs at both sides in air. The reflection and transmission coefficients at the interfaces of the analytical model of Fig. 3 are given by

rs,j,j+1=Zj+1ZjZj+1+Zj,
ts,j,j+1=2Zj+1Zj+1+Zj,
rp,j,j+1=ZjZj+1Zj+Zj+1,
tp,j,j+1=2ZjZj+Zj+1,
where the host PC, the antireflection PC, and the air structure are labeled j=2, j=1,3, and j=0,4, respectively. The transverse impedances Zj and wavevectors (kx,j,ky) are obtained in each polarization by the PWE method, and then the reflection and transmission coefficients at each interface are obtained using Eqs. (3a)–(3d). In addition, using the transfer matrix method, we calculate the reflection and transmission coefficients of the entire structure and plot the transmission coefficients in Figs. 1(e) and 1(f). The analytical results (pink dots) agree excellently with numerical results (blue solid lines). Therefore, we verify the underlying mechanism of our structure of Fig. 1(c). When θ>70° for p polarization at fd=500THz [Fig. 1(e)], the analytical result has two peaks (pink dots), while the numerical result has four peaks (blue solid line). This discrepancy is due to the assumption of the infinite PCs in our analysis; the angular range, where the transmission peaks occur in the finite PC of Fig. 1(c) in the numerical result, falls in the bandgap in the infinite PC in the analytical result.

 figure: Fig. 3.

Fig. 3. Analytical model of the heterostructure of Fig. 1(c). Wavevectors and transverse impedances are obtained by the PWE method. The lengths of the equivalent media are tar=(m1)(tA+tB,ar)2=2149nm and th=(NAB2m+1)ah2=10,336nm, respectively.

Download Full Size | PDF

We have used a number of NAB=97 layers in the entire structure. The number of layers can be reduced; only the number of layers in the host PC is reduced, while that in the antireflection PCs (m=15) is maintained. Reducing the number of layers gives gradual slope in the angular response.

In conclusion, we have explored a 1D PC that operates in the extremely small wavevector regime and, therefore, has an angular transmission window around the normal direction. Based on the analysis of the locally dependent impedance, the interfaces of an antireflection PC with the host PC and the air structure have been selected at the center of the thickness of a material film. This allowed us to use the general method of impedance matching. Consequently, the reflection at the PC–air interface was sufficiently suppressed. In addition, numerical results have revealed that our PC structure provides a high-throughput wide angular transmission window. The current nanofabrication technology for multilayers can be used for our structure. Therefore, our results nicely extend the range of potential applications of 1D PCs. For example, the signal to interference and noise ratio in optical sensor systems and communication systems can be enhanced by suppressing noise and unwanted signals coming from other directions. Based on our antireflection mechanism, one can optimize a desirable angular and frequency response for desired applications.

REFERENCES

1. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Modeling the Flow of Light, 2nd ed. (Princeton University, 2008).

2. H. A. Macleod, Thin-Film Optical Filters, 4th ed. (CRC Press, 2010).

3. B. E. Nelson, M. Gerken, D. A. B. Miller, R. Piestun, C.-C. Lin, and J. S. Harris Jr., Opt. Lett. 25, 1502 (2000). [CrossRef]  

4. R. Ozaki, T. Matsui, M. Ozaki, and K. Yoshimoto, Appl. Phys. Lett. 82, 3593 (2003). [CrossRef]  

5. K.-T. Lee, C. Ji, D. Banerjee, and L. J. Guo, Laser Photon. Rev. 9, 354 (2015). [CrossRef]  

6. B. J. Lee, C. J. Fu, and Z. M. Zhang, Appl. Phys. Lett. 87, 071904 (2005). [CrossRef]  

7. A. P. Raman, M. A. Anoma, L. Zhu, E. Rephaeli, and S. Fan, Nature 515, 540 (2014). [CrossRef]  

8. L. Maigyte and K. Staliunas, Appl. Phys. Rev. 2, 011102 (2015). [CrossRef]  

9. Y. Shen, C. W. Hsu, Y. X. Yeng, J. D. Joannopoulos, and M. Soljacic, Appl. Phys. Rev. 3, 011103 (2016). [CrossRef]  

10. Y. Shen, D. Ye, I. Celanovic, S. G. Johnson, J. D. Joannopoulos, and M. Soljacic, Science 343, 1499 (2014). [CrossRef]  

11. R. E. Hamam, I. Celanovic, and M. Soljacic, Phys. Rev. A 83, 035806 (2011). [CrossRef]  

12. Y. Shen, D. Ye, L. Wang, I. Celanovic, L. Ran, J. D. Joannopoulos, and M. Soljacic, Phys. Rev. B 90, 125422 (2014). [CrossRef]  

13. A. Alu, G. D. Aguanno, N. Mattiucci, and M. J. Bloemer, Phys. Rev. Lett. 106, 123902 (2011). [CrossRef]  

14. C. Niclass, M. Soga, H. Matsubara, S. Kato, and M. Kagami, IEEE J. Solid-State Circuits 48, 559 (2013). [CrossRef]  

15. T. Komine and M. Nakagawa, IEEE Trans. Consum. Electron. 50, 100 (2004). [CrossRef]  

16. Y. Fink, J. N. Winn, S. Fan, C. Chen, J. Michel, J. D. Joannopoulos, and E. L. Thomas, Science 282, 1679 (1998). [CrossRef]  

17. B. Shi, Z. M. Jiang, and X. Wang, Opt. Lett. 26, 1194 (2001). [CrossRef]  

18. G. Liang, P. Han, and H. Wang, Opt. Lett. 29, 192 (2004). [CrossRef]  

19. K.-Y. Xu, X. Zheng, C.-L. Li, and W.-L. She, Phys. Rev. E 71, 066604 (2005). [CrossRef]  

20. H. Nemec, P. Kuzel, L. Duvillaret, A. Pashkin, M. Dressel, and M. T. Sebastian, Opt. Lett. 30, 549 (2005). [CrossRef]  

21. Y. Takeda, H. Iizuka, S. Mizuno, K. Hasegawa, T. Ichikawa, H. Ito, T. Kajino, A. Ichiki, and T. Motohiro, J. Appl. Phys. 116, 014501 (2014). [CrossRef]  

22. H. Iizuka and N. Engheta, Phys. Rev. B 90, 115412 (2014). [CrossRef]  

23. K. M. Ho, C. T. Chan, and C. M. Soukoulis, Phys. Rev. Lett. 65, 3152 (1990). [CrossRef]  

24. W. C. Chew, Waves and Fields in Inhomogeneous Media (Van Nostrand Reinhold, 1990).

25. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, 2nd ed. (Wiley, 2007).

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 (3)

Fig. 1.
Fig. 1. (a) Dispersion diagram of a 1D PC in s and p polarizations. Extended modes exist in the shaded regions. The oblique black dashed lines represent the design angle of θ d = 6 ° and light lines, respectively. The horizontal pink solid line represents the design frequency of f d = 500 THz . The diagram is enlarged around f d in (b). (c) Configuration of a PC having angular selectivity. (d) Two-dimensional map of transmission coefficients obtained numerically by using Eqs. (1a)–(1d) and the transfer matrix method. Angular dependence of transmission coefficients at (e) 500 THz [horizontal pink solid line in (d)] and at (f) 515 THz [horizontal blue dashed line in (d)]. [blue solid lines, numerical results, Eqs. (1a)–(1d); pink dots, analytical results, Eqs. (3a)–(3d)] (parameters: n A = 2.4 , n B = 1.45 , t A = 214 nm , t B , h = 90 nm , a h = t A + t B , h = 304 nm , t B , ar = 93 nm , N A B = 97 , m = 15 ).
Fig. 2.
Fig. 2. Transverse impedance of the host PC at f d = 500 THz . The transverse impedance is normalized with respect to the free-space intrinsic impedance Z 0 . (a) Locally dependent distribution of the real part (blue solid line) and the imaginary part (pink dashed line) along the x axis at normal incidence. (b) Angular dependence of the real part of the transverse impedance Z h of the host PC (blue solid line) at the location F . Z ar 2 / Z air , which comes from the antireflection condition, is plotted as the green dashed line. The imaginary parts of the impedances Z h , Z ar , and Z air are zero (not presented).
Fig. 3.
Fig. 3. Analytical model of the heterostructure of Fig. 1(c). Wavevectors and transverse impedances are obtained by the PWE method. The lengths of the equivalent media are t ar = ( m 1 ) ( t A + t B , ar ) 2 = 2149 nm and t h = ( N A B 2 m + 1 ) a h 2 = 10,336 nm , respectively.

Equations (10)

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

r s , j , j + 1 = k x , j k x , j + 1 k x , j + k x , j + 1 ,
t s , j , j + 1 = 2 k x , j k x , j + k x , j + 1 ,
r p , j , j + 1 = k x , j n j 2 k x , j + 1 n j + 1 2 k x , j n j 2 + k x , j + 1 n j + 1 2 ,
t p , j , j + 1 = 2 k x , j n j 2 k x , j n j 2 + k x , j + 1 n j + 1 2 ,
Z s ( x ) = E z ( x ) H y ( x ) ,
Z p ( x ) = E y ( x ) H z ( x ) .
r s , j , j + 1 = Z j + 1 Z j Z j + 1 + Z j ,
t s , j , j + 1 = 2 Z j + 1 Z j + 1 + Z j ,
r p , j , j + 1 = Z j Z j + 1 Z j + Z j + 1 ,
t p , j , j + 1 = 2 Z j Z j + Z j + 1 ,
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.