计算化学公社

标题: 求助关于molpro2012版说明书中例子的相关问题 [打印本页]

作者
Author: 慢半拍的小陈 时间: 2021-11-11 10:58
标题: 求助关于molpro2012版说明书中例子的相关问题
本帖最后由慢半拍的小陈于 2021-12-28 15:10 编辑

  本人是molpro的初学者，最近在弄清楚molpro2012说明书中的例子。通过搜索论坛和阅读资料，已解决一些问题，还有一些问题因缺乏相应理论知识无法解决。所以来请教各位老师。  1.关于wf卡的定义，主要是不明白如何定义对称性。手册中有一个氧气的例子：”有说到单占据轨道1.6（xz）和1.7（yz）的对称性乘积是4（xy）,因此总波函的对称性是4”我主要是不理解对单占据轨道对称性为何定义为1.6和1.7。
  2.关于occ卡的定义。我有了解到先弄清楚分子对应的点群，然后通过高斯做一个hf计算，在布局分析里找到不可约表示的部分，数点群对应的不可约表示的个数。在molpro2012手册的第四章定义的甲醛的例子occ为512，在第七章定义的甲醛occ为723。为何同一物质occ设置不同？后面我又用高斯对氧气分子做了hf计算。氧气是D2h点群，然后我根据上面描述的方法定义氧气的occ卡，我在输出文件里无法找到氧气点群对应的不可约表示。输出文件附图。
  3.关于closed卡的定义。我还未了解到这个卡的定义，我在论坛里有搜到关于活性空间的选取的论文，是否阅读完就能知道如何定义closed卡。

作者
Author: sobereva 时间: 2021-11-11 17:02
无化学背景，也没量子化学经验，就用molpro是严重不恰当的。本来molpro就是给量子化学领域内行人用的，你的情况如同自动挡的车还没开熟就开手动挡，不费劲才怪。没有对基本原理的理解，光看一堆例子也没用，到头来不能正确举一反三，最后算出一堆错误的数据还是发表不了。

应该想清楚到底有没有必要用molpro，是否明显更容易用的Gaussian或ORCA就足够解决问题。

没有培训专门讲molpro的，除了开发者以前弄的一些workshop。即便北京科音高级量子化学培训班也不会专门讲molpro，本身这东西卖得就太贵，用免费的ORCA、OpenMolcas等程序绝大多数情况都能代替。

作者
Author: zjxitcc 时间: 2021-11-11 17:13
本帖最后由 zjxitcc 于 2021-11-11 17:20 编辑

无化学背景的硕士生用这个是有点被boss坑的，总得组里有会molpro的人带吧，没会的人带的话赶紧弃坑，节约时间，珍惜生命；就算是科班出身的博士，没找到老司机询问的话，独立学习molpro也是很艰难的。molpro是一个你若不懂量子化学知识就根本用不了的软件，这与Gaussian有明显不同。

不如说出你的研究问题，大家可以推荐你其他的研究方式，不是非得molpro，“柔性扫描，单点，频率计算，最后算出能量”这些Gaussian擅长得飞起

作者
Author: 慢半拍的小陈 时间: 2021-11-11 20:45

sobereva 发表于 2021-11-11 17:02
无化学背景，也没量子化学经验，就用molpro是严重不恰当的。本来molpro就是给量子化学领域内行人用的，你的 ...

好的，谢谢老师。我还是想弄清楚说明书里的例子，老师能对我的上述提问解答一下吗？我花了挺多时间的，看书看资料搜论坛什么的。谢谢老师，打扰了！

作者
Author: 慢半拍的小陈 时间: 2021-11-11 20:48

zjxitcc 发表于 2021-11-11 17:13
无化学背景的硕士生用这个是有点被boss坑的，总得组里有会molpro的人带吧，没会的人带的话赶紧弃坑，节约时 ...

我要算烷基与氧气生成过氧烷基这一条路径，最后画出它的能垒图。单参考态方法行不通，需要多参考态方法来计算。

作者
Author: zjxitcc 时间: 2021-11-11 21:19
本帖最后由 zjxitcc 于 2021-11-11 21:41 编辑

慢半拍的小陈发表于 2021-11-11 20:48
我要算烷基与氧气生成过氧烷基这一条路径，最后画出它的能垒图。单参考态方法行不通，需要多参考态方法来 ...

