计算化学公社

标题: 高斯计算过渡态，结果有一个虚频，但感觉模式不对 [打印本页]

作者
Author: 鱼叉君 时间: 2017-11-27 19:14
标题: 高斯计算过渡态，结果有一个虚频，但感觉模式不对
体系是Cu-salen（催化剂，Cu周围有2个O2个N配位）、O2和β-异佛尔酮，想做的过渡态是Cu-salen结合了O2，然后O和β-异佛尔酮上某C成键，即一个（Cu-salen）-O-O-（β-异佛尔酮）的结构，最终能得到氧化异佛尔酮。过渡态结构如图1所示。
设置的时候加长了O-C的键长以及C-H的键长（考虑到O连接上去的话H可能会离开），所以我预期的是能看到O-C和C-H的伸缩振动。目前这次结果是正常结束，查看频率是发现有一个虚频：-25.67，如图2。

查看动画，感觉是个摇摆振动，感觉整个β-异佛尔酮都在摆动。那这个应该不是过渡态吧？

那我应该继续调整初猜结构？调整初猜结构后容易发生502不收敛错，想问下大家对调整初猜结构有啥建议吗？谢谢！~

作者
Author: sobereva 时间: 2017-11-27 20:16
这么点虚频几乎不可能是过渡态

作者
Author: laoman 时间: 2017-11-27 20:40
如果能从文献找到类似的体系，按文献的结果来设置初始结构，没有的话就按下面这个帖子介绍的柔性扫描来确定初始结构：
http://bbs.keinsci.com/forum.php ... &extra=page%3D1
（在“如何通过柔性扫描产生较好的过渡态初猜” 这一节）

作者
Author: 鱼叉君 时间: 2017-11-30 12:43

sobereva 发表于 2017-11-27 20:16
这么点虚频几乎不可能是过渡态

谢谢Sob老师。我另外一个结果得到的唯一虚频是-134.17，能看出来是Cu-O、O-O（如图1 ）的伸缩振动，那这个是过渡态吗？
我用这个构型算了IRC前后态（g16版本算的），都是L502 Convergence criterion not met。
前态：
#p irc=(forward,calcfc,maxpoints=500,stepsize=10,maxcyc=30,lqa) b3lyp/genecp scrf=(solvent=pyridine) SCF=vshift=400
0 2
……
Cu 0
lanl2dz
****
C H O N 0
6-31g*
****

Cu
lanl2dz

错误：Convergence criterion not met。

作者
Author: sobereva 时间: 2017-11-30 14:38

鱼叉君发表于 2017-11-30 12:43
谢谢Sob老师。我另外一个结果得到的唯一虚频是-134.17，能看出来是Cu-O、O-O（如图1 ）的伸缩振动，那这 ...

是不是过渡态看虚频方向，我不清楚你当前研究的反应，若虚频方向和预期的相符，能有8成把握说是

走IRC是最靠谱检验方法，不收敛看
解决SCF不收敛问题的方法
http://sobereva.com/61
不要写maxcyc=30这种多余关键词

作者
Author: YYHHPP 时间: 2024-2-26 22:32
您好，想问一下您这篇关于salen-Cu催化反应的文章发表了吗，想学习一下。

作者
Author: wxyhgk 时间: 2024-2-27 09:17
发结构

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