Liouville Skin Effect in Open Quantum Systems

• Liouville skin effect is the boundary localization of eigenmodes of the Liouvillian superoperator governing open-system dynamics.
• It reveals how non-normal spectral geometry and asymmetric mode localization induce distinct relaxation behaviors under open boundary conditions.
• Models like dissipative SSH chains and asymmetric spin exchange illustrate the practical impact on relaxation times and steady-state properties.

Liouville skin effect, more commonly termed the Liouvillian skin effect (LSE), is the boundary localization of eigenmodes of the Liouvillian superoperator that governs open-system density-matrix dynamics. In contrast to the non-Hermitian skin effect of Hamiltonian eigenstates, the relevant objects are operator-valued modes in Liouville space, and the primary consequences appear in relaxation, transient response, and approach to the steady state rather than in stationary wave mechanics alone (Longhi, 20 Jan 2026, Haga et al., 2020).

1. Definition and conceptual scope

In the standard formulation, an open quantum system evolves according to a Lindblad master equation,

$\dot\rho=\mathcal{L}[\rho],$

and the LSE refers to the situation in which right and left eigenmodes of $\mathcal{L}$ become exponentially localized at opposite boundaries under open boundary conditions. This is the direct Liouville-space analogue of the non-Hermitian skin effect, but with three distinctions emphasized repeatedly in the literature: the generator is a non-Hermitian superoperator rather than a Hamiltonian, the eigenmodes are operators rather than wavefunctions, and the physical consequences are encoded in dissipative dynamics with a steady state rather than in energy spectra alone (Longhi, 20 Jan 2026).

This distinction is especially important because localization now refers to matrix elements of eigenoperators, typically concentrated near specific spatial regions in the doubled space of indices $(n,m)$. In representative one-dimensional chains with non-reciprocal dissipation, right eigenmodes localize near one edge while the dual left eigenmodes localize near the opposite edge, producing strong biorthogonal asymmetry and boundary sensitivity (Longhi, 20 Jan 2026).

A distinct operational usage also exists. In quadratic Lindbladian SSH chains with uniform local coupling to the environment, the steady state remains spatially uniform and the Dirac-density profiles of individual normal modes show no static boundary localization. In that setting, the “Liouvillian skin effect” is identified with chiral damping: a highly asymmetric boundary-propagating relaxation pattern generated by interference among many Liouvillian normal modes rather than by boundary-localized single-mode densities (Zhou et al., 2021). This definitional split is one of the central conceptual nuances in the subject.

2. Liouville-space spectral framework

The spectral description uses right and left Liouvillian eigenmodes,

$$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$

with biorthogonality

$(L_\alpha|R_\beta)=\delta_{\alpha\beta}.$

For a stable Markovian problem, $\lambda_1=0$ labels the steady state, while all other eigenvalues have negative real parts and determine decay toward stationarity. Expanding an initial state in this basis,

$$\rho_i=\rho_E+\sum_{\alpha>1} c_\alpha R_\alpha,\qquad c_\alpha=\mathrm{Tr}(L_\alpha^\dagger \rho_i),$$

gives

$$\rho(t)=\rho_E+\sum_{\alpha>1} c_\alpha e^{\lambda_\alpha t}R_\alpha,$$

so relaxation depends not only on spectral gaps but also on overlap coefficients fixed by the geometry of left and right eigenmodes (Longhi, 20 Jan 2026).

The non-normal character of the Liouvillian is therefore central. When

$[\mathcal{L},\mathcal{L}^\dagger]\neq 0,$

right eigenvectors are not mutually orthogonal, left eigenvectors are not adjoints of the right ones, and the conditioning of the eigenbasis can become poor. In LSE regimes this non-normality is amplified by opposite-edge localization of left and right modes, so two initial states with the same spectral gap can exhibit radically different effective relaxation because their projections onto the slow mode differ exponentially (Longhi, 20 Jan 2026).

For quadratic fermionic Lindbladians, this structure can often be made explicit by third quantization. In the exactly solvable dissipative SSH chain, the Liouvillian in the even-parity sector is written as

$$\hat{\mathcal{L}}_+ = -\sum_{m=1}^{2n}\beta_m\, b'_m b_m,$$

