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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

嘯宇天涯

木蟲 (正式寫手)

應(yīng)助: 26 (小學(xué)生)
金幣: 3697
散金: 69
紅花: 6
帖子: 498
在線: 1098小時
蟲號: 2050401
注冊: 2012-10-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

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

你體系79個原子，你為什么設(shè)置磁矩會那么多？

贊一下

回復(fù)此樓

11樓2013-11-21 18:14:23

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

anewtry

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時
蟲號: 2072989
注冊: 2012-10-19
專業(yè): 無機材料化學(xué)

引用回帖:

11樓: Originally posted by 嘯宇天涯 at 2013-11-21 18:14:23
你體系79個原子，你為什么設(shè)置磁矩會那么多？

什么意思？我就是設(shè)置了79個磁矩啊，我算的是自旋軌道耦合，每個原子的磁矩給出要用三個值

贊一下

回復(fù)此樓

12樓2013-11-21 19:32:14

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

sxjsn1

金蟲 (小有名氣)

veni vedi veci

應(yīng)助: 5 (幼兒園)
金幣: 1108
散金: 165
紅花: 4
帖子: 238
在線: 101小時
蟲號: 941839
注冊: 2010-01-13
性別: GG
專業(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è)置的ISPIN=2磁矩，磁矩設(shè)置是MAGMOM= 48*0.5 8*5 8*5 16*0.5。您能說說用VASP做自旋軌道耦合具體的步驟嗎？我看說明手冊做的，先在自旋情況下優(yōu)化結(jié)構(gòu)，然后進行自洽計算，得到CHGCAR和WAVE ...

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

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

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

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

然后，你再用這個CHAGCAR計算band和DOS。

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

贊一下(2人)

回復(fù)此樓

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

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

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

anewtry

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時
蟲號: 2072989
注冊: 2012-10-19
專業(yè): 無機材料化學(xué)

引用回帖:

13樓: Originally posted by sxjsn1 at 2013-11-21 19:45:57
你這有點亂。。。先自洽計算得到晶格結(jié)構(gòu)（第一次計算）

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

這時INCAR中LSORBIT=.T（開自旋極化），此時的磁矩，每三個數(shù)據(jù)為一 ...

你是說，在自洽計算時加入自旋軌道耦合，得到CHGCAR，然后后面計算能帶和態(tài)密度時，直接用得到的CHGCAR，INCAR中不再設(shè)MAGMOM值了是吧？那我在自洽計算時，還設(shè)置ISPIN=2嗎？還是直接設(shè)置MAGMOM值

贊一下

回復(fù)此樓

14樓2013-11-21 20:10:35

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

anewtry

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時
蟲號: 2072989
注冊: 2012-10-19
專業(yè): 無機材料化學(xué)

引用回帖:

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

我把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
還是提示相同的錯誤

贊一下

回復(fù)此樓

15樓2013-11-21 20:16:26

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

sxjsn1

金蟲 (小有名氣)

veni vedi veci

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

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

★ ★
identation: 金幣+2, 應(yīng)助指數(shù)+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 ...

你這個磁矩，不應(yīng)該是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這個我記不清了，你這個順序沒錯吧？

贊一下

回復(fù)此樓

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

16樓2013-11-21 21:46:13

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

sxjsn1

金蟲 (小有名氣)

veni vedi veci

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

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

★ ★ ★ ★ ★ ★ ★
identation: 金幣+2, 謝謝交流 2013-11-23 10:21:56
anewtry: 金幣+5 2013-12-02 10:13:59

還有，我記得你前面說過，是48+16+16個原子，你就是這樣設(shè)置原子數(shù)目也不對啊

贊一下

回復(fù)此樓

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

17樓2013-11-21 21:50:38

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

sxjsn1

金蟲 (小有名氣)

veni vedi veci

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

引用回帖:

17樓: Originally posted by sxjsn1 at 2013-11-21 21:50:38
還有，我記得你前面說過，是48+16+16個原子，你就是這樣設(shè)置原子數(shù)目也不對啊

看錯了，是79個，尷尬

贊一下

回復(fù)此樓

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

18樓2013-11-21 21:51:35

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

anewtry

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1152.5
散金: 3
帖子: 155
在線: 104.8小時
蟲號: 2072989
注冊: 2012-10-19
專業(yè): 無機材料化學(xué)

引用回帖:

18樓: Originally posted by sxjsn1 at 2013-11-21 21:51:35
看錯了，是79個，尷尬...

嗯，是79個，47+8+8+16.
我昨天又看了下說明書，上面說非線性的磁性，第一步先進行非磁性的結(jié)構(gòu)優(yōu)化，得到CHGCAR和WAVECAR。然后讀得到的CHGCAR和WAVECAR，并加旋軌耦合的參數(shù)進行自洽，最后再算非自洽的band和dos.
我一開始使用ISPIN=2優(yōu)化的，所以我覺得不太對，就想是不是先ISPIN=1 優(yōu)化下，再試試。
還有就算是，說明書上給MAGMOM是000 000 001 這樣給的。我在想我算的前面47 和后面16個磁矩為0 的MAGMOM能直接47*000 和16*000嗎？請指教，謝謝

贊一下

回復(fù)此樓

19樓2013-11-22 08:21:00

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

sxjsn1

金蟲 (小有名氣)

veni vedi veci

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

★ ★
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.
我昨天又看了下說明書，上面說非線性的磁性，第一步先進行非磁性的結(jié)構(gòu)優(yōu)化，得到CHGCAR和WAVECAR。然后讀得到的CHGCAR和WAVECAR，并加旋軌耦合的參數(shù)進行自洽，最后再算非自洽的band和do ...

我就是想說，應(yīng)該是47*0 0 0和16*0 0 0，你這又有0，又有0 0 5，格式都不一樣。。。
還有，讀取CHGCAR和WAVCAR進行自洽？你確定？除了算dos和band，其它應(yīng)該都不需要讀CHGCAR了。。。而且，你讀取CHGCAR后，INCAR設(shè)置磁矩，兩者肯定沖突啊。。。

贊一下

回復(fù)此樓

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

20樓2013-11-22 21:36:11

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

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

返回列表

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

[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 295求調(diào)劑 +4	一志愿京區(qū)211 2026-03-18	6/300	2026-03-20 23:41 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 260求調(diào)劑 +3	朱芷琳 2026-03-20	3/150	2026-03-20 20:35 by 學(xué)員8dgXkO
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[基金申請] 學(xué)校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	8/400	2026-03-20 15:58 by babero
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 0703化學(xué) 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike