[论文解读] A (P)HMC algorithm for N_f = 2 + 1 + 1 flavours of twisted mass fermions

T. Chiarappa, R. Frezzotti+1 more|ArXiv.org|Sep 27, 2005

Quantum Chromodynamics and Particle Interactions被引用 4

一句话总结

本文提出了一种用于$N_f = 2+1+1$扭曲质量费米子的(P)HMC算法，采用偶奇预条件化、质量位移预条件化及多时间尺度方法。该方法引入了对重夸克区中狄拉克矩阵平方根的有理函数逼近，通过伪费米子场和重加权技术，实现了非简并奇异夸克与粲夸克的高效且稳定的模拟。

ABSTRACT

We present a detailed design of a (P)HMC simulation algorithm for N_f = 2 + 1 + 1 maximally twisted Wilson quark flavours. The algorithm retains even/odd and mass-shift preconditionings combined with multiple Molecular Dynamics time scales for both the light mass degenerate, u and d, quarks and the heavy mass non-degenerate, s and c, quarks. Various non-standard aspects of the algorithm are discussed, among which those connected to the use of a polynomial approximation for the inverse (square root of the squared) Dirac matrix in the s and c quark sector.

研究动机与目标

• 开发一种稳定且高效的(P)HMC算法，用于$N_f = 2+1+1$扭曲质量费米子，涵盖轻夸克与重夸克区。
• 解决通过双味伪费米子场表示法模拟非简并重夸克（奇异与粲夸克）的挑战。
• 通过狄拉克矩阵平方根的多项式逼近，确保费米子行列式计算的数值稳定性和准确性。
• 结合多时间尺度与预条件化技术，提升分子动力学更新的性能。

提出的方法

• 该算法使用双味伪费米子场$\Phi$表示$h$-夸克狄拉克算符$\tilde{Q}_{2\times 2}$的行列式，从而实现PHMC方法。
• 采用多项式$P_n(s) \approx 1/\sqrt{s}$逼近矩阵$\hat{S} = \tilde{Q}_h \tilde{Q}_h^\dagger$的平方根倒数，确保正定性与计算稳定性。
• 多项式逼近构造为$P_{n,s_L}(\hat{S}) = \mathcal{N} \prod_{k=1}^n (\hat{S} - z_k)$，其根关于复共轭对称，以维持正定性。
• 使用第二个高精度多项式$\tilde{P}_{\tilde{n},s_L}(\hat{S})$生成伪费米子场并用于重加权，相对误差$\tilde{\delta}_{IR} \leq 10^{-6}$，当$\tilde{n} = 620$且$s_L = 10^{-4}$时成立。
• 该方法结合偶奇预条件化与质量位移预条件化，并采用轻夸克与重夸克区的多时间步分子动力学更新。
• 算法采用64位浮点数运算，通过稳定递归算法（如Clenshaw递推）控制多项式计算中的舍入误差。

实验结果

研究问题

• RQ1 如何设计一种(P)HMC算法，以高效模拟$N_f = 2+1+1$扭曲质量费米子中的非简并重夸克？
• RQ2 在重夸克区中，狄拉克矩阵平方根的倒数最优多项式逼近策略为何？
• RQ3 在$\hat{S}$谱中存在低本征值的情况下，如何维持计算的数值稳定性和准确性？
• RQ4 多时间尺度与预条件化对HMC算法性能与接受率有何影响？
• RQ5 如何利用多项式逼近高精度生成所需的伪费米子场？

主要发现

• 当多项式次数为$n = 162$时，相对拟合误差满足$|R_{n,s_L}(s)| \leq 0.03$，$s \in [10^{-4}, 1]$，且以$\delta_{IR} = |R_{n,s_L}(10^{-4})|$作为保守误差界。
• 第二个高精度多项式（次数为$\tilde{n} = 620$）在相同区间内实现相对误差$\tilde{\delta}_{IR} \leq 10^{-6}$，支持精确的伪费米子生成与重加权。
• 采用双味伪费米子场$\Phi$可将非简并$h$-夸克狄拉克算符的行列式表示为$\det[\tilde{Q}_{2\times 2}] \propto \int \mathcal{D}\Phi \, e^{-\Phi^\dagger P(\hat{S}) \Phi}$，从而实现PHMC。
• 该算法在64位浮点数运算下具备数值稳定性，通过Clenshaw递推实现多项式稳定递归计算。
• 该方法支持多时间尺度与预条件化（偶奇与质量位移），显著提升了分子动力学更新的效率。
• 该方法实现了对$N_f = 2+1+1$扭曲质量费米子的精确模拟，误差可控，其有效性通过关键低本征值区域中极低的相对拟合误差得到验证。