---
title: Leakage-Aware Randomized Benchmarking
url: https://www.emergentmind.com/topics/leakage-aware-randomized-benchmarking-rb
type: topic
---

# Leakage-Aware Randomized Benchmarking

Leakage-aware randomized benchmarking (RB) comprises randomized benchmarking protocols that explicitly treat population transfer outside a designated computational subspace rather than assuming that all errors remain internal to that subspace. In the cited literature, leakage is modeled by embedding the computational space in a larger Hilbert space, modifying the RB experiment and its fitting model so that decay parameters can estimate survival, leakage, seepage, and subspace fidelity even when standard RB would misinterpret the data. Across the main protocol families, the central technical theme is that leakage changes the sequence-length dependence from the standard single-exponential form to single-, bi-, multi-, or matrix-exponential decays, with the precise form determined by the subspace structure, the gate set, and the noise model [1412.4126, 1704.03081, 2010.07974, 2502.00154].

## 1. Subspace structure and leakage metrics

Leakage-aware RB starts from an explicit decomposition of the physical Hilbert space into a computational part and a leakage part. The notation varies across papers, but the structure is the same: \(\mathcal{H} = \mathcal{H}_1 \oplus \mathcal{H}_2\), \(X = X_1 \oplus X_2\), or \(\mathcal{H} = \mathcal{H}_C \oplus \mathcal{H}_L\). Leakage is the unwanted transfer of population from the computational subspace into the leakage subspace, while seepage denotes population flowing back into the computational subspace [1412.4126, 1704.03081, 2502.00154].

A basic quantity is the survival rate of a state \(\rho\) under a channel \(\mathcal{E}\) with respect to \(\mathcal{H}_1\),
\[
s(\rho|\mathcal{E},\mathcal{H}_1) =
\frac{\mathrm{Tr}[\mathcal{P}_{\mathcal{H}_1}\mathcal{E}(\rho)]}{\mathrm{Tr}[\rho]}.
\]
Averaging over states in \(\mathcal{H}_1\) gives
\[
s(\mathcal{E},\mathcal{H}_1)=
\mathrm{Tr}[\mathcal{P}_{\mathcal{H}_1}\mathcal{E}(d_{\mathcal{H}_1}^{-1}\mathcal{P}_{\mathcal{H}_1})],
\]
and for the whole space this average survival is denoted \(s_{\rm inc.}(\mathcal{E})\). In the same framework, incoherent leakage is
\[
l_{\rm inc.}(\mathcal{E}) = 1-s_{\rm inc.}(\mathcal{E}),
\]
while coherent leakage is
\[
l_{\rm coh.}(\mathcal{E}) = s_{\rm inc.}(\mathcal{E}) - s_{\rm coh.}(\mathcal{E}),
\]
with
\[
s_{\rm coh.}(\mathcal{E}) =
\mathrm{Tr}\left[\mathcal{I}_{\mathcal{H}_1}\mathcal{E}\left(\frac{1}{d_{\mathcal{H}_1}}\mathcal{I}_{\mathcal{H}_1}\right)\right]
+
\mathrm{Tr}\left[\mathcal{I}_{\mathcal{H}_2}\mathcal{E}\left(\frac{1}{d_{\mathcal{H}_2}}\mathcal{I}_{\mathcal{H}_2}\right)\right].
\]
These definitions separate irreversible population loss from reversible population exchange [1412.4126].

A complementary formulation introduces state leakage
\[
L(\rho)=\operatorname{Tr}[{}_2\rho]=1-\operatorname{Tr}[{}_1\rho],
\]
together with channel leakage and seepage rates
\[
L_1(E)=L\!\left(E\left(\frac{{}_1}{d_1}\right)\right), \qquad
L_2(E)=1-L\!\left(E\left(\frac{{}_2}{d_2}\right)\right).
\]
In this notation, \(L_1\) is the average population lost from \(X_1\) and \(L_2\) is the average population returning from \(X_2\) [1704.03081].

