【摘要】 Stoner 判据 I·D(EF)>1 怎么用于 SAA 磁性初筛？本文说明 DOS 筛查、自旋 DFT、多初始磁矩和临界体系复核方法。

科研人必备干货｜全套科研知识库

常用科研资料汇总，点击领取即可获取。

 

Stoner 判据 I·D(EF) > 1 可以作为单原子合金磁性倾向的初步筛查工具，其中 D(EF) 表示费米能级附近的态密度。它适合帮助决定哪些体系应优先做自旋 DFT，但不能替代最终的自旋计算和多初始磁矩复核。

 

为什么要先看磁性倾向

对 Cu、Ag、Au 等宿主表面中的过渡金属单原子，局域电子结构可能产生显著磁矩。若直接沿用非自旋设置，可能在优化、吸附能和反应能学计算开始时就选错电子基态。

磁性筛查的目标不是给体系贴上绝对标签，而是识别需要进一步核验的候选体系。

 

四步初筛流程

1. 运行非自旋 DOS

先对结构进行非自旋计算，检查费米能级附近的态密度，并记录体系、参数和收敛状态。DOS 图的能量零点和展宽设置应保持一致，便于不同体系比较。

2. 应用 Stoner 判据

对候选体系估算 I·D(EF)。当结果明显大于 1 时，可将其列入优先自旋计算清单；当结果接近临界值时，不宜仅凭判据下结论。

页面报告该判据在所讨论分类任务中的准确率约为 90%。这个比例属于特定研究范围的结果，不应直接当作所有 SAA 体系的普适准确率。

3. 设置多个初始磁矩

对候选结构使用多组初始磁矩计算，比较最终磁矩、总能和收敛状态。不同初值收敛到不同磁态时，应进一步判断哪一个状态更稳定、计算是否充分收敛。

4. 做 spin/non-spin 对照

在结构和参数可比的前提下，比较自旋与非自旋总能。页面讨论的部分体系能差接近 2 eV，说明默认关自旋可能造成明显基态偏差。

 

磁性与吸附能的联系

自旋极化会改变局域轨道占据、d 带中心、杂化和电荷重排，因此可能改变吸附能。页面给出的代表性对照中，H 吸附能变化最大约 0.33 eV，CO 体系最大约 1.55 eV。

分析吸附能时要同时记录吸附位点、参考态、结构是否重新优化、磁态和能量符号。不能把特定体系的差值直接作为所有材料的统一修正。

 

临界体系的复核清单

  • 初始磁矩是否覆盖合理范围；

  • 自旋与非自旋计算是否使用可比结构和参数；

  • 电子与离子收敛是否达到目标；

  • slab 层数、偶极修正和 SOC 是否影响结论；

  • 磁矩、总能、DOS 和吸附能是否相互对应。

科学指南针·唯理计算页面展示专业 DFT 计算服务和计算方案咨询。需要讨论自旋初筛时，可准备 DOS、元素组合、表面模型和当前计算设置。

论文信息：Zhang, S.; Riviere, C. I.; Montemore, M. M. “Magnetic Moments in Single-Atom Alloys: Trends, Mechanisms, and Implications for Tunable Adsorption.” ACS Physical Chemistry Au (2026). DOI: 10.1021/acsphyschemau.6c00086。

 

常见问题

I·D(EF)>1 就一定有磁性吗？

它更适合作为磁性倾向的初筛线索，不能替代自旋 DFT。临界值附近的体系尤其需要多初始磁矩、总能和最终磁矩的联合判断。

初始磁矩设置一个值够吗？

对可能存在多磁态的体系，一个初始值可能不足以找到更稳定的收敛状态。建议根据元素和局域位点设置多组初始磁矩，并比较最终状态。

自旋计算为什么会改变吸附能？

自旋会改变电子占据、d 带中心、吸附质与表面的杂化和电荷重排。最终变化大小取决于体系、吸附位点、磁态和结构优化条件，不能用单一固定数值概括。

配套阅读：SAA 自旋计算与基态判断 和 单原子合金 DFT 计算方案。