## Polarization-controlled non-Hermitian metasurfaces for ultrasensitive terahertz sensing

Xintong Shi  , Hai Lin 1 1, * , Tingting Liu 2,3 , Yun Shen 4, * , Rongxin Tang  , Le Li  , 4 4

Junyi Zhang 1, 4 , Yanjie Wu 1, 5 , Shouxin Duan  , Chenhui Zhao  , Shuyuan Xiao 4 4 2, 3, *

- 1 *
- College of Physics Science and Technology, Central China Normal University, Wuhan 430079, China

2 School of Information Engineering, Nanchang University, Nanchang 330031, China 3 Institute for Advanced Study, Nanchang University, Nanchang 330031, China 4 Institute of Space Science and Technology, Nanchang University, Nanchang 330031, China 5 School of Electronic Science and Engineering, University of Electronic Science and Technology of

- China, Chengdu, 611731 China

*Corresponding author(s). E-mail(s): linhai@mail.ccnu.edu.cn; syxiao@ncu.edu.cn shenyun@ncu.edu.cn;

## Abstract

Non-Hermitian systems offer significant advantages in sensor design, especially at the exceptional points.

However, the extreme sensitivity near these points poses great challenges due to fabrication errors and system  noises,  which  degrade  sensing  performance. To  address  this,  we  introduce  a  novel  approach leveraging the polarization degrees of freedom in non-Hermitian systems. In this work, we establish a direct relation between the incident polarization and the transmission phase of a coupled metasurface system  and  achieve  the  polarization-controlled  phase  singularity  even  post-fabrication.  The  incident polarization angle can be utilized as a sensing index, which enables indirect and accurate measurement. The  theoretical  approach  is  experimentally  validated  using  a  general  design  of  THz  non-Hermitian metasurface sensors. Our method enhances robustness and sensitivity, opening new avenues for practical applications in ultra-sensitive sensing.

## Introduction

Any natural, non-isolated physical systems are intrinsically non-Hermitian, because they continually exchange  energy  with  the  surrounding  environment 1,2 .  The  non-Hermitian  phenomena  have  been observed across various fields such as optics, acoustics, microwaves, and electronics 3-9 . The intriguing characteristics of these phenomena have inspired numerous exciting applications 10-13 , for example, nonHermitian state-based topological lasers 14 ,  unclonable physical keys 15 ,  and  more flexible wavefront manipulation  and  topological  engineering  capabilities 16,17 .  Additionally,  leveraging  the  eigenvalue splitting near exceptional points (EPs ） 18-20 , which are among the most renowned exotic states in nonHermitian systems, can facilitate the development of high-performance sensors with diverse photonic structures like optical resonators, microcavity, and metasurfaces 21-29 . Among them, metasurfaces are highly favored due to their unprecedently flexible and versatile design strategies 30-37 , and in particular provide feasible solutions in previously inaccessible THz regime. Their significantly enhanced local field 38-40 , combined with the unique characteristics of THz waves, present tremendous advantages in detecting subtle biological signals 41-43 .

Despite  EPs'  superior  performance,  a  critical  issue  has  been  consistently  overlooked  in  the  existing design strategies for previous non-Hermitian sensors. The ultra-sensitivity near the EPs, offering both benefits and drawbacks. Eigenvalue variations under a perturbation are not a good measure of the overall sensor  performance  near  an  EP 44 .  The  factors  causing  eigenvalue  splitting  include  not  only  the perturbation from the analytes but also fabrication errors, which would lead to a significant degradation in sensing performance. Considering the geometric parameters are fixed once the structures are fabricated, some studies have introduced gain media to compensate for uncertainty 45-47 . However, such approach is not suitable for periodic structures like metasurfaces, let alone the lack of such gain media in the terahertz regime.  This  issue  appears  to  be  a  deadlock,  severely  obstructing  the  practical  application  of  nonHermitian metasurface sensors.

In this work, we resort to the polarization degree of freedom to achieve controllable singularity in nonHermitian systems for sensing. We theoretically construct an effective Hamiltonian 48 by incorporating the coupling efficiency loss which comes from the mismatch between the linearly polarized light and the excited eigenmodes of the structure. By solving its eigenvalues, we derive the explicit relation between the polarization of incident light and the transmission phase of coupled system, and observe the phase singularity in the polarization space. This phase singularity exhibits high tolerance to fabrication errors. Based on this, the polarization of incident light, which can be easily tuned for phase singularity even after  fabrication,  is  utilized  as  a  sensing  index.  Finally,  the  theoretical  method  is  experimentally demonstrated  using  a  general  design  of  THz  metasurface  sensors  and  measured  with  the  THz  time domain  spectrometer.  The  proposed  method  explores  the  new  degree  of  freedom  in  non-Hermitian physics, and open an avenue for robust ultra-sensitive sensing.

## Results

## Theoretical derivations

In previous studies, non-Hermitian metasurfaces have typically relied on circular polarization or strictly linear  polarization.  Experimental  measurements  of  the  linear  polarization  or  synthesized  circular polarization usually require strict  horizontal  and vertical  polarization  to evaluate  co-polarization  and cross-polarization. Here we introduce a sensing method that incorporates polarization degrees of freedom.

Figure 1 shows the schematic illustration of the proposed non-Hermitian metasurface sensor and its working  mechanism.  The  unit  cell  of  the  metasurface  sensor  is  constructed  by  two  coupled  and orthogonal  resonators.  The resonators possess similar frequencies but different intrinsic rates, respectively. A linearly polarized light is incident upon the metasurface with a controllable polarization angle 𝜃𝜃 , and a detector collects signals corresponding to the horizonal and vertical polarizations aligned with two resonators, respectively. By rotating the first polarizer to control the incident polarization angle, the  excitation  strength  of  two  resonators  can  be  precisely  tuned.  Due  to  the  complementary  relation between the incident polarization angle and the orthogonal resonators, the transmissions in horizontal and vertical directions exhibit a distinct behavior. One exhibits phase anomaly and the other does not. For the anormal one, we also observe that the coupled system is expected to exhibit a polarizationcontrolled phase singularity of transmission (which is extremely sensitive to any subtle change in the surrounding environment) at a certain  angle.  Consequently,  the  sensing  strategy  here,  constructs  the method to characterize the analyte using mechanical control for incident polarization. The subfigures illustrate a phenomenon where the polarity of group delay reverses, providing a clearer representation of phase singularity. On either side of the critical polarization angle, the calculation of the winding number for  the  transmission  spectrum  reveals  that  the  transmission  signals  corresponding  to  the  incident polarizations exhibit different topological invariants. This discrepancy serves as the origin of the phase singularity. Next, we proceed to construct the equivalent Hamiltonian.

Fig 1. The conceptual illustration of the proposed non-Hermitian metasurface sensor and working mechanism. The unit cell consists of two orthogonal resonators. A pair of polarizers, functioning as a polarizer and an analyzer respectively, controls the polarization of the incident light and analyzes the polarization of transmission into two directions along the resonators (horizontal and vertical directions). The phase changes from lagging to leading, the corresponding group delay exhibits a transition, and the topological invariant of transmission also changes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00579v3-figure-d6509b70-a1ba-4070-9973-8c8c0cc83825.jpg)

