24小时热门版块排行榜    

查看: 1173  |  回复: 2

杨寻

铁杆木虫 (著名写手)

[求助] pwscf计算费米面出错

1. 自恰计算 INPUT:
&control
calculation='scf'
restart_mode='from_scratch',
pseudo_dir = './',
outdir='./'
prefix='cu'
/

&system
ibrav = 2, celldm(1) =6.73, nat= 1, ntyp= 1,
ecutwfc = 25.0, ecutrho = 300.0
occupations='smearing', smearing='gaussian', degauss=0.02
/

&electrons
diagonalization='david'
conv_thr = 1.0e-8
mixing_beta = 0.7

/
ATOMIC_SPECIES
Cu 63.55 Cu.pz-d-rrkjus.UPF
ATOMIC_POSITIONS
Cu 0.0 0.0 0.0
K_POINTS (automatic)
8 8 8 0 0 0

这一步正常

2. 得到较密K点 并且 非自恰计算: INPUT

&control
calculation='nscf'
prefix='cu'
pseudo_dir='./'
outdir='./'
/

&system
ibrav=2
celldm(1)=6.73
nat=1
ntyp=1
ecutwfc=25.0
ecutrho=300
nbnd=8
/

&electrons
conv_thr=1.0e-8
mixing_beta=0.7
/

ATOMIC_SPECIES
Cu 63.55 Cu.pz-d-rrkjus.UPF

ATOMIC_POSITIONS
Cu 0.0 0.0 0.0

K_POINTS
   216
    0.000000    0.000000    0.000000  1.00
   -0.200000    0.200000   -0.200000  1.00
   -0.400000    0.400000   -0.400000  1.00
....
....

执行命令 pw.x < cu.fs.in > cu.fs.out
出现如下错误:


     Program PWSCF v.4.3.2      starts on  9Apr2014 at 22:29:17

     This program is part of the open-source Quantum ESPRESSO suite
     for quantum simulation of materials; please cite
         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
          URL http://www.quantum-espresso.org",
     in publications or presentations arising from this work. More details at
     http://www.quantum-espresso.org/ ... ng_Quantum-ESPRESSO

     Parallel version (MPI), running on     8 processors
     R & G space division:  proc/pool =    8

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
     Waiting for input...
     Reading input from stdin
               file Cu.pz-d-rrkjus.UPF: wavefunction(s)  3D renormalized

     Atomic positions and unit cell read from directory:
     ./cu.save/


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

  *            1  *           19  /=           59
  *            2  *           19  /=           59
  *            3  *           19  /=           60
  *            4  *           18  /=           60
  *            5  *           17  /=           60
  *            6  *           19  /=           59
  *            7  *           19  /=           59
  *            8  *           19  /=           59

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from  fft_dlay_set  : error #         7
      inconsistent number of sticks
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

google也找不到原因,请教大家帮忙 谢谢。。
计算费米面的方法参见:
http://valenhou.blog.edu.cn/home ... g&id=138335
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

确定google过了?
我google出来的,仅供参考:
http://www.democritos.it/piperma ... January/001876.html
http://qe-forge.org/pipermail/pw_forum/2007-September/081723.html
我相信你能看懂。
看了以后要是还不知道怎么做,我们再讨论讨论
2楼2014-04-29 08:58:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

火力悠92

新虫 (初入文坛)

lz,qe初学者,想问一下,产生计算费米面的k点的脚本kvecs_FS.x是qe自带的脚本么,表示没有找到

发自小木虫Android客户端
3楼2017-04-20 13:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 杨寻 的主题更新
信息提示
请填处理意见