24小时热门版块排行榜

>论坛更新日志 (4903)
>博后之家 (1070)
>虫友互识 (485)
>导师招生 (434)
>论文投稿 (277)
>文献求助 (248)
>学术会议 (213)
>休闲灌水 (169)
>功能材料 (116)
>招聘信息布告栏 (105)
>硕博家园 (77)
>考研 (71)
>考博 (70)
>基金申请 (58)
>教师之家 (53)
>攻关文献(高奖励) (43)

返回列表

dwd0826

铁杆木虫 (著名写手)

应助: 7 (幼儿园)
金币: 6550.9
帖子: 1329
在线: 427小时
虫号: 1100436

[交流] 计算Fe的超晶胞晶格常数与能量

刚计算了Fe54超晶胞的晶格常数与能量，通过状态方程拟合，得到的V0＝614.403(即a0＝8.5013)，total enthalpy＝-443.214eV，参考某篇文献中的a0＝8.6001，total enthalpy＝-443.295eV。可以看到，能量较为符合，但最后得到的a0值有将近0.1A的差异，是否存在问题？
Pseudopotential：GGA-PAW91；以下是INCAR：
SYSTEM=Fe54
ENCUT=520
ISTART=0;ICHARG=2
VOSKOWN=1
IALGO=48
ISPIN=2
ISMEAR=-5
PREC=Accurate
请各位指教！

回复此楼

» 猜你喜欢

售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有3人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有4人回复
面上没中，邀请各位路过的虫友分析一下分数已经有13人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有8人回复
两块石头已经有8人回复
初秋的晨风已经有6人回复
学科评审组评审是指会评吗？已经有9人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有14人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有7人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有17人回复

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

求助castep计算中：晶胞总能量E 与晶胞体积V 的关系图如何得到，看内容已经有24人回复
晶格常数和能量是这种形式？已经有3人回复
【求助】晶格常数大了以后 VASP的计算值准确吗已经有5人回复
【求助成功】晶格优化的结果很诡异已经有13人回复
【求助成功】如何用VASP计算反铁磁已经有14人回复
【原创】如何计算中性杂质的形成能已经有66人回复
【转帖】如何用VASP计算晶格常数已经有8人回复

1楼 2011-12-26 15:24:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

WDD880227

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
fzx2008(金币+1): 谢谢交流 2011-12-26 16:08:10
dwd0826(金币+1): 感谢解答！ 2011-12-26 19:00:09

ISPIN=2
想要计算自旋的话只设置ISPIN=2貌似是不够的，还需要MAGMOM和NUPDOWN吧~~

赞一下(2人)

回复此楼

2楼2011-12-26 16:04:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dwd0826

铁杆木虫 (著名写手)

应助: 7 (幼儿园)
金币: 6550.9
帖子: 1329
在线: 427小时
虫号: 1100436

引用回帖:

楼: Originally posted by WDD880227 at 2011-12-26 16:04:01:
ISPIN=2
想要计算自旋的话只设置ISPIN=2貌似是不够的，还需要MAGMOM和NUPDOWN吧~~

在计算纯Fe体系时MAGMOM应该默认为NIONS*1吧？另外，关于NUODOWN，手册上描述为:"If no value is set(or NUODOWN=-1)a full relaxation will be performed. this is also the default."指的是默认时不固定磁矩吗？是否在这里也应该是不固定磁矩，通过静态计算结果拟合后得到基态能量与磁矩？

赞一下

回复此楼

3楼2011-12-26 19:12:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