计算化学公社

 找回密码 Forget password
 注册 Register

molpro中CAS方法活化轨道输入的问题

查看数: 12855 | 评论数: 14 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2016-8-27 17:54

正文摘要:

各位老师,初学molpro,我想问一下关于molpro的活化空间输入问题,我现在看过Atkins分子对称性那章关于不可约表示,也看了http://www.pci.tu-bs.de/aggericke/PC4e/Kap_V/AO_H2O.htm这一篇关于水分子分子轨道不可约 ...

回复 Reply

SongChao 发表于 Post on 2016-9-2 21:38:36
sobereva 发表于 2016-9-2 20:23
合理。
MO总带有点离域性,不可能MO恰好对应那个AO,总会有其它AO掺入一些。把等值面数值调大可以把关键 ...

好的,非常感谢!
sobereva 发表于 Post on 2016-9-2 20:23:19
SongChao 发表于 2016-9-2 15:34
Sob老师:
     我看完懂了关于不可约表示的东西,重复一个一篇文献中的phenylnitrene分子得到的结构也 ...

合理。
MO总带有点离域性,不可能MO恰好对应那个AO,总会有其它AO掺入一些。把等值面数值调大可以把关键区域突出,而不显示数值较小的区域,这样更有助于指认MO的主要特征。
SongChao 发表于 Post on 2016-8-31 13:45:59
遗忘爱神 发表于 2016-8-31 12:21
molpro手册里有一张就是讲对称性输入的,苯应该是symmetry=x,y,z这样输入吧(如果没记错的话)

谢谢,我找到manual里面输入对称性的说明了
遗忘爱神 发表于 Post on 2016-8-31 12:21:02
SongChao 发表于 2016-8-30 18:13
谢谢您回复
b3lyp/cc-pvtz算出来 一个苯 结果他给认成了C1对称性,所有轨道都是A
这样选活化轨道就都是 ...

molpro手册里有一张就是讲对称性输入的,苯应该是symmetry=x,y,z这样输入吧(如果没记错的话)
SongChao 发表于 Post on 2016-8-31 09:36:54
sobereva 发表于 2016-8-31 01:09
电子态的不可约表示就是各个轨道不可约表示的直积
看过科顿的《群论在化学中的应用》就自然明白了。看书 ...

谢谢sob老师指点,我去看一看cotton的书。
sobereva 发表于 Post on 2016-8-31 01:09:00
SongChao 发表于 2016-8-30 16:59
感谢您的回复。
我懂了,我刚查看了一下用Gaussian算完的out文件,确实会给出轨道对称性的不可约表示。
...

电子态的不可约表示就是各个轨道不可约表示的直积
看过科顿的《群论在化学中的应用》就自然明白了。看书,不要看文章。文章里没有你要的知识,要看书。
sobereva 发表于 Post on 2016-8-31 01:06:58
SongChao 发表于 2016-8-30 18:13
谢谢您回复
b3lyp/cc-pvtz算出来 一个苯 结果他给认成了C1对称性,所有轨道都是A
这样选活化轨道就都是 ...


什么程序?高斯的话初始结构没有做对称化,没判断出D6h点群。
甲醛也是初始结构不满足C2v点群。在gview里先做对称化再保存输入文件。
SongChao 发表于 Post on 2016-8-30 18:13:23
遗忘爱神 发表于 2016-8-30 16:37
你可以从molpro的输出文件里读分子轨道的对称性

谢谢您回复
b3lyp/cc-pvtz算出来 一个苯 结果他给认成了C1对称性,所有轨道都是A
这样选活化轨道就都是A对称性了 这样怎么选呢?
SongChao 发表于 Post on 2016-8-30 16:59:38
sobereva 发表于 2016-8-30 15:14
取的不是原子轨道,就是分子轨道。
只不过分子体系中,有些分子轨道恰好或近似对应于分子轨道,这从轨 ...

感谢您的回复。
我懂了,我刚查看了一下用Gaussian算完的out文件,确实会给出轨道对称性的不可约表示。
但是我一般算完的分子,用凑活的基组算完即使很像是C2v的对称性,键长键角也会有细微的不同,比如算一个甲醛分子,他最后轨道的对称性都是A'和A'',我觉得他应该是一个C2v的对称性,但被划到了Cs对称性,这种情况应该按Cs对称性输入cas轨道吗?
还有如果输入基态结构计算激发态,不知道激发态对称性归为哪类点群,应该按初始结构输入轨道对称性吗?

关于波函数对称性,我只知道闭壳层输入时是A1,sob老师能否推荐几个如何确定开壳层波函数对称性的文章。

谢谢sob老师
遗忘爱神 发表于 Post on 2016-8-30 16:37:11
你可以从molpro的输出文件里读分子轨道的对称性
sobereva 发表于 Post on 2016-8-30 15:14:49
SongChao 发表于 2016-8-30 12:50
感谢sob老师回复。
我没有查到关于波函数对称性公式少一点的文章,能否求一些推荐,告诉我怎么划分开壳 ...


取的不是原子轨道,就是分子轨道。
只不过分子体系中,有些分子轨道恰好或近似对应于分子轨道,这从轨道形状上一看就知道。
不知你说的对称性公式少一点的文章指什么。
不管开壳层闭壳层,HF计算完了都会有轨道不可约表示输出。之后做CASSCF用的只是其中一种自旋的轨道,或者UNO轨道(UHF的无自旋自然轨道)。
SongChao 发表于 Post on 2016-8-30 12:50:38
sobereva 发表于 2016-8-28 02:42
1 波函数不可约表示跟体系净电荷数没直接关系。闭壳层状态所有轨道都是双占据的,此时必定是A1。
2 分子轨 ...

感谢sob老师回复。
我没有查到关于波函数对称性公式少一点的文章,能否求一些推荐,告诉我怎么划分开壳层轨道对称性呢?
另外我看到的关于活化空间选取建议的文章http://bbs.keinsci.com/forum.php?mod=viewthread&tid=126 这个帖子里的几个,说的都是建议选取的原子轨道,molpro例子里面输入的解释也是说选取了原子轨道,我应该怎样确定组成分子轨道的是否是我想要的原子轨道呢?
sobereva 发表于 Post on 2016-8-28 02:42:29
1 波函数不可约表示跟体系净电荷数没直接关系。闭壳层状态所有轨道都是双占据的,此时必定是A1。
2 分子轨道。
3 你可以先用高斯算一下,看看轨道的不可约表示和能级顺序,可视化也比较方便,如《使用Multiwfn观看分子轨道》(http://sobereva.com/269)。跟NBO毛关系也没有。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 21:39 , Processed in 0.167256 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list