[论文解读] ana_cont: Python package for analytic continuation

Josef Kaufmann, Karsten Held|arXiv (Cornell University)|May 24, 2021

Advanced Chemical Physics Studies被引用 8

一句话总结

本论文提出 ana_cont，一个用于费米子和玻色子格林函数解析延拓的 Python 软件包，采用帕德近似和最大熵方法。该工具可实现从虚时或马苏巴拉频率到实频率的稳健、灵活且用户友好的解析延拓，支持图形界面和程序化访问，适用于可复现的工作流程及准确谱函数重构的超参数优化。

ABSTRACT

We present the Python package ana_cont for the analytic continuation of fermionic and bosonic many-body Green's functions by means of either the Pade approximants or the maximum entropy method. The determination of hyperparameters and the implementation are described in detail. The code is publicly available on GitHub, where also documentation and learning resources are provided.

研究动机与目标

• 提供一个可靠、开源的工具，用于将量子多体关联函数从虚时或马苏巴拉频率解析延拓至实频率。
• 实现帕德近似和最大熵方法，具备稳健的超参数优化和数值稳定性。
• 为非专家（通过图形界面）和高级用户（通过库接口）提供对费米子和玻色子系统进行解析延拓的支持。
• 通过文档化的教程、测试用例以及与科学 Python 生态系统的集成，实现可复现且可扩展的工作流程。
• 通过解决解析延拓中对噪声的敏感性和参数依赖性问题，促进准确的谱函数重构。

提出的方法

• 该软件包使用 Cython 编译的核心函数实现帕德插值，以高效地在复平面上对格林函数进行有理逼近。
• 对于最大熵方法，采用贝叶斯推断，结合默认模型和超参数 α 与 b，以优化谱函数的后验概率。
• 该方法基于最大熵原理，选择与噪声马苏巴拉数据最一致的最可能谱函数。
• 引入预模糊化处理以提升最大熵方法的稳定性，尤其适用于噪声较大或病态条件的数据。
• 代码支持费米子和玻色子格林函数，对马苏巴拉频率和虚时关联函数有专门的处理机制。
• 用户可通过图形用户界面快速分析，或通过 Python 库接口以程序化方式实现自定义工作流。

实验结果

研究问题

• RQ1 如何使量子多体格林函数的解析延拓既准确又对非专家友好？
• RQ2 在存在噪声和有限数据分辨率的情况下，最大熵方法的最优超参数（α, b）是什么？
• RQ3 帕德近似与最大熵方法在从噪声马苏巴拉数据重构谱函数方面有何比较优势？
• RQ4 预模糊化在多大程度上提升了最大熵解析延拓的稳定性和准确性？
• RQ5 能否开发一个统一、可扩展的 Python 软件包，支持费米子和玻色子系统，并具备一致且文档化的使用流程？

主要发现

• ana_cont 软件包成功实现了帕德和最大熵方法的解析延拓，具备高度的数值稳定性和可复现性。
• 通过系统的 χ² 最小化，确定了最大熵方法的最优超参数 α 和 b，结果在论文图中可视化展示。
• 预模糊化显著提升了最大熵方法的收敛性和鲁棒性，尤其在噪声数据下表现突出，如 SrVO₃ 案例研究所示。
• 图形界面使非专家可仅通过极少代码完成解析延拓，而库接口则允许高级用户自定义工作流。
• 该软件包支持费米子和玻色子格林函数，与标准量子多体数据格式和数值方法完全兼容。
• 测试用例和教程（包括 SrVO₃ 的详细示例）展示了该软件包从噪声马苏巴拉数据中恢复已知谱特征的能力。