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

Phase distribution and circular dichroism switchable terahertz chiral metasurface

Open Access Open Access

Abstract

Chiral metasurfaces have many applications in the terahertz (THz) band, but they still lack modulation flexibility and functionality expansion. This paper presents a terahertz chiral metasurface with switchable phase distribution and switchable circular dichroism (CD). The metasurface unit consists of a metallic inner ring embedded in vanadium oxide and a vanadium oxide outer ring, state switching by thermal control of vanadium oxide and a change in the frequency of the incident wave. Based on the switchable phase distribution, we designed a focusing vortex beam generator with adjustable focal lengths through simulation. Based on the switching CD capability, we simulate its mode switching in near-field imaging using numerical simulation, and innovatively propose an optical encryption method. Utilizing the chiral property, we also designed dual-channel switchable holographic imaging in the same frequency band, which combined with the state change of VO2 can realize a total of 4 holograms switching. Our proposed metasurface is expected to provide new ideas for the study of optical encryption and wavefront modulation of dynamics.

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

1. Introduction

Metasurfaces are man-made two-dimensional planar materials that can be tuned to control electromagnetic waves by adjusting their microstructure [13]. Chirality is the inability of an object to coincide with its mirror image [4]. However, with the development of chiral metasurfaces, symmetry-breaking structures can also cause different responses to different circularly polarized wave, giving rise to phenomena such as circular dichroism (CD) [5] and asymmetric transmission (AT) [6]. These structures are widely used in biomolecular detection, and CD spectroscopy [7,8]. Highly efficient chiral metasurfaces can strongly absorb one of the circularly polarized waves while transmitting (reflecting) the opposite wave, thus enabling independent modulation of LCP or RCP waves. Currently, the applications of chiral metasurfaces are also extended to wavefront manipulation, such as holographic imaging [9,10], vortex beam generation [1113], and beam splitting [14], and many interesting functions have been developed.

Devices with multifunctionality have great potential in future terahertz communication systems. However, most metasurfaces lack active tunability, which affects their functional expansion. Tunable materials such as graphene [15,16], vanadium dioxide [17,18] provide functional expansion of metasurfaces. Some studies [19,20] have designed graphene-based chiral metasurfaces that can dynamically adjust the amplitude of the metasurface, but only achieve on/off control of a single function. Utilizing the thermal phase transition properties of vanadium dioxide, works [21,22] have also achieved switchable wavefront manipulation between reflective and transmissive modes, but were only able to switch between transmissive and reflective for a single function. Niu et al [23] utilized a double-open ring structure combined with $\textrm{V}{\textrm{O}_\textrm{2}}$ to achieve a switchable phase distribution at the same frequency. Li et al [24] used a similar structure and designed a reconfigurable metasurface based on photosensitive silicon and $\textrm{V}{\textrm{O}_\textrm{2}}$, which can realize switchable phase distributions in three frequency bands, but it is difficult to practically applied because it requires corresponding stimulation of each unit. These two studies [23,24], despite realizing multifunctional switching, are limited to switching the phase distribution, and there is no way to broaden the scope of more functionalities, such as circular dichroism (CD).

In this paper, a switchable chiral metasurface based on $\textrm{V}{\textrm{O}_\textrm{2}}$, which can utilize the thermal phase transition properties of $\textrm{V}{\textrm{O}_\textrm{2}}$ and the electromagnetic response of the metasurface units at different frequencies to switch the already aligned phase distributions, and at the same time switch the circular dichroism (CD). Recently, the vortex phase distribution of the metasurfaces has been experimented with using simple metasources and used for photonic routing and topological excitation [25,26], achieving efficient excitation in the microwave band. However, there is still a lack of stable vortex beam sources in the terahertz band and a lack of tunability. We utilize the property that the proposed metasurface has a switchable arbitrary phase distribution to design a terahertz-band focused vortex beam generator, which can be adjusted to different focal lengths by temperature changes, providing a new idea for tunable vortex beams in the terahertz band. Subsequently, the switchable CD was used to develop a near-field imaging function, where the switching of two different images was achieved by thermal control of $\textrm{V}{\textrm{O}_\textrm{2}}$, which is a reversible process, and to simulate its application in optical encryption. Using the CD of the metasurface and its mirror structure, we realize polarization-multiplexed holography in the same frequency band, which, combined with the state switching of $\textrm{V}{\textrm{O}_\textrm{2}}$, achieves a total of four holograms switching, As shown in Fig. 1. This work is expected to be applicable to a new generation of devices such as optical encryption and message multiplexing communications.

 figure: Fig. 1.

Fig. 1. (a) Chiral metasurface structure. (b) Focused vortex beam switching by phase distribution at different temperatures. (c) Near-field image switching using CD at different temperatures. (d) Switchable Polarization Multiplexing Holography.

Download Full Size | PDF

2. Theoretical analysis and structural design

Using the classical Jones matrix method [27,28], it is possible to analyze the reflected waves generated by terahertz waves incident on the metasurface from a direction perpendicular to the metasurface. Because the terahertz waves excited by the transmitter are linearly polarized, the reflected electric field generated by a beam of linearly polarized electromagnetic waves irradiating the metasurface is represented as follows:

$$\left( {\begin{array}{{c}} {{E}_{x}^{r}}\\ {{E}_{y}^{r}} \end{array}} \right){=R}\left( {\begin{array}{{c}} {{E}_{x}^{i}}\\ {{E}_{y}^{i}} \end{array}} \right){=}\left( {\begin{array}{cc} {{{r}_{{xx}}}}&{{{r}_{{xy}}}}\\ {{{r}_{{yx}}}}&{{{r}_{{yx}}}} \end{array}} \right)\left( {\begin{array}{{c}} {{E}_{x}^{i}}\\ {{E}_{y}^{i}} \end{array}} \right)$$
where $\textrm{E}_\textrm{x}^\textrm{r}$, ${E}_{y}^{r}$ represents the electric field vector of the reflected waves, x and y represent the x-polarized and y-polarized waves, R denotes the reflectance transmission matrix of the incident wave polarized by the metasurface line, and ${E}_{x}^{i}$, ${E}_{y}^{i}$ represents the electric field vector of the incident waves.

