[论文解读] Controlling Residual Chiral Symmetry Breaking in Domain Wall Fermion Simulations

Dwight Renfrew, Thomas Blum+3 more|arXiv (Cornell University)|Feb 16, 2009

Quantum Chromodynamics and Particle Interactions参考文献 6被引用 15

一句话总结

该论文提出了一种基于行列式的权重因子，利用威尔逊-狄拉克费米子行列式中的手征性扭曲质量，以抑制域墙费米子（DWF）格点QCD模拟中的残留手征对称性自发破缺。通过调节质量参数，该方法显著降低了残留质量 $m_{\rm res}$，同时保持了足够的拓扑隧穿能力，从而在强耦合区域实现了更大的物理体积和对手征对称性破坏的更好控制。

ABSTRACT

At stronger gauge-field couplings, the domain wall fermion (DWF) residual mass, a measure of chiral symmetry breaking, grows rapidly. This measure is largely due to near zero fermion eigenmodes of logarithm of the 4D transfer matrix along the fifth dimension, and these eigenmodes increase rapidly at strong coupling. To suppress these eigenmodes, we have added to the DWF path integral a multiplicative weighting factor consisting of a ratio of determinants of Wilson-Dirac fermions having a chirally twisted mass with a large negative real component and a small imaginary chiral component. Numerical results show that this weighting factor with an appropriate choice of twisted masses significantly suppresses the residual mass while allowing adequate topological tunneling.

研究动机与目标

• 降低域墙费米子模拟中的残留质量 $m_{\rm res}$，其在强规范耦合下迅速增大，原因在于四维转移矩阵 $T$ 的近零特征值。
• 在蒙特卡罗模拟中保持足够的拓扑隧穿能力，这需要访问具有近零狄拉克特征值的规范配置。
• 开发一种计算上可行的权重因子，通过抑制厄米特威尔逊-狄拉克算符的近零特征值来代理抑制 $H_4 = -\log(T)$ 的近零特征值。
• 在强耦合下实现更大的物理格点体积（例如 $L \approx 4.5$ fm），并控制 $m_{\rm res} \approx 0.002$，对应 $1/a \approx 1.4$ GeV。

提出的方法

• 在DWF路径积分中引入一个乘法权重因子，由具有手征性扭曲质量的威尔逊-狄拉克费米子行列式之比构成：一个具有大负实质量，另一个具有小虚部手征分量。
• 使用厄米特威尔逊-狄拉克算符 $\gamma_5 D_{\mathcal{W}}(-M)$ 的平方行列式作为代理，以抑制驱动 $m_{\rm res}$ 的 $H_4 = -\log(T)$ 的近零特征值。
• 采用分子动力学（MD）演化结合该权重因子，生成排斥力，使系统远离包含近零特征值的规范场构型。
• 调节扭曲质量参数（$\epsilon_f$, $\epsilon_b$）以平衡 $m_{\rm res}$ 的抑制与拓扑荷隧穿的保持。
• 通过数值评估不同参数集下的 $m_{\rm res}$ 和拓扑荷 $\langle |\nu_{\rm top}| \rangle$ 来评估方法的有效性。
• 在零温与非零温条件下比较结果，以确认 $m_{\rm res}$ 对热效应的独立性，并验证在不同耦合强度下的抑制效果。

实验结果

研究问题

• RQ1 具有手征性扭曲质量的行列式之比是否能有效抑制域墙费米子模拟中的残留质量 $m_{\rm res}$？
• RQ2 所提出的权重因子是否能保持足够的拓扑隧穿能力，以确保规范场拓扑的正确采样？
• RQ3 在平衡 $m_{\rm res}$ 抑制与拓扑荷保持方面，扭曲质量参数的最优范围是什么？
• RQ4 厄米特威尔逊-狄拉克算符的近零特征值抑制与 $m_{\rm res}$ 抑制及拓扑荷之间的相关性如何？
• RQ5 该方法是否能实现在 $1/a \approx 1.4$ GeV 下 $32^3 \times 64$ 格点上 $m_{\rm res} \approx 0.002$ 的模拟，对应 $L \approx 4.5$ fm？

主要发现

• 该方法在参数集 $\epsilon_f/\epsilon_b = 0.005/0.50$ 下对 $m_{\rm res}$ 的抑制效果最佳，与未使用任何权重因子的模拟相比，显著降低了其值。
• 在 $\epsilon_f/\epsilon_b = 0.005/0.50$ 时，拓扑荷 $\langle |\nu_{\rm top}| \rangle$ 几乎被完全消除，表明对关键隧穿过程的极接近零的特征值实现了强烈抑制。
• 在强耦合区域，$m_{\rm res}$ 的抑制效果最显著，相对抑制程度在QCD相变温度以下增强，在相变点附近减弱。
• 数值结果确认 $m_{\rm res}$ 几乎与温度无关，表现为在 $\beta = 1.75$ 下，$16^4 \times 32$ 与 $16^3 \times 8 \times 32$ 格点上的值保持一致。
• 特征值范围 $|\lambda| \lesssim 0.1-0.5$ 被识别为 $m_{\rm res}$ 抑制的关键区域，而 $|\lambda| \lesssim 0.01-0.05$ 则主导拓扑荷的抑制。
• 结果表明，在 $1/a \approx 1.4$ GeV 下，$32^3$ 格点上实现 $m_{\rm res} \approx 0.002$ 且 $L \approx 4.5$ fm 是可行的，从而可实现比以往更大的物理体积。