24小时热门版块排行榜

返回列表

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
请各位指教！

回复此楼

» 猜你喜欢

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

求助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的回帖