dx0620

木蟲 (正式寫手)

應(yīng)助: 14 (小學(xué)生)
金幣: 4053.6
散金: 5
紅花: 6
帖子: 393
在線: 910小時(shí)
蟲號(hào): 1303311
注冊(cè): 2011-05-22
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] vasp 考慮自旋軌道耦合計(jì)算能帶的過程？已有2人參與

現(xiàn)在算的體系需要考慮自旋軌道耦合（SOC）對(duì)能帶、態(tài)密度的影響，那請(qǐng)問
1）具體的步驟及INCAR、KPOINTS的設(shè)定是如何的？
2）下面是我認(rèn)為的順序：
   1、正常結(jié)構(gòu)優(yōu)化
   2、正常自洽計(jì)算得到電荷密度（1、2兩步與一般計(jì)算能帶過程一樣吧？�。�
   3、加上SOC計(jì)算能帶、態(tài)密度
3）每一步（尤其是自洽和計(jì)算能帶兩步）的KPOINTS如何選取？

我的計(jì)算經(jīng)常出現(xiàn)如下問題：
====================================================
magnetization density of overlapping atoms calculated
entering main loop
   N    E                   dE          d eps    ncg    rms       rms(c)
  WARNING: DENMP: can't reach specified precision
  Number of Electrons is NELECT = 44.9843750000000
DAV: 1 -0.607409496817+115 -0.60741+115 -0.41430+115  4680 0.151+114
  WARNING: DENMP: can't reach specified precision
  Number of Electrons is NELECT = 45.0312500000000
DAV: 2 -0.762479432673+116 -0.70174+116 -0.75190+116  7368 0.225+113
WARNING: Sub-Space-Matrix is not hermitian in DAV          6
WARNING: Sub-Space-Matrix is not hermitian in DAV          6
  9.401988633190671E+099
WARNING: Sub-Space-Matrix is not hermitian in DAV          6
  9.401988633190671E+099
WARNING: Sub-Space-Matrix is not hermitian in DAV          6
====================================================

回復(fù)此樓

» 猜你喜歡

各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩已經(jīng)有12人回復(fù)
359求調(diào)劑已經(jīng)有7人回復(fù)
070300化學(xué)專業(yè)279調(diào)劑已經(jīng)有10人回復(fù)
環(huán)境工程297分求調(diào)劑一志愿杭高院已經(jīng)有9人回復(fù)
一志愿武理材料工程302調(diào)劑環(huán)化或化工已經(jīng)有9人回復(fù)
309分085801求調(diào)劑已經(jīng)有5人回復(fù)
一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩已經(jīng)有5人回復(fù)
材料調(diào)劑已經(jīng)有4人回復(fù)
材料調(diào)劑已經(jīng)有3人回復(fù)
311（085601）求調(diào)劑已經(jīng)有12人回復(fù)

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

vasp計(jì)算能帶時(shí)考慮自旋軌道耦合的情況？已經(jīng)有7人回復(fù)
vasp計(jì)算時(shí)考慮自旋軌道耦合的設(shè)置已經(jīng)有9人回復(fù)
vasp中有關(guān)自旋已經(jīng)有31人回復(fù)
vasp計(jì)算磁性的時(shí)候 OSICAR中有總磁矩的大小，但是OUTCAR中每個(gè)原子的磁矩均顯示為0 已經(jīng)有8人回復(fù)
自旋軌道耦合后怎么算態(tài)密度已經(jīng)有4人回復(fù)
自旋軌道耦合第一性原理計(jì)算已經(jīng)有11人回復(fù)
VASP計(jì)算-考慮自旋極化已經(jīng)有3人回復(fù)
VASP如何讀取分子軌道？已經(jīng)有18人回復(fù)
vasp 中 NELECT的設(shè)置已經(jīng)有5人回復(fù)
求助vasp中的自旋軌道耦合已經(jīng)有3人回復(fù)
ADF計(jì)算自旋軌道耦合已經(jīng)有7人回復(fù)
vasp計(jì)算磁矩的困惑已經(jīng)有10人回復(fù)
【求助】vasp 離子馳豫與電子自旋的關(guān)系已經(jīng)有11人回復(fù)
【求助】VASP考慮自旋軌道耦合的話，如何編譯已經(jīng)有17人回復(fù)
【求助】用vasp計(jì)算磁矩的問題已經(jīng)有13人回復(fù)
【求助】Molpro如何執(zhí)行旋軌耦合計(jì)算如何執(zhí)行態(tài)平均計(jì)算? 已經(jīng)有10人回復(fù)

