24小时热门版块排行榜    

查看: 3385  |  回复: 5

jmlv

金虫 (著名写手)

[求助] VASP参数NUPDOWN如何使用? 已有1人参与

问题1:什么情况下需使用NUPDOWN?
问题2:启用NUPDOWN之后,还需设置ISPIN=2吗?
问题3:难道是设置一系列NUPDOWN值,取体系能量最低时的NUPDOWN值,也就是对NUPDOWN进行优化吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jmlv: 金币+10, ★★★很有帮助, 感谢大侠解惑。 2016-03-10 08:24:07
问题1:什么情况下需使用NUPDOWN?
答: 當你確定你的系統是有磁性,且你知道磁矩大小時

问题2:启用NUPDOWN之后,还需设置ISPIN=2吗?
答: 要。考慮 spin-polarization 的計算 (ISPIN) 才會用到 NUPDOWN

问题3:难道是设置一系列NUPDOWN值,取体系能量最低时的NUPDOWN值,也就是对NUPDOWN进行优化吗?
答: 你可以這樣做,但我相信很少人會這麼做


做磁性系統有兩個參數很好用
一個是 MAGMOM,你可以給各別原子的磁矩大小,然後VASP透過自洽計算會對磁矩做最佳化。MAGMOM 的預設值為 1,也就是每顆原子磁矩一開始是 1。
另一個就是 NUPDOWN,這參數是 spin up - spin down 的差值。如上面所講,VASP 會透過自洽對磁矩做最佳化,磁矩最佳化跟結構最佳化的過程是差不多的,也就是說路途上是有可能會陷入某個 local minimum,導致結果出現非基態的磁矩。因此,NUPDOWN 的作用就是為了避開這現像的發生。

當然,如果你做的是未知系統,磁矩多少你根本不曉得,或是你不知要怎麼去推測,那~ 是有可能會做像你問題三所講的事....
2楼2016-03-10 06:34:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmlv

金虫 (著名写手)

引用回帖:
2楼: Originally posted by jpchou at 2016-03-10 06:34:48
问题1:什么情况下需使用NUPDOWN?
答: 當你確定你的系統是有磁性,且你知道磁矩大小時

问题2:启用NUPDOWN之后,还需设置ISPIN=2吗?
答: 要。考慮 spin-polarization 的計算 (ISPIN) 才會用到 NUPDOWN

问 ...

大侠貌似是台湾同胞。根据你的回复,我追加提问两个问题。
问题4:一个作业中,设置NUPDOWN值并设ISPIN=2,VASP还将优化NUPDOWN值吗?
问题5:我计算的是金属氧化物的表面。我因为不知道我体系的磁矩,才打算用设置UPDOWN来推测体系的磁矩。理论上,体系能量最低时的NUPDOWN值,与仅设置ISPIN=2优化后的磁矩,应该是相等的,是吗?
3楼2016-03-10 08:31:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hy1987007

木虫 (著名写手)

引用回帖:
2楼: Originally posted by jpchou at 2016-03-10 06:34:48
问题1:什么情况下需使用NUPDOWN?
答: 當你確定你的系統是有磁性,且你知道磁矩大小時

问题2:启用NUPDOWN之后,还需设置ISPIN=2吗?
答: 要。考慮 spin-polarization 的計算 (ISPIN) 才會用到 NUPDOWN

问 ...

学习了,谢谢前辈分享!
4楼2016-03-10 12:57:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2016-03-11 20:18:49
jmlv: 金币+10, ★★★★★最佳答案, 谢谢大侠答疑! 2016-05-04 11:25:47
引用回帖:
3楼: Originally posted by jmlv at 2016-03-10 08:31:37
大侠貌似是台湾同胞。根据你的回复,我追加提问两个问题。
问题4:一个作业中,设置NUPDOWN值并设ISPIN=2,VASP还将优化NUPDOWN值吗?
问题5:我计算的是金属氧化物的表面。我因为不知道我体系的磁矩,才打算用设 ...

问题4:一个作业中,设置NUPDOWN值并设ISPIN=2,VASP还将优化NUPDOWN值吗?
答: 若 INCAR 內有出現 NUPDOWN,但 ISPIN = 1 (或是忘了填)
     那~ VASP 還是做 spin non-polarization 計算
     也就是 VASP 會直接忽略掉 NUPDOWN


问题5:我计算的是金属氧化物的表面。我因为不知道我体系的磁矩,才打算用设置UPDOWN来推测体系的磁矩。理论上,体系能量最低时的NUPDOWN值,与仅设置ISPIN=2优化后的磁矩,应该是相等的,是吗?
答: 請思考我回答 1 2 3 題後,之後寫的那些東西
     "理論上" , 只有 ISPIN=2 是可以算出系統能量最低(global minimum)
     的 NUPDOWN 值,NUPDOWN 的值也就是所謂的磁矩
     對大多數系統的確是可以的
     但,某些系統是 "有可能" 發生計算停在某個非最低能量的磁矩的態 (local minimum)
     所以,我只能說,要小心檢查計算結果
     做計算最擔心也是最為人垢病的情況就是所謂的
      garbage in, garbage out
     所以檢查結果合理性是理論計算人的基本要則之一
5楼2016-03-10 16:17:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
送红花一朵
赞!
6楼2017-03-21 09:34:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jmlv 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-11 09:20 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 3/150 2026-11-11 09:07 by 1Nb0IVlmV6v0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-11 08:34 by 1Nb0IVlmV6v0
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-11 06:57 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 9/450 2026-11-11 04:54 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 04:32 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +4 yexuqing 2026-11-04 6/300 2026-11-10 11:13 by susizheng
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
信息提示
请填处理意见