计算化学公社

标题: 同一个晶体用cp2k先优化再做频率分析，得到的电子能量为什么不同？ [打印本页]

作者
Author: wanlichuan 时间: 2023-11-22 18:41
标题: 同一个晶体用cp2k先优化再做频率分析，得到的电子能量为什么不同？

同一个晶体用cp2k先优化再做频率分析，用的方法是一样的，用CELL_OPT得到的电子能量和用VIBRATIONAL_ANALYSIS得到的电子能量应该相同，但是我得到的并不相同，请大家帮忙看看问题出在哪儿：

CELL_OPT结果：

ENERGY| Total FORCE_EVAL ( QS ) energy[a.u.]: -4498.470708911696420

VIBRATIONAL_ANALYSIS结果：

Minimum Structure - Energy and Forces:

VIB| Totalenergy: -4498.126056127895

主要参数都是一样的，比如：

计算水平：PBEsol/ 6-31G*

CUTOFF 800

REL_CUTOFF 55

EPS_SCF 1.0E-10

SCHEMEMONKHORST-PACK 4 2 1

输入文件和输出文件已经放在附件中（主要原子坐标删去）。

请指教。

作者
Author: 吞木木 时间: 2023-11-23 21:12
个人觉得这俩不是一种能量，应该不是你说的只有电子能量。
官网给出的能量的解释

作者
Author: sobereva 时间: 2023-11-24 07:19
振动分析没带D3，优化带了

作者
Author: wanlichuan 时间: 2023-11-24 13:12

吞木木发表于 2023-11-23 21:12
个人觉得这俩不是一种能量，应该不是你说的只有电子能量。
官网给出的能量的解释

对，是有一点差别，但差别不会这么大。多谢回复。

作者
Author: wanlichuan 时间: 2023-11-24 13:13

sobereva 发表于 2023-11-24 07:19
振动分析没带D3，优化带了

哦，是的，当时还是没严格保持相同算法，忽略了D3。谢谢社长。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)