高温度分辨率??Fe穆斯堡尔谱学精细表征铁矿物的采集与分析

《ChemPhysChem》：Collection and Analysis of High Temperature-Resolution 57Fe M?ssbauer Spectroscopy for Detailed Iron Mineral Characterisation

　　

??Fe穆斯堡尔谱学（M?ssbauer spectroscopy）是确定含铁材料氧化态和磁性质的有力技术。由于磁性质依赖于温度，穆斯堡尔谱（M?ssbauer spectra）通常需要在多个温度下采集，尤其是在磁性阻塞温度（magnetic blocking temperature）等特征转变点附近。穆斯堡尔数据采集时间缓慢，有时需要超过24小时才能达到高信噪比（signal-to-noise ratio，SNR），这限制了可采集谱图的数量。本研究提出一种替代方法，即测量具有高温度分辨率（如1 K温度步长）但较低信噪比的谱图。研究人员以不同结晶度的水铁矿（ferrihydrite）为例，在不显著增加总采集时间的情况下揭示了更详细的矿物学性质。超精细参数（hyperfine parameters）的温度依赖性在定性分析中被清晰揭示，而在温度和能量域中同时进行的定量分析则揭示了纯水铁矿和混合相样品中颗粒结晶度分布的细节。该方法将改进温度依赖性矿物性质的测量，能够分析阻塞温度分布，改善超精细参数在温度剖面中的插值，并有助于区分混合相样品中重叠的谱图。

高温度分辨率??Fe穆斯堡尔谱学：铁矿物精细表征的方法创新与矿物学应用解读

一、研究背景与问题

穆斯堡尔谱学（M?ssbauer spectroscopy）是基于穆斯堡尔活性同位素的无声反冲γ射线发射与共振吸收的物理化学分析技术，能够精确分析目标同位素的电子态、氧化态和原子间相互作用。其中??Fe是最受关注和研究最广泛的穆斯堡尔活性同位素，已在材料科学、凝聚态物理、应用物理和化学等领域得到广泛应用，并在化学工程和地球与行星科学等方向持续增长。

??Fe穆斯堡尔谱学的重要应用之一是测定含铁材料的磁性质。铁离子的电子自旋构型决定了材料的磁性质，而这些性质强烈依赖于温度——热能会逐步破坏磁有序。在高温下净磁场平均为零，谱图呈现双峰（doublet）或单峰（singlet）特征；在低温下则呈现六线峰（sextet），根据内磁场矢量与电场梯度主轴之间的夹角产生六条或八条共振线。两种状态之间的转变发生在阻塞温度（blocking temperature，T_B）处，即样品中半数Fe的弛豫时间超过测量特征时间的温度，操作上定义为谱图中半数呈现六线峰的温度。T_B通常低于材料的磁有序温度，且依赖于磁畴的形状和尺寸，因此是分析材料结晶度的有用参数。此外，对含有混合相（各相分别处于其T_B之上或之下）的样品进行分析，可用于区分具有相似超精细参数（hyperfine parameters）或谱图重叠而无法精确拟合的物相。

然而，利用穆斯堡尔谱学确定T_B通常需要采集T_B上下的多个谱图。对于具有显著粒径分布的矿物相，在T_B附近测量的谱图包含塌缩六线峰和双峰的混合物，可以从有序和无序位点的比例随温度的变化中插值得到T_B。但在大多数实际情况下，现有方法无法精确分析平均T_B或其分布。传统穆斯堡尔谱数据采集时间缓慢，高信噪比（signal-to-noise ratio，SNR）谱图有时需要超过24小时，严重限制了可采集的温度点数，导致温度分辨率不足。这一瓶颈制约了对矿物磁性质、结晶度和混合相组成的精细表征。

二、研究内容与意义

针对上述问题，研究人员提出了一种时间高效的高温度分辨率穆斯堡尔温度剖面分析方法。该方法利用超精细参数的温度依赖性来补偿较低信噪比测量的不足，在定性和定量分析中同时考虑所有温度下的谱图。核心思路在于：根据信噪比与测量循环次数的标度关系，信噪比的增长与√N成正比（N为测量循环次数），因此测量时间加倍仅带来41%的信噪比提升，十倍测量时间也仅带来216%的提升——大部分谱学信息在测量的初始阶段就已获得。研究人员将节省的计数时间用于提高温度分辨率，以1 K温度间隔代替常规的10 K间隔，在相同总采集时间内获得更丰富的温度依赖性信息。

