Accurate and efficient simulation of photoemission spectroscopy via Kohn-Sham scattering states

本文引入了一种高效的第一性原理框架，该框架将光电子态计算为 Kohn-Sham 散射解，从而实现对角分辨光电子能谱（ARPES）准确、透明且广泛兼容的模拟，并通过与石墨烯和 WSe$_2$ 实验数据的极佳一致性证明了这一点。

要点

想象一下，你正试图拍摄一张完美的城市夜景照片。你不仅想看到建筑物的轮廓，还想精确地观察光线是如何在窗户上反射的，阴影是如何投射的，以及城市的布局如何改变从你特定相机角度看过去的光影效果。

在物理学世界中，**角分辨光电子能谱（ARPES）**就像是那台相机。科学家用高能光子（光）照射材料，将电子从材料中撞击出来。通过测量这些飞行的电子的速度和方向，他们可以绘制出材料内部的“城市规划图”——即其电子结构。

然而，这里有一个陷阱。你得到的照片不仅仅是城市的直接图像，它还是城市布局以及光线在到达相机之前如何反射、撞击墙壁并产生干涉的复杂混合体。长期以来，在计算机上模拟这一过程就像是在解一个巨大的、缠绕在一起的绳结。现有的方法要么过于僵化（仅适用于特定类型的材料），要么速度太慢，无法满足现代复杂实验的实际需求。

全新的“相机镜头”方法
本文介绍了一种高效模拟这些照片的新方法。作者 Gian Parusa 及其团队开发了一种将逃逸电子视为波浪拍打海岸的方法。

他们没有使用只能针对特定材料的复杂专用软件，而是构建了一个能够与大多数科学家已经在使用的标准“蓝图”（计算机代码）协同工作的工具。他们的方法求解了一个特定的数学问题（Kohn-Sham 方程），并采用了特殊的边界条件，告诉计算机：“想象这些电子正从材料中向着真空空间奔跑而去。”

为什么这种方法更好？
可以这样理解：

• 旧方法： 就像是试图通过手工制作每一块砖头来盖房子。这行得通，但速度很慢，而且你无法轻易更改设计。
• 这种新方法： 就像是使用一套高质量的预制墙体系统，它可以适配任何标准的房屋设计。它快速、灵活，让你在还没盖好房子之前，就能看清光线是如何照射在墙上的。

机器中的“幽灵”：伪势（Pseudopotentials）
在这些模拟中，最大的障碍之一是处理沉重的原子核（原子核及内层电子）。为了节省计算能力，科学家经常使用“伪势”——它们就像是代表重原子的简化面具，无需计算每一个微小的细节。

团队测试了这些“面具”是否足以准确预测高速电子的散射情况。他们发现：

1. 简单的面具对于许多材料都表现良好，前提是该面具具有高质量。
2. 然而，对于重原子（如 WSe2 中的钨），面具需要包含“深层秘密”（半芯态/semicore states）。 如果遗漏了这些，模拟会对“阴影”产生错误的判断，导致照片失真。这就像戴着一副遮住了眼睛却忘了遮住耳朵的面具；你能看见，但你会错过那些会改变你对世界反应的关键声音信号。

证据：石墨烯与 WSe2
为了证明他们的方法有效，他们模拟了两种材料：

1. 石墨烯（单层碳）： 他们预测了光模式（称为圆二色性）的表现。他们的模拟与现实世界的实验完美契合，甚至预测到了其他方法错过的微妙“节点线”（信号消失的地方）。
2. WSe2（体相晶体）： 他们展示了在面具中包含那些“深层秘密”（半芯态）对于获得正确模式是至关重要的。如果没有它们，模拟结果看起来就像是真实实验的一个模糊且错误的变体。

底线
这篇论文不仅提供了一种更快的数学计算方式，它还提供了一个更清晰的窗口，让我们观察光与物质是如何相互作用的。通过计算逃逸电子精确的“飞行路径”，科学家现在可以：

• 理解为什么实验中会出现某些特定的图案。
• 区分材料的真实本质与测量过程造成的“光学幻觉”。
• 使用标准且广泛可用的计算机工具，研究复杂的材料，甚至是处于运动状态的材料（例如被激光泵浦的材料）。

简而言之，他们为研究人员提供了一个更锐利、更灵活的镜头，用以观察固体内部隐形的电子世界。

技术摘要

技术摘要：通过 Kohn-Sham 散射态实现精确且高效的光电子能谱模拟

问题陈述
角分辨光电子能谱（ARPES）是绘制固体电子结构结构的主要工具，然而在本质材料属性（初态）与实验观测光谱之间存在显著差距。观测到的强度不仅取决于初始 Bloch 态，还关键取决于最终光电子态，而后者受实验几何结构、光极化、原子内及原子间散射以及晶体势场内多重散射的影响。虽然“一步法”（one-step）光电子理论提供了一个严谨的框架，但现有的实现方式——特别是成熟的 Korringa-Kohn–Rostoker (KKR) 多重散射方法——通常与专门的电子结构形式化方法绑定。这限制了它们与广泛使用的平面波或实空间密度泛函理论（DFT）代码的集成，阻碍了先进交换相关泛函和现代工作流的应用。此外，对于描述高能光电子态而言，对计算至关重要的伪势（PP）在准确性方面仍是一个开放性问题。

