计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12849|回复 Reply: 14
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

各位老师,初学molpro我想问一下关于molpro的活化空间输入问题,我现在看过Atkins分子对称性那章关于不可约表示,也看了http://www.pci.tu-bs.de/aggericke/PC4e/Kap_V/AO_H2O.htm
这一篇关于水分子分子轨道不可约表示的讨论。我看了molpro manual里的例子https://www.molpro.net/info/2015.1/doc/quickstart/node18.html还是有点不太懂。
occ,7,2,3              !7a1, 2b1, 3b2 occupied orbitals
wf,16,1,0}             !16 electrons, Symmetry 1 (A1),singlet
1.      这个里面波函数对称性的1我看到网上有英文的介绍说只要选的电子数等于原子核的电荷数这个就可以写1,这意思是说分子是中性的波函数对称性就可以写1吗?
2.      occ指出的轨道是原子轨道还是分子轨道?molpro这个例子里面说取的原子轨道,但是我看到一篇使用cas方法计算phenylnitrene的文章里说活化空间选择了“7个pi和pi*轨道和一个氮原子的p轨道”,这是什么情况?
3.      如果选取的是分子轨道的话,怎么画出多原子分子的分子轨道图来确定对称性的不可约表示?使用NBO分析之类的方法让软件画出来再选择吗?
问题比较浅陋,还请老师赐教

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2016-8-28 02:42:29 | 只看该作者 Only view this author
1 波函数不可约表示跟体系净电荷数没直接关系。闭壳层状态所有轨道都是双占据的,此时必定是A1。
2 分子轨道。
3 你可以先用高斯算一下,看看轨道的不可约表示和能级顺序,可视化也比较方便,如《使用Multiwfn观看分子轨道》(http://sobereva.com/269)。跟NBO毛关系也没有。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2016-8-30 12:50:38 | 只看该作者 Only view this author
sobereva 发表于 2016-8-28 02:42
1 波函数不可约表示跟体系净电荷数没直接关系。闭壳层状态所有轨道都是双占据的,此时必定是A1。
2 分子轨 ...

感谢sob老师回复。
我没有查到关于波函数对称性公式少一点的文章,能否求一些推荐,告诉我怎么划分开壳层轨道对称性呢?
另外我看到的关于活化空间选取建议的文章http://bbs.keinsci.com/forum.php?mod=viewthread&tid=126 这个帖子里的几个,说的都是建议选取的原子轨道,molpro例子里面输入的解释也是说选取了原子轨道,我应该怎样确定组成分子轨道的是否是我想要的原子轨道呢?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

4#
发表于 Post on 2016-8-30 15:14:49 | 只看该作者 Only view this author
SongChao 发表于 2016-8-30 12:50
感谢sob老师回复。
我没有查到关于波函数对称性公式少一点的文章,能否求一些推荐,告诉我怎么划分开壳 ...


取的不是原子轨道,就是分子轨道。
只不过分子体系中,有些分子轨道恰好或近似对应于分子轨道,这从轨道形状上一看就知道。
不知你说的对称性公式少一点的文章指什么。
不管开壳层闭壳层,HF计算完了都会有轨道不可约表示输出。之后做CASSCF用的只是其中一种自旋的轨道,或者UNO轨道(UHF的无自旋自然轨道)。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

45

帖子

0

威望

503

eV
积分
548

Level 4 (黑子)

5#
发表于 Post on 2016-8-30 16:37:11 | 只看该作者 Only view this author
你可以从molpro的输出文件里读分子轨道的对称性

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2016-8-30 16:59:38 | 只看该作者 Only view this author
sobereva 发表于 2016-8-30 15:14
取的不是原子轨道,就是分子轨道。
只不过分子体系中,有些分子轨道恰好或近似对应于分子轨道,这从轨 ...

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

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

谢谢sob老师

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2016-8-30 18:13:23 | 只看该作者 Only view this author
遗忘爱神 发表于 2016-8-30 16:37
你可以从molpro的输出文件里读分子轨道的对称性

谢谢您回复
b3lyp/cc-pvtz算出来 一个苯 结果他给认成了C1对称性,所有轨道都是A
这样选活化轨道就都是A对称性了 这样怎么选呢?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

8#
发表于 Post on 2016-8-31 01:06:58 | 只看该作者 Only view this author
SongChao 发表于 2016-8-30 18:13
谢谢您回复
b3lyp/cc-pvtz算出来 一个苯 结果他给认成了C1对称性,所有轨道都是A
这样选活化轨道就都是 ...


什么程序?高斯的话初始结构没有做对称化,没判断出D6h点群。
甲醛也是初始结构不满足C2v点群。在gview里先做对称化再保存输入文件。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

9#
发表于 Post on 2016-8-31 01:09:00 | 只看该作者 Only view this author
SongChao 发表于 2016-8-30 16:59
感谢您的回复。
我懂了,我刚查看了一下用Gaussian算完的out文件,确实会给出轨道对称性的不可约表示。
...

电子态的不可约表示就是各个轨道不可约表示的直积
看过科顿的《群论在化学中的应用》就自然明白了。看书,不要看文章。文章里没有你要的知识,要看书。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2016-8-31 09:36:54 | 只看该作者 Only view this author
sobereva 发表于 2016-8-31 01:09
电子态的不可约表示就是各个轨道不可约表示的直积
看过科顿的《群论在化学中的应用》就自然明白了。看书 ...

谢谢sob老师指点,我去看一看cotton的书。

45

帖子

0

威望

503

eV
积分
548

Level 4 (黑子)

11#
发表于 Post on 2016-8-31 12:21:02 | 只看该作者 Only view this author
SongChao 发表于 2016-8-30 18:13
谢谢您回复
b3lyp/cc-pvtz算出来 一个苯 结果他给认成了C1对称性,所有轨道都是A
这样选活化轨道就都是 ...

molpro手册里有一张就是讲对称性输入的,苯应该是symmetry=x,y,z这样输入吧(如果没记错的话)

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2016-8-31 13:45:59 | 只看该作者 Only view this author
遗忘爱神 发表于 2016-8-31 12:21
molpro手册里有一张就是讲对称性输入的,苯应该是symmetry=x,y,z这样输入吧(如果没记错的话)

谢谢,我找到manual里面输入对称性的说明了

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2016-9-2 15:34:18 | 只看该作者 Only view this author
sobereva 发表于 2016-8-31 01:06
什么程序?高斯的话初始结构没有做对称化,没判断出D6h点群。
甲醛也是初始结构不满足C2v点群。在gvie ...

Sob老师:
     我看完懂了关于不可约表示的东西,重复一个一篇文献中的phenylnitrene分子得到的结构也和文献中的结果吻合,但是对于所说的原子轨道还是有一些疑惑。
文献中提到选用的活化空间是苯的7个pi和pi*轨道加上一个“in-plane p AO on nitrogen”,我用Gaussian算完用multiwfn看了所有轨道,很容易找到了最低的7条pi和pi*轨道,但是“in-plane p AO on nitrogen”我只找到了一个比较像的分子轨道,如图一,这是默认isovalue=0.05的图。我觉得这是比较像文献中描述的轨道,然后我把isovalue跳到0.125之后平面就更符合“in-plane p AO on nitrogen”的描述了,如图2。我用这条轨道输入活化空间得到结果和文献吻合,但是不知道我这样调isovalue看轨道像什么原子轨道的方法是不是正确。我还是对这个原子轨道有点不理解。

QQ截图20160902152832.jpg (304.86 KB, 下载次数 Times of downloads: 87)

图1

图1

QQ截图20160902153101.jpg (209.1 KB, 下载次数 Times of downloads: 75)

图2

图2

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

14#
发表于 Post on 2016-9-2 20:23:19 | 只看该作者 Only view this author
SongChao 发表于 2016-9-2 15:34
Sob老师:
     我看完懂了关于不可约表示的东西,重复一个一篇文献中的phenylnitrene分子得到的结构也 ...

合理。
MO总带有点离域性,不可能MO恰好对应那个AO,总会有其它AO掺入一些。把等值面数值调大可以把关键区域突出,而不显示数值较小的区域,这样更有助于指认MO的主要特征。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

28

帖子

0

威望

184

eV
积分
212

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2016-9-2 21:38:36 | 只看该作者 Only view this author
sobereva 发表于 2016-9-2 20:23
合理。
MO总带有点离域性,不可能MO恰好对应那个AO,总会有其它AO掺入一些。把等值面数值调大可以把关键 ...

好的,非常感谢!

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 17:55 , Processed in 0.181956 second(s), 23 queries , Gzip On.

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