氦原子散射中金属表面电子-声子耦合涉及的有效声子频率和质量

《The Journal of Physical Chemistry Letters》：The Effective Phonon Frequency and Mass Involved in Electron–Phonon Coupling at Metal Surfaces from Helium Atom Scattering

　　

在这项工作中，研究人员提出了一种经典的德拜-沃勒（Debye–Waller，DW）理论的新公式，用于描述He原子从金属表面的散射。探针原子被第一表面原子层正上方的表面自由电子密度排斥，晶格振动产生的表面电荷密度调制通过电子-声子（electron–phonon，el-ph）耦合发生。在德拜-沃勒指数中，表面投影声子谱被简化为一个负责总电子-声子耦合的有效单一声子分支。该分支的特征在于有效表面声子频率和有效表面原子质量，这些参数从高和低表面温度极限中提取。无论考虑声学还是光学声子分支，可以分别假设表面的德拜（Debye）或爱因斯坦（Einstein）模型，从而在有效参数之间形成简单的关系。根据Lapujoulade等人报道的HAS数据，爱因斯坦模型对于铜表面（110）、（113）、（001）和（111）效果良好，而对于干净和氧化的Nb(001)表面，最近由Arias、Sibener等人报道的HAS数据分析表明，表面声学声子在表面电子-声子相互作用中占主导地位。

**论文解读：氦原子散射揭示金属表面电子-声子耦合中的有效声子频率与质量**

**研究背景与问题**

氦原子散射（Helium Atom Scattering，HAS）是研究固体表面动力学的重要实验技术，其弹性散射强度随温度变化的衰减由德拜-沃勒（Debye–Waller，DW）因子描述。经典DW理论基于He原子与表面原子直接碰撞的假设，适用于绝缘体或半导体表面。然而，对于金属表面，He原子在距第一原子层上方几?处即被表面自由电子密度排斥，晶格振动引起的表面电荷密度调制仅通过电子-声子（electron–phonon，el-ph）耦合实现。因此，经典DW理论不再适用，直接提取el-ph耦合信息需要新的理论框架。此外，已有研究表明，HAS的DW指数温度依赖性与表面el-ph耦合常数（质量增强因子）直接相关，但如何将表面投影声子谱简化为可解析的物理参数（如有效频率和有效质量）仍是一个关键问题。本研究旨在提出一种新的DW理论公式，通过引入有效单一声子分支，从实验DW指数的高温和低温极限中提取有效表面声子频率和有效原子质量，从而定量刻画金属表面的el-ph耦合。

**研究内容与结论**

研究人员提出了经典DW理论的新公式，将表面投影声子谱简化为一个有效单一声子分支，该分支完全由有效频率^*ω^*和有效质量M^*表征。根据声子分支类型，可采用德拜模型（描述声学分支）或爱因斯坦模型（描述光学分支），两者通过简单关系关联。利用Lapujoulade等人（1980, 1983）报道的铜（Cu）表面（110）、（113）、（001）和（111）的HAS数据，以及Arias、Sibener等人（2022-2025）报道的干净和氧化Nb(001)表面的HAS数据，研究人员对DW指数进行双曲余切（coth）拟合，提取了有效参数。结果表明：对于铜表面，爱因斯坦模型与实验数据高度吻合，有效频率对应于表面光学声子分支（如S4、S2、S_e等），有效质量小于铜原子质量（63.55 a.u.），例如Cu(001)的M_E^*≈32.1 a.u.（约0.5M_Cu），Cu(110)的M_E^*≈47.0 a.u.（约0.74M_Cu），Cu(113)的M_E^*≈58.4 a.u.（接近M_Cu），Cu(111)的M_E^*≈26.1 a.u.（约0.41M_Cu），这些数值与表面原子反相运动（光学模式）的约化质量一致。对于Nb(001)表面，德拜模型更适用，有效频率?ω_D^*=19 meV，有效质量M_D^*=380 a.u.（约4.1倍Nb原子质量92.9 a.u.），对应表面剪切垂直（shear-vertical，SV）声学分支（瑞利波）；而爱因斯坦模型得到的有效频率（?ω_E^*=9.5 meV）则对应于表面纵向声学（surface longitudinal acoustic，SLA）分支，但DFT计算表明SV分支贡献了DW指数的主要部分（52.4%），因此德拜模型更优。氧化后的(3×1)-O/Nb(001)表面同样符合德拜模型，有效频率?ω_D^*=22 meV，有效质量M_D^*=107 a.u.（接近NbO分子质量），表明氧的吸附稳定了表面并改变了声子谱。该研究的重要意义在于：提供了一种直接从HAS实验数据提取表面el-ph耦合参数（有效频率和质量）的系统方法，并揭示了铜表面与铌表面el-ph耦合机制的本质差异——前者由光学声子主导，后者由声学声子主导，这为理解金属表面超导性、表面粗糙化转变等提供了新的实验和理论工具。论文发表在《The Journal of Physical Chemistry Letters》。

**主要关键技术方法**

研究人员主要采用了以下关键技术方法：（1）基于HAS实验的弹性镜面散射强度随温度变化的数据，数据来源包括Lapujoulade等人（1980, 1983）对Cu(110)、(113)、(001)、(111)表面的测量，以及Arias、Sibener等人（2022-2025）对干净和氧化Nb(001)表面的测量。实验数据通过双曲余切（coth）函数拟合，得到DW指数2W(T)在高温极限（线性渐近线）和零温极限（2W(0)）的值，从而利用公式（9a）和（9b）提取有效爱因斯坦频率?ω_E^*和有效质量M_E^*，或通过德拜模型公式（10）得到等效德拜参数。（2）为校正He原子在表面吸引势中动能增加，采用了Beeby修正（公式11），将入射动量垂直分量替换为有效值。（3）理论计算方面，使用了嵌入原子方法（Embedded Atom Method，EAM）、多体经验势（Many-Body Empirical Potential，MBEP）、密度泛函微扰理论（Density Functional Perturbation Theory，DFPT）以及力常数模型（Force-Constant Model，FCI）等计算表面声子色散曲线，与实验提取的有效参数进行对比。

**研究结果**

**Cu(110)表面**：通过HAS数据（入射波矢k_i=6.35 ?^–1，极角θ_i=20°）的coth拟合，得到有效爱因斯坦频率?ω_E^*=22.1 meV，有效质量M_E^*=47.0 a.u.。该频率对应表面光学声子分支S4，该分支在整个表面布里渊区（BZ）内具有爱因斯坦类平坦色散；德拜模型给出的有效声学频率?ω_D^*≈44.2 meV，对应异常纵向共振S3分支，但有效质量M_D^*≈35.25 a.u.，表明光学模式主导el-ph耦合。

**Cu(113)表面**：HAS数据（相同入射条件）得到?ω_E^*=19.1 meV，M_E^*=58.4 a.u.。有效频率对应BZ中心的光学模式S4和BZ边界台阶边缘模式S_e（准一维声子分支），后者沿台阶原子行局域。有效质量接近铜原子质量，与边缘模式沿单原子链传播的特性一致。

**Cu(001)表面**：两组不同入射极角（55.5°和31.8°）的HAS数据均得到一致的?ω_E^*=23.8 meV，M_E^*≈32.1 a.u.。有效频率对应BZ边界间隙表面模式S2，该模式在第一和第二原子层中具有离子纵向位移，产生强el-ph耦合。有效质量约为铜原子质量的一半，与fcc(001)表面双原子层反相运动（光学模式）的约化质量一致。

**Cu(111)表面**：由于低温数据不足，无法直接coth拟合，研究人员利用DFPT计算（Chis等人，2008）识别出爱因斯坦类S2表面光学共振（第二层剪切垂直极化）和S3异常纵向共振（BZ边界）作为主要el-ph贡献。通过将?ω_E^*与S2分支关联（25.0 meV），得到M_E^*=26.1 a.u.，约为铜原子质量的0.41倍。这一数值与表面三原子层单元（trimer）中反相振动的理论约化质量（0.382M）吻合，验证了模型的有效性。

