[论文解读] Spanning of Topological sectors, charge and susceptibility with naive Wilson fermions

Abhishek Chowdhury, Asit K. De+5 more|arXiv (Cornell University)|Nov 8, 2011

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

一句话总结

本研究利用HMC与DDHMC算法，在β = 5.6和5.8下对两种简并的外尔威尔逊费米子进行格点QCD模拟，研究了拓扑电荷与拓扑敏感度。结果明确显示，随着夸克质量减小，拓扑敏感度被抑制，证实了手征微扰理论的预测，并解决了以往使用未改进费米子模拟中的长期模糊性。

ABSTRACT

We study the topological charge and the topological susceptibility in lattice QCD with two degenerate flavors of naive Wilson fermions at two values of lattice spacings and different volumes, for a range of quark masses. Configurations are generated with DDHMC/HMC algorithms and smoothened with HYP smearing. We present integrated autocorrelation time for both topological charge and topological susceptibility at the two lattice spacing values studied. The spanning of different topological sectors as a function of the hopping parameter kappa is presented. The expected chiral behaviour of the topological susceptibility (including finite volume dependence) is observed.

研究动机与目标

• 解决长期以来在格点QCD模拟中关于使用未改进威尔逊费米子时，随着夸克质量减小，拓扑敏感度是否被抑制的模糊性。
• 研究不同跃迁参数（κ）和格点体积下，拓扑规范层的跨越行为。
• 评估HMC与DDHMC算法在使用非改进威尔逊费米子时，访问拓扑规范层的效率与准确性。
• 通过在多个格点间距与体积下进行系统模拟，验证拓扑敏感度的预期手征行为，包括有限体积效应。
• 为未来比较提供精确的自关联时间估计，作为拓扑电荷与敏感度的基准结果。

提出的方法

• 使用HMC与DDHMC算法，基于未改进的威尔逊费米子和规范作用量，生成了格点配置，涉及两种简并的夸克味。
• 通过索莫默参数确定格点间距：β = 5.6时为0.077 fm，β = 5.8时为0.061 fm。
• 拓扑电荷通过格点近似处理的陈-西蒙斯电流计算，该方法源自DeGrand、Hasenfratz与Kovacs对SU(3)的处理，并在MILC代码中实现。
• 采用20步优化系数（α = 0.75，α₂ = 0.6，α₃ = 0.3）的HYP投影，以抑制紫外格点效应，并改善拓扑电荷的整数实现。
• 通过蒙特卡罗轨迹的分箱统计分析，计算了拓扑电荷与敏感度的积分自关联时间。
• 通过轴矢量Ward恒等式非微扰地确定夸克质量，拓扑敏感度定义为χ = (⟨Q²⟩ − ⟨Q⟩²)/V。

实验结果

研究问题

• RQ1 在使用非改进威尔逊费米子的格点QCD中，拓扑敏感度是否如手征微扰理论所预测的那样，随夸克质量减小而被抑制？
• RQ2 格点间距、体积与夸克质量如何影响使用未改进费米子模拟中拓扑规范层的跨越行为？
• RQ3 HYP投影在多大程度上改善了拓扑电荷的整数实现，相较于APE投影？
• RQ4 在不同κ与β值下，HMC与DDHMC模拟中拓扑电荷与敏感度的积分自关联时间是多少？
• RQ5 观测到的拓扑敏感度行为能否与手征Ward恒等式及一阶手征微扰理论预测相一致？

主要发现

• 在β = 5.6与β = 5.8下，拓扑敏感度均表现出随夸克质量减小而明显抑制，证实了手征微扰理论的预测。
• 在β = 5.6与体积32³×64下，当夸克质量从100 MeV降至15 MeV（MS̄在2 GeV下）时，拓扑敏感度从~1.2×10⁻⁴降至~0.3×10⁻⁴ GeV⁴。
• 在β = 5.6下，拓扑电荷的积分自关联时间τ_int^Q在30至590之间变化，最大值出现在临界κ附近（例如，κ = 0.15825，β = 5.6，24³×48时为590(48)）。
• HYP投影显著减少了拓扑电荷的涨落，并使数值更接近整数，β = 5.6、κ = 0.15825下，20步HYP投影后98%的配置显示Q = 0。
• 在β = 5.8下，κ = 0.15475时观察到拓扑冻结的微弱迹象，τ_int^Q = 211(30)，τ_int^χ = 109(14)，表明拓扑规范层内混合缓慢。
• 拓扑电荷分布近似为高斯分布，⟨Q⟩ ≈ 0，且⟨Q²⟩与χV一致，支持在大多数配置中不存在大尺度的拓扑冻结。