计算化学公社

标题: cp2k形成能(标准摩尔生成焓）计算结果和实验值有明显差异如何解决 [打印本页]

作者
Author: WangSiTian 时间: 2025-11-15 14:11
标题: cp2k形成能(标准摩尔生成焓）计算结果和实验值有明显差异如何解决
本帖最后由 WangSiTian 于 2025-11-15 18:20 编辑

各位老师大家好，我打算高通量计算分子形成能，在CP2K使用二氧化碳作为一个示例，我使用的O2和CO2均为在MS中建模，C是在COD上找了个石墨结构。
均采用PBE泛函和DZVP-MOLOPT-SR-GTH基组进行结构优化和频率分析，所上传的为优化后的结构,而后在shermo上进行热力学量输出，ZPE矫正因子设定为0.9806。
因为石墨中的晶胞中含有两个碳原子，最后将它的焓除以二，拿二氧化碳减去氧气和单个碳的焓即-37.7480766+31.9208677+11.381193/2=-0.1366124a.u.=-358.67588352kJ/mol,而二氧化碳的标准摩尔生成焓约为-393kJ/mol,这中间的差异是由于什么形成的呢，是因为ZPE矫正因子还是因为基组泛函的选取呢，或者是石墨这种周期性结构也要设成三重态呢，如何才能减小这种差距，希望老师捞捞菜鸡

作者
Author: sobereva 时间: 2025-11-16 05:47
直接拿那么小的C的原胞做振动分析结合Shermo算热力学数据根本不靠谱，缺乏基本计算常识
(, 下载次数 Times of downloads: 2)

北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里面“热力学量的计算”部分有详细的介绍和很多完整细致的例子，正经学一遍才能避免胡算瞎算

作者
Author: WangSiTian 时间: 2025-11-16 11:26

sobereva 发表于 2025-11-16 05:47
直接拿那么小的C的原胞做振动分析结合Shermo算热力学数据根本不靠谱，缺乏基本计算常识

好的，谢谢卢老师，我也打算报名12月份的这个课，目前刚报名了初级量子化学培训班

作者
Author: zjxitcc 时间: 2025-11-16 13:28
你这个问题看起来更适合使用量子化学软件进行孤立体系的计算

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