Theoretical Study on All-Dielectric Elliptic Cross Metasurface Sensor Governed by Bound States in the Continuum

Abstract

The appearance of all-dielectric micro–nano photonic devices constructed from high refractive index dielectric materials offers a low-loss platform for the manipulation of electromagnetic waves. The manipulation of electromagnetic waves by all-dielectric metasurfaces reveals unprecedented potential, such as focusing electromagnetic waves and generating structured light. Recent advances in dielectric metasurfaces are associated with bound states in the continuum, which can be described as non-radiative eigen modes above the light cone supported by metasurfaces. Here, we propose an all-dielectric metasurface composed of elliptic cross pillars arranged periodically and verify that the displacement distance of a single elliptic pillar can control the strength of the light–matter interaction. Specifically, when the elliptic cross pillar is C₄ symmetric, the quality factor of the metasurface at the Γ point is infinite, also called the bound states in the continuum. Once the C₄ symmetry is broken by moving a single elliptic pillar, the corresponding metasurface engenders mode leakage; however, the large quality factor still exists, which is called the quasi-bound states in the continuum. Then, it is verified by simulation that the designed metasurface is sensitive to the refractive index change of the surrounding medium, indicating that it can be applied for refractive index sensing. Moreover, combined with the specific frequency and the refractive index variation of the medium around the metasurface, the information encryption transmission can be realized effectively. Therefore, we envisage that the designed all-dielectric elliptic cross metasurface can promote the development of miniaturized photon sensors and information encoders due to its sensitivity.

1. Introduction

The development of metamaterial photonics has been greatly promoted since Pendry et al. proposed a composite structure composed of high refractive index dielectric cylinders and verified that the structure has effective negative permeability [1]. The typical metamaterial is the metal split-ring resonator operating in the microwave band. However, due to the inherent ohmic loss of metals, potential applications of metal metamaterials are hindered [2,3]. Excitingly, the proposal of all-dielectric metasurfaces based on high refractive index dielectric materials provides an excellent platform for the exploration of lowloss subwavelength photonic devices [4]. In recent years, all-dielectric metasurfaces have been widely harnessed in many fields, such as beam focusing [5,6,7,8,9], holography [10,11,12], imaging [13,14,15], structured light generators [16,17,18,19], etc.

The function of restricting electromagnetic waves possessed by numerous metasurfaces can be attributed to strong resonance, whose corresponding quality factor (Q factor) is an important indicator to measure the strength of light–matter interaction. Generally, tailoring the shape, size and periodicity of all dielectric meta-atoms can increase the Q factor and enhance the local electromagnetic field, which can be employed in a wide range of realms, such as sensing [20,21,22], optical switch [23], and high-order harmonic generation [24]. Recently, the proposal of metasurfaces supporting bound states in the continuum (BICs) provides a new scheme for further enhancing the light–matter interaction. Specifically, under the incidence of electromagnetic waves, some specific modes generated by metasurfaces above the light cone carry infinite Q factors, which cannot be excited by radiation from the far field [25,26]. However, for practical optical devices, due to the influence of material absorption and fabricating technology, they have finite Q factors and narrow resonance linewidths, called Quasi BICs. The general strategy to achieve Quasi BICs is to break the structural symmetry of the resonant unit [27,28,29,30]. In other words, rotating the resonant unit in a specific modality or adding (reducing) a part of the initial resonant unit can generate Quasi BICs with sharp Fano lines under the incidence of electromagnetic waves [31,32,33,34]. For example, Sinclair et al. proposed that the removal of highly symmetrical resonant units can induce coupling between orthogonal modes, thereby generating Fano resonance with a high-quality Q factor. Furthermore, they verified the feasibility of the scheme through experiments [31]. Tittl et al. broke the C₂ symmetry of the resonant unit by changing the tilt angle of a pair of elliptic pillars, and produced Quasi BICs which can be exploited for highly sensitive molecular sensors [33]. In addition, by using a metasurface that supports a symmetry-protected Quasi BIC, the near-field imaging can be realized [34].

Here, we propose an all-dielectric metasurface consisting of elliptic cross silicon pillars with C₄ symmetry arranged periodically. By calculating the eigenmodes and transmission spectra excited by different modes of electromagnetic waves in C₄ symmetry, it is verified that the metasurface has an infinite Q factor at the Γ point, that is, the metasurface supports BIC. After introducing an asymmetric parameter to break C₄ symmetry, the metasurface can generate sharp Fano resonance under the incidence of electromagnetic waves. The relationship between the radiative Q factor and the asymmetric parameters manifests that the BICs supported by the metasurface are symmetry protected BICs. Once different refractive index mediums are introduced into the gap of the elliptic cross silicon pillars, the metasurface exhibits a sensitivity to the change of the refractive index. Moreover, based on the high Q properties of Quasi BICs, information encoding can be realized effectively by selecting specific meta-atoms to construct a metasurface. By changing the refractive index of the medium around the metasurface, information encryption of near-field imaging is realized, which effectively verifies the feasibility of information encoding. These above schemes provide a new working platform for real-time liquid refractive index detection and information encoding.

2. Results and Discussion

Here, in order to show the excitation effect of symmetry on BICs supported by the metasurface, we propose an all-silicon elliptic cross metasurface (refractive index n₂ = 3.45) placed on a substrate composed of silica (refractive index n₁ = 1.45), whose corresponding schematic diagram is illustrated in Figure 1. It can be seen from Figure 1 that the thickness of the metasurface is T = 150 μm, and the thickness of the corresponding substrate is h = 50 μm. The period of the meta-atom is p = 300 μm, and the corresponding major and minor semi-axes of a single elliptic pillar are L = 120 μm and W = 30 μm, respectively. When the values of a₁ and a₂ are equal, the meta-atom is C₄ symmetric. Moving a single elliptic pillar makes a₁ ≠ a₂, which can break the C₄ symmetry of the meta-atom. Next, we will analyze the BICs supported by the designed metasurface under the incidence of TM and TE modes, respectively.

