| 查看: 1555 | 回复: 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 如上述，真有这么三个能量值，我该看哪一个能量值来确定它是不是收敛？求大神赐教 |
回复此楼» 收录本帖的淘帖专辑推荐
 计算 |
» 猜你喜欢
271求调剂
已经有30人回复
-
材料类284调剂
已经有50人回复
-
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp 几何优化在先？参数优化在先？
已经有17人回复
- castep几何优化一直不收敛
已经有17人回复
- castep计算单点能的几何优化问题
已经有12人回复
- 晶体模型几何优化后，所得结果反而不如未优化的模型
已经有19人回复
- 几何优化
已经有7人回复
- 表面分子吸附的几何优化
已经有5人回复
- CASTEP几何优化问题
已经有9人回复
- MS几何优化好的超胞，算能量的时候为什么得不到结果
已经有6人回复
- 关于CASTEP的几何优化的两个疑惑。希望大家一起讨论讨论
已经有6人回复
- Gaussian几何优化后的键长不合理
已经有41人回复
- 利用CASTEP对建立的晶胞进行几何优化的问题 （急！！！）
已经有13人回复
- 【求助】咨询：高斯优化中局域能量最小值和全局最小值的问题
已经有6人回复
- 【求助】几何优化问题
已经有25人回复
- 【求助】Casteps几何优化中，遇到这样的提示通常该怎没办？
已经有3人回复
- 【求助】castep做几何优化时其参数如何设置？
已经有19人回复
- 【求助】还是收敛测试和几何优化的问题，又来麻烦大家了...
已经有3人回复
- 【求助】Gulp中的几何优化，MD等和Dmol3的中的有什么区别？
已经有11人回复
- 【求助】紧急求助几何优化后变形
已经有7人回复
2楼2012-12-03 21:37:28
fanqy1991
金虫 (著名写手)
- 应助: 32 (小学生)
- 金币: 7838.1
- 散金: 399
- 红花: 6
- 帖子: 1008
- 在线: 511.5小时
- 虫号: 1752732
- 注册: 2012-04-12
- 专业: 凝聚态物性 II ：电子结构
|
是的，都收敛， +-----------+-----------------+-----------------+------------+-----+ <-- 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
fanqy1991
金虫 (著名写手)
- 应助: 32 (小学生)
- 金币: 7838.1
- 散金: 399
- 红花: 6
- 帖子: 1008
- 在线: 511.5小时
- 虫号: 1752732
- 注册: 2012-04-12
- 专业: 凝聚态物性 II ：电子结构
4楼2012-12-04 12:46:24
5楼2012-12-05 21:56:31
fanqy1991
金虫 (著名写手)
- 应助: 32 (小学生)
- 金币: 7838.1
- 散金: 399
- 红花: 6
- 帖子: 1008
- 在线: 511.5小时
- 虫号: 1752732
- 注册: 2012-04-12
- 专业: 凝聚态物性 II ：电子结构
6楼2012-12-08 00:10:34
sesea7
银虫 (小有名气)
- 应助: 13 (小学生)
- 金币: 705.4
- 散金: 4
- 帖子: 56
- 在线: 26.3小时
- 虫号: 1965351
- 注册: 2012-08-31
- 性别: GG
- 专业: 理论环境化学
【答案】应助回帖
★ ★
liliangfang: 金币+2, 谢谢指教 2012-12-08 18:07:25
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
