---
title: Optimized PEC+QED Hybrid Mitigation
url: https://www.emergentmind.com/topics/optimized-pec-qed
type: topic
---

# Optimized PEC+QED Hybrid Mitigation

Searching arXiv for recent papers on probabilistic error cancellation combined with quantum error detection.
arxiv_search(query="probabilistic error cancellation quantum error detection PEC QED arXiv", max_results=10, sort_by="submittedDate")
Searching arXiv for recent papers on optimized PEC+QED and related co-design approaches.
Optimized PEC+QED denotes a class of hybrid error-mitigation protocols that combine probabilistic error cancellation (PEC) with quantum error detection (QED) and post-selection. In this setting, QED filters away detectable fault branches, producing an accepted logical channel that is weaker than the original physical noise, while PEC is applied only to the residual accepted noise rather than to the full channel [2605.12149]. Recent work has treated this as a co-design problem over code choice, detection spacing, truncation order, noise characterization, and symmetry-measurement configuration, because the reduction in PEC sampling overhead must be balanced against rejection and measurement costs introduced by QED [2604.19871][2607.01072].

## 1. Protocol architecture and operational setting

The canonical feedback-free QED+PEC protocol encodes \(n-2\) logical qubits into \(n\) physical qubits using the \([n,n-2,2]\) Iceberg code, with stabilizers
\[
S_X=X^{\otimes n},\qquad S_Z=Z^{\otimes n},\qquad \Pi=\tfrac12(I+S_X)\,\tfrac12(I+S_Z).
\]
The logical circuit is divided into \(M=L/T\) blocks, each containing \(T\) layers of noisy Clifford gates. After each block, all stabilizers are measured and the run is post-selected on the \(+1\) outcome; failed shots are discarded. On the surviving branch, one obtains a reduced logical channel \(\Lambda_{\rm acc}\), and PEC is then applied to undo that channel up to truncated order \(K\). No mid-circuit recovery is required [2605.12149].

A complementary formulation uses the Sparse-Pauli-Lindblad noise model. A noisy layer is written as
\[
\Lambda(\rho)=\bigcirc_{k\in\mathcal K}\bigl(w_k\,\rho +(1-w_k)\,E_k\,\rho\,E_k\bigr),
\]
with inverse quasi-channel
\[
\Lambda^{-1}(\rho)=\gamma\;\bigcirc_{k\in\mathcal K}\bigl(w_k\,\rho -(1-w_k)\,E_k\,\rho\,E_k\bigr),\qquad
\gamma=\prod_{k\in\mathcal K}(2w_k-1)^{-1}.
\]
Within PEC+QED, symmetry measurements are used to detect a subset of errors, and PEC is restricted to the undetectable subset. In the ideal-check limit, the total one-layer sampling cost becomes
\[
C_{\rm PEC+QED}=\frac{(\gamma^{\rm undet})^2}{p_{\rm no\text{-}det}},
\]
so the principal optimization target is the trade-off between shrinking \(\gamma^{\rm undet}\) and maintaining a favorable acceptance probability [2607.01072].

## 2. Accepted logical channel and perturbative inverse construction

Let \(\mathcal K\) index all elementary Pauli fault locations in one block, each with probability \(w_i\). Propagating \(P_i\) to the block end gives \(\widetilde P_i\). For a subset \(I\subseteq\mathcal K\),
\[
w_I=\prod_{i\in I}w_i,\qquad
\pi_I=w_I\prod_{j\not\in I}(1-w_j),\qquad
\mathcal P_I(\rho)=\widetilde P_I\,\rho\,\widetilde P_I^\dagger.
\]
The subset survives post-selection iff \([\widetilde P_I,S_X]=[\widetilde P_I,S_Z]=0\). Defining
\[
\mathcal K_{\rm pass}=\{\,I:[\widetilde P_I,S_X]=[\widetilde P_I,S_Z]=0\},
\]
the exact unnormalized accepted channel and success probability are
\[
\overline\Lambda(\rho)=\sum_{I\in\mathcal K_{\rm pass}}\pi_I\,\mathcal P_I(\rho),\qquad
p_{\rm succ}=\sum_{I\in\mathcal K_{\rm pass}}\pi_I,
\]
and the normalized accepted logical channel is
\[
\Lambda_{\rm acc}(\rho)=\Lambda_{\rm res}(\rho):=\frac{\overline\Lambda(\rho)}{p_{\rm succ}}.
\]
\(\Lambda_{\rm res}\) is the residual channel that PEC must approximate and invert [2605.12149].