where the rapidities $\mathcal{L}$0 control damping through $\mathcal{L}$1, and the Liouvillian gap is

$\mathcal{L}$2

In that model the rapidities are obtained from a non-Hermitian SSH damping matrix, making the Liouville-space skin effect exactly solvable (Yang et al., 2022).

3. Canonical mechanisms and exactly solvable realizations

A minimal mechanism is non-reciprocal dissipation. In a bosonic tight-binding chain with coherent hopping $\mathcal{L}$3, asymmetric incoherent hoppings $\mathcal{L}$4, and tunable end-to-end coherent coupling $\mathcal{L}$5, the open-chain limit $\mathcal{L}$6 exhibits a Liouvillian skin effect. In the classical limit $\mathcal{L}$7, the stationary distribution is

$\mathcal{L}$8

and the localization length scales as

$\mathcal{L}$9

For $(n,m)$0, right Liouvillian eigenvectors localize near the right edge and left eigenvectors near the left edge; for $(n,m)$1, the roles reverse (Longhi, 20 Jan 2026).

A second canonical realization is the exactly solvable bond-dissipative SSH chain. There the damping matrix factorizes as

$(n,m)$2

with $(n,m)$3 a non-Hermitian SSH Hamiltonian with asymmetric hoppings, so that the rapidities are

$(n,m)$4

Periodic and open boundary spectra are related by a complex momentum shift, and the normal master modes inherit the same sensitivity to boundary conditions that characterizes the non-Hermitian SSH skin effect (Yang et al., 2022).

An exactly solvable many-body construction arises from asymmetric dissipative spin exchange. In the diagonal Liouville-space sector, the effective Liouvillian reduces to a non-Hermitian XXZ chain with asymmetric hopping. Under open boundary conditions, an imaginary gauge transformation maps the bulk problem to a Hermitian XXZ chain while leaving an exponential weight in the wavefunction, producing a skin-localized steady state; in the large-$(n,m)$5 limit this approaches a domain-wall configuration with all up-spins packed at one edge (Mao et al., 2024).

These models show that the microscopic origin of the LSE need not be unique. The effect can emerge from asymmetric incoherent hopping, bond-dissipative lattice fermion models, or integrable dissipative spin chains, provided the Liouvillian develops a non-Hermitian spectral geometry with boundary-sensitive eigenmodes.

4. Relaxation geometry, non-normality, and protocol dependence

A foundational result is that Liouvillian skin localization changes the relation between the Liouvillian gap and the actual longest relaxation time. In the prototypical asymmetric-hopping bosonic chain, the maximal relaxation time over initial states and local observables is

$(n,m)$6

with $(n,m)$7 the system size, $(n,m)$8 the Liouvillian gap, and $(n,m)$9 the localization length. The consequence is that $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$0 can diverge as $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$1 without any gap closing (Haga et al., 2020).

The mechanism is geometric rather than purely spectral. Because the slowest right and left eigenmodes are localized near opposite edges, their overlap is exponentially small. The coefficient of the slow mode in the time evolution can therefore become exponentially large or exponentially small depending on the initial condition, even though the asymptotic decay rate $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$2 is unchanged. This is the sense in which LSE makes relaxation prefactor dominated rather than gap dominated (Haga et al., 2020).

This non-normal geometry can also be exploited constructively. In the quantum Pontus–Mpemba protocol, the same initial state $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$3 is relaxed toward the same stationary state by two routes: a direct protocol with $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$4, and a two-step protocol in which a short coherent boundary coupling first swaps the excitation approximately to $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$5, after which the same Liouvillian is restored. The two protocols share the same asymptotic decay rate set by $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$6, but the slow-mode overlap changes by

$$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$7

so the two-step protocol can reach a practical threshold in trace distance or Hilbert–Schmidt distance much earlier. The effect disappears at $$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$8, exactly when the skin effect vanishes (Longhi, 20 Jan 2026).

5. Boundary conditions, disorder, and generalized reciprocity

Sensitivity to boundary conditions is one of the defining signatures of the LSE. In the exactly solvable dissipative XXZ chain, generalized boundary conditions distinguish counter-flow from co-flow edge hopping. Counter-flow hopping is exponentially suppressed after the imaginary-gauge transformation and leaves the LSE intact in the thermodynamic limit, whereas any finite co-flow boundary hopping is exponentially amplified and destroys the skin effect, making the spectrum and eigenfunctions effectively periodic-boundary-like at large system size (Mao et al., 2024).

Disorder introduces a different regime. In a model with globally reciprocal Liouvillian dynamics and locally asymmetric incoherent hopping,

$$\mathcal{L}|R_\alpha)=\lambda_\alpha |R_\alpha),\qquad (L_\alpha|\mathcal{L}=\lambda_\alpha (L_\alpha|,$$9

the steady state no longer exhibits boundary accumulation. Instead, it shows disorder-dependent, erratic localization at realization-specific bulk positions, and in the incoherent-hopping-dominated regime the transport is not ballistic or diffusive but Sinai-type subdiffusive, with spreading dramatically slower than ordinary diffusion in symmetric stochastic lattices (Longhi, 16 Feb 2026). This establishes a sharp contrast with globally reciprocal Hamiltonian problems, where the suppression of the conventional skin effect can coexist with ballistic transport.

Another route away from conventional LSE is interchain coupling of opposite skin directions. In a two-chain Lindbladian whose uncoupled legs localize toward opposite edges, increasing interchain coupling gradually delocalizes the skin steady state. In the single-particle sector the Liouvillian gap scaling crosses from

$(L_\alpha|R_\beta)=\delta_{\alpha\beta}.$0

to

$(L_\alpha|R_\beta)=\delta_{\alpha\beta}.$1

and in the many-body problem the steady-state entanglement scaling changes from an area law to a logarithmic law. In the thermodynamic limit, even arbitrarily small coupling can induce dramatic spectral and steady-state changes; this singular boundary sensitivity is termed the critical Liouvillian skin effect (Feng et al., 2024).

6. Extensions, platforms, and current directions

The Markovian Lindblad setting is not the endpoint of the subject. In a tilted chain coupled collectively to a bosonic bath and treated with hierarchical equations of motion, non-Markovian memory produces a thick skin effect: boundary-localized modes persist but broaden into the bulk. The same analysis identifies coherent skin modes with complex eigenvalues, oscillatory relaxation, and cross-site coherence that appears only when counter-rotating terms are retained. In that regime, both the dominant skin mode and the steady-state coherence remain robust against additional dephasing noise (Kuo et al., 2024).

In magnetic systems, a microscopic spin chain coupled to a shared bosonic reservoir clarifies when the Liouvillian viewpoint reduces to an effective non-Hermitian Hamiltonian and when it does not. In the dilute-magnon, zero-temperature limit, the Liouvillian spectrum is completely determined by a non-Hermitian magnon Hamiltonian, so NHSE and LSE coincide. Beyond that limit, however, the effective non-Hermitian Hamiltonian fails to capture relaxation times, pumping effects, and the conditions under which boundary accumulation is experimentally visible before decay erases it (Li et al., 2023).

Finite-temperature electronic realizations have also been identified. A two-dimensional electron system on a substrate with spin–orbit coupling, transverse magnetic field, and energy dissipation exhibits both $(L_\alpha|R_\beta)=\delta_{\alpha\beta}.$2 and $(L_\alpha|R_\beta)=\delta_{\alpha\beta}.$3 Liouvillian skin effects. Their temperature dependence is controlled by the band-splitting scale induced by spin–orbit coupling and magnetic field, and the $(L_\alpha|R_\beta)=\delta_{\alpha\beta}.$4 case yields charge accumulation under quench dynamics. In contrast to previously reported LSEs with system-size-dependent slowdown, the relaxation time there is independent of system size because the localization length is scale-free, analogous to non-Hermitian critical skin effects (Shigedomi et al., 23 May 2025).

Taken together, these developments show that “Liouville skin effect” now denotes a family of boundary-sensitive non-Hermitian phenomena in open quantum systems. In its strict form it is the exponential localization of Liouvillian eigenmodes; in broader operational usage it can also denote chiral, interference-driven transient response in Liouville space. Across both usages, the recurring themes are non-normal spectral geometry, strong boundary sensitivity, and a breakdown of the expectation that the Liouvillian gap alone determines dissipative timescales.