Influence of ligand field and correlation on the electronic structure of NiO and CoO from DFT+DMFT calculations

本研究采用电荷自洽的 DFT+DMFT 计算，探讨了晶体结构、配体场以及变化的关联强度（包括氧 2p 轨道关联）如何影响顺磁性 NiO 和 CoO 的电子结构与谱函数。

要点

想象一个由微小、充满活力的舞者（电子）在舞台上移动的世界。在大多数材料中，这些舞者像音乐会上的观众一样自由移动。但在被称为过渡金属氧化物（如氧化镍和氧化钴）的特殊材料中，舞者们拥挤不堪且彼此高度敏感，以至于它们开始像一个紧密团结的群体那样行动。如果一个舞者移动，其他所有人都会立即做出反应。这被称为“强关联”，它使得这些材料表现出棘手而迷人的行为——有时表现为绝缘体（阻断电流），而非导体。

您提供的这篇论文深入探讨了这两种特定材料中这些舞者的确切行为：氧化镍 (NiO) 和 氧化钴 (CoO)。研究人员旨在弄清楚这些材料为何表现出如此行为，以及如何利用计算机模拟准确预测它们的行为。

以下是他们发现的分解，使用了简单的类比：

1. 问题：“糟糕的地图”

科学家有一种标准的材料映射方法，称为 DFT（密度泛函理论）。将 DFT 想象成一辆在开阔高速公路上通常表现出色的 GPS 导航仪。然而，对于这些拥挤的“舞蹈”材料，GPS 却失灵了。它预测舞者可以自由移动（金属态），但实际上，它们却被困在原地（绝缘态）。这张地图之所以错误，是因为它忽略了舞者们不断相互碰撞和反应的事实。

2. 解决方案：带有“自我修正”的更好 GPS

为了解决这个问题，作者使用了一种超级先进的方法，称为 DFT+DMFT。

• DMFT（动力学平均场理论）： 想象在模拟中加入一位“编舞家”，他观察每一个舞者，并实时计算他们如何对邻居做出反应。这解决了“金属与绝缘体”的问题。
• SIC（自相互作用修正）： 研究人员意识到，舞者不仅仅是金属原子（镍或钴）；它们周围的氧原子也是舞蹈的一部分。他们为氧舞者添加了一条特殊的“自我修正”规则，以确保模拟不会重复计算它们的能量。这种组合（DFT+DMFT+SIC）创造了一张与真实世界实验几乎完美匹配的地图。

3. 两个舞池：八面体 vs. 四面体

研究人员在两种不同的“舞池”（晶体结构）上测试了这些材料：

• 岩盐结构 (RS)： 舞者们排列成 八面体（像一个有 6 个顶点的球体）。这是这些材料稳定且自然的形状。
• 闪锌矿结构 (ZB)： 舞者们排列成 四面体（像一个有 4 个顶点的金字塔）。这是一种不稳定的人造形状，研究人员特意创造它，只是为了观察舞者们会对不同的环境做出何种反应。

发现： 就像舞者在圆形舞池和方形舞池中感觉不同一样，电子会根据舞池的形状以不同的方式分裂其能级。

• 在 八面体 (RS) 舞池中，能级以一种方式分裂。
• 在 四面体 (ZB) 舞池中，能级发生翻转。
  研究证实，八面体 舞池更稳定，因为舞者（电子）可以安顿在能量更低、更舒适的位置。四面体 舞池迫使太多舞者进入“反键”（不舒服、高能量）的位置，导致结构不稳定。这解释了为什么氧化镍自然存在于八面体形状中，而不存在于四面体形状中。

4. “卫星”与“带隙”

主要目标之一是测量 带隙——即让舞者开始移动所需的能量跳跃。

• 实验： 真实世界的实验显示带隙约为 5 到 6 电子伏特 (eV)。
• 旧模拟： 如果没有特殊的修正，计算机预测的带隙很小（约 2.5–3 eV），这是错误的。
• 新模拟： 通过对氧原子使用“自我修正” (SIC)，研究人员的计算机模型预测带隙为 5.1 eV。这与现实世界完美匹配！

