Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}_2$ from lattice QCD at physical quark masses

该论文利用物理夸克质量下的格点量子色动力学计算，在大型动量有效理论框架下首次提取了质子同位旋矢量螺旋度部分子分布函数（PDF）及 twist-3 矩 $\tilde{d}_2$，并通过高阶微扰匹配与重求和技术确定了 $\overline{\rm MS}$ 方案下 $x\in[0.25,0.75]$ 区域的 PDF 及其相关物理量。

要点

这篇论文就像是一次**“质子内部结构的超级 X 光扫描”**，科学家们试图用超级计算机（格点量子色动力学）来解开一个困扰物理学界几十年的谜题：质子的自旋（Spin）到底是从哪里来的？

为了让你更容易理解，我们可以把质子想象成一个高速旋转的“宇宙陀螺”，而构成它的夸克（Quarks）和胶子（Gluons）就是陀螺里的“小零件”。

以下是这篇论文的通俗解读：

1. 核心谜题：陀螺的旋转力从何而来？

早在几十年前，科学家发现质子内部夸克的自旋加起来，只占质子总自旋的很小一部分（大约 30%）。剩下的 70% 去哪了？是胶子在旋转？还是夸克在内部乱跑产生的“轨道运动”？这就是著名的**“质子自旋之谜”**。

这篇论文的目标，就是要把质子内部那些**“带着自旋方向的夸克”（也就是所谓的“螺旋度部分子分布函数”，PDF）画出来。这就好比我们要给质子内部画一张“自旋地图”**，看看在质子内部不同位置（用 $x$ 表示，代表夸克携带动量的比例），夸克是顺着质子旋转方向转，还是反着转。

2. 科学家的新工具：LaMET（大动量有效理论）

在传统的实验室里，我们无法直接“看到”质子内部，因为夸克被强力锁在里面，而且我们只能在“静止”的参考系里做实验，但质子内部的东西是沿着“光锥”（光速方向）运动的。

这就好比你想看清一辆飞驰的赛车内部，但你只能站在路边看。

• 传统方法：很难直接看清。
• 新方法（LaMET）：科学家发明了一种“时间机器”般的数学技巧。他们让质子在我们的超级计算机里**“加速”**到接近光速（虽然实际上是在格点上模拟高动量状态）。
• 比喻：这就好比你把那个静止的陀螺突然加速到超音速。在高速状态下，原本复杂的内部运动变得简单了，我们可以用一套“翻译器”（微扰匹配），把我们在计算机里看到的“高速快照”（准分布函数），翻译回我们想要的“真实地图”（光锥分布函数）。

3. 这次做了什么？（两大成果）

成果一：绘制了更清晰的“自旋地图”

科学家在物理质量的夸克（最真实的模拟环境）上，计算了质子内部夸克的自旋分布。

• 过程：他们把质子加速到不同的速度（0, 0.25, 1.02, 1.53 GeV），就像给陀螺拍不同角度的高速照片。
• 处理：照片里有很多噪点（数学上的发散），他们用了高级的“滤镜”（重整化方案，特别是混合方案和重求和技巧）去除了噪点。
• 结果：他们成功画出了质子中间区域（$x$ 在 0.25 到 0.75 之间）的自旋分布图。
• 发现：他们发现，在这个中间区域，夸克顺着质子旋转方向的比例，比之前全球科学家通过实验数据拟合出来的结果要更高一些。这就像大家原本以为陀螺中间是“空”的，结果发现其实转得挺欢。

成果二：首次测量了“扭曲力”（Twist-3 Moment $\tilde{d}_2$）

这是论文最亮眼的创新点。

• 背景：除了自旋，夸克在质子内部还会受到一种复杂的力，就像在拥挤的舞池里跳舞，不仅要转圈，还要被周围的人（胶子）推来推去。这种**“夸克 - 胶子之间的相互作用力”**被称为“色洛伦兹力”。
• 比喻：想象质子是一个拥挤的舞池。夸克是舞者，胶子是周围的推挤者。$\tilde{d}_2$ 这个数值，衡量的就是**“舞者平均受到的侧向推力有多大”**。
• 突破：以前很难直接算出这个力，因为数学上太复杂（会有无穷大的项）。这次，科学家利用一种新的数学技巧（算符乘积展开 OPE），第一次在 MS 方案（一种标准的数学语言）下，直接从短距离的矩阵元素中提取了这个力。
• 结果：算出来的数值非常小（接近于零）。这意味着，在质子内部，这种复杂的“侧向推力”在平均意义上被强烈抑制了。就像舞池里的推挤力虽然存在，但大家跳得很有序，平均下来并没有把谁推倒。

4. 为什么这很重要？

• 验证理论：这是第一次从第一性原理（纯数学计算，不依赖实验拟合）出发，直接算出了这个复杂的力，验证了我们对强相互作用的理解。
• 填补空白：之前的实验很难精确测量这个力，现在有了理论上的“标尺”，未来的实验（如电子离子对撞机 EIC）可以更有针对性地去验证。
• 技术突破：他们展示了如何在物理质量的格点上，同时处理“自旋”和“复杂力”的计算，为未来更精确的质子结构研究铺平了道路。

总结

这篇论文就像是用超级计算机给质子做了一次**“高精度 CT 扫描”**。

1. 他们不仅看清了质子内部夸克**“怎么转”**（自旋分布），还发现中间区域的旋转比预想的更活跃。
2. 他们第一次算出了夸克在内部受到的**“侧向推力”**（$\tilde{d}_2$），发现这个力其实很小，说明质子内部虽然拥挤，但秩序井然。

这让我们离彻底解开“质子自旋之谜”又近了一大步，也展示了人类利用数学和超级计算探索微观世界的能力。

技术摘要

这是一份关于论文《Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}_2$ from lattice QCD at physical quark masses》（物理夸克质量下格点 QCD 计算的质子同位旋矢量螺旋度部分子分布函数及 twist-3 矩 $\tilde{d}_2$）的详细技术总结。

1. 研究背景与问题 (Problem)

• 质子自旋之谜与螺旋度 PDF： 质子自旋中夸克自旋的贡献（由螺旋度部分子分布函数 $g_1(x)$ 描述）是理解质子结构的关键。虽然实验（如 RHIC, COMPASS）提供了约束，但在小 $x$ 区域和高精度下仍存在不确定性。格点 QCD 作为第一性原理计算方法，对于验证和补充实验数据至关重要。
• 计算挑战：
  • 物理夸克质量与连续极限： 许多之前的格点计算要么在重夸克质量下进行，要么未外推至连续极限。同时实现物理夸克质量和精细晶格间距极具挑战性。
  • 非微扰定义： 螺旋度 PDF 定义在光锥上，无法直接在欧几里得时空的格点上计算。
  • 高阶修正与重整化： 需要次次领头阶（NNLO）的匹配精度以及处理线性发散和重整化子（renormalon）模糊性，以控制理论误差。
  • Twist-3 矩 $\tilde{d}_2$： 这是一个描述夸克 - 胶子关联（平均色洛伦兹力）的 twist-3 矩。由于 twist-3 算符混合了低维算符且存在 $1/a$ 发散，直接计算非常困难。此前缺乏在 $\overline{\text{MS}}$ 方案下直接从非局域关联函数提取 $\tilde{d}_2$ 的格点结果。

2. 方法论 (Methodology)

该研究采用了大动量有效理论 (LaMET) 框架，结合短距离因子化 (SDF) 方法，在物理夸克质量下进行了计算。

• 格点设置：

  • 系综： 使用 HotQCD 合作组生成的 $N_f=2+1$ 味 HISQ 作用量系综，物理夸克质量（$m_\pi \approx 140$ MeV）。
  • 参数： 晶格间距 $a = 0.076$ fm，体积 $64^3 \times 64$，共 350 个构型。
  • 动量： 质子态被Boost 到 $P_z = \{0, 0.25, 1.02, 1.53\}$ GeV。
  • 技术优化： 使用多网格求解器 (QUDA) 和所有模式平均 (AMA) 技术提高信噪比；使用库仑规范动量平滑 (momentum smearing) 改善 Boost 态的信号。

• 矩阵元提取：

  • 计算包含 Wilson 线的非局域夸克双线性算符的三点函数。
  • 提取同位旋矢量（$u-d$）的螺旋度通道（$\tilde{h}_1$）和 twist-3 横向通道（$\tilde{h}_T$）的裸矩阵元。
  • 通过两态拟合（two-state fit）消除激发态污染。

• 重整化与匹配：

  • 混合重整化方案 (Hybrid Scheme)： 在短距离使用比值方案，长距离减去线性发散。
  • 重整化子重求和 (Renormalon Resummation)： 引入线性重整化子重求和（LRR）来正则化线性发散中的红外模糊性，确保微扰匹配核与重整化方案一致。
  • 微扰匹配： 将重整化后的准 PDF (quasi-PDF) 通过微扰匹配核转换到 $\overline{\text{MS}}$ 方案的光锥 PDF。匹配计算达到了 NNLO 精度，并包含了重整化群 (RG) 和阈值 (Threshold) 重求和，以处理 $x \to 0$ 和 $x \to 1$ 区域的大对数项。

• Twist-3 矩 $\tilde{d}_2$ 的提取：

  • 构建 RG 不变的比值 $M_T$，利用算符乘积展开 (OPE)。
  • 基于 Ref. [94] 发展的 NLO 因子化框架，构造减除项以分离出真正的 twist-3 贡献。
  • 通过分析短距离矩阵元对 Ioffe 时间 $\lambda$ 的依赖关系（$\lambda^2$ 项），直接提取 $\tilde{d}_2$。

• 端点区域约束：

  • 由于 LaMET 在 $x \to 0, 1$ 区域失效，结合 SDF 方法，利用坐标空间矩阵元的短距离行为，通过参数化模型（$x^a(1-x)^b$）拟合端点区域，从而获得全 $x$ 范围的 PDF。

3. 关键贡献 (Key Contributions)

1. 物理点与高精度计算： 首次在物理夸克质量、精细晶格（$a=0.076$ fm）下，结合 NNLO 微扰匹配和重求和技术，计算了质子的同位旋矢量螺旋度 PDF。
2. 首次提取 $\overline{\text{MS}}$ 方案下的 $\tilde{d}_2$： 首次从格点 QCD 的非局域关联函数中，在 $\overline{\text{MS}}$ 方案下直接提取了 twist-3 矩 $\tilde{d}_2$。这避免了传统局域算符方法中复杂的算符混合和方案转换问题。
3. 理论框架的完善： 系统性地应用了线性重整化子重求和（LRR）和 NNLO 匹配，显著降低了微扰展开的不确定性，并明确了 LaMET 方法在中间 $x$ 区域（$x \in [0.25, 0.75]$）的适用范围。
4. SDF 与 LaMET 的结合： 展示了如何利用 SDF 约束端点行为，弥补 LaMET 在大动量展开失效区域的不足，提供了完整的 $x$ 依赖 PDF。

4. 主要结果 (Results)

• 螺旋度 PDF 矩：

  • 提取了 $\mu = 2$ GeV 下的矩：$\langle x \rangle^{u-d} = 0.291(7)$，$\langle x^2 \rangle^{u-d} = 0.101(10)$。
  • 这些矩略高于 JAM22 全局拟合的结果。

• Twist-3 矩 $\tilde{d}_2$：

  • 提取结果：$\tilde{d}_2^{u-d}(2 \text{ GeV}) = 0.0024(46)$。
  • 物理意义： 结果在统计上与零一致，表明在同位旋矢量通道中，最低阶的手征偶 twist-3 夸克 - 胶子关联矩受到强烈抑制。这与唯象学约束和极化深度非弹性散射（DIS）的提取结果一致。

• 螺旋度 PDF $x$ 分布：

  • 在 $x \in [0.25, 0.75]$ 区域，PDF 分布呈现比全球拟合（如 JAM22, BDSSV24）更高的值（特别是在 $x \in [0.4, 0.7]$）。
  • 端点参数化拟合得到 $g_1^+$ 的 $b \approx 1.68(87)$，$g_1^-$ 的 $b \approx 0.15(22)$（后者误差较大）。
  • 与库仑规范（无 Wilson 线）的格点计算结果相比，在 $x < 0.6$ 区域一致性良好，但在 $x > 0.6$ 区域存在差异，这可能源于当前动量较低（$P_z \approx 1.5$ GeV）导致的幂次修正。

5. 意义与展望 (Significance)

• 第一性原理验证： 该工作为质子自旋结构提供了独立于实验的第一性原理约束，特别是在中间 $x$ 区域，有助于解决全球拟合中的不确定性。
• Twist-3 物理的新窗口： 成功提取 $\tilde{d}_2$ 证明了利用非局域算符和短距离因子化研究高阶 twist 物理的可行性，为未来研究更复杂的夸克 - 胶子关联开辟了道路。
• 方法论进步： 结合 NNLO 匹配、重求和及混合重整化方案，展示了格点 QCD 在计算 PDF 方面已达到极高的理论精度。
• 未来方向： 为了进一步减小系统误差（特别是幂次修正和端点区域的不确定性），未来需要在更高动量（$P_z$）下进行计算，并增加统计量和控制激发态。此外，需要更精细地量化端点建模带来的模型误差。

总结： 这是一项具有里程碑意义的格点 QCD 工作，它不仅提供了高精度的质子螺旋度 PDF 数据，还首次实现了对关键 twist-3 矩 $\tilde{d}_2$ 的 $\overline{\text{MS}}$ 方案提取，极大地推动了强子结构从第一性原理角度的理解。