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

Highly efficient tunable broadband terahertz polarizers based on a graphene metasurface

Open Access Open Access

Abstract

Ultrathin tunable broadband terahertz transmission mode linear-to-circular polarizers are proposed and numerically validated. Each structure consists of a dielectric substrate with slotted graphene patches on one or both sides of the substrate. Both polarizers yield tunable broadband conversion with large transmission coefficients in terahertz frequency range. Compared to the previous polarizers, the proposed structure with a single-layer of graphene patches exhibits a 120% larger fractional bandwidth and 60% larger transmittance. The polarizer with two layers of graphene patches features a 138% larger fractional bandwidth and 25% larger transmittance, compared to the previous polarizers, with excellent axial ratio. The tuning frequency range for the proposed single- and double-layer structures is (2.4 to 4.9) and (1.8 to 5.8) THz, respectively. The tuning is achieved by varying the graphene Fermi energy (chemical potential) from 0.2 to 1 eV, by the application of a DC bias voltage. Moreover, the polarization conversion performance is well maintained under oblique incidence condition, over a wide range of incident angles up to 70°. The size of the square-shaped unit cell is λ0/10, much smaller than the free-space wavelength of the incident waves, λ0, while the thickness is λ0/22. With the above features, the present work offers a further step in developing controllable polarization converters with potential applications in imaging, sensing, and communications.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

To date, several types of electromagnetic polarization converters based on metamaterials and metasurfaces have been reported [111]. Practical implementation of a polarizer can be a challenging task as the frequency band of interest may alter. However, graphene has been used in many tunable devices because of its attractive properties such as flexibility, strong interact with light at THz and infrared frequencies, and its tunable surface conductivity [1217].

According to the above features, recently, various tunable polarizers have been introduced [1844]. Most of these polarizers operate within narrow bandwidths with low transmission efficiencies. For instance, H. Cheng et al. designed tunable infrared polarization converters, based on graphene metasurfaces with cross-shaped unit cells, but the bandwidth was limited [40]. T. Guo et al. designed a broadband polarizer graphene metasurface with a unit cell comprising a rectangular graphene patch [30]; however, the transmission was reduced. Therefore, a tunable broadband polarization converter of high transmission efficiency is an on-going challenging task.

In the present work, two graphene-based tunable broadband polarizers are proposed. Each structure consists of a dielectric substrate with slotted graphene patches, on one or both sides of the substrate. The structure with graphene patches on one side of the substrate is called the single-layer, and the one using two layers of graphene patches on both sides of the substrate is called the double-layer. It is numerically shown the proposed metasurface polarizers yield tunable broadband conversion with large transmission coefficients at terahertz frequencies.

Compared to the previous polarizers, the single-layer structure exhibits a 120% larger fractional bandwidth and 60% larger transmittance, and the double-layer structure features a 138% larger fractional bandwidth and 25% larger transmittance, with an excellent axial ratio (AR). AR is defined as the ratio of the magnitude of the major and minor axes of the polarization ellipse of the transmitted or reflected waves. AR is always a positive real number, often expressed in the logarithmic scale in decibels (dBs). For a linear to circular polarizer, it is desirable to achieve small AR over a broad frequency range of interest. The frequency range over which AR is less than 3 dBs is considered as the bandwidth of the polarizer.

The frequency and polarization states can be dynamically tuned by varying the graphene chemical potential (Fermi energy), using a variable DC bias voltage. The overall operation frequency over which the axial ratio is less than 3 dB, for the proposed single- and double-layer structures are (2.4 to 4.9) THz and (1.8 to 5.8) THz, respectively. In both cases, the wideband operation is well maintained up to an incident angle of 70°.

In Section 2, the proposed converters and their basic features are introduced. In Section 3, a rigorous discussion of the simulation results and comparisons with the previous polarizers are presented. Finally, the conclusions are presented in Section 4.

2. Structure and numerical simulations

Figure 1 shows the schematic of the proposed metasurface made of an array of slotted rectangular graphene patches on a dielectric substrate.

 figure: Fig. 1.

Fig. 1. (a) Schematic of the metasurface made of an array of graphene patches. (b) unit cell.

Download Full Size | PDF

The patches may be placed on one or both sides of the substrate. The surface conductivity of graphene is the sum of the intraband and interband terms. In the low terahertz frequency region, the interband term is ignored and the intraband conductivity is modeled as the Drude expression [45],

$$\sigma = \displaystyle{{-je^2E_F} \over {\pi {\mathrm \hbar }\left( {\omega -j2\Gamma } \right)}}$$
where ${{\boldsymbol E}_{\boldsymbol F}}{\boldsymbol \; \; }$ is the Fermi energy of graphene, ${\boldsymbol \varGamma } = {\boldsymbol ev}_{\boldsymbol F}^2/{\boldsymbol \mu }{{\boldsymbol E}_{\boldsymbol F}}{\boldsymbol \; }$ is the electron scattering rate, ${{\boldsymbol v}_{\boldsymbol F}} \approx {\boldsymbol c}/$300 is the Fermi velocity, ${\boldsymbol \mu }$=1 m2/Vs is the measured DC mobility [30], ${\boldsymbol c}$ the speed of light in free space, and ${\boldsymbol e}$ the electron charge. The relaxation time ${\boldsymbol \tau } = 1/{\boldsymbol \varGamma }$ is, approximately, ${\boldsymbol \tau }$=1 ps. After some optimizations, the structural parameters, with respect to Fig. 1(b), are chosen as: ${\boldsymbol p}$=7.6 µm, ${\boldsymbol y}$=250 nm, L1 = 6.5 µm, L2 = 5.2 µm, ${\boldsymbol t}$=1.5 µm, ${{\boldsymbol \varepsilon }_{\boldsymbol r}}$=2.25, and ${{\boldsymbol E}_{\boldsymbol F}}$=0.95 eV. The parameters are chosen based on simulation, using commercial software package CST, where, the excitation source is considered to be a linearly polarized plane wave of x- or y- polarization. The transmittance, or transmission coefficients amplitude, is defined as:
$${T_{ij}} = \left|{\frac{{{\boldsymbol E}_j^{Trans}}}{{{\boldsymbol E}_i^{Inc}}}} \right|,\,({i.j = x.y} )$$
where ${\boldsymbol E}_{\boldsymbol i}^{{\boldsymbol Inc}}$ is the ${\boldsymbol x}$- or ${\boldsymbol y}$-polarized incident waves, and ${\boldsymbol E}_{\boldsymbol j}^{{\boldsymbol Trans}}$ is the ${\boldsymbol x}$- or ${\boldsymbol y}$-component of the transmitted waves, respectively [30]. The phase terms are defined as
$${\mathrm{\Phi }_{ij}} = arg\left( {\frac{{{\boldsymbol E}_j^{Trans}}}{{{\boldsymbol E}_i^{Inc}}}} \right)$$

3. Results and discussions

3.1 Single-layer structure

For the single-layer structure, the amplitude and phase of the transmission coefficients for the x- and y-polarizations are shown in Fig. 2, for different values of W. The phase difference $\mathrm{\Phi }\Delta = {\mathrm{\Phi }_{\textrm{xx}}} - {\mathrm{\Phi }_{\textrm{yy}}}$ is also shown in Fig. 2(b). The cross-polarization conversion has been negligible. As observed from the Fig. 2, the transmission resonance is shifted when the converter is illuminated by x- or y-polarized waves. The two distinct resonances arise from the asymmetry between the length and the width of the grapheme patches. For W = 0 ($y$=0), the structure simplifies to the one in [30]. In this case, the transmittance, for both polarizations are equal, ${T_{xx}} = {T_{yy}} = 0.55$, at 4.75 THz, while the phase difference is exactly 90°. Interestingly, the 90° phase difference is obtained in a broad frequency range, from 4.3 to 5.1 THz, as shown in Fig. 2(b) [30]. By increasing W, the transmittance can be increases by 60%, from ${T_{xx}} = {T_{yy}} = 0.55,$ for W = 0, to ${T_{xx}} = {T_{yy}} = 0.88,$ for W = 3000 nm. Also, the frequency range in which the phase difference is constant and equal to 90° increases. This behavior can be interpreted through the resonance condition. In general, the resonant frequencies of the periodic graphene patch structure have been demonstrated to be affected by the length of the slot L and the Fermi energy ${E_F}$ following the relationship [17,30],

$${f_r} \propto \; \sqrt {\frac{{{E_F}}}{L}} $$

 figure: Fig. 2.

Fig. 2. (a) Transmission coefficients and (b) phase difference for the x- and y-polarized normally incident waves for the proposed single-layer structure with different values of W, when p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and ${E_F} = 0.95\textrm{}eV$ and comparing them with [30]. The unit of W is nm.

Download Full Size | PDF

According to this relation, by increasing the length of the slot W along the x direction, the resonant frequency associated with ${T_{xx}}$ graph decreases, and as a result, the magnitude and phase of ${T_{xx}}$ shift to the left. Therefore, the intersection of ${T_{xx}}$ with ${T_{yy}}$ curves occurs at a large value, also, the frequency range in which the phase difference is constant increases.

Next, the axial ratio of the transmitted wave is calculated for the proposed single-layer structure. The general formula for calculating AR is [46]:

$$AR = \sqrt {\frac{{E_x^2 + E_y^2 + \sqrt {E_x^4 + E_y^4 + 2E_x^2E_y^2\cos ({2({{\Phi _x} - {\Phi _y}} )} )} }}{{E_x^2 + E_y^2 - \sqrt {E_x^4 + E_y^4 + 2E_x^2E_y^2\cos ({2({{\Phi _x} - {\Phi _y}} )} )} }}} $$
when ${E_x}$ and ${E_y}$ are the amplitudes of the transverse electric field components, and ${\mathrm{\Phi }_x}$ and ${\mathrm{\Phi }_y}$ are the phases of the individual components. In the present work, as in [30], the incident wave is chosen as:
$$\begin{array}{l} \; |{{\boldsymbol E}_i^{Inc}} |= {E_{0i}}\; ;\; \; \; ({i = x.y} )\\ \mathrm{\Phi }_i^{Inc} = {\varphi _{0i}}\; ;\; \; \; ({i = x.y} )\end{array}$$
where ${E_{0i}}$ and ${\varphi _{0i}}$ are given values. In order to obtain the axial ratio of the transmitted wave, the amplitude and phase of the transmitted wave components are obtained as:
$$\begin{array}{l} |{{\boldsymbol E}_i^{Trans}} |= {T_{ii}}\; |{{\boldsymbol E}_i^{Inc}} |\; ;\; \; \; ({i = x.y} )\\ \mathrm{\Phi }_i^{Trans} = {\mathrm{\Phi }_{ii}} + \mathrm{\Phi }_i^{Inc}\; ;\; \; \; ({i = x.y} )\end{array}$$

By placing Eq. (6) in Eq. (7), and Eq. (7) in Eq. (5), AR is obtained as:

