Non-Hermitian degeneracies of Lindblad generators (Liouvillian exceptional points) can induce non-exponential relaxation and higher-order poles in dynamical response functions. A collective spin coupled to a polarized Markovian bath exhibits an \emph{exceptional spectral phase} in which defective Liouvillian modes imprint super-Lorentzian features in frequency-resolved spectra. We compute the emission spectrum via the Liouvillian resolvent, identify symmetry-sector selection rules, and demonstrate that exceptional-point signatures are strongly state-dependent: they are suppressed in steady-state fluorescence yet become unambiguous for generic (infinite-temperature or random) initial states. Our results provide an experimentally accessible spectroscopic diagnostic of many-body Liouvillian exceptional phases and clarify when steady-state emission can (and cannot) reveal them.
Non-Hermitian degeneracies of Lindblad generators (Liouvillian exceptional points) can induce non-exponential relaxation and higher-order poles in dynamical response functions. A collective spin coupled to a polarized Markovian bath exhibits an \emph{exceptional spectral phase} in which defective Liouvillian modes imprint super-Lorentzian features in frequency-resolved spectra. We compute the emission spectrum via the Liouvillian resolvent, identify symmetry-sector selection rules, and demonstrate that exceptional-point signatures are strongly state-dependent: they are suppressed in steady-state fluorescence yet become unambiguous for generic (infinite-temperature or random) initial states. Our results provide an experimentally accessible spectroscopic diagnostic of many-body Liouvillian exceptional phases and clarify when steady-state emission can (and cannot) reveal them.
Usually, when investigating exceptional points (EPs) of an open Markovian bosonic system, one deals with spectral degeneracies of a non-Hermitian Hamiltonian (NHH), which can correctly describe the system dynamics only in the semiclassical regime. A recently proposed quantum Liouvillian framework [Minganti et al., Phys. Rev. A 100, 062131 (2019)] enables the complete determination of the dynamical properties of such systems and their EPs (referred to as Liouvillian EPs, or LEPs) in the quantum regime by taking into account the effects of quantum jumps, which are ignored in the NHH formalism. Moreover, the symmetry and eigenfrequency spectrum of the NHH become a part of much larger Liouvillian eigenspace. As such, the EPs of an NHH form a subspace of the LEPs. Here we show that once an NHH of a dissipative linear bosonic system exhibits an EP of a certain finite order n, it immediately implies that the corresponding LEP can become of any higher order $m\ensuremath{\ge}n$ defined in the infinite Hilbert space. These higher-order LEPs can be identified by the coherence and spectral functions at the steady state. The coherence functions can offer a convenient tool to probe extreme system sensitivity to external perturbations in the vicinity of higher-order LEPs. As an example, we study a linear bosonic system of a bimodal cavity with incoherent mode coupling to reveal its higher-order LEPs; particularly, of second and third order via first- and second-order coherence functions, respectively. Accordingly, these LEPs can be additionally revealed by squared and cubic Lorentzian spectral lineshapes in the power and intensity-fluctuation spectra. Moreover, we demonstrate that these EPs can also be associated with spontaneous parity-time (PT) and anti-PT-symmetry breaking in the system studied. These symmetries can be switched in the output fields (the so-called supermodes) of an additional linear coupler with a properly chosen coupling strength. Thus, we show that the initial loss-loss dynamics for the supermodes can be equivalent to the balanced gain-loss evolution.
We study a model of a quantum collective spin weakly coupled to a spin-polarized Markovian environment and find that the spectrum is divided into two regions that we name normal and exceptional Liouvillian spectral phases. In the thermodynamic limit, the exceptional spectral phase displays the unique property of being made up exclusively of second-order exceptional points. As a consequence, the evolution of any initial density matrix populating this region is slowed down and cannot be described by a linear combination of exponential decays. This phase is separated from the normal one by a critical line in which the density of Liouvillian eigenvalues diverges, a phenomenon analogous to that of excited-state quantum phase transitions observed in some closed quantum systems. In the limit of no bath polarization, this criticality is transferred onto the steady state, implying a dissipative quantum phase transition and the formation of a boundary time crystal.
Spectral degeneracies in Liouvillian generators of dissipative dynamics generically occur as exceptional points, where the corresponding non-Hermitian operator becomes nondiagonalizable. Steady states, i.e., zero modes of Liouvillians, are considered a fundamental exception to this rule since a no-go theorem excludes nondiagonalizable degeneracies there. Here, we demonstrate that the crucial issue of diverging timescales in dissipative state preparation is largely tantamount to an asymptotic approach toward the forbidden scenario of an exceptional steady state in the thermodynamic limit. With case studies ranging from NP-complete satisfiability problems encoded in a quantum master equation to the dissipative preparation of a symmetry protected topological phase, we reveal the close relation between the computational complexity of the problem at hand, and the finite size scaling toward the exceptional steady state, exemplifying both exponential and polynomial scaling. Formally treating the weight W of quantum jumps in the Lindblad master equation as a parameter, we show that exceptional steady states at the physical value W=1 may be understood as a critical point hallmarking the onset of dynamical instability.
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.
Non-invertible symmetries extend conventional symmetries beyond group-like transformations and play an important role in quantum field theory and many-body physics. However, a general framework of non-invertible symmetries in open quantum systems is still lacking. Here, we formulate a general Liouville-space duality framework that extends Hilbert-space non-invertible dualities to open quantum systems. The resulting defect relates Liouvillians with dual Hamiltonian and dissipative parameters and becomes an exact symmetry at self-duality. We apply this framework to the periodic transverse-field Ising chain with two collective Lindblad channels exchanged by Kramers--Wannier (KW) duality. A strong Z2 symmetry divides Liouville space into four parity sectors, while the defect acts nontrivially only in the even-even sector. At self-duality, the KW defect imposes a momentum-resolved selection rule on Liouvillian Jordan blocks. Modes with different Liouville momenta or KW-defect eigenvalues cannot form a common exceptional Jordan block. For a three-site chain, a global Rényi-2 KW asymmetry dynamically identifies the self-dual conditions. At the same time, the spectral structure exhibits both a semisimple degeneracy and a self-dual exceptional point with a perturbation-selective response. These results extend non-invertible dualities to Markovian quantum dynamics and establish a symmetry-based selection rule for Liouvillian Jordan structures.
Exceptional points (EPs), indicative of parity-time (PT) symmetry breaking, play a central role in non-Hermitian physics, yet most studies begin from deliberately engineered effective Hamiltonians whose parameters are tuned to exhibit exceptional behavior. In realistic open quantum systems, however, dynamics are governed by Lindblad superoperators whose spectral structure is high-dimensional, symmetry-constrained, and not obviously reducible to minimal non-Hermitian models. A general framework for discovering exceptional dynamics directly from microscopic dissipative models has been lacking. Here we introduce a symmetry-resolved approach for identifying and characterizing exceptional points directly from the full Liouvillian generator. Correlated dissipation induces graph symmetries that decompose Liouville space into low-dimensional invariant sectors, within which minimal non-Hermitian blocks govern the onset of EPs and PT-breaking behavior. We further introduce a numerical diagnostic - the exceptional-point strength E - based on eigenvector conditioning, which quantifies proximity to defective dynamics without requiring analytic reduction. Applied to tight-binding models with correlated dephasing and relaxation, the method reproduces analytically predicted exceptional seams and reveals universal scaling of E near EP manifolds. More broadly, the framework enables systematic discovery of hidden exceptional structure in complex or high-dimensional open systems and is naturally compatible with matrix-free and tensor-network implementations for scalable many-body applications.