AAPPS bulletin

Research and Review

Encircling the Liouvillian exceptional points: a brief review

writerKonghao Sun & Wei Yi

Vol.34 (Aug) 2024 | Article no.22 2024

Abstract

Exceptional points are the branch-point singularities of non-Hermitian Hamiltonians and have rich consequences in open-system dynamics. While the exceptional points and their critical phenomena are widely studied in the non-Hermitian settings without quantum jumps, they also emerge in open quantum systems depicted by the Lindblad master equations, wherein they are identified as the degeneracies in the Liouvillian eigenspectrum. These Liouvillian exceptional points often have distinct properties compared to their counterparts in non-Hermitian Hamiltonians, leading to fundamental modifications of the steady states or the steady-state-approaching dynamics. Since the Liouvillian exceptional points widely exist in quantum systems such as the atomic vapors, superconducting qubits, and ultracold ions and atoms, they have received increasing amount of attention of late. Here, we present a brief review on an important aspect of the dynamic consequence of Liouvillian exceptional points, namely the chiral state transfer induced by the parametric encircling the Liouvillian exceptional points. Our review focuses on the theoretical description and experimental observation of the phenomena in atomic systems that are experimentally accessible. We also discuss the ongoing effort to unveil the collective dynamic phenomena close to the Liouvillian exceptional points, as a consequence of the many-body effects therein. Formally, these phenomena are the quantum-many-body counterparts to those in classical open systems with nonlinearity, but hold intriguing new potentials for quantum applications.

1 Introduction

In a typical non-Hermitian system, the dynamics is effectively driven by a non-Hermitian Hamiltonian [1]. As a common feature of a wide class of non-Hermitian matrices, multiple eigenvectors and eigenvalues of a non-Hermitian Hamiltonian can simultaneously coalesce at certain critical points in the parameter space. These critical points, known as the exceptional points (EPs), correspond to the so-called branch-point singularities in the eigenspectrum, where the Hamiltonian cannot be diagonalized [2,3,4]. As such, the EPs are fundamentally different from conventional degeneracies in Hermitian systems. In recent years, the EPs have been widely discussed in the context of non-Hermitian models with the parity-time symmetry [5,6,7,8,9,10,11,12,13], where they appear as the demarcation between a spectral region with entirely real eigenvalues (dubbed the parity-time unbroken regime), and one without (dubbed the parity-time broken regime). Being the transition point between these spectrally and dynamically distinct regions, the EPs have been extensively studied in the context of parity-time symmetric models in connection with phenomena such as power oscillations [14,15,16], directional transport [17, 18], and lasing [19,20,21,22]. However, EPs also exist in systems without the parity-time symmetry and are therefore more general [1, 10]. Thanks to the singular nature of EPs, a system exhibits many interesting behaviors in their vicinity, including universal criticality [23,24,25,26], non-reciprocal dynamics [27,28,29,30,31,32,33,34,35,36,37,38], enhanced entanglement generation [39, 40], and strong sensitivity to external perturbations [41,42,43]. These properties have been confirmed in a wide range of classical and quantum mechanical systems, including optics and photonics [29, 30], optomechanics [44, 45] acoustics [46,47,48], atomic gases [33, 49], and superconducting qubits [40, 50, 51]. Understanding EPs and the related phenomena create opportunities for developing useful applications in terms of sensing and topological transport.

In the quantum regime, EPs, and non-Hermitian physics in general, are mostly discussed in the context of conditional dynamics of an open quantum system. Specifically, the full quantum dynamics of an open system is governed by the Lindblad master equation under the Markov approximation [52]. In the quantum trajectory picture, the density matrix evolution can be unraveled as an ensemble of quantum trajectories [53]. The quantum state in each trajectory undergoes evolution driven by an effective non-Hermitian Hamiltonian, but is interrupted by stochastic quantum jumps. It follows that, when considering the transient dynamics or by selecting trajectories without quantum jumps (hence conditioned), the dynamics is governed purely by the non-Hermitian Hamiltonian [1]. While such a scheme can be experimentally implemented through post selection, the practice is limited to single-particle systems, or non-interacting ones where quantum statistics is not essential.

On the other hand, the Lindblad master equation, complete with full quantum jump processes, describes the evolution of the density matrix under the Liouvillian superoperator. Since the Liouvillian itself can be represented as a non-Hermitian matrix acting on a vectorized density matrix, the EPs can also arise in the Liouvillian spectrum [54,55,56,57], leading to unique dynamics in the full quantum dynamics of the master equation. Along this line of thinking, the Liouvillian EPs are defined and studied in recent theoretical and experimental works [50, 5154, 56,57,58,59,60]. In Fig. 1, we illustrate two recent experiments where LEPs were engineered in superconducting qubits and trapped ions, respectively. It is found that the Liouvillian EPs have distinct impact on the open-system dynamics, compared to their counterparts in non-Hermitian Hamiltonians. More importantly, since the application of the quantum master equation is not limited to non-interacting systems or transient dynamics, the Liouvillian EPs can become relevant in open systems with many-body correlations [61,62,63,64]. Thus, understanding Liouvillian EPs in a generic open quantum system exemplifies the advancing frontier of non-Hermitian physics, where the wealth of non-Hermitian phenomena is no longer restrained to the classical or semi-classical regime, but lends meaningful insights to quantum many-body processes.

Fig. 1
figure 1

a Schematics of LEP engineering in the superconducting qubit experiment [51]. b The corresponding Liouvillian exceptional structure in the parameter space of Ref. [51]. c Level scheme for generating LEP in the trapped ion experiment [59]. b The corresponding LEP in the parameter space of [59]. Here, (a and b) are adapted from Ref. [51], copyright © 2022 by the American Physical Society. And (c and d) are adapted from Ref. [59], copyright © by the American Physical Society


In the following, we provide a brief review on the study of Liouvillian EPs, focusing on the experimentally relevant atomic systems. We first review conventional EPs in non-Hermitian Hamiltonians, using a simple two-level system as an example. Then, by adding different quantum jump operators, we discuss how the density-matrix dynamics admits a non-Hermitian description as well as the emergence of Liouvillian EPs in the framework of Lindblad master equations. As the main content of the review, we give some concrete examples, all experimentally related, to the role of Liouvillian EPs in open-system dynamics, and show how non-Hermitian physics emerge within the framework of master equations [33, 50, 51, 58, 65]. Finally, we discuss a recent experiment on the collective topological state transfer near a Liouvillian EP, where both non-Hermiticity and many-body effects play a key role.

2 Hamiltonian EP and Liouvillian EP

To illustrate the connection and difference between the EPs in non-Hermitian Hamiltonians and those of the Liouvillians, we first consider an exemplary two-level model. The Hamiltonian in the corresonding spin basis can be written as

\(\begin{aligned} H=J\sigma _x-i\frac{\Gamma }{2}\sigma _z, \end{aligned}\)
(1)

where J and \(\Gamma\) are positive coefficients. The Hamiltonian possesses the parity-time symmetry, with \((PT)H(PT)^{-1}=H\), where the parity operator \(P=\sigma _x\), and the time-reversal operator T corresponds to complex conjugation. The eigenvalues of the Hamiltonian are straightforward to calculate

\(\begin{aligned} E_{\pm }=\pm \sqrt{J^2-\frac{\Gamma ^2}{4}}. \end{aligned}\)
(2)

Here, the eigenvalues are real for \(J\ge \Gamma /2\) and imaginary for \(J<\Gamma /2\). Clearly, \(J=\Gamma /2\) corresponds to the aforementioned EP, separating the parity-time unbroken and broken regimes. To understand the criticality of the EP, we examine the eigenvectors for \(J\ge \Gamma /2\), with

\(\begin{aligned} |\psi _{\pm }\rangle =\frac{\sqrt{2}}{2J} \left( \begin{array}{c} -i\frac{\Gamma }{2}\pm \sqrt{J^2-\frac{1}{4}\Gamma ^2}\\ J \end{array}\right) . \end{aligned}\)
(3)

Notably, the eigenvectors coalesce at the EP \(J=\Gamma /2\), where the non-Hermitian matrix cannot be diagonalized. The ill-conditioned matrix leads to many intriguing observable phenomena, and have stimulated much research interest.

Likewise, exceptional structures also exist in the Liouvillian eigenspectrum. Based on Hamiltonian (1), we consider the following Lindblad master equation (we set \(\hbar\)=1)

\(\begin{aligned} \dot{\rho }=\mathcal {L}\rho =-i[H,\rho ]+L\rho L^\dagger -\frac{1}{2}L^\dagger L\rho -\frac{1}{2}\rho L^\dagger L, \end{aligned}\)
(4)

where \(H=J\sigma _x\), and the quantum-jump operator \(L=\sqrt{\Gamma } \sigma _-\). Note that when we drop the recycling term \(L\rho L^\dagger\), the remaining equation of motion describes the evolution of the density matrix by a non-Hermitian effective Hamiltonian \(H_{\text {eff}}=H-i\frac{\Gamma }{2}L^\dagger L\), which reproduces Hamiltonian (1). As we discuss in the following section, neglecting the recycling term is referred to as the no-jump condition, also known as the semi-classical limit or the post-selection condition. Though with limitations, it offers a convenient and practical route, connecting the dynamics of quantum open systems with those under a non-Hermitian effective Hamiltonian.

More relevant to the discussion here, we vectorize the density matrix

\(\begin{aligned} \rho =\left[ \begin{array}{cc} \rho _{11} & \rho _{12}\\ \rho _{21} & \rho _{22} \end{array}\right] \rightarrow \left( \begin{array}{c} \rho _{11}\\ \rho _{12}\\ \rho _{21}\\ \rho _{22} \end{array}\right) , \end{aligned}\)
(5)

and write the Liouvillian superoperator in a non-Hermitian matrix form

\(\begin{aligned} \mathcal {L}= \left[ \begin{array}{cccc} 0 & iJ & -iJ & \Gamma \\ iJ & -\frac{\Gamma }{2} & 0 & -iJ \\ -iJ & 0 & -\frac{\Gamma }{2} & iJ\\ 0 & -iJ & iJ & -\Gamma \\ \end{array}\right] . \end{aligned}\)
(6)

Importantly, now that the Liouvillian is represented by a non-Hermitian matrix, it can host exceptional structures. The exceptional points in the Liouvillian eigenspectrum are referred to as the Liouvillian exceptional points (LEPs), in contrast to the Hamiltonian exceptional points (HEPs). Specifically, the Liouvillian eigenspectrum \(\lambda\) is calculated through \(\mathcal {L}\rho _\lambda =\lambda \rho _\lambda\), with the corresponding Liouvillian eigenstates \(\rho _{\lambda }\). For the Liouvillian in (6), we then have

\(\begin{aligned} \lambda _1=0,\,\lambda _2=-\frac{\Gamma }{2},\,\lambda _{3,4}=\frac{1}{4}\left( -3\Gamma \pm \sqrt{\Gamma ^2-64J^2}\right) . \end{aligned}\)
(7)

While the eigenvalue \(\lambda _1=0\) corresponds to the steady-state solution, \(\lambda _{3,4}\) coalesce at \(\Gamma =8J\) and \(\lambda _{2,3}\) coalesce at \(J=0\), representing second-order LEPs. These LEPs lie away from the steady-state solution, suggesting that their impact is in the intermediate time scales, before the system relaxes to the steady state at long times. Note that, equivalent to the treatment above, one may also write down the corresponding optical Bloch equations, and extract information of the LEPs from the coefficient matrix [55].

Fig. 2
figure 2

The Liouvillian exceptional structure of Eq. (8) on the \(\delta\)J plane. The green lines are the second-order Liouvillian exceptional lines, the yellow stars indicate the third-order LEP. The red point is the second-order HEP of the corresponding non-Hermitian Hamiltonian under the no-jump condition. We take \(\Gamma =1\) as the unit of energy for our calculations


More generally, the connection and distinction between the LEP and HEP have been systematically discussed in Refs. [56, 57]. Although LEPs and HEPs have essentially different properties, an open system can have exceptional structures that are the direct correspondence to a HEP of the non-Hermitian Hamiltonian under the no-jump condition. This is the case with the Liouvillian above. The Liouvillian eigenspectrum features an exceptional point, but shifted in parameters compared to those of the HEP of the non-Hermitian Hamiltonian.

Furthermore, when a detuning term \(\frac{\delta }{2}\sigma _z\) is added to the Hermitian Hamiltonian H, the matrix form of the Liouvillian operator becomes

\(\begin{aligned} \mathcal {L}= \left[ \begin{array}{cccc} 0 & iJ & -iJ & \Gamma \\ iJ & -\frac{\Gamma }{2}-i\delta & 0 & -iJ \\ -iJ & 0 & -\frac{\Gamma }{2}+i\delta & iJ\\ 0 & -iJ & iJ & -\Gamma \\ \end{array}\right] . \end{aligned}\)
(8)

The corresponding Liouvillian eigenspectrum shows exceptional lines on the \(\delta\)J plane, ending at two third-order exceptional points as illustrated in Fig. 2. As such, the HEP of the non-Hermitian Hamiltonian develops into an exceptional structure in the Liouvillian eigenspectrum, consisting of exceptional lines and higher-order exceptional points.

3 Dynamic consequence of EPs

One outstanding feature of the EPs is the sensitivity of eigenvalues to external perturbations. For instance, by adding a weak perturbative term \(\epsilon \sigma _x\) to Hamiltonian (1), the eigenvalues at the EP become

\(\begin{aligned} E_{\pm }=\pm \sqrt{\epsilon (2J+\epsilon )}, \end{aligned}\)
(9)

