24小时热门版块排行榜    

查看: 1980  |  回复: 6

fanqy1991

金虫 (著名写手)

[求助] 几何优化有三个能量值 我该看哪一个?

Final energy, E             =  -12266.31912109     eV
Final free energy (E-TS)    =  -12266.33460092     eV
(energies not corrected for finite basis set)

NB est. 0K energy (E-0.5TS)      =  -12266.32686101     eV
如上述,真有这么三个能量值,我该看哪一个能量值来确定它是不是收敛?求大神赐教
回复此楼

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

计算

» 猜你喜欢

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

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

anlen0615

银虫 (小有名气)

【答案】应助回帖

★
感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-12-06 20:06:41
确定是不是收敛  主要看.castep文件中BFGS是不是都是YES
   要求体系最后的能量  要在优化好的模型基础上算Energy  得出的能量那才是有效的!
厚积而薄发!
2楼2012-12-03 21:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanqy1991

金虫 (著名写手)

引用回帖:
2楼: Originally posted by anlen0615 at 2012-12-03 21:37:28
确定是不是收敛  主要看.castep文件中BFGS是不是都是YES
   要求体系最后的能量  要在优化好的模型基础上算Energy  得出的能量那才是有效的!

是的,都收敛,

+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
| Parameter |      value      |    tolerance    |    units   | OK? | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
|  dE/ion   |   7.491404E-008 |   5.000000E-006 |         eV | Yes | <-- BFGS
|  |F|max   |   8.616090E-003 |   1.000000E-002 |       eV/A | Yes | <-- BFGS
|  |dR|max  |   6.978532E-005 |   5.000000E-004 |          A | Yes | <-- BFGS
|   Smax    |   6.914192E-003 |   2.000000E-002 |        GPa | Yes | <-- BFGS
+-----------+-----------------+-----------------+------------+-----+ <-- BFGS
但是我不知道该看哪一个能量值,您能帮我明确指出来吗?
3楼2012-12-04 12:44:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanqy1991

金虫 (著名写手)

引用回帖:
2楼: Originally posted by anlen0615 at 2012-12-03 21:37:28
确定是不是收敛  主要看.castep文件中BFGS是不是都是YES
   要求体系最后的能量  要在优化好的模型基础上算Energy  得出的能量那才是有效的!

就是上面三个能量值我该看哪一个?不是说要确定参数的时候要能量间隔小于0.001eV/atom吗?
4楼2012-12-04 12:46:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anlen0615

银虫 (小有名气)

【答案】应助回帖

★
sunyang1988: 金币+1, 谢谢交流 2012-12-06 20:06:51
引用回帖:
4楼: Originally posted by fanqy1991 at 2012-12-04 12:46:24
就是上面三个能量值我该看哪一个?不是说要确定参数的时候要能量间隔小于0.001eV/atom吗?...

不是看优化里面的能量    而是在优化后的结构基础上再算一下Energy
那里面就只会有一个能量!   那才是有效的!
厚积而薄发!
5楼2012-12-05 21:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanqy1991

金虫 (著名写手)

引用回帖:
5楼: Originally posted by anlen0615 at 2012-12-05 21:56:31
不是看优化里面的能量    而是在优化后的结构基础上再算一下Energy
那里面就只会有一个能量!   那才是有效的!...

啊?真的?
6楼2012-12-08 00:10:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sesea7

银虫 (小有名气)

【答案】应助回帖

★ ★
liliangfang: 金币+2, 谢谢指教 2012-12-08 18:07:25
lz是用castep计算的吧,收敛了,是有三个能量值。具体参见help文件,我复制粘贴一部分:
The CASTEP output is as follows. The line

Final energy, E             =  -861.8315373731     eV

gives the Kohn-Sham energy, EKS, which ignores the entropy. However, it is the total Mermin free energy that is minimized by CASTEP and reported in the SCF cycles in the line

Final free energy (E-TS)    =  -861.8315374521     eV

If you were actually interested in the electronic states at some reasonable nonzero temperature, you would smear the Fermi level with a thermal (i.e., Fermi-Dirac) smearing method and the free energy reported by CASTEP would be the value you require.

However, it is often the T = 0 K value that is required and the smearing is necessary to improve the convergence of the numerical algorithm.  Because the dependence of the free energy and the Kohn-Sham energy on the smearing width is the same (to first order), averaging the two values gives a much better approximation to the T = 0 K result:

NB est. 0K energy (E-0.5TS) =  -861.8315374126     eV
喵星人
7楼2012-12-08 10:38:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanqy1991 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
信息提示
请填处理意见