24小时热门版块排行榜    

查看: 3747  |  回复: 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分析方法 +6 majunge000 2026-08-10 6/300 2026-08-12 20:43 by JChen2024
[基金申请] 不应该看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 云上清扬
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看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
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 确定了,国自然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
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +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 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于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 +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
请填处理意见