Locally Scaled Self-Interaction Corrected Energy Functionals with Complex Optimal Orbitals

本文提出了一种全变分局部尺度自相互作用修正能量泛函，该泛函利用复最优轨道和基于动能密度的尺度因子，在不同电子密度机制下动态调整修正量，从而提高了对原子、分子及固态系统的预测能力。

要点

大局观：修复有缺陷的地图

想象一下，你正试图通过一张地图来在城市中导航。长期以来，科学家们一直使用一种特定类型的地图（称为 KS-DFT）来预测原子和分子的行为。这张地图非常有用且快速，但它有一个著名的缺陷：它存在“自相互作用误差”。

类比：
把电子想象成一个在拥挤房间里行走的人。在现实中，一个人并不会撞到自己。然而，这张旧地图错误地计算出这个人正在撞向自己，从而创造了一个由其自身重量和存在感构成的虚假“幽灵”。这个幽灵干扰了对原子间结合强度或移动它们所需能量的计算。

前人的尝试：“一刀切”的修复方案

科学家们意识到这个“幽灵”问题需要解决。他们发明了一种被称为自相互作用校正 (SIC) 的修正方法。

• 完全修复 (SIC)： 想象一下告诉地图：“把幽灵完全删掉。”如果房间里只有一个人（单个电子），这会运作得非常完美。地图变得无懈可疵。
• 半数修复 (1/2 SIC)： 但当房间里有很多人（许多重叠的电子）时，完全删除幽灵会导致地图过度反应，产生偏差。因此，科学家们尝试只删除一半的幽灵。这在某些方面（如分子间的粘附程度）表现良好，但在其他方面（如原子在激发态或远离彼此时的行为）却表现不佳。

问题在于，科学家们必须做出选择：要么使用“完全修复”（适用于单电子，但不利于人群），要么使用“半数修复”（适用于人群，但不利于单电子）。他们无法两者兼得。

新的解决方案：一个“智能调光开关”

这篇论文介绍了一种名为局部缩放自相互作用校正 (LSSIC) 的新方法。

类比：
与其使用一个全局开关，为整个房间统一开启或关闭（或开启一半）幽灵校正，作者构建了一个智能调光开关，它能根据你在房间里的位置自动进行调节。

• 在孤立区域（低密度）： 如果一个电子是孤身一人的（比如氢离子中的单个电子），调光器会将校正强度调至 100%。幽灵被完全移除，从而得到完美的结果。
• 在拥挤区域（高密度）： 如果电子聚集在一起并发生重叠，调光器会将校正强度调低甚至关闭。这可以防止地图过度校正，从而避免出现异常情况。

这个“调光器”由一个数学函数（称为 $z(r)$）控制，该函数通过观察电子的“交通密度”来运作。它准确地知道何时应用完全修复，以及何时应该克制。

秘密武器：“复数”轨道

论文还提到了使用**“复数最优轨道 (Complex Optimal Orbitals)”**。
类比：
想象电子不仅仅是在直线上行走，而是在进行 3D 螺旋运动。以前的地图试图将这种 3D 螺旋简化为 2D 线条以简化数学计算，但这丢失了一些细节。新方法则拥抱了完整的 3D 螺旋（即“复数”特性）。这使得“智能调光器”能够更清晰地观察交通模式，并以更高的精度调整校正。

他们测试了什么？

作者在几种场景下测试了这种新的“智能地图”：

1. 单电子（氢离子）：
   • 结果： 新方法运作完美。它正确预测了单个电子的行为，就像旧的“完全修复”一样准确，但没有副作用。
2. 单个原子（碳、氮、氧）：
   • 结果： 新方法在预测获取额外电子所需的能量（电子亲和能）方面表现出色。在预测移除电子的难度（电离能）方面，虽然没有达到革命性的突破，但依然非常精确。
3. 分子（原子对）：
   • 结果： 当两个原子结合（如两个碳原子或两个氮原子）时，新方法非常准确地预测了化学键的强度和距离。它的表现通常优于“半数修复”，并避免了“完全修复”带来的误差。

核心结论

这篇论文为科学家用于模拟化学和材料的工具提供了一次重大升级。通过创建一个结合了复数轨道（3D 螺旋）的局部缩放函数（智能调光器），他们构建了一种既能：

• 在电子单独存在时完美修复“幽灵”误差；
• 又在电子拥挤时不会过度校正；
• 且适用于单原子、分子及固体材料的方法。

这就像是从一个强迫你在两条糟糕路线中做选择的地图，升级到了一个能根据你遇到的每种具体交通状况自动寻找最佳路线的 GPS。

技术摘要

技术摘要：基于复数最优轨道的局部定标自相互作用校正能量泛函

