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

Broadband of linear-to-linear and double-band of linear-to-circular polarization converter based on a graphene sheet with a π-shaped hollow array

Open Access Open Access

Abstract

In this manuscript, we propose a tunable multifunctional reflection polarization converter based on a graphene sheet, which is decorated by a π-shaped carved-hollow array. Our design can achieve the linear-to-linear (LTL) polarization conversion with 3.58 THz bandwidth and two linear-to-circular (LTC) polarization conversion bands with opposite handedness, which originates from the excitation of the three graphene surface plasmons (GSPs) modes. The LTL polarization converter can be switched to the LTC polarization converter by only changing the medium height. The performances of this multifunctional reflection polarization converter under different geometric parameters, incident angles, and electron scattering time of graphene are investigated. Besides, it is found that the broadband of LTL and double bands of LTC have the blue shift with increasing the Fermi energy (Ef) of graphene. Especially, the conversion between line-to-left-circular and line-to-right-circular polarization converter can be achieved by changing Ef. All the simulation results are conducted by finite element method, and our design has potential applications in integrated terahertz systems.

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

1. Introduction

Terahertz wave manipulation is a technology with great potential, which is of great significance in the research and application of wireless communication, space exploration and other fields [14]. In recent years, manipulating the polarization state of terahertz waves has become a hot research in high resolution imaging, terahertz communication, remote sensing and other mid-infrared applications [57]. The traditional polarization converter is usually based on the birefringence effect generated by grating [8] and dichroic crystal [9] based on natural birefringent materials. The realization of polarization manipulation requires phase delay between two orthogonal components [10,11]. However, it is difficult for traditional devices to achieve nanometer size because of their large structure and long distance to achieve phase accumulation [12,13]. In this case, the metasurface provides a possible way for a highly integrated polarizer. As a kind of artificial subwavelength structure, metasurface makes up for the shortage of the traditional photoelectric element and becomes a hot research [1418]. Surface plasmon can manipulate the amplitude, phase and polarization [1921] of the electromagnetic wave in nanoscale and enhance the interaction between light and matter. This may provide scientific basis for the highly integrated miniaturization of nanoscale optoelectronic devices and the application of metamaterial design [22,23]. Due to its ultra-thin shape and efficient operation, several metasurface structures have been proposed to control the polarization of the terahertz electromagnetic waves such as linear-to-circular (LTC) [2428], linear-to-linear (LTL) [2931] and circular-to-circular [32,33] polarization conversion. Metasurface devices have developed a unique method to manipulate the polarization state of terahertz waves. However, most of these polarized operating devices are not dynamically adjustable.

Graphene is a single two-dimensional plane of carbon atoms packed tightly together. Novoselov et al. prepared graphene by mechanical stripping in 2004 [34] and found that graphene has excellent photoelectric properties. Graphene supports low loss surface plasmons in the mid-infrared region [3537], which provides a possible way to manipulate the amplitude and phase of electromagnetic wave. The thickness of the graphene is ultra-thin and the electrical properties of graphene can be adjusted by doping concentration and bias voltage [3842]. So far, the graphene metasurface has been extensively studied for tunable polarization control [38,44,45]. But, the multifunctional polarization converter, only based on one simple design of graphene metasurface, has been rarely reported.

Multifunctional polarization converters can achieve the switch between the LTL and LTC polarization conversion by changing the parameters of the metasurface [43,46,47]. Yao et al. proposed a multifunctional polarization converter based on the cross double-ellipse graphene patches, but the bandwidths of LTL and LTC polarization converters are both smaller than 1THz [43]. Guo et al. proposed the L-shaped and rectangular graphene patches to achieve LTL and LTC polarization conversion, respectively [46]. The LTC and LTL polarization converters were realized by the single-layer and double-layer graphene metasurfaces, respectively [47]. Our design can achieve broadband LTL polarization conversion and two LTC polarization conversion bands with opposite handedness by only changing the medium height.

In this paper, a tunable multifunctional reflective polarization converter is proposed, which is composed of a graphene sheet decorated by a π-shaped carved-hollow array. The LTL polarization conversion with 3.58 THz bandwidth can be achieved, and the polarization conversion ratio is larger than 80%. By carefully choosing the medium height, the broad band of LTL is changed to the two bands of LTC polarization conversion, including linear to left-handed and right-handed circular polarization conversion simultaneously. The multifunctional reflection polarization converter originates from the excitation of three graphene surface plasmons (GSPs) modes. The influences of incident angle, electron scattering time of graphene and geometric parameter over LTL and LTC polarization converter are investigated. In addition, the broad band of LTL and double bands of LTC all have blue shifts with increasing the Fermi energy (Ef) of graphene, and the band of line-to-left-circular-polarization (LTLCP) can be dynamically tuned to the one of line-to-right-circular-polarization (LTRCP), by changing Ef. All the simulation results based on finite element method verify the design of this paper. The broadband effect for LTL polarization converter has been reported a lot [2931], but our proposed graphene metasurface can realize the switch from broadband LTL to LTC polarization converter just by changing the medium height and without changing the shape of metasurface, which will improve the practicability and adjustability of this design. So the proposed polarization converter has promising potential on the miniaturization and integration of the terahertz system.

2. Design and theories

The structural diagram of proposed polarization converter is shown in Fig. 1(a). It consists of a graphene sheet decorated with the periodical π-shaped hollows. The graphene sheet is placed over an intermediate silica layer, which is supported by Au substrate to improve the reflection. One unit of the proposed structure is shown in Fig. 1(b). The lengths L1 and L2 of the lateral and vertical arms are 505 and 665 nm, respectively. The width of arms and the distance between the two vertical arms are d(=60 nm) and g(=110 nm), respectively. In our simulation, the periods (Px and Py) along the x and y directions are both set as 990 nm. Our proposed graphene metasurface could be fabricated by the existing preparation methods. For example, Helium Ion Beam direct writing technique could be employed to pattern graphene in nanometer-scale graphene metasurface under 5nm extreme precision [48,49], which meets our demand. The incident angle is defined as θ as shown in Fig. 1(c).

 figure: Fig. 1.

