| 查看: 3715 | 回复: 2 | ||
[求助]
单个原子结合能求解，到底指的哪一个？？？ 已有1人参与
|
|
小弟要计算过渡金属原子吸附在二维材料表面的问题，根据要求需要计算结合能，binding energy. Eb = E(TM) + E(graphene) − E(graphene + TM), where E(TM), E(graphene) and E(graphene + TM) are the energies of the isolated TM atom, the pristine graphene and the TM atom adsorbed graphene, respectively，TM为过渡金属。 现在单个过渡金属原子的能量，假设其为Co原子，小弟算出来3种： 1. 不开自旋， ISPIN=1, POSCAR为 1.0 8 0 0 0 8 0 0 0 8 1 direct 0 0 0 KPOINTS为 0 Monkhost pack 1 1 1 0 0 0 POTCAR用的PBE-PAW INCAR中ISMEAR=0,,,,,,,,,这样优化完结果为： free TOTEN= -0.45900580 energy without entropy=-0.37450127 energy(sigma->0)=-0.41675353 2. 开自旋，ISPIN=2, POSCAR晶格常数为 8 0 0， 9 0 0 ，10 0 0，其他设置与1中情况一样，优化完结果为 free TOTEN=-1.75919912， energy without entropy=-1.72248815 energy(sigma->0)=-1.74084364， 3.开自旋，这次我用的是VASP官网的O原子的例子，INCAR中只写了:SYSTEM=Co,ISMEAR=0,ISPIN=2,POSCAR晶格常数为 8 0 0， 9 0 0 ，10 0 0算出来结果free TOTEN=-1.84363826energy without entropy=-energy-1.67698926(sigma->0)=-1.76031376， 这3种情况应该以哪一种为准，希望广大虫友给予解答，由于ISMEAR=0，所以我总能取得energy without entropy=-0.37450127，不知对不对？？？？？？ 还是应该取free TOTEN的值，这3个能量到底什么区别，再一个，这些能量为什么都是负值？？还有就是 E(graphene)和 E(graphene + TM)我取得是结构优化完后OUTCAR中的energy without entropy值，不知对不对？？？还是应该取静态自洽完后的能量值？？？另外形成能和结合能一样吗？？？问题有点多？大家帮个忙。 |
回复此楼» 猜你喜欢
面上的水很深
已经有16人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
国社科通讯评审
已经有4人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2017-05-25 16:44:01
【答案】应助回帖
|
这3种情况应该以哪一种为准，希望广大虫友给予解答，由于ISMEAR=0，所以我总能取得energy without entropy=-0.37450127，不知对不对？？？？？？ 还是应该取free TOTEN的值，这3个能量到底什么区别，再一个，这些能量为什么都是负值？？ 能量负值没问题，自旋肯定要开 还有就是 E(graphene)和 E(graphene + TM)我取得是结构优化完后OUTCAR中的energy without entropy值，不知对不对？？？还是应该取静态自洽完后的能量值？？？另外形成能和结合能一样吗？？？问题有点多？大家帮个忙。 取OSZICAR里的能量即可 结合能与形成能不一样，多看文献 |
3楼2018-01-30 18:02:22
10