Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Multi-functional high-efficiency reflective polarization converter based on an ultra-thin graphene metasurface in the THz band

Open Access Open Access

Abstract

In this study, an ultra-thin reflective metasurface is proposed for polarization conversion in the terahertz band. Each unit cell of metasurface is composed of graphene ribbons lying diagonally on silicon substrate. A reflective metal is also placed at the bottom of the structure. Our polarization converter works as a linear polarization converter (LPC) and linear to circular polarization converter (LTC-PC) by variation of the chemical potential of graphene, which can actively be changed by chemical doping or electrical bias of the graphene. The working bandwidth of LPC changes by adjusting the chemical potential of the graphene. The LPC structure has more than 99% polarization conversion ratio in the frequency range of 0.83-0.92 THz, even by changing the angle of incident wave up to 45°, the results are still acceptable. The LTC-PC has less than 3dB axial ratio (AR) in the frequency range of 0.6-0.67 THz for left-handed circularly polarized (LHCP) waves and 0.72-0.97 THz for right-handed circularly polarized (RHCP) waves. To verify the simulation results, an equivalent circuit model based on the structure performance is proposed. Equivalent circuit model results agree very well with the simulation results. Due to the fabrication feasibility, ultra-thin thickness, incident angle insensitive, and high efficiency, our structure has great potential in state-of-the-art technologies such as imaging, sensing, communication, and other optical applications.

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

1. Introduction

Terahertz (THz) region, the frequency range of 0.1 THz to 10 THz, has been unexplored for many decades. Recently, THz regime has attracted a great interest in detection, imaging, sensing, spectroscopy, and security devices [14]. Lack of natural materials which can interact with THz waves has made metasurfaces a good choice for practical terahertz devices [5]. Metasurface, a two-dimensional artificial material which can be an appropriate replacement for bulky metamaterials, consists of periodic subwavelength metal or dielectric patterns. By changing shape, size, and material of the unit cells, different optical properties, such as anomalous optical reflection and negative refractive indices, are achievable [6,7]. Graphene, a single layer of carbon atoms arranged in a honeycomb lattice, shows unique carrier mobility, mechanical, electrical, thermal, electromagnetic, and optical properties [810]. In spite of having a single atom thickness, graphene exhibits strong interaction with incident wave and supports surface plasmon polaritons (SPPs) with long propagation length from THz to mid-infrared (MIR), therefore, it is a good choice to be integrated into metamaterials and metasurfaces [9,1114]. Controlling the chemical doping and electrical bias of graphene lead to a tunable surface conductivity of graphene which affects its electrical and optical response [7,12]. In addition, graphene is transparent in THz region [9,11]. Recently, graphene has been used in many optical and electrical devices such as electro-optic modulators, absorbers, couplers, polarization converters, solar cells, and transistors [12,1417]. Periodic structure enhances the energy coupled between incident wave and SPPs of graphene [9]. The graphene metasurfaces introduce actively tunable applications which their properties can change not only by size and shape, but also by graphene conductivity. This means that we would have a reconfigurable structure which does not need to be refabricated.

Many practical devices are sensitive to polarization, therefore, manipulation of the polarization state of electromagnetic waves has got great importance [18]. Polarization converters (PCs) are capable of using in imaging, sensing, detection and communication applications including modulators and switches [8,1922]. Conventional polarization converters use birefringent crystals or optical gratings to change the polarization state of the waves [18,22,23]. These PCs are not appropriate for THz regime because of their bulky structures, low efficiency and narrow bandwidth [6,8]. Recently, metamaterials and metasurfaces have suppressed this limitation and many structures have been proposed based on metamaterials and metasurfaces [14,18,2428]. Metamaterials are still bulky and are not appropriate to be used in the state-of-the-art applications. Although metasurfaces eliminate this problem, many of them have narrow bandwidth, low efficiency, complicated structure which cannot be easily fabricated and are not tunable [14,24,27,28].

In this study, we propose a high-efficiency linear polarization converter (LPC) with graphene metasurface for THz regime. It is shown that by changing the graphene chemical potential, $\mu _c$, of the proposed LPC structure, linear to circular (LTC) polarization conversion is also achievable. To verify the proposed structure performance and to find a fast method to get the analysis results, an equivalent circuit model is proposed. Results of equivalent circuit model (ECM) and simulations are in very good agreement. Our proposed structure is ultra-thin and can be easily fabricated. By using graphene in the structure, the optical response and application of the metasurface can be varied by changing the graphene chemical potential, ${\mu _c}$, without re-fabricating the structure.

The remaining of this paper is organized as follows: In Section 2, the structure design and analysis method are described. The results are presented and discussed in Section 3. In Section 4 an equivalent circuit model is proposed for the structure. Finally, the paper is concluded in Section 5.

2. Structure design and analysis method

