Gaussian 16 A.03 WIN64优化乙烯分子点群出现降级，应该如何修正

Haike_tyz · 发表于 Post on 2021-9-23 16:45:00

使用GView6.0.16绘制乙烯分子，送G16 A.03 WIN64按如下关键词计算
# opt b3lyp/6-31g geom=connectivity
得到的分子显示属于点群 C2h，然而根据理论知识乙烯分子应当是C2v，同学用G09计算结果显示C1（更离谱了）
使用Gview去看C-H键键长，H-C=C键键角，H-C-H键角均无问题，log文件底部数据也显示分子偶极矩为000（这里那位同学的G09输出文件虽然键长键角也没问题，但出现了偶极矩，令人迷惑的问题增加了）
想知道，是哪里出了问题？
附件为我的计算结果log文件

wzkchem5 · 发表于 Post on 2021-9-23 17:29:54

菜单Edit -> Symmetrize，把判断对称性的阈值调低即可。
另外纠正几个问题：（1）乙烯的点群是D2h；（2）计算必须加极化函数，不能用6-31G，至少要用6-31G*

Haike_tyz · 发表于 Post on 2021-9-23 18:22:56

wzkchem5 发表于 2021-9-23 17:29
菜单Edit -> Symmetrize，把判断对称性的阈值调低即可。
另外纠正几个问题：（1）乙烯的点群是D2h；（2） ...

谢谢老师，了解错误了
对以上我有一个疑问，非环分子，尤其是小分子添加极化函数的作用是什么呢？

wzkchem5 · 发表于 Post on 2021-9-23 19:33:13

Haike_tyz 发表于 2021-9-23 11:22
谢谢老师，了解错误了
对以上我有一个疑问，非环分子，尤其是小分子添加极化函数的作用是什么呢？

极化函数的作用是：
（1）让成键轨道有更多的变形自由度，可以更加精细地调整自己的形状，来让能量最小化；
（2）对于post-HF计算，还有一个额外作用：提供角向相关能。
只要是分子，都必须至少对非氢原子加极化函数。原子原则上不用加，但是原子本来计算量就小，没必要特意为原子去掉极化函数，而且你但凡要和分子的能量求差，原子和分子的理论级别必须统一。这和是否非环、是否是小分子没有任何关系，不知道你为什么单单对“非环分子，尤其是小分子”添加极化函数的作用有疑问

liyuanhe211 · 发表于 Post on 2021-9-23 22:31:37

没有G16W没法尝试，但高斯的默认symmetry判定过于严格，优化的收敛限不够低就判定不了高点群，可以用IOp(2/17=N)和IOp(2/18=N)来调整symmetry的判据严格程度。

sobereva · 发表于 Post on 2021-9-24 09:41:55

Haike_tyz 发表于 2021-9-23 18:22
谢谢老师，了解错误了
对以上我有一个疑问，非环分子，尤其是小分子添加极化函数的作用是什么呢？

好好把此文看了，要不然做的计算大多数都可能是垃圾数据
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

Haike_tyz · 发表于 Post on 2021-9-25 20:06:58

谢谢各位老师，我再去补补基组的基础

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

[新手求助] Gaussian 16 A.03 WIN64优化乙烯分子点群出现降级，应该如何修正