---
title: Liouvillian Exceptional Point
url: https://www.emergentmind.com/topics/liouvillian-exceptional-point
type: topic
---

# Liouvillian Exceptional Point

Searching arXiv for recent papers on Liouvillian exceptional points to ground the article in published work.
arXiv search query: "Liouvillian exceptional point"
A Liouvillian exceptional point is a non-diagonalizable spectral degeneracy of the Liouvillian superoperator that generates open-system dynamics. In operator form, the density matrix obeys a master equation $\dot\rho=\mathcal L[\rho]$, and an LEP occurs when two or more Liouvillian eigenvalues coalesce simultaneously with their corresponding eigenoperators, so that the generator develops a Jordan block. Because $\mathcal L$ includes both coherent evolution and quantum-jump recycling terms, LEPs are generally inequivalent to exceptional points of an effective non-Hermitian Hamiltonian, and can differ in existence, location, order, and topology. Recent work places LEPs at the intersection of non-Hermitian spectral theory, dissipative quantum dynamics, open-system thermodynamics, and topological mode conversion, with extensions from standard Lindblad generators to non-Markovian pseudomode constructions and even discrete completely positive trace-preserving circuits [2002.11620] [2512.06311] [2510.10629].

## 1. Formal definition and algebraic structure

For a time-independent Markovian open quantum system, the Liouvillian is defined through the Lindblad equation
\[
\frac{d\rho}{dt}
=
\mathcal L[\rho]
=
-\,i[H,\rho]
+
\sum_\mu
\Bigl(
L_\mu \rho L_\mu^\dagger
-
\frac12\{L_\mu^\dagger L_\mu,\rho\}
\Bigr).
\]
After vectorization, $\rho\mapsto |\rho\rangle\rangle$, the Liouvillian becomes a non-Hermitian matrix acting on Liouville space. One then solves
\[
\mathcal L |r_\mu\rangle\rangle = \lambda_\mu |r_\mu\rangle\rangle,
\qquad
\langle\langle \ell_\mu|\mathcal L = \lambda_\mu \langle\langle \ell_\mu|,
\]
with biorthogonal left and right eigenmodes. The complex eigenvalues determine decay rates and oscillation frequencies, while the zero mode is the steady state when it is unique [2602.01375] [2002.11620].

An LEP of order two occurs when $\lambda_1=\lambda_2\equiv\lambda_{\rm EP}$ and the corresponding eigenoperators coalesce, so that the algebraic multiplicity exceeds the geometric multiplicity. Equivalently, the relevant invariant subspace is described by a Jordan block
\[
\mathcal J_2(\lambda_{\rm EP})
=
\begin{pmatrix}
\lambda_{\rm EP} & 1\\
0 & \lambda_{\rm EP}
\end{pmatrix},
\]
or, in rank language, the kernel dimension drops below the algebraic multiplicity. In matrix formulations this is detected by simultaneous eigenvalue coincidence and defectiveness, often through a vanishing discriminant together with a rank condition [2512.06311] [2401.14993] [2605.12867].

The dynamical hallmark of defectiveness is the appearance of polynomial prefactors multiplying exponentials. At a second-order Jordan block one obtains terms of the form $t\,e^{\lambda_{\rm EP} t}$; at higher order, higher-degree polynomials appear. This structure is the Liouvillian analogue of critical damping in non-Hermitian dynamics and underlies the anomalous relaxation signatures associated with LEPs [2101.11553] [2409.08100].

## 2. Relation to non-Hermitian Hamiltonian exceptional points

The effective non-Hermitian Hamiltonian
\[
H_{\rm eff}=H-\frac{i}{2}\sum_\mu L_\mu^\dagger L_\mu
\]
describes no-jump evolution. Its exceptional points are Hamiltonian exceptional points, defined purely at the conditional-wavefunction level. LEPs instead belong to the full Liouvillian and therefore retain the recycling terms $L_\mu \rho L_\mu^\dagger$. This enlargement from Hilbert space to operator space is not a minor technicality: the Liouvillian spectrum is generically richer, and the corresponding singularities need not track those of $H_{\rm eff}$ [2002.11620] [2408.11435].

A useful interpolation between the two descriptions is provided by the hybrid-Liouvillian formalism,
\[
\mathcal L_\eta[\rho]
=
-\,i[H_{\rm eff},\rho]
+
\eta\sum_\mu L_\mu \rho L_\mu^\dagger,
\qquad
\eta\in[0,1],
\]
or equivalently by the jump-weight parameter $q$ in $\mathcal L'(q)$. At $\eta=0$ one recovers the no-jump non-Hermitian evolution; at $\eta=1$ one obtains the full Liouvillian. Varying $\eta$ or $q$ shows explicitly that quantum jumps can move, split, remove, or create spectral degeneracies, and can also change their order [2002.11620] [2506.02902].

