Terahertz Sensor Study Based on Spoof Surface Plasmon Polaritons

+e spoof surface plasmon polaritons (SSPPs) structure can be used as a sensor in THz region for the biosensing. +e accuracy of resonance and amplitude for sensor is very important for biosensing. +e momentum matching of SSPPs determines the resonance position and the gap distance determines the amplitude. For the biomolecular sensing, the sample is positioned between the prism base and the SSPPs structure. +e momentum matching condition at the current study does not consider the effect of sample refractive index and the resonance position has a significant error. Here the correction is made to the momentum matching condition which considers the effect of the sample refractive index. A comparative study of surface plasmon resonance (SPR) sensing performance based on frequency and angle variations shows that the sensing sensitivity for frequency region is superior to that of angle region; in the meanwhile, as an application of biosensors, we have detected different types of brain lesions in the frequency range. Furthermore, the reflection amplitude is related to gap size between the prism and SSPPs.+e relationship of gap size and reflection amplitude is studied. By using the relationship between gap size and reflection amplitude, the amplitudes at different frequencies or incident angles for different refractivities have the same reflection dips compared to the other published results. +e simulation is performed and the results proved the theory.


Introduction
In recent years, basic research and development based on surface plasmon polaritons (SPPs) structures and devices have attracted tremendous research enthusiasm, especially in the field of sensing. SPPs are collective oscillations of electron-electromagnetic coupling at a metal/dielectric interface [1]. eir unique properties make them have a wide range of practical applications, including the guidance and manipulation of light at the nanometer scale, biomacromolecules detection, enhanced transmission of light fields through subwavelength aperture, and the absorption and attenuation of electromagnetic (EM) waves [2][3][4][5][6]. e sensors based on surface plasmon resonance (SPR) are one of the most successful applications of SPR at optical frequencies [7]. Many important biological molecules have rotational, translational, and vibrational responses, which uniquely occur in the THz band, and THz radiation is nonionizing radiation, which enables terahertz sensing to have broad application prospects in the field of biology sensing [8]. However, the PEC properties of metallic materials in the THz region make the conventional SPPs unable to be supported at the metal/dielectric interface [9]. Fortunately, Pendry et al. demonstrated theoretically that this problem can be solved by the concept of THz structural surfaces [10,11], that is, through the modification of the metal surface, to support a terahertz frequency mimicking SPPs, called the spoof SPPs (SSPPs). e plasma frequency of SSPPs can be easily designed through the geometric parameters of the corrugated surfaces, allowing SSPPs to be used in sensing applications from gigahertz (GHz) to midinfrared frequencies.
We note that SSPPs sensing has been well discussed in the THz range [12][13][14][15][16][17][18][19][20][21][22][23]. For example, Ng et al. studied the relationship between spectral shift and refractive index of SPPs in 1D periodic grooves by terahertz time-domain spectroscopy (THz-TDS) [12]. Yao and Zhong proposed a high-order mode of spoof SPPs for periodic metal grooves for ultrasensitive terahertz refractive index sensing [14,15]. Zhang et al. analyzed the angle dependence of the resonance angle on the refractive index change [16]. Li et al.designed the THz wave switches based on K-R and Otto configurations using the principle of SPR refractive index sensing [20,21]. Recently, Niknam et al. and others have introduced a new type of double-corrugated ultrasensitive resonance sensor [22].
Generally, most of the SSPPs-based applications in the above literatures are the sensing of the refractive index of filled fluid; however, firstly, the momentum matching condition does not consider the influence of the sample refractive index, which will result in a significant error in the resonance position and reduce the sensitivity of the sensor. Obviously, it is highly undesirable for biomolecular detection. For example, the research work of Bray et al. [24] proved that there are an estimated 18.1 million new cancer cases (17 million excluding non-melanoma skin cancer) and 9.6 million cancer deaths (9.5 million excluding non-melanoma skin cancer). erefore, the accuracy of the resonance frequency of the SPR biosensor is crucial to improve the sensitivity of the sensor to detect these diseases earlier. Secondly, almost all designs are only used to study the frequency dependence or angle dependence, and there is no comparison between the performances of frequency sensing and angle sensing.
irdly, since the SSPPs wave decays exponentially in the direction perpendicular to the interface, the coupling gap distance will vary with the filling refractive index, but, in the existing literature [12,14,15], the sensing of different refractive indexes (RI) is carried out according to the same coupling gap. As the RI of the filling medium increases, the optimal coupling regime is broken, resulting in increasing amplitude and wider absorption peak of reflectivity, which limits the possibility of the RI sensing over a wider range.
In this paper, considering the effect of the sample refractive index (RI), we first modified the momentum matching condition, which will help to improve the accuracy and sensitivity of the "refractive index sensing." en, we verified the correctness of this modification by comparing the resonant positions of the dispersion curve obtained by CST Microwave Studio and the reflection spectrum calculated by COMSOL software. Furthermore, the optimal coupling gap for different refractive indexes is fitted, and the reflectivity is recalculated based on the fitting gap values. Without the shift of the resonance position, the amplitudes at different frequencies or incident angles for different RI have the same reflection dips. Finally, as an application platform for SSPPs biosensing, we carried out sensing detection on brain cancerous cells. e resonance frequency shifts of normal cells and cancer cells infiltrating the gap and groove were calculated by finite element method (FEM).

