氧取代诱导的UC1?xOx结构、力学及电子演化：面向TRISO UCO燃料应用的第一性原理研究

《Annals of Nuclear Energy》：Oxygen substitution-induced structural, mechanical, and electronic evolution in UC1?xOx: A first-principles study for TRISO UCO fuel applications

　　

采用第一性原理方法，系统研究了氧取代对无序铀碳氧化物固溶体UC1?xOx在低氧含量范围内的结构、热力学、电子、磁性与力学性能的影响。研究将密度泛函理论（DFT）中的GGA+U框架与特殊准随机结构（SQS）方法相结合，以模拟阴离子亚晶格上碳/氧的取代无序，氧含量分别为0、4、8、12和16 at.%。同时采用占据矩阵控制（OMC）方案，确保所有组分下铀5f电子态的一致收敛。计算结果表明，在研究的全部氧含量范围内，UC1?xOx均保持岩盐（NaCl型）结构，未发生对称性降低。晶格参数呈现非单调变化：在4 at.%氧含量时出现晶格收缩，归因于较短且较强U–O键的形成；而在8–16 at.%范围内晶格参数单调增大，反映了累积晶格应变与电荷再分布效应。形成焓随氧含量增加而单调降低，表明氧取代在热力学上稳定碳化物晶格，且U–O键的离子性增强。电子结构分析显示，随着氧含量增加，铀磁矩从近乎零（本征UC）单调增大至约0.4–0.5 μB，费米能级处态密度逐渐降低，表明铀5f电子局域化增强及U–O杂化加强。在最高氧含量下出现极小的带隙特征，被解释为电子性格的交叉而非明确的金属-绝缘体转变。弹性性质计算表明，氧取代使剪切常数C44及剪切模量G显著增大，而纵向弹性常数C11基本保持不变，体模量B仅略有变化，泊松比ν随氧含量增加而降低。结合晶体轨道Hamilton布居（COHP）分析，研究人员确认U–O键的键合强度显著高于U–C键，且随氧含量增加U–O键贡献持续增强。这些结果表明，氧取代增强了碳化物晶格的剪切刚度和键方向性，同时保持了体积稳定性。该研究建立了低氧含量区间内无序UC1?xOx组分-结构-性能之间的定量关系，为TRISO UCO燃料在辐照条件下碳/氧再分布所导致的碳化物相演化提供了原子尺度的参考基准，对燃料性能预测和辐照行为建模具有重要参考价值。

论文解读

一、研究背景与科学问题

先进核能系统对燃料材料提出了日益严苛的要求，包括高燃耗、高温、长辐照寿命及瞬态事故工况下的可靠性。高温气冷堆（HTGR）中，三层各向同性（TRISO）包覆颗粒燃料因其优异的裂变产物包容能力和结构稳定性，已成为先进反应堆的主要燃料概念。TRISO燃料核心由铀氧碳化物（UCO）组成，这是一种碳化物-氧化物混合体系，旨在辐照过程中维持受控的氧势。与纯UO₂燃料芯核相比，UCO组分同时包含碳化铀和二氧化铀相，可有效缓冲氧势、抑制有害气体生成并增强与SiC层的相容性。

在辐照前，UCO芯核通常呈现UO₂、α-UC₂和UC多相共存的非均匀微观结构。随着燃耗推进，α-UC₂相因热力学不稳定而逐渐分解，剩余的碳化物和氧化物区域向非化学计量固溶体转变。氧和碳在各自宿主晶格中的不对称溶解度导致氧取代碳化物（UC_1?xO_x）和碳取代氧化物（UO_2?xC_x）相的形成。热力学模型和CALPHAD评估预测，UC中氧溶解度在辐照相关温度下可达约12–16 at.%，而UO₂中的碳溶解度相对有限。这种不对称性驱动碳化物富集区域在燃耗过程中逐渐富氧，导致碳化物相发生组分演化。