In order to analyze circularly polarized waves, it is necessary to establish a relationship between circular polarization and line polarization. Here, a fundamental mode matrix U is introduced [29]:

$${U}= \displaystyle{{\rm 1} \over {\sqrt {\rm 2} }}\left( \begin{matrix} {\rm 1} & {\rm 1} \\ {\rm i} & {{\rm -i}} \end{matrix}\right)$$

According to the principle that two linearly polarized waves can form a circularly polarized wave, the transmission matrix R of the linearly polarized wave in Eq. (1) can be converted into the transmission matrix R’ of the circularly polarized wave using the fundamental mode matrix U:

$$\begin{array}{cc} {{R^{\prime}\;\ =\ \;\ }\left( {\begin{array}{{cc}} {{{r}_{\textrm{RR}}}}&{{{r}_{\textrm{RL}}}}\\ {{{r}_{\textrm{LR}}}}&{{{r}_{\textrm{LL}}}} \end{array}} \right)\textrm{=}{{U}^{\; {{-}1}}}{RU = }{{U}^{{{-}1}}}\left( {\begin{array}{{cc}} {{{r}_{{xx}}}}&{{{r}_{{xy}}}}\\ {{{r}_{{yx}}}}&{{{r}_{{yy}}}} \end{array}} \right){U}}\\ {{ =}\frac{{1}}{{2}}\left( {\begin{array}{{cc}} {{{r}_{{xx}}}{+}{{r}_{{yy}}}\textrm{+i(}{{r}_{{xy}}}\textrm{{-}}{{r}_{{yx}}}{)}}&{{{r}_{{xx}}}\textrm{{-}}{{r}_{{yy}}}{{-} \textrm i(}{{r}_{{xy}}}\textrm{+}{{r}_{{yx}}}{)}}\\ {{{r}_{{xx}}}{{-}}{{r}_{{yy}}}{+ \textrm i(}{{r}_{{xy}}}{ + }{{r}_{{yx}}}{)}}&{{{r}_{{xx}}}{+ }{{r}_{{yy}}}{{-}\textrm i(}{{r}_{{xy}}}{ {-}}{{r}_{{yx}}}{)}} \end{array}} \right)} \end{array}$$

Here, L and R denote the LCP and RCP waves, respectively, and ${{r}_{{ij}}}$ represents the transmission coefficient of the i-polarized reflected electric field component of the j-polarized incident electric field of amplitude 1, ${i,j} \in { [x,y,R,L]}$.

The difference in the reflection coefficients between LCP and RCP waves mainly stems from the spin-selective absorption of CP waves by the metasurface. Using the above reflection Jones matrix, the absorption of LCP and RCP waves by the designed metasurface can be expressed as follows:

$$\begin{aligned}{{{A}_{1}}{= 1{-} |}{{r}_{{RL}}}{{|}^{2}}{{-}|}{{r}_{{LL}}}{{|}^{2}}}\\ {{{A}_{2}}{= 1{-}|}{{r}_{{LR}}}{{|}^{2}}{{-}|}{{r}_{{RR}}}{{|}^{2}}} \end{aligned}$$

${{A}_{1}}$ in Eq. (4) denotes the absorption of the metasurface for LCP waves, while ${{A}_{2}}$ denotes the absorption of the metasurface for RCP waves. Circular dichroism (CD) can be characterized by the differential absorption of the metasurface for LCP and RCP waves:

$${CD=}{{A}_{2}}{{-}}{{A}_{1}}$$

In this work, the designed metasurface structural unit is shown in Fig. 2(a) and consists of three layers of structure. The bottom layer is a gold substrate with a thickness of 2 µm. The intermediate layer is a polyimide dielectric material with a thickness of 30µm and a dielectric constant of 3.5 + 0.00945i [30]. The top layer is made of two metal split rings, the outer ring is made of $\textrm{V}{\textrm{O}_\textrm{2}}$ material with a thickness of 1.5 µm, and the inner ring is made of metal and embedded $\textrm{V}{\textrm{O}_\textrm{2}}$ with a thickness of 1 µm. The other parameters are P = 115 µm, R1 = 32.5 µm, R2 = 52.5 µm, W1 = 5 µm, W2 = 12.5 µm, θ1 = 40°, θ2 = 4°. As the temperature changes from 300 K to the phase transition temperature of 350 K, the change in $\textrm{V}{\textrm{O}_\textrm{2}}$ dielectric constant can be similarly described by the Drude model [31]:

