Liouvillian Exceptional Points in Open Quantum Systems

• Liouvillian exceptional points are non-diagonalizable spectral degeneracies where eigenvalues and eigenmatrices coalesce into high-dimensional Jordan blocks.
• They control transitions between dynamical regimes, leading to anomalous decay behaviors and enhanced sensitivity in quantum measurements.
• Experimental approaches like quantum process tomography and spectroscopic analysis reveal LEPs, advancing quantum sensing and dissipative control strategies.

A Liouvillian exceptional point (LEP) is a non-Hermitian non-diagonalizability in the spectrum of the Liouvillian superoperator that governs Lindbladian open-system quantum dynamics. Unlike Hamiltonian exceptional points (HEPs), which pertain to non-Hermitian Hamiltonians neglecting quantum jumps, LEPs intrinsically describe the full dissipative and measurement-including evolution—quantum-jump processes, non-reciprocal relaxation, and full density-matrix dynamics. LEPs control the transition between dynamical regimes (e.g., underdamped/overdamped, oscillatory/monotonic decay), underlie critical sensitivity features in measurement and response, and enable robust construction of non-Hermitian quantum protocols impossible within the pure-state Hamiltonian paradigm. LEPs are characterized mathematically by a coalescence of multiple eigenvalues and associated right and left eigenmatrices of the Liouvillian, leading to higher-dimensional Jordan blocks, altered power-law splittings of spectral lines, anomalous time evolution, and topological phenomena under parameter variation. Recent research has established LEPs as a ubiquitous and experimentally accessible feature of both few- and many-body open quantum systems, with rigorous connections to system control, quantum sensing, and nonequilibrium phase transitions.

1. Formal Definition and Spectral Structure

The Lindblad master equation describes the evolution of a density operator $\rho$ for a generic Markovian open quantum system: $$\frac{d}{dt}\rho = \mathcal L[\rho] = -i[H,\rho] + \sum_k \gamma_k \left( L_k \rho L_k^\dagger - \frac{1}{2}\{ L_k^\dagger L_k, \rho \} \right)$$ where $H$ is the system Hamiltonian, $L_k$ are the quantum jump operators, and $\gamma_k$ their rates.

The Liouvillian $\mathcal L$ is viewed as a superoperator acting on the operator space (Liouville space), typically via vectorization $\rho \mapsto |\rho\rangle\rangle$. The right and left spectral problems are

$$\mathcal L |R_n\rangle = \lambda_n |R_n\rangle, \quad \langle L_n| \mathcal L = \lambda_n \langle L_n|$$

The spectrum $\{\lambda_n\}$ consists of a unique zero eigenvalue for the steady state, and a collection of generally complex eigenvalues for decaying and oscillating modes.

A Liouvillian exceptional point of order $m$ arises at parameter values where $m$ eigenvalues $\lambda_i$ coalesce, and the corresponding right and left eigenmatrices also merge, reducing the geometric multiplicity to one and yielding a non-trivial $m \times m$ Jordan block; that is, the Liouvillian is nondiagonalizable and only similar to a block-upper-triangular form. Explicitly, for a parameter $p$,

$$\mathcal L(p_{\mathrm{EP}}) \longrightarrow \text{Jordan block: } \begin{pmatrix} \lambda_{EP}&1&0&\dots&0\ 0&\lambda_{EP}&1&\dots&0\ \vdots &&\ddots&\ddots&\ 0&\ldots&0&\lambda_{EP} \end{pmatrix}$$

A necessary and sufficient algebraic condition is that the first $m-1$ derivatives of the characteristic polynomial vanish at the degenerate eigenvalue $\lambda_{\mathrm{EP}}$: $$P(\lambda_{\mathrm{EP}}) = 0, \quad P'(\lambda_{\mathrm{EP}}) = 0, \quad \dots, \quad P^{(m-1)}(\lambda_{\mathrm{EP}}) = 0, \quad P^{(m)}(\lambda_{\mathrm{EP}}) \neq 0$$ This non-diagonalizability underlies anomalous system dynamics and is the spectral signature distinguishing LEPs from diabolical points, where only eigenvalues—itself, not eigenmatrices—coincide (Abo et al., 2024, Wu et al., 1 Dec 2025).