The same 2017 framework also defines coherent leakage through off-diagonal blocks between computational and leakage sectors:
\[
CL(\rho)=\|P_C(\rho)\|_1=\|\rho-P_I(\rho)\|_1,
\]
with the bound
\[
CL(\rho)\le 2\sqrt{p_l(1-p_l)},
\]
where \(p_l=L(\rho)\). Channel-level coherent leakage and seepage obey analogous bounds,
\[
CL_1(E)\le 2\sqrt{L_1(1-L_1)}, \qquad
CL_2(E)\le 2\sqrt{L_2(1-L_2)}.
\]
This makes explicit that population transfer and coherence between subspaces are distinct aspects of leakage [1704.03081].

A more recent parameterization emphasizes the computational sector directly. For a channel \(\Lambda\), the computational population is
\[
t[\Lambda]=\frac{1}{d_C}\mathrm{Tr}[\mathbb{I}_C\Lambda(\mathbb{I}_C)],
\]
the leakage rate is \(\tau=1-t\), the depolarizing parameter is
\[
r[\Lambda] := \frac{1}{d_C^2-1}\sum_{i=1}^{d_C^2-1}\mathrm{Tr}[P_{C,i}\Lambda(P_{C,i})],
\]
and the computational error is \(\lambda=t-r\). The average gate fidelity on computational inputs is then
\[
F[\Lambda]=\frac{(d_C-1)r[\Lambda]+t[\Lambda]}{d_C}.
\]
This formulation is designed for fidelity estimation in the presence of leakage rather than leakage-rate estimation alone [2502.00154].

## 2. Modified decay laws and leakage signatures

A defining result of leakage-aware RB is that leakage alters the functional form of the RB decay. In an extended qutrit model, leakage errors lead to a sum of exponentials rather than a single exponential. With initial state and measurement restricted to the qubit subspace, and with an expanded single-qubit Clifford group that includes both \(\pm\pi/2\) and \(\pm 3\pi/2\) rotations, the average fidelity decay takes the form
\[
F(m)=C_1\alpha^m + C_2\mathcal{A}_{55}^m + C_3.
\]
For unital leakage errors, \(\mathcal{A}_{51}=0\), so the asymptotic fidelity is \(F(\infty)=1/3\), whereas in the absence of leakage it saturates at \(1/2\). The change in asymptote is therefore a diagnostic of leakage in that model [1308.2928].

A scalable protocol for explicitly estimating leakage rates was developed by modifying the RB randomization so that gates are sampled from a group
\[
\mathcal{G}=\{g=v\oplus(\pm w)\},
\]
where \(v\) and \(w\) are unitary 1-designs on \(\mathcal{H}_1\) and \(\mathcal{H}_2\), respectively. The experiment prepares \(|0\rangle\) in \(\mathcal{H}_1\), applies a random sequence, and measures the probability of being in \(\mathcal{H}_1\) after \(m\) steps. The resulting averaged decay depends on the type of leakage [1412.4126].

For incoherent leakage, with \(\mathcal{H}_1=\mathcal{H}\), the sequence average follows
\[
\mathds{E}_{\mathbf{k}}p_{\mathbf{k}}
=
A\, s_{\rm inc.}^{m-1}(\mathcal{E}) + O(m\epsilon),
\]
where \(A\) is SPAM-dependent and \(\epsilon\) quantifies gate dependence. For coherent leakage, with \(\mathcal{H}_1\subsetneq \mathcal{H}\), the decay is bi-exponential,
\[
\mathds{E}_{\mathbf{k}}p_{\mathbf{k}}
=
B\lambda_+^{m-1}+C\lambda_-^{m-1}+O(m\epsilon),
\]
and
\[
s_{\rm coh.}(\mathcal{E})=\lambda_+ + \lambda_-.
\]
The distinction between a single exponential and a two-exponential decay is the operational signature that separates incoherent leakage from coherent population transfer [1412.4126].

This directly contradicts a common simplification in which leakage is treated as a small correction to ordinary RB. The cited works instead show that the structure of the decay itself can be altered, including non-scalar matrix-exponential behavior in more general settings [1308.2928, 2010.07974].

## 3. Clifford-based leakage benchmarking protocols

