Photonic crystals (PCs) are artificial microstructures with periodic refractive index.[
Chinese Physics B, Volume. 29, Issue 10, (2020)
Thermal tunable one-dimensional photonic crystals containing phase change material
To obtain the adjustable photonic crystals (PCs), we numerically investigate one-dimensional (1D) PCs with alternating VO2 and SiO2 layers through transfer matrix method. The dispersion relation agrees well with the transmittance obtained by the finite element calculation. Tunable band gaps are achieved with the thermal stimuli of VO2, which has two crystal structures. The monoclinic crystal structure VO2 (R) at low temperature exhibits insulating property, and the high temperature square rutile structure VO2 (M) presents metal state. Concretely, the bandwidth is getting narrower and red shift occurs with the higher temperature in VO2 (R)/SiO2 PCs structure. Based on the phase change characteristics of VO2, we can flexibly adjust the original structure as VO2 (R)/VO2 (M)/SiO2. By increasing the phase ratio of VO2 (R) to VO2 (M), the band gap width gradually becomes wider and blue shift occurs. The discrete layers of gradient composites on the dispersion of 1D PCs are also investigated, which enhances the feasibility in practical operation. Thus, our proposed thermal modulation PCs structure paves a new way to realize thermal tunable optical filters and sensors.
1. Introduction
Photonic crystals (PCs) are artificial microstructures with periodic refractive index.[
The emergence of liquid crystal,[
In this paper, we have proposed 1D PCs structure with alternating VO2 and SiO2 layers to achieve a tunable band gap under thermal stimuli. The photonic band gap is numerical calculated with the transfer matrix method and agrees well with the finite element method (FEM). The band gap position and bandwidth of the 1D PCs are fully explored with the phase change of VO2 from insulation phase to metal phase. Two sets of the proposed PCs structures are calculated for comparison with different portions of the phase states.
2. Theoretical methods
VO2 has two crystal structures.[
Figure 1.Real part of the effective dielectric constant of VO2 at different temperature in THz region.
3. Results and discussion
Figure 2 shows the 3D configuration and the cross section view of the proposed PCs structure, which is composed of alternating VO2 (R) and SiO2. The incident light is in TE mode along the z-axis. N denotes the period and is taken as 10. The lattice constant d is 2 μm. The thickness ratio of VO2 (R) to SiO2 is a : b = 3 : 2. The refractive index of SiO2 is taken as 1.5 and the dielectric constant of VO2 is given by Eq. (1).
Figure 2.The 3D configuration and the cross section view of the VO2 (R)/SiO2 PCs structure.
The dispersion relationship and transmittance of VO2 (R)/SiO2 PCs are illustrated in Fig. 3. It is obviously found that three photonic band gaps exist in the range of 0–125 THz in Fig. 3(a). The first band gap locates at 26.1–37.1 THz. The second bang gap is in the region between 56.5 THz and 73.1 THz. The third band gap is from 91.5 THz to 106.6 THz. As a result, the electromagnetic wave is not allowed to propagate through the VO2 (R)/SiO2 PCs in this circumstance. The corresponding transmittance of the structure is numerically shown in Fig. 3(b). It is clearly found that the transmissivity is almost zero in the photonic band gap and achieves a good self-consistent with the band gap characteristics of the PCs, verifying the correctness of our results. Besides, the transmission spectrum range is not limited to 0–125 THz. It can be extended to high frequency region with careful chosen thickness ratio. The working frequency region can be extended to 125–200 THz while the relevant parameters are taken as d = 2 μm and the thickness ratio of VO2 (R) to SiO2 is a : b = 3 : 7. Thus, with a smaller thickness ratio between VO2 (R) and SiO2, the transmission spectrum will have a significant blue shift.
Figure 3.(a) The dispersion relation of VO2 (R)/SiO2 PCs structure obtained with transfer matrix method. (b) The transmission spectrum of the structure obtained with numerical simulation.
The optical properties of VO2 are closely related to the external temperature. Based on the phase transition characteristics of VO2, the influences of temperature on the dispersion relationship are explored in Fig. 4. A transition from the insulating phase to the metal phase occurs at the critical temperature 341 K. In this process, there is a red shift in the dispersion relationship, and the band gap bandwidth is gradually decreasing. The first and second band gaps disappear when VO2 is converted to metal phase at 350 K. With the phase change of VO2 from insulating to metal phase, it is difficult for the low frequency electromagnetic waves to pass through the metallic interface due to multiple scattering overlapping. Thus, the photon band gap at low frequency will gradually disappear.
Figure 4.The influence of temperature on the dispersion relation of VO2 (R)/SiO2 PCs in the range of 0–125 THz.
To enrich the tunable properties of PCs, we make further adjustments to the VO2 (R)/SiO2 PCs structure. Now the VO2 is divided into two parts (VO2 (R)/VO2 (M)). Figure 5(a) depicts 3D configuration and cross section view of the structure, which is composed of three alternating layers VO2 (R)/VO2 (M)/SiO2. a1 and a2 are the thicknesses of VO2 (R) and VO2 (M). The parameters are chosen as follows: d = 2 μm, a : b = 3 : 2, and N = 10. Figure 5(b) illustrates the dispersion relation and transmission spectrum with the thickness ratio a1 : a2 = 5 : 1. Compared to Fig. 3, the dispersion curve demonstrates a red shift, and the bandwidth of each band gap is getting smaller. Now the band gaps are in the ranges of 23.8–34.2 THz, 51.7–67.9 THz, and 89.5–95.8 THz respectively. In addition, the band gap will disappear in the higher frequency region. The band gap distribution in Fig. 5(b) is consistent with the transmittance in Fig. 5(c). It is revealed that the red shift of the dispersion relation and the decrease of the band gap are caused by the emergence of metal phase VO2 (M). For the combined layer of VO2 (R) and VO2 (M), it can be regarded as a composite material layer. Its effective dielectric constant is obtained through series capacitance method.[
Figure 5.(a) The 3D configuration and cross section view of the VO2 (R)/VO2 (M)/SiO2 PCs structure. (b) The dispersion curve obtained with transfer matrix method. (c) The transmission spectrum of VO2 (R)/VO2 (M)/SiO2 PCs structure obtained with numerical calculation.
To verify the influence of VO2 (M) on the dispersion relationship, the trend of the dispersion relationship through different thickness ratio of VO2 (R) to VO2 (M) is calculated in Fig. 6. The thickness ratio of a1 : a2 is from 0 to 1.0 with a thicker layer of VO2 (R) in Fig. 6(c) and it is getting larger than 1.0 with a thicker layer of VO2 (M) in Fig. 6(a). It can be intuitively found that the dispersion behavior of VO2 (R) and VO2 (M) is changing with different ratios. As the thickness of the VO2 (M) continues to decrease, the dispersion relationship demonstrates a blue shift and the width of each band gap gradually increases. According to the series capacitance method, the effective dielectric constant of the composite layer increases with the decrease of the thickness of the metal phase. Consequently, the dielectric constant of the composite layer is significantly different from SiO2, leading to significant band shift of PCs. Besides, the position and band gap width of the first band gap are extracted in Figs. 6(b) and 6(d). In Fig. 6(b), when the proportion of VO2 (M) is very low, the growth trend of both is relatively flat. When the thickness of the VO2 (M) is larger than that of the VO2 (R), the blue shift and growth trend are obvious in Fig. 6(d). After further calculation, we extract the position distributions of the second and the third band gaps with different ratios in Fig. 7. The results denote that the variations of the second and the third band gaps are the same as that of the first band gap. Namely, as the thickness of the metal phase decreases, the dispersion relationship blue shifts and the band gaps are getting smaller.
Figure 6.(a) and (c) The influence of different proportions of VO2 (R) and VO2 (M),
Figure 7.(a) and (b) The position and band gap of the second band with different
Finally, the number of layers on the dispersion of 1D PCs is considered in Fig. 8. VO2 is divided into n layers (n > 1), the first and the last layers are set at 303 K (VO2 (R)) and 343 K (VO2 (M)). As n increases, the temperature of each layer is evenly divided. The division is listed as follows: 2 layers, T = 303 K, 343 K; 3 layers, T = 303 K, 323 K, 343 K; 4 layers, T = 303 K, 318 K, 333 K, 343 K; 5 layers: T = 303 K, 313 K, 323 K, 333 K, 343 K. Figure 8(a) presents the dispersion relationship when VO2 is divided into different layers. When the temperature division becomes more intensive, the dispersion relationship undergoes a slight blue shift and the band gap width does not change significantly. Figure 8(b) shows the position and bandwidth of the first band gap. The band gap position is moving into higher frequency region with a larger n. The band gap width is slightly changed. Thus, to initiate the gradient temperature distribution, we can divide it into several separated parts with concrete thermal stimuli, which enhances the feasibility in practical operation.
Figure 8.(a) The dispersion curve of 1D PCs with different number of layers. (b) The position and the first band width with different layers of VO2.
4. Conclusion
In summary, theoretical calculation of 1D PCs consisting of alternating vanadium oxide and dielectric layer with the transfer matrix method is conducted in this work. Our results of the dispersion relation coincide well with the numerical calculation of transmittance. The position and width of the band gap can be dynamically tuned with the different phase ratios of VO2 (R) to VO2 (M). Thermal tunable photonic band gap can be achieved with different portions of the metallic phase of VO2. In addition, to initiate the gradient temperature distribution, we can divide it into several separated parts with concrete thermal stimuli, which enhances the feasibility in practical operation. Therefore, this work will guide some potential applications in the thermal tunable optical filters and sensors.
[1] H T Jiang, H Chen, H Q Li, Y W Zhang, J Zi, S Y Zhu. Phys. Rev. E, 69(2004).
[2] S Colodrero, M Ocana, H Míguez. Langmuir, 24, 4430(2008).
[3] E Yablonovitch. Phys. Rev. Lett., 58, 2059(1987).
[4] S John. Phys. Rev. Lett., 58, 2486(1987).
[5] M Notomi, E Kuramochi, H Taniyama. Opt. Express, 16(2008).
[6] M Fujita, S Takahashi, Y Tanaka, T Asano, S Noda. Science, 308, 1296(2005).
[7] S Noda, K Tomoda, N Yamamoto, A Chutinan. Science, 289, 604(2000).
[8] F Benabid, J C Knight, G Antonopoulos, P Russell, J St. Science, 298, 399(2002).
[9] S P Li, H J Liu, Q B Sun, N Huang. IEEE Photonic Tech. L, 27, 752(2015).
[10] W Y Wang, Y X Cui, Y R He, X Y Lin, X M Tian, T Ji, S L He. Opt. Lett., 39, 331(2014).
[11] X Zhu, X D Yang, X Wang. Prog Electromagn Res., 67, 103(2017).
[12] E Yablonovitch, T J Gmitter, K M Leung. Phys. Rev. Lett., 67, 2295(1991).
[13] G K Jamshidi, F Moslemi. Appl. Opt., 56, 4146(2017).
[14] A H Aly, H A Elsayed, A A Ameen, S H Mohamed. Int. J. Mod. Phys. B, 31(2017).
[15] I V Timofeev, D N Maksimov, A F Sadreev. Phys. Rev. B, 97(2018).
[16] P F Zhao, B Li, Z H Tang, Y Gao, H M Tian, H L Chen. Smart Mater. Struct., 28(2019).
[17] Y F Zhang, C C Chan, J Sun. Sens. Actuat. A-Phys., 157, 276(2010).
[18] K Busch, S John. Phys. Rev. Lett., 83, 967(1999).
[19] P Snapp, P Kang, J Leem, S W Nam. Adv. Funct. Mater., 29(2019).
[20] Y J Ke, I Balin, N Wang, Q Lu, A Y Tok, T J White, S Magdassi, I Abdulhalim, Y Long. Acs. Appl. Mater. Inter., 8(2016).
[21] C Z Fan, J Q Wang, S M Zhu, J N He, P Ding. J. Opt., 15(2013).
[22] Q Liu, S G Li, J S Li, C Dou, X Y Wang, G Y Wang, M Shi. J. Lightwave. Technol., 34, 2484(2016).
[23] F Ghasemi, S R Entezar, S Razi. Phys. Lett. A, 383, 2551(2019).
[24] F Fan, Y Hou, Z W Jiang, X H Wang, S J Chang. Appl. Opt., 51, 4589(2012).
[25] J R Liang, P Li, X L Song, L W Zhou. Appl. Phys. A, 123, 794(2017).
[26] R Wang, W Y Yang, S Gao, X J Ju, P F Zhu, B Li, Q Li. J. Mater. Chem. C, 7, 8185(2019).
[27] M Suzuki. Phys. Rev. B, 31, 2957(1985).
[28] M Lin, Z B Ouyang, J Xu, G X Qiu. Opt. Express, 17, 5861(2009).
[29] F Bréchet, J Marcou, D Pagnoux, P Roy. Opt. Fiber. Technol., 6, 181(2000).
[30] X D Xiao, H L Cheng, G P Dong, Y G Liu, L H Chen, L Miao, G Xu. CrystEngComm, 15, 1095(2013).
[31] H T Kim, B G Chae, D H Youn, S L Maeng, G Kim, K Y Kang, Y S Lim. New J. Phys., 6, 52(2004).
[32] B J Kim, Y W Lee, S Choi, J W Lim, S J Yun, H T Kim. Phys. Rev. B, 77(2008).
[33] P U Jepsen, B M Fischer, A Thoman, H Helm. Phys. Rev. B, 74(2006).
[34] T L Cocker, L V Titova, S Fourmaux, H C Bandulet, D Brassard, J C Kieffer, M A E Khakani, F A Hegmann. Appl. Phys. Lett., 97(2010).
[35] C Z Fan, J Q Wang, J N He, P Ding, E J Liang. Chin. Phys. B, 22(2013).
Get Citation
Copy Citation Text
Yuanlin Jia, Peiwen Ren, Chunzhen Fan. Thermal tunable one-dimensional photonic crystals containing phase change material[J]. Chinese Physics B, 2020, 29(10):
Received: May. 10, 2020
Accepted: --
Published Online: Apr. 21, 2021
The Author Email: Chunzhen Fan (chunzhen@zzu.edu.cn)