Theory and Configuration of SPR-Based Sensing
We used the conventional Otto prism-coupling configuration to excite SSPPs [25], as shown in Figure 1, which consists of a coupling prism and a metal grating with a onedimensional groove structure, where the period is p, the groove depth is h, the groove width is a, and the gap distance between the surface of the metal grating and the prism base is g. In order to excite SSPPs on the metal grating, a prism with a high permittivity ε p is used to increase the wave vector for phase matching. e sample to be measured with the permittivity ε d is filled in area I and area II. e merits of such a sensing configuration are that the metal groove can be regarded as a microfluidic channel, which significantly reduces the sample consumption and makes the terahertz sensing of biological molecules possible. Furthermore, since the sample is in a relatively closed channel, the disturbance is reduced to a certain extent, so that the reaction occurring in the microfluidic channel can be monitored more accurately in real time. e collimated p-polarized terahertz beam is injected into the prism and refracted at its base at the internal incident angle θ int . When the momentum matching condition is satisfied, the evanescent wave formed at the base of the prism interacts with the sample in the gap and groove, and the SSPPs wave with exponential decay will be excited on the structure surface of the metal groove. At the frequency for the excitation of SSPPs, the reflectivity will show a sharp dip, known as the SPR phenomenon. erefore, the sample filled in the microfluidic channel can be decided by measuring the frequency shifts or angle shifts of reflection spectrum due to the variation of sample RI. In order to couple the evanescent wave into SSPPs, θ int must be greater than the critical angle of total reflection. e wave vector of incident light in the prism is k 2 P � k 2 ⊥ + k 2 ‖ , where k ‖ and k ⊥ represent the parallel and vertical components of the wave vector, respectively. Note that the conventional excitation of SPR requires that the parallel component of two wave vectors of the evanescent wave, k ‖ , and that of the SSPPs, k SSPPs , satisfy the momentum matching condition [14,15,23,25,26]: However, formula (1) is only suitable for the case where regions Ι and ΙΙ in Figure 1 are air (nd � 1); and when the refractive index nd of the filling medium in regions Ι and ΙΙ is n d >1, formula (1) cannot reflect the effect of the medium on the momentum matching condition. e consequence of using formula (1) will lead to the deviation of the SPR resonance point. erefore, in order to fully consider the effect of the filling medium in the microfluidic channel on the momentum matching condition, formula (1) is revised as formula (2) in this paper: where n p � �� ε p , n d � �� ε d √ , and k 0 � ω/c is the wave vector in the vacuum. In addition, under the limit of a, P ≪ λ, according to the effective medium theory and considering the effect of metal loss [11], the dispersion relationship of SSPPs on the surface of the grooved metal can be expressed as follows [27]: where ε d is the permittivity of the filled medium, δ is the skin depth of the metal, and ε m is the complex permittivity of the metal. According to equation (3), the dispersion relation of SSPPs is determined by the subwavelength structure. In other words, the structure of SSPPs has excellent design flexibility, because the effective plasma frequency can be designed within the desired frequency region, which is very desirable for biosensing.

