24小时热门版块排行榜    

查看: 1283  |  回复: 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的回帖

火力悠92

新虫 (初入文坛)

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

发自小木虫Android客户端
3楼2017-04-20 13:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上复盘~欢迎讨论 +8 晴天加油 2026-08-26 9/450 2026-08-27 14:14 by gltch
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +6 苏东坡二世 2026-08-23 6/300 2026-08-27 13:57 by zhangjianl
[基金申请] 青基如何破局 +3 lyz123lyz 2026-08-27 4/200 2026-08-27 13:08 by kulium
[基金申请] 基金未中,这种答复是模板吗? +4 zhaosm1982 2026-08-27 5/250 2026-08-27 12:17 by zhaosm1982
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 面上合作单位盖章 +3 ssyjh 2026-08-27 3/150 2026-08-27 09:59 by 阿拉贡aragon
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
请填处理意见