splitting by a small amount \(\sim \epsilon ^{\frac{1}{2}}\). More generally, at a higher-order EP, where more than two eigenvectors and eigenvalues coalesce, the splitting is of the order \(\epsilon ^{\frac{1}{n}}\). Such a sensitivity is the basis for discussions of EP-enhanced sensing, which have been extensively reported and reviewed in recent years [10, 11, 41,42,43].

Here, we focus on a general dynamic consequence, which derives from the complex spectral topology in the parameter space near an EP. Take, for instance, a parameterized non-Hermitian model based on Eq. (1)

\(\begin{aligned} H(t)=J(t)\sigma _x-[\Omega (t)+i\frac{\Gamma }{2}]\sigma _z, \end{aligned}\)
(10)

where \(J(t)=\Gamma /2+r\cos (\omega t)\), \(\Omega (t)=r\sin (\omega t)\), and t is understood as a parameter but would eventually parameterize the time dependence of the variables. The eigenvalues are

\(\begin{aligned} E_{\pm }=\pm \sqrt{r^2+\Gamma r e^{i\omega t}}, \end{aligned}\)
(11)

which, in the parameter space of \((J,\Omega )\), lead to the Riemann sheets illustrated in Fig. 3. Here, an EP exists at \((J=\Gamma /2, \Omega =0)\), which represents the endpoint of the branch cut along \(\Omega =0\) with \(J\le 0\). As such, the EP is usually referred to as the branch-cut singularity of the non-Hermitian matrix.

Such a geometry has interesting dynamic consequences. Consider the evolution of an eigenstate under H(t) where t now represents time. Intuitively, when the rate of time variation is slow compared to the real energy gap along the path (roughly with \(\frac{2\pi }{\omega }\text {min}(|E_+-E_-|)\gg 1\)), the time-evolved state should adiabatically follow the instantaneous right eigenstate of H(t)

\(\begin{aligned} |\psi _{\pm }\rangle \propto \left( \begin{array}{c} -[\Omega (t)+i\frac{\Gamma }{2}]\pm \sqrt{J(t)^2+(\Omega (t)+i\frac{\Gamma }{2})^2}\\ J(t) \end{array}\right) . \end{aligned}\)
(12)

We track such a process with the trajectory \(\overline{E}(t)\) on the Riemann sheet, with

\(\begin{aligned} \overline{E}(t)=\frac{\sum _{i=\pm }|\langle \chi _i(t)|\psi (t)\rangle |^2 E_i(t)}{\sum _{i=\pm }|\langle \chi _i(t)|\psi (t)\rangle |^2}, \end{aligned}\)
(13)

where \(|\chi _i(t)\rangle\) is the left eigenstate of the instantaneous Hamiltonian, defined as \(H^{\dagger }(t)|\chi _i(t)\rangle =E^*_i(t)|\chi _i(t)\rangle\). Importantly, when the trajectory crosses the branch cut, the band indices (±) of the right eigenstate switch, due to the sign change in the front of the square root in Eq. (12). Correspondingly, assuming the system is initialized in an eigenstate, under one full cycle of the parameter change, it undergoes an eigenstate switching under the adiabatic condition. This is illustrated in Fig. 3a, where the parameters change in a counterclockwise direction along the trajectory.

Fig. 3
figure 3

Trajectories along the Riemann surface corresponding to the non-Hermitian model (10). The Riemann surface is colored red (blue), indicating the gain (loss) of the eigenvalues. Parameters change in (a) counterclockwise and (b) clockwise directions end up with different states after one cycle. For numerical calculations, we take the parameters \(\Gamma =1\), \(r=0.1\), and \(T=100\)


However, the adiabatic condition does not always hold. As illustrated in Fig. 3b, when the parameters change in a clockwise direction, the system comes back to the original eigenstate, due to a non-adiabatic jump that switches the eigenstate one more time along the path. Interestingly, such a non-adiabatic jump occurs even if the condition \(\frac{2\pi }{\omega }\text {min}(|E_+-E_-|)\gg 1\) is satisfied. Such a phenomenon derives from the non-Hermiticity of the system. Specifically, expanding the time-evolved state onto the basis of the instantaneous right eigenstates [31]

\(\begin{aligned} |\psi (t)\rangle =\sum \limits _{n=\pm } c_n(t)|\psi _n(t)\rangle , \end{aligned}\)
(14)

the coefficients satisfy

\(\begin{aligned} i\hbar \frac{\partial }{\partial t}c_n=\left[E_n(t)-i\hbar \left\langle \chi _n\left|\frac{\partial }{\partial t}\right|\psi _n\right\rangle \right]c_n-i\hbar \left\langle \chi _n\left|\frac{\partial }{\partial t}\right|\psi _{-n}\right\rangle c_{-n}. \end{aligned}\)
(15)

The last term on the right-hand side describes the non-adiabatic transition between different eigenstates. The term can have significant impact when an unoccupied eigenstate has an eigenvalue with a positive imaginary component (or simply a less negative imaginary component), leading to an exponentially increased likelihood of a non-adiabatic jump into the said state. As a result, depending on the initial state and the encircling direction, the system undergoes an eigenstate switch only along one direction of the parameter change. This is referred to as the chiral state transfer in the literature [10, 11]. While the process of the parameter change is known as the EP encircling, the path of the parameter change does not have to encircle the EP: the chiral state transfer can occur when the path lies close to the exceptional structure. Exceptions can arise under special circumstances, where the chiral state transfer can occur on trajectories far way from any EPs [32, 33].

4 Implementing non-Hermitician Hamiltonian in quantum systems

The chiral state transfer under the Hamiltonian EP encircling has been experimentally demonstrated in a variety of physical systems [33,34,35,36,37,38, 50, 51]. A recent experiment further establishes the phenomenon in an ultracold gas of fermions [33]. While in all cases, the dynamics are driven by a non-Hermitian effective Hamiltonian, in quantum systems, the no-jump condition leading to the non-Hermiticity is typically imposed through post selection.

Fig. 4
figure 4

a The level scheme in Ref. [33], where states in the \(^1S_0\) and \(^3P_1\) manifolds of ultracold \(^{173}\)Yb atoms are used. b Level scheme for a simplified four-level model


Take the cold-atom experiment for example. Ultracold \(^{173}\)Yb atoms are prepared in the \(|m_F=3/2\rangle\) and \(|m_F=5/2\rangle\) states of the ground-state \(^1S_0\) (\(F=5/2\)) manifold. The two states are coupled through a two-photon Raman process, and the state \(|m_F=3/2\) is further laser coupled to an electronically excited state in the \(^3P_1\) (\(F=7/2\)) manifold. This introduces a controlled atom loss to the system, due to the finite linewidth of the excited state. The level scheme is illustrated in Fig. 4a. Neglecting for now the spontaneous emission back to the \(\{|m_F=3/2\rangle , |m_F=5/2\rangle \}\) subspace, one assumes that the atoms undergo the spontaneous emission process are lost from the system. Under such a condition, the dynamics of the remaining atoms are governed by a non-Hermitian effective Hamiltonian

