利用Dmol3进行茂金属催化剂链转移过渡态优化和能垒计算

葛立恒super

      各位老师好，


      我现在刚接触到MS软件，想利用MS软件的Dmol3模块来计算茂金属催化剂在烯烃聚合过程中链转移的反应能垒。


      为了计算方便，目前反应物模型和生成模型如下：其中反应物模型我是在金属中心上连了一个丙基，同时在金属中心周围放了一个乙烯分子，从而来模拟烯烃分子配位到金属中心上；生成物模型是经过链转移，丙基上的β-H断裂跑到乙烯分子上同时金属中心与丙基断裂形成丙烯分子，乙烯分子的双键中的一个键断裂，一半连氢原子一半连金属中心。我对上述反应物模型和生成模型先进行了结构优化，Vibrational Analysis后发现反应物模型有多个虚频，生成物模型也是类似情况。接下来我应该怎么消除这些虚频呢？我看了一些B站视频，主要是提高运行精度，提高SCF tolerance 的精度这两个方法。我现在正在尝试这两个方式。但是我有一个猜想，有没有可能我的反应物模型和生成物模型本就是不稳定的，就应该是存在多个虚频，但是如何确定具体几个虚频呢。


     希望有这方面的老师能够解答这些疑惑，恳请老师们能够说的详细一点。因为是刚开始接触MS软件和理论计算，如果有不专业的地方，还请各位老师批评指正。麻烦各位老师了。