$$\mathrm{\varepsilon (\omega )\;\ =\ \;\ }{\mathrm{\varepsilon }_\infty }\textrm{{-}}\frac{{{\mathrm{\omega }_{p}}^{2}\mathrm{(\sigma )}}}{{\textrm{(}{\mathrm{\omega }^\textrm{2}}{\ +\ i\gamma \omega )}}}$$
Where the damping frequency ${\gamma \;\ =\ \;\ 5}\mathrm{.73\ \times 1}{\textrm{0}^{\textrm{13}}}\textrm{rad/s}$, the dielectric constant at infinite frequency ${{\varepsilon }_\infty }\textrm{=12}$ the plasma frequency of $\textrm{V}{\textrm{O}_\textrm{2}}$ can be written as ${{\omega}_{p}}^\textrm{2}{(\sigma)\;\ =\ \;\ (\sigma /}{\mathrm{\sigma }_\textrm{0}}\textrm{)}{\mathrm{\omega }_{p}}^\textrm{2}\textrm{(}{\mathrm{\sigma }_\textrm{0}}\textrm{)}$, where ${{\sigma }_{\textrm{0}}}\mathrm{\ =\ \;\ 3\ \times 1}{\textrm{0}^\textrm{5}}\textrm{S/m}$ and ${{\omega }_{p}}({{{\sigma }_{0}}} )\textrm{= 1}\mathrm{.4\ \times 1}{\textrm{0}^{\textrm{15}}}\textrm{rad/s}$. Thus, the conductivities of $\textrm{V}{\textrm{O}_\textrm{2}}$ for different phase transitions $\mathrm{\sigma \;\ =\ \;\ 200\;\ S/m}$ and $\mathrm{\sigma \;\ =\ \;\ 2\ \times 1}{\textrm{0}^\textrm{5}}\textrm{S/m}$, corresponding to the insulating and metallic states [32].

 figure: Fig. 2.

Fig. 2. (a) Schematic of the designed metasurface structured unit. (b) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: metasurface unit. (c) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: co-polarized reflection coefficient. (d) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: metasurface unit. (e) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: co-polarized reflection coefficient.

Download Full Size | PDF

The proposed metasurfaces are also easy to prepare in practice by utilizing the photolithography process. The general preparation steps are as follows [33]. First, a layer of polyimide was deposited on the metal, and then a VO2 film and photoresist were successively deposited, followed by patterning of the VO 2 film using photolithography with a mask. After a series of standard photolithography processes such as exposure and development, the VO 2 pattern was left on the polyimide layer. The same technique can be used for the preparation of metal inner rings.

Calculation of the electromagnetic response of metasurfaces using the finite difference in time domain (FDTD) approach. Figures 3(a) and (b) show the reflection amplitude spectra of the metasurface structure at different temperatures respectively ${{r}_{{ij}}}$. It is obvious from the figure that at different temperatures, $\textrm{V}{\textrm{O}_\textrm{2}}$ has different frequency band responses when it is in different states. When $\textrm{V}{\textrm{O}_\textrm{2}}$ is metallic, the inner ring is equivalent to being shorted and only the outer ring of $\textrm{V}{\textrm{O}_\textrm{2}}$ is operating, as shown in Fig. 3(c), which exhibits a high polarization selectivity at 0.68 THz, ${CD=}{{A}_{2}}\textrm{{-}}{\textrm{A}_\textrm{1}}\textrm{=0}\textrm{.6}$. In contrast, when ${CD=}{{A}_{2}}\textrm{-}{{A}_\textrm{1}}\textrm{=0}\textrm{.6}$ is in the insulated state, as shown in Fig. 3(d), only the metal-open inner ring is in operation and exhibits great polarization selectivity at 0.87 THz, $\textrm{CD=}{\textrm{A}_\textrm{2}}\textrm{{-}}{\textrm{A}_\textrm{1}}\textrm{= 0}\textrm{.62}$. This means that the response of circularly converted dichroism in two different frequency bands can be switched by the thermal phase transition of $\textrm{V}{\textrm{O}_\textrm{2}}$, and temperature control is the easiest to achieve compared to electric, optical, and magnetic modulation.

 figure: Fig. 3.

Fig. 3. (a) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: amplitude spectra. (b) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: amplitude spectra. (c) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Absorption spectra. (d) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: absorption spectra.

Download Full Size | PDF

The Pancharatnam-Berry (PB) theory can often be used to manipulate the phase of CP waves [34,35]. When the metasurface structure is rotated around the z-axis by a certain angle θ, the incident CP wave is co-polarized by the action of the metasurface, and the reflected wave carries an additional phase delay of ±2θ. In this section, the phase gradient of the chiral metasurface structure is simulated when VO2 is in different states, as shown in Fig. 4(a). When $\textrm{V}{\textrm{O}_\textrm{2}}$ is in insulated state, the inner ring rotates every 30° from 0° to 180°. From the figure, we can see that the phase changes linearly, covering 2π, and the amplitudes are all kept constant at 0.8. When $\textrm{V}{\textrm{O}_\textrm{2}}$ is in the metallic state, as shown in Fig. 4(b), it also satisfies the phase 2π coverage, while the amplitude remains constant at 0.78. We selected six units with rotation angles (0°, 30°, 60°, 90°, 120°, 150°) for subsequent phase alignment. We also simulated the change in CD values for different rotations of the inner and outer rings, see Fig. S1 in the Supplement 1, the CD values do not change much as the angle changes.

 figure: Fig. 4.

Fig. 4. (a) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: Phase gradient formed by rotating open metal ring. (b) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Phase gradient formed by rotating $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring.

Download Full Size | PDF

Here, we analyze the generation of switchable CD and switchable phase distributions. As shown in Fig. 5(a) and (b), it can be seen that when the embedded $\textrm{V}{\textrm{O}_\textrm{2}}$ is in the metallic state, the opening of the inner metallic ring is shorted to form a complete circle, thus, the inner ring has almost no electric field distribution. Only the outer ring is operating, with a strong electric field response at 0.68 THz. As shown in Fig. 5(b) and (c), when the embedded $\textrm{V}{\textrm{O}_\textrm{2}}$ is in the insulated state, the reflectance factor is close to 0.9 at 0.87 THz, which indicates that the metal inner ring is not shorted. This further illustrates that by exploiting the thermal phase change properties of $\textrm{V}{\textrm{O}_\textrm{2}}$, switchable CD and switchable phase distributions can be achieved.

 figure: Fig. 5.

