过度态不收敛，请教一下各位老师

土拔鼠

各位老师好

体系是110个原子的团簇模型 体系里只有C H O N，我的整个反应文献报道是有两步基元反应，我在找第一步基元反应的过渡态

在优化好结构之后 用关键词 pm6d3 opt=(calcall,gdiis,ts,noeigen) 并添加了隐式水模型，想在低级别下先找一下过渡态作为初始构想在较高级别下再找一下 但是能量一直不收敛 四个指标离收敛标准一直有一到两个数量级

我查了一下卢老师的博文 好像内坐标体系对团簇模型的收敛可能会难一点 所以我改了关键词 opt=(calcall,ts,noeigen,cartesian)

改变之后能量收敛的比较快 结构的变化过程比较符合我设想的过渡态结构 但突然有后期有几步一能量和受力突然升高，变得震荡起来，如图

请问老师 这种情况是否是我的初猜有问题（但是我觉得目前已经很接近我的过渡态，正好在反应的底物与产物之间） 还是我的关键词需要改 有必要继续等下去吗

zjxitcc

等下去是没必要了。如果突跃前几个点结构合理的话，可以提取一个结构出来作为初猜，重新找过渡态。

我个人的经验是PM6（哪怕加色散、溶剂）不是很靠谱···也就有机小分子的简单反应管用。分子一大、自由度一多感觉就不行了。
110个原子，初猜结构构建得好一点，用g16算，然后任务分两步，
（1）PM7 freq 隐式溶剂
（2）opt=(rcfc,ts,noeigen) freq B3LYP/6-31G(d) em=GD3 隐式溶剂 g09default
感觉这种更靠谱。

土拔鼠

zjxitcc 发表于 2019-8-24 16:24
等下去是没必要了。如果突跃前几个点结构合理的话，可以提取一个结构出来作为初猜，重新找过渡态。

我个 ...

感谢 我看看用B3LYPD3/6-31G*试试 我发现后面的结构能量一直在升高 我只有G09 所以用不了PM7

sobereva

图快的话建议Gaussian+xtb的组合代替PM6/7，对于某些情况比PM6/7明显更可靠
将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/421（http://bbs.keinsci.com/thread-10106-1-1.html）

如果继续用PM6D3，也可以取之前优化过程中能量较低，而且已经基本靠谱的结构继续优化过渡态，但是再带上maxstep=5，避免中途跑飞

土拔鼠

sobereva 发表于 2019-8-25 00:30
图快的话建议Gaussian+xtb的组合代替PM6/7，对于某些情况比PM6/7明显更可靠
将Gaussian与Grimme的xtb程序 ...

谢谢社长