It is worth emphasizing that the BICs supported by a metasurface have a strong interaction with infinite Q factors and zero leakage. Electromagnetic waves of a particular frequency coexist with extended states in the radiative continuum, while remaining perfectly localized in the absence of radiation. Therefore, BICs with radiation leakage can be obtained by adjusting some geometric or structural parameters of the metasurface. One of the more common ones is to break the symmetry of the structure and obtain the symmetry-protected quasi-BICs. Thus, in order to verify that the designed metasurface can support symmetry-protected BICs, COMSOL Multiphysics software is harnessed to build a single meta-atom model and calculate the eigenmodes (see Methods). Figure 2a shows two band diagrams of the TM-like modes in the case of C₄ symmetry. Among them, the radiation Q factors of TM-like 1 and TM-like 2 modes are illustrated in Figure 2b,c, respectively. The frequency eigenvalue of quasi-BIC solved in COMSOL software (COMSOL Multiphysics 5.6, Burlington, MA, USA) is a complex number, expressed as ω-iγ, with the real part representing the frequency center and the imaginary part representing the radiation rate. Therefore, the Q value can be obtained based on the above results. Here, the equation for calculating Q factor is: Q = ω/2γ, where ω represents the resonant frequency and γ represents the corresponding radiation rates. Apparently, these two modes have an infinite Q factor at Γ point. The insets show the electric field distribution of the corresponding |E_z|: when the frequencies are 0.611 and 0.641 THz, the electric field at the Γ point is confined to the elliptic cross pillar. The electric field of TM-like 1 mode is concentrated at the center of the elliptic cross pillar, while the electric field of TM-like 2 mode is mainly distributed at the boundary. Then, the C₄ symmetry is broken by moving a single elliptic pillar of the elliptic cross pillar along the x-direction. Here, the asymmetric parameter is defined as: α = |(a₁ − a₂)/(a₁ + a₂)|. Furthermore, the movement distance of a single elliptic pillar along the x-direction is: d ≈ |a₁ − a₂|/2. As shown in Figure 2d, there is an inverse relationship between the radiation Q factor of TM-like 1 mode and the square of the asymmetric parameter α, which is consistent with reference [32]. Figure 2e shows the radiation Q factor of TM-like 2 mode corresponding to different asymmetric parameters. It can be concluded that the radiation Q factor of such mode is basically inversely proportional to the square of the asymmetric parameter α. Figure 2f,h shows the electric field distribution of the TM-like 1 mode in the xoz plane with the asymmetric parameters α = 0 and α = 0.8. Of particular emphasis, the absolute value of electric field is considered here. It can be seen that in the case of C₄ symmetry, TM-like 1 mode is bound inside the entire elliptic cross pillar, and after breaking the C₄ symmetry, TM-like 1 mode will leak into the far field. However, part of the electric field is still distributed inside the elliptic cross pillar, and such a mode is called a Quasi BIC 1. Similarly, we obtained the electric field distribution of TM-like 2 mode in xoz plane with asymmetric parameters α = 0 and α = 0.8 through simulation, as shown in Figure 2g,i. The results exhibit that the TM-like 2 mode is confined inside the dielectric pillar with the asymmetric parameter α = 0, while the strong leakage occurs with the asymmetric parameter α = 0.8. The above results demonstrate that the TM-like 1 and TM-like 2 modes supported by the designed metasurface in the case of C₄ symmetry are BIC 1 and BIC 2 at Γ points, respectively, and these two BICs belong to the symmetry-protected BICs.

Next, CST Studio Suite software (CST STUDIO SUITE 2019, Computer Simulation Technology AG, Darmstadt, Germany) is used to build a single elliptic cross pillar (see Methods). We set TM mode input at the incident port and TM mode output at the outgoing port. The transmitted amplitude of TM mode at different frequencies in Figure 3a is obtained by employing a frequency domain solver. Here, the thickness of the substrate is maintained at 50 μm. It is worth noting that when the movement distance d of a single elliptic pillar is not zero, Fano resonances will occur around 0.61 and 0.64 THz, respectively. When d increases from 10 μm to 60 μm, the lower resonant frequency decreases by 4.9 GHz, which means that the increase in d will lead to a blue shift in the resonant frequency. The linewidth of the Fano resonance decreases with the decrease in d, and disappears in the case of d = 0 μm, again confirming the presence of BIC 1 and BIC 2 in Figure 2b,c. In addition, when the movement distance of a single elliptical pillar is d = 20 μm, the amplitude changes of Fano resonances generated near 0.61 and 0.64 THz is more intense, which may be applied to narrowband bandpass filtering. With this, the transmitted amplitude of the TM mode with different thicknesses of substrate h is further studied. Here, the control variate methods are harnessed to maintain the movement distance of the elliptic pillar as d = 20 μm. The transmitted amplitude of TM mode is calculated when the thicknesses of substrate are h = 50 μm, 70 μm, 90 μm, 110 μm, 130 μm and 150 μm, as shown in Figure 3b. It can be seen from Figure 3b that with the increase in the thickness of the substrate h, the Fano resonance has a redshift. Meanwhile, with the change of thickness of the substrate h, the maximum amplitude variation in Fano resonance near 0.61 THz is inapparent, while the maximum amplitude variation in Fano resonance near 0.64 THz is more obvious. Therefore, the change of thickness of the substrate h will cause the frequency shift of the transmitted spectrum of TM mode.

