[论文解读] Lattice Calculation of the Connected Hadronic Light-by-Light Contribution to the Muon Anomalous Magnetic Moment

Luchang Jin, Thomas Blum+4 more|arXiv (Cornell University)|Sep 28, 2015

Particle physics theoretical and experimental studies参考文献 7被引用 7

一句话总结

本文提出了一种新颖的格点QCD方法，以高精度计算强子光子-光子自能（cHLbL）对μ子反常磁矩的贡献。通过使用精确的光子传播子和基于矩的方法，相比以往方法，显著降低了统计误差和计算成本，得到在$32^3 \times 64$格点上、具有物理类似π介子和μ子质量下的结果$ F_2 = (0.1054 \pm 0.0054)(\alpha/\pi)^3 $。

ABSTRACT

The anomalous magnetic moment of muon, $g-2$, is a very precisely measured quantity. However, the current measurement disagrees with standard model by about 3 standard deviations. Hadronic vacuum polarization and hadronic light by light are the two types of processes that contribute most to the theoretical uncertainty. I will describe how lattice methods are well-suited to provide a first-principle's result for the hadronic light by light contribution, the various numerical strategies that are presently being used to evaluate it, our current results and the important remaining challenges which must be overcome.

研究动机与目标

• 通过使用第一性原理的格点QCD计算连接的强子光子-光子自能（cHLbL）贡献，解决μ子反常磁矩中较大的理论不确定性。
• 克服以往格点方法在计算cHLbL振幅时存在的噪声、系统误差和高计算成本等局限性。
• 开发并验证一种基于精确光子传播子和矩方法的新评估策略，实现在有限体积中直接计算$ F_2(q^2=0) $。
• 相比早期方法，实现更高的精度和更低的计算成本，使模拟更接近物理动力学条件。
• 为未来在物理π介子质量及更大体积下的计算奠定基础，包括有限体积效应和真空泡图的处理。

提出的方法

• 使用精确解析表达式和傅里叶变换计算光子传播子，避免随机噪声，消除了低阶噪声源。
• 采用基于矩的方法，直接在$ q^2 = 0 $处评估cHLbL振幅，避免外推，降低统计不确定性。
• 将矩方法应用于分离短程和长程贡献：通过直接采样夸克传播子对计算短程部分，通过加权采样$ 1/|r|^4 \exp(-0.01|r|) $分布计算长程部分。
• 使用AMA（代数多重网格加速）技术加速夸克传播子的求解，尤其适用于大体积和细网格模拟。
• 在Breit参考系中计算μ子自能，使初态和末态μ子能量相等，以确保动量守恒并正确投影到单粒子态。
• 利用光子-μ子顶点函数的有限体积矩阵元提取$ F_2(0) $，即反常磁矩贡献，有限体积修正可控。

实验结果

研究问题

• RQ1 是否可以不依赖随机规范场，仅通过格点QCD高精度计算μ子g-2的连接强子光子-光子自能贡献？
• RQ2 与以往格点策略相比，使用精确光子传播子和矩方法如何降低统计误差和计算成本？
• RQ3 在物理π介子和μ子质量下，cHLbL贡献的大小是多少？有限体积效应如何随体积变化？
• RQ4 该新方法能否扩展至物理π介子质量和更大格点，同时保持精度和效率？
• RQ5 cHLbL振幅的有限体积误差结构是什么？能否系统性地控制？

主要发现

• 新方法显著降低了统计误差和计算成本：在$ 24^3 \times 64 $格点上，$ m_\pi = 333\,\text{MeV} $时，$ N_{\text{prop}} = 26,568 $，相比之前方法的$ 81,000 $大幅减少。
• 在$ 32^3 \times 64 $格点上，$ m_\pi = 171\,\text{MeV} $，$ m_\mu = 134\,\text{MeV} $时，结果为$ F_2 = (0.1054 \pm 0.0054)(\alpha/\pi)^3 $，总计算成本为13.2 BG/Q机架天。
• 有限体积误差被证实按$ \mathcal{O}(1/L^2) $缩放，支持该方法在无限体积外推中的可靠性。
• 该方法可直接计算$ F_2(q^2=0) $，无需外推，提高了精度并减少了系统性不确定度。
• 结果与模型估计值$ (0.08 \pm 0.02)(\alpha/\pi)^3 $一致，尽管是在非物理质量和有限体积下计算的。
• 该方法现正应用于物理π介子质量下的$ 48^3 \times 5.5\,\text{fm} $格点模拟，计划引入有限体积修正和真空泡图。