The schematic view of the proposed three-dimensional reflective polarization converter (PC) structure and its unit cell are shown in Fig. 1. Our proposed structure is composed of three layers, which from bottom to top are gold layer, silicon substrate and graphene layer, respectively. The bottom layer is a reflective gold which reflects all the incident wave and plays an important role in obtaining efficient polarization conversion. The silicon substrate with refractive index of 3.4 is deposited over the gold layer. Graphene layer is composed of a diagonal graphene ribbon in the middle of the unit cell and two small graphene triangles at corners of the unit cell as shown in Fig. 1(b). These two small triangles make the graphene ribbons continuous in the metasurface, so the PC can be easily fabricated. We believe that the fabrication of our proposed structure is feasible and straightforward. Many of metasurfaces have spirals, double rings, ginkgo-leafs, and other complicated geometrical patterns which are not easy for lithography. On the other hand, the graphene ribbons have simple and common structure in comparison with those complicated structures. To pattern a large-area graphene, the etching method has been widely employed by using a mask consisting of a resist pattern defined by electron beam lithography, pre-deposited nanowires, or self-assembled block copolymer (BCP) nanopatterns [29]. Graphene ribbons can also be written by lithography without a mask [29]. Also, there are plenty of methods for graphene ribbons fabrication [29].

 figure: Fig. 1.

Fig. 1. (a). The schematic view of the proposed three dimensional linear polarization converter (LPC) and linear to circular polarization converter (LTC- PC) and (b) the structure unit cell.

Download Full Size | PDF

Graphene is modeled as a two-dimensional surface and its surface conductivity is calculated by Kubo formula, which is consisted of inter-band and intra-band conductivities [30]:

$${{\sigma}}={{\sigma}_{intra}} + {{\sigma}_{inter}}$$

Intra-band conductivity can be evaluated as:

$${{\sigma}_{intra}}=\frac{-j{e^2}{{k}_{B}}T}{\pi{\hbar^2}(\omega-j2\Gamma)} \left[ \frac{{{\mu}_{c}}}{{{k}_{B}}T} +2 \ln(e^ \frac{ -{{\mu}_{c}}}{{{k}_{B}}T} +1)\right]$$
where $ \omega$, $ e$, $ \hbar$, $ k_B$, $\Gamma$, $ \mu _c$, and $ T$ are respectively angular frequency, electron charge, reduced Planck’s constant, Boltzmann constant, scattering rate, graphene chemical potential, and room temperature.

In special conditions, when $ {k_B T \ll |\mu _c|, \hbar \omega }$ inter-band can be approximated as [30]:

$${{\sigma}_{inter}}=\frac{-j{e^2}}{4\pi \hbar} \ln \frac{2|{\mu}_{c}|-(\omega-j2\Gamma)\hbar}{2|{\mu}_{c}|+(\omega-j2\Gamma)\hbar}$$

Some parameters of equations and the optimum dimensions for maximum efficiency of our proposed structure of Fig. 1 are given in Table 1. As can be seen, the optimum thickness of substrate is large, which would require a large bias voltage, $V_g$, to change the $\mu _c$, for excitation of graphene ribbons. As shown in Fig. 1, to address this issue and get a practical structure with applicable bias voltage, a thin layer of a transparent conductive oxide (TCO) can be applied between the ${t}_{s}$ and rest of the substrate, which due to the extreme thinness of the TCO, its effect in simulation is negligible [31,32]. Indium tin oxide (ITO), fluorine doped tin oxide (FTO), doped zinc oxide or networks of polymers such as polysilicon can be used as TCO.

Tables Icon

Table 1. The optimum dimensions for maximum efficiency of our proposed structure of Fig. 1; and the values of some parameters of Eqs. (1)–(3)

The relation between ${\mu }_{c}$ and $V_g$ can be approximated to the below expression [31,3335] :

$${{\mu}_{c}}\simeq\hbar {\nu}_{f} \sqrt{\frac{\pi \epsilon_{r_{Si}}\epsilon_{0} {V}_{g}}{e{t}_{s}}}$$
where ${\epsilon }_{r_{Si}}$, ${\epsilon }_{0}$, and $\nu _f$ are the relative permittivity of substrate, permittivity of vacuum, and the Fermi velocity ($1.1 \times 10^6$ m/s in graphene), respectively. Therefore, the maximum of bias voltage to bias graphene ribbons for ${t}_{s}$=2 ${\mu } m$ will be 47 V which is a reasonable value.

The PC structure is analyzed by the finite difference time domain (FDTD) method. The boundary conditions are considered periodic in $ x$ and $ y$ directions and perfectly matched layer (PML) in $ z$ axis. The meshes in $ x$ and $ y$ directions on the graphene layer are considered 0.5 $\mu$m. The PC is illuminated by an $ x$-polarized normal plane wave from the top. When the PC is illuminated by an $ x$-polarized wave, a plasmonic resonance is excited along the graphene ribbon which leads to $ x$ and $ y$-polarized reflected wave.

To describe the reflected wave and its polarization state, we introduce the Stokes Parameters as follows [3638]:

$${{I}_{n}}={|{R}_{xx}|^2}+{|{R}_{yx}|^2}$$
$${{Q}_{n}}={|{R}_{xx}|^2}-{|{R}_{yx}|^2}$$
$${{U}_{n}}=2{|{R}_{xx}|}{|{R}_{yx}|}cos(\Delta\phi)$$
$${{V}_{n}}=2{|{R}_{xx}|}{|{R}_{yx}|}sin(\Delta\phi)$$
where $ {R_{xx}}$ and $ {R_{yx}}$ are respectively the amplitudes of the co-polarized and cross-polarized reflection coefficients when the structure is illuminated by a normal $ x$-polarized wave. ${\Delta \phi }$ is the difference of corresponding phases $ {\phi _{xx}}$ and $ {\phi _{yx}}$. The Stokes Parameters $ {I_n}$, $ {Q_n}$, $ {U_n}$, and $ {V_n}$ stand respectively for total reflection, the horizontal and vertical linear polarization state of the reflection, the linear +45 or -45 polarization state of the reflection, and circular polarization state. The subscript $ n$ denotes that the Stokes Parameters are for reflections, i.e. are normalized to incident wave. The axial ratio (AR) is defined as [36]:
$${AR}=10log(\tan \beta)$$
where $ \beta$ is the ellipticity angle defined as [36]:
$${\beta}=\frac{1}{2}{\sin^{{-}1}\left(\frac{{V}_{n}}{{I}_{n}}\right)}$$

The axial ratio (AR) is used for evaluating the circular polarization performance. The analyzed wave is circularly polarized if the axial ratio is less than 3dB.

3. Results and discussion

Figure 2(a) illustrates the spectra of the reflection coefficients when the LPC structure is illuminated by a normal $ x$-polarized wave with central frequency of 0.83 THz (360 $\mu$m) amplitudes in $ z$- direction. As shown, most of the reflection is cross-polarized, while co-polarized reflected wave is negligible in the frequency range of 0.78-0.97 THz. $ {Q_n}$ parameter of Eq. (6) is depicted in Fig. 2(b). In the frequency range of 0.78 THz to 0.97 THz, $ {Q_n}$ is approximately equal to -1, which stands for linearly $ y$-polarized reflected wave, and when the reflected wave is mostly $ x$-polarized, $ {Q_n}$ is around +1. To approve LPC’s performance, polarization conversion ratio (PCR) is calculated as [7,39]:

$${PCR}=\frac{{{R}_{yx}^2}}{{{R}_{xx}^2}+{{R}_{yx}^2}}$$

 figure: Fig. 2.

Fig. 2. The spectra of (a) Co-polarized and cross-polarized reflection coefficient amplitudes, (b) Qn parameter, and (c) the PCR of the proposed structure of Fig. 1, as a linear polarization converter (LPC).

Download Full Size | PDF

As shown in Fig. 2(c) PCR is more than 0.9 in the frequency range of 0.78-0.97 THz (308-384 $\mu$m), which implies that linear polarization is properly converted to its cross polarized wave, which the device performs as a linear polarization converter (LPC).

Figure 3(a) demonstrates the spectra of the reflection coefficients amplitudes for different height of the structure. The related PCR spectra are depicted in Fig. 3(b). As shown, PCR is not very ideal for $ {h}$=20 $\mu$m, but for $ {h}$=23 $\mu$m or 26 $\mu$m PCR is ideal with a little difference in bandwidth. Hence, $ {h}$=23 $\mu$m is chosen in this study.

 figure: Fig. 3.

Fig. 3. The spectra of (a) Reflection coefficient amplitude for three different heights of the substrate and (b) the related PCR of the structure of Fig. 1.

Download Full Size | PDF

In practical applications, the transmitter and receiver block each other in normal incidence. Hence, oblique incidence of plane wave must be considered. For this purpose, we have simulated LPC structure for different angles of incidence of $ {\theta }$ = $15^{\circ }$, $30^{\circ }$, and $45^{\circ }$. The related PCRs are depicted in Fig. 4. As shown, even for the wave with $45^{\circ }$ incident angle PCR is more than 90% in a wide frequency range of 0.77 THz to 0.97 THz.

 figure: Fig. 4.

Fig. 4. The spectra of PCR of structure of Fig. 1 for different angle of incidence.

Download Full Size | PDF

The spectra of PCR for different $ {\mu _c}$ of graphene is demonstrated in Fig. 5. As shown, altering $ {\mu _c}$ from 0.4 eV to 0.7 eV, changes the LPC’s working frequency range.

 figure: Fig. 5.

Fig. 5. The LPC spectra of PCR of structure of Fig. 1 for different ${\mu _c}$ of the graphene.

Download Full Size | PDF

The cross-polarized reflection coefficient amplitude versus $ {\mu _c}$ at frequency of 0.83 THz (360 $\mu m$) is depicted in Fig. 6. As illustrated, for $ {\mu _c}$ more than 0.25, the structure behaves as linear polarization converter. Therefore, our proposed LPC is an active metasurface and one can find the best performance for the desired frequency range just by adjusting $ {\mu _c}$. Figure 6 indicates that when $ {\mu _c}$ changes in the range of 0.14-0.17 eV, the reflection coefficient amplitude of the cross-polarized wave is around 0.7 (3 dB), which means that the amplitude of the reflected $ {y}$-polarized wave is equal to the $ {x}$-polarized one. Thus, a circularly-polarized wave is achievable.

Figure 7 illustrates the spectra of the reflection coefficients amplitude and phases of the co-polarized and cross-polarized electrical fields reflected in $ {z}$ direction when the graphene chemical potential is set as 0.152 meV and the structure is illuminated by a normal $ {x}$-polarized wave with central frequency of 0.83 THz. As shown in Fig. 7(a), normalized reflection coefficient amplitudes of cross-polarized and co-polarized reflected waves at the frequency range of 0.61 THz to 0.67 THz (wavelength range of 457-492 $\mu m$) and 0.72-0.97 THz (wavelength range of 308-413 $\mu m$) are between 0.5 and 0.8 which can be considered approximately equal. As depicted in Fig. 7(b), in these ranges the cross-polarized and co-polarized reflected waves have $\pm 90^\circ$ phase difference. The equal amplitudes and $\pm 90^\circ$ phase differences approve that the reflected wave is a circularly polarized wave.

 figure: Fig. 6.

Fig. 6. Variation of the cross-polarized reflection coefficient amplitude of structure of Fig. 1 vs. chemical potential, $ {\mu _c}$, at the central frequency of 0.83 THz.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. The spectra of (a) amplitude and (b) phases of the reflection coefficients of the structure of Fig. 1.

Download Full Size | PDF

Figure 8(a) illustrates the spectra of AR in desired frequency range. As shown, AR is less than 3dB in the frequency range of 0.6-0.67 THz (wavelength range of 457-492 $\mu m$) and 0.72-0.97 THz (wavelength range of 308-413 $\mu m$), which indicates that LTC polarization conversion has properly happened. The ellipticity of reflected wave is calculated as ${{V}_{n}}/{{I}_{n}}$ which is depicted in Fig. 8(b). As shown in Fig. 8(b), in the frequency range of 0.72-0.97 THz ellipticity is around +1, which indicates that the reflected wave is right-handed circularly polarized (RHCP) wave. Ellipticity is around -1 in the range of 0.6 to 0.67 which stands for a left-handed circularly polarized (LHCP) wave.

 figure: Fig. 8.

Fig. 8. The spectra of (a) axial ratio and (b) the ellipticity for proposed LTC-PC structure of Fig. 1.

Download Full Size | PDF

Based on the results of the simulations, the proposed structure of Fig. 1 can be used as linear polarization converter (LPC) and also linear to circular polarization converter (LTC-PC) by changing the graphene chemical potential, $\mu _c$. Essential dimensions of the LPC and LTC-PC are given in Table 1. The dimensions have been optimized for obtaining maximum efficiency. By setting $\mu _c$ as 0.152 eV, our proposed structure works as LTC-PC which converts linearly polarized incident wave to a circularly polarized wave. For $\mu _c$ more than 0.25 eV the structure works as LPC and converts the linearly polarized incident wave ($ {x}$- or $ {y}$-polarized) to its cross polarized wave ($ {y}$- or $ {x}$-polarized). The best performance of LPC is achievable at $\mu _c$ = 0.5 eV. We have investigated that our proposed PC structure is more efficient and feasible to fabricate, compared to other similar structures reported in Table 2.

Tables Icon

Table 2. Comparison of our proposed structure and some other similar structures

4. Equivalent circuit model (ECM)

To verify simulation results, due to the perforxmance of the PC structure of Fig. 1, we have proposed an equivalent circuit model (ECM). This equivalent circuit is shown in Fig. 9. Impedance of the silicon substrate is considered as ${{Z}_{Si}}={{\eta }_{0}}/{\sqrt {{\epsilon }_{{r}_{Si}}}}$, where ${{\eta }_{0}}\cong 120\pi$ which is the free space wave impedance. Gold layer behaves as a perfect electric conductor (PEC) because of its short skin depth in THz range, and therefore, its impedance value is about zero. As it can be seen in Fig. 1(a), graphene ribbons are lying on silicon substrate with 45-degree angle with $ {x}$ and $ {y}$ axis. The structure is illuminated by an $ {x}$-polarized normal plane wave which can be decomposed into $ {x^\prime }$ and $ {y^\prime }$ axis as:

$$\textbf{E}_{\textbf{inc}}={E}_{0x}\textbf{a}_{\textbf{x}}=\frac{\sqrt{2}}{2}{E}_{0x}(\textbf{a}_{\textbf{x}^{\prime}}+\textbf{a}_{\textbf{y}^{\prime}})$$

 figure: Fig. 9.

Fig. 9. The equivalent circuit model for the structure of Fig. 1, in (a) $x^{\prime }$ and (b) $y^{\prime }$ directions.

Download Full Size | PDF

We have considered two equivalent circuit models for $ {x^\prime }$ and $ {y^\prime }$ axis, individually. In $ {y^\prime }$ direction, periodic arrays of graphene ribbons (PAGR) can be modeled with parallel R-L-C circuits. The relative values of the resistance, inductance, and capacitance of the R-L-C are obtained as [41]:

$${R}_{GR}=\frac{P}{{{S}_{1}^2}}\frac{\pi {\hbar^2}}{{e^2}\mu_{c}\tau}$$
$${L}_{GR}=\frac{P}{{{S}_{1}^2}}\frac{\pi {\hbar^2}}{{e^2}\mu_{c}}$$
$${C}_{GR}=\frac{{{S}_{1}^2}}{P}\frac{2 \epsilon_{eff}}{{q}_{1}}$$
where subscript $ GR$ stands for graphene, ${{{S}_{1}^2}}={8a}/{9}$, in which $ a$ is the width of the ribbons; and $q_1$ is the eigenvalue of the equation which is given in Table 1 of [41], $q_1$ is the first eigenvalue of the equation governing the current on the graphene sheet. $P$ and $\tau$ are respectively the period of the unit cell, and the carrier relaxation time. The effective permittivity of graphene’s surrounding area is calculated by ${{\epsilon }_{eff}}={({\epsilon }_{0}+{\epsilon }_{Si})}/{2}$, in which $\epsilon _0$ and $\epsilon _{Si}$ are the permittivity of the air and silicon, respectively. The relative permittivity of silicon is 3.4. The impedance of each graphene ribbon can be calculated as:
$${Z}_{G1}={{R}_{GR}}+j{\left(\omega{{L}_{GR}}-\frac{1}{\omega{{C}_{GR}}}\right)}$$

In $ x^\prime$ direction, graphene is modeled as two parallel inductances, each representing the magnetic energy stored in the upper and lower media, as [41]:

$${L}_{i}=\frac{{{\eta}_{0}}P}{c\pi}\ln{\left[\csc{\left(\frac{\pi a}{2P}\right)}\right]}{\left(1-\alpha+\frac{\alpha}{\sqrt{1-{({\lambda}_{ci}/{\lambda})^2}}}\right)}, {i=1,2}$$
where $c$ and $\lambda$ are respectively speed of light and wavelength of incident wave. The incident wavelength is considered as 360 $\mu m$, ${\lambda }_{ci}={n}_{i}{(1+\sin {\theta _{i}})}P$, and ${\alpha }\approx \sqrt {1-[(P-a)/P]^4}$, $n_1$ and $n_2$ are refractive indices of media on top and bottom of graphene, which are equal to $n_0$ and $n_{Si}$, respectively. $\theta _1$ and $\theta _2$ are angle of incident and reflected wave respectively, both are considered to be zero. In this condition, the impedance of PAGR is calculated as [42]:
$${Z}_{G2}= j\omega\frac{L_1L_2}{L_1+L_2}+\frac{P}{\sigma a}$$
where $\sigma$ is the surface conductivity of graphene. By neglecting the loss of the structure, the equivalent input impedance of silicon, terminated with gold sheet can be calculated as:
$${Z}_{t}= Z_{Si}\frac{Z_{Au}+jZ_{Si}\tan(\beta_{Si} h)}{Z_{Si}+jZ_{Au}\tan(\beta_{Si} h)}$$
where $\beta _{Si}$ is the propagation constant of the THz wave in the silicon substrate. The height of the silicon substrate, $h$, is assumed to be 23$\mu m$. The equivalent impedances $Z_{T1}$ and $Z_{T2}$ depicted in Fig. 9(a) and 9(b), are calculated by:
$${Z}_{T1}= \frac{Z_{G1}Z_{t}}{Z_{G1}+Z_{t}}$$
$${Z}_{T2}= \frac{Z_{G2}Z_{t}}{Z_{G2}+Z_{t}}$$
Using these total equivalent impedances, $Z_{T1}$ and $Z_{T2}$, the reflection coefficient of the structure in $ x^\prime$ and $ y^\prime$ directions can be obtained by:
$${R}_{y^{\prime}}= \frac{Z_{T1}-\eta_{0}}{Z_{T1}+\eta_{0}}$$
$${R}_{x^{\prime}}= \frac{Z_{T2}-\eta_{0}}{Z_{T2}+\eta_{0}}$$

The reflection coefficients are calculated by using equivalent circuit model and FDTD method which are illustrated in Fig. 10(a) and 10(b), respectively, for LPC and LTC-PC structure. The results of these two methods are in good agreement. As depicted in Fig. 10(a), equivalent circuit model of cross polarized reflection coefficient is approximately equal compared with the negligible co-polarized reflected wave, in the frequency range of 0.81-1 THz. Thus, we can either realize that linear polarization conversion is done properly using equivalent circuit model. Also, reflection coefficient amplitudes obtained by equivalent circuit model in Fig. 10(b), show that amplitude of co-polarized and cross-polarized reflection coefficients are approximately the same, in frequency range of 0.75-0.97 THz. This means that linear to circular polarization conversion could be done in this range. Equivalent circuit model and FDTD results have slight difference in both LPC and LTC-PC structures, which have occurred due to the approximations in ECM and different analysis methods used in ECM and FDTD methods.

 figure: Fig. 10.

Fig. 10. Reflection coefficient amplitudes calculated by equivalent circuit model theory (ECM) and the FDTD method for (a). LPC structure and (b). LTC-PC structure of Fig. 1.

Download Full Size | PDF

5. Conclusion

