Abstract
In recent years, metasurface-based focusing elements have gradually become an indispensable type of terahertz lenses. However, the meta-lens often suffers from chromatic aberration due to the intrinsic dispersion of each element, especially in the broadband application scenarios. In this paper, we design and demonstrate a silicon-based achromatic meta-lens working from 0.6 to 1.0 THz, which is polarization insensitive because of the adopted symmetrical structures. The simulated focal length and the full width at half maximum (FWHM) of the foci at different frequencies prove the achromatic behavior of our meta-lens compared with the chromatic counterpart. We also show that the focus shift incongruence of our design originates from the transmission amplitude distribution of the meta-lens. This article not only provides an achromatic planar lens working at terahertz domain but also reveals the importance of the amplitude distribution in the achromatic metasurface design.
1. Introduction
In recent years, the terahertz frequency domain has ushered in extensive and in-depth development. Various terahertz techniques and systems have been proposed for sensing, imaging, and communication scenarios, showing promising application prospects [1]. In almost all terahertz systems, the lens is an important component that can directly determine the performance of the system [2, 3]. Traditional terahertz lenses manipulate the output wavefront through phase accumulation along the propagation path via thickness variation, which is bound to a bulky size and curved shape, damaging the compactness of the lens. The metasurface, a 2D version metamaterial composed of an array of artificial subwavelength structures with the ability to engineer the phase, amplitude, and polarization features, has provided a new paradigm for the design of terahertz lenses. From the early metal gradient index terahertz lens to the all-dielectric meta-lens, which has attracted wide attention in recent years, metasurface-based focusing elements have gradually become an indispensable type of terahertz lenses. Compared with traditional terahertz lenses, terahertz meta-lenses composed of graphene [4, 5], silicon [6, 7], metals, etc. [8–14] can control the amplitude and/or phase distributions of the terahertz wavefront at quite a high freedom level to achieve novel focusing functions [15, 16]. Moreover, most meta-lenses have a planar shape and wavelength order thickness with potential outstanding integration capability.
However, compared with the silicon-based traditional terahertz lens, the meta-lens still faces the challenge of large chromatic aberration introduced by the mode dispersion of each meta-atom. For example, as the most mature terahertz technology, the time-domain spectroscopy uses a broadband terahertz pulse to analyze the spectrum of various materials, which has an natural demand for achromatic lenses. Recently in the visible domain, a huge line of meta-lens research has been dedicated to achieving the dispersion engineering. For instance, Khorasaninejad et al. [17] carried out a series of scientific researches on the achromatic meta-lens with a reflection mode. Chen et al. [18, 19] designed several achromatic meta-lenses with anisotropic nanostructures in the visible domain. Shrestha et al. [20] proposed an achromatic meta-lens working from 1,200 to 1,650 nm. Wang et al. [21, 22] demonstrated an achromatic meta-lens through integrated resonances for imaging. Despite those inspiring advances in the visible light domain, terahertz metasurfaces are critically different from visible light metasurfaces in terms of raw materials and structure design. Only sporadic study of achromatic meta-lenses has been conducted in the terahertz domain; Cheng et al. [23] proposed an achromatic lens based on silicon metasurface with a high aspect ratio. Though the achromatic behavior was demonstrated, the efficiency of their terahertz meta-lens is somehow limited, and the dispersion engineering capability is polarization sensitive due to the adopted geometrical phase modulation method.
In this paper, we propose and investigate a silicon based meta-lens with a focal length of 11 mm and a diameter of 6.4 mm. Etched silicon pillars with 190 µm height are adopted as the unit cells to accumulate the required transmissive phase distribution. By changing the cross section of the pillars, alteration of the linear dispersion slop of each unit cell is achieved, which satisfies the requirement of an achromatic lens. The numerical investigation shows that our meta-lens has an achromatic bandwidth ranging from 0.6 to 1.0 THz that is polarization insensitive due to the symmetry of the cross section geometry. The FWHMs of the foci at different frequencies are smaller than those of the chromatic counterpart. We also analyze the deviation between the simulated focal length and the designed value by the Rayleigh-Sommerfeld integration. This study has provided an effective method for solving the chromatic aberration problem in the terahertz lens design and paved the way for the research on terahertz achromatic devices that require dispersion engineering.
2. Design and Unit Cell Construction
The designed meta-lens is illustrated by Figure 1A. The desired phase distribution of a meta-lens working for a monochromatic wave can be expressed aswhere ω is the wave frequency, c is the speed of light in vacuum, f is the focal length, and r and represent the positions of the unit cell and the reference cell on the meta-lens, respectively. represents the required optical path difference to realize a focusing effect. However, for an achromatic meta-lens working in the frequency range , should also be equivalent at different frequencies because the focal point is fixed at different frequencies. [24]. With as the minimal frequency in the design process, the required phase distribution for an achromatic meta-lens is
FIGURE 1
The unit cells for the achromatic meta-lens are thus expected to have a phase distribution similar to the chromatic lens designed for a monochromatic wave and a linear dispersed phase whose dispersion slope is determined by . To ensure the achromatic effect, the phase generated by the unit cells should be proportional to the frequency in the working frequency range.
The first step in the construction of the unit cells is to select the right raw material by using the commercial software CST Microwave StudioTM and applying the periodic boundary condition with normal incidence wave. In the visible light domain, the substrate and the resonate cavity are often made of different materials in order to utilize the dispersion feature of the Mie scattering of a closed cavity. However, in the terahertz band, due to the ease of the fabrication and the consideration of absorption, dielectric metasurfaces prefer to etch structures directly onto the semi-insulating silicon substrates. As shown in Figure 2A, when both the substrate and the structure are made of silicon, a larger variation range of the phase and dispersion can be obtained by subwavelength pillar structure, compared to the fused silica structure with silicon substrate. Conversely, when the unit cell is made of a silicon pillar with fused silica substrate, its phase and dispersion have nearly the same ranges as the all silicon unit cells’, while the larger absorption and material dispersion of silica may cause problems in the design process.
FIGURE 2
Considering the aspect ratio in the dry etching process, we set the height of the silicon pillars as 190 µm, and the fabrication accessibility of our selection is supported by our previous work [25]. The periodicity of the unit cell was set as 130 µm, and our meta-lens has 49 unit cells in total with a diameter of 6.4 mm. A focal length of 13.5 mm and the target frequency ranging from 0.6 to 1.0 THz were determined as well for the followed unit cell construction process. Inspired by the previous work of Shrestha et al [20], we appropriately changed the cross-section geometry of the silicon pillars while retaining the symmetry to hold the polarization-independent characteristics of the structures to engineer the phase and dispersion simultaneously. The specific cross section shapes are presented in Figure 1B, which are square ring, crisscross and hollow regular octagon, respectively. As depicted in Figures 2B,C, the three unit cells have different phases and amplitudes, while the phase has a quasi-linear dispersion and the transmission is relatively high in the target frequency range, indicating the silicon pillar structure proposed here has the capability to realize the desired phase range and liner dispersion. Besides, we found that adding a thin layer of pure gold on the dielectric silicon pillars could further endow the freedom level of the phase control. When the structure size was the same, the phase range of the metal-attached structure was significantly increased (Figure 3A). Hence we adopt this method in the construction of the reference unit cell. Next, by changing the cross section shape and scales in the simulation, we constructed a unit cell base and drew its phase-dispersion library in Figure 3B. We found that for all the investigated unit cells, the phase decreased linearly as the frequency increased from 0.6 to 1.0 THz, and the dispersion slope differs according to the specific cross section shape and geometric parameters. Finally, taking the lens edge as the reference position in the design, we calculated the phase distribution at and the dispersion slop for each position r. Then, different shapes of unit cells whose phase and dispersion slop satisfy the required at as well as the dispersion slop were chosen to form the metalens. The whole sketch of our metalens is depicted in Figure 4D. The adopted cross section shape and the geometric parameters are listed in Table 1. It can be seen in Figure 3C that the phase and dispersion of the selected unit cells fit the theoretical design quite well.
FIGURE 3
FIGURE 4
TABLE 1
| r/µm | Outer diameter a/µm | Inner diameter b/µm | r/µm | Outer diameter a/µm | Inner diameter b/µm |
|---|---|---|---|---|---|
| 0 | 130 | 61 | 1,690 | 130 | 73 |
| 130 | 130 | 61 | 1,820 | 130 | 78 |
| 260 | 130 | 62 | 1,950 | 130 | 84 |
| 390 | 130 | 63 | 2,080 | 130 | 89 |
| 520 | 130 | 66 | 2,210 | 130 | 93 |
| 650 | 130 | 62 | 2,340 | 130 | 97 |
| 780 | 130 | 58 | 2,470 | 130 | 102 |
| 910 | 130 | 56 | 2,600 | 130 | 106 |
| 1,040 | 130 | 42 | 2,730 | 45 | 31 |
| 1,170 | 130 | 49 | 2,860 | 40 | 27 |
| 1,300 | 130 | 56 | 2,990 | 33 | 20 |
| 1,430 | 130 | 62 | 3,120 | 60 | 24 |
| 1,560 | 130 | 68 |
The parameters of the selected structures and the corresponding position on the meta-lens.
3. Numerical Demonstration and Analysis
The performance of our achromatic meta-lens was also evaluated using the CST Microwave Studio. A full meta-lens was constructed by the selected unit cells with the center of the meta-lens set at the origin of the simulation environment, the optical axis of which was set as the z-axis. With the open boundary conditions and a plane wave incidence with linear polarization, we obtained the electric field distributions of the achromatic focusing effects by setting three monitors at 0.6, 0.8, and 1.0 THz, respectively, and the cut of the distributions are shown in Figure 4. To inspect the achromatic feature of our meta-lens, a control meta-lens designed only for 0.6 THz and composed of pure crisscross pillars with other conditions remaining the same was investigated as well. It is quite clear that both the achromatic meta-lens and the control show focusing capabilities. By analyzing the maximal value of the electric field intensity, we verified that the foci of the achromatic meta-lens at 0.6, 0.8 and 1.0 THz lay at z = 11.20, 11.80 and 11.85 mm, respectively, which only has an maximal offset of 650 µm (Figure 4A). In contrast, the foci of the control meta-lens at the corresponding frequencies lay at z = 11.20, 11.80 and 13.30 mm, respectively, which has a much larger focal shifts of 2 mm as Figure 4B showed. Figure 4C depicts the evolution of the focal length as frequencies changes for the achromatic meta-lens and the control. The focal length of the achromatic meta-lens varies less than that of the control. We also found that the maximal electric field intensity point deviated 6 mm from the foci at 0.8 and 1.0 THz. Our achromatic meta-lens design is successful since it demonstrated a clearer focusing behavior and a more stable focal lengths over the whole interested frequency band than the control.
Figure 5A depicts the intensity distributions of our achromatic meta-lens at the focal planes () for the monitored frequencies, respectively. It can be seen that the focus spot shows a perfect circle shape and a typical Gaussian intensity distribution. Fringes around the spot are a common phenomenon in diffraction focusing devices caused by the limited unit cell periodicity. The horizontal cuts along the x-axis for the achromatic meta-lens and the control meta-lens are shown in Figures 5B,C. The FWHM of the achromatic meta-lens at 0.6, 0.8 and 1.0 THz are 1.628, 1.277, and 1.063 mm, respectively, yet the FWHM of the control at 0.6, 0.8, and 1.0 THz are 1.621, 1.382 and 1.892 mm, respectively. The evident reduction of the FWHM is a typical feature for achromatic lens, indicating a higher spatial resolution in broadband imaging or microscopy applications which is a well-known advantage of achromatic lens. The polarization dependence of the achromatic meta-lens was manifested by cut of the electric field intensities as shown in Figure 6 It is clear that the distributions of the electric field are almost the same, with only small variances compared with the electric field distribution shown in Figure 5, which confirms the symmetry of each unit cell.
FIGURE 5
FIGURE 6
4. Discussion
Although the simulation results have verified the feasibility of the designed achromatic meta-lens, the numerically obtained focal length is around 11.2 mm, which deviates 2.3 mm from the designed focal length of 13.5 mm. In addition, the focal lengths at 0.6 and 1.0 THz also differ by 650 µm. We believe this focal shift and the residual chromatic aberration are due to two reasons. The first one is the hexahedral mesh method in the simulation, which divides the space along Cartesian coordinates and will thus cause a calculation error because of the mismatch between the mesh grid and the tilted edge of the octagon pillars. The second reason is that focal shift is a prevalent phenomenon of the lens design [26], and we analyzed it from the perspective of the reversible light path. By assuming a point terahertz source located at , it is easy to obtain the resulted amplitude and the phase distributions on the meta-lens. Though the phase distribution has been well satisfied in the whole interested bandwidth, the required amplitude variance, which equals , was not considered in our numerical investigation. To investigate the influence of this amplitude distribution on the focal shift, we numerically testified an achromatic meta-lens with and without the amplitude distribution by the Rayleigh-Sommerfeld diffraction formula, respectively [27]. As shown in Figure 7, without the amplitude distribution, the calculated focal lengths of the achromatic meta-lens are in accordance with the simulated results, verifying our calculation program. While by adding the right amplitude distribution to the constructed achromatic meta-lens, the foci of this ‘ideal’ achromatic meta-lens at 0.6, 0.8, and 1.0 THz now all move to , which is highly consistent with our design and shows even smaller chromatic features. This comparison inspires the researchers to pay more attention to the amplitude distribution in the design of achromatic meta-lenses, especially for meta-lenses with a large focal length.
FIGURE 7
5. Conclusion
In summary, we designed and numerically investigated an achromatic meta-lens composed of subwavelength silicon pillars working at terahertz frequencies. The achromatic behavior of the meta-lens was verified by the numerical simulation, which showed that its focal length changed by 650 µm when frequency ranged from 0.6 to 1.0 THz. The variation of the focal length is merely . The FWHM of the achromatic meta-lens foci ranges from 1.063 to 1.628 mm, indicating a higher spatial resolution compared with that of the chromatic counterpart. The incongruence of the focal lengths with the design was well explained by adding the correct amplitude distribution in the Rayleigh-Sommerfeld diffraction calculation. The achromatic meta-lens proposed here not only provides an achromatic planar lens working at terahertz domain but also reveals the importance of the amplitude distribution in achromatic metasurface design.
Funding
National Natural Science Foundation of China (61307125, 61420106006, 61422509, 61575141, 61622505, 61975143); the US National Science Foundation (ECCS-1232081).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
YG and JG discussed the design principle. YG and RJ completed the unit cell construction. YG performed the CST simulations, the matlab calculations, and the data analysis. YG wrote the paper, and JG modified the writing. All the authors commented on the paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2020.606693/full#supplementary-material.
References
1.
ZhangXCXuJ. Introduction to THz wave photonics, Vol. 29. New York, NY: Springer (2010).
2.
CooperKBDenglerRJLlombartNThomasBChattopadhyayGSiegelPH. Thz imaging radar for standoff personnel screening. IEEE Transactions on Terahertz Science and Technology (2011) 1:169–82. 10.1109/TTHZ.2011.2159556
3.
TrichopoulosGCMosbackerHLBurdetteDSertelK. A broadband focal plane array camera for real-time thz imaging applications. IEEE Trans Antenn Propag (2013) 61:1733–40. 10.1109/TAP.2013.2242829
4.
LiuWHuBHuangZGuanHLiHWangXet alGraphene-enabled electrically controlled terahertz meta-lens. Photon Res (2018) 6:703–8. 10.1364/PRJ.6.000703
5.
HuangZHuBLiuWLiuJWangY. Dynamical tuning of terahertz meta-lens assisted by graphene. JOSA B (2017) 34:1848–54. 10.1364/JOSAB.34.001848
6.
YuQGuJYangQZhangYLiYTianZet alAll-dielectric meta-lens designed for photoconductive terahertz antennas. IEEE Photonics J (2017) 9:1–9. 10.1109/JPHOT.2017.2721503
7.
JiaDTianYMaWGongXYuJZhaoGet alTransmissive terahertz metalens with full phase control based on a dielectric metasurface. Optic Lett (2017) 42:4494–7. 10.1364/OL.42.004494
8.
YangQGuJWangDZhangXTianZOuyangCet alEfficient flat metasurface lens for terahertz imaging. Optic Express (2014) 22:25931–9. 10.1364/OE.22.025931
9.
YangQGuJXuYZhangXLiYOuyangCet alBroadband and robust metalens with nonlinear phase profiles for efficient terahertz wave control. Advanced Optical Materials (2017) 5:1601084. 10.1002/adom.201601084
10.
LuoJYuHSongMZhangZ. Highly efficient wavefront manipulation in terahertz based on plasmonic gradient metasurfaces. Optic Lett (2014) 39:2229–31. 10.1364/OL.39.002229
11.
HuDMorenoGWangXHeJChahadihAXieZet alDispersion characteristic of ultrathin terahertz planar lenses based on metasurface. Optic Commun (2014) 322:164–8. 10.1016/j.optcom.2014.02.017
12.
HuDWangXFengSYeJSunWKanQet alUltrathin terahertz planar elements. Adv Opt Mater (2013) 1:186–91. 10.1002/adom.201200044
13.
JiangXYYeJSHeJWWangXKHuDFengSFet alAn ultrathin terahertz lens with axial long focal depth based on metasurfaces. Optic Express (2013) 21:30030–8. 10.1364/OE.21.030030
14.
WangSWangXKanQYeJFengSSunWet alSpin-selected focusing and imaging based on metasurface lens. Optic Express (2015) 23:26434–41. 10.1364/OE.23.026434
15.
ChangCCHeadlandDAbbottDWithayachumnankulWChenHT. Demonstration of a highly efficient terahertz flat lens employing tri-layer metasurfaces. Opt Lett (2017) 42:1867–70. 10.1364/OL.42.001867
16.
JiangXChenHLiZYuanHCaoLLuoZet alAll-dielectric metalens for terahertz wave imaging. Optic Express (2018) 26:14132–42. 10.1364/OE.26.014132
17.
KhorasaninejadMShiZZhuAYChenWTSanjeevVZaidiAet alAchromatic metalens over 60 nm bandwidth in the visible and metalens with reverse chromatic dispersion. Nano Letters (2017) 17:1819–24. 10.1021/acs.nanolett.6b05137
18.
ChenWTZhuAYSanjeevVKhorasaninejadMShiZLeeEet alA broadband achromatic metalens for focusing and imaging in the visible. Nat Nanotechnol (2018) 13:220–6. 10.1038/s41565-017-0034-6
19.
ChenWTZhuAYSislerJBharwaniZCapassoF. A broadband achromatic polarization-insensitive metalens consisting of anisotropic nanostructures. Nat Commun (2019) 10:1–7. 10.1038/s41467-019-08305-y
20.
ShresthaSOvervigACLuMSteinAYuN. Broadband achromatic dielectric metalenses. Light Sci Appl (2018) 7:1–11. 10.1038/s41377-018-0078-x
21.
WangSWuPCSuVCLaiYCChuCHChenJWet alBroadband achromatic optical metasurface devices. Nat Commun (2017) 8:187. 10.1038/s41467-017-00166-7
22.
WangSWuPCSuVCLaiYCChenMKKuoHYet alA broadband achromatic metalens in the visible. Nat Nanotechnol (2018) 13:227–32. 10.1038/s41565-017-0052-4
23.
ChengQMaMYuDShenZXieJWangJet alBroadband achromatic metalens in terahertz regime. Sci Bull (2019) 64:1525–31. 10.1016/j.scib.2019.08.004
24.
AietaFKatsMGenevetPCapassoF. Applied optics. Multiwavelength achromatic metasurfaces by dispersive phase compensation. Science (2015) 347:1342–5. 10.1126/science.aaa2494
25.
XuYLiQZhangXWeiMXuQWangQet alSpin-decoupled multifunctional metasurface for asymmetric polarization generation. ACS Photonics (2019) 6:2933–41. 10.1021/acsphotonics.9b01047
26.
LiYWolfE. Focal shift in focused truncated Gaussian beams. Optic Commun (1982) 42:151–6. 10.1016/0030-4018(82)90128-6
27.
LuckeRL. Rayleigh–sommerfeld diffraction and Poisson’s spot. Eur J Phys (2006) 27:193–204. 10.1088/0143-0807/27/2/002
Summary
Keywords
achromatic, meta-lens, polarization independent, terahertz, metasurfaces
Citation
Gao Y, Gu J, Jia R, Tian Z, Ouyang C, Han J and Zhang W (2020) Polarization Independent Achromatic Meta-Lens Designed for the Terahertz Domain. Front. Phys. 8:606693. doi: 10.3389/fphy.2020.606693
Received
15 September 2020
Accepted
16 November 2020
Published
21 December 2020
Volume
8 - 2020
Edited by
Longqing Cong, University of Pennsylvania, United States
Reviewed by
Lin Chen, University of Shanghai for Science and Technology, China
M. Ren, Nankai University, China
Updates
Copyright
© 2020 Gao, Gu, Jia, Tian, Ouyang, Han and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianqiang Gu, gjq@tju.edu.cn; Weili Zhang, weili.zhang@okstate.edu
This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.