• Photonics Research
  • Vol. 8, Issue 6, 830 (2020)
Bingshuang Yao1, Xiaofei Zang1、2、*, Zhen Li1, Lin Chen1、2, Jingya Xie1、2, Yiming Zhu1、2、3, and Songlin Zhuang1
Author Affiliations
  • 1Terahertz Technology Innovation Research Institute, Terahertz Spectrum and Imaging Technology Cooperative Innovation Center, Shanghai Key Laboratory of Modern Optical System, University of Shanghai for Science and Technology, Shanghai 200093, China
  • 2Shanghai Institute of Intelligent Science and Technology, Tongji University, Shanghai 200092, China
  • 3e-mail: ymzhu@usst.edu.cn
  • show less
    DOI: 10.1364/PRJ.387672 Cite this Article Set citation alerts
    Bingshuang Yao, Xiaofei Zang, Zhen Li, Lin Chen, Jingya Xie, Yiming Zhu, Songlin Zhuang. Dual-layered metasurfaces for asymmetric focusing[J]. Photonics Research, 2020, 8(6): 830 Copy Citation Text show less

    Abstract

    Asymmetric transmission, defined as the difference between the forward and backward transmission, enables a plethora of applications for on-chip integration and telecommunications. However, the traditional method for asymmetric transmission is to control the propagation direction of the waves, hindering further applications. Metasurfaces, a kind of two-dimensional metamaterials, have shown an unprecedented ability to manipulate the propagation direction, phase, and polarization of electromagnetic waves. Here we propose and experimentally demonstrate a metasurface-based directional device consisting of a geometric metasurface with spatially rotated microrods and metallic gratings, which can simultaneously control the phase, polarization, and propagation direction of waves, resulting in asymmetric focusing in the terahertz region. These dual-layered metasurfaces for asymmetric focusing can work in a wide bandwidth ranging from 0.6 to 1.1 THz. The flexible and robust approach for designing broadband asymmetric focusing may open a new avenue for compact devices with potential applications in encryption, information processing, and communication.

    1. INTRODUCTION

    Asymmetric transmission is defined as the difference in the total transmission between the forward and backward directions, exhibiting great potential applications such as directionally sensitive beam splitting [1,2], multiplexing [3], communications [4], and optical interconnection [5]. A typical method for asymmetric transmission is to break the time-reversal symmetry using an external magnetic field [6], time-varying component [7,8], or nonlinear materials [9], leading to a topological-protected edge mode (nonreciprocal transmission). Although this nonreciprocal edge mode is immune to backscattering, it suffers from complicated materials, external biases, and bulk configurations. Based on the conventional passive and linear materials (chiral metamaterials [1014], asymmetric metallic gratings [15], traditional polarization devices [16], and photonic crystals [17]), the reciprocal asymmetric transmission is realized by manipulating the propagation direction of electromagnetic (EM) waves.

    Metasurfaces, two-dimensional counterparts of metamaterials, can tailor the propagation direction, phase, and polarization of EM waves, providing a platform to manipulate asymmetric transmission in multiple degrees of freedom, i.e., polarization and phase modulation. Benefiting from the superior capability in the local manipulation of the wavefront of EM waves, a plethora of applications such as metalenses [1826], wave plates [2730], vortex beam convertors [3133], and holograms [3439] have been demonstrated. Nevertheless, devices based on a single-layered metasurface show symmetric transmission, i.e., identical transmission between the forward and backward directions. Recently, multilayered plasmonic metasurfaces have been proposed to demonstrate the asymmetric transmission. For example, Frese et al. [40] and Chen et al. [41] designed dual/triple-layered plasmonic metasurfaces to realize an asymmetric meta-hologram and directional Janus, respectively. Sun et al. [42] reported asymmetric dual-frequency meta-holograms based on a triple-layered transmissive metasurface. These metasurface-based asymmetric devices work in a narrow band because of their resonant unit cells. In addition, study on metasurface-based asymmetric focusing is limited. Here we propose and experimentally demonstrate the broadband asymmetric focusing in the terahertz (THz) region triggered by simultaneous manipulation of the propagation direction, phase, and polarization of THz waves. This directional device for asymmetric focusing consists of a geometric metasurface (with a nonresonant unit cell) and metallic gratings. Unlike previously demonstrated resonant metasurfaces with different structured unit cells, the geometric metasurface consists of identical microrods, and the manipulation of phase, polarization, and energy flux is realized by rotating the microrods (to generate geometric phase). Our designed dual-layered metasurfaces for asymmetric focusing can work in a wide bandwidth ranging from 0.6 to 1.1 THz. The approach for designing directional meta-devices may provide a robust platform in developing future high-performance devices that can manipulate EM waves in multiple degrees of freedom.

    2. DESIGN AND METHOD

    Figure 1 shows the schematic of an asymmetric focusing device. Upon the illumination of x-polarized THz waves in the forward direction, a focal point with orthogonal polarization is observed after the dual-layered metasurfaces. In contrast, for backward incidence, the focal point is not generated, leading to asymmetric focusing. This directional device consists of a geometric metasurface with spatially rotated microrods and metallic gratings, which act as a linear polarizer that can transmit x-polarized THz waves and reflect y-polarized THz waves. The geometric metasurface is designed to simultaneously control the phase (for focusing) and polarization (for polarization rotation), generating a focal point with orthogonal polarization relative to incident THz waves. Since a linearly polarized THz beam consists of left-handed circularly polarized (LCP) and right-handed circularly polarized (RCP) components with the same amplitude, two opposite phase profiles [see Eq. (1)] should be required to focus both the LCP and RCP components: where f is the focal length and λ is wavelength. Therefore, the focused LCP and RCP components [induced by phase modulations in Eq. (1)] with the same amplitude can be combined into a linearly polarized focal spot. Another bulk phase α that can be generated by rotating each microrod with an angle of α/2 should be introduced into the geometric metasurface to rotate the polarization orientation of the focal spot. Therefore, the desired phase profile of the geometric metasurface is governed by where α=90° is the polarization rotation angle of the focal spot.

    Schematic of asymmetric focusing. Under the illumination of x-polarized THz waves in the forward direction, a y-polarized focal spot is observed, while the focal spot is not generated for backward x-polarized incidence.

    Figure 1.Schematic of asymmetric focusing. Under the illumination of x-polarized THz waves in the forward direction, a y-polarized focal spot is observed, while the focal spot is not generated for backward x-polarized incidence.

    To manipulate the incident THz waves based on geometric phase, each microrod should be designed as a quasi-perfect half-wave plate. The optimized structure parameters of a gold microrod are as follows: the width and length of the rods are w1=50  μm and l=160  μm [see Fig. 2(a1)], and the period is 170 μm. Under the illumination of circularly polarized THz waves, the transmitted THz waves will contain two parts: one is the co-polarized THz waves, and the other is the cross-polarized THz waves. When the microrod is rotated counterclockwise/clockwise with an angle θ, the cross-polarized THz waves will be added with an additional phase delay of ±2θ as shown in Fig. 2(a2). Figure 2(b1) shows the schematic of metallic gratings where the long side of each grating is along the y axis. The width and period of these metallic gratings are 20 μm and 40 μm, respectively. As shown in Fig. 2(b2), the transmission for the incident x-polarized THz waves is much higher than that of the incident y-polarized THz waves. By introducing a geometric metasurface (with gold microrods) and metallic (gold) gratings [on the other side of the polyimide film (with permittivity of ε=3.5+0.035i)], asymmetric transmission is realized [see Figs. 2(c1) and 2(d2)]. For forward x-polarized incidence, the transmission is more than 82%, ranging from 0.6 to 1.1 THz, while it is lower than 52% under the the illumination of x-polarized THz waves from the backward direction, demonstrating an asymmetric transmission that can work in a broad bandwidth.

    Design of dual-layered metasurfaces: (a1) and (a2) schematic and the corresponding geometric phase of the microrod; (b1) and (b2) schematic and the corresponding transmission spectra of the metallic gratings; (c1) and (c2) schematic and the corresponding transmission spectra of the metasurface combined with metallic gratings; (d1) and (d2) optical images of the metasurface and metallic gratings.

    Figure 2.Design of dual-layered metasurfaces: (a1) and (a2) schematic and the corresponding geometric phase of the microrod; (b1) and (b2) schematic and the corresponding transmission spectra of the metallic gratings; (c1) and (c2) schematic and the corresponding transmission spectra of the metasurface combined with metallic gratings; (d1) and (d2) optical images of the metasurface and metallic gratings.

    3. RESULTS

    Figure 3 shows numerical simulations of asymmetric focusing. The size of this final metasurface composed of 80×80 microrods is 1.36  cm×1.36  cm. Based on the finite-difference time-domain (FDTD) method, the electric field distributions of asymmetric focusing are calculated and shown in Figs. 3(a1)3(f2). For simplicity, we calculated (and measured) electric field distributions at 0.6, 0.85, and 1.1 THz to show the broadband characteristic of our directional devices. Figure 3(a1) shows the electric field distribution of x-polarized transmission for forward THz waves at f=0.6  THz. The polarization of one focal spot observed at z=2.0  mm is along the y axis, which is orthogonal to the polarization of the incident THz waves, while the focal spot [see Fig. 3(b1)] is not generated for backward incidence, leading to asymmetric focusing. Here the mechanism of the asymmetric focusing can be understood as follows: for forward incidence, the x-polarized THz waves are almost completely transmitted through the metallic gratings and simultaneously interact with the metasurface (microrods with different orientations). Therefore, the transmitted x-polarized THz waves are partly converted into y-polarized THz waves and simultaneously focused into a focal spot. In contrast, for backward incidence, a portion of the waves interact with the metasurface and are converted into y-polarized THz waves. Then the converted waves are reflected completely after interacting with metallic gratings. The other part of nonconverted (x-polarized) THz waves can pass through metallic gratings. However, the nonconverted THz waves are not modulated by the metasurface and thus cannot be focused into a focal spot for backward waves. The detailed mechanism of asymmetric focusing under the illumination of x-polarized (and y-polarized) THz waves in the forward and backward directions is given in Appendix A. Figures 3(a2) and 3(b2) show electric field distributions in the xy plane (z=±2.0  mm). For forward transmission, a circularly shaped focal spot is observed at the real focal plane (z=2.0  mm), while it cannot be generated (at z=2.0  mm) for the backward case [see Figs. 3(a2) and 3(b2)]. When f=0.85 THz (1.1 THz), a focal point located at z=3.5  mm (z=5.5  mm) with polarization along the y axis is observed for forward incidence as shown in Figs. 3(c1) and 3(e1). In comparison with Figs. 3(a1), 3(c1), and 3(e1), the distance between the focal spot and the device is getting longer and longer due to the chromatic aberration. The focal spot cannot be generated for backward incidence [see Figs. 3(d1) and 3(f1)]. The corresponding electric field distributions in the xy plane are shown in Figs. 3(c2)3(f2), demonstrating asymmetric focusing between the forward and backward directions. The calculated efficiencies of the directional device at these three frequencies are 23.03%, 26.7%, and 21.34%, respectively. The calculated electric field (|Ex|2) distributions and the performances (broadband property, size, and excitation ratio of the focal spot) of this directional device are given in Appendix B and Appendix C, respectively. In addition, the focusing properties under the illumination of THz waves with different incident angles are shown in Appendix D.

    (a1)–(f1) Numerical simulation of electric field distributions in the x−z plane under the illumination of x-polarized THz waves in the forward/backward direction at 0.6, 0.85, and 1.1 THz; (a2)–(f2) the corresponding electric field distributions in the x−y plane.

    Figure 3.(a1)–(f1) Numerical simulation of electric field distributions in the xz plane under the illumination of x-polarized THz waves in the forward/backward direction at 0.6, 0.85, and 1.1 THz; (a2)–(f2) the corresponding electric field distributions in the xy plane.

    In experiment, we utilize near-field terahertz microscopy (NSTM) [26,43] to detect field distributions. The position of the THz source and THz tip are fixed while the sample is rotated along the y axis to measure the field distributions under the illumination of THz waves in the forward and backward directions. For forward x-polarized incidence (with working frequency of 0.6 THz), a focal point is generated with polarization along the y axis (at z=2.0  mm) as shown in Fig. 4(a1). In contrast, the focal spot is not generated for backward incidence [see Figs. 4(b1)], resulting in asymmetric focusing. Figures 4(c1)4(f1) show the measured electric field distributions at 0.85 THz and 1.1 THz, respectively. Under the illumination of x-polarized THz waves in the forward direction, a focal spot located at 3.5 mm (or 5.5 mm) is observed as shown in Figs. 4(c1) and 4(e1), while this focal spot cannot be generated for backward incidence [see Figs. 4(d1) and 4(f1)]. The measured electric field distributions in the xy plane at z=±2.0  mm, ±3.5  mm, ±5.5  mm are shown in Figs. 4(a2)4(f2), and they also reveal the asymmetric functionality. The measured efficiencies of the directional device at f=0.6, 0.85, and 1.1 THz are 20.6%, 25.3%, and 20.3%, respectively. The experimental results are matched well with the numerical simulations, except for a slight discrepancy, which can be attributed to the fabrication error and nonperfect half-wave plate of each microrod. In addition, the signal-to-noise ratio of the THz source (in the experiment) decreases with the increase of the frequency, and thus the focal spot is distorted in comparison with the case in simulation. The measured electric field (|Ex|2) distributions and the performances (size and excitation ratio of the focal spot) of the directional device are given in Appendix B and Appendix C, respectively.

    (a1)–(f1) The measured electric field distributions in the x−z plane under the illumination of x-polarized THz waves in the forward/backward direction at 0.6, 0.85, 1.1 THz; (a2)–(f2) the corresponding electric field distributions in the x−y plane.

    Figure 4.(a1)–(f1) The measured electric field distributions in the xz plane under the illumination of x-polarized THz waves in the forward/backward direction at 0.6, 0.85, 1.1 THz; (a2)–(f2) the corresponding electric field distributions in the xy plane.

    Our approach can not only manipulate asymmetric transmission with one focal point but also enable the capability to harness the multiple focal points with functionality of asymmetric transmission. We numerically demonstrate the asymmetric transmission with longitudinal and transverse multiple focal spots as shown in Fig. 5. The corresponding phase profile requirements of the longitudinal and transverse multiple focal spots are written as where φLCP1/φRCP1 and φLCP2/φRCP2 have the same formula as shown in Eq. (1) but different focal length [f1 (focal length) for φLCP1 and φRCP1, while f2 (focal length) for φLCP2 and φRCP2]. φLCP3(x,y)=φRCP3(x,y)=2πλ[(xxi)2+y2+f2f] and φLCP4(x,y)=φRCP4(x,y)=2πλ[(xxj)2+y2+f2f]. The designed parameters are xi=xj=1  mm, f1=4  mm, f2=6  mm, and f=4  mm. Figure 5(a) shows the electric field (|Ey|2) distribution of the asymmetric transmission of longitudinal multiple focal spots upon the illumination of x-polarized THz waves (at 0.85 THz) from the forward direction. Two separated focal points located at z=3.9  mm and z=5.9  mm are observed along the longitudinal direction. Here the discrepancy in focus spot position between theoretical design and numerical simulation results from the nonperfect half-wave plate of the microrods. When the x-polarized THz waves interact with the geometric metasurface (with the phase profile shown in ΦT), two transversely distributed focal points located at (x=1  mm, z=3.9  mm) and (x=1  mm, z=3.9  mm) are generated as shown in Fig. 5(c). These two focal points are not generated for backward incidence as shown in Figs. 5(b) and 5(d). The calculated efficiencies of the asymmetric transmission with longitudinal and transverse focal spots are 23.93% and 25.35%. The calculated electric field (|Ex|2) distributions and the performances (excitation ratio of the focal spot) of the directional device are given in Appendix B and Appendix C, respectively.

    Numerical simulation of asymmetric transmission with (a) and (b) longitudinal and (c) and (d) transversal multiple focal spots.

    Figure 5.Numerical simulation of asymmetric transmission with (a) and (b) longitudinal and (c) and (d) transversal multiple focal spots.

    4. DISCUSSION AND CONCLUSION

    The approach can be utilized to design ultrathin and ultraflat geometric metasurfaces by simultaneously controlling the polarization, phase, and transmission direction in all spatial dimensions, leading to unprecedented capabilities in developing multifunctional components and ultrathin systems. Typically, the asymmetric functionality of our designed metasurfaces can be traditionally realized by putting together a lens, a half-wave plate, and a linear polarizer. However, the conventional devices should typically be thick enough and exhibit certain surface topography to ensure the gradual phase accumulation (along the propagation direction) to realize the desired wave-manipulating functionalities. The thickness of these traditional THz elements is up to the order of centimeters, resulting in a large and bulky functional device. Since device miniaturization and system integration are two continuing trends for next-generation integration-optics applications, our designed directional device opens a new avenue to address the aforementioned issues. The thickness of the metasurface and metallic gratings is 0.15 μm, while it is 30 μm for the spacer (polyimide). The total thickness of the device is 30.3 μm, which is an order of 0.09λ0 (λ0 is the central working wavelength of the designed device). Therefore, the uniqueness and robustness of our proposed approach may provide a platform for designing future multifunctional and ultracompact devices. The sample fabrication is given in Appendix E. In addition, the efficiency plays an important role for device performance. The efficiency of the directional device is more than 20% [in simulation (see Fig. 10)] ranging from 0.6 to 1.1 THz. For the traditional devices, the efficiency of the THz lens (half-wave plate and linear polarizer) is about 52% (75% and 90%). Therefore, the total efficiency of the asymmetric focusing (realized by traditional devices) is about 35%. In comparison with the traditional bulk devices, our designed directional device not only shows a comparable performance in the efficiency but also enables the ability to design ultrathin devices.

    From the practical point of view, some other distinct properties of the designed directional device are of great interest. Unlike traditional asymmetric transmission, which only controls the propagation direction of EM waves, the developed directional device can simultaneously control the polarization, phase, and transmission direction, leading to a directional metalens. As is well known, lenses can be applied for imaging. Recently, various types of metalenses such as dual-polarity plasmonic metalenses [18,19], achromatic metalenses [23,44], metalens arrays [45], and polarization-dependent metalenses [26,46] have been developed for imaging. However, all of the previous metalenses are symmetric in the forward and backward directions. The developed directional metalens can be applied for unidirectional imaging, which means that we can grab information/imaging from the “enemy” side and avoid our own information being detected. As a result, the designed ultrathin device can have an immediate impact on information security and cloaking. Moreover, the proposed approach can be extended to design asymmetric transmission with multiple focal points, while a traditional lens generates only one focal point. Therefore, the flexible design will enable the directional metalens with unprecedented functions that are impossible to achieve by traditional lenses. Due to the simplicity and robustness of the method, our approach can be easily extended into the visible and IR frequency regions by designing the metasurface and metallic gratings with proper sizes, i.e., the same structure with reduced size. While remarkable progress has been achieved, great challenges, i.e., further enhancement of the focusing efficiency and asymmetric focusing under oblique incidence, exist that provide new opportunities for future studies.

    In conclusion, we have proposed and experimentally demonstrated an approach to realizing asymmetric focusing in the THz region. Based on a geometric metasurface and polarization-selective gratings, a focal spot is observed under the illumination of x-polarized THz waves in the forward direction, while this focal spot is not generated for backward incidence. This metasurface-based directional device works in a broad bandwidth ranging from 0.6 to 1.1 THz. The robust approach can be also extended to design asymmetric transmission with multiple focal points. The presented dual metasurfaces are a powerful tool that can not only flexibly control the wavefront of EM but also open a new window to manipulate the propagation direction of EM waves.

    APPENDIX A: PRINCIPLE OF ASYMMETRIC FOCUSING

    Metasurface with Linearly Polarized Focusing and Polarization Rotation

    We assume a (transmission-type) metasurface with the long axis of each metallic microrod along the x axis. Under the normal illumination of a linearly polarized (x-polarized) THz beam, the polarization of the transmitted beam is the same as the incident THz waves. The Jones vector of the transmitted beam can be written as [10]=12(12[1i]+12[1?i]).If each of the microrods is rotated counterclockwise with an angle of α2, the transmitted beam will contain two parts. One part is a polarization-rotated (converted) beam in which the polarization is rotated counterclockwise with an angle (α). The other part is a co-polarized (nonconverted) beam without polarization rotation. It can be understood as follows: under the illumination of a linearly polarized THz beam, the LCP component is partially converted into an RCP THz beam with an abrupt phase delay of α, while the RCP component is partially converted into LCP THz waves with an additional phase delay of ?α. The transmitted electric field for the converted LCP and RCP beams can be written as ELCP/RCP=η(λ)e?iα[1±i],where η(λ) is the conversion efficiency. Both of these converted (LCP/RCP) components can be combined to a linearly polarized THz beam because of the same transmission intensity. The electric field is expressed as follows: Econ=η(λ)e?iα[1i]+η(λ)eiα[1?i]=η(λ)[cos?αsin?α].As a consequence, the polarization of the converted (linearly polarized) THz beam is rotated counterclockwise with an angle of α, which is realized by counterclockwise rotation of each microrod with an angle of α2.

    To enable the functionality of focusing, two opposite phase profiles are introduced to focus both the LCP and RCP components of the converted THz waves. One phase profile is φ=φR(x,y)=2πλ(x2+y2+f2?|f|) (f is the focal length), and the other phase profile is φL(x,y)=?φR(x,y)=?2πλ(x2+y2+f2?|f|). Here the first phase profile can focus (defocus) the RCP (LCP) component of the converted THz waves, while the second phase profile can focus (defocus) the LCP (RCP) component of the converted THz waves. The Jones vector of the converted (polarization-rotated) beam under the modulation of the aforementioned two opposite phase profiles can be given as [21] η(λ)2{12[1i]exp(?iα)exp(iφL(x,y))+12[1?i]exp(iα)exp(?iφL(x,y))}+η(λ)2{12[1i]exp(?iα)exp(iφR(x,y))+12[1?i]exp(iα)exp(?iφR(x,y))}.The above equation can also be expressed as follows: η(λ)exp(iφL)[cos?αsin?α]+η(λ)exp(iφR)[cos?αsin?α].Accordingly, the required phase distribution for our proposed metalens is governed by Φ(x,y)=arg{exp[i(α+φL(x,y))]+exp[i(α+φR(x,y))]}.It should be noted that α is the rotation angle of polarization, while φL(x,y) and φR(x,y) are the required phase profiles for focusing the LCP and RCP components. Therefore, under the illumination of a linearly polarized THz beam, it can be focused into a linearly polarized focal point with the corresponding polarization rotated with an angle of α. Here we neglect the coefficient of conversion efficiency because it does not affect the total phase profile as shown in Eq.?(6).

    Metasurface Combined with Metallic Gratings for Asymmetric Focusing

    Figure?6 shows the schematics [see Figs.?6(a1)6(c1)], transmission, and reflection of a microrod, metallic gratings, and a unit cell of the directional device under the illumination of linearly polarized THz waves. Figures?6(a2)6(a4) illustrate the co-polarized, cross-polarized, and total transmission/reflection of the microrod under the illumination of linearly polarized THz waves. It should be noted that the microrod enables functionality of a quasi-perfect half-wave plate, and the intersection angle between the long axis of the microrod and the x axis is 45°. Therefore, under the illumination of x-polarized (or y-polarized) THz waves, the transmitted/reflected waves contain two parts: one has the same polarization as the incident THz waves, and the other has the cross-polarization. As shown in Fig.?6(a1), the single-layered metasurface (microrods) shows perfect symmetry to the incident THz waves in the forward and backward directions. As a consequence, the transmission (or reflection) of the microrod upon the illumination of linearly polarized THz waves in the forward direction is the same as that for backward incidence as shown in Figs.?6(a2)6(a4). For metallic gratings [see Fig.?6(b1)], the co-polarized transmission of x-polarized THz waves for forward and backward incidence is more than 86%, ranging from 0.6 to 1.1?THz, while the co-polarized reflection in this case is less than 8% as shown in Fig.?6(b2). For y-polarized incidence, the THz wave is reflected completely [Ryy?(forward)=Ryy?(backward)=100% in Fig.?6(b2)] without transmission. In addition, the cross-polarized transmission and reflection are all near zero [see Fig.?6(b3)]. Therefore, the designed metallic gratings with a long axis along the y axis can block y-polarized THz waves and transmit x-polarized THz waves as shown in Fig.?6(b4).

    (a1)–(c1) Schematics of a microrod, metallic gratings, and a unit cell of the directional device. The intersection angle between the long axis of the microrod and the x axis is 45°, while the long axis of gratings is along the x axis. (a2)–(a4) The co-polarized/cross-polarized/total transmission and reflection of the microrod under the illumination of linearly polarized THz waves. (b2)–(b4) The co-polarized/cross-polarized/total transmission and reflection of the metallic gratings under the illumination of linearly polarized THz waves. (c2)–(c4) The co-polarized/cross-polarized/total transmission and reflection of the unit cell of the directional device under the illumination of linearly polarized THz waves. Tij (Rij) is the transmission (reflection) of the i-polarized THz waves under the illumination of j-polarized THz waves (i,j=x,y). Ti (Ri,i=x,y) is the total transmission (total reflection) under the illumination of i-polarized THz waves.

    Figure 6.(a1)–(c1) Schematics of a microrod, metallic gratings, and a unit cell of the directional device. The intersection angle between the long axis of the microrod and the x axis is 45°, while the long axis of gratings is along the x axis. (a2)–(a4) The co-polarized/cross-polarized/total transmission and reflection of the microrod under the illumination of linearly polarized THz waves. (b2)–(b4) The co-polarized/cross-polarized/total transmission and reflection of the metallic gratings under the illumination of linearly polarized THz waves. (c2)–(c4) The co-polarized/cross-polarized/total transmission and reflection of the unit cell of the directional device under the illumination of linearly polarized THz waves. Tij(Rij) is the transmission (reflection) of the i-polarized THz waves under the illumination of j-polarized THz waves (i,j=x,y). Ti(Ri,i=x,y) is the total transmission (total reflection) under the illumination of i-polarized THz waves.

    (a1)–(f1) Numerical simulation of electric field distributions in the x−z plane under the illumination of y-polarized THz waves in the forward direction at 0.6, 0.85, and 1.1 THz; (a2)–(f2) the calculated electric field distributions for backward incidence.

    Figure 7.(a1)–(f1) Numerical simulation of electric field distributions in the xz plane under the illumination of y-polarized THz waves in the forward direction at 0.6, 0.85, and 1.1 THz; (a2)–(f2) the calculated electric field distributions for backward incidence.

    Figure?7 shows the calculated electric field distributions of the directional device under the illumination of y-polarized THz waves in the forward and backward directions. As shown in Figs.?7(a1) and 7(b1), the focal spot cannot be generated for forward incidence. In contrast, an x-polarized focal point is observed for backward incidence [see Figs.?7(a2) and 7(b2)], indicating that a portion of the y-polarized THz waves (for backward incidence) are converted into x-polarized THz waves and simultaneously focused into a focal spot. For f=0.85??THz (1.1?THz), an x-polarized focal point can also be observed for backward incidence, while it is not generated for forward incidence [see Figs.?7(c1)7(f1) and 7(c2)7(f2)]. In comparison with Figs.?7(b2), 7(d2), and 7(f2), the distance between the focal spot and the device is getting longer and longer, which can be attributed to the chromatic aberration.

    APPENDIX B: ELECTRIC FIELD (Ex) DISTRIBUTIONS OF THE DIRECTIONAL DEVICE

    Figure?8 shows the calculated and measured electric field (|Ex|2) distributions of the directional device (that can unidirectionally generate a y-polarized focal point) under the illumination of x-polarized THz waves in the forward and backward directions. As shown in Figs.?8(a1), 8(c1), and 8(e1), no focal spot can be observed under the illumination of x-polarized THz waves in the forward direction. In experiment, the corresponding electric field (|Ex|2) distributions in Figs.?8(b1), 8(d1), and 8(f1) also demonstrate that there is no (x-polarized) focal spot. For backward incidence, the designed directional device cannot modulate nonconverted THz waves into a focal spot as shown in Figs.?8(a2)8(f2).

    The calculated and measured electric field distributions in the x−z plane under the illumination of the x-polarized THz waves from the (a1)–(f1) forward and (a2)–(f2) backward directions at 0.6, 0.85, and 1.1 THz.

    Figure 8.The calculated and measured electric field distributions in the xz plane under the illumination of the x-polarized THz waves from the (a1)–(f1) forward and (a2)–(f2) backward directions at 0.6, 0.85, and 1.1 THz.

    The calculated electric field distributions in the x−z plane under the illumination of the x-polarized THz waves in the (a1), (b1) forward and (a2), (b2) backward directions at 0.85 THz.

    Figure 9.The calculated electric field distributions in the xz plane under the illumination of the x-polarized THz waves in the (a1), (b1) forward and (a2), (b2) backward directions at 0.85 THz.

    Calculated efficiency of the directional device under the illumination of x-polarized THz waves in the forward direction.

    Figure 10.Calculated efficiency of the directional device under the illumination of x-polarized THz waves in the forward direction.

    APPENDIX C: PERFORMANCE OF ASYMMETRIC FOCUSING

    Broadband Property of Asymmetric Focusing

    Figure?10 shows the calculated efficiency of our designed device under the illumination of x-polarized THz waves in the forward direction. The efficiency is more than 20% ranging from 0.6 to 1.1?THz, demonstrating a broadband directional device.

    Figure?11 shows the schematic of multiple transmissions from the dual-layered metasurfaces. For forward x-polarized incidence, the total output THz waves are the superposition of all the transmitted THz waves passing through the microrods. The transmitted THz waves are expressed as follows: Eout=Eout1+Eout2+Eout3+?

    Schematic of multiple transmissions from the dual-layered metasurfaces.

    Figure 11.Schematic of multiple transmissions from the dual-layered metasurfaces.

    Comparison of the numerical (blue curves) and experimental (red curves) focusing properties: (a)–(c) the corresponding electric field distributions at x=0 in the focal plane.

    Figure 12.Comparison of the numerical (blue curves) and experimental (red curves) focusing properties: (a)–(c) the corresponding electric field distributions at x=0 in the focal plane.

    Schematics for the extinction ratio defined as (a) TEy/TEx and (b) TEy1/TEy2.

    Figure 13.Schematics for the extinction ratio defined as (a) TEy/TEx and (b) TEy1/TEy2.

    Extinction Ratio Between |Ey|2 and |Ex|2

    Frequency (THz)0.60.8511.1
    Sim. extinction ratio1.22:13.79:12.27:1
    Exp. extinction ratio1.6:11.7:11.6:1

    For traditional asymmetric transmission, the extinction ratio is defined as the ratio of the transmission coefficient for forward incidence to that for backward incidence. Therefore, the extinction ratio of the asymmetric focusing in our paper can also be analogously defined as the ratio of the Ey component for forward x-polarized incidence to that for backward incidence [40] [extinction ratio=TEy1/TEy2 and |Einc1|2(forward)=|Einc2|2(backward)] as shown in Fig.?13(b). Here TEy1 and TEy2 are defined as the transmission of the Ey component under the illumination of x-polarized THz waves in the forward and backward directions. In this case, the extinction ratios at 0.6, 0.85, and 1.1?THz are 36.6:1, 81.5:1, and 76.2:1 as shown in Table?4. In the experiment, the corresponding extinction ratios are 17:1, 71:1, and 15.5:1 at 0.6, 0.85, and 1.1?THz.

    Extinction Ratio Between the Forward and Backward Directions

    Frequency (THz)0.60.851.1
    Sim. extinction ratio36.6:181.5:176:1
    Exp. extinction ratio17:171:115.5:1

    The extinction ratios (defined as the ratio of the transmission coefficient between Ex and Ey components under the illumination of x-polarized THz waves in the forward direction) of the directional device with two unidirectional focal points are shown in Table?5. They are 5.4:1 and 6:1, respectively. For the extinction ratio, defined as the ratio of the transmission coefficient (|Ey|2) for forward waves to that for backward waves, they are 82:1 and 38:1 as shown in Table?5.

    Extinction Ratio for the Directional Device with Two Focal Spots

    Extinction RatioLongitudinalTransversal
    |Ey|2/|Ex|25.4:16:1
    |Ey1|2/|Ey2|282:138:1

    APPENDIX D: PERFORMANCES OF THE DIRECTIONAL DEVICE UNDER THE ILLUMINATION OF THz WAVES WITH DIFFERENT INCIDENT ANGLES

    Figure?14 shows the calculated electric field (|Ey|2) distributions of the directional device under the illumination of x-polarized THz waves with different incident angles in the forward and backward directions. (The corresponding electric field distributions at f=0.6??THz and 1.1?THz show the same trend as in Fig.?14.) For incident angles ranging from 10° to 30°, the focal spot is more and more?deviated from the center of the device and shifted to the blow with the increase of incident angle [see Figs.?14(a1)14(c1)]. When the incident THz waves with a tilted wavefront (nonzero incident angle) interact with the directional device, the original position for generating the constructive interference (at x=0 and z=3.5??mm) is shifted downward, and thus the focal spot is?deviated to the blow. Since the original geometric phase with respect to normally incident THz waves shows circular symmetry, it cannot keep the circular symmetry under the illumination of x-polarized THz waves with nonzero incident angles. Therefore, the deviated focal spots will inevitably be distorted as shown in Figs.?14(a1)14(c1). For backward incidence, no focal spot can be observed [see Figs.?14(a2)14(c2)].

    (a1)–(c1) Calculated electric field (|Ey|2) distributions in the x−z plane under the illumination of x-polarized THz waves (with different incident angles) in the forward directions at 0.85 THz; (a2)–(c2) the calculated electric field distributions for backward incidence. Insets show the schematics for the incident THz waves with a tilted wavefront.

    Figure 14.(a1)–(c1) Calculated electric field (|Ey|2) distributions in the xz plane under the illumination of x-polarized THz waves (with different incident angles) in the forward directions at 0.85 THz; (a2)–(c2) the calculated electric field distributions for backward incidence. Insets show the schematics for the incident THz waves with a tilted wavefront.

    APPENDIX E: SAMPLE FABRICATION

    We fabricated the directional device based on traditional photolithography. The AZP5214 photoresist (with thickness of 1.6?μm) was spun (with spin speed of 3000 r/min) onto one side of the polyimide film. Then a mask was utilized for exposure processing. After the metal coating and ultrasonic striping, gold microrods were formed, resulting in a metasurface. The gold gratings can also be fabricated on the other side of the polyimide film by repeating the same processing.

    References

    [1] S. Cakmakyapan, H. Caglayan, A. E. Serebryannikov, E. Ozbay. Experimental validation of strong directional selectivity in nonsymmetric metallic gratings with a subwavelength slit. Appl. Phys. Lett., 98, 051103(2011).

    [2] D. Dai, Z. Wang, J. E. Bowers. Ultrashort broadband polarization beam splitter based on an asymmetrical directional coupler. Opt. Lett., 36, 2590-2592(2011).

    [3] S. Cakmakyapan, A. E. Serebryannikov, H. Caglayan, E. Ozbay. Spoof-plasmon relevant one-way collimation and multiplexing at beaming from a slit in metallic grating. Opt. Express, 20, 26636-26648(2012).

    [4] D. F. Tang, C. Wang, W. K. Pan, M. H. Li, J. F. Dong. Broad dual-band asymmetric transmission of circular polarized waves in near-infrared communication band. Opt. Express, 25, 11329-11339(2017).

    [5] A. Cicek, M. B. Yucel, O. A. Kaya, B. Ulug. Refraction-based photonic crystal diode. Opt. Lett., 37, 2937-2939(2012).

    [6] D. L. Sounas, C. Caloz. Electromagnetic nonreciprocity and gyrotropy of graphene. Appl. Phys. Lett., 98, 021911(2011).

    [7] A. Shaltout, A. Kildishev, V. Shalaev. Time-varying metasurfaces and Lorentz non-reciprocity. Opt. Mater. Express, 5, 2459-2467(2015).

    [8] D. L. Sounas, A. Alu. Non-reciprocal photonics based on time modulation. Nat. Photonics, 11, 774-783(2017).

    [9] A. M. Mahmoud, A. R. Davoyan, N. Engheta. All-passive nonreciprocal metastructure. Nat. Commun., 6, 8359(2015).

    [10] V. A. Fedotov, P. L. Mladyonov, S. L. Prosvirnin, A. V. Rogacheva, Y. Chen, N. I. Zheludev. Asymmetric propagation of electromagnetic waves through a planar chiral structure. Phys. Rev. Lett., 97, 167401(2006).

    [11] C. Huang, Y. Feng, J. Zhao, Z. Wang, T. Jiang. Asymmetric electromagnetic wave transmission of linear polarization via polarization conversion through chiral metamaterial structures. Phys. Rev. B, 85, 195131(2012).

    [12] Z. Li, M. Gokkavas, E. Ozbay. Manipulation of asymmetric transmission in planar chiral nanostructures by anisotropic loss. Adv. Opt. Mater., 1, 482-488(2013).

    [13] C. Menzel, C. Helgert, C. Rockstuhl, E.-B. Kley, A. Tünnermann, T. Pertsch, F. Lederer. Asymmetric transmission of linearly polarized light at optical metamaterials. Phys. Rev. Lett., 104, 253902(2010).

    [14] V. A. Fedotov, A. S. Schwanecke, N. I. Zheludev, V. V. Khardikov, S. L. Prosvirnin. Asymmetric transmission of light and enantiomerically sensitive plasmon resonance in planar chiral nanostructures. Nano Lett., 7, 1996-1999(2007).

    [15] M. Stolarek, D. Yavorskiy, R. Kotynski, C. J. Z. Rodriguez, J. Lusakowski, T. Szoplik. Asymmetric transmission of terahertz radiation through a double grating. Opt. Lett., 38, 839-841(2013).

    [16] M. Mutlu, A. E. Akosman, A. E. Serebryannikov, E. Ozbay. Diodelike asymmetric transmission of linearly polarized waves using magnetoelectric coupling and electromagnetic wave tunneling. Phys. Rev. Lett., 108, 213905(2012).

    [17] H. Kurt, D. Yilmaz, A. E. Akosman, E. Ozbay. Asymmetric light propagation in chirped photonic crystal waveguides. Opt. Express, 20, 20635-20646(2012).

    [18] X. Chen, L. Huang, H. Mühlenbernd, G. Li, B. Bai, Q. Tan, G. Jin, C. Qiu, S. Zhang, T. Zentgraf. Dual-polarity plasmonic metalens for visible light. Nat. Commun., 3, 1198(2012).

    [19] X. Chen, M. Chen, M. Q. Mehmood, D. Wen, F. Yue, C. Qiu, S. Zhang. Longitudinal multifocimetalens for circularly polarized light. Adv. Opt. Mater., 3, 1201-1206(2015).

    [20] Z. Zhang, D. Wen, C. Zhang, M. Chen, W. Wang, S. Chen, X. Chen. Multifunctional light sword metasurface lens. ACS Photon., 5, 1794-1799(2018).

    [21] X. Zang, H. Ding, Y. Intaravanne, L. Chen, Y. Peng, Q. Ke, A. V. Balakin, A. P. Shkurinov, X. Chen, Y. Zhu, S. Zhuang. A multi-foci metalens with polarization-rotated focal points. Laser Photon. Rev., 13, 1900182(2019).

    [22] A. Arbabi, Y. Horie, A. J. Ball, A. Faraon. Subwavelength-thick lenses with high numerical apertures and large efficiency based on high-contrast transmitarrays. Nat. Commun., 6, 7069(2015).

    [23] S. Wang, P. Wu, V. Su, Y. Lai, M. Chen, H. Kuo, B. Chen, Y. Chen, T. Huang, J. Wang, R. Lin, C. Kuan, T. Li, Z. Wang, S. Zhu, D. Tsai. A broadband achromatic metalens in the visible. Nat. Nanotechnol., 13, 227-232(2018).

    [24] E. Arbabi, A. Arbabi, S. M. Kamali, Y. Horie, A. Faraon. Multiwavelength polarization-insensitive lenses based on dielectric metasurfaces with meta-molecules. Optica, 3, 628-633(2016).

    [25] N. Davidson, A. Friesem, E. Hasman. Holographic axilens: high resolution and long focal depth. Opt. Lett., 16, 523-525(1991).

    [26] X. Zang, W. Xu, M. Gu, B. Yao, L. Chen, Y. Peng, J. Xie, A. V. Balakin, A. P. Shkurinov, Y. Zhu, S. Zhuang. Polarization-insensitive metalens with extended focal depth and longitudinal high-tolerance imaging. Adv. Opt. Mater., 8, 1901342(2020).

    [27] N. Yu, F. Aieta, P. Genevet, M. A. Kats, Z. Gaburro, F. Capasso. A broadband, background-free quarter-wave plate based on plasmonic metasurfaces. Nano Lett., 12, 6328-6333(2012).

    [28] B. Yang, W. M. Ye, X. D. Yuan, Z. H. Zhu, C. Zeng. Design of ultrathin plasmonic quarter-wave plate based on period coupling. Opt. Lett., 38, 679-681(2013).

    [29] R. Fan, Y. Zhou, X. Ren, R. Peng, S. Jiang, D. Xu, X. Xiong, X. Huang, M. Wang. Freely tunable broadband polarization rotator for terahertz waves. Adv. Mater., 27, 1201-1206(2015).

    [30] S. Jiang, X. Xiong, Y. Hu, Y. Hu, G. Ma, R. Peng, C. Sun, M. Wang. Controlling the polarization state of light with a dispersion-free metastructure. Phys. Rev. X., 4, 021026(2014).

    [31] E. Karimi, S. A. Schulz, I. D. Leon, H. Qassim, J. Upham, R. W. Boyd. Generating optical orbital angular momentum at visible wavelengths using a plasmonic metasurface. Light Sci. Appl., 3, e167(2014).

    [32] F. Yue, D. Wen, C. Zhang, B. D. Gerardot, W. Wang, S. Zhang, X. Chen. Multichannel polarization-controllable superpositions of orbital angular momentum states. Adv. Mater., 29, 1603838(2017).

    [33] X. Zang, Y. Zhu, C. Mao, W. Xu, H. Ding, J. Xie, Q. Cheng, L. Chen, Y. Peng, Q. Hu, M. Gu, S. Zhuang. Manipulating terahertz plasmonic vortex based on geometric and dynamic phase. Adv. Opt. Mater., 7, 1801328(2018).

    [34] D. Wen, D. F. Yue, G. Li, G. Zheng, K. Chan, S. Chen, M. Chen, K. F. Li, P. W. H. Wong, K. W. Cheah, E. Y. B. Pun, S. Zhang, X. Chen. Helicity multiplexed broadband metasurface holograms. Nat. Commun., 6, 8241(2015).

    [35] G. Zheng, H. Mühlenbernd, M. Kenney, G. Li, T. Zentgraf, S. Zhang. Metasurface holograms reaching 80% efficiency. Nat. Nanotechnol., 10, 308-312(2015).

    [36] B. Wang, F. Dong, Q. T. Li, D. Yang, C. Sun, J. Chen, Z. Song, L. Xu, W. Chu, Y. F. Xiao, Q. Gong, Y. Li. Visible-frequency dielectric metasurfaces for multiwavelength achromatic and highly dispersive holograms. Nano Lett., 16, 5235-5240(2016).

    [37] X. Zang, F. Dong, F. Yue, C. Zhang, L. Xu, Z. Song, M. Chen, P. Chen, G. S. Buller, Y. Zhu, S. Zhuang, W. Chu, S. Zhang, X. Chen. Polarization encoded color image embedded in a dielectric metasurface. Adv. Mater., 30, 1707499(2018).

    [38] L. Jin, Z. Dong, S. Mei, Y. Yu, Z. Wei, Z. Pan, S. Rezaei, X. Li, A. I. Kuznetsov, Y. S. Kivshar, J. K. W. Yang, C. W. Qiu. Noninterleaved metasurface for (26−1) spin- and wavelength-encoded holograms. Nano Lett., 18, 8016-8024(2018).

    [39] Y. W. Huang, W. T. Chen, W. Tsai, P. Wu, C. Wang, G. Sun, D. P. Tsai. Aluminum plasmonic multicolor meta-hologram. Nano. Lett., 15, 3122-3127(2015).

    [40] D. Frese, Q. Wei, Y. Wang, L. Huang, T. Zentgraf. Nonreciprocal asymmetric polarization encryption by layered plasmonicmetasurfaces. Nano. Lett., 19, 3976-3980(2019).

    [41] K. Chen, G. Ding, G. Hu, Z. Jin, J. Zhao, Y. Feng, T. Jiang, A. Alu, C. W. Qiu. Directional Janus metasurface. Adv. Mater., 32, 1906352(2019).

    [42] Q. Sun, Z. Zhang, Y. Huang, X. Ma, M. Pu, Y. Guo, X. Li, X. Luo. Asymmetric transmission and wavefrontmanipulation toward dual-frequency meta-holograms. ACS Photon., 6, 1541-1546(2019).

    [43] X. Zang, C. Miao, X. Guo, G. You, H. Yang, L. Chen, Y. Zhu, S. Zhuang. Polarization-controlled terahertz superfocusing. Appl. Phys. Lett., 113, 071102(2018).

    [44] W. T. Chen, A. Y. Zhu, V. Sanjeev, M. Khorasaninejad, Z. J. Shi, E. Lee, F. Capasso. A broadband achromatic metalens for focusing and imaging in the visible. Nat. Nanotechnol., 13, 220-226(2019).

    [45] R. J. Lin, V.-C. Su, S. Wang, M. K. Chen, T. L. Chung, Y. H. Chen, H. Y. Kuo, J.-W. Chen, J. Chen, Y.-T. Huang, J.-H. Wang, C. H. Chu, P. C. Wu, T. Li, Z. Wang, S. Zhu, D. P. Tsai. Achromatic metalens array for full-colour light-field imaging. Nat. Nanotechnol., 14, 227-231(2018).

    [46] W. T. Chen, A. Y. Zhu, J. Sisler, Z. Bharwani, F. Capasso. A broadband achromatic polarization-insensitive metalens consisting of anisotropic nanostructures. Nat. Commun., 10, 355(2019).

    [47] F. Yue, C. Zhang, X. Zang, D. Wei, B. D. Gerardot, S. Zhang, X. Chen. High-resolution grayscale image hidden in a laser beam. Light Sci. Appl., 7, 17129(2018).

    Bingshuang Yao, Xiaofei Zang, Zhen Li, Lin Chen, Jingya Xie, Yiming Zhu, Songlin Zhuang. Dual-layered metasurfaces for asymmetric focusing[J]. Photonics Research, 2020, 8(6): 830
    Download Citation