Concrete models illustrate all of these possibilities. A driven-lossy qubit with $\Gamma_-=\sqrt{\gamma_-}\sigma_-$, $\Gamma_x=\sqrt{\gamma_x}\sigma_x$, and $\Gamma_y=\sqrt{\gamma_y}\sigma_y$ has a diagonal effective Hamiltonian whose eigenvalues never coalesce, hence no HEPs, yet its Liouvillian exhibits second-order LEPs at $\gamma_x=\gamma_y\pm\omega$ [2401.14993]. In an effective atomic-vapor model, the no-jump sector can host an EP3, while turning on repopulation reduces the degeneracy to an EP2 or removes it entirely, thereby showing that the physical Liouvillian singularity can have a lower order than its NHH precursor [2506.02902]. These results establish that LEPs are not merely re-labeled HEPs, but singularities of the full quantum dissipative generator.

## 3. Dynamical consequences: relaxation, critical damping, and thermodynamic speedup

The long-time approach to stationarity is governed by the Liouvillian spectral gap. In a standard ordering,
\[
\Delta=-\Re \lambda_1
\quad\text{or}\quad
\Delta=-\max_{\alpha\neq 0}\Re\lambda_\alpha,
\]
depending on convention. Several exactly solvable or reduced models show that tuning to an LEP can maximize this gap and therefore minimize the characteristic relaxation time. In a three-level atomic model with coherent drive $\Omega$ and decay $\gamma$, the LEP occurs at $\gamma=2\Omega$, where the static gap reaches $\Delta_{\rm LEP}=\gamma/4=\Omega/2$, and Floquet modulation of the dissipation can enlarge the effective gap beyond the static limit [2305.12745]. In a trapped-ion quantum-battery model, the reduced slow Liouvillian block develops an EP when the cubic discriminant vanishes, and the spectral gap $\Delta$ reaches a local maximum exactly at that defective degeneracy, minimizing the charging time $\tau_s(\epsilon)\simeq (1/\Delta)\ln(1/\epsilon)$ [2605.12867].

Higher-order LEPs generate critical-damping-type dynamics. In the two-qubit thermal machine analyzed by Khandelwal, Brunner, and Haack, a third-order LEP occurs when $\Delta\Gamma^2=4g^2$. At that point three eigenvalues merge at $\bar\lambda=-\Gamma/2$, the corresponding eigenmatrices coalesce, and the exact solution acquires both $t\,e^{\bar\lambda t}$ and $t^2 e^{\bar\lambda t}$ terms. The resulting relaxation is the open-system analogue of a critically damped oscillator and is asymptotically faster than the overdamped nonoscillatory alternatives considered in that model [2101.11553]. Closely related behavior survives beyond the conventional master-equation regime: in a dissipative double quantum dot, the exact Heisenberg-equation generator and the corresponding master-equation Liouvillian share the same square-root structure and the same critical coupling $g_{\rm EP}=|\Gamma_1-\Gamma_2|/4$, and the critical-damping signature persists beyond Born-Markov validity [2409.08100].

These spectral effects have direct thermodynamic and control implications. Encircling an LEP in a single-ion quantum heat engine yields positive net work, with larger net work obtained closer to the LEP, while LEP conditions in sideband and EIT cooling reproduce analytically optimal detuning and Rabi-frequency conditions for fast cooling [2302.13450] [2305.12745]. The common mechanism is dynamical rather than static: spectral reorganization of slow Liouvillian modes changes the route to the steady state.

## 4. Encircling, chirality, and non-Markovian topology

Parameter-space encircling of LEPs generalizes the branch-point physics of non-Hermitian Hamiltonians to open-system superoperators. For a time-dependent Liouvillian $\mathcal L(t)$, the propagator is
\[
V(t)=\mathcal T\exp\!\left[\int_0^t \mathcal L(t')\,dt'\right],
\qquad
\rho(t)=V(t)\rho(0).
\]
Near a second-order LEP, the instantaneous spectral sheets are joined by a square-root branch cut, and slow encircling can induce directional mode conversion among decay modes. In Liouvillian settings, however, this chirality is typically transient because the steady state eventually dominates the dynamics [2408.11435].

This transient nature was made explicit in single-photon interferometric simulations of open-system dynamics. In that work, clockwise and counterclockwise encircling of a Liouvillian EP produce different final states only for intermediate encircling times. When dephasing opens a nonzero Liouvillian gap, sufficiently long protocols drive both directions back to the same instantaneous steady-state sheet, and the chirality vanishes. The experiment further reported a scaling collapse of the chirality $C$ as a function of $\gamma_2 T^{1/\nu}$, with fitted exponent $\nu\approx 1.7221$ [2501.10349]. This result sharpens an important distinction: LEP encircling can control transient decay-mode populations even when it does not alter the ultimate steady state.

A qualitatively different topological structure arises in the non-Markovian pseudomode construction of a qubit coupled to a reservoir with memory. There the extended Liouvillian
\[
\mathcal L\rho
=
i\bigl[\rho H_{NH}^\dagger-H_{NH}\rho\bigr]
+
\kappa b\rho b^\dagger
\]
hosts two coinciding LEP2s at $(\Omega=\kappa/4,\Delta=0)$. For a closed loop
\[
\Omega(k)=\frac{\kappa}{2}+r\cos k,\qquad
\Delta(k)=r\sin k,\qquad
r>\frac{\kappa}{4},
\]
the generalized winding number
\[
\mathcal W_j
=
\frac{1}{2m\pi}
\int_0^{2m\pi}
dk\,
\frac{d}{dk}\arg[\lambda_j(k)-\lambda_{\rm EP}]
\]
takes two distinct values on a single encircling path:
\[
\mathcal W_{1}=\mathcal W_{3}=+\frac12,
\qquad
\mathcal W_{2}=\mathcal W_{4}=-\frac12.
\]
This “hybrid” invariant was measured in a superconducting circuit where a transmon was sideband-coupled to a decaying resonator acting as a Lorentzian reservoir. Repeating the experiment at large $\kappa$ removed the distinction, so the dual half-integer topology was identified as a purely non-Markovian effect [2512.06311].

## 5. Higher-order LEPs, asymptotic exceptional steady states, and geometric methods

LEPs are not restricted to simple EP2 scenarios. In a driven two-level Liouvillian with decay and dephasing,
\[
\dot\rho
=
-\,i[H,\rho]
+
\gamma_0 L_{|g\rangle\langle e|}[\rho]
+
\gamma_\phi L_{|e\rangle\langle e|}[\rho],
\]
the $4\times4$ Liouvillian matrix has one steady eigenvalue and three nontrivial modes. A third-order LEP occurs when
\[
\frac{\Delta}{\Omega}=\pm \frac{1}{\sqrt 8},
\qquad
\frac{\gamma}{\Omega}
=
\frac{\sqrt{27/2}}{|1-2\alpha|},
\qquad
\alpha=\frac{\gamma_0}{\gamma},
\]
so that all three nonzero modes coalesce. Because the pure-decay and pure-dephasing Liouvillians do not commute, the LEP position moves with the mixing parameter $\alpha$, and at $\alpha=1/2$ the EP is pushed to infinity [2512.01217]. This is a direct illustration of the sensitivity of high-order LEPs to the structure, not just the magnitude, of dissipation.

In linear bosonic and continuous-variable settings, the possible order can be much larger. For dissipative linear bosonic systems, an $n$th-order HEP of the effective non-Hermitian Hamiltonian implies LEPs of order $m\ge n$ in the Liouvillian, and steady-state coherence functions diagnose them through non-exponential time dependence and higher-power Lorentzian lineshapes: squared Lorentzians in power spectra and cubic Lorentzians in intensity-fluctuation spectra were identified for second- and third-order signatures in a bimodal cavity with incoherent coupling [2006.03557]. For a damped harmonic oscillator in a generic environment, the LEP occurs when the modified oscillator frequency vanishes, $\Omega=0$, corresponding to critical damping in Caldeira-Leggett and Markovian Hu-Paz-Zhang equations or to an infinite-effective-mass limit in a modified Kossakowski-Lindblad equation. At the EP, each fixed-$N$ sector exhibits an $(N+1)$-fold degeneracy and a single Jordan block of order $N+1$ [2304.05792].

Exactly solvable fermionic models show a complementary many-body route to high-order LEPs. In a class of quadratic open fermion models treated by third quantization, the Liouvillian develops an EP whose order approaches system size and appears as a quasisteady state with gapless spectrum. Adding many-body quantum-jump perturbations lifts that degeneracy and opens a finite Liouvillian gap with fractional power-law scaling whose exponent depends sensitively on the perturbation form [2602.00486]. Discrete-time open dynamics can also host LEPs: in two-qubit brickwork CPTP circuits, the stroboscopic transfer matrix develops a Jordan block at a surface defined by the vanishing of an explicit square root $Q$, and the resulting time-domain signal acquires an $n\,\mu_0^n$ factor in the defective subspace [2510.10629].

The geometry of LEP splitting has also been systematized algebraically. Newton polygons and tropical geometry encode the Puiseux exponents of eigenvalue branches near a defective point. A polygon edge of slope $-1/n$ signals the leading $\epsilon^{1/n}$ splitting of an $n$th-order EP, while multiple lower faces reveal anisotropic or hybrid splittings, such as coexistence of linear and cube-root scaling or square-root splitting with frozen modes, depending on the perturbation direction [2510.08156]. At the opposite end of the spectrum, zero modes obey a no-go theorem in finite dimension: steady-state eigenspaces at physical jump strength are semisimple. Even so, dissipative many-body systems can approach an “asymptotic exceptional steady state” in the thermodynamic limit, with $\Delta(N,1)\to 0$ and normalized overlap $|\nu|^2\to 1$ between the steady state and first metastable mode, exhibiting exponential scaling in a dissipative 3SAT solver and polynomial scaling in AKLT-type symmetry-protected topological state preparation [2504.02937].

## 6. Diagnostics, experimental platforms, and conceptual boundaries

LEPs are diagnosed by reconstructing or probing the Liouvillian rather than only the conditional Hamiltonian. Quantum process tomography reconstructs the short-time transfer matrix $S\simeq I+dt\,\mathcal L$ and hence the experimental Liouvillian; resolvent spectroscopy identifies higher-order poles and super-Lorentzian lineshapes; coherence functions reveal Jordan-block-induced polynomial decay; and full state tomography of reduced or extended systems can track the complex Liouvillian eigenvalues along parameter loops [2401.14993] [2602.01375].

| Platform | Diagnostic | Key LEP feature |
|---|---|---|
| IBM superconducting processor | Quantum process tomography | LEPs without any HEP [2401.14993] |
| Superconducting qubit + decaying resonator | Joint-state tomography in extended space | Non-Markovian twofold LEP2 with $\pm 1/2$ winding [2512.06311] |
| Single trapped ion | Time-resolved tomography and parameter scans | Heat-engine enhancement and third-order LEPs [2302.13450] [2512.01217] |
| Collective spin with polarized bath | Frequency-resolved resolvent spectroscopy | Super-Lorentzian response and state-dependent visibility [2602.01375] |

Spectroscopy introduces additional subtleties. In the collective-spin model, the emission spectrum
\[
S(\omega)
=
\frac{1}{\pi}\Re\,\mathrm{Tr}\!\left[J_-(i\omega-L)^{-1}(\rho_{\rm ss}J_+)\right]
\]
contains both simple and second-order poles near an EP, leading to a super-Lorentzian correction
\[
S(\omega)\simeq
\frac{1}{\pi}
\left[
a\frac{1}{\Delta^2+\gamma^2}
+
b\frac{\gamma^2-\Delta^2}{(\Delta^2+\gamma^2)^2}
\right].
\]
However, weak $U(1)$ symmetry sectors impose selection rules, and the overlap of the source state with the Jordan chain can strongly suppress the $b$ term. As a result, steady-state fluorescence may conceal an EP that becomes unambiguous for infinite-temperature or random initial states [2602.01375]. This state dependence is a recurrent theme in LEP diagnostics.

The conceptual status of LEPs outside time-local descriptions remains unsettled. One line of work argues that exact non-Markovian dynamics, as described by NEGF Dyson equations or non-Markovian master equations with memory kernels, does not admit a fixed finite-dimensional Liouvillian matrix and therefore does not allow the introduction of exceptional points in the strict spectral sense [2401.04011]. Another line constructs non-Markovian LEPs in an extended Liouvillian acting on system-plus-pseudomode space and reports experimentally accessible non-Markovian topology in that enlarged generator [2512.06311]. The literature therefore contains two distinct viewpoints: LEPs as singularities of a time-local generator, and LEPs as singularities of an extended effective generator that embeds memory degrees of freedom. This suggests that the meaning of “non-Markovian LEP” depends on what is taken to be the fundamental dynamical operator.

Across these settings, the central theme is stable: LEPs are spectral singularities of dissipative quantum generators, not merely of no-jump Hamiltonians. They reorganize relaxation channels, modify spectral gaps, produce polynomial decay, reshape topological encircling phenomena, and, in higher-order or many-body settings, connect to asymptotic criticality and computational complexity. Their modern study now spans Markovian and pseudomode descriptions, bosonic and fermionic many-body systems, continuous and discrete time, and both tomographic and spectroscopic diagnostics [2408.11435].

Source: https://www.emergentmind.com/topics/liouvillian-exceptional-point