$$AR = \sqrt {\frac{{\left[ {T_{xx}^2 + T_{yy}^2 + \sqrt {T_{xx}^4 + T_{yy}^4 + 2T_{xx}^2T_{yy}^2\cos ({2\Delta \Phi } )} } \right]}}{{\left[ {T_{xx}^2 + T_{yy}^2 - \sqrt {T_{xx}^4 + T_{yy}^4 + 2T_{xx}^2T_{yy}^2\cos ({2\; \Delta \Phi } )} } \right]}}} $$
where $\mathrm{\Phi }\Delta = {\mathrm{\Phi }_{\textrm{xx}}} - {\mathrm{\Phi }_{\textrm{yy}}}$.

By inserting the amplitude and phase of the transmission coefficients of Fig. 2 in Eq. (8), the axial ratio of the transmitted wave is obtained, as shown in Fig. 3. It is observed that, we can increase the bandwidth, frequency range over which AR < 3 dB, from 0.8 to 1.5 THz by increasing W, from 0 to 3000 nm. Table 1 shows the center frequency, ${f_0}\; $, transmittance, ${T_{ij}},\; ({i.j = x.y} )$ at ${f_0}\; ,\; $ phase difference, ΔΦ, bandwidth, $\varDelta f$, the fractional bandwidth (FBW), axial ratio, AR, and the averaged axial ratio in the 3-dB bandwidth (AAR), in terms of W. According to the table, increasing W, the center frequency decreases while the transmittance and bandwidth increase. However, increasing W, the axial ratio within the 3-dB bandwidth, AAR, increases slightly, as shown in Fig. 3. This is because the phase difference in this range decreased significantly as W increases, as observed from Fig. 2(b). Increasing W, the frequency range in which the phase difference is near 90° increases, but the phase difference in this range decreases and this leads to an increase of the axial ratio in the 3-dB bandwidth. This indicates that there is a trade-off between the 3-dB bandwidths and the averaged axial ratio in the 3-dB bandwidth, AAR.

 figure: Fig. 3.

Fig. 3. Axial ratio for the proposed single-layer structure with different values of W when p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and ${E_F} = 0.95\textrm{}eV$ and comparing them with [30]. The unit of W is nm.

Download Full Size | PDF

The operation of the polarizer is highly adjustable via the graphene Fermi energy ${E_F}$. The Fermi energy can be changed either by chemical doping or by the application of a bias voltage, ${V_B}$ [30,47]. The relation between ${E_F}\textrm{}$ and ${V_B}\textrm{}$ is given in [30,4852].

The axial ratio for different values of ${E_F}$, 0.4, 0.6, 0.8 and 1 eV, is calculated and the results are shown in Fig. 4. It is observed that, by varying ${E_F}$, one may effectively vary ${f_0}\; $ and the frequency range over which the axial ratio is less than 3 dB. Table 2 shows the variation of the operational frequency range and the corresponding bandwidth, in terms of the graphene Fermi energy. According to the table, by decreasing the Fermi energy, ${E_F}$, the center frequency is decreased, but the axial ratio increased. We consider the frequency range in which the axial ratio is less than 3 dB from 2.4 THz, for ${E_F} = 0.4\textrm{eV}$, to 4.9 THz, for ${E_F} = 1\textrm{eV}$, as the overall operational bandwidth with tuning. The fractional bandwidth, corresponding to the above frequency range, 2.4 to 4.9 THz, is 68%.

Tables Icon

Table 1. Central frequency, ${{\boldsymbol f}_0}$, transmittance at ${{\boldsymbol f}_0}$, phase difference at ${{\boldsymbol f}_0}$, bandwidth, ${\boldsymbol \Delta f}$, fractional bandwidth (FBW), and averaged axial ratio (AAR) over the bandwidths, for different values of ${\boldsymbol W\; }$ of the proposed single-layer structure. Comparison between Case 1 [30] and Case 4 is made in the fifth row. The parameters of the proposed structure are p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${{\boldsymbol \varepsilon }_{\boldsymbol r}}$=2.25, and ${{\boldsymbol E}_{\boldsymbol F}}$= 0.95 eV.

 figure: Fig. 4.

Fig. 4. Axial ratio of the proposed single-layer structure calculated for different values of graphene Fermi energy, ${E_F} = $ 0.4, 0.6, 0.8 and 1 eV.

Download Full Size | PDF

Tables Icon

Table 2. Center frequency ${{\boldsymbol f}_0}$, 3-dB bandwidth $\Delta {f}$, fractional bandwidth FBW and axial ratio AR at different values of Fermi energy, ${{\boldsymbol E}_{\boldsymbol F}} = $1, 0.8, 0.6 and 0.4 eV for the proposed single-layer structure

In order to apply the electrostatic bias voltage to the graphene patches, narrow metallic ribbons (or transparent electrodes) can be placed under the substrate [30,48], as a DC ground plane. One end of the bias voltage ${V_B}\textrm{}$ is connected to the graphene plate and the other end is connected to the DC ground plate. The coplanar graphene patches are also interconnected with narrow metallic ribbons [4951]. The thin wires are short circuited in the DC mode, and open circuited in the AC mode.

3.2 Double-layer structure

In order to increase the phase difference, and reduce the axial ratio, we used two identical graphene patch arrays on both sides of the dielectric substrate. The results are shown in Fig. 5, Fig. 6, and Table 3. The results of [30] are added for comparison. In Table 3, the center frequency ${f_0}$, transmittance at ${f_0}$, phase difference at ${f_0}$, bandwidth, and the average axial ratio over the bandwidth are shown. Comparison with [30] is added in the third row. Comparing the double-layer and single-layer structures, it is seen that by increasing the number of graphene layers, the transmission coefficient decreases by 22% (from 88% to 69%), but $\mathrm{\Phi }\Delta $ and bandwidth increase, and the averaged axial ratio (AAR) over the bandwidth decreases by 48% (from 2.3 to 1.2 dB). According to Fig. 5, Fig. 6, and Table 3, the bandwidth and transmission coefficient of the double-layer structure compared to [30], increase by 138% (from 16.8% to 40%) and 25% (from 55% to 69%), respectively, but the averaged axial ratio (AAR) increased slightly from 1.5 to 1.2 dB.

 figure: Fig. 5.

Fig. 5. (a) Transmission coefficients amplitude, for the x- and y-polarized normally incident waves, Insets: the z-component of the electric field distribution at the two resonances. (b) phase difference for the proposed double-layer structure, for $W = $2500 nm, p = 7.6 µm, y = 250 nm, L1 = 7.2 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and ${E_F} = 0.95\textrm{}eV$ and comparing them with [30].

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Axial ratio for the proposed double-layer structure, with $W$=2500 nm, p = 7.6 µm, y = 250 nm, L1 = 7.2 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, and ${E_F}$= 0.95 eV and comparing it with [30].

Download Full Size | PDF

Tables Icon

Table 3. Center frequency, ${{\boldsymbol f}_0}$, transmittance ${{\boldsymbol T}_{{xx}}} = {{\boldsymbol T}_{{yy}}}$ at ${{\boldsymbol f}_0}$, phase difference, ΔΦ, bandwidth, ${\boldsymbol \Delta f}$, fractional bandwidth (FBW), and averaged axial ratio in the passband (AAR), for the proposed double-layer structure with ${\boldsymbol W}$=2500 nm, ${\boldsymbol p}$=7.6 µm, ${\boldsymbol y}$=250 nm, ${{\boldsymbol L}_1}$=7.2 µm, L2 = 5.2 µm, ${\boldsymbol t}$=1.5 µm, ${{\boldsymbol \varepsilon }_{\boldsymbol r}}$=2.25, and ${{\boldsymbol E}_{\boldsymbol F}}$=0.95 eV and comparing them with [30]

The z-component of the electric field distribution at the two resonance modes are calculated and shown in the inset of Fig. 5(a). Dipolar resonant modes are excited along the edges of the graphene patches with different orientations, depending on the polarization of the incident wave. In this structure, W = 2500 nm and the other parameters are L1 = 7.2 µm, L2 = 5.2 µm, y = 250 nm, t = 1.5 µm, ${\varepsilon _r}$=2.25, and ${E_F}$=0.95 eV.

The transmission amplitude ratio ${{\boldsymbol T}_{xx}}$/${{\boldsymbol T}_{yy}}$ for the double-layer structure is shown in Fig. 7 for ${L_1}$= 6500 and 7200 nm, when $p$=7.6 µm, $y$=250 nm, ${L_2}$=5.2 µm, $t$=1.5 µm, ${\varepsilon _r} = 2.25$, $W$= 2.5 µm, and ${E_F}$= 0.95 eV. It is observed from the figure that increasing L1, the slope of the ${{\boldsymbol T}_{xx}}$/${{\boldsymbol T}_{yy}}$ curve decreases by 70% about the center frequency of f0= 4.75 THz, where, ΔΦ is close to 90°. Therefore, a wideband linear to circular polarization conversion is achieved with a great accuracy.

 figure: Fig. 7.

Fig. 7. Transmission amplitude ratio for the double-layer structure for L1 = 6500 and 7200 nm, with p = 7.6 µm, y = 250 nm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, W = 2.5 µm, and ${E_F}$=0.95 eV.

Download Full Size | PDF

Another feature of the proposed polarizer is its robustness at oblique incidence. The transmittance and phase difference at oblique incident angles from 0° to 70° is shown in Fig. 8. As observed from the figure, the transmittance and phase difference do not vary much in a broad frequency range around 4.75 THz, when the incident angle is varied from 0° to 50°. Figure 9 shows that at the incident angle of 70°, the axial ratio is less than 3 dB at 4.75 THz. This indicates that even at large incidence angles, the conversion is performed well.

 figure: Fig. 8.

Fig. 8. (a) Transmission coefficients amplitude for the x- and y-polarized incident waves and (b) phase difference, under oblique incidence, for angles θ=0 to 70°, by steps of 10°, for the double-layer structure, with p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, W = 2.5 µm, and ${E_F}$=0.95 eV.

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. Axial ratio under oblique incidence, for incident angles of θ=0, to 70°, for the double-layer structure, when p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, W = 2.5 µm and ${E_F}$=0.95 eV.

Download Full Size | PDF

The phase difference of the proposed ultrathin polarization converter can be tuned by controlling the graphene Fermi energy, ${E_F}$, as shown in Fig. 10. It is observed that the proposed structure can be tuned to operate in a broad range of frequency, and its operational bandwidth blueshifts as the Fermi energy increases. In addition, the polarization state of the transmitted wave can be dynamically tuned by varying a DC bias voltage. These features can be interpreted through the relationship (4), according to which, the resonant frequency is directly related to ${E_F}$. Therefore, increasing ${E_F}$, $f_r^{{T_{xx}}}$ and $f_r^{{T_{yy}}}\; $ increase, and as a result, ${T_{xx}}$ and ${T_{yy}}$ transmission diagrams and their phase difference shift towards higher frequencies.

 figure: Fig. 10.

Fig. 10. Calculated phase difference $\Delta \Phi = {\Phi _{xx}} - {\Phi _{yy}}\textrm{}$ of the circular polarization converter as a function of the operation frequency and Fermi energy for the double-layer structure. The parameters are p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and W = 2.5 µm.

Download Full Size | PDF

