乙氧基G16优化后的结构出现-800多的虚频，G09D就没有虚频

倪宇晴

用G16对乙氧基进行优化和振动分析，结果出现了-800多的虚频。但是用G09D算同一个输入文件就没有虚频。
请问各位老师，是我的输入文件有问题么？
以下是我的输入文件：
%nprocshared=28
%mem=6GB
%chk=2.chk
# opt freq M062X def2TZVP

Title Card Required

0 2
O                 -1.25762700   -0.36381700    0.00000500
C                 -0.17980700    0.47571600    0.00000500
H                 -0.30322100    1.15874200   -0.86509300
H                 -0.30319500    1.15883600    0.86502500
C                  1.19019500   -0.19461900    0.00000000
H                  1.98770800    0.55461800   -0.00015400
H                  1.30876700   -0.82400200    0.88515600
H                  1.30863200   -0.82423700   -0.88500500

Mikasa

G16和G09默认的积分格点不一致，G16默认ultrafine，G09默认fine

倪宇晴

Mikasa 发表于 2018-3-9 10:58
G16和G09默认的积分格点不一致，G16默认ultrafine，G09默认fine

但这是最小点的优化啊。之前只遇到过100以下的虚频，现在是-800多，都变成过渡态了。请问老师，这也是积分格点导致的么。

Mikasa

倪宇晴 发表于 2018-3-9 11:20
但这是最小点的优化啊。之前只遇到过100以下的虚频，现在是-800多，都变成过渡态了。请问老师，这也是积 ...

对于明尼苏达泛函而言，积分格点影响很大。看这个博文http://sobereva.com/278

liyuanhe211

观察虚频是不是破坏对称性方向的虚频，如果是的话应该在输入文件中破坏对称性。如果是这个情况，那可能是因为G16数值精度高，噪声小，对称性得以保持；而G09有较大的数值噪声，恰好把错误的对称性破坏掉了。

倪宇晴

Mikasa 发表于 2018-3-9 11:55
对于明尼苏达泛函而言，积分格点影响很大。看这个博文http://sobereva.com/278

谢谢老师的帮助。我去看一下博文

倪宇晴

liyuanhe211 发表于 2018-3-9 12:22
观察虚频是不是破坏对称性方向的虚频，如果是的话应该在输入文件中破坏对称性。如果是这个情况，那可能是因 ...

谢谢老师。我去鼓捣一下。

倪宇晴

liyuanhe211 发表于 2018-3-9 12:22
观察虚频是不是破坏对称性方向的虚频，如果是的话应该在输入文件中破坏对称性。如果是这个情况，那可能是因 ...

老师，虚频的方向大致是CH2处的两个CH键一个变长一个变短的运动。
在输入文件里加nosymm之后还是有800多的虚频。
但是加了int=fine之后，虚频就没了。

liyuanhe211

倪宇晴 发表于 2018-3-9 23:33
老师，虚频的方向大致是CH2处的两个CH键一个变长一个变短的运动。
在输入文件里加nosymm之后还是有800多 ...

nosymm 根本不是我说的破坏初猜对称性，nosymm跟本没有这个作用。我是指让初猜的几何结构不具有(近似成立的)对称元素，例如故意改变几何变量使它扭曲。关于nosymm阅读：
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297

倪宇晴

liyuanhe211 发表于 2018-3-11 02:30
nosymm 根本不是我说的破坏初猜对称性，nosymm跟本没有这个作用。我是指让初猜的几何结构不具有(近似成立 ...

。。是这样啊。谢谢老师。
我试了几次，虚频越来越小了，我再继续试试。