24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1573  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 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 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 药学求调剂 +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
信息提示
请填处理意见