The frequency range in which the axial ratio is less than 3 dB may be adjusted by changing the Fermi energy, as shown in Fig. 11. For each value of ${E_F}$, the linear to circular polarization conversion occurs at five frequencies. According to Fig. 11, by decreasing the Fermi energy, the central frequency decreases, but the axial ratio increases. Table 4 shows the variation of center frequency, ${f_0}$, the 3-dB bandwidth, $\Delta f$, and the operational frequency range in which the axial ratio is less than 3 dB, for several values of $\textrm{}{E_F}$. According to the table, the operational frequency range blueshifts with ${E_F}$. The 3-dB bandwidth varies from 0.9 THz, for ${E_F}$=0.2 eV, to 1.9 THz, for ${E_F}$=1 eV; the center frequency varies from 2.25 THz, for ${E_F}$=0.2 eV, to 4.85 THz, for ${E_F}$=1 eV. The minimum fractional bandwidth of 39.2% is obtained for ${E_F}$=1 eV. The overall operation bandwidth with tuning is (1.8 to 5.8), for which an overall tuning fractional bandwidth may be calculated to be 105%.

 figure: Fig. 11.

Fig. 11. Axial ratio of the double-layer structure calculated for different values of graphene Fermi energy, ${E_F} = $ 0.2, 0.4, 0.6, 0.8 and 1 eV.

Download Full Size | PDF

Tables Icon

Table 4. Center frequency ${{\boldsymbol f}_0}$, 3-dB bandwidth $\Delta {\boldsymbol f}$, fractional bandwidth FBW and axial ratio AR at different values of Fermi energy, ${{\boldsymbol E}_{\boldsymbol F}} = $1, 0.8, 0.6, 0.4 and 0.2 eV for the double-layer structure

A comparison of the two polarizers of the present work with the previous ones is given in Table 5. The transmission and bandwidth of the present polarizers are well improved compared to the previous ones. Meanwhile, the proposed structures have simple geometries.

Tables Icon

Table 5. Comparison of the proposed structures with previous transmission-mode polarization converters, considering the center frequency ${{\boldsymbol f}_0}$, transmittance at ${{\boldsymbol f}_0}$, fractional bandwidth (FBW), axial ratio (AR), structural complexity, and dependence on the incident angle (θ)

It is worth mentioning that the application of the bias voltage to the graphene in the double-layer structure is similar to that of the single-layer. In both cases, the co-planar graphene patches are interconnected with narrow metallic ribbons (transparent electrodes). The ribbons are short circuited in the DC mode, and open circuited in the AC mode. In the double-layer structure, the electrostatic bias voltage is applied between the two planes.

4. Conclusions

We numerically demonstrated two broadband tunable terahertz linear-to-circular polarization converters, in transmission mode, based on graphene metasurface. Each converter comprised a dielectric substrate with slotted graphene patches on one side (single-layer) or both sides (double-layer) of the substrate. Compared to the previous polarizers, the single-layer structure exhibited 120% larger fractional bandwidth and 60% larger transmittance, and the double-layer converter showed a 138% larger fractional bandwidth and 25% larger transmittance, with excellent axial ratio. Both structures are dynamically tuned via the graphene Fermi energy (chemical potential), by the application of a DC bias voltage. The tunable operational frequency range over which the axial ratio is less than 3 dB, for the single-layer polarizer is (2.4 to 4.9) THz, and for the double-layer structures is (1.8 to 5.8) THz, achieved by varying the graphene Fermi energy from 0.2 to 1 eV. The fractional bandwidth, corresponding to the above frequency ranges, is 68% and 105%, for the single-layer and double-layer structures, respectively. The 3-dB bandwidth of operation blueshifts as the Fermi level is increased, while the fractional bandwidth remains between 37 and 45 percent, for each value of Fermi level. The polarization conversion performance is well maintained under oblique incidence, for a wide range of incident angles from 0° to 70°, indicating robust operation in practice. The unit cell size is ${\lambda _0}/10,\; $ much smaller than the free-space wavelength of the incident waves, $\textrm{}$ while the thickness is ${\lambda _0}/22$.

Funding

Shiraz University.

Disclosures

The authors declared that there are no conflicts of interests.

Data Availability

Simulation data, information, and results are available upon request.

References

1. N. K. Grady, J. E. Heyes, D. R. Chowdhury, Y. Zeng, M. T. Reiten, A. K. Azad, A. J. Taylor, D. A. Dalvit, and H.-T. Chen, “Terahertz metamaterials for linear polarization conversion and anomalous refraction,” Science 340(6138), 1304–1307 (2013). [CrossRef]  

2. L. Cong, N. Xu, J. Gu, R. Singh, J. Han, and W. Zhang, “Highly flexible broadband terahertz metamaterial quarter-wave plate,” Laser Photonics Rev. 8(4), 626–632 (2014). [CrossRef]  

3. S. Bhattacharyya, S. Ghosh, and K. V. Srivastava, “A wideband cross polarization conversion using metasurface,” Radio Sci. 52(11), 1395–1404 (2017). [CrossRef]  

4. Z. Li, R. Zhao, T. Koschny, M. Kafesaki, K. B. Alici, E. Colak, H. Caglayan, E. Ozbay, and C. Soukoulis, “Chiral metamaterials with negative refractive index based on four “U” split ring resonators,” Appl. Phys. Lett. 97(8), 081901 (2010). [CrossRef]  

5. X. Huang, D. Yang, and H. Yang, “Multiple-band reflective polarization converter using U-shaped metamaterial,” J. Appl. Phys. 115(10), 103505 (2014). [CrossRef]  

