魔法态蒸馏

魔法态蒸馏（magic state distillation，MSD）是容错架构里非 Clifford 门的补给线：把物理设备能制备的噪声魔法态（|T⟩ 或 |H⟩ 型叠加态）作为原料，用理想 Clifford 操作加测量做后选择提纯。Bravyi 与 Kitaev 2004 年证明只要原料沿"魔法方向"的极化超过阈值，纯魔法态即可蒸馏出来，从而通用计算成立；其 15-to-1 协议（15 个噪声态换 1 个好态、输出误差约 35ε³）至今仍是架构基准，Meier–Eastin–Knill 的 10-to-2 与 Bravyi–Haah 的三正交码族把开销指数从 2.47 压到趋近 1.6。

魔法态与 Bloch 八面体

Clifford 群把单比特 Bloch 球的八面体对称群映到自身，因此稳定子态（Clifford 轨道里的 $|0⟩,|+⟩$ 等）之外的方向分两类：$|T⟩=T|+⟩$ 型（8 个，对应 $120^{\circ }$ 旋转轴）与 $|H⟩=H|+⟩$ 型（12 个，对应 $180^{\circ }$ 轴）。Gottesman–Knill 定理说 Clifford 操作可经典模拟，所以这些非稳定子”魔法”方向正是通用性的最小缺口Bravyi & Kitaev 2005 (PRA)。

图:单比特稳定子态在 Bloch 球上构成内接八面体的顶点;T 型(120° 轴)与 H 型(180° 轴)魔法态占据对称轴上八面体之外的位置(来源:QE296)

蒸馏的可行性有边界：$T$ 方向投影超过 $2F_T^2-1\approx 0.655$（$F_T=0.910$）、$H$ 方向超过 $2F_H^2-1\approx 0.718$（$F_H\approx 0.9272$）即可启动；Reichardt 后来证明 $T$ 型的普适阈值为 $(1+1/\sqrt{2})/2\approx 0.854$，而 Campbell–Browe 证明 $|T⟩$ 在稳定子边界附近不可蒸馏（no-go）。

五比特码：5ε² 的热身

对角化于 $A$ 基的输入 $\rho =(1-ϵ)|T_0⟩⟨T_0|+ϵ|T_1⟩⟨T_1|$，在 $[[5,1,3]]$ 完美码上测 $Z$ 综合征、校正、再测 $A$ 综合征（非平凡即弃）、解码，得 ${ϵ}_{\mathrm{out}}\approx 5{ϵ}^2$，阈值 ${ϵ}_0=(1-\sqrt{3/7})/2\approx 0.1732$；成功概率从 $1/6$（$ϵ=0$）单调降到 $1/16$，产率约 $1/30$；递归给出 $ϵ\to \frac{1}{5}(5ϵ{)}^{2^k}$Bravyi & Kitaev 2005 (PRA)。平方压制已经可用，但三次压制的 15-to-1 每比特成本更低。

15-to-1：码与横向门

码从刺穿 Reed–Muller 码构造。记 $[f]\in {\mathbb{F}}_2^{15}$ 为布尔函数 $f$ 在除 $0000$ 外全部输入上的取值表，取两个线性空间

$${\mathcal{L}}_1=\mathrm{span}\{[x_1],[x_2],[x_3],[x_4]\},{\mathcal{L}}_2={\mathcal{L}}_1+\mathrm{span}\{[x_i{x}_j]{\}}_{i<j},$$

即一阶/二阶刺穿 RM 码（$\dim {\mathcal{L}}_1=4$、$\dim {\mathcal{L}}_2=10$；${\mathcal{L}}_1$ 是 $[15,5,8]$ 单纯形码的对偶）。考虑 $\mathrm{CSS}({\sigma }^z,{\mathcal{L}}_2;A,{\mathcal{L}}_1)$，其中 $A=({\sigma }^x+{\sigma }^y)/\sqrt{2}$ 是非 Pauli 的 Clifford 算符。关键引理链Bravyi & Kitaev 2005 (PRA)：