2. Physical Manifestations and Dynamical Consequences

Crossing an LEP produces a qualitative transition in system dynamics. For second-order LEPs, the eigenvalues exhibit a square-root branch-point splitting in parameters, leading to a transition from purely real (overdamped/exponential) to complex-conjugate (underdamped/oscillatory or spiraling) decay rates (Abo et al., 2024, Zhou et al., 2023). For third- or higher-order LEPs, more complex (e.g., cubic-root, fractional) splittings appear (Wu et al., 1 Dec 2025, Xu et al., 31 Jan 2026).

At the LEP, the Liouvillian propagator acquires additional polynomial-in-time prefactors multiplying the exponential decay: $$e^{\lambda_{EP} t} (A_0 + A_1 t + \cdots + A_{m-1} t^{m-1})$$ This directly results from the structure of the Jordan block. Such behavior gives rise to critical damping in open quantum systems—e.g., maximal aperiodic (non-oscillatory) decay with the fastest possible relaxation rate (Zhou et al., 2023, Khandelwal et al., 2021).

In spectral observables, LEPs manifest as higher-order poles in the resolvent $(z-\mathcal L)^{-1}$, leading to super-Lorentzian line shapes in emission or response spectra, distinguishable from ordinary Lorentzian profiles by their slower power-law decay in frequency and the presence of flat-top or non-Lorentzian wings (Molina, 1 Feb 2026, Arkhipov et al., 2020).

LEPs also underlie chiral state transfer and non-reciprocal dynamical responses when the system is driven along a closed loop in parameter space around the exceptional point—this yields nontrivial topological exchange of dynamical eigenmodes, not realizable in conventional Hermitian systems (Sun et al., 2024, Zhang et al., 6 Dec 2025, Gao et al., 17 Jan 2025).

3. Distinction from Hamiltonian Exceptional Points and Hybrid Liouvillian Formalism

Hamiltonian exceptional points (HEPs) describe coalescences of eigenvalues and eigenvectors of effective non-Hermitian Hamiltonians, which govern conditional evolution in the absence of quantum jumps—semiclassical or post-selected scenarios. Liouvillian exceptional points, in contrast, are defined for the full, quantum-jump-inclusive Liouvillian (Minganti et al., 2020, Kopciuch et al., 3 Jun 2025).

The spectrum, location, and even order of EPs in the full Liouvillian can differ dramatically from their Hamiltonian counterparts. For example, systems can exhibit a LEP where no HEP exists or the order of coalescence is greater (or less) than in the Hamiltonian picture due to the enlarged operator space (Abo et al., 2024, Wu et al., 1 Dec 2025, Arkhipov et al., 2020). Explicitly, quantum jumps can "lift" a higher-order HEP into multiple lower-order LEPs, shift their positions in parameter space, or even destroy the correspondence altogether (Kopciuch et al., 3 Jun 2025).

The hybrid Liouvillian is a parametric family

$\mathcal L_\eta = \mathcal L' + \eta \mathcal J$

interpolating between the no-jump (pure Hamiltonian) limit ($\eta=0$) and full Lindblad (quantum jumps, $\eta=1$) evolution, useful for tracing EP evolution with detector efficiency or postselection. This formalism enables continuous control and connects measurement and decoherence back-action to the nature and location of exceptional points (Minganti et al., 2020).

4. Experimental Detection and Tomographic Reconstruction

LEPs can be observed directly by reconstructing the Liouvillian via quantum process tomography (QPT) (Abo et al., 2024). The superoperator is reconstructed experimentally by preparing a basis set of input states, evolving them with the quantum process for a short time, and measuring probabilities in various bases. A short-time propagator $S = 1 + dt\,\mathcal L$ is extracted, from which the full matrix $\mathcal L$ is determined and diagonalized to find the spectrum and Jordan block structure.

