## Abstract

In this work, we propose that one-dimensional ultratransparent dielectric photonic crystals with wide-angle impedance matching and shifted elliptical equal frequency contours are promising candidate materials for illusion optics. The shift of the equal frequency contour does not affect the refractive behaviors, but enables a new degree of freedom in phase modulation. With such ultratransparent photonic crystals, we demonstrate some applications in illusion optics, including creating illusions of a different-sized scatterer and a shifted source with opposite phase. Such ultratransparent dielectric photonic crystals may establish a feasible platform for illusion optics devices at optical frequencies.

© 2017 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Illusionary effects like mirages have fascinated people for centuries and have inspired myths, novels, and films. With the development of metamaterials [1–4] and transformation optics (TO) [5–26], more fascinating illusionary effects have been proposed, including invisibility cloaking [12], transforming the images of one object into another [13, 14], opening virtual holes in a wall [13, 15], superscattering and gateway [16–22], super absorbers [24], illusion of shifted sources [25], etc. The creation of such intriguing illusionary effects requires materials with unusual properties such as negative index or anisotropy in both permittivity and permeability, which usually can only be realized by metamaterials [1–4] with the units in the deep subwavelength scale. Due to the difficulty in realizing metamaterials with the ideal TO parameters, “reduced parameters” that sacrifice the property of omnidirectional impedance matching have been generally adopted [7, 19–21].

Very recently, by using photonic crystals (PhCs) with the units at the wavelength scale, ultratransparency media with omnidirectional impedance matching and the ability of forming aberration-free virtual images have been realized, which provides another candidate materials for TO devices [27]. Interestingly, the equal frequency contours (EFCs) of such ultratransparent PhCs are not centered the Brillouin zone center, but shifted towards the $X$ point. In other words, such ultratransparent media are nonlocal or spatially dispersive, beyond the description of local TO media. Recently, we find that such ultratransparency effect can be achieved for both polarizations in a broad range of frequencies by using one-dimensional (1D) dielectric PhCs [27, 28] for a relatively narrow range of incident angles.

In this work, we propose and demonstrate the application of such 1D ultratransparent PhCs with wide-angle impedance matching and shifted elliptical EFCs in TO, and more specifically, in the creation of optical illusions, which is also known as illusion optics. We demonstrate two examples: one creates the illusion of a different-sized scatterer, and the other creates the illusion of a shifted source with an opposite phase. We show that the shift of the EFCs of the ultratransparent PhCs does not affect the refractive behaviors, but provides a new degree of freedom for additional phase modulation. With such 1D ultratransparent PhCs, TO and illusion optics become more feasible in the optical frequency regime.

## 2. TO media and ultratransparent PhCs

To begin with, we consider a TO medium obtained by stretching the coordinate along the $x$ direction in the background medium of air. For transverse electric (TE) polarization with electric fields polarized in the $z$ direction, the parameters of the TO medium are written as [23],

where $s$ is the stretching ratio. ${\epsilon}_{z}$, ${\mu}_{x}$ and ${\mu}_{y}$ are, respectively, the $z$ component of relative permittivity tensor, $x$ and $y$ components of relative permeability tensor.Thus, the dispersion of the TO medium, i.e. $\frac{{k}_{x}^{2}}{{\mu}_{y}}+\frac{{k}_{y}^{2}}{{\mu}_{x}}={\epsilon}_{z}{k}_{0}^{2}$, can be rewritten as ${s}^{2}{k}_{x}^{2}+{k}_{y}^{2}={k}_{0}^{2}$, indicating that the EFC is an ellipse with the same height as the EFC of air in the ${k}_{y}$ direction, as shown in Fig. 1(a). Here, ${k}_{x}$ and ${k}_{y}$ are the $x$ and $y$ components of wave vectors in the TO medium, respectively. ${k}_{0}$ is the wave number in air. As an example, the stretching ratio is set to be larger than unity, i.e., $s>1$. Thus, the EFC of the TO medium (blue lines) is slimmer than the EFC of air (gray lines) in Fig. 1(a). Thus, the coordinate mesh is looser in the $x$ direction in the TO medium region, as illustrated in Fig. 1(a).

If a “shifted” EFC of the same shape as that of the TO medium can be obtained in a 1D ultratransparent PhC [see Fig. 1(b)], then the refraction of a slab of such 1D ultratransparent PhC would be the same as the refraction of a slab of TO media, because the shape of the EFC determines the refractive behaviors. The wide-angle impedance matching of ultratransparent PhCs ensures almost zero reflection and total transmission. Therefore, the functionality of such a slab of 1D ultratransparent PhC can be applied to replace a slab of the TO media. The only difference is that there may exist additional phases in the ultratransparent PhCs due to the shift of the EFC.

Through a trial-and-error iterative optimal algorithm, wide-angle impedance-matched 1D ultratransparent PhCs with the same shaped EFCs as that of the TO media can be obtained [28]. The design of such 1D ultratransparent PhC turns out to be surprisingly easy, due to the large degrees of freedoms in engineering the nonlocal effective parameters in PhCs [27]. In Fig. 2, we show a specific example of 1D ultratransparent PhCs. The designed PhC is composed of symmetric unit cells, which consist of dielectric components I and II with a lattice constant of $a$ [see Fig. 2(a)]. The relative permittivities $\epsilon $ and the thicknesses $d$ of the two components are presented in Fig. 2(a). In Fig. 2(b), by using software COMSOL Multiphysics based on finite-element methods, we plot the EFC of the PhC (red solid lines) at the normalized frequency $fa/c=0.402$ for TE polarization. $f$ and $c$ are the eigen-frequency and light velocity in vacuum, respectively. Evidently, the EFC can be approximately regarded as a part of an ellipse located at the *X* point. Compared with the EFC of the TO medium with $\left\{{\epsilon}_{z},{\mu}_{x},{\mu}_{y}\right\}=\left\{0.5,2,0.5\right\}$ at the same frequency [blue dashed lines in Fig. 2(b)], we see that the EFC of the PhC has the same height in the ${k}_{y}$ direction, which satisfies the requirement proposed above. In addition, in Fig. 2(c), we plot the transmittance through the 1D PhC slab composed of $n$ unit cells as the function of the incident angle at $fa/c=0.402$. We observe almost perfect transmission under the incident angle of 70° irrespective of the number of unit cells, indicating a wide-angle impedance matching effect. Thus, such a 1D ultratransparent PhC can approximately work as the TO medium with $\left\{{\epsilon}_{z},{\mu}_{x},{\mu}_{y}\right\}=\left\{0.5,2,0.5\right\}$, whose incident angle-dependent transmittance is denoted by the purple lines in Fig. 2(c).

To further verify the above results and show the effect of the shift of the EFC, with the software COMSOL Multiphysics, we simulate the wave propagation through the TO medium slab with a thickness of $5a$ (upper inset) and a PhC slab with 5 unit cells (lower inset) in Fig. 2(d). The waves are incident from the left port boundary, and the upper and lower boundaries are set to be periodic boundary condition in the simulation setup. Simulation results show perfect transmission under the incident angle of 45°. Interestingly, we notice a $\pi $ phase difference in the transmission waves, which actually is caused by the shift of the PhC’s EFC. As the center of the PhC’s EFC is located at the $X$ point, an addition $\pi $ phase difference will be imposed on waves when passing through each unit cell. As a consequence, when the total number of the unit cells in the PhC slab is odd, a $\pi $ phase difference will be seen when compared with the TO medium. On the other hand, if the total number is even, there will be no phase difference, which is confirmed in the simulation results in Fig. 2(e). In Fig. 2(e), the thickness of the TO medium slab is $6a$ (upper inset), and the number of unit cells of the PhC slab is 6 (lower inset). Comparing the transmitted waves under the incident angle of 45°, we can see identical phases, in contrary to the case with an odd number of unit cells shown in Fig. 2(d).

## 3. Illusion devices via ultratransparent PhCs

Here, with the 1D ultratransparent PhCs, we can realize applications in illusion optics. As an example, we consider a cylindrical model in Fig. 3. Here, we start with an air cylinder with a radius of ${R}_{1}$, and assume a virtual core with a radius of ${R}_{2}$, as illustrated in Fig. 3(a). Then, we change the radius of the core to be ${R}_{3}$ through transforming the coordinate in the radial direction based on the method of TO, as illustrated in Fig. 3(b). Thus, the parameters of the shell will be gradient functions of the radius [8, 17], that is, ${\mu}_{r}=1/{\mu}_{\theta}=1+\frac{{R}_{1}}{r}\frac{{R}_{2}-{R}_{3}}{{R}_{1}-{R}_{2}}$ and ${\epsilon}_{z}=\frac{{R}_{1}-{R}_{2}}{r}\frac{r\left({R}_{1}-{R}_{2}\right)+{R}_{1}\left({R}_{2}-{R}_{3}\right)}{{\left({R}_{1}-{R}_{3}\right)}^{2}}$ with ${R}_{1}\ge r\ge {R}_{3}$. And the parameters of the core are ${\mu}_{r}={\mu}_{\theta}=1$ and ${\epsilon}_{z}={\left({R}_{2}/{R}_{3}\right)}^{2}$. Next, we discretize the gradient shell into four uniform layers, as denoted by A, B, C and D in Fig. 3(c). The parameters of the four layers are determined roughly by the arithmetic mean of the graded parameters, as presented in Fig. 3(e). Finally, we design four corresponding 1D ultratransparent PhCs to replace the four uniform layers, as illustrated in Fig. 3(d). The relative permittivities and thickness of the dielectric components I and II of the four 1D PhCs are presented in Fig. 3(e). All the PhCs have the same lattice constant of $a$. The numbers of unit cells in layers A, B, C and D are 4, 2, 2 and 1, respectively. The working frequency is chosen as $fa/c=0.3$.

