Investigation of Multiple High Quality-Factor Fano Resonances in Asymmetric Nanopillar Arrays for Optical Sensing

: A novel asymmetric all-dielectric metasurface supporting multiple Fano resonances with high quality-factor through the excitation of quasi-bound states in the continuum is theoretically investigated. It is demonstrated that two resonances in the near-infrared wavelength are excited by the symmetry-protected bound state in the continuum, which can be transformed into the electric dipole and the toroidal dipole quasi-BIC resonance with high quality-factor by breaking the symmetry of metasurface. Moreover, the sensing properties based on different liquid refractive indexes are researched theoretically. The results show that the maximum quality-factor of the Fano resonance peak is 8422, and the sensitivity can reach 402 nm/RIU, with a maximum figure of merit of 2400 RIU − 1 . This research is believed to further promote the development of optical sensing and nonlinear optics.


Introduction
As a two-dimensional metamaterial, metasurfaces can be artificially engineered into subwavelength structures to obtain electromagnetic properties not found in natural materials [1,2].By optimizing the structural parameters of the metasurface, high quality-factor (Q-factor) Fano resonances are achieved, which have important roles in the fields of optical refractive index sensing [3,4], optical switches [5], nonlinear optics [6,7], and slow light [8].The high Q-factor Fano resonance not only enhances the light-matter interactions but also has a high sensitivity to tiny environmental changes [9,10].
High Q-factor Fano resonances are obtained in many ways.Mohammadi et al. presented a photonic crystal ring resonator including three different photonic crystal sensors with Q-factors up to 5365 [11].Zhou et al. designed a guided-mode resonance sensor with a shallow subwavelength structure with a Q-factor up to 8000 [12].All-dielectric metasurfaces offer greater freedom and design flexibility, allowing for more diversified functionality [13].Different optical functions can be achieved by precisely designing nanostructures that allow precise modulation of several parameters of light, such as polarization, wavelength, phase, and amplitude [14].At the same time, the all-dielectric metasurface is capable of wide-bandwidth functionality, allowing simultaneous operation of resonances in multiple frequency ranges.In addition, all-dielectric metasurface sensors have better tunability.Wavelength selectivity of the sensor is achieved by adjusting the geometry of the structure, material properties or electromagnetic excitation.
In recent years, silicon nitride has been an important member of the silicon-based family for making resonators [15,16].Yang et al. designed a refractive index sensor based on bilayer silicon nitride with a Q-factor of 7605 [17].Guo et al. proposed an all-dielectric metasurface based on a silicon nitride substrate, where the Q-factor of the resonance peak can reach 54757 [18].In addition, silicon nitride is widely used in a variety of resonators (drum resonators [19], micro-ring resonators [20], waveguide-coupled resonators [21,22]), oscillators [23], and semiconductor lasers [24]).Silicon nitride-based optics can reach Qfactors of the order of 10 6 in most cases, with the highest reaching 2.6 × 10 8 .Although silicon nitride resonators have extremely high Q-factors, all-dielectric metasurface sensors composed of silicon and silicon dioxide are still very necessary because silicon and silicon dioxide are easily manufactured using current manufacturing techniques [25].The alldielectric metasurface composed of silicon and silica can excite multiple Fano resonance peaks simultaneously.More resonant peaks are better suited for multi-channel sensing applications [26].
The all-dielectric metasurface is believed to be a great way to obtain high Q-factor Fano resonance due to its greater freedom, design flexibility, and compatibility with complementary metal oxide semiconductor (CMOS) processes [27].The excitation of Fano resonance in metasurfaces has evolved from a single resonance peak to multiple resonance peaks with the development of sensing technology [28,29].Multi-Fano resonances have potential applications in the fields of multi-wavelength sensing and multi-channel biosensing [30].Yang et al. proposed an all-dielectric metasurfaces consisting of an array of square nanopores, which excited four Fano resonances by introducing nanopore asymmetry [31].Wang et al. investigated an all-dielectric metasurface consisting of an asymmetric nanocylindrical dimer with five Fano resonance peaks that were excited [28].
An efficient way for generating high Q-factor Fano resonances in all-dielectric metasurfaces is through the utilization of bound states in the continuum (BIC) [32,33].The BIC is a non-radiative state characterized by a resonance frequency embedded in the continuum spectrum of radiative modes in the surrounding space.The coupling of the resonant mode to the surrounding external channels is prohibited when symmetry is maintained, and this type of BIC is referred to as a symmetry-protected BIC [34,35], which is considered to be an ideal BIC arising from an infinitely high Q-factor and extremely narrow resonance width [36,37].However, this ideal BIC cannot be observed because of the absence of farfield radiation and the inability to be excited by incident waves.It is necessary to create a radiation channel to the outside by breaking the symmetry of the result or by oblique incidence of the light source [38].Through this method, the symmetry-protected BIC is switched to the quasi-BIC mode so that the Q-factor of the resonance becomes finite [39,40].
In this paper, an all-dielectric metasurface that can excite multiple Fano resonance peaks in the near-infrared spectral region is designed and simulated.The metasurface is consisting of five rectangular silicon blocks periodically placed on a silica-based substrate.Changing the length of two of the rectangular silicon blocks breaks the symmetry of the structure, and thus, four Fano resonance peaks are excited.The Q-factors of the four Fano resonance peaks all reach 10 3 , and the highest can reach 8422, the sensitivity is up to 402 nm/RIU, and the figure of merit (FOM) can reach up to 2400 RIU −1 .The preparation process of the structure is investigated as a way to demonstrate the realizability of the structure.It is shown that the structure has potential applications in the field of refractive index sensors.

