考虑受阻转子时，对柔性扫描计算的疑问

400__LUX

      本人学习MESS的过程中，遇到考虑低频受阻转子的问题，故基于文献所给结构和计算级别，对HCOOCH2分子的C-O键进行关于二面角的柔性扫描，尝试复现论文中的扫描结果，最终扫描结果和计算结果如下：





      在考虑受阻转子时，转动组的对称数Symmetry（CH2）为2，则在考虑输入Potential时，需要取360°/(2*10°)=18个点，但扫描结果显示在CH2旋转360°的过程中两个峰具有明显不同，参考其他帖子中的讨论（参考帖），我想问的问题如下：1.该扫描结果的不同是否是因为柔性扫描“滞后性”的结果；2.从比对结果上来看，考虑输入Potential的18个点时是应该以第一个180°为主，直接舍弃后面的180°吗？；3.类似地，知道了转动组的对称数Symmetry=n后，得到相应的取点数360°/(n*10°)=m，是否只需要计算m个点，而不用进行360°的扫描；4.对于涉及到T1诊断值较大的体系，说明其具有一定的多参考态特征，是否应该用多参考态方法进行柔性扫描？
     另外想问一个与本题不太相关的问题：5.我所研究的体系的TS虽然T1诊断值较大，但IRC过程能够正确地连接反应复合物和产物，是否还需要考虑用多参考态方法进行重新优化？

     烦请各位大佬帮忙解答，谢谢。

wzkchem5

我不太确定这个扫描曲线不连续问题怎么解决最好，但“IRC过程能够正确地连接反应复合物和产物”和是否有必要用多参考态方法之间没有任何关系。另外多参考态性质不要用T1诊断值判断，T1诊断值衡量的是HF轨道对于post-HF计算而言有多差，不是直接衡量多参考态性质的，说T1诊断值是衡量多参考态性质的容易被审稿人喷。应当结合FOD、MP2自然轨道占据数等综合判断。如果确实有明显的多参考态性质，可以考虑用CASPT2优化、做扫描，这么小的体系用CASPT2完全算得动。

400__LUX

wzkchem5 发表于 2024-7-2 22:00
我不太确定这个扫描曲线不连续问题怎么解决最好，但“IRC过程能够正确地连接反应复合物和产物”和是否有必 ...

好的，谢谢老师。目前我还未开始接触多参考态计算，看到一篇文献提到了TS结构和IRC的问题，提问题4时正好想到了，可能是我理解上存在错误，我在回头好好看看文献，学习一下多参考态理论。