ATAT钴自由能：统一参照、合并与作图

本篇完成同一Co实例的最后一步：所有项使用每原子单位，静态能量只加入一次。

1 先明确每一列是什么#

cd ATAT-Co
cat reference/model-summary.json
head -n 6 reference/Co-common-free-energy.csv

本例静态参照选同一固定几何两原子FM晶胞的E0/2，即−7.0189525 eV/原子。振动项保留零点能；电子项是原生felec修正；磁性项以同一FM态为零点。

F_total(T) = E_static_FM/atom
           + F_vib(T)/atom
           + DeltaF_elec(T)/atom
           + DeltaF_mag(T)/atom

不能把AFM的完整总能量再加进来，也不能在已经包含振动零点能的曲线上再加一次零点能。

2 原生mkteci怎样接文件#

cd reference/combined
cat Trange.in
cat eci.out
head -n 8 fvib.eci
mkteci fvib.eci felec.eci fmag.eci
cat teci.out
cd ../..

本例只用一个常数能量系数来演示组合，不是假设已经拟合出任意合金的温度簇展开。Trange.in为：

1000 11

表示到1000 K，共11个温度点；本例对应0、100、…、1000 K。eci.out放一项静态参照，各修正文件按温度块依次排列，原生mkteci输出总和。

3 一条短命令重新走完已有数据后处理#

bash postprocess.sh
cat work/Co-common-free-energy.csv

在包根目录运行。脚本将已提供的输入复制到新建work，依次运行fitsvsl → svsl → felec → icamag → mkteci，最终得到CSV。已有work时换一个解压目录，保留原来的结果。

这条命令不启动VASP。若通过submit.sh生成自己的新DFT结果，应先把新位移力、对应结构、自旋总DOS和两种磁态能量更新为同一套数据，再执行后处理；不要把fresh的新能量与reference中的旧力和DOS混合。

4 作图与实际结果#

bash plot.sh
ls plots

绘图读取reference的本站实算数据，生成结构、力拟合、DOS和自由能图片；自己的后处理总表在work/Co-common-free-energy.csv。

左图分开显示各项，右图比较直接相加与原生mkteci结果。图中的主曲线是固定体积，不能与准谐最优体积的振动项混用。

原胞DOS的静态E0为−7.0164067 eV/原子，与两原子FM参照相差约2.546 meV/原子，反映本示例有限网格的采样差异；这里只把原胞DOS的热修正接到选定参照，不把两套静态能量重复相加。用于精细材料比较时，应继续收敛网格和超胞。

最终你会得到：每个温度的静态参照、振动项、电子修正、有限磁性修正和总自由能。最近邻弹簧、FM背景与两格点磁性近似都应保留在结果说明中。