• Photonics Research
  • Vol. 9, Issue 4, 567 (2021)
Jie Li1、†, Chenglong Zheng1、†, Guocui Wang2、3, Jitao Li1, Hongliang Zhao1, Yue Yang1, Zhang Zhang1, Maosheng Yang4, Liang Wu1, Jining Li1, Yating Zhang1、5、*, Yan Zhang2, and Jianquan Yao1、6、*
Author Affiliations
  • 1Key Laboratory of Opto-Electronics Information Technology (Tianjin University), Ministry of Education, School of Precision Instruments and Opto-Electronics Engineering, Tianjin University, Tianjin 300072, China
  • 2Beijing Key Laboratory for Metamaterials and Devices, Key Laboratory of Terahertz Optoelectronics, Ministry of Education, and Beijing Advanced Innovation Center for Imaging Technology, Department of Physics, Capital Normal University, Beijing 100048, China
  • 3Beijing Engineering Research Center for Mixed Reality and Advanced Display, School of Optics and Photonics, Beijing Institute of Technology, Beijing 100081, China
  • 4School of Mechanical Engineering, Jiangsu University, Zhenjiang 225009, China
  • 5e-mail: yating@tju.edu.cn
  • 6e-mail: jqyao@tju.edu.cn
  • show less
    DOI: 10.1364/PRJ.415547 Cite this Article Set citation alerts
    Jie Li, Chenglong Zheng, Guocui Wang, Jitao Li, Hongliang Zhao, Yue Yang, Zhang Zhang, Maosheng Yang, Liang Wu, Jining Li, Yating Zhang, Yan Zhang, Jianquan Yao, "Circular dichroism-like response of terahertz wave caused by phase manipulation via all-silicon metasurface," Photonics Res. 9, 567 (2021) Copy Citation Text show less

    Abstract

    Chiral metasurfaces based on asymmetric meta-atoms have achieved artificial circular dichroism (CD), spin-dependent wavefront control, near-field imaging, and other spin-related electromagnetic control. In this paper, we propose and experimentally verify a scheme for achieving high-efficiency chiral response similar to CD of terahertz (THz) wave via phase manipulation. By introducing the geometric phase and dynamic phase in an all-silicon metasurface, the spin-decoupled terahertz transmission is obtained. The giant circular dichroism-like effect in the transmission spectrum is observed by using a random phase distribution for one of the circular polarization components. More importantly, the effect can be adjusted when we change the area of the metasurface illuminated by an incident terahertz beam. In addition, we also demonstrate the spin-dependent arbitrary wavefront control of the transmitted terahertz wave, in which one of the circularly polarized components is scattered, while the other forms a focused vortex beam. Simulated and experimental results show that this method provides a new idea for spin selective control of THz waves.

    1. INTRODUCTION

    Chirality is a special geometric asymmetry in nature, which exists widely in microscopic matter and macroscopic objects [1]. The chiral structure of the optical medium may trigger a chiral optical response, which can be explained as different refractive indices for left- and right-handed circularly polarized (LCP and RCP) light [2]. The difference between real and imaginary parts of the refractive index is manifested as optical activity and circular dichroism (CD), that is, the difference in the propagation speed or absorption of the two components. CD has become one of the important expressions of the molecular conformation of biological macromolecules, DNA, and many drug molecules, which is widely used in biology, medicine, chemistry, and other fields [3,4]. However, the inherent CD in natural materials is extremely weak, which requires a high-precision testing system [5].

    Metamaterials are new artificial structures composed of subwavelength meta-atoms, which have achieved a variety of electromagnetic manipulation and novel physical phenomena. Metasurfaces are two-dimensional metamaterials that have realized efficient electromagnetic control on extremely small longitudinal scales [6]. Both metamaterials and metasurfaces have achieved broadband or giant chiral optical phenomena, including CD [7,8], chiral wavefront control [9], circular polarization imaging or detection [10,11], and tunable chiral response [12]. In addition, both metasurfaces with or without chirality have been reported for greatly enhancing the CD of chiral molecules [1315]. In the past 20 years, various chiral phenomena have been observed in metasurfaces includes intrinsic chirality caused by asymmetric meta-atoms and extrinsic chirality caused by achiral structures under oblique incidence [16,17]. The realized functions include spin-selective absorption or reflection [18,19], polarization rotation, and so on [2027]. However, these implementations require careful design of asymmetric units or large incident angles, and most of them can only work in a narrow frequency band or cannot achieve arbitrary wavefront control.

    In this paper, we propose a new method to achieve chiral optical response based on phase manipulation rather than asymmetry meta-atoms. Simultaneously introducing geometric phase and dynamic phase into an all-silicon metasurface, spin-decoupled phase control of the transmitted terahertz wave is realized. The randomly arranged phase of one of the circularly polarized components causes giant circular dichroism-like effect in transmission. The CD can be controlled by the size of the incident terahertz beam. Not only that, but we also demonstrate another metasurface that achieves arbitrary wavefront control of the circularly polarized wave which is not scattered. This scheme provides a new idea for spin-dependent terahertz wave control.

    2. RESULTS AND DISCUSSION

    Schematic diagram of the metasurface which can independently control the phase of two circularly polarized incident waves. The cells are cross-shaped elliptical columns; both the cells and the substrate are made of high-resistance silicon.

    Figure 1.Schematic diagram of the metasurface which can independently control the phase of two circularly polarized incident waves. The cells are cross-shaped elliptical columns; both the cells and the substrate are made of high-resistance silicon.

    In order to achieve the above functions, design of the metasurface units with spin decoupling function is necessary, so the dynamic phase and geometric phase (Pancharatnam–Berry phase, P-B phase) need to be introduced at the same time. The P-B phase means an additional phase in circularly polarized electromagnetic waves determined by the geometric orientation of the meta-atoms, which appears in anisotropic structures. The P-B phase is spin-dependent, which means the two circular polarization components always obtain opposite phases, while the dynamic phase is determined by the geometry of the meta-atom. Assuming that the dynamic phases of the anisotropic element for x and y polarization components are φx and φy, the geometric phase can be obtained by rotating the element around the z axis by an angle of θ. To achieve the independently controllable circular polarization phase, the dynamic phase and rotation angle of the units need to meet the following equations [28]: {φx=(φRCP+φLCP)/2φy=(φRCP+φLCP)/2πθ=(φRCPφLCP)/4,where φRCP and φLCP represent the circular polarization phase response of the metasurface units.

    Simulation and results selection of the geometric parameters of meta-atoms. (a) Schematic diagram of the units. (b)–(e) The transmission amplitudes and phase shifts of the silicon cylinder for x- and y-polarized terahertz waves with different values of Lx and Ly. (f), (g) The selected 64 sets of Ly and Lx values, in which the phases are divided into eight values at intervals of 45 deg.

    Figure 2.Simulation and results selection of the geometric parameters of meta-atoms. (a) Schematic diagram of the units. (b)–(e) The transmission amplitudes and phase shifts of the silicon cylinder for x- and y-polarized terahertz waves with different values of Lx and Ly. (f), (g) The selected 64 sets of Ly and Lx values, in which the phases are divided into eight values at intervals of 45 deg.

    Phase design of the metasurface with circular dichroism-like effect consists of 140×140 units. (a) The phase matrices consist of equal elements or pseudorandom numbers (submatrix with equal elements). (b) The phase distributions corresponding to the phase matrices (in π/4). (c)–(f) Far-field patterns of the metasurfaces with phase distributions of a constant or a random number matrix.

    Figure 3.Phase design of the metasurface with circular dichroism-like effect consists of 140×140 units. (a) The phase matrices consist of equal elements or pseudorandom numbers (submatrix with equal elements). (b) The phase distributions corresponding to the phase matrices (in π/4). (c)–(f) Far-field patterns of the metasurfaces with phase distributions of a constant or a random number matrix.

    Figure 3(c) is the far-field pattern of the transmitted wave with a phase matrix of Φ0, where the incident wave only undergoes a phase shift after passing through the metasurface, and the transmission direction is almost unchanged. According to Figs. 3(d)–3(f), after the incident wave passes through the random phase matrix (Φ1,Φ2,Φ3), the transmitted wave is scattered in almost all directions, and the energy will be dispersed in a large angle range. When monitoring the transmitted wave at a certain point in the far field, only a small amount of energy will be observed. Then a CD-like response of the transmitted wave is obtained. In order to obtain a higher scattering efficiency, we compared the far-field patterns of several phase matrices with different repetition rates, that is, the phase matrix is composed of 1×1, 2×2, and 4×4 identical elements. It can be found that a larger number of repetitions cause a smaller scattering range and efficiency in the far field. Among them, Fig. 3(d) has the highest scattering efficiency.

    Morphology characterization and transmission spectrum of the metasurface. (a) Optical photograph and scanning electron microscope (SEM) images of sample 1. (b) Simulated two-dimensional electric field and phase distributions of the transmitted LCP and RCP waves. (c) Terahertz polarization measurement time-domain spectroscopy system. (d) Simulation results of the TCD spectrum for metasurfaces with different numbers of units. (e)–(j) Experiment results of the linear and circular polarization transmission coefficients and TCD spectra.

    Figure 4.Morphology characterization and transmission spectrum of the metasurface. (a) Optical photograph and scanning electron microscope (SEM) images of sample 1. (b) Simulated two-dimensional electric field and phase distributions of the transmitted LCP and RCP waves. (c) Terahertz polarization measurement time-domain spectroscopy system. (d) Simulation results of the TCD spectrum for metasurfaces with different numbers of units. (e)–(j) Experiment results of the linear and circular polarization transmission coefficients and TCD spectra.

    Figure 4(c) shows the experiment setup for the transmission spectrum measurement of sample 1. We placed four wire grid polarizers in a typical terahertz time-domain spectroscopy system to form a terahertz polarization test system [30]. Among them, the polarizers P1 and P3 are placed on the side close to the emitter and the receiver, respectively, to ensure that the polarization states of the transmitting and receiving waves remain unchanged during the measurement. In addition, P2 and P4 are rotated by ±45  deg (observed along the propagation direction of the terahertz beam) as the x and y coordinate axes of the measurement, as shown by the red and blue arrows in Fig. 4(c). Using the system in the figure, the four transmission coefficients for linear polarization components (txx, tyx, txy, and tyy) can be obtained. Then the transmission matrix of circularly polarized waves can be calculated as [31] Tcir=(trrtrltlrtll)=12(txx+tyy+i(txytyx)txxtyyi(txy+tyx)txxtyy+i(txy+tyx)txx+tyyi(txytyx)).We define the transmission difference of the two circular polarization components measured in a certain direction as TCD, which can be calculated by [19] TCD=TRTL=(|trr|2+|tlr|2)(|tll|2+|trl|2).Since our metasurface is a non-periodic structure, different numbers of cells may produce different circular polarization responses, thereby achieving tunable CD. In order to verify this function, we place the sample at different positions of the focused terahertz beam during the measurement, so that different numbers of units are illuminated by the terahertz wave. Before that, we first simulated the transmission coefficients of two metasurfaces with 60×60 elements and 120×120 elements, respectively, in which we placed the electric field probe directly in front of the center of the metasurface. Then we calculated the TCD spectrum in the range of 0.5–2.5 THz using Eq. (4), as shown in Fig. 4(d). It can be found that the maximum value of TCD and the curve width are slightly different in the two cases. The red curve reaches a maximum value of about 0.4 around 1.3 THz, and the peak range is narrow. The maximum value of the black curve is about 0.3, while the curve width is larger. It is worth mentioning that due to the reflection loss of high-resistance silicon, the maximum initial transmission coefficient is about 0.7. The maximum value of TCD is as high as about 0.4, showing a polarization conversion efficiency of about 80%. In order to experimentally verify the circular dichroism of the sample, we use the experiment setup shown in Fig. 4(c) to perform the polarization spectrum measurement of sample 1. First, we place the sample at the focal plane (position A, d=0  cm) of the off-axis parabolic lens, and the measured linear and circular polarization spectra are shown in Figs. 4(e)–4(g). It can be seen in the circularly polarized transmission spectrum that the values of the two co-polarized components (trr and tll) are relatively small. For the cross-polarization components, trl is as high as about 0.55, but tlr is close to zero. Therefore, the calculated TCD value in Fig. 4(g) reaches about 0.3 near 1.4 THz. Then, we moved the sample about 4 cm along the incident direction and placed it in position B for a similar measurement. The results can be seen in Figs. 4(h)–4(j). The diameter of the terahertz beam at position B (d=4  cm) is larger than that at position A, so there will be more units in the sample working. It can be found that the TCD curve in Fig. 4(j) is narrower than that in Fig. 4(g), and the maximum value reaches 0.4 near 1.4 THz. This change rule is consistent with the simulation results in Fig. 4(d), and the tunable CD is obtained. The slight deviation of the working frequency may be caused by the processing deviation of the sample.

    Spin-selective wavefront control of transmitted terahertz wave using the proposed metasurface. (a) Optical photograph and SEM images of sample 2. (b) Experimental setup of sample 2. (c)–(j) Simulations and measured results of the two-dimensional electric field and phase distributions.

    Figure 5.Spin-selective wavefront control of transmitted terahertz wave using the proposed metasurface. (a) Optical photograph and SEM images of sample 2. (b) Experimental setup of sample 2. (c)–(j) Simulations and measured results of the two-dimensional electric field and phase distributions.

    We use the terahertz imaging system in Fig. 5(b) to measure the electric field and phase of the transmitted wave. The beam output from the amplified femtosecond laser system (the pulse width is 50 fs, repetition frequency is 1 kHz, and center wavelength is 800 nm) is divided into two paths, which are used as the pump beam and probe beam for terahertz imaging. The terahertz wave generated by the ZnTe crystal on the left is vertically polarized (y axis). After passing through the metasurface, the THz wave is irradiated on the crystal on the right and the probe laser beam experiences an electro-optic modulation effect. Then the modulated laser transmits the terahertz electric field image to the CCD camera through the optical imaging system (lenses L2 and L3, quarter-wave plate, and Wollaston prism). During the measurement, the incident and detected terahertz polarization is changed by rotating the sample or the optical polarizer and the half-wave plate at the pump end, and the terahertz electric fields Exx, Eyx, Eyy, Exy of the four linear polarization components are observed; then the circularly polarized terahertz electric field of the sample can be obtained by Eq. (3).

    Before measuring the sample 2, we simulated the transmitted electric field of it with commercial software. Figures 5(c) and 5(d) show the electric fields of the RCP and LCP waves at the focal plane, and Figs. 5(e) and 5(f) show the phase distributions. It can be seen that the RCP component is a vortex beam, while the LCP component is a messy scattered light field. The experimental results at z=6000  μm are shown in Figs. 5(g)–5(j), with a working frequency of 1.38 THz. It can be seen that the measured electric field and phase are basically consistent with the simulation results. It should be pointed out that due to the limitation of the aspect ratio in the etching process, we select the designed etching depth as 150 μm, which results in our optimal design working frequency of about 1.3 THz. More importantly, the actual operating frequency is 1.38 THz due to the processing deviation. The terahertz imaging system we use shows a low signal-to-noise ratio (SNR) at 1.38 THz, which led to the deviation in the results of Fig. 5. Fortunately, we believe that the SNR is still acceptable. This may also show that our scheme has good fabrication tolerance.

    3. CONCLUSION

    In summary, we propose a new method to achieve giant chiral response of transmitted terahertz waves based on phase control. By introducing geometric phase and dynamic phase at the same time, we designed a spin-decoupled all-silicon terahertz metasurface, which can set different phase matrices for LCP and RCP incident waves. We analyzed the required conditions of the spin-decoupled metasurface elements, and then simulated and selected their geometric parameters. Using the equivalent phase matrix and the random number phase matrix, we achieve the CD-like effect of the transmitted terahertz wave. Simulated and measured results show that the maximum polarization conversion efficiency is about 80%, and the maximum CD value is 0.4. The CD value can be adjusted by the relative position of the sample in the focused terahertz beam. In addition, we also show that this method is used to achieve spin-dependent arbitrary wavefront control. Measurement results of the two-dimensional terahertz electric field are consistent with the simulation results. The proposed method does not require complex asymmetric chiral units, but only needs to set the phase matrix of the metasurface, which is expected to be used in the design of new chiral devices for terahertz wave.

    References

    [1] L. D. Barron. True and false chirality and parity violation. Chem. Phys. Lett., 123, 423-427(1986).

    [2] J. Mun, M. Kim, Y. Yang, T. Badloe, J. Ni, Y. Chen, C. W. Qiu, J. Rho. Electromagnetic chirality: from fundamentals to nontraditional chiroptical phenomena. Light Sci. Appl., 9, 139(2020).

    [3] Z. Cao, H. Gao, M. Qiu, W. Jin, S. Deng, K. Wong, D. Lei. Chirality transfer from sub-nanometer biochemical molecules to sub-micrometer plasmonic metastructures: physiochemical mechanisms, biosensing, and bioimaging opportunities. Adv. Mater., 32, 1907151(2020).

    [4] Y. Tang, A. E. Cohen. Optical chirality and its interaction with matter. Phys. Rev. Lett., 104, 163901(2010).

    [5] Y. Tang, A. E. Cohen. Enhanced enantioselectivity in excitation of chiral molecules by superchiral light. Science, 332, 333-336(2011).

    [6] N. Yu, P. Genevet, M. A. Kats, F. Aieta, J. Tetienne, F. Capasso, Z. Gaburro. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science, 334, 333-337(2011).

    [7] J. K. Gansel, M. Thiel, M. S. Rill, M. Decker, K. Bade, V. Saile, G. Freymann, S. Linden, M. Wegener. Gold helix photonic metamaterial as broadband circular polarizer. Science, 325, 1513-1515(2009).

    [8] A. Kuzyk, R. Schreiber, Z. Fan, G. Pardatscher, E. Roller, A. Hogele, F. Simmel, A. Govorov, T. Liedl. DNA-based self-assembly of chiral plasmonic nanostructures with tailored optical response. Nature, 483, 311-314(2012).

    [9] H. Xu, G. Hu, Y. Li, L. Han, J. Zhao, Y. Sun, F. Yuan, G. Wang, Z. Jiang, X. Ling, T. Cui, C. Qiu. Interference-assisted kaleidoscopic meta-plexer for arbitrary spin-wavefront manipulation. Light Sci. Appl., 8, 3(2019).

    [10] W. Li, Z. J. Coppens, L. V. Besteiro, W. Wang, A. O. Govorov, J. Valentine. Circularly polarized light detection with hot electrons in chiral plasmonic metamaterials. Nat. Commun., 6, 8379(2015).

    [11] L. Kang, S. P. Rodrigues, M. Taghinejad, S. Lan, K. Lee, Y. Liu, D. H. Werner, A. Urbas, W. Cai. Preserving spin states upon reflection: linear and nonlinear responses of a chiral meta-mirror. Nano Lett., 17, 7102-7109(2017).

    [12] T. Kim, S. Oh, H. Kim, H. Park, O. Hess, B. Min, S. Zhang. Electrical access to critical coupling of circularly polarized waves in graphene chiral metamaterials. Sci. Adv., 3, e1701377(2017).

    [13] J. Guirado, M. Svedendahl, J. Puigdollers. Enhanced chiral sensing with dielectric nanoresonators. Nano Lett., 20, 585-591(2020).

    [14] L. D. Barron, N. Gadegaard, M. Kadodwala. Ultrasensitive detection and characterization of biomolecules using superchiral fields. Nat. Nanotechnol., 5, 783-787(2010).

    [15] W. Zhang, T. Wu, R. Wang, X. Zhang. Amplification of the molecular chiroptical effect by low-loss dielectric nanoantennas. Nanoscale, 9, 5701-5707(2017).

    [16] E. Plum, V. A. Fedotov, N. I. Zheludev. Metamaterials: optical activity without chirality. Phys. Rev. Lett., 102, 113902(2009).

    [17] I. Sersic, M. A. Haar, F. B. Arango, A. F. Koenderink. Ubiquity of optical activity in planar metamaterial scatterers. Phys. Rev. Lett., 108, 223903(2012).

    [18] L. Jing, Z. Wang, R. Maturi, B. Zheng, H. Wang, Y. Yang, L. Shen, R. Hao, W. Yin, E. Li, H. Chen. Gradient chiral metamirrors for spin-selective anomalous reflection. Laser Photonics Rev., 11, 1700115(2017).

    [19] Z. Wang, H. Jia, K. Yao, W. Cai, H. Chen, Y. Liu. Circular dichroism metamirrors with near-perfect extinction. ACS Photonics, 3, 2096-2101(2016).

    [20] F. Neubrech, M. Hentschel, N. Liu. Reconfigurable plasmonic chirality: fundamentals and applications. Adv. Mater., 32, 1905640(2020).

    [21] D. C. Hooper, A. G. Mark, C. Kuppe, J. T. Collins, P. Fischer, V. K. Valev. Strong rotational anisotropies affect nonlinear chiral metamaterials. Adv. Mater., 29, 1605110(2017).

    [22] C. Niu, Z. Wang, J. Zhao, L. Du, N. Liu, Y. Liu, X. Li. Photonic heterostructures for spin-flipped beam splitting. Phys. Rev. Appl., 12, 044009(2019).

    [23] J. Fan, Y. Cheng, B. He. High-efficiency ultrathin terahertz geometric metasurface for full-space wavefront manipulation at two frequencies. J. Phys. D, 54, 115101(2021).

    [24] Y. Cheng, F. Chen, H. Luo. Multi-band giant circular dichroism based on conjugated bilayer twisted-semicircle nanostructure at optical frequency. Phys. Lett. A, 384, 126398(2020).

    [25] Y. Cheng, J. Fan, H. Luo, F. Chen. Dual-band and high-efficiency circular polarization convertor based on anisotropic metamaterial. IEEE Access, 8, 7615-7621(2020).

    [26] J. Fan, Y. Cheng. Broadband high-efficiency cross-polarization conversion and multi-functional wavefront manipulation based on chiral structure metasurface for terahertz wave. J. Phys. D, 53, 025109(2020).

    [27] Y. Li, Y. Li, L. Chen, M. Hong. Reflection tuning via destructive interference in metasurface. Opto-Electron. Eng., 44, 313-318(2017).

    [28] J. P. B. Mueller, N. A. Rubin, R. C. Devlin, B. Groever, F. Capasso. Metasurface polarization optics: independent phase control of arbitrary orthogonal states of polarization. Phys. Rev. Lett., 118, 113901(2017).

    [29] T. J. Cui, M. Q. Qi, X. Wan, J. Zhao, Q. Cheng. Coding metamaterials, digital metamaterials and programmable metamaterials. Light Sci. Appl., 3, e218(2014).

    [30] W. Choi, G. Cheng, Z. Huang, S. Zhang, T. B. Norris, N. A. Kotov. Terahertz circular dichroism spectroscopy of biomaterials enabled by kirigami polarization modulators. Nat. Mater., 18, 820-826(2019).

    [31] C. Menzel, C. Rockstuhl, F. Lederer. Advanced Jones calculus for the classification of periodic metamaterials. Phys. Rev. A, 82, 053811(2010).

    Jie Li, Chenglong Zheng, Guocui Wang, Jitao Li, Hongliang Zhao, Yue Yang, Zhang Zhang, Maosheng Yang, Liang Wu, Jining Li, Yating Zhang, Yan Zhang, Jianquan Yao, "Circular dichroism-like response of terahertz wave caused by phase manipulation via all-silicon metasurface," Photonics Res. 9, 567 (2021)
    Download Citation