Arbitrary order exceptional point induced by photonic spin–orbit interaction in coupled resonators

Many novel properties of non-Hermitian systems are found at or near the exceptional points—branch points of complex energy surfaces at which eigenvalues and eigenvectors coalesce. In particular, higher-order exceptional points can result in optical structures that are ultrasensitive to external perturbations. Here we show that an arbitrary order exceptional point can be achieved in a simple system consisting of identical resonators placed near a waveguide. Unidirectional coupling between any two chiral dipolar states of the resonators mediated by the waveguide mode leads to the exceptional point, which is protected by the transverse spin–momentum locking of the guided wave and is independent of the positions of the resonators. Various analytic response functions of the resonators at the exceptional points are experimentally manifested in the microwave regime. The enhancement of sensitivity to external perturbations near the exceptional point is also numerically and analytically demonstrated.

I f the Hamiltonian of a system takes the form of a Hermitian matrix, it has only real eigenvalues accompanied by a set of orthogonal eigenvectors. If the Hamiltonian is non-Hermitian, its eigenvalues are complex in general. By tuning the amount of non-Hermiticity, one can make two eigenvalues coalesce at one point, which is called an exceptional point (EP) [1][2][3] . In contrast to the degeneracy in a Hermitian system (a diabolic point), an EP is a branch point in the complex energy surface where the Hamiltonian matrix becomes defective. EPs were found in non-Hermitian systems possessing parity-time (PT) symmetry, where the simultaneous increase of gain and loss can induce symmetry breaking, characterised by a change of frequency spectrum from being purely real to being complex 4,5 . EPs can also exist in passive systems without PT symmetry as a result of the unbalanced competition between loss contrast and coupling strength of two states 6 . The phase transition associated with EPs leads to numerous counterintuitive phenomena such as loss (pump)-induced revival (depressing) of lasing 6,7 , optical isolation 8 , unidirectional transport of light 4,[9][10][11] , topological energy transfer 12 and others [13][14][15][16][17][18][19][20][21] . Recently, higher-order EPs have also attracted a lot of attention [22][23][24][25][26][27] . A higher-order EP is a result of coalescence of multiple eigenvalues 22 , which normally requires a delicate variation of parameters in a larger parameter space and is much more difficult to achieve 25,28 . For a simple two-level non-Hermitian system with a Hamil- , an EP can be obtained under the condition of (a − b) 2 + 4cd = 0. This can be achieved when Re(a) = Re(b), a − b is purely imaginary and c ¼ d Ã 4,6-9 or when a − b is real and c = d is purely imaginary [28][29][30][31] . In the former case, the non-Hermiticity comes from the asymmetric loss/gain of the two identical states, whereas in the latter case the non-Hermiticity appears in the coupling of two states mediated by some open channels. However, EPs can also be achieved when a = b and either c = 0 or d = 0, corresponding to a unidirectional coupling between two states. EPs of this kind can be achieved in systems with structural asymmetry 32,33 .
In this paper, we propose a different scenario to achieve unidirectional couplings. Our systems consist of a set of identical dipolar resonators sitting on a strip waveguide. The unidirectional coupling between any two resonators arises from the evanescent couplings between the chiral dipolar modes of the resonators and the waveguide mode. The transverse spin-momentum locking of the guided wave excludes the backward couplings among the resonators. EPs achieved in this manner can be of arbitrary order, depending only on the number of resonators. Although the positions of the resonators do not affect the order of an EP, they do determine the spectral response functions of the resonators. We have studied such systems theoretically using both coupled mode theory (CMT) 34,35 and full-wave simulations. An experimental demonstration of the theory is also carried out in a system consisting of multiple high-dielectric spheres on a strip waveguide in the microwave regime.

