计算化学公社

标题: 两个简单底物的pka计算和预期不符的求助 [打印本页]

作者
Author: taciturn__ 时间: 2025-3-25 23:44
标题: 两个简单底物的pka计算和预期不符的求助
最近在用sob老师量化中级班的方法算pka，级别稍微改了下，是CBSQB3，M052x/def2svp对应的是E_gas的能量，M052x/def2svp-smd(acetonitrile)对应的是E_smd的能量，得到的结果按课上的公式：
ΔG=(E_A_CBSQB3-E_AH_CBSQB3)+(E_A_smd-E_A_gas)-(E_HA_smd-E_HA_gas)-270.29
pka=ΔG*0.733

得到结果如图。算的两个底物一个是质子化的苯胺N上的H，一个是质子化的氧化吡啶O上的H，但是结果算出来和实验值10附近都差的很多
(, 下载次数 Times of downloads: 22) (, 下载次数 Times of downloads: 19)

不知道问题出在哪？

作者
Author: wzkchem5 时间: 2025-3-26 08:30
这个H+溶解自由能数据确定是乙腈里的，不是其他溶剂里的是吧？
另外可以考虑做显式溶剂化

作者
Author: imasen 时间: 2025-3-26 09:42
苯胺的共轭酸的pKa怎么会在10附近呢，是不是看到的是pKb呢？

作者
Author: taciturn__ 时间: 2025-3-26 10:37

wzkchem5 发表于 2025-3-26 08:30
这个H+溶解自由能数据确定是乙腈里的，不是其他溶剂里的是吧？
另外可以考虑做显式溶剂化

哦谢谢王老师！是这个我没考虑到
更正后的结果不做显式溶剂也至少合理了

作者
Author: taciturn__ 时间: 2025-3-26 10:37

imasen 发表于 2025-3-26 09:42
苯胺的共轭酸的pKa怎么会在10附近呢，是不是看到的是pKb呢？

我在ibond上查的结果就是这样，应该是乙腈的原因

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