vasp计算高熵合金FeCoCrNi

一、项目背景与科学问题

高熵合金以多种主元素近等摩尔比例合金化为设计理念，打破了传统合金”一种主元加微量合金化”的框架，展现出高强度、优异的热稳定性与抗辐照性能。FeCoCrNi是面心立方结构高熵合金的模型体系，其单相固溶体是否稳定、各元素混合焓的正负与大小如何，直接关系到合金能否在服役温度下保持单相组织而不发生偏聚或析出。本项目运用第一性原理VASP软件计算FeCoCrNi及其子三元系的混合焓，并结合构型熵评估其相稳定性，从热力学根源上回答高熵合金”熵稳定”的微观机制。

二、计算方法与理论背景

高熵合金相稳定性的判据是Gibbs自由能G=H-TS。混合焓H_mix反映不同原子间的化学亲和力：若异种原子结合更牢固，H_mix为负，倾向于形成有序相或金属间化合物；若同种原子更易聚集，H_mix为正，倾向发生偏聚。单相固溶体得以存在的关键在于高温下构型熵项-TS_conf足够大，能够补偿正的混合焓。对等摩尔n元合金，构型熵为S_conf=kB·ln(n)，FeCoCrNi四元等摩尔体系的S_conf=kB·ln4，在1000 K时-TS_conf约为-0.119 eV每原子。本项目用特殊准随机结构（SQS）在有限超胞内近似无规原子分布，再以DFT精确计算体系总能，从而得到混合焓。

三、计算设置与模型构建

构建了FeCoCr、FeCoNi、FeCrNi、CoCrNi四个三元子体系与FeCoCrNi四元体系的SQS超胞，每超胞含32至108个原子以保证随机性足够好。VASP计算采用PAW方法描述离子-电子相互作用，PBE交换关联泛函，平面波截断能取400 eV。考虑Fe、Co、Ni、Cr均为磁性元素，所有计算均开启自旋极化（ISPIN=2），初始终磁矩按铁磁构型设置。倒空间k点采样按标准精度设置，结构优化收敛后计算总能，混合焓定义为合金总能减去各纯元素参考总能之差，再按原子数平均。

四、结果分析与案例图解读

图2的左侧面板（a）给出各体系混合焓的柱状图，单位为meV每原子。可见FeCoCr、FeCoNi、FeCrNi的混合焓为正，分别为21、12与33 meV每原子，说明这些三元组合中同种原子有一定聚集倾向；CoCrNi则呈现约-15 meV每原子的负混合焓，表明Co与Cr、Ni之间亲和力较强。四元FeCoCrNi的混合焓为6 meV每原子，绝对值明显小于各子体系的最大值，反映出多元合金化对局部原子应力的平均化效应。图中紫色虚线标出1000 K时构型熵贡献-TS_conf=-119 meV每原子，其绝对值约为混合焓的二十倍，说明在该温度下构型熵完全可以压制正的混合焓，从而稳定单相面心立方固溶体——这正是高熵合金”熵稳定”的直接数值证据。

图2的右侧面板（b）给出FeCoCrNi自旋分辨的态密度。多数自旋通道与少数自旋通道的态密度在费米能级附近出现明显的不对称分布，说明合金存在显著的磁矩贡献，Fe与Co的局域磁矩使费米能级附近态密度呈现分裂特征。d电子带在费米能级附近保持较高的态密度，表明该合金具有良好的金属导电性；同时各元素的d轨道相互重叠形成宽广的d带，有利于固溶强化与良好的塑性。自旋分辨态密度还显示少数自旋态在费米面附近占据较高，为后续讨论磁性对相稳定性的影响提供了线索。

两面板结合，左面板从能量角度量化了混合焓与熵的竞争，右面板揭示了电子结构层面的成键与磁性特征，共同完成FeCoCrNi高熵合金相稳定性的第一性原理评估。

五、结论与展望

本项目用VASP结合SQS模型计算了FeCoCrNi及其子体系的混合焓，四元系混合焓仅6 meV每原子，1000 K时-TS_conf达-119 meV每原子，构型熵占主导，单相面心立方结构稳定。案例图以混合焓柱状图与自旋态密度双面板呈现。后续可拓展：(1)对更复杂五元、六元体系开展SQS高通量计算，筛选新高熵合金；(2)引入HSE杂化泛函考察磁性对电子结构的修正；(3)结合声子计算评估晶格振动熵与有限温度自由能。