版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

下一頁(yè)

anewtry

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時(shí)
蟲(chóng)號(hào): 2072989
注冊(cè): 2012-10-19
專(zhuān)業(yè): 無(wú)機(jī)材料化學(xué)

[求助] VASP自旋軌道耦合

各位前輩，我在用VASP做自旋軌道耦合時(shí)，提示錯(cuò)誤：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

初次算這個(gè)，請(qǐng)各位賜教。謝謝。

回復(fù)此樓

1樓 2013-11-20 20:22:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

jpchou

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +20
1ST強(qiáng)帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時(shí)
蟲(chóng)號(hào): 585915
注冊(cè): 2008-08-08
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

★ ★
identation: 金幣+2, 謝謝交流 2013-11-23 10:20:26

MAGMOM 是給 initial moment
有多少原子就要給多少數(shù)字
但這個(gè)參數(shù)是有長(zhǎng)度限制
所以原子太多~ MAGMOM 寫(xiě)得太長(zhǎng)的話
就會(huì)出問(wèn)題
最簡(jiǎn)單的解決之道就是減少數(shù)字的出現(xiàn)
比方你有一百顆原子
每個(gè) initial moment 都是 0 的話
那就寫(xiě)成 MAGMOM = 100*0
這樣就可以有 100 個(gè) 0 了

贊一下(2人)

回復(fù)此樓

4樓2013-11-21 06:17:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sxjsn1

金蟲(chóng) (小有名氣)

veni vedi veci

應(yīng)助: 5 (幼兒園)
金幣: 1108
散金: 165
紅花: 4
帖子: 238
在線: 101小時(shí)
蟲(chóng)號(hào): 941839
注冊(cè): 2010-01-13
性別: GG
專(zhuān)業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

★ ★
identation: 金幣+2, 謝謝交流 2013-11-23 10:21:30

引用回帖:

10樓: Originally posted by anewtry at 2013-11-21 17:08:54
哦，我之前自洽時(shí)直接設(shè)置的ISPIN=2磁矩，磁矩設(shè)置是MAGMOM= 48*0.5 8*5 8*5 16*0.5。您能說(shuō)說(shuō)用VASP做自旋軌道耦合具體的步驟嗎？我看說(shuō)明手冊(cè)做的，先在自旋情況下優(yōu)化結(jié)構(gòu)，然后進(jìn)行自洽計(jì)算，得到CHGCAR和WAVE ...

你這有點(diǎn)亂。。。先自洽計(jì)算得到晶格結(jié)構(gòu)（第一次計(jì)算）

然后使用優(yōu)化后的POSCAR，新的INCAR,KPONTS和POTCAR重新計(jì)算（第二次計(jì)算）。

這時(shí)INCAR中LSORBIT=.T（開(kāi)自旋極化），此時(shí)的磁矩，每三個(gè)數(shù)據(jù)為一個(gè)原子三個(gè)軸方向的磁矩分量。比如，MAGMOM=0 0 5，表示一個(gè)原子x方向磁矩0，y方向磁矩0，z方向磁矩5。

設(shè)置好后，進(jìn)行計(jì)算，得到CHAGCAR.

然后，你再用這個(gè)CHAGCAR計(jì)算band和DOS。

看你主樓的INCAR，ICHARG=11,這種情況下是調(diào)用已有的CHAGCAR文件。你如果想產(chǎn)生新的CHAGCAR文件，這個(gè)參數(shù)需刪掉。

贊一下(2人)

回復(fù)此樓

夜半酒醒人去后，更持紅燭賞殘花

13樓2013-11-21 19:45:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲(chóng)號(hào): 1007127
注冊(cè): 2010-04-27
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
identation: 金幣+2, 謝謝交流 2013-11-23 10:20:35

MAGMOM設(shè)置時(shí)，要續(xù)行的話，用 \符號(hào)（反斜杠符號(hào)）。

贊一下(1人)

回復(fù)此樓

5樓2013-11-21 08:00:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +20
1ST強(qiáng)帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時(shí)
蟲(chóng)號(hào): 585915
注冊(cè): 2008-08-08
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+3, 謝謝指導(dǎo) 2013-11-23 17:01:50
anewtry: 金幣+3, ★★★很有幫助, 真的很感謝。 2013-12-02 10:10:04

原子數(shù)量 47+8+8+6=79
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

這樣的設(shè)置是錯(cuò)誤的
因?yàn)殚_(kāi)啟了 spin-orbit calculation
"每個(gè)"原子的 MOMENT 都有三個(gè)方向
所以正確寫(xiě)法是

MAGMOM=47*0 47*0 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 16*0 16*0
或
MAGMOM=141*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 48*0

MAGMOM 這行包頭包尾總共最多就是 255 個(gè)字元
你用換行符號(hào) " \ " 也是沒(méi)有救的
如果你需要超過(guò)
就真的只能重新編輯 VASP了

贊一下(1人)

回復(fù)此樓

22樓2013-11-23 15:46:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

zuocuiping

木蟲(chóng) (職業(yè)作家)

應(yīng)助: 14 (小學(xué)生)
金幣: 4893.9
散金: 966
紅花: 60
帖子: 4024
在線: 1597.4小時(shí)
蟲(chóng)號(hào): 890422
注冊(cè): 2009-11-01
性別: GG
專(zhuān)業(yè): 半導(dǎo)體光電子器件

可能要重新編譯一下vasp
你可以在小木蟲(chóng)上搜搜

贊一下

回復(fù)此樓

2樓2013-11-20 21:08:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anewtry

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時(shí)
蟲(chóng)號(hào): 2072989
注冊(cè): 2012-10-19
專(zhuān)業(yè): 無(wú)機(jī)材料化學(xué)

你是說(shuō)MAGMOM設(shè)置時(shí)給的數(shù)值太多，所以要編譯VASP嗎？之前也搜過(guò)，但不會(huì)編譯。有沒(méi)有可能不是這個(gè)原因呢

贊一下

回復(fù)此樓

3樓2013-11-20 21:11:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anewtry

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時(shí)
蟲(chóng)號(hào): 2072989
注冊(cè): 2012-10-19
專(zhuān)業(yè): 無(wú)機(jī)材料化學(xué)

引用回帖:

4樓: Originally posted by jpchou at 2013-11-21 06:17:49
MAGMOM 是給 initial moment
有多少原子就要給多少數(shù)字
但這個(gè)參數(shù)是有長(zhǎng)度限制
所以原子太多~ MAGMOM 寫(xiě)得太長(zhǎng)的話
就會(huì)出問(wèn)題
最簡(jiǎn)單的解決之道就是減少數(shù)字的出現(xiàn)
比方你有一百顆原子
每個(gè) initial momen ...

我按照您的提示改了下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=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
但是還是提示錯(cuò)誤：Error reading item 'MAGMOM' from file INCAR.
Error code was IERR=0 ... . Found N=  111 data.
而且正好寫(xiě)完一行，不用換行。

贊一下

回復(fù)此樓

6樓2013-11-21 09:05:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jpchou

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +20
1ST強(qiáng)帖: 5
應(yīng)助: 231 (大學(xué)生)
金幣: 2342
散金: 4367
紅花: 103
帖子: 1061
在線: 757.3小時(shí)
蟲(chóng)號(hào): 585915
注冊(cè): 2008-08-08
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

引用回帖:

6樓: Originally posted by anewtry at 2013-11-21 09:05:40
我按照您的提示改了下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 - ...

你有全部共有幾顆原子

贊一下

回復(fù)此樓

7樓2013-11-21 15:20:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anewtry

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時(shí)
蟲(chóng)號(hào): 2072989
注冊(cè): 2012-10-19
專(zhuān)業(yè): 無(wú)機(jī)材料化學(xué)

我體系里一共有79個(gè)原子。前面48個(gè)設(shè)置初始磁矩為0，最后16個(gè)原子設(shè)置初始磁矩為0.中間有16個(gè)原子設(shè)置的磁矩不為0.這是我預(yù)想的，是不是INCAR里磁矩還是設(shè)置的不對(duì)呀？

贊一下

回復(fù)此樓

8樓2013-11-21 15:47:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sxjsn1

金蟲(chóng) (小有名氣)

veni vedi veci

應(yīng)助: 5 (幼兒園)
金幣: 1108
散金: 165
紅花: 4
帖子: 238
在線: 101小時(shí)
蟲(chóng)號(hào): 941839
注冊(cè): 2010-01-13
性別: GG
專(zhuān)業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

【答案】應(yīng)助回帖

★ ★
感謝參與，應(yīng)助指數(shù) +1
identation: 金幣+2, 謝謝交流 2013-11-23 10:20:58

引用回帖:

8樓: Originally posted by anewtry at 2013-11-21 15:47:01
我體系里一共有79個(gè)原子。前面48個(gè)設(shè)置初始磁矩為0，最后16個(gè)原子設(shè)置初始磁矩為0.中間有16個(gè)原子設(shè)置的磁矩不為0.這是我預(yù)想的，是不是INCAR里磁矩還是設(shè)置的不對(duì)呀？

你這是自洽計(jì)算完了算能帶吧？這一步的話磁矩是從CHAGCAR里讀取，INCAR中不要設(shè)置了，不然會(huì)沖突。
我猜你初始優(yōu)化晶格計(jì)算時(shí)，磁矩就是這么設(shè)置沒(méi)出問(wèn)題吧。

贊一下

回復(fù)此樓

夜半酒醒人去后，更持紅燭賞殘花

9樓2013-11-21 16:26:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

anewtry

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時(shí)
蟲(chóng)號(hào): 2072989
注冊(cè): 2012-10-19
專(zhuān)業(yè): 無(wú)機(jī)材料化學(xué)

引用回帖:

9樓: Originally posted by sxjsn1 at 2013-11-21 16:26:28
你這是自洽計(jì)算完了算能帶吧？這一步的話磁矩是從CHAGCAR里讀取，INCAR中不要設(shè)置了，不然會(huì)沖突。
我猜你初始優(yōu)化晶格計(jì)算時(shí)，磁矩就是這么設(shè)置沒(méi)出問(wèn)題吧。...

哦，我之前自洽時(shí)直接設(shè)置的ISPIN=2磁矩，磁矩設(shè)置是MAGMOM= 48*0.5 8*5 8*5 16*0.5。您能說(shuō)說(shuō)用VASP做自旋軌道耦合具體的步驟嗎？我看說(shuō)明手冊(cè)做的，先在自旋情況下優(yōu)化結(jié)構(gòu)，然后進(jìn)行自洽計(jì)算，得到CHGCAR和WAVECAR，然后加入LSORBIT=.TRUE.
SAXIS=0 0 1
NBANDS=784
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
這些項(xiàng)，進(jìn)行非自洽的計(jì)算。不知道這樣對(duì)不對(duì)？

贊一下

回復(fù)此樓

10樓2013-11-21 17:08:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 anewtry 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑 +5	蔣昌鵬qtj 2026-04-01	5/250	2026-04-02 08:10 by fxue1114
[考研] 289求調(diào)劑 +21	新時(shí)代材料 2026-03-27	23/1150	2026-04-01 22:42 by peike
[考研] 考研調(diào)劑 +11	Amber00 2026-03-31	11/550	2026-04-01 11:32 by wangjy2002
[考研] 0855機(jī)械初試280求調(diào)劑 +3	kazenotori 2026-03-31	3/150	2026-04-01 10:08 by JourneyLucky
[考研] 求調(diào)劑：一志愿：南京大學(xué) 專(zhuān)業(yè)：0705 總分320 ，本科985，四六級(jí)已過(guò) +3	lfy760306 2026-03-31	3/150	2026-04-01 01:57 by Creta
[考研] 070300化學(xué)專(zhuān)業(yè)279調(diào)劑 +10	哈哈哈^_^ 2026-03-31	10/500	2026-03-31 23:13 by liu823948201
[考研] 080200學(xué)碩，機(jī)械工程專(zhuān)業(yè)277分，求帶走！ +4	瓶子PZ 2026-03-31	4/200	2026-03-31 20:16 by vgtyfty
[考研] 085602 307分求調(diào)劑 +10	不知道叫什么！ 2026-03-26	10/500	2026-03-31 19:53 by Dyhoer
[考研] 物理學(xué)調(diào)劑 +4	小羊36 2026-03-30	4/200	2026-03-31 16:16 by lishahe
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-31	4/200	2026-03-31 11:56 by jp9609
[考研] 22408 359分調(diào)劑 +4	Qshers 2026-03-27	8/400	2026-03-31 08:53 by Qshers
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4	XQTJZ 2026-03-30	4/200	2026-03-31 04:10 by fmesaito
[考研] 370求調(diào)劑 +3	080700調(diào)劑 2026-03-30	3/150	2026-03-31 01:09 by A_Zhe
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[有機(jī)交流] 考研調(diào)劑 +8	watb 2026-03-26	8/400	2026-03-30 18:40 by 544594351
[考研] 材料專(zhuān)碩 085600求調(diào)劑 +7	BBQ233 2026-03-30	7/350	2026-03-30 17:44 by oooqiao
[考研] 考研調(diào)劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息

24小時(shí)熱門(mén)版塊排行榜

[求助] VASP自旋軌道耦合

» 收錄本帖的淘帖專(zhuān)輯推薦

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