24小时热门版块排行榜    

查看: 2727  |  回复: 12

MVP41

木虫 (小有名气)

[求助] 有关castep能量计算,请高手赐教!

我用castep算了几个优化后的分子的能量,可是有一个最后算完没有Convergence.xcd和Energies.xcd这两个格式文件,但是完成时却显示的complete successfully,这样的结果,最后的算出来的能量值能用吗,还有为什么会出现这种情况,请各位高手不灵赐教,拜谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lifeycc

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
御剑江湖: 金币+2, 谢谢~~ 2012-05-08 20:00:30
MVP41: 金币+4, ★★★很有帮助, 非常感谢 2012-05-09 08:54:10
BFGS: finished iteration    49 with enthalpy= -4.99445386E+003 eV
  
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    tolerance    |    units   | OK? | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
|  dE/ion   |   3.206652E-006 |   2.000000E-005 |         eV | Yes | <-- BFGS
|  |F|max   |   4.732415E-002 |   5.000000E-002 |       eV/A | Yes | <-- BFGS
|  |dR|max  |   1.426954E-003 |   2.000000E-003 |          A | Yes | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
  
BFGS: Geometry optimization completed successfully.
从底部往上看第一个这个表,三个yes就说明收敛了
superman爱喝摩卡
5楼2012-05-08 18:35:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MVP41

木虫 (小有名气)

引用回帖:
5楼: Originally posted by lifeycc at 2012-05-08 18:35:58:
BFGS: finished iteration    49 with enthalpy= -4.99445386E+003 eV
  
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    toleran ...

谢谢您的指点,突然发现我的结构基本上就没有收敛的,悲剧了,|  dE/ion   |   9.387614E-005 |   1.000000E-005 |         eV | No  | <-- BFGS
|  |F|max   |   2.180862E-001 |   3.000000E-002 |       eV/A | No  | <-- BFGS
|  |dR|max  |   2.390368E-003 |   1.000000E-003 |          A | No  | <-- BFGS
|   Smax    |   1.191771E-002 |   5.000000E-002 |        GPa | Yes | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
  
BFGS: WARNING - Geometry optimization failed to converge after        100 steps
很多都是这样的,还有全是NO的,只有个别的是全是yes,是不是这样的结构就根本不能用啊,但是每次算完时都是complete successfully,还有就是请教一下怎样修改啊?本以为体系快要算完了,哎,无比失落啊!
6楼2012-05-09 08:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MVP41

木虫 (小有名气)

引用回帖:
5楼: Originally posted by lifeycc at 2012-05-08 18:35:58:
BFGS: finished iteration    49 with enthalpy= -4.99445386E+003 eV
  
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    toleran ...

不好意思再麻烦您一下,小生还有一问,我用castep计算能量,怎样判断算完的能量时正确可用的?再次拜谢!
7楼2012-05-09 09:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

chaizhm

木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+1, 鼓励交流 2012-05-07 13:46:54
MVP41: 金币+1, 有帮助 2012-05-09 08:53:40
有用 我也遇到这样的情况,一般原子比较少就这样,具体原因我也不知道
壁立万仞,无欲则刚
2楼2012-05-07 13:23:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lifeycc

铁虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
chaizhm: 金币+1, 谢谢~ 2012-05-08 22:44:19
一般情况下只要*.castep文件中最后是3个YES,表示收敛了,能量就可以用的,应该是通讯问题那两个图没有download下来,没关系的。
superman爱喝摩卡
3楼2012-05-08 16:31:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MVP41

木虫 (小有名气)

引用回帖:
3楼: Originally posted by lifeycc at 2012-05-08 16:31:17:
一般情况下只要*.castep文件中最后是3个YES,表示收敛了,能量就可以用的,应该是通讯问题那两个图没有download下来,没关系的。

没见到过有三个yes呀,最后是这些啊
Initialisation time =      3.69 s
Calculation time    = 162618.61 s
Finalisation time   =     29.31 s
Total time          = 162651.62 s
Peak Memory Use     = 619688 kB
不过好像有这些就是收敛的!我还是想知道一下您所说的三个yes,请赐教!
4楼2012-05-08 16:56:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lifeycc

铁虫 (小有名气)

引用回帖:
6楼: Originally posted by MVP41 at 2012-05-09 08:58:32:
谢谢您的指点,突然发现我的结构基本上就没有收敛的,悲剧了,|  dE/ion   |   9.387614E-005 |   1.000000E-005 |         eV | No  | <-- BFGS
|  |F|max   |   2.180862E-001 |   3.000000E-002 |       ...

在castep提交任务里的More选项中有步数选择,一般默认100步,你可以把步数调大就行,建议都调500步。你可以在算完的的那个构型的基础上再提交一次任务,继续优化,达到收敛就行。
superman爱喝摩卡
8楼2012-05-09 09:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lifeycc

铁虫 (小有名气)

引用回帖:
7楼: Originally posted by MVP41 at 2012-05-09 09:32:27:
不好意思再麻烦您一下,小生还有一问,我用castep计算能量,怎样判断算完的能量时正确可用的?再次拜谢!

只要收敛了,然后把最后优化得到的构型算个能量,在castep计算中选择energy那一项,得到下面的数据:
Final energy, E             =  -6589.629949027     eV
Final free energy (E-TS)    =  -6589.711410673     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -6589.670679850     eV
然后看你需要,统一选用其中一个能量,我一般用NB est. 0K energy (E-0.5TS)
superman爱喝摩卡
9楼2012-05-09 09:46:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

MVP41

木虫 (小有名气)

引用回帖:
9楼: Originally posted by lifeycc at 2012-05-09 09:46:48:
只要收敛了,然后把最后优化得到的构型算个能量,在castep计算中选择energy那一项,得到下面的数据:
Final energy, E             =  -6589.629949027     eV
Final free energy (E-TS)    =  -6589.71141067 ...

非常感谢,看来虽然我的结构不收敛,但是能量计算还是没问题,呵呵,不过既然结构都不收敛,那能量也不能用了,做了这么长的工作还要重来啊,伤不起啊!
10楼2012-05-09 10:11:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 MVP41 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗?20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见