Simulation Results and Discussion
In order to verify the correctness of formula (2), we use two different numerical simulation pieces of software to compare the results. at is, the dispersion curve of SSPPs wave propagating on the grooved metal at different frequencies is solved by using the eigenmode solver of CST microwave studio software and compared with the reflectivity obtained by the multiphysics field simulation software COMSOL.
Firstly, the dispersion curves of SSPPs supported by 1D periodic rectangular grooves filled with different refractive indexes are analyzed, as shown in Figure 2. e size of periodic unit cell structure is as follows: a � 20 um, h � 30 um, and P � 60 um. e boundary condition is set as follows: the x direction is the periodic boundary, the y direction is set as magnetic wall, and the z direction is set as electrical wall. e fluid refractive indexes corresponding to each dispersion curve of SSPPs in Figure 2 are nd1 � 1, nd2 � 1.025, nd3 � 1.05, nd4 � 1.075, and nd5 � 1.1. e refractive index of the prism is n p � 1.5. In Figure 2(a), the air line (Line_Air) is the dispersion curve of incident light in the air, while the prism lines (Line_np/nd) are the dispersion curves of parallel wave vectors in the prism when the refractive index of the filling sample is nd. In the low-frequency band, the SSPPs dispersion curve is close to the air line but deviates from the air line when approaching the edge of the first Brillouin region, and the degree of this deviation will increase with the increase of nd. e resonant frequency of SPR sensing is determined by the dispersion curve of incident light in the prism and the dispersion curve of SSPPs supported by the metal surface.
at is, when k SSPPs � k ∕ ∕ � k 0 (n p /n d )sin(θ int ), the prism lines for different refractive indexes will intersect the dispersion curves of SSPPs, and each intersection point represents the SPR frequency of different refractive indexes nd. As shown in Figure 2(b), considering formula (2), the resonance points corresponding to the refractive indexes nd1, nd2, nd3, nd4, and nd5 are points A, B, C, D, and E, respectively. Moreover, the prism line (Line_np/nd) will become closer to the air line (Line_Air) as nd gradually increases. Obviously, compared with the resonance points A, B1, C1, D1, and E1 obtained by formula (1), the frequency shift Δf res caused by the variation ∆nd of refractive index in formula (2) will be larger. From the perspective of sensing, formula (2) itself implies that its real sensitivity is higher. However, the influence of refractive index nd on momentum matching condition is not considered in the existing literature [13-16, 19, 20, 23, 25, 26]. In view of the dispersion curve, its sensitivity is lower than that in the real case.
Next, in order to verify the above numerical discussion and theoretical analysis and perform comparison with the dispersion curve calculated by CST, the reflectivity of SPR sensor when filledwith different and was calculated by using multiphysical field simulation software COMSOL based on finite element method (FEM). e schematic diagram of the simulation unit cell of the periodic structure is shown in Figure 3(a), where a = 20 um, h = 30 um, P = 60 um, and g = 140 um. e boundary condition in the propagation direction of SSPPs is set to the Floquet period boundary, and the k vector of the Floquet period is set to . TM polarization wave is incident to the port on the top of the prism at an internal incident angle θ int � 60°, and the propagation constant of the incident port is set to β � k ⊥ � k 0 (n p /n d )cos(θ int ). In the simulation, the permittivity of the metal is ε m � − 1.2 × 10 5 − j7.2 × 10 5 . In order to ensure the convergence and accuracy of the simulation results, the finite element mesh size of each domain is extremely refined. Figure 3(b) shows the reflection spectrum of SPR sensor with different nd as a function of incident frequency, where the range of refractive index is consistent with the previous dispersion curve, from 1 to 1.1, and the interval is 0.025. On the one hand, it is clear that there is a minimum value of  Figure 1: e Otto sensing configuration of the SSPPs excitation in the attenuated total reflection regime using a high-permittivity prism and a periodic metal grating. e width, depth, and period of the groove are a, h, and P, respectively. e gap distance between the prism base and the surface of the metal grating is g; k SPP is the wave vector of the SSPPs wave. θ int is the internal incident angle at the bottom of the prism. e SPR reflectivity is close to zero, and the off-SPR reflectivity is close to one. reflectivity in the reflection spectrum of Figure 3(b), which, respectively, corresponds to the resonance frequency of SSPPs excited at different nd. On the other hand, the resonance point of the reflectivity curve for the different nd in Figure 3(b) is very consistent with the resonance points (A, B, C, D, and E) of the dispersion curve in Figure 2(b), and, with the increase of nd, the resonance frequency gradually red-shifts, which also verifies the necessity and correctness of modifying formula (1) to formula (2).
Furthermore, in order to compare the angle dependence and frequency dependence of SPR sensing at the sensitivity level, taking the resonance frequencies in the dispersion curves of Figure 2 as the incident frequencies, we solved the relationship between the reflectance curves and the incident  angles, as shown in Figure 3(c). It is clear that the reflectivity at each resonance angle θ res also has a distinct dip, but the resonance angle is proportional to the refractive index. e detailed values of the resonance frequency f res and resonance angle θ res in Figure 3(b) and 3(c) are shown in Table 1.