Structure Design
The schematics of the designed metasurface structure are shown in Figure 1. Figure 1a shows the array layout of the structure, which is periodically distributed on a substrate with the material SiO 2 .The structure of a single cell is shown in Figure 1b, where the labeled P x and P y denote the period in the x and y directions, respectively, and h is the height of the silicon block.Figure 1c shows the top view of the structure, the length of the two rectangular blocks on the left (L 1 ) is 280 nm, the length of the two rectangular blocks on the right (L 2 ) is 170 nm, the length of the rectangular blocks in the center (L 3 ) is 660 nm, the width of the rectangular blocks (w) is 150 nm, and the distance between the rectangular blocks (g) is 100 nm.The structure is asymmetrical along the y-axis since L 1 ̸ = L 2 .δ = L 1 − L 2 as an asymmetry parameter is used to explore the effect of asymmetry on the excitation of the Fano resonance peaks.
with the material SiO2.The structure of a single cell is shown in Figure 1b, where the labeled P and P denote the period in the x and y directions, respectively, and h is the height of the silicon block.Figure 1c shows the top view of the structure, the length of the two rectangular blocks on the left (L ) is 280 nm, the length of the two rectangular blocks on the right (L ) is 170 nm, the length of the rectangular blocks in the center (L ) is 660 nm, the width of the rectangular blocks (w) is 150 nm, and the distance between the rectangular blocks (g) is 100 nm.The structure is asymmetrical along the y-axis since L ≠ L .δ = L − L as an asymmetry parameter is used to explore the effect of asymmetry on the excitation of the Fano resonance peaks.The system of Maxwell's equations is transformed into finite difference equations by the finite difference in time domain (FDTD) method, enabling the simulations to be conducted.With this method, the distribution information of the electromagnetic field can be calculated when the boundary conditions and initial values are known, which greatly reduces the difficulty of the calculation.In the simulation, it is necessary to set the appropriate mesh size.The resolution of the mesh needs to be detailed enough to capture the details and rapidly changing features that appear in the simulation.Excessive mesh size leads to a decrease in the accuracy of the simulation results, as tiny structures and details cannot be adequately represented.If the mesh size is too small, it will lead to an increase in the computational effort, thus decreasing the computational efficiency.Therefore, the mesh was set to 8 nm × 8 nm × 8 nm during the simulation.In order to simulate the behavior of a system in three dimensions, boundary conditions need to be set.In this case, the X and Y directions are set to periodic boundary conditions, which means the simulation accurately represents the behavior of the system across multiple repetitions.On the other hand, the Z direction is set to Perfectly Matched Layer (PML).This is a type of boundary condition that effectively absorbs outgoing waves in the system, preventing reflections and ensuring accurate results.By combining periodic boundary conditions in the X and Y directions with a PML boundary condition in the Z direction, a comprehensive and reliable simulation can be achieved.The y-polarized plane wave perpendicular to the z-axis as a light source is used in the simulations.The plane wave is used in the simulation process.Plane waves propagate through space at a constant velocity without changing shape or decaying.Optical parameters of silicon and silicon dioxide used in the simulations are referenced in the Palik manual.