The key complication is that post-selection correlates accepted fault branches through stabilizer-commutation constraints, so the sparse Pauli-Lindblad factorization underlying bare PEC no longer applies directly. The inverse is therefore constructed perturbatively. With total noise weight
\[
W=\sum_{i\in\mathcal K}w_i=O(c\,n\,T),
\]
one truncates the Bernoulli expansion to subsets \(|I|\le K\), forms the truncated accepted channel \(\hat\Lambda_{\rm acc}=\mathrm{id}+R_K\), and then uses the degree-\(K\) Neumann inverse
\[
\Lambda_{\rm inv}^{(K)}=\operatorname{Trunc}_{\le K}\!\Bigl[\sum_{r=0}^K(-1)^r\,R_K^{\circ r}\Bigr]
=\sum_P c_P^{(K)}\,\mathcal P .
\]
Sampling from \(\{c_P^{(K)}\}\) implements \(\Lambda_{\rm inv}^{(K)}\) up to \(O(W^{K+1})\) error. Only subsets of size \(\le K\) are enumerated, so classical preprocessing scales as
\[
\sum_{r=0}^K\binom{m}{r}=O(m^K)
\]
with \(m=|\mathcal K|\), and each branch needs \(O(n)\) symplectic updates, giving \(O(m^K\,n)\) cost per block rather than \(2^m\). For \(K\ge1\),
\[
\gamma_K=1+O(W),\qquad p_{K,\rm succ}=1-O(W),\qquad C_K=\frac{\gamma_K^2}{p_{K,\rm succ}}=1+O(W),
\]
while the one-block error satisfies
\[
\|\Lambda_{\rm inv}^{(K)}\circ\Lambda_{\rm acc}-\mathrm{id}\|_\diamond=O(W^{K+1}).
\]
Over multiple blocks, the error accumulates at most additively and the variance factor multiplies blockwise [2605.12149].

## 3. Optimization variables: detection interval, truncation order, and the discrete-Zeno trade-off

A central optimization variable is the QED interval: how often detection cycles are inserted. One efficiency condition is obtained by dividing a circuit of total depth \(M\) into \(N=M/L\) mitigable units of \(L\) logical layers. If \(\epsilon_p\) is the per-layer physical error, \(\epsilon_\ell(L)\) is the post-selected logical error per unit, and \(p_{rej}(L)\) is the rejection probability, break-even requires
\[
4(\epsilon_p-\epsilon_\ell(L))>\frac{p_{rej}(L)}{L},
\qquad\text{equivalently}\qquad
\frac{\epsilon_p-\epsilon_\ell(L)}{p_{rej}(L)}>\frac14.
\]
This criterion makes the interval itself an architectural knob rather than a fixed code-level choice [2604.19871].

In explicit circuit-level studies, the canonical strategy \(L=1\) does not break even in the tested codes. For depolarizing two-qubit error \(p_2=10^{-3}\) with single-qubit \(p_1=p_2/10\), Stim-based simulations gave
\[
(\epsilon_p-\epsilon_\ell(1))/p_{rej}(1)<0.25
\]
for the \([[6,4,2]]\) and \([[12,10,2]]\) Iceberg codes and the \([[9,1,3]]\) surface code. For the \([[6,4,2]]\) Iceberg code, \(L_{\rm thresh}=4\) and \(L_{\rm opt}\sim345\); for \([[12,10,2]]\), break-even begins at \(L_{\rm thresh}=2\); for the \([[9,1,3]]\) surface code, break-even was not observed at any \(L\) up to hundreds of layers [2604.19871].

A different optimization picture appears in the ideal-stabilizer-measurement analysis of the \([n,n-2,2]\) Iceberg code. There, a toy model gives
\[
\ln C_{\rm QED+PEC}
=\gamma T\!\Bigl(1+\tfrac4N\Bigr)
+\gamma^2 T\tau\!\Bigl(\tfrac2N-\tfrac7{N^2}\Bigr)
+O((\gamma\tau)^2),
\]
so for \(N>4\), making \(\tau\) smaller reduces \(C\). The paper therefore recommends the smallest practical detection interval \(T=1\) when syndrome extraction is cheap, and first order \(K=1\) as the minimal nontrivial choice for distance-2 codes. Higher \(K\) yields residual error \(O(W^{K+1})\) per block but increases preprocessing as \(O(m^K)\) and the quasiprobability norm roughly as \((1+\|W\|)^K\) [2605.12149]. This suggests interval optimality is architecture- and noise-model-dependent: cheap, idealized detection and explicit syndrome-extraction overhead lead to different optima.

## 4. First-cycle transients and steady-state extraction

A major obstacle to naive PEC+QED is that repeated QED cycles violate the position-independent noise assumption used in standard PEC. Transition-matrix analysis shows that the post-selected evolution contains one fast mode of order \(O(p)\) that decays in \(O(1)\) cycles and slow modes of order \(1-O(p)\) that govern steady-state behavior. Physically, the first QED cycle injects leakage states at \(O(p)\) but cannot remove any, whereas later cycles reach a balance of detection and injection [2604.19871].

If one characterizes a gate by tomography of “prepare + gate + QED” without removing this first-cycle transient, the learned model embeds that fast leakage component. The result is that naive PEC+QED can degrade accuracy below the QED-only baseline. The proposed remedy is steady-state extraction (SSE), defined through
\[
S_{G,\mathrm{steady}}=S_G^{\mathrm{concat}}\circ S_{\mathrm{init}}^{-1},
\]
where \(S_{\mathrm{init}}\) is the superoperator of “prepare logical + first QED cycle” and \(S_G^{\mathrm{concat}}\) is the concatenated superoperator “\(S_{\mathrm{init}}\circ\) (ideal \(G\) plus its QED).” Operationally, one tomographically learns \(S_{\mathrm{init}}\), learns \(S_G^{\mathrm{concat}}\) for each logical gate, and computes \(S_{G,\mathrm{steady}}\) by matrix inversion [2604.19871].

On the \([[4,2,2]]\) code, Hamiltonian-level QuTiP simulations showed that SSE alone reduces per-cycle prediction error by \(\sim39.7\times\) versus naive single-cycle tomography, and in end-to-end PEC studies reduces the observable bias by up to \(10.2\times\) below the QED-only baseline for depolarizing, dephasing, and amplitude-damping noise. Within optimized PEC+QED, SSE therefore functions as a characterization protocol that isolates the steady-state channel actually seen by long runs, rather than as an additional mitigation layer [2604.19871].

## 5. Optimization over symmetry measurements

When QED is implemented by measuring symmetries, the choice of which symmetries to measure becomes a classical optimization problem. Each candidate symmetry \(S_i\) detects a subset \(D_i\) of error generators and incurs mitigation overhead \(\mu_i\) from the circuit that measures it. With error weights \(\lambda_k\), one objective is
\[
\min_{\mathcal C\subseteq\mathcal S}
\underbrace{\sum_{k\notin\mathcal D(\mathcal C)}\lambda_k}_{\text{undetected errors}}
+
\underbrace{\sum_{i\in\mathcal C}\mu_i}_{\text{symmetry circuits}},
\qquad
\mathcal D(\mathcal C)=\bigcup_{i\in\mathcal C}D_i .
\]
This is a weighted set-cover problem with element penalty. A greedy approximation is used to generate a candidate pool, followed by subset enumeration once idling-aware costs are included [2607.01072].

