和频产生（SFG）光谱与神经网络势（NNP）分子动力学联合解析水-Fe3O4界面反应性羟基网络及水取向失配

《Journal of Colloid and Interface Science》：Towards decoding the structure of water-magnetite interface

　　

联合振动光谱与机器学习驱动的分子动力学（Molecular Dynamics, MD）探测磁铁矿界面水结构。模拟揭示纳秒尺度上的反应性表面羟基网络。界面OH物种被归属到和频产生（Sum Frequency Generation, SFG）光谱中的光谱特征。解释理论与实验之间的水取向失配。结果挑战既有表面模型，并对替代水-磁铁矿终止结构施加严格约束。

论文主体解读

研究背景与问题：铁氧化物是地球丰度较高的矿物，具有多种晶相和氧化态，其中磁铁矿（Fe₃O₄）在地质风化、电化学、大气化学、催化反应、太阳能水分解和生物医学等领域具有重要作用。由于Fe₃O₄常与含水溶液接触，水-磁铁矿界面受到广泛研究，尤其是Fe₃O₄的稳定(001)和(111)表面。既有研究多基于超高真空（Ultrahigh Vacuum, UHV）条件和密度泛函理论（Density Functional Theory, DFT）计算，增进了对热力学稳定性、吸附、溶解和钝化趋势的理解。然而，大多数水-磁铁矿相互作用研究在UHV下进行，观察到部分解离水分子、预吸附表面羟基和有序铁羟基氧化物相；磁铁矿与宏观量水接触时的实验行为尚未被探索。计算研究也多限于经典模拟或短时间DFT模拟。此前研究人员开发了Behler-Parrinello架构的神经网络势（Neural Network Potential, NNP），用于暴露亚表面阳离子空位（Subsurface Cation Vacancy, SCV）重构表面的水-磁铁矿界面，并提供了解离机制和水在表面扩散的见解。本工作结合SFG光谱与NNP基MD模拟，解析磁铁矿-水界面的分子结构，识别界面OH物种并计算SFG响应，以解释实验光谱特征，并考察实验与模拟之间水取向失配的结构含义。

关键技术方法（不超过250字）：实验采用表面特异的相位分辨SFG光谱，在SSP偏振模式下测量O-H和O-D伸缩区；样品包括射频磁控溅射（Radio Frequency sputtering, RF-sputtering）制备的约50 nm Fe₃O₄/CaF₂薄膜，并用拉曼光谱（Raman spectroscopy）和X射线光电子能谱（X-ray Photoelectron Spectroscopy, XPS）表征，另使用外延生长的55 nm Fe₃O₄/MgO。计算采用NNP，通过n2p2与LAMMPS接口进行MD。模拟体系为含1152个O和800个Fe的磁铁矿板与2176个水分子接触，包含两个SCV界面，300 K，0.5 fs时间步；统计轨迹29 ns，光谱轨迹为四条独立1.3 ns轨迹。分析包括OH物种识别、取向概率密度、氢键长度和寿命自相关，并计算Im[χ⁽²⁾]谱。

研究结果

Experiment：研究人员通过实验与模拟联合分析界面水结构。实验部分首先制备并表征磁铁矿膜，然后进行SFG测量，最后从谱学数据推断界面水取向与氢键状态。

Fe₃O₄ layer preparation and characterization：约50 nm磁铁矿膜由RF-sputtering制备，拉曼光谱在不同位置相似，表明膜均匀。观察到约685 cm^?1的A_1g模式和约1343 cm^?1的双声子散射模式，约330 cm^?1的T_2g模式与CaF₂衬底320 cm^?1模式重叠。XPS Fe 2p区域与α-Fe₂O₃参考比较，Fe₃O₄中Fe²⁺和Fe³⁺混合导致Fe 2p_3/2和Fe 2p_1/2峰较赤铁矿略向低结合能移动，719 eV附近Fe³⁺相关振激卫星表明表面氧化。

SFG spectroscopy：在O-H伸缩区3000–3750 cm^?1测量磁铁矿/水界面，并用D₂O谱归一化。菲涅尔校正对整体谱形影响显著，并减弱3500 cm^?1以上未校正谱的宽贡献。校正后，H₂O在约3100 cm^?1处SFG比值大于1，在3300 cm^?1以上小于1；pH 11谱在3180 cm^?1以下大于1，以上小于1。相位分辨Im[χ⁽²⁾]显示，H₂O拟合包含一个低波数负带和两个高波数正带，分别代表较强和较弱氢键物种；负或正号指示O-H中H远离或朝向界面。pH 11下高波数带符号相反，中波数带幅度显著增加，符合带电界面特征。磁铁矿等电点约6.5–7，pH 11下表面可能带负电；盐筛选实验表明pH 7近中性，pH 11带电。OD区测量与OH区一致，MgO上外延Fe₃O₄层结果类似。

Insights into the water structure：H₂O谱在低波数近零，3030 cm^?1以上为正，提示OH基团以H朝向界面。pH 11谱在约3400 cm^?1以下为正，表明O-H基团以H朝向磁铁矿表面，以上有小的负信号。用非共振贡献和三个复洛伦兹（Lorentzian）峰拟合共振χ⁽²⁾，再现实验谱形。H₂O的共振虚部含一个低波数负带和两个高波数正带；pH 11中高波数带符号相反且中波数带增强，是带电界面附近水的典型特征。

Simulation：SCV表面直接暴露两层水分子，最外层为八面体铁和氧行，第二层为四面体铁行，每个表面单元缺失一个铁。水可在该位点解离，一个H结合表面氧，形成“解离OH”；剩余羟基称“remaining OH”，吸附于八面体铁；解离位点上方强相互作用的水称“water on dissociation”。分析第一水层至距表面3.8 ?。16个表面单元中平均有15.1个解离OH和15.0个remaining OH；13.4个位点发生一个解离，1.7个无解离，0.9个有两个解离；water on dissociation为15.2个。八面体水上每个单元3.2个，四面体水上3.1个。water on dissociation的取向与解离数目相关。

The orientation of the OH species at the magnetite interface：解离OH的OH键以cos θ=0.71（θ=45°）并垂直于铁行取向，φ=π/4。remaining OH主要垂直于铁行、略离表面，cos θ=0.25（θ=75.5°）；偏离垂直方向伴随θ增加，源于附近水氢原子的排斥作用。平行铁行的中间态来自水分子间氢转移。water on dissociation两个OH键形成夹角105°，一个指向表面，cos θ=?0.67（θ=132°），另一个平行且背离表面，cos θ=0.4（θ=66.5°）；偶极居中于cos θ=?0.22（θ=102.5°）。偶极向下或向上峰对应无解离或与两个解离OH形成氢键的情形。八面体铁上水的OH键均背离表面，峰在平行铁行略上方，cos θ=0.14（θ=82°），约0.34的OH键位于该峰，其余分布于峰周围。四面体铁上的水更加无序，偶极优先使H朝向表面。

Bond length analysis：OH键长与氢键长度联合概率显示，体相水中共价OH键长与氢键长度反相关。该关系在解离OH与water on dissociation之间、以及water on dissociation与表面氧之间更明显。remaining OH的相关较弱。体相水平均OH键长0.99 ?，氢键1.97 ?；解离OH平均OH键长1.02 ?，氢键1.65 ?；water on dissociation平均OH键长1.03 ?，氢键1.57 ?；remaining OH平均OH键长0.98 ?，氢键2.26 ?，说明其向四面体铁上方水提供很弱氢键。

Stability of the water-magnetite interface：通过指示函数自相关分析物种寿命。解离OH寿命约0.6 ns，主导过程为解离质子经解离位点上方水在相邻表面氧之间跳跃；remaining OH寿命约0.5 ns，质子转移发生在八面体铁位点吸附水之间；water on dissociation更稳定，寿命约1.9 ns，涉及不同构型转变和与周围水交换。界面反应性高，但寿命长于计算SFG所需皮秒尺度。

