计算化学公社

标题: Gaussian 16 A.03 WIN64优化乙烯分子点群出现降级，应该如何修正 [打印本页]

作者
Author: Haike_tyz 时间: 2021-9-23 16:45
标题: Gaussian 16 A.03 WIN64优化乙烯分子点群出现降级，应该如何修正
使用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文件

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

作者
Author: Haike_tyz 时间: 2021-9-23 18:22

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

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

作者
Author: wzkchem5 时间: 2021-9-23 19:33

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

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

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

作者
Author: sobereva 时间: 2021-9-24 09:41

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

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

作者
Author: Haike_tyz 时间: 2021-9-25 20:06
谢谢各位老师，我再去补补基组的基础

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