For  an  isolated  resonator,  when  the  direction  of  incident  light  aligns  with  that  of  its  eigenmode,  it experiences the strongest excitation and harvests the highest energy, and when the incident polarization directions are orthogonal, it acquires the lowest energy. In a more general case, if the incident polarization directions form an angle, it would lead to losses between the resonator and the incident field. For a coupled system composed of two resonators, the situation becomes more complex as the resonators can exchange energy with each other. If the two resonators are orthogonal, the losses that come from their coupling with the incident light would change in opposite ways. We bring these polarization-induced losses  into  consideration  (for  more  information  please  see  supplementary  Note  SN1).  The  effective Hamiltonian of the coupled system can be written as

$$H = \begin{bmatrix} \omega _ { 1 } - i \gamma _ { 1 } - i \Gamma _ { 1 } & \kappa e ^ { i \beta } \\ \kappa e ^ { i \beta } & \omega _ { 2 } - i \gamma _ { 2 } - i \Gamma _ { 2 } \end{bmatrix},$$

where 𝜔𝜔 1 2 , are the frequency of the resonators, 𝛾𝛾 1 2 , are the intrinsic loss of the resonators, and Γ1 2 , are the  polarization-induced  losses.  The  coupling  strength  in  coupled  system  is  denoted  as 𝜅𝜅𝑒𝑒 𝑖𝑖𝑖𝑖 . The complex frequencies are denoted Ω1 and Ω2 , and the corresponding eigenvalues can be denoted as

$$\lambda _ { \mp } & = \bar { \Omega } \pm \sqrt { \Delta \Omega ^ { 2 } / 4 + \kappa ^ { 2 } e ^ { 2 i \beta } },$$

Ω1+Ω2

where Ω � = 2 , and ΔΩ = Ω1 - Ω2 . It is obvious that the two eigenvalues coalescing at EP with 𝜆𝜆 ∓ = Ω � when  the  equation  satisfying  with  the  discriminant  of  a  quadratic  equation ΔΩ 2 /4 + 𝜅𝜅 2 𝑒𝑒 2𝑖𝑖𝑖𝑖 = 0 .

While the physics of EPs may appear straightforward, it is often challenging to reach the point in practice because  EPs  require  quite  stringent  structural  parameters.  Therefore,  in  the  following,  we  introduce tuneable losses by leveraging the polarization degree of freedom, which allows us to dynamically control the coupled system response.

Therefore, the eigenvalues are re-written as

$$\lambda _ { \mp } \, = \, \overline { \Omega } \pm \sqrt { k ^ { 2 } c o s 2 \beta + \Delta \omega ^ { 2 } / 4 - ( \Delta \gamma + \Delta \Gamma ) ^ { 2 } / 4 - i \kappa ^ { 2 } s i n 2 \beta - i \Delta \omega ( \Delta \gamma + \Delta \Gamma ) } \,, \quad \ \ ( 3 )$$

where Δ𝑖𝑖 is the difference between polarization-induced losses of the resonators. Considering the trend of  this  difference,  the  system  responses  should  have  an  extremum  (for  more  information  please  see supplementary note SN1 and Figure S3). It is obvious to find the zero occurs at Ω � -�ΔΩ 2 /4 + 𝜅𝜅 2 𝑒𝑒 2𝑖𝑖𝑖𝑖 and Ω � + �ΔΩ 2 /4 + 𝜅𝜅 2 𝑒𝑒 2𝑖𝑖𝑖𝑖 .

