24小时热门版块排行榜    

查看: 478  |  回复: 2
当前主题已经存档。

michelle8812

新虫 (小有名气)

[交流] [求助] MS中GULP模块计算表面能出错

使用MS的GULP模块计算Y3AlO12晶体(001)面的表面能,建好表面模型后优化计算,系统提示:

Start of surface optimisation :

  Cycle:      0 Energy:     -2439.461143  Gnorm:      3.619983  CPU:    2.672
  ** Hessian calculated **
  Cycle:      1 Energy:     -2789.807186  Gnorm:      2.999379  CPU:   10.391
  ** Hessian calculated **
  Cycle:      2 Energy:     -3268.548968  Gnorm:      2.141373  CPU:   19.500
  ** Hessian calculated **
  Cycle:      3 Energy:     -3496.883530  Gnorm:      2.090267  CPU:   28.875
  ** Hessian calculated **
  Cycle:      4 Energy:     -3962.867941  Gnorm:      1.027086  CPU:   38.156
  ** Hessian calculated **
  Cycle:      5 Energy:     -4309.378260  Gnorm:      0.708665  CPU:   47.344
  ** Hessian calculated **
  Cycle:      6 Energy:     -4310.303948  Gnorm:      0.710238  CPU:   56.625
  ** Hessian calculated **
  Cycle:      7 Energy:     -4310.359825  Gnorm:      0.710336  CPU:   66.016
  ** Hessian calculated **
  Cycle:      8 Energy:     -4311.057138  Gnorm:      0.711578  CPU:   75.125
  ** Hessian calculated **
  Cycle:      9 Energy:     -4311.057138  Gnorm:      0.711578  CPU:   84.266
  ** Hessian calculated **


  **** Conditions for a minimum have not been satisfied. However ****
  **** no lower point can be found - treat results with caution  ****
  **** unless gradient norm is small (less than 0.1)             ****


请问可能是什么原因导致出错呢?是表面模型建的不对吗?最后这三行字是什么意思?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★
michelle8812(金币+3,VIP+0):谢谢,我不太清楚是不是创建表面模型得的时候要满足那些要求,只是去了一个单位晶胞,我再好好研究下了
这是说,计算没有能够很好地收敛。你要小心处理以上数据,要求gnorm(即gradient norm,梯度模) 小于0.1才好。

出现这个情况,一般是因为模型(如势函数)不大好,或者初始构型不很好。

你可以这样:
1。将目前得到的优化构型做为初始构型接着算,或者稍微把这个构型变化一点点。
或者:2. 改变优化算法。
或者:3.改变你的势函数。

你可以单独下载一个GULP程序进行计算,灵活性比在MS中计算要大得多。
2楼2009-01-15 10:20:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

michelle8812

新虫 (小有名气)

谢谢
请问如何修改势函数或是如何改变优化算法
3楼2009-01-15 11:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 michelle8812 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 06:53 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 17/850 2026-11-01 06:48 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见