Results and Discussion
The transmission spectrum of the structure in the liquid refractive index of 1.35 state is simulated, while the background refractive index is set to n = 1.35 .When L = The system of Maxwell's equations is transformed into finite difference equations by the finite difference in time domain (FDTD) method, enabling the simulations to be conducted.With this method, the distribution information of the electromagnetic field can be calculated when the boundary conditions and initial values are known, which greatly reduces the difficulty of the calculation.In the simulation, it is necessary to set the appropriate mesh size.The resolution of the mesh needs to be detailed enough to capture the details and rapidly changing features that appear in the simulation.Excessive mesh size leads to a decrease in the accuracy of the simulation results, as tiny structures and details cannot be adequately represented.If the mesh size is too small, it will lead to an increase in the computational effort, thus decreasing the computational efficiency.Therefore, the mesh was set to 8 nm × 8 nm × 8 nm during the simulation.In order to simulate the behavior of a system in three dimensions, boundary conditions need to be set.In this case, the X and Y directions are set to periodic boundary conditions, which means the simulation accurately represents the behavior of the system across multiple repetitions.On the other hand, the Z direction is set to Perfectly Matched Layer (PML).This is a type of boundary condition that effectively absorbs outgoing waves in the system, preventing reflections and ensuring accurate results.By combining periodic boundary conditions in the X and Y directions with a PML boundary condition in the Z direction, a comprehensive and reliable simulation can be achieved.The y-polarized plane wave perpendicular to the z-axis as a light source is used in the simulations.The plane wave is used in the simulation process.Plane waves propagate through space at a constant velocity without changing shape or decaying.Optical parameters of silicon and silicon dioxide used in the simulations are referenced in the Palik manual.