Two cases of TM-like mode in elliptic cross pillar are discussed and analyzed above. Next, we will further analyze the TE-like mode. Similar to the method adopted in Figure 2a, the band diagram of TE-like mode under C₄ symmetry is obtained through COMSOL Multiphysics software (see Methods, COMSOL Multiphysics 5.6, Burlington, MA, USA), as shown in Figure 4a. Figure 4b shows the corresponding radiation Q factor, while the illustration is the magnetic field distribution (|H_z|) at Γ point (BIC 3). It can be seen from the illustration that when the frequency is 0.646 THz, the magnetic field |H_z| is mainly concentrated in the center of the elliptic cross pillar. Figure 4c shows the relationship curve between the radiation Q factor and the asymmetric parameter α after moving a single elliptic pillar: similar to the results in Figure 2c, the radiation Q factor is basically inversely proportional to the square of the asymmetric parameter α. In addition, Figure 4d,e show the |H| distributions of the TE-like mode on the xoz plane with the asymmetry parameter α = 0 and α = 0.8, respectively. Clearly, the magnetic field is bound in the dielectric pillar with the asymmetry parameter α = 0, while strong magnetic field leakage occurred with the asymmetry parameter α = 0.8. In short, BIC 3 at Γ point (0.646 THz) in Figure 4b is a symmetry-protected BIC.

Further, a single elliptic cross pillar is constructed using CST Studio Suite software with the frequency range of 0.6–0.7 THz. We set TE mode input at the incident port and TE mode output at the outgoing port. The transmitted amplitude curve of TE mode in Figure 5a is obtained by employing the frequency domain solver. Here, the thickness of the base is maintained as h = 50 μm. Once the movement distance d of the elliptic pillar is increased, the linewidth of Fano resonance gradually increases, and the transmitted spectrum has a slight blue shift. Apparently, when d = 0 μm, the linewidth of Fano resonance in the transmitted spectrum disappears, which verifies that the designed elliptic cross pillar can support BIC 3 under C₄ symmetry. In addition, in the case of d ≠ 0 μm, the Fano resonance in the transmitted spectrum is called Quasi BIC 3. The movement distance d of the elliptic pillar is set as 20 μm, and the transmitted amplitude of TE mode with different thickness h has been calculated (see Figure 5b). Figure 5b manifests that with the increase in h, Fano resonance in transmitted amplitude curve of TE mode has a blue shift, while the corresponding linewidth has hardly changed. Therefore, it can be determined that the change of h mainly affects the frequency of Fano resonance, but hardly change the linewidth of Fano resonance.

To evaluate the sensing performance of the designed all-dielectric elliptic cross metasurface, the refractive index sensing capability of the metasurface is quantified. Concretely, a medium with different refractive index n is filled in the gap of the metasurface to simulate the situation of immersing the metasurface into a solution mixture with different refractive index. As shown in Figure 6a, the thickness of the substrate is set as h = 50 μm, and the movement distance of a single elliptic pillar is d = 20 μm, then the transmitted spectra at n = 1.1, 1.15, 1.2, 1.25, 1.3 and 1.35 are obtained through simulation. When the refractive index is n = 1.1, the transmission curve of the designed metasurface under the incidence of TM mode engenders Fano resonance at 0.593 THz (Quasi BIC 1) and 0.631 THz (Quasi BIC 2). In addition, with the increase in the refractive index n, the two Fano resonances in the transmitted curve have red shifted, and the linewidths of Quasi BIC 1 and Quasi BIC 2 exist slightly changed. To further measure the performance of the sensor, the sensitivity (S = |Δf/Δn|) is calculated according to the frequency offset of Quasi BIC 1 and Quasi BIC 2 in the transmitted curve with different refractive index n in Figure 6a, as plotted in Figure 6b. Among them, Δn is the difference between adjacent refractive indices (in Figure 6a), and Δf is the frequency shift of the resonant frequency in response to the refractive index change Δn. The blue curve and the red curve represent the sensitivity of Quasi BIC 1 and Quasi BIC 2, respectively. It is not difficult to find that the sensitivity of Quasi BIC 1 is always higher than 85 GHz/RIU, and the sensitivity of Quasi BIC 2 is in the range of 65–77 GHz/RIU. Therefore, under the incidence TM mode incidence, the sensing performance of Quasi BIC 1 generated by the metasurface is better than that of Quasi BIC 2.

Furthermore, the simulated results of the metasurface filled with different media under the illumination of TE mode are obtained through CST Studio Suite software (see Figure 6c); with the increase in the refractive index n of the filled media, the resonant frequency at Quasi BIC 3 has a significant red shift, and the resonant linewidth increases accordingly. As the refractive index n increases from 1.1 to 1.35, the corresponding frequency of Quasi BIC 3 moves from 0.635 to 0.599 THz. Figure 6d shows the sensitivity after the change of refractive index n. According to Figure 6d, it can be determined that under the incidence of TE mode, when the difference between adjacent refractive indices remains the same, increasing the refractive index of the filled medium can yield greater sensitivity. In addition, the sensitivity of excited Quasi BIC 3 after the incidence of TE mode is always greater than 120 GHz/RIU, signifying that Quasi BIC 3 has better sensing performance than excited Quasi BIC 1 and 2 after the incidence of TM mode. Wang et al. ‘s work mentions that sensitivity on the order of 10² GHz/RIU can be used for sensing [35], so it can be concluded that when the thickness of substrate is h = 50 μm and the movement distance of a single elliptic pillar is d = 20 μm, the constructed metasurface is sensitive to the refractive index of the surrounding medium, which may be suitable for sensing. It is worth emphasizing that the sensor with three peaks can improve the sensing accuracy compared with the sensor with only one peak under two incident modes. Apart from this, to facilitate the comparison of the sensing performance of sensors applied to different wavebands, Figure of Merit (FOM) is proposed, which is defined as [36]: FOM = S/FWHM, where FWHM is full width at half maxima of resonant peak. It is worth emphasizing that the greater the value of FOM, indicating that the sensor has better sensing performance. Thus, according to the above formula, the FOM of the three Quasi BICs in Figure 6b,d are calculated (see Figure 6e,f). It can be found that the FOM of Quasi BIC 1 is higher than that of the other two Quasi BICs, and the corresponding FOM reached a maximum of 93 RIU⁻¹. Interestingly, with the increase in the refractive index n, the sensitivity and FOM of the three Quasi BICs all illustrated an increasing trend, which verified that the interaction between the medium around the metasurface and the electric field is enhanced with the increase in the refractive index n. Meanwhile, performance comparison between such sensor and some other recently developed sensors [37,38,39,40,41,42] is listed in Table 1. Comparatively, the sensing performance of the designed metasurface is not optimal, which may be attributed to the characteristics of the material itself [43].

