计算化学公社

标题: 甲烷分子：molprol里将Td点群转换为C2v点群 [打印本页]

作者
Author: zazzzzzz 时间: 2022-1-5 10:37
标题: 甲烷分子：molprol里将Td点群转换为C2v点群
本帖最后由 zazzzzzz 于 2022-1-5 17:29 编辑

大家好，我是个molpro新手，现在想用molpro计算得到甲烷分子的轨道信息以及molden文件。先在高斯里面得到轨道信息 (, 下载次数 Times of downloads: 49)
在molpro里转换的时候一直出现错误，我猜测是坐标的问题，想转换为C2V点群，但是改变坐标出现过D2点群和CS点群，然后就不知道咋改
molpro输入文件和出错页面如下：

(, 下载次数 Times of downloads: 46)

(, 下载次数 Times of downloads: 32)

作者
Author: jiangning198511 时间: 2022-1-5 16:22
问题中“molpro转换是什么意思”？程序报错的信息是什么？

作者
Author: zazzzzzz 时间: 2022-1-5 17:37

jiangning198511 发表于 2022-1-5 16:22
问题中“molpro转换是什么意思”？程序报错的信息是什么？

谢谢您的回复！
我需要得到甲烷分子的轨道信息和molden文件，后续用到的程序只识别阿贝尔群，甲烷分子是Td点群，需要转化为C2v点群，不知道molpro的输入文件应该怎么写，关于我自己写的输入文件和报错信息已经上传。（报错信息我猜应该是我的坐标问题，导致提示对称性为Cs点群的X对称元素，但是不知道怎么改）

作者
Author: fatpig 时间: 2022-1-5 17:59
如果你在Gaussian里已经得到了 Td 的结构，那放在molpro里直接默认就应该是C2v。

把 symmetry, x, y 这一行删掉呢？

作者
Author: zazzzzzz 时间: 2022-1-5 21:34

fatpig 发表于 2022-1-5 17:59
如果你在Gaussian里已经得到了 Td 的结构，那放在molpro里直接默认就应该是C2v。

把 symmetry, x, y 这 ...

尝试删过了，出来的结果虽然不报错，但不是我想要的点群数据

作者
Author: jiangning198511 时间: 2022-1-6 11:39
你确定你的结构文件书写是正确的吗？元素后面两个逗号是正确的吗？实在不行把用手册里面的例子一点一点的修改

作者
Author: fatpig 时间: 2022-1-6 16:55
是molpro本身的问题，似乎要转成这种写法才行。

我是在GaussView里添加虚原子x，令x-CH4限定在C2v。

写输入文件时x原子可删可不删，都是C2v。

作者
Author: zazzzzzz 时间: 2022-1-6 17:19

fatpig 发表于 2022-1-6 16:55
是molpro本身的问题，似乎要转成这种写法才行。

我是在GaussView里添加虚原子x，令x-CH4限定在C2v。

好的，谢谢您，很有用，我也试了用z坐标，molpro没有出错，但是轨道没有对上，A1对称性相差一个，所以准备用您的办法！再次感谢！

作者
Author: youyno 时间: 2022-1-6 22:27
这种没有办法的，Molpro和Molcas是一样的只能识别阿贝尔群，而Td点群的最高阿贝尔群就是C2v，所以只能按照C2v处理。

作者
Author: zazzzzzz 时间: 2022-1-7 10:08

youyno 发表于 2022-1-6 22:27
这种没有办法的，Molpro和Molcas是一样的只能识别阿贝尔群，而Td点群的最高阿贝尔群就是C2v，所以只能按照C ...

谢谢您的回复，我尝试转换到C2v，但是轨道不简并

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