[论文解读] Manual for the Flexible DM-NRG code

Örs Legeza, Cătălin Paşcu Moca+3 more|ArXiv.org|Sep 18, 2008

Quantum many-body systems参考文献 4被引用 23

一句话总结

本论文提出了一种灵活的开源 DM-NRG 代码，用于使用密度矩阵数值重整化群方法求解量子杂质模型，支持任意 U(1)、SU(2) 和 Z₂ 对称性。该代码通过用户定义的输入文件，实现对谱函数、期望值和有限尺寸能谱的精确计算，输出结果结构化，便于使用 sfb 和 es 等配套工具进行详细分析。

ABSTRACT

Quantum impurity models describe interactions between some local degrees of freedom and a continuum of non-interacting fermionic or bosonic states. The investigation of quantum impurity models is a starting point towards the understanding of more complex strongly correlated systems, but quantum impurity models also provide the description of various correlated mesoscopic structures, biological and chemical processes, atomic physics and describe phenomena such as dissipation or dephasing. Prototypes of these models are the Anderson impurity model, or the single- and multi-channel Kondo models. The numerical renormalization group method (NRG) proposed by Wilson in mid 70's has been used in its original form for a longtime as one of the most accurate and powerful methods to deal with quatum impurity problems. Recently, a number of new developments took place: First, a spectral sum-conserving density matrix NRG approach (DM-NRG) has been developed, which has also been generalized for non-Abelian symmetries. In this manual we introduce some of the basic concepts of the NRG method and present recently developed Flexible DM-NRG code. This code uses user-defined non-Abelian symmetries dynamically, computes spectral functions, expectation values of local operators for user-defined impurity models. The code can also use a uniform density of states as well as a user-defined density of states. The current version of the code assumes fermionic bath's and it uses any number of U(1), SU(2) charge SU(2) or Z(2) symmetries. The Flexible DM-NRG code can be downloaded from http://www.phy.bme.hu/~dmnrg

研究动机与目标

• 为量子杂质模型提供一种全面、可扩展且用户友好的 DM-NRG 方法实现。
• 支持任意阿贝尔与非阿贝尔对称性（U(1)、SU(2)、Z₂），以提升计算效率与精度。
• 实现谱函数、格林函数及静态算符期望值的高精度计算。
• 提供工具用于分析有限尺寸能谱，并提取临界指数或共形场论数据。
• 通过结构化输入文件、版本控制与在线文档，促进结果可重现性与开发便利性。

提出的方法

• 代码采用灵活的输入文件格式，用于定义模型参数、对称性、哈密顿量与算符。
• 采用威尔逊 NRG 方法结合密度矩阵截断方案（DM-NRG），以提升收敛性与精度。
• 通过对称性通过用户定义的区块与量子数标签动态处理，实现高效的态空间约化。
• 谱函数通过带展宽技术的迭代对角化方法计算，输出按迭代步骤存储。
• 代码支持平坦与用户自定义的 bath 状态密度，允许能量依赖的谱函数计算。
• 提供后处理工具（sfb、es），用于生成展宽谱函数并分析有限尺寸能级结构。

实验结果

研究问题

• RQ1 如何以灵活、可扩展且用户可配置的方式实现 DM-NRG 方法，以适用于任意量子杂质模型？
• RQ2 在 NRG 框架中，如何最优地集成非阿贝尔对称性（如 SU(2) 自旋或电荷对称性）以提升精度？
• RQ3 如何可靠地计算并分解谱函数与格林函数，覆盖 NRG 的偶数与奇数次迭代？
• RQ4 在高对称性与大态空间条件下运行 DM-NRG 模拟的关键计算与内存需求是什么？
• RQ5 如何分析有限尺寸能谱，以提取共形场论数据或费米液体参数？

主要发现

• DM-NRG 代码通过结合偶数与奇数 NRG 迭代的数据，成功实现了高分辨率的谱函数计算。
• 该代码支持任意 U(1)、SU(2) 与 Z₂ 对称性，显著降低了计算成本并提升了精度。
• 工具 sfb 基于 input.dat 中的输入参数，生成展宽谱函数与希尔伯特变换格林函数。
• 工具 es 可对有限尺寸能谱进行详细分析，支持识别共形场论数据与无关算符维数。
• 代码将静态算符与谱算符的期望值分别存储于独立文件中，且偶数与奇数次迭代的数据分开保存。
• m_star 参数标记了态截断开始的迭代步数，该参数被保留以确保谱函数生成的一致性。