[论文解读] Exact solutions of the QCD evolution equations using Monte Carlo method

S. Jadach, M. Skrzypek|ArXiv.org|Dec 26, 2003

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

一句话总结

本文提出了一种基于马尔可夫蒙特卡洛方法的精确、高精度数值求解QCD演化方程的方法，实现了在广泛 $Q$ 和 $x$ 范围内约 0.1% 的精度。结果表明，现代CPU算力使得基于蒙特卡洛方法的精确部分子分布演化成为可行，同时提出了对发展非马尔可夫蒙特卡洛算法的挑战，该算法可内建完整的QCD演化以处理初态辐射，而无需依赖预计算的PDF表格。

ABSTRACT

We present the exact and precise (~0.1%) numerical solution of the QCD evolution equations for the parton distributions in a wide range of $Q$ and $x$ using Monte Carlo (MC) method, which relies on the so-called Markovian algorithm. We point out certain advantages of such a method with respect to the existing non-MC methods. We also formulate a challenge of constructing non-Markovian MC algorithm for the evolution equations for the initial state QCD radiation with tagging the type and $x$ of the exiting parton. This seems to be within the reach of the presently available computer CPUs and the sophistication of the MC techniques.

研究动机与目标

• 测试利用现代CPU算力通过蒙特卡洛方法精确求解QCD演化方程的实际可行性。
• 证明马尔可夫蒙特卡洛方法可在广泛 $Q$ 和 $x$ 范围内实现约 0.1% 的精度求解演化方程。
• 挑战开发非马尔可夫蒙特卡洛算法，以整合初态辐射的完整QCD演化，同时标记部分子类型和 $x$-分数。
• 倡导将精确蒙特卡洛演化集成到未来的部分子喷注生成器中，以实现LHC和DIS数据的同时拟合。
• 探索用动态演化分布替代预计算PDF表格在事件生成器中的可行性。

提出的方法

• 将马尔可夫蒙特卡洛方法应用于求解QCD演化方程，将其解释为在 $t = \ln Q$ 上的随机过程，包含 $t$ 方向的随机步长和部分子衰变。
• 该方法使用演化方程的积分形式，重写为一系列被积函数为正的迭代积分，可解释为具有适当归一化转移概率的马尔可夫过程。
• 通过分裂函数和衰变速率定义转移概率 $\omega(t_i, K_i | t_{i-1}, K_{i-1})$，确保归一化和概率一致性。
• 该算法通过在 $t$ 和部分子态空间中的随机游走模拟部分子喷注，最终在目标 $Q$ 尺度处设置停止规则。
• 该方法避免了基于网格的离散化，转而依赖演化路径的随机采样以精确计算 $D_I(t)$。
• 该方法被扩展用于考虑非马尔可夫算法处理初态辐射，其中在高 $Q$ 处标记部分子类型和 $x$，采用暴力或优化采样。

实验结果

研究问题

• RQ1 马尔可夫蒙特卡洛方法是否能在广泛 $Q$ 和 $x$ 范围内实现QCD演化方程的精确且高精度（约0.1%）数值求解？
• RQ2 与传统的非蒙特卡洛网格求解器相比，使用蒙特卡洛方法在实际应用中有哪些优势？
• RQ3 如今是否可行构建一种非马尔可夫蒙特卡洛算法，以包含初态辐射的完整QCD演化，并对部分子类型和 $x$-分数进行标记？
• RQ4 现代CPU算力和蒙特卡洛技术是否足以支持一种非马尔可夫蒙特卡洛算法，实现从低 $Q$ 到高 $Q$ 的PDF演化，而无需依赖预计算表格？
• RQ5 将精确蒙特卡洛演化集成到部分子喷注生成器中，对同时拟合LHC和DIS数据有何影响？

主要发现

• 马尔可夫蒙特卡洛方法在 $Q$ 从 1 GeV 到 1 TeV 的范围内求解QCD演化方程时，实现了约 0.1% 的精度。
• 从 $Q = 1$ GeV 到 $Q = 1$ TeV 的质子喷注中，平均部分子发射次数约为 20 次，表明存在显著的辐射。
• 平均 $q,\bar{q} \to G$ 和 $G \to q,\bar{q}$ 转换次数较低，介于 0.5 到 1 之间，表明这些过程可预先制表。
• 研究表明，由于过去15年CPU算力提升了10,000倍，标准马尔可夫蒙特卡洛方法在初态辐射中效率低的问题已不再是根本障碍。
• 作者认为，具备内建演化的非马尔可夫蒙特卡洛算法如今已可行，可实现更精确、自洽的事件生成器。
• 作者提出一项挑战：开发一种蒙特卡洛算法，使用NLL $\overline{\text{MS}}$ 核函数，从 $Q_{\min}$ 到 $Q_{\max}$ 演化PDF，固定 $Q_{\max}$ 处的部分子类型和 $x$，且无需依赖预计算的PDF表格。