说明:本文华算科技主要介绍高熵效应、晶格畸变效应、迟滞扩散效应和鸡尾酒效应各自对应的计算对象,并讨论四类判断成立所需的相态、温度、局域结构与迁移机制条件。



高熵材料通常含有多个占比较高的组元。理想随机固溶体的构型熵写作 ΔSconf = −RΣciln ci;等原子比、N 个组元时为 RlnN。这个数值描述理想占位组合的统计规模,它通过自由能中的 −TΔSconf 项参与相竞争。有限温度下的实际构型熵还受短程有序、偏聚和相分离约束,组元数量本身只给出理想上限。
“高熵四大效应”把多组元体系里常见的相稳定、局域位移、原子迁移和性能协同概括为四个名称。用于计算时,四个名称应分别改写成可比较的量:相自由能与有序参数对应高熵效应,键长或原子位移分布对应晶格畸变,示踪扩散系数与迁移能垒对应迟滞扩散,局域位点的电子结构和反应能分布对应鸡尾酒效应。每一项都带有参照态,脱离参照态便无法判定效应幅度。
图1中的占位构型来自有限温度蒙特卡洛采样。同一名义成分随温度降低依次出现部分 B2 有序和 BCC/B2 相混合,说明构型数多与保持单相属于两项不同判断。若异种原子配对显著降低焓,化学有序或分相带来的能量收益可以抵消熵项。高熵效应的证据应来自不同候选相的自由能或有序转变温度,而非一张等原子比配方表。
第一性原理计算给出构型能、局域松弛能和电子态;簇展开或机器学习势把这些能量扩展到大构型空间;蒙特卡洛追踪化学有序与相占比;分子动力学处理热振动和扩散。计算方法与效应名称之间没有一一绑定,同一效应常由两类以上可观测量约束。例如晶格位移统计依赖结构弛豫,温度响应还涉及动力学采样。
四类效应之间存在物理耦合。局域化学有序会改写键长、缺陷能和电子态,进而改变位移分布、迁移路径与吸附位点。计算报告若把四项并列成材料标签,往往会遗漏共同的微观变量。较明确的表述应指出材料处于何种相、采用何种温度或磁态、比较哪一个参照体系,以及观测量变化有多大。



相稳定的直接对象是吉布斯自由能 G = H − TS。构型熵随组元数增加,混合焓、振动熵、磁熵和电子自由能同样随成分与相结构变化。高熵效应成立的范围由自由能交叉温度限定;某个固溶相在高温占优,冷却后仍可能进入有序相或多相区。零温 DFT 构型能可筛选焓驱动趋势,却不直接给出有限温度的相比例。
元素分辨的 B2 有序参数在高温接近零,降温后逐渐偏离零;增加 Al 还会推高有序起始温度。有序参数提供了比“单相或多相”更细的构型证据,因为宏观晶格仍可保持 BCC,而两个子晶格的占位偏好已经发生变化。此时把理想 RlnN 全部计入自由能会高估无序程度,构型熵应由采样到的占位概率或热力学积分修正。
相竞争计算还受局域松弛控制。不同原子半径和键合偏好使原子离开理想格点,弛豫后的构型能排序可能改变,随之移动有序转变和相分解温度。固定理想格点的哈密顿量覆盖化学占位,弛豫模型才包含位移自由度。两类模型给出差异时,差值本身就是晶格松弛参与相稳定的量化结果。
热容峰对应构型能涨落增强,代表体系发生有序转变或结构重排。峰位随模型是否允许局域松弛而移动,表明熵稳定温区必须与同一能量模型下的焓项配套。有限超胞会限制长周期有序和相分离尺度,计算温区还应与合成、退火和测量温度对应;否则,高温随机固溶体的判断会错误外推至室温状态。
相稳定判断可写成明确的证据组合:候选相自由能差、短程有序参数、热容异常和相占比随温度的变化。只有构型熵公式时,结论限于理想随机模型,无法覆盖有序、偏聚、析出或非等原子比体系。材料名称中的“高熵”描述成分设计类别,热力学高熵效应则对应特定温度下可计算的自由能贡献。



多种元素共占一个晶格时,原子半径、价电子和键强差异会产生分布化的局域平衡位置。结构优化后,单个原子相对理想格点的位移、最近邻键长分布和局域应变张量都可描述静态畸变。晶格常数的平均值不承载局域位移分布,所以平均晶胞看似规则时,原子尺度仍可能存在宽广的键长与应力分布。
常用量之一是均方原子位移 MSAD = N−1Σ|ri − ri,0|2,其平方根与长度同量纲。参照位置 ri,0 可取相同晶格常数下的理想格点。参照晶格、超胞形状和是否放松体积会直接改变 MSAD,跨论文比较时应核对定义一致。零温优化所得为静态位移,有限温度轨迹还包含热振动,两者可通过时间平均位置或低温外推分开。
MSAD1/2 把大量原子位移压缩为一个尺度,适合比较不同成分和计算方法。单一均值仍会掩盖元素分辨差异,完整判断还应读取各元素位移概率、键长分布尾部及局域应力。较大的 MSAD 指向较宽的离位尺度,却不单独指定畸变来自尺寸失配、化学键偏好还是局域有序。弛豫能 ΔErelax 可补充衡量原子离开理想格点后释放的能量。
允许局域松弛后,部分原子对的短程有序参数明显偏离未松弛结果,说明位移场会反馈到化学占位偏好。尺寸差较大的原子对若能通过局部位移降低应变能,其配对概率可能增加;某些化学上有利的配对还可产生特定键长。把畸变效应写成“原子种类多导致晶格严重扭曲”会丢失这种双向作用。
计算上的证据强度可按信息量排列:弛豫前后能量差确认位移自由度降低了能量,元素分辨位移与键长分布确定畸变集中在哪些局域环境,多个独立构型给出统计波动。单个 SQS 构型只代表一种有限尺寸近似,其 MSAD 若随构型种子变化显著,平均值应附带离散范围,避免把偶然局域排布写成材料固有尺度。



