版塊導(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í)熱門版塊排行榜

返回列表

下一頁(yè)

本帖產(chǎn)生 1 個(gè) 1ST強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

langren24

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 127.5
散金: 11
帖子: 193
在線: 96.4小時(shí)
蟲號(hào): 1136177
注冊(cè): 2010-11-01
專業(yè): 新型信息器件

[求助] spin-orbit計(jì)算求助！

最近用spin-orbit計(jì)算體系的磁性總是出錯(cuò)，現(xiàn)將計(jì)算過程中的一些問題和疑惑貼出來(lái)，求大蝦解答：
1. 用spin-orbit與spin-polarizedINCAR文件的不同有哪些？除了需要加
LSORBIT=.TRUE.
SAXIS =X Y Z
NBANDS變?yōu)楣簿€的兩倍，MAGMON也變?yōu)閤yz方向之外
需不需要加
GGA_COMPAT=.FALSE.
ISYM =0
呢？（這是在VASP手冊(cè)上看到的）
2.自旋量子軸SAXIS =X Y Z 應(yīng)該怎么設(shè)置啊？應(yīng)該與體系晶格結(jié)構(gòu)有關(guān)，軸的設(shè)置應(yīng)該根據(jù)實(shí)空間晶格設(shè)置呢，還是根據(jù)倒空間晶格設(shè)置呢，求舉例！
3.根據(jù)VASP手冊(cè)上提供的方法，計(jì)算spin-orbit之前應(yīng)該先計(jì)算共線情況下的WAVECAR和CHGCAR,然后再開始spin-orbit的計(jì)算。
但是我在計(jì)算的時(shí)候?qū)胫坝?jì)算的WAVECAR和CHGCAR，計(jì)算總會(huì)出錯(cuò)，不知道是什么原因？是不是應(yīng)該計(jì)算非共線情況下的WAVECAR和CHGCAR，然后再計(jì)算spin-orbit呢？
4.計(jì)算spin-orbit需不需要使用超胞呢？（摻雜或者不摻雜的情況下）
使用超胞的計(jì)算結(jié)果該怎么分析呢？
5.計(jì)算spin-orbit對(duì)贗勢(shì)（譬如GGA-PBE等）有沒有要求，對(duì)算法（譬如ALGO=Normal、ALGO=Damped等），對(duì)計(jì)算精度等有沒有要求？

[ Last edited by langren24 on 2011-12-13 at 15:28 ]

回復(fù)此樓

1樓 2011-12-13 15:24:44

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

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

zhang668

木蟲 (正式寫手)

1ST強(qiáng)帖: 3
應(yīng)助: 11 (小學(xué)生)
貴賓: 0.061
金幣: 1735
散金: 1294
紅花: 11
帖子: 843
在線: 324.9小時(shí)
蟲號(hào): 733242
注冊(cè): 2009-03-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
youzhizhe(金幣+5, 1ST強(qiáng)帖+1): 謝謝交流。授予EPI一枚。計(jì)入專家考核存檔。 2011-12-13 20:30:15
langren24(金幣+2): 2011-12-14 14:13:18

1. ISYM =0在打開SO時(shí)要加，用GGA交換泛函時(shí)需要用GGA_COMPAT=.FALSE.改善精度。
2. SAXIS和MAGMOM有兩種常用設(shè)置，手冊(cè)上有介紹，默認(rèn)是0 0 1方向，是相對(duì)實(shí)空間笛卡爾坐標(biāo)系的。一種是SAXIX為0 0 1方向，此時(shí)MAGMOM為任意方向，另外一種是MAGMOM均為 0 0 z或0 0 -z方向，SAXIX為任意方向，第二種允許讀取先前存在的共線或非共線電荷密度和波函數(shù)，且精度更高。
3. 導(dǎo)入之前計(jì)算的WAVECAR和CHGCAR出錯(cuò)原因是KPOINTS在SO打開時(shí)對(duì)稱性發(fā)生變化。你可以用SO生成的IBZKPT,然后拷貝為KPOINTS,重新計(jì)算WAVECAR和CHGCAR，然后再讀取就可以了。
4. 不需要用超胞！
5. SO只支持PAW贗勢(shì)，對(duì)精度要求高，能量精度EDIFF一般為10^-7eV.

贊一下(4人)

回復(fù)此樓

3樓2011-12-13 20:17:54

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

普通回帖

fzx2008

榮譽(yù)版主 (著名寫手)

1ST強(qiáng)帖: 9
應(yīng)助: 272 (大學(xué)生)
貴賓: 1.045
金幣: 21814.1
散金: 710
紅花: 33
沙發(fā): 1
帖子: 1808
在線: 1601小時(shí)
蟲號(hào): 621341
注冊(cè): 2008-10-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

關(guān)注！

回復(fù)此樓

2樓2011-12-13 15:32:44

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

雪狼1210

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 919.7
散金: 43
帖子: 128
在線: 155.9小時(shí)
蟲號(hào): 939826
注冊(cè): 2010-01-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★
感謝參與，應(yīng)助指數(shù) +1
fzx2008(金幣-1): 非應(yīng)助請(qǐng)不要點(diǎn)應(yīng)助回帖 2011-12-13 22:04:47

SO的強(qiáng)度怎么計(jì)算呢？有的文獻(xiàn)中討論SO強(qiáng)度變化（0%至100%），這個(gè)怎么實(shí)現(xiàn)？還請(qǐng)高手指點(diǎn)

贊一下(1人)

回復(fù)此樓

4樓2011-12-13 21:02:05

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

spatial

榮譽(yù)版主 (職業(yè)作家)

應(yīng)助: 1 (幼兒園)
貴賓: 17.95
金幣: 2915.7
帖子: 3871
在線: 44.7小時(shí)
蟲號(hào): 46929
注冊(cè): 2004-05-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

