求助Y(OMe)3几何优化出现虚频，有18个负值

Novice · 发表于 Post on 2020-12-16 16:02:16

wzkchem5 发表于 2020-12-15 20:21
我是说Y-O-C是弯的，不是说Y是突出的。
建议找个教材自学金属有机化学再来做计算吧，这种分子的结构都是 ...

他这个结构挺有意思的，我试着用B3LYP/（6-31g**+SDD)、PBE0和MN15L优化了一下，发现最终结构都是向着图中C-O-Y成一直线的平面三角形跑去。更有意思的是，xtb的gfn2优化能给出我前面贴的那个非直线三角锥坐标，但是对过渡金属配合物可能效果还好点的gfn1优化却不收敛出错！感兴趣可以试试。

wzkchem5 · 发表于 Post on 2020-12-16 16:21:33

Novice 发表于 2020-12-16 16:02
他这个结构挺有意思的，我试着用B3LYP/（6-31g**+SDD)、PBE0和MN15L优化了一下，发现最终结构都是向着图 ...

你能重复出来他的虚频吗？

zjxitcc · 发表于 Post on 2020-12-16 16:30:48

wzkchem5 发表于 2020-12-16 16:21
你能重复出来他的虚频吗？

题主写的基组和赝势不对，例如一些H原子连基函数都是0个，结果不可信

不能与之对比

Novice · 发表于 Post on 2020-12-16 16:48:15

本帖最后由 Novice 于 2020-12-16 16:54 编辑

wzkchem5 发表于 2020-12-16 16:21
你能重复出来他的虚频吗？

哈哈，你说完我还专门去算了下频率，只有一个虚频（在-25左右），是沿着一个甲基转动的方向的。我觉得他这个结构很可能是势能面太平缓，需要减小步长。当然我没做过配合物，对Y也不熟悉，不清楚这个结构本身是不是化学稳定或者合理的。
楼上说的我也注意到了，我算的时候改成正确的输入了。

paramecium86 · 发表于 Post on 2020-12-16 17:01:42

主要问题还是在H的基组没定义。我试了试改了之后直接算还是会有虚频。破坏一点对称性加上calcall就可以解决了。

Novice · 发表于 Post on 2020-12-16 17:09:55

paramecium86 发表于 2020-12-16 17:01
主要问题还是在H的基组没定义。我试了试改了之后直接算还是会有虚频。破坏一点对称性加上calcall就可以解 ...

你这个最后一步的频率计算位移提示还是很大，倒是没有虚频了，肯定是势能面十分平缓导致的。哈哈，这个结构还是有点意思的。一眼看过去感觉是不太合理。

paramecium86 · 发表于 Post on 2020-12-16 17:33:27

Novice 发表于 2020-12-16 17:09
你这个最后一步的频率计算位移提示还是很大，倒是没有虚频了，肯定是势能面十分平缓导致的。哈哈，这个结 ...

嗯是的。不过还是怀疑这个基组不太行。我试了试用我们组以前算这个族的金属有机物比较常用6-311G** /SDD 搭配去算。加上calcall，收敛也正常了。

wzkchem5 · 发表于 Post on 2020-12-16 18:24:34

zjxitcc 发表于 2020-12-16 16:30
题主写的基组和赝势不对，例如一些H原子连基函数都是0个，结果不可信不能与之对比

好吧，原来问题在这里。。。

wzkchem5 · 发表于 Post on 2020-12-17 08:30:25

抱歉啊楼主，我之前说Y-O-C应该是弯的，是从你算出来有一大堆虚频反推出真实结构应该比你优化的结构的对称性低，属于事后诸葛亮。。。我昨天用更高精度的计算（Y用triple zeta，用ZORA描述相对论效应）验证了，确实基本是直的，Y-O-C键角174度。事后推测可能是和O 2p到Y 4d的pi键有关。没想到你算出来虚频另有原因，结果我鲁莽地从虚频得到你结构不合理的错误结论，是我不对。

z472673893 · 发表于 Post on 2020-12-17 10:27:05

wzkchem5 发表于 2020-12-17 08:30
抱歉啊楼主，我之前说Y-O-C应该是弯的，是从你算出来有一大堆虚频反推出真实结构应该比你优化的结构的对称 ...

哈哈，谢谢大家帮助，我主要是做实验的，这几天忙着做实验，然后也是想自己学一下Gaussian，不过还是自己太粗心了，重复了十几次，都没注意没标定H，，一直以为Gaussian很简单，还是要好好学哈哈

z472673893 · 发表于 Post on 2020-12-17 10:28:09

sobereva 发表于 2020-12-16 02:27
此文明确说了
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（h ...

谢谢Sob老师，找到问题关键了。

z472673893 · 发表于 Post on 2020-12-17 10:28:47

zjxitcc 发表于 2020-12-16 16:30
题主写的基组和赝势不对，例如一些H原子连基函数都是0个，结果不可信不能与之对比

终于找到问题关键了，，

z472673893 · 发表于 Post on 2020-12-17 10:29:45

paramecium86 发表于 2020-12-16 17:01
主要问题还是在H的基组没定义。我试了试改了之后直接算还是会有虚频。破坏一点对称性加上calcall就可以解 ...

谢谢帮助

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 求助Y(OMe)3几何优化出现虚频，有18个负值

评分 Rate

评分 Rate

浏览过的版块