| 6 | 1/1 | 返回列表 |
| 查看: 822 | 回復(fù): 5 | ||||
mouwat鐵蟲(chóng) (初入文壇)
|
[求助]
The eigenvector number1 has zero norm
|
|
在計(jì)算graphene時(shí) 不知道為何會(huì)在這一部分算不過(guò)去 由於我是由單層BN的input改寫(xiě)的 應(yīng)該只需要變更lattice constant吧 但卻出現(xiàn)了虛數(shù)無(wú)法計(jì)算的部分 可以請(qǐng)問(wèn)是哪裡出問(wèn)題了呢 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 |
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +174 |

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

鐵蟲(chóng) (初入文壇)
鐵蟲(chóng) (初入文壇)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +8 | 15188958825 2026-03-25 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 085600,專(zhuān)業(yè)課化工原理,321分求調(diào)劑 +4 | 大饞小子 2026-03-28 | 4/200 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
[考研] 320分,材料與化工專(zhuān)業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 329求調(diào)劑 +6 | 星野? 2026-03-26 | 6/300 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研]
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|