24小时热门版块排行榜    

查看: 1803  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:24 by jHcrHMhbKCFV
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:23 by jHcrHMhbKCFV
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 19/950 2026-08-02 16:21 by jHcrHMhbKCFV
[基金申请] 面上提前没消息,有中的吗 +9 archvillain 2026-08-02 11/550 2026-08-02 15:44 by alexon
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 8/400 2026-08-02 15:01 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 14:58 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 13:34 by jHcrHMhbKCFV
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
信息提示
请填处理意见