扩展哈伯德链中条件算符格林函数方法谱截断的精确对角化基准

《Physica B: Condensed Matter》：Exact-diagonalization benchmark of spectral truncation relevant to the conditional-operator Green-function method in extended Hubbard chains

　　

•我们使用条件算符格林函数方法对长度至L=6的扩展哈伯德链进行了基准测试，证实了尽管Fock空间维度以4^L增长，但数量级为多项式增长的极点数即可捕捉主要的谱权重（N90 ?w/dim ? 从L=2到L=6由0.19降至0.0027）。•当L>2时，基于权重的截断被证明显著比能量截断更高效——在L=4时，按权重排序仅需6个极点即可达到90%的谱权重，而按能量排序需要11个极点——这一差异归因于基态Fock空间参与比的增长。•提供了L=2,4,6的开放哈伯德链在t=1, U=4, V=1, ??=2下G??↑(ω)的精确莱曼分解，包括L=4的完整Fock基和条件算符表，作为权威参考。•祖巴列夫/哈伯德-I平均场解耦在这些参数下被证明完全遗漏了关联物理：由于??=0，其对G??↑(ω)的预测与非相互作用链完全相同，而精确光谱在L=4时有38个极点，在L=6时有237个极点。•我们概述了条件算符方法的实用截断方案：基态应按矩阵元权重而非仅按能量分母进行排序，权重分布由基态的构型振幅决定，而非由扇区级能量排序决定。•完整的、经过验证的Python和Mathematica精确对角化代码作为补充材料提供。引言范围与免责声明。本文未展示任何可标度的格林函数求解器。本文是一篇精确对角化基准论文，旨在确定促使条件算符方法的谱结构能否在二聚体之外得以保持。所有L≤6的结果均通过完整稠密对角化获得；L=8的结果使用稀疏Lanczos方法（ARPACK）获得，如补充代码中所实现的。向更大系统尺寸的推广以及实用条件算符截断实现方案的构建留待未来工作。特别是，虽然基准测试表明在所研究的系统尺寸范围内，相关的莱曼贡献数相对于完整Fock空间是小的，但四个系统尺寸不足以确定N_w90(L)的渐近标度形式。Bowen、Williams和Mancini引入了一种基于Masson定理的格林函数方法，其中留数矩阵元在条件（位移）算符Φ^(k)_{α→β}=I_β c_k I_α的基中展开，其中I_α是到Fock态|α?上的幂等投影子[1], [2], [3]。Masson定理在其规定的数学条件下提供了收敛性，包括在所适用哈密顿量类别中任意强耦合情况[1], [2]。截断近似解是各自良好定义的且可系统改进的，这与祖巴列夫解耦[4]或任意连分数截断形成对比。有限截断在一般情况下并非精确的；精确终止可在相关动力学子空间闭合时发生，但在此处研究的较长链中未假设任何此类闭合。一项伴随的二聚体计算在具体层面上确立了这一形式体系：它枚举了两点扩展哈伯德二聚体的所有32个条件湮灭算符及其约旦-温格相位，获得了精确基态和G??↑(ω)的精确莱曼分解，并表明仅保留|λ|最小的两个极点即可捕捉90%的谱权重，其余10%在下一壳层恢复[5]。该工作还基于跳跃连通性给出了数量级估计，表明在能量截断Λ_max下存活的条件算符数量可能随链长L多项式增长，而完整Fock空间以4^L增长，并将明确的长链计算标记为直接的下一步。二聚体结果出现在该特定哈密顿量、参数和可观测量下的六维半填充扇区中，且对于该参数组，主要谱权重恰好与最低|λ|极点重合。先验上并不明显两极点、90%集中于两极点的图像是该方法的普遍特征还是最小可能系统尺寸的人为产物。本文直接回答了这一问题，通过对相同参数下L=4,6,8的开放哈伯德链进行精确对角化，并报告了二聚体直觉成立之处以及需要修正之处。Masson/Feenberg条件算符方法与其他几种非微扰的关联电子格林函数方法并列。在一维情况下，DMRG是在大得多的系统尺寸下一维哈伯德系统基态和动力学性质最准确的既定数值方法之一[6], [7]。Mancini和Avella的复合算符方法[8]从在运动方程下闭合的复合场基构建格林函数，与本文明使用的条件算符方法共享投影技术理念，但在基的选择和闭合方面有所不同。最近，Lira等人[9]发展了关联系统格林函数的相关投影子方法，为截断运动方程层级中的挑战提供了独立的视角。他们的方法与本文明使用的条件算符方法共享投影技术理念，但在基的选择和闭合方面有所不同。与他们结果的比较突显了在尝试大规模实现之前对小系统进行截断方案与精确对角化基准测试的重要性。少量位点数的有限哈伯德系统也可以在超冷原子光晶格中实现和探测[10], [11]。章节摘要哈密顿量我们研究的扩展哈伯德链为H = -t Σ_{σ=↑,↓} Σ_{i=1}^{L-1} (c?_{iσ} c_{i+1,σ} + h.c.) + U Σ_{i=1}^L n_{i↑} n_{i↓} + V Σ_{i=1}^{L-1} n_i n_{i+1} - μ Σ_{i=1}^L n_i，其中n_i = n_{i↑} + n_{i↓}，采用开放边界条件，并使用典范约旦-温格轨道排序1↑ < 1↓ < 2↑ < 2↓ < ? < L↑ < L↓。对于L=2，这正好是二聚体计算[5]的哈密顿量；以下结果使用相同参数t=1, U=4, V=1, μ=2（能量以t为单位），使得L=2的情况重现该计算报告的結果作为正确性检查（第2.3节）。精确对角化与莱曼表示因为条件算符枚举与截断标度在展示谱结果之前，我们提供条件算符在简单能量截断下存活数量的单独组合估计。这与第4节报告的莱曼极点数不同；它计算的是Liouville空间基中的算符，而非谱极点。以下推导建立了固定能量截断下条件算符基的多项式标度。对于固定能量截断Λ，条件湮灭算符Φ^(α→β)的数量表2总结了L=2,4,6,8下G??↑(ω)的精确谱分解。比较了两种截断排序：按权重排序的数值分辨莱曼贡献（对目标精度而言是物理相关的排序），以及按|λ|递增排序的莱曼能量排序（类似于二聚体计算中用于L=2截断分析的能量截断排序）。对于L=2，两个主要莱曼贡献的总权重为9/20 + 9/20 = 9/10 = 0.90，因此它们满足与一阶祖巴列夫解耦参考的比较二聚体计算中与一阶祖巴列夫解耦的比较促使了此处所研究较长链的相应基准测试。本节按照二聚体计算的规范将同样的比较推广到L=4,6,8。该比较仅作为谱分辨能力的基准；并非对两种方法计算复杂度或形式收敛性质的比较。N_w90的增长缓慢且亚指数在所研究的四个系统尺寸范围内，捕捉90%谱权重所需的数值分辨莱曼贡献数相对于完整Fock空间维度仍然很小：2→6→11→22，对应16→256→4096→65536。比值N_w90/dim ? 从L=2时的0.125降至L=4时的0.023、L=6时的0.0027和L=8时的0.000336。这与所研究系统尺寸范围内的亚指数增长一致。然而，精确的渐近标度形式无法从四个系统尺寸确定；结果摘要在匹配已报告二聚体基准的参数下，对扩展哈伯德链进行至L=8的精确对角化，表明：1.捕捉90%谱权重所需的数值分辨莱曼贡献数随L=2,4,6,8分别为2→6→11→22增长，而完整Fock空间维度从16增至256、4096和65536。所要求的贡献数在所研究的四个系统尺寸范围内始终相对于完整Fock空间维度很小，CRediT作者贡献声明Jay D. Mancini：原始论文撰写、软件、项目管理、方法学、调查、概念化。数据与代码可用性本文报告的所有莱曼极点和权重的完整数据集，连同用于生成它们的Mathematica和Python实现，已存档于Zenodo：J. D. Mancini, Exact-Diagonalization Benchmark of Spectral Truncation Relevant to the Conditional-Operator Green-Function Method in Extended Hubbard Chains, Zenodo (2026)。版本DOI：https://doi.org/10.5281/zenodo.22800578概念DOI：https://doi.org/10.5281/zenodo.22800577许可证：CC-BY-4.0。存档包含两个文件：•利益冲突声明作者声明他们没有已知的竞争性经济利益或个人关系可能对本论文报告的工作产生影响了影响。Jay D. Mancini