Results
Exceptional points in a model system. We first consider a model system consisting of two identical spheres located near a strip waveguide that forms a circle, as shown in Fig. 1a. The purpose of studying this model system is twofold. It not only proves that the model system in the next subsection, where EP appears naturally, is just a limiting case of the model system in Fig. 1a, but also uncovers the critical role of unidirectional coupling in achieving arbitrary order EPs. In cylindrical coordinates (r,θ,z), the waveguide's geometry is defined by R ≤ r ≤ R + t and −w/2 ≤ z ≤ w/2, with t being the waveguide's thickness and w being the waveguide's width. The waveguide has a rectangular cross-section and is made of gold. The positions of the two spheres are denoted by two azimuthal angles θ 1 = 3π/2 and θ 2 = π/2. The permittivities of the spheres are described by the Drude model, which gives an electric dipole resonance at about 300 THz.
The strip waveguide shown in Fig. 1a supports guided plasmon modes. The dispersion relation of the lowest-order mode is shown in Fig. 1b by the solid blue line, where n eff = k plas /k 0 is the effective mode index, with k plas being the mode propagation constant and k 0 being the vacuum wave number. The corresponding mode electric field is shown in the inset of Fig. 1b, where the colour denotes the field magnitude and the arrows indicate field directions. We note that the electric field is well confined around the waveguide, facilitating long-range coupling between any two dipole resonators through the guided plasmon 36 . Because the plasmon field (i.e. evanescent electromagnetic field) in Fig. 1a has both transverse (E r ) and longitudinal (E θ ) electric field components, it carries intrinsic spin angular momentum (SAM) in the z direction, i.e. the plasmon field is in general elliptically polarised 37 . The SAM is determined by the material and geometrical properties of the waveguide. It has been shown that such a propagating plasmon can be mapped into a quantum spin Hall system with the property of spin-momentum locking (plasmon carrying opposite SAM propagating in opposite directions) 38,39 .
Each sphere in Fig. 1a supports two degenerate electric dipole modes alongx andŷ directions, respectively. A superposition of the two modes gives rise to a pair of chiral dipole modes p ¼ pê ± , whereê ± ¼x ± iŷ ð Þ= ffiffi ffi 2 p is the circular basis of the spheres [40][41][42][43] . Such chiral modes should not be confused with the eigenmodes of a sphere made of chiral materials, which has intrinsic chirality 44 . Below, we will focus on the interaction between such chiral dipole modes mediated by the waveguide mode. In the presence of the strip waveguide, the resonant chiral dipole field asymmetrically couples to the propagating plasmon of the waveguide 45 . We carefully designed the waveguide (tuning material parameters and cross-sectional geometry) so that the SAM of the chiral dipole field matches the SAM of the guided plasmon. Because of the locking between the SAM and the linear orbital momentum (k), the excited plasmon propagates in one direction only. Such a process can be regarded as the conversion of light's SAM to linear orbital momentum (i.e. photonic spin-orbit interaction 45 ), which can produce many interesting applications/phenomena such as chiral photonic waveguides 46 and lateral optical forces 44,47,48 . The resulting coupling between the two spheres is unidirectional, i.e. either clockwise or counter-clockwise, as shown in Fig. 1a, depending on the SAM of the incident light. The system can be described by an effective Hamiltonian Here, ω 0 is the dipole resonance frequency of the isolated spheres, γ 1,2 and γ c denote the material loss and radiation loss, respectively, of the spheres and κ 21 (κ 12 ) denotes the coupling parameter from Sphere 1 (2) to Sphere 2 (1), as shown in Fig. 1a for the case under left-handed circular polarisation (LCP) excitation. Here, the coupling parameters have complex values in general, which are in contrast to those systems where the mode couplings are realised through near-field hopping 6,49 and are usually real-valued. The corresponding eigenfrequencies of Eq. are , where Γ = (γ 1 + γ 2 + 2γ c )/2 and Δγ = γ 2 − γ 1 denotes the loss contrast of the two spheres. Using CMT, the amplitudes of the electric chiral dipole ARTICLE NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-019-08826-6 moments can be determined as (see Methods) where A in is the incident field amplitude.
In general, the Hamiltonian in Eq. (1) is non-Hermitian even if the system has no loss. The non-Hermiticity here originates from the unidirectional coupling which leads to κ 12 ≠ (κ 21 ) Ã . We note that this does not imply a broken reciprocity since the Hamiltonian only describes a subsystem corresponding to the coupling of one chiral mode. The other chiral mode (with opposite chirality) couples via the plasmon propagating in the opposite direction. The two scenarios have equal coupling strength and therefore the system is reciprocal 40,50 . The value of ω ± can be tuned by varying the product κ 12 κ 21 , which can be expressed as κ 12 κ 21 ¼ Àiκ 0 12 e iϕ 12 À Á Àiκ 0 21 e iϕ 21 À Á ¼ Àκ 0 12 κ 0 21 e i ϕ 12 þϕ 21 ð Þ51 . Here, κ 0 12 and κ 0 21 are the positive real-valued coupling strengths and ϕ 12 + ϕ 21 accounts for the total dynamic phase accumulated by the guided plasmon. An EP can be achieved when κ 12 κ 21 = (Δγ/4) 2 .
Since the loss contrast Δγ has a real value, this requires κ 12 κ 21 to be real and positive, therefore, ϕ 12 + ϕ 21 = (2m + 1)π, m = 0,1,..., which corresponds to destructive interference of the plasmon. Such a condition is easily fulfilled by tuning the radius R of the waveguide. Figure 1c shows the plasmon amplitude when the circular waveguide in Fig. 1a is excited by a point source. The peaks correspond to constructive interferences similar to the Mie resonances associated with dielectric spheres/disks. The dips marked by the red arrows correspond to destructive interferences, where κ 12 κ 21 ¼ κ 0 12 κ 0 21 . The value of the product κ 0 12 κ 0 21 can be controlled by tuning the loss of the waveguide. Therefore, it is possible to achieve an EP when κ 12 To demonstrate such a possibility, we did full-wave simulations of the system shown in Fig. 1a with R = 200 nm under the incidence of a plane wave of LCP, in which case the left-handed chiral dipole modes were excited and coupled with each other. The electric dipole moments of the spheres were calculated and fitted to the CMT expressions in Eq. (2), from which we obtained the values of γ 1,2 , γ c and κ 12 , κ 21 . In the simulations, we only introduced loss into the lower-half waveguide (θ ∈ [π/2,3π/2]) shown in Fig. 1a to suppress κ 12 κ 21 (see Methods for details). When the introduced waveguide loss γ wg was increased, an EP emerged at κ 12 κ 21 = (Δγ/ 4) 2 , where the real parts of the eigenfrequencies merge and their imaginary parts bifurcate. Note that γ wg contributes little to the phase of κ 12 κ 21 (see Supplementary Figure 1). Figure 1d shows such a scenario for different values of Δγ. As Δγ is reduced, the EP moves to the right with an enlarged splitting of Re (ω ± ) and reduced splitting of Im(ω ± ). In the limit of Δγ → 0, i.e. approaching the limit of identical spheres, the EP will be reached when κ 12 → 0 (γ wg is sufficiently large so that the coupling from p 2 to p 1 is completely suppressed). In this case, the system Hamiltonian is reduced to a Jordan block form (2) reduces to Different from Eq. (2), Eq. (3) exhibits a single pole but of different orders for p 1 and p 2 . The presence of a second-order pole in p 2 is a signature of EPs 51,52 . It is noted that the unidirectional coupling renders p 1 unperturbed, whereas it contributes to p 2 . The interference between the first-order and second-order terms in p 2 can give rise to an asymmetric spectrum different from the Lorentzian response in p 1 53 . Below, we will focus on EPs in the case of unidirectional coupling.
Exceptional point in coupled-dipole-resonator model system with mutual coupling. a Two spherical dipole resonators (p 1 , p 2 ) placed near a circular plasmonic waveguide. Their chiral dipole resonance modes couple with each other through the unidirectionally propagating plasmon under spin-momentum locking. b Effective mode index of the guided plasmon. The inset shows the magnitude (colour) and direction (arrows) of the mode electric field. c Amplitude of the guided plasmon as a function of radius R under the excitation of a point source near the circular waveguide. The peaks correspond to constructive interferences of the plasmon and the dips (marked by the red arrows) correspond to destructive interferences. d Real and imaginary parts of the eigenfrequencies as a function of the loss γ wg of the lower-half waveguide (θ ∈ [π/2,3π/2]) obtained by coupled mode theory. Δγ is the loss difference of the two spheres Exceptional points achieved with unidirectional coupling. At the EP of the case where Δγ → 0, the coupling from p 2 to p 1 is suppressed by the lossy waveguide, and thus the model system shown in Fig. 1a can be simplified to that shown in Fig. 2a, where the two spheres are placed on a straight waveguide and are separated by a distance d. The guided plasmon propagates unidirectionally under spin-momentum locking and is absorbed at the end of the waveguide under an absorbing boundary condition. The absorption of the plasmon prevents reflection at the end of the waveguide, and the coupling between spheres is unidirectional (i.e. κ 12 = 0). We carried out full-wave simulations for such a configuration. The simulation details and system parameters are given in the Methods section. The electric dipole moments of the two spheres are shown by the blue (p 1 ) and red (p 2 ) symbols in Fig. 2b-k for d = 2.0λ plas (Fig. 2b, c), d = 2.2λ plas (Fig. 2d, e), d = 2.4λ plas (Fig. 2f, g), d = 2.6λ plas (Fig. 2h, i), and d = 2.8λ plas (Fig. 2j, k). Here, λ plas = 2π/k plas is the wavelength of the guided plasmon. For d = 3.0λ plas , the spectra of two dipole moments are identical to those of d = 2.0λ plas , exhibiting a periodic pattern with λ plas being the period. We note that p 1 has a typical Lorentzian shape and remains unchanged when d changes, whereas p 2 changes dramatically because of the interference effect under unidirectional coupling. The CMT results are also shown in Fig. 2b-k as solid lines, which match well with the simulation results. In Fig. 2c, p 2 shows a dip at the resonance frequency, which can be understood from Eq. (3). At ω = ω 0 , Eq. (3) gives  Fig. 2l. When d is an integer multiple of λ plas , κ 21 has no phase contribution from the guided plasmon and we have κ 21 ¼ À iκ 0 21 51 , and therefore Z(d) has a minimum value of Γ=2 À κ 0 21 , which explains the dip of p 2 in Fig. 2c. When d is a half-integer multiple of λ plas , Z(d) has a maximum value of Γ=2 þ κ 0 21 , which explains the enhancement of p 2 in Fig. 2g. Note that we have ignored the absorption of the waveguide in Fig. 2a for simplicity, but the results would be qualitatively the same if absorption is also considered.
Higher-order exceptional points. The realisation of higher-order EPs is hotly sought after, not only for the interesting physics, which goes beyond the square-root singularity, but also for their application potential such as for ultrasensitive sensors 32,54-57 . In conventional non-Hermitian systems, the realisation of higherorder EPs is usually very difficult as it involves complicated design and fine-tuning of many parameters 22,25,26 . We show that EPs of arbitrary order can be readily achieved in our system by introducing more identical resonators on the waveguide. For N identical spheres on the waveguide with unidirectional couplings, as shown in Fig. 3a, all of the coupling parameters κ ij , i < j vanish so that the effective Hamiltonian can be expressed as The Hamiltonian matrix in Eq. (4) represents an EP of order N. It is defective, having a single eigenvalueω 0 ¼ ω 0 À iΓ=2 with algebraic multiplicity N, i.e. the N eigenvalues coalesce. Thus, the order of the EP can be increased simply by putting more resonators on the waveguide. The coupling constant κ ij will only affect the forms of the spectral response functions of the resonators. Here, we consider an example with five equally spaced spheres with a separation d = λ plas , as shown in Fig. 3a. Figure 3b shows the electric dipole moments of the five spheres under LCP excitation. The solid lines denote the results obtained from CMT with the Hamiltonian in Eq. (4). The symbols denote the results of full-wave simulations. Excellent agreement between the two is found. It is noted that only p 1 shows a Lorentzian shape, whereas the rest are suppressed because of interference. Therefore, at the resonance frequency, the dipole moment as a function of the resonator's position has a monotonic trend, as shown in Fig. 3c, where each point denotes the dipole moment magnitude of the corresponding sphere. Under LCP excitation, the dipole moment gradually decreases from left to right, whereas this trend is reversed under right-handed circular polarisation (RCP). Such a phenomenon can be understood based on CMT. Because the period of the array is d = λ plas , all of the coupling parameters in Eq. (4) have zero phase, i.e. κ ij ¼ κ 21 ¼ Àiκ 0 21 ; 8i>j. Thus, at ω = ω 0 , one can easily derive the expressions for the resonant dipole moments using CMT (Supplementary Note 1) as . Because both κ 0 21 and Γ take positive real values and Γ À 2κ 0 21 À Á =Γ <1 according to the results in Fig. 2l, we have p i < p i−1 , which indicates a monotonic decrease of the dipole moments and the decrease is reversed for the RCP excitation by mirror symmetry. Such a decrease can also be understood as a kind of cascade effect under destructive interferences between the spheres, where the coupling field from p i suppresses the field of the spheres on its right side. It is also possible to have constructive interference for the spheres, such that the sphere at the end has much larger emission power than an isolated sphere (see Supplementary Figure 2). The mechanism revealed here enables spin-controlled dipole array emission and energy accumulation in the lateral direction.
Higher-order EPs are characterised by an enhanced sensitivity to external perturbations. In the case of small perturbations, the eigenfrequency splitting induced near an EP of order N is of the order ffiffi ffi ϵ N p , where ϵ denotes the strength of the small perturbation. Because the order of EP in our system is directly determined by the number of spheres N, the corresponding sensitivity can be easily controlled. To simplify the numerical calculation, we demonstrate the higher-order sensitivity enhancement in our system by considering a two-dimensional (2D) version of our model system, i.e. five parallel cylinders on a 2D metal waveguide, as shown in Fig. 4a. The cylinders are equally spaced with a distance of d = 5λ plas . The waveguide is truncated on the right side with a distance of 5λ plas between its edge and the last cylinder. An absorption boundary condition is applied on the left edge. Because the waveguide truncation on the right edge acts as a mirror, a propagating LCP plasmon from left to right will be totally reflected by the truncated edge. Because the reflected plasmon is of RCP, it will also excite the RCP modes of the five cylinders. Thus, the system involves 10 unidirectionally coupled modes with five LCP modes travelling to the right and five RCP modes travelling to the left and is at an EP of order 10.
To demonstrate the enhanced sensitivity of order ffiffi ffi ϵ 10 p , perturbation is introduced by inserting a small parallel cylinder with an elliptical cross-section between the first cylinder and the left edge. The strength of perturbation ϵ is proportional to the reflection caused by the small elliptical cylinder, which can be controlled by the amount of loss in the small cylinder (see Methods). Figure 4b, with results calculated numerically using COMSOL, shows the spectra of total electric dipole moment strength of the system at different values of ϵ. We notice that as ϵ increases, more resonance peaks appear in the spectrum, and the splitting between peaks increases as well. When ϵ is large enough, a total of 10 peaks are conspicuous, corresponding to 10 eigenvalues, as expected. In Fig. 4c, we plot the splitting of the eigenvalues Δω n = |ω n − ω 0 | as a function of ϵ for the first and the last peaks shown in Fig. 4b in log-log scales. For small ϵ values, the results show a straight line with a slope close to 1/10, i.e. ffiffi ffi ϵ 10 p behaviour for an EP of order 10.
The above system can be described by the following Hamiltonian: where the two N × N block matrices H L and H R denote the unidirectionally coupled N LCP modes and N RCP modes, respectively. They both have the form shown in Eq. (4), with N = 5 and κ ij ¼ Àiκ 0 21 . The block matrix κ RL denotes the one-way conversion from LCP modes to RCP modes occurring at the right edge of the waveguide and κ LR denotes the excitation of rightgoing LCP modes at the perturbation cylinder on the left side of the waveguide caused by reflection. These two block matrices have the forms κ RL ij ¼ Àiκ 0 21 and κ LR ij ¼ Àiϵκ 0 21 , respectively. Note that the corresponding eigenvectors have the form of  Figure 4d shows the real parts of the normalised eigenvalues ω n =κ 0 21 as a function of perturbation strength ϵ, which exhibits the eigenvalue splitting phenomenon shown in Fig. 4b. Note that ω 5 and ω 6 have identical real parts. In the numerical simulations, all 10 eigenvalues can be observed because the spacing between cylinders is not exactly an integer multiple of plasmon wavelength. In Fig. 4e, we plot the splitting of Δω n =κ 0 21 as a function of ϵ in log-log scales, where only n = 1,2,3,4 is considered because of the symmetric splitting of the eigenvalues. Enhanced sensitivity of order ffiffi ffi ϵ 10 p can be clearly seen in the small ϵ region, where all of the splitting follows a trend with a slope of 1/10. When the perturbation is large, only Δω 1 =κ 0 21 maintains a 1/10 slope, as shown by the zoom-in inset in Fig. 4e, in accordance with the numerical results in Fig. 4c. The trends of the other eigenvalue splitting deviate from the 1/10 slope because of the contributions of high-order perturbation terms.
Experimental demonstration at microwave frequencies. The physics discussed above also apply to the microwave regime. We note that the magnetic field of light can also carry SAM. Consider a dielectric sphere which supports two orthogonal magnetic dipole modes alongx andŷ directions, respectively. A superposition of the two modes gives rise to a pair of chiral magnetic dipole modes m ¼ mê ± , one of which is shown by the inset in the lower right corner of Fig. 5a. A simple dielectric waveguide can support guided wave with evanescent magnetic field carrying spin. The combination of the dielectric spheres and dielectric waveguide then corresponds to the magnetic version of the plasmonic model system and can be easily realised at microwave frequencies. We carried out an experimental demonstration of the properties associated with the EPs in our system using the setup shown in Fig. 5a. We used high-dielectric spheres as the resonators, which have a magnetic dipole resonance near 9 GHz. The waveguide was also made of a dielectric material, and its fundamental mode served to guide the wave in the considered frequency region. The two ends of the waveguide were inserted into absorbing materials to suppress the reflected waves. The chiral magnetic dipole mode was excited using a circularly polarised microwave incident along the z direction. To measure the amplitude of the magnetic dipole moments, a probe was used to detect the magnetic field above the spheres (see Fig. 5a). Figure 5b shows the magnitude of the magnetic field near the left and right spheres, respectively, under LCP excitation. The results of RCP excitation have similar features and are shown in Supplementary Figure 3. The symbols here denote the experimental results and the solid lines denote the analytical fitting results based on Eq. (3). The measurements were done at different coupling distances d. We note that as the coupling distance changes from d = 50 mm to d = 60 mm, the left dipole resonance remains unchanged, whereas the other dipole undergoes dramatic variations, in accordance with the numerical results shown in Fig. 2b-k. The position of the peak for the left sphere (marked by a dashed line) is a constant. The interference feature of the right sphere gradually changes as d increases. This can be seen from the positions of the maximum and minimum values near the resonance frequency ω 0 , shown in Fig. 5c, which undergoes a red-shift when d increases. Note that there is a switch of the maximum at d = 2λ wg = 55 mm, with λ wg being the wavelength of the guide wave. The experimental results match the analytical CMT results very well, as shown in Fig. 5c. The asymmetry of the spectra, as well as the evolution of the maximum/ minimum for the right sphere, can be explained by Eq. (2) with κ 12 = 0. The numerator of p 2 has a factor of Γ/2 − iκ 21 − iΔω, where Δω = ω − ω 0 . Because κ 21 is a complex number, it is obvious that p 2 is asymmetric with respect to Δω = 0 except for the special case of Re(κ 21 ) = 0. This can be seen from Fig. 5d, where the circular symbols denote the values of function Z(d) = Γ/2 − iκ 21 on the complex plane, obtained from the fitting results in Fig. 5b. For any case with Re κ 21 ð Þ≠0, the values of |Γ/2 − iκ 21 + iΔω| and |Γ/2 − iκ 21 − iΔω| (denoted by the length of the red arrows) are not equal, which indicates that p 2 has a shape that is asymmetric with respect to centre frequency ω 0 . It is the looping of the complex coupling parameter κ 21 that accounts for the shift of the minimum/ maximum in Fig. 5b, c. We note that similar phenomena also exist in the plasmonic model system shown in Fig. 2a except that the spectra asymmetry of p 2 is weak due to a smaller coupling strength κ 0 21 compared with Γ/2 (corresponding to a smaller circle in Fig. 2l). The experiment confirmed the mechanism of unidirectional coupling between identical resonators which can be described by the Hamiltonian in Eq. (4) (where an EP naturally arises). It also demonstrated the interference effect between the firstorder and second-order poles in Eq. (3). We emphasize that the second-order pole emerges in the Green's function as a direct result of the coalescence of two eigenvalues at the EP 52,58 . The interference effect at the EP is different from the ordinary interference of two resonators without an EP. In the latter case, the spectra of both resonators should be asymmetric in general and depend on the coupling. Here, as shown in Fig. 5b, the spectrum of Sphere 1 is approximately symmetric and does not depend on the value of coupling distance d, while the spectrum of Sphere 2 is asymmetric and depends on the value of d.
We also performed measurements for the higher-order EP, which manifested spin-controlled dipole array emission. Because of the finite length of the waveguide, we put four spheres on the waveguide with equal spacings of λ wg and measured the resonant magnetic field near each sphere under LCP and RCP excitations. The results are shown by solid circles in Fig. 6. In accordance with the prediction in Fig. 3c, the magnetic field has a monotonic trend determined by the SAM of the incident microwave wave. The slight asymmetry between LCP and RCP results is owing to the imperfection of the microwave source wavefront. We note that the field dependence in Figs interference effect in presence of higher-order poles at the EPs. We also note that, different from the dipole moment, the evanescent waves coupled from the spheres to the waveguide interference constructively since d = λ plas , which leads to a growing field magnitude of the waveguide (Supplementary Figure 4).

Discussion
For a physical system described by a Hamiltonian involving many physical parameters, it typically requires tedious tuning of multiple parameters to achieve a higher-order EP. Worse still, in a real system, the parameters may be correlated, making their tuning even more complicated. This issue renders a higher-order EP rather difficult to achieve. One feasible way is to have a physical mechanism that guarantees all of the parameters fulfil the condition of the EPs. In our system, the mechanism is the unidirectional coupling of identical resonators induced by the locking between the light SAM and linear orbital momentum. With such a mechanism, the phenomenon associated with arbitrary order EPs can be realised simply by expanding the system to include more spheres. The higher-order EPs demonstrated here should have potential applications in sensing. In such applications, resonators with high-quality factors are needed to detect small perturbations, although it is not necessary for a proof-ofconcept demonstration of the EPs here. In the model system with unidirectional coupling between identical spheres, we have assumed a lossless plasmonic waveguide. When loss is taken into account, the coupling strength between spheres will be reduced, and the coupling parameters in the Hamiltonian will have different magnitudes in general. However, the Hamiltonian will remain in the same form (i.e. a triangular matrix), and the discussed physics about EPs still hold. In this case, one has to assemble the spheres closely to obtain strong EP phenomena.
In summary, we have demonstrated that the spin-orbit interaction of light can induce unidirectional coupling between chiral dipole modes of identical spherical resonators through a waveguide if the chiral modes are selectively excited using circularly polarised light. The unidirectional coupling leads to arbitrary order EPs depending on the number of resonators on the waveguide. Various analytic response functions of the resonators at the EPs have been experimentally manifested. The discussed mechanism also provides an easy way to control the excitation of the coupled resonance modes using the SAM of light. The proposed configuration can serve as a convenient platform for studying light-matter interactions at EPs and can be integrated into a photonic circuit for on-chip applications.