One major line of work modifies standard Clifford RB so that leakage and seepage can be estimated together with fidelity. In leakage randomized benchmarking (LRB) for a Clifford gate set, the protocol prepares \(\rho_0=|0\rangle\langle 0|\) in \(X_1\), applies a random sequence of \(m\) Clifford gates together with a final recovery Clifford, and measures computational-subspace populations. The leakage-population decay is fit to
\[
p_{_1}(m)=A+B\lambda_1^m,
\]
with
\[
L_1(E)=(1-A)(1-\lambda_1), \qquad
L_2(E)=A(1-\lambda_1).
\]
When leakage is weak, the approximations
\[
A \approx \frac{L_2}{L_1+L_2}, \qquad
B \approx \frac{L_1}{L_1+L_2}, \qquad
\lambda_1=1-L_1-L_2
\]
connect the fit directly to leakage and seepage. A second fit,
\[
p_0(m)=A_0+B_0\lambda_1^m + C_0\lambda_2^m,
\]
yields the average gate fidelity
\[
\overline{F}=\frac{1}{d_1}\big[(d_1-1)\lambda_2 + 1 - L_1\big].
\]
This protocol is intended to estimate \(L_1\), \(L_2\), and \(\overline{F}\) simultaneously for a Clifford gate set [1704.03081].

A different but related construction extends the Clifford group itself to randomize phases on leakage levels:
\[
\mathcal{C}_N^\ast = \mathcal{C}_N \times \left\{{}_2 \oplus (\pm1)\right\}^{\otimes N}.
\]
The purpose of this extension is to destroy residual coherence in the non-computational subspace so that the effective dynamics reduce to populations on a \(3^N\)-dimensional diagonal subspace. After twirling, the error channel becomes a real, non-negative \(3^N\times 3^N\) Markov matrix, and the sequence fidelity has the generic form
\[
\Phi_y=\sum_i a_i\lambda_i^y.
\]
This formulation admits arbitrary gate-dependence in the error channels and does not require perturbative gate-independence assumptions [1505.00580].

The same protocol analyzes the variance of sequence fidelities and states that, to leading order,
\[
\mathrm{Var}(\Phi_y)\propto y e^{-\kappa y}.
\]
It is also compatible with interleaved randomized benchmarking and can be expanded to benchmarking of arbitrary gates. For an interleaved gate \(V\), the extracted gate error is
\[
\varepsilon_V=\varepsilon_{\mathcal{C}^\ast\times V}-\varepsilon_{\mathcal{C}^\ast}.
\]
The setting is described as relevant for superconducting transmon qubits, among other systems [1505.00580].

Taken together, these Clifford-based approaches establish two complementary viewpoints: one emphasizes direct estimation of leakage and seepage rates together with subspace fidelity, while the other emphasizes a complete RB protocol whose effective twirl converts leakage dynamics into a classical Markov problem [1704.03081, 1505.00580].

## 4. General RB theory and representation-theoretic extensions

A general framework for randomized benchmarking later placed leakage-aware protocols inside a broader theory of RB. In this framework, the implementation map \(\phi\) that associates a gate to a quantum channel is not restricted to the computational subspace and may be trace-non-increasing or non-trace-preserving, thereby encompassing leakage and loss. The core structural statement is that, under realistic conditions, the RB signal is well described by
\[
p(m)\approx \sum_{\lambda\in\Lambda}\mathrm{Tr}(A_\lambda M_\lambda^m),
\]
a linear combination of matrix exponential decays. In leakage-aware scenarios, additional subspaces appear, so there can be multiple decay channels associated with leakage-space population and return processes [2010.07974].

The same framework states that accurate interpretation of these decay rates requires Markovianity and time-independence, closeness on average in diamond norm to a reference representation,
\[
\frac{1}{|\mathcal{G}|}\sum_{g\in\mathcal{G}}\|\omega(g)-\phi(g)\|_\diamond \le \delta,
\]
with \(\delta \le 1/9\), and uniform or well-controlled sampling over the group. When these conditions hold, dominant decay rates can be used as figures of merit even in the presence of leakage [2010.07974].

A notable technical contribution of this framework is post-processing that isolates decay components without requiring physical inversion gates. For a subrepresentation \(\lambda\), the filtered data are
\[
k_\lambda(m):=\frac{1}{|G|}\sum_{g\in G}\alpha_\lambda(g)p(m,g),
\]
where \(\alpha_\lambda\) is a filter function. The same paper recommends MUSIC/ESPRIT and related spectral-estimation methods for extracting multiple decay rates, which is particularly relevant when leakage produces matrix-exponential or heavily mixed exponential behavior [2010.07974].

