Liouvillian Spectral Resolution

• Liouvillian spectral resolution is the explicit decomposition of the Liouvillian superoperator into eigenvalues, right/left eigenmodes, and biorthogonal projectors, crucial for understanding open system dynamics.
• It employs operator-theoretic and algebraic methods, including Jordan block analysis and Floquet techniques, to address phase transitions and topological properties in nonequilibrium systems.
• This spectral approach provides actionable insights into relaxation processes, transport phenomena, and the emergence of steady states in complex quantum and classical dissipative systems.

A Liouvillian spectral resolution refers to the explicit decomposition of a Liouvillian superoperator—governing quantum or classical open-system dynamics—into its spectral components: eigenvalues, right/left eigenmodes, and biorthogonal projectors. This construction is central to the analysis of nonequilibrium steady states, dissipative phase transitions, topological phases in open systems, and the dynamical properties of Lindbladian and more general non-Hermitian evolutions. Both operator-theoretic and algebraic approaches to Liouvillian spectrum are present in the literature, with technical frameworks suited to Markovian master equations, spectral theory of differential operators, and infinite-dimensional rigged spaces, depending on context.

1. Liouvillian Superoperator Formalism and Spectral Structure

Given a quantum system with density matrix $\rho$, the general generator of Markovian evolution is the time-dependent (or time-independent) Liouvillian superoperator $\mathcal{L}$, most commonly in the Lindblad (GKLS) form:

$$\dot{\rho} = \mathcal{L}[\rho] = -i[H, \rho] + \sum_\mu \gamma_\mu \left(L_\mu \rho L_\mu^\dagger - \frac12 \{L_\mu^\dagger L_\mu,\rho\} \right),$$

where $H$ is the system Hamiltonian and $L_\mu$ are jump (collapse) operators.

The Liouvillian acts linearly, typically on the space of trace-class operators or, after vectorization, as a non-Hermitian matrix on an enlarged Hilbert (or Liouville) space. Its spectral resolution takes the general form

$\mathcal{L} = \sum_n \lambda_n P_n,$

with right eigenmodes $\mathcal{L}[\rho_n] = \lambda_n \rho_n$, left eigenmodes $$\mathcal{L}^\dagger[\sigma_n] = \lambda_n^* \sigma_n$$, and biorthogonality $\mathrm{Tr}(\sigma_n^\dagger \rho_m) = \delta_{nm}$ (Minganti et al., 2021). The spectrum is, in general, complex and may display nontrivial Jordan block (defective) structure at exceptional points (Molina, 1 Feb 2026, Kopciuch et al., 3 Jun 2025).

2. Spectral Decomposition, Projectors, and Jordan Structure

In finite-dimensional systems or compact operators, the Liouvillian admits complete biorthogonal spectral resolutions:

• If diagonalizable:

$$\mathcal{L} = \sum_i \lambda_i \Pi_i, \qquad R(z) = (z-\mathcal{L})^{-1} = \sum_i \Pi_i (z-\lambda_i)^{-1}$$

with projectors $\Pi_i$ and meromorphic resolvent $R(z)$ (Molina, 1 Feb 2026).

• In the presence of exceptional points (defective eigenvalues), the Jordan structure manifests as higher-order poles in the resolvent:

$$R(z) = \sum_i \sum_{k=1}^{r_i} \Pi_{i,k} (z-\lambda_i)^{-k},$$

with $r_i$ denoting the order of the largest Jordan block at $\lambda_i$, and generalized projectors $\Pi_{i,k}$ built from nilpotent chains. This leads to non-exponential, polynomial prefactors (e.g., $t^{r_i-1}e^{\lambda_i t}$) in dynamical evolution and response (Molina, 1 Feb 2026).

Biorthogonal completeness and the explicit action of projectors on observables enable precise tracking of each eigenmode's contribution to physical quantities (Honda et al., 2010).

3. Liouvillian Spectral Theory in Periodic, Topological, and Driven Systems

Liouvillian spectral resolution is essential for analyzing topological and dynamical phases in non-equilibrium open systems subjected to periodic driving (Floquet systems). In "Liouvillian topology and non-reciprocal dynamics in open Floquet chains" (Koch et al., 20 Nov 2025), the full time-periodic Liouvillian $\mathcal{L}(t)$ is discretized into a Floquet propagator $F$:

$$F = U(T,0) = e^{\mathcal{L}_3 T/3} e^{\mathcal{L}_2 T/3} e^{\mathcal{L}_1 T/3},$$

acting on vectorized density matrices. The Floquet-Liouvillian spectral decomposition yields eigenvalues $\eta_a(k)$ and eigenvectors $\lvert R_a(k) \rangle\rangle$, $\langle\langle L_a(k)\rvert$ for each momentum block, and the full stroboscopic dynamics is constructed as

$$F(k) = \sum_a \eta_a(k) \lvert R_a(k)\rangle\rangle \langle\langle L_a(k)\rvert.$$

This enables the definition of Liouvillian spectral winding numbers for each band:

$$W(\eta_a) = -\frac{i}{2\pi} \int_{-\pi}^{\pi} \eta_a(k)^{-1} \partial_k \eta_a(k) dk,$$

and a system-wide winding $W_S$ characterizing topological transport properties and their relation to observable dynamics and boundary effects (Koch et al., 20 Nov 2025).

