计算化学公社

标题: 求助，如何计算反应的中间体和过渡态的焓 [打印本页]

作者
Author: zeral 时间: 2021-12-8 11:25
标题: 求助，如何计算反应的中间体和过渡态的焓

根据文章，我想计算反应物以及过渡态的焓，输入文件如下，

[attach]42548[/attach]

文件中未加频率校正因子，焓的计算方法为添加了校正因子后算opt和freq得到焓的热力学校正量A，然后算出溶剂的溶解自由能B，然后提高精度再算一次单点C，体系的焓为A+B+C,

问题1：这个计算焓的方法是否正确（参考Sob老师,但是不知道是否理解有误）

问题2：不加校正因子得到的焓的热力学校正量是否可用

问题3：如果不加溶剂，体系的焓直接就是A+C吗，用不用再加气相下的能量

问题4：文献中只给出了以下数据（DFT/B3LYP,6-311++G*/Lanl2dz,358K），不知道他怎么得到的焓，我得到的焓和他计算得到的不同，是否我的更准确。

文中方法和反应的焓的变化在下图中

(, 下载次数 Times of downloads: 8) (, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 7)

作者
Author: zeral 时间: 2021-12-8 11:31
输入文件因操作失误没有发上去，现在补充。%nprocshared=36
%mem=64GB
%chk=B-11.chk
# opt=(calcfc,noeigentest,tight) freq b3lyp/genecp Int=UltraFine em=GD3BJ scrf=(solvent=water) temperature=358

作者
Author: sobereva 时间: 2021-12-9 06:39

zeral 发表于 2021-12-8 11:31
输入文件因操作失误没有发上去，现在补充。%nprocshared=36
%mem=64GB
%chk=B-11.chk

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

作者
Author: sobereva 时间: 2021-12-9 06:43
1 不是严格的溶剂下的焓，毕竟没法靠隐式溶剂模型得到严格的溶解焓，只能这么凑合近似一下

2 加了校正因子更好，不加也能接受

3 C本来不就是气相下算的能量？何谓用不用再加气相下的能量

4 此文的计算级别很烂，泛函和基组都选择不理想，不要效仿。焓怎么算这里说了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

作者
Author: zeral 时间: 2021-12-9 10:07

sobereva 发表于 2021-12-9 06:43
1 不是严格的溶剂下的焓，毕竟没法靠隐式溶剂模型得到严格的溶解焓，只能这么凑合近似一下

2 加了校正因 ...

好的，感谢Sob的解答，本来很注意发帖了，但是第一次发还是犯错了，下次会更认真检查

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