Magic State Factory Catalog

[\![210,\,40,\,3]\!] factory

\mathrm{T}_{0}\cdot \mathrm{T}_{1}\cdot \mathrm{T}_{2}\cdot \mathrm{T}_{3}\cdot \mathrm{T}_{4}\cdot \mathrm{T}_{5}\cdot \mathrm{T}_{6}\cdot \mathrm{T}_{7}\cdot \mathrm{T}_{8}\cdot \mathrm{T}_{9}\cdot \mathrm{T}_{10}\cdot \mathrm{T}_{11}\cdot \mathrm{T}_{12}\cdot \mathrm{T}_{13}\cdot \mathrm{T}_{14}\cdot \mathrm{T}_{15}\cdot \mathrm{T}_{16}\cdot \mathrm{T}_{17}\cdot \mathrm{T}_{18}\cdot \mathrm{T}_{19}\cdot \mathrm{T}_{20}\cdot \mathrm{T}_{21}\cdot \mathrm{T}_{22}\cdot \mathrm{T}_{23}\cdot \mathrm{T}_{}\cdots

distance, proved here
d=3
error coefficient
A_{3}=228
\gamma_\rho
1.5094
V_{\mathrm{ex}}
40
wires
N=53
Clifford correction
none with T powers
found by
AI search pure T

The circuit

The 53\times210 matrix. Rows are wires and columns are rotations. A 1 means the wire takes part in the rotation, and the number beside each wire is its row weight. Scroll sideways to see every column.

  • output wires
  • check wires, measured and postselected