Fig. 1. (a) Schematic diagram of multifunctional polarization converter, in which yellow and gray regions represent gold and silicon dioxide layer, respectively. (b) One unit of (a). Geometric parameters include: L1 = 505 nm, L2=665 nm, d=60 nm, g=110 nm and Px= Py = 990 nm. Besides, the incident angle is defined as θ as shown in Fig. 1(c).

Download Full Size | PDF

In this paper, all the simulation results are obtained by finite element method (FEM) via frequency domain solver in COMSOL Multiphysics. The Floquet boundary condition is employed along the x and y directions. The graphene thickness (Δ) is chosen as 1 nm, and the graphene permittivity ε is as follows:

$$\varepsilon = 1\textrm{ + }\frac{{i{\sigma _g}{\eta _0}}}{{{k_0}\Delta }}$$
where η0 (≈ 377 Ω) is the impedance of air and k0 is the wave vector of the incident light in air. The optical conductivity of graphene σg can be calculated as:
$$\begin{aligned} \mathop \sigma \nolimits_g &= \frac{{\mathop {ie}\nolimits^2 \mathop E\nolimits_f }}{{\pi {\hbar ^2}({\mathop {\omega + i\tau }\nolimits^{ - 1} } )}} + \frac{{\mathop {ie}\nolimits^2 }}{{4\pi \hbar }}\ln \left[ {\frac{{\mathop {2E}\nolimits_f - (\omega + i{\tau^{ - 1}})\hbar }}{{\mathop {2E}\nolimits_f + (\omega + i{\tau^{ - 1}})\hbar }}} \right] \\ &+ \frac{{\mathop {ie}\nolimits^2 {k_B}T}}{{\pi {\hbar ^2}({\omega + i{\tau^{ - 1}}} )}}\ln \left[ {\exp \left( { - \frac{{{E_f}}}{{{k_B}T}}} \right) + 1} \right] \end{aligned}$$

Here, Ef is Fermi energy of graphene, T(=300K) is temperature, ω is frequency of the incident light, e is electronic charge, τ is electron scattering time, ℏ is reduced Planck’s constant, and kB is Boltzmann constant. For dynamically controlling the electromagnetic property of graphene, a thin transparent ion-gel layer can be covered at top of the structure so that each cell could be electrically biased through transparent electrodes. If large Fermi energy is required, the chemical doping is employed before electrical stimuli, and the targeted Fermi energy is obtained. Assume that the dielectric constant of silica is 2.1. The Drude model is used to simulate the metal layer. The typical carrier concentration was 4.11×1027 m−3 and the scattering time was 20 fs.

The incident and reflected electric fields can be calculated by the following formula:

$$\left( \begin{array}{l} E_A^r\\ E_B^r \end{array} \right) = r\left( \begin{array}{l} E_A^i\\ E_B^i \end{array} \right)$$

where EAi and EBi represent the electric field of incident light in the +45° and −45° to positive x-direction respectively. EAr and EBr represent the electric field of the reflected light in the +45° and −45° directions, respectively. The general reflection matrix of the reflective graphene metasurface can be expressed as:

$$r = \left( \begin{array}{ll} {r_{AA}} &{r_{AB}}\\ {r_{BA}} &{r_{BB}} \end{array} \right)$$
where ${r_{jk}} = E_j^r/E_k^i$ (j, k = A, B)) represents the reflection coefficient. In this manuscript, the polarization angle of incident beam is 45° respect to the x axis, so rAA and rBA are the co-polarization and cross-polarization reflection coefficients, and ${\varphi _{AA}}$ and ${\varphi _{BA}}$ are the corresponding phases, respectively. Hence, the polarization conversion ratio (PCR) can be described as [19]:
$$PCR = \frac{{{R_{BA}}}}{{{R_{AA}} + {R_{BA}}}}$$

Here, we define co-polarization reflection RAA and cross-polarization reflection RBA as:

$${R_{AA}} = {|{{r_{AA}}} |^2},{R_{BA}} = {|{{r_{BA}}} |^2}$$

Here, the proposed design is believed to achieve cross polarization conversion when PCR is over 80%.

To verify the polarization state of the reflected wave, the Stokes parameters [50] are introduced as follows:

$$\begin{aligned} {S_0} &= {|{{r_{AA}}} |^2} + {|{{r_{BA}}} |^2}\\ {S_1} &= {|{{r_{AA}}} |^2} - {|{{r_{BA}}} |^2}\\ {S_2} &= 2|{{r_{AA}}} ||{{r_{BA}}} |\cos ({\Delta \varphi } )\\ {S_3} &= 2|{{r_{AA}}} ||{{r_{BA}}} |\sin ({\Delta \varphi } )\end{aligned}$$
where S0 presents the intensity of output reflected waves and S3 demonstrates the circular polarized component. When $|{{{r}_{{AA}}}} | = |{{{r}_{{BA}}}} |$, and ${\Delta} \; \varphi \; ( = {\varphi_{BA}} - {\varphi_{AA}})= {\pi}/2 + {k}{\pi}$ (k is an integer), the LTC polarization converter can be realized. Here, we define the ellipticity χ as χ = S3/S0 to indicate the performance of the LTC polarization converter. When χ is equal to 1, the reflected wave is the perfect left-handed circularly polarized wave, and when χ is equal to −1, the reflected wave is the perfect right-handed circularly polarized wave.

3. Broadband of LTL polarization converter

