[论文解读] Testing a generalized cooling procedure in the complex Langevin simulation of chiral Random Matrix Theory

Keitaro Nagata, Jun Nishimura+1 more|arXiv (Cornell University)|Nov 27, 2015

Quantum many-body systems参考文献 18被引用 5

一句话总结

本文在手征性随机矩阵理论（chRMT）的复 Langevin 方法（CLM）中提出了一种广义的规范冷却程序，利用基于 Dirac 算符的范数来抑制轻夸克和近零特征值引起的奇异漂移。该方法通过迭代冷却减少原点附近的特征值密度，在标准 CLM 失效的轻夸克区域成功重现了精确结果，展示了在符号问题区域中改进的收敛性和稳定性。

ABSTRACT

The complex Langevin method has been attracting much attention as a solution to the sign problem since the method was shown to work in finite density QCD in the deconfined phase by using the so-called gauge cooling procedure. Whether it works also in the confined phase with light quarks is still an open question, though. In order to shed light on this question, we apply the method to the chiral Random Matrix Theory, which describes the epsilon regime of finite density QCD. Earlier works reported that a naive implementation of the method fails to reproduce the known exact results and that the problem can be solved by choosing a suitable coordinate. In this work we stick to the naive implementation, and show that a generalized gauge cooling procedure can be used to avoid the problem.

研究动机与目标

• 为低温下具有轻夸克的有限密度 QCD 中的费米子奇异漂移问题提供解决方案。
• 克服在轻夸克情况下标准 CLM 在 chRMT 中失效的问题，此时由于 $D+m$ 的特征值接近零，漂移项发散。
• 开发一种不依赖坐标变换的规范冷却程序，同时保持笛卡尔变量。
• 测试基于 Dirac 算符的范数是否能有效控制动力学并确保收敛到精确结果。

提出的方法

• 引入两种新的规范冷却范数：$\mathcal{N}_1 = \mathrm{Tr}[(D+m)(D+m)^\dagger]$ 和 $\mathcal{N}_2 = \mathrm{Tr}[(D+m)^\dagger(D+m)]^{-1}$，对靠近零的特征值分布敏感。
• 在每次 Langevin 步骤后应用冷却程序，通过最小化所选范数的复化 $\mathrm{GL}(N,\mathbb{C})$ 变换来实现。
• 使用基于梯度的优化方法确定变换的大小和方向，每步重复 10 次以确保有效降低。
• 在 $N=30$，$\sqrt{N}\mu=2$，$N_f=2$，$\tilde{m}=mN$ 的 chRMT 中实现该方法，使用 $\Delta t = 5 \times 10^{-5}$，并舍弃 20,000 步热化。
• 在热化后每 200 步测量一次可观测量（手征性凝聚、重子密度），并与精确解析解进行比较。
• 使用 $r=0$ 控制 $\mathcal{N}_1$ 的冷却强度，$r=0.01$ 控制 $\mathcal{N}_2$ 的冷却强度，$\xi=300$，$n_{\rm ev}=2$ 用于 $\mathcal{N}_2$。

实验结果

研究问题

• RQ1 基于 Dirac 算符特征值结构的广义规范冷却程序是否能稳定 chRMT 中轻夸克区域的复 Langevin 方法？
• RQ2 使用 $\mathcal{N}_1$ 或 $\mathcal{N}_2$ 范数的冷却是否能有效抑制由 $D+m$ 特征值接近零引起的危险漂移项？
• RQ3 采用新冷却程序的 CLM 是否能在小 $mN$ 区域重现手征性凝聚和重子数密度的精确结果？
• RQ4 冷却范数的选择如何影响复 Langevin 模拟中 $D+m$ 的特征值分布？
• RQ5 该方法是否可推广至无需坐标重参数化的有限密度 QCD？

主要发现

• 对于 $\tilde{m} = mN \lesssim 10$，朴素的复 Langevin 方法失效，手征性凝聚和重子数密度与精确结果存在显著偏差。
• 采用 $\mathcal{N}_1$ 和 $\mathcal{N}_2$ 范数的新规范冷却程序，CLM 在整个 $\tilde{m}$ 范围内成功重现了精确结果，包括轻夸克区域。
• 在两种冷却程序下，$D+m$ 在原点附近的特征值分布均显著被抑制，其中 $\mathcal{N}_1$ 在实轴方向压缩了分布，$\mathcal{N}_2$ 则将特征值从原点排斥。
• 手征性凝聚（一个全纯可观测量）无论使用何种冷却范数，均收敛至精确值，表明在正确稳定下方法具有一致性和普遍性。
• 即使在 $\tilde{m} = 3$ 时，朴素 CLM 完全失效，该方法仍保持稳定和准确，显著扩展了 CLM 的适用范围。
• 冷却程序有效控制了动力学变量在虚轴方向的长距离漂移，防止了复 Langevin 理论基础的崩溃。