The above results have confirmed that the resonance in the transmission spectrum of the designed element at TM (TE) incidence shifts with the change of refractive index of the medium around the element atoms. Here, we take advantage of this property to encode the information based on near-field imaging [34]. As shown in Figure 7a, three representative meta-atoms are carefully selected, denoted as A, B and C respectively. Thereinto, A is the unit with C₄ symmetry, B is the unit formed after moving a single elliptic pillar by 20 μm, and C is the unit formed after rotating B by 90° counterclockwise. For efficient numerical coding, the transmission amplitudes below 0.45 are defined as 0, and the transmission amplitudes above 0.45 are defined as 1. Figure 7b plots the numerical codes of these three meta-atoms at frequencies of 0.607, 0.636 and 0.645 THz under TM incidence when the refractive index of the surrounding medium is n = 1. Meanwhile, Figure 7c,d display the numerical results when the refractive index of surrounding medium is 1.1 and 1.3 respectively. It can be seen that as n changes, the numerical information of these meta-atoms at the same frequency changes. Furthermore, Figure 7e–g shows the results of the selected atoms under TE incidence when the refractive index of the surrounding medium is 1, 1.1 and 1.3. With these, we construct a metasurface to show the function of information coding (see Figure 7h), which is composed of 18 × 18 units. It can be seen that the area in the black box is composed of meta-atom B, the area in the yellow box is composed of meta-atom C, and the other areas are composed of meta-atom A.

Then, the near-field intensity distribution of the metasurface designed in the air under the normal incidence of TM wave is obtained by simulation, as shown in Figure 8a (where TM incidence is considered as an example). It is obvious that the designed metasurface has the same intensity distribution when the frequencies are 0.607 and 0.636 THz: some regions with stronger energy form a pattern of ‘-’, while others have weaker energy. In addition, when the frequency is 0.645 THz, some regions with stronger energy form a pattern of ‘|’. Figure 8b depicts the intensity distribution at three frequencies under TM incidence when n = 1.1: when the frequency is 0.636 THz, the corresponding intensity distribution carries obvious patterns, while no obvious patterns are found at the frequencies of 0.607 and 0.645 THz. Similarly, when the refractive index of the surrounding medium is adjusted to 1.3, the designed metasurface has different intensity distributions at the above three frequencies. It can be found that when the frequency is 0.607 THz, the intensity distribution generated has a ‘|’ pattern, while the intensity distribution of the other two frequencies has no obvious pattern. In conclusion, the metasurface constructed based on three meta-atoms can generate specific near-field images under TM and TE incidences. Meanwhile, by changing the refractive index of the medium around the metasurface, different intensity distributions can be generated at the same frequency, which may be used for information encryption.

It is common knowledge that ordinary resonance generally has a wide working frequency band, which will cause that even if the central frequency is unknown during decryption, the result close to the real information may be obtained in a larger bandwidth near the central frequency. The advantages of a Quasi-BIC mode are high Q factor and narrow working frequency band, so it is necessary to know the monitoring frequency accurately in advance for decryption, which enhances the difficulty of decryption. With this, we proposed two simple application conceptions for the mentioned information encryption. First, encoding information encryption: by defining ‘-’ as 2, ‘|’ as 1, and none as 0, a ternary cipher compilation system can be obtained. This system can encrypt information through frequency, polarization direction and refractive index. For example, defines a code 3611 TM, it corresponds to 0.636 TH_Z, 1.1 refractive index and TM polarization (such information can be designed to be more secretive). Then the 3611 TM encrypted information can be converted to the near-field field strength of ‘|’, that is, 1. Suppose that a letter A needs to be transmitted, and the binary code of A is ‘01000001′, which is converted to hexadecimal ‘2102′. For a unique 16 × 16 lattice, it can be compiled through 3610 TM, 3611 TM, 4513 TM, 0710 TM (where TM incidence is considered as an example). The advantage of this design is that the same information can be obtained in many ways to enhance deciphering difficulty. Second is the graphics encryption: by placing independent regions similar to the B and C regions at different angles, a basic graph can be constructed with ‘-’, and the near-field graph can only be decrypted with specific polarization, frequency and refractive index values, while the wrong polarization, frequency and refractive index information cannot obtain the correct graph. Hence, such above coding encryption scheme can be employed as a supplement to hardware encryption equipment and software encryption, improving the confidentiality of data.

3. Conclusions

In conclusion, we propose and verify that the all-dielectric metasurface constructed by elliptic cross pillars with C₄ symmetry can support three different symmetry-protected BICs. Under the excitation of the TM mode, two Fano resonances can appear by adjusting the movement distance of a single elliptic pillar, and the linewidths of these two resonances gradually decrease and disappear with the decrease in movement distance, accompanied by a slight frequency shift. The change in thickness of the substrate mainly affects the shift of resonance frequency. Similarly, under the excitation of the TE mode, the change of movement distance mainly changes the resonant linewidth of Quasi BIC 3, while the thickness of the substrate mainly affects the corresponding frequency of Quasi BIC 3. After introducing different refractive index media between adjacent elliptic cross pillars, the transmitted spectrum generated by the metasurface exists as a significant shift. Sensitivity is introduced to quantify the sensitive ability of the metasurface to the change of the refractive index. The simulated results show that the sensitivity of the designed metasurface under the incidence of the TE mode is greater than 120 GHz/RIU, indicating that this design scheme may promote the development of real-time liquid refractive index detection and biosensors. In addition, information can be encoded effectively with high encryption by combining the specific frequency with the refractive index of the medium around the metasurface.