In the first example, we set the core to be a perfect electric conductor (PEC) to achieve illusionary scattering effect. For comparison, in Fig. 4(a), we first simulate the wave scattering by a PEC cylinder with a radius of $R=4a$ under the illumination of a TE polarized Gaussian beam in the air background. Then, in Fig. 4(b), we coat the PEC cylinder with the ideally gradient TO medium shell, as that in Fig. 3(b). The simulation results show that the scattering pattern is changed into that of a relative larger cylinder, as shown in Fig. 4(d). In Fig. 4(c), the 1D ultratransparent PhCs are used to replace the TO media, showing almost the same scattering pattern. As a result, any observer will see an illusion of a larger cylinder from scattered waves.

Such an illusionary effect can be understood from the theory of TO. Since the TO medium shell is obtained by compressing an air shell in the radial direction by varying the inner radius from $6.5a$ to $4a$, the model consisting of a PEC core with $R=4a$ and the TO medium shell is equivalent to the model consisting of a PEC core with $R=6.5a$ and an air shell. As the background medium is air, the wave scattering in Figs. 4(b) and 4(c) will be the same as the wave scattering by a PEC cylinder with $R=6.5a$, as confirmed by the simulation results in Fig. 4(d).

In the second example, we place a TE polarized monopolar point source inside the core to show the source illusion effect. The relative permittivity of the dielectric core is 2.64. In Fig. 5(a), the point source has an offset distance ${d}^{\text{'}}$ ($=1.6a$) from the center of the core, and the shell is composed of the 1D ultratransparent PhCs proposed above. From the field distribution, we see relatively standard cylindrical radiated waves in air due to the wide-angle impedance matching effect of the PhC shell. Moreover, one interesting phenomenon is that by back-tracking the radiation waves, the point source should have an offset distance of ${d}^{\text{'}\text{'}}=2.6a$, as displayed in Fig. 5(b). This means that the PhC shell enables a real source to resemble an illusion source off its original position. By using the TO theory, the relationship of the offset distance of the real and illusion source can be derived as,

Equation (2) reveals that the offset distance of the illusion source actually can be further controlled by varying the inner radius of the PhC shell. Another interesting observation is that the phase of the illusion source is opposite to that of the real source. This is induced by a total odd number of the unit cells of the ultratransparent PhCs. We note that such a source illusion effect has been demonstrated in [25] by using anisotropic negative-index materials. Interestingly, here only dielectric materials are required, which largely simplifies the design and fabrication.

Finally, it is worth noting the performances of the illusion devices consisting of ultratransparent PhCs can be further improved by making more discretized layers to achieve quasi-continuous parameters. The number of discretized layers is dependent on the size of the device. In fact, the above results reveal that four discretized layers are sufficient to obtain relatively good performances for most illusions devices.

## 4. Conclusions

In summary, we have demonstrated that 1D ultratransparent PhCs with shifted elliptical EFC and wide-angle impedance matching can approximately work as the TO media to realize optical illusionary effects. The shift of the EFC does not affect the refractive behaviors, but provide a new degree of freedom in phase modulation. Ultratransparent PhCs may provide a feasible platform for realization of more illusionary effects, especially in optical frequencies. We believe that there are diverse experimental methods to realize such illusionary effects. In the microwave frequency regime, a single material index with various filling ratios of holes is probably the simplest approach. In the optical frequency regime, the simplest method may be the thin-film fabrication, where multiple materials can be deposited alternatively to form the multiple film structure.

## Funding

National Natural Science Foundation of China (No. 11374224, 11574226, 11704271); Natural Science Foundation of Jiangsu Province (No. BK20170326); Natural Science Foundation for Colleges and Universities in Jiangsu Province of China (No. 17KJB140019); Jiangsu Planned Projects for Postdoctoral Research Funds (1701181B); Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

## References and links

**1. **J. B. Pendry, A. J. Holden, W. J. Stewart, and I. Youngs, “Extremely low frequency plasmons in metallic mesostructures,” Phys. Rev. Lett. **76**(25), 4773–4776 (1996). [PubMed]

**2. **J. B. Pendry, A. J. Holden, D. J. Robbins, and W. J. Stewart, “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE T. Microw. Theory **47**, 2075–2084 (1999).

**3. **D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett. **84**(18), 4184–4187 (2000). [PubMed]

**4. **R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science **292**(5514), 77–79 (2001). [PubMed]

**5. **J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science **312**(5781), 1780–1782 (2006). [PubMed]

**6. **U. Leonhardt, “Optical conformal mapping,” Science **312**(5781), 1777–1780 (2006). [PubMed]

**7. **D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial electromagnetic cloak at microwave frequencies,” Science **314**(5801), 977–980 (2006). [PubMed]

**8. **M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, D. R. Smith, and J. B. Pendry, “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s equations,” Photon. Nanostructures **6**, 87–95 (2008).

**9. **J. Li and J. B. Pendry, “Hiding under the carpet: A new strategy for cloaking,” Phys. Rev. Lett. **101**(20), 203901 (2008). [PubMed]

**10. **X. Ni, Z. J. Wong, M. Mrejen, Y. Wang, and X. Zhang, “An ultrathin invisibility skin cloak for visible light,” Science **349**(6254), 1310–1314 (2015). [PubMed]

**11. **Y. Lai, J. Ng, H. Chen, Z. Zhang, and C. T. Chan, “Illusion optics,” Front. Phys. China **5**, 308–318 (2010).

**12. **Y. Lai, H. Chen, Z. Q. Zhang, and C. T. Chan, “Complementary media invisibility cloak that cloaks objects at a distance outside the cloaking shell,” Phys. Rev. Lett. **102**(9), 093901 (2009). [PubMed]

**13. **Y. Lai, J. Ng, H. Chen, D. Han, J. Xiao, Z. Q. Zhang, and C. T. Chan, “Illusion optics: The optical transformation of an object into another object,” Phys. Rev. Lett. **102**(25), 253902 (2009). [PubMed]

**14. **J. Pendry, “All smoke and metamaterials,” Nature **460**, 579–580 (2009).

**15. **C. Shen, J. Xu, N. X. Fang, and Y. Jing, “Anisotropic complementary acoustic metamaterial for canceling out aberrating layers,” Phys. Rev. X **4**, 041033 (2014).

**16. **T. Yang, H. Chen, X. Luo, and H. Ma, “Superscatterer: enhancement of scattering with complementary media,” Opt. Express **16**(22), 18545–18550 (2008). [PubMed]

**17. **H. Chen, X. Zhang, X. Luo, H. Ma, and C. T. Chan, “Reshaping the perfect electrical conductor cylinder arbitrarily,” New J. Phys. **10**, 113016 (2008).

**18. **X. Luo, T. Yang, Y. Gu, H. Chen, and H. Ma, “Conceal an entrance by means of superscatterer,” Appl. Phys. Lett. **94**, 223513 (2009).

**19. **C. Li, X. Meng, X. Liu, F. Li, G. Fang, H. Chen, and C. T. Chan, “Experimental realization of a circuit-based broadband illusion-optics analogue,” Phys. Rev. Lett. **105**(23), 233906 (2010). [PubMed]

**20. **H. Xu, G. Wang, K. Ma, and T. J. Cui, “Superscatterer illusions without using complementary media,” Adv. Opt. Mater. **2**, 572–580 (2014).

**21. **G. D. Bai, F. Yang, W. X. Jiang, Z. L. Mei, and T. J. Cui, “Realization of a broadband electromagnetic gateway at microwave frequencies,” Appl. Phys. Lett. **107**, 153503 (2015).

**22. **C. Liu, P. Bai, and Y. Lai, “Sound-impenetrable holes in water based on acoustic complementary medium,” EPL **115**, 58002 (2016).

**23. **Y. Lai, H. Zheng, Z. Zhang, and C. T. Chan, “Manipulating sources using transformation optics with ‘folded geometry’,” J. Opt. **13**, 024009 (2011).

**24. **J. Ng, H. Chen, and C. T. Chan, “Metamaterial frequency-selective superabsorber,” Opt. Lett. **34**(5), 644–646 (2009). [PubMed]

**25. **J. Yi, P. Tichit, S. N. Burokur, and A. de Lustrac, “Illusion optics: Optically transforming the nature and the location of electromagnetic emissions,” J. Appl. Phys. **117**, 084903 (2015).

**26. **B. Zheng, H. A. Madni, R. Hao, X. Zhang, X. Liu, E. Li, and H. Chen, “Concealing arbitrary objects remotely with multi-folded transformation optics,” Light Sci. Appl. **5**, e16177 (2016).

**27. **J. Luo, Y. Yang, Z. Yao, W. Lu, B. Hou, Z. H. Hang, C. T. Chan, and Y. Lai, “Ultratransparent media and transformation optics with shifted spatial dispersions,” Phys. Rev. Lett. **117**(22), 223901 (2016). [PubMed]

**28. **Z. Yao, J. Luo, and Y. Lai, “Photonic crystals with broadband, wide-angle, and polarization-insensitive transparency,” Opt. Lett. **41**(21), 5106–5109 (2016). [PubMed]