From Figures 3(b) and 3(c)
, it can be seen that as the refractive index of the filled medium deviates from the initial value (nd � 1) corresponding to the optimal coupling gap (g �140 um), the absorption peak of the reflection spectrum begins to broaden and become shallower. is variation of absorption peak can be attributed to the increase of filling refractive index, which not only affects the dispersion curve of SSPPs but also affects the penetration depth of evanescent wave from the prism base into medium. e comparison in optics index between the prism and its contacting sample will be less due to the increase in the refractive index, which will cause a more extension of evanescent wave into the sample gap, and this expansion will lead to an extra reflection on the surface of the metal grating.
at is, the approximately optimal sample gap adopted is only applicable to the structure with the specific initial nd. Hence, when the coupling gap departs from its optimal value g opt , the optimal coupling state will be broken. e more the coupling gap deviates from g opt , the more thoroughly the optimal coupling state broken is. e broken coupling leads to a poorer coupling efficiency and further appears on the reflection spectrum as a broader, shallower reflectivity dip. erefore, in order to study the influence of the gap on SSPPs, we first analyzed the reflectivity of different coupling gaps at nd � 1.025, as shown in Figure 4. Figure 4(a) represents reflectivity as a function of frequency variation, while Figure 4(b) represents reflectivity as a function of incident angle variation. It was apparent that the coupling gap only affects the coupling efficiency of SSPPs, which is reflected in the amplitude of the reflectivity, but has no effect on f res and θ res . at is, for a certain nd, there is an optimal coupling gap g opt , which corresponds to the reflectance curve with the smallest and narrowest absorption peak. For example, the optimal coupling gap when nd � 1.025 is g opt � 162 um. e electric field of SSPPs with different nd is different, which is the fundamental reason why the coupling gap affects the coupling efficiency of SSPPs. e insets in Figure 4(a) and 4(b) from left to right are the surface electric field distribution for different coupling gaps g � 130 um, g � 140 um, and g � 162 um, when the refractive index is 1.025. According to the inset, the electric field for the optimal coupling gap is the strongest, and the larger the coupling gap deviates from the optimal coupling gap, the weaker the electric field is. For example, the maximum values of the electric field for g � 162 um, g � 140 um, and g � 130 um in the insets of According to the above analysis, the sample gap has an optimal value for the highest coupling efficiency with a given refractive index. Here, the optimal sample coupling gaps were obtained by comparing the peak values of the reflectivity at the resonance points with different gap distances adopted (not shown in paper). rough several simulations, the fitting curve of the optimal sample coupling gap for different nd was obtained, as shown in Figure 5. Obviously, with the increase of nd, g increases exponentially. is is because the optimal sample coupling gap is governed by the composite effect of the confinement property of SSPPs and the extension of evanescent wave from the prism base into the sample. Firstly, the complete expression of the evanescent wave is where F(x) � e j(ωt− β 1 x) and α 1 � k 0 n p ������������������ � (sin(θ int )) 2 − (n d /n p ) 2 . Secondly, the transverse electric fields (E 2 ) of spoof SPPs mode decay exponentially away from the corrugated metal surface and can be expressed as [28,29] where α 2 � k 0 ����������������� (n p (sin(θ int ))) 2 − n 2 d . e extension of the evanescent wave is defined as the attenuation length of electromagnetic field to 1/e, which is expressed as follows: Similarly, the confinement of SSPP waves, defined as attenuation length of electromagnetic field to 1/e, is expressed as follows: where k 0 � ω/c and k SSPPs is the wave vector of SSPPs obtained by equation (3). According to the resonance frequencies for the different refractive indexes in Table 1, the curve of L1 + L2 can be obtained by calculation, as shown in Figure 5, and compared with the optimal coupling gap. On the one hand, the electromagnetic field of SSPPs decays exponentially in the direction perpendicular to the metal surface, and the attenuation will slow down with the increase of nd [30]; on the other hand, for an incident wave of a certain frequency, the evanescent wave at the prism base is a bound state, and the depth of its electromagnetic field extending into the gap is limited, but this limitation will be weakened with the increase of nd. ese two aspects will yield a small amount of the "extended" electromagnetic field approaching the edges of the metal groove, which lead to an extra reflection, resulting in a poorer coupling efficiency than the optimal coupling regime. In addition, since the curve trends of L1 + L2 and g opt are consistent and L1 + L2 can be calculated by formulae (8) and (9), in order to show and accurately obtain the existence of an optimum prismmetal gap, we also have fitted the relationship between L1 + L2 and g opt , which can be expressed as International Journal of Antennas and Propagation  where G(L1, L2) is the optimal coupling gap, which is a function of L1 and L2, and alpha and n are the fitting coefficients of the structure. α � 0.042 and n � 1.788. According to the data in Figure 5, under the condition that each nd is matched with its optimal coupling gap, we simulated the reflection spectra of different refractive indexes as a function of resonance frequency and resonance angle, and the results are illustrated in Figures 6(a) and 6(b), respectively. Obviously, since gap only affects the amplitude of the reflectivity and has no effect on the resonance point (f res and θ res ), when the coupling gaps of SSPPs with different nd are optimal, the resonance frequencies of SPR in   Figure 6 are all less than 0.05%, almost close to 0, and the absorption peak is sharper, and the width of the absorption peak is narrower, representing a higher Q-factor and higher sensitivity.
Here, the mutual verification of resonant frequencies obtained by CST and COMSOL and the comparison of frequency dependence and angle dependence of SPR sensing are carried out, respectively, as shown in Table 1. It is clear that the SPR resonance frequencies obtained by the two numerical methods are basically the same, which verifies the correctness of the correction of formula (2) in our previous theoretical analysis. Furthermore, according to the previous theoretical analysis, compared to formula (1), the revised formula (2) itself implies the fact that the sensitivity is higher and up to 3.57 THz/RIU (refractive index unit), which is much more sensitive than in previous work. For example, the sensitivities in literature [14,15,23]  In order to more intuitively reflect the comparison of the two numerical results in Table 1, we draw the fitting curves of the relationship between the resonance frequency and resonance angle and the refractive index of the filling fluid. As shown in Figure 7, the two fitting curves of the resonance frequency are very close, and, according to the fitting curve, we can easily obtain the resonance frequencies for different nd, and, combining formulae (6) and (7), we can calculate the curve of L1 + L2, so as to obtain the approximate optimal coupling gaps based on the proportional relationship of the vertical axis, which provides a new method for determining the optimal sample coupling gap. Moreover, the absolute value 0.74 rad/RIU of the slope of the fitting curve of resonance angle is much less than the slope value 3.57 THz/RIU of resonance frequency, which indicates that the sensing performance of SPR based on frequency variation is better than that of angle variation.
Unfortunately, millions of people die from cancer every year. erefore, the detection of highly sensitive cancerous cells has important biological significance [8,24]. Here, as an application of our proposed biosensing platform based on the improved formula, this article ends with the detection of brain lesions.
at is, the resonance frequency shifts of normal cells and cancer cells infiltrating the gap and groove were calculated by finite element method. Here the refractive index (RI) of the brain cancer cells is shown in Table 2 [31,32]. For a solution with 100% distilled water, the refractive index (RI) is 1.333, which acts as an external    International Journal of Antennas and Propagation reference, and here the CSF with RI of 1.3333 is assumed as an internal reference. According to formula (3), the plasma frequency can be designed by redesigning the geometric parameters of the corrugated metal, so as to effectively excite SSPPs on the surface of periodic structure when the cancer cells and normal cells are infiltrated into the gaps and grooves. It can be seen from Figure 2(a) that the dispersion curve of SSPPs will deviate further apart from the light line (Line_Air) as the refractive index increases, but the parallel component k ∕ ∕ of the wave vector in the prism is closer to the light line. Note that when nd increases to a certain value (of course, since SSPPs are a slow wave, i.e., k SSPPs /k 0 > n d , the theoretically detectable sample maximum refractive index is n d < ��������� n p sin(θ int )), there will be no intersections between the parallel component of the wave vector in the prism and the dispersion curve of SSPPs.
at is, SSPPs cannot be effectively excited. Here, when the filling sample is the brain cancer cells, in order to effectively excite SSPPs, we redesigned the geometric dimensions of the metal grooves as follows: a � 20 um, h � 50 um, and P � 60 um. e permittivity of the prism is ε p � n 2 p � 11.63. As mentioned in the previous theory and simulation analysis, the sensing sensitivity of SPR in the frequency domain is better than that in the angle domain, so next we will mainly conduct the reflection spectra of brain lesions in the frequency region. Here, we carry out simulation analysis for the two cases, θ int � 45°and θ int � 60°, and when the metal groove and the gap space under the prism were infiltrated with brain cancer cells, the reflectivity at different incident frequencies is calculated by COMSOL, respectively, as illustrated in Figures 8(a) and 8(b); the reflectance spectrum of the proposed sensing device manifests that the SPR frequency of cancer cells will red-shift as the refractive index increases. Note that the reflection spectrums in Figure 8 are calculated under the optimal coupling gaps corresponding to different cells, and the determination of g opt is similar to the process discussed above. In addition, in order to more clearly show the change of resonance frequency in Figure 8, we draw the resonance positions of brain lesions at different incident angles in Figure 9. It can be seen from Figure 9 that the resonant frequency corresponding to 45 degrees is lower than that corresponding to 60 degrees. is is because the larger the incident angle, the larger the wave vector k of the SSPPs, and the smaller the frequency of the SPR. e sensitivity is a key factor in evaluating the efficiency of a sensor. For the sensitivity of the refractive index sensor, we define it as the frequency shift value divided by the refractive index change of the biological sample: where Δf res � f res (cancerous cell) − f res (normal cell) and Δn d � n d (cancerous cell) − n d (normal cell). Besides, the figure of merits (FOM) is another significant indicator of sensor performance, which can be defined as where Δf 1/2 is the spectral full width at half maximum of the reflection dip. Obviously, the sensitivity when the incident angle is 45 degrees is better than that when the incident angle is 60 degrees. e electric field confinement ability of the SSPPs on the surface of the periodic groove structure at two incident angles is the main reason for this difference in sensing performance. Here, the CSF are considered as a reference and we calculated the electric field intensity normE � �������������������������� � abs(Ex) 2 + abs(Ey) 2 + abs(Ez) 2 , as shown in Figure 10. It is clear that the electric field for 45 degrees is stronger than that for 60 degrees; that is, the interaction between the biological sample and SSPPs at 45 degrees is stronger than that at 60 degrees, which leads the SSPPs to be more sensitive to changes in the dielectric environment near the groove surface.
In addition, in order to more clearly describe the detection effect of different degrees of brain lesions, we plotted the relationship between S and FOM with brain cancer cells in Figures 11(a) and 11(b). It can be seen intuitively from Figure 11 that this structure has excellent performance in detecting different degrees of brain lesions, and, with the deterioration of cancerization, the detected S value and FOM value will be larger.
Finally, Table 4 compares the sensitivity and FOM value of proposed structure in this research with the previous designed sensors operated in the THz region. e results show that the performance of the biosensor based on the modified momentum matching formula proposed in this paper is far superior to previous work. Owing to its high sensing performance above, it could be believed that the proposed SSPPs sensor based improved momentum matching condition would provide effective guides for achieving high-sensitivity RI biosensor in THz region.      Figure 11: e sensitivity (a) and FOM (b) for different brain cancerous cells at different incident angles.

