 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 731 | 回復(fù): 2 | |||
燦爛的幸福新蟲 (小有名氣)
|
[求助]
超晶胞電聲耦合計算的問題
|
|
我用MS構(gòu)建了一個超晶胞(已經(jīng)運行過例子，結(jié)果正確)，然后準(zhǔn)備計算它的電聲耦合，但是運行第四步時出現(xiàn)如下錯誤： 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)價值貼推薦，對您同樣有幫助:
- QE軟件計算電聲子耦合，最后一步T_c出現(xiàn)NaN，什么原因呢？
已經(jīng)有6人回復(fù)
- 怎么補(bǔ)算電聲子耦合？
已經(jīng)有3人回復(fù)
- 電聲耦合計算中l(wèi)ambd.x的一些問題
已經(jīng)有10人回復(fù)
- MS 能計算超導(dǎo)嗎
已經(jīng)有6人回復(fù)
燦爛的幸福
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時
- 蟲號: 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 請各路大神指點迷津，謝謝！ |
2樓2014-04-14 19:30:30
燦爛的幸福
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 906.3
- 帖子: 93
- 在線: 101.5小時
- 蟲號: 1335565
- 注冊: 2011-07-01
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
不要沉，自己頂下。這個問題很難嗎？ 3樓2014-04-21 10:12:48
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 燦爛的幸福 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 張zz111 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
|
[考研] 283求調(diào)劑（080500） +4 | A child 2026-03-27 | 4/200 |
|
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +3 | Zuhui0306 2026-03-25 | 3/150 |
|
|
|
[考研] 317求調(diào)劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
|
[考研] 085600材料與化工306 +10 | z1z2z3879 2026-03-21 | 11/550 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +5 | Sxy112 2026-03-21 | 7/350 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|