[论文解读] Flow-based density of states for complex actions

Jan M. Pawlowski, Julian M. Urban|arXiv (Cornell University)|Mar 2, 2022

Theoretical and Computational Physics参考文献 56被引用 6

一句话总结

该论文提出一种基于流的新型方法，可直接计算格点场论中复作用量下的态密度（DoS），避免了传统方法中对对数导数进行数值积分时易出错的问题。通过利用归一化流实现直接密度估计，该方法成功定位了李-杨零点，并在零维、一维和二维标量场论中准确复现了DoS，展现出更高的鲁棒性与显式归一化控制能力。

ABSTRACT

Emerging sampling algorithms based on normalizing flows have the potential to solve ergodicity problems in lattice calculations. Furthermore, it has been noted that flows can be used to compute thermodynamic quantities which are difficult to access with traditional methods. This suggests that they are also applicable to the density-of-states approach to complex action problems. In particular, flow-based sampling may be used to compute the density directly, in contradistinction to the conventional strategy of reconstructing it via measuring and integrating the derivative of its logarithm. By circumventing this procedure, the accumulation of errors from the numerical integration is avoided completely and the overall normalization factor can be determined explicitly. In this proof-of-principle study, we demonstrate our method in the context of two-component scalar field theory where the $O(2)$ symmetry is explicitly broken by an imaginary external field. First, we concentrate on the zero-dimensional case which can be solved exactly. We show that with our method, the Lee-Yang zeroes of the associated partition function can be successfully located. Subsequently, we confirm that the flow-based approach correctly reproduces the density computed with conventional methods in one- and two-dimensional models.

研究动机与目标

• 为解决格点场论中复作用量的符号问题，通过归一化流直接计算态密度（DoS）。
• 克服传统方法中通过数值积分对数导数重构DoS时固有的误差累积问题。
• 证明基于流的采样可实现显式归一化，并准确计算配分函数和磁化等热力学量。
• 在精确可解的零维模型上验证该方法，并将其扩展至具显式O(2)对称性破缺的一维和二维格点场论。
• 探索基于流的方法在解决遍历性问题（如拓扑冻结和临界慢化）方面的潜力。

提出的方法

• 该方法利用归一化流建模给定虚部作用量生成器 $ c = X(\phi) $ 条件下的场构型条件分布，实现直接密度估计。
• 流结构采用16个仿射耦合层，耦合参数 $ s $ 和 $ t $ 依赖于上下文 $ c $，并通过全连接或卷积神经网络进行条件化。
• 应用变量变换公式，直接从流的可逆变换计算概率密度 $ \rho(c) $，避免对 $ \partial_c \log \rho(c) $ 进行数值积分。
• 训练过程中使用Adam优化器，学习率为 $ 10^{-3} $，批量大小为10k（0D）或1k（1D/2D），最多进行50k次梯度更新。
• 通过在每个 $ c $-值处采样10k至1000万组构型（$ \Delta c = 0.01 $）来评估DoS，并与基于HMC的传统DoS重构结果进行对比。
• 该方法显式计算配分函数 $ Z = \int dc \, \rho(c) e^{-i h c} $，从而可直接获取磁化等可观测量及李-杨零点。

实验结果

研究问题

• RQ1 能否利用归一化流在复作用量格点场论中直接计算态密度，从而避免对对数导数进行数值积分？
• RQ2 基于流的DoS方法在精确可解的零维模型中是否能准确复现已知结果，包括李-杨零点？
• RQ3 与传统MCMC方法相比，该基于流的方法在一维和二维复作用量标量场论中的表现如何？
• RQ4 该方法能否通过显式采样配置空间中欠采样的区域，解决如拓扑冻结等遍历性问题？
• RQ5 显式归一化与无误差密度估计对磁化、配分函数等热力学可观测量的准确性有何影响？

主要发现

• 基于流的方法在零维模型中成功定位了配分函数的李-杨零点，证实其能准确捕捉非解析行为。
• 在一维和二维模型中，基于流的DoS与传统MCMC重构的DoS高度一致，验证了其可靠性。
• 该方法通过直接计算 $ \rho(c) $ 避免了数值积分带来的误差累积，实现了密度的显式归一化。
• 该方法准确计算了磁化随外加虚部磁场的变化关系，揭示了与相变一致的不连续性。
• 基于流的采样在所有测试的 $ c $-值范围内均获得稳定结果，包括传统重加权因符号问题而失效的区域。
• 该方法在高方差区域表现出强鲁棒性，为复作用量问题提供了一种可行的MCMC替代方案。