A key discovery is the existence of quantum-jump-induced topological phases, which occur only at specific regimes of the Lindblad jump rate $\gamma$ and have no counterpart in non-Hermitian Hamiltonian approaches; the full Liouvillian spectrum encodes both coherent and jump-induced bands, including new "outlier" spectral branches whose winding is distinct (Koch et al., 20 Nov 2025).

4. Liouvillian Spectrum, Nonequilibrium Steady States, and Transport

Open-system steady states (NESS) and transport are directly governed by the Liouvillian's spectrum and eigenvectors:

• The unique NESS corresponds to the zero eigenmode $\mathcal{L}[\rho_\mathrm{NESS}] = 0$, projected out from the full spectrum by contour integration of the resolvent around $z=0$ (Tanaka et al., 2010).
• In chain models coupled to baths, the full Liouvillian spectral decomposition allows explicit construction of NESS, incorporates spatial bath-system correlations, and separates Landauer-type transport (first-order spatial correlations, single-bath connection) from true many-body, non-Landauer conductivity (second-order correlations) (Tanaka et al., 2010).
• For quadratic fermionic and bosonic systems, advanced "third-quantization" techniques reduce the Liouvillian spectrum to eigenvalues of effective one-body matrices, enabling an exact mapping to the spectrum of single-particle rapidities and efficient computation of the relaxation gap (Pouranvari, 26 Nov 2025, Wang et al., 2024, Yuan et al., 2020).

Table: Liouvillian Spectral Features in Various Physical Contexts

System Class	Spectrum Structure	Physical Implications
Periodically driven	Floquet bands, windings	Topological transport, skin effect
Quadratic fermion chain	Rapidity matrix, level statistics	Transport/localization diagnostics
Atomic vapor	Liouvillian/Hamiltonian EPs, jumps	Spectral singularities, quantum jumps
Harmonic oscillator	Discrete, ladder structure	Mode-resolved damping, expectation flow

5. Spectral Singularities, Collapse, and Dissipative Phase Transitions

Liouvillian spectral theory has established the existence of non-Hermitian degeneracies (exceptional points), spectral collapse, and non-trivial criticality in dissipative quantum systems:

• Liouvillian exceptional points (LEPs) of order $n$ correspond to $n$-fold degenerate, non-diagonalizable eigenvalues, leading to non-exponential, power-law temporal response and higher-order poles in frequency-resolved spectra (Kopciuch et al., 3 Jun 2025, Molina, 1 Feb 2026). Experimental signatures include the emergence of super-Lorentzian line shapes and state-dependent visibility in emission spectra.
• Liouvillian spectral collapse describes the non-analytic closure of the real parts of infinitely many decay modes at a dissipative phase transition ($\lambda_j\to0$ for $j\ge1$), leading to anomalously slow relaxation, emergent multistability, and hysteresis even in the absence of symmetry breaking (Minganti et al., 2021).
• In quadratic or exactly solvable Lindbladians, the precise spectrum (gaps, bands, exceptional points) prescribes the dynamical landscape, relaxation timescales, and phase behavior (Wang et al., 2024, Yuan et al., 2020).

6. Algebraic, Analytic, and Algorithmic Aspects of Liouvillian Spectral Sets

For algebraic differential operators and non-quantum settings, the "Liouvillian spectral set" $\Sigma$ refers to the locus in parameter space where a linear ODE admits explicit (Liouvillian) integrable solutions. This set is shown to be a countable union of non-accumulating algebraic varieties, with precise structural theorems available for polynomial coefficient ODEs via Kovacic's algorithm and the Asymptotic Iteration Method (Acosta-Humánez et al., 2019).

In the analytic domain, classical Liouville transformations and Delsarte transmutation operators provide constructive bridges between Sturm–Liouville operators and canonical (Schrödinger) forms, enabling high-precision numerical and analytic spectral resolution and linking formal power series (SPPS) representations to eigenfunction expansions (Kravchenko et al., 2014, Jin, 2017). The spectral data (eigenvalues/eigenfunctions/spectral measures) thus recovered directly inform probabilistic and fractal geometric properties of associated diffusions.

7. Methodologies for Liouvillian Spectral Computation and Applications

Methodological advances in Liouvillian spectral resolution include:

• Block-diagonalization and vectorization for translationally invariant and Floquet systems (Koch et al., 20 Nov 2025).
• Third-quantization and Majorana covariance formalism for exactly solving quadratic (Gaussian) models (Pouranvari, 26 Nov 2025, Wang et al., 2024).
• Explicit biorthogonal construction of right/left eigenoperators and associated rigged spaces to capture the full spectral content in infinite-dimensional settings, e.g., for the quantum harmonic oscillator (Honda et al., 2010).
• Machine-learning-based variational approaches (e.g., RBM ansatz) for accessing Liouvillian gaps and dominant decay modes in high-dimensional many-body systems, benchmarked against Bethe-ansatz results when possible (Yuan et al., 2020).
• Algorithmic recursion based on closed hierarchies of correlation functions, extending Wick's theorem, to reconstruct the entire spectrum for quadratic and certain quartic Lindbladians (Wang et al., 2024).

These approaches enable the full exploitation of Liouvillian spectral data for diagnostics of unidirectional transport, localization, phase transitions, topological phenomena, and nonequilibrium response in open quantum systems.