24小时热门版块排行榜    

查看: 3746  |  回复: 2

billvasp

新虫 (小有名气)

[求助] 单个原子结合能求解,到底指的哪一个??? 已有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值,不知对不对???还是应该取静态自洽完后的能量值???另外形成能和结合能一样吗???问题有点多?大家帮个忙。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

770968942DX

新虫 (初入文坛)

请问这个问题解决了吗?我最近也遇到了这个问题
2楼2017-05-25 16:44:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

前田庆次

木虫 (著名写手)

【答案】应助回帖

这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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 billvasp 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 21/1050 2026-08-12 16:08 by 云上清扬
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
信息提示
请填处理意见