[论文解读] Improved Lieb-Oxford bound on the indirect and exchange energies

Mathieu Lewin, Élliott H. Lieb+1 more|arXiv (Cornell University)|Mar 23, 2022

Advanced Chemical Physics Studies被引用 4

一句话总结

该论文通过引入基于最优传输对偶性和数值优化的新型变分方法，改进了多体量子系统中交换能与间接能量的Lieb-Oxford界。作者推导出Lieb-Oxford常数的新上界1.5765，显著优于此前最佳估计值1.64，并为Hartree-Fock态建立了更紧的上界1.25，证明交换能远低于完整的间接库仑能。

ABSTRACT

The Lieb-Oxford inequality provides a lower bound on the Coulomb energy of a classical system of $N$ identical charges only in terms of their one-particle density. We prove here a new estimate on the best constant in this inequality. Numerical evaluation provides the value 1.58, which is a significant improvement to the previously known value 1.64. The best constant has recently been shown to be larger than 1.44. In a second part, we prove that the constant can be reduced to 1.25 when the inequality is restricted to Hartree-Fock states. This is the first proof that the exchange term is always much lower than the full indirect Coulomb energy.

研究动机与目标

• 改进已知的Lieb-Oxford常数上界，该常数约束了密度泛函理论中的间接库仑能。
• 严格建立在Hartree-Fock态限制下交换能的更紧上界。
• 开发一种结合最优传输对偶性与离散测度逼近的数值稳定框架，用于计算Lieb-Oxford常数。
• 通过系统性地细化离散化参数，验证数值方法的收敛性与准确性。

提出的方法

• 作者采用多边际最优传输问题的对偶公式，将间接能量表示为对偶势与耦合测度的函数。
• 他们引入基于一对测度μ和ν的Lieb-Oxford常数的变分表示，并定义一个能量泛函I(μ,ν)，用以从上方界定该常数。
• 通过同心球面表示均匀球形测度，对问题进行离散化，其中径向密度在半径R内以M个点的网格进行近似。
• 在GPU集群上并行应用Broyden-Fletcher-Goldfarb-Shanno（BFGS）优化算法，最小化I(μ,ν)关于离散化测度系数的函数。
• 使用迭代算法改进格林函数的近似，从而提升常数的上界估计。
• 通过改变M（网格点数）、R（截断半径）和K（同心球面数）来测试数值收敛性，结果稳定在10−4精度内。

实验结果

研究问题

• RQ1 对于具有相同粒子的N体系统，Lieb-Oxford常数的最紧可能上界是什么？
• RQ2 当限制在Hartree-Fock态时，能否显著改进Lieb-Oxford界，已知一般情况下该常数的下界为1.44？
• RQ3 使用有限维径向测度离散化，能否在数值上高精度逼近最优耦合测度对(μ,ν)？
• RQ4 离散化参数M、R和K的选择在多大程度上影响计算上界的结果收敛性与准确性？
• RQ5 该数值方法是否足够稳健，能够识别出非平凡的、非均匀的测度对，其结果优于均匀球体情况？

主要发现

• 作者为Lieb-Oxford常数建立了1.5765的新上界，显著优于此前最佳估计值1.64。
• 数值评估在10−4精度内可靠收敛，且随着M和R值增加，结果趋于稳定。
• 通过BFGS优化获得的最优测度对(μ,ν)是非均匀的：μ在单位球内稍靠内处具有平滑支持，而ν由三个位于球面上的狄拉克型分量组成，包括一个位于半径1处的分量。
• 对于Hartree-Fock态，交换能的上界为常数1.25，严格小于一般情况的下界1.44，证明交换能远小于完整的间接库仑能。
• 方法在M和R上均表现出收敛性，随着R增大，上界减小；随着M增加，上界上升，证实了数值近似方案的有效性。
• 利用均匀球体情况下Ψμν的精确解析形式，实现了高精度基准测试，验证了数值框架的可靠性。