In this study, we have proposed a multi-functional ultra-thin high efficiency reflective polarization converter (PC), using active graphene metasurface in terahertz (THz) band. Our proposed structure can work as linear polarization converter (LPC) and linear to circular polarization converter (LTC-PC) by changing the chemical potential of graphene. Both PCs are simulated by finite difference time domain method, and also equivalent circuit models are proposed for both structures. Results of these two methods were in a good agreement. A perfect linear polarization conversion was achieved with polarization conversion ratio (PCR) higher than 99% in the frequency range of 0.83–0.92 THz. LTC-PC is achieved in two different wavelength ranges. In the frequency range of 0.6 to 0.67 left-handed circularly polarized (LHCP) wave is achieved and in the frequency range of 0.72-0.97 THz the right-handed circularly polarized (RHCP) wave is achieved. Our proposed active PC structure is composed of graphene nanoribbons over silicon substrate which its fabrication is simple. This property in addition to its ultra-thin thickness and high efficiency are highly desirable properties in THz band to give this PC structure great application potentials in imaging, sensing, and communications.

Disclosures

The authors declare that there is no conflict of interest.

Data availability

No data were generated or analyzed in the presented research.

References

1. C. Jansen, S. Wietzke, O. Peters, M. Scheller, N. Vieweg, M. Salhi, N. Krumbholz, C. Jördens, T. Hochrein, and M. Koch, “Terahertz imaging: applications and perspectives,” Appl. Opt. 49(19), E48–E57 (2010). [CrossRef]  

2. P. U. Jepsen, R. H. Jacobsen, and S. Keiding, “Generation and detection of terahertz pulses from biased semiconductor antennas,” J. Opt. Soc. Am. B 13(11), 2424–2436 (1996). [CrossRef]  

3. L. Qi, “Broadband and polarization independent terahertz metamaterial filters using metal-dielectric-metal complementary ring structure,” J. Opt. Soc. Korea 20(2), 263–268 (2016). [CrossRef]  

4. L. Ye, F. Zeng, Y. Zhang, and Q. H. Liu, “Composite graphene-metal microstructures for enhanced multiband absorption covering the entire terahertz range,” Carbon 148, 317–325 (2019). [CrossRef]  

5. G. Deng, J. Yang, and Z. Yin, “Broadband terahertz metamaterial absorber based on tantalum nitride,” Appl. Opt. 56(9), 2449–2454 (2017). [CrossRef]  

6. H.-T. Chen, A. J. Taylor, and N. Yu, “A review of metasurfaces: physics and applications,” Rep. Prog. Phys. 79(7), 076401 (2016). [CrossRef]  

7. V. S. Yadav, S. K. Ghosh, S. Bhattacharyya, and S. Das, “Graphene-based metasurface for a tunable broadband terahertz cross-polarization converter over a wide angle of incidence,” Appl. Opt. 57(29), 8720–8726 (2018). [CrossRef]  

8. A. Grigorenko, M. Polini, and K. Novoselov, “Graphene plasmonics,” Nat. Photonics 6(11), 749–758 (2012). [CrossRef]  

9. M. Rahmanzadeh, H. Rajabalipanah, and A. Abdolali, “Multilayer graphene-based metasurfaces: robust design method for extremely broadband, wide-angle, and polarization-insensitive terahertz absorbers,” Appl. Opt. 57(4), 959–968 (2018). [CrossRef]  

10. L. Ye, K. Sui, Y. Zhang, and Q. H. Liu, “Broadband optical waveguide modulators based on strongly coupled hybrid graphene and metal nanoribbons for near-infrared applications,” Nanoscale 11(7), 3229–3239 (2019). [CrossRef]  

11. C. Jördens, M. Scheller, M. Wichmann, M. Mikulics, K. Wiesauer, and M. Koch, “Terahertz birefringence for orientation analysis,” Appl. Opt. 48(11), 2037–2044 (2009). [CrossRef]  

12. R. Amin, Z. Ma, R. Maiti, S. Khan, J. B. Khurgin, H. Dalir, and V. J. Sorger, “Attojoule-efficient graphene optical modulators,” Appl. Opt. 57(18), D130–D140 (2018). [CrossRef]  

13. R.-M. Gao, Z.-C. Xu, C.-F. Ding, and J.-Q. Yao, “Intensity-modulating graphene metamaterial for multiband terahertz absorption,” Appl. Opt. 55(8), 1929–1933 (2016). [CrossRef]  

14. Z. Liu and B. Bai, “Ultra-thin and high-efficiency graphene metasurface for tunable terahertz wave manipulation,” Opt. Express 25(8), 8584–8592 (2017). [CrossRef]  

15. Z. Li, L. Bai, X. Li, E. Gu, L. Niu, and X. Zhang, “U-shaped micro-ring graphene electro-optic modulator,” Opt. Commun. 428, 200–205 (2018). [CrossRef]  

16. K. Wang, M. M. Elahi, L. Wang, K. M. Habib, T. Taniguchi, K. Watanabe, J. Hone, A. W. Ghosh, G.-H. Lee, and P. Kim, “Graphene transistor based on tunable dirac fermion optics,” Proc. Natl. Acad. Sci. 116(14), 6575–6579 (2019). [CrossRef]  

17. L. Ye, K. Sui, and H. Feng, “High-efficiency couplers for graphene surface plasmon polaritons in the mid-infrared region,” Opt. Lett. 45(2), 264–267 (2020). [CrossRef]  

18. A. Suhail, G. Pan, D. Jenkins, and K. Islam, “Improved efficiency of graphene/si schottky junction solar cell based on back contact structure and duv treatment,” Carbon 129, 520–526 (2018). [CrossRef]  

19. J. Yang, Z. Zhu, J. Zhang, C. Guo, W. Xu, K. Liu, X. Yuan, and S. Qin, “Broadband terahertz absorber based on multi-band continuous plasmon resonances in geometrically gradient dielectric-loaded graphene plasmon structure,” Sci. Rep. 8(1), 3239 (2018). [CrossRef]  

20. S. Luo, B. Li, A. Yu, J. Gao, X. Wang, and D. Zuo, “Broadband tunable terahertz polarization converter based on graphene metamaterial,” Opt. Commun. 413, 184–189 (2018). [CrossRef]  

21. X. Zhao, F. Boussaid, A. Bermak, and V. G. Chigrinov, “High-resolution thin “guest-host” micropolarizer arrays for visible imaging polarimetry,” Opt. Express 19(6), 5565–5573 (2011). [CrossRef]  

22. F. Zeng, L. Ye, L. Li, Z. Wang, W. Zhao, and Y. Zhang, “Tunable mid-infrared dual-band and broadband cross-polarization converters based on u-shaped graphene metamaterials,” Opt. Express 27(23), 33826–33839 (2019). [CrossRef]  

23. T. Sun, J. Kim, J. M. Yuk, A. Zettl, F. Wang, and C. Chang-Hasnain, “Surface-normal electro-optic spatial light modulator using graphene integrated on a high-contrast grating resonator,” Opt. Express 24(23), 26035–26043 (2016). [CrossRef]  

24. V. Sharma, D. Madaan, and A. Kapoor, “Extremely short length and high extinction ratio plasmonic te mode pass polarizer,” IEEE Photonics Technol. Lett. 29(7), 559–562 (2017). [CrossRef]  

25. W. Mo, X. Wei, K. Wang, Y. Li, and J. Liu, “Ultrathin flexible terahertz polarization converter based on metasurfaces,” Opt. Express 24(12), 13621–13627 (2016). [CrossRef]  

26. X. Gao, X. Han, W.-P. Cao, H. O. Li, H. F. Ma, and T. J. Cui, “Ultrawideband and high-efficiency linear polarization converter based on double v-shaped metasurface,” IEEE Trans. Antennas Propag. 63(8), 3522–3530 (2015). [CrossRef]  

27. B. Lin, B. Wang, W. Meng, X. Da, W. Li, Y. Fang, and Z. Zhu, “Dual-band high-efficiency polarization converter using an anisotropic metasurface,” J. Appl. Phys. 119(18), 183103 (2016). [CrossRef]  

28. H. Zhang, Y. Liu, Z. Liu, X. Liu, G. Liu, G. Fu, J. Wang, and Y. Shen, “Multi-functional polarization conversion manipulation via graphene-based metasurface reflectors,” Opt. Express 29(1), 70–81 (2021). [CrossRef]  

29. Z. Geng, B. Hähnlein, R. Granzner, M. Auge, A. A. Lebedev, V. Y. Davydov, M. Kittler, J. Pezoldt, and F. Schwierz, “Graphene nanoribbons for electronic devices,” Ann. Phys. 529(11), 1700033 (2017). [CrossRef]  

30. G. W. Hanson, “Dyadic green’s functions and guided surface waves for a surface conductivity model of graphene,” J. Appl. Phys. 103(6), 064302 (2008). [CrossRef]  

31. S. Guan, J. Cheng, T. Chen, and S. Chang, “Bi-functional polarization conversion in hybrid graphene-dielectric metasurfaces,” Opt. Lett. 44(23), 5683–5686 (2019). [CrossRef]  

32. L. Ye, X. Chen, C. Zhu, W. Li, and Y. Zhang, “Switchable broadband terahertz spatial modulators based on patterned graphene and vanadium dioxide,” Opt. Express 28(23), 33948–33958 (2020). [CrossRef]  

33. Y. Zhang, Y. Feng, B. Zhu, J. Zhao, and T. Jiang, “Graphene based tunable metamaterial absorber and polarization modulation in terahertz frequency,” Opt. Express 22(19), 22743–22752 (2014). [CrossRef]  

34. S. Quader, J. Zhang, M. R. Akram, and W. Zhu, “Graphene-based high-efficiency broadband tunable linear-to-circular polarization converter for terahertz waves,” IEEE J. Sel. Top. Quantum Electron. 26(5), 1–8 (2020). [CrossRef]  

35. T. Guo and C. Argyropoulos, “Broadband polarizers based on graphene metasurfaces,” Opt. Lett. 41(23), 5592–5595 (2016). [CrossRef]  

36. Y. Jiang, L. Wang, J. Wang, C. N. Akwuruoha, and W. Cao, “Ultra-wideband high-efficiency reflective linear-to-circular polarization converter based on metasurface at terahertz frequencies,” Opt. Express 25(22), 27616–27623 (2017). [CrossRef]  

