应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP自旋轨道耦合
272/3返回列表上一页123下一页
查看: 4757  |  回复: 26
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

啸宇天涯

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +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文件，这个参数需删掉。
赞一下(2人) 回复此楼
夜半酒醒人去后，更持红烛赏残花
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 的主题更新
272/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