应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » error when running pw.x
81/1返回列表
查看: 1868  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

provels

金虫 (正式写手)

[求助] error when running pw.x

Subspace diagonalization in iterative solution of the eigenvalue problem:
     a serial algorithm will be used


     Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min           0       0      0                    0        0       0
     Max           1       1      1                    5        5       5
     Sum           1       1      1                    5        5       5


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine allocate_fft (1):
     wrong ngm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

不知道问题在哪  错误不明确  
QE 新手 求教

在LOG FILE 中，
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 4
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 1
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 3
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 2
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 6
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 7
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 5
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 1


多谢！！！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-08-16 05:46:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
编译完后测试了没有？
把输入文件贴上来看看。
赞一下 回复此楼
2楼2013-08-16 14:38:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

provels

金虫 (正式写手)
编译没问题

# self-consistent calculation
#
cat > alas.scf.in << EOF
&control
    calculation='scf'
    restart_mode='from_scratch',
    tstress = .true.
    tprnfor = .true.
    prefix='alas',
    pseudo_dir = '/home/x_weisu/noback/weiwei/yunguo/raman/mos',
    outdir='/home/x_weisu/noback/weiwei/yunguo/raman/mos'
/
&system
    nosym = .false.,
    ibrav=  4, celldm(1)=0.23372 ,celldm(3) =5.03366, nat=  3, ntyp= 2,
    ecutwfc = 37
    occupations='smearing',
    smearing='mp'
    degauss=0.00367
/
&electrons
    conv_thr =  1.0d-8
    mixing_beta = 0.7
/
ATOMIC_SPECIES
  S   32.06  S.pbe-n-rrkjus_psl.0.1.UPF
  Mo  95.96  Mo.pbe-spn-rrkjus_psl.0.2.UPF
ATOMIC_POSITIONS
S  0.9999999933631344  0.0000000035897103  0.4023766108613526
S  0.9999999964301295  0.0000000066521953  0.5976233891367776
Mo 0.6666666768734000  0.3333333230914300  0.5000000000018697
K_POINTS {automatic}
6 6 1 0 0 0
EOF

mpprun /software/testing/espresso/5.0.2/build01/bin/pw.x < alas.scf.in > alas.scf.out
赞一下 回复此楼
3楼2013-08-16 18:14:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
原来是要计算二硫化钼，那就好办多了。

这一行：
“ ibrav=  4, celldm(1)=0.23372 ,celldm(3) =5.03366, nat=  3, ntyp= 2,”

二硫化钼晶格长度大概是3.16 angstrom，换算成bohr大概是5.97，你的celldm(1)错了。
赞一下 回复此楼
4楼2013-08-16 18:59:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
provels: 金币+5 2013-08-16 20:54:10
还有“ATOMIC_POSITIONS”改成"ATOMIC_POSITIONS {crystal}"。
ATOMIC_POSITIONS后面是有选项的，默认的是alat，这个选项对应直角坐标而不是分数坐标。

算二硫化钼要加展宽？我记得不加展宽也能收敛。
赞一下(1人) 回复此楼
5楼2013-08-16 19:03:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

provels

金虫 (正式写手)
引用回帖:
4楼: Originally posted by liqizuiyang at 2013-08-16 18:59:16
原来是要计算二硫化钼，那就好办多了。

这一行：
“ ibrav=  4, celldm(1)=0.23372 ,celldm(3) =5.03366, nat=  3, ntyp= 2,”

二硫化钼晶格长度大概是3.16 angstrom，换算成bohr大概是5.97，你的celldm(1)错 ...

帮我看看哪有问题？？？

phonon at Gamma point

####
phonons of AlAs at Gamma
&inputph
  tr2_ph=1.0d-12,
  prefix='alas',
  trans=.true.,
  epsil=.true.,
  lraman=.true.,
  elop=.true.,
  amass(1)=32.06,
  amass(2)=95.96,
  outdir='./',
  fildyn='alas.dynG',
  fildrho='alas.drho',
/
0.0 0.0 0.0
###########

多谢！！！

log file:
_________________________________________________________
Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW
     Min          22      22      7                 1782     1782     309
     Max          23      23      8                 1799     1799     318
     Sum         367     367    121                28719    28719    5021


     Check: negative/imaginary core charge=   -0.000014    0.000000

     negative rho (up, down):  0.119E-02 0.000E+00

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     Error in routine phq_readin (1):

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
赞一下 回复此楼
6楼2013-08-16 21:13:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
引用回帖:
6楼: Originally posted by provels at 2013-08-16 21:13:11
帮我看看哪有问题？？？

phonon at Gamma point

####
phonons of AlAs at Gamma
&inputph
  tr2_ph=1.0d-12,
  prefix='alas',
  trans=.true.,
  epsil=.true.,
  lraman=.true.,
  elop=.tru ...

不好意思，我没算过声子谱，所以不知道哪里出了问题。
赞一下 回复此楼
7楼2013-08-16 23:19:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)
看起来你的输入文件是用QE中附带的例子改的。

输入文件最好自己手写，因为改已有的输入文件难保不会出现忘记修改某个参数的情况。可能算声子谱时出错是由文件名出错导致的。
赞一下 回复此楼
8楼2013-08-17 09:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 provels 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年也是没消息就是没中么 +3 vito刘 2026-06-13 3/150 2026-06-14 06:40 by zzahkj
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[找工作] 售SCI文章，我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] 囍 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[公派出国] T^T +5 q1zg6x4lf3 2026-06-12 5/250 2026-06-13 18:53 by lkmi35jdzt
[硕博家园] O(∩_∩)O +5 jk31rb6vyz 2026-06-12 7/350 2026-06-13 18:45 by lkmi35jdzt
[公派出国] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:33 by lkmi35jdzt
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o&gt;_&lt;o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考博] T^T +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 14:45 by 80j16ujcu8
[考研] O(∩_∩)O +3 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 14:42 by 80j16ujcu8
[考博] T^T +3 m584r2z4j3 2026-06-11 3/150 2026-06-13 12:04 by 80j16ujcu8
[找工作] 售SCI文章，我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[找工作] *囍 +3 m584r2z4j3 2026-06-12 3/150 2026-06-13 08:29 by t4hz8376wy
[论文投稿] →_→ +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 08:13 by t4hz8376wy
[基金申请] 面上项目没有好文章就没希望了吗？ +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
关闭
请填处理意见
关闭