有机体系用UDFT和UCCSD(T)还是靠谱的

作者
Author: 慢半拍的小陈 时间: 2021-11-12 08:24

zjxitcc 发表于 2021-11-11 21:19
有机体系用UDFT和UCCSD(T)还是靠谱的

好的，谢谢老师。老师能帮我解决一下上述molpro卡的定义问题吗？打扰了。

作者
Author: biogon 时间: 2021-11-12 09:10
本帖最后由 biogon 于 2021-11-12 09:12 编辑

慢半拍的小陈发表于 2021-11-11 20:48
我要算烷基与氧气生成过氧烷基这一条路径，最后画出它的能垒图。单参考态方法行不通，需要多参考态方法来 ...

这事用单参考并非不能做，初学者上来就是多参考也是很离谱的，更别说用molpro做

作者
Author: zjxitcc 时间: 2021-11-12 09:29
本帖最后由 zjxitcc 于 2021-11-12 09:31 编辑

慢半拍的小陈发表于 2021-11-12 08:24
好的，谢谢老师。老师能帮我解决一下上述molpro卡的定义问题吗？打扰了。

我不用对称性，即使我用molpro，我也会用MOKIT调用molpro自动做多参考态计算。我甚至觉得molpro手册里很多说明难以阅读，一些技巧过时了。

作者
Author: 慢半拍的小陈 时间: 2021-11-12 09:31

biogon 发表于 2021-11-12 09:10
这事用单参考并非不能做，初学者上来就是多参考也是很离谱的，更别说用molpro做

好的谢谢老师。那老师能帮我看看如何解决molpro卡的定义问题吗？麻烦了~

作者
Author: 慢半拍的小陈 时间: 2021-11-12 09:33

zjxitcc 发表于 2021-11-12 09:29
我不用对称性，即使我用molpro，我也会用MOKIT调用molpro自动做多参考态计算。我甚至觉得molpro手册里很 ...

好的谢谢老师~

作者
Author: biogon 时间: 2021-11-12 10:28

慢半拍的小陈发表于 2021-11-12 09:31
好的谢谢老师。那老师能帮我看看如何解决molpro卡的定义问题吗？麻烦了~

你要是非要用molpro的话就按9楼说的用mokit自动算

作者
Author: zjxitcc 时间: 2021-11-12 10:48
molpro今年也有workshop，就在这两天（昨天好像故障了，改成了今天），你感兴趣可以去听。不过听完这个也只能有个大概的了解，没法学会做手册里的那种古董式计算。
Molpro Workshop开始时间为：下午2:00
会议链接如下：
会议主题：Molpro Workshop
https://meeting.tencent.com/dm/IlLm86RkjvyH
会议 ID：746 554 006

作者
Author: 慢半拍的小陈 时间: 2021-11-12 10:50

zjxitcc 发表于 2021-11-12 10:48
molpro今年也有workshop，就在这两天（昨天好像故障了，改成了今天），你感兴趣可以去听。不过听完这个也只 ...

非常感谢老师~

作者
Author: hebrewsnabla 时间: 2021-11-12 11:09

慢半拍的小陈发表于 2021-11-11 20:48
我要算烷基与氧气生成过氧烷基这一条路径，最后画出它的能垒图。单参考态方法行不通，需要多参考态方法来 ...

这个反应里面的物质没啥对称性啊，不用对称性就行了。

作者
Author: 慢半拍的小陈 时间: 2021-11-12 11:21

hebrewsnabla 发表于 2021-11-12 11:09
这个反应里面的物质没啥对称性啊，不用对称性就行了。

不用对称性的话是直接定义为1吗？我看到有一个过氧乙酸甲酯与氧气反应的输入文件，在wf卡中定义为69，1，1

作者
Author: hebrewsnabla 时间: 2021-11-12 11:41

慢半拍的小陈发表于 2021-11-12 11:21
不用对称性的话是直接定义为1吗？我看到有一个过氧乙酸甲酯与氧气反应的输入文件，在wf卡中定义为69，1， ...

可以用mokit自动产生molpro输入文件

作者
Author: 慢半拍的小陈 时间: 2021-11-13 10:46
有人能解决我的问题吗？顶一下。

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