Results and Discussion
The transmission spectrum of the structure in the liquid refractive index of 1.35 state is simulated, while the background refractive index is set to n = 1.35.When L 1 = L 2 (δ = 0 nm), the transmission spectrum is obtained as shown by the black solid line in Figure 2a.Two Fano resonance peaks marked as R1 and R2 are observed at the wavelength of 1123 nm and 1163 nm, respectively.The asymmetric metasurface transmission spectrum is shown in Figure 2b.Two new Fano resonance peaks, labeled R3 and R4, are innovatively excited at 1202 nm and 1221 nm when the length of L 2 is changed to 170 nm (δ = 110 nm).As the symmetry of the metasurface is broken, radiation channels are established, allowing energy to leak outward, resulting in a symmetry-protected BIC converting to a quasi-BIC.
of 1123 nm and 1163 nm, respectively.The asymmetric metasurface transmission spectrum is shown in Figure 2b.Two new Fano resonance peaks, labeled R3 and R4, are innovatively excited at 1202 nm and 1221 nm when the length of L is changed to 170 nm (δ = 110 nm ).As the symmetry of the metasurface is broken, radiation channels are established, allowing energy to leak outward, resulting in a symmetry-protected BIC converting to a quasi-BIC.In addition, the excited Fano resonance peaks of the periodic all-dielectric metasurface can be analyzed and fitted by the Fano model with Equation (1) [36]: where  ,  , and  are constant real numbers,  is the resonant frequency, and  is the damping loss.The Fano curves fitted by Equation ( 1) are shown as red and green dashed lines in Figure 2a.The fitting results are in general agreement with the simulation results.Equation ( 2) is extracted by Equation ( 1), and the Q-factor of the Fano resonance can be obtained by Equation (2).
The The relationship between the Q-factor and the degree of asymmetry is discussed.α = ∆/ as the degree of asymmetry, where  is the area of the silicon block when symmetric and ∆ is the area of the reduced silicon block when asymmetric.As ∆S decreases, the Qfactor is larger.After observing several sets of data, the Q-factor of the excited Fano resonance is found to be related to the degree of asymmetry as  ∝  , as shown in Figure 3a.In addition, the excited Fano resonance peaks of the periodic all-dielectric metasurface can be analyzed and fitted by the Fano model with Equation (1) [36]: where a 1 , a 2 , and b are constant real numbers, ω 0 is the resonant frequency, and γ is the damping loss.The Fano curves fitted by Equation ( 1) are shown as red and green dashed lines in Figure 2a.The fitting results are in general agreement with the simulation results.Equation ( 2) is extracted by Equation ( 1), and the Q-factor of the Fano resonance can be obtained by Equation (2).
The relationship between the Q-factor and the degree of asymmetry is discussed.α = ∆S/S as the degree of asymmetry, where S is the area of the silicon block when symmetric and ∆S is the area of the reduced silicon block when asymmetric.As ∆S decreases, the Q-factor is larger.After observing several sets of data, the Q-factor of the excited Fano resonance is found to be related to the degree of asymmetry as Q ∝ α −2 , as shown in Figure 3a.
In order to further investigate the nature of the multiple resonance peaks generated by this metasurface, the electromagnetic field distributions of the resonance peaks under asymmetry are simulated, and their near-field distribution is shown in Figure 4.For the resonance peak R1, the magnetic field forms two clockwise loops in the x-z plane, and the electric field is linear in the x-y plane, which is a typical electric dipole (ED) feature.The resonance mode of R3 is also considered to be ED, but its direction is opposite to that of R1.The resonance mode of R2 can be seen in the figure as magnetic dipoles (MDs), in which the current forms three loops with the same direction in the y-z plane and the magnetic field is linearly polarized in the x-y plane.Two current loops in opposite directions can be observed in the x-y plane of R4, and the resulting magnetic field creates a large loop in the x-z plane, which is typical of toroidal dipole (TD) resonance.In order to further investigate the nature of the multiple resonance peaks generated by this metasurface, the electromagnetic field distributions of the resonance peaks under asymmetry are simulated, and their near-field distribution is shown in Figure 4.For the resonance peak R1, the magnetic field forms two clockwise loops in the x-z plane, and the electric field is linear in the x-y plane, which is a typical electric dipole (ED) feature.The resonance mode of R3 is also considered to be ED, but its direction is opposite to that of R1.The resonance mode of R2 can be seen in the figure as magnetic dipoles (MDs), in which the current forms three loops with the same direction in the y-z plane and the magnetic field is linearly polarized in the x-y plane.Two current loops in opposite directions can be observed in the x-y plane of R4, and the resulting magnetic field creates a large loop in the x-z plane, which is typical of toroidal dipole (TD) resonance.The characteristics of the transmission spectrum for different geometrical parameters are also discussed to obtain the optimal parameter, and the simulation results are shown in Figure 5. Figure 5a shows the transmission spectrums of the silicon block with a high between 180  In order to further investigate the nature of the multiple resonance peaks generated by this metasurface, the electromagnetic field distributions of the resonance peaks under asymmetry are simulated, and their near-field distribution is shown in Figure 4.For the resonance peak R1, the magnetic field forms two clockwise loops in the x-z plane, and the electric field is linear in the x-y plane, which is a typical electric dipole (ED) feature.The resonance mode of R3 is also considered to be ED, but its direction is opposite to that of R1.The resonance mode of R2 can be seen in the figure as magnetic dipoles (MDs), in which the current forms three loops with the same direction in the y-z plane and the magnetic field is linearly polarized in the x-y plane.Two current loops in opposite directions can be observed in the x-y plane of R4, and the resulting magnetic field creates a large loop in the x-z plane, which is typical of toroidal dipole (TD) resonance.The characteristics of the transmission spectrum for different geometrical parameters are also discussed to obtain the optimal parameter, and the simulation results are shown in Figure 5. Figure 5a shows the transmission spectrums of the silicon block with a high between 180 The characteristics of the transmission spectrum for different geometrical parameters are also discussed to obtain the optimal parameter, and the simulation results are shown in Figure 5. Figure 5a shows the transmission spectrums of the silicon block with a high between 180 nm and 220 nm, and it is clear that the resonance peaks are sensitive because the four resonance peaks red-shift as h increases, especially R1.The transmission spectrums of the silicon block with a wide w between 130 nm and 170 nm are presented in Figure 5b, where all the four resonance peaks are red-shifted with w increasing.In addition, the linewidth of R1 decreases as w increases, and the modulation depth of R4 becomes deeper as w increases.The transmission spectrum corresponding to the change in length of the intermediate silicon block (L 3 ) is displayed in Figure 5c.The L 3 produces little effect on the transmission spectrum of the metasurface.Figure 5d parades the transmission spectra at different periods.The Fano resonance peaks show different degrees of red-shift as the period increases.However, compared with h and w, this red-shift is slight, while the modulation depth of R4 undergoes a little change.
as w increases.The transmission spectrum corresponding to the change in length of the intermediate silicon block (L ) is displayed in Figure 5c.The L produces little effect on the transmission spectrum of the metasurface.Figure 5d parades the transmission spectra at different periods.The Fano resonance peaks show different degrees of red-shift as the period increases.However, compared with h and w , this red-shift is slight, while the modulation depth of R4 undergoes a little change.In order to investigate the influence of light absorption and scattering caused by liquid and materials on the transmission spectrum of the structure, different extinction coefficients (k) are set for symmetric and asymmetric cases, respectively.The simulation results are shown in Figure 6.The linewidth of the resonance peak increases with increasing extinction coefficient; however, the modulation depth and Q-factor decrease.The sensitivity of the resonance peaks to the extinction coefficient varies, with R4 being the most sensitive and R1 the least sensitive.What this result presents is that the same extinction coefficient will have different effects on the Fano resonance formed by different resonant modes.These simulations demonstrated that the presence of a large intrinsic absorption in the liquid and materials used would cause the resonance to attenuate or even disappear.In order to investigate the influence of light absorption and scattering caused by liquid and materials on the transmission spectrum of the structure, different extinction coefficients (k) are set for symmetric and asymmetric cases, respectively.The simulation results are shown in Figure 6.The linewidth of the resonance peak increases with increasing extinction coefficient; however, the modulation depth and Q-factor decrease.The sensitivity of the resonance peaks to the extinction coefficient varies, with R4 being the most sensitive and R1 the least sensitive.What this result presents is that the same extinction coefficient will have different effects on the Fano resonance formed by different resonant modes.These simulations demonstrated that the presence of a large intrinsic absorption in the liquid and materials used would cause the resonance to attenuate or even disappear.The proposed metastructure is shown to be investigated as a refractive index sensor due to its advantages of excitation of high Q-factor Fano resonance and local field enhancement.The potential application of metasurfaces in the field of refractive index sensors is verified by analyzing the changing patterns of transmission spectra at different ambient refractive indexes.Figure 7 show the transmission spectra of R1, R2, R3, and R4 at ambient refractive index (n) in steps of 0.01 from 1.33 to 1.37, respectively.As can be seen from Figure 7, the Fano resonance peaks are significantly red-shifted despite the little change in refractive index.The sensing performance of the metasurface is believed to be determined by two metrics, the sensitivity S and the FOM.The sensitivity is calculated by The proposed metastructure is shown to be investigated as a refractive index sensor due to its advantages of excitation of high Q-factor Fano resonance and local field enhancement.The potential application of metasurfaces in the field of refractive index sensors is verified by analyzing the changing patterns of transmission spectra at different ambient refractive indexes.Figure 7 show the transmission spectra of R1, R2, R3, and R4 at ambient refractive index (n) in steps of 0.01 from 1.33 to 1.37, respectively.As can be seen from Figure 7, the Fano resonance peaks are significantly red-shifted despite the little change in refractive index.The sensing performance of the metasurface is believed to be determined by two metrics, the sensitivity S and the FOM.The sensitivity is calculated by Equation (3) [41]: in which ∆λ is the offset of the resonance peak and ∆n is the difference in refractive index.
Refractive index unit (RIU) represents the ratio of the speed of light relative to the speed of light in a vacuum as it enters the medium from a vacuum.The sensitivities of the Fano resonance peaks are 129, 333, 403, 348 by calculation.Figure 7b,d show the wavelengths corresponding to the four resonance peaks at different refractive indexes, and the slope of the curve corresponding to each resonance peak is the sensitivity of that resonance peak.The FOM is calculated by Equation ( 4) [42]: in which S is the sensitivity and FWHM is the full width at half peak of the resonance peak.
The FOM corresponding to the four resonance peaks is calculated to be 123, 701, 448, and 2400, respectively.In addition, the proposed metastructure is compared with the existing work for number of resonance peaks, sensitivity, and FOM, as shown in Table 1.The results show that the metastructure proposed in this paper can generate more resonance peaks with higher sensitivity and FOM, which is more suitable for multichannel sensing applications.In addition, the proposed metastructure is compared with the existing work for number of resonance peaks, sensitivity, and FOM, as shown in Table 1.The results show that the metastructure proposed in this paper can generate more resonance peaks with higher sensitivity and FOM, which is more suitable for multichannel sensing applications.The proposed sensors are mainly based on silicon on insulator (SOI) fabrication.SOI device processes are compatible with standard silicon CMOS processes, which are easy and inexpensive to fabricate, making them ideal for large-scale semiconductor chip integration.The preparation process of the proposed metastructure is shown in Figure 8.The device preparation procedure consists of the following steps: cleaning of the SOI substrate, lowpressure chemical vapor deposition (LPCVD), spin coating of the resist, electron beam lithography (EBL), development, inductively coupled plasma (ICP) etching and removal of the resist.

Conclusions
In summary, we proposed and theoretically demonstrated a structure-asymmetric all-dielectric metasurface based on asymmetric nanopillar arrays.Four sharp and highspectral-contrast Fano resonances are excited with the best Q-factor of 8422.The results of the electromagnetic field distribution explicate that the resonances at 1107 nm and 1202 nm are governed by the ED mode, while the resonances at 1129 nm and 1221 nm are governed by the MDs mode and TD mode, respectively.In addition, the relationship between the different extinction coefficients (K) and transmission spectrum indicates that the resonances are attenuated or even disappear, occurring when there is a large intrinsic absorption in the liquid.Finally, the sensing characteristics based on different liquid refractive indices are analyzed to certify that the designed all-dielectric metasurface is employed as a refractive index sensor with the maximum sensitivity, and the FOM value is 402 nm/RIU and 2400 RIU −1 , respectively.This dielectric metasurface with multiple sharp resonances offers a good platform for multichannel sensing as well as optical modulation.

Figure 1 .
Figure 1.Schematic of the metasurface structure consisting of five rectangular columns.(a) Structural layout of the designed metasurface; (b) diagram of the structure of a single cell; (c) top view of a single cell.

Figure 1 .
Figure 1.Schematic of the metasurface structure consisting of five rectangular columns.(a) Structural layout of the designed metasurface; (b) diagram of the structure of a single cell; (c) top view of a single cell.

