Threshold theorem: below p_T, concatenated circuits simulate C to distance ε with polylog overhead
OpenQAlgorithms.GottesmanIntro.threshold_theoremThere exist an code with fault-tolerant gadgets of a common depth (Definitions 2–5) and a threshold with the following property. For every with there are constants and such that for every ideal circuit (Definition 6: preparation, , , CNOT, wait and measurement locations, no input) and every there is a level for which the -fold concatenated fault-tolerant simulation satisfies
- and , where is the number of locations of ;
- for every local stochastic error model on (Definition 12) with error bounds , ,
The output of is decoded classically level by level. The threshold theorem shows that arbitrarily long quantum computations can be performed reliably with polylogarithmic overhead once the physical error rate is below a constant.
Formalization Note The fault-tolerant circuit is pinned to the concatenated circuit of the proof (p.42); an unconstrained circuit with free classical computation could sample the output of classically. is encoded as , which is positive also for .
import Definitions.Def_QAlgorithms_Defs_FTCircuit import Definitions.Def_QAlgorithms_GottesmanIntro_Ch05b_Model set_option autoImplicit false
namespace QAlgorithms.GottesmanIntro
/-- Gottesman intro p.41, Theorem 10 (threshold theorem): "There is a threshold error rate `p_T`.
Suppose we have a local stochastic error model with `p_i ≤ p < p_T`. Then for any ideal circuit
`C`, and any `ε > 0`, there exists a fault-tolerant circuit `C′` which, when it undergoes the
error model, produces an output which has statistical distance at most `ε` from the output of
`C`. `C′` has a number of qubits and a number of timesteps which are at most `polylog(|C|/ε)`
times bigger than the number of qubits and timesteps in `C`, where `|C|` is the number of
locations in `C`."
There are a fault-tolerant protocol `P` (an `[[n, 1, 2t + 1]]` code with gadgets of a common
depth) and a threshold `p_T > 0` such that for every `0 ≤ p < p_T` there are constants `K > 0` and
`c` with: for every ideal circuit `C` and every `ε > 0` there is a concatenation level `L` such
that `C′ = C^{(L)}` (the `L`-fold fault-tolerant simulation, p.42) has at most
`K (log₂(2 + |C|/ε))^c` times as many qubits and time steps as `C`, and for every local
stochastic error model on `C′` with bounds `p_i ≤ p` (`p_i < 1`, Def. 12) the decoded output
distribution of `C′` is within total variation distance `ε` of the output distribution of `C`. -/
theorem threshold_theorem :
∃ P : FTProtocol, P.IsFaultTolerant ∧ P.HasUniformDepth ∧
∃ pT : ℝ, 0 < pT ∧ ∀ p : ℝ, 0 ≤ p → p < pT →
∃ (K : ℝ) (c : ℕ), 0 < K ∧
∀ (N : ℕ) (C : LocCircuit N), C.IsIdealCircuit → ∀ ε : ℝ, 0 < ε →
∃ L : ℕ,
((concatCircuit P L C).numQubits : ℝ) ≤
K * Real.logb 2 (2 + (C.size : ℝ) / ε) ^ c * (C.numQubits : ℝ) ∧
((concatCircuit P L C).numSteps : ℝ) ≤
K * Real.logb 2 (2 + (C.size : ℝ) / ε) ^ c * (C.numSteps : ℝ) ∧
∀ (M : LocalStochasticModel (N * (P.n + P.A) ^ L)) (b : ℕ × ℕ → ℝ),
(∀ ℓ ∈ (concatCircuit P L C).locs, b ℓ < 1) →
(∀ ℓ ∈ (concatCircuit P L C).locs, b ℓ ≤ p) →
IsLocalStochastic (concatCircuit P L C) M b →
tvDist (pushDist (concatDecode P L C) (noisyOutDist (concatCircuit P L C) M))
C.idealOutDist ≤ ε := by
sorry
end QAlgorithms.GottesmanIntro
Read-back
What the Lean code literally says, in plain math · claude-opus-5-5
Notions used.
A protocol consists of:
- numbers ;
- a code ;
- an encoder and a decoder ;
- gadget circuits for the location kinds prep, , , CNOT, meas and wait;
- an EC circuit ;
- a classical .
" is fault tolerant" means all of the following:
- , and all operators on qubits are detectable on .
- is an isometry onto .
- is a Kraus channel decoding to for all in the span of weight- Pauli strings.
- All gadgets are well formed with the stated live data wires.
- The Meas, Prep A/B, Gate A/B and EC A/B superoperator identities hold under Pauli faults. Each is stated with the filters (projection onto Pauli-weight- errors of ) whenever fault count plus filter parameters is .
"Uniform depth" means all gadgets have the same number of time steps.
Concatenation. is the level- concatenation, on wires: and . replaces each step by a gadget phase followed by an EC phase on the blocks of non-measured wires. The decoder applies the one-level logical record decoding times, from the top level down. For a circuit:
- is its wire count;
- is its number of time steps;
- is its number of locations.
Distributions. For a local stochastic model (described below):
The ideal distribution is the output distribution with no errors and a trivial one-dimensional environment. The push-forward is , and the total variation distance is
A model is local stochastic for a circuit with bound when:
- its initial environment is a density matrix;
- its weights are nonnegative, vanish off subsets of the locations, and sum to over those subsets;
- each with is a channel on the location's wires environment;
- for every .
Statement. There exists a fault-tolerant protocol with uniform depth, and a real , such that the following holds.
For every with there exist a real and , depending only on , with the following property. For every , every ideal circuit on wires (well formed from and to no live wires, all conditions constantly true) and every , there is a level with:
Moreover, for every local stochastic model on wires with any bound satisfying and at every location of :
Remarks on edge cases.
- The qubit count of is by definition . For the qubit inequality reads .
- The constant inside the logarithm makes the logarithm at least .
- is permitted.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.