In this section, the graphene sheet is placed over an intermediate silica layer (the thickness of silica layer m = 3600 nm). A tunable broadband LTL polarization converter in mid-infrared region can be achieved when Fermi energy (Ef) of graphene is considered to be 1 eV. The incident beam with polarization angle +45° (the angle between polarization angle and positive x-direction) is converted to reflected beam with polarization angle −45° under normal incidence. The results of PCR, co-polarization (RAA) and cross-polarization (RBA) reflections are shown in Fig. 2. There are three dips in the co-polarization (RAA) reflection spectrum, which are located at 14.02, 15.78 and 17.36 THz respectively. And there are three cross-polarization (RBA) reflection peaks around the same frequencies, respectively. The PCRs reach the three peaks at 14.02, 15.78 and 17.36 THz, whose values are 91.16%, 99.78% and 97.69%, respectively. In the PCR spectrum from 13.94 to 17.52 THz, PCRs are over 80%, proving that the structure can realize LTL polarization conversion with the bandwidth of 3.58 THz, yielding a fractional bandwidth (the ratio of bandwidth to central frequency) about 23% with reference to the central frequency of 15.73 THz. Thus, a broadband LTL polarization converter is realized in reflective mode.

 figure: Fig. 2.

Fig. 2. The simulation results of (a) co-polarization and cross-polarization reflection (RAA and RBA), and (b) PCR are plotted in black, blue and black solid line, respectively.

Download Full Size | PDF

We attempt to further explain the physical mechanism of broadband LTL polarization converter, and find that there are three eigenmodes related to the physical origin. The incident polarization light (polarization angle +45°) can be divided into two orthogonal polarization light components Ex and Ey, whose polarization direction is along the x-axis and y-axis respectively. Rxx is the co-polarization reflection under the incident Ex and Ryy is the co-polarization reflection under the incident Ey. In Fig. 3(a), the two reflectance spectra (drawn with blue and green solid lines, respectively) show three dips at 13.84, 16.00 and 17.58 THz which are very close to three PCR peaks (14.02, 15.78 and 17.36 THz). The three dips originate from three eigenmodes, which are excited by these two components (Ex and Ey). Ex excited one mode at 16.00 THz and Ey excited two modes at 13.84 and 17.58 THz. As shown in Figs. 3(b), (c) and (d), the three modes show symmetric (16.00THz) and antisymmetric (13.84 and 17.58 THz) magnetic field distributions about the y-axis, which originate from the GSPs resonances, respectively. In addition, the phase delay (PD) between Rxx and Ryy in the reflected wave occurs in the three resonant modes, which is approximately 180° at 14.02, 15.78 and 17.36 THz respectively, shown in Fig. 3(a). Thus, the superposition of two reflected orthogonal components produces a polarization rotation of 90° at 14.02, 15.78 and 17.36 THz. Thus, the three PCR peaks happen at these three frequencies as shown in Fig. 2(b), and the broadband LTL converter with 3.58 THz bandwidth is achieved.

 figure: Fig. 3.

Fig. 3. (a) The co-polarized reflection spectrum under x-polarized Ex and y-polarized Ey incidence, respectively. The simulation results of Rxx and Ryy are drawn with blue and green solid lines. The dotted red line represents the PD between Rxx and Ryy. Magnetic field distributions (Hz) of 13.84 THz(b) and 17.58 THz (d) under Ey incidence, and 16.00 THz (c) under Ex incidence, respectively.

Download Full Size | PDF

Next, we study the broadband performance of LTL polarization converter with different geometric parameters L1, L2 and d. When L1 increases from 465 to 545 nm, other geometric parameters remain unchanged, and PCRs are shown in Fig. 4(a). The third PCR peak corresponding to antisymmetric GSPs mode has a significant red shift, the first and second PCR peaks are almost unchanged, and the bandwidth of LTL polarization converter narrows with increasing L1. Besides, the LTL polarization converter shows the bad broadband performance (PCR<80%) when L1 is smaller than 485 nm. In Fig. 4(b), when L2 increases from 625 to 705 nm, the frequency of first PCR peak has a red shift from 14.74 to 13.36 THz, and other two PCR peaks are almost unchanged. Thus, the bandwidth of LTL polarization converter broadens with increasing L2. In Fig. 4(c), the positions of first and third PCR peak are almost unchanged, but the second PCR peak corresponding to symmetric GSPs mode has a blue shift from 15.10 to 16.86 THz with the increase of d. The bandwidth of LTL polarization converter broadens with increasing d, and LTL polarization converter shows bad broadband performance when d is larger than 60 nm.

 figure: Fig. 4.

Fig. 4. Simulated PCR of the LTL polarization converter under normal incidence with different (a) length L1, (b) length L2 and (c) width d of the π-shaped hollow array. Other parameters are the same with those of Fig. 2.

Download Full Size | PDF

Figure 5 shows the broadband performance of LTL polarization converter under different incident angles from 0° to 70°, and other parameters are the same with those of Fig. 2. As the incident angle increases from 0° to 40°, the operating frequency range of the LTL polarization converter (PCR is over 80%) remains from 13.94 to 17.52 THz and the fractional bandwidth is unchanged. The resonant frequencies of three GSPs modes are insensitive to the incident angle. The simulation results clearly show that there is no significant degradation of PCR performance in the range from 0° to 40°. When the incident angle is larger than 40°, the broadband performance turns down, because the interaction between graphene and electromagnetism weakens the intensity of the reflected lights due to the large incident angle.

 figure: Fig. 5.

Fig. 5. Simulation of the PCR under different incident angles from 0° to 70°. Other parameters are the same with those of Fig. 2.

Download Full Size | PDF

