Directional dissipation can concentrate relaxation modes in space, but connecting this accumulation to controlled microscopic processes in correlated matter requires separating virtual charge motion from spin coherence. We study this connection in a half-filled Hubbard chain with reciprocal hopping and number-conserving, direction-selective returns of virtual doublon--hole defects. Eliminating charge defects and then Mott coherences yields asymmetric spin-exchange rates under an explicit separation of timescales, with full short-chain dynamics supporting the reduction on finite exchange times. Balanced random returns concentrate the stationary population and several low-lying right modes at sample-selected positions rather than a predetermined edge, producing random Liouvillian skin localization despite zero end-to-end logarithmic rate bias. The accumulated rate imbalance determines an exact finite-density hard-core stationary state and a one-particle activated relaxation scale governed by random barriers. This connection between locally calibrated return rates, many-body stationary weights and slow relaxation offers a means of controlling transport in constrained open quantum matter without changing the reciprocal Hamiltonian.
Directional dissipation can concentrate relaxation modes in space, but connecting this accumulation to controlled microscopic processes in correlated matter requires separating virtual charge motion from spin coherence. We study this connection in a half-filled Hubbard chain with reciprocal hopping and number-conserving, direction-selective returns of virtual doublon--hole defects. Eliminating charge defects and then Mott coherences yields asymmetric spin-exchange rates under an explicit separation of timescales, with full short-chain dynamics supporting the reduction on finite exchange times. Balanced random returns concentrate the stationary population and several low-lying right modes at sample-selected positions rather than a predetermined edge, producing random Liouvillian skin localization despite zero end-to-end logarithmic rate bias. The accumulated rate imbalance determines an exact finite-density hard-core stationary state and a one-particle activated relaxation scale governed by random barriers. This connection between locally calibrated return rates, many-body stationary weights and slow relaxation offers a means of controlling transport in constrained open quantum matter without changing the reciprocal Hamiltonian.
We unveil a quantum Pontus-Mpemba effect enabled by the Liouvillian skin effect in a dissipative tight-binding chain with asymmetric incoherent hopping and coherent boundary coupling. The skin effect, induced by non-reciprocal dissipation, localizes relaxation modes near the system boundaries and gives rise to non-orthogonal spectral geometry. While such non-normality is often linked to slow relaxation, we show that it can instead accelerate relaxation through a two-step protocol - realizing a quantum Pontus-Mpemba effect. Specifically, we consider a one-dimensional open chain with coherent hopping J, asymmetric incoherent hoppings JR=JL, and a controllable end-to-end coupling ε. For ε=0, the system exhibits the Liouvillian skin effect, with left and right eigenmodes localized at opposite edges. We compare two relaxation protocols toward the same stationary state: (i) a direct relaxation with ε=0, and (ii) a two-step (Pontus) protocol where a brief coherent evolution transfers the excitation across the lattice before relaxation. Although both share the same asymptotic decay rate, the two-step protocol relaxes significantly faster due to its reduced overlap with the slow boundary-localized Liouvillian mode. The effect disappears when JR=JL, i.e., when the skin effect vanishes. Our results reveal a clear connection between boundary-induced non-normality and protocol-dependent relaxation acceleration, suggesting new routes for controlling dissipation and transient dynamics in open quantum systems.
We develop a microscopic theory for how slow Liouvillian sectors are selected in an open correlated lattice. The starting point is not a postulated non-Hermitian band, but a local interacting resonance between an on-site doublon and a branch-resolved nearest-neighbor bond. This resonance defines a composite shell orbital whose doublon weight controls reservoir visibility and whose mixed doublon-bond character controls shell mobility. Projecting the microscopic hopping onto the selected shell yields a branch-selective dimerized channel. In the dilute regime, a boundary doublon-loss channel yields an exponentially slow edge-memory pole through a Zeno-type return. At the shell-critical point, the edge pole is replaced by a near-zero standing-wave doublet with an algebraic coherent spacing. At finite shell filling, the same local shell becomes density dressed. A number-conserving phase-locking jump removes a bright mismatch sector, leaving defects as the asymptotic slow variables and producing a diffusive finite-size gap. We derive the local shell, the projected branch topology, the edge-memory law, the shell-critical doublet, the density-dressed shell Hamiltonian, and the defect generator within one Schur-projection framework. The resulting mechanism identifies the reservoir-engineered fast block as the selector of the observable slow sector, while the microscopic parent shell remains fixed.
Liouvillian skin effects are usually associated with dissipation whose rates favor unidirectional flow. In this work, we show that directional edge localization can emerge even when population exchange between sites is exactly balanced, with a single tunable phase parameter α deciding their direction and a dual winding structure deciding their nature. Although this number-conserving bath is known for eigenstate selection, we demonstrate here that its phase controls a steady current, which Hamiltonian hopping converts to directional edge localization. Further, we study the spectral winding structure by independently twisting the left (ket) and right (bra) vector indices of the density matrix, revealing two distinct windings whose relationship is dictated by an antiunitary mirror symmetry. Nearest-neighbor jump operators break this symmetry: both windings coincide, giving rise to a common spectral winding W0(α)=±1, producing a unipolar skin effect in which the density accumulates at a specific edge. However, next nearest neighbor jump preserve it. The two windings are forced to be opposite, and the ket and bra are pushed to opposite edges, giving a bipolar skin effect of coherences present in eigenmodes, without edge localization of the density in the steady state. Remarkably, in a dimerized chain the phase-controlled density skin effect gives the steady-state zero modes a size-independent weight at the skin edge, accompanied by a kink at the topological transition. In addition, the relaxation time grows linearly with the system size, however at specific phases α where the jump operators are rendered Hermitian, edge-localized eigenmodes with skin character coexist with diffusive relaxation.
Non-Hermitian systems with globally reciprocal couplings -- such as the Hatano-Nelson model with stochastic imaginary gauge fields -- avoid the conventional non-Hermitian skin effect, displaying erratic bulk localization while retaining ballistic transport. An open question is whether similar behavior arises when non-reciprocity originates at the Liouvillian level rather than from an effective non-Hermitian Hamiltonian obtained via post-selection. Here, a lattice model with globally reciprocal Liouvillian dynamics and locally asymmetric incoherent hopping is investigated, a disordered setting in which Liouvillian-specific effects have remained largely unexplored. While the steady state again shows disorder-dependent, erratic localization without boundary accumulation, {\color{black}excitations in the incoherent-hopping regime spread via {\em Sinai-type subdiffusion}, dramatically slower than ordinary diffusion in symmetric stochastic lattices.} {\color{black}This highlights that the genuinely distinct Liouvillian signature is the coexistence of global reciprocity with ultra-slow, disorder-induced subdiffusive transport, rather than the erratic localization itself.} {\color{black}These results reveal a fundamental distinction between globally reciprocal Hamiltonian and Liouvillian systems: in both cases the skin effect is suppressed, but only in Liouvillian dynamics erratic skin localization can coexist with subdiffusive transport
Reservoir engineering can stabilize states inaccessible to unitary dynamics. Directed particle-conserving dissipation creates Lindbladian skin states, where Pauli exclusion turns edge accumulation into a many-body density imbalance. In a monitored fermion chain with tunable hopping range, we identify, within a Gaussian trajectory approximation, two finite-size scaling regimes: short-range hopping is consistent with complete skin accumulation and area-law entanglement, whereas sufficiently long-range hopping produces a finite bulk tail and effective algebraic sub-volume-law entanglement. Dissipation and coherent hopping thus jointly control skin localization and quantum entanglement, highlighting their close interconnection.
Nonreciprocal relaxation matrices can have skin-localized right eigenmodes, but their imprint on a mixed steady state is not fixed by the density profile alone. We develop an exact steady-state theory for number-conserving Gaussian fermion chains and show that the dominant natural orbital of the correlation matrix provides a mode-resolved diagnostic of hidden skin order. The steady-state correlator admits a biorthogonal decomposition in terms of the left and right eigenmodes of the relaxation matrix X and the source matrix Y. This formula separates three ingredients: slow rapidity denominators, source loading by left eigenmodes, and real-space geometry from right eigenmodes. For a local pump, the pump position is read by the left modes, whereas the selected profile is drawn by the right modes. In a single-slow-mode regime, the dominant natural orbital locks to the Euclidean-normalized slow right mode. The density can follow the same boundary trend, but it is a less selective incoherent sum over occupied natural orbitals. We verify this selection law in a nonreciprocal Hatano--Nelson chain and show that, in a nonreciprocal SSH chain, the selected natural orbital crosses over from a topological edge candidate to a slow bulk-skin candidate. These results identify natural-orbital locking as a steady-state diagnostic of nonreciprocal localization in Gaussian open fermion chains.