计算化学公社

标题: 从只有一个虚频的初猜开始寻找过渡态，为什么收敛不了而且终点没有了虚频？ [打印本页]

作者
Author: 星斗如盘 时间: 2020-5-17 14:26
标题: 从只有一个虚频的初猜开始寻找过渡态，为什么收敛不了而且终点没有了虚频？
我为了寻找双芳基镍络合物转金属的过渡态，首先做了一个限制性优化，文件freez1。优化之后频率分析只有一个虚频，应该是对应着我需要的过渡态，用它作为初猜进行过渡态优化。但是一直无法收敛，而别结构逐渐偏离过渡态的结构，我取了最后一帧做了一个频率分析，发现虚频消失了。过渡态优化不是在保留虚频的情况下优化吗？为什么会出现这种情况？本人新手，第一次求助，这个结构已经能够优化两个星期了，希望大家有时间看一下。

作者
Author: liyuanhe211 时间: 2020-5-17 17:58
本帖最后由 liyuanhe211 于 2020-5-17 17:59 编辑

某结构是过渡态的条件有两个：有且只有一个虚频并且（无冻结的前提下）【梯度为0】。另外不要过于迷信限制性优化/柔性扫描法找过渡态。能找到固然好，用这办法找不到也没什么奇怪的，还是要进一步调整初猜。

作者
Author: 星斗如盘 时间: 2020-5-17 18:06

liyuanhe211 发表于 2020-5-17 17:58
某结构是过渡态的条件有两个：有且只有一个虚频并且（无冻结的前提下）【梯度为0】。另外不要过于迷信限制 ...

初猜做振动分析有一个虚频，但是优化之后就没有了

作者
Author: sobereva 时间: 2020-5-17 22:50

星斗如盘发表于 2020-5-17 18:06
初猜做振动分析有一个虚频，但是优化之后就没有了

没什么新鲜的，初始结构满足这个条件又不代表优化过程中必定能保持这个条件，更何况你之前还对结构做了冻结，引入了人为因素

结合calcall或者recalc=x再试

写5d关键词毫无意义，何况genecp的时候默认就是用的5d

作者
Author: 星斗如盘 时间: 2020-5-18 09:16
本帖最后由星斗如盘于 2020-5-18 09:34 编辑

sobereva 发表于 2020-5-17 22:50
没什么新鲜的，初始结构满足这个条件又不代表优化过程中必定能保持这个条件，更何况你之前还对结构做了冻 ...

好的，谢谢社长，我尝试一下。我想在补充问一下，如果使用calcall，因为计算的特别慢，还有必要使用maxstep之类的关键字吗？

作者
Author: shenchaoren 时间: 2020-5-18 18:46
楼主的基组搭配不合理.既然Ni已经用了赝势，那么原子序数更大的Br为什么不用赝势而要用全电子基组6-31G(d)?建议两个Br上也用赝势

作者
Author: shenchaoren 时间: 2020-5-18 18:52

星斗如盘发表于 2020-5-18 09:16
好的，谢谢社长，我尝试一下。我想在补充问一下，如果使用calcall，因为计算的特别慢，还有必要使用maxst ...

http://sobereva.com/164 楼主可以看看社长写的

作者
Author: sobereva 时间: 2020-5-18 23:38

星斗如盘发表于 2020-5-18 09:16
好的，谢谢社长，我尝试一下。我想在补充问一下，如果使用calcall，因为计算的特别慢，还有必要使用maxst ...

先尝试maxstep=10结合recalc=3，实在不行再calcall结合maxstep=10

作者
Author: 星斗如盘 时间: 2020-5-20 08:58

shenchaoren 发表于 2020-5-18 18:46
楼主的基组搭配不合理.既然Ni已经用了赝势，那么原子序数更大的Br为什么不用赝势而要用全电子基组6-31G(d)? ...

Br用什么赝势比较合理呢？和镍用同一种可以吗？

作者
Author: shenchaoren 时间: 2020-5-20 09:19

星斗如盘发表于 2020-5-20 08:58
Br用什么赝势比较合理呢？和镍用同一种可以吗？

http://sobereva.com/373 看社长写的《谈谈赝势基组的选用》里面的 “2 什么时候该用赝势基组？”和“3 常用的几种赝势基组”，社长也特意标红了，“Lanl2DZ曾被用得非常广泛，即便是现在还用得很多。但笔者强烈不推荐用Lanl2DZ，因为在如今看来其精度明显偏低。”。里面也写了关于计算配合物时赝势基组与全电子基组的匹配问题

作者
Author: 星斗如盘 时间: 2020-5-20 21:44

shenchaoren 发表于 2020-5-20 09:19
**** 本内容被作者隐藏 ****

我使用的G09貌似不支持社长推荐的Lanl2TZ，显示无法识别Lanl2TZ。。。

作者
Author: 星斗如盘 时间: 2020-5-20 21:49

shenchaoren 发表于 2020-5-20 09:19
**** 本内容被作者隐藏 ****

想再问您一个小问题，我定义了赝势基组，但是在输入计算方法和基组的时候错误的输入了b3lyp/6-311g(d)，这样后面计算的结果是在赝势基组下计算的，还是全电子基组6-311g(d)下计算的呀？

作者
Author: sobereva 时间: 2020-5-21 07:43

星斗如盘发表于 2020-5-20 21:49
想再问您一个小问题，我定义了赝势基组，但是在输入计算方法和基组的时候错误的输入了b3lyp/6-311g(d)， ...

你若没写genecp显然还是用的6-311G*算的

作者
Author: sobereva 时间: 2020-5-21 07:43

星斗如盘发表于 2020-5-20 21:44
我使用的G09貌似不支持社长推荐的Lanl2TZ，显示无法识别Lanl2TZ。。。

显然得自定义基组，仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

作者
Author: 星斗如盘 时间: 2020-5-21 08:54

sobereva 发表于 2020-5-21 07:43
你若没写genecp显然还是用的6-311G*算的

谢谢社长，我查看了一下输入文件，确实是这样，按照6-311*计算的

作者
Author: 星斗如盘 时间: 2020-5-21 08:55

sobereva 发表于 2020-5-21 07:43
显然得自定义基组，仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/ ...

好的，社长，我试一试

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