\(\begin{aligned} H_{\text {eff}}=\frac{\delta }{2}\sigma _z-\Omega \sigma _x-i\frac{\Gamma }{4}(1-\sigma _z), \end{aligned}\)
(16)

where \(\Omega\) and \(\delta\) are respectively the Rabi frequency and detuning of the Raman process, and \(\Gamma\) is the laser-induced loss rate. For the convenience of discussion, we label \(|1\rangle =|m_F=5/2\rangle\) and \(|2\rangle =|m_F=3/2\rangle\) and define the Pauli operators such that \(\sigma _z=|1\rangle \langle 1|-|2\rangle \langle 2|\). We also label the excited state used for the laser-induced loss as \(|3\rangle\), as illustrated in Fig. 4b.

From the perspective of quantum open systems, the non-Hermitian effective Hamiltonian above can be derived by imposing a no-jump condition on the quantum master equation

\(\begin{aligned} \dot{\rho }=-i[H,\rho ]+L_{1}\rho L_{1}^{\dagger }-\frac{1}{2}L_{1}^\dagger L_{1}\rho -\frac{1}{2}\rho L_{1}^\dagger L_{1}, \end{aligned}\)
(17)

where \(H=\frac{\delta }{2}\sigma _z-\Omega \sigma _x\), and \(L_1=\sqrt{\Gamma }|4\rangle \langle 2|\) describes the laser-induced loss from state \(|2\rangle\) to a bystander state \(|4\rangle\). The no-jump condition corresponds to neglecting the recycling term \(L_1\rho L_1^\dagger\), which applies when one focuses only on the subspace \(\{|1\rangle ,|2\rangle \}\). Equivalently, from all possible time-evolved states, one selects out only those with no support on \(|4\rangle\). In practice, this is natural if one assumes that atoms in \(|4\rangle\) are no longer trapped while only those remain trapped are detected. Since the atoms that remain in the \(\{|1\rangle , |2\rangle \}\) subspace necessarily have not gone through the quantum jump (spontaneous emission) to state \(|4\rangle\), the non-Hermitian description is therefore always valid within the subspace. The exponential decrease in the modulus square of the wave function directly corresponds to the exponential decay of the total atom number within the subspace. It is worth pointing out that such a natural implementation of post selection (by detection) only applies when the atoms are non-interacting. Even as the recycling term mixes density matrices from sectors with different atom numbers, the overall density matrix in each sector is similar in structure, featuring a direct product of those of individual atoms. This is not the case in the presence of interactions, under which density matrices in different atom-number sectors are coupled through the interplay of dissipation and interaction.

Under the condition above, the master equation reduces to the non-Hermitian description (16), which, having the same structure as Hamiltonian (10), forms the basis of the observed EP encircling in Ref. [33]. Therein, as the system’s parameters vary along a closed path, the internal atomic states flip but only one way around. The parameters are chosen such that all the trapped atoms encircle a common EP in a similar fashion, leading to a collective chiral state flip.

In another more commonly discussed scenario, the jump operator is replaced by \(L^\prime _1=\sqrt{\Gamma }|1\rangle \langle 2|\), meaning the spontaneous emission from state \(|2\rangle\) to \(|1\rangle\). Such a decay channel can be engineered in the current setup, for instance, by coupling \(|2\rangle\) to an excited state which decays back to state \(|1\rangle\). One can still recover the non-Hermitian description by following the no-jump condition. But as the final state of the spontaneous emission remains in the \(\{|1\rangle , |2\rangle \}\) subspace, in practice, the no-jump condition becomes more difficult to satisfy for an ensemble of atoms. Specifically, the timescale within which the non-Hermitian Hamiltonian reins decreases exponentially with the total atom number, as any single jump would leave the corresponding atom in a mixed state and beyond the non-Hermitian description. The accumulation of these events would effectively heat up the system, and eventually make the non-Hermitian description invalid even on an approximate level. This is why such a scheme is typically adopted for single- or few-qubit systems, but difficult to implement in atom gases.

5 Encircling the LEP

However, the reality is more complicated than the simple picture above. While the non-Hermiticity is introduced through the laser-assisted atom loss, the quantum jump operator \(L_1\) is not the only spontaneous emission process that can happen. When the atoms are in the excited \(^3P_1\) state, there is a finite probability for it to decay back to the original state (\(|2\rangle\) in our convention). After adiabatic elimination of the excited state, the process gives rise to an additional jump operator \(L_2=\sqrt{\gamma }|2\rangle \langle 2|\), which corresponds to the dephasing between states \(|1\rangle\) and \(|2\rangle\), due to the spontaneous emission from the excited state back to \(|2\rangle\).

Fig. 5
figure 5

a Liouvillian eigenspectra and exceptional structure on the \(\delta -J\) plane of (18). The green lines are the second-order Liouvillian exceptional lines and the yellow stars indicate the third-order LEP. b, c Trajectories against the Liouvillian eigenspectra of (18) in the adiabatic limit for (b) clockwise and (c) counterclockwise rotations. de Trajectories against the Liouvillian eigenspectra for an intermediate encircling time. The encircling path is as follows: \(\delta (t)=0.5\sin {(\pm 2\pi t/T+2\pi /3)}\), and \(J(t)=0.5+0.5\cos {(\pm 2\pi t/T+2\pi /3)}\). The eigenspectra are colored red (blue), indicating the quasi-steady (excited) nature of the corresponding Liouvillian eigenstates. We take \(T=150\) for (a and b) and \(T=10000\) for (and d). Other parameters are \(\Gamma =1/20\), \(\gamma =1/100\)


The dephasing process further restricts the applicability of the non-Hermitian Hamiltonian. As discussed in the previous section, under \(L_1\), it is sufficient to implement a no-jump condition by detecting (post selecting) atoms that remain in the trap (and hence in the \(\{|1\rangle , |2\rangle \}\) subspace). With the addition of \(L_2\), atoms that remain in the trap do not necessarily satisfy the no-jump condition. In fact, the condition requires the complete absence of spontaneous emission from \(|3\rangle\) to \(|2\rangle\), which limits the validity of the non-Hermitian Hamiltonian to transient dynamics. And, similar to the case with the jump operator \(L^\prime _1\), the timescale of this transient dynamics decreases exponentially with increasing atom number.

Nevertheless, the exceptional points do have consequences in the open-system dynamics. This can be revealed by analyzing the Liouvillian eigenspectrum and studying the LEP encircling. Specifically, we write the Liouvillian superoperator in a non-Hermitian matrix form