1樓 2014-11-17 17:03:38

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

whywhy01

金蟲 (小有名氣)

應(yīng)助: 35 (小學(xué)生)
金幣: 725.8
紅花: 4
帖子: 113
在線: 46小時(shí)
蟲號(hào): 1515708
注冊(cè): 2011-11-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★
感謝參與，應(yīng)助指數(shù) +1
dx0620(ljw4010代發(fā)): 金幣+2, 謝謝指導(dǎo)解答！ 2014-11-17 18:54:56

第2步的時(shí)候應(yīng)該設(shè)定NBANDS和LMAXMIX，我看后面有NELECT的提示，我猜可能是NBANDS的問題。第2步的NBANDS應(yīng)該是第1步的兩倍

贊一下

回復(fù)此樓

2樓2014-11-17 17:17:16

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

dx0620

木蟲 (正式寫手)

應(yīng)助: 14 (小學(xué)生)
金幣: 4053.6
散金: 5
紅花: 6
帖子: 393
在線: 910小時(shí)
蟲號(hào): 1303311
注冊(cè): 2011-05-22
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

2樓: Originally posted by whywhy01 at 2014-11-17 17:17:16
第2步的時(shí)候應(yīng)該設(shè)定NBANDS和LMAXMIX，我看后面有NELECT的提示，我猜可能是NBANDS的問題。第2步的NBANDS應(yīng)該是第1步的兩倍

手冊(cè)上寫的不是第3步是第2步的兩倍嗎？
第2步設(shè)定LMAXMIX了，這個(gè)就是一個(gè)數(shù)吧？如果兩種元素如何設(shè)置？

贊一下

回復(fù)此樓

3樓2014-11-17 17:27:00

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

whywhy01

金蟲 (小有名氣)

應(yīng)助: 35 (小學(xué)生)
金幣: 725.8
紅花: 4
帖子: 113
在線: 46小時(shí)
蟲號(hào): 1515708
注冊(cè): 2011-11-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by dx0620 at 2014-11-17 17:27:00
手冊(cè)上寫的不是第3步是第2步的兩倍嗎？
第2步設(shè)定LMAXMIX了，這個(gè)就是一個(gè)數(shù)吧？如果兩種元素如何設(shè)置？...

對(duì)的我記錯(cuò)了....選LMAXMIX比較大的那個(gè)數(shù)

贊一下

回復(fù)此樓

4樓2014-11-17 19:26:14

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

zh11111

銀蟲 (正式寫手)

---

應(yīng)助: 2 (幼兒園)
金幣: 338
散金: 1051
紅花: 6
帖子: 396
在線: 372.5小時(shí)
蟲號(hào): 2656941
注冊(cè): 2013-09-15
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

我在進(jìn)行第三步soc計(jì)算時(shí)INCAR為SYSTEM
PREC=High
ENCUT=520
EDIFF=1e-5
ISIF=3
ISYM=0
NSW=0
ISMEAR=-5
ISTART=1
ICHARG=11
LWAVE=TRUE
LCHARG=TRUE
NPAR=4
EDIFFG=-1e-2
ISPIN=2
LORBIT=11
LMAXMIX=4
LSORBIT=T
LORBMOM=T
SAXIS=1 0 0
NBANDS=192
NGXF=120
NGYF=54
NGZF=126
結(jié)果顯示錯(cuò)誤：LDA part: xc-table for Ceperly-Alder, standard interpolation
found WAVECAR, reading the header
  number of k-points has changed, file: 39 present: 125
  trying to continue reading WAVECAR, but it might fail
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(          1 )
reading WAVECAR
reading wavefunctions of collinear run, up
ERROR: while reading WAVECAR, plane wave coefficients changed    15319
   15285
這個(gè)問題您知道如何改進(jìn)嗎？

贊一下

