---
title: Quantum Subspace Verification
url: https://www.emergentmind.com/topics/quantum-subspace-verification
type: topic
---

# Quantum Subspace Verification

Quantum subspace verification is the task of deciding whether a quantum state lies in a target subspace \(V\subseteq H\), rather than verifying equality with a single target pure state. In the current literature, the target subspace may be a stabilizer code space, a genuinely entangled subspace, a symmetry sector, a computational subspace induced by an encoding, or a constraint-satisfying subspace. The subject has developed into a distinct framework that combines hypothesis testing, operator design, and restricted-measurement analysis, with particular emphasis on local Pauli measurements, adaptive LOCC protocols, and experimentally realistic certification of noisy mixed states [2409.19699][2410.12551][2509.01281].

## 1. Definition and statistical formulation

A standard formulation fixes a target subspace \(\mathcal V\subseteq \mathcal H\) with projector
\[
\Pi := \sum_j |\psi_j\rangle\langle \psi_j|,
\]
where \(\{|\psi_j\rangle\}_j\) is an orthonormal basis of \(\mathcal V\). The verification problem is then posed as a binary distinction between a good case, in which every tested sample satisfies \(\operatorname{Tr}[\Pi \sigma_i]=1\), and a bad case, in which every sample satisfies \(\operatorname{Tr}[\Pi \sigma_i]\le 1-\epsilon\) for some tolerance \(\epsilon>0\). Pure-state verification is recovered as the special case \(\dim V=1\) [2409.19699][2509.01281].

A verification strategy is specified by a family of POVM test operators \(M\in\mathcal M\), sampled with probabilities \(\mu(M)\), and summarized by the verification operator
\[
\Omega=\sum_{M\in\mathcal M}\mu(M)M.
\]
The acceptance probability on \(\sigma\) is \(\operatorname{Tr}[\Omega\sigma]\). For a valid subspace verifier, perfect completeness requires
\[
\operatorname{Tr}[\Omega \sigma]=1,\qquad \forall \sigma\in \mathscr D(\mathcal V),
\]
equivalently \(\operatorname{Tr}[\Omega\Pi]=\operatorname{rank}(\Pi)\). Soundness is governed by the restriction of \(\Omega\) to the orthogonal complement,
\[
\tilde\Omega:=(1-\Pi)\Omega(1-\Pi),
\]
with worst-case bad-case acceptance
\[
p(\Omega)=1-\bigl(1-\lambda_{\max}(\tilde\Omega)\bigr)\epsilon.
\]
Defining the spectral gap
\[
\nu(\Omega):=1-\lambda_{\max}(\tilde\Omega),
\]
one obtains
\[
\Pr\{\text{accept}\mid \sigma_1,\dots,\sigma_N\}\le (1-\nu(\Omega)\epsilon)^N,
\qquad
N(\Omega)\ge \frac{1}{\nu(\Omega)}\cdot \frac{1}{\epsilon}\ln\frac{1}{\delta}.
\]
Within this framework, maximizing \(\nu(\Omega)\) minimizes the number of copies needed for confidence \(1-\delta\) [2409.19699].

A complementary formulation emphasizes subspace fidelity. For a projector \(\mathcal P\) onto \(\mathcal V\), the subspace fidelity is
\[
F(\mathcal P,\rho)=\operatorname{Tr}(\mathcal P\rho),
\qquad
\epsilon_\rho=1-F(\mathcal P,\rho).
\]
In the hypothesis-testing version for quantum error-correction subspaces, the good device satisfies \(F(\mathcal P,\rho)\ge 1-\tau\epsilon\) and the bad device satisfies \(F(\mathcal P,\rho)\le 1-\epsilon\), where \(0<\tau<1\). If
\[
\tau\Delta_{\max}(\Omega)<\Delta_{\min}(\Omega),
\]
with \(\Delta_{\min}\) and \(\Delta_{\max}\) the two relevant spectral gaps outside the target subspace, then the required number of rounds is expressed through a Bernoulli KL divergence and scales as \(O(1/\epsilon)\) for small \(\epsilon\). If the verification operator has maximum eigenvalue \(\lambda<1\), the sample complexity reverts to \(O(1/\epsilon^2)\), reflecting the loss of perfect acceptance for ideal code states [2410.12551].

## 2. Stabilizer-code subspaces and local-measurement protocols

