通过联合中子衍射和计算技术精修α-FeOOH中的氢位置

《Dalton Transactions》：Refining hydrogen positions in α-FeOOH through combined neutron diffraction and computational techniques

　　

研究人员通过中子衍射研究了铁锈的关键组分针铁矿α-FeOOH的氢位置和磁结构。利用不可约表示和磁空间群方法，分析了空间群Pnma、磁波矢q_m = (0, 0, 0) r.l.u.下所有对称允许的磁结构。磁矩沿b轴排列形成反铁磁自旋排布，这一结果被第一性原理计算所复现。准确确定氢位置对于理解α-FeOOH中CO₂催化还原的机制至关重要。通过衍射和计算精确确定了这些位置，突显了将这两种方法联合用于非氘代化合物的有效性。

**论文解读：联合中子衍射与第一性原理计算精修α-FeOOH中的氢位置**

**1. 研究背景、现存问题与研究动机**

氢是材料科学中极具重要性的元素，其小原子半径使其易嵌入多种无机材料的晶格，显著改变材料的物理化学性质。针铁矿（α-FeOOH）是自然界最普遍的铁氧氢氧化物矿物之一，也是铁锈的主要成分，其晶体结构、磁性质和表面化学已在矿物学、固态化学和环境科学中得到广泛研究。近年来，α-FeOOH被发现可作为光化学还原CO₂为甲酸（HCOOH）的催化剂（在[Ru(bpy)₃]²⁺作为氧化还原光敏剂条件下），且其产HCOOH的速率和选择性显著高于α-Fe₂O₃，这表明α-FeOOH中的羟基（OH^-）基团起着关键作用。理论计算和实验均证实α-FeOOH表面OH基团在CO₂吸附及质子耦合电子转移反应中扮演重要角色。因此，精确获取α-FeOOH体相中的氢原子位置和氢键几何结构，对于约束表面质子环境、模拟表面终止以及质子迁移步骤至关重要，即使催化反应发生在表面。

然而，X射线衍射因氢原子电子密度极低而难以定位氢。中子衍射虽能克服该局限，但氢的非相干散射截面大且散射长度为负（?3.7390 fm），导致衍射信号复杂，尤其对于非氘代样品，精准定位氢仍具挑战。此前对α-FeOOH的中子衍射研究虽通过标准Rietveld精修确定过氢位置，但氢的精确度依然存疑。因此，本研究的动机是：通过联合中子衍射和第一性原理计算，以互补方式精确确定α-FeOOH中氢的体相位置，并同时测定其磁结构，为理解CO₂催化还原机制提供关键结构信息。

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

本研究采用商业α-FeOOH粉末（Kojundo Chemical Lab. Co., Ltd, Japan, 99%+）作为样品，在日本研究反应堆（JRR-3）的HERMES高分辨率中子粉末衍射仪上采集数据（波长λ = 1.34239(7) ?，温度范围60–485 K）。使用FullProf Suite进行Rietveld精修，采用pseudo-Voigt峰形函数和Chebyshev多项式背景。第一性原理计算采用投影缀加波方法（PAW），基于VASP代码，采用PBEsol型交换-关联势，截断能550 eV，并采用DFT+U方法（U = 5.3 eV）处理Fe 3d电子的强关联效应。通过CLUPAN代码识别对称独立的反铁磁构型，共计算了58种自旋构型。晶体结构完全弛豫至原子间残余力小于0.02 eV ?^-1。

**3. 研究结果**

**确定磁转变温度与磁结构**
研究人员首先通过变温中子粉末衍射（64.8(1) K和430.1(9) K）确定α-FeOOH的磁转变温度。从200和101衍射峰的积分强度温度依赖性拟合得到Néel温度T_N = 382(2) K。利用不可约表示（irrep）分析（Kovalev符号），在空间群Pnma、磁波矢q_m = (0, 0, 0) r.l.u.下，识别出8个一维不可约表示对应的12个基矢（BV），并列出所有最大磁空间群（MSG）。通过比较R_mag因子（64.8(1) K数据），最佳拟合为MSG Pnma′ (#62.445)，对应Γ₈的ψ₁₂（R_mag = 8.53%），磁矩为3.64(4) μ_B per Fe³⁺ site。该磁结构为沿b轴排列的简单反铁磁结构，磁单胞与晶胞相同。这一结果与以往在Pbnm设定下的报道一致（通过轴变换等价）。

**氢位置的精修与温度依赖性**
在顺磁相430.1(9) K下，衍射图谱符合空间群Pnma的α-FeOOH结构。Rietveld精修成功获得O2–H键长随温度的变化（图2(d)），其在0.96 ?（水中的O–H键长）附近略有增大（1.008(7) ? at 64.8 K），整个温度范围变化<5%，未观察到质子有序-无序转变或平均结构对称性变化。还评估了H?O1距离、O1?O2供体-受体距离和O2–H?O1角随温度的变化（图S1），以描述氢键几何的演化。

**计算与实验对比**
第一性原理计算获得的最稳定反铁磁构型的磁化密度分布（Δρ^↑↓）显示Fe–Fe自旋沿[100]方向反铁磁、沿[010]方向铁磁，与实验精修磁结构一致。优化晶格常数a = 9.885, b = 3.010, c = 4.542 ?，与64.8(1) K实验值差异<1%。O2–H计算键长为1.02253 ?，与实验值1.008(7) ?相差1.5%。氢键几何参数（H?O1、O1?O2、∠O2–H?O1）的计算值与实验值一致，给出近乎线性的氢键排布。投影态密度（PDOS）计算显示α-FeOOH带隙为2.16 eV，价带顶和导带底主要由Fe(3d)和O(2p)轨道构成，OH基团的分子轨道位于较深能级。

**方法有效性验证**
鉴于实验与计算的近乎完全一致，研究人员证明对非氘代样品联合使用中子衍射和第一性原理计算可有效确定氢位置。对于本高对称性化合物，精修相对直接；但对于更复杂结构或更多Wyckoff位置的化合物，氘代可能仍属必要。成功精修的氢位置为理解CO₂还原和析氧机制提供了最小必要信息。

**4. 讨论与结论总结**

研究人员指出，表面板模型通常基于体相结构构建，表面O位点的初始质子化模式继承自体相OH取向和氢键拓扑。如果体相H位置约束不佳，则需假设多种可能的氢键网络和表面羟基构型，这将导致对CO₂还原中吸附几何和质子耦合电子转移能的不同预测。本精修提供了明确的体相参考（包括温度依赖的氢键几何），可作为反应建模的一致起点，并用于评估缺陷或表面终止如何扰动质子环境。未来方向包括氢动力学的激发研究以及CO₂气氛下的原位测量。

**研究结论**（原文翻译）：
综上所述，研究人员精修了针铁矿α-FeOOH的磁结构和氢位置。基于不可约表示和磁空间群分析，在零场下利用粉末样品成功精修了磁结构。通过Rietveld精修和第一性原理计算准确确定了氢位置。实验与计算结果的一致性突显了同时使用这两种方法的有效性，即便对于非氘代样品也是如此。该研究成果发表在《Dalton Transactions》期刊上。