计算化学公社

标题: [已解决]关于Molpro 无法自动确定对称性的问题 [打印本页]

作者
Author: Finale 时间: 2020-6-11 10:58
标题: [已解决]关于Molpro 无法自动确定对称性的问题
本帖最后由 Finale 于 2020-6-11 13:36 编辑

各位老师，最近开始学习 molpro，目前计算了一个23个原子的体系，在Gasussian中算单点的时候，Gaussian可以确定该分子为C2v点群，但是在 molpro 将其确定为 C1，这就带来很多问题比如，OCC 中的分子轨道不可约表示，我是使用高斯算了一个HF来确定的，两个程序给出的点群不同，OCC出错在所难免
而且，我在 molpro 中使用 symmetry,X,Y 也就是将分子的对称性确定为 C2v，但molpro 不认啊，它不同意，所以总是出现如第一张图所示的错误

假如我关掉对称性，又会出现第二张图的错误，正常来说，括号内是不同不可约表示的轨道占据数依次排列（比如 2，1，1 之类的），
但这里molpro将这些不同不可约表示的轨道占据数全部加起来了（括号内的数），我反复试验后发现，这是由于molpro将分子认定为 C1 所导致的（因为C1没有不同的不可约表示，所以就将所有的轨道数加成一个数，不再依次排列）

我使用的是笛卡尔坐标，难道使用笛卡尔坐标molpro就不能自动确定分子对称性了吗？但为什么我另一个D2h对称性的体系就可以自动确定呢？同样也是笛卡尔坐标
还是说我这个所谓的 C2v 体系对称性并不是C2v？但该体系在Gaussian中就被确定为 C2v对称性

ps：1、附件是输入文件
2、该体系是 103 个电子，多重度=2，但是由于我使用了赝势基组，所以这里电子数我写成了 93

谢谢各位老师不吝赐教！

已解决：使用 Z-矩阵坐标就可以了

作者
Author: 叫个啥名字呢 时间: 2020-6-11 18:13
难道是笛卡尔坐标和Z-矩阵坐标的精度不一样吗？

作者
Author: Peter_zhong 时间: 2020-6-11 18:56
请问各位大佬和楼主：1、某体系怎么判断它是xx点群？比如楼主说的 C2v点群，肉眼观察得出还是通过什么软件？xx点群命名方式也不是很懂 2、知道了是什么点群在高斯里面怎么添加关键字，利用对称大大缩短计算时间呢？谢谢大家~

作者
Author: sobereva 时间: 2020-6-11 22:07

Peter_zhong 发表于 2020-6-11 18:56
请问各位大佬和楼主：1、某体系怎么判断它是xx点群？比如楼主说的 C2v点群，肉眼观察得出还是通过什么软件 ...

根据结构化学知识肉眼就能判断
也可以把结构文件载入Multiwfn程序，载入后屏幕上会直接显示出判断的点群。用gview等程序也可以判断

仔细看
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297（http://bbs.keinsci.com/thread-1355-1-1.html）

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