计算化学公社

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

[综合交流] HOMO轨道全占,HOMO-1轨道半占是否合理

[复制链接 Copy URL]

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

开壳层阴离子体系,在Cs对称下计算得到分子HOMO轨道全占,HOMO-1轨道半占,尝试了不同的方法:UHF,UBP96,UPBE1PBE,UM06L,CCSD,都是一样的趋势图中8a''为半占SOMO轨道能量较7a‘(HOMO)全占轨道低,光电子能谱结果表明目前HOMO的beta电子被激发,与目前能级顺序一致,也就是HOMO上beta电子离去,形成三重激发态,但目前结果电子占据不满足能量最低原理。
如果用RO方法可以避免这个问题,但ROPBE1PBE能量较UPBE1PBE高3kJ/mol,且需要从HOMO-1而不是HOMO激发电子。
该体系中性平衡结构为Cinfv,三重态,也就是现在的8a''和7a‘在中性结构中为简并pi轨道,均为alpha全占。可以看出,添加一个电子后简并pi轨道发生JT变形,能级分裂为8a''和7a‘轨道,导致8a''的alpha轨道能量低于全占7a‘轨道。
请教大家这样的结果合理吗?

感谢耐心帮助!

level.jpg (45.39 KB, 下载次数 Times of downloads: 51)

level

level

1664

帖子

5

威望

4772

eV
积分
6536

Level 6 (一方通行)

喵星人

2#
发表于 Post on 2020-9-22 14:16:56 | 只看该作者 Only view this author
先用stable检查波函数稳定性,另外确认不同自旋多重度的可能性

4279

帖子

4

威望

9477

eV
积分
13836

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2020-9-22 14:18:39 | 只看该作者 Only view this author
还是要多交待一点计算细节,例如你这个CCSD结果,别是看成了UHF轨道能量去了。UCCSD本身是没有轨道能级之说的,高斯fchk文件里默认存的是UHF轨道能量。要看UCCSD自然轨道得另算。
自动做多参考态计算的程序MOKIT

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2020-9-22 14:40:46 | 只看该作者 Only view this author
zjxitcc 发表于 2020-9-22 14:18
还是要多交待一点计算细节,例如你这个CCSD结果,别是看成了UHF轨道能量去了。UCCSD本身是没有轨道能级之说 ...

非常感谢耐心解答,确实都是用fchk的轨道,我再算一下CCSD的自然轨道。

4279

帖子

4

威望

9477

eV
积分
13836

Level 6 (一方通行)

MOKIT开发者

5#
发表于 Post on 2020-9-22 14:57:44 | 只看该作者 Only view this author
wzh 发表于 2020-9-22 14:40
非常感谢耐心解答,确实都是用fchk的轨道,我再算一下CCSD的自然轨道。

各种UDFT和UHF的波函数稳定性也要注意。
自动做多参考态计算的程序MOKIT

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-9-22 15:15:48 | 只看该作者 Only view this author
喵星大佬 发表于 2020-9-22 14:16
先用stable检查波函数稳定性,另外确认不同自旋多重度的可能性

谢谢耐心解答,波函数是稳定的,不同多重度也测试了

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2020-9-22 15:20:04 | 只看该作者 Only view this author
本帖最后由 wzh 于 2020-9-23 12:08 编辑
zjxitcc 发表于 2020-9-22 14:57
各种UDFT和UHF的波函数稳定性也要注意。

还有个问题,光电子能谱应该是从非定域轨道激发吧?

4279

帖子

4

威望

9477

eV
积分
13836

Level 6 (一方通行)

MOKIT开发者

8#
发表于 Post on 2020-9-22 15:22:29 | 只看该作者 Only view this author
wzh 发表于 2020-9-22 15:20
还有个问题,自然轨道是定域的,光电子能谱应该是从非定域轨道激发吧?

自然轨道不一定是局域的,没有任何限制,算出来是啥就是啥。例如多并苯体系的自然轨道是十分离域的。“光电子能谱应该是从非定域轨道激发”你这种表述已经用化学直觉默认是正则轨道图像了。
自动做多参考态计算的程序MOKIT

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2020-9-23 09:29:04 | 只看该作者 Only view this author
本帖最后由 wzh 于 2020-9-23 22:47 编辑
zjxitcc 发表于 2020-9-22 15:22
自然轨道不一定是局域的,没有任何限制,算出来是啥就是啥。例如多并苯体系的自然轨道是十分离域的。“光 ...

谢谢您,又引伸出一个问题如CH4的正则轨道是满足Td对称性的,也就是有一个三重简并轨道,而用自然键轨道则是4重简并的。二者轨道能级也不同,如何考虑激发态的计算呢?

4279

帖子

4

威望

9477

eV
积分
13836

Level 6 (一方通行)

MOKIT开发者

10#
发表于 Post on 2020-9-23 10:12:55 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2020-9-23 10:16 编辑
wzh 发表于 2020-9-23 09:29
谢谢您,又引伸出一个问题如CH4的正则轨道是满足Td对称性的,也就是有一个三重简并轨道,而用自然轨道 ...

自然轨道没有能级,只有占据数,用GaussView打开自然轨道的fchk文件看轨道,电子上下箭头旁的数字都是>=0的小数,那是占据数。你说能级,可能你看错了。激发态计算本来就不一定非得考虑对称性,那只是习惯对小分子用对称性的一类人的做法而已。
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

9740

eV
积分
21936

Level 6 (一方通行)

11#
发表于 Post on 2020-9-23 11:05:35 | 只看该作者 Only view this author
wzh 发表于 2020-9-23 09:29
谢谢您,又引伸出一个问题如CH4的正则轨道是满足Td对称性的,也就是有一个三重简并轨道,而用自然轨道 ...

你把自然轨道和局域轨道搞混了吧,CH4的自然轨道也是一个单独的轨道和一组三重简并的轨道
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

1万

帖子

0

威望

9740

eV
积分
21936

Level 6 (一方通行)

12#
发表于 Post on 2020-9-23 11:08:34 | 只看该作者 Only view this author
完全正常。Kohn-Sham理论只要求alpha、beta电子分别符合aufbau规则,不要求两种自旋放在一起仍然符合aufbau规则。所以只要HOMO-1对应的空轨道在HOMO之上,就没问题。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
wzh + 3 赞!

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2020-9-23 11:21:06 | 只看该作者 Only view this author
wzkchem5 发表于 2020-9-23 11:08
完全正常。Kohn-Sham理论只要求alpha、beta电子分别符合aufbau规则,不要求两种自旋放在一起仍然符合aufbau ...

非常感谢

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

14#
 楼主 Author| 发表于 Post on 2020-9-23 11:21:43 | 只看该作者 Only view this author
zjxitcc 发表于 2020-9-23 10:12
自然轨道没有能级,只有占据数,用GaussView打开自然轨道的fchk文件看轨道,电子上下箭头旁的数字都是>=0 ...

非常感谢

17

帖子

0

威望

40

eV
积分
57

Level 2 能力者

15#
 楼主 Author| 发表于 Post on 2020-9-23 12:00:43 | 只看该作者 Only view this author
本帖最后由 wzh 于 2020-9-23 12:02 编辑
wzkchem5 发表于 2020-9-23 11:05
你把自然轨道和局域轨道搞混了吧,CH4的自然轨道也是一个单独的轨道和一组三重简并的轨道

http://mu****g.com/t-4370956-1”木虫不显示了
这个帖子是算的自然键轨道,四重简并的。

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 18:33 , Processed in 0.219485 second(s), 26 queries , Gzip On.

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