We discuss the tunable performance of the LTL polarization converter by changing Fermi energy (Ef). The Fermi energy (Ef) of graphene can be dynamically regulated by changing chemically doped or electronically gated [39], showing the great potential of tunable devices. As the Fermi energy (Ef) decreases from 1 to 0.84 eV, the ideal wideband LTL effect with 3.58 THz bandwidth is almost maintained. The resonant frequency positions have red shifts from 14.02, 15.78 and 17.36 to 12.84, 14.18 and 15.90 THz, and the corresponding PCR values decrease from 91.16%, 99.78% and 97.69% to 80.04%, 99.40% and 94.96%. The operating frequency range of the LTL polarization converter (PCR is over 80%) has a red shift from the frequency range (13.94 to 17.52 THz) to one (12.78 to 16.02THz). When the Fermi energy (Ef) decreases from 0.84 to 0.7 eV, the performance of wideband LTL turns down.

The red shift behavior of the GSPs’ resonant frequency with decreasing Fermi energy can be interpreted by the following way. The wave vector of surface plasmon is satisfied with the equation

$$K{}_{spp} = \frac{{\hbar {\omega ^2}}}{{2{\alpha _0}{E_f}c}} \propto \frac{1}{{{L_g}}},{\alpha _0} = \frac{{{e^2}}}{{\hbar c}}$$

α0 is the fine structure constant. The resonant frequency f can be written as

$$f = \frac{\omega }{{2\pi }} \propto \sqrt {\frac{{{\alpha _0}c{E_f}}}{{2{\pi ^2}\hbar {L_g}}}} \propto \sqrt {\frac{{{E_f}}}{{{L_g}}}}$$

Lg is the resonance characteristic length of graphene carved-hollows. Equation (8) shows that decreasing Fermi energy (Ef) can cause red shift behavior. The dynamically tunable characteristic makes our proposed structure more practical than metal integrated metasurface polarization converter, which can be achieved by adjusting Fermi energy rather than geometric parameters.

Secondly, we change the electron scattering time τ to investigate the influence on the performances of the LTL polarization converter. Previous studies have shown that electron scattering time τ can vary over a wide range of ps or sub-ps [51]. Figure 6(b) shows PCR at different electron scattering time τ. It is obvious that as the electron scattering time τ decreases from 1 to 0.2 ps, the frequency position of the PCR peak is almost unchanged, but the PCR values of the three peaks will gradually decrease. When the electron scattering time τ is less than 0.6 ps, the PCR values will drop rapidly, and the broadband performance of LTL converter turns down. This is because a decrease in the scattering time will increase the loss, resulting in a decrease in the amplitude of the reflected wave.

 figure: Fig. 6.

Fig. 6. Simulation of the PCR under different (a) Fermi energy (Ef) varying from 1 to 0.7 eV and (b) electron scattering time τ varying from 1 to 0.2ps. Other parameters are the same with those of Fig. 2.

Download Full Size | PDF

As shown in Fig. 7, we find that the resonant frequencies of the symmetric GSPs mode blue shift and those of two antisymmetric GSPs modes are almost unchanged, when increasing the medium height m. Besides, as the medium height increases from 3600 to 4500 nm, the PCR at the resonant frequencies of the two antisymmetric GSPs modes gradually decreases, which leads to the possible emergence of the double LTC polarization conversion bands. In the next part, we will continue a further discussion about the double LTC bands.

 figure: Fig. 7.

Fig. 7. Simulation of the PCR when the medium height m varying from 3600 to 4500 nm. Other parameters are the same with those of Fig. 2.

Download Full Size | PDF

4. Double LTC polarization bands with opposite handedness

Based on the results of Fig. 7, when we set the medium height as m( = 4485 nm), and other parameters are as shown in Fig. 2, the perfect line-to-right-circular-polarization (LTRCP) and line-to-left-circular-polarization (LTLCP) conversion are simultaneously achieved. We use ellipticity χ and axial ratio to measure the performance of LTC polarization converter. Figure 8 shows the performance of LTC polarization converter under normal incident angle. Excited by the incident THz wave polarized at 45° with respect to x-direction, our proposed structure converts linear to right-handed and left-handed circular polarized waves. Figure 8(a) shows co-polarization (RAA) and cross-polarization (RBA) reflections in black and red lines, respectively. As shown in Fig. 8(b), ${\varphi _{BA}}$, ${\varphi _{AA}}$ and $\mathrm{\Delta }\varphi $($= {\varphi _{\textrm{BA}}} - {\varphi _{\textrm{AA}}}$) are plotted in black dashed line, blue dashed line and red solid line, respectively. RAA, and RBA are equal at 14.88 and 17.32 THz, while $\mathrm{\Delta }\varphi $ is equal to −90° and 90°, respectively. So, the perfect right-hand and left-hand circular polarization conversion are simultaneously formed at 14.88 and 17.32 THz in the reflection mode, respectively. Figure 8(c) shows that the band of LTRCP conversion is obtained at the lower frequency range, and the band of LTLCP conversion could be realized at the higher frequency range. Results indicate that the two bands of LTRCP and LTLCP operate in the frequency range of 13.82-16.08 THz and 16.84-17.74 THz with an axial ratio ($\textrm{ = 10log}|{{\textrm{r}_{\textrm{BA}}}\textrm{/}{\textrm{r}_{\textrm{AA}}}} |$) of less than 3 dB, respectively. The fractional bandwidths reach 15.1% and 5.2%, respectively. Besides, the ellipticity χ is also calculated in Fig. 8(c) to measure the bandwidth of the two bands of LTC polarization conversion. It can be found that χ is smaller than −0.95 from 14.24 to 15.74 THz, which indicates that LTRCP conversion is achieved and the fractional bandwidth is 10.0%. χ is −1 at 14.88 THz, which indicates that perfect LTRCP conversion is achieved. From 17.02 to 17.48 THz, χ is larger than 0.95, which indicates that LTLCP is achieved and the fractional bandwidth is 2.7%. χ is 1 at 17.32 THz, which indicates that perfect LTLCP is achieved. Obviously, the bandwidth defined by χ is a little smaller than the one defined with axial ratio.

 figure: Fig. 8.

