《RSC Advances》：Pressure- and temperature-driven local structural stability and 5f/d hybridization in a metastable U–5.28 wt% Nb alloyOpen Access
编辑推荐：
亚稳态 U–Nb 合金在接近 6 wt% Nb 时，属于富 U 锕系金属材料，其中 Nb 取代、晶格畸变、振动稳定性以及 U-5f 电子退局域化之间存在强耦合。尽管高压实验和模拟已阐明 U–6Nb 及有序 U–Nb 化合物的关键特征，但对于稀溶性 U–5.28
亚稳态 U–Nb 合金在接近 6 wt% Nb 时，属于富 U 锕系金属材料，其中 Nb 取代、晶格畸变、振动稳定性以及 U-5f 电子退局域化之间存在强耦合。尽管高压实验和模拟已阐明 U–6Nb 及有序 U–Nb 化合物的关键特征，但对于稀溶性 U–5.28 wt% Nb 在压力和温度扰动下的原子级稳定性响应仍知之甚少。本研究将低焓结构搜索、声子谱、投影态密度（projected density of states, PDOS）以及从头算分子动力学（ab initio molecular dynamics, AIMD）导出的有限温度构型相结合，揭示压缩与加热如何重组局域稳定性。压力筛选结构在 50 GPa 至 300 GPa 区间呈现回滞的单斜-正交-单斜相变序列。在 50 GPa、100 GPa、250 GPa 和 300 GPa 条件下获得单斜晶胞，而正交晶胞出现在 150 GPa 和 200 GPa。这种非单调的对称性演化表明，中等压缩促进更有效的堆积，而更强的压缩则恢复局域单斜畸变的能量作用。声子计算支持这些压力筛选结构为动力学稳定的局域极小值，并显示整体压力诱导的模硬化效应。有限温度 AIMD 弛豫在 773 K 时得到低对称性构型，在 1023 K 时得到动力学不稳定构型，在 1273 K 时则得到更高对称性的动力学稳定构型。PDOS 分析显示，在压缩和加热条件下，U-5f 衍生态展宽以及 U-5f/U-6d/Nb-4d 杂化增强。这些结果为亚稳态 U–5.28 wt% Nb 在极端条件下的结构选择、声子响应与 5f/d 杂化之间的关联建立了原子级稳定性框架。
**论文解读：亚稳态U–5.28 wt% Nb合金的压力与温度驱动局域结构稳定性与5f/d杂化机制**
**一、研究背景与科学问题**
铀铌（U–Nb）合金是一类重要的富铀金属材料，Nb的添加能够稳定γ衍生及α″相关亚稳态、调控马氏体相变路径，并显著影响其力学与腐蚀性能。在5–6 wt% Nb成分范围内，合金可同时保留低对称性α″相关结构与γ衍生亚稳态，其分解、有序化及应变协调行为强烈依赖于热历史、冷却速率、微观组织和应力状态。这种亚稳特性赋予U–Nb合金独特性能组合，但也使其结构响应对外界压力和温度极为敏感。
压力与温度扰动U–Nb合金中不同但相互耦合的自由度：压缩减少体积、改变U–U和U–Nb配位并增强费米能级附近轨道重叠；加热增强原子振动、加速局域重排并改变静态畸变与高对称构型之间的平衡。已有实验和理论研究提供了重要参考：基于U–6Nb的原位高压高温中子衍射表明亚稳α″结构在显著压缩下仍可存续，而加热促进相变与分解行为；热力学与动力学评估显示U–Nb稳定性取决于成分相关的自由能和扩散系数；计算研究方面，Bao等人结合第一性原理与智能算法预测了0–200 GPa范围内U–Nb化合物性质，结构预测研究识别了时效U–Nb合金中的有序终态相和压力稳定的U–Nb金属间化合物，近期原子模拟揭示U–6Nb的γ → γ0 → α″马氏体相变通过局域原子改组与剪切而非简单晶体学转变完成，机器学习势模拟澄清了α-U基体的晶格动力学与高温稳定性。
然而，对于稀溶性U–5.28 wt% Nb这一工程合金成分，仍存在若干悬而未决的问题：其一，多数高压结构预测研究聚焦化学计量比化合物或有序合金，而工程合金为稀溶性替位固溶体；其二，有限温度结构响应常从纯U相区或后处理衍射结果推断，高温下合金构型的动力学稳定性缺乏直接解析；其三，压力诱导的U-5f态演化通常与结构稳定性分开讨论，而5f电子的局域-巡游倾向是锕系成键的核心问题。因此，需要建立一种统一描述，将局域结构、声子响应与费米能级附近电子态联系起来。
**二、研究内容与技术方法**
本研究以U–5.28 wt% Nb为模型亚稳富U合金，系统考察压力与温度如何重塑局域结构稳定性。在50–300 GPa压力范围内通过晶体结构搜索确定低焓结构，利用声子谱评估其动力学稳定性，分析费米能级附近U与Nb衍生态的演化；同时基于AIMD获得773 K、1023 K和1273 K的有限温度构型，并采用相同的结构-声子-电子框架进行检验。主要技术方法包括：采用CALYPSO粒子群优化算法进行晶体结构搜索（种群规模30、代数50、禁用对称性约束、最大晶胞维度3×3×3、重复结构筛选阈值1%）；基于Vienna Ab initio Simulation Package（VASP）的第一性原理计算，使用投影缀加平面波方法和Perdew–Burke–Ernzerhof广义梯度近似，平面波截断能520 eV，布里渊区采用Monkhorst–Pack网格（倒空间间距0.05 ??1），电子自洽收敛标准10?? eV/原子，结构优化力收敛标准0.02 eV ??1，所有计算均启用自旋极化；声子谱采用有限位移法计算；PDOS用于评估近费米能级各轨道贡献；AIMD在正则系综下进行，时间步长1 fs，总时长100 ps，前20 ps作为平衡段。研究未使用特定实验样本队列，全部基于第一性原理计算。
**三、主要研究结果**
**3.1 压力诱导的结构序列**
结构搜索在50 GPa、100 GPa、150 GPa、200 GPa、250 GPa和300 GPa六个代表性压力点均收敛至低焓构型。结构搜索揭示非单调的压力响应：50 GPa和100 GPa为单斜结构（分别属于Cm空间群和8原子Pm晶胞），150 GPa和200 GPa为更高对称性的正交结构（分别为Immm和Pmmm），250 GPa和300 GPa重新出现单斜畸变（分别为P2/m和Pm）。这一单斜-正交-单斜（monoclinic–orthorhombic–monoclinic, M–O–M）序列表明，理解U–5.28 wt% Nb的压力响应不能仅视为简单的密度效应，而是堆积效率与局域畸变能量竞争的结果。形成焓计算显示所有压力下ΔH
f均为正值（0.111–0.299 eV/原子），证实该稀溶体相对分解为纯α-U和bcc-Nb在热力学上是亚稳的，与已知的富U U–Nb合金亚稳特性一致。
**3.2 压缩条件下的动力学稳定性**
所有压力筛选结构沿高对称路径计算的声子谱均无显著虚频，证实这些低焓候选结构在相应压力下是动力学稳定的局域极小值。此外，声子频率整体随压力增加向高频偏移，表明压力导致原子间力常数硬化。这一结果与密度增加和费米能级附近更宽的轨道杂化程度一致，并将M–O–M序列的物理图像进一步明确：正交结构与单斜结构均在其压力区间内对应动力学稳定的局域极小值，而非瞬时结构。
**3.3 高压下的电子结构**
PDOS计算显示费米能级附近态密度主要由U-5f、U-6d、Nb-4d和Nb-5s贡献，且费米能级处有限态密度确认所有压力选择结构均具有金属特性。随压力增加，U-5f衍生峰展宽并与U-6d及Nb-4d态重叠增强，这是5f巡游性增强和U–Nb态耦合加强的关键电子特征。高压下同时出现的U衍生态和Nb衍生态展宽与更强的轨道重叠、更大的力常数和声子硬化一致，为M–O–M序列提供了结构-声子-电子关联的清晰物理解释。
**3.4 AIMD揭示的温度驱动结构响应**
AIMD模拟选取三个特征温度：773 K对应纯铀α相稳定区、1023 K位于α→β转变区附近、1273 K对应γ相高温稳定区。弛豫后结构显示：773 K构型保持低对称单斜晶胞（C2/m）；1023 K构型演化为四方结构（P4?m2），但该构型声子谱存在显著虚频，揭示出动力学不稳定的变换敏感区；1273 K构型弛豫为高对称立方结构（P2/m3?），其中Nb以单个替位缺陷形式存在。密度从773 K的17.89 g cm?3降至1273 K的16.87 g cm?3，与热膨胀和静态畸变成分的减少一致。MSD和RDF分析表明1273 K合金在模拟时间尺度内保持晶态固体特征，无液态扩散行为。这一温度序列识别出合金从低对称局域极小值出发，经由软模变换敏感区间，最终到达高对称动力学稳定构型的原子级演化路径。
**3.5 高温电子结构**
AIMD导出构型的PDOS计算显示，773 K和1023 K时U-5f/U-6d与Nb-4d/Nb-5s在费米能级附近仍为主要贡献者；1273 K高对称构型表现出更强的U–Nb态杂化。高温电子响应可描述为热膨胀、原子重排与电子再杂化的耦合：高温降低低对称畸变的持续性，而费米能级附近U和Nb态继续定义金属键合网络。声子与PDOS联合结果表明，1273 K弛豫构型的稳定性同时来源于更高对称性和更强的U–Nb电子混合。
**四、讨论与结论**
**4.1 压力诱导的局域稳定极小值重排**
核心结果是U–5.28 wt% Nb中压力依赖的低焓动力学稳定结构候选的回滞序列。压缩通过改变配位、单斜畸变和U–Nb电子杂化来改变局域极小值的层级结构，而非单纯驱动合金向更高对称性发展。证据层级在三个描述符间内部自洽：焓排序识别候选结构极小值、声子谱确认动力学稳定性、PDOS趋势将结构序列与压力诱导的5f/d展宽和U–Nb杂化关联。所有压力构型的正形成焓进一步澄清了稀溶体的亚稳程度，为M–O–M对称性演化趋势提供热力学基础。
**4.2 亚稳U–Nb中温度辅助的对称性演化**
温度计算将AIMD导出结构置于U–Nb合金亚稳性图景中。773 K与1273 K衍生构型弛豫至动力学稳定的局域极小值，而1023 K构型动力学不稳定，识别出一个变换敏感窗口而非平滑的稳定中间结构序列。这一模式与晶格畸变、成分再分布与局域应变耦合控制的相演化一致。温度依赖结构应理解为原子级稳定性指示器，长时间扩散、化学分解、沉淀与介观尺度应力协调是自然的后续拓展方向。
**4.3 压力与温度驱动稳定性的电子起源**
PDOS结果提供了连接压力与温度观测的电子线索。压力下U-5f衍生态展宽并与U-6d和Nb-4d态重叠增强，将M–O–M结构序列与增强的金属键合和5f巡游性联系起来；高温下1273 K高对称构型也显示更强的U–Nb重叠，表明电子再杂化同时伴随压缩驱动和温度辅助的稳定化过程。整体物理图景为：U–5.28 wt% Nb通过高效堆积、低对称应变弛豫和U-5f/Nb-4d杂化之间的竞争选择局域结构。
**五、研究结论**
本研究建立了亚稳态U–5.28 wt% Nb中压力与温度驱动稳定性的原子级结构-声子-电子框架，主要结论如下：(1) 50 GPa至300 GPa压力筛选低焓结构呈现非单调的单斜-正交-单斜序列，表明压缩改变局域配位和晶格畸变而非简单推动合金向更高对称性演化；(2) 声子计算显示压力筛选结构无显著虚频，证实其局域动力学稳定性，压力诱导的声子模硬化与密度增加、力常数增强和近费米轨道重叠增强一致；(3) AIMD导出构型在773 K呈低对称结构，1273 K为高对称动力学稳定构型，1023 K构型存在显著虚声子模，识别出两弛豫态之间的变换敏感动力学不稳定性；(4) PDOS分析表明U-5f/U-6d和Nb-4d/Nb-5s态主导费米能级附近金属键合，压力展宽U-5f衍生态并增强U–Nb杂化，1273 K弛豫构型亦显示更强的U–Nb电子混合，增强的5f/d杂化是压力诱导序列与温度辅助对称性响应的关键电子因素。该工作通过聚焦局域畸变与U-5f/Nb-4d杂化的耦合作用，区分了稀溶性富U固溶体行为与此前研究的U–6Nb实验、有序U–Nb相及化学计量比金属间化合物，为理解亚稳锕系合金在极端条件下的稳定性提供了统一框架。
