超越Vegard定律：Gd2O3–UO2燃料短程结构非均匀性的原子尺度机制及其对热物理性能的影响

《ACTA MATERIALIA》：Beyond Vegard’s Law: Short-Range Structural Heterogeneity in Gd?O?–UO? Fuel and Its Impact on Thermophysical Properties

　　

氧化钆（Gd₂O₃）被广泛用作二氧化铀（UO₂）燃料中的可燃中子吸收剂（burnable neutron absorber），但支配其对热物理性能影响的原子尺度机制仍不明确。在此，研究人员通过在Gd掺杂UO₂（Gd质量分数6%–15%）中整合X射线衍射（X-ray diffraction, XRD）、扫描电子显微镜/能谱（scanning electron microscopy/energy-dispersive X-ray spectroscopy, SEM/EDS）、基于同步辐射的扩展X射线吸收精细结构谱（Extended X-ray Absorption Fine Structure, EXAFS）、共振超声谱（resonant ultrasound spectroscopy, RUS）以及热物理性能测量，建立了直接的结构-性能关系。XRD结果证实了单一相萤石固溶体（fluorite solid solution）的形成，其晶格线性收缩并遵循Vegard定律；而EXAFS则揭示了局部畸变、氧空位（oxygen vacancies）以及缺乏XRD检测所需长程结晶学序的类bixbyite纳米尺度环境。微观结构分析表明Gd分布均匀，无次生相或微裂纹，且较高的掺杂水平会降低晶粒尺寸并增加孔隙率。热物理测量表明，由于声子-缺陷散射（phonon–defect scattering）增强，Gd掺入显著降低了热扩散率（thermal diffusivity）与热导率（thermal conductivity），且在较低温度下尤为显著。共振超声谱进一步揭示弹性刚度（elastic stiffness）随Gd含量增加而逐渐降低，与缺陷诱导的晶格软化（lattice softening）一致。热阻分析表明，杂质散射（impurity scattering）在Gd掺杂UO₂中占主导地位，而声子-声子Umklapp散射（phonon–phonon Umklapp scattering）贡献随掺杂剂浓度增加而受抑制。这些结果表明，Gd掺杂UO₂在整体上是单相的但在局部是非均匀的，且这种短程结构复杂性支配着热导率行为。研究结果为基于原子尺度信息的燃料性能模型（atomistically informed fuel performance models）提供了关键实验依据，并强调了解析局部原子环境在评估掺杂核燃料中的重要性。

一、研究背景与意义

核能为全球提供了可靠的低碳能源。二氧化铀（UO₂）燃料驱动着目前全球90%以上的民用核反应堆，提升其效率与安全性是核工业的核心任务。为抑制辐照初期反应性、实现更高燃耗，可燃中子吸收剂（burnable absorber, BA）被添加至核燃料中。氧化钆（Gd₂O₃）因¹⁵⁵Gd和¹⁵⁷Gd的中子吸收截面大、燃耗速率与²³⁵U相当、且无非长寿命寄生副产物，成为极具吸引力的可燃吸收剂燃料（burnable absorber fuel, BAF）材料，已在压水堆（pressurized water reactor, PWR）和沸水堆（boiling water reactor, BWR）中应用。

尽管已有研究成功制备了高致密度Gd掺杂UO₂燃料并考察了Gd在其中的行为，若干关键问题仍未解决。例如Hirai等人曾报道UO₂燃料中出现微裂纹，但其潜在机制未被阐明；传统认知认为Gd与UO₂形成均匀固溶体，但Gd的局部结构可能因偏聚、畸变或短程有序而偏离理想取代模型。Gd³⁺作为大尺寸三价阳离子掺入UO₂晶格会引入电荷补偿缺陷，深刻影响燃料密度、缺陷化学、热力学与热输运行为。热导率是评估核燃料安全高效运行的关键参数，已有实验表明Gd浓度升高可使热导率下降高达50%，且低温下降更为显著，但声子散射的微观结构起源仍缺乏原子尺度的直接实验证据。

针对上述问题，研究人员在Gd掺杂UO₂（Gd质量分数6%–15%）中开展了多尺度结构表征与热物理性能研究，建立了明确的原子尺度结构-性能关系。核心结论为：Gd掺杂UO₂在长程尺度上是单一萤石相并遵循Vegard定律发生线性晶格收缩，但在原子尺度上存在短程结构非均匀性，即Gd³⁺取代诱发局部畸变、氧空位和类bixbyite纳米尺度局域环境；正是这种短程复杂性增强了缺陷散射，导致热导率显著下降、弹性刚度降低，并重塑了热阻中杂质散射与声子-声子Umklapp散射（Umklapp scattering）的相对贡献。该论文发表于《ACTA MATERIALIA》，为基于原子尺度信息的燃料性能模拟框架（如BISON、CENTIPEDE）提供了关键实验依据。

二、关键研究方法

研究人员将depleted UO₂粉末与Gd₂O₃按目标Gd质量分数6.3%、10%和14.9%称量，经高能球磨混合并添加0.25 wt.% EBS（Ethylene Bis-Stearamide）烧结助剂，在2023.15 K、Ar-H₂气氛下烧结制备芯块（另制备3%–14.9% Gd系列样品专用于XRD）。主要技术包括：采用Bruker D2 Phaser X射线衍射仪（Cu Kα₁辐射）采集数据并用GSAS-II软件进行Rietveld精修；用Zeiss EVO扫描电子显微镜（SEM）结合能谱（EDS）分析微观结构与元素分布；在布鲁克海文国家实验室国家同步光源II（NSLS-II）6-BM光束线上于10 K采集U和Gd的LIII边X射线吸收近边结构（XANES）与扩展X射线吸收精细结构谱（EXAFS）；采用激光闪射法测热扩散率、差示扫描量热法（DSC）测比热容量、热膨胀仪测线膨胀以计算热导率；利用共振超声谱（RUS）结合Rayleigh-Ritz反演提取弹性常数。

