| 查看: 2973 | 回复: 1 | ||
andy101199银虫 (正式写手)
|
[求助]
QE计算声电耦合最后lambda.x出现问题
|
|
QE计算声电耦合在进行最后一步ambda.x出现问题，之前步骤都正常结束。 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 错误提示 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文件，但是感觉明显是错误的 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 ********* 请教各位大神，这个问题如何解决？ |
回复此楼» 猜你喜欢
2026年WR青拔进展
已经有5人回复
-
国自然申请五篇代表作大比拼，感觉这个是最重要的
已经有10人回复
-
无聊看看时间戳打发时间
已经有4人回复
-
中！中！中！
已经有9人回复
-
咨询
已经有3人回复
-
职称论文投稿
已经有9人回复
-
基于自然哲学类比的风化壳型稀土矿
已经有14人回复
-
评委有多少概率知道其他专家手中有哪些人的本子？
已经有6人回复
-
求推荐期刊，重谢
已经有3人回复
-
青A35岁以下通知答辩了吗
已经有4人回复
2楼2018-12-19 22:32:34