回復(fù)此樓

---

5樓2015-05-26 07:37:39

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

jlmxtyy

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2273.5
散金: 1029
紅花: 2
帖子: 382
在線: 103小時(shí)
蟲號(hào): 246059
注冊(cè): 2006-04-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

5樓: Originally posted by zh11111 at 2015-05-26 07:37:39
我在進(jìn)行第三步soc計(jì)算時(shí)INCAR為SYSTEM
PREC=High
ENCUT=520
EDIFF=1e-5
ISIF=3
ISYM=0
NSW=0
ISMEAR=-5
ISTART=1
ICHARG=11
LWAVE=TRUE
LCHARG=TRUE
NPAR=4
EDIFFG=-1e-2
ISPIN=2
LORBIT=11
LMAXM ...

ISIF=3?你這個(gè)是要優(yōu)化還是干啥��？
另外算能帶的時(shí)候，可以不用之前的WAVECAR，只用CHGCAR,就不會(huì)有錯(cuò)誤提示

贊一下

回復(fù)此樓

6樓2016-06-13 16:56:38

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

wobure

新蟲 (小有名氣)

應(yīng)助: 14 (小學(xué)生)
金幣: 657.2
散金: 100
紅花: 5
帖子: 177
在線: 84.6小時(shí)
蟲號(hào): 3268517
注冊(cè): 2014-06-11
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

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

不想麻煩的話，第二步就開自旋軌道耦合計(jì)算，K點(diǎn)和第一步一樣，第三步用第二步的Charge計(jì)算能帶，k點(diǎn)line mode，自定義都可以。讀入標(biāo)準(zhǔn)dft的共線計(jì)算，個(gè)人認(rèn)為只是為了增加收斂速度而已。個(gè)人經(jīng)驗(yàn)，僅供參考！

贊一下

回復(fù)此樓

7樓2016-06-15 13:49:01

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

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

[考研] 288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行 +22	lllllos 2026-03-30	22/1100	2026-03-31 22:16 by 天倌賜福
[考研] 0856材料與化工調(diào)劑，339 +8	10213207 2026-03-31	8/400	2026-03-31 22:08 by lbsjt
[考研] 080200專業(yè)277分，求帶走！ +4	瓶子PZ 2026-03-31	4/200	2026-03-31 20:16 by vgtyfty
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-30	11/550	2026-03-31 16:34 by 邱gl
[考研] 362求調(diào)劑 +9	西南交材料專碩3 2026-03-31	9/450	2026-03-31 16:25 by w蟲蟲123
[考研] 物理學(xué)調(diào)劑 +4	小羊36 2026-03-30	4/200	2026-03-31 16:16 by lishahe
[考研] 化學(xué)0703 調(diào)劑 306分一志愿211 +10	26要上岸 2026-03-28	10/500	2026-03-31 16:04 by 記事本2026
[考研] 367求調(diào)劑 +7	芋泥啵啵… 2026-03-28	7/350	2026-03-31 14:55 by 不吃魚的貓
[考研] 323分食品與營(yíng)養(yǎng)調(diào)劑 +3	嘿ooo 2026-03-31	3/150	2026-03-31 09:38 by longlotian
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3	小圓圈圈ooo 2026-03-30	3/150	2026-03-30 21:05 by 余震yz
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5	天賦帶到THU 2026-03-29	6/300	2026-03-30 20:49 by dick_runner
[考研] 279求調(diào)劑 +12	j的立方 2026-03-29	12/600	2026-03-30 20:30 by dick_runner
[考研] 292求調(diào)劑 +13	是妍子也是研子 2026-03-30	13/650	2026-03-30 18:01 by 小徐0109
[考研] 材料專碩 085600求調(diào)劑 +7	BBQ233 2026-03-30	7/350	2026-03-30 17:44 by oooqiao
[考研] 337求調(diào)劑 +6	《樹》 2026-03-29	6/300	2026-03-30 10:15 by herarysara
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8	mikasa的圍巾 2026-03-28	8/400	2026-03-29 12:48 by 無際的草原
[考研] 343求調(diào)劑 +6	愛羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無際的草原
[考研] 356求調(diào)劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7	小小甜面團(tuán) 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry