24小时热门版块排行榜    

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

xiao_ma523

铁杆木虫 (著名写手)

[求助] MS中CASTEP计算结构优化时不收敛问题,请各位指教。。。

刚开始学习使用 ms castep ,小白一只,在进行结构优化计算的时候,出现如下问题,烦请各位前辈们指教。。。。谢谢。。。
(计算函数选取为GGA-PBE;计算参数选择分别为:k-points:20×20×20;SCF设置:SCF tolerance:5.0e-7eV/atom,Max. SCF sysles: 100,Electronic minimizer:All B ands/EDFT
 http://

------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -1.05658475E+003  4.29446436E+001                        40.51  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
      1  -1.66908169E+003  2.06199975E+001   3.06248473E+002     201.90  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
      2  -1.71727432E+003  1.54132235E+001   2.40963134E+001     368.41  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
      3  -1.73512023E+003  1.31966544E+001   8.92295366E+000     505.40  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
      4  -1.74004485E+003  1.24138206E+001   2.46230895E+000     636.64  <-- SCF
      5  -1.74135250E+003  1.21899463E+001   6.53828871E-001     773.64  <-- SCF
      6  -1.74205719E+003  1.20787843E+001   3.52340973E-001     907.32  <-- SCF
      7  -1.74259344E+003  1.19953030E+001   2.68127227E-001    1037.62  <-- SCF
      8  -1.74306614E+003  1.19286041E+001   2.36350991E-001    1177.77  <-- SCF
      9  -1.74358796E+003  1.18488409E+001   2.60909149E-001    1298.85  <-- SCF
     10  -1.74403744E+003  1.17938519E+001   2.24738643E-001    1421.06  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     11  -1.74455876E+003  1.17365970E+001   2.60662432E-001    1547.33  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     12  -1.74497637E+003  1.16890613E+001   2.08802683E-001    1672.85  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     13  -1.74535476E+003  1.16477635E+001   1.89198766E-001    1784.51  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     14  -1.74569639E+003  1.16017743E+001   1.70815013E-001    1896.44  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     15  -1.74603497E+003  1.15495874E+001   1.69286593E-001    2027.35  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     16  -1.74633929E+003  1.15088101E+001   1.52159230E-001    2140.21  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     17  -1.74664029E+003  1.14673428E+001   1.50503130E-001    2250.68  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     18  -1.74691967E+003  1.14291301E+001   1.39687296E-001    2383.77  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     19  -1.74716002E+003  1.13983485E+001   1.20175874E-001    2524.37  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     20  -1.74738725E+003  1.13643370E+001   1.13616354E-001    2637.47  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     21  -1.74764344E+003  1.13346137E+001   1.28091921E-001    2761.41  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     22  -1.74783115E+003  1.13112076E+001   9.38565737E-002    2880.58  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     23  -1.74801534E+003  1.12912355E+001   9.20978199E-002    3019.49  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
     24  -1.74815772E+003  1.12827309E+001   7.11860533E-002    3166.94  <-- SCF
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
     25  -1.74831435E+003  1.12683472E+001   7.83167816E-002    3295.82  <-- SCF
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
     26  -1.74846036E+003  1.12568548E+001   7.30063648E-002    3429.77  <-- SCF
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016
     27  -1.74859810E+003  1.12399413E+001   6.88709270E-002    3565.08  <-- SCF
     28  -1.74872556E+003  1.12249860E+001   6.37300369E-002    3706.03  <-- SCF
     29  -1.74883372E+003  1.12156147E+001   5.40795198E-002    3822.78  <-- SCF
*Warning* predicted gradient is positive



主要问题就是两个warning怎么才能解决:
*Warning* predicted gradient is positive
*Warning* step too small in line minimisation; setting to
  2.220446049250313E-016

[ Last edited by xiao_ma523 on 2013-1-25 at 10:22 ]
回复此楼

» 猜你喜欢

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

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

磊俊陈

木虫 (小有名气)

【答案】应助回帖

★ ★
xiao_ma523(franch代发): 金币+2, 谢谢回帖交流,, 2013-04-22 20:38:15
楼主的K点设置的确实有点多了,先用默认参数试试,成功之后再往精度高了调
以一颗学习的心态面对一切未知的事物!
10楼2013-04-17 19:52:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

xw2008

铜虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2013-01-26 15:22:31
k-points:20×20×20设的挺多的 你算的是什么材料啊?最好有把生成的.castep文件,才好判断。
Max. SCF sysles先把设大点,试试。
2楼2013-01-25 13:34:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoupeng87

版主 (文学泰斗)

优秀版主优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
。。。。。你这个怎么整的这么复杂的
3楼2013-01-25 14:17:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao_ma523

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by xw2008 at 2013-01-25 13:34:52
k-points:20×20×20设的挺多的 你算的是什么材料啊?最好有把生成的.castep文件,才好判断。
Max. SCF sysles先把设大点,试试。

上边的问题就是在.castep里边截出来的。。。算的是Ti基多元固溶体,采用的是虚晶近似。我试一下,可是Max. SCF cycles要调高多少比较好呢。。。
4楼2013-01-25 14:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 20:40 by jyNB2vrbdUpR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 20:31 by jyNB2vrbdUpR
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 uFOAtb1CDngh 2026-09-02 9/450 2026-09-03 20:20 by jyNB2vrbdUpR
[考博] 找导师 +7 yuanjiabao 2026-08-29 8/400 2026-09-03 20:15 by yinnelly
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ero8OE6tv9cu 2026-09-02 5/250 2026-09-03 20:09 by jyNB2vrbdUpR
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ero8OE6tv9cu 2026-09-02 9/450 2026-09-03 20:07 by jyNB2vrbdUpR
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +3 jackleilei 2026-09-03 3/150 2026-09-03 19:13 by god_tian
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 18:50 by jyNB2vrbdUpR
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 rFEsUBKXRll0 2026-09-02 6/300 2026-09-03 16:49 by T0rGB46095mJ
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 16:38 by T0rGB46095mJ
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 10/500 2026-09-03 16:27 by T0rGB46095mJ
[基金申请] 国社科又开始会评了,不知道这次命运如何 +10 雨打竹帘 2026-08-30 14/700 2026-09-03 13:17 by qsd10086
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
请填处理意见