- Boxes
- definitions
- Ellipses
- theorems and lemmas
- Blue border
- the statement of this result is ready to be formalized; all prerequisites are done
- Orange border
- the statement of this result is not ready to be formalized; the blueprint needs more work
- Blue background
- the proof of this result is ready to be formalized; all prerequisites are done
- Green border
- the statement of this result is formalized
- Green background
- the proof of this result is formalized
- Dark green background
- the proof of this result and all its ancestors are formalized
- Dark green border
- this is in Mathlib
Under the hypotheses of the Fault Tolerance Theorem, for any spacetime fault \(F\) with \(\mathrm{weight}(F) {\lt} d\):
That is, \(F\) either triggers a detector (is detectable) or does not affect logical information (is a stabilizer).
Under the hypotheses of the Fault Tolerance Theorem, for any \(t \in \mathbb {N}\) with \(2t + 1 \leq d\), the code can tolerate \(t\) faults:
In particular, the code can correct up to \(\lfloor (d-1)/2\rfloor \) faults.
For any fault pattern \(p\) with empty syndrome and preserving logical information, there exists a list of generators \(\mathrm{gens}\) such that:
every generator in \(\mathrm{gens}\) has empty syndrome,
every generator in \(\mathrm{gens}\) preserves logical information, and
the \(\mathbb {Z}_2\)-sum of their logical effects is \(0\).
An abstract predicate capturing whether a fault affects the logical information. Given a logical effect predicate \(\ell : \mathrm{SpacetimeFault} \to \mathrm{Prop}\) and a fault \(f\):
This captures changes to the gauging measurement outcome or introduction of logical Pauli errors.
Given base measurement outcomes and a spacetime fault \(f\), the faulted outcomes are defined by:
Time-faults flip the corresponding measurement outcomes.
Given base measurement outcomes and a spacetime fault, the faulted outcomes are obtained by flipping the classical measurement outcome at each location where a time-fault is active:
A Bivariate Bicycle (BB) code with parameters \(\ell , m\) is specified by two polynomials \(A, B \in \mathbb {F}_2[M]\). The structure consists of:
\(\mathrm{polyA} : \mathbb {F}_2[M]\) — the first polynomial \(A\),
\(\mathrm{polyB} : \mathbb {F}_2[M]\) — the second polynomial \(B\).
The X parity check matrix is
an \(\ell m \times 2\ell m\) matrix over \(\mathbb {F}_2\) with rows indexed by \(M\) (X checks) and columns indexed by \(M \oplus M\) (left \(\oplus \) right qubits), formed by horizontally concatenating the matrix representations of \(A\) and \(B\).
The Z parity check matrix is
an \(\ell m \times 2\ell m\) matrix over \(\mathbb {F}_2\) with rows indexed by \(M\) (Z checks) and columns indexed by \(M \oplus M\) (left \(\oplus \) right qubits), formed by horizontally concatenating the matrix representations of \(B^\top \) and \(A^\top \).
The support of a Pauli operator \(X(p, q)\) as a binary vector over \(M \oplus M\) (left \(\oplus \) right qubits) is defined as:
The support of a Pauli operator \(Z(p, q)\) as a binary vector over \(M \oplus M\) (left \(\oplus \) right qubits) is defined as:
The group algebra \(\mathbb {F}_2[x, y]/(x^\ell - 1, y^m - 1)\) is defined as
i.e., the set of finitely-supported functions \(M \to \mathbb {F}_2\) with convolution multiplication.
The monomial set is defined as
This is an additive abelian group where addition represents multiplication of monomials.
The matrix representation of a group algebra element \(p \in \mathbb {F}_2[M]\) is the \(|M| \times |M|\) matrix over \(\mathbb {F}_2\) whose \((\alpha , \beta )\)-entry is \(p(\alpha - \beta )\):
This gives a circulant-like structure on the product group.
The transpose of a group algebra element \(p \in \mathbb {F}_2[M]\) is defined as \(p^\top = p(x^{-1}, y^{-1})\), computed by applying the monomial transpose to the domain:
This maps each monomial \((a, b)\) to \((-a, -b)\) while preserving coefficients.
A structure capturing the full boundary condition convention:
a perfect boundary assumption,
an error correction configuration,
a proof that the configuration provides full protection,
a flag indicating this is a simplification (default: true).
The theoretical convention can be replaced by realistic conditions: given a boundary condition convention and a realistic boundary context, the realistic context is returned unchanged (since in practice, the surrounding computation determines the conditions).
A set of logical operators \(\mathcal{L}\) satisfies the bounded overlap condition with bound \(k\) if for every qubit \(i\), the overlap degree is at most \(k\):
This condition is required to maintain the LDPC property during code deformation.
The characteristic vector map as a type equivalence \(\mathrm{Finset}(\alpha ) \simeq \mathrm{BinaryVector}(\alpha )\), with forward map \(\chi \) and inverse map \(\mathrm{fromBinaryVector}\). Theleft and right inverses are given by the two round-trip properties.
The characteristic vector map is a linear equivalence
with forward map \(\chi = \mathrm{characteristicVector}\), inverse map \(\mathrm{fromBinaryVector}\), additivity \(\chi (S + T) = \chi (S) + \chi (T)\), and scalar compatibility \(\chi (c \cdot S) = c \cdot \chi (S)\).
The three types of checks in the deformed code:
gaussLaw: The Gauss’s law operators \(A_v\) (X-type on vertices and incident edges),
deformed: The deformed stabilizer checks \(\tilde{s}_j\) (from the code deformation),
flux: The flux operators \(B_p\) (Z-type on cycle edges).
The Cheeger constant of a simple graph \(G = (V, E)\) is defined as
where \(\partial S\) is the edge boundary of \(S\). The infimum is taken over all nonempty subsets \(S\) of \(V\) satisfying \(2|S| \leq |V|\).
The circuit phases of the gauging circuit are defined as the inductive type with five constructors:
InitializeEdges: Initialize edge qubits in \(\ket {0}\).
FirstEntangling: Apply the first entangling circuit \(\prod _v \prod _{e \ni v} \mathrm{CX}_{v \to e}\).
MeasureXVertices: Measure \(X_v\) on all vertices.
SecondEntangling: Apply the second entangling circuit (identical to the first).
MeasureZEdges: Measure \(Z_e\) on all edges and discard.
The natural ordering of circuit phases is given by the function \(\mathrm{toNat} : \mathrm{CircuitPhase} \to \mathbb {N}\) defined by:
The cleaned operator \(\bar{L}\) obtained from \(L\) by multiplying by \(\prod _{v \in S} A_v\) is defined by:
where \(\triangle \) denotes the symmetric difference.
A bridge edge configuration consists of:
A finset of bridge connections \(\mathrm{bridges} \subseteq \mathrm{Fin}(W_1) \times \mathrm{Fin}(W_2)\), specifying which pairs of vertices to connect between the two systems.
A proof that at least one bridge edge exists (\(\mathrm{bridges}\) is nonempty).
The Cohen hypergraph \(\mathrm{CohenHypergraph}(W, \mathrm{numChecks})\) is a hypergraph with vertex set \(\mathrm{Fin}(W)\) and hyperedge set \(\mathrm{Fin}(\mathrm{numChecks})\). The vertices represent the \(W\) qubits in \(\mathrm{supp}(L)\), and the hyperedges are the restricted Z-type checks.
The Cohen kernel property for a Cohen hypergraph \(HG\) states that the operator kernel of \(HG\) is exactly \(\{ 0, L\} \):
This is the defining property of the Cohen scheme: the restriction of Z-type checks to \(\mathrm{supp}(L)\) has kernel spanned by the logical operator \(L\) (i.e., \(L\) is irreducible).
The layered Cohen construction builds a hypergraph from the original Cohen hypergraph \(HG\). The incidence function is defined by:
For a chain edge \(\mathrm{chain}(l, i)\): the incidence set is \(\{ (l, i),\, (l+1, i)\} \).
For a hypergraph copy \(\mathrm{hypergraphCopy}(\ell , c)\): the incidence set is the image of \(HG\)’s incidence of check \(c\) under the embedding \(v \mapsto (\ell , v)\).
The layered hyperedge type combines two kinds of edges:
Chain edge \(\mathrm{chain}(l, i)\): connects vertex \((l, i)\) to \((l+1, i)\) for \(l \in \mathrm{Fin}(d)\) and \(i \in \mathrm{Fin}(W)\).
Hypergraph copy \(\mathrm{hypergraphCopy}(\ell , c)\): a copy of check \(c\) in layer \(\ell \), for \(\ell \in \mathrm{Fin}(d+1)\) and \(c \in \mathrm{Fin}(\mathrm{numChecks})\).
The product hyperedge type has three constructors:
\(\mathrm{left}(h)\): a check from the first ancilla system, indexed by \(h \in \mathrm{Fin}(nc_1)\).
\(\mathrm{right}(h)\): a check from the second ancilla system, indexed by \(h \in \mathrm{Fin}(nc_2)\).
\(\mathrm{bridge}(i, j)\): a bridge edge connecting vertex \(i \in \mathrm{Fin}(W_1)\) to vertex \(j \in \mathrm{Fin}(W_2)\).
The product measurement hypergraph combines two hypergraphs \(HG_1\) and \(HG_2\) with bridge edges. Its incidence function is:
For \(\mathrm{left}(h)\): the image of \(HG_1.\mathrm{incidence}(h)\) under \(\mathrm{inl}\).
For \(\mathrm{right}(h)\): the image of \(HG_2.\mathrm{incidence}(h)\) under \(\mathrm{inr}\).
For \(\mathrm{bridge}(i, j)\): the set \(\{ \mathrm{inl}(i), \mathrm{inr}(j)\} \).
Let \(C\) be a stabilizer code. A finite set of logical operators \(\mathcal{L} \subseteq \mathrm{LogicalOp}(C)\) satisfies the commutativity condition if every pair of operators in \(\mathcal{L}\) is Pauli-compatible:
The surrounding computation context for a gauging measurement, consisting of:
the position of the gauging measurement (beginning, middle, end, or isolated),
the total number of logical operations,
a Boolean indicating whether surrounding operations provide additional error correction,
the effective number of EC rounds provided by surrounding operations.
The syndrome of a spacetime fault \(f\) with respect to a detector collection \(\mathrm{DC}\) and base outcomes is the set of detectors violated by \(f\):
A CSS code (Calderbank–Shor–Steane code) is a structure consisting of:
A positive integer \(n\) (the number of physical qubits), with \(n {\gt} 0\).
A number \(n_X\) of \(X\)-type check generators and \(n_Z\) of \(Z\)-type check generators.
For each \(i \in \mathrm{Fin}(n_X)\), a nonempty support set \(\mathrm{xCheckSupport}(i) \subseteq \mathrm{Fin}(n)\) indicating the qubits on which \(X\) acts.
For each \(j \in \mathrm{Fin}(n_Z)\), a nonempty support set \(\mathrm{zCheckSupport}(j) \subseteq \mathrm{Fin}(n)\) indicating the qubits on which \(Z\) acts.
The CSS orthogonality condition: for all \(i \in \mathrm{Fin}(n_X)\) and \(j \in \mathrm{Fin}(n_Z)\),
\[ |\mathrm{xCheckSupport}(i) \cap \mathrm{zCheckSupport}(j)| \equiv 0 \pmod{2}. \]
The initialization hypergraph of a CSS code \(C\) is the hypergraph with vertex set \(\mathrm{Fin}(n)\) (physical qubits) and hyperedge set \(\mathrm{Fin}(n_Z)\) (one hyperedge per \(Z\)-type check), where the incidence function maps each \(Z\)-check index \(j\) to the support \(\mathrm{zCheckSupport}(j)\).
The CSS initialization vertex type for a CSS code \(C\) is the disjoint union of physical qubits and dummy vertices:
where:
\(\mathrm{qubit}(i)\) for \(i \in \mathrm{Fin}(n)\) represents physical qubit \(i\),
\(\mathrm{dummy}(j)\) for \(j \in \mathrm{Fin}(n_X)\) represents the dummy vertex for the \(j\)-th \(X\)-check.
The pairwise \(XX\) operator support for qubit index \(i \in \mathrm{Fin}(n)\) is the function \(\mathrm{SteaneVertex}(n) \to \mathbb {Z}/2\mathbb {Z}\) given by:
This represents the \(X_i^{\mathrm{data}} \otimes X_i^{\mathrm{ancilla}}\) operator acting on matching qubit pairs.
The readout support function assigns \(1 \in \mathbb {Z}/2\mathbb {Z}\) to ancilla vertices and \(0\) to data vertices:
This represents the fact that the readout step measures \(Z\) on ancilla qubits only.
The Steane gauging steps are the three phases of Steane-style measurement via gauging:
State preparation: Initialize ancilla code block via CSS gauging.
Entangling measurement: Pairwise \(XX\) gauging between data and ancilla blocks.
Readout: \(Z\) measurement on ancilla qubits (ungauging).
The Steane vertex type for \(n\) qubits is the disjoint union of a data block and an ancilla block:
where \(\mathrm{data}(i)\) represents qubit \(i\) in the data block and \(\mathrm{ancilla}(i)\) represents qubit \(i\) in the ancilla block.
The cycle degree of an edge \(e\) with respect to a set of cycles \(C\) (each given as a finset of edges) is the number of cycles in the generating set that contain \(e\):
A cycle-sparsification of a graph \(G\) with cycles \(C\) is a structure consisting of:
a base graph \(G\) with its chosen generating set of cycles (of type \(\mathrm{GraphWithCycles}\; V\; E\; C\)),
a number of additional layers \(R\) (total layers \(= R + 1\)),
a cycle-degree bound \(c \in \mathbb {N}\),
a cycle-layer assignment \(\mathrm{cycleAssignment} : C \to \mathrm{Fin}(R+1)\),
for each cycle, a list of vertices representing the cycle in order,
the cycle-degree constraint: for every edge \(e \in E\), \(\mathrm{cycleDegree}(\mathrm{cycles}, e) \leq c\).
The embedding \(V \hookrightarrow \mathrm{SparsifiedVertex}(S)\) that sends each vertex \(v\) to its layer-\(0\) copy \(\mathrm{ofOriginal}(v)\). This is injective: if \(\mathrm{ofOriginal}(v) = \mathrm{ofOriginal}(w)\), then \(v = w\) by the second component of the equality.
An inter-layer edge connects a vertex to its copy in an adjacent layer. Formally, \((\ell v_1, \ell v_2)\) is an inter-layer edge if:
An intra-layer edge in the sparsified graph connects two vertices in the same layer that are adjacent in the original graph. Formally, \((\ell v_1, \ell v_2)\) is an intra-layer edge if:
A triangulation edge \((\ell v_1, \ell v_2)\) is an edge added to cellulate a cycle. Formally:
The adjacency relation in the sparsified graph is defined as: two layered vertices \(\ell v_1, \ell v_2\) are adjacent if \(\ell v_1 \neq \ell v_2\) and at least one of the following holds:
\((\ell v_1, \ell v_2)\) is an intra-layer edge,
\((\ell v_1, \ell v_2)\) is an inter-layer edge,
\((\ell v_1, \ell v_2)\) is a triangulation edge.
The set of triangulation edge pairs for a cycle \(c\) in layer \(\ell \) is:
If the cycle \(c\) is assigned to layer \(\ell \): the image of \(\mathrm{triangulationEdgePairs}(\mathrm{cycleVertices}(c))\) under the embedding into layer \(\ell \).
Otherwise: the empty list.
A logical operator of the deformed code extends a deformed Pauli operator with:
Nontrivial: \(L\) is not the identity.
Cocycle: \(S_X^E\) is a \(1\)-cocycle (commutes with all flux operators \(B_p\)).
Cleaned vertex nontrivial: for any \(S\) with \(S_X^E = \partial S\), the cleaned vertex restriction \(\bar{L}|_V\) is nontrivial. This follows from \(L'\) being a logical (not a stabilizer): if \(\bar{L}|_V\) were trivial, then \(\bar{L}\) would be pure edge \(Z\)-support, hence a stabilizer, contradicting that \(L'\) is a logical.
Cleaned is original logical: for any \(S\) with \(S_X^E = \partial S\), the cleaned vertex restriction is a logical of the original code.
A Pauli operator on the deformed system (vertices \(+\) edges) is a structure
where \(S_X^V, S_Z^V \subseteq V\) are the \(X\)- and \(Z\)-supports on vertex qubits, \(S_X^E, S_Z^E \subseteq E\) are the \(X\)- and \(Z\)-supports on edge qubits, and \(\phi \in \mathbb {Z}/4\mathbb {Z}\) is a phase.
A detector for types \(V\), \(E\), \(M\) (with decidable equality) is a structure consisting of:
\(\mathrm{events} : \mathrm{Finset}(\mathrm{DetectorEvent}(V, E, M))\) — the set of events in the detector,
\(\mathrm{expectedParity} \in \mathbb {Z}/2\mathbb {Z}\) — the expected parity in fault-free execution (\(0\) for \(+1\), \(1\) for \(-1\)).
Given an outcome assignment and a time fault \(f\), the faulted outcomes are:
The combination of two detectors \(D_1\) and \(D_2\) is defined by taking the symmetric difference of their event sets and the XOR (sum mod \(2\)) of their expected parities:
This corresponds to multiplying the parity checks.
A structure capturing the deterministic outcome property of measuring \(X\) on \(\ket {+}\). It consists of:
A state \(\mathrm{is\_ plus\_ state}\) of type \(\mathrm{InitialState}\),
A proof that \(\mathrm{is\_ plus\_ state} = \mathrm{plus}\),
An outcome of type \(\mathrm{GraphMeasurementOutcome}\),
A proof that \(\mathrm{outcome} = \mathrm{plus}\).
The set of dummy vertices in the subdivided graph is \(\mathrm{dummyVerticesOfSubdivision}'(G) := \{ \mathrm{inr}(e) \mid e \in E(G)\} \subseteq \texttt{SubdividedVertex'}(V, \mathrm{Sym2}(V))\), defined as the image of \(G.\mathrm{edgeFinset}\) under \(\mathrm{inr}\).
The edge boundary \(\partial S\) of a vertex set \(S\) in a simple graph \(G = (V, E)\) is the set of edges with exactly one endpoint in \(S\):
Formally, \(\partial S\) is computed by filtering the edge set of \(G\) to those edges that cross the boundary of \(S\).
Given a vertex set \(S \subseteq V\) and an edge \(e = \{ u, v\} \in \mathrm{Sym}_2(V)\), we say that \(e\) crosses the boundary of \(S\) if exactly one of its endpoints lies in \(S\), i.e., \((u \in S \land v \notin S) \lor (u \notin S \land v \in S)\).
Effective code parameters under a given round configuration, consisting of:
the original code distance \(d\),
the effective distance (which may be reduced with fewer rounds),
a Boolean indicating whether a fault-tolerance threshold exists,
the round configuration used.
Given a simple graph \(G = (V, E)\) with decidable adjacency, the entangling circuit \(\prod _v \prod _{e \ni v} \mathrm{CX}_{v \to e}\) is a structure consisting of:
the underlying graph \(G\),
a gate specification: for each vertex \(v \in V\) and each edge \(e\) in the incidence set of \(v\), a \(\mathrm{CX}\) gate with control \(\mathrm{inl}(v) \in V \oplus \mathrm{Sym}_2(V)\) and target \(\mathrm{inr}(e) \in V \oplus \mathrm{Sym}_2(V)\).
A structure capturing the configuration of error correction rounds around the gauging measurement. It consists of:
\(d\): the code distance (a positive natural number),
roundsBefore: the number of EC rounds before gauging,
roundsAfter: the number of EC rounds after gauging.
An abstract error process in the fault-tolerant procedure, consisting of:
weight: the total number of faults,
involvesGauging: whether the process involves the gauging measurement,
involvesInitialBoundary: whether the process involves the initial boundary,
involvesFinalBoundary: whether the process involves the final boundary.
An error process spans to the boundary if it involves the gauging measurement and at least one boundary (initial or final):
The exactness condition for a graph with cycles \(G\) consists of:
Cycles valid: every cycle \(c \in C\) is a valid cycle edge set (every vertex has even degree in \(c\)).
Cuts generate: every element of \(\ker (\delta _2)\) lies in \(\operatorname {im}(\delta )\), i.e., every \(1\)-cocycle is a coboundary.
The expander existence specification is the proposition that for all \(W \ge 2\), there exists a simple graph \(G\) on \(\mathrm{Fin}(W)\) with decidable adjacency such that:
\(G\) has bounded degree: \(\exists \, d \in \mathbb {N}\) such that \(\forall \, v,\; \deg _G(v) \le d\);
\(G\) is a strong expander: \(\mathrm{IsStrongExpander}(G)\) holds (i.e., the Cheeger constant \(h(G) \ge 1\)).
This is a cited result from the expander graph literature (probabilistic method on random regular graphs, or explicit constructions such as Ramanujan graphs or Margulis–Gabber–Galil graphs).
A fault tolerance configuration is a structure bundling all preconditions needed to establish \(d_{ST} = d\). It consists of:
A code distance \(d \in \mathbb {N}\) with \(d {\gt} 0\).
An initial time \(t_i\) and final time \(t_o\) with \(t_i {\lt} t_o\) (the deformation interval is nonempty).
The condition \((t_o - t_i) \geq d\) (sufficient measurement rounds).
A Cheeger constant \(h(G) \in \mathbb {R}\) with \(h(G) \geq 1\) (strong expansion).
A spacetime fault \(F\) has a spacetime decomposition with respect to a detector collection \(DC\), base outcomes, logical effect predicate, and configuration \(\mathrm{cfg}\) if there exist spacetime faults \(F_{\mathrm{space}}\) and \(F_{\mathrm{time}}\) such that:
\(F\) is equivalent modulo stabilizers to \(F_{\mathrm{space}} \cdot F_{\mathrm{time}}\),
\(F_{\mathrm{space}}\) is a pure space fault at a single time step \(t_i\), and
\(F_{\mathrm{time}}\) is a pure time fault.
The lower bound case enumeration for a spacetime fault \(F\) consists of two cases:
Case 1 (timeNontrivial): There exists a pure time fault \(F_{\mathrm{time}}\) with \(\mathrm{weight}(F_{\mathrm{time}}) \geq \mathrm{numRounds}\) and \(\mathrm{weight}(F) \geq \mathrm{weight}(F_{\mathrm{time}})\).
Case 2 (spaceLogical): There exists a pure space fault \(F_{\mathrm{space}}\) at time \(t_i\) with \(\mathrm{weight}(F_{\mathrm{space}}) \geq d\) and \(\mathrm{weight}(F) \geq \mathrm{weight}(F_{\mathrm{space}})\).
An original code logical operator consists of:
A time step \(t\) of application (should be \({\lt} t_i\) or \({\gt} t_o\)).
A map \(\mathrm{vertexPaulis} : V \to \mathrm{PauliType}\) assigning Pauli operators to each vertex qubit.
A weight \(w \in \mathbb {N}\).
The conversion of an original logical operator \(L\) to a spacetime fault assigns:
Space errors: For a vertex qubit \(v\) at time \(t\), the error is \(L.\mathrm{vertexPaulis}(v)\) if \(t = L.\mathrm{time}\), and \(I\) otherwise. For edge qubits, the error is always \(I\).
Time errors: All time errors are \(\mathrm{false}\) (no measurement errors).
The type \(\mathrm{Finset}(\alpha )\) forms an additive commutative group with symmetric difference as addition, \(\emptyset \) as the zero element, and identity as negation (since every element is its own inverse in \(\mathbb {Z}/2\mathbb {Z}\)). The group axioms hold because:
Associativity: \((S \mathbin {\triangle } T) \mathbin {\triangle } U = S \mathbin {\triangle } (T \mathbin {\triangle } U)\).
Identity: \(\emptyset \mathbin {\triangle } S = S\) and \(S \mathbin {\triangle } \emptyset = S\).
Inverse: \(S \mathbin {\triangle } S = \emptyset \).
Commutativity: \(S \mathbin {\triangle } T = T \mathbin {\triangle } S\).
The type \(\mathrm{Finset}(\alpha )\) forms a module (vector space) over \(\mathbb {Z}/2\mathbb {Z}\), where addition is symmetric difference and scalar multiplication is as defined above. The module axioms are verified by case analysis on elements of \(\mathbb {Z}/2\mathbb {Z}\) (which are either \(0\) or \(1\)).
A summary structure capturing the complete trade-off for flux measurement frequency. It records whether \(B_p\) can be measured less often, whether it can be inferred entirely, whether fault distance scaling is preserved, whether detector cells are large, whether a threshold is expected, and whether the strategy is useful for small instances.
The Freedman–Hastings decongestion specification asserts that there exists a constant \(C {\gt} 0\) such that for all \(W \ge 2\) and \(d \ge 1\) (representing any constant-degree-\(d\) graph on \(W\) vertices), there exist natural numbers \(R\) and \(\mathrm{cycleWeightBound}\) satisfying:
\(R \le C \cdot (\log _2 W)^2 + C\) (i.e., \(R = O(\log ^2 W)\));
\(\mathrm{cycleWeightBound} \le 3\) (all generating cycles have weight at most \(3\), achieved by triangulation).
This is a cited result from Freedman–Hastings requiring topological methods.
The complete gauging circuit for a simple graph \(G = (V,E)\) is a structure consisting of:
the underlying graph \(G\),
a first entangling circuit \(\mathcal{E}_1\) of type \(\mathrm{EntanglingCircuit}(G)\),
an \(X\)-measurement specification on the vertices \(V\),
a second entangling circuit \(\mathcal{E}_2\) of type \(\mathrm{EntanglingCircuit}(G)\),
a \(Z\)-measurement specification on the edge set \(\mathrm{Sym}_2(V)\),
subject to the following conditions:
\(\mathcal{E}_1 = \mathcal{E}_2\) (the two entangling circuits are identical),
the \(X\)-measurement covers all vertices: \(\mathrm{qubits} = V\),
the \(X\)-measurement keeps all qubits: \(\mathrm{keep}(v) = \mathrm{true}\) for all \(v\),
the \(Z\)-measurement covers all edges: \(\mathrm{qubits} = E_G\),
the \(Z\)-measurement discards all qubits: \(\mathrm{keep}(e) = \mathrm{false}\) for all \(e\).
A circuit qubit \(q\) is an edge qubit if it has the form \(\mathrm{edge}(e)\) for some \(e\). Formally,
A gauging graph convention specifies how a connected graph \(G\) relates to the logical operator \(L\) being measured. It consists of:
A vertex type \(V\) with a \(\mathrm{Fintype}\) instance and decidable equality,
A simple graph structure on \(V\) with decidable adjacency,
A proof that the graph is connected,
A finite set \(\mathrm{logicalSupport} \subseteq V\) representing the support of \(L\),
A proof that \(\mathrm{logicalSupport}\) is nonempty (i.e., \(L\) is nontrivial).
This captures the conventions: (1) vertices of \(G\) are identified with qubits in \(\operatorname {supp}(L)\), (2) each edge \(e \in E_G\) gets an auxiliary gauge qubit in \(\ket {0}\), (3) dummy vertices (\(V_G \setminus \operatorname {supp}(L)\)) get auxiliary qubits in \(\ket {+}\), and (4) dummy vertices have no effect on the measurement outcome.
The effective logical support when gauging \(L\) with dummy vertices is the entire vertex set \(V_G\). This captures the convention that we gauge the operator \(L \cdot \prod _{v \in \mathrm{dummy}} X_v\):
The Three Desiderata for Gauging Graph \(G\).
A graph \(G\) (equipped with a cycle structure) satisfies the gauging graph desiderata with respect to a collection of \(Z\)-type supports \(\mathcal{Z}\), path bound \(k\), and cycle bound \(w\) if the following three conditions hold simultaneously:
Short paths: \(\mathrm{ShortPaths}(G.\mathrm{graph}, \mathcal{Z}, k)\) — there exist paths of length at most \(k\) between any two vertices in the same \(Z\)-type support set.
Sufficient expansion: \(\mathrm{SufficientExpansion}(G.\mathrm{graph})\) — the Cheeger constant satisfies \(h(G) \geq 1\).
Low-weight cycles: \(\mathrm{LowWeightCycles}(G, w)\) — all generating cycles have weight at most \(w\).
The informal justifications from the remark are:
Short paths \(\Longrightarrow \) deformed checks have bounded weight.
\(h(G) \geq 1\) \(\Longrightarrow \) deformed code distance \(\geq \) original code distance.
Low-weight cycles \(\Longrightarrow \) \(B_p\) flux operators have bounded weight.
All three together \(\Longrightarrow \) constant qubit overhead.
Given \(z \in (\mathbb {Z}/2\mathbb {Z})^E\) in the image of \(\delta \) (i.e., there exists \(c\) with \(\delta c = z\)), the byproduct cochain \(c' = \mathrm{byproductCochain}(G, z)\) is a chosen \(0\)-cochain satisfying \(\delta c' = z\), obtained by the axiom of choice.
For a \(0\)-cochain \(c \in (\mathbb {Z}/2\mathbb {Z})^V\) and outcomes \(\varepsilon \), define
This represents \(\prod _{v : c_v = 1} \varepsilon _v^{c_v}\) in additive notation.
The Gauss law measurement outcomes assign to each vertex \(v \in V\) a value in \(\mathbb {Z}/2\mathbb {Z}\), where \(0\) represents the outcome \(+1\) and \(1\) represents the outcome \(-1\). Formally, \(\mathrm{GaussLawOutcomes}(V) = V \to \mathbb {Z}/2\mathbb {Z}\).
The measured outcome \(\sigma \) is the product of all Gauss law outcomes, represented in \(\mathbb {Z}/2\mathbb {Z}\) as
Here \(\sigma = 0\) means \(+1\) (even number of \(-1\) outcomes) and \(\sigma = 1\) means \(-1\) (odd number).
A configuration of key time points in the gauging procedure, consisting of:
\(t_{\mathrm{start}}\): the start of the procedure,
\(t_{\mathrm{initial}}\): the time of the initial gauging code deformation,
\(t_{\mathrm{final}}\): the time of the final ungauging code deformation,
subject to the ordering constraint \(t_{\mathrm{start}} \le t_{\mathrm{initial}} \le t_{\mathrm{final}}\).
The four directions of generalization for the gauging procedure form an inductive type with constructors:
finiteGroupReps: Any finite group representation with tensor product factorization.
nonPauliOperators: Non-Pauli operators (which can produce magic states).
quditSystems: Qudit systems with \(d {\gt} 2\) levels per site.
nonabelianGroups: Nonabelian groups (with a local-vs-global caveat).
A graph with cycles is a structure \((G, V, E, C)\) where:
\(V\), \(E\), \(C\) are finite types with decidable equality,
\(G\) is a simple graph on \(V\) with decidable adjacency,
\(\mathrm{edgeEndpoints} : E \to V \times V\) assigns to each edge its two endpoints,
for each edge \(e\), the endpoints are adjacent: \(G.\mathrm{Adj}(\mathrm{edgeEndpoints}(e)_1, \mathrm{edgeEndpoints}(e)_2)\),
edges are symmetric: adjacency holds in both directions,
\(\mathrm{cycles} : C \to \mathrm{Finset}(E)\) assigns to each cycle index its set of edges.
The standard basis vectors are defined as:
\(\mathrm{basisV}(v) = \mathbf{e}_v \in \mathbb {Z}_2^V\): the vector with \(1\) at position \(v\) and \(0\) elsewhere,
\(\mathrm{basisE}(e) = \mathbf{e}_e \in \mathbb {Z}_2^E\): the vector with \(1\) at position \(e\) and \(0\) elsewhere,
\(\mathrm{basisC}(c) = \mathbf{e}_c \in \mathbb {Z}_2^C\): the vector with \(1\) at position \(c\) and \(0\) elsewhere.
Formally, \(\mathrm{basisV}(v) = \pi _{\mathrm{single}}(v, 1)\) (and similarly for edges and cycles).
The second boundary map \(\partial _2 : \mathbb {Z}_2^C \to \mathbb {Z}_2^E\) is the \(\mathbb {Z}_2\)-linear map defined by:
where \(\partial _2 c\) is the characteristic vector of edges in cycle \(c\).
The boundary map \(\partial : \mathbb {Z}_2^E \to \mathbb {Z}_2^V\) is the \(\mathbb {Z}_2\)-linear map defined by:
where \(\partial e = \mathbf{e}_v + \mathbf{e}_{v'}\) for edge \(e = \{ v, v'\} \).
The second coboundary map \(\delta _2 : \mathbb {Z}_2^E \to \mathbb {Z}_2^C\) is the \(\mathbb {Z}_2\)-linear map defined by:
where \(\delta _2 e\) is the characteristic vector of cycles containing \(e\).
A deformable Pauli operator on a graph \(G = (V, E, C)\) is a structure consisting of:
\(S_X^V \subseteq V\): the \(X\)-type support on vertices,
\(S_Z^V \subseteq V\): the \(Z\)-type support on vertices,
\(S_X^E \subseteq E\): the \(X\)-type support on edges,
\(S_Z^E \subseteq E\): the \(Z\)-type support on edges,
\(\sigma \in \mathbb {Z}_4\): the global phase (\(0 = +1\), \(1 = +i\), \(2 = -1\), \(3 = -i\)),
The deformability condition: \(|S_Z^V| \equiv 0 \pmod{2}\).
The deformability condition is equivalent to \(P\) commuting with the logical operator \(L = \prod _{v \in V} X_v\).
The symplectic form between the deformed operator \(\tilde{P}\) and the Gauss law operator \(A_v\) is defined as the count of anticommuting pairs:
where \(\mathrm{inc}(v)\) denotes the set of edges incident to \(v\).
The deformed check of a stabilizer check \(s\) along an edge-path \(\gamma \) is defined as
which is an alias for the DeformedOperator construction applied to a stabilizer check.
A valid deformed code setup for a graph \(G\) and code parameters \(\mathrm{code}\) consists of:
The original stabilizer checks: \(\mathrm{checks} : \mathrm{Fin}(n-k) \to \mathrm{StabilizerCheck}(G)\).
Valid deforming paths for each check: \(\mathrm{paths}(j)\) is an edge path for each \(j\).
Original checks have no Z-support on edges: \(\forall j,\; \mathrm{checks}(j).\mathrm{zSupportOnE} = \emptyset \).
Paths satisfy the boundary condition: \(\forall j,\; \partial (\mathrm{paths}(j)) = \mathcal{S}_Z(\mathrm{checks}(j)) \cap V_G\).
All cycles in \(C\) are valid (geometric constraint).
\(|V| \geq 1\).
Euler’s formula holds: \(|C| = |E| - |V| + 1\).
Graph is connected: \(|V| \leq |E|\).
The original code has at least one logical qubit: \(k \geq 1\).
Given a deformable Pauli operator \(P\) on graph \(G\) and an edge-path \(\gamma \subseteq E\), the deformed operator \(\tilde{P} = P \cdot \prod _{e \in \gamma } Z_e\) is defined as the deformable Pauli operator with:
\(S_X^V(\tilde{P}) = S_X^V(P)\) (X-support on vertices unchanged),
\(S_Z^V(\tilde{P}) = S_Z^V(P)\) (Z-support on vertices unchanged),
\(S_X^E(\tilde{P}) = S_X^E(P)\) (X-support on edges unchanged),
\(S_Z^E(\tilde{P}) = S_Z^E(P) \oplus \gamma \) (Z-support on edges is the symmetric difference with \(\gamma \)),
\(\sigma (\tilde{P}) = \sigma (P)\) (phase unchanged).
The deformability condition \(|S_Z^V(\tilde{P})| \equiv 0 \pmod{2}\) is inherited from \(P\).
The boundary of an edge-path \(\gamma \) is defined as \(\partial \gamma = \partial _1(\boldsymbol {\gamma })\), where \(\partial _1\) is the boundary map and \(\boldsymbol {\gamma }\) is the binary vector representation of \(\gamma \). This gives a vector in \(\mathbb {Z}_2^V\) whose value at each vertex \(v\) counts (modulo 2) the number of edges in \(\gamma \) incident to \(v\).
Given a graph \(G\) and an edge-path \(\gamma \subseteq E\), the edge-path vector is the binary vector \(\boldsymbol {\gamma } \in \mathbb {Z}_2^E\) defined by
The symplectic form between two flux operators. For \(Z\)-type operators:
Formally:
The edge support of flux operator \(B_p\): \(1\) at edges in cycle \(p\), \(0\) elsewhere. This represents \(B_p = \prod _{e \in p} Z_e\). For a graph \(G\) and cycle \(p\):
i.e., the image of cycle \(p\) under the second boundary map.
The vertex support of flux operator \(B_p\): zero everywhere. Since \(B_p\) is a \(Z\)-type operator, it has no \(X\) component on vertices. For any graph \(G\) and cycle \(p\):
A gauging measurement for a graph \(G\) consists of a measurement outcome function \(\varepsilon : V \to \mathbb {Z}/2\mathbb {Z}\), recording the outcome \(\varepsilon _v \in \{ +1, -1\} \) for each Gauss law operator \(A_v\) (where \(0\) encodes \(+1\) and \(1\) encodes \(-1\)).
The symplectic form between two Gauss law operators \(A_v\) and \(A_w\) is defined as
For \(X\)-type operators, this equals \(|\emptyset | + |\emptyset | = 0\).
The edge support of Gauss law operator \(A_v\) is \(1\) at edges incident to \(v\) and \(0\) elsewhere. Formally,
the coboundary of vertex \(v\). This represents \(\prod _{e \ni v} X_e\) in the product \(A_v = X_v \prod _{e \ni v} X_e\).
The vertex support of Gauss law operator \(A_v\) is a binary vector with \(1\) at position \(v\) and \(0\) elsewhere. Formally,
the basis vector at \(v\). This represents \(X_v\) in the product \(A_v = X_v \prod _{e \ni v} X_e\).
A graph \(G\) has valid paths if for any even-cardinality subset \(S \subseteq V\), there exists an edge-path \(\gamma \) such that \(\gamma \) is a valid deforming path for \(S\):
This follows from graph connectivity: for two vertices, take any path between them; for larger even-cardinality sets, pair up vertices and connect each pair.
The number of independent cycles in a connected graph according to Euler’s formula. For a connected graph with vertex set \(V\) and edge set \(E\):
We represent this as an integer to handle the subtraction correctly:
The initial stabilizer condition: each edge qubit in state \(\lvert 0\rangle \) is stabilized by \(Z_e\). Represented as \(Z\lvert 0\rangle = +\lvert 0\rangle \), so the measurement outcome is \(0\) (for \(+1\)) in \(\mathbb {Z}/2\mathbb {Z}\):
An edge-path \(\gamma \) is a minimum-weight valid deforming path for \(S_Z^V\) if:
\(\gamma \) is a valid deforming path: \(\partial \gamma = \mathbf{s}_{S_Z^V}\), and
\(\gamma \) has minimum weight among all valid paths: for all valid paths \(\gamma '\), \(w(\gamma ) \leq w(\gamma ')\).
A cycle \(p\) is valid if every vertex has \(0\) or \(2\) incident edges from the cycle:
An edge-path \(\gamma \) is a valid deforming path for a \(Z\)-support \(S \subseteq V\) if the boundary of \(\gamma \) equals the \(Z\)-support vector:
where \(\mathbf{s}_S\) is the binary vector representation of \(S\). Equivalently, \(\partial \gamma = S\) as subsets of \(V\) (via their characteristic functions in \(\mathbb {Z}_2^V\)).
An original code state encodes the input assumption that the original code state \(|\psi \rangle \) is in the \(+1\) eigenspace of all stabilizer generators. This is the defining property of stabilizer codes: the code space is the simultaneous \(+1\) eigenspace of all stabilizer generators.
A stabilizer check from the original code is a DeformablePauliOperator on the graph \(G\). It carries the additional semantic meaning that it represents a generator of the original code’s stabilizer group. The key inherited property is that \(|\mathcal{S}_Z \cap V_G|\) is even (the check commutes with \(L\)).
We define the binary vector spaces:
\(\mathbb {Z}_2^V := V \to \mathbb {Z}_2\) (binary vectors over vertices),
\(\mathbb {Z}_2^E := E \to \mathbb {Z}_2\) (binary vectors over edges),
\(\mathbb {Z}_2^C := C \to \mathbb {Z}_2\) (binary vectors over cycles).
Each carries a natural \(\mathbb {Z}_2\)-module structure.
Given a graph \(G = (V, E, C)\) and a Pauli operator with \(Z\)-type support \(S_Z \subseteq V\), the \(Z\)-type support restricted to the graph vertices is \(S_Z \cap V_G\). Since \(S_Z\) is already given as a subset of \(V\), this is simply \(S_Z\) itself.
The concrete gauging graph for \(\bar{X}_\alpha \) in the Gross code, with parameters:
\(|V| = 12\) vertices,
\(18\) initial edges from \(Z\)-check connectivity,
\(4\) expansion edges,
\(4\) redundant cycles,
maximum new degree \(6\),
maximum affected degree \(7\).
The structure recording the parameters of the gauging graph for \(\bar{X}_\alpha \) in the Gross code. It consists of:
numVertices : the number of vertices (monomials in \(f\)),
numInitialEdges : the number of initial edges from \(Z\)-check connectivity,
numExpansionEdges : the number of additional expansion edges,
numRedundantCycles : the number of redundant cycles from the BB code’s \(Z\)-check structure,
maxNewDegree : the maximum Tanner graph degree for new elements,
maxAffectedDegree : the maximum Tanner graph degree for affected existing elements.
The \(12\) monomials appearing in \(f\), represented as elements of \(\mathbb {Z}/12\mathbb {Z} \times \mathbb {Z}/6\mathbb {Z}\):
These correspond to the monomials of \(f = 1 + x + x^2 + x^3 + x^6 + x^7 + x^8 + x^9 + (x + x^5 + x^7 + x^{11})y^3\).
The number of \(Z\) checks adjacent to \(\bar{X}_\alpha \):
These arise from the connectivity condition: vertices \(\gamma , \delta \in f\) are connected if \(\gamma = B_i^T B_j \delta \) for some \(i, j \in \{ 1,2,3\} \).
A half-integer time representation, defined as an inductive type with two constructors:
\(\mathtt{integer}(n)\) represents the integer time step \(n \in \mathbb {Z}\),
\(\mathtt{halfInteger}(n)\) represents the half-integer time step \(n + \tfrac {1}{2}\) for \(n \in \mathbb {Z}\).
The predecessor function on half-integer times, retreating by \(\tfrac {1}{2}\):
That is, \(n \mapsto n - \tfrac {1}{2}\) and \(n + \tfrac {1}{2} \mapsto n\).
The successor function on half-integer times, advancing by \(\tfrac {1}{2}\):
That is, \(n \mapsto n + \tfrac {1}{2}\) and \(n + \tfrac {1}{2} \mapsto n + 1\).
A class asserting that a graph \(G\) satisfies the Freedman–Hastings decongestion bound. This states:
There exists a layer bound \(\mathrm{lb}\) with \(\mathrm{lb}.\mathrm{numLayers} \leq (\log _2 |V| + 1)^2 \cdot (\mathrm{maxDegree} + 1)\).
The cycle-degree bound \(c {\gt} 0\).
The graph has bounded degree: for all \(v \in V\), \(|G.\mathrm{incidentEdges}(v)| \leq \mathrm{maxDegree}\).
The predicate \(\mathrm{hasLogicalFault}(\mathit{DC}, \mathit{logicalEffect})\) holds if there exists at least one spacetime fault \(F\) such that \(F\) is a spacetime logical fault:
A hypergraph on vertices \(V\) with hyperedges indexed by \(H\) is a structure consisting of an incidence function
mapping each hyperedge to the set of its incident vertices. This generalizes a graph, where each edge connects exactly two vertices, to the case where each hyperedge may connect an arbitrary subset of vertices.
The hyperedge support of \(A_v\) is a binary vector on \(H\) with \(1\) at each hyperedge containing \(v\):
The vertex support of the generalized Gauss law operator \(A_v\) is a binary vector on \(V\) with \(1\) at position \(v\) and \(0\) elsewhere, representing the \(X_v\) factor:
The generalized coboundary map \(\delta : \mathbb {Z}_2^V \to \mathbb {Z}_2^H\) for the hypergraph is the \(\mathbb {Z}_2\)-linear map defined by
This is the transpose of the incidence matrix.
The incidence entry of a hypergraph \(\mathrm{HG}\) at hyperedge \(h\) and vertex \(v\) is defined as
where the values are taken in \(\mathbb {Z}/2\mathbb {Z}\).
Given a \(\mathrm{GraphWithCycles}\) \(G\) on vertices \(V\), edges \(E\), and cycles \(C\), we construct a hypergraph on vertices \(V\) with hyperedges indexed by \(E\) by setting
for each edge \(e \in E\). This embeds ordinary graphs into the hypergraph framework.
The operator kernel of a hypergraph \(\mathrm{HG}\) is the kernel of the parity-check map:
This is a \(\mathbb {Z}_2\)-submodule of \(\mathbb {Z}_2^V\) that characterizes the abelian group of commuting \(X\)-type operators.
The parity-check map \(H_Z : \mathbb {Z}_2^V \to \mathbb {Z}_2^H\) of a hypergraph \(\mathrm{HG}\) is the \(\mathbb {Z}_2\)-linear map defined by
for each hyperedge \(h \in H\). The kernel of this map characterizes the abelian group of commuting \(X\)-type operators.
The symplectic inner product on vertex qubits for two Pauli operators with \(X\)-type supports \(xS_1, xS_2\) and \(Z\)-type supports \(zS_1, zS_2\) is
The complete trade-off for the inferred flux measurement strategy:
\(B_p\) can be measured less often: true,
\(B_p\) can be inferred: true,
Fault distance preserved: true,
Large detector cells: true,
Threshold expected: false,
Useful for small instances: true.
An initialization event for vertex type \(V\) and edge type \(E\) consists of:
A qubit location \(\mathrm{qubit} : \mathrm{QubitLoc}(V, E)\),
A time step \(\mathrm{time} : \mathrm{TimeStep}\),
An expected parity \(\mathrm{expectedParity} \in \mathbb {Z}/2\mathbb {Z}\).
An initialization fault over qubit types \(V, E\) is a structure consisting of:
qubit : \(\texttt{QubitLoc}(V, E)\) — the qubit being initialized,
time : \(\mathbb {N}\) — the time step of initialization,
effectiveError : \(\texttt{PauliType}\) — the effective Pauli error applied after perfect initialization.
A faulty initialization of state \(|\psi \rangle \) is modeled as perfect initialization followed by an immediate space-fault.
An initialization fault \(f\) is converted to an equivalent space-fault:
The space-fault occurs at the same qubit and time as the initialization.
A stabilizer code \(C\) is Low-Density Parity-Check (LDPC) if there exist constants bounding:
Maximum check weight \(w_{\max } \in \mathbb {N}\): for every generator \(s \in C.\texttt{stabilizers}.\texttt{generators}\), \(|\operatorname {supp}(s)| \leq w_{\max }\).
Maximum qubit degree \(\Delta _{\max } \in \mathbb {N}\): for every qubit \(i \in \operatorname {Fin}(C.n)\), the number of generators whose support contains \(i\) is at most \(\Delta _{\max }\).
A spacetime fault \(F\) is a logical fault if it has empty syndrome but is non-trivial:
A fault \(f\) is a spacetime logical fault with respect to a detector collection \(\mathrm{DC}\), base outcomes, and logical effect predicate \(\ell \) if:
\(f\) has empty syndrome: \(\mathrm{hasEmptySyndrome}(\mathrm{DC}, \mathrm{base}, f)\), and
\(f\) affects logical information: \(\mathrm{affectsLogicalInfo}(\ell , f)\).
A fault \(f\) is a spacetime stabilizer with respect to a detector collection \(\mathrm{DC}\), base outcomes, and logical effect predicate \(\ell \) if:
\(f\) has empty syndrome: \(\mathrm{hasEmptySyndrome}(\mathrm{DC}, \mathrm{base}, f)\), and
\(f\) preserves logical information: \(\mathrm{preservesLogicalInfo}(\ell , f)\).
A spacetime fault \(F\) is a spacetime stabilizer if it has empty syndrome and is trivial:
The adjacency relation for the combined graph used in generalized lattice surgery: two copies of a graph \(G\) on vertex set \(V\) with a set of bridge edges. Two vertices \(x, y \in \mathrm{Bool} \times V\) are adjacent if:
\(x_1 = y_1\) and \(G.\mathrm{Adj}\; x_2\; y_2\) (intra-copy adjacency), or
\(x_1 \neq y_1\) and \(((x_2, y_2) \in \mathrm{bridges}\) or \((y_2, x_2) \in \mathrm{bridges})\) (cross-copy bridge).
Grid graph adjacency on \(\mathrm{Fin}(\mathrm{rows}) \times \mathrm{Fin}(\mathrm{gap\_ width} + 2)\): two vertices \(x, y\) are adjacent if they are unit-distance neighbors, i.e.,
The adjacency relation for the ladder graph on \(n\) rungs. Two vertices \(x, y \in \mathrm{Bool} \times \mathrm{Fin}(n)\) are adjacent if:
\(x_1 \neq y_1\) and \(x_2 = y_2\) (horizontal rung), or
\(x_1 = y_1\) and \(x_2 + 1 = y_2\) (vertical rail going up), or
\(x_1 = y_1\) and \(y_2 + 1 = x_2\) (vertical rail going down).
The edge type for the ladder graph on \(n\) rungs is an inductive type with three constructors:
\(\mathrm{rung}(i)\) for \(i \in \mathrm{Fin}(n)\): horizontal rung connecting \((\mathrm{false}, i)\) to \((\mathrm{true}, i)\).
\(\mathrm{leftRail}(i)\) for \(i \in \mathrm{Fin}(n-1)\): left vertical rail connecting \((\mathrm{false}, i)\) to \((\mathrm{false}, i+1)\).
\(\mathrm{rightRail}(i)\) for \(i \in \mathrm{Fin}(n-1)\): right vertical rail connecting \((\mathrm{true}, i)\) to \((\mathrm{true}, i+1)\).
The endpoint function for each ladder edge:
\(\mathrm{rung}(i) \mapsto ((\mathrm{false}, i), (\mathrm{true}, i))\),
\(\mathrm{leftRail}(i) \mapsto ((\mathrm{false}, i), (\mathrm{false}, i+1))\),
\(\mathrm{rightRail}(i) \mapsto ((\mathrm{true}, i), (\mathrm{true}, i+1))\).
The ladder graph on \(n \geq 2\) rungs, built as a \(\mathrm{GraphWithCycles}\) instance to connect to the gauging framework (Gauss’s law operators from Def 2, flux operators from Def 3). The graph is the ladder graph, with edge endpoints and cycle structure as defined above.
Relaxed expansion: the condition \(|\partial S| \geq |S|\) is required only for subsets \(S \subseteq V\) satisfying \(0 {\lt} |S|\), \(2|S| \leq |V|\), and \((S \cap \mathrm{logicalSupport}) \neq \emptyset \). Formally,
A layer bound for a graph \(G\) with cycle-degree bound \(c\) consists of:
a number of layers \(R \in \mathbb {N}\),
a proof that a valid cycle-layer assignment exists with \(R\) layers and cycle-degree bound \(c\).
This represents the value \(R_G^c\) from the paper.
A layered vertex in a graph with \(R+1\) layers is a pair \((i, v)\) where \(i \in \{ 0, 1, \ldots , R\} \) is the layer index and \(v \in V\) is the original vertex. Layer \(0\) contains the original vertices; layers \(1, \ldots , R\) contain dummy copies.
Formally, given a type \(V\) and a natural number \(R\), a layered vertex is a structure consisting of:
a layer index \(\mathrm{layer} \in \mathrm{Fin}(R+1)\),
an original vertex \(\mathrm{vertex} \in V\).
Alogical effect predicate \(\ell \) is group-like if:
The identity preserves logical info: \(\neg \ell (1)\).
Closure under multiplication: for all \(f, g\), if \(\neg \ell (f)\) and \(\neg \ell (g)\), then \(\neg \ell (f \cdot g)\).
Closure under inverse: for all \(f\), if \(\neg \ell (f)\), then \(\neg \ell (f^{-1})\).
Given a detector collection \(\mathit{DC}\), base outcomes, a logical effect predicate, and a set of time steps \(\mathit{times}\), the set of logical fault weights is
A logical operator for a stabilizer code \(C\) is a structure consisting of:
an underlying Pauli operator \(\texttt{op} : \texttt{PauliOp}(C.n)\),
a proof that \(\texttt{op}\) commutes with all stabilizers,
a proof that \(\texttt{op}\) is not in the stabilizer group.
Desideratum 3: Low-Weight Cycles Property.
Let \(G\) be a graph with a chosen collection of cycles (a \(\mathrm{GraphWithCycles}\) structure over vertex type \(V\), edge type \(E\), and cycle type \(C\)), and let \(w \in \mathbb {N}\) be a constant weight bound.
The low-weight cycles property \(\mathrm{LowWeightCycles}(G, w)\) holds if every cycle \(c \in C\) has at most \(w\) edges:
This ensures that the \(B_p\) flux operators \(B_p = \prod _{e \in p} Z_e\) have bounded weight.
Measurement frequency options for stabilizer checks:
everyRound: Measured in every error correction round (standard approach),
sparse: Measured less often than every round,
inferred: Never measured directly; inferred from other data (initialization and readout).
A measurement schedule assigns a measurement frequency to each check type. It consists of:
gaussLawFreq: the frequency for Gauss law operators \(A_v\),
deformedFreq: the frequency for deformed checks \(\tilde{s}_j\),
fluxFreq: the frequency for flux operators \(B_p\).
Given a simple graph \(G = (V,E)\), the default gauging circuit \(\mathrm{mkGaugingCircuit}(G)\) is the gauging circuit where:
both entangling circuits use \(\mathrm{CX}\) gates with \(\mathrm{control} = \mathrm{inl}(v)\) and \(\mathrm{target} = \mathrm{inr}(e)\) for each vertex-edge incidence pair,
the \(X\)-measurement covers all vertices with \(\mathrm{keep} = \mathrm{true}\),
the \(Z\)-measurement covers all edges of \(G\) with \(\mathrm{keep} = \mathrm{false}\).
All required conditions (identical entangling circuits, full coverage, keep/discard properties) hold by reflexivity.
The original code distance structure captures:
A positive integer \(d {\gt} 0\) (the distance).
A predicate \(\mathrm{isLogical}\) identifying logical operators of the original code.
The distance property: for every nontrivial logical operator \(\mathrm{op}\),
\[ \mathrm{op}.\mathrm{weight} \; \ge \; d. \]
The set of original vertices in the subdivided graph is \(\mathrm{originalVertices}'(G) := \{ \mathrm{inl}(v) \mid v \in V\} \subseteq \texttt{SubdividedVertex'}(V, \mathrm{Sym2}(V))\), defined as the image of \(\mathrm{Finset.univ}\) under \(\mathrm{inl}\).
An OverheadBreakdown provides a step-by-step breakdown of overhead sources:
\(W\): weight of the logical operator,
\(\mathtt{matchingEdges} \leq W\): edges from the matching step,
\(\mathtt{expansionEdges} \leq 3W\): edges from the expansion step (constant-degree expanders),
\(\mathtt{layers}\): number of sparsification layers, satisfying \(\mathtt{layers} \leq C \cdot (\log _2 W)^2 + C\) for all \(C {\gt} 0\).
The qubit overhead structure for a given configuration, setting:
The overlap degree of a qubit \(i\) with respect to a set of logical operators \(\mathcal{L}\) is the number of operators in \(\mathcal{L}\) whose support contains \(i\):
A set of logical operators \(\mathcal{L}\) has pairwise disjoint support if every pair of distinct operators has disjoint support:
A set of logical operators is parallelizable if it satisfies both conditions for parallel gauging measurement. A ParallelizableLogicals structure for a stabilizer code \(C\) consists of:
A finite set \(\mathcal{L}\) of logical operators (the operators to measure in parallel),
An overlap bound \(k \in \mathbb {N}\) (must be constant to maintain the LDPC property),
A proof that \(\mathcal{L}\) satisfies the commutativity condition,
A proof that \(\mathcal{L}\) satisfies the bounded overlap condition with bound \(k\).
Given a set of logical operators \(\mathcal{L}\) with pairwise disjoint support, we construct a ParallelizableLogicals structure with overlap bound \(1\), using Lemma 292 for the commutativity condition and Lemma 293 for the bounded overlap condition.
Two Pauli operators \(P, Q \in \mathrm{PauliOp}(n)\) are Pauli-compatible if for every qubit \(i \in \mathrm{Fin}(n)\), their Pauli types at position \(i\) are compatible:
A Pauli operator on \(n\) qubits is a structure consisting of:
a function \(\texttt{paulis} : \operatorname {Fin}(n) \to \texttt{StabPauliType}\) assigning a single-qubit Pauli type to each qubit, and
a global phase \(\texttt{phase} \in \mathbb {Z}/4\mathbb {Z}\), where \(0 = +1\), \(1 = +i\), \(2 = -1\), \(3 = -i\).
For a Pauli operator \(P\) on \(n\) qubits and a set \(S \subseteq \mathrm{Fin}(n)\), the anticommuting positions with pure \(X\) is
Two Pauli operators \(P\) and \(Q\) on \(n\) qubits commute (ignoring global phase) if and only if the number of positions where their single-qubit Pauli types anticommute is even:
Given a finite set \(S \subseteq \operatorname {Fin}(n)\), the pure \(X\)-type operator \(\operatorname {pureX}(S)\) is the Pauli operator with
and phase \(0\).
The X-Z decomposition of a Pauli operator \(P\) on \(n\) qubits is a structure witnessing that:
The X-type support of \(P\) is exactly \(\{ i \mid (\mathrm{xzExponents}(P_i))_1 = 1\} \).
The Z-type support of \(P\) is exactly \(\{ i \mid (\mathrm{xzExponents}(P_i))_2 = 1\} \).
This captures the decomposition \(P = i^{\sigma } \prod _{v \in \mathcal{S}_X} X_v \prod _{v \in \mathcal{S}_Z} Z_v\).
Every single-qubit Pauli type can be written as \(X^a Z^b\) for \(a, b \in \{ 0,1\} \) (ignoring phase). The X-Z exponents map assigns:
The multiplication of Pauli types is defined by the standard Pauli multiplication table (ignoring global phases):
\(I\) is the identity: \(I \cdot p = p \cdot I = p\) for all \(p\).
Each Pauli is self-inverse: \(X^2 = Y^2 = Z^2 = I\).
Products of distinct non-identity Paulis: \(XY = Z\), \(YX = Z\), \(YZ = X\), \(ZY = X\), \(ZX = Y\), \(XZ = Y\).
A structure encoding the Perfect Boundary Assumption: the boundary measurements are error-free. It consists of:
initialQuality: the quality of initial boundary measurements (default: perfect),
finalQuality: the quality of final boundary measurements (default: perfect).
A summary structure capturing all aspects of the practical measurement rounds trade-off:
the theoretical requirement (number of rounds),
whether this is a worst-case guarantee,
whether practical implementations may use fewer rounds,
whether the choice affects distance,
whether the choice is context-dependent.
The set of pure space-faults forms a subgroup of \(\texttt{SpacetimeFault}(V, E, M)\):
This is verified to be closed under multiplication (since \(\text{false} \oplus \text{false} = \text{false}\)), contains the identity (which has all time errors \(\text{false}\)), and is closed under inversion (since the inverse preserves time errors).
The set of pure time-faults forms a subgroup of \(\texttt{SpacetimeFault}(V, E, M)\):
This is verified to be closed under multiplication (since \(I \cdot I = I\)), contains the identity (which has all space errors \(I\)), and is closed under inversion (since \(I^{-1} = I\)).
A QubitOverhead structure records the components of the qubit overhead in the gauging measurement construction. Each edge qubit is an auxiliary qubit initialized in \(|0\rangle \). It consists of:
\(W\): the weight of the logical operator (number of qubits in its support),
\(\mathtt{baseEdges}\): the number of edges in the base graph \(G_0\),
\(\mathtt{numLayers}\): the number of layers added for cycle sparsification,
\(\mathtt{cellulationEdges}\): the number of cellulation (triangulation) edges added.
An inductive type classifying the different types of qubits involved in the gauging measurement protocol:
LogicalSupport: Original code qubits in the support of \(L\).
EdgeQubit: Auxiliary gauge qubits on edges, initialized in \(\ket {0}\).
DummyQubit: Auxiliary qubits for dummy vertices, initialized in \(\ket {+}\).
An inductive type representing the possible initial quantum states for qubits:
zero: The \(\ket {0}\) state (computational basis zero).
plus: The \(\ket {+}\) state, i.e., \((\ket {0} + \ket {1})/\sqrt{2}\).
encoded: Unspecified state (for logical qubits, determined by the encoding).
The function assigning an initial state to each qubit type:
\(\mathrm{LogicalSupport} \mapsto \mathrm{encoded}\) (state determined by the code encoding),
\(\mathrm{EdgeQubit} \mapsto \mathrm{zero}\) (initialized in \(\ket {0}\)),
\(\mathrm{DummyQubit} \mapsto \mathrm{plus}\) (initialized in \(\ket {+}\)).
A round requirement specifies how many error correction rounds are needed before and after gauging. It consists of:
a code distance \(d \in \mathbb {N}\) with \(d {\gt} 0\),
the number of rounds before gauging, \(r_{\mathrm{before}} \in \mathbb {N}\),
the number of rounds after gauging, \(r_{\mathrm{after}} \in \mathbb {N}\),
a requirement type (theoretical or practical).
The edge type for the Shor-style gauging graph with path connectivity on dummy vertices is an inductive type with two constructors:
\(\mathrm{cross}(i)\) for \(i \in \mathrm{Fin}(W)\): a cross edge connecting \((\mathtt{false}, i)\) to \((\mathtt{true}, i)\),
\(\mathrm{dummyPath}(i)\) for \(i \in \mathrm{Fin}(W-1)\): a path edge connecting \((\mathtt{true}, i)\) to \((\mathtt{true}, i+1)\).
The edge type for the Shor-style gauging graph with star connectivity on dummy vertices is an inductive type with two constructors:
\(\mathrm{cross}(i)\) for \(i \in \mathrm{Fin}(W)\): a cross edge connecting \((\mathtt{false}, i)\) to \((\mathtt{true}, i)\),
\(\mathrm{dummyStar}(i)\) for \(i \in \mathrm{Fin}(W-1)\): a star edge connecting \((\mathtt{true}, 0)\) to \((\mathtt{true}, i+1)\).
Two vertices \(x, y \in \mathrm{ShorVertex}(W)\) are adjacent in the path-connectivity graph if one of the following holds:
\(x_1 \neq y_1\) and \(x_2 = y_2\) (cross edge between support and dummy partner),
\(x_1 = \mathtt{true}\), \(y_1 = \mathtt{true}\), and \(x_2 + 1 = y_2\) (forward path edge),
\(x_1 = \mathtt{true}\), \(y_1 = \mathtt{true}\), and \(y_2 + 1 = x_2\) (backward path edge).
The Shor-style path graph as a\(\mathrm{GraphWithCycles}\) instance with cycle type \(\mathrm{Empty}\) (since the graph is a tree with \(0\) independent cycles). The instance consists of:
Graph: \(\mathrm{shorPathGraph}(W)\),
Edge endpoints: \(\mathrm{shorPathEdgeEndpoints}(W)\),
Edge adjacency and symmetry from the established theorems,
Cycles: the empty function from \(\mathrm{Empty}\) (no cycles).
Two vertices \(x, y \in \mathrm{ShorVertex}(W)\) are adjacent in the star-connectivity graph if one of the following holds:
\(x_1 \neq y_1\) and \(x_2 = y_2\) (cross edge),
\(x_1 = y_1 = \mathtt{true}\), \(x_2 = 0\), and \(y_2 {\gt} 0\) (outgoing star edge from center),
\(x_1 = y_1 = \mathtt{true}\), \(y_2 = 0\), and \(x_2 {\gt} 0\) (incoming star edge to center).
Given a positive integer \(W = |\operatorname {supp}(L)|\), the vertex type for the Shor-style gauging graph is \(\mathrm{Bool} \times \mathrm{Fin}(W)\), where:
\((\mathtt{false}, i)\) represents the \(i\)-th support vertex (qubit \(i\) in \(\operatorname {supp}(L)\)),
\((\mathtt{true}, i)\) represents the \(i\)-th dummy vertex \(u_i\).
Desideratum 1: Short Paths Property.
Let \(G = (V, E)\) be a simple graph, let \(\mathcal{Z}\) be a collection of \(Z\)-type support sets (each a finite subset of \(V\)), and let \(k \in \mathbb {N}\) be a constant path length bound.
The short paths property \(\mathrm{ShortPaths}(G, \mathcal{Z}, k)\) holds if: for every \(Z\)-type support set \(S \in \mathcal{Z}\) and every pair of vertices \(u, v \in S\), there exists a walk from \(u\) to \(v\) in \(G\) of length at most \(k\):
This ensures that deformed checks have bounded weight.
A space-fault over qubit types \(V, E\) is a structure consisting of:
qubit : \(\texttt{QubitLoc}(V, E)\) — the qubit where the error occurs,
time : \(\mathbb {N}\) — the time step when the error occurs,
pauliType : \(\texttt{PauliType}\) — the type of Pauli error (\(I\), \(X\), \(Y\), or \(Z\)).
Given two space-faults \(f\) and \(g\) at the same qubit and time (i.e., \(f.\texttt{qubit} = g.\texttt{qubit}\) and \(f.\texttt{time} = g.\texttt{time}\)), their composition is:
The set of bulk detectors during deformation at time \(t\) (used for \(t_i {\lt} t {\lt} t_o\)) is the union:
The set of bulk detectors for original checks at time \(t\) (used for \(t {\lt} t_i\) and \(t {\gt} t_o\)) is
The time classification of a detector:
\(\mathrm{bulk}\): repeated measurement of the same observable at times \(t - \tfrac {1}{2}\) and \(t + \tfrac {1}{2}\),
\(\mathrm{initialBoundary}\): initialization at \(t_i - \tfrac {1}{2}\), first measurement at \(t_i + \tfrac {1}{2}\),
\(\mathrm{finalBoundary}\): last measurement at \(t_o - \tfrac {1}{2}\), readout at \(t_o + \tfrac {1}{2}\).
The final boundary detectors at \(t = t_o\) are:
A gauging region is a structure specifying the time boundaries of the gauging measurement procedure, consisting of:
A start time \(t_i\) (the start of code deformation),
An end time \(t_o\) (the end of code deformation),
The validity condition \(t_i {\lt} t_o\), ensuring the deformation has positive duration.
The initial boundary detectors at \(t = t_i\) are:
The type of operator involved in a detector, given as an inductive type with constructors:
\(\mathrm{originalCheck}(j)\): the original stabilizer check \(s_j\),
\(\mathrm{gaussLaw}(v)\): the Gauss law operator \(A_v\),
\(\mathrm{flux}(p)\): the flux operator \(B_p\),
\(\mathrm{deformedCheck}(j)\): the deformed check \(\tilde{s}_j\),
\(\mathrm{edgeZ}(e)\): a single-qubit \(Z\) measurement on edge \(e\).
A cleaning sequence is a list of Pauli pair moves that collectively move all Pauli errors to a target time \(t_{\mathrm{ref}}\). The combined effect is:
Space errors: net effect is to relocate all Paulis to \(t_{\mathrm{ref}}\).
Time errors: XOR of all induced measurement faults.
The combined measurement errors of a cleaning sequence is the function \(M \to \mathrm{TimeStep} \to \mathrm{Bool}\) obtained by folding over the list of moves: for each measurement \(m\) and time \(t\), the result is the XOR of all induced measurement faults from moves whose source time or successor time equals \(t\).
Two spacetime faults \(F\) and \(G\) are equivalent modulo stabilizers (with respect to a detector collection \(\mathrm{DC}\), base outcomes, and logical effect predicate) if there exists a spacetime stabilizer \(S\) such that \(F = G \cdot S\).
A spacetime fault \(F\) is a pure space fault at a single time step \(t\) if:
All space errors occur only at time \(t\): for all qubits \(q\) and times \(t' \neq t\), the space error \(F.\mathrm{spaceErrors}(q, t') = I\).
There are no time errors: for all measurements \(m\) and times \(t'\), \(F.\mathrm{timeErrors}(m, t') = \mathrm{false}\).
A Pauli pair move captures a single operation of moving a Pauli error from time \(t\) to \(t+1\) using a Pauli pair stabilizer. It consists of:
A qubit location (vertex or edge).
The source time step (fromTime).
The Pauli type (\(X\), \(Y\), or \(Z\)).
The induced measurement faults on anticommuting checks.
A spacetime fault over qubit types \(V, E\) and measurement type \(M\) is a structure consisting of:
\(\texttt{spaceErrors} : \texttt{QubitLoc}(V, E) \to \mathbb {N} \to \texttt{PauliType}\) — the Pauli error at each (qubit, time) location; identity means no error.
\(\texttt{timeErrors} : M \to \mathbb {N} \to \texttt{Bool}\) — whether each (measurement, time) has a flipped outcome; \(\text{false}\) means no error.
This functional representation allows efficient pointwise composition.
The composition of spacetime faults \(f\) and \(g\) is defined pointwise:
Given an initialization fault \(f\), the corresponding spacetime fault is obtained by first converting \(f\) to its equivalent space-fault via \(\operatorname {toSpaceFault}\), then embedding it via \(\operatorname {fromSpaceFault}\).
Given a single space-fault \(f\), the corresponding spacetime fault has:
and \(\texttt{timeErrors}(m, t) = \text{false}\) for all \(m, t\).
Given a single time-fault \(f\), the corresponding spacetime fault has \(\texttt{spaceErrors}(q, t) = I\) for all \(q, t\), and:
The inverse of a spacetime fault \(f\) is:
since Pauli operators are self-inverse and XOR is self-inverse.
The space-only component of a spacetime fault \(f\) retains the space errors and sets all time errors to \(\text{false}\):
The time-only component of a spacetime fault \(f\) sets all space errors to \(I\) and retains the time errors:
The spacetime fault-distance \(d_{\mathrm{ST}}\) is defined as:
where \(W = \mathrm{logicalFaultWeights}(\mathit{DC}, \mathit{times})\). When logical faults exist, the minimum is computed using the well-foundedness of \({\lt}\) on \(\mathbb {N}\) (via \(\mathrm{WellFounded.min}\)). The value \(0\) serves as a sentinel when no logical faults exist.
A cleaning preserves parity structure, for a fault \(F\) and its cleaned version \(F_{\mathrm{cleaned}}\), asserts:
Equivalence: There exists a spacetime stabilizer \(S\) such that \(F_{\mathrm{cleaned}} = F \cdot S\).
Parity preservation: For each spatial position \(q\), there exists \(k \in \mathbb {N}\) such that for any set of times, the number of non-identity Pauli faults in \(F_{\mathrm{cleaned}}\) plus the number in \(F\) at position \(q\) equals \(2k\) (i.e., has constant parity).
A fault distance configuration is a structure consisting of:
A code distance \(d \in \mathbb {N}\) with \(d {\gt} 0\),
An initial time \(t_i\) and final time \(t_o\) with \(t_i {\lt} t_o\) (the deformation interval is nonempty),
The condition \((t_o - t_i) \geq d\) (sufficient measurement rounds),
A Cheeger constant \(h(G) \in \mathbb {R}\) with \(h(G) \geq 1\).
An original logical operator applied at a specific time consists of:
A time step \(t\) (intended to be outside the deformation region),
A Pauli assignment \(\mathrm{vertexPaulis} : V \to \mathrm{PauliType}\) on each vertex qubit,
A weight \(w \in \mathbb {N}\) equal to \(|\{ v : \mathrm{vertexPaulis}(v) \neq I\} |\).
Given an original logical operator \(L\) applied at time \(t\), define the corresponding spacetime fault:
Space errors: at qubit \(q\) and time \(t'\), the error is \(L.\mathrm{vertexPaulis}(v)\) if \(q\) is a vertex \(v\) and \(t' = t\), and is \(I\) otherwise (including all edge qubits).
Time errors: identically false (no measurement errors).
An inductive type capturing the decomposition from the Spacetime Decoupling Lemma (Lem 6) and the resulting case analysis. For a spacetime logical fault \(F\):
Case 1 (timeNontrivial): There exist \(F_{\mathrm{space}}\), \(F_{\mathrm{time}}\), and a stabilizer \(S\) with \(F = F_{\mathrm{space}} \cdot F_{\mathrm{time}} \cdot S\), where \(F_{\mathrm{time}}\) is a nontrivial spacetime logical fault that spans the interval, and \(|F| \geq |F_{\mathrm{time}}|\).
Case 2 (spaceOnly): There exists \(F_{\mathrm{space}}\) which is a pure space fault at time \(t_i\) and a stabilizer \(S\) with \(F = F_{\mathrm{space}} \cdot S\), where \(|F_{\mathrm{space}}| \geq d\) and \(|F| \geq |F_{\mathrm{space}}|\).
A structure asserting that a time fault \(F_{\mathrm{time}}\) spans the interval \([t_i, t_o)\):
For every \(t\) with \(t_i \leq t {\lt} t_o\), there exists a measurement \(m\) such that \(F_{\mathrm{time}}\) has a time error at \((m, t)\).
\(F_{\mathrm{time}}\) is a pure time fault (no space errors).
A detector effect models the effect of a generator’s faults on a detector. It is a structure with three \(\mathbb {Z}_2\)-valued fields:
\(\mathtt{pauliFlip\_ t} \in \mathbb {Z}_2\): flip from a Pauli at time \(t\) affecting the measurement at \(t + \tfrac {1}{2}\),
\(\mathtt{pauliFlip\_ t1} \in \mathbb {Z}_2\): flip from a Pauli at time \(t+1\) affecting the measurement at \(t + \tfrac {3}{2}\),
\(\mathtt{measFault} \in \mathbb {Z}_2\): flip from a measurement fault at \(t + \tfrac {1}{2}\).
The net effect on detector \(c^t\) (comparing measurements at \(t - \tfrac {1}{2}\) and \(t + \tfrac {1}{2}\)) is defined as
The measurement at \(t - \tfrac {1}{2}\) is unaffected, while the recorded outcome at \(t + \tfrac {1}{2}\) is flipped by both the Pauli at time \(t\) and the measurement fault.
The net effect on detector \(c^{t+1}\) (comparing measurements at \(t + \tfrac {1}{2}\) and \(t + \tfrac {3}{2}\)) is
A generator \(g\) has empty syndrome if, for every relevant detector, the net \(\mathbb {Z}_2\) effect is zero. Specifically:
For spaceStabilizer: \(\mathrm{netEffect}_{c^t} = 0\) and \(\mathrm{netEffect}_{c^{t+1}} = 0\) using the space stabilizer effect;
For Pauli pair generators (pauliPairOriginal, vertexXPair, vertexZPair, edgeXPair, edgeZPair): both anticommuting and commuting detector effects vanish at \(c^t\) and \(c^{t+1}\);
For initFaultPlusXe: \(\mathrm{netEffect}_{c^t} = 0\);
For boundary \(Z_e\) pair generators: both \(c^t\) and \(c^{t+1}\) effects vanish;
For finalBoundaryXePair: \(\mathrm{netEffect}_{c^t} = 0\);
For finalBoundaryBareZe: \(\mathrm{netEffect}_{c^t} = 0\).
The function \(\mathrm{generatorToPattern}\) converts each generator type to an explicit fault pattern. For example:
spaceStabilizer: empty Pauli faults, empty measurement faults;
pauliPairOriginal(\(q\), \(P\), \(t\), \(S\)): Pauli faults \(\{ (q, t, P), (q, t{+}1, P)\} \), measurement faults \(\{ (s_j, t) : j \in S\} \);
vertexXPair(\(v\), \(t\), \(S\)): Pauli faults \(\{ (\mathrm{vertex}(v), t, X), (\mathrm{vertex}(v), t{+}1, X)\} \) with deformed check measurement faults;
initFaultPlusXe(\(e\), \(t\)): single Pauli fault \(\{ (\mathrm{edge}(e), t, X)\} \);
and analogously for all other generator types.
The predicate \(\mathrm{generatorValidInRegion}(\mathrm{gen}, \mathrm{region})\) specifies which time regions a generator is valid in:
spaceStabilizer with region \(r\): valid only in region \(r\);
pauliPairOriginal: valid in beforeGauging or afterGauging;
vertexXPair, vertexZPair, edgeXPair: valid in duringGauging or initialBoundary;
edgeZPair: valid in duringGauging;
initFaultPlusXe, initialBoundaryZePair: valid in initialBoundary;
finalBoundaryXePair, finalBoundaryBareZe, finalBoundaryZePair: valid in finalBoundary.
The logical effect of a generator \(g\) is a value in \(\mathbb {Z}_2\) (\(0\) = trivial, \(1\) = nontrivial). For every generator type, \(\mathrm{logicalEffect}(g) = 0\):
Stabilizers act as the identity on the code space;
Pauli pairs cancel: \(P \cdot P = I\);
Initialization fault plus \(X\) cancels;
Boundary generators act on qubits being initialized or discarded.
A fault pattern \(p\) has empty syndrome if for every detector \((c, t)\):
A fault pattern \(p\) is generated by the listed generators if there exists a list of generators whose patterns product to \(p\):
The net Pauli parity of a fault pattern \(p\) at qubit \(q\) is the number of Pauli faults at \(q\) modulo \(2\):
A value of \(0\) means an even number of Paulis (they cancel: \(P \cdot P = I\)), while \(1\) means a nontrivial Pauli remains.
A fault pattern \(p\) preserves logical information if either:
For every qubit \(q\), \(\mathrm{netPauliParity}(p, q) = 0\) (all Paulis cancel pairwise), or
\(p.\mathrm{pauliFaults} = \emptyset \) (the pattern is a pure space stabilizer with no explicit Pauli faults).
A spacetime stabilizer generator is one of the following types:
For \(t {\lt} t_i\) and \(t {\gt} t_o\) (original code):
spaceStabilizer: a check operator \(s_j\) (or \(\tilde{s}_j\), \(A_v\), \(B_p\)) at time \(t\);
pauliPairOriginal: a Pauli \(P\) at time \(t\) and \(P\) at time \(t+1\), together with measurement faults on all anticommuting checks \(s_j\).
For \(t_i {\lt} t {\lt} t_o\) (deformed code):
vertexXPair: \(X_v\) at \(t\) and \(t+1\) with measurement faults on anticommuting \(\tilde{s}_j\);
vertexZPair: \(Z_v\) at \(t\) and \(t+1\) with measurement faults on \(A_v\) and anticommuting \(\tilde{s}_j\);
edgeXPair: \(X_e\) at \(t\) and \(t+1\) with measurement faults on \(B_p\) (for \(p \ni e\)) and anticommuting \(\tilde{s}_j\);
edgeZPair: \(Z_e\) at \(t\) and \(t+1\) with measurement faults on \(A_v\) (for \(v \in e\)).
For \(t = t_i\) (initial boundary):
initFaultPlusXe: initialization fault on \(|0\rangle _e\) plus \(X_e\) at time \(t\);
initialBoundaryZePair: \(Z_e\) at \(t+1\) with \(A_v\) measurement faults for \(v \in e\).
For \(t = t_o\) (final boundary):
finalBoundaryXePair: \(X_e\) at \(t\) with \(Z_e\) measurement fault;
finalBoundaryBareZe: bare \(Z_e\) at \(t\) (since \(Z|0\rangle = |0\rangle \));
finalBoundaryZePair: \(Z_e\) at \(t-1\) with \(A_v\) measurement faults for \(v \in e\).
A time region in the gauging measurement procedure is one of:
beforeGauging: \(t {\lt} t_i\) (original code),
duringGauging: \(t_i {\lt} t {\lt} t_o\) (deformed code),
afterGauging: \(t {\gt} t_o\) (original code),
initialBoundary: \(t = t_i\),
finalBoundary: \(t = t_o\).
Under the conditions that the logical effect predicate is group-like and the syndrome computation respects the group structure, the spacetime stabilizers form a subgroup of \(\mathrm{SpacetimeFault}(V, E, M)\).
The carrier is the spacetime stabilizer set. Closure under multiplication, the identity element, and closure under inverses are verified using the respective group-like and syndrome homomorphism properties.
A space-time tradeoff for gauging measurement consists of:
A code distance \(d \in \mathbb {N}\),
A tradeoff parameter \(m \in \mathbb {N}\) with \(m {\gt} 0\),
A proof that \(m \mid d\).
One performs \(2m - 1\) measurements of equivalent logical operators in parallel, for \(d/m\) rounds, and takes a majority vote to determine the classical outcome.
An \([[n, k, d]]\) stabilizer code is a structure consisting of:
the number of physical qubits \(n \in \mathbb {N}\),
the number of logical qubits \(k \in \mathbb {N}\),
the code distance \(d \in \mathbb {N}\),
a stabilizer group on \(n\) qubits,
a constraint \(n \geq k\),
a positivity constraint \(d {\gt} 0\).
A Pauli operator \(P\) is in the stabilizer group of code \(C\) if there exists a finite set of generators \(\texttt{gens} \subseteq C.\texttt{stabilizers}.\texttt{generators}\) such that:
if \(\texttt{gens} = \emptyset \) then \(P\) equals the identity, and
\(P\) is one of the generators, or \(P\) is the identity.
(This is a simplified definition; a full treatment would define the group generated by the generators.)
A stabilizer group on \(n\) qubits is a structure consisting of:
a set of generating Pauli operators \(\texttt{generators} \subseteq \texttt{PauliOp}(n)\),
a proof that the generating set is finite,
a proof that all generators mutually commute: for all \(s_1, s_2 \in \texttt{generators}\), \(\operatorname {commutes}(s_1, s_2)\).
Two Pauli types \(p, q \in \mathrm{StabPauliType}\) are compatible if they do not conflict on the same qubit. A conflict occurs when one is purely \(X\)-type and the other is purely \(Z\)-type. Formally,
Multiplication of single-qubit Pauli types (ignoring global phase) is defined by the rules:
\(I \cdot p = p\) and \(p \cdot I = p\) for all \(p\),
\(X \cdot X = Y \cdot Y = Z \cdot Z = I\),
\(X \cdot Y = Z\), \(Y \cdot X = Z\), \(Y \cdot Z = X\), \(Z \cdot Y = X\), \(Z \cdot X = Y\), \(X \cdot Z = Y\).
Given a vertex type \(V\) and an edge type \(E\), the subdivided vertex type is the disjoint sum \(V \oplus E\), written \(\texttt{SubdividedVertex'}(V, E) := V \sqcup E\). Original vertices are embedded via \(\mathrm{inl}\) and dummy vertices (one per edge) via \(\mathrm{inr}\).
Let \(G\) be a simple graph on \(V\). The subdivided graph \(\mathrm{subdivideGraph}'(G)\) is the simple graph on \(\texttt{SubdividedVertex'}(V, \mathrm{Sym2}(V))\) with adjacency given by \(\mathrm{subdivisionAdj}'(G)\). The symmetry and irreflexivity axioms are verified as follows:
Symmetry: We case-split on \(x\) and \(y\). If \(x = \mathrm{inl}(v_1)\) and \(y = \mathrm{inr}(e)\), then the adjacency hypothesis directly gives the reverse. The cases \(\mathrm{inl}\)–\(\mathrm{inl}\) and \(\mathrm{inr}\)–\(\mathrm{inr}\) are contradictory by definition.
Irreflexivity: For any \(x\), we case-split; in both cases (\(\mathrm{inl}\) or \(\mathrm{inr}\)), the adjacency relation simplifies to \(\mathrm{False}\).
Let \(G\) be a simple graph on \(V\) with decidable adjacency. The subdivision adjacency on \(\texttt{SubdividedVertex'}(V, \mathrm{Sym2}(V))\) is defined by:
That is, a vertex \(\mathrm{inl}(v)\) is adjacent to \(\mathrm{inr}(e)\) if and only if \(v\) is an endpoint of the edge \(e\). There are no \(\mathrm{inl}\)–\(\mathrm{inl}\) or \(\mathrm{inr}\)–\(\mathrm{inr}\) adjacencies.
Desideratum 2: Sufficient Expansion.
A simple graph \(G = (V, E)\) has sufficient expansion if it is a strong expander, i.e., its Cheeger constant satisfies \(h(G) \geq 1\). Formally, \(\mathrm{SufficientExpansion}(G)\) is defined to be the property \(\mathrm{IsStrongExpander}(G)\).
Given a detector collection \(\mathrm{DC}\), base outcomes, and a spacetime fault \(F\), the syndrome is:
The syndrome computation respects the group structure if:
The identity has empty syndrome.
For all \(f, g\), if both have empty syndrome, then \(f \cdot g\) has empty syndrome.
For all \(f\), if \(f\) has empty syndrome, then \(f^{-1}\) has empty syndrome.
A syndrome linearity condition for a detector collection \(\mathrm{DC}\) and base outcomes consists of two properties:
Local dependence: For all spacetime faults \(f, g\) and all detectors \(D \in \mathrm{DC}.\mathrm{detectors}\), if \(f\) and \(g\) agree on all time-errors at \(D\)’s measurement events, then \(D\) is violated by \(f\) if and only if it is violated by \(g\).
Parity determines violation: For all faults \(f\) and detectors \(D \in \mathrm{DC}.\mathrm{detectors}\), detector \(D\) is violated if and only if the number of measurement events of \(D\) at which \(f\) has a time-error is odd:
\[ D.\mathrm{isViolated}(\mathrm{applyFaultToOutcomes'}(\mathrm{base}, f)) \iff |\{ e \in D.\mathrm{measEvents} \mid f.\mathrm{timeErrors}(e.\mathrm{measurement}, e.\mathrm{time})\} | \equiv 1 \pmod{2}. \]
The syndrome as a function to \(\mathbb {Z}/2\mathbb {Z}\): for each detector \(D\),
A time-fault (measurement error) over measurement type \(M\) is a structure consisting of:
measurement : \(M\) — the measurement that is affected,
time : \(\mathbb {N}\) — the time step when the measurement occurs,
isFlipped : \(\texttt{Bool}\) — whether the classical outcome is flipped (\(\text{true}\) means error).
Given two time-faults \(f\) and \(g\) at the same measurement and time, their composition is:
where \(\oplus \) denotes XOR (exclusive or).
The \(A_v\) measurement fault chain for check \(m\) and deformation interval \(I\) is the spacetime fault defined by:
\(\mathrm{spaceErrors}(q, t) = I\) (identity) for all qubits \(q\) and times \(t\),
\(\mathrm{timeErrors}(m', t) = \mathrm{true}\) if and only if \(m' = m\) and \(t_i \leq t {\lt} t_o\).
This represents measurement faults on a single \(A_v\) check at all times in the interval \([t_i, t_o)\).
A fault \(F\) has a gauging logical effect on interval \(I\) if there exists a check \(m\) such that the number of time errors of \(F\) on \(m\) over the interval \([t_i, t_o)\) has odd cardinality:
This captures the idea that flipping an odd number of \(A_v\) outcomes changes \(\sigma = \prod _v \varepsilon _v\).
The set of pure time logical faults for a detector collection \(DC\), base outcomes, and logical effect predicate is:
The pure time spacetime fault distance is the minimum weight over all pure time logical faults, defined using the well-founded minimum on \(\mathbb {N}\):
A spacetime stabilizer generator is one of three types:
Initial: “init fault \(+\) \(X_e\) at \(t_i\)”, parameterized by an edge.
Propagator: “\(X_e\) at \(t\), \(X_e\) at \(t+1\), measurement faults between”, parameterized by an edge and time.
Final: “\(X_e\) at \(t_o\) \(+\) \(Z_e\) readout fault”, parameterized by an edge.
The canonical decomposition of a Type 2 fault string \(s\) uses:
Initial generator: \(\mathrm{initial}(s.\mathrm{edge})\),
Propagator at index \(i\): \(\mathrm{propagator}(s.\mathrm{edge},\, s.I.t_i + i)\),
Final generator: \(\mathrm{final}(s.\mathrm{edge})\).
A Type 2 decomposition of a Type 2 fault string \(s\) consists of:
An initial generator (of type \(\mathrm{initial}(s.\mathrm{edge})\)),
A sequence of \(s.I.\mathrm{numRounds} - 1\) propagator generators (each of type \(\mathrm{propagator}(s.\mathrm{edge}, t)\) for some \(t\)),
A final generator (of type \(\mathrm{final}(s.\mathrm{edge})\)).
A fault \(F\) violates comparison detector \(D = (c, t)\) if the parity of the fault count at \((c, t)\) differs from the parity at \((c, t-1)\):
When \(t = 0\), the previous count is taken to be \(0\).
Given a list of cycle vertices \([v_0, v_1, \ldots , v_{n-1}]\), the triangulation edge pairs are the diagonal edges that triangulate the cycle using a zigzag pattern. For a cycle of length \(n \geq 4\), this produces \(n - 3\) triangulation edges following the pattern:
Even-indexed steps connect from the low index to the high end, and odd-indexed steps connect from the high end to the next low index. For cycles of length \({\lt} 4\), the result is the empty list.
The trivial extension of an original code operator to the deformed system. Given an operator \(\mathrm{op}\) with support only on vertices, the trivial extension \(\mathrm{op} \otimes I_E\) is defined as a deformed Pauli operator with:
\(X\)-support on vertices: \(\mathrm{op}.\mathrm{xSupport}\),
\(Z\)-support on vertices: \(\mathrm{op}.\mathrm{zSupport}\),
\(X\)-support on edges: \(\emptyset \),
\(Z\)-support on edges: \(\emptyset \),
phase: \(0\).
A trivially extendable logical for a graph \(G\) with cycles and an original code distance specification consists of:
An original logical operator \(\mathrm{op}\),
A proof that \(\mathrm{op}\) is a logical of the original code (\(\mathrm{code}.\mathrm{isLogical}(\mathrm{op})\)),
A proof that \(\mathrm{op}\) is nontrivial (\(\mathrm{op}.\mathrm{isNontrivial}\)),
A proof that \(\mathrm{op}\) has the minimum weight \(d\) (\(\mathrm{op}.\mathrm{weight} = \mathrm{code}.d\)),
A proof that the trivial extension commutes with all Gauss law operators \(A_v\): for all vertices \(v\), \(|\mathrm{op}.\mathrm{zSupport} \cap \{ v\} | \equiv 0 \pmod{2}\).
The worst-case construction specification is the proposition:
More precisely, assuming expander existence and the Freedman–Hastings decongestion lemma, for every \(W \ge 2\) there exist types \(V\), \(E\), \(C\) (with decidable equality and finiteness) and a \(\mathrm{GraphWithCycles}\) structure \(G\) on \((V, E, C)\) together with a family of \(Z\)-type supports \(\mathrm{zTypeSupports}\), such that:
Desideratum 1 (Short paths): There exists \(\kappa \in \mathbb {N}\) such that \(\mathrm{ShortPathsProperty}(G.\mathrm{graph},\, \mathrm{zTypeSupports},\, \kappa )\) holds.
Desideratum 2 (Sufficient expansion): \(\mathrm{SufficientExpansion}(G.\mathrm{graph})\) holds.
Desideratum 3 (Low-weight cycles): \(\mathrm{LowWeightCyclesProperty}(G, 3)\) holds (all generating cycles have weight \(\le 3\)).
Overhead bound: There exists \(C \in \mathbb {N}\) with \(|V| \le W \cdot (C \cdot (\log _2 W)^2 + C + 1)\), i.e., \(O(W \log ^2 W)\) vertices.
An \(X\)-measurement specification for a qubit type \(Q\) consists of:
a finite set \(\mathcal{Q} \subseteq Q\) of qubits to measure in the Pauli-\(X\) basis,
a function \(\mathrm{keep} : Q \to \mathrm{Bool}\) indicating whether to keep each qubit after measurement.
An \(X\)-type logical operator for a stabilizer code \(C\) extends a logical operator with the additional property that the underlying Pauli operator is purely \(X\)-type. This captures the convention
Given a stabilizer code \(C\) and a finite set \(S \subseteq \operatorname {Fin}(C.n)\), together with proofs that \(\operatorname {pureX}(S)\) commutes with all stabilizers and is not in the stabilizer group, we construct an \(X\)-type logical operator with underlying operator \(\operatorname {pureX}(S)\). The purely \(X\)-type property follows from the lemma that \(\operatorname {pureX}(S)\) is purely \(X\)-type.
A zigzag state tracks the progress of building triangulation edges for a cycle. It consists of:
\(\mathrm{low} : \mathbb {N}\), the current low index,
\(\mathrm{high} : \mathbb {N}\), the current high index (exclusive),
\(\mathrm{fromLow} : \mathrm{Bool}\), indicating whether the next edge starts from the low index.
One step of the zigzag edge generation. Given a list of cycle vertices of length \(n\), and a zigzag state \(s\) with \(s.\mathrm{low} {\lt} n\), \(s.\mathrm{high} \leq n\), and \(s.\mathrm{low} + 2 {\lt} s.\mathrm{high}\), the step produces an edge pair \((v_1, v_2)\) and a new zigzag state:
If \(\mathrm{fromLow}\) is true: the edge connects \(v_{\mathrm{low}}\) to \(v_{\mathrm{high}-1}\), and the next state has \(\mathrm{high}' = \mathrm{high} - 1\) and \(\mathrm{fromLow}' = \mathrm{false}\).
If \(\mathrm{fromLow}\) is false: the edge connects \(v_{\mathrm{high}-1}\) to \(v_{\mathrm{low}+1}\), and the next state has \(\mathrm{low}' = \mathrm{low} + 1\) and \(\mathrm{fromLow}' = \mathrm{true}\).
A \(Z\)-measurement specification for a qubit type \(Q\) consists of:
a finite set \(\mathcal{Q} \subseteq Q\) of qubits to measure in the Pauli-\(Z\) basis,
a function \(\mathrm{keep} : Q \to \mathrm{Bool}\) indicating whether to keep each qubit after measurement.
If a time fault \(f\) has \(f.\mathrm{isFlipped} = \mathrm{true}\), and a measurement event \(e\) has \(e.\mathrm{measurement} = f.\mathrm{measurement}\) and \(e.\mathrm{time} = f.\mathrm{time}\), then
Let \(G\) be a graph with cycles, let \(w \in \mathbb {N}\), and suppose the low-weight cycles property \(\mathrm{LowWeightCycles}(G, w)\) holds. Then for every cycle \(c \in C\), the flux operator \(B_c = \prod _{e \in c} Z_e\) has weight at most \(w\):
This is immediate from the definition of flux operators.
The sum (product) of edge supports of two flux operators corresponds to the symmetric difference. For any edge \(e\):
A fault \(f\) has empty syndrome if and only if every detector \(D\) in the collection is satisfied by the faulted outcomes:
Let \(F\) be a spacetime fault, \(I\) a deformation interval, and \(m\) a check. If no interior comparison detector on \(m\) is violated, and there exists \(t_0 \in [t_i, t_o)\) with \(\mathrm{Odd}(\mathrm{timeFaultCountAt}(F, m, t_0))\), then for all \(t \in [t_i, t_o)\), \(\mathrm{timeFaultCountAt}(F, m, t) {\gt} 0\).
Let \(F\) be a spacetime fault and \(I\) a deformation interval. If no comparison detector \((m, t)\) is violated for any interior time \(t_i {\lt} t {\lt} t_o\), then for all \(t_1, t_2 \in [t_i, t_o)\):
The trivial extension of an operator with empty \(X\)-support on edges is a \(1\)-cocycle. Every cycle has even (zero) intersection with the empty edge set:
A trivially extended nontrivial operator is nontrivial in the deformed code. If \(\mathrm{op}\) is nontrivial (i.e., \(\mathrm{op}.\mathrm{isNontrivial}\) holds), then \(\mathrm{trivialExtension}(\mathrm{op}).\mathrm{isNontrivial}\) holds.
Let \(\textit{conv}\) be a boundary condition convention. For every error process \(\textit{ep}\) such that \(\texttt{spansToBoundary}(\textit{ep}) = \texttt{true}\), under the propagation model (spanning to the initial boundary requires weight \(\geq \texttt{roundsBefore} + 1\), and similarly for the final boundary), we have
Let \(\textit{conv}\) be a boundary condition convention with code distance \(d\). For every error process \(\textit{ep}\) with \(\texttt{weight} \leq d\), under the propagation model, the error cannot span to any boundary:
This captures the fact that boundary quality is irrelevant for low-weight errors.
Let \(\textit{ec}\) be an error correction configuration providing full protection, and let \(\textit{ep}\) be an error process with \(\texttt{weight} {\gt} 0\) that spans to a boundary. Suppose the propagation model holds: if the process involves gauging and the initial boundary then \(\texttt{weight} \geq \texttt{roundsBefore} + 1\), and similarly for the final boundary. Then
The cardinality of a finset \(S\) can be recovered from its characteristic vector:
where \(\mathrm{val} : \mathbb {Z}/2\mathbb {Z} \to \mathbb {N}\) is the distance-minimal representative.
For finsets \(S, T \subseteq \alpha \),
Addition of characteristic vectors corresponds to symmetric difference. This is the fundamental property that allows linear algebra over \(\mathbb {Z}/2\mathbb {Z}\) to encode set operations.
When \(h(G) {\gt} 1\), the bound is still \(d^* \geq d\) (not better). Specifically, for any deformed logical operator \(L\):
The \(\min (h(G), 1)\) factor caps the improvement at \(d\).
When \(0 {\lt} h(G) {\lt} 1\), the lower bound from Lemma 2 gives \(|L| \geq h(G) \cdot d\), and furthermore \(h(G) \cdot d {\lt} d\), confirming that the distance is reduced. Formally:
Part 1 Main Result: When \(h(G) = 1\), the deformed code distance exactly equals the original code distance:
The equality follows from matching upper and lower bounds: the lower bound \(d^* \geq d\) comes from Lemma 2, and the upper bound \(d^* \leq d\) comes from the existence of a trivially extended logical operator of weight exactly \(d\) in the set of deformed logicals.
When the Cheeger constant is exactly 1, the deformed code distance satisfies \(d^* \geq d\) and there exists a deformed logical with weight exactly \(d\) (the trivial extension of an original minimum-weight logical). That is:
Main Optimality Theorem: \(h(G) = 1\) is optimal for code distance preservation. Given a graph \(G\) with exactness condition, Cheeger constant \(h(G) {\gt} 0\), original code distance \(d\), and a nonempty set of deformed logicals satisfying the Cheeger and weight bound hypotheses:
When \(h(G) \geq 1\): \(d^* \geq d\) (distance preserved).
When \(h(G) {\lt} 1\): \(d^* \geq h(G) \cdot d\) (distance reduced by factor \(h(G)\)).
In all cases: \(d^* \geq \min (h(G), 1) \cdot d\).
Given the exactness condition (i.e., \(\ker (\delta _2) = \operatorname {im}(\delta )\)), every \(1\)-cocycle \(S_X^E\) is a coboundary: there exists a vertex set \(S \subseteq V\) such that
An edge support set \(S_X^E\) is a \(1\)-cocycle if and only if \(|S_X^E \cap p| \equiv 0 \pmod{2}\) for all cycles \(p \in C\):
For good qLDPC codes with \(d = \Theta (n)\) and \(W = O(n)\), we have \(d {\gt} (\log _2 W)^2\), so the gauging measurement provides a significant improvement. Specifically, for \(W \geq 4\) and \(d {\gt} (\log _2 W)^2\), there exists \(c {\gt} 0\) such that
If a gauging measurement occurs in the middle of a large computation with surrounding error correction, then for any constant \(c\) with \(0 {\lt} c \leq d\), the constant round configuration satisfies: it uses constant rounds \(c\), and both the practical rounds before and after are positive.
Assuming \(\mathrm{ExpanderExistenceSpec}\) and \(\mathrm{FreedmanHastingsSpec}\), for every \(W \ge 2\):
Step 1 property: For any graph \(G\) on \(\mathrm{Fin}(W)\) and any adjacent vertices \(u, v\) with \(G.\mathrm{Adj}(u,v)\), there exists a walk of length \(1\) from \(u\) to \(v\).
Overhead arithmetic: For all \(R, C \in \mathbb {N}\) with \(R \le C \cdot (\log _2 W)^2 + C\), we have \(\mathrm{vertexCountWithLayers}(W, R) \le W \cdot (C \cdot (\log _2 W)^2 + C + 1)\).
For every \(X\)-check index \(i \in \mathrm{Fin}(n_X)\), the \(X\)-check vector \(\mathbf{x}_i\) lies in the operator kernel of the initialization hypergraph:
Equivalently, for all \(j \in \mathrm{Fin}(n_Z)\):
The sum of all pairwise \(XX\) operator support vectors equals the all-ones vector:
This represents the fact that the product of all \(XX\) operators gives the logical operator \(L = \prod _v X_v\).
Inter-layer edges only connect consecutive layers: if \((\ell v_1, \ell v_2)\) is an inter-layer edge, then
If \(|F|_{\mathit{times}} {\lt} d_{\mathrm{ST}}\), then either \(F\) does not have empty syndrome (it is detectable) or \(F\) does not affect logical information (it is a stabilizer):
Let \(D\) be a detector and let \(f : M \to \mathrm{TimeStep} \to \mathbb {Z}/2\mathbb {Z}\) be an expected measurement outcome function. If \(D\) is consistent in the sense that
then \(D\) is satisfied under the fault-free outcomes \(\mathrm{faultFreeOutcomes}(f)\).
Let \(D\) be a detector and \(\mathrm{faultedMeas}\) a finite set of faulted measurement events. Suppose:
For every \(e \in D.\mathrm{measEvents}\) with \(e \in \mathrm{faultedMeas}\): \(\mathrm{outcomes}(e.\mathrm{measurement}, e.\mathrm{time}) = 1\),
For every \(e \in D.\mathrm{measEvents}\) with \(e \notin \mathrm{faultedMeas}\): \(\mathrm{outcomes}(e.\mathrm{measurement}, e.\mathrm{time}) = 0\),
\(D.\mathrm{initParity} = 0\),
\(D.\mathrm{expectedParity} = 0\).
Then \(D\) is violated if and only if the fault syndrome equals \(1\):
For codes supporting many disjoint logical representatives, parallel logical gates are possible. Formally, if \(\mathcal{L}\) has pairwise disjoint support, then there exists a \(P : \mathrm{ParallelizableLogicals}(C)\) with \(P.\mathrm{logicals} = \mathcal{L}\) and \(P.\mathrm{overlapBound} = 1\).
The dummy vertex for edge \(\{ u, v\} \) is adjacent to both endpoints in the subdivided graph: if \(G.\mathrm{Adj}(u, v)\), then
If a fault has empty syndrome, then it is either a spacetime stabilizer or a logical fault:
Every fault with empty syndrome is either a spacetime stabilizer or a spacetime logical fault: if \(\mathrm{hasEmptySyndrome}(\mathrm{DC}, \mathrm{base}, f)\), then
A Pauli error at integer time \(t\) affects the state between measurements at \(t - \tfrac {1}{2}\) and \(t + \tfrac {1}{2}\). Formally:
Let \(d {\gt} 0\) be the code distance. If a standard configuration (all checks measured every round) has the fault distance scaling property \(\mathrm{faultDistance} \geq d\), and a sparse configuration (with standard primary checks for \(A_v\) and \(\tilde{s}_j\)) has the same fault distance value, then the sparse configuration also has the scaling property.
Main Theorem (Fault Tolerance). Under the conditions:
\(h(G) \geq 1\) (Cheeger constant at least 1),
\((t_o - t_i) \geq d\) (sufficient measurement rounds),
the spacetime fault-distance equals exactly \(d\):
Assuming that every spacetime logical fault \(F\) admits a lower bound case decomposition, and there exists a weight-\(d\) logical fault, then
Under the hypotheses of the Fault Tolerance Theorem,
Given an original logical operator \(L_{\mathrm{orig}}\) with \(|\{ v \mid L_{\mathrm{orig}}.\mathrm{vertexPaulis}(v) \neq I\} | = d\), applied at time \(L_{\mathrm{orig}}.\mathrm{time} {\lt} t_i\), and assuming \(L_{\mathrm{orig}}.\mathrm{time} \in \mathrm{times}\) and the resulting spacetime fault is a logical fault, there exists a spacetime fault \(F\) that is a spacetime logical fault with \(\mathrm{weight}(F, \mathrm{times}) = d\).
Given that the logical effect is group-like, the syndrome is a group homomorphism, and cleaning exists, any spacetime logical fault \(F\) decomposes as \(F \sim F_{\mathrm{space}} \cdot F_{\mathrm{time}}\), where \(F_{\mathrm{space}}\) is a pure space fault at time \(t_i\) and \(F_{\mathrm{time}}\) is a pure time fault.
For all initialization data init and readout data readout, if all edges are initialized in \(|0\rangle \) (i.e., init.allInitializedInZero = true), then there exists a cycle membership function such that the flux value is well-defined and can be computed from the readout without direct measurement.
Flux \(B_p\) can be inferred if we have: (1) knowledge that edges were initialized in \(|0\rangle \), and (2) final \(Z\)-measurement outcomes on all edges in the cycle. For any readout and cycle membership, the flux value computation is well-defined.
For any gauging circuit \(C\) and any qubit \(q \in \mathrm{allQubits}(C)\), either
That is, every qubit in the circuit is either a vertex qubit or an edge qubit; no additional ancilla qubits are required.
The five phases execute in strict order:
The main theorem of this remark: dummy vertices have no effect on the gauging measurement outcome. For any finite set of dummy vertices, the product of their contributions is \(1\):
For a connected graph \(G\) and any \(z \in (\mathbb {Z}/2\mathbb {Z})^E\), if \(c'\) satisfies \(\delta c' = z\), then for all \(c \in (\mathbb {Z}/2\mathbb {Z})^V\):
That is, the fiber \(\{ c : \delta c = z\} \) has exactly two elements: \(c'\) and \(c' + \mathbf{1}\).
For a connected graph \(G\), outcomes \(\varepsilon \), and \(c'\) with \(\delta c' = z\), let \(c_0 = c'\) and \(c_1 = c' + \mathbf{1}\). Then:
\(\varepsilon (c_0) = \varepsilon (c')\) and \(\varepsilon (c_1) = \varepsilon (c') + \sigma \),
\(X_V(c_0) = X_V(c')\) and \(X_V(c_1) = X_V(c') + L\).
This shows the sum over the fiber \(\{ c : \delta c = z\} \) has exactly two terms whose contributions combine as \(\varepsilon (c') X_V(c')(I + \sigma L)\).
(Gauging Measurement Equivalence.) Let \(G\) be a connected graph with cycles, let \(\varepsilon \) be the Gauss law measurement outcomes, and let \(z \in (\mathbb {Z}/2\mathbb {Z})^E\) be in the image of \(\delta \). Let \(c' = \mathrm{byproductCochain}(G, z)\) and \(\sigma = \sigma (\varepsilon )\). Then:
Fiber structure: For all \(c\), \(\delta c = z \iff (c = c' \text{ or } c = c' + \mathbf{1})\).
Phase relation: \(\varepsilon (c' + \mathbf{1}) = \varepsilon (c') + \sigma \).
Operator relation: \(X_V(c' + \mathbf{1}) = X_V(c') + L\).
Projector properties: \(\sigma + \sigma = 0\) and \(L + L = 0\) (i.e., \(\sigma ^2 = 1\) and \(L^2 = I\)).
Idempotence: \(\sigma \cdot \sigma = \sigma \).
Together, these establish that the gauging procedure is equivalent to projective measurement of \(L\) with eigenvalue \(\sigma \), up to the byproduct operator \(X_V(c')\).
On the \((-\sigma )\)-eigenspace of \(L\) where \(L|\psi _{-\sigma }\rangle = -\sigma |\psi _{-\sigma }\rangle \), the projector \(\frac{1}{2}(I + \sigma L)\) annihilates the state. The key algebraic property is again \(\sigma \cdot \sigma = \sigma \).
On the \(\sigma \)-eigenspace of \(L\) where \(L|\psi _\sigma \rangle = \sigma |\psi _\sigma \rangle \), the projector \(\frac{1}{2}(I + \sigma L)\) acts as the identity. The key algebraic property is \(\sigma \cdot \sigma = \sigma \) in \(\mathbb {Z}/2\mathbb {Z}\).
For a commutative group \(G\), any function \(\mathrm{charges} : \mathrm{Fin}(n) \to G\), and any two permutations \(\sigma _1, \sigma _2\) of \(\mathrm{Fin}(n)\),
For a commutative group \(G\) and any function \(\mathrm{charges} : \mathrm{Fin}(n) \to G\) and any permutation \(\sigma \) of \(\mathrm{Fin}(n)\),
This is the key property that allows local measurements to determine the global charge.
For a nonabelian group \(G\) (i.e., one in which there exist \(a, b \in G\) with \(ab \neq ba\)), there exists a function \(\mathrm{charges} : \mathrm{Fin}(2) \to G\) such that
This is why, for nonabelian groups, local measurements do not fix a definite global charge.
For a (not necessarily abelian) group \(G\), if there exist elements \(a, b \in G\) such that \(a \cdot b \neq b \cdot a\), then there exist elements \(g_1, g_2 \in G\) whose products differ based on order: \(g_1 \cdot g_2 \neq g_2 \cdot g_1\).
For any \(f \in \mathbb {Z}_2^C\) and \(g \in \mathbb {Z}_2^E\),
That is, \(\langle \partial _2 f, g \rangle = \langle f, \delta _2 g \rangle \), so \(\delta _2 = \partial _2^T\).
For any \(f \in \mathbb {Z}_2^E\) and \(g \in \mathbb {Z}_2^V\),
That is, \(\langle \partial f, g \rangle = \langle f, \delta g \rangle \), so \(\delta = \partial ^T\).
For a valid cycle \(p\) and any vertex \(v\), the number of cycle edges incident to \(v\) is even:
(If \(v \notin p\): \(0\) edges incident; if \(v \in p\): exactly \(2\) edges incident.)
The number of edges in a cycle incident to any vertex is even (0 or 2). For a cycle \(p\) and vertex \(v\), if
then
If \(P\) is a deformable operator with no \(Z\)-support on edges (\(S_Z^E(P) = \emptyset \)) and \(\gamma \) is a valid deforming path for \(S_Z^V(P)\), then
That is, \(\tilde{P}\) commutes with all Gauss law operators.
If \(\gamma \) is a valid deforming path for \(S_Z^V\) (so that \(\partial \gamma = \mathbf{s}_{S_Z^V}\)), then the deformed operator commutes with the Gauss law operator \(A_v\) for every vertex \(v\):
Let \(s\) be a stabilizer check with no Z-support on edges (\(s.\mathrm{zSupportOnE} = \emptyset \)), and let \(\gamma \) be a valid deforming path with \(\partial \gamma = \mathcal{S}_{Z}(s) \cap V_G\). Then the boundary condition ensures anticommutations cancel:
Given a stabilizer check \(s\) with \(s.\mathrm{zSupportOnE} = \emptyset \) and a valid deforming path \(\gamma \), the deformed check commutes with all Gauss law operators:
Given a stabilizer check \(s\) with \(s.\mathrm{zSupportOnE} = \emptyset \) and a valid deforming path \(\gamma \) (i.e., \(\partial \gamma = \mathcal{S}_Z \cap V_G\)), the deformed check commutes with the Gauss law operator \(A_v\) for every vertex \(v\):
Composing two deformations with paths \(\gamma _1\) and \(\gamma _2\) is equivalent to a single deformation with their symmetric difference:
The deformed check’s edge \(Z\)-support is the symmetric difference of the original support and the path’s edge \(Z\)-operator support:
The deformed check \(\widetilde{s}_j\) has eigenvalue \(+1\) on the initial state \(|\psi \rangle |0\rangle ^{\otimes E}\):
For a stabilizer check \(s\) in Set \(\mathcal{S}\) with a non-empty valid path \(\gamma \) and \(s.\mathrm{zSupportOnE} = \emptyset \), the deformed check’s edge support differs from the original:
Given a stabilizer check \(s\) with \(s.\mathrm{zSupportOnE} = \emptyset \) and a valid deforming path \(\gamma \) with \(\partial \gamma = \mathcal{S}_Z \cap V_G\), the following hold:
The deformed check commutes with all Gauss law operators: for all \(v \in V\),
\[ \mathrm{deformed\_ gaussLaw\_ symplectic\_ simple}(G, s.\mathrm{zSupportOnV}, \gamma , v) \bmod 2 = 0. \]The deformed check preserves deformability:
\[ |(\mathrm{DeformedCheck}(G, s, \gamma )).\mathrm{zSupportOnV}| \bmod 2 = 0. \]
Given a valid deformed code setup, the following all hold:
Part 1: All \(A_v\) become stabilizers via measurement:
\(\forall v \in V,\; \mathrm{correctedOutcomes}(G, m)(v) = 0\),
\(\forall v, w \in V,\; \omega (A_v, A_w) = 0\).
Part 2: All \(B_p\) are stabilizers:
\(\forall p \in C,\; \sum _{e \in \mathrm{cycles}(p)} \mathrm{initialEdgeStabilizerOutcome}(e) = 0\),
\(\forall p \in C,\; \forall v \in V,\; \omega (B_p, A_v) \bmod 2 = 0\).
Part 3: All \(\widetilde{s}_j\) are stabilizers:
\(\forall j,\; \forall v \in V,\; \omega _{\mathrm{simple}}(\widetilde{s}_j, A_v) \bmod 2 = 0\),
\(\forall j,\; \mathrm{original\_ eigenvalue} + \sum _{e \in \gamma _j} \mathrm{initialEdgeStabilizerOutcome}(e) = 0\).
Part 4: The dimension is \(2^{k-1}\):
\[ \mathrm{total\_ qubits'} - \mathrm{total\_ stabilizers'} = k - 1. \]
The deformed code is an \([[n + |E|, k - 1, d']]\) code. Specifically, given \(k \geq 1\):
New qubit count: \(\mathrm{total\_ qubits'} = n + |E|\).
New logical qubit count: \(k - 1\).
Dimension formula: \(2^{(\mathrm{total\_ qubits'} - \mathrm{total\_ stabilizers'})} = 2^{k-1}\).
For a deformable operator \(P\), edge-path \(\gamma \), and edge \(e\),
That is, the binary vector representation of the deformed operator’s edge \(Z\)-support is the \(\mathbb {Z}_2\)-sum of \(P\)’s edge \(Z\)-support vector and \(\gamma \)’s characteristic vector.
Flux operators commute with Gauss law operators. Given the cycle validity hypothesis that each vertex has an even number of incident edges from the cycle:
For a valid cycle \(p\), \(B_p\) commutes with all \(A_v\). The symplectic form \(\omega (B_p, A_v)\) counts Z-X anticommutations, which equals the number of cycle edges incident to \(v\). This is always even:
In the computational basis, \(\lvert 0\rangle \) is a \(+1\) eigenstate of \(Z\). The eigenvalue of \(B_p\) on \(\lvert 0\rangle ^{\otimes E}\) is \((+1)^{|p|} = +1\). In \(\mathbb {Z}/2\mathbb {Z}\), the phase contribution from \(Z\) operators on \(\lvert 0\rangle \) states is \(0\):
\(B_p\) is Hermitian with eigenvalues \(\pm 1\). This is represented by \(B_p^2 = I\), which in \(\mathbb {Z}/2\mathbb {Z}\) is: \(2 \cdot \operatorname {support} = 0\). Formally, for all edges \(e\):
The initial eigenvalue of \(B_p\) is \(+1\): since \(Z|0\rangle = |0\rangle \), the eigenvalue is \(+1\) on the initial state \(|0\rangle ^{\otimes E}\). In \(\mathbb {Z}/2\mathbb {Z}\):
The flux operator support is exactly the characteristic vector of the cycle. For all edges \(e\):
Due to the constraint \(\prod _v A_v = L\), the number of independent \(A_v\) is \(|V| - 1\). However, \(L\) itself now becomes a stabilizer (it is the product of all \(A_v\) outcomes). Given \(|V| \geq 1\):
\(\mathrm{gaussLaw\_ independent\_ count}(G) = |V| - 1\),
\(\mathrm{gaussLaw\_ product\_ vertexSupport}(G) = v \mapsto 1\).
For any vertex \(v\) and coordinate \(w\),
This represents \(A_v^2 = I\). If \(A_v|\psi \rangle = \lambda |\psi \rangle \) and \(A_v^2 = I\), then \(\lambda ^2 = 1\), so \(\lambda \in \{ -1, +1\} \).
There exists a generator that is a linear combination of the others. Specifically, there exists \(v_0 \in V\) such that for all \(w \in V\),
The XOR of all measurement outcomes determines the \(L\) eigenvalue:
That is, the product is \(0\) (meaning \(L\) eigenvalue \(+1\)) if and only if the number of \(-1\) outcomes is even.
The Gauss law operators \(A_v\) form an obstruction: if \(P\) anticommutes with \(L\), then for any edge-path \(\gamma \), either there exists a vertex \(v\) at which the deformed Gauss law symplectic condition fails (i.e., the symplectic inner product is not even), or \(\gamma \) is not a valid deforming path:
The product constraint \(\prod _v A_v = L\) holds:
Vertex support: all \(1\)s (i.e., \(\mathrm{gaussLaw\_ product\_ vertexSupport}(G) = v \mapsto 1\), which equals the support of \(L = \prod _v X_v\)).
Edge support: all \(0\)s (each edge appears twice, \(X_e^2 = I\)), i.e., \(\mathrm{gaussLaw\_ product\_ edgeSupport}(G) = 0\).
The product of all Gauss law operators equals \(L\):
Specifically:
Vertex support: \(\sum _v \operatorname {vertexSupport}(A_v) = \mathbf{1}\) (the support of \(L = \prod _v X_v\)).
Edge support: \(\sum _v \operatorname {edgeSupport}(A_v) = \mathbf{0}\) (since \(X_e^2 = I\), edge terms cancel).
The Gauss law operator \(A_v\) corresponds to the coboundary of its vertex support:
This shows \(A_v = X_v \cdot \prod _{e \ni v} X_e\) in the binary vector representation.
If any valid deforming path exists for \(S_Z^V\), then a minimum-weight valid deforming path exists:
Complete formalization of Remark 12: There is no deformed version of a Pauli operator \(P\) that does not commute with the logical \(L\). Specifically:
If \(P\) anticommutes with \(L\), then \(|\mathcal{S}_Z \cap V_G|\) is odd.
A path boundary \(\partial \gamma \) always has even cardinality.
Therefore no edge-path \(\gamma \) with \(\partial \gamma = \mathcal{S}_Z \cap V_G\) exists.
Hence \(P\) has no well-defined deformed version.
If a Pauli operator \(P\) anticommutes with \(L\) (i.e., \(|\mathcal{S}_Z \cap V_G|\) is odd), then no valid deforming path exists:
This formalizes the remark’s central claim: there is no deformed version of a Pauli operator that does not commute with the logical \(L\).
Given a graph \(G\) and a gauging measurement \(m\), the following three properties hold:
After measurement and corrections, all \(A_v\) have eigenvalue \(+1\): \(\forall v \in V,\; \mathrm{correctedOutcomes}(G, m)(v) = 0\).
All \(A_v\) mutually commute: \(\forall v, w \in V,\; \omega (A_v, A_w) = 0\).
Each \(A_v\) is self-inverse: \(\forall v, w \in V,\; 2 \cdot \mathrm{gaussLawOperator\_ vertexSupport}(G, v, w) = 0\).
Given that \(p\) is a valid cycle (each vertex has \(0\) or \(2\) incident edges from \(p\)), the following hold:
\(B_p\) has initial eigenvalue \(+1\) on \(|0\rangle ^{\otimes E}\).
\(B_p\) commutes with all measured operators \(A_v\).
Therefore \(B_p\) remains a stabilizer after the gauging measurement.
Given:
\(s_j\) is a stabilizer from the original code (no Z-support on edges),
\(\gamma \) is a valid deforming path with \(\partial \gamma = \mathcal{S}_{Z,j} \cap V_G\),
the original code state \(|\psi \rangle \) is a \(+1\) eigenstate of \(s_j\) (input assumption),
then:
\(\widetilde{s}_j\) commutes with all \(A_v\) (by the boundary condition).
\(\widetilde{s}_j\) has eigenvalue \(+1\) on \(|\psi \rangle |0\rangle ^{\otimes E}\).
Therefore \(\widetilde{s}_j\) is a stabilizer of the deformed code.
Every stabilizer check falls into exactly one of two categories:
\(s\) is in Set \(\mathcal{C}\) and \(\mathrm{DeformedCheck}(G, s, \emptyset ) = s\), or
\(s\) is in Set \(\mathcal{S}\) and \(\neg \, \mathrm{IsValidDeformingPath}(G, s.\mathrm{zSupportOnV}, \emptyset )\).
The logical structure of the argument consists of four parts, all holding simultaneously:
\(\mathrm{anticommutesWithL}(\mathcal{S}_Z) \implies |\mathcal{S}_Z| \bmod 2 = 1\).
\(\forall \gamma ,\; |\mathrm{boundarySupport}(G,\gamma )| \bmod 2 = 0\).
\(\forall \gamma ,\; \mathrm{IsValidDeformingPath}(G,\mathcal{S}_Z,\gamma ) \implies \mathrm{boundarySupport}(G,\gamma ) = \mathcal{S}_Z\).
\(\mathrm{anticommutesWithL}(\mathcal{S}_Z) \implies \nexists \gamma ,\; \mathrm{IsValidDeformingPath}(G,\mathcal{S}_Z,\gamma )\).
For a stabilizer check \(s\) in Set \(\mathcal{C}\) with \(s.\mathrm{zSupportOnE} = \emptyset \), and for any vertex \(v\):
If \(P\) anticommutes with \(L\), then for every edge-path \(\gamma \), \(\gamma \) is not a valid deforming path for \(\mathcal{S}_Z\). This captures the paper’s claim: “there is no way to multiply such a \(P\) with stabilizers \(Z_e\) and \(s_j\) to make it commute with all the Gauss’s law operators \(A_v\).”
For any monomial \(\alpha \), the \(\bar{Z}'_\alpha \) operator has zero left polynomial and the \(\bar{X}_\alpha \) operator has zero right polynomial:
For any monomial \(\alpha \), we have:
For any monomial \(\beta \), we have:
A fault has empty syndrome if and only if all detectors in the collection are satisfied:
For any two vertices \(v_1, v_2 \in V\), the generalized Gauss law operators \(A_{v_1}\) and \(A_{v_2}\) commute. Specifically, the symplectic inner product of their supports is zero on both vertex and hyperedge qubits:
If \(t_{\mathrm{initial}} - t_{\mathrm{start}} = \texttt{roundsBefore}\) for an error correction configuration \(\textit{ec}\), then the initial boundary separation equals \(\texttt{roundsBefore}\):
A fault is a logical fault if and only if it has empty syndrome and is non-trivial:
A fault \(f\) with empty syndrome is a spacetime logical fault if and only if it is not a spacetime stabilizer:
When the total vertex count exceeds the logical support cardinality, there exist vertices outside the support (i.e., unconstrained by relaxed expansion):
The sum of all Gauss’s law vertex supports on the ladder is the all-ones vector:
This is an instantiation of the general property from the Gauss’s law operator definition.
Combining: \(\prod _v A_v = L\) on the ladder, i.e., the vertex support equals the all-ones vector and the edge support equals zero:
Relaxed expansion is strictly weaker than full expansion: subsets \(S\) disjoint from the logical support are unconstrained. If \(\mathrm{Disjoint}(S, \mathrm{logicalSupport})\), then \((S \cap \mathrm{logicalSupport})\) is not nonempty.
When the logical support equals the full vertex set \(V\), relaxed expansion implies full expansion, since every nonempty subset intersects the full vertex set:
After the split step, gauge qubits are discarded and only original code qubits remain. For surface code blocks of size \(\mathrm{rows} \times \mathrm{cols}\) with \(\mathrm{rows}, \mathrm{cols} \geq 2\):
Let \(\textit{ec}\) be an error correction configuration providing full protection, and let \(\textit{ep}\) be an error process with \(\texttt{weight} \leq d\). Under the propagation model, \(\textit{ep}\) does not span to any boundary:
A measurement at time \(t + \tfrac {1}{2}\) occurs between qubit states at times \(t\) and \(t + 1\). Formally, for any measurement time step \(t\):
A detector \(D \in \mathrm{DC}.\mathrm{detectors}\) is in the syndrome if and only if it is violated by the fault:
Let \(c {\lt} d\). If a mid-computation implementation uses \(c\) rounds before gauging and an isolated implementation uses \(d\) rounds before gauging (with the same code distance), then the mid-computation implementation uses strictly fewer rounds before gauging.
For graphs satisfying the Freedman–Hastings bound, there exists a valid layer assignment with \(R = O(\log ^2 n)\) where \(n = |V|\). Specifically, there exists \(R \in \mathbb {N}\) such that:
(Corollary 1: Overhead Bound.) For any LDPC code with parameters \(W \geq 2\), \(d \geq 1\), and Freedman–Hastings constant \(C_{\mathrm{FH}} {\gt} 0\), there exists \(c {\gt} 0\) and a qubit overhead structure \(Q\) with \(Q.W = W\) such that
If \(R \le C \cdot (\log _2 W)^2 + C\), then
This is the proven arithmetic consequence of the Freedman–Hastings specification, showing \(O(W \log ^2 W)\) overhead.
Gauging measurement can be applied to multiple logical operators in parallel if and only if both conditions are satisfied. Formally, for a stabilizer code \(C\), a finite set \(\mathcal{L}\) of logical operators, and a bound \(k\):
If a Pauli operator \(P\) commutes with an \(X\)-type logical operator \(L\) whose support is the full set \(\mathrm{Fin}(n)\), then the Z-type support of \(P\) has even cardinality:
If a Pauli operator \(P\) commutes with an \(X\)-type logical operator \(L = \prod _{v \in \mathrm{supp}(L)} X_v\) of a stabilizer code \(C\), then the Z-type support of \(P\) restricted to \(\mathrm{supp}(L)\) has even cardinality:
For any overhead configuration \(\mathrm{cfg}\) with parameters \(W \geq 2\), \(d \geq 1\), and \(C_{\mathrm{FH}} {\gt} 0\),
where \(c(d, C_{\mathrm{FH}}) = 2d + C_{\mathrm{FH}}(d+1) + C_{\mathrm{FH}}\).
The CX commutation relation \(\mathrm{CX}^\dagger (X \otimes I) \mathrm{CX} = X \otimes X\) gives the Gauss’s law operator. Formally, for each \(i \in \mathrm{Fin}(W)\), both \((\mathtt{false}, i) \sim (\mathtt{true}, i)\) and \((\mathtt{true}, i) \sim (\mathtt{false}, i)\) hold in \(\mathrm{shorPathGraph}(W)\).
The \(Z\) stabilizers of the GHZ state on \(W\) qubits are \(Z_i Z_{i+1}\) for \(i = 0, \ldots , W-2\). These correspond exactly to the \(Z\) operators on the dummy path edges: for each \(i \in \mathrm{Fin}(W-1)\), there exists an edge \(e\) that is a dummy edge with endpoints \((\mathtt{true}, i)\) and \((\mathtt{true}, i+1)\).
The full Shor-style/gauging correspondence consists of five properties:
There are \(W\) cross edges,
There are \(W - 1\) dummy path edges,
Each cross edge \(\mathrm{cross}(i)\) has endpoints \((\mathtt{false}, i)\) and \((\mathtt{true}, i)\),
Each dummy path edge \(\mathrm{dummyPath}(i)\) has endpoints \((\mathtt{true}, i)\) and \((\mathtt{true}, i+1)\),
The total edge count is \(2W - 1\).
Combined result: the product of all Gauss’s law operators on the Shor path graph yields \(L\). That is, the vertex support product is the all-ones vector and the edge support product is zero:
For a small code instance with very high-weight flux checks (\(\texttt{maxCheckWeight} {\gt} 2d\)), the trade-off is favorable: the code instance remains small (the benefit of avoiding high-weight measurements outweighs the lack of threshold for small instances).
A space fault at time \(t\) (with \(t {\gt} 0\)) is between measurements at \(t - \tfrac {1}{2}\) and \(t + \tfrac {1}{2}\):
For any detector configuration and time \(t\), and any check index \(j {\lt} n_{\mathrm{checks}}\), the elementary detector \(\langle \mathrm{originalCheck}(j), t, \mathrm{bulk} \rangle \) belongs to the set of bulk original check detectors.
The elementary detector parities are all zero in the error-free case. Specifically, the following five conditions hold simultaneously:
Bulk detectors: For all \(m\), \(\mathrm{bulkDetectorParity}(m, m) = 0\).
Initial \(B_p\): \(\mathrm{xorParity}(0, 0) = 0\).
Initial \(\tilde{s}_j\): For all \(m_{s_j}\), with \(m_{\tilde{s}} = m_{s_j} + 0\), \(\mathrm{xorParity}(m_{s_j}, m_{\tilde{s}}) = 0\).
Final \(B_p\): For all \(m_{B_p}, m_{Z_e}\), if \(m_{B_p} = m_{Z_e}\) then \(\mathrm{xorParity}(m_{B_p}, m_{Z_e}) = 0\).
Final \(\tilde{s}_j\): For all \(m_{\tilde{s}}, m_{s_j}, m_{Z_\gamma }\), if \(m_{\tilde{s}} = m_{s_j} + m_{Z_\gamma }\) then \(m_{\tilde{s}} + m_{s_j} + m_{Z_\gamma } = 0\).
Since \(B_p = \prod _{e \in p} Z_e\) by definition, measuring \(B_p\) and measuring all \(Z_e\) individually then taking the product (XOR in \(\mathbb {Z}/2\mathbb {Z}\)) give the same result. If \(m_{B_p} = m_{Z_e\text{-product}}\), then \(\mathrm{xorParity}(m_{B_p}, m_{Z_e\text{-product}}) = 0\).
For any gauging region \(R\) and any time \(t\), exactly one of the following holds:
That is, the five time regions partition all times.
The time regions are mutually exclusive. Specifically, for any time \(t\):
\(\neg (\mathrm{isBefore}(t) \land \mathrm{isStart}(t))\),
\(\neg (\mathrm{isBefore}(t) \land \mathrm{isDuring}(t))\),
\(\neg (\mathrm{isStart}(t) \land \mathrm{isDuring}(t))\),
\(\neg (\mathrm{isDuring}(t) \land \mathrm{isEnd}(t))\),
\(\neg (\mathrm{isDuring}(t) \land \mathrm{isAfter}(t))\).
At \(t = t_i\), edge qubits are initialized in \(|0\rangle \). Since \(|0\rangle \) is a \(+1\) eigenstate of \(Z\), we have \(B_p = \prod _{e \in p} Z_e\) with eigenvalue \(\prod _{e \in p}(+1) = +1\). The detector at \(t_i\) compares the initialization outcome (\(+1\), encoded as \(0\)) with the \(B_p\) measurement (\(+1\), encoded as \(0\)). The parity is \(\mathrm{xorParity}(0, 0) = 0\).
At \(t_i - \tfrac {1}{2}\), we measure \(s_j\) with outcome \(m_{s_j}\). At \(t_i + \tfrac {1}{2}\), we measure \(\tilde{s}_j = s_j \cdot Z_\gamma \). Since \(Z_\gamma |0\rangle = |0\rangle \) (eigenvalue \(+1\), encoded as \(0\)), we have \(m_{\tilde{s}} = m_{s_j} + 0 = m_{s_j}\). Therefore \(\mathrm{xorParity}(m_{s_j}, m_{\tilde{s}}) = 0\).
The XOR of consecutive bulk detector parities telescopes to the endpoints. For a sequence of outcomes \(a_0, a_1, \ldots , a_{n+1}\):
This is because in \(\mathbb {Z}/2\mathbb {Z}\), \((a_0 + a_1) + (a_1 + a_2) + \cdots + (a_n + a_{n+1}) = a_0 + a_{n+1}\), since middle terms appear twice and cancel.
The elementary detectors form a generating set for the fault-tolerant gauging measurement procedure. This is established through eight conditions:
Part 1 — Verification:
For all \(m\), \(\mathrm{bulkDetectorParity}(m, m) = 0\).
\(\mathrm{xorParity}(0, 0) = 0\).
For all \(m_{s_j}\), \(\mathrm{xorParity}(m_{s_j}, m_{s_j} + 0) = 0\).
For all \(m_{B_p}, m_{Z_e}\) with \(m_{B_p} = m_{Z_e}\), \(\mathrm{xorParity}(m_{B_p}, m_{Z_e}) = 0\).
For all \(m_{\tilde{s}}, m_{s_j}, m_{Z_\gamma }\) with \(m_{\tilde{s}} = m_{s_j} + m_{Z_\gamma }\), \(m_{\tilde{s}} + m_{s_j} + m_{Z_\gamma } = 0\).
Part 2 — Completeness:
For all \(t {\lt} t_i\) and all check indices \(j\), there exists a bulk detector for \(s_j\) at \(t\).
For all cycle indices \(p\) and check indices \(j\), initial boundary detectors exist.
For all cycle indices \(p\) and check indices \(j\), final boundary detectors exist.
Let the logical effect be group-like. If \(F_{\mathrm{space}} \cdot F_{\mathrm{time}}\) affects the logical information, \(F_{\mathrm{space}}\) has empty syndrome, and \(F_{\mathrm{time}}\) is a spacetime stabilizer, then \(F_{\mathrm{space}}\) is a spacetime logical fault.
Let \(\mathrm{DC}\) be a detector collection, let the logical effect be group-like, and let the syndrome be a group homomorphism. Let \(I = [t_i, t_o]\) be a gauging interval, and let \(F\) be a spacetime logical fault.
Suppose that for any such fault, there exists a cleaning stabilizer \(S_{\mathrm{clean}}\) (built from Pauli pair stabilizers as in Lemma 4) such that \(S_{\mathrm{clean}}\) is a spacetime stabilizer and the space errors of \(F \cdot S_{\mathrm{clean}}\) are concentrated at time \(t_i\).
Then there exist spacetime faults \(F_{\mathrm{space}}\) and \(F_{\mathrm{time}}\) such that:
\(F\) is equivalent to \(F_{\mathrm{space}} \cdot F_{\mathrm{time}}\) modulo stabilizers.
\(F_{\mathrm{space}}\) is a pure space fault at the single time step \(t_i\).
\(F_{\mathrm{time}}\) is a pure time fault (only measurement errors).
Let the logical effect be group-like. If \(F_{\mathrm{space}} \cdot F_{\mathrm{time}}\) affects the logical information, \(F_{\mathrm{time}}\) has empty syndrome, and \(F_{\mathrm{space}}\) is a spacetime stabilizer, then \(F_{\mathrm{time}}\) is a spacetime logical fault.
We have \(d_{\mathrm{ST}} {\gt} 0\) if and only if logical faults exist and every spacetime logical fault has positive weight:
If \(F\) is a spacetime logical fault, \(F\) decomposes as \(F = F_{\mathrm{space}} \cdot F_{\mathrm{time}} \cdot S\) for some stabilizer \(S\), \(F_{\mathrm{time}}\) is a nontrivial spacetime logical fault that spans the interval \([t_i, t_o)\), and \(|F| \geq |F_{\mathrm{time}}|\), then \(|F| \geq d\).
If \(F\) is a spacetime logical fault equivalent to a pure space fault \(F_{\mathrm{space}}\) at time \(t_i\) (modulo a stabilizer \(S\)), \(|F_{\mathrm{space}}| \geq d\) (from the space distance bound with \(h(G) \geq 1\)), and \(|F| \geq |F_{\mathrm{space}}|\) (cleaning preserves weight), then \(|F| \geq d\).
Under the conditions:
\(h(G) \geq 1\) (strong expansion),
\((t_o - t_i) \geq d\) (sufficient measurement rounds),
the spacetime fault-distance equals exactly \(d\):
There exists a spacetime logical fault of weight exactly \(d\). Specifically, given an original code logical operator \(L_{\mathrm{orig}}\) of weight \(d\) applied at a time outside the deformation region (either \(t {\lt} t_i\) or \(t {\gt} t_o\)) that produces a spacetime logical fault, the resulting spacetime fault has weight \(d\).
For any fault pattern \(p\) satisfying \(\mathrm{hasEmptySyndrome}(p)\) and \(\mathrm{preservesLogical}(p)\), there exists a list of spacetime stabilizer generators \(\mathrm{gens}\) such that every generator in the list has empty syndrome and preserves logical information:
For any qubit \(q\), Pauli type \(P\), time step \(t\), positive integer \(k\), and function assigning anticommuting checks, there exist \(k\) spacetime stabilizer generators (each a pauliPairOriginal). That is:
Each factor \((P_{t+i}, P_{t+i+1})\) is a pauliPairOriginal generator at adjacent time steps.
The net Pauli effect of \(k\) adjacent pairs is \(P\) at the first time and \(P\) at the last time: each intermediate \(P\) appears twice and cancels (\(P \cdot P = I\)). Formally, \(\mathrm{foldl}(\lambda \, \mathrm{acc}\, \_ .\; \mathrm{acc} + 2,\; 0,\; [0, \ldots , k{-}1]) = 2k\).
\(\mathrm{netEffect}_{c^{t+1}}(\text{pauliPairEffect\_ anticommuting}) = 0\).
The recorded outcome at \(t + \tfrac {1}{2}\) equals base (as above), and the physical state at \(t + \tfrac {3}{2}\) is \(\text{base} + 1 + 1 = \text{base}\) (since \(P\) at \(t\) and \(P\) at \(t+1\) give \(P^2 = I\)).
\(\mathrm{netEffect}_{c^t}(\text{pauliPairEffect\_ anticommuting}) = 0\).
The measurement at \(t - \tfrac {1}{2}\) is at the base value (unaffected), and the recorded outcome at \(t + \tfrac {1}{2}\) is \(\text{base} + 1 + 1 = \text{base}\) (the Pauli flip and measurement fault cancel).
If every generator \(g\) in a list \(\mathrm{gens}\) satisfies \(\mathrm{logicalEffect}(g) = 0\), then the \(\mathbb {Z}_2\)-sum of their logical effects is \(0\):
For any fault pattern satisfying the local spacetime stabilizer conditions (empty syndrome and preserving logical information), there exists a decomposition into the listed generators. This shows the generators form a generating set spanning all local stabilizers.
The listed generators form a generating set of local spacetime stabilizers for the fault-tolerant gauging measurement procedure:
Empty Syndrome: For every generator \(g\), \(\mathrm{generatorHasEmptySyndrome}(g)\) holds.
Preserves Logical: For every generator \(g\), \(\mathrm{logicalEffect}(g) = 0\).
No fault can be both a spacetime stabilizer and a spacetime logical fault:
A fault cannot be simultaneously a spacetime stabilizer and a logical fault:
The union of spacetime stabilizers and spacetime logical faults equals the set of empty-syndrome faults:
For a graph \(G\) on vertices \(V\) and any two adjacent vertices \(u, v \in V\) with \(G.\mathrm{Adj}(u,v)\), there exists a walk \(w\) from \(u\) to \(v\) of length exactly \(1\):
This captures the Step 1 property: matched pairs have a direct edge, giving path length 1.
The subdivided graph is bipartite: every edge connects an original vertex to a dummy vertex. Formally, for all \(x, y \in \texttt{SubdividedVertex'}(V, \mathrm{Sym2}(V))\) with \(\mathrm{subdivisionAdj}'(G)(x,y)\), either
or
For each original edge \(e = \{ u, v\} \) in \(G\) (i.e., \(G.\mathrm{Adj}(u, v)\)), the subdivided graph contains both adjacencies:
Under the syndrome linearity condition, the syndrome is \(\mathbb {Z}_2\)-linear with respect to fault composition:
Under the detector structure assumptions (that detectors correspond to interior comparison detectors), three results hold simultaneously:
The \(A_v\) chain is a pure time logical fault.
The \(A_v\) chain has weight exactly \(I.\mathrm{numRounds} = t_o - t_i\).
For all pure time logical faults \(F\): if no interior comparison detectors are violated and faults exist in the interval, then \(I.\mathrm{numRounds} \leq F.\mathrm{weight}(\mathrm{times})\).
Therefore the pure time fault distance is exactly \(t_o - t_i\).
Any pure time spacetime logical fault \(F\) satisfies:
provided the interval rounds are contained in \(\mathrm{times}\), no interior comparison detectors are violated, and faults exist in the interval.
Every Type 2 fault string \(s\) has a decomposition \(d\) such that all propagator generators are of the form \(\mathrm{propagator}(s.\mathrm{edge}, t)\), the initial generator is \(\mathrm{initial}(s.\mathrm{edge})\), and the final generator is \(\mathrm{final}(s.\mathrm{edge})\).
Let \(F\) be a pure time fault, \(I\) a deformation interval, and \(\mathrm{times}\) a set of time steps containing all interval rounds. Assume no interior comparison detector is violated for any check \(m\), and there exists some check \(m\) and time \(t_0 \in [t_i, t_o)\) with odd fault count. Then:
If the code can tolerate weight-\(t\) faults and \(|F|_{\mathit{times}} \leq t\), then either \(F\) is detectable or \(F\) is a stabilizer:
The deformed code distance \(d^*\) is at most \(d\) when the original code has a logical operator that can be trivially extended. Specifically, if \(\ell \) is a trivially extendable logical, then
For any vertex set \(S \subseteq V\) and edge \(e \in E\),
where \(\delta \) is the coboundary map and \(\mathbf{1}_S\) is the characteristic function of \(S\) over \(\mathbb {Z}/2\mathbb {Z}\).
Given:
\(h_G \ge 0\), \(d {\gt} 0\),
\(\mathrm{cleanedWeight} \ge d\),
\(\mathrm{cleaningSetSize} \le \mathrm{cleanedWeight}\),
\(\mathrm{boundarySize} \ge h_G \cdot \mathrm{cleaningSetSize}\),
then
For an \(X\)-type logical operator \(L\) over a stabilizer code \(C\), the Pauli content of \(L\) on each qubit matches \(\operatorname {pureX}\) applied to its support. That is,
This is the formal statement that \(L = \prod _{v \in \operatorname {supp}(L)} X_v\) (ignoring global phase).