Dual-Channel Mid-Infrared Toroidal Metasurfaces for Wavefront Modulation and Imaging Applications

In this paper, we propose a dual-channel mid-infrared toroidal metasurface that consists of split equilateral triangular rings. The electromagnetic responses are analyzed by the finite-difference-time-domain (FDTD) method and temporal coupled-mode theory (TCMT). The results show that one channel of the metasurface is insensitive to the polarization angle of the incident light and temperature, while the other channel is sensitive. The reflectance and resonance wavelength can be manipulated by the polarization angle and temperature independently. Based on such a mechanism, we propose metasurfaces for two-bit programmable imaging and thermal imaging. The metasurfaces are believed to have potential applications in information processing and thermal radiation manipulation.


Introduction
Composed of subwavelength antennas, metasurfaces are widely used to manipulate the polarization, phase, and amplitude of electromagnetic waves because metasurfaces have many advantages such as the small size and little loss [1]. Lots of devices based on metasurfaces have been proposed, e.g., absorbers [2], invisibility cloaks [3], special beam generators [4], programmable elements [5], and thermal imaging devices [6]. The design of programmable metasurfaces is an effective approach to suppressing the crosstalk [7][8][9][10][11], and thermal imaging metasurfaces can be used for temperature measurement [12], nightvision [13], and heat detection [14]. For example, Shang et al. designed a reconfigurable thermal metasurface to achieve arbitrary objects' thermograms [15]. Alipour et al. proposed a bolometer sensor with high resistance and sensitivity [16]. Liu et al. illustrated a metalliquid-crystal-metal metasurface to realize dynamic thermal camouflage [17]. With these advantages, metasurfaces are believed to have great potential applications.
Metasurfaces are mainly divided into sensitive and insensitive types. Sensitive metasurfaces can work as sensors [18,19], bolometers, and photodetectors. The tunability can be realized by introducing external excitation methods including electric [20], optical [21][22][23], mechanical controlling [24], and thermal tuning [25,26]. As for the insensitive type, they usually maintain a relatively stable performance to environmental changes. Recently, in order to meet various practical applications, metasurfaces with higher flexibility have become a research hotspot. Although both sensitive and insensitive metasurfaces have been proposed and extensively studied in recent years, metasurfaces with the ability to switch between sensitivity and insensitivity are rarely reported.
In this paper, a dual-channel mid-infrared toroidal metasurface that is composed of split equilateral triangular ring resonators is proposed. The reflectance spectrum of the metasurface is studied by both the finite-difference-time-domain (FDTD) method and temporal coupled-mode theory (TCMT). The results illustrate that the toroidal dipole metasurface has two channels. Channel A is insensitive to the polarization angle of the incident light and temperature, while channel B is sensitive to the polarization angle and temperature. Therefore, a two-bit programmable metasurface is achieved based on different sensitivities of the two channels to polarization angles. Furthermore, another metasurface is proposed for thermal imaging based on the different sensitivities of the two channels to temperature. The proposed metasurfaces offer alternative approaches to software-controlled digital metasurfaces and the photonic thermotropic.

