求助，如何计算反应的中间体和过渡态的焓

zeral · 发表于 Post on 2021-12-8 11:25:02

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

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

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

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

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

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

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

zeral · 发表于 Post on 2021-12-8 11:31:02

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

sobereva · 发表于 Post on 2021-12-9 06:39:06

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

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

sobereva · 发表于 Post on 2021-12-9 06:43:15

1 不是严格的溶剂下的焓，毕竟没法靠隐式溶剂模型得到严格的溶解焓，只能这么凑合近似一下

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

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

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

zeral · 发表于 Post on 2021-12-9 10:07:47

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

2 加了校正因 ...

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 求助，如何计算反应的中间体和过渡态的焓