wire
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200
out 0 (61)··1·1····1·1···1·1·1···1···1·1·····11····1111···········11····1·1···1··1·1··1···1·1········1·1·····1111··········1···1·1···1···1·1··1···1·1······1··1111········1·1····11···1···11···1····11·····1·1····111·······
out 1 (45)···············1·······1·1·1·1·····11····1111······················11·······1·1·1·1········1·1·····1111··················1·1········11··1·1······1··1111··················1·1········1·1··11·····1·1····111·······
out 2 (89)·1·1···11··1····11111····1·1·1·1·1·1·1···11·1···1·······111··11···1···11·1·1·111······1···111···1·11····1···111·1·111··1····1··1111··1··1·1·1·····1·11·1··1·····11111·11·1·1···111··1·111····1··111··1·1····1··1·1
out 3 (63)····1···1··1·1·1···11····1···1······1·····111··1·1··1····1···1··1·1·1····11···1···1··········1······111··1·1·1···1··1··1···1·····11··1····1······1···111·1·1·1····1···1·11··1····11····1···1·······1····111··11··1
out 4 (81)·1·11··1·1·1···11·1·1··········1·1·1·1·1·1·1·1··1···1···1·········1·1·1····1·1·11·1···1···111·····1111111···1·1·11·1·1·1···11···1·1·········1·····1·1·1·1·1··1··11··1····1·11··1····1·1·1·11·1··111····111111··1··
out 5 (47)····1···1········1··1······1········1·1···1·11·1····1····1···1·········1··1·····1············1·1····1·11·1···1···1··1··········1··1·····1········1···1·111···1····1···1·········1·1·······1········11···1·11·1···1
out 6 (31)····1····1·····1·······1···1·········1·····1·1···········1····1·····1·······1···1·············1······1·1·········1···1·····1········1···1·········1···1·1·········1····1····1········1····1··············1·1······
out 7 (45)····1····1·····1·······1···1·······11····11·11···········1····1·····1·······1···1··········111·····11·11·········1···1·····1········1···1·······11··11·11·········1····1····1········1····1······111····1·11······
out 8 (105)·1·1···111·1···11·111··1·1·1···1·1·111···1·1···111··1···11···1····1·1····1·11111··1···1···1111··1·11111·11·11·1·1··111·1···11·1·111·11··1···1····11·1·1··111·1··111·1·1··1·11····1··11111··1·1··1111·1·1111·1111··
out 9 (99)·1·11··11···1···11111····1·····1·1··11······11·111··1···1····11·111···11·1·1·1111·1···1···1·11··1·1·1··111·11·1·11·11···1···1··1111··1······1····11····11111·1··11··1·1111·1···111··1·111·11·1··1·11·1·11··11111··
out 10 (83)·1·1···1···1···111111·····11·1·1·1···1····1·11··1·······11····1···1·1·11·1·1·1·1······1···1·1···1·1····11···111·1··1···1···11··1111····11·1·1·····1··1·11·1·····111·1··1·1·11··111··1·1·11···1··1·1··1·1···11··1·1
out 11 (91)·1·1···111·1···111111··11····1·1·1···1····11·1··1·······11···1····1·1·11·1·11··11·····1···1·1···1·1··1111···111·1··111·1···11··1111·1·····1·1·····1··11·1·1·····111·1·1··1·11··111··11··1·1··1··1·1··1·1·1111··1·1
out 12 (97)·1·1···1···1···111111······111·1·1·111···1111··11·······11····1···1·1·11·1·1·1·1·1····1···1111··1·11·1··11··111·1··1···1···11··1111·····111·1····11·1111·11·····111·1··1·1·11··111··1·1·1····1··1111·1·1·1··11·1·1
out 13 (31)···········1·····1··1······1········1·····1····1·······1········1······1··1·····1············1······1····1·····1·······1·······1··1·····1········1···1···1·····1········1·······1·1·······1········1····1····1····
out 14 (29)·····1··1·····1········1·····1······1·····1··1···············1··············1·····1··········1······1··1············1·····1·········1·····1······1···1··1···········1·1····1·········1·····1·······1····1··1······
out 15 (55)········1··1···1·1·11··1·1·1··············11·1···1···········1··1···1··1·11·1·1·1····1··············11·1···1········1··1···1···1·11·11··1····1·······11·1··1··········1·1···1···111··1·1··1·1···········11·1··1···
out 16 (93)·1··1··11······1111·1····1···1·1····11····1·11·111······1····11·1·1·1·11···1·1111·····1·1·1·11··1·1····111·1111··1·11······11··11·1··1····1·1····11··1·11111····1···1·1111·11··11···1·111·1··1··1·11·1·1···11111·1
out 17 (83)·1·1···111·1····1·11···1·····1·1·1···1····111···11······11···1····1···1··11111·11·····1···1·1···1····1··1··1111·1··111·1····1···11··1·····1·1·····11·111··11····111·1·1··1·1···1·11·111·1·1··1··1·1··1···1··1·11·1
out 18 (63)·1······11·1·····1··1··1·1···1······1·····111··1·1··1·1······11·1······1··1·1·1···1··········1······111··1·1·11·····11·1·······1··1·11····1······1···111·1·1·1········111·······1·1··1·1···1·······1····111··11··1
out 19 (63)····1···1······1·1·11····1·1········1·····1111···1··1····1···1······11·1·11···1·1············1······1111···1·1···1··1······1·1·1·11··1··1········1···1111··1·1····1···1·····1·1·111····1··1········1····1111··1··1
out 20 (97)·1·1···11······11111·····1·1···111·1·1···11·····11······11···11·1·1·1·11·111·111··1···11··111···1·11··1·1··1111·1··11······11··111···1··1···1·····1·11····11····111·1·1111·11··1111·1·111··1·11·111··1·1··1·1·11·1
out 21 (63)1···1···1··1···1·1·1·····1·1·1············11·1·1·1··11···1···1··1···1··1·1····1·1·1·················11·1·1·1·1···1··1··1···1···1·1···1··1·1··········11·11·1·11···1···1·1···1···11·····1··11············11·1·11··1
out 22 (83)·1·11··1·1·1····11111··1···1·1·1·1··11····1··11·1·······1·········1···11·1·111·1······1···1·11··1·1···11····111·11·1·1·1····1··1111·1···1·1·1····11··1··1·1·····11··1····1·1···111··111·1····1··1·11·1·1··11···1·1
out 23 (45)····1···1······1···11··1·····1·····1·····1·····1···11····1···1······1····11·1·····1········1·······1·····1···1···1··1······1·····11·1·····1·········1····1··11····1···1·····1····11··1·····1·····1···········1··11
out 24 (97)·1·1···111·1···111111····1·1···1·1·111···11····11·······11···1····1·1·11·1·1·11···1···1···1111··1·11··1·11··111·1··111·1···11··1111··11·1···1····11·11···11·····111·1·1··1·11··111··1·11···1·1··1111·1·1··1·11·1·1
out 25 (31)······1·1··············1·1·1·1············1·1··············1·1··············1·1·1·1·················1·1·············1···············11··1·1··········1·1·············11··············1·1··11············1·1·······
out 26 (89)·1·1···1··1····1111········1···1·1·111···11····11·······11····111·1·1·11··11·1·1··1···1···1111··1·11··1·11··111·1··1··1····11··11·······1···1····11·11···11·····111·1··111·11··11·1·1·1·1··1·1··1111·1·1··1·11·1·1
out 27 (31)···············1·1·1···1················1··1·1···1··················1··1·1··1····················1···1·1···1···············1···1·1··1·················1·1··1················1···11···1················1··1·1··1···
out 28 (61)·········1·1···1·1·····1·1···1·····11····1111··1·1············1·1···1··1····1·1···1···1····1·1·····1111··1·1·········1·1···1···1····11····1···1··1··1111·1·1···········11···1···1····1·1···1·1···1·1····111··11···
out 29 (47)········1··1········1····1·1·1······1·····1·1··1····1········1··1·········1···1·1·1·······1··1······1·1··1···1······1··1··········1··1··1·1····1·1···1·1·1···1········1·1·········1····1··11····1··1····1·1··1···1
out 30 (75)·1·11··1·1·1····1······1···1·1·1·1··11····111···1·······1·········1···1·1·1111·1······1···1·11··1·1··1··1···111·11·1·1·1····1·······1···1·1·1····11··111··1·····11··1····1·1···1··1·111·1····1··1·11·1·1·1··1··1·1
out 31 (45)····1···1······1····11···1···1·····1·····1···1···1··1····1···1······1·····1···1···1········1·······1···1···1·1···1··1······1······11·1····1·········1···1··1·1····1···1·····1·····11···1···1·····1·········1··1··1
out 32 (47)····1···1······1···1··1··1·1········1······1···1·1··1····1···1······1····1·1··1·1············1·······1···1·1·1···1··1······1·····1···1··1········1····1··1·1·1····1···1·····1····1··1··1··1········1·····1···11··1
out 33 (31)····1···1·······1······1···1········1·······11···········1···1··············1···1············1········11·········1··1·······1·······1···1········1·····11·········1···1······1·······1····1········1······11······
out 34 (45)·······1·1·············1·1·········11····11111·1············1·1·············1·1············1·1·····11111·1·········1·1··············11···········1··111111·············1·············1·1·········1·1····1111·1····
out 35 (83)·1·11··1········111·1··1·1·····1·11··1······1··11···1···1·····1·1·1···11···111111·1···1··11·1···1·1·1···11··1·1·11·1········1··11·1·11······1·····1····1·11··1··11··1··111·1···11···11111·11·1·11·1··1·11···11·1··
out 36 (47)····1······1·····1·······1·1·11·····1·····1··1·1····1····1······1······1······1·1·11·········1······1··1·1···1···1·····1·······1·····1··1·11·····1···1··11···1····1·····1·······1······1··11·······1····1··1·1···1
out 37 (47)····1···1······1···11··1···1········1······1·1···11······1···1······1····11·1···1············1·······1·1···11····1··1······1·····11·1···1········1····1·1··1······1···1·····1····11··1····1········1·····1·1··11··
out 38 (71)····1···11·1·····1·11····1·1·1·1····1·····1111·1····1····1···11·1······1·11···1·1·1·1········1······1111·1···1···1··11·1·······1·11··1··1·1·1····1···11111···1····1···111·······111····1··11·······1····1111·1···1
out 39 (89)·1·11··1·1·1····11111··1·······1·1·1·1···111·1··········1·········1···11·1·111·11·1···1···111···1·11·1111·1·111·11·1·1·1····1··1111·1·······1·····1·111·1·······11··1····1·1···111··111·1·11·1··111··1·1·1111··1·1
check 40 (122)·1·····1··1·1··11·11·1···11·1111·1·1·1······11··1··1·1·11111·11·1·11·111··1111·11····111·11·11·1111111··11·1111····1··1·1··11···11·1·1·1·1111·····1····11·1·1·11111111111111··111·1·111·1·1·11111·11111111··1111·1
check 41 (86)········11·11·1·111·11·1·1·11·······11·······1··1··11··11·1·····1·1·1·····111··11·1··11··11······11··111·1··········11·11·1·1··11·1111··11·······11·····1·1·11·11··11···11·11·····1·11··1·1111·11·····11·111·1····
check 42 (84)···11···111·1·11·11···········11···1·1····1·11·1·1·1······1··11·1·1·111·····1111·····11·····11·····1·1·1·1··11··11··111·1·11···11··········11·····1··1·111·11······11·1111·11·11·····1111···11····11·····1·1·1·1·1
check 43 (84)···1···1·1111··1···111·······1···1·1·1···1·11··111·11······1·1····1·1······11111··1····1·11·11·····1·1111···1···1··1·1111··1·····111······1·······1·1·11·11111···1···11··1··1·······11111··1··111·11·····1111··1··
check 44 (64)··············1·111·11···················1··1··111·1···············1111··1·11111···············1111·11·11···11············1·1··11·11················1··1·1111·············1·1·11·1··11111···········111111·11··1·1
check 45 (94)···1····11·1···1···1···1·11·1111·1··11···11·11··11·11········11·1···1····1··1·1·1····111·1111··1111··1··11··1···1···11·1···1·····1··11·1·1111····11·11·11·1111···1····111···1····1···1·1··1·1111111·1111·1··11·1··
check 46 (126)111··111··111111·11···1·1··11·1111111111111·1111111··11·1··11··1111·1··11··1·1··11·11··111·1·11111·11·1111111··1111·11···1··111··1111111··1··1111··1··1·····11·1·111··11··1··111·111·1·111·111··1·1····1·1······11
check 47 (86)··1111··1····11·1··1··111··············1···1··1··11·1·1·····111·1·····1·111·1··11·1·1·1·1·1·1·1···1·111···11··1·1··111·1······1·111·1·1·1·1·1·111·11·111··11····111·1·1··1·1·1···1··111··············1···1··1·11·1
check 48 (116)·111·1·1···1111·1··11·11···111·1·1111···11···1·1·········1····1··1····1·11··11·1·1····1··111111·11111·11111·11111·1111·1··11···11111111·1···1111··111·111·1·11·111111·1··1··11·111·11111···11····11··1··11·11··11·
check 49 (86)····11····111111·1111··11·····11···1·1···11··1·1·11······1·····1111·1··1111·11·····11······1··1····11··1·1·11···1·1·······11·11··11··11······11·11····1111··11···1·1·······11·11·11····11···11····11·····111·1··11
check 50 (100)·1·11·1·111111··11111·111··1·1·11·11······1··11···1·····1··111·1·1····11·1··1··1······1111111·1·1·1···11··1··11·11·1········1··11··1·11·1·1·1··1111··1···11111··11··1··111·1···11·111·11·····1·····1·1·1··1·11111·
check 51 (98)·1·11··1·1·1···111····1·11·····1·1·1·1111··11·11·11·1···1·········1·1·111·1···111·1···1···111··1·11·11···111··1·11·1·1·1·11··11·11111·1·····1·····1111··1·1·1···11··1····11··11··1·1111···11·1··111··1····111··111
check 52 (86)·········1·1···1·1·····1·1·11·11111··1111·1···1·1·1·1·········1·1···1··1····1·1·11·11··111····1111··1···1·1·11·······1·1···1···1····11·1··1··11111·11·1111·11··········11···1···1····1·1·1·111··1111···1·111·11·1·

The gate is read off the output rows. An output row of odd weight gives a T, a pair of output rows with odd overlap gives a CS, and an odd triple overlap gives a CCZ. The circuit is a valid factory because every such parity that touches a check row is even.

Figures of merit

inputs n210Noisy T states per run.
outputs k40
distance d3Exact, proved and witnessed here.
error coefficient A_{3}228Faults of weight 3 in the circuit shown that pass every check and damage the output. The output error is A_{3}\,p^{3} to leading order.
wires N53
checks r13
V_{\mathrm{ex}}40T states per run. CCZ counts as two.
\gamma_\rho = \log(n/V_{\mathrm{ex}})/\log d1.5094Overhead per extractable T state. Lower is better.
\gamma = \log(n/k)/\log d1.5094Overhead per output qubit.
\gamma_T = \log(n/T)/\log dnot computedOverhead per unit of T-count.
T-countnot computedMinimal T-count of the output gate. Not computable at this width.
phase-polynomial degreenot computedMinimised over changes of output basis.
effective width40Rank of the output rows modulo the check rows.

Output gate

40 T0 CS0 CCZon 40 output wires

As gates: \mathrm{T}_{0}\cdot \mathrm{T}_{1}\cdot \mathrm{T}_{2}\cdot \mathrm{T}_{3}\cdot \mathrm{T}_{4}\cdot \mathrm{T}_{5}\cdot \mathrm{T}_{6}\cdot \mathrm{T}_{7}\cdot \mathrm{T}_{8}\cdot \mathrm{T}_{9}\cdot \mathrm{T}_{10}\cdot \mathrm{T}_{11}\cdot \mathrm{T}_{12}\cdot \mathrm{T}_{13}\cdot \mathrm{T}_{14}\cdot \mathrm{T}_{15}\cdot \mathrm{T}_{16}\cdot \mathrm{T}_{17}\cdot \mathrm{T}_{18}\cdot \mathrm{T}_{19}\cdot \mathrm{T}_{20}\cdot \mathrm{T}_{21}\cdot \mathrm{T}_{22}\cdot \mathrm{T}_{23}\cdot \mathrm{T}_{24}\cdot \mathrm{T}_{25}\cdot \mathrm{T}_{26}\cdot \mathrm{T}_{27}\cdot \mathrm{T}_{28}\cdot \mathrm{T}_{29}\cdot \mathrm{T}_{30}\cdot \mathrm{T}_{31}\cdot \mathrm{T}_{32}\cdot \mathrm{T}_{33}\cdot \mathrm{T}_{34}\cdot \mathrm{T}_{35}\cdot \mathrm{T}_{36}\cdot \mathrm{T}_{37}\cdot \mathrm{T}_{38}\cdot \mathrm{T}_{39}

As monomials: 0+1+2+3+4+5+6+7+8+9+10+11+12+13+14+15+16+17+18+19+20+21+22+23+24+25+26+27+28+29+30+31+32+33+34+35+36+37+38+39 (output wires per term, separated by +)

This row stands for every gate that a CNOT circuit and S, Z and CZ gates on the outputs can reach from this one. They all prepare the same magic state. The wires are shown in the S_k-canonical frame.

Clifford correction

With every rotation a \mathrm{T}=\mathrm{diag}(1,e^{i\pi/4}), the rotations deposit the gate above times a diagonal Clifford. Undo it after the rotations and before the checks are measured, with 44 single-wire gates and 349 \mathrm{CZ} gates: \mathrm{Z} on output 0, \mathrm{Z} on output 1, \mathrm{S} on output 3, \mathrm{S} on output 5, \mathrm{S} on output 6, \mathrm{Z} on output 7, \mathrm{S}^\dagger on output 9, \mathrm{S}^\dagger on output 10, \mathrm{S}^\dagger on output 11, \mathrm{S} on output 13, \mathrm{Z} on output 14, \mathrm{S} on output 15, \mathrm{Z} on output 16, \mathrm{S}^\dagger on output 17, \mathrm{S} on output 18, \mathrm{S} on output 19, \mathrm{S} on output 21, \mathrm{S}^\dagger on output 22, \mathrm{Z} on output 23, \mathrm{S} on output 25, \mathrm{S} on output 27, \mathrm{Z} on output 28, \mathrm{S} on output 29, \mathrm{S}^\dagger on output 30 and 369 more (listed below). The circuit's action is then exactly the gate above, and an error-free run passes every check with certainty.

All 393 correction gates

\mathrm{Z} on output 0, \mathrm{Z} on output 1, \mathrm{S} on output 3, \mathrm{S} on output 5, \mathrm{S} on output 6, \mathrm{Z} on output 7, \mathrm{S}^\dagger on output 9, \mathrm{S}^\dagger on output 10, \mathrm{S}^\dagger on output 11, \mathrm{S} on output 13, \mathrm{Z} on output 14, \mathrm{S} on output 15, \mathrm{Z} on output 16, \mathrm{S}^\dagger on output 17, \mathrm{S} on output 18, \mathrm{S} on output 19, \mathrm{S} on output 21, \mathrm{S}^\dagger on output 22, \mathrm{Z} on output 23, \mathrm{S} on output 25, \mathrm{S} on output 27, \mathrm{Z} on output 28, \mathrm{S} on output 29, \mathrm{S}^\dagger on output 30, \mathrm{Z} on output 31, \mathrm{S} on output 32, \mathrm{S} on output 33, \mathrm{Z} on output 34, \mathrm{S}^\dagger on output 35, \mathrm{S} on output 36, \mathrm{S} on output 37, \mathrm{S} on output 38, \mathrm{S}^\dagger on check 40, \mathrm{S} on check 41, \mathrm{Z} on check 42, \mathrm{Z} on check 43, \mathrm{S} on check 45, \mathrm{S} on check 46, \mathrm{S} on check 47, \mathrm{Z} on check 48, \mathrm{S} on check 49, \mathrm{Z} on check 50, \mathrm{S}^\dagger on check 51, \mathrm{S} on check 52, \mathrm{CZ} on output 0, output 1, \mathrm{CZ} on output 0, output 2, \mathrm{CZ} on output 0, output 4, \mathrm{CZ} on output 0, output 7, \mathrm{CZ} on output 0, output 8, \mathrm{CZ} on output 0, output 12, \mathrm{CZ} on output 0, output 20, \mathrm{CZ} on output 0, output 23, \mathrm{CZ} on output 0, output 24, \mathrm{CZ} on output 0, output 26, \mathrm{CZ} on output 0, output 28, \mathrm{CZ} on output 0, output 31, \mathrm{CZ} on output 0, output 34, \mathrm{CZ} on output 0, output 39, \mathrm{CZ} on output 0, check 40, \mathrm{CZ} on output 0, check 42, \mathrm{CZ} on output 0, check 43, \mathrm{CZ} on output 0, check 48, \mathrm{CZ} on output 0, check 51, \mathrm{CZ} on output 0, check 52, \mathrm{CZ} on output 1, output 2, \mathrm{CZ} on output 1, output 4, \mathrm{CZ} on output 1, output 7, \mathrm{CZ} on output 1, output 8, \mathrm{CZ} on output 1, output 12, \mathrm{CZ} on output 1, output 20, \mathrm{CZ} on output 1, output 23, \mathrm{CZ} on output 1, output 24, \mathrm{CZ} on output 1, output 26, \mathrm{CZ} on output 1, output 28, \mathrm{CZ} on output 1, output 31, \mathrm{CZ} on output 1, output 34, \mathrm{CZ} on output 1, output 39, \mathrm{CZ} on output 1, check 40, \mathrm{CZ} on output 1, check 42, \mathrm{CZ} on output 1, check 43, \mathrm{CZ} on output 1, check 48, \mathrm{CZ} on output 1, check 51, \mathrm{CZ} on output 1, check 52, \mathrm{CZ} on output 2, output 4, \mathrm{CZ} on output 2, output 7, \mathrm{CZ} on output 2, output 8, \mathrm{CZ} on output 2, output 12, \mathrm{CZ} on output 2, output 14, \mathrm{CZ} on output 2, output 16, \mathrm{CZ} on output 2, output 20, \mathrm{CZ} on output 2, output 23, \mathrm{CZ} on output 2, output 24, \mathrm{CZ} on output 2, output 26, \mathrm{CZ} on output 2, output 28, \mathrm{CZ} on output 2, output 31, \mathrm{CZ} on output 2, output 34, \mathrm{CZ} on output 2, output 39, \mathrm{CZ} on output 2, check 40, \mathrm{CZ} on output 2, check 41, \mathrm{CZ} on output 2, check 45, \mathrm{CZ} on output 2, check 46, \mathrm{CZ} on output 2, check 47, \mathrm{CZ} on output 2, check 49, \mathrm{CZ} on output 2, check 51, \mathrm{CZ} on output 2, check 52, \mathrm{CZ} on output 4, output 7, \mathrm{CZ} on output 4, output 8, \mathrm{CZ} on output 4, output 12, \mathrm{CZ} on output 4, output 14, \mathrm{CZ} on output 4, output 16, \mathrm{CZ} on output 4, output 20, \mathrm{CZ} on output 4, output 23, \mathrm{CZ} on output 4, output 24, \mathrm{CZ} on output 4, output 26, \mathrm{CZ} on output 4, output 28, \mathrm{CZ} on output 4, output 31, \mathrm{CZ} on output 4, output 34, \mathrm{CZ} on output 4, output 39, \mathrm{CZ} on output 4, check 40, \mathrm{CZ} on output 4, check 41, \mathrm{CZ} on output 4, check 45, \mathrm{CZ} on output 4, check 46, \mathrm{CZ} on output 4, check 47, \mathrm{CZ} on output 4, check 49, \mathrm{CZ} on output 4, check 51, \mathrm{CZ} on output 4, check 52, \mathrm{CZ} on output 7, output 8, \mathrm{CZ} on output 7, output 12, \mathrm{CZ} on output 7, output 20, \mathrm{CZ} on output 7, output 23, \mathrm{CZ} on output 7, output 24, \mathrm{CZ} on output 7, output 26, \mathrm{CZ} on output 7, output 28, \mathrm{CZ} on output 7, output 31, \mathrm{CZ} on output 7, output 34, \mathrm{CZ} on output 7, output 39, \mathrm{CZ} on output 7, check 40, \mathrm{CZ} on output 7, check 42, \mathrm{CZ} on output 7, check 43, \mathrm{CZ} on output 7, check 48, \mathrm{CZ} on output 7, check 51, \mathrm{CZ} on output 7, check 52, \mathrm{CZ} on output 8, output 12, \mathrm{CZ} on output 8, output 14, \mathrm{CZ} on output 8, output 16, \mathrm{CZ} on output 8, output 20, \mathrm{CZ} on output 8, output 23, \mathrm{CZ} on output 8, output 24, \mathrm{CZ} on output 8, output 26, \mathrm{CZ} on output 8, output 28, \mathrm{CZ} on output 8, output 31, \mathrm{CZ} on output 8, output 34, \mathrm{CZ} on output 8, output 39, \mathrm{CZ} on output 8, check 40, \mathrm{CZ} on output 8, check 41, \mathrm{CZ} on output 8, check 45, \mathrm{CZ} on output 8, check 46, \mathrm{CZ} on output 8, check 47, \mathrm{CZ} on output 8, check 49, \mathrm{CZ} on output 8, check 51, \mathrm{CZ} on output 8, check 52, \mathrm{CZ} on output 9, output 14, \mathrm{CZ} on output 9, output 16, \mathrm{CZ} on output 9, check 41, \mathrm{CZ} on output 9, check 42, \mathrm{CZ} on output 9, check 43, \mathrm{CZ} on output 9, check 45, \mathrm{CZ} on output 9, check 46, \mathrm{CZ} on output 9, check 47, \mathrm{CZ} on output 9, check 48, \mathrm{CZ} on output 9, check 49, \mathrm{CZ} on output 10, output 14, \mathrm{CZ} on output 10, output 16, \mathrm{CZ} on output 10, check 41, \mathrm{CZ} on output 10, check 42, \mathrm{CZ} on output 10, check 43, \mathrm{CZ} on output 10, check 45, \mathrm{CZ} on output 10, check 46, \mathrm{CZ} on output 10, check 47, \mathrm{CZ} on output 10, check 48, \mathrm{CZ} on output 10, check 49, \mathrm{CZ} on output 11, output 14, \mathrm{CZ} on output 11, output 16, \mathrm{CZ} on output 11, check 41, \mathrm{CZ} on output 11, check 42, \mathrm{CZ} on output 11, check 43, \mathrm{CZ} on output 11, check 45, \mathrm{CZ} on output 11, check 46, \mathrm{CZ} on output 11, check 47, \mathrm{CZ} on output 11, check 48, \mathrm{CZ} on output 11, check 49, \mathrm{CZ} on output 12, output 14, \mathrm{CZ} on output 12, output 16, \mathrm{CZ} on output 12, output 20, \mathrm{CZ} on output 12, output 23, \mathrm{CZ} on output 12, output 24, \mathrm{CZ} on output 12, output 26, \mathrm{CZ} on output 12, output 28, \mathrm{CZ} on output 12, output 31, \mathrm{CZ} on output 12, output 34, \mathrm{CZ} on output 12, output 39, \mathrm{CZ} on output 12, check 40, \mathrm{CZ} on output 12, check 41, \mathrm{CZ} on output 12, check 45, \mathrm{CZ} on output 12, check 46, \mathrm{CZ} on output 12, check 47, \mathrm{CZ} on output 12, check 49, \mathrm{CZ} on output 12, check 51, \mathrm{CZ} on output 12, check 52, \mathrm{CZ} on output 14, output 16, \mathrm{CZ} on output 14, output 17, \mathrm{CZ} on output 14, output 20, \mathrm{CZ} on output 14, output 22, \mathrm{CZ} on output 14, output 24, \mathrm{CZ} on output 14, output 26, \mathrm{CZ} on output 14, output 30, \mathrm{CZ} on output 14, output 35, \mathrm{CZ} on output 14, output 39, \mathrm{CZ} on output 14, check 40, \mathrm{CZ} on output 14, check 41, \mathrm{CZ} on output 14, check 42, \mathrm{CZ} on output 14, check 46, \mathrm{CZ} on output 14, check 49, \mathrm{CZ} on output 14, check 50, \mathrm{CZ} on output 14, check 51, \mathrm{CZ} on output 16, output 17, \mathrm{CZ} on output 16, output 20, \mathrm{CZ} on output 16, output 22, \mathrm{CZ} on output 16, output 24, \mathrm{CZ} on output 16, output 26, \mathrm{CZ} on output 16, output 30, \mathrm{CZ} on output 16, output 35, \mathrm{CZ} on output 16, output 39, \mathrm{CZ} on output 16, check 40, \mathrm{CZ} on output 16, check 41, \mathrm{CZ} on output 16, check 42, \mathrm{CZ} on output 16, check 46, \mathrm{CZ} on output 16, check 49, \mathrm{CZ} on output 16, check 50, \mathrm{CZ} on output 16, check 51, \mathrm{CZ} on output 17, check 41, \mathrm{CZ} on output 17, check 42, \mathrm{CZ} on output 17, check 43, \mathrm{CZ} on output 17, check 45, \mathrm{CZ} on output 17, check 46, \mathrm{CZ} on output 17, check 47, \mathrm{CZ} on output 17, check 48, \mathrm{CZ} on output 17, check 49, \mathrm{CZ} on output 20, output 23, \mathrm{CZ} on output 20, output 24, \mathrm{CZ} on output 20, output 26, \mathrm{CZ} on output 20, output 28, \mathrm{CZ} on output 20, output 31, \mathrm{CZ} on output 20, output 34, \mathrm{CZ} on output 20, output 39, \mathrm{CZ} on output 20, check 40, \mathrm{CZ} on output 20, check 41, \mathrm{CZ} on output 20, check 45, \mathrm{CZ} on output 20, check 46, \mathrm{CZ} on output 20, check 47, \mathrm{CZ} on output 20, check 49, \mathrm{CZ} on output 20, check 51, \mathrm{CZ} on output 20, check 52, \mathrm{CZ} on output 22, check 41, \mathrm{CZ} on output 22, check 42, \mathrm{CZ} on output 22, check 43, \mathrm{CZ} on output 22, check 45, \mathrm{CZ} on output 22, check 46, \mathrm{CZ} on output 22, check 47, \mathrm{CZ} on output 22, check 48, \mathrm{CZ} on output 22, check 49, \mathrm{CZ} on output 23, output 24, \mathrm{CZ} on output 23, output 26, \mathrm{CZ} on output 23, output 28, \mathrm{CZ} on output 23, output 31, \mathrm{CZ} on output 23, output 34, \mathrm{CZ} on output 23, output 39, \mathrm{CZ} on output 23, check 40, \mathrm{CZ} on output 23, check 42, \mathrm{CZ} on output 23, check 43, \mathrm{CZ} on output 23, check 48, \mathrm{CZ} on output 23, check 51, \mathrm{CZ} on output 23, check 52, \mathrm{CZ} on output 24, output 26, \mathrm{CZ} on output 24, output 28, \mathrm{CZ} on output 24, output 31, \mathrm{CZ} on output 24, output 34, \mathrm{CZ} on output 24, output 39, \mathrm{CZ} on output 24, check 40, \mathrm{CZ} on output 24, check 41, \mathrm{CZ} on output 24, check 45, \mathrm{CZ} on output 24, check 46, \mathrm{CZ} on output 24, check 47, \mathrm{CZ} on output 24, check 49, \mathrm{CZ} on output 24, check 51, \mathrm{CZ} on output 24, check 52, \mathrm{CZ} on output 26, output 28, \mathrm{CZ} on output 26, output 31, \mathrm{CZ} on output 26, output 34, \mathrm{CZ} on output 26, output 39, \mathrm{CZ} on output 26, check 40, \mathrm{CZ} on output 26, check 41, \mathrm{CZ} on output 26, check 45, \mathrm{CZ} on output 26, check 46, \mathrm{CZ} on output 26, check 47, \mathrm{CZ} on output 26, check 49, \mathrm{CZ} on output 26, check 51, \mathrm{CZ} on output 26, check 52, \mathrm{CZ} on output 28, output 31, \mathrm{CZ} on output 28, output 34, \mathrm{CZ} on output 28, output 39, \mathrm{CZ} on output 28, check 40, \mathrm{CZ} on output 28, check 42, \mathrm{CZ} on output 28, check 43, \mathrm{CZ} on output 28, check 48, \mathrm{CZ} on output 28, check 51, \mathrm{CZ} on output 28, check 52, \mathrm{CZ} on output 30, check 41, \mathrm{CZ} on output 30, check 42, \mathrm{CZ} on output 30, check 43, \mathrm{CZ} on output 30, check 45, \mathrm{CZ} on output 30, check 46, \mathrm{CZ} on output 30, check 47, \mathrm{CZ} on output 30, check 48, \mathrm{CZ} on output 30, check 49, \mathrm{CZ} on output 31, output 34, \mathrm{CZ} on output 31, output 39, \mathrm{CZ} on output 31, check 40, \mathrm{CZ} on output 31, check 42, \mathrm{CZ} on output 31, check 43, \mathrm{CZ} on output 31, check 48, \mathrm{CZ} on output 31, check 51, \mathrm{CZ} on output 31, check 52, \mathrm{CZ} on output 34, output 39, \mathrm{CZ} on output 34, check 40, \mathrm{CZ} on output 34, check 42, \mathrm{CZ} on output 34, check 43, \mathrm{CZ} on output 34, check 48, \mathrm{CZ} on output 34, check 51, \mathrm{CZ} on output 34, check 52, \mathrm{CZ} on output 35, check 41, \mathrm{CZ} on output 35, check 42, \mathrm{CZ} on output 35, check 43, \mathrm{CZ} on output 35, check 45, \mathrm{CZ} on output 35, check 46, \mathrm{CZ} on output 35, check 47, \mathrm{CZ} on output 35, check 48, \mathrm{CZ} on output 35, check 49, \mathrm{CZ} on output 39, check 40, \mathrm{CZ} on output 39, check 41, \mathrm{CZ} on output 39, check 45, \mathrm{CZ} on output 39, check 46, \mathrm{CZ} on output 39, check 47, \mathrm{CZ} on output 39, check 49, \mathrm{CZ} on output 39, check 51, \mathrm{CZ} on output 39, check 52, \mathrm{CZ} on check 40, check 41, \mathrm{CZ} on check 40, check 45, \mathrm{CZ} on check 40, check 46, \mathrm{CZ} on check 40, check 47, \mathrm{CZ} on check 40, check 49, \mathrm{CZ} on check 40, check 51, \mathrm{CZ} on check 40, check 52, \mathrm{CZ} on check 41, check 42, \mathrm{CZ} on check 41, check 46, \mathrm{CZ} on check 41, check 49, \mathrm{CZ} on check 41, check 50, \mathrm{CZ} on check 41, check 51, \mathrm{CZ} on check 42, check 43, \mathrm{CZ} on check 42, check 46, \mathrm{CZ} on check 42, check 48, \mathrm{CZ} on check 42, check 49, \mathrm{CZ} on check 42, check 50, \mathrm{CZ} on check 42, check 52, \mathrm{CZ} on check 43, check 45, \mathrm{CZ} on check 43, check 47, \mathrm{CZ} on check 43, check 50, \mathrm{CZ} on check 43, check 52, \mathrm{CZ} on check 45, check 47, \mathrm{CZ} on check 45, check 48, \mathrm{CZ} on check 45, check 50, \mathrm{CZ} on check 45, check 51, \mathrm{CZ} on check 46, check 49, \mathrm{CZ} on check 46, check 50, \mathrm{CZ} on check 46, check 51, \mathrm{CZ} on check 47, check 48, \mathrm{CZ} on check 47, check 50, \mathrm{CZ} on check 47, check 51, \mathrm{CZ} on check 48, check 50, \mathrm{CZ} on check 48, check 52, \mathrm{CZ} on check 49, check 50, \mathrm{CZ} on check 49, check 51, \mathrm{CZ} on check 51, check 52