\(\begin{aligned} \mathcal {L}= \left[ \begin{array}{cccc} 0 & iJ & -iJ & 0 \\ iJ & -i\delta -\frac{\Gamma +\gamma }{2} & 0 & -iJ \\ -iJ & 0 & i\delta -\frac{\Gamma +\gamma }{2} & iJ\\ 0 & -iJ & iJ & -\Gamma \\ \end{array}\right] , \end{aligned}\)
(18)

where we enforce the no-jump condition for operator \(L_1\). Physically, this corresponds to detecting atoms that remain in the trap, in the presence of dephasing dictated by \(L_2\). Similar to the example in Sect. 2, the Liouvillian eigenspectrum \(\lambda\) is calculated through \(\mathcal {L}\rho _\lambda =\lambda \rho _\lambda\). But a steady-state solution with \(\lambda =0\) is no longer present, due to the post selection process which does not reserve the trace of the density matrix. We show the Liouvillian eigenspectrum in Fig. 5, where we identify the high-lying sheet as the quasi-steady-state solution. The quasi-steady state (with their \(\text {Re}\lambda <0\)) would become steady state (with \(\lambda =0\)) when the strength of \(L_1\) is continuously turned off. More importantly, an exceptional structure is identified in the eigenspectrum below the quasi-steady state, see for instance Fig. 5a. More concretely, with the addition of \(L_2\), the HEP of the original non-Hermitian Hamiltonian develops into an exceptional structure in the Liouvillian eigenspectrum, featuring exceptional lines ending at higher-order LEPs, with a structure similar to that in Fig. 2.

The encircling dynamics of similar LEPs was first experimentally studied in superconducting qubits [50, 51], where chiral behaviors were reported to persist for adiabatic encircling. However, a relevant question is, since there are no exceptional structures in the steady-state subspace of these systems, the long-time dynamics should still be dominated by the steady state, and not by the exceptional structure.

The question is resolved in Ref. [60], where chiral state transfer in the presence of LEP is studied in the context of the cold-atom experiment. First, since adiabaticity corresponds to the requirement that the rate of parameter change be much smaller than the instantaneous Liouvillian gap, the system would remain in the (quasi)-steady state of the Liouvillian. It follows that, since the LEP structure lies in the excited eigenstates of the Liouvillian, it has no bearing on the adiabatic dynamics. Indeed, as demonstrated in Ref. [60], in the long-time limit, the system always returns to the initial state, regardless of the encircling direction. This is illustrated in Fig. 5a, b.

By contrast, a chiral state transfer is observed in the intermediate regime, in the density-matrix evolution driven by the Liouvillian superoperator. This is illustrated in Fig. 5c, d. Qualitatively, such a chiral state transfer can be understood as the manifestation of the LEP encircling, which is exactly relevant at intermediate times. We note that, the observed LEP encircling in superconducting qubits also occur at intermediate time scales [50, 51, 60]. Similarly, a recent demonstration of the LEP-facilitated quantum heat engine is also based on the dynamic impact of LEP at intermediate timescales [58].

6 Steady-state LEP

Exceptional structures can also emerge in the steady-state subspace of a quantum open system. In this case, the EP-related chiral state transfer occurs in the long-time limit, even as the system adiabatically follows the steady state. The bifurcation of the steady-state solutions at the LEP corresponds to multi-stability, which is commonly observed in non-linear systems or quantum many-body systems under the mean-field approximation [65, 66].

Fig. 6
figure 6

a Rydberg level scheme. b Steady-state solutions of \(n_R=\rho _{22}\) from the optical Bloch equations under \(\Omega =2\). c, d Trajectories on the landscape of the steady-state Rydberg population. The green lines are the second-order Liouvillian exceptional lines and the yellow star is the third-order LEP. The red (blue) region indicates the stable (unstable) steady state. The encircling path is as follows: \(\Omega (t)=3.85+1.477\sin {(\pm 2\pi t/T+\phi _0)}\), \(\Delta (t)=-5.6+1.477\cos {(\pm 2\pi t/T+\phi _0)}\), and \(\phi _0=-\arctan {(9/4)}\). Here \(T=50000\), and other parameters are: \(\gamma =1\), and \((N-1)V=-11\)


For instance, consider the following optical Bloch equations

\(\begin{aligned} \dot{\rho }_{22} & =-\Omega {\text {Im}} \rho _{21}-\gamma \rho _{22}, \end{aligned}\)
(19)

\(\begin{aligned} \dot{\rho }_{2 1} & =i\left[ \Delta -(N-1)V \rho _{22}\right] \rho _{21}-\frac{\gamma }{2} \rho _{21}+i \Omega \left( \rho _{22}-\frac{1}{2}\right) , \end{aligned}\)
(20)

which describe the dynamics of a driven-dissipative Rydberg gas of N atoms. As illustrated in Fig. 6a, the states \(|1\rangle\) and \(|2\rangle\) correspond to the ground and Rydberg states, respectively. The density-dependent detuning \(\Delta -(N-1)V \rho _{22}\) originates from the Rydberg interactions under the mean-field approximation [67,68,69,70,71], while \(\Omega\) and \(\Delta\) are the Rabi frequency and detuning of the Rydberg coupling laser.

The steady-state solutions of the optical Bloch equations are obtained by setting \(\dot{\rho }=0\) and solving the resulting algebraic equations. With the non-linear detuning, the system admits either one or three steady-state solutions, depending on the parameters. As shown in Fig. 6b, in the regime with three solutions, there are an unstable solution and two stable ones. The two stable states are both many-body steady states featuring high and low Rydberg excitations, respectively. This is the well-known bistability which gives rise to the hysteresis in the light transmission of Rydberg gases in the electromagnetically induced transparency (EIT) measurements.

Crucially, new insights can be obtained from the corresponding Liouvillian eigenspectrum of the optical Bloch equations. As shown in Fig. 6c, d, the boundary of the bistable region marks the coalescence of the unstable and stable steady-state solutions. They can hence be identified as the second-order Liouvillian exceptional lines, which terminate at a third-order exceptional point where all three steady-state solutions merge. Here, only the eigenstates of the Liouvillian superoperator coalesce, while their eigenvalues remain zero due to the steady-state nature. This is different from the EPs (HEP or LEP) in linear systems where both eigenstates and eigenvalues coalesce. Furthermore, since the exceptional structure in the steady states arises from non-linearity, they do not reside on the Riemann surfaces, as in the case of conventional EPs. Nevertheless, the overall topology of the exceptional structure is similar, leading to similar dynamics features as discussed below. In the literature, such exceptional structures consisting of EP lines and higher-order EPs are also called exceptional nexus [72, 73].