6. N. Kanda, K. Konishi, and M. Kuwata-Gonokami, “Terahertz wave polarization rotation with double layered metal grating of complimentary chiral patterns,” Opt. Express 15(18), 11117–11125 (2007). [CrossRef]  

7. X.-x. Zheng, Z.-y. Xiao, and X.-y. Ling, “Broadband and efficient reflective polarization converter based on a three-dimensional metamaterial,” Opt. Quantum Electron. 48(10), 461 (2016). [CrossRef]  

8. M. I. Shalaev, J. Sun, A. Tsukernik, A. Pandey, K. Nikolskiy, and N. M. Litchinitser, “High-efficiency all-dielectric metasurfaces for ultracompact beam manipulation in transmission mode,” Nano Lett. 15(9), 6261–6266 (2015). [CrossRef]  

9. S. Zhang, J. Zhou, Y.-S. Park, J. Rho, R. Singh, S. Nam, A. K. Azad, H.-T. Chen, X. Yin, and A. J. Taylor, “Photoinduced handedness switching in terahertz chiral metamolecules,” Nat. Commun. 3(1), 942 (2012). [CrossRef]  

10. Y. Yu, F. Xiao, I. D. Rukhlenko, and W. Zhu, “High-efficiency ultra-thin polarization converter based on planar anisotropic transmissive metasurface,” AEU - Int. J. Electron. Commun. 118, 153141 (2020). [CrossRef]  

11. S. Sarkar and B. Gupta, “Multiband ultrathin chiral metasurface for polarization conversion, symmetric and asymmetric transmission,” AEU - Int. J. Electron. Commun. 142, 154009 (2021). [CrossRef]  

12. A. Farmani, M. Miri, and M. H. Sheikhi, “Analytical modeling of highly tunable giant lateral shift in total reflection of light beams from a graphene containing structure,” Opt. Commun. 391, 68–76 (2017). [CrossRef]  

13. D. L. Sounas and C. Caloz, “Electromagnetic nonreciprocity and gyrotropy of graphene,” Appl. Phys. Lett. 98(2), 021911 (2011). [CrossRef]  

14. D. L. Sounas and C. Caloz, “Gyrotropy and nonreciprocity of graphene for microwave applications,” IEEE Trans. Microwave Theory Tech. 60(4), 901–914 (2012). [CrossRef]  

15. J. Ding, B. Arigong, H. Ren, J. Shao, M. Zhou, Y. Lin, and H. Zhang, “Mid-infrared tunable dual-frequency cross polarization converters using graphene-based L-shaped nanoslot array,” Plasmonics 10(2), 351–356 (2015). [CrossRef]  

16. Y. Zhao and A. Alù, “Manipulating light polarization with ultrathin plasmonic metasurfaces,” Phys. Rev. B 84(20), 205428 (2011). [CrossRef]  

17. Z. Hamzavi-Zarghani, A. Yahaghi, and L. Matekovits, “Electrically tunable mantle cloaking utilizing graphene metasurface for oblique incidence,” AEU - Int. J. Electron. Commun. 116, 153080 (2020). [CrossRef]  

18. Z. Hamzavi-Zarghani, A. Yahaghi, L. Matekovits, and I. Peter, “Tunable polarization converter based on graphene metasurfaces,” in 2018 IEEE Radio and Antenna Days of the Indian Ocean (RADIO), 2018: IEEE, pp. 1–2.

19. J. Zhu, S. Li, L. Deng, C. Zhang, Y. Yang, and H. Zhu, “Broadband tunable terahertz polarization converter based on a sinusoidally-slotted graphene metamaterial,” Opt. Mater. Express 8(5), 1164–1173 (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. M. Chen, W. Sun, J. Cai, L. Chang, and X. Xiao, “Frequency-tunable mid-infrared cross polarization converters based on graphene metasurface,” Plasmonics 12(3), 699–705 (2017). [CrossRef]  

22. X. Yu, X. Gao, W. Qiao, L. Wen, and W. Yang, “Broadband tunable polarization converter realized by graphene-based metamaterial,” IEEE Photonics Technol. Lett. 28(21), 2399–2402 (2016). [CrossRef]  

23. M. Chen, X. Xiao, L. Chang, C. Wang, and D. Zhao, “High-efficiency and multi-frequency polarization converters based on graphene metasurface with twisting double L-shaped unit structure array,” Opt. Commun. 394, 50–55 (2017). [CrossRef]  

24. 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]  

25. L. Zeng, T. Huang, G.-B. Liu, and H.-F. Zhang, “A tunable ultra-broadband linear-to-circular polarization converter containing the graphene,” Opt. Commun. 436, 7–13 (2019). [CrossRef]  

26. J. Li, P. Yu, H. Cheng, W. Liu, Z. Li, B. Xie, S. Chen, and J. Tian, “Optical polarization encoding using graphene-loaded plasmonic metasurfaces,” Adv. Opt. Mater. 4(1), 91–98 (2016). [CrossRef]  

27. Y. Zhang, Y. Feng, B. Zhu, J. Zhao, and T. Jiang, “Switchable quarter-wave plate with graphene based metamaterial for broadband terahertz wave manipulation,” Opt. Express 23(21), 27230–27239 (2015). [CrossRef]  

28. J. Peng, Z. Zhu, J. Zhang, X. Yuan, and S. Qin, “Tunable terahertz half-wave plate based on hybridization effect in coupled graphene nanodisks,” Appl. Phys. Express 9(5), 055102 (2016). [CrossRef]  

29. C. Yang, Y. Luo, J. Guo, Y. Pu, D. He, Y. Jiang, J. Xu, and Z. Liu, “Wideband tunable mid-infrared cross polarization converter using rectangle-shape perforated graphene,” Opt. Express 24(15), 16913–16922 (2016). [CrossRef]  

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