Materials and Methods
The proposed metasurface is composed of periodical arraying unit cells and the unit cell is depicted in Figure 1a. A calcium fluoride (CaF 2 ) substrate is fabricated by chemical mechanical polishing. Then, a split equilateral triangular silicon (Si) ring is etched on the substrate [27] and the gap of the ring is filled with GST [28]. In addition, the ring is immersed in air. The thicknesses of the split ring and the gap are D = D 1 = 0.5 µm. There is an X-Y plane view of the unit cell in Figure 1b. The lattice constant of the unit cell is K = 3 µm. The outer side length of the triangular ring is L = 1.74 µm, whereas the inner side length is L 1 = 0.87 µm. The length and width of the gap are H = 0.7 µm and W = 0.27 µm, respectively. In Figure 1c, a linearly polarized incident plane wave illuminates vertically on the unit cell and a thermal field is applied on the unit cell. In other words, the unit cell is heated by an external heat source. As is shown in Figure 1d, the incidence angle θ of the incident plane wave is defined in the X-Z plane. The incidence angle θ varies between 0 • and 60 • and the initial incidence angle θ is set as 0 • . Optical responses of the metasurface are simulated by FDTD. Periodic boundary conditions are chosen in the Xand Ydirections and perfectly matched layers are set in the Zdirection, and this simulation model is semi-infinite in the Zdirection. The refractive indexes of Si and CaF 2 are obtained from Ref. [29].
metasurfaces have been proposed and extensively studied in recent years, metasurfaces with the ability to switch between sensitivity and insensitivity are rarely reported.
In this paper, a dual-channel mid-infrared toroidal metasurface that is composed of split equilateral triangular ring resonators is proposed. The reflectance spectrum of the metasurface is studied by both the finite-difference-time-domain (FDTD) method and temporal coupled-mode theory (TCMT). The results illustrate that the toroidal dipole metasurface has two channels. Channel A is insensitive to the polarization angle of the incident light and temperature, while channel B is sensitive to the polarization angle and temperature. Therefore, a two-bit programmable metasurface is achieved based on different sensitivities of the two channels to polarization angles. Furthermore, another metasurface is proposed for thermal imaging based on the different sensitivities of the two channels to temperature. The proposed metasurfaces offer alternative approaches to software-controlled digital metasurfaces and the photonic thermotropic.

Materials and Methods
The proposed metasurface is composed of periodical arraying unit cells and the unit cell is depicted in Figure 1a. A calcium fluoride (CaF2) substrate is fabricated by chemical mechanical polishing. Then, a split equilateral triangular silicon (Si) ring is etched on the substrate [27] and the gap of the ring is filled with GST [28]. In addition, the ring is immersed in air. The thicknesses of the split ring and the gap are D = D1 = 0.5 μm. There is an X-Y plane view of the unit cell in Figure 1b. The lattice constant of the unit cell is K = 3 μm. The outer side length of the triangular ring is L = 1.74 μm, whereas the inner side length is L1 = 0.87 μm. The length and width of the gap are H = 0.7 μm and W = 0.27 μm, respectively. In Figure 1c, a linearly polarized incident plane wave illuminates vertically on the unit cell and a thermal field is applied on the unit cell. In other words, the unit cell is heated by an external heat source. As is shown in Figure 1d, the incidence angle θ of the incident plane wave is defined in the X-Z plane. The incidence angle θ varies between 0° and 60° and the initial incidence angle θ is set as 0°. Optical responses of the metasurface are simulated by FDTD. Periodic boundary conditions are chosen in the X-and Ydirections and perfectly matched layers are set in the Z-direction, and this simulation model is semi-infinite in the Z-direction. The refractive indexes of Si and CaF2 are obtained from Ref. [29].  Affected by an external thermal field [30], GST is transformed from amorphous to crystalline. Permittivities of GST at amorphous (ε aGST ) and crystalline (ε cGST ) states are Nanomaterials 2022, 12, 3300 3 of 12 derived from [31]. The permittivity of GST (ε GST ) related with the crystalline fraction m (0 ≤ m ≤ 1) is described as [32]: In Figure 2, the temperature is set as 100 • and GST is in an amorphous state. The simulation results in Figure 2 are obtained when a Y-direction-linearly polarized incident plane wave illuminates vertically on the metasurface. The structural parameters in Figure 2 are the same as those in Figure 1. In Figure 2a, reflectance spectra of the metasurface are shown and the resonance wavelengths are at 3.51 µm and 3.82 µm, respectively. To further understand the resonances in the metasurface, TCMT is used to theoretically analyze the interactions between the metasurface and the incident light. In TCMT, normalized resonant amplitudes a 1 and a 2 as functions of the input amplitudes (S 1+ and S 2+ ) and the output amplitudes (S 1− and S 2− ) are expressed as [33,34]: , and B 2 = 0 1 1 0 . w 1 and w 2 represent the resonance frequencies; γ i1 , γ i2 , γ e1 , and γ e2 are the intrinsic decay rates and the external decay rates of different resonant modes, respectively. The coupling coefficient κ illustrates the coupling between two resonant modes. B 1 shows the coupling matrix between two resonant modes and the incident light. B T 1 is the transpose matrix of B 1 . B 2 is the direct scattering matrix. When the relationships between B 1 , B 2 , and γ e are satisfied as B 2 B * 1 = −B 1 and B * 1 B 1 = 2γ e , combined with Equations (2)-(4), the reflectance R can be concluded as: In Equation (4), because Si and CaF 2 have quite small optical loss [29], the absorption is small enough to ignore compared with reflectance. The calculation results are depicted with the black dotted line in Figure 2a, which is in good agreement with the simulation results.
In Figure 2b, Cartesian multipole moments are shown to measure the contributions to the optical resonances of the metasurface. The electric dipole P; magnetic dipole M; toroidal dipole T; electric quadrupole Q e ; magnetic quadrupole Q m ; scattering powers of each multipole moment I P , I M , I T , I Q e , I Q m ; reflectance R are obtained as follows [35][36][37][38][39]: where E(t) denotes the electric field of a metasurface unit cell extracted from FDTD simulation; ε d represents the relative permittivity of vacuum; c expresses the speed of light in vacuum; Cartesian coordinates α, β = x, y; µ 0 is the magnetic permeability of vacuum; k 0 is the wave vector; E 0 is the electric field of the incident light; subscript ( . . . ) // denotes the projection of the multipole moments into the X-Y plane; the direction of k 0 is the propagation of the incident light. Equation (11) illustrates the relationship between reflectance R, reflection coefficient r, and multipole moments. Here, we mainly study the contributions of P, M, T, Q e , and Q m , and higher-order multipoles are neglected. As shown in Figure 2b, two obvious toroidal dipole peaks are observed. The positions well-coincide with the reflection peaks in Figure 2a. As a result, the reflection peaks are mostly attributed to the toroidal dipoles. The far-field distributions at the resonance wavelengths are depicted in Figure 2c. The normalized electric field intensity distributions of the unit cell in the X-Y plane at the resonance wavelengths are shown in Figure 2d,e when Z = 0. The normalized electric field intensity distributions of the unit cell in the Y-Z plane at resonance wavelengths are shown in Figure 2f,g when X = 1.5 µm.