三、研究结果

3.1 相与晶格参数演变
XRD与Rietveld精修表明，所有Gd掺杂样品均仅指标定至萤石结构（Fm3?m），未检出Gd₂O₃多晶型或次生相衍射峰；随Gd含量增加，主衍射峰系统性向高2θ方向偏移，晶格参数线性减小，符合Vegard定律，证实Gd³⁺在U位取代，经验关系为a(?)=5.4693?0.00225x。

3.2 微观结构表征
SEM显示随Gd掺杂量增加晶粒尺寸减小；孔隙主要分布于晶界，6%和10% Gd样品中可见晶内小孔隙；与Hirai (1990)的报道不同，未观察到微裂纹；EDS未检测到晶内或晶界处的次生相，证实Gd均匀分布。

3.3 XANES与EXAFS结果
Gd LIII边XANES显示Gd保持+3价态，但Gd掺杂改变了其局域配位对称性（白线强度及边后特征增强）；U LIII边XANES白线强度随Gd含量增加而降低，15% Gd样品出现轻微蓝移，与缺陷团簇特征相符。EXAFS拟合表明Gd–O配位壳层分裂为2.26–2.43 ?多个距离，且Gd–O无序度大于U–O；Gd/U壳层解析出三个独立原子间距离（3.7 ?、3.80 ?、4.2 ?），表明萤石取代构型与类bixbyite局域构型共存；仅当Gd达15%时U–O壳层出现可分辨分裂；德拜-瓦洛因子随Gd含量系统性增大；EXAFS推导的晶格参数与XRD结果一致。

3.4 热物理性能演变
比热容量随Gd浓度升高而增大；热扩散率在低温段显著降低（15% Gd最低），归因于晶格缺陷的声子散射；热膨胀行为与UO₂总体一致，仅10%和15% Gd在高温段略有增大；热导率随Gd含量显著下降且低温更明显，采用λ(T)=1/(A+BT)模型拟合，A为杂质散射项、B为声子-声子Umklapp散射项。

3.5 共振超声谱弹性性能
未掺杂UO₂的弹性常数C₁₁和C₄₄与文献一致；随Gd₂O₃含量增加，C₁₁、C₄₄及推导的体积模量（bulk modulus, B）和杨氏模量（Young's modulus, E）均系统性下降，表明缺陷诱导的晶格软化；弹性模量降低意味着声子群速度下降，与增强的缺陷散射共同导致热导率退化。

四、讨论部分总结

讨论指出，Fm3?m（萤石）向Ia3?（bixbyite）的相变本质上是阴阳离子子晶格均参与的对称性破缺位移型转变：阳离子子晶格中U原子沿三轴发生微小位移，使Gd³⁺得以被容纳于局域Ia3?晶格；阴离子子晶格中氧原子由8c Wyckoff位置分裂为48e不等价位点，可优先与Gd³⁺配位。基于未掺杂UO₂晶格参数估算O/M比约为2.004，研究人员提出在还原性烧结条件下，电荷中性通过围绕U⁵⁺原子形成缺陷团簇实现，缺陷团簇呈立方体八面体（cuboctahedron）构型并产生一长一短两组U–O键距，这与EXAFS中氧壳层分裂及XANES白线强度降低相互印证。RUS弹性软化进一步支持异价取代与缺陷存在的判断。针对Hirai报道的微裂纹，研究人员认为其可能源于立方（C）→单斜（B）相变引起的体积收缩，而本工作未观察到单斜Gd₂O₃。EXAFS额外特征可借立方bixbyite结构配位环境描述，但强调这些类bixbyite构型局限于短程而不具长程一致性。与Cr掺杂促进晶界扩散和晶粒生长不同，Gd掺杂的缺陷主要局域于晶格内部。热阻分析显示杂质散射项随Gd浓度升高而增大，Umklapp散射项下降，表明Gd掺杂剂破坏了长程声子-声子相互作用，且结构紊乱可能降低总体声子电阻。最终，研究提出Gd掺杂UO₂应被描述为"整体萤石固溶体+局域短程非均匀性"的复合图像，这种局部结构二元性是缺陷浓度增加、弹性软化和声子散射增强的机制根源，将其纳入模拟框架有望提升热导率退化、弹性软化及长程缺陷演化的预测能力。

五、研究结论

研究得出以下结论：（1）Gd主要以萤石晶格形式被容纳，产生符合Vegard定律的系统性晶格收缩；（2）XANES和EXAFS揭示取代引入了复杂局域环境，包括氧空位、缺陷团簇和纳米尺度类bixbyite结构基序，这些局域畸变仅在短程尺度上关联，在介观距离上失去关联，从而保持了XRD所观测的平均萤石结构，同时改变了XAS所探测的局域原子环境；（3）SEM/EDS显示Gd均匀分布，无次生相或微裂纹，且Gd浓度升高伴随晶粒尺寸略微减小；（4）Gd掺入显著降低热扩散率与热导率，低温至中温区间尤为显著；（5）RUS揭示弹性刚度（C₁₁、C₄₄、B、E）随Gd含量增加而系统性降低，与缺陷诱导的晶格软化和缺陷形成相符；（6）热阻分析表明杂质散射占主导地位，声子-声子Umklapp散射随Gd含量升高而逐步受抑制；（7）所建立的结构-热物理关系为改进基于原子尺度信息的燃料性能模型提供了物理基础，凸显了同时表征平均晶格与局域原子环境的必要性。