On platforms supporting only unitary operations—such as IBM Quantum devices—general completely positive maps are implemented by embedding the system in a larger Hilbert space, using ancilla qubits and the Choi-Jamiołkowski isomorphism. Controlled-unitary circuits simulate the effect of non-unitary evolution, with experimental details including error mitigation and dynamical decoupling for robustness (Abo et al., 2024).

Signatures such as square-root closing of relaxation gaps, maximal spectral broadening, and coalescence of eigenmatrices (measured through overlap $|\langle L_1|R_2\rangle|$ approaches $1$ near the LEP) provide unambiguous identification, distinguishing LEPs from diabolical points (degeneracies without eigenmode coalescence).

In emission and spectroscopy, higher-order spectral poles are extracted via line-shape fits to Lorentzian and super-Lorentzian forms, with the weight $r=|b|/(|a|+|b|)$ used as an experimental diagnostic of an LEP (Molina, 1 Feb 2026).

5. Higher-Order and Many-Body Liouvillian Exceptional Points

Liouvillian superoperators can exhibit higher-order exceptional points far beyond second order, especially in systems with large Hilbert space dimension (e.g., multi-level atoms, collective spins, bosonic and fermionic chains). In linear bosonic systems, any $n$-th order exceptional point of an effective non-Hermitian Hamiltonian always generates an infinite hierarchy of higher-order Liouvillian EPs—detectable at arbitrarily high orders in the hierarchy of operator moments (through $g^{(2)}$, $g^{(3)}$, and higher correlators) (Arkhipov et al., 2020, Xu et al., 31 Jan 2026).

The presence of higher-order LEPs in fermionic chains leads to a gapless Liouvillian spectrum and algebraic (power-law) relaxation toward steady state, with the order of the Jordan block scaling with system size. Perturbations (e.g., multi-body quantum-jump processes) split the giant Jordan block, opening a finite Liouvillian gap with universal fractional scaling laws, $\Delta \sim z^{1/(d+1)}$ for a $d+1$-st order EP (Xu et al., 31 Jan 2026).

The generalized eigenvectors at an $N$-th order LEP allow freedom in basis choice under triangular similarity transformations, but the long-time dynamics universally exhibits secular polynomial prefactors $t^k e^{\lambda t}$, with $k$ up to $N-1$ (Tay, 2023, Xu et al., 31 Jan 2026).

6. Topological and Chiral Phenomena

Parametric encircling of an LEP in multidimensional parameter space results in nontrivial topological effects: eigenmodes are non-reciprocally transferred after one or more loops—chiral state transfer—a dynamical effect without Hermitian analogue and diagnostic of non-Hermitian topology (Sun et al., 2024, Gao et al., 17 Jan 2025). In the open quantum context, such chiral dynamics persists within timescales set by Liouvillian gaps; long-time evolution erases chirality due to relaxation to the unique steady state (Gao et al., 17 Jan 2025).

In collective or many-body systems, encircling can implement nontrivial braiding of eigenmodes, and multiple windings or hybrid winding numbers may arise in the presence of multiple coalescences or non-Markovian dynamics (Zhang et al., 6 Dec 2025). These signatures are experimentally accessible via joint quantum-state tomography and spectroscopic monitoring, with topological invariants extracted from the winding of complex eigenvalues in parameter space.

7. Applications: Control, Sensing, and Quantum Thermodynamics

LEPs are a resource for enhancing quantum device performance:

• Quantum Sensing: Enhanced parameter sensitivity appears near higher-order LEPs, with responses scaling as $(\theta-\theta_{EP})^{1/m}$ for an order-$m$ EP, providing increased quantum Fisher information (Wu et al., 1 Dec 2025, Popkov et al., 12 Oct 2025).
• Dissipative State Preparation: LEPs yield maximal Liouvillian gaps, enabling rapid convergence to the steady state—key for dissipative engineering of pure or entangled quantum states (Zhou et al., 2023).
• Quantum Heat Engines: Operation at, or staggered across, LEPs enhances efficiency and power extraction in quantum Otto engines, allowing for optimal work cycles by leveraging the boundary between oscillatory and monotonic decay (Zhang et al., 2022).
• Stroboscopic Quantum Circuits: LEPs and their sensing properties arise in discrete-time (brickwork) CPTP circuits realized on modern quantum computers, bridging the gap between continuous and digital open-system dynamics (Popkov et al., 12 Oct 2025).

8. Experimental Realizations and Outlook

LEPs have been observed in a wide range of platforms, including alkali vapor magnetometers (Kopciuch et al., 3 Jun 2025), superconducting transmon qubits (Chen et al., 2021), trapped ions (Wu et al., 1 Dec 2025, Zhang et al., 2022), collective atomic ensembles, and photonic and spin-based quantum simulators (Molina, 1 Feb 2026, Abo et al., 2024, Arkhipov et al., 2020). Quantum process tomography, time-resolved state tomography, dynamical encircling, and noise spectroscopy are principal methods for their identification in experiment.

Prospects for future work include non-Markovian generalizations (where the LEP may relate to effective Liouvillian descriptions after mapping to pseudomodes or extended system-bath Hilbert spaces) (Zhang et al., 6 Dec 2025), higher-order and many-body LEPs in large-scale quantum devices (Xu et al., 31 Jan 2026), and systematic design of LEPs with desired order and topology via algebraic and tropical-geometric protocols (P et al., 9 Oct 2025). These developments are expected to impact precision metrology, dissipative quantum information processing, and the engineering of nonequilibrium phases.

---

References:

• Experimental Liouvillian exceptional points in a quantum system without Hamiltonian singularities (Abo et al., 2024)
• Hybrid-Liouvillian formalism connecting exceptional points of non-Hermitian Hamiltonians and Liouvillians via postselection of quantum trajectories (Minganti et al., 2020)
• Accelerating relaxation through Liouvillian exceptional point (Zhou et al., 2023)
• Spectroscopic Signatures of a Liouvillian Exceptional Spectral Phase in a Collective Spin (Molina, 1 Feb 2026)
• Experimental Witness of Quantum Jump Induced High-Order Liouvillian Exceptional Points (Wu et al., 1 Dec 2025)
• Encircling the Liouvillian exceptional points: a brief review (Sun et al., 2024)
• Photonic chiral state transfer near the Liouvillian exceptional point (Gao et al., 17 Jan 2025)
• Liouvillian exceptional points of any order in dissipative linear bosonic systems (Arkhipov et al., 2020)
• Characterizing Liouvillian Exceptional Points Through Newton Polygons and Tropical Geometry (P et al., 9 Oct 2025)
• Signatures of Liouvillian exceptional points in a quantum thermal machine (Khandelwal et al., 2021)
• Decoherence Induced Exceptional Points in a Dissipative Superconducting Qubit (Chen et al., 2021)
• Lindblad dynamics of open multi-mode bosonic systems: Algebra of bilinear superoperators, exceptional points and speed of evolution (Gaidash et al., 2024)
• Dynamical Control of Quantum Heat Engines Using Exceptional Points (Zhang et al., 2022)
• Exploring the topology induced by non-Markovian Liouvillian exceptional points (Zhang et al., 6 Dec 2025)
• Emergent Liouvillian exceptional points from exact principles (Khandelwal et al., 2024)
• Liouvillian and Hamiltonian exceptional points of atomic vapors: The spectral signatures of quantum jumps (Kopciuch et al., 3 Jun 2025)
• Liouvillian exceptional points in continuous variable system (Tay, 2023)
• Exactly solvable higher-order Liouvillian exceptional points in dissipative fermionic systems (Xu et al., 31 Jan 2026)
• Liouvillian Exceptional Points in Quantum Brickwork Circuits (Popkov et al., 12 Oct 2025)
• Liouvillian exceptional points of an open driven two-level system (Seshadri et al., 2024)