Fig. 5. (a) Variation of reflection coefficient of cp wave passing through metasurface when VO2 is in metallic and insulating states, respectively. (b) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Electric field distribution on the metasurface. (c) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: Electric field distribution on the metasurface.

Download Full Size | PDF

3. Results and discussions

Our proposed metasurface unit can switch the already aligned phase distributions when $\textrm{V}{\textrm{O}_\textrm{2}}$ is in different states, and also switch the CD. To verify its potential, we design and simulate a focused vortex beam generator with switchable focal length function and a near-field imaging device with switchable imaging modes, and innovatively propose an optical encryption method. In Fig. s2 of the Supplement 1, we also implemented switchable vortex beam and near-field imaging. In order to enhance the potential applications of this metasurface, we have also designed dual-channel switchable holographic imaging in the same frequency band.

3.1 Focused vortex beam generator with switchable focal lengths

To obtain a vortex beam carrying OAM, when lining up the metasurface units, the units at different positions must satisfy the corresponding phase distribution in Eq. (7). Where l is the number of topological charges generating the OAM and (x, y) is the given coordinates in space:

$${\mathrm{\varphi }_\textrm{1}}{(x,y)\; = \; l\; } \cdot \textrm{arctan(}\frac{\textrm{y}}{\textrm{x}}\textrm{)}$$

The need to control the focusing of the vortex beam to a certain point, which facilitates the reception of the communication system. At this point, the units at different positions must satisfy the corresponding focusing phase distribution in Eq. (8). where f is the focal length and λ is the wavelength.

$${\mathrm{\varphi }_\textrm{2}}{(x,\;\ y)\;\ =\ \;\ (2\pi } \cdot \textrm{(}\sqrt {{\textrm{f}^{\textrm{2}}}\textrm{+(}{\textrm{x}^\textrm{2}}\textrm{+}{\textrm{y}^\textrm{2}}\textrm{)}} \textrm{{-}}\mathrm{\mid }{f}\mathrm{\mid }\mathrm{))\;\ /\;\ \lambda }$$

To focus the desired vortex beam to a certain point, the phase variation of the unit structure can superimpose the phase distribution of Eq. (7) and (8). According to the PB phase principle, the metasurface cell rotates at an angle of 1/2 of the corresponding phase. Therefore, the phase distribution of the unit structure at the (x, y) position can be converted to the corresponding rotation angle, satisfying Eq. (9) [36].

$$\begin{array}{{c}} {\theta (x,y) =\frac{\textrm{1}}{\textrm{2}}{\mathrm{\varphi }_\textrm{1}}{(x, y) + }\frac{\textrm{1}}{\textrm{2}}{\mathrm{\varphi }_\textrm{2}}{(x,y)}}\\ {\textrm{ = }\frac{\textrm{1}}{\textrm{2}}{l} \cdot \textrm{arctan(}\frac{{y}}{{x}}{)+}\frac{{1}}{\textrm{2}}\mathrm{(2\pi } \cdot \textrm{(}\sqrt {{{f}^{{2}}}\textrm{ + (}{{x}^{2}}\textrm{+}{\textrm{y}^\textrm{2}}\textrm{)}} {{-}\ | f\ |))\;\ /\;\ \lambda }} \end{array}$$

Based on the above principal analysis of phase distribution, we used 25 × 25 metasurface units for two switchable phase arrangements. Figures 6(a) and (b) show the phase distribution of the vortex beam and the phase distribution of the focused beam, respectively, and Fig. 6(c) represents the phase distribution of the focused vortex beam after the superposition of the vortex beam and the focused beam. To realize the switchable focal lengths of the focused vortex beams, the phase arrangement of the focused vortex beams at different focal lengths for each of the inner and outer rings of the metasurface was determined, as shown in Fig. 6(d) and (e). When the temperature is 300 K and the incident wave frequency is 0.87 THz, $\textrm{V}{\textrm{O}_\textrm{2}}$ is in an insulated state. Set the focal length to 1600 µm and the topological charge number to 1. When the temperature is 350 K and the incident wave frequency is 0.68 THz, $\textrm{V}{\textrm{O}_\textrm{2}}$ is in a metallic state. Set the focal length to 900 µm and the topological charge to 1.

 figure: Fig. 6.

Fig. 6. (a) Phase distribution of the vortex beam. (b) Phase distribution of the focused beam. (c) Phase distribution of the superposition of the vortex and focused beams. (d) Schematic of the metasurface structure. (e) Enlargement of the localized structure.

Download Full Size | PDF

Figure 7 shows the simulation results when $\textrm{V}{\textrm{O}_\textrm{2}}$ is insulated. As shown in Fig. 7(a), it can be seen that the beam has a good focusing effect. The energy amplitude of the beam has a hollow circle in the center, which is the zero of the electric field. This result proves that a focused vortex beam is produced. As shown in Fig. 7(b), it can be seen that the electric field is strongest at a distance of z = 1530 µm, indicating that this is the focal point, close to the preset focal length value. Figure 7(c) shows the phase distribution of the generated vortex beam. Mode purity is often used to characterize the quality of a vortex beam. Mode purity can be obtained by comparing the power of the main mode with that of all modes. As shown in Fig. 7(d), the mode purity at this point is 0.84. Figure 7(e) shows the distribution of the electric field along the X-axis corresponding to the focal point at a distance of z = 1530 µm. It can be seen that when the points x = 203 µm and x = -193 µm approach 1, the electric field is the maximum value. The full width half maximum (FWHM) of the focal spot center is 613 µm.

 figure: Fig. 7.

Fig. 7. $\textrm{V}{\textrm{O}_\textrm{2}}$ in the insulating state:(a) The electric field in the X-Y plane. (b) The electric field in the X-Z plane. (c) Phase of the X-Y plane. (d) Mode purity. (e) The electric field in the X-axis (z = 1530 µm).