然而，尽管UC_1?xO_x在辐照演化的UCO芯核中扮演核心角色，其原子尺度的本征性质仍缺乏系统表征。现有原子尺度研究多集中于化学计量UC或孤立缺陷，而非随组分变化的固溶体。此外，许多计算研究采用有序超胞和周期性氧排列，无法准确捕捉取代氧掺入的固有统计无序性。铀5f电子部分填充且位于费米能级附近，在锕系碳化物中介导键合相互作用，使UC处于巡游金属与局域关联电子行为之间的中间状态。标准广义梯度近似（GGA）可能低估在位库仑相互作用并使5f态部分退局域化。DFT+U框架中的在位库仑修正可增强电子局域化和结构性质的描述，但缺乏系统评估氧取代如何影响辐照相关组分下5f态演化的研究。

取代无序与关联5f电子行为之间的相互作用进一步增加了UC_1?xO_x原子尺度建模的复杂性。氧掺入降低局域对称性并创造化学上不同的U–C和U–O环境，增加了多种低能电子构型存在的可能性。因此，在统一的DFT+U框架内，同时解决结构、热力学、电子和力学性质的综合组分依赖第一性原理研究仍属空白。低氧含量下与TRISO芯核演化最相关的氧浓度-晶格畸变-5f电子再分布-键合特征-弹性响应之间的定量关系尚未建立。

二、研究开展内容

本研究通过系统性的第一性原理研究填补了这一空白，采用DFT+U结合特殊准随机结构（SQS）方法模拟阴离子亚晶格上的取代C/O无序。考虑0、4、8、12和16 at.%的氧浓度，对应TRISO UCO燃料芯核中辐照驱动氧再分布相关的低氧范围。研究分析了晶格参数、形成焓、铀5f电子态、磁矩和弹性性质，以阐明氧掺入如何改变碳化物富集的UC晶格。该工作的新颖性在于建立了低氧区间内无序UC_1?xO_x的组分依赖原子结构-性能关系，为未来缺陷形成、裂变产物掺入和辐照诱导组分演化研究提供了基线参考。

三、关键技术方法

研究采用维也纳第一性原理模拟软件包（VASP 6.5.0）中的投影缀加平面波（PAW）方法进行密度泛函理论（DFT）计算。交换关联能采用Perdew-Burke-Ernzerhof（PBE）形式的广义梯度近似（GGA）处理。铀的局域库仑相互作用采用Dudarev形式主义中的有效Hubbard参数U_eff=3 eV描述。构建2×2×2面心立方NaCl型UC超胞，共64个原子（32个铀原子和32个阴离子位点）。氧原子取代碳亚晶格上的碳原子，对应实际氧含量为0、4.69、7.81、12.50和15.63 at.%，按名义浓度报告为0、4、8、12和16 at.%。特殊准随机结构（SQS）方法通过Alloy Theoretic Automated Toolkit（ATAT）生成，使相关函数尽可能接近理想无限随机合金。占据矩阵控制（OMC）方案用于所有含氧组分，初始占据矩阵采用单轨道占据模式（0100000），在自洽迭代中约束前30个电子步以稳定局域化5f构型。弹性常数通过应力-应变方法计算，体模量、剪切模量和杨氏模量分别通过Voigt-Reuss-Hill近似推导。对12 at.%氧含量组分进行了64原子与216原子超胞的收敛性测试，确认64原子SQS模型在结构、热力学和磁性趋势分析中已充分收敛。

四、研究结果

4.1 结构、体积与热力学响应

SQS模拟表明，UC_1?xO_x在整个研究组分范围内保持岩盐结构，直至16 at.%氧含量均无对称性降低的证据。晶格参数随氧含量呈显著非单调变化：4 at.%氧含量时晶格参数降至4.981 ?，相对于本征UC（5.032 ?）发生收缩，归因于较短且较强的U–O键取代U–C键带来的局域键合增强。在8–16 at.%氧含量范围内，晶格参数单调增大，16 at.%时达到5.093 ?，表明从局域键缩短主导转向长程晶格相互作用主导的机制转变。键长统计证实，4 at.%时U–C和U–O键长分别降至2.496和2.493 ?，16 at.%时增至2.547和2.543 ?。Bader电荷分析显示，氧亚晶格上的总电荷积累从4 at.%时的?2.780 e增加至16 at.%时的?14.620 e，证实U–O键的离子贡献随氧掺入增强。形成焓计算以α-铀金属、石墨碳和分子氧（O₂）为参考态，本征UC形成焓为?0.86 eV/atom，随氧含量增加单调降低至16 at.%时的?1.13 eV/atom，表明氧取代在热力学上稳定碳化物晶格。这一趋势与UO₂中碳掺入导致晶格参数单调增大的行为形成鲜明对比，反映了碳化物和氧化物主导燃料中根本不同的键合环境。