Or no Clifford gate at all: run \mathrm{T}^3 on rotations 2, 9, 21, 25, 27, 35, 40, 45, 57, 59, 60, 61, 64, 65, 69, 71, 73, 74, 77, 79, 82, 83, 86, 87 and 55 more (listed below); \mathrm{T}^5 on rotations 1, 10, 11, 13, 14, 15, 18, 29, 38, 41, 111; \mathrm{T}^\dagger on rotations 0, 5, 7, 8, 20, 22, 30, 31, 33, 34, 39, 55, 80, 113, and \mathrm{T} on the rest. A power of \mathrm{T} costs one magic state, like \mathrm{T} itself. For a code this is a transversal \mathrm{T} gate, \mathrm{T}^j on each qubit.

All 79 rotations at T³

2, 9, 21, 25, 27, 35, 40, 45, 57, 59, 60, 61, 64, 65, 69, 71, 73, 74, 77, 79, 82, 83, 86, 87, 88, 89, 91, 93, 95, 96, 97, 100, 101, 106, 108, 110, 115, 116, 118, 124, 126, 128, 134, 137, 139, 142, 144, 146, 148, 151, 152, 153, 156, 163, 164, 166, 167, 168, 169, 170, 176, 177, 183, 184, 185, 187, 188, 189, 191, 192, 195, 197, 198, 199, 200, 201, 202, 204, 205

Rotations are numbered by matrix column, from 0. The correction comes from the columns by the weight expansion of Bravyi and Haah (2012). verify_catalog.py re-derives it, checks the corrected circuit's phase input by input against the gate, and on small circuits simulates it gate by gate.

Convert the output with CNOT and S

A CNOT circuit on the outputs, followed by S, Z and CZ gates, turns this factory's gate into other gates that prepare the same magic state. Type a target gate to check whether it is reachable. If it is, you get the circuit, checked on every basis state.

Or list every gate it reaches. The list is built in your browser and saved as a text file, one gate per line.

Python code that writes the same list
import json, urllib.request

NAMES = {1: "T", 2: "CS", 3: "CCZ"}


def equivalent_gates(k, terms, limit=1_000_000):
    """Every gate a CNOT circuit and S, Z, CZ gates reach from this one.

    A gate is its set of terms mod Clifford: {a} is T, {a,b} is CS and
    {a,b,c} is CCZ on those outputs. CNOT c->t sends x_t to
    x_c + x_t - 2 x_c x_t, so a term on t but not c adds the term with t
    replaced by c, and, below CCZ, the term with c added.
    """
    start = frozenset(frozenset(m) for m in terms)
    seen, queue = {start}, [start]
    for gate in queue:
        for c in range(k):
            for t in range(k):
                if c == t:
                    continue
                new = set(gate)
                for m in gate:
                    if t in m and c not in m:
                        new ^= {(m - {t}) | {c}}
                        if len(m) < 3:
                            new ^= {m | {c}}
                new = frozenset(new)
                if new not in seen:
                    seen.add(new)
                    queue.append(new)
                    if len(seen) > limit:
                        raise RuntimeError(f"more than {limit} equivalent gates")
    return seen


