Spectrally Separable Lindbladians

• Spectrally separable Lindbladians are GKLS generators whose Liouvillian spectra decompose into distinct, invariant sectors via mechanisms like dissipative gaps, block-triangular structures, and direct-integral decompositions.
• The separation between peripheral and decaying spectra ensures precise asymptotic dynamics and rapid convergence to steady states, as well as robust Poisson eigenvalue statistics in certain regimes.
• Advanced techniques such as charge-ordered triangularity, operator-space fragmentation, and additive rapidity constructions provide practical tools for analyzing and exactly solving quadratic dissipative models.

Spectrally separable Lindbladians are Gorini–Kossakowski–Sudarshan–Lindblad generators whose spectra admit an exact decomposition into simpler spectral pieces. In the literature, this phrase is used in several closely related but non-identical senses: separation of peripheral and strictly decaying spectrum by a dissipative gap; block-triangular Liouville-space structure whose diagonal blocks determine the full eigenvalue set; direct-integral decomposition over quasi-momentum; decomposition into invariant operator-space sectors; separation of inter-orbit and intra-orbit dynamics under fast control; and additive many-body spectra built from single-mode rapidities (Albert, 2018, Wang et al., 19 Mar 2026, Klausen et al., 2022, Essler et al., 2020, Rooney et al., 2016, Honda et al., 2010, Wang et al., 2024). In all of these usages, the common theme is that the Liouvillian eigenproblem becomes reducible to lower-dimensional, fiberwise, sectorwise, or modewise data.

1. GKLS framework and the range of meanings

A bounded Lindbladian in standard form acts on density matrices as

$$\frac{d\rho}{dt}=\mathcal{L}(\rho)=-i[H,\rho]+\sum_k\Big(L_k\rho L_k^\dagger-\tfrac12\{L_k^\dagger L_k,\rho\}\Big),$$

with adjoint

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$

The dynamics preserves complete positivity, trace, and Hermiticity, and the Lindbladian spectrum satisfies $\operatorname{Re}(\lambda)\le 0$, with complex eigenvalues occurring in conjugate pairs (Albert, 2018).

Within this common framework, separability is realized in different ways. One line of work defines it by the existence of an asymptotic projection onto a steady-state or peripheral manifold separated from the decaying spectrum by a nonzero dissipative gap (Albert, 2018). Another defines it by a charge-ordered block-triangular Liouville matrix in which recycling terms are strictly off-diagonal, so that the full spectrum is exactly the spectrum of the no-jump effective Hamiltonian superoperator (Wang et al., 19 Mar 2026). In translation-invariant single-particle problems, separability means a direct integral

$$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$

with finite-range Laurent fibers (Klausen et al., 2022). In fragmented one-dimensional many-body models, it means a block decomposition

$$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$

into exponentially many invariant operator sectors (Essler et al., 2020). Under fast Hamiltonian controllability, it refers instead to the separation of eigenvalue dynamics from flag dynamics on unitary orbits (Rooney et al., 2016). In quadratic exactly solvable models, it becomes additive factorization of Liouvillian eigenvalues into dissipation and phase or rapidity quantum numbers (Honda et al., 2010, Wang et al., 2024).

The coexistence of these usages is a genuine feature of the subject rather than a notational accident. The term consistently signals that the Liouvillian spectral problem is reducible, but the structural source of reducibility depends on the physical setting.

2. Peripheral spectrum, dissipative gap, and the steady-state manifold

For Lindbladians with multiple steady states, the central object is the asymptotic subspace

$$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$

where $P$ is idempotent, $P^2=P$. The asymptotic evolution has the form

$$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$

with $H_\infty$ the residual peripheral Hamiltonian. When $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$0, $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$1 is a stationary steady state (Albert, 2018).

In this setting, spectral separability means a decomposition of the spectrum into peripheral eigenvalues on the imaginary axis and strictly decaying eigenvalues with negative real part, separated by a nonzero dissipative gap

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$2

The gap controls exponential convergence to the asymptotic manifold. The paper further states that all Jordan blocks of $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$3 with pure imaginary eigenvalues are diagonal, so the peripheral part is semisimple and supports only blockwise unitary motion (Albert, 2018).

