24小时热门版块排行榜    

查看: 4720  |  回复: 26
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

anewtry

金虫 (小有名气)

[求助] VASP自旋轨道耦合

各位前辈,我在用VASP做自旋轨道耦合时,提示错误:Error reading item 'MAGMOM' from file INCAR.
Error code was IERR=0 ... . Found N=  124 data.
我的INCAR如下:
SYSTEM=La2CoMnO6
ENCUT=550
ISTART=1
ICHARG=11
ISMEAR=-5
PREC=low
ISPIN=2
#IBRION=2
#POTIM=0.2
#ISIF=3
NELM=300
#NSW=200
LORBIT=11
LDAU=.TRUE.
LDAUL=-1  2 2 -1
LDAUU=0.5 5 5 0.5
LDAUJ=0  1  1  0
#LAECHG=.TRUE.
RWIGS=0.820 1.323 1.302 1.535
LSORBIT=.TRUE.
SAXIS=0 0 1
NBANDS=784
#LMAXMIX=1 4 4 1
MAGMOM=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

初次算这个,请各位赐教。谢谢。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

★ ★
identation: 金币+2, 谢谢交流 2013-11-23 10:21:30
引用回帖:
10楼: Originally posted by anewtry at 2013-11-21 17:08:54
哦,我之前自洽时直接设置的ISPIN=2磁矩,磁矩设置是MAGMOM= 48*0.5 8*5 8*5 16*0.5。您能说说用VASP做自旋轨道耦合具体的步骤吗?我看说明手册做的,先在自旋情况下优化结构,然后进行自洽计算,得到CHGCAR和WAVE ...

你这有点乱。。。  先自洽计算得到晶格结构(第一次计算)

然后使用优化后的POSCAR,新的INCAR,KPONTS和POTCAR重新计算(第二次计算)。

这时INCAR中LSORBIT=.T(开自旋极化), 此时的磁矩,每三个数据为一个原子三个轴方向的磁矩分量。 比如,MAGMOM=0 0 5,表示一个原子x方向磁矩0,y方向磁矩0,z方向磁矩5。

设置好后,进行计算,得到CHAGCAR.

然后,你再用这个CHAGCAR计算band和DOS。

看你主楼的INCAR,ICHARG=11,这种情况下是调用已有的CHAGCAR文件。你如果想产生新的CHAGCAR文件,这个参数需删掉。
夜半酒醒人去后,更持红烛赏残花
13楼2013-11-21 19:45:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 27 个回答

zuocuiping

木虫 (职业作家)

可能要重新编译一下vasp
你可以在小木虫上搜搜
2楼2013-11-20 21:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anewtry

金虫 (小有名气)

你是说MAGMOM设置时给的数值太多,所以要编译VASP吗?之前也搜过,但不会编译。有没有可能不是这个原因呢
3楼2013-11-20 21:11:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★
identation: 金币+2, 谢谢交流 2013-11-23 10:20:26
MAGMOM 是給 initial moment
有多少原子就要給多少數字
但這個參數是有長度限制
所以原子太多~ MAGMOM 寫得太長的話
就會出問題
最簡單的解決之道就是減少數字的出現
比方你有一百顆原子
每個 initial moment 都是 0 的話
那就寫成 MAGMOM = 100*0
這樣就可以有 100 個 0 了
4楼2013-11-21 06:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见