24小时热门版块排行榜    

查看: 5067  |  回复: 8

王梦ariel

新虫 (小有名气)

[求助] castep几何优化的收敛性问题。Max.force 一直不收敛,其他的参数已经收敛

最近在做有机物在晶面上的吸附计算,但是设了构型之后,吸附构型的几何优化,一直都不能收敛。下面是我的收敛的结果图和参数设置图,请问下怎么设置能够加快我的收敛速度呀!而且计算中出现了一些warning,应该怎么修改,才能找到吸附的稳定结构。

castep几何优化的收敛性问题。Max.force 一直不收敛,其他的参数已经收敛
收敛趋势图.GIF


castep几何优化的收敛性问题。Max.force 一直不收敛,其他的参数已经收敛-1
参数设置.GIF


castep几何优化的收敛性问题。Max.force 一直不收敛,其他的参数已经收敛-2
warning.GIF
回复此楼

» 猜你喜欢

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

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★ ★
gmy1990: 金币+3 2014-10-03 00:23:58
试试增加k或是增加一下smearing,k不充分的条件下不考虑smearing导致的结果就是力的计算不可靠,从而导致几何优化收敛困难。当然,截断能不充分也有可能导致问题,所以建议做测试。

如果其他参数都合理,那么只好不采用BFGS了,不过我对阻尼动力学实在不太熟悉。如果只是初始结构不好也不妨试试Dmol3做个预先优化,收敛后再用CASTEP优化。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2014-10-02 23:42:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王梦ariel

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 卡开发发 at 2014-10-02 23:42:05
试试增加k或是增加一下smearing,k不充分的条件下不考虑smearing导致的结果就是力的计算不可靠,从而导致几何优化收敛困难。当然,截断能不充分也有可能导致问题,所以建议做测试。

如果其他参数都合理,那么只好 ...

请问收敛性测试也是用几何优化的算法,参数从低到高做吗?我的体系比较大,计算太慢了,下面是我的一个几何优化的收敛图,从前面看,三个收敛标准都是往收敛的方向走的,但是中途又开始增加和震荡,这个参数都只是设到medium,已经走了20天左右了,这个应该怎么办呀!初始的构型也是在discover里做了一下。不知道测试应该怎么测?随便一算都是好几天不收敛。
castep几何优化的收敛性问题。Max.force 一直不收敛,其他的参数已经收敛-3
收敛图.GIF

3楼2014-10-03 10:06:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by 王梦ariel at 2014-10-03 10:06:15
请问收敛性测试也是用几何优化的算法,参数从低到高做吗?我的体系比较大,计算太慢了,下面是我的一个几何优化的收敛图,从前面看,三个收敛标准都是往收敛的方向走的,但是中途又开始增加和震荡,这个参数都只是 ...

不妨把结构截图给出来,以及*.castep文件,否则我没办法识别参数以及迭代的信息。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2014-10-03 10:41:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王梦ariel

新虫 (小有名气)

引用回帖:
4楼: Originally posted by 卡开发发 at 2014-10-03 10:41:28
不妨把结构截图给出来,以及*.castep文件,否则我没办法识别参数以及迭代的信息。...

这是我的结构截图和我的castep文件。其中结构是FeOOH的100晶面,上下表面采用H和OH饱和,吸附分子为CH3-CONOH这样一个含有氧肟酸基团的小分子,谢谢你帮助。希望你能看出不合理的地方。
castep几何优化的收敛性问题。Max.force 一直不收敛,其他的参数已经收敛-4
结构截图.GIF

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 3.castep
  • 2014-10-03 18:03:35, 11.68 M
5楼2014-10-03 18:04:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur


gmy1990: 金币+1 2014-10-05 18:51:58
引用回帖:
5楼: Originally posted by 王梦ariel at 2014-10-03 18:04:07
这是我的结构截图和我的castep文件。其中结构是FeOOH的100晶面,上下表面采用H和OH饱和,吸附分子为CH3-CONOH这样一个含有氧肟酸基团的小分子,谢谢你帮助。希望你能看出不合理的地方。

结构截图.GIF
...

初始结构最好别用力场去处理,有时候反倒是糟糕。另外你可以看看xtd文件,看看结构优化的走势,手工做调整。另外,单独的基体是否能够结构优化成功?
不一定挂在论坛,计算问题问题欢迎留言。
6楼2014-10-03 19:51:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王梦ariel

新虫 (小有名气)

引用回帖:
6楼: Originally posted by 卡开发发 at 2014-10-03 19:51:28
初始结构最好别用力场去处理,有时候反倒是糟糕。另外你可以看看xtd文件,看看结构优化的走势,手工做调整。另外,单独的基体是否能够结构优化成功?...

几何优化的话就会有力的收敛条件这一项,那要用别的模块做初始结构吗?
7楼2014-10-04 20:02:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王梦ariel

新虫 (小有名气)

引用回帖:
6楼: Originally posted by 卡开发发 at 2014-10-03 19:51:28
初始结构最好别用力场去处理,有时候反倒是糟糕。另外你可以看看xtd文件,看看结构优化的走势,手工做调整。另外,单独的基体是否能够结构优化成功?...

几何优化就会有力场的收敛条件,难道要用其他的模块做优化吗?
8楼2014-10-04 20:03:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

★ ★
gmy1990: 金币+2 2014-10-05 18:51:09
引用回帖:
7楼: Originally posted by 王梦ariel at 2014-10-04 20:02:07
几何优化的话就会有力的收敛条件这一项,那要用别的模块做初始结构吗?...

力收敛和使用分子力场优化完全不是一回事,没必要用其他模块做预先的优化,如果你觉得初始结构不好可以降低计算精度作为预优化。
不一定挂在论坛,计算问题问题欢迎留言。
9楼2014-10-04 20:06:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 王梦ariel 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见