尖端增强拉曼光谱（TERS）是一种强大的技术，可用于对表面束缚体系的振动模式和化学特性进行成像表征。然而，过去大多数的理论模拟往往将系统视为孤立模型，忽略了诸如金属表面等基底支撑的影响。本文的研究表明，这种忽略会导致模拟结果与实验测量数据产生偏差。为此，作者开发了一种新的有限场周期公式化方法，用于第一性原理模拟TERS光谱，能够处理扩展体系，从而解决了上述问题。
研究结果显示，使用团簇模型计算的缺陷二硫化钼（MoS2）单层TERS图像，与考虑基底周期性后计算得到的图像存在显著的定性差异。对于一个可与实验直接对比的体系——镁(II)卟啉（MgP）吸附在银（100）表面（MgP/Ag(100)），这些模拟对于解释TERS强度图案的空间变化至关重要，并帮助研究人员揭示了TERS光谱学的基本原理。研究阐明并解释了为什么表面相互作用对垂直平面（out-of-plane）振动模式图像的影响远大于面内（in-plane）模式，从而为未来解释更复杂体系的TERS图像提供了重要工具。
方法学进展
作者提出的新方法是对先前工作的延伸，能够从第一性原理出发，在密度泛函理论（DFT）模拟中完全考虑等离激元近场的原子描述以及扩展金属表面的周期性。该方法的技术细节在文中及其支持信息中详述，并被应用于研究缺陷和吸附物的TERS信号。
缺陷二硫化钼单层的TERS成像
作为新方法的首次应用，研究人员模拟了单层MoS2中与缺陷相关的振动模式的TERS图像，特别是最常见的硫单空位。在完美晶格中，对应于A1′振动的信号在397 cm-1处占主导地位，其TERS图像在系统内表现出非常均匀的强度。然而，在含有4%硫空位浓度（5×5超胞）的缺陷体系中，由于对称性降低，在A1′峰（396.3 cm-1）的右侧出现了一个低强度的肩峰（400.4 cm-1，即D带）。两种缺陷诱导振动的TERS图像都显示出以空位为中心的特征图案：A1′模式在缺陷位点中心呈现较低强度区域，周围环绕着一个明显具有C3对称性的高强度环；而D带则呈现出以空位为中心的三叶草状对称性。这些结果强调了在周期性系统（而非有限尺寸的薄片模型）中模拟的重要性，以获得正确的图像对称性。
镁(II)卟啉/Ag(100)体系的TERS成像分析
为了验证方法的可靠性并揭示TERS图像图案的基本原理，研究者系统地分析了MgP/Ag(100)体系。该分子在吸附后基本保持其平面D4h对称性。研究者选取了三个代表性振动模式进行分析：193.4 cm-1的A2u模式（主要由Mg原子面外振动主导）、1359.9 cm-1的B1g模式（吡咯环的非对称呼吸模式）以及3180.3 cm-1的A2g模式（氢原子非对称伸缩振动）。
表面效应与屏蔽作用
模拟揭示，先前常用的气相近似（孤立分子模型）在某些情况下（如A2g模式）可能给出与实验接近的图像，但这并不通用。例如，对于A2u模式，气相模拟在中心Mg原子处呈现一个显著的强度最小值，而实验图像中此处却是最高强度区。当模型显式包含原子级Ag(100)周期性表面后，所有模式的模拟结果与实验图像的一致性均显著提升。A2g模式保持了良好的一致性，B1g模式的峰位与实验吻合，而A2u模式则成功再现了中心的高强度区域。
进一步的深入分析表明，这种图像形状的显著变化源于金属基底对分子振动过程中诱导偶极的电荷屏蔽（screening）效应。在气相中，A2u模式的TERS振幅在整个分子范围内几乎全部为正，导致图像中Mg原子处为极小值。而在表面吸附的情况下，振幅变为负值，从而使图像中的中心区域转变为极大值。这一符号反转主要是由Ag(100)基底对振动过程中电子密度响应的屏蔽引起的，并且是全局性的，与近场分布细节无关。研究指出，只有那些在散射方向（z方向）具有非零拉曼强度的模式（如面外A2u模式）才会表现出这种显著的符号变化效应。
TERS强度的原子贡献分解
为了更深刻地理解TERS强度增强图案的来源，研究者对TERS强度进行了原子级分解。TERS振幅（Azz）可表达为对所有原子贡献之和。因此，TERS强度（Izz）可以分解为自项（self terms）和交叉项（cross terms）。自项恒为非负，而交叉项则可为负。
对于MgP的A2u模式分析发现，在气相中，强烈的交叉项几乎完全抵消了自项，导致尽管Mg原子在该模式中运动幅度最大，但其中心处的总强度却很小。而在表面上，交叉项的幅度显著减小，导致两者抵消不完全，从而在原点处观察到非零的强度峰。这表明，连接不同原子的交叉项可以导致TERS图像出现超越振动态几何形状的非平凡图案，挑战了“TERS强度较大处原子运动幅度较大”的常见假设。
结论与展望
综上所述，本研究通过严格的基准测试和实际体系模拟证明，将表面近似为团簇或完全忽略（气相计算）会导致包含人为假象的图像，对周期性的正确处理是准确建模的关键。对MgP/Ag(100)体系的模拟与实验直接对比，加深了人们对TERS图像形状起源的理解，揭示了其对分子-表面距离的依赖性，以及基底屏蔽效应的关键影响。
文中所发展的方法适用于与基底有强相互作用的多种体系，无论基底电子性质如何。未来的工作可以朝多个方向扩展，例如考虑更大的针尖模型、处理共振TERS信号、包括再辐射增强机制以及超越基态DFT描述的现象等。这项工作所实现的理解，得益于文中提出的新型第一性原理计算框架。该方法已集成在FHI-aims软件中，并结合自动化工作流程，可供学界广泛使用，为研究任意复杂度的体系开辟了新的道路，能够有力地支撑和解释当前TERS实验，并揭示这些测量背后的新物理机制。
本文主要方法基于有限场微扰理论，将含时问题转化为受扰动的静态哈密顿量问题，并利用密度泛函微扰理论（DFPT）或有限差分法计算依赖于针尖位置的局域极化率张量zz分量。该方法依赖于几个合理假设：针尖与分子之间无电子密度重叠、外加电场强度处于线性响应区、非共振拉曼散射条件。通过自洽求解包含针尖近场耦合的哈密顿量，计算散射子系统的诱导偶极矩，进而通过有限差分获得极化率，最终在谐振近似下生成TERS图像。所有计算基于PBE泛函，并使用Tkatchenko–Scheffler色散修正。
