[修改于 7年4个月前 - 2017/08/19 21:58:03]
引用 ne555:那个，肟基不是应该包括那个碳原子嘛。
肟基（=NOH）许多人都比较陌生
首先鼓励楼主发帖，但有一些存在的问题我必须讲一讲。（原文献没有看（其实我手头也有一篇计算pKa的综述，不过我也没有看））
我个人觉得使用DFT方法计算应该算“从头算”，半经验这种说法感觉不太合适。
本质上来讲，计算pKa就是要得到电离过程的准确的Gibbs自由能变。因此应当选择足够准确的理论水平进行计算。
楼主使用的泛函B3LYP应该没太大问题（不过计算有机物热力学量更推荐使用M06-2X），但是基组的使用很不合适。6-31+G这个基组没有极化函数，很不可取。通常而言，基组中的极化函数是不能缺少的，而弥散函数则在几何优化及频率计算时可以不加。为了得到准确的热力学量，还应当在高理论水平下计算单点能（你的文中提到没有算氢离子是因为该物种无电子，所以没有电子能，热力学校正部分只存在平动能\(3/2 RT\)）。
所以我建议的计算的理论水平是，在M062X/defSVP水平下优化和频率分析（使用相应的校正因子），在B2PLYP-D3/def2TZVP水平下计算单点能（如果计算资源吃紧，那就用M062X/def2TZVP）。
关于溶剂模型，目前比较推荐使用SMD模型。而IEFPCM比Gaussian默认的PCM模型还差，尽量别用吧23333
其实我想说的是，很多文献使用的计算方法并不一定合理。除非是一定要重复文献，自己动手算的时候可以自己考虑一下有没有更合理的理论方法。
前面提到的那个综述是doi: 10.1002/wcms.1218，有兴趣可以看一看。
引用 last-order:半经验法我也觉得不准确，文献中也没有这么说，但毕竟需要实验数据，所以我也拿不准，本来就像看看大家的意见
首先鼓励楼主发帖，但有一些存在的问题我必须讲一讲。（原文献没有看（其实我手头也有一篇计算pKa的综述，不过我也没有看））
我个人觉得使用DFT方法计算应该算“从头算”，半经验这种说法感觉不太合适。
本质上来讲，计算pKa就是要得到电离过程的准确的Gibbs自由能变。因此应当选择足够准确的理论水平进行计算。
楼主使用的泛函B3LYP应该没太大问题（不过计算有机物热力学量更推荐使用M06-2X），但是基组的使用很不合适。6-31+G这个基组没有极化函数，很不可取。通常而言，基组中的极化函数是不能缺少的，而弥散函数则在几何优化及频率计算时可以不加。为了得到准确的热力学量，还应当在高理论水平下计算单点能（你的文中提到没有算氢离子是因为该物种无电子，所以没有电子能，热力学校正部分只存在平动能\\(3/2 RT\\)）。
所以我建议的计算的理论水平是，在M062X/defSVP水平下优化和频率分析（使用相应的校正因子），在B2PLYP-D3/def2TZVP水平下计算单点能（如果计算资源吃紧，那就用M062X/def2TZVP）。
关于溶剂模型，目前比较推荐使用SMD模型。而IEFPCM比Gaussian默认的PCM模型还差，尽量别用吧23333
其实我想说的是，很多文献使用的计算方法并不一定合理。除非是一定要重复文献，自己动手算的时候可以自己考虑一下有没有更合理的理论方法。
前面提到的那个综述是doi: 10.1002/wcms.1218，有兴趣可以看一看。
引用 三硝基二甲苯:
高质量发帖啊。鼓励一下楼主（没KCB随便说说，手动滑稽）引用 ne555:那个，肟基不是应该包括那个碳原子嘛。
肟基（=NOH）许多人都比较陌生
我不是很清楚pKa = ΔGaq/(RT ln 10) ) = ΔGaq/(2.303RT)这个式子是怎么来的，还希望楼主能给细讲一下。
这个半经验法还是蛮不错的，是不是对所有羟基氢都适用？
引用 ne555:明白了引用 三硝基二甲苯:看了一下，确实有碳，闹笑话了
高质量发帖啊。鼓励一下楼主（没KCB随便说说，手动滑稽）引用 ne555:那个，肟基不是应该包括那个碳原子嘛。
肟基（=NOH）许多人都比较陌生
我不是很清楚pKa = ΔGaq/(RT ln 10) ) = ΔGaq/(2.303RT)这个式子是怎么来的，还希望楼主能给细讲一下。
这个半经验法还是蛮不错的，是不是对所有羟基氢都适用？
ΔG=-RTlnK
pK=-lgK
另外还是叫从头算法吧
不同种类化合物的羟基需要不同的常数，比如醇和羧酸，可以看看帖子最后的文献
| 时段 | 个数 |
|---|---|
| {{f.startingTime}}点 - {{f.endTime}}点 | {{f.fileCount}} |
200字以内，仅用于支线交流，主线讨论请采用回复功能。