24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1379  |  回复: 4

dwd0826

铁杆木虫 (著名写手)


[交流] 计算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
请各位指教!
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+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楼2011-12-26 16:04:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dwd0826

铁杆木虫 (著名写手)


引用回帖:
: 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的回帖

lxbcsrc

铜虫 (初入文坛)


★ ★ ★ ★ ★ ★ ★ ★ ★
dwd0826: 金币+9, thx! 2012-12-29 10:32:18
好的
5楼2012-12-29 10:31:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2011-12-28 12:41   回复  
顶一下,感谢分享
相关版块跳转 我要订阅楼主 dwd0826 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见