计算化学公社

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

[量化理论] 如何判断α MO和β MO是否属于同一轨道?

[复制链接 Copy URL]

210

帖子

1

威望

4360

eV
积分
4590

Level 6 (一方通行)

各位老师好,本人在用Gaussian做三线态双自由基的计算,用的是UB3LYP/6-31G(d)。输出文件用GaussView查看分子轨道,其中有的α MO和β MO从形状和能量上不好判断是否属于同一个轨道,请问这种情况该如何处理?比如:如何判断哪个α MO和β MO同属于HOMO?谢谢解答。
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

4290

帖子

4

威望

9556

eV
积分
13926

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2021-1-29 15:50:19 | 只看该作者 Only view this author
那你做三重态RO计算不就好了
自动做多参考态计算的程序MOKIT

210

帖子

1

威望

4360

eV
积分
4590

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2021-1-29 16:03:38 | 只看该作者 Only view this author
zjxitcc 发表于 2021-1-29 16:50
那你做三重态RO计算不就好了

感谢回复。
RO也算过,目前只是想知道有什么方法可以判断U的结果中α MO和β MO轨道归属的问题。
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

4290

帖子

4

威望

9556

eV
积分
13926

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2021-1-29 16:15:33 | 只看该作者 Only view this author
wangzhe 发表于 2021-1-29 16:03
感谢回复。
RO也算过,目前只是想知道有什么方法可以判断U的结果中α MO和β MO轨道归属的问题。

只有那些芯轨道能用肉眼看出来,alpha/beta轨道长得几乎一样,这些轨道在RHF/ROHF下也有,没有歧义。其他轨道(价轨道)若强行指认,都是存在歧义的,即这个人看起来像,那个人又说不像。因此不建议强行指认。

除非你算出来自旋是严格三重态(即几乎等于三重态ROHF),否则既然是非限制性的方法,就注定一些轨道的alpha/beta空间部分长得不一样,无法归属,这是方法决定的。
自动做多参考态计算的程序MOKIT

210

帖子

1

威望

4360

eV
积分
4590

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2021-1-29 16:29:41 | 只看该作者 Only view this author
zjxitcc 发表于 2021-1-29 17:15
只有那些芯轨道能用肉眼看出来,alpha/beta轨道长得几乎一样,这些轨道在RHF/ROHF下也有,没有歧义。其他 ...

感谢回复。
因为最近在做SOMO-HOMO conversion相关的计算,文献中用的都是U来计算的。如果用RO的话,轨道能量上没法判断是否具有SOMO-HOMO conversion的性质,因为RO计算的轨道都是SOMO能量最高。文献中的分子也计算过了,肉眼判断alpha/beta轨道长得一摸一样,因此我觉得可能文献作者是通过肉眼判断。
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

4290

帖子

4

威望

9556

eV
积分
13926

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2021-1-29 16:50:00 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2021-1-29 16:51 编辑
wangzhe 发表于 2021-1-29 16:29
感谢回复。
因为最近在做SOMO-HOMO conversion相关的计算,文献中用的都是U来计算的。如果用RO的话,轨 ...

哦你是在SOMO-HOMO conversion这种体系里分辨是否空间轨道一致啊,那有办法,你根据alpha/beta分子轨道间的重叠积分(即空间部分的重叠积分)矩阵元判断,数值越接近1.0的表明它俩长得越像。你可以在Multiwfn里找找有没这个功能。
自己写代码变换出空间部分的重叠积分也行,用别的软件也行。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
wangzhe + 5 谢谢

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

210

帖子

1

威望

4360

eV
积分
4590

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2021-1-29 16:54:22 | 只看该作者 Only view this author
zjxitcc 发表于 2021-1-29 17:50
哦你是在SOMO-HOMO conversion这种体系里分辨是否空间轨道一致啊,那有办法,你根据alpha/beta分子轨道间 ...

太感谢了!我去看一下,谢谢!
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

8#
发表于 Post on 2021-1-29 18:50:21 | 只看该作者 Only view this author
Multiwfn里直接就有计算alpha和beta轨道间重叠积分的功能,在手册3.100.5节介绍了。载入含有非限制性波函数的文件,然后主功能100里选子功能5再选1,可以得到诸如下面的输出,可见alpha和beta匹配程度最高的轨道对的序号都告诉你了,重叠积分数值越接近1说明匹配程度越高



北京科音自然科学研究中心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

210

帖子

1

威望

4360

eV
积分
4590

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2021-1-29 20:25:48 | 只看该作者 Only view this author
sobereva 发表于 2021-1-29 19:50
Multiwfn里直接就有计算alpha和beta轨道间重叠积分的功能,在手册3.100.5节介绍了。载入含有非限制性波函数 ...

感谢sob老师回复。刚刚使用Multiwfn成功计算了重叠积分,也验证了SOMO-HOMO conversion,感谢老师开发这么好用的工具。另外,有一个问题想请教一下:如下图中这种,alpha MO的78和79都与beta的78有最大重叠积分,但是beta MO的77在结果中没有显示,这种情况该怎么理解?理论学的不太好,希望老师能给予解答,谢谢。

ovlint.png (70.33 KB, 下载次数 Times of downloads: 94)

ovlint.png
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

10#
发表于 Post on 2021-1-30 00:06:12 来自手机 | 只看该作者 Only view this author
wangzhe 发表于 2021-1-29 20:25
感谢sob老师回复。刚刚使用Multiwfn成功计算了重叠积分,也验证了SOMO-HOMO conversion,感谢老师开发这 ...

我感觉这有点为了套概念而去做分析。SOMO-HOMO conversion这个概念能说明或解释什么样的实验现象? 另外,这个概念所依存的前提首先是单电子近似,其次需要考虑的一个beta虚轨道—而虚轨道对单电子近似下的实际波函数无任何贡献,因而对体系所有性质无任何作用。更何况,根据你的结果,UHF下的alpha-beta配对已经相当模糊了,说服力不够。(另外我不认为应该考虑ROHF,毕竟UHF使得体系基态能量最低,波函数更接近真实波函数—变分原理。更何况在ROHF框架下,SOMO-HOMO conversion这个概念根本就不成立)。还是那句话,SOMO-HOMO conversion有什么重要意义或深刻的insights吗?

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

11#
发表于 Post on 2021-1-30 00:11:39 | 只看该作者 Only view this author
wangzhe 发表于 2021-1-29 20:25
感谢sob老师回复。刚刚使用Multiwfn成功计算了重叠积分,也验证了SOMO-HOMO conversion,感谢老师开发这 ...

Multiwfn也给出了完整的alpha和beta之间的重叠矩阵ovlpmat.txt,看看对应beta 77的第77列中哪些行(alpha 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

210

帖子

1

威望

4360

eV
积分
4590

Level 6 (一方通行)

12#
 楼主 Author| 发表于 Post on 2021-1-30 10:19:50 | 只看该作者 Only view this author
sobereva 发表于 2021-1-30 01:11
Multiwfn也给出了完整的alpha和beta之间的重叠矩阵ovlpmat.txt,看看对应beta 77的第77列中哪些行(alpha ...

多谢老师回复,大変勉強になりました。
Ph.D., Asst. Prof. at UOsaka.
Homepage: wongzit.github.io

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 21:28 , Processed in 0.414906 second(s), 24 queries , Gzip On.

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