Methods
Coupled mode theory for two dipole resonators. The rate equations for the two chiral dipole modes can be written as 34,35 where a i is the mode amplitude. The above equations can be re-written as where Under the excitation of a harmonic wave with a in = A in e −iωt , the steady-state solutions for a i = A i e −iωt obtained from Eqs. (6) and (7) are where Γ = (2γ c + γ 1 + γ 2 )/2. The magnitude of the induced dipole moments for each sphere are p 1 = |A 1 | and p 2 = |A 2 |. In the case of unidirectional coupling, e.g. κ 12 = 0 and κ 21 ≠ 0, we obtain Eq. (3).
Numerical simulation. For the plasmonic model systems shown in Figs. 1a and 2a, the spheres have diameters of 10 nm. The distance between the spheres and the waveguide is h p = 17 nm. The permittivities of the spheres are described by the Drude model ε i ¼ 1 À ω 2 p = ω 2 þ iωγ i À Á , where ω p = 3.3 × 10 15 rad s −1 and the sub-script 'i' denotes the ith sphere. For the system in Fig. 1a, we set γ 1 = 1.0 × 10 13 rad s −1 and γ 2 = γ 1 + Δγ with Δγ being the loss difference in Fig. 1d. For the system in Fig. 2a, we set γ 1 = γ 2 = 1.0 × 10 13 rad s −1 . The waveguide in Fig. 1a has a cross-section of w × t = 20 nm × 5 nm and forms a circle of radius R = 200 nm. We set ε wg = −50 for the upper half of the waveguide and ε wg = −50 + 50iγ wg (1 − 2|θ − π|/π),θ ∈ [π/2,3π/2] for the lower half. The waveguide in Fig. 2a has the same cross-section and ε wg = −50. In the 2D model system in Fig. 4a, the cylinder is made of metal and coated with a dielectric layer to guarantee a good dipole approximation. The permittivity of the metal is described by the Drude model ε ¼ 1 À ω 2 p = ω 2 þ iωγ ð Þwith ω p = 3.57 × 10 15 rad s −1 and γ = 3 × 10 12 rad s −1 . The shell has a relative permittivity of 9. The cylinders have an inner radius of 5 nm and an outer radius of 6 nm. The waveguide has a relative permittivity of ε wg = −10. The distance between the cylinders and the waveguide is h p = 25 nm. The perturbation cylinder is 20 nm above the waveguide and has a permittivity of −2 + iδ. δ controls the perturbation strength ϵ. The value of ϵ was determined by calculating the reflection coefficient of the perturbation cylinder. All of the numerical simulations were performed using COMSOL Multiphysics (www.comsol.com).
Microwave measurement. The spheres have diameters of 5.56 mm and are made of ZrO 2 with a relative permittivity of ε d = 33.7 + 0.06i. The waveguide has a crosssection of w × t = 9.9 mm × 3.0 mm and was cut from a commercial ECCOSTOCK dielectric slab with a relative permittivity of ε wg = 9 + 0.02i. The spheres were separated from the waveguide by a 4-mm-thick foam spacer. The probe for detecting magnetic field has a split-ring shape with a radius of 1.5 mm. The detecting point is located 1 mm above the spheres. The ends of the waveguide were inserted into foam absorbers to realise the absorbing boundary condition.