Light transfer transitions beyond higher-order exceptional points in parity-time and anti-parity-time symmetric waveguide arrays

We propose two non-Hermitian arrays consisting of $N=2l+1$ waveguides and exhibiting parity-time ($\mathcal{PT}$) or anti-$\mathcal{PT}$ symmetry for investigating light transfer dynamics based on $N$th-order exceptional points (EPs). The $\mathcal{PT}$-symmetric array supports two $N$th-order EPs separating an unbroken and a broken phase with real and imaginary eignvalues, respectively. Light transfer dynamics in this array exhibits radically different behaviors, i.e. a unidirectional oscillation behavior in the unbroken phase, an edge-towards localization behavior in the broken phase, and a center-towards localization behavior just at $N$th-order EPs. The anti-$\mathcal{PT}$-symmetric array supports also two $N$th-order EPs separating an unbroken and a broken phase, which refer however to imaginary and real eigenvalues, respectively. Accordingly, light transfer dynamics in this array exhibits a center-towards localization behavior in the unbroken phase and an origin-centered oscillation behavior in the broken phase. These nontrivial light transfer behaviors and their controlled transitions are not viable for otherwise split lower-order EPs and depend on the underlying $SU(2)$ symmetry of spin-$l$ matrices.

Hinging on the mathematical isomorphism between the Schrödinger equation and the Maxwell paraxial wave equation, investigations on non-Hermitian physics have been extended to and are flourishing in optics and photonics in view of the great convenience on experimental realizations [45][46][47][48][49][50]. Typical optical platforms for investigating non-Hermitian phenomena include arrays of resonators and waveguides described by engineered mode-coupling matrices with complex diagonal or offdiagonal elements [51][52][53]. It is found, for instance, light transport in two-mode PT -symmetric waveguides with balanced gain and loss changes drastically from a reciprocal to an nonreciprocal behavior across the 2ndorder EPs of mode-coupling matrices [45,46]. More interesting non-Hermitian characteristics and applications, such as richer topological photonics [54], enhanced spontaneous emission [55], improved sensitivity of spectra and dynamics [56][57][58], and delayed sudden death of entanglement [59], have also been found around higher-order EPs of mode-coupling matrices, e.g. in PT -symmetric waveguide arrays discussed in terms of underlying SU (2) symmetry [42]. Meanwhile, passive optical systems exhibiting anti-PT symmetry have spurred a rapidly growing interest in that they hold unique non-Hermitian features [60][61][62][63] yet without the need of gain [64][65][66] limiting realistic implementations and applications. In this regard, multi-mode anti-PT -symmetric waveguides deserve further studies on higher-orders EPs and nontrivial light transport.
In this work we consider two non-Hermitian arrays consisting of N = 2l + 1 waveguides in the regimes of PT symmetry and anti-PT symmetry, respectively, for achieving two highest N th-order EPs and examining a few nontrivial light transfer behaviors. Mode-coupling matrices with imaginary diagonal (sub-diagonal) elements are constructed here for describing the PT -FIG. 1: A PT -symmetric array of 2l + 1 single-mode waveguides arranged (extending) along the x-axis (z axis). Each mode an exhibits an imaginary propagation constant −in∆ indicating balanced gain (n < 0) and loss (n > 0) for ∆ > 0 while two real coupling constants Γn and Γn−1 (Γn+1 and Γn) with its neighbors an+1 and an−1, respectively, for n > 0 (n < 0). symmetric (anti-PT -symmetric) waveguide array on the basis of two spin-l matrices [43,44] embodying the underlying SU (2) symmetry. Different from the ballistic (breathing-like) light transport in a Hermitian waveguide array, the light transport in our PT -symmetric waveguide array can be switched between a unidirectional oscillation behavior in the unbroken phase and an edgetowards localization behavior (related to the skin effect [67]) in the broken phase through a center-towards localization behavior at one N th-order EP. Such a sensitive transition of light transfer behaviors becomes impossible when both N th-order EPs are split into a series of lower-order EPs by modifying somehow the corresponding mode-coupling matrix. Compared to the PT -symmetric waveguide array with balanced gain and loss, the anti-PT -symmetric waveguide array with passive couplings [65,[68][69][70] supports another transition between distinct light transfer behaviors around one N thorder EP: a center-towards localization behavior even in the unbroken phase and an origin-centered oscillation behavior in the broken phase.

II. PT -SYMMETRIC WAVEGUIDE ARRAY
In this section, we start by considering in Fig. 1 an array consisting of N = 2l + 1 single-mode waveguides labeled by a n for n ∈ {−l, l}. Our goal is to examine nontrivial light transfer behaviors based on N th-order EPs by constructing a PT -symmetric mode-coupling matrix for this N -mode waveguide array with underlying SU (2) symmetry. Generally speaking, N th-order EPs are viable to be observed via a fine tuning of propagation and coupling constants, though may be easily split into lowerorder EPs even in case of a slightly broken N th-order SU (2) symmetry. In principle, two N th-order EPs will be found at α = ±β for an N × N non-Hermitian matrix H non = αS i + iβS j (i = x and j = z or vice versa), being S i (S j ) the Cartesian operator along the i-axis (j-axis) of a spin-l system [44]. Such non-Hermitian matrices have been adopted to describe various quantum or classical systems, including the Bose-Hubbard model [43], ultracold Bose gases [71], and waveguide arrays [58].
With above considerations, we construct the PTsymmetric mode-coupling matrix with imaginary diagonal elements H n,n = −in∆ for n ∈ {−l, l} and real sub-diagonal elements H n±1,n = H n,n±1 = Γ n restricted by Γ n = Γ (l ± n + 1)(l ∓ n) for n ∈ ±{0, l − 1}. With respect to the waveguide array as illustrated in Fig. 1 [46], imaginary diagonal elements refer to gain or loss constants adjusted, e.g. by a temperature gradient [72,73], while real sub-diagonal elements refer to adjacent couplings adjusted, e.g. by the waveguide separations [74,75]. To be more specific, the waveguide array described by Eq. involutory permutation matrix while the time reversal operator T indicates complex conjugation with H † = T H T T = T HT . We plot real parts in Fig. 2(a) while imaginary parts in Fig. 2(b) of eigenvalues λ n of PT -symmetric modecoupling matrix H in Eq. (1) with l = 7. It is easy to see that all fifteen eigenvalues coalesce at ∆/Γ = ±2 with λ n ≡ 0, yielding thus two highest 15th-order EPs marked by Υ l and Υ r . It is more important that the PT -symmetric waveguide array can be engineered on demand to work in the unbroken phase with purely real eigenvalues in the case of |∆/Γ| < 2, while in the broken phase with purely imaginary eigenvalues in the case of |∆/Γ| > 2. To clarify the importance of choosing proper sub-diagonal elements, we then make similar plots in Fig. 2(c) and Fig. 2(d) for another PT -symmetric mode-coupling matrix where Γ n = Γ (l ± n + 1)(l ∓ n) is replaced by Γ n = Γ for n ∈ ±{0, l − 1}. It is found that each 15th-order EP degrades into a dozen of 2nd-order EPs and a single 3rd-order EP. To be more specific, points Υ l and Υ r are split into , respectively, in a way of exact mirror symmetry with respect to ∆/Γ = 0. This then results in a very narrow unbroken phase between points A ± l and A ± r ; twelve partially broken phases between points G ± l and A ± l as well as between points A ± r and G ± r ; two fully broken phases at the left of points G ± l or at the right of points G ± r . As a light beam is incident upon one or more of the N waveguides, its longitudinal propagation along the z-axis and transverse transfer between adjacent waveguides evolve according to a Schrödinger-like equation i∂ z ψ = Hψ [45,46]. Here an N × 1 column vector ψ = (a −l , ..., a 0 , ..., a l ) T has been introduced to describe the collective excitation of this PT -symmetric waveguide array, being a n an excitation amplitude of the nth mode. Using the PT-symmetric mode-coupling matrix H in Eq. (1), it is straightforward to expand the Schrödingerlike equation into restricted by the boundary conditions n ∈ {−l, l} and Γ −l−1 = Γ l+1 = 0. We examine with Eq. (3) evolutions of mode excitations characterized by |a n (z)| 2 / n |a n (z)| 2 in the PTsymmetric waveguide array for a single-site initial excitation a n (0) = δ n,0 . Plotted with ∆/Γ = 1 in the unbroken phase, Fig. 3(a) shows that the incident light first transfers unidirectionally toward waveguides of larger gain constants until reaching the left edge, then oscillates symmetrically between the left and right edges with a period δz = 2π/ √ 4Γ 2 − ∆ 2 . This is different from the (breathing-like) bidirectional Bloch oscillations observed in corresponding Hermitian waveguide arrays [72] because mode-coupling matrix H exhibits real eigenval- ues and linearly gradient gain/loss constants. In other words, the unidirectional oscillation occurs only if an excitation transfer from the central (n = 0) waveguide is increasingly evanescent toward the right (n > 0) waveguides while increasingly amplified toward the left (n < 0) waveguides in an exactly asymmetric manner. Plotted with ∆/Γ = 2 at one N th-order EP Υ r , Fig. 3(b) shows that mode excitations are localized around the central waveguide with no oscillations after a single reflection from the left edge. This may be attributed to a strong centripetal force arising from the highest-degree degeneracy of all eigenvalues λ n ≡ 0 and promises a robust directional output of mode excitations from the interface of middle waveguides with loss and gain respectively. Plotted with ∆/Γ = 4 in the broken phase, Fig. 3(c) shows that mode excitations suffer neither oscillations nor reflections and are localized near the left edge instead, similar to the skin effect. This promises another robust directional output of mode excitations clinging to the left edge and can be understood by considering that complex eigenvalues of mode-coupling matrix H lead to evanescent evolutions in general, leaving only amplified supermodes on waveguides of largest gain [45]. Similar light transfer behaviors have been discussed in [45], but it is worth noting that mode excitations in our work are confined in a much smaller number of waveguides, yielding thus much tighter light transfer behaviors facilitating potential applications, applications, as spin-l matrices S x and S z are used to construct mode-coupling matrix H.
Light transfer behaviors in the PT -symmetric waveg-uide array will become quite different, provided that each N th-order EP is split into a series of lower-order EPs as shown in Fig. 2(c) and Fig. 2(d). This is examined in Fig. 3(d) for an unbroken phase between 2nd-order EPs A ± l and A ± r ; in Fig. 3(e) at 2nd-order EPs E ± l separating two partially broken phases; in Fig. 3(f) for a fully broken phase at the right of 2nd-order EPs G ± r . We can see that evolutions of mode excitations vary from an asymmetric (breathing-like) bidirectional oscillation behavior, through an imperfect left-towards localization behavior, to a perfect left-towards localization behavior when ∆/Γ is increased to across one by one the 2nd and 3rd-order EPs. This nontrivial transition of light transfer behaviors depends critically on the order (degree of degeneracy) and the number of EPs. It is worth noting that the center-towards localization behaviors observed in Fig. 3(b) cannot be observed again in the absence of highest N th-order EPs.
Each N th-order EP in Fig. 2(a) and Fig. 2(b) for the PT -symmetric waveguide array indicates also a sensitive transition with respect to distinct light transfer behaviors. This is justified in Fig. 4 by examining |a 0 (z)| 2 /Σ n |a n (z)| 2 for a few typical ratios of ∆/Γ and another single-site initial excitation a n = δ n,7 . It is easy to see that the central waveguide exhibits a robust mode excitation with neither periodic oscillations nor a fast decay for ∆/Γ = 2 referring to one N th-order EP Υ r . Adjusting ∆/Γ to enter the unbroken or broken phase beyond this N th-order EP, we find instead periodic oscillations or a fast decay of mode excitation in the central waveguide. Hence it is viable to attain a robust output of mode excitation from the central waveguide or realize a sensitive switching between distinct light transfer behaviors in the PT -symmetric waveguide array by accurately modulating ∆ or Γ.

III. ANTI-PT -SYMMETRIC WAVEGUIDE ARRAY
In this section, we try to construct an anti-PTsymmetric mode-coupling matrix so as to further investigate highest N th-order EPs and nontrivial light transfer behaviors in an alternative array consisting of N = 2l + 1 single-mode waveguides. This can be done by rotating both S x and S z in the PT -symmetric mode-coupling matrix H around S y , yielding thus S x → S z and S z → −S x . The resultant anti-PT -symmetric mode-coupling matrix is which reserves the underlying SU (2) symmetry on one hand and satisfies the anti-commutation relation {PT , H} = 0 on the other hand. Obviously, the 2l + 1 diagonal elements H n,n become real so that all waveguide modes a n suffer neither loss nor gain, while the 4l sub-diagonal elements H n,n±1 become imaginary so that waveguide modes a n and a n±1 are passively coupled. However, passive couplings are usually not realistic, hence one question arises: is it viable to realize modecoupling matrix H by modifying the waveguide array as shown in Fig. 1?
We plot in Fig. 6(a) and Fig. 6(b), respectively, real and imaginary parts of eigenvalues λ n of the anti-PTsymmetric matrix H in Eq. (4) with l = 7. In comparison to that observed in Fig. 2(a) and Fig. 2(b), the common feature is that all fifteen eigenvalues coalesce again at ∆/Γ = ±2 with λ n ≡ 0 to yield two highest 15th-order EPs Υ l and Υ r ; the difference lies in that the unbroken and broken phases refer instead to purely imaginary and purely real eigenvalues, respectively [68]. To be more specific, the unbroken (broken) phase is attained with a smaller (larger) absolute ratio of imaginary propagation constant −in∆ and real coupling constant Γ n for H in Eq. (1), while with a smaller (larger) absolute ratio of real propagation constant −n∆ and imaginary coupling constant −iΓ n for H in Eq. (4). Therefore, distinct light transfer behaviors are expected to occur in the PT -symmetric and anti-PT -symmetric waveguide arrays even if they are engineered to work with the same parameters as discussed below.
We then examine with Eq. (8) evolutions of mode excitations characterized by |a n (z)| 2 / n |a n (z)| 2 in the anti-PT -symmetric waveguide array for different singlesite initial excitations a n (0) = δ n,−7 , a n (0) = δ n,0 , and a n (0) = δ n,7 as shown in Fig. 7. In the unbroken phase with ∆/Γ = 1, we can see from Figs. 7(a-c) that mode excitations always converge to and finally leave from a few middle waveguides without exhibiting oscillations and independent of initial single-site excitations. Such a robust center-towards localization behavior is clearly different from the left-edge-towards localization behavior observed in the broken phase of an PT -symmetric waveguide array as shown in Fig. 3(c). The underlying physics for this center-towards localization behavior should be that all adjacent waveguides are passively coupled with the gradually reduced imaginary coupling constants −iΓ n = −iΓ −n from both edges to the center. In the broken phase with ∆/Γ = 3, we can see from Figs. 7(d-f) that mode excitations oscillate with an identical fixed period determined by the ratio ∆/Γ but of different centers identical to origins of single-site excitations. This more localized origin-centered oscillation may be attributed to a larger ratio ∆/Γ = 3 and is also very different from that observed in the unbroken phase of an PT -symmetric waveguide array as shown in Fig. 3(a).
To gain a more intuitive insight into the transition between different light transfer behaviors, we focus on the central (n = 0) waveguide by examining |a 0 (z)| 2 /Σ n |a n (z)| 2 in Fig. 8 for different values of ∆/Γ while a fixed initial excitation a n = δ n,0 . As the anti-PTsymmetric waveguide array works at and below its N thorder EP Υ r , consistent with that observed in Fig. 7(b) but different from that observed in Fig. 4, the central mode excitation quickly decays to a steady value ∼ 0.5 independent of the ratio ∆/Γ. As the anti-PT -symmetric waveguide array works above the N th-order EP Υ r , consistent with that observed in Fig. 7(e) but different from that observed in Fig. 4, the central mode excitation starts to oscillate with a period δz = 2π/ √ ∆ 2 − 4Γ 2 tunable in an arbitrarily large range (unattainable yet in the PTsymmetric waveguide array) since ∆/2Γ may approach either unit or infinity. Hence, it is viable to attain a controlled output of mode excitation from the central waveguide by adjusting the ratio ∆/Γ, which may serve as an efficient scheme of high-precision measurement on relevant parameters determining ∆ or Γ.

IV. CONCLUSIONS
In summary, we have studied nontrivial light transfer dynamics by working at and beyond one of two N thorder EPs, referring to highest-degree degeneracies of all N eigenvalues of relevant mode-coupling matrices, in two non-Hermitian waveguide arrays exhibiting a PT or an anti-PT symmetry. For the PT -symmetric waveguide array with balanced gain and loss, light transfer dynamics is viable to switch between a unidirectional oscillation behavior in the unbroken phases, through a center-towards localization behavior at one N th-order EP, to an edgetowards localization behavior in the broken PT phases. For the anti-PT -symmetric waveguide array with passive couplings, light transfer dynamics is very different again in the broken and unbroken phases, hence supporting the transition between an origin-centered oscillation behavior and a center-towards localization behavior across one N th-order EP. These drastic transitions of light transfer behaviors are expected to facilitate advanced photonic applications like sensitive optical switching and high-precision measurement, while providing also new perspectives for further theoretical studies.