他们还观察了一个“卫星峰”（数据中的高能隆起）。他们发现，虽然他们的模型改进了主带隙，但卫星峰仍然有点难以精确定位，但整体图景比以前清晰得多。

5. 镍 vs. 钴：“一步之差”

镍和钴是元素周期表上的邻居。镍比钴多一个电子。

• 氧化镍： 拥有特定数量的舞者，形成了一个非常稳定的高能“张 - 李”态（一种特殊的、紧密结合的舞者对）。
• 氧化钴： 少了一个舞者。这产生了一个“空穴”（一个空位）。研究人员发现，这个缺失的舞者使得钴材料的行为略有不同，在能带边缘产生了一个更尖锐、更强烈的峰。这就像从一个拥挤的房间里移除一个人；剩下的人会发生移动，并对空出的空间做出更强烈的反应。

总结

简而言之，这篇论文是关于构建一个 极其精确的计算机模拟，以描绘电子在氧化镍和氧化钴中如何舞蹈。

1. 他们证明，必须考虑 氧原子（而不仅仅是金属）才能得到正确的答案。
2. 他们表明，晶体形状（八面体与四面体）会改变电子分裂其能量的方式，解释了为什么某些形状是稳定的，而另一些则不是。
3. 他们实现了与真实世界实验的 近乎完美的匹配，证明他们的方法是理解这些复杂材料的可靠工具。

作者得出结论，这种方法对于希望设计更好材料的科学家来说是一个强大的工具，这些材料可用于 催化剂（用于分解水或制造燃料）和 能量存储，因为理解电子的“舞蹈”是释放其潜力的关键。

技术摘要

技术摘要：DFT+DMFT 计算中配体场与关联对 NiO 和 CoO 电子结构的影响

问题陈述
过渡金属氧化物（TMOs）如 NiO 和 CoO 表现出复杂的电子行为，这是由电子关联、配体场、晶体结构和杂化之间的相互作用所驱动的。虽然采用局域或广义梯度近似（LDA/GGA）的标准密度泛函理论（DFT）无法重现这些强关联系统的绝缘基态和正确的激发谱，但需要更先进的方法。具体而言，需要系统地理解晶体结构（岩盐型与闪锌矿型）、成分（Ni 与 Co）以及在过渡金属（TM）$3d$ 和氧 $2p$ 轨道中关联的处理方式如何影响谱函数和带隙。以往的研究往往专注于特定的相互作用参数，或者忽略了氧配体中的关联，而这对准确描述电荷转移绝缘体至关重要。

方法论
作者采用密度泛函理论与动力学平均场理论（DFT+DMFT）的电荷自洽组合，并辅以针对氧轨道的自相互作用修正（SIC）方案（DFT+sicDMFT）。

• 计算框架：DFT 部分利用混合基赝势代码（MBPP）和模守恒赝势。关联子空间定义为以 TM 离子为中心的类原子 $d$ 轨道，并在截断的 Kohn-Sham 基组中进行分解。
• 关联处理：DMFT 杂质问题使用连续时间量子蒙特卡洛（QMC）求解器 CTHYB（TRIQS 包）求解。TM $3d$ 子空间内的电子 - 电子相互作用采用旋转不变的 Slater 哈密顿量建模，包含 Hubbard 参数 $U$ 和 Hund 交换 $J$。
  • 对于 NiO：$U = 10$ eV，$J = 1$ eV。
  • 对于 CoO：$U = 9$ eV（也测试了 $U=6$ eV），$J = 1$ eV（也测试了 $J=0.9$ eV）。
• 配体关联：为了考虑与 TM 态杂化的 O($2p$) 轨道中的关联，作者对氧赝势应用了 SIC 方案。参数设定为 $\alpha=0.8$，对 O($2s$) 进行完全修正（$w_l=1.0$），对 O($2p$) 进行 80% 修正（$w_l=0.8$）。
• 研究体系：分析了顺磁态 NiO 和 CoO 在立方岩盐（RS）和闪锌矿（ZB）结构中的情况。虽然 RS 是稳定相，但 ZB 被作为模型系统研究，以隔离配体场效应（八面体配位与四面体配位）。
• 分析：实频率谱函数通过解析延拓（最大熵方法）导出。研究聚焦于总谱函数、元素投影谱函数（TM 和 O）以及局域轨道分辨谱函数，特别分析了重心（COGs）和带隙特征。

