极低能垒反应IRC中forward和reverse关键词路径相同的问题

lonepairs

各位老师好，我之前复现一个分步Diels-Alder体系，TS-2-N(虚频156.5i，对应分步反应的第二步C6-C12成键)跑IRC时，使用forward和reverse关键词，均跑向INT而非产物。
TS-2-N的两个方向IRC中都出现校正步不收敛(经验不足，步长设置不合理)；
其中reverse(19steps)方向键长增加，forward(11steps)方向键长稍微减小后增加；两个方向能量均单调下降；调整stepsize=4后，IRC跑了更多点(C6-C12键长从2.47A到2.9A，高于INT-N的2.8A)，也出现上述趋势。
文献中对这一步的总结是"almost barrierless"，计算发现自由能垒<2.0kcal/mol，无论是否采用准谐振近似。
计算级别：M06-2X/6-31G* for opt freq IRC; M06-2X/6-311G** for SP; CPCM(dichloromethane); Truhlar’s quasiharmonic correction；复现中计算唯一区别是使用G16而非文献中的G09。
想请问：1.IRC此现象的原因是什么，有没有解决方法？2.这些结果能否验证或否定TS-2-N？3.类似的低能垒过渡态TS-2-X(文献键长2.53A)一直找不到，在2.5A-2.56A以及2.4A-2.6A扫描结果能量单调递增，最后附了一次优化的结果。请问各位老师有没有好的建议？








TS-2-X:

sobereva

在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）
上文第二节说过IRC两个方向相同的问题

mmyy321

你好，我前段时间也在复现这篇文章，并且遇到了和你相同的问题，我当时TS-2-X的过渡态没有找出来（现在也没找出来，主要是后面没再去找）。对于ts-2-n的irc，我当时用#p M062X/6-31G(d) IRC=(calcfc,maxpoints=100,maxcycle=300,stepsize=4,recalc=10) scrf=(cpcm,solvent=dichloromethane) 的指令，并且成功了，你可以也尝试一下。另外如果你ts-2-x成功找出来，可以告诉我方法或者指令吗，谢谢。祝顺利！