Unlike the previously discussed LEPs in the excited states, the exceptional structure in the steady-state subspace has significant impact on the long-time dynamics. An outstanding example is the chiral state transfer, as recently demonstrated in a thermal Rydberg gas [65]. When the parameters are slowly modulated around the exceptional structure in a closed loop, the final state either switches or remains unchanged, depending on the modulation direction. As illustrated in Fig. 6c, d, along one direction, the system adiabatically follows the steady state, whereas along the other direction, a jump is unavoidable as the trajectory crosses the exceptional line. Based on numerical simulations, a set of sufficient conditions for the chiral state transfer is summarized [65]: (i) the initial stable steady state should be prepared in the bistable region; (ii) the two nearest (with respect to the initial position) intersection points between the trajectory and the exceptional lines should be on either side of the third-order LEP. Under these conditions, the chiral state transfer can still be observed, even as the trajectory crosses the exceptional lines multiple times. Furthermore, since the chiral state transfer only occurs when the trajectories are traversed sufficiently slowly, the exceptional landscape in the steady state becomes important. The chirality disappears for sufficiently fast parameter changes [65]. An interesting feature of the system is that both the exceptional structure and the chiral state transfer are subject to the tuning of many-body parameters. In thermal Rydberg atoms, these parameters include additional microwave fields that couple different Rydberg states, or the temperature which affects the density of the thermal gas. These possibilities pave the way for interesting schemes of quantum control.

Note that similar exceptional structures have been reported in non-linear non-Hermitian systems [72, 74, 75]. The exceptional structure therein derives from the classical non-linearity, with the simultaneous coalescence of eigenstates and eigenvalues. It is therefore fundamentally different from the Liouvillian exceptional structure in the steady-state manifold. Nevertheless, chiral state transfer can also be observed based on the non-linear HEPs [66].

7 Summary

Exceptional structures and the corresponding exceptional dynamics occur in a wide range of settings. While they have attracted significant interest in classical non-Hermitian models, recent studies have revealed their relevance and impact in quantum open systems, where exceptional structures emerge in the Liouvillian eigenspectrum and affect the steady-state approaching dynamics. Our brief review focuses on the chiral state transfer, a particular dynamic consequence, near the Liouvillian EPs. Based on these studies, it would be interesting to explore other exceptional features such as the EP-enhanced sensitivity and criticality near the Liouvillian EPs in quantum open settings. This is particularly intriguing in physical platforms such the driven-dissipative Rydberg gases, where the development of novel quantum control and sensing schemes would also benefit the on-going explorations therein for quantum information and computation. It is also interesting to study the generation and dynamic consequences of higher-order LEPs, wherein the complicated exceptional landscape can lead to richer encircling possibilities but have rarely been explored in quantum open systems. Finally, given the recent interest in the non-Abelian braid topology in non-Hermitian multi-band systems [76, 77], it is tempting to investigate and exploit similar features in quantum open systems, based on our understanding of the LEPs.

Availability of data and materials

Codes and numerical results of this review are available from the corresponding author upon reasonable request.