Conclusions
In summary, we have numerically investigated the performance of THz SSPPs-based sensing of refractive index changes under the premise that the filling sample refractive index will affect momentum matching condition. Firstly, we modify the dispersion formula of the parallel component of the wave vector in the prism and verify the correctness of the formula from the two aspects of dispersion curve and reflectivity. e results show that reflection amplitude depends on the gap distance. For a given filled sample with a certain nd, there is an optimal coupling gap for minimal reflection amplitude. We fitted the relationship curve between the optimal coupling gap and the filling nd. en, the reflectivity is calculated on the basis of adjusting the optimal coupling gap distance, so that the coupling of SSPPs with different nd is in the optimal regime; that is, the absorption peak of the reflectivity for each nd is very sharp, and the minimum value of reflectivity is below 0.05%, which realizes the high-sensitivity detection of specific analytes. e SPRbased sensing performance of the frequency region is superior to that of the angle region, and the sensing results of filling fluids with different nd show that the sensitivity of the improved formula (2) is as high as 3.57 THz/RIU and the corresponding FOM is 3966.
In addition, as a sensing application platform, we conducted the brain lesions sensing with different types of cancerous cells. e numerical results show that the maximum sensitivity is about 2 THz/RIU, and the corresponding FOM is 3337, which is significantly better than the previously designed sensors that work in the terahertz region.
Moreover, since the response of the SSPPs can be engineered by designing its geometrical sizes, the sensing for biological samples can be performed in the desired frequency range, which is crucial for biological and security sensing based on the terahertz microchannel.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare no conflicts of Interest.