**清洁Nb(001)表面**：基于Méndez等人（2025）的ab initio HAS计算，对DW指数进行拟合得到德拜模型参数?ω_D^*=19 meV，M_D^*=380 a.u.；爱因斯坦模型参数?ω_E^*=9.5 meV，M_E^*=506 a.u.。德拜频率对应BZ边界的声学SV（瑞利波）分支的线性外推，而爱因斯坦频率对应平坦的SLA分支。DFT计算表明SV分支贡献了DW指数52.4%，SLA仅贡献7.4%，因此德拜模型更合适。有效质量远大于Nb原子质量，反映了声学模式位移场指数衰减且涉及多层原子的特征。

**氧化(3×1)-O/Nb(001)表面**：类似分析得到德拜模型参数?ω_D^*=22 meV，M_D^*=107 a.u.，对应SLA分支的线性拟合，且有效质量接近NbO分子质量，表明氧的引入提高了SLA分支频率，稳定了表面结构。

**总结与讨论**

研究人员在讨论中指出，对于金属表面HAS，He原子不直接碰撞表面原子，而是通过表面自由电子密度相互作用，因此DW指数的主要贡献来自具有最大el-ph耦合的声子模式。铜表面普遍存在光学声子模式，其反相运动导致表面电荷密度振荡最大，故爱因斯坦模型适用；而铌表面缺乏局域光学模式，声学声子（尤其是瑞利波）主导el-ph耦合，故德拜模型更优。经典德拜模型在铜表面给出过大的有效频率（远高于体铜德拜频率），而爱因斯坦模型与实验及理论计算高度一致。对于Nb(001)表面，SLA分支的异常软化（低于体LA带边缘）是el-ph耦合增强的直接表现，该软化在氧化后因氧的稳定作用而减弱。研究结论部分翻译如下：总之，已证明在金属表面的HAS中，He原子不直接与表面原子碰撞，而是通过表面自由电子密度相互作用，DW指数的主要贡献来自具有最大el-ph耦合的表面声子。通常这些是表面光学声子，其表面原胞内原子的反相位移导致最大的表面电子密度振荡。具有单个有效光学频率和小于金属原子质量的有效质量的爱因斯坦模型，为解释金属表面HAS的DW因子温度依赖性提供了适当的近似。这已在四种不同的铜表面上得到证实，从实验DW指数温度依赖性推导出的有效爱因斯坦频率与最相关的表面光学分支非常吻合。另一方面，自HAS光谱学早期以来文献中采用的基于德拜模型的经典DW因子近似，对于铜表面给出了过大的德拜有效频率，远高于体铜的已知值。尽管如此，有趣的是，对应于ω_D^*的声学模式在长波极限下的斜率正是普遍存在的异常纵向声学共振的斜率，这种共振作为表面增强el-ph耦合的效应，几乎在所有金属表面都能通过HAS发现。尽管有这种对应关系，将金属表面DW指数温度依赖性与具有大el-ph耦合的爱因斯坦类光学分支联系起来要合理得多。这是否也为Lapujoulade等人在Cu(110)和邻位面上观察到的500 K以上He原子散射表面粗糙化转变（其中λ_HAS显著增加）提供了一种机制，值得进一步研究。干净和氧化铌(001)表面的情况提供了与铜表面互补的启发性场景。在铌（常压下最好的元素超导体）中，el-ph耦合远大于铜，但在上述两个表面上，由于缺乏（或弱）表面局域光学声子分支，其耦合显著降低。对于这两个表面，el-ph耦合主要由表面声学分支支持，因此更好地用经典德拜近似建模。此外，应用于清洁Nb(001)表面的爱因斯坦模型，其中SLA分支在布里渊区大部分区域异常软化且平坦，表明该分支在铌表面el-ph相互作用中具有相关作用。