Wavefront Modulation
The impacts on reflectance spectra caused by the changes in geometrical parameters are shown in Figure 3. The split triangular ring can be seen as an L-C circuit of which the inductance and capacitance are inversely proportional to the resonance wavelengths of the metasurface [26]. Figure 3a shows W changing from 0.25 µm to 0.29 µm and the other structural parameters are the same as the parameters in. Because W represents the position occupied by the current flow, increasing W leads to a decrease in inductance of the gap of the ring and an obvious redshift in the right resonance wavelength. In Figure 3b, the two resonance wavelengths increase when H varies from 0.2 µm to 1.2 µm with the other parameters unchanged. An apparent redshift occurs in the right resonance wavelength with increasing reflectance. The gap can be considered as a capacitor with both sides accumulating negative and positive charges. As for a capacitor, an increment in Nanomaterials 2022, 12, 3300 5 of 12 the distance between the two electrodes leads to the growth of charging and discharging abilities. With the increase in H, the capacitance is reduced and the wavelength of the reflection peak is increased. Figure 3c illustrates that increasing L from 1.74 µm to 2.08 µm with the other parameters unchanged causes the reduction in the inductance of the ring and the redshifts in both resonance wavelengths with growing reflectance. As for Figure 3d, the growth of L 1 from 0.87 µm to 1.05 µm with the other parameters unchanged brings blueshifts in both resonance wavelengths without an apparent change in reflectance. This may be attributed to the growth of the induced current density of the ring.  Figure 4 shows a further study on the metasurface that has the same structural parameters as Figure 1. In Figure 4a, the reflectance spectra of the metasurface are evaluated at different polarization angles. The left peak is considered as channel A and the right peak is considered as channel B. When the polarization angle is 0°, it means that the electric field of a linearly polarized incident plane wave is parallel to the Y-direction, and when the polarization angle is 90°, the electric field is parallel to the X-direction. Increasing the polarization angle from 0° to 90°, from the Y-direction to the X-direction, the reflectance in channel A is invariant but the reflectance in channel B is reduced. Figure  4b shows the results calculated by TCMT, which has the same tendency with Figure 4a. Figure 4c shows normalized far-field scattering patterns of channel A with different polarization angles. Figure 4d shows the normalized far-field scattering patterns of channel B with different polarization angles. The reflected light is mainly limited in the ranges of 60° to 120° and 240° to 300° and the reflectance decreases with the change in polarization angle from 0° to 90°. With growth of the crystalline fraction of GST in the gap, the left resonance wavelength is fixed and the right has a redshift as shown in Figure 4e. Increasing the crystalline fraction causes growth of the effective permittivity of the gap. Growth of the effective permittivity can be reflected as the increase in W and H of the gap in Figure 3a,b. Therefore, changing the crystalline fraction has similar impacts on the reflectance spectrum of the metasurface as changing W and H. As is shown in Figure 4f, calculated by TCMT, the results have the same tendency as those in Figure 4e. In Figure  4g,h, the reflectance of the left and right resonances reduce when the crystalline fraction changes from 0 to 1. From Figure 4a to Figure 4h, it is concluded that channel A is insensitive to the change in polarization angle and crystalline fraction of GST. However, channel B is sensitive and the reflectance and resonance wavelength of channel B can be manipulated by the polarization angle and the crystalline fraction of GST independently. In Figure 4i, with the change in incidence angle θ from 0° to 60°, the two resonance wavelengths are fixed and the reflectance at the two resonance wavelengths is decreased.  Figure 4 shows a further study on the metasurface that has the same structural parameters as Figure 1. In Figure 4a, the reflectance spectra of the metasurface are evaluated at different polarization angles. The left peak is considered as channel A and the right peak is considered as channel B. When the polarization angle is 0 • , it means that the electric field of a linearly polarized incident plane wave is parallel to the Y-direction, and when the polarization angle is 90 • , the electric field is parallel to the X-direction. Increasing the polarization angle from 0 • to 90 • , from the Y-direction to the X-direction, the reflectance in channel A is invariant but the reflectance in channel B is reduced. Figure 4b shows the results calculated by TCMT, which has the same tendency with Figure 4a. Figure 4c shows normalized far-field scattering patterns of channel A with different polarization angles. Figure 4d shows the normalized far-field scattering patterns of channel B with different polarization angles. The reflected light is mainly limited in the ranges of 60 • to 120 • and 240 • to 300 • and the reflectance decreases with the change in polarization angle from 0 • to 90 • . With growth of the crystalline fraction of GST in the gap, the left resonance wavelength is fixed and the right has a redshift as shown in Figure 4e. Increasing the crystalline fraction causes growth of the effective permittivity of the gap. Growth of the effective permittivity can be reflected as the increase in W and H of the gap in Figure 3a,b. Therefore, changing the crystalline fraction has similar impacts on the reflectance spectrum of the metasurface as changing W and H. As is shown in Figure 4f, calculated by TCMT, the results have the same tendency as those in Figure 4e. In Figure 4g,h, the reflectance of the left and right resonances reduce when the crystalline fraction changes from 0 to 1. From Figure 4a to Figure 4h, it is concluded that channel A is insensitive to the change in polarization angle and crystalline fraction of GST. However, channel B is sensitive and the reflectance and Nanomaterials 2022, 12, 3300 6 of 12 resonance wavelength of channel B can be manipulated by the polarization angle and the crystalline fraction of GST independently. In Figure 4i, with the change in incidence angle θ from 0 • to 60 • , the two resonance wavelengths are fixed and the reflectance at the two resonance wavelengths is decreased. Figure 4j shows that the calculated results have the same tendency as those in Figure 4i. In Figure 4k,l, the reflectance of the left and right resonance wavelengths diminishes with the growth of incident angle θ from 0 • to 60 • . (i) reflectance spectra with different incidence angle θ; (j) reflectance spectra with different incidence angle θ calculated by TCMT; (k) far-field scattering patterns in X-Y plane of channel A with different incidence angle θ; (l) far-field scattering patterns in X-Y plane of channel B with different incidence angle θ.
Because the reflection peaks are mainly attributed to toroidal dipoles in Figure 2b and toroidal dipoles are calculated by current density in Equation (8), the surface current distributions are introduced to illustrate the physical mechanism of Figure 4a. The surface current distributions of the unit cell with different channels and polarizations are shown in Figure 5. The surface current distributions are evaluated at Z = 0. The color scale represents the normalized current density. White arrows in Figure 5 represent the directions of the electric field of the unit cell. According to Equation (5), the current has the same direction as the electric field. As a result, white arrows in the ring could represent directions of the currents. In channel A, when the unit cell is illuminated by X-directionlinearly polarized (XLP) light, surface currents are almost localized in two legs of the equilateral triangular ring and flow from positive (+) to negative (−). There are pathways 1 and 2 from "+" to "−". When the X-direction-linearly polarized light turns counterclockwise to Y-direction-linearly polarized (YLP) light by 90 degrees, currents exist in two legs and the base of the triangular ring. The positions of "+" and "−" turn counterclockwise by 90 degrees, and there are still two pathways 1 and 2 rotated by 90 degrees following the change in polarized light. In channel A, currents flow from "+" to "−" in two pathways, and the two pathways turn by 90 degrees with the change in the incident light from the X-direction to Y-direction. The two pathways are varied with the change in incident light, while the total amount of current inside the ring is almost (i) reflectance spectra with different incidence angle θ; (j) reflectance spectra with different incidence angle θ calculated by TCMT; (k) far-field scattering patterns in X-Y plane of channel A with different incidence angle θ; (l) far-field scattering patterns in X-Y plane of channel B with different incidence angle θ.
Because the reflection peaks are mainly attributed to toroidal dipoles in Figure 2b and toroidal dipoles are calculated by current density in Equation (8), the surface current distributions are introduced to illustrate the physical mechanism of Figure 4a. The surface current distributions of the unit cell with different channels and polarizations are shown in Figure 5. The surface current distributions are evaluated at Z = 0. The color scale represents the normalized current density. White arrows in Figure 5 represent the directions of the electric field of the unit cell. According to Equation (5), the current has the same direction as the electric field. As a result, white arrows in the ring could represent directions of the currents. In channel A, when the unit cell is illuminated by X-direction-linearly polarized (XLP) light, surface currents are almost localized in two legs of the equilateral triangular ring and flow from positive (+) to negative (−). There are pathways 1 and 2 from "+" to "−". When the X-direction-linearly polarized light turns counterclockwise to Y-directionlinearly polarized (YLP) light by 90 degrees, currents exist in two legs and the base of Nanomaterials 2022, 12, 3300 7 of 12 the triangular ring. The positions of "+" and "−" turn counterclockwise by 90 degrees, and there are still two pathways 1 and 2 rotated by 90 degrees following the change in polarized light. In channel A, currents flow from "+" to "−" in two pathways, and the two pathways turn by 90 degrees with the change in the incident light from the X-direction to Y-direction. The two pathways are varied with the change in incident light, while the total amount of current inside the ring is almost unchanged. This guarantees a relatively stable reflectance and insensitivity to polarization angles in channel A. As for channel B, with X-direction-linearly polarized light, few currents exist in the triangular ring. The direction of the electric field is almost unchanged inside the triangular ring. However, with Y-direction-linearly polarized light, currents emerge and are mainly sustained in the base of the ring. Currents flow counterclockwise in the ring. In channel B, with the incident light rotated from the X-direction to Y-direction, currents inside the ring emerge. This leads to an increase in the reflectance and sensitivity to polarization angles in channel B. unchanged. This guarantees a relatively stable reflectance and insensitivity to polarization angles in channel A. As for channel B, with X-direction-linearly polarized light, few currents exist in the triangular ring. The direction of the electric field is almost unchanged inside the triangular ring. However, with Y-direction-linearly polarized light, currents emerge and are mainly sustained in the base of the ring. Currents flow counterclockwise in the ring. In channel B, with the incident light rotated from the X-direction to Ydirection, currents inside the ring emerge. This leads to an increase in the reflectance and sensitivity to polarization angles in channel B.