Download Full Size | PDF

Figure 8 shows the simulation results when $\textrm{V}{\textrm{O}_\textrm{2}}$ is in a metallic state and the $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring works independently. Similar to Fig. 7(a), Fig. 8(a) also demonstrates the generation of a focused vortex beam. As shown in Fig. 8(b), the maximum electric field intensity is observed at the focal point z = 820 µm, which is close to the present value. Figure 8(c) shows the phase distribution corresponding to the vortex beam. As shown in Fig. 8(d), the mode purity is 0.88 at this point. Figure 8(e) shows the distribution of the electric field on the X-axis corresponding to the focal point at a distance of z = 820 µm. It can be seen that the electric field is maximized when the points x = 186 µm and x = -183 µm approach 1. The full width half maximum (FWHM) of the focal spot center is 605 µm.

 figure: Fig. 8.

Fig. 8. $\textrm{V}{\textrm{O}_\textrm{2}}$ in the metallic state: (a) The electric field in the X-Y plane. (b) The electric field in the X-Z plane. (c) Phase of the X-Y plane. (d) Mode purity. (e) The electric field in the X-axis (z = 820 µm).

Download Full Size | PDF

When the $\textrm{V}{\textrm{O}_\textrm{2}}$ is in different states, the designed focused vortex generator has excellent performance at 0.68 and 0.87 THz. In addition, the designed metasurface can also realize more functional switching combinations in addition to focusing, such as beam splitting and topological charge number switching, with great regulatory freedom.

3.2 Dual-band switchable near-field imaging

The designed chiral metasurface has high-efficiency CD in two different frequency bands, which can be used for near-field imaging applications without considering the PB phase. According to the previous analysis and verification, the designed unit array composed of metasurface units has great potential in imaging display. Because the structure is chiral, the reverse CD can be obtained by mirroring. As shown in Fig. 9(a) and (b), arrange the unit structure of forward CD and reverse CD as 25 × 25 units. Among them, the $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring and its mirror image structure are arranged into the letter “H”, and the metal inner ring and its mirror image structure are arranged into the letter “N”, and then the pattern is displayed by detecting the reflection fields of the LCP and RCP.

 figure: Fig. 9.

Fig. 9. (a) The $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring and its mirror image structure are arranged into the letter “H”. (b) The metal inner ring and its mirror image structure are arranged into the letter “N”. (c) Switching of “H” and “N” at RCP wave incidence. (d) Switching of “H” and “N” at LCP wave incidence.

Download Full Size | PDF

As shown in Fig. 9(c) and (d), a detector is placed 200µm away from the array element to observe the image formed by the co-polarized reflection waves after the incident waves are RCP and LCP. The images of the letter “H” and the letter “N” can be displayed, and the imaging effect is better. Switching of the two-letter images is achieved by changing the temperature. When the temperature is 300 K and $\textrm{V}{\textrm{O}_\textrm{2}}$ is in the insulated state, the metasurface responds strongly at 0.87 THz, at which time the letter “N” arranged by the metal inner ring is imaged. This means that we encode the original image and the encrypted image to the metasurface according to a specific amplitude configuration. At room temperature and a specific incident wave frequency, the reflected wave field is the original image. Only by changing both the temperature and the frequency of the incident wave can the encrypted image appear, and the entire switching process is reversible.

We also innovatively propose an optical encryption method. As shown in Fig. 10(a), the entire structure consisting of a 30 × 10 metasurface units. The entire array is divided into 12 sub-regions, and the same 5 × 5 units is adopted in each sub-region. The initial and mirror structure of the inner and outer rings are arranged into each sub-region according to the design requirements. The reflected strong electric field is marked as “1”, and the weak electric field is marked as “0”. As can be seen from the results in Fig. 10(c), when the incident frequency of the LCP wave is 0.87 THz, when the temperature is 300 K, and $\textrm{V}{\textrm{O}_\textrm{2}}$ is in an insulated state, the near-field number is displayed as “101010/010101”. When the incident frequency of LCP wave is 0.68 THz, the temperature is 350 K, and $\textrm{V}{\textrm{O}_\textrm{2}}$ is in the metallic state, the near-field number is displayed as “110101/010010”. According to this result, we can express the information that needs to be transmitted by binary coding, in which the metal inner ring is responsible for storing the binary code of the original information, and the $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring is responsible for storing the binary code of the encrypted information. As shown in Fig. 10(c), the reflected electric field of the inner metal ring is the valid information of “101010/010101”, while the reflected electric field of the outer ring of $\textrm{V}{\textrm{O}_\textrm{2}}$ is the invalid information after encryption of “110101/010010”. The original information and the encrypted information can be set independently and do not interfere with each other. The encryption/decryption process can be realized by changing the temperature and frequency of the incident CP wave. It is believed that the design of the unit structure combined with better processing technology can provide higher security in the encryption/decryption and transmission of information.

 figure: Fig. 10.

Fig. 10. (a) Arrays are made using a 30 × 10 metasurface unit structure and are divided into 12 subregions. (b) Structure of the metasurface unit corresponding to each subregion. (c) Reflected electric field intensity distribution at different temperatures under incident LCP wave. (d) Reflected electric field intensity distribution at different temperatures under incident RCP wave.

Download Full Size | PDF

3.3 Dual-channel switchable holographic imaging at the same frequency

From Fig. 9(ab) in Subsection 3.2, it can be seen that the initial and mirror structures of the metasurface unit have opposite electromagnetic responses for the incidence of LCP waves. This is a property of chiral structures, where the initial structure responds only to the LCP and the mirror structure responds only to the RCP. Using this feature, we can realize the holographic imaging function of dual-channel switching under the same frequency. From the previous section, it can be seen that the inner and outer rings of the metasurface have independent electromagnetic responses to two different frequency bands in the two states of $\textrm{V}{\textrm{O}_\textrm{2}}$, which indicates that the same-frequency dual-channel switchable holographic imaging can be established in each state, and that a total of four holographic image switching can be achieved. This also requires the phase gradient of the mirror structure to be established as well, see Figure S3 in the Supplement 1.

We constructed holograms of the letters “H”, “T”, “L” and “A”. We would like to get “H” and “T” for LCP and RCP channels respectively at 0.87 THz, and “L “ and “A”. The phase distribution of metasurface holographic imaging is generally obtained by an algorithm for phase retrieval, and this design utilizes the classical G-S algorithm, which uses the Fresnel-Kirchhoff diffraction formula:

$$\begin{matrix}{{U(x,y)\; = }\displaystyle{{\rm 1} \over {{j\lambda }}}\mathop \int\nolimits_{{\rm -}\infty }^\infty \mathop \int\nolimits_{{\rm -}\infty }^\infty {U(}{x}_{0}{\rm ,}{y}_{0}{\rm )}\displaystyle{{{exp(jkr)}} \over {r}}{K(\theta )d}{x}_{\rm 0}{d}{y}_{\rm 0}} \end{matrix} $$

We use the partitioned G-S algorithm to obtain the desired phase distribution of the electric field on the metasurface so that the cell distributions of the LCP and RCP channels are arranged alternately. The flow of the partitioned G-S algorithm is shown in Fig. 11, where the initial random phase and the partitioned amplitude distribution M1 (M2) are input. After Fourier transform, the amplitude distribution of the computed image is evaluated. If the image quality is poor, the amplitude of the target image is combined with the computed phase distribution to perform an inverse Fourier transform to obtain a new input. Iterate repeatedly until a more like phase distribution of the electric field is obtained. The phase distributions of the respective cells of the initial and mirror structures were finally obtained.

 figure: Fig. 11.

Fig. 11. Schematic flow of partitioned iterative GS algorithm with two channels “LCP” and “RCP”.

Download Full Size | PDF

The simulation results are shown in Fig. 12, and it can be seen that the crosstalk between the two channels is small. As can be seen from Fig. 11(bc), when the temperature is increased and VO2 is in the metallic state, the outer ring responds, and the initial and specular structures of the outer ring are arranged alternately, which allows switching between “H” and “T” holographic imaging at a frequency of 0.68 THz by simply changing the incident circular polarization state. From Fig. 11, it can be seen that when the temperature decreases, the inner loop response, at 0.87 THz, is able to accomplish the switching between “L” and “A”. When the frequency and temperature are changed simultaneously, the holographic imaging map formed by the outer ring and the holographic imaging map formed by the inner ring can also be switched. This greatly improves the switching freedom of this metasurface.

 figure: Fig. 12.

Fig. 12. (a) Schematic diagram of the metasurface. (b) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Outer ring independent work. (c) Letter switching between “H” and “L” at the same frequency. (d) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: Inner ring independent work. (e) Letter switching between “H” and “L” at the same frequency.

Download Full Size | PDF

4. Conclusion

In conclusion, we designed a terahertz chiral metasurface with a switchable phase distribution and switchable circular dichroism (CD). By simulation, the inner and outer rings of this structure have independent phase distributions (covering 2π) and circular dichroism (CD) at frequencies of 0.68 and 0.87 THz. Temperature control of the phase change of $\textrm{V}{\textrm{O}_\textrm{2}}$ and changing the frequency of the incident wave enable different phase arrangements and circular dichroism (CD) switching. Utilizing this property, we design focus-switchable vortex beam generators and near-field imaging switching, and present their application in optical encryption. Using the partitioned phase design, we realized switchable holographic imaging in the same frequency band, which combined with the state change of VO2 can realize the switching of four holograms. The proposed metasurface structure is expected to provide a noteworthy perspective for research in the fields of optical encryption and dynamic wavefront modulation.

Funding

National Key Research and Development Program of China (2022YFA1203500); National Natural Science Foundation of China (12274424, 61988102, 62064004, 62174046, 62274048); Open Fund for Innovation and Entrepreneurship of College Students (RSXH20231165808X, RSXH20231165817X).

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.

Supplemental document

See Supplement 1 for supporting content.

References

1. L. Kang, S. Rodrigues, M. Taghinejad, et al., “Preserving Spin States upon Reflection: Linear and Nonlinear Responses of a Chiral Meta-Mirror,” Nano Lett. 17(11), 7102–7109 (2017). [CrossRef]  

2. Q. He, S. Sun, S. Xiao, et al., “High-efficiency metasurfaces: Principles, realizations, and applications,” Adv. Opt. Mater. 6(19), 1800415 (2018). [CrossRef]  

3. G. Rubio, J. Mosquera, V. Kumar, et al., “Micelle-directed chiral seeded growth on anisotropic gold nanocrystals,” Science 368(6498), 1472–1477 (2020). [CrossRef]  

4. J. Hwang, M. Song, B. Park, et al., “Electro-tunable optical diode based on photonic bandgap liquid-crystal heterojunctions,” Nat. Mater. 4(5), 383–387 (2005). [CrossRef]  

5. T. Cao, Y. Li, X. Zhang, et al., “Theoretical study of tunable chirality from graphene integrated achiral metasurfaces,” Photonics Res. 5(5), 441–449 (2017). [CrossRef]  

6. N. Parappurath, F. Alpeggiani, L. Kuipers, et al., “The origin and limit of asymmetric transmission in chiral resonators,” ACS Photonics 4(4), 884–890 (2017). [CrossRef]  

7. E. Philip, M. Güngördü, S. Pal, et al., “Review on Polarization Selective Terahertz Metamaterials: from Chiral Metamaterials to Stereometamaterials,” J. Infrared, Millimeter, Terahertz Waves 38(9), 1047–1066 (2017). [CrossRef]  

8. S. Park, S. Cha, G. Shin, et al., “Sensing viruses using terahertz nano-gap metamaterials,” Biomed. Opt. Express 8(8), 3551–3558 (2017). [CrossRef]  

9. Q. Wang, E. Plum, Q. Yang, et al., “Reflective chiral meta-holography: multiplexing holograms for circularly polarized waves,” Light: Sci. Appl. 7(1), 25–227 (2018). [CrossRef]  

10. W. Pan, Z. Wang, Y. Chen, et al., “High-efficiency generation of far-field spin-polarized wavefronts via designer surface wave metasurfaces,” Nanophotonics 11(9), 2025–2036 (2022). [CrossRef]  

11. J. Li, Y. Zhang, J. Li, et al., “Amplitude modulation of anomalously reflected terahertz beams using all-optical active Pancharatnam–Berry coding metasurfaces,” Nanoscale 11(12), 5746–5753 (2019). [CrossRef]  

12. H. Xu, G. Hu, Y. Li, et al., “Interference-assisted kaleidoscopic meta-plexer for arbitrary spin-wavefront manipulation,” Light: Sci. Appl. 8(1), 3–62 (2019). [CrossRef]  

13. Z. Zhou and Z. Song, “Terahertz mode switching of spin reflection and vortex beams based on graphene metasurfaces,” Opt. Laser Technol. 153, 108278 (2022). [CrossRef]  

14. J. Li and C. Zhou, “Multifunctional reflective dielectric metasurface in the terahertz region,” Opt. Express 28(15), 22679–22689 (2020). [CrossRef]  

15. X. Tao, L. Qi, T. Fu, et al., “A tunable dual-band asymmetric transmission metasurface with strong circular dichroism in the terahertz communication band,” Opt. Laser Technol. 150, 107932 (2022). [CrossRef]  

16. H. Rajabalipanah, K. Rouhi, A. Abdolali, et al., “Real-time terahertz meta-cryptography using polarization-multiplexed graphene-based computer-generated holograms,” Nanophotonics 9(9), 2861–2877 (2020). [CrossRef]  

17. T. Lv, Y. Li, C. Qin, et al., “Versatile polarization manipulation in vanadium dioxide-integrated terahertz metamaterial,” Opt. Express 30(4), 5439–5449 (2022). [CrossRef]  

18. F. Baranzadeh and N. Nozhat, “High Performance Plasmonic Nano-Biosensor Based on Tunable Ultra-Narrowband Perfect Absorber Utilizing Liquid Crystal,” Plasmonics 16(1), 253–262 (2021). [CrossRef]  

19. J. Li, J. Li, C. Zheng, et al., “Dynamic control of reflective chiral terahertz metasurface with a new application developing in full grayscale near field imaging,” Carbon 172, 189–199 (2021). [CrossRef]  

20. H. Liu, S. Duan, C. Chen, et al., “Graphene-enabled chiral metasurface for terahertz wavefront manipulation and multiplexing holographic imaging,” Opt. Mater. 147, 114654 (2024). [CrossRef]  

21. X. Qiao, A. Wang, W. Xiao, et al., “A simple structured and homo-frequency point mode switchable THz full-space metasurface based on temperature-controlled vanadium dioxide,” Opt. Laser Technol. 170, 110233 (2024). [CrossRef]  

22. L. Yang and J. Li, “Terahertz vortex beam generator carrying orbital angular momentum in both transmission and reflection spaces,” Opt. Express 30(20), 36960–36972 (2022). [CrossRef]  

23. J. Niu, C. Li, Q. Yao, et al., “Switchable metasurface with double phase distributions for wavefront manipulation,” Phys. Scr. 98(8), 085517 (2023). [CrossRef]  

24. J. Li, J. Ren, R. Xiong, et al., “Reconfigurable three multi-mode terahertz metasurface,” Opt. Express 31(22), 35583–35593 (2023). [CrossRef]  

25. Z. Guo, Y. Long, H. Chen, et al., “Anomalous unidirectional excitation of high-k hyperbolic modes using all-electric metasources,” Adv. Photon. 3(03), 036001 (2021). [CrossRef]  

26. X. Wang, Z. Guo, H. Chen, et al., “Unique Huygens-Fresnel electromagnetic transportation of chiral Dirac wavelet in topological photonic crystal,” Nat. Commun. 14(1), 3040 (2023). [CrossRef]  

27. J. Li, J. Li, C. Zheng, et al., “Lossless dielectric metasurface with giant intrinsic chirality for terahertz wave,” Opt. Express 29(18), 28329–28337 (2021). [CrossRef]  

28. C. Menzel, C. Rockstuhl, F. Lederer, et al., “Advanced Jones calculus for the classification of periodic metamaterials,” Phys. Rev. A 82(5), 053811 (2010). [CrossRef]  

29. Z. Wang, H. Jia, K. Yao, et al., “Circular Dichroism Metamirrors with Near-Perfect Extinction,” ACS Photonics 3(11), 2096–2101 (2016). [CrossRef]  

30. F. Hu, L. Wang, B. Quan, et al., “Design of a polarization insensitive multiband terahertz metamaterial absorber,” J. Phys. D: Appl. Phys. 46(19), 195103 (2013). [CrossRef]  

31. Q. Wen, H. Zhang, Q. Yang, et al., “Terahertz metamaterials with VO2 cut-wires for thermal tunability,” Appl. Phys. Lett. 97(2), 021111 (2010). [CrossRef]  

