Predicting core-level X-ray photoemission spectra of oxide surfaces from first principles -- a case study for SnO$_2$

本文提出了一种基于第一性原理的 Z+1 方法，用于预测各种 SnO$_2$(110) 表面终止面及缺陷态的内壳层 X 射线光电子能谱，证明了含有吸附物的还原表面之计算谱图与实验测量结果高度吻合，并成功区分了不同的表面化学环境。

要点

想象一下，你正试图仅通过观察一个物体的影子来推测它的成分。这本质上就是科学家在研究像二氧化锡（$SnO_2$）——一种用于气体传感器和透明电子器件的物质——时所使用的技术：X射线光电子能谱（XPS）。

在XPS中，科学家向材料发射X射线，以将电子从原子中撞击出来。通过测量撞击出这些电子所需的能量（“结合能”），他们可以了解表面上有哪些种类的原子以及它们的行为方式。然而，这里有一个大问题：现实世界的表面是杂乱无章的。它们有缺失的原子、多余的原子以及粘附在上面的分子。这就像是在一个拥挤且雾气弥漫的房间里，仅凭听到的声音来识别一个特定的人；这些声音（或者说在这种情况下是光谱峰）会发生重叠并变得令人困惑。

问题：嘈杂的房间

多年来，科学家们一直在争论当二氧化锡暴露在氧气中时，其表面究竟发生了什么。有人认为氧分子粘附在表面并抓取了电子；也有人认为表面存在“孔洞”（空位），即原子缺失的地方，而氧填充了这些孔洞。

麻烦在于，实验数据（即“影子”）对于不同的情景看起来非常相似。如果没有一张清晰的地图，很难知道哪种理论才是正确的。

解决方案：第一性原理地图

本文的作者利用计算机模拟创建了一张“地图”，用以预测不同表面条件下XPS“影子”应该呈现出的精确样子。他们并非仅仅靠猜测；他们是从底层开始，利用量子物理学构建了一个数字模型。

为了使数学计算更简单、更稳定，他们使用了一个被称为 $Z+1$ 方法 的巧妙技巧。

• 类比： 想象一下，如果你想观察从一个氧原子中移除一个特定电子会发生什么。与其尝试去计算留下“空穴”后那混乱的物理过程，他们干脆假装这个氧原子已经被替换成了一个氟原子（氟原子多出一个质子）。
• 为什么有效： 这就像是用一个尺寸略有不同但能完美契合的齿轮，去替换掉钟表里一个损坏的齿轮，从而让时钟能够继续走动，以便你测量时间。这使他们能够精确计算能级，而不会导致计算机崩溃。

侦探工作：测试不同的表面

团队构建了五种不同状态下的二氧化锡表面数字模型，并预测了它们的XPS“影子”会是什么样子：

1. 完美的表面（化学计量比表面）： 一个洁净、平衡的表面。
   • 预测： 该表面会在低能量处显示出一个由“桥联”氧原子（像桥一样坐落在上方的原子）引起的奇怪额外凸起。
2. “完全还原”的表面： 一个许多氧原子缺失（产生空位）的表面。
   • 预测： 该表面会产生一个非常平滑、对称的峰。
3. “愈合”的表面： 带有吸附氧分子（$O_2$）或水分子（$H_2O$）的还原表面。
   • 预测： 这些表面会在高能级处显示出一个新的“肩峰”或凸起。

结论：匹配线索

研究人员将他们的数字预测与其他人（Kucharski及其同事）进行的真实实验进行了对比。

• 暴露于氧气之前： 真实的实验数据显示出一个平滑、对称的峰。这与**“完全还原”**模型完美吻合。这意味着科学家们观察到的表面实际上充满了缺失的氧原子（空位），而不是一个完美的表面。
• 暴露于氧气之后： 当真实表面暴露在氧气中时，光谱的高能端出现了一个新的凸起。
  • 计算机模型显示，无论是**吸附的氧分子（$O_2$）还是羟基（OH）**都会产生这个高能端的凸起。
  • 作者得出结论，表面的“愈合”不仅仅是氧填充了空位；更有可能是氧分子粘附在表面或形成了OH基团，这产生了那个特定的高能信号。

核心要点

该论文声称，通过使用这种特定的计算机方法（$Z+1$），他们可以准确预测复杂、杂乱的表面在XPS光谱中应该呈现的样子。

他们发现，那个“杂乱”的表面（充满空位）在数据中看起来反而最干净，而“洁净”的表面看起来却很杂乱。此外，引入氧气后看到的额外信号，很可能是由于氧分子或OH基团粘附在表面，而非仅仅是简单的空位被填充。

