zhang668

木蟲 (正式寫手)

[交流] 【問題已解決】vasp計算自旋軌道耦合讀問題

如題，先進行關(guān)閉自旋軌道耦合的共線的scf計算后，得到CHGCAR和WAVCAR，然后設(shè)置ICHARGE=11，計算非自洽時候，總出現(xiàn)錯誤，提示為：
running on 64 nodes
distr:  one band on 1 nodes, 64 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  O  Fe Bi
POSCAR found :  3 types and 10 ions
LDA part: xc-table for Ceperly-Alder, standard interpolation
found WAVECAR, reading the header
  number of k-points has changed, file: 19 present: 125
  trying to continue reading WAVECAR, but it might fail
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(          1 )
reading WAVECAR
reading wavefunctions of collinear run, up
reading wavefunctions of collinear run, down
the WAVECAR file was read sucessfully
charge-density read from file: BiFeO3
magnetization density read from file 1
LAPACK: Routine ZPOTRF failed!          1
LAPACK: Routine ZPOTRF failed!          1
LAPACK: Routine ZPOTRF failed!          1
LAPACK: Routine ZPOTRF failed!          1

就是K點數(shù)變化，出現(xiàn)問題。不管第一步有沒設(shè)置ISYM=0,第二步的計算的K點數(shù)都和第一步不一樣，所以總是出現(xiàn)問題。只有第二步計算不讀取波函數(shù)和電荷密度，計算才可以進行下去。

[ Last edited by zhang668 on 2010-12-2 at 22:21 ]

回復(fù)此樓

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

vasp計算資料

有機

自旋軌道耦合

» 猜你喜歡

283分材料與化工求調(diào)劑已經(jīng)有12人回復(fù)
349求調(diào)劑已經(jīng)有7人回復(fù)
296材料專碩求調(diào)劑已經(jīng)有11人回復(fù)
一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè) 已經(jīng)有7人回復(fù)
材料340分調(diào)劑已經(jīng)有8人回復(fù)
316求調(diào)劑已經(jīng)有13人回復(fù)
321求調(diào)劑已經(jīng)有11人回復(fù)
315求調(diào)劑已經(jīng)有6人回復(fù)
330分求調(diào)劑已經(jīng)有13人回復(fù)
化工調(diào)劑303分，過四級已經(jīng)有9人回復(fù)

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

關(guān)于VASP中MAGMON的磁矩設(shè)置問題已經(jīng)有6人回復(fù)
vasp計算磁性的時候 OSICAR中有總磁矩的大小，但是OUTCAR中每個原子的磁矩均顯示為0 已經(jīng)有8人回復(fù)
怎么補算電聲子耦合？已經(jīng)有3人回復(fù)
VASP自旋軌道耦合已經(jīng)有26人回復(fù)
NBANDS 的設(shè)置問題已經(jīng)有6人回復(fù)
自旋軌道耦合第一性原理計算已經(jīng)有11人回復(fù)
VASP 運行報錯已經(jīng)有11人回復(fù)
vasp自旋極化計算時，初始設(shè)置自旋能否不被優(yōu)化已經(jīng)有6人回復(fù)
VASP如何讀取分子軌道？已經(jīng)有18人回復(fù)
VASP+SOC后的能帶計算已經(jīng)有21人回復(fù)
求助vasp中的自旋軌道耦合已經(jīng)有3人回復(fù)
ADF計算自旋軌道耦合已經(jīng)有7人回復(fù)
VASP計算：如何處理波函數(shù)占用空間太大的問題已經(jīng)有5人回復(fù)
請問MS里怎么考慮自旋軌道耦合作用？已經(jīng)有9人回復(fù)
【求助】LAPACK: Routine ZPOTRF failed! 已經(jīng)有5人回復(fù)
【求助】vasp 離子馳豫與電子自旋的關(guān)系已經(jīng)有11人回復(fù)
【求助】VASP考慮自旋軌道耦合的話，如何編譯已經(jīng)有17人回復(fù)
【求助成功】charge-density read from file:unknown 是怎么回事呢? 已經(jīng)有3人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-12-01 23:20:53

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

