[论文解读] Improved estimator for numerical renormalization group calculations of the self-energy

Fabian B. Kugler|arXiv (Cornell University)|Feb 8, 2022

Quantum and electron transport phenomena参考文献 46被引用 22

一句话总结

本文提出了一种用于数值重整化群（NRG）计算的新自能估计算法，通过结合一阶与二阶运动方程，消除了退格林函数虚部中的伪影。该方法确保了高频极限下的精确归一化，抑制了低能区的非物理振荡与波动，并在更少保留态数下实现收敛，显著提升了多轨道DMFT计算的精度与效率。

ABSTRACT

We present a new estimator for the self-energy based on a combination of two equations of motion and discuss its benefits for numerical renormalization group (NRG) calculations. In challenging regimes, NRG results from the standard estimator, a ratio of two correlators, often suffer from artifacts: the imaginary part of the retarded self-energy is not properly normalized and, at low energies, overshoots to unphysical values and displays wiggles. We show that the new estimator resolves the artifacts in these properties as they can be determined directly from the imaginary parts of auxiliary correlators and do not involve real parts obtained by Kramers-Kronig transform. Furthermore, we find that the new estimator yields converged results with reduced numerical effort (for a lower number of kept states) and thus is highly valuable when applying NRG to multiorbital systems. Our analysis is targeted at NRG treatments of quantum impurity models, especially those arising within dynamical mean-field theory, but most results can be straightforwardly generalized to other impurity or cluster solvers.

研究动机与目标

• 解决NRG自能计算中持续存在的伪影问题，如退格林函数虚部中的非物理振荡与波动。
• 解决标准NRG实现中实部高频渐近行为的归一化问题，该问题可能被违反。
• 通过在NRG计算中减少保留态数，降低精确自能计算的数值成本。
• 为量子点模型中的自能估计提供一种稳健且可推广的框架，尤其适用于DMFT和多轨道体系。
• 提升NRG作为强关联材料与晶格模型中点系问题求解器的可靠性与效率。

提出的方法

• 提出一种基于一阶与二阶运动方程联合应用的新自能估计算法，避免依赖实部的Kramers–Kronig变换。
• 推导出涉及三个关联函数（G, F, I）的公式，直接从辅助关联函数确定自能的虚部。
• 通过构造确保实部自能的高频渐近行为精确满足，利用fdm NRG中的谱权守恒。
• 直接使用辅助关联函数F与I的虚部，避免通过Kramers–Kronig变换获取实部，后者易受数值不准确影响。
• 以保持对称性与数值稳定性的方式实现新估计算法，尤其在矩阵值情形下。
• 在多个模型中验证该方法：单轨道Anderson模型，以及在一、二、三轨道Hubbard模型中的DMFT设置。

实验结果

研究问题

• RQ1 能否推导出一种新的自能估计算法，避免标准NRG计算中观察到的ImΣν在低能区的非物理振荡与波动？
• RQ2 新估计算法是否能在具有挑战性的多轨道区域中，精确保持ReΣ高频渐近行为的归一化？
• RQ3 与标准估计算法相比，新方法是否能在更少保留态数下实现自能的收敛结果？
• RQ4 在粒子--hole对称与强关联区域，特别是在低能区，新估计算法表现如何？
• RQ5 新估计算法在多轨道体系中矩阵值关联函数的推广程度如何？

主要发现

• 新估计算法ΣIFG消除了低能区ImΣν中的非物理振荡，而标准NRG估计算法ΣFG中曾存在此类现象。
• 自能的虚部−ImΣIFGν呈现出清晰、非负的抛物线形状，其顶点精确位于原点，符合因果律与费米液体理论。
• 在相同精度要求下，新估计算法收敛所需的保留态数显著更少——例如，Nkp = 500即可实现ImΣIFGν的收敛，而ΣFG需Nkp = 2000才能接近收敛。
• 新估计算法精确保持了ReΣ的高频渐近行为，而标准方法可能违反此解析约束。
• ImΣIFG(0)的值相比ΣFG提升了数个数量级，在粒子-hole对称情况下趋近于精确值零。
• 该方法通过在较低Nkp下实现高精度结果，显著降低了数值计算量，使其在多轨道DMFT应用中极为高效。