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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP自旋轨道耦合
272/3返回列表上一页123下一页
查看: 5238  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 87/4350 2026-08-18 06:35 by 医学老男孩
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