|
本帖最后由 qingmang 于 2025-2-3 09:42 编辑
这是一个自学且初学者的计算经验总结贴,经验来源于自己的计算经历与师兄的一些建议,主要为了以后自己回想起来这些经验,会随缘更新。欢迎大家讨论,同时也欢迎各位前辈,同学批评指正(轻喷)
…………………………………………………………………………………………………分……割……线…………………………………………………………………………………………………………………………………………………………………
过渡态找不到,虚频变成甲基旋转或者骨架变动:
0、考虑反应是否存在很明显的活性位点(比如质子,孤对电子)通过加入其他成分(酸、碱、溶剂)来惰化
1、利用扫描,刚性扫描比较快速,一般是扫描中能量最低点(毕竟不大可能直接契合那个虚频的对应的方向),但要关注势能曲线不够平滑或者比较特殊的点,
2、在可能结构附近精细(减少步长)扫描,看看有没有合适的原子位置(扫描后的鞍点)
3、减少计算循环的步长(MexStep = 3),可以让某种二阶不稳定的结构找到可能位置(鞍点和鞍点之间亦有稳定性差别)
中间体找不到,出现甲基旋转或者骨架变动的虚频:
0、此种情况一般是势能面上两个点(鞍点与局部最小点)离得很近导致的
1、积极利用recalc,Maxstep(如:recalc=5,maxstep=3)计算二阶导,可以将两个势能面很接近的位点(鞍点,局部最小点)区分
|
|