1 Introduction
With the significant advancements in terahertz (THz) science and technology, terahertz waves have become indispensable in radar [
1], communication [
2,
3], and imaging [
4,
5]. These applications demand highly efficient optical devices, including lenses, wave plates, switches, and polarization rotators. Among these devices, polarization converters are crucial for controlling polarization and are typically constructed using optical gratings [
6], birefringent crystals [
7,
8] or anisotropic materials. However, traditional converters suffer from drawbacks such as large size, susceptibility to damage, and integration challenges [
9,
10]. Moreover, the conventional natural materials chosen for traditional converters have significant limitations in electromagnetic response. Consequently, researchers are considering the use of metamaterials as substitutes for traditional materials.
Metamaterials [
11,
12], as artificial composite materials, have gained substantial attention due to their unique electromagnetic properties, offering alternative possibilities for controlling the polarization of terahertz waves [
13–
17]. In 2018, Sun et al. demonstrated a transmissive metal polarization rotator with a double-layer structure, which, relative to pre-existing technologies, exhibited reduced loss and improved polarization conversion performance [
18]. Cheng et al. further utilized chiral metamaterials to design a versatile terahertz polarization converter, achieving linear-to-circular polarization conversion at 1.14 and 1.34 THz, as well as cross-linear polarization conversion between 2.19 and 2.47 THz [
19]. In 2020, Zi et al. proposed a transmissive dual-function terahertz wave plate based on all-dielectric metamaterials, enabling two different polarization conversions at 1.01 THz [
20]. These designs operate in the transmissive mode, thereby to some extent avoiding the associated drawbacks of traditional converters and facilitating easy integration. However, their effective working bandwidth remains relatively narrow. Therefore, achieving a transmissive broadband device with high polarization conversion efficiency is important for development of control of the polarization state of terahertz waves.
In this paper, we present a broadband all-silicon transmissive multi-functional polarization converter with high efficiency, employing cross-shaped microstructures. This design enhances the effective electromagnetic response, thereby improving the device efficiency. In comparison to previously reported converters with limited working bandwidth, our proposed design offers a significantly broader operational range. Numerical simulations confirm the device's capability to achieve polarization conversion for both cross-linear and linear-to-circular polarized waves in the terahertz frequency range. The proposed device holds potential applications in terahertz spectroscopy, imaging, and communications.
2 Theoretical model
The unit cell of this polarization convertor is composed of a cross-shaped patch on top of a dielectric substrate (schematic in Fig. 1b). When a beam of terahertz waves passes through this device, its polarization direction will change. Here,
Ei =
Eixe
−jkz is the incident polarized wave amplitude,
Eix is incident wave amplitude,
k is the wave number, and
Eix represents the amplitude of the incident wave along the
x-axis [
21]. Then the transmitted wave can be expressed
where Etx and φx represent the amplitude and additional phase of the transmitted wave, respectively.
If the cross-shaped microstructure in
u-v coordination rotates through angle
θ (
θ is the included angle between the cross-shaped minor axis and the
x-axis, as shown in Fig. 1b), the incident and transmitted wave amplitude in the
x-direction can be expressed as [
22]
The cross-shaped microstructure does not change θ, so for simplicity we usually consider that: φx = φu and φy = φv. The transmitted waves are then given by
Then the polarization state of transmitted waves is related to the rotated angle and the phase difference between the two orthogonal coordinate axes. The linear polarization direction of transmitted wave will rotate π/2 as long as the two polarization components along the
x- and
y-axis produce a phase difference of π [
23],
where
λ is the wavelength of incident wave,
nu and
nv are the refractive indices in
u- and
v-directions, respectively, and
h is the thickness of cross-shaped microstructure. Since
, the modulation of polarization state is related to permittivity [
24], which is a unique and significant characteristic of metamaterials. We can handle the polarization direction through artificially changing the shape, size and distribution of the unit cell [
25,
26]. It can be seen from Fig. 1b that the high-resistance silicon and air are arranged periodically. Therefore, the effective permittivity calculation formulas in the
u-direction and
v-direction can be expressed as
where ε1 and ε2 are the permittivities of high-resistance silicon and air, respectively, εeff is the effective permittivity, and f is the duty cycle of high resistance silicon. By appropriately adjusting the proportion of the high-resistance silicon and changing the refractive index of the convertor, the polarization direction of terahertz waves can be deflected toward the target angle.
From Eq. (4), the transmitted waves of x- and y-polarized components are also related to the rotated angle θ. When Δφ = φy − φx = π, the transmitted wave amplitudes are Et = −(cos 2θ + sin 2θ)eiφxejkz. Therefore, when θ = 45°, the x-polarized incident waves can be completely turned into y-polarized waves. According to Eq. (4), we set θ = 45° and the polarization direction of incident waves is then along the x-axis. The transmittance in x- and y-polarized directions are tx = Ex/Ei, ty = Ey/Ei. The polarization conversion rate (PCR) is defined as PCR = |ty|2/(|tx|2 + |ty|2). Here, Ex and Ey represent the output electrical field intensity of the terahertz waves polarized along x- and y-axis, respectively, Ei represents the electrical field intensity of incident terahertz waves.
To verify the design, the simulation is carried out in CST MWS, where the material is high-resistance silicon with ε = 11.9 and μ = 1.0. In the simulation, the periodic boundary condition in the x–y plane is applied to simulate the periodic array state. Meanwhile, the open condition along the z-direction is set to match the practical incident and transmitted terahertz wave paths. For the light source setup, we use a linearly polarized wave along the x-direction incident along the − z-axis to the metamaterial.
3 Results and discussion
3.1 Influence of parameters on device performance
There are several factors that affect the device performance: length, width, height, period, and substrate thickness. By varying these factors, we can obtain the influence of structural parameters on PCR, and the optimal structural parameters.
Table 1 shows the parameter scanning range of the cross-shaped microstructure. Through the parameter scanning, the optimal solution of PCR and the effective working bandwidth are obtained.
Figure 2 illustrates the simulated PCR distributions for various structural parameters. (Here, the red area indicates that the device has high polarization conversion efficiency in this frequency band). The red areas represent PCR ≥ 80%. These area are concentrated on the range of 1.00–2.5 THz. Figure 2a–e show the impact of the cross-shaped parameters on the effective working bandwidth and PCR. When a = 30 μm, c = 30 μm, the device performance is the best. Figure 2b, d, e indicate that the effective working bandwidth increases first and then drops off, with the maximum at b = 80 μm, d = 40 μm, h = 180 μm. As shown in Fig. 2f, the period has great influence on device efficiency. With the period increasing, the effective working frequency range continuously reduces, and reaches the best value at P = 80 μm. By analyzing the influence of microstructure parameters on the polarization conversion, the optimum parameter combination of the cross-shaped microstructure is as follows: a = 30 μm; b = 80 μm; c = 30 μm; d = 40 μm; h = 180 μm; P = 80 μm.
3.2 Function of broadband half-wave plate
According to the optimal structural parameters obtained by parameter scanning, the microstructure is then simulated. Figure 3 shows the transmittance and corresponding PCR of the device at the optimal parameters as determined by simulation. It can be understood that in the frequency range of 1.00–2.32 THz, PCR is greater than 80%, representing a very broad effective working bandwidth. In addition, Table 2 gives the performance of the proposed structure in comparsion with previous works. Compared with previously reported designs, our design for cross-polarization conversion has simple structure, wide operating bandwidth, and works well at wide incidence angles.
3.3 Function as a quarter-wave plate
In addition, this device can be used not only as a half-wave plate, but also as a quarter-wave plate. The quarter wave plate can convert linearly polarized waves into circularly or elliptically polarized waves. The expression of the quarter-wave plate can be expressed as
When the two orthogonal components of the transmitted waves are equal, the linearly polarized waves can be converted into circularly polarized waves. The Stokes formula describing ellipticity, χ = S1/S2, where S1 and S2 are the Stokes parameters, is used to evaluate the circular polarization effect. Here, S1 = 2|tx||ty|sinφ, S2 =|tx|2 + |ty|2, and φ is the phase difference. From Fig. 4, when the x-polarized incident waves impinge on the device, the transmitted waves are completely converted into left-handed circularly polarized waves in 0.85 THz.
In Fig. 4, a
y-polarized wave is incident into the device, then it will be converted into right-handed circularly polarized wave at 0.85 THz. Additionally, we use PB (Pancharatnam-Berry) phase [
32,
33] to change the polarization state of transmitted waves. By adjusting the spatial transformation of this device, it can also achieve a relatively effective modulation effect.
We also consider the angular dependence of the proposed polarization convertor. (The substrate thickness is set to 10 μm here). The numerically resolved transmittance and PCR are depicted in Fig. 4. As shown in Fig. 4a, b, the transmittance decreases with the increase of incident angle. While the polarization conversion performance (PCR) is minimally affected (Fig. 5c). These results show that the structure has good robustness, which is of great significance for practical applications.
Regarding the application situation of the proposed convertor, the influence of substrate thickness on device performance is analyzed. As shown in Fig. 4d–f, the red ranges represent that the transmittance and PCR are better than white areas. The redder the area, the higher the transmittance and PCR. Thus, the device is more easily integrated by reducing the thickness of substrate.
3.4 Optimization of polarization convertor
As shown in Fig. 5a–c, the incident angle seriously affects transmittance in
x- and
y-directions. With the angle increases, the transmittance decreases, which means that a lot of energy is lost in transmission. Therefore, it's important to improve the energy efficiency of the device. We plate an anti-reflection film on the surface of the device to increase the transmittance [
33]. Since the device adopts a metamaterial structure, the coating method is not applicable. A “moth-eye” (ME) structure is proposed to improve transmittance [
34–
36]. The ME structure is a periodic and regular arrangement of protrusions, as shown in Fig. 5a–b. It is always added at the bottom of the microstructure. Due to the protrusions on the surface, the reflected radiation is reduced, that is, the reflectivity is diminished and the transmittance is enhanced. Figure 5c shows the effect of the device on transmittance of terahertz waves after adding the ME structure. The ME structure with height
l = 25 μm and diameter
D = 35 μm have good influence on the performance of the device.
4 Conclusion
We present a new all-silicon metamaterial featuring a cross-shaped microstructure capable of simultaneously achieving cross-linear and linear-to-circular polarization conversion in the terahertz frequency range. By analyzing the scanning results, we determine the optimal structural parameters. Our device achieves a broadband conversion range of approximately 1.32 THz and demonstrates high efficiency, reaching up to 80% for orthogonal linear polarization conversion. Furthermore, we investigate the device's response at various incident angles and find that it exhibits no angular dependence. Additionally, we observe that the thickness of the substrate has minimal impact on the device's performance, confirming its excellent integration potential. To further enhance the transmittance of terahertz waves, we incorporate a magneto-electric (ME) structure. Overall, our device serves as a valuable reference for the design and theoretical analysis of transmissive broadband polarization converters, particularly in integrated terahertz systems.