4.2 电子与磁性质

本征UC的铀磁矩接近0 μB，符合Pauli顺磁体特征，表明5f电子以巡游性为主。Hill判据表明UC的U–U间距接近3.4–3.6 ?的临界阈值，使该化合物处于磁不稳定性边缘。随氧掺入，铀磁矩单调增大，最高氧含量时达到约0.4–0.5 μB，反映U–O键较U–C键更强所诱导的5f局域化增强。自旋轨道耦合（SOC）效应对UC电子结构的影响主要体现在5f流形的精细结构上，使未占据5f态分裂为j=5/2和j=7/2分量，但不会诱导金属-绝缘体转变。随氧浓度增加，O-2p态在约?6至?3 eV能量范围内出现并与U-5f态杂化，费米能级处态密度单调降低。最高氧含量时出现极小的带隙特征，但因带隙极小且对k点采样、展宽宽度和有限超胞尺寸敏感，被解释为电子性格交叉而非真正的金属-绝缘体转变。

4.3 力学性质

弹性常数计算显示，Hubbard U参数对纵向弹性常数C₁₁影响极小（311–318 GPa，约2%变化），但对剪切常数C₄₄影响显著，从U=1 eV时的36.6 GPa增大至U=3 eV时的69.9 GPa，增加近90%。体模量基本不变（189–194 GPa），而剪切模量和杨氏模量分别增大47%和40%。SOC的加入使体模量从约193 GPa降至约152 GPa，接近实验值范围（156–158 GPa），但不会改变U驱动剪切硬化的定性趋势。由于包含SOC的全弹性张量计算在大型SQS超胞中计算成本高昂，且本研究主要目标是分析氧掺入诱导的相对力学趋势，后续UC_1?xO_x弹性计算在非SOC DFT+U框架下进行。氧掺入对弹性张量的影响表现为：C₁₁基本不变，C₁₂逐渐减小，C₄₄持续增强，剪切模量相应增大，体模量仅略有变化，泊松比随氧浓度增加而降低。COHP分析表明，U–O键的ICOHP值（4 at.%时为?3.510 eV/bond）比U–C键（?2.220 eV/bond）更负，且随氧含量增加U–O键数从12增至60，加权U–O贡献从?42.120增至?237.000 eV，证实氧掺入引入强U–O键合相互作用，同时削弱和减少U–C键合网络。

五、结论与讨论

本研究建立了低氧含量区间内无序UC_1?xO_x的定量组分-结构-性能关系。氧取代在热力学上稳定碳化物晶格，增强铀5f态局域化，提高剪切刚度并增加键方向性，同时保持纵向刚度和体积稳定性。从工程角度看，这些原子尺度趋势与辐照TRISO UCO燃料芯核中碳化物富集区域直接相关：燃耗过程中碳和氧的再分布可产生氧取代的UC_1?xO_x区域。岩盐框架的保持表明碳化物相在低氧区间具有结构鲁棒性，而非单调的晶格响应表明氧再分布可能影响局域体积应变和尺寸稳定性。热力学稳定性的增强和U–O键贡献的增加表明氧掺入可改变缺陷容纳能力和裂变产物相互作用行为。剪切相关刚度的增加和泊松比的降低表明富氧UC区域可能展现增强的键方向性和改变的局域应力响应。该研究建立的定量关系为UCO TRISO燃料芯核的缺陷形成、裂变产物掺入和辐照诱导微观结构演化研究提供了有用的原子尺度参考。论文发表在《Annals of Nuclear Energy》。