简而言之，他们构建了一个可靠的翻译器，将X射线数据中那些令人困惑的“噪声”转化为关于表面原子层级实际发生情况的清晰故事。这有助于科学家们停止猜测，转而确切地了解这些材料上的化学环境。

技术摘要

问题陈述
X射线光电子能谱（XPS）是分析氧化物表面化学性质的关键技术，然而其解释过程仍然具有挑战性。真实的表面通常具有复杂的终止面、重构、吸附物和缺陷，每种因素都可能留下潜在重叠的光谱指纹。特别对于二氧化锡（SnO2）——一种广泛用于传感器和催化领域的n型半导体——控制表面化学（例如在氧气暴露下的电阻变化）的微观机制仍存在争议。传统的解释依赖于离子吸附理论或氧空位浓度的变化，但由于缺乏针对无数种可能吸附物构型的准确参考光谱，将实验O 1s能谱中的特定峰归属于特定的化学环境变得十分困难。

方法论
作者提出了一种从第一性原理出发预测氧化物表面XPS光谱的方法，特别侧重于SnO2 (110)。该方法的核心在于利用 Z+1 近似法 计算核心电子结合能（CEBE）位移。在该方案中，具有核心空位的原子核被替换为一个增加了一个质子的原子核（例如，将一个O原子替换为F原子），以在不违反泡利不相容原理的情况下模拟核心空位。

关键方法组成部分包括：

• 建模： 使用包含最多7层原子的金红石型SnO2 (110) 周期性平板模型，并设置至少50 Å的真空层。
• 研究体系： 包括化学计量比表面、完全还原表面（所有桥联氧均被移除）、以及具有各种氧空位（面内、亚桥联、亚面内）的表面，以及带有吸附O2、OH和H2O的完全还原表面。
• 计算细节： 计算使用全电子FHI-aims代码，利用PBE泛函进行结构弛豫和Z+1计算，并使用B3LYP计算态密度。通过缩放的零阶相对论近似（ZORA）捕捉相对论效应。
• 光谱构建： 基于计算出的CEBE位移，通过高斯函数的加权和生成模拟的O 1s光谱，并纳入光电子逃逸概率（平均自由程 $\lambda = 2$ Å）和展宽（$\sigma = 0.9$ eV）。
• 验证： 通过将Z+1方法与完整的 $\Delta$-自洽场（$\Delta$-SCF）方法以及基态Kohn-Sham本征值进行对比，以评估其准确性以及末态效应的作用。

主要结果

1. 方法验证： 通过Z+1方法获得的CEBE位移与计算成本更高的 $\Delta$SCF 方法表现出极佳的一致性，差异仅为0.02–0.05 eV。
2. 化学计量比表面与还原表面：
   • 化学计量比表面表现出一个额外的低结合能峰，这归因于桥联氧原子；由于其配位数较低且具有较高的负Mulliken电荷，这些原子具有较大的负CEBE位移（约 -1.30 eV）。
   • 完全还原表面产生了一个高度对称的O 1s光谱，与近期测得的还原表面（在接触氧气前）的实验测量结果相匹配。该表面具有金属特性，末态屏蔽效应显著降低了相对于基态计算的CEBE位移。
3. 还原表面的吸附物：
   • O2吸附： 在完全还原表面上吸附的O2分子会产生相对于主峰的正CEBE位移（更高的结合能）。最外层O原子的位移为0.44 eV，而桥联位置的O原子位移为1.49 eV。研究发现O2分子从表面接受了两个电子。
   • OH吸附： 吸附的OH基团同样产生高结合能特征（位移约为 +1.44 eV）。
   • H2O吸附： 吸附的水分子导致显著更高的结合能位移（约为 +4.1 eV）。
4. 与实验对比： 计算出的暴露于氧气（含有O2或OH吸附物）的还原表面光谱显示出一个高能肩峰，这与Kucharski等人观察到的实验结果一致。研究表明，实验中观察到的高能特征很可能是由吸附的氧分子（O2）或OH基团引起的，而非氧空位。

意义与主张
论文声称，所提出的基于Z+1的方法是一种强大的通用工具，可用于支持实验XPS光谱的分析。通过将第一性原理预测与测量值直接对比，作者为氧化物表面氧环境的化学本质提供了微观见解。

具体而言，这项工作挑战了将某些XPS峰传统地归因为“氧空位”的做法。作者认为，通常被标记为源自氧空位的峰，实际上可能起源于“额外的氧”，即通过与空位相互作用而形成的OH或桥联氧物种。研究结论指出，该方法可以简单地应用于其他表面，以解决复杂的氧化物表面化学中的歧义，并为缺乏实验参考资料的谱学数据解释提供可靠的参考。