Assuming that the coupled system satisfying the conditions of zero, Ω � = �ΔΩ 2 /4 + 𝜅𝜅 2 𝑒𝑒 2𝑖𝑖𝑖𝑖 . Under such condition, the Magnitude 𝑀𝑀𝜆𝜆 = 0 , and the phase expression is,

$$\varphi _ { \lambda _ { - } } = \arg ( \lambda _ { - } ) = \arctan \frac { - \omega _ { 1 } ( \Gamma _ { 2 } + \gamma _ { 1 } ) - \omega _ { 2 } ( \Gamma _ { 1 } + \gamma _ { 2 } ) - \kappa ^ { 2 } s i n 2 \beta } { \omega _ { 1 } \omega _ { 2 } - ( \gamma _ { 1 } + \Gamma _ { 1 } ) ( \gamma _ { 2 } + \Gamma _ { 2 } ) - \kappa ^ { 2 } c o s 2 \beta } }.$$

The  phase  expression  has  a  singularity  point  when  the  denominator  resembles  a  critical  coupling condition ( 1 -𝛼𝛼 2 -𝜅𝜅 2 = 0 ). If the two resonators share the same eigenfrequency, the total losses equal the energy exchanged between the two resonators. Both of the real and imaginary parts 𝜆𝜆 -are  zero, which means the non-transmission appears along the direction of this resonator.

Figure 2a illustrates the parameter space of the coupled system under tight coupling ( 𝛽𝛽 = 0 ) as functions of frequency difference Δ𝜔𝜔 and incident polarization angle 𝜃𝜃 , as indicated by the equation 2. The upper and lower panels show the real and imaginary parts of eigenvalues, respectively. The real eigenvalues not  only  separate  as  the  frequency  difference  increases,  but  are  also  influenced  by  the  incident polarization, with the effect on the real eigenvalues diminishing as the frequency offset between the two resonators  increases.  The  ring  appears  on  the 𝜃𝜃 axis  when  the  two  resonators  share  the  same eigenfrequency ( Δ𝜔𝜔 = 0 ).

The imaginary part of eigenvalue is mainly controlled by the incident polarization. When the polarization angle is distant  from  the quadrant center (45 ° ),  the  imaginary  eigenvalues  separate. Otherwise,  the imaginary  eigenvalues  degenerate.  When  the  frequency  offset  increases,  the  imaginary  eigenvalues exhibit only an intersection point. Otherwise, they degenerate into a line along the 𝜃𝜃 axis.

Fig. 2 The theoretical calculations (a). The parameter space of coupled system with phase difference Δ𝜔𝜔 and incident polarization angle 𝜃𝜃 . The upper panel and down panel are real part and imaginary part of eigenvalues of the effective Hamiltonian. The two parts of eigenvalues form a closed loop and a line close to critical polarization angle, respectively. (b-c). The Magnitudes of two eigenvalues and the phases of two eigenvalues. The results are calculated by the equations (2) and (4). The parameters used for the calculations are 𝛾𝛾 1 = 0.2 , 𝛾𝛾 2 = 0.13 , and 𝜅𝜅𝑒𝑒 𝑖𝑖𝑖𝑖 = 0.03 𝑒𝑒 𝑖𝑖0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00579v3-figure-219c3789-211a-4e29-9bb7-6084bb3801ea.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00579v3-figure-3952b7ff-600a-4cba-abd9-45ccfbe385fa.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00579v3-figure-044ce8a8-161a-4fac-be6b-ce781ffce420.jpg)

In this work, we concentrate on the system's phase singularity. Figures 2b and 2c show the calculated magnitude and phase of eigenvalues based on equation 4. The magnitude of 𝜆𝜆 -goes to the minimum, approaching to zero (see the legend of 2b), and the phase exhibits a singularity, when Δ𝜔𝜔 equals zero. The phase changes from lagging to leading when the incident polarization crosses the critical angle. The magnitude of 𝜆𝜆 + is much larger than 𝜆𝜆 -, and the phase is trivial.

It is worth noting that the magnitude and phase in Figures 2b and 2c correspond to the eigenvalues of the effective Hamiltonian, representing the effects of the Hamiltonian operator on the incident field. These values do not include the propagating phase and therefore represent the intrinsic (non-spectral) phase of the effective Hamiltonian rather than the directly measured phase in the spectra. In addition, even if there is a certain offset between the eigenfrequencies of two resonators, the phase singularity can still exist, albeit appearing at different incident polarization angles. It is important to note that the phase singularity is influenced by the incident polarization, thus establishing a polarization-controlled singularity system.

