24小时热门版块排行榜    

查看: 1330  |  回复: 2

Dirac12345

新虫 (小有名气)

[交流] QE nscf 计算错误 已有2人参与

各位好 小弟计算 nscf 准备 用作下一步的dos计算
小弟有两个问题:
一是 nscf计算时 nbnd需要手动设置吗? 还是用默认值即可
二是不明白为什么出现了以下错误,附小弟的输入文件。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #        44
     from iosys : error #         1
     occupations tetrahedralnot implemented
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@^@     from iosys : error #         1
     occupations tetrahedralnot implemented


Error in routine iosys (1):
     occupations tetrahedralnot implemented


===================
&control
    title = unitcell
    calculation = 'nscf',
    restart_mode='from_scratch',
    prefix='unitcell_vdw_scf',   
    pseudo_dir = '/auto/rcf-proj2
    outdir= './',
!    forc_conv_thr = 1.0D-3,  !1.
!    etot_conv_thr = 1.0D-4, !1.0
    nstep= 1,
/

&system
    ibrav= 0,
    celldm(1)=1.8897,
    nat= 48,
    ntyp= 5,
    ecutwfc = 40.D0,
    ecutrho   = 300.D0, !for norm
    occupations = 'tetrahedral',
!    smearing = 'gauss', !gauss i
!    degauss = 0.02,
    input_dft = 'vdw-DF2'  !or vd
   nbnd = 150       !for dos and
/

&electrons
!    mixing_mode = 'plain'      !
!    mixing_beta = 0.1D0        !
    conv_thr =  1.0d-8           
    electron_maxstep = 100
/

&IONS !use only calculation = rel
  ion_dynamics = 'bfgs'         
  pot_extrapolation = "second_ord
  wfc_extrapolation = "second_ord
  trust_radius_max = 0.2D0,
  trust_radius_ini = 0.1D0,
==================

谢谢
回复此楼

» 猜你喜欢

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

都是斯坦惹的祸
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
ljw4010: 金币+1, 谢谢交流! 2015-04-08 22:45:54
nscf不能写这个:
occupations = 'tetrahedra'
2楼2015-04-08 13:57:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
因为拼写错误啊,是tetrahedra,不是tetrahedral。
另外,ibrav=0时需要指定所有六个celldm(i),或者用CELL_PARAMETERS来定义晶格。
顺便,nbnd需要指定多少,跟自己想要看多大的导带有关。当然,nbnd大了自然计算时间就要长一些了~
思想重于技巧,内涵重于表象
3楼2015-04-19 21:11:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Dirac12345 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600 材料与化工 329分求调剂 +8 Mr. Z 2026-03-25 8/400 2026-03-25 23:00 by peike
[考研] 材料与化工328分调剂 +6 。,。,。,。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 289材料与化工(085600)B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 一志愿重庆大学085700资源与环境,总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 求调剂 +7 Auroracx 2026-03-22 7/350 2026-03-22 12:38 by 素颜倾城1988
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
信息提示
请填处理意见