计算化学公社

标题: 乙氧基G16优化后的结构出现-800多的虚频，G09D就没有虚频 [打印本页]

作者
Author: 倪宇晴 时间: 2018-3-9 09:47
标题: 乙氧基G16优化后的结构出现-800多的虚频，G09D就没有虚频
本帖最后由倪宇晴于 2018-3-9 10:03 编辑

用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

作者
Author: Mikasa 时间: 2018-3-9 10:58
G16和G09默认的积分格点不一致，G16默认ultrafine，G09默认fine

作者
Author: 倪宇晴 时间: 2018-3-9 11:20

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

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

作者
Author: Mikasa 时间: 2018-3-9 11:55

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

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

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

作者
Author: 倪宇晴 时间: 2018-3-9 14:05

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

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

作者
Author: 倪宇晴 时间: 2018-3-9 14:07

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

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

作者
Author: 倪宇晴 时间: 2018-3-9 23:33

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

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

作者
Author: liyuanhe211 时间: 2018-3-11 02:30

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

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

作者
Author: 倪宇晴 时间: 2018-3-11 10:05

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

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

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