24小时热门版块排行榜    

查看: 2488  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cj4566

木虫 (正式写手)

[求助] pwscf计算integrated charge出错及不收敛

我在对氧化锌吸附染料的体系进行relax优化的时候,出现下面提示:
WARNING: integrated charge=  1826.40000000, expected=  1827.00000000

每一次iteration后面都有这warning,但是正常计算直到完成100次iteration,最后显示没收敛,大家帮我看看啊
End of self-consistent calculation
convergence NOT achieved after 100 iterations: stopping

我的输入文件:
&CONTROL
                 calculation  = 'relax',
                 restart_mode = 'from_scratch',
                 prefix       = "ZnO-533-dye.rx",
                 pseudo_dir   = "/pkg/suse11/quantum/5.0/pseudo",
                 outdir       = "/home/n7388683/QE/tmp",                 
/
&SYSTEM
                       ibrav = 8,
                   celldm(1) = 30.77419871526555,
                   celldm(2) = 0.9614369051274179,
                   celldm(3) = 2.76327909118821,
                         nat = 244,
                        ntyp = 6,
                     ecutwfc = 25,
                     ecutrho = 250,
/
&ELECTRONS
                 conv_thr    = 1.D-6,
                 mixing_beta = 0.6,
/
&IONS
                ion_dynamics = 'bfgs',           
/
ATOMIC_SPECIES
   Zn   65.3799972534  Zn.pbe-van.UPF
    O   15.9989995956  O.pbe-van_ak.UPF
    C   12.0109996796  C.pbe-van_bm.UPF
    H    1.0080000162  H.pbe-van_ak.UPF
    N   14.0069999695  N.pbe-van_ak.UPF
    S   32.0600013733  S.pbe-van_bm.UPF
ATOMIC_POSITIONS crystal
。。。
K_POINTS automatic
  1 1 1   0 0 0

体系见下图:
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
9楼: Originally posted by frank_zhan at 2012-07-23 08:54:43
可是试着把ecutwfc和ecutrho提高一些,比如ecutwfc=50,ecutrho=500.然后mix_beta=0.3,试试。

好的  谢谢啊,我在原来的输入文件里添加了occupations = 'smearing', degauss = 0.02D0, 没有提示电子数错误,现在正在算,完了再提高ecut试试
10楼2012-07-23 09:22:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

huazhorg

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-07-21 16:53:28
cj4566: 金币+5 2012-07-23 07:03:56
cj4566: 金币+5 2012-07-23 07:13:02
你grep Total outfile看看力的变化趋势
2楼2012-07-21 16:48:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
2楼: Originally posted by huazhorg at 2012-07-21 16:48:54
你grep Total outfile看看力的变化趋势

grep Total outfile是怎么用的啊   电子数不对为什么要看力的变化呢??
3楼2012-07-22 14:10:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)


sunyang1988: 金币+1, 谢谢交流 2012-07-22 16:16:00
引用回帖:
3楼: Originally posted by cj4566 at 2012-07-22 14:10:50
grep Total outfile是怎么用的啊   电子数不对为什么要看力的变化呢??...

在命令行下敲入:grep Total outfile(你的输出文件名)。既然是优化,一般情况下力会越变越小,直到收敛到指定的阈值之内。看力才知道有没有收敛
4楼2012-07-22 14:33:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +14 于海海海海 2026-03-19 14/700 2026-03-19 19:40 by maocaozhuxi
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 求调剂 +3 Mqqqqqq 2026-03-19 3/150 2026-03-19 14:11 by peike
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
信息提示
请填处理意见