 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 727 | 回復(fù): 2 | ||
燦爛的幸福新蟲(chóng) (小有名氣)
|
[求助]
超晶胞電聲耦合計(jì)算的問(wèn)題
|
|
我用MS構(gòu)建了一個(gè)超晶胞(已經(jīng)運(yùn)行過(guò)例子，結(jié)果正確)，然后準(zhǔn)備計(jì)算它的電聲耦合，但是運(yùn)行第四步時(shí)出現(xiàn)如下錯(cuò)誤： symdm9.F90:213:ERROR Informations are missing in the DDB. The dynamical matrix number 3 cannot be built. since no blok with wavevector 0.00 0.50 0.00 has been found. Action: add the required blok in the DDB , or modify your input file. 還請(qǐng)各位不吝賜教。此處貼上我的第一步輸入文件telphon_1.in： ndtset 5 # # DATASET 1 : make ground state wavefunctions and density # tolwfr1 1.0d-20 nline1 8 # This is to expedite the convergence of higher-lying bands rfphon1 0 # for DS1 do _not_ do perturbation nqpt1 0 # for DS1 do _not_ do perturbation prtwf1 1 # need GS wavefunctions for further runs getwfk1 0 # enforce calculation of forces at each SCF step optforces 1 qpt2 0.0 0.0 0.0 qpt3 0.5 0.0 0.0 qpt4 0.5 0.5 0.0 # # DS5 get DDK perturbation (stored in GKK matrix element files) # tolwfr5 1.0d-14 qpt5 0 0 0 rfphon5 0 rfelfd5 2 prtwf5 0 # # general data for all phonon calculations: # rfatpol 1 8 rfdir 1 1 1 rfphon 1 prtwf 0 # for response function runs, do not keep first order wavefunctions tolvrs 1.0e-5 getwfk 1 nqpt 1 prepgkk 1 # force all perturbations to be calculated for q-points considered prtgkk 1 # print out GKK files containing electron-phonon coupling # # Common data for GS and perturbation runs # # # the kpoint grid is minimalistic to keep the calculation # manageable. # ngkpt 12 12 12 kptopt 3 # keep all k-points. Also for ground state, because of eventual transport calculations # # use a centered grid for the kpoints: obligatory for el-ph for the moment # nshiftk 1 shiftk 0.0 0.0 0.0 # # as is the kinetic energy cutoff # ecut 30 acell 3*13.472 rprim 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 nband 80 # # include metallic occupation function with a small smearing # occopt 7 tsmear 0.001 natom 8 typat 1 1 1 1 1 1 1 1 xred 0.00 0.00 0.00 0.50 0.00 0.00 0.00 0.50 0.00 0.50 0.50 0.00 0.00 0.00 0.50 0.50 0.00 0.50 0.00 0.50 0.50 0.50 0.50 0.50 nstep 800 ntypat 1 znucl 29 |
回復(fù)此樓» 猜你喜歡
285求調(diào)劑
已經(jīng)有6人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
291求調(diào)劑
已經(jīng)有6人回復(fù)
-
315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
277跪求調(diào)劑
已經(jīng)有8人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
295材料工程專(zhuān)碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
083000學(xué)碩274求調(diào)劑
已經(jīng)有4人回復(fù)
-
275求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- QE軟件計(jì)算電聲子耦合，最后一步T_c出現(xiàn)NaN，什么原因呢？
已經(jīng)有6人回復(fù)
- 怎么補(bǔ)算電聲子耦合？
已經(jīng)有3人回復(fù)
- 電聲耦合計(jì)算中l(wèi)ambd.x的一些問(wèn)題
已經(jīng)有10人回復(fù)
- MS 能計(jì)算超導(dǎo)嗎
已經(jīng)有6人回復(fù)
燦爛的幸福
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時(shí)
- 蟲(chóng)號(hào): 1335565
- 注冊(cè): 2011-07-01
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
再補(bǔ)充第四步的輸入文件telphon_4.in： # turn on calculation of the electron-phonon quantities elphflag 1 # Path in reciprocal space along which the phonon linewidths # and band structure will be calculated nqpath 7 qpath 0.0 0.0 0.0 1/2 1/2 0.0 1 1 1 1/2 1/2 1/2 1/2 1/2 0.0 1/2 3/4 1/4 1/2 1/2 1/2 # Coulomb pseudopotential parameter mustar 0.136 # Minimalistic qpoint grid ngqpt 2 2 2 # impose acoustic sum rule in a symmetric way asr 2 dipdip 0 # bravais lattice necessary brav 2 # qpt grid nqshft 1 q1shft 0.0 0.0 0.0 # ifcflag 1 ifcana 1 # ifc for all atoms? natifc 0 atifc 1 2 3 # print dielectric matrix with freq dependence dieflag 0 # print out eigenvectors and symmetrize dyn matrix eivec 1 #Wavevector list number 1 (Reduced coordinates and normalization factor) nph1l 1 qph1l 0.00000000E+00 0.00000000E+00 0.00000000E+00 1 # This line added when defaults were changed (v5.3) to keep the previous, old behaviour symdynmat 0 請(qǐng)各路大神指點(diǎn)迷津，謝謝！ |
2樓2014-04-14 19:30:30
燦爛的幸福
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時(shí)
- 蟲(chóng)號(hào): 1335565
- 注冊(cè): 2011-07-01
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
不要沉，自己頂下。這個(gè)問(wèn)題很難嗎？ 3樓2014-04-21 10:12:48
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 燦爛的幸福 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 315分求調(diào)劑 +6 | 26考研上岸版26 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 295材料工程專(zhuān)碩求調(diào)劑 +3 | 1428151015 2026-03-27 | 3/150 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 308求調(diào)劑 +6 | 墨墨漠 2026-03-27 | 6/300 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書(shū)讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 材料專(zhuān)碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +4 | 要好好無(wú)聊 2026-03-21 | 4/200 |
|