Liouvillian Skin Effect in Open Quantum Systems

• Liouvillian Skin Effect (LSE) is a non-Hermitian phenomenon in open quantum systems characterized by boundary-localized eigenmodes under open boundary conditions.
• It arises from point-gap topology, resulting in anomalously slow, boundary-limited relaxation that scales with both system size and localization length.
• LSE offers practical insights for experimental quantum setups, enabling controlled state-preparation and novel dissipative dynamics in many-body systems.

The Liouvillian Skin Effect (LSE) is a non-Hermitian phenomenon manifesting in open quantum systems governed by Lindblad master equations. Under open boundary conditions, a macroscopic fraction of the non-equilibrium (right) eigenmodes of the Liouvillian superoperator become exponentially localized at one edge of the system. This boundary accumulation leads to anomalously slow, boundary-limited relaxation, which cannot be deduced simply from the bulk spectral gap. The LSE arises from point-gap topology in the complex spectrum of the Liouvillian and is robustly connected to irreversible transport, topological invariants, and extreme spectral sensitivity to the choice of boundary conditions (Haga et al., 2020, Hamanaka et al., 2023, Yang et al., 2022).

1. Mathematical Foundation and Mechanism

Consider a general driven-dissipative quantum system, with density matrix ρ evolving under the Lindblad equation: $$\frac{d\rho}{dt} = \mathcal{L}[\rho] = -i[H, \rho] + \sum_\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 quantum jump operators describing system-bath coupling. The Liouvillian superoperator $\mathcal{L}$ is generically non-Hermitian due to dissipative terms, and acts linearly on the $D^2$-dimensional space of operators.

The eigenvalue problem for $\mathcal{L}$ is defined by

$$\mathcal{L}(\rho^R_j) = \lambda_j \rho^R_j,\quad \mathcal{L}^\dagger(\rho^L_j) = \lambda_j^* \rho^L_j,$$

with right and left eigenmodes biorthogonal as $$\mathrm{Tr}[(\rho^L_j)^\dagger \rho^R_k] = \delta_{jk}$$. If $|\mathrm{Re}\,\lambda_0|=0$ corresponds to the unique steady state $\rho_\mathrm{ss}$ and $H$0 is the Liouvillian spectral gap, then the relaxation rate is often naively expected to be $H$1.

The LSE occurs if, under open boundary conditions (OBC), the right eigenmodes $H$2 with $H$3 take the asymptotic form

$H$4

i.e., are exponentially localized at a system boundary, with localization length $H$5 (Haga et al., 2020). This behavior renders the system extremely sensitive to boundary conditions: under periodic boundary conditions (PBC), the spectrum forms closed loops in the complex plane; under OBC, spectral loops collapse, and eigenmodes pile at the edge.

The LSE is intimately related to the non-Hermitian skin effect (NHSE) of effective non-Hermitian Hamiltonians, but can persist even in regimes where the latter fails, especially due to the presence of quantum jumps and in many-body or correlated systems (Hamanaka et al., 2023, Feng et al., 2023).

2. Topological Diagnosis and Boundary Condition Sensitivity

The emergence of the LSE is topologically protected by a point-gap winding number in the complex Liouvillian spectrum. For a block-diagonal (Bloch-type) Liouvillian $H$6 under PBC, one defines

$H$7

where $H$8 is a generalized boundary twist. A nonzero $H$9 indicates a spectral loop encircling $L_\alpha$0, implying that OBC must induce a distinct eigenvalue flow and mode localization (skin effect) (Hamanaka et al., 2023, Feng et al., 2024).

Boundary conditions thus crucially determine the physical manifestation of the LSE:

• Under OBC, localized skin modes and boundary-sensitive steady states appear.
• Under PBC, the modes are extended and the LSE vanishes.

Fragility or robustness of the LSE to boundary perturbations is system-dependent. For instance, in exactly solvable many-body models, the LSE survives counter-flow boundary hopping but is destabilized by even infinitesimal co-flow boundary hopping (Mao et al., 2024).

3. Anomalous Relaxation Scaling

Unlike conventional dissipative systems, LSE systems exhibit a longest relaxation time $L_\alpha$1 determined not merely by the Liouvillian gap $L_\alpha$2, but also by the localization length $L_\alpha$3 and system size $L_\alpha$4: $L_\alpha$5 (Haga et al., 2020). For fixed $L_\alpha$6, $L_\alpha$7 diverges linearly with $L_\alpha$8 if $L_\alpha$9 is finite. This is due to the exponentially small Hilbert-Schmidt overlap between a localized left mode and an initial state at the opposite boundary; information or particles must propagate a macroscopic distance to excite the slowest decaying mode pinned at the boundary.

Numerically, this scaling is observed in asymmetric dissipative hopping chains, where the steady profile is $\mathcal{L}$0 and the time for the boundary occupation to relax grows with $\mathcal{L}$1 (Haga et al., 2020). In interacting systems, the same mechanism is operative, but now the localization length depends on the imaginary part of complex-valued interactions $\mathcal{L}$2 (Hamanaka et al., 2023).

More exotic scaling, such as $\mathcal{L}$3 in systems with gradient non-reciprocal hopping, is found when the standard LSE paradigm breaks down (Wang et al., 2023).

4. Physical Manifestations and Experimental Proposals

The LSE produces highly asymmetric dynamics:

• Directional, boundary-focused accumulation of probability or particles in transient and steady-state evolution.
• "Cutoff phenomena" in large systems, where no significant decay occurs until a system-size dependent time, after which a sudden relaxation ensues (Feng et al., 2023).
• Boundary-sensitive entanglement features, such as transitions from area-law to log-law scaling of steady-state entropy when delocalization occurs via interchain coupling (Feng et al., 2024).

Experimental analogs are proposed in engineered cold atom setups (with tuneable boundary dissipation), trapped-ion sideband cooling, and mesoscopic transport in magnetic multilayers (Li et al., 2023, Cai et al., 2024). The engineering of additional dissipative channels can increase Liouvillian gaps, optimizing state-preparation and cooling efficiency by leveraging LSE-based boundary transport (Cai et al., 2024).

In some platforms, the interplay of disorder and local reciprocity suppresses skin accumulation, but can yield erratic, sample-dependent bulk localization and Sinai-type subdiffusive transport, distinguishing Liouvillian dynamics from conventional non-Hermitian Hamiltonian models (Longhi, 16 Feb 2026).

5. LSE in Many-Body and Correlated Systems

The LSE persists in many-body integrable quantum chains, including XXZ-type spin models and interacting lattice fermions with dissipative two-body loss (Hamanaka et al., 2023, Mao et al., 2024). In these cases, the LSE is established analytically via Bethe ansatz, with eigenfunctions displaying boundary-pinned exponential profiles. The LSE's fragility and persistence are manifest under generalized boundary conditions.

In correlated systems, complex interactions open point-gaps in the Liouvillian spectrum, endowing it with a topological invariant and producing interaction-induced skin modes. These phenomena manifest as transient, boundary-accumulated local densities and are absent in the periodic setting where topological invariants vanish (Hamanaka et al., 2023).

Non-Markovian generalizations using hierarchical equations of motion reveal that memory effects can "thicken" the skin modes and generate robust cross-site coherence, leading to novel relaxation oscillations while preserving the characteristic $\mathcal{L}$4 scaling (Kuo et al., 2024).

6. Connections to Broader Non-Hermitian Phenomena and Topology

The LSE generalizes the NHSE from single-particle Hamiltonians to full Liouville dynamics of open systems, incorporating the impact of quantum jumps. It is distinguished from the NHSE by the structure of open quantum dynamics, many-body interplay, and richer topology:

• LSE supports both $\mathcal{L}$5 (winding number) and $\mathcal{L}$6 (Kramers-paired) topological protection, as seen in two-dimensional electron systems with spin-orbit coupling and magnetic fields (Shigedomi et al., 23 May 2025).
• LSE can enable critical delocalization transitions—i.e., a "critical LSE"—where even an infinitesimal coupling between subsystems radically alters the spectral scaling and spatial localization, shifting from boundary-pinned to extended steady states (Feng et al., 2024).

Further, the spatial mode structure of the LSE provides a pathway for the design and engineering of protocol-dependent relaxation phenomena, such as the quantum Mpemba and Pontus–Mpemba effects, where the overlap of initial conditions with skin modes enables anomalously fast or engineered relaxation behaviors (Longhi, 20 Jan 2026, Sun et al., 22 Jan 2026).

7. Summary Table: Key Formulas and Scaling Relations

Phenomenon / Parameter	Formula / Scaling	Reference
Liouvillian gap	$\mathcal{L}$7	(Haga et al., 2020)
Skin mode localization	$\mathcal{L}$8	(Haga et al., 2020)
Longest relaxation time	$\mathcal{L}$9	(Haga et al., 2020)
Topological invariant (winding)	$D^2$0	(Hamanaka et al., 2023)
Critical LSE gap scaling	$D^2$1; delocalized steady state	(Feng et al., 2024)
LSE relaxation (gradient)	$D^2$2 (anomalous, fast)	(Wang et al., 2023)

The Liouvillian skin effect is thus a robust mechanism by which non-Hermitian topology and irreversible dynamics produce boundary-localized eigenmodes, fundamentally modifying relaxation, entanglement, and transport behavior in open quantum many-body systems. It provides both a diagnostic and a functional tool for controlling and understanding quantum dissipative architectures, beyond conventional bulk spectral analyses.