在设计一个具有ESIPT分子的时候,先模拟了下,发现一个问题,基态IRC过程看起来有点怪,如下图:1为优化后E点,3为优化后的K点,原以为2为K点,结果优化完后跑回E点,最终找到了TS在4位置,跑出来的IRC曲线,就看起来很怪,怀疑是不是找错了。因为E和K的能量相差0.054hartree,算了下有34Kcal/mol,太大了。请大家指点,谢谢。(G16: b3lyp/Def2SVP IRC(calcfc,maxpoints=100,LQA)) (, 下载次数 Times of downloads: 2)