醋酸根几何优化出现问题

LPY12345

各位老师好，我今天优化一个醋酸根负离子，但是得到的结构十分诡异，O-C-O的键角明显太小，而且换其他基组泛函包括加弥散都不能解决这个问题，而且我检查了输入文件也看不哪里有问题，希望大家解惑。

%chk=OAc.chk
#p opt freq scrf=(smd,solvent=toluene) empiricaldispersion=gd3bj temperature=383.15 pbe1pbe/def2svp

OAc

1 1
O                  4.36267743    0.44239041    3.20051174
C                  4.72171155    0.41336643    4.36507497
C                  6.16998931    0.56755605    4.76836303
H                  6.77065939    0.86281749    3.89935003
H                  6.26998194    1.31043574    5.57222240
H                  6.54012906   -0.39187234    5.16138082
O                  3.95284003    0.24596128    5.40431207

axuanyan

醋酸根电荷是-1吧

zjxitcc

电荷写错了。另外，阴离子必须要加弥散函数，而且你这基组对于这么小的分子太小了。用def2TZVPD, ma-TZVP, 6-311+G(d)等为宜。

LPY12345

多谢各位老师回复，我献丑了，我居然一直以为Gaussian里输入电荷不需要加正负号。。。。。。。。

LPY12345

zjxitcc 发表于 2020-12-1 21:04
电荷写错了。另外，阴离子必须要加弥散函数，而且你这基组对于这么小的分子太小了。用def2TZVPD, ma-TZVP,  ...

老师，这个醋酸根我不是单独算的，这个醋酸根实际上是给一个配合物大体系做能量补偿用的。而大体系的几何优化过渡金属用def2tzvp.，其他元素用def2svp。因此醋酸根也得跟着用def2svp，我的问题是如果我用自定义基组把醋酸根的氧定义成def2svpD，醋酸根其它元素还用def2svp，这样是否在最后核算能量的时候是否也可以做到合理

zjxitcc

LPY12345 发表于 2020-12-1 23:08
老师，这个醋酸根我不是单独算的，这个醋酸根实际上是给一个配合物大体系做能量补偿用的。而大体系的几何 ...

你这情况比较合理的方案：醋酸根中的C, H: def2SVP; O: def2SVPD。无论是单独算醋酸根，还是算复合物（该复合物中含醋酸根），都得这样。这就意味着，如果你的复合物如果一开始没给氧用def2SVPD，现在用了，结构就得再优化。

LPY12345

zjxitcc 发表于 2020-12-1 23:17
你这情况比较合理的方案：醋酸根中的C, H: def2SVP; O: def2SVPD。无论是单独算醋酸根，还是算复合物（该 ...

好的，多谢老师指点