24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1572  |  回复: 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的回帖

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的回帖
查看全部 7 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 22408 312求调剂 +24 门路摸摸 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
信息提示
请填处理意见