应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » pwscf计算integrated charge出错及不收敛
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 2513  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

体系见下图：
回复此楼

» 猜你喜欢
1楼 2012-07-21 14:28:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

cj4566

木虫 (正式写手)
继续求助啊，高人指点帮帮忙啊
赞一下 回复此楼
5楼2012-07-22 14:36:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)
引用回帖:
4楼: Originally posted by huazhorg at 2012-07-22 14:33:25
在命令行下敲入：grep Total outfile(你的输出文件名）。既然是优化，一般情况下力会越变越小，直到收敛到指定的阈值之内。看力才知道有没有收敛...

还是不会啊，计算是在远程机群上运行的，我在自己windows电脑上只看到out文件啊，而且我没有计算力啊
赞一下 回复此楼
6楼2012-07-23 07:03:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)
引用回帖:
4楼: Originally posted by huazhorg at 2012-07-22 14:33:25
在命令行下敲入：grep Total outfile(你的输出文件名）。既然是优化，一般情况下力会越变越小，直到收敛到指定的阈值之内。看力才知道有没有收敛...

我加了occupation现在没有那错误了，正在算，也看到force了，谢谢啊
赞一下 回复此楼
7楼2012-07-23 07:12:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)
引用回帖:
2楼: Originally posted by huazhorg at 2012-07-21 16:48:54
你grep Total outfile看看力的变化趋势

再请教一下啊，我看输出文件，一直在进行Self-consistent Calculation，每次自洽要迭代iteration十几步，然后给出k点，能量，力和坐标信息，请问这样的自洽计算到什么时候会结束啊，是依据哪个指标啊，谢谢！
赞一下 回复此楼
8楼2012-07-23 07:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)
可是试着把ecutwfc和ecutrho提高一些，比如ecutwfc=50,ecutrho=500.然后mix_beta=0.3,试试。
赞一下 回复此楼
frank_zhan,zy2zhan@gmail.com
9楼2012-07-23 08:54:43
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 cj4566 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +15 121. 2026-04-04 15/750 2026-04-07 02:14 by BruceLiu320
[考研] 304求调剂 +7 素年祭语 2026-04-06 10/500 2026-04-06 23:31 by 素年祭语
[考研] 338求调剂 +4 wxygxsaaaaa 2026-04-06 4/200 2026-04-06 17:14 by 尚水阁主
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 化学0703-一志愿211-338分求调剂 +8 vants 2026-04-05 8/400 2026-04-06 06:17 by houyaoxu
[考研] 生物工程求调剂 +6 喜欢还是不甘心 2026-04-05 6/300 2026-04-05 10:28 by 唐沐儿
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 085601一志愿北理325分求调剂 +6 找调剂，， 2026-04-02 6/300 2026-04-03 22:20 by 柟�?
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +5 redcom 2026-04-01 6/300 2026-04-03 15:15 by Howard28
[考研] 338求调剂，一志愿能源动力，外语是日语203 +5 zzz，，r 2026-04-02 5/250 2026-04-03 09:45 by 蓝云思雨
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 一志愿武汉理工0856，初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 08工科，295，接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
信息提示
关闭
请填处理意见
关闭