收藏切换
Ultrasensitive Quantum Sensors Based on High-Order Exceptional Bound States
收藏切换
PDF
Shaohui Liu, Tian Chen*, Deyuan Zou, Xiangdong Zhang*
Research. Vol 9 Article ID 1091
Less
收藏切换
Research. Vol 9 Article ID 1091
• RESEARCH ARTICLE •
Ultrasensitive Quantum Sensors Based on High-Order Exceptional Bound States
Full
Shaohui Liu, Tian Chen*, Deyuan Zou, Xiangdong Zhang*
Affiliations
  • Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements of Ministry of Education, School of Physics, Beijing Institute of Technology, Beijing 100081, China.
About Author:

† These authors contributed equally to this work.

Published: 2026-02-20 doi: 10.34133/research.1091
Outline
收藏切换

High-precision sensors are of fundamental importance in modern society and technology. Although there are various schemes for the construction of sensors relying on different physical mechanisms, obtaining sensors with higher levels of sensitivity and stronger robustness has always been expected. In particular, non-Hermitian quantum sensors have recently attracted substantial attention due to their unique properties. So far, 2 types of non-Hermitian sensors based on exceptional points and topological zero modes have been realized. Here, high-order exceptional bound states with robust properties are constructed for the first time. Based on these states, we propose theoretically and demonstrate experimentally another new type of non-Hermitian quantum sensors. Such sensors not only are robust against disorders but also have unprecedented sensitivity. Their sensing performance can display the improvement of many orders of magnitude over the previous non-Hermitian sensors. Furthermore, we design and fabricate such sensors based on circuit networks. Taking weak magnetic field detection as an example, we also experimentally demonstrate their sensing capabilities. Our work opens up new avenues for the development of highly sensitive sensors, which have a wide range of applications in various fields.

