计算化学公社

标题: 过度态不收敛，请教一下各位老师 [打印本页]

作者
Author: 土拔鼠 时间: 2019-8-24 11:27
标题: 过度态不收敛，请教一下各位老师
各位老师好

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

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

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

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

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

作者
Author: zjxitcc 时间: 2019-8-24 16:24
本帖最后由 zjxitcc 于 2019-8-24 16:26 编辑

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

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

作者
Author: 土拔鼠 时间: 2019-8-24 22:50

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

我个 ...

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

作者
Author: sobereva 时间: 2019-8-25 00:30
图快的话建议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，避免中途跑飞

作者
Author: 土拔鼠 时间: 2019-8-25 12:33

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

谢谢社长

欢迎光临计算化学公社 (http://bbs.keinsci.com/)