37. D. H. Goldstein, Polarized light (CRC, 2017).

38. M. Barkabian, N. Dalvand, H. Zandi, and N. Granpayeh, “Terahertz linear to circular polarization converter based on reflective metasurface,” Sci. Iran. pp. – (2021).

39. Q. Zheng, C. Guo, G. A. Vandenbosch, P. Yuan, and J. Ding, “Ultra-broadband and high-efficiency reflective polarization rotator based on fractal metasurface with multiple plasmon resonances,” Opt. Commun. 449, 73–78 (2019). [CrossRef]  

40. J. Huang, T. Fu, H. Li, Z. Shou, and X. Gao, “A reconfigurable terahertz polarization converter based on metal–graphene hybrid metasurface,” Chin. Opt. Lett. 18(1), 013102 (2020). [CrossRef]  

41. A. Khavasi, “Design of ultra-broadband graphene absorber using circuit theory,” J. Opt. Soc. Am. B 32(9), 1941–1946 (2015). [CrossRef]  

42. A. Khavasi, “Ultra-sharp transmission resonances in periodic arrays of graphene ribbons in te polarization,” J. Lightwave Technol. 34(3), 1020–1024 (2016). [CrossRef]  

Data availability

No data were generated or analyzed in the presented research.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. (a). The schematic view of the proposed three dimensional linear polarization converter (LPC) and linear to circular polarization converter (LTC- PC) and (b) the structure unit cell.
Fig. 2.
Fig. 2. The spectra of (a) Co-polarized and cross-polarized reflection coefficient amplitudes, (b) Qn parameter, and (c) the PCR of the proposed structure of Fig. 1, as a linear polarization converter (LPC).
Fig. 3.
Fig. 3. The spectra of (a) Reflection coefficient amplitude for three different heights of the substrate and (b) the related PCR of the structure of Fig. 1.
Fig. 4.
Fig. 4. The spectra of PCR of structure of Fig. 1 for different angle of incidence.
Fig. 5.
Fig. 5. The LPC spectra of PCR of structure of Fig. 1 for different ${\mu _c}$ of the graphene.
Fig. 6.
Fig. 6. Variation of the cross-polarized reflection coefficient amplitude of structure of Fig. 1 vs. chemical potential, $ {\mu _c}$ , at the central frequency of 0.83 THz.
Fig. 7.
Fig. 7. The spectra of (a) amplitude and (b) phases of the reflection coefficients of the structure of Fig. 1.
Fig. 8.
Fig. 8. The spectra of (a) axial ratio and (b) the ellipticity for proposed LTC-PC structure of Fig. 1.
Fig. 9.
Fig. 9. The equivalent circuit model for the structure of Fig. 1, in (a) $x^{\prime }$ and (b) $y^{\prime }$ directions.
Fig. 10.
Fig. 10. Reflection coefficient amplitudes calculated by equivalent circuit model theory (ECM) and the FDTD method for (a). LPC structure and (b). LTC-PC structure of Fig. 1.

Tables (2)

Tables Icon

Table 1. The optimum dimensions for maximum efficiency of our proposed structure of Fig. 1; and the values of some parameters of Eqs. (1)–(3)

Tables Icon

Table 2. Comparison of our proposed structure and some other similar structures

Equations (23)

Equations on this page are rendered with MathJax. Learn more.

σ = σ i n t r a + σ i n t e r
σ i n t r a = j e 2 k B T π 2 ( ω j 2 Γ ) [ μ c k B T + 2 ln ( e μ c k B T + 1 ) ]
σ i n t e r = j e 2 4 π ln 2 | μ c | ( ω j 2 Γ ) 2 | μ c | + ( ω j 2 Γ )
μ c ν f π ϵ r S i ϵ 0 V g e t s
I n = | R x x | 2 + | R y x | 2
Q n = | R x x | 2 | R y x | 2
U n = 2 | R x x | | R y x | c o s ( Δ ϕ )
V n = 2 | R x x | | R y x | s i n ( Δ ϕ )
A R = 10 l o g ( tan β )
β = 1 2 sin 1 ( V n I n )
P C R = R y x 2 R x x 2 + R y x 2
E inc = E 0 x a x = 2 2 E 0 x ( a x + a y )
R G R = P S 1 2 π 2 e 2 μ c τ
L G R = P S 1 2 π 2 e 2 μ c
C G R = S 1 2 P 2 ϵ e f f q 1
Z G 1 = R G R + j ( ω L G R 1 ω C G R )
L i = η 0 P c π ln [ csc ( π a 2 P ) ] ( 1 α + α 1 ( λ c i / λ ) 2 ) , i = 1 , 2
Z G 2 = j ω L 1 L 2 L 1 + L 2 + P σ a
Z t = Z S i Z A u + j Z S i tan ( β S i h ) Z S i + j Z A u tan ( β S i h )
Z T 1 = Z G 1 Z t Z G 1 + Z t
Z T 2 = Z G 2 Z t Z G 2 + Z t
R y = Z T 1 η 0 Z T 1 + η 0
R x = Z T 2 η 0 Z T 2 + η 0
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.