Of late, considerable attention has been paid to the transformation method,12.–3 because it provides an intuitive and direct way of designing interesting materials that can control electromagnetic,12.–3 acoustic,4 and elastic5 wave propagation at will. Many interesting devices, such as invisible cloaks,2,3 have been designed based on transformation method. This method has also been explored to improve optical information processing. Li et al.6 designed a Fourier lens that has wider input spatial frequency bandwidth by applying transformation optics to the conventional graded index (GRIN) lens. Their result, however, depended on an anisotropic material, which can usually be realized with local resonant mechanism or metamaterial technology,2 and is challenging for low-loss and broadband applications. For example, the lens constructed by Li et al.6 can only work in a particular wavelength range. Furthermore, this anisotropic material also limits the application of the lens to the fractional Fourier transform (FRFT),7 which is widely used in image processing and other information systems.8 To indicate this limit, we consider the two-dimensional (2D) device for transverse electric (TE) or transverse magnetic (TM) waves. A conventional GRIN lens can have the refractive index distribution as , where is the radial distance from the optical axis and and are the GRIN lens parameters.9 Under the paraxial approximation, if the lens thickness , the signal at the output facet becomes a Fourier transformed one of the input signal.9 According to the optically interpreted FRFT, the -order FRFT can be physically defined as the functional form of the output signal at from the input facet,7 which is located at the inner region of the lens. As the anisotropic unit block is constructed from different material layers,6 the system is also challenging for small-scale fabrication and thus the inner region of the lens is discrete to some degree. This discreteness will influence the signal detection precision in the inner region. Therefore, to promote the broadband and low-loss practical application of optical FRFT, where the traditional Fourier transform is viewed as a special case when , an isotropic material is in great demand for such a transformation lens.
To achieve the isotropic transformation media, one can use the numerical method based on inverse Laplace’s equation with sliding boundaries.10 For further application, it is better to get the transformation material in analytical form, such as conformal mapping in the complex plane.3 Here we propose a potentially more flexible alternative of this design, with which the space mapping can be obtained step by step directed by straightforward geometrical senses. To this end, the deformation view of the transformation optics is introduced, where the transformed material parameters can be expressed in a geometrical way11 as1) is established in the local principal system of the deformation. Instead of directly adjusting the explicit co-ordinate transformation expression, Eq. (1) suggests that one can adjust the principal stretches, which have very clear geometrical meaning, to manipulate the transformed material parameters.
Accordingly, our design should be divided into two steps, as shown in Fig. 1, by which the input facet of the conventional GRIN lens is compressed while the output facet remains the same.6 First, a rectangle ABCD is bent to a sector with the side length and the output facet size , both of which remain unchanged, where superscript indicates the first-step transformation. At each point, there is only rescaling in the direction during the deformation, and this stretch can be easily obtained by comparing the arc length with its original length , thus the principal stretches are1), this transformation will result in the anisotropic medium, as was shown in Ref. 6. Now, let us carry out the second-step transformation , in the polar co-ordinate system shown in Fig. 1(b) and 1(c), where the continuous function is to be determined. This transformation will lead to principal stretches 2) and (3) into Eq. (4), one has 5) as 1), will be , , , and . Thus, the isotropic refractive index for TE or TM waves is , where is the refractive index of the conventional rectangle-shaped GRIN lens.
To validate the transformed material parameters, a numerical simulation was performed using the software COMSOL Multiphysics, where we set and , , with , i.e., . The enhancement factor of the designed lens is the compression factor of the input boundary and can be obtained from Eqs. (5) and (6). In this example, . The simulation results are shown in Fig. 2, where the wavelength of the input TE wave is 0.15 μm, and the wave is impinging on the input facet with aperture width , therefore the rectangular aperture function is optically Fourier transformed to sinc function at the output facet. For the designed lens, we set , shown in Fig. 2(b); for the conventional GRIN lens, we set , shown in Fig. 2(a). Owing to paraxial approximation, the widths of the output signals that are well matched with the ideal ones are about 1.6 μm. It is clear that a narrower input signal in the designed lens can obtain the same effective output signal, shown in Fig. 2(e) as that in the conventional GRIN lens in Fig. 2(d) with a -times wider input signal; thus, a -times enhancement for the input bandwidth was achieved.6 Furthermore, we set also equal to 0.3 μm for the conventional GRIN lens in Fig. 2(c) and checked its output signal , where denotes spatial frequency and .9 It can be shown that in the conventional lens if , where the sinc function becomes zero as seen in Fig. 2(f), the corresponding in the designed lens is about , shown in Fig. 2(e), thus a greater amount (2.38 times) of spatial frequency is revealed in the latter.
The designed lens itself is an FRFT device because of its semigroup property.7 To obtain its order FRFT, the detection line is usually located at from the input facet. However, the FRFT result will be different from that of the conventional one with the same order , because this detection location is not the transformation counterpart of the original one. The points at line in Fig. 1(a) are transformed to points at line , not the line in Fig. 1(c), as shown in Eq. (6). Thus, according to the idea of transformation optics, to obtain the same FRFT results in the designed lens as those of the conventional one for the same order , the detection location must be moved from the line to . To validate the adjustment effects, FRFTs with , 0.5, and 0.75 were checked, respectively, in the same simulation environment as Fig. 2(a) and 2(b), and the results are shown in Fig. 3, where the compressed widths of the FRFT results in the designed lens have already been normalized for convenience of comparison. Only the results near the optical axis are considered according to the paraxial approximation. It is clear that the designed lens can provide almost the same FRFT results as that of the conventional lens for the same order. More demonstrations of the optical FRFT in conventional GRIN lens can be found in Ref. 12.
To conclude, keeping the principal stretches equal in the space mapping, an isotropic material can be designed for the Fourier lens that has an enhanced input spatial frequency bandwidth. The isotropic material, having important features like broadband and low loss, is also less discrete, and that helps the practical application of FRFT. More discussions on the deformation view of the transformation method can be found in Refs. 1314.–15. The clear geometrical sense provided by this view can help in optimizing more potential transformation media designs.
This work was supported by the National Natural Science Foundation of China (11172037) and Excellent Young Scholars Research Fund of Beijing Institute of Technology (2011YR0509).
A. GreenleafM. LassasG. Uhlmann, “On non-uniqueness for Calderon’s inverse problem,” Math. Res. Lett. 10(5–6), 685–693 (2003).1073-2780Google Scholar
J. HuZ. ChangG. K. Hu, “Approximate method for controlling solid elastic waves by transformation media,” Phys. Rev. B 84(20), 201101(R) (2011).PRBMDO1098-0121Google Scholar
D. MendlovicH. M. Ozaktas, “Fractional Fourier transforms and their optical implement: I,” J. Opt. Soc. Am. A 10(9), 1875–1881 (1993).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.10.001875Google Scholar
A. Yariv, Optical Electronics, 3rd ed., Holt Reinhart, New York (1985).Google Scholar
Z. Changet al., “Design method for quasi-isotropic transformation materials based on inverse Laplace's equation with sliding boundaries,” Opt. Express 18(6), 6089–6096 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.006089Google Scholar
H. M. OzaktasD. Mendlovic, “Fourier transforms of fractional order and their optical interpretation,” Opt. Commun. 101(3–4), 163–169 (1993).OPCOB80030-4018http://dx.doi.org/10.1016/0030-4018(93)90359-DGoogle Scholar
J. HuX. LiuG. K. Hu, “Constraint condition on transformation relation for generalized acoustics,” Wave Motion 50(2), 170–179 (2013).WAMOD90165-2125http://dx.doi.org/10.1016/j.wavemoti.2012.08.004Google Scholar
J. HuX. Y. Lu, “Determining the full transformation relations in the transformation method,” Appl. Phys. A 109(4), 971–]977 (2012).APAMFC0947-8396http://dx.doi.org/10.1007/s00339-012-7370-zGoogle Scholar