计算化学公社

标题: 求问能从Gaussian计算得到的IRC直接提取反应能垒和隧穿宽度信息吗？ [打印本页]

作者
Author: mxr21 时间: 2022-1-7 21:22
标题: 求问能从Gaussian计算得到的IRC直接提取反应能垒和隧穿宽度信息吗？
(1)通过Gaussian计算，其中maxpoints=50,stepsize=20得到了IRC曲线，两端平缓，可以直接读取点的坐标计算反应能垒（过渡态纵坐标-反应物纵坐标），隧穿宽度（生成物横坐标-反应物横坐标）吗？(2)如果可以那么反应物和生成物该怎么寻找（看到公社中帖子提到对IRC端点进行极小值能量优化，但是这样优化出来了就找不到对应的横坐标了，无法计算隧穿宽度）

作者
Author: snljty 时间: 2022-1-7 23:11
本帖最后由 snljty 于 2022-1-8 08:26 编辑

能垒不能。再平缓也没收敛到极小值，差距往往还很明显。另一方面，通常应该对过渡态和底物产物更高级别计算的高精度单点能，与优化和频率分析的相对较低级别的自由能校正项相加。

作者
Author: sobereva 时间: 2022-1-8 06:24
如果两端已经基本跑平了，跑出来的IRC部分就已经足够定义隧穿宽度了

作者
Author: mxr21 时间: 2022-1-8 10:16

sobereva 发表于 2022-1-8 06:24
如果两端已经基本跑平了，跑出来的IRC部分就已经足够定义隧穿宽度了

Bell35模型可以用来描述深度隧穿，在文献中，计算虚频的公式如下，其中m定义为发生隧穿的原子的质量（H的质量），但是通过量纲计算，我发现m的单位为（kn*mol），请问这种情况下如何对m取值呢？
计算过程：
根据量纲相等准则：量纲(RT)=量纲（h/2pi*v*)——J/mol=J*s*v——v=1/mol*s——m=kg*mol
(通过其它文献给出的临界温度Tc和虚频，反推得到m应该在10——19量级）

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