Consequently,  this  relation  allows  us  to  establish  a  sensing  function  with  respect  to  the  incident polarization angle. The optical response of the coupled system would change after a small amount of analyte is added onto the metasurface. In equation 4, the resonance frequencies 𝜔𝜔 1 2 , (main) and intrinsic losses 𝛾𝛾 1 2 , (very  small  in  the  THz  regime)  are  influenced  by  the  perturbations,  but  the  polarizationinduced losses Γ1 2 , depend on the incident field. Then the equation 𝜔𝜔 1′𝜔𝜔 2′ - 𝜅𝜅 2 𝑐𝑐𝑐𝑐𝑐𝑐 2 𝛽𝛽 = ( 𝛾𝛾 1 + Γ1 )( 𝛾𝛾 2 + Γ2 ) determines the new incident angle. Hence, it bridges the two aspects of polarization and analyte, allowing us to characterize the perturbation in the polarization space, which links the electromagnetic parameters of the metasurface with external mechanical parameter.

## Design and simulation results

We  employ  the  finite  integration  technique  (FIT)  to  validate  the  proposed  theoretical  method.  The metasurface sensor consists of simple and classical structures, namely split ring resonators (SRRs) and cut  wires  (CWs).  The  selection  of  these  simple  structures  ensures  manageable  responses  while maintaining sufficient coupling effects. The parameters of the structures, including the SRRs and CWs, are estimated using the ruler equation 49 . By combining these two structures into a single unit, a tightly coupled  system  is  constructed.  Figure  3a  provides  a  schematic  illustration  of  the  metasurface configuration. We  conduct  simulations  to  obtain  the  phase  response  of  the  metasurface  without  any analytes under different incident polarizations. The designed parameters of the metasurface sensor are presented in Table 1(for more information please see supplement information SN2 and SN3), and the corresponding  results  are  plotted  in  Figure  3b.  The  upper  panel  shows  the  transmission  phase  of horizontal  direction,  and  the  lower  panel  shows  the  transmission  phase  of  vertical  direction.  The responses have an approximate frequency offset of 0.11 THz. The responses are consistent with the theoretical calculations: one exhibits a phase singularity, while the other is trivial.

Table 1 Parameters of the non-Hermitian THz metasurface sensor

| Parameters                                                          | Symbols   |   Values ( μ m ) |
|---------------------------------------------------------------------|-----------|------------------|
| Length of SRR in the x- direction                                   | L 1x      |             60   |
| Length of SRR in the y- direction                                   | L 1y      |             50   |
| Length of CWs                                                       | L 2       |             70   |
| Width of SRR                                                        | W 1       |             10   |
| Width of CWs                                                        | W 2       |             10   |
| Gap of SRR                                                          | G         |             10   |
| Distance between SRR and CWs in the x- direction                    | d x       |             42.5 |
| Displacement of SRR from the center of unit cell in the y direction | d y       |             20   |

We proceed with the phase singularity. It can be observed that under different linearly polarized excitation, the transmitted phase lags on one side of the singularity point, and upon crossing the singularity, on the other side, the lagging phase transitions to leading. This behaviour is consistent with the results shown in Figure 2. What distinguishes this from the calculations in Figure 2 is that here we have the real phase of the resonators, including the propagation phase. The phase singularity appears at 42 ° and 0.827 THz. The phase variation becomes pronounced near the singularity point. We also investigate the topological origin used winding number for this phase singularity, as included in supplementary information SN5.

Retrieving the phase singularity becomes challenging when considering the influence of propagation phase,  and  it  is  also  difficult  to  characterize  it  as  a  single  point.  From  an  engineering  perspective, analyzing the differential form of the phase is much more convenient. When the extremum of the group delay ( 𝑑𝑑𝜑𝜑 / 𝑑𝑑𝜔𝜔 ) undergoes a transition, it indicates that the phase has crossed the singularity. We also illustrate the corresponding group delays of the two directions in the inset of Figure 3b. In other words, the extremum of the group delay corresponds well to the phase singularity. Moreover, since the derivative of the propagation phase with respect to frequency is constant, its effect is directly eliminated in this process.

Fig. 3 Simulation results. (a). The schematic of the non-Hermitian THz metasurface sensor. (b). The transmission phase (group delay) of the metasurface without loading analyte, (c). The relation between the incident polarization angle ( x -axis) and refractive index ( y -axis), with group delay indicated by color. (d). The three refractive indexes are extracted from Figure 3c to illustrate the results. (e) and (f) are the group delay with varying the thickness of analytes, with the refractive index are 1.2 and 1.4, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00579v3-figure-0d742e2a-25e5-47a3-8393-b594a68955f2.jpg)

In the FIT simulations, a 20 μ m dielectric layer is introduced onto the metasurface sensor to analog the analyte 50 .  The  range  of  analyte  refractive  index  considered  is  from  1  to  1.2,  while  the  incident polarization angle is varied between 40 and 50 degrees. The simulation results are presented in Figure 3c,  where  each  pixel  corresponds  to  once  simulated  group  delay  (for  more  information  please  see supplement information SN4). The horizontal axis represents the incident polarization angle, while the vertical axis represents the refractive index of the dielectric layer. The color scale indicates the extreme values of the group delay. Notably, a distinct transition is observed, indicating a relation between the two variables characterized by a series of singularities. Beside the transition region, the group delay exhibits different symbols, with one side being positive and the other side being negative. This transition signifies that  within  the  given  range,  each  value  of  the  refractive  index  corresponds  to  a  unique  incident polarization  angle. Thus,  the  result  demonstrates  the  existence  of  a  specific  correlation  between  the refractive index and the incident polarization angle within a defined range.