Methods: Eigenmode solution: After building the model, periodic boundary conditions are applied along the x- and y-directions. Such a model is symmetrical up and down, so we can have two perfect sets of energy bands: TE-like and TM-like. TE-like and TM-like are derived from TE and TM modes. Because the designed meta-atom does not have continuous translational symmetry in the z-direction, there is no real TE or TM mode, but it can be based on the symmetry plane of the meta-atom in the z-direction Γ. The direction of the point electric field and magnetic field defines the entire energy band as TE-like or TM-like, but the above premise is that the entire atomic element is mirror symmetric in the z-direction. If the meta-atom is located in the z-direction without mirror symmetry, such as air on one side and substrate on the other, the energy band obtained is neither perfectly TE-like or TM-like. Here, the model we choose is with the atoms as all air, so it has perfect two sets of energy bands. This requires us to add a perfect magnetic conductor (PMC) or perfect electric conductor (PEC) on the high symmetry plane, and at the same time, we can simplify the model—only need to calculate half of the region. The boundary condition set on its lower surface is PMC (TE-like) or PEC (TM-like). It should be noted that the thickness of the designed atomic element should be reduced by half. Then, different eigen modes and corresponding Q factors can be calculated by the eigen mode solver.

Spectrum solution: In CST Studio Suite software, the x-direction and y-direction are set as unit cell, the z-direction is set as open boundary, and the frequency range is set as 0.56–0.67 THz.

Author Contributions

Conceptualization, methodology, software, writing—original draft: H.C.; data curation, validation: X.Y.; writing—review & editing, supervision: L.M. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Table

Figure 1. Schematic diagram of all-dielectric metasurface supporting BICs. The metasurface is composed of an elliptic cross silicon pillar, with a thickness of T = 150 μm. The substrate is made of silica, and the corresponding thickness is h = 50 μm. The illustration is a top view of a single atomic element. The period is p = 300 μm, and the corresponding major and minor semi-axes of a single elliptic pillar are L = 120 μm and W = 30 μm, respectively.

Figure 2. (a) Band diagram of TM-like mode in an elliptic cross pillar with C4 symmetry. The illustration shows a k-space diagram, and Γ point shows a highly symmetric point. (b,c) The radiation Q factor of TM-like 1 and TM-like 2 modes in an elliptic cross pillar with C4 symmetry, and the illustration shows the electric field distribution at Γ point. (d,e) The correlation curve of radiation Q factor of TM-like 1 (TM-like 2) mode and asymmetric parameter α. (f,g) When the asymmetric parameters meet the condition: α = 0, the electric field distributions of TM-like 1 and TM-like 2 modes at xoz plane. (h,i) When the asymmetric parameters meet the condition: α = 0.8, the electric field distributions of TM-like 1 and TM-like 2 modes at xoz plane.

Figure 3. (a) When the movement distances of a single elliptic pillar are d = 0 μm, 10 μm, 20 μm, 30 μm, 40 μm, 50 μm and 60 μm, the transmission amplitude under the incidence of TM mode. The thickness of the substrate is h = 50 μm. (b) When the thicknesses of the substrate are h = 50 μm, 70 μm, 90 μm, 110 μm, 130 μm and 150 μm, the transmission amplitude under the incidence of TM mode. The movement distance of the single elliptic pillar is d = 20 μm.

Figure 4. (a) Band diagram of transverse electric-like (TE-like) mode in an elliptic cross pillar with C4 symmetry. (b) The radiation Q factor of TE-like mode in an elliptic cross pillar with C4 symmetry, and the illustration shows the electric field distribution at Γ point. (c) The correlation curve of radiation Q factor of TE-like mode and asymmetric parameter α. (d) When the asymmetric parameters meet the condition: α = 0, the magnetic field distribution of TE-like mode. (e) When the asymmetric parameters meet the condition: α = 0.8, the magnetic field distribution of TE-like mode.

Figure 5. (a) When the movement distances of a single elliptic pillar are d = 0 μm, 10 μm, 20 μm, 30 μm, 40 μm, 50 μm and 60 μm, the transmission amplitude under the incidence of TE mode. The thickness of the substrate is h = 50 μm. (b) When the thicknesses of the substrate are h = 50 μm, 70 μm, 90 μm, 110 μm, 130 μm and 150 μm, the transmission amplitude under the incidence of TE mode. The movement distance of the single elliptic pillar is d = 20 μm.

Figure 6. (a) The relation between the transmitted amplitude and frequency generated by the meta-atom under the incidence of TM mode. Here, n represents the refractive index of the material in the gap between adjacent elliptic cross pillars. (b) The sensitivities of Quasi BIC 1 and 2 in TM-like mode is calculated by taking the refractive index difference between two adjacent refractive indexes as a unit. (c) The relation curve between transmitted amplitude and frequency generated by meta-atom under the illumination of TE mode. (d) Sensitivity of Quasi BIC 3 engendered by meta-atom under the incidence of TE mode. (e,f) Figure of merit (FOM) of the three Quasi BICs under the incidence of TM and TE modes.

Figure 7. (a) The three representative meta-atoms harnessed to encode information are denoted as A, B and C. (b–d) When the refractive index n = 1, 1.1 and 1.3, the number of codes of the meta-atom at different frequencies under TM incidence. (e–g) When the refractive index n = 1, 1.1 and 1.3, the number of codes of the meta-atom at different frequencies under TE incidence. (h) The structure diagram of the constructed information encoder.

