24小时热门版块排行榜    

查看: 1777  |  回复: 7

jmlv

金虫 (著名写手)

[求助] 请教:VASP如何得到Ni(111)表面的正确磁矩(magnetic moment)? 已有1人参与

我构建6x4x4的Ni(111)表面,真空层厚度15埃。欲得到正确的磁矩,我提交8个作业,作业之间的唯一不同是INCAR文件里MAGMOM的数值不同,第一个是MAGMOM=96*0.9,第二个是MAGMOM=96*1.0,一直到第八个MAGMOM=96*1.7。
优化完毕之后,能量最低,磁矩最大的那个,就是我们想要的。我这样优化磁矩的做法,是正确的吗?多谢。
回复此楼

» 猜你喜欢

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

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

wmy8802217

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
jmlv: 金币+4, ★★★很有帮助, 多谢。问题是,咱们不知道用哪个初始的MAMMOM值,才能是总能量最低呀。 2015-06-25 09:24:21
有必要么。。我怎么觉得设置一个就能收敛到能量最低的那个磁矩
凝聚态理论与计算物理
2楼2015-06-24 15:16:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

引用回帖:
2楼: Originally posted by wmy8802217 at 2015-06-24 15:16:08
有必要么。。我怎么觉得设置一个就能收敛到能量最低的那个磁矩

难道不在INCAR文件里明设MAGMOM值,就用MAGMOM的初始值?
3楼2015-06-25 09:25:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by jmlv at 2015-06-25 09:25:01
难道不在INCAR文件里明设MAGMOM值,就用MAGMOM的初始值?...

你设的只是初始值,最终系统会收敛到能量最低的那个磁矩大小的。。你不是设了那么多个,算完你看看最终磁矩是不是一样。
凝聚态理论与计算物理
4楼2015-06-25 09:29:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

引用回帖:
4楼: Originally posted by wmy8802217 at 2015-06-25 09:29:52
你设的只是初始值,最终系统会收敛到能量最低的那个磁矩大小的。。你不是设了那么多个,算完你看看最终磁矩是不是一样。...

我的计算显示,VASP没有那么智能。收敛后的磁矩和总能量,与MAGMOM初始值直接关联。如下表所示。
MAGMOM              mag_mom                Etot(eV)
96×0.9                78.53                -158.964
96×1.0                78.70                -158.897
96×1.1                77.51                -159.439
96×1.2                81.22                -171.965
96×1.3                76.75                -173.028               
96×1.4                90.06                -196.535               
96×1.5                94.39                -196.464               
96×1.6                80.52                -182.063               
96×1.7                80.52                -182.063
96×1.8                80.50                -182.064
96×1.9                80.53                -182.065
96×2.0                80.52                -182.063
5楼2015-06-25 16:10:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
jmlv: 金币+4, ★★★很有帮助, 我也咨询过业内同事同行,他们在做磁性单金属的表面计算时,貌似都在INCAR文件里明设MAGMOM值,只是设定ISPIN=2。得到的诸如吸附能等能量值,很接近文献报道的值。也就是说,对磁性单金属表面,在INCAR文件里设置MAGMOM的意义不大。 2015-06-26 08:28:15
引用回帖:
5楼: Originally posted by jmlv at 2015-06-25 16:10:10
我的计算显示,VASP没有那么智能。收敛后的磁矩和总能量,与MAGMOM初始值直接关联。如下表所示。
MAGMOM              mag_mom                Etot(eV)
96×0.9                78.53                -158.964
96×1.0                78.70                -158.897
96×1.1                77.51                -159.4 ...

据我师兄说,一般设的比收敛的大比较容易收敛到合适的磁矩,按这个说法,你后面几个应该算是比较合适的结果
凝聚态理论与计算物理
6楼2015-06-25 16:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

引用回帖:
6楼: Originally posted by wmy8802217 at 2015-06-25 16:45:46
据我师兄说,一般设的比收敛的大比较容易收敛到合适的磁矩,按这个说法,你后面几个应该算是比较合适的结果...

前帖“貌似都在INCAR文件里明设MAGMOM值”应该是“貌似都不在INCAR文件里明设MAGMOM值”。
7楼2015-06-26 08:29:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wmy8802217

木虫 (正式写手)

引用回帖:
7楼: Originally posted by jmlv at 2015-06-26 08:29:05
前帖“貌似都在INCAR文件里明设MAGMOM值”应该是“貌似都不在INCAR文件里明设MAGMOM值”。...

但是根据你测试的结果,还是有区别的
凝聚态理论与计算物理
8楼2015-06-26 09:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jmlv 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
信息提示
请填处理意见