扩散速率由缺陷形成、局域迁移能垒和迁移路径的空间连通程度共同控制。多组元随机环境会把单一迁移能垒展开成分布,空位、间隙原子或交换机制还会采样不同的元素组合。迟滞扩散是给定温度与机制下的相对速率判断,参照物可以是纯金属、低组元合金或相同晶格的成分序列。缺少参照扩散系数时,“扩散迟滞”只剩定性命名。
示踪扩散系数 Di 从元素 i 的长时间均方位移斜率得到,也可由缺陷浓度、跃迁频率和相关因子组成。Arrhenius 斜率给出表观活化能,它混合了主导缺陷和路径随温度的变化。某一元素的 Di 降低时,其余组元的 Dj 仍可上升或保持不变,多组元体系常出现分量扩散系数分叉,最快迁移组元还可能改变缺陷哑铃或空位邻域的组成。
Ni–Fe 连续成分结果显示,扩散系数随 Fe 含量呈非单调变化,等原子比位置并无普遍的最小值。迁移哑铃的元素组成同步变化,反映局域缺陷态选择了不同路径。化学无序度与扩散速率之间缺少固定单调关系,因为提高组分复杂度既可能引入深能量陷阱,也可能增加低能迁移路径。具体结果由能量分布和路径连通程度决定。
两类局域缺陷态即使只相差形成能或迁移能中的一项,便会产生不同的成分—扩散曲线。低能位点数量达到连通阈值后,长程扩散可转入新的渗流网络;高能垒位点若绕行成本较低,对宏观 D 的抑制有限。迁移能垒直方图没有记录路径连通关系,动力学蒙特卡洛、长时间分子动力学或显式跳跃网络可补足空间相关信息。
扩散判断应报告晶相、缺陷类型、成分、温区和比较基准。低温下某类陷阱主导的迟滞,在高温跨越多条路径后可能减弱;辐照产生的间隙扩散与热平衡空位扩散也对应不同能量网络。“高熵合金扩散普遍缓慢”超出了单一机制计算能够支持的范围,更具体的结论应写成某元素、某缺陷机制和某温区内 D 或表观活化能的相对变化。



鸡尾酒效应用来描述多种元素共同存在后出现的非线性性能响应。计算中没有名为“鸡尾酒”的独立标量,证据通常来自局域配位、应变、电荷转移和轨道杂化共同造成的性质分布。可预测的协同应能追溯到具体位点或相邻位点组合,并与单元素、低组元或打乱局域结构的参照模型比较。只列出各元素功能会把协同作用误写成成分清单。
在电催化高熵合金表面,同一晶面含有大量邻域组合。局域电负性、d 带占据和应变会改变 H*、OH* 等中间体的吸附自由能,水解离与氢脱附因而可能由不同金属位点承担。位点分工的计算对象是反应步骤分辨的能垒和吸附自由能,而非把某个平均 d 带中心同时解释全部反应步骤。覆盖度升高后,吸附物之间的相互作用还会重新排列优势位点。
高熵合金的多位点密度提供了连续调节吸附能的空间,也扩大了不利位点的数量。真正的协同表现为某些邻域降低一个基元步骤的自由能代价,另一些邻域接续处理后续中间体,且整体反应路径的最大上坡自由能减小。位点能量分布的宽度只描述异质程度,活性还取决于分布中心、尾部占比、位点连通和反应物到达概率。
在 FeCoNiMnRu 案例中,计算把水解离与 H* 吸附分别投影到不同金属位点:Co 位点的水解离能垒为 0.34 eV,Ru 位点的 ΔGH* 接近热中性的 −0.07 eV。两项能量对应串联反应中的不同任务,原位拉曼中的金属—氧和 Ru—H 信号则提供工作电位下的中间体线索。计算负责位点归属与能量比较,光谱负责检验界面上是否出现相应吸附物。
这类位点归属仍受表面偏析、溶剂、电势和吸附覆盖度影响。真空零覆盖度模型给出的优势位点,在工作条件下可能因电极双电层或邻位吸附而重排。鸡尾酒效应的判断范围应与所建表面构型和反应环境一致,动态重构后出现的新配位环境应重新计算能垒,并用原位信号确认活性位点是否仍然存在。
