[论文解读] Stochastic calculation of the QCD Dirac operator spectrum with Mobius domain-wall fermion

Guido Cossu, Hidenori Fukaya+3 more|arXiv (Cornell University)|Jan 5, 2016

Quantum Chromodynamics and Particle Interactions被引用 3

一句话总结

本文提出了一种基于切比雪夫滤波与特征值计数的随机方法，用于在莫比乌斯域墙费米子框架下计算量子色动力学（QCD）狄拉克算符的完整谱。该方法实现了整个特征值范围内精确的单次测量谱分析，通过一阶圈级有效场论拟合，得到的轴向矢量凝聚能为 $[\Sigma(2~\text{GeV})]^{1/3} = 260.0(1.7)~\text{MeV}$。

ABSTRACT

We calculate the spectral function of the QCD Dirac operator using the four-dimensional effective operator constructed from the Mobius domain-wall implementation. We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number. The spectrum in the entire eigenvalue range is obtained with a single set of measurements. Results on 2+1-flavor ensembles with Mobius domain-wall sea quarks at lattice spacing ~ 0.08 fm are shown.

研究动机与目标

• 开发一种在格点QCD中计算QCD狄拉克算符完整谱的计算高效方法。
• 通过实现随机、灵活的谱估计，克服传统特征值求解器在大体积模拟中的局限性。
• 利用Banks-Casher关系与有效场论，非微扰地提取轴向矢量凝聚能。
• 在 $a \sim 0.08$ fm 的 2+1 条味莫比乌斯域墙费米子组态上验证该方法，其海夸克质量取值不同。

提出的方法

• 利用切比雪夫多项式滤波，随机估计狄拉克算符 $D^\dagger D$ 在给定区间 $[a,b]$ 内的特征值数量。
• 采用包含 $N_v$ 个高斯随机向量的随机估计器，计算 $\bar{n}[a,b] = \frac{1}{N_v} \sum_{k=1}^{N_v} \xi_k^\dagger h(A) \xi_k$，其中 $h(A)$ 充当谱滤波器。
• 应用切比雪夫展开 $h(x) = \sum_{j=0}^{p} g_j^p \gamma_j T_j(x)$，并引入杰克逊阻尼以抑制振荡并改善收敛性。
• 利用递推关系 $T_j(x) = 2xT_{j-1}(x) - T_{j-2}(x)$ 实现高效多项式求值，降低计算成本。
• 使用莫比乌斯域墙费米子的4维有效狄拉克算符，其 $D^{(4)\dagger}D^{(4)}$ 的特征值被限制在 $[0,1]$ 范围内，从而在适中多项式阶数下实现高精度滤波。
• 通过一阶圈级 $N_f=2$ 有效场论实现轴向矢量凝聚能的轴向外推，从零附近谱密度中提取 $\Sigma$。

实验结果

研究问题

• RQ1 切比雪夫滤波与随机估计相结合，能否为格点QCD中完整狄拉克谱的计算提供一种高效且精确的方法？
• RQ2 在 $p$-区域中，谱密度 $\rho(\lambda)$ 如何依赖于海夸克质量 $m_{ud}$？
• RQ3 一阶圈级有效场论在多大程度上能描述零附近 $\rho(\lambda)$ 的格点数据？
• RQ4 利用非微扰格点数据从谱密度中提取的轴向矢量凝聚能 $\Sigma$ 的值是多少？

主要发现

• 该方法通过单次测量集合成功计算出完整的狄拉克谱，采用随机切比雪夫滤波。
• 谱函数 $\rho(\lambda)$ 明显依赖于 $m_{ud}$，对于较重的海夸克，零附近出现峰值，这是由于近零模被抑制程度降低所致。
• 轴向矢量凝聚能提取为 $[\Sigma(2~\text{GeV})]^{1/3} = 260.0(1.7)~\text{MeV}$，与Banks-Casher关系一致。
• 对于 $\delta=0.01$，一阶圈级有效场论在 $\lambda=0$ 附近对数据的描述良好，$\chi^2/\text{dof} = 1.13$。
• 切比雪夫近似带来的系统误差低于统计误差，$\delta=0.01$ 和 $\delta=0.005$ 时的误差分别为 0.8% 和 1.5%。
• 2+1 条味有效场论拟合无法再现数据，可能是因为奇异夸克质量过重，导致一阶近似失效。