24小时热门版块排行榜    

查看: 1349  |  回复: 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 的主题更新
信息提示
请填处理意见