• Lemma 1：${\mathcal{L}}_1$ 中向量重量 $\equiv 0(\mathrm{mod}8)$；${\mathcal{L}}_2$ 全偶重；${\mathcal{L}}_1^{\perp }={\mathcal{L}}_2\oplus [1]$，${\mathcal{L}}_2^{\perp }={\mathcal{L}}_1\oplus [1]$；任意 $u,v\in {\mathcal{L}}_1$ 有 $|u\cdot v|\equiv 0(\mathrm{mod}4)$。
• Lemma 2（横向门）：在 $Z$ 综合征平凡的子空间上，$A(u)$ 与 Pauli 算符 ${\sigma }^x(u)$ 的作用重合（相位因子 $e^{i\pi |u|/4-i\pi |u\cdot v|/2}=1$ 正是由模 8/模 4 条件保证）——于是非 Pauli 稳定子可以用 Pauli 替身测量，整个综合征提取留在 ${\mathcal{O}}_{\mathrm{ideal}}$（理想 Clifford）内。等价地：逐位作用 $W^{\otimes 15}$（$W{\sigma }^x{W}^{†}=A$，$W\notin$ Clifford）保持码空间并在其上实现非 Clifford 逻辑门。后世约定（Bravyi–Haah 等）把这个门直接叫横向 $T$：同一 $[[15,1,3]]$ 码上 ${⨂}_j{T}^{u_j}$ 是逻辑 $T$。

蒸馏子程序六步：①测 ${\sigma }^z([x_j])$、${\sigma }^z([x_i{x}_j])$ 得 $Z$ 综合征 $\mu$；②找 $w$ 使 $(w,v)=\mu (v)$；③施加校正 $A(w{)}^{†}$；④测 ${\sigma }^x([x_j])$ 得 $A$ 综合征 $\eta$；⑤$\eta \ne 0$ 即宣布失败；⑥Clifford 解码到单比特。容错性来自横向性：单物理故障在 15 个输入上至多造成单比特错误，天然不级联。

15-to-1：误差多项式的完整推导

输入 ${\rho }^{\otimes 15}={\sum }_u{ϵ}^{|u|}(1-ϵ{)}^{15-|u|}|A_u⟩⟨A_u|$ 后，接受概率与输出误差都是重量枚举多项式（weight enumerator）的取值：

$$p_s=W_{{\mathcal{L}}_1^{\perp }}(ϵ,1-ϵ),{ϵ}_{\mathrm{out}}=\frac{W_{{\mathcal{L}}_2}(ϵ,1-ϵ)}{W_{{\mathcal{L}}_1^{\perp }}(ϵ,1-ϵ)},$$

其中 $W_{\mathcal{L}}(x,y)={\sum }_{u\in \mathcal{L}}{x}^{n-|u|}{y}^{|u|}$。用 MacWilliams 恒等式

$$W_{\mathcal{L}}(x,y)=\frac{1}{|{\mathcal{L}}^{\perp }|}{W}_{{\mathcal{L}}^{\perp }}(x+y,x-y),$$

加上 ${\mathcal{L}}_2^{\perp }={\mathcal{L}}_1\oplus [1]$ 与 ${\mathcal{L}}_1$ 全偶重，两式都化到 ${\mathcal{L}}_1$ 的枚举多项式上——而单纯形码的枚举多项式简单到一眼可读：

$$W_{{\mathcal{L}}_1}(x,y)=x^{15}+15x^7{y}^8.$$

代入即得闭式（原文式 35–36）：

$$p_s=\frac{1+15(1-2ϵ{)}^8}{16},$$ $${ϵ}_{\mathrm{out}}=\frac{1-15(1-2ϵ{)}^7+15(1-2ϵ{)}^8-(1-2ϵ{)}^{15}}{2(1+15(1-2ϵ{)}^8)}.$$

小 $ϵ$ 展开验证 ${ϵ}^0,{ϵ}^1,{ϵ}^2$ 项逐项相消、${ϵ}^3$ 系数为 $1120/32$，给出著名的

$${ϵ}_{\mathrm{out}}\approx 35{ϵ}^3,$$

——系数 35 完全由 $x^{15}+15x^7{y}^8$ 的结构决定。数值解 ${ϵ}_{\mathrm{out}}(ϵ)=ϵ$ 给阈值 ${ϵ}_0\approx 0.141$。递归 $k$ 层后

$${ϵ}_{\mathrm{out}}(n,ϵ)\sim (\sqrt{35}ϵ{)}^{n^{\xi }},\xi =\frac{1}{{\log }_{3}15}\approx 0.404,$$

即蒸馏到精度 $1/ϵ$ 需要输入态数 $n=O({\log }^{\gamma }(1/ϵ))$、$\gamma ={\log }_{3}15\approx 2.47$。

图:15-to-1 协议的输出误差 ε_out(ε) 与成功概率 p_s(ε)——ε_out 曲线与对角线的交点即阈值 ε₀≈0.141(来源:QE296)

记号警示：Bravyi–Kitaev 原文 Sec. VI 用上述 $A$ 基蒸馏的是 $H$ 型态（$|A_0⟩=e^{i\pi /8}H{K}^{†}|H⟩$）；Bravyi–Haah 之后的文献把同一 $[[15,1,3]]$ 码与横向 $T$ 用于 $T$ 型态 $|A⟩=T|+⟩$，输入输出误差关系同为 $35p^3$。两文的”A 基”字母同名不同义，读原始文献时务必消歧。

改进史：MEK 10-to-2 与 Bravyi–Haah 三正交码

MEK 10-to-2：Meier、Eastin 与 Knill 用 $[[4,2,2]]$ 检错码蒸馏 $|H⟩$ 态——该码稳定子 $\{XXXX,ZZZZ\}$，横向 $H^{\otimes 4}={\bar{H}}_1{\bar{H}}_2\cdot {\mathrm{SWAP}}_{12}$。4 个受控-Hadamard 门加 8 个资源态完成编码的 ${\bar{H}}_1{\bar{H}}_2$ 测量，解码后测综合征，两者皆平凡才接受；10 个输入换 2 个输出，$e(p)=u(p)/a(p)$，其中 $u(p)=9p^2-56p^3+\cdots +32p^8$、$a(p)=1-10p+58p^2-\cdots +64p^8$Meier, Eastin & Knill 2013 (QIC)。Knill 早年的 7 比特 Steane 码版本与之等价（Reichardt 证明）；Reichardt 的 7-to-1 把 $|H⟩$ 阈值推到 0.146。

图:10-to-2(MEK)与 15-to-1 的条件误差 e(p) 对比——MEK 每轮产率更高,两者级联组合显著降低总开销(来源:QE301)

Bravyi–Haah 三正交码：二进制矩阵 $G$ 若任意两行、三行支撑集交集均为偶（$|g^a\cdot g^b|\equiv |g^a\cdot g^b\cdot g^c|\equiv 0\mathrm{mod}2$）则称三正交；奇重行张出 $k$ 个逻辑比特，码 $\mathrm{CSS}(X,G_0;Z,G^{\perp })$ 上 ${⨂}_j{T}^{u_j}$ 经 Clifford 修正（只用 $\Lambda (Z)$ 与 $S$）化为逻辑 $T$Bravyi & Haah 2012 (PRA)。距离 2、码率 1/3 的码族给出 $3k+8\to k$ 协议（$p\to O(p^2)$），开销指数 $\gamma ={\log }_2((3k+8)/k)\to {\log }_{2}3\approx 1.585$；并构造出首个 $[[49,1,5]]$ 横向 $T$ 码。成本表（输入误差 $10^{-2}$、目标 $10^{-12}$）：15-24-36 协议序列需 $C=187.9$ 个输入态，对比纯 15-to-1/10-to-2 组合的 $418$——约 $2.2\times$ 缩减；目标 $10^{-30}$ 时为 1301 对 3483。任何蒸馏协议 $\gamma \ge 1$ 的猜想也在此文论证。

图:三正交码蒸馏子程序——编码 |+⟩ 态、逐位消耗 |A⟩ 完成横向 T、测 X 型稳定子后选择、解码输出 k 个提纯态(来源:QE298)

2025–2026 年的新方向继续压常数：unfolding 协议把三正交码折叠出更优距离-码率权衡2025 (arXiv)，零级 CCZ 蒸馏瞄准直接产出 CCZ 资源态2026 (arXiv)。另一条与蒸馏竞争的路线——在硬件上直接”培育”魔法态——见magic-state-injection；两者在资源账本中的位置见ft-resource-overhead。

延伸阅读

• magic-state-injection
• teleportation-gates
• transversal-gates
• eastin-knill
• quantum-reed-muller
• ft-resource-overhead