24小时热门版块排行榜    

查看: 4757  |  回复: 26

啸宇天涯

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
identation: 金币+2, 谢谢交流 2013-11-23 10:21:16
你体系79个原子,你为什么设置磁矩会那么多?
11楼2013-11-21 18:14:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anewtry

金虫 (小有名气)

引用回帖:
11楼: Originally posted by 啸宇天涯 at 2013-11-21 18:14:23
你体系79个原子,你为什么设置磁矩会那么多?

什么意思?我就是设置了79个磁矩啊,我算的是自旋轨道耦合,每个原子的磁矩给出要用三个值
12楼2013-11-21 19:32:14
已阅   回复此楼   关注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的回帖

anewtry

金虫 (小有名气)

引用回帖:
13楼: Originally posted by sxjsn1 at 2013-11-21 19:45:57
你这有点乱。。。  先自洽计算得到晶格结构(第一次计算)

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

这时INCAR中LSORBIT=.T(开自旋极化), 此时的磁矩,每三个数据为一 ...

你是说,在自洽计算时加入自旋轨道耦合,得到CHGCAR,然后后面计算能带和态密度时,直接用得到的CHGCAR,INCAR中不再设MAGMOM值了是吧?那我在自洽计算时,还设置ISPIN=2吗?还是直接设置MAGMOM值
14楼2013-11-21 20:10:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anewtry

金虫 (小有名气)

引用回帖:
9楼: Originally posted by sxjsn1 at 2013-11-21 16:26:28
你这是自洽计算完了算能带吧? 这一步的话磁矩是从CHAGCAR里读取,INCAR中不要设置了,不然会冲突。
我猜你初始优化晶格计算时,磁矩就是这么设置没出问题吧。...

我把INCAR改成了:SYSTEM=La2CoMnO6
ENCUT=550
ISTART=0
ICHARG=2
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=47*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 16*0
还是提示相同的错误
15楼2013-11-21 20:16:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

【答案】应助回帖

★ ★
identation: 金币+2, 应助指数+1, 谢谢交流 2013-11-23 10:22:11
引用回帖:
15楼: Originally posted by anewtry at 2013-11-21 20:16:26
我把INCAR改成了:SYSTEM=La2CoMnO6
ENCUT=550
ISTART=0
ICHARG=2
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 ...

你这个磁矩,不应该是47*0 0 0, 16*0 0 5, 16*0 0 0 吗? 为什么写成47*0, 16* 0 0 5, 16*0呢?
另外,磁矩是47*0 0 0 还是0 0 0*47这个我记不清了,你这个顺序没错吧?
夜半酒醒人去后,更持红烛赏残花
16楼2013-11-21 21:46:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
identation: 金币+2, 谢谢交流 2013-11-23 10:21:56
anewtry: 金币+5 2013-12-02 10:13:59
还有,我记得你前面说过,是48+16+16个原子,你就是这样设置原子数目也不对啊
夜半酒醒人去后,更持红烛赏残花
17楼2013-11-21 21:50:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

引用回帖:
17楼: Originally posted by sxjsn1 at 2013-11-21 21:50:38
还有,我记得你前面说过,是48+16+16个原子,你就是这样设置原子数目也不对啊

看错了,是79个,尴尬
夜半酒醒人去后,更持红烛赏残花
18楼2013-11-21 21:51:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anewtry

金虫 (小有名气)

引用回帖:
18楼: Originally posted by sxjsn1 at 2013-11-21 21:51:35
看错了,是79个,尴尬...

嗯,是79个,47+8+8+16.
我昨天又看了下说明书,上面说非线性的磁性,第一步先进行非磁性的结构优化,得到CHGCAR和WAVECAR。然后读得到的CHGCAR和WAVECAR,并加旋轨耦合的参数进行自洽,最后再算非自洽的band和dos.
我一开始使用ISPIN=2优化的,所以我觉得不太对,就想是不是先ISPIN=1 优化下,再试试。
还有就算是,说明书上给MAGMOM是000 000 001 这样给的。我在想我算的前面47 和后面16个磁矩为0 的MAGMOM能直接47*000 和16*000吗?请指教,谢谢
19楼2013-11-22 08:21:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjsn1

金虫 (小有名气)

veni vedi veci

★ ★
identation: 金币+2, 谢谢交流 2013-11-23 10:22:27
引用回帖:
19楼: Originally posted by anewtry at 2013-11-22 08:21:00
嗯,是79个,47+8+8+16.
我昨天又看了下说明书,上面说非线性的磁性,第一步先进行非磁性的结构优化,得到CHGCAR和WAVECAR。然后读得到的CHGCAR和WAVECAR,并加旋轨耦合的参数进行自洽,最后再算非自洽的band和do ...

我就是想说,应该是47*0 0 0和16*0 0 0, 你这又有0,又有0 0 5,格式都不一样。。。
还有,读取CHGCAR和WAVCAR进行自洽? 你确定? 除了算dos和band,其它应该都不需要读CHGCAR了。。。 而且,你读取CHGCAR后,INCAR设置磁矩,两者肯定冲突啊。。。
夜半酒醒人去后,更持红烛赏残花
20楼2013-11-22 21:36:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anewtry 的主题更新
信息提示
请填处理意见