Figure 8. (a–c) The near-field intensity distributions of the designed information encoder at TM incidence, when the refractive index of the surrounding medium is 1, 1.1 and 1.3.

Table 1

Comparison of sensor performance in recent years.

Reference	Freguency (THz)	Material	S (GHz/RIU)	FOM
[37]	0.5–1.1	Silicon	77	11.1
[38]	0.2–1.2	Gold	105	7.51
[39]	0.15–0.9	Aluminium	139.2	-
[40]	0.4–1.6	Graphene	177.7	59.3
[41]	555–568	Silicon nitride	-	163.5
[42]	0.2–2.1	LiTaO3	438	515
This work	0.5–0.7	Silicon	147	93

References

1. O′Brien, S.; Pendry, J.B. Photonic band-gap effects and magnetic activity in dielectric composites. J. Phys.-Condens. Mat.; 2002; 14, pp. 4035-4044. [DOI: https://dx.doi.org/10.1088/0953-8984/14/15/317]

2. Yang, Q.L.; Chen, X.Y.; Xu, Q.; Tian, C.X.; Xu, Y.H.; Cong, L.Q.; Zhang, X.Q.; Li, Y.F.; Zhang, C.H.; Zhang, X.X. et al. Broadband terahertz rotator with an all-dielectric metasurface. Photonics Res.; 2018; 6, pp. 1056-1061. [DOI: https://dx.doi.org/10.1364/PRJ.6.001056]

3. Deng, C.Z.; Ho, Y.L.; Clark, J.K.; Yatsui, T.; Delaunay, J.J. Light Switching with a Metal-Free Chiral-Sensitive Metasurface at Telecommunication Wavelengths. ACS Photonics; 2020; 7, pp. 2915-2922. [DOI: https://dx.doi.org/10.1021/acsphotonics.0c01377]

4. Li, J.; Yue, Z.; Li, J.; Zheng, C.; Wang, S.; Li, M.; Zhang, Y.T.; Zhang, Y.; Yao, J. Diverse terahertz wavefront manipulations empowered by the spatially interleaved metasurfaces. Sci. China-Inf. Sci.; 2023; 66, 132301. [DOI: https://dx.doi.org/10.1007/s11432-022-3499-4]

5. Yue, Z.; Zheng, C.L.; Li, J.; Li, J.T.; Liu, J.Y.; Wang, G.C.; Chen, M.Y.; Xu, H.; Tan, Q.; Zhang, H.J. et al. A dual band spin-selective transmission metasurface and its wavefront manipulation. Nanoscale; 2021; 13, pp. 10898-10905. [DOI: https://dx.doi.org/10.1039/D1NR02624K]

6. Li, J.; Yue, Z.; Li, J.; Zheng, C.; Liu, J.; Yang, F.; Li, H.; Zhang, Y.; Zhang, Y.; Yao, J. Wavefront-controllable all-silicon terahertz meta-polarizer. Sci. China Mater.; 2023; 66, pp. 300-308. [DOI: https://dx.doi.org/10.1007/s40843-022-2126-0]

7. Li, J.T.; Wang, G.C.; Yue, Z.; Liu, J.Y.; Li, J.; Zheng, C.L.; Zhang, Y.T.; Zhang, Y.; Yao, J.Q. Dynamic phase assembled terahertz metalens for reversible conversion between linear polarization and arbitrary circular polarization. Opto-Electron. Adv.; 2022; 5, 210062-1. [DOI: https://dx.doi.org/10.29026/oea.2022.210062]

8. Yue, Z.; Li, J.; Li, J.; Zheng, C.; Liu, J.; Wang, G.; Xu, H.; Chen, M.; Zhang, Y.; Zhang, Y. et al. Terahertz metasurface zone plates with arbitrary polarizations to a fixed polarization conversion. Opto-Electron. Sci.; 2022; 1, 210014. [DOI: https://dx.doi.org/10.29026/oes.2022.210014]

9. Yue, Z.; Li, J.T.; Liu, J.Y.; Li, J.; Zheng, C.L.; Wang, G.C.; Xu, H.; Chen, M.Y.; Zhang, Y.T.; Zhang, Y. et al. Versatile Polarization Conversion and Wavefront Shaping Based on Fully Phase-Modulated Metasurface with Complex Amplitude Modulation. Adv. Opt. Mater.; 2022; 10, 2200733. [DOI: https://dx.doi.org/10.1002/adom.202200733]

10. Wang, Q.; Plum, E.; Yang, Q.L.; Zhang, X.Q.; Xu, Q.; Xu, Y.H.; Han, J.G.; Zhang, W.L. Reflective chiral meta-holography: Multiplexing holograms for circularly polarized waves. Light-Sci. Appl.; 2018; 7, 25. [DOI: https://dx.doi.org/10.1038/s41377-018-0019-8]

11. Zheng, G.X.; Zhou, N.; Deng, L.G.; Li, G.F.; Tao, J.; Li, Z.L. Full-space metasurface holograms in the visible range. Opt. Express; 2021; 29, pp. 2920-2930. [DOI: https://dx.doi.org/10.1364/OE.417202] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33726478]

12. Zhang, X.; Li, X.; Zhou, H.Q.; Wei, Q.S.; Geng, G.Z.; Li, J.J.; Li, X.W.; Wang, Y.T.; Huang, L.L. Multifocal Plane Display Based on Dual Polarity Stereoscopic Metasurface. Adv. Funct. Mater.; 2022; 32, 2209460. [DOI: https://dx.doi.org/10.1002/adfm.202209460]

13. Slobozhanyuk, A.P.; Poddubny, A.N.; Raaijmakers, A.J.E.; van den Berg, C.A.T.; Kozachenko, A.V.; Dubrovina, I.A.; Melchakova, I.V.; Kivshar, Y.S.; Belov, P.A. Enhancement of Magnetic Resonance Imaging with Metasurfaces. Adv. Mater.; 2016; 28, pp. 1832-1838. [DOI: https://dx.doi.org/10.1002/adma.201504270] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26754827]

14. Zhao, F.; Shen, Z.C.; Wang, D.C.; Xu, B.J.; Chen, X.N.; Yang, Y.M. Synthetic aperture metalens. Photonics Res.; 2021; 9, pp. 2388-2397. [DOI: https://dx.doi.org/10.1364/PRJ.440185]

15. Zhou, J.X.; Zhao, J.X.; Wu, Q.Y.; Chen, C.F.; Lei, M.; Chen, G.H.; Tian, F.L.; Liu, Z.W. Nonlinear Computational Edge Detection Metalens. Adv. Funct. Mater.; 2022; 32, 2204734. [DOI: https://dx.doi.org/10.1002/adfm.202204734]

16. Li, J.; Li, J.T.; Yue, Z.; Zheng, C.L.; Wang, G.C.; Liu, J.Y.; Xu, H.; Song, C.Y.; Yang, F.; Li, H. et al. Structured Vector Field Manipulation of Terahertz Wave along the Propagation Direction Based on Dielectric Metasurfaces. Laser Photonics Rev.; 2022; 16, 2200325. [DOI: https://dx.doi.org/10.1002/lpor.202200325]

17. Yue, Z.; Liu, J.Y.; Li, J.T.; Li, J.; Zheng, C.L.; Chen, M.Y.; Xu, H.; Zhang, H.J.; Yang, F.; Zhang, Y.T. et al. Vector beam generation based on spin-decoupling metasurface zone plate. Appl. Phys. Lett.; 2022; 120, 191704. [DOI: https://dx.doi.org/10.1063/5.0093708]

18. Zheng, C.L.; Li, J.; Liu, J.Y.; Li, J.T.; Yue, Z.; Li, H.; Yang, F.; Zhang, Y.T.; Zhang, Y.; Yao, J.Q. Creating Longitudinally Varying Vector Vortex Beams with an All-Dielectric Metasurface. Laser Photonics Rev.; 2022; 16, 2200236. [DOI: https://dx.doi.org/10.1002/lpor.202200236]

19. Wang, Q.H.; Ni, P.N.; Xie, Y.Y.; Kan, Q.; Chen, P.P.; Fu, P.; Deng, J.; Jin, T.L.; Chen, H.D.; Lee, H.W.H. et al. On-Chip Generation of Structured Light Based on Metasurface Optoelectronic Integration. Laser Photonics Rev.; 2021; 15, 2000385. [DOI: https://dx.doi.org/10.1002/lpor.202000385]

20. Chen, Y.; Zhao, C.; Zhang, Y.Z.; Qiu, C.W. Integrated Molar Chiral Sensing Based on High-Q Metasurface. Nano Lett.; 2020; 20, pp. 8696-8703. [DOI: https://dx.doi.org/10.1021/acs.nanolett.0c03506]

21. Hsiao, H.H.; Hsu, Y.C.; Liu, A.Y.; Hsieh, J.C.; Lin, Y.H. Integrated Molar Chiral Sensing Based on High-Q Metasurface. Ultrasensitive Refractive Index Sensing Based on the Quasi-Bound States in the Continuum of All-Dielectric Metasurfaces. Adv. Opt. Mater.; 2022; 10, 2200812. [DOI: https://dx.doi.org/10.1002/adom.202200812]

22. Maksimov, D.N.; Gerasimov, V.S.; Romano, S.; Polyutov, S.P. Refractive index sensing with optical bound states in the continuum. Opt. Express; 2020; 28, pp. 38907-38916. [DOI: https://dx.doi.org/10.1364/OE.411749] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33379449]

23. Huang, C.; Zhang, C.; Xiao, S.M.; Wang, Y.H.; Fan, Y.B.; Liu, Y.L.; Zhang, N.; Qu, G.Y.; Ji, H.J.; Han, J.C. et al. Ultrafast control of vortex microlasers. Science; 2020; 367, 1018. [DOI: https://dx.doi.org/10.1126/science.aba4597]

24. Bernhardt, N.; Koshelev, K.; White, S.J.U.; Meng, K.W.C.; Froch, J.E.; Kim, S.; Tran, T.T.; Choi, D.Y.; Kivshar, Y.; Solntsev, A.S. Quasi-BIC Resonant Enhancement of Second-Harmonic Generation in WS2 Monolayers. Nano Lett.; 2020; 20, pp. 5309-5314. [DOI: https://dx.doi.org/10.1021/acs.nanolett.0c01603] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32530635]

25. Li, J.T.; Li, J.; Zheng, C.L.; Yue, Z.; Wang, S.L.; Li, M.Y.; Zhao, H.L.; Zhang, Y.T.; Yao, J.Q. Free switch between bound states in the continuum (BIC) and quasi-BIC supported by graphene-metal terahertz metasurfaces. Carbon; 2021; 182, pp. 506-515. [DOI: https://dx.doi.org/10.1016/j.carbon.2021.06.037]

26. Li, J.T.; Yue, Z.; Li, J.; Zheng, C.L.; Zhang, Y.T.; Yao, J.Q. Ultra-narrowband terahertz circular dichroism driven by planar metasurface supporting chiral quasi bound states in continuum. Opt. Laser Technol.; 2023; 161, 109173. [DOI: https://dx.doi.org/10.1016/j.optlastec.2023.109173]

27. Hsu, C.W.; Zhen, B.; Stone, A.D.; Joannopoulos, J.D.; Soljacic, M. Bound states in the continuum. Nat. Rev. Mater.; 2016; 1, 16048. [DOI: https://dx.doi.org/10.1038/natrevmats.2016.48]

28. Cong, L.Q.; Singh, R. Symmetry-Protected Dual Bound States in the Continuum in Metamaterials. Adv. Opt. Mater.; 2019; 7, 1900383. [DOI: https://dx.doi.org/10.1002/adom.201900383]

29. Wang, J.; Kuhne, J.; Karamanos, T.; Rockstuhl, C.; Maier, S.A.; Tittl, A. All-Dielectric Crescent Metasurface Sensor Driven by Bound States in the Continuum. Adv. Funct. Mater.; 2021; 31, 2104652. [DOI: https://dx.doi.org/10.1002/adfm.202104652]

30. Hu, Y.Z.; Tong, M.Y.; Hu, S.Y.; He, W.B.; Cheng, X.A.; Jiang, T. Spatiotemporal Lineshape Tailoring in BIC-Mediated Reconfigurable Metamaterials. Adv. Funct. Mater.; 2022; 32, 2203680. [DOI: https://dx.doi.org/10.1002/adfm.202203680]

31. Campione, S.; Liu, S.; Basilio, L.I.; Warne, L.K.; Langston, W.L.; Luk, T.S.; Wendt, J.R.; Reno, J.L.; Keeler, G.A.; Brener, I. et al. Broken Symmetry Dielectric Resonators for High Quality Factor Fano Metasurfaces. ACS Photonics; 2016; 3, pp. 2362-2367. [DOI: https://dx.doi.org/10.1021/acsphotonics.6b00556]

32. Koshelev, K.; Lepeshov, S.; Liu, M.K.; Bogdanov, A.; Kivshar, Y. Asymmetric Metasurfaces with High-Q Resonances Governed by Bound States in the Continuum. Phys. Rev. Lett.; 2018; 121, 193903. [DOI: https://dx.doi.org/10.1103/PhysRevLett.121.193903] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30468599]

33. Yesilkoy, F.; Arvelo, E.R.; Jahani, Y.; Liu, M.K.; Tittl, A.; Cevher, V.; Kivshar, Y.; Altug, H. Ultrasensitive hyperspectral imaging and biodetection enabled by dielectric metasurfaces. Nat. Photonics; 2019; 13, pp. 390-396. [DOI: https://dx.doi.org/10.1038/s41566-019-0394-6]

34. Li, J.T.; Li, J.; Zheng, C.L.; Yue, Z.; Yang, D.Y.; Wang, S.L.; Li, M.Y.; Zhang, Y.T.; Yao, J.Q. Spectral amplitude modulation and dynamic near-field displaying of all-silicon terahertz metasurfaces supporting bound states in the continuum. Appl. Phys. Lett.; 2021; 119, 241105. [DOI: https://dx.doi.org/10.1063/5.0067937]

35. Wang, Y.L.; Han, Z.H.; Du, Y.; Qin, J.Y. Ultrasensitive terahertz sensing with high-Q toroidal dipole resonance governed by bound states in the continuum in all-dielectric metasurface. Nanophotonics; 2021; 10, pp. 1295-1307. [DOI: https://dx.doi.org/10.1515/nanoph-2020-0582]

36. Sherry, L.J.; Chang, S.-H.; Schatz, G.C.; Van Duyne, R.P.; Wiley, B.J.; Xia, Y. Localized surface plasmon resonance spectroscopy of single silver nanocubes. Nano Lett.; 2005; 5, pp. 2034-2038. [DOI: https://dx.doi.org/10.1021/nl0515753]

37. Zhong, Y.; Du, L.; Liu, Q.; Zhu, L.; Meng, K.; Zou, Y.; Zhang, B. Ultrasensitive specific sensor based on all-dielectric metasurfaces in the terahertz range. RSC Adv.; 2020; 10, pp. 33018-33025. [DOI: https://dx.doi.org/10.1039/D0RA06463G]

38. Xie, Q.; Dong, G.-X.; Wang, B.X.; Huang, W.Q. High-Q fano resonance in terahertz frequency based on an asymmetric metamaterial resonator. Nanoscale Res. Lett.; 2018; 13, 294. [DOI: https://dx.doi.org/10.1186/s11671-018-2677-0]

39. Yahiaoui, R.; Tan, S.; Cong, L.; Singh, R.; Yan, F.; Zhang, W. Multispectral terahertz sensing with highly flexible ultrathin metamaterial absorber. J. Appl. Phys.; 2015; 118, 083103. [DOI: https://dx.doi.org/10.1063/1.4929449]

40. He, X.; Zhang, Q.; Lu, G.; Ying, G.; Wu, F.; Jiang, J. Tunable ultrasensitive terahertz sensor based on complementary graphene metamaterials. RSC Adv.; 2016; 6, pp. 52212-52218. [DOI: https://dx.doi.org/10.1039/C5RA21974D]

41. Schiattarella, C.; Sanità, G.; Alulema, B.G.; Lanzio, V.; Cabrini, S.; Lamberti, A.; Rendina, I.; Mocella, V.; Zito, G.; Romano, S. High-Q photonic aptasensor based on avoided crossing bound states in the continuum and trace detection of ochratoxin A. Biosens. Bioelectron. X; 2022; 12, 100262. [DOI: https://dx.doi.org/10.1016/j.biosx.2022.100262]

42. Chen, X.; Fan, W. Ultrahigh-Q toroidal dipole resonance in all-dielectric metamaterials for terahertz sensing. Opt. Lett.; 2019; 44, pp. 5876-5879. [DOI: https://dx.doi.org/10.1364/OL.44.005876] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31774802]

43. Hu, S.; Liu, D.; Yang, H.L.; Wang, H.X.; Wang, Y. Staggered H-shaped metamaterial based on electromagnetically induced transparency effect and its refractive index sensing performance. Opt. Commun.; 2019; 450, pp. 202-207. [DOI: https://dx.doi.org/10.1016/j.optcom.2019.05.038]