| 查看: 3601 | 回复: 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值，不知对不对？？？还是应该取静态自洽完后的能量值？？？另外形成能和结合能一样吗？？？问题有点多？大家帮个忙。 |
回复此楼» 猜你喜欢
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
假如你的研究生提出不合理要求
已经有6人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
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