A parallel development extended character randomized benchmarking to non-multiplicity-free groups and derived a new leakage RB protocol for more general groups of gates. The key relaxation is that the group need only preserve the computational-leakage subspace splitting \(\mathcal{H}=\mathcal{H}_1\oplus\mathcal{H}_2\), rather than factor into independently sampled 2-designs on the computational and leakage sectors. In Liouville form, the leakage and seepage rates are
\[
L=\frac{1}{d_1}\langle\langle \mathbb{1}_2|\hat{\Lambda}|\mathbb{1}_1\rangle\rangle, \qquad
S=\frac{1}{d_2}\langle\langle \mathbb{1}_1|\hat{\Lambda}|\mathbb{1}_2\rangle\rangle,
\]
and the subspace average fidelity is
\[
F_{\Lambda,1}=\int d\psi_1 \langle\langle \psi_1|\hat{\Lambda}|\psi_1\rangle\rangle.
\]
For the trivial irrep, the decay is fit in a form that yields
\[
L=(1-B)(1-\lambda), \qquad S=B(1-\lambda),
\]
and more explicitly
\[
S_0(N)=\frac{S}{L+S}+\frac{L}{L+S}(1-L-S)^{N+1}.
\]
Under additional assumptions, the computational-subspace fidelity can be extracted as
\[
F_{\Lambda,1} = \frac{(d_1^2-1)\lambda_{1\perp} + (d_1 + 1)(1-L)}{d_1^2 + d_1}.
\]
This representation-theoretic extension broadens leakage RB to groups that do not satisfy the multiplicity-free restriction of earlier character RB [2011.00007].

## 5. Regime-specific methods and universal-gate-set variants

Recent work has emphasized that there is no single leakage-aware fit model that is optimal across all experimental regimes. One 2025 treatment develops several methods for measuring fidelity with randomized benchmarking in the presence of leakage errors, explicitly organized by error regime and by whether leakage detection or post-selection is available. It also states that, when leakage is substantial, standard randomized benchmarking can cause a large overestimation of fidelity [2502.00154].

Two leakage-detection strategies are identified there: explicitly addressing leakage states and ancilla- or gadget-based leakage detection. The same work gives regime-dependent survival-probability formulas, including short-sequence linearization, an EXP-LIN model when computational error dominates, double-exponential fitting in the no-seepage regime, and basis-randomized fitting for pure population transfer [2502.00154].

| Regime | Method | Representative fit |
|---|---|---|
| Short sequences | Comp. survival | \(\bar{p}(\ell)=1-\ell(1-F)\) |
| Comp. error dominant | EXP-LIN | \(\frac{d_C-1}{d_C}(1-\lambda-\ell\tau)(1-\lambda)^{\ell-1}+\frac{1-\ell\tau}{d_C}\) |
| Comp. error dominant | LPS | \(\frac{d_C-1}{d_C}(1-\lambda)^\ell+\frac{1}{d_C}\) |
| No seepage | 2EXP | \(\frac{d_C-1}{d_C}r^\ell+\frac{1}{d_C}t^\ell\) |
| No seepage | LPS | \(\frac{d_C-1}{d_C}\left(\frac{r}{t}\right)^\ell+\frac{1}{d_C}\) |
| Population transfer | Basis randomization | \(\frac{d_C-1}{d_C}r^\ell+\frac{1}{d_C}\) |

The same source associates retention models with the post-selected cases, namely \(p_\mathrm{retention}(\ell)=1-\ell\tau\) in the computational-error-dominant regime and \(p_\mathrm{retention}(\ell)=t^\ell\) in the no-seepage regime [2502.00154].

