再次感谢老师解答 |
Susan.say 发表于 2024-12-16 20:42 这是个误导性/诱导性问题,不能直接回答 是/不是。“将这个irc里面的过渡态和反应物拿出来分别优化”这个操作和描述是不对的,体现你没有理解你在算什么。因为IRC的最高点就是当前计算级别的过渡态,无需结构优化;IRC的反应物方向 得到的确实是近似的反应物结构,需要结构优化。 当前计算级别下,既然已经找到过渡态,那必然有过渡态的电子能量高于反应物电子能量。你需要计算Gibbs自由能,如果过渡态Gibbs自由能 低于等于 反应物Gibbs自由能,那可以说 当前反应在当前计算级别下是无能垒反应。 如果直接使用当前计算级别的过渡态和反应物结构,用大基组计算电子能量,加上上述计算的热力学校正量,得到大基组下近似的Gibbs自由能,若过渡态Gibbs自由能 低于等于 反应物Gibbs自由能,那同样可以说 当前反应在当前计算策略下是无能垒反应。 不是你描述的那些判定方式。 |
zjxitcc 发表于 2024-12-16 19:18 谢谢老师解答,请问一下,我将这个irc里面的过渡态和反应物拿出来分别优化并用更高级别基组计算能量,如果得出来的过渡态电子能量高于反应物电子能量,可以判定这个反应是无垒反应吗?再次感谢老师 |
王二葛 发表于 2024-12-16 20:11 谢谢老师解答,请问一下,我将这个irc里面的过渡态和反应物拿出来分别优化并用更高级别基组计算能量,如果得出来的过渡态电子能量高于反应物电子能量,可以判定这个反应是无垒反应吗?再次感谢老师 |
首先,请放大你关心的 TS 右侧,能算出虚频 = 沿着虚频方向能量下降 其次,把步长减小 最后是 IRC 的通用注意事项,见 2 楼最后一段 |
zjxitcc 发表于 2024-12-16 19:18 十分感谢老师解答· |
计算上 无能垒反应 一般说的是找不到过渡态,或者找到了过渡态但其Gibbs自由能 低于 反应物Gibbs自由能。而如果找到了过渡态,那么过渡态电子能量一定高于反应物电子能量,也一定高于产物电子能量,否则就是用户搞错了。当你说到无能垒反应时,一定要清楚这点,也要清楚 不能把Gibbs自由能与电子能量混用,要明确指明阐述的物理量。 换句话说,你的问题不能简单地回答为 是/不是无能垒反应。建议你先把IRC计算做得更合理一些,例如forward/reverse方向分成两个任务提交;其中一边不能只跑6个点(那样啥也说明不了),步数多一点。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-14 21:58 , Processed in 0.193845 second(s), 25 queries , Gzip On.