A分子存在下B与C反应过渡态优化收敛和寻找问题

JCenter

我的体系中A和B分子可以形成复合物，B和C可以发生亲核取代反应。因此，我研究A存在下B和C的反应能垒的变化（前4周期CHON和Cl原子构成），而这需要寻找过渡态。我的体系共143个原子，体系很大。我首先在b3lyp-d3(bj)/6-31g**结合iefpcm溶剂模型进行了三元体系下的B和C的反应成键的柔性扫描，通过电子能量寻找到了能量最高点，并作为过渡态初猜（图2）。然后在同级别下进行ts寻找，即opt(TS,calcfc,noeigen)。结果发现能量一直在上下波动，无法收敛（图1）。我尝试了1.qst3,效果很差，直接远离了过渡态结构。2.在计算了hessian矩阵后，使用opt(TS,ReadFC,recalc=5,noeigen)，暂时还未收敛，感觉够呛收敛。希望有无大佬和道友给些实用的小技巧和建议让其收敛，感谢大家。
（PS：我目前想到的另一个方法，对柔性扫描上曲线的能量最高点周围的点均做自由能计算，找到自由能最高值，将其所在的结构作为过渡态，不知道是否可取）

JCenter

sai77 发表于 2024-10-16 08:49
可以试试简化结构，跟反应位点离得比较远的原子先拿掉
从扫描曲线看势能面比较平缓，可以试试减小步长上限 ...

非常感谢

sai77

可以试试简化结构，跟反应位点离得比较远的原子先拿掉
从扫描曲线看势能面比较平缓，可以试试减小步长上限
扫描曲线的最高点只能作为初猜，直接作为过渡态肯定不行