Imagine you're a mechanic listening to an engine through a stethoscope — you hear one combined rumble, but you need to know which cylinder is misfiring and whether it's a timing issue or a fuel issue. That's the problem with randomized benchmarking (RB). It gives you one number — a decay rate — that summarizes the aggregate noisiness of your quantum gates. But when different gates have different errors (which they always do), that single number becomes maddeningly opaque. This paper hands you the decomposition. The committed claim: for any low-error gate set forming a unitary 2-design, the RB decay parameter splits cleanly into a linear contribution from Pauli-stochastic errors and a gauge-invariant quadratic contribution from Hamiltonian (coherent) errors. Pauli-correlation and active generators only show up at higher order. This isn't just a perturbative expansion — it's an interpretability theorem. The decay rate is now a readable function of physically meaningful error coordinates. The method works by factorizing each noisy gate as the ideal gate followed by exp(small error generator), then perturbing the so-called L-matrix that governs the RB signal decay. The perturbation theory is in the error generator magnitude — think of it as a Taylor expansion where the small parameter is how noisy your gates are. For the workhorse single-qubit gate set {Xπ/2, Yπ/2}, the authors go all the way to an explicit formula in three first-order gauge-invariant (FOGI) coordinates: two intrinsic over-rotations and one relational axis-tilt error. These are the physically tunable knobs a hardware engineer actually cares about. The payoff is immediate for two RB variants that have confused practitioners. Direct RB (DRB) uses random circuit layers rather than Clifford group elements, and its decay rate depends on the sampling distribution over those layers — the paper's formula explains exactly why. Compiled RB (CRB), which compiles random Clifford gates into native gates, has a compilation-dependent decay rate — again, the formula shows how different compilation strategies weight the underlying FOGI errors differently. These weren't mysterious before; people knew the dependence existed. But now there's an equation, not a hand-wave. The architecture here is pure analytical perturbation theory on representation-theoretic objects — no numerics, no simulation, no hardware. The L-matrix lives in the space of superoperators, and the key insight is that the RB protocol's averaging over random sequences projects this matrix onto a subspace where the perturbation theory becomes tractable. The gauge-invariance of the result is critical: it means the formula doesn't depend on arbitrary decomposition choices that plague other error characterizations. Integrity is strong for a theory paper — the results are mathematical derivations, not empirical claims. The main vulnerability is the perturbative regime assumption: errors must be small enough for the linear+quadratic truncation to be accurate. The paper doesn't numerically benchmark how small is small enough, which is the obvious next step. There's no code, no simulation, no experimental validation — this is a theorem paper, and the proofs are the validation. The practical import is that hardware teams running RB on their processors can now back out which physical error mechanisms dominate from the decay rates they already collect, especially if they run DRB with multiple sampling distributions or CRB with multiple compilers. The formula converts a diagnostic ritual into a diagnostic instrument.