A different 2023 development addresses universal gate sets and multi-qubit systems by proposing leakage randomized benchmarking (LRB) and interleaved LRB (iLRB) that require only access to the computational Pauli group and do not require operations that twirl the leakage subspace. The formalism maps leakage dynamics to a classical Markov chain over invariant subspaces. For the noise channel \(\Lambda\), the average leakage and seepage are
\[
L_{ave}(\Lambda)=\operatorname{Tr}[{\Pi}_l\Lambda(\widetilde{\Pi}_c)], \qquad
S_{ave}(\Lambda)=\operatorname{Tr}[{\Pi}_c\Lambda(\widetilde{\Pi}_l)],
\]
and the averaged sequence outcome is
\[
p_{\Pi_c}(m)=((\hat{\Pi}_c|Q^{m-1}|\tilde{\rho}_0)),
\]
where \(Q\) is the transition matrix of the induced Markov process [2304.07884].

In the single-site leakage model, \(Q\) reduces to an \((n+1)\times(n+1)\) matrix,
\[
Q=
\begin{pmatrix}
1-\sum_i p_i & 2q_1 & \cdots & 2q_n \\
p_1 & 1-2q_1 & \cdots & 0 \\
\vdots & \vdots & \ddots & \vdots \\
p_n & 0 & \cdots & 1-2q_n
\end{pmatrix},
\]
while in the cross-talk-free independent leakage model the computational-subspace survival factorizes over qubits. For a target gate benchmarked via iLRB, random Paulis are interleaved with the gate of interest, and for two-qubit examples such as iSWAP and CZ the paper gives analytic relations between fitted decay constants and leakage/seepage rates. This extends leakage-aware benchmarking beyond Clifford-only settings and into generic \(n\)-site gates under stated noise assumptions [2304.07884].

## 6. Numerical validation, practical significance, and limitations

The principal motivation for leakage-aware RB is that leakage can invalidate the interpretation of ordinary RB data. Multiple papers therefore evaluate robustness numerically. In the 2014 survival-rate protocol, simulations with weakly gate-dependent Pauli noise gave \(r^2=0.9991\) for incoherent leakage, and simulations of a three-level shelving system gave \(r^2=0.9904\) for coherent leakage, both in good agreement with the predicted exponential forms [1412.4126].

Earlier analysis of realistic RB error models concluded that, in almost all cases considered, benchmarking provides better than a factor-of-two estimate of average error rate, including leakage to higher levels, provided the leakage-aware sum-of-exponentials model is used when appropriate. The same study stated that only a small number of trials are required for high-confidence estimation if the standard error of the fidelity measurements is small, and reported reliable extraction at modest sampling such as \(K \sim 100\) in the leakage examples [1308.2928].

The 2023 LRB/iLRB framework reported simulations showing robustness to strong depolarizing preparation noise and measurement bit-flip errors, as well as accurate recovery of theoretical leakage and seepage rates for two-qubit gates under realistic device-inspired models. The 2025 regime-based study numerically demonstrated the methods for two-qubit randomized benchmarking and then applied them to previously shared data from Quantinuum systems. In that case, leakage-aware methods generally estimated higher (worse) 2Q infidelity than simple or legacy RB analysis by more than one standard deviation, although the absolute difference remained small; the same analysis noted that the short-sequence method underestimates infidelity when longer sequences are used because the linear approximation breaks down [2304.07884, 2502.00154].

Several limitations recur across the literature. Many derivations assume Markovian and time-independent noise, and some variants additionally assume gate-independence or specific commutation properties of the target-gate noise [2010.07974, 2304.07884]. More general leakage patterns can produce multi-exponential or matrix-exponential signals that are difficult to fit directly, which is why filtering and spectral-estimation methods such as MUSIC/ESPRIT were proposed in the general framework [2010.07974]. Another recurrent point is that \(L_1+L_2\) alone is insufficient for full characterization: thermal relaxation, erasure, and unitary leakage can share the same sum while having different physical behavior and different implications for error correction [1704.03081].

Within these constraints, leakage-aware RB has developed from modified Clifford twirls and survival-rate protocols into a family of subspace-sensitive benchmarking methods. The common conclusion across the cited work is that leakage must be treated as a first-class error mechanism rather than as a perturbation of ordinary RB fits, because it changes both the operational figures of merit and the mathematical form of the benchmarking signal [1412.4126, 1505.00580, 2011.00007, 2502.00154].

Source: https://www.emergentmind.com/topics/leakage-aware-randomized-benchmarking-rb