[论文解读] Strong-coupling formula for momentum-dependent susceptibilities in dynamical mean-field theory

Junya Otsuki, Kazuyoshi Yoshimi+2 more|arXiv (Cornell University)|Dec 28, 2018

Physics of Superconductivity and Magnetism参考文献 38被引用 28

一句话总结

该论文提出了一种强耦合近似（SCL公式），用于在动力学平均场理论（DMFT）中计算动量依赖的磁化率，通过利用局域关联简化贝特-萨皮特方程，显著降低了计算成本。该方法在弱耦合到强耦合范围内准确捕捉了如动能交换和RKKY相互作用等物理现象，实现了多轨道体系中自旋和轨道磁化率的高效计算，并提供了三种实用的近似层级。

ABSTRACT

Computing momentum-dependent susceptibilities in the dynamical mean-field theory (DMFT) requires solving the Bethe-Salpeter equation, which demands large computational cost. Exploiting the strong-coupling feature of local fluctuations, we derive a simplified formula that can be solved at a considerably lower cost. The validity and the physical meaning of the formula are confirmed by deriving the effective intersite interactions in the strong-coupling limit, such as the kinetic exchange and RKKY interactions. Furthermore, numerical calculations for single-orbital and multiorbital models demonstrate surprisingly wider applicability including weak-coupling region. Based on this formula, we propose three levels of practical approximations that can be chosen depending on complexity of problems. Simpler evaluations of spin and orbital susceptibilities in multiorbital systems thus become possible within DMFT.

研究动机与目标

• 开发一种计算高效的动量依赖磁化率计算方法，以解决由于贝特-萨皮特方程导致的计算成本过高问题。
• 利用局域涨落的强耦合物理特性，推导出一种简化公式，保留关键物理相互作用，如动能交换和RKKY相互作用。
• 将该方法扩展至多轨道体系，并为实际材料提供可扩展、实用的近似方案。
• 在广泛耦合范围内验证该方法的有效性，包括其在弱耦合区域出人意料的准确性。

提出的方法

• 通过利用局域顶点函数的强耦合特性，对贝特-萨皮特方程进行简化，使其形式类似于RPA，但输入为完全关联的局域磁化率。
• 提出SCL公式：χ(q) ≈ χloc + χloc [X0(q) - X0,loc] χloc，其中χloc为局域磁化率，X0(q)为非相互作用双粒子格林函数。
• 通过验证其正确再现已知的强耦合物理现象（如 Hubbard 模型中的动能交换和周期性 Anderson 模型中的 RKKY 相互作用），验证了公式的有效性。
• 提出三种实用的近似层级：SCL1（仅使用局域χ）、SCL2（使用动量依赖的X0）、SCL3（对u0(iω)采用两极点近似以高效估算χloc）。
• 采用连续时间量子蒙特卡洛和 Hubbard-I 近似作为杂质求解器，结合频率截断和拟合程序以保证数值稳定性。
• 对主导模式函数u0(iω)采用两极点近似，以实现多轨道体系中χloc的高效计算。

实验结果

研究问题

• RQ1 能否从DMFT中的贝特-萨皮特方程导出一种简化公式，以大幅降低计算成本，同时保持关键物理特性？
• RQ2 所提出的SCL公式是否能准确描述Hubbard模型和周期性Anderson模型中已知的强耦合相互作用（如动能交换和RKKY）？
• RQ3 SCL公式在强耦合区域之外的表现如何，特别是在弱耦合和多轨道体系中？
• RQ4 该方法能否系统地分层近似为多个精度和复杂度逐步提升的层级，以适用于实际材料模拟？

主要发现

• SCL公式在从弱耦合到强耦合的整个U/W范围内均能提供准确的动量依赖磁化率，这与预期相反，即此类近似仅在强耦合极限下成立。
• 该方法正确再现了Hubbard模型中的动能交换相互作用和周期性Anderson模型中的RKKY相互作用，证实了其物理有效性。
• 对双轨道模型的数值结果表明，SCL3近似（采用u0(iω)的两极点拟合）在计算成本显著降低的同时，与完整贝特-萨皮特解保持良好一致。
• 对u0(iω)的两极点近似成功捕捉了主导模式结构，拟合参数E± ≈ E±,min，验证了在近似中使用最小激发能的合理性。
• SCL3方案通过χloc ∝ s0|u0|2避免了直接计算s0，从而在无需求解完整杂质问题的情况下实现了高效计算。
• 该公式近似满足两粒子求和规则，并支持将χloc作为可靠输入使用，尽管在弱耦合区域其并非精确值。