同样做CO在Au表面的吸附能,为什么师姐做的总比你好? | 您所在的位置:网站首页 › cif计算公式是什么 › 同样做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)表面的吸附能了,小伙伴们学会了吗,赶紧动手试一试。 本文所有内容文字、图片和音视频资料,版权均属科学指南针网站所有,任何媒体、网站或个人未经本网协议授权不得以链接、转贴、截图等任何方式转载。 来源:同样做CO在Au表面的吸附能,为什么师姐做的总比你好? |
CopyRight 2018-2019 实验室设备网 版权所有 |