问题陈述
Kohn–Sham 密度泛函理论 (KS-DFT) 是电子结构计算的标准方法，但存在自相互作用误差 (SIE)，这种误差源于自库仑相互作用与自交换相关相互作用的相互抵消不完全。该误差导致在描述局域化电子态、里德堡激发、过渡金属磁性和化学键能方面出现显著的不准确性。虽然 Perdew-Zunger 自相互作用校正 (PZ-SIC) 可以逐轨道消除 SIE，但其完全变分实现往往会在轨道重叠的高密度区域过度校正能量。相反，通过 $1/2$ 因子进行经验定标的校正方法 ($1/2$SIC) 虽然改善了许多性质，但无法恢复准确描述里德堡态和局域电荷分布所需的 $-1/r$ 电势依赖性。挑战在于开发一种框架，能够在必要时（例如单电子极限或孤立轨道）应用全 SIC，同时在离域区域减少校正以避免过度校正，且所有这些都要在能够容纳复数最优轨道的完全变分形式体系内实现。

方法论
作者在开源代码 GPAW 中实现了一种完全变分的局部定标自相互作用校正 (LSSIC) 能量泛函。其核心创新在于引入了一个空间依赖的局部定标函数 $z_\sigma(\mathbf{r})$，它取代了 PZ-SIC 形式中的全局定标因子 ($a$)。

• 等轨道指示器： 定标函数源自基于动能密度的等轨道指示器。它被构建在 0 到 1 之间，将 SIC 校正从离域区域（均匀电子气极限）的 0 缩放到单轨道或孤立轨道区域的 1。
• 复数最优轨道： 本工作的关键理论进展是重新推导了 von Weizsäcker 和 Kohn-Sham 动能密度，以处理复数值的最优轨道。局部定标函数 $z_\sigma(\mathbf{r})$ 包含一个依赖于复数相位梯度（通过电流密度 $\mathbf{J}_\sigma(\mathbf{r})$ 表示）的项。
• 变分实现： 通过直接优化方法对所有电子自由度进行能量最小化 ($dE/d\psi^*_k = 0$)。该实现利用了投影缀加平面波 (PAW) 方法的幺正变体，并将波函数视为复数。
• 近似： 作者将 LSSIC 泛函应用于伪密度，并对芯电子采用冻结芯近似。

主要贡献

1. 理论公式： 本文提出了第一个基于复数最优轨道的完全变分 LSSIC 泛函实现。这包括推导了包含复数相位梯度的广义等轨道指示器（公式 5 和 6）。
2. 数值实现： 该方法已在 GPAW 代码中实现，允许在 LSSIC 框架下对复数轨道进行自洽优化。
3. 连接基态与激发态： 该框架旨在弥合需要全 SIC（用于激发态和局域电荷）的方法与需要降低 SIC（用于多电子系统的基态能量学）的方法之间的差距。

结果
作者针对若干体系，将 LSSIC 方法与标准 PBE、$1/2$SIC 以及全 SIC (PZ-SIC) 计算进行了评估：

• 氢二离子阳离子 ($H_2^+$)： 作为单电子系统，LSSIC 正确恢复了精确的 SIC 极限 ($z \approx 1$)。它给出了正确的解离极限和键能，与全 SIC 结果一致，而 PBE 失效且 $1/2$SIC 低估了结合能。
• 原子电子亲和力 (EA) 和电离能 (IE)： 对于 C、N 和 O 原子，包含相位项的变分 LSSIC 提供了与实验电子亲和力最好的吻合度，显著优于 PBE 和 $1/2$SIC。对于电离能，LSSIC 通常与 PBE 匹配或略有改进，但并未在所有元素上一致地优于全 SIC。复数相位项的引入略微改善了 N 和 O 的 IE 预测。
• 二聚体键能与键长： 对于 $C_2$、$N_2$ 和 $O_2$，LSSIC 泛函预测的键能和键长与实验值偏差在百分之几以内。
  • 对于 $C_2$ 和 $O_2$，LSSIC（无论是否包含相位项）产生的键能比 PBE、$1/2$SIC 或全 SIC 更接近实验值。
  • 对于 $N_2$，LSSIC 相比 $1/2$SIC 略微高估了键能，但势能面形状保持相似。
  • “单次”（One-shot）LSSIC（即定标函数固定自 PBE 计算）在 $C_2$ 和 $O_2$ 键能预测上表现出惊人的良好效果，但在整体上逊于完全变分方法。

意义与主张
论文声称 LSSIC 框架代表了完全变分 SIC 泛函的一个重要里程碑。通过引入能够自然适应电子密度和轨道重叠的局部定标函数，该方法在严格需要全 SIC 的区域（如单电子极限）保留了全 SIC 的优势，同时减轻了在高密度离域区域的过度校正。作者指出，这种方法提高了 DFT 对原子和分子性质的预测精度，特别是电子亲和力，并提供了一个适用于原子、分子和固态系统的通用框架。他们强调，该方法是改进异质催化计算的重要一步，因为在这些领域，既需要准确描述分子吸附质，也需要准确描述固体基底。这项工作并不声称具有超越所有现有方法的普适优越性，而是证明了 LSSIC 提供了一个平衡的选择，在特定语境下可以优于标准的半局域泛函和固定定标方案。