[论文解读] The {\it victory} project v1.0: an efficient parquet equations solver

Gang Li, Anna Kauch+2 more|arXiv (Cornell University)|Aug 24, 2017

Physics of Superconductivity and Magnetism参考文献 11被引用 4

一句话总结

本论文提出 Victory v1.0，一种用于 Hubbard 模型的抛物线方程的高性能 Fortran 实现，可自洽求解单粒子与双粒子关联函数。其采用两级核近似与 Matsubara 频率空间的高频调节，确保收敛性与准确性，并通过特征值分析揭示相变信号。

ABSTRACT

{\it Victory}, i.e. \underline{vi}enna \underline{c}omputational \underline{to}ol deposito\underline{ry}, is a collection of numerical tools for solving the parquet equations for the Hubbard model and similar many body problems. The parquet formalism is a self-consistent theory at both the single- and two-particle levels, and can thus describe individual fermions as well as their collective behavior on equal footing. This is essential for the understanding of various emergent phases and their transitions in many-body systems, in particular for cases in which a single-particle description fails. Our implementation of {\it victory} is in modern Fortran and it fully respects the structure of various vertex functions in both momentum and Matsubara frequency space. We found the latter to be crucial for the convergence of the parquet equations, as well as for the correct determination of various physical observables. In this release, we thoroughly explain the program structure and the controlled approximations to efficiently solve the parquet equations, i.e. the two-level kernel approximation and the high-frequency regulation.

研究动机与目标

• 开发一种在强关联电子系统背景下数值高效且可扩展的抛物线方程求解器。
• 解决在 Matsubara 频率空间中依赖于三个动量与三个频率的大型双粒子顶点函数带来的计算瓶颈。
• 通过恰当处理频率空间中的开放边界条件，确保收敛性与物理准确性。
• 通过计算不同通道中 Bethe-Salpeter 方程的特征值，实现对相变的探测。
• 提供一个并行化、开源的工具，用于研究 Hubbard 模型中涌现的反铁磁性与超导性等现象。

提出的方法

• 该程序实现了完整的抛物线形式，自洽求解耦合的 Bethe-Salpeter 方程、Schwinger-Dyson 方程与抛物线方程。
• 采用两级核近似，通过将双粒子顶点函数近似为变量更少的函数，降低其复杂度。
• 应用高频调节以建模顶点函数的渐近频率依赖性，确保在有限频率区间内的收敛性。
• 使用 MPI 进行并行化，采用集体通信（一对一与全对全）以最小化点对点消息通信，提升可扩展性。
• 在数值实现中显式施加了平移不变性、SU(2) 自旋对称性与时间反演对称性。
• 程序计算不同通道（磁性、配对、电荷）中 Bethe-Salpeter 核的最大特征值，以探测即将发生的相变。

实验结果

研究问题

• RQ1 如何高效求解 Hubbard 模型在动量空间与 Matsubara 频率空间中大型双粒子顶点函数的抛物线方程？
• RQ2 在有限频率与动量网格下求解抛物线方程时，需要哪些数值近似以确保收敛性与物理准确性？
• RQ3 如何利用 Bethe-Salpeter 方程的特征值探测并表征强关联体系中的相变？
• RQ4 高频正则化在 Matsubara 频率空间中稳定抛物线方程求解过程中起到何种作用？
• RQ5 该程序的并行化策略在最多 6×6 的集群规模下，跨多个核心的扩展性能如何？

主要发现

• 两级核近似与高频正则化对于在有限频率区间内实现抛物线方程的收敛至关重要。
• 在半满的正方晶格上，磁性通道的特征值 λ 随温度降低而单调增加并趋近于 1，表明出现反铁磁不稳定性。
• 在配对通道中，高温下 s 波配对占主导，而低温下 d 波配对（特别是 $d_{x^2-y^2}$）成为主导。
• 磁性通道的最大特征值出现在 ${\bf q} = (\pi, \pi)$，与正方晶格上的反铁磁序一致。
• 该程序成功计算了 Bethe-Salpeter 方程的特征函数，可通过动量空间波函数的模式识别出配对对称性。
• Victory v1.0 已可实现对 6×6 集群上任意电子数的 Hubbard 模型的研究，即使未经进一步优化亦可行。