[论文解读] A self-consistent ground-state formulation of the first-principles Hubbard U parameter validated on one-electron self-interaction error

Glenn Moynihan, Gilberto Teobaldi+1 more|arXiv (Cornell University)|Apr 26, 2017

Advanced Chemical Physics Studies被引用 5

一句话总结

本文提出了一种自洽的、基于基态的从头算Hubbard U参数形式，采用完全弛豫的约束密度，实现了精确且计算高效的DFT+U计算。结果表明，通过U^(2)自洽方案，DFT+U能精确校正H₂⁺中的单电子自相互作用误差，验证了该方法在从头算一致性及不同U值系统间能量可比性方面的有效性。

ABSTRACT

In electronic structure methods based on the correction of approximate density-functional theory (DFT) for systematic inaccuracies, Hubbard $U$ parameters may be used to quantify and amend the self-interaction errors ascribed to selected subspaces. Here, in order to enable the accurate, computationally convenient calculation of $U$ by means of DFT algorithms that locate the ground-state by direct total-energy minimization, we introduce a reformulation of the successful linear-response method for $U$ in terms of the fully-relaxed constrained ground-state density. Defining $U$ as an implicit functional of the ground-state density implies the comparability of DFT + Hubbard $U$ (DFT+$U$) total-energies, and related properties, as external parameters such as ionic positions are varied together with their corresponding first-principles $U$ values. Our approach provides a framework in which to address the partially unresolved question of self-consistency over $U$, for which plausible schemes have been proposed, and to precisely define the energy associated with subspace many-body self-interaction error. We demonstrate that DFT+$U$ precisely corrects the total energy for self-interaction error under ideal conditions, but only if a simple self-consistency condition is applied. Such parameters also promote to first-principles a recently proposed DFT+$U$ based method for enforcing Koopmans' theorem.

研究动机与目标

• 解决DFT+U计算中Hubbard U参数自洽性问题的未解难题。
• 提出一种仅依赖基态密度的从头算Hubbard U，确保不同U值系统间能量可比性。
• 在直接总能量最小化框架（包括线性标度求解器）中，实现精确且计算便捷的DFT+U计算。
• 验证该方法在理想条件下校正单电子自相互作用误差的能力。
• 为通过从头算U参数实现DFT+U中的Koopmans定理提供基础。

提出的方法

• 将线性响应方法中的Hubbard U重新表述为完全弛豫的约束基态密度形式，使U成为基态密度的隐式泛函。
• 通过总能量最小化定义U为基态性质，确保与依赖基态优化的DFT求解器的一致性。
• 应用U^(2)自洽准则，即在保持基态密度一致的前提下，对总能量关于U进行最小化。
• 采用旋转不变的DFT+U能量项：E_U = ∑_{I,σ} (U^I / 2) Tr[n^Iσ − n^Iσ n^Iσ]，其中n^Iσ为投影密度矩阵。
• 利用约束密度泛函理论计算约束密度对微扰的响应，实现U计算而无需激发态信息。
• 在H₂⁺离子上验证该方法，该体系为自相互作用误差纯粹为单电子且可解析处理的最小体系。

实验结果

研究问题

• RQ1 是否能完全在DFT的基态密度框架内构建自洽的Hubbard U参数，确保不同系统间能量可比性？
• RQ2 U^(2)自洽方案是否能产生一种从头算的、物理一致的Hubbard U，从而在总能量中校正自相互作用误差？
• RQ3 该形式的DFT+U能否在H₂⁺等最小体系中精确校正单电子自相互作用误差？
• RQ4 与依赖激发态或非基态量的传统线性响应方法相比，该方法有何优势？
• RQ5 DFT+U能量表达式中的线性项在确保整数电子填充时正确行为方面起什么作用？

主要发现

• 所提出的自洽基态U形式确保DFT+U总能量在不同U值系统间具有可比性，因为U和密度均为基态性质。
• U^(2)自洽方案在H₂⁺离子中实现了对单电子自相互作用误差的精确校正，理想条件下DFT+U总能量与精确结果完全一致。
• 该方法验证了理论预期：当正确U值与自洽方案结合时，DFT+U可完全校正自相互作用误差。
• 该形式为DFT+U₁+U₂方法提供了从头算基础，该方法将自相互作用误差校正扩展至最高占据本征值，近似实现Koopmans定理。
• 本研究表明，标准线性响应方法中的U并非基态性质，而所提出的重构形式通过将U锚定于基态密度解决了这一问题。
• 该方法简化了自洽方案的分析，并识别出U^(2)是基于基态密度的变分线性响应DFT+U的最优准则。