研究人员以不同结晶度的合成水铁矿（ferrihydrite）为案例验证了该方法。高温度分辨率数据在定性分析中清晰揭示了超精细参数的温度依赖性，在定量分析中通过同时在温度和能量域拟合谱图，揭示了纯水铁矿和混合相样品中颗粒结晶度分布细节。该方法能够分析阻塞温度分布，改善超精细参数在温度剖面中的插值，并有助于区分混合相样品中重叠的谱图。该方法对理解铁氧化物/氢氧化物矿物的结晶度、粒径分布和磁性转变行为具有重要价值，在环境地球科学、土壤科学和环境矿物学等领域具有直接应用前景。

三、关键技术方法

研究人员采用了以下主要技术方法：第一，自动化穆斯堡尔谱采集协议，通过自定义Python代码控制WSS-10谱仪（WissEL，德国），实时计算信噪比（SNR，最强共振峰与谱图两端噪声标准差之比并经谱面积归一化），达到预设质量目标后自动降温并重新采集；第二，扩展Voigt基拟合（Extended Voigt-Based Fitting，xVBF）线形模型，用于拟合单个温度谱图及温度剖面的静态弛豫分析；第三，xVBF-Wickman修正模型，在xVBF基础上引入与磁弛豫时间τ成比例的内四线洛伦兹线宽展宽；第四，基于自旋翻转随机两相模型的Blume–Tjon动态弛豫模型，通过SYNCmoss软件实现；第五，温度依赖性参数建模，包括中心位移（centre shift，δ）按二阶多普勒位移计算、超精细场（hyperfine field，B_hf）按布里渊函数（Brillouin function）结合Bean–Rodbell模型计算、谱面积按Lamb–M?ssbauer因子归一化；第六，所有温度谱图同时拟合的定量分析协议，基于Scipy optimize（Scipy v1.16.3）包实现。样本为??Fe富集的合成水铁矿样品，包括6线水铁矿（6L-Fh）、2线水铁矿（2L-Fh）和低结晶度2线水铁矿（2Lc-Fh，在2线水铁矿合成过程中悬浮10 g/L土壤制备），以及??Fe富集针铁矿（goethite）与水铁矿的混合样品。

四、研究结果

高温度分辨率与低温度分辨率剖面的对比

研究人员以??Fe富集2L-Fh为对象，对比了相同总采集时间下的高、低温度分辨率剖面。低分辨率剖面以10 K间隔、每温度30分钟采集；高分辨率剖面以1 K间隔、每温度3分钟采集，总采集时间相同。将高分辨率数据以能量-温度矩阵形式绘制，用颜色标度表示相对吸收强度。10 K分辨率在此例中足以将T_B估计到几K以内，因为T_B分布较宽。但当包含双峰和六线峰混合的谱图数量太少时，拟合精度下降。低温度分辨率数据仅有少量谱图包含六线峰和双峰的混合物，且难以在缺乏相邻温度谱图上下文的情况下明确分配面积，对T_B估计有显著影响。

定性分析揭示样品间磁性质差异

研究人员对三种水铁矿样品（6L-Fh、2L-Fh和2Lc-Fh）进行了高温度分辨率测量。6L-Fh具有较高的T_B，与较高的结晶度一致，且T_B分布较大，与较大的粒径范围一致。相比之下，2Lc-Fh的双峰持续到低于主磁有序区域（约40 K）的温度，表明其颗粒结晶度分布可能不是高斯分布，或形成了具有不同结晶度的多个水铁矿相。这种不规则分布可能是部分颗粒在晶体形成过程中与土壤组分相互作用所致。

混合相样品中低丰度相的鉴别

研究人员将16%??Fe富集针铁矿与84%??Fe富集水铁矿的混合物以高温度分辨率测量。在能量-温度矩阵中可以清晰区分两种物相，即使单个温度谱图的信噪比较低、无法在单谱中识别低丰度相，矩阵视图仍能揭示其存在。多个T_B值的观察有助于揭示样品中存在多少物相，从而为单温度谱图的分析提供依据。图像处理技术也可用于自动识别谱图特征。

三种拟合模型的定量比较

研究人员对2L-Fh在5–140 K温度范围和?12至+12 mm/s速度范围内同时拟合所有135个温度谱图，比较了三种线形模型：xVBF静态模型、xVBF-Wickman修正模型和Blume–Tjon动态弛豫模型。三种模型均基于T_B分布叠加多个子谱来模拟部分有序六线峰的不对称展宽。xVBF模型给出T_B=64.23 K，xVBF-Wickman模型给出62.97 K，与单谱logistic拟合的T_B（67.25 K）最为接近，且这两种模型具有最低的χ2（分别为2.78和2.37）和最低的全局NRMSE（分别为2.70%和2.49%），能最好地再现5 K、110 K和T_B附近的谱图特征。Blume–Tjon模型则无法捕捉饱和六线峰形状、塌缩六线峰宽度和T_B附近双峰间距等关键特征，导致T_B被高估至78.61 K，χ2高达13.58，全局NRMSE为5.73%。三种模型拟合的T_B分布标准差（σ(T_B)）为9–12 K，表明2L-Fh具有较大的T_B分布。

拟合参数与模型局限

三种模型拟合的超精细参数在关键参考温度下与既往研究中Voigt基拟合和xVBF拟合的结果相当。所有模型中超精细场温度依赖性参数ν（约0.3）拟合值偏低，可能表明B_hf从0 K饱和值到T_B处零场的转变较为平缓，与水铁矿纳米颗粒表面效应一致。但残差模式表明B_hf随温度的梯度被低估。德拜温度（Debye temperature，θ_D）的拟合值（117–168 K）低于铁蛋白等类似材料，受纳米颗粒边缘效应、低质量反冲和厚吸收体饱和效应影响，且对超精细参数的温度依赖性变化非常敏感。三种模型均无法完全捕捉T_B附近部分塌缩谱图的温度依赖性，可能因为未考虑水铁矿的不规则磁行为。

计算成本与实用考量

定性分析的计算成本可忽略不计，定量静态弛豫模型在标准工作站（Intel Core i7，8核，32 GB RAM）上各运行不到2小时，而动态弛豫模型需要15小时计算。可通过减少迭代次数、降低数值积分分辨率、缩小温度范围等方式降低计算成本。高温度分辨率剖面通常替代具有较多高信噪比谱图的低温度分辨率剖面。采集时间应根据分析目的和样品性质逐案决定：??Fe富集样品低信噪比谱图采集3分钟，天然同位素丰度样品（??Fe含量2.119%）通常每温度至少需要30分钟。用户需在信噪比、采集时间、温度范围和温度分辨率之间进行权衡。

五、讨论与结论

研究人员指出，拟合参数反映的是纳米颗粒的功能性质而非块体材料性质，可用于理解矿物结晶度等性质，但各参数值可能受多种物理因素耦合影响。与单能量穆斯堡尔谱拟合常需使用物理上不合理的参数不同，基于温度的拟合更依赖于模型对样品物理特性的适用性。不能在所有温度下捕捉关键材料性质的模型会产生拟合妥协，这些模型不足同样影响单温度谱图的拟合。水铁矿的磁阻塞行为仍需进一步研究。塌缩谱图的建模仍是争论中的问题，弛豫模型的选择取决于样品性质。随着穆斯堡尔谱学和计算能力的日益普及，温度剖面采集和分析方法的创新将不断推进。

研究结论为：高温度分辨率穆斯堡尔谱学测量通过重新平衡信噪比与温度分辨率之间的权衡，为含铁材料的磁性质提供了独特视角。无论定性还是定量分析，高温度分辨率数据都能提供关于材料结晶度的额外清晰度，并可能揭示复杂矿物混合样品中的低丰度相。自动化协议实现了时间高效的测量，在不增加仪器时间需求的情况下提高了温度分辨率。作为一种严格表征含铁材料和矿物混合物的工具，高温度分辨率穆斯堡尔谱学有望为多个科学学科的研究项目提供独特价值。