The simulated response spectra：计算Im[χ⁽²⁾]，对应实验SSP几何。第一水层总谱由约3200 cm^?1负峰主导，来自八面体水；四面体水在高波数产生窄正峰，约3500 cm^?1形成净表面取向响应。低于3000 cm^?1出现解离OH和water on dissociation向下OH峰，强氢键且平均取向相反，贡献符号相反。water on dissociation第二个H产生约3300 cm^?1负峰。约3600 cm^?1弱负贡献可能来自构型b的water on dissociation；remaining OH在相近波数有弱峰。加入更多水层后谱逐渐向正Im[χ⁽²⁾]移动，表明第一层外水仍取向有序并相干贡献。

Discussion：核心结果是实验推断的界面水净取向与模拟预测相反。实验从SFG符号和强度推断OH偶极朝向磁铁矿表面占优，而模拟预测OH基团H优先朝向体相。若模拟可靠，当前SCV终止可能不代表SFG实验样品。原文提出支持包括：表面X射线衍射显示Fe₃O₄(001)在mbar水蒸气下重构解除并表面粗化；液相接触可广泛羟基化并形成羟基氧化物相；宏观水可能引起Fe溶解、羟基化或重构，溶解氧可缓慢氧化。模拟重现关键光谱活性OH种类：低波数负信号来自解离OH与解离位点上方水的强相互作用，正贡献来自四面体铁行上方水。因此差异更多来自相对丰度和取向权重，而非单种物种错误。可能结构原因包括八面体铁吸附位点部分耗尽、类羟基氧化物相、较低解离度或更多解离位点。pH 11增强低波数强度与OH丰富且带负电界面一致；八面体吸附羟基朝体相产生高频特征，实验高频附加特征支持该解释。

Conclusions：在本工作中，研究人员将表面特异SFG光谱与反应性NNP基MD模拟结合，解析水-磁铁矿界面水和羟基物种的分子结构。界面水的振动响应不仅由吸附水和羟基决定，还由其相对种群、取向和氢键环境决定。结合实验与模拟，研究人员表明这些微观特征可直接联系到测量非线性光学响应的符号、光谱位置和强度。主要实验发现是水-磁铁矿界面在OH/OD伸缩区具有显著且pH依赖的SFG响应。近中性pH下盐筛选表明表面电荷诱导水有序的贡献可忽略，可直接与计算界面响应比较；pH 11下增强的低波数SFG强度与带负电、更羟基化界面一致。这些观察表明磁铁矿处水结构对表面质子化和电荷高度敏感，二者是胶体稳定性、腐蚀、环境反应性和铁氧化物催化过程中的关键变量。模拟提供光谱响应的分子解释：NNP-MD轨迹揭示异质且反应性的界面羟基网络，由解离表面OH、吸附于八面体铁位点的羟基、解离位点结合水及四面体行上方水组成。这些物种在取向、氢键强度和寿命上差异显著。计算谱将低波数响应归属于解离位点强氢键OH物种，高波数贡献来自较弱氢键羟基和界面水。与早期UHV、DFT和经典模拟研究相比，本工作将液体水条件下测得的表面敏感光谱与原子分辨的反应性界面构型直接联系。关键结果是实验推断的界面水净取向与模拟SCV终止Fe₃O₄(001)表面预测存在系统差异。该差异不意味着单个模拟物种非物理，而是表明其相对种群或取向权重与实验界面不同。这种不一致具有化学信息，约束了磁铁矿可能的水相终止。结果提示真实界面可能含有更大比例以H朝向表面取向的水分子，可能由于八面体铁位点部分耗尽、类羟基氧化物基序或水解离度低于理想SCV模型。未来工作应构建和测试替代水相磁铁矿终止，包括部分羟基化、缺铁和类羟基氧化物表面，并将反应性模拟扩展到显式长程静电和其他表面终止。实验上，系统测量pH、离子强度、同位素组成和表面制备将有助于分离局部界面OH特征与扩散层贡献。