[论文解读] A stochastic approach to the reconstruction of spectral functions in lattice QCD

Hai-Tao Shu, Heng-Tong Ding+3 more|arXiv (Cornell University)|Oct 10, 2015

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

一句话总结

本文提出了一种随机优化方法（SOM），用于在无需依赖贝叶斯推断或默认模型的情况下，从格点QCD中的欧几里得相关函数重建谱函数。通过将谱函数参数化为一组随机更新的矩形框，并利用模拟退火法最小化似然函数，SOM成功重建了包含共振态、输运峰和连续谱成分的模拟谱函数，其结果与最大熵方法（MEM）在真实底粲强子相关函数上的表现相当。

ABSTRACT

We present a Stochastic Optimization Method (SOM) for the reconstruction of the spectral functions (SPFs) from Euclidean correlation functions. In this approach the SPF is parameterized as a sum of randomly distributed boxes. By varying the width, location and height of the boxes stochastically an optimal SPF can be obtained. Using this approach we reproduce mock SPFs fairly well, which contain sharp resonance peaks, transport peaks and continuum spectra. We also analyzed the charmonium correlators obtained from $N_τ$=96, 48, 32 lattices using SOM and found similar conclusion on the dissociation temperatures of charmonium ground states as that obtained using the Maximum Entropy Method.

研究动机与目标

• 开发一种无需贝叶斯方法的谱函数重建方法，用于从格点QCD相关函数中进行重建。
• 解决从欧几里得空间到闵可夫斯基空间的解析继续这一不适定逆问题。
• 在包含复杂谱结构（如尖锐共振态、输运峰和连续谱贡献）的模拟数据上测试该方法的鲁棒性。
• 将SOM的结果与真实格点数据（$N_{\tau} = 96, 48, 32$）下最大熵方法（MEM）的结果进行比较。
• 提供一种可靠的MEM替代方案，避免对默认模型或先验信息的依赖。

提出的方法

• 谱函数被参数化为$K$个随机分布的矩形框之和，其宽度、高度和中心位置可变。
• 定义似然函数$Q$，用于比较重建的相关函数$\widetilde{G}(\tau)$与平均格点相关函数$\overline{G}(\tau)$，其中使用了来自配置平均的协方差矩阵$C$。
• 通过随机应用基本更新操作（移动、缩放、合并或分裂）来探索解空间。
• 采用模拟退火法最小化$Q$，冷却计划通过$\delta Q/T$控制接受概率。
• 通过平均$L=100$个独立解来获得最终的谱函数，以减少噪声和波动。
• 通过$\widetilde{G}(\tau) = \int_0^\infty \frac{d\omega}{2\pi} \widetilde{\rho}(\omega) K(\omega,\tau,T)$计算重建的相关函数，其中核$K$编码了热场理论边界条件。

实验结果

研究问题

• RQ1 是否能够通过一种随机优化方法，在不依赖贝叶斯先验的情况下，重建包含尖锐共振态、输运峰和连续谱成分的谱函数？
• RQ2 SOM在具有复杂谱结构的模拟数据上的表现与最大熵方法（MEM）相比如何？
• RQ3 SOM是否能够可靠地从$N_{\tau} = 96, 48, 32$的真实格点QCD相关函数中提取谱函数？
• RQ4 对于矢量通道中的底粲强子，SOM与MEM在峰位置和高度上的结果有何差异？
• RQ5 连续谱成分的存在是否显著影响SOM中尖锐共振峰的重建？

主要发现

• SOM能够以高保真度成功重建包含尖锐共振态、输运峰和连续谱成分的模拟谱函数。
• 对于截距较大的模拟数据（截距为17），该方法能较好地再现连续谱，但由于噪声和高$\omega$区域的盒分布偏差，对截距的重建存在困难。
• 在真实$N_{\tau} = 96, 48, 32$的底粲强子相关函数上，SOM成功再现了$T = 0.73T_c$处的输运峰，并在$T > T_c$时显示出共振态抑制现象。
• 在$T = 0.73T_c$时，SOM计算得到的底粲强子第一峰位置为3.44 GeV，与空间相关函数测得的3.472 GeV屏蔽质量非常接近。
• SOM得到的峰位置（3.44 GeV）与MEM结果（3.48 GeV）非常接近，但SOM的峰高更大（1.75 vs. 0.98）。
• 对100个解的平均有效抑制了波动，改善了在噪声存在下对连续谱和共振特征的重建。