[论文解读] Towards thermodynamics with $N_f=2+1+1$ twisted mass quarks

Florian Bürger, Grit Hotzel+3 more|arXiv (Cornell University)|Nov 7, 2013

High-Energy Particle Collisions Research参考文献 14被引用 4

一句话总结

本研究首次报告了使用 $N_f=2+1+1$ 动态扭曲质量夸克的格点QCD热力学结果，采用固定尺度方法以最小化 $T=0$ 模拟。与 $N_f=2$ 模拟相比，发现手征对称性交叉转变温度约降低 20 MeV，且迹异常的胶子部分表现出更高、更宽的峰，表明第二代夸克对状态方程有显著影响。

ABSTRACT

We present preliminary results achieved within a recently started project dealing with QCD thermodynamics in the presence of a fully dynamical second quark family. We are employing the Wilson twisted mass discretization. To reduce the amount of zero temperature simulations and the cost of analysis we have chosen the fixed-scale approach. We show a variety of basic thermodynamic observables for temperatures ranging from 158 to 633 MeV. Simulations were performed for three lattice spacings below 0.1 fm each and at a single value of the pion mass which allows a comparison with previously obtained $N_f=2$ results. We determine the chiral crossover temperature from the bare chiral susceptibility and show results for the gauge part of the trace anomaly.

研究动机与目标

• 使用威尔逊扭曲质量费米子研究包含一个完全动态的第二代夸克（粲夸克与奇异夸克）的QCD热力学。
• 通过采用固定尺度方法降低计算成本，最小化 $T=0$ 模拟，并实现与现有 $T=0$ 组合的比较。
• 利用有限温格点模拟确定 $N_f=2+1+1$ 夸克存在下的手征对称性交叉转变温度和状态方程。
• 评估截断效应，并在多个格点间距和温度范围内验证结果。

提出的方法

• 轻夸克与重夸克采用威尔逊扭曲质量费米子作用量，胶子部分使用伊瓦萨基规范作用量。
• 在固定逆温度 $\beta$ 下进行有限温模拟，改变 $N_\tau$，对应温度范围为 158 至 633 MeV。
• 固定尺度方法使用来自ETMC的 $T=0$ 组合，格点间距分别为 $a \approx 0.086$、0.078 和 0.061 fm，且经过调节的π介子质量约为 400 MeV。
• 通过 $T=0$ 时在距离 $r_0$ 处的静态 $\bar{Q}Q$ 格势对 Polyakov 圈和手征凝聚进行重正化，并减去 $T=0$ 时的发散项。
• 通过 $I_g = a \frac{d\beta}{da} \langle S_g \rangle_{\text{sub}}$ 计算迹异常的胶子部分，其中 $\beta$-函数通过拟合 $\beta$ 与 $m_\rho$ 和 $m_\pi$ 确定。
• 通过改变拟合函数估计 $\beta$-函数中的系统误差，并将误差传播至 $I_g$。

实验结果

研究问题

• RQ1 在QCD热力学中，包含一个动力学第二代夸克（$N_f=2+1+1$）如何影响手征对称性交叉转变温度？
• RQ2 与 $N_f=2$ 模拟相比，第二代夸克对状态方程（尤其是迹异常）有何影响？
• RQ3 截断效应在多大程度上影响热力学可观测量，特别是在交叉区域附近？
• RQ4 使用 $T=0$ 组合的固定尺度方法在多大程度上能实现高精度的有限温研究，同时最小化 $T=0$ 模拟？

主要发现

• 基于非连通手征敏感度峰值，$N_f=2+1+1$ 模拟中手征对称性交叉转变温度估计比 $N_f=2$ 情况低约 20 MeV。
• 与 $N_f=2$ 相比，$N_f=2+1+1$ 中迹异常的胶子部分 $I_g$ 显示出显著更高且更宽的峰，表明偏离理想气体行为更明显。
• 在 $N_\tau \geq 8$ 时，交叉区域附近的截断效应较小，但在高温区域（$T \gtrsim 350$ MeV）变得显著，尤其在更细格点上。
• 重正化的 Polyakov 圈在交叉区域单调增加，最粗格点的数据在相变区域与其它数据吻合良好。
• $\beta$-函数通过拟合 $m_\rho$ 和 $m_\pi$ 成功获得，$\chi^2/\text{dof} \approx 1.0$，系统误差通过改变拟合函数估计。
• 该方法成功计算了迹异常且误差可控，结果在交叉区域内各格点间距间保持一致。