---
title: 'OSQMBS: Open-System Quantum Many-Body Scars'
url: https://www.emergentmind.com/topics/open-system-quantum-many-body-scars-osqmbs
type: topic
---

# OSQMBS: Open-System Quantum Many-Body Scars

Open-System Quantum Many-Body Scars (OSQMBS) denotes a family of nonergodic structures in open, effectively open, or open-system-dual quantum many-body dynamics that extend the scar paradigm beyond closed Hermitian Hamiltonians. In recent arXiv literature, the term is used in several non-identical but related senses: an effective open-system reformulation of kinetically constrained scar dynamics that leads to “Grand Canonical-like Thermalization” on an energy–quasiparticle-number plane [2603.06075]; exact Lindbladian embeddings of scar towers into decoherence-free subspaces with persistent oscillations for generic initial states [2304.08515]; a commutant-algebra definition of stationary pure states coexisting with a typical infinite-temperature stationary state on the orthogonal complement, together with asymptotic variants [2509.18023]; and a detailed-balance duality in which special Lindbladian eigenoperators become exact scars of a closed bilayer Hermitian Hamiltonian, notably a thermofield double state [2304.03155]. Experimentally, related work on single scars in superconducting-qubit platforms is relevant because it develops non-revival diagnostics under controlled imperfections, even though it does not formulate a Lindblad scar theory [2410.14613].

## 1. Conceptual scope and relation to closed-system scars

Quantum many-body scars in closed systems are atypical excited eigenstates embedded in otherwise thermal spectra, often associated with ETH violation, low entanglement, and, in tower constructions, coherent revivals from special initial states. In the open-system setting, the same phenomenology is not transferred in a unique way. One branch of the literature treats scars as exact or asymptotically long-lived structures of a Lindbladian itself; another uses an effective open-system representation to reinterpret constrained unitary dynamics; a third identifies special open-system eigenoperators whose doubled-space images are exact scar states of a Hermitian bilayer model [2509.18023].

A recurring misconception is that OSQMBS must mean a conventional scar that merely survives weak decoherence. The current literature is broader and more structured. In one formulation, OSQMBS are exact stationary pure states selected by an unconventional strong-symmetry structure, with all orthogonal initial states converging to an infinite-temperature state on the complement [2509.18023]. In another, a scar tower is embedded into a decoherence-free subspace (DFS), so that the asymptotic dynamics is not stationary but coherently oscillatory inside the DFS [2304.08515]. In yet another, the “open-system” content is an exact non-Markovian representation of closed constrained dynamics, from which a revised ETH and an emergent quasiparticle number arise [2603.06075].

A common motif is the coexistence of an atypical sector with a thermal or infinite-temperature sector. This suggests a unifying interpretation, although the literature does not yet provide a single universally accepted definition. What is shared is that the exceptional sector is isolated relative to a thermalizing bulk, but the mathematical realization differs: stationary density matrices, nondecaying Liouvillian eigenoperators, operator-space eigenmodes, or low-density sectors in a generalized thermodynamic landscape.

## 2. Effective open-system formulation in kinetically constrained scar models

A direct open-system reinterpretation of scarred constrained dynamics is developed for a one-dimensional spin chain of length \(D\) with periodic boundary conditions and on-site spin \(j\ge 1\) [2603.06075]. The unconstrained dynamics is
\[
H_0=\sum_{k=1}^{D}s_k^x .
\]
Kinetic constraints forbid neighboring sites from occupying
\[
|x\rangle_{k,k+1}=|j\rangle_k\otimes |-j\rangle_{k+1},
\]
and define the constrained Hilbert space through
\[
\hat P=\prod_{k=1}^{D}\left(\mathbb I-|x\rangle_{k,k+1}\langle x|\right).
\]
The constrained Hamiltonian is obtained by projection,
\[
\hat P H_0\hat P=H\hat P=\hat P H.
\]

The same dynamics can be represented by a non-Hermitian effective generator together with two compensating dissipative channels. The resulting equation is written as
\[
\partial_t \rho = \mathcal{L}(\rho) \equiv -i H_N \rho + i \rho H_N^{\dagger} + \mathcal{J}(\rho),
\]
with a jump term \(\mathcal J\) that annihilates all states in the constrained subspace because \(\pi_k\hat P=\hat P\pi_k=0\). Restricted to \(\mathcal H\), the construction reproduces the original unitary dynamics exactly. In this sense, the dynamics is “open-system” only effectively: the constrained closed system is represented as if it were coupled to an auxiliary environment. One dissipation rate is negative, \(\gamma_2<0\), so the effective dynamics is explicitly non-Markovian and incorporates information backflow [2603.06075].

This representation motivates an emergent quasiparticle number. For \(\rho(t)\in\mathcal H\), the two jump probabilities are equal, \(P_+=P_-\), and the information-exchange rate is proportional to \(\mathrm{Tr}(\rho\hat N)\), where
\[
\hat N=\sum_{k=1}^{D} M_k^\dagger M_k=\sum_{k=1}^{D}|y\rangle\langle y|_{k,k+1}.
\]
The operator \(\hat N\) counts local \(|y\rangle\)-patterns and measures how strongly the state participates in the virtual transitions implementing the blockade. It is neither exactly conserved nor universal; it is emergent and model-specific. For an eigenstate \(|E_i\rangle\), the quasiparticle number is
\[
N_i=\langle E_i|\hat N|E_i\rangle.
\]
Numerically, scar states have anomalously low \(N_i\), which is interpreted as reduced information exchange with the auxiliary environment [2603.06075].

The construction was studied numerically for \(D=10\), \(j=1\), where the model maps to the PXP model on a chain of doubled length. Two special initial states are emphasized:
\[
|\psi_1(0)\rangle=|j,j,\dots,j\rangle,
\qquad
|\psi_2(0)\rangle=\frac{1}{\sqrt D}\sum_{k=1}^{D}|1,1,\dots,0_k,\dots,1\rangle .
\]
The first exhibits strong revivals, while both have large overlap with scarred eigenstates [2603.06075].

## 3. Revised ETH, generalized density of states, and coherence diagnostics

Within the effective open-system description, the conventional ETH is replaced by a two-parameter, “grand canonical-like” version [2603.06075]. The diagonal ETH ansatz becomes
\[
\langle E_i | \hat{O} | E_i \rangle \simeq \mathcal{O}(\mathcal{E},\mathcal{N}),
\]
with \(\mathcal E=E_i\) and \(\mathcal N=N_i\). For an initial state \(|\psi(0)\rangle=\sum_i c_i|E_i\rangle\), the long-time average is
\[
\bar O \simeq \mathcal O(\mathcal E,\mathcal N)+\mathbb O(\Delta_E^2)+\mathbb O(\Delta_N^2),
\]
where \(\mathcal E=\sum_i |c_i|^2 E_i\) and \(\mathcal N=\sum_i |c_i|^2 N_i\). The associated effective ensemble is
\[
\rho(\beta,\mu)\propto \sum_i e^{-\beta(E_i-\mu N_i)}|E_i\rangle\langle E_i|.
\]
Because \(\hat N\) is not conserved, this is not a fundamental thermodynamic grand canonical ensemble; it is an eigenstate-weighted effective ensemble over \((E_i,N_i)\). Nonetheless, it predicts long-time averages and reduced density matrices more accurately than the canonical ensemble for the scarred states studied numerically [2603.06075].

The off-diagonal ETH is generalized as
\[
\langle E_i | \hat{O} | E_{i'} \rangle_{i\neq i'}
=
\Omega(\mathcal{E},\mathcal{N})^{-1/2}
\, f(\mathcal{E},\mathcal{N},\omega,\nu)\, R_{i,i'} ,
\]
so the suppression factor is no longer a one-dimensional density of states \(\Omega(\mathcal E)\), but a generalized density of states on the \((\mathcal E,\mathcal N)\) plane. The smoothed DOS is defined by Gaussian kernel smoothing over eigenstate coordinates \((E_i,N_i)\). Low-density regions in this plane are predicted to support relatively large off-diagonal matrix elements, and therefore anomalously large temporal fluctuations.

To make the off-diagonal analysis less observable-dependent, the paper introduces the cross coherence purity (CCP),
\[
\mathcal T_{i,i'}\equiv \mathrm{Tr}\!\left(\rho_A^{i,i'\dagger}\rho_A^{i,i'}\right),
\]
where \(\rho_A^{i,i'}=\mathrm{Tr}_{\bar A}(|E_i\rangle\langle E_{i'}|)\) is the reduced cross-density matrix on subsystem \(A\). For \(i=i'\), CCP reduces to the usual subsystem purity. For \(i\neq i'\), it quantifies the coherence retained between the two eigenstates on \(A\). A rigorous bound implies that CCP controls the typical size of all local off-diagonal matrix elements on \(A\), and the authors infer
\[
\mathcal T_{i,i'}\propto \Omega(\mathcal E,\mathcal N)^{-1}.
\]
They stress that this is not an exact theorem, since the CCP inequality is only an upper bound. Numerically, however, the inverse averaged CCP tracks the generalized DOS closely once data are binned in the \((\mathcal E,\mathcal N)\) plane and restricted to fixed \((\omega,\nu)\) windows [2603.06075].

This leads to a revised fluctuation estimate,
\[
\Delta_t^2 O
=
\sum_{i\neq i'} |c_i|^2 |c_{i'}|^2 |O_{i,i'}|^2
\propto
\sum_{i\neq i'} |c_i|^2 |c_{i'}|^2 \, \Omega(\mathcal E,\mathcal N)^{-1},
\]
and an alternative CCP-based estimate,
\[
\Delta_t O \sim \sqrt{\sum_{i\neq i'} |c_i|^2 |c_{i'}|^2 \mathcal T_{i,i'}}.
\]
Across all translation-invariant initial states tested, both estimates scale proportionally with the observed long-time fluctuation amplitudes. The anomalous fluctuations and quasi-periodic dynamics of scarred states are therefore recast as consequences of low-DOS regions in the generalized \((E,N)\) landscape rather than as isolated failures of thermalization [2603.06075].

The same low-DOS picture is connected to the approximate spectrum-generating algebra (SGA),
\[
[H,Q^\pm]=\pm Q^\pm+\hat R,
\]
with small ladder deviation for scar states quantified by parameters \(r_i^\pm\). Numerically, scar states have much smaller \(r_i^\pm\) than thermal states, and the enhancement is attributed to a low generalized DOS that amplifies the relevant denominator \(d_i^\pm\). In this interpretation, the SGA is not a competing explanation but a structure amplified in low-DOS sectors [2603.06075].

## 4. Lindbladian embedding into decoherence-free scar sectors

A distinct OSQMBS program constructs Lindblad dynamics whose asymptotic nondecaying sector is exactly the scar sector of a many-body Hamiltonian [2304.08515]. The starting point is the Lindblad equation
\[
\frac{d\rho}{dt} = -i[H,\rho] + \gamma\sum_j \left( 2L_j\rho L_j^\dagger-\{L_j^\dagger L_j,\rho\} \right).
\]
If a scar tower \(\{|S_n\rangle\}\) satisfies
\[
L_j|S_n\rangle=0,
\qquad
H|S_n\rangle=E_n|S_n\rangle,
\]
then the Liouvillian contains exact nondecaying eigenoperators
\[
\mathcal L(|S_n\rangle\langle S_m|)=-i(E_n-E_m)|S_n\rangle\langle S_m|.
\]
When the tower is equally spaced, the Liouvillian inherits a uniformly spaced set of purely imaginary eigenvalues, so observables display persistent periodic oscillations for generic initial states, including mixed states, after transient relaxation into the DFS [2304.08515].

The engineered jumps are taken as
\[
L_j=V_jP_j,
\]
where \(P_j\) are local projectors annihilating the whole scar tower and \(V_j\) are local operators chosen so that the dark space is not accidentally enlarged. The Hamiltonian is essential: it preserves the desired scar subspace while rotating unwanted states out of the null space of the dissipators, after which those unwanted states decay. This mechanism differs from ordinary dark-state engineering, where the focus is usually on a stationary dark state rather than a coherently oscillating dark tower.

The tower structure is encoded by a restricted spectrum-generating algebra on the scar subspace \(W\),
\[
([H,Q^\dagger]-\omega Q^\dagger)W=0,
\qquad
|S_n\rangle=(Q^\dagger)^n|S_0\rangle.
\]
To find projectors annihilating the whole tower, the paper compresses the tower into
\[
|S(\beta)\rangle = e^{\beta Q^\dagger}|S_0\rangle = \sum_{n=0}\frac{\beta^n}{n!}|S_n\rangle,
\]
and then uses an MPS-based construction to determine local annihilating projectors. This “compressed-MPS” procedure is the technical backbone of the general framework [2304.08515].

The paper demonstrates the construction in four representative models: a Dicke-state toy model, a spin-1 \(XY\) model, the AKLT model, and a domain-wall preserving spin-\(\tfrac12\) model. These examples illustrate three different situations: exact coincidence between the dissipator null space and the scar subspace; a larger null space that is filtered by the Hamiltonian; and the need for additional higher-locality projectors to remove unwanted exact dark states. The AKLT example is especially important because it shows that two-local projectors can leave irrational unwanted modes, requiring three-local refinement [2304.08515].

A practically relevant extension is dissipative preparation. If a strong symmetry \(A\) commutes with both \(H\) and all \(L_j\), and distinct scar states carry distinct eigenvalues of \(A\), then each symmetry sector evolves independently, and a simple product state in the appropriate sector relaxes to the target scar state. The paper further proposes a digital quantum simulation protocol with ancilla qubits that are repeatedly reset. A Hamiltonian step, a system–ancilla coupling
\[
H_{\rm coup}^j=\sqrt{2\gamma}(L_j\tau_j^+ + L_j^\dagger\tau_j^-),
\]
and an ancilla reset reproduce the Lindblad evolution to order \(O(\delta t^2)\). Superconducting qubits, Rydberg arrays, and trapped ions are identified as candidate platforms [2304.08515].

## 5. Formal Lindbladian theory: exact stationary OSQMBS and asymptotic OSQMBS

A later work undertakes the problem of formally introducing quantum many-body scarring for Markovian open systems governed by the Lindblad equation and proposes a precise definition of OSQMBS in terms of stationary pure states and thermal complements [2509.18023]. In that formulation, a pure state \(\hat\rho_\psi=\ket{\psi}\bra{\psi}\) is an OSQMBS iff
\[
\mathcal{L}(\hat\rho_\psi)=0
\]
and, for every initial state \(\hat\rho(0)\) such that
\[
\mathrm{Tr}[\hat\rho_\psi \hat\rho(0)] = 0,
\]
one has
\[
\lim_{t\to+\infty} e^{\mathcal{L}t}\hat\rho(0)=\frac{\mathds{\hat 1}-\hat\rho_\psi}{d-1},
\]
with \(d=\dim(\mathcal H)\). The scar is therefore an exact stationary pure state, while every orthogonal initial state thermalizes to the infinite-temperature state on its orthogonal complement [2509.18023].

The theory is built from a bond algebra \(\mathcal A\) generated by local Hamiltonian terms and Hermitian jump operators, and its commutant \(\mathcal C\). The Hilbert space decomposes as
\[
\mathcal{H}=\bigoplus_\lambda \left(\mathcal H_\lambda^{(\mathcal A)}\otimes \mathcal H_\lambda^{(\mathcal C)}\right).
\]
OSQMBS arise when the commutant has an unconventional strong-symmetry structure: a one-dimensional scar block plus a single irreducible thermal block on the orthogonal complement. In the single-scar case, the stationary manifold is
\[
\hat\rho_{ss}(p)= p\,\ket{\psi}\bra{\psi} +(1-p)\frac{\mathds{\hat 1}-\ket{\psi}\bra{\psi}}{d-1},
\qquad 0\le p\le 1,
\]
with \(p=\mathrm{Tr}[\ket{\psi}\bra{\psi}\hat\rho(0)]\). This shifts the scar notion away from revival physics and toward the stationary-state manifold of the Liouvillian [2509.18023].

The framework also treats towers of scars. Depending on the commutant structure, the scar manifold can consist of multiple one-dimensional sectors or a higher-dimensional scar block. The paper studies benchmark spin-\(\tfrac12\) and spin-1 constructions in which the atypical stationary states are ferromagnets or closed-system scar states adapted to the open-system commutant structure.

A further technical contribution concerns coherences between a scar sector and the thermal sector. For \(\hat\Pi_s=\ket{\psi}\bra{\psi}\) and \(\hat\Pi_{th}=1-\hat\Pi_s\), the off-diagonal components evolve under an effective non-Hermitian operator
\[
\hat H_{\mathrm{eff}}
=
-i\sum_\alpha J_\alpha(\hat h_\alpha-\epsilon_\alpha)
-\frac12\sum_j \gamma_j(\hat l_j-\lambda_j)^2
=
-i\hat H_1-\frac12\hat H_2,
\]
and decay exponentially when \(\hat H_2\) is gapped on the thermal block. For dephasing examples with isolated ferromagnetic scars, the paper derives explicit lower bounds on the decay rate of scar–thermal coherences [2509.18023].

The same work introduces asymptotic OSQMBS (AOSQMBS): states that are not exact stationary states, do relax to the typical infinite-temperature state, but display anomalously large relaxation times. The proposed diagnostic is fidelity scaling. For one class of Lindbladians,
\[
\mathcal F(t)=f(t^2/L^2),
\]
while for another,
\[
\mathcal F(t)=g(t/L^2).
\]
These Ansätze are motivated by short-time expansions and are supported numerically for spin-1 tower models at modest system sizes. The paper emphasizes that the strongest exact results concern stationary-state structure and coherence sectors, whereas AOSQMBS are characterized operationally through scaling and finite-size benchmarks rather than by a single formal definition of the same precision as the exact OSQMBS definition [2509.18023].

## 6. Detailed-balance duality and thermofield double scars

A complementary route to OSQMBS is based on an exact duality between open many-body systems satisfying detailed balance and closed bilayer systems with a self-adjoint Hamiltonian [2304.03155]. The open-system dynamics is written in the Heisenberg picture as
\[
\partial_t O_t = \mathcal L^\dagger O_t \equiv i[H,O_t] + \mathcal D^\dagger O_t,
\]
with Lindblad operators \(L_j\) obeying
\[
[\beta H_0,L_j] = -\beta L_j,
\qquad
[H_0,H]=0.
\]
After a temperature-dependent similarity transformation and a Wick rotation \(g=i\gamma\), the vectorized generator becomes a Hermitian bilayer Hamiltonian \(\mathcal H_\beta\) acting on \(\mathscr H\otimes \mathscr H\) [2304.03155].

The operator-to-vector map is
\[
|O\rangle\rangle_\beta
=
\sum_{m,n} O_{mn}\,
e^{-\frac{\beta}{4}H_0}|m\rangle
\otimes
e^{-\frac{\beta}{4}H_0^T}|n\rangle.
\]
Under this map, the identity operator is sent to the thermofield double (TFD) state,
\[
|\mathbb 1\rangle\rangle_\beta
=
\sum_n e^{-\beta E_n/2}|n\rangle\otimes |n\rangle,
\]
and because \(\mathcal L^\dagger \mathbb 1=0\), one has
\[
\mathcal H_\beta |\mathbb 1\rangle\rangle_\beta = 0.
\]
The TFD is thus an exact bulk eigenstate of the dual Hermitian bilayer Hamiltonian and is identified as a quantum many-body scar [2304.03155].

The entanglement structure is a central feature. Tracing out one layer yields a Gibbs state,
\[
\rho_{\rm red}=\frac{e^{-\beta H_0}}{Z(\beta)},
\]
so the interlayer entanglement entropy is
\[
S_{\rm ent} = \beta \langle H_0\rangle_\beta + \ln Z(\beta).
\]
At \(\beta\to 0\), the TFD reduces to the maximally entangled EPR scar; at \(\beta\to\infty\), the interlayer entanglement is minimized or vanishes for a unique ground state. This gives an exact scar with temperature-tunable entanglement, a feature that is unusual in the conventional scar literature [2304.03155].

More generally, explicit Lindbladian eigenoperators \(Q\) satisfying
\[
\mathcal L^\dagger Q = -\Gamma Q
\]
have expectation values
\[
\langle Q\rangle_t=e^{-\Gamma t}\langle Q \rangle_0
\]
for arbitrary initial states. Under the duality, these eigenoperators generate additional exact scar states or scar towers of the bilayer Hamiltonian. The paper provides broad model classes, including dissipative \(XY\)-type spin models, dephasing-like \(\sigma^z\) dissipators, bond dissipators, and generalized Hubbard models [2304.03155].

This duality-based construction is foundational for OSQMBS, but it comes with a clear caveat: the scar lives in the doubled Hilbert space of the dual bilayer system, not in the original single-layer open system. The open-system analogue is therefore most naturally described as an “operator scar,” while the scar nomenclature is literal on the Hermitian bilayer side.

## 7. Experimental diagnostics, misconceptions, and current limitations

Experimental work on single scars in fixed-frequency, fixed-coupling superconducting transmons is closely related to OSQMBS because it develops diagnostics adapted to imperfect, noisy, and stroboscopic settings, even though it is not a formal open-system theory [2410.14613]. The proposal implements single-scar parent Hamiltonians using trotterized cross-resonance interactions and studies a fully \(x\)-polarized product scar and a cluster-state scar. Because a single scar does not produce tower-like revival phenomenology, the paper proposes alternative signatures: persistence of nonthermal local observables, persistence of low entanglement, comparison with locally deformed nearby states, sensitivity to controlled coefficient randomization,
\[
c_{i,j}^\alpha \to (1+r\,u_{i,j}^\alpha)c_{i,j}^\alpha,
\]
and sensitivity to Trotter resolution. Exact Floquet analysis for \(L=12\) finds one mode with overlap \(0.9993\) with the \(x\)-polarized scar and \(0.9992\) with the cluster scar; under controlled protocol distortion the scar remains noticeably more stable than its deformed counterpart up to about \(r\lesssim 0.05\). The paper is explicit that this is not Lindblad decoherence, but a controlled coherent distortion of the protocol [2410.14613].

The present OSQMBS literature contains several nontrivial limitations. In the effective constrained-dynamics approach, the quasiparticle number \(\hat N\) is emergent, model-dependent, and not conserved, while the open-system representation is non-Markovian because one jump rate is negative; most numerics are for \(D=10\), \(j=1\) [2603.06075]. In the commutant-based formal theory, the strongest results assume Hermitian jump operators and focus on exact stationary-state structure rather than a full Liouvillian spectral theory; the AOSQMBS evidence is finite-size [2509.18023]. In the DFS-embedding program, local projectors annihilating the whole tower must be found model by model, unwanted dark states can persist, and higher-locality corrections may be needed, as in the AKLT construction [2304.08515]. In the duality framework, Hermiticity of the bilayer Hamiltonian requires detailed balance, and the scar is directly realized only in the doubled space [2304.03155].

These caveats also delimit the main conceptual controversies. One is terminological: whether OSQMBS should refer only to intrinsic Lindbladian stationary or nondecaying modes, or also to effective open-system reinterpretations of closed constrained dynamics. Another concerns universality: current constructions are exact and illuminating, but strongly model-specific. A further open issue is extension to genuinely dissipative scar systems in which environment-induced decoherence competes with scar coherence rather than merely protecting, dualizing, or reparametrizing it. The literature repeatedly points to this unresolved direction, along with scaling to larger systems, relaxation beyond Hermitian jumps, and the status of quasiparticle-number constructions beyond constrained models [2603.06075].

Source: https://www.emergentmind.com/topics/open-system-quantum-many-body-scars-osqmbs