Shaohui Liu, Tian Chen, Deyuan Zou, Xiangdong Zhang. Ultrasensitive Quantum Sensors Based on High-Order Exceptional Bound States[J]. Research, 2026 , 9 (2) : 1091 - . DOI: 10.34133/research.1091
Physical properties in many real systems, such as the presence of loss and/or gain, should be described by non-Hermitian Hamiltonians [1]. Effective non-Hermitian Hamiltonians can induce novel topological phases and unique phenomena [2–7], which do not exist in the Hermitian system. For example, the exceptional points (EPs) in non-Hermitian systems can exhibit special spectral degeneracies, at which 2 or more eigenvalues are degenerated and the related eigenvectors coalesce. The rich physics at and around EPs has been elucidated in recent years, especially the various designs of sensors related to the EP [8–35]. While the EP sensors have shown remarkable potential, their practical implementation requires precise tuning to operate at the EP, and they face the inherent challenge of noise amplification. In parallel, non-Hermitian topological states provide a fundamentally different and more robust framework for non-Hermitian sensing.
The non-Hermitian skin effect (NHSE), namely, that the majority of eigenstates of a non-Hermitian operator are localized at boundaries, has also been discovered as another manifestation of the non-Hermitian topology. Broader implications of NHSE have been under investigation [36–48], and the effect has been observed in experiments [11,49–54]. Recent studies have shown that the energy shift of a midgap non-Hermitian topological state in response to boundary perturbations can be exponentially amplified with increasing system size [55]. This effect has been experimentally demonstrated using topolectrical circuits to detect minor variations in electrical parameters [53,56,57]. Although non-Hermitian topological sensors demonstrate promising application prospects, their detectable areas are rather localized, which are limited to the edges of systems in common [53–57].
Recently, a new class of robust states known as exceptional bound (EB) states, which are distinct from the well-known topological and non-Hermitian skin boundary states, has been found [58–60]. However, these revealed that states only correspond to the case with the second-order EPs. The question is whether there exist higher-order EB states corresponding to high-order EPs, which remains unknown. If the higher-order EB states exist, would they be helpful for achieving ultrasensitive sensing?
In this work, we first demonstrate that the high-order EB states, a new class of robust quantum states, can be constructed in exceptional systems. Compared to the second-order case, high-order EB states exhibit more robust against disorders. Especially, the high-order EB states with flat spatial distributions can be realized. Based on such states, we construct a novel class of quantum sensors to achieve ultrasensitive sensing across the entire surface of the system. Their sensing performances can display the improvement of many orders of magnitude when comparing with the previous non-Hermitian sensors. Furthermore, we design and fabricate corresponding circuit network and integrated circuit sensors to experimentally observe these phenomena. Taking weak magnetic field detection as an example, we experimentally demonstrate the sensing capabilities of such circuit sensors. Thus, our work opens up the exciting possibility for realizing robust sensors with unprecedented sensitivity.
Consider a generic Nth-order non-Hermitian Hamiltonian H, this system exhibits simultaneous coalescence of both eigenvalues and eigenvectors, characteristic of the Nth-order EP. The corresponding biorthogonal basis is constructed from the left and right eigenvectors of H, HψiR=ϵiψiR and H†ψiL=ϵi∗ψiL, where ϵi denotes the eigenvalue of H. These eigenvectors satisfy the biorthogonality condition ψiLψjR=δij. Starting from these left and right eigenvectors of the second-order Hamiltonian, EB states embodying the nonzero entanglement entropy have been reported [58–60]. Those EB states are special boundary modes arising from the asymmetric truncation of long-range correlations and defectiveness near an EP. Their corresponding truncated projector exhibits anomalous eigenvalues beyond the usual bounds of 0 and 1. Here, we provide the construction of EB states for any Nth-order Hamiltonian as
Pk=12I+Hϵ1=∣ψ1Rk〉〈ψ1Lk∣+∑i≠1N121+ϵiϵ1|ψiRk〉〈ψiLk|
where the first term in the right part is the same as the construction of EB states in the second-order Hamiltonian, and the biorthogonal projection operator for the high-order EP system is introduced by weighting 121+ϵiϵ1 with its corresponding eigenvalues.
The Hamiltonian P can be obtained through Fourier transform,
Pl1αl2β=∑k′Pαβke−ik′·l1−l2,
where the element refers to the coupling from the sublattice α of unit cell l1 to the sublattice β of unit cell l2. The value of α (β) is chosen within 1N for the N-order EB system. In addition, to avoid the singularity caused by the degenerate at the EP, we introduce a small offset ∆ in the momentum space as k′=k+∆ in the Fourier transform.
In Fig. 1A, a 2-dimensional (2D) real space is considered. The scale of system is expressed as L=Lx1×Lx2, where Lx1 (Lx2) is the scale along the x1x2 direction. The 2D momentum space is spanned by k=kx1kx2, which is related to x=x1x2 in 2D real space by the Fourier transform. Here, l1=x1l1x2l1 refers to the unit cell, and x1l1 [x2l1] represents the real space coordinate x1 (x2) of l1. The quantity l2 is the same as above. The orange balls in Fig. 1A represent unit cells, and each cell consists of N sublattice sites that are represented by gray balls as shown in Fig. 1B. The intracell couplings (a1 and b1) between those adjacent sublattices are nonreciprocal, which correspond to Pl1αl1α+1=a1 and Pl1αl1α−1=b1, respectively.
Figure 1C shows couplings along the x1 direction. The red solid line represents the nearest-neighbor (NN) intercell couplings that the sublattice sites are coupled with staggered sites by nonreciprocal couplings (a2 and b2). Here, Pl1αl2β=a2 and Pl2βl1α=b2 with x1l1+1=x1l2 and x2l1=x2l2 in l1 and l2. Similarly, other long-range couplings can be obtained by the similar rules. The nonreciprocal couplings are al=1Lx1∑kx1′fkx1−1/2eikx1′l−1 and bl=1Lx1∑kx1′fkx11/2eikx1′l−1. The term fkx1=21−coskx1B is taken, and B represents the mode when the system gradually reaches EP. The selection of fkx1 is that it possesses the most local Fourier coefficients [54]. Here, B=4 is taken, which represents the minimal exponent ensuring maximally identical coefficients al. In such a case, all of the al are assigned identical values that a1=a2=…=aLx1=a=12Lx1∆−4 when considering ∆ is a small offset, and for the couplings bl, b0=3,b1=−2,b2=12,b3=…=bLx1≈0. In Supplementary Information S1, we show the different forms of fkx1 and related coupling coefficients.
In Fig. 1D, the NN intercell couplings in the x2 direction are shown. The couplings are nonreciprocal (y0andy1) at the edge of the x2 direction, Pl1αl2α=y0, direction that which satisfies x2l1=1andx2l2=2,orx2l1=Lx2andx2l2=Lx2−1. The reciprocal couplingsare trivial at the bulk of the x2 Pl1αil2αi=y1wherex1l1=x1l2andx2l1+1=x2l2,orx2l1−1=x2l22≤x2l1≤Lx2−1. We engineer nonreciprocal couplings at the edge of x2direction to achieve a uniform coupling structure in real space along that direction. This prevents any undesired edge-induced modulation, thereby flattening the spatial distribution of the EB states and ensuring the flat spatial distribution across the entire 2D plane.
After the above construction, we further explore the properties of P. Setting Lx2=1yields a 1D Nth−order EB system (x1 direction chain in Fig. 1A). Numerical diagonalization of P gives eigenvalues and eigenstates for arbitrary N (details in Supplementary Information S2). Figure 1E shows the density of states (DOS) for the 1D fifth-order EB system with Lx1=300 and ∆=0.001. Besides the trivial states confined within [0, 1], 4 isolated EB eigenvalues emerge far outside this range, marked by red (first), blue (second), green (third), and purple (fourth). Each kind corresponds to multiple eigenvalues and eigenstates as shown in Fig. 1F and G. In Fig. 1F, colored dots (e.g., 8 dots for the first order) match the DOS peaks in Fig. 1E. Notably, higher-order EB eigenvalues (E) exhibit larger separations, with the first-order EB states being the most isolated—enhancing their robustness, unlike lower-order states. Figure 1G displays the whole state distributions of the EB system with respect to unit cell index. Remarkably, the first EB state exhibits a well-defined flat spatial distribution.
For Lx2≠1, the system becomes 2D EB system with eigenvalues and eigenstates obtained by numerical diagonalization (details in Supplementary Information S2). As shown in Fig. 1H with Lx1=60,Lx2=10, and ∆=0.001, the 2D system exhibits the characteristic similar to the 1D case, with the first EB states showing flat spatial profiles (red) versus the heavily fluctuating trivial states (gray).
Building upon these 1D and 2D results, we emphasize that our construction method can generalize to arbitrary dimensions. Using eigenvalues and eigenstates of Nth−order EP system H0, the approach remains valid via the Fourier transform in higher-dimensional spaces (Supplementary Information S3).
We demonstrate the robustness of higher-order EB states following the EB spectral flow under spatial truncation [41]. The real-space cut can lead to a new non-Hermitian Pcut=RPR, where R indicates the cut operation. Figure 2A shows the calculated results for the 1D EB system Pcut, where the absolute eigenvalues Ecut of Pcut evolve as the truncation cut xcut varying from 2 to 300 (we fix the untruncated region xLxR=1xcut). This EB spectral flow is reminiscent of the spectral flow of topological states as flux is threaded. The existence of this spectral flow implies that the EB state that is well-defined at a large truncated system position xcut is adiabatically connected to a qualitatively similar EB state at a very small truncated system with xcut∝O1 when the size of lattice approximates a continuum. Three red curves in Fig. 2A means the change of 8 eigenvalues of the first EB states with the cut position. Even if there are 3 different absolute values for these first EB states, they are all distributed flat on the spatial positions. In the right part of Fig. 2A, the spatial lattice profiles of the first EB states are shown at xcut=100 and 200, respectively. When comparing these 2 distributions, it is found that the same flat spatial profile is maintained, which means that the first EB states have good robustness against spatial truncation.
The robustness of these EB states can be further demonstrated by analyzing the response of system against hopping disorder. We introduce disorder by multiplying all coupling terms by a random factor 1+rR, where rR is a random variable uniformly distributed over the interval −μμ. In this way, the element Pij changes to Pij1+rR. The range of disorder μ is chosen to be μ=0.01,0.02,0.05, and 0.1. As shown in Fig. 2B, at zero disorder, the eigenvalues E of the first and second EB states are plotted by red and blue dots. The details about the third and fourth EB states and the bulk states are indicated in the black box (top right of Fig. 2B). It is shown that with the increase of μ, the fluctuation near the initial eigenvalue also increases. However, it is obvious that the eigenvalues of the first and second EB states are separated from the other eigenvalues, and the first EB state is further separated. The corresponding details of the first and second EB states are shown in the middle right and bottom right of Fig. 2B. It can be seen that the first and second EB states are still robust even though the largest disorder exists (μ=0.1).
The similar robustness can also be found for the 2D EB non-Hermitian system. Details are shown in Supplementary Information S4. By considering that these flat distribution EB states have good robustness against spatial truncation and hopping disorder, we can construct ultrasensitive robust quantum sensors.
In the following, we analyze the sensitivity of the constructed EB systems to external perturbations. For the 1D system, we consider the measurand Γ shown in Fig. 3A. When considering the connection between the first sublattice site of unit-cell l1 (indexed by state l11) and the last sublattice site of unit-cell l2 (indexed by state l2N), the Hamiltonian of the external perturbations can be expressed as ∆Pl1,l2=Γl11〈l2N+l2N〈l11. Thus, the measurand of connecting all the unit cells of the 1D chain with perturbation strength Γ in Fig. 3A is expressed as ∆P1D=∑l1∑l2≠l1∆Pl1,l2. By performing numerical diagonalization, we obtain the eigenvalues EΓ of the first EB states and the shift ∆E=EΓ−E0, where E0 is the eigenvalues of the first EB states without perturbations.
In Fig. 3B, we plot numerical results of ∆E as a function of Γ with various a, which are represented by the circles. It is found that the curves of the relation between ∆E and Γ exhibit 2 distinct regimes. To better analyze the changing characteristics of ∆E, we perform systematic fitting for these numerical results. In fact, ∆E depends on several variables, including Γ, a, N, and L. Thus, through numerical calculations, we first take determined values for a and L, and discuss the functional relation of ∆E on Γ by numerical calculations. It is found that ∆E∝Γ at small ∆E regime and a power-law scaling ∆E∝Γ1/N at large ∆E regime. Then, with Γ and L fixed, we observe ∆E∝aN−1/2 at small ∆E and ∆E∝aN−1/N at large ∆E. Finally, keeping Γ and a constant, we obtain the size scaling ∆E∝Lθ1N−1+1, where θ1 is determined by all of the nonreciprocal couplings and θ1≈0.66 provides excellent fitting results. These findings lead us to the comprehensive scaling relations:
∆E∝ΓaN−1/2Lx1θ11−N+1small∆EΓ1/NaN−1/NLx1θ11−N+1large∆E.
The further detailed numerical fitting procedures and the first-order effect using perturbation theory are provided in Supplementary Information S5. The results in Eq. (3) also clearly demonstrates that the behavior of ∆E is governed by a and N. This can be attributed to 2 key factors: First, a can take significantly large values since a=12Lx1∆−4, where ∆ is a small offset. Second, the exponential dependence on N further amplifies the effect of a. For example, when taking Lx1=10 and Γ=10−15 with the same measurand via couplings of connecting all the unit cells, the shift of EB system is ∆E≈10−2, and the result of non-Hermitian Su–Schrieffer–Heeger (SSH) sensor [55] is ∆E≈10−14. Therefore, our sensing performance can display the improvement of 12 orders of magnitude when comparing with the non-Hermitian SSH sensor. When taking Lx1=20andΓ=10−15, our sensing performance can give the 10 orders of magnitude improvements.
The solid lines and dashed lines in Fig. 3B represent the fitted curves for 2 regimes, respectively. The agreements between numerical results and fitting curves are very well for a wide range of Γ. Furthermore, the sensitivity S=∂∆E∂Γ of the system can be calculated by the fitting expressions of ∆E. Compared with the recent proposed non-Hermitian SSH sensor [55], in which only edge states are used as the sensor unit, our proposed sensors can exhibit many advantages in sensing objects because all unit cells in the structure can be used as the sensor unit.
Similarly, we consider the measurand of connecting all the unit cells of the 2D EB system with perturbation strength Γ, which has been shown in Fig. 3C. Following the same fitting procedures, we can also establish the fitting relations that ∆E∝Γ,∆E∝aN−1/2, and ∆E∝Lθ21−N+1=Lx1×Lx2θ21−N+1 at small ∆E regime, while ∆E∝Γ1/N, ∆E∝aN−1/N, and ∆E∝Lθ21−N+1 at large ∆E regime. The results also lead us to the comprehensive scaling relations:
∆E∝ΓaN−1/2Lθ21−N+1small∆EΓ1/NaN−1/NLθ21−N+1large∆E,
where θ2 represents the parameter for the 2D case. As revealed in Fig. 3D, good agreements between the fitting results and the numerical results are observed for the energy shift ∆E with perturbation Γ with N=5. Further details of the fitting process are described in Supplementary Information S5. The results are shown in Fig. 3D, where the circles denote the numerical results and lines denote the fitted results. Because all unit cells in the 2D structure can be used as the sensor unit, our system demonstrates superior performance ability. For example, when we take a perturbation Γ=10−15, the corresponding ∆E≈106 and S≈1021 with Lx1=100, Lx2=8, and a=106. This means that a very small perturbations can cause significant change of ∆E. That is to say, even minimal changes can be perceived by such a sensor.
Recent studies have shown that the signal-to-noise ratio (SNR) performance of EP-based sensors remains a controversial issue [61–65]. In these works, the SNR is defined as the standard deviation of the measured frequency shift, which is also employed in our study. The SNR at the EP is proportional to perturbations, which limits the detection accuracy enhancement of EP-based sensors. In contrast, our EB system exhibits stable SNR performance. This stability originates from the system’s operation at the EB state, where the energy is isolated and insensitive to common noise sources. The detailed analysis is provided in Supplementary Information S6, where we simulate the SNR under structural disorder and Gaussian noise, and show that the SNR can be stable even for strong nonreciprocal couplings, validating the robustness of the measurement.
Now, we explore how to experimentally construct the sensors mentioned above. Recent investigations have shown that classical electric circuits can be used to simulate various topological physics [49,52,53,56,59,66–72]. In the following, we explore the design of circuit networks to realize such sensors.
Our designed circuit sensor for the 2D fifth-order Hamiltonian is schematically shown in Fig. 4A. The circuit in Fig. 4A includes orange spheres with Lx1×Lx2=4×4, which correspond to the unit cells in Fig. 1A. Each orange sphere includes 5 nodes (gray spheres in Fig. 4B) corresponding to the sites within the unit cells. The connections between 5 nodes for the nonreciprocal couplings are achieved by capacitors (CU and CA) and buffer, which has been shown in Fig. 4B. The purpose of the grounding inductances LU and LA is to eliminate the on-site energy offsets in the system. However, it should be noted that the offset at each lattice site is distinct, so these differing offsets cannot be fully eliminated merely by treating them as a translated spectrum.
For the above circuit of unit cell, it can be described by the Laplacian J1ω, where the elements are Jl1αil1αi+1=−iωCU and Jl1αil1αi−1=−iωCA, and satisfy CU:CA=a1:b1. Correspondingly, appropriate values of the grounding components are also selected to make CU=ω02LU−1 and CA=ω02LA−1, with ω0 denoting the resonant frequency. Their specific nonreciprocal circuit connection diagrams for the NN, second NN, and third NN couplings are shown in Fig. 4C. The components are taken as CB:ω2LC−1:CD=b2:b3:b4. As shown in Fig. 4D, the nonreciprocal coupling along the x2 direction is realized by the combination of 2 capacitors CY and one buffer, and the other reciprocal couplings is performed by the capacitor CY. The Laplacian Jω for the circuit network above corresponds to the Hamiltonian P. Detailed demonstration is given in Topolectric circuit design. The implementation for the circuit is provided in Supplementary Information S7.
Now, we explore the sensor effect of the circuit by tuning the resistor R as the perturbation, ∆J=R−1∑l1∑l2≠l1l1α1〈l2αN+l2αN〈l1α1, and the impedance of the resistor 1/R corresponds to the perturbation strength Γ. We take the connection between the first node of the l1 unit (indexed by the state l1α1) and the last node of the l2 unit (indexed by the state l2αN). By connecting all the units with resistors, we observe the changes of the eigenvalues with the perturbation, which can be directly observed by measuring the impedance response [66]. The red line in Fig. 4E shows the frequency of the EB flat state without the perturbation. The other colored lines represent the experimental measured spectra of the impedance response with various perturbations. It is observed clearly that the shift of resonant frequency corresponds to the peak impedance shifts. Details are given in The impedance measurement on the circuit.
The black dots in Fig. 4F show the measured results for the frequency shifts of impedance peaks ∆f as a function of the perturbation (1/R), while the blue line is the fitted result. The agreements between the experimental and fitted results are very well. With the choice of electric components, the frequency shift of impedance satisfies with the perturbation as ∆f∝1/R1/5, which is the same to the results in the lattice model in the large eigenvalue-shift case. In Fig. 4G, we provide the sensitivity results based on the resonance frequency. The black dots represent the experimental results, while the orange line is the fitted sensitivity that is described as Scir∝1/R4/5. It is also the same to the results in the lattice model. Furthermore, the intrinsic robustness of EB states consequently ensures the stable sensing performance, with detailed analysis provided in Supplementary Information S8. In addition, due to the intrinsic robustness of the long-range coupling design and effective suppression of circuit noise through precision components and active compensation, the SNR of the designed circuit remains very stable. Detailed noise analysis is provided in Supplementary Information S9.
As a practical application example, now we explore how to use our designed circuit sensor to measure magnetic fields. In general, the magnetic field can have an impact on a region, and therefore, the signals on the region can be accumulated to be detected. Not like the non-Hermitian SSH chain sensor [55] whose working areas are just 2 ends of chain, our proposed EB circuit sensor can capture the signal on the whole surface of circuit network. Figure 5A shows the schematic diagram of the experimental device, including the magnetic field generation, the magnetoresistance chip, the designed EB circuit, and the impedance measurement part.
The weak magnetic field to be measured is generated by the Helmholtz coil, and the required strength of the weak magnetic field can be accurately controlled by adjusting the high-precision current source connected to the coil. The magnetoresistance chip is made of a thin layer of nickel–iron film and arranged as a resistor. In our experiment, many magnetoresistance chips serve as the perturbation resistors. After applying different external magnetic fields B,the resistanceR also changes in a proportional relationship due to the magnetoresistance effect, realizing the transformation from the measured physical quantity (magnetic field strength) to the circuit quantity. In addition, in order to eliminate the influence of the geomagnetic field, we compensate by biasing the current to generate a reverse magnetic field. Then, due to the high sensitivity and strong amplification effect of our sensor circuit, the weak change of the resistance causes a significant shift of the impedance peak in the circuit, which corresponds to the significant shift of the eigenvalue of the EB state under the perturbation. Finally, the shift of the impedance peak can be measured by connecting the EB circuit and impedance analyzer through positive and negative electrodes. Details are provided in Supplementary Information S10.
Figure 5B displays the impedance spectra at the first node of the first unit cell under magnetic fields of the magnetic field B=10,102,103,and 104nT (purple, yellow, orange, and blue lines, respectively). The rationale for selecting this node is that the distribution of the EB states at this location is flat, which accurately reflects the flat characteristics of the EB states. It can be clearly seen that the impedance spectrum measured under different external magnetic fields varies with the frequency, and the specific relationship between ∆f and B is shown in Fig. 5C. Figure 5D also shows the results of relationship between sensitivity and B. It is found that the experimental results based on the magnetic field are in good agreement with the theoretical results. Therefore, the EB sensor circuit we designed has been successfully applied to the ultrasensitive detection of weak magnetic fields. Little difference between points and lines comes from the external interference magnetic fields and loss of electric components. Moreover, the noise of resistor is also considered. In Supplementary Information S9, we have studied the effect of different levels of noise, the robustness of the sensing mechanism ensures performance stability under different fluctuations.
At end, we would like to emphasize that the designed circuit sensors are based on printed circuit board (PCB) platforms. In fact, the corresponding on-chip ultrasensitive sensing systems can also be realized by using a complementary metal-oxide semiconductor (CMOS) process technology [53,73]. The chips using CMOS technology can show the advantages as it is easily configurable, the parasitic effect is reduced enormously, and the operating frequency of the mode can also be increased. Thus, the sensitivity of circuit sensors based on CMOS technology can be improved by orders of magnitude compared to those based on PCB platform. That is, the detection limit of the magnetic field in Fig. 5 is nT. As a comparison, when the CMOS technology is used, weak magnetic signals as small as fT are expected to be detected.
Here, we have demonstrated magnetic field sensing using the magnetoresistor. Indeed, the core sensing mechanism is not limited to magnetic fields, and the same platform can be directly adapted to sense any physical quantity convertible into a change in circuit, such as displacement, angle, or liquid level, by connecting corresponding capacitive transducers. In addition, our results reveal that the circuit serves as a powerful and flexible platform for realizing the targeted quantum Hamiltonian. This enables precise access to the relevant eigenstates of the system, thereby facilitating experimental demonstration of the sensing enhancement predicted by quantum theory.
The sensing mechanism proposed in this work can be further analyzed by comparing it with other advanced sensing schemes, particularly EP-based sensors, topological quantum sensors, and conventional quantum sensors. EP-based sensors exploit the spectral singularity at non-Hermitian degeneracies, and the frequency splitting for an nth−order EP is generically proportional to the nth root of perturbation strength Γ. However, their practical utility is hindered by the amplification of noise due to extreme non-orthogonality of modes (divergent Petermann factor) [9,10]. For the non-Hermitian topological chain sensors, the energy shift of a boundary mode scales exponentially with system size. However, its exponential sensing performance often requires the measurand connecting between the 2 ends of chain [55]. For the conventional quantum sensor [74], although it shows unprecedented sensitivity and spatial resolution, the decoherence time is very short, which leads to increased control difficulty. As a comparison, our proposed sensor shows the stable SNR, and the detected areas of sensor are not limited to the ends of system. In addition, our sensor does not suffer from the decoherence and is able to continue working for long times.
In this work, we have theoretically disclosed a new class of robust quantum states and higher-order EB states in exceptional systems. These quantum states can be engineered by tuning the nonreciprocal couplings of the system. Especially, we can design quantum states with flat spatial distributions for both 1D and 2D systems. Compared to the second-order EB states, our designed high-order EB states exhibit more robust against disorders. Based on such states, we have constructed a novel class of quantum sensors. We have also demonstrated that such sensors not only have unprecedented sensitivity but also are robust against disorders. They have been proven to have significant advantages compared to existing quantum sensors. Furthermore, we have found that higher-order EB states can also exist in classical systems. We have designed and fabricated corresponding circuit network and integrated circuit sensors to experimentally observe the phenomena. Taking magnetic field detection as an example, we have demonstrated the effectiveness of our designed circuit sensor. Therefore, our proposed new quantum sensing scheme is expected to have widely potential applications in various fields.
Topological circuits are an emerging synthetic topological material platform that allows any tight-binding model of any coupling range or dimension to be realized. In such circuits, the behavior of the system is governed by fundamental electrical laws, among which Kirchhoff's law plays a central role. Specifically, based on Kirchhoff’s law, the voltage and current at the node have the relation:
Im=∑n≠mYmnVm−Vn+wmVm,
where Im and Vm are the input current and electrical voltage at node m. According to the conservation of current, Ix is equal to the total current flowing out of node m toward all other nodes n linked by the admittance Ymn, plus the current flowing from node m to the ground through the path with the admittance wm. For a single capacitive, resistive, and inductive component, the admittance are taken as iωC,1/R,and1/iωL, respectively.
In the circuit, the critical design consideration involves the realization of nonreciprocal coupling mechanisms. Here, the fundamental nonreciprocal connection is shown in Fig. 4B. The black triangle represents the buffer, which play the role of voltage follower. By connecting capacitor CU in series with the buffer, the other capacitance CA in parallel is added. Based on the characteristic that buffer can block the input current while keeping the output voltage stable, we get the following equation by carrying out Kirchhoff’s law on the circuit node 1 (left) and node 2 (right):
I1=iωCAV1−V2I2=iωCAV2−V1+iωCUV2−V1.
Equation (6) can be reexpressed in a matrix form as:
I1I2=iωCA−iωCAiωCU+CAiωCU+CAV1V2.
It can be seen that the conductance matrix is a non-Hermitian matrix, where 2 off-diagonal elements of CA and CU+CA represent effective values of the connecting capacitor between 2 nodes, and the nonreciprocal capacitance with the value being CA and CU+CA can be achieved.
Similarly, other nonreciprocals can be realized in circuit and schematically depicted in Fig. 4C and D. They are all realized by buffers, capacitors, and inductors, and the corresponding grounding components are also different. In the x2 direction, the circuits for realizing 2 kinds of couplings are shown in Fig. 4D. The nonreciprocal coupling is realized by the combination of 2 capacitors CY and one buffer, and the other reciprocal couplings is performed by the capacitor CY. For such a circuit network above, we can theoretically derive its circuit Laplacian Jω. In our study, the currents flowing into each node of circuit can be written as:
I1=∑j=03iωCUV1−V2+5j+4iωCYV21I2=∑j=03iωCUV2−V3+5j+iωCAV2−V1+1iωLBV2−V6+iωCCV2−V11+4iωLU+1iωLA+iωCB+1iωLCV2+iωCYV22I3=∑j=03iωCUV3−V4+5j+iωCAV3−V2+1iωLBV3−V7+iωCCV3−V12+4iωLU+1iωLA+iωCB+1iωLCV3+iωCYV23…I80=iωCCV80−V69+1iωLBV80−V74+iωCAV80−V79+1iωLA+iωCB+1iωLCV80+2iωCYV60
Equation (8) can be recast into a matrix form:
I1I2I3⋮I80=iωY1−CU0⋯0−CAY2−CU⋯00−CAY3⋯0⋮⋮⋮⋱⋮000⋯Y80V1V2V3⋮V80,
where Y1=4CU−4ω2LU,Y2=4CU+CA−1ω2LB+CC−4ω2LU−1ω2LA+CB−1ω2LC,…,Y80=CA−1ω2LB+CC+1iωLA+iωCB+1iωLC are the diagonal terms that are determined by all the components connected to the corresponding node and the grounding components. The circuit Laplacian can be obtained as:
Jω=iωY1−CU0⋯0−CAY2−CU⋯00−CAY3⋯0⋮⋮⋮⋱⋮000⋯Y80.
If the parameters are chosen as CU:CA:−ω02LB−1:CC:CY=a:b1:b2:b3:y0,ω=ω0 and CU=1ω02LU,CA=1ω02LA,CB=1ω02LB,CC=1ω02LC, i.e. Y1=Y2=…=Y80=0, the circuit Laplacian Jω has a good correspondence to the lattice model Hamiltonian at the resonance frequency ω0 that Jω0∝P.
In the circuit, we usually observe the characteristics of the system through the impedance response Zabω=Va−VbI, which is the ratio of voltage between nodes a and b due to the magnitude of the current that enters at a and leaves at b. Here, the 2-point impedance of node a and b reads:
Zab=∑nφn,a−φn,b2Ecir,n,
where Ecir,n and φn,a−φn,b denote the eigenvalue and the difference between the amplitudes of the nth admittance eigenmode of circuit Laplacian J. Based on the above expression of the impedance, it can be observed that Zab becomes significantly large when there exists a finite density of nontrivial eigenmodes with small eigenvalues Ecir,n. For each eigenvalue Ecir,nω, there corresponds a resonant frequency ωn~ at which Ecir,nωn~ attains a very small value ε. Consequently, this gives rise to a dominant term in the impedance, leading to the approximate relation Zωn~≈1εφn,12. Given the existence of multiple eigenvalues Ecir,n, a series of distinct resonant frequencies ωn~ can be identified. Furthermore, any shift in the eigenvalues directly corresponds to a shift in these resonant frequencies. Therefore, for the resonance frequency ω0 of the circuit Laplacian J without perturbation, there is a peak in impedance spectrum Zω, and the frequency shifts of this impedance peak represent the energy shift of the first EB mode.
1.
Moiseyev N. Non-hermitian quantum mechanics. Cambridge: Cambridge University Press; 2011.
2.
Feng L, El-Ganainy R, Ge L. Non-hermitian photonics based on parity-time symmetry. Nat Photonics. 2017;11(12):752–762.
3.
El-Ganainy R, Makris KG, Khajavikhan M, Musslimani ZH, Rotter S, Christodoulides DN. Non-hermitian physics and pt symmetry. Nat Phys. 2018;14(1):11–19.
4.
Miri MA, Alù A. Exceptional points in optics and photonics. Science. 2019;363(6422):Article eaar7709.
5.
Ozdemir SK, Rotter S, Nori F, Yang L. Parity-time symmetry and exceptional points in photonics. Nat Mater. 2019;18(8):783–798.
6.
Bergholtz EJ, Budich JC, Kunst FK. Exceptional topology of non-hermitian systems. Rev Mod Phys. 2021;93(1):Article 015005.
7.
Ding K, Fang C, Ma GC. Non-hermitian topology and exceptional-point geometries. Nat Rev Phys. 2022;4(12):745–760.
8.
Wiersig J. Enhancing the sensitivity of frequency and energy splitting detection by using exceptional points: Application to microcavity sensors for single-particle detection. Phys Rev Lett. 2014;112(20):Article 203901.
9.
Chen W, Kaya Ozdemir S, Zhao G, Wiersig J, Yang L. Exceptional points enhance sensing in an optical microcavity. Nature. 2017;548(7666):192–196.
10.
Hodaei H, Hassan AU, Wittek S, Garcia-Gracia H, El-Ganainy R, Christodoulides DN, Khajavikhan M. Enhanced sensitivity at higher-order exceptional points. Nature. 2017;548(7666):187–191.
11.
Xiao L, Zhan X, Bian ZH, Wang KK, Zhang X, Wang XP, Li J, Mochizuki K, Kim D, Kawakami N, et al.Observation of topological edge states in parity-time-symmetric quantum walks. Nat Phys. 2017;13(11):1117–1123.
12.
Chen P-Y, Sakhdari M, Hajizadegan M, Cui Q, Cheng MM-C, El-Ganainy R, Alù A. Generalized parity–time symmetry condition for enhanced sensor telemetry. Nat Electron. 2018;1(5):297–304.
13.
Zhao H, Chen Z, Zhao R, Feng L. Exceptional point engineered glass slide for microscopic thermal mapping. Nat Commun. 2018;9(1):1764.
14.
Chen PY, El-Ganainy R. Exceptional points enhance wireless readout. Nat Electron. 2019;2(8):323–324.
15.
Dong ZY, Li ZP, Yang FY, Qiu CW, Ho JS. Sensitive readout of implantable microsensors using a wireless system locked to an exceptional point. Nat Electron. 2019;2(8):335–342.
16.
Hokmabadi MP, Schumer A, Christodoulides DN, Khajavikhan M. Non-hermitian ring laser gyroscopes with enhanced sagnac sensitivity. Nature. 2019;576(7785):70–74.
17.
Wang S, Hou B, Lu W, Chen Y, Zhang ZQ, Chan CT. Arbitrary order exceptional point induced by photonic spin-orbit interaction in coupled resonators. Nat Commun. 2019;10(1):832.
18.
Xiao Z, Li H, Kottos T, Alù A. Enhanced sensing and nondegraded thermal noise performance based on pt-symmetric electronic circuits with a sixth-order exceptional point. Phys Rev Lett. 2019;123(21):Article 213901.
19.
Zhang M, Sweeney W, Hsu CW, Yang L, Stone AD, Jiang L. Quantum noise theory of exceptional point amplifying sensors. Phys Rev Lett. 2019;123(18):Article 180501.
20.
Park J-H, Ndao A, Cai W, Hsu L, Kodigala A, Lepetit T, Lo Y-H, Kanté B. Symmetry-breaking-induced plasmonic exceptional points and nanoscale sensing. Nat Phys. 2020;16(4):462–468.
21.
Xu H, Mason D, Jiang L, Harris JG. Topological energy transfer in an optomechanical system with exceptional points. Nature. 2016;537(7618):80–83.
22.
Hassan AU, Zhen B, Soljacic M, Khajavikhan M, Christodoulides DN. Dynamically encircling exceptional points: Exact evolution and polarization state conversion. Phys Rev Lett. 2017;118(9):Article 093002.
23.
Hu W, Wang H, Shum PP, Chong YD. Exceptional points in a non-hermitian topological pump. Phys Rev B. 2017;95(18):Article 184306.
24.
Doppler J, Mailybaev AA, Bohm J, Kuhl U, Girschik A, Libisch F, Milburn TJ, Rabl P, Moiseyev N, Rotter S. Dynamically encircling an exceptional point for asymmetric mode switching. Nature. 2016;537(7618):76–79.
25.
Yoon JW, Choi Y, Hahn C, Kim G, Song SH, Yang KY, Lee JY, Kim Y, Lee CS, Shin JK, et al.Time-asymmetric loop around an exceptional point over the full optical communications band. Nature. 2018;562(7725):86–90.
26.
Zhang XL, Wang SB, Hou B, Chan CT. Dynamically encircling exceptional points: In situ control of encircling loops and the role of the starting point. Phys Rev X. 2018;8(2):Article 021066.
27.
Naghiloo M, Abbasi M, Joglekar YN, Murch KW. Quantum state tomography across the exceptional point in a single dissipative qubit. Nat Phys. 2019;15(12):1232–1236.
28.
Liu W, Wu Y, Duan CK, Rong X, Du J. Dynamically encircling an exceptional point in a real quantum system. Phys Rev Lett. 2021;126(17):Article 170506.
29.
Abbasi M, Chen W, Naghiloo M, Joglekar YN, Murch KW. Topological quantum state control through exceptional-point proximity. Phys Rev Lett. 2022;128(16):Article 160401.
30.
Nasari H, Lopez-Galmiche G, Lopez-Aviles HE, Schumer A, Hassan AU, Zhong Q, Rotter S, LiKamWa P, Christodoulides DN, Khajavikhan M. Observation of chiral state transfer without encircling an exceptional point. Nature. 2022;605(7909):256–261.
31.
Schumer A, Liu YGN, Leshin J, Ding L, Alahmadi Y, Hassan AU, Nasari H, Rotter S, Christodoulides DN, LiKamWa P, et al.Topological modes in a laser cavity through exceptional state transfer. Science. 2022;375(6583):884–888.
32.
Feng L, Wong ZJ, Ma RM, Wang Y, Zhang X. Single-mode laser by parity-time symmetry breaking. Science. 2014;346(6212):972–975.
33.
Ding K, Ma GC, Xiao M, Zhang ZQ, Chan CT. Emergence, coalescence, and topological properties of multiple exceptional points and their experimental realization. Phys Rev X. 2016;6(2):Article 021007.
34.
Bai K, Lin C, Liu T, Li J-Z, Lyu X, Xiao M. Nonlinear chiral-like state transfer realized with a minimal set of parameters. Nat Commun. 2025;16(1):5844.
35.
Fang L, Bai K, Guo C, Liu T-R, Li J-Z, Xiao M. Exceptional features in nonlinear hermitian systems. Phys Rev B. 2025;111(16):L161102.
36.
Lee TE. Anomalous edge state in a non-hermitian lattice. Phys Rev Lett. 2016;116(13):Article 133903.
37.
Li Y, Peng YG, Han L, Miri MA, Li W, Xiao M, Zhu XF, Zhao J, Alu A, Fan S, et al.Anti-parity-time symmetry in diffusive systems. Science. 2019;364(6436):170–173.
38.
Zhang K, Yang Z, Fang C. Correspondence between winding numbers and skin modes in non-hermitian systems. Phys Rev Lett. 2020;125(12):Article 126402.
39.
Yao S, Wang Z. Edge states and topological invariants of non-hermitian systems. Phys Rev Lett. 2018;121(8):Article 086803.
40.
Palacios LS, Tchoumakov S, Guix M, Pagonabarraga I, Sanchez S, Grushin AG. Guided accumulation of active particles by topological design of a second-order skin effect. Nat Commun. 2021;12(1):4691.
41.
Borgnia DS, Kruchkov AJ, Slager RJ. Non-hermitian boundary modes and topology. Phys Rev Lett. 2020;124(5):Article 056802.
42.
Okuma N, Kawabata K, Shiozaki K, Sato M. Topological origin of non-hermitian skin effects. Phys Rev Lett. 2020;124(8):Article 086801.
43.
Gong ZP, Ashida Y, Kawabata K, Takasan K, Higashikawa S, Ueda M. Topological phases of non-hermitian systems. Phys Rev X. 2018;8(3):Article 031079.
44.
Kunst FK, Edvardsson E, Budich JC, Bergholtz EJ. Biorthogonal bulk-boundary correspondence in non-hermitian systems. Phys Rev Lett. 2018;121(2):Article 026808.
45.
Yokomizo K, Murakami S. Non-bloch band theory of non-hermitian systems. Phys Rev Lett. 2019;123(6):Article 066404.
46.
Kawabata K, Shiozaki K, Ueda M, Sato M. Symmetry and topology in non-hermitian physics. Phys Rev X. 2019;9(4):Article 041015.
47.
Wang HY, Song F, Wang Z. Amoeba formulation of non-bloch band theory in arbitrary dimensions. Phys Rev X. 2024;14(2):Article 021011.
48.
Kawabata K, Numasawa T, Ryu S. Entanglement phase transition induced by the non-hermitian skin effect. Phys Rev X. 2023;13(2):Article 021007.
49.
Zou D, Chen T, He W, Bao J, Lee CH, Sun H, Zhang X. Observation of hybrid higher-order skin-topological effect in non-hermitian topolectrical circuits. Nat Commun. 2021;12(1):7201.
50.
Ghatak A, Brandenbourger M, van Wezel J, Coulais C. Observation of non-hermitian topology and its bulk-edge correspondence in an active mechanical metamaterial. Proc Natl Acad Sci USA. 2020;117(47):29561–29568.
51.
Wang B, Chen T, Zhang XD. Observation of novel robust edge states in dissipative non-hermitian quantum walks. Laser Photonics Rev. 2020;14(7):2000092.
52.
Helbig T, Hofmann T, Imhof S, Abdelghany M, Kiessling T, Molenkamp LW, Lee CH, Szameit A, Greiter M, Thomale R. Generalized bulk–boundary correspondence in non-hermitian topolectrical circuits. Nat Phys. 2020;16(7):747–750.
53.
Deng W, Zhu W, Chen T, Sun H, Zhang X. Ultrasensitive integrated circuit sensors based on high-order non-hermitian topological physics. Sci Adv. 2024;10(38):eadp6905.
54.
Xie X, Ma F, Rui WB, Dong Z, Du Y, Xie W, Zhao YX, Chen H, Gao F, Xue H. Non-hermitian dirac cones with valley-dependent lifetimes. Nat Commun. 2025;16(1):1627.
55.
Budich JC, Bergholtz EJ. Non-hermitian topological sensors. Phys Rev Lett. 2020;125(18):Article 180403.
56.
Yuan H, Zhang WX, Zhou ZL, Wang WL, Pan NQ, Feng Y, Sun HJ, Zhang XD. Non-hermitian topolectrical circuit sensor with high sensitivity. Adv Sci. 2023;10(19):2301128.
57.
Zhou X, Zhang W, Cao W, Zhang X. Non-hermitian floquet topological sensors for ultrasensitive detection of dynamic signals. Phys Rev Lett. 2025;135(10):Article 106601.
58.
Lee CH. Exceptional bound states and negative entanglement entropy. Phys Rev Lett. 2022;128(1):Article 010402.
59.
Zou D, Chen T, Meng H, Ang YS, Zhang X, Lee CH. Experimental observation of exceptional bound states in a classical circuit network. Sci Bull. 2024;69(14):2194–2204.
60.
Yang M, Lee CH. Beyond the non-hermitian skin effect: Scaling-controlled topology from exceptional-bound bands. arXiv. 2025. https://doi.org/10.48550/arXiv.2510.06338.
61.
Lai YH, Lu YK, Suh MG, Yuan Z, Vahala K. Observation of the exceptional-point-enhanced sagnac effect. Nature. 2019;576(7785):65–69.
62.
Lau HK, Clerk AA. Fundamental limits and non-reciprocal approaches in non-hermitian quantum sensing. Nat Commun. 2018;9(1):4320.
63.
Wang H, Lai YH, Yuan Z, Suh MG, Vahala K. Petermann-factor sensitivity limit near an exceptional point in a brillouin ring laser gyroscope. Nat Commun. 2020;11(1):1610.
64.
Wiersig J. Prospects and fundamental limits in exceptional point-based sensing. Nat Commun. 2020;11(1):2454.
65.
Bai K, Fang L, Liu T-R, Li J-Z, Wan D, Xiao M. Nonlinearity-enabled higher-order exceptional singularities with ultra-enhanced signal-to-noise ratio. Natl Sci Rev. 2023;10(7):nwac259.
66.
Lee CH, Imhof S, Berger C, Bayer F, Brehm J, Molenkamp LW, Kiessling T, Thomale R. Topolectrical circuits. Commun Phys. 2018;1(1):39.
67.
Imhof S, Berger C, Bayer F, Brehm J, Molenkamp LW, Kiessling T, Schindler F, Lee CH, Greiter M, Neupert T, et al.Topolectrical-circuit realization of topological corner modes. Nat Phys. 2018;14(9):925–929.
68.
Lee CH, Sutrisno A, Hofmann T, Helbig T, Liu Y, Ang YS, Ang LK, Zhang X, Greiter M, Thomale R. Imaging nodal knots in momentum space through topolectrical circuits. Nat Commun. 2020;11(1):4385.
69.
Pan N, Chen T, Sun H, Zhang X. Electric-circuit realization of fast quantum search. Research. 2021;2021:Article 9793071.
70.
Zhang J, Chen T, Deng W, Tong X, Zhang X. Quantum 2-player games and realizations with circuits. Research. 2024;7:Article 0480.
71.
Liu S, Gao W, Zhang Q, Ma S, Zhang L, Liu C, Xiang YJ, Cui TJ, Zhang S. Topologically protected edge state in two-dimensional Su–Schrieffer–Heeger circuit. Research. 2019;2019:Article 8609875.
72.
Luo K, Yu R, Weng H. Topological nodal states in circuit lattice. Research. 2018;2018:Article 6793752.
73.
Cao W, Wang C, Chen W, Hu S, Wang H, Yang L, Zhang X. Fully integrated parity-time-symmetric electronics. Nat Nanotechnol. 2022;17(3):262–268.
74.
Du J, Shi F, Kong X, Jelezko F, Wrachtrup J. Single-molecule scale magnetic resonance spectroscopy using quantum diamond sensors. Rev Mod Phys. 2024;96(2):Article 025001.
Year 2026 volume 9 Issue 2
PDF
102
55
Cite this Article
BibTeX
Article Info
doi: 10.34133/research.1091
  • Receive Date:2025-10-18
  • Online Date:2026-10-16
  • Published:2026-02-20
Article Data
Affiliations
History
  • Received:2025-10-18
  • Revised:2025-12-16
  • Accepted:2025-12-22
Affiliations
    Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements of Ministry of Education, School of Physics, Beijing Institute of Technology, Beijing 100081, China.

Corresponding:

* Address correspondence to: (T.C.);
* Address correspondence to: (X.Z.)
References
Share
https://castjournals.cast.org.cn/joweb/research/EN/10.34133/research.1091
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)
属
Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT