| 查看: 3748 | 回复: 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值，不知对不对？？？还是应该取静态自洽完后的能量值？？？另外形成能和结合能一样吗？？？问题有点多？大家帮个忙。 |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