主要结果

1. SIC 和相互作用强度（$U$）的作用：

   • 对于 RS-CoO，将 $U$ 从 6 eV 增加到 9 eV，使低能 $d^7\underline{L}$ 峰（其中 $\underline{L}$ 表示配体空穴）从 $-1.0$ eV 移至 $-1.6$ eV。
   • 应用 O($2p$)-SIC 使该峰向相反方向移动至 $-0.6$ eV。$U=9$ eV 与 O($2p$)-SIC 的组合产生的谱形和峰位（O($2p$) 带为 $-3.1$ eV）与 UPS 实验数据高度吻合。
   • 至关重要的是，只有当 $U$ 足够大时，显著的带隙才会打开。包含 O($2p$)-SIC 将未占据的上 Hubbard 带（UHB）移至更高能量，导致 RS-CoO 的计算带隙为 5.1 eV。该值落在 5–6 eV 的实验范围内，相较于之前未包含 O($2p$)-SIC 的 DFT+DMFT 研究（报告的带隙低至 2.5–3.8 eV）有了显著改进。

2. 结构和配体场效应（RS 与 ZB）：

   • 计算证实了晶体场（CF）理论的预测。在 RS（八面体）中，$t_{2g}$ 态的能量低于 $e_g$。在 ZB（四面体）中，能级顺序反转（$e$ 低于 $t_2$）。
   • 通过积分谱函数进行的占据分析证实了高自旋构型：RS-NiO 为 $t_{2g}^6 e_g^2$，RS-CoO 为 $t_{2g}^5 e_g^2$。对于 ZB 结构，占据数为 $e^4 t_2^4$（NiO）和 $e^4 t_2^3$（CoO）。
   • 发现 ZB 结构由于反键态的占据而较不稳定。在 ZB 中，$t_2$ 轨道是反键的；CoO 和 NiO 中分别存在 3 个和 4 个电子占据这些状态，导致其结构相对于 RS 相变得不稳定。

3. 谱特征与屏蔽：

   • 与 NiO 相比，CoO 谱表现出更尖锐的峰和增加的低能谱权重，这归因于 $d$ 壳层中额外的空穴。
   • CoO 价带边缘附近的尖锐峰被解释为 $d^7\underline{Z}$ 态（其中 $\underline{Z}$ 表示 Zhang-Rice 空穴），类似于 Li 掺杂 NiO 中讨论的 $d^8\underline{Z}$ 行为。
   • 电荷转移能量（$\Delta$）定义为占据的 O($2p$) 态与未占据的 TM($3d$) 态之间的能量差，根据结构和元素的不同计算值为 4.3–4.9 eV，显示出与实验趋势一致的层级关系（NiO < CoO）。

意义与主张
本文主张 DFT+sicDMFT 方法为描述强关联 TMOs 的电子结构提供了一个稳健的框架。通过 SIC 显式地包含氧 $2p$ 轨道中的关联效应，该方法成功重现了顺磁态 NiO 和 CoO 的实验激发谱和带隙值，而静态平均场方法（如 DFT+U）或标准 DFT+DMFT 往往无法准确捕捉这些特征。

该研究表明：

• Hubbard $U$ 参数与配体关联之间的相互作用对于与实验的定量一致至关重要。
• 该方法可以系统地分离晶体结构（配体场）和成分对谱函数的影响。
• 计算出的谱特征，如带隙和电荷转移能量，可作为可靠的描述符，只要电子结构被准确建模，即可支持 TMOs 在催化和能量转换应用中的优化。

作者保持了谦逊的态度，指出虽然 ZB-NiO 不是稳定的体相，但将其作为模型系统研究允许对配体场效应进行系统比较。该工作基于先前关于 NiO 的发现，验证了该方法在 CoO 上的适用性，并将分析扩展到了结构变化。