24小时热门版块排行榜    

查看: 3337  |  回复: 3

小小鑫000

铁虫 (小有名气)

[求助] 会用VASP计算单个原子能量的高手进来! 已有2人参与

各位前辈大家好:
过渡金属Ni"(n=1一13)团簇"团簇中每个原子的平均束缚能EB采用如下公式计算:
EB=(NxE atom一Ecluster)/n     ;Eamot和Ecluster分别代表单个Ni原子的自由能和所研究原子团簇的自由能。
这是我从文献上看到的,我想计算一下单个Ni 原子的自由能。但不知道怎么算,于是从侯老师的手册上得知了两种计算原子能量的方法。
他说单个原子的计算有两个目的:1) 检验赝势的好坏;2) 对称性被破坏后自旋极化情况下的原子基态能量,对结合能进行修正。
然后就第一个目的,他INCAR设置如下:
SYSTEM = Al: atom
ENCUT = 250.00 eV   
NELMDL = 5
ISMEAR = 0; SIGMA=0.1
就第二个目的,INCAR设置为:
SYSTEM = Al: atom
ENCUT = 250  
ISYM = 0 ! no symmetry
ISPIN = 2 ! allow for spin polarisation
VOSKOWN = 1    ! this is important, in particular for GGA
ISMEAR = 0 ! Gaussian smearing, otherwise negative occupancies
SIGMA = 0.1 ! intermid. smearing width
AMIX = 0.2 ! mixing set manually
BMIX = 0.0001
NELM = 20 ! 20 electronic steps
ICHARG = 1
我不明白上面的意思,为什要分两种情况,主要是不知道我所要算的Ni原子的自由能到底该按上面哪种情况设置,还是采用别的设置方法。我这方面真是一塌糊涂,还望前辈细心解答,不胜感激!
回复此楼

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

好好研究,天天向上!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mywai520

铁杆木虫 (著名写手)


【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
小小鑫000: 金币+5 2014-05-04 09:31:33
若是你不做自旋极化直接自洽计算得到总能即可,若是自旋极化加上自旋项,然后自洽得到总能即可。
2楼2014-05-03 23:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by mywai520 at 2014-05-03 23:38:18
若是你不做自旋极化直接自洽计算得到总能即可,若是自旋极化加上自旋项,然后自洽得到总能即可。

谢谢前辈,前辈的意思是说无自旋按照第一种设置和有自旋按第二种方法设置INCAR 吗?

[ 发自手机版 http://muchong.com/3g ]
好好研究,天天向上!
3楼2014-05-04 09:35:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
小小鑫000: 金币+5 2014-05-05 10:06:45
?? ???} ???? PAPER ?????????????
????????v??????^???????????????e (reference)
?@???????]?п??^???????????????????????????


????W·????????????v?????????????W?
http://wenku.baidu.com/view/a98cfc2c4b73f242336c5ff3
???@?v?x???????? (???O?@???v?x??_???????????n???v?x)
???????????????????f (??????F?????`??)
???????^??  "  1) ???????????" ????]???P?S...... ????? ??
????N??  "?????????????????Ni??????????????????????????????" ?? ?


?????}????μ????????????`??
???`?????????????????????
?F????? & ????w?????? (?@?e????????????????....????^????w...)
???????????F?????
???
?? "??????"  ???????w???????????
(???????????С??)

??~???????v??????
?? "???? INCAR" ???????w???????????
4楼2014-05-05 03:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小鑫000 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +69 AI8RaGaPaSCR 2026-08-07 72/3600 2026-08-14 07:45 by u5Hy7ZmCpaka
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 29/1450 2026-08-14 07:42 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +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 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
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 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
信息提示
请填处理意见