To illustrate  this  relation  clearly,  three  specific  refractive  index  parameters  (1.05,  1.1,  and  1.15)  are selected from Figure 3c. Figure 3d depicts the corresponding group delay values for these refractive indexes. The behavior of the three lines aligns with the theoretical predictions. Notably, at the incident polarization angles of 43°, 45°, and 46°, the group delays undergo a transition from positive to negative, revealing the presence of phase singularities. In essence, these findings demonstrate that the incident polarization angles can effectively represent the refractive index.

However, it is important to note that the refractive index of analytes in different concentrations is not the only  factor  affecting  terahertz  biosensing.  The  effective  refractive  index  varies  depending  on  the properties of the substrate and the attached material. Additionally, terahertz detection often involves the removal  of  moisture,  further  complicating  the  use  of  refractive  index  alone  for  accurate  terahertz biosensing. To address this, it is necessary to simulate different concentrations of analytes by employing dielectric blocks with varying thickness.

Here, we simulate the metasurfaces coated with a dielectric layer of constant refractive index, with the thickness ranging from 0 to 20 μ m . Two types of dielectric layers with distinct refractive indexes are used to emulate different analytes. The corresponding results are presented in Figure 3e and Figure 3f, revealing that the sensitivity to the thickness parameter is non-linear. While the thickness variation of the dielectric layer is linear, the corresponding incident polarization angle exhibits a non-linear response. Moreover,  as  the  dielectric  layer  becomes  thicker,  the  incident  polarization  angle  undergoes  slower changes.  Additionally,  the  incident  polarization  angle  varies  differently  depending  on  the  type  of dielectric  layer,  demonstrating  the  sensor's  capability  to  specifically  distinguish  between  different analytes. Furthermore, the system's response is more sensitive when a thinner dielectric layer is employed, and  the  curves  corresponding  to  different  refractive  indexes  (representing  different  analytes)  exhibit distinct patterns. Based on this principle, we attribute the sensor with the capability to accurately identify different analytes.

## Experimental validations

Finally, we also conduct experimental demonstrations of the metasurface sensors to validate the proposed method. The analyte used in the experiments is aflatoxin B1 with varying concentrations. During the measurement process, a switchable linearly polarized terahertz wave is illuminated onto the metasurface sensor at a specific polarization angle 𝜃𝜃 .  The  incident  polarization  angle is  precisely  controlled  by a rotating polarizer. The transmitted signal is then separated into the horizontal axis using a polarizer and subsequently detected by a receiver.

Additionally, we implement the key measure to enhance measurement precision. To minimize the need for frequent adjustments of the polarizer angle, multiple metasurfaces are fabricated, each coated and dried with aflatoxin B1 at different concentrations. This process plays a crucial role in improving the accuracy of measurement. The calibration procedure involved setting the first polarizer to have the same polarization as the photoconductor antenna, while adjusting the polarization of the second polarizer to be  orthogonal  to  the  first.  The  second  polarizer  is  carefully  aligned  to  achieve  calibration  when  the transmission amplitude reaches its minimum. The angle of the first polarizer is then set by selecting points at 5-degree intervals around 45 degrees. After each setting, the unload sensor and other sensors with varying aflatoxin B1 concentrations are sequentially placed on the light path, and a set of responses is measured.

Fig. 4. Experiment results. (a). The fabricated sample of metasurface sensor. (b). The magnified imaging of metasurface sensors coated with aflatoxin B1 of different concentrations, (c). The group delay calculated using a set of 5×5 measurements, revealing the relation between incident polarization angle and concentration of analytes. (d-h). The fitted group delay based on measurement data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00579v3-figure-491c80cd-2ced-4d5d-bb30-f9b59ab684b1.jpg)

Figure 4a presents a photograph of the fabricated sensor. Figure 4b illustrates magnified photographs of the metasurface sensors coated with different concentrations of aflatoxin B1 (unloaded, 0.001, 0.01, 0.1, and 1 μ g/mL , respectively). In these images, the gray region represents the silicon substrate, the white region corresponds to the aluminum structure, and the black dots represent the analytes. The processed group delay as shown in Figure 4c (the raw data pleases see supplementary information SN 6), the same method depicts the group delay between incident polarization and concentration of aflatoxin B1. The results have verified the relation, it is as similar as the anticipation of theoretical and simulation results.

The phase transition of the unload sensor is observed between 40 and 45 degrees. For sensors with concentrations of about 0.001 μ g/mL , 0.01μ g/ mL , and 0.1 μ g/mL , the phase transitions appear between 45 and 50 degrees. Sensors with a concentration of 1 μ g/mL of aflatoxin B1 exhibit a transition between 50 and 55 degrees. Based on the proposed theory, we can solve the singularity point of each concentration.

