| 3 | 1/1 | 返回列表 |
| 查看: 331 | 回復(fù): 2 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者瀟瀟的KK將贈(zèng)送您 5 個(gè)金幣 | ||
瀟瀟的KK鐵蟲 (初入文壇)
|
[求助]
QE計(jì)算電聲耦合的時(shí)候報(bào)錯(cuò)Error in routine lambda (100)wrong or too many modes 已有1人參與
|
|
|
各位老師,最近用QE用DFPT計(jì)算體系的電聲耦合的時(shí)候發(fā)現(xiàn)lambda.x時(shí)出現(xiàn)Error in routine lambda (100):wrong # or too many modes。體系的原子數(shù)一共有50個(gè),用如下的計(jì)算其他的比這個(gè)原子少的都沒(méi)有問(wèn)題。請(qǐng)各位老師指點(diǎn)指點(diǎn) 50 0.12 1 8 0.000000000 0.000000000 0.000000000 1 0.000002612 0.000000000 -0.203199324 1 0.000000000 -0.286852146 0.000000000 1 0.000002612 -0.286852146 -0.203199324 1 -0.500000000 0.000000000 -0.000005312 1 -0.499997388 0.000000000 -0.203204636 1 -0.500000000 -0.286852146 -0.000005312 1 -0.499997388 -0.286852146 -0.203204636 1 elph_dir/elph.inp_lambda.1 elph_dir/elph.inp_lambda.2 elph_dir/elph.inp_lambda.3 elph_dir/elph.inp_lambda.4 elph_dir/elph.inp_lambda.5 elph_dir/elph.inp_lambda.6 elph_dir/elph.inp_lambda.7 elph_dir/elph.inp_lambda.8 0.10 lambda.in |
木蟲 (正式寫手)
|
這個(gè) Error in routine lambda (100): wrong # or too many modes 在 QE 的 lambda.x 里是一個(gè)非常典型、而且和“體系大”強(qiáng)相關(guān)的問(wèn)題。對(duì) 50 原子體系,下面任何一條都會(huì)觸發(fā)你這個(gè)錯(cuò)誤: (A) 有 q 點(diǎn) 聲子計(jì)算沒(méi)有完全結(jié)束 ph.x 在某個(gè) q 點(diǎn):超時(shí), SCF 沒(méi)收斂, 內(nèi)存不夠被 kill, 但你仍然生成了 elph.inp_lambda.i 結(jié)果:該 q 點(diǎn)的模式數(shù) < 150 lambda.x 的 “wrong # or too many modes” 基本可以確定是:某個(gè) q 點(diǎn)的 elph.inp_lambda 文件中聲子模式數(shù) ≠ 3×Nat(你的體系是 150),通常由 ph.x 未完整結(jié)束、SCF 不收斂或混用 elph 文件導(dǎo)致。逐個(gè)檢查 elph 文件中的模式數(shù)即可定位。 |
木蟲 (正式寫手)
|
排查步驟(強(qiáng)烈建議) Step 1:檢查每個(gè) elph 文件的模式數(shù) 對(duì)每個(gè)文件: grep -i "mode" elph_dir/elph.inp_lambda.1 | wc -l 或者更粗暴一點(diǎn): grep -i "omega" elph_dir/elph.inp_lambda.1 | wc -l 必須每個(gè)都是 150 如果你看到: 147 144 0 或明顯不同, 那就是問(wèn)題源頭。 Step 2:檢查 ph.x 輸出是否完整 在對(duì)應(yīng)的 ph.out 里搜索: grep -i "Diagonalizing the dynamical matrix" ph.out 你應(yīng)該看到 每個(gè) q 點(diǎn)都完成了對(duì)角化。 如果某個(gè) q 點(diǎn):沒(méi)有這行 或后面沒(méi)有列出 150 個(gè)頻率 那個(gè) q 點(diǎn)的 elph 文件 必壞 Step 3:確認(rèn)沒(méi)有“混用 elph 文件” 檢查時(shí)間戳: ls -l elph_dir/elph.inp_lambda.* 如果發(fā)現(xiàn):某幾個(gè)文件明顯比其他舊, 或來(lái)自不同計(jì)算 刪掉全部 elph 文件,重新跑一遍 ph + elph Step 4:確認(rèn)聲子計(jì)算參數(shù) 對(duì)大體系,強(qiáng)烈建議: &inputph ldisp = .true. elph = .true. tr2_ph = 1.0d-14 epsil = .false. / 并且: 不要用過(guò)低的 tr2_ph, 確保 ecutwfc / ecutrho 足夠 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +3 | 土土小蟲 2026-03-03 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 學(xué)碩材料275調(diào)劑 +4 | 路三三 2026-03-03 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研]
材料325求調(diào)劑
30+3
|
mariusuki 2026-03-02 | 7/350 |
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考研] 哈工大計(jì)算機(jī)劉劼團(tuán)隊(duì)招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
[考研] 材料學(xué)調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
[基金申請(qǐng)]
|
Doma 2026-03-01 | 7/350 |
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|