The operator space admits a four-corners decomposition using the orthogonal projection $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$4 onto the maximal non-decaying subspace of $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$5 and $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$6:

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$7

Population in the decaying block vanishes asymptotically,

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$8

and the asymptotic subspace is block-diagonal in this decomposition. In important cases with multiple steady states and genuine decay, the asymptotic projection is explicitly

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$9

where $\operatorname{Re}(\lambda)\le 0$0 is the inverse restricted to the decaying block. This formula makes precise how initial weight in $\operatorname{Re}(\lambda)\le 0$1 feeds into the non-decaying sector while coherences in $\operatorname{Re}(\lambda)\le 0$2 decay away (Albert, 2018).

A second structural statement is the duality between Schrödinger and Heisenberg fixed points. If $\operatorname{Re}(\lambda)\le 0$3 spans the asymptotic manifold and $\operatorname{Re}(\lambda)\le 0$4 satisfies $\operatorname{Re}(\lambda)\le 0$5, then the number of steady-state basis elements equals the number of linearly independent conserved quantities and one can normalize them so that $\operatorname{Re}(\lambda)\le 0$6. The same framework includes decoherence-free subspaces and noiseless subsystems. The paper emphasizes, however, that a Noether-type implication fails: symmetry block decomposition does not guarantee conserved quantities in every block, so $\operatorname{Re}(\lambda)\le 0$7 and $\operatorname{Re}(\lambda)\le 0$8 must be determined explicitly (Albert, 2018).

3. Charge-ordered triangularity, recycling, and spectral statistics

A different notion of spectral separability arises from the no-jump/recycling split

$\operatorname{Re}(\lambda)\le 0$9

with

$$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$0

After vectorization,

$$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$1

so the no-jump part acts independently on ket and bra legs while recycling couples them rung by rung (Wang et al., 19 Mar 2026).

The defining conditions are structural and algebraic. There must exist a physical $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$2 charge $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$3 with $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$4, so the no-jump Liouville superoperator commutes with the total Liouville charge $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$5. In addition, each jump operator must change the physical charge by a fixed sign, so the recycling term shifts $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$6 only in one direction, for example $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$7 for lowering jumps. In a basis ordered by $$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$8, the full Lindbladian is then block upper- or lower-triangular, with diagonal blocks given by the charge-resolved no-jump superoperator and strictly off-diagonal recycling blocks (Wang et al., 19 Mar 2026).

The immediate consequence is

$$\mathcal{J}^\uparrow(\mathcal{L})=\int_{-\pi}^{\pi}{}^\oplus M(k)\,dk, \qquad \sigma(\mathcal{L})=\bigcup_k \sigma(M(k)),$$9

with eigenvalues

$$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$0

where $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$1 are eigenvalues of $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$2. Recycling mixes eigenvectors across charge sectors but does not change the eigenvalue set. The paper explicitly stresses that this is not a symmetry decomposition of $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$3: the spectrum is separable, while eigenvectors remain entangled across sectors (Wang et al., 19 Mar 2026).

This structure yields robust Poisson eigenvalue statistics for the full Lindbladian. The stated mechanism is the decoupled-difference structure of $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$4 in the complex plane together with the fact that strictly block-off-diagonal recycling does not generate additional level repulsion. In uniform damping, with $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$5 and $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$6, the eigenvalues take the banded form

$$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$7

so real parts are equally spaced vertical shifts indexed by total Liouville magnetization $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$8. The paper reports that spectrally separable Lindbladians then show Poisson statistics even when the associated $$O=\bigoplus_\alpha O_\alpha, \qquad \mathrm{spec}(\mathcal{L})=\bigcup_\alpha \mathrm{spec}(\mathcal{L}|_{O_\alpha}),$$9 has chaotic spectral correlations (Wang et al., 19 Mar 2026).

The limits of this mechanism are equally sharp. Dephasing jumps $$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$0 conserve charge but do not act one-sidedly on $$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$1; they induce inter-leg couplings within the same charge sector and can restore chaos in interacting models. Mixed gains and losses or boundary driving likewise destroy triangularity by making the recycling action two-sided (Wang et al., 19 Mar 2026).

4. Direct-integral separability in translation-invariant lattice models

