合金团簇上振动红移、电子结构与水离解之间的相关性：一项探索性计算研究

《Molecules》：Correlations Between Vibrational Red Shifts, Electronic Structure, and Water Dissociation on Alloy Clusters: An Exploratory Computational Study

　　

理解水分子的活化机理对于通过水分解实现催化水离解具有根本重要性。在本研究中，研究人员采用第一性原理分子动力学（First-Principles Molecular Dynamics, FPMD）模拟研究了H₂O在88个合金团簇上的催化离解。结果表明，吸附态H₂O的耦合ν₁和ν₃伸缩模式中心的红移越大，离解温度越低。研究人员观察到一种经验相关性，即当吸附能大小与HOMO–LUMO能隙（最高占据分子轨道–最低未占分子轨道能隙，HOMO–LUMO gap）相当时，离解更容易发生；这被提出为一种有待激发态方法进一步验证的假说。此外，前驱体与中间态的前沿分子轨道比较表明，与离解的H原子重叠增加的前沿轨道数目与离解温度呈负相关。这些发现为孤立气相团簇上水离解的活化机理提供了原子尺度的见解，并建立了评估此类模型体系本征反应性的光谱描述符。研究人员提醒，将这些发现推广到实际催化水分解时，需要考虑催化剂载体、溶剂和实际反应条件等其他因素。

水离解是水分解制氢过程的核心基元步骤，其活化机理的理解对高效催化剂设计至关重要。氢能被视为解决全球变暖、环境污染和化石能源危机的潜在清洁能源载体，但水分解反应的吉布斯自由能变（ΔG° = 237.13 kJ/mol）为正值，动力学势垒显著，通常无法自发进行。水的三种振动模式中，对称伸缩（ν₁）和反对称伸缩（ν₃）通过Darling–Dennison共振强耦合，而弯曲模式（ν₂）与伸缩模式耦合较弱。根据Polanyi规则，水离解属“晚期”势垒反应，振动激发可有效促进反应进行。尽管已有大量理论和实验研究关注水在金属表面的离解，但只有少数研究将吸附水分子的红外光谱化学位移作为反应性描述符。研究人员此前在Cu(10?n)Ptn和Cu(9?n)Pdn合金催化剂上引入了ν₁/ν₃强度加权中心的概念，发现其红移程度与离解势垒相关。本研究系统考察了88个合金团簇上吸附水的红外振动化学位移、电子结构与离解温度之间的关联，旨在检验该描述符的普适性并发展预测性理解。研究采用密度泛函理论（DFT）计算和第一性原理分子动力学（FPMD）模拟，揭示了振动红移、HOMO–LUMO能隙匹配以及前沿轨道参与度对水离解活性的综合影响。该工作发表于《Molecules》杂志，为模型合金团簇上的水活化机理提供了原子尺度认识，并提出了可用于评估本征催化活性的光谱描述符。

研究人员采用的关键技术方法包括：基于BIOVIA DMol3程序（2019）的密度泛函理论（DFT）计算，使用广义梯度近似（GGA）下的Perdew–Burke–Erzerhof（PBE）泛函、相对论半芯赝势（DSPP）和双数值加极化（DNP）基组；所有团簇均带?1电荷，以负电荷孤立水分子H₂O^?为振动参考态；通过优化几何结构计算红外光谱，定义振动描述符ΔIR为吸附水ν₁/ν₃强度加权中心相对参考态的频移；采用NVE系综的FPMD模拟热离解，温度范围50–1500 K，离解判据为O–H距离超过1.5 ?；对部分体系执行三条独立轨迹以评估统计变异性。

**3.1 催化体系的识别与构建**
研究人员考察了元素周期表前75个单原子在正、中、负三种电荷态下吸附H₂O的ν₁/ν₃中心红移，发现负电荷态红移显著大于其他电荷态，且随元素性质呈周期性变化；惰性元素在负电荷态因易于向水分子转移电子而表现出高活性。以负电Cu原子红移615 cm^?1为参考，构建了24个双原子、30个三原子和34个四原子合金团簇，主要包含Pt、Na及红移超过Cu的元素。在88个初始团簇中，84个成功优化并用于后续计算，4个体系（NaPt-H₂O、CdPt-H₂O、Ni₂Pt-H₂O、Ni₃-H₂O）因几何优化不收敛而排除。

**3.2 离解温度对红外振动光谱的依赖**
通过FPMD模拟获得各团簇上H₂O的离解温度T_d，并将其与ΔIR关联。结果显示，ΔIR小于300 cm^?1的团簇离解温度不超过200 K；ΔIR在300–700 cm^?1之间（除Sr₃Pt和CuPt外）均表现出较高活性，离解温度不高于600 K；ΔIR在700–1016 cm^?1范围内催化活性差异显著，而更大ΔIR对应体系未观察到明显活性。这表明ν₁/ν₃中心红移是描述水催化离解活性的有效光谱描述符，与Polanyi规则和SVP模型定性一致。

**3.3 吸附能与HOMO–LUMO能隙的能量匹配**
研究人员分析了Kohn–Sham HOMO–LUMO能隙（E_g）与吸附能（E_a）对离解温度的影响。数据表明，当E_g + E_a接近零时，离解温度往往较低。统计相关分析显示，排除异常点前T_d与|E_g + E_a|的Pearson相关系数为0.78；按离群值判据排除SrRe₃和FPt后，相关系数提升至0.90。当|E_g + E_a|超过0.25 eV时，反应通常不利。该结果提示吸附能释放与电子激发能匹配可能以共振样机制促进离解，但直接证据需依赖非绝热动力学或激发态计算。

**3.4 前沿分子轨道的作用**
通过比较前驱体和中间态的HOMO-1、HOMO、LUMO，发现大多数体系的前沿轨道构型在反应前后基本一致，仅少数发生能级顺序调整。参与反应的前沿轨道定义为在离解过程中与H原子重叠增加的轨道，其数目与离解温度呈反比。Mulliken电荷分析进一步支持这一判断，例如CuPt仅LUMO参与且无能量匹配，离解温度高于1500 K；SrPt满足能量匹配条件，且HOMO和LUMO均参与H离解，离解温度较低；Sr₂Re通过HOMO-1和LUMO的能级变化参与催化，降低了离解温度。该结果说明前线轨道参与程度是影响水活化的重要电子结构因素。

讨论部分强调，上述相关性均基于平衡态基态性质，FPMD模拟从优化后的吸附态出发，未直接模拟吸附能释放后的能量再分配过程；因此提出的共振机制应视为需进一步验证的初步假说，明确验证还需非绝热分子动力学、过渡态计算或自由能分析。研究结论可翻译为：振动描述符ΔIR可作为评价水分子催化活化的光谱描述符，但全面理解活化机理还需考虑前沿轨道参与度和熵效应。结果暗示吸附能与HOMO–LUMO能隙接近可能与催化活性相关，此类共振样电子效应需非绝热或激发态计算确认。这些发现为模型合金团簇上的水离解提供了理论基础，并提出可能指导未来催化剂设计的光谱描述符；但计算基于孤立气相阴离子团簇，推广至实际负载催化剂时需考虑载体、溶剂和真实反应条件等因素。