Figure 2 .
Figure 2. Simulated transmission spectra of the metasurface.(a) Transmission spectrum of symmetric structure and equation fitting plot; (b) transmission spectrum of asymmetric structure.

Figure 2 .
Figure 2. Simulated transmission spectra of the metasurface.(a) Transmission spectrum of symmetric structure and equation fitting plot; (b) transmission spectrum of asymmetric structure.

Figure 3 .
Figure 3. (a) Plot of the Q-factor of the excited Fano resonance versus the degree of asymmetry α; (b) schematic representation of ∆ and .

Figure 4 .
Figure 4. Near−field distribution of the electromagnetic field corresponding to the four Fano resonance peaks in the asymmetric case.The electric field vector distributions are shown as white arrows and the magnetic field vector distributions are shown as black arrows.

Figure 3 .
Figure 3. (a) Plot of the Q-factor of the excited Fano resonance versus the degree of asymmetry α; (b) schematic representation of ∆S and S.

Figure 3 .
Figure 3. (a) Plot of the Q-factor of the excited Fano resonance versus the degree of asymmetry α; (b) schematic representation of ∆ and .

Figure 4 .
Figure 4. Near−field distribution of the electromagnetic field corresponding to the four Fano resonance peaks in the asymmetric case.The electric field vector distributions are shown as white arrows and the magnetic field vector distributions are shown as black arrows.

