 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 732 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
燦爛的幸福新蟲 (小有名氣)
|
[求助]
超晶胞電聲耦合計(jì)算的問題
|
|
|
我用MS構(gòu)建了一個(gè)超晶胞(已經(jīng)運(yùn)行過例子，結(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. 還請各位不吝賜教。此處貼上我的第一步輸入文件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ù)此樓» 猜你喜歡
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有5人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有9人回復(fù)
-
085602 化工專碩 338分 求調(diào)劑
已經(jīng)有9人回復(fù)
-
352分 化工與材料
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有3人回復(fù)
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有4人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
275求調(diào)劑
已經(jīng)有10人回復(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的一些問題
已經(jīng)有10人回復(fù)
- MS 能計(jì)算超導(dǎo)嗎
已經(jīng)有6人回復(fù)
燦爛的幸福
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時(shí)
- 蟲號(hào): 1335565
- 注冊: 2011-07-01
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
不要沉，自己頂下。這個(gè)問題很難嗎？ 
3樓2014-04-21 10:12:48
燦爛的幸福
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時(shí)
- 蟲號(hào): 1335565
- 注冊: 2011-07-01
- 性別: GG
- 專業(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 請各路大神指點(diǎn)迷津，謝謝！ |
2樓2014-04-14 19:30:30
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +5 | 辰沐5211314 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 286求調(diào)劑 +8 | PolarBear11 2026-03-26 | 8/400 |
|
|
|
[考研]
 一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
+7
|
18419759900 2026-03-25 | 8/400 |
|
|
|
[考博] 26申博 +3 | 加油沖啊！ 2026-03-26 | 3/150 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級(jí)已過，有科研經(jīng)歷，求調(diào)劑！ +3 | Zuhui0306 2026-03-25 | 3/150 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研]
085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
|
[考研] 調(diào)劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|