Liouvillian Eigenstates in Open Systems

• Liouvillian eigenstates are generalized eigenoperators that enable a complete spectral decomposition of open-system dynamics.
• They reveal key phenomena such as exceptional points, flat bands, and skin effects through biorthogonal structures.
• Their rigorous analysis advances control over dissipation, decoherence, and topological transitions in quantum and classical settings.

Liouvillian eigenstates are generalized eigenoperators of the Liouvillian superoperator, which governs the evolution of density matrices in both classical and quantum dynamical systems. Their rigorous characterization, structure, and physical significance have become central in diverse settings, from classical Hilbert-space mechanics to open quantum systems, non-Hermitian many-body settings, and frameworks with rich symmetry and topological properties.

1. Definition and General Properties

The Liouvillian operator $\mathcal{L}$ acts on the space of density operators $\rho$ via

$\frac{d\rho}{dt} = \mathcal{L}[\rho].$

In quantum settings with Markovian open dynamics, $\mathcal{L}$ typically takes the Lindblad form: $$\mathcal{L}[\rho] = -i[H, \rho] + \sum_\alpha \gamma_\alpha \left(L_\alpha \rho L_\alpha^\dagger - \frac{1}{2}\{L_\alpha^\dagger L_\alpha, \rho\}\right),$$ where $H$ is the system Hamiltonian and $L_\alpha$ are jump operators.

Liouvillian eigenstates are operators $r_n$ (right) and $\ell_n$ (left) satisfying

$$\mathcal{L}[r_n] = \lambda_n r_n, \quad \mathcal{L}^\dagger[\ell_n] = \lambda_n^* \ell_n,$$

with $\lambda_n$ complex, encoding both oscillatory (Re $\lambda_n$) and dissipative (Im $\lambda_n$) dynamics (Richter et al., 24 Nov 2025).

Biorthogonality and completeness are imposed via

$$\mathrm{Tr}(\ell_m^\dagger r_n) = \delta_{mn}, \quad \sum_n |r_n)(l_n| = \mathbb{I},$$

with operator-ket notation $|r_n)$.

Physical states (in particular steady states) correspond to the subspace of right eigenoperators with $\lambda_n = 0$. For non-Hermitian $\mathcal{L}$, the full spectrum may exhibit degeneracies, exceptional points, and biorthogonal structure (Richter et al., 24 Nov 2025, Honda et al., 2010, Tay, 2023).

2. Construction and Structure in Classical and Quantum Systems

2.1. Koopman–von Neumann (KvN) Classical Formalism

In the KvN Hilbert-space embedding of classical mechanics, Liouvillian eigenstates $\psi(q,p)$ solve

$$\hat{L} \psi(q, p) = i\hbar \{ H, \psi \} = \lambda \psi(q,p),$$

where $\{, \}$ denotes the Poisson bracket (Amin et al., 11 Dec 2025). The spectrum can be made discrete or continuous depending on the gauge and the choice of separation variables, with product-form eigenfunctions for canonical ensembles and time–energy conjugacy manifesting in the $(\tau, H)$ representation. KvN eigenstates admit Hilbert-space orthonormality and support superpositions after fixing the associated gauge.

2.2. Lindbladian Systems and Quantum Markovian Dynamics

For Lindblad dynamics, Liouvillian eigenstates form a biorthogonal basis of operator-space. The spectral decomposition allows arbitrary $\rho$ (including pure-state projectors from quantum trajectories) to be expanded as

$$\rho = \sum_n c_n r_n, \qquad c_n = \mathrm{Tr}(\ell_n^\dagger \rho)$$

with time evolution $c_n(t) = e^{\lambda_n t} c_n(0)$ (Richter et al., 24 Nov 2025).

Right eigenstates encode the decay modes of the system, with the steady state $\rho_\mathrm{ss}$ corresponding to the zero mode. The non-Hermitian nature of $\mathcal{L}$ leads to a nontrivial biorthogonal structure, including gauge redundancies (the $r_n$ and $\ell_n$ may be independently normalized up to nonzero complex scalars).

2.3. Flat Bands, Skin Effects, and Nontrivial Spectra

Liouvillian eigenstates reflect special structural phenomena:

• Flat Bands: In systems with “Liouvillian flat bands,” all modes in a band have the same eigenvalue, leading to degenerate relaxation rates and dynamically localized “compact localized normal master modes” (CLNMMs) that are strictly supported at certain sites (Liu et al., 2023).
• Liouvillian Skin Effect: Non-Hermitian Liouvillians may exhibit the macroscopic localization of many eigenstates at a boundary (skin modes), with localization lengths and boundary sensitivity controlled by spectral topology and boundary conditions (Mao et al., 2024, Shigedomi et al., 23 May 2025). Loss of the effect is triggered by specific boundary hoppings or temperature.
• Nontrivial Spectra: Complex eigenvalues that do not arise from simple Hamiltonian differences may appear, as in the T-type quantum dot model, arising from correlated bra–ket (two-particle) dynamics and captured via effective two-body Hamiltonian mappings (Nakano et al., 2010).

3. Algebraic and Analytical Structure

3.1. Explicit Eigenbasis: Harmonic Oscillator Case

For quantum harmonic oscillators, the explicit ladder of Liouvillian eigenstates is constructed with biorthogonal right and left eigenoperators, labeled by dissipation and phase-rotation quantum numbers $(j, k)$. The right eigenoperators can be written as

$$\rho^{(R)}_{j, k} = \mathcal{N}_{j, k} \sum_{m=0}^{l} (-1)^m \frac{(k+m)! (l-m)!}{m! l!} (a^\dagger)^{k+m} e^{-\beta \omega a^\dagger a} a^m,$$

with spectrum

$$\Lambda_{j, k} = -(\gamma - \gamma') j - i \omega k,$$

providing a complete resolution with well-defined ladder operators generating the full basis (Honda et al., 2010, Tay, 2020).

3.2. Quadratic Bosonic Systems, Exceptional and Diabolical Points

For generic quadratic Hamiltonians, the Liouvillian eigenproblem can be mapped to the diagonalization of a dynamical matrix $A$ in operator moment space, with eigenvalues $\lambda_i$ and associated eigenoperators $O_i$ explicitly constructed from the eigenvectors $v_i$ of $A$ (Jr et al., 2022). Coalescence of eigenvalues and eigenvectors produces exceptional points (EPs) with nontrivial Jordan block structure, fully classified for both single-mode and two-mode systems (Tay, 2023, Jr et al., 2022). These EPs govern non-exponential dynamics and can be inherited, genuine, induced, or hidden according to system reductions or coupling structure.

4. Topological, Symmetry, and Many-Body Aspects

4.1. Spectral Topology and Skin Modes

Liouvillian spectra can support nontrivial winding numbers (e.g., $\mathbb{Z}$ or $\mathbb{Z}_2$), directly linking bulk topology to the presence and robustness of edge-localized eigenstates (Liouvillian skin modes) (Shigedomi et al., 23 May 2025, Koch et al., 20 Nov 2025). For open Floquet chains, stroboscopic Floquet-Liouvillians display band and system winding numbers, with quantum jumps producing "jump-induced" topological phases and robust edge transport. In open XXZ chains, the Liouvillian in the diagonal sector can be solved exactly by Bethe ansatz, supporting analytically tractable skin modes, whose localization length and boundary sensitivity are determined by the nonreciprocal hopping and the precise form of boundary conditions (Mao et al., 2024).

4.2. Quasi-Spin Symmetry and Multiplet Structure

In nonlinear oscillator systems, e.g., the squeeze-driven Kerr oscillator, Liouvillian eigenstates reorganize into $su(2)\otimes su(2)^*$ multiplets at integer detuning/frequency ratios. The spectral structure forms double-ellipsoidal shells in the complex plane, with degeneracies controlled by the multiplet structure and relaxation times by the Liouvillian spectral gap (Iachello et al., 2024).

4.3. Symmetry-Induced Degeneracies

System symmetries induce Liouvillian block diagonalization, producing multiple steady states in distinct symmetry sectors. When strong symmetries are present, true steady-state subspaces can be constructed using projection operators or by orthonormalizing the kernel of $\mathcal{L}$ via Gram–Schmidt and enforcing positivity. Large deviation techniques allow the extraction of extremal current-carrying steady states in nonequilibrium settings (Thingna et al., 2021).

5. Dynamical and Statistical Implications

Expansion in the Liouvillian eigenbasis enables detailed analysis of both trajectory and ensemble dynamics:

• Any state can be decomposed in terms of decay eigenmodes, with projectors onto these modes evolving as $e^{\lambda_n t}$.
• "Quasiprobability" weights $p_n = c_n d_n$ for trajectory-eigenstate overlaps encode the localization/delocalization of individual quantum trajectories in Liouvillian space (Richter et al., 24 Nov 2025), with statistical properties—such as the inverse participation ratio (IPR) and center of mass—revealing how trajectories explore the decay spectrum and correlating with the purity of the steady state.
• In degenerate or topologically nontrivial cases, dynamics can be dominated by slowest-decaying or edge-localized eigenstates, with observable effects such as unidirectional transport, dynamical charge accumulation at boundaries, or steady-state localization under the Liouvillian skin effect (Shigedomi et al., 23 May 2025, Koch et al., 20 Nov 2025, Mao et al., 2024, Liu et al., 2023).

6. Schematic Comparisons of Key Structural Phenomena

Physical Setting	Liouvillian Spectrum/Eigenstates	Special Phenomena
Harmonic Oscillator (Lindblad)	Discrete, biorthogonal ladder, explicit forms	Rigged Hilbert space, ladder superoperators, spectral completeness
Quadratic Bosonic (General)	Mapped to moments/dynamical matrix eigenproblem	Exceptional/diabolical points, Jordan blocks, hidden degeneracies
Koopman–von Neumann (Classical)	Continuous/discrete via gauge, product-form	Canonical ensembles, uncertainty relations, coherent superpositions
Many-Body NH Chains (Skin Effect)	Macroscopic boundary-localized modes (OBC)	Sensitivity to boundary, Bethe ansatz solvability, skin fragility
Floquet/Open Chains (Topology)	Bands, winding numbers, outlier modes, bulk-edge	Jump-induced topological transition, NH skin, unidirectional pumping
Squeeze-Driven Kerr Oscillator	$su(2)\otimes su(2)^*$ multiplets, ellipsoids	Quasi-spin symmetry, long-lived modes, double-ellipsoidal structure

7. Summary and Outlook

Liouvillian eigenstates provide a complete, generically biorthogonal basis for the spectral analysis of evolution generators in both classical and quantum dynamical systems. Their explicit structure encodes all physical dissipation, decoherence, and oscillatory dynamics, with intricate mathematical features—exceptional points, topological winding, boundary sensitivity, and multiplet degeneracies—emerging in many-body, non-Hermitian, and nonlinear contexts. The ability to construct and resolve Liouvillian eigenstates underpins fundamental advances in the understanding of relaxation processes, nonequilibrium steady states, quantum trajectory statistics, and quantum-to-classical connections (Amin et al., 11 Dec 2025, Richter et al., 24 Nov 2025, Koch et al., 20 Nov 2025, Nakano et al., 2010, Honda et al., 2010, Tay, 2023). Ongoing research continues to probe their relevance for phase transitions, dissipative engineering, and the control of open-system quantum information.