24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1656  |  回复: 7

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


多谢!!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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,这个选项对应直角坐标而不是分数坐标。

算二硫化钼要加展宽?我记得不加展宽也能收敛。
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 的主题更新
信息提示
请填处理意见