| 5 | 1/1 | 返回列表 |
| 查看: 824 | 回復: 5 | |||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||
[求助]
The eigenvector number1 has zero norm
|
|||||
|
在計算graphene時 不知道為何會在這一部分算不過去 由於我是由單層BN的input改寫的 應該只需要變更lattice constant吧 但卻出現了虛數無法計算的部分 可以請問是哪裡出問題了呢 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 |
專家顧問 (著名寫手)
![]() |
專家經驗: +174 |
嘗試多次,附件下載不下來~LZ的DataSet 1有沒有很好的收斂?我測試了一下樓主的input,發(fā)現DataSet 1的時候就非常難以收斂。 仔細觀察了一下LZ的input,發(fā)現至少有兩處不妥。 1,對于xred的有效數字問題。為了更有效準確地找到對稱性,Abinit一般要求至少10位有效數字,而LZ明顯只給了6位。比較簡單的解決方法顯然是寫為分數: xred 1/3 2/3 1/2 2/3 1/3 1/2; 2,對于Graphene類的東西,由于是模擬2D的物質,所以k點在z方向上為1,也即:ngkpt 12 12 1。 經過此改動后,DataSet 1和DataSet 2的收斂是沒有問題了。 但是DataSet 3的收斂還是有問題的(準確說就是無法收斂),現在在嘗試解決RF的收斂問題。 如果DataSet 3無法有效收斂,在達到額定步數后,程序依然會認為計算完成,并把結果用于下一步的計算,自然在DataSet 4會出現NaN~ |

專家顧問 (著名寫手)
![]() |
專家經驗: +174 |

|
DATASET 4 計算開始會NAN 我附上了output檔 感謝你的回應 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 總分293求調劑 +6 | 加一一九 2026-03-25 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經歷 +6 | cccchenso 2026-03-29 | 6/300 |
|
|
[考研] 一志愿鄭州大學,080500學碩,總分317分求調劑 +8 | 舉個栗子oi 2026-03-24 | 9/450 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 312,生物學求調劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 復試調劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 291求調劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 308求調劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 一志愿華東理工大學081700,初試分數271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 286求調劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 材料調劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 調劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 中國科學院深圳先進技術研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 考研一志愿蘇州大學初始315(英一)求調劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 086003食品工程求調劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 300分,材料,求調劑,英一數二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 一志愿山東大學藥學學碩求調劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 336化工調劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研] 石河子大學(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|