For single-particle translation-invariant Lindbladians on $$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$2, spectral separability is realized as a momentum-space direct integral. The Hilbert space is $$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$3, the Hamiltonian is the tight-binding Laplacian

$$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$4

and the jump operators are finite-range translates of a seed $$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$5,

$$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$6

Translation covariance then implies

$$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$7

On Hilbert–Schmidt space, vectorization, a conditional shift, and a partial Fourier transform yield an isometric isomorphism

$$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$8

under which

$$\mathrm{As}(\mathcal{H})\equiv P\,\mathrm{Op}(\mathcal{H}), \qquad P\equiv \lim_{t\to\infty} e^{t\mathcal{L}},$$9

with Laurent fiber $P$0 and finite-rank jump term $P$1 (Klausen et al., 2022).

Norm continuity in $P$2 gives the spectrum identity

$P$3

This is the defining notion of separability in that work. The Laurent part has symbol calculus, so $P$4 is the range of its symbol and is purely essential; finite-rank fiber perturbations can only add point spectrum. For rank-one fibers,

$P$5

reducing jump-induced eigenvalues to a scalar resolvent equation (Klausen et al., 2022).

The paper proves further spectral properties under the standing assumptions. The spectrum is purely approximate point, so there is no residual spectrum:

$P$6

If $P$7, then either the generator is gapless at $P$8 or $P$9 is infinite-dimensional. For finite-volume periodic approximants, one obtains exact block diagonalization over discrete momenta and Hausdorff convergence of periodic finite-volume spectra to the infinite-volume spectrum; periodic boundary conditions are essential because free boundaries can show strong deviations associated with the non-Hermitian skin effect (Klausen et al., 2022).

Canonical examples are worked out explicitly. For local dephasing, $P^2=P$0 with coupling $P^2=P$1,

$P^2=P$2

and the full spectrum is

$P^2=P$3

This model is gapless at $P^2=P$4 (Klausen et al., 2022).

5. Operator-space fragmentation and integrable sector decompositions

In one-dimensional many-body systems, spectral separability can arise from operator-space fragmentation. Here the operator space $P^2=P$5 decomposes as

$P^2=P$6

or equivalently, after vectorization, $P^2=P$7 with invariant subspaces $P^2=P$8. The mechanism is a complete set of mutually commuting local projectors on $P^2=P$9 that commute with $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$0 (Essler et al., 2020).

For the spin-$$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$1 ASEP family,

$$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$2

the Liouvillian commutes with onsite rank-one projectors $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$3 and $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$4, while $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$5 projects onto the two-dimensional local diagonal span. The number of invariant sectors is therefore $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$6. A sector with defects at sites $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$7 decomposes into independent segments between neighboring defects, and on long enough segments the projected Liouvillian is similar to an open XXZ chain. In the defect-free diagonal sector, the projected Liouvillian is the classical ASEP generator

$$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$8

which is integrable. A local similarity transform maps it to an XXZ chain with a diagonal twist at the boundary (Essler et al., 2020).

The sectorwise spectral content is correspondingly explicit. In the diagonal sector, the ground-state energy is $$\rho(t)\xrightarrow{t\to\infty}\rho_\infty = e^{-iH_\infty t}\,P(\rho(0))\,e^{iH_\infty t},$$9 with degeneracy $H_\infty$0, given by equal-weight states over fixed magnetization sectors. In fully fragmented sectors, any product state is an eigenstate; two special product states, $H_\infty$1 and $H_\infty$2, are annihilated by $H_\infty$3. The total ground-state degeneracy is therefore $H_\infty$4. The Liouvillian gap in the diagonal sector scales as $H_\infty$5, while on any open-segment defect sector one has the bound

$H_\infty$6

Thus the slow modes are concentrated in the diagonal ASEP sector, and late-time dynamics is dominated by that block (Essler et al., 2020).

The construction extends to arbitrary local dimension. For $H_\infty$7-level systems with jumps

$H_\infty$8

the number of sectors becomes $H_\infty$9. In the fully diagonal sector the projected Liouvillian reduces to the SU($$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$00)-invariant Sutherland chain, while generic defect sectors yield open SU($$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$01) chains. The significance of this framework is that separability is not restricted to steady subspaces: it furnishes a complete block decomposition of the full dissipative dynamics into integrable pieces (Essler et al., 2020).

6. Fast-control separation of eigenvalue and flag dynamics

A control-theoretic notion of spectral separability appears in finite-dimensional Lindblad systems with fast Hamiltonian controllability. Writing

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$02

with unbounded piecewise-continuous controls, one may make intra-orbit unitary motion much faster than the dissipative time scale set by the bounded Lindblad superoperator. The density operator is decomposed as

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$03

where the eigenvalue data describe the inter-orbit motion and the projectors, or complete flag $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$04 when eigenvalues are simple, describe the intra-orbit motion (Rooney et al., 2016).

The spectral drift obeys a classical rate equation controlled by the flag:

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$05

with

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$06

Thus $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$07 is a classical Markov generator whose coefficients depend on the chosen flag. After projection to the trace-free simplex coordinates $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$08, the dynamics becomes

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$09

Given a spectral/flag trajectory, one can reconstruct a Hamiltonian $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$10 that implements the required flag motion and corrects for dissipative couplings between eigenspaces (Rooney et al., 2016).

In this usage, spectral separability does not mean that the spectrum of $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$11 decomposes into invariant blocks independent of control. It means that inter-orbit spectral evolution and intra-orbit flag evolution can be treated separately, with the flag acting as an effective control parameter. The paper develops local controllability criteria in this reduced picture. Strong local controllability holds if

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$12

At the completely mixed state $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$13, the accessible velocities are the convex hull of the $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$14 permutations of the eigenvalues of

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$15

These $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$16 permutation flags provide a natural finite control set and are described as highly effective for low-purity orbits near $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$17 (Rooney et al., 2016).

This framework clarifies a frequent ambiguity in the term. Here the separated objects are not spectral subspaces of the Liouvillian itself, but the slow eigenvalue dynamics of the state and the fast flag dynamics on unitary orbits.

7. Additive rapidities, exact solvability, and correlator-based reconstruction

The most explicit additive realization of spectral separability appears in exactly solvable quadratic models. For the damped harmonic oscillator with Hamiltonian $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$18 and jump operators $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$19, $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$20, the full Liouvillian spectrum is

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$21

The real part is controlled by the non-negative dissipation quantum number $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$22 and the imaginary part by the phase-rotation quantum number $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$23. Ladder superoperators $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$24 and $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$25 generate the spectral towers from the thermal stationary state, with $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$26 raising both $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$27 and $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$28 by one and $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$29 raising $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$30 while lowering $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$31 by one (Honda et al., 2010).

This additive organization generalizes to quadratic fermionic Liouvillians with quadratic Hamiltonian

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$32

linear jump operators

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$33

and diagonalizable normal-form block $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$34. The Liouvillian can be written in terms of Majorana superoperators and brought to canonical form

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$35

with parity-independent rapidities $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$36. The spectrum then becomes

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$37

so even and odd parity sectors consist of all such sums with even or odd occupation parity, respectively. This is spectral separability in a strict additive sense: many-body eigenvalues are linear sums of single-mode rapidities (Wang et al., 2024).

The same work derives closed recursion equations for all even-order correlation functions,

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$38

together with a matrix Wick theorem for Gaussian initial states. After projection onto fully antisymmetric index combinations, the reduced blocks $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$39 have eigenvalues

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$40

reproducing the additive spectral construction from correlator dynamics (Wang et al., 2024).

The additive picture has a sharp boundary. For quartic Liouvillians with quadratic dissipation, closure of the two-point sector holds if and only if

$$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$41

Under these conditions, one still obtains block-triangular recursion for even-order correlators and constructs the even-parity spectrum as the union of eigenvalues of reduced blocks $$\mathcal{L}^\dagger(X)=-i[H,X]+\sum_k\Big(L_k^\dagger X L_k-\tfrac12\{L_k^\dagger L_k,X\}\Big).$$42. The paper states, however, that strict additivity into sums of single-mode rapidities generally fails once the extra quadratic-dissipation couplings are present. In this sense, blockwise spectral construction survives, but the stronger form of spectral separability does not (Wang et al., 2024).

Taken together, these exactly solvable models show the strongest available form of the concept: the Liouvillian spectrum is not merely reducible sector by sector, but can be generated algebraically from a small set of elementary modes or rapidities.