24小时热门版块排行榜    

查看: 2289  |  回复: 12

小小鑫000

铁虫 (小有名气)

[求助] 急!请教高手用vasp 优化Ni 团簇遇到的问题~ 已有1人参与

各位前辈,大家好:
最近用vasp 优化Ni13团簇结构,还优化将氢气吸附到团簇的某个原子上的结构,但是我有一个疑问,就是据说Ni 是磁性材料,那我在优化这个结构时,INCAR 中需不需要加ISPIN 和MAGMOM 两个条件限制,加这两个参数和不加,对我最后优化的结构有没有影响?

[ 发自手机版 http://muchong.com/3g ]
回复此楼

» 本帖@通知

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

好好研究,天天向上!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
小小鑫000: 金币+5, ★★★★★最佳答案 2014-06-22 16:23:35
fzx2008: 金币+2, 谢谢指导 2014-06-22 21:55:57
ISPIN 需要
MAGMOM 視情況~可加可不加

不加這兩個參數
你計算出?淼木褪菦]有磁性的 Ni13
對最後優化的結構應該會些許有些影響
2楼2014-06-22 04:21:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by jpchou at 2014-06-22 04:21:23
ISPIN 需要
MAGMOM 視情況~可加可不加

不加這兩個參數
你計算出?淼木褪菦]有磁性的 Ni13
對最後優化的結構應該會些許有些影響

多谢前辈指导,我在优化完这些个团簇结构时,需要知道其总能量,从网上得知需要再进行一次静态计算,我的INCAR设置如下,我老感觉自己的能量算的不准确,前辈帮我看看INCAR的参数设置是否合理,还有ISPIN需要加吧?MAGMOM 也不需要加了吧?谢谢前辈了!
SYSTEM=Ni
ISTART=0
ICHARG=2
PREC=Accurate
ENCUT=450
ISMEAR=0
SIGMA=0.1
ISPIN=2
好好研究,天天向上!
3楼2014-06-22 16:30:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

??л???????????????????Щ???????????????????????????????????????????ξ?????????INCAR?????????????о????????????????????????????INCAR?????????????????????ISPIN???????MAGMOM ????????????лл??????
SYSTEM=Ni
ISTART=0
ICHARG=2
PREC=Accurate
ENCUT=450
ISMEAR=0
SIGMA=0.1
ISPIN=2



"??о????????????????" ??? ????N??????@????X????
???????c
?? PREC ?{???c
EDIFF ?? EDIFFG ?{???c
???????????
INCAR ????????????
4楼2014-06-22 18:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by jpchou at 2014-06-22 18:08:30
多谢前辈指导,我在优化完这些个团簇结构时,需要知道其总能量,从网上得知需要再进行一次静态计算,我的INCAR设置如下,我老感觉自己的能量算的不准确,前辈帮我看看INCAR的参数设置是否合理,还有ISPIN需要加吧? ...

前辈,我语言表达不清楚,还望见谅,我的意思是我的INCAR里也不需要加MAGMOM这个参数了吧,不加这个参数对最后的总能量大小不影响吧?
我还有第二的问题就是:这个总能量我该取哪个值最好,是TOTEN那个还是energy without entropy还是energy(sigma趋向于零)这个值,因为我算出来的TOTEN和energy without entropy的值相差很大,不知道哪个是总能量值?
好好研究,天天向上!
5楼2014-06-22 20:06:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

???????:
5?: Originally posted by СС??000 at 2014-06-22 20:06:18
???????????????????????????????????????INCAR??????????MAGMOM????????????????????????????????????С???????
????е???????????????????????????????????TOTEN???????energy ...

????? MAGMOM ?@??????????????
??????o ISPIN=2??VASP ???? spiin-polarization ???
VASP ???A?O???????? moment = 1
??^ SCF ??VASP ??????????? moment
??r??Y????o?????C????????_??
??? moment space ?? energy space ???????
?????? global minimum ?? local minimum

????f
???O????o ISPIN=2???]?нo MAGMOM
??????y (13 ?w Ni) ???? moment ???? 13
????????????~ "????" ????? 3
????~ Ni13 ?????? moment=1 (???????....)
?@???~??_???o ISPIN=2 ??????
???o MAGMOM=13*0.76
???_??? moment ????? 1
???~ ???????????? moment=1 ??Y??


moment ???~????????????
?????????? moment ~ ????????????

?? TOTEN ??? energy without entropy ?
?@??????
?@Щ????????????"??????"
???~??????????? formation energy ?@?????
??~ ?????????
(formation energy ??? Ni bulk energy
Ni bulk ????? TOTEN?? ?? Ni13 ???? TOTEN?? ?@????"???"
6楼2014-06-23 02:09:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)

引用回帖:
6楼: Originally posted by jpchou at 2014-06-23 02:09:53
你知道 MAGMOM 這個參數的用途嬤?
當你只給 ISPIN=2,VASP 會做 spiin-polarization 計算
VASP 會預設所有原子的 moment = 1
透過 SCF 後,VASP 能得到所有原子 moment
此時的結果是無法保證完全是正確的
因 ...

谢谢前辈的详细解答~前辈的意思就是说MAGMOM这个参数不影响结构优化的结果,但会影响最后的总能量,是这样的吧?

我最终是要算吸附能,吸附能的计算公式中是要算几个部分的总能,这个总能也是取哪个能量都行吧?

[ 发自手机版 http://muchong.com/3g ]
好好研究,天天向上!
7楼2014-06-23 16:16:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

???????:
7?: Originally posted by СС??000 at 2014-06-23 16:16:17
лл????????????????????????MAGMOM???????????????????????????????????????????????????

?????????????????????????????????????????????????????????????????????а? ...

MAGMOM ??????Y??   "??"   ??
?? "?????" ...
?????????y?????
??????????y
???????

??t??
??? ENTOT ?? energy without entropy ??????
????????~???
??~~???
??????@?????????????
???? "TOTEN??energy without entropy????????" ~ ????С??~???l????N??????
???^~??f
?????^??"?????" ??????? "??"?  ??? eV?
????ǔ???~??????? meV ~~ ???????]??...........
8楼2014-06-23 17:30:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)

引用回帖:
8楼: Originally posted by jpchou at 2014-06-23 17:30:07
MAGMOM 不是對結構   "不"   影響
是 "影響不大" ...
對大多數系統影響不大
磁性?姷南到y
影響就會很大

原則上
你取 ENTOT 或 energy without entropy 都可以
每一項都取一致~就好
...

前辈,我的有的相差0.4ev了,这个算不算大呀?正不正常?
好好研究,天天向上!
9楼2014-06-23 18:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小鑫000

铁虫 (小有名气)

引用回帖:
8楼: Originally posted by jpchou at 2014-06-23 17:30:07
MAGMOM 不是對結構   "不"   影響
是 "影響不大" ...
對大多數系統影響不大
磁性?姷南到y
影響就會很大

原則上
你取 ENTOT 或 energy without entropy 都可以
每一項都取一致~就好
...

前辈,我在静态计算Ni4的能量时,第一次只在INCAR中加入了ISPIN=2这个参数,没有加MAGMOM,也就是让每个的磁矩都默认为1;第二次,我还是加了ISPIN=2,这时令MAGMOM=4*2;第三次也加了ISPIN=2,这时令MAGMOM=4*0.75,结果很奇怪,这三次最后OSICAR里的总能量相同,磁矩也都相同,都是MAG=3.6897,这是怎么回事呀?前辈?
好好研究,天天向上!
10楼2014-06-23 18:47:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小小鑫000 的主题更新
信息提示
请填处理意见