计算化学公社

 找回密码 Forget password
 注册 Register
Views: 209|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] cp2k形成能(标准摩尔生成焓)计算结果和实验值有明显差异如何解决

[复制链接 Copy URL]

5

帖子

0

威望

73

eV
积分
78

Level 2 能力者

本帖最后由 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矫正因子还是因为基组泛函的选取呢,或者是石墨这种周期性结构也要设成三重态呢,如何才能减小这种差距,希望老师捞捞菜鸡

202511151359318877..png (23.47 KB, 下载次数 Times of downloads: 0)

202511151359318877..png

202511151358262801..png (23.53 KB, 下载次数 Times of downloads: 1)

202511151358262801..png

202511151356412875..png (23.85 KB, 下载次数 Times of downloads: 0)

202511151356412875..png

C.pdb

1.1 KB, 下载次数 Times of downloads: 0

Cfeq.inp

3.58 KB, 下载次数 Times of downloads: 3

Cfeq.out

108.26 KB, 下载次数 Times of downloads: 0

CO2.pdb

4.02 KB, 下载次数 Times of downloads: 0

CO2feq.inp

4.25 KB, 下载次数 Times of downloads: 0

CO2feq.out

73.01 KB, 下载次数 Times of downloads: 0

O2.pdb

3.09 KB, 下载次数 Times of downloads: 0

O2feq.inp

4.09 KB, 下载次数 Times of downloads: 2

O2feq.out

59.88 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2025-11-16 05:47:33 | 只看该作者 Only view this author
直接拿那么小的C的原胞做振动分析结合Shermo算热力学数据根本不靠谱,缺乏基本计算常识


北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)里面“热力学量的计算”部分有详细的介绍和很多完整细致的例子,正经学一遍才能避免胡算瞎算
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

5

帖子

0

威望

73

eV
积分
78

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2025-11-16 11:26:56 | 只看该作者 Only view this author
sobereva 发表于 2025-11-16 05:47
直接拿那么小的C的原胞做振动分析结合Shermo算热力学数据根本不靠谱,缺乏基本计算常识

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

4279

帖子

4

威望

9473

eV
积分
13832

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2025-11-16 13:28:10 | 只看该作者 Only view this author
你这个问题看起来更适合使用量子化学软件进行孤立体系的计算
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 00:01 , Processed in 0.165673 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list