[论文解读] What can Mott insulators teach us about density-functional theory (and vice versa)?

Neemias Alves de Lima, L. N. Oliveira+1 more|arXiv (Cornell University)|May 27, 2002

Physics of Superconductivity and Magnetism被引用 40

一句话总结

该论文表明，通过在交换-关联泛函中引入导数不连续性，基于Bethe Ansatz的局域密度近似（BA-LDA）在密度泛函理论中能正确预测一维 Hubbard 模型中的莫特绝缘体带隙。与传统泛函不同，BA-LDA 在所有关联强度下均能精确给出莫特带隙，揭示了关联诱导的不连续性在低维系统中主导于能带结构带隙。

ABSTRACT

We study the Mott insulating phase of the one-dimensional Hubbard model using a local-density approximation (LDA) that is based on the Bethe Ansatz (BA). Unlike conventional functionals, the BA-LDA has an explicit derivative discontinuity. We demonstrate that as a consequence of this discontinuity the BA-LDA yields the correct Mott gap, independently of the strength of the correlations. A convenient analytical formula for the Mott gap in the thermodynamic limit is also derived. We find that in one-dimensional quantum systems the contribution of the discontinuity to the full gap is more important than that of the band-structure gap, and discuss some consequences this finding has for electronic-structure calculations.

研究动机与目标

• 研究密度泛函理论（DFT）在一维 Hubbard 模型莫特绝缘相中的适用性。
• 解决标准 LDA 泛函因缺乏导数不连续性而无法描述莫特带隙的问题。
• 评估基于 Bete Ansatz 的局域密度近似（BA-LDA）泛函，该泛函显式包含不连续性，是否能准确描述莫特带隙。
• 比较 Kohn-Sham 带隙与 xc 不连续性在决定低维系统中总多体带隙中的贡献。
• 评估这些发现对量子线和纳米管中第一性原理电子结构计算的启示。

提出的方法

• 本研究采用基于 Bete Ansatz 的局域密度近似（BA-LDA）作为交换-关联（xc）泛函，其来源于无限一维 Hubbard 模型的精确解。
• 利用 Bete Ansatz 解对 xc 能量进行参数化，并显式处理整数粒子数处的导数不连续性。
• 使用精确的多体表达式 Δ = E(N+1) + E(N-1) - 2E(N) 计算总能隙，该表达式在系综 DFT 中有效。
• 通过 Δ = ΔKS + Δxc 将能隙分解为 Kohn-Sham（KS）和 xc 不连续性贡献，其中 Δxc 由 xc 能量的泛函导数计算得出。
• 在有限体系（如 150 根 site 的链）上进行数值计算，并将结果外推至热力学极限。
• 推导出在 U→0 和 U→∞ 极限下能隙的渐近解析表达式，并与 BA-LDA 结果进行比较。

实验结果

研究问题

• RQ1 具有导数不连续性的 DFT 泛函能否正确描述一维 Hubbard 模型中的莫特带隙？
• RQ2 Kohn-Sham 带隙与 xc 不连续性在决定总多体带隙中的贡献如何比较？
• RQ3 基于精确 Bete Ansatz 解的 BA-LDA 泛函是否能在所有相互作用强度下准确再现精确的莫特带隙？
• RQ4 标准 LDA 泛函为何无法描述莫特带隙，而导数不连续性在克服这一失败中起到何种作用？
• RQ5 xc 不连续性对准一维体系（如量子线和纳米管）中电子结构计算有何启示？

主要发现

• BA-LDA 泛函在所有相互作用强度 U 下均能正确再现一维 Hubbard 模型中的精确莫特带隙，包括非解析的 U→0 极限。
• 莫特带隙主要由 xc 不连续性 Δxc 主导，其贡献超过 Kohn-Sham 带隙 ΔKS，尤其在热力学极限下更为显著。
• 在有限体系中，xc 带隙近似按 (1/L)² 缩放，而 KS 带隙按 1/L 线性缩放，表明其具有不同的物理起源。
• 标准 LDA 泛函因缺乏不连续性，导致 Δxc = 0，因此即使在均匀极限下精确，也无法捕捉莫特带隙。
• 解析公式 Δ = U + 4cos(π/β(U)) 准确捕捉了在 U→∞（Δ≈U−4）和 U→0（Δ≈(8/π)√U exp(−2π/U)）极限下的渐近行为。
• 传统 LDA 无法再现带隙的原因并非 LDA 近似本身，而是 xc 泛函中缺乏导数不连续性。