| 3 | 1/1 | 返回列表 |
| 查看: 367 | 回復(fù): 2 | |||
| 【懸賞金幣】回答本帖問題,作者瀟瀟的KK將贈送您 5 個金幣 | |||
瀟瀟的KK鐵蟲 (初入文壇)
|
[求助]
QE計算電聲耦合的時候報錯Error in routine lambda (100)wrong or too many modes 已有1人參與
|
||
|
各位老師,最近用QE用DFPT計算體系的電聲耦合的時候發(fā)現(xiàn)lambda.x時出現(xiàn)Error in routine lambda (100):wrong # or too many modes。體系的原子數(shù)一共有50個,用如下的計算其他的比這個原子少的都沒有問題。請各位老師指點(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 |
木蟲 (正式寫手)
|
這個 Error in routine lambda (100): wrong # or too many modes 在 QE 的 lambda.x 里是一個非常典型、而且和“體系大”強(qiáng)相關(guān)的問題。對 50 原子體系,下面任何一條都會觸發(fā)你這個錯誤: (A) 有 q 點(diǎn) 聲子計算沒有完全結(jié)束 ph.x 在某個 q 點(diǎn):超時, SCF 沒收斂, 內(nèi)存不夠被 kill, 但你仍然生成了 elph.inp_lambda.i 結(jié)果:該 q 點(diǎn)的模式數(shù) < 150 lambda.x 的 “wrong # or too many modes” 基本可以確定是:某個 q 點(diǎn)的 elph.inp_lambda 文件中聲子模式數(shù) ≠ 3×Nat(你的體系是 150),通常由 ph.x 未完整結(jié)束、SCF 不收斂或混用 elph 文件導(dǎo)致。逐個檢查 elph 文件中的模式數(shù)即可定位。 |
木蟲 (正式寫手)
|
排查步驟(強(qiáng)烈建議) Step 1:檢查每個 elph 文件的模式數(shù) 對每個文件: grep -i "mode" elph_dir/elph.inp_lambda.1 | wc -l 或者更粗暴一點(diǎn): grep -i "omega" elph_dir/elph.inp_lambda.1 | wc -l 必須每個都是 150 如果你看到: 147 144 0 或明顯不同, 那就是問題源頭。 Step 2:檢查 ph.x 輸出是否完整 在對應(yīng)的 ph.out 里搜索: grep -i "Diagonalizing the dynamical matrix" ph.out 你應(yīng)該看到 每個 q 點(diǎn)都完成了對角化。 如果某個 q 點(diǎn):沒有這行 或后面沒有列出 150 個頻率 那個 q 點(diǎn)的 elph 文件 必壞 Step 3:確認(rèn)沒有“混用 elph 文件” 檢查時間戳: ls -l elph_dir/elph.inp_lambda.* 如果發(fā)現(xiàn):某幾個文件明顯比其他舊, 或來自不同計算 刪掉全部 elph 文件,重新跑一遍 ph + elph Step 4:確認(rèn)聲子計算參數(shù) 對大體系,強(qiáng)烈建議: &inputph ldisp = .true. elph = .true. tr2_ph = 1.0d-14 epsil = .false. / 并且: 不要用過低的 tr2_ph, 確保 ecutwfc / ecutrho 足夠 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 6/300 |
|
|---|---|---|---|---|
|
[基金申請] 提交后的基金本子,已讓學(xué)校撤回了,可否換口子提交 +3 | dut_pfx 2026-03-10 | 3/150 |
|
|
[考研] 312求調(diào)劑 +5 | 陌宸希 2026-03-10 | 5/250 |
|
|
[考研] 307求調(diào)劑 +6 | 超級伊昂大王 2026-03-10 | 6/300 |
|
|
[考研] 調(diào)劑 +5 | 呵唔哦豁 2026-03-10 | 5/250 |
|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 8/400 |
|
|
[考研] 求調(diào)劑! +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
[考研] 一志愿湖師大化學(xué)289求調(diào)劑 +5 | XMCMM3.14159 2026-03-10 | 5/250 |
|
|
[考研] 278求調(diào)劑 +8 | Gale1314 2026-03-06 | 8/400 |
|
|
[考研] 物理學(xué)求調(diào)劑 +4 | sx宋 2026-03-05 | 4/200 |
|
|
[考研] 085600材料與化工,一志愿廣州985,求調(diào)劑 +15 | qqyyaill 2026-03-05 | 15/750 |
|
|
[考研] 一志愿山東大學(xué),總分327,英語二79,有論文,有競賽,已過四六級 +3 | 木木目目1 2026-03-09 | 3/150 |
|
|
[考研] 085701環(huán)境工程專碩求調(diào)劑 +5 | 幼稚園大班打手 2026-03-08 | 5/250 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
[考研] 276求調(diào)劑 本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4 | 路lyh123 2026-03-06 | 6/300 |
|
|
[考研] 求調(diào)劑 +4 | 呼呼?~+123456 2026-03-06 | 4/200 |
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
[考研] 328求調(diào)劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|