32. Z. Zhang, Q. Xie, L. Guo, et al., “Dual-controlled tunable dual-band and ultra-broadband coherent perfect absorber in the THz range,” Opt. Express 30(17), 30832–30844 (2022). [CrossRef]  

33. S. Zhu, S. Bo, G. Guo, et al., “Terahertz Metasurfaces for Thermally Controlled Optical Encryption,” Laser Photonics Rev. 17(10), 2300233 (2023). [CrossRef]  

34. H. Xu, G. Wang, T. Cai, et al., “Tunable Pancharatnam-Berry metasurface for dynamical and high-efficiency anomalous reflection,” Opt. Express 24(24), 27836–27848 (2016). [CrossRef]  

35. R. Mao, G. Wang, T. Cai, et al., “Ultra-thin and high-efficiency full-space Pancharatnam-Berry metasurface,” Opt. Express 28(21), 31216–31225 (2020). [CrossRef]  

36. R. Nie, C. He, R. Zhang, et al., “Vanadium dioxide-based terahertz metasurfaces for manipulating wavefronts with switchable polarization,” Opt. Laser Technol. 159, 109010 (2023). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Phase distribution and circular dichroism switchable terahertz chiral metasurface:Supplemental Document

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) Chiral metasurface structure. (b) Focused vortex beam switching by phase distribution at different temperatures. (c) Near-field image switching using CD at different temperatures. (d) Switchable Polarization Multiplexing Holography.
Fig. 2.
Fig. 2. (a) Schematic of the designed metasurface structured unit. (b) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: metasurface unit. (c) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: co-polarized reflection coefficient. (d) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: metasurface unit. (e) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: co-polarized reflection coefficient.
Fig. 3.
Fig. 3. (a) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: amplitude spectra. (b) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: amplitude spectra. (c) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Absorption spectra. (d) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: absorption spectra.
Fig. 4.
Fig. 4. (a) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: Phase gradient formed by rotating open metal ring. (b) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Phase gradient formed by rotating $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring.
Fig. 5.
Fig. 5. (a) Variation of reflection coefficient of cp wave passing through metasurface when VO2 is in metallic and insulating states, respectively. (b) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Electric field distribution on the metasurface. (c) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: Electric field distribution on the metasurface.
Fig. 6.
Fig. 6. (a) Phase distribution of the vortex beam. (b) Phase distribution of the focused beam. (c) Phase distribution of the superposition of the vortex and focused beams. (d) Schematic of the metasurface structure. (e) Enlargement of the localized structure.
Fig. 7.
Fig. 7. $\textrm{V}{\textrm{O}_\textrm{2}}$ in the insulating state:(a) The electric field in the X-Y plane. (b) The electric field in the X-Z plane. (c) Phase of the X-Y plane. (d) Mode purity. (e) The electric field in the X-axis (z = 1530 µm).
Fig. 8.
Fig. 8. $\textrm{V}{\textrm{O}_\textrm{2}}$ in the metallic state: (a) The electric field in the X-Y plane. (b) The electric field in the X-Z plane. (c) Phase of the X-Y plane. (d) Mode purity. (e) The electric field in the X-axis (z = 820 µm).
Fig. 9.
Fig. 9. (a) The $\textrm{V}{\textrm{O}_\textrm{2}}$ outer ring and its mirror image structure are arranged into the letter “H”. (b) The metal inner ring and its mirror image structure are arranged into the letter “N”. (c) Switching of “H” and “N” at RCP wave incidence. (d) Switching of “H” and “N” at LCP wave incidence.
Fig. 10.
Fig. 10. (a) Arrays are made using a 30 × 10 metasurface unit structure and are divided into 12 subregions. (b) Structure of the metasurface unit corresponding to each subregion. (c) Reflected electric field intensity distribution at different temperatures under incident LCP wave. (d) Reflected electric field intensity distribution at different temperatures under incident RCP wave.
Fig. 11.
Fig. 11. Schematic flow of partitioned iterative GS algorithm with two channels “LCP” and “RCP”.
Fig. 12.
Fig. 12. (a) Schematic diagram of the metasurface. (b) Metallic $\textrm{V}{\textrm{O}_\textrm{2}}$: Outer ring independent work. (c) Letter switching between “H” and “L” at the same frequency. (d) Insulating $\textrm{V}{\textrm{O}_\textrm{2}}$: Inner ring independent work. (e) Letter switching between “H” and “L” at the same frequency.

Equations (10)

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

( E x r E y r ) = R ( E x i E y i ) = ( r x x r x y r y x r y x ) ( E x i E y i )
U = 1 2 ( 1 1 i i )
R   =     ( r RR r RL r LR r LL ) = U 1 R U = U 1 ( r x x r x y r y x r y y ) U = 1 2 ( r x x + r y y +i( r x y {-} r y x ) r x x {-} r y y i ( r x y + r y x ) r x x r y y + i ( r x y + r y x ) r x x + r y y i ( r x y r y x ) )
A 1 = 1 | r R L | 2 | r L L | 2 A 2 = 1 | r L R | 2 | r R R | 2
C D = A 2 A 1
ε ( ω )   =     ε {-} ω p 2 ( σ ) ( ω 2   +   i γ ω )
φ 1 ( x , y ) = l arctan( y x )
φ 2 ( x ,   y )   =     ( 2 π ( f 2 +( x 2 + y 2 ) {-} f ) )   /   λ
θ ( x , y ) = 1 2 φ 1 ( x , y ) + 1 2 φ 2 ( x , y )  =  1 2 l arctan( y x ) + 1 2 ( 2 π ( f 2  + ( x 2 + y 2 )   | f   | ) )   /   λ
U ( x , y ) = 1 j λ U ( x 0 , y 0 ) e x p ( j k r ) r K ( θ ) d x 0 d y 0
Select as filters


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