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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP自旋轨道耦合
272/3返回列表上一页123下一页
查看: 4886  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中科技大学071000，求调剂 +3 沿岸有贝壳6 2026-03-21 3/150 2026-03-21 10:35 by 暮云清寒
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +6 墨墨漠 2026-03-20 6/300 2026-03-21 09:23 by 天青等烟雨
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 一志愿中国石油大学（华东） 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸！ 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[论文投稿] 申请回稿延期一个月，编辑同意了。但系统上的时间没变，给编辑又写邮件了，没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
关闭
请填处理意见
关闭