Lattice Quantum Chromodynamics for Quantum Simulations

技术摘要：用于量子模拟的晶格量子色动力学

问题陈述
量子色动力学（QCD）是强相互作用的基础理论，然而其非微扰机制——对于理解强子结构和高密度物质至关重要——在计算上仍具有挑战性。标准的晶格 QCD 方法依赖于欧几里得路径积分和重要性采样，但在存在“符号问题”（sign problem）的情况下（例如在实时动力学或包含 $\theta$ 角时），这些方法会失效。虽然哈密顿量表述为实时模拟提供了一条路径，但随着晶格体积增加，希尔伯特空间的指数级增长使得大规模系统的经典模拟变得难以实现。因此，迫切需要一种将晶格规范场论映射到量子计算机上的框架，特别是要解决规范自由度和费米子自由度的编码、物理（规范不变）态的构建，以及高效时间演化算法的开发问题。

方法论
作者开发了一个用于 $SU(N_c)$ 晶格规范场论（含夸克）量子模拟的综合框架，重点研究二维和三维空间中的 $N_c=3$（QCD）。该方法建立在以下支柱之上：

1. 哈密顿量表述与离散化：
   本工作采用了 Kogut-Susskind 哈密顿量表述。作者同时考虑了交错费米子（staggered fermions）和威尔逊费米子（Wilson fermions），以解决费米子倍增问题。该理论在连续时间的立方空间晶格上进行离散化。哈密顿量被分解为电场项（$H_E$）、磁场项（$H_B$）、$\theta$ 项（$H_\theta$）、动能项（$H_{kin}$）和质量项（$H_{mass}$）。$\theta$ 项作为与拓扑电荷密度的耦合项被引入，这在哈密顿量中引入了复系数。

2. 表示基底与物理希尔伯特空间：
   一个核心贡献是在**表示基底（representation basis）**中构建物理希尔伯特空间。不同于直接将群元素映射到量子比特，规范自由度（位于链路/link上）和费米子自由度（位于格点/site上）是使用 $SU(N_c)$ 的不可约表示（irreps）进行编码的。

   • 链路希尔伯特空间： 链路由状态 $|(R, r_R, r_L)\rangle$ 描述，其中 $R$ 是一个不可约表示，$r_{R,L}$ 是左右半链路的基向量。该基底使二次卡西米尔算符（quadratic Casimir operators）对角化。
   • 格点希尔伯特空间： 费米子通过占据数基底进行处理，随后利用 Clebsch-Gordan 系数（CGCs）将其转化为表示基底。这使得粒子和反粒子态可以根据其在 $SU(N_c)$ 下的变换性质进行分类。
   • 规范不变性： 物理态被构建为由所有在晶格格点处相遇的不可约表示（链路与费米子）通过张量积形成的单态（singlets/trivial representations）。作者明确定义了构造这些单态所需的“格点顺序”（site order）约定，以确保满足高斯定律约束。

3. 截断与编码：
   由于 $SU(N_c)$ 的希尔伯特空间是无限维的，作者实施了截断方案。这包括对规范链路中最高权重（$n_1$）或二次卡西米尔特征值（$C_2$）的截断，以及对费米子颜色占据数的限制。物理态通过**局部编码（local encoding）**方案编码进量子比特，其中量子比特代表特定的不可约表示和单态多重性指标。这种方法被认为在处理大型希尔伯特空间时比直接编码更具量子比特效率。

4. 量子模拟算法：
   作者提出了用于时间演化的 Trotterization 算法。时间演化算符 $e^{-iHt}$ 被近似为对应于分解后的哈密顿量项的酉算符乘积。

   • 门构造： 普拉克特（plaquette）算符和链路算符的演化被映射为 Givens 旋转（双能级酉算符）和相位移动。作者详细说明了用于这些操作的量子电路构造，包括处理费米子反对易关系所需的 Jordan-Wigner 链（相位翻转）。
   • 矩阵元： 文献提供了物理基底下所有哈密顿量算符的显式公式，这些公式表现为涉及 CGCs 的“格点因子”乘积。这些公式对于计算所需量子门的旋转角度和相位至关重要。

主要贡献

• **含 $\theta$ 角的 $SU(3)$框架：** 本文展示了在哈密顿量方法下，在量子计算机上模拟三维空间中带有$\theta$ 角的晶格 QCD 的首次构架与模拟。
• 表示基底形式体系： 它提供了一种使用不可约表示和单态 CGCs 构建物理希尔伯特空间的严谨构造方法，为规模化替代群元素基底提供了途径。
• 算法细节： 该工作详细阐述了用于 Trotter 化演化的特定量子电路，包括将复杂的哈密顿量项（如 $\theta$ 项和 Wilson 动能项）分解为基本门的过程。
• 矩阵元计算： 作者推导并列表了普拉克特、$\theta$ 和费米子动能算符的矩阵元，为实现模拟电路提供了必要的数据。

结果
作者在小型晶格（最多为 $2\times2$ 周期性晶格和 $2\times2\times2$ 开边界立方体）上进行了无噪声模拟，使用了高达 32 个量子比特的经典模拟器（Qiskit Aer 和 CUDA-Q）。模拟表明：

• 希尔伯特空间截断效应： 分析了不同的规范和费米子截断如何影响观测物理量，如平均电场（$\langle \Pi^2 \rangle$）和颜色占据数（$\langle N_q \rangle$）。结果表明，强耦合抑制了高占据态，使得适度的截断对于交错费米子是有效的，而威尔逊费米子则表现出对截断更高的敏感性。
• 参数依赖性： 研究了观测量随费米子质量和 Wilson 参数 $r$ 变化的关系。研究发现，较大的质量会抑制向费米子丰富态的振荡。
• $\theta$ 角动力学： 改变 $\theta$ 角的模拟揭示了强耦合机制下的非周期行为，且正负 $\theta$ 具有截然不同的动力学特征。
• 强子态与弦动力学： 作者演示了“带电派子（charged pion）”态的制备以及弦态（包括回路弦和费米子弦）的演化。他们观察到了弦断裂（string breaking）动力学，即通量管衰变为费米子-反费米子对，并追踪了弦断裂配置与未断裂配置的概率。
• 重子化学势： 通过在单个立方体上进行精确对角化，作者计算了添加一个重子所需的能量代价，这作为研究高密度物质的原理验证。

意义与主张
本文旨在为基础物理学的实际量子模拟开发奠定基础。作者声称，其工作为设计哈密顿量表述下的晶格 QCD 量子模拟提供了必要的“要素”。

• 概念验证： 其主要意义在于证明了复杂的 $SU(3)$晶格规范场论（包括$\theta$ 角和费米子物质）可以在量子硬件（此处为模拟硬件）上进行编码和演化。
• 扩展路径： 局部编码和表示基底被呈现为迈向更大规模模拟的可行路径，相比全局编码可能需要更少的量子比特。
• 未来展望： 作者谦虚地指出，目前的结果具有探索性，由于晶格尺寸较小、处于强耦合机制以及使用了激进的截断，尚未对真实的实时 QCD 做出严肃的物理结论。他们指出，若要达到连续极限并模拟具有物理意义的机制，需要改进的哈密顿量（反项/counterterms）、更好的态制备算法以及容错硬件。这项工作被定位为未来模拟致密物质（如中子星）和轴子-QCD 动力学的垫脚石。