方法论
作者引入了一种第一性原理框架，通过直接求解带有散射边界条件的 Kohn-Sham (KS) 方程来计算最终光电子态。

• 形式化方法： 该方法假设系统为二维周期性系统，并采用波函数的面内 Laue 表示。通过施加渐近边界条件（即在 $z \to +\infty$ 时为出射波，在 $z \to -\infty$ 时为入射/驻波），将 KS 方程从特征值问题（用于束缚态）转化为散射态的线性方程。
• 与 Lippmann-Schwinger 的等价性： 作者证明，通过这些边界条件求解 KS 方程在形式上等价于 Lippmann-Schwinger (LS) 公式。然而，直接的 KS 方法允许使用源自迭代对角化的高效预处理策略，从而比直接求解 LS 积分方程具有显著更快的收敛速度。
• 实现： 该方法在标准的平面波和实空间 DFT 代码中实现。它利用离散的 $z$ 网格，并采用高阶离散化（最高达 $\mathcal{O}(h^5)$）和迭代线性求解器。
• 伪势分析： 本研究的一个核心组成部分是对伪势的系统评估。作者将全电子（AE）计算与广义正规守恒伪势（PP）进行了比较，特别研究了包含半芯态（semicore states）和非局部势对于准确重现高能散射属性的必要性。

主要贡献

1. 高效框架： 本文提出了一种与标准 DFT 工作流（平面波和实空间）完全兼容的方法，使得使用先进泛函和大规模模拟成为可能，而这些在以往专门的一步法代码中是难以实现的。
2. 显式波函数： 不同于仅产生强度的算法，该方法提供了显式的光电子波函数。这使得对矩阵元效应、多重散射以及将本质初态属性与外在末态效应进行解耦的研究变得透明且可行。
3. 伪势验证： 研究表明，高质量的伪势（特别是具有优化非局部部分并包含半芯态的伪势）可以准确描述真空紫外（VUV）到极紫外（XUV）波段的光电子态，使该方法具备常规使用的实用性。

结果
该框架通过对两个系统的圆偏振角分辨光电子能谱（CD-ARPES）模拟得到了验证：

• 单层石墨烯： 将能量积分后的圆偏振二色性角分布（CDAD）理论结果与实验数据（文献 [61]）进行了对比。该方法成功重现了节线（nodal lines）、整体 CDAD 纹理以及 75 eV 处的 CDAD 符号反转。全电子法与伪势法的比较显示，只要伪势包含非局部项，两者在 100 eV 以下表现出极佳的一致性。研究发现非局部部分对于重现散射属性至关重要，尽管其影响随光子能量和角动量通道的变化而不同（例如，在 52 eV 处对碳的 $d$ 类偏波贡献的影响可以忽略不计）。
• 体相 2H-WSe$_2$： 将在 CELIA 设施获取的新实验数据与理论模拟进行了对比。结果表明，在伪势中包含半芯态（W 的 $5s, 5p$ 和 Se 的 $3d$）是必不可少的。使用含半芯态伪势的模拟在所有 $K$ 点上均实现了与实验谷集成 CDAD 的定量一致；而未使用半芯态的模拟则表现出显著偏差，特别是在 $K$ 点 1 和 2 处。这证实了准确的原子位点周围散射相位对于捕捉光电子发射的干涉本质是必要的。

意义与主张
作者声称，这项工作为定量 ARPES 建模建立了一条稳健且易于获取的路径。通过弥合广泛使用的从头算平台与光电子理论之间的鸿沟，该方法为轨道纹理、多体效应和非平衡态现象的系统研究打开了大门。

• 多功能性： 显式计算矩阵元允许对任意实验几何结构和光极化进行后处理。
• 可扩展性： 该方法可以与格林函数技术（如 $GW$）结合，以处理超越标准 DFT 的关联效应。它还被定位为模拟时间分辨 ARPES（trARPES）在泵浦-探测场景下的基础，为目前局限于标准 DFT 的 KKR 基 trARPES 实现提供了一种潜在替代方案。
• 局限性： 作者谦逊地指出，处理具有极大胞元的表面可能具有挑战性，且真正的半无限系统可能仍需要像 KKR 这样的格林函数形式。此外，在极高光子能量下，核态的贡献可能需要通过伪波函数重建全电子波函数。

总之，本文提供了一种计算高效、在形式上等价于 Lippmann-Schwinger 公式的方法，该方法能与现代 DFT 代码无缝集成，并通过高精度的实验 CD-ARPES 数据得到了验证。