2
mon in non-Hermitian systems, but the acquired phases
can be largely enhanced by encircling DPs or EPs [37–39].
Moreover, DPs and EPs are useful in testing and classify-
ing phases and phase transitions [40,41]. For example, a
Liouvillian spectral collapse in the standard Scully-Lamb
laser model occurs at a quantum DP [42,43].
Recent studies on EPs have also shown that by ex-
ploiting a nontrivial topology in the vicinity of EPs in
the energy spectrum can lead to a swap-state effect,
where the initial state does not come back to itself af-
ter a round trip around an EP. Such phenomenon has
been predicted theoretically [44,45] and observed exper-
imentally in [21,37,46–48], while performing ‘static’,
i.e., independent, measurements at various locations in
the system parameter space. However, when encircling
an EP dynamically, another intriguing effect can be in-
voked; namely, a chiral mode behavior, such that a start-
ing state, after a full winding period, can eventually re-
turn to itself [49–52]. The latter effect stems from the
breakdown of the adiabatic theorem in non-Hermitian
systems [49,53]. This asymmetric mode switching phe-
nomenon has also been experimentally confirmed in var-
ious platforms [38,54–58]. A number of studies have
demonstrated the practical feasibility to observe the chi-
ral light behavior on a pure quantum level [59] and even
in a so-called hybrid mode [60], where by exploiting var-
ious measurement protocols, one can switch between the
system dynamics described by a quantum Liouvillian and
the corresponding classical-like effective NHH.
Other works, both theoretical [61] and experimen-
tal [62], have pointed that a crucial ingredient in detect-
ing a dynamical flip-state asymmetry is the very curved
topology near EPs. In other words, it is not necessary
to wind around EPs in order to observe such phenom-
ena. However, the dynamical contours must be in a close
proximity to EPs [61].
More recently, much effort is put on studying the be-
havior of modes while encircling high-order or multiple
EPs in a parameter space of multimode systems. In-
deed, the presence of high-order or multiple low-order
EPs in a system spectrum, along with the non-Hermitian
breakdown of adiabaticity, can impose a substantial diffi-
culty to manipulate the mode-switching behavior on de-
mand [52,63,64]. That is, a system may end up only
in a few states out of many regardless of the encircling
direction and winding number.
In this work we demonstrate that dynamically wind-
ing around exceptional curves (ECs), whose trajecto-
ries can additionally cross diabolic curves (DCs), pro-
vides a feasible route to realize a programmable multi-
mode switch with controlled mode chirality. We use a
four-mode parity-time (PT )-symmetric bosonic system,
which is governed by an effective NHH, as an exemplary
platform to demonstrate this programmable switch. At
the crossing of EC and DC a new type of a spectral sin-
gularity is formed, referred to as diabolically degener-
ate exceptional points (DDEPs) [65]. By exploiting the
presence of DDEPs in dynamical loops of the system pa-
rameter space, one can restore the swap-state symmetry,
which breaks down in two-mode non-Hermitian systems.
This implies that the initial state can eventually return to
itself after a state flip in a double cycle. In other words,
the interplay between the topologies of EPs and DPs en-
ables one to restore (impose) mode symmetry (asymme-
try) on demand. These results are valid also for purely
dissipative systems (i.e., loss only systems without gain)
and can be extended to arbitrary multimode systems.
RESULTS
Theory
We start from the construction of a four-mode NHH,
possessing both exceptional and diabolic degeneracies.
For this, we follow the procedure described in [65], where
one can construct a matrix, whose spectrum is a combi-
nation of the spectra of two other matrices by exploiting
Kronecker sum properties. Namely, by taking two PT -
symmetric matrices
M1=i∆k
k−i∆, M2=0g
g0,(1)
one can form a PT -symmetric 4×4non-Hermitian matrix
H=M1⊗I+I⊗M2,(2)
where Iis the 2×2identity matrix. Explicitly, the matrix
Hreads
H=
i∆g k 0
g i∆ 0 k
k0−i∆g
0k g −i∆
.(3)
The symbols in Eq. (3) can have various physical mean-
ings, but in our context they may denote, e.g., cou-
pling (g, k) and dissipation (∆) strengths in a pho-
tonic system (see the text below). The PT -symmetry
operator is expressed via the parity operator P=
antidiag[1,1,1,1] and the time-reversal operator T, thus,
implying PT HPT −1=H. The matrix Hcan be related
to a linear four-mode NHH operator ˆ
H, written in the
mode representation, i.e.,
ˆ
H=Xˆa†
jHˆak,
where ˆai(ˆa†
i) are the annihilation (creation) operators
of bosonic modes i= 1,...,4. Such an NHH can be
associated, e.g., with a system of four coupled cavities or
waveguides (see Fig. 1a). A similar scheme, based on two
lossy and two amplified subsystems, has been proposed