References

  1. Y. Ashida, Z. Gong, M. Ueda, Non-Hermitian Physics, Adv. Phys. 69, 249 (2020)

  2. T. Kato, Perturbation Theory of Linear Operators (Springer, 1966)

  3. M.V. Berry, D.H.J. O’Dell, Diffraction by volume gratings with imaginary potentials, J. Phys. Math. Gen. 31, 2093 (1998)

  4. W.D. Heiss, Phases of wave functions and level repulsion, Eur. Phys. J. D 7, 1 (1999)

  5. C.M. Bender, S. Boettcher, Real spectra in non-Hermitian Hamiltonians having PT symmetry, Phys. Rev. Lett. 80, 5243 (1998)

  6. C.M. Bender, Making sense of non-Hermitian Hamiltonians, Rep. Prog. Phys. 70, 947 (2007)

  7. V.K. Konotop, J. Yang, D.A. Zezyulin, Nonlinear waves in 𝒫𝒯-symmetric systems, Rev. Mod. Phys. 88, 035002 (2016)

  8. L. Feng, R. El-Ganainy, L. Ge, Non-Hermitian photonics based on parity–time symmetry, Nat. Photonics 11, 752 (2017)

  9. R. El-Ganainy, K.G. Makris, M. Kharjavikhan, Z.H. Mussilimani, S. Rotter, D.N. Christodoulides, Non-Hermitian physics and PT symmetry, Nat. Phys. 14, 11 (2018)

  10. M.-A. Miri, A. Alù, Exceptional points in optics and photonics, Science 363, eaar7709 (2019)

  11. ŞK. Özdemir, S. Rotter, F. Nori, L. Yang, Parity–time symmetry and exceptional points in photonics, Nat. Mat. 18, 783 (2019)

  12. K. Ding, C. Fang, G.-C. Ma, Non-Hermitian topology and exceptional-point geometries, Nat. Rev. Phys. 4, 745 (2022)

  13. A.-D. Li, H. Wei, M. Cotrufo, W.J. Chen, S. Mann, X. Ni, B.-C. Xu, J.-F. Chen, J. Wang, S.-H. Fan, C.-W. Qiu, A. Alù, L. Chen, Exceptional points and non-Hermitian photonics at the nanoscale, Nat. Nanotech. 18, 706 (2023)

  14. K.G. Makris, R. El-Ganainy, D.N. Christodoulides, Z.H. Musslimani, Beam Dynamics in 𝒫⁢𝒯 Symmetric Optical Lattices, Phys. Rev. Lett. 100, 103904 (2008)

  15. S. Klaiman, U. Günther, N. Moiseyev, Visualization of Branch Points in 𝒫⁢𝒯-Symmetric Waveguides, Phys. Rev. Lett. 101, 080402 (2008)

  16. C.E. Rüter, K.G. Makris, R. El-Ganainy, D.N. Christodoulides, M. Segev, D. Kip, Observation of parity–time symmetry in optics, Nat. Phys. 6, 192 (2010)

  17. B. Peng, ŞK. Özdemir, F.-C. Lei, F. Monifi, M. Gianfreda, G.-L. Long, S.-H. Fan, F. Nori, C.M. Bender, L. Yang, Parity–time-symmetric whispering-gallery microcavities, Nat. Phys. 10, 394 (2014)

  18. L. Chang, X. Jiang, S. Hua, C. Yang, J. Wen, L. Jiang, G. Wang, M. Xiao, Parity–time symmetry and variable optical isolation in active–passive-coupled microresonators, Nat. Photon. 8, 524 (2014)

  19. M. Brandstetter, M. Liertzer, C. Deutsch, P. Klang, J. Schöberl, H.E. Tüerci, G. Strasser, K. Unterrainer, S. Rotter, Reversing the pump dependence of a laser at an exceptional point, Nat. Commun. 5, 4034 (2014)

  20. B. Peng, ŞK. Özdemir, S. Rotter, H. Yilmaz, M. Liertzer, F. Monifi, C.M. Bender, F. Nori, L. Yang, Loss-induced suppression and revival of lasing, Science 346, 328 (2014)

  21. H. Hodaei, M.-A. Miri, M. Heinrich, D. N. Christodoulides, M. Khajavikhan, Loss-induced suppression and revival of lasing, Science 346, 975 (2014)

  22. L. Feng, Z.J. Wong, R.-M. Ma, Y. Wang, X. Zhang, Single-mode laser by parity-time symmetry breaking, Science 346, 972 (2014)

  23. K. Kawabata, Y. Ashida, M. Ueda, Information retrieval and criticality in parity-time-symmetric systems, Phys. Rev. Lett. 119, 190401 (2017)

  24. B. Dóra, M. Heyl, R. Moessner, The Kibble-Zurek mechanism at exceptional points, Nat. Commun. 10, 2254 (2019)

  25. L. Xiao, K. Wang, X. Zhan, Z. Bian, K. Kawabata, M. Ueda, W. Yi, P. Xue, Observation of critical phenomena in parity-time-symmetric quantum dynamics, Phys. Rev. Lett. 123, 230401 (2019)

  26. L. Xiao, D. Qu, K. Wang, H.-W. Li, J.-Y. Dai, B. Dóra, M. Heyl, R. Moessner, W. Yi, P. Xue, Non-Hermitian Kibble-Zurek mechanism with tunable complexity in single-photon interferometry, Phys. Rev. X Quantum 2, 020313 (2021)

  27. O. Latinne, N.J. Kylstra, M. Dörr, J. Purvis, M. Terau-Dunseath, C.J. Joachain, P.G. Burke, C.J. Noble, Laser-induced degeneracies involving autoionizing states in complex atoms, Phys. Rev. Lett. 74, 46 (1995)

  28. R. Lefebvre, O. Atabek, M. Šindelka, N. Moiseyev, Resonance coalescence in molecular photodissociation, Phys. Rev. Lett. 103, 123003 (2009)

  29. R. Uzdin, A. Mailybaev, N. Moiseyev, On the observability and asymmetry of adiabatic state flips generated by exceptional points, J. Phys. A: Math. Theor. 44, 435302 (2011)

  30. M.V. Berry, R. Uzdin, Slow non-Hermitian cycling: exact solutions and the Stokes phenomenon, J. Phys. A: Math. Theor. 44, 435303 (2011)

  31. T.J. Milburn, J. Doppler, C.A. Holmes, S. Portolan, S. Rotter, P. Rabl, General dynamical description of quasi-adiabatically encircling exceptional points, Phys. Rev. A 92, 052124 (2015)

  32. H. Nasari, G. Lopez-Galmiche, H.E. Lopez-Aviles, A. Schumer, A.U. Hassan, Q. Zhong, S. Rotter, P. LiKamWa, D.N. Christodoulides, M. Khajavikhan, Observation of chiral state transfer without encircling an exceptional point, Nature 605, 256–261 (2022)

  33. Z.-J. Ren, D. Liu, E.-T. Zhao, C.-D. He, K.K. Pak, J. Li, G.-B. Jo, Chiral control of quantum states in non-Hermitian spin–orbit-coupled fermions, Nat. Phys. 18, 385 (2022)

  34. J. Doppler, A.A. Mailybaev, J. Böhm, U. Kuhl, A. Girschik, F. Libisch, T.J. Milburn, P. Rabl, N. Moiseyev, S. Rotter, Dynamically encircling an exceptional point for asymmetric mode switching, Nature 537, 76 (2016)

  35. X.-L. Zhang, S.-B. Wang, B. Hou, C.T. Chan, Dynamically Encircling Exceptional Points: In situ Control of Encircling Loops and the Role of the Starting Point, Phys. Rev. X 8, 021066 (2018)

  36. Q.-H. Song, M. Odeh, J. Zúñiga-Pérez, B. Kanté, P. Genevet, Plasmonic topological metasurface by encircling an exceptional point, Science 373, 6559 (2021)

  37. L.-Y. Ding, K.-Y. Shi, Y.-X. Wang, Q.-X. Zhang, C.-H. Zhu, L.-D. Zhang, J.-Q. Yi, S.-N. Zhang, X. Zhang, W. Zhang, Information retrieval and eigenstate coalescence in a non-Hermitian quantum system with anti-𝒫𝒯 symmetry, Phys. Rev. A 105, L010204 (2022)

  38. W.-Q. Liu, Y. Wu, C.-K. Duan, X. Rong, J.-F. Du, Dynamically encircling an exceptional point in a real quantum system, Phys. Rev. Lett. 126, 170506 (2021)

  39. Z.-Z. Li, W. Chen, M. Abbasi, K.W. Murch, K.B. Whaley, Speeding up entanglement generation by proximity to higher-order exceptional points, Phys. Rev. Lett. 131, 100202 (2023)

  40. P.-R. Han, F. Wu, X.-J. Huang, H.-Z. Wu, C.-L. Zou, W. Yi, M. Zhang, H. Li, K. Xu, D. Zheng, H. Fan, J. Wen, Z.-B. Yang, S.-B. Zheng, Exceptional entanglement phenomena: Non-Hermiticity meeting nonclassicality, Phys. Rev. Lett. 131, 260201 (2023)

  41. W.-J. Chen, ŞK. Özdemir, G.-M. Zhao, J. Wiersig, L. Yang, Exceptional points enhance sensing in an optical microcavity, Nature 548, 192 (2017)

  42. Y.-H. Lai, Y.-K. Lu, M.-G. Suh, Z.-Q. Yuan, K. Vahala, Observation of the exceptional-point-enhanced Sagnac effect, Nature 576, 65 (2019)

  43. M.-Z. Zhang, W. Sweeney, C.W. Hsu, L. Yang, A.D. Stone, L. Jiang, Quantum noise theory of exceptional point amplifying sensors, Phys. Rev. Lett. 123, 180501 (2019)

  44. H. Xu, D. Mason, L. Jiang, J.G.E. Harris, Topological energy transfer in an optomechanical system with exceptional points, Nature 537, 80 (2016)

  45. H. Jing, S. Özdemir, X.-Y. Lü, J. Zhang, L. Yang, F. Nori, 𝒫𝒯-Symmetric Phonon Laser, Phys. Rev. Lett. 113, 053604 (2014)

  46. X. Zhu, H. Ramezani, C. Shi, J. Zhu, X. Zhang, 𝒫⁢𝒯-Symmetric Acoustics, Phys. Rev. X 4, 031042 (2014)

  47. R. Fleury, D. Sounas, A. Alù, An invisible acoustic sensor based on parity-time symmetry, Nat. Commun. 6, 5905 (2015)

  48. K. Ding, G. Ma, M. Xiao, Z.Q. Zhang, C.T. Chan, Emergence, coalescence, and topological properties of multiple exceptional points and their experimental realization, Phys. Rev. X 6, 021007 (2016)

  49. J. Li, A.K. Harter, J. Liu, L. de Melo, Y.N. Joglekar, L. Luo, Observation of parity-time symmetry breaking transitions in a dissipative Floquet system of ultracold atoms, Nat. Commun. 10, 855 (2019)

  50. W.-J. Chen, M. Abbasi, Yogesh N. Joglekar, K. W. Murch, Quantum jumps in the non-Hermitian dynamics of a superconducting qubit, Phys. Rev. Lett. 127, 140504 (2021)

  51. W.-J. Chen, M. Abbasi, B. Ha, S. Erdamar, Y.N. Joglekar, K.W. Murch, Decoherence-induced exceptional points in a dissipative superconducting qubit, Phys. Rev. Lett. 128, 110402 (2022)

  52. H. P. Breuer, F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2002)

  53. M.B. Plenio, P.L. Knight, The quantum-jump approach to dissipative dynamics in quantum optics, Rev. Mod. Phys. 70, 101 (1998)

  54. P. Kumar, K. Snizhko, Y. Gefen, Near-unit efficiency of chiral state conversion via hybrid-Liouvillian dynamics, Phys. Rev. A 104, L050405 (2021)

  55. M. Am-Shallem, R. Kosloff, N. Moiseyev, Exceptional points for parameter estimation in open quantum systems: Analysis of the Bloch equations, New J. Phys. 17, 113036 (2015)

  56. F. Minganti, A. Miranowicz, R.W. Chhajlany, F. Nori, Quantum exceptional points of non-Hermitian Hamiltonians and Liouvillians: The effects of quantum jumps, Phys. Rev. A 100, 062131 (2019)

  57. F. Minganti, A. Miranowicz, R.W. Chhajlany, I.I. Arkhipov, F. Nori, Hybrid-Liouvillian formalism connecting exceptional points of non-Hermitian Hamiltonians and Liouvillians via postselection of quantum trajectories, Phys. Rev. A 101, 062112 (2020)

  58. J.-W. Zhang, J.-Q. Zhang, G.-Y. Ding, J.-C. Li, J.-T. Bu, B. Wang, L.-L. Yan, S.-L. Su, L. Chen, F. Nori, ŞK. Özdemir, F. Zhou, H. Jing, M. Feng, Dynamical control of quantum heat engines using exceptional points, Nat. Commun. 13, 6225 (2022)

  59. J.-T. Bu, J.-Q. Zhang, G.-Y. Ding, J.-C. Li, J.-W. Zhang, B. Wang, W.-Q. Ding, W.-F. Yuan, L. Chen, ŞK. Özdemir, F. Zhou, H. Jing, M. Feng, Enhancement of quantum heat engine by encircling a Liouvillian exceptional point, Phys. Rev. Lett. 130, 110402 (2023)

  60. K. Sun, W. Yi, Chiral state transfer under dephasing, Phys. Rev. A 108, 013302 (2023)

  61. F. Song, S. Yao, Z. Wang, Non-Hermitian skin effect and chiral damping in open quantum systems, Phys. Rev. Lett. 123, 170401 (2019)

  62. Y.-N. Zhou, L. Mao, H. Zhai, Rényi entropy dynamics and Lindblad spectrum for open quantum systems, Phys. Rev. Research 3, 043060 (2021)

  63. T. Li, Y.-S. Zhang, W. Yi, Engineering dissipative quasicrystals, Phys. Rev. B 105, 125111 (2022)

  64. H. Li, H. Wu, W. Zheng, W. Yi, Many-body non-Hermitian skin effect under dynamic gauge coupling, Phys. Rev. Research 5, 033173 (2023)

  65. C. Xie, K. Sun, K.-D. Wu, C.-F. Li, G.-C. Guo, W. Yi, G.-Y. Xiang, Chiral switching of many-body steady states in a dissipative Rydberg gas, arXiv:2402.02779

  66. Z. Li, C. Li, G. Xu, W. Chen, Z. Xiong, H. Jing, J.-S. Ho, C.-W. Qiu, Synergetic positivity of loss and noise in nonlinear non-Hermitian resonators, Sci. Adv. 9, 27 (2023)

  67. C. Carr, R. Ritter, C.G. Wade, C.S. Adams, K.J. Weatherill, Nonequilibrium phase transition in a dilute Rydberg ensemble, Phys. Rev. Lett. 111, 113901 (2013)

  68. N. Šibalić, C.G. Wade, C.S. Adams, K.J. Weatherill, T. Pohl, Driven-dissipative many-body systems with mixed power-law interactions: Bistabilities and temperature-driven nonequilibrium phase transitions, Phys. Rev. A 94, 011401(R) (2016)

  69. N.R. de Melo, C.G. Wade, N. Šibalić, J.M. Kondo, C.S. Adams, K.J. Weatherill, Intrinsic optical bistability in a strongly driven Rydberg ensemble, Phys. Rev. A 93, 063863 (2016)

  70. D.-S. Ding, H. Busche, B.-S. Shi, G.-C. Guo, C.S. Adams, Phase diagram and self-organizing dynamics in a thermal ensemble of strongly interacting Rydberg atoms, Phys. Rev. X 10, 021023 (2020)

  71. K. Wadenpfuhl, C.S. Adams, Emergence of synchronization in a driven-dissipative hot Rydberg vapor, Phys. Rev. Lett. 131, 143002 (2023)

  72. K. Bai, L. Fang, T.-R. Liu, J.-Z. Li, D.-D. Wan, M. Xiao, Nonlinearity-enabled higher-order exceptional singularities with ultra-enhanced signal-to-noise ratio, Natl. Sci. Rev. 10, 7 (2022)

  73. X. Zhou, X.-J. Ren, D.-B. Xiao, J.-Q. Zhang, R. Huang, Z.-P. Li, X.-P. Sun, X.-Z. Wu, C.-W. Qiu, F. Nori, H. Jing, Higher-order singularities in phase-tracked electromechanical oscillators, Nat. Commun. 14, 7944 (2023)

  74. S. Assawaworrarit, X.-F. Yu, S.-H. Fan, Robust wireless power transfer using a nonlinear parity–time-symmetric circuit, Nature 546, 387 (2017)

  75. H.-W. Wang, S. Assawaworrarit, S.-H. Fan, Dynamics for encircling an exceptional point in a nonlinear non-Hermitian system, Opt. Lett. 44, 3 (2019)

  76. K. Wang, A. Dutt, C.C. Wojcik, S. Fan, Topological complex-energy braiding of non-Hermitian bands, Nature 598, 59 (2021)

  77. Y.S.S. Patil, J. Höller, P.A. Henry, C. Guria, Y. Zhang, L. Jiang, N. Kralj, N. Read, J.G.E. Harris, Measuring the knot of non-Hermitian degeneracies and non-commuting braids, Nature 607, 271 (2022)

Acknowledgements

See the “Funding” section.

Funding

This research is supported by the National Natural Science Foundation of China (Grant No. 12374479).

Author information

Authors and Affiliations

Contributions

K. S. performed theoretical derivations and numerical calculations. W. Y. supervised the project. Both authors contributed to the writing of the manuscript.

Corresponding author

Correspondence to Wei Yi.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

[Source: https://link.springer.com/article/10.1007/s43673-024-00129-3]