Figure 4 .
Figure 4. Near−field distribution of the electromagnetic field corresponding to the four Fano resonance peaks in the asymmetric case.The electric field vector distributions are shown as white arrows and the magnetic field vector distributions are shown as black arrows.

Figure 5 .
Figure 5. Transmission spectra of the metasurface with varying geometrical parameters.(a) Height of the silicon block (h), (b) width of the silicon block (w), (c) length of the middle silicon block (L ), and (d) period of the metasurface (P).

Figure 5 .
Figure 5. Transmission spectra of the metasurface with varying geometrical parameters.(a) Height of the silicon block (h), (b) width of the silicon block (w), (c) length of the middle silicon block (L 3 ), and (d) period of the metasurface (P).

Figure 6 .
Figure 6.(a) Transmission spectra of symmetric structures with different extinction coefficients, (b)Transmission spectra of asymmetric structures with different extinction coefficients.

Figure 6 .
Figure 6.(a) Transmission spectra of symmetric structures with different extinction coefficients, (b)Transmission spectra of asymmetric structures with different extinction coefficients.

Figure 7 .
Figure 7. Transmission spectra of metasurfaces at different refractive indexes.(a) Transmission spectra of R1 and R2 at different refractive indexes; (b) wavelength plots corresponding to R1 and R2 at different refractive indexes; (c) transmission spectra of R3 and R4 at different refractive indexes; (d) wavelength plots corresponding to R3 and R4 at different refractive indexes.

Figure 7 .
Figure 7. Transmission spectra of metasurfaces at different refractive indexes.(a) Transmission spectra of R1 and R2 at different refractive indexes; (b) wavelength plots corresponding to R1 and R2 at different refractive indexes; (c) transmission spectra of R3 and R4 at different refractive indexes; (d) wavelength plots corresponding to R3 and R4 at different refractive indexes.

Figure 8 .
Figure 8.The fabrication process flow chart of the metastructure.

Table 1 .
The number of resonance peaks, sensitivity, and FOM of this sensor are compared with other types.

Table 1 .
The number of resonance peaks, sensitivity, and FOM of this sensor are compared with other types.
a not reported.