1. Introduction
Chirality is an inherent feature of chiral molecules. That is, an object cannot coincide with its mirror image through translation, turnover and other in-plane operations. Similar to our left and right hands, they mirror each other and cannot overlap. Chiral objects and their mirror images are called enantiomers. In the organic world, there are many chiral molecules, such as protein, DNA, sugar molecules, viruses, amino acids and liquid crystals, but there are few chiral structures in nature. Later, it is found that chiral materials can control the spin of electrons to a certain extent [
1,
2,
3]. In addition, thanks to the development of metamaterials, the control of the spin of photons has been realized [
4,
5,
6,
7,
8,
9,
10]. At present, the study of the interaction between circularly polarized light (CPL) and matter has been involved in many wavebands, such as microwave [
11], terahertz [
12], infrared [
13] and visible [
14]. It is well known that chiral metamaterials show different electromagnetic responses such as asymmetric transmission (AT) and circular dichroism (CD) when they interact with CPL with different spins. In other words, chiral metamaterials have selective transmission, reflection and absorption of incident CPL with different spins.
In recent years, metamaterials with spin selectivity have played a very important role in many fields. For example, spin-selective beam deflectors [
15], polarization-sensitive chiral circularly polarized light detectors [
16], circularly polarized light multiplexed holograms [
17] and spin- and wavelength-sensitive gray imaging [
18]. However, the optical effect of metamaterials is difficult to change after manufacturing, which limits the practical application of functional metamaterials. Hence, actively tunable metamaterials are gaining increasing attention. At present, there are two design methods that can make metamaterials have the property of active regulation. The first method is to change the geometry of metamaterials. Hu Tao et al. proposed a reconfigurable anisotropic metamaterial at terahertz frequency, which can bend the structure through thermal stimulation to change the structural response [
19]. Yuan Hsing Fu et al. designed a micromechanically reconfigurable metamaterial. Dynamic tuning can be realized by changing the coupling distance between cells through bidirectional micromechanical actuators [
20]. The second method is to add controllable materials such as silicon, graphene and phase change materials into metamaterials. Lei Kang et al. presented a nonlinear chiral meta-mirror based on silicon, which can be used for an all-optical picosecond optical polarization switch. The chiral optical response of the structure is dynamically adjusted by the change of carriers in silicon by pump light [
21]. Chuanchuan Ding et al. proposed an absorber integrating a periodic Ge
2Sb
2Te
5 (GST) resonator. Through the phase transition of the GST resonator from an amorphous state to a crystalline state, a change in CD intensity in the range of 0.03 to 0.75 can be realized [
22].
In contrast to the past, some researchers have shifted their attention from tuning the chiral response intensity to switching the chiral response. Meng Liu et al. took advantage of vanadium dioxide (VO
2) to change the chirality of metamaterial resonators and realized a change in CD intensity from −75% to +55% [
23]. Xiaoqing Luo et al. proposed a dynamically reversible chiral metamaterial by combining vanadium dioxide and Babinet’s principle, which can achieve a reversible CD tuning range of ±0.5 at 0.97 THz [
24]. Yongkai Wang et al. realized a reversible chiral response in the terahertz band by making use of achiral dielectric and graphene bands. Chiral response switching can be realized by reversing the Fermi level of the graphene band [
25]. From the perspective of practical application, the design of a dynamically reversible metamaterial with strong circular dichroism is more in line with the application requirements.
Here, we theoretically investigate a switchable chiral metasurface based on a phase change material GST, which can switch between a right-handed circularly polarized (RCP) mirror and a left-handed circularly polarized (LCP) mirror. The metasurface is composed of asymmetric split-ring resonators (ASRRs) arranged periodically. The simulated results show that when the conductivity of GST σ is equal to 0 S/m, ASRR_1 will reflect incident right-handed circularly polarized light (RCPL) and absorb incident left-handed circularly polarized light (LCPL) at 0.76 THz. As σ rises to 3 × 105 S/m, the response of ASRR_1 to circularly polarized light will be reversed. That is, it reflects the incident LCPL and absorbs the incident RCPL at 0.66 THz. In addition, we also simulate ASRR_2, the mirror structure of ASRR_1. Then, the effect of conductivity on the properties of ASRR_1 was also studied. In addition, the physical principle of ASRR_1 was explained. Finally, we arrange ASRR_1 into a phase continuous metasurface. Before and after the GST phase transition, it can anomalously reflect the incident LCPL and RCPL out of the whole metasurface at different angles.
2. Materials and Methods
The performance of the switchable metasurface was numerically simulated by a commercial full-wave numerical software computer simulation technology (CST) Microwave Studio. The frequency-domain solver is selected as the simulated method. On a single periodic structural element, the boundary conditions in the x and y directions are set as unit cells, and the boundary condition in the z direction is open. When x- and y-polarized plane waves are incident, due to the existence of polarization conversion, we obtain the following reflection coefficients of linearly polarized light:
rxx,
ryx,
ryy and
rxy. Then, the reflection coefficient of CPL can be obtained by Equation (1). Here,
x,
y,
r and
l in the subscript represent
x-polarized light,
y-polarized light, right-handed circularly polarized light and left-handed circularly polarized light, respectively.
Rij represents the reflection coefficient of
i-polarized light reflected by the structure under the incidence of
j-polarized light,
i,
j∊{
x,
y,
r,
l}. It should be pointed out that we will observe circularly polarized light facing the propagation direction of electromagnetic waves. At a fixed position, if the polarization direction rotates clockwise with time, it is RCPL; in contrast, it is LCPL.
Figure 1a,b shows a schematic diagram of the unit structure of ASRR_1. Here, a continuous uniform metal film
rr with a thickness of 200 nm is used as the reflecting layer. The top of each unit is an ASRR composed of metal and phase change material GST. A 50 μm thick polyimide (PI) dielectric layer separates the bottom metal film from the top ASRR. In the simulation, the metal structures were made of aluminum (Al) with a conductivity of 3.56 × 10
7 S/m. The dielectric constant and loss tangent of PI are 2.93 and 0.044, respectively. The thickness, outer radius R and linewidth w of ASRR_1 are 200 nm, 65 μm and 20 μm, respectively. The angle of gaps α and β are 20° and 40°. φ is the angle at which ASRR_1 rotates around the z-axis. The period of the unit structure P is 160 μm.
Figure 1c shows the schematic of the electromagnetic responses of a gradient chiral metamirror, which is obtained by arranging ASRR_1 into a phase continuous metasurface.
3. Results
First, we simulate the reflection spectra of ASRR_1 and ASRR_2, which are mirror images of each other.
Figure 2a–d shows the unit structure diagram and the reflection spectrum results of ASRR_1. As shown in
Figure 2b, when the conductivity of GST σ is 0 S/m, ASRR_1 will efficiently absorb the incident LCPL at a frequency of 0.76 THz, while it will strongly reflect the incident RCPL. That is,
Rrr = 0.76 and
Rll = 0. At this time, ASRR_1 is called RCP mirror. As σ rises to 3 × 10
5 S/m, it can be seen from
Figure 2c that ASRR_1 is switched to the LCP mirror. That is, it can efficiently absorb the incident RCPL and strongly reflect the incident LCPL. The reflectivity is
Rrr = 0.02 and
Rll = 0.66. Since GST after the phase transition increases the electrical length of the whole structure, the working frequency of ASRR_1 shifts red to 0.66 THz. The quantitative calculation of circular dichroism is CD =
Rrr −
Rll. In
Figure 2d, before the phase transition, the circular dichroism of ASRR_1 is CD = 76%. After the phase transition, the circular dichroism of ASRR_1 is CD = −64%. This also shows that switching between the RCP mirror and LCP mirror is realized by using the phase transition property of GST.
Figure 2e–h shows the unit structure diagram and the reflection spectrum results of ASRR_2. As ASRR_2 and ASRR_1 are mirror image relationships in structure, their reflection spectrum and CD spectrum also show a reversal effect. That is, ASRR_2 is an LCP mirror before the GST phase transition and will switch to an RCP mirror after the GST phase transition. Its circular dichroism changes from CD = −76% to CD = 64%.
Next, there are several reflection spectra of ASRR_1 with obvious changes under different conductivities in
Figure 3. We have explained that the reflection effects of ASRR_1 and ASRR_2, which are mirror structures to each other, are reversed to each other in
Figure 2. Therefore, only the simulation results of ASRR_1 will be discussed in future research. With increasing conductivity, the reflection valley of LCPL at 0.76 THz is gradually passivated, and the reflectivity increases. Moreover, the reflectivity difference between the RCPL and LCPL will also gradually decrease. As σ = 6 × 10
4 S/m, the reflection valley of LCPL basically disappears. At this time, reflection valley of RCPL appeared at 0.66 THz. Moreover, the reflection valley of RCPL will gradually sharpen, and the reflectivity will decrease with increasing conductivity. At the same time, the reflectivity difference between the RCPL and LCPL also increases. When σ is 3 × 10
5 S/m, ASRR_1 completes the switching from the RCP mirror to the LCP mirror.
To explain the physical mechanism of this switchable chiral metasurface, the surface current distribution of ASRR_1 under different conductivities is simulated, as shown in
Figure 4. Here, we mark the current oscillation direction along the two arms of ASRR_1 with a white arrow. When the conductivity of GST is 0 S/m, ASRR_1 shows the effect of absorption of LCPL and reflection of RCPL at 0.76 THz. As seen from
Figure 4a, the surface current oscillates in a single direction on the two arms, which shows that the incident LCPL excites a magnetic dipole resonance in ASRR_1. This resonant mode can bind the incident electromagnetic wave energy inside the structure, and the lossy PI layer will slowly absorb the energy. Therefore, LCPL is highly absorbed after incident on ASRR_1. In
Figure 4b, the surface current oscillates in different directions along the two arms of the resonant ring, indicating that the incident RCPL excites an electric dipole with low loss in ASRR_1. Thus, the metasurface at this time presents the performance of reflecting RCP incident light. As σ is 3 × 10
5 S/m, ASRR_1 shows the effect of reflecting LCPL and absorbing RCPL at 0.66 THz. According to the current distribution in
Figure 4c, the incident LCPL excited in the structure is an electric dipole; thus, LCPL is reflected. Similarly, by analyzing
Figure 4d, it can be seen that the magnetic dipole that will cause high absorption is excited after the RCPL is incident on ASRR_1. Through the analysis, we can clearly understand that the different resonance modes excited by incidence lead to the different responses of the chiral metasurface to circularly polarized light with different spins under different conductivities.
According to the Pancharatnam–Berry phase [
26], when a structure rotates φ around the z-axis, the phase of the reflected light will change by 2φ, while the amplitude of the reflected light remains unchanged. Then, the reflection spectrum of ASRR_1 with different orientations φ is simulated. As shown in
Figure 5a, RCPL is incident to ASRR_1. The conductivity of GST is σ = 0 S/m. Regardless of the value of φ, the reflection remains almost unchanged and approximately equals 0.76. When φ changes from 0 to π at an interval of π/18, the phase of reflected RCPL changes between −π and π, that is, the phase change range is 2π. The operating frequency is 0.76 THz.
Figure 5b shows that ASRR_1 with σ = 3 × 10
5 S/m is irradiated by LCPL. Similarly, the amplitude of reflected LCPL, which is approximately 0.66, will not change with different φ. The phase of reflected LCPL depends on 2φ and changes in the range of 2π. The operating frequency here is 0.66 THz.
In the generalized Snell’s law, we can control the direction of the emitted light by designing a linear varied phase shift
along the surface of a metasurface. For reflected light, the equation of generalized Snell’s law is defined as [
27]:
where
and
represent the reflection angle and incident angle, respectively;
and
represent the refractive indices of the output and input media, which are equal to 1;
is the speed of light in vacuum;
corresponds to the operating frequency of the metasurface; and
stands for a constant phase gradient. PB phase is a method to accurately control
. That is
[
28]. φ here is the rotation angle of ASRR_1 around the z-axis we mentioned above. The sign plus or minus is determined by RCPL and LCPL. As seen from
Figure 5, ASRR_1 can change the phase of its reflected light by rotating around the z-axis, and the phase size can cover the whole phase period 2π. Therefore, we chose six different rotation angles: 0, π/6, π/3, π/2, 2π/3 and 5π/6. The six small units are arranged in a row to form a superperiodic unit cell and thus a chiral metasurface with linear phase change that covers the whole 2π range. Part of the structure of the chiral metasurface is shown in
Figure 1c. The period in the x direction is P
x = 960 μm, and the period in the y direction is P
y = 160 μm. Similarly, we use the simulation software CST to simulate the reflected field distribution of the metasurface.
Figure 6a,b show the results of LCPL and RCPL incident on the chiral metasurface with a conductivity of 0 S/m, respectively. It can be clearly seen that the RCPL incident vertically has obvious deflection after being reflected by the metasurface, that is, anomalous reflection. Because the metasurface under this conductivity will efficiently absorb the incident LCPL, the reflected electric field is almost invisible. In addition, we simulate the reflected electric field with a conductivity of 3 × 10
5 S/m. As shown in
Figure 6c,d, the results correspond to the LCP incidence and RCP incidence, respectively. When LCPL is incident vertically, it will be anomalously reflected off the chiral metasurface. However, the incident RCPL will be absorbed. The uniform wavefront in
Figure 6b,c shows that the amplitudes of circularly polarized light reflected by the six structures are almost equal, which is consistent with the results obtained in
Figure 5. Comparing the two figures, we can see that the angles of the chiral metasurface anomalously reflected under the two conductivities are different. There are two reasons for the difference: the first is the different working frequencies, and the second is the reflected light with different spins. Regardless of the conductivity, the linear phase distribution of the metasurface takes π/3 as the gradient. Thus, according to Equation (2), the two deflection angles are 24.3° and −28.3°.