GEM: An implementation of the ghost-Gutzwiller approximation for simulating interacting quantum systems

技术摘要：“GEM：一种用于模拟相互作用量子系统的 ghost-Gutzwiller 近似实现”

问题陈述
强关联电子模型的非微扰计算面临着物理保真度、计算成本与可获取观测范围之间的根本权衡。虽然动力学平均场理论（DMFT）为处理强局部关联提供了严谨的框架，但由于其需要重复求解涉及频率相关格林函数的相互作用量子杂质问题，随着轨道数增加、温度降低或研究对称性破缺相时，其计算成本变得难以承受。相反，传统的 Gutzwiller 近似（GA）虽然计算效率高，但无法重现非相干的高能哈伯德带（Hubbard bands）。ghost-Gutzwiller 近似（ghost-GA）通过增加变分辅助空间来恢复高能光谱特征和具有 DMFT 特性的总能量，提供了一个折中方案，且仅需适度扩大辅助空间。然而，在此项工作之前，尚无开源软件实现了基于自由能泛函的有限温度 ghost-GA 公式，该公式统一了零温 ghost-GA、其有限温度扩展以及在无限辅助态极限下的 DMFT。

方法论
本文介绍了 GEM（Ghost Embedding Method），这是一个实现了量子嵌入格式下 ghost-Gutzwiller 近似的开源 Python 软件包。该方法基于动力学自由能泛函，其驻点条件通过有效哈密顿量的静态或热期望值而非频率相关的格林函数来表达。

关键理论与算法组件包括：

• Ghost 嵌入格式： 该方法为每个物理轨道引入 $B$ 个辅助“ghost”浴模式。自能和杂化函数通过涉及矩阵 $\{R_i, \Lambda_i, D_i, \Lambda^c_i\}$ 的极点展开进行参数化。
  • $B=1$ 回归标准的零温 GA。
  • $B \to \infty$ 回归 DMFT 极限。
  • 中间值（如 $B=3$）通过系统性改进 GA，以较低的 DMFT 成本捕捉高能特征和热力学性质。
• 自洽循环： 解法通过满足鞍点方程（公式 17–20）获得，该方程使三个辅助问题之间的热一阶密度矩阵相匹配：
  1. 代表晶格的准粒子哈密顿量（$\hat{H}_{qp}$）。
  2. 代表杂质的相互作用嵌入哈密顿量（$\hat{H}_{emb}$）。
  3. 用于拟合的非相互作用嵌入哈密顿量（$\hat{H}_{0,emb}$）。
• 有限温度处理： 与需要求解频率相关量的 DMFT 不同，GEM 的自洽性依赖于匹配静态热密度矩阵。对于有限温度，变分参数的更新涉及一个非线性最小二乘拟合程序（密度矩阵拟合），通过近端惩罚项（proximal penalty）来处理规范自由度和低温下的数值平坦性问题。
• 软件架构： GEM 采用模块化设计，将晶格（准粒子）问题与局部片段（嵌入）问题分离。它包含一个内部精确对角化求解器（SimpleED），并为外部求解器（如 DMRG、张量网络）提供通用接口（gemSolver）。它支持基于 MPI 的 $k$ 点求和与对称性扇区并行化。

核心贡献

1. 开源实现： GEM 是第一个公开可用的、维护中的开源软件包，实现了基于自由能泛函的有限温度 ghost-GA 公式。
2. 系统性可提升性： 该软件包允许用户通过改变 ghost 轨道的数量 ($B$) 来系统地调节精度，从而架起廉价的静态 Gutzwiller 方法与完全动力学 DMFT 之间的桥梁。
3. 计算效率： 通过避免频率相关的杂质求解器并依赖热密度矩阵，GEM 为研究多轨道晶格模型、对称性破缺相和有限温度热力学提供了一个高效的框架。
4. 互操作性： 其设计与 TRIQS 生态系统集成，用于模型构建，并允许轻松引入外部杂质求解器。

结果与验证
作者通过三个原型示例验证了实现效果：

1. 贝特晶格热力学： 单轨道哈伯德模型（$U/D=2.0$）的有限温度计算表明，总能量、熵和双占据率随温度平滑收敛。结果显示，$B=3$ 的计算结果与有限温度预期一致且准确。
2. 方格子磁相图： 代码计算了方格子哈伯德模型的反铁磁序参量。使用 $B=3$ 时，该方法正确捕捉了 BCS-BEC 交叉行为，并重现了与既定 DMFT 文献（参考文献 [42]）一致的临界温度 ($T_c$) 和相边界。
3. 多轨道 Hund 物理： 贝特晶格上三轨道 Hubbard-Kanamori 模型的零温计算展示了准粒子权重 ($Z$) 随相互作用强度的演化。结果表明，将 $B$ 从 1 增加到 3 显著改善了对 Hund 物理（$Z$ 的重整化）的描述，使其与 DMFT 基准（参考文献 [44, 45]）更加吻合。

意义与主张
论文声称 GEM 在量子嵌入方法层级中占据了一个至关重要的中间地带。它为研究静态和热观测量的“低成本嵌入方法”提供了支持，在较小的浴尺寸（$B=3$）下，其与 DMFT 的一致性已经非常出色，且无需承担完整的动力学杂质求解的计算开销。此外，收敛的 ghost-GA 解可以作为 DMFT 计算的高质量起点。

作者将 GEM 定位为开发或应用量子嵌入方法的科研人员的工具，能够实现跨模型、求解器、温度和近似水平的可重复比较。该软件旨在具有可扩展性，未来计划与 DMRG/MPS 求解器对接，引入迟滞相互作用，并连接至 ab initio 电子结构代码（例如 DFT+ghost-GA）以进行特定材料的计算。这项工作强调，GEM 使 ghost 嵌入方法在模型计算、基准测试和方法论开发方面变得触手可及，填补了缺乏开源有限温度 ghost-GA 实现所留下的空白。