In fact, it is not necessary to precisely measure each of group delay under different incident polarization angles. The location of the phase singularity can be determined through multiple measurements on an experimental platform that has not been specially modified. The mathematical form of the phase related to the incident polarization has been established in equation 4. The constant parameter can be determined by switching the incident polarization and measuring the corresponding responses without loading any analyte. Subsequently, when a perturbation is introduced into the system, the new polarization angle associated with the singularity can be determined. Hence, the location of phase singularity can be fitted based on the measured data, when the analytes are added onto the metasurface. Then, the entire group delay with the incident polarization from 40 to 60 degrees are depicted in Figure 4d-h. It is clearly that each concentration of aflatoxin B1 matches one phase singularity. The result corroborates the theory, the singularity on the curve is clear. This method mainly adds precision of measurement for phase singularity. Employing  this  method  helps  to  avoid  unstable  measurements  at  the  singularity,  thereby  improving measurement  accuracy.  Furthermore,  compared  to  direct  measurements,  responses  obtained  from locations far away from the singularity exhibit greater stability and are less susceptible to interference. This leads to more accurate results and allows for solution derivation based on a set of measurement data, thereby reducing the impact of measurement errors on the final outcomes.

Table 2 shows the concentration of aflatoxin B1, corresponding incident polarization angle. Considering the define of sensitivity of traditional metasurface sensors 50-55 , we are not able to take comparison of sensitivity with traditional methods directly. Compared to our previous metasurface sensor 56 , the sensing indicator changes distinctly when the analyte is in trace levels. More importantly, compared to the change of amplitude and frequency, the incident polarization angle is an independent mechanical quantity which can be amplified by mechanical methods and geometrical optics. Currently, the measurement is also limited by the spectrometer, which is a challenging. However, we believe that with breakthroughs in technologies, the advantages of this method will surpass other sensors.

Table2 The concentration and polarization

| Concentration ( μ g /mL)   |     0 |   0.001 |   0.01 |   0.1 |     1 |
|----------------------------|-------|---------|--------|-------|-------|
| Polarization （ degree ）  | 43.31 |   46.72 |  47.54 | 49.78 | 50.67 |

## Discussion

In  summary,  we  have  proposed  a  novel  THz  sensing  method  based  on  polarization-control  phase singularity  to address  the  current  challenges  faced  by  non-Hermitian  metasurface  sensors.  We demonstrate  the  second-order  non-Hermitian  system  using  a  metasurface  and  controllable  linearly polarized  THz  wave.  By  introducing  the  concept  of  polarization-controlled  losses,  we  establish  a connection between two different types of physical quantities, the polarization of the incident light and the phase singularities of the coupled metasurface system. This method is not only well-suited in THz regime, but can also be extended to other frequencies including microwave (for more results please see supplementary information SN 7 and the attached video file), infrared, and visible. Furthermore, we develop a sensor characterization method based on the phase expression of the system's eigenvalues, using  the  polarization  rotation  angle  as  an  indicator.  We  also  propose  a  method  to  locate  phase singularities,  which  circumvents  challenges  associated  with  directly  measuring  the  system  response, significantly enhancing the sensing performance. Note that without sacrificing generality, here we have opted for  classical  meta-structures  in  designing  the  sensor  for  the  proof  of  concept,  and  the  sensing performance can be further improved through the integration of advanced designs, nonlinear effects, and higher-order coupling. We aspire for this sensing method to augment the application potential of nonHermitian THz metasurface sensors, and in particular, enhance the stability of non-Hermitian sensor measurements.

## References

- 1. Colom, R. et  al. Crossing  of  the  Branch  Cut:  The  Topological  Origin  of  a  Universal  2 πPhase Retardation in Non-Hermitian Metasurfaces. Laser Photonics Rev. 17 , 2200976 (2023).
- 2. Li,  C. et  al. Experimental  Demonstration  of  Controllable  PT  and Anti-  PT  Coupling  in  a  NonHermitian Metamaterial. Phys. Rev. Lett. 132 , 156601 (2024).
- 3. Ergoktas,  M.  S. et al. Topological  engineering  of  terahertz  light  using  electrically  tunable exceptional point singularities. Science 376 , 184-188 (2022).
- 4. He, T. et al. Scattering exceptional point in the visible. Light Sci. Appl. 12 , 229 (2023).
- 5. Park, S. H. et al. Observation  of  an  exceptional  point  in a non-Hermitian  metasurface.

Nanophotonics 9 , 1031-1039 (2020).

- 6. Rao, J. W. et al. Controlling Microwaves in Non-Hermitian Metamaterials. Phys. Rev. Appl. 15 , L021003 (2021).
- 7. Baek, S. et al. Non-Hermitian chiral degeneracy of gated graphene metasurfaces. Light Sci. Appl. 12 , 87 (2023).
- 8. Lawrence, M. et al. Manifestation of P T Symmetry Breaking in Polarization Space with Terahertz Metasurfaces. Phys. Rev. Lett. 113 , 093901 (2014).
- 9. Wang, D. et al. Superconductive PT-symmetry phase transition in metasurfaces. Appl. Phys. Lett. 110 , 021104 (2017).
- 10. Song, W. et al. Breakup and Recovery of Topological Zero Modes in Finite Non-Hermitian Optical Lattices. Phys. Rev. Lett. 123 , 165701 (2019).
- 11. Liang, C., Tang, Y., Xu, A.-N. &amp; Liu, Y.-C. Observation of Exceptional Points in Thermal Atomic Ensembles. Phys. Rev. Lett. 130 , 263601 (2023).
- 12. Chen, W., Kaya Özdemir, Ş., Zhao, G., Wiersig, J. &amp; Yang, L. Exceptional points enhance sensing in an optical microcavity. Nature 548 , 192-196 (2017).
- 13. De Carlo, M., De Leonardis, F., Soref, R. A., Colatorti, L. &amp; Passaro, V. M. N. Non-Hermitian Sensing in Photonics and Electronics: A Review. Sensors 22 , 3977 (2022).
- 14. Schumer, A. et al. Topological modes in a laser cavity through exceptional state transfer. Science 375 , 884-888 (2022).
- 15. Yang, M. et  al. Electromagnetically unclonable functions generated by non-Hermitian absorber-

