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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » 有关castep能量计算，请高手赐教！
131/2返回列表12下一页
查看: 2244  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

MVP41

木虫 (小有名气)

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

我用castep算了几个优化后的分子的能量，可是有一个最后算完没有Convergence.xcd和Energies.xcd这两个格式文件，但是完成时却显示的complete successfully，这样的结果，最后的算出来的能量值能用吗，还有为什么会出现这种情况，请各位高手不灵赐教，拜谢！
回复此楼

» 收录本帖的淘帖专辑推荐

晶体模型计算分析 xuexijisuan

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-05-07 10:55:40
已阅   回复此楼   关注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就说明收敛了
赞一下(1人) 回复此楼
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，还有就是请教一下怎样修改啊？本以为体系快要算完了，哎，无比失落啊！
赞一下(1人) 回复此楼
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计算能量，怎样判断算完的能量时正确可用的？再次拜谢！
赞一下(1人) 回复此楼
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 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