Experimental realization of non-Abelian permutations in a three-state non-Hermitian system

Abstract Eigenstates of a non-Hermitian system exist on complex Riemannian manifolds, with multiple sheets connecting at branch cuts and exceptional points (EPs). These eigenstates can evolve across different sheets—a process that naturally corresponds to state permutation. Here, we report the first experimental realization of non-Abelian permutations in a three-state non-Hermitian system. Our approach relies on the stroboscopic encircling of two different exceptional arcs (EAs), which are smooth trajectories of order-2 EPs appearing from the coalescence of two adjacent states. The non-Abelian characteristics are confirmed by encircling the EAs in opposite sequences. A total of five non-trivial permutations are experimentally realized, which together comprise a non-Abelian group. Our approach provides a reliable way of investigating non-Abelian state permutations and the related exotic winding effects in non-Hermitian systems.


INTRODUCTION
Permutation is a process of both fundamental and practical importance. For example, one way to distinguish fermions from bosons is to consider the exchange of the wave functions of two or more identical particles. Permutations of multiple states can emerge as the phenomenon of multi-state geometric phases [1,2]. They are generally non-commutative and can therefore be mapped to non-Abelian groups. This perspective suggests the possibility of emulating non-Abelian permutations by the parallel transport of three or more degenerate states. However, despite notable attempts in the fields of optics [3], cold atoms [4] and other topological systems [5,6], its realization remains a considerable experimental challenge, with the excitation and manipulation of multiple degenerate but coupled modes posing a major obstacle.
Recent advances in non-Hermitian physics have sparked the development of many intriguing applications related to optics and other classical waves [7,8]. Although non-Hermitian systems can be straightforwardly constructed from Hermitian systems by the inclusion of loss and/or gain or non-reciprocal hopping, they possess unique char-acteristics that are not found in their Hermitian counterparts. Perhaps the most notable distinction is that the eigenvalues are generally complex numbers. This simple fact permits the existence of multiple eigenvalue Riemann sheets connected at branch cuts [9][10][11]. The endpoints of the branch cuts are branch-point singularities known as exceptional points (EPs). Encircling an EP inevitably crosses one or multiple branch cuts-a process that causes the eigenstates to be exchanged and can even produce fractional winding numbers [11][12][13][14][15][16][17]. These fascinating behaviors, which are useful for topological energy transfer [11] and asymmetric mode switching [16] applications, have a topological origin: a non-Hermitian system lives on a complex Riemannian manifold that naturally permits state permutations. Hence, non-Hermitian systems offer a new vantage point for the study of state permutations. Recent theoretical investigations suggest that the encircling of multiple order-2 EPs or higherorder EPs is non-Abelian in character [13,[18][19][20][21] and can give rise to a myriad of exotic winding effects [16,[22][23][24][25]. These findings are consistent with the group theory point of view, as at least three degrees of freedom are required for non-Abelian processes to emerge. However, experimental confirmations of these proposals are lacking.
In this work, we theoretically investigate and experimentally realize the non-Abelian permutations of three states in a non-Hermitian system. By embedding the system's Riemannian manifolds in a 3D parameter space, two exceptional arcs (EAs), smooth trajectories of order-2 EPs, are found. As we will show, encircling them induces a unique permutation of the eigenstates. Five distinct types of state permutations are realized by encircling the EAs individually or sequentially. These five permutations, together with an identity element, holistically form a dihedral group of degree three, called the D 3 group, which can be used to describe the symmetry operations on an equilateral triangle, as shown in Fig. 1. All five permutations in Fig. 1a and the equivalent permutations in Fig. 1b are experimentally realized via a stroboscopic approach [17,23,[26][27][28] in acoustic experiments. We further show that the permutation operations are described by 3 × 3 unitary matrices, also known as U (3) non-Abelian Berry phases (NABPs), which connect the three-state evolutions on the system's complex Riemannian manifold.

EAs in a three-state non-Hermitian system
We begin with an exceptional nexus (EX) that emerges in a three-state non-Hermitian Hamiltonian H = (ω 0 + i γ 0 )I + H EP , where ω 0 + i γ 0 denotes the complex onsite energy and H EP determines the core physics and has the following form: (1) H EP lives on a 3D parameter space spanned by (η, ζ, ξ ) ∈ R 3 . There is also another parameter g , which for the convenience of discussion is not regarded as a separate dimension. Here, all coefficients are normalized by κ (where κ < 0), which is the hopping coefficient between neighboring sites. A ternary cavity system can be used to experimentally realize the Hamiltonian in acoustics, as shown in Fig. 2a. The second-order cavity mode is chosen as the onsite resonance mode. The parameters η and ξ represent detuning to onsite resonant frequencies, while i ζ and i g are introduced as losses. Figure 2b shows the three eigenmode profiles from a lower frequency (State 1) to a higher frequency (State 3) in the absence of non-Hermiticity. More details of the experimental set-up are given in Section III of the supplementary information. When g = 0, an EX exists at (η, ζ, ξ ) = (0, 0, 0) , which is an order-3 EP that connects to four EAs [23], each of which is a trajectory of order-2 EPs formed by two of the three eigenstates of Equation (1). These three eigenstates constitute a Hilbert space, which can be figuratively referred to as a fiber, at each parametric point (η, ξ, ζ ), thus forming fiber bundles that stick to the base manifold in the parameter space. The non-Hermiticity of the system means that the three eigenvalue Riemannian sheets connect at branch cuts, which naturally allows the exchange of states by encircling the EPs. Hence, each EA can be characterized by the two surrounding eigenstates in the permutation. As shown in Fig. 1, the permutations μ 1 and μ 3 constitute two generating operations of the D 3 group and the other elements of the D 3 group can be generated by ordered operations of μ 1 and μ 3 , i.e. The identity element is not of interest here, since it generates no changes. The issue of how to realize two EAs that possess the μ 1 and μ 3 types of permutation is therefore crucial for the demonstration of non-Abelian permutations.
In order to achieve this, we introduce the second term in Equation (1). When g = 0, the EX splits into two order-2 EPs in the ζ ξ plane at η = 0. In this way, the four EAs converging at the EX become a pair of smooth EAs. Figure 2c shows the two EAs, denoted α and β, for g = 0.61. We note that the way in which the EAs connect is dependent on the sign of g (see Sections I and II of the supplementary information for details) and we focus on the case with positive values of g in the main text. This configuration allows us to trace the evolution of states around the EAs, making it suitable for analysing the non-Abelian permutation of states that is the focus of this work.

Two generating permutations by encircling an EA
We first demonstrate two generating operations, μ 1 and μ 3 , that exchange two of the three states. To facilitate the discussion, we order the eigenstates based on the real parts of the eigenfrequencies at the starting point of the loop. We set η = 0.33, which is depicted as a light-green plane in Fig. 2c. The ζ ξ plane intersects with both EAs at two EPs, as shown by the red dots on the eigenvalue Riemann surface (real part) in Fig. 2d. The purple loop encircles EAα, which is formed by the coalescence of States 2 and 3 at η = 0.33. Hence, one complete cycle must  cross the branch cut once, resulting in the swapping of the two states, and consequently the operation μ 1 : 123 → 132 is realized. Likewise, it is straightforward to see that μ 3 , which encloses EA-β, exchanges States 1 and 2, i.e. μ 3 : 123 → 213. These permutations are experimentally observed via a stroboscopic approach. The parameters of the acoustic system are tuned to the specific values defined by the chosen loop. To achieve this, a Green's function method is used to determine the experimental parameters at each parametric point from the measured pressure response spectra (the details of this process are presented in Sections IV and IX of the supplementary information). The complex eigenfrequencies are then obtained by using the above parameters from the Green's function method and their real parts are plotted as the open circles in Fig. 3b and e for μ 1 and μ 3 , respectively. The solid lines in the figure show the theoretical results and their colors share the same notation as in Fig. 1. The measured eigenvalues are schematically labeled on the Riemann surfaces in Fig. 3a and d, which clearly delineate the evolutions associated with μ 1 and μ 3 . The salient feature that two states exchange at the branch cuts is clearly seen, and thus Fig. 3b and e agree well with our expectation.
Next, the evolutions of eigenfunctions are also obtained experimentally by measuring the acoustic field profile in all three cavities (see Section V of the supplementary information for details). The results indeed show the swapping of eigenfunctions across the branch cut where the real parts of the eigenvalues cross. In Fig. 3c and f, we plot the representative eigenfunctions at five chosen points along the encircling path and state exchanges are observed. The results shown in Fig. 3c can be taken as an example. We see that at starting Point I, State 2 (shown in green) has a large amplitude at Site A (the middle site). As the system is driven along the μ 1 path, the amplitude at Site A gradually diminishes, while that at Site C increases. At the last two points, State 2 at Point IV smoothly connects to State 3 at Point V, as a direct consequence of crossing the branch cut. Likewise, State 3 at Point IV (shown by the red lines) connects to State 2 at Point V. Meanwhile, State 1 remains almost unchanged throughout the evolution. Upon the completion of one closed cycle, the final outcome is the exchange of States 2 and 3.
We further remark that as the parameters change, the eigenfunctions of the three states also vary. It is therefore crucial to correctly identify how the eigenfunctions evolve along the parametric points, especially in the vicinity of the branch cut where the state exchange takes place. We examine the inner products for all the neighboring states, i.e. | ψ L i,l +1 |ψ R j,l | 2 , where |ψ R j,l is the right eigenfunction of the j th state at the parametric point l and ψ L i,l +1 | is the left eigenfunction of the i th state at point l + 1, where i, j = 1, 2, 3. The two neighboring eigenfunctions that yield an inner product close to unity are connected by parallel transport [29]. This procedure was performed for all states at all the parametric points presented in our work.
The state permutation induced by μ 1 can be captured by a U (3) NABP [1] (see Section VI of the supplementary information for details). Using the eigenvectors of H EP as a basis, the NABP for μ 1 is: From Equation (2), we can further obtain a multiband Berry phase as: This phase factor can be observed as a π -phase difference between State 2 at Points I (shown in green) and V (shown in red) in Fig. 3c. These results are consistent with the knowledge that an order-2 EP possesses a fractional winding number of 1/2 and the fact that encircling the EP twice restores both states with a Berry phase of π . The μ 3 -induced state permutation can also be seen by tracing the eigenfunction evolutions in Fig. 3f and its corresponding NABP is: which also yields a Berry phase of μ 3 = −π from Equation (3). Although μ 1 and μ 3 are the same, the two NABPs U μ 1 and U μ 3 are different and they do not commute.
The three-state permutations are achieved by sequentially encircling both EA-α and EA-β. Without loss of generality, we can anchor the two loops μ 1 and μ 3 at a common vantage point P η, ζ, ξ = (0.33, 0, 0), as depicted by the black hexagon in Fig. 2c. The point P is also the starting and end point of the encircling. In Fig. 4a-c, the μ 3 operation is executed first by encircling EA-β, which swaps States 1 and 2. The μ 1 operation is then carried out by encircling EA-α, thus exchanging the new States 2 and 3. The net result is the swapping of all three states, as defined by ρ 1 . The ρ 2 operation is also experimentally achieved by first encircling α and then β, as shown in Fig. 4d-f. The two experimental outcomes, i.e. the mode profiles at the parametric Point VII in Fig. 4c and f, are clearly distinct, thus unambiguously validating the non-Abelian characteristics. Again, we can summarize the three-state permutations with the NABPs as: It is clear that U ρ 1 = U ρ 2 , although the Berry phases in both cases are ρ 1 = ρ 2 = 0 [mod(2π )], as can be identified from the mode profiles at the parametric Points I and VII in Fig. 4c and f. This verifies the non-Abelian character by encircling different types of EAs.

Multiple permutations of encircling EAs
We have already demonstrated four of the five non-trivial permutations depicted in Fig. 1a. The remaining operation is μ 2 : 123 → 321, which exchanges States 1 and 3. Following the rule used to label the EAs, we would expect that EPs exist near to the EX that correspond to the permutation μ 2 . This can be attained by shifting the encircling loop of EA-α to η = 0, as shown in Fig. 2c. At first sight, it seems counterintuitive that μ 2 can exist in our system, since the hopping between Sites B and C is zero in Equation (1). To see how μ 2 emerges, we first note that η represents onsite detuning in Sites B    and C, thus letting η cross zero causes the inversion of the lowest and highest frequency modes (States 1 and 3). At η = 0, the two order-2 EPs (EP-α and EP-β in Fig. 2e) are linked by a branch cut that is parallel to the ζ axis, which connects the lowest and highest frequency sheets. Hence, an evolution that follows the blue loop in Fig. 2e exchanges States 1 and 3 and leaves State 2 unchanged, thus realizing μ 2 .
The μ 2 operation is also experimentally realized using our acoustic system. The results for the eigenvalues and eigenfunctions are shown in Fig. 3h and i, respectively, where the exchange of States 1 and 3 can clearly be seen. We have also computed the corresponding NABP: and μ 2 = −π . We further remark that, as an element in D 3 , . This indicates that the permutation μ 2 can be treated as the operation of encircling the EAs multiple times in our non-Hermitian system. To show this, we can shift the position of the blue loop in Fig. 2c slightly to η = 0.055, so that it transverses three different branch cuts, with each traversal exchanging two states. These results are presented in Section VII of the supplementary information. Since μ 2 completes the D 3 group here, all other operations that encircle the EAs in Fig. 2c multiple times must be equivalent to the single operation shown in Figs 3 and 4.

DISCUSSION AND CONCLUSION
A common practice for characterizing topological manifolds is to consider equivalence classes of loops, in which winding numbers play a vital role. Non-Hermitian topology can be characterized by the eigenvalue winding number, sometimes called the eigenvalue vorticity or discriminant number [28,30,31], which is often considered to be sufficient to reveal the topological structure of the complex Riemann surfaces. However, our results show that eigenvalues are not directly associated with state permutations. Even the eigenvector winding numbers underlain by the Berry phase do not contain explicit information on state permutations. The state permutations and their non-Abelian characteristics are disclosed either by tracing the parallel transport of all three states or by computing the NABP matrix. Hence, the EAs and their interactions constitute the non-Hermitian counterparts of the knot and link structures of nodal lines in Hermitian band structures [19,20,32]. A question naturally arises as to how the winding numbers relate to the non-Abelian permutations demonstrated in this work. To illustrate this, we recall that the two processes defined by ρ 1 and ρ 2 yield identical Berry phases ρ 1 = ρ 2 = 0 [mod(2π )], which can be regarded as the same eigenvector winding number. We have numerically confirmed that the eigenvalue winding numbers for ρ 1 and ρ 2 are also identical in these two cases and are consistent with their Berry phases. As discussed above, the evolutions ρ 1 and ρ 2 are equivalent to performing both μ 1 and μ 3 in opposite orders. However, the two concatenated loops μ 1 and μ 3 are equivalent to the larger loop encircling both EA-α and EA-β (see Section VIII of the supplementary information). When this loop is followed, three complete cycles are needed to restore all three states, which gives rise to a fractional winding number of 2/3 [23,33]. In other words, one complete parametric cycle following ρ 1 and ρ 2 does not recover all the states. It follows that the states after one cycle are dependent on the states at the starting point. This is the reason for the non-Abelian outcomes demonstrated in our work.
In summary, we have successfully demonstrated that all the non-trivial operations comprising the D 3 group can be realized by encircling EAs in a three-state non-Hermitian system. Our work builds on recent developments in non-Hermitian physics that have introduced a kaleidoscope of EP structures with distinct topological characteristics. Experimentally, our studies are based on the stroboscopic approach so that the non-adiabatic transitions typically encountered in dynamic evolutions can be avoided [16,[34][35][36]. Our work and the methodology can be extended to study knot and link structures formed by different EAs [37][38][39][40]. The combined strength of these theoretical developments and experimental techniques in non-Hermitian physics, in conjunction with the rich arsenal of non-Abelian theories, will open new avenues to the discovery of exotic phenomena and the development of rich applications in a diversity of fields. For example, non-Abelian permutations around multiple EAs provide additional degrees of freedom to manipulate wave propagation [16] and on-chip energy transfer [11]. Relating to our work are several recent studies proposing a new class of anyonic-parity-time symmetric systems [41,42] that can benefit applications such as lasers [43]. On the other hand, the existence and evolutions of multiple EPs in a multi-parameter phase space give rise to rich opportunities of more sophisticated usage of EPs, which may benefit applications such as sensors [44,45], absorbers [46,47], scattering control [48], etc.

DATA AVAILABILITY
Data are available upon request to the corresponding authors.

SUPPLEMENTARY DATA
Supplementary data are available at NSR online.