31. H. Cheng, S. Chen, P. Yu, J. Li, B. Xie, Z. Li, and J. Tian, “Dynamically tunable broadband mid-infrared cross polarization converter based on graphene metamaterial,” Appl. Phys. Lett. 103(22), 223102 (2013). [CrossRef]  

32. A. Forouzmand and A. B. Yakovlev, “Electromagnetic cloaking of a finite conducting wedge with a nanostructured graphene metasurface,” IEEE Trans. Antennas Propag. 63(5), 2191–2202 (2015). [CrossRef]  

33. H. Cheng, S. Chen, P. Yu, W. Liu, Z. Li, J. Li, B. Xie, and J. Tian, “Refraction: Dynamically Tunable Broadband Infrared Anomalous Refraction Based on Graphene Metasurfaces (Advanced Optical Materials 12/2015),” Adv. Opt. Mater. 3(12), 1743 (2015). [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. M. Sajjad, X. Kong, S. Liu, A. Ahmed, S. U. Rahman, and Q. Wang, “Graphene-based THz tunable ultra-wideband polarization converter,” Phys. Lett. A 384(23), 126567 (2020). [CrossRef]  

36. R. Zhang, B. You, S. Wang, K. Han, X. Shen, and W. Wang, “Broadband and switchable terahertz polarization converter based on graphene metasurfaces,” Opt. Express 29(16), 24804–24815 (2021). [CrossRef]  

37. 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]  

38. M. Masyukov, A. Vozianova, A. Grebenchukov, K. Gubaidullina, A. Zaitsev, and M. Khodzitsky, “Optically tunable terahertz chiral metasurface based on multi-layered graphene,” Sci. Rep. 10(1), 3157 (2020). [CrossRef]  

39. B. Bakhtiari and H. Oraizi, “Tunable Terahertz Polarization Converter Based on Graphene Metasurfaces,” in 2020 14th European Conference on Antennas and Propagation (EuCAP), 2020: IEEE, pp. 1–4.

40. H. Cheng, S. Chen, P. Yu, J. Li, L. Deng, and J. Tian, “Mid-infrared tunable optical polarization converter composed of asymmetric graphene nanocrosses,” Opt. Lett. 38(9), 1567–1569 (2013). [CrossRef]  

41. H.-F. Zhang, L. Zeng, G.-B. Liu, and T. Huang, “Tunable linear-to-circular polarization converter using the graphene transmissive metasurface,” IEEE Access 7, 158634–158642 (2019). [CrossRef]  

42. X. Zhang, H. Ye, Y. Zhao, and H. Zhang, “Linear-to-Circular Polarization Converter with Adjustable Bandwidth Realized by the Graphene Transmissive Metasurface,” Plasmonics 17(3), 1079–1089 (2022). [CrossRef]  

43. M. Barkabian, N. Dalvand, H. Zandi, and N. Granpayeh, “Terahertz linear to circular polarization converter based on reflective metasurface,” Scientia Iranica (2021).

44. X. Qi, J. Zou, C. Li, J. Zhang, C. Guo, and Z. Zhu, “Graphene-based electrically controlled terahertz polarization switching between a quarter-wave plate and half-wave plate,” Opt. Express 28(26), 39430–39442 (2020). [CrossRef]  

45. P.-Y. Chen, C. Argyropoulos, and A. Alu, “Terahertz antenna phase shifters using integrally-gated graphene transmission-lines,” IEEE Trans. Antennas Propag. 61(4), 1528–1537 (2013). [CrossRef]  

46. C. Balanis, “Antenna theory: analysis and design (4-th ed.). John Wiley & Sons, Inc,” 2016.

47. B. Wu, H. M. Tuncer, M. Naeem, B. Yang, M. T. Cole, W. I. Milne, and Y. Hao, “Experimental demonstration of a transparent graphene millimetre wave absorber with 28% fractional bandwidth at 140 GHz,” Sci. Rep. 4(1), 4130 (2015). [CrossRef]  

48. A. Dolatabady and N. Granpayeh, “Manipulation of the Faraday rotation by graphene metasurfaces,” J. Magn. Magn. Mater. 469, 231–235 (2019). [CrossRef]  

49. M. Huang, Y. Cheng, Z. Cheng, H. Chen, X. Mao, and R. Gong, “Based on graphene tunable dual-band terahertz metamaterial absorber with wide-angle,” Opt. Commun. 415, 194–201 (2018). [CrossRef]  

50. A. Ottomaniello, “Giant Faraday rotation in a hybrid graphene-split ring resonators metasurface with magneto-electric tunability,” 2016.

51. S. Barzegar-Parizi and A. Ebrahimi, “Ultrathin, polarization-insensitive multi-band absorbers based on graphene metasurface with THz sensing application,” J. Opt. Soc. Am. B 37(8), 2372–2381 (2020). [CrossRef]  

52. I. Mazraeh-fard and A. Alighanbari, “Equivalent circuit model for analysis and design of graphene based tunable,” Appl. Opt. 61(19), 5760–5768 (2022). [CrossRef]  

53. S. K. Ghosh, S. Bhattacharyya, and S. Das, “Graphene-based metasurface for wideband linear to circular polarization conversion,” in 2020 URSI Regional Conference on Radio Science (URSI-RCRS), 2020: IEEE, pp. 1–4.

Data Availability

Simulation data, information, and results are available upon request.

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 (11)

Fig. 1.
Fig. 1. (a) Schematic of the metasurface made of an array of graphene patches. (b) unit cell.
Fig. 2.
Fig. 2. (a) Transmission coefficients and (b) phase difference for the x- and y-polarized normally incident waves for the proposed single-layer structure with different values of W, when p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and ${E_F} = 0.95\textrm{}eV$ and comparing them with [30]. The unit of W is nm.
Fig. 3.
Fig. 3. Axial ratio for the proposed single-layer structure with different values of W when p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and ${E_F} = 0.95\textrm{}eV$ and comparing them with [30]. The unit of W is nm.
Fig. 4.
Fig. 4. Axial ratio of the proposed single-layer structure calculated for different values of graphene Fermi energy, ${E_F} = $ 0.4, 0.6, 0.8 and 1 eV.
Fig. 5.
Fig. 5. (a) Transmission coefficients amplitude, for the x- and y-polarized normally incident waves, Insets: the z-component of the electric field distribution at the two resonances. (b) phase difference for the proposed double-layer structure, for $W = $2500 nm, p = 7.6 µm, y = 250 nm, L1 = 7.2 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and ${E_F} = 0.95\textrm{}eV$ and comparing them with [30].
Fig. 6.
Fig. 6. Axial ratio for the proposed double-layer structure, with $W$=2500 nm, p = 7.6 µm, y = 250 nm, L1 = 7.2 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, and ${E_F}$= 0.95 eV and comparing it with [30].
Fig. 7.
Fig. 7. Transmission amplitude ratio for the double-layer structure for L1 = 6500 and 7200 nm, with p = 7.6 µm, y = 250 nm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, W = 2.5 µm, and ${E_F}$=0.95 eV.
Fig. 8.
Fig. 8. (a) Transmission coefficients amplitude for the x- and y-polarized incident waves and (b) phase difference, under oblique incidence, for angles θ=0 to 70°, by steps of 10°, for the double-layer structure, with p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, W = 2.5 µm, and ${E_F}$=0.95 eV.
Fig. 9.
Fig. 9. Axial ratio under oblique incidence, for incident angles of θ=0, to 70°, for the double-layer structure, when p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r}$=2.25, W = 2.5 µm and ${E_F}$=0.95 eV.
Fig. 10.
Fig. 10. Calculated phase difference $\Delta \Phi = {\Phi _{xx}} - {\Phi _{yy}}\textrm{}$ of the circular polarization converter as a function of the operation frequency and Fermi energy for the double-layer structure. The parameters are p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ${\varepsilon _r} = 2.25$, and W = 2.5 µm.
Fig. 11.
Fig. 11. Axial ratio of the double-layer structure calculated for different values of graphene Fermi energy, ${E_F} = $ 0.2, 0.4, 0.6, 0.8 and 1 eV.

Tables (5)

Tables Icon

Table 1. Central frequency, f 0 , transmittance at f 0 , phase difference at f 0 , bandwidth, Δ f , fractional bandwidth (FBW), and averaged axial ratio (AAR) over the bandwidths, for different values of W of the proposed single-layer structure. Comparison between Case 1 [30] and Case 4 is made in the fifth row. The parameters of the proposed structure are p = 7.6 µm, y = 250 nm, L1 = 6.5 µm, L2 = 5.2 µm, t = 1.5 µm, ε r =2.25, and E F = 0.95 eV.

Tables Icon

Table 2. Center frequency f 0 , 3-dB bandwidth Δ f , fractional bandwidth FBW and axial ratio AR at different values of Fermi energy, E F = 1, 0.8, 0.6 and 0.4 eV for the proposed single-layer structure

Tables Icon

Table 3. Center frequency, f 0 , transmittance T x x = T y y at f 0 , phase difference, ΔΦ, bandwidth, Δ f , fractional bandwidth (FBW), and averaged axial ratio in the passband (AAR), for the proposed double-layer structure with W =2500 nm, p =7.6 µm, y =250 nm, L 1 =7.2 µm, L2 = 5.2 µm, t =1.5 µm, ε r =2.25, and E F =0.95 eV and comparing them with [30]

Tables Icon

Table 4. Center frequency f 0 , 3-dB bandwidth Δ f , fractional bandwidth FBW and axial ratio AR at different values of Fermi energy, E F = 1, 0.8, 0.6, 0.4 and 0.2 eV for the double-layer structure

Tables Icon

Table 5. Comparison of the proposed structures with previous transmission-mode polarization converters, considering the center frequency f 0 , transmittance at f 0 , fractional bandwidth (FBW), axial ratio (AR), structural complexity, and dependence on the incident angle (θ)

Equations (8)

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

σ = j e 2 E F π ( ω j 2 Γ )
T i j = | E j T r a n s E i I n c | , ( i . j = x . y )
Φ i j = a r g ( E j T r a n s E i I n c )
f r E F L
A R = E x 2 + E y 2 + E x 4 + E y 4 + 2 E x 2 E y 2 cos ( 2 ( Φ x Φ y ) ) E x 2 + E y 2 E x 4 + E y 4 + 2 E x 2 E y 2 cos ( 2 ( Φ x Φ y ) )
| E i I n c | = E 0 i ; ( i = x . y ) Φ i I n c = φ 0 i ; ( i = x . y )
| E i T r a n s | = T i i | E i I n c | ; ( i = x . y ) Φ i T r a n s = Φ i i + Φ i I n c ; ( i = x . y )
A R = [ T x x 2 + T y y 2 + T x x 4 + T y y 4 + 2 T x x 2 T y y 2 cos ( 2 Δ Φ ) ] [ T x x 2 + T y y 2 T x x 4 + T y y 4 + 2 T x x 2 T y y 2 cos ( 2 Δ Φ ) ]
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.