γ-Al2O3表面晶面依赖的MoO3前驱体吸附能景观：对MoS2基加氢脱硫催化剂的可能启示

《Computational Materials Science》：Facet-dependent adsorption energy landscapes of the MoO3 precursor on γ-Al2O3: Possible implications for MoS2-based hydrodesulfurization catalysts

　　

采用密度泛函理论（DFT）映射方法，研究了MoO₃前驱体在非尖晶石结构γ-Al₂O₃的脱羟基化(001)面和真实羟基化(110)面上的锚定行为。通过自动化二维空间映射策略对势能面（PES）进行全面采样，以最大限度减少前驱体初始放置带来的偏倚。结合吸附能、几何约束和电荷密度差（CDD）分析，建立了理解表面晶体学如何调控吸附能景观和Mo前驱体空间定位的第一性原理框架。结果表明，脱羟基化(001)面呈现相对连续的吸附能景观，前驱体通过双齿锚定模式与配位不饱和表面原子作用，形成较短的Mo–O键（1.98 ?），并伴随定域化电荷转移（纯电子相互作用能为?489 kJ/mol）；而羟基化(110)面由于致密羟基覆盖层的空间位阻效应（能垒超过250 kJ/mol），吸附景观高度受限，前驱体仅能容纳于特定吸附域（Sites 12和16），其稳定化机制转变为由一个Mo–O配位和三个氢键构成的多点相互作用网络，纯电子相互作用能为?380 kJ/mol。键价分析表明，(001)面和(110)面上Mo–O_surf键的有效键级分别为0.82 v.u.和0.63 v.u.，与纯电子相互作用能的变化趋势一致。研究揭示，γ-Al₂O₃载体作为活性结构组分而非被动载体，通过晶面特异性的几何和电子约束引导前驱体的局域化分布，为理解和控制MoS₂基加氢处理催化剂中活性相的分散提供了原子尺度依据。

加氢处理（HDT）是现代石油炼制的核心催化过程，用于脱除硫、氮等杂原子以生产清洁运输燃料并减少环境污染。工业上广泛使用的催化剂为负载于γ-氧化铝（γ-Al₂O₃）上的二硫化钼（MoS₂），通常以钴或镍为助剂。经典rim-edge模型以及扫描隧道显微镜（STM）和理论模型均证实，催化性能与活性相的形貌性质密切相关——催化反应主要发生在配位不饱和的边缘位点而非基面，因此控制MoS₂纳米片的分散度、长度和堆叠程度是最大化催化活性的必要条件。γ-Al₂O₃载体并非被动物理承载者。Kn?zinger和Ratnasamy提出的经典经验模型强调了表面羟基的作用，预测载体的晶体学和化学性质主动引导活性相的成核和最终形貌。多晶γ-Al₂O₃暴露多种晶面，其中(110)和(001)面占暴露面积的绝大部分。值得注意的是，Pigeon等人通过勃姆石拓扑转变详细阐述的现代晶体学重新指派已更新了历史命名，将传统(100)面重新定义为(001)面。表面科学实验和光谱研究已证实催化剂制备过程中存在显著的结构敏感性。Bara等人利用模型平面氧化铝基底的开创性研究提供了前驱体-载体相互作用高度依赖晶面的实验证据：在(110)等效晶面上Mo物种被阻碍完全硫化而截留在非活性氧化态；而目前定义为(001)的晶面通常为活性相的聚集和硫化提供更有利的环境。然而，这些不同的宏观行为表明，支配这些相互作用的确切物理化学规律需要更深入的基础研究。第一性原理计算方法，特别是密度泛函理论（DFT），已被广泛用于研究催化剂浸渍和焙烧早期阶段钼前驱体（通常以MoO₃建模）初始吸附的原子尺度机制。早期理论研究为Mo在氧化铝表面的吸附提供了有价值的见解。然而，许多后续研究，包括Yu等人对α-MoO₃/γ-Al₂O₃复合材料的研究，对所有氧化铝晶面均采用理想化脱羟基模型。这些模型虽推进了对前驱体吸附的理解，但忽略了实际催化剂制备条件下不同γ-Al₂O₃表面所特有的不同羟基化状态。完全脱羟基氧化铝表面的广泛使用往往不能代表标准催化剂制备条件。γ-Al₂O₃表面复杂的氢键网络和振动性质已被充分表征，凸显了采用真实结构模型的必要性。非尖晶石模型能更忠实地再现局部配位环境。在此结构框架基础上，Ferreira等人报道的从头算热力学相图揭示了各γ-Al₂O₃晶面不同的水合行为：(001)面在相对较低的焙烧温度下即可脱水，暴露配位不饱和的铝和氧位点；而(110)面在相同条件下保持致密且热力学稳定的羟基覆盖层。因此，忽略这些晶面依赖的羟基化状态可能导致对界面相互作用和空间约束的描述不完整。为解决这一空白，研究人员引入系统性DFT映射方案，研究MoO₃前驱体在非尖晶石γ-Al₂O₃脱羟基化(001)面和真实羟基化(110)面上的锚定机制。为确保势能面（PES）的全面采样并减少与前驱体初始放置相关的偏倚，对两个晶面均采用自动化二维空间映射方法。通过结合吸附能、几何约束和电荷密度差（CDD）分析，建立了理解表面晶体学如何调控吸附能景观和Mo前驱体空间定位的第一性原理框架。结果进一步表明，(110)面上的羟基覆盖层将锚定机制从定域双齿配位转变为多点氢键网络，揭示了MoO₃前驱体在γ-Al₂O₃上晶面依赖吸附行为的原子尺度起源及其对MoS₂基HDT催化剂制备的启示。

研究人员开展研究用到的主要关键技术方法包括：密度泛函理论（DFT）计算，采用Quantum ESPRESSO软件包，广义梯度近似（GGA）结合Perdew–Burke–Ernzerhof（PBE）泛函，并应用Grimme的D3半经验色散校正；采用投影增强波（PAW）方法描述电子-离子相互作用；基于Ferreira等人验证的非尖晶石结构模型构建γ-Al₂O₃的(001)和(110)晶面平板模型，羟基化状态依据从头算热力学相图指派；采用自动化二维空间映射协议对势能面进行系统采样（(001)面56个初始构型，(110)面81个初始构型）；通过电荷密度差（CDD）分析和冻结碎片相互作用能计算分离电子相互作用贡献；采用Pauling–Brown键价描述符评估有效键级。样本队列来源为基于文献报道的非尖晶石γ-Al₂O₃结构模型和从头算热力学相图，不涉及实验样本队列。

研究结果方面：3.1 脱羟基化(001)面上的吸附能景观——通过对56个系统分布的初始构型进行几何优化，结构收敛为10个不同的吸附极小值。最稳定区域由两个近简并构型（Sites 3和6）组成，以Site 3为零能量参考。收敛结构的相对能量跨度约105.0 kJ/mol，近40%的初始构型弛豫至这两个最低能量极小值，表明存在明确的能量盆地。研究通过自上而下能量等高线图和全局最小值构型分析得出，能量差异主要源于前驱体取向和局部配位变化而非载体显著形变。邻近吸附位点Site 3和Site 4（间距小于1 ?）的能量差约56 kJ/mol，Site 3中MoO₃优先结合表面氧原子最大化Lewis酸碱相互作用同时保持局部骨架完整，而Site 4涉及直接配位表面铝阳离子诱导晶格位移，表明Site 4的稳定性降低与局部结构畸变相关。3.2 羟基化(110)面上的吸附能景观——81个初始构型中51个成功收敛为16个不同吸附极小值，30个未满足收敛标准。研究通过空间约束分析和能量分布映射得出，未收敛并非方法学假象而是(110)面严重空间约束的直接物理表现：29个收敛结构已占据>250 kJ/mol的高度排斥能量平台，空间上环绕30个未收敛网格点，证实未收敛点位于同一空间排斥区的最陡排斥峰处，51个收敛构型完全捕获了(110)吸附景观上所有物理可及的盆地。吸附能景观分布超过270 kJ/mol，远宽于(001)面。Boltzmann加权分析（500 °C，RT≈6.4 kJ/mol）表明仅4个结构位于最低能量稳定窗口（0–10 kJ/mol），中间能量范围（30–65 kJ/mol）的构型贡献低于1%。研究通过结构分析得出，能量最优态中Mo原子直接配位四面体铝位点（Al^IV）上的羟基并受相邻氢键网络稳定，前驱体须定位表面质子密度中的空间裂隙才能穿透；亚稳态中Mo与八面体铝位点（Al^VI）上的氧作用但局部几何强制次优取向；高能排斥态中致密质子层限制进入，迫使分子与羟基发生Pauli排斥。3.3 载体效应对绝对锚定强度和空间约束的影响——最稳定锚定构型在(001)和(110)面上的绝对吸附能分别为21和24 kJ/mol，相似的温和吸热性与高温焙烧促进前驱体铺展和与载体表面作用一致。(110)面上最低能量吸附位点（Sites 12和16）空间邻近（间距1.7 ?），限制多个MoO₃单元同时容纳；(001)面上Sites 3和6间距2.4 ?（跨周期边界6.0 ?），位阻约束较小允许更广泛的吸附构型。结构参数表明(001)面Mo–O_surf键长1.98 ?，(110)面为2.08 ?，后者通过延伸氢键网络（1.71 ?、1.82 ?和2.11 ?）稳定前驱体。3.4 电子洞察——冻结碎片相互作用能表明(001)面E_int=?489 kJ/mol，CDD显示高度定域电荷重分布，双齿锚定基序将电荷转移限制在局域空间域；(110)面E_int=?380 kJ/mol，CDD呈更弥散拓扑模式。通过局部去质子化控制计算（E_int从?380增强至?439 kJ/mol），研究定量分离出局部O–H质子化贡献59 kJ/mol（约15%）的电子弱化，其余50 kJ/mol差异源于本征晶体学位点几何（(001)面相邻不饱和Al阳离子可实现双齿螯合基序）。3.5 构效关系与竞争吸附——通过Pauling–Brown键价描述符s=exp[(R₀?d)/B]评估有效键级，(001)面短Mo–O_surf键（1.98 ?）对应s=0.82 v.u.，(110)面伸长键（2.08 ?）降至s=0.63 v.u.，约23%的键级降低与约22%的纯电子相互作用能降幅定量吻合。研究通过竞争吸附分析得出，(110)面结合严重依赖羟基网络中的特定表面裂隙，竞争物种较易破坏这些弥散相互作用网络，而(001)面刚性双齿键可能提供更强的抗竞争置换能力。

讨论部分总结：研究人员将纯电子相互作用能、界面电荷转移拓扑和初级Mo–O_surf连接的有效键级相关联，建立了前驱体锚定稳定性的定量构效关系。电子和键级参数可作为前驱体锚定的稳健理论描述符，这与数据驱动催化科学中利用键级、电荷分布和局部配位等电子和结构描述符桥接微观吸附机制与宏观催化性能的趋势一致。竞争吸附方面，实际制备或反应环境中Mo前驱体须与残余水或硫化剂（如H₂S和SO₂）等配位物种竞争锚定位点。由于(110)面结合严重依赖羟基网络中的特定表面裂隙，竞争物种可能较易破坏这些弥散相互作用网络，而(001)面刚性双齿键可能提供更强的抗竞争置换能力。活性位点对Mo前驱体相对于竞争剂的选择性将由热力学稳定性（氢键网络强度）和动力学可及性（穿透羟基层空间约束）之间的微妙平衡决定。研究结论翻译：本工作采用密度泛函理论（DFT）映射方法研究了MoO₃前驱体在γ-Al₂O₃脱羟基化(001)面和羟基化(110)面上的晶面依赖锚定行为。第一性原理结果表明，吸附能景观和前驱体的潜在空间分布强烈受载体表面晶体学和化学性质的调控。在脱羟基化(001)面上，MoO₃吸附呈现相对连续的能量景观，不同构型主要由前驱体取向和局部表面弛豫决定。前驱体与配位不饱和表面原子作用，形成定域双齿锚定基序，包括较短的Mo–O相互作用（1.98 ?）和前驱体端氧与表面铝位点的额外相互作用。CDD分析揭示与强前驱体-载体相互作用（?489 kJ/mol）相关的定域电荷重分布模式。结合多个吸附构型的可及性，该表面环境可能为MoS₂活性相形成中涉及的前驱体再分布和后续聚集过程提供充足的空间灵活性。相比之下，羟基化(110)面因致密表面羟基层的存在呈现高度受限的吸附景观。羟基覆盖层产生显著空间位阻，形成能量惩罚超过250 kJ/mol的不利区域，将前驱体容纳限制在特定吸附域（Sites 12和16）。虽然绝对吸附能与(001)面相当，但稳定化机制根本不同。在(110)面上，前驱体与载体保持直接Mo–O配位，但表面羟基的参与改变了该相互作用的电子性质，并通过分布式氢键网络引入额外稳定化。CDD分析揭示，由一个Mo–O配位和三个氢键组成的四点相互作用基序提供强界面稳定化，同时限制可及吸附空间和相邻Mo物种的空间排列。总体而言，这些结果表明晶面依赖的前驱体行为不仅由绝对吸附强度决定，还由有效键级、电荷转移拓扑和空间相互作用网络决定，这些因素定义了可及吸附构型及其对竞争吸附的敏感性。因此，γ-Al₂O₃载体作为活性结构组分而非被动载体，通过晶面特异性的几何和电子约束引导前驱体局域化。这些原子尺度见解和构效描述符为理解和控制MoS₂基加氢处理催化剂中活性相分散提供了基础性指导。