emitter. Sci. Adv. 9 , eadg7481 (2023).

- 16. Song, Q., Odeh, M., Zúñiga-Pérez, J., Kanté, B. &amp; Genevet, P. Plasmonic topological metasurface by encircling an exceptional point. Science 373 , 1133-1137 (2021).
- 17. Yang, Z. et al. Creating pairs of exceptional points for arbitrary polarization control: asymmetric vectorial wavefront modulation. Nat. Commun. 15 , 232 (2024).
- 18. Özdemir, Ş. K., Rotter, S., Nori, F. &amp; Yang, L. Parity -time symmetry and exceptional points in photonics. Nat. Mater. 18 , 783-798 (2019).
- 19. Ashida, Y., Gong, Z. &amp; Ueda, M. Non-Hermitian physics. Adv. Phys. 69 , 249-435 (2020).
- 20. De Carlo, M. (INVITED) Exceptional points of parity-time- and anti-parity-time-symmetric devices for refractive index and absorption-based sensing. Results Opt. 2 , 100052 (2021).
- 21. Wiersig, J. Sensors operating at exceptional points: General theory. Phys. Rev. A 93 , 033809 (2016).
- 22. Wiersig, J. Review of exceptional point-based sensors. Photonics Res. 8 , 1457 (2020).
- 23. Huang,  R. et al. Exceptional  Photon  Blockade:  Engineering  Photon  Blockade  with  Chiral Exceptional Points. Laser Photonics Rev. 16 , 2100430 (2022).
- 24. Jin, B. et al. High-Performance Terahertz Sensing at Exceptional Points in a Bilayer Structure. Adv. Theory Simul. 1 , 1800070 (2018).
- 25. Park, J.-H. et al. Symmetry-breaking-induced plasmonic exceptional points and nanoscale sensing. Nat. Phys. 16 , 462-468 (2020).
- 26. Sakhdari, M., Farhat, M. &amp; Chen, P.-Y. PT-symmetric metasurfaces: wave manipulation and sensing

using singular points. New J. Phys. 19 , 065002 (2017).

- 27. Yuan, H. et  al. Non-Hermitian Topolectrical  Circuit  Sensor  with  High  Sensitivity. Adv.  Sci. 10 , 2301128 (2023).
- 28. Xu, J. et al. Experiment on parity-time-symmetry phase transition in non-Hermitian metasurfaces. J. Opt. Soc. Am. B 39 , 1847 (2022).
- 29. Li, Y. et  al. Bifunctional  sensing based on an exceptional point with bilayer metasurfaces. Opt. Express 31 , 492 (2023).
- 30. Chen, Y., Hu, J., Yin, S., Zhang, W. &amp; Huang, W. Bimodal Absorber Frequencies Shift Induced by the Coupling of Bright and Dark Modes. Materials 17 , 3379 (2024).
- 31.  Yin, S. et al. Phase mismatch induced suppression of eigenmode resonance in terahertz metasurfaces. Opt. Express 30 , 19176 (2022).
- 32. Huang,  W. et  al. Robust  and  broadband  integrated  terahertz  coupler  conducted  with  adiabatic following. New J. Phys. 21 , 113004 (2019).
- 33. Liu, T. et  al. Edge Detection Imaging by Quasi-Bound States in the Continuum. Nano Lett. 24 , 14466-14474 (2024).
- 34. Huang, W. et al. Inverse Design of Metamaterial Bound States in the Continuum Device via Coupled Mode Theory. IEEE Trans. Microw. Theory Tech. 72 , 1520-1528 (2024).
- 35. Yin, S. et al. Coupling-enabled chirality in terahertz metasurfaces. Nanophotonics 12 , 1317-1326 (2023).
- 36. Singh, S. et al. Advances in Metasurfaces: Topology, Chirality, Patterning, and Time Modulation.

IEEE Antennas Propag. Mag. 64 , 51-62 (2022).

