Jones 多项式把结图的拓扑信息编码为一个 Laurent 多项式:逐项展开会指数爆炸,精确计算在大多数求值点上落在 -hard 的难度量级。然而在特定的单位根处,同一个数值可以改写成 braid group 酉表示的矩阵元或 Markov trace——而“估计一个酉矩阵元”正是量子线路(配合 Hadamard test)的母语。本词条从零开始推导这条链路:先讲清结不变量为什么难算,再把结代数化为 braid word 与 Temperley—Lieb 代数,然后在单位根上构造 path model 酉表示,得到 Aharonov—Jones—Landau(AJL)的多项式时间量子算法;最后讨论两处“完备性”边界——plat closure 的加性近似是 BQP-complete、trace closure 的归一化迹估计是 DQC1-complete——以及把 Jones “范畴化”之后的 Khovanov 同调在 2025 年出现的量子算法与它的适用条件。
为什么值得关心? 结不变量是区分不同结与链环的第一道筛子:若两个链环的不变量不同,它们一定不同痕。判定“两个结是否等价”(结问题)在经典上是可判定的,但已知的通用判定过程(如 Haken 的 normal surface 理论)复杂度高得难以使用;而不变量这一侧,精确计算又普遍 -hard。量子算法给出的第三条路是特定归一化尺度下的加性近似——它既不是精确值,也不是相对误差,理解“这到底买到了什么”是本词条贯穿始终的主线,也是本词条 index 页反复强调的原则:一个 -hard 精确量不因存在某种加性量子估计就变得普遍易算。
历史脉络(与文末参考文献对应)。 结不变量一侧:Reidemeister 于 1926 年证明同一链环的任意两张结图可由三类局部移动互变,这给出了“图上计算不变量”的合法性;Alexander 于 1928 年给出第一个多项式不变量(他的辫表示定理则早在 1923 年,见第 2.2 节);Jones 于 1984 年在完全不同的语境(von Neumann 代数的子因子理论)中发现了以他命名的新多项式,并因此获得 Fields 奖;Temperley 与 Lieb 早在 1971 年研究冰模型传递矩阵时就写下了后来以他们命名的代数,Kauffman 于 1987 年用它给出 Jones 多项式的括号演算——本词条第 3 节的整套计算就来自这条路线。物理一侧:Witten 于 1989 年指出 Jones 多项式是 Chern—Simons 拓扑量子场论中 Wilson 圈的期望值,Reshetikhin—Turaev 随后给出了严格的数学构造;Freedman、Kitaev、Larsen 与 Wang 在 2000 年前后证明了这类 TQFT 与量子计算可以互相模拟——这既预示了拓扑量子计算(用 anyon 的编织做[通用量子门](../fundamentals/universal-gate-set)),也隐含了“逼近 Jones 多项式与量子计算一样难”的完备性结论。算法一侧:Aharonov—Jones—Landau 于 2005 年(STOC 2006)避开 TQFT、直接用 Temperley—Lieb 代数的 path model 表示,写出了显式的多项式时间量子算法(Zoo 2、4、41、42 对应这条线);Shor 与 Jordan 随后证明某些 trace closure 的 Jones 近似在 DQC1 模型中完备(Zoo 83);HOMFLY 与 quantum-double 型链环不变量有类似的表示论算法(Zoo 93、174);与 Tutte/Potts 配分函数的联系见 Zoo 3 及本词条 配分函数量子算法 词条。2025 年,Schmidhuber 等人给出了第一个 Khovanov 同调的量子算法(Zoo 510—511),其核心是把同调秩的计算化为 Hodge Laplacian 近零谱质量的估计。
前置知识。 本词条假定读者已完成第 1—8 章:量子态与测量、Hadamard 门、Grover 搜索、相位估计与振幅估计、基本的线性代数(内积、投影、酉算子)。Hadamard test 与 Hoeffding 不等式会在这里自足地推导,不需要回头翻;但振幅估计只引用结论“ 次相干调用给出加性 估计”。拓扑一侧不要求任何先验知识:结、链环、braid、Reidemeister 移动都会就地定义。
1. 问题从哪里来:结、链环与 Jones 多项式
1.1 结与链环
结 (knot) 是三维空间 (或三维球面 )中一条自身不相交的闭曲线;链环 (link) 是若干条互不相交的闭曲线的并。两个链环等价(同痕,ambient isotopic),如果其中一个可以通过连续变形(曲线始终不自交、互不穿透)变成另一个。最简单的结是未打结的圆周,称为平凡结 (unknot);最简单的非平凡结是三叶结 (trefoil),最少需要三个交叉才能画出。
要把链环交给算法,必须先拍一张二维照片:取一个投影方向,把链环投影到平面上,并在每个交叉处记录哪条线在上方。这样得到的四价平面图称为结图 (knot diagram);一张结图的规模用它包含的交叉数 衡量。同一张图里“谁上谁下”是投影的真实信息,而“交叉的位置与连法”则依赖于投影方向——Reidemeister 定理(1926)把这些自由度完全刻画清楚:
定理 (Reidemeister)。 同一个链环的任意两张结图,可以通过反复施加以下三类局部移动互相变换;反之,能如此互变的两张图代表同一个链环。
- R1:增删一个“卷”(一段线自己跨过自己形成的单个交叉);
- R2:把一条压在另一条上的两平行线段整体滑走(同时增删两个交叉);
- R3:三段线之间的三角形滑动(交叉数不变,只改变哪个交叉在左边)。
这条定理的意义是算法性的:任何“逐交叉计算的图泛函”只要在 R1、R2、R3 下不变,就是链环的不变量。后文的 Kauffman 括号正是这样被逐条验证的。
1.2 Kauffman 括号与 Jones 多项式
Kauffman 括号 (Kauffman bracket) 是定义在(无向)结图 上、取值于 的函数,由三条规则完全确定:
- 归一化:(无交叉的圆周取值 );
- 脱圈规则:,其中 , 表示 与一个不相交的圆圈之并;
- 展开规则(skein 关系):在每个交叉处, 即把该交叉的四个端点两两重连:一种重连方式是“两条竖线”(平直型),另一种是“上弧接上弧、下弧接下弧”(转折型)。哪种平滑配 、哪种配 ,由交叉处上方线段的旋向决定;不同教材的约定可能恰好交换 与 (因此交换了链环的两种手性),本词条固定一种约定,结论不受影响。
反复应用规则 3,可以把任何结图化到没有交叉的情形,再用规则 1、2 读出 的幂——所以括号总等于 的幂的整系数组合,是一个 Laurent 多项式。但它是链环不变量吗?逐条对照 Reidemeister 移动(第 3.3 节将把这些验证代数化)会发现:
- R2、R3 下不变:局部计算恰好化为 Temperley—Lieb 代数的恒等式(第 3 节的主题);
- R1 下会差一个因子:单个卷的括号是 而不是 (第 8 节将实际算出)。
修补办法是经典的:给结图 定义拧数 (writhe) 为全部交叉的带号总和(每个交叉按“右手拧”记 、“左手拧”记 ),然后带上一个补偿因子。定义
一个卷会同时让 变化 、让 乘上 (这正是 R1 失败的形态),两者恰好相消;R2 不改变 且括号不变;R3 两者都不变。因此 在全部三类移动下不变,是链环不变量。以 换元,就得到 Jones 多项式:
它也等价地由 Jones 原始的 skein 关系 (配合 )定义,其中 是只在同一个交叉处取三种不同处理的图。对多分量链环,(半整数幂是正常的);对结(单分量)则是 的整幂 Laurent 多项式。两个基本值值得记住:,二分量平凡链环 (第 8 节会顺带验证后者)。
1.3 经典算法能做到什么程度
精确计算的指数瓶颈。 按定义计算括号: 个交叉每个有两种平滑方式,展开成一棵深度为 的二叉树,叶子数 ,每个叶子是若干不相交圆圈之并、其括号是 的幂。这个“skein 树”算法的时间是 :关于交叉数指数。动态规划可以省掉重复子树,但对一般输入仍然指数。
不能指望多项式精确算法。 Jones 多项式与图论的 Tutte 多项式血脉相连:对交错结,Jones 在 处的值可以用平面图的 Tutte 多项式在相应点读出(媒介构造见本词条 配分函数量子算法 第 5 节)。而 Jaeger—Vertigan—Welsh 的一族结果表明,Tutte 多项式在大多数固定点上的精确求值是 -hard 的;相应地,Jones 多项式在大多数固定求值点上的精确计算也被证明(或普遍相信,视具体点集而定)落在同样的难度量级。只有少数特殊点(如 等退化情形)已知有多项式时间经典算法。本词条反复出现的告诫在这里原样适用:-hard 的精确量不会被任何通用多项式算法(经典或量子)解决。
近似怎么办? 经典侧对这类拓扑量缺少通用的随机近似格式(对比:铁磁 Potts 配分函数有 FPRAS,但那是“正权重、可混合”的结构带来的)。量子侧的机会在于:把不变量改写成酉算子的矩阵元,矩阵元可以用 Hadamard test 以 次采样(或振幅估计的 次)估计到加性 。注意这个承诺的形状:加性、相对于某个自然尺度。若该尺度本身随输入指数增长,加性近似并不给出相对误差,更不给出精确系数——第 5.4 节将把这句话量化。
1.4 本词条路线图
- 第 2 节:把链环输入代数化为 braid word,区分 trace closure 与 plat closure 两种读出方式;
- 第 3 节:Temperley—Lieb 代数与 Kauffman 的字典“交叉 平滑的线性组合”,验证 确实是 braid 群表示;
- 第 4 节:在单位根上构造 path model,使表示变酉,并逐项解释 的电路复杂度;
- 第 5 节:plat closure 算法(Hadamard test + 采样复杂度 + 误差尺度);
- 第 6 节:trace closure 算法(Markov trace、提纯、DQC1);
- 第 7 节:为什么这些近似问题对量子计算是完备的;
- 第 8 节:一个从头算到尾的数值例子( 上的 Hopf 链环与三叶结);
- 第 9 节:推广到 HOMFLY 与 Khovanov 同调,以及 2025 年的 Hodge Laplacian 量子算法。
2. 把结代数化:braid 群与两种闭合
2.1 Braid 群
直观上,一个 braid(辫) 是 股竖直的线:顶部有 个端点、底部有 个端点,每股从顶部某点沿竖直方向(允许左右绕行、永不回头)到达底部某点,两股交叉时记录谁在上。Artin braid 群 是这些辫在“端点固定、连续形变”意义下的等价类构成的群,乘法为竖直拼接(先做上段再做下段)。它由生成元 给出,其中 表示第 股与第 股的一次正交叉(第 股从上方跨过), 表示对应的负交叉。定义关系为:
第二条即著名的辫关系 (braid relation) :两股在两处交叉后“谁在内谁在外”的两种顺序给出同一个辫。注意 是无限群( 的幂互不相同),而把每个 映为对换 得到满同态 :辫比置换记录了更多信息——正是“交叉怎么绕”的那部分信息。
指数和 (exponent sum) 是 的同态:,其中 。第 1.2 节的 writhe 对“全部股向下定向”的辫图恰好等于 ——这个量稍后将进入归一化因子。
2.2 Alexander 定理与两种闭合
链环与辫由闭合 (closure) 操作连接起来。
定理 (Alexander, 1923)。 任何有向链环都同痕于某个辫的trace closure。
Trace closure(迹闭合):把 股辫的第 个顶端与第 个底端用一条不交叉的弧连接(共 条弧,全部绕在辫的外侧),得到一个闭链环。例如 2 股辫 的 trace closure 是 型环面链环: 是 Hopf 链环、 是三叶结(第 8 节的例子正是它们)。
Plat closure(plat 闭合):要求股数 为偶数,把相邻端点成对封闭——顶部把 接 、 接 、……,底部同样。这相当于在辫的上下各加一排“杯”与“帽”。并非每个链环都直接是某个辫的 plat closure,但任何链环都适当加厚后可用偶数股的 plat 表示;对本词条的算法而言,输入就限定为“已经写成辫”的情形。
为什么要区分两种闭合? 因为它们对应两种完全不同的读出方式:
- plat closure 矩阵元。杯和帽选定了一对确定的初态/末态(后文的 ),链环的不变量正比于辫表示算子在该态上的矩阵元 。矩阵元用一次 Hadamard test 就能采样——这是 AJL 算法的读出,也是 BQP 完备性的载体。
- trace closure 迹。闭合把“第 股接回第 股”翻译成“对股指标求和”,对应表示算子的(加权)迹。归一化迹是 DQC1 模型的母语(第 6.3 节),这是 Shor—Jordan 结果的载体。
输入表示与复杂度参数。 本词条的输入是长度 的 braid word
连同股数 与单位根参数 (第 4 节)。输入大小为 。对比第 1.3 节以交叉数 计的结图:braid word 的长度 正是闭合图的交叉数,所以“关于 指数”的 skein 树与“关于 指数”是同一件事。同章 三维流形不变量 词条对“输入表示决定复杂度”有同款讨论,可参照。
3. Temperley—Lieb 代数:从交叉到平滑
3.1 图示演算
Temperley—Lieb (TL) 代数 的元素是“没有交叉的配对图”的形式线性组合:图有 个顶部端点与 个底部端点,两两配对相连,连线互不交叉、不出边框。两个图的乘法是竖直堆叠:上图的底端接下图的顶端,然后压平。堆叠可能出现封闭的圈,规则是每个圈抹去并乘上环值 。
对每个位置 定义生成元 :把第 与第 个顶端用一条弧相连、第 与第 个底端用一条弧相连,其余股竖直穿过。直接数圈可以验证三条定义性关系:
逐条解释为什么成立:
- :把两个 堆叠,中间形成一个封闭小圈(内部那对弧上下相接围成),抹去后剩下的图恰好还是 ;所以 等于“ 乘上一个圈的值 ”。
- :三层的堆叠中,外侧两股各绕一个“挂钩”后回到原位,中间没有封闭圈( 与 位置的竖直股被 的弧包住但不闭合),压平后恰是 。这是三条中最需要画图核对的一条,留作习题 1 让读者把连通分支数出来。
- 远交换律: 与 触及的端点集 与 不相交(),堆叠次序显然无关。
这套演算与第 1.2 节的括号规则严丝合缝:TL 的图就是括号展开后的“无交叉中间态”,环值 就是脱圈规则里的 。
3.2 Braid 表示:把交叉展开成平滑
Kauffman 的关键观察:辫的每个交叉都可以就地展开成“平直型 + 转折型”两种平滑的线性组合,而平直型正是恒等图 、转折型正是 。于是定义表示候选
(约定提醒:哪种平滑配 取决于正交叉的画法约定,另一约定恰好交换 与 ,对应镜像链环;本词条固定此式。)由于 只依赖 附近两个位置,这个公式在任意 下同时定义。因为 是标量(中心元),多项式展开时可以随意合并同类项——后文的全部计算只用到这一点和三条 TL 关系。
它真的是表示吗? 这不是自动的:把群生成元的像随便写成代数元素,定义关系未必被保持。下一小节逐条检验,并发现一个漂亮的事实:三条辫关系全部成立,当且仅当环值恰好是 ——与括号的拓扑要求一模一样。
3.3 表示定理
Lemma 1(逆元). 对任意 , 与 互逆,当且仅当 。
证明。直接相乘(用 ):
右端等于 当且仅当系数 。Q.E.D.
Lemma 2(远交换律). 对 成立,对任意 。
证明。此时 (TL 远交换律),而 是 与 的线性组合,故两组乘法可交换。Q.E.D.
Lemma 3(辫关系). 成立,当且仅当 。
证明。记 、。展开左边( 是标量,八项分别为从每个因子取 或取生成元的所有组合;用 、):
(逐项来源:;两次“”各贡献 ;“”贡献 ;“”与“”各贡献 倍的 ;“”贡献 ;“”贡献 ,最后一步用了 。)右边由对称性()为
两式相减:
由于 (两个不同的图),辫关系成立当且仅当 ,两边乘 即 。Q.E.D.
把三个引理合起来:
Theorem 4(Kauffman). 取环值 时, 是群同态(把 映为上述可逆元素);并且此时括号在 R2 移动下的不变性(Lemma 1 的恒等式)与 R3 移动下的不变性(Lemma 3 的恒等式)逐条成立,R1 的亏损恰好是 因子,被 补偿。
同一个条件 出现了两次:一次保证代数(这是表示),一次保证拓扑(括号在 R2 下不变)。第 4.1 节将看到它第三次出现——保证量子(表示可取为酉)。这不是巧合:环值 同时是“圈抹去后的代数权重”、“拓扑自洽的参数”、“酉性的参数”这三个角色的唯一交点。
闭合后的读出。对 trace closure,把 TL 元素 ( 为配对图)闭合:每个图 闭合后是若干不相交圆圈之并,其括号是 的幂。由线性性,
其中 表示 的 trace closure。第 8 节的例子将把这条“数圈”规则一步步算出来。
4. 单位根与 path model:让表示变酉
4.1 酉性引理
量子线路只能实现酉算子,而 一般既不酉也不可逆(可逆性已被 Theorem 4 安排好)。现在把 实现为一个 Hilbert 空间上的算子,并问:何时 酉?
Lemma 5(酉性). 设 是 Hilbert 空间上的自伴算子()且满足 ,。取 。则
是酉算子,当且仅当 。
证明。由 有 ,故 。于是(用自伴性与 )
若 ,则 ,且同理 , 酉。反之,若 酉则 ,即 ;取 的表示(后文的 path model 即是)即得 。Q.E.D.
于是策略清楚了:在单位根上取 ,使 为实数,并把 实现为自伴算子。具体地,固定整数 (后文称 level),取
( 是 次单位根;分支选取的差别只影响整体相位与手性约定。)此时 。注意一个即将变得关键的数值事实:
4.2 图与 Perron—Frobenius 权重
取 个顶点排成一条线、相邻相连的路径图 (Dynkin 记号),顶点标号为 。它的邻接矩阵最大特征值(谱半径)是 ,对应的正特征向量分量为
这两句话可以直接验证:顶点 的两个邻居是 (端点 与 各只有一个邻居,下式同样成立),于是特征向量方程就是和差化积恒等式
即向量 被邻接矩阵乘后整体放大 倍。两个端点只有一个邻居,方程退化为 ,即 ——二倍角公式,同样成立(顶点 处由 化归为同一等式)。由 Perron—Frobenius 定理,正特征向量属于谱半径,故谱半径确为 。 时 、谱半径 ,与 一致。
就是原文文献中的 Perron—Frobenius 权重,它将出现在两个地方:构造使 自伴的内积(本节),以及 Markov trace 的终点权重(第 6 节)。
4.3 Path model:路径基与量子比特编码
Hilbert 空间。 path model 的态空间由路径标记:路径是图 上长度为 的顶点序列
即从顶点 出发、每步向相邻顶点移动、永不越界的行走。全体这样的路径张成一个有限维复向量空间,取路径为基(内积由 权重按文献标准方式修正后使 自伴——修正的具体形状不影响本词条后续内容,我们只用到“ 自伴、满足 TL 关系”这一事实,其存在性的构造性证明见 AJL 文)。
量子比特编码。 每一步只有“向上 ()”与“向下 ()”两种选择,正好编码为一个量子比特; 步用 个量子比特。全体 个比特串中,越界(触到 或 )的是非法路径,被排除在合法子空间之外;由于每一步的合法性只依赖“当前高度”,合法子空间可以用逐位的受控约束标记。对固定 ,合法路径数目按 量级增长(受谱半径控制,这是随机游走的标准渐近),严格小于 但对变化 仍是指数——这正是“经典上指数大、量子上 个比特就装得下”的典型处境。
生成元的作用。 只触及路径的第 与 步附近的局部形状:若这两步构成“先上后下”或“先下后上”的转折, 把转折改写(翻转或加权保留),矩阵元由 Perron—Frobenius 权重的比值 给出;若两步同向(单调穿过), 的作用为零。这一局部性的代数后果是:TL 的三条关系(第 3.1 节)逐条成立——它们的验证就是把比值化简,背后正是上节的和差化积恒等式。对本词条而言,关键结论是:
在单位根 ()处,配以 权重内积,每个 是作用在常数个相邻步寄存器上的酉算子。
(这是原文文献的表述:适当单位根上,path-model 内积使这些算子酉。第 4.1 节的 Lemma 5 说明了“为什么恰好在这个参数上”。)
Plat 态。 plat closure 的杯—帽结构对应一个特殊初态 :其支集集中在“路径在相邻高度间规则往返”的路径上(每个杯位强迫一次上—下往返)。这样的态由位置的奇偶性确定,是浅层电路可制备的乘积型态。
4.4 电路复杂度:逐项清单
把 braid word 编译成量子电路 ,按从右到左逐个执行 。总规模 的每个因子各有着落:
- (交叉数):每个字母是一个局部酉门,门数关于 线性。这是 AJL 算法的复杂度陈述:运行时间关于 多项式。
- (股数):门 作用在第 步寄存器(各 或常数个比特,视编码而定),但“第 个”需要寻址:用交换网络把目标寄存器移到一起、作用、再移回,每次门的布线代价 (或用寻址门 ),总计 量级。此外总寄存器数为 个步比特。
- (level):门的矩阵元是 型比值与 的幂,都是可经典预计算的 个常数(实现时取足够精度的定点近似);合法性判定需要知道“高度是否在 内”,涉及 个比特。
- 辅助:Hadamard test 需要额外 个比特(第 5.2 节)。
所以:电路大小 、量子比特数 。注意 既可作为输入(算法对 多项式),也可固定为常数(第 7 节的完备性恰恰取固定的 )。
5. Plat closure:估计一个矩阵元
5.1 公式的形状
对偶数股的 braid word ,其 plat closure 的 Jones 值在单位根 处可以写成
其中三个符号的含义:
- 是第 4 节编译出的电路;
- 是 cup-cap 初态(第 4.3 节);
- 是已知的、经典多项式时间可计算的归一化因子,它收集三类修正:辫的指数和(writhe)带来的 framing 因子 (plat 闭合时还包含股方向反转引起的符号修正,AJL 文中有显式公式)、闭合产生的圈数对应的 的幂、以及态 的归一化。
这句话的算法含义是:**量子部分只需要输出一个复数 的近似;所有拓扑换算都在经典侧完成。**剩下的任务是把“估计矩阵元”做快、做准。
5.2 Hadamard test 的完整推导
记 (,因为它是酉算子单位矢间的内积——Cauchy—Schwarz)。Hadamard test 用一个辅助比特(ancilla)与受控- 把 的实部与虚部变成辅助比特的测量概率。推导如下。
第一步,制备辅助比特 (一个 Hadamard 门作用于 即得),目标寄存器制备 ,联合态为
第二步,执行受控-(辅助比特为 时作用 ):
(这里用了 、;展开出的四项中,辅助比特指标相同的两项因 保留,指标不同的两项因 消去。)保留的两项是
即辅助比特的 期望恰是 ,从而测得 的概率为 。
第四步,测虚部:对辅助比特测量 。用 、(即 、):
(末步用了 。)实践上 基测量等价于在受控门后给辅助比特加一个相位门再做 基测量。于是
两个概率各用一组重复测量即可还原复数 。
5.3 采样复杂度:Hoeffding 与振幅估计
单次 Hadamard test 给出一个 Bernoulli 随机变量(参数 )。用 次独立重复的经验频率 估计 ,Hoeffding 不等式给出
要常数置信度(如 )需 ,即 ,也就是 。取 (这样 ),得实部到加性 需要
次电路运行;虚部同阶。每次运行执行一遍 电路(代价 ),所以总复杂度
振幅估计 (amplitude estimation) 把精度因子开根号: 次相干调用给出同样加性精度(已证;其机制是把“成功概率 ”编成 Grover 型迭代的旋转角 ,,再用相位估计读角)。代价有二:需要实现受控- 的幂 (把电路重复 次,几何级数求和后总门数仍是 );且不含测量地相干复用中间态。于是
两种读出的取舍就是“ 次含测量运行”与“ 次相干运行”之间的取舍;对 两者都是多项式,对指数小的 两者都指数——精度的多项式性是硬边界。
5.4 误差尺度:加性近似买到了什么
把矩阵元的估计 ()代回第 5.1 节的公式,Jones 值的误差是
这个不等式是本词条所有“保留条款”的技术根源,值得把每一层说透:
- ** 可以关于 指数大。**它包含 的幂(每多一个闭合圈贡献一个 , 的模是常数但幂次随 增长)与 framing 因子。于是“矩阵元的 加性近似”换算到 Jones 值是“尺度 上的加性近似”——若 指数大,绝对误差 对典型的 Jones 值可能毫无信息量。
- **它不是相对误差。**相对误差要求 ;对 的链环,由 可知这要求 ,而典型矩阵元振幅约 、随 指数小——相对逼近一般需要指数采样。
- **它更不是精确系数。**Jones 是 Laurent 多项式,精确读出某个系数需要把误差压到多项式间距以内;一般求值点上这保持 -hard(第 1.3 节)。
一句话总结:**量子算法给出的是“相对自然归一化 的加性近似”,这是一个明确定义、有内容、但范围有限的承诺。**夸大它(说成“量子计算机算出了 Jones 多项式”)与贬低它(说成“毫无意义”)都不对——第 7 节将看到,恰恰是这个受限的近似问题,对量子计算本身是完备的。
6. Trace closure:Markov trace 与 DQC1
6.1 为什么 trace closure 需要加权迹
回忆 Markov 定理:两个辫 、 的 trace closure 同痕,当且仅当它们可通过两类移动互变:共轭(,同群内)与稳定化(,加一股打一次交叉)。要在表示论一侧得到链环的泛函,必须找到一个对这两类移动不变(差已知因子)的线性泛函:
- 共轭不变 它必须是迹(trace 性质 );
- 稳定化不变 它在“多闭合一股”时的行为被固定,这唯一确定了各股指标上的权重——计算的结果正是 Perron—Frobenius 权重。
在 path model 的路径基上,这个 Markov trace 写成对路径终点的加权归一化迹:
求和跑遍合法路径 ,权重 只依赖终点 所在顶点(),归一化成概率分布。Jones 值等于这个迹乘上一个已知的经典归一化(包含 writhe 因子与 的幂)。直觉:trace closure 把“第 股接回第 股”翻译成“对股指标求和”,而路径的终点记录的正是“这 股组合成的拓扑荷落在哪个 sector”——加权求和就是按 sector 配平。
6.2 两条实现路线与一个参数平衡
路线一:提纯 (purification)。 只依赖终点寄存器,故加权叠加态
是一个可高效制备的态:权重只依赖终点,故 是一个经典高效可采样的分布——先按 ( 为终点在 的合法路径数,可动态规划经典计算)抽终点 、再均匀抽一条到 的路径——相干版沿同一结构制备即可。在它上面做一次 Hadamard test(第 5.2 节,把 换成 ),期望值就是
中间一步用了路径基的正交归一性 。于是一次 Hadamard test 直接采样 Markov trace 本身,第 5.3 节的采样复杂度原样适用:,或振幅估计的 。
**路线二:经典抽样 + 逐个 Hadamard test。**也可以先用经典随机数按 抽路径 ,再对每个 做 Hadamard test 估计 ,最后平均。这不需要制备 ,代价是误差来源变成两个。
**两阶段误差的平衡(展示求解过程)。**设抽出 个路径、每个做 次 Hadamard 测量。估计量的误差有两个独立来源:(i) 对 的 Monte-Carlo 抽样,贡献标准差 ,其中 是各 sector 对角元在 下的方差;(ii) 每个路径的 Bernoulli 测量噪声,贡献 ( 为常数,由第 5.3 节 Hoeffding 的换算决定)。总误差近似为两者平方和的平方根,为简单起见按最坏情形相加:
把 代入,得 。对 求导置零:
代回得 ,即最优总误差 ,且最优分配让两项误差同阶(各占一半)——这就是“取参数使两项同阶”的平衡原则的一个具体实例:样本该在“多抽几个路径”与“每个路径测准一点”之间按 分配。若 很小(各 sector 对角元几乎相同),就少抽路径、每个多测;反之亦然。注意两个阶段只能换来 :预算被劈成两半,每个阶段只分到 次资源、各贡献 的误差。对比之下,路线一(提纯)把两个误差源合并成一个——经典混合的方差被相干地吸收进 ,只剩量子测量噪声一项,误差回到 (即 采样)。这就是提纯优于两阶段抽样的定量理由。
6.3 DQC1:归一化迹的母语
DQC1(deterministic quantum computation with one clean qubit,单干净比特量子计算)是 Knill 与 Laflamme 于 1998 年提出的计算模型: 个完全混合的量子比特加一个干净比特,允许对全体做任意酉操作后测量干净比特。它介于经典与 BQP 之间(,一般认为严格包含于 BQP),最著名的展示是它能在多项式时间内估计任意酉算子的归一化迹 ——推导只需三行。
初始态(干净比特 ,其余最大混合):
对干净比特做 ,再执行受控-,然后只看干净比特的约化态。 之后 。受控- 的共轭作用逐块计算:对角块 与 分别映到 与 ( 酉);非对角块 、。于是
对混合寄存器求迹(,),记 :
(最后两个等号: 等。)整个过程没有制备任何纯的 比特态——归一化迹本来就是混合态可观测量。
这正是 Markov trace 的形状:一个“对角矩阵元的加权平均”。Shor 与 Jordan 由此证明,某些 Jones trace closure 近似(在论文指定的单位根、closure 与加性误差尺度下)是 DQC1-complete 的:它既能在 DQC1 中完成(上述电路加上对路径 sector 的经典随机化,把权重 并入抽步骤),又足以模拟 DQC1 中的任意计算。注意保留条款:是“某些”——完备性结论绑定于特定的根与误差归一化,不是说所有 Jones 迹都同样难;也不是说 DQC1 能做 plat closure 的完备版本(那属于 BQP,第 7 节)。
7. 为什么 BQP-hard:拓扑量子计算视角
第 4、5 节给出了“plat closure 的加性近似 BQP”。反方向更强:同一个近似问题对 BQP 是完备的。证明思路来自拓扑量子计算,两个要件:
要件一:编织就是计算(表示的稠密性)。在 维 TQFT 描述的任意子系统中, 个任意子的编织(braid)作用于系统的融合空间,给出正是 braid 群的(射影)酉表示——与本词条的 同族。对合适的 level(例如 ,对应 Fibonacci 型任意子理论),该表示的像在相应的酉群中稠密:任意目标酉门都能被某个编织词逼近。这是 Freedman—Larsen—Wang 一线的稠密性定理;配合 Solovay—Kitaev 型编译(把“稠密”升级为“多项式长度词逼近到 精度”),得到:任何量子电路都可以被经典多项式时间编译成一个多项式长度的 braid word 。(AJL 文的表述即:其算法在五次单位根 处解决的是 BQP-complete 问题,依据正是 Freedman 等人的结果。)
要件二:电路的接受振幅被编码为辫的矩阵元。编译可以安排得使任意电路 的接受振幅(或接受概率)与 的表示矩阵元相差一个已知的归一化因子:
直觉: 扮演电路的 初态,编织 扮演电路演化,读出矩阵元扮演测量振幅——拓扑系统“物理上就是”一台量子计算机,其振幅天然是拓扑不变量的值。
**拼起来的归约。**假设存在(经典或量子)算法 能在 AJL 指定的归一化尺度内加性逼近 。给定任意 BQP 电路 :经典多项式时间编译出 ,运行 ,乘上已知尺度读出 的近似——精度取到足以区分“接受概率 ”与“”即可(中间留有常数间隙,编译误差可以做得远小于间隙)。于是 解决任意 BQP 问题,即该近似问题 BQP-hard;结合第 5 节的算法,它是 BQP-complete。
保留条款(务必逐条记牢)。被证明完备的是以下三者同时成立的问题版本:
- 特定单位根: 取使表示稠密的那一档(如 ),不是“对所有根”;
- 特定 closure:plat closure 的矩阵元形式(trace closure 对应 DQC1 侧);
- 特定加性误差尺度:误差相对于指定的(可指数大的)归一化 声明,这是问题定义的一部分,不是算法偷懒。
因此该结论不表示“Jones 多项式被量子计算机攻克了”:精确值与相对误差在第 5.4 节的意义下仍然困难(一般点 -hard),也不表示所有求值点同等困难(特殊点有多项式经典算法,第 1.3 节)。同章 三维流形不变量 词条对 Turaev—Viro 不变量的 BQP-complete 结论有平行的讨论,可对照阅读。
8. 手算例子: 上的 Hopf 链环与三叶结
本节把前文的机制在一个能完全手算的设定里走一遍:两股辫、。所有代数只用三条 TL 关系与“数圈”规则,量子算法逼近的正是我们一步步算出的这些数。
**设定。**取 ,故 ,(则 )。环值
注意 恰是 的谱半径(第 4.2 节),且 、。两股()的合法路径只有
所以 path model 的态空间是二维的——本例中量子电路 就是一个 酉矩阵(作用在这两个基矢张成的空间上)。
第一步:验证酉性条件。,故 Lemma 5 的条件 成立, 酉;Theorem 4 同时保证它是 的表示( 只有一个生成元,辫关系自动成立,但引理 1 的可逆性条件已核对)。把 这个关键数值实际算出来,看它确实落在单位圆上(提出公因子 ,再用 与 ):
而 本身模长为 : 的两个“候选本征方向”( 取本征值 或 )上的本征值 与 都在单位圆上——这正是 Lemma 5 在本例中的数值体现。
第二步:单个卷(校准坐标)。。把两项分别 trace-close:恒等图(两条竖线)闭合后是两个不相交圆圈(每条竖线与自己的闭合弧成一圈),括号 ;(上弧接上弧、下弧接下弧)闭合后是一个圆圈,括号 。所以
这正是第 1.2 节预言的单卷行为(),与“补偿因子 ”的机制吻合:带因子后 ,而 是平凡结,。校准通过。
**第三步:Hopf 链环()。**用 展开( 是标量):
闭合读出(恒等图 、 图 ):
(第三个等号代入 ,末步展开平方。)与文献中 Hopf 链环的括号值一致。换到 Jones:闭合图有 个正交叉,
在根值处:,故 ;而
顺带一提,二分量平凡链环的值是 (环值亲自现身),而 Hopf 的镜像给出 :这个根值足以区分 Hopf 与平凡链环、以及两种手性。
**第四步:三叶结()。**同法展开(二项式型求和,):
闭合:
逐项代入 并展开:;;。求和( 的系数 , 的系数 ):
与文献中三叶结的括号值完全一致。Jones 值(,):
(每一步只用 型的指数加法。)在根值处:
镜像三叶结(把 与 互换的约定)在同一根处取值 。
第五步:算法视角复盘。量子算法在这个例子上做什么?它不展开上面的代数,而是:
- 编译: 作为二维空间上的酉矩阵 由三个局部门(每个对应一个 )实现;
- 制备 cup 态 (或提纯态 ,权重 );
- Hadamard test 测 与虚部, 次运行得加性 ;
- 乘上经典可算的 (含 、圈数幂与态归一化),得到 的 加性近似。
我们上面用纯经典的手算交叉验证了第 4 步的目标值。规模增大时,手算的项数指数增长( 的 TL 展开有指数多项),而量子电路仍是 个门的乘积——这就是“指数经典、多项式量子”差距的具体形状。同时注意本例的 是常数阶,加性近似足以区分 (相差 );对一般输入,第 5.4 节的告诫仍然全部有效。
9. 更远的一步:HOMFLY 与 Khovanov 同调
9.1 Hecke 代数与 HOMFLY
Jones 多项式有一个双变量推广 HOMFLY 多项式 (以发现者姓名首字母命名),它同时单变量化为 Jones 与 Alexander 多项式。表示论上,把 Temperley—Lieb 代数换成 Hecke 代数(对称群代数的 -形变;TL 代数是它的一个商),同一套“辫 代数表示 矩阵元/迹”的机器原样运转:HOMFLY 在适当参数点上同样可以被量子算法做加性近似(Zoo 93 一线的工作)。再往外,quantum double(Drinfeld double)构造给出另一族链环不变量(Zoo 174),其算法面貌与本词条 三维流形不变量 词条的 Turaev—Viro/Reshetikhin—Turaev 状态和同源。这一族推广的共同结构是:“结的拓扑 某个代数的表示论 酉矩阵元”三层字典,本词条第 3—5 节建立的正是其中最简单的一层。
9.2 范畴化:Khovanov 同调
2000 年,Khovanov 对 Jones 多项式做了范畴化 (categorification):构造一族分次阿贝尔群——Khovanov 同调 ——使其分次 Euler 特征标恢复 Jones 多项式(差一个已知的移位与线性替换),而同调群本身携带严格更多的信息。构造的形状:
- 链群:对结图的每个交叉选一种平滑(共 种选择,称为 Kauffman 状态),每个状态按圈数与平滑类型获得双分次 ,生成分次自由阿贝尔群。链空间维数 是 量级——指数大。
- 边缘算子:(合并/分裂圈的标准映射配上符号安排),满足 (符号安排的相容性需要验证,这是构造的技术核心)。
- 同调:
闭链模边缘链。分次 Euler 特征标 (把维数放进生成函数)恢复 Jones。
信息量的差距是实质性的:存在 Jones 多项式相同而 Khovanov 同调不同的链环对;Khovanov 同调能识别平凡结(同调的秩达到平凡结的最小值就蕴含链环平凡——这是已被证明的定理;而 Jones 多项式是否也能识别平凡结,即 是否蕴含 平凡,至今是未解决问题)。因此“算 Khovanov 同调”是严格更强的算法问题。复杂性一侧,精确计算当然至少和 Jones 一样难;2025 年的工作(Zoo 510—511)进一步把秩的加性近似按精度分了档:越来越精确的加性近似分别是 DQC1-hard、BQP-hard、-hard——精度档位每提高一级,难度上升一个世界,这与第 5.4 节“加性 vs 相对 vs 精确”的阶梯精确呼应。
9.3 Hodge Laplacian 方案与它的条件
2025 年的量子算法(arXiv:2501.12378)把同调秩的计算化为一个谱问题。核心是Hodge Laplacian
其中 是边缘算子关于链空间标准内积的伴随。
Lemma 6(调和代表). ,且其维数等于同调维数(Betti 数 / Khovanov 秩)。
证明。设 。计算(用 自伴、 与 半正定):
(第二个等号:,这是伴随的定义;第三项同理。)范数平方非负且和为零,故 (闭)且 (余闭)。反之闭且余闭的 显然被 消灭。于是 = 闭且余闭的链(调和链)。每个同调类 有唯一的调和代表(Hodge 分解 ),映射 是同调群到调和链空间的同构,故维数相等。Q.E.D.
算法的形状。于是“算 ”变成“算 的零空间维数”。由于链空间维数指数大,不能直接对角化;方案是估计 的近零谱质量——零空间占全空间的比例:
- Gibbs / 相位估计提取:制备 Hodge Laplacian 的 Gibbs 态 (或对 做相位估计),测量“能量近零”的概率,它近似 与全维数之比。Gibbs 态制备与相位估计分别关联(Gibbs 态制备)与(相位估计)的工具,谱隙决定所需分辨精度。
- 预热化 (pre-thermalization):若 Betti 数远小于链空间维数(同调极小、链群巨大——一般情形正是如此),近零谱质量指数小,直接采样测不到。预热化步骤先把态集中到近零/低能子空间附近再做 Gibbs/相位估计提取,专门克服“秩 维数”的失衡——这是该工作相对先前量子同调算法的关键改进。
**适用条件(务必保留的条款)。**论文声明的高效性依赖两个承诺,并有数值与解析证据支持、但不是无条件定理:
- 可热化: 的 Gibbs 态能在多项式时间内制备(典型 Gibbs 制备的障碍——谱隙、温度、与环境的耦合设计——这里全部内嵌于“thermalize”一词,见 Gibbs 态词条的讨论);
- 谱隙:零空间(严格说零特征值簇)与其余谱之间有 的间隙,否则相位估计无法在多项式时间内把“近零”与“非零”分开。
没有这些条件时:链空间指数大、同调秩可能极小,近零质量被指数压制,任何基于谱估计的方案都需要指数资源——不能宣称普适的指数加速。与第 5.4 节 Jones 的处境对照:一边是“归一化尺度的加性近似”,一边是“热化与谱隙承诺下的秩估计”,两者都以“精确说明自己买到了什么”为底线,这正是本词条导言的原则在本词条的两次落地。
10. 小结与习题
小结。
- 结图的 Kauffman 括号由三条规则定义(归一化、脱圈 、交叉展开),经 writhe 补偿因子 后给出 Jones 多项式;R2/R3 不变性恰好化为 Temperley—Lieb 代数的恒等式。
- 链环经 Alexander 定理化成 braid word(长度 、股数 ); 是 braid 群表示当且仅当 (Lemma 1—3 的计算),这也是括号的拓扑自洽条件与酉性条件(Lemma 5)——同一个等式出现三次。
- 单位根 、 处, 与路径图 的谱半径重合,Perron—Frobenius 权重 (和差化积恒等式)构造出使 自伴的路径表示;braid word 编译为 的局部酉电路(因子 的来源逐项清点于第 4.4 节)。
- Plat closure = 矩阵元:Hadamard test 以 采样(Hoeffding)或振幅估计 给出 的加性 估计;Jones 值的误差是 ——加性、非相对、非精确,且 可指数大(第 5.4 节)。
- Trace closure = Markov trace:终点权重 的加权归一化迹(Markov 定理迫使的形状);可经提纯态一次 Hadamard test 估计,或两阶段抽样(第 6.2 节给出 的平衡解);归一化迹是 DQC1 的母语,某些 trace-closure 近似 DQC1-complete(Shor—Jordan)。
- plat closure 的加性近似在适当固定的 (如 ,Fibonacci 型稠密表示)下 BQP-complete(AJL + Freedman—Larsen—Wang + Solovay—Kitaev 编译);完备性绑定特定根、closure 与误差尺度。
- Khovanov 同调范畴化 Jones(分次 Euler 特征标恢复之、信息严格更多);2025 年方案把秩计算化为 Hodge Laplacian 的近零谱质量( 同调,Lemma 6),以预热化克服“秩 维数”,高效性依赖多项式热化与谱隙承诺;秩的加性近似按精度分为 DQC1-hard / BQP-hard / -hard 三档。
习题。
- **(TL 图示演算)**画图验证第 3.1 节的三条 TL 关系:(a) 对 ,数出堆叠中封闭圈的个数并说明为何恰好是一个;(b) 对 ,追踪三条外侧股的连通情况,说明中间为何没有产生封闭圈;(c) 用 (a)(b) 复核第 3.3 节 Lemma 3 展开中 与 的每次使用。
- **(计算)**沿第 8 节的方法计算 (提示:,逐项写出)与 ,并求 trace closure( 环面链环)的 Jones 值在 处的显式复数。再用 与 skein 关系 核对你的答案。
- **(采样置信度)**第 5.3 节用 Hoeffding 不等式得到常数置信度下 。若要求置信度 ,写出 关于 与 的完整表达式,并解释为什么对 的依赖只是对数级的(联系振幅估计中“失败重试”的标准论证)。
- **(两阶段平衡)**在第 6.2 节的两阶段方案中,设 、、总预算 次电路运行。求最优的 与相应的总误差;若把预算提高为 ,误差改善多少?由此说明“把样本在两个阶段间按 分配、使两项误差同阶”的含义。
- **(误差尺度)**设某链环的 、而 。(a) 要让加性估计足以判断 的符号, 需多小?采样复杂度随 如何增长?(b) 由此解释为什么“加性近似 BQP-complete”与“相对误差近似可能需要指数资源”并不矛盾,并说明这与第 5.4 节、第 9.3 节两处保留条款的关系。
- **(概念辨析)**分别给出以下三个量能区分、而前一个不能区分的链环例子或论证方向:(a) Jones 多项式 vs Khovanov 分次 Euler 特征标(它们何时给出相同信息?);(b) Euler 特征标 vs Khovanov 同调群(提示:同一组分次维数可以给出相同 Euler 特征标而群结构不同;再想想“识别平凡结”的性质落在哪一侧);(c) 结合第 9.2 节的三档硬度阶梯,说明“把秩估计到越来越准”为什么在每个档位上落进不同的复杂度类。
参考文献与 Zoo 覆盖
- Zoo 2、4、41、42:Aharonov—Jones—Landau 与 TQFT universality/BQP-hardness。
- Zoo 83、93、174:DQC1 trace closure、HOMFLY 与 quantum-double link invariants。
- Zoo 3:Tutte/Potts 推广。
- Zoo 510—511:cohomology complexity 与 Khovanov homology 量子算法。
参考资料
- 本词条整理自《量子计算算法教程》原文:ch13-topology-statistical-physics/knot-invariants
- 内容遵循 CC BY-NC-SA 4.0 许可协议