同样做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)表面的吸附能了，小伙伴们学会了吗，赶紧动手试一试。

本文所有内容文字、图片和音视频资料，版权均属科学指南针网站所有，任何媒体、网站或个人未经本网协议授权不得以链接、转贴、截图等任何方式转载。

更多科研干货教程，可以点击下面链接获取哦~

模拟计算