- 37. Wu, Y. et al. Extremely asymmetric absorption and reflection near the exceptional point for threedimensional metamaterial. J. Phys. Appl. Phys. 57 , 365101 (2024).
- 38. Shi, S.-F. et al. Optimizing Broadband Terahertz Modulation with Hybrid Graphene/Metasurface Structures. Nano Lett. 15 , 372-377 (2015).
- 39. Melik-Gaykazyan,  E. et al. From  Fano  to  Quasi-BIC  Resonances  in  Individual  Dielectric Nanoantennas. Nano Lett. 21 , 1765-1771 (2021).
- 40. Liang, Y. et al. Ultra-broadband long-wave infrared metasurface absorber based on Peano fractal curve. Results Phys. 33 , 105169 (2022).
- 41. Ho, L., Pepper, M. &amp; Taday, P. Signatures and fingerprints. Nat. Photonics 2 , 541-543 (2008).
- 42. Shi, X. &amp; Han, Z. Enhanced terahertz fingerprint detection with ultrahigh sensitivity using the cavity defect modes. Sci. Rep. 7 , 13147 (2017).
- 43. Lee, S. H. et al. Detection and discrimination of SARS-CoV-2 spike protein-derived peptides using THz metamaterials. Biosens. Bioelectron. 202 , 113981 (2022).
- 44. Duggan, R., Mann, S. A. &amp; Alù, A. Limitations of Sensing at an Exceptional Point. ACS Photonics 9 , 1554-1566 (2022).
- 45. Lv, X., Lin, Q., Qiu, W., Guan, H. &amp; Lu, H. Non-Hermitian Optical Tunable System Based on Lithium Niobate Coupling Resonator. IEEE Photonics J. 14 , 1-5 (2022).
- 46. Hodaei,  H. et  al. Enhanced  sensitivity  at  higher-order  exceptional  points. Nature 548 ,  187-191 (2017).
- 47. Wang, H., Lai, Y.-H., Yuan, Z., Suh, M.-G. &amp; Vahala, K. Petermann-factor sensitivity limit near an exceptional point in a Brillouin ring laser gyroscope. Nat. Commun. 11 , 1610 (2020).
- 48. Rao, Z. et al. Braiding reflectionless states in non-Hermitian magnonics. Nat. Phys. (2024)
- 49. Yin, S. et al. Ruler equation for precisely tailoring the resonance frequency of terahertz U-shaped metamaterials. J. Opt. 21 , 025101 (2019).
- 50. Zhang, J. et al. Highly sensitive detection of malignant glioma cells using metamaterial-inspired THz  biosensor  based  on  electromagnetically  induced  transparency. Biosens.  Bioelectron. 185 , 113241 (2021).
- 51. Evlyukhin, A. B. et al. Detuned Electrical Dipoles for Plasmonic Sensing. Nano Lett. 10 ,  45714577 (2010).
- 52. Singh, R. et al. Ultrasensitive terahertz sensing with highQ Fano resonances in metasurfaces. Appl. Phys. Lett. 105 , 171101 (2014).
- 53. Wang, Y., Han, Z., Du, Y. &amp; Qin, J. Ultrasensitive terahertz sensing with highQ toroidal dipole resonance governed by bound states in the continuum in all-dielectric metasurface. Nanophotonics 10 , 1295-1307 (2021).
- 54. Beruete, M. &amp; Jáuregui-López, I. Terahertz Sensing Based on Metasurfaces. Adv. Opt. Mater. 8 , 1900721 (2020).
- 55. Sun, R., Li, W., Meng, T. &amp; Zhao, G. Design and optimization of terahertz metamaterial sensor with high sensing performance. Opt. Commun. 494 , 127051 (2021).
- 56. Xu, J. et al. Terahertz toroidal dipole metamaterial sensors for detection of aflatoxin B1. Chin. Phys.

B 33 , 038701 (2024).

## Declarations

Funding : This work was supported by:

Self-determined  research  funds  of  CCNU  from  the  colleges  basic  research  and  operation  of  MOE (CCNU23CG016)

Fundamental Research Funds for the Central Universities (CCNU24JCPT011)

National Natural Science Foundation of China (Grants No. 12364045, No. 12264028, and No. 12304420),

Natural Science Foundation of Jiangxi Province (Grants No. 20232BAB201040 and No. 20232BAB211025),

Young Elite Scientists Sponsorship Program by JXAST (Grant No. 2023QT11, No. 2025QT04).

## Author contributions:

Conceptualization: Hai Lin, Xintong Shi

Resources: Yun Shen, Hai Lin, Rongxin Tang

Investigation: Xintong Shi, Yanjie Wu, Le Li, Chenhui Zhao, Shouxin Duan

Visualization: Xintong Shi, Junyi Zhang, Shuyuan Xiao

Supervision: Hai Lin, Yun Shen, Shuyuan Xiao

Writing-original draft: Xintong Shi

Writing-review &amp; editing: Hai Lin, Tingting Liu, Shuyuan Xiao

Competing interests: All other authors declare they have no competing interests

Data and materials availability: All data are available in the main text or the supplementary materials.
Xintong Shi Hai Lin Tingting Liu Yun Shen Rongxin Tang Le Li Junyi Zhang Yanjie Wu Shouxin Duan Chenhui Zhao Shuyuan Xiao
Polarization-controlled non-Hermitian metasurfaces for ultra-sensitive terahertz sensing
Non-Hermitian systems offer significant advantages in sensor design,
especially at the exceptional points. However, the extreme sensitivity near
these points poses great challenges due to fabrication errors and system
noises, which degrade sensing performance. To address this, we introduce a
novel approach leveraging the polarization degrees of freedom in non-Hermitian
systems. In this work, we establish a direct relation between the incident
polarization and the transmission phase of a coupled metasurface system and
achieve the polarization-controlled phase singularity even post-fabrication.
The incident polarization angle can be utilized as a sensing index, which
enables indirect and accurate measurement. The theoretical approach is
experimentally validated using a general design of THz non-Hermitian
metasurface sensors. Our method enhances robustness and sensitivity, opening
new avenues for practical applications in ultra-sensitive sensing.