甲烷分子：molprol里将Td点群转换为C2v点群

zazzzzzz · 发表于 Post on 2022-1-5 10:37:21

本帖最后由 zazzzzzz 于 2022-1-5 17:29 编辑

大家好，我是个molpro新手，现在想用molpro计算得到甲烷分子的轨道信息以及molden文件。先在高斯里面得到轨道信息

在molpro里转换的时候一直出现错误，我猜测是坐标的问题，想转换为C2V点群，但是改变坐标出现过D2点群和CS点群，然后就不知道咋改
molpro输入文件和出错页面如下：

jiangning198511 · 发表于 Post on 2022-1-5 16:22:32

问题中“molpro转换是什么意思”？程序报错的信息是什么？

zazzzzzz · 发表于 Post on 2022-1-5 17:37:10

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

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

fatpig · 发表于 Post on 2022-1-5 17:59:55

如果你在Gaussian里已经得到了 Td 的结构，那放在molpro里直接默认就应该是C2v。

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

zazzzzzz · 发表于 Post on 2022-1-5 21:34:51

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

把 symmetry, x, y 这 ...

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

jiangning198511 · 发表于 Post on 2022-1-6 11:39:07

你确定你的结构文件书写是正确的吗？元素后面两个逗号是正确的吗？实在不行把用手册里面的例子一点一点的修改

fatpig · 发表于 Post on 2022-1-6 16:55:33

是molpro本身的问题，似乎要转成这种写法才行。

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

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

zazzzzzz · 发表于 Post on 2022-1-6 17:19:06

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

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

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

youyno · 发表于 Post on 2022-1-6 22:27:44

这种没有办法的，Molpro和Molcas是一样的只能识别阿贝尔群，而Td点群的最高阿贝尔群就是C2v，所以只能按照C2v处理。

zazzzzzz · 发表于 Post on 2022-1-7 10:08:30

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

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

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

[Molpro] 甲烷分子：molprol里将Td点群转换为C2v点群