同一个晶体用cp2k先优化再做频率分析，得到的电子能量为什么不同？

wanlichuan

同一个晶体用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

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

请指教。

吞木木

个人觉得这俩不是一种能量，应该不是你说的只有电子能量。
官网给出的能量的解释

sobereva

振动分析没带D3，优化带了

wanlichuan

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

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

wanlichuan

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

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