Noisy symmetry measurements introduce two distinct penalties. Errors during the check can create false detections, increasing variance without reducing bias, while undetected double-errors from pairs of detectable errors can survive post-selection and increase bias. The framework therefore discourages symmetries with large \(\mu_i\) unless they cover sufficiently large \(\lambda_k\)-weight. This is the sense in which the method optimizes QED for PEC, rather than simply appending post-selection to an existing symmetry set [2607.01072].

For GHZ-state preparation, the full stabilizer group contains \(n-1\) \(ZZ\) generators and the global \(X^{\otimes n}\). Because measuring \(X^{\otimes n}\) costs \(n\) CNOTs, the optimizer excludes that global check and instead selects nonlocal weight-2 stabilizers. For \(n=50\), a 4-element configuration detects \(\sim60\%\) of the total error weight with only \(\sim4\) CNOTs of check cost. At fixed shot count \(10^6\), optimized PEC+QED yields up to \(10\times\) lower total-square-error than PEC for \(n\sim50\). For the generalized superfast encoded Fermi-Hubbard model, PEC+QED improves observable estimation on a \(2\times2\) lattice, and for larger systems mitigation overheads can be reduced by measuring only subsets of stabilizers; for \(m=5,6\), the optimal subsets omit \(17\)–\(22\%\) of stabilizers [2607.01072].

## 6. Benchmark regimes, empirical performance, and relation to physical-layer hybrids

The current literature reports gains in three distinct regimes. In logical GHZ-state preparation with the \([n,n-2,2]\) Iceberg code under circuit-level depolarizing noise and ideal stabilizer measurements, first-order QED+PEC reaches \(n=200\) physical qubits and lowers sampling overhead by three to four orders of magnitude relative to standard PEC while maintaining \(F\simeq0.956\). In the same setting, the bare PEC sample-variance factor is \(\sim10^7\), whereas QED+PEC gives \(\sim10^3\) [2605.12149].

In a more hardware-oriented co-design study, a \([[6,4,2]]\) Iceberg code was used for 4-vertex MaxCut QAOA at depths \(p=4,6,8\) under depolarizing noise \(p_2=10^{-3}\), \(p_1=p_2/10\), with two QED checks per circuit and \(10^5\) attempted shots. PEC+QED achieved \(2\)–\(11\times\) lower absolute error and up to \(31\times\) lower mean-squared error versus PEC on physical qubits. Across \(p=4,6,8\), the reported quasiprobability norms were \(1.053,1.080,1.108\) for PEC only and \(1.005,1.006,1.006\) for PEC+QED, with acceptance rates between \(90.2\%\) and \(93.9\%\) [2604.19871].

A related line of work pushes PEC and QED to the physical layer inside a logical code without modifying the decoder. In that framework, any linear QEM method can be integrated into the physical layer because QEC is itself a linear quantum map. For code-capacity memory, a combined PEC+QED construction cancels all error weights \(k\le \omega+1\), so the leading logical failure becomes \(O(p^{\omega+2})\), which is described as increasing the effective code distance by \(2\). Simulations on repetition and rotated surface codes showed that a distance-3 code with physical-level PEC achieves logical error rates lower than or similar to a distance-5 unmitigated code while using \(40\%\) and \(64\%\) fewer qubits, respectively [2601.18384].

Two recurrent limitations define the practical boundary of optimized PEC+QED. First, noisy symmetry or stabilizer measurements can negate the advantage: readout-only flips mainly increase post-selection cost, while noisy GHZ-assisted global stabilizer extraction can remove the advantage entirely [2605.12149]. Second, optimization is code- and circuit-specific: high-rate distance-2 Iceberg codes and carefully spaced detection cycles can satisfy the relevant cost conditions, whereas low-rate codes may not [2604.19871]. A plausible implication is that optimized PEC+QED is best viewed not as a single protocol, but as a family of co-designed mitigation architectures in which post-selection reshapes the effective channel that PEC must invert.

Source: https://www.emergentmind.com/topics/optimized-pec-qed