Two-Bit Programmable Imaging
In Figure 4a, the reflectance spectra of the metasurface can be tuned with the change in the polarization angle. In Figure 6a, the metasurface is proposed for digital imaging according to Figure 4a. This metasurface is an 8 × 8 array of the unit cells with K = 2 μm. Quadrants 2 and 4 are formed by rotating quadrants 1 and 3 counterclockwise by 90 degrees. Figure 6b shows the magnetic field distributions when the metasurface is illuminated by X-and Y-polarized incident light. In channel A, it is shown that strong magnetic field distributions exist in all quadrants. As for channel B, strong magnetic distributions exist in quadrants 1 and 3 with X-direction-linearly polarized light, while the strong magnetic distributions exist in quadrants 2 and 4 with Y-direction-linearly polarized light. The quadrants with weak and strong magnetic field distributions can be coded as "0" and "1", respectively. In Figure 6b, channel A shows a stable pattern as "1111" and channel B expresses the digital images as "1010" and "0101".

Two-Bit Programmable Imaging
In Figure 4a, the reflectance spectra of the metasurface can be tuned with the change in the polarization angle. In Figure 6a, the metasurface is proposed for digital imaging according to Figure 4a. This metasurface is an 8 × 8 array of the unit cells with K = 2 µm. Quadrants 2 and 4 are formed by rotating quadrants 1 and 3 counterclockwise by 90 degrees. Figure 6b shows the magnetic field distributions when the metasurface is illuminated by Xand Y-polarized incident light. In channel A, it is shown that strong magnetic field distributions exist in all quadrants. As for channel B, strong magnetic distributions exist in quadrants 1 and 3 with X-direction-linearly polarized light, while the strong magnetic distributions exist in quadrants 2 and 4 with Y-direction-linearly polarized light. The quadrants with weak and strong magnetic field distributions can be coded as "0" and "1", respectively. In Figure 6b, channel A shows a stable pattern as "1111" and channel B expresses the digital images as "1010" and "0101". To quantitatively assess the quality of the generated images, imaging efficiency  is introduced. i P represents the reflected power of each quadrant and r P denotes the reflected power of the metasurface.  is expressed as: The imaging efficiency in Figure 6b can be calculated. In channel A, four quadrants all have  = 25% with both X-and Y-polarization, while in channel B,  of zone "0" and "1" are 14% and 36%, respectively.

