| 5 | 1/1 | 返回列表 |
| 查看: 823 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
mouwat鐵蟲 (初入文壇)
|
[求助]
The eigenvector number1 has zero norm
|
|||
|
在計(jì)算graphene時(shí) 不知道為何會(huì)在這一部分算不過去 由於我是由單層BN的input改寫的 應(yīng)該只需要變更lattice constant吧 但卻出現(xiàn)了虛數(shù)無法計(jì)算的部分 可以請(qǐng)問是哪裡出問題了呢 para.err為 0: 0: fxphas_seq.F90:137:ERROR 0: fxphas_seq: 0: The eigenvector number1 has zero norm. tnlo_2.in檔如下 #Definition of lattice parameters #******************************** acell 4.6686384981 4.6686384981 30 rprim 1 0 0 -0.5 0.8660254038 0.0 0 0 1 #Definition of atomic positions and types #************************************* natom 2 ntypat 1 znucl 6 typat 1 1 xred 0.33333 0.66667 0.5 0.66667 0.33333 0.5 ixc 7 #Parameters of the SCF cycles #**************************** iscf 7 nstep 200 #Plane wave basis and k-point grid #********************************* ecut 50 # ecutsm 0.5 dilatmx 1.05 ngkpt 12 12 2 nshiftk 1 shiftk 0.0 0.0 0.5 ndtset 5 jdtset 1 2 3 4 5 #DATASET1 : scf calculation: GS WF in the BZ #******************************************** prtden1 1 kptopt1 1 toldfe1 1.0d-12 #DATASET2 : non scf calculation: GS WF in the whole BZ #***************************************************** getden2 1 kptopt2 2 iscf2 -2 getwfk2 1 tolwfr2 1.0d-22 nband2 4 #DATASET3 : ddk on the MP grid of k-points #***************************************** getwfk3 2 rfdir3 1 1 1 rfelfd3 2 tolwfr3 1.0d-22 nband3 4 kptopt3 2 #DATASET4 : ddE on the MP grid of k-points #***************************************** prtden4 1 getwfk4 2 getddk4 3 rfdir4 1 1 1 rfphon4 1 rfatpol4 1 2 tolvrs4 1.0d-12 rfelfd4 3 nband4 4 kptopt4 2 prepanl4 1 #DATASET5 : 3DTE calculation #*************************** getden5 1 get1den5 4 getwfk5 2 get1wf5 4 nband5 4 kptopt5 2 optdriver5 5 rf1elfd5 1 rf1phon5 1 rf1atpol5 1 2 rf1dir5 1 1 1 rf2elfd5 1 rf2dir5 1 1 1 rf3elfd5 1 rf3dir5 1 1 1 nbdbuf 0 |
鐵蟲 (初入文壇)
|
DATASET 4 計(jì)算開始會(huì)NAN 我附上了output檔 感謝你的回應(yīng) |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |
嘗試多次,附件下載不下來~LZ的DataSet 1有沒有很好的收斂?我測(cè)試了一下樓主的input,發(fā)現(xiàn)DataSet 1的時(shí)候就非常難以收斂。 仔細(xì)觀察了一下LZ的input,發(fā)現(xiàn)至少有兩處不妥。 1,對(duì)于xred的有效數(shù)字問題。為了更有效準(zhǔn)確地找到對(duì)稱性,Abinit一般要求至少10位有效數(shù)字,而LZ明顯只給了6位。比較簡單的解決方法顯然是寫為分?jǐn)?shù): xred 1/3 2/3 1/2 2/3 1/3 1/2; 2,對(duì)于Graphene類的東西,由于是模擬2D的物質(zhì),所以k點(diǎn)在z方向上為1,也即:ngkpt 12 12 1。 經(jīng)過此改動(dòng)后,DataSet 1和DataSet 2的收斂是沒有問題了。 但是DataSet 3的收斂還是有問題的(準(zhǔn)確說就是無法收斂),現(xiàn)在在嘗試解決RF的收斂問題。 如果DataSet 3無法有效收斂,在達(dá)到額定步數(shù)后,程序依然會(huì)認(rèn)為計(jì)算完成,并把結(jié)果用于下一步的計(jì)算,自然在DataSet 4會(huì)出現(xiàn)NaN~ |

鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +7 | Koxui 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 5/250 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 打過很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|