The most developed concrete setting is the verification of stabilizer code spaces. For an \([\![n,k,d]\!]\) stabilizer code, the code space is the \(+1\) eigenspace of a stabilizer group \(\mathcal S_k\subseteq\mathcal P_n\), generated by \(n-k\) independent stabilizer generators \(\mathcal G_k=\{S_1,\dots,S_{n-k}\}\). Verifying the stabilizer subspace is equivalent to certifying the logical qubits encoded by the code. The central practical constraint is that the protocols are required to use only local Pauli measurements, be non-adaptive, and work on mixed states [2409.19699].

Two universal stabilizer strategies are immediate. Strategy I samples a non-identity stabilizer \(P\in\mathcal S_k\setminus\{1\}\) uniformly, measures it, and accepts on outcome \(+1\). Its verification operator is
\[
\Omega_{\rm I}:=\frac{1}{2^{n-k}-1}\sum_{P\in \mathcal S_k\setminus\{1\}} P^+,
\qquad
P^+:=\frac{1+P}{2},
\]
with spectral gap
\[
\nu(\Omega_{\rm I})=\frac{2^{n-k-1}}{2^{n-k}-1},
\]
so that
\[
N(\Omega_{\rm I})\approx 2\,\frac{1}{\epsilon}\ln\frac{1}{\delta}.
\]
This is near-global optimal in copy complexity but requires \(2^{n-k}-1\) measurement settings. Strategy II samples a generator \(S\in\mathcal G_k\) uniformly and uses
\[
\Omega_{\rm II}:=\frac{1}{n-k}\sum_{S\in\mathcal G_k}S^+,
\qquad
\nu(\Omega_{\rm II})=\frac{1}{n-k},
\]
hence
\[
N(\Omega_{\rm II})=(n-k)\frac{1}{\epsilon}\ln\frac{1}{\delta}.
\]
It uses only \(n-k\) settings but needs more copies [2409.19699].

Additional code structure yields sharper protocols. For graph code subspaces, a derived graph \(G'\) on the generators is covered by independent sets \(\mathscr A=\{A_1,\dots,A_m\}\), giving
\[
P_\ell=\prod_{a\in A_\ell}\frac{1+S'_a}{2},
\qquad
\Omega_{\mathscr G}=\frac{1}{m}\sum_{\ell=1}^m P_\ell,
\qquad
\nu(\Omega_{\mathscr G})=\frac{1}{m}.
\]
Since \(m\le n-k\), the coloring strategy uses no more and often fewer settings than generator-based verification. For CSS codes with \(H_XH_Z^T=0\), the \(X\)- and \(Z\)-checks factorize, and the XZ strategy chooses \(P\in\{X,Z\}\) uniformly, measures \(P^{\otimes n}\), and accepts iff all checks of that type return \(+1\). Its verification operator is
\[
\Omega_{\rm CSS}^{\rm XZ}
=\frac{1}{2}\left(
\prod_{\mathbf c_X}\frac{1+X^{\mathbf c_X}}{2}
+
\prod_{\mathbf c_Z}\frac{1+Z^{\mathbf c_Z}}{2}
\right),
\]
with
\[
\nu(\Omega_{\rm CSS}^{\rm XZ})=\frac{1}{2},
\qquad
N(\Omega_{\rm CSS}^{\rm XZ})=2\,\frac{1}{\epsilon}\ln\frac{1}{\delta}.
\]
For dual-containing CSS codes with \(H_X=H_Z\), the XYZ strategy samples uniformly from \(\{X,Y,Z\}\), achieves
\[
\nu(\Omega_{\rm DC}^{\rm XYZ})=\frac{2}{3},
\qquad
N(\Omega_{\rm DC}^{\rm XYZ})=\frac{3}{2}\,\frac{1}{\epsilon}\ln\frac{1}{\delta},
\]
and reduces the copy cost by \(25\%\) relative to the XZ strategy [2409.19699].

A closely related framework for error-correction benchmarking defines verification operators for stabilizer codes and QLDPC codes in terms of subspace fidelity rather than perfect-membership testing. For stabilizer codes, measuring all stabilizers gives \(\Omega_{\text{all}}=(I+\mathcal P)/2\) with \(\Delta_{\min}=\Delta_{\max}=1/2\), whereas measuring only generators gives \(\Delta_{\min}=1/m\), \(\Delta_{\max}=1\), where \(m=n-k\). A chromatic construction based on the bit-wise commutativity graph \(G_{\mathcal S}\) yields
\[
\Omega_{\text{chr}(\mathcal S)}
=
\frac{1}{\chi(G_{\mathcal S})}
\sum_{u=1}^{\chi(G_{\mathcal S})}
\prod_{v_i\in\mathcal I_u}\frac{I+S_i}{2},
\]
with \(\Delta_{\min}=1/\chi(G_{\mathcal S})\). For CSS codes, including surface codes, toric codes, and hypergraph product codes, \(\chi^*=2\), so both the setting number and the sample complexity become independent of \(m\) up to constants. For generic QLDPC codes with local projector support size \(s\), the support graph has degree at most \(s^2\), which yields at most \(s^2\) measurement settings in the chromatic construction [2410.12551].

| Subspace family | Measurement settings | Copy/sample complexity |
|---|---:|---:|
| Global optimum | 1 | \(\frac{1}{\epsilon}\ln\frac{1}{\delta}\) |
| Stabilizer code, full stabilizer group | \(2^{n-k}-1\) | \(\approx 2\,\frac{1}{\epsilon}\ln\frac{1}{\delta}\) |
| Stabilizer code, generators | \(n-k\) | \((n-k)\frac{1}{\epsilon}\ln\frac{1}{\delta}\) |
| Graph code coloring | \(m\le n-k\) | \(m\,\frac{1}{\epsilon}\ln\frac{1}{\delta}\) |
| CSS code, XZ strategy | 2 | \(2\,\frac{1}{\epsilon}\ln\frac{1}{\delta}\) |
| Dual-containing CSS code, XYZ strategy | 3 | \(\frac{3}{2}\,\frac{1}{\epsilon}\ln\frac{1}{\delta}\) |

These results establish a recurring pattern: additional algebraic structure decreases both the number of measurement settings and the required number of copies, and several local-measurement protocols approach the global optimum under realistic experimental constraints [2409.19699][2410.12551].

## 3. Entangled subspaces and local-verifiability limits

Quantum subspace verification is not restricted to code spaces. A prominent example is the three-qubit GHZ-W genuinely entangled subspace
\[
\mathcal V_3=\mathrm{span}\{|\mathrm{GHZ}\rangle,|W\rangle\},
\]
with
\[
|\mathrm{GHZ}\rangle=\frac{|000\rangle+|111\rangle}{\sqrt2},
\qquad
|W\rangle=\frac{|001\rangle+|010\rangle+|100\rangle}{\sqrt3}.
\]
For this target, efficient verification requires adaptive local measurements and one-way classical communication. The adaptive framework measures one qubit first, then chooses a conditional two-qubit test on the remaining pair. The general form is
\[
M_P=P^+M_P^+ + P^-M_P^-,
\qquad P\in\{X,Z\}.
\]
Two concrete protocols are given: an XZ strategy and a rotation strategy [2412.19540].

The XZ strategy uses one \(Z\)-based adaptive test and three \(X\)-based adaptive tests, one for each choice of first measured qubit. Its verification operator is
\[
\Omega_{\rm XZ}
=
\mu(Z)M_Z+\frac{1}{3}\mu(X)\sum_{i=1}^3 M_{X,i}.
\]
Numerical optimization yields
\[
\nu(\Omega_{\rm XZ})\approx 0.262
\quad\text{at}\quad
\mu(Z)\approx 0.424,
\]
hence
\[
N(\Omega_{\rm XZ})\approx 3.817\times \frac{1}{\epsilon}\ln\frac{1}{\delta}.
\]
The rotation strategy exploits permutation symmetry and local phase rotations \(R_{2\pi/3}\) and \(R_{4\pi/3}\). Averaging over three qubit choices and six rotated variants produces a ten-test family with exact spectral gap
\[
\nu(\Omega_\mu)=\min\left\{\frac{47}{80}\mu(X),\,1-\frac{11}{15}\mu(X)\right\},
\]
maximized at
\[
\mu^\star(X)=\frac{240}{317},
\qquad
\nu(\Omega_{\mu^\star})=\frac{141}{317}\approx 0.445.
\]
The resulting sample complexity is
\[
N(\Omega_{\mu^\star})=\frac{317}{141}\times \frac{1}{\epsilon}\ln\frac{1}{\delta}
\approx 2.248\times \frac{1}{\epsilon}\ln\frac{1}{\delta}.
\]
This is significantly closer to the global optimum than the XZ strategy while remaining local and adaptive [2412.19540].

The same work also identifies an intrinsic limitation of local subspace verification through a complete classification of two-dimensional two-qubit subspaces. Some are unverifiable by LOCC because their orthogonal complement contains only one product state. For such a subspace, the best possible local test has zero spectral gap,
\[
\Omega_u=1-|\tau\rangle\langle\tau|,
\qquad
\nu(\Omega_u)=0.
\]
Other subspaces are verifiable with two tests, and perfectly verifiable subspaces admit a single local test operator. A common misconception is therefore that subspace verification simply generalizes state verification without changing locality properties. The classification shows that local certifiability can fail at the level of subspaces even when the constituent states are individually meaningful entanglement resources [2412.19540].

## 4. Restricted measurements, visibility, and duality with data hiding

A more abstract line of work analyzes subspace verification through restricted-measurement distinguishability. For a measurement class \(\mathbb M\) and nontrivial proper subspace \(V\subsetneq H\), the central quantity is the \(\epsilon\)-visibility
\[
\gamma_{V,\mathbb M(\epsilon)}
=
\frac{1}{\epsilon}
\sup_{\Omega\in\mathbf M}
\min_{\substack{\rho\in\mathsf D(H):\,\rho\subseteq V\\
\sigma\in\mathsf D(H):\,\sigma\Pi_V\le 1-\epsilon}}
\Omega(\rho-\sigma),
\]
which also admits the minimax form
\[
\gamma_{V,\mathbb M(\epsilon)}
=
\min_{\substack{\rho\subseteq V\\
\sigma\Pi_V\le 1-\epsilon}}
\frac{1}{2\epsilon}\|\rho-\sigma\|_{\mathbb M}.
\]
For universal strategies satisfying \(\Omega\Pi_V\propto\Pi_V\), the limiting quantity \(\hat\gamma_{V,\mathbb M}\) is obtained as \(\epsilon\to 0\) when \(\mathbb M\) is informationally complete. The operational content is that sample complexity for sequential measurements lies between a lower bound scaling like \(1/\gamma_{V,\mathbb M(\epsilon)}\) and an upper bound scaling like \(1/\gamma_{V,\mathbb M(\epsilon)}^2\), up to logarithms and constants [2509.01281].

The same framework introduces the \(\epsilon\)-distinguishability ratio \(\mu_{\rho,\mathbb M(\epsilon)}\) for quantum data hiding and proves a duality with subspace verification. For a subspace \(V\),
\[
\mu_{V,\mathbb M(\epsilon)}:=\min_{\rho\subseteq V}\mu_{\rho,\mathbb M(\epsilon)}
\le
\gamma_{V,\mathbb M(\epsilon)}.
\]
More strongly, for informationally complete \(\mathbb M\) and rank bound \(r\),
\[
\min_{\dim V\le r}\hat\gamma_{V,\mathbb M}
=
\mu_{r,\mathbb M},
\qquad
\mu_{r,\mathbb M}:=
\min_{\operatorname{rank}(\rho)\le r}\hat\mu_{\rho,\mathbb M}.
\]
Thus the hardest subspaces to verify and the most secure states for data hiding are the same extremal objects. This establishes a geometric correspondence between verification hardness and hiding security under restricted measurements [2509.01281].

The theory yields explicit protocols and bounds. Under non-adaptive single-qudit local measurements on \((\mathbb C^d)^{\otimes n}\),
\[
\hat\gamma_{V,\mathbb{LO}_n}\ge \frac{1}{2\sqrt{18^n\,\dim V}},
\]
implying
\[
O\!\left(18^n\,\dim(V)\,\frac{\log(1/\delta)}{\epsilon^2}\right)
\]
samples for arbitrary \(V\). For the symmetric subspace \(\vee_n\mathbb C^d\),
\[
\hat\gamma_{\vee_n\mathbb C^d,\mathbb{LO}_n}\ge
\frac{d}{\dim(\vee_n\mathbb C^d)},
\qquad
\dim(\vee_n\mathbb C^d)=\binom{n+d-1}{d-1}.
\]
The associated protocol applies Haar-random local unitaries, measures in the computational basis, and accepts iff all outcomes are equal, implementing
\[
\Omega
=
\sum_{i=0}^{d-1}\int (u|i\rangle\langle i|u^\dagger)^{\otimes n}\,du
=
\frac{d}{\dim(\vee_n\mathbb C^d)}\Pi_{\vee_n\mathbb C^d}.
\]
For qubits this gives \(\hat\gamma\ge 2/(n+1)\). Under PPT measurements,
\[
\hat\gamma_{V,\mathbb{PPT}_n}\ge \frac{1}{2\sqrt{\dim V}},
\]
so any multipartite pure state can be verified with constant sample complexity independent of \(n\) and \(d\), and the universal PPT strategy
\[
\Omega_\Phi=\Phi+\frac13(I-\Phi)
\]
achieves \(\hat\gamma_{\Phi,\mathbb{PPT}_n}\ge 2/3\), with \(2/3\) proved optimal. For 4-design POVMs,
\[
\hat\gamma_{V,\mathbb{DES}_4}\ge \frac{1}{6\sqrt{\dim V}},
\]
again implying \(O(\dim(V)\log(1/\delta)/\epsilon^2)\)-type sample complexity. Stabilizer POVMs give analogous \(O(\dim(V)\log(1/\delta)/\epsilon^2)\) protocols [2509.01281].

## 5. Verification as correction, mitigation, and reduced-space control

In several applications, subspace verification is not only a certification primitive but also the front end of correction or mitigation. Quantum subspace correction constructs stabilizer-like operators for a constraint-satisfying subspace and then uses syndrome extraction to detect and repair violations. For Independent Set on a graph \(\mathcal G=(V,E)\), each edge \(e=(i,j)\) has local legality subspace
\[
\mathcal V_e=\text{span}\{|00\rangle,|01\rangle,|10\rangle\},
\qquad
\bar{\mathcal V}_e=\text{span}\{|11\rangle\},
\]
with stabilizer
\[
\hat S_{(i,j)}=\frac12(I+Z_i+Z_j-Z_iZ_j).
\]
More generally, for a Boolean constraint \(c\), the stabilizer-like operator is
\[
\hat S_c=
\sum_{b\in\mathcal V_c}P_0^{b_0}\cdots P_k^{b_k}
-
\sum_{b\in\bar{\mathcal V}_c}P_0^{b_0}\cdots P_k^{b_k},
\]
and ancilla-assisted syndrome extraction is implemented by
\[
\hat G_c=
\left(\sum_{b\in\mathcal V_c}P_0^{b_0}\cdots P_k^{b_k}\right)I_{a_c}
+
\left(\sum_{b\in\bar{\mathcal V}_c}P_0^{b_0}\cdots P_k^{b_k}\right)X_{a_c}.
\]
The operational loop is explicit: measure syndromes, infer violations, apply recovery, and repeat until the target subspace is reached. In the Independent Set example, the violating region is reset locally, and the protocol becomes a quantum analogue of partial rejection sampling, yielding exact uniform or weighted independent-set distributions when it terminates [2310.20191].

In near-term fermionic simulation, the relevant target is a stabilizer-symmetric computational subspace defined by the fermion-to-qubit encoding together with optional conserved spin-parity symmetries. The projector is written as
\[
\mathsf{M}_{\mathbb S}
=
\prod_{\mathsf S_i\in\mathbb S}\frac{\mathsf I+\mathsf S_i}{2}.
\]
Symmetry verification and post-selection discard shots whose propagated noise anticommutes with at least one stabilizer, while Subspace Noise Tailoring applies probabilistic error cancellation only to the undetectable errors. If \(\mathbb U_k\) denotes the undetectable errors at layer \(k\), the restricted inverse noise map is
\[
\mathcal N_k^{-1}[\bullet]
=
\left(1+\sum_{i\,:\,\mathsf P_i\in\mathbb U_k}p_i^{(k)}\right)\bullet
-
\sum_{i\,:\,\mathsf P_i\in\mathbb U_k}p_i^{(k)}\mathsf P_i\bullet\mathsf P_i.
\]
The residual bias scales as
\[
\mathrm{Bias}[\mathsf O_{\mathrm{est.}}]\sim \mathcal O(\lambda'^2),
\]
while the method preserves the low cost of subspace filtering. The efficacy depends strongly on the encoding: local encodings such as LE, PA, and HX provide many local stabilizers and large detectable-noise fractions, unlike JW, whose global stabilizer is less suitable for scalable parity-check verification [2503.11785].

A further operational use appears in contextual-subspace VQE for the Kagome Heisenberg model. There the reduced Hilbert space is defined by a commuting set \(\mathcal S\) of exact or approximate stabilizers, with projector
\[
\mathbb P=
\frac{1}{2^K}\prod_{S\in\mathcal S}(\mathbbm 1+\nu_S S),
\]
followed by Clifford rotation, projection, and partial trace:
\[
H_{\mathrm{CS}}
=
\operatorname{Tr}_{\mathcal I}\!\left(\mathbb P' RHR^\dagger \mathbb P'\right).
\]
Approximate stabilizers are ranked by a weighted commutation score \(w_n\), and DMRG wavefunctions are used to bias the stabilizer choice toward the desired low-energy sector. The reduced 5-qubit Hamiltonian deliberately retains a \(\mathbb Z_2\) symmetry \(\sigma_z^{(0)}\), enabling hardware symmetry verification by post-selecting on the \(-1\) eigensector of \(\sigma_z^{(0)}\). The reported mitigation stack REM+SV+ZNE reduces the final error ratio from about \(1.210\%\) to \(0.019\%\), which suggests that retaining a verifiable symmetry sector can materially improve reduced-space computations [2506.12391].

## 6. Formal extensions and neighboring verification paradigms

Subspace verification also appears as a logical primitive in quantum program verification. In measurement-based verification of quantum Markov chains, earlier subspace-based temporal logic is treated as a special case of a broader measurement-based linear-time temporal logic. A subspace proposition corresponds to a projection \(P_X\) onto \(X\subseteq\mathcal H\), with
\[
\rho=|\psi\rangle\langle\psi|\models (P_X,[1,1])
\iff
(P_X\rho)=1
\iff
|\psi\rangle\in X.
\]
The quantum Markov chain model \(\mathcal M=(\mathcal H,\mathcal E,\rho_0)\) is then analyzed through spectral properties of the super-operator matrix
\[
M_{\mathcal E}=\sum_k E_k\otimes E_k^*,
\]
and approximate verification is reduced to \(\omega\)-regular model checking via periodic stability and symbolic neighborhoods of the asymptotic cycle. This line of work does not implement subspace verification as a laboratory certification protocol; rather, it embeds subspace membership into a temporal-logic semantics for quantum dynamics [2405.05825].

A distinct but related generalization lifts symmetry verification from states to channels. Symmetric channel verification assumes a known symmetry operator \(S\) commuting with the ideal channel and inserts coherent symmetry-sector tests around the noisy implementation. With symmetry projectors \(\{\Pi_i\}\), the phase-tagging operator is
\[
V_S=\sum_j \exp\!\left[\frac{2\pi i}{2^m}j\right]\Pi_j,
\]
and the detection supermap becomes
\[
(\Theta^{\mathrm{det}_S}(\mathcal U_{\mathcal N}))(\cdot)
=
\sum_{ij}\Pi_i\mathcal U_{\mathcal N}(\Pi_i\cdot\Pi_j)\Pi_j.
\]
This can be interpreted as a channel-level analogue of subspace or symmetry verification, because it enforces the symmetry-preserving block structure of the ideal channel rather than the symmetry of a specific input state. At the same time, it is not a generic replacement for quantum subspace verification: it requires known channel symmetries, targets only symmetry-breaking noise, and in the virtual implementation purifies expectation values rather than directly output states. In Clifford-only settings under Pauli symmetry, the corresponding detection and correction criteria are proved optimal for the analyzed Hamiltonian-simulation scenario [2503.13114].

Taken together, these developments define quantum subspace verification as a family of verification tasks rather than a single protocol. The common core is the certification of membership in a known subspace or sector under restricted measurements. The main divergences concern what is being verified—states, code spaces, entangled subspaces, constraint subspaces, computational subspaces, temporal trajectories, or symmetric channels—and what follows from successful verification: statistical acceptance, fidelity bounds, rejection of bad shots, local recovery, reduced-space simulation, or channel purification. The strongest current results are concentrated in structured settings, especially stabilizer and CSS code spaces, whereas entangled-subspace verification and local-verifiability obstructions indicate that no comparably uniform local theory exists for arbitrary target subspaces [2409.19699][2412.19540][2509.01281].

Source: https://www.emergentmind.com/topics/quantum-subspace-verification