def name(gate, k):
    sep = "," if k > 10 else ""
    terms = sorted((sorted(m) for m in gate), key=lambda m: (len(m), m))
    return "·".join(NAMES[len(m)] + sep.join(map(str, m)) for m in terms) or "I"


url = "https://shubhamj810.github.io/magic-state-factory-repo/data/factories/210.40.3.a.json"
record = json.load(urllib.request.urlopen(url))
k = record["parameters"]["k"]
gates = equivalent_gates(k, record["gate"]["monomials"])
lines = sorted(name(g, k) for g in gates)
print(len(lines), "equivalent gates")
with open("equivalent-gates-210.40.3.a.txt", "w") as out:
    out.write("\n".join(lines) + "\n")

Distance

d=3, exact. No fault of weight below 3 is both undetectable and damaging. One of weight 3 is.

Exactly 228 faults of weight 3 are undetectable and damaging, so A_{3}=228. If each rotation fails independently with probability p, an accepted run has a damaged output with probability A_{3}\,p^{3}+O(p^{4}).

A damaging fault of weight 3 sits on columns [0,21,75].

A fault is a set of faulty columns. It is undetectable when their check parts cancel, and damaging when their output parts do not.

Original source

S. P. Jain, A. Wills, and S. Singh, "Symmetry and AI-assisted discovery of magic-state factories," arXiv:2610.06535 (2026). arxiv.org/abs/2610.06535

Provenance

a paired-column contraction of an existing triorthogonal code, in the pure-T distillation-exponent release (gamma = log(n/k)/log d, n < 1000); a verified witness, not a maximum -- only the release's own frontier points are the lowest gamma it knows

Found by: AI search (gamma frontier): contraction of an existing code

Discovery: AI search

1 recorded source
[
 {
  "regime": "AI search (gamma frontier): contraction of an existing code",
  "file": "gamma_frontier_release/protocols/paired_column_contraction_210_40_3_1aea0b9276aa.json",
  "label": "paired_column_contraction_210_40_3_1aea0b9276aa",
  "N": 53,
  "d": 3,
  "d_is_exact": true,
  "origin": "Existing triorthogonal code with physical paired-column contraction",
  "provenance": "gamma frontier release record paired_column_contraction_210_40_3_1aea0b9276aa (source catalogue id paired_column_contraction_210_40_3_1aea0b9276aa); construction: Existing triorthogonal code with physical paired-column contraction; matrix 1aea0b9276aa33de...; the release states d=3 (exact), gamma=1.5093842420185073 (exact)",
  "notes": "matrix gamma_frontier_release/matrices/1aea0b9276aa33dedf466f4a3e1ee1574f2d168355470150f67ab844fb8b0c49.matrix"
 }
]

Use it in code

Load this factory and re-derive every number from its circuit with the catalogue's own verifier.

# git clone https://github.com/shubhamj810/magic-state-factory-repo
# cd magic-state-factory-repo/master_catalog
import json, urllib.request
import verify_catalog as VC

url = "https://shubhamj810.github.io/magic-state-factory-repo/data/factories/210.40.3.a.json"
record = json.load(urllib.request.urlopen(url))
columns = record["circuit"]["columns"]
k, N = record["parameters"]["k"], record["parameters"]["N"]

facts, problems = VC.derive(columns, k, N)    # gate, T-count, degree, ...
distance = VC.measure_distance(columns, k, N)  # exhaustive fault search
print(facts["gate_human"], distance["d_exact"] or distance["d_at_least"],
      problems or "verified")

Raw record

Show the JSON this page is built from

Source: master_catalog.json, factories[850].