[论文解读] kMap.py: A Python program for simulation and data analysis in photoemission tomography

Dominik Brandstetter, Xiaosheng Yang+3 more|arXiv (Cornell University)|Sep 28, 2020

Electron and X-Ray Spectroscopy Techniques参考文献 51被引用 22

一句话总结

kMap.py 是一款基于 Python 的开源软件，采用 PyQt 图形用户界面（GUI），利用密度泛函理论（DFT）计算的轨道密度，通过平面波近似和快速傅里叶变换（FFT）模拟并分析分子轨道的光电子动量图（k-图）。该工具可实现模拟与实验角度分辨光电子能谱（ARPES）数据的直接一对一比较，支持自动化优化，高精度提取分子取向、轨道权重和几何参数。

ABSTRACT

For organic molecules adsorbed as well-oriented ultra-thin films on metallic surfaces, angle-resolved photoemission spectroscopy has evolved into a technique called photoemission tomography (PT). By approximating the final state of the photoemitted electron as a free electron, PT uses the angular dependence of the photocurrent, a so-called momentum map or k-map, and interprets it as the Fourier transform of the initial state's molecular orbital, thereby gains insights into the geometric and electronic structure of organic/metal interfaces. In this contribution, we present kMap.py which is a Python program that enables the user, via a PyQt-based graphical user interface, to simulate photoemission momentum maps of molecular orbitals and to perform a one-to-one comparison between simulation and experiment. Based on the plane wave approximation for the final state, simulated momentum maps are computed numerically from a fast Fourier transform of real space molecular orbital distributions, which are used as program input and taken from density functional calculations. The program allows the user to vary a number of simulation parameters such as the final state kinetic energy, the molecular orientation or the polarization state of the incident light field. Moreover, also experimental photoemission data can be loaded into the program enabling a direct visual comparison as well as an automatic optimization procedure to determine structural parameters of the molecules or weights of molecular orbitals contributions. With an increasing number of experimental groups employing photoemission tomography to study adsorbate layers, we expect kMap.py to serve as an ideal analysis software to further extend the applicability of PT.

研究动机与目标

• 开发一款用户友好、开源的软件工具，用于模拟和分析有机/金属界面中的光电子动量图。
• 实现分子轨道模拟与实验 ARPES 数据的直接、可视化比较。
• 提供自动化优化程序，从实验动量图中提取分子取向、轨道权重和结构参数。
• 通过提供灵活、易访问的模拟平台，支持光电子断层成像在研究分子吸附层中的日益广泛应用。

提出的方法

• 使用最终态的平面波近似模拟动量图。
• 通过密度泛函理论（DFT）计算得到的实空间分子轨道分布的快速傅里叶变换（FFT）数值计算模拟动量图。
• 提供基于 PyQt 的图形用户界面（GUI），支持交互式模拟和数据加载。
• 允许用户调整关键参数，如最终态动能、分子取向和入射光偏振态。
• 支持加载实验 ARPES 数据，实现与模拟结果的直接可视化比较，并采用最小二乘法进行自动化拟合。
• 基于最小化模拟与实验动量图之间 χ² 差异的优化程序，确定轨道权重和结构参数。

实验结果

研究问题

• RQ1 如何高效地模拟分子轨道的动量图，并与实验 ARPES 数据进行比较？
• RQ2 在实验动量图中，通过自动化拟合可多大程度上提取分子取向和轨道权重？
• RQ3 kMap.py 是否能准确从实验 k-图中重构分子倾斜角等几何参数？
• RQ4 平面波近似在模拟复杂有机分子在金属表面的动量图时表现如何？
• RQ5 kMap.py 在将重叠的谱带特征解卷积为单个轨道贡献方面具有多大潜力？

主要发现

• kMap.py 成功再现了在 Ag(111) 表面的 pentacene 实验动量图，最佳拟合分子倾斜角为 ϑ ≈ 20°，与先前研究结果一致。
• 该程序准确地将 PTCDA/Ag(110) 的 M3 发射特征解卷积为四个紧密相邻的分子轨道（C、D、E、F）的贡献，得到轨道投影的态密度。
• 通过拟合获得的轨道权重显示出明显的动能依赖性，证实了解卷积过程的有效性。
• 该软件可在极少用户输入下精确确定分子取向和轨道贡献，表现出高度的可重复性和准确性。
• 基于 FFT 的模拟实现支持快速、实时的可视化与优化，显著提升了分析效率。
• 未来扩展方向，如实空间轨道重构的相位恢复以及超越平面波最终态模型的构建，被识别为可行且前景广阔。