收藏，下一步準(zhǔn)備so去

贊一下

回復(fù)此樓

http://spatial.blogchina.com/

5樓2011-12-14 12:41:05

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

langren24

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 127.5
散金: 11
帖子: 193
在線: 96.4小時(shí)
蟲號(hào): 1136177
注冊(cè): 2010-11-01
專業(yè): 新型信息器件

引用回帖:

3樓: Originally posted by zhang668 at 2011-12-13 20:17:54:
1. ISYM =0在打開SO時(shí)要加，用GGA交換泛函時(shí)需要用GGA_COMPAT=.FALSE.改善精度。
2. SAXIS和MAGMOM有兩種常用設(shè)置，手冊(cè)上有介紹，默認(rèn)是0 0 1方向，是相對(duì)實(shí)空間笛卡爾坐標(biāo)系的。一種是SAXIX為0 0 1方向，此時(shí)M ...

spin-orbit對(duì)算法（譬如ALGO=Normal、ALGO=Fast、ALGO=Damped等）有沒有要求呢？我在計(jì)算的時(shí)候換不同的算法總會(huì)出現(xiàn)不收斂的情況，用什么樣的算法最合適呢？

贊一下

回復(fù)此樓

6樓2011-12-14 14:24:52

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

langren24

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 127.5
散金: 11
帖子: 193
在線: 96.4小時(shí)
蟲號(hào): 1136177
注冊(cè): 2010-11-01
專業(yè): 新型信息器件

引用回帖:

非常感謝您的解答，我受益匪淺，還有幾個(gè)問題，都寫在一起，請(qǐng)大蝦解答：
1.spin-orbit對(duì)算法（譬如ALGO=Normal、ALGO=Fast、ALGO=Damped等）有沒有要求呢？我在計(jì)算的時(shí)候換不同的算法總會(huì)出現(xiàn)不收斂的情況，用什么樣的算法最合適呢？
2.spin-orbit用LDA+U或者LSDA+U有什么不同嗎，或者哪種方法更適合s-o？
3.用spin-orbit計(jì)算MAE的話，對(duì)于MAGMOM和SAXIS的設(shè)置用第一種還是第二種設(shè)置合適呢，因?yàn)橛?jì)算MAE需要改變SAXIX的方向？

贊一下

回復(fù)此樓

7樓2011-12-14 14:52:02

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

zhang668

木蟲 (正式寫手)

1ST強(qiáng)帖: 3
應(yīng)助: 11 (小學(xué)生)
貴賓: 0.061
金幣: 1735
散金: 1294
紅花: 11
帖子: 843
在線: 324.9小時(shí)
蟲號(hào): 733242
注冊(cè): 2009-03-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
sunyang1988(金幣+2): 謝謝指點(diǎn) 2011-12-14 17:17:35
ben_ladeng: 專家考核存檔 2011-12-15 17:41:52

1. SO對(duì)算法應(yīng)該沒要求，因?yàn)槲襂ALGO 38（DAV）或AlGO=Fast都算成功過。
2. SO肯定要打開自旋的，是LSDA+U。
3. 第二種，改變SAXIS相當(dāng)于改變磁場(chǎng)方向，以此來(lái)計(jì)算磁各向異性能。

贊一下(2人)

回復(fù)此樓

8樓2011-12-14 16:08:28

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

雪狼1210

銀蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 919.7
散金: 43
帖子: 128
在線: 155.9小時(shí)
蟲號(hào): 939826
注冊(cè): 2010-01-09
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

4樓: Originally posted by 雪狼1210 at 2011-12-13 21:02:05:
SO的強(qiáng)度怎么計(jì)算呢？有的文獻(xiàn)中討論SO強(qiáng)度變化（0%至100%），這個(gè)怎么實(shí)現(xiàn)？還請(qǐng)高手指點(diǎn)

不好意思，點(diǎn)錯(cuò)了，我的問題可否指點(diǎn)一下

贊一下

回復(fù)此樓

9樓2011-12-14 16:22:04

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

zhang668

木蟲 (正式寫手)

1ST強(qiáng)帖: 3
應(yīng)助: 11 (小學(xué)生)
貴賓: 0.061
金幣: 1735
散金: 1294
紅花: 11
帖子: 843
在線: 324.9小時(shí)
蟲號(hào): 733242
注冊(cè): 2009-03-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

ben_ladeng: 專家考核存檔 2011-12-15 17:42:41

引用回帖:

9樓: Originally posted by 雪狼1210 at 2011-12-14 16:22:04:
不好意思，點(diǎn)錯(cuò)了，我的問題可否指點(diǎn)一下

你是指認(rèn)為調(diào)節(jié)參數(shù)來(lái)改變SO強(qiáng)度？這個(gè)跟勢(shì)能V（r）有關(guān)。我沒見到這么調(diào)節(jié)的

贊一下(1人)

回復(fù)此樓

10樓2011-12-14 16:49:34

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

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

返回列表

下一頁(yè)

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

[考研] 0703化學(xué)調(diào)劑 +11	妮妮ninicgb 2026-03-15	15/750	2026-03-21 19:15 by ColorlessPI
[考研] 求調(diào)劑 +4	要好好無(wú)聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 306求0703調(diào)劑一志愿華中師范 +5	紙魚ly 2026-03-21	5/250	2026-03-21 17:11 by 學(xué)員8dgXkO
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 299求調(diào)劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	5/250	2026-03-20 22:11 by 云游重陽(yáng)
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[考研] 生物學(xué)調(diào)劑招人�。�！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu

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

[求助] spin-orbit計(jì)算求助！

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

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

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

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

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

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

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