版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

andy101199

銀蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 7296
散金: 1160
紅花: 7
帖子: 433
在線: 456.6小時(shí)
蟲號(hào): 1326650
注冊(cè): 2011-06-19
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] QE計(jì)算聲電耦合最后lambda.x出現(xiàn)問題

QE計(jì)算聲電耦合在進(jìn)行最后一步ambda.x出現(xiàn)問題，之前步驟都正常結(jié)束。

lambda.in如下（ Program PWSCF v.5.4.0）：
75 0.02 1 ! emax (something more than highest phonon mode in THz), degauss, smearing method
8          ! Number of q-points
0.000000000 0.000000000 0.000000000  1.00
0.000000000 0.000000000  -0.136976104  1.00
0.000000000 0.288675135 0.000000000  6.00
0.000000000 0.288675135  -0.136976104  6.00
0.000000000  -0.577350269 0.000000000  3.00
0.000000000  -0.577350269  -0.136976104  3.00
0.250000000 0.433012702 0.000000000  6.00
0.250000000 0.433012702  -0.136976104  6.00
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

錯(cuò)誤提示
At line 60 of file lambda.f90 (unit = 5, file = 'stdin')
Fortran runtime error: End of file

查看了lambda.f90中第60行附近代碼：

58    real*8 mustar, omegalog(20), Tc, x

60    read(5,*) emax, degaussq, ngaussq
61 deltae=emax/(nex-1)
62    read(5,*) nks
63       if (nks.gt.npk) call errore('lambda','too many q-points',npk)
64    sum=0.d0
65    do ik=1,nks
  66    read(5,*) q(1,ik), q(2,ik), q(3,ik), wk(ik)
  67    sum = sum + wk(ik)
  68    end do
  69    do ik=1,nks
  70    wk(ik)=wk(ik)/sum
  71       end do

且能輸出lambda.out文件，但是感覺明顯是錯(cuò)誤的
lambda       omega_log       T_c
  -0.15572    2477.084       73200.645
  -0.14600    2490.699       88108.999
  -0.14491    2488.983       89922.337
  -0.14131    2486.484       96426.060
  -0.13638    2486.396       *********
  -0.13187    2491.111       *********
  -0.12821    2495.033       *********
  -0.12502    2501.643       *********
  -0.12261    2504.104       *********
  -0.12064    2512.136       *********
  -0.11938    2516.256       *********
  -0.11828    2520.488       *********
  -0.11777    2523.697       *********
  -0.11731    2528.947       *********
  -0.11708    2531.623       *********
  -0.11712    2534.020       *********
  -0.11707    2536.620       *********
  -0.11729    2538.828       *********
  -0.11730    2541.920       *********
  -0.11756    2543.034       *********
  -0.11773    2545.929       *********
  -0.11785    2545.691       *********
  -0.11799    2546.964       *********
  -0.11818    2550.281       *********
  -0.11844    2550.324       *********
  -0.11855    2551.999       *********
  -0.11878    2551.735       *********
  -0.11905    2553.195       *********
  -0.11936    2554.038       *********
  -0.11976    2556.991       *********

請(qǐng)教各位大神，這個(gè)問題如何解決？

回復(fù)此樓

» 猜你喜歡

271材料工程求調(diào)劑已經(jīng)有6人回復(fù)
材料專碩274一志愿陜西師范大學(xué)求調(diào)劑已經(jīng)有8人回復(fù)
288求調(diào)劑已經(jīng)有9人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有4人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有10人回復(fù)
321求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有6人回復(fù)
266求調(diào)劑已經(jīng)有10人回復(fù)
本人考085602 化學(xué)工程專碩已經(jīng)有19人回復(fù)

1樓 2016-08-11 10:01:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

NickyHCQ

木蟲 (正式寫手)

應(yīng)助: 12 (小學(xué)生)
金幣: 2004
散金: 55
紅花: 8
帖子: 642
在線: 754小時(shí)
蟲號(hào): 4546266
注冊(cè): 2016-03-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

不知樓主解決了沒，我也遇到了這個(gè)問題，但是out是空的

贊一下

回復(fù)此樓

2樓2018-12-19 22:32:34