byin

銅蟲 (正式寫手)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時
蟲號: 948127

★
zhang668(金幣+1):謝謝參與

你好，我也遇到了相類似的問題，提示K點不一致。請問應(yīng)該如何解決？

謝謝！

贊一下(1人)

回復(fù)此樓

2樓2011-07-19 11:34:33

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

deligent

禁蟲 (初入文壇)

★
zhang668(金幣+1):謝謝參與

本帖內(nèi)容被屏蔽

3樓2011-07-19 15:50:30

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

byin

銅蟲 (正式寫手)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時
蟲號: 948127

引用回帖:

Originally posted by deligent at 2011-07-19 15:50:30:
可能是NBANDS設(shè)置的問題，NBANDS = 2 * number of bands of collinear run。

設(shè)置了，但是還不行。

贊一下

回復(fù)此樓

4樓2011-07-24 22:17:19

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

whw19850730

鐵桿木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 5875.8
帖子: 64
在線: 121.3小時
蟲號: 1048864

★ ★
zhang668(金幣+1):謝謝參與
bingmou(金幣+1): 感謝參與 2011-07-31 17:55:16

引用回帖:

Originally posted by byin at 2011-07-24 22:17:19:
設(shè)置了，但是還不行。

第一步scf計算也開啟自詡軌道耦合和非共線，這樣產(chǎn)生的電荷密度讀入就沒錯了，說明書說的有問題。

贊一下(2人)

回復(fù)此樓

5樓2011-07-30 23:43:00

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

byin

銅蟲 (正式寫手)

1ST強帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 245
帖子: 395
在線: 174.4小時
蟲號: 948127

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

5樓: Originally posted by whw19850730 at 2011-07-30 23:43:00:
第一步scf計算也開啟自詡軌道耦合和非共線，這樣產(chǎn)生的電荷密度讀入就沒錯了，說明書說的有問題。

你的意思是直接一步計算，對嗎？

istart=0
icharge=2

自旋軌道耦合開

直接這么算？

贊一下(1人)

回復(fù)此樓

6樓2011-08-01 00:20:05

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

whw19850730

鐵桿木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 5875.8
帖子: 64
在線: 121.3小時
蟲號: 1048864

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

6樓: Originally posted by byin at 2011-08-01 00:20:05:
你的意思是直接一步計算，對嗎？

istart=0
icharge=2

自旋軌道耦合開

直接這么算？

是的讀入非共線的電荷密度就不會有錯了。你可以測試一下

贊一下(1人)

回復(fù)此樓

7樓2011-08-01 10:34:05

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

sunyb83

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 130.4
帖子: 14
在線: 2.9小時
蟲號: 468982

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

ISYM=0
關(guān)閉對稱性操作。當(dāng)啟用自旋軌道耦合時，不同磁矩指向情況下體系的對稱性存在差異，因此k點也會相應(yīng)變化。
我用的辦法是從共線計算開始就關(guān)掉對稱性操作，然后用IBZKPT 代替KPOINTS，保證各個方向計算時使用相同的K點。

贊一下(1人)

回復(fù)此樓

8樓2011-12-06 17:23:47

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

lastfreedom

木蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 4209.4
帖子: 99
在線: 658.9小時
蟲號: 1087212

★
小木蟲(金幣+0.5):給個紅包，謝謝回帖

引用回帖:

8樓: Originally posted by sunyb83 at 2011-12-06 17:23:47:
ISYM=0
關(guān)閉對稱性操作。當(dāng)啟用自旋軌道耦合時，不同磁矩指向情況下體系的對稱性存在差異，因此k點也會相應(yīng)變化。
我用的辦法是從共線計算開始就關(guān)掉對稱性操作，然后用IBZKPT 代替KPOINTS，保證各個方向計算時 ...

