楼主您好，我想请问一下最后一步直接离去一个氢是怎么做的？ 直接去掉H然后算的能量，一步到位的么？还是说会经历一个能垒的过程？

升腾 发表于 2018-3-8 18:24
谢谢您的回复，然后单位是在图下标注的,截图时没有标记，△E指势垒，单位kcal/mol;△H指反应热，单位kcal ...

确实不应当这样做，解释不出来，要是反应30°C，你用默认25°C，那还可以解释一个温度很接近影响不大，反应温度很高很低就不容易解释了。倒不如直接向审稿人认错改正，如果保留了chk文件，可以很快的用freqchk工具计算出其他温度下的热力学量。

用△E表示活化能是太好的，用△E≠好些。

能量数值最多一位小数就可以了。根本达不到0.1kcal/mol精度。

如果最后两步是这样，根据你的△H，后两步都是很大的正值，这如何解释呢？

liyuanhe211 发表于 2018-3-8 12:26
反应是多少温度就在多少温度下算热力学量。

图中几处不太妥当，比如没有单位，未标明数值是反应能变还是 ...

谢谢您的回复，然后单位是在图下标注的,截图时没有标记，△E指势垒，单位kcal/mol;△H指反应热，单位kcal/mol;最后一步的计算根据参考的众多文献步骤没有问题。然后现在我的问题其实就是这个温度必须选择一个定值计算，但是当时是使用默认的，突然问为什么我就蒙了。

反应是多少温度就在多少温度下算热力学量。

图中几处不太妥当，比如没有单位，未标明数值是反应能变还是活化能，有意义往往是自由能而不是焓。

最后一步芳构化的活化能感觉太高，违背直觉，可能不太应当在模型中直接离去一个氢原子，而是体系里找个物种做个氢转移可能活化能低一些。