同样做CO在Au表面的吸附能，为什么师姐做的总比你好？
小编之前给大家分享了吸附位点的探究和吸附结构的优化，很多小伙伴就很好奇吸附能是不是就可以直接从结果文件中得到。大家注意不要有这种想法，吸附能的计算要比这个稍微复杂一些，今天就和小编一起来学习计算CO吸附在Au(111)表面的吸附能。
首先，我们先来看一下吸附能的公式
Eads=EAu(111)+CO–EAu(111)–ECO
从上面这个公式我们可以看到，如果计算CO吸附在Au(111)表面的吸附能，我们需要得到的数据是E Au(111)+CO，E Au(111)，ECO的能量。
现在明确了计算目的，那么我们就接着上次的教程继续往下进行（Au(111)+CO的输入文件上期教程中已经提供，然后还需要准备Au(111)和CO的输入文件）
1.准备INCAR、POSCAR、KPOINTS、POTCAR输入文件
INCAR，POSCAR，KPOINTS，POTCAR准备就绪,POSCAR，KPOINTS，POTCAR就不多做赘述了，这里只展示INCAR文件。
【ps:注意】为了节省计算时间，这里在优化晶胞时没有打开LWAVE和LCHARG；ISIF=3是为了优化Au的离子位置、晶胞常数和体积。
2.查看结构是否收敛
优化完成之后，可以输入grep reached OUTCAR命令来查看结构是否收敛，如果我们看到reached required accuracy - stopping structural energy minimisation，就表明我们的优化任务已经完成了。
3.获取能量
现在就是获取能量，然后计算吸附能了。（在之前，可以把优化后的结构导出来看一下，即CONTCAR文件，这个是我们优化后的结构文件，然后用VESTA打开就可以看到结构了，也可有转化成.cif文件，用MS打开）
OUTCAR会展示出好几个能量，而我们一般计算用到的是energy(sigma->0) = 后面的那个，这和OSZICAR中E0后面的那个数值是一样的，感兴趣的小伙伴可以去翻一下这两个文件进行对比。
小编在这里使用的的提取命令是grep without OUTCAR | tail -n 1
4.计算吸附能
获得的能量总结如下：（小编保留了两位小数，大家是具体情况而定）
Au(111)+CO的能量为：-76.58 eV
Au(111)表面的能量为：-62.01 eV
CO的能量为：-14.80 eV
根据上面的公式计算得到CO分子的吸附能：Eads = 0.23 eV
到这里大家就已经得到了CO吸附在Au(111)表面的吸附能了，小伙伴们学会了吗，赶紧动手试一试。
END
往期推荐
点“在看”你最好看