24小时热门版块排行榜

返回列表

cj4566

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 17
散金: 24
红花: 3
帖子: 473
在线: 399.8小时
虫号: 356133
注册: 2007-04-26
性别: GG
专业: 无机非金属类光电信息与功

[求助] 非自洽计算中断报错

我用pwscf对ZnO表面模型进行非自洽计算的时候总是中断，大家帮忙看看啊，谢谢了！
错误信息：
   ethr =  3.47E-10,  avg # of iterations = 61.6
   Message from routine efermit:
   internal error, cannot bracket Ef
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%    from tweights : error #       1
   bad Fermi energy
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   stopping ...

输入文件：
&CONTROL
            calculation  = 'nscf',
            pseudo_dir = "/pkg/suse11/quantum/5.0/pseudo",
            outdir    = "/home/n7388683/QE/tmp",
            prefix    = "100scf",
/
&SYSTEM
                     ibrav = 8,
               celldm(1) = 12.42306308591684,
               celldm(2) = 1.61271676300578,
               celldm(3) = 2.859066169759659,
                     nat = 32,
                     ntyp = 2,
                     nbnd = 30,
               occupations = 'tetrahedra',
                  ecutwfc = 30,
                  ecutrho = 300,
/
&ELECTRONS
         diagonalization = 'david'
/
ATOMIC_SPECIES
Zn 65.3799972534  Zn.pbe-van.UPF
O 15.9989995956  O.pbe-van_ak.UPF
ATOMIC_POSITIONS crystal
O       0.263888864 0.344244423  -0.001685542
O       0.013888845 0.343515874 0.148543251
O       0.013888801 0.093520173 0.053390015
O       0.263888855 0.094247127 0.203632594
O       0.763888929 0.344244426  -0.001685670
O       0.513888956 0.343516022 0.148543308
O       0.513888890 0.093520319 0.053390042
O       0.763889048 0.094247130 0.203632690
O       0.263888731 0.844247062  -0.001684524
O       0.013888855 0.843518205 0.148542338
O       0.013888867 0.593519332 0.053387801
O       0.263888857 0.594246737 0.203631467
O       0.763889008 0.844246982  -0.001684655
O       0.513888910 0.843518345 0.148542355
O       0.513888853 0.593519414 0.053387765
O       0.763888942 0.594246802 0.203631536
Zn    0.263888921 0.173444211 0.024427317
Zn    0.013888927 0.159770140 0.166488906
Zn    0.013888878 0.409759797 0.035464366
Zn    0.263888824 0.423464137 0.177485655
Zn    0.763888787 0.173444233 0.024427220
Zn    0.513889002 0.159770379 0.166489335
Zn    0.513888941 0.409759858 0.035464244
Zn    0.763889020 0.423464160 0.177485798
Zn    0.263888817 0.673445736 0.024425680
Zn    0.013888949 0.659772028 0.166489391
Zn    0.013888826 0.909759911 0.035467283
Zn    0.263888943 0.923464754 0.177487167
Zn    0.763888925 0.673445781 0.024425793
Zn    0.513888892 0.659772030 0.166489391
Zn    0.513888940 0.909760022 0.035467025
Zn    0.763888851 0.923464777 0.177487253
K_POINTS automatic
  10 10 1 0 0 0

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

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

1楼 2012-08-10 16:34:17

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

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
御剑江湖: 金币+2, 谢谢 2012-08-10 21:35:47
cj4566: 金币+10, ★★★很有帮助 2012-08-13 06:51:04

你的计算scf的时候没有出现类似的问题？二者差别应该不大的，除非你的scf计算很不准确。
原则上，要解决这个问题，需要修改程序的。我以前就干过，不过我的体系特别，Fermi能级很高。你这个体系Feimi能级难不成也很大？

赞一下

回复此楼

2楼2012-08-10 18:39:32

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

goldenfisher

金虫 (著名写手)

1ST强帖: 3
应助: 55 (初中生)
贵宾: 0.153
金币: 5779.1
散金: 2908
红花: 21
帖子: 1552
在线: 351.5小时
虫号: 117764
注册: 2005-11-26
性别: GG
专业: 原子和分子物理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
御剑江湖: 金币+1, 谢谢 2012-08-10 21:35:57
cj4566: 金币+10, ★★★很有帮助 2012-08-13 06:51:11

我看到了，你的nbnd才30，你在搞笑。

赞一下

回复此楼

3楼2012-08-10 18:40:13

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

cj4566

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 17
散金: 24
红花: 3
帖子: 473
在线: 399.8小时
虫号: 356133
注册: 2007-04-26
性别: GG
专业: 无机非金属类光电信息与功

引用回帖:

2楼: Originally posted by goldenfisher at 2012-08-10 18:39:32
你的计算scf的时候没有出现类似的问题？二者差别应该不大的，除非你的scf计算很不准确。
原则上，要解决这个问题，需要修改程序的。我以前就干过，不过我的体系特别，Fermi能级很高。你这个体系Feimi能级难不成也很 ...

scf计算正常的，得到的费米能级是-1.3505 ev，这算是正常吧。

赞一下

回复此楼

4楼2012-08-13 06:48:39

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