Fig. 8. (a) Simulation results of the co-polarization reflection (RAA), cross-polarization reflection (RBA). (b) Phase of φAA and φBA together with their phase difference Δφ. (c) Ellipticity and axial ratio. The medium height m is 4485 nm, and other parameters are the same with those of Fig. 2.

Download Full Size | PDF

We further explore the influence of different geometric parameters L1, L2 and d on the performances of two bands of LTC polarization converter with opposite handedness, including ellipticity χ, phase difference Δφ, and the axial ratio. As shown in Figs. 9(a)-(c), the band of LTLCP conversion disappears when L1 is larger than 505 nm, because χ is smaller than 0.95. The band of LTLCP (χ>0.95) narrows from the frequency range (17.06 to 17.80 THz) to one (17.02 to 17.48 THz), and the band with 3 dB bandwidth of LTLCP narrows from frequency range (16.90 to 18.00 THz) to one (16.84 to 17.74 THz), with increasing L1 from 485 to 505 nm. Besides, the bandwidth of LTRCP is almost unchanged. It is because that there is a red shift of the third GSPs mode, and the first and second GSPs modes are almost unchanged with the increase of L1 from 465 to 545 nm (shown as Fig. 4(a)). It can be found that the band of LTRCP broadens, and the one of LTLCP is almost unchanged with increasing L2, as shown in Figs. 9(d)-(f). The band with |χ|> 0.95 of LTRCP broadens from the frequency range (14.94 to 15.92THz) to one (13.62 to 15.62THz), and the band with 3 dB bandwidth of LTRCP broadens from frequency range (14.52 to 16.22THz) to one (13.20 to 15.98THz), with increasing L2 from 625 to 705 nm. It is because that there is a red shift of the first GSPs mode, and the second and third GSPs mode are almost unchanged with the increase of L2 from 625 to 705 nm (shown as Fig. 4(b)). As shown in Figs. 9(g)-(i), it is found that the band of LTRCP broadens, and the one of LTLCP narrows with increasing d. The band with |χ|> 0.95 of LTRCP broadens from the frequency range (14.50 to 14.88 THz) to one (14.00 to 16.72 THz), and the one of LTLCP narrows from frequency range (16.64 to 17.20 THz) to one (17.54 to 17.74 THz), with increasing d from 50 to 80 nm. In the similar way, the band with 3 dB bandwidth of LTRCP broadens from frequency range (13.74 to 15.44 THz) to one (13.78 to 16.86 THz), and the one of LTLCP narrows from frequency range (16.48 to 17.52 THz) to one (17.40 to 17.94 THz). It is because that there is a blue shift of the second GSPs mode, and the first and third GSPs modes are almost unchanged with the increase of d from 50 to 80 nm (shown as Fig. 4(c)).

 figure: Fig. 9.

Fig. 9. Simulated (a), (d), (g) ellipticity, (b), (e), (h) phase difference Δφ, and (c), (f), (i) the axial ratio of the LTC polarization converter under different L1, L2 and d, respectively. Other parameters are the same with those of Fig. 8.

Download Full Size | PDF

Figure 10 shows performances of the double LTC polarization bands with opposite handedness under different incident angles from 0° to 70°, and other parameters are the same with those of Fig. 8. Figure 10(a) shows the band with |χ|> 0.95 of LTRCP narrows from the frequency range (14.24 to 15.74 THz) to one (14.44 to 15.68 THz), and the one of LTLCP narrows from frequency range (17.02 to 17.48 THz) to one (17.04 to 17.42 THz), with the incident angle increasing from 0° to 30°. When the incident angle is larger than 30°, the double bands of LTC polarization are destroyed due to |χ|< 0.95. Figure 10(b) shows that the band with 3 dB bandwidth of LTRCP narrows from frequency range (13.82 to 16.08 THz) to one (14.28 to 15.58 THz), and the one of LTLCP narrows from frequency range (16.84 to 17.74 THz) to one (16.84 to 17.44 THz), with the increase of incident angle from 0° to 40°. When the incident angle is larger than 40°, the double LTC polarization bands with 3 dB bandwidth are destroyed.

 figure: Fig. 10.

Fig. 10. (a) The ellipticity and (b) axial ratio of the proposed polarization converter under the incident angle ranging from 0°to 70°. Other parameters are the same with those of Fig. 8.

Download Full Size | PDF

With the decrease of Fermi energy (Ef), LTRCP and LTLCP bands both have an obvious red shift, and it is because the three GSPs’ modes have a red shift with decreasing Fermi energy based on Eq. (9). As shown in Fig. 11(a), the frequency range of LTRCP band with |χ|> 0.95 covers from 11.84 to 15.74 THz, the one of LTLCP covers from 14.02 to 17.48 THz with |χ|>0.95 with the Fermi energy (Ef) decreasing from 1 to 0.7 eV. In a similar way, as shown in Fig. 11(b), the 3dB band of LTRCP covers from 11.54 to 16.08 THz, the one of LTLCP covers from 14.02 to 17.74 THz. Thus, by changing the Fermi energy, the two broad bands of LTC are also achieved. Especially, the simulation results demonstrate the conversion between LTLCP and LTRCP can be achieved by changing Ef. For example, when Ef=0.7eV, the working frequency 14.1 THz is located at the LTLCP band (χ> 0.95 and axial ratio<3dB); when Ef=1eV, the working frequency 14.1 THz is located at the LTRCP band (χ<−0.95 and axial ratio<3dB).

 figure: Fig. 11.

Fig. 11. Simulation of the (a) ellipticity, and (b) the axial ratio of the LTC polarization converter under different Fermi energy (Ef) varying from 1 to 0.7 eV. Other parameters are the same with those of Fig. 8.

Download Full Size | PDF

With the decrease of electron scattering time τ, the performances (ellipticity and the axial ratio) of the LTC polarization converter are investigated as shown in Fig. 12. It is obvious that as the electron scattering time τ decreases from 1 to 0.2 ps, the two LTC bands have no frequency shift when other parameters are the same with those of Fig. 8. The bandwidths for the two LTC bands are almost unchanged, and the perfect right-hand and left-hand circular polarization always happens at 14.88 and 17.32 THz, respectively. Three peaks of axial ratio decrease with decreasing electron scattering time, which originates from increased free-carrier loss in graphene. When the electron scattering time τ is smaller than 0.2 ps, the performance of double bands of LTC polarization conversion turns down.

 figure: Fig. 12.

Fig. 12. Simulation of the (a) ellipticity, and (b) the axial ratio of the LTC polarization converter under different electron scattering time τ varying from 1.0 to 0.2ps. Other parameters are the same with those of Fig. 8.

Download Full Size | PDF

5. Conclusion

In summary, we have proposed a novel design of a tunable multifunctional polarization converter, which consists of a single graphene layer decorated with a π-shaped carved-hollow array. The broadband LTL and double bands of LTC polarization converter with opposite handedness can be realized respectively, by carefully choosing the medium height. The polarization conversions, including broadband LTL and double bands of LTC, are all related to the excitation of three GSPs modes. By increasing the Fermi energy, the switch between LTRCP and LTLCP can be realized, and the bands of polarization conversions all blue shift. In addition, the performance of the multifunctional polarization converter is also investigated under different geometrical parameters, incident angles and electron scattering time. This design provides a feasible method for the tunable multifunctional polarization control of terahertz field.

Funding

National Natural Science Foundation of China (1148081606193050); National College Students Innovation and Entrepreneurship Training Program (202010295062); The Key Research and Development Program of Jiangsu Province (BE2020756).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. B. Ferguson and X. C. Zhang, “Materials for terahertz science and technology,” Nat. Mater. 1(1), 26–33 (2002). [CrossRef]  

2. F. Zangeneh-Nejad and R. Safian, “A graphene-based THz ring resonator for label-free sensing,” IEEE Sens. J. 16(11), 4338–4344 (2016). [CrossRef]  

3. J. P. Fan and Y. Z. Cheng, “Broadband high-efficiency cross-polarization conversion and multi-functional wavefront manipulation based on chiral structure metasurface for terahertz wave,” J. Phys. D: Appl. Phys. 53(2), 025109 (2020). [CrossRef]  

4. W. L. Chan, M. L. Moravec, R. G. Baraniuk, and D. M. Mittleman, “Terahertz imaging with compressed sensing and phase retrieval,” Opt. Lett. 33(9), 974–976 (2008). [CrossRef]  

5. X. J. 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]  

6. I. A. Gayduchenko, G. E. Fedorov, M. V. Moskotin, D. I. Yagodkin, S. V. Seliverstov, G. N. Goltsman, A. Yu Kuntsevich, M. G. Rybin, E. D. Obraztsova, V. G. Leiman, M. S. Shur, T. Otsuji, and V. I. Ryzhii, “Manifestation of plasmonic response in the detection of sub-terahertz radiation by graphene-based devices,” Nanotechnology 29(24), 245204 (2018). [CrossRef]  

7. C. P. Huang, Q. J. Wang, X. G. Yin, Y. Zhang, J. Q. Li, and Y. Y. Zhu, “Break through the limitation of Malus’ Law with plasmonic polarizers,” Adv. Opt. Mater. 2(8), 723–728 (2014). [CrossRef]  

8. R. H. Fan, Y. Zhou, X. P. Ren, R. W. Peng, S. C. Jiang, D. H. Xu, X. Xiong, X. R. Huang, and M. Wang, “Freely tunable broadband polarization rotator for terahertz waves,” Adv. Mater. 27(7), 1201–1206 (2015). [CrossRef]  

9. S. Nersisyan, N. Tabiryan, D. M. Steeves, and B. R. Kimball, “Axial polarizers based on dichroic liquid crystals,” J. Appl. Phys. 108(3), 033101 (2010). [CrossRef]  

10. C. Y. Chen, T. R. Tsai, C. L. Pan, and R. P. Pan, “Room temperature terahertz phase shifter based on magnetically controlled birefringence in liquid crystals,” Appl. Phys. Lett. 83(22), 4497–4499 (2003). [CrossRef]  

11. L. Zhang, H. Zhong, C. Deng, C. Zhang, and Y. Zhao, “Characterization of birefringent material using polarization-controlled terahertz spectroscopy,” Opt. Express 18(19), 20491–20497 (2010). [CrossRef]  

12. C. C. Chang, Z. X. Zhao, D. F. Li, A. J. Taylor, S. H. Fan, and H. T. Chen, “Broadband linear-to-circular polarization conversion enabled by birefringent off-resonance reflective metasurfaces,” Phys. Rev. Lett. 123(23), 237401 (2019). [CrossRef]  

13. Y. M. Yang, W. Y. Wang, P. Moitra, I. I. Kravchenko, D. P. Briggs, and J. Valentine, “Dielectric meta-reflectarray for broadband linear polarization conversion and optical vortex generation,” Nano Lett. 14(3), 1394–1399 (2014). [CrossRef]  

14. L. L. Huang, X. Z. Chen, H. Muhlenbernd, G. X. Li, B. F. Bai, Q. F. Tan, G. F. Jin, T. Zentgraf, and S. Zhang, “Dispersionless phase discontinuities for controlling light propagation,” Nano Lett. 12(11), 5750–5755 (2012). [CrossRef]  

15. S. X. Xia, X. Zhai, Y. Huang, J. Q. Liu, L. L. Wang, and S. C. Wen, “Multi-band perfect plasmonic absorptions using rectangular graphene gratings,” Opt. Lett. 42(15), 3052–3055 (2017). [CrossRef]  

16. J. Y. Liu, Z. C. Li, W. W. Liu, H. Cheng, S. Q. Chen, and J. G. Tian, “High-Efficiency mutual dual-band asymmetric transmission of circularly polarized waves with few-layer anisotropic metasurfaces,” Adv. Opt. Mater. 4(12), 2028–2034 (2016). [CrossRef]  

17. L. Q. Jing, Z. J. Wang, Y. H. Yang, L. Shen, B. Zheng, F. Gao, H. P. Wang, E. P. Li, and H. S. Chen, “Diodelike spin-orbit interactions of light in chiral metasurfaces,” IEEE Trans. Antennas Propag. 66(12), 7148–7155 (2018). [CrossRef]  

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

19. J. M. Hao, Y. Yuan, L. X. Ran, T. Jiang, J. A. Kong, C. T. Chan, and L. Zhou, “Manipulating electromagnetic wave polarizations by anisotropic metamaterials,” Phys. Rev. Lett. 99(6), 063908 (2007). [CrossRef]  

20. S. J. Li, Y. B. Li, L. Zhang, Z. J. Luo, B. W. Han, R. Q. Li, X. Y. Cao, Q. Cheng, and T. J. Cui, “Programmable controls to scattering properties of a radiation array,” Laser Photonics Rev. 15(2), 2000449 (2021). [CrossRef]  

21. S. J. Li, Y. B. Li, H. Li, Z. X. Wang, C. Zhang, Z. X. Guo, R. Q. Li, X. Y. Cao, Q. Cheng, and T. J. Cui, “A thin self-feeding janus metasurface for manipulating incident waves and emitting radiation waves simultaneously,” Ann. Phys. 532(5), 2000020 (2020). [CrossRef]  

22. W. L. Barnes, A. Dereux, and T. W. Ebbesen, “Surface plasmon subwavelength optics,” Nature 424(6950), 824–830 (2003). [CrossRef]  

23. L. Billings, “Metamaterial world,” Nature 500(7461), 138–140 (2013). [CrossRef]  

24. X. F. Zang, S. J. Liu, H. H. Gong, Y. J. Wang, and Y. M. Zhu, “Dual-band superposition induced broadband terahertz linear-to-circular polarization converter,” J. Opt. Soc. Am. B 35(4), 950–957 (2018). [CrossRef]  

25. Y. Jiang, H. Zhao, and L. Wang, “Broadband linear-to-circular polarization converter based on phosphorene metamaterial,” Opt. Mater. Express 9(5), 2088 (2019). [CrossRef]  

26. B. Q. Lin, J. L. Wu, and X. Y. Da, “A linear-to-circular polarization converter based on a second-order band-pass frequency selective surface,” Appl. Phys. A: Mater. Sci. Process. 123(1), 43 (2017). [CrossRef]  

27. Z. Y. Li, S. J. Li, B. W. Han, G. S. Huang, Z. X. Guo, and X. Y. Cao, “Quad-Band transmissive metasurface with linear to dual-circular polarization conversion simultaneously,” Adv. Theory Simul. 4(8), 2100117 (2021). [CrossRef]  

28. B. W. Han, S. J. Li, Z. Y. Li, G. S. Huang, J. H. Tian, and X. Y. Cao, “Asymmetric transmission for dual-circularly and linearly polarized waves based on a chiral metasurface,” Opt. Express 29(13), 19643–19654 (2021). [CrossRef]  

29. J. Xu, R. Q. Li, J. Qin, S. Y. Wang, and T. C. Han, “Ultra-broadband wide-angle linear polarization converter based on H-shaped metasurface,” Opt. Express 26(16), 20913–20919 (2018). [CrossRef]  

30. Y. Li, Q. S. Cao, and Y. Wang, “A wideband multifunctional multilayer switchable linear polarization metasurface,” Antennas Wirel. Propag. Lett. 17(7), 1314–1318 (2018). [CrossRef]  

31. J. Xu, R. Q. Li, J. Qin, S. Y. Wang, and T. C. Han, “Ultra-broadband linear polarization converter based on anisotropic metasurface,” Opt. Express 26(20), 26235–26241 (2018). [CrossRef]  

32. J. K. Gansel, M. Latzel, A. Frölich, J. Kaschke, M. Thiel, and M. Wegener, “Tapered gold-helix metamaterials as improved circular polarizers,” Appl. Phys. Lett. 100(10), 101109 (2012). [CrossRef]  

33. L. Wu, Z. Y. Yang, Y. Z. Cheng, M. Zhao, R. Z. Gong, Y. Zheng, J. A. Duan, and X. H. Yuan, “Giant asymmetric transmission of circular polarization in layer-by-layer chiral metamaterials,” Appl. Phys. Lett. 103(2), 021903 (2013). [CrossRef]  

34. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, “Electric field effect in atomically thin carbon films,” Science 306(5696), 666–669 (2004). [CrossRef]  

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

36. C. H. Gan, H. S. Chu, and E. P. Li, “Synthesis of highly confined surface plasmon modes with doped graphene sheets in the midinfrared and terahertz frequencies,” Phys. Rev. B 85(12), 125431 (2012). [CrossRef]  

37. F. Bonaccorso, Z. Sun, T. Hasan, and A. C. Ferrari, “Graphene photonics and optoelectronics,” Nat. Photonics 4(9), 611–622 (2010). [CrossRef]  

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

39. G. W. Hanson, “Dyadic Green’s functions for an anisotropic, non-local model of biased graphene,” IEEE Trans. Antennas Propag. 56(3), 747–757 (2008). [CrossRef]  

40. A. K. Geim, “Graphene: status and prospects,” Science 324(5934), 1530–1534 (2009). [CrossRef]  

41. K. F. Mak, L. Ju, F. Wang, and T. F. Heinz, “Optical spectroscopy of graphene: from the far infrared to the ultraviolet,” Solid State Commun. 152(15), 1341–1349 (2012). [CrossRef]  

42. Y. Jiang, H. D. Zhang, J. Wang, C. N. Gao, J. Wang, and W. P. Cao, “Design and performance of a terahertz absorber based on patterned graphene,” Opt. Lett. 43(17), 4296–4299 (2018). [CrossRef]  

43. Z. F. Yao, T. T. Wei, Y. K. Wang, M. J. Lu, C. Y. Zhang, and L. L. Zhang, “Tunable multifunctional reflection polarizer based on a graphene metasurface,” Appl. Opt. 58(13), 3570–3574 (2019). [CrossRef]  

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

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

46. M. Sajjad, X. K. Kong, S. B. 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]  

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

48. S. X. He, R. Tian, W. Wu, W. D. Li, and D. Q. Wang, “Helium-ion-beam nanofabrication: extreme processes and applications,” Int. J. Extreme Manuf. 3(1), 012001 (2021). [CrossRef]  

49. A. N. Abbas, G. Liu, B. L. Liu, L. Y. Zhang, H. Liu, D. Ohlberg, W. Wu, and C. W. Zhou, “Patterning, characterization, and chemical sensing applications of graphene nanoribbon arrays down to 5 nm using helium ion beam lithography,” ACS Nano 8(2), 1538–1546 (2014). [CrossRef]  

50. L. Cong, W. Cao, X. Q. Zhang, Z. Tian, J. Q. Gu, R. J. Singh, J. G. Han, and W. L. Zhang, “A perfect metamaterial polarization rotator,” Appl. Phys. Lett. 103(17), 171107 (2013). [CrossRef]  

51. P. Neugebauer, M. Orlita, C. Faugeras, A. L. Barra, and M. Potemski, “How perfect can graphene be?” Phys. Rev. Lett. 103(13), 136403 (2009). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable 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 (12)

Fig. 1.
Fig. 1. (a) Schematic diagram of multifunctional polarization converter, in which yellow and gray regions represent gold and silicon dioxide layer, respectively. (b) One unit of (a). Geometric parameters include: L1 = 505 nm, L2=665 nm, d=60 nm, g=110 nm and Px= Py = 990 nm. Besides, the incident angle is defined as θ as shown in Fig. 1(c).
Fig. 2.
Fig. 2. The simulation results of (a) co-polarization and cross-polarization reflection (RAA and RBA), and (b) PCR are plotted in black, blue and black solid line, respectively.
Fig. 3.
Fig. 3. (a) The co-polarized reflection spectrum under x-polarized Ex and y-polarized Ey incidence, respectively. The simulation results of Rxx and Ryy are drawn with blue and green solid lines. The dotted red line represents the PD between Rxx and Ryy. Magnetic field distributions (Hz) of 13.84 THz(b) and 17.58 THz (d) under Ey incidence, and 16.00 THz (c) under Ex incidence, respectively.
Fig. 4.
Fig. 4. Simulated PCR of the LTL polarization converter under normal incidence with different (a) length L1, (b) length L2 and (c) width d of the π-shaped hollow array. Other parameters are the same with those of Fig. 2.
Fig. 5.
Fig. 5. Simulation of the PCR under different incident angles from 0° to 70°. Other parameters are the same with those of Fig. 2.
Fig. 6.
Fig. 6. Simulation of the PCR under different (a) Fermi energy (Ef) varying from 1 to 0.7 eV and (b) electron scattering time τ varying from 1 to 0.2ps. Other parameters are the same with those of Fig. 2.
Fig. 7.
Fig. 7. Simulation of the PCR when the medium height m varying from 3600 to 4500 nm. Other parameters are the same with those of Fig. 2.
Fig. 8.
Fig. 8. (a) Simulation results of the co-polarization reflection (RAA), cross-polarization reflection (RBA). (b) Phase of φAA and φBA together with their phase difference Δφ. (c) Ellipticity and axial ratio. The medium height m is 4485 nm, and other parameters are the same with those of Fig. 2.
Fig. 9.
Fig. 9. Simulated (a), (d), (g) ellipticity, (b), (e), (h) phase difference Δφ, and (c), (f), (i) the axial ratio of the LTC polarization converter under different L1, L2 and d, respectively. Other parameters are the same with those of Fig. 8.
Fig. 10.
Fig. 10. (a) The ellipticity and (b) axial ratio of the proposed polarization converter under the incident angle ranging from 0°to 70°. Other parameters are the same with those of Fig. 8.
Fig. 11.
Fig. 11. Simulation of the (a) ellipticity, and (b) the axial ratio of the LTC polarization converter under different Fermi energy (Ef) varying from 1 to 0.7 eV. Other parameters are the same with those of Fig. 8.
Fig. 12.
Fig. 12. Simulation of the (a) ellipticity, and (b) the axial ratio of the LTC polarization converter under different electron scattering time τ varying from 1.0 to 0.2ps. Other parameters are the same with those of Fig. 8.

Equations (9)

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

ε=1 + iσgη0k0Δ
σg=ie2Efπ2(ω+iτ1)+ie24πln[2Ef(ω+iτ1)2Ef+(ω+iτ1)]+ie2kBTπ2(ω+iτ1)ln[exp(EfkBT)+1]
(EArEBr)=r(EAiEBi)
r=(rAArABrBArBB)
PCR=RBARAA+RBA
RAA=|rAA|2,RBA=|rBA|2
S0=|rAA|2+|rBA|2S1=|rAA|2|rBA|2S2=2|rAA||rBA|cos(Δφ)S3=2|rAA||rBA|sin(Δφ)
Kspp=ω22α0Efc1Lg,α0=e2c
f=ω2πα0cEf2π2LgEfLg
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.