24小时热门版块排行榜    

查看: 1688  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +24 医学老男孩 2026-08-13 50/2500 2026-08-15 06:38 by 学员8dgXkO
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 04:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 4/200 2026-08-15 01:08 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 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
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见