【摘要】 小编之前给大家分享了吸附位点的探究和吸附结构的优化，很多小伙伴就很好奇吸附能是不是就可以直接从结果文件中得到。

小编之前给大家分享了吸附位点的探究和吸附结构的优化，很多小伙伴就很好奇吸附能是不是就可以直接从结果文件中得到。大家注意不要有这种想法，吸附能的计算要比这个稍微复杂一些，今天就和小编一起来学习计算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)表面的吸附能了，小伙伴们学会了吗，赶紧动手试一试。

 

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