請教：這樣計算和relax的時候就加soc，費的時間是不是會少點，對大體系呢

贊一下(1人)

回復(fù)此樓

9樓2011-12-06 17:42:07

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

ffzx

金蟲 (正式寫手)

應(yīng)助: 53 (初中生)
金幣: 1266.9
帖子: 455
在線: 158.9小時
蟲號: 1074368

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

1943759樓: Originally posted by lastfreedom at 2011-12-06 17:42:07
請教：這樣計算和relax的時候就加soc，費的時間是不是會少點，對大體系呢...

大體系本來就慢，直接加SOC relax豈不更慢？

贊一下

回復(fù)此樓

10樓2012-06-02 20:58:16

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

國花1989

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 768.8
帖子: 109
在線: 56.8小時
蟲號: 2001124

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

樓主你好，你的問題解決了嗎？？我的自旋軌道耦合沒有加進去，是不是參數(shù)設(shè)置錯啦？？

贊一下

回復(fù)此樓

11樓2013-06-09 08:37:13

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

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

[考研] 318求調(diào)劑，計算材料方向 +10	吸喵有害笙命 2026-04-01	11/550	2026-04-02 16:29 by oooqiao
[考研] 農(nóng)學(xué)考研求調(diào)劑 +3	dkdkxm 2026-04-01	3/150	2026-04-02 16:04 by wangjagri
[考研] 318求調(diào)劑 +3	篤行致遠. 2026-03-31	4/200	2026-04-02 15:56 by Jaylen.
[考研] 285求調(diào)劑 +14	AZMK 2026-04-02	14/700	2026-04-02 15:54 by 上九天攬月（好?/a>
[考研] 302求調(diào)劑一志愿華中師范大學(xué) +5	小江小江江江 2026-04-02	5/250	2026-04-02 15:47 by zikuzi1853
[考研] 288求調(diào)劑一志愿哈工大材料與化工 +30	洛神哥哥 2026-03-31	32/1600	2026-04-02 14:49 by 賴春艷
[考研] 環(huán)境285分，過六級，求調(diào)劑 +6	xhr12 2026-04-02	6/300	2026-04-02 13:38 by sun200805022
[考研] 329求調(diào)劑 +8	miaodesi 2026-04-02	10/500	2026-04-02 10:24 by 不吃魚的貓
[考研] 322求調(diào)劑 +5	熹僖XX 2026-03-31	6/300	2026-04-02 10:08 by 求調(diào)劑zz
[考研] 279求調(diào)劑 +6	學(xué)而思兮知 2026-04-01	6/300	2026-04-02 09:16 by vgtyfty
[考研] 一志愿北京科技，085601總分305求調(diào)劑 +9	半生瓜！ 2026-04-01	11/550	2026-04-02 08:28 by Wang200018
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +10	腳滑的守法公民 2026-03-27	10/500	2026-04-01 17:23 by pies112
[考研] 0856，材料與化工321分求調(diào)劑 +13	大饞小子 2026-03-27	14/700	2026-04-01 15:58 by yanflower7133
[考研] 土木304求調(diào)劑 +5	頂級擦擦 2026-03-31	5/250	2026-04-01 08:15 by fdcxdystjk￥
[考研] 299求調(diào)劑 +8	嗯嗯嗯嗯2 2026-03-27	8/400	2026-03-31 18:23 by lizhi8172
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。 +10	XX小鄧 2026-03-29	10/500	2026-03-31 16:46 by 不吃魚的貓
[考研] 英一數(shù)一總分334求調(diào)劑 +4	陳陽坤 2026-03-31	4/200	2026-03-31 14:22 by 記事本2026
[考研] 一志愿大連理工大學(xué)，機械工程學(xué)碩，341 +3	西瓜田的守望者 2026-03-30	3/150	2026-03-31 11:08 by asdfzly
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3	小圓圈圈ooo 2026-03-30	3/150	2026-03-30 21:05 by 余震yz
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs