| 查看: 4918 | 回復(fù): 26 | ||||||
anewtry金蟲 (小有名氣)
|
[求助]
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è),請各位賜教。謝謝。 |
第一性原理 | 第一性原理和電化學(xué) | 自旋軌道耦合 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
你這有點(diǎn)亂。。。 先自洽計(jì)算得到晶格結(jié)構(gòu)(第一次計(jì)算) 然后使用優(yōu)化后的POSCAR,新的INCAR,KPONTS和POTCAR重新計(jì)算(第二次計(jì)算)。 這時(shí)INCAR中LSORBIT=.T(開自旋極化), 此時(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ù)需刪掉。 |

至尊木蟲 (職業(yè)作家)
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
原子數(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)殚_啟了 spin-orbit calculation "每個(gè)"原子的 MOMENT 都有三個(gè)方向 所以正確寫法是 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) " \ " 也是沒有救的 如果你需要超過 就真的只能重新編輯 VASP了 |
木蟲 (職業(yè)作家)
金蟲 (小有名氣)
金蟲 (小有名氣)
|
我按照您的提示改了下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. 而且正好寫完一行,不用換行。 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
金蟲 (小有名氣)

金蟲 (小有名氣)
|
哦,我之前自洽時(shí)直接設(shè)置的ISPIN=2磁矩,磁矩設(shè)置是MAGMOM= 48*0.5 8*5 8*5 16*0.5。您能說說用VASP做自旋軌道耦合具體的步驟嗎?我看說明手冊做的,先在自旋情況下優(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ì)算。不知道這樣對不對? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +5 | 傅文秋 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 312求調(diào)劑 +3 | 小小墨123 2026-04-02 | 4/200 |
|
|
[考研] 材料調(diào)劑 +6 | 一樣YWY 2026-04-02 | 6/300 |
|
|
[考研] 一志愿同濟(jì)大學(xué)323分(080500)求調(diào)劑 +6 | yikeniu 2026-04-01 | 6/300 |
|
|
[考研] 學(xué)碩化學(xué)工程與技術(shù),一志愿中國海洋大學(xué)320+求調(diào)劑 +8 | 披星河 2026-04-02 | 8/400 |
|
|
[考研] 348求調(diào)劑 +6 | 吳彥祖24k 2026-04-02 | 6/300 |
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +8 | pipiver 2026-03-30 | 8/400 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +7 | Ong3 2026-03-27 | 7/350 |
|
|
[考研] 材料專碩322分 +11 | 哈哈哈吼吼吼哈 2026-04-01 | 11/550 |
|
|
[考研] 343求調(diào)劑085601 +4 | 要努力學(xué)習(xí)x 2026-03-29 | 5/250 |
|
|
[考研] 一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩,331分求調(diào)劑 +5 | 蔣昌鵬qtj 2026-04-01 | 5/250 |
|
|
[考研] 070300化學(xué)279求調(diào)劑 +15 | 哈哈哈^_^ 2026-03-31 | 17/850 |
|
|
[考研] 0817化工學(xué)碩調(diào)劑 +11 | 努力上岸中! 2026-03-31 | 11/550 |
|
|
[碩博家園] 考研調(diào)劑 +5 | 駱駝男人 2026-04-01 | 5/250 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +15 | minht 2026-03-29 | 16/800 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +6 | 小小機(jī)靈蟲 2026-03-29 | 12/600 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考研] 323分 食品與營養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|