Thermal Imaging
According to Figure 4c, the reflectance spectra of the metasurface are varied with change in the crystalline fraction of GST in the gap. As the crystalline fraction of GST is in direct proportion to the external temperature, we design a metasurface that can be used for thermal imaging. The structure has an array of 10 × 10 unit cells with K = 3 μm, as shown in Figure 7a. In Figure 7b, the crystalline fraction in each unit cell of the metasurface is shown. The crystalline fraction changes from 0 to 1 in a step of 0.25. The crystalline fraction distribution in Figure 7b is used to simulate thermal distribution on the metasurface. To quantitatively assess the quality of the generated images, imaging efficiency η is introduced. P i represents the reflected power of each quadrant and P r denotes the reflected power of the metasurface. η is expressed as: The imaging efficiency in Figure 6b can be calculated. In channel A, four quadrants all have η = 25% with both Xand Y-polarization, while in channel B, η of zone "0" and "1" are 14% and 36%, respectively.

Thermal Imaging
According to Figure 4c, the reflectance spectra of the metasurface are varied with change in the crystalline fraction of GST in the gap. As the crystalline fraction of GST is in direct proportion to the external temperature, we design a metasurface that can be used for thermal imaging. The structure has an array of 10 × 10 unit cells with K = 3 µm, as shown in Figure 7a. In Figure 7b, the crystalline fraction in each unit cell of the metasurface is shown. The crystalline fraction changes from 0 to 1 in a step of 0.25. The crystalline fraction distribution in Figure 7b is used to simulate thermal distribution on the metasurface.
In Figure 7b, the thermal imaging metasurface is divided into five parts according to different crystalline fractions. Reflectance spectra of five parts are shown in Figure 8a. For example, the black dotted line in Figure 8a denotes the reflectance spectrum of the part in which the crystalline fraction is 0.0 in Figure 7b. In channel A, the reflection peaks are almost localized at the wavelength λ 0 . However, in channel B, the reflection peaks increase from λ 1 to λ 5 with the change in crystalline fraction. In Figure 8b, for channel A, the magnetic field distribution at reflection peak λ 0 is shown and the imaging efficiencies of crystalline fractions 0.00, 0.25, 0.50, 0.75, and 1.00 are 23%, 23%, 29%, 10%, and 15%, respectively. In Figure 8c-g, for channel B, strong magnetic distributions at reflection peaks λ 1 , λ 2 , λ 3 , λ 4 , and λ 5 are mainly localized in the field of crystalline fractions of 0.00, 0.25, 0.50, 0.75, and 1.00, respectively. The maximal imaging efficiencies at reflection peaks λ 1 , λ 2 , λ 3 , λ 4 , and λ 5 are 40%, 37%, 40%, 37%, and 46%, respectively. Each reflection peak could represent a crystalline fraction and a temperature. In this way, thermal imaging is realized because the thermal distribution can be reflected with magnetic field distributions at different reflection peaks in channel B. In other words, the thermal distribution on the metasurface can be illustrated by the reflected light distributions. In Figure 7b, the thermal imaging metasurface is divided into five parts according to different crystalline fractions. Reflectance spectra of five parts are shown in Figure 8a. For example, the black dotted line in Figure 8a denotes the reflectance spectrum of the part in which the crystalline fraction is 0.0 in Figure 7b. In channel A, the reflection peaks are almost localized at the wavelength λ0. However, in channel B, the reflection peaks increase from λ1 to λ5 with the change in crystalline fraction. In Figure 8b, for channel A, the magnetic field distribution at reflection peak λ0 is shown and the imaging efficiencies of crystalline fractions 0.00, 0.25, 0.50, 0.75, and 1.00 are 23%, 23%, 29%, 10%, and 15%, respectively. In Figure 8c-g, for channel B, strong magnetic distributions at reflection peaks λ1, λ2, λ3, λ4, and λ5 are mainly localized in the field of crystalline fractions of 0.00, 0.25, 0.50, 0.75, and 1.00, respectively. The maximal imaging efficiencies at reflection peaks λ1, λ2, λ3, λ4, and λ5 are 40%, 37%, 40%, 37%, and 46%, respectively. Each reflection peak could represent a crystalline fraction and a temperature. In this way, thermal imaging is realized because the thermal distribution can be reflected with magnetic field distributions at different reflection peaks in channel B. In other words, the thermal distribution on the metasurface can be illustrated by the reflected light distributions. In Table 1, we collect the previous works and compare the results with our work to show the advantages of the proposed structure in terms of wavelength range, functions, and imaging applications. Our work can be manipulated with dual channels, polarization angle, and temperature. In addition, two-bit coding and thermal imaging metasurfaces are achieved based on multiple manipulations. In other words, our designed metasurface In Table 1, we collect the previous works and compare the results with our work to show the advantages of the proposed structure in terms of wavelength range, functions, and imaging applications. Our work can be manipulated with dual channels, polarization angle, and temperature. In addition, two-bit coding and thermal imaging metasurfaces are achieved based on multiple manipulations. In other words, our designed metasurface is more tunable and can be designed for thermal imaging compared with others. Table 1. Comparison of our work with analogous structures.

Conclusions
In conclusion, a dual-channel mid-infrared toroidal metasurface is designed with split equilateral triangular ring resonators. Calculated by TCMT, the toroidal dipole metasurface has two channels of which one is sensitive and another is insensitive to the polarization angle of incident light and temperature. A coding metasurface is realized with Xor Ypolarization expressed as "0" or "1" in digital mode. Another metasurface is proposed to achieve thermal imaging that can represent the external thermal distribution with a reflected wavelength distribution. The proposed metasurfaces offer alternative approaches to many information processing and energy-harvesting applications.