计算化学公社

标题: 关于cp2k做两个小分子的吸附的能量差值问题 [打印本页]

作者
Author: liuxp 时间: 2025-4-28 16:05
标题: 关于cp2k做两个小分子的吸附的能量差值问题
目前催化方向做了两个小分子的吸附，结果如图所示，两个有机物的在结构上没有区别，仅仅是由单键变成了双键，两个小分子通过shermo计算出来自由能的差值在14 eV ,根据文献这两个有机物的差值在2eV 左右，目前计算能量的级别是PBE(D3)/TZV2P 算出来差值依然很大，目前有如下几个问题
1.基组已经用到 TZV2P-MOLOPT-GTH的级别，是否需要加到 TZV2PX-MOLOPT-GTH （目前这个方案已提交超算）
2.对于这两个有机物的优化仅限于在MS里边的dmol模块画完以后0维体系的优化，不确定是不是最优的结构
3.由于吸附物是有机物不确定在设置上是否需要其他的变动

作者
Author: sobereva 时间: 2025-4-28 16:31
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread-38911-1-1.html）

作者
Author: liuxp 时间: 2025-4-28 17:04

sobereva 发表于 2025-4-28 16:31
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678（http://bbs.keinsci.com/thread ...

sob 老师并没有复现文献，算别的都没问题，只是算有机物，两个有机物差了单纯的一根键的问题，现在算出来两个自由能差距太大，不知道问题出在哪

作者
Author: Dempey 时间: 2025-4-28 17:15
你的1.png比2.png多一个氢，组成不一样不能比较绝对值，你可以比较吸附能的差异

作者
Author: sobereva 时间: 2025-4-28 17:22

liuxp 发表于 2025-4-28 17:04
sob 老师并没有复现文献，算别的都没问题，只是算有机物，两个有机物差了单纯的一根键的问题，现在算出 ...

量子化学计算根本不看有没有键、形式键级设成几，因此这种描述没有任何意义。看下文
谈谈原子间是否成键的判断问题
http://sobereva.com/414（http://bbs.keinsci.com/thread-9840-1-1.html）

作者
Author: liuxp 时间: 2025-4-29 10:13

Dempey 发表于 2025-4-28 17:15
你的1.png比2.png多一个氢，组成不一样不能比较绝对值，你可以比较吸附能的差异

谢谢刚发现这个问题

作者
Author: liuxp 时间: 2025-4-29 10:15

sobereva 发表于 2025-4-28 17:22
量子化学计算根本不看有没有键、形式键级设成几，因此这种描述没有任何意义。看下文
谈谈原子间是否成键 ...

好的谢谢sob老师

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